熱サイクル試験 — CAE用語解説

その通りです。定常状態を仮定する定常熱伝導とは異なり、熱サイクルは時間変化が本質です。支配方程式は非定常（過渡）熱伝導方程式で、以下のように表されます。

はい、境界条件または発熱量Qを時間関数として定義します。例えば、自動車のパワーエレクトロニクスモジュールの試験規格「AQG 324」では、-40℃から+125℃の間で温度を繰り返し変化させます。シミュレーションでは、モデル表面の熱伝達率を非常に大きな値（例: 10000 W/m²K）に設定し、周囲温度を時間関数

良い着眼点です。無視できない場合がほとんどです。例えば、はんだのヤング率は室温で約40 GPaですが、125℃では約20 GPaと半減します。熱伝導率や比熱も変化します。高精度なシミュレーションでは、材料特性を温度依存性としてテーブル入力する必要があります。Ansysの材料ライブラリには、サーマルサイクリング解析用に、SnAgCuはんだなどの温度依存データが既に登録されているものもあります。

温度そのものよりも、温度変化に伴って生じる「熱応力」とそれによる「疲労破壊」が評価対象です。異なる材料が接合されている部分、例えば半導体チップと基板を接合するはんだ部や、ワイヤーボンディング部で、熱膨張係数（CTE）の差により応力が発生します。CAEでは、熱解析で得た温度分布を構造解析に読み込み、熱応力を計算します。その繰り返しによる疲労寿命を、Coffin-Manson則などのモデルを用いて予測します。

時間領域を離散化する「時間積分法」を用います。代表的なのは陰解法の一種である「後退差分法（Backward Euler）」や「クランク・ニコルソン法」です。これらは無条件安定なので、熱サイクル解析のように長時間の計算に向いています。例えば、時刻

それは状況によりますが、一般にダメです。熱がモデル内を伝わる速度（熱拡散率）と要素サイズを考慮する必要があります。粗いメッシュで大きな

熱伝導のみを考えるのであれば、可能です。1サイクル終了時の温度分布を、次のサイクルの初期条件として設定します。しかし、熱応力と塑性変形を考慮する「熱-構造連成解析」では話が変わります。1サイクル目で生じた塑性ひずみは2サイクル目に引き継がれ、累積していきます。このため、寿命予測では複数サイクルを連続して計算するか、安定した応力-ひずみループが形成されるまでの数サイクルを計算し、その結果を基に寿命を外挿する手法が取られます。

まずは「試験規格」の確認です。対象製品や業界によって適用される規格が異なります。自動車電子部品なら「AQG 324」や「JESD22-A104」、民生品なら「JIS C 60068-2-14」などです。規格にサイクル条件（温度範囲、昇降温レート、保持時間）、許容差が定義されているので、それをシミュレーションの入力条件とします。例えば、AQG 324の温度プロファイルは、実際の試験装置の制御ロジックを反映して、単純な三角波ではなく台形波になることが多いです。

熱応力・ひずみが集中し、破壊の起点となる「インターフェース」部分です。具体的には、半導体チップと基板間のはんだ接合部、ワイヤーボンディングの足元、モールド樹脂とリードフレームの界面などです。これらの層は物理的に薄い（はんだは数十～百µm）ため、少なくとも3層以上の要素厚さでメッシュを切る必要があります。Abaqus/CAEでは「ボルトメッシュ」機能を使って、これらの薄層領域に構造化メッシュを生成するのが一般的です。

実験データがなくても行える検証手順があります。1) **収束性検証**: メッシュサイズと時間ステップ幅を段階的に細かくし、関心量（最大応力、特定点の温度履歴）の変化が無視できるレベルになることを確認する。2) **エネルギー収支チェック**: ソルバーの出力ログで、系への流入熱量と内部エネルギー増加・流出熱量のバランスが取れているか確認する。3) **極限ケースでの比較**: 熱伝達率を無限大に設定した場合、表面温度が周囲温度に瞬時に追従するか、あるいは断熱条件で温度が一定になるかなど、理論的に予測できる単純ケースでソルバーの挙動をテストする。

Coffin-Manson則は

もちろんAbaqusでも可能です。手法のアプローチが少し異なります。Ansysでは「Transient Thermal」で温度分布を求め、その結果を「Static Structural」に読み込んで熱応力を計算する、いわゆる「順次連成解析」が標準的です。一方、Abaqusでは「Coupled temperature-displacement」解析ステップを利用して、熱伝導方程式と力の平衡方程式を完全連成で解くこともできます。特に、接触条件が温度に依存する場合や、発熱がひずみ速度に依存するような問題では、Abaqusの連成解析が有利です。

COMSOLは、熱サイクルに伴う複雑な物理を直接的にモデル化するのに適しています。例えば、はんだ接合部の「クリープ」挙動は温度と応力に強く依存し、時間とともに緩和します。COMSOLの「Nonlinear Structural Materials Module」を使えば、温度依存のアネーラルド・パワー則クリープモデルなどを熱-構造連成解析に直接組み込めます。また、ユーザーが支配方程式をカスタマイズできるため、学術研究や新しい疲労モデルの開発プロトタイプとしてよく用いられます。

はい、その代表格が「Ansys Sherlock」です。これは汎用FEMソルバーではなく、PCB（プリント基板）の設計データ（ODB++やIPC-2581形式）を直接読み込み、部品配置、配線パターン、ビア情報から自動的に簡易FEMモデルを生成します。材料ライブラリには電子部品用の材料特性が豊富で、JEDECやIPCの規格に基づく熱サイクル条件を適用し、はんだ接合部の疲労寿命を迅速に見積もることができます。設計の初期段階で多くのケーススタディを行う「スクリーニング」用途に強みがあります。

理論的には可能ですが、実務的なハードルは高いです。Code_Asterでは熱-構造連成解析が可能ですが、温度依存の材料特性や複雑なクリープモデルの入力、後処理での疲労寿命計算の自動化などには、独自のスクリプト（Python）をかなり書く必要があります。また、商用ソフトのような直感的なGUIや、電子部品向けに最適化された材料ライブラリ、検証済みの疲労モデルは基本的にありません。研究目的や、ソルバーそのもののアルゴリズムを学ぶには良いですが、製品開発のスケジュール内で信頼性の高い結果を出すには、専門知識と時間が要求されます。

最も多い原因は、「熱膨張の基準温度（Reference Temperature）」の設定ミスです。熱応力は、ある基準温度からの温度変化

材料特性の温度依存性が急激に変化するポイントでよく起こる問題です。例えば、樹脂材料の比熱や熱伝導率のデータが、ある温度を境に不連続に変化するテーブル入力になっていると、ソルバーが収束に苦労します。対策としては、1) 材料データをスムージングする、2) 時間ステップ幅を一時的にさらに小さくする（ソルバーの自動制御で可能）、3) 発散する温度付近で材料特性を一定値に近似する、などがあります。また、相転移や潜熱を考慮している場合も、同様の現象が起きるので注意が必要です。

これは薄層要素に生じる「ロッキング」と呼ばれる数値的な問題です。特に完全積分要素を使っていると、せん断変形を過度に剛に評価してしまい、異常に高い応力が計算されます。対策は以下の通りです。1) **要素タイプの変更**: 減積分要素（AbaqusのC3D8Rなど）を使用する。2) **要素形状の改善**: 薄層部分の要素アスペクト比をできるだけ1に近づける（現実的には難しいが、10以下を目指す）。3) **特殊要素の使用**: インターフェース専用のコヒーシブ要素やコンタクト要素を用いて薄層をモデル化する。Ansysでは「Interior Bonded」接触設定で代用できる場合もあります。

寿命が桁で違う場合は、根本的なモデル化か材料定数に問題がある可能性が高いです。チェックすべき点は: 1) **使用した疲労モデルの適用限界**: Coffin-Manson則は塑性ひずみが支配的な場合のモデルです。弾性ひずみが主の場合は、応力ベースのモデル（Basquin則）を検討する。2) **平均応力効果の考慮**: 同じひずみ幅でも、引張平均応力がある