2. 热分析基础与材料参数

热结构耦合分析的可靠性，首先取决于热分析本身是否被正确建立。热分析要回答的核心问题只有一个：温度场如何分布，以及它随时间怎样变化。而温度场一旦确定，热应变、热应力、热变形才有计算的基础。因此本节先厘清三种基本传热方式，再说明四类关键材料参数随温度变化的定义方法。

一、三种基本传热方式

热传导发生在固体内部或相互接触的物体之间，依靠微观粒子的热运动传递能量，不需要物质宏观流动。其数学描述为傅里叶定律：热流密度与温度梯度成正比，比例系数就是导热系数 k（W/(m·K)）。一维稳态形式为 q = -k·dT/dx，负号表示热量自发从高温流向低温。在 ANSYS 中，热传导由单元自身的导热矩阵体现，只要赋予材料导热系数，实体单元即可自动参与传导计算。

热对流发生在固体表面与流体之间，依靠流体的宏观运动带走或带来热量。牛顿冷却定律给出 q = h·(T_s − T_f)，其中 h 为对流换热系数（W/(m²·K)），T_s 为壁面温度，T_f 为流体参考温度。对流是典型的边界条件，在 ANSYS 中通过表面效应单元或直接施加 Convection 载荷实现。h 的取值对结果影响极大：自然空气对流约 5～25，强制风冷可达 25～250，水冷则常在 500～10000 量级，工程中应优先依据经验关联式或实测确定，而不是随意取值。

热辐射通过电磁波传递能量，不需要介质，且与绝对温度的四次方成正比。净辐射热流可写为 q = εσ(T_s⁴ − T_sur⁴)，其中 ε 为发射率（0～1），σ 为斯特藩-玻尔兹曼常数 5.67×10⁻⁸ W/(m²·K⁴)。辐射在高温工况（如数百摄氏度以上）不可忽略，且具有强烈的非线性，会显著增加求解迭代次数。ANSYS 中通过 SURF152 等表面单元或 Radiosity 求解器处理，需要定义发射率和视角因子。

二、材料参数随温度变化的定义方法

热结构耦合分析中，材料属性往往不是常数。金属的导热系数、比热容、弹性模量、热膨胀系数都会随温度明显变化。若全部按室温常数处理，高温工况下的热应力和变形可能产生百分之十几甚至更大的偏差。因此需要掌握在 ANSYS 中定义温度相关属性的方法。

在 Engineering Data 中，每项材料属性都可以从 Constant 切换为 Tabular，输入“温度—属性值”数据对。例如某不锈钢的弹性模量可定义为：20 ℃ 时 195 GPa，200 ℃ 时 185 GPa，400 ℃ 时 175 GPa，600 ℃ 时 165 GPa。求解器会在相邻数据点之间线性插值，超出范围则按端点值外推。因此表格的温度范围必须覆盖实际工况温度区间，否则外推会带来不可控误差。

四类参数的作用各不相同：导热系数 k 决定温度场本身，是热分析的核心输入；比热容 c 与密度 ρ 共同决定热惯性，瞬态分析中 ρc 越大，温度响应越迟缓；热膨胀系数 α 把温度变化转换为热应变 ε_th = α·ΔT，是热应力的直接来源；弹性模量 E 决定材料对热应变的约束刚度，E 越高，同样约束条件下的热应力越大。四者中任一随温度变化，都会改变最终的热应力结果。

需要特别注意热膨胀系数的参考温度。ANSYS 中热应变按 α·(T − T_ref) 计算，T_ref 默认常为 0 ℃ 或用户指定值。若参考温度设置不当，即使温度场正确，热应变也会整体偏移。工程上通常把参考温度设为装配或应力自由状态的温度，并在全文和模型中保持一致。

上图的逻辑是：三种传热方式共同决定温度场，而温度场又通过随温度变化的材料参数反馈到热应变与热应力计算中。理解这条链路，就能明白为什么材料参数表必须覆盖工况温度范围。

三、实操要点与常见问题

在 Workbench 中建立温度相关材料时，建议先在 Engineering Data 中完成表格录入，再回到 Mechanical 中检查材料分配是否正确。对于瞬态热分析，密度和比热容缺一不可；对于稳态热分析，比热容和密度不参与计算，但若后续要做热应力，仍需保证弹性模量和热膨胀系数完整。

常见问题有三类。第一，参数表温度范围不足，导致高温区外推失真，解决方法是把表格上限设到工况最高温度以上 10%～20%。第二，热膨胀系数参考温度与装配温度不一致，造成整体热应变偏移，应在材料或分析设置中显式指定 T_ref。第三，忽略辐射或对流系数取值随意，使温度场偏差过大，建议先用简化手算估算各传热路径的热阻，再决定边界条件的详细程度。

掌握以上内容后，就可以进入几何建模与网格划分环节。热分析的质量，很大程度上在材料参数和边界条件阶段就已经决定。