CN104408220A - 一种改进的轮齿加载接触分析方法 - Google Patents
一种改进的轮齿加载接触分析方法 Download PDFInfo
- Publication number
- CN104408220A CN104408220A CN201410523305.0A CN201410523305A CN104408220A CN 104408220 A CN104408220 A CN 104408220A CN 201410523305 A CN201410523305 A CN 201410523305A CN 104408220 A CN104408220 A CN 104408220A
- Authority
- CN
- China
- Prior art keywords
- tooth
- contact
- sigma
- flank
- formula
- 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
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/80—Technologies aiming to reduce greenhouse gasses emissions common to all road transportation technologies
- Y02T10/82—Elements for improving aerodynamics
Landscapes
- Gears, Cams (AREA)
Abstract
本发明公开了一种改进的轮齿加载接触分析方法。该方法基于变形协调条件,通过几何计算得到齿轮的齿面方程和法向量方程;并在此基础上进行无载荷的轮齿接触分析,得到齿轮啮合的接触印痕、传动误差,齿面的主曲率、主方向,相切齿面的诱导法曲率,从而建立起接触椭圆大小、齿面压缩量与载荷三者之间的关系;结合Hertz理论公式和柔度系数法,考虑安装误差、支撑误差的影响和接触柔度系数的非线性,求出齿轮的载荷分布;根据齿间载荷分配结果,采用Hertz公式计算小大齿轮某时刻下的齿面压力分布和接触应力分布,利用三维有限元模型进行受力计算得到可靠的弯曲应力。轮齿加载接触分析方法对齿轮学的设计、计算和应用有重要意义。
Description
技术领域
本发明属于齿轮结构设计分析领域,具体地说,涉及一种改进的轮齿加载接触分析方法。
背景技术
轮齿加载接触分析(LTCA)方法在国内外已引起了广泛的重视,但大量的研究还停留在轮齿接触分析的基础上,如何在齿轮研究领域中建立好几何设计与力学分析的桥梁,特别是获得齿轮在拟真实工况下的工作性能,对齿轮学的设计、计算和发展具有重要意义。
现有公开的技术文献“齿轮轮齿承载接触分析(LTCA)的模型和方法”(机械传动,1998,22(2),p1-3)中,基于柔度系数法和变形协调条件提出了一种轮齿加载接触分析(LTCA)方法,用于计算轮齿接触应力和弯曲应力,但该方法忽略了接触柔度系数的非线性影响,以及将接触区域载荷近似处理为线接触载荷来替代实际的区域载荷,当实际接触区域不是细长形区域,该方法的计算精度将变得更差,故不具有普遍性和较好的计算精度;其次,该方法将线接触载荷离散为一些点上的集中载荷,当离散点取得较多时,计算变量较多,计算复杂,对算法的收敛性要求高。
发明内容
为了避免传统TCA方法忽略轮齿承载接触变形和弯曲变形的影响,能够真实有效地反应轮齿在加载情况下的齿间载荷分布情况,齿面载荷分布情况,传动误差以及精确的计算轮齿接触应力和弯曲应力,本发明提出一种改进的轮齿加载接触分析方法。该方法基于变形协调条件,结合Hertz理论公式和柔度系数法,考虑安装误差、支撑误差的影响和接触柔度系数的非线性,求出齿轮的载荷分布;在此基础上,利用三维有限元模型进行受力计算得到可靠的弯曲应力,利用Hertz理论进行了接触应力的计算。
本发明解决其技术问题所采用的技术方案是:一种改进的轮齿加载接触分析方法,其特征在于包括以下步骤:
步骤1.推导及获得齿面方程、法向量;
根据微分几何与啮合原理,求出一对小大齿轮的齿面方程及法向量,其表示为:在与小齿轮固联的坐标系S1下,小齿轮相邻三个轮齿的齿面方程及法向量为r11(u1,v1),n11(u1,v1),r12(u1,v1),n12(u1,v1)和r13(u1,v1),n13(u1,v1);在与大齿轮固联的坐标系S2下,与小齿轮各轮齿相啮合的大齿轮齿面的齿面方程及法向量为r21(u2,v2),n21(u2,v2),r22(u2,v2),n22(u2,v2)和r23(u2,v2),n23(u2,v2);
步骤2.进行无载荷的轮齿接触分析;
①计算接触印痕和传动误差
以小齿轮三轮齿的中间齿为基准面,分别计算进行多对齿面的无载荷轮齿接触分析,得到在小齿轮取不同转角φ1时,多对齿面副的接触印痕和传动误差,其具体如下:
建立一个装配坐标系Sf,分别将小齿轮、大齿轮的某对啮合轮齿的齿面方程r1j(u1,v1),r2j(u2,v2)及其法向量n1j(u1,v1),n2j(u2,v2)转换到装配坐标系Sf下,建立轮齿接触分析方程组;
给定一系列的小齿轮转角φ1,通过上述方程组可求出另外五个未知量:u1,v1,u2,v2和大齿轮转角φ2;
将上述求得变量分别代入小大齿轮的齿面方程r1j(u1,v1),r2j(u2,v2),则可得到某对齿面上一系列的接触点,即接触路径;
将(1)式求解参数代入传动误差公式则可得到接触点处的传动误差,并绘制不同轮齿接触对的传动误差曲线;
②推导齿面的主曲率、主方向
已知,一个齿面方程为r(u,v),则单位法矢n=(ru×rv)/|ru×rv|,式中,
曲面r(u,v)上某点P(x,y)的法曲率表示如下:
式中,E,F,G为曲面的第一基本量,L,M,N为曲面的第二基本量,du/dv为点P(x,y)处的切线方向,其中,F=ru·rv,L=-n·ruu,M=n·ruv,N=n·rvv;
由(2)式知,kn随μ=du/dv而变化,其最大值和最小值成为点P(x,y)的主曲率,对应切线方向即为主方向;
kn是μ的函数,对μ求导得
当dkn/dμ=0时,kn取极值,有
(knE-L)μ+(knF-M)=0 (4)
联立式(2)与(4)得到
(knF-M)μ+(knG-N)=0 (5)
联立式(4)与(5),消除μ,则有
求解式(6)则可得到主曲率;
联立式(4)与(5),消除kn,则有
(EM-FL)μ2+(EN-GL)μ+(FN-GM)=0 (7)
求解式(7)则可得到主方向;
根据罗德里克定理,判断与主曲率对应的主方向,其关系式为:
根据式(7)求出的主方向参数μ,可确定主方向的单位矢量:
③相切齿面的诱导法曲率,求得接触椭圆的长短轴方向
根据步骤②求得的主曲率和主方向可得到相切齿面的诱导法曲率,求得接触椭圆的长短轴方向,其过程如下:
规定小齿轮齿面上的第一主方向到大齿轮齿面上的第一主方向的有向角为
已知曲面的主曲率及主方向,根据法曲率的欧拉公式表达出任意方向的法曲率
式中, i=1,2;
两曲面在切点沿同一方向的法曲率之差,即该方向的诱导法曲率,表示为
将σ(2)=σ(1)-σ(12)代入式(12),化简得
是σ(1)的函数,对式(13)求导,令得到诱导法曲率取极值时对应的σ(1),其化简后表示为:
式中,
由此,可以进一步得到接触椭圆的长短轴方向:
④建立单个接触位置上接触椭圆大小、齿面压缩量与载荷的关系
由上述步骤①、②、③,求得某时刻接触点位置、法向量、主曲率、主方向、接触椭圆长短轴方向,建立某时刻接触点位置上接触椭圆大小、齿面压缩量与载荷的关系,具体如下:
根据Hertz接触理论,可得到两轮齿压缩量δc与接触椭圆大小之间的关系
式中
同时,可得到两轮齿压缩量δc与载荷Pi的关系
式中,v1,v2,E1,E2分别为小大齿轮的泊松比和弹性模量;Re为等效半径,Re=(R′R″)1/2,1/R′=1/R1′+1/R2′,1/R″=1/R1″+1/R2″,Ri′和Ri″分别为小大齿轮在椭圆中心位置的主曲率半径;F2(e)函数可看作是对于椭圆偏心率的修正因子,其相对于相对曲率比值(R′/R″)的变化曲线已知;
接触椭圆大小、齿面压缩量与载荷的关系由式(16)、(17)确定;
步骤3.计算载荷分配;
通过步骤1,步骤2,得到某时刻小大齿轮副的多齿面接触情况,即某时刻各个齿面的接触位置、传动误差,同时建立单个接触位置上接触椭圆大小、齿面压缩量与载荷的关系,下面依据变形协调条件可求出各个接触位置上的载荷分量;
首先,根据轮齿无载荷接触分析得到的多对齿面的传动误差,可判断小齿轮在转动某转角φ1时,小大齿轮的接触点对数m;
设各接触齿面对上分担的载荷分别为P1,P2,……Pm,则根据变形协调条件和总载荷不变得到方程组:
式中,RGi为第i对齿轮接触点距转动中心的转动半径;αti表示载荷与切平面之间的夹角;Δθi为加载传动误差,T为小大齿轮传递总扭矩;
已知,加载传动误差Δθi由三部分组成,表示为
Δθi(Pi)=Δθio(Pi)+Δθic(Pi)+Δθib(Pi) (19)
式中,Δθio(Pi)为无加载传动误差,Δθic(Pi)为接触变形引起的加载传动误差,Δθib(Pi)为轮齿弯曲变形引起的加载传动误差;
下面对三种传动误差分量进行求解:
①求解Δθio(Pi),根据上述无载荷轮齿接触分析,可得到
②求解Δθic(Pi),上述按Hertz接触理论建立起两轮齿压缩量δic与载荷Pi的关系即δic=f(Pi),则有
③求解Δθib(Pi),首先,分别建立大小齿轮的三齿三维有限元网格模型,进行齿轮各个齿面的弯曲柔度系数矩阵提取,其方法如下:
对齿轮的三齿三维有限元模型添加两侧面的位移约束和底面位移约束,下面,分别给工作齿面上的各个节点施加1N的法向力,提取齿面上所有节点的位移,若某个工作齿面上有n个节点,则提取出的柔度系数矩阵为[Cf1]n×n×3;
对齿轮的三齿三维有限元模型添加两侧面的位移约束和底面位移约束,同时将各个非工作齿面添加位移约束,下面分别给工作齿面上的各个节点施加1N的法向力,提取齿面上所有节点的位移,若某个工作齿面上有n个节点,则提取出的柔度系数矩阵为[Cf2]n×n×3;
则齿轮的弯曲柔度系数矩阵为
[Cf]=[Cf1]-[Cf2] (22)
根据Hertz接触理论可得,一般外形曲面的表面压力为
将接触椭圆区域离散化,建立分割区域内载荷对接触点的弯曲弹性变形量的微分表达式
dδbi=p(ρ,θ)·cf(ρ,θ)·ρdρdθ (24)
由此,建立一个分割区域内载荷对接触点的弯曲弹性变形量的近似表达式
Δδbi=p(ρ,θ)·cf(ρ,θ)·ρ·Δρ·Δθ (25)
其中,ρ,θ是通过将接触椭圆在极坐标下的表示,且有
故ρ的积分区间为[0,1],θ的积分区间为[0,2π];cf(ρ,θ)为分割区域类的柔度系数,通过对区域范围内或周围的节点柔度系数插值而言,插值选取节点越多,cf(ρ,θ)越精确;
最后,对接触椭圆内所有分割区域内载荷对接触点的弯曲弹性变形量进行求和,则可得到一个齿面上由弯曲引起的接触点位置的位移量:
由弯曲变形引起的传动误差分量为
式中,δbi1,δbi2分别代表小大齿轮一对齿轮副在接触点位置的轮齿弯曲变形量;
联立式(18)~(28),进行求解,则可得到齿间载荷分布P1,P2,……Pm;
再次计算计入弯曲变形和接触变形后的传动误差,判断小齿轮在转动某转角φ1时,小大齿轮的接触点对数k,若k>m,则将k赋予m,重新进行步骤3,若k=m,则进行下一步;
步骤4.进行强度计算;
按照上述分析计算得到的齿间载荷分配结果,则可根据Hertz理论公式计算小大齿轮各对轮齿副在小齿轮转动某转角φ1时的齿面压力分布和接触应力分布;
根据上述建立小、大齿轮三维有限元模型,将计算好的齿面压力分布等效施加在接触区域内的节点上,采用现有有限元软件进行弯曲强度计算,则可得到轮齿接触产生的弯曲应力分布。
有益效果
本发明提出的一种改进的轮齿加载接触分析方法,通过几何计算得到齿轮的齿面方程和法向量方程;在此基础上进行无载荷的轮齿接触分析,得到齿轮啮合的接触印痕、传动误差,齿面的主曲率、主方向,相切齿面的诱导法曲率,从而建立接触椭圆大小、齿面压缩量与载荷三者之间的关系;通过上述得到的传动误差、啮合印痕、以及载荷与接触椭圆大小、齿面压缩量之间的关系,联系变形协调方程、Hertz方程,则可求解一对齿轮副的齿间载荷分布;最终,根据齿间载荷分配结果,采用Hertz公式计算小大齿轮某时刻下的齿面压力分布和接触应力分布,采用现有有限元软件计算得到轮齿接触产生的弯曲应力分布。该方法考虑了安装误差、支撑误差的影响和接触柔度系数的非线性,能够精确和简便的计算轮齿接触应力和弯曲应力,对齿轮学的设计、计算和发展有重要意义。
附图说明
下面结合附图和实施方式对本发明一种改进的轮齿加载接触分析方法作进一步详细说明。
图1为相切齿面的诱导法曲率示意图。
图2a为非工作齿面无位移约束的有限元计算模型示意图。
图2b为添加非工作齿面位移约束的有限元计算模型示意图。
图3为接触椭圆区域的分割示意图。
具体实施方式
本实施例轮齿加载接触分析方法是以一对三齿几何模型进行计算分析。
第一步,推导及获得齿面方程、法向量;
根据微分几何与啮合原理,求出一对小大齿轮的齿面方程,其表示方式为:
小齿轮中间啮合齿面在与小齿轮固联的坐标系S1下的齿面方程及法向量为r11(u1,v1)和n11(u1,v1),根据坐标旋转变换可得到小齿轮中间齿左右轮齿的齿面方程及其法向量为r12(u1,v1),n12(u1,v1)和r13(u1,v1),n13(u1,v1);
同理,在与大齿轮固联的坐标系S2下,与小齿轮各轮齿相啮合的大齿轮齿面的齿面方程及法向量为r21(u2,v2),n21(u2,v2)、r22(u2,v2),n22(u2,v2)和r23(u2,v2),n23(u2,v2);
第二步,进行无载荷的轮齿接触分析;
①计算接触印痕和传动误差
以小齿轮三轮齿的中间齿为基准轮齿,分别计算进行多对齿面的无载荷轮齿接触分析,得到在小齿轮取不同转角φ1时,多对齿面副的接触印痕和传动误差,其具体如下:
建立一个装配坐标系Sf,分别将小、大齿轮的某对啮合轮齿的齿面方程r1j(u1,v1),r2j(u2,v2)及其法向量n1j(u1,v1),n2j(u2,v2)转换到装配坐标系Sf下,建立轮齿接触分析方程组;
给定一系列的小齿轮转角φ1,通过上述方程组可求出另外五个未知量:u1,v1,u2,v2和大齿轮转角φ2;
将上述求得变量分别代入小大齿轮的齿面方程r1j(u1,v1),r2j(u2,v2),则可得到某对齿面上一系列的接触点,即接触路径;将上述求得变量分别代入小大齿轮的法向量方程n1j(u1,v1),n2j(u2,v2),则可得到某对齿面上对应接触点位置的法向量;
通过传动误差公式则可得到接触点处的传动误差,并绘制不同轮齿接触对的传动误差曲线;
②推导齿面的主曲率、主方向
已知,一个齿面方程表示为r(u,v),则单位法矢n=(ru×rv)/|ru×rv|,式中,
曲面r(u,v)上某点P(x,y)的法曲率表示如下:
式中,E,F,G为曲面的第一基本量,L,M,N为曲面的第二基本量,du/dv为点P(x,y)处的切线方向;其中,F=ru·rv,L=-n·ruu,M=n·ruv,N=n·rvv;
由(2)式知,kn随μ=du/dv而变化,其最大值和最小值成为点P(x,y)的主曲率,对应切线方向即为主方向;
kn是μ的函数,对μ求导得
当dkn/dμ=0时,kn取极值,有
(knE-L)μ+(knF-M)=0 (4)
联立式(2)与(4)得到
(knF-M)μ+(knG-N)=0 (5)
联立式(4)与(5),消除μ,则有
求解式(6)则可得到主曲率;
联立式(4)与(5),消除kn,则有
(EM-FL)μ2+(EN-GL)μ+(FN-GM)=0 (7)
求解式(7)则可得到主方向;
根据罗德里克定理,判断与主曲率对应的主方向,其关系式为:
根据式(7)求出的主方向参数μ,可确定主方向的单位矢量:
下面,用分别代表主方向的单位矢量、主曲率,上标(i)=(1)、(2)分别代表小、大齿轮,下标j=I、II分别代表第一主曲率、第二主曲率及其对应的主方向,下面计算并规定小齿轮的第一主曲率为法曲率的最大值,大齿轮的第一主曲率为法曲率的最小值;
③相切齿面的诱导法曲率,求得接触椭圆的长短轴方向
根据步骤②求得的主曲率和主方向可得到相切齿面的诱导法曲率,求得接触椭圆的长短轴方向,其过程如下:
参阅图1,Σ1、Σ2为空间相切的小、大齿面,切点为P,在切点处两齿面的单位法矢n的方向相同,σ(i)分别为主方向到任意方向的有向角;
规定小齿轮齿面上的第一主方向到大齿轮齿面上的第一主方向的有向角为
已知曲面的主曲率及主方向,根据法曲率的欧拉公式表达出任意方向的法曲率
式中, i=1,2;
两曲面在切点沿同一方向的法曲率之差称为该方向的诱导法曲率,表示为
将σ(2)=σ(1)-σ(12)代入式(12),化简得
是σ(1)的函数,对式(13)求导,令得到诱导法曲率取极值时对应的σ(1),其化简后表示为:
式中,
由此,可进一步得到接触椭圆的长短轴方向:
④建立单个接触位置上接触椭圆大小、齿面压缩量与载荷的关系
由上述步骤①、②、③,求得某时刻接触点位置、法向量、主曲率、主方向、接触椭圆长短轴方向,建立某时刻接触点位置上接触椭圆大小、齿面压缩量与载荷的关系,具体如下:
根据Hertz接触理论,可得到两轮齿压缩量δc与接触椭圆大小之间的关系
式中
同时,可得到两轮齿压缩量δc与载荷Pi的关系
式中,v1,v2,E1,E2分别为小大齿轮的泊松比和弹性模量;Re为等效半径,Re=(R′R″)1/2,1/R′=1/R1′+1/R2′,1/R″=1/R1″+1/R2″,Ri′和Ri″分别为小大齿轮在椭圆中心位置的主曲率半径;F2(e)函数可作为是对于椭圆偏心率的修正因子,其相对于相对曲率比值(R′/R″)的变化曲线已知;
至此,接触椭圆大小、齿面压缩量与载荷的关系由式(16)、(17)确定;
第三步,计算载荷分配;
通过第一步,第二步,得到某时刻小大齿轮副的多齿面接触情况,即某时刻各个齿面的接触位置、传动误差,同时建立单个接触位置上接触椭圆大小、齿面压缩量与载荷的关系,下面依据变形协调条件则可求出各个接触位置上的载荷分量;
首先,根据上述轮齿无载荷接触分析得到多对齿面的传动误差,可判断小齿轮在转动某转角φ1时,小大齿轮的接触点对数m;
假设,各接触齿面对上分担的载荷分别为P1,P2,……Pm,则根据变形协调条件和总载荷不变得到方程组:
式中,RGi为第i对齿轮接触点距转动中心的转动半径,αti为载荷与切平面之间的夹角,Δθi为加载传动误差,T为小大齿轮传递总扭矩;
已知,加载传动误差Δθi由三部分组成,表示为
Δθi(Pi)=Δθio(Pi)+Δθic(Pi)+Δθib(Pi) (19)
式中,Δθio(Pi)为无加载传动误差,Δθic(Pi)为接触变形引起的加载传动误差,Δθib(Pi)为轮齿弯曲变形引起的加载传动误差;
下面对这三种传动误差分量进行求解:
①求解Δθio(Pi),根据上述无载荷轮齿接触分析,可得到
②求解Δθic(Pi),上述按Hertz接触理论建立起两轮齿压缩量δic与载荷Pi的关系,即δic=f(Pi),则有
③求解Δθib(Pi),首先,分别建立大小齿轮的三齿三维有限元网格模型,进行齿轮各个齿面的弯曲柔度系数矩阵提取,其方法如下:
参阅图2a,对齿轮的三齿三维有限元模型添加两侧面的位移约束和底面位移约束,分别给工作齿面上的各个节点施加1N的法向力,提取齿面上所有节点的位移,若某个工作齿面上有n个节点,则提取出的柔度系数矩阵为[Cf1]n×n×3;
参阅图2b,对齿轮的三齿三维有限元模型添加两侧面的位移约束和底面位移约束,同时将各个非工作齿面添加位移约束,分别给工作齿面上的各个节点施加1N的法向力,提取齿面上所有节点的位移,若某个工作齿面上有n个节点,则提取出的柔度系数矩阵为[Cf2]n×n×3;
则齿轮的弯曲柔度系数矩阵为
[Cf]=[Cf1]-[Cf2] (22)
根据Hertz接触理论可得,一般外形曲面的表面压力为
参阅图3,将接触椭圆区域离散化,建立分割区域内载荷对接触点的弯曲弹性变形量的微分表达式
dδbi=p(ρ,θ)·cf(ρ,θ)·ρdρdθ (24)
由此,建立一个分割区域内载荷对接触点的弯曲弹性变形量的近似表达式
Δδbi=p(ρ,θ)·cf(ρ,θ)·ρ·Δρ·Δθ (25)
其中,ρ,θ是通过将接触椭圆在极坐标下的表示,且有
故ρ的积分区间为[0,1],θ的积分区间为[0,2π];cf(ρ,θ)为分割区域类的柔度系数,通过对区域范围内或周围的节点柔度系数插值而言,插值选取节点越多,cf(ρ,θ)越精确;
最后,对接触椭圆内所有分割区域内载荷对接触点的弯曲弹性变形量进行求和,则可得到一个齿面上由弯曲引起的接触点位置的位移量:
由弯曲变形引起的传动误差分量为
式中,δbi1,δbi2分别代表小大齿轮一对齿轮副在接触点位置的轮齿弯曲变形量;
联立式(18)~(28),进行求解,则可得到齿间载荷分布P1,P2,……Pm;
再次计算计入弯曲变形和接触变形后的传动误差,判断小齿轮在转动某转角φ1时,小大齿轮的接触点对数k,若k>m,则将k赋予m,重新进行第三步,若k=m,则进行下一步;
第四步,进行强度计算;
按照上述分析计算得到的齿间载荷分配结果,则根据Hertz理论公式计算小大齿轮各对轮齿副在小齿轮转动某转角φ1时的齿面压力分布和接触应力分布;
根据上述建立小、大齿轮三维有限元模型,将计算好的齿面压力分布等效施加在接触区域内或者含接触区域附近的节点上,采用现有有限元软件ANSYS或ABAQUS进行弯曲强度计算,则可分别得到小大轮齿承载接触产生的弯曲应力分布。
本实施例所计算轮齿弯曲应力和接触应力的方法,考虑安装误差、支撑误差的影响和接触柔度系数的非线性,同时又尽量保证接触区域的几何形状,能较好的获得齿轮在拟真实工况下的工作性能。
Claims (1)
1.一种改进的轮齿加载接触分析方法,其特征在于包括以下步骤:
步骤1.推导及获得齿面方程、法向量;
根据微分几何与啮合原理,求出一对小大齿轮的齿面方程及法向量,其表示为:在与小齿轮固联的坐标系S1下,小齿轮相邻三个轮齿的齿面方程及法向量为r11(u1,v1),n11(u1,v1),r12(u1,v1),n12(u1,v1)和r13(u1,v1),n13(u1,v1);在与大齿轮固联的坐标系S2下,与小齿轮各轮齿相啮合的大齿轮齿面的齿面方程及法向量为r21(u2,v2),n21(u2,v2),r22(u2,v2),n22(u2,v2)和r23(u2,v2),n23(u2,v2);
步骤2.进行无载荷的轮齿接触分析;
①计算接触印痕和传动误差
以小齿轮三轮齿的中间齿为基准面,分别计算进行多对齿面的无载荷轮齿接触分析,得到在小齿轮取不同转角φ1时,多对齿面副的接触印痕和传动误差,其具体如下:
建立一个装配坐标系Sf,分别将小齿轮、大齿轮的某对啮合轮齿的齿面方程r1j(u1,v1),r2j(u2,v2)及其法向量n1j(u1,v1),n2j(u2,v2)转换到装配坐标系Sf下,建立轮齿接触分析方程组;
给定一系列的小齿轮转角φ1,通过上述方程组可求出另外五个未知量:u1,v1,u2,v2和大齿轮转角φ2;
将上述求得变量分别代入小大齿轮的齿面方程r1j(u1,v1),r2j(u2,v2),则可得到某对齿面上一系列的接触点,即接触路径;
将(1)式求解参数代入传动误差公式则可得到接触点处的传动误差,并绘制不同轮齿接触对的传动误差曲线;
②推导齿面的主曲率、主方向
已知,一个齿面方程为r(u,v),则单位法矢n=(ru×rv)/|ru×rv|,式中,
曲面r(u,v)上某点P(x,y)的法曲率表示如下:
式中,E,F,G为曲面的第一基本量,L,M,N为曲面的第二基本量,du/dv为点P(x,y)处的切线方向,其中,F=ru·rv,L=-n·ruu,M=n·ruv,N=n·rvv;
由(2)式知,kn随μ=du/dv而变化,其最大值和最小值成为点P(x,y)的主曲率,对应切线方向即为主方向;
kn是μ的函数,对μ求导得
当dkn/dμ=0时,kn取极值,有
(knE-L)μ+(knF-M)=0 (4)
联立式(2)与(4)得到
(knF-M)μ+(knG-N)=0 (5)
联立式(4)与(5),消除μ,则有
求解式(6)则可得到主曲率;
联立式(4)与(5),消除kn,则有
(EM-FL)μ2+(EN-GL)μ+(FN-GM)=0 (7)
求解式(7)则可得到主方向;
根据罗德里克定理,判断与主曲率对应的主方向,其关系式为:
根据式(7)求出的主方向参数μ,可确定主方向的单位矢量:
③相切齿面的诱导法曲率,求得接触椭圆的长短轴方向
根据步骤②求得的主曲率和主方向可得到相切齿面的诱导法曲率,求得接触椭圆的长短轴方向,其过程如下:
规定小齿轮齿面上的第一主方向到大齿轮齿面上的第一主方向的有向角为
已知曲面的主曲率及主方向,根据法曲率的欧拉公式表达出任意方向的法曲率
式中, i=1,2;
两曲面在切点沿同一方向的法曲率之差,即该方向的诱导法曲率,表示为
将σ(2)=σ(1)-σ(12)代入式(12),化简得
是σ(1)的函数,对式(13)求导,令得到诱导法曲率取极值时对应的σ(1),其化简后表示为:
式中,
由此,可以进一步得到接触椭圆的长短轴方向:
④建立单个接触位置上接触椭圆大小、齿面压缩量与载荷的关系
由上述步骤①、②、③,求得某时刻接触点位置、法向量、主曲率、主方向、接触椭圆长短轴方向,建立某时刻接触点位置上接触椭圆大小、齿面压缩量与载荷的关系,具体如下:
根据Hertz接触理论,可得到两轮齿压缩量δc与接触椭圆大小之间的关系
式中
同时,可得到两轮齿压缩量δc与载荷Pi的关系
式中,v1,v2,E1,E2分别为小大齿轮的泊松比和弹性模量;Re为等效半径,Re=(R′R″)1/2,1/R′=1/R1′+1/R2′,1/R″=1/R1″+1/R2″,Ri′和Ri″分别为小大齿轮在椭圆中心位置的主曲率半径;F2(e)函数可看作是对于椭圆偏心率的修正因子,其相对于相对曲率比值(R′R″)的变化曲线已知;
接触椭圆大小、齿面压缩量与载荷的关系由式(16)、(17)确定;
步骤3.计算载荷分配;
通过步骤1,步骤2,得到某时刻小大齿轮副的多齿面接触情况,即某时刻各个齿面的接触位置、传动误差,同时建立单个接触位置上接触椭圆大小、齿面压缩量与载荷的关系,下面依据变形协调条件可求出各个接触位置上的载荷分量;
首先,根据轮齿无载荷接触分析得到的多对齿面的传动误差,可判断小齿轮在转动某转角φ1时,小大齿轮的接触点对数m;
设各接触齿面对上分担的载荷分别为P1,P2,……Pm,则根据变形协调条件和总载荷不变得到方程组:
式中,RGi为第i对齿轮接触点距转动中心的转动半径;αti表示载荷与切平面之间的夹角;△θi为加载传动误差,T为小大齿轮传递总扭矩;
已知,加载传动误差△θi由三部分组成,表示为
△θi(Pi)=△θio(Pi)+△θic(Pi)+△θib(Pi) (19)
式中,△θio(Pi)为无加载传动误差,△θic(Pi)为接触变形引起的加载传动误差,△θib(Pi)为轮齿弯曲变形引起的加载传动误差;
下面对三种传动误差分量进行求解:
①求解△θio(Pi),根据上述无载荷轮齿接触分析,可得到
②求解△θic(Pi),上述按Hertz接触理论建立起两轮齿压缩量δic与载荷Pi的关系即δic=f(Pi),则有
③求解△θib(Pi),首先,分别建立大小齿轮的三齿三维有限元网格模型,进行齿轮各个齿面的弯曲柔度系数矩阵提取,其方法如下:
对齿轮的三齿三维有限元模型添加两侧面的位移约束和底面位移约束,下面,分别给工作齿面上的各个节点施加1N的法向力,提取齿面上所有节点的位移,若某个工作齿面上有n个节点,则提取出的柔度系数矩阵为[Cf1]n×n×3;
对齿轮的三齿三维有限元模型添加两侧面的位移约束和底面位移约束,同时将各个非工作齿面添加位移约束,下面分别给工作齿面上的各个节点施加1N的法向力,提取齿面上所有节点的位移,若某个工作齿面上有n个节点,则提取出的柔度系数矩阵为[Cf2]n×n×3;
则齿轮的弯曲柔度系数矩阵为
[Cf]=[Cf1]-[Cf2] (22)
根据Hertz接触理论可得,一般外形曲面的表面压力为
将接触椭圆区域离散化,建立分割区域内载荷对接触点的弯曲弹性变形量的微分表达式
dδbi=p(ρ,θ)·cf(ρ,θ)·ρdρdθ (24)
由此,建立一个分割区域内载荷对接触点的弯曲弹性变形量的近似表达式
△δbi=p(ρ,θ)·cf(ρ,θ)·ρ·△ρ·△θ (25)
其中,ρ,θ是通过将接触椭圆在极坐标下的表示,且有
故ρ的积分区间为[0,1],θ的积分区间为[0,2π];cf(ρ,θ)为分割区域类的柔度系数,通过对区域范围内或周围的节点柔度系数插值而言,插值选取节点越多,cf(ρ,θ)越精确;
最后,对接触椭圆内所有分割区域内载荷对接触点的弯曲弹性变形量进行求和,则可得到一个齿面上由弯曲引起的接触点位置的位移量:
由弯曲变形引起的传动误差分量为
式中,δbi1,δbi2分别代表小大齿轮一对齿轮副在接触点位置的轮齿弯曲变形量;
联立式(18)~(28),进行求解,则可得到齿间载荷分布P1,P2,……Pm;
再次计算计入弯曲变形和接触变形后的传动误差,判断小齿轮在转动某转角φ1时,小大齿轮的接触点对数k,若k>m,则将k赋予m,重新进行步骤3,若k=m,则进行下一步;
步骤4.进行强度计算;
按照上述分析计算得到的齿间载荷分配结果,则可根据Hertz理论公式计算小大齿轮各对轮齿副在小齿轮转动某转角φ1时的齿面压力分布和接触应力分布;
根据上述建立小、大齿轮三维有限元模型,将计算好的齿面压力分布等效施加在接触区域内的节点上,采用现有有限元软件进行弯曲强度计算,则可得到轮齿接触产生的弯曲应力分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410523305.0A CN104408220B (zh) | 2014-10-08 | 2014-10-08 | 一种改进的轮齿加载接触分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410523305.0A CN104408220B (zh) | 2014-10-08 | 2014-10-08 | 一种改进的轮齿加载接触分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104408220A true CN104408220A (zh) | 2015-03-11 |
CN104408220B CN104408220B (zh) | 2017-08-25 |
Family
ID=52645851
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410523305.0A Active CN104408220B (zh) | 2014-10-08 | 2014-10-08 | 一种改进的轮齿加载接触分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104408220B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105787149A (zh) * | 2015-12-27 | 2016-07-20 | 北京航空航天大学 | 一种由弧齿锥齿轮传动***轴上功率谱向齿面应力谱精确转换的方法 |
CN108115217A (zh) * | 2017-12-13 | 2018-06-05 | 长安大学 | 一种基于高阶传动误差的螺旋锥齿轮加工方法 |
CN108253094A (zh) * | 2018-01-17 | 2018-07-06 | 河南科技大学 | 一种rv减速器摆线针轮承载啮合印痕确定方法及装置 |
CN109190324A (zh) * | 2018-11-12 | 2019-01-11 | 中南大学 | 一种弧齿锥齿轮含误差的齿面载荷接触变形数值计算方法 |
CN109376456A (zh) * | 2018-11-12 | 2019-02-22 | 中南大学 | 弧齿锥齿轮含安装误差的齿面载荷接触性能数值计算方法 |
CN109376455A (zh) * | 2018-11-12 | 2019-02-22 | 中南大学 | 一种弧齿锥齿轮载荷弹性接触变形数值计算方法 |
CN109446711A (zh) * | 2018-11-12 | 2019-03-08 | 中南大学 | 弧齿锥齿轮含安装误差的齿面载荷传动误差数值计算方法 |
CN109446710A (zh) * | 2018-11-12 | 2019-03-08 | 中南大学 | 一种弧齿锥齿轮载荷传动误差数值计算方法 |
CN109492307A (zh) * | 2018-11-12 | 2019-03-19 | 中南大学 | 一种弧齿锥齿轮齿面载荷接触性能参数的数值计算方法 |
CN110457820A (zh) * | 2019-08-12 | 2019-11-15 | 重庆大学 | 一种确定齿轮副接触椭圆的方法 |
CN110807281A (zh) * | 2019-10-30 | 2020-02-18 | 重庆大学 | 一种考虑材料夹杂物和残余应力的齿轮psn曲线绘制方法 |
CN113408082A (zh) * | 2021-08-20 | 2021-09-17 | 宁波东力传动设备有限公司 | 一种外啮合直齿圆柱齿轮副动态啮合力的计算方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107577876B (zh) * | 2017-09-07 | 2020-06-23 | 清华大学 | 一种螺旋锥齿轮齿面加载性能多目标优化方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101476614A (zh) * | 2009-01-20 | 2009-07-08 | 河南科技大学 | 直线-渐开线齿廓内啮合圆柱齿轮副 |
US20130085722A1 (en) * | 2011-09-27 | 2013-04-04 | Romax Technology Limited | Rotating machines |
CN104021303A (zh) * | 2014-06-19 | 2014-09-03 | 中国北方车辆研究所 | 一种基于实测载荷谱的弧齿锥齿轮修形齿面的修形方法 |
-
2014
- 2014-10-08 CN CN201410523305.0A patent/CN104408220B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101476614A (zh) * | 2009-01-20 | 2009-07-08 | 河南科技大学 | 直线-渐开线齿廓内啮合圆柱齿轮副 |
US20130085722A1 (en) * | 2011-09-27 | 2013-04-04 | Romax Technology Limited | Rotating machines |
CN104021303A (zh) * | 2014-06-19 | 2014-09-03 | 中国北方车辆研究所 | 一种基于实测载荷谱的弧齿锥齿轮修形齿面的修形方法 |
Non-Patent Citations (3)
Title |
---|
屈文涛: "双圆弧齿轮的静态和动态接触模型", 《机械设计》 * |
方宗德: "齿轮轮齿承载接触分析_LTCA_的模型和方法", 《机械传动》 * |
李永祥: "基于ANSYS的直齿面齿轮的承载接触分析", 《机械科学与技术》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105787149A (zh) * | 2015-12-27 | 2016-07-20 | 北京航空航天大学 | 一种由弧齿锥齿轮传动***轴上功率谱向齿面应力谱精确转换的方法 |
CN105787149B (zh) * | 2015-12-27 | 2019-03-08 | 北京航空航天大学 | 一种由弧齿锥齿轮传动***轴上功率谱向齿面应力谱精确转换的方法 |
CN108115217A (zh) * | 2017-12-13 | 2018-06-05 | 长安大学 | 一种基于高阶传动误差的螺旋锥齿轮加工方法 |
CN108253094A (zh) * | 2018-01-17 | 2018-07-06 | 河南科技大学 | 一种rv减速器摆线针轮承载啮合印痕确定方法及装置 |
CN108253094B (zh) * | 2018-01-17 | 2019-08-20 | 河南科技大学 | 一种rv减速器摆线针轮承载啮合印痕确定方法及装置 |
CN109376455A (zh) * | 2018-11-12 | 2019-02-22 | 中南大学 | 一种弧齿锥齿轮载荷弹性接触变形数值计算方法 |
CN109376455B (zh) * | 2018-11-12 | 2023-04-18 | 中南大学 | 一种弧齿锥齿轮载荷弹性接触变形数值计算方法 |
CN109446711A (zh) * | 2018-11-12 | 2019-03-08 | 中南大学 | 弧齿锥齿轮含安装误差的齿面载荷传动误差数值计算方法 |
CN109446710A (zh) * | 2018-11-12 | 2019-03-08 | 中南大学 | 一种弧齿锥齿轮载荷传动误差数值计算方法 |
CN109492307A (zh) * | 2018-11-12 | 2019-03-19 | 中南大学 | 一种弧齿锥齿轮齿面载荷接触性能参数的数值计算方法 |
CN109190324A (zh) * | 2018-11-12 | 2019-01-11 | 中南大学 | 一种弧齿锥齿轮含误差的齿面载荷接触变形数值计算方法 |
CN109376456B (zh) * | 2018-11-12 | 2023-04-18 | 中南大学 | 弧齿锥齿轮含安装误差的齿面载荷接触性能数值计算方法 |
CN109446711B (zh) * | 2018-11-12 | 2023-04-18 | 中南大学 | 弧齿锥齿轮含安装误差的齿面载荷传动误差数值计算方法 |
CN109376456A (zh) * | 2018-11-12 | 2019-02-22 | 中南大学 | 弧齿锥齿轮含安装误差的齿面载荷接触性能数值计算方法 |
CN109446710B (zh) * | 2018-11-12 | 2023-04-18 | 中南大学 | 一种弧齿锥齿轮载荷传动误差数值计算方法 |
CN109492307B (zh) * | 2018-11-12 | 2023-04-07 | 中南大学 | 一种弧齿锥齿轮齿面载荷接触性能参数的数值计算方法 |
CN109190324B (zh) * | 2018-11-12 | 2023-04-18 | 中南大学 | 一种弧齿锥齿轮含误差的齿面载荷接触变形数值计算方法 |
CN110457820B (zh) * | 2019-08-12 | 2022-11-18 | 重庆大学 | 一种确定齿轮副接触椭圆的方法 |
CN110457820A (zh) * | 2019-08-12 | 2019-11-15 | 重庆大学 | 一种确定齿轮副接触椭圆的方法 |
CN110807281A (zh) * | 2019-10-30 | 2020-02-18 | 重庆大学 | 一种考虑材料夹杂物和残余应力的齿轮psn曲线绘制方法 |
CN110807281B (zh) * | 2019-10-30 | 2024-01-26 | 重庆大学 | 一种考虑材料夹杂物和残余应力的齿轮psn曲线绘制方法 |
CN113408082A (zh) * | 2021-08-20 | 2021-09-17 | 宁波东力传动设备有限公司 | 一种外啮合直齿圆柱齿轮副动态啮合力的计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104408220B (zh) | 2017-08-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104408220A (zh) | 一种改进的轮齿加载接触分析方法 | |
CN106844818A (zh) | 基于粗糙表面的直齿轮三维接触刚度计算方法 | |
CN108052760A (zh) | 一种齿轮副非线性动力学计算方法 | |
CN101770538B (zh) | 含损伤性单齿故障圆柱直齿轮啮合刚度仿真分析方法 | |
CN104462637B (zh) | 齿轮齿条式变比转向器齿轮副齿条齿廓的设计方法 | |
CN103615501B (zh) | 一种小型谐波减速器及其优化设计方法 | |
Liu et al. | Transmission principle and geometrical model of eccentric face gear | |
CN106570255A (zh) | 一种基于行人保护的负泊松比吸能结构的优化方法 | |
CN107696048A (zh) | 一种具有力位传感的刚柔耦合机械手指及其感知方法 | |
CN104281730B (zh) | 一种大转动变形的板壳结构动响应的有限元分析方法 | |
CN109726520A (zh) | 考虑复杂基体与裂纹扩展路径的直齿轮啮合刚度计算方法 | |
CN105972184A (zh) | 一种基于速度瞬心法的摆线轮齿廓方程设计方法 | |
CN109885976B (zh) | 一种旋转叶片位移场反演重构方法及其*** | |
CN103235846A (zh) | 柔性线缆的实时装配仿真方法及装置 | |
CN107763173A (zh) | 一种基于有限元分析的斜齿轮时变啮合刚度计算方法 | |
CN104091065A (zh) | 一种求解浅水问题模拟间断水流数值的方法 | |
CN102645200B (zh) | 基于双面啮合的齿轮齿距偏差测量方法 | |
CN202659811U (zh) | 一种具有傅里叶函数节曲线的非圆齿轮副 | |
CN105787149A (zh) | 一种由弧齿锥齿轮传动***轴上功率谱向齿面应力谱精确转换的方法 | |
CN104615878A (zh) | 一种计算带有边缘接触的齿间载荷分配方法 | |
Li et al. | Analytical calculation of mesh stiffness for spiral bevel gears with an improved global tooth deformation model | |
CN105427375A (zh) | 一种三维路面轮廓生成的方法及*** | |
CN104142424B (zh) | 一种准确测量非正弦电压信号的方法及其*** | |
CN110020463A (zh) | 一种计算旋转柔性曲梁动力学响应的仿真方法 | |
CN103411546B (zh) | 钢结构三维精度的检验方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |