一、疲劳分析概述：到底什么是疲劳？

各位好，我是老张。做风电结构分析十几年了，今天咱们聊聊疲劳分析。

很多人刚接触风电时，第一反应是：「风机不是算强度吗？算完静强度不就够了？」

嗯，这个问题我当年也问过。直到有一次，我看到一个塔筒在运行三年后，焊缝处出现了肉眼可见的裂纹——而它的静强度校核是合格的。那一刻我才真正明白：强度够，不代表寿命够。

1.1 什么是疲劳？

说白了，疲劳就是材料在反复加载下，慢慢「累坏了」。

你想想看，一根铁丝，你用力拉一次，它可能不断。但你来回弯折几十次，它就断了。这就是疲劳。

风机的情况更复杂。叶片每转一圈，塔筒就受一次风载；阵风来了，载荷更大；停机启动，又是一次冲击。这些载荷反复作用，日积月累，材料内部就会产生微裂纹，慢慢扩展，最终断裂。

1.2 为什么风机必须分析疲劳？

原因很简单：风机设计寿命通常是20年，甚至25年。在这期间，叶片要旋转几千万到上亿次。

我做过一个项目，业主问：「能不能把塔筒壁厚减薄一点？省点成本。」我说可以，但疲劳寿命会从20年降到12年。业主沉默了。后来还是按原设计做了。

所以，疲劳分析不是「要不要做」的问题，而是「怎么做才能保证20年不出事」的问题。

1.3 三个核心概念

要搞懂疲劳分析，有三个概念必须吃透。我一个个讲。

1.3.1 S-N曲线

S-N曲线，全称是应力-寿命曲线。它描述的是：一个材料在某个应力水平下，能承受多少次循环。

横轴是循环次数N（对数坐标），纵轴是应力幅值S（线性或对数坐标）。曲线越往下，说明材料越「抗疲劳」。

举个例子：某钢材的S-N曲线显示，在200MPa应力下，能承受10^6次循环；在100MPa下，能承受10^8次循环。应力越低，寿命越长。

应力幅值 (MPa)	循环次数 N	说明
250	5×10^5	高应力，短寿命
200	1×10^6	典型设计点
150	5×10^6	中等应力
100	1×10^8	低应力，长寿命

1.3.2 Miner线性累积损伤

实际工况中，风机承受的载荷不是恒定的。今天刮大风，明天小风，后天停机。每个载荷水平都会产生一定的损伤。

Miner法则说：总损伤等于各载荷水平损伤之和。当总损伤达到1时，结构就失效了。

公式很简单：

D = Σ (ni / Ni)

其中：
D  = 总损伤
ni = 第i个应力水平下的实际循环次数
Ni = 第i个应力水平下允许的最大循环次数（从S-N曲线查得）

举个例子：假设某焊缝在200MPa下允许10^6次循环，在150MPa下允许5×10^6次循环。实际运行中，200MPa出现了20万次，150MPa出现了100万次。那么：

D = 200000/1000000 + 1000000/5000000
  = 0.2 + 0.2
  = 0.4

损伤0.4，说明还有60%的寿命余量。当D=1时，就该检修或更换了。

1.3.3 雨流计数法

这个名词听起来挺玄乎，其实它就是把乱七八糟的载荷时间历程，拆解成一个个完整的应力循环。

风机实际测到的载荷数据，是像心电图一样上下波动的曲线。雨流法就像「数雨滴」一样，把每个波峰波谷配对，统计出每个循环的幅值和均值。

为什么叫「雨流」？因为算法模拟了雨水从屋顶流下的路径——从波峰开始，沿着波形往下流，遇到更低的波谷就「滴落」，形成一个循环。

代码实现其实不复杂，我贴一个简单的Python示例：

def rainflow_counting(data):
    """
    简单的雨流计数法实现
    data: 载荷时间序列（列表）
    返回: 循环列表，每个元素为 (幅值, 均值)
    """
    cycles = []
    # 提取波峰波谷
    peaks_valleys = []
    for i in range(1, len(data)-1):
        if (data[i] >= data[i-1] and data[i] >= data[i+1]) or \
           (data[i] <= data[i-1] and data[i] <= data[i+1]):
            peaks_valleys.append(data[i])
    
    # 雨流计数核心逻辑（简化版）
    # 实际工程中建议使用现成库，如 rainflow 包
    # 这里仅展示思路
    stack = []
    for pv in peaks_valleys:
        stack.append(pv)
        while len(stack) >= 3:
            # 判断是否形成循环
            if abs(stack[-1] - stack[-2]) <= abs(stack[-2] - stack[-3]):
                cycle_range = abs(stack[-2] - stack[-1])
                cycle_mean = (stack[-2] + stack[-1]) / 2
                cycles.append((cycle_range/2, cycle_mean))
                # 移除已配对的点
                stack.pop(-2)
            else:
                break
    return cycles

1.4 知识体系总览

为了让你更直观地理解这三个概念的关系，我画了一张图：

这张图把三个概念串起来了。你记住一句话：雨流法是「拆」、S-N曲线是「查」、Miner法则是「加」。三步走完，疲劳寿命就出来了。

1.5 本章小结

好了，这一章就讲这些。总结一下：

• 疲劳是反复加载导致的渐进式损伤，风机必须分析是因为它要转20年。
• S-N曲线告诉你材料在不同应力下能扛多少次。
• Miner法则把不同载荷的损伤加起来，超过1就失效。
• 雨流计数法把乱七八糟的载荷数据变成一个个完整的循环。

下一章，我会带你走一遍完整的疲劳分析流程，从载荷数据到寿命评估，每一步该怎么做。咱们到时候见。

---

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