CN111059255B - 一种双圆弧谐波减速器齿面磨损的计算方法 - Google Patents
一种双圆弧谐波减速器齿面磨损的计算方法 Download PDFInfo
- Publication number
- CN111059255B CN111059255B CN201911371974.XA CN201911371974A CN111059255B CN 111059255 B CN111059255 B CN 111059255B CN 201911371974 A CN201911371974 A CN 201911371974A CN 111059255 B CN111059255 B CN 111059255B
- Authority
- CN
- China
- Prior art keywords
- wear
- tooth
- tooth profile
- gear
- arc
- 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
Images
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
- F16H57/00—General details of gearing
- F16H57/01—Monitoring wear or stress of gearing elements, e.g. for triggering maintenance
-
- 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
- F16H57/01—Monitoring wear or stress of gearing elements, e.g. for triggering maintenance
- F16H2057/012—Monitoring wear or stress of gearing elements, e.g. for triggering maintenance of gearings
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Gears, Cams (AREA)
Abstract
本发明公开了一种双圆弧谐波减速器齿面磨损的计算方法,首先,基于原始参数建立双圆弧谐波减速器柔轮—刚轮的齿廓方程;其次,确定齿形参数、接触正压力和磨损系数等参数;基于改进的Archard磨损模型,计算柔轮—刚轮齿面节点k在一个磨损周期内的磨损量;重复磨损循环直至该节点总的磨损量达到预设阈值时,齿面的接触正压力和齿形参数因齿面磨损而改变;进而,重新构建柔轮‑刚轮齿廓方程,重新计算接触正压力;将齿廓变化m次后节点k的磨损量求和,直至磨损量达到最大许用磨损量。本发明提出了双圆弧谐波减速器齿面磨损的计算方法,为提高谐波减速器传动精度和寿命预测的准确性提供了理论基础。
Description
技术领域
本发明属于齿轮磨损计算技术领域,具体涉及一种双圆弧谐波减速器齿轮磨损的计算方法。
背景技术
谐波减速器的结构包括波发生器、柔轮、刚轮。其中柔轮是柔性外齿轮,刚轮是刚性内齿圈。谐波减速器具有结构紧凑、体积小、质量轻、传动比与承载能力大、传动精度高等优点。相对于渐开线齿形的谐波减速器,双圆弧齿形的柔轮与刚轮在啮合过程中存在共轭运动及“双共轭”啮合区间,这有效减少了干涉,并对扭转刚度和传动精度的提高有重要作用。
国内外学者针对齿面磨损特性展开了一系列研究,多数研究者大多采用Archard磨损公式,认为磨损系数、接触载荷及滑动距离是影响齿面磨损的主要因素。Osman等建立了宽齿面和斜齿面的Archard磨损模型,预测齿面磨损沿齿宽及啮合方向的分布并研究了齿面磨损和啮合偏差对齿面载荷分布的影响;等基于Archard磨损公式建立了内啮合齿轮副的磨损模型,并通过数值仿真和试验的方式分析了沿啮合线方向齿面磨损的分布。张俊等基于Hertz理论和Archard磨损公式建立了渐开线直齿准静态磨损模型,并数值仿真了磨损的磨损量。
以上研究表明,现有的齿面磨损模型多用于渐开线齿轮,缺乏双圆弧齿轮磨损模型及机理研究。谐波减速器接触齿对较多,具有较大的啮合区,啮合齿对间的载荷分布规律难以确定,并且齿面磨损还会受到柔轮变形及啮合区润滑特性的影响,因此,其磨损机理本身具有复杂性。目前对柔轮—刚轮齿面磨损失效机理的研究非常有限,严重制约着谐波减速器运行可靠性的提高。
发明内容
本发明要解决的技术问题在于,针对现有技术的上述缺陷,提供一种双圆弧谐波减速器齿面磨损的计算方法。
本发明解决其技术问题所采用的技术方案是:构造一种双圆弧谐波减速器齿面磨损的计算方法,包括步骤:
(1)基于待检测的双圆弧谐波减速器,根据齿轮的基本参数构建出柔轮的齿廓方程,并得出共轭的刚轮齿廓方程;
(2)基于柔轮和刚轮的齿廓方程、接触模型、有限元模型和磨损实验对比方法分别得出接触正压力q、磨损系数k;
(3)基于Archard磨损公式,引入滑动系数推导出与接触正压力相关的磨损公式;
其中,基于齿面节点k在一个磨损周期内的磨损量,循环计算个磨损周期内齿面节点k的磨损量当磨损量达到预设阈值εq后,修正齿形参数和接触正压力q;将齿廓变化m次后节点k的磨损量求和,直至磨损量hr,g(k)达到最大许用磨损量εt,结束循环,以对双圆弧谐波减速器齿轮的寿命进预测。
其中,停止循环后,根据齿面各点的磨损量,构建出磨损后的柔轮-刚轮齿廓方程。
其中,谐波减速器采用特殊的双圆弧齿形,双圆弧齿形采用两段圆弧和一段直线构成,为了考虑柔轮弹性变形对共轭运动的影响,故用包络法求解刚轮的齿廓方程。
其中,在Archard磨损公式中引入滑动系数,磨损深度h为接触正压力q的函数,磨损公式为:
其中,k为磨损系数,q为接触正压力,bR为齿宽,H材料硬度,ε为滑动系数。
区别于现有技术,本发明的双圆弧谐波减速器齿面磨损的计算方法通过构建双圆弧谐波减速器柔轮和刚轮的齿廓方程,计算接触正压力和磨损系数,基于Archard磨损公式,引入滑动系数推导出与接触正压力相关的磨损公式,并计算齿面节点k在一个磨损周期内的磨损量,实现对双圆弧谐波减速器齿面磨损的计算。通过本发明,能够精确计算双圆弧谐波减速器的齿轮齿面磨损,避免柔轮变形及啮合区润滑特性的影响,提高谐波减速器运行的可靠性。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1为本发明提供的一种双圆弧谐波减速器齿面磨损的计算方法中双圆弧齿轮磨损计算流程图。
图2为本发明提供的一种双圆弧谐波减速器齿面磨损的计算方法中双圆弧齿廓示意图。
图3为本发明提供的一种双圆弧谐波减速器齿面磨损的计算方法中柔轮接触点及旋转角度示意图。
图4为本发明提供的一种双圆弧谐波减速器齿面磨损的计算方法中磨损前后齿廓线示意图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现对照附图详细说明本发明的具体实施方式。
如图1所示本发明的计算方法如下:
双圆弧谐波减速器柔轮、刚轮齿廓方程的建立:谐波减速器采用特殊的双圆弧齿形,双圆弧齿形采用两端圆弧和一段直线构成,如图2所示为双圆弧齿廓示意图。根据齿面的基本参数构建出柔轮的齿廓方程,为了考虑柔轮弹性变形对共轭运动的影响,故用包络法求解刚轮的齿廓方程;其中柔轮的齿廓方程中AB段圆弧方程可表示为:
其中ρa为柔轮齿廓AB段圆弧半径,α2为齿廓AB段圆弧内任一点对应的圆心角,la为柔轮AB段齿廓圆心偏移量,ea为AB段齿廓圆心移距量,hf为齿根高。t为柔轮壁厚;
柔轮的齿廓方程中BC段直线方程可表示为:
其中Sa为柔轮齿廓AB段节线齿厚,Wf为柔轮齿廓CD段节线齿槽宽,lf为柔轮齿廓CD段圆心偏移量,ρf为柔轮齿廓CD段圆弧半径。
其中b代表齿廓方程中BC段直线方程的截距,k代表齿廓方程中BC段直线方程的斜率。
柔轮的齿廓方程中CD段圆弧方程可表示为:
其中α1为齿廓CD段圆弧内任一点对应的圆心角,ef为CD段齿廓圆心移距量。
在保证传动比不变的前提下,以柔轮齿廓方程为基础,通过包络法求解柔轮和刚轮的共轭齿廓,从而得到对应的刚轮齿廓,柔轮坐标系变换到刚轮坐标系的坐标变换公式为:
从而得到刚轮齿廓的表达式:
刚轮齿廓方程中齿根圆弧段的共轭方程为:
刚轮齿廓方程中直线段的共轭方程为:
刚轮齿廓方程中齿顶圆弧段的共轭方程为:
基于齿形参数、接触模型、有限元模型和磨损实验对比分别得出柔轮圆弧半径(ρa、ρf)、接触正压力q、磨损系数k;
由于双圆弧特殊的齿形是由两端圆弧和一段直线段构成,所以其齿面接触正压力会随接触位置的不同而发生变化。如图3所示为柔轮受力分析图,齿面接触正压力q的方程可表示为:
其中,AB、CD是双圆弧齿面的圆弧段,BC代表直线段,T为双圆弧谐波减速器柔轮-刚轮承受的力矩,θ为不同接触点对应的角度,θbc为齿廓BC直线段的倾角。
磨损系数k通过有限元仿真结合磨损试验来确定。一定工况下,采用与柔轮—刚轮啮合副相同材质的摩擦副进行摩擦磨损试验,得到磨损量随磨损周期的实际变化曲线;
建立与摩擦磨损试验样件材料属性及工况相同的三维有限元模型,设定不同磨损系数,得到不同磨损系数下,磨损量随磨损周期的仿真变化曲线;对比仿真结果与试验结果,两条曲线基本吻合的条件下所对应的磨损系数值即为所需磨损系数k的预估值。
在Archard磨损公式中引入滑动系数,磨损深度h为接触正压力q的函数,磨损公式为:
其中,k为磨损系数,q为接触正压力,bR为齿宽,H材料硬度,ε为滑动系数。在此公式中可以看出,磨损深度h只与变量接触正压力q有关。
滑动系数的定义:
当齿面节点k的循环磨损量达到该节点预设阈值εq时,双圆弧齿面的齿形会因磨损发生变化,其接触正压力q也会随之改变。此时,需要修正参数再进行磨损循环过程。
将齿廓变化m次后节点k的磨损量求和:
齿面磨损量hr,g(k)达到最大许用磨损阈值εt后,停止循环。
停止循环后,根据齿面各点的磨损量,构建出磨损后的柔轮-刚轮齿廓方程。
随着磨损的进行,齿廓形状会发生改变,对齿廓线修正的方法如图4示,磨损后齿廓方程为:
其中未磨损齿面上节点k的坐标是xk,yk,沿齿面法向磨损量达到hk后,节点k的坐标变为x′k,y′k,上式中αk是反映斜率的角,以下就称之为斜率角。
利用修正后的节点坐标可以插值得到齿廓线的三次样条函数,进而利用样条函数计算新的斜率角αk,由几何关系可知αk与接触点对应的啮合角θ互余,可得出新的θ角度,用来计算新的接触正压力p和半径ρa,ρf,进行下一工作时间段的磨损计算。拟合出不同接触点经过磨损后的柔轮-刚轮齿面轮廓。
上面结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护之内。
Claims (4)
1.一种双圆弧谐波减速器齿面磨损的计算方法,其特征在于,包括以下步骤:
(1)基于待检测的双圆弧齿形谐波减速器,根据双圆弧齿形齿轮的基本参数构建出柔轮的齿廓方程,并得出共轭的刚轮齿廓方程;
双圆弧齿形齿轮包括AB、CD代表的双圆弧齿面的圆弧段,BC代表的直线段;其柔轮齿廓方程包括AB段圆弧方程、BC段直线方程和CD段圆弧方程;
柔轮的齿廓方程中AB段圆弧方程表示为:
其中ρa为柔轮齿廓AB段圆弧半径,α2为齿廓AB段圆弧内任一点对应的圆心角,la为柔轮AB段齿廓圆心偏移量,ea为AB段齿廓圆心移距量,hf为齿根高,t为柔轮壁厚;
柔轮的齿廓方程中BC段直线方程表示为:
其中Sa为柔轮齿廓AB段节线齿厚,Wf为柔轮齿廓CD段节线齿槽宽,lf为柔轮齿廓CD段圆心偏移量,ρf为柔轮齿廓CD段圆弧半径;
其中b代表齿廓方程中BC段直线方程的截距,k代表齿廓方程中BC段直线方程的斜率;
柔轮的齿廓方程中CD段圆弧方程表示为:
其中α1为齿廓CD段圆弧内任一点对应的圆心角,ef为CD段齿廓圆心移距量;
在保证传动比不变的前提下,以柔轮齿廓方程为基础,通过包络法求解柔轮和刚轮的共轭齿廓,从而得到对应的刚轮齿廓,柔轮坐标系变换到刚轮坐标系的坐标变换公式为:
从而得到刚轮齿廓的表达式:
(2)基于柔轮和刚轮的齿廓方程、接触模型、有限元模型和磨损实验对比方法分别得出接触正压力q、磨损系数k;
其中,接触正压力q公式表示为:
磨损系数k通过有限元仿真结合磨损试验来确定,采用与柔轮—刚轮啮合副相同材质的摩擦副进行摩擦磨损试验,得到磨损量随磨损周期的实际变化曲线;
建立与摩擦磨损试验样件材料属性及工况相同的三维有限元模型,设定不同磨损系数,得到不同磨损系数下,磨损量随磨损周期的仿真变化曲线;对比仿真结果与试验结果,两条曲线基本吻合的条件下所对应的磨损系数值即为所需磨损系数k的预估值;
(3)基于Archard磨损公式,引入滑动系数推导出与接触正压力相关的磨损公式;
磨损公式为:
其中,k为磨损系数,q为接触正压力,bR为齿宽,H为材料硬度,ε为滑动系数;
滑动系数的定义:
3.根据权利要求1所述的双圆弧谐波减速器齿面磨损的计算方法,其特征在于,停止循环后,根据齿面各点的磨损量,构建出磨损后的柔轮-刚轮齿廓方程。
4.根据权利要求1所述的一种双圆弧谐波减速器齿面磨损的计算方法,其特征在于,谐波减速器采用特殊的双圆弧齿形,双圆弧齿形采用两段圆弧和一段直线构成,为了考虑柔轮弹性变形对共轭运动的影响,故用包络法求解刚轮的齿廓方程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911371974.XA CN111059255B (zh) | 2019-12-27 | 2019-12-27 | 一种双圆弧谐波减速器齿面磨损的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911371974.XA CN111059255B (zh) | 2019-12-27 | 2019-12-27 | 一种双圆弧谐波减速器齿面磨损的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111059255A CN111059255A (zh) | 2020-04-24 |
CN111059255B true CN111059255B (zh) | 2021-08-20 |
Family
ID=70302707
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911371974.XA Active CN111059255B (zh) | 2019-12-27 | 2019-12-27 | 一种双圆弧谐波减速器齿面磨损的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111059255B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111810608B (zh) * | 2020-06-05 | 2023-05-09 | 金齿传动科技(大连)有限公司 | 具有连续啮合齿形的谐波齿轮的齿形设计方法 |
CN111737833B (zh) * | 2020-06-15 | 2024-05-07 | 大连理工大学 | 具有连续啮合齿形的谐波齿轮滑动系数的计算方法 |
CN112883519B (zh) * | 2021-03-19 | 2023-03-14 | 河南科技大学 | 一种摆线轮寿命预测方法及装置 |
CN113408075B (zh) * | 2021-07-05 | 2022-06-21 | 清华大学深圳国际研究生院 | 一种谐波减速器刚轮齿形优化方法 |
CN113722843B (zh) * | 2021-08-01 | 2023-04-04 | 北京工业大学 | 一种谐波减速器柔轮滚齿加工齿面残余高度计算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE69723291T2 (de) * | 1997-10-02 | 2004-04-22 | Harmonic Drive Systems Inc. | Elastisches zahnrad |
CN106845046A (zh) * | 2017-04-14 | 2017-06-13 | 东北大学 | 直齿圆柱齿轮磨损量计算方法 |
CN109871652A (zh) * | 2019-03-14 | 2019-06-11 | 东北大学 | 一种基于动态啮合力的齿轮副磨损量预测方法 |
CN110162924A (zh) * | 2019-06-03 | 2019-08-23 | 珠海格力电器股份有限公司 | 一种谐波传动机构以及谐波减速器 |
CN110188409A (zh) * | 2019-05-09 | 2019-08-30 | 东南大学 | 扭簧加载消隙齿轮齿面磨损量计算模型 |
-
2019
- 2019-12-27 CN CN201911371974.XA patent/CN111059255B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE69723291T2 (de) * | 1997-10-02 | 2004-04-22 | Harmonic Drive Systems Inc. | Elastisches zahnrad |
CN106845046A (zh) * | 2017-04-14 | 2017-06-13 | 东北大学 | 直齿圆柱齿轮磨损量计算方法 |
CN109871652A (zh) * | 2019-03-14 | 2019-06-11 | 东北大学 | 一种基于动态啮合力的齿轮副磨损量预测方法 |
CN110188409A (zh) * | 2019-05-09 | 2019-08-30 | 东南大学 | 扭簧加载消隙齿轮齿面磨损量计算模型 |
CN110162924A (zh) * | 2019-06-03 | 2019-08-23 | 珠海格力电器股份有限公司 | 一种谐波传动机构以及谐波减速器 |
Also Published As
Publication number | Publication date |
---|---|
CN111059255A (zh) | 2020-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111059255B (zh) | 一种双圆弧谐波减速器齿面磨损的计算方法 | |
CN109707822B (zh) | 一种应用于机器人、探测器的小模数谐波传动啮合齿形的设计方法 | |
Jyothirmai et al. | A finite element approach to bending, contact and fatigue stress distribution in helical gear systems | |
Brauer | Analytical geometry of straight conical involute gears | |
Danieli et al. | New developments in variable radius gears using constant pressure angle teeth | |
CN106845046B (zh) | 直齿圆柱齿轮磨损量计算方法 | |
Palmer et al. | Evaluation of methods for calculating effects of tip relief on transmission error, noise and stress in loaded spur gears | |
Kacalak et al. | Innovative design of non-backlash worm gear drives | |
Liu et al. | Investigation on contact and bending stress of face-hobbed and face-milled hypoid gear | |
CN111299982A (zh) | 一种农机齿轮修形加工工艺 | |
CN110374870B (zh) | 一种双螺杆压缩机螺杆转子间隙布置方法 | |
Inoue et al. | Derivation of path of contact and tooth flank modification by minimizing transmission error on face gear | |
CN115221655B (zh) | 一种斜齿非圆齿轮的传动精度评价方法 | |
Wink et al. | Investigation of tooth contact deviations from the plane of action and their effects on gear transmission error | |
Li et al. | A Loaded Analysis Method for RV Cycloidal-pin Transmission Based on the Minimum Energy Principle. | |
Liang et al. | Generation principle and meshing characteristics of conjugate-curve circular arc gears | |
Vishwakarma et al. | Finite element analysis of helical gear using three-dimensional cad model | |
Ambade et al. | Involute gear tooth contact and bending stress analysis | |
Liang et al. | Theoretical and experimental investigations on parallel-axis gear transmission with tubular meshing surfaces | |
Liang et al. | Calculation method of sliding ratios for conjugate-curve gear pair and its application | |
Pabiszczak et al. | Kinematic analysis and design of eccentric rolling transmission | |
Chen et al. | Contact stress of modified curvilinear gears | |
Yun et al. | Load distribution analysis for the plastic helical gear meshing with steel worm | |
Samya et al. | Analysis of profile shift factor's effect on bending stress of spur gears using the finite element method | |
CN109902437A (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 |