谐波齿轮的侧隙规律研究与有限元模型仿真

第２９卷第６期
中国机械工程
V o l ．２９㊀N o ．６２０１８年３月
C H I N A M E C HA N I C A LE N G I N E E R I N G
p p
．６５９Ｇ６６６谐波齿轮的侧隙规律研究与有限元模型仿真
杨朋朋１,２㊀陈晓霞１,２㊀邢静忠１,２㊀姚云鹏１,２
１．天津工业大学机械工程学院,天津,３００３８７
２．天津市现代机电装备技术重点实验室,天津,３００３８７
摘要:为更真实地反映谐波齿轮的侧隙分布,改进侧隙算法中柔轮齿根的定位方式,建立坐标变换下的齿廓方程以代替原有侧隙算法中的齿厚方程,提出基于周向位移定位和弧长定位的侧隙计算方法,建立含渐开线齿廓的平面齿圈实体有限元模型,获得空载啮合状态下侧隙分布.将理论计算的侧隙值与有限元模型计算的侧隙值比较发现,两者所得结果一致性较高.同时为了揭示侧隙偏差的来源,获取了有限元模型柔轮中性层的径向位移㊁周向位移和法线转角,并求解了周向位置极角.与理论算法结果比较发现,柔轮齿根周向位移偏差是引起侧隙偏差的主要因素
.
关键词:谐波齿轮;齿根定位;侧隙算法;有限元模型中图分类号:T H １３２
D O I :１０．３９６９/j ．i s s n ．１００４ １３２X．２０１８．０６．００５开放科学(资源服务)标识码(O S I D )
:R e s e a r c ho nB a c k l a s hR e g u l a r i t y a
n dF i n i t eE l e m e n t S i m u l a t i o no f H a r m o n i cG e a r s
Y A N GP e n g p e n g １,
２㊀C H E N X i a o x i a １,
２㊀X I N GJ i n g z h o n g １,
２㊀Y A O Y u n p e n g
１,
２
１．S c h o o l o fM e c h a n i c a l E n g i n e e r i n g ,T i a n j i nP o l y t e c h n i cU n i v e r s i t y ,T i a n j
i n ,３００３８７２．T i a n j i nK e y L a b o r a t o r y o fM o d e r n M e c h a t r o n i c sE q u i p m e n tT e c h n o l o g y ,T i a n j
i n ,３００３８７A b s t r a c t :I no r d e r t or e f l e c t t h eb a c k l a s hd i s t r i b u t i o n so fh a r m o n i c g e a r s m o r ea c c u r a t e l y
,t h e t o o t h r o o t p o s i t i o n i n g m e t h o d s o f f l e x s p l i n ew e r e i m p r o v e d i n t h eb a c k l a s ha l g o r i t h m ,a n d t h e t o o t h p r o f i l e e q u a t i o n su n d e r t h ec o o r d i n a t e t r a n s f o r m a t i o n w a se s t a b l i s h e d i n s t e a do f t h e t o o t ht h i c k n e s s e q u a t i o n s i n t h e o r i g i n a l b a c k l a s ha l g o r i t h m．B o t ho f t h em e t h o d s f o r c a l c u l a t i n g t h eb a c k l a s h sb a s e d o n c i r c u m f e r e n t i a l d i s p l a c e m e n t p o s i t i o n i n g a n da r c l e n g t h p o s i t i o n i n g w e r e p r o p
o s e d ．T h e f i n i t ee l e Ｇm e n tm o d e l o f p l a n a r t o o t h r i n g
sw i t h i n v o l u t e t o o t h p r o f i l ew a s e s t a b l i s h e d ,a n d t h e b a c k l a s hd i s t r i Ｇb u t i o n s u n d e r n o Ｇl o a d c o n d i t i o n sw a s o b t a i n e d ．C o m p a r i n g
t h e t h e o r e t i c a l b a c k l a s h sw i t h t h a t f r o mt h e f i n i t e e l e m e n tm o d e l ,t h e i r c o n s i s t e n c y
i s e x c e l l e n t ．I n o r d e r t o r e v e a l t h e c a u s e s o f t h e b a c k l a s h d e v i a Ｇt i o n s ,t h e r a d i a l d i s p l a c e m e n t ,t h e c i r c u m f e r e n t i a l d i s p l a c e m e n t s a n d t h e n o r m a l r o t a t i o n a n g
l e s o f t h e f l e x s p l i n en e u t r a l l a y
e r s i n t h e
f i n i t e e l e m e n tm o d e lw e r eo b t a i n e d ,a n d t h e c i r c u m f e r e n t i a l p o l a r a n Ｇ
g l e sw e r e s o l v e d t o c o m p a r ew i t
h t h e t h e o r e t
i c a l a l g o r i t h m．T h e c i r c u m f e r e n t i a l d i s p l a c e m e n t s o f t h e f l e x s p l i n e t e e t ha r e t h em a j
o r f a c t o r o f t h eb a c k l a s hd e v i a t i o n s ．K e y w
o r d s :h a r m o n i c g e a r ;t o o t h r o o t p o s i t i o n i n g ;b a c k l a s ha l g o r i t h m ;f i n i t e e l e m e n tm o d e l 收稿日期:２０１７０３２０
基金项目:国家自然科学基金资助项目(５１５７５３９０);天津市应用基础与前沿技术研究一般项目(１４J C Y B J C １９２００
)０㊀引言
谐波齿轮传动技术是在２０世纪５０年代末发
展起来的一种新型机械传动技术.谐波齿轮具有运动精度高㊁质量和体积小㊁传动比大㊁参与啮合齿数多和传动效率高等独特优点,广泛应用于机
器人㊁航空航天㊁光学仪器等领域[
１
].在设计时保证齿间合理的侧隙是保持谐波齿轮这些优越性能的前提,侧隙不足,就会在过载时导致干涉.因此
如何获得齿间合理的侧隙,保证轮齿不发生干涉,
是一个很重要的课题[
２Ｇ４
].伊万诺夫[２]
对齿廓干涉条件以及侧隙理论计
算方法进行了研究.沈允文等[３
]对谐波齿轮传动
的啮合干涉问题以及侧隙计算问题建立了理论方
法.辛洪兵等[５]建立了初始啮合侧隙计算的数学
模型,研究了谐波齿轮初始啮合侧隙的变化规律.
殷燕[６]进行了零侧隙渐开线谐波齿轮传动的参数
优化设计及有限元分析,以零侧隙为目标函数,完
成了啮合参数的优化.C H E N 等[７
Ｇ８]提出了双圆弧公切线齿廓谐波齿轮的侧隙计算方法,并依据
９５６ Copyright ©博看网. All Rights Reserved.
齿间侧隙进行了干涉检查.刘邓辉等[９]考虑柔轮筒体锥度变形进行了空间齿廓设计,并求解了空
间齿廓的侧隙分布.以上的研究虽有涉及侧隙计
算,但并未考虑柔轮齿根定位方式与柔轮中性层
的位移情况对周向侧隙(下文涉及周向侧隙以柔
轮齿顶点为计算点)的影响.
关于谐波齿轮仿真模型,国内外已有多人进
行过相关研究.董惠敏等[１０]建立了谐波齿轮传动中柔轮在空载和负载时板壳的有限元分析模
型.刘文芝等[１１]以杯形柔轮为例,建立了柔轮啮合的仿真实体模型,用三维弹性接触有限元法计
算和分析了承载柔轮齿圈和筒体的应力大小及分
布规律.O S T A P S K I等[１２]提出了一种复杂形状薄壁壳结构的弹性变形问题的几何非线性壳理论
求解方法,并通过有限元方法进行了计算.付军
锋等[１３]建立了柔轮的三维实体有限元分析模型,并对柔轮模型和波发生器模型在接触条件下进行了有限元分析.以上涉及的柔轮齿圈模型多为当量厚度的壳单元模型.
本文在文献[２Ｇ３]谐波齿轮柔轮齿根定位方式与侧隙计算方法基础上,建立了四滚轮波发生器作用下能够准确表达齿廓信息的平面齿圈有限元模型,并将上述理论计算侧隙结果与有限元模型计算结果进行比对,分析两者之间的偏差.通过改进柔轮齿根定位方式和建立坐标变换下的齿廓方程,提出基于周向位移定位的侧隙计算方法与基于弧长定位的侧隙计算方法.同时为揭示侧隙偏差的来源,获取了有限元模型中性层的径向位移㊁周向位移和法线转角,并求解了周向位置极角来与理论公式计算结果进行比较.
１㊀问题的引出
１．１㊀现有文献侧隙算法的分析
柔轮齿廓与刚轮齿廓的侧隙分布情况是齿轮
啮合性能评价的重要指标.在理论方法中,柔轮
齿根的定位方式和柔轮中性层的径向位移㊁周向
位移,以及齿对称轴线相对于矢径的转角(以下简
称法线转角)是影响侧隙结果的重要因素.目前
主要有两种侧隙计算方法.
(１)文献[２]计算的侧隙为在波发生器旋转状态下,初始状态位于柔轮变形短轴处的一对柔轮齿廓与刚轮齿廓在不同时刻下的周向侧隙.周向侧隙j a g定义为
j a g＝|v a g|－(s a g＋s y b)/２(１)式中,s a g㊁s y b分别为柔轮齿顶圆处齿厚与对应刚轮处齿厚;v a g为柔轮齿对称线和齿顶圆交点与此点到刚轮齿对
称线交点的距离.
(２)文献[３]计算的侧隙为与柔轮变形长轴不同夹角的任意啮合位置柔轮齿顶与刚轮齿廓间的周向侧隙.该算法以变形前的柔轮弧长等于变形后的柔轮弧长为齿根定位方式,建立了柔轮与刚轮的齿廓参数方程,最后将柔轮齿顶到刚轮齿廓间的最短距离作为齿廓间的周向侧隙.
为验证理论方法的有效性,本文建立准确表达渐开线齿廓的平面齿圈有限元模型,对以上两种理论侧隙算法的结果进行仿真验证.
１．２㊀有限元实体模型的建立
设置齿轮模数m＝０．２,柔轮齿数z g＝１４０,刚轮齿数z b＝１４２,柔轮变位系数x g＝２．１３,刚轮变位系数x b＝１．９２５,基准齿形角α＝２０ʎ,径向变形系数w∗０＝１,波数U＝２,柔轮壁厚δ＝０．３mm,齿顶高系数h∗a＝１,径向间隙系数c∗＝０．３５.
考虑到结构的对称性,基于以上参数建立含齿圈的１/４有限元模型,柔轮㊁刚轮和波发生器均选用P l a n e１８３单元.在建立柔轮的齿廓时,利用二分法求解从柔轮齿顶到齿根６个均布的齿廓参数u,基于齿廓方程得到６个不同的关键点坐标,将其连线表达为柔轮齿廓.刚轮齿廓建模过程与此相同.然后定义波发生器外表面与柔轮内壁的接触关系,接触对间选用刚体柔体的面面接触单元,以波发生器的上半圆为目标面,定义T a r g e１６９目标单元;以柔轮的内壁为柔性接触面,定义C o n t a１７２接触单元.在柔轮变形长轴区和短轴区施加对称位移边界条件,选用大变形求解选项.图１为经过后处理得到的柔轮齿圈径向位移U X .
图１㊀有限元模型径向位移分布
F i g．１㊀R a d i a l d i s p l a c e m e n t d i s t r i b u t i o no f f i n i t e
e l e m e n tm o d e l
在有限元模型中求解柔轮齿顶侧隙,首先需获取柔轮齿顶的节点坐标,继而寻找与该点极径相等的刚轮齿廓上的节点;判断柔轮和刚轮齿顶点是否位于啮合区间内,即设柔轮齿顶节点坐标A(x a１,y a１),刚轮齿顶节点坐标B(x a２,y a２),若
０６６
中国机械工程第２９卷第６期２０１８年３月下半月Copyright©博看网. All Rights Reserved.
判断满足
(x ２a １＋y ２
a １)ɤ
(x ２a ２＋y ２
a ２)
则齿轮顶点处于啮合区间内(
图２),柔轮和刚轮齿顶点间的距离即为齿廓的侧隙.
图２㊀柔轮与刚轮啮合状态图
F i g ．２㊀E n g a g e m e n t d i a g r a mo f f l e x s p
l i n e a n d c i r c u l a r s p
l i n e １．３㊀有限元模型结果验证
图３的纵坐标表示柔轮与刚轮齿廓间的周向侧隙j t ,
横坐标表示柔轮齿圈部各啮合齿的位置角度φ.在φɪ(－１０ʎ,３０ʎ)区间内,文献[２]和文献[３]两种算法结果与有限元模型计算结果差距较小;而在φɪ(３０ʎ,６０ʎ
)区间内,三者计算结果差距都较大,其中文献[２]
结果与有限元模型计算结果偏差最大.产生偏差的原因分析如下.
图３㊀理论计算方法与有限元模型计算方法得到的
周向侧隙对比
F i g ．３㊀C o m p a r i s o nd i a g r a mo f t h e c i r c u m f e r e n t i a l b a c k l a s hb e t w e e n t h e t h e o r e t i c a l a l g
o r i t h ma n d t h e f i n i t e e l e m e n tm o d e l
(１)文献[２
]算法计算侧隙结果与有限元模型计算结果偏差相对较大,主要是定位柔轮齿根和侧隙计算位置出现偏差.求解柔轮齿顶点的周向坐标v a g :
v a g ＝v ＋(r a g －r )θ－(r a g ＋w )φb －v b
(２
)式中,v ㊁w 分别为柔轮齿根运动的周向位移和径向位移;r a g ㊁
r 分别为柔轮齿顶圆半径与分度圆半径;θ为轮齿转动引起的法线转角;φb 为柔轮与刚轮相对转过的角度;
v b 为变形后刚轮的周向位移(忽略不计)
.柔轮齿顶圆变形后径向坐标:
w a g ＝(r a g ＋w )c o s φb －r －w b
(３
)式中,w b 为变形后刚轮的径向位移(
忽略不计).式(３)求解径向坐标w a g 时,
未考虑法线转角引起的轮齿转动;同时式(１)中计算的并不是真正的柔轮齿顶点,而是柔轮齿对称线与齿顶圆的交点,这些简化也会造成侧隙偏差.
(２
)如图４所示,在柔轮齿和刚轮齿的横剖面内,设柔轮齿坐标系S f {o f ,x f ,y f }原点位于柔轮中性层上,y f 与柔轮齿对称线重合,柔轮变形长轴与y b 轴重合.设中性层半径为r m ,原始曲线矢径为ρ,轴y b 与柔轮输出端矢径o b o １的夹角
为φ,与柔轮变形端矢径o b o f 的夹角为φ１[
３
],当波发生器从y b 轴顺时针旋转时,φ值为正.
图４㊀柔轮与刚轮齿廓侧隙图
F i g ．４㊀B a c k l a s hd i a g r a mo f t h e f l e x s p
l i n e t o o t h p r o f i l e a n d t h e c i r c u l a r s p
l i n e t o o t h p r o f i l e 文献[３
]算法与有限元模型计算结果出现偏差的原因是在进行柔轮齿根的定位时,考虑到求
解积分公式的复杂性,将弧长公式[１４
]近似取
φ１＝φ＋v
/r m (４
)２㊀齿根定位方式与侧隙算法的改进
２．１㊀基于周向位移定位的侧隙算法
针对文献[２]中柔轮齿根的定位方式,将周向位移简化为切向位移,即用o ３o f 代替o ２o f ,利用变形端柔轮齿根转过角度来定位(见图５中Δo w o ３o f ).这种方式定位虽也有近似,但是θv 足
够小,此种定位方式的侧隙偏差要小于文献[２]计算的侧隙偏差.同时利用坐标变换方法建立柔轮与刚轮的齿廓方程,求解柔轮齿顶的侧隙.利用此种定位方式的侧隙计算方法称为基于周向位移定位的侧隙算法(简称位移法)
.图５为轮齿对称线初始位置与柔轮变形长轴
相差π/２个相位的柔轮齿F 与刚轮齿B 啮合示意图.设刚轮坐标系S b {o b ,x b ,y b }
固定,柔轮变形端坐标系为S f {o f ,x f ,y f }
,波发生器坐标系为
１６６ 谐波齿轮的侧隙规律研究与有限元模型仿真
杨朋朋㊀陈晓霞㊀邢静忠等Copyright ©博看网. All Rights Reserved.
S w {o w ,x w ,y w }
.定义输出端矢径o b o ２与柔轮变形长轴y w 的夹角φ为自变量,与刚轮坐标系x b 的夹角为φg ;定义柔轮变形端矢径o b o f 与柔轮变形长轴y w 的夹角为φ１,柔轮变形端矢径o b o f 与输出端矢径o b o ２的夹角为θv (
定义为周向位置极角).柔轮变形长轴y w 与y b 轴的夹角为φw ,柔
轮齿廓对称线y f 与刚轮坐标轴x b 的夹角为ψ;刚
轮齿顶圆和齿根圆半径分别是r a b 和r f b .当波发生器从y b 轴顺时针旋转时,φ值为正
.
图５㊀谐波齿轮传动的几何关系
F i g ．５㊀
G e o m e t r i c a l r e l a t i o n s h i p o
f h a r m o n i c d r i v e 谐波齿轮运动转换关系如下:波发生器从初
始位置转动到任意角φw ,柔轮输出端相对刚轮转过角度φg ,柔轮齿根产生径向位移w 和周向位移
v 时,柔轮变形端齿根由o １运动到o f ,
柔轮齿对称线相对于齿根矢径转过角度μ.当φ＝０ʎ时,柔轮变形长轴y w 与轴y f 重合,柔轮齿与刚轮齿处于完全啮入状态.当φ＝９０ʎ时,轴x w ㊁柔轮齿对称线y f 与刚轮齿对称线x b 重合,柔轮齿与刚轮齿处于齿顶对齿顶的完全脱开状态,且有
φg ＝
φw i b h g ＝π－２φ
２(i b
h g ＋１
)(５)φ１＝φg －θ
v (６)ψ＝θ
v －φg ＋μ(７)θv ＝a
r c s i n [v /(w ＋r m )](８)ρ＝r m ＋w
(９
)式中,i b
h g 为刚轮b 固定时的传动比.
(１
)确定柔轮齿廓方程.先给出在坐标系S w {o w ,x w ,y w }中位于柔轮变形长轴的渐开线柔
轮右齿廓的参数方程[
３
]:x １＝r １[－s i n (u a １－θ１)＋u a １c o s α０c o s (u a １－θ１＋α０)
]y １＝r
１[c o s (u a １－θ１)＋u a １c o s α０s i n (u a １－θ１＋α０)]}
(１０
)式中,u a １为渐开线柔轮齿顶参数值;θ１为柔轮齿分度圆
齿厚所对中心角之半;α０为基准齿形角;r １为柔轮分度圆半径.
通过坐标变换得到在坐标系S b {o b ,x b ,y b }
下的柔轮齿廓方程:
S １＝x a １y a １１éëêêêùûúúú＝R １
R ２x １y １１é
ëêêêùû
úúú(１１
)R １＝c o s (π/２)－s i n (π/２)－r m s i n (π/２)c o s (π/２)０００１éëêêêùû
úúúR ２＝c o s ψ－s i n ψρc o s φ１s i n ψc o s ψρs i n φ１００１éëêêêùû
ú
ú
ú基于上述柔轮齿廓方程得到柔轮齿顶坐标为
M １(x a １,y a
１).(２
)确定刚轮齿廓方程.给出在坐标系S w {o w ,x w ,y w }中位于柔轮变形长轴的渐开线
刚轮齿槽右齿廓的参数方程[
３]
:x ２＝r ２[
－s i n (u M ２－θ２)＋u M ２c o s α０c o s (u M ２－θ２＋α０)]y ２＝r
２[c o s (u M ２－θ２)＋u M ２c o s α０s i n (u M ２－θ２＋α０)
]üþ
ýï
ï
ïï(
１２)式中,u M ２为渐开线刚轮齿廓上对应点处的参数值;θ２为
刚轮齿分度圆齿厚所对中心角的１/２;r ２为刚轮分度圆半径.
通过坐标变换得到坐标系S b {o b ,x b ,y b }下
的刚轮齿廓方程:
S ２＝x M ２y M ２１
éëêêêùû
úúú＝R ３x ２y ２１éëêêêùûú
ú
ú
(１３
)R ３＝c o s φ２－s i n φ２０s i n φ２c o s φ２０００１éëêêêùû
ú
ú
úφ２＝π
/２－π/z b (１４)式中,φ２为齿槽对称线位于y b 轴的刚轮左齿廓旋转到齿对称线位于x b 轴所运动过的角度.
基于上述刚轮齿廓方程可得到与柔轮齿顶点
极径相等的点M ２(x M ２,y M ２)
.(３
)四滚轮波发生器作用下柔轮中性层的径向位移[２]
为
w ＝w ０(C c o s φ＋φs i n βs i n φ－π４)
C －
π４㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀０ɤφɤβ
w ０[D s i n φ＋(π２－φ)c o s βc o s φ－π４]C －π
４
β＜φɤπ
/２ìîíïï
ïïïïïï
ïïïï(１５
)
２６６ 中国机械工程第２９卷第６期２０１８年３月下半月Copyright ©博看网. All Rights Reserved.
C ＝s i n β＋(π/２－β)c o s β
D ＝c o s β＋βs i n β
式中,β为四滚轮波发生器与变形长轴的夹角.
假定中线不伸长,得到周向位移:
v ＝－ʏ
w d φ
(１６)法线转角:
μ＝a
r c t a n (１ρd
ρd φ
)(１７
)(４)由图５可见,周向侧隙j t 定义为点M １与
M ２之间的周向距离,
即j
t ＝(x M ２－x a １)２＋(y a
１－y M ２)２
(１８
)２．２㊀基于弧长定位的侧隙算法
文献[３]对基于弧度定位的侧隙算法有全面计算说明,由图４可见,基于弧A o f 等于圆弧
B o １,
有r m φ＝
ʏ
φ１
０
ρ２＋ρ
２
d φ１
(１９
)法线转角
μ＝a
r c t a n ρ
ρ
(２０
)但文献[３
]对以上公式进行了近似计算(见式(４
)),而本文通过数值求解方法来逼近准确的φ１[１４
],
即将这种定位柔轮齿根的侧隙算法称为基于弧长定位的侧隙算法(以下简称弧长法)
.(１
)柔轮齿廓方程[１５
]为
x a １＝r １{s i n [ψ－(
u a １－θ１)]＋u a １c o s α０ ㊀c o s [ψ－(u a １－θ１＋α０)]}＋ρs i n φ１－r m s i n ψ(２１)y a １＝r １{c o s [ψ－(u a １－θ１)]－u a １c o s α０ ㊀s i n [ψ－(u a １－θ１＋α０)]}＋ρc o s φ１－r m c o s ψ
(２２)ψ＝φ１＋μ
(２
)刚轮齿廓方程[１５]为x M ２＝r ２{s i n [φ２－(u M ２－θ２)]＋u M ２ c o s α０c o s [φ２－(
u M ２－θ２＋α０)]}(２３)y M ２＝r ２{c o s [φ２－(
u M ２－θ２)]－u M ２ c o s α０s i n [φ２－(
u M ２－θ２＋α０)]}(２４
)(３
)四滚轮波发生器作用下柔轮中性层的径向位移和法线转角由式(１５)㊁式(１７)确定.(４)由图４可见,周向侧隙j t 定义为点M １与
点M ２之间的周向距离,即
j
t ＝(x M ２－x a １)２＋(y a
１－y M ２)２
(２５
)２．３㊀有限元模型结果验证
本文给出了两个从不同角度计算侧隙的改进方法,下面借助有限元模型对改进后的算法进行验证.
图６为１．２节案例参数下,
四种不同的理论侧隙算法与有限元模型结果得到的周向侧隙比较图.图６显示,位移法为改进文献[２
]的侧隙计算方法,位移法更加吻合有限元模型计算结果
;弧长法为改进文献[３]中的侧隙计算方法,可以看到弧长法也更加吻合有限元模型计算结果;其中弧长法侧隙曲线最优.
图６㊀有限元模型验证结果
F i g
．６㊀T h e v a l i d a t i o n r e s u l t s o f f i n i t e e l e m e n tm o d e l ３㊀柔轮齿根定位方式对周向侧隙的影响
３．１㊀理论算法与有限元模型的周向侧隙比较
图７为啮合区间内不同啮合位置的周向侧隙图.图７显示:在φ＝０ʎ长轴区附近,有限元模型计算结果与两种理论算法计算结果一致.弧长法与有限元模型结果比较:在φɪ(－１０ʎ,０ʎ)区间,弧长法结果偏大,在φɪ(０ʎ,１２ʎ)区间,弧长法结果偏小;当φȡ１２ʎ时,弧长法结果偏大;且φ＝３７ʎ时偏差最大,为０．７２μm ;在φ＝５６ʎ左右侧隙值相
等.位移法与有限元模型计算结果比较如下:在
φɪ(
－１０ʎ,０ʎ)区间,位移法结果偏大;在φɪ(０ʎ,１３ʎ)区间,弧长法结果偏小;在φɪ(１３ʎ,２６ʎ)区间,位移法偏大;在φȡ２６ʎ时,位移法结果偏小;且在φ＝５５ʎ左右偏差最大,为３．１６μm .
图７㊀理论算法与有限元模型的周向侧隙对比图F i g ．７㊀C o m p a r i s o nd i a g r a mo f t h e c i r c u m f e r e n t i a l b a c k l a s hb e t w e e n t h e t h e o r e t i c a l a l g
o r i t h ma n d t h e f i n i t e e l e m e n tm o d e l
３．２㊀柔轮中性层变形与柔轮齿根定位方式比较
文中侧隙的两种理论算法都是基于小变形假定和柔轮中性层不伸长假定.但是研究显示,不同的波发生器作用下,柔轮中性层会有不同程度
３６６ 谐波齿轮的侧隙规律研究与有限元模型仿真
杨朋朋㊀陈晓霞㊀邢静忠等Copyright ©博看网. All Rights Reserved.
的伸长[１６],因此这必然会引起理论算法与有限元模型仿真结果的差异.柔轮中性层的变形和柔轮齿根定位方式是影响侧隙大小的主要因素.图８为理论算法的柔轮中性层的径向位移㊁周向位移和法线转角与有限元模型计算结果的差值s .
图８㊀柔轮齿根变形位置偏差
F i g．８㊀D e v i a t i o no f d e f o r m a t i o n p o s i t i o no f t h e
f l e x s p l i n e t o o t h r o o t
图８显示:在φɪ(０ʎ,４９ʎ)区间,理论算法径向位移偏小;在φɪ(４９ʎ,９０ʎ)区间,径向位移偏大;并在φ＝９０ʎ偏差达到最大.理论算法的周向位移结果偏大,在φ＝５９ʎ左右偏差最大为３．１μm.在φɪ(０ʎ,２８ʎ)区间,理论算法的法线转角偏大;在φɪ(２８ʎ,９０ʎ)区间,法线转角偏小.理论算法在求解径向位移时,以小变形假定为前提,但本算例中的径向位移量偏大,最终也导致周向位移与法线转角出现偏差.
柔轮齿根周向位置极角定义为θv＝φ１－φ,周向位置极角偏差s１为两种理论算法与有限元模型计算的周向位置极角的差,它主要反映柔轮齿根定位方式的不同.图９显示:与有限元模型相比,弧长法的周向位置极角在φɪ(０ʎ,５０ʎ)时偏大,且在φ＝３７ʎ时偏差最大,为３．２８ˑ１０－５r a d;在φɪ(５０ʎ,９０ʎ)时,周向位置极角偏小.位移法的周向位置极角一直偏大,在φ＝５９ʎ左右偏差最大,为２．２４ˑ１０－４r a d.
３．３㊀弧长法周向侧隙分析
弧长法与有限元模型的比较如图７与图９所示,在φɪ(０ʎ,１２ʎ)时,弧长法侧隙结果偏小,是因为周向位置极角偏大,柔轮齿根定位偏右;在φȡ１２ʎ时,弧长法侧隙结果偏大,是因为弧长法的周向位置极角渐渐地出现偏大趋势,即柔轮齿根定位渐渐偏左,且法线转角偏小.图９中φ＝３７ʎ周向位置极角偏差最大为３．２８ˑ１０－５r a d,图７中φ＝３７ʎ时,周向侧隙偏差最大为０．７２μm,由于弧长法计算过程未用到周向位移,故用周向位置极角来解释周向侧隙偏差
.
图９㊀柔轮齿根定位方式偏差
F i g．９㊀D e v i a t i o no f l o c a t i o nm e t h o d s o f t h e
f l e x s p l i n e t o o t h r o o t
３．４㊀位移法周向侧隙分析
位移法与有限元模型的比较如下:在φɪ(０ʎ,１６ʎ)时,由于周向位移偏差和法线转角偏差都比较小,故其两个因素都有所抵消,在这区间内位移法侧隙结果与有限元模型结果虽有偏差,但比较接近.在φȡ１６ʎ时,位移法侧隙结果偏小,这是因为理论计算方法的周向位移比有限元模型计算结果偏大,同时周向位置极角偏大,导致柔轮齿根定位偏右.图９中,在φ＝５９ʎ左右周向位置极角偏差最大为２．２４ˑ１０－４r a d;图８中,在φ＝５９ʎ左右周向位移偏差最大,为３．１μm;图７中,在φ＝５５ʎ左右周向侧隙偏差最大,为３．１６μm.由此可见位移法中周向位移与周向位置极角偏差可同步解释周向侧隙偏差规律.
综上所述,侧隙出现偏差的根本原因仍是未给出准确的柔轮中性层变形位置,即未严格遵守小变形假定,导致径向位移㊁周向位移和法线转角偏差,这会对不同的柔轮齿根定位方式产生影响,进而影响周向侧隙.虽然理论算法与有限元模型侧隙存在偏差,但是在实际谐波齿轮加工生产过程中,以上偏差仅为５级加工精度齿廓总偏差的１/１０[１７].４㊀结论
(１)基于小变形假定的理论,获取了有限元模型柔轮中性层的径向位移㊁周向位移和法线转角,并求解了周向位置极角.与理论计算结果比较发现,理论侧隙与有限元模型计算的侧隙结果的差异主要是由周向位移引起的.
(２)与有限元模型计算结果相比,在柔轮与刚轮将脱离啮合区域,基于弧长定位的侧隙算法计算结果偏大,基于周向位移定位的侧隙算法计算结果偏小.相对来说,基于弧长定位的侧隙算法定位更加准确.因此选择合理的柔轮齿根定位方式可以提高侧隙计算准确性.
４６６
中国机械工程第２９卷第６期２０１８年３月下半月Copyright©博看网. All Rights Reserved.
参考文献:
[１]㊀谢金瑞．国内外谐波传动的应用和发展[J]．光学精密工程,１９７９(４):２２Ｇ３１．
X I E J i n r u i．T h e A p p l i c a t i o n a n d D e v e l o p m e n to f
H a r m o n i cD r i v ea t H o m ea n d A b r o a d[J]．O p t i c s
a n dP r e c i s i o nE n g i n e e r i n g,１９７９(４):２２Ｇ３１．[２]㊀伊万诺夫M H．谐波齿轮传动[M]．沈允文,译．北京:国防工业出版社,１９８７．
I W A N N O V M H．T h e H a r m o n i cD r i v e[M]．S H E N
Y u n w e n,T r a n s．B e i j i n g:D e f e n s e I n d u s t r y P r e s s,１９８７．[３]㊀沈允文,叶庆泰．谐波齿轮传动的理论和设计[M]．北京:机械工业出版社,１９８５．
S H E N Y u n w e n,Y E Q i n g t a i．T h e o r y a n dD e s i g no f
H a r m o n i cD r i v e[M]．B e i j i n g:M e c h a n i c a l I n d u s t r y
P r e s s,１９８５．
[４]㊀辛洪兵．双圆弧谐波齿轮传动基本齿廓设计[J]．中国机械工程,２０１１,２２(６):６５６Ｇ６６２．
X I N H o n g b i n g．D e s i g n f o rB a s i cR a c ko fH a r m o n i c
D r i v ew i t hD o u b l eＧc i r c u l eＧa r cT o o t hP r o f i l e[J]．C h iＧ
n aM e c h a n i c a l E n g i n e e r i n g,２０１１,２２(６):６５６Ｇ６６２．[５]㊀辛洪兵,陈菲．谐波齿轮传动初始啮合侧隙变化规律的研究[J]．试验技术与试验机,２００２,４２(１/２):１３Ｇ１４．
X I N H o n g b i n g,C H E N F e i．S t u d y o nt h eC h a n g e
L a wo f t h e I n i t i a lM e s h i n g B a c k l a s h o f t h eH a r m o nＧ
i cG e a rD r i v e[J]．T e s tT e c h n o l o g y a n dT e s t i n g M aＧ
c h i n e,２００２,４２(１/２):１３Ｇ１４．
[６]㊀殷燕．零侧隙渐开线谐波齿轮传动的参数优化设计及有限元分析[D]．秦皇岛:燕山大学,２０１０．
Y I N Y a n．P a r a m e t e rO p t i m i z a t i o nD e s i g n a n dF i n i t e
E l e m e n tA n a l y s i s o f Z e r oＧl a t e r a l s p a c e I n v o l u t eH a rＧ
m o n i c G e a r D r i v i n g[D]．Q i n h u a n g d a o:Y a n s h a n
U n i v e r s i t y,２０１０．
[７]㊀陈晓霞,林树忠,邢静忠,等．圆弧齿廓谐波齿轮侧隙及干涉检查仿真[J]．计算机集成制造系统,２０１１,１７(３):６４３Ｇ６４８．
C H E N X i a o x i a,L I N S h u z h o n g,X I N GJ i n g z h o n g,
e t a l．S i m u l a t i o no nG e a rB a c k l a s ha n dI n t e r
f e r e n c e
C h e c ko f H a r m o n i c
D r i v e w i t h C i r c u l a rＧa r c T e e t h
P r o f i l e[J]．C o m p u t e r I n t e g r a t e d M a n u f a c t u r i n g
S y s t e m s,２０１１,１７(３):６４３Ｇ６４８．
[８]㊀C H E N X i a o x i a,L I U Y u s h e n g,X I N G J i n g z h o n g,
e t a l．T h eP a r a m e t r i cD e s i g no fD o u b l eＧc i r c u l a rＧa r c
T o o t hP r o f i l ea n dI t sI n f l u e n c eo nt h e F u n c t i o n a l
B a c k l a s h o f H a r m o n i c D r i v e[J]．M e c h a n i s m a n d
M a c h i n eT h e o r y,２０１４,７３(２):１Ｇ２４．[９]㊀刘邓辉,邢静忠,陈晓霞．渐开线谐波齿轮的空间齿廓设计及仿真分析[J]．计算机集成制造系统,２０１５,２１(３):７０９Ｇ７１５．
L I U D e n g h u i,X I N G J i n g z h o n g,C H E N X i a o x i a．S p a t i a lT o o t hP r o f i l eD e s i g na n dS i m u l a t i o nA n a l yＧ
s i so f H a r m o n i cD r i v ew i t hI n v o l u t eT o o t hP r o f i l e
[J]．C o m p u t e rI n t e g r a t e d M a n u f a c t u r i n g S y s t e m s,２０１５,２１(３):７０９Ｇ７１５．
[１０]㊀董惠敏,张晓青．基于实验建模的谐波齿轮传动柔轮的有限元分析研究[J]．机械传动,２００１,２５
(２):１６Ｇ１９．
D O N G H u i m i n,Z HA N G X i a o q i n g．F i n i t e
E l eＧ
m e n tA n a l y s i so f t h eF l e x i b l eo f H a r m o n i cD r i v e
w i t hE x p e r i m e n t a l M o d e l[J]．M e c h a n i c a lT r a n sＧ
m i s s i o n,２００１,２５(２):１６Ｇ１９．
[１１]㊀刘文芝,张乃仁,张春林,等．谐波齿轮传动中杯形柔轮的有限元计算与分析[J]．机械工程学报,
２００６,４２(４):５２Ｇ５７．
L I U W e n z h i,Z HA N G N a i r e n,Z HA N G C h u n l i n,
e t a l．T h eF i n i t eE l e m e n tC a l c u l a t i o na n dA n a l y s i s
o nt h e C o l u m nＧs h a p e d F l e x S p l i n eo f H a r m o n i c
D r i v e[J]．J o u r n a l o f M e c h a n i c a l
E n g i n e e r i n g,
２００６,４２(４):５２Ｇ５７．
[１２]㊀O S T A P S K IW,MU K HAI．S t r e s sS t a t eA n a l y s i s o fH a r m o n i cD r i v eE l e m e n t sb y F E M[J]．B u l l e t i n
o f t h eP o l i s hA c a d e m y o fS c i e n c e s:T e c h n i c a l S c iＧ
e n c e s,２００７,５５(１):１１５Ｇ１２３．
[１３]㊀付军峰,董海军,沈允文．谐波齿轮传动中柔轮应力的有限元分析[J]．中国机械工程,２００７,１８
(１８):２２１０Ｇ２２１４．
F UJ u n f e n g,D O N
G
H a i j u n,S H E N Y u n w e n．T h e
F i n i t eE l e m e n tA n a l y s i s o f t h e S t r e s s o f F l e x i b l e i n
H a r m o n i cD r i v e[J]．C h i n a M e c h a n i c a lE n g i n e e rＧ
i n g,２００７,１８(１８):２２１０Ｇ２２１４．
[１４]㊀陈晓霞,林树忠,邢静忠．谐波齿轮传动中基于柔轮装配变形的共轭精确算法[J]．中国机械工程,
２０１０,２１(１７):２０５３Ｇ２０５７．
C H E N X i a o x i a,L I NS h u z h o n g,X I N GJ i n g z h o n g．
E x a c t C o n j u g a t e A l g o r i t h m B a s e d o n A s s e m b l y
D e f o r m a t i o no fF l e x s p l i n e i n H a r m o n i cD r i v e s[J]．
C h i n a M e c h a n i c a l E n g i n e e r i n g,２０１０,２１(１７):
２０５３Ｇ２０５７．
[１５]㊀沈允文．渐开线谐波齿轮传动的几何计算[J]．齿轮,１９８６,１０(６):５０Ｇ５５．
S H E N Y u n w e n．G e o m e t r i cC a l c u l a t i o no f I n v o l u t e
H a r m o n i cG e a rD r i v e[J]．G e a r,１９８６,１０(６):５０Ｇ
５５．
[１６]㊀陈晓霞,刘玉生,邢静忠,等．谐波齿轮中柔轮中性层的伸缩变形规律[J]．机械工程学报,２０１４,
５０(２１):１８９Ｇ１９６．
C H E N X i a o x i a,L I U Y u s h e n g,X I N GJ i n g z h o n g,
e t a l．N e u t r a lL i n eS t r e t c ho fF l e x s p l i n ei n H a rＧ
m o n i cD r i v e r[J]．J o u r n a l o fM e c h a n i c a lE n g i n e e rＧ
i n g,２０１４,５０(２１):１８９Ｇ１９６．
５６６
谐波齿轮的侧隙规律研究与有限元模型仿真 杨朋朋㊀陈晓霞㊀邢静忠等
Copyright©博看网. All Rights Reserved.。