RU EN
lrg rst sml

Авторизация

Запомнить

OS

35%Windows XP
34%Windows 7
8%Windows NT
6%Mac OS X
5%Linux
floors_1.pngSvetotechnika N3 2005
Главная arrow Публикации arrow К теоретическому анализу эмпирической яркости фасадов
К теоретическому анализу эмпирической яркости фасадов Печать E-mail
Автор Д.В.Бахарев, И.А.Зимнович (Светотехника, 2008, №3, с.10-17)   
22:07:2008 г.

Недавняя дискуссия (см. "Светотехника", 2006, №№1-3), при всем разнообразии мнений, выразила, однако, общую неудовлетворенность состоянием отечественной методики расчета ЕО. В ее заключение [1] редколлегия журнала рекомендовала теоретически проанализировать учитывающие отраженную составляющую КЕО эмпирические коэффициенты СП 23-102-2003 современными методами компьютерного моделирования.

Достоверность расчета КЕО прежде всего определяется яркостью светящих объектов внешней среды (ОВС), составляющих источник ЕО помещений. Поэтому в первую очередь проанализируем среднюю эмпирическую яркость фасада b_{\small{\cyr f}} противостоящего здания для основной расчетной схемы №1 СП 23-102-2003, на основании которой определяются яркости b_{\small{\cyr f{j}} для всех остальных расчетных схем застройки и, видимо, значения эмпирических коэффициентов K_{\small{\cyr ZD}} и r_0.

В эмпирическую методику расчета ЕО учет яркости фасада был впервые введен в 1962 г в СНиП П-А 8-62 в виде 1/10 значения КЕО в расчетной точке помещения от участка неба, закрываемого противостоящим зданием. В СНиП П-А 8-72 он был дифференцирован коэффициентом K, который согласно табл.6 СНиП возрастал от 0,12 до 0,25 при увеличении коэффициента отражения отделочного материала фасада от 25 до 60%. В СНиП П-4-79 коэффициент K заменили коэффициентом R, косвенно учитывающим зависимости b_\small{\cyr f} от геометрических параметров и расстояния между противостоящими зданиями с помощью индексов z_{\small 1} в плане и z_{\small 2} в разрезе расчетной схемы застройки и помещения. В явном виде табличная функция b_{\small{\cyr f}}\left( \rho_{\small{\cyr f}}, a/H_{\small p}, l/a \right), где \rho_{\small{\cyr f}} - средневзвешенный коэффициент отражения фасада, a/H_{\small p} - отношение длины a противостоящего здания к его расчетной высоте, l/a - отношение расстояния между зданиями к длине a, появилась в МГСН 2.06-99. В табл. В.2 СП 23-102-2003 она представлена в более подробной градации \rho_{\small{\cyr f}}.

В настоящее время проверить компьютерными методами эти результаты полувекового пути, пройденного развитием эмпирической методики расчета ЕО, не представляет затруднений. Почти все программы компьютерной графики содержат сейчас стандартные модели неба МКО при различном состоянии облачности, включая принятое в СНиП облачное небо [2], и довольно точно рассчитывают яркость ОВС любой конфигурации при диффузном, зеркальном и даже двунаправленном отражении. Излагаемые ниже результаты анализа b_{\small{\cyr f{j}} получены при диффузном отражении методом статистического моделирования в программе Radiance 3.8. Программа имеет открытый доступ к исходному коду, подробное теоретическое обоснование и описание [3] и признана наиболее достоверным инструментом компьютерного моделирования и визуализации световых полей.

Image
Рис.1. Зависимости эмпирической и теоретической яркости фасадов под небом МКО от разрыва между зданиями при разных коэффициентах отражения и пропорциях фасадов и процентные расхождения между ними

На рис.1 дано сравнение данных табл.В.2 с расчетами, выполненными Radiance по схеме №1 для двух одинаковых параллельно стоящих зданий при коэффициенте отражения поверхности земли \rho_{\small{\cyr f}} = 0.2. Как видно на рисунке при \rho_{\small{\cyr f}} = 0,2 - 0,6 эмпирическая яркость фасадов оказалась до 40-48% выше теоретической. На графиках процентных расхождений не показаны уходящие в бесконечность завышения b_{\small{\cyr f{j}}} в табличном интервале l/a = "0,25 и менее". Ясно, что при сближении зданий яркость не может оставаться постоянной, а устремляется к 0. Поэтому для указанного интервала следовало особо оговорить необходимость интерполяций яркости между b_{\small{\cyr f{j}}} = 0 при l/a = 0 и табличным значением b_{\small{\cyr f{j}}} при l/a = 0,25. Отсутствием этого разъяснения к табл.В.2 можно, видимо, объяснить обнаруженное в [4] трехкратное завышение КЕО программой Солярис при разрыве 0,5Hp. Для принятых в [4] параметров зданий (a = 62 м, a/H_{\small p} = 16 м, a/H_{\small p} = 3.86, l/a = 0.13) искомое значение b_{\small{\cyr f{j}} находится в середине интервала l/a = "0,25 и менее".

Исправленные на рис.1 графики расхождений эмпирической и теоретической яркости, где ошибки в интервале l/a = 0 - 0,25 стремятся к нулю, имеют довольно упорядоченный экстремальный характер с заметно выраженными зависимостями от \rho_{\small{\cyr f}}, a/H_{\small p} и l/a. Для их анализа необходимо знать, каким образом были получены данные табл.В.2. К сожалению, обнаружить в литературе соответствующие сведения об экспериментальной установке, моделях зданий и методах и средствах измерения и преобразования их результатов в значения b_{\small{\cyr f{j}}} табл.В.2, не удалось. Пришлось предположить, что эти значения получены моделированием под 4,5-метровым куполом НИИСФ, и довольствоваться его описанием и фотографиями полувековой давности [5].

Image
Рис.2. Ортогональные проекции компьютерной модели купола НИИСФ и облачного неба МКО в расчетной точке M фасада здания (слева) и структурированные изображения полей яркости купола и земного диска (справа)

По данным [5] в Radiance была разработана статистическая модель сферического купола 1950 г, замкнутого в уровне центра сферы сооруженным позднее "земным" диском с кольцевой щелью и регулируемыми экранами (рис.2). Параметры диска, щели, экранов и современные характеристики 24 прожекторов, зеркальных ламп и люминесцентных светильников в литературе отсутствуют. В компьютерной модели установлены 24 светильника с асимметричной КСС, создающих результирующую яркость купола, близкую к распределению облачного неба МКО. Асимметричная КСС исключает засветку моделей светильниками в пределах конуса GZ с "земным" диском в основании и вершиной в "зените" купола. Согласно [5] коэффициент отражения купола принят равным 0,8, "земли" - 0,2. В компьютерной модели кольцевая щель для светильников шириною 0,5 м не отражает свет.

На рис.2. показаны синтезированные Radiance гладкие и структурированные ортогональные изображения купола и диска. Здесь отчетливо видна неравномерность яркости купола на стыках светильников у основания полусферы. В вертикальном разрезе купола расхождение с табл.В.1 коэффициента q не превышают 5%. Эти расхождения не существенны, поскольку яркость купола совпадает с яркостью неба только в центре О полусферы. Кажущаяся, видимая нами небесная полусфера в действительности является сегментом шара с радиусом, превышающим земной (6371 км) на высоту нижней границы облачного слоя в 0,1 - 2 км. Высота (стрела) сегмента зависит от типа облаков, формирующих облачный слой, а его радиус (видимый горизонт) - от высоты наблюдения. Для глаза наблюдателя, расположенного на высоте 1,6 - 25 м от поверхности Земли, он составляет 5 - 20 км. При меньшей метеорологической дальности видимости горизонт исчезает в атмосферной дымке. В городе он обычно закрыт окружающей застройкой.

Отсутствие геометрического подобия купола земной атмосфере влечет за собою нарушение фотометрического подобия. Как видно на рис.2, видимая из точки M яркость неба не зависит от ее положения, в то время как видимая яркость купола меняется тем больше, чем больше точка M смещается от центра О купола. Например, в секторе NML она определятся светимостью участка NQL, почти половина которого располагается ниже горизонта, где в натурных условиях находится земля. Яркость земной поверхности под открытым облачным небом МКО зависит только от ее коэффициента отражения. Яркость "земного" диска в основании купола радиально снижается от 0,21 в центре до 0,15 на краю. Создать равномерную освещенность диска светящим куполом, имитирующим яркость облачного неба, видимо, принципиально невозможно [6].

Существенно отличается также 80%-ное отражение поверхности купола от альбедо облачного покрова. Яркость облачного слоя формируется в его толщи в зависимости от поглощения и рассеяния света капельным и кристаллическим водным и пылевым аэрозолем. Альбедо нижней границы облаков обычно не превышает 0,5 [7]. Стандартное небо МКО, установленное для среднего альбедо города, не предполагает его взаимодействия с отдельными зданиями в городской застройке. В модельном эксперименте исключить здания из процесса многократных отражений света и, следовательно, взаимодействия с куполом невозможно.

Image
Рис.3. Синтезированные Radiance фотографии модели зданий под небом МКО (слева) и под куполом НИИСФ (справа) и структурированные поля яркости неба, купола, земли и фасадов зданий (внизу)

На рис.3 показаны синтезированные Radiance фотографии модели двух параллельных зданий под куполом и небом МКО, сделанные с высоты человека (1,7 м), стоящего на "земном" диске в 3-х метрах от его центра, в котором пересекаются диагонали прямоугольного участка земли между зданиями. Они наглядно иллюстрируют отмеченные нарушения подобия экспериментальной установки натуре и порождаемые этим расхождения в яркости фасадов зданий. Представленное на рис.4. сравнение модельной и натурной средней яркости фасадов показывает, что эти расхождения не очень велики. Даже при сравнительно крупном масштабе модели зданий (высота зданий - 1 м, ширина - 0, 15 м) они не превышают 15%. Качественный характер расхождений довольно очевиден.

Image
Рис.4. Зависимости теоретической яркости фасадов под небом и под куполом от разрыва между зданиями при разных коэффициентах отражения и пропорциях фасадов и процентные расхождения между ними

Как видно на рис.4, расхождения возрастают по мере увеличения разрыва между зданиями и их длины. При этом различие между b_{\small{\cyr f{j}}} высоких и протяженных фасадов под небом уменьшается скорее, чем под куполом, поскольку при открытом горизонте натурное значение b_{\small{\cyr f{j}}} мало зависит от пропорций фасада a/H_{\small p} и в основном определяется \rho_\small{\cyr f} и \rho_\small{\cyr z}. Тогда точное значение прямой составляющей КЕО фасада от неба МКО составляет 39,6% [8], а отражение от земли с \rho_\small{\cyr z} = 0,2 добавляет к нему не более 10%. Поэтому при открытом горизонте и \rho_\small{\cyr f} = 0,2; 0,4 и 0,6, b_\small{\cyr f{j}} не может превышать 0,1; 0,2 и 0,3. Соответственно, Radiance дала под небом b_\small{\cyr f{j}} = 0,096; 0,193 и 0,289, а под куполом - 0,106; 0,214 и 0.321, т.е. на 11% выше. При неизменном масштабе модели проследить в эксперименте динамику расхождений при увеличении расстояния между зданиями невозможно, поскольку модель выходит за пределы диска "земли" и "врезается" в купол. На рис.4. этому соответствует случай a/H_{\small p} = 4 при l/a > 1. Здесь штриховой линией показано предполагаемое продолжение графика b_\small{\cyr f{j}}., и анализировать вызванное полным нарушением подобия резкое снижение расхождений не имеет смысла.

Image
Рис.5. Поля яркости фасада с пропорциями a/H_p = 0,25 под небом (левая симметричная половина фасада) и куполом (правая половина) при различных относительных расстояниях l/a(l/H_p) между зданиями и \rho_{\cyr f} = 0,4
 
Image
Рис.6. Поля яркости фасада с пропорциями a/H_p = 1,0 под небом (левая симметричная половина фасада) и куполом (правая половина) при различных относительных расстояниях l/a(l/H_p) между зданиями и \rho_{\cyr f} = 0,4
 
Image
Рис.7. Поля яркости фасада с пропорциями a/H_p = 4,0 под небом (левая симметричная половина фасада) и куполом (правая половина) при различных относительных расстояниях l/a(l/H_p) между зданиями и \rho_{\cyr f} = 0,4

Характер расхождений модели с натурой наглядно иллюстрирует приведенное на рис. 5-7 сравнение симметричных полей результирующей яркости фасадов. Здесь видно, как увеличивается засветка фасадов расположенной ниже горизонта зоной купола по мере приближения к боковым кромкам и, в особенности, к верхней кромке фасада при разных пропорциях и расстоянии между фасадами. При a/H_\small{p} = 4 и l/a > 1 (см. рис.7) поле яркости фасада под куполом качественно искажается за счет входа модели зданий в зону подсветки купола и соответствующего затенения купола и диска "земли".

Image
Рис.8. Зависимости эмпирической и теоретической яркости фасадов под куполом от разрыва между зданиями при разных коэффициентах отражения и пропорциях фасадов и процентные расхождения между ними

Как видно на рис.8, расхождения табличных значений b_\small{\cyr f} с компьютерной моделью купола при l/a = 0,25 и 2,0 достигают 26%, т.е. в 1,5 раза превышают расхождения b_\small{\cyr f{j}} между куполом и небом. Причины этих расхождений очевидны на графиках рис.8. При l/a = 0,25 и a/H_{\small p} = 0,25 и 0,5, а для \rho_\small{\cyr f} = 0,2, даже при a/H_{\small p} = 1, эмпирические значения b_\small{\cyr f{j}} совпадают, чего заведомо не может быть в натуре, поскольку с уменьшением расстояния между фасадами b_\small{\cyr f{j}} должно обязательно снижаться. Здесь 26%-ные расхождения вызваны, видимо, очень грубым округлением до сотых долей и без того малых значений b_\small{\cyr f{j}}. Ясно также, что в интервале l/a > 1 графики b_\small{\cyr f{j}} для разных a/H_\small{p} не могут быть параллельными. При увеличении расстояния между фасадами они должны сходиться к предельным значениям b_\small{\cyr f{j}} при открытом горизонте, которые, как было отмечено выше, мало зависят от a/H_\small{p} и в основном определяются \rho_\small{\cyr f} и \rho_\small{\cyr z}.

На рис.1, 4 и 8 графики b_\small{\cyr f{j}} и процентных расхождений построены в принятых в табл. В.2 "двуотносительных" координатах a/H_{\small p} и l/a. Это позволяет компактно построить таблицы и графики, но зрительно искажает действительные зависимости b_\small{\cyr f{j}} от пропорций и расстояний между зданиями. Как видно на графиках рис.9, построенных в обычных относительных координатах l /H_\small{p} и a/H_\small{p}, интервалы l/a в действительности смещаются по натуральной, сильно вытянутой до значения в 8 высот здания оси l /H_\small{p}. В градостроительстве разрывы более 4H_\small{p} встречаются редко. Практический интерес представляют разрывы ≤2H_{\small{p}. В этом интервале расхождения табличных значений b_\small{\cyr f{j}} с компьютерной моделью купола не выходят за пределы ±5-20%, т.е. обычной российской точности измерений прошлого века [2, стр.52].

Результаты сравнения подтвердили принятую для анализа предпосылку, что табл.В.2 была составлена на основе измерений яркости под куполом НИИСФ на моделях одинаковых параллельно противостоящих зданий или, быть может, имитирующих их фасады плоских щитов, увеличивающих расхождения с натурой. Значительные, превышающие 30% расхождения с яркостью фасадов под небом МКО в интервале разрывов ≤2H_\small{p} следует, видимо, объяснить тем, что в табличные значения b_\small{\cyr f{j}} не были внесены поправки на систематическое, доходящее до 15% завышение яркости, порождаемое нарушением подобия экспериментальной установки натуре.

Image
Рис.9. Натуральные зависимости эмпирической и теоретической яркости фасадов под куполом от разрыва l/H_p между зданиями при \rho_{\cyr f} = 0,4 и разных пропорциях фасадов

Как уже отмечалось в [9], задачи расчета ЕО не выходят за рамки классической теории светового поля. Теоретическая проверка табл.В.2 еще раз подтвердила, что примитивное экспериментирование под куполом не может дать ничего другого, кроме разнообразных ошибок и недоразумений. Не будем здесь повторять замечания [10] по методике преобразования b_\small{\cyr f{j}} из табл.В.2 в яркость фасадов для расчетных схем №2-6, на которые не был дан ответ в дискуссии. Ограничимся лишь демонстрацией ее непригодности для определения b_\small{\cyr f{j}} даже по основной расчетной схеме №1.

Image
Рис.10. Поля яркости земли (слева) и протяженного и торцевого фасадов взаимно перпендикулярных зданий (справа) при их действительном положении (вверху) и ошибочной замене параллельными зданиями согласно указаниям СП 23-102-2003 (внизу). Стрелками показана правильная замена зданий

Рассмотрим два одинаковых взаимно перпендикулярных здания с относительными размерами 1х1х0,25 и \rho_\small{\cyr f} = 0,4, освещаемых небом МКО. Рассчитанные и визуализированные Radiance поля яркости земли и параллельных фасадов зданий показаны на рис.10. При проверке соответствия взаиморасположения зданий нормативным требованиям по ЕО любое из зданий "А" и "В" может стать расчетным, а расчетное помещение может занимать любую позицию на их фасадах.

По методике СП 23-102-2003 расчетное здание на схеме №1 заменяется расчетным помещением (см.рис.В.1) и яркость видимого из расчетной точки через окно участка параллельного помещению фасада противостоящего здания определяется по табл.В.2 в зависимости от \rho_\small{\cyr f}, пропорций a/H_\small{p} и расстояния l/a до этого здания. Если окно расчетного помещения расположено в торце здания "В", то для противостоящего здания "A" a^A/H^A_\small{p} = 1, l^A/a^A = 0,25 и согласно табл.В.2 b^A_\small{\cyr f} = 0,1. Если расчетное помещение находится в здании "А", то для здания "В" a^B/H^B_\small{p} = 0,25, l^B/a^B = 1 и по табл.В.2 b^B_\small{\cyr f} = 0,15. Как видно на рис.10, в действительности, наоборот, b^A_\small{\cyr f} = 0,15, т.е. 1,5 раза выше, а b^B_\small{\cyr f} = 0,1, в 1,5 раза ниже найденных по табл.В.2 значений.

Согласно оптической теории светового поля яркость поверхностей ОВС является размытым изображением неба, проецируемого сложной диафрагмой контура затеняющих его зданий на поверхности земли и фасадов. На земле оно зрительно воспринимается "теневой ямой", образованной сложением рассеянных и размытых теневых изображений зданий. На рис.10 хорошо видны следы плоскостей пространственного рисунка изображений, в которых происходит излом изолиний поля [11]. Яркость протяженного фасада здания "А" воспринимается как теневое изображение здания "В" с четким рисунком, разделяющим изображение на зону торца и зоны ребер, в которых складываются рассеянные изображения торца и протяженных фасадов здания "В". Рассеянное изображение светящего объекта несет более полную информацию о нем, чем четкое зеркальное изображение. В четком изображении здания "В", близком по структуре к фронтальной перспективе, его протяженные фасады не могут быть видимыми. Поэтому нельзя заменять фасады здания плоскими щитами или "условными затеняющими плоскостями", как это предписывается в СП 23-102-2003. На торце здания "В", как в осколке зеркала, умещается только часть пространственного теневого изображения протяженного фасада здания "А". Четкое воспроизведение на рис.10 структуры полей свидетельствует о высокой точности статистической оценки яркости программой Radiance. Средние яркости первого отражения на фасадах совпали с вычисленными программой Lara детерминированным методом значениями с точностью 10-3.

Таким образом, поля яркости фасадов противостоящих зданий есть их взаимные оптические изображения. Поэтому при одинаковых значениях \rho_\small{\cyr f} противостоящий диффузный отражатель можно заменить отражением объекта, а при разных \rho_\small{\cyr f} его среднюю яркость скорректировать пропорционально отношению \rho_\small{\cyr f} отражателя к среднему \rho_\small{\cyr f} обоих зданий. Здесь очевидна аналогия с зеркалом, средняя яркость которого определяется яркостью фона и глядящего в зеркало субъекта. Однако это не дано видеть слепому, пренебрегающему оптической (зримой) теорией, эмпиризму. На рис.10 видно, что замена противостоящих зданий расчетными дает на них искомые значения b_\small{\cyr f{j}} всего на 2-4% меньше рассчитанных Radiance.

Возлагать, однако, большие надежды на замену отражателя отражаемым изображением не следует, поскольку она справедлива только для схем с двумя зданиями. В реальной застройке яркость фасада определяется обычно изображением нескольких зданий, заменить которое изображением одного расчетного здания уже невозможно. Кроме того, в плотной застройке значения КЕО, рассчитанные по среднему значению яркости фасада b_\small{\cyr f{j}}, является заведомо завышенными, поскольку яркость светящих через окно в расчетную точку участков фасадов, выделенных на рис.10 прямоугольным контуром, на нижних этажах значительно меньше b_\small{\cyr f{j}}.

Ошибки, порождаемые слепой выборкой из табл.В.2 значений b_\small{\cyr f{j}}, зависят от соотношения размеров противостоящего и расчетного зданий и могут достигать ±2,5-кратной величины. Занижения b_\small{\cyr f{j}}, видимо, компенсируются эмпирическими поправками K_\small{\cyr{ZD}j} и K_\small{\cyr{ZD}0}, "учитывающими изменение внутренней отраженной составляющей КЕО в помещении при наличии противостоящих зданий" и "при полном закрытии небосвода зданиями, видимыми из расчетной точки" и определяемыми по формуле (В.5) и по табл.В.6 СП 23-102-2003. Согласно табл.В.6 для расчетных точек в середине и в глубине помещений значения K_\small{\cyr{ZD}0} = 1,0 - 2,07, а из (В.5) следует, что K_\small{\cyr{ZD}j} всегда больше 1. Дальнейший совместный анализ K_\small{\cyr{ZD}0} и r_0 покажет, в какой мере множитель K_\small{\cyr{ZD}j} > 1, способен компенсировать завышения b_\small{\cyr f{j}} и к каким последствиям это приводит в расчете КЕО по эмпирической методике СП 23-102-2003.

Задачу расчета поля яркости ОВС в застройке с достаточной для практических целей точностью можно считать инвариантной к расчету КЕО в помещениях. Если это поле не рассчитано предварительно, то едва ли возможен достоверный расчет ЕО помещений. Как было продемонстрировано выше, компьютерные расчеты полей яркости ОВС не представляют сейчас затруднений. Остается только сожалеть, что бюджетные средства впустую уходят на реанимацию безнадежно устаревшей ручной методики расчета КЕО, а не на разработку отечественного программного обеспечения инженерных расчетов ЕО. Это положение необходимо срочно изменить, иначе отставание российской строительной светотехники от мирового уровня и запросов отечественной практики, давно перешедшей на компьютерные методы проектирования, станет необратимым.

Литература

  1. Заключение о дискуссии "О нормах и методах расчета естественного освещения и инсоляции в зданиях" // Светотехника, 2006, №3 , с.65-66.
  2. Справочная книга по светотехнике, 3 изд. Москва: "Знак". 2006. стр. 972 (см. раздел 8, параграф 12, Компьютерные методы расчета ЕО).
  3. Ward Larson, G., R. Shakespeare, Rendering with Radiance: The Art and Science of Lighting Visualization. San Francisco. CA: Morgan Kaufman Publishers. 1998. р. 644.
  4. Зимнович И.А. Выступление в дискуссии по инсоляции и ЕО // Светотехника, 2006, №3 , с.64.
  5. Гусев Н.М., Естественное освещение зданий, М.: Госстройиздат, 1961, с.170.
  6. Бахарев Д.В. Ларисы отражающих систем//Светотехника, 1994. № 1. с. 15 - 19.
  7. Фейгельсон Е.М., Радиационные процессы в слоистообразных облаках, Институт физики атмосферы АН СССР, М.: Изд."Наука". 1964. с.230.
  8. Бахарев Д.В. Освещенность от неба конхоидальной яркости // Светотехника. 1989. № 9. с. 12 - 13.
  9. Бахарев Д.В., Зимнович И.А. О светопропускании окон // Светотехника. 2007. №5. с. 4 - 9.
  10. Бахарев Д.В. О методике расчета естественного освещения // Светотехника. 2006. № 1. с. 57 - 59.
  11. Бахарев Д.В. О структуре световых полей // Светотехника. 2005. № 3. с.40-44.
Обсудить новость в форуме. (0 сообщений)
 
« Николай Алексеевич Рынин   О визуализации спектральной модели безоблачного неба и солнца »