Python 齿轮振动信号模型构建

---

前言

---

前面的文章已从齿轮动力学机理上揭示了齿轮振动特性，为更直观的表述齿轮振动信号特征，现从齿轮振动信号模型角度来进一步分析其振动特征，以便为后续齿轮的状态监测提供参考。

一、齿轮振动信号模型

由前面动力学分析可知，齿轮的振动信号可以表示成齿轮啮合频率及其倍频的正弦信号叠加而成，故障等其他激励可以看做是对载波信号的调频、调幅，由此可得齿轮振动信号数学形式：

式中，K为阶次、a为调幅分量、fm为啮合频率、b为调频分量、theta为相位。当只考虑一阶分量时，齿轮振动信号数学模型为：

根据上述数学模型，我们便可以仿真出齿轮在各类工况下的振动特性。

二、仿真代码

1.Demo

代码如下（示例）：

import numpy as np
import matplotlib.pyplot as plt
import FFT

w = 5
z = 30
fs = 1024
fsw = 5
time = 1
f = w * z
t = np.linspace(0, time - 1 / fs, int(time * fs))
x = (1 + 1 * np.sin(2 * np.pi * 20 * t)) * np.sin(2 * np.pi * f * t)
amp, fre = FFT(x, 1024, fs)

plt.subplot(2, 1, 1)
plt.plot(t, x)
plt.ylabel('Amplitude')
plt.xlabel('time')
plt.subplot(2, 1, 2)
plt.plot(fre, amp)
plt.ylabel('Amplitude')
plt.xlabel('Frequency')
plt.show()

2.模型输出

上图为齿轮健康状态下的振动信号时域和频域变化特征，在健康状态下，时域信号为平稳变化的正弦信号，频域特征为齿轮的啮合频率。

上图为齿轮某个故障下的振动信号时域和频域变化特征，在该故障下，时域信号存在调幅情况，频域中的特征频率除了啮合频率，还在啮合频率两侧出现了边带。

---

总结

上述只是对齿轮部件振动特点做了一个简单分析，实际中的情况远比模型复杂，分析过程也更加严谨，下图是比较常见的齿轮频谱示例图(来源网络侵删),后续将逐步拆解分析齿轮各故障特点。

PS：欢迎各位交流，后续有啥想实现的信号处理功能，请在下方评论区留言，或者关注公众号：不说话上代码