📚 课程目录第4章雨流计数法处理随机载荷谱

4. 雨流计数法处理随机载荷谱

随机载荷谱的雨流计数处理

实测载荷时间历程通常是不规则的随机波形，不能直接用于疲劳损伤计算。雨流计数法将随机载荷历程分解为一系列完整的应力循环（闭合滞回环）和半循环，并统计其均值与幅值的二维分布，从而得到可用于Miner累积损伤计算的载荷谱。本节以一段实测随机载荷为例，完整演示雨流计数的实施过程。

载荷数据与预处理

假设从某构件危险点采集到一段应变-时间历程，经名义应力换算后得到载荷序列（单位：MPa），采样点如下（示例数据）：

载荷序列 S = [120, 80, 150, 60, 130, 90, 170, 70, 110, 50, 140, 100, 160, 75, 125]

雨流计数的前提是：序列必须由峰、谷交替组成。若原始数据存在连续上升或连续下降的中间点，需先进行峰谷提取，只保留极值点。上例已满足交替条件，可直接计数。

雨流计数规则

雨流计数法的核心规则可概括为：

将载荷历程旋转90°，使时间轴竖直向下，载荷轴水平，雨水沿屋檐（载荷峰谷）向下流动。
雨水从每个峰（或谷）的内侧开始流动，当遇到比起点更高的峰（或更低的谷）时停止。
若雨水流到对面屋檐，则形成一个完整循环；否则形成半循环。
完整循环的幅值为该循环的应力范围的一半，均值为循环上下限的平均值。

计数结果（示例）：

循环编号	类型	起点	终点	幅值 Sa (MPa)	均值 Sm (MPa)
1	完整循环	80	150	35	115
2	完整循环	60	130	35	95
3	完整循环	90	170	40	130
4	完整循环	70	110	20	90
5	完整循环	50	140	45	95
6	完整循环	100	160	30	130
7	半循环	120	80	20	100
8	半循环	75	125	25	100

均值-幅值二维分布统计

将上述循环按幅值和均值进行区间划分，统计频次，形成二维载荷谱。例如幅值区间取10 MPa，均值区间取20 MPa：

幅值区间 (MPa)	均值区间 (MPa)	频次
20~30	80~100	1
20~30	100~120	1
30~40	100~120	1
30~40	120~140	1
40~50	80~100	1
40~50	120~140	1
20~30	90~110	1
20~30	90~110	1

该二维分布即为后续Miner累积损伤计算的输入载荷谱。每个循环的损伤贡献为 n_i / N_i，其中 N_i 由该幅值和均值对应的S-N曲线确定。

实施要点：雨流计数前必须确保载荷序列峰谷交替；完整循环和半循环需分别记录；半循环在后续损伤计算中通常按半个循环处理，或与其它半循环配对成完整循环。

结果判读与后续使用

从二维分布可以看出：本例中完整循环的幅值主要集中在20~50 MPa，均值分布在80~140 MPa之间；半循环幅值较小（20~25 MPa），对损伤贡献有限。该载荷谱可直接用于Miner线性累积损伤计算：对每个循环，根据其幅值和均值查S-N曲线得到许用循环数 N_i，计算损伤 D_i = n_i / N_i，总损伤 D = Σ D_i。当 D ≥ 1 时判定疲劳失效。

注意：雨流计数得到的半循环不能直接与完整循环等同处理。工程上通常将半循环两两配对成完整循环，或按0.5个循环计入损伤。若半循环数量为奇数，最后一个半循环需单独处理或保守估计。