COMSOL偶极子天线散射体仿真：远场方向图畸变与参数分析

做射频仿真的人，十有八九会遇到这么个情况：天线单独仿真时指标漂亮得不行，一放进整机里就立马翻脸。增益掉了、方向图歪了、谐振点飘了，怎么看怎么别扭。最典型的困扰之一，就是天线旁边多了个金属散射体之后，原本对称的偶极子远场方向图变得不再对称，主瓣被拉扯，后瓣忽大忽小。这篇内容我拿 COMSOL Multiphysics 里一个带散射体的偶极子天线模型来完整拆一遍，重点放在远场图的物理含义、二维方向图的提取方法，以及散射体距离变化时方向图怎么跟着变。适合正在用 COMSOL RF 模块做天线仿真、但又被远场后处理坑过的人，尤其建议刚开始接触散射体分析的同行先通读一遍。

1. 模型设计思路：为什么要在偶极子旁边放散射体

1.1 基础模型选型与关键尺寸确认

我用的模型是 2.45 GHz 的半波偶极子天线，工作波长 λ 约为 122.45 mm。两臂各取 0.24λ，总长度约 58.8 mm，半径取 0.005λ，约 0.61 mm。不含散射体时，这个偶极子的谐振点是比较准的，S11 在 2.45 GHz 附近能压到 -20 dB 以下。

选择半波偶极子有两个理由。第一，它的远场方向图理论上非常清楚，在自由空间里是标准的面包圈形，E 面方向图是“8”字形，H 面方向图是圆形，拿来做散射体影响分析时，任何方向图畸变都能直接归因到散射体上，不会掺杂天线本身的设计问题。第二，偶极子对周围的金属物体非常敏感，散射体稍微靠近一点，馈电点的输入阻抗就开始剧烈变化，这正好方便我们观察谐振偏移和匹配恶化。

散射体我用了两种几何做了对比，一种是矩形金属贴片，尺寸 0.5λ × 0.25λ，厚度设为 1 mm，另一种是直径 0.2λ 的金属圆柱。实际跑下来，矩形贴片产生的方向图畸变更明显，因为它的迎风面积更大，二次辐射更强，本文后面讨论的仿真结果都以矩形贴片为准。

1.2 求解域和边界条件的搭配：SBC 还是 PML

求解域我是用一个空气球包住天线和散射体，这个球体的半径至少取 0.5λ，大约 61 mm。外层边界有两种常见选择：散射边界条件（SBC）和完美匹配层（PML）。

如果你刚接触 COMSOL，先记住一个结论：只关心天线在自由空间的辐射方向图和 S11 时，SBC 比 PML 更省事。SBC 理论上能把入射到边界上的波近似吸收掉，反射系数在正入射时比较低，但在大角度斜入射时会变差。偶极子方向图比较全向，边界不同位置入射角差异很大，SBC 会引起一定的非物理反射，尤其远场计算包含低电平后瓣时，这个反射会污染结果。

PML 的吸收性能更好，对入射角不敏感，但需要额外画一层厚度约为 λ/4 到 λ/2 的层域，网格数量上去之后计算量明显增加。我的经验是：做散射体参数扫描时要跑很多组工况，用 SBC 把单次求解时间压下来；最后确定两三个关键工况时，再换成 PML 做核算。千万不要一开始就无脑上 PML，等参数扫描跑完可能天都亮了。

关于远场计算有一个非常重要的边界设置原则：COMSOL 里 Far-Field Domain 节点必须和吸收边界配合使用，远场计算是基于边界上的等效电流和磁流积分完成的。如果你选 SBC，就把 SBC 所在的球面设成远场积分边界；如果你选 PML，PML 层外面不需要再算远场，通常是把空气球和 PML 的交界面作为远场积分面。这里最容易犯的错误是，把远场积分面选在 PML 内部。PML 内部是坐标拉伸后的非物理区域，场分布已经做过变换，从那里积分得到的远场方向图完全不对，而且多数情况下图形形状诡异、数值忽大忽小，一眼看不出来哪里错了。

2. 散射体影响方向图的物理机理拆解

2.1 感应电流与二次辐射：方向图畸变的根源

金属散射体在偶极子产生的入射电磁场中，表面会被感应出传导电流，这个感应电流本身又向外辐射电磁波，相当于一个次生辐射源。偶极子辐射的直接波，和散射体感应电流产生的二次辐射波，在空间任意观察点叠加，形成干涉。两者之间的相位差主要由路径差决定，路径差又由散射体相对偶极子的方位和距离决定。

用等效二源干涉模型来理解就非常直观。把偶极子看成源点 O，散射体等效成一个位于 S 点的次级源。如果观察点 P 在空间某个方向，直接波从 O 到 P 的路径长度，和二次辐射从 O 到 S 再到 P 的路径长度不一样，两者的相位差大约是：

Δφ ≈ (2π/λ) × (OS + SP - OP)

这个相位差随风向角变化，因此不同方向上的叠加结果有的增强、有的减弱，方向图就不再是原来干净的“8”字形。更为关键的是，散射体通常不是各向同性的点源，它的感应电流分布和自身形状有关，矩形贴片还会激发出明显的方向性辐射，所以实际方向图的不对称程度比点源模型预测得还要复杂。

我在仿真里做过一组对比：当矩形贴片放置在偶极子正侧方、距离很近时，E 面方向图在贴片那一侧出现了很大的凹陷，并且另一侧冒出了新的瓣。这就是典型的二次辐射干涉结果——某些方向上直接波和二次辐射相位接近同相，场强增强；另一些方向接近反相，场强被抵消。

2.2 用参数扫描量化散射体的几何影响

为了把散射体的影响说清楚，我用 COMSOL 的参数化扫描功能，把矩形贴片到偶极子轴线的距离 d 设成扫描参数，分别计算了 d = 0.25λ、0.5λ、0.75λ、1.0λ 四种情况，同时保留一组无散射体结果作为基准。

扫参之前需要把几何里的距离参数和网格设置绑定好。这里有个技巧，空气域尺寸不要跟着 d 变，固定不变，只移动散射体位置。否则每次改动几何时网格会重新划分，不同工况之间的数值差异就掺杂了网格不一致带来的误差，方向图的细微变化就很难归因到物理机制上。

仿真下来，几个典型结果如下表所示。注意这些数值是我这组尺寸参数下的具体结果，放到别的模型里绝对值会变，但趋势是一致的。

可以清楚看到，散射体靠近后，最大增益逐渐下降，后向电平升高，谐振频率向低频偏移，S11 最差点也在变差。原因是散射体和偶极子之间存在互耦，相当于在天线输入端额外引入了一个与频率相关的阻抗，改变了天线的输入阻抗特性，所以谐振点产生偏移。

这个表格里的“后向电平”是指 180° 方向上的增益相对主瓣最大增益的差值。无散射体时偶极子的 H 面全向性很好，E 面是双向对称的，“后向”其实就是主瓣的反方向。加上散射体后，主瓣方向本身还会偏移，所以对比后向电平时一定要先重新找主瓣方向，否则会得到一个没什么意义的数据。

3. 远场方向图的计算概念与二维绘制实战

3.1 远场节点设置中容易被忽略的三个点

在 COMSOL 中，远场计算不是直接用求解器算出来的，而是把近场结果通过 Stratton-Chu 积分公式变换得来的，也就是等效源法。求解器在空气域内部算出每个网格点上的电磁场，远场后处理把这些近场值在预先定义的外边界面上积分，推出无穷远处的辐射场。理解这一点很重要，因为它决定了我们为什么要刻意保证积分边界的场分布准确性。

实际建模时，我用的是“电磁波，频域”物理场接口，在模型树下依次添加：

1. 一个Lumped Port边界条件，设置在偶极子两臂之间的馈电缝隙上，阻抗设为 50 Ω。
2. 一个Far-Field Domain节点，作用域选为最外层的吸收边界或空气域外表面。
3. 一个Far-Field Calculation的派生值，用于提取特定方向上的远场值。

Far-Field Domain 节点有三个参数要特别留意：

• 电场分量参考方向：决定远场结果里 Ex、Ey、Ez 的坐标系，通常用球坐标的 theta 和 phi 方向来做方向图，这样最直观。
• 远场积分边界的选择：必须是一个封闭面，不能漏选面，也不能选中内部金属表面。金属表面属于边界条件域，不算远场积分的一部分。
• 表达式中的单位：COMSOL 默认远场电场单位是 V/m，增益单位是 dBi。后处理画增益方向图时，要选对表达式，我习惯直接用Gain或Gain dBi，避免自己手动算对数，省得单位换算出错。

提示：如果你同时选中了 SBC 外表面和内部某个介质分界面作为远场域，COMSOL 会做多次积分然后叠加，结果几乎是必错的。选的时候只保留最外层的那一个封闭面。

3.2 二维方向图的绘制方法：从三维方向图到截面

COMSOL 里画二维方向图，入口通常在“派生值”菜单下的“二维远场图”。这里的关键是，你要清楚自己画的是哪个截面的方向图。

偶极子天线有两个标准截面：

• E 面：包含天线电流方向的平面。对沿 z 轴放置的偶极子来说，E 面就是 xz 平面或 yz 平面，方向图呈 8 字形。
• H 面：垂直于天线电流方向的平面，也就是 xy 平面，方向图是圆形。

在三维远场图的基础上，二维方向图等于把观察球面上某个大圆上的增益值单独提出来画成极坐标或直角坐标曲线。比如我要看 E 面方向图，就在二维远场图设置里把 phi 固定为 0°，theta 从 0° 扫到 360°，横轴是 theta，纵轴是 gain；或者反过来，固定 theta 看 phi 变化，那就是 H 面方向图。

我经常看到新手在这里犯迷糊：用二维远场图里的“极坐标”选项，结果画出来一个上下不对称、但感觉哪里不对的图，原因往往是 theta 的起始角度设错。COMSOL 默认 theta 从 0° 到 180°，对于偶极子的 E 面，0° 到 180° 只能画出一半方向图（上半平面），另一半需要把 theta 范围设到 360°，或者在设置里勾选“对称”相关的选项。如果你只画了 0 到 180°，另一半完全缺失，看起来像是主瓣被砍掉了一大半。

下面是一个典型的二维 E 面方向图绘制步骤，基于我自己的模型：

1. 在结果节点下新建“二维远场图”。
2. 把“数据”来源选成“来自父节点”或直接选中你的远场计算节点。
3. 设置坐标类型为“极坐标”。
4. theta 起点设 0°，终点设 360°，步长可以先用 1° 或 0.5°，看曲线平滑度再调。
5. 如果只看某个固定 phi 截面，在“选择”栏里指定 phi 角，比如 0° 对应 xz 平面，90° 对应 yz 平面。
6. 横轴（方向轴）选 theta，纵轴选Gain dBi或abs(emw.Efar)/max(abs(emw.Efar))。

3.3 三维远场图和二维截面的对应关系

三维远场方向图看起来直观，但你需要从三维图里读出具体方向上的数值时，就会觉得它远不如二维方向图实用。三维方向图是一张彩色球面，增益大小用颜色和 z 轴高度同时表示；二维方向图则是把某个截面的增益画成一条曲线，可以直接读出半功率波束宽度、前后比、副瓣电平这些量化指标。

我在做完散射体参数扫描后，通常是先把每组的二维 E 面方向图叠加到同一个坐标系里，纵轴增益设为相对值（归一化到自身最大值）。这样能非常清楚地看到散射体在哪个角度范围内把主瓣压低了，在哪个角度范围产生新瓣，甚至能看到零陷位置的移动规律，比逐张翻三维图高效得多。

COMSOL 里把多条曲线放到同一张图里，方法是在同一个“二维绘图组”节点下新建多个“二维远场图”子节点，每个子节点对应一个参数扫描解。如果你的扫描变量名是d_scat，在子节点的“参数选择”里分别勾选不同的扫描参数值即可。曲线颜色记得手动区分，不然几根线混在一起很难分辨。

4. 实操过程中最常见的四个问题与排查思路

4.1 远场方向图不光滑，出现毛刺

这个现象在低增益区域尤其明显。方向图曲线在 -20 dBi 以下的部分往往会有锯齿状抖动，看起来像收敛了，又像没收敛。

排查顺序是：

1. 先看网格。远场积分面对应区域的最大单元尺寸要控制在 λ/8 以内，理想情况下靠近天线的局部网格做到 λ/12 甚至更细。不要只盯着天线附近的网格，散射体周围和远场积分面附近的网格同样重要。因为散射体的感应电流分布直接决定了二次辐射的准确性，如果散射体表面网格太粗，表面的感应电流会有数值震荡，表现到远场上就是方向图局部跳动。
2. 再看边界条件。如果你用的是 SBC，尝试把空气域半径加大到 0.75λ 或 1λ，看看毛刺是不是明显降低。如果是 SBC 反射造成的，加大边界距离通常能改善。
3. 检查是否选中了非封闭的远场域。有时候几何布尔操作会在内部留下一些碎片面，被自动选进 Far-Field Domain，这种碎片面的积分贡献是个错误源。

4.2 S11 仿真值和手算谐振长度对不上

偶极子臂长取 0.48λ 时，理想情况下谐振频率应该很接近 2.45 GHz，但加上散射体后谐振点经常往低频偏移。这个不是仿真 bug，而是互耦带来的必然结果。

散射体靠近时，相当于在天线附近引入了一个寄生单元，它让天线的等效电长度变长，所以谐振频率向下漂。如果你需要让加上散射体后的天线仍然谐振在目标频率，最直接的办法是在参数扫描里把偶极子臂长也设成变量，做一个二维扫描：一个变量是散射体距离，一个变量是臂长，然后找出每个散射体距离下最优的臂长修正量。

需要注意的是，集总端口（Lumped Port）的积分线方向不能搞反。偶极子两臂之间的电压方向决定了电流方向，如果积分线方向反了，端口阻抗的虚部会非常奇怪，S11 曲线直接不对。检查方法是看电场矢量箭头方向是否从一臂指向另一臂，如果反了就翻转积分线方向。

4.3 参数扫描计算量太大，一个晚上跑不完

散射体参数扫描本质上是多工况求解，每个工况都要重新划分网格（因为几何变了），这是最主要的耗时来源。如果我用的是三维全波模型，每个工况网格量容易超过 200 万自由度，在普通工作站上单工况求解可能要十几分钟，扫十个点就要两三个小时。

我常用的提速办法是按顺序来：

• 先用二维轴对称或 2D 模型做预扫描。如果天线系统绕某个轴近似旋转对称，散射体也在这个对称轴上，可以考虑二维模型。但散射体在侧方时破坏对称性，二维模型用不了，只能全波三维。
• 网格策略上，把扫描范围内的几何变化区单独设置一个“扫掠”网格。散射体移动时只重画它附近的活动区域，其他区域网格保持不变。COMSOL 的自动重新划分网格可能会导致整体重画，但手动控制几何序列里的网格域可以部分避免。
• 设置好求解器的自适应频率。如果扫描频点较多，先用粗频率步长（比如 100 MHz）找到谐振点的大致位置，再在谐振点附近加密频点，不必全频段都跑细扫。
• 合理利用对称面。如果没有散射体，偶极子本身关于 xz 和 yz 面对称，可以只建四分之一模型。但加散射体后对称性被破坏，这招基本作废。如果散射体是对称放置的，比如在偶极子两侧各放一块完全相同的贴片，那么对称面和偶极子本身的对称面重合，仍然可以减半建模。

4.4 二维方向图主瓣方向漂移导致的误判

散射体靠近后，主瓣方向不再老老实实待在 90° 方向，而是会向远离散射体的那一侧偏转。我在 d = 0.25λ 那组数据里，主瓣方向从 90° 偏到了大约 102°，偏移量不小。所以做前后比或者副瓣电平时，不能直接拿 90° 作为主瓣方向去查值，否则会把主瓣边缘当成最大值，所有方向性指标全部算错。

正确做法是在后处理里用“最大增益”探针先找到主瓣方向，再用该方向做归一化。如果散射体很强，还可能出现主瓣分裂成两个瓣的情况，这时候只报最大增益会掩盖方向图畸变的真相。建议除了最大增益之外，再加一个“主瓣方向附近 ±20° 内增益波动”的指标，用来评判方向图恶化程度。

5. 几个可以让结果更接近实际工程的经验补丁

5.1 散射体的材料属性不要只盯着理想导体

电磁仿真里最常见的一个陷阱是，默认金属就是完美导体（PEC），但工程中的散热片、支架、结构件大多是铝合金，甚至是镀镍钢，它们的电导率有限。对于 2.45 GHz 这个频段，铝的电导率大约是 3.56 × 10^7 S/m，趋肤深度不到 1.6 μm，所以损耗不算大，在很多场景下用 PEC 近似误差可以接受。

但如果散射体表面做了导电氧化处理，或者表面粗糙度比较大，等效表面电阻会明显升高，感应电流会减弱，导致二次辐射强度下降，方向图畸变程度也会轻微减弱。追求精度时可以把散射体设置成阻抗边界条件，输入实际电导率和相对磁导率。我用过某型号铝合金材料参数仿完之后，Gain 的变化量也就是零点几个 dBi，方向图形状几乎不变，所以先判断你的应用对精度要求多高，再决定要不要舍简求繁。

5.2 加散射体后不要忘记重新优化匹配

很多人做天线加结构件仿真，盯着方向图看半天，却忘了看馈电端匹配。我的经验是：散射体不仅影响方向图，对 S11 的影响往往比方向图更致命。d = 0.25λ 时，我的偶极子原本在 2.45 GHz 匹配得很好，加上贴片后 S11 最差值只剩 -9.8 dB，也就是有大约 10% 以上的功率被反射掉，实际辐射效率自然跟着下降。

这种情况下，单纯调整天线臂长只能把谐振点拉回来，但匹配深度不一定够。比较实用的方法是加一个简单的 Γ 型匹配网络，或者在 COMSOL 里做一个集总元件参数扫描，把并联电感和串联电容的值扫出来，找到让 S11 恢复到 -15 dB 以下的组合。这个优化过程和天线单独设计时差别不大，但要把散射体带来的互耦效应一起纳入优化目标。

5.3 二维方向图与三维远场结果交叉验证

单看二维方向图容易漏掉一些信息。我习惯用三维远场图做一个快速交叉验证：先看三维图的整体形状，观察有没有在非 E 面/H 面的方向上出现异常大瓣。如果三维图上某个斜向方向出现了一个二维截面完全没覆盖到的强瓣，说明散射体激起了较强的交叉极化辐射。这时候必须额外画一个 co-pol 和 cross-pol 的远场方向图对比，确认极化纯度是否满足系统指标。COMSOL 里可以在远场表达式中直接写出emw.Efarphi和emw.Efartheta，分别对应两个正交极化分量，画在同一张图里就能清楚看到交叉极化电平。

在我仿的那组带矩形贴片的模型里，d = 0.5λ 工况下交叉极化分量大概比主极化低 18 dB 左右，不算严重，但方向图副瓣区域交叉极化明显抬升，会对某些对极化纯度敏感的系统造成影响。如果你做的是圆极化天线或者分集天线，这一步绝对不能省。

5.4 导出数据做二次分析的小技巧

COMSOL 的绘图节点可以直接导出文本格式的曲线数据。我通常把每一条二维方向图曲线导成 CSV，然后在外部工具里做比较和后处理。导出时注意两点：第一，角度分辨率要足够高，建议 theta 步长取 0.5° 以下，否则后期算半功率波束宽度时插值误差会很大；第二，导出的增益数值本身已经是 dBi，不需要再手动取对数，但你在外部工具里做归一化时建议用线性值归一化后再转回 dB，而不是直接对 dB 值做减法，两者在数值上有细微差异，容易导致曲线之间出现奇怪的整体平移。

用 MATLAB 或者 Python 的 matplotlib 都能很方便地把多条 CSV 曲线画在一起。我自己经常用 Python 里的小脚本做批量处理，把无散射体和不同 d 值的方向图曲线自动叠加，标注主瓣方向和前后比，一步到位，节省大量重复操作。

6. 把散射体效应当设计工具，而不是单纯的问题源

最后分享一个我个人在多次实操之后转变的思路。散射体并不总是坏事，熟练的射频工程师反而会利用金属支架或者地板来改善特定方向的覆盖。比如在某些手持设备里，天线旁边的金属边框会对远场方向图产生整形作用，如果控制好位置和尺寸，可以压低向人体方向的辐射，同时增强朝向自由空间一侧的增益。这本质上和我上面仿真的散射体效应是同一个物理过程，只是目标从“抑制畸变”变成了“定向整形”。

所以每次做带散射体的仿真，我都会把参数扫描当成一次方向图整形能力探索，先看辐射机理，再谈指标。这样跑仿真时的心态会完全不同，不焦虑，反而能发现很多设计上的可能性。如果你也正在处理天线被周围结构器件“带偏”的问题，建议按本文的模型搭一遍，亲手看看不同散射距离下方向图的演化过程，你会对远场图有更直觉的理解。