有限元软件ANSYS在传输线特性阻抗计算中的应用

摘要: 具有复杂截面的传输线的特性阻抗计算一直是领域内的热点问题,尽管已有相关的 解析解法,但是当截面形状不规则时,数值方法无疑是更好的选择。本文介绍了如何采用 有限元方法(FEM)计算均匀无耗传输线的特性阻抗,着重探讨了如何在大型通用有限元软 件ANSYS上实现传输线特性阻抗计算的数值模拟,并以一种TEM传输室的设计计算为例,应用 ANSYS软件计算各段结构的截面特性阻抗,指导其各部分结构的确定,从而达到优化TEM传输 室各设计参数的目的。测试表明计算的结果与实传输线阻抗值相吻合,这种方法大大的方便 了特性阻抗的确定并提高了特性阻抗计算的精度。
关键词:特性阻抗;有限元;传输线;ANSYS
Application of Finite Element Software ANSYS to Calculate Characteristic
Impedance of Transmission Line
LIU Bin,LI Gang
(EMC Lab,Department of Mechanical Engineering, Southeast Unive rsity, Nanjing,210096, China)
Abstract: The calculation of characteristic impedance of transmission line with complicate d section is a hot spot in field Although some analytic method is available, numeric method is a better choice when the section is very irregular. This article presents how to calculate characteristic impedance of transmissi on line based on the Finite Element Method (FEM) The numerical simulation of c haracteristic impedance of transmission line, which is implemented by using ANSY S FEM software, is focused especially here Taking the design of a sort of TEM as example,The ANSYS FEM software is used to determine the dimension of se ction of the cell, in order to optimize the design parameters Subsequent testi ng inosculates with the result from ANSYS simulation.The application indicates that the method have a good achievement for improving the design of TEM cell
Keywords:characteristic impedance; finite element; transmiss ion line;ANSYS
1 引言
   有限元方法是当今数值计算领域应用最广泛、最成熟的一种计算方法,其优点是通用性强、 精度高。随着有限元理论的深入研究,各种通用软件也应运而生,ANSYS以 其强大的前后处理功能,可靠的计算精度成为了业界的佼佼者,充分应用他可以使研究者 的注意力集中到工程设计中来。遗憾的是目前其在电磁领域的应用明显弱于力学、结构等学 科。
   本文结合广泛应用于电磁兼容性试验TEM传输室和GTEM传输室的设计计算,应用ANSYS在电磁领域的计算模块,计算传输线特性阻抗,优化并简化TEM传输室的设计,并据此在ANSYS二次 开发出相应的模块。特殊截面传输线特性阻抗的计算一直是热点问题,目前普遍采用的解决方法通常有2类:第一类是解析类的保角变换法,通过保角变换可以把复杂的边界转 换成比较简单的边界,然后通过比较简单的办法进一步求得解析解,但是对于不规则的截面 形状找到合适的变换函数有时几乎是不可能的,因此这种方法实际上只能解决一些特殊截面 问题;第二类是近似的数值解法,以有限元法、有限差分法、边界元法为代表。
2 有限元法计算特性阻抗的基本原理
    由传输线理论可知,无耗均匀传输线特性阻抗Z0是表征长线固有特性的参数,仅决定 于传输线的界面形状、尺寸、周围介质而与频率无关。
其中:L0为传输线单位长度的电感量,C0为单位长度的电容量。
电磁波在真空中的传播速度(光速)C由式(2)确定:
由式(3)可知,特性阻抗的确定关键是其分布电容的确定。
用有限元法求解分布电容的思路:以数值方法确定截面的电场分布,再根据区域内能量与场 强的关系式:
得到区域内的电场能,以下是这一思路的基本推导过程:
    对于以TEM为主模传输的大部分传输线来说,其内部场分部的求解可以采用二维准静态 的方法处理,在二维平面场的情况下,区域内电场分布用Lap lace方程第一类边界条件描述 为:
    由变分原理可知上述边值问题等价于在满足条件φc=f1的函数类中 ,求φ以使下面的泛函取得最小值:
   不难得出,变分原理在此处的物理意义为:当用边值问题求解静电场电位时,在求解区 域中当电磁场稳定后电场能量达到最小。
进行场域的剖分、插值,对式(6)进行取极值得到刚度矩阵:
式(8)即为有限元法的代数方程组,[k]为其系数矩阵,[k]=(kij) (i,j=1,2, 3,…,n0)。
对方程组的求解即可得到场域位函数Φ(x,y),根据式(9)容易得到区域内的电场解:
至此电场分量的求解完成。将式(4)应用到这里就可得到区域内的电场能量,如式(10 ):
其中:Ex,Ey是单元I上电场的对应分量,Si为对应场域的面积。
求得电场能量后,再由式(11)即可得出截面分布电容:
其中:ΔV为两导体间电势差。
3 分析过程
3.1建立模型
   对于一个给定的传输线截面来说,一般可以分为2种情况:一是截面为封闭系统,同轴线结 构及其各种形式的变种,如微带线、GTEM小室等,如图1所示。第二种是两导体间的介质空 间具有开放的边界,如平行结构双线,如图2所示。两种情况建模时要分别对待。对第一 种情况只需在内外导体间加电位差,划分网格计算即可;而对第二种情况则要在开放的边界 外再设置无限边界单元,然后计算区域内的电场能量。
    此处以TEM小室内传输截面为例属于第一种情况如图3所示。
    根据前面的分析可知,提取截面的电容只要用到ANSYS的静电场分析模块,依次进行如下步 骤:
    生成几何图形→设置材料属性、添加二维静电场单元PLANE121→设置内导体及侧板局部细化 参数→进行网格划分进行上述操作后即得到如图4所示的计算模型。
3.2加载边界条件及求解
   对于用静电场模块提取电容的加载很简单,只需在相应导体上施加电压载荷即可,在这里将 相当于外导体的各边施加电压0 V,内导体各边施加电压1 V,这样即在内外导体间施加了1 V的电势差,以上工作完成后求解即可。图5为求解后模型的能量分布。
3.3截面电容的提取
   ANSYS强大的后处理可以使用户方便地调用任何结果参数,结果的获得可以通过命令或者GUI 方式,在后处理程序中按照以下步骤:
    定义单元列表→提取单元总能量→在参数菜单里提取电容。结果界面如图6所示。
    需要注意的是,当所分析截面为如图2所示的情况时应首先用二维无限单元(如ANSYS程序中的 单元110)对外围阴影表示的无限区域进行网格划分,无限区域的范围随不同模型的具体情况 而变。
4结论
    本例中,复杂截面传输的特性阻抗得到了比较理想的控制,用Tektronix公司的1502B时域 反射计测试设计阻抗的结果同样验证了方法的可行性。
   由上述分析可知,利用ANSYS计算传输线特性阻抗方便快捷,相比经验公式他能提供更准确 的数据,又可省去繁琐的编程过程。
   ANSY还提供了宏(MACRO)、参数化设计语言(APDL)、用户界面设计语言(UIDL)和用户 可编程特性(UPFs)几种二次开发工具,利用这些二次开发工具可以容易的开发出针对特性 阻抗计算的模块,在ANSYS中特性阻抗的计算方法是:在后处理中,定义能量单元表,并相 加各单元的能量求出总能量,然后按式(11)求出分布电容,进而求得特性阻抗。

 

返回Ansys文章专题列表>>>