ESD(静電気放電)シミュレーション
理論と物理
ESDとは何か
静電気放電のシミュレーションって何をモデル化するんですか? バチッとくるあれですか?
そう、まさにあの「バチッ」だよ。ただし、人間が感じるのは3kV以上だけど、半導体ICは数百Vで壊れる。しかもESDの放電パルスは立ち上がり時間がわずか0.7〜1ナノ秒で、ピーク電流は数アンペアから数十アンペアに達する。雷の1000倍速い衝撃波がICの酸化膜を数ナノメートルの経路で絶縁破壊させるんだ。
えっ、ナノ秒って…。そんな速い現象をシミュレーションで再現できるんですか?
できる。人体モデル(HBM)で4kVのESD電流は立ち上がり1nsでピーク7.5Aになる。この超高速パルスがPCBパターンを伝搬して、ICの入力保護回路に到達するまでの電圧・電流分布を解析するのがESDシミュレーションだ。IEC 61000-4-2のレベル4(±8kV接触放電)をクリアするには、ESDパスのインピーダンスを3D FEMで最適化する必要がある。
具体的にはどういう場面で使うんですか?
典型的なのはスマートフォンのUSB-Cコネクタだね。ユーザーが帯電した状態でケーブルを挿すと、コネクタのシェル→GNDプレーン→ICの経路でESD電流が流れる。この経路のインピーダンスが高いと、ICのVSS-VDDピン間に規格外の電圧がかかって誤動作やラッチアップが起きる。出荷前にシミュレーションで電流密度マップを確認しておけば、量産後のリコールを防げるわけだ。
ESD試験モデル(HBM・MM・CDM・IEC)
HBMとかIECとか、いろんなモデルがあるって聞いたんですけど、どう違うんですか?
ESD試験には対象レベルに応じて複数のモデルがある。それぞれ回路定数が違うから、電流波形の形が全然違うんだ。
| モデル | 規格 | C [pF] | R [Ω] | 立ち上がり | 対象レベル |
|---|---|---|---|---|---|
| HBM | ANSI/ESDA/JEDEC JS-001 | 100 | 1,500 | 2〜10 ns | デバイス(ICチップ) |
| MM | EIAJ ED-4701 | 200 | 0(実質 <10) | <1 ns | デバイス(振動的波形) |
| CDM | ANSI/ESDA/JEDEC JS-002 | デバイス依存 | — | 100〜400 ps | デバイス(実装状態) |
| IEC | IEC 61000-4-2 | 150 | 330 | 0.7〜1 ns | システム(製品全体) |
CDMの立ち上がりが100ピコ秒って…。シミュレーションで再現するのすごく大変そうですね。
その通り。CDMは時間刻みをピコ秒オーダーにする必要があるから計算コストが桁違いに高い。だから実務ではまずIECモデルでシステムレベルの合否判定をして、NGならHBMモデルでICの入力保護回路を詳細に解析する、という2段階アプローチが一般的だよ。
支配方程式
ESD解析の「土台」になる方程式って何ですか?
基本はマクスウェル方程式だ。ESDパルスは周波数成分がDC〜数GHzに広がるから、準静的近似は使えない。フルウェーブ(全波)解析が必要になる。
マクスウェル方程式(時間領域)
$$\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}$$ $$\nabla \times \mathbf{H} = \mathbf{J} + \frac{\partial \mathbf{D}}{\partial t}$$ $$\nabla \cdot \mathbf{D} = \rho_v$$ $$\nabla \cdot \mathbf{B} = 0$$ESDの場合、特に重要なのは変位電流項 $\partial \mathbf{D}/\partial t$ だ。ESDパルスの $di/dt$ が数GA/sに達するから、磁場の変動も激しく、EMI(電磁干渉)として周囲の回路にノイズを撒き散らす。この「撒き散らし」を正確に追うには、マクスウェル方程式を省略なしで解かなきゃいけない。
HBM電流波形の数理モデル
HBMの波形って、どんな数式で表現するんですか?
HBMの等価回路は100pFのキャパシタを1.5kΩの抵抗を通して放電するRCモデルだ。ただし、放電経路には寄生インダクタンス $L_p$ があるから、実際の波形は二重指数関数になる。
HBM電流波形(二重指数関数モデル)
$$i_{\text{HBM}}(t) = \frac{I_p}{1 - k}\left(e^{-t/\tau_1} - e^{-t/\tau_2}\right)$$ここで:
- $I_p$:ピーク電流 [A](4kV HBMで約2.67A)
- $\tau_1 = RC = 150\,\text{ns}$:減衰時定数
- $\tau_2 = L_p / R \approx 2\text{--}10\,\text{ns}$:立ち上がり時定数
- $k = \tau_2 / \tau_1$:時定数比
$\tau_1$ と $\tau_2$ の差が大きいから、急激に立ち上がってゆっくり減衰する形になるんですね。
その通り。あとESDが蓄えるエネルギーも重要だ。
ESD蓄積エネルギー
$$E_{\text{ESD}} = \frac{1}{2}CV^2$$例えば IEC モデル($C = 150\,\text{pF}$, $V = 8\,\text{kV}$)では:
$$E_{\text{ESD}} = \frac{1}{2} \times 150 \times 10^{-12} \times (8000)^2 = 4.8\,\text{mJ}$$4.8mJは微小に見えるけど、これがナノ秒で放出されると瞬時電力は数キロワットになる。ICの入力保護回路(ESDクランプ)がこのエネルギーを安全に吸収できるかどうかが設計の勝負どころだ。
IEC 61000-4-2 電流源モデル
IECの波形はHBMとどう違うんですか?
IEC 61000-4-2の波形は2つのピークを持つ独特な形をしている。最初の鋭いスパイクと、その後のブロードなハンプだ。規格ではこの波形を4つのパラメータで定義している。
| パラメータ | レベル4(±8kV接触) | 許容範囲 |
|---|---|---|
| 初期ピーク電流 $I_1$ | 30 A | ±15% |
| $t = 30\,\text{ns}$ での電流 $I_2$ | 16 A | ±30% |
| $t = 60\,\text{ns}$ での電流 $I_3$ | 8 A | ±30% |
| 立ち上がり時間 $t_r$ | 0.7〜1 ns | — |
8kVで初期ピーク30A! HBMの約10倍じゃないですか。
そう。IECはシステムレベルの試験だから、人体の手→筐体→PCB→ICという現実的な放電経路全体を模擬している。シミュレーションでは、この波形を電流源として放電ポイントに注入し、PCB上の応答を解析する。等価回路としては次のような2段RLC回路で近似できる。
IEC波形の等価回路近似
$$i_{\text{IEC}}(t) = I_1 \cdot A_1 \cdot e^{-t/\tau_a}\sin(\omega_1 t) + I_2 \cdot A_2 \cdot \left(e^{-t/\tau_b} - e^{-t/\tau_c}\right)$$第1項:初期スパイク($\tau_a \approx 2\,\text{ns}$, $\omega_1$ は寄生Lに依存)
第2項:ブロードハンプ($\tau_b \approx 20\,\text{ns}$, $\tau_c \approx 2\,\text{ns}$)
グラウンドバウンスの物理
グラウンドバウンスってよく聞くんですけど、ESDでなぜ問題になるんですか?
グラウンドプレーンは「完全なゼロ電位の面」じゃない。実際にはプレーンにもインダクタンス $L_{\text{gnd}}$ がある。ESD電流が流れると、このインダクタンスによって局所的な電位上昇が起きる。これがグラウンドバウンスだ。
グラウンドバウンス電圧
$$V_{\text{bounce}} = L_{\text{gnd}} \cdot \frac{di}{dt}$$IECレベル4の初期スパイクで $di/dt$ を概算すると:
$$\frac{di}{dt} \approx \frac{30\,\text{A}}{1\,\text{ns}} = 30 \times 10^9\,\text{A/s} = 30\,\text{GA/s}$$もし $L_{\text{gnd}} = 1\,\text{nH}$(ビア1本分程度)なら:
$$V_{\text{bounce}} = 1 \times 10^{-9} \times 30 \times 10^9 = 30\,\text{V}$$ビア1本で30V!? グラウンドがそんなに暴れるんですか…。
だからESD対策ではGNDビアを放電点の直近に複数本配置して $L_{\text{gnd}}$ を下げることが鉄則なんだ。例えばビアを4本並列にすれば $L_{\text{gnd}}$ は約1/4になり、バウンス電圧も7.5Vに下がる。シミュレーションでは、GNDプレーンの電位分布をカラーマップで可視化して、バウンスが大きい「ホットスポット」を特定するのが定番の使い方だよ。
TVSクランプ電圧モデル
TVSダイオードの「クランプ電圧」って何ですか? ESD対策で必ず出てきますよね。
TVS(Transient Voltage Suppressor)はESD電流を吸収するための保護素子だ。正常時は高インピーダンスで信号に影響しないけど、ESDが印加されるとブレークダウンして低インピーダンスになり、電圧を一定値に「クランプ」する。
TVSクランプ電圧の関係式
$$V_{\text{clamp}} = V_{\text{BR}} + I_{\text{ESD}} \cdot R_{\text{dyn}}$$ここで:
- $V_{\text{BR}}$:ブレークダウン電圧(TVSの設計値、例:5V品なら6.4V typ.)
- $R_{\text{dyn}}$:動的抵抗(0.1〜2Ω程度。小さいほど高性能)
- $I_{\text{ESD}}$:TVSに流れるESD電流 [A]
例:$V_{\text{BR}} = 6.4\,\text{V}$, $R_{\text{dyn}} = 0.5\,\Omega$, $I_{\text{ESD}} = 16\,\text{A}$ の場合
$$V_{\text{clamp}} = 6.4 + 16 \times 0.5 = 14.4\,\text{V}$$ICの絶対最大定格が例えば7Vだったら、14.4Vだとアウトですよね…。
まさにそこがESDシミュレーションの核心だ。クランプ電圧が高すぎるなら、$R_{\text{dyn}}$ が小さいTVSに変更するか、TVSの配置を放電点に近づけて配線のインダクタンスによる電圧上昇分を減らす。シミュレーションでは、TVSのSPICEモデルと3D電磁場解析を連成させて、ICピンでの実効クランプ電圧を算出するのが正攻法だ。
もう1つ重要なのが、TVS〜ICピン間の配線インダクタンス $L_{\text{trace}}$ の影響だ。
ICピンでの実効電圧
$$V_{\text{IC}} = V_{\text{clamp}} + L_{\text{trace}} \cdot \frac{di}{dt}$$$L_{\text{trace}} = 2\,\text{nH}$(約3mm配線相当)、$di/dt = 30\,\text{GA/s}$ の場合:
$$V_{\text{IC}} = 14.4 + 2 \times 10^{-9} \times 30 \times 10^9 = 14.4 + 60 = 74.4\,\text{V}$$74V!? たった3mmの配線で60Vも上乗せされるんですか!
だからESD対策の世界では「TVSはICの隣に置け、1mmでも近く」が鉄則なんだ。ここまで来ると、回路図上の接続だけでは対策できない。PCBレイアウトの3D電磁場解析が不可欠になる理由がわかるだろう?
ESDの放電時間は1ナノ秒——雷より1000倍速い衝撃
冬にドアノブに触れて「バチッ」とくる静電気は約3〜5kV。この放電が起きている時間はわずか数ナノ秒で、雷の放電時間(数マイクロ秒)の約1000分の1だ。1960年代、IBMのエンジニアが半導体工場で不良率が冬に急増する原因を調査し、静電気が原因であることを突き止めた。そこで考案されたのが人体モデル(HBM)——人体を100pFのコンデンサと1.5kΩの抵抗でモデル化するという、今なお現役の規格である。もし人体の静電容量がもっと大きかったら、ESD対策のコストは跳ね上がっていただろう。
各項の物理的意味
- HBM二重指数関数の第1項 $e^{-t/\tau_1}$:RCの放電減衰を表す。$\tau_1 = RC = 100\,\text{pF} \times 1.5\,\text{k}\Omega = 150\,\text{ns}$ で、放電パルスの「尾」の長さを決める。この時定数が短いほどエネルギーが速く放出され、ICへの熱ストレスが増大する。
- HBM二重指数関数の第2項 $e^{-t/\tau_2}$:放電路の寄生インダクタンス $L_p$ による立ち上がりを表す。$\tau_2 = L_p / R$ で、$L_p$ が小さいほど立ち上がりが急峻になる。CDMでは $L_p$ が極めて小さいため、100ps台の超高速立ち上がりとなる。
- グラウンドバウンス $V = L \cdot di/dt$:ファラデーの電磁誘導法則の回路版。自己インダクタンスによる逆起電力であり、ESD電流の $di/dt$ が数十GA/sに達するため、わずかなnHのインダクタンスでも数十Vの電圧が発生する。
- TVS動的抵抗 $R_{\text{dyn}}$:TVSのI-V特性のブレークダウン領域における微分抵抗。データシートの「クランプ電圧 at $I_{pp}$」から逆算可能。小信号用TVSは $R_{\text{dyn}}$ が大きいため、大電流ESDでは過大なクランプ電圧になることがある。
仮定条件と適用限界
- HBM・IEC波形モデルは集中定数近似:放電路長がESD波形の波長($c / f_{\text{max}} \approx 0.1\,\text{m}$)より十分短い場合に有効。PCBが大型の場合は分布定数モデルが必要
- TVSの動的抵抗 $R_{\text{dyn}}$ は定常近似:実際にはナノ秒領域でのターンオン遅延(0.5〜1ns)があり、その間TVSは高インピーダンス。この遅延中のオーバーシュートが問題になることがある
- GNDプレーンのインダクタンス推定はDC近似:GHz帯域では表皮効果や共振モードが効き、単純な $L \cdot di/dt$ では不十分。フルウェーブ3D解析が必要
- 線形媒質仮定:アーク放電によるプラズマ生成は非線形現象であり、放電ギャップ近傍では追加モデルが必要
次元解析と単位系
| 物理量 | SI単位 | ESD特有の注意点 |
|---|---|---|
| ESD電圧 | V | kV表記が多い。8kV = 8000V。入力ミスに注意 |
| ピーク電流 | A | IEC 8kVで30A。HBM 4kVで2.67A |
| 立ち上がり時間 | s | ns表記。1ns = 10⁻⁹s。周波数換算 $f \approx 0.35/t_r$ |
| 寄生L | H | nH表記。PCBビア1本 ≈ 0.5〜1nH |
| 寄生C | F | pF表記。TVSの寄生容量は0.1〜数pF(高速信号で問題) |
| $di/dt$ | A/s | GA/s表記。IEC初期スパイクで30GA/s程度 |
数値解法と実装
FDTD法によるESD過渡解析
ESD解析で一番よく使われる数値手法って何ですか?
ESD解析の主力はFDTD法(時間領域有限差分法)だ。理由は3つある。
- ESDパルスは広帯域(DC〜数GHz)なので、1回の時間領域シミュレーションで全周波数成分の応答が得られる
- 陽解法なので大規模3Dモデルでもメモリ効率がいい
- CST Studio Suiteの「Transient Solver」がFDTDベースで、ESD専用テンプレートを持っている
FDTDの時間刻みはどう決めるんですか?
FDTD法の安定条件はCFL(Courant-Friedrichs-Lewy)条件で決まる。
CFL安定条件(3D)
$$\Delta t \leq \frac{1}{c \cdot \sqrt{\frac{1}{\Delta x^2} + \frac{1}{\Delta y^2} + \frac{1}{\Delta z^2}}}$$等間隔メッシュ($\Delta x = \Delta y = \Delta z = \Delta$)の場合:
$$\Delta t \leq \frac{\Delta}{c\sqrt{3}}$$ESD解析では最高周波数 $f_{\max} \approx 0.35 / t_r \approx 5\,\text{GHz}$($t_r = 0.7\,\text{ns}$ の場合)。
$\lambda_{\min} = c / f_{\max} \approx 60\,\text{mm}$。メッシュは $\lambda_{\min} / 20 = 3\,\text{mm}$ 以下が目安。
3mmメッシュだと、PCBの微細パターン(0.15mmピッチとか)は全然足りないですよね?
いい指摘だ。だからCST Studio SuiteではPBA(Perfect Boundary Approximation)やTST(Thin Sheet Technology)を使って、メッシュ境界上の薄い導体を正確に扱う。また、ESD電流が集中する経路(TVS周辺、GNDビア列など)だけ局所的にメッシュを細分化するサブグリッド技術も重要だよ。
3D FEMによるESDパス最適化
FDTDじゃなくてFEMを使う場面もあるんですか?
ある。特にESD保護回路の3Dインピーダンス最適化ではFEMが威力を発揮する。FEMは非構造メッシュ(四面体)が使えるから、複雑な形状のコネクタやBGAパッケージの内部構造を忠実にモデル化できる。Ansys HFSSの時間領域解析やCOMSOL Multiphysicsが代表的だ。
FEMでESD経路のインピーダンスを求める場合、弱形式は次のようになる。
電磁場FEMの弱形式(辺要素 / Nedelec要素)
$$\int_{\Omega}\left(\frac{1}{\mu}\nabla \times \mathbf{N}_i \cdot \nabla \times \mathbf{E} + \epsilon\,\mathbf{N}_i \cdot \frac{\partial^2 \mathbf{E}}{\partial t^2} + \sigma\,\mathbf{N}_i \cdot \frac{\partial \mathbf{E}}{\partial t}\right)d\Omega = -\int_{\Omega}\mathbf{N}_i \cdot \frac{\partial \mathbf{J}_s}{\partial t}\,d\Omega$$$\mathbf{N}_i$:辺要素の形状関数(接線成分の連続性を自動保証)
$\mathbf{J}_s$:ESD電流源項
辺要素って、電磁場専用の要素なんですか?
そう。通常の節点要素でベクトル場(電場・磁場)を扱うと「スプリアスモード」という偽の固有モードが出てしまう。辺要素(Nedelec要素)は要素の辺に沿った接線成分を自由度とするから、この問題を根本的に回避できる。高周波電磁場解析では事実上の標準だ。
SPICE連成解析
TVSダイオードの非線形特性も一緒に解くんですよね? それってどうやるんですか?
TVSの非線形I-V特性はマクスウェル方程式の中では直接扱えない。そこで3D電磁場解析(FDTD/FEM)とSPICE回路シミュレータを連成させるアプローチをとる。CST Studio Suiteでは「Schematic」機能で、3DモデルのポートにSPICEサブ回路(TVSモデル)を接続できる。
| 解析対象 | 手法 | 役割 |
|---|---|---|
| PCBパターン・ビア・筐体 | 3D FDTD / FEM | 電磁場分布・電流密度マップ |
| TVSダイオード | SPICEモデル | 非線形クランプ特性 |
| IC入力保護回路 | IBIS / SPICEモデル | IC端子の実効電圧波形 |
| ESD電流源 | 波形定義(IEC準拠) | 放電ポイントへの電流注入 |
かなり複雑な連成解析になるんですね。計算時間はどのくらいかかるんですか?
一般的なスマートフォン基板(8層、100mm×50mm)で、IEC 8kVの100nsパルスを解析する場合、メッシュ数は数百万セル、計算時間はGPUアクセラレーション使用で30分〜2時間程度だ。GPUなしだと半日コースになることもある。だから放電ポイントを絞り込んでから走らせるのが実務的だね。
メッシュ要件と時間刻み
ESD解析特有のメッシュの注意点ってありますか?
構造解析のメッシュとは考え方がかなり違う。ESD解析では以下が重要だ。
| 項目 | 推奨値 | 理由 |
|---|---|---|
| 最小メッシュサイズ | $\leq \lambda_{\min}/20$ | 波長分解能の確保(5GHzで3mm以下) |
| 導体表面メッシュ | $\leq \delta/3$(表皮深さの1/3) | 表皮効果の正確な解像 |
| PCBビア周辺 | ビア径の1/3以下 | 電流集中部の解像 |
| TVSパッド周辺 | 0.1mm以下 | クランプ動作の正確な再現 |
| 時間刻み | CFL条件 × 0.9 | 安定性マージン |
| シミュレーション時間 | $\geq 5\tau_1$(HBM: 750ns) | パルス減衰の十分な観測 |
表皮深さとメッシュの関係
ESD電流の周波数成分は数GHzに達する。銅の表皮深さは $\delta = \sqrt{2/(\omega\mu\sigma)}$ で、1GHzでは約2μm、5GHzでは約0.9μmだ。PCBの銅箔厚が35μmだとすると、ESD電流は表面のごく薄い層にしか流れない。これを正確に解くには導体表面に極めて細かいメッシュが必要になるが、3D FDTDでは非現実的。そこでCST Studio Suiteでは「Thin Sheet」近似を使い、導体をメッシュの面として扱うことで計算コストを大幅に削減している。
実践ガイド
ESD解析フロー
実際にESD解析をやるとしたら、どういう手順で進めるんですか?
標準的なフローは以下の5ステップだ。
- 放電ポイントの特定:IEC 61000-4-2の試験規格に基づき、ユーザーがアクセス可能な金属部(コネクタ、ネジ頭、筐体の隙間)を洗い出す
- 3Dモデル構築:PCBのGerberデータ + 筐体CADを取り込み、メッシュ化。ESD電流が通過しない遠方部品は省略して計算コストを削減
- 電流源の設定:IEC 61000-4-2のレベル4波形(±8kV接触)を放電ポイントに注入。GNDリターンパスも忘れずに定義
- 保護回路のモデリング:TVSダイオード、フェライトビーズ、バイパスコンデンサのSPICEモデルをポートに接続
- 結果評価:ICピンの電圧波形が絶対最大定格以内か、GNDバウンスのホットスポットはどこか、電流密度が許容値を超えるパターンはないかを確認
PCBモデリングのポイント
PCBをモデル化するとき、何が一番難しいですか?
一番の落とし穴はGNDプレーンのスリット(切り欠き)だ。PCB設計者がアナログ/デジタルのGNDを分離するために入れたスリットが、ESD電流の帰還経路を遮断してしまうケースが非常に多い。シミュレーションで電流密度マップを見ると、スリットの端で電流が集中して迂回している様子が一目でわかる。
GNDは一枚板がベストってことですか?
ESD耐性の観点だけなら、そうだ。ただし、ノイズ分離の要求もあるから、実際にはスリットの上をESD電流が横切らないレイアウトにする。具体的には:
- TVSダイオードのGNDパッドは、スリットのない側のプレーンに接続する
- コネクタのシェルGNDとICのVSSが同じプレーン島にあることを確認する
- スリットを跨ぐ信号線にはブリッジコンデンサ(100pF程度)を追加する
TVSダイオード配置の最適化
TVSの配置って、回路図上で繋がってればいいんじゃないんですか?
それが最大の誤解だ。回路図では同じ「GND」でも、PCB上ではインダクタンスが全然違う。TVS配置の3原則を教えよう。
- 最短経路の原則:ESD電流は「放電ポイント → TVS → GND」の最短ループを通る。このループ面積を最小化する。具体的には、TVSのGNDビアはTVSパッドの直下に配置
- 分流の原則:TVSは信号線とICの間ではなく、放電ポイントとICの間に配置する。ESD電流の大部分がTVSを通ってGNDに逃げ、ICに流入する分を最小化する
- 独立GNDの原則:TVSのGNDビアは、ICのGNDビアとは別に独立して設ける。共有するとGNDバウンスがICに直接伝わる
なるほど…。シミュレーションなしで最適配置を見つけるのは難しそうですね。
その通り。特に高密度実装のスマホ基板では、TVSの位置を1mm動かすだけでICピンの電圧が数十V変わることがある。だからこそパラメトリックスタディ——TVSの位置を数パターン試して、ICピン電圧が最小になる配置を探す——が重要なんだ。
よくある失敗と対策
初心者がやりがちな失敗ってありますか?
| 失敗パターン | 原因 | 対策 |
|---|---|---|
| TVSを付けたのにESD試験NG | TVSのGNDパスのインダクタンスが高い | GNDビアをTVS直下に最低3本配置 |
| ICは壊れないが誤動作する | GNDバウンスによるロジックレベル変動 | GNDプレーンのスリット除去、デカップリングC追加 |
| シミュレーション結果と実測が合わない | 筐体のGNDリターンパスが未モデル化 | 筐体・シャーシGNDを含めた3Dモデル構築 |
| 計算が発散する | CFL条件違反、またはポートインピーダンス不整合 | メッシュ細分化、ポート定義の見直し |
| ESD電流が実測の3倍になる | IEC波形パラメータの入力ミス(kV→V換算忘れ) | 入力単位の一覧表を作成し必ずクロスチェック |
「ESD試験でNGが出たら基板を見ろ」——現場の格言
ESD試験でICが誤動作するとき、原因の約7割はPCB上のグラウンドプレーン設計にあると言われている。放電電流が帰還経路を探して「迷走」し、コントローラのリセットピンに乗る事例は量産ライン調査でもよく見られる。ある車載ECUメーカーでは、GNDプレーンに幅0.3mmのスリットが1本入っていただけで、IEC 4kVで全数NGという事態に陥った。シミュレーションで電流密度マップを確認すると、スリットの端で電流密度が周囲の50倍に集中していることが判明。スリットを0402サイズの100pFチップコンデンサ1個で橋渡ししただけで、8kVまでクリアした。
初心者が陥りやすい落とし穴
「TVSのデータシートに"ESD保護 ±15kV"って書いてあるから大丈夫」——これが最も危険な思い込みだ。データシートのESD耐性はTVS単体の値であり、PCBに実装した状態での保護レベルではない。実装後のクランプ電圧は、TVS単体のクランプ電圧 + 配線インダクタンスによる電圧上昇 + GNDバウンス分の合算になる。データシートの値を鵜呑みにして基板を量産してしまい、IEC試験で全数NGという悲劇は珍しくない。
GNDリターンパスの考え方
ESD電流は「行き」だけでなく「帰り」もある。ESDガンから放電された電流は、被試験機のGNDを通ってESDガンのGNDケーブルに戻る。このリターンパスのインピーダンスが高いと、被試験機のGND全体が浮き上がる。シミュレーションでは、被試験機をグラウンドプレーン(GRP: Ground Reference Plane)の上に配置し、筐体のGNDケーブル接続点まで含めた完全なモデルを構築することが、実測との一致度を上げるコツだ。
ソフトウェア比較
ESD対応ツール比較
ESDシミュレーションに使えるソフトって、どんなものがありますか?
ESD解析に対応する主要ツールを比較してみよう。
| ツール名 | 主な手法 | ESD特有の機能 | SPICE連成 | 強み |
|---|---|---|---|---|
| CST Studio Suite (Dassault SIMULIA) | FDTD / FIT | IEC波形テンプレート、ESD Current Path表示 | ○ | ESD専用テンプレートが最も充実。時間領域解析の定番 |
| Ansys HFSS (Ansys Inc.) | FEM(時間/周波数) | Transient解析対応、SIwave連携 | ○(Circuit) | PCBインピーダンス解析との統合。SI/PI解析からの拡張が容易 |
| Ansys SIwave (Ansys Inc.) | MoM / FEM | ESD電流密度マップ、PDN解析 | ○ | PCBレイアウトデータ直接読み込み。GNDバウンス可視化に強い |
| COMSOL Multiphysics (COMSOL AB) | FEM | カスタム波形定義 | △(外部連携) | マルチフィジクス(熱連成など)。研究用途に柔軟 |
| Keysight EMPro (Keysight) | FDTD / FEM | ESD Current Injectionテンプレート | ○(ADS連携) | ADS/高周波回路設計との統合ワークフロー |
選定の指針
結局、初めてESD解析をやるならどれがおすすめですか?
用途別に分けるとこうなる:
- PCB + 筐体の3D解析(IEC準拠試験のバーチャル再現)→ CST Studio Suite 一択。ESD専用テンプレートがあり、学習コストが最も低い
- PCBレイアウトのGNDバウンス評価(Gerberデータから直接解析)→ Ansys SIwave。レイアウト変更のフィードバックが速い
- ICの入力保護回路の詳細設計(デバイスレベル、CDM対策)→ Cadence Spectre + Ansys RedHawk。SPICE精度のトランジスタレベル解析が必要
- 研究用途(新しい保護構造の提案、論文用)→ COMSOL。自由度が高く、カスタム物理の追加が容易
予算が限られている場合はどうすればいいですか?
オープンソースのopenEMS(FDTD法)でもESD波形の注入と時間領域解析は可能だ。ただし、GUIやESD専用テンプレートがないので、波形定義やポスト処理をスクリプトで書く必要がある。学習用途やプロトタイプ検証には十分使える。また、商用ツールの多くは無料の学生版・評価版を提供しているので、まずはそこから始めるのも手だよ。
ESD解析ツール選定——「電磁場屋」vs「回路屋」の流儀
ESDシミュレーションの世界では「電磁場屋」と「回路屋」のアプローチが合流する珍しい分野だ。電磁場屋はCST Studio SuiteでPCB全体を3Dフルウェーブ解析し、電流密度マップから問題箇所を特定する。回路屋はSPICEモデルでTVSの非線形特性を精密にモデル化し、ICピンの電圧波形から合否を判定する。2020年代に入り、両者の連成(co-simulation)が成熟し、「3D電磁場でPCBを解き、ポートにSPICEモデルを接続する」ワークフローが標準になりつつある。どちらか一方だけでは不十分な時代になった。
先端技術
チップ-パッケージ-ボード一体ESD解析
ESD解析の最先端ってどんな感じなんですか?
最大のトレンドはチップ-パッケージ-ボード(CPB)一体解析だ。従来はICチップ(CDMレベル)とPCB(IECレベル)を別々に解析していたが、最新のアプローチではICのESD保護回路のSPICEモデル + BGAパッケージの3Dモデル + PCBの3Dモデルを統合して解析する。
パッケージの中身まで解析するんですか! めちゃくちゃ計算量が多そうですね。
そう。だからCPB一体解析ではマルチスケール手法が使われる。チップレベル(μmオーダー)はSPICEで、パッケージレベル(mmオーダー)はFEMのSパラメータモデルで、ボードレベル(cmオーダー)はFDTDで、それぞれ最適な手法で解いて連結する。Ansysの「Chip-Package-System」ワークフローやSynopsysのRaptor QCMが代表的なプラットフォームだ。
機械学習によるESD耐性予測
機械学習を使ったESD解析って出てきてますか?
研究レベルでは活発だ。主に2つのアプローチがある。
- サロゲートモデル:TVS配置パラメータ(位置、GNDビア数、配線幅など)を入力、ICピン電圧を出力とするニューラルネットワークを3Dシミュレーション結果で学習させる。パラメトリックスタディが数千倍高速化される
- ESD故障予測:過去のESD試験データ(合否、故障モード、PCBレイアウト特徴量)から、新設計のESD耐性をレイアウト段階で予測する。IPC(国際電子回路協会)が2024年にガイドラインを発行した
将来的にはシミュレーションなしでESD耐性が予測できるようになるんですか?
完全な置き換えは難しいだろうけど、設計初期段階のスクリーニング(明らかにNGな設計を早期に弾く)には有効だ。「レイアウトDRC(Design Rule Check)にESDルールを組み込む」方向で、PCB設計ツールへの統合が進んでいる。KiCadやAltium Designerのプラグインとして使える日も遠くないかもしれない。
トラブルシューティング
ESD解析のデバッグ手順
シミュレーション結果が実測と合わないとき、どこから手を付ければいいですか?
ESD解析のデバッグには定番の手順がある。順番を守ることが大事だ。
- 入力波形の確認:注入しているESD電流波形がIEC規格の4パラメータ($I_1$, $I_2$, $I_3$, $t_r$)を満たしているか。単位ミス(kV→V)は最も多いバグ
- GNDリターンの確認:ESD電流の帰還経路がモデルに含まれているか。「グラウンドリフト」が起きていないか(モニターポイントを追加)
- TVSモデルの検証:TVSのSPICEモデルが正しくブレークダウンしているか。I-V特性をプロットしてデータシートと比較
- メッシュ収束性:ESD電流集中部のメッシュを2倍に細分化して、ICピン電圧が10%以内で変化しないことを確認
- 筐体の簡略化度合い:筐体のGND接続を省略していないか。金属筐体がある場合、その接地インピーダンスが結果に大きく影響する
よくあるエラーと対策
ソルバーのエラーメッセージで困ったことがあるんですが…
| エラー / 症状 | 原因 | 対策 |
|---|---|---|
| CST: "Time domain solver diverged" | CFL条件違反、PEC面の微小ギャップ | メッシュの最小セルを確認。PEC間のギャップが1セル以上あるか確認 |
| CST: "Port impedance mismatch warning" | 離散ポートのインピーダンスが実回路と不一致 | ポートを50Ωから実際のTVSのオン抵抗(0.5Ω程度)に変更 |
| HFSS: "Mesh too coarse for transient" | 最高周波数に対してメッシュが粗い | 最大周波数を $0.35/t_r$ に設定し、適応メッシュを再実行 |
| ICピン電圧が異常に高い(数百V) | TVSモデルがブレークダウンしていない | TVSのSPICEモデルの極性を確認。双方向TVSか片方向か |
| ESD電流の積分値がゼロにならない | GNDリターンパスが開放 | ESD電流源のリターン端子がGNDに接続されているか確認 |
| シミュレーション時間が異常に長い | パルス減衰が不十分でソルバーが停止しない | エネルギー減衰基準を-30dBに設定。またはシミュレーション時間を手動で $5\tau_1$ に制限 |
いやぁ、ESD解析って本当に奥が深いですね。でも先生の説明のおかげで、何をモデル化して何に注意すべきか、かなり見えてきました!
ESD対策は「回路設計」「レイアウト設計」「3D電磁場解析」の3つが交差する分野だ。どれか1つでも欠けると最適解にたどり着けない。まずは簡単なモデル(USB-Cコネクタ + TVS + IC)で一連のフローを回してみるのが一番の勉強になるよ。実際に手を動かして、電流密度マップを見たときの「おお!」という感動が、理解を加速させてくれる。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——PCB全体ではなく、TVS1個 + 信号線1本 + GNDプレーンだけの「最小ESDモデル」で現象を再現する
- 1つだけ変えて再実行——メッシュ密度だけ変える、TVSモデルだけ変える、というように対照実験の原則を守る
- 物理に立ち返る——ESD電流の総量が $Q = CV$ と一致しているか、エネルギー保存が成り立っているかを確認。基本が合っていなければ入力データの根本的なミスを疑う
なった
詳しく
報告