modeling and - 程序员联合开发网

1、超宽带天线的同轴线馈电结构建模与模拟研究摘要 同轴线馈电是超宽带天线常用的馈电方式。文中运用时域有限差分法（FDTD），对同轴馈电的单极子天线的电磁特性进行时域仿真，基本思路是把同轴馈电的天线计算分为激励网格和天线网格，并通过反射场/总场的分离边界把入射场加入到天线网格中。首先对算法的原理介绍及其在仿真计算中的具体实现作介绍，同时给出了阶梯近似对同轴线馈电和单极振子天线建模的具体算法，然后实例分析了同轴馈电单极子天线的若干性能，对比粗网格模拟和细网格模拟单极振子天线所产生的不同。用FDTD法对天线的仿真分析，可用作天线设计的一种快速、经济的辅助手段。关键词 时域有限差分法 阶梯近似网格剖分 二

2、阶MUR吸收边界条件 单极子天线 MODELING AND SIMULATION OF COAXIAL FEED LINE STRUCTURES FOR ULTRA-WIDE BAND ANTENNAAbstract:Coaxial feed structures are widely used in ultra-wide band antennas . This paper modeled the characteristic of the monopole antenna feeded by coaxial line by FDTD in the time-domiain,which sh

3、owes that . Firstly, it introduced the theory of the arithmetic and the particularly realization in the calculation; then it described the use in the time-domain; finally it analysed several characteristics of the monopole antenna. The arithmetic used in the microstrip antenna is also a quick and ec

4、onomical way to design the antenna.Key words:FDTD;2nd order MUR absorbing-boundary; staircasing technique; monopole antenna1.引言 超宽带天线是针对天线收发信号的相对带宽而言的：当信号的带宽与中心频率之比小于1%时称为窄带，带宽与中心频率之比在1%与25%之间称为宽带，带宽与中心频率之比大于25%时称为超宽带。超宽带天线就是用来发射和接收超宽带信号的天线装置。而天线的辐射场要靠源来激发，如何设置符合实际的激励源，是计算天线辐射特性的关键之一。源的设置方式要与天线的馈电方式

5、相一致，以保证天线的辐射，图形与实际一致。不同天线有不同的馈电方式，相应的激励源的设置也有所不同。同轴线是以TEM模为主模的传输线，可以传输超宽频带信号，由于其本身的一些优点，同轴线馈电是超宽带天线和高功率微波常用的馈电方式。 1966年KSYee创立了计算电磁场的时域有限差分法(FDTD)1，它是一种计算时变电磁场有效的数值方法。它把含时间变量的Maxwell旋度方程转化为差分形式，并在时间轴上依次推进的求解空间电磁场的数值计算方法。其已在计算电磁场的各个领域得到了广泛应用，尤其是在天线分析。作为一种电磁场的数值计算方法，时域有限差分法在计算天线的特性上具有一些很突出的优点：作为一种瞬态方法

6、，在脉冲波的激励下，FDTD方法的一次计算结果经Fourier变换后可获得丰富的频域信息；适合模拟各种复杂电磁结构，用FDTD的离散空间网格可以较精确的模拟天线(阵)的实际结构；易于得到计算空间场的暂态分布情况，这既便于定性理解其工作的物理过程，又便于得到供定量分析的有关电参量；它所需要的计算机内存和CUP时间与网格单元成正比，并且不需要矩阵求逆，明显优于传统的矩量法。本文首先从FDTD几个关键步骤入手，从理论上加以分析说明。并用这些措施对单极子天线进行分析计算。给出电磁计算的仿真结果。因此FDTD为天线的设计提供一种快速、经济的辅助手段。下面我们将用时域有限差分法对超宽带天线同轴馈电结构实例

7、进行研究，即对同轴馈电的单极子天线进行建模与模拟研究,天线结构如图1所示。图 1 金属平板上圆柱天线模型2.天线建模与仿真 本文对图1中天线进行建模和模拟研究的基本思路是把同轴馈电的天线计算分为激励网格和天线网格，并通过反射场/总场的分离边界把入射场加入到天线网格中。这种思想将激励设置划分出来成为一个单独的网格空间(激励空间)，在此网格空间只存在与天线馈线相同的传输线，无任何其它结构，而所研究的天线结构处于另一个网格空间之内(天线结构空间)，见图2。激励空间的作用是迭代产生天线馈线入射波场，为使激励空间仅有入射波，传输线两端用吸收边界端接，在本文中吸收边界采二阶MUR吸收边界。 下面首先介绍时

8、域有限差分基本原理，然后讨论激励网格和天线网格的建模，同时给出同轴线的激励源设置方法，最后实例分析了同轴馈电单极子天线的若干性能，对比粗网格模拟和细网格模拟单极振子天线所产生的不同点。图2.1 天线网格与入射场网格的连接示意图2.1时域有限差分法基本原理麦克斯韦方程及其差分格式麦克斯韦方程的积分方程 （2-1-1） （2-1-2） （2-1-3） （2-1-4）它们相应的微分方程为 （2-1-5） （2-1-6） （2-1-7） （2-1-8）其中E为电场强度，单位为伏特米()；D为电通量密度，单位为 ()；H为磁场强度，单位为安培米()；B为磁通量密度，单位为 ()；J为电流密度，单位为()

9、；为磁流密度，单位为 ()。各向同性线性介质中的本构关系为D=E B=H J=E （2-1-9）其中表示介质介电系数，单位为法拉米(Fm)；表示磁导系数，单位为亨利米(Hm)；表示电导率，单位为西门子米(Sm)；表示导磁率，单位为欧姆米(m)。和分别为介质的电损耗和磁损耗。真空中0，0，以及=8.85F/m=H/m在直角坐标系中，麦克斯韦旋度方程（2-1-5）、（2-1-6）式写为 （2-1-10） （2-1-11） （2-1-12） （2-1-13） （2-1-14） （2-1-15）以上六个方程是FDTD离散差分的基础。 首先，将问题空间沿三个坐标轴向分成很多网格单元，用公、却和酝分别表示

10、在x、y和Z坐标方向的网格空间步长，用t表示时间步长，任意一个空间和时间的函数可表示为: （2-1-16）然后用中心差分式来表示函数对空间和时间的偏导数，这种差分式实质上是一种蛙跳法，具有二阶精度。为了实现空间坐标的差分计算，并考虑到电磁场在空间互相正交和铰链的关系，在FDTD离散中电场和磁场各节点空间排布如图2.2，这就是著名的Yee元胞。图 2.2 FDTD离散中的Yee元胞Yee元胞中E、H个分量空间节点与时间步取值的整数和半整数约定如表2-1所示。表2-1 Yee元胞中E、H各分量节点位置电磁场分量 空间分量取样时间轴t取样X坐标Y坐标Z坐标E节点 jk nkijH节点i jk下面直接

11、给出直角坐标中三维情形下的FDTD形式：设观察点（x,y,z）为的节点，即，以及时刻，于是（2-1-10）离散为 （2-1-17）式中（2-1-18）（2-1-19）同样（2-1-11）离散为（2-1-20）式中、可由（2-1-18）、（2-1-19）推得。同样（2-1-12）离散为 （2-1-21）式中、可由（2-1-18）、（2-1-19）推得。同样（2-1-13）离散为 （2-1-22）式中 （2-1-23） （2-1-24）同样（2-1-14）离散为 （2-1-25）式中、可由（2-1-23）、（2-1-24）推得。同样（2-1-11）离散为 （2-1-26）式中、可由（2-1-23）

12、、（2-1-24）推得。2.1.2数值稳定性 由于FDTD方程只是原Maxwell旋度方程的一种近似，在计算中存在误差。同时，由于FDTD方法是一个迭代过程，因此它的数值稳定性至关重要。（1）数值色散问题与空间步长 （2-1-27）（2）数值稳定性与时间步长 （2-1-28）界面上的介质参数采用介质边界算术平均值条件： （2-1-29） （2-1-30） （2-1-31） （2-1-31）对于理想金属导体界面，其切向电场为零。在Yee的差分格式中，已经给出电导率为的导体中的FDTD差分格式。其中，对电流项采用的是一般的中心差分格式。对于良导体，通常，因此，那么方程右边第一项系数为负数，接近于-

13、1，这将导致时间步进趋于不稳定。一般来说，为使解稳定，电流项的差分格式应介于中心平均近似和前向近似之间。Luebbers就曾对用代替中心平均近似使其差分方程对良导体稳定。此外，还可以采用指数差分格式。下面采用这两种差分格式离散麦克斯韦方程的旋度方程：（1）前向近似差分格式以（2-1-10）为例设观察点（x,y,z）为的节点，即，以及时刻， （2-1-32）令，有 （2-1-33）于是（2-1-10）离散为（2-1-34）式中 （2-1-35） （2-1-36）上式中标号。同样（2-1-11）离散为 （2-1-37）式中标号。同样（2-1-12）离散为 （2-1-38）式中标号。（2）指数差分格

14、式仍以（2-1-10）式为例，首先将（2-1-10）式变为方程 （2-1-39）解微分方程（2-1-39）式：令，有式中，整理得： （2-1-40）所以，即 （2-1-41）式中，上式中标号。同样（2-1-11）离散为 （2-1-42）式中标号。同样（2-1-12）离散为 （2-1-43）式中标号。用FDTD分析电磁散射、辐射等开发或者半开放性质问题时,受计算机内存容量限制,不可能直接对无限的空间进行计算,因此必须在截断处设置适当的吸收边界条件,以便用有限网格空间模拟开发的无限空间.目前对吸收边界条件的比较系统和深入的研究,主要是沿着两个方向进行的： 一是通过波动方程的因子分解获得单行波方程并

15、取近似来建立吸收边界条件；二是在边界上引入吸收材料,电磁波在无反射地进入吸收材料后被衰减掉,如PML.下面介绍几种应用较为广泛的吸收边界条件.（1）Mur吸收边界条件一阶Mur吸收边界条件： （2-1-44）二阶Mur吸收边界条件： （2-1-45）图2.3 二阶吸收边界条件所涉及的10个节点 图2.4 同轴线阶梯模拟2.2同轴线馈电单极子天线建模实例同轴线的矩形网格阶梯模拟对于具有复杂几何外形的目标，首先将其按几何外形的特点分解成几个部件，然后对各个部件分别建立几何外形及尺寸的描述文件并进行FDTD剖分。虽然同轴线的几何外形并不复杂，但这种思想仍然可取。实际上同轴线就是有几个圆柱体叠合拼接而

16、成，因此在对其建模时可以将其拆分为对多个圆柱体的建模，然后对重叠部分别处理。同轴线的矩形网格阶梯模拟如图2.4所示，模拟的原则是:以网格的中心点为准，如果网格中心点位于位于同轴线的内外导体之间，此网格即为空气网格，否则此网格为导体网格。具体算法见附录1。同轴线与同轴线馈电单极子建模基本思想基本一致，这里不再赘述。 同轴线激励源的最简单设置是在内外导体之间的某个位置加入激励电场，但由于同轴线中的主模是TEM模，如果加入的电场与TEM模不完全匹配，则在激励源附近就有相应的高次模存在，经过一段距离的传输之后，高次模逐渐衰减至可以忽略不计的程度，传输线上只剩余主模，此时的场分量才能够作为同轴线的入射波

17、加入到所要计算的同轴馈电的天线计算区域中。如果同轴线上强迫激励平面所加入的电场与TEM模基本匹配，则激励源附近产生的高次模就很小，经过很短一段距离的传输，高次模就可衰减至忽略不计的程度，同轴线上的场就是我们需要的主模场。 激励源的时间函数选用Gauss脉冲。如图2.5轴线上TEM模电场只有方向分量，而FDTD网格所用的坐标系是直角坐标系，因此强迫场只能通过激励源平面电场Ex、Ey加入，为使强迫场与同轴线主模电场一致，Ex与Ey合成场的方向应是同轴线半径方向，并且幅度与所在位置的半径成反比，具体的计算公式如下: （2-2-1） （2-2-2） （2-2-3） （2-2-4）图2.5同轴线场矢量分解示意图由此知道，激励源设置算法，见附录2。 天线结构空间的计算网格划分为总场区和反射场区。如图2.1在总场区中，所有网格点用Maxwel