4、参数扫描收敛失败的表现形式：不收敛、伪收敛、振荡、发散、结果跳变、仿真中断等典型现象
在天线设计与电磁仿真中，参数扫描（Parametric Sweep）是评估设计余量、优化几何尺寸、分析频带响应的核心手段。然而，当扫描变量（如贴片长度、馈线宽度、介质厚度）跨越某些临界值时，求解器可能无法在预设迭代次数内达到稳定解，或即使“看似收敛”却给出物理上不可信的S参数。本章系统梳理六种典型失效现象——不收敛、伪收敛、振荡、发散、结果跳变、仿真中断，从数值机理、外在表现、判断路径到修复策略进行全链路拆解，帮助工程师在仿真排障中快速定位根因。
4.1 不收敛：迭代无法达到预设精度
核心原理：不收敛是指求解器在达到最大迭代次数（如HFSS中的Maximum Number of Passes）或能量误差阈值（Delta S）之前，无法使S参数变化量降至设定值以下。对于有限元法（FEM），每次自适应网格加密后，端口S参数应逐渐趋于稳定；若网格加密后S参数仍剧烈变化，或残差曲线呈水平趋势，则判定为不收敛。常见诱因包括：模型存在极薄介质层导致网格长宽比过大、端口尺寸与波长不匹配、材料损耗正切设置异常、或求解频率设置过高导致电尺寸过大。
判断方法：观察收敛曲线（Convergence Curve）和网格统计信息。若Delta S曲线在多次迭代后仍高于0.02，且网格数量已超过内存限制，则需人工干预。具体操作流程如下：
- 检查求解设置中的最大迭代次数是否过小（如小于6次）；
- 查看每次迭代的Delta S值，若连续3次迭代变化量均大于0.05，则判定为不收敛；
- 检查模型是否存在小于λ/1000的几何特征（如薄铜箔、细缝）；
- 尝试将求解频率降低至中心频率的1.2倍，观察收敛性是否改善。
某2.4GHz贴片天线，介质基板厚度0.254mm（Rogers RO4350B），初始仿真设置最大迭代10次，Delta S=0.02。运行至第8次迭代时，Delta S仍为0.038，且网格数量已达120万。分析发现，基板厚度与贴片宽度之比为1:150，导致Z方向网格过度细化。修复措施：将基板厚度方向网格层数手动设为3层，并启用“Thin Sheet”近似处理铜箔，最终在第6次迭代达到Delta S=0.015，收敛成功。
4.2 伪收敛：结果看似稳定但物理失真
核心原理：伪收敛是指求解器在数值上达到收敛判据（如Delta S<0.02），但结果与理论值或实测值严重偏离。其本质是网格剖分未能捕捉到关键电磁场分布——例如天线馈电点附近的电流奇异性、介质基板中的表面波模式、或辐射边界附近的倏逝波。伪收敛比不收敛更危险，因为它会给出“可信”的S参数曲线，误导后续优化决策。
判断方法：不能仅依赖收敛曲线，需结合场图验证和能量守恒检查。具体步骤：
- 查看天线方向图是否出现非物理的栅瓣或零陷；
- 检查端口反射系数S11是否在谐振频率处出现“假谐振”（即带宽极窄的尖峰）；
- 计算辐射效率，若超过100%或低于50%（对无耗介质），则存在伪收敛；
- 对比不同网格密度下的结果，若S11曲线差异大于0.5dB，则判定为伪收敛。
| 检查项 | 正常范围 | 伪收敛特征 | 验证工具 |
|---|---|---|---|
| S11谐振深度 | -15~-30dB | <-40dB且带宽极窄 | Smith圆图 |
| 辐射效率 | 70%~95% | >100%或<50% | 场计算器 |
| 方向图对称性 | 主瓣对称 | 出现非对称旁瓣 | 3D方向图 |
| 网格收敛性 | 网格加密后S11变化<0.1dB | 变化>0.5dB | 参数扫描对比 |
4.3 振荡与发散：迭代过程中的数值失稳
核心原理：振荡表现为Delta S曲线在迭代过程中反复波动（如0.02→0.05→0.03→0.06），无法单调下降；发散则表现为Delta S持续增大或出现NaN/Inf值。两者均源于数值不稳定，常见于：模型包含高Q值谐振结构（如介质谐振器）、端口激励与模式不匹配、或材料参数存在负值（如错误设置负介电常数）。振荡的物理本质是网格加密后，局部场强变化导致S参数在相邻迭代间跳变；发散则通常意味着矩阵求解失败，可能是由于网格严重畸变或边界条件冲突。
判断与修复流程：
- 若Delta S曲线呈周期性波动，检查是否启用了“Interpolating Sweep”且频率点过密（建议每频段不超过201点）；
- 若出现发散，首先检查模型是否存在重叠面或未连接的几何体（使用Modeler→Boolean→Unite检查）；
- 尝试将求解器从“Direct”切换为“Iterative”，并设置合适的收敛容差（如1e-6）；
- 对于高Q结构，启用“Coupled”模式或增加“Radiation Boundary”到模型边界的距离（至少λ/4）。
某介质谐振器天线（DRA）工作在5.8GHz，介电常数εr=38，仿真时在第3次迭代后Delta S突变为NaN。检查发现，DRA与微带馈线之间的空气缝隙宽度为0.01mm，导致网格剖分时产生负雅可比矩阵。修复措施：将缝隙宽度改为0.05mm（仍满足加工公差），并启用“Small Edge”处理，仿真在第5次迭代后正常收敛，S11曲线与实测偏差小于0.8dB。
4.4 结果跳变：扫描参数微小变化引起响应突变
核心原理：结果跳变是指扫描变量（如贴片长度L）从L=10mm变化到L=10.1mm时，S11从-25dB突变至-3dB，而理论预期应为平滑变化。这种跳变通常源于模式切换——当天线工作频率接近高次模的截止频率时，主模与高次模之间发生能量转移。另一种常见原因是端口参考面位置不当，导致S参数相位跳变。
判断方法：绘制S11随扫描变量的变化曲线，若出现“台阶状”或“尖峰状”跳变，需进一步分析：
- 在跳变点附近加密扫描步长（如从0.1mm改为0.01mm），确认是否为连续过渡；
- 查看该频点处的电场分布，判断是否出现高次模（如TM11模）；
- 检查端口激励是否包含多个模式（Wave Port的Modes设置），若仅设置1个模式，可能遗漏模式耦合。
4.5 仿真中断：求解器异常终止
核心原理：仿真中断包括软件崩溃、内存溢出、许可证超时或网格生成失败。其中内存溢出最常见于电大尺寸天线阵列或包含精细结构的模型。中断的根因往往在网格剖分阶段——当模型包含极小特征（如过孔间距0.1mm）时，四面体网格数量可能超过1000万，导致内存不足。
应对策略：
- 启用“Domain Decomposition”方法，将模型分割为多个子域并行求解；
- 使用“Mesh Fusion”技术，对不同区域采用不同网格密度；
- 若模型包含周期性结构，使用“Unit Cell”边界条件简化；
- 检查日志文件中的错误代码，如“Out of Memory”或“Mesh Generation Failed”。
| 中断类型 | 典型日志信息 | 修复优先级 | 推荐措施 |
|---|---|---|---|
| 内存溢出 | Error: insufficient memory | 高 | 启用64位求解器，增加RAM |
| 网格失败 | Mesh generation failed at element 12345 | 高 | 检查几何重叠，简化倒角 |
| 许可证超时 | License checkout timed out | 中 | 缩短单次仿真时间，分批扫描 |
| 求解器崩溃 | Solver terminated unexpectedly | 低 | 更新软件版本，检查GPU驱动 |
4.6 综合排障流程与复盘
面对参数扫描收敛失败，建议遵循“先模型后求解器、先几何后数值”的排查顺序。首先检查模型几何是否存在重叠、缝隙、尖角等奇异特征；其次验证材料参数和边界条件是否物理合理；然后调整网格剖分策略；最后才修改求解器设置。下表总结了六种现象与修复措施的对应关系：
| 现象 | 首要检查项 | 次要检查项 | 推荐修复动作 |
|---|---|---|---|
| 不收敛 | 最小特征尺寸 | 迭代次数设置 | 简化薄层，增加迭代上限 |
| 伪收敛 | 场图验证 | 能量守恒 | 加密网格，启用自适应 |
| 振荡 | 高Q结构 | 频率点密度 | 启用Coupled模式 |
| 发散 | 几何重叠 | 材料参数 | 修复几何，检查εr/μr |
| 结果跳变 | 模式切换 | 端口模式数 | 细化扫描步长 |
| 仿真中断 | 内存资源 | 网格数量 | 域分解，简化模型 |