熱-構造連成設計（Thermo-Mechanical Design）

いや、熱応力の評価はほんの一側面に過ぎない。もっと大事なのは熱変形による寸法精度の崩れだ。例えば半導体露光装置のステージは0.1nm精度が求められるんだけど、温度変動がたった0.01°Cあるだけで数nmの変形が出る。

そう。だからインバー（Fe-Ni低膨張合金、CTE ≈ 1.2×10⁻⁶/K）のような特殊材料の選定、冷却流路の最適配置、拘束条件の設計まで含めて「熱-構造連成設計」と呼ぶ。単に $\sigma = E \alpha \Delta T$ を計算して終わりじゃない。

代表的なのはこの5分野だ：

数字で見ると怖さがわかるよ。例えば長さ1mの鋼材（CTE ≈ 12×10⁻⁶/K）の温度が10°C上がると：

0.12mm——普通の機械設計なら問題ないけど、工作機械の加工精度がμm級だと致命的だ。しかも温度分布が一様じゃないから、単純な膨張じゃなく「反り」や「ねじれ」が出る。これが熱変形の本当の怖さだ。

その通り。だから設計段階で「温度分布の非一様性をいかに小さくするか」と「残った非一様性にどう耐えるか」を両面から考える。前者が冷却設計、後者が材料選定と拘束設計だ。

そう。例えば電子基板では、シリコンチップ（CTE ≈ 2.6×10⁻⁶/K）がFR-4基板（CTE ≈ 14×10⁻⁶/K）の上にはんだ付けされている。温度サイクルでCTE差がそのまま界面せん断応力になるんだ。

ここで $\gamma$ ははんだのせん断ひずみ、$L$ はチップの半幅、$h_{solder}$ ははんだ高さ。CTE差が大きいほど、チップが大きいほど、はんだが薄いほどせん断ひずみは大きくなる。これがBGAはんだクラックの根本原因だ。

鋭いね。そしてタービンブレードでも同じ問題がある。Ni超合金の上にTBC（熱遮蔽コーティング、セラミック系、CTE ≈ 10×10⁻⁶/K）を吹き付けるんだけど、Ni超合金のCTEは約13×10⁻⁶/K。この差が繰り返しの起動停止でTBCの剥離を引き起こす。

まず熱伝導方程式（過渡）：

ここで $\rho$ は密度、$c_p$ は比熱、$k$ は熱伝導率、$Q$ は内部発熱。そして構造の平衡方程式：

連成のキモは構成則だ。全ひずみから熱ひずみを引いて弾性ひずみを求める：

基本はそう。ただし材料の弾性率 $E(T)$ や降伏応力 $\sigma_y(T)$ が温度依存する場合は、構造側の挙動が温度に依存するので「片方向連成」でも温度の影響が正確に反映される。さらに構造変形が接触状態を変えて熱抵抗が変わる場合（例：ボルト締結部の接触熱伝導）は双方向連成が必要だ。

それが大きな落とし穴だ。過渡状態のほうが定常より応力が大きいケースは非常に多い。蒸気タービンのロータを例にしよう。起動時に高温蒸気が外表面に当たると、外面はすぐ温まるが内部はまだ冷たい。このとき外面に圧縮、内面に引張の大きな温度勾配が生じる。

定常に達すると温度勾配が緩和されて応力は下がる。つまり起動時と停止時の過渡状態こそが最厳しい条件になる。実際に蒸気タービンメーカーは起動パターン（cold start / warm start / hot start）ごとの過渡解析を必ず行って、許容昇温レート（°C/min）を決めている。

だから熱-構造連成設計ではライフサイクル全体（起動→定常→停止→再起動）の過渡解析が必須なんだ。これが単なる「熱応力計算」と本質的に違うポイントだよ。

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

典型例は鍛造や摩擦攪拌接合（FSW）だ。塑性変形による発熱 $Q_{plastic} = \eta \cdot \sigma : \dot{\varepsilon}^p$（$\eta$ は塑性仕事の熱変換係数、通常0.9）が温度を上げ、温度が材料特性を変えて変形挙動が変わる——完全にフィードバックしている。Abaqusの `*Coupled Temperature-Displacement` がこの強連成を扱う代表的なステップだ。

構造の有限要素方程式は：

ここで熱荷重ベクトル $\{\mathbf{F}_{th}\}$ は各要素の熱ひずみから求める：

$\mathbf{B}$ はひずみ-変位マトリクス、$\mathbf{D}$ は弾性マトリクス。つまり「温度差 → 熱ひずみ → 等価節点力」という変換で、熱の効果が機械的な荷重と同じ形式で構造方程式の右辺に入る。同じソルバーで解けるから実装が楽なんだ。

熱伝導方程式のFEM離散化で後退Euler法を使うと：

$[\mathbf{C}]$ は熱容量マトリクス、$[\mathbf{K}_T]$ は熱伝導マトリクス。後退Euler（陰解法）は無条件安定だから時間刻みを大きくできるけど、過渡的な温度ピークを捕まえるには十分小さい時間刻みが必要だ。目安は：

これは実務で本当に重要なポイントだ。熱解析ではメッシュが粗くても収束するけど、構造では応力集中部に細かいメッシュが必要になる。メッシュが非一致の場合、温度データの転写（マッピング）が必要で：

Ansys Workbenchなら `Imported Load` で自動マッピング、Abaqusなら `*TEMPERATURE, FILE=job_thermal` で熱結果ファイルを直接読み込める。ただし転写時のエネルギー保存性を必ず確認すること。転写誤差が大きいと、ないはずの熱応力が出たり、あるはずの応力が消えたりする。

熱-構造連成設計の標準的なワークフローは6ステップだ：

その通り。よくある失敗は「定常運転だけ解析して、起動時の過渡応力を見落とす」パターン。蒸気タービンの例で言えば、cold start 時の昇温レートが設計の制約条件になることが多い。最も過酷な瞬間がいつなのかをまず特定する。

CTEだけじゃ不十分だ。熱変形感度は複数の材料特性で決まる。よく使われる指標を紹介しよう：

具体例で比べてみよう：

いいところに気づいたね。アルミは「膨張するけど温度が均一化しやすい」から、温度の非一様性（反り・ねじれ）は小さくなりやすい。一方インバーは膨張しないけど熱伝導が悪いから温度ムラが残りやすい。どちらが有利かは用途による。だからFEMで具体的な温度分布と変形を計算する必要がある。

冷却設計と拘束設計はセットで考える。ポイントは3つ：

剛体の6自由度（3並進 + 3回転）を最小限の点で拘束する方法だ。1面3点で面直方向を拘束（3自由度消去）、1辺2点で面内回転と1並進を拘束（2自由度消去）、1点で残りの1並進を拘束。こうすると熱膨張に対して余計な拘束力が発生しない。過拘束（例えば4隅を完全固定）すると、膨張しようとする材料が拘束されて不必要な熱応力が発生する。

当然だ。熱-構造連成設計の最終目標は「設計寿命を満足するか」だからね。特に温度サイクルが繰り返される環境では熱機械疲労（TMF: Thermo-Mechanical Fatigue）が支配的な破壊モードになる。

はんだ接合の疲労寿命は Coffin-Manson 則で評価する：

$\Delta \gamma$ はせん断ひずみ範囲、$C$ と $n$ は材料定数。Sn-Ag-Cu系はんだなら $n \approx 1.9$, $C \approx 50$ 程度。つまりひずみ範囲を半分にすると寿命が約4倍になる。CTE不整合を減らすアンダーフィル（樹脂充填）の効果はまさにこの $\Delta \gamma$ の低減だ。

主要ツールの熱-構造連成機能を比較しよう：

機能面では大差ないから、用途と既存環境で決めるのが実務的だ：

それと「チームに使える人がいるか」も重要。高機能なソフトでも使いこなせなければ意味がない。まずは順次連成（片方向）で感度をつかんで、必要に応じて強連成に切り替える——この段階的アプローチが現場では最も効率的だ。

入る。むしろ熱-構造連成トポロジー最適化はホットトピックだ。目的関数の典型例：

$C_{mech}$ は構造コンプライアンス（剛性の逆数）、$C_{th}$ は熱コンプライアンス（温度分布の均一性指標）。$w_1$, $w_2$ は重み。例えばヒートシンクの設計で「熱抵抗を最小化しつつ、取り付けボルト荷重に耐える剛性を確保する」という多目的最適化を行える。

実際に出てくる。有機的な樹状構造になることが多い。ただし製造制約（抜き勾配、最小板厚）を入れないと、金型で作れない形状になるから要注意。3Dプリンタとの組み合わせが本命だね。

最先端の例を紹介しよう。半導体露光装置では、実機のセンサー（温度、変位）データをリアルタイムで縮退モデル（ROM: Reduced Order Model）に入力し、構造全体の熱変形を推定して補正する。FEMで作った精密モデルを数千の固有モード和に縮退させて、計算時間をms級に短縮しているんだ。

そう。「設計段階のCAE」が「運用段階の制御」とつながる——これがデジタルツインの本質だ。ガスタービンでもGE やSiemens Energyが同様のアプローチで運転中のブレード温度と寿命を推定している。

むしろAMこそ熱-構造連成設計の最前線だ。レーザー粉末床溶融（LPBF）では、レーザーが溶融池を作りながら移動する過程で：

これを予測するAMプロセスシミュレーション（Ansys Additive, Simufact Additive, Amphyon等）が急速に発展している。数千層の積層を固有ひずみ法で効率的に計算し、ベースプレート取り外し後の変形を事前に予測する。トポロジー最適化 → AMプロセスシミュレーション → 後処理設計という一気通貫ワークフローが現実的になりつつある。

構造解析のトラブルに加えて、連成特有のものがある。よくある5つを紹介しよう：

まさにそう。Abaqusのデフォルト $T_{ref}$ は0°C、Ansysはモデルに依存する。「自由膨張ケース」（全拘束を外して温度だけ与える）で応力がゼロになることを最初に確認する——これが熱-構造連成解析のスモークテストだ。

結果の品質を保証するために、以下のチェックを毎回行うこと：

うん。熱-構造連成設計は「温度を構造に渡すだけ」じゃなく、材料選定・冷却設計・拘束設計・疲労評価まで含むトータルな設計方法論だ。FEMはその中の強力なツールだけど、最終的には「何のために計算するか」——寸法精度なのか、寿命なのか、安全率なのか——を常に意識して進めてほしい。