CN106989157B - 基于运动学原理的渐开线圆柱齿轮双面啮合测量仿真方法 - Google Patents
基于运动学原理的渐开线圆柱齿轮双面啮合测量仿真方法 Download PDFInfo
- Publication number
- CN106989157B CN106989157B CN201710180993.9A CN201710180993A CN106989157B CN 106989157 B CN106989157 B CN 106989157B CN 201710180993 A CN201710180993 A CN 201710180993A CN 106989157 B CN106989157 B CN 106989157B
- Authority
- CN
- China
- Prior art keywords
- gear
- tooth
- involute
- coordinate system
- equation
- 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
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
- F16H55/00—Elements with teeth or friction surfaces for conveying motion; Worms, pulleys or sheaves for gearing mechanisms
- F16H55/02—Toothed members; Worms
- F16H55/08—Profiling
- F16H55/0806—Involute profile
-
- 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
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Mechanical Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
基于运动学原理的渐开线圆柱齿轮双面啮合测量仿真方法,属于精密测试技术与仪器领域。该方法是通过建立齿轮模型,进而求解齿轮的双面啮合模型,得到在笛卡尔坐标系下的齿轮双面啮合的啮合线表达式,将齿轮的基圆半径误差、齿槽半角误差引入啮合线表达式中,得到齿轮在双面啮合过程中的径向、切向的变化量曲线。啮合线方程是在笛卡尔坐标系下求得,所以得到的是啮合线在径向、切向的参数方程。
Description
技术领域
本发明涉及精密测试技术与仪器领域,具体是一种基于渐开线齿轮双面啮合数学模型的齿轮双面啮合误差仿真分析方法。
背景技术
齿轮传动机构在航空航天、仪器仪表、汽车工业、精密机械等各个领域都有广泛的应用。齿轮精度、质量直接影响机器设备的效率、噪声、运动精度和使用寿命。因此,对齿轮的检验就显得尤为重要。
对齿轮进行检验的仪器有很多,包括齿轮单面啮合检查仪、齿轮双面啮合检查仪、齿轮测量中心等。其中对于中小模数齿轮,大多采用齿轮双面啮合检查仪,其具有对环境要求不高、结构简单、操作方便等优点。
通过理论分析对齿轮双面啮合过程进行仿真的方法中,传统的方法是通过齿轮建模,求解出齿轮双面啮合过程中,中心距变化与齿轮误差的关系,从而得出齿轮双面啮合的仿真误差曲线。该方法所得到的误差曲线为齿轮中心距变化曲线,即为齿轮径向变动的误差曲线,对于齿轮在双面啮合过程中的切向变动量并没有考虑。
发明内容
为解决背景技术中提出的问题,本发明提出了一种新的齿轮双面啮合误差分析方法。该方法是通过建立齿轮模型,进而求解齿轮的双面啮合模型,得到在笛卡尔坐标系下的齿轮双面啮合的啮合线表达式,将齿轮的基圆半径误差、齿槽半角误差引入啮合线表达式中,得到齿轮在双面啮合过程中的径向、切向的变化量曲线。
本发明采用的技术方案为基于运动学原理的渐开线圆柱齿轮双面啮合测量仿真方法,该方法的实现过程如下:
1)在笛卡尔坐标系下,写出齿轮一个齿面的渐开线表达式,再通过坐标变化,得到齿轮所有渐开线齿面的表达式,进而得到齿轮的数学模型;
2)在建立齿轮数学模型的基础上,根据齿轮啮合原理——齿轮啮合过程中,在相互包络齿廓的接触点处,相对运动的速度矢量应当垂直于齿廓的法线矢量,求解出齿轮双面啮合方程;
3)将齿轮的渐开线表达式与齿轮双面啮合方程结合,得到双面啮合的啮合线方程;
4)由于该啮合线方程是在笛卡尔坐标系下求得,所以得到的是啮合线在径向、切向的参数方程。该参数方程的变量中含有齿轮的基本参数,由误差理论推导出啮合线在径向、切向的变化量与齿轮参数变化量之间的关系,从而得到齿轮双面啮合的径向、切向误差变化曲线。
5)有误差的齿轮双面啮合过程主要分为两个部分——渐开线啮合部分和顶刃啮合部分,最终的双啮误差分析仿真为两部分的结合。
由于采用了上述推导过程,本发明的有益效果是:
根据齿轮基本参数的变化量,带入到推导出的齿轮双面啮合误差曲线函数中,可以得到齿轮双面啮合径向及切向变化曲线。
附图说明
图1齿轮的数学模型
图2齿轮的双面啮合数学模型
图3一个周节角内基圆半径变化-10μm对径向偏差的影响
图4一个周节角内基圆半径变化-10μm对切向偏差的影响
具体实施方式
为使本方法的实现过程更加清晰,以下结合附图及分析公式对本方法做进一步说明。
在笛卡尔坐标系下,求解出齿轮1一个渐开线齿面的表达式:
为展开角,θ为渐开线角,α为压力角,ξ为齿槽半角,rb1为齿轮1的基圆半径,xR1、yR1分别表示齿轮1第一个渐开线右齿面上任意一点在笛卡尔坐标系下的横坐标、纵坐标。
上述表达式表示第一个齿的右齿面方程,当齿轮1有z个齿时,则其第i个齿的右齿面方程可通过该表达式乘以旋转矩阵:
z表示齿轮1的的齿数
即齿轮1第i个齿的右齿面的渐开线方程为:
xRi、yRi分别表示齿轮1第i个渐开线右齿面上任意一点在笛卡尔坐标系下的横坐标、纵坐标。
对于齿轮1左齿面而言,只需求取齿轮1右齿面第一段渐开线关于x轴对称的曲线,即可得到第z个齿的左齿面渐开线方程:
xLz、yLz分别表示齿轮1第z个渐开线右齿面上任意一点在笛卡尔坐标系下的横坐标、纵坐标。
上式表示齿轮1第z个齿的左齿面,要将其转化为齿轮1第i个齿的左齿面只需乘以旋转矩阵:
即齿轮1第i个齿的左齿面的渐开线方程为:
xLi、yLi分别表示齿轮1第i个渐开线右齿面上任意一点在笛卡尔坐标系下的横坐标、纵坐标。
将左右齿面渐开线模型联立,即可得到齿轮渐开线数学模型,参照图1。
在所建立的齿轮渐开线数学模型基础上,建立齿轮双面啮合数学模型,改模型的建立过程如下:
如图2中所示齿轮1和齿轮2,以齿轮中心为坐标原点,建立坐标系。其中齿轮1所在的坐标系C1[O1;x1,y1,z1]为被测齿轮的可动坐标系,实线坐标系为坐标系C1初始位置,虚线坐标系为坐标系C1随齿轮1旋转后的位置;齿轮2所在的坐标C2[O2;x2,y2,z2]为标准齿轮的可动坐标系,实线坐标系为坐标系C1初始位置,虚线坐标系为坐标系C1随齿轮1旋转后的位置;两齿轮分度圆交点为固定坐标系的原点O,其参考坐标系为C[O;x,y,z]。
结合齿轮啮合原理——齿轮啮合过程中,在相互包络齿廓的接触点处,相对运动的速度矢量应当垂直于齿廓的法线矢量,对双面啮合数学模型进行求解,解出齿轮双面啮合方程:
齿轮1第i个齿的右齿面啮合方程:
表示齿轮1旋转的角度,a表示齿轮1与齿轮2的理论中心距,i21表示齿轮2与齿轮1的转动比。
齿轮1第i个齿的左齿面啮合方程:
将上述的啮合线方程与推导出的齿轮渐开线方程联立,及可得到齿轮双面啮合的啮合线表达式。啮合线上任意一点的坐标(x,y)与rb1,ξ四个量有关,即:
当这四个量存在偏差时,理论啮合线与实际啮合线的位置就会发生变化,从而引起测量上的误差。考虑四个量中有两个量存在偏差,分别为Δrb1,Δξ,理论的啮合点为Q(x,y),变化之后的啮合点为Q′,则Q′的坐标为(x′,y′),则:
x′=x+Δx
y′=y+Δy
由误差理论可知:
将参数带入上述表达式,可得:
左齿面:
右齿面:
当被测齿轮存在基圆半径偏差时,会出现顶刃啮合的现象。在顶刃啮合期间,被测齿轮的齿顶与标准齿轮的齿面相接触,不符合啮合原理。此时,对于基圆半径偏差引起的啮合线的变化,需另作计算。
具有基圆半径误差为Δrb1的被测齿轮在啮合时,由顶刃啮合引起的啮合线增量可按下式计算:
ΔL表示啮合线方向增量,表示齿轮1从顶刃啮合开始到结束齿轮转过的角度,Δx表示径向变化量,Δy表示切向变化量。
通过Matlab将渐开线啮合与顶刃啮合的结果进行处理,即可得到基圆半径存在定值误差时,齿轮双面啮合在径向、切向的误差变化曲线。
基于上述误差分析,本发明给出了如下的一组齿轮数据进行仿真实验:
名称 | 被测齿轮(齿轮1) | 标准齿轮(齿轮2) |
模数m | 2.5 | 2.5 |
齿数z | 45 | 70 |
齿顶高系数 | 1 | 1 |
齿根高系数 | 1.25 | 1.25 |
压力角α | 20° | 20° |
Claims (2)
1.基于运动学原理的渐开线圆柱齿轮双面啮合测量仿真方法,其特征在于:该方法的实现过程如下,
1)在笛卡尔坐标系下,写出齿轮一个齿面的渐开线表达式,再通过坐标变化,得到齿轮所有渐开线齿面的表达式,进而得到齿轮的数学模型;
2)在建立齿轮数学模型的基础上,根据齿轮啮合原理——齿轮啮合过程中,在相互包络齿廓的接触点处,相对运动的速度矢量应当垂直于齿廓的法线矢量,求解出齿轮双面啮合方程;
3)将齿轮的渐开线表达式与齿轮双面啮合方程结合,得到双面啮合的啮合线方程;
4)由于该啮合线方程是在笛卡尔坐标系下求得,所以得到的是啮合线在径向、切向的参数方程;该参数方程的变量中含有齿轮的基本参数,由误差理论推导出啮合线在径向、切向的变化量与齿轮参数变化量之间的关系,从而得到齿轮双面啮合的径向、切向误差变化曲线;
5)有误差的齿轮双面啮合过程主要分为两个部分——渐开线啮合部分和顶刃啮合部分,最终的双啮误差分析仿真为两部分的结合。
2.根据权利要求1所述的基于运动学原理的渐开线圆柱齿轮双面啮合测量仿真方法,其特征在于:
在笛卡尔坐标系下,求解出齿轮1一个渐开线齿面的表达式:
为展开角,θ为渐开线角,α为压力角,ξ为齿槽半角,rb1为齿轮1的基圆半径,xR1、yR1分别表示齿轮1第一个渐开线右齿面上任意一点在笛卡尔坐标系下的横坐标、纵坐标;
上述表达式表示第一个齿的右齿面方程,当齿轮1有z个齿时,则其第i个齿的右齿面方程可通过该表达式乘以旋转矩阵:
z表示齿轮1的的齿数;
即齿轮1第i个齿的右齿面的渐开线方程为:
xRi、yRi分别表示齿轮1第i个渐开线右齿面上任意一点在笛卡尔坐标系下的横坐标、纵坐标;
对于齿轮1左齿面而言,只需求取齿轮1右齿面第一段渐开线关于x轴对称的曲线,即可得到第z个齿的左齿面渐开线方程:
xLz、yLz分别表示齿轮1第z个渐开线右齿面上任意一点在笛卡尔坐标系下的横坐标、纵坐标;
上式表示齿轮1第z个齿的左齿面,要将其转化为齿轮1第i个齿的左齿面只需乘以旋转矩阵:
即齿轮1第i个齿的左齿面的渐开线方程为:
xLi、yLi分别表示齿轮1第i个渐开线右齿面上任意一点在笛卡尔坐标系下的横坐标、纵坐标;
将左右齿面渐开线模型联立,即可得到齿轮渐开线数学模型;
在所建立的齿轮渐开线数学模型基础上,建立齿轮双面啮合数学模型,该模型的建立过程如下:
齿轮1和齿轮2,以齿轮中心为坐标原点,建立坐标系;其中齿轮1所在的坐标系C1[O1;x1,y1,z1]为被测齿轮的可动坐标系,实线坐标系为坐标系C1初始位置,虚线坐标系为坐标系C1随齿轮1旋转后的位置;齿轮2所在的坐标C2[O2;x2,y2,z2]为标准齿轮的可动坐标系,实线坐标系为坐标系C1初始位置,虚线坐标系为坐标系C1随齿轮1旋转后的位置;两齿轮分度圆交点为固定坐标系的原点O,其参考坐标系为C[O;x,y,z];
结合齿轮啮合原理——齿轮啮合过程中,在相互包络齿廓的接触点处,相对运动的速度矢量应当垂直于齿廓的法线矢量,对双面啮合数学模型进行求解,解出齿轮双面啮合方程:
齿轮1第i个齿的右齿面啮合方程:
表示齿轮1旋转的角度,a表示齿轮1与齿轮2的理论中心距,i21表示齿轮2与齿轮1的转动比;
齿轮1第i个齿的左齿面啮合方程:
将左齿面啮合方程与推导出的左齿面的渐开线方程联立,将右齿面啮合方程与推导出的右齿面的渐开线方程联立即可得到齿轮双面啮合的啮合线表达式;啮合线上任意一点的坐标(x,y)与rb1,ξ四个量有关,即:
当这四个量存在偏差时,理论啮合线与实际啮合线的位置就会发生变化,从而引起测量上的误差;考虑四个量中有两个量存在偏差,分别为Δrb1,Δξ,理论的啮合点为Q(x,y),变化之后的啮合点为Q′,则Q′的坐标为(x′,y′),则:
x′=x+Δx
y′=y+Δy
由误差理论可知:
将参数带入上述表达式,可得:
左齿面:
右齿面:
当被测齿轮存在基圆半径偏差时,会出现顶刃啮合的现象;在顶刃啮合期间,被测齿轮的齿顶与标准齿轮的齿面相接触,不符合啮合原理;此时,对于基圆半径偏差引起的啮合线的变化,需另作计算;
具有基圆半径误差为Δrb1的被测齿轮在啮合时,由顶刃啮合引起的啮合线增量可按下式计算:
ΔL表示啮合线方向增量,表示齿轮1从顶刃啮合开始到结束齿轮转过的角度,Δx表示径向变化量,Δy表示切向变化量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710180993.9A CN106989157B (zh) | 2017-03-24 | 2017-03-24 | 基于运动学原理的渐开线圆柱齿轮双面啮合测量仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710180993.9A CN106989157B (zh) | 2017-03-24 | 2017-03-24 | 基于运动学原理的渐开线圆柱齿轮双面啮合测量仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106989157A CN106989157A (zh) | 2017-07-28 |
CN106989157B true CN106989157B (zh) | 2019-06-14 |
Family
ID=59413362
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710180993.9A Active CN106989157B (zh) | 2017-03-24 | 2017-03-24 | 基于运动学原理的渐开线圆柱齿轮双面啮合测量仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106989157B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108256174B (zh) * | 2017-12-27 | 2021-05-11 | 西安建筑科技大学 | 用于分析剃削过程对齿形质量影响的剃齿啮合建模方法 |
CN109751386B (zh) * | 2019-02-14 | 2021-09-21 | 重庆模源齿轮有限公司 | 一种反渐开线齿轮传动啮合的设计计算方法 |
CN110657986B (zh) * | 2019-10-11 | 2021-07-02 | 北京工业大学 | 齿轮双面啮合测量中测量力引入轮齿变形计算方法 |
CN111950178B (zh) * | 2020-07-22 | 2022-08-16 | 中国第一汽车股份有限公司 | 一种基于HyperWorks软件的齿轮自动化加载方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5662008A (en) * | 1993-08-30 | 1997-09-02 | Teijin Seiki Boston, Inc. | Extended contact harmonic drive devices |
SI22423B (sl) * | 2006-12-15 | 2011-10-28 | Gorazd Hlebanja | Parabolična valjasta polžasta zobniška dvojica |
CN100567935C (zh) * | 2008-06-27 | 2009-12-09 | 北京工业大学 | 一种齿轮整体误差测量装置及方法 |
CN104565285A (zh) * | 2013-10-14 | 2015-04-29 | 空中客车营运有限公司 | 齿轮构造方法和数字设备 |
CN105240486B (zh) * | 2015-09-24 | 2017-11-10 | 无锡职业技术学院 | 一种弧齿锥齿轮的滑动系数及重合度的计算方法 |
-
2017
- 2017-03-24 CN CN201710180993.9A patent/CN106989157B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106989157A (zh) | 2017-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106989157B (zh) | 基于运动学原理的渐开线圆柱齿轮双面啮合测量仿真方法 | |
CN102636097B (zh) | 基于双面啮合的齿轮齿廓偏差测量方法 | |
CN108645322B (zh) | 安装误差影响下的渐开线圆柱直齿轮齿距偏差评定方法 | |
CN101915667B (zh) | 齿轮副整体误差测量技术及其测量方法 | |
CN112257192B (zh) | 一种基于实测齿廓偏差的一齿切向综合偏差计算方法 | |
CN102706306B (zh) | 渐开线齿轮波度样板 | |
CN108645323A (zh) | 安装误差影响下的渐开线圆柱直齿轮齿廓偏差评定方法 | |
CN102645200B (zh) | 基于双面啮合的齿轮齿距偏差测量方法 | |
CN108730480A (zh) | 一种rv减速器的摆线轮及其齿廓逆向主动修形方法 | |
CN104990707B (zh) | 一种圆柱齿轮参数精确反求方法 | |
CN103162959A (zh) | 多功能齿轮样板 | |
CN105300691A (zh) | 基于最佳安装距的锥齿轮传动误差测量方法 | |
CN106980724A (zh) | 一种基于复杂齿面网格规划原则的齿面精确模型截面放样构建方法 | |
CN108253094A (zh) | 一种rv减速器摆线针轮承载啮合印痕确定方法及装置 | |
CN104728390A (zh) | 基于圆锥面渐开线的圆锥齿轮齿形设计方法 | |
CN111666643A (zh) | 一种用于确定复杂齿面接触性能的方法 | |
Yu et al. | A virtual measurement method of the transmission error based on point clouds of the gear | |
CN106326575B (zh) | 一种基于坐标法的非圆齿轮检测方法 | |
CN115221655B (zh) | 一种斜齿非圆齿轮的传动精度评价方法 | |
CN105241415A (zh) | 一种接触式渐开线蜗杆齿形的测量方法 | |
CN203365034U (zh) | 一种带齿轮转盘轴承装配及摩擦力矩测量机 | |
CN204807010U (zh) | 一种螺旋锥齿轮齿面三坐标测量的辅助量具 | |
CN111859576A (zh) | 机器人用rv减速器含间隙机构传动误差计算方法 | |
CN109489599B (zh) | 一种自校准的渐开线齿轮副振动位移及动态传动误差测量方法 | |
Wang et al. | Positioning accuracy prediction and tolerance allocation of cycloid reducers based on a BDTCA model |
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 |