一种非平稳非高斯振动损伤等效试验谱制定方法技术

本发明专利技术提供了一种非平稳非高斯振动损伤等效试验谱制定方法，通过计算实测非平稳非高斯信号的峭度；再利用幅值调制法对实测非平稳非高斯信号进行建模；最终求解振动损伤等效试验谱。相较于现有技术，本发明专利技术可准确模拟非平稳非高斯振动损伤等效试验谱，提高振动试验精度，提高产品可靠性和疲劳寿命的评估精度，从而提高企业的经济效益。而提高企业的经济效益。而提高企业的经济效益。

【技术实现步骤摘要】
一种非平稳非高斯振动损伤等效试验谱制定方法


[0001]本专利技术涉及机械振动领域中非平稳非高斯信号的振动损伤等效试验谱的制定方法。

技术介绍

[0002]一直以来，随机振动台都是利用功率谱密度(PSD,Power Spectra l Dens ity)来模拟产品经历的真实的随机振动环境。事实上，只有当实测随机振动信号是平稳的且服从高斯分布时，PSD才能完全描述该类振动环境。研究表明，峭度大于3的实测非平稳非高斯信号会导致产品更快地发生疲劳失效。为了在实验室准确评估非平稳非高斯随机振动环境给产品带来的振动损伤，疲劳损伤谱(Fatigue Damage Spectrum,FDS)被广泛用于制定振动损伤等效试验谱。为了减小误差，非平稳非高斯信号的FDS需要在时域进行求解。然后根据损伤等效原理，由时域FDS反推得到用PSD来表征的振动损伤等效试验谱，从而评估产品的使用寿命和可靠性。
[0003]然而，现有技术是利用时域法直接计算非平稳非高斯信号的FDS谱，当数据量较多时，该计算方法效率较低，导致产品的疲劳寿命评估时间过长。此外，当利用公式(10)反推振动损伤等效试验谱时，如果输入PSD在每个单自由度的半功率带宽内为非平直谱，得到的振动损伤等效试验谱将会引入较大误差，从而导致产品的疲劳损伤评估不准确，产品质量无法得到保证。

技术实现思路

[0004]本专利技术旨在提供一种非平稳非高斯振动损伤等效试验谱制定方法。
[0005]本申请的一种非平稳非高斯振动损伤等效试验谱制定方法，包括如下步骤：r/>[0006](1)计算实测非平稳非高斯信号的峭度
[0007]对于一个实测非平稳非高斯信号y(t)，其高阶原点矩可用于计算峭度：
[0008][0009]其中，μ为均值，当均值为0时：
[0010][0011][0012]其中，σ为均方根值，K为峭度；
[0013](2)利用幅值调制法对实测非平稳非高斯信号进行建模；
[0014]幅值调制法是用一个高斯信号g(t)乘以一个缓慢变化的幅值调制信号w(t)来模拟一个非平稳非高斯信号y(t)：
[0015]y(t)＝g(t)
·
w(t) (14)
[0016]假设高斯与非高斯信号的均值为0，且高斯信号与幅值调制信号相互独立，则非平稳非高斯信号的峭度为：
[0017][0018]其中，K
w
为幅值调制信号w(t)的峭度，幅值调制信号w(t)可用Weibull分布来描述，其概率密度函数为：
[0019][0020]其中，k为形状参数，λ是尺度参数；
[0021]其高阶原点矩为：
[0022][0023]由公式(15)和(17)可知，非平稳非高斯信号的峭度为：
[0024][0025]由公式(13)和(18)可知，当计算出实测信号的峭度后，可直接求解出形状参数k；
[0026]对于实测非平稳非高斯信号，其幅值调制信号变化缓慢，在每一个特定的时间段内可认为是常数，对于任意一个时间段t
i
，假设幅值调制信号的幅值为c
i
，该时间段与整个时域信号长度的比值为r
i
，那么每一个单自由度系统的疲劳损伤可表示为：
[0027]D
y
(t
i
)＝r
i
c
ib
D
g
(t
i
) (19)
[0028]其中，D
g
(t
i
)为高斯信号g(t)对应的疲劳损伤；
[0029]总损伤可表示为：
[0030][0031]其中，m为时间段数，T为总时域信号长度；
[0032]为了确保y(t)和g(t)的均方根值相同，需要对幅值调制信号的均方根值进行归一化，因此得到y(t)和g(t)的疲劳损伤比值为：
[0033][0034](3)求解振动损伤等效试验谱；
[0035]非平稳非高斯信号的损伤等效PSD谱可表示为：
[0036][0037]其中，D
y
(f
n
)指的是非高斯信号激励下每个共振频率f
n
上的疲劳损伤；
[0038]和实测信号具有相同均方根值的高斯信号g(t)对应的PSD谱为:
[0039][0040]其中，D
g
(f
n
)指的是高斯信号激励下每个共振频率f
n
上的疲劳损伤；
[0041]根据公式(21)
‑
(23)，两个PSD谱的比值为：
[0042][0043]因此，计算原信号的PSD谱，然后将实测数据的PSD谱和公式(24)得到的比值相乘即可得到非平稳非高斯振动损伤等效试验谱。
[0044]有益效果：本专利技术可准确模拟非平稳非高斯振动损伤等效试验谱，提高振动试验精度，提高产品可靠性和疲劳寿命的评估精度，从而提高企业的经济效益。
附图说明
[0045]图1是FDS的计算流程图；
[0046]图2是现场实测非平稳非高斯信号图；
[0047]图3是非平稳非高斯信号的原PSD和振动损伤等效PSD试验谱图；
[0048]图4是FDS对比图。
具体实施方式
[0049]下面详细描述本专利技术的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本专利技术，而不能理解为对本专利技术的限制。
[0050]在本专利技术的描述中，需要理解的是，术语“中心”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”、“逆时针”、“顺时针”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本专利技术和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本专利技术的限制。在本专利技术的描述中，除非另有说明，“多个”的含义是两个或两个以上。
[0051]在本专利技术的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接，可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本专利技术中的具体含义。
[0052]实施例
[0053]为了便于相关技术人员了理解现有技术，本实施例首先将介绍现有技术存在的缺陷以及具体缘由，内容如下：
[0054]现有技术利用时域计数法得到非平稳非高斯信号的应力循环分布，接着利用不同应力水平的循环数，结合S
本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种非平稳非高斯振动损伤等效试验谱制定方法，其特征在于，包括如下步骤：(1)计算实测非平稳非高斯信号的峭度对于一个实测非平稳非高斯信号y(t)，其高阶原点矩可用于计算峭度：其中，μ为均值，当均值为0时：其中，μ为均值，当均值为0时：其中，σ为均方根值，K为峭度；(2)利用幅值调制法对实测非平稳非高斯信号进行建模；幅值调制法是用一个高斯信号g(t)乘以一个缓慢变化的幅值调制信号w(t)来模拟一个非平稳非高斯信号y(t)：y(t)＝g(t)
·
w(t) (14)假设高斯与非高斯信号的均值为0，且高斯信号与幅值调制信号相互独立，则非平稳非高斯信号的峭度为：其中，K
w
为幅值调制信号w(t)的峭度；幅值调制信号w(t)可用Weibul l分布来描述，其概率密度函数为：其中，k为形状参数，λ是尺度参数；其高阶原点矩为：由公式(15)和(17)可知，非平稳非高斯信号的峭度为：由公式(13)和(18)可知，当计算出实测信号的峭度后，可直接求解出形状参数k；对于实测非平稳非高斯信号，其幅值调制信号变化缓慢，在每一个特定的时间段内可认为是常数，对于任意一个时间段t
i
，假设幅值调制信号的幅值为c
i
，该时间段与整个时域
信号长度的比值为r
i
，那么每...

【专利技术属性】
技术研发人员：徐飞，张森，
申请(专利权)人：盐城工学院，
类型：发明
国别省市：