电磁波前沿问题仿真工具解析.doc

计算机水平的迅猛发展，为我们开拓了许多前沿研究的方法。进行前沿研究时，特别是像高频电磁学（包含微波、光学等范围）这种物理机理较为复杂的情形，传统只依靠实验来进行研究的方法，已经力不从心了。于是近十几年来，兴起了一种流行的研究方法，计算机辅助工程（CAE），用以辅助和验证我们的理论和实验模型，并形成了一个电磁学的分支计算电磁学。然而，大多数运用计算机仿真手段进行研究的人员，本质的工作还在于我们所研究的物理现象本身，除了少数人能够精通编程（或者说有时间专门研究计算电磁学）之外，大部分研究人员还是选择使用现有的仿真工具来进行研究。主流的工具包括OptiFDTD、Rsoft Photonics Suite、Lumerical FDTD Solutions、XFDTD、COMSOL Multiphysics、ANSOFT（HFSS）、CST等。这些主流软件使用的算法，可以归结为以下几种：时域有限差分法（FDTD）、有限元法（FEM）以及混合法（也就是各类数值计算方法的组合）。下面就来大概介绍和对比一下这些软件的一些功能和特点。OptiFDTD、Rsoft Photonics Suite、Lumerical FDTD Solutions、XFDTD这几个软件，从名字上可以一眼看出，他们都是使用FDTD算法的。FDTD的思想是将求解域划分成许多网格，把随时间变化的Maxwell微分方程转化为每个网格结点的有限差分方程，每个结点的电(磁)场仅与其相邻网格上的磁(电)场以及上一时间步长的场量有关。选取合适的边界周期条件和初始条件，求解四维数值解，这四维是指空间的三维加上时间维度。从算法的原理可以看到，这一类的软件，最大的特点就是能够一次计算完成扫频。FDTD计算的是时域波形的演变，而一个有一定波形的时域信号，它在频域上一定不是单频的，而是会有一定的谱宽。这样的最大好处便是一次可以计算复杂时域的波形。用来计算电磁波的传输时非常的方便。另一方面，由于FDTD算法是基于时域的，也就是说，计算时需要求解电磁波传输一段距离，这就可能带来一些困扰。比如，在进行光纤的模场分析时，选取光纤的一个截面进行分析。这个时候不需要关心电磁波在光纤传输方向的传输状态，但FDTD需要计算电磁波在这个方向上传播一段距离才能得出在截面上稳定的状态（往往要计算到较长距离才能达到这种状态），这无形中增加了不少的计算量。另一个重要的一点，在涉及到多尺度的问题时（微纳光学器件里经常会遇到），电磁波通过小尺度的求解域时，需要对时间进行细致的划分，这将会导致在大尺度求解域的划分，也必须按照这个标准（不均匀度不能超过一定限度）来进行，最终导致网格数太大而使得计算不可行。在做仿真时也需要考虑到这点。COMSOL Multiphysics这款软件是一个纯有限元的软件。最早听说这个软件和看这个软件的介绍时，觉得只是一个做多场耦合的通用软件，未必能够解决前沿电磁波的问题。但深入使用之后才发现这款软件的强大。COMSOL在求解电磁波问题时，是基于频域来求解的，采用的是一种谐波近似的方法。很多时候我们想计算电磁波的传输，但是其实我们并不关心这个电磁波的时域波形有多复杂，这种情况下，我们就可以用谐波近似，把电磁波近似正弦波。正弦波在频域上就是一个单频，因此COMSOL做这种计算很快。另外还有一种计算，就是模场分析。COMSOL直接在频域求解亥姆霍兹方程（波动方程的频域形式）即可。COMSOL最早是MATLAB的一个求解偏微分方程（PDEs）的工具箱，因此传承了MATLAB平台性、开放性的特征，能够对绝大多数的物理问题，使用数学符号来描述。而其本质，也就是求解各种形式的偏微分方程。比如，某种非线性很强的材料，介电常数与电场强度的平方成正比，在COMSOL里，可以直接将这个材料的介电常数直接写成电场的平方的表达式就可以了，求解时，会将这个表达式代入基本方程中介电常数的项，软件直接求解最后这个形式的方程。这对于求解非线性的问题是一个很大的亮点。在这个求解与解有关的材料系数问题里，COMSOL直接避免了一个时间先后的问题：是解先影响了材料，还是材料先影响了解？关于材料问题，目前研究得比较热门的是电磁波与金属的相互作用问题。金属介质的模型不准，或者当遇到金属介质的时候，收敛缓慢。金属的材料特性与普通介电介质太不一样。这种材料的不连续当然使得电磁场在界面上变化很陡峭。FDTD在计算时需要进行时间推演，场变得快则时间步就得更小，计算量大，收敛比较慢。在COMSOL没有这个问题，快变没有太大的影响，一个正弦波，变得再快，频率再高，也只是求解单频而已。前面讲到的多尺度问题，也是COMSOL的强项。得益于有限元的算法，COMSOL处理复杂结构模型时，能够将网格的不均匀度控制到一个很大的范围。也就是在细微处能够划得很密而在需要粗略处理的地方划得很粗。这样整个求解域的网格数就能控制到一个很有效的水平，从而节省计算量。最后再说说Ansoft和CST，这两个软件是也是在电磁场计算中非常有名的软件。Ansoft使用的是有限元和矩量法相结合的方法；CST是一个专门求解电磁问题的软件包，有FDTD、矩量法、传输线法等方法（没有有限元法FEM）。这两个软件，更多的是面向于理论比较成熟的工程应用，模块比较固化，对实际的问题只需对号入座便可，能够方便地提高学习效率和操作性。对于前沿理论的研究来说，灵活性和开放性就稍显不足，用户发挥的余地较小。以上就是我对目前一些主流计算电磁学软件应用于前沿电磁学研究的一些粗浅看法，限于学识和经验，不能理解和表达透彻，希望各位能够及时提点和修正。ANSYS 12.0电磁学有限元分析从入门到精通张倩, 胡仁喜, 康士廷等编著机械工业出版社 2010O441.4-39/062中文图书ANSYS 10.0电磁学有限元分析实例指导教程孔明礼, 胡仁喜, 崔海蓉等编著机械工业出版社 2007O441-39/048中文图书ANSYS 10.0工程电磁分析技术与实例详解阎照文等编著中国水利水电出版社 2006O441.4-39/847中文图书Ansoft工程电磁场有限元分析刘国强, 赵凌志, 蒋继娅编著电子工业出版社 2005O441.4/840中文图书时域有限差分法FDTD Method高本庆编著国防工业出版社 1995.3O241.3/717中文图书HFSS电磁仿真设计应用详解李明洋编著人民邮电出版社 2010O441.4-39/148中文图书FEKO 5.4电磁场分析技术与实例详解阎照文, 苏东林, 袁晓梅编著中国水利水电出版社 2009O441.4/847.1中文图书电磁散射与辐射问题仿真理论与方法(mom-po)童创明编著西北工业大学出版社 2010O441.4/764.1中文图书电磁兼容问题的计算机模拟与仿真技术周佩白 . 等 著中国电力出版社 2006TN03/666中文图书信号完整性仿真分析方法阎照文等编著中国水利水电出版社 2011中文图书EMC电磁兼容设计与测试案例分析.第2版郑军奇编电子工业出版社 2010TN03/881.01中文图书工程电磁场数值分析与综合谢德馨 .等编著机械工业出版社 2009O441.4/862中文图书电磁场与波分析中半解析法的理论方法与应用盛剑霓等著科学出版社 2006O441.4/161中文图书电磁场分析中的应用教学余恬, 雷虹编著北京邮电大学出版社 2009O441.4/682.3中文图书电磁场理论形式逻辑分析杨本洛著上海交通大学出版社 2009O441.4-53/218中文图书电气工程电磁场数值分析颜威利 . 等 著机械工业出版社 2005TM153/816中文图书工程电磁场数值计算.第2版倪光正, 杨仕友, 邱捷等编著机械工业出版社 2010.4中文图书工程电磁场数值计算倪光正 . 等 编著机械工业出版社 2004TB13/631中文图书计算电磁学的数值方法吕英华编著清华大学出版社 2006TM15/416中文图书电磁场数值计算与MATLAB实现何红雨编著华中科技大学出版社 2004O441.4/601中文图书微波与高速电路中的电磁场理论及其数值方法李征帆, 曹毅著科学出版社 2002TN710/164.1中文图书课程名称：ANSYS电磁场分析时间: 5天 上午9:00-12:00 下午：1:00-4:30 第一天 低频电磁场1 电磁场仿真概述2 二维静磁学3 实例：简单直流致动器4 模型边界条件5 分析单位6 BH曲线输入7 实例：轴对称直流致动器8 永磁体9 耦合及约束方程10 练习;第二天 二维交流和瞬态分析1 交流模拟概念2 二维轴对称有限元模型3 实例：电机槽内导体4 实数解与虚数解5 实例：带圆环的交流致动器6 电路单元应用7 2D变压器应用实例8 瞬态分析9 致动器-采用跃变电压第三天 三维电磁模拟1 三维模拟单元类型（标量势单元，单元边列式，磁矢量式）2 三维应用3 三维直流衔接磁铁4 三维直流致动器实例TEAM20标准检查程序第四天 耦合场模拟1 结构-磁耦合2 热-电磁3 实例分析：三维汇流排变形第五天 高频电磁场分析1 电磁场与电磁波介绍2 电磁场理论3 麦克斯韦方程组，波方程，电磁场边界条件，格林函数，等效面与近/远场,许瓦兹反射定律与镜像法，泊印亭定理，唯一性定理，叠加定理，散射矩阵4 高频电磁场分析的有限元方法5 麦克斯韦方程弱形式，总场与散射场，腔体与波导结构，基函数，边界条件与PML，激励源，参数提取6 ANSYS高频电磁场仿真工具7 导波结构8 谐振腔9 电磁波传播10 电磁波辐射11 电磁波散射光的电磁说说明光在本质上是电磁波的理论。电磁辐射不仅与光相同，并且其反射、折射以及偏振之性质也相同）由麦克斯韦的理论研究表明，空间电磁场是以光速传播。这一结论已被赫兹的实验证实。麦克斯韦，在1865年得出了结论：光是一种电磁现象。按照麦克斯韦的理论c/v=( * ) 式中c为真空中的光速。为在介电常数为和导磁系数为的媒质中的光速，因为c/v=n(折射率),所有n=( * ) 这个关系式给出了物质的光学常数，电学常数和磁学常数之间的关系。当时从上述的公式中看不出n应随着光的波长而改变，因而无法解释光的色散现象。后来罗仑兹在1896年创立了电子论，从这一理论看，介电常数是依赖于电磁场的频率，即依赖于波长而变的，从而搞清了光的色散现象。光的电磁理论能够说明光的传播、干涉、衍射、散射、偏振等许多现象，但不能解释光与物质相互作用中的能量量子化转换的性质，所以还需要近代的量子理论来补充。 光的微粒说关于光的本性的一种学说。17世纪曾为牛顿等所提倡。这种学说认为光由光源发出的微粒、它从光源沿直线行进至被照物，因此可以想像为一束由发光体射向被照物的高速微粒。这学说很直观地解释了光的直进及反射折射等现象，曾被普遍接受；直到19世纪初光的干涉等现象发现后，才被波动说所推翻。1905年提出光是一种具有粒子性的实物（光子）。但这观念并不摒弃光具有波动性质。这种关于光的波粒二象性的认识，是量子理论的基础。 光的波动说关于光的本性的一种学说。第一位提出光的波动说的是与牛顿同时代的荷兰人惠更斯。他在17世纪创立了光的波动学说，与光的微粒学说相对立。他认为光是一种波动，由发光体引起，和声一样依靠媒质来传播。这种学说直到19世纪初当光的干涉和衍射现象被发现后才得到广泛承认。19世纪后