电磁波作用PCB上微带线的电磁干扰时域分析方法技术

本发明专利技术提供了一种电磁波作用PCB上微带线的电磁干扰时域分析方法，首先，采用传输线方程构建电磁波作用PCB上微带线的电磁耦合模型，传输线方程的建模精度取决于微带线的单位长度分布参数和等效分布源项。因此，根据微带线的结构尺寸，结合经验公式，计算得到微带线的单位长度分布参数。然后，采用时域有限差分(FDTD)方法结合亚网格技术对PCB的基片和接地板进行建模，计算得到微带线的激励场，并引入到传输线方程作为等效分布源项。建立好传输线方程之后，采用FDTD方法的中心差分格式进行离散，从而迭代求解得到微带线和端接电路上的瞬态响应。本发明专利技术避免了对微带线结构的直接建模，可以有效改善传统全波算法在模拟PCB上微带线电磁耦合时剖分所需网格量较大而导致计算效率低下的问题。

Time domain analysis method of electromagnetic interference of microstrip line on PCB by electromagnetic wave

【技术实现步骤摘要】
电磁波作用PCB上微带线的电磁干扰时域分析方法
本专利技术涉及PCB板上微带线的电磁干扰分析方法，提出高效的时域混合算法，适用于电磁波对微带线电磁干扰的快速模拟与分析。
技术介绍
随着4G、5G通信技术的快速发展，使得通信、电子等领域使用的电子设备集成度越来越高。电子设备中的电路将半导体器件和微波、毫米波芯片等集成到PCB板上，各器件和芯片之间通过微带线进行数据传输和通信。对于PCB板电路而言，能够承受的功率比较小，在空间强电磁辐射源的作用下，将通过微带线耦合产生电磁干扰信号，进而对PCB板上的元器件造成干扰甚至是破坏。为此，需要分析电磁波对PCB板上微带线的电磁干扰特性，进而提出有针对性的干扰抑制方法。考虑到PCB板电路结构复杂，直接进行微带线的电磁干扰测试比较困难。因此，需采用数值算法进行模拟，以寻求微带线的电磁耦合规律。目前，分析微带线的电磁干扰最常用的方法是全波算法，例如，时域有限差分(FDTD)方法、矩量法(MOM)，等等。其中，FDTD方法是一种应用广泛的时域算法，其非常适用于宽频带信号在空间传播、辐射和耦合的数值模拟。采用FDTD方法对微带线进行网格剖分，然后按照时间步进的方式即可迭代求解得到电磁波在微带线上耦合产生的瞬态响应。但是，FDTD方法用于微带线电磁耦合分析的劣势在于需要对微带线直接进行网格剖分，由于微带线的结构比较精细，必然会造成网格需剖分得很细而导致网格量增加，造成计算资源的浪费。因此，迫切需要研究一种高效的时域混合算法，能够在避免对微带线直接建模的前提下，快速计算得到电磁波在微带线和端接电路上耦合产生的瞬态响应，并能保证与全波算法可比拟的计算精度。
技术实现思路
本专利技术所要解决的技术问题，就是针对现有技术处理微带线电磁干扰时耗费大量计算资源的缺点，提供一种高效的时域混合算法，能够快速地模拟和分析电磁波对微带线的电磁干扰问题。本专利技术解决所述技术问题，采用的技术方案是，电磁波作用PCB上微带线的电磁干扰时域分析方法，包括以下步骤：采用传输线方程，构建电磁波作用PCB上微带线的电磁耦合模型；结合经验公式，根据微带线的结构参数，提取微带线的单位长度分布参数；对移除了微带线的PCB板结构进行建模，计算微带线周围空间的电磁场分布，作为等效分布源项引入到传输线方程；采用FDTD方法对传输线方程进行差分离散，迭代求解得到微带线和端接负载上的电压和电流瞬态响应。本专利技术将FDTD方法的时域全波模拟的优势与传输线方程高效建立场线耦合模型的特点结合起来，并引入相应的网格技术，形成一种高效的数值计算方法，该方法首先根据传输线方程的思想，建立电磁波对微带线的电磁耦合模型。通过经验公式，根据微带线的结构参数，提取微带线的单位长度分布电感和电容参数，然后采用FDTD方法结合亚网格技术对移除了微带线的PCB板结构进行建模，仿真得到微带线周围的切向电场分量和垂直电场分量，从而引入到传输线方程作为等效分布电压源和电流源项。最后，采用FDTD的中心差分格式对传输线方程进行离散，迭代求解得到微带线和端接负载上的瞬态电压和电流响应。本专利技术的有益效果是，采用传输线方程建立电磁波对微带线的电磁耦合模型，有效地避免了对微带线精细结构的直接建模，并且应用亚网格技术对PCB板电路进行建模，提高了FDTD方法模拟微带线激励场的计算效率。附图说明为了更清楚地说明本专利技术实施例，下面将对实施例中所需要使用的附图和表格作简单地介绍。图1为本专利技术的流程图；图2为PCB板上两根微带线的结构示意图；图3为微带线激励场的计算示意图；图4为时空分离亚网格示意图；图5为空间亚网格界面磁场Hf插值示意图；图6为空间亚网格界面电场Ef的修正示意图；图7为时域混合算法与FDTD方法的仿真结果对比图。具体实施方式下面结合附图和实施例，详细描述本专利技术的技术方案。本实施例以电磁波作用两根微带线的电磁干扰为例进行说明。电磁波作用两根微带线的电磁耦合模型如图2所示，包括接地板1、基片2、微带线3，4、负载5，6，7，8、电磁波9。接地板1的尺寸为Lc×Wc。基片2的厚度为h，选用材料的相对介电常数为εr。微带线3和微带线4为矩形金属条，金属条长度为l，横截面宽度为w，厚度为t。负载5、负载6、负载7和负载8为电阻，电阻的阻值可自行定义。电磁波9为外部电磁干扰源，可以任意角度和极化方向照射微带线。Lc、Wc、w、t、l和h的具体参数值可以自行设定。如图1所示，本专利技术的实现流程包括如下：步骤1，根据微带线的宽度w、厚度t、与间距d和基片的厚度h与材料的介电常数，采用经验公式计算微带线的自电感Lt、自电容Ct、互电感Lm和互电容Cm，进而获得微带线的单位长度分布电感和电容参数。具体的经验公式为：μ0、ε0、εr分别表示自由空间磁导率、自由空间介电常数和基片材料的相对介电常数。步骤2，如图3所示，微带线与接地板之间的间距一般小于空间入射波的最小波长，在计算微带线的电磁耦合时，可以忽略微带线的电磁辐射效应。因此，对移除微带线的PCB板进行建模，并应用亚网格进行剖分，即对PCB板采用细网格进行剖分，其他区域仍采用粗网格进行剖分。对于亚网格技术而言，保持粗细网格界面上场分量的连续性是保证其良好后时间稳定性的前提。因此，本专利技术应用时空分离亚网格技术，将粗细网格界面分解为时间和空间两层亚网格界面，如图4所示。对时间亚网格界面的亚网格区磁场Hf由粗网格磁场H通过Taylor级数展开求得，表示为其中A＝(Hn+1-Hn-1)/2和B＝Hn+1-2Hn+Hn-1，v分别取为1/3,2/3或1，n表示时间步。空间亚网格界面的细网格磁场h由亚网格区的粗网格磁场Hf按照如图5所示的线性插值方式得到。线性插值的公式为h＝αHf1/3+(1-α)Hf2/3，α分别取为1或2，Hf1和Hf2分别表示相邻两个亚网格区的粗网格磁场。在一次细网格时间步进之后，如图6所示，根据细网格电场e修正空间亚网格界面上的粗网格电场Ef，修正公式为：Ef＝e1/9+2e2/9+e3/3+2e4/9+e5/9。在一次粗网格时间步进之后，将亚网格界面的电场Ef赋值给同位置的E以修正时间亚网格界面上的粗网格电场。最后，采用FDTD方法进行模拟，获得微带线沿线的切向电场分量以及与接地板之间的垂直电场分量。步骤3，建立适用于微带线电磁耦合分析的传输线方程，其形式为：L和C分别为微带线的单位长度电感和电容分布参数，V(y,t)和I(y,t)分别为微带线上的电压和电流，VF(y,t)和IF(y,t)分别为等效分布电压源和电流源项，其公式为：ET(y,t)和EL(y,t)由空间电磁场计算得到，可表示为h表示线缆距离地面的高度，和分别为微带线沿线的入射电场分量和垂直于线缆的入射电场分量。(x,y,z,t)本文档来自技高网...

【技术保护点】
1.电磁波作用PCB上微带线的电磁干扰时域分析方法，其特征在于，包括以下步骤：/n采用传输线方程，构建电磁波作用PCB上微带线的电磁耦合模型；/n结合经验公式，根据微带线的结构参数，提取微带线的单位长度分布参数；/n对移除了微带线的PCB板结构进行建模，计算微带线周围空间的电磁场分布，作为等效分布源项引入到传输线方程；/n采用FDTD方法对传输线方程进行差分离散，迭代求解得到微带线和端接负载上的电压和电流瞬态响应。/n

【技术特征摘要】
1.电磁波作用PCB上微带线的电磁干扰时域分析方法，其特征在于，包括以下步骤：
采用传输线方程，构建电磁波作用PCB上微带线的电磁耦合模型；
结合经验公式，根据微带线的结构参数，提取微带线的单位长度分布参数；
对移除了微带线的PCB板结构进行建模，计算微带线周围空间的电磁场分布，作为等效分布源项引入到传输线方程；
采用FDTD方法对传输线方程进行差分离散，迭代求解得到微带线和端接负载上的电压和电流瞬态响应。


2.根据权利要求1所述电磁波作用PCB上微带线的电磁干扰时域分析方法，其特征在于：所述微带线的单位长度分布参数包括电感参数和电容参数。


3.根据权利要求1或2所述电磁波作用PCB上微带线的电磁干扰时域分析方法，其特征在于：所述经验公式为：












Lt表示自电感，Ct表示自电容，Lm表示互电感，Cm表示互电容，l表示微带线长度，w表示微带线横截面宽度，t表示微带线厚度，μ0、ε0、εr分别表示自由空间磁...

【专利技术属性】
技术研发人员：叶志红，吴小林，周健健，苟丹，石艳超，汝梦祖，
申请(专利权)人：重庆邮电大学，
类型：发明
国别省市：重庆;50