CST 软件仿真超材料特性

2025-10-09
广东
本文字数：2500 字
阅读完需：约 8 分钟

简要介绍

人工电磁材料或称之为超材料，已被广泛应用于天线、透镜、隐身斗篷以及微波无源器件的研究设计中。具体可分为 DSG、EBG、FSS 等等。高阻抗表面(HIS, high impedance surfaces)是广泛应用的一种人工电磁材料。下面以 HIS 为例分析使用 CST 软件仿真超材料的各个参数。

高阻抗表面具有独特的电磁特性：

一方面，当平面波入射时（kx^2+ky^2 ≤ k0^2，kz 为实数，其中，kx，ky，kz 分别代

表 x, y, z 方向上的波数），HIS 的反射相位为零，这就是同相位反射特性，与 PMC 特性相同。与之对应的是 PEC，反射相位为 180°。反射特性如下图所示。

另一方面，由于 HIS 的周期性结构特征，能支持多种空间谐波。按照 HIS 的具体结构的不同，在表面波传播特性方面可以划分为有禁带的 HIS 和无禁带的 HIS 两种类型。有禁带的 HIS 表现出电磁带隙特性：当表面波入射时（kx^2+ky^2 ≥ k0^2，kz 为虚数，其中，kx，ky，kz 分别为 x, y, z 的波数），HIS 呈现出频率带隙，在这一带隙之内，任意极化方向的波均不能传播。HIS 的电磁带隙特点在抑制表面波、减少天线单元间互耦以及提高天线的增益等方面有着重要应用。此外，无带隙的 HIS 由于能够支持表面波的传播，因而被广泛应用于需要对表面波的传播进行控制的透镜、全息天线以及电磁隐身等方面[1]。

同向反射区 CST 软件仿真

1、建模：

如上图所示，我们建立最简单的一个晶胞结构。以笔者为例，上层片为 2mm，介质采用 Rogers5880，厚度 1.5mm，介质长 6.42mm，宽 4.9mm，最底层为完整地。

2、参数设置

（1）设置工作频率

Simulation->Frequency

包含中心频率即可，为较少仿真时间，也可缩小频率范围。

（2）设置边界条件

Simulation->Boundaries

由于我们只建模了一个晶胞，实际为在 X-Y 平面内的无穷大结构，于是边界条件中，在 X、Y 方向我们采用 unit cell 结构，在 Z 方向采用 open。

（3）设置求解器

Home->Simulation->Frequency Domain

将导航树下 ports 端口中的 Zmax 和 Zmin 分别双击，选择 mode 为 2，这样我们可以仿真两个模式，即 TE 和 TM 模式的传播特性。

其次我们更改 Zmax 的 distance 到结构的上表面，我们需要的是该结构上表面的同向反射特性，Zmin 不需要此要求。

最后，我们打开频域求解器，Source 采用 Zmax 激励，模式采用 All，Start 即可。

（4）得到同向反射区结果

仿真结束后，我们选择 1D result 下的 S-parameters 选择 SZmax(1),Zmin(1)和 SZmax（2）,Zmin(2)，接着选择

1D plot->Plot Type->Phase，即可得到下图，可以发现在超表面上表面的同向反射特性。

色散图仿真

关于布里渊区、色散图为何物，读者自行阅读文献。

考虑一下入射到周期性表面上的平面波。对于表面波来说，入射角等于 90°（与表面平行）。对于此值，无法使用平面波响应研究周期性结构上的波传播。通带和带隙在频谱中的确切位置只能通过沿简约布里渊区轮廓的表面波的色散关系（即计算本征模的谐振频率）来获得。

对上文叙述的结构，继续进行色散曲线仿真。

1、边界条件

如下图所示，本征模求解器不支持 unit cell、不支持开放边界，我们采用 periodic 来表示周期阵列。

选择 Simulation->Background，将 Upper Z distance 设置为 10 倍的介质厚度。

CST 软件仿真中很多实用经验表明，在背景材料和边界设置中，关键参数是介质基板上方空气空间的高度，以仿真结构上的自由空间（在“背景设置”中，该空间为“上部 Z 距离”），必须在晶胞上放置一个高度约为介质板厚度许多倍的空气盒。基于大量计算机仿真，并将获得的结果与分析和实验考虑因素进行比较，确定了大约十倍于基板厚度的正确选择，即上部 Z 距离=10*h。

接着将边界条件 Simulation->Boundaries 设置为上文所述。

模型顶部和底部的边界应定义为电导体或 PEC。

最后，我们选择 Boundaries 内的第三项 Phase Shift，新建两个变量 phase_x 和 phase_y，将参数“相位”分配给周期性边界，以便可以在参数扫描中使用相移，以测量该结构的频率色散行为，即用来扫参布里渊区。

2、色散图查看方式

（1）由用户定义监视器，并且可以构建自定义的监视器。用户定义的参数扫描监视，将第一模式的组速度、相速度和色散图添加到“1D result”文件夹中的导航树中。同样，可以在“1D result”文件夹中看到功率流和阻抗。在参数扫描窗口中，选择“编辑”以查看或修改源代码。

（2）数据后处理模板：

色散图是基于“相位”的，预先定义参数的“频率”与空间相位变化的关系图，可以通过以下逐步过程获得：

调用 TBP，然后选择“2D and 3D Field Results->3D Eigenmode Result”

模式选择 mode1 即可。

笔者采用（2）方式。

3、本征模求解器设置

求解器选择本征模求解器，method 采用 JDM，mode 选择 1。

接着选择 par.Sweep，选择 New Seq.->New Par，选择 phase_x 从 0-180 读扫描，样本值 19，即以 10°为步进扫描。

当然你也可以降低精度提高速度，样本值采用 7，扫描步进 30°。

接着 Start 即可。

4、得到色散图

布里渊区如下如所示。

沿着图中箭头所指方向遍历所有可能的传播矢，即Γ->X->M->Z->Γ，就可以求解出整个周期性结构的色散特性。令 phase_y=0°，以 10°为间隔从 0°递增到 180°，实现Γ->X；令 phase_x=180，以 10°为间隔从 0°递增到 180°，实现 X->M；令 phase_y=180°，以 10°为间隔从 180°递减到 0°，实现 M->Z；令 phase_x=0°，以 10°为间隔从 180°递减到 0°，实现 Z->Γ。

具体步骤为，首先按照上文扫描 phase_x=0°-180°，完成换成 phase_y=0-180°。接着导出数据。设置 phase_y=180°，扫描 phase_x=180°-0°，完成后换成扫描 phase_y=180°-0°，即可得到全部数据。

将数据导入 origin 绘制。同时由公式 f=c*k/(2*pi)，我们同样绘制参考光线，可以得到该超材料的色散图。