[go: up one dir, main page]

CN111912373A - A method for measuring tooth profile deviation using a roughness profiler - Google Patents

A method for measuring tooth profile deviation using a roughness profiler Download PDF

Info

Publication number
CN111912373A
CN111912373A CN202010670501.6A CN202010670501A CN111912373A CN 111912373 A CN111912373 A CN 111912373A CN 202010670501 A CN202010670501 A CN 202010670501A CN 111912373 A CN111912373 A CN 111912373A
Authority
CN
China
Prior art keywords
tooth profile
involute
deviation
gear
parameters
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
CN202010670501.6A
Other languages
Chinese (zh)
Other versions
CN111912373B (en
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.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN202010670501.6A priority Critical patent/CN111912373B/en
Publication of CN111912373A publication Critical patent/CN111912373A/en
Application granted granted Critical
Publication of CN111912373B publication Critical patent/CN111912373B/en
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
    • G01B21/20Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring contours or curvatures, e.g. determining profile
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/2416Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures of gears

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

The invention discloses a tooth profile deviation measuring method by using a roughness profiler, which comprises the steps of obtaining tooth profile data of an involute cylindrical gear on the roughness profiler, establishing an involute tooth profile model according to parameters of a measured gear, solving fitting parameters by using an optimization solution idea through constructing a least square objective function of original measured data and the involute tooth profile model in a normal direction, obtaining an original measured data curve and the involute tooth profile model after orthogonal distance fitting, further calculating to obtain the tooth profile deviation of any point in the involute normal direction, and finally evaluating and calculating the deviation value to obtain the tooth profile deviation and the precision grade defined by national standards.

Description

一种利用粗糙度轮廓仪的齿廓偏差测量方法A method for measuring tooth profile deviation using a roughness profiler

技术领域technical field

本发明属于精密测量领域,更具体的涉及渐开线圆柱齿轮齿廓偏差的测量与评定,尤其涉及一种渐开线齿廓测量数据的处理方法和一种渐开线齿廓偏差数据的评定计算方法。The invention belongs to the field of precision measurement, and more particularly relates to the measurement and evaluation of involute cylindrical gear tooth profile deviation, in particular to a processing method of involute tooth profile measurement data and an evaluation of involute tooth profile deviation data calculation method.

背景技术Background technique

齿轮是一种传动件,轮齿齿面的质量对齿轮传动的性能如传动误差、承载能力、振动噪声等有直接的影响。齿廓偏差是齿轮精度的重要参数,需要通过对齿轮齿廓信息进行评价得到。因此,对于齿轮齿廓信息的获取与处理显得极为重要。Gear is a kind of transmission part, and the quality of gear tooth surface has a direct influence on the performance of gear transmission, such as transmission error, bearing capacity, vibration and noise, etc. Tooth profile deviation is an important parameter of gear accuracy, which needs to be obtained by evaluating gear tooth profile information. Therefore, the acquisition and processing of gear tooth profile information is extremely important.

为了测量齿廓偏差,通常的做法是根据渐开线形成原理,通过机械或电子展成的方法形成理论渐开线,然后由测头记录实际齿廓跟理论渐开线的偏差,典型的仪器是齿轮测量中心。此外,不同于展成测量原理的三坐标测量机、光学测量仪器在齿轮测量中也有应用。In order to measure the tooth profile deviation, the usual practice is to form a theoretical involute by mechanical or electronic generation according to the principle of involute formation, and then record the deviation between the actual tooth profile and the theoretical involute by a probe. A typical instrument is the gear measurement center. In addition, CMM and optical measuring instruments, which are different from the principle of generative measurement, are also used in gear measurement.

粗糙度轮廓仪是一种接触式综合测量仪器,既可用于检测工件的二维形位误差,又被广泛的应用于表面粗糙度、波纹度、原始轮廓等微观轮廓参数的检测。Roughness profiler is a contact-type comprehensive measuring instrument, which can not only be used to detect the two-dimensional shape and position error of the workpiece, but also be widely used in the detection of surface roughness, waviness, original contour and other micro-profile parameters.

发明内容SUMMARY OF THE INVENTION

本发明提供了一种渐开线齿廓测量数据的处理方法,根据被测齿轮的参数建立渐开线齿廓模型,通过构造原始测量数据与渐开线齿廓模型在法向上的最小二乘目标函数,利用最优化求解的思想求解拟合参数,得到拟合后的原始测量数据曲线与渐开线齿廓模型。The invention provides a method for processing involute tooth profile measurement data. An involute tooth profile model is established according to the parameters of the measured gear. The objective function is to solve the fitting parameters by using the idea of optimal solution, and obtain the original measurement data curve and involute tooth profile model after fitting.

本发明提供了一种渐开线齿廓偏差数据的评定计算方法,根据求解得到的拟合参数,计算得到渐开线法向上的齿廓任意点偏差,将该偏差值进行评定计算,得到国家标准定义的齿廓偏差与精度等级。The invention provides an evaluation and calculation method for involute tooth profile deviation data. According to the obtained fitting parameters, the deviation of any point of the tooth profile in the normal direction of the involute is calculated, and the deviation value is evaluated and calculated to obtain the national Standard-defined profile deviation and accuracy class.

本发明通过在粗糙度轮廓仪上获取渐开线圆柱齿轮齿廓数据,根据被测齿轮的参数建立渐开线齿廓模型,通过构造原始测量数据与渐开线齿廓模型在法向上的最小二乘目标函数,利用最优化求解的思想求解拟合参数,得到拟合后的原始测量数据曲线与渐开线齿廓模型,进而计算得到渐开线法向上的齿廓任意点偏差值,最后将该偏差值进行评定计算,得到国家标准定义的齿廓偏差与精度等级。The invention obtains the tooth profile data of the involute cylindrical gear on the roughness profiler, establishes the involute tooth profile model according to the parameters of the measured gear, and constructs the minimum difference between the original measurement data and the involute tooth profile model in the normal direction. The square objective function is used to solve the fitting parameters by using the idea of optimal solution, and the fitted original measurement data curve and the involute tooth profile model are obtained, and then the deviation value of the tooth profile at any point on the involute normal direction is calculated. The deviation value is evaluated and calculated to obtain the tooth profile deviation and accuracy grade defined by the national standard.

附图说明Description of drawings

图1为本发明的三维结构等轴测示意图。FIG. 1 is a schematic isometric view of a three-dimensional structure of the present invention.

图2为本发明的结构侧视图。Figure 2 is a side view of the structure of the present invention.

图3为本发明的结构正视图。Figure 3 is a front view of the structure of the present invention.

图4为本发明的轮齿最佳测量位置示意图。FIG. 4 is a schematic diagram of the optimal measurement position of the gear teeth of the present invention.

图5为本发明的渐开线直角坐标示意图。5 is a schematic diagram of an involute rectangular coordinate of the present invention.

图6为本发明的预处理示意图。FIG. 6 is a schematic diagram of the preprocessing of the present invention.

图7为本发明的最小距离点示意图。FIG. 7 is a schematic diagram of the minimum distance point of the present invention.

图8为本发明的正交距离拟合结果示意图。FIG. 8 is a schematic diagram of an orthogonal distance fitting result of the present invention.

图9为本发明的渐开线齿廓偏差计算与评定示意图。FIG. 9 is a schematic diagram of the calculation and evaluation of the involute profile deviation of the present invention.

具体实施方式Detailed ways

本发明提供了一种在粗糙度轮廓仪上获取渐开线圆柱齿轮齿廓数据的方法。如图1所示,粗糙度轮廓仪包括立柱1、驱动箱2、测头系统3、轴向移动台4、底座5组成。如图2所示,渐开线齿轮夹具6包括电机7、联轴器8、主轴轴系9、圆光栅系统10、被测渐开线齿轮11。测头系统3由测量触针和传感器测杆组成,连接在驱动箱2上,使其可在X轴平面内(切向轴)沿被测齿面运动,立柱1可使测头系统3和驱动箱2在Z轴(垂直轴)做上下运动,渐开线齿轮夹具6实现对被测渐开线齿轮的装夹和X轴平面内的旋转,轴向移动台4实现渐开线齿轮夹具6在Y轴(径向轴)上的运动,使被测渐开线齿轮11在其齿宽方向上做测量位置调整,圆光栅系统10安装在渐开线齿轮夹具6上,用于检测被测渐开线齿轮11每次测量后的旋转角度。主轴轴系9左端通过联轴器8连接在电机7上,右端连接被测渐开线齿轮11,被测渐开线齿轮11的一端通过轴肩进行限位,另一端通过螺杆、螺母和轴套进行轴向的夹紧固定,主轴轴系9通过轴承支撑固定在轴向移动台4上。其中,电机轴、圆光栅、被测渐开线齿轮三者同步转动。The invention provides a method for acquiring tooth profile data of an involute cylindrical gear on a roughness profiler. As shown in Figure 1, the roughness profiler includes a column 1, a drive box 2, a probe system 3, an axial moving table 4, and a base 5. As shown in FIG. 2 , the involute gear fixture 6 includes a motor 7 , a coupling 8 , a main shaft shaft 9 , a circular grating system 10 , and an involute gear 11 to be measured. The probe system 3 is composed of a measuring stylus and a sensor rod, which is connected to the drive box 2 so that it can move along the measured tooth surface in the X-axis plane (tangential axis). The column 1 can make the probe system 3 and the The drive box 2 moves up and down on the Z axis (vertical axis), the involute gear clamp 6 realizes the clamping of the involute gear to be measured and the rotation in the X axis plane, and the axial moving table 4 realizes the involute gear clamp. 6 The movement on the Y axis (radial axis) makes the measured involute gear 11 adjust the measurement position in the direction of its tooth width, and the circular grating system 10 is installed on the involute gear fixture 6 to detect the Measure the rotation angle of the involute gear 11 after each measurement. The left end of the main shaft shaft 9 is connected to the motor 7 through the coupling 8, and the right end is connected to the measured involute gear 11. One end of the measured involute gear 11 is limited by the shaft shoulder, and the other end is limited by the screw, nut and shaft. The sleeve is clamped and fixed in the axial direction, and the main shaft shaft 9 is supported and fixed on the axial moving table 4 through the bearing. Among them, the motor shaft, the circular grating, and the measured involute gear rotate synchronously.

为了准确、全面地获取被测渐开线齿轮的渐开线齿廓,被测渐开线齿轮在安装后,需要进行测量位置的精确定位。待被测渐开线齿轮装夹完毕后,通过电机驱动被测轮齿旋转至最佳测量位置状态,使得被测轮齿齿面的齿根位置上的齿根成形点Ff与齿顶位置上的齿顶成形点Fa处在一个水平线上,如图4所示,使触针测头系统的测量变动范围降至最小,减少渐开线齿廓测量中的非线性误差。In order to obtain the involute tooth profile of the measured involute gear accurately and comprehensively, it is necessary to precisely locate the measurement position after the measured involute gear is installed. After the involute gear to be tested is clamped, the gear teeth to be tested are driven to rotate to the optimal measurement position by the motor, so that the root forming point F f on the tooth root position of the tooth surface of the tested gear and the position of the tooth tip are formed. The top forming point Fa is on a horizontal line, as shown in Figure 4, which minimizes the measurement variation range of the stylus probe system and reduces the nonlinear error in the measurement of the involute tooth profile.

吗待被测渐开线齿轮装夹完毕、齿面测量位置调整完毕后,进行齿廓数据的测量。测量的实施过程如下:第一步,控制轴向移动台,调整测头在轮齿齿廓的轴向中心位置上;第二步,将测头移动到被测轮齿的齿根成形点Ff附近,并继续朝齿根方向移动0.05mm,操作轮廓仪使测头自当前位置测量至超过齿顶成形点Fa附近0.05mm处,完成被测渐开线齿轮的齿廓测量,并将测量数据存储在计算机中。After the measured involute gear is clamped and the tooth surface measurement position is adjusted, the tooth profile data is measured. The implementation process of the measurement is as follows: the first step is to control the axial moving table, and adjust the probe to the axial center position of the tooth profile; the second step, move the probe to the tooth root forming point F of the measured tooth near f , and continue to move 0.05mm in the direction of the tooth root, operate the profiler to make the probe measure from the current position to 0.05mm beyond the tooth tip forming point Fa , complete the measurement of the tooth profile of the involute gear to be measured, and measure the tooth profile of the involute gear. The measurement data is stored in the computer.

本发明提供了一种渐开线齿廓测量数据的处理方法。第一步,根据被测渐开线齿轮的参数,基于渐开线展成法使用参数方程建模。渐开线是发生线绕基圆作纯滚动时,发生线上任意一点的轨迹。渐开线直角坐标如图5所示,根据此定义可得渐开线上任意一点K的滚动角uk为:The invention provides a method for processing involute tooth profile measurement data. In the first step, according to the parameters of the involute gear under test, a parametric equation is used for modeling based on the involute generation method. The involute is the trajectory of any point on the generating line when the generating line rolls purely around the base circle. The involute rectangular coordinates are shown in Figure 5. According to this definition, the rolling angle uk of any point K on the involute can be obtained as:

Figure BDA0002582112380000031
Figure BDA0002582112380000031

式中:θk为K点的展角,αk为K点的压力角。In the formula: θ k is the spread angle of point K, and α k is the pressure angle of point K.

渐开线齿廓参数方程表示为:The involute profile parameter equation is expressed as:

Figure BDA0002582112380000032
Figure BDA0002582112380000032

式中:rb为基圆半径。Where: r b is the radius of the base circle.

设最佳测量位置相对于式(2)中渐开线齿廓参数方程建模初始位置的旋转角为

Figure BDA0002582112380000039
通过对式(2)中渐开线齿廓参数方程进行旋转,最终得到最佳测量位置的渐开线齿廓理论模型:Set the rotation angle of the optimal measurement position relative to the initial position of the involute profile parameter equation modeling in Eq. (2) as
Figure BDA0002582112380000039
By rotating the involute tooth profile parameter equation in formula (2), the theoretical model of the involute tooth profile for the optimal measurement position is finally obtained:

Figure BDA0002582112380000033
Figure BDA0002582112380000033

式(3)中滚动角u的最小值umin与最大值umax分别取值于渐开线的起点与终点处。对于渐开线圆柱齿轮,渐开线的起点F点即为过渡曲线与渐开线齿廓的交点:In formula (3), the minimum value u min and the maximum value u max of the rolling angle u are respectively taken at the starting point and the ending point of the involute. For an involute cylindrical gear, the starting point F of the involute is the intersection of the transition curve and the involute tooth profile:

Figure BDA0002582112380000034
Figure BDA0002582112380000034

Figure BDA0002582112380000035
Figure BDA0002582112380000035

式中:αF为F点的压力角,αt为端面压力角,

Figure BDA0002582112380000036
为齿顶高系数,x为变位系数,z为齿数,dF为F点的直径,db为基圆直径。where α F is the pressure angle at point F, α t is the end face pressure angle,
Figure BDA0002582112380000036
is the addendum height coefficient, x is the displacement coefficient, z is the number of teeth, d F is the diameter of point F, and db is the diameter of the base circle.

若忽略齿顶倒角,渐开线的终点在齿顶圆上。因此,滚动角u的最小值umin与最大值umax的计算公式为:If the addendum chamfer is ignored, the end point of the involute is on the addendum circle. Therefore, the calculation formula of the minimum value u min and the maximum value u max of the roll angle u is:

Figure BDA0002582112380000037
Figure BDA0002582112380000037

Figure BDA0002582112380000038
Figure BDA0002582112380000038

式中:da为齿顶圆直径。Where: d a is the diameter of the addendum circle.

第二步,测量数据预处理。为了保证算法的稳定性,同时提高算法的收敛速度,预处理使渐开线齿廓模型与原始测量数据曲线在正交距离拟合前快速尽可能地靠近。The second step is measurement data preprocessing. In order to ensure the stability of the algorithm and improve the convergence speed of the algorithm, the preprocessing makes the involute tooth profile model and the original measurement data curve as fast as possible before the orthogonal distance fitting.

由轮廓仪对齿轮齿廓测量所得的原始测量数据不具有位置信息,只具有长度、形状信息。预处理在x方向上,将原始测量数据与渐开线齿廓模型两者的中点位置平移对齐。在y方向上,通过计算渐开线齿廓模型与原始测量数据在y方向上两者的数据均值差,将原始测量数据与渐开线齿廓模型对齐,预处理后的原始测量数据为:The original measurement data obtained by measuring the gear tooth profile by the profiler has no position information, but only length and shape information. Preprocessing In the x-direction, the raw measurement data and the midpoint position of the involute tooth profile model are aligned in translation. In the y direction, by calculating the mean difference between the involute tooth profile model and the original measurement data in the y direction, the original measurement data is aligned with the involute tooth profile model, and the preprocessed original measurement data is:

Figure BDA0002582112380000041
Figure BDA0002582112380000041

如图6所示,通过预处理,最大可能地使得渐开线齿廓模型与原始测量数据曲线的位置更加接近。As shown in Fig. 6, through preprocessing, the position of the involute tooth profile model and the original measurement data curve can be made closer to the greatest possible extent.

第三步,最优化求解。正交距离拟合算法(Orthogonal Distance FittingAlgorithm)需要求解预处理后的测量数据曲线与渐开线齿廓模型在法向上的最小二乘解,通过构造目标函数利用最优化求解的思想,求解拟合参数,得到预处理后的测量数据曲线与渐开线齿廓模型的正交距离拟合结果。The third step is to optimize the solution. The Orthogonal Distance Fitting Algorithm needs to solve the least squares solution of the preprocessed measurement data curve and the involute tooth profile model in the normal direction. By constructing the objective function, the idea of optimization is used to solve the fitting parameters to obtain the fitting result of the orthogonal distance between the preprocessed measurement data curve and the involute tooth profile model.

首先利用列文伯格-马夸尔特算法(Levenberg-Marquardt Algorithm),求解预处理后的每个测量数据点对应渐开线齿廓模型上的最小距离点,得到所有测量点对应渐开线齿廓模型上所有最小距离点的位置参数u。如图4所示,Mi为预处理后测量数据中的任意点(i=1,2,...,n,n为测量数据总个数),Ti为Mi在渐开线齿廓模型上的最小距离点,即正交距离对应点,x(u)为渐开线齿廓模型,di即为Mi与最小距离点Ti的法向距离。First, use the Levenberg-Marquardt Algorithm to solve the minimum distance point on the involute tooth profile model corresponding to each measurement data point after preprocessing, and obtain all measurement points corresponding to the involute curve The position parameter u of all minimum distance points on the tooth profile model. As shown in Figure 4, M i is any point in the pre-processed measurement data (i=1, 2,..., n, n is the total number of measurement data), and T i is M i in the involute tooth The minimum distance point on the profile model, that is, the corresponding point of the orthogonal distance, x(u) is the involute tooth profile model, and d i is the normal distance between Mi and the minimum distance point Ti .

设法向距离为:The managed distance is:

Figure BDA0002582112380000042
Figure BDA0002582112380000042

求解最小距离点Ti的位置参数u即可转换为对目标函数D(u)极值的求解:Solving the position parameter u of the minimum distance point T i can be converted into solving the extreme value of the objective function D(u):

Figure BDA0002582112380000043
Figure BDA0002582112380000043

L-M算法为高斯-牛顿迭代法(Gauss-Newtoniteration method)的改进,引入阻尼系数λ使得迭代具有更大的收敛区间,由此得到迭代公式:The L-M algorithm is an improvement of the Gauss-Newton iteration method. The damping coefficient λ is introduced to make the iteration have a larger convergence interval, and the iteration formula is obtained:

Figure BDA0002582112380000044
Figure BDA0002582112380000044

L-M算法实现所有预处理后的测量数据点对渐开线齿廓模型上最小距离点位置参数u求解,所求解的所有位置参数u用作高斯-牛顿迭代法所有位置参数u′的迭代初值,同时加入旋转参数

Figure BDA0002582112380000045
和平移参数x0、y0,通过高斯-牛顿迭代法最优化求解,实现预处理后的测量数据曲线与渐开线齿廓模型的正交距离拟合,并求解出相关拟合参数(旋转参数
Figure BDA0002582112380000046
平移参数x0、y0,所有位置参数u′):The LM algorithm realizes that all preprocessed measurement data points solve the position parameter u of the minimum distance point on the involute tooth profile model, and all the solved position parameters u are used as the initial iterative value of all position parameters u' of the Gauss-Newton iteration method. , while adding the rotation parameter
Figure BDA0002582112380000045
and translation parameters x 0 , y 0 are optimized and solved by Gauss-Newton iterative method to realize the orthogonal distance fitting between the preprocessed measurement data curve and the involute tooth profile model, and solve the relevant fitting parameters (rotation). parameter
Figure BDA0002582112380000046
Translation parameters x 0 , y 0 , all positional parameters u′):

Figure BDA0002582112380000047
Figure BDA0002582112380000047

设旋转和平移处理后的最小距离对应点为:Let the minimum distance corresponding point after rotation and translation processing be:

T′=R-1M+X0 (13)T'=R -1 M+X 0 (13)

基于正交距离拟合的思想,需要求解旋转和平移处理后的渐开线齿廓模型与预处理后的测量数据之间最小距离的平方和最小化,且预处理后的测量数据给定点与渐开线齿廓模型之间的每个距离也应最小化。旋转、平移处理后的渐开线齿廓模型与预处理后的测量数据曲线的距离平方和为:Based on the idea of orthogonal distance fitting, it is necessary to minimize the sum of the squares of the minimum distances between the involute tooth profile model after rotation and translation processing and the preprocessed measurement data, and the given point of the preprocessed measurement data is equal to Each distance between the involute tooth profile models should also be minimized. The sum of the squares of the distances between the involute tooth profile model after rotation and translation processing and the preprocessed measurement data curve is:

Figure BDA0002582112380000051
Figure BDA0002582112380000051

一阶必要条件为:The first-order necessary conditions are:

Figure BDA0002582112380000052
Figure BDA0002582112380000052

基于高斯-牛顿迭代法的思想,可得迭代公式:Based on the idea of Gauss-Newton iteration method, the iterative formula can be obtained:

J|kΔb=(M-T′)|k,bk+1=bk+αΔb (16)J| k Δb=(MT′)| k , b k+1 =b k +αΔb (16)

迭代公式(16)的展开形式为:The expanded form of iterative formula (16) is:

Figure BDA0002582112380000053
Figure BDA0002582112380000053

式(17)是关于Δb的线性超定方程组,需求最小二乘解。Equation (17) is a system of linear overdetermined equations for Δb that requires a least-squares solution.

设定停止条件为:Set the stop condition as:

|bk+1-bk|<ε (18)|b k+1 -b k |<ε (18)

对于需要求解的拟合参数b,最终求解结果中的

Figure BDA0002582112380000054
决定着拟合的精度,进而影响着齿廓偏差计算的结果。由此,迭代的停止条件为:For the fitting parameter b that needs to be solved, in the final solution result
Figure BDA0002582112380000054
It determines the accuracy of fitting, and then affects the result of tooth profile deviation calculation. Thus, the stopping condition of the iteration is:

Figure BDA0002582112380000055
Figure BDA0002582112380000055

Figure BDA0002582112380000056
Figure BDA0002582112380000056

Figure BDA0002582112380000057
Figure BDA0002582112380000057

求解出的拟合参数b包含正交距离拟合的最佳旋转参数

Figure BDA0002582112380000058
平移参数x0、y0和所有位置参数u′。为了保留测量齿廓在整个齿轮上的相对位置,最终的拟合结果对渐开线齿廓模型旋转
Figure BDA0002582112380000059
对原始测量数据平移x0、y0,最终得到渐开线齿廓模型与预处理后的测量数据的正交距离拟合结果,如图8所示(偏差放大)。The solved fitting parameter b contains the optimal rotation parameter for the orthogonal distance fitting
Figure BDA0002582112380000058
Translation parameters x 0 , y 0 and all positional parameters u'. In order to preserve the relative position of the measured tooth profile on the entire gear, the final fitting result is rotated on the involute tooth profile model
Figure BDA0002582112380000059
The original measurement data is shifted by x 0 and y 0 , and finally the orthogonal distance fitting result between the involute tooth profile model and the preprocessed measurement data is obtained, as shown in Figure 8 (bias magnification).

本发明提供了一种渐开线齿廓偏差数据的评定计算方法。第一步,根据ODF算法得到的所有位置参数u计算得到ODF处理后的齿廓法向偏差:The invention provides a method for evaluating and calculating involute tooth profile deviation data. The first step is to calculate the normal deviation of the tooth profile after ODF processing according to all the position parameters u obtained by the ODF algorithm:

Figure BDA00025821123800000510
Figure BDA00025821123800000510

第二步,根据式(22)计算得到ODF处理后的齿廓法向偏差结果,基于圆柱齿轮精度制现行国家标准GB/T10095.1-2008中的定义进行计算与评价。其中,齿廓偏差计算根据圆柱齿轮精度制现行国家标准GB/T10095.1-2008中有关齿廓偏差的规定:齿廓总偏差(Fα)表示在计值范围Lα内,包容实际齿廓迹线的两条设计齿廓迹线间的距离;齿廓形状偏差

Figure BDA0002582112380000061
表示在计值范围Lα内,包容实际齿廓迹线的,与平均齿廓迹线完全相同的两条迹线间的距离,且两条曲线与平均齿廓迹线的距离为常数;齿廓倾斜偏差
Figure BDA0002582112380000062
表示在计值范围Lα内,两端与平均齿廓迹线相交的两条设计齿廓迹线间的距离,如图9所示。In the second step, the normal deviation result of the tooth profile after ODF processing is calculated according to formula (22), and the calculation and evaluation are carried out based on the definition in the current national standard GB/T10095.1-2008 for the precision manufacturing of cylindrical gears. Among them, the tooth profile deviation is calculated according to the regulations on tooth profile deviation in the current national standard GB/T10095.1-2008 for the precision of cylindrical gears: the total tooth profile deviation (F α ) is expressed in the calculation range L α , including the actual tooth profile The distance between the two design tooth profile traces of the trace; tooth profile shape deviation
Figure BDA0002582112380000061
Represents the distance between two traces that are exactly the same as the average tooth profile trace within the value range L α , including the actual tooth profile trace, and the distance between the two curves and the average tooth profile trace is constant; Profile tilt deviation
Figure BDA0002582112380000062
Indicates the distance between two designed tooth profile traces whose ends intersect with the average tooth profile trace within the value range L α , as shown in Figure 9.

第三步,根据计算结果,基于圆柱齿轮精度制现行国家标准GB/T10095.1-2008查表得到齿廓总偏差(Fα)、齿廓形状偏差

Figure BDA0002582112380000063
齿廓倾斜偏差
Figure BDA0002582112380000064
三个项目的单项精度,完成对被测渐开线齿轮的精度评价。In the third step, according to the calculation results, based on the precision of cylindrical gears, the current national standard GB/T10095.1-2008 is made to look up the table to obtain the total tooth profile deviation (F α ), tooth profile shape deviation
Figure BDA0002582112380000063
Tooth profile inclination deviation
Figure BDA0002582112380000064
The single precision of the three items completes the precision evaluation of the measured involute gear.

Claims (2)

1. A tooth profile deviation measuring method using a roughness profiler is characterized in that: establishing an involute tooth profile model according to parameters of a measured gear, and solving fitting parameters by using an optimization solution idea through constructing a least square objective function of original measurement data and the involute tooth profile model in a normal direction to obtain a fitted original measurement data curve and the involute tooth profile model;
and calculating to obtain the deviation of any point of the tooth profile in the normal direction of the involute according to the fitting parameters obtained by solving, and evaluating and calculating the deviation value to obtain the tooth profile deviation and the precision grade defined by the national standard.
2. The method of tooth profile deviation measurement using a roughness profiler as set forth in claim 1, wherein: the method comprises the steps of obtaining tooth profile data of an involute cylindrical gear on a roughness profiler, establishing an involute tooth profile model according to parameters of a measured gear, constructing a least square objective function of original measured data and the involute tooth profile model in a normal direction, solving fitting parameters by using an optimization solving idea to obtain a fitted original measured data curve and the involute tooth profile model, further calculating to obtain tooth profile arbitrary point deviation in the involute normal direction, and finally evaluating and calculating the deviation value to obtain tooth profile deviation and precision grade defined by national standards.
CN202010670501.6A 2020-07-13 2020-07-13 A method for measuring tooth profile deviation using a roughness profiler Active CN111912373B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010670501.6A CN111912373B (en) 2020-07-13 2020-07-13 A method for measuring tooth profile deviation using a roughness profiler

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010670501.6A CN111912373B (en) 2020-07-13 2020-07-13 A method for measuring tooth profile deviation using a roughness profiler

Publications (2)

Publication Number Publication Date
CN111912373A true CN111912373A (en) 2020-11-10
CN111912373B CN111912373B (en) 2021-11-12

Family

ID=73227933

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010670501.6A Active CN111912373B (en) 2020-07-13 2020-07-13 A method for measuring tooth profile deviation using a roughness profiler

Country Status (1)

Country Link
CN (1) CN111912373B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113283025A (en) * 2021-05-11 2021-08-20 北京理工大学 Involute tooth profile error modeling method containing system error
CN113566772A (en) * 2021-07-19 2021-10-29 北京工业大学 Local tooth surface positioning method based on coordinate measurement
CN115031678A (en) * 2022-06-13 2022-09-09 北京工业大学 Noise gear screening method based on tooth profile waviness information and gear transmission error information
CN115143922A (en) * 2022-06-19 2022-10-04 北京工业大学 Gear tooth surface roughness measurement and evaluation method for transmission performance analysis
CN116429047A (en) * 2023-05-04 2023-07-14 扬州保来得科技实业有限公司 Gear profile measuring and evaluating method
CN116518907A (en) * 2023-05-05 2023-08-01 北京工业大学 A gear tooth profile roughness measurement method, system, equipment and medium
CN119268637A (en) * 2024-12-06 2025-01-07 昆山奥德鲁自动化技术有限公司 A waviness detection device
CN119293396A (en) * 2024-10-14 2025-01-10 昆明理工大学 A method for measuring arbitrary tooth thickness deviation considering installation eccentricity

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0131537A1 (en) * 1983-06-09 1985-01-16 BHS-Höfler Maschinenbau GmbH Device for checking an involute gear tooth profile
EP0510431A2 (en) * 1991-04-22 1992-10-28 Toyota Jidosha Kabushiki Kaisha Method and apparatus for evaluating gear motion characteristics, based on tooth profile deflection differentiated by rotation angle of the gear
CN102216726A (en) * 2009-04-24 2011-10-12 株式会社东京技术 Method of measuring an involute gear tooth profile
CN107588742A (en) * 2017-10-25 2018-01-16 北京工业大学 A kind of cylindrical gear profile bias measurement method based on line-structured light
CN110455244A (en) * 2019-08-28 2019-11-15 北京工业大学 A Method of Measuring Pitch Deviation of Cylindrical Gears Using Roughness Profiler
CN111060061A (en) * 2020-01-13 2020-04-24 北京工业大学 Involute gear sample plate tooth profile measuring and correcting method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0131537A1 (en) * 1983-06-09 1985-01-16 BHS-Höfler Maschinenbau GmbH Device for checking an involute gear tooth profile
EP0510431A2 (en) * 1991-04-22 1992-10-28 Toyota Jidosha Kabushiki Kaisha Method and apparatus for evaluating gear motion characteristics, based on tooth profile deflection differentiated by rotation angle of the gear
CN102216726A (en) * 2009-04-24 2011-10-12 株式会社东京技术 Method of measuring an involute gear tooth profile
CN107588742A (en) * 2017-10-25 2018-01-16 北京工业大学 A kind of cylindrical gear profile bias measurement method based on line-structured light
CN110455244A (en) * 2019-08-28 2019-11-15 北京工业大学 A Method of Measuring Pitch Deviation of Cylindrical Gears Using Roughness Profiler
CN111060061A (en) * 2020-01-13 2020-04-24 北京工业大学 Involute gear sample plate tooth profile measuring and correcting method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
林虎等: "基于正交距离回归齿面的齿轮误差评定", 《光学精密工程》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113283025B (en) * 2021-05-11 2022-09-06 北京理工大学 A Modeling Method of Involute Profile Error Including System Error
CN113283025A (en) * 2021-05-11 2021-08-20 北京理工大学 Involute tooth profile error modeling method containing system error
CN113566772A (en) * 2021-07-19 2021-10-29 北京工业大学 Local tooth surface positioning method based on coordinate measurement
CN115031678A (en) * 2022-06-13 2022-09-09 北京工业大学 Noise gear screening method based on tooth profile waviness information and gear transmission error information
CN115143922A (en) * 2022-06-19 2022-10-04 北京工业大学 Gear tooth surface roughness measurement and evaluation method for transmission performance analysis
CN115143922B (en) * 2022-06-19 2023-09-26 北京工业大学 Gear tooth surface roughness measurement and evaluation method for transmission performance analysis
CN116429047B (en) * 2023-05-04 2023-10-13 扬州保来得科技实业有限公司 Gear profile measuring and evaluating method
CN116429047A (en) * 2023-05-04 2023-07-14 扬州保来得科技实业有限公司 Gear profile measuring and evaluating method
CN116518907A (en) * 2023-05-05 2023-08-01 北京工业大学 A gear tooth profile roughness measurement method, system, equipment and medium
CN116518907B (en) * 2023-05-05 2025-07-25 北京工业大学 Gear tooth profile roughness measurement method, system, equipment and medium
CN119293396A (en) * 2024-10-14 2025-01-10 昆明理工大学 A method for measuring arbitrary tooth thickness deviation considering installation eccentricity
CN119293396B (en) * 2024-10-14 2025-05-13 昆明理工大学 Method for measuring arbitrary tooth thickness deviation by considering installation eccentricity
CN119268637A (en) * 2024-12-06 2025-01-07 昆山奥德鲁自动化技术有限公司 A waviness detection device

Also Published As

Publication number Publication date
CN111912373B (en) 2021-11-12

Similar Documents

Publication Publication Date Title
CN111912373A (en) A method for measuring tooth profile deviation using a roughness profiler
CN105423946B (en) Axle journal axle center measuring device based on laser displacement sensor and measurement scaling method
CN108050946B (en) A method for measuring gear tooth thickness based on line structured light
CN107588742A (en) A kind of cylindrical gear profile bias measurement method based on line-structured light
CN111060061B (en) A method for measuring and correcting tooth profile of an involute gear template
CN108871229A (en) A kind of measurement method of ball nut spiral interior rollaway nest curved surface and outer diameter
CN110440723A (en) A kind of abnormally-structured part negative camber measuring device for surface roughness and measurement method
CN115307571B (en) A Pose Calibration Part and Calibration Method for a Planar Line Laser Sensor
LIN et al. Tooth profile deviation measurement of cylindrical gears based on a roughness profilometer
CN106796095B (en) Method for operating a coordinate measuring device, coordinate measuring device and computer program
Garratt A new stylus instrument with a wide dynamic range for use in surface metrology
CN115540730A (en) A Coordinate Measurement System and Method for High Steep or Deep Concave Complex Surface
CN216846033U (en) Inner wall measurement system based on deep sagittal workpiece
Hill et al. Computer Aided Error Analysis Of Three-dimensional Precision Surface Mapping System
CN100491895C (en) Three-coordinate system calibration tester
CN112344863B (en) Wall thickness detection method for free-form surface revolving body workpiece
Zapico et al. Extrinsic calibration of a conoscopic holography system integrated in a CMM
CN113310439A (en) Method and system for detecting surface roughness of spherical shell part based on white light confocal sensor
CN115752294B (en) Method for measuring three-dimensional surface profile of large complex shaft of aero-engine
CN111060056A (en) Reconstruction device and reconstruction method for accurately reconstructing parallel contour
CN205138446U (en) Axle journal axle center measuring device based on laser displacement sensor
CN115307572B (en) A fixed-width planar line laser sensor pose calibration device and calibration method
CN115307570B (en) U-groove planar line laser sensor pose calibration piece and calibration method
Zhang et al. The development of cylindrical coordinate measuring machines
CN113310455B (en) Method for detecting eccentricity of compression molding processing double-sided mirror

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