高周波焼入シミュレーション

誘導焼入れの物理と理論基礎

表皮効果が全てを決める

誘導焼入れは、コイルの交流磁場がワーク表層に誘起する渦電流のジュール熱で表面だけを急速加熱し、直後の急冷でマルテンサイト硬化層を作るプロセスです。加熱の空間分布を支配するのが表皮効果で、電流密度は表面から浸透深さ \( \delta \) で指数減衰します。

$$ \delta = \sqrt{\frac{2\rho}{\mu \omega}} = \sqrt{\frac{\rho}{\pi \mu f}} $$

鋼の室温では \( f = 10 \) kHzで \( \delta \) は0.1 mm台、\( f = 1 \) kHzで1 mm弱——周波数の選択がそのまま硬化深さの設計変数です。ここに誘導加熱特有のツイストが加わります：温度が上がると抵抗率 \( \rho \) は増え、キュリー点（鋼で約770 °C）で強磁性が消えて \( \mu \) が激減、浸透深さが一気に数倍へ跳ねます。つまり加熱中に「加熱される場所」自体が動的に変わる——この強い非線形が、誘導焼入れシミュレーションを単なる熱解析と分ける核心です。

電磁・熱・組織の3物理連成

電源周波数と加熱時間のスケール差が6桁あるため、電磁場は周波数領域（時間調和）解析で1周期平均の発熱密度を求め、それを熱解析のソースにする分離が標準です。温度が変わるたびに物性（ρ, μ）を更新して電磁場を解き直す弱連成反復で全体が構成されます。

焼入れ側の物理——何をもって「硬化」とするか

硬化の条件は2段です。①加熱でオーステナイト化温度（Ac3、昇温速度依存で上方シフト）を超えること、②冷却速度が臨界冷却速度を上回りマルテンサイト変態（Ms点以下でKoistinen-Marburger則）すること。したがってシミュレーションの成果物は温度のコンターではなく、「オーステナイト化した領域∩急冷できた領域」＝硬化層パターンと硬さ分布です。冷却はスプレー/浸漬の沸騰熱伝達（温度依存の激しい沸騰曲線）が支配し、ここも強非線形です。

数値解析の手法

電磁-熱連成の実装

1. 電磁場（調和解析） — コイル電流（または電源電力）を与え、渦電流と発熱密度分布を計算。メッシュは表皮深さ内に最低2〜3要素（キュリー前の薄いδ基準）が鉄則
2. 熱解析（過渡） — 発熱密度をソースに温度を進める
3. 物性更新と再解析 — 温度上昇でρ・μ（B-H曲線）を更新し電磁場を解き直す。更新間隔は温度変化（特にキュリー点通過）で制御
4. 冷却・変態解析 — 沸騰熱伝達で急冷し、CCT/TTT情報に基づく変態モデル（拡散型はJMAK系、マルテンサイトはKM式）で組織分率と硬さを計算

収束と精度の急所

• キュリー点通過 — μの急落で発熱分布が激変。この区間は物性更新を細かく（温度刻み10〜20 °C）しないと加熱履歴を外す
• B-H非線形 — 表面近傍は磁気飽和領域で動作する。線形μ近似は発熱を過大評価。実測B-H曲線（温度依存）が精度の土台
• 移動（走査焼入れ） — コイルまたはワークが動く場合、位置ごとの調和解析＋移動熱源化、または移動座標系での定常化で扱う
• 電源のモデル化 — 実機は定電力/定電流制御＋整合回路。コイル電流一定の仮定が実機と合わない場合、電源制御をループに含める

冷却熱伝達の校正

スプレー冷却の熱伝達係数は表面温度依存（膜沸騰→遷移沸騰→核沸騰）で1桁以上変化し、文献値のばらつきも大きい——冷却側hは加熱側と独立に校正すべきパラメータです。実務では、試験片の埋込み熱電対による冷却曲線からhの温度依存カーブを逆同定します（リフローのh同定と同じ発想の、より非線形が強い版）。

実務適用の指針

標準ワークフロー

1. 材料データ整備 — 温度依存のρ・B-H曲線・熱物性、CCT線図、（変形まで見るなら）高温力学物性。ここが最大の準備工数
2. 加熱の校正 — 実機条件で試験片を加熱し、放射温度計/熱電対の温度履歴に電磁-熱モデルを合わせる（実効電流・効率の同定）
3. 冷却の校正 — 冷却曲線から沸騰hカーブを同定
4. 硬化パターン予測と検証 — 断面マクロ・硬さ分布（ビッカース走査）と比較。硬化深さ±10〜15%程度が校正済みモデルの相場
5. 条件探索 — 周波数・電力・加熱時間・コイル形状・冷却条件のスタディで狙いの硬化層・変形最小の条件を絞る

シミュレーションが効く典型場面

• コイル設計 — 段付き軸・歯車のような形状で、エッジの過熱（角部は磁束集中で焼き割れリスク）と谷部の未硬化を両立させるコイル形状・カップリング距離の検討。試作コイルの削減効果が大きい
• 新製品の条件出し — 類似品からの条件流用が効かない形状変更時の初期条件絞り込み
• 変形・割れ対策 — 変態膨張＋熱応力による変形・残留応力の予測。焼割れ（引張残留応力×急冷）のリスク部位特定

実測との突き合わせの作法

検証の基準データは「硬さ分布の断面マップ」が最強です。温度履歴の点計測より情報量が多く、加熱（オーステナイト化域）と冷却（臨界冷却速度達成域）の両方の妥当性を同時に検査できます。硬化深さだけでなくパターンの形状（エッジの回り込み・谷部の途切れ）まで合っているかを見ること——深さは合うがパターンが違う場合、電磁場側（コイルモデル・カップリング距離）の誤りを示唆します。

主要ツールの対応状況

ツール比較

選定の考え方

①電磁側の精度要求（コイル最適化までやるなら電磁専業級のB-H取り扱い）、②組織・硬さ・変形までの一貫性（DANTE級の材料モデルが要るか）、③使用者（解析専任かコイル設計現場か）で決めます。共通の注意として、ツールに内蔵の材料データベースは自社鋼種と一致しないことが多く、少なくともB-H曲線とCCT情報は自社データ整備が精度の前提になります。

先端動向

コイル形状の自動最適化

硬化パターンを目標に、コイル形状・巻き数・カップリング距離をパラメトリック/トポロジー的に最適化する研究・実装が進んでいます。電磁-熱連成1ケースが分オーダーに収まるようになったことで、実用的な反復回数の最適化が可能になりました。3Dプリントによる複雑コイル製作と組み合わせ、「計算で設計したコイルをそのまま作る」流れが現実になりつつあります。

デジタルツイン・プロセス制御との統合

電源の実測波形・赤外温度計測をリアルタイムでモデルと同化し、ばらつき（材料ロット・コイル摩耗・冷却剤劣化）を補償する制御研究が進行中です。物理モデルはROMやサロゲートに縮約してミリ秒応答にし、全数品質保証（各ワークの推定硬化深さをログ化）へつなげる構想が産業側の駆動力です。

組織・特性予測の高度化

急速短時間加熱に特有の非平衡オーステナイト化（Ac3の昇温速度依存・炭化物の溶け残り）を扱う速度論モデル、旧オーステナイト粒径と硬さ・疲労強度の関係、残留応力を含めた疲労寿命予測までの一貫連成が研究の前線です。「硬化深さ」から「部品性能（疲労限）」へ、予測のゴールが一段深くなる方向です。

トラブル対応

症状別の原因と対策

外れたとき、どの物理から疑うか

関連記事：産業熱プロセスの記事一覧、アーク溶接シミュレーション（同じ移動熱源・相変態の親戚）、リフロー炉シミュレーション。