载荷谱、雨流计数、ncode中雨流计数的实现

本文雨流计数原理解释部分来自：https://wenku.baidu.com/view/814883e8998fcc22bcd10d69.html

在进行疲劳分析时，通常使用电阻应变片来测量部件关键位置上的交变载荷。加载点所承受的应力或载荷实际上是一个连续随机变化过程。把载荷循环中峰值和谷值随时间变化的过程，简称“应力一时间历程”或“载荷一时间历程”。

循环计数法的编谱步骤：①测量实际的载荷历程。②对实际的载荷历程作离散化处理，只记取峰值和谷值并作适当的压缩，压缩后的载荷历程消去了一些小量值的或变化不大的载荷。③使用循环计数法计算峰值或幅值。④对计数结果进行统计，计算出各载荷区间内所具有的频次、频率、累积频次和累积频率等统计量。⑤绘制载荷峰值或幅值变化的频率直方图和累积频次曲线，确定设计用的典型载荷谱。

雨流计数法，简称雨流法。它根据载荷历程得到全部的载荷循环，分别计算出全循环的幅值，并根据这些幅值得到不同幅值区间内所具有的频次，绘制出频次直方图。
（1）计数原理如图4-9c所示，载荷历程形同一座高层建筑物，雨点依次由上向下流动，根据雨点向下流动的轨迹确定出载荷循环，并计算出每个循环的幅值大小。每个载荷循环如用于疲劳寿命计算时，就对应于一个应力循环。
（2）计数规则采用雨流法时应遵守4项规则：①重新安排载荷历程，以最高峰值或最低谷值为雨流的起点（视二者的绝对值哪一个更大而定）。②雨流依次从每个峰值或谷值的内侧往下流，在下一个峰值或谷值处落下，直到对面有一个比开始时的峰值更大或谷值更小的值时停止。③当雨流遇到来自上面屋顶流下的雨流时即行停止。④取出所有的全循环，并记录下各自的幅值和均值。

 

把图a的载荷历程顺时针旋转90°，得到图c。根据上述规则，对这个载荷历程进行一次雨流计数，共计8个雨流：①a→b→b'→d点，然后下落。②b→c→d点的对应处，由于d的谷值比b为低。③c→b，遇到来自上面的雨流abb'd。bc与c