1、计算泄漏同轴线的传播与辐射特性一、题目A22 计算如图所示泄漏同轴线的传播与辐射特性同轴线外导体沿轴向连续开槽二、FDTD 方法分析泄漏同轴线FDTD 的主要思想是把 Maxwell 方程在空间、时间上离散化，用差分方程代替一阶偏微分方程，求解差分方程组，从而得出各网格单元的场值。2.1 泄漏同轴电缆简介泄漏同轴电缆（Leaky Coaxial Cable）是遵循特定的电磁场理论，沿着同轴电缆的外导体周期性或非周期性开槽而形成的。信号在该电缆中传输时，能把电磁能量的一部分按要求从特殊开槽口以电磁波的形式辐射到周围外部空间，因此泄漏电缆能实现与外部空间的全方位双工通信。利用泄漏同轴电缆既具有传输
2、线的特性又具有无线电发射天线的性质，人们近年来广泛地利用它来解决无线电无法传输信息区域的信息传输问题。本文研究的是同轴线外导体沿轴向连续开槽的问题。2.2 Maxwell 方程在介质基片和空气中，Maxwell 旋度方程为：（1）JtDH（2）mtBE其中， 为电场强度，单位为 V/m, 为电通量密度，单位为 C/m2, 为磁场强度，单位为 Wb/ , EDH2m为电流密度，单位为 , 为磁流密度，单位为 。J2/AmJ/V各向同性介质中本构关系为：， ， ， （3）EDHBEJHm其中， 表示介电系数，单位为 , 表示磁导系数，单位为 ，/Fm/表示电导率，单位为 ， 表示导磁率，单位为 。
3、S和 分别为介质的电损耗和磁损耗。m在真空中(4)12078.5/4mFmH2.3 圆柱坐标系中差分格式的建立在圆柱坐标系中， （1） 、 （2）式可写为：（5）()1zrrrzrzzHEttErt以及（6）()1()1zrmrrzrzmzEHtt柱坐标下的 Yee 元胞如下图所示：图 1 柱坐标下的 Yee 元胞由图 1 可见，每一个磁场分量由四个电场分量环绕；同样，每一个电场分量由四个磁场分量环绕。这种电磁场分量的空间取样方式不仅符合法拉第感应定律和安培环路定律的自然结构，而且这种电磁场各分量的空间相对位置也适合麦克斯韦方程的差分计算，能够恰当的描述电磁场的传播特性。此外，电场和磁场在时间
4、上顺序上交替抽样，抽样时间间隔彼此相差半个时间步，使麦克斯韦旋度方程离散以后构成显式差分方程，从而可以在时间上迭代求解，而不需要进行矩阵求逆运算。因而，由给定相应电磁问题的初值及边界条件，利用 FDTD 方法就可以逐步地求得以后各个时刻空间电磁场的分布。取差分近似后 (5) 式 的 FDTD 离散式为rE(7)1 11/21/2 1(,)(,)(,)(,)22(,)(,)(nnnnrr rrnnz zijkijkEijkEijktHijHijir A(7)式整理后得：(8)1/21/211/21/2(,)(,)(,)(,)(,)(,)nnz zn nr rnnijkHijkttEijkEijk
5、irHijHijkz 同理可得：（9）1/21/211/21/2 1(,)(,)(,)(,),(,)nnr rn nnnz zijkHijkttEijkEijkzHHijr （1011/ /2()(,)(,)1 2(,)(,)221/,niijkiijkttnEijkEijkzz rnnijkrri）以上为电场的 FDTD 递推公式，同样可得到磁场的 FDTD 离散式为（111/2 1/2 11(,)(,)2222(,)(,)(,),nnzzn nmr r mn EijkEijkt tHijkHijkrEijkEz ）（12）1/2 1/2 11(,)(,)2222(,)(,)(,),nnrz
6、n nmmnzz EijkEijkt tHijkHijkEijkEr （131/2 1/2 11(),)(,)2222(,)(,)(,),12nnn nmz z mnnrr iEjkiEjkt tij ijkrEijkij ）2.4 数值稳定性分析FDTD 方法是以一组有限差分方程来代替麦克斯韦旋度方程，即以差分方程组的解来代替原来电磁场偏微分方程的解。只有离散后差分方程组的解是收敛和稳定的，这种代替才是有意义的。收敛性是指当离散时间隔趋于零时，差分方程的解在空间任意一点和任意一时刻都一致地趋于原方程的解。稳定性是指寻求一种离散间隔所满足的条件，在此条件下差分方程的数值解与原方程的严格解之间的
7、差为有界。在圆柱坐标系下，对空间步长和时间步长的要求比在直角坐标系下更为严格。沿同一径向，在角度步长相同时，横向长度单元是不同的，因此，为了保证数值收敛性，要取最小的值来计算。因为不同径向距离的横向单元长度为 ，但是考虑到电场和磁场相差半个时间步长，可取最靠近轴向的网格ir的横向单元长度，即取 ，此时时间步长的取值为 （14）1/2i 2221()(tcrz由（14）式可知，在轴向 ，在计算中会出现奇点，同样， 在轴向也存在奇点。0,0irirH这些奇点不是 Maxwell 方程本身所固有的，而是由微分方程转变为差分方程后带来的，可以用安培定律去除，因为 （15）csHdlEdStAA取上式的
8、积分得： （16）200(,)(,)zrethrztt式中， ，将（16）式转换为差分格式为：0/r（17）1 1/24(,)(,)(,)2nnnzztejkejkhjkr 有了 后， 及其他格式将自动满足。z 0,zj2.5 吸收边界条件：由于天线是一个无限空间的问题，而数值计算只能在有限空间进行，因此在计算空间的边界时必须设置吸收边界来模拟无限空间，使场满足辐射条件。为了让这种有限空间与无限空间等效，需对有限空间的周围边界面做特殊处理，使得向边界面进行的波在边界处保持“外向行进”的特征、无明显反射现象，并且不会使内部空间的场产生畸变。本文采 Mur 吸收边界条件。在圆柱坐标系中，三维的二阶
9、 Mur 吸收边界条件的差分格式为：1） ：0（18）1 11(,)(,)2(0,)(,)3(0,)(,)40,15,6,7)8nnnnnnnnnnjkcjkcjkjkcjkjkjjj j 其中, ， ， ，1/204ctt2()/04cctt， ， ，222( )/08cctz234/0256/cz， 。2117()/4tt 2118()/tt2） ：max0N(19)110 00000 1(,)(,)2(,)(,)3(1,)(,1)4 5,6,7nnnnnnjkcjkcNjkjkcNjkjkc 8,其中， ， ，21104ctt2()/04cctt， ， ，222( )/08tz23/25
10、6/0cz， 。2117()/4cctt2118()/04ctt3） ：min0（2011 11(,0)(,)2(,0)(,)3(,0)(,)4,15,60,7)8nnnnnnnnnnikcikcikikcikikiiii ）4） ：max0L（21110 00000 100(,)(,)2(,)(,)3(1,)(1,)415, 6,7nnnnnnnnnikcikciLkikciLkikiii 08,ni）式（20）和（21）具有相同的系数：， ， ，210ctt21()/0cctt22()/0ctz， ， ，23()/0424()/4256/c， 。217()/cctt815） ：min0z（
11、221 11(,)(,1)(,0)(,)3(,0)(,)4,15,6,78nnnnnnnnnjcijcijijcijijjjijijici ）6） ：max0zM（2310000001(,)(,1)2(,)(,)3(1,)(1,)4 5, 6,7nnnnnnijcijcijMijcijMijijijij 8,ij）式（22）和（23）具有相同的系数：， ， ，210ctzvt21()/0cctzt22()/0ct， ， ，23()/0424()/256/c， 。217()/0cctzt812.6 激励源：采用 Gauss 脉冲激励， （ 24）20()exp()/iEttT式中 是出现最大值的
12、时间，T 与脉冲宽度宽度有关，应适当选择 和 T，以保证源的初始条件，即0t 0t。对式（24）进行傅里叶变换得 （25）()iE2()exp()iEff本文取 ， 。1.67es03tT2.7 近远场外推：2.7.1 基本思想由于 FDTD 方法只能计算空间有限区域的电磁场，要获得远区的辐射场则必须应用等效原理，如图 2所示，在源的周围引入虚拟界面 S，而 S 外是自由空间。令 S 之外存在原来的场，在 S 之内的场为零，且到处是自由空间，为支持这样的场，S 面上必须存在等效面电流 和面磁流 。sJmJ（a）原问题 （b）等效问题图 2 等效原理示意图,smJnHEn三维自由空间的格林函数为
13、： （26） exp(|)(,)4jkrGr其中 分别为观察点和源点位置失量，,r取远区近似得： ，其中 为 方向单位矢量，|rreAr|r于是， （26）式近似为 （27） xp()(,)e()4rjkGjA（28） (),),mArJrdVFrJGdV将式（27）代人式（28）得： （29） exp()()(exp)4)rsmrsjkrArJjkdsjFjrA电流和磁流的辐射场为 （30）11()EAFjAj jHAF在远区， （30）式中的算子 可以用 来代替，因此： （31）()jk ()kEjkjAjHjFj将（31）式右端写成球坐标分量形式，从而可以得到：（32）EjkFA（33）
14、Hjkj2.7.2 三维瞬态场的外推如果 FDTD 计算时入射波是时域脉冲的形式，就需要考虑如何将时域近场数据外推到计算区域以外的情况，对于远区场，由式（33）可得到： （34）()EjkZAF下面将频域远区场外推表达式（34）通过傅里叶变换转换到时域。为此在频率中令：（35） exp()()(exp)4)rsmrsjkrWjkArjJjkdsjUjFj jr A其中， ，记其傅里叶变换为： /,(),)(,mkcJrJ（36） )exp)(,)(,(smmsjrtrjtdjtJjt根据傅里叶变换性质，有：（37） (,)(,)exp(),smmsjrtjJrjtdjtjjt以及 （38） (
15、,)(,)exp()()exp()(,)(,)()()()rsrmmsejrtJjkrjrjtdcjtrjjjtAA由（37） 、 （38）式得（36）式的傅里叶变换为： （39） 1()(,)4()(,)rsrmsewtjtdsrctcutjttA所以，将（34） 、 （35）式并作傅里叶变换，得到时域中的关系式为： （40）)(euZwt其中， 是电场 的傅里叶变换。()etE三、FDTD 方法计算泄漏同轴线的传播与辐射特性3.1 用 FDTD 方法求解泄漏同轴电缆传播和辐射场的步骤：1、将时域 Maxwell 的旋度方程展成其坐标分量式，本文采用圆柱坐标系，用中心有限差分格式代替各场分量
16、对空间、时间的微分，得到 FDTD 的基本方程。2、确定空间网格基本单元尺寸 ，单元尺寸大小可以相等，也可以不等，为了减小数值色散，,rz应满足 ， ，本文中，取 。min() /10rz3、选择时间步长 ，根据前面的讨论计算得到： ，t 2221()(tcrz本文近似取 。/2trc4、沿网格空间的外表面设置吸收边界条件。5、选用和设置吸收边界条件。6、确定运算的总时间步数，对于时域计算，应使激励源主要频谱部分的信号达到稳定状态，特别应使低频谱信号达到稳定状态。7、计算边界处的电磁场，利用等效原理和近远场外推的方法，求得辐射场。8、如图 2 所示，利用 A，B 两个横截面上的电磁场分布，计算
17、出相应的功率，从而可推导得到传播常数。图 2 泄漏同轴电缆的侧面图注：计算 A、B 两个横截面处的电磁场，是为了避免在激励附近高次模的影响。3.2 FDTD 方法的 MATLAB 编程思路：1、按所传输波频率大小以及所仿真物理模型大小合理划分网格，一般取所传输波最高频率所对应波长的1/10-1/20，还要结合模型自身尺寸来确定空间步长。2、按Courant稳定条件确定时间步长和迭代循环次数N 。3、赋予各点相对应的物理参数，包括介电常数、磁导率、电导率等。这里需要注意的是，当网格处于两种不同介质交界面时，计算该处电场，取介电常数，电导率为在两种介质中的平均值。4、计算电磁场迭代时需要用到的系数，5、循环开始，计算电场(或磁场 )，加入激励源，计算磁场 (或电场)，设置吸收边界。6、用 MATLAB 语句显示场图。四、HFSS 仿真结果泄漏同轴线参数：内半径 r1=1cm, 外半径 r2=2cm, f=3G, 介质层相对介电常数 ，沿着轴向连1.25r续开槽，长度 L=40cm,轴向宽度 30。1.远区辐射方向图图 3 3D 辐射图
