机械元动作单元多故障耦合的动力学建模与特性分析

doi:10.3969/j.issn.1004-132X.2026.09.007

摘要/Abstract

摘要：

针对机械元动作单元多故障并发时特征模糊导致诊断控制困难的问题，提出了一种基于多故障耦合建模分析故障特性的方法。以蜗杆转动元动作单元为例，提取关键故障，基于分段线性弹簧阻尼模型模拟支座松动，基于附加综合不对中激励力模拟联轴器不对中故障。融合滚动轴承动力学模型建立蜗杆动力学模型，集成上述故障模型，建立蜗杆转动元动作单元多故障耦合动力学模型。采用数值分析法求解蜗杆在正常、单一松动故障、单一不对中故障及耦合故障状态下的加速度响应并进行频谱分析。根据现有实验平台，采集并分析上述4种状态的实际加速度时域信号并进行频谱分析。仿真与实验结果对比验证了所建多故障耦合模型的准确性。

关键词: 元动作单元, 多故障耦合, 数值分析, 故障特征提取

Abstract:

To address the problem that feature ambiguity caused by concurrent multiple faults in mechanical meta-action units leads to difficulties in diagnosis and control， a method based on multi-fault coupling modeling was proposed to analyze fault characteristics. First， taking the worm rotation meta-action unit as an example， key faults were extracted. Looseness of the support was simulated based on a piecewise linear spring-damper model. The coupling misalignment fault was simulated based on an additional combined misalignment excitation force. Second， integrating a rolling bearing dynamic model， the worm dynamic model was established. The above fault models were integrated to establish a multi-fault coupling dynamic model for the worm rotation meta-action unit. Then， the acceleration responses of the worm under normal conditions， single looseness fault， single misalignment fault， and the coupled fault state were solved using numerical analysis method. Spectral analyses were performed on these responses. Finally， based on the existing experimental platform， actual acceleration time-domain signals under the above four states were collected and analyzed. Spectral analyses were also performed on the actual signals. A comparison between simulation and experimental results verifies the accuracy of the established multi-fault coupling model.

Key words: meta-action unit, multi-fault coupling, numerical analysis, fault feature extraction

中图分类号:

TH17

葛红玉, 孙佳瑞, 李艺, 杨满芝, 高扬. 机械元动作单元多故障耦合的动力学建模与特性分析[J]. 中国机械工程, 2026, 37(9): 2141-2152.

GE Hongyu, SUN Jiarui, LI Yi, YANG Manzhi, GAO Yang. Dynamic Modeling and Characteristic Analyses of Multi-Fault Coupling in Mechanical Meta-action Unit[J]. China Mechanical Engineering, 2026, 37(9): 2141-2152.

导出引用管理器 EndNote|Ris|BibTeX

链接本文: https://www.cmemo.org.cn/CN/10.3969/j.issn.1004-132X.2026.09.007

https://www.cmemo.org.cn/CN/Y2026/V37/I9/2141

图/表 21

图1 蜗杆转动元动作单元结构示意图

Fig.1 Structural diagram of worm rotation meta-action unit

图2 蜗杆转动元动作单元试验台

Fig.2 Worm rotation meta-action unit test bench

表 1 故障频率统计表

Tab.1 Fault frequency statistics table

故障事件	年度发生频次	发生频率/%	主要伴随事件（发生频次）	并发频率/%
滚动轴承支座松动故障	20	28.60	联轴器不对中（14次）	70.00
联轴器不对中故障	18	25.70	滚动轴承支座松动（14次）	77.80
蜗杆出现裂纹	6	8.60	联轴器不对中（2次）	33.30
滚动轴承滚道划伤	6	8.60	滚动轴承支座松动（4次）	66.70
电机调速失灵	4	5.70	联轴器梅花键弹性形变（3次）	75.00
联轴器梅花键弹性形变	4	5.70	电机调速失灵（3次）	75.00
滚动轴承滚动体破损	4	5.70	滚动轴承滚道划伤（1次）	25.00
轴承装配间隙过大	4	5.70	滚动轴承滚道划伤（1次）	25.00
平键定位孔磨损	3	4.30	联轴器梅花键弹性形变（2次）	66.70
蜗杆腐蚀	1	1.40

图3 蜗杆转动元动作单元故障原理示意图

Fig.3 Schematic diagram of fault principle in worm rotation meta-action unit

图4 滚动轴承赫兹接触模型简化图

Fig.4 Simplified diagram of rolling bearing Hertz contact model

表 2 仿真参数

Tab.2 Simulation parameters

参数	数值
蜗杆密度ρ/（kg·m $- 3$ ）	7850
弹性模量E/GPa	200
滚动轴承等效质量m_bl1，m_bl2/kg	0.045
基础支撑等效质量m_s/kg	20
联轴器质量m_c/kg	2
基础支撑等效刚度k₃/（N·m $- 1$ ）	$1 × 105$
轴承赫兹接触刚度k_cb/（N·m $- 1.5$ ）	$2 × 108$
轴承间隙r₀/m	$0.015$
轴承外径D/mm	35
轴承内径d/mm	15
滚子数目	8
松动螺栓等效刚度k₁/（N·m $- 1$ ）	$1 × 104$
未松动螺栓等效刚度k₂/（N·m $- 1$ ）	$5 × 104$
松动螺栓等效阻尼c₁/（N·s·m $- 1$ ）	$10$
未松动螺栓等效阻尼c₂/（N·s·m $- 1$ ）	$100$
基础等效阻尼c₃/（N·s·m $- 1$ ）	$500$
支座松动间隙δ/m	0.005

表 2 仿真参数

Tab.2 Simulation parameters

参数	数值
蜗杆密度ρ/（kg·m $- 3$ ）	7850
弹性模量E/GPa	200
滚动轴承等效质量m_bl1，m_bl2/kg	0.045
基础支撑等效质量m_s/kg	20
联轴器质量m_c/kg	2
基础支撑等效刚度k₃/（N·m $- 1$ ）	$1 × 105$
轴承赫兹接触刚度k_cb/（N·m $- 1.5$ ）	$2 × 108$
轴承间隙r₀/m	$0.015$
轴承外径D/mm	35
轴承内径d/mm	15
滚子数目	8
松动螺栓等效刚度k₁/（N·m $- 1$ ）	$1 × 104$
未松动螺栓等效刚度k₂/（N·m $- 1$ ）	$5 × 104$
松动螺栓等效阻尼c₁/（N·s·m $- 1$ ）	$10$
未松动螺栓等效阻尼c₂/（N·s·m $- 1$ ）	$100$
基础等效阻尼c₃/（N·s·m $- 1$ ）	$500$
支座松动间隙δ/m	0.005

表 3 梁单元参数

Tab.3 Beam element parameters

编号	1	2	3	4	5	6	7
长度/mm	65	20	15	15	15	15	15
直径/mm	40	20	25	20	20	20	20
编号	8	9	10	11	12	13	14
长度/mm	15	15	15	15	15	15	15
直径/mm	20	20	20	20	20	20	25

图5 临界转速坎贝尔图

Fig.5 Critical speed Campbell diagram

图6 不同转速下关键节点加速度响应幅值

Fig.6 Acceleration response amplitudes of key nodes at different speeds

图7 正常状态下时域图与频谱图

Fig.7 Time domain and frequency domain diagrams under normal conditions

图8 不对中故障下时域图与频谱图

Fig.8 Time domain and frequency domain diagrams under misalignment fault

图9 松动故障下时域图与频谱图

Fig.9 Time domain and frequency domain diagrams under looseness fault

图10 不对中与松动故障耦合下时域图与频谱图

Fig.10 Time domain and frequency domain diagrams under coupled fault （misalignment + looseness）

表4 不同故障下节点P1仿真加速度响应的特征频率

Tab.4 Characteristic frequencies of node P1 simulation acceleration response under different faults

故障类型	Y方向加速度基频幅值/（m·s $- 2$ ）	Z方向加速度2倍频幅值/（m·s $- 2$ ）	Y方向加速度基频幅值/（m·s $- 2$ ）	Z方向加速度2倍频幅值/（m·s $- 2$ ）	主要频率/Hz
仅不对中	0.83	0.37	0.80	0.35	27.3，54.7
仅松动	0.64	0.01	0.64	0.01	27.3
不对中+松动	0.42	1.21	0.46	1.23	27.3，54.7，110.3，137.7
正常状况	0.60	0	0.61	0	27.3

表4 不同故障下节点P1仿真加速度响应的特征频率

Tab.4 Characteristic frequencies of node P1 simulation acceleration response under different faults

故障类型	Y方向加速度基频幅值/（m·s $- 2$ ）	Z方向加速度2倍频幅值/（m·s $- 2$ ）	Y方向加速度基频幅值/（m·s $- 2$ ）	Z方向加速度2倍频幅值/（m·s $- 2$ ）	主要频率/Hz
仅不对中	0.83	0.37	0.80	0.35	27.3，54.7
仅松动	0.64	0.01	0.64	0.01	27.3
不对中+松动	0.42	1.21	0.46	1.23	27.3，54.7，110.3，137.7
正常状况	0.60	0	0.61	0	27.3

图11 联轴器不对中故障实验

Fig.11 Misalignment fault experiment of coupling

图12 松动故障实验

Fig.12 Looseness fault experiment

图13 正常状态下实验时域图与频谱图

Fig.13 Experimental time domain and frequency domain diagrams under normal conditions

图14 不对中故障下实验时域图与频谱图

Fig.14 Experimental time domain and frequency domain diagrams under misalignment fault

图15 松动故障下实验时域图与频谱图

Fig.15 Experimental time domain and frequency domain diagrams under looseness fault

图16 不对中与松动故障耦合下实验时域图与频谱图

Fig.16 Experimental time domain and frequency domain diagrams under coupled fault （misalignment + looseness）

表 5 不同故障下实验加速度响应的特征频率

Tab.5 Characteristic frequencies of experimental acceleration response under different faults

故障类型	Y方向加速度基频幅值/（m·s $- 2$ ）	Z方向加速度2倍频幅值/（m·s $- 2$ ）	Y方向加速度基频幅值/（m·s $- 2$ ）	Z方向加速度2倍频幅值/（m·s $- 2$ ）	主要频率/Hz
仅不对中	0.80	0.28	0.80	0.33	27，55，110
仅松动	0.59	0.01	0.58	0.01	28
不对中+ 松动	0.33	1.37	0.34	1.40	28，54，110，137
正常状态	0.53	0	0.54	0	28

表 5 不同故障下实验加速度响应的特征频率

Tab.5 Characteristic frequencies of experimental acceleration response under different faults

故障类型	Y方向加速度基频幅值/（m·s $- 2$ ）	Z方向加速度2倍频幅值/（m·s $- 2$ ）	Y方向加速度基频幅值/（m·s $- 2$ ）	Z方向加速度2倍频幅值/（m·s $- 2$ ）	主要频率/Hz
仅不对中	0.80	0.28	0.80	0.33	27，55，110
仅松动	0.59	0.01	0.58	0.01	28
不对中+ 松动	0.33	1.37	0.34	1.40	28，54，110，137
正常状态	0.53	0	0.54	0	28

参考文献 13

[1]	李健，张根保，冉琰，等. 基于自由度约束的元动作单元误差传递路径识别方法［J］. 中国机械工程， 2023， 34（2）： 135-141.
	LI Jian， ZHANG Genbao， RAN Yan， et al. Error Transfer Path Identification Method of Meta-action Unit Based on DOF Constraints［J］. China Mechanical Engineering， 2023， 34（2）： 135-141.
[2]	陈果. 含不平衡-碰摩-基础松动耦合故障的转子-滚动轴承系统非线性动力响应分析［J］. 振动与冲击， 2008， 27（9）： 100-104.
	CHEN Guo. Nonlinear Dynamic Response Analysis of Rotor-ball Bearing System Including Unbalance-rubbing-looseness Coupled Faults［J］. Journal of Vibration and Shock， 2008， 27（9）： 100-104.
[3]	甄满，田拥胜，孙涛，等. 具有不对中故障的双跨转子系统非线性动力学特性［J］. 机械工程学报， 2020， 56（16）： 109-117.
	ZHEN Man， TIAN Yongsheng， SUN Tao， et al. Nonlinear Dynamics of Two-span Rotor-bearing System with Flexible Coupling Misalignment［J］. Journal of Mechanical Engineering， 2020， 56（16）： 109-117.
[4]	张伦，胡蛟，杨翼. 直升机主减速器锥齿轮系动力学建模与故障响应分析［J］. 国防科技大学学报， 2025， 47（5）： 236-245.
	ZHANG Lun， HU Jiao， YANG Yi. Bevel Gear Dynamic Modeling and Fault Response Analysis of Helicopter Main Reducer［J］. Journal of National University of Defense Technology， 2025， 47（5）： 236-245.
[5]	殷帆丽. 双转子-机匣耦合系统碰摩故障动力学特性研究［D］. 沈阳：东北大学， 2017.
	YIN Fanli. Research on Dynamic Characteristics of Dual Rotor-casing Coupling System with Rubbing Fault［D］. Shenyang： Northeastern University， 2017.
[6]	ZHENG Yang， WU Ke， WANG Luchao， et al. Structural Damage Assessment and Failure Mode Analysis for Cross-fault Submarine Tunnels［J］. Engineering Failure Analysis， 2024， 157： 107853.
[7]	CHEN Guowang， TANG Qian， LI Heng， et al. Rigid-flexible Coupling Dynamic-assisted Imbalanced Fault Diagnosis for Helicopter Tail Transmission System［J］. Measurement， 2025， 242： 115934.
[8]	TANG Hesheng， REN Yan， XIANG Jiawei， et al. Numerical and Experimental Analysis of Rotor-bearing System for Axial Piston Pump with Misalignment-Rubbing Coupling Fault［J］. Journal of Sound and Vibration， 2023， 559： 117786.
[9]	ZHAO Zhiyao， LI Xuanhao， CUI Xiaoyu， et al. Fault Coupling Analysis and Reliability Assessment of Actuation System Based on Bond Graph Model［J］. Applied Sciences， 2023， 13（13）： 7462.
[10]	LI Junyuan， GAO Peng， TONG Bang， et al. Revealing the Mechanism of Pack Ceiling Failure Induced by Thermal Runaway in NCM Batteries： a Coupled Multiphase Fluid-structure Interaction Model for Electric Vehicles［J］. eTransportation， 2024， 20： 100335.
[11]	张立志，曹成军，张伟，等. 基于场路耦合法的10 kV干式空心串联电抗器匝间绝缘故障电气特征量研究［J］. 电力科学与技术学报， 2025， 40（2）： 265-275.
	ZHANG Lizhi， CAO Chengjun， ZHANG Wei， et al. Study on Electrical Features of Inter-turn Insulation Fault of 10 kV Dry-type Air-core Series Reactor by Field-circuit Coupling Method［J］. Journal of Electric Power Science and Technology， 2025， 40（2）： 265-275.
[12]	王继康，孙敏，陈建恩. 刚度特征对振子链能量传递非互易性的影响［J］. 振动与冲击， 2025， 44（16）： 101-107.
	WANG Jikang， SUN Min， CHEN Jian’en. Effect of Stiffness Characteristics on Non-reciprocity of Energy Transfer in Oscillator Chain［J］. Journal of Vibration and Shock， 2025， 44（16）： 101-107.
[13]	李雪萍. 数控转台基于元动作的运动可靠性分析及建模方法研究［D］. 重庆：重庆大学， 2022.
	LI Xueping. Research on Motion Reliability Analysis and Modeling Method of NC Turntable Based on Meta-action［D］. Chongqing： Chongqing University， 2022.
	下转第 2176 页）

[1]	葛红玉, 赵展, 郭安祥, 孙佳瑞. 基于SABO-VMD的数控机床元动作单元故障可诊断性评价[J]. 中国机械工程, 2025, 36(08): 1774-1783.
[2]	王治超1, 刘超1, 冉琰2, 陈一凡3, 蒋东翔1, 张根保2, 3. 考虑分布和传递模型的机械系统可靠性分配方法[J]. 中国机械工程, 2025, 36(06): 1143-1150.
[3]	黄襄茂. RV减速机摆线齿轮接触应力数值分析方法与试验[J]. 中国机械工程, 2023, 34(24): 3015-3023.
[4]	李健, 张根保, 冉琰, 王勇勤. 基于自由度约束的元动作单元误差传递路径识别方法[J]. 中国机械工程, 2023, 34(02): 135-141.
[5]	张亢, 田泽宇, 陈向民, 廖力达, 吴家腾. 多通道多分量分解方法在变转速工况齿轮故障特征提取中的应用[J]. 中国机械工程, 2022, 33(20): 2483-2491.
[6]	徐洁, 俞树荣, 严如奇, 丁雪兴, 王世鹏, 丁俊华. 柔性端面气膜密封流场分析及密封特性研究[J]. 中国机械工程, 2022, 33(06): 656-663.
[7]	唐贵基;李楠楠;王晓龙. 综合改进奇异谱分解和奇异值分解的齿轮故障特征提取方法[J]. 中国机械工程, 2020, 31(24): 2988-2996.
[8]	高富东;王德心;王海东;曹建平. 基于射流冲击作用的舰机适配性数值分析与优化[J]. 中国机械工程, 2020, 31(12): 1425-1436.
[9]	王志永;杜伟涛. 螺旋锥齿轮铣齿机主轴\\|轴承系统的动力学分析[J]. 中国机械工程, 2019, 30(18): 2211-2216,2223.
[10]	周宏明, 张衡, 张祥雷, 彭姣娇, 张翔翔. 钻头冷镦参数建模及成形工艺优化[J]. 中国机械工程, 2016, 27(21): 2962-2966.
[11]	孟宗, 胡猛, 谷伟明, 赵东方. 基于LMD多尺度熵和概率神经网络的滚动轴承故障诊断方法[J]. 中国机械工程, 2016, 27(04): 433-437.
[12]	张国涛, 尹延国. 基于叶顶曲率半径变化的变量叶片泵叶片-定子副润滑数值分析[J]. 中国机械工程, 2015, 26(18): 2481-2485,2490.
[13]	屈小章, 韩旭, 毕仁贵, 谭艳 . 基于Kriging近似模型的轨道主冷叶轮多目标遗传优化[J]. 中国机械工程, 2015, 26(14): 1938-1943.
[14]	刘玉杰, 赵建社, 干为民, 祈璐. 轻量化阶梯孔结构电解加工试验研究[J]. 中国机械工程, 2015, 26(11): 1429-1433,1515.
[15]	陈仁祥, 汤宝平, 苏祖强. Hilbert-Huang变换分析仪及其在机械故障特征提取中的应用[J]. 中国机械工程, 2013, 24(13): 1755-1760.