基于TDIEFDTD混合算法电磁环境效应分析.docVIP

基于TDIEFDTD混合算法电磁环境效应分析 　　摘 要： 电磁环境效应是在现代航空、航海、航天以及地面系统设计中必须考虑的因素。针对带有线缆结构的复杂系统的电磁环境效应问题，发展了一种时域积分方程（TDIE）和时域有限差分（FDTD）的混合方法，采用时域平面波（PWTD）加速的TDIE方法求解场问题，高阶FDTD（2，4）方法求解线缆响应。数值算例表明了方法的有效性。 　　关键词： 电磁环境效应； 线缆； 时域积分方程； 时域有限差分； 时域平面波 　　中图分类号： TN01?34 文献标识码： A 文章编号： 1004?373X（2013）21?0157?06 　　0 引 言 　　随着数字设备和集成电路的广泛应用，现代电子信息系统对雷电、高功率脉冲武器等产生的强电磁环境极为敏感。这种瞬态电磁场会通过天线、孔缝或外部电缆，耦合到信息设备机箱内部线缆网络从而对终端电路产生电磁干扰（EMI），甚至造成损伤。 　　在系统设计阶段就对其进行电磁环境效应分析无疑是一种有效手段。但对于典型的系统级分析会面临多尺度问题，电子信息系统的载体一般是电大或超电大尺寸，系统的线缆网络在一维上（传输方向）电大、另一维（横截面）电小。 　　为了解决极其复杂的电磁环境效应分析问题，20世纪70年代Baum提出了电磁拓扑理论[1]（Electromagnetic Topology，EMT），其基本思想是将复杂系统分割为相互联系的空间，各空间单独求解，再将结果进行综合，并在此基础上发展了BLT（Baum?Liu?Tesche）方程方法[2]，但EMT理论人为地对系统进行分割，分析误差较大。传统的数值计算方法，如矩量法[3]（Method of Moment，MoM）、有限元法（Finite Element Method，FEM） [4]、时域有限差分方法（Finite Difference Time Domain，FDTD）[5]、传输线矩阵法（Transmission Line Matrix，TLM）[6]、时域积分方程法[7]等等，以及其加速、改进算法等[8?12]，解决多尺度问题、线缆加载问题仍有较大难度。一种可行的解决方案是将多种方法进行混合，已有的代表性的成果包括：Ba?c?采用FFT加速的TDIE方法求解带线缆结构飞机的电磁环境效应问题[13]；Xie采用FDTD方法和SPICE等效电路模型求解线缆网络的响应[14?15]； Liu采用FDTD方法和MNA方法混合求解带线缆和非线性负载结构的系统响应[16?18]。另外，一些电磁场商业软件（如CST、HFSS等）也纷纷推出复杂系统的电磁兼容和电磁环境效应仿真解决方案。 　　本文发展了一种基于TDIE?FDTD的混合方法，采用PWTD加速的TDIE方法求解场问题，高阶FDTD（2，4）方法求解线缆响应，结合两种方法的优势对系统的电磁环境效应进行有效预测，数值算例表明了方法的有效性。 　　1 PWTD加速的TDIE方法求解场问题 　　基于Maxwell方程组，金属问题的时域电场积分方程表达式为： 　　[n×n×Eir，t=n×n×μ04πs1R?Jr，τ?tds-14πε0?sds-∞τ??Jr，tRdt r，r∈s] 　　式中：[Eir，t]表示入射场；[Jr，t]为表面电流；[n×n]表示取矢量的表面负切向分量；[?]作用于场[r]上；[?]作用于源[r]上；[s]表示电流源分布表面；[τ=t-Rc]表示时间延迟。 　　传统的求解方法是时间递推方法： 　　式中：[Ij]为第[j]时刻电流系数向量；[Vj]为第[j]时刻激励向量；[Zi]为阻抗矩阵。 　　MOT方法的主要计算量在于式右端的矩阵矢量乘积求和运算，它表示先前各时刻表面电流源对当前时刻各场点的贡献，其运算量为[O（N2s）]。由于上述求和运算需要进行[Nt]个时间步，经典MOT算法的总计算量为[O（NtN2s）]。实际的计算实体（飞机、舰船）基本都是电大尺寸问题，计算量和存储量都十分巨大，因此必须采用加速方法来进行。 　　时域平面波（PWTD）方法的机理在于瞬态远场的平面波展开，通过这种展开可以减少式右端的求和运算。对于空间中某点的散射场，其贡献来源可分为两部分：一部分由该场点所在组的近场组[NFP（α）]中的源所产生；另一部分由远场组[FFP（α）]中的源产生。因此，利用PWTD加速MOT算法求解电磁散射的迭代公式可表达如下： 　　[Z0Ij=Vj-α∈NFP（α）l=0j-1Zααj-lIαl- α′∈FFP（α）n∈α′（n）fm（r），LeJn（r，t）t=lMtΔt+iΔt] （3） 　　式中：[Le·]表示电场积分算子；[Mt]为基本子信号离散序列的点数，[i]为整数，且[0≤i