浙江大学学报(工学版), 2019, 53(12): 2289-2297 doi: 10.3785/j.issn.1008-973X.2019.12.005

机械与能源工程

制造误差对谐波齿轮应力的影响规律

陈茜,, 李俊阳,, 王家序, 蒋倩倩, 唐挺

Influence law of manufacturing error on harmonic gear stress

CHEN Qian,, LI Jun-yang,, WANG Jia-xu, JIANG Qian-qian, TANG Ting

通讯作者: 李俊阳,男,副教授. orcid.org/0000-0001-8317-6228. E-mail: lijunyang1982@sina.com

收稿日期: 2019-03-2  

Received: 2019-03-2  

作者简介 About authors

陈茜(1995—),女,硕士生,从事谐波减速器齿廓研究.orcid.org/0000-0002-1353-2378.E-mail:957772748@qq.com , E-mail:957772748@qq.com

摘要

为研究谐波齿轮各种制造误差对柔轮齿面应力与磨损情况的影响,根据响应面思想,设计实验,采用四因素响应面模型,分析柔轮跨距值(M值)偏差、刚轮M值偏差、波发生器长半轴偏差及短半轴偏差,分别建立相应谐波齿轮模型,并进行有限元仿真分析. 实验结果表明,凸轮长半轴偏差对柔轮齿上应力敏感度最大,柔轮、刚轮M值偏差次之,凸轮短半轴偏差敏感度最小;针对在制造加工时可能出现的刚轮与柔轮齿廓不同偏差情况提出4种误差补偿方案,即当柔轮、刚轮M值均为负偏差或分别为正偏差、负偏差时减小凸轮长半轴,当柔轮、刚轮M值均为正偏差或分别为负偏差、正偏差时增大凸轮长半轴,以此改善啮合情况并提高精度保持性.

关键词: 谐波齿轮传动 ; 双圆弧齿廓 ; 制造误差 ; 有限元 ; 响应面法

Abstract

The experiment was designed based on the response surface theory to study the influence about various manufacturing errors of harmonic gears on the tooth surface stress and wear of the flexspline. A four-factor response surface model was adopted, which includes the span value (M-value) deviation of the flexspine, the M-value deviation of the circular spline, the long half axis deviation, and the short half axis deviation of the wave generator. Then corresponding harmonic gear models were established for finite element simulation analysis, respectively. The experimental results show that the cam long half axis deviation is most sensitive to the stress on the flexspline tooth, the M-value deviation of the flexspline and the M-value deviation of the circular spline less, and the cam short axis deviation sensitivity the least. Four kinds of error compensation schemes were proposed for different deviations of the gears of the circular spline and the flexspline, which may occur during the manufacturing process, to improve the meshing condition and the accuracy retention. That is, the long half axis of the cam can be reduced when the M-values of the flexspline and the circular spline are both negative deviation, or positive and negative deviation respectively; the long half axis of the cam can be increased when the M-values of the flexspline and the circular spline are both positive deviation, or negative and positive deviation respectively.

Keywords: harmonic gear transmission ; double-circular-arc tooth profile ; manufacturing error ; finite element ; response surface method

PDF (1467KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

陈茜, 李俊阳, 王家序, 蒋倩倩, 唐挺. 制造误差对谐波齿轮应力的影响规律. 浙江大学学报(工学版)[J], 2019, 53(12): 2289-2297 doi:10.3785/j.issn.1008-973X.2019.12.005

CHEN Qian, LI Jun-yang, WANG Jia-xu, JIANG Qian-qian, TANG Ting. Influence law of manufacturing error on harmonic gear stress. Journal of Zhejiang University(Engineering Science)[J], 2019, 53(12): 2289-2297 doi:10.3785/j.issn.1008-973X.2019.12.005

随着航空航天[1-2]、机器人[3]、医疗器械、交通运输等领域的快速发展,其中的关键部件之一谐波减速器逐渐成为近年来的研究重点,人们也对其提出了更严苛的要求. 区别于传统刚性构件的传动方式,谐波传动通过构建可控制弹性变形来实现机械运动的传递,具有结构紧凑、传动平稳、效率高、轻量化等优点. 由柔轮、刚轮、波发生器三大部件组成的谐波减速器最主要的失效形式是柔轮的疲劳断裂[4],因此,柔轮的应力与变形规律是当前研究重点. 为减小啮合应力,增加柔轮疲劳寿命,许多学者针对齿廓形状参数[5-6]、结构参数[7-8]等因素进行了大量研究. 王家序等[9]提出的三维双圆弧空间齿廓,提高了啮合性能. 杨勇等[10]针对双圆弧谐波齿廓传动,提出一种增大“双共轭”啮合区间的求解方法,优化了齿廓参数,提高了扭转刚度及传动精度. 随着计算机技术的快速提高,有限元法逐渐成为研究谐波传动应力变化的主要手段. 李奇[11]采用有限元法,分析柔轮结构不同参数,发现柔轮径向变形对柔轮应力以及疲劳寿命影响最大,壁厚和长径比次之. 张雷等[12]基于有限元法与响应面法,综合分析了柔轮圆角、凸轮结构等参数对柔轮刚度和内壁、齿圈等的应力影响规律,但是未考虑齿形的影响. 然而,在实际加工过程中提高齿廓的加工制造精度是减小啮合应力、提高啮合性能的最主要手段. 刚轮与柔轮的齿厚偏差、凸轮长短轴偏差是能够通过检测来量化的4个重要参数. 但是,如何通过优化上述可检验项目上、下偏差取值来控制柔轮的应力,进行实验探究及定量分析,目前仍然没有文献涉及.

本文以无公切线式双圆弧齿廓为研究对象,考虑引起谐波传动齿廓变化的三大重要构件制造误差,包括刚轮与柔轮的齿厚误差和凸轮的长短轴加工误差,通过响应面法与有限元法揭示谐波齿轮制造误差对柔轮应力分布情况的影响规律,以期得到生产制造中对柔轮应力影响最小的误差范围,减少因高加工精度带来的高成本.

1. 双圆弧谐波传动齿廓设计

由于双圆弧齿廓有更大的啮合弧长,目前谐波传动齿廓研究热点为公切线双圆弧齿廓[13-14]. 本次研究采取的无公切线双圆弧齿廓相较于公切线双圆弧齿廓,其共轭区间更大且设计流程更简单,更便于生产加工. 如图1所示为局部坐标系下的双圆弧基本齿形. 该柔轮局部坐标系S1是以柔轮轮齿的对称轴为Y1轴,以Y1轴和柔轮中性层曲线的交点O1为坐标原点,以过O1点垂直于Y1轴的直线为X1轴. 该无公切线双圆弧齿廓包括三段相切圆弧,分别为柔轮齿凸齿廓AB、凹齿廓BC、齿根圆过渡齿廓CD. 如表1所示为双圆弧齿廓的基本参数含义.

图 1

图 1   双圆弧齿廓坐标系

Fig.1   Tooth profile of flexspline with double-circular-arc tooth profile


表 1   双圆弧齿廓参数

Tab.1  Parameters of double-circular-arc tooth profile

符号 意义 符号 意义
${h^{\rm{*}}_{\rm{a}}}$ 齿顶高系数 $\gamma $ 公切线倾角
${h^{\rm{*}}_{\rm{f}}}$ 齿根高系数 $t$ 齿根壁厚
$h$ 全齿高 ${X_{\rm{a}}}$ 凸齿圆心移距量
${\rho _{\rm{a}}}$ 凸圆弧齿廓半径 ${Y_{\rm{a}}}$ 凸齿圆心偏移量
${\rho _{\rm{f}}}$ 凹圆弧齿廓半径 ${X_{\rm{f}}}$ 凹齿圆心移距量
${Y_{\rm{f}}}$ 凹齿圆心偏移量 m 柔轮模数
kt 齿厚比

新窗口打开| 下载CSV


根据图1,以从齿顶A处起的齿廓弧长u为参数,r为右侧齿廓的矢径, $ {n}$ 为其对应法向量.

1)在AB段上:

$ \begin{split} & {{{r}}_1} = \left[ {\begin{array}{*{20}{c}} {{\rho _{\rm{a}}}\cos \;(\theta - {u / {{\rho _{\rm{a}}}}}) + {x_{{\rm{oa}}}}} \\ {{\rho _{\rm{a}}}\sin \;(\theta - {u / {{\rho _{\rm{a}}}}}) + {y_{{\rm{oa}}}}} \\ 1 \end{array}} \right], \\ & {{{n}}_1} = \left[ {\begin{array}{*{20}{c}} {\cos \;(\theta - {u / {{\rho _{\rm{a}}}}})} \\ {\sin \;(\theta - {u / {{\rho _{\rm{a}}}}})} \\ 1 \end{array}} \right]. \end{split} $

式中:

$\left.\begin{split} & u \in (0,{\rm{ }}{l_1}),\;{x_{{\rm{oa}}}} = - {l_{\rm{a}}}, \\ & {y_{{\rm{oa}}}} = {h_{\rm{ }}}_{\rm{f}} + t - {X_{\rm{a}}}, \\ & {l_1} = {\rho _{\rm{a}}}(\theta - \gamma ), \\ & \theta = \arcsin \left[ {{{({h_{\rm{a}}} + {X_{\rm{a}}})} / {{\rho _{\rm{a}}}}}} \right]. \end{split} \right\}$

2)在BC段上:

$ \begin{split} & {{{r}}_2} = \left[ {\begin{array}{*{20}{c}} {{x_{{\rm{of}}}} - {\rho _{\rm{f}}}\cos \;[\gamma + ({{u - {l_{\rm{1}}}}})/{{{\rho _{\rm{f}}}}}]} \\ {{y_{{\rm{of}}}} - {\rho _{\rm{f}}}\sin \;[\gamma + ({{u - {l_{\rm{1}}}}})/{{{\rho _{\rm{f}}}}}]} \\ 1 \end{array}} \right], \\ & {{{n}}_2} = \left[ {\begin{array}{*{20}{c}} {\cos\; [\gamma + ({{u - {l_1}}})/{{{\rho _{\rm{f}}}}}]} \\ {\sin\; [\gamma + ({{u - {l_1}}})/{{{\rho _{\rm{f}}}}}]} \\ 1 \end{array}} \right]. \end{split} $

式中:

$\left.\begin{split} & u \in ({l_1},{\rm{ }}{l_2}), \\ & {x_{{\rm{of}}}} = ({\rho _{\rm{a}}} + {\rho _{\rm{f}}})\cos \gamma + {h_{\rm{l}}}\tan \gamma - {l_{\rm{a}}}, \\ & {y_{{\rm{of}}}} = {h_{\rm{f}}} + t + {X_{\rm{f}}}, \\ & {l_2} = {l_1} + {\rho _{\rm{f}}}\left\{ {{\rm{arcsin}}\left[ {{{({X_{\rm{f}}} + {h_{\rm{f}}})} /{{\rho _{\rm{f}}}}}} \right] - \gamma } \right\}. \end{split} \right\}$

首先建立啮合基本方程并采用包络法设计刚轮的理论共轭齿廓. 采用数值离散思想,将各段的弧长参数u离散成s个点,再将对应齿廓段点的rn代入啮合基本方程,即可求得该点刚柔轮共轭转角集,最后通过坐标变换求得与柔轮齿廓共轭的刚轮理论齿廓.

2. 谐波齿轮制造误差分析

2.1. 齿厚偏差

按照谐波减速器的建构组成,将制造误差分为柔轮、刚轮和凸轮制造误差. 柔轮、刚轮的制造误差主要以传动侧隙的大小影响谐波齿轮传动工作性能,为本实验主要研究对象. 侧隙较大会产生齿间冲击,影响齿轮传动的平稳性. 一般传递运动的精密谐波齿轮传动要求其最小侧隙为零. 侧隙是在谐波齿轮传动中柔轮与刚轮啮合处的间隙,其大小以圆周侧隙jt表示,设K1XK1,YK1)为柔轮齿廓上某点的坐标,以rK 为半径作圆弧与相邻刚轮齿廓相交,可得到对应刚轮齿廓上的交点K2(XK2,YK2),则有

$\left.\begin{split} & {j_{\rm{t}}} = {\left[ {{{\left( {{X_{K2}} - {X_{K1}}} \right)}^2} + {{\left( {{Y_{K1}} - {Y_{K2}}} \right)}^2}} \right]^{{\rm{1/2}}}}, \\ & {r_K} = {\left( {X_{K1}^2 + Y_{K1}^2} \right)^{{\rm{1/2}}}}. \end{split} \right\}$

当谐波齿轮承受负载时,柔轮产生弹性扭转变形,并使轮齿产生附加的切向位移,如果侧隙过小,将会引起变形干涉,故将柔轮因弹性变形而引起的侧隙减小量je作为分析因素之一,计算公式为

${j_{\rm{e}}} \approx {{Tb}}/({{d_1^2\delta G}}).$

式中:T为作用的扭矩,b为柔轮齿宽;d1为柔轮分度圆直径,δ为柔轮壁厚,G为柔轮剪切模量.

传动侧隙的大小取决于齿厚偏差ΔEsn,其指分度圆柱面上法向齿厚实际值与公称值之差. Esns为齿厚上偏差,Esni为齿厚下偏差. 如图2所示为齿厚偏差示意图. 在工程实际中一般采用测量跨距值(M值)的方式来度量齿厚偏差. 对于圆弧齿轮,量柱与齿面接触为一直线,测量方便. 该方法计算简便、测量准确、数值稳定,且不需要测量基准面,消除了齿顶圆误差、孔径误差和齿圈径向跳动的影响,适用于谐波齿轮的测量.

图 2

图 2   齿厚偏差示意图

Fig.2   Diagram of tooth thickness deviation


2.1.1. 柔轮理论M值的计算原理

柔轮为外啮合齿轮,设点E为柔轮齿凸齿廓的圆心,其坐标为E(xEyE),点F为量柱的圆心,点H为点EY轴上的投影. 量柱与凸齿廓接触相切于点G,显然点EFG在同一条直线上. 如图3所示为柔轮齿廓理论M值计算原理图,图中,z为柔轮齿数,α为柔轮半齿角,dp为量柱直径,dM为测量跨距值时的实际读数(量柱直径与跨距值之和),则有

图 3

图 3   柔轮齿廓理论M值计算原理图

Fig.3   Schematic diagram of theoretical M-value for tooth profile of flexspline


$ M = {d_{\rm{M}}} + {d_{\rm{p}}},\;\angle FOG{\rm{ = }}{\text{π}} /z,\;\angle EOH = \alpha \tan \;({{{x_E}} / {{y_E}}}); $

${l_{OE}} = \sqrt {x_E^2 + y_E^2} ;$

${l_{EF}} = {l_{FG}} + {l_{EH}} = {d_{\rm{p}}}/2 + {\rho _{\rm{f}}}.$

由三角形余弦定理可知:

${{l^2_{EF}}} = {l^2_{OE}} + {l^2_{OF}} - 2 \times {l_{OE}} \times {l_{OF}} \times \cos\; (\angle EOF).$

由式(7) ~ (10)可得柔轮齿廓的理论M值.

2.1.2. 刚轮理论M值的计算原理

刚轮为内啮合齿轮,选择不同直径大小的量柱,量柱会与刚轮不同的齿廓段相切.

1)相切于凸齿廓的刚轮理论M值的计算.

当选择较大直径的量柱时,量柱与刚轮凸齿廓相切,如图4(a)所示为相切于凸齿廓的刚轮齿廓理论跨距值计算原理图,则有

图 4

图 4   刚轮齿廓理论M值计算原理图

Fig.4   Schematic diagram of theoretical M-value for tooth profile of circular spline


$ \angle EOF{\rm{ = }}\alpha \tan \;({{{x_E}} / {{y_E}}}),\;M = {d_{\rm{M}}} - {d_{\rm{p}}}. $

由式(8)~(11)可求得lOF值,从而得到相切于凸齿廓的刚轮齿廓理论M值.

2)相切于凹齿廓的刚轮理论M值的计算.

选择较小直径的量柱,则会与刚轮相切于凹齿齿廓,如图4(b)所示为相切于凹齿廓的刚轮齿廓理论M值计算原理图,则有

${l_{EF}} = {l_{EG}} - {l_{FG}} = {\rho _{\rm{f}}} - {d_{\rm{p}}}/2.$

由式(8)、(10)~(12),可得到相切于凹齿廓的刚轮齿廓理论M值. 通过比较理论M值和实际测量M值,即可确定制造精度.

2.2. 凸轮长短轴偏差

波发生器迫使柔轮作弹性波动变形,决定了谐波齿轮传动的运动原始曲线,故波发生器的制造误差对谐波传动的精度影响显著. 凸轮长短轴偏差对标准椭圆凸轮制造误差影响显著,直接影响柔轮径向变形量,改变了柔轮筒体张角大小,进而影响谐波齿轮应力状态,因此将凸轮长短轴偏差作为分析因素之一.

设椭圆波发生器长短轴半径分别为ab,柔性轴承厚度为t′,波发生器最大径向变形量为w0′,柔轮内径为de,则上述参数存在以下关系:

$\left. \begin{array}{l} a = 1/2{d_{\rm{e}}} - {t'} + {w_0}', \\ b = 1/2{d_{\rm{e}}} - {t'} - {w_0}'. \\ \end{array} \right\}$

3. 采用响应面法分析制造误差对柔轮应力变化规律的影响

3.1. 响应面法的外切中心复合设计

响应面分析法(response surface methodology, RSM)通过选取合适的因素水平,设计较少的实验次数并拟合获得因素与响应值之间良好的函数关系,是寻求最优参数的一种解决多因素影响问题的数学统计方法. 中心复合设计(central composite designs, CCDs)因其以最小实验次数提供较多关于实验因素与误差的信息与更好的灵活性等优点,在工程、工业等领域应用颇广. CCDs拟合的二阶响应曲面模型[15]

$y = {\beta _0} + \sum\limits_{i = 1}^k {{\beta _i}{x_i}} + \sum\limits_{i<j}^{} {\sum\limits_{}^n {{\beta _{ij}}{x_i}{x_j} + } } \sum\limits_{i = 1}^k {{\beta _{ii}}{x_i}^2} + \varepsilon .$

式中:y为响应变量,β0为常量,βiβijβii为各项系数,xi为各项自变量,ε为误差项. 模型中共有1+2k+k(k−1)/2个参数,由于k≥1,至少有3个参数,即因素至少要有3个水平.

CCDs又分为外切中心复合设计(central composite circumscribed designs,CCC)、嵌套中心复合设计(central composite inscribed designs,CCI)、面心立方设计(central composite face-centred designs,CCF). 有以下关系式:

$\eta = \sqrt[4]{F},\;F = {2^k}.$

式中:k为因素个数,η为轴向距离参数,F为析因设计点的数目.

对3种模型的拟合参数比较发现,CCC估计的一致精度最高,模型参数估计最有效且具有可旋转性,故选择采用CCC模型进行本次实验设计.

3.2. 基于外切中心复合设计的多因素多水平实验

为研究谐波减速器各个制造误差对柔轮齿上应力的影响,拟合出各制造误差与柔轮齿上观测点应力的函数关系模型,通过Design Expert响应面法实验设计软件设计多因素多水平实验. 其中,研究因素如下:柔轮M值偏差δ1,刚轮M值偏差δ2,凸轮长半轴偏差δ3,凸轮短半轴偏差δ4,共4个因素,即k=4,通过式(15)得到η=2,F=16.

由于目前没有可查阅的双圆弧柔性齿轮精度等级标准,本实验参考小模数渐开线圆柱齿轮精度及参数表,柔轮采用六级制造精度,刚轮采用七级精度. 综合采集到工厂生产的CSF25-120谐波减速器样机制造误差数据,计算出样机的刚轮与柔轮M值,以及凸轮长短半轴误差取值范围. 如表2所示为最终实验的因素范围与水平. 通过对不同因素水平合理组合设计得到30组实验.

表 2   谐波齿轮制造误差因素水平表

Tab.2  Factors and levels of harmonic gear manufacturing error

水平编码 δ1/mm δ2/mm δ3/mm δ4/mm
−2 −0.050 −0.050 −0.050 −0.050
−1 −0.025 −0.025 −0.025 −0.025
0 0 0 0 0
1 0.025 0.025 0.025 0.025
2 0.050 0.050 0.050 0.050

新窗口打开| 下载CSV


3.3. 谐波齿轮传动响应面模型建立

3.3.1. 谐波齿轮传动有限元分析

由于样机制作成本较高,且有限元法与理论计算结果较为一致[16],目前有限元分析(finite element analysis,FEA)被广泛用于研究谐波齿轮传动的柔轮应力研究.

本次实验建模以哈默纳科(Harmonic Drive,HD)公司的CSF-25-120为实例,传动比i=120,柔轮齿数z1=240,刚轮齿数z2=242,柔轮内壁直径de=61.32 mm. 采用压平设计的标准椭圆凸轮,其理论径向变形量w0'=0.261 mm,计算出刚轮理论M值为63.354 9 mm,柔轮理论M值为64.081 7 mm.

建立刚轮、柔轮、波发生器的三维模型,进行网格划分,为提高计算精确性与收敛性,采用八节点六面体网格. 对刚轮与柔轮轮齿细化处理,如图5所示为谐波齿轮有限元模型网格划分. 如表3所示为刚轮、柔轮及波发生器的材料性能参数. 表3中,E为弹性模量,μ为泊松比,ρ为密度.

图 5

图 5   谐波齿轮有限元分析(FEA)模型网格划分

Fig.5   Meshing of harmonic gear finite element analysis (FEA) model


表 3   谐波齿轮材料物理性能参数

Tab.3  Material physical performance parameters of harmonic gear

部件 材料 E/GPa μ ρ/(kg·m−3)
柔轮 30CrMnSiA 196 0.3 7 750
刚轮 45 210 0.269 7 850
波发生器 45 210 0.269 7 850

新窗口打开| 下载CSV


为减少工作量并缩短求解时间,对模型作部分结构上的简化. 将波发生器简化为椭圆刚体,去掉柔轮杯底螺栓孔,刚轮只取有齿部分并简化为内圆柱齿圈. 因加载对柔轮齿应力影响很小[16],不考虑加载,设置以下2个分析步骤.

1)在柔轮法兰盘面建立耦合点,对该点加上固定约束,波发生器椭圆分为上、下两半,分别设置上、下位移约束,使波发生器外表面分别与柔轮内壁进行挤压接触,将其撑开.

2)在刚轮外圆表面建立耦合点,并对该点添加位移约束,将刚轮装进柔轮,与柔轮装配,刚轮与柔轮齿廓进行接触.

图6所示为以上2个分析步骤完成后的柔轮齿圈应力云. 由图6可知,柔轮齿上最大应力位置在啮入区附近的齿圈处.

图 6

图 6   柔轮齿圈应力云图

Fig.6   Stress nephogram of flexpline ring gear


3.3.2. 响应面模型的建立

采集有限元仿真实验各组数据,基于响应面思想拟合得到柔轮M值偏差δ1、刚轮M值偏差δ2、凸轮长半轴偏差δ3、凸轮短半轴偏差δ4四因素对于柔轮齿上固定观测点应力σ的响应面模型如下:

对该方程进行回归显著性检验,P值大于F值的概率小于0.001,表明该模型具有统计学意义,显著性表现为显著. 对各因素进行敏感度分析,一次项系数显示各因素的重要程度为δ3>δ2>δ1>δ4,二次项系数显示各因素的重要程度为δ3>δ2>δ4>δ1. 综合分析可知,δ3对该模型的影响程度最大,δ2δ1次之,δ4的影响程度最小. 即凸轮长半轴偏差对柔轮齿上应力敏感度最大,刚轮M值偏差、柔轮M值偏差次之,凸轮短半轴偏差敏感度最小. 因此在实际生产过程中,为达到产品高性能及低成本的目标,需优先保证凸轮长半轴加工精度,再保证其他因素的精度.

为探究各因素对柔轮应力的具体影响规律,根据得到的数学模型,对上述4个因素分别进行定量分析.

3.4. 基于响应面模型的单因素对柔轮应力影响
3.4.1. 柔轮M值偏差对柔轮齿上应力的影响

δ2=0.03 mm、δ3=0、δ4=0时,根据实验数据结果,绘制出的柔轮M值偏差δ1对柔轮齿上观测点应力σ的影响曲线如图7所示. 根据图7,在给定范围内,柔轮齿上应力随柔轮M值偏差增大递增. 这是由于柔轮M值偏大使柔轮与刚轮啮合位置不正常,最大啮合深度增加,齿根发生过度啮合致使应力增大. 此外柔轮挠度随壁厚增大而降低,最终应力表现为增加. 根据柔轮疲劳强度评价准则体系,应力在300~1500 MPa为正常啮合范围,则柔轮M值偏差合理范围(−0.01 mm,0 mm).

图 7

图 7   柔轮M值偏差对柔轮齿上观测点应力影响曲线

Fig.7   Effect of M-value deviation of flexspline on stress of flexspline tooth observation point


3.4.2. 刚轮M值偏差对柔轮齿上应力的影响

δ1=−0.03 mm、δ3=0、δ4=0时,如图8所示为刚轮M值偏差δ2对柔轮齿上应力影响曲线. 根据图8,柔轮齿上应力随刚轮M值偏差增大递减. 这是因为当刚轮M值偏大使其齿廓外移,接触齿面随啮合深度减小而减小,导致柔轮轮齿上应力递减. 刚轮M值偏差的合理范围为(0.005 mm,0.030 mm).

图 8

图 8   刚轮M值偏差对柔轮齿上观测点应力影响曲线

Fig.8   Effect of M-value deviation of circular spline on stress of flexspline tooth observation point


3.4.3. 凸轮长半轴偏差与短半轴偏差对柔轮齿上应力的影响

1)当δ1=−0.03 mm、δ2=0、δ4=0时,如图9所示为绘制出的凸轮长半轴偏差δ3对柔轮齿上应力影响曲线. 由图9知,凸轮长半轴偏差对柔轮齿上应力影响显著,柔轮齿上应力随凸轮长半轴偏差增大递增. 该曲线基于数学模型推导得到,故曲线中出现了应力为负的情况. 在工程实际中不考虑此部分的取值. 凸轮长半轴偏差最佳范围为(−0.04 mm,−0.02 mm).

图 9

图 9   凸轮长半轴偏差对柔轮齿上观测点应力影响曲线

Fig.9   Effect of cam long half axis deviation on stress of flexspline tooth observation point


2)当δ1=0、δ2=−0.03 mm、δ3=−0.05 mm时,如图10所示为凸轮短半轴偏差δ4对柔轮齿上应力影响曲线. 根据图10,当δ4∈(−0.05 mm,−0.02 mm)时,柔轮齿上应力随δ4增大略减;当δ4∈(−0.02 mm,0.05 mm)时,柔轮齿上应力随δ4增大而略增. 故制造加工凸轮时可控制凸轮短半轴偏差在(−0.05 mm,0.05 mm),保证在较低加工成本下得到良好啮合与较小应力波动.

图 10

图 10   凸轮短半轴偏差对柔轮齿上观测点应力影响曲线

Fig.10   Effect of cam short half axis deviation on stress of flexspline tooth observation point


波发生器的长短轴半径是按理论假设的线性关系推导出波发生器内侧与柔轮内壁接触处的最大径向变形量再计算得到,而研究表明,理论设计的波发生器在装配后未能使柔轮达到理论变形效果,柔轮沿轴向的径向变形量不符合线性关系,会造成啮入区干涉,引起应力集中[9]. 故凸轮长短半轴和柔轮最大径向变形量密切相关,由式(13)可得凸轮径向变形量w0′=(ab)/2,凸轮长半轴增大或短半轴减小都会使实际径向变形量超过理论径向变形量,实际啮合深度增大. 虽然提高了承载能力[16],但柔轮齿上应力剧增. 本实验按柔轮理论最大径向变形量设计凸轮理论长短轴半径且未修型,故应力偏大.

综上分析可知,柔轮应力随柔轮M值减小而减小,随刚轮M值增大而减小. 如图11所示为δ1δ2对柔轮齿上应力的三维响应;当柔轮M值偏大时,减小凸轮长半轴可使柔轮齿上应力不剧增;当柔轮M值偏小时,增加凸轮长半轴可使柔轮齿上应力不骤减. 如图12所示为δ1δ3对柔轮齿上应力的三维响应;此外,减小凸轮长半轴可减小刚轮M值偏小时增大的柔轮应力. 如图13所示为δ2δ3对柔轮齿上应力的三维响应.

图 11

图 11   柔轮与刚轮M值偏差对柔轮齿上观测点应力的响应面3D图

Fig.11   Response surface of M-value deviation of flexspline and circular spline on stress of flexspline tooth observation point


图 12

图 12   柔轮M值偏差与凸轮长半轴偏差偏差对柔轮齿上观测点应力的响应面3D图

Fig.12   Response surface of M-value deviation of flexspline and cam long half axis deviation on stress of flexspline tooth observation point


图 13

图 13   刚轮M值偏差与凸轮长半轴偏差偏差对柔轮齿上观测点应力的响应面3D图

Fig.13   Response surface of M-value deviation of circular spline and cam long half axis deviation on stress of flexspline tooth observation point


3.5. 不同参数制造误差补偿方案设计

柔轮的传动精度决定着谐波减速器的传动精度. 在谐波减速器各部件的生产加工过程中,无可避免会出现刚轮、柔轮齿廓加工误差,从而引起柔轮齿上应力集中,影响啮合性能. 一般来说,刚轮与柔轮齿加工难度及生产成本较高,而凸轮形状规则,加工难度和生产成本较低. 故在刚轮、柔轮已完成生产加工后,可通过调整凸轮来控制刚轮与柔轮齿廓加工误差导致的应力变化. 现根据实际生产中四种不同情形,设计出以下4种应力补偿方案. 如图14所示为凸轮的4种调整方案曲线. 图中,Δ1Δ2Δ3Δ4为4种方案凸轮长半轴的调整量.

图 14

图 14   4种方案的凸轮长半轴调整曲线

Fig.14   Adjustment curves of cam long half axis by four kinds of solutions


1)方案一:当柔轮M值、刚轮M值为负偏差,可采用减小凸轮长半轴的方式以减小柔轮应力. 例如当柔轮M值偏差δ1=−0.02mm,刚轮M值偏差δ2=−0.02 mm,柔轮齿上应力σ与凸轮长半轴调整量Δ1关系曲线如图14(a)所示.

2)方案二:当柔轮M值偏差、刚轮M值偏差为正偏差,可增大凸轮长半轴以改善啮合不足. 例如当柔轮M值偏差δ1=0.02 mm,刚轮M值偏差δ2=0.05 mm,柔轮齿上应力σ与凸轮长半轴调整量Δ2关系曲线如图14(b)所示.

3)方案三:当柔轮M值偏差为负偏差,刚轮M值偏差为正偏差,可采取增大凸轮长半轴的方式来改善啮合不足。例如柔轮M值偏差δ1=−0.03 mm,刚轮M值偏差δ2=0.03 mm,柔轮齿上应力σ与凸轮长半轴调整量Δ3关系曲线如图14(c)所示。

4)方案四:当柔轮M值偏差为正偏差,刚轮M值偏差为负偏差,可采取减小凸轮长半轴的方式以减小柔轮应力。例如当柔轮M值偏差δ1=0.02 mm,刚轮M值偏差δ2=−0.02 mm,柔轮齿上观测点应力σ与凸轮长半轴调整量Δ4关系曲线如图14(d)所示。

4. 结 论

(1)在刚轮与凸轮的共同作用下,柔轮上最大应力出现在柔轮齿面啮入区位置处.

(2)凸轮长半轴偏差对柔轮应力影响最大,柔轮M值偏差、刚轮M值偏差次之,凸轮短半轴偏差最小. 若凸轮长半轴或柔轮M值增大,柔轮齿上应力随之增大;若刚轮M值增大,柔轮齿上应力减小;柔轮齿上应力随凸轮短半轴先减小后增大,但变化很小.

(3)在谐波减速器的加工制造过程中,当刚轮与柔轮出现齿廓偏差时,可通过控制凸轮长半轴偏差来使谐波减速器轮齿啮合情况不受较大影响,改善啮合情况,提高精度保持性,降低生产成本和废品率.

参考文献

丰飞, 王炜, 唐丽娜, 等

空间高精度谐波减速器的应用及其发展趋势

[J]. 机械传动, 2014, 38 (10): 98- 107

[本文引用: 1]

FENG Fei, WANG Wei, TANG Li-na, et al

Application and development trends of the space harmonic reducer with high precision

[J]. Journal of Mechanical Transmission, 2014, 38 (10): 98- 107

[本文引用: 1]

陶孟仑, 陈阳鹏, 陈定方, 等

谐波减速器测试技术研究现状及展望

[J]. 机械传动, 2018, 42 (7): 175- 180

[本文引用: 1]

TAO Meng-lun, CHEN Yang-peng, CHEN Ding-fang, et al

Research present status and outlook of harmonic reducer testing technology

[J]. Journal of Mechanical Transmission, 2018, 42 (7): 175- 180

[本文引用: 1]

付小月, 王伟, 王雷, 等

变刚度关节驱动器动力学特性的分析与研究

[J]. 机器人, 2017, 39 (4): 466- 473

[本文引用: 1]

FU Xiao-yue, WANG Wei, WANG Lei, et al

Analysis and study on dynamical characteristics of variable stiffness joint actuator

[J]. Robot, 2017, 39 (4): 466- 473

[本文引用: 1]

夏田, 杨世勇, 王亚峰

谐波减速器柔轮疲劳分析

[J]. 机械传动, 2017, 41 (10): 133- 135

[本文引用: 1]

XIA Tian, YANG Shi-yong, WANG Ya-feng

Fatigue analysis of the harmonic reducer flexspline

[J]. Journal of Mechanical Transmission, 2017, 41 (10): 133- 135

[本文引用: 1]

DENNIS L, NELSON A, ANDRE’S T

Statistical analysis of the influence of tooth geometry in the performance of a harmonic drive

[J]. Journal of the Brazilian Society of Mechanical Sciences and Engineering, 2015, 37: 723

DOI:10.1007/s40430-014-0197-0      [本文引用: 1]

ISHIKAWA S, KIYOSAWA Y. Flexing contact type gear drive of non-profile-shifted two-circular-arc composite tooth profile: United States, US 5458023[P]. 1995–10-17.

[本文引用: 1]

赵建虎. 谐波传动机构柔轮的应力分布及寿命特性分析[D]. 哈尔滨: 哈尔滨工业大学, 2013: 26-34.

[本文引用: 1]

ZHAO Jian-hu, Analysis of stress distribution and life characteristics of flexspline in harmonic drive [D]. Harbin: Harbin Institute of Technology, 2013: 26-34.

[本文引用: 1]

董惠敏, 李德举, 齐书学

基于正交试验和有限元分析的谐波传动柔轮杯体结构优化

[J]. 机械传动, 2013, 37 (1): 34- 38

[本文引用: 1]

DONG Hui-min, LI De-ju, QI Shu-xue

Structural optimization of the flexspline cup in harmonic drives based on orthogonal test and finite element analysis

[J]. Journal of Mechanical Transmission, 2013, 37 (1): 34- 38

[本文引用: 1]

王家序, 周祥祥, 李俊阳, 等

杯形柔轮谐波传动三维双圆弧齿廓设计

[J]. 浙江大学学报: 工学版, 2016, 50 (4): 616- 624

[本文引用: 2]

WANG Jia-xu, ZHOU Xiang-xiang, LI Jun-yang, et al

Three dimensional profile design of cup harmonic drive with double-circular-arc common-tangent tooth profile

[J]. Journal of Zhejiang University: Engineering Science, 2016, 50 (4): 616- 624

[本文引用: 2]

杨勇, 王家序, 周青华, 等

双圆弧谐波齿轮传动柔轮齿廓参数的优化设计

[J]. 四川大学学报: 工程科学版, 2016, 48 (1): 186- 193

[本文引用: 1]

YANG Yong, WANG Jia-xu, ZHOU Qing-hua, et al

Optimization design for flexspline tooth profile parameters of double-circular-arc harmonic drives

[J]. Journal of Sichuan University: Engineering Science Edition, 2016, 48 (1): 186- 193

[本文引用: 1]

李奇. 不同参数对谐波减速器柔轮动态特性的影响[D]. 重庆: 重庆大学, 2016. 26-30.

[本文引用: 1]

LI Qi. Influence of different parameters to the dynamic characteristics of the flexspline of harmonic reducer [D]. Chongqing: Chongqing University, 2016.

[本文引用: 1]

张雷, 杨伟超, 蒋倩倩, 等

谐波减速器核心部件结构应力分析

[J]. 机械传动, 2018, 42 (12): 114- 117

[本文引用: 1]

ZHANG Lei, YANG Wei-chao, JIANG Qian-qian, et al

Structural stress analysis of the core components of harmonic reducer

[J]. Journal of Mechanical Transmission, 2018, 42 (12): 114- 117

[本文引用: 1]

王家序, 刘彪, 周祥祥, 等

双圆弧谐波齿轮传动齿廓设计与参数分析

[J]. 四川大学学报: 工程科学版, 2016, 48 (3): 164- 170

[本文引用: 1]

WANG Jia-xu, LIU Biao, ZHOU Xiang-xiang, et al

Double-circular-arc tooth profile design and parametric analysis of harmonic drive

[J]. Journal of Sichuan University: Engineering Science Edition, 2016, 48 (3): 164- 170

[本文引用: 1]

王家序, 周祥祥, 李俊阳, 等

不同共轭原理的双圆弧齿廓谐波齿轮传动分析

[J]. 四川大学学报: 工程科学版, 2015, 47 (5): 160- 166

[本文引用: 1]

WANG. Jia-xu, ZHOU Xiang-xiang, LI Jun-yang, et al

Double-circular-arc tooth profile of harmonic drive analysis based on different conjugate principle

[J]. Journal of Sichuan University: Engineering Science Edition, 2015, 47 (5): 160- 166

[本文引用: 1]

张志红, 何桢, 郭伟

在响应曲面方法中三类中心复合设计的比较研究

[J]. 沈阳航空航天大学学报, 2007, 24 (1): 87- 91

DOI:10.3969/j.issn.2095-1248.2007.01.026      [本文引用: 1]

ZHANG Zhi-hong, HE Zhen, GUO Wei

A comparative study of three central composite designs in response surface methodology

[J]. Journal of Shenyang Institute of Aeronautical Engineering, 2007, 24 (1): 87- 91

DOI:10.3969/j.issn.2095-1248.2007.01.026      [本文引用: 1]

张立勇, 刘新猛, 王长路, 等

径向变形量对谐波减速器啮合特性及柔轮应力的影响分析

[J]. 机械传动, 2017, 41 (9): 166- 169

[本文引用: 3]

ZHANG Li-yong, LIU Xin-meng, WANG Chang-lu, et al

Influence of radial deformation on stress of flexspline and meshing characteristic of harmonic reducer

[J]. Journal of Mechanical Transmission, 2017, 41 (9): 166- 169

[本文引用: 3]

/