フーリエの法則 — トラブルシューティング
より充実した内容を fourier-law.html でご覧いただけます。
よくあるトラブルと対策
定常熱伝導解析で困ったとき、まず何を見ればいいですか？
よくある問題を整理しよう。
1. 温度が発散する・非現実的な値になる
原因: 境界条件の不足。定常熱伝導には最低1つの温度拘束（Dirichlet条件）または対流条件（Robin条件）が必要だ。全面にNeumann条件（熱流束）だけ指定すると、温度の基準が決まらず解が一意に定まらない。
対策: モデルのどこかに温度の基準点を設ける。物理的に温度が既知の面を見直す。
それは構造解析でいう剛体変位の問題と同じですね。
まさにそうだ。構造で拘束なしだと剛体移動するのと同様に、熱で温度拘束なしだと温度が浮動する。
2. エネルギー収支が合わない
原因: メッシュ不整合、接触面のギャップ、材料定義の誤り。
対策: 反力の合計値（Ansysなら FSUM, Abaqusなら Output > Energy）を確認。流入熱量と流出熱量の差が全体の1%未満であることを検証する。大きなずれがある場合、ノードのマージ漏れやTied Contactの設定ミスを疑う。
3. 温度分布が予想と大きく異なる
解析は回ったけど結果が実測と合わないときは？
チェックリストはこうだ。
| チェック項目 | よくあるミス |
|---|---|
| 単位系 | W/(m K) と W/(mm K) の混同（1000倍の差） |
| 熱伝達係数h | 自然対流で5〜10 W/(m2K) のところに100を入れている |
| 接触熱抵抗 | TIMやボルト締結部のギャップコンダクタンス未設定 |
| 材料テーブル | 温度依存テーブルの範囲外で外挿されている |
| 放射 | 高温部品で放射を無視している |
特に単位系の混同は致命的だ。Ansys MechanicalはSI単位系（m, kg, s）とmm単位系（mm, t, s）でkの値が変わる。mm系ではk=200がmW/(mm K)=0.2 W/(mm K)になる点に注意。
4. メッシュ収束しない
定常熱伝導で温度のメッシュ収束が悪い場合、特異点の近傍（角部、材料不連続点）を疑う。熱流束は温度勾配の微分なので、温度が収束していても熱流束は収束しにくい。
実務的にはどうすればいいですか？
温度そのものに着目すれば十分なケースが多い。熱流束の局所値が必要な場合は、パスに沿った分布を確認し、特異点から離れた位置で評価する。
異方性材料でのフーリエ則適用
炭素繊維強化プラスチック（CFRP）は繊維方向k≈7 W/m·K、垂直方向k≈0.7 W/m·Kと10倍の異方性を持つ。等方性を仮定して解析すると温度分布の最大誤差が30%に達することがあり、テンソル形式の熱伝導率行列の入力が必須だ。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——フーリエの法則の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告