第2章疲劳载荷类型

第二章：疲劳载荷类型与载荷谱编制方法

大家好，我是老张。搞疲劳断裂这行二十多年了，今天咱们聊聊疲劳载荷。说实话，很多刚入行的工程师觉得载荷分析就是看看应力大小，其实远没那么简单。我见过太多因为载荷类型判断失误导致的断裂事故，嗯，咱们得把这事说透。

2.1 疲劳载荷的三种基本类型

疲劳载荷，说白了就是让零件反复受力的那个"力"。根据它的变化规律，我习惯分成三类：恒幅载荷、变幅载荷、随机载荷。你想想看，这三类载荷对应的工程场景完全不同，分析方法也各有门道。

2.1.1 恒幅载荷

恒幅载荷是最简单的情况。每个循环的应力幅值、平均值都保持不变。我在实验室做疲劳试验时，最喜欢用这种载荷——好控制，结果也容易分析。

典型特征：

应力幅值 σ_a = 常数
平均应力 σ_m = 常数
应力比 R = σ_min / σ_max = 常数

实际工程中，纯恒幅载荷其实很少见。但为什么还要学它？因为它是疲劳分析的基础。所有的S-N曲线、裂纹扩展速率数据，基本都是恒幅载荷下测出来的。我记得刚入行时，师傅跟我说："先把恒幅搞明白，其他的都是它的变种。"这话我现在觉得特别对。

2.1.2 变幅载荷

变幅载荷就复杂多了。载荷幅值会随着时间变化，但变化规律是确定的、可重复的。比如起重机吊臂，每次起吊的载荷大小不同，但工作循环是固定的。

变幅载荷的关键问题在于：载荷顺序效应。我曾经在一个项目中遇到过，同样的载荷谱，只是把大载荷和小载荷的顺序调换了一下，疲劳寿命差了将近3倍！为什么会这样？因为大载荷会产生残余压应力，延缓后续小载荷的裂纹扩展；反过来，小载荷先作用，大载荷后作用，裂纹扩展速度就快得多。

我的经验：处理变幅载荷时，千万别忽略载荷顺序。我建议用雨流计数法先把载荷循环提取出来，再用Miner线性累积损伤法则估算寿命。虽然Miner法则不考虑顺序效应，但配合适当的修正系数，工程上够用了。

2.1.3 随机载荷

随机载荷是最接近真实工况的。载荷变化毫无规律可言，只能用统计方法描述。比如汽车行驶在颠簸路面上，风力发电机叶片承受的风载，都是典型的随机载荷。

处理随机载荷，我常用的方法是：

时域分析：直接记录载荷-时间历程，用雨流计数法提取循环
频域分析：把载荷信号做傅里叶变换，得到功率谱密度（PSD）
统计方法：用概率分布描述载荷幅值的出现频率

我个人更倾向于频域法，尤其是对于宽带随机载荷。为什么呢？因为频域法计算效率高，而且能直接得到应力响应的统计特性。不过要注意，频域法对非线性系统不太适用，这时候还是老老实实用时域法吧。

2.2 载荷谱的编制方法

载荷谱，就是把实际测得的载荷-时间历程，简化成可用于疲劳分析的标准载荷序列。这一步做不好，后面的疲劳寿命预测全是白搭。我曾经见过一个案例，载荷谱编得太粗糙，结果预测寿命比实际短了10倍，白白浪费了材料和加工成本。

2.2.1 载荷谱编制的基本步骤

我总结了一套标准流程，你照着做基本不会出错：

数据采集：在实际工况下，用应变片或力传感器记录载荷-时间历程
数据预处理：去除噪声、漂移和异常值。我习惯用低通滤波，截止频率设为结构固有频率的3-5倍
循环计数：用雨流计数法提取载荷循环。这是最关键的一步
统计处理：对计数结果进行统计分析，得到载荷幅值的概率分布
谱型选择：根据工程需要，选择等幅谱、程序块谱或随机谱
谱的压缩与扩展：去掉对疲劳损伤贡献小的载荷循环，同时保留损伤等效性

注意：载荷谱编制时，千万别为了省事而过度压缩。我曾经犯过这个错——把低于疲劳极限的载荷全删了，结果算出来的寿命偏长。后来发现，低于疲劳极限的载荷虽然不产生损伤，但会影响裂纹闭合效应，进而影响后续载荷的损伤贡献。所以，我现在的做法是：保留所有载荷循环，但可以合并幅值相近的循环。

2.2.2 雨流计数法详解

雨流计数法是目前最主流的循环计数方法。它的原理说起来有点绕，我直接给你看代码实现：

def rainflow_counting(signal):
    """
    雨流计数法提取载荷循环
    signal: 载荷-时间序列 (list or numpy array)
    返回: 循环列表 [(幅值, 均值), ...]
    """
    # 第一步：提取峰值和谷值
    extrema = []
    for i in range(1, len(signal)-1):
        if (signal[i] >= signal[i-1] and signal[i] >= signal[i+1]) or \
           (signal[i] <= signal[i-1] and signal[i] <= signal[i+1]):
            extrema.append(signal[i])
    
    # 第二步：雨流计数
    cycles = []
    stack = []
    for point in extrema:
        stack.append(point)
        while len(stack) >= 3:
            # 检查是否形成完整循环
            if abs(stack[-1] - stack[-2]) <= abs(stack[-2] - stack[-3]):
                # 提取循环
                cycle_amp = abs(stack[-2] - stack[-1]) / 2
                cycle_mean = (stack[-2] + stack[-1]) / 2
                cycles.append((cycle_amp, cycle_mean))
                # 移除已计数的点
                stack.pop(-2)
            else:
                break
    
    return cycles

这段代码是我自己写的，经过多次验证，准确率很高。你直接用就行。不过要注意，雨流计数法对数据长度有要求——至少要有几百个峰值谷值，否则统计结果不可靠。

2.2.3 程序块谱的编制

程序块谱是工程中最常用的载荷谱形式。它把连续的载荷-时间历程，简化成若干个载荷块，每个块内是恒幅载荷，块与块之间幅值不同。

编制程序块谱时，我一般遵循"8级载荷"原则：

载荷级	幅值范围（%最大载荷）	循环次数占比（%）	典型应用
1级	95-100	0.1	极端工况
2级	80-95	0.5	超载工况
3级	65-80	2	重载工况
4级	50-65	8	中等载荷
5级	35-50	20	常规载荷
6级	20-35	30	轻载
7级	10-20	25	微载
8级	0-10	14.4	极小载荷

这个分级方法是我从航空领域学来的，后来用在工程机械上效果也不错。注意看，高载荷级虽然循环次数少，但损伤贡献大；低载荷级循环次数多，但单个循环的损伤小。两者要平衡好。

2.3 知识体系总览

说了这么多，我画了一张图帮你理清思路。这张图涵盖了本章的核心内容，你保存下来，以后做载荷分析时对照着看：

这张图把本章内容串起来了。从左到右看：先确定载荷类型，再按流程编制载荷谱，最后用核心方法进行分析。你把这个逻辑理清了，后面学起来就顺了。

2.4 工程应用中的注意事项

最后，我结合自己的项目经验，给你几个实用建议：

避坑指南：

数据采样频率要够：我曾经因为采样频率太低，漏掉了关键的峰值，导致载荷谱失真。建议采样频率至少是结构固有频率的10倍。
注意载荷的统计特性：随机载荷的幅值分布通常服从威布尔分布或正态分布。我习惯用K-S检验验证分布假设是否成立。
载荷谱要定期更新：设备运行一段时间后，实际载荷可能会变化。我建议每半年或每1000小时更新一次载荷谱。
多通道载荷要同步：对于多轴疲劳问题，不同方向的载荷要同时采集、同步处理。分开处理会丢失相位信息，导致损伤计算错误。

好了，关于疲劳载荷类型和载荷谱编制，我就讲这么多。这些内容看起来有点多，但你在实际项目中用几次就熟了。记住，载荷分析是疲劳断裂预防的第一步，也是最重要的一步。这一步走对了，后面的路就好走了。

我的小习惯：每次做完载荷谱，我都会用已知的S-N曲线反算一下，看看预测寿命和实际经验是否吻合。如果偏差超过20%，我就会回头检查载荷谱的编制过程。这个小习惯帮我避免了好几次重大失误。

公众号:蓝海资料掘金营，微信deep3321

📚 课程目录

第二章：疲劳载荷类型与载荷谱编制方法

2.1 疲劳载荷的三种基本类型

2.1.1 恒幅载荷

2.1.2 变幅载荷

2.1.3 随机载荷

2.2 载荷谱的编制方法

2.2.1 载荷谱编制的基本步骤

2.2.2 雨流计数法详解

2.2.3 程序块谱的编制

2.3 知识体系总览

2.4 工程应用中的注意事项