Application of JEMS-FDTD in EMP regional propagation simulation
-
摘要: 电磁脉冲区域传播数值模拟是电磁环境效应分析的重要环节, 面临空间尺度巨大、多辐射源、复杂地貌等技术挑战。本文介绍了三维时域全波电磁模拟并行软件JEMS-FDTD的研制进展。在大规模并行计算FDTD方法的基础上结合自适应网格技术, 研制了一种适应于大区域电磁脉冲传播的高效时域全波计算方法。并基于该技术实现了千km2级城市电磁脉冲区域传播的数值模拟, 获取了全空间的时域电磁场信息, 验证了软件在电磁脉冲区域传播仿真应用中的可行性。Abstract: Numerical simulation of electromagnetic pluses regional propagation is an important link in the analysis of electromagnetic environment effect, faced with technical challenges such as huge spatial scale, multiple radiation sources and complex geomorphology.This paper demonstrates a massively parallel 3 Dfull-wave electromagnetic field simulation program software JEMS-FDTD and its development progress.Based on the large-scale parallel FDTD method and adaptive mesh refinement technology, an efficient full-wave time-domain computation method for large-area electromagnetic pulse propagation is developed.The numerical simulation of regional propagation of electromagnetic pulse in 1000 square kilometers city has been realized.This work verifies the feasibility of the software in the simulation of electromagnetic pluses regional propagation.
-
Key words:
- FDTD /
- EMP /
- regional propagation /
- adaptive mesh refinement /
- numerical simulation
-
具有更宽的工作带宽和更广的扫描范围的天线,增强了雷达的抗干扰能力、杂波抑制能力,同时保证了雷达在大范围内进行探测、跟踪和精确定位目标。然而,阵列单元之间的相互耦合会影响阵列扫描性能甚至带来破坏性的扫描盲点,使设计变得困难。传统的超宽带天线阵为保持阵列的性能,往往采用一些方法来抑制或者补偿单元间耦合,例如在单元间加入电磁带隙结构[1],或者超材料[2],但这些方法往往带来一些负面影响,如增益下降、方向图畸变、增加设计难度等。
紧耦合偶极子阵列(TCDA) 最早由Munk于2003年提出[3]。与传统阵列不同,TCDA通过减少相邻单元之间的距离来增加耦合,使得带宽有所拓宽。常见紧耦合单元结构有末端重叠型[4-5]及“交指”型[6]。
为了实现宽带和宽角扫描,需要在阵列孔径上方放置一层或者多层电介质板作为宽角匹配层(WAIM)。文献[7]中阵列上方引入两层电介质板,可实现9∶1阻抗带宽,扫描角度达60°。但是必须通过外部180°混合器来抑制偶极子相邻馈线之间的共模谐振。文献[8]采用同轴平衡馈线,并加载两层电介质板,可实现4倍阻抗带宽,但扫描角仅50°。文献[9]中的阵列引入了一种新型的寄生覆层,可以与偶极子巴伦共基板印刷,能够实现E面70°扫描,然而18段阶梯阻抗变换器用于馈电,结构复杂,加工精度非常高。文献[10]提出一种球杯型紧耦合结构,并加载单层电介质板,但仅能实现X波段最大70°扫描。为了减少输入阻抗便于匹配,文献[11-14]采用双偶极子形式,Wilkinson功分器馈电,但会导致E平面大角度扫描时电阻损失过大,扫描角限制在45°。
基于Munk的概念,本文提出了一种新型的TCDA设计。摒弃传统电介质层,选用频选表面(FSS) 作为宽角匹配层;巴伦选用微带线到共面平行双线过渡结构实现宽带馈电,结构简单。仿真结果表明,该阵列天线可实现3倍频内E面80°、H面45°的扫描角,且有源驻波比小于3。
1. 天线单元设计
1.1 TCDA等效电路
该文所提出的TCDA采用偶极子双面印刷于介质板,单元边缘部分相互重叠以增加耦合电容,重叠面积为w×l。图 1给出理想馈源馈电的TCDA模型及等效电路,其中dE, dH分别表示E面、H面单元间距;电介质层均由传输线表示,其特性由每个相应层内的传播Floquet模式确定。Ldipole表示偶极子自感,由偶极子的长度及宽度决定;Ccouple表示单元间耦合电容,取决于重叠面积,两者共同作用于抵消由地平面引入的电抗分量,以进一步扩宽带宽[15],表达式如下
Ldipole =μ04πdE×ln(2dEle) (1) 式中:le=l/4,为平面偶极子的等效宽度,μ0为自由空间磁导率;而Cp的计算采用平行板电容器模型,有
Ccouple =Cp1Cp2/(Cp1+Cp2) (2) 式中::Cp1=Cp2=ε0εr1S4πkd,k=8.99×109,S=w×l/2,ε0为自由空间介电常数,εr1为辐射层介电常数,d为介质板厚度。
利用射频软件ADS建立相应的等效电路,通过算法优化,进而初步确定参数值。等效电路结构简单,优化效率非常高。由于实际馈电方式引入的自身阻抗及金属结构都会影响电场分布,因此电路仿真结果不是非常准确,但仍可作为初步设计的参考。
1.2 宽带巴伦设计
对于TCDA单元,当E面尺寸dE等同于H面尺寸dH时,天线的输入阻抗约为200 Ω。如上文所述,为了便于阻抗匹配,很多TCDA设计采用双偶极子结构,即减小E面尺寸使得dE=0.5dH,此时天线输入阻抗降至100 Ω;为了减少端口数量,采用Wilkinson功分器同时给两个偶极子馈电。但是隔离电阻的存在,严重限制了E面大角度扫描。TCDA中dE=dH时,宽带巴伦须实现50~200 Ω的阻抗变换。由于辐射层为偶极子,采用同轴线直接馈电会导致电流分布不均,方向图畸变,因此巴伦在实现宽带阻抗变换的同时还需实现平衡不平衡转换。
巴伦结构如图 2(b)所示,介质材料为Arlon AD1000 (tm) (介电常数为10.2)。微带线宽带渐变,用于实现阻抗变换;扇形贴片φ=60°可将微带线与地之间垂直电场转换为平行双线之间的水平电场;平行双线实现电流不均匀分布到均匀分布,给偶极子平衡馈电。优化得到最佳参数值见表 1,S参数仿真结果如图 2(c)所示。可以看出,2~10 GHz频带内,该结构反射系数S11<-15 dB,传输系数S21>-3 dB,说明这种微带线到平行双线的过渡结构在5倍频带内实现良好的馈电匹配及平衡不平衡转换。
表 1 优化参数值Table 1. Optimized parameter value(mm) dE dH w l W2 W3 W4 W5 l2 l3 l4 Wb 14.3 14.3 3.4 7.7 0.55 1.25 3 3 5 18 3.2 8 1.3 天线单元结构及参数
根据等效电路初步确定的参数在全波仿真软件HFSS中建立TCDA模型。模型中引入优化后的宽带巴伦,完整TCDA结构包含宽角匹配层、偶极子辐射层、渐变巴伦、电阻型FSS以及金属反射板,如图 2(a)所示。除巴伦外,其余介质材料均为Rogers RO4003 (介电常数为3.55)。
对引入实际馈电方式的模型进行参数微调优化,得到最佳参数值见表 1及图 2。其中辐射层dE=dH=14.3 mm,w=3.4 mm,l=7.7 mm。电阻型FSS介质层厚度为0.6 mm,贴片长度l1=13.7 mm,宽度w1=1.7 mm,方阻为10 Ω/sq。对电阻型FSS进行仿真,将辐射层与金属反射板作为参考面,仿真得到谐振频率为7.5 GHz,而工作频带为2~6 GHz,谐振频率在工作频带之外。电阻型FSS具有一定损耗,因此工作频带内可以吸收部分反射波,有源驻波比得到降低,以牺牲辐射效率为代价来换取带宽的提高。
单元间的互耦使得阵列单元输入阻抗会随着扫描角度发生变化。Magill和Wheeler提出在阵列口径上方放置单层或多层介质板,在阵列参考位置上反射最接近纯电纳,可减小总的阻抗变化,进而起到宽角匹配效果[16]。所设计的TCDA采用频选表面取代传统电介质板,图 3给出等效电路。FSS介质层可等效为L-C级联电路;自由空间由传输线表示,其特性由传播的Floquet模式确定。FSS贴片谐振频率正比于单元结构的周长,因此该结构工作频率低于其谐振频率,旨在为沿z轴传播的平面波提供并联电容。这以与电介质覆层相同的方式降低了介质的特征阻抗,并且增加的电容减轻了接地平面引入的电感分量。从等效传输线观点看,宽角匹配层相当于阻抗渐变线,将天线与真空中的阻抗进行阻抗变化,而这种频选表面比传统电介质层增加了多个设计自由度,使得阻抗随角度变化更加稳定,因此匹配效果更好。
对于紧耦合阵列而言,E面耦合程度高于H面,因此E面扫描特性好,而H面较差。所加载的宽角匹配层设计时主要针对H面,以实现主平面的宽角扫描。频选表面介质材料选定为Rogers RO4003 (介电常数为3.55)。对天线扫描特性影响较大的参数优化后数值为:FSS与辐射层的距离gfss=1 mm;FSS介质层厚度hfss=1 mm;贴片长度lfss=8.58 mm及宽度wfss=1.2 mm。
为了验证FSS的性能,以H面45°扫描为参考基准,对比仿真加载电介质层的TCDA。除宽角匹配层外两者其余结构参数完全一致,周期边界下仿真无限大阵列。电介质层变量为介电常数εd,介质板厚度hd。εd分别取1.7,2.2,3.55,不同εd下优化hd,得到最优组合。图 4给出最优组合下H面45°扫描时的仿真结果。可以看出,当采用单层电介质板时,为了实现良好的宽角匹配,介电常数越低,所需介质板越厚,天线整体剖面越高。采用高介电常数覆层时,介质板厚度虽然降低,但Floquet模式更容易耦合到表面波,高阶模式将占主导地位而导致失谐。因此采用FSS作为宽角匹配层,可以极大地降低天线剖面及成本,同时还能实现更大的扫描角度。
2. 仿真结果与分析
在上节分析基础上,确定天线单元最终参数见表 1,使用电磁仿真软件HFSS对天线特性进行仿真与分析。图 5给出周期边界下无限阵列仿真结果。当天线进行波束扫描时,其有源驻波系数会随着扫描角发生变化。E面可以实现2~6 GHz工作频段内高达80°扫描,且有源驻波系数小于3,高频处匹配较差,因此驻波较低频处增大。H面有源驻波系数变化形式一致,随扫描角的增大而增大,可以实现3倍频内45°扫描。
图 6给出无限大阵列仿真得到的法向可实现增益,天线口径面积A一定时,增益随频率升高而增大。由于阵列电阻型FSS介质损耗及馈电部分损耗,低频处共面极化增益与理想最大增益差距略大,中心频率及高频处一致性较好。同时该阵列轴向交叉极化隔离度大于25 dB。
紧耦合天线阵由大量紧密排列的偶极子天线构成,为了加快仿真速度和节省计算资源,建立如图 7所示无限×有限阵列。有限维由20个元素组成,而无限维采用周期性边界条件,以达到近似全阵列仿真效果。阵列上方采用理想匹配层的边界条件,可以有效吸收阵列辐射。分别取4,5,6 GHz三个频点下对于20×∞单元的行阵列与列阵列做波束扫描分析。波束激励幅度按照均匀直线阵设置,相位差即为
Δφ=2πDfsinθscan /c0 (3) 式中:c0为自由空间波速,D为单元间距,f为工作频率,θscan为扫描角度。
传统相控阵增益会随着扫描角的增大而降低,并且单元间互耦使得有源反射系数增加。当主瓣扫描至60°时,增益下降3 dB;扫描角超过60°时,增益急剧下降;扫描至80°时,增益下降10 dB,严重限制了超宽角扫描[17]。
为了与传统相控阵相比较,图 8给出所提出的紧耦合天线阵相扫时归一化实际增益方向图,各频点处相扫实际增益均与法向实际增益相归一化。E面60°扫描时,归一化实际增益下降在3 dB以内,而当扫描角达到80°时,主瓣峰值处增益下降6 dB以内。H面最大角度45°扫描时,增益下降在2 dB以内。对比传统相控阵,紧耦合相控阵在实现超宽角扫描时具有极大优势。
同时在扫描波束范围内均未出现盲区,共面极化与交叉极化相差约20 dB,保持了良好的交叉极化隔离。对于接近法向的扫描角,旁瓣低于-13 dB,但对于较大的扫描角,旁瓣略有增加。对比三个频率E面及H面相扫方向图,可以看出更高频率时,阵列尺寸与波长相比变大,波束宽度变窄,方向图指向性更好。而当频率低于4 GHz时,由于阵列电尺寸小,辐射特性差,大角度扫描时波束峰值偏移较严重,所以未给出低频阵列方向图。综上所述,在无限×有限阵列,天线辐射特性较好,扫描时指向性良好,可以适用于具有宽带宽角扫描特性相控阵系统中。
3. 结论
本文基于紧耦合原理,提出一种宽带宽角扫描紧耦合偶极子阵列天线单元。结合等效电路,设计参数在电路中直接优化,最后在全波模拟求解器HFSS中进行微调,提高优化效率。频选表面取代传统电介质板用以改善宽角扫描时单元输入阻抗的变换,大大降低天线剖面、质量及成本。仿真表明天线可实现3倍频内E面80°、H面45°的扫描角,且有源驻波系数小于3。并仿真20×∞阵列进行验证,相扫过程中,交叉极化始终低于-20 dB,旁瓣低于-13 dB。所提出的TCDA设计满足宽带宽角扫描特性,且结构简单紧凑,易于加工制作。
致谢: 感谢吴鸿超老师在天线设计方面给予的指导和帮助。 -
-
[1] 刘尚合. 武器装备的电磁环境效应及其发展趋势[J]. 装备学院学报, 2005, 16(1): 1-6. doi: 10.3783/j.issn.1673-0127.2005.01.001Liu Shanghe. Effect of electromagnetic environment to weaponry and its trend of development. Journal of the Academy of Equipment Command & Amp, 2005, 16(1): 1-6 doi: 10.3783/j.issn.1673-0127.2005.01.001 [2] 孙国至, 刘尚合, 陈京平, 等. 战场电磁环境效应对信息化战争的影响[J]. 军事运筹与系统工程, 2006, 20(3): 43-47. doi: 10.3969/j.issn.1672-8211.2006.03.010Sun Guozhi, Liu Shanghe, Chen Jingping, et al. The effect of battlefield electromagnetic environment on information warfare. Military Operations Research and Systems Engineering, 2006, 20(3): 43-47 doi: 10.3969/j.issn.1672-8211.2006.03.010 [3] 张青洪, 廖成, 盛楠, 等. 森林环境电波传播抛物方程模型的改进研究[J]. 物理学报, 2013, 62(20): 189-195. https://www.cnki.com.cn/Article/CJFDTOTAL-WLXB201320025.htmZhang Qinghong, Liao Cheng, Sheng Nan, et al. Improved study on parabolic equation model for radio wave propagation in forest. Acta Physica Sinica, 2013, 62(20): 000189-195 https://www.cnki.com.cn/Article/CJFDTOTAL-WLXB201320025.htm [4] 申荣. 小区域的电磁环境评价研究[D]. 西安: 西安电子科技大学, 2006.Shen Rong. Study on evaluation of electromagnetic environment in small area. Xi'an: Xidian University, 2006 [5] 刘艳梅. 基于射线追踪法的电磁环境仿真系统实现的关键技术研究[D]. 北京: 北京邮电大学, 2013.Liu Yanmei. Research on key technologies for realization of electromagnetic environment simulation system based on ray tracing method. Beijing: Beijing University of Posts and Telecommunications, 2013 [6] Li Hanyu, Zhou Haijing, Liu Yang, et al. Massively parallel FDTD program JEMS-FDTD and its applications in platform coupling simulation[C]//IEEE International Symposium on Electromagnetic Compatibility (EMC Europe). 2014: 229-233. [7] 李瀚宇, 周海京, 廖成. 时域全波电磁计算程序JEMS-FDTD在复杂电磁环境研究中的应用[J]. 强激光与粒子束, 2014, 26: 073213. doi: 10.11884/HPLPB201426.073213Li Hanyu, Zhou Haijing, Liao Cheng. Application of JEMS-FDTD in complicated electromagnetic environment study. High Power Laser and Particle Beams, 2014, 26: 073213 doi: 10.11884/HPLPB201426.073213 [8] 鲍献丰, 李瀚宇, 伍月千, 等. JEMS-FDTD软件在飞机HIRF仿真中的应用[J]. 强激光与粒子束, 2017, 29: 103204. doi: 10.11884/HPLPB201729.170175Bao Xianfeng, Li Hanyu, Wu Yueqian, Application of JEMS-FDTD in HIRF simulation on aircraft. High Power Laser and Particle Beams, 2017, 29: 103204 doi: 10.11884/HPLPB201729.170175 [9] 徐小文, 莫则尧, 刘青凯, 等. 自适应结构网格上扩散方程隐式时间积分算法及其应用[J]. 计算物理, 2012, 29(5): 684-692. doi: 10.3969/j.issn.1001-246X.2012.05.008Xu Xiaowen, Mo Zeyao, Liu Qingkai, et al. An implicit time-integration algorithm for diffusion equations with structured AMR and applications. Chinese Journal of Computational Physics, 2012, 29(5): 684-692 doi: 10.3969/j.issn.1001-246X.2012.05.008 [10] 莫则尧. 并行自适应结构网格应用支撑软件框架研制进展[C]//中国工程物理研究院科技年报. 2009: 61-62.Mo Zeyao. Development of jparallel adaptive structured mesh applications infrastructure//Annual Report of Science and Technology of China Academy of Engineering Physics, 2009: 61-62 -