对于COMSOL——雷达散射截面案例的分析

简介：

本教学案例演示背景场在电磁散射问题中的应用，描述的是被雷达检测到的船身，同样的技术可用于模拟孤立物体受来自远处的电磁波影响的任何情况。例如，将模型尺寸缩小几个数量级，同样常见的应用是等离激元共振纳米颗粒。本例除了设置背景场并在一定的入射角范围内对其进行扫描之外，还介绍了如何计算远场和雷达散射截面(RCS)。

模型定义：

本例计算船身与雷达发射器的入射场之间的相互作用。将发射器视为距离船身足够远,这样我们可以将入射场视为平面波,因而可以从模型几何结构中排除发射器,仅分析船身及其周围环境。
虽然二维下的建模过程与三维类似,但本例使用二维模式可以实现快速设置和求解。为了着重阐述概念,我们特意采用非常简单的几何结构(见图 1)。

图一：模型几何结构。船身长14米，被空气包围-界面位于水面之上。

几何结构中的内圆表示船身周围的空气，外圆是完美匹配层（PML）。为了得到比较好的结果,内圆需要完全包围船身,还必须延伸到波长的相当大一部分,及其外部倏逝波的特征长度。但在实践中,为了最大限度地减少所需时间和内存使用率,本应避免将圆绘制得过大。出于本例的目的,为了得到近场的清晰视图,内圆的尺寸超出了必要大小。外圆半径不重要,只要满足 PML 网格划分的需要(即 5-6 个网格单元)即可。雷达的背景电磁场由其面外电场分量描述:

在该方程中：

j是虚数单位；

k0 = 2f/c是真空中的波数；

f = 100MHz是频率，是模型的设定值，即条件；

是射入角，以 = 0 对应x轴正方向的入射波，在模型中即为从右向左沿x轴负方向传播。

建模操作说明：

从文件菜单窗口中点击新建

在不偏离预期实验效果的情况下，对三维的模型进行简化，只取其在海面（二维平面）上的截面，我们将入射场视为平面场，仅分析船身及周围的环境。使用二维模型简替三维模型话可以实现快速设置和求解，其主要目的归根结底还是为了着重阐述概念。

新建

在新建窗口中，点击模型向导

模型向导

1在模型向导窗口中，单击二维

2在选择物理场树中选择射频>电磁波，频域(emw)。

3单击添加。

4单击→研究。

5在选择研究树中选择一般研究>频域。 

6单击√完成。

研究1

步骤1:频域

1在模型开发器窗口的研究1节点下，单击步骤1:频域.

2在频域的设置窗口中，定位到研究设置栏。

3在频率文本框中键入"100 [MHz]”。

这一步是为了设置全局的雷达波频率为100MHz

全局定义

参数1 

1在模型开发器窗口的全局定义节点下，单击参数1。

2在参数的设置窗口中，定位到参数栏。

3在表中输入以下设置:

几何1

圆1 (c1) 

1在几何1工具栏中右键，单击圆。 

2在圆的设置窗口中，定位到大小和形状栏。

3在半径文本框中键入"r0+3[m]”。 

4单击以展开层栏。在表中输入以下设置:

圆弧1 (ca1) 

1在几何工具栏中右键，然后选择圆弧。 

2在圆弧的设置窗口中，定位到属性栏。 

3从指定列表中选择端点和起始角度。 

4定位到起点栏。在x文本框中键入"6”。

5定位到终点栏。在x文本框中键入"4",在y文本框中键入"2"。

圆弧2 (ca2) 

1在几何工具栏中右键，选择圆弧。

2在圆弧的设置窗口中，定位到属性栏。 

3从指定列表中选择端点和起始角度。 

4定位到起点栏。在x文本框中键入"4"。在y文本框中键入“-2”。 

5定位到终点栏。在x文本框中键入"6"。在y文本框中键入"0”。

6定位到角度栏。在起始角度文本框中键入“270”。 

线段1 (Is1) 

1在几何工具栏中右键，选择线段。 

2在线段的设置窗口中，定位到起点栏。 

3从指定列表中选择坐标。 

4定位到终点栏。从指定列表中选择坐标。

5定位到起点栏。在×文本框中键入“-4"。在y文本框中键入“-2”。 

6定位到终点栏。在x文本框中键入"4”。在y文本框中键入“-2”。 

线段2 (ls2) 

1在几何工具栏中右键，然后选择线段。 

2在线段的设置窗口中，定位到起点栏。 

3从指定列表中选择坐标。 

4定位到终点栏。从指定列表中选择坐标。 

5定位到起点栏。在x文本框中键入“-4”。在y文本框中键入"2”。 

6定位到终点栏。在x文本框中键入“4”。在y文本框中键入“2”。 

二次贝塞尔曲线1(gb1) 

1在几何工具栏中右键，然后选择二次贝塞尔曲线。 

2在二次贝塞尔曲线的设置窗口中，定位到控制点栏。

3在第1行中，将×设为-4并将y设为2。

4在第2行中，将×设为-6并将y设为2。 

5在第3行中，将x设为-8。 

二次贝塞尔曲线2 (gb2) 

1在几何工具栏中右键，然后选择二次贝塞尔曲线。 

2在二次贝塞尔曲线的设置窗口中，定位到控制点栏。 

3在第1行中，将x设为-8。 

4在第2行中，将x设为-6并将y设为-2。 

5在第3行中，将x设为-4并将y设为-2。 

此处与原案例pdf文件不同，照原pdf文件去做，结果是形同二次贝塞尔曲线1。

转换为实体1 (csol1) 

1在几何工具栏中右键，然后选择转换>转换为实体。 

2选择“对象” cal、ca2、Is1、Is2、qb1和qb2。

定义

显式1 

1在定义工具栏中右键，单击选择>显式。 

2选择“域”4。 

3在显式的设置窗口中，定位到输出实体栏。 

4从输出实体列表中选择相邻边界。  

在本例中，我们会多次用到船身的边界条件。把整个零散的部分组合为一个显式，可以方便我们进行后续一系列处理。

完美匹配层1(pmll) 

1在定义工具栏中右键，单击完美匹配层。 

2选择“域”1、2、5和6。 

3在完美匹配层的设置窗口中，定位到几何栏。

4从类型列表中选择圆柱型。 

添加材料

1在主屏幕工具栏中，右键从库中添加材料。 

2转到添加材料窗口。

3在模型树中选择内置材料>Air。 

4右键选中后，单击添加到“组件1（comp1）”。 

5在模型树中选择内置材料>Aluminum。

6右键选中后，单击添加到“组件1（comp1）”。 

材料Aluminum (mat2) 

1在材料的设置窗口中，定位到几何实体选择栏。 

2从几何实体层列表中选择边界。

3从选择列表中选择显式1。

电磁波，频域(EMW) 

1在模型开发器窗口的组件1 (comp1)节点下，单击电磁波，频域(emw)。 

2选择“域”1-3、5和6。 

为了简化模型，也为了节省时间和内存，使用边界条件来表示船身，令船内部的场为0，省去计算量和额外干扰。

3在电磁波，频域的设置窗口中，定位到公式栏。 

4从列表中选择散射场。 

5将Eb矢量指定为：

1[V/m]*exp(j*emw.k0*(x*cos(phi)+y*sin(phi)))

在这个表达式中，真空中波的传播常数emw.k0由物理场接口自动提供。

6定位到分量栏。从求解的电场分量列表中选择面外矢量。 

阻抗边界条件1 

1在物理场工具栏中右键，然后选择阻抗边界条件。 

2在阻抗边界条件的设置窗口中，定位到边界选择栏。

3从选择列表中选择显式1。

远场域1 

在物理场工具栏中右键，然后选择远场域。 

网格1 

在模型开发器窗口的组件1 (comp1)节点下，右键单击网格1并选择全部构建。 

研究1 

步骤1:频域

1在模型开发器窗口的研究1节点下，单击步骤1:频域。 

2在频域的设置窗口中，单击以展开研究扩展栏。 

3选中辅助扫描复选框。 

4单击+添加。 

5在表中输入以下设置:

range（a，b，c）函数表示以b为步长，在（a，b）区间内的求解模型，属于枚举函数的一种。设置越小步长的角数可以得到越精细的RCS图形，但相应的会线性延长求解时间和模型文件的大小。

6 在主屏幕工具栏中单击=计算

结果

选择 1

1 在模型开发器窗口中展开电场 (emw)节点。

2 右键单击表面并选择选择。

3 选择 “域” 3。

4 在电场 (emw) 工具栏中单击绘制。

此时的图像为入射角等于360度的总电场模。后一步是作者以30度为示范，通过改变参数值（phi（deg））来展示如何查阅和绘制不同入射角下的电场模图形。

电场 (emw)

1 在模型开发器窗口的结果节点下，单击电场 (emw)。

2 在二维绘图组的设置窗口中，定位到数据栏。

3 从参数值 (phi (deg)) 列表中选择 30。

4 在电场 (emw) 工具栏中单击绘制。

如此，我们得到的图示便是phi（入射角）为30度时的场。如果我们想要别的角度，便按照以上方法，更改参数值（phi（deg））到对应的值即可。

表面

1 在模型开发器窗口中，单击表面。

2 在表面的设置窗口中，单击表达式栏右上角的替换表达式。从菜单中选择组件 1 (comp1)> 电磁波，频域 > 电 > 电场 - V/m>emw.Ez - 电场， z 分量。

3 在电场 (emw) 工具栏中单击绘制。

4 回到表面节点，单击表达式栏右上角的替换表达式。从菜单中选择组件 1 (comp1)> 电磁波，频域 > 电 >相对电场 - V/m>emw.relEz - 相对电场， z 分量。

5 在电场 (emw) 工具栏中单击绘制。

此时的图形表示30度入射角下船身产生的相对电场的瞬时值。

由于相对场是观测场与背景场之差，因此在尾流中以及总场因反射而增强的区域，其大小都会增大。如果我们绘制相对场的绝对值，这种趋势会更加明显。

6 定位到表达式栏。在表达式文本框中键入 “abs(comp1.emw.relEz)”。

此表达式中的abs（）为绝对值，括号内的函数为我们刚刚设置的替换表达式。

7 在电场 (emw) 工具栏中单击绘制。

此时的图形表示30度入射角下船身产生的相对电场的大小。

二维远场 (emw)

1 在模型开发器窗口的结果节点下，单击二维远场 (emw)。

2 在极坐标绘图组的设置窗口中，定位到数据栏。

3 从参数选择 (phi) 列表中选择来自列表。

4 从参数值 (phi (deg)) 列表中选择 30。

辐射方向图 1

1 在模型开发器窗口中展开二维远场 (emw) 节点，然后单击辐射方向图 1。

2 在辐射方向图的设置窗口中，定位到计算栏。

3 找到角度子栏。在角数文本框中键入 “360”。

4 在二维远场 (emw) 工具栏中单击绘制。

该极坐标图是远场随角度的变化图，用于测定反射场如何向不同方向辐射。

极坐标绘图组 3

右键模型开发器窗口中的结果节点，选择极坐标绘图组。

全局 1

1 右键单击极坐标绘图组 3 并选择全局。

2 在全局的设置窗口中，定位到 r 轴数据栏。

3 在表中输入以下设置：

10*log10(at2(r0*cos(phi),r0*sin(phi),emw.bRCS2D))

4 定位到角数据栏。从参数列表中选择表达式

5 在表达式文本框中键入 “phi”。

6 单击以展开标题栏。从标题类型列表中选择手动。

7 在标题文本区中键入 “单位长度 RCS，单位为分贝”。

8 在极坐标绘图组 3 工具栏中单击绘制。