3、振动测量与信号分析：传感器选型、数据采集、时域与频域分析、FFT与模态参数识别
振动测量是振动控制的“眼睛”与“耳朵”。若传感器选型不当、数据采集链路存在噪声或分析算法误用，后续的模态分析、减振设计与状态评估都将失去可靠基础。本章围绕“传感器—采集—分析—参数识别”这一完整链路展开，重点解决工程中三个核心问题：如何选对传感器并正确安装？如何从原始时域信号中提取有效的频域特征？如何从频响函数中识别模态参数并指导结构修改？
3.1 传感器选型：从原理到工程决策
核心原理：振动传感器将机械运动（位移、速度、加速度）转换为电信号。工程中最常用的是压电式加速度计，其利用压电晶体的正压电效应——当晶体受到外力变形时，在极化面上产生与加速度成正比的电荷。压电式传感器具有频响宽（0.1Hz~50kHz）、动态范围大（可达±1000g）、无需外部电源（电荷模式）等优点，但低频响应受限且不能测量静态加速度。压阻式加速度计则基于MEMS硅微桥的压阻效应，可测0Hz静态加速度，适合低频和冲击测量，但温漂较大。电容式加速度计精度高、稳定性好，多用于惯性导航和低频精密测量。
选型判断流程：工程选型需按以下步骤逐项筛选：
- 确定测量对象：是振动还是相对振动？是轴承座、基础还是柔性结构？
- 估算频率范围：旋转机械通常关注1kHz以内，齿轮箱需到10kHz，超声应用需到50kHz以上。
- 估算幅值范围：根据经验公式或历史数据，预估最大加速度峰值，留出30%余量。
- 检查环境约束：温度、湿度、电磁干扰、安装空间、线缆长度等。
- 质量效应校验：传感器质量应小于被测结构等效质量的10%，否则会改变结构固有频率。
| 类型 | 测量量 | 频响范围 | 典型灵敏度 | 适用场景 | 局限性 |
|---|---|---|---|---|---|
| 压电式加速度计 | 加速度 | 0.1Hz~50kHz | 10~100 mV/g | 通用机械、结构、冲击 | 不能测静态，低频差 |
| 压阻式加速度计 | 加速度 | 0~5kHz | 0.5~5 mV/g | 冲击、低频、跌落测试 | 温漂大，灵敏度低 |
| 电容式加速度计 | 加速度 | 0~1kHz | 100~1000 mV/g | 低频精密、惯性导航 | 量程小，价格高 |
| 速度传感器（磁电式） | 速度 | 10Hz~2kHz | 20~100 mV/(mm/s) | 旋转机械轴振 | 体积大，低频受限 |
| 电涡流位移传感器 | 位移 | 0~10kHz | 8 V/mm | 轴相对振动、间隙 | 需被测面为导体 |
3.2 数据采集：从模拟到数字的关键链路
核心原理：数据采集系统将传感器输出的模拟信号转换为数字信号，其性能直接决定分析精度。采样定理（Nyquist）要求采样率至少为最高分析频率的2倍，工程上通常取2.56倍以留出抗混叠滤波器的过渡带。抗混叠滤波器是模拟低通滤波器，在采样前滤除高于奈奎斯特频率的成分，否则高频信号会折叠到低频段产生虚假峰值。ADC位数决定动态范围：16位ADC的理论动态范围为96dB，24位可达144dB，但实际有效位数（ENOB）通常比标称值低2~3位。
操作流程：正确设置采集参数是保证数据质量的前提：
// 典型模态测试采集参数设置示例
采样率 fs = 2.56 × f_max = 2.56 × 2000 = 5120 Hz
采样点数 N = 2^14 = 16384
频率分辨率 Δf = fs / N = 5120 / 16384 = 0.3125 Hz
采样时间 T = N / fs = 3.2 s
平均次数 = 20（线性平均）
窗函数 = Hanning（稳态激励）/ Force窗（锤击激励）
触发方式 = 自由触发（稳态）/ 力锤触发（冲击）通道一致性校验：多通道测量时，各通道间的幅值偏差应小于±1%，相位偏差小于±1°。可通过背靠背校准或并联同一信号源进行校验。对于模态测试，通道间的时间延迟会导致相位误差，进而影响振型识别精度，因此应使用同步采样模式。
3.3 时域与频域分析：从波形到频谱的特征提取
核心原理：时域分析直接观察振动波形随时间的变化，可提取峰值、峰峰值、RMS值、波形因子、峭度等统计特征。RMS值反映振动能量，适用于稳态振动评估；峭度对冲击敏感，适用于轴承早期故障检测。频域分析通过傅里叶变换将时域信号分解为不同频率成分的幅值和相位，功率谱密度（PSD）表示单位频率带宽内的功率分布，是随机振动分析的核心工具。
FFT关键参数与窗函数选择：FFT将N点时域样本变换为N/2条有效谱线。窗函数用于抑制频谱泄漏——当信号频率不是频率分辨率的整数倍时，能量会泄漏到相邻频带。Hanning窗适用于稳态正弦或窄带信号，幅值精度约±1.5dB；FlatTop窗幅值精度可达±0.01dB但主瓣宽；矩形窗适用于瞬态冲击信号；Exponential窗用于锤击测试中衰减缓慢的结构响应。
| 窗函数 | 主瓣宽度 | 旁瓣衰减 | 幅值精度 | 适用场景 |
|---|---|---|---|---|
| 矩形窗 | 2Δf | -13dB | ±0.1dB | 瞬态冲击、校准 |
| Hanning窗 | 4Δf | -31dB | ±1.5dB | 稳态正弦、窄带随机 |
| Hamming窗 | 4Δf | -43dB | ±1.8dB | 语音、窄带信号 |
| FlatTop窗 | 8Δf | -44dB | ±0.01dB | 幅值精确测量 |
| Exponential窗 | 可变 | 可变 | 取决于衰减率 | 锤击法模态测试 |
频响函数（FRF）与相干函数：FRF定义为输出响应与输入激励之比，H(f) = S_xy(f) / S_xx(f)，其中S_xy为互功率谱，S_xx为自功率谱。FRF包含幅值和相位信息，是模态参数识别的基础。相干函数γ² = |S_xy|² / (S_xx · S_yy) 表示输出中由输入引起的比例，取值0~1。相干值低于0.8时，说明存在噪声干扰、非线性或泄漏，测量结果不可信。
背景：某化工厂离心泵（转速2970rpm，叶片数5）轴承座垂直方向振动超标，RMS值从2.1mm/s升至7.8mm/s。
测量方案：选用ICP型加速度计（灵敏度100mV/g，量程±50g），磁座安装于轴承座顶部。采样率fs=5120Hz，分析频率2000Hz，Hanning窗，线性平均20次。
时域分析：波形显示明显的周期性冲击，冲击间隔约20ms（对应50Hz），峰值达4.2g，峭度从3.1升至8.7，提示存在滚动体故障。
频域分析：FFT频谱中，1倍频（49.5Hz）幅值0.8g，2倍频（99Hz）幅值0.3g，但出现明显的边带族：以49.5Hz为中心，间隔为3.2Hz的边带（对应保持架故障频率BPFO=3.2Hz）。
结论：诊断为轴承外圈故障。拆检确认外圈滚道存在剥落。更换轴承后，RMS降至1.8mm/s，峭度恢复至3.0。
复盘：若仅看RMS值只能判断振动超标，但无法定位故障源；时域峭度提示冲击特征，频域边带分析精确定位故障元件。三者结合形成完整诊断链。
3.4 模态参数识别：从FRF到固有频率、阻尼与振型
核心原理：模态参数识别是从实测FRF中提取结构的固有频率、阻尼比和振型。单自由度系统FRF的幅值峰值对应固有频率，半功率带宽法可估算阻尼比：ζ = (f₂ - f₁) / (2fₙ)，其中f₁、f₂为幅值下降3dB处的频率。多自由度系统需采用更复杂的识别算法，如峰值拾取法（PP）、最小二乘复指数法（LSCE）、PolyMAX等。
操作流程：典型锤击法模态测试步骤如下：
- 结构准备：弹性悬挂（软弹簧或橡皮绳）模拟自由边界，悬挂频率应低于第一阶固有频率的1/5。
- 激励点选择：避开节点位置，选择刚度较大处，避免双击。
- 响应点布置：按结构几何形状均匀布置，每阶模态至少需要3~5个测点。
- 数据采集：力锤激励，Exponential窗处理响应信号，Force窗处理力信号，平均5~10次。
- FRF计算：采用H₁估计器（假设噪声主要在输出端），检查相干函数。
- 参数提取：使用PolyMAX或LSCE算法，通过稳态图筛选物理极点。
- 振型动画与MAC校验：计算模态置信准则（MAC），MAC值应大于0.8（同阶）且小于0.3（异阶）。
| 识别方法 | 适用场景 | 优点 | 局限性 |
|---|---|---|---|
| 峰值拾取法（PP） | 小阻尼、模态稀疏 | 简单快速 | 频率分辨率受限，阻尼精度差 |
| 频域分解法（FDD） | 环境激励、运行模态 | 抗噪声强 | 无法识别重模态 |
| LSCE（时域） | 脉冲响应清晰 | 精度高 | 对噪声敏感，需定阶 |
| PolyMAX（频域） | 通用，含重模态 | 稳态图清晰，鲁棒性好 | 计算量大 |
总结复盘：振动测量与信号分析是一个环环相扣的工程链路。传感器选型决定了测量数据的物理意义和可信度；数据采集参数设置决定了频率分辨率和动态范围；时域与频域分析提供了从波形特征到频谱结构的全面视角；模态参数识别则将测试数据转化为结构动力学模型。工程实践中，应始终遵循“先校验、后测量、再分析”的原则，对每个环节进行质量检查：传感器校准证书是否在有效期内？采集通道是否通过一致性校验？FRF相干函数是否满足要求？MAC矩阵是否对角占优？只有每个环节都可靠，最终识别出的模态参数才能用于指导结构修改和减振设计。建议工程师建立标准化的测量流程检查表，将上述要点固化为操作规范，减少人为误差。