CN112036049A - 一种实际工况下斜齿轮副时变啮合刚度的快速计算方法 - Google Patents
一种实际工况下斜齿轮副时变啮合刚度的快速计算方法 Download PDFInfo
- Publication number
- CN112036049A CN112036049A CN202010966081.6A CN202010966081A CN112036049A CN 112036049 A CN112036049 A CN 112036049A CN 202010966081 A CN202010966081 A CN 202010966081A CN 112036049 A CN112036049 A CN 112036049A
- Authority
- CN
- China
- Prior art keywords
- meshing
- gear
- tooth
- pair
- calculating
- 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
- 238000004364 calculation method Methods 0.000 title claims description 16
- 238000000034 method Methods 0.000 claims abstract description 23
- 150000001875 compounds Chemical class 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 2
- 238000006073 displacement reaction Methods 0.000 abstract description 2
- 230000005540 biological transmission Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000001133 acceleration Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 238000005312 nonlinear dynamic Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007334 memory performance Effects 0.000 description 1
Images
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
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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)
- Gears, Cams (AREA)
- Gear Transmission (AREA)
Abstract
本发明的一种实际工况下斜齿轮副时变啮合刚度的快速计算方法,首先根据具体工况下斜齿轮重合度及某时刻该斜齿轮副的啮合位置,计算该时刻斜齿轮中参与啮合的轮齿上的接触线长度;从而将此时刻斜齿轮副中参与啮合的各个轮齿沿齿宽方向分解为厚度相同的无数个独立直齿切片;并计算该时刻斜齿轮副端面直齿切片的啮合半径;从而考虑齿轮副变位系数,重合度等,计算此刻斜齿轮副中,同时,对参与啮合的各个轮齿的所有直齿切片啮合副啮合刚度进行积分求和,计算此刻斜齿轮副各个啮合轮齿的啮合刚度及斜齿轮副的总啮合刚度;最后依此类推,计算实际工况下,各个离散时刻斜齿轮副的啮合刚度,从而快速计算出实际工况下该斜齿轮副的时变啮合刚度。
Description
技术领域
本发明属于齿轮传动系统动力学分析领域,特别涉及实际工况下斜齿轮副时变啮合刚度的快速计算方法。
技术背景
斜齿轮传动平稳、承载能力强且传动效率高,因此在汽车、航空等各个工业领域被大量应用,目前针对齿轮的研究也很多,但在工程应用中,还存在很多问题,比如齿轮传动的振动噪声问题。时变啮合刚度激励作为齿轮传动振动噪声的主要内部激励源,大量学者对其进行了研究,但仍然难以满足实际工程应用,比如实际工况下,齿轮系统的输入转速和扭矩均为时变,每个时刻齿轮副输入转速和扭矩均不相同,参与啮合齿轮的齿数时变,每个啮合轮齿上的啮合位置时变,啮合时刻的啮合力也时刻发生变化,而当前以有限元法及数值分析方法为代表的斜齿轮时变啮合刚度无法满足实际工况下时变啮合刚度的快速计算。其中有限元法效率低,难以针对实际时变转速和扭矩工况下斜齿轮副的时变啮合刚度进行快速求解;解析法也主要针对恒定输入扭矩和转速下直齿轮或斜齿轮时变啮合刚度的快速计算,忽略了实际工况下输入扭矩和转速的时变性。
发明内容
为解决现有技术存在的上述问题,本发明提供一种实际工况下斜齿轮副时变啮合刚度的快速计算方法,用以解决现有斜齿轮副时变啮合刚度存在的不足,能够实现快速计算实际工况下斜齿轮副的时变啮合刚度,为分析实际工况下齿轮系统的非线性动力学特性提供了基础。
一种实际工况下斜齿轮副时变啮合刚度的快速计算方法,其特征在于,包括以下步骤:
步骤S1:获取实际工况下斜齿轮副初始啮合位置、输入转速、输入扭矩及斜齿轮副的基本参数;
步骤S2:根据S1中斜齿轮副的基本参数,计算斜齿轮副的重合度以及已旋转的角度,并对比判断斜齿轮轴向重合度和端面重合度的大小;
步骤S3:根据S1、S2中斜齿轮副参数、重合度及旋转角度,判断此刻参与啮合的齿轮对数,计算该时刻每个参与啮合的轮齿对上的接触线长度;并将斜齿轮离散为无数个直齿切片,判断参与啮合的直齿切片的啮合半径范围,基于啮合半径,计算直齿切片的刚度;
步骤S4:根据S3中直齿切片的刚度计算斜齿轮副中参与啮合的各个轮齿对的啮合刚度;
步骤S5:根据S4中此时刻斜齿轮参与啮合的各个轮齿对的啮合刚度,计算斜齿轮副的啮合刚度。
步骤S6:根据S5中该时刻斜齿轮副的啮合刚度,依此类推,计算实际工况下,各个离散时刻斜齿轮副的啮合刚度,即可快速计算出实际工况下该斜齿轮副的时变啮合刚度。
进一步的,步骤S3包括如下步骤:
步骤S31:判断此刻参与啮合的轮齿对数;
步骤S32:计算此时刻各啮合轮齿对上啮合线的长度;
步骤S33:将斜齿轮沿齿宽方向离散为无数个直齿切片;
步骤S34:计算此时刻参与啮合的直齿切片的啮合半径范围;
步骤S35:计算此时刻参与啮合的轮齿的综合刚度;
进一步的,步骤S31中参与啮合的轮齿对数判断公式为:
式中,M为斜齿轮副中参与啮合的轮齿对数;εβ为斜齿轮端面重合度和轴向重合度中数值较小的一个;εα为两个重合度中的另外一个;εγ为斜齿轮的整体重合度;zp为主动齿轮的齿数;Θ为主动齿轮旋转的角度,计算公式如下:Θ=ωptn
式中,ωp为主动齿轮旋转速度;tn为工况开始至此刻的时间;X为所选基准轮齿在所有啮合轮齿中的序号;zM为所有其他参与啮合的轮齿的序号;Pbt为基圆端面齿距,计算公式为:
Pbt=πmn/cosβ
式中,εmin=min(εε,εβ)和εmax=max(εε,εβ)分别为端面重合度和轴向重合度中较小和较大值,Lmax该斜齿轮副上的最长接触线长度,计算公式为:
进一步的,步骤S34中此时刻参与啮合的主动齿轮和被动齿轮的直齿切片的啮合半径范围为:
式中主动齿轮和被动齿轮上直齿切片啮合半径的为:
(1)当εα<εβ时,主动齿轮:
从动齿轮:
(2)当εα>εβ时,主动齿轮:
从动齿轮:
上述式中:
αt=arctan(tanα/cosβ)
式中,N1N2为轮齿副的作用线;B1B2为轮齿副实际啮合线;N1B1为作用线起点到实际啮合线起点的距离;N1B2为作用线起点到实际啮合线终点的距离;Rap、Rag分别为主动和被动齿轮的齿顶圆半径;OpB1、OgB1分别为主动齿轮Op和从动齿轮中心Og到B1的距离;Rp、Rg分别为主动和被动齿轮的分度圆半径;αt为齿轮端面压力角;α为分度圆压力角;β为齿轮螺旋角。
进一步的,步骤S35中此时刻参与啮合的主动齿轮和被动齿轮的直齿切片的啮合刚度计算公式为:
Kji(rji)=(A0+A1Xj)+(A2+A3Xj)(Rji-rj)m/(1+Xj)
其中:
式中,Kji(rji)为第i个(i=1,2...,N)主/被动轮齿(j=p,g)的端面直齿切片刚度;Xj为主被动轮齿的齿顶修正系数,m为法向模数;Rji为此刻第i个(i=1,2...,N)主被动轮齿(j=p,g)端面啮合半径;rj、zj根据j=p,j=g分别表示为主动齿轮p、被动齿轮g的节圆半径和齿数。
结合S34中此时刻的主被动齿轮上参与啮合的轮齿切片的啮合半径范围,计算参与啮合的各个轮齿的刚度为:
进一步的,步骤S4包括以下两个步骤:
步骤S41:计算此刻斜齿轮上各个啮合轮齿对的啮合刚度;
步骤S42:计算此刻斜齿轮副的啮合刚度。
进一步的,在步骤S41中:
根据主动轮齿和被动轮齿啮合的串联关系,计算斜齿轮上参与啮合的各个轮齿对的啮合刚度:
其中,Kpi和Kgi分别为第i个(i=1,2...,N)轮齿中主动齿轮的刚度,C(Fi)为载荷分配系数。
进一步的,在步骤S42中:
根据斜齿轮上参与啮合的各个轮齿对的啮合刚度,依据各对啮合轮齿对的并联关系,则斜齿轮副的啮合刚度为:
式中,λx(x=1,2,…,N)表示各齿轮副的载荷分配比,C(λxF)为对应的载荷系数。
本发明的有益效果是:
本发明的一种实际工况下斜齿轮副时变啮合刚度的快速计算方法,首先根据具体工况下斜齿轮重合度及某时刻该斜齿轮副的啮合位置,计算该时刻斜齿轮中参与啮合的轮齿上的接触线长度;从而将此时刻斜齿轮副中参与啮合的各个轮齿沿齿宽方向分解为厚度相同的无数个独立直齿切片;并计算该时刻斜齿轮副端面直齿切片的啮合半径;从而考虑齿轮副变位系数,重合度等,计算此刻斜齿轮副中,同时,对参与啮合的各个轮齿的所有直齿切片啮合副啮合刚度进行积分求和,计算此刻斜齿轮副各个啮合轮齿的啮合刚度及斜齿轮副的总啮合刚度;最后依此类推,计算实际工况下,各个离散时刻斜齿轮副的啮合刚度,从而快速计算出实际工况下该斜齿轮副的时变啮合刚度。
本发明根据啮合点半径计算该时刻直齿切片刚度,并判断直齿切片是否参与啮合以及参与啮合的直齿切片的啮合半径范围,通过积分求和快速计算某时刻斜齿轮副的综合啮合刚度,使得离散时刻下快速计算各瞬时齿轮副的啮合刚度成为现实,从而为快速计算实际时变工况下斜齿轮副的时变啮合刚度成为现实。
综上,本发明的一种实际工况下斜齿轮副时变啮合刚度的快速计算方法,能够实现快速计算实际工况下斜齿轮副的时变啮合刚度,为分析实际工况下齿轮系统的非线性动力学特性提供了基础。
附图说明
图1为端面直齿切片啮合简图。
图2为端面直齿切片啮合过程简图。
图3为加速工况下的输入转速和扭矩情况。
图4为斜齿轮作用面简图。
图5为斜齿轮啮合过程原理图。
图6为斜齿轮啮合过程接触线变化图。
图7为实际加速工况下斜齿轮副四个时刻点刚度图。
具体实施方案
为更好的解释本发明内容,下面结合附图,通过具体实施计算案例,对本发明进行详细说明。
本发明的一种实际工况下斜齿轮副时变啮合刚度的快速计算方法,将斜齿轮离散为无数个直齿切片,根据实际工况下各个时刻斜齿轮的啮合情况及接触线长度,判断各轮齿上所有直齿切片的啮合半径范围,通过积分求和以及主被动轮齿的并联关系,计算此刻各轮齿的啮合刚度,然后串联求和各轮齿刚度即得到该时刻此斜齿轮副的时变啮合刚度,以此类推,得到实际工况下,斜齿轮副的时变啮合刚度。
步骤S1:获取实际工况下斜齿轮副初始啮合位置、输入转速、输入扭矩及斜齿轮副齿数、模数、压力角、齿宽、螺旋角、修正系数等齿轮基本参数,如表1所示。假设,工况开始时,某一齿轮齿顶刚好啮入,则假设此轮齿为此时刻参与啮合的第X个轮齿。
表1斜齿轮副基本参数
步骤S2:根据S1中斜齿轮副的基本参数,计算斜齿轮副的重合度以及已旋转的角度,并对比判断斜齿轮轴向重合度和端面重合度的大小;
根据步骤S1中齿轮的基本参数,计算齿轮的重合度为:
εα=[zp(tanαap-tanαt')+zg(tanαag-tanαt')]/2π
εβ=Bsinβ/(πmn)
εγ=εα+εβ
αaj=arccos(Rbj/Raj)
αt′=arctan(tanα/cosβ)
上述式中,εα、εβ分别为端面重合度和轴向重合度;B为齿宽;β为螺旋角;mn为法向模数;αaj为齿轮的齿顶圆压力角(当j=p时代表主动齿轮,当j=g时代表从动齿轮);Raj和Rbj分别为齿顶圆和基圆半径(当j=p时代表主动齿轮,当j=g时代表从动齿轮);α为法向压力角。根据表1中的参数,计算齿轮的重合度为:εα=1.7572、εβ=1.6863和εγ=3.4435。此时的轴向重合度εα大于端面重合度εβ。
步骤S31:判断此刻参与啮合的轮齿对数;
式中,M为斜齿轮副中参与啮合的轮齿对数;Θ为主动齿轮旋转的角度,由S1及工况图(如图3所示)可知,随意取一个点,假设tn=0.38,则:Θ=ωptn=0.002
M=floor(εγ)+1=4
即此时有四个轮齿参与啮合,所选基准轮齿刚啮入,则此基准轮齿前还有三个齿轮参与啮合;啮合情况如图1、4和图5所示,共有7个进程(ab、bc、cd、de、ef、fg、gh),第一个轮齿即将啮出时第四个轮齿才开始啮入。
步骤S32:计算此时刻各啮合轮齿对上啮合线的长度;
根据S1中的重合度和S31中参与啮合的轮齿个数以及所选基准轮齿在此刻参与啮合的轮齿中的序号(为第4个参与啮合的轮齿),则此情况下的最大接触线为:
则计算各个轮齿的啮合线的长度,先计算各轮齿的旋转角度为:
根据7个进程,各轮齿的旋转角度,结合接触线长度的计算公式,则各轮齿的接触线长度根据所在的进程范围,计算公式为:
L(Θ1)=2.1421×10-6
L(Θ2)=0.0123
L(Θ3)=0.0177
L(Θ4)=0.0054
L(Θ5)=0
L(Θ6)=0
L(Θ7)=0
上述七个轮齿,由于判断5、6、7齿不在啮合范围内,则此时刻这三个轮齿未参与啮合,所以接触线长度为0.
步骤S33:将斜齿轮沿齿宽方向离散为无数个直齿切片;
步骤S34:计算此时刻参与啮合的直齿切片的啮合半径范围;
步骤S34中此时刻参与啮合的主动齿轮和被动齿轮的直齿切片的啮合半径范围,计算公式如下:
上式中,RpMi中i=1,2,…,7代表参与啮合的七个啮合进程(即一个轮齿从啮入到啮出共经历了三个三齿啮合区(bc、ed、fg),四个四齿啮合区(ab、cd、ef、gh))。同理,则被动齿轮为:
结合此时εα=1.7572、εβ=1.6863,则:
主动齿轮上参与啮合的前四个轮齿的啮合半径范围为:
Rp1∈(RpM-min,RpM-max)=(0.0302,0.0376)
Rp2∈(RpM-min,RpM-max)=(0.0302,0.0416)
Rp3∈(RpM-min,RpM-max)=(0.0338,0.0388)
Rp4∈(RpM-min,RpM-max)=(0.0338,0.0429)
从动齿轮上参与啮合的四个轮齿的啮合半径范围为:
Rg1∈(RgM-min,RgM-max)=(0.0383,0.043)
Rg2∈(RgM-min,RgM-max)=(0.0391,0.043)
Rg3∈(RgM-min,RgM-max)=(0.0383,0.0393)
Rg4∈(RpM-min,RpM-max)=(0.0387,0.0393)
步骤S35:计算此时刻参与啮合的轮齿的综合刚度;
基于参与啮合的主动齿轮和被动齿轮的直齿切片的啮合刚度计算公式为:Kji(rji)=(A0+A1Xj)+(A2+A3Xj)(Rji-rj)m/(1+Xj)
其中:
则主动齿轮和被动齿轮参数为:A0p=30.5004,A0g=31.7262,A1p=26.0268,A1g=24.9635,A2p=-17.4520,A2g=-18.3166,A3p=-21.8459,A3g=-21.7940。
结合S34中此时刻的主被动齿轮上参与啮合的轮齿切片的啮合半径范围,计算参与啮合的各个轮齿的刚度,其中各主动齿轮轮齿的刚度为:
各从动齿轮轮齿的刚度为:
步骤S41:计算此刻斜齿轮上各个啮合轮齿对的啮合刚度;
根据主动轮齿和被动轮齿啮合的串联关系,计算斜齿轮上参与啮合的各个轮齿对的啮合刚度:
步骤S42:计算此刻斜齿轮副的啮合刚度:
根据斜齿轮上参与啮合的各个轮齿对的啮合刚度,依据各对啮合轮齿对的并联关系,则斜齿轮副的啮合刚度为:
式中,λx(x=1,2,…,N)表示各齿轮副的载荷分配比,C(λxF)为对应的载荷系数。
进一步的,根据以上斜齿轮啮合刚度的计算方法,即可得到各离散时刻斜齿轮的啮合刚度,从而得到此工况下斜齿轮副的时变啮合刚度,本实施例共选择5个时间点,除了上述第一个点的刚度外,其余各四个时间点分别为t1=10s、t1=11s、t1=198s、t1=200s,对应的主被动轮齿刚度分别为:
齿轮副综合啮合刚度分别为:
Kpair(t=10)=1.728046×108、Kpair(t=11)=1.785617×108
Kpair(t=198)=1.671155×108、Kpair(t=200)=1.73723×108
同理,可计算其他时刻下齿轮啮合刚度。
本发明实施例得到的全油门工况下斜齿轮副主动轮啮合刚度如图7。实际加速工况下的输入转速和扭矩情况如图3,该工况下共分为20000个离散时刻点,此斜齿轮副的时变啮合刚度如图7所示,用32G内存性能的电脑,3300秒即计算出该工况下的时变啮合刚度,证明本发明能快速分析计算出实际工况下斜齿轮副的时变啮合刚度。
以上结合附图对本发明的实施例的技术方案进行完整描述,需要说明的是所描述的实施例仅仅是本发明的一部分实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
Claims (9)
1.一种实际工况下斜齿轮副时变啮合刚度的快速计算方法,其特征在于,包括以下步骤:
步骤S1:获取实际工况下斜齿轮副初始啮合位置、输入转速、输入扭矩及斜齿轮副的基本参数;
步骤S2:根据S1中斜齿轮副的基本参数,计算斜齿轮副的重合度以及已旋转的角度,并对比判断斜齿轮轴向重合度和端面重合度的大小;
步骤S3:根据S1、S2中斜齿轮副参数、重合度及旋转角度,判断此刻参与啮合的齿轮对数,计算该时刻每个参与啮合的轮齿对上的接触线长度;并将斜齿轮离散为无数个直齿切片,判断参与啮合的直齿切片的啮合半径范围,基于啮合半径,计算直齿切片的刚度;
步骤S4:根据S3中直齿切片的刚度计算斜齿轮副中参与啮合的各个轮齿对的啮合刚度;
步骤S5:根据S4中此时刻斜齿轮参与啮合的各个轮齿对的啮合刚度,计算斜齿轮副的啮合刚度。
步骤S6:根据S5中该时刻斜齿轮副的啮合刚度,依此类推,计算实际工况下,各个离散时刻斜齿轮副的啮合刚度,即可快速计算出实际工况下该斜齿轮副的时变啮合刚度。
2.根据权利要求1所述的实际工况下斜齿轮副时变啮合刚度的快速计算方法,其特征在于,步骤S3包括如下步骤:
步骤S31:判断此刻参与啮合的轮齿对数;
步骤S32:计算此时刻各啮合轮齿对上啮合线的长度;
步骤S33:将斜齿轮沿齿宽方向离散为无数个直齿切片;
步骤S34:计算此时刻参与啮合的直齿切片的啮合半径范围;
步骤S35:计算此时刻参与啮合的轮齿的综合刚度。
5.根据权利要求4所述的实际工况下斜齿轮副时变啮合刚度的快速计算方法,其特征在于,步骤S34中此时刻参与啮合的主动齿轮和被动齿轮的直齿切片的啮合半径范围为:
式中主动齿轮和被动齿轮上直齿切片啮合半径的为:
(1)当εα<εβ时,主动齿轮:
从动齿轮:
(2)当εα>εβ时,主动齿轮:
从动齿轮:
上述式中:
αt=arctan(tanα/cosβ)
式中,N1N2为轮齿副的作用线;B1B2为轮齿副实际啮合线;N1B1为作用线起点到实际啮合线起点的距离;N1B2为作用线起点到实际啮合线终点的距离;Rap、Rag分别为主动和被动齿轮的齿顶圆半径;OpB1、OgB1分别为主动齿轮Op和从动齿轮中心Og到B1的距离;Rp、Rg分别为主动和被动齿轮的分度圆半径;αt为齿轮端面压力角;α为分度圆压力角;β为齿轮螺旋角。
6.根据权利要求5所述的实际工况下斜齿轮副时变啮合刚度的快速计算方法,其特征在于,,步骤S35中此时刻参与啮合的主动齿轮和被动齿轮的直齿切片的啮合刚度计算公式为:
Kji(rji)=(A0+A1Xj)+(A2+A3Xj)(Rji-rj)m/(1+Xj)
其中:
式中,Kji(rji)为第i个(i=1,2...,N)主/被动轮齿(j=p,g)的端面直齿切片刚度;Xj为主被动轮齿的齿顶修正系数,m为法向模数;Rji为此刻第i个(i=1,2...,N)主被动轮齿(j=p,g)端面啮合半径;rj、zj根据j=p,j=g分别表示为主动齿轮p、被动齿轮g的节圆半径和齿数。
结合S34中此时刻的主被动齿轮上参与啮合的轮齿切片的啮合半径范围,计算参与啮合的各个轮齿的刚度为:
7.根据权利要求2所述的实际工况下斜齿轮副时变啮合刚度的快速计算方法,其特征在于,步骤S4包括以下两个步骤:
步骤S41:计算此刻斜齿轮上各个啮合轮齿对的啮合刚度;
步骤S42:计算此刻斜齿轮副的啮合刚度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010966081.6A CN112036049B (zh) | 2020-09-15 | 2020-09-15 | 一种实际工况下斜齿轮副时变啮合刚度的快速计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010966081.6A CN112036049B (zh) | 2020-09-15 | 2020-09-15 | 一种实际工况下斜齿轮副时变啮合刚度的快速计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112036049A true CN112036049A (zh) | 2020-12-04 |
CN112036049B CN112036049B (zh) | 2024-04-23 |
Family
ID=73589352
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010966081.6A Active CN112036049B (zh) | 2020-09-15 | 2020-09-15 | 一种实际工况下斜齿轮副时变啮合刚度的快速计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112036049B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113051681A (zh) * | 2021-03-25 | 2021-06-29 | 天津职业技术师范大学(中国职业培训指导教师进修中心) | 一种齿轮副啮合刚度计算方法及终端设备 |
CN113051682A (zh) * | 2021-03-25 | 2021-06-29 | 天津职业技术师范大学(中国职业培训指导教师进修中心) | 斜齿圆柱齿轮副热弹啮合刚度计算方法和装置 |
CN113378312A (zh) * | 2021-05-25 | 2021-09-10 | 武汉理工大学 | 一种考虑基圆与齿根圆相对位置的斜齿轮时变啮合刚度计算方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008175694A (ja) * | 2007-01-18 | 2008-07-31 | Fuji Heavy Ind Ltd | 歯車対の評価装置、評価プログラム、及びこれを用いた歯車対の評価方法 |
JP2010174983A (ja) * | 2009-01-29 | 2010-08-12 | Ogino Kogyo Kk | 揺動型歯車装置 |
JP2012043347A (ja) * | 2010-08-23 | 2012-03-01 | Ricoh Co Ltd | 歯車設計支援方法、歯車設計支援プログラムを記録した記録媒体及び歯車設計支援装置 |
CN105784360A (zh) * | 2016-05-12 | 2016-07-20 | 重庆长安汽车股份有限公司 | 一种基于啮合接触线长度变化确定齿轮啮合动刚度的方法 |
CN107391876A (zh) * | 2017-08-11 | 2017-11-24 | 东北大学 | 斜齿轮副时变啮合刚度计算方法 |
CN107798200A (zh) * | 2017-11-10 | 2018-03-13 | 西安电子科技大学 | 一种考虑轴向变形的斜齿圆柱齿轮时变啮合刚度计算方法 |
-
2020
- 2020-09-15 CN CN202010966081.6A patent/CN112036049B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008175694A (ja) * | 2007-01-18 | 2008-07-31 | Fuji Heavy Ind Ltd | 歯車対の評価装置、評価プログラム、及びこれを用いた歯車対の評価方法 |
JP2010174983A (ja) * | 2009-01-29 | 2010-08-12 | Ogino Kogyo Kk | 揺動型歯車装置 |
JP2012043347A (ja) * | 2010-08-23 | 2012-03-01 | Ricoh Co Ltd | 歯車設計支援方法、歯車設計支援プログラムを記録した記録媒体及び歯車設計支援装置 |
CN105784360A (zh) * | 2016-05-12 | 2016-07-20 | 重庆长安汽车股份有限公司 | 一种基于啮合接触线长度变化确定齿轮啮合动刚度的方法 |
CN107391876A (zh) * | 2017-08-11 | 2017-11-24 | 东北大学 | 斜齿轮副时变啮合刚度计算方法 |
CN107798200A (zh) * | 2017-11-10 | 2018-03-13 | 西安电子科技大学 | 一种考虑轴向变形的斜齿圆柱齿轮时变啮合刚度计算方法 |
Non-Patent Citations (2)
Title |
---|
FANG GUO ET AL: "The statistical analysis of the dynamic performance of a gear system considering random manufacturing errors under different levels of machining precision", JOURNAL OF MULTI-BODY DYNAMICS, vol. 234, no. 1, 19 July 2019 (2019-07-19), pages 3 - 18 * |
刘文;李锐;张晋红;林腾蛟;杨云;: "斜齿轮时变啮合刚度算法修正及影响因素研究", 湖南大学学报(自然科学版), no. 02, 25 February 2018 (2018-02-25), pages 1 - 10 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113051681A (zh) * | 2021-03-25 | 2021-06-29 | 天津职业技术师范大学(中国职业培训指导教师进修中心) | 一种齿轮副啮合刚度计算方法及终端设备 |
CN113051682A (zh) * | 2021-03-25 | 2021-06-29 | 天津职业技术师范大学(中国职业培训指导教师进修中心) | 斜齿圆柱齿轮副热弹啮合刚度计算方法和装置 |
CN113051682B (zh) * | 2021-03-25 | 2022-08-12 | 天津职业技术师范大学(中国职业培训指导教师进修中心) | 斜齿圆柱齿轮副热弹啮合刚度计算方法和装置 |
CN113378312A (zh) * | 2021-05-25 | 2021-09-10 | 武汉理工大学 | 一种考虑基圆与齿根圆相对位置的斜齿轮时变啮合刚度计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112036049B (zh) | 2024-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112036049A (zh) | 一种实际工况下斜齿轮副时变啮合刚度的快速计算方法 | |
JP4429390B2 (ja) | 歯車装置システム | |
CN1013217B (zh) | 环形齿轮/小齿轮设计 | |
US6709356B2 (en) | Vehicle transmission with progressive stepping | |
CN107391876A (zh) | 斜齿轮副时变啮合刚度计算方法 | |
EP1561972A3 (en) | Speed change device | |
CN111062103B (zh) | 一种优化参数的冶金起重机减速器齿轮寿命评估方法 | |
DE102010047143A1 (de) | Planetengetriebe, Planetendifferenzial und Getriebe mit dem Planetendifferenzial | |
CN113051681B (zh) | 一种齿轮副啮合刚度计算方法及终端设备 | |
CN109271702B (zh) | 工程车发动机与液力变矩器逆向匹配的优化方法 | |
CN202690900U (zh) | 一种新型齿形的齿轮偏心传动机构 | |
CN106996436B (zh) | 齿轮机构 | |
EP1445511A1 (en) | Limited slip differential of variable gear ratio type | |
CN115809527A (zh) | 基于切片法的斜齿轮时变啮合刚度及传递误差计算方法 | |
EP1371848A1 (en) | Gear pump with spline function generated gear profile | |
KR20090059107A (ko) | 다단 감속 기어 | |
CN113868755B (zh) | 基于基节误差的斜齿轮啮合刚度计算方法 | |
CN113378312B (zh) | 一种斜齿轮时变啮合刚度计算方法 | |
CN108662116A (zh) | 一种新型双圆弧谐波齿形 | |
US20050166697A1 (en) | Gearing element and motor vehicle gear unit | |
EP1760363A2 (de) | Antriebsbaugruppe in Form eines Winkelgetriebes für eine Arbeitsmaschine | |
KR20220078788A (ko) | 사이클로이드 기어의 설계 방법 | |
KR100965102B1 (ko) | 2개의 제어수단과 기어결합체를 이용한 변속장치 | |
CN112507485A (zh) | 基于切片耦合理论的斜齿轮时变啮合刚度分析方法 | |
CN112395713A (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 |