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 谐波。若要求线性相位、不产生波形时移,应选用零相位滤波(如对信号正反各滤一次),避免峰值位置偏移影响雨流计数。判断依据是:滤波后时域波形毛刺明显减少,频谱中工频尖峰被压到本底附近,而关注频段幅值基本不变。
处理顺序与参数选择
五类处理并非随意排列。合理的顺序是:先去趋势,再剔除异常值,最后滤波去噪。原因是去趋势依赖全局拟合,若先滤波会把漂移与真实低频混在一起难以分离;异常值若不先剔除,会污染滤波器的状态,使尖峰能量扩散到邻近点。参数选择上,去趋势阶数一般取 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 | 工频尖峰压至本底 | 带宽过宽损伤邻近频段 |
| 零相位处理 | 正反各滤一次 | 峰谷位置无时移 | 单次滤波引入相位滞后 |
完成上述清洗后,建议做一次量化对比:统计处理前后信号的标准差、峰值因子,以及用相同雨流计数参数得到的循环总数和最大幅值。若循环总数下降主要来自小幅值噪声循环的消失,而大幅值循环基本保留,说明清洗是成功的;若大幅值循环也明显减少,则说明滤波或去趋势过度,应回到参数调整环节。只有通过这一对比,才能确认清洗后的应变-时间信号可以可靠地进入载荷谱构建与疲劳寿命预测流程。