熱機械疲労（TMF）解析

理論と物理

TMFとは何か

先生、TMFって普通の疲労と何が違うんですか？ 「熱」と「機械」の両方が入ってるのはなんとなく分かるんですけど…

いい質問だ。通常の疲労試験——例えば引張圧縮の繰り返し——は室温一定で行うよね。これを等温疲労（Isothermal Fatigue）と呼ぶ。一方、TMFでは温度サイクルと機械荷重が同時に作用する。ガスタービンブレードを想像してみてくれ。

ジェットエンジンの中のあの翼ですよね？ めちゃくちゃ高温になるやつ。

そう。起動時に室温から1000〜1200°Cまで急加熱され、同時に回転による遠心力で引張荷重がかかる。停止すれば冷却されて荷重も抜ける。このサイクルが何千回も繰り返されるわけだ。問題は、温度が変わると材料の降伏応力・弾性率・クリープ特性がすべて変わるということ。室温の疲労データだけで寿命を予測すると、現実と1桁以上ずれることがある。

え、1桁もずれるんですか？ それはさすがにまずいですね…

だからTMFという独立した評価体系が必要なんだ。TMFが問題になる典型的な部品を挙げると：

• ガスタービンブレード・ベーン：起動-停止サイクルで0〜1200°C
• 自動車エンジン排気マニホールド：冷間始動→高温定常→冷却で200〜900°C
• 発電用蒸気タービンロータ：起動時の熱過渡で表面と内部に温度差
• ブレーキディスク：制動時の摩擦加熱と冷却の繰り返し
• 原子力配管のサーマルストライピング：異温度流体の混合による温度揺動

In-phase と Out-of-phase

TMFの教科書を見ると「IP」とか「OP」って出てくるんですけど、これって何ですか？

TMFの最も重要な概念だ。温度と機械ひずみの位相関係によって、破壊モードが根本的に変わる。

• In-phase（IP, 同相）：引張ひずみが最大のときに温度も最高。高温・引張の組み合わせなのでクリープ損傷が支配的になる。粒界に沿ったき裂が内部から進展するパターンが多い。
• Out-of-phase（OP, 逆相）：圧縮ひずみが最大のときに温度が最高。高温で酸化層が形成され、冷却時に引張に転じると酸化層が割れて表面き裂が発生する。酸化損傷が支配的。

同じ材料・同じひずみ範囲でも、位相が違うだけでき裂の形態が変わるんですか？

その通り。例えばNi基超合金IN738LCの場合、同じ機械ひずみ範囲0.6%でもIP-TMFの寿命はOP-TMFの2〜5倍になることがある。IPではクリープ空孔が粒界に集積してゆっくり破壊するのに対し、OPでは酸化-き裂開口のサイクルが速く進む。実部品では位相が必ずしも0°（IP）や180°（OP）ではなく、中間的な角度を取ることも多い。これを菱形サイクル（Diamond cycle）と呼ぶんだ。

実際のガスタービンブレードだと、どの部位がIPでどの部位がOPになるんですか？

いい着眼点だ。ブレードの前縁は起動時に高温ガスに直接曝されて先に温度が上がり、翼内部はまだ冷たい。そのため前縁は圧縮ひずみ+高温→OP-TMF条件になりやすい。一方、ディスクのボア（内径部）は遠心力による引張ひずみと定常運転時の温度上昇が同期するため、IP-TMF条件になる。つまり同じエンジン内でもIPとOPが混在するんだ。

損傷メカニズムの重畳

普通の疲労だと繰り返し荷重でき裂が入る、というシンプルな話ですけど、TMFだともっと複雑なんですね？

TMFでは最低3つの損傷メカニズムが同時進行する：

1. 純疲労損傷（Fatigue）：機械ひずみの繰り返しによる転位蓄積→き裂核生成。低サイクル疲労（LCF）のメカニズムと同じだが、温度によって材料のサイクリック応答が変化する。
2. 酸化損傷（Oxidation）：高温保持中に表面酸化層が成長し、冷却時のひずみ不適合で酸化層にき裂が入る。Ni基合金ではAl₂O₃やCr₂O₃の保護膜が破壊されると酸化が加速する。
3. クリープ損傷（Creep）：高温・引張応力下で粒界すべりと空孔拡散が進む。長時間保持で損傷が蓄積し、粒界割れに至る。

この3つが線形加算で効くわけではなく、相互作用する。酸化で脆くなった粒界にクリープ空孔が集まれば、単独よりもはるかに速く破壊が進む。

支配方程式とひずみ分解

TMFの寿命予測ってどんな方程式を使うんですか？

TMFのすべてはひずみの分解から始まる。全ひずみ $\varepsilon_{total}$ から熱ひずみ $\varepsilon_{th}$ を差し引いて、寿命に直結する機械ひずみを取り出す：

ここで $\alpha(T)$ は温度依存の線膨張係数だ。機械ひずみはさらに弾性成分と非弾性成分に分かれる：

弾性・塑性・クリープの3つに分けるんですね。で、寿命はどうやって求めるんですか？

最も基本的なアプローチは修正Manson-Coffin式だ。等温疲労のManson-Coffin式を温度依存にしたもの：

各パラメータの意味はこうだ：

• $\sigma'_f(T)$：温度依存の疲労強度係数
• $b(T)$：疲労強度指数（通常 $-0.05$ 〜 $-0.12$）
• $\varepsilon'_f(T)$：疲労延性係数
• $c(T)$：疲労延性指数（通常 $-0.5$ 〜 $-0.7$）
• $E(T)$：温度依存のヤング率

ただし、この式だけでは酸化やクリープの影響を考慮できない。これが修正Manson-Coffin式の最大の限界だ。

Neu-Sehitogluモデル

酸化もクリープも入れた、もっと精密なモデルってあるんですか？

TMF寿命予測のゴールドスタンダードとされるのがNeu-Sehitogluモデルだ。1989年にイリノイ大学のNeuとSehitogluが提案した。損傷を3つのメカニズムに分離し、それぞれ独立に計算して足し合わせる：

それぞれの項を見ていこう。

1. 疲労項 $N_f^{fat}$：修正Manson-Coffin式そのもの。最高温度 $T_{max}$ での材料パラメータを使う場合が多い。

2. 酸化項 $N_f^{ox}$：酸化層の成長速度と、き裂開口のタイミングを考慮する。酸化層厚さ $h_{ox}$ は放物線則で成長する：

ここで $Q_{ox}$ は酸化の活性化エネルギー、$D_0$ は頻度因子、$R$ はガス定数だ。酸化損傷はOP条件で特に大きくなる——高温で酸化層が育ち、冷却時の引張で割れるから。

3. クリープ項 $N_f^{cr}$：高温保持中のクリープひずみ蓄積を積分する。Norton則を仮定すると：

クリープ損傷はIP条件で深刻になる。引張ひずみ最大時に温度が高く、クリープ変形が加速されるからだ。このように、Neu-SehitogluモデルはIPではクリープ項が支配的、OPでは酸化項が支配的という実験事実を自然に再現できる。これが1本の式（修正Manson-Coffin）では不可能だったことなんだ。

なるほど…損傷を分けて足すから、IPとOPで自動的に効く項が変わるわけですね。頭いい！

数値解法と実装

温度依存構成則

TMF解析をFEMでやるとき、材料モデルはどうするんですか？ 温度が変わると全部変わっちゃうんですよね？

TMF解析の成否は構成則（材料モデル）の精度でほぼ決まると言っていい。最低限必要なのは：

• 温度依存の弾塑性モデル：降伏応力、硬化則（等方硬化・移動硬化）のパラメータを複数温度で同定。Chabocheモデル（非線形移動硬化）がTMFでは事実上の標準。
• クリープ則：Norton則やノートン-ベイリー時間硬化則。高温保持の影響を再現する。
• 統合モデル（Unified Constitutive Model）：Bodner-Partom、Walker、Chaboche-Lemaitreなど。塑性とクリープを区別せず、粘塑性ひずみとして一体的に扱う。TMF向けには統合モデルが理想的だが、パラメータ同定が大変。

温度ごとにパラメータが違うってことは、データベースが膨大になりませんか？

そうなんだ。例えばChabocheモデルでパラメータが8個、それを6温度水準で同定すると48個のデータが要る。しかもサイクリック応力-ひずみ試験を各温度で実施する必要がある。航空エンジンメーカーでは専用の材料データベースを何十年もかけて構築している。これがTMF解析の最大の参入障壁と言ってもいい。

FEMによる熱-構造連成

FEMで温度と応力を同時に計算するんですよね？ 連成の方法って1つじゃないんですか？

大きく2つのアプローチがある：

1. 逐次連成（弱連成）：まず熱解析で温度場の時刻歴を求め、それを構造解析に温度荷重としてマッピングする。TMFの多くのケースではこれで十分だ。なぜなら、構造変形が温度場にフィードバックする効果は通常小さいから。

2. 完全連成（強連成）：熱と構造を同時に解く。塑性仕事による発熱（Taylor-Quinney係数）が温度場に影響する場合や、接触面の熱伝達が変形状態に依存する場合に必要。ブレーキディスクのTMFがこの例だ。

逐次連成の解析フローをもう少し具体的に言うと：

1. 熱伝導解析で起動→定常→停止の温度場時刻歴を計算（$T(\mathbf{x}, t)$）
2. 温度場を構造メッシュに転写
3. 各時刻ステップで $\Delta\varepsilon_{th} = \alpha(T)\Delta T$ を温度荷重として付与
4. 非線形構造解析（温度依存弾塑性＋クリープ）を実行
5. 安定化サイクルの応力-ひずみ-温度履歴を抽出 → 寿命予測モデルに入力

疲労ポスト処理

FEMの結果から寿命を出すところは、FEMソルバーの中でやるんですか？

現場では疲労専用のポスト処理ソフトを使うのが一般的だ。FEMソルバーが出力する応力・ひずみ・温度の時刻歴データを読み込み、各節点でTMF寿命を計算してコンター図を表示する。代表的なツール：

• fe-safe（Dassault Systemes）：Abaqusと密連携。TMF専用モジュールあり。Neu-Sehitogluモデル対応。
• nCode DesignLife（HBK/Siemens）：マルチソルバー対応。SWT・SRP・Manson-Coffin等を選択可能。
• FEMFAT（ECS）：自動車OEMで広く採用。高速処理。

なぜFEM内で完結しないかというと、疲労モデルは材料・環境・表面処理に応じて頻繁に更新が必要で、専用ソフトのほうがデータベース管理と手法切り替えの柔軟性が高いからだ。

サイクル外挿法

ガスタービンの寿命って数万サイクルですよね？ FEMで全サイクル計算したら何年かかるんですか…？

現実的に不可能だ。だからサイクル外挿法を使う：

• 安定化サイクル法：非線形FEMで数サイクル（5〜20サイクル程度）を計算し、応力-ひずみループが安定化（シェイクダウン）した時点のループを代表サイクルとして使う。Ni基超合金では3〜5サイクルで安定することが多い。
• 直接サイクルジャンプ法：AbaqusのDirect Cyclic機能。Fourier級数展開で安定化サイクルを直接求める。収束が速い。
• サブモデリング：全体モデル（グローバル）で温度・変位を求め、き裂危険部位をカットバウンダリで切り出して局所的にメッシュを細分化する。計算コストを1/10〜1/100に削減できる。

実践ガイド

TMF解析ワークフロー

実際にTMF解析をやるとき、最初の一歩から教えてください！

ステップバイステップで説明しよう：

Step 1. 運転条件の整理

• 起動→定常→停止→冷間保持の時間プロファイルを定義
• ガス温度・圧力・回転数の時刻歴を収集
• ミッションミックス（異なる運転パターンの組み合わせ）を決定

Step 2. 熱解析

• CFD（外部流れ）から熱伝達係数 $h$ と流体温度 $T_{\infty}$ を取得
• 内部冷却流路がある場合は流路側のHTCも設定
• 過渡熱伝導解析で温度場の時刻歴 $T(\mathbf{x}, t)$ を求める

Step 3. 構造解析

• 温度場を構造モデルにマッピング
• 遠心力・ガス圧力などの機械荷重を設定
• 温度依存の弾塑性＋クリープ構成則で非線形解析
• 安定化サイクル（5〜20サイクル）の応力-ひずみ履歴を抽出

Step 4. 寿命評価

• 疲労ポスト処理ソフトに応力/ひずみ/温度履歴を入力
• Neu-Sehitoglu等のTMFモデルで寿命マップを計算
• 最低寿命部位と安全率を報告

材料データの取得と管理

先生、TMF用の材料データってどこで手に入るんですか？ 論文を読んでも値がバラバラで…

これはTMF解析で最も苦労するポイントだ。信頼できるデータソース：

• NIMS（物質・材料研究機構）クリープデータシート：日本が誇る世界最大のクリープ試験データベース。無料でアクセス可能。
• MMPDS/CMH-17：航空宇宙材料の公式データハンドブック。米国規格。
• エンジンメーカー内部データベース：GE、P&W、Rolls-Royceなどは数十年分のTMF試験データを蓄積。社外秘。
• 自社試験：材料ロットのばらつきを考慮するなら、使用ロットからの試験が最も確実。

注意：論文の値はロット・試験片形状・試験条件が異なるため、そのまま使うと危険。特にクリープパラメータは試験温度範囲と応力範囲に敏感で、外挿すると大きく外れることがある。

メッシュ戦略

TMFのメッシュって、通常の構造解析と何か違いがありますか？

温度勾配が急な領域でメッシュが粗いと、ひずみ評価が大きくずれる。具体的な指針：

• 表面から3要素以内で温度変化の90%を解像する密度が必要。特にOP-TMFでは表面層の温度勾配が寿命を支配する。
• ブレード前縁・冷却孔周りには最低5層の境界層メッシュ。要素サイズは0.1〜0.3mm程度。
• フィレット・ノッチ部：応力集中と温度集中が重なる。六面体2次要素（C3D20R）が理想だが、複雑形状では四面体2次（C3D10）も許容。
• メッシュ収束性：3水準以上のメッシュ密度で安定化サイクルの機械ひずみ範囲が±5%以内に収束するか確認。

よくある失敗と対策

初心者がTMF解析でやりがちな失敗パターンってありますか？ 事前に知っておきたいです。

一番怖いのは「等温疲労データだけで済ませて安全だと思い込む」パターンですね…

その通り。等温LCFデータで1万サイクルの寿命が出ても、OP-TMF条件では2000サイクルで割れることがある。「等温≠TMF」——これをチーム全体で共有することが何より重要だ。

ソフトウェア比較

FEMソルバー比較

TMF解析に使えるFEMソルバーってどんなものがありますか？

TMFには温度依存の非線形弾塑性＋クリープ解析が必須だから、対応できるソルバーは限られる：

Abaqusの「Direct Cyclic」って便利そうですね。何が嬉しいんですか？

通常は安定化サイクルを得るために1サイクルずつ時刻歴解析を回す必要がある。Direct Cyclic法ではFourier級数展開を使って安定化後のループを直接求めるから、計算サイクル数が1/5〜1/10で済む。ただし、大きな塑性変形を伴うケースでは収束しにくいことがあるので注意が必要だ。

疲労寿命専用ツール

疲労専用ツールのfe-safeとnCodeはどう違うんですか？

先端技術

結晶塑性によるTMF

タービンブレードって単結晶合金を使ってますよね？ 普通の連続体モデルで大丈夫なんですか？

鋭い指摘だ。Ni基単結晶超合金（CMSX-4、PWA1484等）は結晶方位によって弾性定数や強度が異なる異方性材料だ。{001}方位と{111}方位では弾性率が2倍近く違う。だから精密なTMF解析には結晶塑性モデル（Crystal Plasticity Finite Element Method, CPFEM）が不可欠になりつつある。

CPFEMでは各すべり系ごとにせん断ひずみ速度をSchmid則で計算し、γ/γ'二相組織のラフト化（高温クリープによる析出相の方向性成長）も取り込める。計算コストは桁違いに増えるが、き裂発生の方位依存性を予測できる。

機械学習による寿命予測

最近は機械学習も使われてるんですか？

2つの方向で研究が進んでいる：

• サロゲートモデル：FEMの入出力関係をニューラルネットワークやGaussian Processで学習し、パラメータスタディを高速化。温度プロファイルの変更ごとにFEMを回さなくて済む。
• PINN（Physics-Informed Neural Network）：支配方程式を損失関数に組み込んだNNで、少数のTMF試験データから寿命予測モデルを構築。データ不足をフィジックスで補う。2024年以降、Ni基合金のTMF寿命予測で±30%精度の報告が出始めている。

ただし、学習データの範囲外（未知の温度条件・新合金）への外挿は依然として危険で、物理ベースモデルの置き換えにはまだ至っていない。現状は「物理モデル＋ML補正」のハイブリッドが現実的だ。

積層造形材のTMF

3DプリントのタービンブレードもTMF解析が必要ですか？

積層造形（AM）材は従来の鋳造材とは異なるミクロ組織（柱状晶・残留気孔・積層方向依存性）を持つため、TMF特性も大きく異なる。GE AviationがLEAPエンジンの燃料ノズルにAM部品を採用して以来、AM材のTMFデータ蓄積が急ピッチで進んでいる。主な課題：

• 残留気孔（特にLack of Fusionタイプ）がき裂核になりやすい
• HIP（熱間等方圧プレス）処理で気孔を潰しても、酸化物介在物が残る場合がある
• 積層方向（Build Direction）による異方性をモデルに取り込む必要がある
• 表面粗さの影響がOP-TMFで顕著（酸化促進）

トラブルシューティング

収束失敗

先生もTMF解析で徹夜デバッグしたことありますか？（笑）

何度もある（笑）。TMFは通常の構造解析より収束が難しい。高温域で材料が急に柔らかくなったり、クリープが急激に進んだりして非線形性が跳ね上がるんだ。よくあるケースと対策をまとめよう：

寿命予測と試験の乖離

解析では5000サイクル持つはずなのに、試験したら1000サイクルで割れた、みたいなことってあるんですか？

TMFでは珍しくない。寿命が5倍ずれた場合のチェックリスト：

1. 酸化の影響を入れ忘れていないか？ — Manson-Coffin式だけでは非保守的になりがち
2. 温度依存パラメータの外挿をしていないか？ — 800°Cまでのデータで1000°Cを予測していたら危険
3. 保持時間（dwell）の効果を考慮したか？ — 定常運転中の高温保持はクリープ損傷を蓄積する
4. 表面状態を反映したか？ — 加工粗さ、コーティング劣化、FOD（異物損傷）の影響
5. 機械ひずみの計算で熱膨張係数は正しいか？ — $\alpha$ の温度依存性が不正確だと、機械ひずみ範囲そのものが狂う
6. 多軸応力場を1軸に換算する手法は適切か？ — von Mises等価ひずみだけでは拘束効果を見落とす

ソルバー固有エラー

ソルバーのエラーメッセージが出たとき、TMF特有の対処法ってありますか？

Abaqus：

• ***WARNING: THE CREEP STRAIN INCREMENT IS TOO LARGE — CETOL（許容クリープひずみ増分）を小さく設定（0.001程度）。時間増分の上限も制限する。
• ***ERROR: EXCESSIVE DISTORTION AT ELEMENT — 高温で降伏応力が極端に低い領域。温度依存材料データの最高温度を確認し、外挿による非物理的な値がないかチェック。
• Direct Cyclic not converging — Fourier項数を増やすか、初期推定を改善（前回解からリスタート）。

Ansys Mechanical：

• Creep ratio exceeds limit — CRPLIM コマンドで最大クリープ比を調整。時間刻みの自動制御パラメータ（AUTOTS）を有効化。
• Negative pivot in element stiffness — 高温域での極端な材料軟化。弾性率の温度テーブルの下限値を確認（0にしない）。

TMF解析の全体像がやっと見えてきました。等温疲労とは全く別世界ですね。

そうだね。TMFは温度・機械荷重・時間の3つが絡む複合問題だから難しいけれど、正しく理解すれば非常にパワフルな設計ツールになる。実務では「まず等温データで当たりをつけ、TMF専用試験で検証し、解析モデルを校正する」という段階的アプローチが有効だ。いきなり完璧なモデルを作ろうとせず、段階的に精度を上げていこう。

「解析が合わない」と思ったら

1. まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
2. 等温・弾性の基準解を作る——最高温度で等温弾塑性解析を実行し、基準寿命を把握する
3. 1つだけ変えて再実行——温度依存性ON/OFF、クリープON/OFF、メッシュ粗密を「1つずつ」変えて影響因子を特定
4. 試験片レベルのベンチマーク——単軸TMF試験のFEMモデルを作り、実験のヒステリシスループと比較する。ここで合わなければ部品モデルの議論は無意味
5. 物理に立ち返る——「高温で引張なのにクリープ変形がゼロ」なら、クリープ則の入力が効いていない可能性大