飞行器座舱RCS可视化计算方法研究

图1　(a)座舱插值拟合前网栅图　(b)座舱插值拟合后网栅图

　　通过C-R样条曲面对目标模型的拟合，实现了复杂目标由型值点构成的多边形粗糙模型到光滑真实模型的过滤，在对目标进行几何描述以后，即可应用图形软件标准接口(OpenGL)［4］和图形加速卡硬件对目标进行显示和消隐，从而在微机上实现GRECO的电磁计算.以下是用图像生成程序显示的座舱模型图2(a)，(b)前者是没插值前的图形，后者是用C-R样条拟合后的图形.

图2　(a)拟合前座舱模型　　(b)拟合后座舱模型

其详细计算实现方法请看参考文献［5］.

三、舱内结构散射分析
　　当雷达信号进入座舱，经多次反射后又传出座舱，其传播方向无法预料，其强度也很难估计.这里，将用以下方法解决该问题.首先应用分层媒质理论［1］求解舱罩的反射系数和透射系数.
　　1.分层媒质的反射和透射系数［1］
　　令垂直于Z轴的多层媒质的层数为L，加之其外部边界，共有L+2种媒质.分层媒质界面垂直于Z轴，则在入射媒质与分层媒质中的电磁场可表示为：

　(1)

式中为沿轴的单位矢量，［Mm］的具体表示式为

　(2)

式中

δm=(2πNm/λ)dmcosθm
Nm=εmμm-j(4πσmμm/ω)　(3)

其中λ为入射平面波波长，ω为角频率，ε、μ和σ分别为媒质的介电常数、导磁率和导电率，为简便起见，定义分层结构的导纳Y为

　(4)

故方程(1)可表示为

　(5)
　(6)

式中［B　C］T定义为分层结构的特征矩阵，且

Y=C/B　(7)

　　实际上，分层媒质的反射系数、传输系数和吸收系数可表示为

　(8)
　(9)
　(10)

　　2.能量分布调制法［6］
　　鉴于雷达波穿过舱罩进入舱内，由舱内散射体散射到舱外空间过程中，至少两次穿透舱罩结构，利用分层媒质散射矩阵和能量、相位加权，考虑到舱内射线物理过程产生的随机性，在舱内后向散射的RCS求解中引入了随机因素.利用随机因数生成程序+应用程序，从而获得良好的统计结果.
　　假设透过舱罩进入座舱的能量为ε，由于舱内人体、头盔、座椅以及仪表框架等物体的散射，该能量被随机地在方位角(φ0,φ0)和俯仰角(θu,θd)范围内散布.相当于以某一能量分布函数F(θ，φ)对均匀扩散情况下的平均能量进行调制加权，使能量分布与实际情况更逼近.F(θ，φ)可视不同机种的情况，通过分析和测试予以确定.在此种情况下，舱内电磁能量密度的分布可表示为：

(θ,φ)=εF(θ,φ)/∫θdθu∫θ0-θ0R2sinθdφdθ
=εF(θ,φ)/［2R2φ0(cosθu-cosθd)］　(11)

因此，在某方向(θ，φ)上由引起的RCS值为：

σ(θ,φ)=lim［4πR2(θ,φ)/|Ei|2］
=2πεF(θ,φ)/［φ0(cosθu-cosθd)］　(12)

考虑到雷达波经过舱罩进入舱内，由舱内物体散射回舱外空间的过程中，两次穿透舱罩结构，势必产生能量损耗，于是

(θ,φ)=εF(θ,φ).β/［2R2φ0(cosθu-cosθd)］
σ(θ,φ)=2πεF(θ,φ).β/［φ0(cosθu-cosθd)］ 　(13)

式中β为衰减因子，且β正比于舱罩透射系数的平方.式中F(θ，φ)必须满足

∫θdθu∫φd-φ0F(θ,φ)sinθdθdφ=2φ0(cosθu-cosθd)　(14)

由于

∫θdθu∫φd-φ0(θ,φ)R2sinθdθdφ=ε

令

则

σ(θ,φ)=4πεβF0(θ,φ)/∫∫F0(θ,φ)sinθdθdφ　(15)

　　式(15)中F0(θ，φ)的选取应根据统计结果确定.例如，对均匀分布而言F0(θ,φ)=1，而对高斯分布和对数分布情况下，分别为FG0(θ，φ)和FL0(θ,φ):

式中ξ和α为分布参数.在确定F0(θ，φ)之后,式(15)可用来解舱内散射对某给定方向的RCS值σ(θ，φ).舱内结构散射以及舱外金属面部分构成总的面效应场Esf，棱边部分则构成总的边缘散射场Esw.

四、座舱总散射场
　　舱内结构散射以及舱外金属面部分构成总的面效应场Esf,棱边部分则构成总的边缘散射场Esw.其RCS值为

　(16)

五、数值结果分析
　　对于带有镀层厚度在40～80Å的有机玻璃(其厚度为10mm，相对介电常数和导电率各为εr=2.62,μr=1)，薄膜的导电率取为5×106(1/Ω.m)，则在垂直入射时电磁反射率与镀层厚度的关系如图3和图4所示.

图3　反射率与镀层厚度的关系(重直极化)

图4　反射率与镀层厚度的关系(水平极化)

　　图5(a)为、10.5GHz、VV极化带导电镀层(100Å)、座舱姿态角为(φ,0,0)时的飞机座舱的理论计算曲线，相应的实验曲线图5(b).

图5　目标：座舱模型，频率：10.5GHz，极化：VV，带导电镀层

　　图6(a)为10.5GHz、HH极化、带导电镀层(100Å)、座舱姿态角为(φ，0，0)时的飞机座舱理论计算曲线，相应的实验曲线图(b)

图6　目标：座舱模型，频率：10.5GHz，极化：HH，带导电镀层

　　从实验和计算结果可以看出，采用镀膜结构的座舱RCS平均值为-15dBSM，两者平均值一致，总的趋势一致，证明了该算法的有效性.另外采用镀膜结构的座舱透入舱的能量较小，再返回舱外的散射能量更小，远小于无镀膜的座舱.利用以上计算方法可给出不同飞行状态，极化，和工作频率的RCS值，它可作为航空飞行器CAD/CAM系统中的一个基本模块，与其它CAD模块一起，进行优化处理，为飞行器设计提供依据.