叶志红 苟 丹 吴小林 周健健
(重庆邮电大学通信与信息工程学院 重庆 400065)
随着通信技术和微电子技术的快速发展,新的电磁频谱不断引入,空间电磁环境变得越来越复杂。空间电磁场作用电子设备,会通过设备中的传输线产生干扰信号,进而对端接电路造成干扰,影响电路的正常工作。因此,模拟和分析空间电磁场对传输线端接复杂电路的电磁耦合效应,可为电子设备的电磁防护提供理论和数据支撑。
针对传输线端接复杂电路的电磁耦合分析方法,国内外学者已开展了大量的研究。首先,最直接的求解方法是全波算法。文献[1]通过推导3维集总元件的FDTD迭代公式,实现了对微带低通滤波器与非线性混合电路的全波分析。但是,该方法需要对微带线结构的滤波器直接建模,而且针对不同的电路,则需重新推导FDTD迭代公式,通用性不强。为此,国内外学者研究了多种场路混合算法,以避免对传输线和复杂电路的直接建模。其中,最具代表性的方法有BLT方程[2-5]、FDTD-SPICE方法[6-10]和FDTD-TL算法[11-13],等等。
传统的BLT方程是一种频域方法,只适用于传输线终端为线性和时不变负载的情况。因此,文献[14,15]提出了时域 BLT方程,用于求解传输线端接非线性负载的瞬态响应。但是,该方法需要用到大量的卷积计算,而且需要求解非线性矩阵方程来获得时域反射系数,计算效率不高。FDTD-SPICE方法是一种经典的全波算法与电路分析方法相结合的场路混合时域算法,其通过FDTD方法[16-18]计算传输线的激励场,并作为附加电压源引入到SPICE软件,联合传输线端接电路构建SPICE等效电路模型,再利用SPICE软件仿真得到电路上的瞬态响应。但是,该算法计算效率不高,是因为其将空间电磁场与电路响应计算分开处理。
FDTD-TL算法是本文前期研究成果,同样是一种时域方法,其首先通过FDTD方法模拟传输线周围空间的电磁场分布,并在每个时间步上引入传输线方程作为等效分布源项。然后,采用FDTD方法的中心差分格式离散传输线方程,求解得到传输线及其端接负载上的瞬态响应。相较于其他算法,该方法的优势在于实现了空间电磁场辐射与传输线瞬态响应的同步计算。但是,目前该方法只能模拟传输线端接线性电路的电磁耦合问题。
因此,本文将基于FDTD-TL算法,结合Ngspice软件求解复杂电路瞬态响应的优势,研究一种高效的时域混合算法,快速模拟传输线端接复杂电路的电磁耦合,并实现空间电磁场辐射与复杂电路瞬态响应的协同计算。
图1 时域混合算法的实施步骤
时域混合算法用于传输线端接复杂电路电磁耦合计算的实施步骤(见图1)为:(1) 通过CAD软件绘制移除传输线之后的理想导电板物理模型,保存为STL格式文件;(2) 将该STL格式文件导入FDTD仿真程序,实现对理想导电板的网格剖分,并利用FDTD方法计算得到传输线周围空间的电磁场分布;(3) 根据传输线的单位长度分布参数,求取传输线的特性阻抗,并替换复杂电路加载到传输线端口;(4) 使用传输线方程,构建电磁波作用传输线的电磁耦合模型,并对传输线方程采用FDTD方法的中心差分格式进行离散,迭代求解得到特性阻抗端口上的电流响应;(5) 在每个时间步上,提取特性阻抗上的电流加载到复杂电路作为激励源,联同复杂电路构建网表文件,然后通过Ngspice软件读取网表文件并仿真得到电路各个元件上的瞬态响应。
下面将从复杂电路激励源提取与电路传导干扰的Ngspice仿真两方面详细介绍时域混合算法的具体实施过程。
理想导电板上多导体传输线端接复杂电路的物理模型,如图2所示。当传输线距导电板的高度小于空间电磁场的最小波长时,传输线的辐射效应可以忽略。此时,空间电磁场作用多导线的电磁耦合可以通过传输线方程描述为
图2 导电板上多导线端接复杂电路的物理模型
图3 传输线的FDTD网格划分
将计算得到的特性阻抗上的电流加载到复杂电路端口作为激励源,如图4所示。使用Ngspice软件对电路的传导干扰进行仿真分析,具体的实施步骤为:首先,将激励源与复杂电路联合建立网表文件并保存为cir格式;然后,在FDTD的每个时间步上采用计算得到的入射电流数据替换cir文件中的源项数据;最后,调用Ngspice软件读取cir文件,从而仿真得到电路各元件上的瞬态电压响应,实现了空间电磁场辐射与电路瞬态响应的协同计算。
图4 电路传导干扰的Ngspice仿真过程
图5 理想导电板上端接复杂集总电路的多导体传输线模型
图6 集总电路中电阻R r2的瞬态电压响应
表1 两种方法计算所需内存和时间的对比
图7 理想导电板上端接复杂非线性电路的多导体传输线模型
图8给出了时域混合算法和CST仿真软件,模拟得到的复杂电路电阻元件 Rl2上的电压响应对比曲线。可以看出,虽然两种方法的计算结果在幅值上存在一些偏差,但是曲线的变化趋势和峰值基本保持一致,验证了时域混合算法处理复杂非线性电路的正确性。
算例3 为了验证该时域混合算法,模拟非平面波作用传输线端接复杂电路电磁耦合的正确性,将算例1的传输线端接复杂集总电路模型放置于屏蔽腔内,如图9所示。屏蔽腔的尺寸为L = 60 cm,W = 20 cm, H = 50 cm,厚度为1 cm。腔体的顶部开有一条矩形缝隙,尺寸为10 cm × 2 cm。入射波同样为高斯脉冲,垂直照射屏蔽腔,波形与算例1的相同。
时域混合算法与CST仿真软件模拟得到的复杂电路中电阻元件 Rg2上的电压响应对比曲线,如图10所示。可以看出,时域混合算法与CST的计算结果能够保持良好的一致性,进一步证明了该时域混合算法在处理复杂电磁环境下传输线端接复杂电路电磁耦合的正确性。
图8 非线性电路电阻R l2的瞬态电压响应
图9 屏蔽腔内端接复杂集总电路的多导体传输线模型
图10 屏蔽腔内集总电路电阻R g2的瞬态电压响应
本文将FDTD方法与传输线方程和Ngspice软件结合起来,研究了一种高效的时域混合算法,能够快速模拟空间电磁场作用传输线端接复杂电路的电磁耦合问题,并实现空间电磁场辐射与复杂电路瞬态响应的协同计算。首先,将传输线端接复杂电路通过特性阻抗进行等效,结合FDTD方法和传输线方程,提取特性阻抗上的电流响应作为复杂电路的激励源。然后,在每个时间步上,将该激励源与复杂电路联合构建网表文件。最后,通过Ngspice软件调用网表文件仿真得到电路各元件的瞬态响应。数值仿真表明,该时域混合算法所得电路元件瞬态电压与CST仿真结果基本一致,且占用内存和计算时间少,验证了算法的正确性和高效性。