使用Rainflow的等效时刻

Question

我正在尝试从负载系列中进行疲劳分析，我想提取等效时刻

number cycles=1000000
years=25

我有一个小时的时间序列，如下所示：

然后我读到Rainflow Analysis是一个非常好的工具，可以从时间历史中提取周期。因此我申请：

%Rainflow moment
dt=time(2)-time(1);
[timeSeriesSig, extt] = sig2ext(timeSeries, dt);
rf = rainflow(timeSeriesSig,extt);

我看了   OUTPUT     rf - 雨流循环：矩阵3xn或5xn取决于输入，       rf（1，:)周期振幅，       rf（2，:)周期平均值，       rf（3，:)周期数（0.5或1.0），       rf（4，:)开始时间（当输入包括dt或extt数据时），       rf（5，:)循环周期（当输入包括dt或extt数据时），

如果我对周期数感兴趣，那么术语rf（3，:)是什么意思？它在向量中仅包含0.5和1。我想获得一个直方图，其中包含每个bin振幅的周期数。感谢

Answer 1

如果您使用File Exchange或其他来源的代码，则有助于链接到您获取它的位置。

此处，rf(3,:)显示幅度和其他输出是指半周期还是全周期。雨流算法首先找到半周期，然后匹配相等幅度的拉伸和压缩半周期以找到完整周期。通常有一些剩余的半周期。

如果您检查同一文件交换提交中包含的rfhist函数，您将看到它们如何处理完整和半幅度。基本上，在创建直方图时，任何给定负载幅度下的两个半周期将被计为一个完整周期。