📚 课程目录第2章循环载荷与应力应变描述

2. 循环载荷与应力应变描述

在结构疲劳寿命仿真中，载荷与应力的描述方式直接决定了寿命计算的准确性。工程结构在实际服役中承受的往往不是恒定载荷，而是幅值、均值不断变化的循环载荷。本节以实测载荷谱和有限元节点应力为对象，梳理应力比、幅值、均值、雨流计数与循环应力应变曲线之间的表达与转换关系。

一、循环应力的基本参数

一个等幅循环应力可以用三个参数完整描述：最大应力 σ_max、最小应力 σ_min，以及由它们派生的应力幅 σ_a、平均应力 σ_m 和应力比 R。定义如下：

应力幅：σ_a = (σ_max − σ_min) / 2，反映循环的波动幅度；
平均应力：σ_m = (σ_max + σ_min) / 2，反映循环的中心位置；
应力比：R = σ_min / σ_max，反映循环的不对称程度。

这里约定拉应力为正、压应力为负，与有限元节点应力的符号约定保持一致。当 R = −1 时为对称循环，σ_m = 0；R = 0 时为脉动循环，σ_min = 0；R > 0 表示循环始终处于同号应力区间。工程中常见的焊接接头疲劳试验常取 R = 0.1 或 R = 0.5，而旋转弯曲构件接近 R = −1。

应力比 R、应力幅 σ_a 与平均应力 σ_m 三者并非独立：已知任意两个即可推出第三个。例如已知 σ_a 和 R，则 σ_m = σ_a·(1+R)/(1−R)。

二、从实测载荷谱到应力谱

实测载荷谱通常以时间序列形式记录载荷—时间历程，例如应变片采集的电压信号经标定后得到载荷序列。要将载荷谱转化为可用于疲劳分析的应力谱，需要经过以下关键环节：

载荷标定与去噪：将采集的电压信号通过标定系数转换为物理载荷，剔除零点漂移和电磁干扰引起的毛刺。
峰谷提取：只保留载荷历程中的极大值和极小值，去除中间单调段，形成峰谷序列。这一步可大幅压缩数据量而不丢失疲劳损伤信息。
雨流计数：将峰谷序列分解为一系列完整的应力循环，每个循环对应一个应力幅和平均应力。
应力换算：通过有限元分析或标定关系，将载荷循环转换为关注节点处的应力循环。

雨流计数法是疲劳分析中应用最广泛的循环计数方法。其基本思想是：将载荷—时间历程旋转 90°，使时间轴垂直向下，载荷轴水平，想象雨水沿屋檐流下，凡遇到比当前峰（谷）更高的峰（更低的谷）则停止，形成一个个闭合的滞回环。每个滞回环对应一个完整循环，其幅值和均值可直接读出。

三、雨流计数的实现与结果表达

以一段简化的载荷序列为例，峰谷值为：+80、−40、+120、−60、+40、−20、+100（单位：MPa）。按雨流计数规则处理，可提取出若干完整循环。实际工程中，雨流计数通常由软件自动完成，但理解其输出结构对后续寿命计算至关重要。

雨流计数输出示例（应力单位：MPa）
循环编号   σ_max   σ_min   σ_a    σ_m    R
1          +80     −40     60     20     −0.50
2          +120    −60     90     30     −0.50
3          +40     −20     30     10     −0.50
4          +100    +40     30     70      0.40

计数完成后，通常将结果整理为应力幅—均值二维分布，或按应力比分组统计频次。对于焊接结构，常采用 R = 0.1 的 S-N 曲线，此时需要将其他应力比下的循环通过 Goodman 或 Gerber 等平均应力修正公式等效到 R = 0.1 或 R = −1 的基准状态。

四、循环应力应变曲线

当应力水平较高、局部进入塑性时，仅用应力描述循环行为不够，需要引入循环应力应变曲线。循环应力应变曲线由一系列稳定循环的应力幅—应变幅数据点拟合得到，通常用 Ramberg-Osgood 关系表达：

ε_a = σ_a / E + (σ_a / K')^1/n'

其中 E 为弹性模量，K' 为循环强度系数，n' 为循环应变硬化指数。该曲线反映了材料在循环加载下的稳定响应，是应变疲劳分析的核心输入。与单调拉伸曲线不同，循环应力应变曲线可能表现为循环硬化、循环软化或循环稳定，取决于材料状态和加载历史。

五、从计数结果到寿命评估的衔接

雨流计数得到的每个循环（σ_a, σ_m）需要进一步处理才能用于寿命计算。对于应力疲劳（高周疲劳），通常先进行平均应力修正，将非零均值循环等效为零均值或指定 R 的循环，再查 S-N 曲线得到该应力幅下的许用循环数 N_i，最后按 Miner 线性累积损伤法则求和：D = Σ(n_i / N_i)。

对于应变疲劳（低周疲劳），则需要将应力幅通过循环应力应变曲线转换为应变幅，再用 Coffin-Manson 关系计算寿命。此时循环应力应变曲线的准确性直接影响应变幅的换算精度，因此 K' 和 n' 的取值应来自材料手册或实测数据，而非随意估计。

雨流计数要求载荷序列为完整的峰谷序列。若直接对原始采样点计数，会因中间单调点产生大量虚假小循环，导致损伤被高估。务必先完成峰谷提取再计数。

六、工程示例：某支架节点应力谱处理

假设某支架在有限元分析中关注节点 1024，实测载荷谱经峰谷提取后得到 8 个峰谷值，雨流计数输出 4 个循环，结果如上表所示。其中循环 1 和循环 2 的 R = −0.5，循环 3 的 R = −0.5，循环 4 的 R = 0.4。若采用 R = 0.1 的 S-N 曲线，需对循环 4 进行 Goodman 修正：

σ_a,eq = σ_a / (1 − σ_m/σ_b)

设材料抗拉强度 σ_b = 600 MPa，则循环 4 的等效应力幅为 30 / (1 − 70/600) ≈ 33.96 MPa。修正后所有循环均以 R = 0.1 的基准进行寿命查表，再按 Miner 法则累积损伤。这一流程将实测载荷谱、有限元节点应力、雨流计数和 S-N 曲线串联为完整的疲劳评估链路。