l 引 言

时域有限差分法(fnite difference time domain，fdtd)是1966年k．s．yee首先提出的一种以maxwell方程为基础的解决电磁场问题的数值计算方法[1]。fdtd算法将maxwell方程中的两个旋度方程直接转化为差分形式，将电磁场进行空间和时问上的离散化，得到电磁场演化的迭代方程组，实现在一定体积内和一段时间上对连续电磁场的数据取样压缩。

应用fdtd算法对某一个具体的目标进行电磁仿真时，要将目标分解和描述成fdtd仿真软件可以处理的数据格式。在运用fdtd算法早期，由于受计算资源和图形显示设备的局限，用于仿真的建模数据多采用文本编辑形式，对于较复杂物体，很难及时发现建模过程中出现的错误，这种错误将被延续到仿真过程中，而专门的可视化建模模块可以使建模的效率和准确性得到提高。为了使仿真程序具有在不同操作系统的通用性，将电磁学模拟软件分为建模及结果处理模块、仿真模块两个部分，建模及结果处理模块专门进行目标物体的几何和电磁学参数的预设及相关数据的预处理，对模拟结果进行可视化显示；仿真模块专门对建模生成的数据文件进行相应的fdtd电磁模拟。

最后，以具体实例展现该建模模块的可视化效果、友好的图形化界面、结果的精确性和可靠性。

2建模算法

2．1建模物体的描述

对目标物体进行电磁学建模就是对目标物体的介质参数、几何参数、边界条件、激励源、输出记录等问题的描 述过程，并进行数据的预处理，生成仿真模块处理所需的数据文件。

(1)首先将建模目标物体中使用到的介质参数进行编号并生成介质列表。

(2)然后对目标物体进行几何参数设定，根据fdtd算法和目标物体自身结构特点，在直角坐标系下将物体分解单个矩形块输入。几何参数的输入有两种方式，一种是预先规定各个方向单个网格的尺寸，输入矩形的坐标以网格为单位；另一种是直接输入目标分解的矩形块的实际尺寸。两种情况都要给出矩形区域的介质编号。

(3)对目标物体根据具体问题对计算空间边界选择相应的边界条件，有一阶和二阶mur、一阶和二阶dispersive，pmc，pec，pml等7种边界条件供选择。

(4)目标物体激励源的设定，建模软件激励源采用软激励形式，将激励源信号加在目标物体网格坐标上。在完成上述工作以后，指定用于记录输出的网格坐标，用于记录仿真数据。

2．2 非均匀网格划分算法

为fdtd仿真程序提供尺寸合适的均匀和非均匀网格是建模的重要任务，在物体边界电磁场变化较大的地方宜采用细网格，而在物体边界电磁场变化比较小的地方宜采用粗网格[2,3]；考虑算法的数值色散的影响，任一网格的尺寸δ≤λmin/12，λmin为关心的上限频率所对应的波长[4]。为了节约计算资源和提高工作效率，对一些较复杂的目标物体采用非均匀网格划分方式。非均匀网格的划分遵守下列规则[5,6]：

(1)目标物体任一子区域的边界在网格线上。

(2)在各个方向上任意两个相邻网格尺寸比值≤1．25。

进行非均匀网格划分，根据网格密度分布方式有："密→疏"、"疏→密"、"密→疏→密"等3种分布方式。这里将盼"密→疏→密"非均匀网格为例介绍网格划分算法。

将目标物体在各个方向的几何参数排序划分成若干不重复的子区域，然后在各个子区按上述网格划分规则进行划分。在该模块中，用户可以预先设定3个方向中各个方向最大网格尺寸、最小网格尺寸和相邻网格尺寸比例，然后系统对3个方向网格进行一次性划分，并对网格划分进行规则检测，对用户进行提示，然后用户可以根据各个子区域情况进行优化。这里以x方向上某一子区域上进行网格划分为例来说明。设x方向某个子区域的长度为l，x方向最小的网格尺寸为xmin，x方向最大网格尺寸xmax，x方向预设相邻网格比例r∈[1，1．25]。网格的尺寸由区域的一边向中间逐渐增大，然后向另一边逐渐减小。

取该区域长度的一半l/2，令△x=xmin，由：

取k=min｛△xrm，xmax｝考虑几种情况：

①若l＇≥1．5k，将l＇均分为2个网格；

②若1．5k>l＇≥k，则l＇=k+△l＇，调整r使△l＇=0：

③若l＇<k，调整r使l＇=0。

这样在①条件下生成2m+2个网格，②条件下生成2m+1个网格，③条件下生成2m个网格。

(2)假如△xrm-1>xmax，找出n1的值，满足,△xrn1－1<xmax(n1<m)，在整个子区域内有2n1个网格满足条件。剩下的部分以xmax为网格大小进行均匀划分，得到：

对式（4）取整，得到n2=int(n2′),考虑两种情况：

①若n2′- n2≤0.5,调整r满足