2. 热力学基础与材料参数
热结构耦合仿真的精度，很大程度上不取决于网格有多密、求解器设置有多复杂，而取决于材料参数是否被正确理解和表达。Abaqus 中的热-力耦合分析需要同时求解温度场和位移场，因此材料定义必须覆盖两类属性：一类决定热量如何传导和储存，另一类决定结构在温度变化下如何变形和受力。本节围绕四个核心参数——热膨胀系数、比热容、导热系数、弹性模量——说明它们的物理含义、在 Abaqus 中的输入方式，以及随温度变化时应当如何处理。
四个参数各自控制什么
热膨胀系数（CTE，符号 α）描述材料在温度变化时产生热应变的能力。在 Abaqus 中通过 *EXPANSION 定义，单位是 1/K。它直接决定热应变的来源：当温度从参考温度 T₀ 变化到 T 时，自由热应变为 ε_th = α(T)(T − T₀)。这里有一个必须明确的约定：Abaqus 的参考温度在 *EXPANSION 中通过 ZERO 参数指定，默认值为 0，实际建模中通常设为装配或初始温度。若参考温度设置错误，即使 α 数值正确，热应变也会整体偏移。
比热容（c，单位 J/(kg·K)）决定材料储存热量的能力，通过 *SPECIFIC HEAT 定义。它影响瞬态热分析中温度随时间上升或下降的速率，也影响耦合分析中温度场的演化过程。对于稳态热分析，比热容不参与计算；但只要涉及时间相关的温度变化，它就必须给出。
导热系数（k，单位 W/(m·K)）描述热量在材料内部传导的难易程度，通过 *CONDUCTIVITY 定义。它决定温度梯度的大小和热流方向，是温度场分布的主导参数。金属的导热系数通常远高于聚合物和陶瓷，这也是同一热载荷下不同材料温度分布差异明显的原因。
弹性模量（E，单位 MPa 或 Pa）与泊松比一起描述材料的力学刚度，通过 *ELASTIC 定义。在热应力计算中，弹性模量决定同样的热应变被约束时会产生多大的应力：σ = E·α·ΔT（一维约束情形）。因此弹性模量随温度下降而降低时，高温区的热应力会明显小于按常温模量估算的结果。
参数随温度变化：为什么必须用表格而不是常数
很多初学者习惯把材料参数填成单一常数，这在温度变化不大的工况下可以接受，但在热结构耦合中往往带来明显误差。以常见的 304 不锈钢为例，其弹性模量从 20 ℃ 的约 195 GPa 降到 600 ℃ 的约 150 GPa，降幅超过 20%；导热系数从约 16 W/(m·K) 升到约 21 W/(m·K)；热膨胀系数从约 1.6×10⁻⁵ /K 升到约 1.9×10⁻⁵ /K。如果全部按常温常数计算，高温区的热应力和温度分布都会偏离真实情况。
在 Abaqus 中，温度相关性通过在每个材料选项的数据行中增加温度列来实现。例如定义随温度变化的弹性模量：
*ELASTIC
195000., 0.30, 20.
185000., 0.30, 200.
175000., 0.30, 400.
150000., 0.30, 600.每一行依次是弹性模量、泊松比、温度。热膨胀系数、比热容、导热系数的写法相同，只是列数和物理量不同。需要注意两点：第一，Abaqus 会在给定温度点之间做线性插值，超出表格范围时保持端点值不变，因此表格应覆盖分析中可能出现的全部温度区间；第二，热膨胀系数在 Abaqus 中默认是“割线”定义，即相对于参考温度的总膨胀量，而不是某个温度点的瞬时切线值，输入时应确认数据来源的定义方式是否一致。
材料本构模型的选择
热结构耦合中的本构模型选择，取决于分析目标是否允许材料进入非弹性阶段。如果温度变化引起的应力始终低于屈服强度，线弹性模型配合温度相关的弹性模量和热膨胀系数就足够，计算量小、收敛稳定，适合大多数热应力评估和热变形分析。
当热应力超过屈服强度，或分析涉及焊接、热处理、成形等伴随明显塑性变形的过程时，必须引入塑性本构。Abaqus 中通过 *PLASTIC 定义屈服应力和塑性应变随温度的变化，通常与 *ELASTIC 的温度表格配合使用。此时需要额外关注：屈服面随温度移动、硬化曲线是否随温度变化、以及是否考虑蠕变。对于长时间高温服役的部件，蠕变往往不可忽略，可通过 *CREEP 定义。
选择本构模型时可以用一个简单判断：先按线弹性做一次热应力估算，若最大等效应力接近或超过材料在该温度下的屈服强度，就应切换到弹塑性模型；若服役温度超过材料熔点的 0.3～0.4 倍且持续时间较长，再考虑蠕变。这个判断不需要精确，但能避免用错模型导致结果失真。
一个可复用的材料参数检查清单
在提交热结构耦合分析前，建议按以下清单逐项核对材料定义，这些字段直接对应 Abaqus 输入文件中的关键字：
| 检查项 | Abaqus 关键字 | 常见问题 |
|---|---|---|
| 参考温度 ZERO 是否与分析初始温度一致 | *EXPANSION, ZERO | 默认 0 导致热应变整体偏移 |
| 热膨胀系数是割线值还是切线值 | *EXPANSION | 与数据来源定义不一致 |
| 温度表格是否覆盖全部工况温度 | 各材料选项 | 超出范围后取端点值，产生偏差 |
| 比热容与密度是否同时给出 | *SPECIFIC HEAT, *DENSITY | 瞬态分析缺少密度无法计算 |
| 弹性模量与泊松比的温度列是否对齐 | *ELASTIC | 温度点不一致导致插值错误 |
| 是否需要在高温段启用塑性或蠕变 | *PLASTIC, *CREEP | 线弹性高估热应力 |
把这张清单当作提交前的固定动作，可以过滤掉大部分因材料参数设置不当引起的异常结果。热结构耦合的难点往往不在求解器，而在于是否清楚地知道每个参数在方程中扮演什么角色、随温度如何变化、以及当前分析是否允许忽略非弹性效应。理解这四点，后续的边界条件、载荷步和结果判读才有可靠的基础。