4. 实测载荷信号预处理

实测应变-时间信号在进入雨流计数和载荷谱编制之前,必须先完成一轮系统性的清洗。原因很直接:应变片采集到的电压信号里,真正代表结构受载的应变分量往往只占一部分,其余还叠加了温度漂移、工频干扰、电磁噪声、导线抖动、A/D 量化误差以及偶发的尖峰跳变。如果把这些成分直接送入后续统计,会凭空增加大量小循环,使雨流计数得到的循环数虚高、幅值分布失真,最终把疲劳寿命预测带向保守或非保守的错误方向。本节围绕去趋势、去噪、异常值剔除、平滑与滤波五类处理,说明它们各自解决什么问题、按什么顺序执行,以及如何判断处理是否合格。

先明确处理对象与参照约定

设采样率为 fs,单位为 Hz,采样间隔 Δt = 1/fs。以某旋转部件应变测试为例,fs = 1000 Hz,单次记录 60 s,共 60000 个采样点,应变以微应变 με 为单位,规定拉应变为正、压应变为负。原始信号记为 x[n],清洗后信号记为 y[n]。所有处理都在时域或频域对同一序列操作,不改变时间轴,也不改变正负号约定,这样后续雨流计数的峰值谷值才能与原始物理过程对应。

去趋势:消除缓慢漂移

去趋势针对的是远低于结构关注频段的缓慢变化,典型来源是温度变化引起的桥路零点漂移、放大器直流偏置漂移。它的表现是信号整体抬高或缓慢倾斜,叠加在真实循环之上。处理方法有两种:一是减去全局均值,适合只有恒定偏置的情况;二是最小二乘拟合一条低阶多项式(通常一阶直线或二阶曲线)并逐点减去,适合存在线性或缓变漂移的情况。以一阶去趋势为例,对 x[n] 拟合 a·n + b,得到 y[n] = x[n] − (a·n + b)。判断依据是:去趋势后信号的长期均值应接近零,首尾段基线不再单调抬升或下沉。需要注意,去趋势的阶数不能过高,否则会把真实的长周期载荷成分一并削掉。

异常值剔除:处理尖峰与丢点

异常值通常由电磁脉冲、接触不良、瞬时短路或采集丢点造成,表现为孤立的极大或极小跳变,宽度往往只有一两个采样点。常用判据是三点中值或滑动中值比较:若某点与邻域中值之差超过阈值 k·σ(σ 为局部标准差,k 常取 3~5),则判为异常。剔除后不能简单置零,而应做插值修复,例如用前后各若干点的线性插值或样条插值填补,保持时间轴连续。判断依据是:修复后局部波形与相邻正常段过渡自然,不再出现单点尖刺,且不引入新的阶跃。

去噪、平滑与滤波:抑制随机与周期干扰

去噪主要针对宽带随机噪声和工频及其谐波干扰。平滑是最简单的去噪手段,滑动平均对白噪声有效,但会衰减高频真实成分并展宽峰值;因此对需要保留峰谷幅值的疲劳信号,滑动平均要慎用,或只用于辅助观察。更合适的是低通滤波与带阻滤波组合:低通滤波器保留结构关注频段,截止频率 fc 应高于最高关注载荷频率而低于噪声主频;带阻滤波器专门压制 50 Hz 工频及其 100 Hz、150 Hz 谐波。若要求线性相位、不产生波形时移,应选用零相位滤波(如对信号正反各滤一次),避免峰值位置偏移影响雨流计数。判断依据是:滤波后时域波形毛刺明显减少,频谱中工频尖峰被压到本底附近,而关注频段幅值基本不变。

实测应变信号预处理链路与效果判据 原始信号 x[n] 应变片电压→με 含漂移/噪声/尖峰 ① 去趋势 减均值 / 低阶多项式 消除温度漂移 ② 异常值剔除 中值判据 k·σ 插值修复尖峰 ③ 滤波去噪 低通 + 工频带阻 零相位避免时移 清洗后信号 y[n] 送入雨流计数与载荷谱编制 处理效果评估:时域 + 频域 + 统计量三重判据 时域观察 • 基线均值≈0,无缓慢倾斜 • 无孤立尖刺与阶跃 • 峰谷位置未发生偏移 • 波形连续、过渡自然 示例:60 s / 1000 Hz / 60000 点 频域观察 • 50/100/150 Hz 尖峰被压制 • 关注频段幅值基本不变 • 高频本底噪声下降 • 无新引入的虚假频率 示例:fc 取关注频率上限的 1.2 倍 统计量对比 • 标准差下降但非归零 • 峰值因子趋于合理 • 雨流循环数不再虚高 • 幅值分布形态稳定 示例:处理前后循环数对比

处理顺序与参数选择

五类处理并非随意排列。合理的顺序是:先去趋势,再剔除异常值,最后滤波去噪。原因是去趋势依赖全局拟合,若先滤波会把漂移与真实低频混在一起难以分离;异常值若不先剔除,会污染滤波器的状态,使尖峰能量扩散到邻近点。参数选择上,去趋势阶数一般取 1 或 2;异常值阈值 k 取 3~5,过小会误删真实峰值,过大则漏检;低通截止频率 fc 应高于结构最高关注载荷频率,通常取该频率的 1.2~2 倍,同时低于工频;带阻中心频率取 50 Hz 及谐波,带宽取 2~5 Hz。所有参数都应记录在预处理日志中,便于复现和追溯。

判断处理是否合格,最实用的办法是保留一份原始信号副本,把处理前后波形叠在一起看:真实循环的峰谷应一一对应,只是毛刺和漂移消失;若峰谷数量或位置明显改变,说明滤波过强或去趋势阶数过高,需要回调参数。

一个可复用的预处理检查清单

处理环节关键参数合格判据常见风险
去趋势阶数 1~2基线均值≈0,无单调倾斜阶数过高削掉真实低频
异常值剔除k=3~5,邻域 5~11 点无孤立尖刺,插值过渡自然阈值过小误删真实峰值
低通滤波fc≈1.2~2 倍关注频率关注频段幅值不变fc 过低衰减真实循环
工频带阻50/100/150 Hz,带宽 2~5 Hz工频尖峰压至本底带宽过宽损伤邻近频段
零相位处理正反各滤一次峰谷位置无时移单次滤波引入相位滞后

完成上述清洗后,建议做一次量化对比:统计处理前后信号的标准差、峰值因子,以及用相同雨流计数参数得到的循环总数和最大幅值。若循环总数下降主要来自小幅值噪声循环的消失,而大幅值循环基本保留,说明清洗是成功的;若大幅值循环也明显减少,则说明滤波或去趋势过度,应回到参数调整环节。只有通过这一对比,才能确认清洗后的应变-时间信号可以可靠地进入载荷谱构建与疲劳寿命预测流程。