CN113124800A - 阿基米德螺旋面蜗轮滚剃加工精度检测方法 - Google Patents

阿基米德螺旋面蜗轮滚剃加工精度检测方法 Download PDF

Info

Publication number
CN113124800A
CN113124800A CN202110426668.2A CN202110426668A CN113124800A CN 113124800 A CN113124800 A CN 113124800A CN 202110426668 A CN202110426668 A CN 202110426668A CN 113124800 A CN113124800 A CN 113124800A
Authority
CN
China
Prior art keywords
archimedes
worm
tooth surface
equation
worm gear
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
Application number
CN202110426668.2A
Other languages
English (en)
Other versions
CN113124800B (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.)
Chongqing University
Original Assignee
Chongqing University
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 Chongqing University filed Critical Chongqing University
Priority to CN202110426668.2A priority Critical patent/CN113124800B/zh
Publication of CN113124800A publication Critical patent/CN113124800A/zh
Application granted granted Critical
Publication of CN113124800B publication Critical patent/CN113124800B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16HGEARING
    • F16H55/00Elements with teeth or friction surfaces for conveying motion; Worms, pulleys or sheaves for gearing mechanisms
    • F16H55/02Toothed members; Worms
    • F16H55/22Toothed members; Worms for transmissions with crossing shafts, especially worms, worm-gears

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Gears, Cams (AREA)

Abstract

本发明公开了一种阿基米德螺旋面蜗轮滚剃加工精度检测方法,包括如下步骤:1)构建阿基米德蜗杆的齿面方程和齿面法矢方程;2)根据阿基米德蜗轮和阿基米德蜗杆之间的啮合关系构建阿基米德蜗轮的齿面方程和齿面法矢方程;3)根据阿基米德蜗轮和阿基米德蜗杆在任意接触点的相对运动速度矢量与法向矢量之间的正交关系,简化阿基米德蜗轮的齿面方程和齿面法矢方程;4)采用点阵式测量方法,在蜗轮齿面上选择若干网格点;利用简化得到的基米德蜗轮的齿面方程和齿面法矢方程分别计算网格点的理论坐标和法向矢量;利用测头测量得到网格点处的测量坐标;5)利用网格点的理论坐标、法向矢量和测量坐标,计算得到阿基米德蜗轮齿面的加工误差。

Description

阿基米德螺旋面蜗轮滚剃加工精度检测方法
技术领域
本发明属于机械误差分析技术领域,具体的为一种阿基米德螺旋面蜗轮加工精度检测方法。
背景技术
蜗轮蜗杆传动是一种以紧凑型、大速比为主要特点的传动类型,常用来传递两交错轴之间的运动和动力,在现代工业中被广泛应用。蜗轮一般采用滚齿机加工,主要有滚齿和飞刀切齿两种方法。制造精密蜗轮时,可在滚齿或切齿后再进行剃齿、珩齿或研齿等精整加工。
滚齿采用基本参数与工作蜗杆相同的蜗轮滚刀,按展成法原理(见齿轮加工)切出齿形。如果采用径向进给法滚齿,则滚刀与工件按Z2/Z1的传动比(Z1为工作蜗杆螺纹头数,Z2为蜗轮齿数)对滚,两者逐渐靠近直到其中心距等于工作蜗杆与蜗轮啮合时的中心距为止。采用切向进给法滚齿时,机床除保证刀具旋转外,还要有轴向进给;同时机床的工作台也要增加相应的附加转动,才能实现展成运动,这就要使用差动链。因此,切向进给法的加工精度一般不如径向进给法,但齿面质量较好,且不会产生根切现象。剃齿蜗轮剃齿刀的基本参数与工作蜗杆相同。蜗轮剃齿一般用滚齿机,可以由剃齿刀带动蜗轮自由剃齿,也可在机床传动链控制下强迫剃齿,剃后齿面质量和精度有所提高。珩齿和研齿蜗轮滚齿或飞刀切齿后,为了提高齿面质量、改善蜗轮与蜗杆啮合时的接触情况,可在滚齿机上珩齿或研齿。珩齿工具是用磨料与塑料、树脂的混合物浇铸在基体上而制成的珩磨蜗杆;研齿时用铸铁制成的研磨蜗杆加研磨剂与蜗轮对研。蜗轮的加工精度对蜗轮蜗杆传动至关重要,但现有技术中缺少用于检测蜗轮加工精度的方法。
发明内容
有鉴于此,本发明的目的在于提供一种阿基米德螺旋面蜗轮滚剃加工精度检测方法,能够检测蜗轮齿面加工精度。
为达到上述目的,本发明提供如下技术方案:
一种阿基米德螺旋面蜗轮滚剃加工精度检测方法,包括如下步骤:
1)根据阿基米德蜗杆的螺旋面端截形和轴截形方程构建阿基米德蜗杆的齿面方程和齿面法矢方程;
2)根据阿基米德蜗轮和阿基米德蜗杆之间的啮合关系,以阿基米德蜗杆的齿面方程和齿面法矢方程分别构建阿基米德蜗轮的齿面方程和齿面法矢方程;
3)根据阿基米德蜗轮和阿基米德蜗杆在任意接触点的相对运动速度矢量与法向矢量之间的正交关系,简化阿基米德蜗轮的齿面方程和齿面法矢方程;
4)采用点阵式测量方法,在蜗轮齿面上选择若干网格点;利用简化得到的基米德蜗轮的齿面方程和齿面法矢方程分别计算网格点的理论坐标和法向矢量;利用测头测量得到网格点处的测量坐标;
5)利用网格点的理论坐标、法向矢量和测量坐标,计算得到阿基米德蜗轮齿面的加工误差。
进一步,所述步骤1)中,阿基米德蜗杆的螺旋面端截形方程为:
Figure BDA0003029778310000021
阿基米德蜗杆的螺旋面轴截形方程为:
Figure BDA0003029778310000022
其中,α为齿形压力角;p为螺旋参数;u为自变量。
进一步,构建得到的阿基米德蜗杆的齿面方程为:
Figure BDA0003029778310000023
构建得到的阿基米德蜗杆的齿面法矢方程为:
Figure BDA0003029778310000024
Figure BDA0003029778310000025
Figure BDA0003029778310000026
其中,θ为参变量,表示端截线绕中心转过的角度值;rb为基圆的半径;
Figure BDA0003029778310000027
分别为阿基米德蜗杆的齿面法线
Figure BDA0003029778310000028
在x1、y1、z1坐标轴上的分量。
进一步,所述步骤2)中,构建得到的阿基米德蜗轮的齿面方程为:
Figure BDA0003029778310000031
构建得到的阿基米德蜗轮的齿面法矢方程为:
Figure BDA0003029778310000032
其中,
Figure BDA0003029778310000033
分别为空间啮合运动时阿基米德蜗杆和阿基米德蜗轮对应转过的角度,且
Figure BDA0003029778310000034
i12为阿基米德蜗杆和阿基米德蜗轮之间的传动比;
Figure BDA0003029778310000035
分别为阿基米德蜗轮的齿面法线
Figure BDA0003029778310000036
在x2、y2、z2坐标轴上的分量;A0为阿基米德蜗杆和阿基米德蜗轮啮合时的中心距;M21为蜗轮固联坐标系S2与蜗杆固联坐标系S1的之间的转换矩阵;且:
M21=M2p·Mp0·M01
其中,M01蜗杆固联坐标系S1与动坐标系S之间的转换矩阵;Mp0为固定坐标系S与蜗轮固联的动坐标系Sp之间的转换矩阵;M2p为与蜗轮固联的动坐标系S2与动坐标系Sp之间的转换矩阵。
进一步,所述步骤3)中,阿基米德蜗轮和阿基米德蜗杆在任意接触点的相对运动速度矢量v12为:
v12=v12xi+v12yj+v12zk
Figure BDA0003029778310000037
Figure BDA0003029778310000038
Figure BDA0003029778310000039
其中:ω1为阿基米德蜗杆绕其轴线的角速度;
由于阿基米德蜗轮和阿基米德蜗杆在任意接触点的相对运动速度矢量与法向矢量之间的正交关系,则有:
Figure BDA00030297783100000310
进一步,简化得到的阿基米德蜗轮的齿面方程为:
Figure BDA0003029778310000041
Figure BDA0003029778310000042
Figure BDA0003029778310000043
x1=u cosαcosθcos(u sinα/p)-u cosαsinθsin(u sinα/p)
y1=u cosαsinθsin(u sinα/p)+u cosαcosθsin(u sinα/p)
z1=pθ
Figure BDA0003029778310000044
Figure BDA0003029778310000045
简化得到的阿基米德蜗轮的齿面法矢方程为:
Figure BDA0003029778310000046
Figure BDA0003029778310000047
Figure BDA0003029778310000048
Figure BDA0003029778310000049
进一步,所述步骤4)中,设齿面上的任意网格点M所处位置的参数为βi、ri、(x2i,y2i,z2i),由此可得出在x2oz2平面内的βi、ri、z2i的关系式,如下:
z2i=ri tanβi
网格点M的坐标满足:
Figure BDA00030297783100000410
其中,Ri为网格点的半径,且满足:
Ri=A0-ri cosβi
由此得到用于计算自变量u、参变量θ和空间啮合运动时阿基米德蜗杆转过的角度
Figure BDA00030297783100000411
的联立方程组:
Figure BDA00030297783100000412
其中,ri为阿基米德蜗轮分度圆半径;βi为M点对应的在阿基米德蜗轮分度圆上转过的角度;z2i为M点在阿基米德蜗轮分度圆上的z坐标;x2i为M点在阿基米德蜗轮分度圆上的x坐标;y2i为M点在阿基米德蜗轮分度圆上的y坐标;Ri为阿基米德蜗轮分度圆半径;n为阿基米德蜗轮齿面法矢量;v12为阿基米德蜗轮相对于蜗杆的运动速度;
将计算得到的自变量u、参变量θ和空间啮合运动时阿基米德蜗杆转过的角度
Figure BDA0003029778310000051
的值代入到阿基米德蜗轮的齿面方程和齿面法矢方程中,求解得到任意网格点的理论坐标和法向矢量。
进一步,所述步骤4)中,由于测头为球体,测头测得的坐标值为测头球心的坐标值,则对测头半径进行补偿的方法为:
(x3,y3,z3)理论=(x'3,y'3,z'3)理论+r·(ηxyz)
其中,(x3,y3,z3)理论、(x'3,y'3,z'3)理论和(ηxyz)分别为蜗轮齿面的同一测量点在测量坐标系内的测头球心理论坐标值、测量点的理论坐标值和单位法矢量;r为测头的半径。
进一步,所述步骤5)中,计算阿基米德蜗轮齿面的加工误差时,需将蜗轮坐标系和测头坐标系进行变换,其变换方程为:
Figure BDA0003029778310000052
Figure BDA0003029778310000053
其中,(R0,T0,Z0)为测头坐标系内的坐标值;(nr0,nt0,nz0)为测头坐标系内的法向矢量;(x2,y2,z2)为蜗轮坐标系内的坐标值;(nx2,ny2,nz2)为蜗轮坐标系内的法向矢量;β为蜗轮坐标系转动的角度,且蜗轮坐标系转动β角度后与测量机坐标系完全重合:
β=arctan(y5,3/x5,3)
其中,y5,3为蜗轮坐标系中的齿面纵坐标;x5,3为蜗轮坐标系中的齿面横坐标。
进一步,针对阿基米德蜗轮齿面上的同一网格点,设利用阿基米德蜗轮的齿面方程和齿面法矢方程计算得到的理论坐标值和理论法矢变换到测头坐标系内后分别为(R,T,Z)和(nr,nt,nz),利用测头测量得到的该网格点的坐标值(R',T',Z'),则齿面测量误差为:
δ=(R'-R)nr+(T'-T)nt+(Z'-Z)nz-(θ'-θ)R'·nt
其中,(R'-R)nr、(T'-T)nt和(Z'-Z)nz表示R、T和Z三个方向上的误差分量,(θ'-θ)R'·nt表示转动角度在T方向的误差分量;θ为阿基米德蜗轮理想情况下转过的角度;θ'为阿基米德蜗轮实际情况下转过的角度。
本发明的有益效果在于:
本发明的阿基米德螺旋面蜗轮加工精度检测方法,通过构建阿基米德蜗杆的齿面方程和齿面法矢方程,利用蜗轮蜗杆之间的啮合关系,从而可以通过蜗杆的齿面方程和齿面法矢方程构建阿基米德蜗轮的齿面方程和齿面法矢方程;从而利用点阵测量方法,根据蜗轮齿面方程和蜗轮齿面法矢方程分别计算网格点的理论坐标和法向矢量,而后与测头测量得到的网格点处的测量坐标进行比较,从而得到蜗轮的加工误差,即本发明的阿基米德螺旋面蜗轮加工精度检测方法能够测量蜗轮的加工齿轮与理论齿面之间的误差,从而得到蜗轮的加工精度。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚,本发明提供如下附图进行说明:
图1为阿基米德蜗杆螺旋面的端截形的示意图;
图2为蜗轮蜗杆空间中的啮合运动关系;
图3为蜗轮蜗杆空间中的相对速度;
图4为测头球心等距面的示意图;
图5为阿基米德螺旋面(ZA型)蜗轮的左齿面图;
图6为ZN齿面截交线的渐开线拟合曲线图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明,以使本领域的技术人员可以更好的理解本发明并能予以实施,但所举实施例不作为对本发明的限定。
本实施例的阿基米德螺旋面蜗轮滚剃加工精度检测方法,包括如下步骤:
1)根据阿基米德蜗杆的螺旋面端截形和轴截形方程构建阿基米德蜗杆的齿面方程和齿面法矢方程。
依据蜗杆类型、蜗杆头数z1、蜗轮齿数、模数m、齿形压力角α、蜗杆分度圆直径d1、螺旋升角γ、齿宽b、螺距、传动比i21、中心距A0、齿顶圆直径da1、齿根圆直径df1,确定阿基米德蜗杆的螺旋面端截形和轴截形,具体的:
阿基米德蜗杆的螺旋面端截形方程为:
Figure BDA0003029778310000061
阿基米德蜗杆的螺旋面轴截形方程为:
Figure BDA0003029778310000062
其中,α为齿形压力角;p为螺旋参数,且p=导程pz/2π;u为自变量。
根据阿基米德蜗杆的螺旋面端截形和轴截形,可知阿基米德蜗杆的齿面端面截形为一条阿基米德螺旋曲线。如图1所示,该阿基米德螺旋曲线以半径rb作为基圆,以圆上点b作为起点的阿基米德螺旋曲线bc与x轴的夹角为0°,阿基米德螺旋曲线上的任一点P处的法线Pa与基圆相切于点a,取角u=∠aob作为自变量,由阿基米德螺旋曲线的特性有Pa=rbu,故可得阿基米德螺旋曲线bc的方程式为:
Figure BDA0003029778310000071
将的阿基米德螺旋曲线bc作以z轴为中心的螺旋运动,转过角度记为θ,则会形成一个螺旋曲面,相当于蜗杆的齿面,故构建得到的阿基米德蜗杆的齿面(右旋)方程r1为:
Figure BDA0003029778310000072
根据蜗轮蜗杆副空间啮合原理可知,求解得到阿基米德蜗轮齿面的方程式,首先需要对阿基米德蜗杆齿面法线
Figure BDA0003029778310000073
的方程式进行求解。由参数方程的求导运算规律对式(3)进行求导计算,可得其齿面法线方程。且
Figure BDA0003029778310000074
分别为阿基米德蜗杆齿面法线
Figure BDA0003029778310000075
在x1、y1、z1坐标轴上的分量,即构建得到的阿基米德蜗杆的齿面法矢方程为:
Figure BDA0003029778310000076
Figure BDA0003029778310000077
Figure BDA0003029778310000078
其中,θ为参变量,表示端截线绕中心转过的角度值;rb为基圆的半径;
Figure BDA0003029778310000079
分别为阿基米德蜗杆的齿面法线
Figure BDA00030297783100000710
在x1、y1、z1坐标轴上的分量。
2)根据阿基米德蜗轮和阿基米德蜗杆之间的啮合关系,以阿基米德蜗杆的齿面方程和齿面法矢方程分别构建阿基米德蜗轮的齿面方程和齿面法矢方程。
蜗轮蜗杆的空间中的啮合运动关系,建立其空间坐标系。如图2所示,坐标系S1(O1-x1,y1,z1)与蜗杆固定连接,坐标轴z1与蜗杆中心线相重合,且蜗杆可绕坐标轴做旋转运动;S2(O2-x2,y2,z2)坐标系与蜗轮固定连接,坐标轴z2与蜗轮中心线重合,且蜗轮可绕坐标轴z2做旋转运动;S(O-x,y,z)与Sp(Op-xp,yp,zp)为两个不动坐标系,其中坐标轴z与蜗杆的旋转中心线z1相重合,坐标轴zp与蜗轮的旋转中心线z2相重合,且空间中z与zp的异面夹角Σ=90°;坐标轴x与xp相重合,OOp为其啮合时的中心距,记OOp=A0
其中,角
Figure BDA0003029778310000081
分别为空间啮合运动时蜗杆和蜗轮对应转过的角度,即在蜗杆和蜗轮在持续啮合时,当蜗杆以坐标轴z1为中心轴线旋转
Figure BDA0003029778310000082
角度时对应的蜗轮以坐标轴z2为中心轴线旋转
Figure BDA0003029778310000083
角度,且
Figure BDA0003029778310000084
之间满足传动比关系
Figure BDA0003029778310000085
i12为阿基米德蜗杆和阿基米德蜗轮之间的传动比。
将蜗杆齿面转换到蜗轮齿面S2(O2-x2,y2,z2),则有:
Figure BDA0003029778310000086
式中,M21为蜗轮固联坐标系S2与蜗杆固联坐标系S1的之间的转换矩阵。
根据空间坐标位置关系,推算出蜗杆坐标系S1到蜗轮坐标系S2的变化关系式为:
Figure BDA0003029778310000087
式中,M01蜗杆固定坐标系S1与动坐标系S的转换矩阵:
Figure BDA0003029778310000088
Mp0为固定坐标系S与蜗轮固联的动坐标系Sp的转换矩阵:
Figure BDA0003029778310000091
M2p为与蜗轮固联的动坐标系Sp与坐标系S2的转换矩阵:
Figure BDA0003029778310000092
所以,通过以上空间坐标矩阵的运算推导,可求得蜗轮固联坐标系S2与蜗杆固联坐标系S1的之间的转换矩阵M21为:
Figure BDA0003029778310000093
综上所示,将蜗杆齿面转换到蜗轮齿面S2(O2-x2,y2,z2),构建得到的阿基米德蜗轮的齿面方程为:
Figure BDA0003029778310000094
构建得到的阿基米德蜗轮的齿面法矢方程为:
Figure BDA0003029778310000097
其中,
Figure BDA0003029778310000095
分别为阿基米德蜗轮的齿面法线
Figure BDA0003029778310000096
在x2、y2、z2坐标轴上的分量;A0为阿基米德蜗杆和阿基米德蜗轮啮合时的中心距;M21为蜗轮固联坐标系S2与蜗杆固联坐标系S1的之间的转换矩阵;
3)根据阿基米德蜗轮和阿基米德蜗杆在任意接触点的相对运动速度矢量与法向矢量之间的正交关系,简化阿基米德蜗轮的齿面方程和齿面法矢方程。
(1)阿基米德蜗轮和阿基米德蜗杆在任意接触点的相对运动速度关系
蜗轮蜗杆传动为空间单自由度的运动形式。蜗杆与蜗轮工作时在沿任意方向上无位移,均只做绕各自中心轴的回转运动。如图3所示,设蜗轮蜗杆副工作时,蜗杆以z1为中心轴定角速度ω1做回转运动,蜗轮以z2为中心轴定角速度ω2做回转运动。
由图3所示的运动关系可知,蜗杆的回转角速度与其沿轴线方向移动速度的矢量式为:
Figure BDA0003029778310000101
式中,ω1表示工件角速度的模;
Figure BDA0003029778310000102
表示蜗杆沿轴线移动速度的模,
Figure BDA0003029778310000103
蜗轮的回转角速度与其沿轴线方向移动速度的矢量式为:
Figure BDA0003029778310000104
式中,ω2表示工件角速度的模;
Figure BDA0003029778310000105
表示蜗杆沿轴线移动速度的模,
Figure BDA0003029778310000106
空间中有一点M,在坐标系S中的坐标值为(x,y,z),且OOp=A0i,则有:
Figure BDA0003029778310000107
M点随蜗杆运动时的速度v1
Figure BDA0003029778310000108
M点随蜗轮运动时的速度v2为:
Figure BDA0003029778310000109
则蜗杆与蜗轮传动中,在其啮合点M处的相对运动速度v12的矢量式为:
Figure BDA00030297783100001010
所以:
Figure BDA00030297783100001011
所以,有如下关系
Figure BDA00030297783100001012
用坐标系S1(O1-x1,y1,z1)的参数表示,根据的坐标系S与S1之间的关系,即:
Figure BDA00030297783100001013
则:
v12=v12xi+v12yj+v12zk (22-a)
利用Σ=90°,
Figure BDA0003029778310000111
等对上式进行简化,得:
Figure BDA0003029778310000112
Figure BDA0003029778310000113
Figure BDA0003029778310000114
(2)阿基米德蜗轮和阿基米德蜗杆之间的啮合关系
根据空间齿轮啮合理论可知,蜗轮蜗杆传动中的任一时刻两个齿面间都为线接触,在其接触线上的切点除存在两齿面的公法平面及公法线nI。两个齿面在传动中总是保持连续的啮合状态,且不发生脱离与干涉情况,则证明在该点上两齿面间的相对运动速度矢量v12与其法线矢量的正交关系,根据此结论推导方程式便可的其啮合方程式。
由其正交关系可知,在接触点处应满足如下方程式:
nI·v12=0 (23)
所以,阿基米德蜗杆蜗轮啮合方程式为:
Figure BDA0003029778310000115
(3)简化得到的阿基米德蜗轮的齿面方程和齿面法矢方程
根据齿轮啮合原理可知,一对蜗轮蜗杆副在整个啮合运动过程中,两齿面始终保持相切关系并各自按照所设定的运动规律工作。故,在运动过程中蜗轮蜗杆副齿面为一对共轭齿面。利用蜗轮固联坐标系S2与蜗杆固联坐标系S1的之间的转换矩阵M21(式(10))将蜗杆齿面方程(式(3))变换到蜗轮坐标系中,再联立啮合方程式(24)并进行化简计算,简化得到的阿基米德蜗轮的齿面方程为::
Figure BDA0003029778310000121
Figure BDA0003029778310000122
Figure BDA0003029778310000123
x1=u cosαcosθcos(u sinα/p)-u cosαsinθsin(u sinα/p)
y1=u cosαsinθsin(u sinα/p)+u cosαcosθsin(u sinα/p)
z1=pθ
Figure BDA0003029778310000124
Figure BDA0003029778310000125
根据两个齿面间的共轭关系可知,在两齿面的公切点上存在公法线,公法线同时垂直于蜗轮齿面与蜗杆齿面,所以将蜗杆的齿面法矢方程式进行坐标变换就能推导出蜗轮齿面的法矢方程式。即将蜗杆齿面法矢方程式(4)通过蜗轮固联坐标系S2与蜗杆固联坐标系S1的之间的转换矩阵M21(式(10))处理后,求得简化得到的阿基米德蜗轮的齿面法矢方程为:
Figure BDA0003029778310000126
Figure BDA0003029778310000127
Figure BDA0003029778310000128
Figure BDA0003029778310000129
4)采用点阵式测量方法,在蜗轮齿面上选择若干网格点;利用简化得到的基米德蜗轮的齿面方程和齿面法矢方程分别计算网格点的理论坐标和法向矢量;利用测头测量得到网格点处的测量坐标。
本实施例采用点阵式测量方法,在蜗轮齿面上选择若干网格点,将蜗轮齿面划分为若干个测量单元;按照一定的间距在齿面的两个方向将齿面划分成网格,对网格点进行理论坐标的计算,再利用测头点对点进行测量,从而将整个齿面的误差情况得到准确的检测和评价。三维的空间中的蜗轮齿面是不规则的复杂曲面,很难在其上直接网格划分并计算三维的网格点坐标及其法矢,所以采取先二维网格划分后转成三维曲面的方法进行网格划分。二维网格划分其实是将三维空间中的蜗轮齿面进行平面投影,用两个变量在两个方向上将齿面划分为网格状。二维网格划分完成后又需把齿面还原为三维曲面形状,还原的方法是将已知的两个坐标项带入齿面方程式求解未知坐标项。在划分网格时还需要注意对齿根与齿顶处缩进值及齿面上下两端处的缩进值得设定,因为有些蜗轮齿面的不规则可能使得被网格划分后的理论齿面网格点落在齿面外侧。
用两个参数对投影面进行二维网格的划分,其中一个为直线与X1轴线的夹角的角度参数β,另一个为竖直方向上从蜗杆齿根到齿顶圆的同心圆的半径参数r。β可根据蜗轮齿宽缩进后进行取值,等分为9等份(本文在齿宽方向上缩进以蜗轮齿宽b2为准,缩进比例取0.15);r从蜗杆齿根到齿顶进行取值,分为5份,中间取蜗轮得分度圆半径(半径方向上缩进比例取0.05,以分度圆半径为基准,分别由齿顶和齿根两端向分度圆处中心进行缩进)。M为齿面上任意一点,设其所处位置的参数为βi、ri、(x2i,y2i,z2i),由此可得出在x2oz2平面内的βi、ri、z2i的关系式,如下:
z2i=ri tanβi (27)
网格点M的坐标满足:
Figure BDA0003029778310000131
其中,Ri为网格点的半径,且满足:
Ri=A0-ri cosβi (29)
网格点计算包括计算网格点的三个坐标值及点的法矢的求解。根据蜗轮齿面方程式可知,方程存在三个未知参数u、θ、
Figure BDA0003029778310000132
故需要三个方程进行联立求解。式(27)和(28)为其中两个方程式,将式(25)中的x2、y2、z2分别带入方程式(27)和(28)中,并联立啮合方程式(24),由此得到用于计算自变量u、参变量θ和空间啮合运动时阿基米德蜗杆转过的角度
Figure BDA0003029778310000133
的联立方程组:
Figure BDA0003029778310000134
其中,ri为阿基米德蜗轮分度圆半径;βi为M点对应的在阿基米德蜗轮分度圆上转过的角度;z2i为M点在阿基米德蜗轮分度圆上的z坐标;x2i为M点在阿基米德蜗轮分度圆上的x坐标;y2i为M点在阿基米德蜗轮分度圆上的y坐标;Ri为阿基米德蜗轮分度圆半径;n为阿基米德蜗轮齿面法矢量;v12为阿基米德蜗轮相对于蜗杆的运动速度;
将计算得到的自变量u、参变量θ和空间啮合运动时阿基米德蜗杆转过的角度
Figure BDA0003029778310000135
的值代入到阿基米德蜗轮的齿面方程和齿面法矢方程中,求解得到任意网格点的理论坐标和法向矢量。
本实施例对蜗轮齿面进行实际测量时选用测头为接触式一维测头。如图4所示,由于测头为一个球体,在空间中占有一定的体积,测头球体中心走过的轨迹坐标,而不是齿面被测点处或与测头相接触的点,所以获得的实测数据不能直接与理论数据进行比较,并且在控制测量运动时首先得理论计算处测头中心实际的位置坐标值。计算蜗轮齿面误差前必须先消除测头半径在对数据采集时的影响。为保证测头能准确测量到目标测量点。本论文选用等距面法对测量中的测头半径做补偿计算。球形测头运动到测量曲面上的指定位置时,测头和齿面的接触点即为理论测量点,而测头在齿面上走过轨迹的所有球心位置会构成了一个与被测齿面的等距曲面,且知距离为测头的半径值。在曲面被测点理论坐标及法矢已知情况下,根据图中的位置关系就能计算出理论的测头位置,继而完成对测头半径的补偿,则对测头半径进行补偿的方法如下:
(x3,y3,z3)理论=(x'3,y'3,z'3)理论+r·(ηxyz) (31)
其中,(x3,y3,z3)理论、(x'3,y'3,z'3)理论和(ηxyz)分别为蜗轮齿面的同一测量点在测量坐标系内的测头球心理论坐标值、测量点的理论坐标值和单位法矢量;r为测头的半径;
测头进行实际测量时,测头的读取数据(x1,y1,z1)实际与理论数据(x1,y1,z1)理论的差值在测量点处法矢量上的投影值就表示了该平面的法向误差值:
Figure BDA0003029778310000141
其中,nδ为读取数据(x1,y1,z1)实际与理论数据(x1,y1,z1)理论的差值;n为法矢量;δ法向为法向误差值。
5)利用网格点的理论坐标、法向矢量和测量坐标,计算得到阿基米德蜗轮齿面的加工误差。
计算阿基米德蜗轮齿面的加工误差时,需将蜗轮坐标系和测头坐标系进行变换,其变换方程为:
Figure BDA0003029778310000142
Figure BDA0003029778310000143
其中,(R0,T0,Z0)为测头坐标系内的坐标值;(nr0,nt0,nz0)为测头坐标系内的法向矢量;(x2,y2,z2)为蜗轮坐标系内的坐标值;(nx2,ny2,nz2)为蜗轮坐标系内的法向矢量;β为蜗轮坐标系转动的角度,且蜗轮坐标系转动β角度后与测量机坐标系完全重合:
β=arctan(y5,3/x5,3) (35)
其中,y5,3为蜗轮坐标系中的齿面纵坐标;x5,3为蜗轮坐标系中的齿面横坐标。
具体的,在CNC齿轮测量中心上测量蜗轮齿面误差时,需要控制测头运动对每一个齿面网格点进行测量。齿轮测量中心的坐标系包含四个坐标轴,R、T、Z和θ,在测量过程中,测头处在T向受力方向,位置在T=0处;在R和Z方向运动来找到指定理论位置,再由回转工作台转动变化θ来进行测量。所以在测量时,R轴和Z轴构成了极坐标测量的直线轴,θ轴为极坐标测量的旋转轴,故齿轮测量中心的极坐标系为S(O-x,y,z)。应用极坐标法测量齿面时,首先需要把测量坐标系下求解得到的理论齿面被测网格点经过坐标转换到测量机极坐标系下的ZOR平面上,成为一个平面上的网格点。对此的具体操作方法是根据理论网格点的坐标与单位法矢,分别计算每一点与ZOR平面内的夹角,将每一点旋转各自的角度转换到ZOR平面上,对其单位法矢也做同样的操作。设齿面上一被测网格点为P0,根据以上转化处理可知该点在测量坐标系下的坐标值与法矢分别为(R0ij,T0ij,Z0ij)和(nr0ij,nt0ij,nz0ij),将其坐标值与法矢一同旋转θ角度,求得测量极坐标系下P0点的坐标值(Rθij,Tθij,Zθij)与法矢(nθij,ntθij,nzθij),其转换关系如下:
Figure BDA0003029778310000151
Figure BDA0003029778310000152
在对网格点的坐标及法矢做以上数据转化时,其实是以Z轴为中心的旋转处理,则可知每一个网格点处有Zθij=0,则可知:
Figure BDA0003029778310000153
具体的,本实施例针对蜗轮齿面上的同一网格点,利用齿面方程计算得到的理论坐标值和理论法矢变换到测头坐标系内后分别为(R,T,Z)和(nr,nt,nz),利用测头测量得到的该网格点的坐标值(R',T',Z'),则齿面测量误差为:
δ=(R'-R)nr+(T'-T)nt+(Z'-Z)nz-(θ'-θ)R'·nt (38)
其中,(R'-R)nr、(T'-T)nt和(Z'-Z)nz表示R、T和Z三个方向上的误差分量,(θ'-θ)R'·nt表示转动角度在T方向的误差分量;θ为阿基米德蜗轮理想情况下转过的角度;θ'为阿基米德蜗轮实际情况下转过的角度。
本实施例的阿基米德螺旋面蜗轮加工精度检测方法,通过构建阿基米德蜗杆的齿面方程和齿面法矢方程,利用蜗轮蜗杆之间的啮合关系,根据蜗杆的齿面方程和齿面法矢方程构建阿基米德蜗轮的齿面方程和齿面法矢方程;从而利用点阵测量方法,能够利用蜗轮齿面方程和蜗轮齿面法矢方程分别计算网格点的理论坐标和法向矢量,而后与测头测量得到的网格点处的测量坐标进行比较,从而得到蜗轮的加工误差,即实施例的阿基米德螺旋面蜗轮加工精度检测方法能够测量蜗轮的加工齿轮与理论齿面之间的误差,从而得到蜗轮的加工精度。
以上所述实施例仅是为充分说明本发明而所举的较佳的实施例,本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换,均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。

Claims (10)

1.一种阿基米德螺旋面蜗轮滚剃加工精度检测方法,其特征在于:包括如下步骤:
1)根据阿基米德蜗杆的螺旋面端截形和轴截形方程构建阿基米德蜗杆的齿面方程和齿面法矢方程;
2)根据阿基米德蜗轮和阿基米德蜗杆之间的啮合关系,以阿基米德蜗杆的齿面方程和齿面法矢方程分别构建阿基米德蜗轮的齿面方程和齿面法矢方程;
3)根据阿基米德蜗轮和阿基米德蜗杆在任意接触点的相对运动速度矢量与法向矢量之间的正交关系,简化阿基米德蜗轮的齿面方程和齿面法矢方程;
4)采用点阵式测量方法,在蜗轮齿面上选择若干网格点;利用简化得到的基米德蜗轮的齿面方程和齿面法矢方程分别计算网格点的理论坐标和法向矢量;利用测头测量得到网格点处的测量坐标;
5)利用网格点的理论坐标、法向矢量和测量坐标,计算得到阿基米德蜗轮齿面的加工误差。
2.根据权利要求1所述的阿基米德螺旋面蜗轮滚剃加工精度检测方法,其特征在于:所述步骤1)中,阿基米德蜗杆的螺旋面端截形方程为:
Figure FDA0003029778300000011
阿基米德蜗杆的螺旋面轴截形方程为:
Figure FDA0003029778300000012
其中,α为齿形压力角;p为螺旋参数;u为自变量。
3.根据权利要求2所述的阿基米德螺旋面蜗轮滚剃加工精度检测方法,其特征在于:构建得到的阿基米德蜗杆的齿面方程为:
Figure FDA0003029778300000013
构建得到的阿基米德蜗杆的齿面法矢方程为:
Figure FDA0003029778300000014
Figure FDA0003029778300000021
Figure FDA0003029778300000022
其中,θ为参变量,表示端截线绕中心转过的角度值;rb为基圆的半径;
Figure FDA0003029778300000023
分别为阿基米德蜗杆的齿面法线
Figure FDA0003029778300000024
在x1、y1、z1坐标轴上的分量。
4.根据权利要求3所述的阿基米德螺旋面蜗轮滚剃加工精度检测方法,其特征在于:所述步骤2)中,构建得到的阿基米德蜗轮的齿面方程为:
Figure FDA0003029778300000025
构建得到的阿基米德蜗轮的齿面法矢方程为:
Figure FDA0003029778300000026
其中,
Figure FDA0003029778300000027
分别为空间啮合运动时阿基米德蜗杆和阿基米德蜗轮对应转过的角度,且
Figure FDA0003029778300000028
i12为阿基米德蜗杆和阿基米德蜗轮之间的传动比;
Figure FDA0003029778300000029
分别为阿基米德蜗轮的齿面法线
Figure FDA00030297783000000210
在x2、y2、z2坐标轴上的分量;A0为阿基米德蜗杆和阿基米德蜗轮啮合时的中心距;M21为蜗轮固联坐标系S2与蜗杆固联坐标系S1的之间的转换矩阵;且:
M21=M2p·Mp0·M01
其中,M01蜗杆固联坐标系S1与动坐标系S之间的转换矩阵;Mp0为固定坐标系S与蜗轮固联的动坐标系Sp之间的转换矩阵;M2p为与蜗轮固联的动坐标系S2与动坐标系Sp之间的转换矩阵。
5.根据权利要求4所述的阿基米德螺旋面蜗轮滚剃加工精度检测方法,其特征在于:所述步骤3)中,阿基米德蜗轮和阿基米德蜗杆在任意接触点的相对运动速度矢量v12为:
v12=v12xi+v12yj+v12zk
Figure FDA00030297783000000211
Figure FDA0003029778300000031
Figure FDA0003029778300000032
其中:ω1为阿基米德蜗杆绕其轴线的角速度;
由于阿基米德蜗轮和阿基米德蜗杆在任意接触点的相对运动速度矢量与法向矢量之间的正交关系,则有:
Figure FDA0003029778300000033
6.根据权利要求5所述的阿基米德螺旋面蜗轮滚剃加工精度检测方法,其特征在于:简化得到的阿基米德蜗轮的齿面方程为:
Figure FDA0003029778300000034
Figure FDA0003029778300000035
Figure FDA0003029778300000036
x1=ucosαcosθcos(usinα/p)-ucosαsinθsin(usinα/p)
y1=ucosαsinθsin(usinα/p)+ucosαcosθsin(usinα/p)
z1=pθ
Figure FDA0003029778300000037
Figure FDA0003029778300000038
简化得到的阿基米德蜗轮的齿面法矢方程为:
Figure FDA0003029778300000039
Figure FDA00030297783000000310
Figure FDA00030297783000000311
Figure FDA00030297783000000312
7.根据权利要求6所述的阿基米德螺旋面蜗轮滚剃加工精度检测方法,其特征在于:所述步骤4)中,设齿面上的任意网格点M所处位置的参数为βi、ri、(x2i,y2i,z2i),由此可得出在x2oz2平面内的βi、ri、z2i的关系式,如下:
z2i=ri tanβi
网格点M的坐标满足:
Figure FDA00030297783000000313
其中,Ri为网格点的半径,且满足:
Ri=A0-ri cosβi
由此得到用于计算自变量u、参变量θ和空间啮合运动时阿基米德蜗杆转过的角度
Figure FDA0003029778300000041
的联立方程组:
Figure FDA0003029778300000042
其中,ri为阿基米德蜗轮分度圆半径;βi为M点对应的在阿基米德蜗轮分度圆上转过的角度;z2i为M点在阿基米德蜗轮分度圆上的z坐标;x2i为M点在阿基米德蜗轮分度圆上的x坐标;y2i为M点在阿基米德蜗轮分度圆上的y坐标;Ri为阿基米德蜗轮分度圆半径;n为阿基米德蜗轮齿面法矢量;v12为阿基米德蜗轮相对于蜗杆的运动速度;
将计算得到的自变量u、参变量θ和空间啮合运动时阿基米德蜗杆转过的角度
Figure FDA0003029778300000045
的值代入到阿基米德蜗轮的齿面方程和齿面法矢方程中,求解得到任意网格点的理论坐标和法向矢量。
8.根据权利要求1所述的阿基米德螺旋面蜗轮滚剃加工精度检测方法,其特征在于:所述步骤4)中,由于测头为球体,测头测得的坐标值为测头球心的坐标值,则对测头半径进行补偿的方法为:
(x3,y3,z3)理论=(x'3,y'3,z'3)理论+r·(ηxyz)
其中,(x3,y3,z3)理论、(x'3,y'3,z'3)理论和(ηxyz)分别为蜗轮齿面的同一测量点在测量坐标系内的测头球心理论坐标值、测量点的理论坐标值和单位法矢量;r为测头的半径。
9.根据权利要求1所述的阿基米德螺旋面蜗轮滚剃加工精度检测方法,其特征在于:所述步骤5)中,计算阿基米德蜗轮齿面的加工误差时,需将蜗轮坐标系和测头坐标系进行变换,其变换方程为:
Figure FDA0003029778300000043
Figure FDA0003029778300000044
其中,(R0,T0,Z0)为测头坐标系内的坐标值;(nr0,nt0,nz0)为测头坐标系内的法向矢量;(x2,y2,z2)为蜗轮坐标系内的坐标值;(nx2,ny2,nz2)为蜗轮坐标系内的法向矢量;β为蜗轮坐标系转动的角度,且蜗轮坐标系转动β角度后与测量机坐标系完全重合:
β=arctan(y5,3/x5,3)
其中,y5,3为蜗轮坐标系中的齿面纵坐标;x5,3为蜗轮坐标系中的齿面横坐标。
10.根据权利要求9所述的阿基米德螺旋面蜗轮滚剃加工精度检测方法,其特征在于:针对阿基米德蜗轮齿面上的同一网格点,设利用阿基米德蜗轮的齿面方程和齿面法矢方程计算得到的理论坐标值和理论法矢变换到测头坐标系内后分别为(R,T,Z)和(nr,nt,nz),利用测头测量得到的该网格点的坐标值(R',T',Z'),则齿面测量误差为:
δ=(R'-R)nr+(T'-T)nt+(Z'-Z)nz-(θ'-θ)R'·nt
其中,(R'-R)nr、(T'-T)nt和(Z'-Z)nz表示R、T和Z三个方向上的误差分量,(θ'-θ)R'·nt表示转动角度在T方向的误差分量;θ为阿基米德蜗轮理想情况下转过的角度;θ'为阿基米德蜗轮实际情况下转过的角度。
CN202110426668.2A 2021-04-20 2021-04-20 阿基米德螺旋面蜗轮滚剃加工精度检测方法 Active CN113124800B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110426668.2A CN113124800B (zh) 2021-04-20 2021-04-20 阿基米德螺旋面蜗轮滚剃加工精度检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110426668.2A CN113124800B (zh) 2021-04-20 2021-04-20 阿基米德螺旋面蜗轮滚剃加工精度检测方法

Publications (2)

Publication Number Publication Date
CN113124800A true CN113124800A (zh) 2021-07-16
CN113124800B CN113124800B (zh) 2022-03-01

Family

ID=76778465

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110426668.2A Active CN113124800B (zh) 2021-04-20 2021-04-20 阿基米德螺旋面蜗轮滚剃加工精度检测方法

Country Status (1)

Country Link
CN (1) CN113124800B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113639682A (zh) * 2021-09-24 2021-11-12 重庆大学 一种面齿轮在机检测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS57106801A (en) * 1980-12-24 1982-07-02 Toyoyama Akira Measuring method of relative error between worm and hob
CN103791869A (zh) * 2014-02-18 2014-05-14 北京工业大学 一种平面二次包络环面滚刀的测量方法
CN109604738A (zh) * 2019-02-01 2019-04-12 淮阴工学院 一种基于尼曼蜗轮数学模型的高效侧刃精铣加工方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS57106801A (en) * 1980-12-24 1982-07-02 Toyoyama Akira Measuring method of relative error between worm and hob
CN103791869A (zh) * 2014-02-18 2014-05-14 北京工业大学 一种平面二次包络环面滚刀的测量方法
CN109604738A (zh) * 2019-02-01 2019-04-12 淮阴工学院 一种基于尼曼蜗轮数学模型的高效侧刃精铣加工方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
宋丹等: "阿基米德涡轮齿面建模与计算", 《西安理工大学学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113639682A (zh) * 2021-09-24 2021-11-12 重庆大学 一种面齿轮在机检测方法

Also Published As

Publication number Publication date
CN113124800B (zh) 2022-03-01

Similar Documents

Publication Publication Date Title
KR101721969B1 (ko) 기어들에서 인벌류트들의 위치 결정을 위한 방법
US11247305B2 (en) Numerical control device and numerical control meihod for machine tool
WO2022135612A1 (zh) 一种渐开螺旋圆柱齿轮的法向啮合齿形测量方法
Wang et al. An efficient honing method for face gear with tooth profile modification
CN112580160B (zh) 一种用于成形磨齿机在机测量系统标定方法
CN113124800B (zh) 阿基米德螺旋面蜗轮滚剃加工精度检测方法
JP2019512684A (ja) ウォームギアの測定
Shen et al. Research on shaving processing of spiroid face gear
CN114091210A (zh) 考虑齿廓修形的展成磨齿齿面扭曲预测模型构建方法
Boral et al. Machining of spur gears using a special milling cutter
KR101200171B1 (ko) 5축 cνc 가공기계를 이용한 드럼 캠의 황삭가공을 위한 후처리방법
Lin et al. A five-axis CNC machining method of orthogonal variable transmission ratio face gear
CN109446710B (zh) 一种弧齿锥齿轮载荷传动误差数值计算方法
CN114754698B (zh) 面齿轮齿面测量点规划及在机测量方法
CN113639682B (zh) 一种面齿轮在机检测方法
CN110645933A (zh) 一种直齿圆锥齿轮齿面参数测量方法
KR101199086B1 (ko) 5축 cνc 가공기계를 이용한 롤러 기어 캠의 가공을 위한 후처리방법
Chlost et al. A new method of the positioning and analysis of the roughness deviation in five-axis milling of external cylindrical gear
CN114769737A (zh) 一种少齿数齿轮成型磨削加工方法
Shih et al. Precision evaluation for cycloidal gears
Tsuji et al. Machining and running test of high-performance face gear set
Piotrowski Hob identification methods
Shih et al. A flank correction methodology for hob sharpening on the five-axis cnc hob sharpening machine
CN113145943B (zh) 一种摆线轮加工用等前角剐齿刀的设计方法
CN114993230B (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