留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

飞机红外辐射空间分布特性仿真研究

童忠诚 王亮 吴俊

童忠诚, 王亮, 吴俊. 飞机红外辐射空间分布特性仿真研究[J]. 红外与激光工程, 2023, 52(1): 20220264. doi: 10.3788/IRLA20220264
引用本文: 童忠诚, 王亮, 吴俊. 飞机红外辐射空间分布特性仿真研究[J]. 红外与激光工程, 2023, 52(1): 20220264. doi: 10.3788/IRLA20220264
Tong Zhongcheng, Wang Liang, Wu Jun. Simulation research on IR radiation space distribution characteristic of fight plane[J]. Infrared and Laser Engineering, 2023, 52(1): 20220264. doi: 10.3788/IRLA20220264
Citation: Tong Zhongcheng, Wang Liang, Wu Jun. Simulation research on IR radiation space distribution characteristic of fight plane[J]. Infrared and Laser Engineering, 2023, 52(1): 20220264. doi: 10.3788/IRLA20220264

飞机红外辐射空间分布特性仿真研究

doi: 10.3788/IRLA20220264
详细信息
    作者简介:

    童忠诚,男,教授,硕士生导师,博士,主要从事光电技术方面的研究

  • 中图分类号: TN219

Simulation research on IR radiation space distribution characteristic of fight plane

  • 摘要: 为深入了解飞机在3~5 μm波段红外辐射的空间分布特性,在建立飞机尾焰红外辐射的椭圆体模型后,利用飞机蒙皮红外辐射模型和尾喷口红外辐射模型,仿真计算了飞机在不同方向上3~5 μm波段的红外辐射,给出了红外辐射的空间分布曲线。计算发现飞机的红外辐射关于机翼平面和机体纵向对称平面对称,在机尾和机头各有4个极值。其中在机尾辐射方向为($ \pm $150°,$ \pm $32°)时红外辐射的极大值达到5 177 W,在机头辐射方向为($ \pm $68°,$ \pm $64°)时红外辐射的极大值达到3 461 W。通过分析飞机红外辐射在机翼平面、机体纵向对称平面和极值平面的分布规律,发现在机头正向红外辐射非常小,但当辐射方向不在机轴正向时,随着辐射方向与机轴方向夹角的增大红外辐射急剧增加。还发现在机翼平面辐射方向投影与机轴夹角越小,在机体纵向对称面上红外辐射出现极值的方向越靠近机轴。同样在机体纵向对称面上辐射方向投影与机轴夹角越小,则在机翼平面上红外辐射出现极值的方向也越靠近机轴。
  • 图  1  机体坐标系

    Figure  1.  Coordinate system of fight plane

    图  2  尾焰红外辐射模型

    Figure  2.  IR radiation model of tail flare

    图  3  尾焰混合区红外辐射模型

    Figure  3.  IR radiation model of tail flare mixed region

    图  4  飞机3~5 μm波段红外辐射的空间分布

    Figure  4.  IR radiation space distribution of fight plane in 3-5 μm band

    图  5  飞机3~5 μm波段机翼平面上的红外辐射分布

    Figure  5.  IR radiation distribution of fight plane in wing plane in 3-5 μm band

    图  6  飞机3~5 μm波段红外辐射纵对称平面上分布

    Figure  6.  IR radiation distribution of fight plane in longitudinal symmetrical plane in 3-5 μm band

  • [1] Chai S J, Tong Z X, Li J X, et al. Infrared radiation characteristics and detection of stealth aeroplane [J]. Fire Control & Command Control, 2014, 39(8): 26-33. (in Chinese)
    [2] Sun W, Wang B. Modeling and simulation of infrared radiation from aerial combat aircraft by fluent[C]//2019 Proceedings of the Conference on Control and its Applications (CT), SIAM, 2019: 25-29.
    [3] Sidonie Lefebvre, Antoine Roblin, Suzanne Varet, et al. A methodological approach for statistical evaluation of aircraft infrared signature [J]. Reliability Engineering & System Safety, 2010, 95(5): 484-493.
    [4] Xu D G, Feng W L, Sang J H. IR characteristics of aircraft aft fuselage skin [J]. Infrared and Laser Engineering, 2013, 42(1): 7-13. (in Chinese)
    [5] Wang X, Gao S L, Jin L, et al. Muiti-band infrared radiation characterization and simulation analysis for aerial target [J]. Acta Photonica Sinica, 2020, 49(5): 0511002. (in Chinese) doi:  10.3788/gzxb20204905.0511002
    [6] Liu Y, Jin Z Q. Simulation of air vehicle motion control models and infrared images [J]. Infrared Technology, 2020, 42(9): 863-872. (in Chinese) doi:  10.3724/SP.J.7102910262
    [7] He P, Wang Y Y, Yue S H. Operating range of the advanced infrared detector for the stealth aircraft [J]. Infrared Technology, 2020, 42(9): 899-904. (in Chinese) doi:  10.3724/SP.J.7102910267
    [8] Kang L Z, Zhao J S, Zhou Q, et al. Research on infrared signature for remotely detection from the nose of aircrafts [J]. Infrared Technology, 2017, 39(4): 365-370. (in Chinese)
    [9] Zhang Y J, Xu Z Y, Li J X. Study on simulation model of aircraft infrared hyperspectral image [J]. Aero Weaponery, 2020, 27(4): 91-96. (in Chinese)
    [10] Tang S W, Liu T Y. Radiation computation of tailpipe based on grid [J]. Electro-Optic Technology Application, 2020, 42(9): 863-872. (in Chinese)
    [11] Wang C Z, Tong Z X, Lu Y L. Study on the airplane's infrared radiation characteristics [J]. Laser & Infrared, 2011, 41(9): 996-1001. (in Chinese)
    [12] Yin X M, Wang L, Zhang X L. Wide band k-distribution model for calculating radiation signal of liquid rocket exhaust plume [J]. Laser & Infrared, 2020, 50(11): 1352-1356. (in Chinese)
    [13] Wang X, Chen H C. Algorithm of infrared signature of aeroplane along the line of sight [J]. Journal of Atmospheric and Environmental Optics, 2008, 3(3): 217-222. (in Chinese)
    [14] Wu X Z, Teng P, Lu Y, et al. Calculation and simulation of infrared radiation about some jet aircraft [J]. Infrared Technology, 2008, 30(12): 727-731. (in Chinese)
    [15] Wang W J, Huang J, Yuan G F, et al. Modeling and application analysis of operating range of air-based infrared system [J]. Optics and Precision Engineering, 2020, 28(6): 1295-1301. (in Chinese) doi:  10.3788/OPE.20202806.1295
    [16] Li Y Z, Sun B, Wang D H. Building model of aeroplane infrared radiant and simulation [J]. Infrared Technology, 2008, 30(5): 252-255. (in Chinese)
    [17] Liu Z Y, Ding F, Ye Q. A practical calculation method of the radiation parameters of particles in rocket plume [J]. Laser & Infrared, 2020, 50(3): 315-318. (in Chinese)
    [18] Zhang H X, Zhang J Q, Yang W, et al. Theoretical calculation of the IR radiation of an aeroplane [J]. Journal of Xidian University, 1997, 24(1): 78-81. (in Chinese)
    [19] Zhu B, Pei Y, Song B F. Estimation method of aircraft infrared radiant intensity based on aerodynamic computation [J]. Fire Control & Command Control, 2011, 36(8): 180-183. (in Chinese)
    [20] Cao X Z, Guo L H. Measurement of the infrared radiation characteristic of the aerobat based on images [J]. Optical Technique, 2007, 33(2): 299-301. (in Chinese)
  • [1] 吕蓉, 牛青林, 王晓冰.  攻角对临近空间飞行器侧喷射流红外辐射特性影响数值模拟 . 红外与激光工程, 2024, 53(1): 20230176-1-20230176-11. doi: 10.3788/IRLA20230176
    [2] 程稳, 吉洪湖, 施小娟, 王浩.  带全遮挡导流支板的二元排气系统与后体组合模型的流动传热与红外特性分析 . 红外与激光工程, 2022, 51(9): 20210965-1-20210965-11. doi: 10.3788/IRLA20210965
    [3] 程稳, 孙啸林, 马姗.  基于假设气体法的燃气辐射特性计算模型 . 红外与激光工程, 2022, 51(7): 20220286-1-20220286-9. doi: 10.3788/IRLA20220286
    [4] 于鲲, 郭彪, 丛明煜.  空间目标临边探测背景红外成像建模与图像仿真 . 红外与激光工程, 2019, 48(9): 904005-0904005(10). doi: 10.3788/IRLA201948.0904005
    [5] 戴聪明, 赵凤美, 刘栋, 唐超礼, 魏合理.  强吸收带卷云大气红外辐射特性研究 . 红外与激光工程, 2018, 47(12): 1204003-1204003(7). doi: 10.3788/IRLA201847.1204003
    [6] 王丁, 吉洪湖.  波瓣混合器排气系统红外特性的模型实验和数值分析 . 红外与激光工程, 2017, 46(1): 104003-0104003(7). doi: 10.3788/IRLA201746.0104003
    [7] 常帅, 高秀伟, 佟首峰, 宋延嵩, 董岩.  红外辐射式干扰弹能量测量装置研究 . 红外与激光工程, 2016, 45(1): 117002-0117002(5). doi: 10.3788/IRLA201645.0117002
    [8] 李建勋, 童中翔, 樊晓光, 吴利荣, 陈超, 王超哲, 贾林通.  红外对抗过程的效能评估方法研究 . 红外与激光工程, 2016, 45(3): 304008-0304008(10). doi: 10.3788/IRLA201645.0304008
    [9] 周金伟, 石志广, 李吉成.  NO 有效带宽的k 分布算法 . 红外与激光工程, 2015, 44(2): 432-437.
    [10] 孙占久, 聂宏, 黄伟.  无人机红外辐射特性计算与分析 . 红外与激光工程, 2014, 43(4): 1037-1046.
    [11] 黄晨, 梁小虎, 王建军, 高昕.  临近空间红外探测像面照度建模与仿真 . 红外与激光工程, 2014, 43(10): 3217-3221.
    [12] 任利锋, 张靖周, 单勇, 刘喜岳.  旋翼气动加热对直升机红外特征分布的影响 . 红外与激光工程, 2014, 43(1): 53-60.
    [13] 徐顶国, 冯维林, 桑建华.  飞行器后机身蒙皮红外辐射特性 . 红外与激光工程, 2013, 42(1): 7-13.
    [14] 徐顶国, 桑建华, 罗明东.  无人机蒙皮红外辐射特征研究 . 红外与激光工程, 2013, 42(4): 880-884.
    [15] 李建勋, 童中翔, 刘万俊, 王超哲, 张志波, 禚真福.  航空发动机红外辐射实验与仿真 . 红外与激光工程, 2013, 42(3): 549-555.
    [16] 成志铎, 李明博, 李健, 常晓权, 刘君.  目标与背景的红外辐射特性仿真方法 . 红外与激光工程, 2013, 42(9): 2336-2340.
    [17] 路远, 冯云松, 凌永顺, 乔亚.  飞行器尾焰红外辐射及其被动测距 . 红外与激光工程, 2013, 42(7): 1660-1664.
    [18] 许爱华, 汪中贤, 于坚, 张正武, 卓家靖, 史建军, 赵京山.  高空高速无人机尾焰红外辐射特性计算研究 . 红外与激光工程, 2012, 41(7): 1700-1776.
    [19] 王明明, 郝颖明, 朱枫, 付双飞, 石坤.  空中目标红外辐射特性计算与实时仿真 . 红外与激光工程, 2012, 41(8): 1979-1984.
    [20] 王 盈, 来霄毅, 黄建明, 魏祥泉.  基于Sinda/G和Matlab的空间目标红外辐射特性分析 . 红外与激光工程, 2012, 41(5): 1113-1118.
  • 加载中
图(6)
计量
  • 文章访问数:  343
  • HTML全文浏览量:  58
  • PDF下载量:  128
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-04-19
  • 修回日期:  2022-06-17
  • 刊出日期:  2023-01-18

飞机红外辐射空间分布特性仿真研究

doi: 10.3788/IRLA20220264
    作者简介:

    童忠诚,男,教授,硕士生导师,博士,主要从事光电技术方面的研究

  • 中图分类号: TN219

摘要: 为深入了解飞机在3~5 μm波段红外辐射的空间分布特性,在建立飞机尾焰红外辐射的椭圆体模型后,利用飞机蒙皮红外辐射模型和尾喷口红外辐射模型,仿真计算了飞机在不同方向上3~5 μm波段的红外辐射,给出了红外辐射的空间分布曲线。计算发现飞机的红外辐射关于机翼平面和机体纵向对称平面对称,在机尾和机头各有4个极值。其中在机尾辐射方向为($ \pm $150°,$ \pm $32°)时红外辐射的极大值达到5 177 W,在机头辐射方向为($ \pm $68°,$ \pm $64°)时红外辐射的极大值达到3 461 W。通过分析飞机红外辐射在机翼平面、机体纵向对称平面和极值平面的分布规律,发现在机头正向红外辐射非常小,但当辐射方向不在机轴正向时,随着辐射方向与机轴方向夹角的增大红外辐射急剧增加。还发现在机翼平面辐射方向投影与机轴夹角越小,在机体纵向对称面上红外辐射出现极值的方向越靠近机轴。同样在机体纵向对称面上辐射方向投影与机轴夹角越小,则在机翼平面上红外辐射出现极值的方向也越靠近机轴。

English Abstract

    • 运用红外探测系统实现对飞机的探测与跟踪成为对空预警探测的发展趋势。目前,国内外对飞机红外辐射的研究较多[1-4],但研究主要集中在飞机在平面内的红外辐射分布,对其空间分布特性鲜有文献报道。天空背景在3~5 um波段红外辐射比较弱,利用飞机3~5 μm波段的红外辐射实现对来袭飞机的逼近告警具有一定的优势,因此,研究飞机在3~5 μm波段的红外辐射空间分布特征,对红外搜索跟踪系统对飞机的识别,以及对飞机红外隐身技术的研究,都具有参考意义。

    • 建立如图1所示的机体坐标系计算飞机任意方向上的红外辐射[5]。平面1为飞机机翼平面,平面2为飞机纵向对称平面,平面1与平面2的交线为机轴,也是坐标系的X轴,机头方向为正。辐射方向与机轴的交点O为坐标原点,Y轴在平面2内且垂直于X轴。辐射方向OI在平面1投影OI1X轴夹角为$ \theta $OI在平面2左侧为正,在平面2右侧为负。辐射方向OI在平面2投影OI2X轴夹角为$ \gamma $OI在平面1上方$ \gamma $为正,在平面1下方$ \gamma $为负。

      图  1  机体坐标系

      Figure 1.  Coordinate system of fight plane

      如果能够计算出在飞机机翼平面OI1方向上的红外辐射强度I1和在飞机纵向对称平面OI2方向上的红外辐射强度I2,则OI方向的红外辐射强度I为:

      $$ I = {I_1}\cos \alpha + {I_2}\cos \beta $$ (1)

      式中:$ \alpha $OIOI1的夹角;$\; \beta $OIOI2的夹角。$ \alpha $$\; \beta $由下式计算:

      $$\begin{split} & \alpha = 2\arcsin \left( { \frac{{\sqrt 2 }}{2}\cos \gamma \sqrt {1 - \cos \theta } } \right) \\& \beta = 2\arcsin \left( { \frac{{\sqrt 2 }}{2}\sqrt {1 - \cos \gamma } } \right) \end{split} $$ (2)

      设计算域的中心位于机翼平面OI1和飞机纵向对称平面OI2的交线上,与机头的距离为飞机长度的一半。计算域的长度为机长和尾焰长度2倍的和,宽度为机翼的宽度,高度为垂直尾翼高度的2倍[6]

    • 飞机的红外辐射包括尾焰红外辐射、尾喷管口红外辐射和气动加热的蒙皮红外辐射[7]

    • 尾焰可分为稳定区和混合区,如图2所示[8]。稳定区内温度恒定。稳定区向外直至尾焰边界为混合区,温度逐渐下降。设尾喷口直径为${D_{\rm P}}$,混合区尾焰最宽处直径为${D_{\rm W}}$,稳定区的长为$ {L_1} $,尾喷口到混合区最宽处的长度为$ {L_2} $。将稳定区看为长轴为$ {L_1} $,短轴为${D_{\rm P}}/2$的半个长轴椭圆体,混合区看为长轴为$ {L_2} + {L_0} $,短轴为${D_{\rm W}}/2$的长轴椭圆体。

      图  2  尾焰红外辐射模型

      Figure 2.  IR radiation model of tail flare

      图2$ {L_0} $为:

      $$ {L_0} = \left( { \sqrt {\frac{{{D_{\rm W}^2}}}{{{D_{\rm W}}^2 - {D_{\rm P}}^2}}} - 1 } \right) \cdot {{{L}}_{\text{2}}} $$ (3)
    • (1)稳定区温度

      稳定区温度[9-10]为尾焰在尾喷口处的温度。尾焰在尾喷口处的温度与尾喷口的温度有关,其计算公式为:

      $$ {T}_{2}={T}_{1}\left( { \frac{{P}_{2}}{{P}_{1}} } \right)^{\tfrac{\gamma _a-1}{\gamma_a}} $$ (4)

      式中:$ {T_2} $为尾焰在尾喷口处温度;$ {T_1} $为尾喷口温度;$ {P_2} $为膨胀后的气体压力;$ {P_1} $为尾喷口内的气体压力;$ \gamma_a $为气体的定压热容量和定容热容量之比。

      (2)混合区温度[11-12]

      将混合区划分成若干个离心率相同的长轴椭圆体,如图3中的a1a2···an。设a1a2···an表面积分别为S1S2···Sn,且认为椭圆体表面上温度相同,分别为T1T2···Tna1a2···an都不是完整的长轴椭圆体,需补充完整,如图中点状虚线所示。对于a1a2···an,补全部分不完全一样,但其表面积相差不多。

      图  3  尾焰混合区红外辐射模型

      Figure 3.  IR radiation model of tail flare mixed region

      设第n个长轴椭圆体表面距喷口的最远距离为$ {l_n} $,则其短轴${{{D}}_n}$和表面积$ {S_n} $为:

      $$ \begin{split} & {{{D}}_n} = \frac{{{l_n} + L{}_0}}{{2(L{}_2 + L{}_0)}} \cdot {D_{\rm W}} \\ & {S_n} = \frac{ {{\pi}} }{2} \cdot {D_n}^2 + \frac{ {{\pi}} }{2} \cdot {D_n} \cdot ({L_0} + {l_n})\arcsin ({e_0})/{e_0} - S' \\ &{e}_{0}=\frac{\sqrt{({L}_{2}+{L}_{0}{)}^{2}-({D}_{\rm W}/2{)}^{2}}}{{L}_{2}+{L}_{0}} \end{split}$$ (5)

      $ S' $为补全部分的表面积,其计算公式为:

      $$ \begin{split} & S' = \frac{ {{\pi}} }{4} \cdot {D_{\rm P}}^2 + \frac{ {{\pi}} }{2} \cdot {D_{\rm P}} \cdot {L_0} \cdot \arcsin ({e_1})/{e_1} \\ &{e}_{1}=\frac{\sqrt{{L}_{0}{}^{2}-({D}_{\rm P}/2{)}^{2}}}{{L}_{0}} \end{split} $$ (6)

      稳定区的表面积$ {A_0} $为:

      $$\begin{split} & {{{A}}_0} = \frac{ {{\pi}} }{4} \cdot {D_{\rm P}}^2 + \frac{ {{\pi}} }{2} \cdot {D_{\rm P}} \cdot {L_1} \cdot \arcsin ({e_2})/{e_2} \\ & {e}_{2}=\frac{\sqrt{{L}_{1}{}^{2}-({D}_{\rm P}/2{)}^{2}}}{{L}_{1}} \end{split} $$ (7)

      a1与稳定区的温差为$ \Delta T $anan−1的温差为$ \Delta {T_n} $,则有:

      $$ \Delta {T_n} = {S_0} \cdot \Delta T/{S_n} {T_n} = {T_0} - \left(\Delta T + \sum\limits_1^{n - 1} {\Delta {T_i}} \right) $$ (8)

      式中:$ {T_n} $为第n个长轴椭圆体的表面温度。

    • 稳定区在平面2上沿OI1方向的投影面积为[13]

      $$ {A_0} = \frac{\pi }{2} \cdot {D_{\rm P}} \cdot {L_1} \cdot \cos \theta $$ (9)

      n个长轴椭圆体在平面2上沿OI1方向的投影面积An为:

      $$ \left\{ {\begin{array}{*{20}{l}} {{A_1} = \left( { \dfrac{ {{\pi}} }{4} \cdot {D_1} \cdot \left( { {L_0} + {l_1} } \right) - \dfrac{ {{\pi}} }{4} \cdot {D_{\rm P}} \cdot \left( { {L_0} + {L_1} } \right) } \right)\cos \theta }& {{{n}} = 1} \\ {{A_n} = \left( { \dfrac{ {{\pi}} }{4} \cdot {D_n} \cdot \left( { {L_0} + {l_n} } \right) - \dfrac{ {{\pi}} }{4} \cdot {D_{n - 1}} \cdot \left( { {L_0} + {l_{n - 1}} } \right) } \right)\cos \theta }& {n \gt 1} \end{array}} \right. $$ (10)

      计算它们在平面1上沿OI2方向的投影面积时用$ \omega $替代$ \theta $,后面类似情形同样处理。

    • 尾喷口可看成灰体,其红外辐射与尾喷口温度${{T}}_{\rm P}$和尾喷口直径有关。尾喷口在平面2上沿OI1方向的投影面积${A_{\rm P}}$[14]

      $$\begin{split} & {A_{\rm P}} = \frac{ {{\pi}} }{4}{D_{\rm P}}^2\cos \theta \\& \frac{ {{\pi}} }{2} \lt \theta \lt {{\pi}} 或- {{\pi}} \lt \theta \lt -\frac{ {{\pi}} }{2} \end{split} $$ (11)
    • 在大气层内蒙皮的温度可用下式计算[15-16]

      $$ {T_{\rm s}} = {T_{\rm a}}(1 + 0.164{V^2}) $$ (12)

      式中:${T_{\rm s}}$为蒙皮温度;${T_{\rm a}}$为周围大气的温度;$ V $为飞机的速度。

    • 设机体上方和下方蒙皮面积相等为${A_{\rm u0}}$,机头方向蒙皮面积为${A_{\rm h0}}$,机体侧面蒙皮面积为${A_{\rm s0}}$,则机体上方和下方蒙皮在平面2上的投影面积和机体侧面在平面1上的投影面积为0。

      其中机头蒙皮在平面2上沿OI1方向的投影面积为:

      $$ {A_{\rm h}} = {A_{\rm h0}}\cos \theta \qquad - \frac{ {{\pi}} }{2} \lt \theta \lt \frac{ {{\pi}} }{2} $$ (13)

      机体侧面蒙皮在平面2上沿OI1方向的投影面积为:

      $$ {A_{\rm s}} = {A_{\rm s0}}\sin \theta $$ (14)

      机体上下蒙皮在平面1上沿OI2方向的投影面积为:

      $$ {A_{\rm u}} = {A_{\rm u0}}\sin \omega $$ (15)
    • 尾焰、尾喷口和蒙皮都可以认为是具有一定发射率的灰体,知道它们的发射率、辐射温度和有效辐射面积,就可以得到飞机的辐射强度为:

      $$ I(\lambda ,T) = \sum {{\varepsilon _n} \cdot {{t{A}}_n} \cdot M(\lambda ,{T_n})} + {\varepsilon _{\rm P}} \cdot {A_{\rm P}} \cdot M(\lambda ,{T_1}) \cdot {{\rm e}^{ - \mu \cdot r}} $$ (16)

      式中:$ {\varepsilon _n} $${{{A}}_n}$$ {T_n} $分别为尾焰或蒙皮相应计算单元的红外发射率、有效辐射面积和温度;$ M(\lambda ,{T_n}) $为普朗克定律;${\varepsilon _{\rm P}}$为尾喷口的红外发射率;$ \mu $为尾焰对尾喷口红外辐射的衰减系数;$ r $为在辐射方向上尾喷口红外辐射在尾焰中的传输距离[16]

      $$ r = ({L_2} + {L_0}) \cdot (1 - {e}{{}_0^2})/(1 - {e}{}_0\cos \theta ) $$ (17)

      公式(16)在计算尾喷口红外辐射时考虑了尾焰对它的衰减[17-18]

    • 根据上述模型计算了飞机3~5 μm波段、不同辐射方向上的红外辐射强度,如图4所示。仿真时所用参数主要根据某重型战斗机巡航时的红外图像选取,部分参数根据文献选取。具体为:${D_{\rm P}}$取0.8 m、${D_{\rm W}}$取1.5 m、$ {L_1} $取1.2 m、$ {L_2} $取4 m、${{{T}}_1}$取800 K、$ \gamma $取1.3、P2/P1取0.5、$ V $取300 m/s、${T_{\rm a}}$取250 K、${A_{\rm u0}}$取106 m2${A_{\rm s0}}$取33 m2${A_{\rm h0}}$取8 m2$ \mu $取0.15、蒙皮的发射率取0.5、尾喷口的发射率取0.9、尾焰的发射率取0.85[18-19]

      图  4  飞机3~5 μm波段红外辐射的空间分布

      Figure 4.  IR radiation space distribution of fight plane in 3-5 μm band

    • 图4(a)为飞机3~5 μm波段红外辐射空间分布的三维图,图(b)~(e)分别为三维图的俯视图、侧视图、前视图和后视图。

      图4(a)可知,飞机的红外辐射关于机翼平面和飞机纵向对称平面对称。在机尾和机头各有4个强度相同的极值。根据仿真数据,分析得到机尾极值为5177 W,辐射方向为($ \pm $150o$ \pm $32°)。前头极值为3461 W,辐射方向为($ \pm $68°,$ \pm $64°)。

      图4(b)和(c)可知,飞机的红外辐射在水平面和纵平面的分布都呈鸭梨状,但在纵平面内分布更扁一些,原因是飞机上下表面蒙皮比侧表面蒙皮面积更大。

      图4(d)和(e)可知,在机尾正向和机头正向方向上,飞机的红外辐射比其邻近方向要小,特别是在机头正向上飞机的红外辐射极小。

      从结果来看,仿真得出的飞机红外辐射空间分布曲线与参考文献[1]、[14]一致,均为鸭梨状。但参考文献[14]得出机体红外辐射极值为4500 W,略小于文中极值5177 W,其原因是文中研究机型为某重型飞机,红外辐射的强度与其他机型略有差异。依据参考文献[20]中提供的等效辐射方程,利用某重型战斗机巡航时的红外图像,反推出该战斗机在(145°,−17°)方向上的辐射为4230 W,仿真得到在该方向上的红外辐射为4 823 W,两者相差约为14%,其主要原因是建模时没有考虑不同部位蒙皮之间的遮蔽效应,特别是尾翼蒙皮对机体蒙皮的遮蔽。

    • 根据仿真数据给出了$ \gamma $=0°全方位、$ \gamma $=32°时机尾方向和$ \gamma $=64°时机头方向红外辐射水平面的分布曲线,如图5(a)~(c)所示。其中给出$ \gamma $=32°时机尾方向和$ \gamma $=64°时机头方向红外辐射水平面分布曲线的原因是根据图4在这两个方向上存在极值。

      图5(a)可知,在$ \gamma $=0°时飞机红外辐射水平面分布曲线呈鸭梨状,其极大值为4543 W,此时$ \theta $$ \pm $138°。当$ \theta = {0} $°时,即在机头正面时红外辐射极低,但当$ \theta \ne {0} $°时红外辐射随着$ \theta $的增大急剧增加,其原因是当辐射方向偏离机轴时,飞机侧面蒙皮和尾焰有效辐射面积急剧增加。当$ \theta = {180} $°时,在机尾方面红外辐射有最小值,为3493 W;且当$ \theta $偏离180°时,红外辐射会增加,但增加趋势不如机头方向。

      图5(b)可知,在$ \gamma $=32°时,机尾红外辐射水平分布的极大值出现在$ \theta = \pm {150} $°的方向上。可见随着$ \gamma $的变大,机尾红外辐射极大值的方向越靠近机轴。同样由图5(c)可知,在$ \gamma $=64°时,机头红外辐射水平分布的极大值出现在$ \theta = \pm {68} $°的方向上。而由5(a)可知,当$ \gamma $=0°时,机头红外辐射极值方向为$ \theta = \pm {90} $°,即当$ \gamma $变大时,机头红外辐射极大值出现的方向也越靠近机轴。

      图  5  飞机3~5 μm波段机翼平面上的红外辐射分布

      Figure 5.  IR radiation distribution of fight plane in wing plane in 3-5 μm band

    • 飞机红外辐射的纵平面分布如图6所示。图6(a)中曲线以$ \gamma $$ \pm $90°为分界线,右侧曲线对应的$ \theta $=0°,左侧曲线对应的$ \theta $=180°。图6(b)和图6(c)分别为$\theta = {150^ \circ}$$\theta = {68^\circ}$时飞机红外辐射的纵平面分布曲线。

      图6(a)可知飞机红外辐射纵平面分布曲线与其在水平面十分相似,其中$ \gamma $$ \pm $40°时,机尾方向极大值为4528 W。

      图6(b)和图6(c)可知,当$ \theta = {150} $°时,机尾红外辐射纵平面分布的极大值出现在$ \gamma = \pm {32} $°的方向上,当$ \theta = {68} $°时机头红外辐射纵平面分布的极大值出现在$ \gamma = \pm {64} $°的方向上,都呈现出在水平面上$ \theta $的方向越靠近机轴,则在纵平面上红外辐射出现极值的方向也越靠近机轴。

      图  6  飞机3~5 μm波段红外辐射纵对称平面上分布

      Figure 6.  IR radiation distribution of fight plane in longitudinal symmetrical plane in 3-5 μm band

    • 文中建立了一种尾焰红外辐射的椭圆体模型,计算得到了飞机3~5 μm波段红外辐射的空间分布曲线。计算结果表明:飞机红外辐射的空间分布关于机翼平面和飞机纵向对称面对称,且在机头和机尾方向各有4个极值。在机头正向红外辐射极低,但辐射方向一旦偏离机轴,红外辐射急剧增加。研究还发现:辐射方向在机翼平面(或机体纵向对称面)上的投影与机轴夹角越小,则在机体纵向对称面(或机翼平面)上红外辐射出现极值的方向越靠近机轴。

参考文献 (20)

目录

    /

    返回文章
    返回