雨流循环计数技术手册
1. 概述
雨流循环计数是随机载荷与变幅载荷疲劳分析的核心技术，其核心功能是将连续、无规律的应力-时间历史，转化为离散的“等效恒幅循环”集 合。该技术通过识别应力变化过程中的峰值与谷值，遵循特定规则提取有效疲劳循环，量化每个循环的应力范围与均值，为后续基于Miner线性累积损伤理论的疲劳寿命评估提供标准化输入。
工程意义：随机载荷下结构的疲劳损伤由无数次不同大小的应力循环叠加导致，直接基于原始时间序列计算损伤效率极低且精度不足。雨流循环计数通过“去冗余、提核心”的方式，保留对疲劳损伤起主导作用的循环信息，是连接载荷采集与疲劳评估的关键桥梁，广泛应用于汽车、航空航天、工程机械、轨道交通等领域。
2. 核心原理与基本概念
2.1 核心物理意义
疲劳损伤的本质是材料在交变应力作用下的微裂纹萌生与扩展，而微裂纹的扩展速率与单次应力循环的“应力范围”（ΔS = Sₘₐₓ - Sₘᵢₙ）直接相关。雨流循环计数的物理依据是“应力-应变滞后回线”——每一次完整的应力循环都会在应力-应变曲线上形成一条滞后回线，回线的面积对应单次循环的塑性功，即单次循环的损伤贡献。
雨流循环计数的本质的就是通过识别应力时间历史中的峰值、谷值，间接捕捉这些滞后回线的数量与特征参数（应力范围、均值）。
2.2 关键基础概念
• 峰值（Peak）：应力时间序列中，大于相邻两个点应力值的点，即满足Sᵢ > Sᵢ₋₁且Sᵢ > Sᵢ₊₁的点； • 谷值（Valley）：应力时间序列中，小于相邻两个点应力值的点，即满足Sᵢ < Sᵢ₋₁且Sᵢ < Sᵢ₊₁的点； • 应力范围（ΔS）：单次应力循环中最大应力与最小应力的差值，ΔS = Sₘₐₓ - Sₘᵢₙ，是疲劳损伤计算的核心参数； • 应力均值（Sₘ）：单次应力循环中最大应力与最小应力的平均值，Sₘ = (Sₘₐₓ + Sₘᵢₙ)/2，影响疲劳损伤的速率（拉应力均值加速损伤）； • 等效恒幅循环：雨流计数提取的、能等效替代原始变幅循环损伤贡献的恒幅循环，其核心参数为应力范围ΔS和均值Sₘ。
2.3 核心判断规则
雨流循环计数的核心是通过连续三个特征点（峰值或谷值）的应力范围对比，判断是否存在有效循环，具体规则如下：
取重排后序列中的任意三个连续点，记为点1、点2、点3，计算两点间的应力范围：
ΔS₁₂ = |S₁ - S₂| （点1到点2的应力范围）
ΔS₂₃ = |S₂ - S₃| （点2到点3的应力范围）
判断逻辑：
• 若ΔS₁₂ ≤ ΔS₂₃：点1→点2构成一个有效等效循环，提取该循环的ΔS₁₂和Sₘ（(S₁+S₂)/2）；提取后删除点1和点2，将剩余点依次连接形成新序列，重复计数过程； • 若ΔS₁₂ > ΔS₂₃：点1→点2不构成有效循环，放弃该组判断，移动至下一组连续三点（点2、点3、点4）重新计算判断，直至序列中仅剩1个点或无法形成三组连续点。
3. 完整计数流程
雨流循环计数需遵循“预处理→重排→计数→收尾”的标准化流程，确保计数结果的准确性与一致性，具体步骤如下：
3.1 步骤1：数据预处理——提取峰值与谷值
原始应力时间历史包含大量连续数据点，需先剔除无意义的中间点，仅保留峰值和谷值，形成“峰-谷-峰-谷”交替的离散序列。
目的：简化后续计数过程，避免中间点干扰循环识别（疲劳损伤仅与峰值和谷值决定的应力范围相关）。
示例：将连续的载荷历史（图1）转换为仅含峰值和谷值的序列（图2），得到点1、点2、点3、点4四个特征点。
图1. 原始连续载荷历史
图2. 提取峰值和谷值后的离散序列
3.2 步骤2：序列重排——以绝对值最大特征点为起点
为确保优先提取应力范围最大的循环（最大循环对疲劳损伤贡献占比最高，若遗漏会导致评估结果偏于危险），需对峰值-谷值序列进行重排：
1. 遍历所有峰值和谷值，找到绝对值最大的点（若峰值和谷值绝对值相等，可任选其一）； 2. 以该最大绝对值点为新序列的起点，将该点之后的所有特征点按原始时间顺序跟随其后，再将该点之前的所有特征点按原始时间顺序接在末尾，形成重排后的新序列； 3. 对重排后的序列重新编号，便于后续计数。
示例：图2中的心4为绝对值最大的峰值，将其移至新序列起点，重排后得到新编号的序列（图3）。
图3. 以最大绝对值点为起点的重排序列
3.3 步骤3：循环计数——逐组判断提取循环
基于重排后的序列，按“连续三点判断”规则逐组识别有效循环，具体过程如下：
1. 取重排后序列的前三个点（点1、点2、点3），计算ΔS₁₂和ΔS₂₃； 2. 判断：ΔS₁₂ ≥ ΔS₂₃（不满足有效循环条件），放弃该组，移动至下一组（点2、点3、点4）； 3. 计算ΔS₂₃和ΔS₃₄，判断：ΔS₂₃ ≤ ΔS₃₄（满足有效循环条件），提取点2→点3的循环，其ΔS=ΔS₂₃，Sₘ=(S₂+S₃)/2； 4. 删除点2和点3，剩余点1、点4、点5（原序列点4）连接形成新序列； 5. 对新序列取前三个点（点1、点4、点5），计算ΔS₁₄和ΔS₄₅，判断：ΔS₁₄ = ΔS₄₅（满足有效循环条件），提取点1→点4的循环； 6. 删除点1和点4，剩余点5无法形成三组连续点，计数过程暂停。
示例：通过上述过程，从原始序列中提取出点2→点3和点1→点4两个有效循环（图4）。
图4. 循环提取后的剩余序列
3.4 步骤4：收尾处理——残余点处理
计数结束后，序列中可能剩余1个或2个点（无法形成三组连续点），需根据工程需求处理：
• 剩余1个点：无意义，直接舍弃； • 剩余2个点：若为“峰-谷”或“谷-峰”组合，可视为一个“半循环”，工程中通常按1个有效循环计入（或根据载荷历程的对称性调整）。
4. 典型示例解析
通过简单与复杂载荷历史两个示例，详细说明计数过程中的关键细节，特别是序列重排的重要性。
4.1 示例1：简单载荷历史（验证重排的必要性）
原始载荷历史转换为峰值-谷值序列后，若不进行重排，直接从原始顺序（点1→点2→点3→点4）开始计数：
• 第一组：点1、点2、点3，ΔS₁₂ > ΔS₂₃，不提取循环； • 第二组：点2、点3、点4，ΔS₂₃ ≤ ΔS₃₄，提取点2→点3循环； • 剩余点1、点4，无法形成三组连续点，计数结束，仅提取1个循环。
而重排后（以最大点4为起点）计数，可额外提取点1→点4的大应力范围循环（对损伤贡献更大），共提取2个循环，避免了关键循环的遗漏。
结论：序列重排是雨流计数的关键步骤，可确保最大应力范围循环被优先提取，保证评估结果的安全性。
4.2 示例2：复杂载荷历史（重排过程解析）
对于包含更多特征点的复杂载荷历史，重排逻辑一致，重点是准确识别最大绝对值点并完成序列重组：
1. 原始连续载荷历史（图5）转换为峰值-谷值序列（图6），得到点1至点11共11个特征点； 2. 遍历发现点11为绝对值最大的载荷点； 3. 重排规则：以点11为起点，跟随点11之后的原始顺序点（无，因点11为原始序列末尾点），再连接点11之前的原始顺序点（点1至点10）； 4. 重排后形成新序列并重新编号（图7），后续计数过程与简单示例一致，逐组判断提取循环。 
图5. 原始复杂连续载荷历史
图6. 提取峰值和谷值后的复杂序列
图7. 以最大绝对值点（点11）为起点的重排序列
5. 工程应用要点
雨流循环计数的结果直接影响疲劳寿命评估精度，工程应用中需关注以下关键问题：
5.1 数据采样率的影响
原始应力时间历史的采样率需足够高，确保不遗漏真实的峰值和谷值（采样率过低可能导致峰值被低估、谷值被高估，进而使应力范围计算偏小，评估结果偏于危险）。
建议：采样率需满足Nyquist定理，即采样频率≥2倍载荷的最高主导频率，一般工程中取5~10倍主导频率。
5.2 序列重排的唯一性
若存在多个绝对值相等的最大特征点（如两个峰值绝对值相同），可任选其一作为起点，重排后的计数结果一致（因最大应力范围循环的数量与位置固定）。
5.3 残余点的工程处理
计数结束后剩余的1~2个点，若为短时间历程（如10分钟内的载荷），可舍弃；若为长时间历程（如1小时以上的载荷），剩余点通常构成“半循环”，工程中按0.5个循环计入，或结合载荷的对称性补充为完整循环。
5.4 与其他计数方法的对比
工程中还有“单峰计数法”“跨零计数法”等循环计数方法，但雨流计数法的优势显著：
6. 结语
雨流循环计数是随机疲劳分析的“数据处理器”，其核心价值在于将无序的随机载荷转化为有序的等效恒幅循环，为疲劳损伤计算提供可靠输入。掌握“预处理→重排→计数”的标准化流程，理解“连续三点应力范围对比”的核心规则，是确保计数结果准确的关键。
工程应用中，需结合载荷特性合理设置采样率，重视序列重排步骤以避免关键循环遗漏，并根据标准规范处理残余点，才能使计数结果真实反映载荷的疲劳损伤特性，为结构的可靠性设计与寿命评估提供有力支撑。
