轴承缺陷频率
$\text{BPFO}=\frac{N_b f_s}{2}\left(1-\frac{B_d}{P_d}\cos\alpha\right)$
模拟轴承缺陷、不平衡和对中不良产生的振动信号。通过FFT频谱和包络分析体验机械状态监测原理,自动计算BPFO、BPFI、BSF和FTF轴承特征频率。
$\text{BPFO}=\frac{N_b f_s}{2}\left(1-\frac{B_d}{P_d}\cos\alpha\right)$
滚动轴承外圈存在缺陷时,滚动体周期性滚过缺陷点,产生冲击振动。其发生的特征频率称为外圈故障频率(BPFO),计算公式由轴承几何尺寸决定:
$$\text{BPFO}=\frac{N_b f_s}{2}\left(1-\frac{B_d}{P_d}\cos\alpha\right)$$其中,$N_b$是滚动体数量,$f_s$是轴的旋转频率(Hz),$B_d$是球径,$P_d$是节圆直径,$\alpha$是接触角。这个频率是诊断外圈剥落、裂纹等故障的关键指纹。
为了量化冲击振动的剧烈程度,工程上常用波峰因数作为早期诊断指标。它是峰值与均方根值(RMS)的比值:
$$C_f = \frac{X_{peak}}{X_{RMS}}$$其中,$X_{peak}$是振动信号的峰值,$X_{RMS}$是有效值。正常机械运行平稳,$C_f$通常在2~4之间。当轴承出现初期损伤时,周期性冲击会使峰值显著增大,而RMS变化不大,导致$C_f$升高到4~10,从而发出预警。
风力发电机组监测:风机主轴轴承位于近百米高空,更换成本极高。通过在线振动监测,分析BPFO、BPFI等频率成分,可以提前数月预警轴承故障,避免灾难性停机和大额损失。
高铁与轨道交通:列车轮对轴承的状态直接关系行车安全。在列车入库时或通过地面监测点采集振动数据,利用包络分析技术,可以有效发现轴承的早期疲勞剥落,保障高速运行的安全。
大型工业压缩机与泵:在石油化工、冶金等行业,这些设备是生产的核心。通过定期采集振动频谱,识别不平衡、不对中以及轴承缺陷,可以制定预测性维护计划,减少非计划停机。
数控机床主轴:机床主轴轴承的精度直接影响加工质量。通过监测其振动特征频率,可以在轴承精度尚未明显劣化前进行维护,确保加工零件的尺寸精度和表面光洁度。
使用本工具时,有几个容易产生误解的地方。首先,人们常认为“固有频率越低(质量越大或弹簧越软)越危险”,但这其实取决于具体情况。虽然低频振动确实更容易在日常中出现,但关键在于“激励力频率”与“固有频率”的接近程度。例如,高速旋转的风扇(激励频率100Hz)对固有频率2Hz的机箱几乎不产生影响。反之,发动机怠速振动(20Hz)与固有频率20Hz的安装部件,即使频率较高,也会发生共振。
其次,请注意本工具未考虑“阻尼”因素。实际结构必然存在阻尼(将振动能量转化为热能等从而消耗振动的效应)。阻尼较大时,共振峰处的响应幅值会大幅低于理论值。本模拟器的频率响应图呈现的是“无阻尼系统”的理想结果,因此在实际工作中,下一步需加入阻尼项重新评估模型。
最后是参数设置的陷阱。当为质量或弹簧常数设置多个“串联”或“并联”参数时,经常容易算错合成值。例如,若两个弹簧在质点间呈“串联”关系,则合成弹簧常数k需通过 $1/k = 1/k_1 + 1/k_2$ 计算。在工具中分别调整k1和k2时,若未牢记此关系,可能导致非预期的刚度分布。
这种多自由度振动分析的思想,实际上潜藏在许多看不见的工程领域底层。声学工程中处理的是空气振动——“声音”。扬声器的振膜或乐器的琴身虽属于连续体,但理解其振动模态(分割振动)时,本工具所学的“多自由度系统模态”有助于建立直观印象。特别是二阶及以上模态出现的“节点”(振动近乎为零的点),直接对应着吉他弦的泛音或扬声器锥盆的分割振动。
另一个领域是控制工程。机器人手臂或精密平台的位置控制中,结构体的固有振动会严重影响控制系统的稳定性。若控制器输出指令信号的频率接近机械结构的固有频率,将引发超指令的过大振动(振铃现象),导致精度下降。为防止此问题,在控制系统设计阶段会作为“基于模型的设计”的一环进行特征值分析,从而安全设定控制带宽。
此外,这也与材料力学·断裂力学相关联。振动导致的循环应力是“疲劳破坏”的主要原因。共振状态下,即使激励力很小也会产生巨大应力,轻易超出材料的耐久极限。将预测最大应力产生位置的“应力分析”与本工具所学的“模态分析”相结合,便可实现预测疲劳寿命的“振动疲劳分析”。
熟悉本工具后,可进一步探索“连续体”振动。例如梁、板等质量与刚度连续分布的物体。此时固有频率有无限多个,但基本思路相通。例如,悬臂梁的一阶固有频率 $f_1$ 可表示为 $$f_1 \approx \frac{1}{2\pi} \cdot \frac{1.875^2}{L^2} \sqrt{\frac{EI}{\rho A}}$$ 其中L为长度,EI为弯曲刚度,ρA为单位长度质量。可对比式中如何体现相当于弹簧常数k与质量m的参数,以加深理解。
在数学层面,将多自由度系统的运动方程改写为状态空间表达,能更好地衔接控制理论与数值仿真。此外,工程中常用的有限元法模态分析,可视为本工具计算向数千至数百万自由度的扩展。推荐的学习路径为:1. 建立多自由度系统运动方程;2. 矩阵特征值问题的数值解法(雅可比法、QR法等);3. 含阻尼的广义特征值问题;4. 有限元法基础。最重要的是,先通过本工具亲手尝试大量参数变化下的响应特性,这将是所有进阶学习的基石。