熱機械疲労（TMF）解析
理论与物理
什么是TMF
老师，TMF和普通疲劳有什么区别？“热”和“机械”都包含在里面我大概能理解，但是…
问得好。通常的疲劳试验——例如拉伸压缩的反复——是在室温恒定下进行的吧。这被称为等温疲劳。另一方面，在TMF中，温度循环和机械载荷同时作用。想象一下燃气涡轮叶片。
是喷气发动机里的那个叶片吧？会变得超高温的那个。
对。启动时从室温急剧加热到1000〜1200°C，同时因旋转产生的离心力作用承受拉伸载荷。停止时则冷却下来，载荷也卸除。这个循环会重复数千次。问题在于，温度变化时，材料的屈服应力、弹性模量、蠕变特性都会改变。如果仅用室温疲劳数据来预测寿命，可能会与现实相差一个数量级以上。
诶，会差一个数量级吗？那确实很糟糕啊…
所以TMF才需要一个独立的评价体系。列举一些TMF成为问题的典型部件：
- 燃气涡轮叶片/静叶：启停循环下0〜1200°C
- 汽车发动机排气歧管：冷启动→高温稳态→冷却，200〜900°C
- 发电用蒸汽涡轮转子：启动时的热瞬态导致表面与内部温差
- 刹车盘：制动时的摩擦加热与冷却的反复
- 核电管道热分层：不同温度流体混合引起的温度波动
同相与反相
看TMF的教科书时会出现“IP”和“OP”，这些是什么？
这是TMF最重要的概念。根据温度与机械应变的相位关系，破坏模式会发生根本性变化。
- 同相：拉伸应变最大时温度也最高。高温与拉伸的组合下，蠕变损伤占主导。裂纹多沿晶界从内部扩展。
- 反相：压缩应变最大时温度最高。高温下形成氧化层，冷却时转为拉伸则氧化层开裂，产生表面裂纹。氧化损伤占主导。
即使是相同的材料、相同的应变范围，仅仅因为相位不同，裂纹形态就会改变吗？
没错。例如镍基超合金IN738LC，在相同的机械应变范围0.6%下，IP-TMF的寿命可能是OP-TMF的2〜5倍。IP条件下，蠕变空穴在晶界聚集导致缓慢破坏；而OP条件下，氧化-裂纹开口的循环进展迅速。在实际部件中，相位不一定正好是0°或180°，也常取中间角度。这被称为菱形循环。
实际的燃气涡轮叶片，哪个部位是IP，哪个部位是OP呢？
很好的着眼点。叶片的前缘在启动时直接暴露于高温燃气，温度先上升，而叶片内部仍较冷。因此前缘容易形成压缩应变+高温→OP-TMF条件。另一方面，轮盘的内径部，离心力引起的拉伸应变与稳态运行时的温度上升同步，因此形成IP-TMF条件。也就是说，在同一发动机内部，IP和OP是共存的。
损伤机制的叠加
普通疲劳就是反复载荷导致裂纹产生，这么简单，但TMF就更复杂了是吧？
TMF中至少有3种损伤机制同时进行：
- 纯疲劳损伤：机械应变反复引起的位错累积→裂纹形核。与低周疲劳机制相同，但材料的循环响应随温度变化。
- 氧化损伤：高温保持期间表面氧化层生长，冷却时的应变失配导致氧化层开裂。镍基合金中，Al₂O₃或Cr₂O₃保护膜破坏后氧化会加速。
- 蠕变损伤：高温、拉伸应力下晶界滑移和空穴扩散进行。长时间保持导致损伤累积，直至晶界开裂。
这三者并非线性叠加，而是相互作用。如果氧化脆化的晶界处聚集了蠕变空穴，破坏会比单独作用快得多。
控制方程与应变分解
TMF的寿命预测用什么方程？
TMF的一切都始于应变的分解。从总应变 $\varepsilon_{total}$ 中减去热应变 $\varepsilon_{th}$，提取出直接关系到寿命的机械应变：
这里 $\alpha(T)$ 是温度依赖的线膨胀系数。机械应变进一步分为弹性成分和非弹性成分：
分成弹性、塑性、蠕变三种啊。那寿命怎么求呢？
最基本的方法是修正Manson-Coffin公式。这是将等温疲劳的Manson-Coffin公式改为温度依赖的形式：
各参数的含义如下：
- $\sigma'_f(T)$：温度依赖的疲劳强度系数
- $b(T)$：疲劳强度指数（通常 $-0.05$ 〜 $-0.12$）
- $\varepsilon'_f(T)$：疲劳延性系数
- $c(T)$：疲劳延性指数（通常 $-0.5$ 〜 $-0.7$）
- $E(T)$：温度依赖的杨氏模量
但是，仅凭这个公式无法考虑氧化和蠕变的影响。这是修正Manson-Coffin公式最大的局限。
Neu-Sehitoglu模型
有没有包含氧化和蠕变的、更精密的模型？
被视为TMF寿命预测黄金标准的是Neu-Sehitoglu模型。由伊利诺伊大学的Neu和Sehitoglu于1989年提出。它将损伤分为3种机制，分别独立计算后相加：
我们来看看每一项。
1. 疲劳项 $N_f^{fat}$：就是修正Manson-Coffin公式本身。通常使用最高温度 $T_{max}$ 下的材料参数。
2. 氧化项 $N_f^{ox}$：考虑氧化层的生长速度和裂纹开口的时机。氧化层厚度 $h_{ox}$ 按抛物线规律生长：
这里 $Q_{ox}$ 是氧化的活化能，$D_0$ 是频率因子，$R$ 是气体常数。氧化损伤在OP条件下尤其显著——因为高温下氧化层生长，冷却时拉伸导致其开裂。
3. 蠕变项 $N_f^{cr}$：对高温保持期间的蠕变应变累积进行积分。假设Norton定律：
蠕变损伤在IP条件下变得严重。因为拉伸应变最大时温度高，蠕变变形加速。这样，Neu-Sehitoglu模型就能自然地再现IP条件下蠕变项占主导，OP条件下氧化项占主导的实验事实。这是单一公式（修正Manson-Coffin）无法做到的。
原来如此…因为把损伤分开来加，所以IP和OP下自动起主导作用的项就不同了。真聪明！
各项的物理意义与单位制
| 参数 | SI单位 | 典型值（镍基超合金） |
|---|---|---|
| $\sigma'_f$ 疲劳强度系数 | MPa | 1200〜2000 MPa |
| $b$ 疲劳强度指数 | 无量纲 | $-0.06$ 〜 $-0.10$ |
| $\varepsilon'_f$ 疲劳延性系数 | 无量纲 | $0.1$ 〜 $1.0$ |
| $c$ 疲劳延性指数 | 无量纲 | $-0.5$ 〜 $-0.8$ |
| $Q_{ox}$ 氧化的活化能 | J/mol | 150〜250 kJ/mol |
| $Q_{cr}$ 蠕变的活化能 | J/mol | 250〜400 kJ/mol |
| $n$ Norton定律的应力指数 | 无量纲 | $3$ 〜 $8$ |
| $\alpha$ 线膨胀系数 | 1/K | $12$ 〜 $16 \times 10^{-6}$ |
与SRP法（应变范围分割法）的区别
SRP法将非弹性应变分割为 PP（塑性-塑性）、CC（蠕变-蠕变）、PC、CP 四种，对每种分别应用独立的Manson-Coffin关系。由NASA Glenn（原Lewis）研究所开发，在航空航天领域广泛应用。Neu-Sehitoglu模型按损伤机制（疲劳/氧化/蠕变）分离，而SRP法按应变产生过程分离——方法不同，但两者都能处理IP/OP的差异。SRP法的弱点是需要大量材料常数，试验成本高。
TMF试验机“一台数亿日元”的理由
TMF的实验评价需要专用试验机，其价格通常在1〜3亿日元一台。普通疲劳试验在室温下进行，而TMF试验需要同时用高频感应加热对试样进行快速加热，并实时同步施加精密的拉伸/压缩载荷。温度控制精度在±2°C以内，应变控制用高温引伸计达到±0.001%——这些都需要在数百至数千个循环中稳定重复的技术。模拟航空发动机涡轮叶片的OP-TMF试验，其寿命可能仅为等温疲劳的1/5〜1/10。正是为了用仿真预测来替代这种高成本的试验数据，TMF分析才具有其存在意义。
数值解法与实现
温度依赖本构关系
用FEM做TMF分析时，材料模型怎么处理？温度一变，所有参数都变了对吧？
可以说，TMF分析的成败几乎取决于本构关系（材料模型）的精度。至少需要：
- 温度依赖的弹塑性模型：屈服应力、硬化法则（各向同性硬化、随动硬化）的参数需要在多个温度下标定。Chaboche模型（非线性随动硬化）实际上是TMF的标准。
- 蠕变法则：Norton定律或Norton-Bailey时间硬化法则。用于再现高温保持的影响。
- 统一本构模型：Bodner-Partom、Walker、Chaboche-Lemaitre等。不区分塑性和蠕变，作为粘塑性应变统一处理。对于TMF，统一本构模型是理想的，但参数标定非常困难。
每个温度的参数都不同，那数据库不会变得非常庞大吗？
是的。例如Chaboche模型有8个参数，如果在6个温度水平下进行标定，就需要48个数据。而且需要在每个温度下进行循环应力-应变试验。航空发动机制造商花费数十年时间构建了专用的材料数据库。这可以说是TMF分析的最大进入壁垒。
FEM热-结构耦合
用FEM同时计算温度和应力对吧？耦合方法不止一种吗？
主要有两种方法：
1. 顺序耦合（弱耦合）：先通过热分析求出温度场的时程，然后将其作为温度载荷映射到结构分析中。对于大多数TMF情况，这已经足够。因为结构变形对温度场的反馈效应通常很小。
2. 完全耦合（强耦合）：同时求解热和结构。当塑性功发热（Taylor-Quinney系数）影响温度场，或接触面的热传导依赖于变形状态时需要。刹车盘的TMF就是例子。
更具体地说，顺序耦合的分析流程是：
- 通过热传导分析计算启动→稳态→停止的温度场时程（$T(\mathbf{x}, t)$）
- 将温度场传递到结构网格
- 在每个时间步施加 $\Delta\varepsilon_{th} = \alpha(T)\Delta T$ 作为温度载荷
- 执行非线性结构分析（温度依赖弹塑性+蠕变）
- 提取稳定循环的应力-应变-温度历程 → 输入到寿命预测模型
疲劳后处理
从FEM结果得出寿命，是在FEM求解器内部完成的吗？
在实际应用中，通常使用专用的疲劳后处理软件。读取FEM求解器输出的应力、应变、温度时程数据，在每个节点计算TMF寿命并显示云图。代表性的工具：
- fe-safe（达索系统）：与Abaqus紧密集成。有TMF专用模块。支持Neu-Sehitoglu模型。
- nCode DesignLife（HBK/西门子）：支持多求解器。可选择SWT、SRP、Manson-Coffin等。
- FEMFAT（ECS）：在汽车OEM中广泛采用。处理速度快。
为什么不在FEM内部完成？因为疲劳模型需要根据材料、环境、表面处理频繁更新，专用软件在数据库管理和方法切换方面更…
なった
詳しく
報告