鉄心飽和解析
理論と物理
鉄心の飽和
先生、鉄心の飽和ってどういう現象ですか？
磁界$H$を増やしても磁束密度$B$がほとんど増えなくなる現象。鉄の磁気モーメント（磁区）が全て揃うと、それ以上の磁化は不可能。
B-Hカーブ
非線形の$B$-$H$関係:
| 領域 | $\mu_r$ | 特徴 |
|---|---|---|
| 初期領域 | 100〜500 | 低磁界。可逆的 |
| 急上昇域 | 1000〜10000 | 磁壁移動が支配 |
| 飽和近傍（膝） | 10〜100 | 磁区回転 |
| 飽和域 | ≈1 | $B \approx \mu_0 H + B_{sat}$ |
飽和すると$\mu_r \approx 1$…空気と同じですか？
そう。飽和した鉄心は磁束を通さなくなる。磁束が漏れ出して効率低下、損失増大、発熱の原因になる。
代表的な飽和磁束密度
まとめ
鉄心飽和の物理——磁区構造がBH曲線の「膝」を作る
電磁鋼板の磁化曲線に現れる「膝（Knee Point）」は、磁区の壁移動が完了し磁化回転だけが残る段階への移行だ。ソフト鋼のBsat（飽和磁束密度）は約2.0〜2.1 T、フェライトは0.3〜0.5 Tで大きく異なる。変圧器やモータの鉄心設計では動作点を膝点以下に保つことが基本だが、軽量化・高出力化の要求でBを膝点付近に設定するケースが増えている。飽和域の非線形BH曲線を正確にFEMに組み込むことが鉄損・磁束分布の精度を左右し、Langrangeエレメント vs Nedelecエレメントの選択が重要になる。
各項の物理的意味
- 電場項 $\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
$\nu = \nu(|\mathbf{B}|)$（$\nu = 1/\mu$）が磁束密度に依存。Newton-Raphson反復で解く:
1. 初期推定（$\mu_r = 1000$等で線形解）
2. $|\mathbf{B}|$からB-Hカーブで$\nu$を更新
3. 接線透磁率$d\nu/d(B^2)$でヤコビアン（接線剛性）を構成
4. 残差がゼロになるまで反復
構造の弾塑性と同じ非線形ソルバーですね。
まさにそう。B-Hカーブが応力-ひずみカーブに対応する。
B-Hカーブの入力
鋼板メーカーのデータシート（磁化曲線）をテーブルで入力。注意:
- 飽和域の延長 — データが$B = 1.8$ Tまでしかない場合、$\mu_0$の傾きで延長
- 補間方法 — スプライン or 線形。急激な折れ曲がりは収束に悪影響
まとめ
非線形磁場解析の反復法——Newton-Raphson法とBH曲線の扱い
鉄心飽和を考慮した磁場解析は非線形問題であり、Newton-Raphson法（NR法）が標準的な反復解法だ。NR法は収束が二次的（誤差の二乗で減少）と速いが、初期値が悪いと発散する欠点がある。BH曲線の「急勾配部分（膝付近）」では透磁率の変化が大きく、ヤコビアン行列の計算精度が収束性に直結する。実用上は線形初期解を与えてからNRに移行する二段階法や、Picard反復（固定点法）との組み合わせで安定性を高める手法がANSYS Maxwellなどの商用ツールでも採用されている。
辺要素（Nedelec要素）
電磁場解析に特化した要素。接線成分の連続性を自動的に保証し、スプリアスモードを排除。3D高周波解析の標準。
節点要素
スカラーポテンシャル定式化に使用。静磁場のスカラーポテンシャル法や静電場解析で有効。
FEM vs BEM（境界要素法）
FEM: 非線形材料・非均質媒質に対応。BEM: 無限領域（開領域問題）を自然に扱える。ハイブリッドFEM-BEMも有効。
非線形収束（磁気飽和）
B-Hカーブの非線形性をニュートン・ラフソン法で処理。残差基準: $||R||/||R_0|| < 10^{-4}$が一般的。
周波数領域解析
時間高調波仮定により定常問題に帰着。複素数演算が必要だが、広帯域特性は時間領域解析で取得。
時間領域の時間刻み
最高周波数成分の1/20以下の時間刻みが必要。暗黙的時間積分ではより大きな刻みも可能だが精度に注意。
周波数領域と時間領域の使い分け
周波数領域解析は「ラジオの特定の周波数に合わせる」ようなもの——1つの周波数での応答を効率的に計算できる。時間領域解析は「全チャンネルを同時に録画する」ようなもの——あらゆる周波数成分を含む過渡現象を再現できるが計算コストが高い。
実践ガイド
実務
モーター、変圧器、リアクトルの鉄心設計。動作点が飽和域に入らないように設計。
チェックリスト
「変圧器が設計より熱い」——局所飽和が引き起こす予期せぬ損失
変圧器やリアクトルの設計で「実機の温度上昇が計算値より10〜15 K高い」問題が起きるとき、鉄心の「局所飽和」が見落とされていることが多い。コーナー部・継ぎ目部では磁束が集中して局所的に飽和し、渦電流損失が急増する。断面平均の磁束密度がBsat以下でも局所的に超えているケースがある。FEM解析で磁束密度のカラーマップを確認すると、コーナー部の高密度領域が一目瞭然だ。対策としてコーナーの積層方向変更（斜め切り）や継ぎ目位置の最適化が有効で、FEMで事前検証してから製作する手順が定着している。
解析フローのたとえ
モータの電磁界解析は「ギターの調律」に近い感覚です。弦の太さ（コイル巻数）とブリッジの位置（磁石配置）を調整して、最も美しい音色（効率の良いトルク特性）を引き出す。1つのパラメータを変えると全体のバランスが変わる——だからパラメトリックスタディが重要なんです。
初心者が陥りやすい落とし穴
「空気領域？ なんで空気をメッシュで切るの？」——初めて電磁界解析に触れた人がほぼ全員抱く疑問です。答えは「磁力線は鉄心の外にも広がるから」。解析領域を鉄心ぎりぎりにすると、行き場を失った磁束が壁に「ぶつかって」反射し、実際にはありえない磁束集中が起きます。部屋が狭すぎてボールが壁に跳ね返りまくる状態を想像してみてください。
境界条件の考え方
遠方の境界条件って地味ですが超重要です。「ここから先は無限に広がる空間」ということを数値的に表現する必要がある。設定を間違えると、まるで「見えない壁」があるかのように磁束が跳ね返されてしまいます。
ソフトウェア比較
ツール
全磁場FEMソルバーが非線形B-Hカーブに対応。JMAG、Maxwell、COMSOL、FEMMの全てで飽和解析が可能。JMAGとMaxwellは電磁鋼板の材料データベースが充実。
鉄心飽和解析ツール——ANSYS Maxwell vs JMAG vs FEMM
鉄心飽和を含む磁場解析の定番ツールはANSYS Maxwell、JMAG、FEMM（オープンソース）だ。Maxwellは大規模3D非線形磁場解析と回路連成が得意で、変圧器・モータの統合設計に使われる。JMAGは日本製電磁鋼板の豊富なBHデータベースと鉄損計算精度で国内自動車・電機業界での採用が圧倒的だ。FEMMは2D軸対称・平面問題に限定されるが無料で、教育・研究・簡易設計に広く使われている。シミュレーション精度の決め手は計算手法よりも「BH曲線データの品質」だという事実は、ツール選定より材料測定に投資すべきことを示唆している。
選定で最も重要な3つの問い
- 「何を解くか」：鉄心飽和解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端
動的飽和——インバータ高調波が引き起こす追加鉄損
インバータ駆動モータではPWM高調波が電流に重畳し、鉄心が「主波数の飽和」に加えて「高調波飽和」を繰り返す。主波と高調波の相互作用で鉄損が最大30〜50%増加する事例が報告されており、正弦波駆動と直接比較すると実機の温度差がはっきり現れる。この「動的飽和鉄損」をFEMで正確に計算するには非線形時間ステップ解析が必要で、計算時間が定常周波数解析の10〜100倍になる。JMAGの「鉄損解析専用ソルバ」やANSYS Maxwellの「Transient Solver with Core Loss」がこの問題への実用的解決策を提供する。
トラブルシューティング
トラブル
「FEMが収束しない」——BH曲線の入力ミスと飽和域の発散
非線形磁場解析でFEMが収束しない原因の定番は「BH曲線の入力データの問題」だ。測定データのノイズで単調増加でない点が混在していると、反復解法中に透磁率が負になり発散する。対策は①単調増加を確認するスムージング処理、②飽和領域を十分高い磁場値まで外挿する（Bs後の直線延長）、③磁場が設定範囲外に出たときの警告設定の確認。JMAGは「BH曲線チェック機能」で単調性違反を入力段階で検出する。さらに解析領域の外側（空気領域）のメッシュが粗すぎても磁束ループが形成されず収束しないため、空気領域のサイズと密度も要チェックだ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——鉄心飽和解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告