有耗介质层上多导体传输线的电磁耦合时域分析方法*

叶志红 张杰 周健健 苟丹

(重庆邮电大学通信与信息工程学院, 重庆400065)

(2019 年 8 月 9日收到; 2020 年 1 月 5日收到修改稿)

目前, 针对空间电磁场作用有耗介质层上传输线的电磁耦合, 仍缺乏有效的数值分析方法. 因此, 本文提出一种高效的时域混合算法, 很好地解决了有耗介质层上传输线电磁耦合建模难的问题. 首先, 对经典传输线方程进行改进, 推导了适用于有耗介质层上多导体传输线电磁耦合分析的修正传输线方程. 然后, 结合时域有限差分方法和相应插值技术, 求解修正传输线方程, 获得多导线及其端接负载上的电压和电流响应, 并实现空间电磁场辐射与多导线瞬态响应的同步计算. 最后, 通过相应计算实例的数值模拟, 与CST软件的仿真结果进行对比, 验证了时域混合算法的正确性和高效性.

1 引言

非金属导电材料具备质量轻和导电性好的特点, 广泛应用于工业和军事等领域作为设备制造材料, 例如飞机蒙皮使用的复合材料就是典型的非金属导电材料. 在飞机等复杂目标中, 大量的传输线铺设在非金属导电材料上, 实现不同设备之间的数据通信. 当这类目标处于复杂电磁环境中时, 空间强电磁干扰源将通过目标内的传输线耦合产生强电流信号, 该信号必然流入传输线端接电路, 对电路上的敏感元件造成干扰或破坏. 因此, 实现对非金属导电材料上传输线的电磁耦合建模, 是分析飞机等电大尺寸目标电磁环境效应的前提.

由传输线理论可知, 在空间电磁场的激励下,传输线与其邻近的非金属导电材料之间将形成电流回路. 为了便于分析, 将非金属导电材料等效为具有一定介电常数和电导率的有限尺寸大小和厚度的有耗介质层. 全波算法是模拟有耗介质层上传输线电磁耦合的最直接方法, 其中时域有限差分(FDTD)方法[1−4]是应用非常广泛的时域算法. 但是, 采用FDTD方法进行模拟, 需要对有耗介质层和传输线精细结构直接建模, 剖分所需网格量较大, 势必造成计算效率低. 因此, 国内外学者基于传输线方程理论, 开展了大量的数值算法研究, 实现传输线电磁耦合的高效计算. 这类算法的核心思想是在避免对传输线直接建模的前提下, 获得与全波算法具有相同精度的计算结果. 其中, 主流的算法有 3类: Beam-Liu-Tesche (BLT)方程[5−11]、FDTD-SPICE (simulation program with integrated circuit emphasis)算法[12−17]和 FDTD-TL(transmission line)算法[18−20]. BLT 方程是将传输线看作管道, 激励源和负载看作节点, 构建节点电压电流经管道传播的关系矩阵, 进而求解矩阵以获得负载上的电压和电流响应. 但是, BLT方程是一种频域方法, 当入射波为宽频带信号时, 计算效率不高. FDTD-SPICE算法是一种时域方法, 其首先使用传输线理论建立传输线的SPICE等效电路模型[21], 然后通过FDTD方法模拟传输线的激励场并引入 SPICE软件作为激励源, 最后, 采用SPICE软件仿真得到传输线端接负载上的电压和电流响应. 但是, 该算法在推导SPICE等效电路模型时需要用到大量的理论推导, 且算法中负载瞬态响应与传输线激励场需分开计算, 因此计算效率不高. FDTD-TL算法是本文的前期研究成果, 首先使用传输线方程构建空间电磁场作用传输线的电磁耦合模型, 然后通过FDTD方法模拟传输线周围空间的电磁场分布, 并在FDTD的每个时间步上引入到传输线方程作为等效分布源项, 最后采用FDTD的中心差分格式离散传输线方程[22], 迭代求解得到传输线和端接负载上的瞬态响应. 相较于其他算法, 该算法实现了空间电磁场辐射与传输线瞬态响应的同步计算. 然而, 这类算法针对的研究对象均是理想地和实际地面上的传输线, 而对于有耗介质层上传输线的研究还未开展, 其原因是现有的传输线方程不适用于有耗介质层上传输线的电磁耦合分析.

因此, 本文首先基于经典传输线方程, 推导了适用于有耗介质层上传输线电磁耦合分析的修正传输线方程. 然后, 结合FDTD-TL算法和插值技术, 提出了一种高效的时域混合算法, 实现有耗介质层上多导体传输线瞬态响应的快速计算. 最后,通过相应计算实例的数值模拟, 验证了时域混合算法的正确性和高效性.

2 时域混合算法理论

2.1 有耗介质层上传输线电磁耦合分析的修正传输线方程推导

时域混合算法的核心是建立适用于有耗介质层上传输线电磁耦合分析的传输线方程. 经典的传输线方程包含两个方程, 称之为第一和第二电报方程, 其均由Maxwell方程组推导得到. 下面以经典传输线方程为原型, 详细介绍修正传输线方程的推导过程.

建立第一电报方程时, 需在传输线与接地面零电势位置之间选取合适的闭合回路. 当接地面为金属地时, 零电势位于金属地表面. 接地面为实际地面时, 零电势为无穷远处. 适用于金属地和实际地面上传输线电磁耦合分析的传输线方程及其推导过程, 已在文献 [23]中给出. 然而, 有耗介质层是具有一定介电常数和电导率的有限厚度的损耗介质, 已有的传输线方程无法满足电磁波作用有耗介质层上传输线的电磁耦合建模需求.

图 1 闭合回路和闭合曲面的选取Fig. 1. Selections of closed loop and surface.

当空间电磁场进入有耗介质层之后, 随着深度的增加而不断衰减, 因此可以将有耗介质层下表面近似看成零电势位置. 按照图1中的闭合回路, 应用电磁感应定律并按照积分形式展开得到:

其中,d和h分别表示有耗介质层的厚度与传输线架设的高度;,和均为总场. 考虑到传输线为良导体, 导体表面切向电场为零, 而且电磁场到达有耗介质层下表面时衰减到近似为零, 即此时, 将 (1)式对y求导, 可得

y处传输线与有耗介质层下表面之间的电压可表示为

基于此, (4)式可以表示为

其中,

建立第二电报方程时, 需要在传输线上选取合适的闭合曲面, 如图1所示. 根据安培环路定律方程两边同时进行闭合面积分,并应用奥氏公式和旋度的散度恒等于零的结论可得

其中,E为总电场,H为总磁场,J为传输线上的电流密度. 将电场E分解为入射电场Einc和散射电场Esca的叠加, 可得

将(10)式代入(9)式, 并将两边对y求导, 可得

(11)式即为修正以后的第二电报方程.

将(5)和(11)式由频域转换为时域, 即可获得适用于有耗介质层上传输线电磁耦合分析的时域传输线方程. 该方程可扩展到多导体传输线的情况, 表示为

修正传输线方程的方程结构形式与经典传输线方程[18]保持一致, 两者最大的区别在于等效分布源项的计算. 修正传输线方程的项表示有耗介质层下表面与传输线位置之间的垂直入射电场分量的沿线积分,为传输线位置的切向入射电场分量与有耗介质层下表面的切向电场分量之差. 换言之, 修正传输线方程考虑了有耗介质层内部的电场对传输线电磁耦合的影响.

2.2 有耗介质层上多导线电磁耦合的FDTD求解

修正传输线方程建立了空间电磁场与有耗介质层上多导线之间的电磁耦合关系, 通过求解传输线方程, 即可获得多导线上的瞬态响应. 需要说明的是, 传输线方程中的传输线单位长度分布参数和等效分布源项均为未知量, 其计算精度决定了传输线方程的建模准确度. 因此, 在求解传输线方程之前, 需要准确计算有耗介质层上多导体传输线的单位长度分布参数和等效分布源项.

2.2.1 多导线单位长度分布参数的计算

多导线单位长度电感参数矩阵L可由经验公式和计算得到, 其中i和j分别表示第i根和j根传输线,和分别表示第i根传输线的自电感以及与第j根传输线的互电感;,和分别表示第i根和第j根传输线的高度以及两根传输线之间的距离, 如图2所示. 电容参数矩阵C由公式计算得到.

图2 多导体传输线的横截面几何结构Fig. 2. Cross section geometry of multi-conductor transmission lines.

2.2.2 多导体传输线等效分布源项的计算

多导体传输线等效分布源项由传输线周围的空间电磁场获得, 而空间电磁场由FDTD方法计算得到[24,25]. 由 (16)和 (17)式可以发现, 等效分布源项只与入射电场有关, 而与多导线的散射电场无关, 这是因为: 一方面, 导线的横截面较小, 散射场较弱; 另一方面, 根据镜像原理, 有耗介质层对入射波的反射会抵消导线的部分散射场. 因此, 采用FDTD方法模拟多导线周围空间电磁场时, 只需对有耗介质层进行网格剖分, 而无需对多导线直接建模. 但是, 多导线的高度和间距可为任意值,即多导线未必落在FDTD网格的棱边上. 因此, 多导线等效分布源项所需的电场分量需要采用插值技术由相邻棱边上的FDTD电场分量计算得到,如图3所示.

图3 多导线沿线和垂直电场分量的插值示意图Fig. 3. Interpolation schemes of the electric fields along and perpendicular to the multi-conductor transmission lines.

将计算得到的传输线激励场在FDTD的每个时间步进上引入到传输线方程作为等效分布源项.建立好传输线方程之后, 采用FDTD方法的中心差分格式进行离散, 获得传输线电压和电流的FDTD迭代求解公式, 从而求解得到多导线及其端接负载上的电压和电流响应. 具体的迭代求解公式可由文献[15]获得.

对于多导线端接负载上的电压V0和VN不满足中心差分格式, 需分别采用前向差分和后向差分进行离散. 假定多导线按照FDTD网格划分成N段, 如图4所示, 负载上的电压迭代公式表示为

图4 传输线的 FDTD网格划分Fig. 4. FDTD grid division of transmission lines.

3 数值仿真与分析

采用时域混合算法对有耗介质层上单导体传输线和多导体传输线的电磁耦合进行数值模拟, 并与商业电磁仿真软件CST的计算结果进行对比,来验证算法的正确性和高效性.

算例1有耗介质层上单导线的电磁耦合模型如图 5所示, 有耗介质层大小为 0.2 m × 0.4 m,厚度为 0.01 m, 相对介电常数为 10, 电导率为20 S/m. 单导线长度为 20 cm, 高度为 1.9 cm, 端接负载分别为50和100 W. 入射波为高斯脉冲垂直照射单导线, 幅度为 1000 V/m, 脉宽为 2 ns. 为了保证计算精度, 时域混合算法选用的网格大小为 5 mm. 在计算空间电磁场分布时, 选用各向异性介质完全匹配层(UPML)截断边界, 入射波距离多导线的高度为4个空间网格大小. 图6给出了时域混合算法与CST微波工作室计算得到的负载R2上的电压响应对比曲线. 可以看出, 两种方法的计算结果振荡周期保持一致, 且幅值吻合度非常高. 表1列出了两种算法计算所需内存和时间的对比, 可以看出, 时域混合算法相较于 CST, 节省了47%左右的计算时间, 是因为时域混合算法无需对单导线直接建模. 这里需要说明的是, CST软件虽然提供线缆工作室模拟传输线的电磁耦合, 但是只适用于接地面为金属体的情况.

图5 有耗介质层上单导线的电磁耦合模型Fig. 5. Coupling model of single transmission line on the lossy dielectric layer.

图6 负载 R2 上的电压响应Fig. 6. Voltages on the load R2 computed by the two methods.

表 1 两种方法计算算例1时所需内存和时间对比Table 1. Memories and computation time needed by the two methods for the first example.

算例2有耗介质层上多导体传输线的电磁耦合模型见图7, 有耗介质层的大小为0.4 m ×0.7 m, 厚度为 0.01 m, 相对介电常数为 10, 电导率为50 S/m. 5根导线平行放置在有耗介质层上,长度为 0.5 m, 高度为 1.1 cm, 间距为 4 mm, 半径为 1 mm. 始端负载R1—R5均为 50 W, 终端负载R6—R10均为100 W. 入射波类型和算法选用的网格大小与算例1的相同.

图7 有耗介质层上多导体传输线的电磁耦合模型Fig. 7. Coupling model of multi-conductor transmission lines on the lossy dielectric layer.

图8 入射波垂直照射下的多导线端接负载的电压响应(a)负载R1上的电压; (b)负载R7上的电压Fig. 8. Voltages on the terminal loads of multi-conductor transmission lines under the condition of ambient wave perpendicular to the multi-conductor transmission lines:(a) Voltages on R1; (b) voltages on R7.

首先, 入射波角度设置为q= 180°,f= 90°和a= 180°, 即垂直照射多导线. 采用时域混合算法与电磁仿真软件CST计算得到负载R1和R7上的电压响应对比曲线, 如图8所示. 可以看出, 两种算法的计算结果基本保持一致.

然后, 考虑入射波斜照射的情况, 将入射角度设置为q= 135°,f= 45°和a= 180°, 同样采用两种方法计算得到负载R1和R7上的电压响应对比曲线, 如图9所示. 可以看出, 在入射波斜照射的情况下, 两种算法的计算结果仍能保证很好的吻合度. 另外, 表2列出了两种算法计算所需内存和时间的对比. 相较于算例1, 时域混合算法比CST节省了更多的计算时间, 是因为导线数量增加, 导致CST剖分所需网格量增多.

图9 入射波斜照射下的多导线端接负载的电压响应(a)负载R1上的电压; (b)负载R7上的电压Fig. 9. Voltages on the terminal loads of multi-conductor transmission lines under the condition of ambient wave oblique to the multi-conductor transmission lines: (a) Voltages on R1; (b) voltages on R7.

表2 两种方法计算算例2时所需内存和时间对比Table 2. Memories and computation time needed by the two methods for the second example.

4 结论

经典传输线方程不适用于电磁波作用有耗介质层上传输线的电磁耦合问题. 因此, 基于经典传输线方程, 推导了适用于有耗介质层上多导体传输线电磁耦合分析的修正传输线方程. 然后, 结合前期研究的FDTD-TL算法和相应的插值技术, 提出了一种高效的时域混合算法, 实现修正传输线方程的快速求解, 获得多导体传输线及其端接负载上的瞬态响应. 该时域混合算法避免了对多导线精细结构的直接建模, 并实现了空间电磁场辐射与多导线瞬态响应的同步计算. 通过相应计算实例的数值模拟, 验证了时域混合算法能够与CST微波工作室的全波仿真保持相同的计算精度. 当使用相同网格大小时, CST软件所需内存和计算时间较时域混合算法有优势, 但由于在时域混合算法中精细的传输线结构无需剖分网格, 因此, 在获得相同精度的条件下, 时域混合算法相较于CST软件具有显著的优势.