蠕变-疲劳相互作用 — 高温结构的寿命评估方法
蠕变-疲劳相互作用的理论基础
概要 — 蠕变和疲劳为什么同时成为问题
蠕变和疲劳同时发生会怎样？ 各自单独的情况教科书有，但两者重叠的情况我想象不出来…
简单说，在高温长时间运行的机器——比如核反应堆配管或蒸汽涡轮机——中，定常运行期间蠕变损伤在积累，每次启停时疲劳损伤也在积累。问题是这两种损伤并不独立。
不独立是什么意思？ 就是说"1+1≠2"？
可以说"1+1≥3"。蠕变在晶界上生成孔洞，疲劳裂纹就会从这些孔洞处加速扩展。反过来，疲劳形成的微裂纹也会加快蠕变断裂。所以分别评估的寿命之和往往比实际寿命过于乐观。
具体是哪些部件会有这个问题？
典型的高温材料包括：
- 核反应堆配管、联箱管：550〜600℃运行超过10万小时，每年启停数十次
- 蒸汽涡轮机转子：转子表面的热梯度引起热疲劳+定常运行时的离心力蠕变
- 燃气涡轮动叶片：1000℃以上，离心力+热循环，最苛刻的工况
- 火力发电锅炉高温蒸汽配管：DSS（每日启停）运用导致疲劳循环增加
明白了，都是"高温"+"重复"的组合。那怎样用数值方法定量评估这种相互作用呢？
线性损伤和法则（Dc+Df≦D）
使用最广泛的是线性损伤和法则。在ASME Section III, Subsection NH中规定，用蠕变损伤分率 $D_c$ 和疲劳损伤分率 $D_f$ 的和来评估：
这里 $D$ 是根据材料规定的许用损伤包络线上的值。例如304不锈钢，交点在 $(D_c, D_f) = (0.3, 0.3)$ 附近，不是简单的直线 $D=1$，而是一条双线性的更严格包络线。
$D_c$ 和 $D_f$ 分别怎么求呢？
疲劳损伤分率 $D_f$ 用Miner法则（Miner's rule）计算：
$n_j$ 是第 $j$ 种循环的实际循环数，$N_{d,j}$ 是该工况下的设计疲劳寿命（从S-N曲线读取，再除以安全系数）。
蠕变损伤分率 $D_c$ 用时间分率法求取——下一节详细讲解。
包络线比直线 $D=1$ 更严格是为了什么？ 是不是和刚才说的"1+1=3"有关？
完全对。蠕变与疲劳相互作用时，晶界孔洞和晶内滑动会相乘放大损伤。所以包络线比原点处的直线向内凹。各种材料的凹陷程度不同：
| 材料 | 包络线交点 $(D_c, D_f)$ | 特点 |
|---|---|---|
| 304SS / 316SS | (0.3, 0.3) | 奥氏体钢，晶界薄弱 |
| 2.25Cr-1Mo钢 | (0.1, 0.1) 左右 | 非常严格的包络线 |
| Alloy 800H | (0.3, 0.3) | NH规定值 |
| Modified 9Cr-1Mo (P91) | (0.1, 0.01) | 近期增补，约束很强 |
时间分率法与蠕变损伤
时间分率法具体怎么计算？
保持时间内的应力历程按时间积分，再除以断裂时间。数学上就是：
$t_{h,j}$ 是第 $j$ 循环的保持时间，$t_d(\sigma, T)$ 是该应力 $\sigma$ 和温度 $T$ 下的许可蠕变破断时间。
但保持时间中应力是变化的吧？ 听说应力松弛会发生…
完全对！ 实际运行中多为应变控制保持。热膨胀产生的应变被约束固定，此时应力随蠕变而随时间衰减。应力随时间的变化满足：
这个常微分方程需要数值求解，获得每时刻的应力 $\sigma(t)$，再代入 $t_d(\sigma, T)$ 进行积分。如果忽视应力松弛而用初始应力当常数，蠕变损伤会被严重高估，这是常见的错误。
哦，所以"保持时间÷破断时间"这种简单算法会出问题…
正是。蠕变速率通常用Norton幂律表示：
$A$ 是材料常数，$n$ 是应力指数（通常3〜8），$Q$ 是活化能，$R$ 是气体常数，$T$ 是绝对温度。温度相关性是指数的，仅差50℃就能导致蠕变速率变化一个数量级。
应变范围分割法（SRP）
线性损伤和以外还有别的方法吗？
有的。NASA在1970年代开发的应变范围分割法（Strain Range Partitioning, SRP）就是一个。它把非弹性应变范围分为4种：
- $\Delta\varepsilon_{pp}$：拉伸塑性 → 压缩塑性（纯疲劳）
- $\Delta\varepsilon_{cc}$：拉伸蠕变 → 压缩蠕变
- $\Delta\varepsilon_{pc}$：拉伸塑性 → 压缩蠕变
- $\Delta\varepsilon_{cp}$：拉伸蠕变 → 压缩塑性（损伤最大）
各分量都有各自的Manson-Coffin型寿命曲线：
$cp$ 损伤最大是为什么？
拉伸蠕变期间，拉伸方向的孔洞沿晶界生成。之后压缩塑性变形虽然急剧，但已开启的孔洞难以闭合。孔洞易于扩展并相连，加速晶界破坏。例如燃气涡轮叶片在运行时因离心力受拉伸蠕变，停止时急冷产生压缩塑性——典型的 $cp$ 支配模式。
SRP在工程上也用吗？
在航天领域（尤其NASA有关涡轮发动机设计）有实绩。但要实验确认4种应变分割成本很高。因此核电行业主要用ASME NH的线性损伤和法则，SRP仅在NASA系统部分应用，形成了市场分工。
延性消耗法与能量法
还有其他方法吗？
有几种。我来整理一下主要的：
| 方法 | 基本概念 | 主要应用领域 |
|---|---|---|
| 延性消耗法 (Ductility Exhaustion) | 用蠕变延性消费率定义损伤：$D_c = \int \dot{\varepsilon}_{cr} / \varepsilon_f^* \, dt$ | 英国R5规范（EDF/核电） |
| 能量法 | 用滞后回线面积（耗散能）评估损伤 | 研究阶段 |
| Chaboche模型 | 连续体损伤力学（CDM），$D$ 作为内部变量时间演化 | 先进研究·部分商用代码 |
| ASME NH (线性损伤和) | $D_c + D_f \leq D$（包络线） | 核电·压力容器（全球标准） |
延性消耗法的 $\varepsilon_f^*$ 是蠕变延性吗？ 从单轴蠕变试验得？
是的。$\varepsilon_f^*$ 是依赖多轴度和温度的蠕变延性，以单轴蠕变破断试验的破断伸长为基础。英国R5方法用多轴度修正：
多轴度越高，蠕变延性越低。相比时间分率法，延性消耗法更符合物理，特别对P91钢等回火马氏体钢，往往精度更好。
ASME规范明确"蠕变-疲劳"的历史
蠕变与疲劳相互作用首次在设计规范中被正式纳入是1970年代的ASME Boiler and Pressure Vessel Code（Section III, Division 1, Subsection NH）。背景是液态金属冷却高速反应堆（LMFBR）不锈钢配管实测数据表明，单用疲劳法则会过于乐观估计。试验表明600℃带保持时间的蠕变-疲劳寿命有时仅为纯疲劳的1/10以下。这一冲击性结论驱动了采用包络线的保守评估方法。
蠕变-疲劳相互作用的数值计算方法
FEM蠕变-疲劳分析流程
CAE实际分析蠕变-疲劳时是怎样的流程？
大体上分4步：
- 热分析：启动→定常→停止的温度瞬态分析。得到温度历程 $T(\mathbf{x}, t)$
- 结构分析：以温度历程为荷载进行弹塑性蠕变分析。得到应力、应变历程
- 损伤参数计算：从应变范围算 $D_f$，从保持时间应力历程算 $D_c$
- 包络线判定：检查 $(D_c, D_f)$ 是否在ASME NH包络线内
实务中，第2和3之间还要讨论"弹性跟随效应（elastic follow-up）"。约束较弱的部位蠕变应变会集中，忽视会造成严重后果。
弹塑性蠕变分析是在通常弹塑性分析上加什么？
关键是应变分解。全应变分成：
弹性 $\varepsilon_{el}$、塑性 $\varepsilon_{pl}$（流动独立）、蠕变 $\varepsilon_{cr}$（时间相关）、热应变 $\varepsilon_{th}$。FEM求解器在各积分点算出各分量增量，更新应力。蠕变增量积分采用隐式格式（后退Euler法）比较稳定。
蠕变本构则的实现
Norton法则以外还有蠕变本构吗？ 求解器怎么选？
主要蠕变本构对比如下：
| 本构 | 数式 | 特点 |
|---|---|---|
| Norton法则 （定常蠕变） | $\dot{\varepsilon}_{cr} = A\sigma^n$ | 最简单，仅定常状态。二次蠕变段 |
| 时间硬化 | $\dot{\varepsilon}_{cr} = A\sigma^n \cdot m \cdot t^{m-1}$ | 含一次蠕变，$t$ 是时间 |
| 应变硬化 | $\dot{\varepsilon}_{cr} = (A\sigma^n)^{1/m} \cdot m \cdot \varepsilon_{cr}^{(m-1)/m}$ | 抗应力变动能力强，工程常用 |
| Theta投影法 | $\varepsilon = \theta_1(1-e^{-\theta_2 t}) + \theta_3(e^{\theta_4 t}-1)$ | 含三次蠕变，外推性好 |
| 统一粘塑性模型 (Chaboche等) | 塑性与蠕变统一 | 可内在表达相互作用 |
工程上应变硬化最均衡，应力变动时行为也合理。Abaqus中 *CREEP, LAW=STRAIN，Ansys中 TB,CREEP,,,,2（应变硬化）设置。
应力松弛与保持时间处理
FEM中怎样处理保持时间的应力松弛？ 时间增分设置得很难似的…
这是实务中最易出错的。保持时间（如1000小时定常运行）分析时：
- 时间增分设置：保持初期应力松弛急剧，需小增分（秒~分钟）。之后可指数扩大增分。自动时间增分（CETOL/CREEP）最安全
- CETOL参数：Abaqus中
*VISCO步用CETOL=0.5e-3左右，能恰当控制蠕变应变增量 - 稳定性：隐式法（后退Euler）即便增分大也稳定。显式法的时间增分约束严苛
1000小时全部FEM求解？ 计算时间太长吧…
不用全解。常见做法是：
- 代表循环分析：开始几循环+定常后1循环，求单循环损伤增量
- 循环外推：损伤增量稳定后，剩余循环线性外推
- 简化法：弹性分析+Neuber法则推弹塑性应变，保持时间仅数值积分应力松弛
第3种简化法在ASME NH附录也推荐，没条件做完整FEM时做这个。
损伤参数的后处理
FEM结果算 $D_c$ 和 $D_f$ 时用什么值？ 节点值还是积分点值？
很关键的问题。答案是积分点（Gauss点）值。节点值平均化平滑，应力集中处损伤易被低估。具体后处理步骤：
- 各积分点每循环提取应变范围 $\Delta\varepsilon_t$（拉伸峰→压缩峰）
- 去弹性部分得非弹性应变范围 $\Delta\varepsilon_{in}$
- 从设计疲劳曲线查 $N_d$，算 $D_f = \sum n_j / N_{d,j}$
- 积分保持时间应力-时间曲线：$D_c = \sum \int dt / t_d(\sigma,T)$
- $(D_c, D_f)$ 在包络线上标点判定
是把全部积分点都算，还是最严苛的地方？
通常找出最高温度部和最大应力范围部进行详评。但蠕变和疲劳最危险的位置可能不同——管内面温高蠕变伤大，外面热梯度应力范围大疲劳伤大。所以候选点要多个。
蠕变-疲劳相互作用的实务应用
应用案例 — 核反应堆、涡轮机、锅炉
实际项目中蠕变-疲劳评估怎样进行？ 具体例子能讲吗？
以火力发电锅炉蒸汽联箱（T管接头处）为例，600℃蒸汽，日启停（DSS）运用场景：
- 运行曲线定义：冷启动（室温→600℃，4小时）、定常（600℃，16小时保持）、停止（600℃→200℃，4小时）
- 热过渡分析：内外表温差最大约80℃，这是热应力源
- 弹塑性蠕变分析：T管内侧角是应力集中点。启动时内侧拉应力→定常时应力松弛→停止时转压应力
- 损伤评估：年250循环×40年=10000循环，算 $D_f$ 和 $D_c$，判定包络线
整个分析要花多久？
通常1周左右：模型建立·边界条件1~2天，热分析几小时，结构分析数小时~1天（网格大小而定），后处理1天。但材料数据收集·验证另算，数据不足要委外试验。
网格、边界条件、荷载循环设置
网格要多细？ 普通强度分析有区别吗？
蠕变-疲劳有特殊要求：
- 厚度方向：最少8~10层。温度和应力梯度都要捕捉准确。4层不够
- 应力集中部：焊缝止端、喷嘴接合、管接头内角局部细分。要素大小≤板厚的1/20
- 要素类型：二阶六面体（20节点）为佳。一阶要素会剪切锁定，应变低估
- 轴对称：配管可用轴对称模型大幅降低成本
荷载循环直接用实际运行数据吗？
实务中代表性运行方式分类：比如
- 冷启动（室温→额定温度）：年10回
- 温启动（200℃→额定温度）：年50回
- 热启动（400℃→额定温度）：年200回
- 负荷变动（额定的50~100%）：年1000回
各类型应变范围不同，都要分别用对应的 $N_d$ 评估，再合算。用雨流计数法把不规则荷载分解成循环是标准做法。
材料数据的获取与注意事项
材料数据（如Norton法则的 $A$, $n$）怎样获得？
主要数据来源如下：
| 数据源 | 内容 | 注意 |
|---|---|---|
| ASME Section II | 许用应力、设计疲劳曲线 | 设计值偏保守，非平均值 |
| NIMS（日本） | 蠕变破断数据、10万小时级 | 免费公开，全球最大蠕变DB |
| ECCC（欧州） | 蠕变数据表 | P91、P92等新钢种强 |
| 自有试验 | 特定炉批特性 | 最可信但高成本 |
关键是Norton参数 $A$, $n$, $Q$ 随温度变化。引用文献值必须确认温度范围。蠕变破断数据服从对数正态分布，平均值与下限值（如-20%应力补正）结果差很大。
常见失败与对策
初学者常犯什么错？ 提前知道好…
| 失败模式 | 原因 | 对策 |
|---|---|---|
| 蠕变损伤极大 | 忽视应力松弛，用初始应力积分 | 应变控制保持时必算应力松弛 |
| 疲劳寿命过于乐观 | 仅弹性分析，忽视非弹性应变 | 弹塑性蠕变分析，含非弹性部 |
| 计算发散 | 蠕变时间增分过大 | 设CETOL（Abaqus）或CREEP RATE CONTROL（Ansys） |
| 温度差50℃结果翻倍 | Norton法则温度依赖（指数） | 先验证热分布精度，检查传热系数 |
| 网格高度依赖 | 应力集中部网格不足 | 3档网格密度确认收敛性 |
| 包络线判定材料错 | 304不锈钢和316不锈钢包络不同 | 严格确认ASME NH的材料指定 |
发电厂锅炉的"启动次数限制"与蠕变-疲劳
火力发电厂锅炉有"累积启动次数上限"规定。这是因为启停时高温蒸汽温度急变，600℃超的蒸汽联箱遭热疲劳。同时定格运行中蠕变常年进展，复合损伤年年积累。某国内大型火力电站，用"定期无损检验（超声/涡流）+蠕变-疲劳CAE"重新评估了设计30年的蒸汽联箱，成功延命10年。电力自由化后长寿命化需求激增，这类残余寿命评估项目急剧增加。
蠕变-疲劳相互作用的软件对比
Abaqus、Ansys、Marc、COMSOL的对应状况
各软件做蠕变-疲劳的优缺点是什么？
| 功能 | Abaqus | Ansys Mechanical | MSC Marc | COMSOL |
|---|---|---|---|---|
| Norton/幂律蠕变 | 标准搭载 | 标准搭载 | 标准搭载 | 标准搭载 |
| 应变硬化 | *CREEP, LAW=STRAIN | TB,CREEP,,,,2 | CREEP MODEL 2 | 用户定义 |
| 时间硬化 | *CREEP, LAW=TIME | TB,CREEP,,,,1 | CREEP MODEL 1 | 用户定义 |
| 统一粘塑性 (Chaboche等) | UMAT支持 | USERMAT支持 | HYPELA2支持 | PDE直接写 |
| 自动时间增分 | CETOL | CREEP RATIO | AUTO STEP | 标准搭载 |
| 应力松弛追踪精度 | 高（Visco step） | 高 | 高 | 普通 |
| 后处理损伤计算 | Python (Abaqus/CAE) | APDL / ACT | Mentat Post | COMSOL Post |
| ASME NH自动评估 | 用户实现 | 用户实现 | 用户实现 | 用户实现 |
实际上任何求解器都不自动做ASME NH包络线判定。所有工具都需要后处理时用Python/MATLAB自写脚本。部分第三方工具（fe-safe、nCode等）支持疲劳自动评估，但蠕变损伤的时间分率积分通常还是手工。
UMAT、USERMAT是什么？ 自己写程序吗？
UMAT（Abaqus）或USERMAT（Ansys）是用户自定义材料子程序。用Fortran或C写，编译链接到求解器。Chaboche粘塑性模型或独特的蠕变-疲劳损伤模型都可实现。
标准蠕变模型（Norton、应变硬化等）够用的话UMAT不需要。但P91钢这种焙火马氏体的复杂蠕变行为——组织劣化伴随蠕变加速——就必须靠UMAT才能表达准确。
选择指南
结局蠕变-疲劳评估最合适的工具是哪个？
"最佳"取决于用途，根据产业实绩：
- Abaqus：高温非线性实绩最丰，UMAT/UEL自由度大。研究到工程都行。核电产业事实标准
- Ansys Mechanical：Workbench环境易上手。APDL自动化也可。电力公司工程师友好
- MSC Marc：大变形·接触强，有垫圈的法兰蠕变等接触耦合问题强
- COMSOL：PDE直接写灵活，研究向。工业应用实绩少
- Code_Aster（开源）：EDF开发，R5法实装。免费但法文技术支持为主
ANSYS Creep与Abaqus UMAT的选择
蠕变-疲劳分析工具选型的关键是"套装蠕变模型"还是"用户子程序（UMAT/USERMAT）"。Ansys Mechanical内置时间硬化·应变硬化蠕变本构，通过GUI即可使用，优点是快速上手。Abaqus的UMAT自由度更大，高端粘塑性模型实现容易。实务差异在"非标合金的蠕变行为"。P91钢（9Cr-1Mo）属焙火马氏体组织，随时间经年劣化，蠕变特性复杂，标准蠕变模型中期精度不足，需UMAT自定义才能达标。
蠕变-疲劳相互作用的前沿研究
最新研究动向
蠕变-疲劳领域今后怎么发展？ 新方法有吗？
几个值得注意的趋势：
- 连续体损伤力学（CDM）模型：损伤变量 $D$ 当内部变量纳入本构方程，蠕变和疲劳损伤统一记述。Lemaitre-Chaboche型代表
- 晶体塑性有限元法（CPFEM）：多晶粒界滑移和孔洞形成用微观尺度直接模拟。成本大但损伤机理深刻
- 机器学习寿命预测：少量短期试验数据推长期寿命。高斯过程回归和神经网络数据外推
- 相场法粒界损伤：粒界孔洞核生、生长、连结用扩散界面模型跟踪，无网格依赖
- 数字孪生：实机传感器数据（温度、应变、振动）实时反馈FEM，动态更新残余寿命
数字孪生很诱人呢！ 实时寿命预报就能提前保养…
是的！TBM（定时保养）→CBM（状态监测保养）转换里，数字孪生是关键。稼働率和安全性都提高。但难点多：
- 实时FEM计算成本高→代理模型（ROM, Reduced Order Model）活用
- 传感数据噪声处理与不确定度定量化（UQ）
- 监管机构是否把数字孪生结果认可为设计根据——制度整备
日本原子力规制委员会对数字孪生很谨慎，但欧洲EDF（法国电力）已在自有核反应堆用Code_Aster基础的在线监控系统。这块10年内肯定大动。
蠕变-疲劳相互作用的故障排除
常见错误与对策
先生，蠕变-疲劳分析最易踩的坑是什么？ 事先知道…
按求解器分常见问题：
Abaqus情况：
- "Time increment required is less than the minimum specified"：蠕变增分发散。
CETOL改小（如5e-4→1e-4），或初始时间增分设小 - "Excessive distortion at integration point"：蠕变应变累积要素畸变过度。检查
NLGEOM=YES，网格细化 - 保持时间结果平坦：蠕变参数太小或温度偏低。检查参数单位（时间vs秒）
Ansys Mechanical情况：
- "Element has become highly distorted"：同上。增加
NSUBST细分荷载步 - 蠕变不起效：确认
TB,CREEP激活温度。温度参考表设置对否 - "Solution not converged"：蠕变非线性强，加
NEQIT（最大迭代），调整CNVTOL
"分析结果与实验不合"时怎办？
蠕变-疲劳特有的调试顺序：
- 先单轴蠕变试验重现：1要素模型+单轴恒应力恒温。验证蠕变参数实现无误。合不上是材料数据问题
- 次单轴疲劳试验重现：应变控制低周疲劳试验的滞后回线。检查循环硬化/软化表达
- 保持疲劳试验重现：含拉伸保持的蠕变-疲劳试验寿命。合不上是损伤评估方法问题
- 最后结构问题上升：逐级加网格密度检验收敛性
从1要素开始的"减法调试"很关键呢。一开始复杂形状就原地踏步…
正是。蠕变-疲劳并联了材料、荷载、网格、时间增分4个因素，1个不合就迷宫化。减法确认本质。特 Norton参数 $A$, $n$, $Q$ 文献间单位制不同（$\sigma$ 是MPa还是Pa，时间是秒还是小时）很常见，1要素测试必须验。