陳巧云 朱向冰 倪 建 陳 瑾
(安徽師范大學物理與電子信息學院,安徽蕪湖241000)
作為第四代新型節(jié)能光源的典型代表,LED具有壽命長、節(jié)能、環(huán)保、顯色性好等其他光源無法比擬的優(yōu)點,而且隨著大功率LED的迅速發(fā)展[1~4],越來越廣泛的應用于各種指示、顯示、背光源、普通照明和城市夜景等眾多領域[5~7]。但是LED光源發(fā)出的光不均勻,因此在對光照有嚴格要求的特殊場合的照明中,必須根據LED本身的發(fā)光特性和預期的光強分布,來設計一個二次配光結構,將LED的光強進行重新分配[8~11]。
配光結構一般有反射面、透鏡兩種,目前LED的二次配光結構大部分采用的是透鏡結構[12~14]。考慮到透鏡有兩個透光面和一定的厚度,對于任意一個曲面,若曲面設計或制作過程中稍有偏差,或透鏡中夾有一點雜質,對光線的折射和能量分布會產生很大的影響,且在實際應用中透鏡對光線有較高的吸收率而造成能量損失。而用反射面作為二次配光結構,理論上只需要一個反射面即可,容易加工且能量損失少,較好地解決了LED均勻照明的問題。本文旨在設計一種反射面來對LED進行二次配光。
傳統(tǒng)的反射面配光結構依賴于有解析解的二次曲面或二次曲面的組合。但對于發(fā)光不均勻的光源,簡單的二次曲面或二次曲面的組合都不能實現(xiàn)特定的光強分布[15~17]。本文采用自由曲面設計的方法[18,19],根據光源的發(fā)光特性和預期的光強分布建立一個微分方程,用數(shù)值求解法求解出這個方程,根據求解結果擬合出自由反射曲面,再利用光學軟件對結果進行模擬。
本文中光源發(fā)出的光線一部分直接照射到照明面上,另一部分經設計的反射面反射到照明面上,最終在照明面上實現(xiàn)一個矩形的均勻照明區(qū)域。利用光的反射定律得到出射光線與入射光線的矢量的關系,再利用能量守恒和照度公式求出反射光線攜帶的能量,最后列出微分方程,通過求解微分方程得出確定自由曲面的一系列點,最終確定所需的自由曲面。
如圖1所示,o′為光源所在的位置,建立坐標系x′o′y′。1為照明面,為了使計算更加方便,設光源o′到照明面1的距離為單位長度1,以o′為圓心,作一個單位球,則單位球與平面1相切于o點。平面2為單位球的平行于平面1的一個切平面。在平面1和平面2內分別作平行于x′,y′的x,y軸和u,v軸。z′軸的原點在o′處。
任一條光線與單位球面相交于p點,則op為單位矢量,記為i。連接o與p點并延長交平面2于(u,v,1)點,設p點坐標為(x′,y′,z′),則i=(x′,y′,z′)。
如圖所示,可知x′/u=y′/v=(z′+1)/2,且x′2+y′2+z′2=1,
可求得:i=(1+ω2/4)(u,v,1-ω2/4),其中ω2=u2+v2。

圖1 建立的坐標系
如圖2所示,I為LED發(fā)出的光線,入射到到反射面上的q點,R為反射光線,反射到照明面上的t點,T為o′到t的矢量,N為反射面上q點處的法線。

圖2 照射到t點的反射光線光路圖
根據非成像光學理論中的反射定律:


其中I=ρi,可求出矢量T的表達式T=T(ρ,u,v)。
由幾何關系可知反射面的法向向量為切平面內兩個矢量的矢乘,即N=Au×Av=(ρi)u×(ρi)v。
由i矢量可求出:

將以上三式帶入反射定律,并定義p=ρu,q=ρv,可求得T=T(u,v,ρ,p,q),則反射光線反射到照明面上的點的坐標為

其中

LED照明面垂直于z軸朝下,本文采用朗伯體LED光源,發(fā)出的光一部分直接照到照明面上,另一部分被反射后照射到照明面上,在照明面上產生一個長和寬分別為a和b的矩形均勻照明區(qū)域。根據能量守恒,光源的輻射通量應等于照明面上的輻射通量。
光源輻射到照明面上的輻射通量為Φ1=,式中第一項為直接照射到照明面上的輻射通量,Ω1為該部分入射光線對應的立體角,第二項為經過反射面反射到照明面上的輻射通量,Ω2為該部分入射光線對應的立體角,μ為反射面的反射系數(shù)。照明面上接收到的輻射通量為Φ2=∫E d S,E為照明面上的平均照度,根據能量守恒有Φ1=Φ2。由于反射面的反射率可以高達95%以上,這里為了簡單起見,不考慮反射能量損失,且郎伯型光源的發(fā)光角度為120度,則Φ1=,有Φ1=Φ2,所以E=πI0/4ab。
照到t點的光有兩束,一束為圖示的T光線,另一束為圖示的I光線的反射光線,則有E=E1+E2,其中E1為T光束在t點產生的照度,E2為B光束產生的照度。
如圖3所示,設點光源的發(fā)光強度為I,被照射面積元d S′對它所張的立體角為dΩ,則照射在d S′上的光通量為dΦ′=I dΩ=I d S′cosθ′,從而照度為E
由此可求出直接照射在照明面上的光束T在t點的照度為

圖3 光源的照度

可求出反射到照明面上的光束R的照度為

該反射過程中,入射光I的強度I與反射光R照到t點的照度E2滿足偏微分方程[20]:

其中

x和y由2.1中的結果給出。入射到反射面上的光線I的強度

E2由2.2中的結果給出。將以上結果帶入偏微分方程,列出偏微分方程的具體形式。
偏微分方程的邊界條件由朗伯型光源的發(fā)光特性和矩形均勻照明區(qū)域決定。首先設置反射面的下邊緣也為一個長寬比例為a∶b的矩形,讓直接照射到照明面上的光為一個矩形光斑。而朗伯型LED光源的發(fā)光角度為120°,因此反射面的上邊緣在一個頂角為120°的圓錐面上。用數(shù)值求解法將上式進行離散化,計算偏微分方程。設置反射面的下邊緣初始點坐標作為初始點條件,結合邊界條件,帶入計算機程序迭代計算方程,求出反射面上一系列的點的坐標。
將求解出的反射面的一系列點輸入建模軟件ProE中,擬合成一個實體反射面,如圖4所示。再將該反射面導入光學仿真軟件Tracepro中,模擬驗證最終的效果,得到照明面上的照度分布圖,如圖5所示。在中間照明區(qū)域,光照的均勻度可達85%以上,能實現(xiàn)比較理想的照明效果。因為在建模過程中,是先輸入點,形成樣條曲線,再由樣條曲線形成最終的反射面,這個過程中采用了平滑的算法,引入了誤差,使最后結果有點失真,但失真程度不大。

圖4 反射面結構圖

圖5 照度分布圖
本文通過對微分方程的數(shù)值求解求出自由反射曲面的面型,依據不同的光源和照明需求得到不同的邊界條件和具體的方程,可自由設計反射曲面。反射面型的配光結構可避免透鏡的加工成本高及對能量吸收率高等缺點。
[1]Nadarajah Narendran,Lei Deng.Richard M.Pysar et.al..Performance characteristics of high-power lightemitting diodes[C].SPIE,2004,5187:267~275
[2]Yimin Gu,Nadarajah Narendran.Jean Paul Freyssinier.White LED performance[C].SPIE,2004,5530,119~124
[3]Hueseyin Murat,Herbert De Smet,Dieter Cuypers et al..Increased lumens per etendue by combining pulsed LEDs[C].SPIE,2005,5740:1~12
[4]David L Evans.High-luminance LEDs replace incandescent lamps in new applications[C].SPIE 1997,3002:142~153
[5]D.G.Pelka,K.Patel.An overview of LED applications for general illumination[C].SPIE,2003,5186:15~26
[6]Tsunemasa Taguchi.Recent progress and future prospect of high-performance near-UV based white LEDs:from ECO lighting to medical application[C].SPIE,2009,7422:74220B
[7]Solomon W.S.Chi,Tzer-Perng Chen,Chuan-Cheng Tu et.al.Multicolor white light-em itting diodes for illumination applications[C].SPIE,2004,5187:161~170
[8]Chao-Heng Chien,Zhi-Peng Chen.An integrated LED reflector for backlight system[C].SPIE,2008,7058:705810
[9]吳仍茂,屠大維,黃志華,趙其杰.LED照明系統(tǒng)的光照均勻性設計[J].光學技術,35(1),74~76
[10]余桂英,金驥,倪曉武,鄭穎君.基于光學擴展量的LED均勻照明反射器的設計[J].光學學報,29(8):2297~2301
[11]余桂英,金驥,朱旭平,林敏.反射式發(fā)光二極管汽車前照燈的優(yōu)化設計[J].中國激光,36(s2):112~116
[12]John C.Bortz,Narkis E.Shatz,and David Pitou.Optimal design of a nonimaging projection lens for use with an LED source and a rectangular target[C].SPIE,2000,4902:130~138
[13]H.Ries,J.A.Muschaweck.Tailoring freeform lenses for illuminations[C].SPIE,2001,4442:43~50
[14]Bill Parkyn,David Pelka.Free-form illumination lenses designed by a pseudo-rectangular lawnmower algorithm[C].SPIE,2006,6338:633803
[15]Harald Ries,Ari Rabl.Edge-ray principle of nonimaging optics[J].J.Opt.Soc.Am.A.11(10):2627~2632
[16]Xutao Sun,Zhenrong Zheng,Xu Liu,Peifu GU.Etendue analysis and measurement of light source with elliptical reflector[J].Displays,27(2):56~61
[17]J.M.Gordon and Harald Ries.Tailored edge-ray concentrators as ideal second stages for Fresnel reflectors[J].Applied Optics,32(13):2242~2251
[18]Harald R.Ries,Roland W inston.Tailored edge-ray reflectors for illumination[J].J.Opt.Soc.Am.A,11(4):1260~1264
[19]Harald Ries,Julius Muschaweck.Tailored freeform optical surfaces[J].J.Opt.Soc.AM.A,19(3):590~595
[20]J.S.Schruben.Formulation of a Reflector-Design Problem for a Lighting Fixture[J].Journal of the Optical Society of America,62(12):1498~1501