. 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
对于介质体散射模型 ,将介质体 (损耗介质 )可以等效为在外形上强加一个等效阻抗 ,从而应用阻抗边界条件来进行计算。对于图 4中两介质体 (介质体 1为自由空间 )的等效阻抗计算公式表达为 [ 13214 ]
式中 :d为介质体 2的厚度 ;k1为自由空间波数 ,具体算式见式 (1)注释 ;ηi(i =1 , 2)为介质体 i的
归一化特征阻抗 ,算式定义为
i/εi。
当介质 2为非理想导体 ,其边界面上的阻抗边界条件可表示为 [15 ]
E-(n1 E)n = ηZ0 n ·H(8) 对于 Ez极化的二维情形 ,可具体表示为
9 Ez = j k0 μ1 Ez (9)9n η2
两式中的符号含义见式 (1)。对于式 (9)的边界条件可在式 (5)及式 (6)置入并进行数值计算。表 2中给出了其缩比模型与原始模型在不同阻抗下的 RCS数值计算值及其相对误差值。
从表 2可以看出来 ,随着阻抗变大 ,介质体外形缩比模型采用导体外形缩比模型测量原则计算得到的初始模型 RCS值与原始模型相比较 ,其相对误差也越来越大。不过当阻抗值在 01 5以下时 ,相对误差还是可以容许的。
表 2 介质体翼型原始模型与缩比模型的 RCS计算值 Table 2 RCS numerical results of different scaling2down medium airfoil models
阻抗 /Ω 11 000 01 01 447 21
模型 /频率 ( GHz) RCS/ m Error/ % RCS/ m Error/ %
1 ∶1/ 5 1 ∶2/ 10 1 ∶3/ 15 1 ∶4/ 20 阻抗 /Ω 01 041 755 01 018 918 -91 385 7 01 014 585 41 789 8 01 011 258 71 848 2 01 316 23 01 047 665 01 022 691 -41789 7 01 016 206 11999 4 01 012 201 21389 601 223 61
模型 /频率 ( GHz) RCS/ m Error/ % RCS/ m Error/ %
1 ∶1/ 5 1 ∶2/ 10 1 ∶3/ 15 1 ∶4/ 20 01 049 173 01 023 577 01 016 618 01 012 44 -41 105 9 11 384 9 11 193 7 01 050 214 01 024 168 01 016 9 01 012 603 -31 74 01 967 86 01 394 31
(7)
η2 + jη1tan k1 d
η1 + jη1tan k1 d
4 结 论
基于二维电磁散射的有限元数值计算原理 ,在验证了计算单元尺度收敛性的工作基础上 ,首先对导体缩比模型的散射规律进行了数值计算验
参 考 文 献
[1] 张卫东 ,冯小云 ,孟秀兰 .国外隐身材料研究进展 [J ].宇航材料工艺 ,2000 (3) :124. Zhang Weidong , Feng Xiaoyun , Meng Xiulan. The ad2 vance of foreign stealth materials[J ]. Aerospace Materials and Technology , 2000 (3) : 124. (in Chinese)
[2] 李江海 ,孙秦 .结构型吸波材料及其结构型式设计研究进展[J ].机械科学与技术 ,2003 ,22 (S2) :1882190. Li Jianghai , Sun Qin. The development of structured wave absorbing materials and configuration design [J ]. Mechanical Science and Technology, 2003,22 (S2) :1882
190. (in Chinese)
[3] 郑天亮 ,张璋 ,王轩 ,等.聚苯胺中空微球的改性与电磁特性[J ].航空学报 ,2007 ,28 (6) :153221535. Zheng Tianliang , Zhang Zhang , Wang Xuan ,et al. Modi2 fying of PANI hollow microspheres and electromagnetic property[J]. Acta Aeronautica et Astronautica sinica, 2007 ,28 (6) :153221535. (in Chinese)
[4] Johansson M , Holloway C L , Kuester E F. Effective elec2 tromagnetic properties of honeycomb composites , and hol2 low2pyramidal and alternation2wedge absorbers[J ]. IEEE Trans Antennas Propag , 2005 ,53 (2) :7282735.
[5] Graglia R D , Wilton D R , Peterson A F. Higher order in2 terpolatory vector bases for computational electromagnet2 ics[J ]. IEEE Trans Antennas Propagat , 1997 , 45 (3) : 3292342.
[6] Ilic’M M , Notaro.B M. Higher order hierarchical curved hexahedra vector finite elements for electromagnetic mod2 eling[J ]. IEEE Trans Microwave Theory Tech , 2003 , 51
(3) :102621033.
[7] Hadjali M , Besbes M , Bouillault F. Modeling of magneto2 mechanical phenomena by using shell elements[J ]. IEEE Trans Magn , 2004 ,40 (2) :5692572.
[8] Stupfel B , Mittra R. Numerical absorbing boundary con2 ditions for the scalar and vector wave equations [J ]. IEEE Trans Antennas Propagat , 1996 ,44 (7) :101521022.
[9] OmarMR. Stabilityofbsorbingboundaryconditions[J].
. 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
IEEE Trans Antennas Propagat , 1999 , 47 (4) :5932599.
[ 10 ] ChatterjeeA,VolakisJL. Conformalabsorbingboundary conditions for 32D problems: derivations and applications
[J]. IEEETransAntennas Propagat,1995,43 (8):8602
866.
[ 11 ] Peterson A F. Absorbing boundary conditions for the vec2 torwaveequation[J]. MicrowaveOpt TechLett,1988,1: 62264. 中国航空网 www.aero.cn 航空翻译 www.aviation.cn 本文链接地址:航空学报08大飞机专刊(120)