基于分段间隙函数的螺旋锥齿轮时变啮合参数分析

本资料包含pdf文件1个，下载需要1积分

Time-varying meshing parameters analysis for spiral bevel gear based on sub-clearance functionLIU Zhi-feng，GUO Chun-hua，YANG Wen-tong，CAI Li-gang，ZHANG Zhi-min(Colege of Mechanical Engineering and Applied Electronics Technology，Beijing University of Technology，Beijing 100124，China)Abstract： Meshing parameters including stiffness and damping are key factors of nonlinear vibration of spiral bevelgears.Considering the time-varying characteristics of meshing stiffness and meshing damping，adopting a second-orderFourier series，vibration equation of a pair of spiral bevel gears was established.Based on the variation of clearancefunction under different meshing states，the effects of meshing stiffness and meshing damping on meshing features of a pairof spiral bevel gears under three meshing states were studied.Using Matlab，the vibration displacement curves of themeshing point were given when mean meshing damping，harmonic meshing damping，major harmonic meshing stiffnessand sub-harmonic meshing stiffness were considered. The results showed that increasing mean meshing damping anddecreasing major harmonic meshing stifness can drop vibration displacement peaks and improve transmission stability。

Key words： spiral bevel gear；time-varying meshing damping；time-varying meshing stiffness；clearance function；vibration displacement curves由于啮合阻尼、啮合刚度的动态变化和制造误差、安装误差、摩擦等原因而产生的间隙的存在，使得螺旋锥齿轮系统是-个具有参数激励的强非线性系统。目前国内外学者针对该动力学特性的研究已经作了大量的工作 ，建立了多种数学分析模型。Wang等 以-对准双曲面齿轮为研究对象建立了数学模型，将啮合参数按基频谐波展开，采用数值法求解，讨论了时变参数、系统参数及间隙非线性对齿轮副动力学特性的影响。Kahraman等 建立了直齿轮副的单自由度间隙非线性模型，将啮合刚度、间隙非线性函数、传动误差和外部激励展成多项谐波形式，用谐波平衡法求解。

基金项 目：国家 自然 科学 基 金 (51075006)；国 家科 技 重 大专 项(2010ZX04001-041)；北京市科委项目(Z11l104054111)收稿日期：2011-11-03 修改稿收到日期：2012-03-29第-作者 刘志峰 男，博士，副教授 ，1973年生通讯作者 郭春华 男，硕士，1985年生Litak等 考虑了轴的挠性，建立了三自由度齿轮模型，指出在二自由度模型上增加-个额外 自由度，齿轮副将显示不同的动态特性。

李润方等考虑了正弦型时变啮合阻尼，并假设与啮合刚度成正比，且相应于轮齿分离的速度也为非线性函数，指出非线性时变阻尼可采用常数阻尼来代替。唐进元等 对-类阻尼和刚度系数均含时变参数的强非线性齿轮系统进行了研究，针对时变阻尼项和刚度项之间的耦合作用使周期解的平均值发生漂移问题提出了-种改进的能量迭代法，并推导出系统主振动的幅频和相频响应方程。陈思雨等 考虑了时变刚度和摩擦时轮齿间隙及载荷参数对齿轮系统冲击动力学响应的影响，并首次研究 了非线性动力学的脱啮冲击。王建平等 建立了考虑动态啮合刚度、传递误差和齿侧间隙等多种非线性因素的直齿圆柱齿轮系统非线性动力学模型，将时变刚度按 5次谐波展开，齿侧间隙按 3次多项式拟合，并用多尺度法分析了齿轮系154 振 动 与 冲 击 2013年第 32卷统组合共振和谐振。张锁怀等 考虑了不平衡质量、主动轴转速波动和齿轮副齿侧间隙，并建立了单级齿轮传动系统的拍击振动模型，推导出了相邻两次碰撞之间的映射关系式。申永军等 M 建立了考虑间隙、时变刚度和静态传递误差的直齿轮副动力学模型，研究了阻尼 比和激励水平对幅频曲线的影响。刘国华等 提出了考虑轮齿弹性振动和单双齿啮合区变化的齿轮机构多体弹性非线性动力学模型。

上述研究主要利用解析、数值和实验的方法，分析齿轮时变刚度、齿侧间隙、静态传递误差对齿轮系统的参数激励、谐波共振、冲击等动力学响应，并未涉及由于间隙函数的分段变化，而产生的时变啮合参数在不同的啮合状态对螺旋锥齿轮啮合特性的影响。因此，本文根据啮合点法向位移与间隙不同关系下间隙函数的变化，考虑啮合阻尼和啮合刚度的时变性，研究平均啮合阻尼、谐波啮合阻尼、谐波啮合刚度在不同啮合状态下对螺旋锥齿轮副传动过程 中啮合特性的影响规律。

1 数学分析模型根据螺旋锥齿轮的传动特点，建立了-对齿轮副图 1 螺旋锥齿轮副啮合模型Fig.1 Meshing model of spiral bevel gear传动模型，如图1所示。假定传动轴和支撑部件均为刚性体，则主、从动齿轮的扭转平衡方程分别为：R C(t)(R 6P-R )R K(t)厂(R -R ) (1)Ig 0 -R C(t)(Rp -R )-R K(t)R -R )-Tg (2)令 R ， ，，扎 ，m Ig，乏惫贝0式(1)，式(2)变为：C(t)( - )K( ) - )F(3)- C(t)( -x )-K-厂( - )-F (4)令 将式(3)与式(4)联立解耦，得到降阶的螺旋锥齿轮副振动平衡方程：C(t)XK( )，(X)F (5)厂∞ . 1/-/b bl厂∞ l/-/。 1(a)无量纲化前 (b)无量纲化后图 2 间隙函数Fig.2 Clearance function2 时变啮合参数、间隙函数的无量纲化螺旋锥齿轮啮合过程中，会出现单齿对啮合及多齿对啮合交替变化的情况，使得啮合刚度始终处在不断的变化中。相对于单齿啮合区，齿轮的啮合综合刚度较小，啮合弹性变形较大；在双齿对啮合区，由于是两对轮齿同时承受载荷，因此齿轮的啮合综合刚度较大，啮合弹性变形较校随着齿轮副的连续运转，轮齿的弹性变形呈周期性变化，引起齿轮副角速度周期性变化，进而导致齿轮副振动。此外在考虑摩擦时，应将常数啮合阻尼描述为非线性时变啮合阻尼，即阻尼项也含时变参数～阻尼和刚度按傅里叶级数展开，并取前两项，可得：C(t)cmc1COS( t 。1)c2COS(20t c2)1K(t)kmklOS( t 1)kzCOS(2/2t 鸵)J(6)令 ， F 丢，力褰C丽m ， C丽1， C 2-k，k1 k2K -k-m l， p -k-m， -k-m (7)表 1 轮齿不同啮合状态下间隙函数变化Tab.1 Variation of dearltnce functionunder diferent meshing state- bOL>卢， ，c>P >将式(6)～式(8)代人式(5)得2[ OCOS( 。1)flcos(2/2t c2)] [Kpcos(t 1)vcos(2/'2t 娩)] )F各物理量意义见表2。

表 2 物理量对照表Tab.2 Table of physical symbol表3 螺旋锥齿轮传动系统参数Tab.3 Transmission system parameters of spiral bevel gear参数 数值齿数比 z /轴交角∑/(。)基本齿廓平均齿形角 (。)齿宽 b/mm参考点螺旋角卢 /(。)参考点法向模数 m /mm刀齿组数齿侧间隙/mm工作输入扭矩 /(N·m)工作转速 (r·min )10/419O2030353130.151 000-2 0002500-30003 齿轮系统振动因素分析令 1： 以 垃O，取 F0.8，下面分别从- 1

3.1 -1< <1情况当-1< <1时，齿轮副扭转方程表示为2[ acos(t)cos(2[2t)] F图3所示为齿轮系统振动的位移 -时间曲线。由图可知：随平均阻尼系数 的增大，系统趋向收敛。以t150为例，当 0.03时，啮合点法向位移为1 765，0.15时，啮合点法向位移为 389。平均阻尼系数的增大，使啮合点法向位移振动幅值有了明显的衰减，齿轮啮合过程中的振动迅速降低，提高了传动平稳性。

图4所示阻尼含有主谐波激励时，啮合点法向位移均方根值变化情况。以平均阻尼系数 0.15为例，当主谐波阻尼系数 OL0时，啮合点法向位移均方根值为148，为-条水平的光滑直线；O/0.O1时，法向位移均方根最小值为 173.34，最大值为294.61；当 0.10时，位移均方根最小值为282.11，最大值为400.76。由于加入了主谐波激励项，啮合点法向位移在整个频率范围内已不再是-条光滑的恒值曲线，而是在多个频率处出现波动，且随频率的增加，其值呈现上升趋势，并且主谐波系数的增加，会使位移均方根值迅速增大。

图5为阻尼含有次谐波项时，啮合点法向位移均方根值变化曲线。以 0.15，O/0.10为例，当卢0。

05时均方根值最小为285.76，最大值为447.51∩见当-1

3.2 -1情况当 -1时，齿轮副系统的扭转方程表示为：)/-2[ O/COS(T2t)flcos(2t)] [，cpcos(/)PCOS(2T2t)]( 1)F图6为 X≤ -1时，只考虑平均阻尼系数 情况下，齿轮副振动特性变化曲线。 0.03时，在 力1和 2处均出现了峰值，其值为3.95； O.05时，仅在1处出现了峰值，其值为2.64； 0.15时，峰值降为 1.18。随着 的增加，出现峰值的频率点由原来的两个减少到-个，且峰值随之减校这表明平均阻尼系数的增加有利于降低振动峰值，减小轮齿变形。

图7所示为X≤-1时，主谐波阻尼系数Ot对齿轮副系统振动影响曲线。以 0.1为例，随着 O/的变化，啮合点法向位移均方根值变化很小，在整个频率范围内最大误差为0.97%。图8所示为 -1时，次谐波阻尼系数 的变化对位移均方根值影响曲线。以 0.1， 0.05为例，在 1频率处，IB0时啮合点法向位移均方根值为1.54， 0.02时其值为1.47，最大误差为4.5%。当 ≤-1时，主谐波阻尼系数和次谐波阻尼系数的变化对啮合点法向位移均方根值几乎没有影响。因此，分析齿轮副振动时将啮合阻尼系数视为常数是合理的。

156 振 动 与 冲 击 2013年第 32卷图3 -1

图4 -1

O.5 1.0 1.5 2.0 2.5f图 10 -1， 影响振动位移曲线( ：0.1， 0.05，卢0.001)Fig.10 5 -1.vibrationdisplacement curves on主谐波刚度系数P对齿轮副动态特性影响如图9。

在平均刚度系数相同的情况下，P从 0增加到 0.05时，仍以线性特性为主；p继续增加到0.2，在 力1处首次出现峰值 1.54；当P0.5时，在 1附近和 2处均出现了峰值， 1附近的两个峰值分别为 3.47和3.64，力2处的峰值为 5.92。随着次谐波啮合刚度系数 的增加， 1处的峰值变化不是很明显，而 2O.5 1.0 1.5 2.O 2.5力 图 8 ≤-1， 影响振动位移曲线(K1，p0.2， 0.04)Fig.8 S -1.vibrationdisplacement CHives on- 0 03 -- O 05 · ：010 - 015 . 上 ≥-O.5 1.0 1.5 2.O 2.5力 图 11 ≥1， 影响振动位移曲线(K1，P0.2， 0.04)Fig.1 1 1.vibrationdisplacement curves on处的峰值逐渐增大。图 10所示，当/2由0增加到 o.20时， 2处的峰值也由5.92增加到 8.60。因此表明引入主谐波和次谐波刚度项时，均会对齿轮副振动特性产生影响。

3.3 X>l情况当 ≥1时，齿轮副系统的扭转方程可以表示为十2[ acos cos(2 五)] [Kpcos(/2t)PCOS(2 )](X-1)F- - -- - -r.... ，. .L... ....。LF，.。。。r.....L..。，。.r，，。。.r..L 0 2 0 8 6 4 2 0 8 6 4 5 5 4 4 4 4 4 3 3 3 l 第4期 刘志峰等：基于分段间隙函数的螺旋锥齿轮时变啮合参数分析 157由图 11可看出： 0.03时，力1处出现峰值，其值为 5.18； 0.05时， 1还是为峰值频率，但值降为4.47；随着 分别增加到0.10和0.15，不再存在峰值，力1处位移均方根值分别为4.12和4.09。

图 12所示为当 ≥1时，主谐波阻尼系数 对齿轮副系统振动影响曲线。以 0.15为例，由图可看出在 1.5频率处， 0时的啮合点法向位移均方根值为 4.26， 0.10时的啮合点法向位移均方根值为 4.06，最大误差为4.7%。当增加次谐波阻尼项时，如图13所示，以 0.15，O/0.1为例， 0.01时的啮合点法向位移均方根值为4.55， 0.03时的啮合点法向位移均方根值为4.81，最大误差为5.4%。由此可看出当 1时，阻尼含时变激励项对啮合点法向位移均方根值影响不明显。因此在这个阶段 ，非线性时变阻尼也可由常数阻尼来代替。

力 图 12 ≥1， 影响振动位移曲线(，c1，P0.2， 0.04)Fig.12 1.vibrationdisplacement curves on o/5·45.25.04.84·4·44·24.03.83.6图14所示为主谐波刚度系数P对齿轮副啮合特性位移影响曲线。不考虑谐波啮合刚度项，即P 0时， 2为啮合点法向位移均方根峰值频率，峰值为3.5。当仅将p分别增加到0.01和0.2时，位移均方根值在整个频率内均随频率 的增大而增大，峰值消失。

但当P增加到0.5时，在 1附近和 2处均出现了峰值，其值分别为 5.47、5.31和5.O7。

图15所示，当含有次谐波刚度项时，以K1，P0.5为例，随次谐波刚度系数 的增加，在 1附近和：2处的峰值均发生了变化。 0时，三个峰值分别为 5.47、5.31和5.07；/20.05时，三个峰值分别为 5。

09和5.44。即随着次谐波刚度系数的增加， 1附近的两个峰值在减小，而 2处的峰值却随着次谐波刚度系数的增加而增大。

力 图 13 1，口影响振动位移曲线(K1，P0.2， O.04)Fig.13 1.vibrationdisplacement curves onI 遭 6I Jl；j j ：J J工- - v0.00- 罟 - v0.05- -e- vO10力图 15 1， 影响振动位移曲线( O.1，d0.05，卢0.001)Fig.15 1.vibration displacement curves on4 结 论本文基于间隙函数的分段特性，考虑啮合阻尼和啮合刚度的时变性，并采用二次谐波展开的形式，讨论了平均啮合阻尼、谐波啮合阻尼和谐波啮合刚度在不同啮合状态下对螺旋锥齿轮副传动过程中动力学特性力 图 14 1，P影响振动位移曲线( O.1， O.05，JBO.001)Fig.14 1.vibrationdisplacement clLrYes on P的影响。确定了时变啮合参数下的螺旋锥齿轮副，在其不同啮合阶段时的啮合特性，为螺旋锥齿轮的减振和降噪提供了理论依据。

(1)-l

(2) ≤-1时，平均啮合阻尼系数的增大抑制了齿轮啮合过程中的振动峰值，使齿轮传动趋于平稳；谐波阻尼影响很小，在这个啮合状态分析齿轮副振动时将啮合阻尼系数视为常数是合理的。谐波刚度项系数的增大会使齿轮副啮合过程出现振动峰值，与此同时，谐波刚度项系数越大，亦会使某些频率点处峰值越大，不利于齿轮副的平稳传动。

(下转第 162页)162 振 动 与 冲 击 2013年第 32卷[10]ZHU Nian-yong.Yu Hong.HUANG Qing-an.Analysis ofnonlinear characteristic of nanomechanical resonators[J]。

Chinese Journal of Electron Devices，2005，28(1)：35-37。

夏晓嫒，李昕欣，王跃林.超薄硅纳米谐振梁的制作及谐振特性的测量[J].纳米技术与精密工程 ，2008，6(1)：1-4。

XIA Xiao-Yuan，LI Xin-xin，WANG Xue-Lin. Fabricationand detection of uhrathin silicon nano-beam resonator[J]。

Nanotechnology and Precision Engineering.2008，6(1)：l-4。

徐临燕 ，栗大超，胡小唐，等.基于原子力显微镜的纳米梁杨氏模 量 的测 量 [J].天津 大学 学 报，2007，40 (7)：816 -820。

(上接第 157页)(3) 1时，增大平均阻尼和降低主谐波刚度减弱了系统的非线性；谐波阻尼的增加对系统的振动特性影响较小，可采用常数阻尼代替非线性时变阻尼。

次谐波刚度系数的增大会产生不确定的情况，会使-些频率处的振动峰值降低，-些频率处振动峰值升高。