热分析收敛失败

好问题。不收敛是指离散化支配方程的残差无法降低到设定的容限值（例如1.0e-6）以下的状态。从物理上讲，这意味着模型内的某处找不到满足热平衡的解，或者计算过程中发散。例如，发热量远远超过散热能力，温度会无限升高，这是尝试计算一个不现实的情况。

非稳态热传导的支配方程源自能量守恒定律。对于体积为V的控制体，根据蓄积项=流入项+生成项，微分形式为：

求解器在每次迭代时监测这个离散化代数方程组的残差范数。在Ansys Mechanical APDL中，用`CNVTOL, HEAT`命令设置这个容限值。

因为非线性性增强了。例如，某种绝缘材料的`k`与温度平方成正比，方程就包含`k(T)`项，变成非线性方程。求解器通过迭代计算来寻找解，但如果初始估计值与实际解相差太远，迭代的解更新量就会发散或振荡。在Abaqus中，对于这类非线性问题，默认进行最多10次迭代（可通过`*CONTROLS, ANALYSIS=DISCONTINUOUS`修改），如果不能在这些迭代内收敛，就判定为失败。

这很重要。在非稳态解析中，在某个时间步长`Δt`内进行"迭代"来求解平衡。如果这个时间步长`Δt`太大，单个步长内的温度变化会变得很急剧，迭代计算就无法收敛。首先应该调整的是"时间增量"，即减小`Δt`。Ansys Fluent的默认时间步长由流动的时间尺度决定，但仅进行热传导时，应该参考材料热扩散率α给出的扩散时间尺度

根本上是不同的。直接法（如MUMPS、PARDISO）用有限次运算"精确地"求解离散化的线性方程组。这里的失败主要由内存不足或奇异矩阵（某些节点的条件未定义等）引起。而迭代法（如共轭梯度法：CG）通过逐次降低残差。这种失败发生在问题的条件数很大（材料特性比例悬殊等）或允许的最大迭代次数太少时。在COMSOL中，定常热分析若使用"分离"求解器（迭代法），前处理器的设置（如几何多重网格）会大大影响收敛性。

它控制的是迭代计算中解的更新量。从迭代i到i+1的温度更新定义为

按以下顺序检查。1. **边界条件缺失**：所有外表面都定义了热边界条件（温度、热流、对流、辐射）吗？未定义的面被视为"绝热"，但有没有无意中被绝热的面？2. **材料特性输入错误**：热导率`k`是否有1.0e-6这样极端的小值，密度`ρ`是否为0？单位制（SI与mm单位制）是否一致？3. **接触热阻的设置**：部件间接触面是否定义了"热接触"？如果没有，热量无法传导，温度会在局部发散。

接下来，质疑一下物理上是否存在解。例如，一个密闭机箱内的CPU发热10W，仅靠机箱外表面的自然对流（传热系数5 W/m²K）散热，此时稳态温度会达到数百℃，材料特性在那个温度下可能会大幅变化。首先，用粗估法计算发热量和散热量。如果发热量`Q_gen`持续大于散热量`Q_loss = h A (T_s - T_amb)`，那么物理上不存在稳态解。需要降低发热量或增强强制冷却（h提高到20 W/m²K以上）。

那个时刻发生了某种"事件"。具体来说：1. **相变（熔化·凝固）**：潜热的吸收·释放使有效热容量变成无穷大，求解器找不到平衡温度。2. **温度相关材料特性的不连续点**：数据表输入间隔太粗，某温度下`k`值突然跳变。3. **接触状态变化**：热接触间隙闭合/张开，系统刚度剧变。对策是在事件发生附近的时间启用自动时间步长控制，并增加最大迭代次数（例如到25次）。

"Analysis Settings"中的以下项目：1. **Solver Type**：从程序控制（通常是直接法）改为"Iterative (PCG)"。对大规模模型有效。2. **Nonlinear Controls**：将"Force Convergence"和"Momentum Convergence"的容限值从默认的0.5%放宽到1%或2%。需要验证结果是否物理合理。3. **Line Search**：开启后会调整更新方向，对强非线性问题有助于收敛。4. **Stabilization**：开启"Energy"稳定，系数从1.0e-3这样的小值开始尝试。这起到数值阻尼的作用。

根本不同。**Standard**使用Newton-Raphson迭代法来求解平衡。收敛失败时，通过`*CONTROLS`参数调整容限残差(`c1`)、最小迭代次数(`i1`)、最大迭代次数(`i2`)。例如`*CONTROLS, PARAMETERS=TIME INCREMENTATION`中设`i1=8, i2=15`增加迭代次数。而**Explicit**使用动态显式法，进行时间积分，本身"收敛"的概念不存在。热传导也是用对角化（lumped）热容量矩阵显式求解。稳定性由稳定时间步长

COMSOL的分离求解器是按物理场顺序逐次求解方程的迭代法（分离法）。收敛失败时：1. **调整求解器序列**：改为"全耦合求解器"。这是求解全部物理场的直接法，更稳定但内存消耗大。2. **调整松弛因子**：使用"高度非线性"设置，或手动将"阻尼系数"降低到0.5以下。3. **改进初值**：使用"辅助扫描"，将某个参数（如发热量）从0逐步增加到最终值，逐步求解。这样每步都能得到好的初值。在包含辐射传热的问题中特别有效。

那是最后手段。如果30次迭代都不收敛，很可能连"朝向解收敛"的趋势都没有。首先看求解器输出（.log或.output文件），检查残差范数的历史。如果单调递减，那增加次数才有意义（Ansys中用`NEQIT`命令）。但如果残差在振荡或停滞，应该怀疑根本原因（边界条件、材料特性、接触），或降低松弛因子。盲目增加迭代次数只会浪费计算时间。

这是严重错误，通常是材料特性输入错误或极端非线性导致数值不稳定。热传导矩阵的对角项应该是正值，代表该节点热流出的难易程度。变成负值原因是：1. **热导率`k`输入了负值**（物理上不可能）。2. **温度相关数据表在外推时`k`变成负数**。3. **辐射的形状因子（View Factor）计算不准确**，产生了实质上的负热阻。首先检查所有材料的`k`在整个温度范围都是正数。在Abaqus中确认`*CONDUCTIVITY`的温度相关数据范围要充分，不要混淆`*EXPANSION`。

会出现。这不仅在显式法中发生，在隐式法的自动时间步长控制中也会发生。求解器为了收敛所需的时间步长`Δt_req`比用户设定的最小时间步长`Δt_min`还要小。原因是：1. **负荷变化急剧**：发热量突然通过阶跃函数开启。对策是将负荷线性升高（用Amplitude定义）。2. **接触状态急变**：间隙闭合的瞬间。适当设置接触刚度。3. **某个单元极小**：热扩散的稳定条件

因为辐射热流与温度的4次方