DOI :10. 13700/j . bh . 1001-5965. 2009. 10. 0152009年10月北京航空航天大学学报
第35卷第10期J o u r n a l o f B e i j i n g U n i v e r s i t y o f A e r o n a u t i c s a n d A s t r o n a u t i c s
O c t o b e r 2009
V o l . 35 No . 10
飞机红外辐射及大气透过率计算方法
毛 峡
胡海勇
黄 康
梁晓庚
(北京航空航天大学电子信息工程学院, 北京100191) (中国空空导弹研究院, 洛阳417009)
摘 要:研究了飞机目标的红外辐射特性, 以及红外辐射在大气中的辐射衰减特性. 根据现有的理论基础以及飞机自身结构特性, 将飞机的红外辐射源分为蒙皮、尾喷口和羽流3部分, 并提出了3个辐射源在不同波段、不同角度下的红外辐射特性计算方法. 根据红外辐射
的大气传输特性, 采用了不同波长、不同弹目距离下的大气透过率简易计算方法, 并将其与L o w t r a n 7进行比较. 计算了飞机在不同波段下, 经大气衰减后, 最终到达红外探测器的辐射强度.
关 键 词:红外辐射; 蒙皮; 尾喷口; 羽流; 大气透过率; 飞机
中图分类号:T N 219
文献标识码:A 文章编号:1001-5965(2009) 10-1228-04
C a l c u l a t i o nm e t h o d f o r a i r p l a n eI Rr a d i a t i o na n d a t m o s p h e r i c t r a n s m i t t a n c e
M a o X i a Hu H a i y o n g Hu a n g K a n g
(S c h o o l o f E l e c t r o n i c s a n dI n f o r m a t i o n E n g i n e e r i n g , B e i j i n gU n i v e r s i t y o f A e r o n a u t i c s a n dA s t r o n a u t i c s , B e i j i n g 100191, C h i n a )
L i a n g X i a o g e n g
(C h i n a A i r b o r n e M i s s i l eA c a d e m i c , L u o y a n g 417009, C h i n a )
A b s t r a c t :T h e i n f r a r e d r a d i a t i o n o f a i r p l a n e t a r g e t s a n d t h e r a d i a t i o n a t t e n u a t i o n o f t h e i n f r a r e d r a d i a t i o n
i n t h e a t m o s p h e r e w e r e r e s e a r c h e d . B a s e d o n t h e e x i s t i n g t h e o r e t i c a l s t u d y a n d t h e s t r u c t u r e o f t h e a i r p l a n e , t h e i n f r a r e d r a d i a t i o n s o u r c e s o f t h e a i r p l a n e w e r e d i v i d e d i n t o t h r e e c a t e g o r i e s :f u s e l a g e , h o t s p o t a n d p l u m e . F u r t h e r m o r e , c a l c u l a t i o n a p p r o a c h e s t o t h e c h a r a c t e r o f i n f r a r e d r a d i a t i o n f o r t h e t h r e e r a d i a t i o n s o u r c e s a t d i f -f e r e n t w a v e b a n d a n d d i f f e r e n t a n g l e w e r e p r o p o s e d . A c c o r d i n g t o i n f r a r e d r a d i a t i o n t r a n s m i t t a n c e i n t h e a t m o s -p h e r e , a c a l c u l a t i o n m e t h o d f o r t h e t r a n s m i t t a n c e o f a t m o s p h e r e w a s p u t f o r w a r d a t d i f f e r e n t w a v e l e n g t h a n d d i f f e r e n t t a r g e t d i s t a n c e , a n d t h e n t h i s m e t h o d w a s c o m p a r e d w i t h L o w t r a n 7. F i n a l l y , t h e r a d i a n t i n t e n s i t y r e -c e i v e d b y t h e I Rd e t e c t o r w a s c o m p u t e d a f t e r a t m o s p h e r i c a t t e n u a t i o n a t d i f f e r e n t w a v e b a n d .
K e y w o r d s :i n f r a r e d r a d i a t i o n ; f u s e l a g e ; h o t s p o t ; p l u m e ; a t m o s p h e r i c t r a n s m i t t a n c e ; a i r p l a n e
外场实验获取飞机目标红外图像的费用昂贵, 因此, 真实模拟飞机目标在光学系统的红外成像成为一个新的研究热点, 而其中飞机的红外辐射特性是其研究基础. 国内外学者对此已做了大量研究, 但尚未对飞机的结构进行建模, 且尾喷口及羽流的特征分析不够明确. 为此, 本文构建了飞机几何结构模型, 推导了飞机羽流辐射亮度的波数计算式, 分析了涡轮喷气发动机轴向温度分布, 考虑了羽流对尾喷口红外辐射的吸收和衰减.
收稿日期:2008-10-08
基金项目:航空科学基金资助项目([1**********])
大气状况对探测飞机目标有着至关重要的影响, 研究大气透过率不但有利于进一步真实模拟飞机目标红外成像, 而且是计算导弹作用距离的重要环节, 对空中作战有利. 目前, 国内外学者对大气透过率进行了大量研究, 且美国已开发了几套大气透过率计算软件, 如L o w t r a n 7
[1]
、MO D T -
R A N 等. 由于军事保密性, 其大量关键技术并未公开. 基于此现状, 本文实现了一种大气传输衰减简易计算模型, 并给出了其与L o w t r a n 7的计算结
) , 女, , m o u k y d c n
第10期 毛 峡等:飞机红外辐射及大气透过率计算方法1229
果对比图. A A 2A s i n θf u s =(b d +w n g ) a s p
观测角, 其大小定义为
θr c c o s [c o s θc o s θ]a s p =aa z e l
(6)
1 蒙皮的红外辐射特性
红外辐射是一种热辐射, 任何绝对温度大于
0K 的物体, 都会向外发射红外辐射, 物体在温度T 下辐射亮度的光谱分布
L (λ[2]
式中, θ为观测线与机身中心轴线的夹角, 又称a s p
(7)
其中, θ为观测线上的方位角; θ为观测线上的a z e l 俯仰角.
根据所求的L (λ) , A , 结合式(2) 、式(3) f u s f u s 两式便可得到蒙皮的辐射强度I . 本文算法为夜f u s 间作战情况, 未考虑太阳光对蒙皮红外辐射特性的影响.
可表示为
(1)
ε11
5
e x p (c T )-2/λ
式中, L (λ) 为光谱辐射亮度; c 1, c 2分别为第一、
第二辐射常数; ε为物体的发射率; λ为有效波长; T 为物体绝对温度, K .
由此, 可推得物体的光谱辐射强度以及辐射强度, 分别表示为
I (λ) =L(λ) A
λ2
2 羽流的红外辐射特性
飞机羽流的形状模型
[3]
如图2所示, 以x 0为
(2) (3)
I =I (λ) d λ
λ1
∫
式中, A 为辐射面积; I (λ) 为光谱辐射强度; I 为辐射强度; [λ1, λ2]为对应的波段.
当飞机在大气中高速飞行时, 飞机蒙皮由于受气动加热导致温度升高, 从而产生相应的红外辐射, 其大小可通过求驻点温度的方法求得. 其中驻点温度可由下式计算:
2
T {1+r [(γ-1) /2]M a }f u s =Ta i r
(4)
界, 可分为初始段以及主段两部分, 其中斜线部分
为核心区, 具有组分压强和温度基本恒定且为最大值的特性. 为计算沿任意方向羽流的辐射, 本文采用轴线上距尾喷口距离x 以及与轴线夹角A 两个参数确定, 具体步骤如下.
式中, T 为蒙皮表面驻点温度; T 为飞机飞行高f u s a i r 度上的大气温度; r 为恢复系数(层流取0. 82, 紊流取0. 87) ; γ为空气定压定容热容量之比(取1. 4) ; M a 为飞机飞行速度的马赫数. 由此, 结合式(1) 、式(4) 便可得到蒙皮的光谱辐射亮
度L (λ) . f u s
在计算蒙皮面积时, 飞机可大体分为机头、机身以及机翼三部分, 如图1所示.
图2 飞机羽流示意图
1) 将任意方向所截取的羽流部分分为温度
和组分压强均匀的层, 如图3所示, 每层的几何厚度分别为l m=1, 2, m , 温度为T m , 组分压强为p m (…) ;
图3 羽流温度和组分压强均匀层示意图
其中T 8) ~式(10) 求得m 及p m 可通过式(
1/24T p p m 2(8)
图1 飞机蒙皮结构图
[1-(1/n) ]22
n (y /y0)
T 0
T a
p βm m =pa (9) (10)
2
图1中, A 分别表示机头、机身、机h d , A b d , A w n g
翼的面积. 因此, 当0°≤θ90°时, 蒙皮的辐射a s p
A o s θ(A 2A s i n θ5)
f u s =Ah d c a s p +b d +w n g ) a s p (时, °
式中, T T 0为尾喷口截面处的羽流温度; a 为大气温度; p , n 为中间变量; p βm a 为一个标准大气压;
.
1230北京航空航天大学学报 2009年
2) 确定红外发射组分C O 的主要发2和H 2O 射带的波数范围;
3) 计算标态下各组分在每层的光学厚度:
u P (11) m m l m
T m
4) 根据光谱吸收系数表, 采用线性插值法,
查得各波数的光谱吸收系数K 为波数; m , ω, ω
5) 计算每层的光谱透射率, 得
τx p [-(K ·u ) (12) m , ω=em , ω] 6) 考虑谱线的碰撞展宽和多普勒展宽效应, 将光学深度(K ·u ) 修正为X m , ωm , ω, 得
τx p (-X m , ω=em , ω)
7) 第m 层的光谱辐射亮度为
L L ττm , ω=-ω, T (m , ω-m -1, ω) m W·c m ·s r , 表示为
3-16
c /π1·ω·10
L ω, T m
e x p (c ·10/T12·ωm )- 8) 计算沿任意方向的辐射亮度:
L p l =∑
m
-1
-1
的光谱辐射亮度L (λ) . h s
尾喷口的辐射面积A 可表示为h s
A h s 0 0°≤θa s p
πR c o s θ 90°≤θ≤180°h s a s p a s p
[4]
式中, R 为尾喷口的半径. 根据所求的L λ) 、h s h s (A 2) 、式(3) 便可得尾喷口的辐射强度h s , 结合式(I . h s
经羽流衰减后的尾喷口的辐射强度I ′为h s
I ′·τh s =Ih s p l
(19)
(13) (14)
4 大气透过率
大气对目标辐射亮度的影响基本表现在辐射衰减, 且可从以下三个方面进行研究. 4. 1 大气气体分子的吸收
研究表明, 物体吸收红外辐射的主要因素是水汽和C O 大气中水汽含量用可降水分w 度量, 2. 定义为底面S a 和长度等于红外系统到目标距离R 的圆柱体内大气含水汽凝结的水层厚度, 表示
为
w =w·H 0·R r
(20)
式中, R 为目标距离; H 为空气相对湿度; w r 0为一
式中, L ω, T 为黑体的光谱辐射亮度的波数表示式, m
(15)
∑L
ω
m , ω
·Δω(16)
式中Δω为波数间隔.
9) 计算羽流的平均透过率:
ττωp l m , ω·Δ∑m ∑m ω
(17)
定温度下、空气相对湿度为100%时, 每千米大气中的可降水分, 可通过查表得到.
此外, 考虑水汽吸收能力和水汽含量随高度的变化规律, 对目标距离R 进行如下修正:R ·e e =R
式中H 为飞机飞行的海拔高度.
-0. 5154H
3 尾喷口的红外辐射特性
涡轮喷气发动机的轴向温度分布如图4所
示, 当飞机工作于无加力状态时, 尾喷口温度主要由燃烧室决定; 在加力状态下, 燃气在加力燃烧室内剧烈燃烧, 致使尾喷口温度及其红外辐射急剧上升
.
(21)
综合式(20) 、式(21) , 得出w , 通过查海平面水平路程上水汽光谱透过率表, 采用外推法或内插法求得水汽光谱透过率τ(λ) . H 2O
同理, 考虑二氧化碳分子的吸收能力和质量随高度的变化, 对目标距离R 进行如下修正:
R ·e e =R
-0. 313H
[5]
(22)
CO (λ) 只与辐射通过的距2分子透过率τC O 2离有关, 根据τλ) 与R C O 2(e 对应的数据表, 可求得τ(λ) . C O 2
4. 2 大气中分子、气溶胶、微粒的散射
大气中传输的辐射亮度, 要受大气分子和微粒散射的影响, 纯粹由散射导致的透过率为
τλ) =e-3. 19·1(图4 涡轮喷气发动机轴向温度分布
q
(23)
式中, D λV 为能见度; 0代表大气的透射性能在可见光区的指定波长.
4. 3 气象条件所引起的衰减
、, 根据发动机的结构特性以及飞机的飞行状态得到尾喷口的温度, 结合式(1) 即可得到尾喷口
第10期 毛 峡等:飞机红外辐射及大气透过率计算方法1231
据米氏理论, 上述粒子将产生非选择性辐射散射. 研究表明, 对于决定与其强度相关的雨、雪的衰减系数10. 6μm 波长得到的经验公式:
0 J , J r s =0
(24) 0. 660. 7
0. 66J +6. 5J J 或J r s r s ≠0
式中, J , J 降雪强度, m m /h. r s 分别为降雨、
由雨、雪的衰减所导致的透过率为
τλ)=ex p (-a ·R ) (25) 2(
综上所述, 可求得大气光谱透过率为
τ(λ) =τλ) τλ) τλ) τλ) a i r H 2O (C O 2(1(2(
(26)
因此, 大气平均透过率表示为
1
ττ(λ) d λa i r a i r
λλ2-1λ
1
λ2
射只能从后方观测到. 在3~5μm 波段, 前向探测时, 红外辐射主要由羽流产生, 后向探测主要由尾喷口产生. 在8~14μm 波段, 羽流基本不产生辐
射, 蒙皮辐射成为主要辐射源. 因此, 在侧向或迎头探测时, 可以利用双波段探测器; 后向探测可主要利用3~5μm 波段的尾喷管和羽流辐射.
在计算大气透过率方面, 假定环境温度为294K , 相对湿度为50%,飞行高度7k m , 能见距离为10k m , 传输距离为垂直路径7k m , 无雨、雪天气. 在此基础上, 采用大量数据进行实验, 分别计算了不同波长、不同弹目距离下的大气透过率, 并将本文的计算结果与L o w t r a n 7的计算结果进行了比较, 通过多次比较, 可不断的对各系数以及查表方法进行修正, 以此得到较为精确的计算方法, 其结果如图7所示.
∫
(27)
5 结果分析
基于上述理论研究, 以某型涡轮喷气飞机为
例, 假定其飞行速度为
1. 5M a , 飞行高度为7k m , 环境温度为
294K , 尾喷口温度为1650K , 视距为0k m . 在此基础上, 计算了不同波段、不同角度下的蒙皮、尾喷口、羽流以及总的红外辐射强度, 其结果如图5和图6所示.
图7 不同波长下的大气透过率
由图7可知, 在0~20μm 波段, 存在多个大
气透过率窗口, 但是其中一些窗口由于太窄, 可忽略不计, 8~14μm 波段窗口是最宽的, 而当波长超过14μm 时, 辐射基本不能透过. 为此, 红外探测器工作在3~5μm 以及8~14μm 有利于检测飞机目标.
图8为图5、图6中总辐射强度经垂直路径7k m 大气衰减后, 在3~5μm 及8~14μm 波段到达地面红外探测器的辐射强度.
图5 飞机在3~5μm 波段的红外辐射强度
图6 飞机在8~14μm 波段的红外辐射强度
由图5和图6可知, 当红外探测器侧向或迎
, 图8 经大气衰减后的辐射强度
(下转第1272页)
1272北京航空航天大学学报 2009年
表5 优化前后各设计变量对比
p l i c a t i o n s[M].B e i j i n g :N a t i o n a lD e f e n s e I n d u s t r yP r e s s ,
β/(°) 1915. 614. 0
2006:48-58(i nC h i n e s e ) [4]M y e r s R H , M o n t g o m e r y dC .R e s p o n s es u r f a c em e t h o d o l o g y :
p r o c e s sa n d p r o d u c t o p t i m i z a t i o nu s i n gd e s i g n e d e x p e r i m e n t s [M ].N e wY o r k :Wi l e y , 1995:1278-1302[5]方开泰, 马长兴. 正交与均匀试验设计[M ]. 北京:科学出版
社, 2001:50-59
F a n gK a i t a i , M a C h a n g x i n g . R e s e a r c h o f o r t h o g o n a l a n du n i f o r m d e s i g n [M ].B e i j i n g :Sc i e n t i f i cP r e s s , 2001:50-59(i nC h i -n e s e ) [6]马东立, 武哲. 飞机翼面类部件的雷达散射截面计算[J ]. 系统工程与电子技术, 1994(6) :33-40
M a D o n g l i , WuZ h e .T h e c a l c u l a t i o no f r a d a r c r o s ss e c t i o nf o r a i r c r a f t a i r f o i l s [J ].Jo u r n a l o f S y s t e m s E n g i n e e r i n ga n dE l e c -t r o n i c s , 1994(6) :33-40(i n C h i n e s e ) [7]苏伟, 高正红, 夏露. 隐身性能约束的多目标气动外形优化设计[J ].空气动力学学报, 2006, 24(1) :137-140
S u We i , G a oZ h e n g h o n g , X i aL u .M u l t i o b j e c t i v eo p t i m i z a t i o n d e s i g no f a e r o d y n a m i cc o n f i g u r a t i o nc o n s t r a i n e db ys t e a l t hp e r -f o r m a n c e [J ].A c t a A e r o d y n a m i c a S i n i c a , 2006, 24(1) :137-140(i nC h i n e s e )
项目初始外形优化方案A 优化方案B
λf r 22. 93. 0
φ/(°) 2020. 727
d b /dm 0. 450. 370. 25
5 结 论
通过对非圆截面弹体进行气动隐身一体化优化设计, 结果表明, 本文所采用的基于二阶响应面
的P a r e t o 遗传算法能够较好地分析非圆截面弹体气动与隐身特性, 同时正确地反映了弹体各个设计变量对其气动及隐身特性的制约及响应关系, 优化结果比较合理. 证明了本方法适用于非圆截面弹体气动隐身一体化优化这类复杂的多目标设计问题.
参考文献(R e f e r e n c e s )
[1]M e n d e n h a l l M R . 战术导弹空气动力学[M ]. 2版. 洪金森译.
北京:宇航出版社, 1999:6-18
M e n d e n h a l l M R .T a c t i c a lm i s s i l e a e r o d y n a m i c s p r e d i c t i o n m e t h o d o g y [M ]. 2n de d .T r a n s l a t e db y H o n g J i n s e n . B e i j i n g :A s t r o n a u t i c a l l y P u b l i s h i n gH o u s e , 1999:6-18(i nC h i n e s e ) [2]C r o s s l e y W A , C o o kAM , F a n j o yDW, e t a l .U s i n gt h e t w o -b r a n c ht o u r n a m e n tg e n e t i ca l g o r i t h m f o rm u l t i o b j e c t i v ed e s i g n
[R ].AI A A -98-1914, 1998[3]崔逊学. 多目标进化算法及其应用[M ]. 北京:国防工业出版社, 2006:48-58
C u i X u n x u e . M u l t i o b j e c t i v e e v o l u t i o n a r y a l g o r i t h m s a n dt h e i r a p -
[8]S u nJ , H a n s e nSG .T r a n s o n i c n o s ew a v e d r a go f a x i s y m m e t r i c
b o d i e s [R ].AI A A -82-0308, 1982[9]何开锋, 钱炜祺, 刘刚, 等. 飞行器气动隐身一体化设计方法
研究[J ]. 空气动力学学报, 2006, 24(2) :169-174
H e K a i f e n g , Q i a nWe i q i , L i u G a n g , e t a l .A i r c r a f t m u l t i -o b j e c -t i v e d e s i g no f a e r o d y n a m i ca n ds t e a l t h yp e r f o r m a n c e [J ].A c t a A e r o d y n a m i c a S i n i c a , 2006, 24(2) :169-174(i n C h i n e s e ) [10]S i g n a l A , L a p i d o t E . T h e a e r o d y n a m i c c h a r a c t e r i s t i c o f c o n f i g -u r a t i o n sh a v i n gb o d i e sw i t hs q u a r er e c t a n g u l a ra n d c i r c u l a r
c r o s s s e c t i o n s a t a M a c h n u m b e r o f 0. 75[R ].AI A A -87-2429, 1987
(编 辑:李 晶)
(上接第1231页)
型[J ].北京航空航天大学学报, 2008, 34(8) :952-955
6 结 论
理论分析和计算机仿真结果表明, 本文对飞机蒙皮、羽流以及尾喷口在不同波段、不同角度下的计算结果的比例关系符合实际探测情况, 大气透过率的计算结果也与L o w t r a n 7的计算结果近似. 因此, 本文所提出的飞机红外辐射以及大气透过率计算模型具有一定的合理性. 本文对飞机红外图像仿真以及目标自动识别领域的研究有着重要的意义. 鉴于研究对象存在的诸多变化因子及不可准确测量参数, 本文计算结果的精确度有待进一步提高.
参考文献(R e f e r e n c e s )
[1]R a o GA , M a h u l i k a r S P .E f f e c t o f a t m o s p h e r i c t r a n s m i s s i o n a n d
r a d i a n c e o n a i r c r a f t i n f r a r e ds i g n a t u r e s [J ].J o u r n a l o f A i r c r a f t , 2005, 42(5) :1046-1054
[2, , M a o X i a , D o n g X u y a n g , D i a oWe i h e . C a l c u l a t i o nm o d e l o f i n -f r a r e dr a d i a t i o n o fc r u i s em i s s i l ea n da t m o s p h e r i ca t t e n u a t i o n [J ].J o u r n a l o f B e i j i n g U n i v e r s i t yo f A e r o n a u t i c s a n dA s t r o n a u -t i c s , 2008, 34(8) :952-955(i nC h i n e s e )
[3]高思莉, 汤心溢. 空中飞行目标尾焰红外辐射的建模与仿真
[J ].光电工程, 2007, 34(8) :25-27
G a o S i l i , T a n g X i n y i . B u i l d i n g m o d e l o f t h e p l u m e r e l e a s e df r o m t h ef l y i n gm a c h i n ea n ds i m u l a t i o n [J ].O p t o -E l e c t r o n i cE n g i -n e e r i n g , 2007, 34(8) :25-27(i nC h i n e s e )
[4]张海兴, 张建奇, 杨威. 飞机红外辐射的理论计算[J ]. 西安
电子科技大学学报, 1997, 24(1) :78-80
Z h a n g H a i x i n g , Z h a n gJ i a n q i , Y a n gWe i .T h e o r e t i c a l c a l c u l a -t i o n o f t h e I Rr a d i a t i o no f a n a e r o p l a n e [J ]. J o u r n a l o f X i d i a n U -n i v e r s i t y , 1997, 24(1) :78-80(i nC h i n e s e )
[5]吴晗平. 红外辐射大气透过率的工程理论方法研究[J ]. 光
学精密工程, 1998, 6(4) :35-42
W uH a n p i n g .R e s e a r c hi n t ot h e o r e t i c a l c a l c u l a t i o nm e t h o do n e n g i n e e r i n g o f t r a n s m i t t a n c e o f i n f r a r e dr a d i a t i o nt h r o u g ha t m o s -p h e r e [J ].Op t i c s a n dP r e c i s i o nE n g i n e e r i n g , 1998, 6(4) :35-42(i nC h i n e s e )
(编:娄
DOI :10. 13700/j . bh . 1001-5965. 2009. 10. 0152009年10月北京航空航天大学学报
第35卷第10期J o u r n a l o f B e i j i n g U n i v e r s i t y o f A e r o n a u t i c s a n d A s t r o n a u t i c s
O c t o b e r 2009
V o l . 35 No . 10
飞机红外辐射及大气透过率计算方法
毛 峡
胡海勇
黄 康
梁晓庚
(北京航空航天大学电子信息工程学院, 北京100191) (中国空空导弹研究院, 洛阳417009)
摘 要:研究了飞机目标的红外辐射特性, 以及红外辐射在大气中的辐射衰减特性. 根据现有的理论基础以及飞机自身结构特性, 将飞机的红外辐射源分为蒙皮、尾喷口和羽流3部分, 并提出了3个辐射源在不同波段、不同角度下的红外辐射特性计算方法. 根据红外辐射
的大气传输特性, 采用了不同波长、不同弹目距离下的大气透过率简易计算方法, 并将其与L o w t r a n 7进行比较. 计算了飞机在不同波段下, 经大气衰减后, 最终到达红外探测器的辐射强度.
关 键 词:红外辐射; 蒙皮; 尾喷口; 羽流; 大气透过率; 飞机
中图分类号:T N 219
文献标识码:A 文章编号:1001-5965(2009) 10-1228-04
C a l c u l a t i o nm e t h o d f o r a i r p l a n eI Rr a d i a t i o na n d a t m o s p h e r i c t r a n s m i t t a n c e
M a o X i a Hu H a i y o n g Hu a n g K a n g
(S c h o o l o f E l e c t r o n i c s a n dI n f o r m a t i o n E n g i n e e r i n g , B e i j i n gU n i v e r s i t y o f A e r o n a u t i c s a n dA s t r o n a u t i c s , B e i j i n g 100191, C h i n a )
L i a n g X i a o g e n g
(C h i n a A i r b o r n e M i s s i l eA c a d e m i c , L u o y a n g 417009, C h i n a )
A b s t r a c t :T h e i n f r a r e d r a d i a t i o n o f a i r p l a n e t a r g e t s a n d t h e r a d i a t i o n a t t e n u a t i o n o f t h e i n f r a r e d r a d i a t i o n
i n t h e a t m o s p h e r e w e r e r e s e a r c h e d . B a s e d o n t h e e x i s t i n g t h e o r e t i c a l s t u d y a n d t h e s t r u c t u r e o f t h e a i r p l a n e , t h e i n f r a r e d r a d i a t i o n s o u r c e s o f t h e a i r p l a n e w e r e d i v i d e d i n t o t h r e e c a t e g o r i e s :f u s e l a g e , h o t s p o t a n d p l u m e . F u r t h e r m o r e , c a l c u l a t i o n a p p r o a c h e s t o t h e c h a r a c t e r o f i n f r a r e d r a d i a t i o n f o r t h e t h r e e r a d i a t i o n s o u r c e s a t d i f -f e r e n t w a v e b a n d a n d d i f f e r e n t a n g l e w e r e p r o p o s e d . A c c o r d i n g t o i n f r a r e d r a d i a t i o n t r a n s m i t t a n c e i n t h e a t m o s -p h e r e , a c a l c u l a t i o n m e t h o d f o r t h e t r a n s m i t t a n c e o f a t m o s p h e r e w a s p u t f o r w a r d a t d i f f e r e n t w a v e l e n g t h a n d d i f f e r e n t t a r g e t d i s t a n c e , a n d t h e n t h i s m e t h o d w a s c o m p a r e d w i t h L o w t r a n 7. F i n a l l y , t h e r a d i a n t i n t e n s i t y r e -c e i v e d b y t h e I Rd e t e c t o r w a s c o m p u t e d a f t e r a t m o s p h e r i c a t t e n u a t i o n a t d i f f e r e n t w a v e b a n d .
K e y w o r d s :i n f r a r e d r a d i a t i o n ; f u s e l a g e ; h o t s p o t ; p l u m e ; a t m o s p h e r i c t r a n s m i t t a n c e ; a i r p l a n e
外场实验获取飞机目标红外图像的费用昂贵, 因此, 真实模拟飞机目标在光学系统的红外成像成为一个新的研究热点, 而其中飞机的红外辐射特性是其研究基础. 国内外学者对此已做了大量研究, 但尚未对飞机的结构进行建模, 且尾喷口及羽流的特征分析不够明确. 为此, 本文构建了飞机几何结构模型, 推导了飞机羽流辐射亮度的波数计算式, 分析了涡轮喷气发动机轴向温度分布, 考虑了羽流对尾喷口红外辐射的吸收和衰减.
收稿日期:2008-10-08
基金项目:航空科学基金资助项目([1**********])
大气状况对探测飞机目标有着至关重要的影响, 研究大气透过率不但有利于进一步真实模拟飞机目标红外成像, 而且是计算导弹作用距离的重要环节, 对空中作战有利. 目前, 国内外学者对大气透过率进行了大量研究, 且美国已开发了几套大气透过率计算软件, 如L o w t r a n 7
[1]
、MO D T -
R A N 等. 由于军事保密性, 其大量关键技术并未公开. 基于此现状, 本文实现了一种大气传输衰减简易计算模型, 并给出了其与L o w t r a n 7的计算结
) , 女, , m o u k y d c n
第10期 毛 峡等:飞机红外辐射及大气透过率计算方法1229
果对比图. A A 2A s i n θf u s =(b d +w n g ) a s p
观测角, 其大小定义为
θr c c o s [c o s θc o s θ]a s p =aa z e l
(6)
1 蒙皮的红外辐射特性
红外辐射是一种热辐射, 任何绝对温度大于
0K 的物体, 都会向外发射红外辐射, 物体在温度T 下辐射亮度的光谱分布
L (λ[2]
式中, θ为观测线与机身中心轴线的夹角, 又称a s p
(7)
其中, θ为观测线上的方位角; θ为观测线上的a z e l 俯仰角.
根据所求的L (λ) , A , 结合式(2) 、式(3) f u s f u s 两式便可得到蒙皮的辐射强度I . 本文算法为夜f u s 间作战情况, 未考虑太阳光对蒙皮红外辐射特性的影响.
可表示为
(1)
ε11
5
e x p (c T )-2/λ
式中, L (λ) 为光谱辐射亮度; c 1, c 2分别为第一、
第二辐射常数; ε为物体的发射率; λ为有效波长; T 为物体绝对温度, K .
由此, 可推得物体的光谱辐射强度以及辐射强度, 分别表示为
I (λ) =L(λ) A
λ2
2 羽流的红外辐射特性
飞机羽流的形状模型
[3]
如图2所示, 以x 0为
(2) (3)
I =I (λ) d λ
λ1
∫
式中, A 为辐射面积; I (λ) 为光谱辐射强度; I 为辐射强度; [λ1, λ2]为对应的波段.
当飞机在大气中高速飞行时, 飞机蒙皮由于受气动加热导致温度升高, 从而产生相应的红外辐射, 其大小可通过求驻点温度的方法求得. 其中驻点温度可由下式计算:
2
T {1+r [(γ-1) /2]M a }f u s =Ta i r
(4)
界, 可分为初始段以及主段两部分, 其中斜线部分
为核心区, 具有组分压强和温度基本恒定且为最大值的特性. 为计算沿任意方向羽流的辐射, 本文采用轴线上距尾喷口距离x 以及与轴线夹角A 两个参数确定, 具体步骤如下.
式中, T 为蒙皮表面驻点温度; T 为飞机飞行高f u s a i r 度上的大气温度; r 为恢复系数(层流取0. 82, 紊流取0. 87) ; γ为空气定压定容热容量之比(取1. 4) ; M a 为飞机飞行速度的马赫数. 由此, 结合式(1) 、式(4) 便可得到蒙皮的光谱辐射亮
度L (λ) . f u s
在计算蒙皮面积时, 飞机可大体分为机头、机身以及机翼三部分, 如图1所示.
图2 飞机羽流示意图
1) 将任意方向所截取的羽流部分分为温度
和组分压强均匀的层, 如图3所示, 每层的几何厚度分别为l m=1, 2, m , 温度为T m , 组分压强为p m (…) ;
图3 羽流温度和组分压强均匀层示意图
其中T 8) ~式(10) 求得m 及p m 可通过式(
1/24T p p m 2(8)
图1 飞机蒙皮结构图
[1-(1/n) ]22
n (y /y0)
T 0
T a
p βm m =pa (9) (10)
2
图1中, A 分别表示机头、机身、机h d , A b d , A w n g
翼的面积. 因此, 当0°≤θ90°时, 蒙皮的辐射a s p
A o s θ(A 2A s i n θ5)
f u s =Ah d c a s p +b d +w n g ) a s p (时, °
式中, T T 0为尾喷口截面处的羽流温度; a 为大气温度; p , n 为中间变量; p βm a 为一个标准大气压;
.
1230北京航空航天大学学报 2009年
2) 确定红外发射组分C O 的主要发2和H 2O 射带的波数范围;
3) 计算标态下各组分在每层的光学厚度:
u P (11) m m l m
T m
4) 根据光谱吸收系数表, 采用线性插值法,
查得各波数的光谱吸收系数K 为波数; m , ω, ω
5) 计算每层的光谱透射率, 得
τx p [-(K ·u ) (12) m , ω=em , ω] 6) 考虑谱线的碰撞展宽和多普勒展宽效应, 将光学深度(K ·u ) 修正为X m , ωm , ω, 得
τx p (-X m , ω=em , ω)
7) 第m 层的光谱辐射亮度为
L L ττm , ω=-ω, T (m , ω-m -1, ω) m W·c m ·s r , 表示为
3-16
c /π1·ω·10
L ω, T m
e x p (c ·10/T12·ωm )- 8) 计算沿任意方向的辐射亮度:
L p l =∑
m
-1
-1
的光谱辐射亮度L (λ) . h s
尾喷口的辐射面积A 可表示为h s
A h s 0 0°≤θa s p
πR c o s θ 90°≤θ≤180°h s a s p a s p
[4]
式中, R 为尾喷口的半径. 根据所求的L λ) 、h s h s (A 2) 、式(3) 便可得尾喷口的辐射强度h s , 结合式(I . h s
经羽流衰减后的尾喷口的辐射强度I ′为h s
I ′·τh s =Ih s p l
(19)
(13) (14)
4 大气透过率
大气对目标辐射亮度的影响基本表现在辐射衰减, 且可从以下三个方面进行研究. 4. 1 大气气体分子的吸收
研究表明, 物体吸收红外辐射的主要因素是水汽和C O 大气中水汽含量用可降水分w 度量, 2. 定义为底面S a 和长度等于红外系统到目标距离R 的圆柱体内大气含水汽凝结的水层厚度, 表示
为
w =w·H 0·R r
(20)
式中, R 为目标距离; H 为空气相对湿度; w r 0为一
式中, L ω, T 为黑体的光谱辐射亮度的波数表示式, m
(15)
∑L
ω
m , ω
·Δω(16)
式中Δω为波数间隔.
9) 计算羽流的平均透过率:
ττωp l m , ω·Δ∑m ∑m ω
(17)
定温度下、空气相对湿度为100%时, 每千米大气中的可降水分, 可通过查表得到.
此外, 考虑水汽吸收能力和水汽含量随高度的变化规律, 对目标距离R 进行如下修正:R ·e e =R
式中H 为飞机飞行的海拔高度.
-0. 5154H
3 尾喷口的红外辐射特性
涡轮喷气发动机的轴向温度分布如图4所
示, 当飞机工作于无加力状态时, 尾喷口温度主要由燃烧室决定; 在加力状态下, 燃气在加力燃烧室内剧烈燃烧, 致使尾喷口温度及其红外辐射急剧上升
.
(21)
综合式(20) 、式(21) , 得出w , 通过查海平面水平路程上水汽光谱透过率表, 采用外推法或内插法求得水汽光谱透过率τ(λ) . H 2O
同理, 考虑二氧化碳分子的吸收能力和质量随高度的变化, 对目标距离R 进行如下修正:
R ·e e =R
-0. 313H
[5]
(22)
CO (λ) 只与辐射通过的距2分子透过率τC O 2离有关, 根据τλ) 与R C O 2(e 对应的数据表, 可求得τ(λ) . C O 2
4. 2 大气中分子、气溶胶、微粒的散射
大气中传输的辐射亮度, 要受大气分子和微粒散射的影响, 纯粹由散射导致的透过率为
τλ) =e-3. 19·1(图4 涡轮喷气发动机轴向温度分布
q
(23)
式中, D λV 为能见度; 0代表大气的透射性能在可见光区的指定波长.
4. 3 气象条件所引起的衰减
、, 根据发动机的结构特性以及飞机的飞行状态得到尾喷口的温度, 结合式(1) 即可得到尾喷口
第10期 毛 峡等:飞机红外辐射及大气透过率计算方法1231
据米氏理论, 上述粒子将产生非选择性辐射散射. 研究表明, 对于决定与其强度相关的雨、雪的衰减系数10. 6μm 波长得到的经验公式:
0 J , J r s =0
(24) 0. 660. 7
0. 66J +6. 5J J 或J r s r s ≠0
式中, J , J 降雪强度, m m /h. r s 分别为降雨、
由雨、雪的衰减所导致的透过率为
τλ)=ex p (-a ·R ) (25) 2(
综上所述, 可求得大气光谱透过率为
τ(λ) =τλ) τλ) τλ) τλ) a i r H 2O (C O 2(1(2(
(26)
因此, 大气平均透过率表示为
1
ττ(λ) d λa i r a i r
λλ2-1λ
1
λ2
射只能从后方观测到. 在3~5μm 波段, 前向探测时, 红外辐射主要由羽流产生, 后向探测主要由尾喷口产生. 在8~14μm 波段, 羽流基本不产生辐
射, 蒙皮辐射成为主要辐射源. 因此, 在侧向或迎头探测时, 可以利用双波段探测器; 后向探测可主要利用3~5μm 波段的尾喷管和羽流辐射.
在计算大气透过率方面, 假定环境温度为294K , 相对湿度为50%,飞行高度7k m , 能见距离为10k m , 传输距离为垂直路径7k m , 无雨、雪天气. 在此基础上, 采用大量数据进行实验, 分别计算了不同波长、不同弹目距离下的大气透过率, 并将本文的计算结果与L o w t r a n 7的计算结果进行了比较, 通过多次比较, 可不断的对各系数以及查表方法进行修正, 以此得到较为精确的计算方法, 其结果如图7所示.
∫
(27)
5 结果分析
基于上述理论研究, 以某型涡轮喷气飞机为
例, 假定其飞行速度为
1. 5M a , 飞行高度为7k m , 环境温度为
294K , 尾喷口温度为1650K , 视距为0k m . 在此基础上, 计算了不同波段、不同角度下的蒙皮、尾喷口、羽流以及总的红外辐射强度, 其结果如图5和图6所示.
图7 不同波长下的大气透过率
由图7可知, 在0~20μm 波段, 存在多个大
气透过率窗口, 但是其中一些窗口由于太窄, 可忽略不计, 8~14μm 波段窗口是最宽的, 而当波长超过14μm 时, 辐射基本不能透过. 为此, 红外探测器工作在3~5μm 以及8~14μm 有利于检测飞机目标.
图8为图5、图6中总辐射强度经垂直路径7k m 大气衰减后, 在3~5μm 及8~14μm 波段到达地面红外探测器的辐射强度.
图5 飞机在3~5μm 波段的红外辐射强度
图6 飞机在8~14μm 波段的红外辐射强度
由图5和图6可知, 当红外探测器侧向或迎
, 图8 经大气衰减后的辐射强度
(下转第1272页)
1272北京航空航天大学学报 2009年
表5 优化前后各设计变量对比
p l i c a t i o n s[M].B e i j i n g :N a t i o n a lD e f e n s e I n d u s t r yP r e s s ,
β/(°) 1915. 614. 0
2006:48-58(i nC h i n e s e ) [4]M y e r s R H , M o n t g o m e r y dC .R e s p o n s es u r f a c em e t h o d o l o g y :
p r o c e s sa n d p r o d u c t o p t i m i z a t i o nu s i n gd e s i g n e d e x p e r i m e n t s [M ].N e wY o r k :Wi l e y , 1995:1278-1302[5]方开泰, 马长兴. 正交与均匀试验设计[M ]. 北京:科学出版
社, 2001:50-59
F a n gK a i t a i , M a C h a n g x i n g . R e s e a r c h o f o r t h o g o n a l a n du n i f o r m d e s i g n [M ].B e i j i n g :Sc i e n t i f i cP r e s s , 2001:50-59(i nC h i -n e s e ) [6]马东立, 武哲. 飞机翼面类部件的雷达散射截面计算[J ]. 系统工程与电子技术, 1994(6) :33-40
M a D o n g l i , WuZ h e .T h e c a l c u l a t i o no f r a d a r c r o s ss e c t i o nf o r a i r c r a f t a i r f o i l s [J ].Jo u r n a l o f S y s t e m s E n g i n e e r i n ga n dE l e c -t r o n i c s , 1994(6) :33-40(i n C h i n e s e ) [7]苏伟, 高正红, 夏露. 隐身性能约束的多目标气动外形优化设计[J ].空气动力学学报, 2006, 24(1) :137-140
S u We i , G a oZ h e n g h o n g , X i aL u .M u l t i o b j e c t i v eo p t i m i z a t i o n d e s i g no f a e r o d y n a m i cc o n f i g u r a t i o nc o n s t r a i n e db ys t e a l t hp e r -f o r m a n c e [J ].A c t a A e r o d y n a m i c a S i n i c a , 2006, 24(1) :137-140(i nC h i n e s e )
项目初始外形优化方案A 优化方案B
λf r 22. 93. 0
φ/(°) 2020. 727
d b /dm 0. 450. 370. 25
5 结 论
通过对非圆截面弹体进行气动隐身一体化优化设计, 结果表明, 本文所采用的基于二阶响应面
的P a r e t o 遗传算法能够较好地分析非圆截面弹体气动与隐身特性, 同时正确地反映了弹体各个设计变量对其气动及隐身特性的制约及响应关系, 优化结果比较合理. 证明了本方法适用于非圆截面弹体气动隐身一体化优化这类复杂的多目标设计问题.
参考文献(R e f e r e n c e s )
[1]M e n d e n h a l l M R . 战术导弹空气动力学[M ]. 2版. 洪金森译.
北京:宇航出版社, 1999:6-18
M e n d e n h a l l M R .T a c t i c a lm i s s i l e a e r o d y n a m i c s p r e d i c t i o n m e t h o d o g y [M ]. 2n de d .T r a n s l a t e db y H o n g J i n s e n . B e i j i n g :A s t r o n a u t i c a l l y P u b l i s h i n gH o u s e , 1999:6-18(i nC h i n e s e ) [2]C r o s s l e y W A , C o o kAM , F a n j o yDW, e t a l .U s i n gt h e t w o -b r a n c ht o u r n a m e n tg e n e t i ca l g o r i t h m f o rm u l t i o b j e c t i v ed e s i g n
[R ].AI A A -98-1914, 1998[3]崔逊学. 多目标进化算法及其应用[M ]. 北京:国防工业出版社, 2006:48-58
C u i X u n x u e . M u l t i o b j e c t i v e e v o l u t i o n a r y a l g o r i t h m s a n dt h e i r a p -
[8]S u nJ , H a n s e nSG .T r a n s o n i c n o s ew a v e d r a go f a x i s y m m e t r i c
b o d i e s [R ].AI A A -82-0308, 1982[9]何开锋, 钱炜祺, 刘刚, 等. 飞行器气动隐身一体化设计方法
研究[J ]. 空气动力学学报, 2006, 24(2) :169-174
H e K a i f e n g , Q i a nWe i q i , L i u G a n g , e t a l .A i r c r a f t m u l t i -o b j e c -t i v e d e s i g no f a e r o d y n a m i ca n ds t e a l t h yp e r f o r m a n c e [J ].A c t a A e r o d y n a m i c a S i n i c a , 2006, 24(2) :169-174(i n C h i n e s e ) [10]S i g n a l A , L a p i d o t E . T h e a e r o d y n a m i c c h a r a c t e r i s t i c o f c o n f i g -u r a t i o n sh a v i n gb o d i e sw i t hs q u a r er e c t a n g u l a ra n d c i r c u l a r
c r o s s s e c t i o n s a t a M a c h n u m b e r o f 0. 75[R ].AI A A -87-2429, 1987
(编 辑:李 晶)
(上接第1231页)
型[J ].北京航空航天大学学报, 2008, 34(8) :952-955
6 结 论
理论分析和计算机仿真结果表明, 本文对飞机蒙皮、羽流以及尾喷口在不同波段、不同角度下的计算结果的比例关系符合实际探测情况, 大气透过率的计算结果也与L o w t r a n 7的计算结果近似. 因此, 本文所提出的飞机红外辐射以及大气透过率计算模型具有一定的合理性. 本文对飞机红外图像仿真以及目标自动识别领域的研究有着重要的意义. 鉴于研究对象存在的诸多变化因子及不可准确测量参数, 本文计算结果的精确度有待进一步提高.
参考文献(R e f e r e n c e s )
[1]R a o GA , M a h u l i k a r S P .E f f e c t o f a t m o s p h e r i c t r a n s m i s s i o n a n d
r a d i a n c e o n a i r c r a f t i n f r a r e ds i g n a t u r e s [J ].J o u r n a l o f A i r c r a f t , 2005, 42(5) :1046-1054
[2, , M a o X i a , D o n g X u y a n g , D i a oWe i h e . C a l c u l a t i o nm o d e l o f i n -f r a r e dr a d i a t i o n o fc r u i s em i s s i l ea n da t m o s p h e r i ca t t e n u a t i o n [J ].J o u r n a l o f B e i j i n g U n i v e r s i t yo f A e r o n a u t i c s a n dA s t r o n a u -t i c s , 2008, 34(8) :952-955(i nC h i n e s e )
[3]高思莉, 汤心溢. 空中飞行目标尾焰红外辐射的建模与仿真
[J ].光电工程, 2007, 34(8) :25-27
G a o S i l i , T a n g X i n y i . B u i l d i n g m o d e l o f t h e p l u m e r e l e a s e df r o m t h ef l y i n gm a c h i n ea n ds i m u l a t i o n [J ].O p t o -E l e c t r o n i cE n g i -n e e r i n g , 2007, 34(8) :25-27(i nC h i n e s e )
[4]张海兴, 张建奇, 杨威. 飞机红外辐射的理论计算[J ]. 西安
电子科技大学学报, 1997, 24(1) :78-80
Z h a n g H a i x i n g , Z h a n gJ i a n q i , Y a n gWe i .T h e o r e t i c a l c a l c u l a -t i o n o f t h e I Rr a d i a t i o no f a n a e r o p l a n e [J ]. J o u r n a l o f X i d i a n U -n i v e r s i t y , 1997, 24(1) :78-80(i nC h i n e s e )
[5]吴晗平. 红外辐射大气透过率的工程理论方法研究[J ]. 光
学精密工程, 1998, 6(4) :35-42
W uH a n p i n g .R e s e a r c hi n t ot h e o r e t i c a l c a l c u l a t i o nm e t h o do n e n g i n e e r i n g o f t r a n s m i t t a n c e o f i n f r a r e dr a d i a t i o nt h r o u g ha t m o s -p h e r e [J ].Op t i c s a n dP r e c i s i o nE n g i n e e r i n g , 1998, 6(4) :35-42(i nC h i n e s e )
(编:娄