简介：近场仿真是天线设计与电磁兼容分析中的关键手段，其结果紧密依赖于人为定义的采样路径。与远场方向图不同，近场数据必须依附于采样线或采样面，仿真精度的上限往往不是求解器，而是采样线的位置、长度与点数编排。通过在设计阶段引入变量驱动几何建模，并正确理解波端口激励的功率归一化机制，可有效避免后处理中常见的3dB系统性误差。基于HFSS的RadField模块和VBScript脚本接口，工程师能够批量生成多条近场采样线，自动完成近场设置与数据导出，显著提升参数扫描与优化效率。该流程适用于800MHz至1.4GHz频段的近场分布评估、电路板辐射诊断及实测探头轨迹对比等场景。借助Python的matplotlib库对导出的CSV数据进行反归一化与三维重建，可快速直观地呈现近场分布规律，使仿真结果与实测数据精准对标。
1. 近场仿真最容易被低估的一步：采样线的布置
近场仿真与远场仿真最大的区别在于，远场结果是一组角度扫描曲线，而近场结果必须依附于人为定义的采样线或采样面。也就是说，仿真精度不只由网格和收敛准则决定，采样线的位置、方向、长度和点数的编排方式，会在后处理阶段直接改变近场分布曲线的形态。HFSS 的近场仿真把“场计算”和“场采样”分成两步：先在全波求解完成后，在一块自定义的线或面上做场强插值；再导出该路径上的电场或磁场分量。这就是为什么很多人在同一模型上得到完全不同的近场曲线——不是求解器算错了，而是采样对象没有按物理需求布置。
这份面向 800M-1400MHz 频段的 HFSS 近场仿真流程，完整覆盖从建模变量、波端口激励、求解设置、近场线布置到脚本化批量导出的全部环节。适合需要把近场分布作为评估指标的天线设计、EMC 分析和电路板辐射诊断场景。与常见的远场方向图评估不同，近场数据经常要和实测探头位置对应，因此采样线的物理坐标必须可复现、可参数化。
2. 先立变量再画模型：近场仿真结果可复现的前提
2.1 设计变量替换硬编码尺寸，后处理返工率直降
打开 HFSS 后不要直接画图，先进入HFSS > Design Properties添加设计变量。这份资料里的模型变量包括 DELX、DELY、h1、h2 这类几何参数，也有控制求解频段的变量。用变量替代数值的好处不只是方便调参，更重要的是脚本化之后，每条采样线的生成语句里可以直接引用变量名，避免在 VBS 脚本里再去改坐标常量。常见做法是先把所有驱动尺寸列成表，比如基板长度、宽度、介质层厚度、馈电位置偏移，再在模型里全部绑定到变量。这样当近场采样面需要从 3mm 抬高到 5mm 时，只需改 zpos 变量，无需重建模型。
变量命名也要按近场采样的习惯来。对于近场线扫描，我会额外定义 zpos（采样面高度）、xspan（X 向扫描长度）、yspan（Y 向扫描长度）、Xnum（X 向采样点数）、Ynum（Y 向采样点数）。这些变量不仅用于画采样线，还会写进 VBS 脚本的 InputBox 默认值中，减少每次重复输入的出错概率。
2.2 用变量驱动基础形状与复制镜像
画模型时，按变量名输入尺寸。例如创建一个矩形，Position 写-DELX/2, -DELY/2, h1，XSize 写DELX，YSize 写DELY。这样后续用Duplicate Mirror做对称复制时，镜像平面位置可以直接引用0mm或h2/2，不会出现镜像后几何错位的问题。HFSS 的建模历史树会保留所有参数，之后点击任意一个基础形状，都能看到完整的变量关系。
对称结构在近场仿真中尤其适合用镜像复制。假设一个差分结构，左侧是激励端口，右侧是被动端，先建好一半，再通过Duplicate Mirror沿全局坐标系 YZ 平面复制出另一半。复制后的对象名会自动带_1后缀，后续做Unite布尔合并时，要手动确认对象列表里包含了两部分几何体。
2.3 求解类型与求解频率对近场结果的影响
在HFSS > Solution Type中，这份资料选择了 Driven Modal 模式。近场仿真用 Driven Modal 还是 Driven Terminal，取决于端口定义方式。如果使用波端口（Wave Port）且关心的是模式电压和模式电流，Driven Modal 更自然；如果关心的是终端网络参数，用 Driven Terminal。近场线输出的场强值本身受求解类型影响不大，但激励功率的归一化基准会随端口定义变化，因此建议在仿真开始前就固定求解类型，中途切换会导致已布置的近场线上下文失效。
求解频率的设置与扫频策略需要先规划。资料中以 0.8GHz 作为单频点导出示例，而完整频段是 800MHz-1400MHz。求解频率可以设为 1.4GHz 以保证网格密度在高频端足够，然后做插值扫频。近场线导出的场分布在宽频带内变化剧烈，尤其是谐振点附近，所以扫频类型我一般用 Interpolating，并指定最大求解频点数为 151 个。这样既保证曲线平滑，又不至于让每个频点都做全波求解。
3. 波端口激励与功率归一化：1/2 因子的来源
3.1 波端口激励幅度与入射功率的换算
HFSS 在波端口激励中默认设置的幅度是峰值场强，不是 RMS 值。近场仿真导出的电场强度也是峰值场强。对于一个时谐场，电场瞬时值呈现正弦波形态，峰值场强 (E_{peak}) 与 RMS 值 (E_{rms}) 之间相差 (\sqrt{2})，但功率与电场平方成正比，所以用峰值场强计算功率时要多出 (1/2) 的因子。
这就是资料中专门提醒“因为 hfss 算出的场强为正弦波波峰，所以输入功率有一个 1/2 的因子”的原因。后处理时如果要用导出的近场电场值反算功率密度，必须按下式修正：
P_avg = 0.5 * Re(E x H*)如果直接在 HFSS 的 Field Overlays 里看场分布，看到的是峰值包络；但如果用 Report 导出近场数据再与实测结果（通常是 RMS 仪器读数）对比，就必须先除以 (\sqrt{2})。乘错这个因子，近场分布曲线的整体电平会差 3.01dB，这是近场数据对比中非常容易出现的系统误差。
3.2 Edit Sources 设置激励幅度与相位
在求解之前，进入HFSS > Fields > Edit Sources检查端口激励。对一个两端口模型，需要确认两个端口的幅度和相位设置：
| 端口 | 幅度 | 相位 | 说明 |
|---|---|---|---|
| Port1 | 1 W（归一化） | 0 deg | 主激励 |
| Port2 | 0 W | 0 deg | 匹配负载 |
这里要注意，HFSS 的波端口默认计算 1W 入射功率下的 S 参数，但近场结果中的场强实际是按端口归一化功率激励计算的。如果 Edit Sources 里把某个端口的幅度改成非归一化值，比如为了模拟实际输入功率 15kW，那么场强结果会按比例缩放。配合热搜里常见的“波导端口 15kW 功率施加”问题，此时直接修改 Edit Sources 中的幅度为对应功率，HFSS 会自动换算场强，不需要手动乘功率系数。
3.3 收敛判据与扫频设置的参数建议
在Analysis Setup中添加求解设置时，需要同时设置收敛条件和扫频。常用参数如下：
| 参数项 | 推荐值 | 说明 |
|---|---|---|
| Solution Frequency | 1.4GHz | 以最高工作频点网格加密 |
| Maximum Number of Passes | 10 | 防止无穷迭代 |
| Maximum Delta S | 0.02 | 相邻两次求解 S 参数变化量 |
| Sweep Type | Interpolating | 宽频带扫频 |
| Frequency Range | 0.8GHz - 1.4GHz | 覆盖完整工作频段 |
| Count of Solution Points | 151 | 频点密度 |
求解完成后，先看收敛曲线再布置近场线。如果 Delta S 在最后一次迭代中仍大于 0.02，说明网格未收敛，此时导出的近场场强在局部区域会出现明显跳变。不要急着画采样线，应检查是否存在窄缝、薄层等需要加密网格的几何结构，或者手动添加 Length Based Mesh Operation。
4. 近场采样线的布置：位置、长度与点数
4.1 为什么近场用线采样而不是面采样
近场面采样适合观察大范围内的场分布，但数据量巨大，且后处理时三维曲面反而掩盖了沿某一方向的细节。实测近场探头通常也是沿一条直线扫描，所以仿真中布置近场线与实测轨�迹对应，才能做直接对比。
在 HFSS 里，近场线是一根无厚度、无材质的 Polyline，它的作用不是参与电磁求解，而是为后处理提供一条插值路径。HFSS 会在这条线上按设定的 NumPts 均匀取点，计算该点的电场和磁场矢量。因此，画线时不需指定材料，也不能把它和任何物体做布尔运算，否则 HFSS 会把它当作三维实体参与网格剖分。
4.2 画线参数与坐标设定
画近场线用Modeler > Polyline。以 X 方向扫描线为例，在 3D Modeler 里将两个端点设为：
- 起点：
-xspan/2, y, zpos - 终点：
xspan/2, y, zpos
其中y是这条线的 Y 向偏移量。对于一组近场线，需要在这个 Y 方向上以dy为间距，逐条生成采样线。
线画好之后，在工程树中选中该线，右键选择Assign Excitation >，但实际上近场线不是激励，而是通过HFSS > Fields > Insert Near Field Setup > Line来绑定。此时需要指定线名、采样点数和归一化距离参数。
4.3 近场线“一一对应”的坑
资料里特别提到“近场线要一一对应”，指的是：在建模器里画的线名（比如lin3mx1、lin3mx2）必须与近场线设置及后续报表中的 Context 一一对应。如果在建模树里把线重命名，而 Near Field Setup 中未同步更新，导出数据时 HFSS 会报“Context not found”。
整个流程涉及三个层面的命名匹配：
- 建模器中 Polyline 对象的 Name
- Insert Near Field Setup 中 Line 参数引用的 Name
- Report 中 Trace 的 Context 参数
三者必须完全一致。因此资料中的 VBS 脚本统一用index = "lin3mx" & cstr(n)的规则生成线名，后续 CreateReport 和 AddTraces 也引用同一变量，这就是规避一一对应错误的标准做法。
5. 用 VBScript 批量生成近场采样线并绑定设置
5.1 HFSS 脚本接口的对象层级
HFSS 的 VBScript 接口以AnsoftHfss.HfssScriptInterface为入口，逐级获取桌面、工程、设计和编辑器。各层级对象关系如下：
Set oAnsoftApp = CreateObject("AnsoftHfss.HfssScriptInterface") Set oDesktop = oAnsoftApp.GetAppDesktop() Set oProject = oDesktop.GetActiveProject() Set oDesign = oProject.GetActiveDesign() Set oEditor = oDesign.SetActiveEditor("3D Modeler") Set oModule = oDesign.GetModule("RadField")对象层级中，oEditor控制几何建模，oModule控制辐射场和近场设置。注意，近场设置模块是RadField，不是Boundary，这一点容易搞混。GetModule("RadField")之后才能调用InsertNearFieldLineSetup。
5.2 批量生成 X 方向近场线
以高度 zpos 处的 X 方向扫描为例，先通过 InputBox 交互输入采样参数，再用循环批量生成多段线：
Dim zpos, Xnum, Ynum, xspan, yspan, dx, dy Dim y() Units = "mm" zpos = 3200 xspan = 4000 yspan = 2100 Xnum = 100 Ynum = 15 Redim y(Ynum) dx = xspan / (Xnum - 1) dy = yspan / (Ynum - 1) For n = 1 To Ynum y(n) = -yspan / 2 + (n - 1) * dy ZValue = cstr(zpos) & Units YValue = cstr(y(n)) & Units XValuemin = cstr(-xspan / 2) & Units XValuemmax = cstr(xspan / 2) & Units index = "lin3mx" & cstr(n) oEditor.CreatePolyline Array( _ "NAME:PolylineParameters", _ "IsPolylineCovered:=", true, _ "IsPolylineClosed:=", false, _ Array("NAME:PolylinePoints", _ Array("NAME:PLPoint", "X:=", XValuemin, "Y:=", YValue, "Z:=", ZValue), _ Array("NAME:PLPoint", "X:=", XValuemmax, "Y:=", YValue, "Z:=", ZValue)), _ Array("NAME:PolylineSegments", _ Array("NAME:PLSegment", "SegmentType:=", "Line", _ "StartIndex:=", 0, "NoOfPoints:=", 2)), _ Array("NAME:PolylineXSection", "XSectionType:=", "None"), _ Array("NAME:Attributes", "Name:=", index, "Flags:=", "", _ "Color:=", "(132 132 193)", "Transparency:=", 0, _ "PartCoordinateSystem:=", "Global", _ "MaterialValue:=", Chr(34) & "vacuum" & Chr(34), _ "SolveInside:=", true)) Next这段脚本的关键参数有三个。IsPolylineCovered必须设为true，否则 HFSS 会把线视为开放曲线，禁止用于近场采样。StartIndex与NoOfPoints是分段定义时的起止索引，单段直线用StartIndex=0, NoOfPoints=2即可。MaterialValue中的vacuum需要带英文双引号，所以用Chr(34)拼接，直接写"vacuum"会被 VBScript 解析成字符串字面量而丢失引号。
5.3 逐线绑定近场设置
生成线之后，逐条调用InsertNearFieldLineSetup：
For n = 1 To Ynum index = "lin3mx" & cstr(n) namepara = "NAME:" & index oModule.InsertNearFieldLineSetup Array( _ namepara, _ "UseCustomRadiationSurface:=", false, _ "Line:=", index, _ "NumPts:=", cstr(Xnum)) Next这里NumPts必须与线长度和采样间距匹配。资料中 X 方向扫描 4000mm、取 100 个点，对应的空间采样间距约为 40.4mm，在 0.8GHz 时波长约 375mm，采样间距约波长的十分之一，能够保留近场分布的主要空间频率。如果采样间距大于四分之一波长，近场曲线的空间波动会被严重混叠，看起来像噪声。
选择合适的采样点数时，一个快速经验公式是：
NumPts >= (线长 / (0.1 * 最短波长)) + 1这样可以保证在最高频率 1.4GHz 下，采样间距仍优于波长的十分之一。
5.4 用 CreateReport 批量导出近场曲线
近场数据导出走ReportSetup模块，而不是RadField：
Set oModule = oDesign.GetModule("ReportSetup") oModule.CreateReport "XY Plot 27", "Near Fields", "Rectangular Plot", _ "Setup1 : Sweep", _ Array("Context:=", "lin3mx1"), _ Array("NormalizedDistance:=", Array("All"), _ Array("NAME:VariableValues", "Freq:=", "0.8GHz")), _ Array("X Component:=", "NormalizedDistance", _ "Y Component:=", Array("NearETotal")), _ Array()Context参数指定了读取哪条近场线，这就是上一章强调的一一对应关系。NormalizedDistance表示横坐标是归一化距离，从 0 到 1 对应线的起点到终点。Freq是当前导出的频点。若要导出全部 15 条 X 方向线的数据，需要对每条线重复调用 CreateReport，或者用循环生成不同名称的报表。
导出报表后，右键选择Export To File，数据会保存为 CSV 格式。CSV 中列出的横坐标是归一化距离，单位是归一化值，不是毫米。做三维曲面图前，必须把横坐标乘以线长，换回物理坐标。
6. 从归一化距离到三维曲面：近场后处理的三个细节
6.1 导出 CSV 的列结构与数据对齐
CSV 导出时，每一行对应一个采样点，列包括NormalizedDistance、NearETotal，以及需勾选的场分量（如NearEX、NearEY、NearEZ、NearHTotal等）。多条线分别导出为多个 CSV 文件，文件命名建议与线名保持一致。手动操作时，容易把 15 条线的数据混在一起，所以我的做法是在导出前先选中对应报表，再右键导出，避免选错。
6.2 横坐标反归一化与坐标轴重建
反归一化的公式很简单：
x_mm = NormalizedDistance * 线长 + 线起点坐标比如第 5 条 X 线，起点是-2000mm，线长4000mm，则 CSV 里NormalizedDistance=0.25对应的物理坐标是-2000 + 0.25 * 4000 = -1000mm。在 Excel 里新增一列做换算：
C2 = A2 * 4000 - 2000这里A2是归一化距离，4000是 X 跨度，-2000是起点 X 坐标。如果扫描线是 Y 方向，就换成对应的 Y 起点与跨度。反归一化之后，15 条线的数据才能以(x, y, |E|)三元组组织成曲面数据。
6.3 用 Python 快速重建三维近场图
数据量大时，Excel 曲面图操作繁琐且容易卡顿，我一般直接用 Python 的 matplotlib 处理。以下脚本读取按线名排列的 CSV 文件，重建二维网格并绘制三维曲面：
import numpy as np import pandas as pd import matplotlib.pyplot as plt from mpl_toolkits.mplot3d import Axes3D lines = 15 xspan = 4000.0 yspan = 2100.0 X, Y, E = [], [], [] for n in range(1, lines + 1): df = pd.read_csv(f"lin3mx{n}.csv") x = df["NormalizedDistance"].to_numpy() * xspan - xspan / 2 y = np.full_like(x, -yspan / 2 + (n - 1) * yspan / (lines - 1)) e = df["NearETotal"].to_numpy() X.extend(x) Y.extend(y) E.extend(e) X = np.array(X).reshape(lines, -1) Y = np.array(Y).reshape(lines, -1) E = np.array(E).reshape(lines, -1) fig = plt.figure(figsize=(10, 6)) ax = fig.add_subplot(111, projection="3d") ax.plot_surface(X, Y, E, cmap="jet", edgecolor="none") ax.set_xlabel("X (mm)") ax.set_ylabel("Y (mm)") ax.set_zlabel("|E| (V/m)") plt.savefig("nearfield_3d.png", dpi=300)plot_surface要求输入二维网格，所以脚本中先用lines和采样点数量做reshape。如果某条线的 CSV 采样点数不一致，可以先用df.shape[0]动态获取点数再做reshape。这个脚本可以直接复用于任意高度 zpos 的近场切片，只要把xspan、yspan和lines换成实际参数。近场数据与实测轨迹对齐时，同样用反归一化后的物理坐标去关联探头位置，不再依赖 HFSS 里的归一化坐标。
本文还有配套的精品资源，点击获取
