瞬态电磁场对双层圆柱腔体耦合效应的三维FDTD数值模拟.pdfVIP

瞬态电磁场对双层圆柱腔体耦合效应的三维FDTD数值模拟 孟萃1’2陈雨生2程建平1刘以农1 (1．清华大学工程物理系，北京，100084’2． 西北核技术研究所，西安，710024) 摘要本文利用时域有限差分方法(FDTD)，对置于屏蔽体内圆柱筒的瞬态电磁脉冲 的耦合效应进行了数值模拟。针对不同半宽和上升前沿的入射波的耦合效应做了比 较，给出了屏蔽体内场分布的基本规律以及屏蔽体外壁上的皮电流的分布规律，并对 此开有小天线窗口屏蔽体的屏蔽性能做了评估。 关键字瞬态电磁场 孔洞 耦合 时域有限差分法 1前言 无论电磁脉冲或高功率微波都与射频电磁场一样是通过一定的耦合途径进入电子系统并产生 影响的。耦合途径可分为辐射耦合和传导耦合。而空间电磁波直接通过壳体的孔洞、缝隙、屏蔽 体泄露耦合至天线就是辐射耦合的一种。电子设备处于金属屏蔽体内能够很好地避开电磁干扰， 但屏蔽体上不可避免地要开有天线窗、散热窗或留有电缆通道，电磁场通过这些孔洞耦合入屏蔽 体的场会导致电气设备发生状态翻转、性能失效甚至彻底的毁坏。而峰值高的瞬态场由于频率成 分丰富，造成损害的可能性就更大，从而加大了抗辐射的难度。典型的瞬态场如高空核电磁脉冲 因其幅值高(几万伏／米)，频带宽对系统的危害极为严重。因此研究瞬态电磁场对各种目标的耦 合规律是必要且重要的。 由于目标的尺寸大小、形状不是规则的，要求耦合场的解析解几乎是不可能的，FDTD由于它 的简单、直观等特性得到广泛应用。时域有限差分法直接把含时间变量的l妊axwell弦度方程在Yee 氏网格空间中转换为差分方程。在这种差分格式中每个网格点上的电场(或磁场)分量仅与它相 邻的磁场(或电场)分量及上一时间步该点的场值有关。时域有限差分法把各类问题都作为初值 问题来处理，使电磁波的时域特性被直接反映出来。 2物理模型 选取瞬态场波形是平面双指数波，其波阻抗为平面波阻抗r／=377Q，表达形式为： E=E。(P一一e一声)，调整a和13的值即可以得到不同上升前沿和不同半宽度的入射波。在计算 中，入射波从连接边界引入，连接边界内的电磁场是入射电磁场与目标散射电磁场的总和，而吸 收边界与连接边界间的场只有散射场。计算的物理模型如图1和图2所示。 图l计算物理模型 图2散射体物理模型 3数学方法 整个物理过程满足无源的麦克斯韦方程组： 726 v×亏：一旦否一厅 (1) v×厅：旦西+．7 (2) V·E=p／Fo (3) ， v．百=0 (4) 用时域有限差分法【2】，采用Yee氏元胞将(1)和(2)式进行差分离散，得到时域差分方程如下： ．盯(f+{，，，k)at 碟Ⅵ砖舭卜粪雠”扣∽+志‘疆粤1· 2￡(f+毒，／，k) 2s(i+专，J，七) 日：n+{(f+号，_，+号，七)一日，1(“告，，一{，七)ⅣyI1+兰“+{，_，，七+号)一Ⅳ≯(f+{，-，，七一号) 弩 ＆ 、 其它方程依此类推。 在两种媒质交界面上的电流密度满足：歹：．i；×(疗，一厅．)，由于理想导体内部场量为零，所 以导体壁上的电流密度为歹；疗。／7，。本计算是在宣角坐标系下进行的，而目标的形状为圆柱形， 因此某时刻导体壁上的电流的数学表述如下：