首页/文章/ 详情

S-N 曲线为什么能算疲劳寿命？从循环应力试验到有限元寿命云图

兵心依旧

7天前浏览296

做疲劳分析时，很多人都会有一个疑问：有限元明明算出来的是应力云图，为什么到了疲劳模块里，却能变成“寿命云图”？一个零件某处应力是 180 MPa，软件为什么能给出 20 万次循环，或者 3 万次循环这样的寿命结果？

这件事的核心，不是有限元在“预测未来”。工程师早已在试验机上测过材料在不同循环应力下能承受多少次。S-N 曲线把这些试验经验整理成应力与寿命之间的关系；有限元再把结构各位置的受力状态逐一送去与这条曲线比较。

理解这条链路，才能看懂疲劳寿命云图到底在表达什么。

疲劳寿命首先是试出来的

设想一根标准试样装在疲劳试验机上。试验机不断施加交变载荷：拉一下、松一下、再拉一下；也可以是拉压交替，或弯曲、扭转交替。试样表面起初看不出异常，循环足够多后，局部可能形成微小裂纹。裂纹扩展到一定程度，试样最终断裂。

每组试验至少要记录两类信息：试样承受的循环应力水平，以及试样达到规定失效状态前经历的循环次数。再换一根试样，降低应力水平，重复多组试验。通常会得到一个很直观的结果：应力越高，断裂越快；应力越低，寿命越长。

概念示意：不同循环应力下的断裂试验结果，构成 S-N 曲线。

把各组试验数据画在坐标中，横轴是循环次数 N，纵轴是应力水平 S，得到的就是 S-N 曲线，也叫应力—寿命曲线或 Wöhler 曲线。它回答的是：在规定的材料状态、载荷形式、环境和应力比条件下，这种材料或结构细节面对给定循环应力，大约能承受多少次循环才达到疲劳失效。

例如，下面只作变量关系示意，不代表任何具体材料数据：当交变应力由较高水平降到较低水平时，曲线上对应的寿命可能从 10^4 次提高到 10^6 次以上。疲劳寿命不是软件凭空算出来的，S-N 曲线已经把应力水平和失效循环次数的关系提前记录下来。

曲线为什么常用对数坐标

疲劳寿命的跨度很大。有些零件几千次循环后就断裂，有些要经历几百万次甚至上亿次循环。普通坐标下，低寿命区的数据很容易挤在一起，难以同时看清不同寿命区间的变化。

因此，S-N 曲线通常把循环次数 N 放在对数坐标上。横轴常写成 10^3、10^4、10^5、10^6、10^7，分别表示一千次、一万次、十万次、一百万次和一千万次循环。

σa = C × N^b

在高周疲劳、循环响应以弹性为主的范围，S-N 曲线常可用上式近似。σa 是应力幅，N 是疲劳寿命，C 和 b 是由试验数据拟合得到的参数。应力幅也可以由一个循环中的最大、最小应力得到：

σa = (σmax − σmin) / 2

概念示意：应力循环的基本量，决定了后续怎样查询 S-N 曲线。

把 Basquin 关系放到双对数坐标中，曲线常可近似为直线，便于试验数据拟合和寿命反查。需要注意的是，S-N 曲线不是材料的“身份证号码”。表面粗糙度、热处理状态、尺寸、缺口、焊缝、腐蚀环境和载荷比不同，曲线都可能明显变化。实际工程中常用的也不一定是母材曲线；焊接接头、螺栓孔和齿根等结构细节，可能需要各自对应的疲劳曲线。

有限元补上了试样与结构之间的那一步

试验已经说明某个循环应力水平大致对应某个寿命，但真实零件不是一根光滑试样。圆角、孔洞、焊趾、台阶、连接件和边界条件会让不同位置承受不同应力，危险位置往往就在局部应力集中的区域。

这正是有限元的作用。对结构施加载荷、约束和材料参数后，有限元会计算各位置的应力响应。以带孔板拉伸为例，远离孔洞的位置可能是较低的名义应力，孔边则会出现更高的局部应力。若载荷循环变化，还需要得到该位置每个循环中的最大、最小应力，并进一步计算应力幅与平均应力：

σa = (σmax − σmin) / 2　　σm = (σmax + σmin) / 2

σm 是平均应力，反映载荷整体偏拉还是偏压。接下来，疲劳分析软件会在每个计算位置重复一条流程：读取应力循环；按选定疲劳准则把应力状态转换成等效疲劳应力；在匹配的 S-N 曲线上查询允许寿命；再把各位置的寿命绘制到模型上。

有限元负责获得局部循环应力，S-N 曲线负责给出该应力对应的寿命区间。

寿命云图不是有限元直接“算出裂纹长了多久”。有限元提供局部应力，S-N 曲线提供应力—寿命关系，二者结合后得到寿命预测结果。

寿命云图中的红色区域表示什么

假如寿命云图的图例单位为 cycles，某处显示 Nf = 8.0 × 10^4，它表示：在当前材料疲劳数据、载荷循环、应力修正方法和疲劳准则下，该位置预计在约 8 万次循环附近达到疲劳失效条件。

图中颜色通常只是寿命大小的相对表示。红色区域对应预测寿命最短的位置，因而应优先检查；蓝色区域寿命更长，疲劳风险相对较低。但红色不等于这个位置一定会最先在真实零件上断裂。

它只说明，在当前模型和计算假设下，这里的疲劳损伤最严重。网格、约束、接触、载荷路径或材料数据有问题时，红区的位置也可能不可信。因此，看到红区后真正该问的是：这里的应力是否具有物理意义，能否用于疲劳评估？

不是所有应力峰值都能直接代入 S-N 曲线

固定约束边缘、尖角、点载荷和接触边缘附近，常会出现很高的有限元应力峰值。理想尖角或刚性约束处，网格越细，最大应力可能继续上升，甚至没有稳定值。这类结果常与应力奇异有关。

若直接把这类峰值拿去查 S-N 曲线，软件会得到很短的寿命，甚至接近零循环。它并不必然表示真实零件会立刻断裂，而是说明模型中的理想化边界使局部应力失去了直接评价的意义。

概念示意：先确认应力是否具有物理意义，才能把它送入疲劳曲线。

不同疲劳问题需要匹配不同的应力处理口径。光滑或一般机械零件可采用名义应力法或局部应力法；焊接结构常采用热点应力法、结构应力法或疲劳等级曲线；明显弹塑性循环时，可能需要应变—寿命法，即 ε-N 方法；多轴非比例载荷下，还可能需要临界平面法等准则。

因此，疲劳分析不是把 von Mises 应力放进软件、点击求解这么简单。S-N 曲线必须和应力定义方式匹配。把焊缝疲劳曲线配普通实体单元的节点峰值，或者把约束边缘的奇异应力直接用于寿命计算，结果通常没有工程意义。

变幅载荷下，寿命怎样累积

真实工况很少只有固定应力幅循环。汽车经过不同路面、风机经历不同风速、桥梁承受不同车流时，载荷水平一直变化。一段时间应力幅很大，另一段时间又比较小。

这类问题通常先做雨流计数，把复杂应力—时间历程分解成不同应力幅、不同出现次数的等效循环。然后计算累积损伤，最常用的是 Palmgren–Miner 线性累积规则：

D = Σ(ni / Ni)

ni 是第 i 类循环实际发生的次数，Ni 是该应力水平从 S-N 曲线查得的失效循环次数，D 是累计疲劳损伤。当 D 接近 1 时，工程上通常认为该位置达到疲劳失效条件。

举例说，某应力水平允许承受 10 万次，而实际已经经历 2 万次，则这一等级已消耗约 20% 的寿命，损伤为 0.2。不同应力等级造成的损伤继续累加，就能得到结构的损伤分布或剩余寿命估计。这也是许多疲劳模块除 Life 云图外，还能输出 Damage 云图和 Safety Factor 云图的原因。

平均应力与红色云图，最后还要判断什么

同样是 100 MPa 的应力幅，一种循环可能是拉压交替，另一种则始终在拉伸状态附近波动，它们的疲劳寿命并不一定相同。疲劳裂纹通常更容易在拉伸阶段张开和扩展，平均拉应力升高时，材料可承受的交变应力幅往往会降低。

因此，实际疲劳计算还要考虑平均应力。常见修正方法包括 Goodman、Gerber 和 Soderberg。以 Goodman 修正为例：

σa,eq = σa / (1 − σm / σu)

σu 是材料抗拉强度，σa,eq 是用于查询 S-N 曲线的等效应力幅。平均拉应力增大时，σa,eq 增大，预测寿命会下降。不同修正模型的适用范围并不完全相同，优先使用与实际载荷比相符的试验数据，通常比套用默认修正更可靠。

同一张静力应力云图并不对应唯一的疲劳寿命。最大应力、最小应力、应力幅、平均应力、循环次数和载荷顺序，都会影响寿命预测。

疲劳模块完成计算，只说明软件在给定输入下得到了结果，不说明结果已经能够用于设计。看到红色区域后，除了平均应力，还要核对载荷工况是否符合实际循环过程，材料 S-N 数据和应力类型是否匹配，热点区域的网格是否合理，以及约束、接触、表面状态、焊接细节和环境影响是否已经处理。

下一次看到模型上一片红色区域时，先问一句：这里的红色，究竟是结构真的危险，还是模型、网格和应力定义共同制造出来的假热点？这个问题，往往比寿命数值本身更值得先分析。

来源：兵哥讲力学

ACT 疲劳断裂汽车焊接裂纹材料试验螺栓

著作权归作者所有，欢迎分享，未经许可，不得转载

首次发布时间：2026-08-19

零基础入门ABAQUS的Python二次开发：找回自学二次开发的畅快之感

¥149

还没有评论