双周期参数系统的受迫振动响应

顾京君，童彤，黄迪山; Gu Jingjun; Tong Tong; Huang Dishan

引 en

双周期参数系统的受迫振动响应^*

顾京君¹
机构：
1. 南通振康机械有限公司,南通 226153
×
，童彤¹
机构：
1. 南通振康机械有限公司,南通 226153
×
，黄迪山²
机构：
2. 上海大学机电工程与自动化学院,上海 200444
×

1. 南通振康机械有限公司,南通 226153；
2. 上海大学机电工程与自动化学院,上海 200444；

Forced Response of Dual Periodic Parametric Vibration System^*

Gu Jingjun¹
Affiliation：
1. Nantong Zhenkang Machinery Co., Ltd, Nantong 226153 , China
×
，Tong Tong¹
Affiliation：
1. Nantong Zhenkang Machinery Co., Ltd, Nantong 226153 , China
×
，Huang Dishan²
Affiliation：
2. School of Mechatronic Engineering and Automation, Shanghai University, Shanghai 200444 , China
×

1. Nantong Zhenkang Machinery Co., Ltd, Nantong 226153 , China；
2. School of Mechatronic Engineering and Automation, Shanghai University, Shanghai 200444 , China；

通讯作者:

黄迪山,E-mail:hdishan@shu.edu.cn

中图分类号:O3

文献标识码:A

文章编号:1672-6553-2023-22(3)-026-008

DOI:10.6052/1672-6553-2023-104

全文
评论
参考文献
出版信息

参考文献 1

李润方,王建军.齿轮系统动力学:振动、冲击、噪声 [M].北京:科学出版社,1997.

查找原文

参考文献 2

王建军,韩勤锴,李其汉.参数振动系统响应的频谱成分及其分布规律 [J].力学学报,2010,42(3):535-540.WANG J J,HAN Q K,LI Q H.Spectral components and their distributions of the response for parametric vibration system [J].Chinese Journal of Theoretical and Applied Mechanics,2010,42(3):535-540.(in Chinese)

查找原文

参考文献 3

SINHA S C,WU D H,JUNEJA V,et al.Analysis of dynamic systems with periodically varying parameters via Chebyshev polynomials [J].Journal of Vibration and Acoustics,1993,115(1):96-102.

查找原文

参考文献 4

KOVACIC I,RAND R,MOHAMED SAH S.Ma-thieu's equation and its generalizations:overview of stability charts and their features [J].Applied Mechanics Reviews,2018,70(2):020802.

查找原文

参考文献 5

HUANG D S,SHAO H X.Computation method for forced vibration response of a multiple DOF parametric system [J].International Journal of Structural Stability and Dynamics,2020,20(11):2050126.

查找原文

参考文献 6

HUANG D S,HONG L,LIU C.Computational technique to free vibration response in a multi-degree of freedom parametric system [J].Mechanical Systems and Signal Processing,2020,142:106777.

查找原文

参考文献 7

GAO N,MEESAP C,WANG S Y,et al.Parametric vibrations and instabilities of an elliptical gear pair [J].Journal of Vibration and Control,2020,26(19/20):1721-1734.

查找原文

参考文献 8

WANG S Y,ZHAO X X,XIA Y,et al.Mechanical-electromagnetic coupling elastic vibration instability of symmetrical three-phase external rotor induction motor [J].Nonlinear Dynamics,2019,97(1):1-20.

查找原文

参考文献 9

GU J J,TONG T,HUANG D S,et al.Study on torsional vibration of a harmonic driver based on time-varying stiffness caused by manufacturing error [J].Journal of Vibroengineering,2021,23(3):619-631.

查找原文

参考文献 10

XU C,WANG Z Z,ZHANG H L,et al.Investigation on mode-coupling parametric vibrations and instability of spillway radial gates under hydrodynamic excitation [J].Applied Mathematical Modelling,2022,106:715-741.

查找原文

参考文献 11

DENG J.Numerical simulation of stability and responses of dynamic systems under parametric excitation [J].Applied Mathematical Modelling,2023,119:648-676.

查找原文

目录contents

摘要 Abstract
关键词 Keywords
1 动力学建模与二重三角级数逼近
1.1 双周期参数振动方程
1.2 双周期参数振动响应解形式
2 受迫振动响应求解与计算
2.1 谐波力矩作用下受迫振动
2.2 振动响应计算
3 双周期参数振动特征
4 逼近计算误差
5 小结
参考文献

摘要

谐波减速器在制造过程中,如果柔轮齿与刚轮齿引入了周节累积误差,则谐波减速器在动力传递中将出现双周期时变扭刚度波动,形成双周期参数振动问题.本文引入基于组合频率的二重三角级数,逼近双周期参数系统受迫振动,应用谐波平衡,得到不含时间变量的谐波系数递推式,形成三维矩阵代数方程；采用矩阵降维法,将三维矩阵代数方程转化为两维矩阵代数方程,实现响应谐波系数的求解.研究结果对双周期参数振动中响应预测问题研究和工程应用具有一定的理论指导意义.

Abstract

During the manufacturing process of harmonic reducers, if the accumulated pitch error is introduced among the flexible gear teeth and the rigid gear teeth, there will be dual period time-varying torsional stiffness fluctuations in the power transmission of harmonic reducers, resulting in dual period parameter vibration problems. In this paper, the double trigonometric series based on the combination frequency is introduced to approximate the forced vibration response of the dual periodic parameter system. By using harmonic balance operation, the recurrence formula of harmonic coefficient without time variable is obtained, and the three-dimensional matrix algebraic equation is formed. Using the matrix dimension reduction, the three-dimensional matrix algebraic equation is transformed into two-dimensional matrix algebraic equation, and the response harmonic coefficients have been solved. The research results have certain theoretical guidance significance for the research and engineering application in the forced response problem of dual periodic parameter vibration system.

关键词

双周期参数系统；受迫振动；二重三角级数；谐波减速器

Keywords

dual periodic parameter system ； forced vibration ； double trigonometric series ； harmonic reducer

引言
谐波减速器广泛应用于各种机器人关节，它由柔轮、刚轮、波发生器三大基本构件组成.波发生器镶套在柔轮内圈，柔轮齿与刚轮齿进行内啮合，柔轮齿数比刚轮齿数少，因此，它是一个少齿差传动机构.其中，波发生器由椭圆凸轮外套柔性滚动轴承组成，基本结构如图1所示.
图1 单刚轮谐波减速器基本结构
Fig.1 Basic structure of single rigid wheel harmonic reducer
在制造过程中，如果柔轮齿与刚轮齿引入了周节累积误差，则谐波减速器在动力传递中将出现双周期时变扭刚度波动，形成双周期参数振动问题.
关于单自由度参数系统振动，学者王建军等^[1，2]利用Sylvester理论和Fourier 级数给出了自由和受迫振动响应的级数解；Sinha^[3]基于雪比多夫多项式，推出逼近级数解；Ivana Kovacic^[4]回顾了几种典型参数激励下的稳定性区域求解.近几年，利用组合频率三角级数法^[5，6]应用于多自由度参数振动分析，确定了自由与受迫振动响应的解析逼近.而大多数文献研究内容在基于Floquet理论上讨论参数振动稳定性问题^[7-11].但是，对于双周期参数振动的响应预测问题，仍然需深入探讨.
本文针对电机、谐波减速器和惯量负载传动问题，建立双周期参数振动方程，引入基于组合频率的二重三角级数，对其响应进行级数逼近解研究.
1 动力学建模与二重三角级数逼近
在电机、谐波减速器和惯量负载传动中，考虑双周期时变扭刚度波动，建立双周期参数振动方程；对于受迫振动求解，由等效动力学模型，提出基于组合频率的二重三角级数逼近.
1.1 双周期参数振动方程
如图2所示的是一个谐波减速器的刚度曲线K₀（t），它可以认为在基础刚度K_m上叠加了两个频率不同的刚度波动，其数学表达为

K_{0} (t) = K_{m} (1 + β_{1} c o s ω_{1} t + β_{2} c o s ω_{2} t)

(1)

式中，ω₁和 ω₂ 均为刚度波动频率，又称参数频率; β₁和 β₂为调制指数.
由于两个扭刚度波动的周期不同，但它们相互接近，导致扭刚度波动曲线呈拍频状（图2）.
图2 双周期时变扭刚度曲线K₀（t）
Fig.2 Dual period time-varying torsional stiffness curve K₀ (t)
在谐波减速器传动中，电机、谐波减速器和惯量负载构成双惯量弹性负载系统，如图3所示.
图3 双惯量弹性负载系统
Fig.3 Elastic loading system with dual inertia
双惯量弹性负载系统的扭振动方程为

J_{1} {\ddot{φ}}_{1} + C_{0} ({\dot{φ}}_{1} - {\dot{φ}}_{2}) + K_{0} (t) (φ_{1} - φ_{2}) = T_{1}

(2a)

J_{2} {\ddot{φ}}_{2} - C_{0} ({\dot{φ}}_{1} - {\dot{φ}}_{2}) - K_{0} (t) (φ_{1} - φ_{2}) = - T_{2}

(2b)

在方程（2）中，形成的刚度矩阵将为半正定，因此，双惯量弹性负载系统存在一个零固有频率.
对等式（2a）乘以J₂/J₁减去（2b）的代数运算，得到不含刚体位移的振动方程.

\begin{matrix} J_{2} ({\ddot{φ}}_{1} - {\ddot{φ}}_{2}) + C_{0} (1 + \frac{J_{2}}{J_{1}}) ({\dot{φ}}_{1} - {\dot{φ}}_{2}) + \\ K_{0} (t) (1 + \frac{J_{2}}{J_{1}}) (φ_{1} - φ_{2}) = \frac{J_{2}}{J_{1}} T_{1} + T_{2} \end{matrix}

(3)

虑K₀（t）是式（1）所描述的双周期变刚度扭弹簧，在方程（3）中，设 $θ = φ_{1} - φ_{2} ， J = J_{2} ， C = (1 + \frac{J_{2}}{J_{1}}) C_{0} ， K = (1 + \frac{J_{2}}{J_{1}}) K_{m} ， T (t) = \frac{J_{2}}{J_{1}} T_{1} ， T_{2} = 0$ ，即在没有外界负载干扰下，刻画一个具有双周期时变扭刚度的双惯量弹性负载系统，又称双周期参数系统，其振动方程简化为

\begin{matrix} J \ddot{θ} + C \dot{θ} + K (1 + β_{1} c o s ω_{1} t + β_{2} c o s ω_{2} t) θ \\ = T (t) \end{matrix}

(4)

式中，θ 为两个惯量之间的扭转角，力矩T（t）= T₀+Tcosω_pt来自于电机驱动.其中，T₀是电机恒力矩; Tcosω_pt是电机力矩波动，ω_p为力矩波动频率，一般它与电机轴转频相同.
1.2 双周期参数振动响应解形式
改写参数振动方程（4）为以下形式

\begin{matrix} J \ddot{θ} + C \dot{θ} + K θ \\ = T (t) - K (β_{1} c o s ω_{1} t + β_{2} c o s ω_{2} t) θ \end{matrix}

(5)

根据等式（5）表达，双周期参数系统的等效动力学模型是一个双调制反馈控制系统，如图4所示.
在双调制反馈控制系统中，振动响应θ（t）同时被频率ω₁和ω₂ 所调制，其频率被裂解，在叠加以后，反馈至系统的输入端.通过连续交错地频率裂解和组合过程，在振动响应θ（t）中，存在一系列组合频率ω+mω₁+nω₂谐波分量（m=-∞，···，-1，0，1，···，∞，n=-∞，···，-1，0，1，······，∞）.因此，双周期参数系统振动响应 θ（t）可以用基于组合频率ω+mω₁+nω₂的二重三角级数加以逼近，即响应的数学表达式（6）.

\begin{matrix} θ (t) = \sum_{m = - \infty}^{\infty} \sum_{n = - \infty}^{\infty} [E_{m, n} e^{j (ω + m ω_{1} + n ω_{2}) t} + \\ F_{m, n} e^{- j (ω + m ω_{1} + n ω_{2}) t}] \end{matrix}

(6)

当外激励力T（t）=0时，对应着双周期参数系统的自由响应，式（6）中的ω为主振荡频率，即ω=ω_s；当力矩T（t）≠0，θ（t）对应着双周期参数系统的受迫振动.如果外激励力T（t）=Tcosω_pt时，则ω=ω_p.
在图4的系统输出中，由于振动响应能量的有限性，能量主要分布在频率ω附近，随着三角级数项m→∞或n→∞，则谐波系数E_m，n→0和F_m，n→0.因此，对于振动响应逼近，可以采用有限项二重三角级数计算替代无限项级数的逼近.
图4 双周期参数系统的等效动力学模型
Fig.4 Equivalent dynamic model of double periodic parameter system
2 受迫振动响应求解与计算
在双周期参数系统受迫振动的二重三角级数逼近中，采用矩阵降维法，实现对受迫振动的谐波系数E_m，n的求解，得到受迫振动响应解.
2.1 谐波力矩作用下受迫振动
若双周期参数系统受迫振动方程（4）为

\begin{matrix} J \ddot{θ} + C \dot{θ} + K (1 + β_{1} c o s ω_{1} t + β_{2} c o s ω_{2} t) θ \\ = T c o s ω_{p} t \end{matrix}

(7)

对于系统的受迫振动，其响应解形式为

\begin{matrix} θ (t) = \frac{1}{2} \sum_{m = - k}^{k} \sum_{n = - k}^{k} [E_{m} e^{j (ω_{p} + m ω_{1} + n ω_{2}) t} + \\ F_{m, n} e^{- j (ω_{p} + m ω_{1} + n ω_{2}) t}] \end{matrix}

(8)

将响应解形式（8）和欧拉公式代入方程（7），对方程两边作谐波平衡，从正复指数 $e^{j (m ω_{1} + n ω_{2}) t}$ 部分，得到关于不含时间变量的谐波系数E_m，n的递推式.
当m=0 和n=0时

\begin{matrix} \frac{1}{2} β_{1} K E_{- 1,0} + \frac{1}{2} β_{2} K E_{0, - 1} + [K - J ω_{p}^{2} + \\ j C ω_{p}] E_{0,0} + \frac{1}{2} β_{1} K E_{1,0} + \frac{1}{2} β_{2} K E_{0,1} = T \end{matrix}

(9)

一般情况下，谐波系数E_m，n的递推式为

\begin{matrix} \frac{1}{2} β_{1} K E_{(m - 1), n} + \frac{1}{2} β_{2} K E_{m, (n - 1)} + [K - \\ J {(ω_{p} + m ω_{1} + n ω_{2})}^{2} + j C (ω_{p} + m ω_{1} + \\ n ω_{2})] E_{m, n} + \frac{1}{2} β_{1} K E_{(m + 1), n} + \\ \frac{1}{2} β_{2} K E_{m, (n + 1)} = 0 \end{matrix}

(10)

其中m=-k，···，-1，0，1，···，k和n=-k，···，-1，0，1，···，k.
这样，（2k+1）×（2k+1）个E_m，n的递推式构成了代数方程（11）.

Z E = P

(11)

式中，Z为2k+1阶三维系数矩阵，E为待求的2k+1阶二维谐波系数矩阵，P为 2k+1阶二维力矩阵.在公式（11）中，考虑了三维系数矩阵Z的阶数m足够大时，即m→∞，谐波系数E_m_+1，_n→0，E_m，n₊₁→0，E_（_m_+1），-_n→0和E_-_m_，-（_n_+1）→0.
为了矩阵Z、E和P降维，对于式（10）中的常数项，记

\begin{matrix} t = \frac{1}{2} β_{1} K \\ r = \frac{1}{2} β_{2} K \\ {\tilde{ω}}_{m, n} = K - J {(ω_{p} + m ω_{1} + n ω_{2})}^{2} + \\ j C (ω_{p} + m ω_{1} + n ω_{2}) \end{matrix}

(12)

引入子矩阵

Ω_{h} = [\begin{matrix} {\tilde{ω}}_{- k h} & t \\ ⋱ & ⋱ & ⋱ \\ t & {\tilde{ω}}_{- 1 h} & t \\ t & {\tilde{ω}}_{0 h} & t \\ t & {\tilde{ω}}_{1 h} & t \\ ⋱ & ⋱ & ⋱ \\ t & {\tilde{ω}}_{k h} \end{matrix}]

(13)

R = [\begin{matrix} r \\ ⋱ \\ r \\ r \\ r \\ ⋱ \\ r \end{matrix}]

(14)

谐波系数子向量E_h

(15)

力矩子向量P₀

(16)

将式（12）代入谐波系数E_m，n递推式（9）和（10）; 取h=-k，···，-1，0，1，···，k，将子矩阵和子向量式（13）至（14）按下标从-k至k依次在平面上排列，由此，把三维代数方程（11）展成为二维代数方程（17）.

(17)

记为

\bar{Z} \bar{E} = \bar{P}

(18)

式中， $\bar{Z}$ 为（2k+1）² 阶二维系数矩阵， $\bar{E}$ 为（2k+1）² 阶谐波系数向量， $\bar{P}$ 为（2k+1）² 阶力矩向量.利用矩阵逆运算，直接从方程（18）可解得向量 $\bar{E}$ .
同理，从负复指数 $e^{- j (m ω_{1} + n ω_{2}) t}$ 部分，可以得到另一组谐波系数F_m，n.其中，谐波系数E_m，n与F_m，n互为共轭.
2.2 振动响应计算
在双惯量弹性系统中，外界激励力矩分两部分，一是恒定力矩T₀，这时ω_p=0; 另一个是力矩波动，力矩波动频率ω_p =ω₁（一般电机力矩波动为恒定力矩的5%）.在不考虑负载力矩情况下，对于双惯量弹性负载系统的受迫振动响应，将按谐波力矩和恒力矩两种情况进行讨论.
在双惯量弹性负载系统（7）中，设惯量J=1，阻尼系数C=2.64，平均刚度K=17410，总扭刚度K（t）曲线如图2所示; 参数频率一ω₁=2πf₁，f₁=10.125Hz，调制指数一β₁=0.06; 参数频率二ω₂=2πf₂，f₂=10Hz，调制指数二β₂=0.07 .
算例1：若驱动谐波力矩（电机驱动力矩波动）波幅T=0.5，激励频率ω_p=ω₁，计算双惯量弹性负载系统的受迫扭振动响应.
取级数项m=-11，···，-1，0，1，···，11，n=-11，···，-1，0，1，···，11.根据所给的动力学参数，由公式（18）计算谐波系数矩阵中元素E_m，n，在频率ω_p附近的部分谐波系数数据列于表1.
表1 算例1中部分谐波系数E_m，n计算值
Table1 Value of partial harmonic coefficient E_{m, n} in Example 1
于是，该双惯量弹性负载系统受迫扭振动响应 θ₁（t）表达为：
$\begin{matrix} θ_{1} (t) = \frac{1}{2} \sum_{m = - 11}^{11} \sum_{n = - 11}^{11} E_{m, n} e^{j (20.25 + 20.25 m + 20 n) π t} + \frac{1}{2} \sum_{m = - 11}^{11} \sum_{n = - 11}^{11} E_{m, n}^{*} e^{- j (20.25 + 20.25 m + 20 n) π t} = \\ \sum_{m = - 11}^{11} \sum_{n = - 11}^{11} R e (E_{m, n}) c o s [20.25 (m + 1) + 20 n] π t - \sum_{m = - 11}^{11} \sum_{n = - 11}^{11} I m (E_{m, n}) s i n [20.25 (m + 1) + 20 n] π t \end{matrix}$
设置时间起点为0，步长0.001秒，时间历程为100秒，根据上述双惯量弹性负载系统受迫扭振动的响应表达，计算该系统受迫振动响应时间历程θ₁（t）和频谱Θ₁（ω），结果如图5所示，其中，振动响应有效值0.0318m·rad.
系统在外界激励力矩作用下，从振动响应频谱Θ₁（ω）可以看到，在各主谐峰（基频ω_p的高阶谐频）的左右侧和直流分量的右侧，存在着密集型的边频族分布.
图5 在谐波力矩驱动下的双周期参数系统受迫扭振动响应
Fig.5 Forced torsional vibration under a harmonic torque
算例2：在恒力矩驱动下，计算双惯量弹性负载系统的受迫扭振动响应.
设恒力矩力幅T₀=10，这时驱动频率ω_p=0，根据所给动力学参数，从式（18）计算谐波系数矩阵中元素E_m，n，部分计算值见表2.
于是，受迫扭振动响应θ₂（t）表达为
$\begin{matrix} θ_{2} (t) = \frac{1}{2} \sum_{m = - 11}^{11} \sum_{n = - 11}^{11} [E_{m, n} e^{j (20.25 m + 20 n) π t} + \\ E_{m, n}^{*} e^{- j (20.25 m + 20 n) π t}] \\ = \sum_{m = - 11}^{11} \sum_{n = - 11}^{11} R e (E_{m, n}) c o s (20.25 m + 20 n) π t - \\ \sum_{m = - 11}^{11} \sum_{n = - 11}^{11} I m (E_{m, n}) s i n (20.25 m + 20 n) π t \end{matrix}$
设置时间起点为0，步长0.001秒，时间历程为100秒，根据上述响应表达，计算该双惯量弹性负载系统受迫扭振动响应，其中振动响应时间历程θ₂（t）、响应频谱Θ₂（ω）结果如图6所示，振动谱特征与例1中的类同.其中，振动响应有效值0.0613m·rad.
与线性系统不同，双周期参数系统在恒力矩驱动下，仍然产生稳态的扭振动响应，而且，恒力矩越大，扭振动响应越大.
表2 算例2中部分谐波系数E_m，n计算值
Table2 Value of partial harmonic coefficient E_{m, n} in Example 2
图6 恒力矩驱动下的扭振动响应
Fig.6 Torsional vibration response under a constant torque
3 双周期参数振动特征
无论在算例1还是在算例2，在受迫振动响应谱中，都存在着密集型的边频分布，形成双周期参数振动响应特有的边频族特征.以算例2为例，对受迫振动响应谱图作局部放大处理，如图7所示.
在直流分量附近，存在5个数量级较大的边频分量，它们分别是0、ω₁-ω₂、2（ω₁-ω₂）、3（ω₁-ω₂）和4（ω₁-ω₂）；在第一个谱族附近，存在10个数量级较大的边频分量，它们分别是5ω₁-4ω₂、4ω₁-3ω₂、3ω₁-2ω₂、2ω₁-ω₂、ω₁、ω₂、2ω₂-ω₁、3ω₂-2ω₁、4ω₂-3ω₁、5ω₂-4ω₁；在第二、三、四个谱族附近，同样存在这种密集型边频族现象.双周期参数系统的受迫振动响应频谱具有丰富的边频，形成组合频率特征的边频族，谐波分量具体数值详见表2，其中边频间隔 Δ=ω₁-ω₂为两个参数频率之差.
图7 恒力矩驱动下的扭振动频谱局部放大图
Fig.7 Local enlarged figure of torsional vibration spectrum under a constant torque
为了验证组合频率特征边频族的客观存在性，搭建谐波减速器试验台，对一批双刚轮谐波减速器试验样件，在试验台架上进行扭振实测.其中的一个试验样件，柔轮齿与刚轮齿都有较大的周节累积加工误差，扭振测试结果如图8所示.在扭振动频谱局部放大图中，组合特征边频谱现象明显.其中，f₁=34Hz，f₂=33.58Hz，边频间隔约为Δ≈0.42Hz.
图8 某双刚轮谐波减速器扭振谱的组合频率特征
Fig.8 Characteristic of combination frequency in torsional vibration spectrum of dual rigid wheel harmonic reducer
实验表明：对于周节累积加工误差谐波减速器扭振，采用双参数振动建模，振动理论谱与实测结果一致性好，反映了动力学建模的有效性.
4 逼近计算误差
为了评估三角级数逼近计算精度，根据公式（7），定义逼近误差为：

\begin{matrix} e (t) = (J \hat{\ddot{θ}} + C \hat{\dot{θ}}) / K + (1 + β_{1} c o s ω_{1} t + \\ β_{2} c o s ω_{2} t) \hat{θ} - T c o s ω_{p} t / K \end{matrix}

(19)

其中， $\hat{θ}$ （t）为响应逼近；e（t）为逼近误差，单位为角位移或位移.
力幅T大小影响逼近误差，所以采用逼近计算误差ε（t）考核计算精度

ε (t) = \frac{e (t)}{| x |_{m a x}}

(20)

算例1的逼近计算误差时间历程ε（t），如图9所示，逼近计算误差ε小于5e-13（参考Runge-Kutta法计算误差小于1e-3）.因此，用二重三角级数逼近双参数系统振动响应，其逼近计算精度足高.
当然，双周期参数振动响应的逼近计算误差与二重三角级数逼近的项数有关，当级数逼近项数越多，逼近计算误差越小，同时，计算机耗时也越多.
图9 算例1中逼近计算误差时间历程ε（t）（k=11）
Fig.9 Time history of approximation error in example 1 ε (t) (k=11)
5 小结
在谐波传动中，如果柔轮齿与刚轮齿同时引入了周节累积误差，则谐波减速器在动力传递中将出现双周期参数振动问题.
基于组合频率的二重三角级数，可以对双周期参数系统受迫振动响应进行逼近，得到不含时间变量的谐波系数递推式，形成三维矩阵代数方程.通过引入中间变量，矩阵重新排列，对三维矩阵方程进行降维计算，实现对双周期参数系统的受迫振动响应谐波系数的求解，算法行之有效，并且准确性好.
通过计算分析可知，双周期参数系统受迫振动响应具有密集型组合频率边频族的特征，其中边频间隔Δ=ω₁-ω₂，为两个参数频率之差.
参考文献
- [1] 李润方,王建军.齿轮系统动力学:振动、冲击、噪声 [M].北京:科学出版社,1997.
- [2] 王建军,韩勤锴,李其汉.参数振动系统响应的频谱成分及其分布规律 [J].力学学报,2010,42(3):535-540.WANG J J,HAN Q K,LI Q H.Spectral components and their distributions of the response for parametric vibration system [J].Chinese Journal of Theoretical and Applied Mechanics,2010,42(3):535-540.(in Chinese)
- [3] SINHA S C,WU D H,JUNEJA V,et al.Analysis of dynamic systems with periodically varying parameters via Chebyshev polynomials [J].Journal of Vibration and Acoustics,1993,115(1):96-102.
- [4] KOVACIC I,RAND R,MOHAMED SAH S.Ma-thieu's equation and its generalizations:overview of stability charts and their features [J].Applied Mechanics Reviews,2018,70(2):020802.
- [5] HUANG D S,SHAO H X.Computation method for forced vibration response of a multiple DOF parametric system [J].International Journal of Structural Stability and Dynamics,2020,20(11):2050126.
- [6] HUANG D S,HONG L,LIU C.Computational technique to free vibration response in a multi-degree of freedom parametric system [J].Mechanical Systems and Signal Processing,2020,142:106777.
- [7] GAO N,MEESAP C,WANG S Y,et al.Parametric vibrations and instabilities of an elliptical gear pair [J].Journal of Vibration and Control,2020,26(19/20):1721-1734.
- [8] WANG S Y,ZHAO X X,XIA Y,et al.Mechanical-electromagnetic coupling elastic vibration instability of symmetrical three-phase external rotor induction motor [J].Nonlinear Dynamics,2019,97(1):1-20.
- [9] GU J J,TONG T,HUANG D S,et al.Study on torsional vibration of a harmonic driver based on time-varying stiffness caused by manufacturing error [J].Journal of Vibroengineering,2021,23(3):619-631.
- [10] XU C,WANG Z Z,ZHANG H L,et al.Investigation on mode-coupling parametric vibrations and instability of spillway radial gates under hydrodynamic excitation [J].Applied Mathematical Modelling,2022,106:715-741.
- [11] DENG J.Numerical simulation of stability and responses of dynamic systems under parametric excitation [J].Applied Mathematical Modelling,2023,119:648-676.

基本信息

中图分类号: O3
文献标识码: A
DOI: 10.6052/1672-6553-2023-104
文章编号: 1672-6553-2023-22(3)-026-008

基金信息

国家自然科学基金（51575330）；
National Natural Science Foundation of China （51575330）；

引用信息

稿件历史

收稿日期: 2023-07-08

参考文献

[1] 李润方,王建军.齿轮系统动力学:振动、冲击、噪声 [M].北京:科学出版社,1997.
[2] 王建军,韩勤锴,李其汉.参数振动系统响应的频谱成分及其分布规律 [J].力学学报,2010,42(3):535-540.WANG J J,HAN Q K,LI Q H.Spectral components and their distributions of the response for parametric vibration system [J].Chinese Journal of Theoretical and Applied Mechanics,2010,42(3):535-540.(in Chinese)
[3] SINHA S C,WU D H,JUNEJA V,et al.Analysis of dynamic systems with periodically varying parameters via Chebyshev polynomials [J].Journal of Vibration and Acoustics,1993,115(1):96-102.
[4] KOVACIC I,RAND R,MOHAMED SAH S.Ma-thieu's equation and its generalizations:overview of stability charts and their features [J].Applied Mechanics Reviews,2018,70(2):020802.
[5] HUANG D S,SHAO H X.Computation method for forced vibration response of a multiple DOF parametric system [J].International Journal of Structural Stability and Dynamics,2020,20(11):2050126.
[6] HUANG D S,HONG L,LIU C.Computational technique to free vibration response in a multi-degree of freedom parametric system [J].Mechanical Systems and Signal Processing,2020,142:106777.
[7] GAO N,MEESAP C,WANG S Y,et al.Parametric vibrations and instabilities of an elliptical gear pair [J].Journal of Vibration and Control,2020,26(19/20):1721-1734.
[8] WANG S Y,ZHAO X X,XIA Y,et al.Mechanical-electromagnetic coupling elastic vibration instability of symmetrical three-phase external rotor induction motor [J].Nonlinear Dynamics,2019,97(1):1-20.
[9] GU J J,TONG T,HUANG D S,et al.Study on torsional vibration of a harmonic driver based on time-varying stiffness caused by manufacturing error [J].Journal of Vibroengineering,2021,23(3):619-631.
[10] XU C,WANG Z Z,ZHANG H L,et al.Investigation on mode-coupling parametric vibrations and instability of spillway radial gates under hydrodynamic excitation [J].Applied Mathematical Modelling,2022,106:715-741.
[11] DENG J.Numerical simulation of stability and responses of dynamic systems under parametric excitation [J].Applied Mathematical Modelling,2023,119:648-676.

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。

双周期参数系统的受迫振动响应*

Forced Response of Dual Periodic Parametric Vibration System*

摘要

Abstract

关键词

Keywords

1 动力学建模与二重三角级数逼近

1.1 双周期参数振动方程

(1)

(2a)

(2b)

(3)

(4)

1.2 双周期参数振动响应解形式

(5)

(6)

2 受迫振动响应求解与计算

2.1 谐波力矩作用下受迫振动

(7)

(8)

(9)

(10)

(11)

(12)

(13)

(14)

(15)

(16)

(17)

(18)

2.2 振动响应计算

3 双周期参数振动特征

4 逼近计算误差

(19)

(20)

5 小结

参考文献

基本信息

基金信息

引用信息

稿件历史

参考文献

微信公众号二维码

手机版网站二维码

双周期参数系统的受迫振动响应^*

Forced Response of Dual Periodic Parametric Vibration System^*