具体实施方式
下面结合附图对本发明方法的实施方式做详细说明。
本发明一种基于圆弧啮合线的大重合度内啮合齿轮齿形设计方法的流程图如图1所示,包括(1)利用非均匀B样条构造内啮合齿轮啮合线;(2)满足B样条啮合线的共轭齿轮齿廓曲线与啮合线的对应关系建立及其啮合性质分析;(3)构造满足三次非均匀B样条啮合线的内啮合齿轮共轭齿廓曲线;(4)以齿廓曲线光顺性要求及共轭齿轮啮合性能最优建立目标函数,对非均匀B样条啮合线进行优化;(5)内啮合齿轮共轭齿廓曲线的干涉检查;(6)构造内外齿轮齿根齿廓到齿根圆的过渡曲线,完成内啮合共轭齿轮整体设计。
该方法具体步骤如下:
步骤一、利用非均匀B样条构造内啮合齿轮啮合线。
在齿轮传动过程中,随着接触点在啮合线上滑动时,在齿轮坐标系上划出的轨迹就是参与啮合的齿轮齿廓曲线。内啮合齿轮的啮合线为两条,分别对应着内齿轮和外齿轮的左右齿廓,一般而言,两条啮合线是对称的。内啮合齿轮的啮合线分布在内齿轮和外齿轮齿顶圆之间,如图2所示,齿轮1和齿轮2是一对内啮合齿轮,齿轮1为外齿轮,齿轮2为内齿轮,设L1和L2为啮合线的两段,L1位于第二象限,L2位于第四象限,为保证啮合线过原点且在原点连续,分段啮合线L1和L2相交于原点,且在原点相切。
对于渐开线内啮合齿轮,其啮合线L1和L2为两条线段,相交于节点Of,且与X轴的夹角为压力角α,它们的端点分别为与内齿轮及中心轮齿顶圆的交点,对于摆线齿内啮合齿轮,其啮合线L1和L2为两段相切圆弧,相切与节点Of,其中L1是齿轮1右齿廓齿根与齿轮2左齿廓齿顶的啮合线,L2为齿轮1右齿廓齿顶与齿轮2左齿廓齿根的啮合线;对于一般情况,内啮合齿轮齿形不仅限于渐开线齿形和摆线齿形,只要满足共轭条件且不发生干涉均可成为内啮合齿轮齿形。B样条曲线具有灵活性和自由性,理论上,B样条曲线可以拟合和逼近任意的几何形状,因此可以利用B样条构造内啮合齿轮啮合线,从而构造满足B样条啮合线的内啮合齿轮齿形曲线。k次B样条曲线方程为:
p(u)为分段多项式,其中,di(i=0,1,2,…,n)为控制顶点,Ni,k(u)(i=0,1,2,…,n)为k次规范B样条基函数,简称B样条,它是由节点矢量U=[u0,u1,…,un+k+1]:u0≤u1≤…≤un+k+1所决定的k次多项式。当节点矢量在可能范围内(两端节点重复度小于等于k+1,内节点重复度小于等于k)任意选取,构造的B样条曲线为非均匀B样条曲线,确定节点矢量的方法有里森费尔德方法和哈特利-贾德方法。
当k=3时,三次B样条曲线可以表示为:
其中,样条基函数Ni,3(u)(i=0,1,2,…,n)为三次多项式,根据B样条基函数的定义可以得到:
[Ni-3,3(u) Ni-2,3(u) Ni-1,3(u) Ni,3(u)]=[u3 u2 u 1]Mi (3)
其中,Mi是由多项式的系数矩阵,把式(3)带入式(2)得到三次B样条曲线又可以表示为:
对于非均匀三次B样条曲线,其在节点p(ui)(i=k+1,5,..,n)处为C2阶连续,在其他内点为C∞阶连续,并且非均匀B样条曲线在首末端节点p(uk)、p(un+1)分别与控制多边形首末边相切。
利用非均匀B样条曲线构造分段啮合线L1、L2,他们的控制顶点分别为:di (L1)(i=0,1,2,…,n)、dij (L2)(j=0,1,2,…,m),为保证分段啮合线L1和L2相交于原点,且在原点相切,根据B样条曲线的端点性质,使控制顶点d0 (L1)=d0 (L2)=(0,0),d2 (L1)、d1 (L1)、d0 (L1)、d1 (L2)、d2 (L2)在同一条直线l0上。dn (L1)、dm L2)分别在齿轮1和齿轮2的齿顶圆上,限制了啮合线的啮合范围。
步骤二、满足B样条啮合线的共轭齿轮齿廓曲线与啮合线的对应关系建立及其啮合性质分析
如图3,齿轮1和齿轮2的节圆内切于点P0,设直角坐标系S1:{O1,X1,Y1},S2:{O2,X2,Y2}分别与齿轮1和齿轮2固连,坐标系Sf:{Of,Xf,Yf}为固定坐标系,且Of为节点,则Of在齿轮啮合作用线上。设分别为齿轮1和齿轮2传动时的角位移,则齿轮1和齿轮2的传动内啮合齿轮的中心距E12=R2-R1,其中R1、R2设分别为齿轮1和齿轮2的节圆半径。B样条啮合线在直角坐标系Sf上的齐次坐标表达式为:
其中,xf,yf为B样条啮合线上的点P在直角坐标系Sf上的坐标值。在齿轮啮合过程中,接触点在啮合线上并沿着啮合线滑动,接触点在直角坐标系S1和S2上滑动的轨迹就是齿轮1和齿轮2的齿廓曲线。啮合线上的共轭点在直角坐标系S1上的轨迹的向量表示为:
其中,xf,yf为啮合点P在直角坐标系Sf上的坐标值。在齿轮啮合过程中,接触点在啮合线上并沿着啮合线滑动,接触点在直角坐标系S1和S2上滑动的轨迹就是齿轮1和齿轮2的齿廓曲线。啮合线上的共轭点在直角坐标系S1上的轨迹的向量表示为:
其中,M1f为直角坐标系Sf到S1的坐标变换矩阵。
啮合线上的共轭点在直角坐标系S2上的轨迹的向量表示为:
其中,M1f为直角坐标系Sf到S1的坐标变换矩阵。
根据平面齿轮啮合理论,共轭齿形曲线在切触点处的公法线与回转中心线O1O2相交并将该线分为两段O2P和O1P,两线段之比值为齿轮传动比。得到的齿轮1的啮合方程为:
对式(10)积分得到:
得到齿轮2的角位移的表达式为:
式(11)、(12)建立起了啮合线上的点与共轭齿廓曲线上接触点的对应关系,把式(11)、(12)带入式(6)、(8)得到齿轮1和齿轮2的共轭齿廓曲线表达式。
在得到内啮合齿轮共轭齿廓曲线之前,根据啮合线与共轭齿廓曲线的对应关系(11)和(12)可以分析构造的内啮合齿轮共轭齿廓曲线的啮合性质,如相对法曲率、相对滑动率。
根据曲线的法曲率关系式得到齿轮1齿廓曲线在点P处的法曲率为:
同样,齿轮2齿廓曲线在点P处的法曲率为:
齿轮1和齿轮2共轭齿廓曲线在啮合点P处的相对法曲率为:
κ12=κ1-κ2 (14)
根据齿轮的齿面方程(6)和(8)得到齿轮1的弧长元素为:
齿轮2齿廓曲线的弧长元素为:
齿轮1和齿轮2在点P处接触,在dt时间内滑动的弧长分别为ds1、ds2,则其相对滑动率定义为:
把方程(15a)和(15b)带入方程(16)得到了齿轮1和齿轮2共轭齿廓曲线在啮合点P处的相对滑动率。
步骤三、构造满足三次非均匀B样条啮合线的内啮合齿轮共轭齿廓曲线
当啮合线为三次非均匀B样条表示时,啮合线方程为分段三次多项式(4),设xf(u)=xi(u),yf(u)=yi(u),u∈[ui,ui+1],它满足啮合方程式(10)、(11)。
由于三次B样条曲线为分段三次多项式,不同定义域区间对应不同的三次多项式表达式,需对u进行分段积分。首先确定u∈[ui,ui+1]区间,然后分别对[uk,uk+1]、[uk+1,uk+2]、…、[ui-1,ui]、[ui,u]进行分段积分,得到的最终分段积分为:
得到的为齿轮1齿廓曲线的齿形点从初始位置到参与啮合所转过的弧度,即齿廓曲线齿形点与啮合点的对应关系,把式(17)带入式(6)得到齿轮1的齿廓曲线方程,由于 通过式(8)可以得到齿轮2的齿廓曲线方程。
同样,把式(17)带入式(13a)、(13b)、(14)和(15a)、(15b)、(16)可以得到齿轮1和齿轮2共轭齿廓曲线的相对法曲率和相对滑动率。
步骤四、内啮合齿轮共轭齿廓的的优化设计
在共轭齿轮齿廓的设计过程中,共轭齿廓曲线的光顺度及啮合特性取决于啮合线的几何性质。利用非均匀B样条曲线构造齿轮啮合线,其形状由控制多边形决定,可以通过对控制顶点位置的调整实现啮合线形状的修改。
根据曲线法曲率表达式,曲线的曲率变化主要受p″(u)的影响,光顺曲线的曲率分布需分段线性的,则p″(u)在节点区间上也应趋于线性变化。根据这个特性,建立代表B样条曲线光顺性的能量函数为:
为约束待修正控制顶点di 0(i=2,3,…,n-1)的扰动幅度,可对控制顶点的扰动幅度增加一个权重,建立的控制顶点扰动幅度目标函数为:
其中,w2i为控制顶点扰动幅度的权重。
相对滑动率是共轭齿廓曲线的重要啮合特性之一,建立关于相对滑动率的目标函数为:
对于齿轮1齿根齿廓与齿轮2齿顶齿廓啮合,有ds1-ds2<0,对于齿轮1齿顶齿廓与齿轮2齿根齿廓啮合,ds1-ds2>0,因此,在对于啮合线L1,|ds1-ds2|=ds1-ds2,对于啮合线L2,|ds1-ds2|=-(ds1-ds2)。把式(15a)、式(15b)带入式(19)并约去得到:
共轭齿廓的相对法曲率可以由式(14)表示,则建立的相对法曲率目标函数为:
根据共轭齿廓曲线的法曲率式(13a)、(13b),当齿轮1齿廓曲线法曲率表达式的分子、分母与齿轮2齿廓曲线法曲率表达式的分子、分母之差最小时,共轭齿廓曲线的相对法曲率会很小,约去公因子得到相对法曲率的目标函数为:
根据式(17)、(18)、(20)和(22),得到的整体目标函数为:
F=w1F1+F2+w3F′3+w4F′4 (24)
其中,w1、w3、w4分别为光顺度、相对滑动率及相对法曲率在目标函数中所占的权重。要使目标函数F最小,可令
这是2(n-2)元一次方程组,未知量为di,x、di,y(i=2,…,n-1)共2(n-2)个,其矩阵表达式为:
Md=C (26)
其中,M为系数矩阵,C为(n-1)×2向量,d的表达式为:
d=[dx dy]T dx=[d2x d3x … d(n-1)x] dy=[d2y d3y … d(n-1)y]
求解以上方程组得到经过优化后的B样条啮合线控制顶点的坐标值。
根据步骤一、步骤二和步骤三重新构造内啮合齿轮共轭齿廓。
步骤五、内啮合齿轮共轭齿廓曲线的干涉检查
若共轭齿轮齿廓曲线上出现奇异点,在该点处将产生轮齿根切。根据齿轮啮合原理,轮齿根切的判定条件为接触点在齿廓曲线上的滑动速度为零,即:
齿轮1轮齿根切的判定条件为:
xf(u)x′f(u)+(-R1+yf(u))y′f(u)=0 (29a)
齿轮2轮齿根切的判定条件为:
xf(u)x′f(u)+(-R2+yf(u))y′f(u)=0 (29b)
把式(4)代入式(29a)、(29b)可以判定满足B样条啮合线的内啮合齿轮共轭齿廓是否发生齿廓干涉。因此,以方法构造的内啮合共轭齿廓是否干涉是由啮合线的决定的,在构造共轭齿廓曲线之前就可以进行判定。鉴于内啮合齿轮的复杂性及B样条曲线的自由性,判断此方法构造的内啮合齿轮是否发生齿顶干涉和径向干涉,只有通过对构造共轭齿廓在每一个转动位置进行判定才能实现。此方法可以获得共轭齿廓上的齿形点坐标,通过在几个关键位置齿形点的转换位置计算及干涉判断可以完成构造共轭齿廓齿顶干涉及径向干涉的判定。
若发生干涉,根据干涉出现的位置可以对B样条啮合线上对应控制点进行调整,重新进行内啮合齿轮共轭齿廓曲线构造并进行干涉检查,直到不发生干涉为止。
步骤六、构造内外齿轮齿根齿廓到齿根圆的过渡曲线
确定内外齿轮轮齿齿根齿廓曲线在端点T0处的切向t,构造与齿根圆和齿根齿廓曲线在端点相切的圆弧过渡曲线T0Tf。,如图4所示。
下面具体举例说明:
以齿数z1=35,z2=50,m=4的内啮合齿轮为例,简要说明本方法的实施方式。其他参数为ha1=4mm,hf1=5mm,ha2=4mm,hf2=5mm。
步骤一、利用非均匀B样条构造内啮合齿轮啮合线。
设计的啮合线为外齿轮1右侧齿廓与内齿轮2左侧齿廓的啮合线,分为L1、L2两段,如表1所示,其控制顶点为di (L1)(i=0,1,2,…,n)、dij (L2)(j=0,1,2,…,m),满足d0 (L1)=d0 (L2)=(0,0),d2 (L1)、d1 (L1)、d0 (L1)、d1 (L2)、d2 (L2)在同一条直线l0上。dn (L1)、dm L2)分别在齿轮1和齿轮2的齿顶圆上。
把控制顶点的坐标带入式(4),则得到了啮合线方程的三次非均匀B样条曲线表示。
表1非均匀B样条啮合线控制顶点
步骤二、满足B样条啮合线的共轭齿轮齿廓曲线与啮合线的对应关系建立及其啮合性质分析
根据共轭原理,得到共轭方程表达式式(11),对共轭方程式(11)积分,得到啮合线上的点与共轭齿廓曲线上接触点的对应关系式(11)及式(12)。
根据齿轮1和齿轮2共轭齿廓曲线在啮合点P处的相对法曲率表达式(14)及相对滑动率表达式(16),得到的共轭齿廓曲线的相对法曲率曲线和相对滑动率曲线如图6、图7所示。
步骤三、满足设计啮合线的共轭齿轮齿廓曲线的构造
对于三次非均匀B样条曲线,其方程表达式为分段三次多项式,需对u进行分段积分。首先确定u∈[ui,ui+1]区间,然后分别对[uk,uk+1]、[uk+1,uk+2]、…、[ui-1,ui]、[ui,u]进行分段积分,得到的最终分段积分为:
得到的为齿轮1齿廓曲线的齿形点从初始位置到参与啮合所转过的弧度,即齿廓曲线齿形点与啮合点的对应关系,把式(30)带入式(6)得到齿轮1的齿廓曲线方程,由于 通过式(8)可以得到齿轮2的齿廓曲线方程。
以上得到的共轭齿廓曲线为中心轮1右侧齿廓与内齿圈2左侧齿廓,如图5(a)所示。
步骤四、内啮合齿轮共轭齿廓的优化设计
在内啮合齿轮共轭齿廓设计过程中,共轭齿廓曲线的光顺度及啮合特性取决于啮合线的几何性质。非均匀B样条啮合线的形状由控制顶点决定,可以通过对控制顶点位置的调整实现内啮合齿轮共轭齿廓的优化设计。
以B样条啮合线光顺性、控制顶点位置扰动、共轭齿廓相对法曲率及相对滑动率为优化目标建立目标函数式(23),通过求解目标函数得到了优化后的B样条啮合线控制顶点的位置坐标为di′(L1)(i=0,1,2,…,n)、dj′(L2)(j=0,1,2,…,m),如表2所示。
为了满足d0′(L1)=d0′(L2)=(0,0),d2′(L1)、d1′(L1)、d0′(L1)、d1′(L2)、d2′(L2)在同一条直线l0上。dn′(L1)、dm′(L2)分别在齿轮1和齿轮2的齿顶圆上,在设置控制顶点扰动时,使这几个点不产生扰动,即坐标值不变。
根据步骤一、步骤二和步骤三重新构造内啮合齿轮共轭齿廓及求解其相对法曲率和相对滑动率。得到了优化后的B样条啮合线及构造的共轭齿廓曲线如图5(b)所示,其相对法曲率和相对滑动率如图6、图7所示。
表2优化后非均匀B样条啮合线控制顶点
步骤五、设计共轭齿廓曲线的干涉检查
对于给定的三次非均匀B样条啮合线分段多项式方程表达式(4),带入中心轮1工作齿廓干涉的判定条件式(29a)和内齿圈2工作齿廓干涉的判定条件式(29b),发现实例中构造的B样条啮合线满足:
则根据该啮合线构造的共轭齿廓曲线不会产生工作齿廓干涉。通过对构造的共轭齿廓转动一定角度检查发现,此方法构造的共轭齿廓不发生齿顶干涉和径向干涉。
步骤六、齿轮工作齿廓到齿根圆的过渡曲线构造
中心齿根齿廓曲线的端点为T10(0.7102,-66.4380),在端点T10的切线方向t(0.1303,0.9915),中心齿轮的齿根圆中心O1(0,0),半径rf1=65mm,通过几何求解,与中心轮齿根圆外切于点Tf1的一段圆弧T10Tf1的中心O10(2.3987,-66.6599),半径r10=1.703mm,则中心 轮齿根齿廓到齿根圆的过渡圆弧方程表达式为:
内齿轮齿根齿廓曲线的端点为T20(-1.4718,-103.8182),在端点T20的切线方向t(-0.4384,-0.8988),外齿轮的齿根圆中心O2(0,0),半径rf2=105mm,通过几何求解,与内齿圈齿根圆外切于点Tf2的一段圆弧T20Tf2的中心O20(-3.2803,-102.9353),半径r20=2.0125mm,得到内齿圈齿根齿廓到齿根圆的过渡圆弧方程表达式为:
得到的中心轮和内齿圈轮齿如图8和图9所示,图8为中心轮和内齿圈单齿形状,分别由齿顶圆、齿廓曲线、过渡圆弧曲线及齿根圆组成,如图9所示中心轮和内齿圈整体齿形及啮合线。
本发明未详细阐述部分属于本领域的公知技术。
以上所述,仅为本发明部分具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本领域的人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。