几种计算微带线特性阻抗方法.docVIP

几种计算微带线特性阻抗方法摘 要： 本文从不同的计算思路和计算方法出发，介绍了几种计算微带线特性阻抗的方法，并对几种方法的特点进行了比较，为解决此类问题提供了一些方法的选择。 关键词： 微带线 特性阻抗 差分法 有限元法 1.前言 微带集成电路具有重量轻、体积小、频带宽、可靠性高、省电、低成本和长寿命等优点。在现代电子设备中得到广泛应用。微带线作为微带集成电路的主要部分，其设计和特性参数的计算受到广泛关注［１］。在设计微带器件时经常遇到计算微带线的特性阻抗的问题，目前分析微带线特性阻抗的方法很多，比如差分法、有限元法、保角变换法、格林函数法，等等［３］，本文以差分法和有限元法为基础介绍几种计算微带线特性阻抗的方法。 2.基本原理 微带线上传播的电磁波可近似看成TEM波，于是它的特性阻抗就能用下面的公式计算： Z==（1） 式中C、L分别为微带线单位长度的电容和电感，v为波在微带线上的传播速度。如假定微带线上不存在介质时单位长度的电容为C，这时线的电感L将不会因为电介质的存在与否而改变。又因介质不存在时线上波的传播速度为光速v，而且 v=（2） 由（2）式可解出L为：L= （3） 将L值代入（1）式即可求出微带线的特性阻抗Z为：Z= （4） 由（4）式可见，求微带线特性阻抗的关键在于求出介质存在和不存在时，微带线上单位长度的电容C和C。［３］根据计算思路，求这些电容的思路有：由微带线上的电位分布求解；由微带线单位长度的总电量推导出；由储藏在微带线上电场内的能量而推导出，等等。根据计算方法，求这些电容可以用差分法，也可以用有限元法。以下我就介绍几种常用的方法。 3.几种计算方法 ①先用有限元法（FEM）求微带线切面上的电位分布，再根据电位分布求出单位长度的微带线上的电容，进而求出微带线的特性阻抗。 假设微带线介质基片的厚度为h，相对介电常数为ε，导带的宽度为W，厚度为t，微带线的电位分布满足拉普拉斯方程：+=0 u|=u（5） 其中L是场域s的边界。 将场域s剖分为许多三角形单元，考虑场域s中任意一个单元e，假设其三个节点分别为i（x，y），j（x，y），m（x，y），如（图一）所示；三个节点的电位分u别为u，u，u，则三角形单元e内任意一点的电位u可以表示为三个节点电位值的线性插值，即： u（x，y）=［N］［u］（6） 其中：［N］=［N，N，N］ （6）式中［N］是形状函数，T表示矩阵的转置，对于三角形单元，形状函数［N］为： N=（a+bx+cy） N=（a+bx+cy）（7） N=（a+bx+cy） 其中Δe为第e个三角形单元的面积，Δe=1 x y1 x y1 x y=（bc-bc）。 a=xy-xy，b=y-y，c=x-x; a=xy-xy，b=y-y，c=x-x; a=xy-xy，b=y-y，c=x-x; 对拉普拉斯方程进行伽辽金有限元分析：在（5）式两边同时乘以［N］，并对三角形单元区域Δe进行面积积分得： ［N］（+）dxdy=0（8） 对（8）式进行部分积分，并将（6）式代入得： （ ）+ ）［u］dxdy-［N］dL=0 （9） （9）式左边第二项积分是对三角形单元Δe的边界L进行的，n表示边界的外法向。 （9）式是对区域s 中任意一个单元Δe进行的，对于整个区域则应对所有的单元叠加得到： + ）［u］dxdy-［N］dL=0（10） 对于方程式（10）左边进行第二次积分时要注意，要分清是内单元还是边界上的单元，对于内部单元，由于相邻两个单元之间的平衡关系，各单元边界积分值叠加为零，所以只有在整个求解域的外边界L上，这一项才须计算，故（10）式可以写成： + ）［u］dxdy-［N］［u］CdL=0（11） 其中表示对所有边界上的单元求和。 （7）式代入（11）式，最后可得：［k］［u］=0（12） 其中［k］称为三角形单元的总系数矩阵，解方程（12）即可以求得所需的电位值，进而求得微带线单位长度电容，代入公式（4）即可求得微带线的特性阻抗。［２］ ②先用差分法求出单位长度微带线上的总电荷量Q，再根据C=求出微带线的单位长度电容，进而求出微带线的特性阻抗。 为求这个Q值，取一个如（图二）所示的任意的包围内导体的环路abcd。这个环路的每个边均通过网格节点的中心并平行于坐标轴，其中一个边上任意点p处的电通量法向分量D 应等于： D=εE=-ε（13） 而D 可用p处相邻的网格节点e、f处的电位u、u表示（如图三），于是可得： D=-ε=-ε（1