2次元定常熱伝導 — トラブルシューティング
この記事は統合版に移行しました
より充実した内容を conduction-2d.html でご覧いただけます。
より充実した内容を conduction-2d.html でご覧いただけます。
2次元定常熱伝導 — トラブルシューティング
よくある問題
2次元熱伝導でハマりがちなポイントを教えてください。
1. 平面モデルと軸対称モデルの取り違え
問題: Ansys PLANE55で軸対称を指定すべきところを平面のままにすると、円筒形状の熱抵抗が全く異なる値になる。平面問題は単位厚さあたり、軸対称は2πr を考慮する。
対策: KEY OPT設定を必ず確認。Abaqusでは要素名が異なる（DC2D4 vs DCAX4）ので間違いにくい。
2. 異方性材料の座標系
異方性材料で座標系を間違えるとどうなりますか？
材料座標系とグローバル座標系が不一致だと、面内kと面直kが入れ替わる。PCB基板で $k_{in}=30$, $k_{through}=0.5$ W/(m K) が逆になったら温度上昇が数十倍変わる。Ansys WorkbenchではCoordinate Systemの割り当てを要素ごとに確認すること。
3. コーナー特異点
直角コーナーの内角部では理論的に温度勾配（熱流束）が無限大になる。メッシュを細かくするほど値が増大し、収束しない。
対策: 局所的な最大熱流束値は使わず、ある程度の領域で平均化した値を設計に用いる。または微小なフィレットを付けて特異性を除去する。
構造解析の応力特異点と同じ対処法ですね。
まさに同じだ。温度特異点は実物ではフィレットや面取りで緩和されるので、解析モデルにも微小R（0.1mm程度）を入れるのが現実的だ。
Coffee Break よもやま話
角部の特異性とメッシュ依存性
2次元形状の内角が鋭い再入角では理論的に熱流束が無限大になる特異解が生じる。FEMでこの部分を細かくメッシュ分割するほど温度勾配が大きくなり、収束しない現象が起きる。Abaqusマニュアル（2018年版）では再入角近傍への特殊要素（クォータポイント要素）の適用を推奨している。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——2次元定常熱伝導の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
この記事の評価
ご回答ありがとうございます！
参考に
なった
なった
もっと
詳しく
詳しく
誤りを
報告
報告