车灯线光源的优化设计38083

1、车灯线光源的优化设计杜琳琳 贺思三 伍微摘要 通过对照度和光强关系的分析，阐释了求解出测试屏上任意点照度的理论可行性，并将模型转化为：max E(L) s.t 其中，为某一定值； 指线源长度为L时，C点的照度；指线源长度为L时，B点的照度；K为某一定值，其物理意义详见报告正文中（2）式。根据照度平方反比定律，可以得出每一点照度的微元表达式。理论上，可以利用积分得出屏上各点的照度，具体数值求解时，用特定剖分下的有限和来近似，再根据叠加原理算出线源整体在屏上各点的照度。对线源长度L进行搜索，直至满足目标函数。随着抛物面网格细化和线源分段间隔加密的精细化，最优线源长度L0的数值结果应趋近于其理论值。

2、搜索出L0后，将测试屏划分为5cm×5cm的小方格（小方格内部各点照度相同），依据叠加原求出各小方格的照度，并由此画出发射光的亮区。结果如下： 1最优线源长度L0为4.1mm。由最优光源长度的搜索算法得到光源长度与抛物面网格细化和线源分段间隔加密的精度有关。2反射光亮区图形参见报告正文中图9。为了切合实际，我们还作出了双车灯时的反射光亮区。最后用3DS MAX仿真了单车灯时的反射光亮区。对设计规范合理性的讨论中，兼顾了测试屏位置、测试点位置和约束条件的合理性，并针对近光、远光的不同特点，从“良好的照明”和“不眩目”两个方面给出了若干合理的约束规则。问题重述安装在汽车头部的车灯的形状为

3、一旋转抛物面，车灯的对称轴水平地指向正前方, 其开口半径36毫米，深度21.6毫米。经过车灯的焦点，在与对称轴相垂直的水平方向，对称地放置一定长度的均匀分布的线光源。要求在某一设计规范标准下确定线光源的长度。 该设计规范在简化后可描述如下。在焦点F正前方25米处的A点放置一测试屏，屏与FA垂直，用以测试车灯的反射光。在屏上过A点引出一条与地面相平行的直线，在该直线A点的同侧取B点和C点，使AC=2AB=2.6米。要求C点的光强度不小于某一额定值（可取为1个单位），B点的光强度不小于该额定值的两倍（只须考虑一次反射）。 请解决下列问题：（1）在满足该设计规范的条件下，计算线光源长度，使线光源的功

4、率最小。（2）对得到的线光源长度，在有标尺的坐标系中画出测试屏上反射光的亮区。（3）讨论该设计规范的合理性。模型假设1不考虑线光源的横截面积；2线光源上的功率均匀分布。3均分线光源足够细密时，每一小段线源可看作空间各方向均匀辐射的点源；4只考虑光的直接照射和一次反射，不考虑二次及其以上的反射；5车灯反射面绝对光滑，不存在漫反射；6测试屏无限大，光源发出的所有光线直接照射或经过反射最后落在测试屏上；7光沿直线传播，不考虑衍射、干涉等现象；术语、符号说明基本术语：光通量发光体每秒钟所发出的光量之总和，记为；发光强度发光体在特定方向单位立体角内所发射的光通量，简称为光强，记为I；照度发光体照射在被照

5、物体单位面积上的光通量，记为E；配光 照明灯具在空间各方向上的发光强度；前照灯线光源与抛物面镜组成的整体。符号： P0 线光源的总功率；L 线光源的长度，0<=L<=60mm；E0 照度的额定值，对应于光强的额定值。问题分析照度与光强的关系r·UD·A厘米量级的前照灯相对于25米外的屏而言，可看作一个点源D。以D为原点建立球坐标系（如下图），设前照灯的配光特性为I(,)。根据照度的平方反比定律1，可得屏上任一点U处的照度： （1）可见， 图1 照度与光强的关系示意图即屏上任意定点处的照度与光强成正比。又由（1）式， （2） 其中，可见，当时，任意定点处的照度与光

6、强近似成正比.模型建立目标函数 由题意，易写出目标函数如下：min P(L) s.t （3） where，为C点光强额定值。（1）、（2）式代入（3）式，可得（3）式的等价表达式如下：min P(L) s.t （4） where，C点照度额定值对（4）式作等价变形，得：max E(L) s.t （5）屏上任意点照度的计算由假设2、3，考虑线源上的微元dL均匀辐射，则：根据式（1），微元dL直接照射到屏上的照度为：dE= （6）微元dL经一次反射后照射到屏上的照度为：dE= （7）(式（7）详细推导参见附录)模型求解1最优光源长度的搜索算法 算法思路：将抛物面网格化，线光源均匀分段，各段光源视为

7、空间各方向均匀辐射的点源，根据式（6）、（7）算出点源在屏上各点的照度，然后根据叠加原理算出线源整体在屏上各点的照度。对线源长度L进行搜索，直至满足目标函数（5）式。抛物面网格化：用平行于口径面的等距离平面将抛物面分为m份，用过中轴的等夹角平面将抛物面分为n份，每份对应的角度为360/n。则整个抛物面被分为m×n个网格小面元。效果如下图所示：图2 抛物面网格化示意图符号标记： M对线源长度搜索的总次数；i 第i次对线源长度进行搜索，i=1,2,M； Li 第i次搜索时给定的线源长度； L第i+1次搜索与第i次搜索相比，线源长度的增量； Ni 第i次搜索时均分线源的段数； j 线源Li

8、的第j段，j=1,2,Ni； di 第i次搜索时的分段间隔，di×Ni = Li ；编程时，搜索线源长度时，使其线性递增，故L与i无关,取为某一定值； 分段间隔di可取某一与i无关的定值，即di=d。jNi1······FdLi 图3 线光源均匀分段示意图（F为抛物面焦点）光源长度搜索算法：step1：数据初始化。给定光源总功率P=1， i=1；线源总长初值L1=0.2mm，线源步进增量L=0.2mm，分段间隔d=0.1mm，搜索总次数M=Li|max÷(L/2)=60/0.1=600,并将抛物面分割为m×

9、;n的网格，记下网格顶点数据；step2：将线源分成Ni= Li÷d段，把每一段都看作点源，并编号j=1,2,Ni；step3：j=1；step4：根据反射定律，计算出点源j所发出的光线经每个网格顶点反射后落在屏上的对应点（注：反射镜面取该点理论切平面，并非网格平面，如下图4；具体方法参见附录），求得每个网格面元在屏上的映射面元。网格面元切平面j映射面元屏图4 网格映射step5：依据照度的平方反比定律和每个映射面元内光通量均分的原则，遍历网格，求得屏上B、C两点的照度EB(j)|i，EC(j)|i；step6：if j=N then go on； else j=j+1； goto

10、step4；step7：利用叠加原理，求得屏上B点的照度EB|i =EB(j)|i。C点类似；step8：if i=M then go on；else i=i+1； goto step2；step9：记额定照度OL|i=min(EB/2K，EC)|i 。若OL|i0 = max(OL|i)，则i0对应的光源长度L0即为最优值。2亮区求解算法按照上面的方法划分抛物面，并将测试屏划分为5cm×5cm的小方格（小方格内部各点照度相同）。计算每一小格的照度，具体算法：step1：光源长度取L0, 线源细化精度为L/2=0.1mm，屏上所有小格的照度为零；step2：将线源分成Ni= L0&#

11、247;d段，把每一段都看作点源，并编号j=1,2,Ni； step3：j=1；step4：遍历网格，计算点源j发射的光线在抛物面的每个网格面元中心点反射后落在屏上的位置和该点的照度；step5：利用叠加原理，求出各个小方格内的照度E(x,y)|j；step6：if j=N then go on； else j=j+1； goto step4；step7：屏上各点照度为：E(x,y) =E(x,y)|j，即为所求。结果及分析光源长度采用先粗搜，后细化的策略：1线源均匀分段间隔取为0.1mm，搜索范围060mm，抛物面21×36、42×72、84×144、168&#

12、215;288、336×576网格化；2线源均匀分段间隔取为0.05mm，搜索范围060mm，重复1.中网格化的工作；3线源均匀分段间隔取为0.01mm，搜索范围06mm（由1、2的搜索结果缩小搜索范围），重复1.中网格化的工作；随着抛物面网格细化和线源分段间隔的加密，求得线源的最佳长度L0见下表1： （单位mm）线元个数 面元个数720295211952480961929603006.05.24.44.24.06006.105.204.604.4030006.005.244.724.424.12 注：“”表示未算。表1 不同细化精度下线源长度的最佳值由表1可见，随着抛物面网格的细化

13、，L0有减小的趋势；随着线源分段间隔的加密，L0有增大的趋势。相对而言，网格细化精度对L0的影响更大。我们取最佳长度为4.1mm。不同搜索精度下B、C两点相对照度如图5、6所示：图5 不同线源长度对应的B、C两点的相对照度（线源均匀分段间隔取为0.01mm）图6 不同线源长度对应的B、C两点的相对照度 （线源均匀分段间隔取为0.01mm）由图5、图6可知，线源长度很短时，B、C点照度几乎为0；随着线源长度的增加，B、C点照度都是先增大，后减小；且B点照度的峰值先出现，C点照度的峰值后出现。定性分析亦可得出相同的结论。抛物面镜测试屏CB1、 线光源的长度趋近于零时，相当于位于焦点处的点光源，此时

14、反射光平行射出，测试屏上的亮区面积等于反射面的口径面积。B点和C点都没有被亮区覆盖，照度都为零，如下图7(a）：F (a)2、 随着光源长度增加，亮区面积逐渐扩大，先覆盖B点，后覆盖C点。长度每增加一点，就会产生新的亮区。光源长度从零开始增加时，亮区最远照射点到中心的距离从36mm逐渐增加，如下图7(b)、7(c)：抛物面镜测试屏CFB(c)抛物面镜测试屏CFB(b) 图7 B、C点照度随线源长度变化而变化趋势的定性分析亮区示意图根据亮区求解算法画出的亮区如下图8、图9所示：图8 不同细化精度下的亮区 图9 在有标尺的坐标系中测试屏上反射光的亮区（网格细化为168×288）由上图可以

15、看出反射光亮区特性如下：在屏上反射光的亮度从中心向四周渐渐减小；亮区关于X轴、Y轴对称；亮区在X=0附近内凹；亮区呈条带状分布；X的范围为 -3.65,3.65，亮区水平方向长为7.30米；Y的范围为-0.95,0.95 ，亮区竖直方向最长为1.90米；X=0处Y的范围是0.75,0.75 ，竖直方向最长为1.50米。图上每一种灰度代表一个亮度范围。从图上可以看到灰度存在一些不连续的点，而实际中屏上的亮度变化应该是连续变化的。这主要是由于抛物面细化精度不够，使得映射面元与屏上网格不匹配造成的。双车灯 设左、右前照灯位置坐标分别为（0，0）、（-1.5，0）。 （单位：米）可得如下图所示的反射光

16、亮区：图10 双车灯时测试屏上的亮区由上图可以看出反射光亮区特性如下：亮区呈条带状分布，亮区主要集中在-4.5<x<36,-1<y<1的矩形条带内；在屏上反射光的亮度从中心向四周渐渐减小，亮区关于X轴、Y轴对称。亮区在X=0附近内凹。当线源长4.1mm，其总功率为时，任选屏上亮区中的几点，得到其直射照度与反射照度如下表2：点位置(1300,0)(2600,0)(0,300)(150,200)(-200,-600)反射照度4.48005e-0067.30968e-0073.63534e-0065.39973e-0065.19025e-006直射照度5.06375e-009

17、5.00345e-0095.27119e-0095.20749e-0095.20683e-009比例884.73146.093689.6631036.92996.816表2 直、反射照度对照表可见，反射照度是直射照度的几百倍，这是由于抛物面镜对光线有汇聚作用，测试屏上光强应主要来源于反射光，故可以不考虑直射的影响。仿真利用3DS MAX仿真线源经过抛物面镜反射后投影到25m以外的测试屏上，效果如下图：图11 3DS MAX仿真图11与计算所得的亮区图9相比较，可以发现相差不大，从另一侧面说明了亮区求解算法的可行性。设计规范的合理性1测试屏位置的合理性 根据4，“配光性能应在前照灯基准中心前25

18、m，.的铅垂配光屏幕上测定”。题给设计规范中测试屏的位置与之相符，再综合考虑司机视野范围的限制、夜间行车速度等因素，认为25m的测试距离是合理的。2测试点位置的合理性据网上资料，汽车左右两车灯间距一般为0.91.5m， 单车道宽度一般为45m，车灯在25 m处的水平照亮范围应大于路面宽度5m，以保证能发现路面任意处的来车。本文外侧测试点C对应的水平照射宽度5.2m>5m, 较合理。另外，我们查到一些汽车前照灯的配光特性，如下图12：图12 一些汽车前照灯的配光特性 可以发现，亮区多分布在水平方向上，因为从行车安全考虑，水平方向的信息更为重要。 题给规范中测试点B、C的位于与线源平行的水平

19、方向上，合理。3约束条件（测试屏上照度）的合理性先引入三个术语4：配光：灯具发射可见光的光度(照度或发光强度等)分布。近光：当车辆前方有其他道路使用者时，不致使对方眩目或有不舒适感所使用的近距离照明光束。远光：当车辆前方无其他道路使用者时，所使用的远距离照明光束。资料4显示，“前照灯的配光应使其近光具有足够的照明和不眩目，远光具有良好的照明”。可见，近光和远光应区别对待。1）前照灯照明光束为远光时，所给设计规范的合理性讨论远光灯的主要目的是发现道路前方较远处的车辆，并使对方知道有车来了。主要因素为“良好的照明”。前面已讨论了B、C点位置的合理性。题给设计规范中给定B、C点照度额定值之比为2，其合理性有待探讨，但从视觉角度出发，观察点偏离轴线的距离加倍时，射入观察者眼中的光强减半，具有一定的合理性。但是仅考虑B、C两点的照度，信息量略显不足。为了强调足够亮这一概念，约束条件中还应加上：