電磁力 — 数値解法と実装
FEMでの力計算手法
FEMで電磁力を計算する方法は複数あるんですか？
主に3つの手法がある。
1. マクスウェル応力テンソル法 — エアギャップ中の積分面で力を計算。積分面の位置に結果が依存するため複数面で平均を取る
2. 仮想仕事法 — 物体を微小変位させたときのエネルギー変化から力を算出。精度が最も高い
3. Arkkio法 — エアギャップの体積積分。回転機のトルク計算の標準手法
JMAGではどの手法が使われますか？
JMAGではノーダルフォース法も実装されている。各節点に作用する電磁力を直接計算し、構造解析の荷重として受け渡せる。Ansys Maxwellでは仮想仕事法がデフォルト。
まとめ
- マクスウェル応力 — 簡便だが積分面依存
- 仮想仕事法 — 高精度だが計算コスト大
- Arkkio法 — 回転機トルクの標準
ファラデー——「数学が苦手だった」天才
電磁誘導の法則を発見したマイケル・ファラデーは、正規の教育を受けておらず、高等数学が使えませんでした。彼は「力線」という直感的なイメージで電磁気現象を理解し、実験で次々と発見をしました。後にマクスウェルがファラデーの直感を数学で定式化したのがマクスウェル方程式です。CAEの数式の裏には、常に「物理的な直感」があることを忘れずに。
離散化手法の詳細解説
空間離散化における手法選択が数値精度・安定性・計算コストに与える影響を詳述する。
辺要素（Nedelec要素）
電磁場解析に特化した要素。接線成分の連続性を自動的に保証し、スプリアスモードを排除。3D高周波解析の標準。
節点要素
スカラーポテンシャル定式化に使用。静磁場のスカラーポテンシャル法や静電場解析で有効。
FEM vs BEM（境界要素法）
FEM: 非線形材料・非均質媒質に対応。BEM: 無限領域（開領域問題）を自然に扱える。ハイブリッドFEM-BEMも有効。
マトリクスソルバーの選定指針
問題規模と特性に応じた最適なソルバー選択のガイドライン。
| ソルバー種別 | 詳細・推奨条件 |
|---|---|
| 直接法 | 小〜中規模の電磁場問題に適する。LU分解が主流。周波数掃引で右辺のみ変更可能な場合に効率的。 |
| 反復法（GMRES、BiCGSTAB） | 大規模3D電磁場問題で必須。前処理にはILU(k)やAMGが有効。複素対称系にはCOCG法も選択肢。 |
| DOF別推奨 | 〜5×10⁴ DOF: 直接法、5×10⁴〜: GMRES+ILU前処理、10⁶ DOF〜: 並列GMRES+AMG |
時間積分法と収束判定
ソルバー内部の制御パラメータと収束判定基準について記述する。
非線形収束（磁気飽和）
B-Hカーブの非線形性をニュートン・ラフソン法で処理。残差基準: $||R||/||R_0|| < 10^{-4}$が一般的。
周波数領域解析
時間高調波仮定により定常問題に帰着。複素数演算が必要だが、広帯域特性は時間領域解析で取得。
時間領域の時間刻み
最高周波数成分の1/20以下の時間刻みが必要。暗黙的時間積分ではより大きな刻みも可能だが精度に注意。
数値解法の直感的理解
辺要素のイメージ
電磁界解析で使われる辺要素（エッジ要素）は「道路に沿った交通量」をイメージすると分かりやすい。通常の節点要素が「交差点の値」を扱うのに対し、辺要素は「辺（道路）に沿った量」を直接扱う。電場や磁場は方向を持つベクトル量であり、辺に沿った成分として自然に表現できるため、スプリアス解（非物理的な偽の解）を自動的に排除できる。
周波数領域と時間領域の使い分け
周波数領域解析は「ラジオの特定の周波数に合わせる」ようなもの——1つの周波数での応答を効率的に計算できる。時間領域解析は「全チャンネルを同時に録画する」ようなもの——あらゆる周波数成分を含む過渡現象を再現できるが計算コストが高い。
電磁界解析の精度と計算コストの両立は永遠の課題です。 — Project NovaSolverは、既存ワークフローの改善を目指す取り組みとして、この問題に向き合っています。
CAEの未来を、実務者と共に考える
Project NovaSolverは、電磁力における実務課題の本質に向き合い、エンジニアリングの現場を支える道具づくりを目指す研究開発プロジェクトです。
プロジェクトの最新情報を見る →