ヒステリシスモデリング
ヒステリシスモデリングの理論基礎
磁気ヒステリシス
先生、ヒステリシスって磁化が磁界に追従しない現象ですよね？
強磁性体に交番磁界を加えるとB-H曲線がループを描く。このループの面積がヒステリシス損に対応する。
主要なモデル：
- Jiles-Atherton (J-A)モデル — 5パラメータ。磁区壁のピニングに基づく微分方程式モデル
- Preisachモデル — ヒステロンの重ね合わせ。任意形状のマイナーループを再現
- Play/Stopヒステロンモデル — 機械的ヒステリシスとの類似性
J-Aモデルが一番よく使われますか？
FEMへの組み込みやすさからJ-Aモデルが多い。ただし高精度にはPreisachモデルが優れる。JMAGではPlay modelを採用している。
まとめ
ヒステリシスの物理——磁壁ピンニングと「磁石の記憶」
磁性材料のヒステリシスループは「磁壁の移動と磁化回転」が外部磁場に遅れて応じる現象だ。磁壁が結晶粒界・介在物・欠陥に「ピンニング」されることで磁化の変化が不可逆的となり、保磁力（Hc）が生じる。ソフト磁性材料（電磁鋼板）は磁壁ピンニングが小さくヒステリシスループ面積（鉄損）が小さく、ハード磁性材料（ネオジム磁石）は大きなピンニングで強い保磁力を持つ。Preisach（1935）はヒステリシスを「無限多数のスイッチング要素の重ね合わせ」として数学的にモデル化し、現代のヒステリシスCAEモデルの礎を築いた。
ヒステリシスモデリングの数値計算手法
FEMへの組み込み
ヒステリシスモデルをFEMにどう実装するんですか？
各ガウス点（積分点）でB-H履歴を追跡する。時間ステップごとに：
1. 仮のBを計算
2. ヒステリシスモデルで対応するHを求める
3. 残差を評価してNewton-Raphson反復
計算コストが大きそうですね。
通常のB-H曲線（単一値）に比べて3〜5倍の計算時間になる。Preisachモデルでは各積分点にヒステロンの履歴データを保持するためメモリも増大する。JMAGではPlay modelにより効率的に実装されている。
まとめ
ヒステリシスモデルの数値実装——プライザッハモデルとJiles-Atherton
CAEでのヒステリシスモデリングには「プライザッハ（Preisach）モデル」と「Jiles-Athertonモデル」が主要な二手法だ。プライザッハモデルは測定した一次磁化曲線から密度関数を逆解析で求め、任意の磁化過程を正確に再現するが計算コストが高い。Jiles-Athertonモデルは5パラメータの微分方程式で磁化を記述し、商用FEMソルバへの組み込みが容易で計算が速い。磁気ヒステリシス損失の高精度予測にはプライザッハモデルが優れるが、モータ設計の実用解析ではJiles-Athertonが広く採用されている。
ヒステリシスモデリングの実務適用
実務での適用
鉄損の精密評価、変圧器の残留磁化解析、磁気シールドの交番磁界応答が主な用途。
実務チェックリスト
「計算鉄損と実測が合わない」——ヒステリシスモデルのキャリブレーション
電磁鋼板の鉄損計算でSimation（Steinmetz式）を使うと、実測との差が10〜30%になることがある。原因はSteinmetz係数がエポキシ処理・打ち抜き加工の影響で変化すること、および磁束密度の高調波成分がモデルに含まれないことだ。精度改善には①実機から切り出したEプシュー試料の測定データを用いたJiles-Athertonパラメータフィッティング、②打ち抜き端面影響の追加損失係数の導入が有効だ。JMAGは「磁気特性入力補助ツール」で測定データからパラメータを自動フィッティングする機能を持つ。
ヒステリシスモデリングのソフトウェア比較
ツール
| ツール | 特徴 |
|---|---|
| JMAG | Play modelベース。鉄損分離にヒステリシスモデルを活用 |
| Ansys Maxwell | J-AモデルとPreisachモデルの両方に対応 |
| COMSOL AC/DC | J-Aモデル内蔵。カスタムモデルも実装可能 |
| MagNet (Simcenter) | Preisachモデル対応。変圧器の突入電流解析に実績 |
ヒステリシス解析ツール——JMAG vs ANSYS Maxwell vs Opera FEM
磁気ヒステリシスのFEM解析ツールとして、JMAG（JMAGスチールライブラリの充実）、ANSYS Maxwell（EM＋回路連成）、Cobham Technical Services Opera FEM（モータ・発電機設計の老舗）が代表格だ。JMAGは日本製電磁鋼板メーカ各社の測定データをライブラリとして提供しており、日本の自動車・家電業界での採用が特に多い。Opera FEMは欧米の大型モータ・変圧器設計で実績があり、複雑なヒステリシスモデルのカスタマイズ性が高い。COMSOLのAC/DCモジュールも最近ヒステリシス損失機能を強化し、新規参入メーカへの訴求が続いている。
ヒステリシスモデリングの先端研究
先端技術
ベクトルヒステリシス——回転磁場下での鉄損予測の最前線
モータのステータコアでは回転磁場（磁束ベクトルが時間とともに方向を変える）が生じ、スカラーヒステリシスモデルでは鉄損を過小評価することが知られている。回転磁場下では鉄損が1軸磁化時より20〜40%増加する実験結果があり、「ベクトルヒステリシスモデル」が必要だ。Mayergoyzのベクトルプライザッハモデルや、Bergqvistのベクトルplay-modelが提案されているが、測定データ取得の困難さから実用化が遅れていた。近年のマルチアクシャル磁気特性測定装置の普及でベクトルモデルのキャリブレーションが現実的になっている。
ヒステリシスモデリングのトラブル対応
トラブル
「ヒステリシスループが閉じない」——数値シミュレーションの一般的な誤り
FEMヒステリシス解析で「ループが閉じない（スパイラル状になる）」現象は、時間ステップが大きすぎるときに起きる典型的な数値誤差だ。Jiles-Athertonモデルは微分方程式を時間積分するため、dHが大きすぎると積分誤差が蓄積してループが螺旋状にドリフトする。対策は①1サイクルあたり少なくとも100時間ステップの設定、②適応ステップサイズ制御（Runge-Kutta法）の使用。もう一つの落とし穴は「初期磁化状態」の設定ミスで、解析開始前に十分なプリサイクルを走らせてループを安定させる必要がある。
なった
詳しく
報告