车灯线光源的优化设计河海大学.doc

车灯线光源设计的数值逼近 丁根宏 韩喜君 河海大学理学院，南京 210098 摘　要 本文通过在抛物面内内接小平面镜的方法，来模拟抛物面反射光线的汇聚作用。统计能反射光线到指定点的小平面镜的个数，算出该点的光强。在满足指定的B、C两点光强大小的条件下，求出使线光源总功率最小的灯丝长度为3.832毫米。 关键词 抛物面 光强 汇聚作用 数值逼近 一、问题背景介绍 安装在汽车头部的车灯的形状为一旋转抛物面，车灯的对称轴水平地指向正前方，其开口半径36毫米，深度21.6毫米。经过车灯的焦点，在与对称轴相垂直的水平方向，对称地放置一定长度的均匀分布的线光源。要求在某一设计规范下确定线光源的长度。 该设计规范简述如下：在焦点F正前方25米处的点放置一测试屏，屏与FA垂直，用以测试车灯的反射光。在屏上过A点作一条与地面相平行的直线，在该直线A点的一侧取B、C两点，使米。要求C点的光强度不小于某一额定值，B点的光强度不小于该额定值的两倍（只须考虑一次反射）。 在满足该设计规范的条件下，计算线光源长度，使线光源的功率最小。 为了求解这一问题，首先建立坐标系，通常由解析方法分析得到线光源上发光点与反射面及测试屏之间的一系列反射方程。对于灯丝上的某一点S，直接求出抛物面上可将点S发出的光线反射到B点或C点的反射点，此处的反射点并非为抽象意义上的点，而为一微元D1。S点发出的光线经反射后，在测试屏上得一微元D2，D2、D1之间有一个面积比，即为雅可比行列式的绝对值，而D1处光强可求，所以D2处的光强可求，最后对灯丝积分，得B点（C点）处的光强。但是，用该解析方法解题有困难，对于抛物面上的某些点，其雅可比行列式的分母为0，无法实现数值积分。 下面通过在抛物面内内接小平面的方法，来模拟抛物面反射光线的汇聚作用。 二、内接小平面镜逼近 1、建立坐标系 如图建立坐标系，以车灯的底部为原点，以FA为z轴，与其垂直的过原点的水平线为y轴，根据右手法则可确定x轴，各轴均以毫米为单位。于是在该坐标系下，车灯的抛物面方程可表示为. 图1 光线反射示意图 图1中，S为线光源上某一点，R为在抛物面上的反射点，M为S关于R处切平面的镜像点。 2、用内接小平面镜逼近抛物面 将x轴按长度等分，将y轴按长度等分，这样就在xoy平面形成若干网格。每个网格是一个矩形，矩形网格中的上侧和右侧的边是实边，另外两条边是虚边，以保持矩形网格划分的完备性。设一矩形的四个顶点分别为P1、P2、P3、P4。这四个平面上的点沿z轴投影到抛物面上，得到四个抛物面上的点Q1，Q2，Q3，Q4。z1、z2、z3、z4均满足方程，容易证明Q1、Q2、Q3、Q4在同一平面上，于是Q1Q2Q3Q4构成一个四边形，记为，此即可看作是在抛物面内部内接的一个小平面镜。在灯丝上取定某一点S后，以点Q1、Q2、Q3、Q4为反射点，反射后在处的测试屏上得四个点T1、T2、T3、T4，这四个点也构成一个四边形T1T2T3T4，不妨设为Г。若四边形Г包含B点(点)，则在平面镜上存在唯一的一个反射点Q，把从点S发出的光线反射到B点(点)。 固定灯丝上的S点，由于抛物面在xoy平面上的投影是一个有限区域，故可以用若干个矩形网格覆盖这一有限区域。再对矩形网格进行遍历，若反射光线经过B点(点)，则对反射该光线的平面镜计数一次，于是得到可将此S点发出的光线反射到B点(点) 的平面镜个数。现将灯丝等分成若干小段，对各小段做一次上面的遍历，就能得到从Si发出的光线反射到B点(点)的平面镜总数Ni，这里Si为长度微元，用计算机模拟计算时，Si代替灯丝上的一小段，但可以将小段长度尽可能取得很小以减小误差。 因为旋转抛物面的大小相对于其到测试屏的距离很小，所以所有经抛物面反射到B点的光线的路长可近似取为，而反射到C点的光线的路长可近似取为. 设灯丝上一点Si处的光强为I1，则经一块平面镜反射到测试屏上B点(点)处的光强为，这里L为光线经过的路长，即为灯丝上一点Si发出的光线经所有平面镜反射后在B点(点)处的光强。于是，照此方法，灯丝上各小段发出的光线经所有的平面镜反射后在B点（C点）所形成的光强可求。 由物理学知，任一光源发出的光强与该光源所消耗的功率成正比，比例系数为常数k. 设灯丝单位长度(每一小段)的功率为w，灯丝长为m个小段，B、C两点的总光强分别为IB、IC，而C点的光强度不小于某一额定值I0. 则原问题可转化为下列数学模型。 （1） 其中，这里NB、NC分别为将m个小段所发出的光线反射到B、C点的平面镜数。由于比例系数k为常数，额定值I0也为常数，不妨设常数k、I0均为1。因此，可把光强与功率等价，约束条件可简化为IC≥1、IB≥2。 下面讨论求