パワエレ・スイッチングEMI解析 — dV/dt起因の伝導・放射ノイズをCAEで予測・対策する
理論と物理
なぜスイッチングでEMIが発生するのか
先生、インバータのスイッチングでEMIが出るのはなぜですか? 直流を交流に変えるだけなのに、そんなにノイズが出るものなんですか?
いい質問だね。ポイントは「電圧の変化の速さ」だ。たとえば最近のSiC MOSFETは、50ns以下でスイッチングが完了する。このときのdV/dtは10kV/μsにも達するんだ。
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規格の対象になる。
dV/dtとスペクトルエンベロープ
スイッチング波形のスペクトルって、定量的にはどうなるんですか? 式で見たいです。
台形波のフーリエ級数展開から、スペクトルのエンベロープは次のように表せる。スイッチング周波数 $f_{sw}$、デューティ比 $D$、立ち上がり時間 $t_r$ とすると:
このスペクトルには2つの折れ点がある:
- 第1折れ点 $f_1 = 1/(\pi D T_{sw})$:ここから -20dB/dec で減衰開始
- 第2折れ点 $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に上がり、高周波ノイズがさらに増える。
え、じゃあスイッチングが速いほどEMIは悪化するってことですか? スイッチング損失は減るのに…
まさにそれがパワエレ設計のジレンマだ。スイッチング損失 $P_{sw} \propto t_r$(速いほど少ない)と、EMI $\propto 1/t_r$(速いほど多い)はトレードオフの関係にある。ゲート抵抗でdV/dtを意図的に遅くする「スルーレート制御」や、ソフトスイッチング回路(ZVS/ZCS)が実務で多用される理由がこれだ。
コモンモード(CM)/ディファレンシャルモード(DM)ノイズ
EMI対策でよく「コモンモード」「ディファレンシャルモード」って聞くんですけど、違いがいまいちピンときません。
水道管にたとえると分かりやすい。DM(ディファレンシャルモード)ノイズは、行きの管と帰りの管を「逆方向に」流れるノイズ電流だ。スイッチング電流のリプル成分が原因で、電源ラインの2線間に差動電圧として現れる。
CM(コモンモード)ノイズは、2線が「同じ方向に」流れて、グランド面(シャーシ)を経由して帰ってくる。発生メカニズムは、MOSFETのdV/dtが、ヒートシンクとの寄生容量 $C_{parasitic}$ を通じてコモンモード電流 $I_{CM} = C_{parasitic} \cdot dV/dt$ を流すことだ。
じゃあ寄生容量が大きいほどCMノイズも大きいんですか?
そのとおり。例えば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チョークコイルを使う。ここを間違えると、「フィルタを追加したのにノイズが減らない」という典型的な失敗になる。
| 項目 | DMノイズ | CMノイズ |
|---|---|---|
| 電流の向き | 2線で逆方向 | 2線で同方向(グランド帰還) |
| 主な発生源 | スイッチング電流リプル | dV/dt × 寄生容量 |
| 支配的帯域 | $f_{sw}$ 〜数MHz | 数MHz〜100MHz以上 |
| フィルタ | Xコンデンサ + DMチョーク | Yコンデンサ + CMチョーク |
| LISNでの検出 | $V_{DM} = (V_1 - V_2)/2$ | $V_{CM} = (V_1 + V_2)/2$ |
支配方程式とマクスウェル方程式
EMI解析って結局、マクスウェル方程式を解いているんですか?
そのとおり。EMI解析の根底にあるのはマクスウェル方程式だ。放射EMIの解析では4つの式すべてが必要になる:
ただし、伝導EMI(150kHz〜30MHz)の解析では、回路理論ベースのアプローチが現実的なことが多い。LISNのインピーダンス $Z_{LISN} = 50\,\Omega$ を介したノイズ端子電圧を回路シミュレーション(SPICEベース)で計算し、周波数領域に変換する。
じゃあ、伝導EMIと放射EMIで解析手法が違うんですね。
その通り。ざっくりまとめると:
- 伝導EMI(150kHz〜30MHz):回路シミュレーション(SPICE + 寄生素子モデル)が主力。等価回路にLISNを接続して端子電圧を計算
- 放射EMI(30MHz〜1GHz+):3D電磁界シミュレーション(FDTD/FEM/MoM)が必須。PCBレイアウトと筐体形状を含むフルモデルを解く
- 両者の境界域(10〜100MHz):回路-電磁界連成(co-simulation)が有効
EMIスペクトルのエンベロープ導出
周期 $T$ の台形波パルスのフーリエ級数は、$n$ 次高調波の振幅が:
$$ |c_n| = \frac{2A\tau}{T} \cdot \frac{\sin(n\pi f_{sw}\tau)}{n\pi f_{sw}\tau} \cdot \frac{\sin(n\pi f_{sw}t_r)}{n\pi f_{sw}t_r} $$ここで $A$ は振幅、$\tau$ はパルス幅、$t_r$ は立ち上がり時間。dB表記にすると3つの領域が現れる:
- $f < f_1 = 1/(\pi\tau)$ の領域:振幅一定(0dB/dec)。低周波では全高調波がフル振幅。
- $f_1 < f < f_2 = 1/(\pi t_r)$ の領域:-20dB/dec で減衰。sinc関数の第1引数が効く。
- $f > f_2$ の領域:-40dB/dec で急減衰。両sinc関数が効く。$t_r$ が短いほどこの折れ点が高周波に移動し、EMIが悪化。
CMノイズの等価回路モデル
コモンモードノイズの等価回路は以下のように構成される:
- ノイズ源:MOSFETのドレイン-ソース電圧 $V_{DS}(t)$ の dV/dt
- 結合パス:ヒートシンクとの寄生容量 $C_p$(典型値:10〜100pF)
- 帰還パス:グランドインピーダンス $Z_{GND}$(配線のインダクタンス $L_{wire}$ を含む)
- 終端:LISN $Z_{LISN} = 50\,\Omega$
CM電圧は $V_{CM} = V_{DS} \cdot \frac{Z_{LISN}\|Z_{GND}}{Z_{LISN}\|Z_{GND} + 1/(j\omega C_p)}$ で概算できる。
EVのインバータ開発で「CISPR 25に落ちた!」——現場で何が起きるか
電気自動車(EV)のモータ駆動インバータは、CISPR 25クラス5をパスしなければ量産に進めない。開発の終盤でEMC試験に落ちると、PCBの再設計やフィルタの追加が必要になり、数か月のスケジュール遅延と数千万円のコストが発生する。ある自動車Tier1サプライヤーでは、SiC MOSFETの採用にあたりdV/dtを8kV/μsから5kV/μsに落としてEMCをパスしたが、スイッチング損失が40%増加し、放熱設計をやり直す羽目になった。「EMIは設計の最初から考えるもの」という教訓は、こうした痛い経験から生まれている。
数値解法と実装
時間領域 vs 周波数領域解析
EMI解析って、時間領域と周波数領域のどちらでやるのが正解ですか?
ケースバイケースだけど、ざっくり言うとこうだ:
- 時間領域解析:スイッチング波形をそのまま入力してFFTでスペクトルを得る。非線形デバイス(MOSFET、ダイオード)のモデルをそのまま使える。広帯域の結果が一度に得られる。ただし時間刻みは最高周波数の1/20以下が必要で、計算コストが高い。
- 周波数領域解析:特定の周波数での応答を効率的に計算。フィルタの挿入損失やS-パラメータの算出に適する。ただしスイッチングの非線形性は扱えず、線形化したモデルが前提。
実務ではどっちを使うことが多いんですか?
多くの場合は「時間領域の回路シミュレーション(SPICEベース)でスイッチング波形を計算し、FFTでスペクトルを求める」というフローだ。放射EMIについては、得られたスペクトルを電磁界ソルバーの励振源として入力し、遠方界を計算する。この2段階アプローチが最も実用的だね。
FDTD・FEM・MoMの使い分け
3D電磁界ソルバーにもいろいろありますよね? FDTD、FEM、MoM…どう使い分けるんですか?
EMI解析で使われる主要な3手法を比較しよう:
| 手法 | 原理 | EMI解析での得意分野 | 弱点 |
|---|---|---|---|
| FDTD | 時空間を格子で離散化、Yeeアルゴリズム | 広帯域の過渡解析。PCB+筐体の放射EMI | 曲面の近似が粗い。PML境界設定が重要 |
| FEM | 変分法で場を近似、辺要素 | 複雑形状。周波数領域のSパラメータ | 広帯域は周波数掃引が必要。計算コスト大 |
| MoM | 積分方程式を離散化 | ケーブルハーネス、アンテナ。開領域 | 体積問題(誘電体内部)は苦手。密行列 |
で、EMI解析ではどれが一番使われてるんですか?
車載EMCのような「PCB + ケーブル + 筐体」のフルモデルではFDTDが主力だ。CST Studio SuiteのTime Domain Solverがその代表。一方、コネクタやフィルタ単体のSパラメータ抽出にはFEMベースのAnsys HFSSが強い。ケーブルハーネスの放射はMoMベースが効率的。最近はハイブリッド手法(FEM-MoM連成やFDTD-MoM連成)も増えている。
回路-電磁界連成シミュレーション
回路シミュレーションと電磁界シミュレーションを一緒にやる「連成解析」ってどういう仕組みですか?
回路-電磁界連成(co-simulation)は、EMI解析で最も実用的かつ強力なアプローチだ。2つの方法がある:
- 弱連成(ルーズカップリング):回路シミュレーションでスイッチング波形を計算 → その結果をポートの励振源として電磁界ソルバーに入力。計算は1方向で、相互作用は考慮しない。
- 強連成(タイトカップリング):回路ソルバーと電磁界ソルバーが時間ステップごとにデータを交換。PCBパターンの寄生成分がスイッチング波形に与える影響もリアルタイムに反映される。Ansys Electronics Desktopの「Transient EM-Circuit co-simulation」が代表例。
強連成のほうが正確なんですよね? なぜ弱連成を使うケースがあるんですか?
強連成は計算コストが桁違いに大きい。例えば100kHzスイッチングで100周期分を解くとして、電磁界の時間刻みが0.5nsなら、2,000,000ステップ必要になる。各ステップで3D電磁界を解くから、1ケースに数日かかることもある。弱連成なら回路側は秒単位、電磁界側は別途最適化できるから、設計初期のパラメトリックスタディには弱連成が実用的だ。
メッシュ戦略と表皮深さ
EMI解析のメッシュで気をつけることってありますか? 構造解析とは違いますよね?
EMI解析のメッシュには独自のルールがある。最も重要なのは表皮深さ(スキンデプス)の考慮だ:
例えば銅($\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層
構造解析だと「応力集中部を細かく」ですけど、EMIだと「導体表面と波長」がポイントなんですね。
FDTDの時間刻みとCFL条件
FDTD法では安定条件(Courant-Friedrichs-Lewy条件)により時間刻み $\Delta t$ が空間メッシュサイズ $\Delta x, \Delta y, \Delta z$ で制約される:
$$ \Delta t \leq \frac{1}{c\sqrt{1/\Delta x^2 + 1/\Delta y^2 + 1/\Delta z^2}} $$メッシュを局所的に細かくすると全体の時間刻みが引きずられて計算コストが激増する。EMI解析では「非均一メッシュ + サブグリッディング」でこれを回避するテクニックが重要だ。
実践ガイド
EMI解析フロー
先生、実際にEMI解析をやるときの手順を教えてください。何から始めればいいですか?
EMI解析は以下の5ステップで進める:
- ノイズ源の特定:スイッチングデバイスのdV/dt、di/dtをSPICEまたはデータシートから取得。デバイスメーカーのSPICEモデル(Level 3以上)を使う
- 伝搬経路のモデリング:PCBレイアウト(Gerberデータから3Dモデル生成)、ケーブルハーネス、筐体のCADモデルを準備
- 寄生素子の抽出:PCBトレースのインダクタンス $L_{trace}$、パッドの容量 $C_{pad}$、ビアのインピーダンスを電磁界解析で抽出
- 回路-電磁界連成解析の実行:スイッチング回路にLISN($50\,\Omega / 50\,\mu$H)を接続し、伝導EMI端子電圧を計算。放射EMIは3Dモデルで遠方界を計算
- 規格リミットとの比較:CISPR 25(車載)、CISPR 32(情報技術機器)等のリミットラインと重ね合わせてマージンを確認
LISNって何ですか? よく出てきますけど。
LISN(Line Impedance Stabilization Network)は、電源ラインのインピーダンスを規格値(通常$50\,\Omega / 50\,\mu$H)に標準化する計測用ネットワークだ。シミュレーションでもLISNを回路モデルに入れないと実測と合わない。CISPR 16-1-2で規定されている回路をそのまま使えばいい。
EMIフィルタ設計と挿入損失
EMIフィルタの設計って、どこから始めればいいですか?
まずフィルタの必要挿入損失(IL: Insertion Loss)を求める。これはフィルタなしのEMIスペクトルとリミットラインの差分にマージン(通常6〜10dB)を加えたものだ:
次にフィルタのトポロジーを選ぶ。最も基本的なのはL-C型で、カットオフ周波数 $f_c$ は:
$f_c$ 以上の帯域で、2次LCフィルタなら-40dB/decの減衰特性を示す。例えば必要ILが60dBで-40dB/decの減衰なら、$f_c$ はノイズのピーク周波数より1.5桁(約30倍)低く設定する必要がある。
理論上はそうですけど、実際のフィルタって理想通りに効かないですよね?
鋭いね。実際のフィルタ部品には寄生成分がある。コンデンサのESL(等価直列インダクタンス)は数nHあり、自己共振周波数(SRF)以上ではコンデンサがインダクタとして振る舞う。10μFのセラミックコンデンサでも、SRFは1〜5MHz程度。つまり5MHz以上ではフィルタが効かなくなる。対策として:
- 小容量のコンデンサ(100pF〜1nF)を高周波用に並列追加
- フェライトビーズで100MHz帯をカバー
- EMC用のフィードスルーコンデンサ(ESLが極めて小さい)を使用
これらの寄生素子をモデルに含めないシミュレーションは実測と大きく乖離する。
| フィルタ種別 | 対象ノイズ | 典型的減衰特性 | 実装上の注意点 |
|---|---|---|---|
| π型LCフィルタ(DM) | DM 150kHz〜30MHz | -60dB/dec | 入出力のコンデンサは低ESL品を選定 |
| CMチョーク + Y-cap | CM 1MHz〜100MHz | -40〜-60dB/dec | チョークの飽和電流とインピーダンス特性を確認 |
| フェライトビーズ | CM/DM 30MHz〜1GHz | 部品依存 | DCバイアス特性(電流でインピーダンス低下)に注意 |
| シールドケーブル | 放射 30MHz〜 | シールド効果 40〜80dB | グランド接続方式(ピグテール vs 360度)が重要 |
PCBレイアウト最適化
フィルタだけじゃなくて、PCBのレイアウト自体もEMIに影響しますよね?
むしろレイアウトのほうが支配的だ。EMIの大原則は「ループ面積を最小化する」こと。ファラデーの法則から、ループが囲む面積 $A$ に比例してループに結合する磁束が増え、放射も増える:
具体的なレイアウト指針は:
- スイッチングループ:MOSFETのドレイン → DCバスコンデンサ → ソースまでのループを最小に。理想は10mm²以下
- グランドプレーン:内層にソリッドグランドを設け、リターン電流パスのインピーダンスを最小化
- ゲートドライブ:ゲート-ソース間のループもEMI源。ツイストペアまたは同軸で引き回す
- 部品配置:EMIフィルタはコネクタ直近に配置。フィルタの入力側と出力側でカップリングしないように物理的に分離
レイアウトの良し悪しをシミュレーションで確認できるんですか?
できる。Ansys SIwaveやCadence Sigrity PowerSIはGerberデータをそのまま読み込んで、PCBの寄生LCR抽出と電流分布の可視化ができる。「どのループにどれだけの高周波電流が流れているか」をカラーマップで確認し、問題箇所を特定する。これがEMI解析の最も実用的な使い方だ。
CISPR 25 / CISPR 32適合
CISPR 25のクラス5ってどのくらい厳しいんですか? 具体的な数値を知りたいです。
CISPR 25は車載電子機器のEMI規格で、クラス1(最も緩い)〜クラス5(最も厳しい)がある。クラス5の伝導エミッション限度値の一例:
| 周波数帯域 | クラス5 限度値(ピーク) | クラス5 限度値(準尖頭値) |
|---|---|---|
| 150kHz 〜 300kHz | 40 dBμV | 30 dBμV |
| 300kHz 〜 30MHz | 34 dBμV | 24 dBμV |
| 30MHz 〜 54MHz(放射) | 14 dBμV/m @ 1m | — |
| 76MHz 〜 108MHz(FM帯) | 6 dBμV/m @ 1m | — |
FM放送帯(76〜108MHz)でわずか6dBμV/m @ 1mという限度値は、事実上「ほとんど何も放射してはいけない」レベルだ。SiCインバータでこれをクリアするには:
- スイッチングループ面積の最小化(前述)
- 多段EMIフィルタ(LCL-π型 + CMチョーク)
- シールド筐体+フィルタコネクタ
- ゲート抵抗によるdV/dt制限(5〜8kV/μs程度に)
これらすべてを電磁界シミュレーションでバーチャルに検証してから試作に進む、というのが現代のEMC開発フローだ。
EMIフィルタの「グランド」問題——接地が悪いと逆効果
EMIフィルタを付けたのにノイズが減らない、むしろ増えた——という経験をした人は少なくない。原因のひとつがグランドラインのインピーダンスだ。Yコンデンサはシャーシグランドに接続するが、グランドラインが長かったり薄かったりすると、コンデンサのESLと共振してノイズを増幅してしまう。実務では「フィルタはシャーシの直近に」「グランドは幅広パターンで最短距離」が鉄則。EMI解析で接地ラインのインダクタンスを含めたフルモデルを組まないと、この共振は予測できない。
「6dBマージン」の重要性
EMC試験は再現性が完璧ではない。同じDUTでも試験環境(ケーブルの引き回し、グランドプレーンとの距離、温度)で数dBの変動が出る。CISPR 25の限度値ぴったりを狙うのは危険で、最低6dB(約2倍)のマージンを確保するのが業界の常識だ。マージンが3dB以下だと、量産品の個体差でNG品が出る確率が跳ね上がる。
ソフトウェア比較
EMI解析ツール比較
EMI解析に使える商用ツールって、結局どれがいいんですか?
EMI解析では「回路シミュレータ」「3D電磁界ソルバー」「PCB寄生抽出ツール」の3種類を組み合わせるのが普通だ。主要な選択肢を比較しよう:
| ツール | 開発元 | 主要手法 | EMI解析での強み |
|---|---|---|---|
| CST Studio Suite | Dassault Systemes | FDTD / FEM / MoM | 時間領域の広帯域EMI解析。車載EMCのデファクトスタンダード |
| Ansys HFSS | Ansys Inc. | FEM(周波数領域) | Sパラメータ、フィルタ設計。Electronics Desktop連携 |
| Ansys SIwave | Ansys Inc. | ハイブリッドFEM-MoM | PCB全層の寄生抽出。パワーインテグリティ |
| Keysight ADS | Keysight | 回路 + EM連成 | RF/マイクロ波帯のフィルタ設計 |
| Cadence Sigrity | Cadence | FEM / FDTD | PCB寄生抽出とPI/SI解析。OrCAD連携 |
| COMSOL Multiphysics | COMSOL AB | FEM | マルチフィジックス。熱-電磁連成 |
| Altium EMC Advisor | Altium | ルールベース | PCB設計段階でのリアルタイムEMCチェック |
無料で使えるツールはありますか? 学生なので予算がなくて…
いくつかあるよ:
- openEMS:オープンソースのFDTDソルバー。MATLABまたはOctaveから操作。基本的なEMC解析は可能
- FEMM:2D FEMベースの電磁場ソルバー。低周波EMIの基礎検討に
- LTspice:Analog Devices提供の無料SPICE。伝導EMIの回路シミュレーションに最適
- KiCad + openEMS連携:PCBレイアウトからFDTDモデルを生成するプラグインがある
学術用途ならAnsys、CST、COMSOLとも学生ライセンスを提供しているから、大学経由で申請するといい。
ツール選定の指針
結局のところ、どういう基準で選べばいいんですか?
EMI解析ツールの選定で重要なのは3つ:
- 対象周波数帯域:伝導EMI(150kHz〜30MHz)だけならSPICE+寄生抽出で十分。放射EMI(30MHz〜GHz)まで必要なら3Dソルバーが必須
- 既存設計フローとの連携:使用中のECAD(Altium、Cadence、Mentor)とのデータ連携がスムーズかどうか。Gerberデータの読み込み、ネットリスト同期、結果のバックアノテーション
- 回路-電磁界連成の成熟度:SPICEモデルとの連成が簡単にできるか。Ansys Electronics Desktopは回路-HFSS-SIwave-Icepakの統合環境が強み。CST Studioはデザインスタジオで回路-3D FDTDの強連成ができる
先端技術
WBGデバイスとEMI課題
SiCやGaNの次世代パワーデバイスが普及すると、EMIの問題はもっと深刻になるんですか?
間違いなくそうだ。WBG(ワイドバンドギャップ)デバイスの特徴を見ると:
| 特性 | Si IGBT | SiC MOSFET | GaN HEMT |
|---|---|---|---|
| スイッチング速度 | 200〜500ns | 20〜100ns | 5〜20ns |
| 典型的dV/dt | 1〜3kV/μs | 5〜15kV/μs | 50〜100kV/μs |
| スペクトル第2折れ点 $f_2$ | 0.6〜1.6MHz | 3〜16MHz | 16〜64MHz |
| EMI影響帯域 | 〜数MHz | 〜数十MHz | 〜数百MHz |
GaN HEMTではdV/dtが100kV/μsに達することもあり、$f_2$が64MHzを超える。これはFM放送帯(76〜108MHz)にEMIが直撃することを意味する。従来のSi IGBTでは問題にならなかった帯域が、WBGデバイスでは最大の課題になる。
じゃあGaNを使うなら、EMI対策をまったく新しいレベルで考えないといけないんですね。
そうだ。具体的には:
- アクティブゲートドライブ:dV/dtをリアルタイムで制御するフィードバック回路。EMIと損失のトレードオフを動的に最適化
- パッケージ内蔵フィルタ:GaNモジュール内にフィルタを内蔵し、外部に高周波ノイズを出さない設計
- PCB一体型放熱:ヒートシンクの寄生容量を極限まで小さくする基板埋め込みデバイス技術
- EMI予測付き設計自動化:レイアウトの段階でEMIを自動予測し、問題箇所をリアルタイムでフィードバック
機械学習によるEMI予測
最近はAIとか機械学習でEMI予測ってできるんですか?
学術研究で急速に進んでいる分野だ。代表的なアプローチは:
- サロゲートモデル:数百〜数千ケースのシミュレーション結果でニューラルネットワークを訓練し、パラメータ変更時のEMIスペクトルをミリ秒で予測。設計空間の探索が桁違いに速くなる
- CNN(畳み込みニューラルネット)によるPCBレイアウト評価:PCBの画像データを入力し、EMI放射の強さを予測。Gerberデータを直接入力する研究もある
- 物理インフォームドニューラルネットワーク(PINN):マクスウェル方程式を損失関数に組み込み、物理法則を満たす予測を保証
- ベイズ最適化:フィルタのLCR値を少ないシミュレーション回数で最適化。EMC試験のpass/failを目的関数にした実験計画法
ただし現時点では、これらはあくまでフルシミュレーションの「加速器」であり、FDTDやFEMを完全に置き換えるものではない。
トラブルシューティング
よくある失敗と対策
EMI解析で初心者がやりがちな失敗ってありますか?
山ほどあるよ。よくある失敗トップ5を挙げよう:
| 失敗パターン | 原因 | 対策 |
|---|---|---|
| シミュレーションと実測で20dB以上の乖離 | 寄生素子(ESL、リードインダクタンス、寄生容量)の未モデル化 | 部品のSPICEモデルを実測Sパラメータベースのものに差し替え。PCBトレースのL/C/Rを電磁界抽出 |
| フィルタ追加でノイズが増加 | グランドラインのインダクタンスとYコンデンサの共振 | フィルタ入出力のグランドパスをフルモデルに含める。インピーダンスアナライザで検証 |
| 放射EMIが特定周波数だけ突出 | ケーブルまたはPCBトレースの共振($\lambda/4$, $\lambda/2$) | 共振周波数から配線長を逆算し、フェライトコアの設置やケーブル長の変更を検討 |
| FDTDの結果が発散 | CFL条件違反(時間刻みが大きすぎ)。または材料パラメータのミス | メッシュサイズと時間刻みの整合を確認。$\Delta t < \Delta x / (c\sqrt{3})$ |
| 計算時間が膨大で実用にならない | 全体を均一な細かいメッシュで切っている | 非均一メッシュ、サブグリッディング。遠方界はNTFF(Near-to-Far-Field変換)で計算 |
シミュレーションと実測の乖離
先生、シミュレーション結果と実測が合わないとき、まず何をチェックすればいいですか?
まず深呼吸してから、以下の順番で確認する:
- LISN接続の確認:シミュレーションモデルにLISNが正しく含まれているか。LISN未接続で計算して「合わない」は初心者あるある
- 寄生容量の見直し:特にMOSFETパッケージの端子間容量($C_{oss}$, $C_{iss}$, $C_{rss}$)とヒートシンクへの寄生容量。これがCMノイズの主因
- ケーブルハーネスのモデル:実際のワイヤハーネスの長さ、束ね方、シャーシからの距離が放射特性を大きく左右する
- 試験環境の再現:EMC暗室のグランドプレーン、DUTの設置高さ(通常5cm)、アンテナ距離(1mまたは3m)をモデルに含めているか
- 検波器の違い:シミュレーションのFFTはピーク値だが、CISPR規格は準尖頭値(QP)検波器を使用。準尖頭値はピーク値より最大13dB低いため、補正が必要
5番目の検波器の違いは盲点でした。ピーク値で余裕があっても、準尖頭値検波のことを忘れてたら楽観的すぎたりもしそうですね。
逆だね。準尖頭値はピーク値より常に低い(同じか低い)ので、ピーク値でNGでもQPだとOKになるケースがある。ただし低い繰り返し率のバースト的ノイズではピーク値とQPの差が大きくなるから、シミュレーションのピーク値で6dBマージンを確保しておくのが安全だ。
「シミュと実測の差」を詰めるのは3回かかる
EMI解析の世界では「1回目のシミュレーションは外れて当然」という暗黙の了解がある。ベンダーのデモでは綺麗な結果が出るが、自社製品に適用するとモデルの精度が足りず大きな乖離が出る。原因は、配線のインダクタンス、パッケージの内部ボンドワイヤ、筐体の隙間からの放射など、最初のモデルでは省略されがちな要素だ。「シミュ→実測→モデル修正」を3〜4サイクル繰り返して初めてシミュレーション結果が信頼できるようになる。最初から完璧を目指すのではなく、反復改善を前提としたスケジュールを組むのが現実的だ。
各ソフトの成り立ちって、結構ドラマチックだったりしますか?
Ansys Maxwell
「Ansys Maxwell」について教えてください!
Ansoft Maxwell。低周波電磁場解析。2008年Ansysに統合。
現在の所属: Ansys Inc.
Ansys HFSS
次はAnsys HFSSの話ですね。どんな内容ですか?
Ansoft Corporationが開発した3D高周波電磁界シミュレータ。2008年にAnsysがAnsoftを買収。
現在の所属: Ansys Inc.
COMSOL Multiphysics
「COMSOL Multiphysics」について教えてください!
1986年スウェーデンで設立。MATLAB連携のFEMLABとして開始、後にCOMSOLに改名。マルチフィジックスに強み。
現在の所属: COMSOL AB
なるほど。じゃあ低周波電磁場解析ができていれば、まずは大丈夫ってことですか?
ファイル形式と相互運用性
異なるソフト間でデータを受け渡しするときの注意点ってありますか?
| フォーマット | 拡張子 | 種別 | 概要 |
|---|---|---|---|
| STEP | .stp/.step | 中立CAD | ISO 10303準拠の3D CADデータ交換フォーマット。形状+PMI対応。 |
| IGES | .igs/.iges | 中立CAD | 初期のCADデータ交換規格。曲面データの互換性に課題あり。STEPへの移行が進む。 |
| VTK | .vtk/.vtu | 可視化 | Visualization Toolkit形式。ParaView等で使用。 |
異なるソルバー間でモデルを変換する際は、要素タイプの対応関係、材料モデルの互換性、荷重・境界条件の表現差異に注意が必要になるんだ。特に高次要素や特殊要素(コヒーシブ要素、ユーザー定義要素等)はソルバー間で直接変換できない場合が多い。
なるほど…フォーマットって一見シンプルだけど、実はすごく奥が深いんですね。
実務上の注意点
教科書には載ってない「現場の知恵」みたいなものってありますか?
メッシュ収束性の確認、境界条件の妥当性検証、材料パラメータの感度分析がすごく大事なんだ。
いやぁ、スイッチングノイズEMI解析って奥が深いですね… でも先生の説明のおかげでだいぶ整理できました!
うん、いい調子だよ! 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。
スイッチングノイズの「コモンモードとノーマルモード」——この区別がEMI対策の入り口
EMI対策を始めるとき、必ず最初に問われるのがノイズの種類だ。ノーマルモードノイズは電源ラインの往復2線を逆向きに流れ、ラインフィルタ(差動チョーク)で対処する。コモンモードノイズは2線が同じ方向に流れてシャーシとの間を流れ、コモンモードチョークとYコンデンサで対策する。スイッチング電源で多いのは、主回路のdV/dtがヒートシンク(シャーシ接続)との寄生容量を通じてコモンモード電流を流すケース——「なぜかシャーシが熱い」の原因がこれだったりする。理論上の発生メカニズムを押さえないと、対策部品を闇雲に追加しても解決しない。
各項の物理的意味
- 電場項 $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$:ファラデーの電磁誘導法則。時間変動する磁束密度が起電力を生じさせる。【日常の例】自転車のダイナモ(発電機)は、磁石を回転させることで近くのコイルに電圧が発生する——磁場が時間的に変化すると電場が誘起されるというこの法則の直接的応用。IHクッキングヒーターも同じ原理で、高周波磁場の変化が鍋底に渦電流を誘起し、ジュール熱で加熱する。
- 磁場項 $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$:アンペア-マクスウェルの法則。電流と変位電流が磁場を生成する。【日常の例】電線に電流を流すと周囲に磁場が生じる——これがアンペアの法則。電磁石はこの原理で動作し、コイルに電流を流して強力な磁場を作る。スマートフォンのスピーカーも、電流→磁場→振動板の力というこの法則の応用。高周波(GHz帯のアンテナ等)では変位電流 $\partial D/\partial t$ が無視できなくなり、電磁波の放射を記述する。
- ガウスの法則 $\nabla \cdot \mathbf{D} = \rho_v$:電荷が電束の発散源であることを示す。【日常の例】下敷きで髪の毛をこすると静電気で髪が逆立つ——帯電した下敷き(電荷)から電気力線が放射状に広がり、軽い髪の毛に力を及ぼす。コンデンサ(キャパシタ)の設計では、電極間の電場分布をこの法則で計算する。ESD(静電気放電)対策もガウスの法則に基づく電場解析が基盤。
- 磁束保存 $\nabla \cdot \mathbf{B} = 0$:磁気単極子が存在しないことを表す。【日常の例】棒磁石を半分に割っても、N極だけ・S極だけの磁石は作れない——必ずN極とS極がペアで存在する。これは磁力線が「始点も終点もない閉じたループ」を描くことを意味する。数値解析では、この条件を満たすためにベクトルポテンシャル $\mathbf{B} = \nabla \times \mathbf{A}$ という定式化を用い、磁束保存を自動的に保証する。
仮定条件と適用限界
- 線形材料仮定:透磁率・誘電率が磁場・電場強度に依存しない(飽和領域では非線形B-Hカーブが必要)
- 準静的近似(低周波):変位電流項を無視可能($\omega \varepsilon \ll \sigma$)。渦電流解析で一般的
- 2D仮定(断面解析):電流方向が一様で、端部効果を無視できる場合に有効
- 等方性仮定:異方性材料(珪素鋼板の圧延方向等)では方向別の特性定義が必要
- 適用外ケース:プラズマ(電離気体)、超伝導体、非線形光学材料では追加の構成則が必要
数値解法と実装
数値手法の詳細
具体的にはどんなアルゴリズムでスイッチングノイズEMI解析を解くんですか?
ここまで聞いて、スイッチングノイズがなぜ重要か、やっと腹落ちしました!
離散化の定式化
形状関数 $N_i$ を用いて未知量を近似:
これを数式で表すとこうなるよ。
基礎方程式の離散形
これを数式で表すとこうなるよ。
うーん、式だけだとピンとこないです… 何を表してるんですか?
連続体の支配方程式を離散化すると、以下の代数方程式系が得られる:
ここで $[K]$ は全体剛性マトリクス(または同等のシステムマトリクス)、$\{u\}$ は未知節点変数ベクトル、$\{F\}$ は外力ベクトルなんだ。
あっ、そういうことか! 連続体の支配方程式をってそういう仕組みだったんですね。
要素技術
「要素技術」って聞いたことはあるんですけど、ちゃんと理解できてないかもしれません…
| 要素タイプ | 次数 | 節点数(3D) | 精度 | 計算コスト |
|---|---|---|---|---|
| 四面体1次 | 線形 | 4 | 低(シアロッキング) | 低 |
| 四面体2次 | 二次 | 10 | 高 | 中 |
| 六面体1次 | 線形 | 8 | 中 | 中 |
| 六面体2次 | 二次 | 20 | 非常に高 | 高 |
| プリズム | 線形/二次 | 6/15 | 中〜高 | 中 |
積分スキーム
積分スキームって、具体的にはどういうことですか?
ここまで聞いて、要素タイプがなぜ重要か、やっと腹落ちしました!
収束性と安定性
収束しなくなったら、まず何をチェックすればいいですか?
収束速度: 二次要素で $O(h^2)$ のオーダーで誤差が減少(滑らかな解の場合)
なるほど…メッシュを細分化って一見シンプルだけど、実はすごく奥が深いんですね。
ソルバー設定の推奨事項
具体的にはどんなアルゴリズムでスイッチングノイズEMI解析を解くんですか?
| パラメータ | 推奨値 | 備考 |
|---|---|---|
| 反復法の収束判定 | $10^{-6}$ | 残差ノルム基準 |
| 前処理手法 | ILU(0) or AMG | 問題規模による |
| 最大反復回数 | 1000 | 非収束時は設定見直し |
| メモリモード | In-core | 可能な限り |
辺要素(Nedelec要素)
電磁場解析に特化した要素。接線成分の連続性を自動的に保証し、スプリアスモードを排除。3D高周波解析の標準。
節点要素
スカラーポテンシャル定式化に使用。静磁場のスカラーポテンシャル法や静電場解析で有効。
FEM vs BEM(境界要素法)
FEM: 非線形材料・非均質媒質に対応。BEM: 無限領域(開領域問題)を自然に扱える。ハイブリッドFEM-BEMも有効。
非線形収束(磁気飽和)
B-Hカーブの非線形性をニュートン・ラフソン法で処理。残差基準: $||R||/||R_0|| < 10^{-4}$が一般的。
周波数領域解析
時間高調波仮定により定常問題に帰着。複素数演算が必要だが、広帯域特性は時間領域解析で取得。
時間領域の時間刻み
最高周波数成分の1/20以下の時間刻みが必要。暗黙的時間積分ではより大きな刻みも可能だが精度に注意。
周波数領域と時間領域の使い分け
周波数領域解析は「ラジオの特定の周波数に合わせる」ようなもの——1つの周波数での応答を効率的に計算できる。時間領域解析は「全チャンネルを同時に録画する」ようなもの——あらゆる周波数成分を含む過渡現象を再現できるが計算コストが高い。
実践ガイド
実践ガイド
先生、「実践ガイド」について教えてください!
スイッチングノイズEMI解析の実務的な解析フローと注意点を解説する。
ここまで聞いて、スイッチングノイズがなぜ重要か、やっと腹落ちしました!
解析フロー
最初の一歩から教えてください! 何から始めればいいですか?
1. 前処理 (Pre-processing)
- CADデータのインポートと形状簡略化
- 材料特性の定義
- メッシュ生成(要素タイプ・サイズの決定)
- 境界条件と荷重条件の設定
2. 求解 (Solving)
- ソルバー設定(解法、収束基準、出力制御)
- ジョブ投入と計算実行
- 収束モニタリング
3. 後処理 (Post-processing)
- 結果の可視化(変位、応力、その他の物理量)
- 結果の検証と妥当性確認
- レポート作成
メッシュ生成のベストプラクティス
メッシュの良し悪しってどうやって判断するんですか?
要素品質指標
「要素品質指標」について教えてください!
| 指標 | 理想値 | 許容範囲 | 影響 |
|---|---|---|---|
| アスペクト比 | 1.0 | < 5.0 | 精度低下 |
| ヤコビアン比 | 1.0 | > 0.3 | 要素退化 |
| ワーピング | 0° | < 15° | 精度低下 |
| スキューネス | 0° | < 45° | 収束性悪化 |
| テーパー比 | 0 | < 0.5 | 精度低下 |
メッシュ密度の決定
メッシュ密度の決定って、具体的にはどういうことですか?
境界条件の設定指針
境界条件って、ここを間違えると全部ダメになるって聞いたんですけど…
あっ、そういうことか! 過拘束に注意ってそういう仕組みだったんですね。
商用ツール別の実装手順
いろんなソフトがあるんですよね? それぞれの特徴を教えてください!
| ツール名 | 開発元/現在 | 主要ファイル形式 |
|---|---|---|
| Ansys Maxwell | Ansys Inc. | .aedt, .maxwell |
| Ansys HFSS | Ansys Inc. | .aedt, .hfss |
| COMSOL Multiphysics | COMSOL AB | .mph |
| CST Studio Suite | Dassault Systèmes SIMULIA | .cst |
Ansys Maxwell
「Ansys Maxwell」について教えてください!
Ansoft Maxwell。低周波電磁場解析。2008年Ansysに統合。
現在の所属: Ansys Inc.
Ansys HFSS
次はAnsys HFSSの話ですね。どんな内容ですか?
Ansoft Corporationが開発した3D高周波電磁界シミュレータ。2008年にAnsysがAnsoftを買収。
現在の所属: Ansys Inc.
先生の説明分かりやすい! ツール名のモヤモヤが晴れました。
よくある失敗と対策
初心者がやりがちな失敗パターンってありますか? 事前に知っておきたいです!
| 症状 | 原因 | 対策 |
|---|---|---|
| 計算が収束しない | メッシュ品質不良、不適切な境界条件 | メッシュ改善、拘束条件見直し |
| 応力が異常に大きい | 応力特異点、メッシュ依存 | 特異点回避、局所メッシュ細分化 |
| 変位が非現実的 | 材料定数誤り、単位系不整合 | 入力データ確認 |
| 計算時間が過大 | 不要な細分化、非効率な解法 | メッシュ最適化、並列計算 |
品質保証チェックリスト
教科書には載ってない「現場の知恵」みたいなものってありますか?
いやぁ、スイッチングノイズEMI解析って奥が深いですね… でも先生の説明のおかげでだいぶ整理できました!
うん、いい調子だよ! 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。
EMIフィルタの「グランド」問題——接地が悪いと逆効果になる現実
EMIフィルタを付けたのにノイズが減らない、むしろ増えた——という経験をした人は少なくない。原因のひとつがグランドラインのインピーダンスだ。コモンモードフィルタのYコンデンサはシャーシグランドに接続するが、グランドラインが長かったり薄かったりすると、コンデンサのリード線インダクタンスと共振してノイズを増幅してしまう。実務では「フィルタはできるだけシャーシの近くに」「グランドはワイドなパターンで短く」が鉄則。解析で接地ラインを含めたフルモデルを組まないと、こういう共振は見えてこない。
解析フローのたとえ
モータの電磁界解析は「ギターの調律」に近い感覚です。弦の太さ(コイル巻数)とブリッジの位置(磁石配置)を調整して、最も美しい音色(効率の良いトルク特性)を引き出す。1つのパラメータを変えると全体のバランスが変わる——だからパラメトリックスタディが重要なんです。
初心者が陥りやすい落とし穴
「空気領域? なんで空気をメッシュで切るの?」——初めて電磁界解析に触れた人がほぼ全員抱く疑問です。答えは「磁力線は鉄心の外にも広がるから」。解析領域を鉄心ぎりぎりにすると、行き場を失った磁束が壁に「ぶつかって」反射し、実際にはありえない磁束集中が起きます。部屋が狭すぎてボールが壁に跳ね返りまくる状態を想像してみてください。
境界条件の考え方
遠方の境界条件って地味ですが超重要です。「ここから先は無限に広がる空間」ということを数値的に表現する必要がある。設定を間違えると、まるで「見えない壁」があるかのように磁束が跳ね返されてしまいます。
ソフトウェア比較
商用ツール比較
いろんなソフトがあるんですよね? それぞれの特徴を教えてください!
スイッチングノイズEMI解析に対応する主要な商用CAEツールの機能比較と、各製品の歴史的背景を詳述する。
ここまで聞いて、スイッチングノイズがなぜ重要か、やっと腹落ちしました!
対応ツール一覧
で、スイッチングノイズEMI解析をやるにはどんなソフトが使えるんですか?
| ツール名 | 開発元/現在 | 主要ファイル形式 |
|---|---|---|
| Ansys Maxwell | Ansys Inc. | .aedt, .maxwell |
| Ansys HFSS | Ansys Inc. | .aedt, .hfss |
| COMSOL Multiphysics | COMSOL AB | .mph |
| CST Studio Suite | Dassault Systèmes SIMULIA | .cst |
Ansys Maxwell
「Ansys Maxwell」について教えてください!
Ansoft Maxwell。低周波電磁場解析。2008年Ansysに統合。
現在の所属: Ansys Inc.
Ansys HFSS
次はAnsys HFSSの話ですね。どんな内容ですか?
Ansoft Corporationが開発した3D高周波電磁界シミュレータ。2008年にAnsysがAnsoftを買収。
現在の所属: Ansys Inc.
COMSOL Multiphysics
「COMSOL Multiphysics」について教えてください!
1986年スウェーデンで設立。MATLAB連携のFEMLABとして開始、後にCOMSOLに改名。マルチフィジックスに強み。
現在の所属: COMSOL AB
CST Studio Suite
CST Studioって、具体的にはどういうことですか?
Computer Simulation Technology (ドイツ) が開発。2016年にDassault Systèmesが買収しSIMULIAに統合。
現在の所属: Dassault Systèmes SIMULIA
先輩が「低周波電磁場解析だけはちゃんとやれ」って言ってた意味が分かりました。
機能比較マトリクス
予算も時間も限られてるんですけど、コスパ最強はどれですか?
| 機能 | Maxwell | HFSS | COMSOL | CST |
|---|---|---|---|---|
| 基本機能 | ○ | ○ | ○ | ○ |
| 高度な機能 | ○ | ○ | ○ | △ |
| 自動化/スクリプト | ○ | ○ | ○ | ○ |
| 並列計算 | ○ | ○ | ○ | ○ |
| GPU対応 | △ | △ | △ | ○ |
変換時のリスク
変換時のリスクって、具体的にはどういうことですか?
あっ、そういうことか! 異なるツール間でのモってそういう仕組みだったんですね。
ライセンス形態
「ライセンス形態」って聞いたことはあるんですけど、ちゃんと理解できてないかもしれません…
| ツール | ライセンス | 特徴 |
|---|---|---|
| 商用FEA | ノードロック/フローティング | 高額だが公式サポート付き |
| OpenFOAM | GPL | 無償だがサポートは有償 |
| COMSOL | ノードロック/フローティング | モジュール単位で購入 |
| Code_Aster | GPL | EDF開発のOSSソルバー |
選定の指針
結局どれを選べばいいか、判断基準を教えてもらえますか?
スイッチングノイズEMI解析のツール選定においては以下を考慮:
いやぁ、スイッチングノイズEMI解析って奥が深いですね… でも先生の説明のおかげでだいぶ整理できました!
うん、いい調子だよ! 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。