CN108253094A - 一种rv减速器摆线针轮承载啮合印痕确定方法及装置 - Google Patents
一种rv减速器摆线针轮承载啮合印痕确定方法及装置 Download PDFInfo
- Publication number
- CN108253094A CN108253094A CN201810045833.8A CN201810045833A CN108253094A CN 108253094 A CN108253094 A CN 108253094A CN 201810045833 A CN201810045833 A CN 201810045833A CN 108253094 A CN108253094 A CN 108253094A
- Authority
- CN
- China
- Prior art keywords
- cycloidal
- tooth
- wheel
- pinwheel
- pin wheel
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000004048 modification Effects 0.000 claims abstract description 17
- 238000012986 modification Methods 0.000 claims abstract description 17
- 241000283984 Rodentia Species 0.000 claims abstract description 10
- 230000009467 reduction Effects 0.000 claims description 20
- 230000033001 locomotion Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 5
- 241000208340 Araliaceae Species 0.000 claims 4
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 4
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 4
- 235000008434 ginseng Nutrition 0.000 claims 4
- 230000005540 biological transmission Effects 0.000 abstract description 26
- 238000013461 design Methods 0.000 abstract description 7
- 238000005516 engineering process Methods 0.000 abstract description 5
- 238000012937 correction Methods 0.000 abstract description 4
- 230000006872 improvement Effects 0.000 abstract description 3
- 210000000515 tooth Anatomy 0.000 description 99
- 238000004458 analytical method Methods 0.000 description 19
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000005457 optimization Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000000306 component Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 125000000174 L-prolyl group Chemical group [H]N1C([H])([H])C([H])([H])C([H])([H])[C@@]1([H])C(*)=O 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003638 chemical reducing agent Substances 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000001226 toe joint Anatomy 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
- F16H—GEARING
- F16H1/00—Toothed gearings for conveying rotary motion
- F16H1/28—Toothed gearings for conveying rotary motion with gears having orbital motion
- F16H1/32—Toothed gearings for conveying rotary motion with gears having orbital motion in which the central axis of the gearing lies inside the periphery of an orbital gear
-
- 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
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
- F16H—GEARING
- F16H57/00—General details of gearing
- F16H2057/0087—Computer aided design [CAD] specially adapted for gearing features ; Analysis of gear systems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mechanical Engineering (AREA)
- Gears, Cams (AREA)
Abstract
本发明涉及一种RV减速器摆线针轮承载啮合印痕确定方法及装置,属于摆线针轮承载传动技术领域。本发明首先获取摆线针轮空载时轮齿啮合点的位置参数和初始啮合间隙;然后根据初始啮合间隙确定摆线针轮承载时可能接触轮齿的齿号,并由此确定轮齿各啮合位置的接触力、变形量及接触变形总能量;最后确定实际啮合齿数和啮合接触范围,根据实际啮合齿对的载荷分配,计算轮齿啮合的变形量,最终得到RV减速器摆线针轮承载啮合印痕。本发明考虑了摆线针轮的相位差、修形设计以及多齿同时啮合的误差均化对RV传动的综合影响,得到啮合印痕负荷真实工况,对RV减速器的承载啮合特性、传动精度及精度保持性的改善、预控和应用有重要意义。
Description
技术领域
本发明涉及一种RV减速器摆线针轮承载啮合印痕确定方法及装置,属于摆线针轮承载传动技术领域。
背景技术
在机器人的四大核心部件中,RV(Rotary Vector)减速器最为关键,RV减速是在摆线针轮传动的基础上发展起来的一种二级减速传动装置,采用共用曲柄轴和中心圆盘支撑的结构形式组成封闭式行星传动,具有传动比大、体积小、重量轻、承载能力强、刚度大等优点。在实际工程应用中,RV减速器在承载工况下的高精传动是机器人准确运动和重复定位的基本要求。但与国外相比,我国RV减速器在空载下传动精度和精度寿命等指标良好,而承载下则出现传动误差大、精度寿命低等共性问题,严重影响了国产机器人的运动准确性和可靠性,经常造成机械手动作不到位,甚至导致严重事故。
摆线针轮行星传动与RV减速器输出轴直接相连,对运动精度指标影响最大,其中,摆线轮的修形设计与摆线针轮的传动误差控制是核心技术。国外经过几十年的努力和突破,目前己成功解决了摆线针轮传动中的修形设计、参数优化及运动精度指标预控等技术难题。由于RV减速器的整个传动过程中不仅存在零部件的弹性变形,而且存在复杂的接触关系,接触位置也会随着变形、载荷、材料等因素的改变而改变,导致难以建立其完善的理论模型。
国内专家学者目前主要采用有限元模拟仿真的方法进行RV传动的啮合特性进行分析,例如名称为《RV减速器综合性能实验与仿真》的论文,该论文给出在RV传动机构特性分析的基础上,利用Pro/E软件建立RV减速器中各零部件三维实体模型,进行运动学仿真,借助有限元分析软件对摆线轮与针齿之间的非线性接触问题进行受力分析,但是其模拟结果和分析方法缺少相应的理论支撑,仿真模拟结果与实际工况差异较大,轮齿的实际承载啮合分析方法也不够完善,对RV传动精度的改善和控制的指导意义不大。另外,RV传动的主要特征之一是多齿同时啮合,其误差均化效应明显有利于提高传动精度和承载能力,但多齿啮合传动系统的载荷分布和实际啮合状况更加复杂,在啮合特性分析时,必须同时考虑RV传动中摆线针轮的相位差、修形设计对承载传动特性和啮合接触的非线性关系以及多齿同时啮合接触的误差均化影响。而采用有限元进行分析时,并没有考虑上述影响,导致得到RV减速器摆线针轮承载啮合印痕与RV摆线针轮传动的实际啮合印痕并不一致。
发明内容
本发明的目的是提供一种RV减速器摆线针轮承载啮合印痕确定方法,以解决目前采用有限元分析所获取的RV减速器摆线针轮承载啮合印痕与实际啮合印痕不一致的问题;本发明还提供了一种RV减速器摆线针轮承载啮合印痕确定装置。
本发明为解决上述技术问题而提供一种RV减速器摆线针轮承载啮合印痕确定方法,包括以下方案,方法方案一:该确定方法包括以下步骤:
1)根据RV减速器摆线针轮的几何参数及其相对运动关系,确定摆线轮和针轮的齿面及法线的矢量表达;
2)根据摆线针轮的不同相位差以及多齿啮合对RV传动关系的影响,确定空载时摆线针轮轮齿初始啮合点位置参数,并根据该初始啮合点位置参数确定初始啮合间隙;
3)根据初始啮合间隙初步确定可能的啮合齿数和啮合接触范围,计算啮合齿对的变形总能量;
4)基于能量最小原理获得实际啮合齿对的载荷分配、实际啮合齿数和啮合接触范围,计算摆线针轮各啮合位置的变形量,最终确定承载时摆线针轮的啮合印痕。
本发明考虑了摆线针轮的相位差影响、修形设计对承载传动特性和啮合接触的非线性关系以及多齿同时啮合的误差均化对RV传动的综合影响,能够精确获得RV减速器摆线针轮行星传动在拟真实工况下的承载传动误差和啮合印痕,对于RV减速器的承载啮合特性、传动精度及精度保持性的改善和控制具有重要意义。
方法方案二:在方法方案一的基础上,所述步骤1)确定的针轮坐标系Sp中的齿面矢量rp和单位法矢量np为:
摆线轮坐标系Sc中的齿面矢量rc和单位法矢量nc为:
其中下角标p指齿面矢量和单位法矢量表示在坐标系Sp中;β和bpx分别为针轮的齿面参数,bpx≤bp,bp为针齿宽;rp为针轮分布圆半径;rrp为针齿半径;i表示第i个针齿号;i=1~zp,zp为针轮齿数;下角标c指齿面矢量和单位法矢量表示在坐标系Sc中;α和bcx分别为摆线轮的齿廓参数,bcx≤bc,bc为摆线轮齿宽;zc为摆线轮齿数;a为偏心距;iH=zp/zc;Δrp为移距修形量,当rp减小时,Δrp为负,反之为正;Δrrp为等距修形量,当rrp减小时,Δrrp为负,反之为正;k1=azp/(rp+Δrp)。
方法方案三:在方法方案一的基础上,所述步骤2)中的空载时摆线针轮轮齿初始啮合点位置参数是根据下式确定:
其中rfp和nfp分别为针轮的齿面矢量单位法矢量表示在与RV减速器箱体固连的坐标系Sf中;其中rfc和nfc分别为摆线轮的齿面矢量单位法矢量表示在与RV减速器箱体固连的坐标系Sf中;β和bpx分别为针轮的齿面参数,bpx≤bp,bp为针齿宽;α和bcx分别为摆线轮的齿廓参数,bcx≤bc,bc为摆线轮齿宽;φ1和φ2分别表示针轮和摆线轮的转角。
方法方案四:在方法方案一的基础上,所述的步骤2)得到的初始啮合间隙为:
其中(xpi,ypi)表示第i个针齿中心点坐标,(xci,yci)为摆线针轮各齿对处公法线与摆线轮齿廓交点的坐标,φ1为针轮的转角,rp为针轮分布圆半径,rrp为针齿半径,zp为针轮齿数。
方法方案五:在方法方案一的基础上,所述步骤3)的实现过程如下:
A.根据初始啮合间隙确定可能接触轮齿的齿号,根据赫兹接触理论建立载荷与变形量之间的非线性方程;
B.基于承载下摆线针轮的多齿啮合理论,由变形协调条件可知各啮合点处的转角误差相同,利用精确啮合间隙、变形量和力臂大小获得转角误差;
C.根据能量守恒原理建立各啮合点处的载荷与输出力矩之间的关系计算出摆线针轮各啮合点处的载荷分配大小及接触变形量,以及该啮合状态下的接触变形总能量TE;
D.计算摆线针轮在不同接触轮齿情况下的接触变形总能量,根据最小能量原理确定得到轮齿的实际承载啮合的齿数及啮合范围。
本发明还提供了一种RV减速器摆线针轮承载啮合印痕确定装置,包括以下方案,装置方案一:该确定装置包括存储器和处理器,以及存储在所述存储器上并在所述处理器上运行的计算机程序,所述处理器与所述存储器相耦合,所述处理器执行所述计算机程序时实现以下步骤:
1)根据RV减速器摆线针轮的几何参数及其相对运动关系,确定摆线轮和针轮的齿面及法线的矢量表达;
2)根据摆线针轮的不同相位差以及多齿啮合对RV传动关系的影响,确定空载时摆线针轮轮齿初始啮合点位置参数,并根据该初始啮合点位置参数确定初始啮合间隙;
3)根据初始啮合间隙初步确定可能的啮合齿数和啮合接触范围,计算啮合齿对的变形总能量;
4)基于能量最小原理获得实际啮合齿对的载荷分配、实际啮合齿数和啮合接触范围,计算摆线针轮各啮合位置的变形量,最终确定承载时摆线针轮的啮合印痕。
装置方案二:在装置方案一的基础上,所述步骤1)确定的针轮坐标系Sp中的齿面矢量rp和单位法矢量np为:
摆线轮坐标系Sc中的齿面矢量rc和单位法矢量nc为:
其中下角标p指齿面矢量和单位法矢量表示在坐标系Sp中;β和bpx分别为针轮的齿面参数,bpx≤bp,bp为针齿宽;rp为针轮分布圆半径;rrp为针齿半径;i表示第i个针齿号;i=1~zp,zp为针轮齿数;下角标c指齿面矢量和单位法矢量表示在坐标系Sc中;α和bcx分别为摆线轮的齿廓参数,bcx≤bc,bc为摆线轮齿宽;zc为摆线轮齿数;a为偏心距;iH=zp/zc;Δrp为移距修形量,当rp减小时,Δrp为负,反之为正;Δrrp为等距修形量,当rrp减小时,Δrrp为负,反之为正;k1=azp/(rp+Δrp)。
装置方案三:在装置方案一的基础上,所述步骤2)中的空载时摆线针轮轮齿初始啮合点位置参数是根据下式确定:
其中rfp和nfp分别为针轮的齿面矢量单位法矢量表示在与RV减速器箱体固连的坐标系Sf中;其中rfc和nfc分别为摆线轮的齿面矢量单位法矢量表示在与RV减速器箱体固连的坐标系Sf中;β和bpx分别为针轮的齿面参数,bpx≤bp,bp为针齿宽;α和bcx分别为摆线轮的齿廓参数,bcx≤bc,bc为摆线轮齿宽;φ1和φ2分别表示针轮和摆线轮的转角。
装置方案四:在装置方案一的基础上,所述的步骤2)得到的初始啮合间隙为:
其中(xpi,ypi)表示第i个针齿中心点坐标,(xci,yci)为摆线针轮各齿对处公法线与摆线轮齿廓交点的坐标,φ1为针轮的转角,rp为针轮分布圆半径,rrp为针齿半径,zp为针轮齿数。
装置方案五:在装置方案一的基础上,所述步骤3)的实现过程如下:
A.根据初始啮合间隙确定可能接触轮齿的齿号,根据赫兹接触理论建立载荷与变形量之间的非线性方程;
B.基于承载下摆线针轮的多齿啮合理论,由变形协调条件可知各啮合点处的转角误差相同,利用精确啮合间隙、变形量和力臂大小获得转角误差;
C.根据能量守恒原理建立各啮合点处的载荷与输出力矩之间的关系计算出摆线针轮各啮合点处的载荷分配大小及接触变形量,以及该啮合状态下的接触变形总能量TE;
D.计算摆线针轮在不同接触轮齿情况下的接触变形总能量,根据最小能量原理确定得到轮齿的实际承载啮合的齿数及啮合范围。
附图说明
图1是RV减速器摆线针轮行星传动机构示意图;
图2是RV减速器摆线针轮行星传动的双啮合运动模型示意图;
图3是RV减速器摆线针轮承载双啮合流程图;
图4是RV减速器摆线针轮空载双啮合传动误差曲线图;
图5是摆线针轮的一个啮合周期内的空载啮合印痕示意图;
图6是摆线针轮传动各齿对处的精确初始啮合间隙曲线图;
图7是RV减速器摆线针轮承载传动误差曲线图;
图8是RV减速器摆线针轮承载传动啮合印痕示意图(某时刻各啮合印痕在摆线同一齿面上表示)。
具体实施方式
下面结合附图对本发明的具体实施方式做进一步的说明。
本发明所针对的RV减速器,包括两个摆线轮和一个针轮,如图1所示,分别为第一摆线轮和第二摆线轮,本发明针对该RV减速器提供了一种RV减速器摆线针轮承载啮合印痕确定方法,该方法首先基于RV减速器摆线针轮行星传动运动关系,得到摆线轮和针轮的齿面及法线的矢量表达;然后构建摆线针轮空载双啮合分析模型,计算得到轮齿啮合点的位置参数,并基于此啮合位置获得其余轮齿的精确初始啮合间隙;再基于变形协调条件和力平衡方程,考虑修形设计对承载传动特性以及啮合接触的非线性关系,建立摆线针轮承载双啮合分析模型,并根据初始啮合间隙初步确定可能接触轮齿的齿号,求解得到轮齿各啮合位置的接触力、变形量及接触变形总能量;最后基于能量最小原理确定实际啮合齿数和啮合接触范围,根据实际啮合齿对的载荷分配,计算轮齿啮合的变形量,最终精确获得承载传动误差曲线和啮合印痕。
本发明所针对的RV减速器包括两个或两个以上的摆线轮,下面以某具体RV减速器的摆线针轮传动为例对本发明的实现过程进行说明,该实例中的RV减速器包括两个摆线轮,分别为第一摆线轮和第二摆线轮,该RV减速器中的摆线针轮几何参数如表1所示。
表1
参数项目 | 摆线轮和针轮 |
针轮分布半径rp/mm | 52.15 |
针齿半径rrp/mm | 2 |
针齿数zp | 40 |
摆线齿数zc | 39 |
偏心距a/mm | 0.9 |
摆线轮齿宽bc/mm | 10 |
针齿宽bp/mm | 22 |
等距修形量Δrrp/mm | 0.01 |
移距修形量Δrp/mm | -0.012 |
1.推导摆线针轮的齿面矢量和单位法线矢量。
根据图2建立的RV减速器摆线针轮行星传动的啮合运动模型,运用微分几何原理,推导得到针轮坐标系Sp中的齿面矢量rp和单位法矢量np:
其中,下角标p表示坐标系Sp;β和bpx分别为针轮的齿廓参数(bpx≤bp,bp为针齿宽);rp为针轮分布圆半径;rrp为针齿半径;i表示第i个针齿号;zp为针轮齿数。
同样推导得到摆线轮坐标系Sc中的齿面矢量rc和单位法矢量nc:
式中,下角标c表示坐标系Sc;α和bcx分别为摆线轮的齿廓参数(bcx≤bc,bc为摆线轮齿宽);zc为摆线轮齿数;a为偏心距;iH=zp/zc;Δrp为移距修形量(当rp减小时,Δrp为负,反之为正);Δrrp为等距修形量(当rrp减小时,Δrrp为负,反之为正);k1=azp/(rp+Δrp)。
根据空间坐标变换原理,将针轮的坐标矢量rp与法矢np、第一摆线轮和第二摆线轮的齿面矢量rc与法矢nc分别表示在与RV减速器箱体固连的坐标系Sf。
式中,φ1、φ2和φ2'分别表示针轮、第一摆线轮和第二摆线轮的转角。
2.建立摆线针轮空载双啮合分析模型,计算单齿啮合点的位置参数。
根据空间坐标转换原理,建立摆线针轮空载双啮合分析模型,该模型由方程组(a)和方程组(b)组成的
其中方程组(a)包含5个独立的非线性方程,6个未知量:φ1、φ2、α1、β1、bpx1和bcx1,给定输入角φ1,方程组(1)可解;方程组(b)也有6个未知量:φ1、φ2'、α2、β2、bpx2和bcx2,φ1已知,故方程组(2)也可解。运用牛顿拉菲森数值优化方法对双方程组进行求解,得到不同φ1转角所对应的第一摆线轮和第二摆线轮与针轮的啮合点的位置参数(φ1,φ2,α1,bcx1,β1,bpx1,N1)和(φ1,φ2',α2,bcx2,β2,bpx2,N2)。其中,N1和N2分别为针轮与第一摆线轮和第二摆线轮啮合的针齿号。
将优化求解得到的啮合点位置参数(φ1,φ2,α1,bcx1,β1,bpx1,N1)和(φ1,φ2',α2,bcx2,β2,bpx2,N2),分别带入方程(3)、(4)和(5),得到第一摆线轮和第二摆线轮与针轮所对应的啮合位置坐标,显示化处理后即可得到摆线针轮的空载啮合印痕。由于空载摆线针轮副是单齿啮合,瞬时啮合印痕为一条直线,将一个啮合周期内的空载啮合印痕表示在摆线轮的同一齿面上,如图5所示。
根据摆线针轮空载双啮合分析模型,给定针轮转角φ1=0,分别令N1=1~zp、N2=1~zp,求解轮齿双啮合分析方程组(6),获得第一摆线轮和第二摆线轮的最小相对啮合转角φ2min和φ2min',对应初始啮合针齿号分别为m和m',将两初始啮合点位置记为(φ1m,φ2m)和(φ1m,φ2m'),在区间[φ1m,φ1m+2π/zp]内,给定针轮转角一定步长Δφ1,循环求解每一针轮转角对应的摆线轮转角φ2和φ2',根据δ=(φ2-φ2m)-iH(φ1-φ1m)和δ'=(φ2'-φ2m')-iH(φ1-φ1m)最终获得空载双啮合传动误差曲线。本实例中的空载双啮合传动误差曲线如图4所示。空载双啮合传动误差曲线由所有齿对共同形成,呈周期性变化,周期为2π/zp。
3.计算各轮齿间的精确初始啮合间隙。
轮齿间的精确初始啮合间隙为摆线针轮实际齿廓处公法线方向上齿廓间的距离,以摆线轮1与针轮啮合为例,由步骤2求解得到的初始啮合点位置参数:针轮转角φ1=φ1m、第一摆线轮转角φ2=φ2m、啮合点处的针齿号m,则第i个针齿中心点Ai坐标可表示为:
摆线针轮各齿对处的公法线方程表示为:
Gi(αi)=ncy·(x-rcx)-ncx·(y-rcy) (8)
式中,αi表示摆线轮齿面位置参数;x和y分别是啮合点坐标;rcx和rcy分别为方程(2)所示摆线轮齿面矢量的坐标分量,ncx和ncy分别为摆线轮法线的坐标分量。
各针齿中心Ai(xpi,ypi)为对应公法线上的一点,将Ai(xpi,ypi)代入公法线方程(8)中,计算公法线与摆线轮齿廓交点Bi所对应的位置参数αi,然后代入方程(2)可得到交点坐标Bi(xci,yci);最后,获得初始啮合间隙表示为:
若求解任一转角下的精确啮合间隙,关键是确定任一转角下轮齿啮合点位置参数。给定针轮转角φ1,由步骤2求得的综合传动误差曲线确定对应的摆线轮转角φ2,根据轮齿双啮合分析方程组(6)获得啮合点的其余参数值,利用方程(7)、(8)和(9),可计算获得任一转角下的精确啮合间隙。本实例中各齿对处初始精确啮合间隙曲线如图6所示。
4.建立摆线针轮承载双啮合分析模型。
在承载工况下,摆线针轮行星传动的转角误差由啮合间隙和接触力所产生的接触变形共同导致,本发明以第一摆线轮和第二摆线轮与针齿的第i对轮齿啮合点处的接触力F1i、F2i为未知量(单位为N),根据赫兹接触理论、变形协调条件、力平衡方程及总能量方程,建立摆线针轮传动的承载双啮合分析模型。
其中,下角标1和2分别表示第一摆线轮和第二摆线轮;b1i和b2i为第i对摆线针轮啮合点处的接触半宽;ρi为第i对摆线针轮啮合点处的当量曲率半径,单位mm,μ1,μ2为分别为针轮和摆线轮的泊松比,μ1=μ2=0.3;E1,E2分别为针轮和摆线轮的弹性模量,E1=E2=206GPa;δ1i和δ2i为第i对摆线针轮啮合点处的弹性变形量,单位mm;d1i和d2i为第i对摆线针轮啮合点处的精确啮合间隙;L1i和L2i为第i对摆线针轮啮合点处的力臂大小,单位mm;s1i和s2i分别为第i对摆线针轮啮合点处的转角,单位rad;Tc为输出转矩Tc=72Nm;TE为第一摆线轮和第二摆线轮与针齿的第i对轮齿啮合的接触变形总能量。
该模型是在不考虑摆线针轮制造误差和安装误差等前提下建立的,摆线针轮多齿同时啮合所产生的误差均化效应使得F1i=F2i,根据步骤3计算得到摆线针轮某一转角下的精确初始啮合间隙,并由啮合间隙初步判断可能接触轮齿的齿号,根据赫兹接触理论建立载荷与变形量之间的非线性方程(式(10)中的方程(b)和(c));基于承载下摆线针轮的多齿啮合理论,由变形协调条件可知各啮合点处的转角误差相同,利用精确啮合间隙、变形量和力臂大小获得转角误差方程式(式(10)中的方程(d)和方程(e));根据能量守恒远离建立各啮合点处的载荷与输出力矩之间的关系(式(10)中方程(f)),由此计算出摆线针轮各啮合点处的载荷分配大小及接触变形量;根据式(10)中方程(g)计算该啮合状态下的接触变形总能量TE。然后,按照此过程可得到摆线针轮在不同接触轮齿情况下的接触变形总能量,最后根据最小能量原理获得轮齿的实际承载啮合齿数及啮合范围等参数。
5.获取承载传动误差曲线。
给定摆线轮转角φ2=0,用承载双啮合接触分析方法计算针轮的最小相对啮合转角φ1min=0.000519rad,进行在(φ1min,0)时刻摆线针轮的承载啮合分析。首先假设有n对齿参与啮合,根据图6各齿对处初始精确啮合间隙曲线,确定假设可能啮合接触的针齿号,以变形能最小为目标,计算出每种假设下的能量,Nmin时为实际啮合情况,确定了啮合齿数为n=17,啮合区间为第23齿到第39齿;将各参数代入式(10)的承载双啮合分析模型中,解非线性方程组获得摆线针轮啮合点处的力Fi(i=1~17),进一步计算出针轮转角φ1=φ1min时的转角误差;针轮转角φ1增加一定步长Δφ1,摆线轮转角φ2=Δφ1×zp/zc,带入式(10)的承载双啮合分析模型中,计算出针轮转角φ1=φ1min+Δφ1时的转角误差。重复上述步骤,获得一系列的转角误差,进而得到摆线针轮的承载传动误差曲线,如图7所示。
6.获取承载啮合印痕。
以初始时刻(φ1min,0)为例,得到的摆线针轮的啮合点位置和各啮合点出的承载接触半宽,如表2所示,进一步可计算出各齿对处一系列瞬时啮合点的位置和接触半宽,显示化处理后即可得到RV摆线针轮承载啮合印痕。
RV传动为多齿同时啮合,在某时刻下有多个摆线轮齿的多处位置存在啮合印痕,为了更清晰的表示啮合印痕的位置,可在摆线轮的同一齿面上表示出所有啮合轮齿的印痕,如图8所示。
表2
同样,该方法也能够对两个以上的摆线轮的RV减速器进行摆线针轮承载啮合印痕的确定,其过程与上述过程相同。
上述方法可以作为一种计算机程序,存储在RV减速器摆线针轮承载啮合印痕确定装置中的存储器中并可在RV减速器摆线针轮承载啮合印痕确定装置中的处理器上运行。
本发明考虑了摆线针轮的相位差的影响及修形设计对承载传动特性以及轮齿接触的非线性关系,同时考虑了RV传动中两个摆线轮与针轮同时接触的误差均化影响,所得到的RV减速器摆线针轮承载传动误差和承载啮合曲线,能够精确反映RV减速器摆线针轮在拟真实工况下的实际啮合状况。
Claims (10)
1.一种RV减速器摆线针轮承载啮合印痕确定方法,其特征在于,该确定方法包括以下步骤:
1)根据RV减速器摆线针轮的几何参数及其相对运动关系,确定摆线轮和针轮的齿面及法线的矢量表达;
2)根据摆线针轮的不同相位差以及多齿啮合对RV传动关系的影响,确定空载时摆线针轮轮齿初始啮合点位置参数,并根据该初始啮合点位置参数确定初始啮合间隙;
3)根据初始啮合间隙初步确定可能的啮合齿数和啮合接触范围,计算啮合齿对的变形总能量;
4)基于能量最小原理获得实际啮合齿对的载荷分配、实际啮合齿数和啮合接触范围,计算摆线针轮各啮合位置的变形量,最终确定承载时摆线针轮的啮合印痕。
2.根据权利要求1所述的RV减速器摆线针轮承载啮合印痕确定方法,其特征在于,所述步骤1)确定的针轮坐标系Sp中的齿面矢量rp和单位法矢量np为:
摆线轮坐标系Sc中的齿面矢量rc和单位法矢量nc为:
其中下角标p指齿面矢量和单位法矢量表示在坐标系Sp中;β和bpx分别为针轮的齿面参数,bpx≤bp,bp为针齿宽;rp为针轮分布圆半径;rrp为针齿半径;i表示第i个针齿号;i=1~zp,zp为针轮齿数;下角标c指齿面矢量和单位法矢量表示在坐标系Sc中;α和bcx分别为摆线轮的齿廓参数,bcx≤bc,bc为摆线轮齿宽;zc为摆线轮齿数;a为偏心距;iH=zp/zc;Δrp为移距修形量,当rp减小时,Δrp为负,反之为正;Δrrp为等距修形量,当rrp减小时,Δrrp为负,反之为正;k1=azp/(rp+Δrp)。
3.根据权利要求1所述的RV减速器摆线针轮承载啮合印痕确定方法,其特征在于,所述步骤2)中的空载时摆线针轮轮齿初始啮合点位置参数是根据下式确定:
其中rfp和nfp分别为针轮的齿面矢量单位法矢量表示在与RV减速器箱体固连的坐标系Sf中;其中rfc和nfc分别为摆线轮的齿面矢量单位法矢量表示在与RV减速器箱体固连的坐标系Sf中;β和bpx分别为针轮的齿面参数,bpx≤bp,bp为针齿宽;α和bcx分别为摆线轮的齿廓参数,bcx≤bc,bc为摆线轮齿宽;φ1和φ2分别表示针轮和摆线轮的转角。
4.根据权利要求1所述的RV减速器摆线针轮承载啮合印痕确定方法,其特征在于,所述的步骤2)得到的初始啮合间隙为:
其中(xpi,ypi)表示第i个针齿中心点坐标,(xci,yci)为摆线针轮各齿对处公法线与摆线轮齿廓交点的坐标,φ1为针轮的转角,rp为针轮分布圆半径,rrp为针齿半径,zp为针轮齿数。
5.根据权利要求1所述的RV减速器摆线针轮承载啮合印痕确定方法,其特征在于,所述步骤3)的实现过程如下:
A.根据初始啮合间隙确定可能接触轮齿的齿号,根据赫兹接触理论建立载荷与变形量之间的非线性方程;
B.基于承载下摆线针轮的多齿啮合理论,由变形协调条件可知各啮合点处的转角误差相同,利用精确啮合间隙、变形量和力臂大小获得转角误差;
C.根据能量守恒原理建立各啮合点处的载荷与输出力矩之间的关系计算出摆线针轮各啮合点处的载荷分配大小及接触变形量,以及该啮合状态下的接触变形总能量TE;
D.计算摆线针轮在不同接触轮齿情况下的接触变形总能量,根据最小能量原理确定得到轮齿的实际承载啮合的齿数及啮合范围。
6.一种RV减速器摆线针轮承载啮合印痕确定装置,其特征在于,该确定装置包括存储器和处理器,以及存储在所述存储器上并在所述处理器上运行的计算机程序,所述处理器与所述存储器相耦合,所述处理器执行所述计算机程序时实现以下步骤:
1)根据RV减速器摆线针轮的几何参数及其相对运动关系,确定摆线轮和针轮的齿面及法线的矢量表达;
2)根据摆线针轮的不同相位差以及多齿啮合对RV传动关系的影响,确定空载时摆线针轮轮齿初始啮合点位置参数,并根据该初始啮合点位置参数确定初始啮合间隙;
3)根据初始啮合间隙初步确定可能的啮合齿数和啮合接触范围,计算啮合齿对的变形总能量;
4)基于能量最小原理获得实际啮合齿对的载荷分配、实际啮合齿数和啮合接触范围,计算摆线针轮各啮合位置的变形量,最终确定承载时摆线针轮的啮合印痕。
7.根据权利要求6所述的RV减速器摆线针轮承载啮合印痕确定装置,其特征在于,所述步骤1)确定的针轮坐标系Sp中的齿面矢量rp和单位法矢量np为:
摆线轮坐标系Sc中的齿面矢量rc和单位法矢量nc为:
其中下角标p指齿面矢量和单位法矢量表示在坐标系Sp中;β和bpx分别为针轮的齿面参数,bpx≤bp,bp为针齿宽;rp为针轮分布圆半径;rrp为针齿半径;i表示第i个针齿号;i=1~zp,zp为针轮齿数;下角标c指齿面矢量和单位法矢量表示在坐标系Sc中;α和bcx分别为摆线轮的齿廓参数,bcx≤bc,bc为摆线轮齿宽;zc为摆线轮齿数;a为偏心距;iH=zp/zc;Δrp为移距修形量,当rp减小时,Δrp为负,反之为正;Δrrp为等距修形量,当rrp减小时,Δrrp为负,反之为正;k1=azp/(rp+Δrp)。
8.根据权利要求6所述的RV减速器摆线针轮承载啮合印痕确定装置,其特征在于,所述步骤2)中的空载时摆线针轮轮齿初始啮合点位置参数是根据下式确定:
其中rfp和nfp分别为针轮的齿面矢量单位法矢量表示在与RV减速器箱体固连的坐标系Sf中;其中rfc和nfc分别为摆线轮的齿面矢量单位法矢量表示在与RV减速器箱体固连的坐标系Sf中;β和bpx分别为针轮的齿面参数,bpx≤bp,bp为针齿宽;α和bcx分别为摆线轮的齿廓参数,bcx≤bc,bc为摆线轮齿宽;φ1和φ2分别表示针轮和摆线轮的转角。
9.根据权利要求6所述的RV减速器摆线针轮承载啮合印痕确定装置,其特征在于,所述的步骤2)得到的初始啮合间隙为:
其中(xpi,ypi)表示第i个针齿中心点坐标,(xci,yci)为摆线针轮各齿对处公法线与摆线轮齿廓交点的坐标,φ1为针轮的转角,rp为针轮分布圆半径,rrp为针齿半径,zp为针轮齿数。
10.根据权利要求6所述的RV减速器摆线针轮承载啮合印痕确定装置,其特征在于,所述步骤3)的实现过程如下:
A.根据初始啮合间隙确定可能接触轮齿的齿号,根据赫兹接触理论建立载荷与变形量之间的非线性方程;
B.基于承载下摆线针轮的多齿啮合理论,由变形协调条件可知各啮合点处的转角误差相同,利用精确啮合间隙、变形量和力臂大小获得转角误差;
C.根据能量守恒原理建立各啮合点处的载荷与输出力矩之间的关系计算出摆线针轮各啮合点处的载荷分配大小及接触变形量,以及该啮合状态下的接触变形总能量TE;
D.计算摆线针轮在不同接触轮齿情况下的接触变形总能量,根据最小能量原理确定得到轮齿的实际承载啮合的齿数及啮合范围。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810045833.8A CN108253094B (zh) | 2018-01-17 | 2018-01-17 | 一种rv减速器摆线针轮承载啮合印痕确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810045833.8A CN108253094B (zh) | 2018-01-17 | 2018-01-17 | 一种rv减速器摆线针轮承载啮合印痕确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108253094A true CN108253094A (zh) | 2018-07-06 |
CN108253094B CN108253094B (zh) | 2019-08-20 |
Family
ID=62740963
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810045833.8A Active CN108253094B (zh) | 2018-01-17 | 2018-01-17 | 一种rv减速器摆线针轮承载啮合印痕确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108253094B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108730480A (zh) * | 2018-07-27 | 2018-11-02 | 河南科技大学 | 一种rv减速器的摆线轮及其齿廓逆向主动修形方法 |
CN111075920A (zh) * | 2020-01-13 | 2020-04-28 | 重庆大学 | 基于fft和润滑影响的rv减速器摆线针轮残余应力求解方法 |
CN111859576A (zh) * | 2020-07-27 | 2020-10-30 | 大连交通大学 | 机器人用rv减速器含间隙机构传动误差计算方法 |
CN112464481A (zh) * | 2020-12-02 | 2021-03-09 | 华北水利水电大学 | 机器人用摆线针轮减速机的动态传动精度数值计算方法 |
CN115659549A (zh) * | 2022-11-14 | 2023-01-31 | 浙江大学 | 一种内曲线液压马达凸轮环曲线修形方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002126844A (ja) * | 2000-08-14 | 2002-05-08 | Mitsubishi Materials Corp | 回転塑性加工の数値シミュレーション方法及び記録媒体及びプログラム |
CN104408220A (zh) * | 2014-10-08 | 2015-03-11 | 西北工业大学 | 一种改进的轮齿加载接触分析方法 |
EP3073396A1 (en) * | 2015-03-24 | 2016-09-28 | Rolls-Royce plc | A method of modelling at least a part of a gas turbine engine |
CN107081678A (zh) * | 2017-07-04 | 2017-08-22 | 河南科技大学 | 一种摆线轮成形磨削砂轮修整方法 |
CN107256282A (zh) * | 2017-05-08 | 2017-10-17 | 华南理工大学 | 一种基于形变量补偿的rv减速器摆线轮齿廓修形方法 |
-
2018
- 2018-01-17 CN CN201810045833.8A patent/CN108253094B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002126844A (ja) * | 2000-08-14 | 2002-05-08 | Mitsubishi Materials Corp | 回転塑性加工の数値シミュレーション方法及び記録媒体及びプログラム |
CN104408220A (zh) * | 2014-10-08 | 2015-03-11 | 西北工业大学 | 一种改进的轮齿加载接触分析方法 |
EP3073396A1 (en) * | 2015-03-24 | 2016-09-28 | Rolls-Royce plc | A method of modelling at least a part of a gas turbine engine |
CN107256282A (zh) * | 2017-05-08 | 2017-10-17 | 华南理工大学 | 一种基于形变量补偿的rv减速器摆线轮齿廓修形方法 |
CN107081678A (zh) * | 2017-07-04 | 2017-08-22 | 河南科技大学 | 一种摆线轮成形磨削砂轮修整方法 |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108730480A (zh) * | 2018-07-27 | 2018-11-02 | 河南科技大学 | 一种rv减速器的摆线轮及其齿廓逆向主动修形方法 |
CN111075920A (zh) * | 2020-01-13 | 2020-04-28 | 重庆大学 | 基于fft和润滑影响的rv减速器摆线针轮残余应力求解方法 |
CN111075920B (zh) * | 2020-01-13 | 2023-02-21 | 重庆大学 | 基于fft和润滑影响的rv减速器摆线针轮残余应力求解方法 |
CN111859576A (zh) * | 2020-07-27 | 2020-10-30 | 大连交通大学 | 机器人用rv减速器含间隙机构传动误差计算方法 |
CN111859576B (zh) * | 2020-07-27 | 2024-02-02 | 大连交通大学 | 机器人用rv减速器含间隙机构传动误差计算方法 |
CN112464481A (zh) * | 2020-12-02 | 2021-03-09 | 华北水利水电大学 | 机器人用摆线针轮减速机的动态传动精度数值计算方法 |
CN112464481B (zh) * | 2020-12-02 | 2024-04-16 | 华北水利水电大学 | 机器人用摆线针轮减速机的动态传动精度数值计算方法 |
CN115659549A (zh) * | 2022-11-14 | 2023-01-31 | 浙江大学 | 一种内曲线液压马达凸轮环曲线修形方法 |
CN115659549B (zh) * | 2022-11-14 | 2023-10-24 | 浙江大学 | 一种内曲线液压马达凸轮环曲线修形方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108253094B (zh) | 2019-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108253094B (zh) | 一种rv减速器摆线针轮承载啮合印痕确定方法及装置 | |
Sheveleva et al. | Algorithms for analysis of meshing and contact of spiral bevel gears | |
Li et al. | A manufacturing error measurement methodology for a rotary vector reducer cycloidal gear based on a gear measuring center | |
CN103105152B (zh) | 一种基于齿轮加工仿真模型的齿轮加工表面形貌分析方法 | |
CN107030698A (zh) | 机器人的逆运动学求解系统 | |
CN102636097B (zh) | 基于双面啮合的齿轮齿廓偏差测量方法 | |
Ding et al. | An innovative determination approach to tooth compliance for spiral bevel and hypoid gears by using double-curved shell model and Rayleigh–Ritz approach | |
Brauer | Transmission error in anti-backlash conical involute gear transmissions: a global–local FE approach | |
Qu et al. | An innovative semi-analytical determination approach to numerical loaded tooth contact analysis (NLTCA) for spiral bevel and hypoid gears | |
Iglesias et al. | Planet position errors in planetary transmission: Effect on load sharing and transmission error | |
CN106227940A (zh) | 一种摆线齿轮的建模方法 | |
CN108730480A (zh) | 一种rv减速器的摆线轮及其齿廓逆向主动修形方法 | |
Liu et al. | Investigation on contact and bending stress of face-hobbed and face-milled hypoid gear | |
Ding et al. | Automatic data-driven operation and optimization of uncertain misalignment by considering mechanical power transmission performances of spiral bevel and hypoid gears | |
CN106989157B (zh) | 基于运动学原理的渐开线圆柱齿轮双面啮合测量仿真方法 | |
Ding et al. | High-performance tooth flank collaborative optimization model for spiral bevel and hypoid gears | |
CN102645200A (zh) | 基于双面啮合的齿轮齿距偏差测量方法 | |
Han et al. | Topology modification method of generating gear grinding based on multi-axis linkage parameter optimization | |
CN111859576A (zh) | 机器人用rv减速器含间隙机构传动误差计算方法 | |
Yang et al. | Tooth surface contact analysis of involute rotate vector reducer based on a finite element linear programming method | |
CN106802989A (zh) | 一种考虑错位量影响的准双曲面齿轮接触计算方法 | |
Niu et al. | Geometrical design of variable ratio tooth profile based on Boolean subtraction operation and a novel modification method | |
Pacana | Development of Bevel Gear Motion Transmission Graphs with FEM | |
Li et al. | STATICS ANALYSIS OF CYLINDRICAL GEAR DRIVE AT ANY MESHING POSITION BASED ON ANSYS. | |
Li et al. | Study on transmission precision of 3k-type gear mechanism considering eccentricity error |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240423 Address after: No. 2 Xingye Third Road, Advanced Manufacturing Cluster Zone, Jianxi District, Luoyang City, Henan Province, 471003 Patentee after: Luoyang kedager Transmission Research Institute Co.,Ltd. Country or region after: China Address before: 471003 No. 48, Xiyuan Road, Jianxi District, Henan, Luoyang Patentee before: HENAN University OF SCIENCE AND TECHNOLOGY Country or region before: China |
|
TR01 | Transfer of patent right |