Jiles-Athertonモデルで磁性材料のB-Hループをアニメーション描画。保磁力・残留磁束密度・ヒステリシス損失をリアルタイム計算します。
Jiles-Athertonモデルの核心は、外部から加えた磁界Hに、磁化M自身が作る磁界を足し合わせた「有効磁界」を考えることです。これにより、磁化の遅れ（ヒステリシス）を表現します。
$$H_e = H + \alpha M$$$H_e$: 有効磁界 [A/m], $H$: 外部印加磁界 [A/m], $M$: 磁化 [A/m], $\alpha$: 分子場係数（磁化の相互作用の強さ）
この有効磁界に対する、ピン留めなどの障害が一切ない理想的な磁化の応答を「無磁滞磁化」$M_{an}$と呼びます。これはランゲヴィン関数で表され、材料が本来持つ磁化の限界値（飽和磁化$M_s$）に近づいていきます。
$$M_{an}= M_s \left[ \coth\left(\frac{H_e}{a}\right) - \frac{a}{H_e}\right]$$$M_{an}$: 無磁滞磁化 [A/m], $M_s$: 飽和磁化 [A/m], $a$: 形状係数。この曲線が、シミュレーターで見られるループの「中心線」のような役割を果たします。
実際の磁束密度Bは、真空の磁束密度と材料の磁化による寄与を足し合わせたものです。ヒステリシス損失は、1サイクルで磁化Mを変化させるのに要した仕事量（ループの面積）として計算されます。
$$B = \mu_0 (H + M), \quad W = \mu_0 \oint H \, dM$$$B$: 磁束密度 [T], $\mu_0$: 真空の透磁率, $W$: ヒステリシス損失 [J/m³]。シミュレーター右下に表示される「Loss」がこの$W$に相当します。
変圧器・インダクタの鉄心：電力損失と発熱を抑えるため、ヒステリシスループが細く（保磁力Hcが小さく）、損失Wが最小の軟磁性材料（例：方向性ケイ素鋼板、アモルファス合金）が選ばれます。シミュレーターの「トランスコア」がこれに近い特性です。
永久磁石（モーター・スピーカー）：一度磁化したらその強さを長く保つ必要があるため、残留磁束密度Brと保磁力Hcが共に大きな硬磁性材料（例：ネオジム磁石、フェライト磁石）が使われます。シミュレーターの「ハード磁石」が該当します。
磁気記録媒体（HDD、磁気テープ）：情報のビットを微細な磁区の向きで記録します。書き込み時には局所的に磁化を反転させ、読み出し時にはその状態が安定に保持される必要があるため、適度な保磁力を持つ材料が研究されます。
非破壊検査・磁気シールド：材料の磁気特性の違いをヒステリシスループで評価し、内部の欠陥を検出したり、外部磁場から精密機器を守るシールド材の設計に活用されます。シミュレーターで材料を切り替える行為は、まさにこの材料選定プロセスに相当します。
このシミュレーターを使い始めるとき、いくつかハマりやすいポイントがあるよ。まず「保磁力Hcが大きい材料＝良い磁石」という短絡的な考え方。確かにハード磁石はHcが大きいけど、用途によって「良い」の基準は全く違うんだ。例えば、スピーカーの磁気回路には、ある程度Hcが大きくて温度変化に強い（減磁しにくい）材料が求められる。一方、モーターのコアはHcが小さい（ループが細い）軟鉄材料が必須で、ここにハード磁石を使ったら発熱だらけで回らない。材料選択は「何を実現したいか」から逆算してね。
次にパラメータ設定の落とし穴。シミュレーター上で「飽和磁化」や「分子場係数α」を自由に変えられるけど、現実の材料はこれらのパラメータが互いに独立じゃないんだ。例えば、フェライトのαをむやみに大きくしすぎると、非物理的なループ形状（例えば、角張りすぎたり）が出ることがある。実務でJ-Aモデルのパラメータを同定する時は、実測のB-Hループデータにフィッティングさせて、一組のパラメータ群として決定するのが鉄則だ。例えば、あるフェライト材では $M_s=3.2\times10^5$ [A/m], $a=50$ [A/m], $α=0.001$ といった具合に、バランスの取れた値のセットになるよ。
最後に、「シミュレーションのループ面積＝そのまま発熱量」と思い込むこと。確かにヒステリシス損失Wはループ面積に比例するけど、実際の機器での発熱は、これに「渦電流損失」が加わることを忘れちゃダメ。特に交流磁界中では周波数が高くなるほど渦電流損失が支配的になる。このツールで「トランスコア」を選んで周波数を上げてもループ形状は変わらない（渦電流を考慮してないから）けど、現実では薄い珪素鋼板を積層する理由は、まさにこの渦電流を抑えるためなんだ。
この磁気ヒステリシスの理解は、CAEの中でも特に電磁界解析（EM解析）の基礎中の基礎だ。具体的には、モーターやトランスのコア損失を精度良く見積もる「鉄損計算」に直結する。FEM（有限要素法）ソフトでモーターを解析する時、材料に設定するB-Hカーブや損失係数は、まさにこのJ-Aモデルで表現されたりするんだ。計算が雑だと、モーターの効率予測が大きく外れて、製品開発で痛い目を見る。
また、電力工学・パワーエレクトロニクス分野とも深く繋がっている。スイッチング電源に使われるトランスやインダクタのコア選定では、使用周波数帯域でのヒステリシス損失が効率を左右する。例えば、100kHzで動作するDC-DCコンバータには、ループが細く高周波特性に優れたフェライトコアが選ばれる。逆に、商用周波数（50/60Hz）の大型変圧器には、方向性電磁鋼板が使われる。このシミュレーターで材料を切り替えると、損失の桁がどう変わるか体感できるはずだ。
さらに応用先として非破壊検査（NDT）も挙げられるよ。磁性体材料の内部応力や疲労劣化は、磁気ヒステリシスループの形状を微妙に変化させる。これを測定して材料状態を評価する技術があるんだ。シミュレーターで「ピン留め係数k」を変えるとループの太さが変わるけど、あれが現実では材料内部の欠陥やひずみの増加に対応すると考えられる。つまり、B-Hループは材料の「健康診断書」にもなり得るんだ。
まず次の一歩としては、「動的ヒステリシス」を学ぶことをオススメする。このシミュレーターは準静的な変化（ゆっくり磁界を変えた時）をモデル化しているけど、実際の交流応用では周波数によるループの広がり（渦電流の影響）や、磁界変化の速度そのものが損失に与える効果を考える必要がある。数学的には、微分項 $dM/dt$ を含む拡張J-Aモデルや、損失を分離して計算する「Steinmetzの経験式」といった世界が待っている。
数学的背景をもう一歩深めたいなら、常微分方程式の数値解法に注目してみよう。J-Aモデルの核心は、磁化Mに関する微分方程式 $$ \frac{dM}{dH} = \frac{M_{an} - M}{k \delta - α(M_{an}-M)} $$ を解いていることだ（ここで$\delta$は磁界Hの増減方向の符号）。シミュレーターが滑らかなアニメーションを描く裏側では、この方程式を「オイラー法」や「ルンゲ・クッタ法」といった数値解法でコマ送り計算している。自分で簡単なスクリプトを書いて再現してみると、モデルの本質が手に取るように分かるよ。
最後に、材料側の視点を加える学習も大切だ。微視的な磁性物理学の基礎を知ると、パラメータの意味がより鮮明になる。例えば「飽和磁化$M_s$」は材料中の原子が持つ磁気モーメントの総和に、「分子場係数α」は磁区を形成させる交換相互作用の強さに、それぞれ対応している。シミュレーターで「フェライト」を選ぶとループが少し四角くなるのは、その微視的な磁化反転メカニズム（磁壁移動ではなく磁化回転が主）が背景にあるんだ。ツールのパラメータと物理現象の結びつきを想像できると、CAEエンジニアとしての応用力が一段階上がるはずだ。