简介《疲劳与断裂力学第2章金属材料的常幅疲劳行为》PPT课件面向学习材料力学、疲劳与断裂力学的高年级本科生与研究生聚焦金属材料在交变载荷下的循环应力应变特性、S-N曲线及疲劳极限等核心概念。课件以图文并茂方式系统讲解交变载荷参数、循环滞回环、Bauschinger效应、循环软/硬化行为以及成组试样法、增级试验法确定循环应力应变曲线的方法并延伸平均应力松弛、棘轮行为与非比例附加硬化等进阶内容。整个资源为1个PPT文件压缩包大小约2.3MB内容覆盖第2章全部小节层次清晰便于按节选学。已有101人学习下载适合需要系统掌握金属材料疲劳行为基础理论、理解S-N曲线工程含义及疲劳寿命估计要点的读者。1. 疲劳数据为什么比强度数据更“不讲道理”很多工程师第一次接触疲劳时都会有一个困惑静强度校核明明通过了安全系数给到了 2 以上构件还是在服役几个月后断了。断口上那圈贝壳纹和放射棱几乎就是疲劳破坏的身份证。而真正让设计人员被动的是破坏应力往往远低于材料的屈服强度有时甚至不到抗拉强度的三分之一。这说明金属材料的常幅疲劳行为本质上不是在跟强度打交道而是在跟“循环次数”打交道。本章讨论的就是这一章 PPT 要解决的核心问题在应力幅恒定的循环载荷下金属材料如何从完好到萌生裂纹、再到失稳扩展以及工程师如何用 S-N 曲线、疲劳极限和累积损伤法则把这些现象变成可计算、可查表、可写进设计规范的东西。这套内容既是断裂力学课程里连接材料性能与结构寿命的桥梁也是工程实践中做寿命评估绕不开的底子。2. 金属材料常幅疲劳行为的数学描述与 S-N 曲线拟合2.1 高周疲劳与低周疲劳的分界常幅载荷该用哪一段常幅疲劳指的是应力幅在整个试验过程中基本保持不变、应力比也固定的载荷形式。最常见的对称循环R -1和脉动循环R 0都属于这个范畴。要描述这种载荷下的材料行为第一件事是搞清楚构件落在低周还是高周区间因为两段的控制参量完全不同。低周疲劳通常指循环寿命 N 10⁴ ~ 10⁵ 次的区间这时交变应力水平高塑性应变在总应变中占比明显控制寿命的是应变幅而不是应力幅。高周疲劳则是 N 10⁵ 次的区间名义应力低于屈服强度材料整体处于弹性范围控制寿命的是应力幅。工程上绝大多数机械构件比如轴、螺栓、齿轮齿根、压力容器接管都落在高周疲劳区间这也是 S-N 曲线最常被使用的场景。Luft 这样说过PPT 课件里如果不先讲清楚这条分界线后面的 S-N 曲线、疲劳极限概念就容易让学生误以为单根曲线能包打天下。2.2 S-N 曲线与 Basquin 方程从经验规律到回归参数常幅疲劳行为的核心表达就是 S-N 曲线横坐标是循环次数 N纵坐标是应力幅 σ_a或最大应力 σ_max。在高周区这条曲线在双对数坐标下近似为一条直线对应 Basquin 方程σ_a σ_f × (2N_f)^b其中 σ_f 是疲劳强度系数量纲与应力相同b 是疲劳强度指数对多数金属材料是负值大致在 -0.05 到 -0.15 之间。把方程两边取对数就能得到线性回归形式这也是目前绝大多数疲劳数据处理软件的内部算法。实际工程中常见的是把 S-N 曲线写成另外两种等价形式应力幅形式 σ_a^m × N C或最大应力形式 σ_max^m × N C。这里的 m 和 C 是材料常数m 一般取 312值越大说明曲线越陡也就是寿命对应力越敏感。下表给出常见金属材料的 Basquin 参数大致范围具体数值应以本单位的疲劳试验报告为准材料类别疲劳强度系数 σ_fMPa疲劳强度指数 b疲劳极限/抗拉强度比结构钢8001800-0.06-0.100.350.50铝合金400800-0.08-0.120.250.40钛合金9001400-0.05-0.080.450.60铸铁300600-0.08-0.150.300.45参数 b 的绝对值越大说明材料对循环载荷越敏感寿命曲线越陡。铝合金的 b 值普遍比钢大这也是铝合金构件在交变载荷下寿命分散性更明显的原因之一。2.3 用 Python 对疲劳试验数据做 Basquin 拟合拿到一组σ_a, N_f试验数据后常见做法是在双对数坐标下做线性回归直接得到 Basquin 参数。下面这段代码用 numpy 的最小二乘法完成拟合并输出斜率和截距对应的材料常数import numpy as np # 试验数据[应力幅 MPa, 循环次数] data np.array([ [350, 1.2e5], [320, 3.4e5], [300, 8.9e5], [280, 2.1e6], [260, 6.8e6], [240, 1.5e7], ]) # 取对数x log10(2*Nf)y log10(sigma_a) x np.log10(2 * data[:, 1]) y np.log10(data[:, 0]) # 一次多项式拟合返回 [斜率, 截距] coeffs np.polyfit(x, y, 1) b coeffs[0] # 疲劳强度指数 b log_sigma_f coeffs[1] # log10(sigma_f) sigma_f 10 ** log_sigma_f # 计算拟合优度 R^2 y_pred np.polyval(coeffs, x) residual y - y_pred ss_res np.sum(residual ** 2) ss_tot np.sum((y - np.mean(y)) ** 2) r_squared 1 - ss_res / ss_tot # 计算标准差对数坐标下单位与 y 相同 std_dev np.std(residual, ddof2) print(fBasquin 参数b {b:.4f}, sigma_f {sigma_f:.1f} MPa) print(f拟合优度 R^2 {r_squared:.4f}) print(f对数标准差 s {std_dev:.4f})这段代码的核心逻辑是先把试验数据转换到对数空间然后用一次多项式拟合得到 Basquin 方程的斜率和截距。这里有几个参数需要特别注意b 是负值代码里输出的是直接拟合结果不必取绝对值σ_f 的物理意义是 2N_f 1 时的虚拟应力幅它通常远高于抗拉强度不代表真实应力只是回归参数而已。R² 低于 0.95 时先检查是不是有数据点落在了低周区或过渡区这两段的点不应该参与同一个线性拟合。对数标准差 s 反映数据的分散性s 越大说明材料疲劳性能越不稳定在做设计时对应力修正系数的取值要有心理准备。3. 疲劳极限的测定升降法与分组法的实操口径3.1 疲劳极限的物理前提裂纹不扩展的应力阈值疲劳极限endurance limit是常幅疲劳行为里最有工程价值的概念它代表材料在给定应力比下能承受无限次循环而不失效的最大应力幅。低碳钢、球墨铸铁等材料在 S-N 曲线上存在明显的水平渐近线这就是疲劳极限存在的直接证据。而铝合金、多数不锈钢等材料在双对数坐标下 S-N 曲线持续下降不存在明确的水平段工程上通常用“条件疲劳极限”替代也就是指定寿命比如 10⁷ 次对应的应力幅。为什么有些材料有疲劳极限、有些没有这个问题涉及位错结构和裂纹扩展门槛值的概念。简单说钢材在循环载荷下会发生循环硬化和应变时效使裂纹尖端的塑性区尺寸趋近于零裂纹无法继续扩展而铝合金没有明显的应变时效效应裂纹尖端始终存在可逆滑移微裂纹可以在低于疲劳极限的应力下缓慢长大。做 PPT 课件时这块内容适合放一张铁素体-珠光体钢和铝合金 S-N 曲线的对比图比文字描述直观得多。3.2 升降法测疲劳极限的具体流程与判定工程上测定疲劳极限的标准方法是升降法staircase method也叫 up-and-down 法。它的工作方式非常工程化先预估一个疲劳极限附近的应力水平 σ₀在某级应力下试验一个试样。如果试样在指定循环次数比如 10⁷ 次之前断裂就把应力降低一个台阶 Δσ 再做下一个试样如果试样越过了指定寿命就把应力升高一个台阶。这样后续试样的应力水平完全由前一个试样的结果决定最终所有试验点会围绕疲劳极限的均值上下波动。试验完成后取断裂和未断裂这两组数据中数量较少的那一组参与统计分析。设这一组有 n 个有效数据点最低应力水平为 σ_minΔσ 为应力台阶则疲劳极限的均值按下式估算σ_e σ_min Δσ × (A/n ± 1/2)式中的 A 是对参与统计的那组数据按“编号 × 出现次数”加权求和的结果正负号的选取规则是用断裂组数据时取加号用未断裂组数据时取减号。台阶 Δσ 一般取预估疲劳极限的 3%5%取得太大结果粗糙取得太小浪费试样且容易受数据分散性干扰。有效试样数量建议不少于 68 个否则统计意义很弱。3.3 疲劳极限数据的 Python 处理升降法试验做完后的数据处理非常机械手算容易出错用 Python 处理可以顺便把均值和标准差一起算出来import numpy as np # 试验记录[应力水平 MPa, 结果] 结果1断裂0未断裂 trials np.array([ [300, 1], [280, 1], [260, 0], [280, 1], [260, 0], [280, 0], [300, 1], [280, 0], ]) # 取数量较少的那一组这里未断裂(0)有4个断裂(1)有4个任取一组 target trials[trials[:, 1] 0] stress_levels target[:, 0].astype(int) # 统计各应力水平出现次数 unique, counts np.unique(stress_levels, return_countsTrue) A np.sum(unique * counts) n len(stress_levels) sigma_min unique.min() delta_sigma 20 # 应力台阶单位 MPa # 未断裂组用减号 sigma_e sigma_min delta_sigma * (A / n - 0.5) print(f应力水平分布{dict(zip(unique, counts))}) print(f疲劳极限均值{sigma_e:.1f} MPa)这段代码的核心逻辑是把升降法试验记录整理成程序能读的格式然后按标准公式计算疲劳极限。注意几个关键点参与统计的那一组必须是数量较少的一组这是升降法数据处理的核心约定台阶 Δσ 的取值不在代码里体现它是在试验设计阶段定的数据处理时直接引用即可当断裂组和未断裂组数量相同时两组分别计算的结果理论上应该接近偏差大说明试验数据分散性异常或者 Δσ 取得不合适。另外这里算出的疲劳极限对应 50% 存活率如果要用于设计需要按材料的标准差和要求的可靠度做进一步修正。4. 平均应力与累计损伤常幅疲劳行为向工程寿命估算的延伸4.1 应力比与平均应力S-N 曲线查不到的值怎么换算标准 S-N 曲线通常是在对称循环R -1条件下测出来的但实际构件承受的载荷很少是纯对称循环。转轴除了弯曲应力还有自重引起的平均应力压力容器内压产生的是脉动循环R 0螺栓联接的载荷更是带大预紧力的非对称循环。当平均应力不为零时直接用对称循环的 S-N 曲线查寿命会偏危险必须进行平均应力修正。最常用的三种修正是 Goodman、Gerber 和 Soderberg 公式。Goodman 是线性关系偏保守工程上用得最多Gerber 是抛物线关系更接近试验数据的趋势但不保守Soderberg 用屈服强度替代抗拉强度做分母比 Goodman 还保守适合脆性材料或不允许出现塑性变形的场合。三种公式的统一形式是把实际应力幅折算成等效应力幅 σ_eq修正方法等效应力幅公式特点Goodmanσ_eq σ_a / (1 - σ_m / σ_b)线性、工程常用Gerberσ_eq σ_a / (1 - (σ_m / σ_b)²)抛物线、贴合数据Soderbergσ_eq σ_a / (1 - σ_m / σ_s)最保守、用屈服强度4.2 等效应力幅的换算与代码实现实际计算时平均应力 σ_m 和应力幅 σ_a 可以从载荷谱中分解出来。对正弦波载荷若最大应力为 σ_max、最小应力为 σ_min则 σ_m (σ_max σ_min) / 2σ_a (σ_max - σ_min) / 2。下面这段代码对一组载荷循环做 Goodman 修正import numpy as np # 载荷谱[最大应力 MPa, 最小应力 MPa] spectrum np.array([ [320, 80], [280, 40], [350, 50], ]) sigma_b 620 # 抗拉强度MPa sigma_s 420 # 屈服强度MPa def goodman_equivalent(sigma_max, sigma_min, sigma_b): sigma_m (sigma_max sigma_min) / 2 sigma_a (sigma_max - sigma_min) / 2 if sigma_m sigma_b: return np.inf return sigma_a / (1 - sigma_m / sigma_b) equivalent_amplitudes [] for smax, smin in spectrum: eq goodman_equivalent(smax, smin, sigma_b) equivalent_amplitudes.append(eq) print(fσ_max{smax:3d} MPa, σ_min{smin:3d} MPa, fσ_eq{eq:6.1f} MPa) # 判断哪些循环是危险的 fatigue_limit 200 # 对称循环疲劳极限MPa risk_cycles [eq for eq in equivalent_amplitudes if eq fatigue_limit] print(f超限循环数{len(risk_cycles)} / {len(spectrum)})这里的核心逻辑是先拆解出平均应力和应力幅再用 Goodman 公式折算成对称循环等效应力幅。需要注意的是当最大应力接近抗拉强度时修正后的等效应力幅会急剧增大程序里给 σ_m ≥ σ_b 的情况返回无穷大是为了避免除零错误实际试验中这种工况意味着静强度可能已经失效疲劳分析不再适用。另外Goodman 修正对压应力状态σ_m 0通常不做修正因为压缩平均应力对疲劳寿命是有利的忽略它相当于偏保守。4.3 Miner 线性累积损伤在常幅谱下的应用实际构件承受的载荷不是单一的常幅循环但工程上常用“雨流计数 各级损伤累加”的办法处理变幅载荷。Miner 法则的核心假设是每一级应力循环消耗的损伤量 Dᵢ nᵢ / Nᵢ其中 nᵢ 是第 i 级应力的实际循环次数Nᵢ 是该级应力对应的疲劳寿命从 S-N 曲线查得。总损伤 D Σ(nᵢ / Nᵢ)当 D ≥ 1 时认为构件发生疲劳破坏。Miner 法则是线性累积损伤的典型代表它不关心各级载荷的作用顺序也不考虑载荷之间的相互作用。对于常幅载荷谱应用非常简单如果某个构件的设计寿命是 2000 小时转速 3000 r/min那总循环次数就是 3.6 × 10⁸ 次。通过雨流计数拆成若干个常幅块再按 Miner 法则累加就能估算寿命。代码跟 4.2 节的等效应力幅计算衔接起来把每个等效幅值对应的 Nᵢ 从拟合好的 S-N 曲线上插值出来整套逻辑就完整了。这里必须提醒一个工程上常见的误用Miner 法则里 D 1 只是名义破坏点大量试验数据表明实际破坏时 D 的取值在 0.3 到 3.0 之间波动。做设计时取 D 1 才判破坏是危险的规范做法是引入损伤容限系数把允许的总损伤限制在 0.5 甚至更低。PPT 课件里如果能放一张 D 值分布直方图学生对这个概念的印象会深很多。5. 把常幅疲劳数据整理成能直接上课的 PPT结构、图表与三个容易出错的地方5.1 PPT 数据页的标准结构一份能直接上课用的常幅疲劳行为 PPT数据页应该遵循从原始记录到结论的完整链路不给学生留“这个数从哪来的”的疑问。我一般会把每一组数据的处理过程拆成四页第一页放原始试验记录表包含应力水平、试样编号、循环次数、断口位置第二页放双对数坐标下的散点图和拟合直线标注 Basquin 方程和 R²第三页放按升降法整理的疲劳极限结果含应力水平序列和计算过程第四页放工程应用示例比如某个实际构件在给定载荷谱下的寿命估算。这四页形成一条完整的证据链比单贴一张 S-N 曲线图有说服力得多。5.2 三个容易出错的地方第一个坑是旋转弯曲疲劳极限和轴向加载疲劳极限混用。旋转弯曲试样的应力梯度比轴向加载试样大同样的材料旋转弯曲测出的疲劳极限通常比轴向拉压低 10%20%。查手册时一定要确认数据对应的加载方式否则会把疲劳极限估高。第二个坑是不标注存活率。教科书上的 S-N 曲线中值曲线对应 50% 存活率设计时至少要用 95% 存活率甚至 99% 存活率对应的 P-S-N 曲线直接拿中值曲线做设计意味着有一半的构件会提前失效。第三个坑是 S-N 曲线的横坐标用线性刻度还是对数刻度。疲劳数据跨越几个数量级线性坐标下低周区的点会挤在一起高周区反而拉得很开看上去好像数据很多实际上信息密度很差。统一用双对数坐标回归拟合也方便。5.3 一键输出带参数标注的 S-N 曲线做 PPT 时最常用到的是把拟合结果直接输出为图片。用 matplotlib 可以一次性完成散点、拟合线、参数标注三个动作import numpy as np import matplotlib.pyplot as plt # 承接 2.3 节的拟合结果 x_fit np.linspace(x.min(), x.max(), 100) y_fit np.polyval(coeffs, x_fit) plt.figure(figsize(8, 6)) plt.scatter(x, y, markero, facecolorsnone, edgecolorssteelblue, label试验数据) plt.plot(x_fit, y_fit, r-, linewidth1.5, labelfBasquin: $\\sigma_a {sigma_f:.0f} (2N_f)^{{{b:.3f}}}$) plt.xlabel($\log_{10}(2N_f)$) plt.ylabel($\log_{10}(\\sigma_a / MPa)$) plt.legend() plt.grid(True, whichboth, linestyle--, alpha0.6) plt.tight_layout() plt.savefig(sn_curve.png, dpi300) plt.show()这段代码把拟合曲线和参数直接标注在图上输出 300 dpi 的 PNG 图片插入 PPT 后清晰度足够。需要注意曲线标注里 b 是负值LaTeX 语法里要把负号包在花括号里才能正常显示。图片保存后建议再检查一下坐标轴范围确保散点和拟合线完整显示不要把边缘的数据点裁掉。本文还有配套的精品资源点击获取