李 源, 张京国, 梁 谦
(中国空空导弹研究院,河南 洛阳471009)
近炸引信虚拟样机是在计算机上建立的表示引信系统(包括引战配合特性)及其使用环境的完整、详细的模型,它能够逼真地反映引信在测试、试验和作战使用中的性能,并且它能够用来评价引信的设计方案和评估产品的性能[1]。由于虚拟样机的以上优点,随着计算机技术的发展,近炸引信虚拟样机设计已经成为引信设计中的一个重要阶段。
作为激光引信的一个重要组成部分,激光引信虚拟样机必然要对光学系统进行建模。传统的光学仿真往往采用商业的光学软件来实现,但对于激光引信而言,其采用的光学系统形式不同于传统光学系统,特别是用来进行光电转换的探测器往往具有不规则面形,而光学软件只能针对光学系统的某个功能进行仿真以及优化,对目标特性以及不规则面形探测器都不能进行仿真。基于此种原因,本文采用统计建模与光学追迹的方法,利用Matlab软件编程对此光学系统进行虚拟样机建模。
发射光学系统子午面视场角α;
发射光学系统弧矢面视场角β;
激光源:量子阱阵列GaAs半导体激光器。
图1所示为光学发射机的组成示意图。
图1 光学发射机的组成示意图
按照发射光学系统要求,对激光器发出的激光束进行整形设计[2],如图2所示。
图2 发射光学系统
激光器发出的激光经过光学系统整形后,在子午方向光束能量P(α)是高斯分布;在弧矢方向由于激光器发光结是阵列分布,且又经光学系统扩束,所以光束能量P(β)在此方向即不是高斯分布,又不是理想的均匀分布。这就给光学发射机的激光光束建模带来了难度。所以,在本文中借助试验参数,运用统计建模法将多个测试结果进行融合,最终得到发射光学系统的能量空间分布。
图3是光学发射机在子午面以及弧矢面的能量分布示意图。
图3 光学发射机在子午面和弧矢面的能量分布示意图
上述情况仅考虑了光束能量的空间分布,而在本文中还要考虑激光束功率随时间的变化。本文所研究的激光引信采用的是脉冲工作体制,直接向半导体激光器注入电流,改变发射波形。激光光束功率P(t)随时间的变化如图4所示。
图4 光学发射机光束能量随时间分布示意图
假设由光学系统光面发射到外部空间的功率为PS,将此功率在子午面与弧矢面内进行离散化,离散化后的单位发射光束与目标反射面的截面称为单位反射面元,如图5、图6所示。
由上图可知,单位面元面积 ATAR=
图5 光学发射机在子午面的能量均分示意图
图6 光学发射机在弧矢面的能量均分示意图
(π/4)R2θαθβ,其中R 为发射机到目标间的距离;θα、θβ分别为单位发射光束在子午面以及弧矢面的角宽度。照射到单位面元上的功率PTAR为
式中:A(α,β,t)为光束能量PS随子午面视场角、弧矢面视场角以及时间变化的系数;τ 为发射机到目标间的大气透射比。
由目标单位反射面元反射到半球空间内的回波功率PSP为
式中:ρ为目标的反射比;θ为探测方向与朗伯面表面法线的夹角。离反射源距离为R 的空间点上的辐照度HR 由下式表示:
接收光学系统子午面视场角α;
接收光学系统弧矢面视场角β;
接收口径L×H。
本文所研究的是以二次曲柱面反射镜为主体形式的光学系统,如图7所示。
图7 接收光学系统光线追迹图
入射光线经过接收光学系统聚焦成像在反射镜的焦平面上,形成与反射镜口径等宽度的一段线像,图8是成像光斑的能量分布图。
为了减小探测器的动态范围,本光学系统采用的是异形探测器。由图7可知,入射光线经过光学系统在探测器上的成像位置与其在接收口径上的入射位置有关,由于所需建模的接收光学系统采用的是不规则面形的探测器,这就使得光学设计软件不能进行精确的成像光斑与探测器相对位置关系的仿真。基于以上原因,此接收光学系统采用Matlab软件进行建模。
图8 接收光学系统成像光斑能量分布图
将引信接收窗口在两正交方向离散化,均分为多个单位接收面元,如图9所示。
图9 引信窗口接收面元与目标反射面元的关系示意图
引信窗口单位接收面元AW上接收到的回波功率PW由下式表示:
将公式(1)、(2)、(3)、(4)、(5)整理后得到
由图7可以看出接收光学系统的入射光线是一个在某一时刻并携带了光束能量的空间矢量→S(PW,t),它是一个多维变量,将此矢量在接收光学系统的子午面、弧矢面投影,得到传统定义的子午面入射角α、弧矢面入射角β。在进行光学系统的建模时,可以分别在这两个投影面上进行光线追迹,得到成像光斑位置。
假设在接收窗口点A 处,且子午面入射角为3°,弧矢面入射角为45°时,经过仿真可以精确的判断入射光线经过光学系统后是否成像在探测器的光敏面上,如图10所示。
如果成像光斑落在探测器的光敏面上,根据PIN 硅光电二极管的光电转换特性[3],可以得到
探测器的输出电流信号及其大小
图10 反射光束经过光学系统后成像在探测器上的位置示意图
式中:A 是探测器的响应度,单位为A/W;A0是接收光学系统总效率。由于现阶段引信的激光回波能量一般为几个μw ~几十个μw,完全在探测器的线性响应范围内,所以就可以将某一时刻所有有效的成像光斑所产生的信号线性迭加,得
由上式可以得到在某一时刻t探测器的输出电流信号It,将此信号在时间轴上展开计算得到I(t),它就是探测器与一个发射脉冲相对应的输出电流随时间变化的曲线。
激光引信的光学系统形式以及所要探测的目标特性都较为复杂,本文采用的这种统计建模与Matlab软件编程相结合的光学系统虚拟样机建模方法合理有效,从最终的仿真结果来看,比较符合物理样机的现实情况,真实的反应了光学系统的特性,满足激光引信虚拟样机的要求。
[1] 张京国,周宗海,刘建新.近炸引信虚拟样机技术研究[J].制导与引信,2006,(3).
[2] 赵为党,杨李茗,等.激光光束整形的设计和研究[J].红外与激光工程,1999,28(1).
[3] 安毓英,曾晓东.光电探测原理[M].西安:西安电子科技大学出版社,2004.