減磁曲線 — 数値解法と実装

カテゴリ: 電磁場解析 | 2026-01-20
demagnetization-curve-method
数値解法の舞台裏

FEMでの減磁曲線の扱い

🧑‍🎓

減磁曲線をFEMにどう組み込むんですか?


🎓

永久磁石は等価電流モデルで扱う。磁石の構成式:


$$ \mathbf{B} = \mu_0 \mu_r \mathbf{H} + \mu_0 \mathbf{M}_0 $$

$\mathbf{M}_0$: 残留磁化ベクトル。減磁曲線の傾きがリコイル透磁率$\mu_r$(NdFeBで約1.05)。


🧑‍🎓

不可逆減磁はどう判定するんですか?


🎓

動作点が減磁曲線の膝点(ニーポイント)を下回ると不可逆減磁が発生する。JMAGやMaxwellでは各要素の動作点を減磁曲線上にプロットして膝点以下の領域を可視化できる。温度分布込みの減磁解析では各要素に温度依存の減磁曲線を割り当てる。


まとめ

🎓
  • 等価電流モデル — $\mathbf{B} = \mu_0(\mathbf{H} + \mathbf{M}_0)$
  • 膝点判定 — 不可逆減磁の評価基準
  • 温度連成 — 各要素に温度対応の減磁曲線を使用

Coffee Break よもやま話

電気自動車モータ開発と電磁界解析

テスラのModel 3のモータは、リラクタンストルクと磁石トルクの両方を使うIPMSM(埋込磁石型同期モータ)。この複雑な磁場分布を最適化するには数千回の電磁界FEA解析が必要です。1回の解析に数分としても、最適化ループ全体では数週間のCPU時間。それでも実機を何十台も試作するよりは圧倒的に速くて安い。

離散化手法の詳細解説

空間離散化における手法選択が数値精度・安定性・計算コストに与える影響を詳述する。

辺要素(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は、既存ワークフローの改善を目指す取り組みとして、この問題に向き合っています。

減磁曲線の実務で感じる課題を教えてください

Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。

実務課題アンケートに回答する →