车灯线光源的优化设

车灯线光源的优化设计四组牛振宇戚亚萍石艳迟摘要论文就传统的旋转抛物面形的汽车前灯的线光源优化设计问题进行了讨论。文章根据物理学中的反射定律、能量守恒定律和光强与能流密度的关系，运用空间解析几何、微积分，建立了最优化模型。再结合Matlab和Fortran,用数值方法对模型进行求解。本文建模的主要思想分为两个部分。首先，求光路：线光源上一点P0(x0,y0,z0)朝某一方向(θ，φ)发出一条光线g在抛物面被反射，入射光线与反射光线服从反射定律，根据空间解析几何和二元函数的微分法则，求得法线方程。接下来求反射线方程。论文运用解析几何知识，用入射线与法线的方向数求出反射线的方向数，再通过反射点P1写出反射线方程。得到反射线方程后，再另其与测试屏平面的方程联立，得到反射线在测试屏上的光点的坐标P3，同时可以算得反射光线与测试屏的夹角的正弦值。然后由光强、能流密度和功率P的关系，把光强I与功率P联系起来，因为它们之间存在线性关系。在空间球坐标系下取一束光束（θ，θ+△θ），（φ，△φ），这一束光的功率为△P光束经抛物面反射后在测试屏上形成面积为△S的光斑，此时，由能量守恒定律，求得该区域的能流密度ω，然后积分便得到了测试屏上任意点P(x,y,z)的光强的积分表达式（未考虑直接照射光线）。结合设计规范的要求（即最优化问题的限制条件），就得到了问题一的最优化数学模型。模型的求解是论文的一大特色，这里采用直角坐标与位于离散的点光源处的空间球面坐标相结合的方法，将复杂的解析几何运算转化为离散的运算，使得求解过程教容易在计算机上实现。而且模型求解过程中，采用一种全新的方法计算光强，在处理复杂的问题时，其极强的可操作性使它可以应用于几乎所有的同类问题。离散化处理时，对光源发出的光线的方向，在球面坐标中对出射方向角进行二维离散，求出离散的点光源在测试屏上的光强分布，并根据光线密度与照度以及光强的关系，确定出所求点的照度，进而得到模型的最优解。对某一线光源长度为2L，将其离散化为长度为△L的若干个点，然后由每一点得到测试屏的一个照度分布数据Di(S)，共有N=2L/△L个照度分布数据，该2L长度的数据组成一个集合Data(2L),称为光源长度为2L时测试屏的相对照度数据。为求得规范标准下的最优解，论文引入一个灯丝功效函数：W（L）W(L)=F(C，L)当F(B，L)2F(C，L)时W(L)=F(B，L)/2当F(B，L)2F(C，L)时功效函数的物理意义表示在达到设计规范（F(C，L)=1,F(B，L)=2）下，所能节省的功率，我们的目的是为了求最小值L，这样的引入更利于了问题的求解，由于这里已将约束条件该为比值形式。反射光照图的求解是通过计算机编程在三维域进行搜索，得到测试屏上的光照轮廓。根据常识以及相关资料表明，相关的图示见论文求解过程。本文完整求解了模型所提的问题，并作了相关讨论，模型求解所得到的最优灯丝长度为4mm，这时B点和C点的光强之比为CBII:=2.49，本文还在模型的假设下，得到了灯丝长度与功效函数（相当于节能比例）的离散函数关系,还得到了测试屏水平轴线Y=0，Z=25.015上固定功率的照度随灯丝长度2L以及X坐标变化的离散函数关系。关键词：汽车前灯线光源旋转抛物面光强一、问题重述安装在汽车头部的车灯的形状为一旋转抛物面，车灯的对称轴水平地指向正前方,其开口半径36毫米，深度21.6毫米。经过车灯的焦点，在与对称轴相垂直的水平方向，对称地放置一定长度的均匀分布的线光源。要求在某一设计规范标准下确定线光源的长度。该设计规范在简化后可描述如下。在焦点F正前方25米处的A点放置一测试屏，屏与FA垂直，用以测试车灯的反射光。在屏上过A点引出一条与地面相平行的直线，在该直线A点的同侧取B点和C点，使AC=2AB=2.6米。要求C点的光强度不小于某一额定值（可取为1个单位），B点的光强度不小于该额定值的两倍（只须考虑一次反射）。请解决下列问题：（1）在满足该设计规范的条件下，计算线光源长度，使线光源的功率最小。（2）对得到的线光源长度，在有标尺的坐标系中画出测试屏上反射光的亮区。（3）讨论该设计规范的合理性。二、问题假设(1)线光源在直射与反射的过程中不计能量损失。(2)不考虑线光源的折射、散射现象。(3)将4mm线光源等效成连续均匀分布的点光源。(4)线光源的总功率为p。(5)不考虑光在大气中的折射现象。三、问题的分析:光强I与能流密度的关系为I=c,c为光速,我们采用无量纲化计算,可以取c=1。而能流密度是垂直通过单位面积的功率,因此可得:dSdPI四、符号说明：直线在xoy平面上的投影线与x轴的夹角；)20(：直线与z轴的夹角；)0(i：由线光源发出的入射线；n：旋转抛物面的内法线向量；l：一次反射线；：反射线与测试屏的夹角；0),,(zyxf：旋转抛物面方程；fI反射光光强dI直射光光强L线光源长度的一半F（X，L）灯丝长为2L时X点的照度（标量）。W（L）长度2L灯丝的功效函数W（L）：当F(B，L)2F(C，L)时W(L)=F(C，L)当F(B，L)2F(C，L)时W(L)=F(B，L)/2功效函数的物理意义表示在达到设计规范（F(C，L)=1,F(B，L)=2）下，所能节省的功率。五、模型建立及求解图一坐标系的建立首先建立空间坐标系，取旋转抛物面的顶点为坐标原点，FA连线方向为Z轴方向，X轴与线光源平行。在此坐标系中，旋转抛物面方程为：0)(350),,(22zyxzyxf(1)测试屏平面方程为：015.25z焦点F及A、B、C点坐标分别)015.0,0,0(F，)015.25,0,0(A，)015.25,0,3.1(B，)015.25,0,6.2(C。根据题意和物理学原理，可得模型：dxdyIILdf）（）（总MinPs.t.2BI1CI问题一就是上述最优化问题中使总P最小的L。求入射线与入射点在线光源上任取一点),,(0000zyxP，该点朝某一方向),(发出一条光线，该入射线方程可表示为：cossinsincossin000zzyyxx(4)其中],[0LLx，00y，015.00z。联立)1(式与)3(式可解得入射线与旋转抛物面的交点),,(1111zyxP，即为入射点。直线与曲面的交点可能会出现两个解的情况，但由于入射光线具有方向性，且从光源到入射点的方向与入射光线的方向相同，故可根据空间向量间的方向关系舍去不合理解。求反射点1P的内法线n方程及0P点关于n的对称点),,(2222zyxP内法线方程:11111zzfyyfxxyx(5)矢量20PP法线垂直:20PP0)1,,(11yxff(6)20PP连线的中点满足1P点的法线方程:10210210222211zzzfyyyfxxxyx(7)对(6),(7)联立,求解得2P的坐标:),,(2222zyxP.求反射线方程121121121zzzzyyyyxxxx(8)求反射线方程与测试屏平面的交点3P和夹角联立（2）式和（7）式，求的反射线方程fl在测试屏上的投射点),,(3333zyxP11212133)(xxxzzzzx11212133)(yyyzzzzy015.253z其中，33,yx与,,0x有关：),,(033xxx，),,(03xyfI与测试屏平面的夹角的正弦值：21321321313)()()(sinzzyyxxzz（9）求测试屏上一点的反射光光强首先将线光源离散化，对于线光源上的一小段),(00xxx，方向为（,），立体角为),(),(。该立体角内的光束经旋转抛物面反射后在测试屏上形成光斑区域ΔS，面积即为ΔS。设这束光束的功率为ΔP,则ΔS上的平均光强：sinSPI（10）从而得到ΔS上的反射光强为LLfdxxII)(同理可求得直射线在测试屏上的光强为dI。线光源的总功率为：dSIILPdf)()(总（11）其中，Σ为测试屏平面。对于固定的光源功率，设灯丝长为2L时在测试屏上X点处的照度为F(X,L).根据物理意义认为，当灯丝长度一定时，光强和照度有线性关系。因此用照度来求解不会影响问题的结果。由于这里已将约束条件该为比值形式，为了简化求解，下面还将用照度来求解。线光源上任一点T，该点发出的到达配光屏的光线有两种：直接照射和反射后到达（这里只考虑一次反射）。球面坐标下点T发出的某一条光线，其方向数为（cos,sinsin,cossin），在θ和φ有一微小变化时光线的方向也会发生变化。现在让θ在0到2л变化，φ在0到л变化，然后取足够小的微小角度dθ和dφ，对点T发出的所有光线进行离散化处理，即在取值范围内，每隔一定间隔取一对θ和φ值，然后对点T在每一确定的方向射出的光线进行如下处理：利用光学基本定律：入射角=反射角，可以根据每一条入射光线求得其反射光线的方向数，进而求得该反射光线与配光屏的交点Mi(x,y,z)，设每个交点是一个点元Di(x,y)。（因为点元由坐标x和y所唯一确定，这里略去了坐标z=25.015。对线光源上确定的一点T发出的光线，如果存在点元Di（x,y），则该点元也唯一对应点T发出的一条光线T（θi,φi），如图二。图二设CD方向数为（zyxrrr,,）,则有：222sinzyxzrrrr（12）离散化处理时，如果对θ和φ的间隔取得足够小，就可以得到配光屏上足够多的点元。根据这些点元，用如下方法求得某一确定点M的照度。（1）在屏上以M为中心确定一个以R为边长的正方形。（2）统计该正方形内点元的个数Nd。（3）点D的照度为：F（D，L）=Nd*sinα（13）由此方法进行编程可以依次求得配光屏上任意点的照度。对线光源上某一定点S，它对应会有配光屏的一个照度分布数据Di（S）。对某一线光源长度为2L，将其离散化为长度为△L的若干个点，然后由每一点得到配光屏的一个照度分布数据)(SDi，共有2L/△L个数据，该2L长度的数据组成一个集合Data(2L),称为光源长度为2L时配光屏的照度数据。对每一个照度数据，求灯丝功效：W（L）。取线光源长度初值L为各个值时，由照度数据Data(2L)得到测试屏密度分布如下：L=0.5mmL=2.0mmL=3.0mmL=4.5mm由以上四图可以知道当线光源长度大于4.5mm之后，虽然B、C点满足要求，但会造成灯泡功率的浪费，即不是最优设计，因此取线光源长度为1mm至9mm进行搜索，求解满足规范的最优设计是合理的。取线光源步长为△L=0.01mm，初始值为2L=1mm进行搜索，得到功效在各个灯丝长度下的函数值。解的结果分析当L=2.0mm时，功效函数达到最大值。即最优灯丝长度为4.0mm，这时B点和C点的光强之比为CBII:=2.49。其图形如0.00E+01.00E-32.00E-33.00E-34.00E-35.00E-3x-线线线线(线)0.00E+01.00E+72.00E+73.00E+74.00E+75.00E+7f(x)图五灯丝长度2L=4.0mm时，功效函数W（L）达到最大图六L=2.000时，反射光在测试屏上所覆盖的区域-3.00-2.00-1.000.001.002.003.00x(m)0.501.001.502.002.503.003.504.00L(mm)图七在测试屏水平轴线Y=0上固定功率的照度随灯丝长度以及坐标变化的灰度图由上图可知，当灯丝功率一定的情况下，如果灯丝长度很小，则照度比较集中，随着灯丝长度的增加，照度逐渐分散。如果灯丝过短，虽然光照强度很大，但是照明区很小，达不到汽车司机所必须的明视区，对驾车会产生危险；如果灯丝过长，则虽然照明区很大，但是光照强度很低，也会对行车安全造成危害，而且对能源也是一种浪费。本模型求得的灯丝长度4.0mm，既满足了汽车行驶中司机对路面