3. 材料疲劳性能数据体系

在结构疲劳寿命仿真中,材料疲劳性能数据是决定分析可信度的第一道关口。无论有限元模型多么精细、载荷谱多么真实,如果材料曲线本身选错、拟合不当或分散性处理失当,寿命预测结果都可能偏离数倍甚至一个数量级。本节围绕四类核心数据——S-N曲线、ε-N曲线、da/dN-ΔK曲线和疲劳极限——梳理其获取、拟合、分散性处理与工程选用准则,帮助建立一套可追溯、可复现的材料数据体系。

一、四类曲线各自回答什么问题

S-N曲线描述应力水平与疲劳寿命的关系,适用于高周疲劳(通常寿命大于10⁴~10⁵次)和应力控制场景,如焊接接头、螺栓、弹簧等。其横轴为疲劳寿命N(对数),纵轴为应力幅σₐ或应力范围Δσ,常用Basquin关系σₐ = σ_f′(2N)^b拟合,其中σ_f′为疲劳强度系数,b为疲劳强度指数。

ε-N曲线描述应变幅与寿命的关系,适用于低周疲劳和应变控制场景,如压力容器开孔区、涡轮盘榫槽等局部屈服部位。它通常用Coffin-Manson关系表达:Δε/2 = σ_f′/E·(2N)^b + ε_f′·(2N)^c,总应变幅由弹性项和塑性项叠加,ε_f′为疲劳延性系数,c为疲劳延性指数。

da/dN-ΔK曲线描述裂纹扩展速率与应力强度因子范围的关系,用于损伤容限分析。曲线分为三个区域:近门槛区、Paris稳定扩展区和快速扩展区。工程上最常用Paris公式da/dN = C(ΔK)^m,其中C、m为材料常数,ΔK为应力强度因子范围。

疲劳极限σ₋₁或Δσ_D则回答“低于多少应力可以认为寿命无限”。对于钢铁材料,S-N曲线在10⁶~10⁷次附近趋于水平;对于铝合金等无明显疲劳极限的材料,工程上常以10⁷次对应的条件疲劳强度代替。

材料疲劳性能数据体系与工程选用路径 数据来源 标准手册 / 试验 / 文献 S-N 曲线 应力控制 · 高周疲劳 σₐ = σ_f′(2N)^b ε-N 曲线 应变控制 · 低周疲劳 Coffin-Manson 关系 da/dN-ΔK 曲线 损伤容限 · 裂纹扩展 Paris: da/dN = C(ΔK)^m 疲劳极限 / 条件疲劳强度 钢铁:10⁶~10⁷次水平渐近线;铝合金:10⁷次条件强度 拟合与分散性处理 最小二乘拟合 对数坐标下线性回归 存活率 P-S-N 50% / 95% / 99% 曲线 分散性修正 标准差 / 变异系数 安全系数折减 按规范或工程经验 工程选用准则 高周疲劳 · 应力控制 选 S-N 曲线 + 疲劳极限 低周疲劳 · 应变控制 选 ε-N 曲线 + 循环应力应变 损伤容限 · 裂纹扩展 选 da/dN-ΔK + 门槛值

二、数据获取:从试验到可用曲线

获取S-N曲线最常用的方法是旋转弯曲或轴向拉压疲劳试验。以某低合金钢为例,在室温、应力比R=-1条件下,取7~10级应力水平,每级做3~5根试样,记录失效循环数。若试样在10⁷次未断,记为越出(run-out),不能直接当作失效点参与回归,需用极大似然法或升降法单独处理。试验后,在双对数坐标下以σₐ为纵轴、N为横轴描点,用最小二乘法拟合Basquin直线。

ε-N曲线通常采用轴向引伸计控制的低周疲劳试验,应变比Rε=-1,应变幅从0.2%到1.5%分多级。每级记录稳定循环的应力幅和塑性应变幅,分别拟合弹性项和塑性项。需要特别注意的是,ε-N曲线必须同时给出循环应力-应变曲线,否则无法把名义应变转换为局部应变。

da/dN-ΔK曲线则通过紧凑拉伸(CT)或中心裂纹拉伸(M(T))试样获得。按ASTM E647,在恒幅载荷下测量裂纹长度a随循环数N的变化,用割线法或七点递增多项式法计算da/dN,再由应力强度因子公式反算ΔK。近门槛区需采用降K程序,避免载荷历史干扰。

三、拟合:把散点变成设计曲线

拟合不是简单画一条穿过点群中心的线。以S-N曲线为例,工程上通常拟合中值曲线(50%存活率)和设计曲线(95%或99%存活率)。若对数寿命服从正态分布,则P-S-N曲线可写为:

lg N_P = lg N_50 - k_P · s
其中:
  lg N_50 —— 中值对数寿命
  s      —— 对数寿命标准差
  k_P    —— 与存活率P对应的单侧容限系数
           P=95% 时 k≈1.645;P=99% 时 k≈2.326

对于da/dN-ΔK曲线,Paris公式的C和m同样存在分散性。工程上常取95%置信上限曲线用于保守寿命评估,即在同一ΔK下取较高的da/dN。若数据覆盖近门槛区和快速扩展区,则应采用三参数或四参数模型,而不是强行用一条Paris直线贯穿全范围。

四、分散性处理:承认材料不是理想均匀体

疲劳数据的分散性来源包括:材料批次差异、热处理波动、表面加工状态、试样几何、试验机误差和加载频率等。处理分散性有三条常用路径:

  • 统计法:收集足够样本,拟合P-S-N曲线,直接给出不同存活率下的设计曲线。
  • 安全系数法:在中值曲线基础上,对寿命或应力乘以折减系数,如寿命除以5~10,或应力除以1.5~2.0。
  • 规范法:直接采用IIW、Eurocode 3、ASME等标准中已考虑分散性的设计S-N曲线,适用于焊接接头等标准化场景。
工程选用时,优先顺序建议为:有同材料同工艺的实测P-S-N曲线 → 标准规范设计曲线 → 手册中值曲线加安全系数 → 相近材料曲线类比。任何类比都必须说明差异来源和保守性判断。

五、工程选用准则:让曲线与工况匹配

选用材料疲劳数据时,至少核对以下五项:材料牌号与热处理状态、应力比R、加载频率、环境温度与介质、表面状态与缺口。例如,手册中40Cr调质钢的S-N曲线通常对应R=-1、室温、抛光试样;若实际构件是R=0.1、有应力集中、表面喷丸,则不能直接套用,需进行平均应力修正(Goodman或Gerber)和缺口修正(Kf)。

对于低周疲劳,若构件局部应变超过屈服应变,应优先使用ε-N曲线而非S-N曲线;对于含裂纹或允许裂纹存在的结构,应使用da/dN-ΔK曲线进行剩余寿命评估,并以门槛值ΔK_th判断裂纹是否扩展。疲劳极限只适用于高周疲劳且应力幅低于门槛值的场景,不能用于低周或裂纹扩展分析。

常见错误:把不同R下的S-N曲线混用而不做平均应力修正;把da/dN曲线外推到门槛区以下;用中值曲线直接做保守设计而不考虑存活率。这些都会导致寿命预测偏危险。

六、一份可执行的数据选用检查表

检查项 关键问题 判断依据
材料状态 牌号、热处理、批次是否与构件一致? 不一致时需评估差异并加安全系数
应力比R 曲线R是否等于实际工况R? 不等时做Goodman/Gerber修正
寿命区间 属于高周还是低周? N>10⁴~10⁵用S-N;否则用ε-N
分散性 是否有P-S-N或安全系数? 设计曲线取95%或99%存活率
裂纹状态 是否允许裂纹存在? 允许则用da/dN-ΔK做损伤容限
环境 温度、介质、频率是否一致? 腐蚀或高温需专用曲线

把上述四类曲线和选用准则落实到仿真流程中,意味着每一次寿命分析都应能回答:这条曲线来自哪里、对应什么R和存活率、做了哪些修正、为什么适用于当前构件。只有数据体系可追溯,疲劳寿命仿真才具备工程可信度。