CN110657986B - 齿轮双面啮合测量中测量力引入轮齿变形计算方法 - Google Patents

齿轮双面啮合测量中测量力引入轮齿变形计算方法 Download PDF

Info

Publication number
CN110657986B
CN110657986B CN201910965151.3A CN201910965151A CN110657986B CN 110657986 B CN110657986 B CN 110657986B CN 201910965151 A CN201910965151 A CN 201910965151A CN 110657986 B CN110657986 B CN 110657986B
Authority
CN
China
Prior art keywords
gear
meshing
tooth
stiffness
deformation
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
Application number
CN201910965151.3A
Other languages
English (en)
Other versions
CN110657986A (zh
Inventor
汤洁
杨炳耀
石照耀
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201910965151.3A priority Critical patent/CN110657986B/zh
Publication of CN110657986A publication Critical patent/CN110657986A/zh
Application granted granted Critical
Publication of CN110657986B publication Critical patent/CN110657986B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/02Gearings; Transmission mechanisms
    • G01M13/021Gearings
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M5/00Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
    • G01M5/0041Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining deflection or stress

Abstract

本发明公开了一种齿轮双面啮合测量中测量力引入轮齿变形计算方法,属于精密测试技术与仪器领域。包括步骤:(1)双面啮合时测量力会使轮齿发生变形,从而使双啮中心距发生变化,进而影响双啮误差曲线。计算沿左右啮合线方向的由测量力引起的轮齿变形,将左右啮合线上的轮齿变形量按啮合相位叠加,从而得到测量力引起的双啮中心距变化量。(2)沿啮合线方向的变形量为测量力沿啮合线方向的分量除以啮合线上的齿轮刚度。(3)单条啮合线上的轮齿刚度为该啮合线上各个接触点刚度的迭加。(4)齿轮在一个接触点处的刚度由三部分组成:轮齿的弯曲刚度、剪切刚度和齿面的接触刚度。(5)弯曲和剪切刚度由悬臂梁模型计算,接触变形由赫兹公式计算。

Description

齿轮双面啮合测量中测量力引入轮齿变形计算方法
技术领域
本发明涉及精密测试技术与仪器领域,具体是一种双面啮合测量中测量力影响下的轮齿变形计算方法。
背景技术
齿轮作为核心机械基础件,其质量和性能决定了装备的性能和可靠性。齿轮产业已经成为机械通用基础件中规模最大的产业。齿轮广泛应用于汽车制造、航空航天、风电设备以及装备制造等领域。齿轮精度、质量直接影响装备的效率、噪声、运动精度和使用寿命。由于齿轮测量技术与齿轮加工技术相辅相成,因此,对齿轮的测量就显得尤为重要。
双面啮合测量以其反映信息好、仪器结构简单、操作维护方便、检验效率高、测量稳定、易于实现机械化和自动化、适合生产现场使用等优点,成为了现代齿轮生产中最终检验的主要测试方法之一。在齿轮工程领域中对齿轮的质量和数量的要求不断提高,齿轮生产现场对齿轮快速测量需求与日俱增,因此齿轮双啮综合测量是齿轮测量技术的发展的主要方向之一。
齿轮双面啮合测量仪通常由弹簧提供测量力,使产品齿轮和测量齿轮作双面啮合。对于小模数齿轮,轮齿厚度较小,在测量力作用下容易产生变形。尤其对于塑料齿轮,测量力对双面啮合测量结果影响较大。齿轮在测量力的作用下产生的变形所引入的误差将会对双面啮合测量仪的测量结果产生一定影响。
目前双面啮合测量仪的测量力多为根据经验给出的推荐值,而在现场的齿轮测量过程中,在不同测量力下,仪器往往会给出不同的测量结果,甚至会给出不同的精度评定等级。但是一般情况下仪器对于不同齿轮并没有给出不同的测量力标准,也没有对测量力引起的误差进行补偿。因此,分析测量力作用下的轮齿变形,以及轮齿变形对双面啮合测量的影响是双面啮合测量的一个重要问题。
本发明将围绕以下三个方面研究双面啮合下的轮齿变形:1) 双面啮合测量下的测量力模型;2)单对齿啮合刚度计算模型;3)双面啮合综合啮合刚度计算模型。
本发明通过推导单对齿啮合时的刚度来计算齿轮双面啮合时的齿轮综合刚度,将齿轮综合刚度代入双面啮合测量力模型来计算双面啮合测量时测量力引起的轮齿变形量。
发明内容
为解决背景技术中提出的问题,本发明提出了一种双面啮合测量时测量力引起的轮齿变形量计算方法。该方法通过推导单对齿啮合时的刚度来计算齿轮双面啮合时的齿轮综合刚度,将齿轮综合刚度代入双面啮合测量力模型来计算双面啮合测量时测量力引起的轮齿变形量。
本发明是采用如下理论推导实现的:
双面啮合测量齿轮与产品齿轮做紧密的无侧隙啮合需要给一定的径向压力以保证齿轮的紧密啮合。径向的压力会使轮齿发生变形,从而使双啮中心距发生变化,进而影响双啮误差曲线。计算沿左右啮合线方向的由测量力引起的轮齿变形,将左右啮合线上的轮齿变形量按啮合相位叠加,从而得到测量力引起的双啮中心距变化量。
双面啮合时轮齿沿径向的变形由轮齿沿左右啮合线方向的变形叠加得到。通过计算沿啮合线方向的由测量力引起的轮齿变形得到测量力引起的双啮中心距变化量。单条啮合线上的轮齿变形量等于弹簧力沿啮合线方向的分量除以轮齿在这条啮合线上的综合刚度。单条啮合线上的轮齿刚度为该啮合线上各个接触点刚度的叠加。
我们通过计算轮齿在给定力的作用下沿啮合线方向的变形量来求得一对齿在接触点处的刚度。轮齿沿啮合线方向的变形有三部分组成,轮齿的弯曲变形、剪切变形和齿面的接触变形。
把轮齿看做非均匀的悬臂梁来计算轮齿的变形。轮齿有效接触部分沿对称轴线被分成一系列长方形微元,每个微元由符号i表示,将每一个长方形微元看做一个均匀的悬臂梁,所以整个轮齿的变形量就是单个微元变形的叠加。使用悬臂梁模型可以计算轮齿的弯曲变形和剪切变形。
一对外齿轮的赫兹接触刚度沿着整个作用线线性化为一个常数,与接触位置无关,轮齿的接触变形可由赫兹公式计算。
齿轮在双面啮合期间在一条啮合线上有多对齿同时参与啮合,一对齿轮在啮合时同时有多对轮齿参与啮合,齿轮的啮合刚度是多对齿同时啮合刚度的叠加。齿轮啮合时单齿啮合时期和双齿啮合时期交替出现,在不同啮合时期轮齿的啮合点数不同,各个啮合点所处的位置也不同。齿轮的综合啮合刚度由双面啮合过程和单齿啮合刚度来确定。
啮合循环周期τm包含双齿啮合周期τd和单齿啮合周期τs,由于齿轮存在误差,啮合过程中还包含顶刃啮合阶段。
在各个啮合阶段将各接触点处的赫兹接触刚度、弯曲刚度和剪切刚度叠加得到轮齿综合刚度,轮齿综合刚度是齿轮转角的函数。因此可以得到轮齿变形量关于齿轮转角的函数。将左右啮合线方向的轮齿变形按啮合相位叠加再乘以换算系数1/2sinα,即得到双面啮合中中心距方向的轮齿变形量。
本发明的优势在于:
计算双面啮合测量中测量力引起的轮齿变形,即测量力引起的中心距变化,可以了解测量力对双面啮合测量结果的影响。
附图说明
图1测量力引起的齿轮变形误差计算模型
图2齿轮双面啮合沿啮合线方向变形
图3轮齿的悬臂梁模型
图4渐开线的极坐标和直角坐标
图5齿轮转角
图6轮齿接触变形
图7齿轮啮合周期
图8产品齿轮带有正基节偏差时的双面啮合过程
图9产品齿轮左右啮合线格子图
图10为单个轮齿的变形量关于齿轮转角的函数
图11为一对轮齿啮合的变形量关于齿轮转角的函数
图12为单个轮齿的刚度变化关于齿轮转角的函数
图13为一对轮齿啮合刚度关于齿轮转角的函数
图14为双面啮合时单侧齿面的啮合刚度
图15为双面啮合时单侧齿面的变形量
图16为双面啮合时齿轮在径向力作用下沿径向的总变形量关于齿轮转角的函数
具体实施方式
以下结合附图及推导公式对本发明做进一步说明。
如图1所示,F为作用在产品齿轮上的径向压力,在双面啮合中此压力由弹簧力提供,δL和δR分别为沿左右两条啮合线方向上的变形量,δ为齿轮在径向的变形量,δ与δL和δR的关系如下:
Figure GDA0003060890810000051
α’为啮合角,δL、δR为沿左右两条啮合线方向上的变形量。
轮齿沿啮合线方向的变形量与轮齿在啮合线方向上的受力和轮齿综合刚度有关,即:
Figure GDA0003060890810000052
Figure GDA0003060890810000061
FL、FR为径向力沿左右啮合线方向的分量,kL、kR为沿啮合线方向的综合刚度。
如图2所示,FL、FR为径向力沿左右啮合线方向的分量,其大小为:
Figure GDA0003060890810000062
F为双面啮合时弹簧力所提供的径向力,FL、FR为径向力沿左右啮合线方向的分量,α’为啮合角。
齿轮转到不同位置时一条啮合线上的接触点个数不同,齿轮在左右两条啮合线上的刚度kL、kR分别为左右啮合线上各个接触点处的刚度的叠加。在图2中,N1N2和N3N4分别为左右两条啮合线。当齿轮转至图中所示位置时,轮齿在左啮合线N1N2上的刚度kL是N1N2 上两个接触点b和d两点的刚度叠加,沿右啮合线N3N4方向的刚度 kR是N3N4上两个接触点a和c两点的刚度叠加。
Figure GDA0003060890810000063
kR为齿轮沿右啮合线方向的综合刚度,ki为右啮合线上各个接触点的刚度,j为右啮合线上接触点的个数。
将用解析法推导齿轮啮合刚度方程,当齿轮转到不同位置时齿轮的接触点不同,所以在不同位置轮齿的刚度也不同,因此目标是将啮合刚度方程表示为齿轮旋转角度的函数。首先需要推导在一对齿上的一个接触点处的啮合刚度方程,然后将一条啮合线上的各个接触点的刚度叠加从而得到轮齿在一条啮合线上的啮合刚度。
在一个接触点处轮齿的刚度为啮合的两个齿轮的刚度的叠加,因此一对轮齿在一个接触点处的刚度ki为:
Figure GDA0003060890810000071
kb1、kb2为分别表示测量齿轮和产品齿轮的轮齿弯曲刚度,ks1、 ks2分别表示测量齿轮和产品齿轮的轮齿剪切刚度,kh为轮齿齿面的接触刚度。
在一个齿轮上,轮齿在受力点处的刚度由三部分组成:轮齿的弯曲刚度、剪切刚度和赫兹接触刚度。三种刚度的计算如下:
Figure GDA0003060890810000072
Figure GDA0003060890810000073
Figure GDA0003060890810000074
δs为轮齿的弯曲变形量,δs为轮齿的剪切变形量,δh为轮齿齿面的接触变形量,Fi为作用在轮齿上的给定作用力。
通过材料力学的方法计算轮齿的弯曲和剪切变形,把轮齿看做非均匀的悬臂梁来计算轮齿的变形,如图3所示。轮齿有效接触部分沿对称轴线被分成一系列长方形微元,若将轮齿分为n个为微元,每个微元由符号i表示,将每一个长方形微元看做一个均匀的悬臂梁,所以整个轮齿的变形量就是n个微元变形的叠加。使用悬臂梁模型可以计算轮齿的弯曲变形和剪切变形。设微元宽度为σi,横截面积为Ai,惯性矩为Ii,li表示微元i与载荷作用点之间沿X轴方向的距离,l0 为基圆与载荷作用点之间沿X轴方向的距离,αx为啮合角。
轮齿的弯曲变形量δb和剪切变形量δs由下式表示
Figure GDA0003060890810000081
Figure GDA0003060890810000082
其中n为轮齿分成的微元数,δb为弯曲变形,δs为剪切变形, Fi为沿啮合线方向的给定作用力,σi为微元部分宽度,
Figure GDA0003060890810000086
为微元两个侧面截面模量的均值,G为轮齿的剪切模量,
Figure GDA0003060890810000083
为微元两个侧面截面积的均值,li表示微元i与载荷作用点之间沿X轴方向的距离,μ为泊松比,E为弹性模量,α'x为载荷方向与垂直方向的夹角。
上式中,
Figure GDA0003060890810000084
G=E/2(1+μ)。
Ii为微元处的截面模量,齿轮在微元处的截面为矩形,因此截面模量为:
Figure GDA0003060890810000085
b为齿宽,hi为微元处的齿厚,Ai为微元处的截面积,可由下式计算:
Ai=bhi (13)
式(10)α'x为载荷方向与y轴的夹角,可由下式计算:
α'x=tanαxB (14)
其中α'x为载荷方向与y轴的夹角,αx为啮合角,εB为齿廓与基圆交点和圆心连线与X轴的夹角。
将α'x代入式(10)可得轮齿的弯曲变形δb与剪切变形δs关于受力点处啮合角αx的函数,即δbx)和δsx)。齿轮转角
Figure GDA0003060890810000091
因此可以得到轮齿变形量关于齿轮
Figure GDA0003060890810000092
的函数,即
Figure GDA0003060890810000093
Figure GDA0003060890810000094
轮齿在啮合中,受力点从齿顶移动到齿根,如图5所示,受力点B2随齿轮转动移动到B1,αx从0变到αa,受力点从齿顶移动到齿根的过程中,齿轮转过的角度为
Figure GDA0003060890810000095
Figure GDA0003060890810000096
Figure GDA0003060890810000097
带入式(7)(8)中,可以得到轮齿的刚度也是齿轮转角
Figure GDA0003060890810000098
的函数,即
Figure GDA0003060890810000099
Figure GDA00030608908100000910
轮齿的弯曲刚度kb为:
Figure GDA00030608908100000911
轮齿的剪切刚度为ks为:
Figure GDA00030608908100000912
一对外齿轮的赫兹接触刚度kh沿着整个作用线线性化为一个常数,与接触位置无关。齿面的接触变形如图6所示,一对轮齿在接触点处由于弹簧力的作用而产生接触变形,变形量可由赫兹公式进行计算,齿面接触变形为:
Figure GDA00030608908100000913
μ1为测量齿轮的泊松比,μ2为产品齿轮的泊松比,E1为测量齿轮的弹性模量,E2为产品齿轮的弹性模量。轮齿的接触刚度kh为:
Figure GDA0003060890810000101
其中kh轮齿的接触刚度,δh为轮齿的接触变形,Fi为沿啮合线方向的给定作用力。
通常齿轮在双面啮合期间在一条啮合线上有多对齿同时参与啮合,因此整个齿轮的啮合刚度是多对齿同时啮合刚度的合成。在不同啮合时期轮齿的啮合点数不同,各个啮合点所处的位置也不同。同时参与啮合的齿对数与一对齿轮的重合度有关,齿轮的重合度要大于 1才能保证连续的传动,一般齿轮的重合度大于1小于2。这里仅讨论重合度大于1且小于2的情况下,齿轮综合时变啮合刚度的合成。
图7表示一个齿轮轮齿的啮合循环周期,啮合循环周期τm包含双齿啮合周期τd和单齿啮合周期τs,可定义比例系数为:
λ=τdm (19)
比例系数λ与齿轮对的重合度ε相关,关系如下
λ=ε-1 (20)
有:0≤t≤λτm时,即从初始接触点到单齿接触最低点,为双齿啮合区;λτm≤t≤τm时,即从单齿接触最低点到单齿接触最高点,为单齿接触区;τm≤t≤(1+λ)τm时,即从单齿啮合最高点到最终接触点,为双齿啮合区。至此,一个轮齿完成从进入啮合到脱出啮合的啮合周期。
对于双面啮合的啮合过程如图8所示,齿轮1为测量齿轮,齿轮2为产品齿轮,齿轮2为从动轮,其在弹簧力F的作用下与齿轮1 做双面啮合。假设产品齿轮即齿轮2,左右齿面带有正基节偏差+Δfpb,它与测量齿轮1的啮合过程如图8所示。
图8a中,齿轮2的齿2右齿面与齿轮1的齿3提前于啮合线外开始接触而作顶刃啮合,顶刃啮合阶段齿轮2的齿2的齿顶在齿轮 1的齿3齿面上进行刮行而非渐开线啮合,此过程也称作齿顶刮行或右顶刮。如图b中所示,直到两者的接触点转至啮合线上后两齿才开始渐开线啮合,此时,齿轮2的左右齿面均为双齿啮合区,在左右啮合线上各有两个接触点,齿轮对共有四个接触点。如图c所示,齿轮 2的齿1左齿面齿顶转至啮合线上,由于左齿面存在正的基节偏差 +Δfpb,齿顶在并未退出啮合,而是在齿面上进行刮行,此过程称为左顶刮。直至转至图d位置左顶刮结束,齿轮2的齿1左齿面退出啮合,此时,齿轮2的左齿面进入单齿啮合区,右齿面为双面啮合区,左啮合线上有一个接触点,右啮合线上有两个接触点,齿轮对共有三个接触点。齿轮转至图e位置时,齿轮2的齿3的左齿面进入啮合,齿轮 2的左齿面进入双齿啮合区,左右啮合线上各有两个接触点,齿轮对共有四个接触点。转至图f位置时,齿轮2的齿1的右齿面将要退出啮合,右齿面进入单齿啮合区,左齿面仍为双齿啮合区,右啮合线上有一个接触点,左啮合线上有两个接触点,齿轮对共有三个接触点。齿轮转至图g位置时,齿轮2的齿3右齿面发生右顶刮,齿轮的双啮过程完成一个周期。
组合这些刚度分量以获得一对外齿轮的整体啮合刚度。对于一对接触比在1和2之间的正齿轮,交替地进行顶刃啮合、一对和两对齿啮合。外齿轮对的顶刃啮合、单齿对啮合和双齿对啮合持续时间表示为齿轮2的旋转角位移,顶刃啮合阶段齿轮2的转角为:
Figure GDA0003060890810000121
Figure GDA0003060890810000122
为齿轮2的转角,Δfpb为齿轮2右齿面的基节偏差,ra2为齿轮2齿顶圆半径,αa2为齿轮2的齿顶圆压力角。由图7可知,齿轮2 在双齿啮合区的转角为:
Figure GDA0003060890810000123
Figure GDA0003060890810000124
为齿轮2在双齿啮合区的转角,ε为齿轮对的重合度,z2为齿轮2齿数,pb2为齿轮2的基圆齿距,rb2为齿轮2的基圆半径。齿轮2在单齿啮合区的转角为:
Figure GDA0003060890810000125
Figure GDA0003060890810000126
为齿轮2在单齿啮合区的转角。取双齿啮合时的综合刚度为 kd,单齿啮合时的综合刚度为ks,顶刃啮合阶段的刚度为ka,则齿轮在右啮合线上的总刚度为:
Figure GDA0003060890810000127
刚度kR为齿轮转角的函数即
Figure GDA0003060890810000128
左啮合线上的刚度与右啮合线上的刚度相差一个相位角γar,左啮线上的齿轮刚度表示为
Figure GDA0003060890810000129
γar为右齿面的齿顶到左齿面齿根的间距角;左啮线上的齿轮刚度kL为:
Figure GDA0003060890810000131
kd为双齿啮合时的综合刚度,ksg为单齿啮合时的综合刚度,ka为顶刃啮合阶段的刚度。对于右齿面来说,单齿啮合区时的综合刚度 ksg即为该接触点的刚度。
Figure GDA0003060890810000132
其中下标1和2分别代表测量齿轮和产品齿轮。即kb1、kb2为分别表示测量齿轮和产品齿轮的弯曲刚度,ks1、ks2分别表示测量齿轮和产品齿轮的剪切刚度,kh为轮齿齿面接触刚度。
对于右齿面来说,顶刃啮合时的综合刚度为ka为顶刃啮合时右齿面两个接触点刚度的叠加,但是顶刃啮合阶段齿轮2的齿顶在齿轮 1的齿面上刮行,所以是齿轮2的齿顶处的刚度与齿轮1齿面不同点处的刚度叠加,其表达式如下
ka=ka1+ka2 (27)
其中ka1为第1个接触点处的刚度,ka2为第2个接触点处的刚度,轮齿在第1个接触点处为顶刃啮合,即产品齿轮的齿顶在测量齿轮的齿面上刮行,所以ka1为:
Figure GDA0003060890810000133
kb1,1、ks1,1为分别表示测量齿轮在第1个啮合点处的轮齿弯曲和刚度,,
Figure GDA0003060890810000134
分别表示产品齿轮在第1个啮合点处的轮齿剪切刚度,在这个阶段产品齿轮的齿顶在测量齿轮的齿面上刮行,即产品齿轮在该阶段始终是齿顶与测量齿轮接触,产品齿轮齿顶处的刚度即齿轮转角取
Figure GDA0003060890810000141
时的刚度,kh1为在第1个啮合点处的轮齿齿面接触刚度。
第二个接触点出的刚度ka2为:
Figure GDA0003060890810000142
kb1,2、kb2,2为分别表示测量齿轮和产品齿轮在第2个啮合点处的轮齿弯曲刚度,ks1,2、ks2,2分别表示测量齿轮和产品齿轮在第i个啮合点处的轮齿剪切刚度,kh,2为在第2个啮合点处的轮齿齿面接触刚度。
双齿啮合时右齿面的综合刚度为kd为双齿啮合时两个接触点刚度的叠加,其表达式如下:
Figure GDA0003060890810000143
kd1和kd2分别为第一个接触点和第二个接触点处的综合刚度,i 代表第i个啮合点,kb1,i、kb2,i为分别表示测量齿轮和产品齿轮在第i 个啮合点处的轮齿弯曲刚度,ks1,i、ks2,i分别表示测量齿轮和产品齿轮在第i个啮合点处的轮齿剪切刚度,kh为在第i个啮合点处的轮齿齿面接触刚度。
将kR代入式3便可求得双面啮合中右齿面沿啮合线方向的变形量δR,变形量是齿轮转角的函数
Figure GDA0003060890810000144
同样的可求得双面啮合时左齿面沿啮合线方向的变形量δL,左齿面变形量同样是齿轮转角的函数但与右齿面相差一个相位角γar,左齿面的变形量可表示为
Figure GDA0003060890810000145
将δL和δR带入式(1)便可求得齿轮双面啮合时在弹簧力作用下的径向变形量关于齿轮转角的函数。
通过Matlab对上述双面啮合测量中测量力影响下的轮齿变形算法进行仿真。以模数为0.5mm的钢制齿轮在10N径向力作用下的变形为例,所用齿轮参数如下:
Figure GDA0003060890810000151
单个轮齿在齿轮转动时力的作用点从齿根移动到齿顶的过程中变形量逐渐变大,刚度逐渐减小。双面啮合时单侧齿面的啮合刚度为顶刃啮合阶段、双齿啮合阶段和单齿啮合阶段的刚度的叠加,双面啮合的径向变形量是左右齿面变形量按相位的叠加。

Claims (1)

1.齿轮双面啮合测量中测量力引入轮齿变形计算方法,其特征在于,步骤如下:
1)轮齿变形量由三部分组成:弯曲变形、剪切变形和接触变形;弯曲和剪切模型由悬臂梁模型计算,接触变形由赫兹公式计算;轮齿的弯曲变形量δb和剪切变形量δs由下式表示
Figure FDA0002949468610000011
Figure FDA0002949468610000012
其中n为轮齿分成的微元数,δb为弯曲变形,δs为剪切变形,Fi为沿啮合线方向的给定作用力,σi为微元部分宽度,
Figure FDA0002949468610000013
为微元两个侧面截面模量的均值,G为轮齿的剪切模量,
Figure FDA0002949468610000014
为微元两个侧面截面积的均值,μ为泊松比,E为弹性模量,α'x为载荷方向与垂直方向的夹角,li表示微元i与载荷作用点之间沿X轴方向的距离;
轮齿的接触变形由赫兹公式计算;
Figure FDA0002949468610000015
δh为接触变形,μ1为测量齿轮的泊松比,μ2为产品齿轮的泊松比,E1为测量齿轮的弹性模量,E2为产品齿轮的弹性模量,b为齿宽;
2)轮齿的刚度由弯曲刚度、剪切刚度和接触刚度组成;轮齿的弯曲刚度kb、剪切刚度ks和接触刚度kh分别为:
Figure FDA0002949468610000016
Figure FDA0002949468610000021
Figure FDA0002949468610000022
δb、δs和δh分别为轮齿的弯曲变形、剪切变形和接触变形,Fi为沿啮合线方向的给定作用力;
3)轮齿的单对齿啮合刚度是由参与啮合的测量齿轮和产品齿轮的刚度叠加;齿轮的单对齿啮合刚度ki为:
Figure FDA0002949468610000023
kb1、kb2为分别表示测量齿轮和产品齿轮的轮齿弯曲刚度,ks1、ks2分别表示测量齿轮和产品齿轮的轮齿剪切刚度,kh为轮齿齿面的接触刚度;
4)齿轮双面啮合时在一条啮合线上存在单齿啮合区、双齿啮合区和顶刃啮合区,单齿啮合区只有一对齿参与啮合,齿轮的啮合刚度等于单对齿啮合刚度;双齿啮合区和顶刃啮合区同时有两对齿参与啮合,齿轮的啮合刚度为两个啮合啮合点的刚度的叠加;
单齿啮合区齿轮啮合刚度为:
Figure FDA0002949468610000024
kb1、kb2为分别表示测量齿轮和产品齿轮的弯曲刚度,ks1、ks2分别表示测量齿轮和产品齿轮的剪切刚度,kh为轮齿齿面接触刚度;
双齿啮合区齿轮啮合刚度为:
Figure FDA0002949468610000031
kd1和kd2分别为第一个啮合点和第二个啮合点处的单对齿刚度,i代表第i个啮合点,kb1,i、kb2,i为分别表示测量齿轮和产品齿轮在第i个啮合点处的轮齿弯曲刚度,ks1,i、ks2,i分别表示测量齿轮和产品齿轮在第i个啮合点处的轮齿剪切刚度,kh,i为在第i个啮合点处的轮齿齿面接触刚度;
在顶刃啮合区齿轮啮合刚度为:
ka=ka1+ka2
其中ka1为第1个啮合点处的刚度,即顶刮处的刚度,ka2为第2个啮合点处的刚度;
Figure FDA0002949468610000032
kb1,1、ks1,1为分别表示测量齿轮在顶刮处的轮齿弯曲和刚度,
Figure FDA0002949468610000033
Figure FDA0002949468610000034
分别表示产品齿轮在顶刮处的轮齿剪切刚度,kh1为在顶刮处的轮齿齿面接触刚度;
第二个啮合点出的刚度ka2为:
Figure FDA0002949468610000035
kb1,2、kb2,2为分别表示测量齿轮和产品齿轮在第2个啮合点处的轮齿弯曲刚度,ks1,2、ks2,2分别表示测量齿轮和产品齿轮在第2个啮合点处的轮齿剪切刚度,kh,2为在第2个啮合点处的轮齿齿面接触刚度;
因此齿轮在右啮合线上的刚度为:
Figure FDA0002949468610000041
刚度kR为齿轮转角的函数即
Figure FDA0002949468610000042
左啮合线上的刚度与右啮合线上的刚度相差一个相位角γar,左啮线上的齿轮刚度表示为
Figure FDA0002949468610000043
γar为右齿面的齿顶到左齿面齿根的间距角;左啮线上的齿轮刚度kL为:
Figure FDA0002949468610000044
5)由一条啮合线上的齿轮刚度求得齿轮沿该啮合线方向的变形量,沿左右啮合线方向的变形量分别为δL和δR
Figure FDA0002949468610000045
Figure FDA0002949468610000046
FL、FR为径向力沿左右啮合线方向的分量,kL、kR为沿啮合线方向的综合刚度;
将左右啮合线上的轮齿变形量按啮合相位叠加,从而得到测量力引起的双啮中心距变化量δ:
Figure FDA0002949468610000047
其中α’为啮合角,δL、δR为沿左右两条啮合线方向上的变形量。
CN201910965151.3A 2019-10-11 2019-10-11 齿轮双面啮合测量中测量力引入轮齿变形计算方法 Active CN110657986B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910965151.3A CN110657986B (zh) 2019-10-11 2019-10-11 齿轮双面啮合测量中测量力引入轮齿变形计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910965151.3A CN110657986B (zh) 2019-10-11 2019-10-11 齿轮双面啮合测量中测量力引入轮齿变形计算方法

Publications (2)

Publication Number Publication Date
CN110657986A CN110657986A (zh) 2020-01-07
CN110657986B true CN110657986B (zh) 2021-07-02

Family

ID=69040660

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910965151.3A Active CN110657986B (zh) 2019-10-11 2019-10-11 齿轮双面啮合测量中测量力引入轮齿变形计算方法

Country Status (1)

Country Link
CN (1) CN110657986B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111442742B (zh) * 2020-04-16 2021-11-02 创新奇智(重庆)科技有限公司 齿轮检测设备及方法
PL435092A1 (pl) * 2020-08-25 2022-02-28 Jmj Synergia Spółka Z Ograniczoną Odpowiedzialnością Układ mechaniczny z łańcuchem obustronnym

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1633303A1 (ru) * 1989-04-19 1991-03-07 Институт проблем надежности и долговечности машин АН БССР Способ определени жесткости зубьев зубчатых колес
CN101770538B (zh) * 2010-01-15 2011-11-09 北京工业大学 含损伤性单齿故障圆柱直齿轮啮合刚度仿真分析方法
KR101289434B1 (ko) * 2011-07-07 2013-07-24 한국기계연구원 기어의 다물체 동역학적 해석을 위한 기어치 모델링 방법 및 이를 이용한 다물체 동역학 해석 방법
CN103577687B (zh) * 2013-09-23 2017-05-03 北京工业大学 一种局部缺陷齿轮啮合刚度的时变特性定量计算方法
CN103629315A (zh) * 2013-12-10 2014-03-12 吉林大学 一种双相多轴传动系统
CN106641105B (zh) * 2017-01-22 2020-01-10 北京工业大学 一种齿轮逆序啮合模型的建立方法
CN106989157B (zh) * 2017-03-24 2019-06-14 北京工业大学 基于运动学原理的渐开线圆柱齿轮双面啮合测量仿真方法
CN107729626B (zh) * 2017-09-26 2020-03-17 西安交通大学 一种基于概率分布的齿轮点蚀模型建模方法
CN109813508A (zh) * 2017-11-21 2019-05-28 宁波午虎机电技术有限公司 机电伺服机构刚度测试装置
CN109101737B (zh) * 2018-08-23 2020-08-18 西南交通大学 一种考虑温度影响的直齿圆柱齿轮时变啮合刚度计算方法
CN109783840B (zh) * 2018-11-30 2022-10-25 华南理工大学 一种直齿圆柱齿轮内啮合齿轮副时变啮合刚度的计算方法

Also Published As

Publication number Publication date
CN110657986A (zh) 2020-01-07

Similar Documents

Publication Publication Date Title
CN110657986B (zh) 齿轮双面啮合测量中测量力引入轮齿变形计算方法
Li Effects of machining errors, assembly errors and tooth modifications on loading capacity, load-sharing ratio and transmission error of a pair of spur gears
Chaari et al. Analytical modelling of spur gear tooth crack and influence on gearmesh stiffness
Palmer et al. Evaluation of methods for calculating effects of tip relief on transmission error, noise and stress in loaded spur gears
CN103577687A (zh) 一种局部缺陷齿轮啮合刚度的时变特性定量计算方法
CN112257192A (zh) 一种基于实测齿廓偏差的一齿切向综合偏差计算方法
CN107247856A (zh) 一种单滚柱包络环面蜗杆副时变啮合刚度解析方法
Iglesias et al. Planet position errors in planetary transmission: Effect on load sharing and transmission error
CN202251810U (zh) 分段曲线齿轮
Huang et al. Influence of system-inherent phase in spur gear multi-state mesh based on nonlinear dynamics
CN106989157B (zh) 基于运动学原理的渐开线圆柱齿轮双面啮合测量仿真方法
Liu et al. A novel method to predict static transmission error for spur gear pair based on accuracy grade
Stein et al. 3D involute gear evaluation-part II: deviations-basic algorithms for modern software validation
Bodzás Analysis of the effect of the addendum modification coefficient for contact surfaces of spur gear
JP5330837B2 (ja) 歯車対の設計装置、及び歯車対
Wu et al. Geometrical design of skew conical involute gear drives in approximate line contact
CN109376455A (zh) 一种弧齿锥齿轮载荷弹性接触变形数值计算方法
CN1194214C (zh) 一种齿轮误差测量方法及其装置
JP2000171238A (ja) 遊星歯車装置の性能予測方法
CN113536483B (zh) 一种渐开线齿轮参数的测算方法的应用
CN113010978A (zh) 一种基于动态仿真的航空直齿轮修形方法
Shen et al. The influence of external spur gear surface wear on the mesh stiffness
Karma et al. Analysis of effects of pressure angle and pitch error variation on transmission error in spur gear using interference volume method
CN113704899B (zh) 基于齿轮点云数据的齿轮副传动误差计算方法
CN113591238B (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