構造-熱連成｜温度-応力の相互作用を高精度解析

2025.08.31

構造-熱連成

構造-熱連成は、温度場の変化が熱膨張や材料物性の温度依存性を通じて応力・ひずみを生み、その変形が逆に熱の伝わり方や境界条件に影響を及ぼす現象を同時に扱う解析である。熱伝導・対流・放射を記述するPDEと、力のつり合いを記述する連立方程式を同一フレームで解く点に特徴がある。小変形弾性域では古典的熱弾性、溶接やAMのような高温・大ひずみでは弾塑性・クリープ・相変態を含む強い非線形が重要となる。産業ではエレクトロニクス実装、タービン翼、ブレーキディスク、航空宇宙の加熱環境などで設計上の中核となる。

Table Of Contents

物理モデルと支配方程式

温度場は熱拡散方程式（フーリエの法則）で記述し、機械場は運動方程式（準静的つり合い）で記述する。熱ひずみはεth=α(T)ΔTで与え、全ひずみを弾性・塑性・熱成分に加算する。物性値E(T), α(T), k(T), cp(T), ρ(T)は温度依存を考慮する。強い連成では、変形に伴う接触開離により接触熱伝達係数や放射の形態係数が変わり、また繊維配向や残留応力が有効熱伝導率を変調する。非定常過程では内熱源（レーザ・ジュール発熱・化学反応）も連成源となる。

連成の分類（弱連成と強連成）

  • 弱連成：熱解析→構造解析の一方向。温度分布を荷重としてマップし、構造側から熱側へのフィードバックは無視する。定常熱応力評価などで有効。
  • 強連成：構造側の状態（接触、変形、相変態）が熱側に影響する双方向。溶接残留応力、L-PBFのレイヤヒーティング、ギャップ放射などで不可欠。
  • 時間スケール：定常/非定常、空間スケール：マクロ/ミクロを整理して支配的機構を同定する。

数値解法（FEM中心）

解法は大きくmonolithic（単一の連立方程式を同時解）とpartitioned（staggeredで交互反復）に分かれる。前者は収束性に優れるが線形方程式が巨大化し、前処理と並列計算が鍵となる。後者は既存ソルバの活用とモデル分割が容易で、Aitken緩和やIQN-ILSなどの界面緩和で収束を安定化する。非線形解ではNewton–Raphson、線形解ではAMGやILU前処理CG/GMRESを用いる。時間離散は熱側が拡散型ゆえ暗示的θ法が堅牢で、構造側はNewmark法などの定式化を用いる。

メッシュと時間積分の注意

  • 熱勾配の急峻部（熱源近傍・薄板端部・接触縁）と応力集中部（フィレット・孔縁）は局所細分化する。熱境界層/機械境界層を同時に解像する。
  • 要素選択は熱側に一次/二次四面体・六面体、構造側にひずみロッキングを抑える低減積分や混合要素を検討する。
  • 非定常はFourier数Fo=αΔt/Δx²を指標に時間刻みを選定し、過度な大刻みは熱応力のピークを見逃すため注意する。

境界条件と接触の取り扱い

対流はh(T)・流速依存、放射はεσ(T⁴−Tenv⁴)を用い、放射の非線形性は反復で扱う。機械境界は拘束・荷重・接触（摩擦）を定義する。接触熱伝達係数は面圧・粗さ・介在流体で変化し、開離時はギャップ放射と微小伝導のみとなる。熱流束や温度履歴を測定データから同化する場合は逆解析（Tikhonov正則化等）を併用する。

材料モデルとデータ整備

金属ではE(T)とσy(T)の低下、α(T)の上昇、k(T)とcp(T)の顕著な温度依存が支配的である。高温ではクリープ（Norton則）、相変態潜熱、熱軟化を考慮する。樹脂・封止材は粘弾性で時間温度換算（WLF）を用いる。複合材は異方性熱膨張・熱伝導のテンソル化が必要で、層間熱抵抗も重要となる。信頼できる物性表の収集と外挿範囲の明示が解析の再現性を左右する。

検証・妥当性確認（V&V）

  • コード検証：熱弾性の解析解（平板の一様加熱、二層バイメタル曲げ）で誤差収束を確認する。
  • 解の検証：メッシュ・時間刻み独立性（GCI等）を実施し、エネルギーバランス（入力熱＝蓄熱＋損失）を点検する。
  • 妥当性：熱電対・IRカメラ・DICの実験で温度・ひずみ履歴を突合し、不確かさ伝播を評価する。

代表的な適用例

  • 半導体実装：BGA/TSVのサイクル試験における熱疲労寿命と反り予測。
  • 溶接・AM：移動熱源に伴う残留応力・歪み、割れ感受性の評価。
  • 動力機器：タービン翼の急速加熱時の熱応力とクリープ損傷。
  • 輸送機器：ブレーキディスクの熱割れ、ロータのサーマルバウ。
  • ガラス・セラミックス：急冷テンパリング時の内部応力分布。

設計・解析ワークフロー

  1. 目的関数と評価指標（最大主応力、塑性ひずみ、変形、寿命）を定義する。
  2. 幾何・メッシュ・BCを整備し、弱連成/強連成の要否を判定する。
  3. 解法（monolithic/partitioned）とソルバ設定、前処理、並列化（CPU/GPU）を選定する。
  4. ポスト処理で温度履歴、熱ひずみ、残留応力、座屈リスクを可視化し、感度解析やDoEで設計変数を最適化する。

実務の勘所とトラブル対策

  • 初期温度・接触状態の与え方で収束が大きく変わる。段階載荷と物性のスムージングで安定化する。
  • 大温度勾配では要素の歪みが出やすく、二次要素や局所リメッシュが有効である。
  • 結果の解釈は温度依存降伏を踏まえた等価応力と主応力の併用が有効で、熱疲労はレンジと保持の両方を見る。
  • モデル簡略化（対称・2D化）と重要部のみの詳細化を併用し、計算資源を最適配分する。

構造-熱連成は、熱設計と強度設計を同一土俵で最適化するための基盤技術である。現象の支配因子（熱源、境界、材料、接触）を階層的に整理し、物性と境界の不確かさを管理しつつ、妥当性の確保と設計最適化を反復することが、信頼性・軽量化・コストの同時達成に直結する。

コメント（β版）