首页 / 专利库 / 齿轮 / 基圆 / 一种考虑齿向误差的直齿轮啮合刚度计算方法

一种考虑齿向误差的直齿轮啮合刚度计算方法

阅读:910发布:2021-12-05

专利汇可以提供一种考虑齿向误差的直齿轮啮合刚度计算方法专利检索,专利查询,专利分析的服务。并且本 发明 提出了一种考虑齿向误差的直 齿轮 啮合 刚度 计算方法,在保证计算效率的同时,提高 直齿轮 啮合刚度计算 精度 。实现步骤为:设置直齿轮副的参数;计算直齿轮副的齿基刚度;获取直齿轮副中每个直齿轮沿齿宽方向的多个切片齿;计算切片齿的总误差、每个切片齿的总刚度以及每个直齿轮中相邻切片齿之间的耦合刚度;构建齿轮副中相互啮合的一对齿的切片耦合模型;推导 串联 弹簧 组中 变形 量间的关系式;推导的 力 与变形的平衡公式;计算切片齿的变形量;计算齿轮副的啮合力;计算齿轮副的时变啮合刚度。本发明考虑了直齿轮存在齿向修形和不对中误差,以及切片齿之间的耦合作用,计算精度高,可用于直齿轮的动态性能分析及优化设计。,下面是一种考虑齿向误差的直齿轮啮合刚度计算方法专利的具体信息内容。

1.一种考虑齿向误差的直齿轮啮合刚度计算方法,其特征在于包括如下步骤:
(1)设置直齿轮副的参数:
设置直齿轮副中主动齿轮p和从动齿轮g的齿宽均为B,模数均为m,主动齿轮p和从动齿轮g的弹性模量分别为Ep和Eg、主动齿轮p和从动齿轮g的齿数分别为z1和z2;
(2)计算直齿轮副的齿基刚度ktf:
计算直齿轮副中主动齿轮p的单齿啮合齿基刚度kfp和从动齿轮g的单齿啮合齿基刚度kfg,并通过kfp和kfg计算直齿轮副的齿基刚度ktf;
(3)获取直齿轮副中每个直齿轮沿齿宽方向的多个切片齿:
将直齿轮副中的主动齿轮p等效为沿齿宽方向的N个宽度为Δl的切片齿
同时将从动齿轮g等效为沿齿宽方向的N个宽度为Δl的切片齿
其中, 表示第i个切片齿, 表示第i个切片齿, N≥2;
(4)计算切片齿 和 的总误差、每个切片齿的总刚度以及每个直齿轮中相邻切片齿之间的耦合刚度:
(4a)计算切片齿 和 的总误差Ei,包括切片齿 的齿向修形误差elp、 的齿向修形误差elg和 与 的不对中误差em,Ei=elp+elg+em;
(4b)通过切片齿 的弯曲刚度 剪切刚度 径向压缩刚度 和赫兹接触
度Δkh,计算 的总刚度 通过切片齿 的弯曲刚度 剪切刚度 径向压缩刚度
和赫兹接触刚度Δkh,计算 的总刚度
(4c)计算切片齿 及其相邻切片齿 之间的耦合刚度 计算切片齿 及其相邻
切片齿 之间的耦合刚度
其中Cc表示耦合系数;
(5)构建齿轮副中相互啮合的一对齿的切片耦合模型S:
(5a)将切片齿 模拟为与自身总刚度 相等的主动齿轮切片齿压缩弹簧 得到由N个并行排布的主动齿轮切片齿压缩弹簧组成的集合 将切片齿
及其相邻的切片齿 之间的耦合模拟为与耦合刚度 相等的主动齿轮耦合压缩弹簧 得到由N-1个主动齿轮耦合压缩弹簧组成的集合
并将 与 通过 连接,得到主动齿轮p的切片耦合模型;
(5b)将切片齿 模拟为与自身总刚度 相等的从动齿轮切片齿压缩弹簧 得到由N个并行排布的从动齿轮切片齿压缩弹簧组成的集合 将切片齿
及其相邻的切片齿 之间的耦合模拟为与耦合刚度 相等的从动齿轮耦合压缩弹簧 得到由N-1个从动齿轮耦合压缩弹簧组成的集合
并将 与 通过 连接,得到从动齿轮g的切片耦合模型;
(5c)将主动齿轮p与从动齿轮g对应切片齿 与 之间的接触模拟为与对应切片齿与 之间的接触刚度Δkh相等的接触压缩弹簧 得到由N个接触压缩弹簧组成的集合(5d)将主动齿轮p的切片耦合模型中的 与从动齿轮g的切片耦合模型中的 通过连接,组成串联弹簧组 并将N个串联弹簧组组和成串联弹簧组集合
(5e)将 以及 和
组合成齿轮副中相互啮合的一对齿的切片耦合模型;
(6)推导串联弹簧组 中变形量间的关系式:
(6a)假设切片耦合模型中串联弹簧组 在Fi的作用下产生变形量为δit=δpi+δgi+δhi,N个串联弹簧组变形量的和即为齿轮副的啮合变形量 其中δpi为 的变形量,δgi为 的变形量,δhi为 的变形量;
(6b)通过Fi、δpi、δgi、δhi和串联弹簧组 中各压缩弹簧的刚度所构建的平衡方程以及δit=δpi+δgi+δhi,推导串联弹簧组 中变形量间的关系式:
(7)推导 的力与变形的平衡公式:
(7a)通过Fi和δpi,以及主动齿轮p的切片耦合模型的主动齿轮切片齿压缩弹簧 和主动齿轮耦合压缩弹簧 的刚度 和 推导主动齿轮切片齿压缩弹簧 的变形量δpi和Fi的平衡方程,得到N个主动齿轮切片齿压缩弹簧的变形量和力的平衡方程:
(7b)将N个主动齿轮切片齿压缩弹簧的变形量和力的平衡方程转换为矩阵形式KcX=F,并将其作为 的力与变形的平衡公式,其中,Kc表示刚度矩阵,X和F分别表示变形矢量和施加的力矢量,
T
X=[δp1,δp2,…,δpi,…,δpN]
F=[F1,F2,…,Fi,…,FN]T
[·]T表示矩阵[·]的转置;
(8)计算齿轮副的时变啮合刚度:
(8a)设施加在齿轮副中相互啮合的一对齿的切片耦合模型S上的外力为Fe,力容忍误差为Fε,并令δ=0;
(8b)判断δ与切片齿 和 的总误差Ei的关系δ>Ei是否成立,若是,则切片齿 与切片齿 啮合 ,将主动齿轮p中 所有处于啮合状态的切片齿组合为集合
否则,则切片齿 与切片齿 不啮合,将主动齿轮p中所有
处于不啮合状态的切片齿 组合为集合 其中, 为主动齿
轮p中第w个处于啮合状态的切片齿,W为处于啮合状态的切片齿的总个数, 为主动齿轮p中第s个处于不啮合状态的切片齿,S为处于不啮合状态的切片齿的总个数,S+W=N;
(8c)当切片齿 与切片齿 啮合时,主动齿轮p中W个处于啮合状态的切片齿的变形量均为 根据距离集合Ap”区域最近的处于啮合状态的切片齿 的变形量计算处于不啮合状态的切片齿 的变形量δ″ps,并将计算得到的 的变形量δpi填补到X=[δp1,δp2,…,δpi,…,δpN]T的对应行,其中δ″ps的计算公式为:
其中,Γj(j+1)为第j个切片齿和第j+1个切片齿之间的变形传递系数;
(8d)计算齿轮副的啮合力Fm,并判断|Fe-Fm|≤Fε是否成立,若是,得到齿轮副的啮合变形量δ和齿对的载荷分布,并执行步骤(8e),否则,令 并执行步骤(8b);
(8e)根据齿轮副啮合变形量δ计算齿轮副的啮合刚度k,并计算一个啮合周期内各啮合位置的啮合刚度,即齿轮副的时变啮合刚度。
2.根据权利要求1所述的一种考虑齿向误差的直齿轮啮合刚度计算方法,其特征在于,步骤(2)中所述的主动直齿轮p的单齿啮合齿基刚度kfp、从动直齿轮g的单齿啮合齿基刚度kfg和直齿轮副的齿基刚度ktf,计算公式分别为:
其中, 为接触点处的压力,ufq为力到齿根的距离,Sfq为齿根宽度,L*,M*,P*,Q*为公式的系数,下标q取值为p和g,分别代表主动齿轮和从动齿轮,λp是主动齿轮p的齿基刚度修正系数,λg是从动齿轮g的齿基刚度修正系数。
3.根据权利要求1所述的一种考虑齿向误差的直齿轮啮合刚度计算方法,其特征在于,步骤(4a)中所述的切片齿 的齿向修形误差elp、 的齿向修形误差elg和 与 的不对中误差em,计算公式分别为:
其中Cβ表示齿向修形量;b0表示齿向修形曲线沿齿宽方向的顶点坐标,bi表示各切片沿齿宽方向的坐标;s表示决定齿向修形曲线形状的弯曲系数,βb表示齿轮基圆螺旋角,θx1、θx2、θy1和θy2分别为主动齿轮p和从动齿轮g的摆动自由度,ψ12为和主动齿轮转动方向有关的参数。
4.根据权利要求1所述的一种考虑齿向误差的直齿轮啮合刚度计算方法,其特征在于,步骤(4b)中所述的 的总刚度 和 的总刚度 计算公式分别为:
5.根据权利要求1所述的一种考虑齿向误差的直齿轮啮合刚度计算方法,其特征在于,步骤(8c)所述的第j个切片齿和第j+1个切片齿之间的变形传递系数,计算公式为:
其中,当下标大于N时,
6.根据权利要求1所述的一种考虑齿向误差的直齿轮啮合刚度计算方法,其特征在于,步骤(8d)中啮合力Fm,计算公式为:
F=KX
其中K为刚度矩阵,当单齿啮合时K=Kc,双齿啮合时
7.根据权利要求1所述的一种考虑齿向误差的直齿轮啮合刚度计算方法,其特征在于,步骤(8e)中齿对啮合刚度ktt和齿轮副啮合刚度k,计算公式分别为:

说明书全文

一种考虑齿向误差的直齿轮啮合刚度计算方法

技术领域

[0001] 本发明属于机械动学技术领域,涉及一种考虑齿向误差的直齿轮啮合刚度计算方法,可用于齿向误差较大的直齿轮的动态性能分析及优化设计。

背景技术

[0002] 齿轮是应用最广泛的旋转机械零件之一,齿轮系统的性能直接决定着旋转机械的性能。近几年来,高端数控装备正不断向着高速、高精度和长寿命的方向发展,并对直齿轮这一关键部件的动力学特性提出了更高的要求。而啮合刚度周期性的变化是传动系统的主要激励形式之一,这种啮合刚度周期性的变化称之为时变啮合刚度,其直接影响着传动过程中的动力学特性。当齿轮存在齿向误差时,会对其啮合刚度产生较大的影响,因此研究考虑齿向误差的直齿轮啮合刚度计算方法很有必要。
[0003] 目前,对于齿轮传动系统时变啮合刚度的计算方法,主要有解析法、有限元法和解析有限元法。解析法运用材料力学原理,具有计算效率高、使用方便等优点,广泛应用于齿轮副时变啮合刚度的计算,但是解析法通常忽略很多因素,如齿向修形和不对中情况,因此在计算精度上会受到影响。有限元法可以模拟实际齿廓,包括齿廓修形、加工误差、安装误差等,在齿轮副时变啮合刚度分析中具有较高的精度,然而,有限元方法需要对不同的齿轮进行重复建模且需要很多的计算资源。解析有限元法是解析法和有限元法相结合的一种计算方法,可快速计算各种齿形下的时变啮合刚度,如齿廓修形、薄腹板齿基、带孔齿基、裂纹等,但该方法需要专业人员进行建模并且对于不同的齿轮需要进行重复建模。
[0004] 从目前的公开的资料来看,对于存在齿向修形和不对中的斜齿轮和直齿轮,通常用采用解析法和切片理论计算齿轮副的啮合刚度。但是目前的方法通常认为切片齿之间是相互独立的,忽略了切片齿之间的耦合。例如申请公布号为CN 107153736 A,名称为“一种修正的考虑鼓向修形的齿轮副啮合特性分析方法”的中国专利申请,公开了一种修正的考虑鼓向修形的齿轮副啮合特性分析方法,该方法为获取齿轮副的基本参数及鼓向修形参数,将齿轮副的轮齿沿齿宽方向分解为N个独立且均匀的薄片齿轮,基于鼓向修形齿轮副的齿廓误差,采用考虑延长啮合影响的齿轮副啮合特性分析方法,计算得到每片薄片齿轮副的时变啮合刚度;通过三维绘图模拟软件建立含鼓向修形的啮合齿轮副三维模型,将三维模型导入到ANSYS软件建立三维有限元接触模型,求解齿轮整个啮合过程中的时变啮合刚度数据。该方法考虑到了非线性接触刚度、有限元修正基体刚度、延长啮合作用的影响,在一定程度上提高了齿轮啮合刚度的计算精度,但其没有考虑薄片齿轮之间的耦合作用对时变啮合刚度的影响,并且没有考虑轴线不对中的情况,因此对计算精度有一定的影响。

发明内容

[0005] 本发明的目的在于克服上述现有技术存在的缺陷,提出了一种考虑齿向误差的直齿轮啮合刚度计算方法,旨在保证计算效率的同时,提高存在齿向误差的直齿轮啮合刚度计算精度。
[0006] 为了实现上述目的,本发明采用的技术方案包括如下步骤:
[0007] (1)设置直齿轮副的参数:
[0008] 设置直齿轮副中主动齿轮p和从动齿轮g的齿宽均为B,模数均为m,主动齿轮p和从动齿轮g的弹性模量分别为Ep和Eg、主动齿轮p和从动齿轮g的齿数分别为z1和z2;
[0009] (2)计算直齿轮副的齿基刚度ktf:
[0010] 计算直齿轮副中主动齿轮p的单齿啮合齿基刚度kfp和从动齿轮g的单齿啮合齿基刚度kfg,并通过kfp和kfg计算直齿轮副的齿基刚度ktf;
[0011] (3)获取直齿轮副中每个直齿轮沿齿宽方向的多个切片齿:
[0012] 将直齿轮副中的主动齿轮p等效为沿齿宽方向的N个宽度为Δl的切片齿同时将从动齿轮g等效为沿齿宽方向的N个宽度为Δl的切片齿
其中, 表示第i个切片齿, 表示第i个切片齿, N≥2;
[0013] (4)计算切片齿 和 的总误差、每个切片齿的总刚度以及每个直齿轮中相邻切片齿之间的耦合刚度:
[0014] (4a)计算切片齿 和 的总误差Ei,包括切片齿 的齿向修形误差elp、 的齿向修形误差elg和 与 的不对中误差em,Ei=elp+elg+em;
[0015] (4b)通过切片齿 的弯曲刚度 剪切刚度 径向压缩刚度 和赫兹接触刚度Δkh,计算 的总刚度 通过切片齿 的弯曲刚度 剪切刚度 径向压缩
刚度 和赫兹接触刚度Δkh,计算 的总刚度
[0016] (4c)计算切片齿 及其相邻切片齿 之间的耦合刚度 计算切片齿 及其相邻切片齿 之间的耦合刚度
[0017]
[0018]
[0019] 其中Cc表示耦合系数;
[0020] (5)构建齿轮副中相互啮合的一对齿的切片耦合模型S:
[0021] (5a)将切片齿 模拟为与自身总刚度 相等的主动齿轮切片齿压缩弹簧 得到由N个并行排布的主动齿轮切片齿压缩弹簧组成的集合 将切片齿 及其相邻的切片齿 之间的耦合模拟为与耦合刚度 相等的主动齿轮耦合压缩弹簧 得到由N-1个主动齿轮耦合压缩弹簧组成的集合
并将 与 通过 连接,得到主动齿轮p的切片耦合模型;
[0022] (5b)将切片齿 模拟为与自身总刚度 相等的从动齿轮切片齿压缩弹簧 得到由N个并行排布的从动齿轮切片齿压缩弹簧组成的集合 将切片齿 及其相邻的切片齿 之间的耦合模拟为与耦合刚度 相等的从动齿轮耦合压缩弹簧 得到由N-1个从动齿轮耦合压缩弹簧组成的集合 并
将 与 通过 连接,得到从动齿轮g的切片耦合模型;
[0023] (5c)将主动齿轮p与从动齿轮g对应切片齿 与 之间的接触模拟为与对应切片齿 与 之间的接触刚度Δkh相等的接触压缩弹簧 得到由N个接触压缩弹簧组成的集合
[0024] (5d)将主动齿轮p的切片耦合模型中的 与从动齿轮g的切片耦合模型中的 通过 连接,组成串联弹簧组 并将N个串联弹簧组组和成串联弹簧组集合
[0025] ( 5 e ) 将 以 及 和组合成齿轮副中相互啮合的一对齿的切片耦合模型;
[0026] (6)推导串联弹簧组 中变形量间的关系式:
[0027] (6a)假设切片耦合模型中串联弹簧组 在力Fi的作用下产生变形量为N个串联弹簧组变形量的和即为齿轮副的啮合变形量 其中δpi
为 的变形量,δgi为 的变形量,δhi为 的变形量;
[0028] (6b)通过Fi、δpi、δgi、δhi和串联弹簧组 中各压缩弹簧的刚度所构建的平衡方程以及 推导串联弹簧组 中变形量间的关系式:
[0029]
[0030] (7)推导 的力与变形的平衡公式:
[0031] (7a)通过Fi和δpi,以及主动齿轮p的切片耦合模型的主动齿轮切片齿压缩弹簧和主动齿轮耦合压缩弹簧 的刚度 和 推导主动齿轮切片齿压缩弹簧 的变形量δpi和Fi的平衡方程,得到N个主动齿轮切片齿压缩弹簧的变形量和力的平衡方程:
[0032]
[0033] (7b)将N个主动齿轮切片齿压缩弹簧的变形量和力的平衡方程转换为矩阵形式KcX=F,并将其作为 的力与变形的平衡公式,其中,Kc表示刚度矩阵,X和F分别表示变形矢量和施加的力矢量,
[0034]
[0035] X=[δp1,δp2,…,δpi,…,δpN]T
[0036] F=[F1,F2,…,Fi,…,FN]T
[0037] [·]T表示矩阵[·]的转置;
[0038] (8)计算齿轮副的时变啮合刚度:
[0039] (8a)设施加在齿轮副中相互啮合的一对齿的切片耦合模型S上的外力为Fe,力容忍误差为Fε,并令δ=0;
[0040] (8b)判断δ与切片齿 和 的总误差Ei的关系δ>Ei是否成立,若是,则切片齿与切片齿 啮合,将主动齿轮p中所有处于啮合状态的切片齿组合为集合否则,则切片齿 与切片齿 不啮合,将主动齿轮p中所有处
于不啮合状态的切片齿 组合为集合 其中, 为主动齿轮
p中第w个处于啮合状态的切片齿,W为处于啮合状态的切片齿的总个数, 为主动齿轮p中第s个处于不啮合状态的切片齿,S为处于不啮合状态的切片齿的总个数,S+W=N;
[0041] (8c)当切片齿 与切片齿 啮合时,主动齿轮p中W个处于啮合状态的切片齿的变形量均为 根据距离集合Ap”区域最近的处于啮合状态的切片齿 的变形量计算处于不啮合状态的切片齿 的变形量δ”ps,并将计算得到的 的变形量δpi填补到X=[δp1,δp2,…,δpi,…,δpN]T的对应行,其中δ”ps的计算公式为:
[0042]
[0043] 其中,Γj(j+1)为第j个切片齿和第j+1个切片齿之间的变形传递系数;
[0044] (8d)计算齿轮副的啮合力Fm,并判断|Fe-Fm|≤Fε是否成立,若是,得到齿轮副的啮合变形量δ和齿对的载荷分布,并执行步骤(8e),否则,令 并执行步骤(8b);
[0045] (8e)根据齿轮副啮合变形量δ计算齿轮副的啮合刚度k,并计算一个啮合周期内各啮合位置的啮合刚度,即齿轮副的时变啮合刚度。
[0046] 本发明与现有技术相比,具有如下优点:
[0047] 本发明通过解析法和切片理论获取齿轮副的时变啮合刚度过程中,在计算切片齿和 的总误差时,考虑了对直齿轮啮合刚度有影响的齿向修形误差和不对中误差;在构建切片耦合模型时,将切片齿和相邻切片齿之间的耦合作用等效为弹簧,充分考虑了主动齿轮相邻切片齿之间的耦合作用和从动齿轮相邻切片齿之间的耦合作用,与现有技术相比,在保证计算效率的同时,有效提高了存在齿向误差时的直齿轮时变啮合刚度计算精度。附图说明
[0048] 图1为本发明的实现流程图
[0049] 图2为本发明实施例的齿基啮合刚度的参数定义图;
[0050] 图3为本发明实施例的切片齿示意图;
[0051] 图4为本发明实施例的不对中误差参数定义图;
[0052] 图5为本发明实施例的主动齿轮p的切片耦合模型示意图;
[0053] 图6为本发明实施例的从动齿轮g的切片耦合模型示意图;
[0054] 图7为本发明实施例的齿轮副中相互啮合的一对齿的切片耦合模型示意图;

具体实施方式

[0055] 下面结合附图和具体实施例,对本发明作进一步详细说明。
[0056] 参照图1.本发明包括如下步骤:
[0057] 步骤1)设置直齿轮副的参数:
[0058] 设置直齿轮副中主动齿轮p和从动齿轮g的齿宽均为B,模数均为m,主动齿轮p和从动齿轮g的弹性模量分别为Ep和Eg、主动齿轮p和从动齿轮g的齿数分别为z1和z2;
[0059] 本实施例的参数取值如表1所示:
[0060] 表1直齿轮副的参数
[0061]
[0062]
[0063] 步骤2)计算直齿轮副的齿基刚度ktf:
[0064] 计算直齿轮副中主动齿轮p的单齿啮合齿基刚度kfp和从动齿轮g的单齿啮合齿基刚度kfg,并通过kfp和kfg计算直齿轮副的齿基刚度ktf;
[0065] 其中主动直齿轮p的单齿啮合齿基刚度kfp、从动直齿轮g的单齿啮合齿基刚度kfg计算公式为:
[0066]
[0067] 其中,计算直齿轮副的齿基刚度所需的参数如图2所示,ufq为力到齿根的距离,Sfq为齿根宽度,, 为接触点处的压力,系数L*、M*、P*和Q*可按照多项式函数得到,X*代表系数L*、M*、P*和Q*,其中hfq=rfq/rintq;Ai、Bi、Ci、Di、Ei和Fi的值见表2,下标q取值为p和g,分别代表主动齿轮和从动齿轮;
[0068] 并通过kfp和kfg计算直齿轮副的齿基刚度ktf,计算公式为:
[0069]
[0070] λp是主动齿轮p的齿基刚度修正系数,λg是从动齿轮g的齿基刚度修正系数。
[0071] 表2 Ai,Bi,Ci,Di,Ei,Fi系数值
[0072]
[0073] 步骤3)获取直齿轮副中每个直齿轮沿齿宽方向的多个切片齿:
[0074] 将直齿轮副中的主动齿轮p等效为沿齿宽方向的N个宽度为Δl的切片齿同时将从动齿轮g等效为沿齿宽方向的N个宽度为Δl的切片齿
其中, 表示第i个切片齿, 表示第i个切片齿, N≥2,本实
施例中N=100;
[0075] 每个直齿轮等效为的切片齿如图3所示,图中Δl为每个切片齿的宽度,z轴方向为齿宽方向。
[0076] 步骤4)计算切片齿 和 的总误差、每个切片齿的总刚度以及每个直齿轮中相邻切片齿之间的耦合刚度:
[0077] 步骤4a)计算切片齿 和 的总误差Ei,包括切片齿 的齿向修形误差elp、 的齿向修形误差elg和 与 的不对中误差em,Ei=elp+elg+em;其中切片齿 的齿向修形误差elp、 的齿向修形误差elg和 与 的不对中误差em,计算公式分别为:
[0078]
[0079]
[0080]
[0081] 其中Cβ表示齿向修形量;b0表示齿向修形曲线沿齿宽方向的顶点坐标,bi表示各切片沿齿宽方向的坐标;s表示决定齿向修形曲线形状的弯曲系数,βb表示齿轮基圆螺旋角,θx1、θx2、θy1和θy2分别为如图4所示的主动齿轮p和从动齿轮的摆动自由度,ψ12可表示为:
[0082]
[0083] 其中α12(0≤α12≤2π)为齿轮副相对位置角,如图4所示,在本实施例中齿向修形量Cβ=20μm,不对中角度θy1=0.05°;
[0084] 步骤4b)通过切片齿 的弯曲刚度 剪切刚度 径向压缩刚度 和赫兹接触刚度Δkhi,计算 的总刚度 通过切片齿 的弯曲刚度 剪切刚度 径向
压缩刚度 和赫兹接触刚度Δkh,计算 的总刚度
[0085] 采用势能法,主动齿轮p和从动齿轮g中的每个切片齿的弯曲刚度 剪切刚度径向压缩刚度 和赫兹接触刚度Δkh具有相同的计算公式,计算公式分别为:
[0086]
[0087]
[0088]
[0089]
[0090]
[0091] 其中 为分度圆压力角, 为接触点处的压力角,rbq为基圆半径,θbq为基齿角的一半, 为齿顶高系数,γ和τ为积分范围,y1q为过渡曲线上任一点处的平坐标,y2q为渐开线上任一点处的水平坐标,Iy1q为过渡曲线上任意位置处的截面惯性矩,Iy2q为渐开线上任意位置处的截面惯性矩,Gq为剪切模量,Ay1q为过度曲线上任意位置处的截面面积,Ay2q为渐开线上任意位置处的截面面积,E为合成的弹性模量,下标q取值为p和g,分别代表主动齿轮和从动齿轮。
[0092] 步骤4c)计算切片齿 及其相邻切片齿 之间的耦合刚度 计算切片齿及其相邻切片齿 之间的耦合刚度
[0093]
[0094]
[0095] 其中Cc表示耦合系数;
[0096] 步骤5)构建齿轮副中相互啮合的一对齿的切片耦合模型S:
[0097] 步骤5a)将切片齿 模拟为与自身总刚度 相等的主动齿轮切片齿压缩弹簧得到由N个并行排布的主动齿轮切片齿压缩弹簧组成的集合 将切片齿 及其相邻的切片齿 之间的耦合模拟为与耦合刚度 相等的主动齿轮耦合压缩弹簧 得到由N-1个主动齿轮耦合压缩弹簧组成的集合
并将 与 通过 连接,得到如图5所示的主动
齿轮p的切片耦合模型;
[0098] 步骤5b)将切片齿 模拟为与自身总刚度 相等的从动齿轮切片齿压缩弹簧得到由N个并行排布的从动齿轮切片齿压缩弹簧组成的集合 将切片齿 及其相邻的切片齿 之间的耦合模拟为与耦合刚度 相等的从动齿轮耦合压缩弹簧 得到由N-1个从动齿轮耦合压缩弹簧组成的集合
并将 与 通过 连接,得到如图6所示的从动齿
轮g的切片耦合模型;
[0099] 步骤5c)将主动齿轮p与从动齿轮g对应切片齿 与 之间的接触模拟为与对应切片齿 与 之间的接触刚度Δkh相等的接触压缩弹簧 得到由N个接触压缩弹簧组成的集合
[0100] 步骤5d)将主动齿轮p的切片耦合模型中的 与从动齿轮g的切片耦合模型中的通过 连接,组成串联弹簧组 并将N个串联弹簧组组和成串联弹簧组集合
[0101] 步骤5e)将 以及 和组合成如图7所示的齿轮副中相互啮合的一对齿的切片耦
合模型;
[0102] 步骤6)推导串联弹簧组 中变形量间的关系式:
[0103] 步骤6a)假设切片耦合模型中串联弹簧组 在力Fi的作用下产生变形量为N个串联弹簧组变形量的和即为齿轮副的啮合变形的 其中δpi为 的变形量,δgi为 的变形量,δhi为 的变形量;
[0104] 步骤6b)通过Fi、δpi、δgi、δhi和串联弹簧组 中各压缩弹簧的刚度所构建的平衡方程 以及 推导串联弹簧组 中变形量间的关系式:
[0105]
[0106] 步骤7)推导 的力与变形的平衡公式;
[0107] 步骤7a)通过Fi和δpi,以及主动齿轮p的切片耦合模型的主动齿轮切片齿压缩弹簧和主动齿轮耦合压缩弹簧 的刚度 和 推导主动齿轮切片齿压缩弹簧 的变形量δpi和Fi的平衡方程,得到N个主动齿轮切片齿压缩弹簧的变形量和力的平衡方程:
[0108]
[0109] 步骤7b)将N个主动齿轮切片齿压缩弹簧的变形量的力与变形的平衡公式和力的平衡方程转换为矩阵形式KcX=F,并将其作为 的力与变形的平衡公式,其中,Kc表示刚度矩阵,X和F分别表示变形矢量和施加的力矢量,
[0110]
[0111] X=[δp1,δp2,…,δpi,…,δpN]T
[0112] F=[F1,F2,…,Fi,…,FN]T
[0113] [·]T表示矩阵[·]的转置;
[0114] 步骤8)计算齿轮副的时变啮合刚度:
[0115] 步骤8a)设施加在齿轮副中相互啮合的一对齿的切片耦合模型S上的外力为Fe,力容忍误差为Fε,并令δ=0;
[0116] 步骤8b)判断δ与切片齿 和 的总误差Ei的关系δ>Ei是否成立,若是,则切片齿与切片齿 啮合,将主动齿轮p中所有处于啮合状态的切片齿组合为集合否则,则切片齿 与切片齿 不啮合,将主动齿轮p中所有处
于不啮合状态的切片齿 组合为集合 其中, 为主动齿轮
p中第w个处于啮合状态的切片齿,W为处于啮合状态的切片齿的总个数, 为主动齿轮p中第s个处于不啮合状态的切片齿,S为处于不啮合状态的切片齿的总个数,S+W=N;
[0117] 步骤8c)当切片齿 与切片齿 啮合时,主动齿轮p中W个处于啮合状态的切片齿的变形量均为 根据距离集合Ap”区域最近的处于啮合状态的切片齿 的变形量计算处于不啮合状态的切片齿 的变形量δ”ps,并将计算得到的 的变形量δpi填补到X=[δp1,δp2,…,δpi,…,δpN]T的对应行,其中δ”ps的计算公式为:
[0118]
[0119] 其中,Γj(j+1)为第j个切片齿和第j+1个切片齿之间的变形传递系数;
[0120] 当切片齿啮合时代表其等效的弹簧处于接触状态,非接触状态下弹簧的变形可用接触状态下弹簧的变形来表示,变形传递系数表示相邻并行排布的主动齿轮切片齿压缩弹簧的变形量之间的关系,变形传递系数Γj(j+1)的计算公式为:
[0121]
[0122] 当下标大于N时,
[0123] 步骤8d)计算齿轮副的啮合力Fm,并判断|Fe-Fm|≤Fε是否成立,若是,得到齿轮副的啮合变形δ和齿对的载荷分布,并执行步骤(8e),否则,令 并执行步骤(8b);
[0124] 其中啮合力Fm为力矢量中各分力Fi之和,计算公式为:
[0125] F=KX
[0126]
[0127] 其中K为刚度矩阵,当单齿啮合时K=Kc,双齿啮合时
[0128] 步骤8e)根据齿轮副啮合变形δ计算齿轮副的啮合刚度k,并计算一个啮合周期内各啮合位置的啮合刚度,即齿轮副的时变啮合刚度;
[0129] 齿对啮合刚度ktt和齿轮副啮合刚度k,计算公式分别为:
[0130]
[0131]
[0132] 以下结合仿真实验,对本发明的技术效果作以说明:
[0133] 1.仿真条件和内容:
[0134] 使用Matlab软件,对本发明计算方法和修正的考虑鼓向修形的齿轮副啮合特性分析方法现有技术的直齿轮啮合刚度计算精度进行对比仿真,其结果如表3所示。
[0135] 2.仿真结果分析:
[0136] 从表3可以看出,本发明在考虑齿向修形和不对中情况下,直齿轮啮合刚度计算精度相对于现有技术有显著提高。
[0137] 表3仿真对比结果
[0138]计算方法 20μm齿向修形量下的精度 0.05°不对中角度下的精度
本发明计算方法 98.7% 98%
现有技术计算方法 82.9% 94.9%
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈