永久磁石解析
理論と物理
永久磁石
先生、永久磁石はFEMでどう扱うんですか？
永久磁石は残留磁束密度$B_r$を持つ材料。外部電流なしで磁界を発生させる。
$$ \mathbf{B} = \mu_0 \mu_r \mathbf{H} + \mu_0 \mathbf{M}_r $$
先生、永久磁石はFEMでどう扱うんですか？
永久磁石は残留磁束密度$B_r$を持つ材料。外部電流なしで磁界を発生させる。
または等価的に:
$\mathbf{M}_r$: 残留磁化、$\mathbf{H}_c$: 保磁力。FEMでは磁石の領域に$B_r$と磁化方向を指定する。
主要な永久磁石材料
| 材料 | $B_r$ [T] | $H_{cJ}$ [kA/m] | $(BH)_{max}$ [kJ/m³] | 用途 |
|---|---|---|---|---|
| NdFeB（焼結） | 1.2〜1.5 | 800〜2500 | 300〜450 | モーター、発電機 |
| SmCo | 0.9〜1.1 | 600〜2000 | 150〜250 | 高温用途 |
| フェライト | 0.3〜0.4 | 200〜400 | 25〜40 | 安価、スピーカー |
| Alnico | 0.7〜1.3 | 40〜160 | 10〜80 | 計器、センサ |
NdFeBが圧倒的に強いですね。
EV/HEVのモーターにはNdFeBが不可欠。ただし希土類元素の価格変動リスクがあり、フェライトモーターの復権も研究されている。
まとめ
ネオジム磁石の誕生——佐川眞人が1982年に変えた電磁機器の世界
現代の高性能モータ・スピーカ・MRIに欠かせないネオジム磁石（Nd₂Fe₁₄B）は1982年に住友特殊金属の佐川眞人博士が発明した。最大エネルギー積（BHmax）がそれ以前のサマリウムコバルト磁石の2倍以上で、EV用モータの「小型高出力化」を一気に可能にした。唯一の弱点は高温での保磁力低下で、120〜150°C以上で急速に減磁する。この「温度特性」をCAEで正確に評価することが、EV・ハイブリッド車のモータ設計において不可欠な作業となっている。
各項の物理的意味
- 電場項 $\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ソルバーはこの変換を内部で自動的に行う。
ソルバー設定
まとめ
永久磁石の減磁解析——動作点がBH曲線の「膝」を越えるかを確認する
永久磁石の減磁解析では、各磁石要素の「動作点（B, H）」が減磁曲線の膝点（Knee Point）より上にあるかを確認する。FEMで全磁石要素のBとHを計算し、減磁曲線上に動作点をプロットすることで脱磁リスクを評価する。高温・大電流・逆磁場の組み合わせが最悪ケースとなる。解析手順は①最悪ケースの電流と温度でFEM計算、②各要素のB-Hを抽出、③温度別の減磁曲線と比較してマージン（余裕）を算出。JMAGとANSYS MaxwellはともにこのDemagnetization解析ワークフローを標準サポートしている。
辺要素（Nedelec要素）
電磁場解析に特化した要素。接線成分の連続性を自動的に保証し、スプリアスモードを排除。3D高周波解析の標準。
節点要素
スカラーポテンシャル定式化に使用。静磁場のスカラーポテンシャル法や静電場解析で有効。
FEM vs BEM（境界要素法）
FEM: 非線形材料・非均質媒質に対応。BEM: 無限領域（開領域問題）を自然に扱える。ハイブリッドFEM-BEMも有効。
非線形収束（磁気飽和）
B-Hカーブの非線形性をニュートン・ラフソン法で処理。残差基準: $||R||/||R_0|| < 10^{-4}$が一般的。
周波数領域解析
時間高調波仮定により定常問題に帰着。複素数演算が必要だが、広帯域特性は時間領域解析で取得。
時間領域の時間刻み
最高周波数成分の1/20以下の時間刻みが必要。暗黙的時間積分ではより大きな刻みも可能だが精度に注意。
周波数領域と時間領域の使い分け
周波数領域解析は「ラジオの特定の周波数に合わせる」ようなもの——1つの周波数での応答を効率的に計算できる。時間領域解析は「全チャンネルを同時に録画する」ようなもの——あらゆる周波数成分を含む過渡現象を再現できるが計算コストが高い。
実践ガイド
実務
IPMモーター、SPMモーター、リニアモーター、磁気カップリング、磁気浮上の設計。
チェックリスト
「量産品で減磁ばらつき」——磁石グレード選定と設計マージンの確保
EV用IPMモータの量産では、同一設計でも個体ごとに減磁マージンが異なるケースがある。原因は磁石のロット間ばらつきと、組み付け時の磁石配置精度の差だ。磁石の公称Brが±3%ばらつくと、最小動作点で設計値より動作温度が5〜10 K低くなる。解決策は①磁石グレードの公差を厳しく調達仕様に記載、②FEM感度解析で「Brが-3%のとき最悪ケースの動作点」を確認してマージンを設定、③組み立てライン検査でリング磁石の残留磁束密度を全数測定する体制の確立。シミュレーションによる公差設計が製造品質の核心だ。
解析フローのたとえ
モータの電磁界解析は「ギターの調律」に近い感覚です。弦の太さ（コイル巻数）とブリッジの位置（磁石配置）を調整して、最も美しい音色（効率の良いトルク特性）を引き出す。1つのパラメータを変えると全体のバランスが変わる——だからパラメトリックスタディが重要なんです。
初心者が陥りやすい落とし穴
「空気領域？ なんで空気をメッシュで切るの？」——初めて電磁界解析に触れた人がほぼ全員抱く疑問です。答えは「磁力線は鉄心の外にも広がるから」。解析領域を鉄心ぎりぎりにすると、行き場を失った磁束が壁に「ぶつかって」反射し、実際にはありえない磁束集中が起きます。部屋が狭すぎてボールが壁に跳ね返りまくる状態を想像してみてください。
境界条件の考え方
遠方の境界条件って地味ですが超重要です。「ここから先は無限に広がる空間」ということを数値的に表現する必要がある。設定を間違えると、まるで「見えない壁」があるかのように磁束が跳ね返されてしまいます。
ソフトウェア比較
ツール
| ツール | 特徴 |
|---|---|
| JMAG | モーターの永久磁石解析の業界標準 |
| Ansys Maxwell | 自動適応メッシュ |
| MotorCAD | 磁気回路ベース。高速パラメトリック |
| COMSOL | マルチフィジックス。熱-磁場連成 |
| FEMM | 2D。無償。初期検討に最適 |
永久磁石解析ツール——JMAG vs ANSYS Maxwell vs Opera
永久磁石モータ解析ツールの比較：JMAGはネオジム磁石の温度依存BH曲線データベース（信越化学・TDKなど国産メーカ対応）と脱磁計算機能が充実しており、国内EV・家電メーカの標準ツール。ANSYS Maxwellは永久磁石の異方性・温度感度をモデル化したPM Material Setupウィザードで設定が容易で、HPC並列で大型モータ全体を高速解析できる。Opera FEM（Cobham）は欧州・北米の大型発電機・リニアモータ設計で実績が厚く、複雑形状永久磁石のカスタム着磁パターン設定に強みがある。
選定で最も重要な3つの問い
- 「何を解くか」：永久磁石解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端
部分減磁の検出——均一減磁と不均一減磁でシミュレーションが変わる
永久磁石の減磁は「均一減磁（全体が一様に低下）」と「部分減磁（局所的に膝点を超える）」で挙動が異なる。均一減磁は線形なのでBrを下げた材料特性で再計算すれば対応できるが、部分減磁は局所的に非線形なBH変化が起きるため、各要素の減磁履歴を追跡する「経路依存解析」が必要だ。JMAGの「脱磁計算」機能は有限要素ごとに減磁状態を管理し、繰り返し動作後の磁石特性変化を累積計算できる。この機能は高調波電流や繰り返し過渡動作が生じるインバータ駆動モータの寿命評価に直結する。
トラブルシューティング
トラブル
「磁石を変えたらモータ特性が変わった」——調達変更時の再検証の重要性
ネオジム磁石サプライヤの変更・グレード変更により、モータ性能が想定外に変化した事例は多い。磁石の「Br・Hcjの温度係数」「磁化方向精度」「寸法公差」がメーカ・グレードによって異なり、同じ公称スペックでも実特性に差がある。量産中のサプライヤ変更では「変更前後の磁石実測BH曲線をFEMに入力して比較計算」することが最低限の検証だ。特に高温特性（120〜150°C）の差は公称スペックに現れにくいため、温度範囲での実測が欠かせない。JMAGのデータベースには複数メーカの磁石BHデータが収録されており比較が容易だ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——永久磁石解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告