CN103234021B - 一种确定克林跟贝尔格锥齿轮时变啮合刚度的方法 - Google Patents
一种确定克林跟贝尔格锥齿轮时变啮合刚度的方法 Download PDFInfo
- Publication number
- CN103234021B CN103234021B CN201310144354.9A CN201310144354A CN103234021B CN 103234021 B CN103234021 B CN 103234021B CN 201310144354 A CN201310144354 A CN 201310144354A CN 103234021 B CN103234021 B CN 103234021B
- Authority
- CN
- China
- Prior art keywords
- crin
- bevel gear
- bel
- lattice
- gear
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
一种确定克林跟贝尔格锥齿轮时变啮合刚度的方法,属于机械设计与制造领域,步骤如下:1、将克林跟贝尔格锥齿轮系统简化处理成为齿轮副的扭转振动系统模型;2、在齿轮副的扭转振动系统模型中引入时变啮合刚度,建立克林跟贝尔格锥齿轮动力学方程;3、根据克林跟贝尔格锥齿轮时变啮合刚度原理,建立克林跟贝尔格锥齿轮多项式函数展开的时变啮合刚度方程式;4、根据克林跟贝尔格锥齿轮多项式函数展开的时变啮合刚度方程式求得其时变啮合刚度值,与基于传统方法建立的克林跟贝尔格锥齿轮时变啮合刚度值进行理论仿真分析和比较。本方法为得到该种齿轮时变啮合刚度值,为克林根贝尔格锥齿轮传动的减振、降噪、平稳传动动力学研究打下基础。
Description
技术领域
本发明属于机械设计与制造领域,涉及到一种确定克林跟贝尔格锥齿轮时变啮合刚度的方法,本方法更加精确的给出克林跟贝尔格锥齿轮时变啮合刚度值。
背景技术
克林根贝尔格锥齿轮具有重合度高、传动平稳、承载能力高等优点。广泛应用于各种机械设备中,如汽车、工程机械、旋翼推进的直升机、机床等。克林根贝尔格锥齿轮的制造精度、质量直接影响这些设备的效率、噪声、运动精度和寿命,因此克林根贝尔格锥齿轮一直受到各国有关专家学者的广泛关注和研究,成为齿轮生产中的关键技术和制高点。为了保证齿轮在传动过程中平稳可靠,减小振动和噪声,必须对齿轮传动系统进行动力学分析。有效的建立齿轮啮合刚度模型是动力学分析的基础。啮合刚度模型是否合理,直接影响了动力学分析的精确度。
国内外许多学者对齿轮误差以及动力学特性进行了较广泛而深入的研究,但多数都是利用有限元方法计算并获得了齿轮时变啮合刚度,他们采用的方法是把刚度简化处理成多阶谐波叠加的形式,这种形式只是保证了刚度的变化频率,但得到的时变啮合刚度与时变啮合刚度的实际特性相差很远。该种确定克林跟贝尔格锥齿轮时变啮合刚度的方法不仅可以简洁的得出齿轮时变啮合刚度值,而且更贴近于时变啮合刚度的实际特性。该方法对齿轮传动的减振、降噪、平稳传动等方面的动力学研究提供了理论依据。
发明内容
本发明的目的是提供一种确定克林跟贝尔格锥齿轮时变啮合刚度的方法,针对航空航天机械传动中用的克林跟贝尔格锥齿轮,提出采用多项式的形式对克林跟贝尔格锥齿轮时变啮合刚度进行计算。并在不同的阻尼、载荷、转速工况下,将多项式函数展开的刚度值和多阶谐波叠加形式的刚度值进行分析对比,得到本方计算方法的时变啮合刚度能更贴近于时变啮合刚度的实际特性。
本发明是采用以下技术手段实现的:
一种确定克林跟贝尔格锥齿轮时变啮合刚度的方法,其包括如下步骤:
S1、将克林跟贝尔格锥齿轮系统简化处理成为齿轮副的扭转振动系统模型;
S2、在齿轮副的扭转振动系统模型中引入时变啮合刚度,建立克林跟贝尔格锥齿轮动力学方程。动力学方程如下:
式中:Ii(i=p,g)为主、被动齿轮的转动惯量;λi(i=p,g)为主、被动齿轮的齿轮方向旋转半径;θi(i=p,g)为主、被动齿轮的角位移;Ti(i=p,g)为主、被动齿轮上的扭矩;e为静态误差,c为啮合阻尼,k为啮合刚度。
动态传递误差δ定义为δ=λpθp-λgθg带入2中动力学方程,并化简得:
上式中,p=δ-e,
其中,引入新的变量Tp=Tg=T,λp=λg=λ,将其代入上面的平衡方程中得:
其中:ξ为阻尼系数,f′zzc为齿轮等级精度影响齿轮副齿频周期误差的公差值;p为接触点沿齿面法向位移;km即为克林跟贝尔格锥齿轮的啮合刚度;cm为克林跟贝尔格锥齿轮的啮合阻尼;
S3、建立克林跟贝尔格锥齿轮多项式函数展开的时变啮合刚度方程式。
根据等高齿锥齿轮啮合刚度的变化规律,齿轮啮合刚度的大小与啮合区域的形状及其位置有直接的关系。一对轮齿接触时,在啮合力的作用下,沿齿面法向的弹性位移p越大,啮合面积也就越大。因此,弹性位移p的值一定程度上对齿轮综合啮合刚度起决定作用。基于以上分析,假设啮合刚度k'm关于沿齿面法向位移p的函数关系具有多项式的形式:
其中,k'i为刚度影响系数,i=1,2,3,...,n。
S4、根据克林跟贝尔格锥齿轮多项式函数展开的时变啮合刚度方程式求得其时变啮合刚度值,与基于传统方法建立的克林跟贝尔格锥齿轮时变啮合刚度值进行理论仿真分析和比较。
本发明的特点在于基于多项式理论,提出了一种确定齿轮时变啮合刚度的方法,通过此方法可以准确地得到克林根贝尔格锥齿轮时变啮合刚度值,使其更贴近于时变啮合刚度的实际特性。发明内容包括四部分。在第一部分中,将克林跟贝尔格锥齿轮系统简化处理成为齿轮副的扭转振动系统模型;在第二部分中,在齿轮副的扭转振动系统模型中引入时变啮合刚度,建立克林跟贝尔格锥齿轮动力学方程,并将方程化简;在第三部分中,根据克林跟贝尔格锥齿轮时变啮合刚度原理,建立克林跟贝尔格锥齿轮多项式函数展开的时变啮合刚度方程式;在第四部分中,根据克林跟贝尔格锥齿轮多项式函数展开的时变啮合刚度方程式求得其时变啮合刚度值,与基于传统方法建立的克林跟贝尔格锥齿轮时变啮合刚度值进行理论仿真分析和比较。
通过下面的描述并结合附图说明,本发明会更加清晰,附图说明用于解释本发明方法及实施例。
附图说明
图1实际齿面上接触区域;
图2两个自由度的等高齿锥齿轮模型;
图3大小齿轮坐标系;
图4不同转速、负载扭矩下最大振动位移百分比×100%(i=2或3);
图5在n=500r/min、T=2200Nm下不同阻尼耗散因子系数下趋于稳定时间;
图6在n=500r/min、T=2200Nm下不同阻尼耗散因子系数下稳定时间差值百分比×100%(i=2或3);
图7为本发明的步骤流程图。
具体实施方式
本发明的具体实施步骤如下:
第一步:将克林跟贝尔格锥齿轮系统简化处理成为齿轮副的扭转振动系统模型;
本实施例以航空用某摆线锥齿轮副为研究对象,其具体参数见表1。克林跟贝尔格锥齿轮副动力学模型如图2所示。在该模型中,假设两齿轮的支撑刚度较大,且不考虑传动轴、支承轴承和箱体等的弹性变形对摆线锥齿轮系统的影响,最终将摆线锥齿轮系统简化处理成为齿轮副的扭转振动系统模型。
表1摆线锥齿轮系统参数
第二步:在齿轮副的扭转振动系统模型中引入时变啮合刚度,建立克林跟贝尔格锥齿轮动力学方程。
图2为克林跟贝尔格锥齿轮扭转振动模型,双自由度动态模型(主动轮p,从动轮g)转动方程如下:
式中:e为静态误差,c为啮合阻尼,k为啮合刚度,λi为啮合点半径,Ii为齿轮的转动惯量,θi为齿轮转动角度,Ti为齿轮的转动力矩。参数如图2所示。动态传递误差定义为并带入上式中:
δ=λpθp-λgθg
式中:p=δ-e,f′zzc为齿轮等级精度影响齿轮副齿频周期误差的公差值。为了简化计算,令Tp=Tg=T,λp=λg=λ上述方程简化为:
其中: (ξ为阻尼系数)。
第三步:建立克林跟贝尔格锥齿轮多项式函数展开的时变啮合刚度方程式;
在某一啮合点处,已知啮合点的法向位移p,根据齿轮的几何形状和该点的几何位置可以确定椭圆的具体尺寸。轮齿啮合刚度随着接触区域的变化而变化,接触区域越大,两个齿面啮合在一起的部分越多,越难以再继续产生变形,即刚度值越大。当齿轮传递的载荷F较小时,接触点在齿面法向方向上产生很小的位移Δp,此时接触椭圆区域也很小,此时啮合刚度k较小;当传递的载荷很大时,形成的接触椭圆区域变大了,啮合刚度k也随之变大。综上分析,可以认为刚度k值由椭圆接触区域的大小和形状决定,而椭圆接触区域又由接触点沿齿面法向的位移p决定。于是,假设单齿啮合刚度k是关于接触点沿齿面法向位移p的函数。
根据等高齿锥齿轮啮合刚度的变化规律,齿轮啮合刚度的大小与啮合区域的形状及其位置有直接的关系。在任一时刻,啮合区域中心点的位置和啮合面积的大小直接决定了啮合刚度。一对轮齿接触时,在啮合力的作用下,沿齿面法向的弹性位移p越大,啮合面积也就越大。因此,弹性位移p的值一定程度上对齿轮综合啮合刚度起决定作用。基于上文分析,假设啮合刚度k'm关于沿齿面法向位移p的函数关系具有多项式的形式
其中,k'i为刚度影响系数,i=1,2,3,...,n。
第四步:根据克林跟贝尔格锥齿轮多项式函数展开的时变啮合刚度方程式求得其时变啮合刚度值,与基于传统方法建立的克林跟贝尔格锥齿轮时变啮合刚度值并进行理论仿真分析和比较;
4.1平均啮合刚度值根据机械振动手册中查得:
k1=9.12×108
4.2利用传统方法计算克林跟贝尔格锥齿轮时变啮合刚度值,典型的刚度模型建立为:
在求解时由于计算量的限制,往往只取第一阶简谐波,建立起来的阻尼cm和静态传递误差e的模型具有类似的形式为:
k2=9.12×108×[1+0.2cos(ωt+φ)]
式中齿轮的啮合频率ω=nπz/30,k0=1,kaj=0.2,k0=9.12×108。
4.3利用多项式函数展开方法计算克林跟贝尔格锥齿轮时变啮合刚度值;
由方程:
中只有前几项系数K'0,K'1,K'2起较大作用,因此计算分析时只研究前几项即可。取K'0=9.12×108,得时变啮合刚度值为:
k3=9.12×108×(Ap2+Bp)
式中A与B为多项式的系数,其数值为
A与B的数值是根据机械振动手册中,齿轮传动中齿轮副的平均啮合刚度km与重合度ε=1.6,轴向重合度ε1=0.6的函数关系图查得的。
4.4K1,K2,K3的时域响应趋于稳定时间比较
由于转速和负载扭矩的大小对齿轮振动趋于稳定时所用的时间影响不大,为了便于分析比较所以取转速n=500r/min、负载扭矩T=2200Nm,比较k3和k2在不同阻尼耗散因子系数下,齿轮振动趋于稳定的时间。表2给出了三种刚度模型在不同阻尼系数下,振动趋于平衡的时间以及趋于平衡时间差值的百分比,图7直观给出三个刚度模型振动平衡时间,并且可以看出k3的平衡时间更贴近于k1,振动趋于平衡的时间差值百分比可以直观地看出刚度模型k3的振动趋于平衡的时间差值百分比要比刚度模型k2要小,尤其是在小阻尼情况下这种情况更加明显。
表2在n=500r/min、T=2200Nm下不同阻尼耗散因子系数下趋于稳定时间的比较
由于克林根贝尔格锥齿轮具有重合度高,且重合度大多不是整数,给建立该种齿轮时变啮合刚度计算模型带来一定的难度。且啮合过程中同时参与啮合的齿对数随时间成周期性变化,因而轮齿的啮合综合刚度是随时间周期变化的。这同样给建立该种齿轮时变啮合刚度计算模型带来一定的难度。而本发明方法能够较为准确的建立该种齿轮时变啮合刚度计算模型,并且使得在小阻尼、重载工况下这种多项式形式的时变啮合刚度模型能更贴近于时变啮合刚度的实际特性。该模型为以后对齿轮传动的减振、降噪、平稳传动等方面的动力学研究打下基础。
通过以上实例分析总结出:本发明方法能够运用于克林根贝尔格锥齿轮的时变啮合刚度的计算,并能快速准确地得出时变啮合刚度值。本发明方法不仅为克林根贝尔格锥齿轮传动的减振、降噪、平稳传动等方面的动力学研究打下基础,且为齿轮的时变啮合刚度计算提供了有效借鉴。
Claims (1)
1.一种确定克林跟贝尔格锥齿轮时变啮合刚度的方法,其特征在于:其具体步骤如下:
S1、将克林跟贝尔格锥齿轮系统简化处理成为齿轮副的扭转振动系统模型;
S2、在齿轮副的扭转振动系统模型中引入时变啮合刚度,建立克林跟贝尔格锥齿轮动力学方程;动力学方程如下:
式中:Ii(i=p,g)为主、被动齿轮的转动惯量;λi(i=p,g)为主、被动齿轮的齿轮方向旋转半径;θi(i=p,g)为主、被动齿轮的角位移;Ti(i=p,g)为主、被动齿轮上的扭矩;e为静态误差,c为啮合阻尼,k为啮合刚度;
动态传递误差δ定义为δ=λpθp-λgθg带入S2中动力学方程,并化简得:
上式中,p=δ-e,
其中,引入新的变量Tp=Tg=T,λp=λg=λ,将其代入上面的平衡方程中得:
其中:ξ为阻尼系数,f′zzc为齿轮等级精度影响齿轮副齿频周期误差的公差值;p为接触点沿齿面法向位移;km即为克林跟贝尔格锥齿轮的啮合刚度;cm为克林跟贝尔格锥齿轮的啮合阻尼;
S3、建立克林跟贝尔格锥齿轮多项式函数展开的时变啮合刚度方程式;
根据等高齿锥齿轮啮合刚度的变化规律,齿轮啮合刚度的大小与啮合区域的形状及其位置有直接的关系;一对轮齿接触时,在啮合力的作用下,沿齿面法向的弹性位移p越大,啮合面积也就越大;因此,弹性位移p的值一定程度上对齿轮综合啮合刚度起决定作用;基于以上分析,假设啮合刚度k'm关于沿齿面法向位移p的函数关系具有多项式的形式:
其中,k'i为刚度影响系数,i=1,2,3,...,n;
S4、根据克林跟贝尔格锥齿轮多项式函数展开的时变啮合刚度方程式求得其时变啮合刚度值,与基于传统方法建立的克林跟贝尔格锥齿轮时变啮合刚度值进行理论仿真分析和比较。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310144354.9A CN103234021B (zh) | 2013-04-23 | 2013-04-23 | 一种确定克林跟贝尔格锥齿轮时变啮合刚度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310144354.9A CN103234021B (zh) | 2013-04-23 | 2013-04-23 | 一种确定克林跟贝尔格锥齿轮时变啮合刚度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103234021A CN103234021A (zh) | 2013-08-07 |
CN103234021B true CN103234021B (zh) | 2015-09-02 |
Family
ID=48882074
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310144354.9A Active CN103234021B (zh) | 2013-04-23 | 2013-04-23 | 一种确定克林跟贝尔格锥齿轮时变啮合刚度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103234021B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105784360A (zh) * | 2016-05-12 | 2016-07-20 | 重庆长安汽车股份有限公司 | 一种基于啮合接触线长度变化确定齿轮啮合动刚度的方法 |
CN107131282B (zh) * | 2017-06-09 | 2019-05-31 | 江苏大学 | 融合齿背接触机理的高速圆柱齿轮动态啮合刚度计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4651588A (en) * | 1986-03-03 | 1987-03-24 | Rouverol William S | Low-excitation gearing |
CN101246083A (zh) * | 2008-03-24 | 2008-08-20 | 西安电子科技大学 | 直齿圆柱齿轮动态啮合刚度的测量方法 |
CN101770538A (zh) * | 2010-01-15 | 2010-07-07 | 北京工业大学 | 含损伤性单齿故障圆柱直齿轮啮合刚度仿真分析方法 |
-
2013
- 2013-04-23 CN CN201310144354.9A patent/CN103234021B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4651588A (en) * | 1986-03-03 | 1987-03-24 | Rouverol William S | Low-excitation gearing |
CN101246083A (zh) * | 2008-03-24 | 2008-08-20 | 西安电子科技大学 | 直齿圆柱齿轮动态啮合刚度的测量方法 |
CN101770538A (zh) * | 2010-01-15 | 2010-07-07 | 北京工业大学 | 含损伤性单齿故障圆柱直齿轮啮合刚度仿真分析方法 |
Non-Patent Citations (3)
Title |
---|
基于有限元法的螺旋锥齿轮啮合刚度计算;唐进元等;《机械工程学报》;20110630;第47卷(第11期);第23-28页 * |
直齿锥齿轮轮齿变形及瞬时啮合刚度;李剑锋等;《山东工业大学学报》;19961231;第26卷(第4期);第451-454页 * |
螺旋锥齿轮啮合刚度及参数振动稳定性研究;王延忠等;《航空动力学报》;20100731;第25卷(第7期);第1664-1668页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103234021A (zh) | 2013-08-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107092750B (zh) | 一种内啮合齿轮轴减速器轮齿损伤故障的非线性动力学建模方法 | |
CN103234021B (zh) | 一种确定克林跟贝尔格锥齿轮时变啮合刚度的方法 | |
Zhou et al. | Vibration and noise radiation characteristics of gear transmission system | |
CN109190227A (zh) | 基于解析-有限元的复杂齿基直齿圆柱齿轮啮合刚度计算方法 | |
Shi et al. | Dynamic analysis of three-dimensional drivetrain system of wind turbine | |
CN104091007A (zh) | 一种少齿差行星减速器动力学仿真分析方法 | |
CN103473479B (zh) | 一种履带车辆传动系统扭转振动数学模型建立方法 | |
CN103345583A (zh) | 一种齿轮传动系统非线性动力学建模方法 | |
Qin et al. | Design and evaluation of two-stage planetary gearbox for special-purpose industrial machinery | |
Guo et al. | A new algorithm to solve meshing-in impact considering the measured pitch error and to investigate its influence on the dynamic characteristics of a gear system | |
CN115203898A (zh) | 一种受内部多源故障激励的齿轮箱耦合动力学建模方法 | |
Ren et al. | Investigation of dynamic load sharing behavior for herringbone planetary gears considering multicoupling manufacturing errors | |
CN103207939B (zh) | 一种基于变传动比摆线齿锥齿轮动力学建模方法 | |
Zhang et al. | Study on Empirical Gear Profile Micro-Modifications for Gear Transmission | |
Nigam et al. | Modelling and structural analysis of planetary geared winch | |
Shan et al. | Analysis of nonlinear dynamic accuracy on RV transmission system | |
CN105784360A (zh) | 一种基于啮合接触线长度变化确定齿轮啮合动刚度的方法 | |
Wang et al. | Bifurcation and chaos analysis of gear system with clearance under different load conditions | |
CN105069209B (zh) | 一种直升机主减速器行星架裂纹故障动力学模型建模方法 | |
Lee et al. | Vibrational characteristics of automotive transmission | |
Yang et al. | Research on Dynamic Characteristics of Confluence Planetary Row under Speed Impact | |
Yu et al. | NVH characteristic prediction of an electric vehicle reducer | |
Rosbi et al. | Evaluation of dynamic absorber to suppress subharmonic nonlinear vibration in car powertrain | |
Rashid et al. | Dynamic behaviour analysis of single stage planetary gearbox under non-stationary operations using experimental test rig | |
CN103198226B (zh) | 一种考虑摩擦的摆线锥齿轮振动特性分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |