CN112464481B - 机器人用摆线针轮减速机的动态传动精度数值计算方法 - Google Patents
机器人用摆线针轮减速机的动态传动精度数值计算方法 Download PDFInfo
- Publication number
- CN112464481B CN112464481B CN202011392955.8A CN202011392955A CN112464481B CN 112464481 B CN112464481 B CN 112464481B CN 202011392955 A CN202011392955 A CN 202011392955A CN 112464481 B CN112464481 B CN 112464481B
- Authority
- CN
- China
- Prior art keywords
- gear
- dimension
- tooth
- displacement
- cycloidal
- 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
- 239000003638 chemical reducing agent Substances 0.000 title claims abstract description 79
- 230000005540 biological transmission Effects 0.000 title claims abstract description 75
- 238000004364 calculation method Methods 0.000 title claims abstract description 21
- 238000006073 displacement reaction Methods 0.000 claims abstract description 185
- 238000000034 method Methods 0.000 claims abstract description 68
- 238000013178 mathematical model Methods 0.000 claims abstract description 27
- 238000004519 manufacturing process Methods 0.000 claims abstract description 12
- 229920001971 elastomer Polymers 0.000 claims description 23
- 239000000806 elastomer Substances 0.000 claims description 23
- 230000003068 static effect Effects 0.000 claims description 16
- 210000003781 tooth socket Anatomy 0.000 claims description 16
- 238000009826 distribution Methods 0.000 claims description 7
- 238000005452 bending Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000010354 integration Effects 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 3
- 238000013440 design planning Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000005312 nonlinear dynamic Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 239000000969 carrier Substances 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 11
- 241001442234 Cosa Species 0.000 description 4
- 244000089409 Erythrina poeppigiana Species 0.000 description 4
- 235000009776 Rathbunia alamosensis Nutrition 0.000 description 4
- 238000003754 machining Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Retarders (AREA)
Abstract
一种机器人用摆线针轮减速机的动态传动精度数值计算方法,首先对其各个零部件的关键误差进行定义并进行测、算来确定关键误差,然后建立动态传动误差的力学模型,然后根据误差力学模型计算在零件接触处所产生的位移,进一步的确定各零件的微位移在轴或轴承以及轮齿啮合处所产生的位移,再根据上述误差和位移确定各零件在各个零件接触处的作用力;建立其动态传动误差数学模型;通过计算,即可确定其输出轴在任意时刻的实际转角θc、动态传动误差Δθc,该方法综合考虑各种误差、和接触间隙、零件接触变形、旋转零件的惯性载荷等非线性因素,更接近该减速机实况,得到的传动误差更接近实际情况,能够更准确的指导该减速机的零部件生产和装配过程。
Description
技术领域
本发明涉及机器人用摆线针轮减速机的技术领域,特别是涉及一种机器人用摆线针轮减速机的动态传动精度数值计算方法。
背景技术
摆线针轮传动是采用摆线轮和针轮啮合的少齿差行星传动,该传动机构是在二十世纪三十年代初在德国首先研制成功,发明者LorenzBraren于1931年在慕尼黑创建了赛古乐股份有限公司,并开始制造和销售摆线针轮减速机。该传动机构具有传动比大、传动效率高、体积小、重量轻、承载能力大、运行平稳、噪音低、寿命长等特点。
由于摆线针轮减速机零部件制造误差、装配误差以及传动过程中温度变形和弹性变形的存在,输入输出传动误差在所难免。传动误差是指输出轴实际转角与理论转角之间的偏差值,是评价摆线针轮减速机传动精度的重要指标。
传动误差是指:当输入轴单向旋转时,输出轴转角的实际值与理论值的差值。对于行星架为输出轴的摆线针轮减速机,若假设输入轴(即太阳轮)输入端匀速转动且其转角为θs、输出轴(即行星架)实际转角为θc、***传动比为iz,则该***的传动误差为:Δθc=θc-θs/iz。传动精度即是前面提到的传动误差,图1为传动误差的波形示意图。
传动精度的确定对于机器人用摆线针轮减速机的零部件生产和装配过程具有现实的指导意义,能够提前确定零部件的优化配合关系,对采用互补性好的零部件搭配出性能更优的减速机具有重大的意义。
发明内容
本发明的目的是为了解决上述问题,提供一种机器人用摆线针轮减速机的动态传动精度数值计算方法,具体来说是提供一种机器人用三曲柄摆线针轮减速机的动态传动精度数值计算方法。
本发明的技术方案是:该方法的基本原理是:由于其传动***中各零件的加工误差、安装误差、轮齿啮合间隙、轴承间隙、零件接触变形、旋转零件的惯性载荷等非线性因素的存在,而使零件的实际质心位置及其转角偏离了理想位置,且其偏离量用线位移、角位移来描述,并将这些线位移、角位移统称为微位移;同时将零件间的轴承或轴支承、轮齿啮合等接触处用弹簧进行等效,并用轴承或轴刚度、轮齿啮合刚度物理量进行描述。通过分析该传动***中各零件在理想位置时的受力状况,依据D'Alembert原理,建立该***的动力学微分方程组,即动态传动误差的数值计算数学模型;采用Wilson-θ法、Newmark法等非线性直接积分法求解微分方程组;然后进行编程,并利用计算机进行计算,确定该减速机输出轴在任意时刻的实际转角θc,进而确定其动态传动误差Δθc。
具体的该方法包括如下步骤:
S1.根据待仿真计算的减速机的设计参数来确定三曲柄摆线针轮减速机的基本结构参数,
具体包括如下参数:太阳轮齿数Zs、行星轮齿数Zp、摆线轮齿数Zb、针轮齿数Zr、渐开线齿轮模数m(mm)、渐开线齿轮压力角α(°)、太阳轮轮齿变位系数xcs、行星轮轮齿变位系数xcp、行星轮齿宽Bz(mm)、摆线轮轮齿宽度Bc(mm)、曲柄轴偏心距e(mm)、针齿销半径rd(mm)、针齿销或针齿壳齿槽分布圆的理论半径rz(mm)、2个摆线轮、3个行星轮、3个曲柄轴;
S2.定义该减速机中的太阳轮、3个行星轮、2个摆线轮、3个曲柄轴、行星架、针齿销、针齿壳及滚子轴承的关键误差,然后通过高精度测量仪来测量,并根据测量结果运算来确定这些关键误差;
上述各个零件关键误差的定义如下:
1)太阳轮的关键误差的定义
太阳轮的关键误差是指:太阳轮的基圆偏心误差(Es,βs),太阳轮装配误差为(As,γs);其中:Es、As分别表示其基圆偏心误差、装配误差的大小,量纲(um);βs、γs分别表示其基圆偏心误差、装配误差的方向,量纲(°);
2)行星轮的关键误差定义
行星轮的关键误差是指:3个行星轮的基圆偏心误差分别为(Epi,βpi),(i=1,2,3);其中:Epi(i=1,2,3)分别表示行星轮1、2、3的基圆偏心误差大小,量纲(um);βpi(i=1,2,3)分别表示行星轮1、2、3的基圆偏心误差方向,量纲(°);
3)三曲柄摆线轮的关键误差的定义
三曲柄摆线轮的关键误差是指:摆线轮1、2上的曲柄轴孔1、2、3的偏心误差(Ehji,βhji)(j=1,2;i=1,2,3);摆线轮1、2上的轮齿齿槽偏差Rjk(j=1,2)、齿距偏差Pjk(j=1,2);其中:Ehji表示曲柄轴孔的偏心误差大小,量纲(um);βhji表示曲柄轴孔的偏心误差方向,量纲(°);Rjk、Pjk分别表示摆线轮齿的齿槽偏差、齿距偏差大小,量纲(um);
4)针齿壳的关键误差的定义
针齿壳的关键误差是指:Zr个针齿齿槽的平均半径误差δcr(um),针齿壳上Zr个齿槽圆心的齿槽偏差(简称针齿齿槽偏差)Rk(k=1,2,…,Zr)、齿距偏差(简称针轮齿距偏差)Pk(k=1,2,…,Zr);其中:Rk、Pk分别表示第k个针齿的齿槽偏差、齿距偏差大小,量纲(um);
5)曲柄轴的关键误差的定义
曲柄轴的关键误差是指:曲柄轴1、2、3上的偏心凸轮1、2的偏心误差(Ecji,βcji)(j=1,2;i=1,2,3);其中:Ecji表示偏心凸轮的偏心误差大小,量纲(um);βcji表示偏心凸轮的偏心误差方向,量纲(°);
6)行星架的关键误差的定义
行星架的关键误差是指:行星架上曲柄轴孔1、2、3的偏心误差(Ecai,βcai)(i=1,2,3),行星架上安装轴承的轴颈偏心误差(Ac,γc);其中:Ecai、Ac表示曲柄轴孔、轴颈的偏心误差大小,量纲(um);βcai、γc表示曲柄轴孔、轴颈的偏心误差方向,量纲(°);
7)针齿销的关键误差的定义
针齿销的关键误差是指:Zr个针齿销的平均直径误差δzr,量纲(um);
8)轴承的关键误差的定义
在三曲柄摆线针轮减速机中,摆线轮1、2与曲柄轴1、2、3间的轴承采用滚子轴承(图12),其轴承间隙分别为δbji(j=1,2;i=1,2,3),量纲(um);
行星架与曲柄轴1、2、3间的轴承采用圆锥滚子轴承,其轴承间隙分别为δxi(i=1,2,3),量纲(um);
行星架与针齿壳间的轴承采用角接触球轴承,其轴承间隙为δca,量纲(um);
S3.建立计算三曲柄摆线针轮减速机动态传动误差的力学模型,具体如下
1)在力学模型中,将零件间的轴承或轴支承、轮齿啮合等接触处用弹簧进行等效,并用轴承或轴刚度、轮齿啮合刚度物理量进行描述;太阳轮与行星轮i(i=1,2,3)间的轮齿啮合刚度ki(N/um)(i=1,2,3)、曲柄轴i(i=1,2,3)与摆线轮j(j=1,2)轴孔间的轴承支承刚度kji(N/um)(j=1,2;i=1,2,3)、曲柄轴i(i=1,2,3)和行星架轴孔间的轴承支承刚度kbi(N/um)(i=1,2,3)、行星架与针齿壳间轴承支承刚度kca(N/um)、摆线轮与针齿间的啮合刚度kjk(N/um)(k=1,2,…,Zr)、太阳轮轴的支承刚度ks(N/um);以上刚度的确定方法为:ki采用轮齿弯曲变形来确定;ks通过计算太阳轮轴的弯曲变形来求得;而kji、kbi、kca、kjk接触刚度则采用Palmgren公式确定。
Palmgren公式为:
其中,
式中ν1,ν2——弹性体1、2的泊松比;
E1,E2——弹性体1、2的弹性模量(MPa);
δ——两弹性体间的接触变形量(mm);
l——两弹性体间的接触长度(mm);
F——两弹性体间的载荷(N)。
若用k0表示两弹性体间的接触刚度(N/mm),则
F=k0δ (2)
依据式(1)、(2)可推导出两弹性体间的接触刚度k0为:
由式(3)可知,两弹性体间的接触刚度k0主要取决于两弹性体的材料、接触长度及其作用载荷;
2)确定力学模型的静、动坐标系;以太阳轮轴(或行星架、针齿壳)的理论中心O为原点,垂直减速机轴线的断面为平面静坐标系xoy。从输入端即太阳轮端开始,靠近太阳轮的摆线轮编号j=1,另一个摆线轮编号j=2;以摆线轮的理论质心Oj为原点,以摆线轮j(j=1,2)的偏心方向为ηj轴,拟定太阳轮轴逆时针旋转,则将垂直于ηj轴并沿其公转方向转动90°为ξj轴,建立摆线轮j的动坐标系ηjOjξj;在动力学建模过程中,取j=1的摆线轮ηj轴在起始位置时与静坐标系x轴方向一致;
3)在力学模型,将各弹簧设置在各零件坐标轴的正方向接触处或轮齿啮合处,且设定使弹簧受拉为正、受压为负;φi表示摆线轮(或行星架)上的曲柄轴轴孔的相对位置,取φi=2π(i-1)/3(i=1,2,3),ψj表示2个摆线轮的理论质心Oj的相对位置,取ψj=(j-1)π(j=1,2);
S4.依据S2确定的三曲柄摆线针轮减速机中各零件制造误差、间隙,并结合S3所建立的力学模型,确定该传动***在轴或轴承支承处、轮齿啮合处所产生的位移;具体如下:
1)确定太阳轮的基圆偏心误差、装配误差分别在其轮齿啮合处、轴支承处所产生的位移。
若太阳轮的基圆偏心误差为(Es,βs),则其在轮齿啮合处所产生的位移为:
esi=Escos(θs+βs-Ai) (4)
式中Ai=θc+φi+π/2-α′,θc=θs/iz
α′——太阳轮和行星轮的啮合角(°);
θc——行星架的理论转角(°)。
若太阳轮装配误差为(As,γs),则其在太阳轮轴支承处x、y方向上所产生的位移分别为:
esx=Ascosγs (5)
esy=Assinγs (6)
2)确定行星轮的基圆偏心误差在其轮齿啮合处所产生的位移;3个行星轮的基圆偏心误差(Epi,βpi)(i=1,2,3)在轮齿啮合处所产生的位移为:
epi=-Epicos(βpi-θp-Ai)(i=1,2,3) (7)
式中θp——行星轮的理论自转角(°);
3)确定摆线轮上的曲柄轴孔偏心误差在其轴承支承处所产生的位移,及其齿槽偏差、齿距偏差在轮齿啮合处所产生的位移;
定义Oji表示摆线轮j(j=1,2)上曲柄轴孔i(i=1,2,3)的圆心,Rdc表示摆线轮上曲柄轴孔圆心所在圆周的半径。
若摆线轮曲柄轴孔的偏心误差为:(Ehji,βhji)(j=1,2;i=1,2,3),则其在轴承支承处静坐标系x、y方向上所产生的位移分别为:
ehxji=-Ehjicos(θc+φi+βhji) (8)
ehyji=-Ehjisin(θc+φi+βhji) (9)
定义摆线轮的齿槽偏差Rjk、齿距偏差Pjk,O表示针轮中心,Oj表示摆线轮j(j=1,2)中心,P表示摆线轮j与针轮的节点;
若摆线轮j上的齿槽偏差为Rjk(j=1,2)、齿距偏差为Pjk(j=1,2),则其在轮齿啮合处所产生的位移分别为:
式中αjk——第k个针齿中心与第j个摆线轮节点连线与ηj轴正方向夹角(°),其中k=1,2,…,Zr;
φdjk——第k个针齿中心与第j个摆线轮中心Oj连线与ηj轴正方向的夹角(°)
4)确定针齿壳上齿槽圆心的齿槽偏差、齿距偏差在摆线轮齿啮合处所产生的位移;
定义针齿壳上齿槽圆心的齿槽偏差(简称针轮齿槽偏差)Rk、齿距偏差(简称针轮齿槽偏差)Pk,其中:O表示针轮中心,Oj表示摆线轮j(j=1,2)中心,P表示摆线轮j与针轮的节点。
若针轮的齿槽偏差为Rk(k=1,2,…,Zr)、齿距偏差为Pk(k=1,2,…,Zr),则其在轮齿啮合处所产生的位移分别为:
eRk=-Rkcos(αjk-φjk) (12)
式中φjk——第k个针齿半径方向与第j个摆线轮ηj轴正方向间的夹角(°);
5)确定曲柄轴偏心凸轮的偏心误差在其轴承支承处所产生的位移;
定义曲柄轴i的偏心凸轮j偏心误差,Oji表示摆线轮j(j=1,2)上曲柄轴孔i(i=1,2,3)的圆心,Opj表示曲柄轴i(i=1,2,3)上与摆线轮j(j=1,2)相对应的曲柄轴回转中心;
若曲柄轴i上偏心凸轮j的偏心误差为(Ecji,βcji)(j=1,2;i=1,2,3),则其在x、y方向上所产生的位移分别为:
ecxji=Ecjicos(θp+ψj+βcji) (14)
ecyji=-Ecjisin(θp+ψj+βcji) (15)
6)确定行星架上曲柄轴孔、轴颈的偏心误差在轴承支承处所产生的位移;
定义行星架曲柄轴孔i的偏心误差,其中Oci表示行星架上曲柄轴孔i(i=1,2,3)的圆心,Tout表示作用在行星架上的负载转矩;
若行星架曲柄轴孔i的偏心误差为(Ecai,βcai)(i=1,2,3),则它在支承处x、y方向上所产生的位移分别为:
ecaxi=-Ecaicos(θc+φi+βcai) (16)
ecayi=-Ecaisin(θc+φi+βcai) (17)
若假设行星架的轴颈偏心误差为(Ac,γc),则它在支承处x、y方向上所产生的位移分别为:
ecx=Accosγc (18)
ecy=Acsinγc (19)
7)确定针齿销的直径误差、针齿壳上的齿槽半径误差在摆线轮齿啮合处形成的间隙,并确定该间隙在摆线轮齿啮合处所产生的位移;
定义O表示针齿销的理论分布圆心或针齿壳上齿槽理论圆心的分布圆圆心,δjk表示摆线轮齿的啮合间隙;
若针齿销的平均直径误差δzr(um)、针齿齿槽的平均半径误差δcr(um),则其在摆线轮齿啮合处所产生的间隙δjk及其所产生的位移ejk分别为:
δjk=δzr+δcr (20)
ejk=-δjk (21)
8)各轴承间隙在其支承处所产生的位移
若摆线轮j曲柄轴孔i处的滚子轴承间隙为δbji(j=1,2;i=1,2,3),则它在支承处所产生的位移为:
ebji=-δbji (22)
若行星架在曲柄轴孔i处的圆锥滚子轴承间隙为δxi(i=1,2,3),则它在支承处所产生的位移为:
exi=-δxi (23)
若行星架与针齿壳间的角接触球轴承间隙为δca,则它在支承处所产生的位移为:
eca=-δca (24)
S5.确定三曲柄摆线针轮减速机中各零件的微位移在轴或轴承以及轮齿啮合处所产生的位移;具体如下:
在分析三曲柄摆线针轮减速机中各零件的微位移在接触处所产生的位移时,假设针齿壳固定不动,行星轮与曲柄轴固结在一起,且太阳轮、带曲柄轴的行星轮、摆线轮及行星架的线位移方向与静坐标系或动坐标系的坐标轴正向相同、角位移方向与设计拟定方向相同;
1)太阳轮的微位移在轮齿啮合处、太阳轮轴支承处所产生的位移;
若太阳轮的微位移为xs、ys、θsa-θs,则其在轮齿啮合处所产生的位移为:
si=xscosAi+yssinAi+Rbs(θsa-θs) (25)
式中Rbs——表示太阳轮的基圆半径(mm);
在轴支承处x、y方向上所产生的位移为:
sx=xs (26)
sy=ys (27)
2)行星轮及曲柄轴的微位移在其轮齿啮合处、摆线轮滚子轴承支承处以及行星架圆锥滚子轴承支承处所产生的位移;
若令行星轮i(i=1,2,3)的微位移为xpi、ypi和θpi-θp,则其在轮齿啮合处所产生的位移为:
spi=-xpicosAi-ypisinAi-Rbp(θpi-θp) (28)
式中Rbp——行星轮的基圆半径(mm);
θpi——行星轮的实际自转角(°)。
在摆线轮滚子轴承支承处x、y方向上所产生的位移为:
sjix=xpi-e(θpi-θp)sin(θp+ψj) (29)
sjiy=ypi-e(θpi-θp)cos(θp+ψj) (30)
式中e——曲柄轴的偏心距(mm)。
在行星架圆锥滚子轴承支承处x、y方向上所产生的位移为:
six=xpi (31)
siy=ypi (32)
3)摆线轮的微位移在其轮齿啮合处、曲柄轴滚子轴承支承处所产生的位移;
设摆线轮j(j=1,2)的微位移分别为ηdj、θdj-θc和θOj-θp,其中ηdj表示摆线轮的线微位移,θdj-θc表示摆线轮的自转角微位移,θOj-θp表示摆线轮的公转角微位移;
则摆线轮j(j=1,2)的微位移在其轮齿啮合处所产生的位移为:
sjk=ηdjcosαjk-Rd(θdj-θc)sinαjk+e(θOj-θp)sinαjk (33)
式中Rd——摆线轮的节圆半径(mm);
在曲柄轴滚子轴承支承处x、y方向上所产生的位移分别为:
scjix=Rdc(θdj-θc)sin(θc+φi)+e(θOj-θp)sin(θp+ψj)-ηdjcos(θp+ψj) (34)
scjiy=-Rdc(θdj-θc)cos(θc+φi)+ηdjsin(θp+ψj)+e(θOj-θp)cos(θp+ψj) (35)
式中Rdc——太阳轮与行星轮间的中心距(mm);
4)行星架的微位移在曲柄轴圆锥滚子轴承支承处、针齿壳角接触球轴承支承处所产生的位移;
设行星架的三个微位移为xca、yca和θca-θc,则其在曲柄轴圆锥滚子轴承支承处x、y方向上所产生的位移分别为:
sbix=-xca+Rdc(θca-θc)sin(θc+φi) (36)
sbiy=-yca-Rdc(θca-θc)cos(θc+φi) (37)
在针齿壳角接触球轴承支承处x、y方向上所产生的位移(参考图23)分别为:
scax=xca (38)
scay=yca (39)
S6.确定三曲柄摆线针轮减速机中,各零件在轴或轴承支承处、轮齿啮合处的作用力;
依据上述步骤S4、S5,确定的各零件制造误差、间隙及微位移在其支承、轮齿啮合处所产生的位移,同时考虑在理想状态时,各零件所传递的理论载荷,从而获得各零件在轴或轴承、轮齿啮合处的作用力。若作用力为正,则表示两物体间的接触力受压,反之受拉
1)太阳轮在其支承处x、y方向上的作用力分别为Fsx(N)、Fsy(N),如公式(40)、(41)所示;
Fsx=ks(sx-esx)=ks(xs-esx) (40)
Fsy=ks(sy-esy)=ks(ys-esy) (41)
2)太阳轮与行星轮i(i=1,2,3)间的轮齿啮合力为Fi(N),如公式(42)所示;
式中Fi0——表示太阳轮与行星轮间的理论轮齿啮合力(N);
若公式(42)中[]内数值≤0,则表示作用力为拉力,说明太阳轮与行星轮间有间隙存在而没有接触,此时取Fi的值为零;
3)摆线轮j(j=1,2)与曲柄轴i(i=1,2,3)接触处x、y方向上的作用力分别为Fjix(N)、Fjiy(N),如公式(43)、(44)所示;
式中Fjix0——表示理想状态时,摆线轮j(j=1,2)与曲柄轴i(i=1,2,3)
间的滚子轴承在x方向上所承受的载荷(N);
式中Fjiy0——表示理想状态时,摆线轮j(j=1,2)与曲柄轴i(i=1,2,3)
间的滚子轴承在y方向上所承受的载荷(N);
公式(43)、(44)中±号是依据摆线轮与曲柄轴接触处的变形量大小进行确定;
4)摆线轮j(j=1,2)轮齿与第k(k=1,2,…,Zr)个针齿间的作用力为Fjk(N),如公式(45)所示。
式中Fjk0——表示理想状态时摆线轮j(j=1,2)的轮齿啮合力(N);
若Fjk≤0,则表示轮齿啮合力为拉力,说明摆线轮轮齿与针齿之间有间隙存在而没有接触,则取Fijk=0;
5)行星架与曲柄轴i(i=1,2,3)接触处x、y方向上的作用力分别为Fcix(N)、Fciy(N),如公式(46)、(47)所示;
式中Fci0——表示理想状态时行星架与曲柄轴i(i=1,2,3)间的作用力(N)。
公式(46)、(47)中±号是依据行星架与曲柄轴接触处变形量的大小进行确定。
6)行星架与针齿壳支承处x、y方向上的作用力分别为Fcax(N)、Fcay(N),如公式(48)、(49)所示。
Fcax=kca(scax-ecx±eca)=kca(xca-ecx±eca) (48)
Fcay=kca(scay-ecy±eca)=kca(yca-ecy±eca) (49)
公式(48)、(49)中±号是依据行星架与针齿壳支承处变形量的大小来确定。
S7.建立计算三曲柄摆线针轮减速机动态传动误差的数学模型;
通过分析三曲柄摆线针轮减速机中各零件任意位置时的受力状况,依据D'Alembert原理,建立该***动态传动误差的数学模型;在建立数学模型时,考虑了摆线轮、行星轮及曲柄轴的科氏(Coriolis)加速度,忽略了传动***中的摩擦阻力。通过整理,得出三曲柄摆线针轮减速机动态传动误差的数学模型如下:
式中ms——太阳轮的质量(kg);
msp——行星轮及曲柄轴的质量和(kg);
mbx——摆线轮的质量(kg);
Jop——行星轮及曲柄轴的转动惯量(kg·m2);
Joj——摆线轮的转动惯量(kg·m2);
ωc——行星架理论角速度(rad/s);
ωp——行星轮自转理论角速度(rad/s)。
上述数学模型可整理成以下矩阵形式:
式中M、C、K——分别为质量阵、阻尼阵及刚度阵,矩阵阶数均为21×21;
X——位移向量,其中X=(xs,ys,θsa,xp1,yp1,θp1,xp2,yp2,θp2,xp3,yp3,
θp3,ηd1,θo1,θd1,ηd2,θo2,θd2,xca,yca,θca)T;
Q——广义力矢量。
由于该传动***中的刚度系数kji、kbi、kca、kjk是位移的函数,因此所建立的三曲柄摆线针轮减速机动态传动误差的数学模型属于非线性动力学模型;
S8.采用非线性Newmark法、Wilson-θ法的直接积分法,求解所建立的三曲柄摆线针轮减速机动态传动误差的数学模型;通过编程,并利用计算机进行数值计算,即可确定该减速机输出轴在任意时刻的实际转角θc、动态传动误差Δθc。
优选的,步骤S2中采用的误差测量仪器具有智能模式识别、自动扫描等功能,并具有测量精度高、灵敏度高、效率高的特点;其主要技术参数的要求:最小读取单位≤0.1um、测量精度不低于1.2+4L/1000(um)(L为被测工件长度,量纲为mm)、工作台平面尺寸不小于450mm×400mm。
优选的,步骤S2中采用高精度三坐标测量仪或其他零件误差测量仪器;计算方法采用最小二乘法理论、数学变换与滤波理论以及优化方法的基本理论与方法,测量数据的通过计算机处理技术或方法,从而确定上述各零件的关键误差。
本发明的有益效果是:该方法首先对减速机的各个零部件的关键误差进行定义并进行测量和计算来确定关键误差,然后建立动态传动误差的力学模型,然后根据误差力学模型计算在轴或轴承支承处、轮齿啮合处所产生的位移,进一步的确定各零件的微位移在轴或轴承以及轮齿啮合处所产生的位移,再根据上述误差和位移确定各零件在各个零件接触处的作用力;建立该减速机的动态传动误差数学模型;通过计算机数值计算,即可确定该减速机输出轴在任意时刻的实际转角θc、动态传动误差Δθc,该方法综合考虑各种误差、和接触间隙、零件接触变形、旋转零件的惯性载荷等非线性因素,更接近该减速机实况,得到的传动误差更接近实际情况,能够更准确的指导该减速机的零部件生产和装配过程。
附图说明
图1是机器人用三曲柄摆线针轮减速机的传动误差的波形示意图;
图2是机器人用三曲柄摆线针轮减速机传动***示意图;
图3是机器人用三曲柄摆线针轮减速机的零件拆解图;
图4是其中的太阳轮轴示意图;
图5是其中的行星齿轮示意图;
图6是三曲柄摆线轮的示意图;
图7是行星架与针齿壳处的角接触球轴承的示意图;
图8是其中针齿壳的示意图;
图9是其中曲柄轴的示意图;
图10-11是其中行星架的示意图;
图12是针齿销的示意图;
图13是摆线轮与曲柄轴处的滚子轴承的示意图;
图14是曲柄轴与行星架处的圆锥滚子轴承的示意图;
图15是该三曲柄摆线针轮减速机动态传动误差的力学模型示意图;
图16是太阳轮和第i个行星轮的基圆偏心误差的示意图;
图17是太阳轮的装配误差的示意图;
图18是摆线轮曲柄轴孔的偏心误差的示意图;
图19是摆线轮的齿槽偏差、齿距偏差的示意图;
图20是针齿壳上齿槽圆心的齿槽偏差、齿距偏差的示意图;
图21是曲柄轴上偏心凸轮的偏心误差的示意图;
图22是行星架曲柄轴孔的偏心误差的示意图;
图23是行星架的轴颈偏心误差的示意图;
图24是针齿销与齿槽间的摆线轮齿啮合间隙的示意图;
具体实施方式
实施例一:参见图1-24,图中一种机器人用摆线针轮减速机的动态传动精度数值计算方法,包括如下步骤:
S1.根据待仿真计算的减速机的设计参数来确定三曲柄摆线针轮减速机的基本结构参数,
具体包括如下参数:太阳轮齿数Zs、行星轮齿数Zp、摆线轮齿数Zb、针轮齿数Zr、渐开线齿轮模数m(mm)、渐开线齿轮压力角α(°)、太阳轮轮齿变位系数xcs、行星轮轮齿变位系数xcp、行星轮齿宽Bz(mm)、摆线轮轮齿宽度Bc(mm)、曲柄轴偏心距e(mm)、针齿销半径rd(mm)、针齿销或针齿壳齿槽分布圆的理论半径rz(mm)、2个摆线轮、3个行星轮、3个曲柄轴;如图2~图3所示。
S2.定义该减速机中的太阳轮、3个行星轮、2个摆线轮、3个曲柄轴、行星架、针齿销、针齿壳及滚子轴承的关键误差,然后通过高精度测量仪来测量,并根据测量结果运算来确定这些关键误差;
上述各个零件关键误差的定义如下:
1)太阳轮的关键误差的定义
参见图4,太阳轮的关键误差是指:太阳轮的基圆偏心误差(Es,βs),太阳轮装配误差为(As,γs);其中:Es、As分别表示其基圆偏心误差、装配误差的大小,量纲(um);βs、γs分别表示其基圆偏心误差、装配误差的方向,量纲(°);
2)行星轮的关键误差定义
参见图5,行星轮的关键误差是指:3个行星轮的基圆偏心误差分别为(Epi,βpi),(i=1,2,3);其中:Epi(i=1,2,3)分别表示行星轮1、2、3的基圆偏心误差大小,量纲(um);βpi(i=1,2,3)分别表示行星轮1、2、3的基圆偏心误差方向,量纲(°);
3)三曲柄摆线轮的关键误差的定义
参见图6,三曲柄摆线轮的关键误差是指:摆线轮1、2上的曲柄轴孔1、2、3的偏心误差(Ehji,βhji)(j=1,2;i=1,2,3);摆线轮1、2上的轮齿齿槽偏差Rjk(j=1,2)、齿距偏差Pjk(j=1,2);其中:Ehji表示曲柄轴孔的偏心误差大小,量纲(um);βhji表示曲柄轴孔的偏心误差方向,量纲(°);Rjk、Pjk分别表示摆线轮齿的齿槽偏差、齿距偏差大小,量纲(um);
4)针齿壳的关键误差的定义
参见图8,针齿壳的关键误差是指:Zr个针齿齿槽的平均半径误差δcr(um),针齿壳上Zr个齿槽圆心的齿槽偏差(简称针齿齿槽偏差)Rk(k=1,2,…,Zr)、齿距偏差(简称针轮齿距偏差)Pk(k=1,2,…,Zr);其中:Rk、Pk分别表示第k个针齿的齿槽偏差、齿距偏差大小,量纲(um);
5)曲柄轴的关键误差的定义
参见图9,曲柄轴的关键误差是指:曲柄轴1、2、3上的偏心凸轮1、2的偏心误差(Ecji,βcji)(j=1,2;i=1,2,3);其中:Ecji表示偏心凸轮的偏心误差大小,量纲(um);βcji表示偏心凸轮的偏心误差方向,量纲(°);
6)行星架的关键误差的定义
参见图10-11,行星架的关键误差是指:行星架上曲柄轴孔1、2、3的偏心误差(Ecai,βcai)(i=1,2,3),行星架上安装轴承的轴颈偏心误差(Ac,γc);其中:Ecai、Ac表示曲柄轴孔、轴颈的偏心误差大小,量纲(um);βcai、γc表示曲柄轴孔、轴颈的偏心误差方向,量纲(°);
7)针齿销的关键误差的定义
参见图12,针齿销的关键误差是指:Zr个针齿销的平均直径误差δzr,量纲(um);
8)轴承的关键误差的定义
在三曲柄摆线针轮减速机中,摆线轮1、2与曲柄轴1、2、3间的轴承采用滚子轴承,参见图13,其轴承间隙分别为δbji(j=1,2;i=1,2,3),量纲(um);
行星架与曲柄轴1、2、3间的轴承采用圆锥滚子轴承,参见图14,,其轴承间隙分别为δxi(i=1,2,3),量纲(um);
行星架与针齿壳间的轴承采用角接触球轴承,参见图7,其轴承间隙为δca,量纲(um);
S3.建立计算三曲柄摆线针轮减速机动态传动误差的力学模型,具体如下
1)参见图15,在力学模型中,将零件间的轴承或轴支承、轮齿啮合等接触处用弹簧进行等效,并用轴承或轴刚度、轮齿啮合刚度物理量进行描述;太阳轮与行星轮i(i=1,2,3)间的轮齿啮合刚度ki(N/um)(i=1,2,3)、曲柄轴i(i=1,2,3)与摆线轮j(j=1,2)轴孔间的轴承支承刚度kji(N/um)(j=1,2;i=1,2,3)、曲柄轴i(i=1,2,3)和行星架轴孔间的轴承支承刚度kbi(N/um)(i=1,2,3)、行星架与针齿壳间轴承支承刚度kca(N/um)、摆线轮与针齿间的啮合刚度kjk(N/um)(k=1,2,…,Zr)、太阳轮轴的支承刚度ks(N/um);以上刚度的确定方法为:ki采用轮齿弯曲变形来确定;ks通过计算太阳轮轴的弯曲变形来求得;而kji、kbi、kca、kjk接触刚度则采用Palmgren公式确定。
Palmgren公式为:
其中,
式中ν1,ν2——弹性体1、2的泊松比;
E1,E2——弹性体1、2的弹性模量(MPa);
δ——两弹性体间的接触变形量(mm);
l——两弹性体间的接触长度(mm);
F——两弹性体间的载荷(N)。
若用k0表示两弹性体间的接触刚度(N/mm),则
F=k0δ (2)
依据式(1)、(2)可推导出两弹性体间的接触刚度k0为:
由式(3)可知,两弹性体间的接触刚度k0主要取决于两弹性体的材料、接触长度及其作用载荷;
2)确定力学模型的静、动坐标系;以太阳轮轴(或行星架、针齿壳)的理论中心O为原点,垂直减速机轴线的断面为平面静坐标系xoy。从输入端即太阳轮端开始,靠近太阳轮的摆线轮编号j=1,另一个摆线轮编号j=2;以摆线轮的理论质心Oj为原点,以摆线轮j(j=1,2)的偏心方向为ηj轴,拟定太阳轮轴逆时针旋转,则将垂直于ηj轴并沿其公转方向转动90°为ξj轴,建立摆线轮j的动坐标系ηjOjξj;在动力学建模过程中,取j=1的摆线轮ηj轴在起始位置时与静坐标系x轴方向一致;
3)在力学模型,将各弹簧设置在各零件坐标轴的正方向接触处或轮齿啮合处,且设定使弹簧受拉为正、受压为负;φi表示摆线轮(或行星架)上的曲柄轴轴孔的相对位置,取φi=2π(i-1)/3(i=1,2,3),ψj表示2个摆线轮的理论质心Oj的相对位置,取ψj=(j-1)π(j=1,2);
S4.依据S2确定的三曲柄摆线针轮减速机中各零件制造误差、间隙,并结合S3所建立的力学模型,图15所示,确定该传动***在轴或轴承支承处、轮齿啮合处所产生的位移;具体如下:
1)确定太阳轮的基圆偏心误差、装配误差分别在其轮齿啮合处、轴支承处所产生的位移。图16是太阳轮和第i个行星轮的基圆偏心误差示意图,图17是太阳轮的装配误差示意图。
若太阳轮的基圆偏心误差为(Es,βs),则其在轮齿啮合处所产生的位移为:
esi=Escos(θs+βs-Ai) (4)
式中Ai=θc+φi+π/2-α′,θc=θs/iz
α′——太阳轮和行星轮的啮合角(°);
θc——行星架的理论转角(°)。
若太阳轮装配误差为(As,γs),则其在太阳轮轴支承处x、y方向上所产生的位移分别为:
esx=Ascosγs (5)
esy=Assinγs (6)
2)确定行星轮的基圆偏心误差在其轮齿啮合处所产生的位移,如图16所示;3个行星轮的基圆偏心误差(Epi,βpi)(i=1,2,3)在轮齿啮合处所产生的位移为:
epi=-Epicos(βpi-θp-Ai)(i=1,2,3) (7)
式中θp——行星轮的理论自转角(°);
3)确定摆线轮上的曲柄轴孔偏心误差在其轴承支承处所产生的位移,及其齿槽偏差、齿距偏差在轮齿啮合处所产生的位移;图18是摆线轮曲柄轴孔的偏心误差示意图;
定义其中Oji表示摆线轮j(j=1,2)上曲柄轴孔i(i=1,2,3)的圆心,Rdc表示摆线轮上曲柄轴孔圆心所在圆周的半径。
若摆线轮曲柄轴孔的偏心误差为:(Ehji,βhji)(j=1,2;i=1,2,3),则其在轴承支承处静坐标系x、y方向上所产生的位移分别为:
ehxji=-Ehjicos(θc+φi+βhji) (8)
ehyji=-Ehjisin(θc+φi+βhji) (9)
图19是摆线轮的齿槽偏差Rjk、齿距偏差Pjk示意图,O表示针轮中心,Oj表示摆线轮j(j=1,2)中心,P表示摆线轮j与针轮的节点;
若摆线轮j上的齿槽偏差为Rjk(j=1,2)、齿距偏差为Pjk(j=1,2),则其在轮齿啮合处所产生的位移分别为:
式中αjk——第k个针齿中心与第j个摆线轮节点连线与ηj轴正方向夹角(°),其中k=1,2,…,Zr;
φdjk——第k个针齿中心与第j个摆线轮中心Oj连线与ηj轴正方向的夹角(°)
4)确定针齿壳上齿槽圆心的齿槽偏差、齿距偏差在摆线轮齿啮合处所产生的位移;
图20是针齿壳上齿槽圆心的齿槽偏差(简称针轮齿槽偏差)Rk、齿距偏差(简称针轮齿槽偏差)Pk示意图,其中:O表示针轮中心,Oj表示摆线轮j(j=1,2)中心,P表示摆线轮j与针轮的节点。
若针轮的齿槽偏差为Rk(k=1,2,…,Zr)、齿距偏差为Pk(k=1,2,…,Zr),则其在轮齿啮合处所产生的位移分别为:
eRk=-Rkcos(αjk-φjk) (12)
式中φjk——第k个针齿半径方向与第j个摆线轮ηj轴正方向间的夹角(°);
5)确定曲柄轴偏心凸轮的偏心误差在其轴承支承处所产生的位移;
图21是曲柄轴i的偏心凸轮j偏心误差示意图,Oji表示摆线轮j(j=1,2)上曲柄轴孔i(i=1,2,3)的圆心,Opj表示曲柄轴i(i=1,2,3)上与摆线轮j(j=1,2)相对应的曲柄轴回转中心;
若曲柄轴i上偏心凸轮j的偏心误差为(Ecji,βcji)(j=1,2;i=1,2,3),则其在x、y方向上所产生的位移分别为:
ecxji=Ecjicos(θp+ψj+βcji) (14)
ecyji=-Ecjisin(θp+ψj+βcji) (15)
6)确定行星架上曲柄轴孔、轴颈的偏心误差在轴承支承处所产生的位移;
图22是行星架曲柄轴孔i的偏心误差示意图,其中Oci表示行星架上曲柄轴孔i(i=1,2,3)的圆心,Tout表示作用在行星架上的负载转矩;
若行星架曲柄轴孔i的偏心误差为(Ecai,βcai)(i=1,2,3),则它在支承处x、y方向上所产生的位移分别为:
ecaxi=-Ecaicos(θc+φi+βcai) (16)
ecayi=-Ecaisin(θc+φi+βcai) (17)
若假设行星架的轴颈偏心误差为(Ac,γc),则它在支承处x、y方向上所产生的位移分别为:
ecx=Accosγc (18)
ecy=Acsinγc (19)
7)确定针齿销的直径误差、针齿壳上的齿槽半径误差在摆线轮齿啮合处形成的间隙,并确定该间隙在摆线轮齿啮合处所产生的位移;
图24是针齿销的直径误差、针齿壳上的齿槽半径误差在摆线轮齿啮合处形成的间隙示意图,若针齿销的平均直径误差δzr(um)、针齿齿槽的平均半径误差δcr(um),则其在摆线轮齿啮合处所产生的间隙δjk及其所产生的位移ejk分别为:
δjk=δzr+δcr (20)
ejk=-δjk (21)
8)各轴承间隙在其支承处所产生的位移
若摆线轮j曲柄轴孔i处的滚子轴承间隙为δbji(j=1,2;i=1,2,3),则它在支承处所产生的位移为:
ebji=-δbji (22)
若行星架在曲柄轴孔i处的圆锥滚子轴承间隙为δxi(i=1,2,3),则它在支承处所产生的位移为:
exi=-δxi (23)
若行星架与针齿壳间的角接触球轴承间隙为δca,则它在支承处所产生的位移为:
eca=-δca (24)
S5.确定三曲柄摆线针轮减速机中各零件的微位移在轴或轴承以及轮齿啮合处所产生的位移;具体如下:
在分析三曲柄摆线针轮减速机中各零件的微位移在接触处所产生的位移时,假设针齿壳固定不动,行星轮与曲柄轴固结在一起,且太阳轮、带曲柄轴的行星轮、摆线轮及行星架的线位移方向与静坐标系或动坐标系的坐标轴正向相同、角位移方向与设计拟定方向相同,即图15所示方向相同;
1)太阳轮的微位移在轮齿啮合处、太阳轮轴支承处所产生的位移;
参考图16,若太阳轮的微位移为xs、ys、θsa-θs,则其在轮齿啮合处所产生的位移为:
si=xscosAi+yssinAi+Rbs(θsa-θs) (25)
式中Rbs——表示太阳轮的基圆半径(mm);
在轴支承处x、y方向上所产生的位移为:
sx=xs (26)
sy=ys (27)
2)行星轮及曲柄轴的微位移在其轮齿啮合处、摆线轮滚子轴承支承处以及行星架圆锥滚子轴承支承处所产生的位移;
参考图16,若令行星轮i(i=1,2,3)的微位移为xpi、ypi和θpi-θp,则其在轮齿啮合处所产生的位移为:
spi=-xpicosAi-ypisinAi-Rbp(θpi-θp) (28)
式中Rbp——行星轮的基圆半径(mm);
θpi——行星轮的实际自转角(°)。
在摆线轮滚子轴承支承处x、y方向上所产生的位移为:
sjix=xpi-e(θpi-θp)sin(θp+ψj) (29)
sjiy=ypi-e(θpi-θp)cos(θp+ψj) (30)
式中e——曲柄轴的偏心距(mm)。
在行星架圆锥滚子轴承支承处x、y方向上所产生的位移为:
six=xpi (31)
siy=ypi (32)
3)摆线轮的微位移在其轮齿啮合处、曲柄轴滚子轴承支承处所产生的位移;
参考图15设摆线轮j(j=1,2)的微位移分别为ηdj、θdj-θc和θOj-θp,其中ηdj表示摆线轮的线微位移,θdj-θc表示摆线轮的自转角微位移,θOj-θp表示摆线轮的公转角微位移;
则摆线轮j(j=1,2)的微位移在其轮齿啮合处所产生的位移为:
sjk=ηdjcosαjk-Rd(θdj-θc)sinαjk+e(θOj-θp)sinαjk (33)
式中Rd——摆线轮的节圆半径(mm);
在曲柄轴滚子轴承支承处x、y方向上所产生的位移分别为:
scjix=Rdc(θdj-θc)sin(θc+φi)+e(θOj-θp)sin(θp+ψj)-ηdjcos(θp+ψj) (34)
scjiy=-Rdc(θdj-θc)cos(θc+φi)+ηdjsin(θp+ψj)+e(θOj-θp)cos(θp+ψj) (35)
式中Rdc——太阳轮与行星轮间的中心距(mm);
4)行星架的微位移在曲柄轴圆锥滚子轴承支承处、针齿壳角接触球轴承支承处所产生的位移;
参考图22,设行星架的三个微位移为xca、yca和θca-θc,则其在曲柄轴圆锥滚子轴承支承处x、y方向上所产生的位移分别为:
sbix=-xca+Rdc(θca-θc)sin(θc+φi) (36)
sbiy=-yca-Rdc(θca-θc)cos(θc+φi) (37)
在针齿壳角接触球轴承支承处x、y方向上所产生的位移(参考图23)分别为:
scax=xca (38)
scay=yca (39)
S6.确定三曲柄摆线针轮减速机中,各零件在轴或轴承支承处、轮齿啮合处的作用力;
依据上述步骤S4、S5,确定的各零件制造误差、间隙及微位移在其支承、轮齿啮合处所产生的位移,同时考虑在理想状态时,各零件所传递的理论载荷,从而获得各零件在轴或轴承、轮齿啮合处的作用力。若作用力为正,则表示两物体间的接触力受压,反之受拉
1)太阳轮在其支承处x、y方向上的作用力分别为Fsx(N)、Fsy(N),如公式(40)、(41)所示;
Fsx=ks(sx-esx)=ks(xs-esx) (40)
Fsy=ks(sy-esy)=ks(ys-esy) (41)
2)太阳轮与行星轮i(i=1,2,3)间的轮齿啮合力为Fi(N),如公式(42)所示;
式中Fi0——表示太阳轮与行星轮间的理论轮齿啮合力(N);
若公式(42)中[]内数值≤0,则表示作用力为拉力,说明太阳轮与行星轮间有间隙存在而没有接触,此时取Fi的值为零;
3)摆线轮j(j=1,2)与曲柄轴i(i=1,2,3)接触处x、y方向上的作用力分别为Fjix(N)、Fjiy(N),如公式(43)、(44)所示;
式中Fjix0——表示理想状态时,摆线轮j(j=1,2)与曲柄轴i(i=1,2,3)
间的滚子轴承在x方向上所承受的载荷(N);
式中Fjiy0——表示理想状态时,摆线轮j(j=1,2)与曲柄轴i(i=1,2,3)
间的滚子轴承在y方向上所承受的载荷(N);
公式(43)、(44)中±号是依据摆线轮与曲柄轴接触处的变形量大小进行确定;
4)摆线轮j(j=1,2)轮齿与第k(k=1,2,…,Zr)个针齿间的作用力为Fjk(N),如公式(45)所示。
式中Fjk0——表示理想状态时摆线轮j(j=1,2)的轮齿啮合力(N);
若Fjk≤0,则表示轮齿啮合力为拉力,说明摆线轮轮齿与针齿之间有间隙存在而没有接触,则取Fijk=0;
5)行星架与曲柄轴i(i=1,2,3)接触处x、y方向上的作用力分别为Fcix(N)、Fciy(N),如公式(46)、(47)所示;
式中Fci0——表示理想状态时行星架与曲柄轴i(i=1,2,3)间的作用力(N)。
/>
公式(46)、(47)中±号是依据行星架与曲柄轴接触处变形量的大小进行确定。
6)行星架与针齿壳支承处x、y方向上的作用力分别为Fcax(N)、Fcay(N),如公式(48)、(49)所示。
Fcax=kca(scax-ecx±eca)=kca(xca-ecx±eca) (48)
Fcay=kca(scay-ecy±eca)=kca(yca-ecy±eca) (49)
公式(48)、(49)中±号是依据行星架与针齿壳支承处变形量的大小来确定。
S7.建立计算三曲柄摆线针轮减速机动态传动误差的数学模型;
通过分析三曲柄摆线针轮减速机中各零件任意位置时的受力状况,依据D'Alembert原理,建立该***动态传动误差的数学模型;在建立数学模型时,考虑了摆线轮、行星轮及曲柄轴的科氏(Coriolis)加速度,忽略了传动***中的摩擦阻力。通过整理,得出三曲柄摆线针轮减速机动态传动误差的数学模型如下:
式中ms——太阳轮的质量(kg);
msp——行星轮及曲柄轴的质量和(kg);
mbx——摆线轮的质量(kg);
Jop——行星轮及曲柄轴的转动惯量(kg·m2);
Joj——摆线轮的转动惯量(kg·m2);
ωc——行星架理论角速度(rad/s);
ωp——行星轮自转理论角速度(rad/s)。
上述数学模型可整理成以下矩阵形式:
式中M、C、K——分别为质量阵、阻尼阵及刚度阵,矩阵阶数均为21×21;
X——位移向量,其中X=(xs,ys,θsa,xp1,yp1,θp1,xp2,yp2,θp2,xp3,yp3,
θp3,ηd1,θo1,θd1,ηd2,θo2,θd2,xca,yca,θca)T;
Q——广义力矢量。
由于该传动***中的刚度系数kji、kbi、kca、kjk是位移的函数,因此所建立的三曲柄摆线针轮减速机动态传动误差的数学模型属于非线性动力学模型;
S8.采用非线性Newmark法、Wilson-θ法的直接积分法,求解所建立的三曲柄摆线针轮减速机动态传动误差的数学模型;通过编程,并利用计算机进行数值计算,即可确定该减速机输出轴在任意时刻的实际转角θc、动态传动误差Δθc,如图1所示。
步骤S2中采用的误差测量仪器具有智能模式识别、自动扫描等功能,并具有测量精度高、灵敏度高、效率高的特点;其主要技术参数的要求:最小读取单位≤0.1um、测量精度不低于1.2+4L/1000(um)(L为被测工件长度,量纲为mm)、工作台平面尺寸不小于450mm×400mm。
步骤S2中采用高精度三坐标测量仪或其他零件误差测量仪器;计算方法采用最小二乘法理论、数学变换与滤波理论以及优化方法的基本理论与方法,测量数据的通过计算机处理技术或方法,从而确定上述各零件的关键误差。
该方法综合考虑零件加工误差、安装误差、轮齿啮合间隙、轴承间隙、零件接触变形、旋转零件的惯性载荷等非线性因素,提出来一种动态传动误差数值计算方法。不但能够对传动误差即传动精度进行仿真测算,特别是能用于仿真计算高速运转或带载荷运转的这类减速机,因为这类减速机的零部件自身的惯性轮,以及输入和输出之间的扭矩产生的力,都会对零部件的运动接触间隙和变形量产生影响,相比单纯的静态仿真来说,该方法结合了工程实际来进行仿真,因此更适用于指导实际生产。
Claims (3)
1.一种机器人用摆线针轮减速机的动态传动精度数值计算方法,包括如下步骤:
S1.根据待仿真计算的减速机的设计参数来确定三曲柄摆线针轮减速机的基本结构参数,
具体包括如下参数:太阳轮齿数Zs、行星轮齿数Zp、摆线轮齿数Zb、针轮齿数Zr、渐开线齿轮模数m,量纲mm、渐开线齿轮压力角α,量纲°、太阳轮轮齿变位系数xcs、行星轮轮齿变位系数xcp、行星轮齿宽Bz,量纲mm、摆线轮轮齿宽度Bc,量纲mm、曲柄轴偏心距e,量纲mm、针齿销半径rd,量纲mm、针齿销或针齿壳齿槽分布圆的理论半径rz,量纲mm、2个摆线轮、3个行星轮、3个曲柄轴;
S2.定义该减速机中的太阳轮、3个行星轮、2个摆线轮、3个曲柄轴、行星架、针齿销、针齿壳及滚子轴承的关键误差,然后通过高精度测量仪来测量,并根据测量结果运算来确定这些关键误差;各个零件关键误差的定义如下:
1)太阳轮的关键误差的定义
太阳轮的关键误差是指:太阳轮的基圆偏心误差(Es,βs),太阳轮装配误差为(As,γs);其中:Es、As分别表示其基圆偏心误差、装配误差的大小,量纲um;βs、γs分别表示其基圆偏心误差、装配误差的方向,量纲°;
2)行星轮的关键误差定义
行星轮的关键误差是指:3个行星轮的基圆偏心误差分别为(Epi,βpi),i=1,2,3;其中:Epi,i=1,2,3分别表示行星轮1、2、3的基圆偏心误差大小,量纲um;βpi,i=1,2,3分别表示行星轮1、2、3的基圆偏心误差方向,量纲°;
3)三曲柄摆线轮的关键误差的定义
三曲柄摆线轮的关键误差是指:摆线轮1、2上的曲柄轴孔1、2、3的偏心误差(Ehji,βhji),j=1,2;i=1,2,3;摆线轮1、2上的轮齿齿槽偏差Rjk,j=1,2、齿距偏差Pjk,j=1,2;其中:Ehji表示曲柄轴孔的偏心误差大小,量纲um;βhji表示曲柄轴孔的偏心误差方向,量纲°;Rjk、Pjk分别表示摆线轮齿的齿槽偏差、齿距偏差大小,量纲um;
4)针齿壳的关键误差的定义
针齿壳的关键误差是指:Zr个针齿齿槽的平均半径误差δcr,量纲um,针齿壳上Zr个齿槽圆心的齿槽偏差,简称针齿齿槽偏差Rk,k=1,2,…,Zr、齿距偏差,简称针轮齿距偏差Pk,k=1,2,…,Zr;其中:Rk、Pk分别表示第k个针齿的齿槽偏差、齿距偏差大小,量纲um;
5)曲柄轴的关键误差的定义
曲柄轴的关键误差是指:曲柄轴1、2、3上的偏心凸轮1、2的偏心误差(Ecji,βcji),j=1,2;i=1,2,3;其中:Ecji表示偏心凸轮的偏心误差大小,量纲um;βcji表示偏心凸轮的偏心误差方向,量纲°;
6)行星架的关键误差的定义
行星架的关键误差是指:行星架上曲柄轴孔1、2、3的偏心误差(Ecai,βcai),i=1,2,3,行星架上安装轴承的轴颈偏心误差(Ac,γc);其中:Ecai、Ac表示曲柄轴孔、轴颈的偏心误差大小,量纲um;βcai、γc表示曲柄轴孔、轴颈的偏心误差方向,量纲°;
7)针齿销的关键误差的定义
针齿销的关键误差是指:Zr个针齿销的平均直径误差δzr,量纲um;
8)轴承的关键误差的定义
在三曲柄摆线针轮减速机中,摆线轮1、2与曲柄轴1、2、3间的轴承采用滚子轴承,其轴承间隙分别为δbji,j=1,2;i=1,2,3,量纲um;
行星架与曲柄轴1、2、3间的轴承采用圆锥滚子轴承,其轴承间隙分别为δxi,i=1,2,3,量纲um;
行星架与针齿壳间的轴承采用角接触球轴承,其轴承间隙为δca,量纲um;
S3.建立计算三曲柄摆线针轮减速机动态传动误差的力学模型,具体如下
1)在力学模型中,将零件间的轴承或轴支承、轮齿啮合等接触处用弹簧进行等效,并用轴承或轴刚度、轮齿啮合刚度物理量进行描述;太阳轮与行星轮i,i=1,2,3间的轮齿啮合刚度ki,量纲N/um,i=1,2,3、曲柄轴i,i=1,2,3与摆线轮j,j=1,2轴孔间的轴承支承刚度kji,量纲N/um,j=1,2;i=1,2,3、曲柄轴i,i=1,2,3和行星架轴孔间的轴承支承刚度kbi,量纲N/um,i=1,2,3、行星架与针齿壳间轴承支承刚度kca,N/um、摆线轮与针齿间的啮合刚度kjk,量纲N/um,k=1,2,…,Zr、太阳轮轴的支承刚度ks,量纲N/um;以上刚度的确定方法为:ki采用轮齿弯曲变形来确定;ks通过计算太阳轮轴的弯曲变形来求得;而kji、kbi、kca、kjk接触刚度则采用Palmgren公式确定;
Palmgren公式为:
其中,
式中ν1,ν2——弹性体1、2的泊松比;
E1,E2——弹性体1、2的弹性模量,量纲MPa;
δ——两弹性体间的接触变形量,量纲mm;
l——两弹性体间的接触长度,量纲mm;
F——两弹性体间的载荷,量纲N;
若用k0表示两弹性体间的接触刚度,量纲N/mm,则
F=k0δ (2)
依据式(1)、(2)可推导出两弹性体间的接触刚度k0为:
由式(3)可知,两弹性体间的接触刚度k0主要取决于两弹性体的材料、接触长度及其作用载荷;
2)确定力学模型的静、动坐标系;以太阳轮轴或行星架、针齿壳的理论中心O为原点,垂直减速机轴线的断面为平面静坐标系xoy;从输入端即太阳轮端开始,靠近太阳轮的摆线轮编号j=1,另一个摆线轮编号j=2;以摆线轮的理论质心Oj为原点,以摆线轮j,j=1,2的偏心方向为ηj轴,拟定太阳轮轴逆时针旋转,则将垂直于ηj轴并沿其公转方向转动90°为ξj轴,建立摆线轮j的动坐标系ηjOjξj;在动力学建模过程中,取j=1的摆线轮ηj轴在起始位置时与静坐标系x轴方向一致;
3)在力学模型,将各弹簧设置在各零件坐标轴的正方向接触处或轮齿啮合处,且设定使弹簧受拉为正、受压为负;φi表示摆线轮或行星架上的曲柄轴轴孔的相对位置,取φi=2π(i-1)/3,i=1,2,3,ψj表示2个摆线轮的理论质心Oj的相对位置,取ψj=(j-1)π,j=1,2;
S4.依据S2确定的三曲柄摆线针轮减速机中各零件制造误差、间隙,并结合S3所建立的力学模型,确定传动***在轴或轴承支承处、轮齿啮合处所产生的位移;具体如下:
1)确定太阳轮的基圆偏心误差、装配误差分别在其轮齿啮合处、轴支承处所产生的位移;
若太阳轮的基圆偏心误差为(Es,βs),则其在轮齿啮合处所产生的位移为:
esi=Es cos(θs+βs-Ai) (4)
式中Ai=θc+φi+π/2-α′,θc=θs/iz
α′——太阳轮和行星轮的啮合角,量纲°;
θc——行星架的理论转角,量纲°;
若太阳轮装配误差为(As,γs),则其在太阳轮轴支承处x、y方向上所产生的位移分别为:
esx=As cosγs (5)
esy=As sinγs (6)
2)确定行星轮的基圆偏心误差在其轮齿啮合处所产生的位移;3个行星轮的基圆偏心误差(Epi,βpi)i=1,2,3在轮齿啮合处所产生的位移为:
epi=-Epi cos(βpi-θp-Ai),i=1,2,3 (7)
式中θp——行星轮的理论自转角,量纲°;
3)确定摆线轮上的曲柄轴孔偏心误差在其轴承支承处所产生的位移,及其齿槽偏差、齿距偏差在轮齿啮合处所产生的位移;
定义Oji表示摆线轮j,j=1,2上曲柄轴孔i,i=1,2,3的圆心,Rdc表示摆线轮上曲柄轴孔圆心所在圆周的半径;
若摆线轮曲柄轴孔的偏心误差为:(Ehji,βhji),j=1,2;i=1,2,3,则其在轴承支承处静坐标系x、y方向上所产生的位移分别为:
ehxji=-Ehji cos(θc+φi+βhji) (8)
ehyji=-Ehji sin(θc+φi+βhji) (9)
定义摆线轮的齿槽偏差Rjk、齿距偏差Pjk,O表示针轮中心,Oj表示摆线轮j,j=1,2中心,P表示摆线轮j与针轮的节点;
若摆线轮j上的齿槽偏差为Rjk,j=1,2、齿距偏差为Pjk,j=1,2,则其在轮齿啮合处所产生的位移分别为:
eRjk=Rjkcos(αjk-φdjk) (10)
ePjk=Pjksin(αjk-φdjk) (11)
式中αjk——第k个针齿中心与第j个摆线轮节点连线与ηj轴正方向夹角,量纲°,其中k=1,2,…,Zr;
φdjk——第k个针齿中心与第j个摆线轮中心Oj连线与ηj轴正方向的夹角,量纲°;
4)确定针齿壳上齿槽圆心的齿槽偏差、齿距偏差在摆线轮齿啮合处所产生的位移;
定义针齿壳上齿槽圆心的齿槽偏差,简称针轮齿槽偏差Rk、齿距偏差,简称针轮齿槽偏差Pk,其中:O表示针轮中心,Oj表示摆线轮j,j=1,2中心,P表示摆线轮j与针轮的节点;
若针轮的齿槽偏差为Rk,k=1,2,…,Zr、齿距偏差为Pk,k=1,2,…,Zr,则其在轮齿啮合处所产生的位移分别为:
eRk=-Rk cos(αjk-φjk) (12)
ePk=-Pk sin(αjk-φjk) (13)
式中φjk——第k个针齿半径方向与第j个摆线轮ηj轴正方向间的夹角,量纲°;
5)确定曲柄轴偏心凸轮的偏心误差在其轴承支承处所产生的位移;
定义曲柄轴i的偏心凸轮j偏心误差,Oji表示摆线轮j,j=1,2上曲柄轴孔i,i=1,2,3的圆心,Opj表示曲柄轴i,i=1,2,3上与摆线轮j,j=1,2相对应的曲柄轴回转中心;
若曲柄轴i上偏心凸轮j的偏心误差为(Ecji,βcji),j=1,2;i=1,2,3,则其在x、y方向上所产生的位移分别为:
ecxji=Ecji cos(θp+ψj+βcji) (14)
ecyji=-Ecji sin(θp+ψj+βcji) (15)
6)确定行星架上曲柄轴孔、轴颈的偏心误差在轴承支承处所产生的位移;
定义行星架曲柄轴孔i的偏心误差,其中Oci表示行星架上曲柄轴孔i,i=1,2,3的圆心,Tout表示作用在行星架上的负载转矩;
若行星架曲柄轴孔i的偏心误差为(Ecai,βcai),i=1,2,3,则它在支承处x、y方向上所产生的位移分别为:
ecaxi=-Ecai cos(θc+φi+βcai) (16)
ecayi=-Ecai sin(θc+φi+βcai) (17)
若假设行星架的轴颈偏心误差为(Ac,γc),则它在支承处x、y方向上所产生的位移分别为:
ecx=Ac cosγc (18)
ecy=Ac sinγc (19)
7)确定针齿销的直径误差、针齿壳上的齿槽半径误差在摆线轮齿啮合处形成的间隙,并确定该间隙在摆线轮齿啮合处所产生的位移;
定义O表示针齿销的理论分布圆心或针齿壳上齿槽理论圆心的分布圆圆心,δjk表示摆线轮齿的啮合间隙;
若针齿销的平均直径误差δzr,量纲um、针齿齿槽的平均半径误差δcr,量纲um,则其在摆线轮齿啮合处所产生的间隙δjk及其所产生的位移ejk分别为:
δjk=δzr+δcr (20)
ejk=-δjk (21)
8)各轴承间隙在其支承处所产生的位移
若摆线轮j曲柄轴孔i处的滚子轴承间隙为δbji,j=1,2;i=1,2,3,则它在支承处所产生的位移为:
ebji=-δbji (22)
若行星架在曲柄轴孔i处的圆锥滚子轴承间隙为δxi,i=1,2,3,则它在支承处所产生的位移为:
exi=-δxi (23)
若行星架与针齿壳间的角接触球轴承间隙为δca,则它在支承处所产生的位移为:
eca=-δca (24)
S5.确定三曲柄摆线针轮减速机中各零件的微位移在轴或轴承以及轮齿啮合处所产生的位移;具体如下:
在分析三曲柄摆线针轮减速机中各零件的微位移在接触处所产生的位移时,假设针齿壳固定不动,行星轮与曲柄轴固结在一起,且太阳轮、带曲柄轴的行星轮、摆线轮及行星架的线位移方向与静坐标系或动坐标系的坐标轴正向相同、角位移方向与设计拟定方向相同;
1)太阳轮的微位移在轮齿啮合处、太阳轮轴支承处所产生的位移;
若太阳轮的微位移为xs、ys、θsa-θs,则其在轮齿啮合处所产生的位移为:
si=xs cos Ai+ys sin Ai+Rbs(θsa-θs) (25)
式中Rbs——表示太阳轮的基圆半径,量纲mm;
在轴支承处x、y方向上所产生的位移为:
sx=xs (26)
sy=ys (27)
2)行星轮及曲柄轴的微位移在其轮齿啮合处、摆线轮滚子轴承支承处以及行星架圆锥滚子轴承支承处所产生的位移;
若令行星轮i(i=1,2,3)的微位移为xpi、ypi和θpi-θp,则其在轮齿啮合处所产生的位移为:
spi=-xpi cos Ai-ypi sin Ai-Rbp(θpi-θp) (28)
式中Rbp——行星轮的基圆半径,量纲mm;
θpi——行星轮的实际自转角,量纲°;
在摆线轮滚子轴承支承处x、y方向上所产生的位移为:
sjix=xpi-e(θpi-θp)sin(θp+ψj) (29)
sjiy=ypi-e(θpi-θp)cos(θp+ψj) (30)
式中e——曲柄轴的偏心距,量纲mm;
在行星架圆锥滚子轴承支承处x、y方向上所产生的位移为:
six=xpi (31)
siy=ypi (32)
3)摆线轮的微位移在其轮齿啮合处、曲柄轴滚子轴承支承处所产生的位移;
设摆线轮j,j=1,2的微位移分别为ηdj、θdj-θc和θOj-θp,其中ηdj表示摆线轮的线微位移,θdj-θc表示摆线轮的自转角微位移,θOj-θp表示摆线轮的公转角微位移;
则摆线轮j,j=1,2的微位移在其轮齿啮合处所产生的位移为:
sjk=ηdjcosαjk-Rd(θdj-θc)sinαjk+e(θOj-θp)sinαjk (33)
式中Rd——摆线轮的节圆半径,量纲mm;
在曲柄轴滚子轴承支承处x、y方向上所产生的位移分别为:
scjix=Rdc(θdj-θc)sin(θc+φi)+e(θOj-θp)sin(θp+ψj)-ηdjcos(θp+ψj) (34)
scjiy=-Rdc(θdj-θc)cos(θc+φi)+ηdjsin(θp+ψj)+e(θOj-θp)cos(θp+ψj) (35)
式中Rdc——太阳轮与行星轮间的中心距,量纲mm;
4)行星架的微位移在曲柄轴圆锥滚子轴承支承处、针齿壳角接触球轴承支承处所产生的位移;
设行星架的三个微位移为xca、yca和θca-θc,则其在曲柄轴圆锥滚子轴承支承处x、y方向上所产生的位移分别为:
sbix=-xca+Rdc(θca-θc)sin(θc+φi) (36)
sbiy=-yca-Rdc(θca-θc)cos(θc+φi) (37)
在针齿壳角接触球轴承支承处x、y方向上所产生的位移分别为:
scax=xca (38)
scay=yca (39)
S6.确定三曲柄摆线针轮减速机中,各零件在轴或轴承支承处、轮齿啮合处的作用力;
依据上述步骤S4、S5,确定的各零件制造误差、间隙及微位移在其支承、轮齿啮合处所产生的位移,同时考虑在理想状态时,各零件所传递的理论载荷,从而获得各零件在轴或轴承、轮齿啮合处的作用力;若作用力为正,则表示两物体间的接触力受压,反之受拉
1)太阳轮在其支承处x、y方向上的作用力分别为Fsx,量纲N、Fsy,量纲N,如公式(40)、(41)所示;
Fsx=ks(sx-esx)=ks(xs-esx) (40)
Fsy=ks(sy-esy)=ks(ys-esy) (41)
2)太阳轮与行星轮i,i=1,2,3间的轮齿啮合力为Fi,量纲N,如公式(42)所示;
式中Fi0——表示太阳轮与行星轮间的理论轮齿啮合力,量纲N;
若公式(42)中[]内数值≤0,则表示作用力为拉力,说明太阳轮与行星轮间有间隙存在而没有接触,此时取Fi的值为零;
3)摆线轮j,j=1,2与曲柄轴i,i=1,2,3接触处x、y方向上的作用力分别为Fjix,量纲N、Fjiy,量纲N,如公式(43)、(44)所示;
式中Fjix0——表示理想状态时,摆线轮j,j=1,2与曲柄轴i,i=1,2,3
间的滚子轴承在x方向上所承受的载荷,量纲N;
式中Fjiy0——表示理想状态时,摆线轮j,j=1,2与曲柄轴i,i=1,2,3
间的滚子轴承在y方向上所承受的载荷,量纲N;
公式(43)、(44)中±号是依据摆线轮与曲柄轴接触处的变形量大小进行确定;
4)摆线轮j,j=1,2轮齿与第k,k=1,2,…,Zr个针齿间的作用力为Fjk,量纲N,如公式(45)所示;
式中Fjk0——表示理想状态时摆线轮j,j=1,2的轮齿啮合力,量纲N;
若Fjk≤0,则表示轮齿啮合力为拉力,说明摆线轮轮齿与针齿之间有间隙存在而没有接触,则取Fijk=0;
5)行星架与曲柄轴i,i=1,2,3接触处x、y方向上的作用力分别为Fcix,量纲N、Fciy,量纲N,如公式(46)、(47)所示;
式中Fci0——表示理想状态时行星架与曲柄轴i,i=1,2,3间的作用力,量纲N;
公式(46)、(47)中±号是依据行星架与曲柄轴接触处变形量的大小进行确定;
6)行星架与针齿壳支承处x、y方向上的作用力分别为Fcax,量纲N、Fcay,量纲N,如公式(48)、(49)所示;
Fcax=kca(scax-ecx±eca)=kca(xca-ecx±eca) (48)
Fcay=kca(scay-ecy±eca)=kca(yca-ecy±eca) (49)
公式(48)、(49)中±号是依据行星架与针齿壳支承处变形量的大小来确定;
S7.建立计算三曲柄摆线针轮减速机动态传动误差的数学模型;
通过分析三曲柄摆线针轮减速机中各零件任意位置时的受力状况,依据D'Alembert原理,建立该动态传动误差的数学模型;在建立数学模型时,考虑了摆线轮、行星轮及曲柄轴的科氏Coriolis加速度,忽略了传动***中的摩擦阻力,通过整理,得出三曲柄摆线针轮减速机动态传动误差的数学模型如下:
式中ms——太阳轮的质量,量纲kg;
msp——行星轮及曲柄轴的质量和,量纲kg;
mbx——摆线轮的质量,量纲kg;
Jop——行星轮及曲柄轴的转动惯量,量纲kg·m2;
Joj——摆线轮的转动惯量,量纲kg·m2;
ωc——行星架理论角速度,量纲rad/s;
ωp——行星轮自转理论角速度,量纲rad/s;
上述数学模型可整理成以下矩阵形式:
式中M、C、K——分别为质量阵、阻尼阵及刚度阵,矩阵阶数均为21×21;
X——位移向量,其中X=(xs,ys,θsa,xp1,yp1,θp1,xp2,yp2,θp2,xp3,yp3,θp3,ηd1,θo1,θd1,ηd2,θo2,θd2,xca,yca,θca)T;
Q——广义力矢量;
由于该传动***中的刚度系数kji、kbi、kca、kjk是位移的函数,因此所建立的三曲柄摆线针轮减速机动态传动误差的数学模型属于非线性动力学模型;
S8.采用非线性Newmark法、Wilson-θ法的直接积分法,求解所建立的三曲柄摆线针轮减速机动态传动误差的数学模型;通过编程,并利用计算机进行数值计算,即可确定该减速机输出轴在任意时刻的实际转角θc、动态传动误差Δθc。
2.根据权利要求1所述的机器人用摆线针轮减速机的动态传动精度数值计算方法,其特征是:步骤S2中采用的误差测量仪器具有智能模式识别、自动扫描等功能,并具有测量精度高、灵敏度高、效率高的特点;其主要技术参数的要求:最小读取单位≤0.1um、测量精度不低于1.2+4L/1000,量纲um,L为被测工件长度,量纲为mm、工作台平面尺寸不小于450mm×400mm。
3.根据权利要求1所述的机器人用摆线针轮减速机的动态传动精度数值计算方法,其特征是:步骤S2中采用高精度三坐标测量仪或其他零件误差测量仪器;计算方法采用最小二乘法理论、数学变换与滤波理论以及优化方法的基本理论与方法,测量数据的通过计算机处理技术或方法,从而确定上述各零件的关键误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011392955.8A CN112464481B (zh) | 2020-12-02 | 2020-12-02 | 机器人用摆线针轮减速机的动态传动精度数值计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011392955.8A CN112464481B (zh) | 2020-12-02 | 2020-12-02 | 机器人用摆线针轮减速机的动态传动精度数值计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112464481A CN112464481A (zh) | 2021-03-09 |
CN112464481B true CN112464481B (zh) | 2024-04-16 |
Family
ID=74806418
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011392955.8A Active CN112464481B (zh) | 2020-12-02 | 2020-12-02 | 机器人用摆线针轮减速机的动态传动精度数值计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112464481B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112857794B (zh) * | 2021-03-11 | 2023-04-14 | 苏州普瑞川传动科技有限公司 | 一种齿轮误差快速检测装置 |
CN113792396B (zh) * | 2021-08-25 | 2023-08-22 | 东北林业大学 | 一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101109436A (zh) * | 2006-07-21 | 2008-01-23 | 姜虹 | 适用于动力传动的增速或减速齿轮副 |
CN101705979A (zh) * | 2009-11-27 | 2010-05-12 | 齐中才 | 摆线针轮减速机 |
CN107458598A (zh) * | 2016-06-03 | 2017-12-12 | 贝尔直升机德事隆公司 | 使用叶片桨距固定的马达模块的矩阵的反扭矩控制 |
CN108253094A (zh) * | 2018-01-17 | 2018-07-06 | 河南科技大学 | 一种rv减速器摆线针轮承载啮合印痕确定方法及装置 |
WO2018236254A1 (ru) * | 2017-06-21 | 2018-12-27 | Анатолий Степанович ТОКАРЬ | Трохоидальное зубчатое зацепление и планетарная передача |
WO2020233195A1 (zh) * | 2019-05-22 | 2020-11-26 | 苏州华震工业机器人减速器有限公司 | 精密控制用内啮合减速机 |
-
2020
- 2020-12-02 CN CN202011392955.8A patent/CN112464481B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101109436A (zh) * | 2006-07-21 | 2008-01-23 | 姜虹 | 适用于动力传动的增速或减速齿轮副 |
CN101705979A (zh) * | 2009-11-27 | 2010-05-12 | 齐中才 | 摆线针轮减速机 |
CN107458598A (zh) * | 2016-06-03 | 2017-12-12 | 贝尔直升机德事隆公司 | 使用叶片桨距固定的马达模块的矩阵的反扭矩控制 |
WO2018236254A1 (ru) * | 2017-06-21 | 2018-12-27 | Анатолий Степанович ТОКАРЬ | Трохоидальное зубчатое зацепление и планетарная передача |
CN108253094A (zh) * | 2018-01-17 | 2018-07-06 | 河南科技大学 | 一种rv减速器摆线针轮承载啮合印痕确定方法及装置 |
WO2020233195A1 (zh) * | 2019-05-22 | 2020-11-26 | 苏州华震工业机器人减速器有限公司 | 精密控制用内啮合减速机 |
Non-Patent Citations (2)
Title |
---|
RV减速器摆线轮实体结构参数化设计技术研究;韩林山 等;煤矿机械;20200215(第02期);全文 * |
误差对RV型减速机传动精度的灵敏度研究;袁昕 等;机电产品开发与创新;20090328(第02期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112464481A (zh) | 2021-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112464481B (zh) | 机器人用摆线针轮减速机的动态传动精度数值计算方法 | |
Xu et al. | Dynamic modeling and contact analysis of a cycloid-pin gear mechanism with a turning arm cylindrical roller bearing | |
Xu et al. | Dynamic modelling and contact analysis of bearing-cycloid-pinwheel transmission mechanisms used in joint rotate vector reducers | |
CN112464486B (zh) | 一种机器人用摆线针轮减速机的传动误差数值仿真方法 | |
Ericson et al. | Experimental measurement of the effects of torque on the dynamic behavior and system parameters of planetary gears | |
Tran et al. | Lost motion analysis of one stage cycloid reducer considering tolerances | |
Huang et al. | A study on loaded tooth contact analysis of a cycloid planetary gear reducer considering friction and bearing roller stiffness | |
CN110008610B (zh) | 摆线轮齿廓分段优化设计方法 | |
Li et al. | Mesh stiffness calculation of cycloid-pin gear pair with tooth profile modification and eccentricity error | |
CN111488681B (zh) | 基于不确定性的斜齿轮副随机动力学建模方法 | |
Li et al. | Analysis of contact mechanical characteristics of flexible parts in harmonic gear reducer | |
CN111444607B (zh) | 一种转子-轴承多源激励非线性***建模方法 | |
CN112539721B (zh) | 机器人用减速机的三曲柄摆线轮关键加工误差测量方法 | |
Pham et al. | Hysteresis curve analysis of a cycloid reducer using non-linear spring with a dead zone | |
CN111859576B (zh) | 机器人用rv减速器含间隙机构传动误差计算方法 | |
CN116341105A (zh) | 多源激励下人字齿行星传动***动力学的建模方法 | |
Xu et al. | Design and dynamic transmission error analysis of a new type of cycloidal-pin reducer with a rotatable output-pin mechanism | |
Zheng et al. | Dynamic model-driven intelligent fault diagnosis method for rotary vector reducers | |
CN111881529B (zh) | 机器人用rv减速器主轴承受力分析计算方法 | |
CN108846196B (zh) | 一种2k-v型精密减速机的传动精度优化设计方法 | |
He et al. | Multi-Tooth Contact Analysis and Tooth Profile Modification Optimization for Cycloid Drives in Industrial Robots | |
Wang et al. | Transmission characteristics of an RV reducer based on ADAMS | |
Li et al. | Design and investigation of a cycloid helical gear drive | |
CN115329476A (zh) | 行星滚柱丝杠螺纹牙动载荷分布计算方法 | |
CN109489599B (zh) | 一种自校准的渐开线齿轮副振动位移及动态传动误差测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |