熱衝撃解析
理论与物理
热冲击
老师，什么是热冲击？
急剧的温度变化所产生的热应力。陶瓷、玻璃、耐火材料是典型。表面急冷时表面产生拉伸应力→裂纹。
热冲击抗力
Kingsley的R因子: $R = \sigma_f(1-\nu)/(E\alpha)$。R越大，耐热冲击性越高。
总结
热冲击的断裂力学：热冲击系数R
热冲击（Thermal Shock）是因急剧温度变化产生的瞬时热应力，是导致陶瓷产生裂纹的主要破坏机制。Hasselman（1969年，GE研发中心）将热冲击抗力定义为R＝σf(1−ν)/(Eα)（第一热冲击抗力，裂纹发生临界温差）。高韧性陶瓷ZrO₂（PSZ）的R约为Al₂O₃的3倍，这是其被用于隔热材料/TBC涂层的原因之一。玻璃（钠钙玻璃）的R约为80°C，急冷实验证明，在100°C水中落下时其破坏概率为50%。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您是否有过急刹车时身体被向前甩出的经历？那种“被带走的感觉”正是惯性力。物体越重越难启动，一旦启动也越难停止。地震时建筑物摇晃，也是因为地面突然移动而建筑物的质量“被落下”。静力分析中此项设为零，这是“因为缓慢施力所以加速度可以忽略”的假设。在冲击载荷或振动问题中绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 和 $\nabla \cdot \sigma$。拉伸弹簧时能感觉到“想要恢复的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——铁棒和橡皮筋，用相同的力拉伸，哪个伸得更长？当然是橡皮筋。这种“难以拉伸的程度”就是杨氏模量 $E$，它决定了刚度。常见的误解：“刚度高＝强度高”是不对的。刚度是“不易变形的程度”，强度是“不易破坏的程度”，是不同的概念。
- 外力项（载荷项）：体积力 $f_b$（重力等）和表面力 $f_s$（压力、接触力等）。可以这样理解——桥上卡车的重量是“作用在整个内部上的力”（体积力），轮胎压路面的力是“只作用在表面上的力”（表面力）。风压、水压、螺栓紧固力…全都是外力。这里容易犯的错误：弄错载荷方向。本想“拉伸”却变成了“压缩”——听起来像笑话，但在3D空间中坐标系发生旋转时确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。试着弹一下吉他的弦。声音会一直持续吗？不，会逐渐变小。这是因为振动能量通过空气阻力或弦的内部摩擦变成了热。汽车的减震器也是同样原理——特意吸收振动能量来改善乘坐舒适性。如果阻尼为零会怎样？建筑物在地震后会一直摇晃不停。实际上不会这样，所以设置适当的阻尼很重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 输入为mm时，载荷·弹性模量也需统一为MPa/N系 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢: 约210 GPa，铝: 约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm系中为tonne/mm³（钢为 = 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm系为N，m系统一为N |
数值解法与实现
热冲击的FEM
1. 非稳态热传导分析 — 随时间变化的温度分布
2. 各时刻的结构分析 — 温度分布→热应力
3. 确定最大应力的时刻 — 温度梯度最大的时候
总结
瞬态热应力分析步骤
热冲击的瞬态分析流程为①设置瞬间冷却或加热的边界条件（表面传热系数h值起主导作用），②用隐式法（Crank-Nicolson法）对非稳态热传导方程进行时间积分，③计算各时刻温度分布产生的热应变+弹性应力。陶瓷由于导热率低，表面与内部会产生数百°C的温差，导致表面产生拉伸·压缩的急剧应力变化。Ansys Transient Thermal→Static Structural的耦合分析可自动化此过程，推荐时间步长小于热扩散时间（L²/α）的1/10。
线性单元（1次单元）
节点间线性插值。计算成本低，但应力精度低。注意剪切锁定（可通过减缩积分或B-bar法缓解）。
二次单元（带中间节点）
可表现曲线变形。应力精度大幅提高，但自由度约增加2~3倍。推荐：应力评估重要时使用。
完全积分 vs 减缩积分
完全积分：有过约束（锁定）风险。减缩积分：有沙漏模式（零能量模式）风险。根据情况选择。
自适应网格
基于误差指标（ZZ估计量等）的自动细化。高效提高应力集中部位的精度。有h法（单元细分）和p法（增加阶次）。
牛顿·拉夫森法
非线性分析的标准方法。每次迭代更新切线刚度矩阵。在收敛半径内具有二次收敛性，但计算成本高。
修正牛顿·拉夫森法
切线刚度矩阵使用初始值或每隔数次迭代更新。每次迭代成本低，但收敛速度为线性。
收敛判定标准
力残差范数: $||R|| / ||F_{ext}|| < \epsilon$（通常 $\epsilon = 10^{-3}$〜$10^{-6}$）。位移增量范数: $||\Delta u|| / ||u|| < \epsilon$。能量范数: $\Delta u \cdot R < \epsilon$
载荷增量法
不一次性施加全部载荷，而是分小步增加。弧长法（Riks法）可超越载荷-位移关系的极值点进行追踪。
直接法 vs 迭代法的比喻
直接法是“用笔算精确解联立方程”的方法——可靠但大规模问题耗时过长。迭代法是“反复猜测逼近正确答案”的方法——最初答案粗糙，但每次迭代精度提高。就像查字典时，从第一页开始顺序查找（直接法）不如预估位置翻开再前后调整（迭代法）来得高效。
网格阶次与精度的关系
1次单元是“用直尺近似曲线”——用直线折线表现，精度有限。2次单元是“柔性曲线”——可表现曲线变化，即使网格密度相同精度也显著提高。但每个单元的计算成本增加，需根据总体的成本效益来判断。
实践指南
实务检查清单
反应堆紧急堆芯冷却的热冲击评估
反应堆紧急堆芯冷却系统（ECCS）启动时，低温冷却水（约20°C）被急速注入高温（约320°C）的反应堆压力容器。这300°C温差的熱冲击（Pressurized Thermal Shock, PTS）会在容器壁产生最大400MPa的瞬态拉伸应力。美国NRC的监管指南1.99Rev.2要求考虑辐照脆化的断裂韧性评估（RTNDT转变温度），西屋公司的HEATH分析代码与3D-FEM结合的评价已成为国际标准。在日本，东芝能源系统与解决方案公司实施同等评估。
分析流程的比喻
分析流程其实和烹饪非常相似。首先采购食材（准备CAD模型），进行预处理（网格生成），开火烹饪（求解器执行），最后装盘（后处理可视化）。这里有个重要问题——烹饪中最容易失败的工序是哪里？其实是“预处理”。网格质量差的话，无论使用多优秀的求解器，结果都会一团糟。
初学者容易陷入的陷阱
您确认了网格收敛性吗？是否认为“计算能运行＝结果正确”？这其实是CAE初学者最容易掉入的陷阱。求解器一定会根据给定的网格返回“一个差不多的答案”。但如果网格太粗糙，这个答案就会与现实相差甚远。至少用3个级别的网格密度确认结果是否稳定——如果忽略这一点，就会陷入“因为是计算机给出的答案所以应该正确”的危险误区。
边界条件的思考方式
边界条件的设置，与考试的“出题”是相同的。如果题目出错了呢？无论计算多么精确，答案都是错的。“这个面真的完全固定吗”“这个载荷真的是均匀分布吗”——正确建模现实的约束条件，其实是整个分析中最重要的步骤。
软件比较
工具
所有FEM求解器都支持热-结构耦合。没有差别。
热冲击分析的求解器方法比较
热冲击分析的方法因求解器而异。ABAQUS/Explicit能用显式法追踪接触·断裂，被用于GE制造的燃气轮机隔热涂层（TBC）剥离分析。ANSYS使用ADPCM（自适应热机械耦合）网格追踪裂纹扩展。MSC Nastran自2019年起提供专门针对1,000次循环以上热疲劳的Thermo-Mechanical Fatigue模块。
选定时最重要的3个问题
- “要解决什么问题”：所需的物理模型·单元类型是否支持热冲击分析。例如，流体方面是否有LES支持，结构方面接触·大变形的支持能力会成为差异点。
- “由谁使用”：初学者团队适合GUI完善的工具，有经验者适合脚本驱动的灵活工具。类似于汽车的自动挡（GUI）和手动挡（脚本）的区别。
- “要扩展到什么程度”：基于未来分析规模扩大（HPC支持）、向其他部门展开、与其他工具联动的前瞻性选择，有助于长期降低成本。
尖端技术
尖端
热冲击的断裂力学：陶瓷急冷实验
陶瓷的热冲击破坏在1955年由Kingery实施的氧化铝试片急冷实验实现了量化。从900℃投入水中急冷，当△T=200℃以上时开始开裂，其临界温差可通过毕渥数与断裂韧性KIC的乘积预测。现代太阳能电池板制造中，使用ABAQUS/Explicit分析急冷时的热冲击，实现了KIC＝2.0 MPa√m的Si3N4基板存活率95%的冷却速度设计并已实用化。
故障排除
故障
传热系数（h值）的设置精度问题
热冲击分析中，对结果影响最大的是表面传热系数h值的设置。急冷实验（例：高温陶瓷投入水中）会发生沸腾传热，核沸腾区域h＝10,000〜50,000 W/m²K，膜沸腾区域h＝200〜400 W/m²K，存在100倍以上的差异。工业陶瓷研究所（JFCC，名古屋）2015年的报告书指出，若假设h值为常数，最大热应力会产生50%以上的偏差。结合实验中的IR温度计测数据与逆热传导分析来辨识h值是提高精度的关键。
“分析结果不符”时的应对
- 首先深呼吸——慌张地随机更改设置，问题会变得更加复杂
- 制作最小再现案例——用最简单的形式再现热冲击分析的问题。“减法式调试”最为高效
- 每次只改变一个变量再执行——同时进行多项更改，会不知道哪项起了作用。与科学实验相同的“对照实验”原则
- 回归物理本质——如果计算结果出现“物体逆重力漂浮”等非物理现象，应怀疑输入数据存在根本性错误
なった
詳しく
報告