📚 课程目录第4章稳态热分析基础

4. 稳态热分析基础

稳态热分析是热结构耦合仿真的起点。它求解的是温度场不再随时间变化的平衡状态：进入物体的热量与离开物体的热量相等，各点温度恒定。Abaqus/Standard 中对应的分析类型为 *HEAT TRANSFER, STEADY STATE，其控制方程由能量守恒与傅里叶导热定律导出。理解方程、边界条件与求解设置，是后续热应力计算能否得到合理结果的前提。

一、稳态热传导方程

三维稳态、无内热源、材料导热各向同性时，温度场 T(x,y,z) 满足拉普拉斯方程：

∂²T/∂x² + ∂²T/∂y² + ∂²T/∂z² = 0

含内热源 q（单位体积生热率，W/m³）时：
∂²T/∂x² + ∂²T/∂y² + ∂²T/∂z² + q/k = 0

其中 k 为导热系数（W/(m·K)）。

若存在内热源，Abaqus 通过 *DFLUX 施加体热流；若为各向异性材料，则用导热张量 k_ij 描述，此时热流方向不再与温度梯度方向一致。稳态问题没有时间项，因此不需要密度与比热，但一旦进入瞬态分析，二者就成为必需参数。

二、四类边界条件

稳态热分析必须有足够的边界条件才能得到唯一解。Abaqus 支持以下四类：

第一类：给定温度（Dirichlet）。用 *BOUNDARY 在节点上指定自由度 11（温度）。例如 NODE, 11, 11, 100. 表示该节点温度固定为 100。它是最强的约束，会直接“钉住”局部温度。
第二类：给定热流（Neumann）。用 *DFLUX 在面或体上施加热流密度（W/m²）或体热源（W/m³）。正号表示热量流入模型，负号表示流出。注意 Abaqus 中面热流默认沿外法线方向为正流出，实际使用时应核对符号约定。
第三类：对流换热。用 *FILM 指定膜层散热系数 h（W/(m²·K)）与环境温度 T_∞，热流 q = h(T − T_∞)。h 依赖流速、几何与流体物性，自然对流通常 5～25，强制风冷可达 25～250。
第四类：辐射换热。用 *RADIATION 指定发射率 ε 与视角系数，遵循斯特藩–玻尔兹曼定律 q = εσ(T⁴ − T_∞⁴)，σ = 5.67×10⁻⁸ W/(m²·K⁴)。辐射在高温（>300 ℃）时不可忽略，且是非线性边界，需要迭代求解。

对流与辐射可以同时施加在同一表面，Abaqus 会将两者热流叠加。若表面同时存在接触，还需考虑接触热导。

三、稳态求解设置

在 Abaqus/CAE 中创建 Heat transfer (Steady state) 分析步，关键设置包括：

求解器：默认使用直接稀疏求解器；若模型含辐射或强非线性，需开启 *STEP, NLGEOM=NO, INC=... 中的自动增量与迭代控制。
收敛容差：温度增量与热流残差默认 1×10⁻⁵，辐射问题可适当放宽到 1×10⁻⁴ 以改善收敛。
初始条件：稳态分析通常不需要初始温度场，但若存在非线性材料或辐射，提供合理初值可显著加速收敛，用 *INITIAL CONDITIONS, TYPE=TEMPERATURE 设置。
单元类型：热传导用 DC3D8（六面体）或 DC3D4（四面体）；若后续要转热应力，应选用带温度自由度的耦合单元如 C3D8T。

一个典型的最小输入片段如下：

*STEP, NAME=STEADY_HEAT
*HEAT TRANSFER, STEADY STATE
1., 1., 1e-5, 1.
*BOUNDARY
FIXED_NODES, 11, 11, 25.
*DFLUX
HEAT_FLUX_FACE, S2, 5000.
*FILM
COOL_FACE, F1, 20., 25.
*RADIATION
HOT_FACE, R1, 0.8, 25.
*END STEP

四、结果后处理与判断

求解完成后，在 Visualization 模块中查看 NT11（节点温度）云图。判断结果是否合理，可从三方面入手：

能量平衡：用 *OUTPUT, HISTORY 输出 HEAT FLUX 与 HEAT ENERGY，检查总输入热流是否等于总输出热流，误差应小于 1%。
温度梯度方向：热流应从高温区指向低温区，若云图出现局部反向梯度，通常是边界条件符号或材料方向定义错误。
网格无关性：将网格加密一倍，若最大温度变化小于 2%，可认为网格足够；否则需继续细化或改用二次单元。

稳态分析中若只施加对流或辐射而没有固定温度节点，刚度矩阵会奇异，求解器报“零主元”错误。务必至少有一个节点温度被约束，或通过 *CONTROLS 启用自动稳定。

把上述四类边界条件与求解设置串起来，稳态热分析的工作流可以概括为：定义材料导热系数 → 划分热传导单元 → 施加至少一个温度约束 → 按需叠加对流、辐射或热流 → 提交稳态求解 → 检查 NT11 云图与能量平衡。只有温度场合理，后续将温度作为预定义场导入热应力分析，才能得到可信的热变形与热应力结果。