气动热辐射对超音速状态下的导引头红外成像的影响分析 下载: 634次
1 引言
21世纪以来世界**革命迅猛发展,信息主导成为制胜关键,而导弹作为现代战争的关键**,具备“发现即摧毁”的能力。在追踪锁定目标这一过程中,导弹中的制导系统扮演着极其重要的角色。
导弹通常由战斗部、弹体结构、动力装置和制导系统组成。导引头是安装在导弹头部的制导装置,其内部安装有红外光学成像系统,通过接收探测目标的红外辐射来获取目标信息。不同于其他的探测方法,红外成像制导技术的抗干扰能力、隐蔽能力及夜间能力更加优秀,然而导引头在大气层中以超音速的速度飞行时其光学窗口会因气动加热提高自身的温度,从而产生强烈的红外辐射,形成辐射干扰[1]且光学窗口作为光学探测系统唯一的探测通道,使得探测器成像对于光学窗口的温度变化尤为敏感,气动热辐射严重时会让光电探测器饱和而不能接收来自目标的辐射,这种效应称为气动热效应,会极大地影响导引头的制导能力[2-4]。
以往在气动热方向上被研究的对象大多是头部圆形整流罩,研究重点在于整流罩的外形特性及头罩折射率场的变化等对成像质量的影响,或是气动热过程中防护罩与光学系统的相互作用,研究方法多为数值仿真[5-10]。本文的研究对象为导引头侧面红外光学窗口而非位于头部的光学窗口,侧面光学窗口的设计可以有效降低冷却系统的设计难度,减少气动热带来的影响,提高红外光电探测器成像的质量[11-12]。研究了光学窗口在导引头超音速飞行状态下温度变化,仿真了受光学窗口干扰的探测面辐照度[13],推导出了信噪比随光学窗口温度和探测距离变化的关系[14],最终分析得到了气动热辐射对红外光学系统成像影响的规律。
2 光学窗口温度场仿真
2.1 导引头模型的建立
导引头模型如
根据模型尺寸在ANSYS软件中建立模型,导引头的几何模型总长为510 mm,模型半径为105 mm,内置有红外探测镜头,与窗口法线方向成65°角,光学窗口材料选用蓝宝石。导引头网格模型如
图 2. 导引头网格模型及其流场网格模型。(a)网格模型;(b)流场网格模型
Fig. 2. Seeker mesh model and its flow field grid model. (a) Mesh model; (b) flow field grid model
2.2 光学窗口温度场仿真结果
当飞行器以较低速度飞行时,气动热导致的光学窗口升温较小,且光学窗口的红外辐射率很低,所以光学窗口在探测器上产生的辐射照度也很低,相对于目标红外辐射可以忽略。故本研究中导引头飞行速度为3 Ma(1 Ma≈340.3 m/s),飞行高度分别为2、11、20 km,飞行攻角为0°。Ansys-Fluent仿真外流场边界条件[15]具体参数如
表 1. 外流场边界条件
Table 1. Boundary conditions of external flow field
|
在大气层中高速飞行的导引头外流场是较为复杂的湍流流场,求解器选择三维双精度耦合瞬态密度基求解器。湍流模型选择专门为高速航天航空设计的spalart-allmaras模型,该模型对于高速飞行器外流场具有优秀的边界层计算能力。外流场气体设置为理想气体,对于高速可压缩流动气体采用描述气体粘度的Sutherland定律。壁面条件设置无滑移,表面粗糙度为0.5[16]。Fluent计算时设置时步数量为100,时间步长为0.1 s,即计算飞行时间10 s内的温度场变化,仿真时对Fluent进行的相关求解设置如
表 2. 光学系统参数
Table 2. Optical system parameters
|
图 3. 导引头飞行10 s不同海拔光学窗口温度场分布。(a)2 km;(b)11 km;(c)20 km
Fig. 3. Temperature map of seeker at different altitudes after flying for 10 s. (a) 2 km; (b) 11 km; (c) 20 km
图 4. 光学窗口最高温度随飞行时间变化情况
Fig. 4. Maximum temperature of optical window varying with time of flight
3 光学窗口热辐射对红外成像系统的影响
3.1 光学窗口热辐射计算
物体温度高于绝对零度时都将辐射出能量,该现象称为热辐射。对于绝对黑体,其光谱辐射出射度
式中:c1为第1辐射系数,c1=3.741832×10-12 W·cm²;c2为第2辐射系数,c2=1.4388 K·cm;λ为波长;T为黑体温度。
为了描述一般不具备绝对黑体特性的物体辐射,引入“辐射发射率”
式中:
窗口辐射亮度为
式中:
则光学窗口在探测系统入瞳上的辐射照度
式中:
3.2 光学窗口辐照度仿真
本研究基于TracePro软件进行探测器光学窗口辐照度仿真,TracePro可以根据不同辐射的特性设置相对应的辐射源、建立光学系统实体模型及探测面[17]。
根据导引头实际情况,选用蓝宝石晶体作为光学窗口辐射源材料,蓝宝石晶体在3~5 μm红外波段有良好的透过率,且在高温下有低辐射率特性,满足红外成像系统的工作要求。本研究中光学系统为透射式中波红外光学系统,光学系统及探测器参数如
在TracePro软件中将发射形式设为灰体辐射,光谱范围设为3~5 μm。当光学窗口作为光源时可以近似看作朗伯光源,遵循朗伯余弦定律,故发光场型采用朗伯体发光场型。依据第2小节导引头温度场的计算结果设置光学窗口的温度及对应的发射率。依据光学窗口温度场仿真结果,计算得到探测面上不同时间下辐照度的仿真结果,导引头在海拔2、11、20 km高度飞行时间10 s时的探测面的辐照度分布如
图 5. 导引头飞行10 s探测器辐照度仿真图。(a)2 km海拔探测器辐照度图;(b)11 km海拔探测器辐照度图;(c)20 km海拔探测器辐照度图
Fig. 5. Simulation diagrams of detector irradiance for 10 s flight of seeker. (a) 2 km altitude detector irradiance map; (b) 11 km detector irradiance map; (c) 20 km detector irradiance map
依据TracePro软件仿真光学窗口的不同温度下探测面的辐照度结果,可得到导引头在不同海拔高度飞行10 s内探测器面辐照度的变化趋势,结果如
图 6. 导引头飞行10 s内探测器接收面最大辐照度变化趋势
Fig. 6. Change trend of maximum irradiance of detector's receiving surface within 10 s of seeker's flight
为了更加直观地看出光学窗口热辐射对成像系统的影响,仿真得到了不同海拔高度下飞行10 s时的失真目标图像,其光学窗口温度分别为300、305、318、376 K,如
3.3 光学窗口热辐射对成像系统信噪比的影响
通过第3.2节的失真目标仿真图像可以看出,随着光学窗口的气动加热温度不断升高,热辐射噪声不断加大,红外成像系统目标图像质量愈加下降,目标信号相比于热辐射噪声不断减小,则可以通过目标信号幅值与热辐射噪声的幅值的比值——信噪比(SNR),来描述红外光学系统的成像质量。
SNR为目标信号输出功率与同时输出的噪声信号的比值[18],即
真实的目标辐射计算是相当复杂的,它包括目标体各种热源及其热传导辐射计算、飞行过程中的气动热辐射计算、太阳及地球辐射照射及反射计算等,需要建立复杂的数学模型来计算模拟,并辅以实际测量数据来验证。为了工程计算方便,有时也将目标等效为一个灰体辐射元,设其温度为
通过辐射出射度可以得到目标的辐射强度和辐射亮度,目标的辐射强度为
目标的辐射亮度为
设探测器的接收系统置于目标辐射轴的垂面上,距目标的相对距离是
式中:
又因为入射到单位接收表面上的辐射功率
式中:
由
则在红外光学系统接收孔径
则光学窗口热辐射产生的红外噪声为
式中:光学系统接收孔径面积
式中:
根据上述公式,整理可以得到
3.4 信噪比结果分析
控制其他变量,仅研究光学窗口温度与目标相对距离特性对信噪比的影响。在信噪比模型中,目标物体选择美军战斗机,等效面积为78 m2,光学系统接收孔径约为380 mm2,探测器面积为58.9 mm2,饱和温度为400 K。经信噪比模型计算可得到光学窗口在不同飞行时刻的不同温度下红外光学系统信噪比变化情况,如
导引头在3 Ma的状态下在2 km海拔飞行10 s时对应窗口温度为376.55 K,信噪比降低到3.23,而在11 km和20 km海拔飞行10 s时对应光学窗口为318.66、305.85 K,信噪比则分别是19.618、31.17。即不同海拔高度对于信噪比影响较大,这是因为海拔越低,受气动热影响越大,光学窗口温度越高,故信噪比越低。当光学窗口温度趋于稳定,信噪比也会随之稳定。当光学温度升至340 K左右时,信噪比降至10以下,温度升至380 K左右时,信噪比下降到3以下,红外成像系统探测器被严重干扰。从
若光学窗口经制冷设备被保持在一个稳定的温度,根据信噪比模型可以得出,距离目标越近,信噪比越高。本研究中当光学窗口实际温度被稳定在300 K时,信噪比与探测距离的关系如
4 结论
主要研究对象是导引头侧面光学窗口,采用有限元分析法对导引头口进行温度场的模拟,得到光学窗口的温度变化情况并通过TracePro软件仿真探测器探测面的辐照度。建立了成像系统信噪比模型,分析得出导引头飞行海拔高度、目标探测距离和光学窗口温度与信噪比之间的关系。研究结果表明:导引头飞行海拔高度很大程度上决定红外成像系统受气动热干扰的程度,当导引头的飞行海拔为11、20 km,在飞行10 s后,光学窗口干扰的探测器探测面最大辐照度相比于2 km海拔降低了83%、89%,信噪比增加了6.07倍、9.65倍。海拔越高,辐照度越低,则成像系统信噪比越高,成像质量越好。当光学窗口温度上升至探测器饱和温度后,信噪比下降到1以下时就应当采用制冷方式降低光学窗口辐射影响来提高红外光学系统成像质量。
[1] Cross E. Window-heating effects on airborne infrared system calibration[J]. Proceedings of SPIE, 1993, 1762: 1-8.
[2] 殷兴良. 气动光学原理[M]. 北京: 中国宇航出版社, 2003.
YinX L. Principle of Aerooptics[M]. Beijing: China Astronautic Publishing House, 2003.
[3] 王惠, 党凡阳, 张荣达, 等. 半球形整流罩厚度对气动热辐射的影响[J]. 红外与激光工程, 2018, 47(12): 153-158.
Wang H, Dang F Y, Zhang R D, et al. Influence of thickness of hemispherical dome on aerodynamic thermal radiation[J]. Infrared and Laser Engineering, 2018, 47(12): 153-158.
[4] 张发强, 张维光, 万文博. 基于光线追迹的红外探测光学系统杂散辐射研究[J]. 红外与激光工程, 2019, 48(9): 0904006.
[5] 史要涛, 王玉雷, 解放, 等. 基于Code V的气动热环境下头罩的光传输仿真[J]. 航空兵器, 2015, 22(5): 41-45.
Shi Y T, Wang Y L, Xie F, et al. Optical transmission simulation of dome in aero-dynamic heating environment based on CodeV[J]. Aero Weaponry, 2015, 22(5): 41-45.
[6] 薛文慧. 气动热环境下共形整流罩热辐射特性研究[D]. 哈尔滨: 哈尔滨工业大学, 2016.
XueW H. Research on the thermal radiation characteristics of conformal dome in the aero-dynamic environment[D]. Harbin: Harbin Institute of Technology, 2016.
[7] 吕银环. 防护罩气动热效应对玫瑰扫描系统的影响仿真分析[J]. 红外与激光工程, 2020, 49(7): 20190520.
[8] 董玉磊, 韦承甫, 刘现魁, 等. 狭缝式自由旋涡气动窗口光学质量测量方法研究[J]. 中国激光, 2021, 48(23): 2304003.
[9] 陈勇, 谢伟明, 路大举, 等. 吊窗尾迹流动非稳态气动光学效应研究[J]. 光学学报, 2020, 40(16): 1601001.
[10] 张丽琴, 费锦东. 高速飞行器成像探测气动光学效应研究(特约)[J]. 红外与激光工程, 2020, 49(6): 20201016.
[11] 易仕和, 丁浩林. 稠密大气中高超声速导引头红外成像面临的机遇、挑战与对策[J]. 现代防御技术, 2020, 48(3): 1-10.
Yi S H, Ding H L. Opportunities, challenges and countermeasures for infrared imaging of hypersonic seeker flying in condensed atmosphere[J]. Modern Defence Technology, 2020, 48(3): 1-10.
[12] 王亚辉, 王强, 张伯川, 等. 红外窗口气动热辐射效应评估方法[J]. 红外与激光工程, 2016, 45(2): 0204001.
[13] 金和, 任栖锋, 陈为. 窗口自发辐射对于红外探测器系统性能的影响[J]. 光电工程, 2016, 43(8): 84-88.
[14] 徐振亚, 付奎生, 祁鸣, 等. 一种用于发射后截获的红外导引头探测距离估算方法[J]. 红外技术, 2020, 42(11): 1095-1102.
[15] 姜波, 齐杏林, 贾波, 等. 基于fluent的弹丸外流场仿真计算[J]. 计算机仿真, 2014, 31(3): 38-40, 45.
Jiang B, Qi X L, Jia B, et al. Simulation research on projectile external flow with fluent[J]. Computer Simulation, 2014, 31(3): 38-40, 45.
[16] 赵璇, 孙智, 吕润民, 等. 不同热边界对超音速机翼气动热数值分析的影响[J]. 航空工程进展, 2019, 10(5): 714-721.
Zhao X, Sun Z, Lü R M, et al. Influence of different thermal boundary on numerical analysis of aerodynamic heat of supersonic wing[J]. Advances in Aeronautical Science and Engineering, 2019, 10(5): 714-721.
[17] 孙成明, 赵飞, 张泽. 基于TracePro的空间目标光学散射特性建模与仿真[J]. 光子学报, 2014, 43(11): 1122003.
[18] 李振. 预警机蒙皮红外特征与成像仿真技术研究[D]. 成都: 电子科技大学, 2016.
LiZ. Research on infrared characteristics and imaging simulation technology for AWACS-skin[D]. Chengdu: University of Electronic Science and Technology of China, 2016.
Article Outline
黄振, 江伦, 胡赫, 张明, 李奇, 宋延嵩, 董科研. 气动热辐射对超音速状态下的导引头红外成像的影响分析[J]. 激光与光电子学进展, 2023, 60(2): 0201001. Zhen Huang, Lun Jiang, He Hu, Ming Zhang, Qi Li, Yansong Song, Keyan Dong. Analysis of Influence of Aerodynamic Heat Radiation on Infrared Imaging of Seeker in Supersonic State[J]. Laser & Optoelectronics Progress, 2023, 60(2): 0201001.