技术摘要：
本发明涉及一种用于加速试验的拖拉机零部件载荷谱快速压缩方法，基于S‑N曲线与实测应力载荷谱，融合短时傅里叶变换的时频分析方法，准确辨识提取高疲劳损伤部分，快速压缩形成用于加速试验的载荷谱。该方法能够综合考虑载荷幅值与频率对零部件材料的耦合影响，最大限  全部
背景技术：
关键零部件可靠性是影响拖拉机田间作业连续规范的重要因素，试验是对拖拉机 零部件设计与制造的重要评价手段。在完成拖拉机任意零部件设计后，进行全工况的疲劳 试验，可以有效验证零部件具有设计中规定的强度和使用寿命。基于实验室的加速测试可 通过试验台复现零部件实际作业的受载情况，同时可以对加速试验加载信号进行压缩，因 此被广泛应用。 用于加速试验的零部件载荷谱是实验室加速测试的基础，其核心是疲劳数据编辑 技术，可以保证在加载时间上短于实测载荷谱长度，同时依然可以保证造成大约一致的疲 劳损伤。机械领域零件载荷谱编辑方法大致分为两种，基于时域以及基于频域的疲劳数据 编辑方法，其基础均为传统的、以S-N曲线为核心的疲劳分析理论，完全忽略作业载荷的频 率信息，只考虑应力幅值对于研究对象疲劳寿命的影响。 但是，拖拉机等农机装备在田间工作时，受限于缺少足够的悬架系统，与汽车的受 载情况存在较大差异，作业时主要受到宽频率段、幅值时变等特点的随机非平稳载荷。针对 工作振动明显，载荷频率段宽以及幅值变化大的农业装备编制用于加速试验的零部件载荷 谱，应综合考虑载荷幅值与频率对材料寿命的耦合影响。因此，提出基于短时傅里叶变换的 拖拉机零部件信号压缩方法，充分考虑拖拉机田间作业载荷特性，从载荷幅值与频率的角 度出发，充分识别并严格保留原始应力信号中的高疲劳损伤部分，通过对比信号前后的疲 劳损伤保留量、压缩比、统计参数、幅值域以及频率域来确定最终的缩减信号。
技术实现要素：
针对现有技术中存在的缺陷，本发明的目的在于提供一种用于加速试验的基于短 时傅里叶变换的拖拉机零部件载荷谱快速压缩方法，该方法能够综合考虑载荷幅值与频率 对零部件材料的耦合影响，最大限度压缩原始载荷信号并且能够保证压缩信号与原始信号 在疲劳损伤量、统计参数(均值、均方根值以及峰度系数)、幅值域穿级计数以及频率域功率 谱密度等方面保持一致，实现相同的加载效果。 为达到以上目的，本发明采取的技术方案是： 一种用于加速试验的拖拉机零部件载荷谱快速压缩方法，包括以下步骤： 1)获取实测的拖拉机零部件原始应力载荷谱； 2)对原始应力载荷谱进行预处理，预处理包括：重采样、消除趋势项、滤波以及去 除奇异点等； 3)求解预处理后的实测应力载荷谱的梯度，获得实测应力载荷谱的功率密度曲 线，求解材料S-N曲线的梯度，获得功率密度-寿命曲线方程； 3 CN 111581715 A 说　明　书 2/4 页 4)对实测应力载荷谱的功率密度曲线进行基于短时傅里叶变换的时频分析，获取 某一时刻单位时间段内的累积循环次数，对功率密度-寿命曲线方程求反函数，获取材料疲 劳失效的极限循环次数； 5)在步骤4)的基础上，对某一时刻单位时间段内的累积循环次数和材料疲劳失效 的极限循环次数做比得到单位时间内的累积功率密度； 6)重复步骤5)计算得到多个单位时间内的累积功率密度，多个单位时间内的累积 功率密度形成累积功率密度分布，然后设置累积功率密度分布的功率密度截止阈值，根据 功率密度截止阈值识别累积功率密度分布的高累积功率密度部分以及低累积功率密度部 分； 7)将高累积功率密度部分以及低累积功率密度部分分别映射到实测应力载荷谱 的时域，将实测应力载荷谱中的低疲劳损伤部分进行删减，提取高疲劳损伤部分； 8)将步骤7)中提取的高疲劳损伤部分进行保持原有载荷时间序列的重连接，形成 编辑后的压缩应力载荷谱； 9)计算压缩应力载荷谱与原始应力载荷谱的疲劳损伤保留量、统计参数，其中疲 劳损伤保留量、统计参数的计算误差均小于10％，并且二者的幅值域穿级计数以及频率域 功率谱密度分布具有一致性，以验证压缩应力载荷谱的准确性，否则，重新设置功率密度截 止阈值，直到验证成功； 10)完成拖拉机零部件应力载荷谱的压缩工作。 进一步地，短时傅里叶变换窗函数包括：矩形窗、三角窗、汉宁窗、汉明窗以及高斯 窗。 进一步地，步骤6)中的功率密度截止阈值与压缩应力载荷谱的应力时间历程长度 以及疲劳损伤保留量有密切的关系。 进一步地，所述功率密度截止阈值的合理选择能够使编辑后的压缩应力载荷谱压 缩比达到最大，疲劳损伤保留量达到最大。 进一步地，步骤9)中所述统计参数包括：均值、均方根值以及峰度系数。 进一步地，所述步骤8)中的重连接要保证原有载荷时间序列，尽可能进行平稳过 渡连接。 本发明与现有技术相比，具有以下优点： 1)基于短时傅里叶变换的拖拉机零部件载荷谱快速压缩方法与只关注幅值影响 的基于损伤保留的编辑方法有根本区别，前者综合考虑载荷幅值与频率对材料的耦合影 响，能够保持原始载荷信号的所有载荷特性。 2)基于短时傅里叶变换的拖拉机零部件载荷谱快速压缩方法与基于损伤保留的 编辑方法相比，前者得到的压缩信号压缩比更大，能够保持在47.5％。 3)基于短时傅里叶变换的拖拉机零部件载荷谱快速压缩方法与基于损伤保留的 编辑方法相比，得到的压缩信号的损伤保留量更大，并且保持在99％左右。 4)基于短时傅里叶变换的拖拉机零部件载荷谱快速压缩方法与基于损伤保留的 编辑方法相比，前者的统计参数(均值、均方根值以及峰度系数)误差更小，远远小于10％。 5)基于短时傅里叶变换的拖拉机零部件载荷谱快速压缩方法与基于损伤保留的 编辑方法相比，前者的幅值域穿级计数结果与原始载荷更接近。 4 CN 111581715 A 说　明　书 3/4 页 6)基于短时傅里叶变换的拖拉机零部件载荷谱快速压缩方法与基于损伤保留的 编辑方法相比，前者的频率域功率谱密度分析结果与原始载荷更接近。 附图说明 本发明有如下附图： 图1是用于加速试验的拖拉机零部件载荷谱快速压缩方法流程框图。 图2是利用累积功率密度截止识别并提取高累积功率密度部分示意图，其中图2 (a)为实测数据的累积功率密度分布及阈值水平示意图，图2(b)为实测数据的累积功率密 度分布辨识及删减结果示意图。