本发明涉及集成电路技术领域，尤其涉及一种差分过孔分析方法。
背景技术：
在多层高速pcb或封装系统中，过孔通常被用来连接不同层的信号线或者连接器件与电源平面和地平面。当高速信号沿着过孔传输时，过孔上垂直方向的电流会在平行板之间激励起向外传播的电磁波，从而使噪声耦合到相邻的信号孔上。作为信号传输路径中最常见的不连续因素，过孔会造成模式转换、串扰和不匹配等si(signalintegrity，信号完整性)问题。同时，过孔还可能激励起平行板结构的谐振模式，通过边缘的辐射造成潜在的emi(electromagneticinterference，电磁干扰)问题或者导致电源分配网络中剧烈的电压波动。所以在对高速pcb或封装系统进行si和emi分析时，准确的过孔建模是至关重要的。
目前大部分过孔的建模工作均是针对传统的过孔结构，即孔柱位于圆形焊盘中心的结构。因此对传统的过孔结构，反焊盘区的电磁场分布可以近似为同轴tem(transverseelectromagnetic，横电磁波)模式。
但是随着pcb或封装结构数据速率的升高，越来越多的差分过孔结构被用来传输高速差分信号。例如，差分过孔结构的反焊盘不再是圆形，该结构的特点是多个过孔共用一个不规则形状的反焊盘。对于这种差分过孔结构，反焊盘区的tem模式不再存在解析解。因此，必须使用数值方法来计算反焊盘区的tem模式。
目前已有的实现方案有两种：一种是利用区域划分的三维/二维fem(finiteelementmethod，有限单元法)求解不规则平行板上的紧凑型过孔结构。该方法首先利用二维静态有限元求解器得到反焊盘区的tem模式，然后将其作为源分别作用在过孔区域和平行板区域。其中过孔区域内复杂的场分布用三维有限元求解，并同时利用二维有限元求解平行板区域的传播模式，但对于平行板区域采用二维有限元法分析时，需要对计算区域进行网格划分，因此该分析方法效率不够高。另一种方法是采用三维fem/边界积分法求解无限大平行板上的紧凑型过孔结构，但是该方法局限于无限大平行板的紧凑型过孔结构，无法适用于实际的有限大不规则形状的平行板结构。
技术实现要素：
本发明的目的在于提供一种差分过孔分析方法，以解决现有的差分过孔分析方法效率不够高的问题。
为达此目的，本发明采用以下技术方案：
一种差分过孔分析方法，应用于含有差分走线的多层高速pcb，包括以下步骤：
s1：对差分过孔结构进行区域分解，分解形成有平行板区域和若干个过孔区域，若干个所述过孔区域与平行板区域之间形成有若干个交界面，所述平行板还包括外部边界；
s2：利用三维有限单元法求解若干个所述过孔区域；
s3：利用边界积分方法求解所述平行板区域；
s4：所述过孔区域与平行板区域s参数级联。
可选的，所述步骤s3包括：
通过边界积分方法求解不同所述过孔区域的端口之间的耦合，该边界包括所有所述交界面和所述外部边界。
可选的，所述步骤s3包括：
所述平行板之间的间距属于电小尺寸，使得所述外部边界可以近似等于pmc边界。
可选的，所述步骤s4包括：
计算出每个所述过孔区域的s参数矩阵后，将所有所述过孔区域的s参数矩阵整合成一个总的s参数矩阵。
可选的，所述步骤s4包括：
所述平行板区域的s参数矩阵与总的s参数矩阵通过所述交界面的水平端口级联起来。
与现有技术相比，本发明实施例具有以下有益效果：
本实施例提供的差分过孔分析方法采用边界积分法提取不规则形状平行板上不同过孔区域之间的耦合特性时，不需要对整个二维计算区域进行网格划分，而只用考虑一维边界，能进一步提升紧凑型过孔结构的设计效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其它的附图。
图1为本发明实施例提供的差分过孔分析方法的流程图。
图2为本发明实施例提供的含有三个过孔区域的不规则形状平行板的结构示意图。
图3为本发明实施例提供的用三维fem求解过孔区域的俯视示意图。
图4为本发明实施例提供的用三维fem求解过孔区域的侧视示意图。
图5为本发明实施例提供的沿着过孔区域与平行板区域的交界面定义水平方向端口的示意图。
图6为本发明实施例提供的平行板区域与过孔区域的s参数级联的示意图。
图7为本发明实施例提供的将过孔区域和平行板区域分为源、接收器和耦合路径的示意图。
具体实施方式
为使得本发明的发明目的、特征、优点能够更加的明显和易懂，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，下面所描述的实施例仅仅是本发明一部分实施例，而非全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
下面结合附图并通过具体实施方式来进一步说明本发明的技术方案。
请参阅图1所示，本实施例提供了一种差分过孔分析方法，应用于含有差分走线的多层高速pcb，包括以下步骤：
s1：对差分过孔结构进行区域分解，分解形成有平行板区域和若干个过孔区域，若干个所述过孔区域与平行板区域之间形成有若干个交界面，所述平行板还包括外部边界；
s2：利用三维有限单元法求解若干个所述过孔区域；
s3：利用边界积分方法求解所述平行板区域；
s4：所述过孔区域与平行板区域s参数级联。
具体的，在步骤s1中，过孔区域与平行板区域形成有若干个交界面，平行板还包括其外部边界。
请参阅图2所示，为了更好的说明本发明的技术方案，本实施例在平行板区域中分解有三个过孔区域，分别表示为a、b和c，该三个过孔区域与平行板区域的交界面分别为γa、γb和γc，平行板的外部边界用γout表示。
本发明的技术方案关键在于，在每个交界面的内部是过孔区域，这些过孔区域可以用三维fem来求解。同时，交界面的外部是平行板区域，可以用边界积分方法来求解出不同的过孔区域间的耦合。
请参阅图3和图4所示，过孔区域中的反焊盘是矩形，反焊盘的两边各有一个短路过孔。
其中，位于中间的两个孔为信号孔，位于两侧的两个孔为短路孔。
w、l、a、b、d、r均为过孔区域的各项尺寸。
平行板之间介质材料的介电常数为εr，损耗角正切为tanδ。
更进一步的，在步骤s2中，在解过孔区域时，其磁场强度h需要满足如下方程：
这里的v是指平行板之间的介质和上下的过孔柱，而s＝st+sb+sh，st、sb和sh分别指过孔上端和下端的反焊盘面以及平行板与过孔区域的交界面。在本例中，st和sb可分别定义为两个tem端口。同时，sh被划分为许多线性的片段。由于平行板之间的间距h为电小尺寸，并且sh与过孔结构有一定距离，因此sh上每个片段的场分布为平行板模式，水平方向的磁场与垂直方向的电场近似为常数。因此，在方程(1)中，
n×h＝jzez(2)
ez＝ezez(3)
请参阅图5所示，对于一个水平方向的端口i，其端口电压和电流可定义为：
这里的wi是指端口的长度。假设在st、sb和sh上分别有nt，nb和nh个端口。采用全波电磁仿真软件来计算这个多模式网络的s参数，可写成如下形式：
这里的v指代过孔区域a、b或c。
更进一步的，在步骤s3中，不同过孔区域的端口之间的耦合可以通过边界积分方法来求解。这里的边界包括内部过孔区域与平行板区域的交界面γin＝γa+γb+γc，和平行板的外部边界γout。根据等效原理，内部过孔区域与平行板区域的交界面γin上存在等效电流与等效磁流。通常情况下，平行板之间的间距h在考虑的频带范围内属于电小尺寸的。因此，平行板的外部边界γout可以近似当做pmc边界。同样根据等效原理，平行板的外部边界γout上的等效磁流或垂直方向的电场也是平行板区域中发出散射波的源。
采用二维空间的格林函数：
则平行板区域中由边界面上的等效电流jz激励产生的电场为：
这里的积分线即内部过孔区域与平行板区域的交界面γin＝γa+γb+γc，该交界面上分布着等效电流和等效磁流。同时，由交界面γin和平行板的外部边界γout上的等效磁流mt＝ezez×n激励的电场可以写为：
这里的
将(8)(9)带入(10)，可推导出联系jz和ez积分方程为：
以上的积分方程可用矩量法来求解。首先将图3中内部过孔区域与平行板区域的交界面γin和平行板的外部边界γout分别划分为nh和ne-nh个片段，即边界和交界面上的总片段数为ne。然后将每个片段上定义的单位脉冲函数作为基函数对未知量jz和ez进行展开。
这里的单位脉冲函数πn(r)定义如下，
jn和en分别是第n个片段上jz和ez的展开系数。通过点匹配的方法可以推导出：
zej＝zme(15)
其中j是nh×1的矢量，e是ne×1的矢量。ze和zm均是mom(矩量法，methodofmoment)矩阵。ze是ne×nh的矩阵，其元素可由式(16)计算，
式(16)中的υ＝1.781072418，zm是ne×ne的矩阵，其元素可(17)由计算，
根据(4)(5)交界面γin的每个片段上电压电流的定义，可以将jz和ez的关系转变为交界面上电压和电流的关系，
vp＝zppip(18)
这里的vp和ip是nh×1的矢量。zpp是nh×nh的矩阵，即平行板阻抗矩阵，表征了平行板上不同过孔区域的端口之间的耦合。
将式(4)(5)带入(15)，可以得到，
其中，vout是(ne-nh)×1的矢量，表示平行板外部边界γout上由每个片段的切向电场ez计算所得的电压矢量。
则平行板阻抗矩阵zpp可由(22)计算得到，
将(22)右边矩阵运算的结果取前nh行可得平行板阻抗矩阵zpp。
更进一步的，在步骤s4中，得到每个过孔区域的s参数矩阵后，可以将所有过孔区域的s参数矩阵整合成一个总的s参数矩阵。整合后的过孔区域的总的s参数矩阵可写成如下形式：
其中，(23)中的矢量元素都是由(6)中对应的矢量元素整合而成的。例如，
对于平行板区域，将平行板阻抗矩阵zpp转换成散射矩阵spp后有如下关系，
bp＝sppap(25)
请参阅图6所示，可以将平行板区域的散射矩阵spp与(6)中过孔区域的总的s参数矩阵通过交界面的水平端口级联起来。
根据电磁场的连续性，图6中的ap、bp与ah、bh有如下关系：
ap＝bh(26)
bp＝ah(27)
由(23)可得，
bh＝shtat+shbab+shhah(28)
将(26)(27)(28)代入(23)可得，
ah＝spp(i-shhspp)-1(shtat+shbab)(29)
将(29)代入(23)，可以消掉ah、bh，从而得到一个新的矩阵，该矩阵表征了过孔区域中所有的顶部端口、底部端口之间的散射关系。
这里的每个子矩阵的表达式如下，
每个过孔区域的s参数矩阵表征了过孔区域内部的连接特性，即过孔区域的顶部端口、底部端口和水平端口之间的连接关系。而该连接特性不受到平行板尺寸和边界条件的影响。这意味着，如果平行板上有若干个相同的过孔结构，则只需要计算一个过孔区域的s参数矩阵即可。该过孔区域的s参数矩阵也不会因为过孔的布局方位的改变而改变。不同过孔区域之间的耦合可以通过平行板区域的阻抗矩阵来表征。
请参阅图7所示，源或接收器的内在特性可以由过孔区域的s参数矩阵表征，而平行板区域的阻抗矩阵表征了一个过孔区域到另一个过孔区域的耦合路径。因此，将平行板区域与过孔区域分开求解有助于过孔位置的优化和过孔之间噪声耦合的减小。
本实施例提供的差分过孔分析方法采用边界积分法提取不规则形状平行板上不同过孔区域之间的耦合特性时，不需要对整个二维计算区域进行网格划分，而只用考虑一维边界，能进一步提升紧凑型过孔结构的设计效率。
以上所述，以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。