车灯线光源的优化设计02A课件

1、数学建模竞赛题选讲(2)第1页，共36页。车灯线光源的优化设计 （2002 年全国大学生数学建模竞赛 A 题） 汽车头部的车灯形状为一旋转抛物面，其对称轴水平地指向正前方，并已知其开口半径为36毫米，深度为21.6毫米. 经过车灯的焦点F，在与对称轴相垂直的水平方向上，对称地放置一定长度的均匀分布的线光源. 要求在某一设计规范标准下确定线光源的长度。该设计规范可简单描述如下： 在焦点F正前方25米处的A点放置一测试屏，屏与FA垂直，用以测试车灯的反射光. 在屏上过A 点引出第2页，共36页。一条与地面相平行的直线，在该直线A点的同侧取B点和C点，使AC=2AB=2.6米. 要求C点的光强度不小

2、于某一额定值（可取为1个单位），B点的光强度不小于该额定值的两倍（只须考虑一次反射）. 我们需要解决的是： 满足该设计规范的条件下，计算线光源长度，使线光源的功率最小 . 得到的线光源长度，在有标尺的坐标系中画出测试屏上反射光的亮区 . 讨论该设计规范的合理性 .第3页，共36页。1、线光源宽度与其长度相比小得多的发光体.2、光通量光源在单位时间内辐射出的光能，用 c 国际规定的眼（对光适应的眼）的灵 c 敏度来估定.3、发光效率光源每消耗1瓦特功率可辐射出的光c 通量，又叫光效.4、照 度单位面积所接受的光通量，单位为勒 c 克斯.一、几个概念第4页，共36页。 首先解决本问题的一个难点：由

3、于问题讨论的主要对象是线光源，而线光源处理起来不太容易，因此想到用微元法将线光源划分为若干小段（微元），每一段可视为一个点光源. 这样就将线光源的问题转化为了点光源的叠加问题.二、问题分析 由于该线光源照射到测试屏上，有两种途径：直射和反射. 因此我们在考虑某点的光强时，必须同时考虑两种情况 . 还需解决的几个关键问题是如何度量光强，光强的额定值为多少，如何建立起功率与光强的关系 .第5页，共36页。三、模型假设1、光在反射时及整个传播过程中，能量均无损耗;2、只计此线光源的长度，其体积忽略不计;3、此线光源均匀分布，各点同性; 且任意一点向整个空间各个方向均匀辐射能量; 4、光强额定值为1勒

4、克斯;5、不考虑光的干涉、衍射等现象;6、光强在某点叠加时，可直接取其代数和;7、只考虑一次反射 .第6页，共36页。四、符号说明 1、 C 焦距 （m ） 2、 p 线光源的功率 （W ） 3、 光通量 （lm 流明 ） 4、 发光效率 （lm/W ） 5、 E 照度 （lx(lm/m2) 勒克斯 （ lx=1 lm/m2） ） 6、E0 照度额定值 （ lx(lm/m2) ） 7、 k 发光功率线密度 （W/m ） 8、 L 线光源长度 （m ）第7页，共36页。五、模型的建立首先，建立如图坐标系如下：第8页，共36页。于是, 推出此旋转抛物面的方程为 设 (x0, y0, z0) 为此抛

5、物面上任意一点，求得过此点的切面方程为 结论: 在旋转抛物面上，每个反射点的法线都与对称 c 轴相交。 从上图，可以看出此旋转抛物面是由曲线绕轴旋转得到的。第9页，共36页。进而得法线方程为 联立对称轴方程得到交点坐标为（0，y0+c，0） 将开口半径及深度等数据代入（1）式，得出焦距 c/2=15mm=0.015（m） 然后代入，推得该线光源的范围为 0.03, 0.03m。 然后代入, y=0.015, 推得该线光源的范围为 c 0.03, 0.03m .第10页，共36页。六、模型的求解 光源功率与照度的关系 可见，在点光源垂直照射的情况下，被照面上的照度跟光源的发光强度成正比，跟表面法

6、线与光线方向之间夹角的余弦成正比，与光源到被照面问的距离的平方成反比 .第一问：计算线光源长度说明：将对B点进行分析，对于C点，所用方法及运算c 过程完全类似 . 如果光源的发光强度为 I，光源与被照面的距离是R,光线方向与被照面法线的夹角为 . 那么，这个面上的照度就是第11页，共36页。由物理知识知，光通量与发光强度之间的关系式为：联立得光源功率为 即可推出 对于功率线密度为k瓦米，发光效率为流明瓦的线光源，在距它 r 米的物体表面，照度为第12页，共36页。 用微元法将线光源转化为点光源 对于线光源来说，利用微元法的思想，将其分为若干小段（l），每一段可视为一个点光源，各自独立向四周辐射

7、光能 . 而在测试屏上的任意一受光点的照度等于各点光源在此处的照度代数和 . 这样，就把对线光源的讨论转化成了对点光源的讨论 .第13页，共36页。 求直射光的强度 线光源上任一点光源 D 直射到测试屏上 B 点处有且只有一条光线. 因为 D 发出的光是以球面波的形式向空间传播的，设该点到 B 点的半径为 r，则在 B 点处的波面面积为 据 式，可得任意一点光源直射到 B 点时该点的照度为：第14页，共36页。=25.034(米)于是得出此线光源的直射光在 B 点辐射出的总照度为：而此时线光源所消耗的功率 p=kL 可以看出此量与距离有关，也就是说，不同的点光源在 B 点所对应的照度值不同.

8、但注意到所有的差别太小，因此在实际计算中，我们对作近似处理，令它为一个常数：第15页，共36页。D 线光源上任一点 N 对应于D点的反射点 G 法线与中轴线的交点 求反射光的强度 为了便于分析，我们将图一抽象出来，并旋转一定角度后得右图：第16页，共36页。 对于线光源上的任意一点D(设坐标为(t，0.015，0),假设从它发出的光线经旋转抛物面上的点N反射后，可以射到B点上. 下面利用向量来推导旋转抛物面上满足此条件的所有反射点N的坐标 . 设反射点 N 的坐标为 (x，y，z)，由以上的证明可知，在旋转抛物面上，每个反射点的法线都与对称轴相交，且交点 G 的坐标为 (0，y+c，0). 又

9、已知点B的坐标为(-1.3，25.015，0)，点D 的坐标为(t，0.015，0). 经分析可知，点N 必须满足以下三个条件：第17页，共36页。1）NG是BND的角平分线，即： BNG= DNG故有 cos BNG =cos DNG第18页，共36页。3）点N在抛物面上，故应满足抛物面方程： 联立方程、，即可求解出与线光源上一点D对应的所有的反射点N .2）点G、B、N、D四点共面，故有从而有第19页，共36页。 但实际上我们并不关心 N 点的具体坐标值，而看重的是线光源上每一点能发出多少光线经抛物面反射后照到 B 或 C 点. 由于光源发出的光能否反射受灯反射面几何尺寸大小的限制，并且B

10、、C点的位置相对灯而言是处在斜前方，所以，N点的位置并不具有对称性 . 对于线光源上不同的点，与其对应的N点数目也不同. 但在线光源上位置十分相邻的两点，可认为它们对应的N点数目相同. 即它们发出的光中有相同数目的光线经抛物面反射后恰好照射到B点(或C点).第20页，共36页。 我们采用在小区域采集离散点的方法，使用Mathematic软件，通过改变微小l,得到线光源上当前D点对应的反射点（即N点）的数目. 以此作为此微小区域 l, l+ l 内各点所对应的N点数目 . 0.000 222-0.000 2000.000 422-0.000 4000.000 622-0.000 6000.000

11、 862-0.000 8000.001 052-0.001 0000.001 252-0.001 2000.001 452-0.001 4000.001 654-0.001 6000.001 856-0.001 8000.002 046-0.002 000t 的取值到B 的反射光的条数到C 的反射光的条数t 的取值到B 的反射光的条数到C 的反射光的条数第21页，共36页。t 的取值到B 的反射光的条数到C 的反射光的条数 t 的取值到B 的反射光的条数到C 的反射光的条数0.002 246-0.002 2000.002 445-0.002 4000.002 645-0.002 6000.00

12、2 846-0.002 8000.003 046-0.003 0000.003 245-0.003 2000.003 445-0.003 4000.003 646-0.003 6000.003 845-0.003 8000.004 045-0.004 0000.004 245-0.004 2000.004 445-0.004 4000.004 645-0.004 6000.004 844-0.004 8000.005 044-0.005 000表中t为点光源x的坐标第22页，共36页。于是，可以得到 如果把整个线光源划成个相等的微小区间，则由上述分析，第个微小区间的点发出的光线经抛物面反射，如

13、果能恰好射到B点（或者C点），那么在B点的照度为第23页，共36页。 求光源长度的优化解至此，已分别求得了E直、E反，于是可以得到通过条件，建立线光源长度的优化模型：第24页，共36页。 代入照度额定值(1勒克斯)、高压毛细汞灯和卤素灯的参数，得到以下结果：材料发光效率（流明/瓦）功率线密度（瓦/毫米）B点照度（勒克斯）C点照度（勒克斯）线光源的长度（毫米）高压毛细汞灯50302.012 11.367 742卤素灯30202.126 71.512 279第一问的问题全部解决.第25页，共36页。 线光源反射光的亮区可以近似的看作线光源上一系列离散的点光源的反射光亮区的叠加. 通过如下方式产生线

14、光源反射光的亮区： 首先我们称光能照射到的区域为亮区. 以常用的高压毛细汞灯为例（即深度= 4.2mm），来绘制出测试屏上反射光的亮区.亮区图的绘制步骤：第二问：反射光亮区的绘制第26页，共36页。1) 在线光源上等间距的选取一系列的点作为点光源;2) 在旋转抛物面上产生均匀排列的网格作为反光点；3) 对于每一点光源，计算它到所有反光点的入射光 的方向向量及其法线方向向量；4) 通过入射光方向向量及其法线方向向量计算反射 光方向向量；5) 通过反射光方向向量及反射点的坐标求出反射光 的解析式，进而求出反射光在测试屏上的投影坐标， 并在此处描点. 其中，由入射光方向向量及其法线方向向量计算反射光

15、方向向量具体过程如下：第27页，共36页。 将入射光方向向量 X及其法线方向向量 Y 均转 c化为单位向量（如下图所示）；任一点光源的向量反射光路图 所以 Z 即为反射光的方向向量. 由几何关系可知： 取入射光方向向量X 的反方向向量X；第28页，共36页。亮区图 测试屏上反射光的亮区图第29页，共36页。 第三问 设计规范合理性的讨论 要解释此设计规范的合理性，就需要解释清楚只选择A、B、C点(即25米、1.3米、2.6米这些距离上的点) 并要求这些受光点处的光强不小于某一额定值的原因这需要从几个方面来阐述 . 从实际着眼，司机在夜晚行车，难免有一定的安全隐患，如前方突然出现一横穿公路的行人

16、，或者前方有影响正常、安全行驶的障碍物等等. 要避免交通事故的发生，就需要在一定远的距离外，司机能看到这些突发事件并及时作出反应，采取措施 ( 如立即刹车等) .第30页，共36页。 而司机在遇到突发事件后总有一定的反应时间，因此刹车总距离(也称停车视距)为司机反应时间距离，刹车系统协调时间距离和刹车距离三者合一. 为了确保安全，就必须保证进入停车视距的突发事件被充分照亮(即该处所接受的光强不小于某一额定值).下面分两种情况来计算停车视距 . 第31页，共36页。 情况一 正前方发生突发事件，司机或者立即刹车，或 者调整方向，绕道而行.司机反应时间距离为 对这两种救急措施，我们认为其停车视距的

17、计算可归为一类，即计算上面定义的刹车总距离.资料显示人的反应时间一般为1s，普通公路上车速一般为60km/h，刹车系统协调时间距离和刹车距离一般为8m.下面来计算停车视距.所以其停车视距为 16.7+8=24.7（m） 但由于刹车系统协调时间距离和刹车距离随路面附着系数不同而不同，所以此种情况下的停车视距一般在2430米之间波动. 因此可以看出保证25米处光强足够大能在一定程度上确保安全性.第32页，共36页。情况二 相邻车道发生突发事件，司机立即刹车. 同时从经济的角度讲，不选取更长的距离作为规范，是因为在已经保证安全性的前提下，避免造成浪费. 以上的讨论从实际、安全、经济多角度出发，说明了

18、该设计规范是合理的 . 基于情况一，我们进一步深入. 通过查资料，我们了解到一般公路的车道宽度为 2.54 米，我们以 3 米为例，这时若有行人处于相邻车道的中线，且与车垂直距离为 25 米，正欲横穿本车道，我们来计算一下：假定行人行走速度为 1.5m/s，其余参数同上. 则易算出行人在快要到达本车道中线时，车也刚刚 驶到，自然会出现事故. 因此必须照亮相邻车道的部分情况来 使司机提前采取应急措施. 针对不同的车道宽度，可算出此值为2.54 米。这说明仅照亮正前方 25 米还不够，必须同时保证左右偏大约 2.6 米的地方也被照亮 .第33页，共36页。 7.1 优缺点评价 1、运用的数学工具简单，模型清楚易懂，可读性好，实用性强 . 2、运用近似计算的思想，将运算过程大大简化，且对结果影响不大，增强了该模型的可运行性和可移植性 . 3、引进向量，使模型统一性好，巧妙地避开了运用初等几何知识分析时需进行的大量复杂的讨论 . 4、引入发光效率，照度等物理概念来准确描述功率与光强的关系，受光强度等量，而且通过查资料得到一般测试光强的额定值等重要参数，实际背景较强 . 5、但在距离求和的运算中，考虑到毫米与米在数量级上差4个单位，因此我们忽略了毫米的微小影