CC BY
59
РАСЧЁТ ОПТИЧЕСКОГО ЭЛЕМЕНТА ДЛЯ ФОРМИРОВАНИЯ ОСЕСИММЕТРИЧНОГО РАСПРЕДЕЛЕНИЯ ОСВЕЩЁННОСТИ НА ПОВЕРХНОСТИ ВРАЩЕНИЯ
Досколович Л Л.1'2, Асланов Э.Р.2 1 Институт систем обработки изображений РАН 2 Самарский государственный аэрокосмический университет имени академика С.П. Королёва, (национальный исследовательский университет)
Аннотация
Разработан метод расчёта осесимметричного оптического элемента для получения заданного распределения освещённости на поверхности вращения. Расчёт поверхностей оптического элемента сведён к интегрированию системы трёх дифференциальных уравнений первого порядка, разрешённых относительно производной. В качестве примеров рассчитаны оптические элементы для формирования равномерного распределения освещённости на сферической поверхности и вдоль отрезка оптической оси.
Ключевые слова: геометрическая оптика, оптический элемент, распределение освещённости, светодиод, светотехника.
Введение
Светотехнические устройства являются важным элементом большого числа технических систем, включающих дорожное, жилое, промышленное освещение, светотехнические системы транспортных средств и т.д. До недавнего времени основными используемыми источниками света являлись лампы накаливания. В настоящее время эволюция светотехнических устройств направлена на переход от ламп к светодиодам. Использование светодиодов в системах освещения требует применения вторичной оптики, направляющей излучённый светодиодом световой поток в заданную область пространства и обеспечивающей формирование в этой области требуемого распределения освещённости. Расчёт вторичной оптики является сложной задачей и даже в случае точечного источника требует решения нелинейного дифференциального уравнения в частных производных второго порядка [1].
Для получения заданного распределения освещённости в области с малым угловым размером (50 - 60° и менее) обычно применяют оптические элементы, использующие эффект полного внутреннего отражения (ПВО) для поворота лучей на большие углы. Такие элементы, как правило, имеют две преломляющие поверхности (внутренняя и внешняя поверхности а, с1 на рис. 1) и одну поверхность, работающую на ПВО (боковая поверхность с на рис. 1). Наличие нескольких рабочих поверхностей в общем случае требует использования сложных оптимизационных и итерационных методов [2 - 7].
В случае точечного (компактного) источника излучения и осесимметричного распределения освещённости возможно аналитическое решение задачи расчёта оптического элемента. В частности, в работах [8, 9] расчёт осесимметричных оптических элементов светодиодов, содержащих поверхность, работающую по принципу ПВО, сведён к численному интегрированию двух дифференциальных уравнений первого порядка, разрешённых относительно производной. Оптические элементы в [8, 9] предназначены для формирования заданного осесиммет-
ричного распределения освещённости в плоскости, перпендикулярной оптической оси.
В данной работе рассмотрена задача расчёта оптического элемента, обеспечивающего формирование заданного распределения освещённости на поверхности, образованной вращением кривой вокруг оптической оси линзы. Указанная задача обобщает результаты работ [8, 9] и включает расчёт оптического элемента для формирования заданного распределения освещённости в плоскости в качестве частного случая. Задача сведена к интегрированию трёх дифференциальных уравнений первого порядка, разрешённых относительно производных. В качестве важного частного случая рассмотрен расчёт оптического элемента для фокусировки в отрезок оптической оси.
1. Расчёт формы профиля линзы
Рассмотрим задачу расчёта оптического элемента для формирования требуемого распределения освещённости на поверхности вращения с образующей р = p(z), zе [zb z2] (рис. 1). Функцию p(z) будем считать монотонной (возрастающей или убывающей), чтобы поверхность не содержала областей геометрической тени.
Считаем, что точечный источник находится в начале координат. Для создания требуемого распределения освещённости на поверхности вращения будем использовать осесимметричный преломляющий оптический элемент с профилем, приведённым на рис. 1. Профиль элемента состоит из частей а, b, с и внешней плоской поверхности d (рис. 1). Рабочими поверхностями элемента являются поверхности a и с. При этом поверхность a является преломляющей, а поверхность с работает на эффекте полного внутреннего отражения (ПВО). Поверхность b является сферической с центром в начале координат, поэтому преломления на ней не происходит.
Рассмотрим расчёт части a профиля элемента (рис. 1). Обозначим ra ф0) длину радиус--вектора точки профиля, где р0 е [0, р0 max] - угловая координата падающего луча (рис. 2).
О р
Рис. 1. Профиль ПВО оптики с верхней плоской поверхностью
В работе [8] для га (р0) получено следующее дифференциальное уравнение:
dra (Р0) dP„
= -r. (b„)
sin (P„ (P„ ))
1/n - cos (b„ +Pa (b„ ))
(1)
где функция ра ф0) задаёт направление преломлённого луча (рис. 2), п - показатель преломления материала элемента.
Рис. 2. Преломление луча на поверхности а и на внешней плоской поверхности d
Для расчёта функции pa (Р0) в (1) используем закон сохранения светового потока. Световой поток, падающий на элемент dS поверхности вращения с образующей p(z), должен быть равен соответствующему световому потоку от источника, излучённому в элемент телесного угла dW0 = 2p sin р0 dp0. Таким образом, можно записать следующее равенство:
I„ (р„)dW„ = 2pI„ (p„) sinp„dp„ =
= E (z) dS = 2pE (z) p (z )yj 1+(dp( z)/dz )2dz,
(2)
где I0 (Р0) - интенсивность источника излучения, E(z), z e [z1, z2] - заданная освещённость на поверхности вращения. Представим уравнение (2) в следующем виде:
dz(р0 )_ I0 (р0 ) sinр0
dP„
E(z)p(zU1 + (dp(z)/dz)
(3)
Отметим, что корректное задание освещённости Е(г) требует выполнения закона сохранения светового потока:
P„m
f I„ (р„ ) sinp„dp„ =
= ÍE (z )p( z \ k-
dp( z) dz
(4)
dz.
Запишем уравнение выходящего из оптического элемента луча в следующем параметрическом виде:
[p(p„ )=p* (р„)-1 sin pa (p„ ), I z(p„ ) = h +tcos pa (p„),
(5)
где ph (p„) - координата точки пересечения луча, преломлённого на части a профиля, с выходной плоскостью z = h, pa (p„) — угловая координата
луча после преломления на выходной плоскости. Согласно закону Снеллиуса,
sin pa (p„ ) = n • sin pa (p„ ) .
Из геометрии задачи (рис. 2) несложно получить ph (p„) в (5) в виде:
Ph (p„ ) = Pa (p„ )-( h - za (p„ )) tg pa (p„ ) , (6)
где pa (p„) = ra (p„) sin p„, za (p„) = ra ф„) cos p„ - координаты пересечения падающего луча c профилем a. Будем считать, что p ф^ в правой части первого уравнения в (5) соответствует радиальной координате образующей поверхности вращения p = p (z) при z = z фа) , где функция z(p„) является решением уравнения (3). Тогда, исключая из первого уравнения в (5) переменную t и подставляя полученное выражение во второе уравнение, получим:
z(p„) = h + [rb (p„)-p(p„)]x
Ф-(n sin pa (p„ ))2
(7)
n sin
in pa (p„)
Продифференцируем уравнение (7) по переменной р0 и получим следующее уравнение:
^ = (p„ ) Щ cos
dp„ ЛР„> ¡3 (p„)
pa(p„),
(8)
где
¡1 (p„ ) = 1 - [n sin pa ]
¡2 (p„ bffr)
d r
"тр sin (pa +p„) + dp„
+ ra cos (pa +p„ )- d^cos pa dp„ _
¡3 (p„ ) = ( z(p„) - h) ncos3 pa + + (h - ra cos pa ) ( ¡1(p„))3/2.
- nJz. sin2pa, (9) 2dp„ Pa' W
„
2
2
dra dz
Производные ——, - в (9) определены в
dPo dPo
dp dp dz
уравнениях (1), (3), а —— = —--. После выполнено dz dPo
ния указанных подстановок правая часть уравнения (9) будет зависеть только от переменной ро и функций ra (b0), z(b0). Таким образом, расчёт части a профиля оптического элемента сводится к интегрированию системы трёх дифференциальных уравнений (1), (3), (8), разрешённых относительно производных. При этом уравнение (3) является независимым от уравнений (1), (8) и может быть решено отдельно. Указанная система может быть проинтегрирована с использованием стандартных численных методов (например, метода Рунге-Кутты).
Рассмотрим расчёт части с оптического элемента (рис. 1), работающей на ПВО. Обозначим rc(ро), р0 е [р0 max, л/2] длину радиус-вектора точки профиля с. В [8, 9] для rc (р0) получено следующее дифференциальное уравнение:
MoL-rc,^^, (10)
где функция рс (р0) задаёт направления луча после ПВО (рис. 1). Для расчёта части с профиля оптического элемента могут быть также записаны дифференциальные уравнения, которые полностью аналогичны уравнениям (3), (8).
Важным частным случаем рассмотренной задачи является расчёт оптического элемента для фокусировки в отрезок оптической оси. В этом случае закон сохранения светового потока (2) принимает вид:
2pIo (ро) sinМро = E(z) dz,
dz (ро) 2л г (р ) . р
"ар" = ~ЕЩ (ро)ипро, (11)
где E(z) - заданная освещённость на отрезке [zb z2] оптической оси. При этом функция p (ро) ° о в уравнениях (7) - (9).
2. Результаты численного моделирования
Разработанный метод расчёта оптических элементов был реализован в среде программирования Matlab [1о].
Были рассчитаны оптические элементы для равномерного освещения полусферы и отрезка [z1, z2] на оптической оси. Уравнение образующей полусферы имеет вид:
Р(*) = VRz-Z>ze [z0,z0 + R],
(12)
где Я - радиус полусферы, 20 - координата центра сферы на оптической оси. Расчёт оптического элемента для освещения полусферы производился при параметрах: Я = 50 мм, 20 = 215 мм. Расчёт оптического элемента для равномерного освещения отрезка оптической оси производился при 21 = 40 мм, 22 = 90 мм.
15
10
Л
а) 0 2 4 6 8 10
0 2 4 6 8 10 12
Рис. 3. Профиль линзы, обеспечивающей равномерное освещение полусферы (а); отрезка оптической оси (б)
При расчёте обоих элементов были использованы следующие параметры: толщина элемента h = 15 мм, показатель преломления материала элемента n = 1,41 (соответствует силикону), угловой размер части a профиля а р0 max = 30°. Источник излучения считался ламбертовским, при этом I0 (р0) = cos р0 в (3). Рассчитанные профили элементов приведены на рис. 3. Для проверки правильности расчётов было проведено моделирование работы рассчитанных оптических элементов в программе для светотехнических расчётов Zemax с использованием метода трассировки лучей [11]. Распределения освещённости, рассчитанные для 1 000 000 лучей, приведены на рис. 4, 5. При этом выбор лучей осуществлялся по методу Соболя [12], который позволяет достичь большей точности расчёта.
Расчётные распределения освещённости на рис. 4 показывают высокую равномерность распределения освещённости на полусфере. Некоторые осцилляции освещённости вдоль образующей полусферы на рис. 4а объясняются погрешностью метода трассировки лучей. Среднеквадратичное отклонение освещённости от среднего значения на поверхности полусферы составляет 6,5%. При этом энергетическая эффективность элемента (доля излучённого светового потока, попавшая на поверхность полусферы) превышает 89%.
Расчётные распределения освещённости на рис. 5 показывают формирование области равномерной освещённости вдоль отрезка оптической оси.
Нормированная освещённость
1,0-
0,8 0,6 0,4 0,2 О
/V л ¡АЛ
V
а)
50 40 30 20 10 О -10 -20 -30 -40 -50
б) -55 -40 -20 0 20 40 55
Рис. 4. Нормированные распределения освещённости, полученные при моделировании в программе Zemax оптического элемента на рис. 3а для равномерного освещения полусферы: освещённость вдоль образующей полусферы (а), полутоновое распределение освещённости на полусфере (б)
Нормированная освещённость 1,01
а)
-10 0 10
Рис. 5. Нормированные распределения освещённости, полученные при моделировании в программе Zemax оптического элемента на рис. 3б для равномерного освещения отрезка оптической оси: освещённость вдоль оптической оси (а), полутоновое распределение освещённости в плоскости, содержащей оптическую ось (б)
Отметим, что координата z на рис. 5 отсчитывается от выходной плоскости оптического элемента, поэтому отрезок располагается при z е [25, 75] мм. Среднеквадратичное отклонение освещённости от среднего значения в области отрезка составляет 8,3%. Световая эффективность оптического элемента (доля излучённого светового потока, попавшая на плоскость регистрации на рис. 5б) также превышает 89%.
При расчётах мы предполагали, что источник является точечным, но в реальных условиях размеры излучающей области светодиодов начинаются от десятых долей миллиметра. Отметим, что предлагаемые в статье элементы чувствительны к размеру источника, особенно это касается линз, которые формируют равномерное освещение отрезка. Так, при использовании источника с прямоугольной областью излучения размерами 0,1 мм х 0,1 мм среднеквадра-тическое отклонение возрастало с 8,3% до 15%.
Заключение
Задача расчёта преломляющего оптического элемента, предназначенного для формирования заданного осесимметричного распределения освещённости на поверхности вращения, сведена к численному интегрированию трёх дифференциальных уравнений первого порядка, разрешённых относительно производных. Рассчитаны оптические элементы для формирования равномерного распределения освещённости на сферической поверхности и вдоль отрезка оптической оси. Результаты моделирования в программе для светотехнических расчётов Zemax показывают высокую равномерность формируемых распределений освещённости и подтверждают правильность предложенного метода расчёта.
Благодарности Работа выполнена при поддержке грантов РФФИ 13-07-97001, 12-07-13113 и ФЦП «Кадры» (соглашение 8027).
Литература
1. Guan, P. On a Monge-Ampere equation arising in geometric optics / P. Guan and X.-J. Wang // J. Diff. Geom. -1998. - Vol. 48. - P. 205-223.
2. Sun, L. Free-form microlens for illumination applications / L. Sun, Sh. Jin and S. Cen // Appl. Opt. - 2009. - Vol. 48. - P. 5520-5527.
3. Jiang, J. Optical design of a freeform TIR lens for LED streetlight / J. Jiang, S. To, W.B. Lee and B. Cheung // Optik. - 2010. - Vol. 121. - P. 1761-1765.
4. Vazquez-Molini, D. High-efficiency light-emitting diode collimator / D. Vazquez-Molini, M. Gonzalez-Montez, A. Alvarez and E. Bernabeu // Opt. Eng. - 2010. - Vol. 49. - P. 123001.
5. Zhao, Sh. Lens design of LED searchlight of high brightness and distant spot / Sh. Zhao, K. Wang, F. Chen, D. Wu and Sh. Liu // J. Opt. Soc. Am. - 2011. - Vol. A28. -P. 815-820.
6. Grabovickic, D. TIR RXI collimator / D. Grabovickic, P. Benitez, and J.C. Minano // Opt. Express. - 2012. -Vol. 20. - P. A51-A61.
7. Kari, Th. Reliability of point source approximations in compact LED lens designs / Th. Kari, J. Gadegaard,
Th. Sondergaard, Th.G. Pedersen and K. Pedersen // Opt. Express. - 2011. - Vol. 19. - P. A1190-A1195.
8. Doskolovich, L. Design of TIR optics generating prescribed irradiance distribution in the circle region / L. Doskolovich, M. Moiseev // J. Opt. Soc. Am. - 2012. - A 29.
- P. 1758-1763.
9. Моисеев, М.А. Расчёт оптических элементов, работающих по принципу полного внутреннего отражения и формирующих заданные световые распределения в круглых областях / М.А. Моисеев, Л.Л. Досколович, Е.В. Бызов, С.В. Кравченко// Компьютерная оптика. -2012. - Т. 36, № 3. - С. 333-340.
10. http://www.mathworks.com/
11. http://www.radiantzemax.com/en/zemax/
12. http://www.radiantzemax.com/kb-en/KnowledgebaseArtic le50217.aspx
References
1. Guan, P. On a Monge-Ampere equation arising in geometric optics / P. Guan and X.-J. Wang // J. Diff. Geom. -1998. - Vol. 48. - P. 205-223.
2. Sun, L. Free-form microlens for illumination applications / L. Sun, Sh. Jin and S. Cen // Appl. Opt. - 2009. - Vol. 48.
- P. 5520-5527.
3. Jiang, J. Optical design of a freeform TIR lens for LED streetlight / J. Jiang, S. To, W.B. Lee and B. Cheung // Optik. - 2010. - Vol. 121. - P. 1761-1765.
4. Vazquez-Molini, D. High-efficiency light-emitting diode collimator / D. Vazquez-Molini, M. Gonzalez-Montez, A. Al-
varez and E. Bernabeu // Opt. Eng. - 2010. - Vol. 49. -P. 123001.
5. Zhao, Sh. Lens design of LED searchlight of high brightness and distant spot / Sh. Zhao, K. Wang, F. Chen, D. Wu and Sh. Liu // J. Opt. Soc. Am. - 2011. - Vol. A28. -P. 815-820.
6. Grabovickic, D. TIR RXI collimator / D. Grabovickic, P. Benitez, and J.C. Minano // Opt. Express. - 2012. -Vol. 20. - P. A51-A61.
7. Kari, Th. Reliability of point source approximations in compact LED lens designs / Th. Kari, J. Gadegaard, Th. Sondergaard, Th.G. Pedersen and K. Pedersen // Opt. Express. - 2011. - Vol. 19. - P. A1190-A1195.
8. Doskolovich, L. Design of TIR optics generating prescribed irradiance distribution in the circle region / L. Do-skolovich, M. Moiseev // J. Opt. Soc. Am. - 2012. - A 29. - P. 1758-1763.
9. Moiseev, M.A. Design of led tir-optics producing prescribed illuminance distributions in circle regions / M.A. Moiseev, L.L. Doskolovich, E.V. Byzov, S.V. Krav-chenko // Computer Optics. - 2012. - V. 36, N 3. -P. 333-340. - (In Russian).
10. http://www.mathworks.com/
11. http://www.radiantzemax.com/en/zemax/
12. http://www.radiantzemax.com/kb-en/KnowledgebaseArtic le50217.aspx
DESIGN OF THE OPTICAL ELEMENT GENERATING AXISYMMETRIC UNIFORM LIGHT DISTRIBUTION ON THE SURFACE OF REVOLUTION
L.L. Doskolovich1,2, E. Aslanov2 1 Image Processing System Institute of the RAS, 2 S.P. Korolyov Samara State Aerospace University (National Research University)
Abstract
We present the computation method for axisymmetric optical element, which generates required uniform light distribution on the surface of revolution. The computation problem is reduced to integration of three first-order differential equations. We consider two examples: a uniformly illuminated hemisphere and the uniformly illuminated segment of optical axis.
Key words: optics design, light distribution, LED, lighting engineering.
Сведения об авторах
Досколович Леонид Леонидович, 1966 года рождения, в 1989 году с отличием окончил Куйбышевский авиационный институт (КуАИ, ныне - Самарский государственный аэрокосмический университет имени академика С.П. Королёва - СГАУ) по специальности «Прикладная математика». Доктор физико-математических наук (2001 год), профессор, работает ведущим научным сотрудником лаборатории дифракционной оптики Института систем обработки изображений РАН (ИСОИ РАН), профессором кафедры технической кибернетики СГАУ. Специалист в области дифракционной оптики, лазерных информационных технологий, нанофотоники.
Страница в интернете: http://www.ipsi.smr.ru/staff/Doskolovich.htm . E-mail: leonid@smr.ru .
Leonid Leonidovich Doskolovich (b. 1966) graduated with honours (1989) from the S.P. Korolyov Kuibyshev Aviation Institute (presently, Samara State Aerospace University named after S.P. Korolyov (SSAU)), majoring in Applied Mathematics. He received his Doctor in Physics & Maths (2001) degrees from Samara State Aerospace University. Leading researcher of the Image Processing Systems Institute of the RAS, professor at SSAU's Technical Cybernetics sub-department. Current research interests include diffractive optics, laser information technologies, nanophotonics.
Асланов Эмиль, 1986 года рождения. В 2009 году с отличием окончил Московский государственный университет имени Н.Э. Баумана по специальности «Экологический мониторинг». Занимался исследованиями в области светодиодной оптики в компании LG с 2008 по 2012. В настоящее время является экспертом в области патентования в Центре интеллектуальной собственности. Является членом международного научного общества SPIE. Область научных интересов: геометрическая оптика, дифракционная оптика. E-mail: EAslanov@sk.ru.
Emil Aslanov (b. 1986) graduated with honors (2009) from Bauman Moscow State Technical University, majoring in Ecological Monitoring. He has been working in LG as an optical R&D engineer 2008-2012. Now he is a patent expert in the Center of Intellectual Property. He is a SPIE member. His current research interests include geometrical optics, mathematical modeling, diffractive optics.
Поступила в редакцию 9января 2013 г.
