クーロンの法則
理論と物理
クーロンの法則
先生、クーロンの法則って電磁気学の出発点ですよね？
2つの点電荷間に働く力を記述する最も基本的な法則だ。
$\varepsilon_0 = 8.854 \times 10^{-12}$ F/m（真空の誘電率）。同符号なら斥力、異符号なら引力。
距離の2乗に反比例…重力と同じ形ですね。
そう。ただし電気力は重力の$10^{36}$倍の強さがある。電子2個の間の電気力と重力を比較すると、電気力が圧倒的に強い。日常で電気力を感じにくいのは、正負の電荷がほぼ完全に打ち消し合っているからだ。
電界（電場）
点電荷 $q$ が作る電界：
$$ \mathbf{E} = \frac{1}{4\pi\varepsilon_0}\frac{q}{r^2}\hat{\mathbf{r}} $$
点電荷 $q$ が作る電界：
電界は「単位電荷あたりの力」。電界が分かれば、任意の電荷 $Q$ に働く力は $\mathbf{F} = Q\mathbf{E}$ で計算できる。
重ね合わせの原理
複数の電荷がある場合、各電荷が作る電界のベクトル和が全体の電界。線形性が成り立つ。
連続分布の場合は積分に置き換わる。これがFEMで静電場を解く基礎。
まとめ
クーロンのトーションバランス——1785年の驚異的な精度
クーロンが逆二乗則を実証したのは1785年。使った道具は「ねじりばかり（トーションバランス）」という、細い金属線のねじれ角で力を測る装置でした。帯電した球の間の力をその微小なねじれ角から読み取るというアナログ手法で、現代の計測器と比べても驚くほど精度の高い結果を出しています。当時はまだオームの法則もアンペールの法則も存在せず、「電気力の強さ」すら定量的に知られていなかった時代。クーロンが地道な実験から距離の2乗に反比例する関係を導いたことで、後の電磁気学の理論体系が成立しました。CAEで当たり前のように使う $F = kq_1q_2/r^2$ の裏には、こんな職人的な実験があったんです。
各項の物理的意味
- 電場項 $\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
クーロンの法則をFEMでどう解くんですか？
直接クーロンの法則を計算するのではなく、ポアソン方程式に帰着させる。
$\phi$: 電位、$\rho_v$: 電荷密度。FEMで$\phi$を求め、$\mathbf{E} = -\nabla\phi$で電界を計算。
ポアソン方程式を解くほうが効率的なんですね。
直接的なクーロン力の計算は$O(N^2)$だが、FEMなら疎行列を解くだけで済む。BEM（境界要素法）も開放空間の問題に有効。
まとめ
点電荷モデルで解ける問題、解けない問題
クーロンの法則の数値実装で一番ハマるのが「点電荷の近傍」。距離 $r \to 0$ にすると力が無限大に発散するので、FEMやBEMでは電荷を有限サイズの帯電体として扱うか、特異点処理を入れる必要があります。実務でよくあるのが「針状電極の先端に電界集中が起きてシミュレーションが発散する」問題で、現場では先端を微小な球でモデル化して特異性を逃がす、という技が使われます。理論式がシンプルでも、数値実装では「点電荷」という理想化がそのまま使えないのが面白いところですね。
辺要素（Nedelec要素）
電磁場解析に特化した要素。接線成分の連続性を自動的に保証し、スプリアスモードを排除。3D高周波解析の標準。
節点要素
スカラーポテンシャル定式化に使用。静磁場のスカラーポテンシャル法や静電場解析で有効。
FEM vs BEM（境界要素法）
FEM: 非線形材料・非均質媒質に対応。BEM: 無限領域（開領域問題）を自然に扱える。ハイブリッドFEM-BEMも有効。
非線形収束（磁気飽和）
B-Hカーブの非線形性をニュートン・ラフソン法で処理。残差基準: $||R||/||R_0|| < 10^{-4}$が一般的。
周波数領域解析
時間高調波仮定により定常問題に帰着。複素数演算が必要だが、広帯域特性は時間領域解析で取得。
時間領域の時間刻み
最高周波数成分の1/20以下の時間刻みが必要。暗黙的時間積分ではより大きな刻みも可能だが精度に注意。
周波数領域と時間領域の使い分け
周波数領域解析は「ラジオの特定の周波数に合わせる」ようなもの——1つの周波数での応答を効率的に計算できる。時間領域解析は「全チャンネルを同時に録画する」ようなもの——あらゆる周波数成分を含む過渡現象を再現できるが計算コストが高い。
実践ガイド
静電場解析の実務
高電圧機器の絶縁設計、半導体のゲート電界、MEMS静電アクチュエータが典型的な適用例。
実務チェックリスト
静電塗装でクーロン力を使い倒す
自動車のボディ塗装に使われる「静電塗装」は、クーロン力を実用化した代表例です。塗料の粒子をマイナスに帯電させ、アース（接地）されたボディ側に向けて吹き付けると、静電引力で塗料が表面に均一に吸着します。塗り残しが出やすい入り組んだ部分にも塗料が回り込むのはクーロン力のおかげで、実際に塗着効率は通常塗装の2倍以上になることもあります。CAE的には電界分布をシミュレーションして「塗料が届きにくいポケット形状」を事前に特定する、という使い方が現場で増えています。
解析フローのたとえ
モータの電磁界解析は「ギターの調律」に近い感覚です。弦の太さ（コイル巻数）とブリッジの位置（磁石配置）を調整して、最も美しい音色（効率の良いトルク特性）を引き出す。1つのパラメータを変えると全体のバランスが変わる——だからパラメトリックスタディが重要なんです。
初心者が陥りやすい落とし穴
「空気領域？ なんで空気をメッシュで切るの？」——初めて電磁界解析に触れた人がほぼ全員抱く疑問です。答えは「磁力線は鉄心の外にも広がるから」。解析領域を鉄心ぎりぎりにすると、行き場を失った磁束が壁に「ぶつかって」反射し、実際にはありえない磁束集中が起きます。部屋が狭すぎてボールが壁に跳ね返りまくる状態を想像してみてください。
境界条件の考え方
遠方の境界条件って地味ですが超重要です。「ここから先は無限に広がる空間」ということを数値的に表現する必要がある。設定を間違えると、まるで「見えない壁」があるかのように磁束が跳ね返されてしまいます。
ソフトウェア比較
ツール
| ツール | 特徴 |
|---|---|
| COMSOL AC/DC | 静電場の標準。マルチフィジックス連成 |
| Ansys Maxwell | 3D静電場。電界分布の可視化 |
| FEMM | 2D静電場。無償。教育・初期検討用 |
| Elmer | オープンソースFEM。静電場対応 |
静電界ソルバーの「無限遠境界」をどう扱うか
商用ツールで静電界解析をするとき、「解析領域の外側をどうするか」という問題があります。理論的にはクーロン力は無限遠まで届くので、解析領域を有限に切り取るには何らかの近似が必要です。ABAQUSやANSYSでは「無限要素（Infinite Element）」を境界に配置して遠方場を吸収する手法が使われます。一方COMSOLは「無限変換」を内部で自動処理するオプションが便利です。ツール選定では「この無限遠処理がどのくらい精度よくできるか」も見ておくと、高精度な静電解析で後悔しません。
選定で最も重要な3つの問い
- 「何を解くか」：クーロンの法則に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端技術
クーロン定数 $k$ と真空の誘電率——実は定義値だった
クーロン定数 $k \approx 8.99 \times 10^9$ N·m²/C² は、実は真空の誘電率 $\varepsilon_0$ と $k = 1/(4\pi\varepsilon_0)$ で結びついています。面白いのは、2019年のSI単位系改定以降、$\varepsilon_0$ は測定値ではなく「定義値」になったこと。電気素量 $e$ が厳密に定義されたことで、クーロン定数も厳密な値に固定されました。つまり今日のCAEで使う静電解析の根拠となる定数は、もはや実験で測るものではなく、人類が「そう決めた」値なんです。物理定数の定義って、ちょっと哲学的だな、と思いませんか？
トラブルシューティング
トラブル
「電荷が動いたら静電界じゃない」——境界の見極め方
クーロン力の計算でハマりやすいのが、「本当に静電界として扱っていいのか」という判断です。例えばプリント基板の配線間の静電力を計算しようとしたとき、実際には交流信号が流れているので厳密には準静電界（時変電磁界）です。ただし周波数が低ければ「静電界近似」でも十分な精度が出る。現場では「動作周波数の波長 $\lambda$ に対して装置サイズが $\lambda/10$ 以下なら静電界近似でOK」という経験則がよく使われます。この判断を間違えると解析モデルの選定ミスになるので、トラブル対応の第一歩として必ず確認するポイントです。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——クーロンの法則の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告