辐射发射分析
辐射发射理论基础
概要 — 什么是辐射发射
辐射发射怎样进行仿真？要把整个机柜网格划分吗？
很好的问题。直白地说，把整个机柜网格化一次求解在现实中不可行。PCB图案是mm以下的细节，但3m法测试距离是3,000mm——尺度相差3个数量级以上。想用单个网格来解这样的问题，会导致数十亿个单元，计算代价无法接受。
那怎样处理呢？
机柜的狭缝和缆线出口是辐射的关键通道。实务中常采用PCB级近场分析、随后将结果传递给机柜级FDTD或MoM进行计算的多阶段方法。3m法的测量关联精度为±6dB的工业标准。
±6dB不是很大的偏差吗？那是4倍的功率差。
眼光敏锐。从电压比看是2倍。在EMC领域，±6dB被认为是"实用可信赖的"精度。因为缆线配置和连接器接触阻抗偏差，量产产品本身就能产生5～10dB的变动。所以工程设计中是按6dB裕度来做的。
辐射发射（Radiated Emission, RE）是指电子设备意外放出的电磁波。CISPR 32（信息设备）或CISPR 25（车载设备）等国际标准规定了30 MHz～6 GHz（因标准而异）频段的限制。如果EMC测试不合格，就无法推向市场。因此在设计阶段进行仿真预测就显得极其重要。
辐射机制 — 差模与共模放射
辐射发射从主板的什么地方产生？是从整个板子均匀放出的吗？
有两条主要路径。差模（DM）放射和共模（CM）放射。
- 差模放射：信号电流和回流在PCB上形成环路，作为微小磁偶极子辐射。辐射强度与频率的4次方成正比，所以高频时钟信号是主要问题来源。
- 共模放射：整条缆线同相带电，作为单极子天线工作。仅数微安共模电流就能超过CISPR限值，EMC试验失败的8成是它引起。这是实际现场的认知。
只有数微安共模电流就会NG？这么微小的电流也能超限？
对。比如长度1m的缆线上流有5微安共模电流，在100 MHz会在3m距离产生约40 dBμV/m的电场强度。CISPR 32 Class B限值是30 dBμV/m，轻松超出。所以准确建模共模电流路径是仿真成败的关键。
支配方程式与远场变换
辐射发射分析的基础是麦克斯韦方程组。
从近场到远场的变换采用等效电磁流的远场变换（Far-Field Transformation）。在闭合曲面 $S$ 上已知电场 $\mathbf{E}_s$ 和磁场 $\mathbf{H}_s$ 时，距离 $r$ 处的远场电场 $\mathbf{E}_{far}$ 可用以下公式计算：
其中 $k = 2\pi/\lambda$ 是波数，$\eta_0 = 377\,\Omega$ 是自由空间阻抗，$\hat{n}$ 是曲面法向矢量。通过这种变换，只需在机柜附近的小计算域内求解，随后可高速计算任意测量距离的电场强度。
惠更斯等效原理（Huygens Equivalent Source）
惠更斯等效原理在多阶段分析的哪个环节用到？
这是"连接各阶段的胶水"。PCB级分析得到的近场被转换成等效电流源输入到机柜级分析中。具体来说，PCB周围的闭合曲面上的电磁场转换成"等效电流源"。
惠更斯面 $S_H$ 上的等效电流源定义为：
$\mathbf{J}_s$ 和 $\mathbf{M}_s$ 在曲面上每点都有定义，这些等效源能够完全复现闭合面外侧的电磁场。这就是等效定理（Love's Equivalence Theorem）。多阶段分析中，内侧计算得到的近场数据通过惠更斯面转换成等效源，传递给外侧分析。
缆线共模辐射模型
前面说共模电流是EMC试验失败的主因。具体的放射模型是什么样的？
Clayton Paul的经典模型很实用。长度 $L$ 的缆线上流有共模电流 $I_{CM}$，在测距离 $d$ 处的电场强度约为：
转换成dBμV/m：
其中 $f$ 是频率 [Hz]，$L$ 是缆线长度 [m]，$I_{CM}$ 是共模电流 [A]，$c$ 是光速，$d$ 是测距离 [m]。当缆线长度为半波长整数倍 $L = n\lambda/2$ 时出现共振，放射达到极大。
原来缆线长度和频率的关系影响放射。那缩短缆线就有利？
完全同意。缩短缆线是EMC对策的基本原则。但这样只是让共振频率偏移，宽带上看可能在别的频率出现新峰值。根本对策是抑制共模电流本身，铁氧体磁芯或缆线出口屏蔽处理是有效措施。
辐射发射数值计算方法
多阶段分析方法
多阶段方法具体怎样划分"阶段"？
典型分为3个阶段。
- 第1阶段 — PCB近场分析：用电路仿真器（SPICE等）或2.5D电磁场求解（Ansys SIwave、Cadence Sigrity等）算出PCB图案电流分布，获得基板上方的近场。
- 第2阶段 — 机柜级分析：把第1阶段的近场作为惠更斯等效源输入机柜模型。用FDTD或MoM计算机柜狭缝、缆线出口处的泄漏。
- 第3阶段 — 远场变换：从机柜外表的惠更斯面数据做远场变换，得到3m/10m处的电场强度。与CISPR限值比对。
阶段间数据交接会不会精度下降？
很好的提问。惠更斯面采样密度不足会丢失高频成分。经验法则是间距不超过λ/10。1 GHz时需30mm间距，3 GHz时需10mm。CST Studio Suite有"System Assembly and Modeling（SAM）"功能自动化阶段交接，实务上广泛应用。
FDTD法（时域有限差分法）
FDTD在Yee网格上将电场和磁场在时间和空间上错开半步进行计算。在辐射发射分析中，一次计算可获得宽频（30 MHz～数GHz）的谱，这是最大优势。
稳定条件（CFL条件）：
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 单元大小 | ≤ λmin/20 | 基于最高频率波长 |
| 吸收边界 | CPML（8～12层） | 比PML更宽带稳定 |
| 激励源 | 高斯脉冲/调制正弦波 | 根据分析带宽选择 |
| 截断条件 | 剩余能量−30 dB | 不足会导致低频误差 |
MoM（矩量法）
MoM和FDTD有什么区别？怎样选择？
MoM只离散导体表面，对于缆线和机柜这类开放空间问题很有利。不需要体积网格，计算量在大空间问题上大幅下降。但非均质材料（介质填充）就不擅长。Altair Feko是代表的MoM求解器，在缆线共模辐射预测上很强。
MoM基本方程（EFIE：电场积分方程）：
其中 $G(\mathbf{r},\mathbf{r}') = e^{-jk|\mathbf{r}-\mathbf{r}'|} / (4\pi|\mathbf{r}-\mathbf{r}'|)$ 是格林函数。MoM用RWG（Rao-Wilton-Glisson）基函展开上式，得到密矩阵方程。大规模问题用多极展开法（MLFMM）把计算量减到 $O(N \log N)$。
高频FEM
FEM对复杂形状和非均质材料有很好的适应性，Ansys HFSS是代表求解器。采用边元素（Nedelec元）离散矢量亥姆霍兹方程。
| 方法 | 擅长领域 | 弱点 | 代表工具 |
|---|---|---|---|
| FDTD | 宽带，瞬态分析 | 曲面离散，共振结构 | CST Studio Suite |
| MoM | 缆线辐射，开放问题 | 介质，体积问题 | Altair Feko |
| FEM | 复杂形状，多物理 | 宽带计算代价大 | Ansys HFSS |
网格要求与计算成本
辐射发射仿真的网格，和结构分析的标准不同吧？
完全不同。结构分析只需细分应力集中区，但电磁场分析的基准是最高频率的波长。要分析到6 GHz，波长50mm——需要单元大小2.5mm以下。机柜300mm×200mm×100mm，就是120×80×40 = 约40万个单元。再加PML区域，接近70万单元。
70万个单元……再加入PCB图案细节就是天文数字。
正是这个理由才需要多阶段方法。PCB图案细度0.1mm级，全部纳入会产生数十亿单元，不现实。PCB用2.5D求解器单独计算，然后以等效源形式传递——这是定式做法。
辐射发射实务应用
分析流程
实际做辐射发射仿真时，第一步是什么？
步骤是这样的。
- 确定噪声源：列出时钟频率及其谐波、电源纹波频率。和CISPR规范限值对照，划定有问题的频段。
- 获取PCB近场：通过电路仿真或2.5D电磁分析得PCB电流分布。有近场扫描仪测量数据最好，这样模型误差更小。
- 建立机柜+缆线模型：从CAD导入机柜形状。狭缝、通气孔、连接器开口要忠实建模。缆线用细线近似表示。
- 执行FDTD/MoM求解：把PCB近场作为惠更斯等效源放在机柜位置，计算整个机柜。
- 远场变换+CISPR对比：算出3m/10m处的电场，和规范限值作对比图。找出裕度不足的频段列为对策候选。
PCB近场分析
PCB近场是仿真还是实测好？
理想是两个都做，但设计初期仿真，试产后用实测更现实。近场扫描仪（EMSCAN RFR、Detectus等）可测基板上方10mm处的磁场分布，作为惠更斯等效源输入机柜分析，关联精度很好。纯仿真会漏掉器件寄生参数和PCB生产偏差。
近场扫描仪很贵吧？有便宜办法吗？
便宜的办法是自制近场探针（改装同轴线）配合频谱仪。定量精度差，但找出噪声源位置足够。还有用OpenEMS算PCB 2D电流分布，导入CST的工作流，论文里有报道。
机柜屏蔽效果评估
金属机柜的屏蔽效果(SE)在无开口时极高。但实际机柜有狭缝、通气孔、连接器开口，大幅降低屏蔽效果。
狭缝大小和屏蔽效果有什么规律吗？
经典规则是：狭缝长度达到波长的一半时共振，屏蔽效果接近零。比如300 MHz的波长1m，50cm狭缝就完全泄漏。实务中狭缝长度控制在最高频波长的1/20以下。若要6 GHz级，狭缝要不超过5mm。
还有一点重要的是多个狭缝的排列效果。通气孔多了，单个狭缝的放射会叠加增强。仿真必须准确反映狭缝数量和位置，否则实测对不上。
3m法/10m法测量关联
仿真结果和实测的"关联"怎样评估？±6dB是什么意思？
一般用FSV（Feature Selective Validation）定量评估，是IEEE标准（IEEE 1597.1）。实务中常把各频率峰值偏差在±6dB内，光谱包络形状一致视为"良好关联"。
关联差的主因通常是这些。
- 缆线配线路径不一致：实机缆线走线和模型不同，共模电流路径改变，结果偏差最大。
- 机柜接触阻抗：盖子和底座螺钉连接的接触电阻未建模。
- 地板面反射：3m法半无响室有金属地板反射。低频（30～200 MHz）如不建模会30dB级大偏差。
- 天线系数标定误差：测量天线的AF标定精度本身±2dB。
常见失败与对策
| 现象 | 原因 | 对策 |
|---|---|---|
| 低频(30～200 MHz)偏差超10dB | 地板反射未考虑 | 模型中加金属地板PEC |
| 特定频率峰值偏差超20dB | 缆线共振和实机不符 | 精确建模缆线长度和配置 |
| 1 GHz以上仿真明显偏低 | FDTD网格太粗/狭缝未建模 | 网格到λ/20以下，狭缝尺寸准确 |
| 全帯宽仿真偏高 | 铁氧体或屏蔽处理未建模 | 加入铁氧体阻抗特性 |
| FDTD计算不稳定 | 违反CFL条件/网格过细 | 用局部细化(subcell)技术 |
"试验场失败延期3个月"——设计仿真的经济价值
EMC试验通常在外部认证实验室进行，预约往往数周到数月后。量产前试验中辐射发射NG——改机柜设计→金模修改→重新试验，3个月延期不罕见。某大家电厂商的估算："一次EMC重试的成本=200万试验费+500万金模+3000万销售机会成本≈4000万元"。设计阶段仿真成本数百万，投资回报率超10倍。
辐射发射软件比较
主流工具对比
做辐射发射仿真用哪个软件好？
主要选择有4个。各有所长，要根据用途选。
| 工具 | 开发商 | 主要方法 | 强项 | EMC特有功能 |
|---|---|---|---|---|
| CST Studio Suite | Dassault Systemes SIMULIA | FDTD + FEM + MoM | 多求解器一体，SAM功能 | EMC模板，缆线自动生成 |
| Ansys HFSS | Ansys Inc. | FEM(自适应网格) | 高精度，SIwave联动 | 近远场变换，CISPR对比后处理 |
| Altair Feko | Altair Engineering | MoM + MLFMM | 缆线辐射，大型系统 | 缆线束求解，车载EMC |
| Keysight EMPro | Keysight Technologies | FEM + FDTD | ADS联动，电路协同 | PCB导入，电路-电磁耦合 |
开源的能做辐射发射分析吗？
OpenEMS（FDTD）开源可用。从MATLAB或Octave定义模型，执行分析。精度不逊商用，但无GUI，建模费时。学生研究或特定问题验证够用。还有gprMax（地探用，但FDTD通用）也有人用做EMC。
选型指南
最后怎样选？预算有限时？
选型的3个关键点。
- 模型复杂度：机柜狭缝多用FDTD（CST），缆线主导用MoM（Feko），复杂介质结构用FEM（HFSS）。
- 既有工作流：已用Ansys结构热分析就选HFSS联动。用3DEXPERIENCE就CST一家。Cadence PCB设计就Sigrity。
- 成本：商用软件年费数百万到千万。小企业或研究室可用OpenEMS+近场扫描仪实测数据的组合，可能更实惠。
CST的历程 — 德国小创业团队成为业界标准
CST（计算机仿真技术）1992年成立于德国达姆施塔特。其FDTD微波工作室在天线和EMC领域迅速占有市场，2016年被达索系统（CATIA/SOLIDWORKS母公司）收购。收购价据传超3亿欧元。现已并入SIMULIA产品族，能与Abaqus（结构）、PowerFLOW（流体）做多物理耦合。EMC仿真市场CST Studio Suite的份额最高。
辐射发射先端研究
机器学习EMC预测
最近流行的AI/机器学习也能用在辐射发射分析上吗？
能，而且研究很活跃。主要有两个思路。
- 代理模型：用FDTD结果作教师数据，训练神经网络从PCB布局参数直接预测放射谱。学完后数秒得预测，对设计搜索和裕度评估加速效果显著。
- PINN（物理信息神经网络）：把麦克斯韦方程纳入损失函数。教师数据少也能做出物理合理的预测。但目前3D宽带辐射问题的实用例还不多。
那以后就不用跑FDTD了？
部分是对。"过去设计的类似品"用代理模型数分钟就能判断Go/NoGo。但全新筐体形状或未知放射机制，还需物理级的全波分析。ML用做筛选，全波用做最终确认的分工应该是现实路线。
数字孪生与实时EMC
汽车和航空工业在构建"EMC数字孪生"，在设计阶段保证整体系统的EMC适合性。数十个车载ECU和数km电线束的辐射发射，用分层多尺度模型预测。CST Studio Suite的SAM和Ansys EMC Plus已支持这一用途。
随着5G/毫米波（24～40 GHz）普及，传统EMC规范未覆盖的频段需要新评估。毫米波下PCB通孔和芯片封装的共振变成新的放射源，封装级电磁分析的地位上升。
辐射发射故障排查
常见错误与解决方案
辐射发射分析中"容易掉坑"的失误有哪些？
EMC特有的陷阱有不少。结构和流体分析完全不同的地方在哪里掉坑，列一下。
1. FDTD计算发散（数值不稳定）
现象：时间步推进中电场值指数增长，最后异常终止。
原因：CFL条件不满足的时间步，或PML内数值不稳定。为了建模微小狭缝，局部网格极细，导致CFL条件要求的时间步全局过小，计算时间爆炸。
对策：用子网格技术（Thin Sheet Approximation等）保持网格均一，同时表现微小结构。CST的"Thin Panel"设定就是这个。
2. 共模电流被低估
现象：仿真合格，实机EMC试验不合格。特别是接缆线后恶化。
原因：缆线共模阻抗和PCB接地结构连接建模不准。连接器脚位分配或屏蔽接地处理未体现在模型。
对策：不简化连接器，精细建模屏蔽接地点。用MoM型缆线模型（Feko等）验证敏感度。
3. 近远场变换的位相误差
现象：远场辐射图形和实测明显不符（非水平偏差）。
原因：惠更斯面采样密度不足或面距离源太近，会丧失倏逝波影响。
对策：惠更斯面距源≥λ/4，面上采样间距≤λ/10。
4. 半无响室地板反射未考虑
现象：30～200 MHz实测偏差超过±10dB。
原因：3m法半无响室有金属地板。反射波与直达波干涉，自由空间计算完全不同。
对策：模型中加PEC完全导体地板。CST用"Electric Boundary"，HFSS用"Infinite Ground Plane"。
天哪，地板反射也得考虑……。自由空间算了结果能用做3m法预测吗？
不行。EMC试验的"3m法"前提就是半无响室，仿真也得复现同样试验环境才有意义。反过来说，自由空间的仿真结果当做"3m法预测"用是错误的。读透试验规范，忠实建模测量条件，是提高关联精度最确实的方法。
细节
错误