看一看有限元软件ANSYS在传输线特性阻抗计算

摘要： 具有复杂截面的传输线的特性阻抗计算1直是领域内的热门问题，虽然已有相干的 解析解法，但是当截面形状不规则时，数值方法无疑是更好的选择。本文介绍了如何采取 有限元方法（FEM）计算均匀无耗传输线的特性阻抗，侧重探讨了如何在大型通用有限元软 件ANSYS上实现传输线特性阻抗计算的数值模拟，并以1种TEM传输室的设计计算为例，利用 ANSYS软件计算各段结构的截面特性阻抗，指导其各部分结构的肯定，从而到达优化TEM传输 室各设计参数的目的。测试表明计算的结果与实传输线阻抗值相吻合，这类方法大大的方便 了特性阻抗的肯定并提高了特性阻抗计算的精度。关键词：特性阻抗；有限元；传输线；ANSYS Application of Finite Element Software ANSYS to Calculate CharacteristicImpedance 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 TEMas 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 indicatesthat the method have a good achievement for improving the design of TEM cellKeywords：characteristic impedance; finite element; transmiss ion line；ANSYS 1 引言有限元方法是当今数值计算领域利用最广泛、最成熟的1种计算方法,其优点是通用性强、 精度高。随着有限元理论的深入研究，各种通用软件也应运而生，ANSYS以 其强大的前后处理功能，可靠的计算精度成了业界的佼佼者，充分利用他可以使研究者 的注意力集中到工程设计中来买的房子没有房照拆迁怎么办。遗憾的是目前其在电磁领域的利用明显弱于力学、结构等学 科。本文结合广泛利用于电磁兼容性实验TEM传输室和GTEM传输室的设计计算，利用ANSYS在电磁领域的计算模块，计算传输线特性阻抗，优化并简化TEM传输室的设计，并据此在ANSYS2次 开发出相应的模块。特殊截面传输线特性阻抗的计算1直是热门问题，目前普遍采取的解决方法通常有2类：第1类是解析类的保角变换法，通过保角变换可以把复杂的边界转 换成比较简单的边界，然后通过比较简单的办法进1步求得解析解，但是对不规则的截面 形状找到合适的变换函数有时几近是不可能的，因此这类方法实际上只能解决1些特殊截面 问题；第2类是近似的数值解法，以有限元法、有限差分法、边界元法为代表。 2 有限元法计算特性阻抗的基本原理由传输线理论可知，无耗均匀传输线特性阻抗Z0是表征长线固有特性的参数，仅决定 于传输线的界面形状、尺寸、周围介质而与频率无关。

其中：L0为传输线单位长度的电感量，C0为单位长度的电容量。电磁波在真空中的传播速度(光速)C由式(2)肯定：

由式(3)可知，特性阻抗的肯定关键是其散布电容的肯定。用有限元法求解散布电容的思路：以数值方法肯定截面的电场散布，再根据区域内能量与场 强的关系式：

得到区域内的电场能，以下是这1思路的基本推导进程：对以TEM为主模传输的大部分传输线来说，其内部场分部的求解可以采取2维准静态 的方法处理，在2维平面场的情况下，区域内电场散布用Lap lace方程第1类边界条件描述 为：

由变分原理可知上述边值问题等价于在满足条件φ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建立模型对1个给定的传输线截面来说，1般可以分为2种情况：1是截面为封闭系统，同轴线结 构及其各种情势的变种，如微带线、GTEM小室等林地违建一定要强拆吗，如图1所示。第2种是两导体间的介质空 间具有开放的边界，如平行结构双线，如图2所示。两种情况建模时要分别对待。对第1 种情况只需在内外导体间加电位差，划分网格计算即可；而对第2种情况则要在开放的边界 外再设置无穷边界单元，然后计算区域内的电场能量。

此处以TEM小室内传输截面为例属于第1种情况如图3所示。

根据前面的分析可知，提取截面的电容只要用到ANSYS的静电场分析模块，顺次进行以下步 骤：生成几何图形→设置材料属性、添加2维静电场单元PLANE121→设置内导体及侧板局部细化 参数→进行网格划分进行上述操作后即得到如图4所示的计算模型。3.2加载边界条件及求解对用静电场模块提取电容的加载很简单，只需在相应导体上施加电压载荷即可，在这里将 相当于外导体的各边施加电压0 V，内导体各边施加电压1 V，这样即在内外导体间施加了1 V的电势差，以上工作完成后求解即可。图5为求解后模型的能量散布。

3.3截面电容的提取ANSYS强大的后处理可以使用户方便地调用任何结果参数，结果的获得可以通过命令或GUI 方式，在后处理程序中依照以下步骤：定义单元列表→提取单元总能量→在参数菜单里提取电容。结果界面如图6所示。需要注意的是,当所分析截面为如图2所示的情况时应首先用2维无穷单元(如ANSYS程序中的 单元110)对外围阴影表示的无穷区域进行网格划分，无穷区域的范围随不同模型的具体情况 而变。 4结论本例中，复杂截面传输的特性阻抗得到了比较理想的控制，用Tektronix公司的1502B时域 反射计测试设计阻抗的结果1样验证了方法的可行性。由上述分析可知，利用ANSYS计算传输线特性阻抗方便快捷，相比经验公式他能提供更准确 的数据，又可省去繁琐的编程进程。ANSY还提供了宏（MACRO）、参数化设计语言（APDL）、用户界面设计语言（UIDL）和用户 可编程特性（UPFs）几种2次开发工具，利用这些2次开发工具可以容易的开发出针对特性 阻抗计算的模块，在ANSYS中特性阻抗的计算方法是：在后处理中，定义能量单元表，并相 加各单元的能量求出总能量，然后按式（11）求出散布电容强拆胜诉后如何赔偿，进而求得特性阻抗。要用到的 基本命令为： 参考文献［1］金建铭.电磁场有限元方法［M］.西安：西安电子科技大学出版社 ，1998.［2］Rankin D F,Novac B M, Smith I R.A Two Dimensional Electro magnetic Electroplastic Modeling Technique［A］.Pulsed PowerConference, 2003［3］王秉中.计算电磁学［M］.北京:科学出版社，2002(end)资讯分类行业动态帮助文档展会专题报道5金人物商家文章