### Python 实现雨流计数法 雨流计数法是一种用于疲劳分析的经典算法，主要用于评估材料在循环载荷下的寿命。尽管Python标准库并未直接提供雨 flow counting 的内置支持，但可以通过自定义函数来实现该算法。 以下是基于已有资料和专业知识构建的一个简单版本的雨流计数法实现： #### 雨流计数法的核心逻辑 雨流计数法的主要目标是从给定的时间序列中提取完整的应力循环，并计算每种幅值对应的次数。其核心步骤包括闭合环路检测以及剩余部分的处理[^1]。 ```python def rainflow(array, min_range=0): """ 计算时间序列中的雨流计数。 参数: array (list): 输入的时间序列数组。 min_range (float): 最小范围阈值，默认为0。 返回: cycles (list of tuples): 循环列表，每个元组表示(均值, 幅度, 次数)。 """ stack = [] cycles = [] for i in range(len(array)): while len(stack) >= 2: # 获取最近两个极值点 top = stack[-1] next_top = stack[-2] # 判断是否形成闭环 if abs(top - array[i]) <= abs(next_top - top): break # 移除中间顶点并记录循环 amplitude = abs((next_top - top) / 2) mean_value = (next_top + top) / 2 if amplitude >= min_range: cycles.append((mean_value, amplitude, 1)) # 更新堆栈 stack.pop() # 将当前点压入堆栈 stack.append(array[i]) # 处理剩余未匹配的部分 while len(stack) >= 3: top = stack[-1] next_top = stack[-2] third_top = stack[-3] amplitude = abs((third_top - next_top) / 2) mean_value = (third_top + next_top) / 2 if amplitude >= min_range: cycles.append((mean_value, amplitude, 1)) stack.pop() stack.pop() return cycles # 示例代码运行 if __name__ == "__main__": data = [1, 2, 1, 3, 0, 2, -1, 1] result = rainflow(data) print("Rainflow Counting Results:") for cycle in result: print(f"Mean: {cycle[0]}, Amplitude: {cycle[1]}, Counts: {cycle[2]}") ``` 上述代码实现了基本的雨流计数功能，输入是一个时间序列 `array`，返回的是一个包含 `(均值, 幅度, 次数)` 的列表。通过调整参数 `min_range` 可以过滤掉幅度较小的循环[^2]。 #### 关键概念解释 - **闭环条件**: 当前点与堆栈顶部两点构成的三角形满足特定几何关系时，则认为形成了一个闭环。 - **最小范围筛选**: 设置 `min_range` 来忽略那些不重要的微小波动，从而减少不必要的计算开销[^3]。 --- ###