常用数值计算方法及仿真软件简介a

1、1.1.1常用数值计算方法自1864年麦克斯韦建立电磁场基本方程以来，电磁波理论与应用的发展已 经过了 100多年的历史。对电磁分布边值问题的求解从图解、 模拟、解析到目前 所采用的数值计算方法，经历了四个过程。解析方法只能解决一些经典问题， 具 体到复杂的实际环境，往往需要通过数值解得到具体环境中的电磁波特性。随着高速和大容量计算机技术的飞速发展，电磁数值计算已经发展成为一门新兴的重 要学科，已提出多种实用有效的求解麦克斯韦方程的数值方法，主要有矩量法(MOM )、有限元法(FEM )、有限积分法(FIT)、和时域有限差分法(FDTD) 等。基于这些数值计算方法开发出了许多优秀的电磁仿真软件

2、。一个好的数值算法可以很接近地模拟出微波器件的特性，这对于工程设计和研究而言，可以避免很多次的“ cut-and-try”(试凑)，节省时间从而提高了效率。求解电磁问题的最终要求就是获得满足实际条件的 Maxwell方程的解，借助 于计算数学中的数值算法能够得到大多数电磁问题的近似解。数值算法的基本思 想就是把连续变量函数离散化，把微分方程化为差分方程；把积分方程化为有限 和的形式，从而建立起收敛的代数方程组，然后利用计算机技术进行求解。目前常见的几种数值分析方法如表 错误！文档中没有指定样式的文字。-1 电磁数值算法分类所示。针对本论文所应用到的方法，下面简要叙述常用的几种 数值方法及相应的

3、商业软件。表错误！文档中没有指定样式的文字。-1电磁数值算法分类数值算法频域方法积分形式矩量法(MOM)微分形式有限差分法(FDM) 有限元法(FEM)时域方法积分形式有限积分法(FIT)微分形式时域有限差分法(FDTD)高频近似绕射几何理论(GTD)一致性绕射(UTD) 一致性渐进(UAT)绕射物理理论(PTD) 绕射谱理论(STD)1.1.1.1有限元法基于有限元方法（FEM）计算电磁问题，其基本构想是将由偏微分方程表征的 连续函数所在的封闭场域划分成有限个小区域，每个小区域用一个选定的近似函数来代替，于是整个场域上的函数被离散化，由此获得一组近似的代数方程，并 联立求解，以获得该场域中函

4、数的近似数值。广义的来说，三维麦克斯韦方程是三维电磁问题的三维支配方程，但是，一般情况下为了方便求解和建模，大多选取由麦克斯韦方程组的前两个旋度方程导 出的电场强度满足矢量亥姆赫兹方程作为支配方程。如An soft HFSS软件川的支配方程为：（I一）一Vx 一g E k；®E=0（错误！文档中没有指定样式的文字。-1）由变分原理，上式的泛函可以写为：f"A彳imIF（E ）=仃:&一个亦E）（亦E）k；grE，E>d0（错误！文档中没有指定屮J样式的文字。-2）将这一个三维问题的泛函通过多面体离散成单元小矩阵，矩形块、四面体和六面体等都可以被选用做基本的离散

5、单元，但是，不同离散单元对于有限元运算的精度、速度和内存需求都有不同。 Ansoft HFSS软件采用四面体作为基本离 散单元，如图 错误！文档中没有指定样式的文字。-1所示，并选用上一世纪80 年代以后才被应用于电磁学中的棱边元作为矢量基函数。假设图 错误！文档中没有指定样式的文字。-1所示的四面体内的未知函数e能够近似为e =ae - bex cey dez（错误！文档中没有指定样式的文字。-3）34#-1四面体单元（错误！文档中没有指定样式的文字。-4）（错误！文档中没有指定样式的文字。-5）（错误！文档中没有指定样式的文字。-6）（错误！文档中没有指定样式的文字。-7）（错误！文档中没

6、有指定样式的文字。-8）（错误！文档中没有指定样式的文字。-9）图错误！文档中没有指定样式的文字。则用四个顶点处的值*匚（i =1,2,3,4）来表示:4炉（x,y,z ）=送 L：（x, y,z 黑ei 二式中插值函数L：（x, y, z为Le x,y,z茶 ae - bex qey dz6V而aie,bie,ce,die有下列等式获得：ae 决 a； ； - a27 - a； - a：；6Vbe =丄 be ；e be;七身 b：:6V1Ce = ；e ce ie-c37 - c： ；6Vede = ；e die ie d； ； d爲 d； ：e6Ve其中#VeX1ey1eZ1eX2ey2

7、eZ2eX3ey3eZ31ey4eZ4(错误！文档中没有指定样式的文字。-10)5#在利用变分原理和离散化方法建立了有限元矩阵方程后，我们就面临着求解以结点值为未知数的矩阵方程。将方程写为Ax二b(错误！文档中没有指定样式的文字。-11)式中系数矩阵A是一个nXn方阵，x是待求解的未知量，b表示已知向量。 求得这个矩阵方程得到问题空间的电磁场解， 并进而求得所需参数，例如散射参 数等。1.1.1.2有限积分法有限积分法(FIT)法早在1977年由托马斯魏兰特教授(Prof. ThomasWeila nd)冋提出，进而成为其后在电磁仿真领域中一个重要算法的基石。由FIT所导出的矩阵方程保持了解析

8、麦克斯韦方程各种固有的特性，女口：电荷守恒性和能量守恒性。解析下的梯度、散度和旋度算子在FIT下具有一一对应的矩阵。这 些矩阵满足解析形式下的算子恒等式。故FIT保证了非常好的数值收敛性。另一个区别于其它算法的关键之处在于 FIT可被用于所有频段的电磁仿真问题中。软 件CST Microwave Studio是由德国CST公司开发的一种商用电磁仿真软件，主 要是是采用这种数值计算方法。CST空间离散化也是建立在 Yee网格基础之上，典型划分方法如图错误！文档中没有指定样式的文字。-2示问:图错误！文档中没有指定样式的文字。-2麦克斯韦积分方程离散化图通过这种离散方法，导出相对应的麦克斯韦网格方

9、程: E dS = -A :：i iB dS:二ACe 二b（错误！文档中没有指定样式的文字。-12）H dS 二(:D一J ) dS一Ch 二d j（错误！文档中没有指定-AA一t样式的文字。-13）B dA = 0VSd =q（错误！文档中没有指定样V式的文字。-14） D dA =QASb =0（错误！文档中没有指定样式的文字。-15）经过这些步骤，将积分方程转化为线形方程组来求解， 得出问题空间的电磁 场量。CST软件包含了四种求解器：瞬态求解器，频域求解器，本征模求解器， 模式分析求解器，都有各自最适合的应用范围。瞬态求解器由于其时域算法，只 需要进行一次计算就可以得到在整个频带内的

10、响应， 该求解器适合于大部分高频 应用领域，对宽带问题优点尤为突出。 对于高谐振结构，例如滤波器，需要求得 本征模式，可以使用本征模求解器，结合模式分析求解器可以得到散射参量。对 结构尺寸远小于最短波长的低频问题，其频域求解器最为有效。1.1.2各算法比较对于众多求解电磁问题的算法和应用软件， 需要针对求解目标的实际情况来 选择合适的方法，这一步骤往往起到事半功倍的效果。 显然，计算机硬件能力是 要考虑的因素，求解问题的大小、复杂程度等对软硬件的要求不同。图错误！文档中没有指定样式的文字。-3给出了多种算法占用CPU计算时间的比较，当 结构简单且电尺寸较小时，需要划分网格数也少，多种算法占用时

11、间相当，但随 着计算问题网格数的增加，矩量法（MOM）占用时间与网格数呈三次方的关系增 加，有限元法（FM）是平方关系，而时域有限积分方法（FITD）则其CPU时间与网 格数几乎呈线形关系。另外，高频近似法适合计算电大导电体的电磁场问题， 对于电小结构及介质 或者其他非金属材料构成的系统则存在困难。 相反，各种频域方法和时域方法则 适合在低频或者谐振频率附近使用，可以用于近场和远场的计算。QT时间乌节直兼平力咸正蔭有限虫分图错误！文档中没有指定样式的文字。-3各算法比较1.1.3常用仿真软件作为目前电磁问题主要分析手段，电磁场数值计算方法为国内外广大工作者 所研究，并且随着这些数值方法研究的日

12、趋成熟， 大量商业化计算软件工具不断 涌现。随着应用开发的深入，其功能越来越强大，使用也越来越方便，这为具体 电磁问题的设计分析提供了极大的方便，也使包括天线在内的微波器件的设计周 期大为缩减。这对工程应用类的研究设计人员来说如虎添翼，可以很快地实现和 验证自己的创新设计思想。由于电磁场仿真软件与其核心的数值计算方法密切相关， 不同的软件其适用 的问题也不同。目前，可供选择的电磁仿真软件种类众多， 每种软件都有自己的 优势和劣势，可以根据软件的特点和设计的不同要求来选择软件。表 错误！文档中没有指定样式的文字。-2列出了目前最为常用的几种电磁 仿真分析软件的功能和应用，为合理选取软件分析实际问

13、题提供依据。 这些软件 的出现，使得微波电磁结构的设计可以在电脑上进行， 不但简化了之前的工作难 度，而且减少了重复加工测试的步骤， 大大降低了元器件的设计成本。 微波工程 师在设计各种器件并使用电磁仿真软件时，应该预先了解各种软件的基本的算 法，适用场合以及软件设置，以便达到最好的仿真效果，提高工作效率。表 错误！文档中没有指定样式的文字。-2常用电磁仿真软件的比较软件名称开发商数值算法功能特点应用领域Adva need Desig n System (ADS)Agile ntMOM可实现包括时域和频域、 数字与模拟、线性与非线 性、噪声等多种仿真分析 手段，并可对设计结果进 行成品率分析与

14、优化射频和微波电路的设 计，通信系统的设计， DSP设计和向量仿真An soft Desig nerAn softMOM米用“按需求解”技术， 将高频电路系统、版图和 电磁场仿真工具无缝地集 成到同一个环境中射频和微波电路的设 计，通信系统的设计， 电路板和模块设计， 部件设计An soft HFSSAn softFEM三维结构电磁场仿真软 件，拥有空前电性能分析 能力的功能强大后处理器天线分析与设计Microwave OfficeAWRMOM采用“场”和“路”两种 分析方法对不同的电路结 构进行仿真，加快了运算 速度射频集成电路、微波 单片集成电路、微带 贴片天线和高速印制 电路的设计XFD

15、TDRemcomFDTD三维全波电磁场仿真无线、微波电路、雷 达散射计算，化学、 光学、陆基警戒雷达 和生物组织仿真Zeland IE3DZela ndMOM可以解决多层介质环境下 三维金属结构的电流分布 问题微波射频电路、多层 印刷电路板、平面微 带天线的分析与设计CST MicrowaveStudioCSTFIT高频三维电磁场仿真软 件,除了主要的时域求解 器模块外，还为某些特殊 应用提供本征模及频域求 解器模块移动通信、无线通信（蓝牙系统）、信号集成和电磁兼容等SonnetSonnetMOM提供面向 3D平面咼频电 路设计系统以及在微波、 毫米波领域和电磁兼容/电磁干扰设计微带匹配网络、

16、微带 电路、微带滤波器、HDI和LTCC的转换本文研究所采用的工具是基于有限元法的An soft HFSS和基于有限积分法的CST Microwave Studio，下面对这两种数值算法进行介绍。1.131基于时域有限积分的软件 CST MICROWAVE STUDIOCST MICROWAVE STUDIO iii是德国 CST（Computer Simulation Technology） 公司推出的高频三维电磁场仿真软件， 广泛应用于移动通信、无线通信（蓝牙系 统）、信号集成和电磁兼容等领域。微波工作室使用简洁，能为用户的高频设计 提供直观的电磁特性。微波工作室除了主要的时域求解器模块外

17、，还为某些特殊应用提供本征模及频域求解器模块。 CAD文件的导入功能及SPICE参量的提取 增强了设计的可能性并缩短了设计时间。另外，由于CST设计工作室的开放性体系结构能为其它仿真软件提供链接，使微波工作室与其它设计环境相集成。另外在CST微波工作室中，还引入了 CST的专有技术理想边界拟合（Perfect Boundary Approximation PBA）。它使得长方形网格中材料的填充形式可以任意（单连通或复连通）。由于此技术，CST软件不但保持了通常FDTD 的快速，而且还使其精度大为提高。即，带 PBA的FIT即快又准。具体地讲，CST MICROWAVE STUDIO采用FITD

18、，即先在时域计算，用一 个宽频谱的激励信号（方波或者高斯波都有）去激励模型，在时域计算然后去反 演到频域。系统的网络参数和场参数基本上是反演后的得到的。特点是可以计算 相当大的带宽结果，而不需要像用HFSS,可能要把大带宽分割后分别仿真。CST 计算过程中，由于没有FEM计算过程中矩阵求逆过程，计算时间和网格数成线 性增长关系，而FEM的是指数增长关系。CST的MWS从4.3版起，开始有了 大小网格嵌套技术，在曲面上细化六面体网格逼进曲面。这是其它FDTD套件 所没有的。 传统的电路仿真软件仿真是快速的， 但是，当考虑趋肤效应损耗和材 料的复杂性时，结果的准确性将受到大幅度的影响。 像CST的

19、3D仿真软件克服 了这种限制， 可以解决任意几何形状的下所建立的麦克斯韦方程， 包括复杂的材 料模式。软件最新版本为 CST 2008。1.1.3.2 基于有限元的软件 Ansoft HFSSAnsoft HFSSi是美国An soft公司推出的三维电磁仿真软件，是世界上第一 个商业化的三维结构电磁场仿真软件， 业界公认的三维电磁场设计和分析的电子 设计工业标准。 它可以分析仿真任意三维无源结构的高频电磁场， 可直接得到特 征阻抗、传播常数、 S 参数及电磁场、辐射场、天线方向图等结果。该软件被广 泛应用于无线和有线通信、计算机、卫星、雷达、半导体和微波集成电路、航空 航天等领域。Ansoft

20、 HFSS 采用自适应网格剖分， ALPS 快速扫频，切向元等专利技术， 集成了工业标准的建模系统， 提供了功能强大、 使用灵活的宏语言， 直观的后处 理器及独有的场计算器，可计算分析显示各种复杂的电磁场，并利用 Optometric 可对任意的参数进行优化和扫描分析。HFSS提供了一简洁直观的用户设计界面、精确自适应的场求解器、拥有空 前电性能分析能力的功能强大后处理器，能计算任意形状三维无源结构的 S参数 和全波电磁场。HFSS软件拥有强大的天线设计功能，它可以计算天线参量，如 增益、方向性、远场方向图剖面、远场 3D图和3dB带宽；绘制极化特性，包括 球形场分量、圆极化场分量、Ludwig第三定义场分量 和轴比。使用HFSS可以 计算：1）基本电磁场数值解和开边界问题，近远场辐射问题。2）端口特征阻抗和传输常数。 3）S 参数和相应端口阻抗的归一化 S 参数。 4）结构的本征