《A车灯线光源最优的设计3数学建模(共19页).doc》由会员分享，可在线阅读，更多相关《A车灯线光源最优的设计3数学建模(共19页).doc（19页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上一问题重述汽车头部车灯的线光源的设计对车灯的照明光强有重要影响。一种可能的设计规范标准为，在焦点F正前方25米处的A点放置一测试屏，屏与FA垂直，用以测试车灯的反射光。在屏上过A点引出一条与地面相平行的直线，在该直线A点的同侧取B点和C点，使AC=2AB=2.6米。要求C点的光强度不小于某一额定值（可取为1个单位），B点的光强度不小于该额定值的两倍（只须考虑一次反射）。假设车灯的反射面是旋转抛物面。车灯的对称轴水平地指向正前方, 其开口半径36毫米，深度21.6毫米。线光源对称地一定长度的均匀放置在经过车灯的焦点，在与对称轴相垂直的水平方向。旋转抛物面及其内部线光源的
2、空间位置如图1.1图1.1 旋转抛物面及其内部线光源的空间位置 在满足该设计规范的条件下，线光源的功率存在一个极小值，此时线光源长度最短。二基本假设1、 线光源没有宽度，发光功率均匀分布。2、 线光源上任意一点发出的光的波阵面为球形。3、 光线在被反射时没有能量损失，也就是说，光线被全反射。4、 光线在空气中传播时没有能量损失，也不会遇到障碍物。5、 不考虑二次反射的情况。6、 光线被测试屏完全吸收，即在测试屏处检测的光强就是照射到此处的光强。三符号说明r 旋转抛物面的开口半径h 旋转抛物面的开口深度f 焦距I 光强度（单位面积的光功率）P 光功率四问题分析和解答4.1模型中基本的数学推导4.
3、1.1旋转抛物面的焦距根据题意，（本文中所有长度单位为mm），。根据公式：焦点附近的线光源可看成许多点光源的线性叠加。一个基本的结论是若点光源在焦点位置，它发出的光经过抛物面的反射形成平行光。4.1.2测试屏反射点坐标公式给定线光源上一个点光源,它与焦点的距离为d,它发出的光线经旋转抛物面上的一点的反射,反射光线与测试屏的交点为,下面的工作是建立与和的关系。本文中，旋转抛物面的对称轴设为z轴正向，顶点在原点，水平面为zy面。则其方程为： （1）显然, 坐标为。设，则坐标满足方程（1） （2）根据空间解析几何，曲面在点的法线l的法向量为= （3）则过点的法线方程为 （4）设关于法线l的对称点P的
4、坐标为，则和P连线的中点在法线上，得： （5）另外：显然，故 （6）联立（5）（6）两式，解得 （7）其中过和两点的反射光线的方程为 （8）它与测试屏的交点为。在测试屏上，将（7）带入（8），得交点坐标 （9）其中由（9）式可见，如果（即问题1中B点和C点的情况），则或这说明对线光源上的任一点光源，对测试屏上过轴与屏的交点并与线光源平行的直线上的点的光强有贡献的抛物曲面上的反射点分两部分：一是线光源与轴线所确定的平面与曲面的交点，即；一是由空间点，满足 。4.2问题一4.2.1逆向模型1) 点光源模型点光源可以看作一个没有大小的发光点。发出的光线均匀的向各个方向发散。2) 点光源的能量模型a)
5、 经过以点光源为球心的任意一个球面的光的总能量相等；b) 由a)和球面积公式可得，点光源照射到任意一点的光强度，与该点到点光源的距离平方成反比。c) 本问题中，可以假设单位长度的线光源功率为P，则点光源的功率大小也为常数P；3) 光强度函数光的强度在物理学中有具体的定义。但是在本问题中，并没有要求、也无法求出光强度的具体数值。又由于点光源的功率已假定为一常数P，于是某点的光强度仅与该点到点光源的距离d的平方成反比，因此，在本文中，可以定义光强度函数I为： (10)4) 线光源模型线光源可以看作是无数点光源紧密排列在一起组成的。以下的讨论中，可以先讨论点光源对B、C点光强度的影响，然后对光强度函
6、数在线光源所在的直线上进行积分，就可以得到线光源在B、C点产生的光强度函数I (11)5) 旋转抛物面模型旋转抛物面可以看作一条抛物线沿着它的对称轴旋转180度所得的曲面。其中任意一条抛物线都可以起到反射镜的作用。6) 坐标系的建立与各物体的位置关系坐标系与旋转抛物面如图4.2.1所示。图4.2.2.1 坐标系与旋转抛物面位置7) B、C点光强度的粗略模型显然，照射到B、C点的光线来自两部分，一部分是线光源的直接照射，另一部分经过旋转抛物面的反射后照到B、C点。于是，B、C点光强度应该是两部分分别在两点引起的光强度的叠加。8) 直射光线在B点引起的光亮度大小很容易证明，线光源上每一个点光源都可
7、以照射到B点。设线光源两端点分别为：和。某点光源坐标，到B点的距离d(x)，则线光源在B点引起的光亮度大小为： (12)同理可以得到该段线光源在C点的光亮度大小。9) 反射光线在B的点引起的光亮度大小a)一些基本关系式线光源上任取一点，旋转抛物面上任取一点，两点确定一条入射光线：方向向量B点坐标，与旋转抛物面上一点确定一条反射光线：方向向量：过抛物面上一点的法线方程：方向向量：由于入射光线、反射光线和法线满足如下两个条件：三线共面；法线平分入射光线与反射光线的夹角。由以上两点可以得到如下关系式：因为点和点分别为线光源和旋转抛物面上的点，所以还要满足如下约束条件：b) 反射的讨论由入射光线、反射
8、光线和法线满足的两个条件可知，a、b向量归一化后的和向量与法线同向。于是得到如下关系式：其中：分别为向量a、b的模；与关系式、联立并讨论得：1. 当时，立体问题转化为平面问题。即由点光源发出的光线只经过水平抛物线镜面反射到达B点。2. 当，用x1，z1表示出y1代入，并舍掉数值较小的项，对得到的方程式讨论得到：对于与之间的点光源，总有两个解。也就是说，这两个点之间的点光源经过除水平抛物线镜面之外的抛物面反射（我们可以称之为立体反射）可以照射到B点，而且有两条光线可以到达B点。我们称起反射作用的点为反射点。同样的原理可以得到，对于C点，这样的临界点为：(-1.56,0,15)和。3. 水平抛物线
9、镜面反射光线的讨论由光的可逆性原理，当线光源上某一点光源发出的光经过水平抛物线镜面的反射，可以到达B点时，由B点出发射向旋转抛物面的光线经过反射也会到达该点光源。如图4.2.2所示建立平面坐标系： 图4.2.2.2则水平抛物线镜面的方程： 旋转抛物面的开口直径为72毫米，而测试屏与旋转抛物面的焦点的距离为25米，在B点看旋转抛物面，张角0.072/250.0029；所以，可以作如下近似：由B点出发射向旋转抛物面的光线可以看作一束平行光线。如图4.2.3所示为一束平行光束经平面反射后的光线：图4.2.2.3设入射光线的斜率为k，则k的近似取值为其中，P点为水平抛物线的焦点，(XB,YB)，(XP
10、,YP)分别为B，P点的坐标。入射光线经平面抛物线镜面反射后，将与线光源所在直线产生一系列交点，这些交点就是能够照射到B点的点光源。计算每个点光源在B点产生的光亮度大小，再对这一系列点光源进行积分，就得到了线光源经平面抛物线镜面反射后在B点产生的光亮度大小I2。4. 立体反射光线的讨论我们已经得到，对于B、C点，经立体反射可以到达的临界点分别为：与、与。由于B或C距离到旋转抛物面上某点的距离与旋转抛物面的线度相比，旋转抛物面的线度很小，可以忽略。因此，我们用B点到焦点的距离来近似表示B点到抛物面上反射点的距离。设为dB，用d0表示某点光源到反射点的距离，于是可以近似得到这样的结论：于是，某点光
11、源到B点光程为：由此，立体反射在B点引起的光亮度函数也就随之确定了。在前述区间和线光源长度区间的公共区间上积分，就可得到该段线光源经立体反射在B点产生的光亮度I3。10) B点的光亮度由直射，水平抛物线镜面反射和立体反射三部分叠加而成。即B点的光亮度：C点同理。11) 搜索模型的确立由此,对于一个给定的灯丝长度，我们已经确定出了B点和C点的光的亮度.在此基础上,用计算机程序实现从一个较小的灯丝长度开始遍历,知道找到一个灯丝长度满足B点光的亮度是C点光的亮度的两倍.12) 搜索结果利用程序搜索的结果为满足条件的最小灯丝长度为3.15毫米.4.2.2 面积模型4.2.2.1模型分析通过测试屏坐标点
12、公式的计算,我们可得到以下结论:(1) 对于线光源上的某一点光源发出的光线经旋转抛物面上任意一个交线圆上的或部分的半圆反射,在测试屏上的曲线近似为一个圆。图4.2.3.1是取的交线圆的部分曲线作为反射点,取作为点光源得到的测试屏上的轨迹。图4.2.3.1（2）线光源上的某点光源在抛物曲面上反射后在测试屏上形成的点集， 的交线圆的部分曲线作为反射点得到的测试屏上的轨迹圆（其半径最大）及其内部的平均光强是其余部分的一半（如图4.2.3.2） 图4.2.3.2设图示的圆内面积为S1，其光强为，其余面积为S2，光强为，点光源的反射光功率为P，则下式成立：（3）线光源上的两个不同的点光源（除焦点外）在测
13、试屏上形成的光斑图形相似，其面积比等于最大轨迹圆半径的比的平方。（4）根据测试屏坐标点公式，可计算出只有当线光源长度大于3.12时才有反射光线经过C点；大于1.56时才有反射光线经过B点。综合上述结论，当线光源长度为某一个值2d时，取线光源在x点附近的微小长度dx。反射光在C点的光强和B点的光强如下其中，q表示单位长度线光源的功率，k为反射部分所占的比例，在d很小时近似为常数，即由结论（3）得其中c为系数常数，为给定位置为x的点光源产生的最大轨迹圆的半径。由此C点和B点的光强分别表示为：其中.下面确定函数。的解析解不可能用初等函数表示，故用数值积分的方法。如图4.2.3.3，距焦点为x的点光源
14、经过的交线圆上的点A映射到测试面上的点,点B映射到测试面上的点，点D映射到测试面上的点（根据对称性，和重合）图4.2.3.3 通过计算和两点的坐标，可确定轨迹圆半径。由数值计算可知，当时，恒有，所以问题一的解为1.56*2=3.12。4.3 问题二4.3.1 遍历搜索法根据4.1的推导，只要给定P以及抛物面上的任意一点就可以确定在测试面上相应的反射点，遍历抛物面上的所有点，求出测试面上的点，即可画出这一点的反射轮廓图，如下图所示：（a=1.65）图4.3.1.1在问题一得到的线光源的每点处，都画出反射轮廓图，迭加在一起边得到总的轮廓图如图4.3.2 图4.3.1.2对所得亮区轮廓图的几点分析：
15、1) 如下图为在光源点(1.56,0,15)处的反射测试图，其中“”所标为光源点通过椭球面的边沿反射后的点轨迹，可以发现，“”所围部分每一处都有两个重复的反射点，而“”外面则有四个重复的反射点，可见，外部光强是内部光强的两倍，同时可以看出“”轨迹是一个圆。2) 参看图4.3.3可以发现，整个轮廓完全由x坐标最大的光源点决定，即光源点偏离焦点越远，反射后的轮廓越大。图4.3.1.34.3.2切片叠加法为了画出测试屏上反射光的亮区，需要对光源经抛物面反射进行模拟。这里涉及到对整个旋转抛物面的搜索。按怎样的方式来遍历抛物面成了求解模型的关键。对于遍历球，我们不难想到有三种方式：1) 以线光源为轴旋转
16、切面，如图（1）。但存在问题是在I区域重复切。我们不作考虑。 图4.3.2.12) 图（2），沿FA轴，垂直FA切抛物面，所得为一个个交线圆。我们给出模型具体求解过程。求解前，先给出以下定义：定义 交线圆：垂直于Z轴的一个面与旋转抛物面的交线形成的圆入射光线、法线、反射光线在同一平面，于是我们在该平面上观察，化简成“求已知一直线关于另一直线对称的直线方程”的一个的数学模型。设线光源上任一点的空间坐标为（x0,y0,z0）,记为P0；在某一切面上的任意点坐标为（x1,y1,z1），记为P1；反射到检测屏上的点为（x2,y2,z2），记为P2。见下图4.3.2.2: 图4.3.2.2空间直线P1P
17、0向量： =（x0-x1,y0-y1,z0-z1），空间直线P1P2向量：（x2-x1,y2-y1,z2-z1），旋转抛物面方程 ：，记 ；过抛物面上P1点的法向量记为。由二次曲面法向量求解公式：（，）（x1,y1,z1）=（2x1,2y1,-60）。归一化：向量、模分别记为M1，M2。M1=| P1P0|=，M2=|P1P2|=。归一化后的向量分别为：，可得出方程式：利用Mathetatica强大的计算功能，可得出关于x2,y2一个表达式：当线光源上任一点P0向面镜照射，在经面上任一点P1反射后在测试屏上的位置p2可确定，记录该点。然后先在交线圆上遍历，再对不同半径的交线圆遍历的叠加、记录，
18、最终可得出反射光的区域，明暗程度也可得出。（主要编程算法见模型分析）。模型分析：圆切片模型的优点是直观上理容易解，理论上也容易推导出表达式（见上）。在Matlab 的编程中，通过对Z1值的控制，即步长的设定，一层层地遍历整个抛物面，相应的精度便可以控制。当Z1值知道，切面与旋转抛物面的交线为一确定圆：。那么，圆方程可表示为极坐标方程：， （，）通过对的控制，遍历整个交线圆。在对旋转抛物面遍历、线光源遍历时，本质上同对交线圆的遍历原理相同，也是通过设置步长来记录各个反射点。整个模型是个离散模型。理论上当步长趋于无穷小时，便可等效为连续模型。实际中并不可行。因为求解模型的算法概括为：点对点点对某一
19、交线圆交线圆叠加成面点到线光源。所编程序为一个三重循环。为追求精度，步长设得愈小，精度愈高，实际操作性能也就越低。举个例子，当三个步长均设为0.1，循环次数达到（72021636）。就是说时间复杂度相当高，为。这在实际中并不可行，建议这个模型做理论上的分析。为验证正确性，我们有选择的遍历了两个交线圆。交线圆半径=如图4.3.2.3，是线光源上在X0=1.83的点在经Z1=21.6,Z1=10处，即开口半径分别为36, 24.4949的两个交线圆反射后在测试屏上所得的点。图4.3.2.4是利用模型一中所画出的图。结果符合地相当好。 图4.3.2.3 图4.3.2.43) 平行于轴且垂直于屏地切，
20、即在图4.3.2.1中的图（1）上下平行切。这个算法其实跟方法2) 原理相同，通过层层遍历得到线光源在整个旋转抛物面的反射光在检测屏上的亮区。这里不再重复。五设计规范的讨论设计规范的合理性：5.1 测试屏放在25米处的合理性讨论首先讨论作为远光灯的测试标准时的合理性。人的平均反应时间t约为0.5秒，汽车在普通公路上行使的速度如果按照45千米/时计算。汽车紧急刹车的加速度为-5.0米/平方秒。当司机在夜间驾驶，突然发现前方有紧急情况，经过0.5秒后反应过来采取紧急刹车措施，直到汽车完全停止前进。整个过程汽车的前进路程为21米。测试屏放在25米处，保证此处的亮度，以便司机及时发现情况，并有足够时间
21、和空间采取措施，这样考虑还是合理的。但是随着公里的发展，路况越来越好，车速也可以越来越快。尤其是在高速公路上，车速一般可达80千米/小时，轿车可以达到100-120千米/小时。在这种情况下，在25米处测试显然已经不够了。按照80千米/小时计算，测试距离应该大于60米，而且对光强的要求也就相应的增加了。还有，如果仍然以AC=2AB=2.6米来测试，显然视角也变得太小，不利于司机发现情况。从这些角度考虑，这个方案的不合理性也就显现出来了。然而，如果这个标准作为汽车的近光灯的测试标准，还是比较合理的。近光灯的作用主要是照亮近处的路面状况以及路牌等,还有一个重要的作用就是照亮前面的车辆,以便司机在驾驶
22、时能够保持适合的车距.由于车辆在同方向行使时,相对车速比较小,因此,25米的距离已经足够.5.2 B、C两点的位置关系与光强度的意义。由问题二的结论得，在25米处，反射光亮区在水平方向的临界点在离对称轴大约2.6米处。视角约为12度。每辆车一般都有左右两个这种灯.所以,实际的视角还与车的宽度有关.车宽按照2米来计算,则这个视角实际上是大约52度.这个角度对于司机发现前方路面的状况已经足够了。按照单车道的宽度为7米计算,总共7.2米的亮区宽度刚好照亮本车道的区域.而且，这个宽度还有一个很好的性质就是不会影响到沿相对方向开过来的车辆的司机的视线这也就是为什么在夜间有车辆对开时,司机一般要关掉远光灯
23、,而只留下近光灯照明的原因.对于测试中要求C点的光强度不小于某一额定值,目的显然是为了保证司机能够看得足够清楚.本文建立的模型基于非常理想的情况：旋转抛物面发生全反射；光在空气中传播没有能量损失。事实上，当考虑损失因子时，灯丝的长度将增加。真正设计时，在考虑镜面介质材料、灯罩等影响时，得到损失因子ki(i，代表不同情况下的损失权重)。本题的情况，若要更高的精确度，可设两个权重系数，分别记k1 ,k2。k1：反射光的能量损失因子； k2：直射光的能量损失因子。根据一二问中建立的模型，我们只要分别在各自的总光强前乘以（1-ki）（ i=1,2）的系数即可。六参考文献：大学物理学 华金龙 王少杰 冯伟国 同济大学出版社应用光学 张以谟 机械工业出版社工科数学分析 张传义 包革军 张彪 哈尔滨工业大学数学系近代物理与高新技术物理基础 陈泽民 清华大学出版社线形代数与空间解析几何 俞南雁 韩瑞珠 周建华 科学出版社Mathematica 4.0实用教程 刘元高 刘耀儒 国防工业出版社Matlab实用教程 苏金明 阮沈勇 电子工业出版社专心-专注-专业