Алгоритм расчёта отражений на основе геометрической модели

УДК 699.814

АЛГОРИТМ РАСЧЕТА ОТРАЖЕНИИ НА ОСНОВЕ ГЕОМЕТРИЧЕСКОЕ МОДЕЛИ

Н.В. Благородова, А.А. Замятин, В.В. Сухомлинова

Ростовский государственный строительный университет

Размер минимальных противопожарных расстояний между зданиями городской и промышленной застройки, а также безопасные разрывы между технологическим оборудованием, излучающим тепло, назначаются действующими нормами проектирования в области пожарной и промышленной безопасности. Рекомендуемые нормативные расстояния справедливы на 100% для расположенных строго параллельно элементов простой геометрической формы. При определении рекомендуемых нормами размеров минимальных противопожарных расстояний между зданиями используется классическая теория передачи тепла излучением. При этом учитывается только действие лучистой энергии. Распространение пожара от теплоты, передаваемой конвекцией, от горящих конструктивных элементов и искр, перебрасываемых ветром и нагретым воздухом, не учитывается. Это объясняется тем, что количество теплоты, передаваемое конвекцией, составляет небольшой процент по сравнению с теплотой, передаваемой излучением.

В расчетах лучистой энергии, на наш взгляд, значительный вклад в общую картину вносят первые отражения. На практике, при решении задач, связанных с расчётами лучистой энергии, целесообразно применять геометрическую модель. Эта модель даёт хорошие результаты в тех случаях, когда можно пренебречь волновыми свойствами излучения. В данной модели не рассматриваются явления интерференции, дифракции и др., связанные с волновыми свойствами излучения [1].

Как известно, падающий и отражённый лучи лежат в одной плоскости с нормалью к поверхности, построенной в точке падения, и составляют с ней равные углы. Рассмотрим алгоритм построения отражённого луча при данных источнике и точке падения на отражающей поверхности.

Пусть задана поверхность £ уравнением Г = Г (и, V) . В пределах рассматриваемого отсека поверхность не имеет особых точек. На поверхности задана точка А, её радиус-вектор равен ЯА = Г(иА, VA), задан источник излучения, находящийся в точке I (рис.1). Вектор нормали к поверхности £ в точке А определим по формуле

дГ(иА,VA) дГ(иА,VAУ

N =

ди Зу

Как уже отмечалось, падающий, отражённый лучи и нормаль к поверхности лежат в одной плоскости. Обозначим эту плоскость через Пі (рис. 1) и будем называть лучевой плоскостью. Вектор нормали лучевой плоскости

N1 =

Сделаем такое преобразование координат, чтобы лучевая плоскость Пі в новой системе координат стала фронтальной плоскостью уровня [2].

Рис. 1. Геометрическая модель отражения

Определим угол, составляемый горизонтальной проекцией вектора нормали П/ с осью Oy - угол у:

cos ax

Y = arctg----------,

cos a

y

где cos ax и cos ay направляющие косинусы вектора Ni по осям Ox и Oy, равные

к

cos a

cos a =

лік2 +12 + m 2 l

У лk2 +12 + m 2 ’

к, l, m - координаты вектора Ni.

Определим угол, составляемый Ni с плоскостью xOy:

S = n-az,

2 z

где az - угол, составляемый Nl с осью Oz, равный az = arccos

m

лік2 +12 + m 2

Повернем систему координат вокруг оси Оі на угол у (рис.2). В новой системе координат вектор нормали плоскости П/ будет лежать в координатной плоскости уОі. Матрица данного преобразования имеет вид:

ґ---- —~ п\

А(у) =

Повернем систему координат вокруг оси Ox на угол 0 (рис.2).

cos у sin Y 0

- sin у cos Y 0

0 0 1

Рис. 2. Углы, составляемые вектором нормали лучевой плоскости с осями и плоскостями проекций

В новой системе координат вектор нормали П/ параллелен оси Оу. Следовательно, в этой системе координат П/ является фронтальной плоскостью уровня. Матрица этого преобразования имеет вид:

'10 0 Л

ад

0 cos 3 - sin 3 0 sin 3 cos З

V " ~ ~ /

Получим матрицу суммарного преобразования, умножив Л

^ cos у sin у 0

Г

\

- sin y cos 3 cosacos 3 - sin 3

- sin y - sin 3 cos y sin 3 cos 3

Найдем координаты точек I, А и вектора N в системе координат, в которой П/ является плоскостью уровня. Координаты найдем по соотношению

где х, у, г - координаты соответствующей точки в исходной системе координат, х’, у’, г’ - ее координаты в новой системе координат. Т.к. П/, после преобразования, является фронтальной плоскости уровня, то все ее точки имеют у ’=сотї и трехмерная задача сводится к двумерной.

Сдвинем начало системы координат в точку А. Координаты точек определим по соотношению

Где х’А, 7’а координаты точки А. Повернем систему координат таким образом, чтобы ось

О”^’” совпала с вектором N.

где С

cos ф - sin ф у sin ф cos ф J

- матрица преобразования, ф = arccos

z

N

ml + _т!

NN

x’’’N , z’’’N - координаты вектора N.

/ В

V \ А т \ X \ L /

А

Рис. 3. Построение отраженного луча

В новой системе координат найдем координаты точки В, которая будет определять, совместно с точкой А, отраженный луч. Её координаты равны

X

B

-X , Ув= Уі , z

fff fff /1 1 \

b= zi •

Выполним обратные преобразования, чтобы определить координаты точки B в исходной системе координат

(хВ ZB) = (xB zB)хC(-ф); (і2)(xB ZB) = (x'B ZB)+(xA za);

(

X

в

Ув Z

в

)= (xB

Ув

zB )хГ-1(^,$),

Г-1

где Г" (у, 3)- матрица обратного преобразования по отношению к (7), равная

'eos у sin y eos 3 sin y sin 3Л

sin y eos y eos 3 eos y sin 3

у 0 - sin 3 eos 3 j

Программный алгоритм построения отраженного луча имеет линейную структуру, которая заключается в последовательности выполнении вычислений и преобразований, описанных выше. Программное обеспечение разрабатывалось в среде ObjectARX в среде AutoCAD [3]. Отражающие поверхности аппроксимировались методом Фергюсона [4,5]. Пример вычислений приведен на рис. 4. На рис. 4 показаны отраженные от поверхности Q лучи, исходящие из источника, расположенного в точке I.

Данный алгоритм может быть использован в качестве модуля различных автоматизированных систем, для расчетов вклада первых отраженных лучей в суммарную энергию излучения, в том случае, если допускается применение геометрической модели для описания процессов отражения.

Q

Рис. 4. Пример расчета

Библиографический список

1. Ландсберг Г.С. Оптика / Г.С. Ландсберг. - Москва: Наука, 2003. - 848 с.

2. Замятин А.В. Алгоритм построения точек пересечения нелинейчатых поверхностей. Электронный журнал №3 2010. Электрон. ресурс. Режим доступа: http://www.ivdon.ru/ (дата обращения 29 апреля 2012).

3. Полищук Н.Н. AutoCAD: разработка приложений, настройка и адаптация / Н.Н. Полищук. -Санкт-Петербург: БХВ - Петербург, 2006. - 992 с.

4. Роджерс Д., Адамс Дж. Математические основы машинной графики / Д. Роджерс, Дж. Адамс. - Москва: Мир, 2001. - 604 с.

5. 5. Замятин А.В., Сухомлинова В.В. Аппроксимация порции поверхности по методу Фергюсона / Известия вузов. Северо-Кавказский регион. Технические науки. Приложение №2, 2006г.