车辆工程毕业设计3车灯线光源的最优设计

A题 车灯线光源的最优设计 参赛队员： 王之元 谷德峰 饶彬 指导老师 ：毛紫阳 学校： 湖南长沙国防科技大学 2 A 题 车灯线光源的最优设计 摘 要 车灯线光源的设计具有很强的实际应用意义。该问题属于单目标规划中的非线性规划问题。本文通过已知条件求出了灯光焦点，以及任一条反射光线的空间解析表达式和对应屏上的坐标位置表达式。然后建立 光子跟踪模型 进行求解。光子跟踪模型的原理是把光线粒 子化，及时跟踪光子的运行方向，最后以单位面积打到屏上的光子数来衡量光照度大小，进而反映光强度在屏上的分布规律。这是一种离散型处理方法，其本质是 计算机模拟 。这个模型中基于不同原理又提出了好几种算法：等 间距光子跟踪算法，改进的等间距跟踪算法，等效立体角跟踪算法和随机方向跟踪算法。每一种算法的原理都不一样，层层递增，一步比一步深入，并分别作图进行比较。另外还结合边界条件讨论了线光源长度的临界值以及B、 C 两点光强随线光源长度变化的规律。 在进一步讨论中，我们分析了光线直射到屏上时的情况，对二次反射的影响也做了分析，进 而证明了问题的合理性。 本文我们得出的结论是： 满足功率最小时的灯丝长度为 4.337mm，第二问的答案见下图，图的大概形状是一个心形。 关键词： 光子跟踪模型 计算机模拟 3 一、 问题重述与分析： 1、 问题重述 车灯线光源的设计是一个非常实际的问题。已知车灯的形状为一旋转抛物面，其开口半径 36 毫米，深度 21.6毫米。 经过车灯的焦点，在与对称轴相垂直的水平方向，对称地放置一定长度的均匀分布的线光源。要求在某一设计规范标准下确定使光源功率最小的线光源长度 （规范化要求略，见原题）；并对得到的线光源长度，在有标尺的坐标系中画出测试屏上反射光的亮区和讨论该设计规范的合理性。 2、 问题分析 显然在线光源单位长度光通量一定的情况时，要使光源功率最小，线光源的长度也应该较小。但线光源的长度太小了，有可能出现 C点的光强度小于额定值；线光源的长度过大，虽然能同时满足 B、 C 两点光强度的要求，但线光源的功率也增大了。我们的目的就是在 B、 C 两点光强度满足题目要求的情况下，求出最小的线光源长度。 另外还要特别注意对“光强度”这一概念的理解，我们认为它和物理学上的“发光强度”是一 致的。按光度学中的定义，发光强度是某一方向上单位立体角内所辐射的光通量大小。一般不是用肉眼可以观察到的，主要的测量仪器是前照灯检测仪（参看中国汽车检测网 前照灯检测）。其构件 一般是采用具有把吸收的光能变成电流的光电池元件，按照前照灯主光轴照射光电池产生电流的比例，来测量前照灯的发光强度 。由于本题中的光源不是点光源，直接求光强度比较困难，我们通过对光照度的测量来近似反映 B、 C 两点光强度的大小。光照度是单位面积上所接受的光通量的大小，用来衡量被照明表面明暗程度的物理量。 二、 模型假设： 1、基本假设 （ 1） 光线通过车 灯的前玻璃时能量无损失。 （ 2） 光线在抛物面进行一次反射时，能量考虑成无损失。 （ 3） 不考虑二次反射。 （ 4） 不考虑光的干涉和衍射现象。 （ 5） 截取线光源上很小的一段 dl ，可以看成是在空间呈均匀辐射的点光源。于是线光源可以看成是无穷多个点光源的叠加。 （ 6） 不考虑线光源对反射光能量的阻挡和吸收，即是说线光源不考虑厚度，反射光线可以毫无影响地穿过线光源区域。 2、 符号说明 h -车灯的深度， h=21.6mm f -车灯的焦距 4 l -线光源的长度 Q -单位长度上线光源的光通量，一个点光源的光通量可以近似记为 Qdl F - 一个点光源的光通量大小 F =Qdl I - 一个点光源的发光强度 ， I F4 三、 模型建立及求解： 定理 1： 灯泡的焦距 f =15mm 证明：我们以光轴的正方向为 x轴，竖直 方向为 z轴，水平方向 (AC方向 )为 y轴，抛物面的顶点为原点建立空间直角坐标系。如图是抛物面在 z=0时的函数图象，设函 数为 pxy 22 则焦距 f =2p ，由抛物线的定义 MNMF ，得 6.212)26.21(26.21 2 ppp 解得 p=30,所以 f =15mm。 # 因此曲面方程： xy 602 ， 6.21,0x 作出的图象见右图。 5 该模型是一个单目标非线性规划问题 Min )(lP 其中 K 是某一要求的额定光强度 , EB ,EC 分别表示线光源在 B,C 两点的光强度 我们的目的就是求出满足 P 最小的 l 光子跟踪模型： 本题中线光源可以看成是许多个点光源的叠加，只有焦点处的那个点光源通过抛物面后的反射光线是平行光束。其它点光源的反射光线都不是平行光束，因此射到屏上的总光线分布比较复杂，难以求出解析表达式。但并不是无规律所寻，给定了 空间的一个点光源 G(15,y1,0) ，（其中 2,21 lly ），再给定抛物面T: xzy 6022 （ )6.210 x 上的一任点 P(x0,y0,z0),则反射光线就可以根据 y1,x0,y0,z0 这四个参数唯一确定了。下面我们就根据光的微观粒子性结合反射光线的曲线分布建立光子跟踪模型。 这一模型的主要思想是把光能粒子化，及时纪录每一个“光子”的运动轨迹，判断最终到达屏上的具体位置。最后把每个小面积区域内的“光子 ”总数进行叠加，就得到了屏上光通量的分布图，进而求出光照度和光强度的分布。例如对于点光源 G(15,y1,0) ，（ 2,21 lly ），设其光通量为 F ，单位时间内辐射出的“光子”数为 M(为提高精度， M 可以取得大一些，如 1000， 000 个 )，则每个光子所带能量为 F /M（一秒钟内） ,且这些“光子”围绕 G点在空间是均匀发射的。一部分“光子”通过灯泡外玻璃壳直接射出（这一部分能量非常小，我们在模型的 进一步讨论中将继续讨论），绝大多数的“光子”通过抛物面反射一次后射出，还有可能出现“光子”在抛物面反射两次，在反射第二次时，我们设想该“光子”就被牢牢的地粘在了抛物面上（这种情况概率也很小，模型改进中将继续讨论）。通过把线光源选取一定步长 dl做为一个点光源，可以得到 N个点光源 (N= dll),kEkEBC2 6 我们把这 N 个点光源在屏上某一小区域内的光子数进行求和 ,就可以求得该小区域内的光照度值 我们先分析一下任一个 点光源 G(15,y1,0) ，（ 2,21 lly ）在抛物面上任一点 P(x0,y0,z0)的反射曲线 . 为了直观的理解我们在平面上做出了大致反射光线图 ,实际上这些点图并不是在一个平面上的（见图 2）。 G是点光源， P是反射点， E是 G点关于过 P点切平面的镜面对称点， PN是反射方向。 由空间解析几何的知识，我们不难得出 E点坐标： x2 15 2 3600 x0 60 y0 2 120 y0 y1 60 z0 23600 4 y0 2 4 z0 2 y2 y1 2 1800 y0 120 x0 y0 4 y0 3 3600 y1 4 y0 z0 2 4 y1 z0 23600 4 y0 2 4 z0 2 z2 2 1800 z0 120 x0 z0 4 y0 2 z0 4 y0 y1 z0 4 z0 33600 4 y0 2 4 z0 2 进而由反射光线 PN 和屏面方程 x=25015 可以求出反射光线打在屏上的坐标K(x3,y3,z3),由 7 Mathematica 计算结果见下表。公式较长，但我 们在用计算机进行模拟时需要用到 y325015y0 y12 1800 y0 120 x0 y0 4 y033600 y1 4y0 z024 y1 z023600 4 y024 z0215 x02 3600 x0 60 y02120 y0 y1 60 z023600 4 y024 z02y0 15 x02 3600 x0 60 y02120 y0 y1 60 z023600 4 y024 z02x0y0 y12 1800 y0 120 x0 y0 4 y033600 y1 4y0 z024 y1 z023600 4 y024 z0215 x02 3600 x0 60 y02120 y0 y1 60 z023600 4 y024 z02,z325015 z02 1800 z0 120 x0 z0 4 y02z0 4 y0 y1 z0 4 z033600 4 y024 z0215 x02 3600 x0 60 y02120 y0 y1 60 z023600 4 y024 z02z0 15 x02 3600 x0 60 y02120 y0 y1 60 z023600 4 y024 z02x0z02 1800 z0 120 x0 z0 4 y02z0 4 y0 y1z0 4 z033600 4 y024 z0215 x02 3600 x0 60 y02120 y0 y1 60 z023600 4 y024 z02,x3 25015 8 数据采样： 下面我们就给出具体做法，首先要进行数据采样，即是说给定一个点光源 G，让它遍历抛物面上的所有点。我们遍历的参数是和 r ( 36,0r ,0,360degree)。 其中 r 是抛物面上的一点到 x 轴的垂直距离，垂心为 xr是以 xr为圆心， r 为半径上的圆转过的角度。则抛物面方程化为sincos602rzryrx（如图 3） 当我们的 ， r 定了时，抛物面上一点也定了。 ， r 分别变动很小的 d,dr 时，光线将在抛物面上扫出一个小面积 dS1， 这个面积非常小，我们可以近似的认为是一个小平面，光线照射就可以看成是镜面反射。我们设想有一个虚拟光源 G, G 是 G 关于 dS1 的镜面对称点 .则从 G 点反射出的光线完全可以等效看成是从 G点直射出的光。设 G点通过 dS1 射到屏上的面积为 dS2,由光度学的知识可以得到 dS2 上的光照度 E dFdS2 Icosd2 其中 dF 是点光源 G 在 dS1 面上的光通量 ,d 是 G到 dS2 的距离， 是从 G射出 9 光线与 x 轴的夹角。给定点光源 G，给定 ， r 和 d,dr 时，光线在屏上的位置和照度就唯一确定了。 下面我们就建立算法来实现遍历过程。 1、 等间距光子跟踪算法 step 1 选定点光源步长 dl,确定点光源个数 N= dll,每个点光源 1秒内通过的光子数为设为 M=100000, step 2 将吸收屏均匀分割成面积为 ds 许多小正方形 ,记下位置和编号 step 3 选定 d,dr ,分别确定 和 r 的遍历次数 , N1 36dr , N2 2d ,每次遍历对应的 dS1 通过的光子数记为相同 step 4 for(I=1;I0 的图象，）。 由图可以看出斜率是单调下降的， x0=0 处斜率最小，但不为 0 ，原因是 y1不等于 0， y1 越大，图象在 y 轴的截距也将大些。由图象的单调性我们可以推知，线光源 长度下界 1l 应满足的临界条件是点光源 G(15,y1,0)在（ 21.6,-36,0）处的反射光线恰好过 C 点，由于只 考虑了（ 21.6,-36,0）的情况，还应该考虑（ 21.6,36,0）点的情况，二者求出的较小一个 2 y1 值就是长度下界 1l 。对于长度上界满足的邻界条件刚好相反。临界条件是点光源 G(15,y1,0)在（ 21.6,36,0）处的反射光线恰好过 C 点 。 经 过 解 方 程 计 算 得 到 1.53881 = l /2 =2P 时 ,C 点光照度必满足 =P; 当 B 点光照度大于等于 C 点光照度 2 倍时 , 制约额定功率 P 大小的必要条件是 C 点光照度 ,因为 ,当 C 点光照度满足 =P 时 ,B 点光照度必满足 =2P; 我们是在功率一定的条件下对线光源长度进行模拟搜索 ,所得数据点经过拟合处理后 ,得到上面的两条曲线 ,当 L 增大时 ;起初 , B 点光照度大于等于 C 点光照度 2 倍 , B 是制约额定功率 P 大小的主要因素 ,经过一段距离后 ,B 点光照度小于等于 C 点光照度 2 倍 , C 成为制约额定功率 P 大小的主要因素 ;综 合分析得L/2=2.1685 倍附近时 ,得到满足在功率一定的条件下的最大额定功率条件 ,而原题恰好是一个逆向思维 ,当额定功率一定时 ,分析求解满足最小功率条件的线光源的长度 ,可见 ,这于固定功率 ,求可满足的最大额定功率是等价的。我们作出了在线光源长度 L=4.337mm 情况下屏上的反射光亮区（见图 11），这即是第二问的答案 四、 模型的进一步讨论： 前面的模型都是不考虑直射的情况得出来的，实际上射到屏上的一部分光还有来自于灯泡直射的部分。因此严格的说光斑的形状 还要大些，但也不是圆形，而应该接近于椭圆型。只是由于直射部分的能量占的比重太小，看起来还是心型起绝对作用。但这也反映了设计规范性的合理问题。 下面我们就来推导一下直射时光照度应满足的规律 设点光源 G(15,y1,0) ，所带光通量为 Qdy1,直射到屏上 AC 方向的坐标为M(25015,y3,0) 则该点光源在 M点的光照度为 E I cos d 2 I 25000y3 y1 2 25000 2 32Qdy 14 25000y3 y1 2 25000 2 32 那么线光源在 M 点的光照度就是这些点光源在该点的积分 (y1 2/,2/ ll ) 19 E Ml2l2 Qdy14 25000y3 y1 2 25000 2 32y 1将这一部分光强也考虑进去，我们得到更加合理的线光源长度值 L=4.337mm 这个值比仅考虑反射的情况要稍微大一些。但也可以看出直射光线的能量相对于 反射光线的能量是很小的一部分，近似计算中完全可以不考虑。 另外还需要考虑 2 次反射的情况，可以根据反射光线与抛物面方程的交点来判别是否出现 2 次反射。如果有两个交点就会出现二次反射。二次反射在实际生活中一般是有能量损失的。为了设计更加合理，我们可以对发生二次反射的光线进行加权处理，比如我们可以设所有发生二次反射的光线能量损失 10%。对本题我们进行模拟求解，得到的 L=4.303mm，这个长度和不考虑二次反射的值完全一样。我们可 以认为几乎没有光线发生二次反射或者说是只有几