车灯线光源的计算.docx
- 文档编号:28472349
- 上传时间:2023-07-14
- 格式:DOCX
- 页数:13
- 大小:316.66KB
车灯线光源的计算.docx
《车灯线光源的计算.docx》由会员分享,可在线阅读,更多相关《车灯线光源的计算.docx(13页珍藏版)》请在冰豆网上搜索。
车灯线光源的计算
第十讲车灯线光源的计算
2002年高教社杯全国大学生数学建模竞赛大专组C题为:
安装在汽车头部的车灯的形状为一旋转抛物面,车灯的对称轴水平地指向正前方,其开口半径36毫米,深度21.6毫米。
经过车灯的焦点,在与对称轴相垂直的水平方向,对称地放置长度为4毫米的线光源,线光源均匀分布。
在焦点F正前方25米处的A点放置一测试屏,屏与FA垂直。
请解决下列问题:
1.计算直射光总功率与反射光总功率之比。
2.计算测试屏上直射光的亮区,在有标尺的坐标系中画出其图形。
3.计算测试屏上反射光的亮区,在有标尺的坐标系中画出其图形(只须考虑一次反射)。
一、建模的基本思路
1.在空间以旋转抛物面的顶点为坐标原点,对称轴为OZ轴,使OX轴平行于线光源,建立空间直角坐标系。
2.解决问题1时,由于灯丝长度与车灯抛物面的焦距相比很小,可以看成点光源。
以光源为球心,以光源到开口圆的最大距离为半径,作一球面,直射光照到的球面面积与反射光照到的球面面积之比,即为功率比。
3.结合图形和对题目的分析,在假设理想的条件下,可得出“运动”点光源上发出的光线经过开口圆的任一个直径在测试屏上直射所成的直线与开口圆的直径有一个恒定的比值。
所以直射光在屏上成的像和开口圆相似的叠加图形。
确定测试屏上的亮区为类似于标准跑道的图形。
4.对线光源上的任一点,在旋转抛物面上任取一点,根据物理光学反射定律,入射角等于反射角,入射光线和反射光线及法线在同一平面内,从而得出反射线的方程,计算出反射线与测试屏的交点,借助计算机编程,将这些点画出来即为测试屏上反射光的亮区。
二、算法与结果
取旋转抛物面的顶点为坐标原点、对称轴为OZ,使OX轴平行于线光源AB,建立空间直角坐标系如图。
则旋转抛物面
,即
x==36cosi
卩=36winf
z=21.6
L.
1直射光与反射光的功率比
按照对题目的分析,将点光源都放置于球心位置。
它直射到球上的面积显然为一球冠的表面积。
将直射面积与反射面积作比。
这样就可将功率比问题转化为球缺与球冠表面积之比。
由球冠表面积公式可得,其功率比=表面积比二球缺与球冠高比,问题就简单化,下
面就转化为球缺与球冠的高的比。
再对图形进行分析,在求高的过程中可以将其转化为平面图形来处理。
我们以焦点为圆心。
以焦点到开口圆的最大距离为半径做球面如图所示:
003(3±1..0
卩=2严刑直射区球冠的高=0「一门亘=门订I—"门
1—cosa
而另一部分球冠高二=丁1二丁,很显然它们之比等于「十」.厂,
即直射光总功率厂与反射光总功率二之比,
P\1—cosa
■1-.••。
可将线光源看做是由一个点光源在兰[线上做上下重复运动。
⑴当点光源移动到B时,即B为点光源。
显然BN为圆的半径。
角亡,在三
角形中已知三边长,用余弦定理可得出】,从而可求出二山值。
Pl—=0.693694346
通过计算可得,二二.im王,」i,为其中一个值。
(2)同理,当点光源移动到M时,结果同上
Pi
(3)当点光源移动到焦点
1—C03CI
Pr1十eosa
=0.654444444
。
21.6-15
F时
=0.180327868
通过对点光源依次移动情形的计算比较,可得它们的总功率比值在
[0.6936948460.694444444],误差很小。
所以,最终我们对这个区间均匀求值。
得出最终结果:
P1
&=0.694069645。
2.直射光亮区
因为光线上任一点光源到开口处的平面K的距离恒为(21.6—15)毫米;到测试屏的距离恒为2500毫米。
则开口处的图形(即圆)应与该点光源在测试屏上的直射图形相似,也为圆,且各个点光源所形成的直射圆的半径相等。
由此可知,从最上边的点光源开始最下边点光源结束,直射区应是这些点光源形成的直射区的圆的迭加相当于测试屏上有半径相等的圆从上而下连续地推移。
所以,形
成直射区图形应为类似于操场的形状,两端各为一个半圆,中间为矩形且矩形的一组边与线光源平行。
在面K中,开口处圆的直径2r,任一点光源在测试屏上圆的直径2R,则
姑_21.6-15
因为2r为36毫米,所以可解得2R=272727.2727毫米,
点光源在最高点时,设L为直射区下限点的|x|,
36十2_21.6-16
L+2=25000
2L=287874.7879毫米。
因为形成图形为对称图形,所以在直射亮区中从最高点到最低点的距离为2L
即287874.7879毫米。
所以中间矩形部分的与线光源平行的那一组边缘线长度为
(2L-2R)=15147.751518毫米。
下图为各个点光源经过开口圆上的点在测试屏上所形成的图形的叠加
下图为直射亮区的边缘线
3反射光亮区
在线光源上任选一点(m,0,15),-'r'-,在反射面上任一点为
异
(咖祖胪就丽)。
此点处的内法向量为(-十如,-曲祖旳,设反射光线与测试屏的交点为(x,y,25015),则反射光线与入射光线以及法线应满足以下条件。
(1)入射角二等于反射角;,则二二---',此时入射光线的单位向量为
彳m—TCQSt.——二卜
I60丿
(m—rcost)2+(rs沁十(15——)2
法向量的单位向量为
—s{—rcost9—F曲说,30}
反射向量的单位向量为
r
x—rco3tfy—r5s^25015——
(x—rcost)2+@—rsrnt)2十(25015—二)
O1.1'
于是,閃呂a=盲■它,凸0昭=b,所以才•存=7T・b,
由此可得方程
“;厂I.i]I.-o
(2)反射光线与法线以及入射光线在同一平面内。
因为三条线有一个公共交点,所以只需其中一条直线的向量与另外两条直线所组成平面的法向量垂直即可。
其中入射光线与法线组成平面的法向量为
_宀。
因为「与】'垂直,所以:
「•】一】,即…°-.八一。
由此可
得方程-讥h
由方程组中的第一个方程可以得到两个y值,其中一个必须舍去,
原则为使得二1一的值较大的一个舍去此处丨,分别为
r3
—cos/+15rcos/-30m
60
-15rsm£-—smZ
60
. J-: ;叮.■■-.■■■: : i -r2-ru^cosf+3025015-— 2I60丿 应用计算机编程(程序如下页)求解反射光的亮区。 在编程时应附加一个限定条件,因为在抛物面的开口圆处的光线部分会反射,也有部分贴紧边缘直射出去,所以每条发射线都会对应地在测试屏上有两个光点,而其中必有一条是直射出去的。 因为直射出去的光线必定比经过反射的光线相对于测试屏更为发散, 所以在求解出来的两个解应取I/的值较小的一个作为反射光线在测试屏上形成的点,满足这所有条件的点才可以组成反射亮区 问题的讨论与分析 对于车灯线光源的优化模型的设计与解答,我们已经尽可能的考虑最一般的情况,使答案得以尽量地接近真实的情况。 在建立普通的数学函数模型基础上,我们借助计算机编程来描述函数模型的轨迹。 我们也注意到,在问题一中对功率的损失理想化,这样做比较方便于建立模型和计算操作。 但事实上光在传播、反射中都必然会有能量损失,只不过在正常情形其损失较小,因而所建立的模型存在不太大的偏差。 问题二中运用到线段的等比关系,使问题直观化,便于计算,并且使得计算结果比较精确,模型建立和求解都相当成功。 问题三中主要是对最一般情况的考虑,使结果更接近于准确值,借助于计算机绘图,得出反射光亮区。 模型的进一步优化 在此题中,主要讨论线光源长短对直射亮区,反射亮区的影响,各种等能量 的光源上,当光源长度变化时,直射区和反射区都会发生变化。 但经过分析光源长度的变化对直射区的影响要小于反射区,所以就主要考虑反射区。 当光源长度增长时,反射区就会变大,自然,光的强度就会弱一点,光就发更发散,而当光源长度减小时,反射区就会减小,光的强度大一点,光更集中。 就汽车而言,当在高速路上行驶时建议使用尽量短的光源,使光线尽可能地集中,有利于极早发现远方物体。 而当是在市区的行驶的车量,由于人流比较密集,所以也应使光线较强,应用短光源。 而在郊区时对光强度要求不高,在达到相同的光强下,用点光源则可以节省能源,对照明时间长且光强要求不同的情况都可以使用长光源。 ■3000- -4000- 灯丝長3毫世 〔*wn卧琨A 60004000 20000 k轴[单位亳米) 2CC0 40006000 2000 1600 1000 600 D -600 100D -1500 ■2000 灯丝长2毫米 ■3000”2000-1000010002000 x轴伸位毫米) 3000 程序如下: forr=[5,10,15,20,25,29,33,36] form=linspace(-2,2,10); fort=linspace(0,2*pi,500); a1=(-1)*(r*r*r*cos(t)/60+15*r*cos(t)-30*m)/(-15*r*sin(t)-r*r*r*sin(t)/60); a2=(24985*r*m*sin(t)-rA3*m*sin(t)/60)/(-15*r*sin(t)-r*r*r*sin(t)/60);a3=(-r*cos(t)*(m-r*cos(t))+r*r*sin(t)*sin(t)+(15-r*r/60)*30)/sqrt((m-r*cos(t))*(m-r*cos(t))+r*r*sin(t)*sin(t)+(15-r*r/60)*(15-r*r/60));a4=a3*a3*(a1*a1+1); a5=a3*a3*(2*a1*a2-2*a1*r*cos(t)-2*r*sin(t));a6=a3*a3*(a2*a2-2*r*a2*cos(t)+r*r+(25015-r*r/60)*(25015-r*r/60));a7=-(a1*r*cos(t)+r*sin(t)); a8=r*r-r*a2*cos(t)+30*(25015-r*r/60); a9=a4-a7*a7; a10=a5-2*a7*a8; a11=a6-a8*a8; y1=(-a10+sqrt(a10*a10-4*a9*a11))/2/a9; x1=a1*y1+a2; y2=(-a10-sqrt(a10*a10-4*a9*a11))/2/a9; x2=a1*y2+a2; kl=x1A2+y1A2-x2A2-y2A2; ifkl<0 plot(x1,y1,'r.') holdonelseplot(x2,y2,'r.')holdonendendendendgridonaxisequal
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 车灯 光源 计算