齿轮系统故障诊断与传递路径分析的Matlab实现

1. 齿轮系统故障诊断与传递路径分析概述

齿轮系统作为机械传动的核心部件，其运行状态直接影响整个设备的可靠性。当齿轮出现磨损、断齿或装配偏差时，振动信号会通过轴系、轴承座等路径传递到机壳表面。传递路径分析(Transfer Path Analysis, TPA)正是通过量化激励源-传递路径-响应点的贡献关系，来定位故障源的技术手段。

在工程实践中，我们常遇到这样的情况：设备振动超标，但仅凭频谱分析难以确定是齿轮本身问题还是传递路径上的其他部件异常。这时TPA就能发挥独特优势——它不仅能告诉你振动能量从哪来，还能量化每条路径的贡献度。比如去年我们处理过一台减速箱案例，表面测得的振动加速度达到7.8m/s²，通过TPA分解后发现齿轮啮合频率处的主要能量其实来自轴承座的共振放大，而非齿轮加工误差。

2. TPA核心原理与齿轮系统适配

2.1 基本数学模型解析

TPA的核心公式可表示为：

code复制P = Σ(Hi × Fi)   (i=1,2,...,n)

其中：

• P：目标点响应（如振动加速度、声压级）
• Hi：第i条路径的传递函数（频响函数）
• Fi：第i条路径的激励力

对于齿轮系统，传递路径通常包括：

1. 齿轮啮合力→轴→轴承→箱体
2. 齿轮动载荷→轴系扭振→联轴器
3. 齿轮偏心激励→结构共振

2.2 齿轮系统特殊考量

与一般结构不同，齿轮系统的传递函数具有时变特性。这是因为：

• 啮合刚度随啮合位置周期性变化
• 齿面磨损会导致传递函数漂移
• 多级齿轮箱存在路径耦合现象

我们在某风电齿轮箱项目中实测发现，二级行星轮系的传递函数在运行3000小时后，200-400Hz频段的幅值增大了15dB，这正是齿面点蚀的典型特征。

3. Matlab实现关键步骤

3.1 数据采集与预处理

matlab复制% 示例：振动信号同步采集
fs = 5120; % 采样频率需大于最高啮合频率的2.56倍
t = 0:1/fs:10; 
accel_bearing = acquireData('CH1'); % 轴承座振动
accel_gearbox = acquireData('CH2'); % 箱体表面振动
force_transducer = acquireData('CH3'); % 力传感