时域电场积分方程在求解架空屏蔽电缆蒙皮电流中的应用.pdfVIP

维普资讯 第 17卷 第 4期 强 激 光 与 粒 子 束 Vo1．17，NO．4 2005年 4月 HIGH POW ER LASER AND PARTICLE BEAMS Apr．，2005 文章编号 ： 1001—4322(2005)04—0637—04 时域 电场积分方程在求解架空 屏蔽电缆蒙皮电流 中的应用 刘顺坤， 陈雨生， 孙蓓云，周 辉， 谢彦召 (西北核技术研究所，陕西 西安 710024) 摘 要： 时域电场积分方程是研究电磁辐射及散射等问题的重要方程。在数值计算中，该方程方法 只需 将散射体进行剖分 ，而不必将剖分推至整个计算域内进行 ，计算效率较高 。将时域电场积分方程方法引入到 电 磁脉冲作用下架高屏蔽电缆蒙皮电流 的计算 中，研究 了电缆 的蒙皮感应电流分布及波形特征。通过与辐射波 电磁脉冲模拟器的实验结果 的比对 ，证明了数值结果的可靠性 。 关键词 ： 时域积分方程 ； 有限差分； 电磁脉冲； 屏蔽电缆 中图分类号 ： O441 文献标识码 ： A 研究 电磁脉冲对架高状态屏蔽电缆的耦合特性的意义在于，此种状态是耦合最为严重的情况 。电磁脉冲 对架高 电缆的耦合特性 ，对于电缆抗 电磁脉冲耦合防护具有重要指导意义 。时域 电场积分方程是研究电磁辐 射及散射等问题的重要方程 ，尤其在处理细线结构时具有较高的计算效率 ，这缘于使用该方法只需将散射体进 行剖分，而不必将剖分推至整个计算域 内进行。相对于电磁脉冲最高频率的波长及 电缆的长度 ，电缆的线径要 小得多，在采用时域数值方法求解 电缆蒙皮电流时，可以将 电缆作为细线处理 ，采用时域 电场积分方程方法求 解 电缆上的蒙皮电流分布。 时域 电场积分方程的求解一般采用时域矩量法 [1]，通过展开函数 (基函数)及实验函数 (权函数)的构造 ，求 得散射体上电流的分布，进而求解散射体 (辐射体)的散射 (辐射)特性。由于电场时域积分方程的核具有高阶 的奇异性，因此 ，采用矩量法求解时其展开函数不但在时间上不能采用脉冲函数 ，而且在空间上也不应采用脉 冲函数。为使方程能应用到有急剧弯曲的结构，还要求采用具有非零导数的函数 。这使得时域 电场积分方程 的矩量法求解相当繁琐 ]。此外 ，矩量法求解时域 电场积分方程还涉及结构矩阵的求解 ，当剖分单元过多时， 算法的时序推进稳定性不能得到保证 。这些原因限制了时域 电场积分方程的应用。 为此 ，提出了一种有限差分求解时域 电场积分方程的方法[3]。首先推出了适合有 限差分求解的积分方程 形式 ，进而对该方程进行差分求解 。该求解方法不需引入展开函数 ，直接对未知电流的时间、空间微分进行差 分离散，进而时间推进求解。这样避免了时域矩量法引人展开函数及实验函数的繁琐，同时也避免了求解大型 结构矩阵。因为方程的核仅包含一阶奇异性，所 以数值稳定性相对较好 。 1 数值模拟方法的建立 [3] 由线电流 l(r，tI)、电荷密度p(r，tI)产生的散射场可以推出适合有限差分求解的细线时域 电场积分方程 的另外一种形式 ；· (r，)一 ；· 一 (；． J(r ds (1) ， 有限差分求解的时域 电场积分方程推导如下。 令F(r，￡)一；．F(r，￡)一f；．；，J(r，，￡，)／sds，，H一 ，G(r，￡)一fH(r，，￡)／sds，，得 47 一 3t 一 一 坚 + C2 (2) “ 8t‘ d5 式中：r，t表示场项，r，，￡表示源项 ，且满足 tI一￡一Ir—r，I