Real-time infrared detection imaging simulation of aircraft considering influence of solar radiation
-
摘要: 针对飞行器红外探测过程的成像仿真,有助于红外隐身效能的评估,进而提高应对红外探测的反侦察能力。建立了飞机红外成像仿真模型,并针对不同的飞行姿态、速度以及不同的探测时间,求解其蒙皮温度与红外辐射亮度,结果表明当飞机从低速到高速运行时,太阳辐射条件与气动加热层温度分布分别成为影响其成像结果的主要因素。建模方法及仿真结果对飞行器的隐身设计具有指导意义,为飞机实时红外成像场景的生成提供理论基础Abstract: The imaging simulation for the infrared detection process of the aircraft contributes to the evaluation of the infrared stealth performance, thereby improving the anti-reconnaissance ability to deal with infrared detection. The aircraft infrared imaging simulation model was established, and the skin temperature and infrared radiation brightness for different flight attitudes, speeds and detection moments were solved. The results show that when the aircraft runs from low speed to high speed, the solar radiation conditions and the temperature distribution of the aerodynamic heating layer become the main factors affecting the imaging results. The modeling method and simulation results proposed have guiding significance for the stealth design of the aircraft, and provide a theoretical basis for the generation of real-time infrared imaging scenes.
-
Keywords:
- infrared detection /
- imaging simulation /
- radiation brightness /
- real-time
-
引言
微光夜视仪是在夜间或低照度条件下,利用星光、月光和大气辉光,通过像增强器,将人眼不易看见的极微弱星光和红外辐射等进行光电转换和图像增强,使之变成人眼容易看到的图像,以实现夜间观察和探测[1]。高品质的微光夜视装备是夜间作战的眼睛和窗口,是打赢一场高新技术局部战争的坚强后盾和有力保障之一[2-5]。
微光夜视仪主要由物镜、像增强器、目镜等三部分组成,其中实现微光夜视功能的核心部件是微光像增强器,它是一种真空光电成像器件。光纤倒像器作为像增强器上的光学输入、输出窗口,是影响其成像品质的核心器件[6-8]。目前国内外生产和应用的微光像增强器,主要是以直径Φ18 mm的口径为主,随着微光夜视技术的快速发展和系统性能的不断提高,大口径光纤倒像器因其可实现宽视场、大视野、远视距的探测和观瞄,在未来信息化和高科技战争中具有广阔的应用前景和发展潜力[9-20]。但是随着光纤倒像器口径的增大,易产生边缘分辨率消失、内部斑点、弯曲畸变等缺陷,现有的制备工艺已不再适用,为此,研究人员做了大量的研究工作,对大口径光纤倒像器的制备工艺进行了改进,最终制备出满足技术要求的大口径光纤倒像器。
本文从大口径光纤倒像器的制备技术出发,开展了大口径光纤倒像器关键工艺技术研究,最终解决了大口径光纤倒像器的制备难点,满足了产品性能指标要求。
1 大口径光纤倒像器制备工艺设计
1.1 大口径光纤倒像器的设计要求及制备难点分析
大口径光纤倒像器的主要技术指标如表1所示。其外径大于45 mm,高度大于50 mm,在尺寸和体积增大的同时,其他性能指标与常规倒像器产品相同[21-22]。
表 1 大口径光纤倒像器设计指标Table 1. Design indexes of optical fiber image inverters with large aperture项目 技术指标 有效区直径Φ/mm >42 纤维丝径/μm ≤6.0 热膨胀系数/℃ (87±2)×10−7 数值孔径 >0.65 透过率/% >65 分辨率/(lp·mm−1) >110 剪切畸变/μm ≤25 弯曲畸变/μm ≤30 像位移/μm ≤50 扭转角度/(°) 180±1 放大率 1.0±0.01 斑点尺寸或鸡丝长度/μm <90 如图1所示,我们将大口径光纤倒像器按其外形特征划分为两部分,即直区(纤维未变形区域,包括底部直区和台阶直区)和扭区(纤维扭转变形区域),其中直区要求尽可能地保证不出现弯曲畸变,受热尽量少;而扭区要求纤维变形趋于一致,受热尽量均匀,这就对制备工艺提出了很高的要求,需要充分利用大口径光纤倒像器的形体特征来完成制备过程。
图2所示为本文研究的大口径光纤倒像器的扭转变形示意图。图3所示为光学纤维扭转拉伸变形示意图[4, 23]。大口径光纤倒像器的扭转成型是制备过程的关键工序,其中50%以上的性能指标都与扭转成型工序直接相关,如对比度、分辨率、透过率、弯曲畸变等。这是由于大口径光纤倒像器毛坯外形体积是常规光纤倒像器所用毛坯体积的13.8倍(常规光纤倒像器所用毛坯直径25 mm、高度25 mm),使得毛坯在扭转成型工艺中在高温阶段的加热时间增加了5倍,最外围光学纤维的拉伸量增加了80%。这一方面加重了光学纤维芯料、皮层料界面的离子扩散和渗透,导致界面上出现析晶甚至乳化现象;另一方面大口径光纤板毛坯在受热扭转成型过程中,扭区的光学纤维会被拉伸,皮层玻璃材料的厚度变薄,距离光纤倒像器中心轴距离越远的边缘光学纤维,其纤维的拉伸量越大,使得边缘光学纤维丝的皮层厚度变得越薄。当皮层厚度小于(0.5~1)λ(λ为入射光波长)时,光线就会穿透皮层进入到相邻的纤维中去,造成漏光、串光。为了防止产生纤维的漏光、串光,需要增加皮层厚度,但是皮层厚度过厚,会使得入射在皮层上的无用光通量增加,增加背景噪声,造成像质下降[24-26]。这两方面的共同作用破坏了芯、皮界面的完整性,使部分光线不满足全反射条件从光学纤维中逸出,从而产生斑点、数值孔径降低、透过率降低、分辨率下降等问题;同时随着毛坯尺寸的加大,加重了毛坯受热不均,导致扭转过程易产生炸裂。因此,扭转成型是影响大口径光纤倒像器传像品质的关键工序。
1.2 大口径光纤板毛坯的制备工艺设计
大口径光纤倒像器是由成千上万根微米级的光学纤维平行规则排列后,经热熔压成型制备成光纤板毛坯,再经过外力扭转毛坯,使输入图像反转180°从输出端传出的一种光纤传像元器件。其中每一根光学纤维都是一个独立的信息传输单元,它是由高折射率的芯料玻璃和低折射率的皮层玻璃构成,光线能够按照全反射原理在光学纤维内部传播。根据纤维光学的理论和光的全反射定理可知[24-26],当一束光从光密介质入射到光疏介质的光滑分界面时,入射光与分界面法线的夹角大于或者等于临界角,入射的光束将全部反射光密介质中。光线在每一次全反射时,在芯皮的界面都存在着光的渗透现象,皮层玻璃的厚度满足下式时才能使得光学纤维理想传光:
$$Y \geqslant \dfrac{{\dfrac{1}{2}\lambda }}{{ \sqrt{{\dfrac{3}{{\sqrt {1 + {{\Bigg(\dfrac{{{\text{π}} R}}{H}\Bigg)}^2}} }}} - \dfrac{3}{4}} - \dfrac{1}{2}}}$$ (1) 式中:Y为扭后拉伸变形的光学纤维丝的皮层厚度;H为有效扭区宽度(即实际扭区宽度);R为光纤倒像器有效区半径;λ为入射光波长。我们将R/H定义为径宽比,径宽比越大,光学纤维的拉伸量越大;径宽比越小,光学纤维的拉伸量越小。
从(1)式可以看出,有效扭区宽度H越小,光学纤维皮层所需厚度就越厚,为了获得高分辨率的大口径光纤倒像器,我们设计扭转变形区域的宽度H在40 mm~45 mm之间,有效区直径为Φ40 mm~Φ42 mm,芯/皮面积比控制在60/40,皮层玻璃管的厚度控制在4.3 mm~4.5 mm,这样制备出的大口径光纤倒像器能够满足传像要求的各项技术指标。继续增加皮层玻璃管厚度不仅会增加玻璃管的成型难度,也会使无效光通量增加,降低了大口径光纤倒像器的透过率和对比度。
1.3 大口径光纤倒像器的扭转成型工艺设计
大口径光纤倒像器的扭转成型工艺是将热熔压成型后的光纤板经过滚圆、切割、端面磨抛、开槽等冷加工工序制备成两端为六方柱中间为圆柱形的毛坯(图4所示为制备大口径光纤倒像器的扭转成型用六方毛坯示意图),然后将毛坯放入扭转成型炉,采用刚性夹具夹持毛坯两端六方柱侧面,加热毛坯中间圆柱部位使其受热达到玻璃的软化温度,然后对刚性夹具施加一定扭矩,使毛坯扭转180°。
大口径光纤倒像器选取的毛坯直径为60 mm、高度为60 mm,其成品的扭区需要控制在40 mm~45 mm之间,这就要求所用的加热炉的加热温区小于45 mm,该区域的温度场要力求均匀,并且加热区与其相邻区域间的温度梯度要大,同时保证处于直区的纤维温度不会因为温差过大而炸裂。常规的扭转成型设备[27]采用静态挂重扭转和单一加热源的方式,使得毛坯在扭转成型过程中上下左右受热不均,且光学纤维拉伸变形使得扭区角度分布极不均匀,扭转角度过度集中于扭区的中间位置,导致无法扭转制备出合格产品。
2 实验
2.1 大口径光纤倒像器的制备过程
如图5所示为大口径光纤倒像器的制备工艺流程图[28-31]。为满足产品性能指标要求,本文采用旋转差速扭转成型工艺和双炉加热系统相结合来制备大口径光纤倒像器[32]。
1) 采用旋转差速扭转成型工艺。在扭转成型炉升温加热过程中,毛坯通过装夹在旋转轴上的扭转杆卡具的夹持沿旋转轴旋转,保证光纤板毛坯均匀受热,不会因为受热不均或者温差过大而产生应力,有效避免了大尺寸光纤倒像器因受热不均或温差过大而产生的炸裂问题。当光纤板毛坯中间受热部位达到扭转成型温度时,通过设置左右扭转杆不同的旋转速度达到差速的目的,从而实现光纤板毛坯的180°扭转成型。
2) 采用双炉加热系统。如图6所示,炉内设置2个以串联方式连接的环形加热体,加热体采用电阻丝加热,该设计通过优化扭区的温度场分布,提升扭区光学纤维的变形均匀性,进而改善光纤倒像器边缘透过率和提高边缘分辨率。
通过旋转差速扭转成型工艺和双炉加热系统设计来扭转大口径光纤倒像器毛坯,使得大口径光纤倒像器的制备难点得以突破,制备出符合技术指标要求的大口径光纤倒像器。
2.2 大口径光纤倒像器分辨率测试
将制备的大口径光纤倒像器加工成成品,通过美军标1951 USAF分辨率标准测试卡进行分辨率测试。图7所示为采用常规加热系统制备的大口径光纤倒像器分辨率测试结果。图8所示为采用旋转差速扭转成型工艺和双炉加热系统设计后制备的分辨率测试结果。
2.3 大口径光纤倒像器切片角度测试
为了进一步验证大口径光纤倒像器制备工艺优化设计前后扭转区域的分布效果,我们通过光纤倒像器扭转成型区域的角度分布变化来量化表征扭转区域光学纤维的变形均匀性。将大口径光纤倒像器成品切割成3.0 mm±0.02 mm厚度的薄片,并在投影仪上通过10倍物镜分别测量所切割薄片的扭转角度。图9所示为大口径光纤倒像器的切片示意图。表2所示为大口径光纤倒像器的切片角度测试数据。图10所示为切片角度测试分布图。
表 2 大口径光纤倒像器切片角度数据Table 2. Angle data of slices of optical fiber image inverters with large aperture参数 1# 2# 3# 4# 5# 6# 7# 8# 9# 10# 11# 12# 13# 14# 15# 备注 距离/mm 3.0 6.4 9.8 13.2 16.6 20 23.4 26.8 30.2 33.6 37.0 40.4 43.8 47.2 50.6 常规对比 扭转角/(°) 0 0.32 1.23 3.9 9.01 16.3 24.25 29.41 27.98 21.77 13.61 6.67 2.97 0.35 0 角度集中量/
(°/mm)0.87 6.54 0.73 距离/mm 3.06 6.48 9.92 13.32 16.78 20.18 23.56 26.96 30.38 33.74 37.18 40.56 43.96 47.34 50.66 双炉系统 扭转角/(°) 0.39 1.77 4.67 9.04 14.26 19.22 23.17 24.82 22.88 17.67 12.71 7.45 3.58 1.11 0.24 角度集中量/
(°/mm)1.19 5.98 1.48 3 实验结果及分析
从图7中可以看出,大口径光纤倒像器的边缘分辨率下降明显,测试结果显示,其在128 lp/mm、144 lp/mm、161 lp/mm的分辨率测试出现“模糊不清”状态。由图8的分辨率测试结果可以看出,采用双炉加热系统设计后制备的大口径光纤倒像器成品,其在128 lp/mm、144 lp/mm、161 lp/mm的分辨率测试清晰可见,对大口径光纤倒像器的边缘分辨率改善效果非常明显。从图10中的大口径光纤倒像器切片角度测试也可以看出,采用差速扭转和双炉加热系统的工艺设计后,切片角度由原来的“尖峰”分布转变为扭区中间角度向两侧平缓下滑趋势,使得扭区的温度场分布更均匀。
为方便描述,将扭转角度与扭转角度所在的扭区宽度的比值定义为角度集中量,即每毫米扭转角度(°/mm),“理想状态”下的角度集中量在180°/45 mm~180°/42 mm之间,即4.0°/mm~4.3°/mm的角度集中量,角度集中量越大,光学纤维的拉伸变形越不均匀。由测试结果看出,采用常规加热系统制备的大口径光纤倒像器成品,其扭区宽度为40.8 mm(2#~14#),中间20.40 mm的区域内(6#~11#)的角度集中量为6.54°/mm,两侧边缘区域的角度集中量分别为0.87°/mm和0.73°/mm,角度分布极不均匀,这是导致边缘分辨率下降的主要原因。采用差速扭转和双炉加热系统的工艺设计后,其扭区宽度为44.28 mm(1#~14#),中间20.42 mm的区域内(5#~10#)的角度集中量为5.98°/mm,两侧边缘区域的角度集中量分别为1.19°/mm和1.48°/mm,角度集中量也逐步向“理想状态”靠近,这是大口径光纤倒像器边缘分辨率改善的主要原因。
4 结论
本文研究了大口径光纤倒像器制备过程中遇到的一系列问题,分析了扭转成型过程中的制备难点,并对生产工艺进行了优化设计。通过旋转差速扭转和双炉加热系统设计,在满足大口径光纤倒像器性能指标参数的同时,能很好地拉宽扭转变形区域的纤维变形的均匀性,使得整个扭转区域光学纤维的变形趋于平缓而不发生锥变,使整个视场内的分辨率和透过率更加均匀,成功解决了大口径光纤倒像器的制备难题,使得大口径光纤倒像器的扭转角度过度集中问题得以缓解,分辨率下降问题得到明显改善。制备出了各项性能指标符合要求的产品并实现了批量供货,对大口径光纤倒像器的制备具有较大的指导意义。
-
表 1 不同地表状态反射率
Table 1 Reflectivity of different surface states
地表状态 裸地 沙漠 草地 森林 雪(紧、洁) 雪(湿、脏) 海面(太阳高度角>25°) 海面(太阳高度角 < 25°) ρg/% 10~25 25~46 15~25 10~20 75~95 25~75 < 10 10~70 表 2 仿真计算参数设置
探测时间 偏航角/(°) 俯仰角/(°) 旋转角/(°) 1 12:00AM 30 0 0 2 14:00PM 0 30 0 3 16:00PM 0 0 45 4 18:00PM 45 45 45 -
[1] 康丽珠, 赵劲松, 周倩, 等.远程迎头探测飞机目标的红外辐射特性研究[J].红外技术, 2017, 39(4): 365-371. http://d.old.wanfangdata.com.cn/Periodical/hwjs201704012 KANG Lizhu, ZHAO Jinsong, ZHOU Qian, et al. Research on infrared signature for remotely detection from the nose of aircrafts[J]. Infrared Technology, 2017, 39(4): 365-371. http://d.old.wanfangdata.com.cn/Periodical/hwjs201704012
[2] 夏新林, 艾青, 任德鹏.飞机蒙皮红外辐射的瞬态温度场分析[J].红外与毫米波学报, 2007, 26(3): 174-177. doi: 10.3321/j.issn:1001-9014.2007.03.004 XIA Xinlin, AI Qing, REN Depeng. Analysis on the transient temperature-fields for infrared radiation of aircraft skin[J]. Journal of Infrared and Millimeter Waves, 2007, 26(3): 174-177. doi: 10.3321/j.issn:1001-9014.2007.03.004
[3] LU J W, WANG Q. Aircraft-skin infrared radiation characteristics modeling and analysis[J]. Chinese Journal of Aeronautics, 2009, 22(5):493-497. doi: 10.1016/S1000-9361(08)60131-4
[4] 卢建, 辛玉林, 陈曾平.飞机目标红外成像建模与仿真[J].电子测量技术, 2010, 33(7): 65-68. http://d.old.wanfangdata.com.cn/Periodical/dzcljs201007018 LU Jian, XIN Yulin, CHEN Zengping. Modeling and simulation for the aeroplane target infrared imaging[J]. Electronic Measurement Technology, 2010, 33(7): 65-68. http://d.old.wanfangdata.com.cn/Periodical/dzcljs201007018
[5] 王科伟, 马超杰, 陈炜, 等.红外成像波门形心跟踪算法的误差分析[J].应用光学, 2009, 30(2): 353-356, 360. doi: 10.3969/j.issn.1002-2082.2009.02.038 WANG Kewei, MA Chaojie, CHEN Wei, et al. Error analysis for gate centroid tracking algorithm of infrared imaging[J]. Journal of Applied Optics, 2009, 30(2): 353-356, 360. doi: 10.3969/j.issn.1002-2082.2009.02.038
[6] 成刚, 方帆, 宁飞, 等.无人机载多光谱侦察效能研究[J].应用光学, 2017, 38(5): 685-688. doi: 10.5768/JAO201738.0501001 CHENG Gang, FANG Fan, NING Fei, et al. Study on performance of multi spectral reconnaissance of unmanned aerial vehicle[J]. Journal of Applied Optics, 2017, 38(5): 685-688. doi: 10.5768/JAO201738.0501001
[7] 韩军, 薛小乐, 王星.光电成像系统超分辨成像技术方法研究[J].应用光学, 2011, 32(1): 54-58. doi: 10.3969/j.issn.1002-2082.2011.01.013 HAN Jun, XUE Xiaole, WANG Xing. Super-resolution imaging technology for optical imaging system[J]. Journal of Applied Optics, 2011, 32(1): 54-58. doi: 10.3969/j.issn.1002-2082.2011.01.013
[8] DOUCHIN N, LATGER J. Extension of the SE-Workbench for the computation of aircraft infrared signatures[J]. SPIE, 2010, 21(6):78340I-1-13. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=CC0210396499
[9] HAYNES A W, GILMORE M A, STROUD C A. Accurate scene modeling using synthetic imagery[J]. Proceedings of SPIE-The International Society for Optical Engineering, 2003, 5075:85-96. http://d.old.wanfangdata.com.cn/NSTLHY/NSTL_HYCC026741931/
[10] GRETZMACHER F M, NYBERG S. Camouflage assessment considering human perception data[J]. Proceedings of SPIE-The International Society for Optical Engineering, 1998, 3375: 58-67. http://d.old.wanfangdata.com.cn/NSTLHY/NSTL_HYCC026118215/
[11] BASTIAANSEN M C, BCKER K B, CLUITMANS P J, et al. CAMEVA: a methodology for computerized evaluation of camouflage effectiveness and estimation of target detectability[J]. Proceedings of SPIE-The International Society for Optical Engineering, 1999, 40(9): 229-238. http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ021753329/
[12] JACOBS P A M. Experimental evaluation of camouflage effectiveness in the thermal infrared[C].[S.l.]: SPIE, 1995.
[13] 刘娟, 龚光红, 韩亮, 等.飞机红外辐射特性建模与仿真[J].红外与激光工程, 2011, 40(7) :1209-1213. doi: 10.3969/j.issn.1007-2276.2011.07.004 LIU Juan, GONG Guanghong, HAN Liang, et al. Modeling and simulation of aircraft infrared radiation characteristics[J]. Infrared and Laser Engineering, 2011, 40(7):1209-1213. doi: 10.3969/j.issn.1007-2276.2011.07.004
[14] 杨风暴.红外物理与技术[M].北京:电子工业出版社, 2014. YANG Fengbao. Infrared physics and technology[M]. Beijing:Publishing House of Electronics Industry, 2014.
[15] 张可.飞行器蒙皮红外辐射图像生成技术研究[D].南京: 南京航空航天大学, 2010. http://www.cnki.com.cn/Article/CJFDTotal-JGHW201103008.htm ZHANG Ke. Research on technology of creating aircraft's skin infrared images[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2010. http://www.cnki.com.cn/Article/CJFDTotal-JGHW201103008.htm
[16] 杨立.红外热像仪测温计算与误差分析[J].红外技术, 1999, 21(4): 20-24. http://www.cnki.com.cn/Article/CJFDTotal-HWJS199904005.htm YANG Li. Calculation and error analysis of temperature measurement using thermal image[J]. Infrared Technology, 1999, 21(4):20-24. http://www.cnki.com.cn/Article/CJFDTotal-HWJS199904005.htm