キャパシタンス解析
理論と物理
概要
先生、静電容量ってFEMでどう計算するんですか？
静電容量は導体間の電荷と電圧の比だ。
FEMで電位のポアソン方程式を解き、導体表面の電荷をガウスの法則で積分して求める。
代表的な解析解
理論式と比較できる形状はありますか？
| 構造 | 容量 |
|---|---|
| 平行平板 | $C = \varepsilon A / d$ |
| 同軸円筒 | $C = 2\pi\varepsilon L / \ln(b/a)$ |
| 同心球殻 | $C = 4\pi\varepsilon ab / (b-a)$ |
| 孤立球導体 | $C = 4\pi\varepsilon_0 a$ |
これらでFEM結果の妥当性を確認してから実問題に進むのが鉄則だ。
容量行列
多導体系ではどう扱うんですか？
$n$ 個の導体間の容量関係は容量行列 $[C]$ で記述される。
対角成分 $C_{ii}$ は自己容量、非対角成分 $C_{ij}$ は相互容量で負の値を取る。PCBの配線間クロストーク評価にはこの容量行列が不可欠だ。Ansys Q3Dでは自動抽出が可能だよ。
ライデン瓶——人類初のコンデンサは1745年のオランダで生まれた
静電容量（キャパシタンス）の概念を体現した人類初のコンデンサが「ライデン瓶」です。ガラス瓶の内外を金属箔で覆い、内側に静電気を蓄えます。1746年にフランクリンがこれを使い、雷が電気であることを証明する凧実験を行いました。現代のキャパシタンス理論 $C = Q/V$ の本質——「電位差に対してどれだけの電荷を蓄えられるか」——は、この瓶の観察から着想されたのです。
各項の物理的意味
- 電場項 $\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仮定（断面解析）：電流方向が一様で、端部効果を無視できる場合に有効
- 等方性仮定：異方性材料（珪素鋼板の圧延方向等）では方向別の特性定義が必要
- 適用外ケース：プラズマ（電離気体）、超伝導体、非線形光学材料では追加の構成則が必要
数値解法と実装
数値解法の詳細
容量行列をFEMでどう計算するんですか？
$n$ 個の導体がある場合、導体 $j$ に1Vを印加し他を0Vにしてポアソン方程式を解く。各導体表面の電荷 $Q_i$ がそのまま $C_{ij}$ になる。これを全導体について繰り返す。
$n$ 回も解析するのは大変ですね。
対称性 $C_{ij} = C_{ji}$ を利用すれば回数を減らせる。Ansys Q3Dではこの操作が完全に自動化されている。
エネルギー法
電荷積分以外の方法はありますか？
静電エネルギーから容量を求めるエネルギー法がある。
COMSOLの「Volume Integration」で計算可能だ。電荷法とエネルギー法の両方で算出し、差が1%以内なら信頼できる。差が大きければメッシュが粗い証拠だ。
フリンジ電界が大きい構造だとメッシュ不足になりやすいですよね。
その通り。配線幅とギャップが同程度の場合、フリンジ容量が全体の30〜50%を占めることもある。$C = \varepsilon A/d$ の平板近似だけでは不十分だ。
「平行平板コンデンサ」が教科書から消えない理由
キャパシタンス解析の基本として必ず登場する平行平板コンデンサの公式 $C = \varepsilon_0 \varepsilon_r A / d$。「こんな単純な形、現実にはない」と思いきや、半導体のゲート酸化膜、PCB基板の電源プレーン、フレキシブル基板の積層間——いずれも本質的に平行平板です。実際の設計では端効果（フリンジング効果）が10〜30%の誤差を生じるため、CAEソルバーによる補正が必須。教科書の公式はスタート地点にすぎません。
辺要素（Nedelec要素）
電磁場解析に特化した要素。接線成分の連続性を自動的に保証し、スプリアスモードを排除。3D高周波解析の標準。
節点要素
スカラーポテンシャル定式化に使用。静磁場のスカラーポテンシャル法や静電場解析で有効。
FEM vs BEM（境界要素法）
FEM: 非線形材料・非均質媒質に対応。BEM: 無限領域（開領域問題）を自然に扱える。ハイブリッドFEM-BEMも有効。
非線形収束（磁気飽和）
B-Hカーブの非線形性をニュートン・ラフソン法で処理。残差基準: $||R||/||R_0|| < 10^{-4}$が一般的。
周波数領域解析
時間高調波仮定により定常問題に帰着。複素数演算が必要だが、広帯域特性は時間領域解析で取得。
時間領域の時間刻み
最高周波数成分の1/20以下の時間刻みが必要。暗黙的時間積分ではより大きな刻みも可能だが精度に注意。
周波数領域と時間領域の使い分け
周波数領域解析は「ラジオの特定の周波数に合わせる」ようなもの——1つの周波数での応答を効率的に計算できる。時間領域解析は「全チャンネルを同時に録画する」ようなもの——あらゆる周波数成分を含む過渡現象を再現できるが計算コストが高い。
実践ガイド
実践ガイド
静電容量解析の実務での使い方を教えてください。
代表的な適用例を挙げよう。
PCB配線の寄生容量
- Ansys Q3DでECADデータ（ODB++/Gerber）をインポート
- 信号線間の相互容量を抽出しSPICEモデルとして出力
- SI解析の等価回路に反映
静電容量型センサ
- COMSOLでセンサ電極間の容量変化を計算
- 被測定物の位置・厚さに対する容量の感度解析
メッシュ設計
容量解析特有のメッシュの注意点はありますか？
| 部位 | メッシュ要件 | 理由 |
|---|---|---|
| 導体間ギャップ | ギャップ幅の1/5以下 | 容量精度に直結 |
| 導体エッジ | 微小要素を配置 | フリンジ容量の捕捉 |
| 誘電体界面 | メッシュ整合 | 電束密度の連続性 |
| 遠方領域 | 粗いメッシュ可 | 容量への寄与小 |
フリンジ容量って結構大きいんですね。
配線のアスペクト比が1に近い場合、フリンジ容量が全体の半分を占める。無限平板近似は危険だ。必ず3D FEMかBEMで正確に計算すること。
タッチパネルのキャパシタンス——指一本で変わる0.1pFを検出する技術
スマートフォンの投影型静電容量式タッチパネルは、格子状の電極間の静電容量変化（約0.1〜1pF）を検出して指の位置を割り出します。これほど微小な変化を毎秒120回以上スキャンするため、電極パターンの設計にはキャパシタンス解析が欠かせません。電極間のギャップ・幅・誘電体層の厚みを最適化しないと「指が太くないと反応しない」「端っこだけ感度が悪い」というQA問題が現場で頻発します。
解析フローのたとえ
モータの電磁界解析は「ギターの調律」に近い感覚です。弦の太さ（コイル巻数）とブリッジの位置（磁石配置）を調整して、最も美しい音色（効率の良いトルク特性）を引き出す。1つのパラメータを変えると全体のバランスが変わる——だからパラメトリックスタディが重要なんです。
初心者が陥りやすい落とし穴
「空気領域？ なんで空気をメッシュで切るの？」——初めて電磁界解析に触れた人がほぼ全員抱く疑問です。答えは「磁力線は鉄心の外にも広がるから」。解析領域を鉄心ぎりぎりにすると、行き場を失った磁束が壁に「ぶつかって」反射し、実際にはありえない磁束集中が起きます。部屋が狭すぎてボールが壁に跳ね返りまくる状態を想像してみてください。
境界条件の考え方
遠方の境界条件って地味ですが超重要です。「ここから先は無限に広がる空間」ということを数値的に表現する必要がある。設定を間違えると、まるで「見えない壁」があるかのように磁束が跳ね返されてしまいます。
ソフトウェア比較
商用ツール比較
容量解析に使えるツールを比較してください。
用途別に整理しよう。
| ツール | 手法 | 得意分野 |
|---|---|---|
| Ansys Q3D | BEM/FEM | PCB/パッケージの寄生容量抽出 |
| Ansys Maxwell | FEM | 汎用3D静電場、大型機器 |
| COMSOL AC/DC | FEM | カスタム物理、センサ設計 |
| Ansys SIwave | 2.5D FEM | PCB基板全体の高速解析 |
| FastCap | BEM | オープンソース、研究用 |
Q3Dの具体的なワークフローは？
1. 3Dモデルをインポート
2. 信号端子とグラウンドを定義
3. 自動メッシュ生成・適応細分化
4. R, L, C, G行列を周波数ごとに抽出
5. SPICEネットリストとして出力
パワーモジュールのバスバー設計では寄生インダクタンスも同時に抽出でき、スイッチングサージ予測に直結する。
COMSOLとの使い分けは？
COMSOLはセンサなど構造-静電連成が必要な場合に強い。Q3Dは多導体系の自動容量抽出に特化しており、ICパッケージやPCBでは圧倒的に効率が良い。
EDAとCAEのキャパシタンス解析——なぜ「二重管理」が起きるのか
半導体・基板設計の現場では、EDA（電子設計自動化）ツールの寄生容量抽出とCAEツールの電磁界解析という二つの流れが並行して走ることがよくあります。EDAは配線ルール主体で高速・自動、CAEは3D形状を精密に解くが時間がかかる。精度が必要なクリティカルな配線だけCAEで検証し、残りはEDAに任せるという「ハイブリッド運用」が実務では一般的です。ツール選定の議論の前に、まずワークフロー設計から始めるのが得策です。
選定で最も重要な3つの問い
- 「何を解くか」：キャパシタンス解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端トピック
静電容量解析の最新研究を教えてください。
いくつかの重要テーマがある。
微細化半導体の寄生容量
3nm/2nmプロセスではゲート間やメタル配線間の寄生容量が素子性能を支配する。FinFETやGAA構造では3Dの電界分布が複雑で、Synopsys RaphaelやCalibre xACTで大規模3D容量抽出が行われている。
可変容量MEMS
MEMSの静電容量スイッチの解析は？
電極間の容量が機械変位で変化する非線形連成問題だ。プルイン電圧は
COMSOLの電気-構造連成で静電力→変形→容量更新を反復して予測する。
量子容量
グラフェンなど2D材料ではキャリア状態密度が有限なため量子容量が発生する。
幾何学的容量と直列に加わり全容量を制限する。TCADやab initio計算との連成で研究が進んでいる。
次世代半導体チップの静電容量——フェムトファラッドの世界
最新の3nmプロセスのトランジスタでは、ゲート静電容量がフェムトファラッド（fF、10⁻¹⁵ F）のオーダーです。これほど小さな容量でも、数十億個のトランジスタが集積されると、チップ全体の寄生容量がパフォーマンスを大きく左右します。設計ツールが自動でキャパシタンス解析を走らせ、配線間の容量結合（クロストーク）を事前に把握するのは今や当然の工程。「作ってみたら遅かった」では許されない時代です。
トラブルシューティング
トラブルシューティング
容量解析のよくあるトラブルを教えてください。
代表的な問題を整理しよう。
1. FEM容量値が理論値と合わない
平行平板で $\varepsilon A/d$ と10%ずれます。
原因: フリンジ電界。理論式は無限平板の近似。Palmer補正式 $C = \varepsilon_0 \frac{A}{d}[1 + \frac{d}{\pi w}(1 + \ln\frac{2\pi w}{d})]$ で比較すること。
2. Q3Dと実測が合わない
寄生容量が20%もずれます。
原因: モデルに含まれていない周辺構造（筐体、基板）の影響。PCB基板のDk値は周波数依存。LCRメータの測定周波数での分布定数効果。
対策: 周辺構造を含めた広域モデルで再計算し、誘電率の±10%ばらつきで感度解析を行う。
3. 容量行列の対角要素が負
容量行列がおかしいです。
Maxwell形式とKirchhoff形式の違いだ。Maxwell形式では非対角成分が負になるのが正常。Q3Dの出力形式を確認すること。
4. 適応メッシュが収束しない
Maxwellの適応パスが何十回も回ります。
対策: エネルギー誤差を1%から2%に緩和。初期メッシュを導体間ギャップに手動集中。微小フィレットを簡略化。
キャパシタンス解析で「値がおかしい」——まず疑うべき3箇所
キャパシタンス解析の結果が設計値と大きく外れるとき、原因の大半は①メッシュが粗すぎて電極端部の電場集中を捉えられていない、②誘電体の比誘電率をDC値で入力したが実際には周波数依存性がある、③フリンジング（端効果）を考慮していない、のどれかです。特に①は「メッシュを2倍に細かくしたら値が10%変わった」という経験をすると身に染みます。収束確認は最低でも2段階のメッシュリファインで行いましょう。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——キャパシタンス解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告