電磁気学の基礎 — マクスウェル方程式からCAE電磁解析まで完全解説
電磁気学は電気と磁気の相互作用を記述する物理学の体系です。現代の電気機器・電子機器・モーター・変圧器・アンテナのすべての設計において電磁解析が不可欠です。本稿では、マクスウェル方程式の物理的意味から、CAEソフトウェアで解くために必要な定式化、実務でのモデリング注意点まで、体系的に解説します。
1. 電磁気の基本量と物質定数
主要な電磁気量の定義
先生、電磁気学の式を見ると $\mathbf{E}$, $\mathbf{B}$, $\mathbf{H}$, $\mathbf{D}$ って4つもベクトルが出てきて混乱するんですが、これって全部必要なんですか？
4つが必要になるのは、物質が電磁場に応答する効果を分離して表現するためだ。ざっくり言うと：
- $\mathbf{E}$（電場）: 電荷に力を与える場。単位 V/m
- $\mathbf{D}$（電束密度）: $\mathbf{E}$ に誘電体の応答を加えたもの。 $\mathbf{D} = \varepsilon \mathbf{E}$
- $\mathbf{B}$（磁束密度）: 磁気力の根本。単位 T（テスラ）。導線に力を与える
- $\mathbf{H}$（磁場）: $\mathbf{B}$ から磁性体の応答を除いたもの。 $\mathbf{B} = \mu \mathbf{H}$
真空中では $\varepsilon = \varepsilon_0$、$\mu = \mu_0$ なので4つの区別があまり意味ない。でも変圧器の鉄心みたいな非線形磁性体では $\mathbf{B}$-$\mathbf{H}$ 曲線が本質的な役割を果たすから、区別が重要になるんだ。
| 記号 | 名称 | 単位 | 関係式（等方性線形物質） | 物理的意味 |
|---|---|---|---|---|
| $\mathbf{E}$ | 電場（電界） | V/m | — | 単位正電荷に働く力の方向・大きさ |
| $\mathbf{D}$ | 電束密度（電気変位） | C/m² | $\mathbf{D} = \varepsilon\mathbf{E} = \varepsilon_r\varepsilon_0\mathbf{E}$ | 自由電荷から発生するフラックス |
| $\mathbf{B}$ | 磁束密度（磁場） | T（= Wb/m²） | $\mathbf{B} = \mu\mathbf{H} = \mu_r\mu_0\mathbf{H}$ | 荷電粒子に働くローレンツ力の源 |
| $\mathbf{H}$ | 磁場（磁界強度） | A/m | — | 自由電流から生成される磁場 |
| $\mathbf{J}$ | 電流密度 | A/m² | $\mathbf{J} = \sigma\mathbf{E}$（オームの法則微分形） | 単位面積を流れる電流 |
| $\rho_e$ | 自由電荷密度 | C/m³ | — | 体積内の自由電荷量 |
物質定数と材料分類
| 物質定数 | 記号 | 真空値 | 材料例（代表値） |
|---|---|---|---|
| 誘電率 | $\varepsilon_0 = 8.854 \times 10^{-12}$ F/m | 比誘電率 $\varepsilon_r = 1$ | 空気≈1, ガラス≈5〜10, シリコン≈11.7, BaTiO₃≈1200 |
| 透磁率 | $\mu_0 = 4\pi \times 10^{-7}$ H/m | 比透磁率 $\mu_r = 1$ | 空気≈1, 珪素鋼板 $\mu_r \approx 1000\sim10000$, フェライト≈100〜5000 |
| 電気伝導率 | $\sigma$ | $\sigma = 0$ | 銅：$5.8 \times 10^7$ S/m, アルミ：$3.5 \times 10^7$ S/m, 珪素鋼板：$\sim2\times10^6$ S/m, シリコン（半導体）：$\sim10^{-3}\sim10^3$ S/m |
珪素鋼板の比透磁率が「1000〜10000」って書いてありますけど、この「比透磁率」が大きいと何がうれしいんですか？
比透磁率 $\mu_r$ が大きいほど「同じ電流で強い磁束を作れる」ということだ。変圧器やモーターの鉄心に珪素鋼板を使うのは、磁束を効率よく集めて閉じ込めるため。空気（$\mu_r = 1$）と比べて珪素鋼板は1万倍以上の磁束を通せる。逆に言うと、エアギャップ（空気層）があると磁気回路の「抵抗」が急増して効率が落ちる。だからモーターのエアギャップを小さくするのが設計の要だ。ただし $\mu_r$ は磁束密度 $B$ が高くなると飽和して急降下する（磁気飽和）のが非線形CAEで厄介な点だよ。
2. マクスウェル方程式
4つの方程式の全体像
電磁気学の全体は、4つのマクスウェル方程式で記述されます。積分形と微分形の両方を示します。
| 名称 | 積分形 | 微分形 | 物理的意味 |
|---|---|---|---|
| ガウスの法則（電場） | $\oint_S \mathbf{D}\cdot d\mathbf{S} = \int_V \rho_e \, dV$ | $\nabla\cdot\mathbf{D} = \rho_e$ | 電荷が電場の「湧き出し源」 |
| ガウスの法則（磁場） | $\oint_S \mathbf{B}\cdot d\mathbf{S} = 0$ | $\nabla\cdot\mathbf{B} = 0$ | 磁気単極子は存在しない |
| ファラデーの法則 | $\oint_C \mathbf{E}\cdot d\mathbf{l} = -\frac{d}{dt}\int_S \mathbf{B}\cdot d\mathbf{S}$ | $\nabla\times\mathbf{E} = -\frac{\partial\mathbf{B}}{\partial t}$ | 磁束変化が起電力を生む（電磁誘導） |
| アンペール-マクスウェルの法則 | $\oint_C \mathbf{H}\cdot d\mathbf{l} = \int_S \left(\mathbf{J} + \frac{\partial\mathbf{D}}{\partial t}\right)\cdot d\mathbf{S}$ | $\nabla\times\mathbf{H} = \mathbf{J} + \frac{\partial\mathbf{D}}{\partial t}$ | 電流と電場変化が磁場を生む |
アンペールの法則に $\partial\mathbf{D}/\partial t$ っていう項がありますけど、これって何ですか？マクスウェルが付け加えたやつですよね？
そう、これが変位電流と呼ばれるマクスウェルの最大の貢献だ。コンデンサを充電するとき、プレート間には実際の電流は流れていないのに「あたかも電流が流れているかのように」磁場が生じる。マクスウェルはこれを「電場の変化率が電流と同等の効果を持つ」と考えて式に加えた。
この変位電流項がなければ電磁波は存在できない。電場の変化が磁場を作り（変位電流）、その磁場の変化が電場を作り（ファラデー）、それが空間を伝播する——これが電磁波だ。光も電波も全部この原理だよ。低周波の静電・静磁解析では $\partial\mathbf{D}/\partial t \approx 0$ と置けるが、高周波（GHz帯など）では絶対に無視できない。
マクスウェル方程式から導かれる波動方程式
真空中（$\rho_e = 0$, $\mathbf{J} = 0$）のマクスウェル方程式から $\mathbf{B}$ を消去すると、電場の波動方程式が得られます：
これは波速 $c = 1/\sqrt{\mu_0\varepsilon_0}$ の波動方程式です。代入すると：
これは光速に一致します。マクスウェルは「光は電磁波である」ことを理論的に示したのです。
3. 静電場解析
ポアソン方程式とラプラス方程式
静電場（$\partial/\partial t = 0$）では $\nabla\times\mathbf{E} = 0$ より、電場をスカラーポテンシャル $\varphi$ で表せます：
ガウスの法則 $\nabla\cdot(\varepsilon\nabla\varphi) = -\rho_e$ に代入するとポアソン方程式：
自由電荷がない領域（$\rho_e = 0$）ではラプラス方程式：
静電場解析って具体的にどんなものを解くんですか？なんとなく「コンデンサの設計」以外は思いつかなくて。
実務でよく出てくるのは：
- 高圧機器の絶縁設計：変電所の碍子（がいし）・高圧ケーブルの接続部・変圧器内部で、電界集中が起きると絶縁破壊につながる。どこで電場が強くなるかを計算する
- 静電気対策：半導体工場のクリーンルームや火薬の製造現場で、静電気放電（ESD）を防ぐために帯電量と電位分布を解析する
- MEMS・センサ設計：静電容量型の圧力センサや加速度センサの感度計算
- ECMフィルム（エレクトレット）：マイクやエアフィルターの帯電状態の最適化
高圧機器の電界集中は、コーナーや尖端部で $E \propto 1/r$ で急激に増大するから、CAEで形状最適化して電界を均一化するのが重要な設計タスクだよ。
境界条件の設定
| 境界条件 | 数式 | 物理的意味 | 例 |
|---|---|---|---|
| ディリクレ（Dirichlet）条件 | $\varphi = \varphi_0$ （既知電位） | 電位固定 | 電極面に印加電圧を指定 |
| ノイマン（Neumann）条件 | $\partial\varphi/\partial n = 0$（自然BC） | 法線方向電場ゼロ（対称面） | モデルの対称面 |
| 界面条件（2媒質境界） | $D_{1n} = D_{2n}$, $E_{1t} = E_{2t}$ | 法線電束密度・接線電場の連続 | 誘電体境界面 |
| 導体表面 | $\mathbf{E}_t = 0$（内部はゼロ） | 導体は等電位面 | 金属電極 |
静電容量の計算
2つの導体間の静電容量 $C$ は、ポテンシャル解析から求めた電場エネルギーで計算できます：
平行平板コンデンサの理論値との比較でFEMの検証ができます：$C = \varepsilon A/d$（$A$：面積, $d$：間隔）。フリンジ（端部）効果がある場合はFEMによる解析が理論値より大きい値を与えます。
4. 静磁場解析
ベクトルポテンシャルと磁束密度
静磁場（$\partial/\partial t = 0$）では $\nabla\cdot\mathbf{B} = 0$ より、磁束密度をベクトルポテンシャル $\mathbf{A}$ で表します：
アンペールの法則 $\nabla\times\mathbf{H} = \mathbf{J}$ に $\mathbf{B} = \mu\mathbf{H}$ と代入し、クーロンゲージ $\nabla\cdot\mathbf{A} = 0$ を適用すると：
均一・等方性の線形磁性体では：
これはポアソン方程式（電場の $\nabla^2\varphi = -\rho/\varepsilon$ に対応）です。2D軸対称問題では $\mathbf{A} = A_z\hat{z}$ とスカラー化でき、計算が大幅に簡略化されます（コイル・モーター・変圧器の2Dモデルで多用）。
磁気回路と磁気抵抗
「磁気回路」って電気回路のアナロジーって聞きました。具体的にどう使うんですか？
電気回路のオームの法則 $V = IR$ に対応して、磁気回路には「磁気オームの法則」がある：
- 起磁力 $\mathcal{F} = NI$ [A]（電圧に対応）← コイルの巻数×電流
- 磁束 $\Phi = \int \mathbf{B}\cdot d\mathbf{S}$ [Wb]（電流に対応）
- 磁気抵抗 $\mathcal{R} = l/(\mu A)$ [A/Wb]（電気抵抗に対応）
$\mathcal{F} = \mathcal{R}\Phi$ という関係式が成り立つ。エアギャップの磁気抵抗は鉄心の数百倍になるから、「なるべく鉄で磁束を閉じる」のが鉄心設計の基本だ。変圧器・インダクタ・モーターの一次設計でよく使われるよ。ただし非線形（磁気飽和）や磁束漏れがある場合はFEM解析に頼る必要がある。
インダクタンスの計算
コイルのインダクタンス $L$ は磁場エネルギー $W_m$ または磁束鎖交数 $\Lambda = N\Phi$ から求めます：
5. 渦電流解析（時変磁場）
誘導起電力と渦電流の発生メカニズム
時変磁場がある場合、ファラデーの法則より導電性材料内に誘導電場が生じ、渦電流（エディカレント）が発生します：
正弦波励磁（角周波数 $\omega$）の場合、ベクトルポテンシャルの支配方程式は：
ここで $j = \sqrt{-1}$ は虚数単位（複素ポテンシャル法）です。この方程式は拡散方程式型で、時定数 $\tau = \mu\sigma L^2$ （$L$：代表長さ）を持ちます。
表皮効果と表皮深さ
「表皮効果」っていう言葉を聞いたんですが、これって交流電流が導体の表面にしか流れなくなる現象ですよね？なんで表面だけになるんですか？
そう。電流が流れると磁場が生じて（アンペール則）、その磁場が変化すると誘導電場が生じる（ファラデー則）。この誘導電場は「もとの電流を打ち消す方向」に働く。この打ち消し効果が導体内部で積み重なると、高周波では中心部の電流が抑制されて表面に集中する、という自己一貫したメカニズムだ。
周波数が高いほど表皮が薄くなる。表皮深さ $\delta$ は：$\delta = \sqrt{2/(\omega\mu\sigma)}$（周波数が4倍なら $\delta$ は1/2）。50 Hzの商用電源で銅の表皮深さは約9 mmだが、10 MHzでは約21 μmしかない。高周波コイルの設計ではこれが重要で、導線を細い撚線（リッツ線）にして表皮効果による損失を減らす工夫をするよ。
電流密度は表面から深さ $y$ の位置で指数関数的に減衰します：
振幅は $|J(y)| = J_0 e^{-y/\delta}$、深さ $\delta$ でほぼ1/e（約37%）まで減衰します。
| 材料 | 電気伝導率 σ [S/m] | 比透磁率 μ_r | 表皮深さ δ（50 Hz） | 表皮深さ δ（1 MHz） |
|---|---|---|---|---|
| 銅 | 5.8×10⁷ | 1 | 9.3 mm | 65.7 μm |
| アルミニウム | 3.5×10⁷ | 1 | 12.0 mm | 84.7 μm |
| 鉄（純鉄） | 1.0×10⁷ | 1000 | 0.71 mm | 5.0 μm |
| 珪素鋼板（3%Si） | 2.0×10⁶ | 5000 | 0.45 mm | 3.2 μm |
| ステンレス（SUS304） | 1.4×10⁶ | 1 | 60 mm | 424 μm |
渦電流損失の計算と実務例
渦電流損失（ジュール熱）の体積密度は：
変圧器鉄心の渦電流損失の近似式（Steinmetz）：
ここで $k_e$ は材料定数、$f$ は周波数、$B_m$ は最大磁束密度、$d$ は積層板の厚みです。鉄心を薄い板（0.1〜0.5 mm）に積層する（積層鋼板）のはこの $d^2$ 依存性を利用した損失低減策です。
誘導加熱（IHクッキングヒーター）も渦電流の原理ですよね？どうして非磁性の鍋では使えないんですか？
IHは20〜100 kHzの交流磁場を鍋底に作用させて渦電流を誘導し、鍋自体をジュール加熱する仕組みだ。非磁性のアルミや銅の鍋でも渦電流は流れるが、電気抵抗が低すぎて大きな電流が流れてもあまり発熱しない、という問題がある。逆に磁性鍋（鉄・ステンレスSUS430）は比透磁率が高いから表皮深さが小さく、表面近くに電流が集中して発熱効率が高い。
CAEでの誘導加熱解析は「電磁場と熱の連成問題」で、電磁場解析でジュール発熱分布を求め、それを熱解析のソースとして温度分布を計算する。Ansys Maxwell + Fluent の連成解析や JMAG などの専用ソルバーが使われるよ。
6. 電磁波と高周波解析
電磁波の基本特性
誘電体中の電磁波の波数ベクトル $\mathbf{k}$ と波長 $\lambda$：
ここで $n = \sqrt{\mu_r\varepsilon_r}$ は屈折率です。導電性媒質中では波数が複素数となり、電磁波は減衰しながら伝播します（表皮効果と同じ現象）。
平面波インピーダンスと反射
媒質の波動インピーダンス：
界面での反射係数（法線入射）：
EMC/EMI解析のポイント
EMCって「電磁両立性」ですよね。CAEでどんな解析をするんですか？
EMC解析の主な目的は2つ：①自分が出すノイズ（EMI：電磁妨害）を規制値以下に抑える、②外来ノイズに対して誤動作しない（EMS：電磁感受性）こと。電気自動車では高電圧パワーモジュールのスイッチング電流が強いEMIを出すから、シールドや筐体の設計が重要だ。
CAE解析では：波長が基板・ケーブルサイズに比べて長い「準静電磁場領域」（〜数MHz）ではFEM/BEM、波長と構造サイズが同程度の「共振領域」（数MHz〜数GHz）ではFEMまたはFDTD（時間領域差分）、波長が短い高周波（GHz以上）はSBR（射線追跡）などを使い分ける。
アンテナ解析の基礎
ダイポールアンテナの放射電力密度（ポインティングベクトル $\mathbf{S} = \mathbf{E}\times\mathbf{H}$）と指向性 $D$：
半波長ダイポールアンテナの指向性 $D = 1.64$（2.15 dBi）が理論値です。CAEで近傍電磁場を計算し、遠方場変換（NF-FF transformation）で放射パターンを求めます。
7. FEMによる電磁解析の実装
スカラーポテンシャル法 vs ベクトルポテンシャル法
| 手法 | 未知変数 | 支配方程式 | 適用場面 | 注意点 |
|---|---|---|---|---|
| スカラーポテンシャル法 | $\varphi$（スカラー） | $\nabla^2\varphi = -\rho/\varepsilon$ | 静電場（電流なし） | 電流領域（$\mathbf{J}\neq 0$）では使えない |
| ベクトルポテンシャル法 | $\mathbf{A}$（ベクトル, 3成分） | $\nabla\times(\frac{1}{\mu}\nabla\times\mathbf{A}) = \mathbf{J}$ | 静磁場・渦電流 | ゲージ条件（クーロンゲージ等）が必要 |
| $T$-$\Omega$ 法 | 電流ポテンシャル $T$ + 磁気スカラーポテンシャル $\Omega$ | — | 渦電流（低コスト） | 導体が単連結の場合のみ使用可 |
| $A$-$V$ 法 | $\mathbf{A}$（磁) + $V$（電位） | 渦電流の完全連成 | 渦電流（一般的） | 自由度が大きい |
エッジ要素（ネデレック要素）の必要性
構造解析のFEMでは通常の節点要素を使いますよね。電磁解析では「エッジ要素」っていう特殊な要素が使われるって聞いたんですが、なんで普通の要素じゃダメなんですか？
ベクトルポテンシャル $\mathbf{A}$ を節点要素で離散化すると、2つの問題が起きる。①「不正な解（spurious mode）」が生じる——電磁的に意味のない数値振動が現れる。②導電体と絶縁体の境界で $\mathbf{B}$ の法線成分連続の条件が自然に満たされない。
エッジ要素（ネデレック要素）は節点ではなく「辺（エッジ）」に自由度を置く。これにより $\nabla\times\mathbf{A}$ の接線成分連続が自動的に保証され、スプリアスモードが生じない。Ansys Maxwell、JMAG、GetDP（オープンソース）などの電磁解析ソルバーは標準でエッジ要素を使っているよ。
ゲージ条件と一意性の確保
$\mathbf{A}$ の定義に任意性がある（ゲージ自由度）ため、解を一意に決めるゲージ条件が必要です：
| ゲージ条件 | 式 | 適用場面 |
|---|---|---|
| クーロンゲージ | $\nabla\cdot\mathbf{A} = 0$ | 静磁場・低周波渦電流 |
| ローレンツゲージ | $\nabla\cdot\mathbf{A} + \mu\varepsilon\partial\varphi/\partial t = 0$ | 波動方程式（電磁波解析） |
| ツリーコツリーゲージ | 有限要素のグラフ理論的手法 | エッジ要素での数値実装 |
8. モーター・アクチュエータのCAE
マクスウェル応力テンソルによるトルク計算
電磁気的な力はマクスウェル応力テンソルの面積分として計算できます：
ロータに働くトルク $M$ は、エアギャップ周囲の閉曲面 $S$ 上での応力テンソルの積分：
実務的なモーターCAEでは、この積分をエアギャップ内の「仮想閉曲面」上で実行します。積分面の位置を変えても同じ結果になることが、計算の検証に使えます。
磁気飽和と非線形B-H曲線
モーターの鉄心を磁気飽和させると何がダメなんですか？効率が落ちる以外に問題があるんですか？
いくつかある。①トルクの非線形化：設計通りのトルクが出なくなる。電流を2倍にしても磁束密度が2倍にならないから。②鉄損の増大：高い磁束密度域では渦電流損失・ヒステリシス損失が急増して発熱する。③コギングトルクの増大：スロットの鉄心が局所的に飽和するとトルクリップルが大きくなり、振動・騒音の原因になる。
CAEでは非線形の B-H 曲線（磁化曲線）データを材料として入力し、Newton-Raphson法などで非線形方程式を反復解法で解く。B-H曲線は材料メーカーのデータシートから取得するか、試験で測定する。初期透磁率・飽和磁束密度・残留磁束密度・保磁力の4つが重要パラメータだ。
コギングトルク解析と低減技術
コギングトルク（無通電時のトルクリップル）はFEMで以下のアプローチで解析します：
- ロータをステップ回転させて各角度でのトルクを計算（回転ステップ：電気角1°程度）
- 角度依存のトルク波形をフーリエ分解して高調波成分を特定
- スキュー（ロータまたはステータのスロットを傾ける）やスロット開口幅最適化で低減
コギングトルク係数の推定式：
ここで $\Phi$ は磁束、$\mathcal{R}$ は磁気抵抗、$\theta$ は回転角です。エアギャップの磁気抵抗変化を小さくすることがコギングトルク低減の本質です。
9. 実務でのよくある失敗
境界条件の設定ミス（無限遠の扱い）
静磁場解析でモデルを作るとき、「解析領域の外は空気」なのにどこかで切らないといけないですよね。その「切り方」を間違えると結果に影響するんですか？
大いに影響する。電磁場は理論的には無限遠まで広がるので、有限のメッシュ領域で打ち切ると境界での「反射」が生じる。対策は3つ：
- 解析領域を十分大きく取る：コイル径の5〜10倍の空気領域を確保。一番シンプルだが計算コスト大
- 無限要素（Infinite Element）：外側の要素内で座標変換して無限遠を表現。Ansys Maxwell等で対応
- 完全一致層（PML: Perfectly Matched Layer）：電磁波解析での完全吸収境界。高周波解析で標準的
外部コイル・アンテナ解析では特に重要。内部磁場解析（変圧器・電動機の鉄心内部）は磁束が閉じているから外部条件の影響が小さくて楽だよ。
メッシュ密度の設計注意点
| 箇所 | 理由 | 推奨メッシュサイズ |
|---|---|---|
| エアギャップ（モーター・変圧器） | 磁気抵抗の大部分を占める。精度が全体に影響 | エアギャップ厚みの方向に最低6〜10層 |
| 鉄心コーナー・エッジ | 電界・磁界集中、飽和が局所的に起きる | コーナー近傍で10倍程度の細かいメッシュ |
| 渦電流解析の表皮深さ内部 | 表皮深さ未満の要素が必要 | 第一層サイズ ≤ $\delta/5$、深さ方向に3〜5層 |
| 高周波（波動解析） | 数値分散を防ぐ | 波長の1/10〜1/20 以下の要素サイズ |
周波数依存性の無視による失敗例
周波数依存性を無視するとどんな失敗になるんですか？具体例を教えてください。
よくある失敗例を3つ挙げると：
- インバータ駆動モーターの鉄損計算：正弦波（50 Hz）の鉄損データで計算すると、数kHz成分を含むPWM波形での実際の鉄損を大幅に過小評価する。高調波を含めた損失計算が必要
- 高周波トランス（スイッチング電源）の漏れインダクタンス：直流やDC抵抗として設計すると、数百kHzでの表皮効果による抵抗増加で損失が予測の数倍になることがある
- GHz帯PCBのコネクタ設計：低周波での静電容量は問題なくても、GHz帯では共振モードが生じてインピーダンス不整合が起きる。EMC試験で初めて発覚するケースが多い
解決策は、設計段階で想定動作周波数帯を明確にして、その帯域をカバーする周波数掃引解析（AC Sweep / Frequency Sweep）を実施すること。