电力电子开关EMI分析——利用CAE预测和应对由dV/dt引起的传导与辐射噪声

问得好。关键在于“电压变化的速度”。比如最近的SiC MOSFET，开关动作在50ns以内就能完成。此时的dV/dt甚至能达到10kV/μs。

从时间轴上看，陡峭的电压变化“几乎是矩形波”，但对其进行傅里叶变换后，会发现开关频率整数倍的谐波一直延伸到MHz频段。例如，上升时间 $t_r = 20$ ns时，转折频率 $f_2 = 1/(\pi t_r) \approx 16$ MHz。在此频率以上也存在成分。

正是如此。PCB的布线图案和线束，当其长度相对于波长不可忽略时，就会作为高效的天线发挥作用。例如，100MHz时波长为3m。如果电缆长度达到 $\lambda/4 = 75$ cm，就会发生谐振。这就是辐射EMI。另一方面，通过电源线传播并影响其他设备的是传导EMI，150kHz〜30MHz频段是CISPR标准的对象。

从梯形波的傅里叶级数展开可知，频谱包络可以表示如下。设开关频率 $f_{sw}$、占空比 $D$、上升时间 $t_r$：

这个频谱有两个转折点：

• 第一转折点 $f_1 = 1/(\pi D T_{sw})$：从此处开始以 -20dB/dec 衰减
• 第二转折点 $f_2 = 1/(\pi t_r)$：从此处开始以 -40dB/dec 急剧衰减

例如，$f_{sw} = 100$ kHz、$D = 0.5$、$t_r = 20$ ns时，$f_1 \approx 64$ kHz、$f_2 \approx 16$ MHz。也就是说，直到16MHz都只以-20dB/dec衰减。如果SiC的$t_r$变为10ns，$f_2$会上升到32MHz，高频噪声会进一步增加。

这正是电力电子设计的困境。开关损耗 $P_{sw} \propto t_r$（越快越少）和EMI $\propto 1/t_r$（越快越多）是此消彼长的关系。在实际应用中，通过栅极电阻有意降低dV/dt的“压摆率控制”以及软开关电路（ZVS/ZCS）被广泛使用，原因就在于此。

用水管来比喻就容易理解了。DM（差模）噪声是去路和回路管道中“反向”流动的噪声电流。由开关电流的纹波成分引起，表现为电源线两线间的差动电压。

CM（共模）噪声是两线“同向”流动，通过接地平面（机壳）返回的噪声。其产生机制是，MOSFET的dV/dt通过其与散热片之间的寄生电容 $C_{parasitic}$ 产生共模电流 $I_{CM} = C_{parasitic} \cdot dV/dt$。

正是如此。例如，SiC MOSFET模块与散热片之间的寄生电容 $C_{parasitic} = 50$ pF，$dV/dt = 10$ kV/μs时，瞬间会产生 $I_{CM} = 50 \times 10^{-12} \times 10 \times 10^{9} = 0.5$ A的共模电流。如果这电流通过约30cm的电缆，就会在100MHz频段高效地辐射出去。

对策也不同。DM噪声使用X电容（线路滤波器）和DM扼流圈。CM噪声使用Y电容（线路与地之间）和CM扼流圈。这里如果搞错了，就会导致“加了滤波器但噪声没减少”这种典型失败。

正是如此。EMI分析的基础是麦克斯韦方程组。分析辐射EMI时需要全部4个方程：

不过，对于传导EMI（150kHz〜30MHz）的分析，基于电路理论的方法往往更现实。通过电路仿真（基于SPICE）计算经由LISN阻抗 $Z_{LISN} = 50\,\Omega$ 的噪声端子电压，再转换到频域。

没错。大致总结如下：

• 传导EMI（150kHz〜30MHz）：电路仿真（SPICE + 寄生元件模型）是主力。将LISN连接到等效电路上计算端子电压。
• 辐射EMI（30MHz〜1GHz+）：3D电磁场仿真（FDTD/FEM/MoM）必不可少。需要求解包含PCB布局和机壳形状的完整模型。
• 两者的交界频段（10〜100MHz）：电路-电磁场联合仿真（co-simulation）有效。

视情况而定，但大致可以这么说：

• 时域分析：直接输入开关波形，通过FFT得到频谱。可以直接使用非线性器件（MOSFET、二极管）的模型。能一次性得到宽频带结果。但时间步长需要小于最高频率的1/20，计算成本高。
• 频域分析：高效计算特定频率下的响应。适用于滤波器的插入损耗或S参数的计算。但无法处理开关的非线性，前提是线性化模型。

多数情况下是“用时域电路仿真（基于SPICE）计算开关波形，再用FFT求频谱”这个流程。对于辐射EMI，则将得到的频谱作为电磁场求解器的激励源输入，计算远场。这种两阶段方法最为实用。

我们来比较一下EMI分析中使用的三种主要方法：

像车载EMC这种“PCB + 线缆 + 机壳”的完整模型，FDTD是主力。CST Studio Suite的Time Domain Solver是其代表。另一方面，对于连接器或滤波器单体的S参数提取，基于FEM的Ansys HFSS很强。线束的辐射用基于MoM的方法效率高。最近混合方法（FEM-MoM耦合或FDTD-MoM耦合）也在增多。

电路-电磁场联合仿真（co-simulation）是EMI分析中最实用且强大的方法。有两种方式：

• 弱耦合（Loose Coupling）：用电路仿真计算开关波形 → 将其结果作为端口的激励源输入到电磁场求解器。计算是单向的，不考虑相互作用。
• 强耦合（Tight Coupling）：电路求解器和电磁场求解器在每个时间步长交换数据。PCB图案的寄生分量对开关波形的影响也能实时反映。Ansys Electronics Desktop的“Transient EM-Circuit co-simulation”是代表例子。

强耦合的计算成本要高好几个数量级。例如，对于100kHz开关频率，求解100个周期，如果电磁场时间步长为0.5ns，则需要2,000,000步。每一步都要解3D电磁场，所以一个案例可能需要数天。弱耦合的话，电路侧只需秒级，电磁场侧可以单独优化，因此对于设计初期的参数化研究，弱耦合更实用。

EMI分析的网格有其独特的规则。最重要的是考虑趋肤深度（Skin Depth）：

例如，铜（$\sigma = 5.8 \times 10^7$ S/m）在频率100MHz时，$\delta \approx 6.6\,\mu$m。电流集中在这个厚度的导体表面，因此网格在趋肤深度内至少需要3〜4层。EMI分析需遵守以下原则：

• 导体表面：单元厚度小于趋肤深度 $\delta$ 的1/3
• 空间：小于最高分析频率波长 $\lambda$ 的1/10（FDTD推荐$\lambda/20$）
• PCB过孔周边：局部网格小于过孔直径的1/4
• 薄电介质层：FR-4厚度方向至少2〜3层

Ansoft Maxwell。低周波電磁場解析。2008年Ansysに統合。

現在の所属: Ansys Inc.

Ansoft Corporationが開発した3D高周波電磁界シミュレータ。2008年にAnsysがAnsoftを買収。

現在の所属: Ansys Inc.

1986年スウェーデンで設立。MATLAB連携のFEMLABとして開始、後にCOMSOLに改名。マルチフィジックスに強み。

現在の所属: COMSOL AB

異なるソルバー間でモデルを変換する際は、要素タイプの対応関係、材料モデルの互換性、荷重・境界条件の表現差異に注意が必要になるんだ。特に高次要素や特殊要素（コヒーシブ要素、ユーザー定義要素等）はソルバー間で直接変換できない場合が多い。

メッシュ収束性の確認、境界条件の妥当性検証、材料パラメータの感度分析がすごく大事なんだ。

• メッシュ依存性の検証: 少なくとも3水準のメッシュ密度で収束性を確認
• 境界条件の妥当性: 物理的に意味のある拘束条件の設定
• 結果の検証: 理論解、実験データ、既知ベンチマーク問題との比較

うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。

形状関数 $N_i$ を用いて未知量を近似:

これを数式で表すとこうなるよ。

これを数式で表すとこうなるよ。

連続体の支配方程式を離散化すると、以下の代数方程式系が得られる:

ここで $[K]$ は全体剛性マトリクス（または同等のシステムマトリクス）、$\{u\}$ は未知節点変数ベクトル、$\{F\}$ は外力ベクトルなんだ。

• 完全積分: 全ての項を正確に積分。剛性過大評価の傾向（ロッキング）
• 低減積分: 積分点数を削減。計算効率向上だが、アワーグラスモード発生のリスク
• 選択的低減積分 (B-bar法): 体積項と偏差項を分離して積分。ロッキング回避

• h-refinement: メッシュを細分化（要素サイズ h を小さく）して精度向上
• p-refinement: 要素の多項式次数を上げて精度向上
• hp-refinement: h と p を同時に最適化

収束速度: 二次要素で $O(h^2)$ のオーダーで誤差が減少（滑らかな解の場合）

スイッチングノイズEMI解析の実務的な解析フローと注意点を解説する。

1. 前処理 (Pre-processing)

• CADデータのインポートと形状簡略化
• 材料特性の定義
• メッシュ生成（要素タイプ・サイズの決定）
• 境界条件と荷重条件の設定

2. 求解 (Solving)

• ソルバー設定（解法、収束基準、出力制御）
• ジョブ投入と計算実行
• 収束モニタリング

3. 後処理 (Post-processing)

• 結果の可視化（変位、応力、その他の物理量）
• 結果の検証と妥当性確認
• レポート作成

• 応力集中部: 最低3層以上の要素を配置
• 応力勾配の大きい領域: 要素サイズを周囲の1/3〜1/5に
• 荷重印加点近傍: 局所細分化
• 遠方領域: 粗いメッシュで計算効率を確保

• 過拘束に注意: 剛体移動の拘束は6自由度のみ
• 対称条件の活用: 計算規模の削減
• 荷重の等価分配: 集中荷重 vs. 分布荷重の選択

Ansoft Maxwell。低周波電磁場解析。2008年Ansysに統合。

現在の所属: Ansys Inc.

Ansoft Corporationが開発した3D高周波電磁界シミュレータ。2008年にAnsysがAnsoftを買収。

現在の所属: Ansys Inc.

• メッシュ収束性を3水準以上で確認したか
• 力の釣り合い（反力合計）を検証したか
• 結果が物理的に妥当な範囲か確認したか
• 既知の理論解またはベンチマーク問題と比較したか

うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。

スイッチングノイズEMI解析に対応する主要な商用CAEツールの機能比較と、各製品の歴史的背景を詳述する。

Ansoft Maxwell。低周波電磁場解析。2008年Ansysに統合。

現在の所属: Ansys Inc.

Ansoft Corporationが開発した3D高周波電磁界シミュレータ。2008年にAnsysがAnsoftを買収。

現在の所属: Ansys Inc.

1986年スウェーデンで設立。MATLAB連携のFEMLABとして開始、後にCOMSOLに改名。マルチフィジックスに強み。

現在の所属: COMSOL AB

Computer Simulation Technology (ドイツ) が開発。2016年にDassault Systèmesが買収しSIMULIAに統合。

現在の所属: Dassault Systèmes SIMULIA

• 要素タイプの非互換: ソルバー固有要素は中立フォーマットで表現不可
• 材料モデルの差異: 同名でも内部実装が異なる場合がある
• 境界条件の再定義: 多くの場合、手動での再設定が必要
• 結果データの比較: 出力変数の定義（節点値 vs. 要素値、積分点値）に差異

スイッチングノイズEMI解析のツール選定においては以下を考慮:

• 解析規模: 数万〜数億DOFへのスケーラビリティ
• 物理モデル: 必要な構成則・要素タイプの対応状況
• ワークフロー: CADとの連携、自動化の容易さ
• コスト: 初期投資 + 年間保守 + 教育コスト
• サポート: 技術サポートの質とレスポンス

うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。

スイッチングノイズEMI解析における最新の研究動向と先進的手法を見ていこう。

• 等幾何解析 (IGA): NURBS基底関数を直接使用し、CAD-CAE間のシームレスな連携を実現
• 粒子法 (SPH, MPM): メッシュフリー手法による大変形・破壊の追跡
• 位相場法 (Phase-Field): 界面の暗示的表現による複雑な界面追跡
• 機械学習支援: サロゲートモデル、物理インフォームドニューラルネットワーク (PINN)

症状: ソルバーが指定反復回数内に収束せず異常終了

考えられる原因:

• メッシュ品質の不足（過度に歪んだ要素）
• 材料パラメータの不適切な設定
• 不適切な初期条件
• 非線形性が強すぎる（荷重ステップの不足）

対策:

• メッシュ品質チェックを実施（アスペクト比、ヤコビアン）
• 材料パラメータの単位系を確認
• 荷重を複数ステップに分割（サブステップ数の増加）
• 収束判定基準の緩和（ただし精度に注意）

症状: 応力/変位/温度等が物理的に非現実的な値

考えられる原因:

• 境界条件の誤設定
• 単位系の混在（SI単位と工学単位の混同）
• 不適切な要素タイプの選択
• 応力特異点の存在

対策:

• 反力の合計を確認（力の釣り合い）
• 単位系の一貫性を確認
• 要素タイプの適切性を再検討
• 特異点除去またはサブモデリング

症状: 計算が想定時間の何倍もかかる

対策:

• メッシュの粗密分布の最適化
• 対称性の活用（1/2, 1/4モデル）
• ソルバー設定の最適化（反復法、前処理の選択）
• 並列計算の活用

症状: Out of Memory エラー

対策:

• アウトオブコア解法の使用
• メッシュ規模の削減
• 64bit版ソルバーの使用確認
• メモリ割り当ての増加

• FATAL 2012: 特異剛性マトリクス → 拘束条件の見直し
• USER WARNING 5291: 要素品質不良 → メッシュ修正
• SYSTEM FATAL 3008: メモリ不足 → MEM設定の調整

• Excessive distortion: 要素の過大変形 → NLGEOM確認、メッシュ改善
• Zero pivot: 拘束不足 → 境界条件追加
• Time increment too small: 収束失敗 → ステップ設定見直し