WO2016110236A1 - 齿轮动态传动误差确定方法 - Google Patents

齿轮动态传动误差确定方法 Download PDF

Info

Publication number
WO2016110236A1
WO2016110236A1 PCT/CN2016/070104 CN2016070104W WO2016110236A1 WO 2016110236 A1 WO2016110236 A1 WO 2016110236A1 CN 2016070104 W CN2016070104 W CN 2016070104W WO 2016110236 A1 WO2016110236 A1 WO 2016110236A1
Authority
WO
WIPO (PCT)
Prior art keywords
relative displacement
time
gear
collision
time interval
Prior art date
Application number
PCT/CN2016/070104
Other languages
English (en)
French (fr)
Inventor
刘富豪
Original Assignee
盐城工学院
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 盐城工学院 filed Critical 盐城工学院
Priority to US15/319,228 priority Critical patent/US10274399B2/en
Publication of WO2016110236A1 publication Critical patent/WO2016110236A1/zh

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00—Testing of machine parts
    • G01M13/02—Gearings; Transmission mechanisms
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00—Testing of machine parts
    • G01M13/02—Gearings; Transmission mechanisms
    • G01M13/021—Gearings

Definitions

  • the invention relates to the field of mechanical manufacturing and control, and in particular to a method for determining a dynamic gear transmission error of a gear considering a flank clearance.
  • gear transmission has the characteristics of smooth transmission, accurate transmission ratio, high efficiency and long service life, making gear transmission the most important and widely used mechanical transmission.
  • High-speed, stable, precision gear transmission is a national key scientific and technological project.
  • Transmission error has been widely accepted and considered to be an excitation source for vibration and noise of the gear system because the transmission error affects the force and speed changes in the gear transmission.
  • the dynamic transmission error of the gear is also a very important part of gear detection. Accurate testing and calculation of transmission error have great significance for gear manufacturing and gear dynamics research.
  • the name is the small-modulus gear transmission error measurement method, which firstly measures the gear, the actuator and the first angle measuring encoder.
  • the first coupling and the first motor are sequentially connected, and then the measuring gear, the flexible coupling, the second goniometer encoder, the second coupling and the second motor are sequentially connected.
  • the measured gear and the measuring gear mesh on one side and are driven by two motors respectively.
  • the angular displacement ⁇ 2 of the measuring gear relative to the second angle measuring encoder, the measured gear shaft angular displacement ⁇ 1 and the measured gear shaft angular displacement ⁇ ' are detected, and the transmission error of the small modulus gear is calculated.
  • the numerical calculation of the gear system transmission error has become an increasingly prominent problem.
  • the calculation method of the transmission error mostly stays in the existing conventional methods, such as: Newmark integration method, Runge-Kutta integration method, Gear method and improved Gear method (Wstiff method and Dstiff method).
  • the above calculation methods have different emphasis on the calculation method of the transmission error. Some methods focus more on the accuracy of calculation; some methods have faster calculation speed.
  • Traditional calculation methods usually cannot perfectly control the distribution of calculation time and calculation accuracy: in the case of high precision requirements, it often brings a huge amount of calculation, resulting in a long calculation time; under the condition of ensuring calculation time, The calculation accuracy is so degraded that the gear system dynamic transmission error cannot be properly described and predicted.
  • the invention mainly solves the technical problem that the measurement error of the transmission error is caused by the flank clearance in the actual measurement process in the prior art, and it is difficult to balance the calculation accuracy and the calculation time in the theoretical simulation calculation of the gear dynamic transmission error.
  • a calculation method that combines the actual measurement with the theoretical simulation calculation and can quickly calculate the gear dynamic transmission error under the condition of ensuring a certain calculation accuracy is proposed.
  • a gear dynamic transmission error determination method comprising the following steps:
  • A. Obtain a parameter value, the input shaft angular velocity measured by the angle encoder and the output shaft resistance torque provided by the magnetic powder brake;
  • step C comparing the collision time interval ⁇ t i and the time step size ⁇ t, if ⁇ t ⁇ ⁇ t i , then proceeds to step D; if ⁇ t ⁇ ⁇ t i , then proceeds to step G;
  • step D Calculate the relative displacement x(n+1) by using the time step ⁇ t as a step size. If the absolute value
  • Step E comparing the time step ⁇ t and the minimum time step t min , if the time step ⁇ t is greater than or equal to the minimum time step t min , proceeding to step D; if the time step ⁇ t is less than the minimum time step t min , then entering Step F;
  • step J Calculate the relative displacement with the minimum time step t min as the step size Then enter step J;
  • step J If the absolute value
  • step J determine whether t(n+1) is greater than M*T f , if yes, the calculation is ended, if not, the step number n is increased by 1 and then returns to step B;
  • ⁇ is a small error
  • M is the number of cycles
  • T f is the excitation period
  • t(n) is the nth time sampling point
  • x(n) is the nth relative displacement
  • the origin above the symbol indicates the derivation of time.
  • One dot represents the first derivative
  • two dots represent the second derivative
  • sign is a sign function.
  • the relative displacement x(n+1) and relative velocity of the resulting M cycles And the time t(n+1) can well reflect the dynamic transmission error of the gear.
  • the accuracy of the gearbox can be evaluated by the resulting transmission error.
  • the values of each parameter are measured by the conventional method or set by the user; when the simulation is calculated, the parameter value is manually set or automatically set by the computer.
  • the parameter values include L, ⁇ , T f , R bg , R bp , ⁇ g , ⁇ p , the mean value of the external load moment, the moment of inertia of the driven wheel, the fluctuation part of the external load moment, ⁇ , t 0 , due to the oil film.
  • the collision time interval is calculated as follows:
  • ⁇ t i takes the smallest positive real number of the two values
  • R bg is the base circle radius of the driven wheel
  • R bp is the base circle radius of the driving wheel
  • T m is the ratio of the mean value of the external load moment to the moment of inertia of the driven wheel
  • T p The ratio of the fluctuation of the external load moment to the moment of inertia of the driven wheel
  • ⁇ is the frequency of the external load moment
  • t 0 is the system start time
  • k is the sum of the meshing stiffness of the front and side of the gear pair contact due to the oil film
  • c is the ratio of the sum of the meshing damping of the front and side of the gear pair contact due to the oil film to the moment of inertia of the driven wheel.
  • ⁇ t i when the collision event interval ⁇ t i is calculated, if ⁇ t i has no positive real number, ⁇ t i ⁇ t is directly determined.
  • step D the specific operation of reducing the time step is to reduce ⁇ t by ⁇ t 2 .
  • step D to reduce the impact time interval specific operation is: to reduce the decrease ⁇ t i ⁇ t i 2.
  • the substantial effect brought by the invention is that the influence of the flank clearance on the transmission error is considered, and under the requirement of ensuring a certain calculation accuracy, the calculation time is minimized, and a more optimal allocation of calculation accuracy and calculation time is realized.
  • Figure 1 is an enlarged view of a gear meshing area of the present invention
  • FIG. 2 is a schematic diagram of a parameter measurement system of the present invention
  • FIG. 3 is a schematic view showing the assembly of a driving wheel and a driven wheel according to a standard center distance according to the present invention
  • Figure 4 is a schematic view showing the judgment of the driving wheel and the driven wheel in the case of the adhesion of the driving wheel;
  • Figure 5 is a schematic view showing a state in which ⁇ t ⁇ ⁇ t i when a collision between teeth and teeth occurs in the present invention
  • Figure 6 is a schematic view of the present invention using ⁇ to make a collision between the teeth again accurate
  • Figure 7 is a flow chart of a determining method of the present invention.
  • FIG. 2 is a system for measuring an actual gear set, including a motor 1, a drive wheel 2, an input shaft 3, an angle encoder 4, a driven wheel 5, an output shaft 6, and a magnetic powder brake 7.
  • the driving wheel 2 and the driven wheel 5 are assembled according to a standard center distance as shown in FIG. 3, and
  • FIG. 1 is an enlarged view of the gear meshing area.
  • the gear transmission error (relative displacement) is the difference between the position of the tooth profile on the passive wheel in the actual meshing position and the position under the ideal condition along the direction of the meshing line, which can be defined as:
  • R bg is the base circle radius of the driven wheel
  • R bp is the base circle radius of the driving wheel
  • ⁇ g is the angular displacement of the driven wheel
  • ⁇ p is the angular displacement of the driving wheel
  • the angular displacement can also be performed by the angle encoder measuring.
  • Angle encoder measure the angular velocity of the input shaft in Figure 2.
  • Magnetic powder brake Provides a resistance torque to the output shaft.
  • the resistance torque is:
  • Ig is the moment of inertia of the large gear
  • k 1 and c 1 are the meshing stiffness and meshing damping of the gear pair contact front due to the oil film, respectively
  • k 2 and c 2 are the contact side of the gear pair due to the oil film, respectively Engagement stiffness and meshing damping; with They are the mean of the external load moment and the moment of the fluctuating part.
  • the transmission error can be expressed as:
  • ⁇ t i takes only the smallest positive real number. If ⁇ t i has no minimum positive real number, ⁇ t ⁇ ⁇ t i is considered.
  • ⁇ t i is compared with ⁇ t to determine what the next time step is. If ⁇ t ⁇ ⁇ t i , the relative displacement x(n+1) is calculated in time step ⁇ t. If
  • ⁇ t i is used as the next step to calculate the next relative displacement and velocity.
  • >L+ ⁇ we have to judge x(n+1). If
  • the formula for decreasing ⁇ t i is as follows:
  • the gear dynamic transmission error determining method of the present invention passes the input shaft 3 angular velocity measured by the angle encoder 4, the output shaft 6 resistance torque provided by the magnetic powder brake 7, the flank clearance, the input shaft 3, and the output.
  • the angular displacement of the shaft 6 and other parameter values calculate the collision time interval, and the obtained relative displacement is compared with the flank clearance until the relative displacement value is within the appropriate interval and meets the requirements, and the obtained result is higher precision.
  • Transmission error introduces the influence of the flank clearance on the transmission error, ensures the accuracy of the determination result of the transmission error, and achieves a good balance between accuracy and time, and is suitable for the calculation calculation and simulation calculation of the actual gear set.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Gear Transmission (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Gears, Cams (AREA)

Abstract

一种齿轮动态传动误差确定方法，通过由角度编码仪（4）测得的输入轴（3）角速度、由磁粉制动器（7）提供的输出轴（6）阻力矩、齿侧间隙、输入轴（3）和输出轴（6）的角位移等参数值计算出碰撞时间间隔，将所得到的相对位移与齿侧间隙进行比较，直至相对位移的值在适当的区间内并且符合要求，所得到的结果即为精度较高的传动误差。通过考虑齿侧间隙对传动误差的影响，保证了传动误差确定结果的准确性，同时使精度和时间得到了良好的平衡，适用于对实际齿轮组的测量计算和仿真模拟的计算。

Description

齿轮动态传动误差确定方法 技术领域
本发明涉及机械制造和控制领域，尤其是涉及一种考虑齿侧间隙的齿轮动态传动误差确定方法。
背景技术
在所有的机械传动中，齿轮传动有着传动平稳、传动比精确、效率高、寿命长等特点，使得齿轮传动成为最重要而且应用最广的机械传动。高速，平稳，精密齿轮传动是国家重点科技攻关项目。传动误差已经被人们普遍的接受和认为是齿轮系统振动和噪声的激励源，因为传动误差影响着齿轮传动中的力和速度的变化。同时，齿轮的动态传动误差也是齿轮检测中十分重要的一项内容，精确测试和计算传动误差对齿轮制造和齿轮动力学研究具有重大意义。
中华人民共和国国家知识产权局于2008年10月01公开了公告号为CN101275881的专利文献，名称是小模数齿轮传动误差测量方法，其先将被测齿轮、带动器、第一测角编码器、第一联轴器和第一电机依次连接，再将测量齿轮、柔性联轴器、第二测角编码器、第二联轴器和第二电机依次连接。被测齿轮和测量齿轮单面啮合，并分别由两个电机进行驱动。检测出测量齿轮相对于第二测角编码器的角位移φ2、测量齿轮轴系角位移φ1和被测齿轮轴系角位移φ′，经过计算得到小模数齿轮的传动误差。但是，由于齿轮的加工误差，安装误差，润滑和修行等因素使得齿廓表面相对于理想齿廓位置的偏移会引起齿侧间隙(图1)，导致齿轮在进入和退出啮合时，啮入和啮出点的位置会偏离理论啮合点，产生线外啮合，使啮合齿面间产生碰撞冲击。这种碰撞冲击对高速齿轮传动的平稳性和精密性起着关键的作用。对于精密传动而言，正反转的时候会有微量空行程，而且有冲击，影响传动精度。显然，齿侧间隙直接影响着 齿轮系统得传动误差。然而，此方案未考虑齿侧间隙对传动误差的影响，会使最终的计算结果与实际值发生较大偏差。
同时，对于齿轮系统传动误差的数值计算成为了一个日益突出的问题。目前，对于传动误差的计算方法，多停留在应用现有的传统方法，如：Newmark积分法，Runge-Kutta积分法，Gear法和改进的Gear法(Wstiff法和Dstiff法)。上述计算方法对于传动误差的计算方法的侧重点各有不同，有些方法更侧重于对计算精度的把握；也有一些方法拥有较快的计算速度。传统的计算方法通常无法完美的控制好计算时间与计算精度的分配：在精度要求较高的情况下往往带来庞大的计算量，导致计算时间的漫长；在保证计算时间的条件下，又会使计算精度下降以致于无法正确的描述及预测齿轮系统动传动误差。
发明内容
本发明主要是解决现有技术在实际测量过程中因齿侧间隙导致传动误差测量不准确的技术问题，及在理论模拟计算齿轮动态传动误差难以平衡计算精度和计算时间的问题。在考虑到齿侧间隙对于传动误差影响的前提下，提出了一种将实际测量和理论模拟计算相结合，并能在保证一定的计算精度情况下快速计算齿轮动态传动误差的计算方法。
本发明针对上述技术问题主要是通过下述技术方案得以解决的：一种齿轮动态传动误差确定方法，包括以下步骤：
A、获取参数值，所述参数值包括由角度编码仪测得的输入轴角速度和由磁粉制动器提供的输出轴阻力矩；
B、计算两个轮齿之间的碰撞时间间隔Δti；
C、比较碰撞时间间隔Δti和时间步长Δt的大小，如果Δt＜Δti，则进入步骤D；如果Δt≥Δti，则进入步骤G；
D、以时间步长Δt为步长来计算相对位移x(n+1)，如果相对位移x(n+1)的绝对值|x(n+1)|小于或等于齿侧间隙L，则进入步骤J；如果相对位移x(n+1)的 绝对值|x(n+1)|大于齿侧间隙L，则缩小时间步长Δt然后进入步骤E；
E、比较时间步长Δt和最小时间步长tmin，如果时间步长Δt大于或等于最小时间步长tmin，则进入步骤D；如果时间步长Δt小于最小时间步长tmin，则进入步骤F；
F、以最小时间步长tmin为步长计算相对位移
Figure PCTCN2016070104-appb-000001
Figure PCTCN2016070104-appb-000002
然后进入步骤J；
G、以碰撞时间间隔Δti为步长来计算相对位移x(n+1)，如果相对位移x(n+1)的绝对值|x(n+1)|小于或等于L-ε，则进入步骤J；如果相对位移x(n+1)的绝对值|x(n+1)|小于或等于L+ε且相对位移的绝对值|x(n+1)|大于L-ε，则进入步骤I；如果相对位移的绝对值|x(n+1)|大于L+ε，则减小碰撞时间间隔Δti，然后进入步骤H；
H、比较碰撞时间间隔Δti和最小时间步长tmin的大小，如果碰撞时间间隔Δti大于或等于最小时间步长，则返回步骤G；如果碰撞时间间隔Δti小于最小时间步长tmin，则取最小时间步长为步长，即t(n+1)＝t(n)+tmin，然后进入步骤I；
I、计算相对位移
Figure PCTCN2016070104-appb-000003
使步长数n增大1，再取碰撞时间常数∈Δt为步长，即t(n+1)＝t(n)+∈Δt，计算相对速度为：
Figure PCTCN2016070104-appb-000004
Figure PCTCN2016070104-appb-000005
计算相对位移为：x(n+1)＝x(n)，然后进入步骤J；
J、判断t(n+1)是否大于M*Tf，如果是，则结束计算，如果否，则步长数n增大1然后返回步骤B；
ε为误差小量，M为周期数，Tf为激励周期，t(n)为第n个时间采样点，x(n)为第n个相对位移，符号上方带原点表示对时间的求导，一个圆点表示一阶导数，两个圆点表示二阶导数，sign为符号函数。
最终得到的M个周期的相对位移x(n+1)、相对速度
Figure PCTCN2016070104-appb-000006
和时间t(n+1)就能良好地反映齿轮的动态传动误差。通过所获得的传动误差，可以对齿轮箱 的精度进行评价。
在对实际齿轮组进行测量时，各项参数值由常规方式测量得到或由用户设定；仿真模拟计算时，参数值由手动设定或者计算机自动设定。参数值包括L、ε、Tf、Rbg、Rbp、θg、θp、外载荷力矩的均值、从动轮转动惯量、为外载荷力矩的波动部分、ω、t0、由于油膜引起的在齿轮副接触正面和侧面的啮合刚度、由于油膜引起的在齿轮副接触正面和侧面的啮合阻尼等。
作为优选，步骤B中，所述碰撞时间间隔计算方式如下：
Figure PCTCN2016070104-appb-000007
式中，
Figure PCTCN2016070104-appb-000008
Figure PCTCN2016070104-appb-000009
Δti取两个值中最小的正实数，Rbg为从动轮的基圆半径，Rbp为主动轮的基圆半径，Tm为外载荷力矩的均值与从动轮转动惯量的比值，Tp为外载荷力矩的波动部分与从动轮转动惯量的比值，ω为外载荷力矩的频率，t0为系统起始时间，k为由于油膜引起的在齿轮副接触正面和侧面的啮合刚度的和与从动轮转动惯量的比值，c为由于油膜引起的在齿轮副接触正面和侧面的啮合阻尼的和与从动轮转动惯量的比值。
作为优选，计算碰撞事件间隔Δti时，如果Δti无正实数，则直接判定Δti＜Δt。
作为优选，步骤D中，缩小时间步长的具体操作为：将Δt减小减少Δt2。
作为优选，步骤D中，减小碰撞时间间隔的具体操作为：将Δti减小减少Δti 2。
本发明带来的实质性效果是，考虑了齿侧间隙对传动误差的影响，在保证一定的计算精度要求之下，尽量减少计算时间，实现计算精度与计算时间的一种更加优化的分配。
附图说明
图1是本发明的一种齿轮啮合区放大图；
图2是本发明的一种参数测量系统示意图；
图3是本发明的一种主动轮和从动轮按标准中心距装配示意图；
图4是本发明的一种主动轮和从动轮粘着连动时的判断情况的示意图；
图5是本发明的一种发生齿与齿之间碰撞时Δt≥Δti的情况的示意图；
图6是本发明的一种利用ε对轮齿之间的碰撞进行再一次精确的情况的示意图；
图7是本发明的一种确定方法流程图；
图中：1电机；2小齿轮；3输入轴；4角度编码仪；5大齿轮；6输出轴；7磁粉制动器。
具体实施方式
下面通过实施例，并结合附图，对本发明的技术方案作进一步具体的说明。
实施例：图2是对实际齿轮组进行测量的系统，包括电机1、主动轮2、输入轴3、角度编码仪4、从动轮5、输出轴6和磁粉制动器7。主动轮2和从动轮5按标准中心距装配如图3所示，图1是齿轮啮合区放大图。
齿轮的传动误差(相对位移)为与传动特性有关的沿啮合线方向度量被动轮上的齿廓在实际啮合时所处位置同理想条件下应处位置之间的差值，可以定义为：
x(t)＝Rbpθp(t)-Rbgθg(t)  (1)
式中，Rbg为从动轮的基圆半径，Rbp为主动轮的基圆半径；θg为从动轮的角位移，θp为主动轮的角位移，角位移也可以通过角度编码仪进行测量。
角度编码仪：测量图2中的输入轴的角速度
Figure PCTCN2016070104-appb-000010
磁粉制动器：给输出轴提供一个阻力矩，该阻力矩为：
Figure PCTCN2016070104-appb-000011
式中，
Figure PCTCN2016070104-appb-000012
和
Figure PCTCN2016070104-appb-000013
分别为外载荷力矩的均值和波动部分的力矩；ω为该阻力矩的频率。
当|x(t)|<L时，根据达朗贝尔原理建立该齿轮系统动力学方程：
Figure PCTCN2016070104-appb-000014
其中，Ig为大齿轮的转动惯量；k1和c1分别为由于油膜引起的在齿轮副接触正面 的啮合刚度和啮合阻尼；k2和c2分别为由于油膜引起的在齿轮副接触侧面的啮合刚度和啮合阻尼；
Figure PCTCN2016070104-appb-000015
和
Figure PCTCN2016070104-appb-000016
分别为外载荷力矩的均值和波动部分的力矩。将公式(3)两边同时除以转动惯量Ig得：
Figure PCTCN2016070104-appb-000017
这里，c＝(c1+c2)/Ig,k＝(k1+k2)/Ig，
Figure PCTCN2016070104-appb-000018
和
Figure PCTCN2016070104-appb-000019
根据公式(4)，可以得到：
Figure PCTCN2016070104-appb-000020
从t0时刻到t时刻，对公式(5)进行积分，可以得到相应的大齿轮的角速度和角位移
Figure PCTCN2016070104-appb-000021
Figure PCTCN2016070104-appb-000022
Figure PCTCN2016070104-appb-000023
将公式(7)与公式(8)带入公式(1)，可得：
Figure PCTCN2016070104-appb-000024
Figure PCTCN2016070104-appb-000025
当|x(t)|＝L时，齿轮系统发生弹性碰撞，此时系统动力学方程为：
Figure PCTCN2016070104-appb-000026
其中，e为碰撞系数，碰撞时间常数∈Δt为非常小的一个正数。在碰撞时，传动误差可以表示为：
Figure PCTCN2016070104-appb-000027
这里，Δti＝ti-t0为碰撞时间间隔；若
Figure PCTCN2016070104-appb-000028
则表示碰撞发生在轮齿的正面，此时，从动轮在主动轮的推动下转动。若
Figure PCTCN2016070104-appb-000029
则表示碰撞发生在轮齿的侧面，此时，主动轮在从动轮的推动下转动。
将公式(9)结合公式(12)得：
Figure PCTCN2016070104-appb-000030
对公式(13)左右两边同时在小量Δti附件进行2阶泰勒公式展开，并忽略比
Figure PCTCN2016070104-appb-000031
高阶的无穷小量，得：
Figure PCTCN2016070104-appb-000032
从上式可以求出：
Figure PCTCN2016070104-appb-000033
式中，
Figure PCTCN2016070104-appb-000034
Figure PCTCN2016070104-appb-000035
这里根据物理意义，Δti只取其中最小的正实数。若Δti无最小正实数，则认为Δt<Δti。
如图7所示，将Δti与Δt进行比较来确定下一个时间步长为多少，如果是Δt<Δti，则以时间步长Δt来计算相对位移x(n+1)。若|x(n+1)|<L，则计算下一时刻的时间Δti。若|x(n+1)|>L，则说明相对位移大于齿侧间隙，此时应该按照公式
Δt＝Δt-Δt2  (16)
缩小时间步长Δt。若公式(16)缩小后的时间步长大于系统所设置的最小时间步长tmin，则重新计算相对位移x(n+1)。反之，系统以最小时间步长tmin为步长计算相对位移x(n+1)。此时系统大小齿轮发生粘动现象(如图4所示)，此时根据其物理意义，可以定义
Figure PCTCN2016070104-appb-000036
如果是Δt≥Δti(如图5所示)，则用Δti作为下一个步长去计算下一次的相对位移与速度。同样，这里我们要对x(n+1)进行判断，如果|x(n+1)|>L+ε，则证明这里的相对位移x(n+1)已经超出了精确范围，如图6所示，这里就需要减小Δti，来保证精确到允许的误差范围内，减小Δti的公式如下：
Figure PCTCN2016070104-appb-000037
以公式(18)中的Δti为步长，重新计算x(n+1)与
Figure PCTCN2016070104-appb-000038
如果此时|x(n+1)|>L+ε，重复公式(18)，直到Δti<tmin发生时，系统将以tmin为步长计算相对位移和相对速度。此时，系统发生碰撞，根据碰撞的实际物理意义，可以令：
Figure PCTCN2016070104-appb-000039
为保障系统在碰撞前后，其相对位移值相等，这里引入一个时间小量，即碰撞时间常数∈Δt，使得系统以∈Δt为步长，即：
t(n+1)＝t(n)+∈Δt  (20)
重新计算一次相对位移和相对速度，其碰撞后的相对速度和相对位移为：
Figure PCTCN2016070104-appb-000040
若不满足公式|x(n+1)|>L+ε，且满足|x(n+1)|>L-ε,系统发生碰撞，此时重复公式(19)，(20)和(21)。
若不满足|x(n+1)|>L+ε,也不满足|x(n+1)|>L-ε，说明系统没有碰撞，此时重新计算位移和速度。
这样不断循环，可以记录下大量的x(n+1)，
Figure PCTCN2016070104-appb-000041
t(n+1)等数据。在其整个运算过程中，若t(n+1)<M*Tf,则把其计算得到的相对位移和相对速度代入公式(15)中重新计算Δti。若t(n+1)>M*Tf则终止运算，即可得到M个周期的相对位移(传动误差)，相对速度(传动误差的速度)和时间的数据。M为激励的周期数，Tf为激励周期，激励由磁粉制动器提供。
综上所述，本发明的齿轮动态传动误差确定方法，通过由角度编码仪4测得的输入轴3角速度、由磁粉制动器7提供的输出轴6阻力矩、齿侧间隙、输入轴3和输出轴6的角位移等参数值计算出碰撞时间间隔，将所得到的相对位移与齿侧间隙进行比较，直至相对位移的值在适当的区间内并且符合要求，所得到的结果即为精度较高的传动误差。本发明引入了齿侧间隙对传动误差的影响，保证了传动误差确定结果的准确性，同时使精度和时间得到了良好的平衡，适用于对实际齿轮组的测量计算和仿真模拟的计算。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
尽管本文较多地使用了相对位移、碰撞时间间隔等术语，但并不排除使用其它术语的可能性。使用这些术语仅仅是为了更方便地描述和解释本发明的本质；把它们解释成任何一种附加的限制都是与本发明精神相违背的。

Claims (5)

  1. 一种齿轮动态传动误差确定方法，其特征在于，包括以下步骤：
    A、获取参数值，所述参数值包括由角度编码仪测得的输入轴角速度和由磁粉制动器提供的输出轴阻力矩；
    B、计算两个轮齿之间的碰撞时间间隔Δti；
    C、比较碰撞时间间隔Δti和时间步长Δt的大小，如果Δt＜Δti，则进入步骤D；如果Δt≥Δti，则进入步骤G；
    D、以时间步长Δt为步长来计算相对位移x(n+1)，如果相对位移x(n+1)的绝对值|x(n+1)|小于或等于齿侧间隙L，则进入步骤J；如果相对位移x(n+1)的绝对值|x(n+1)|大于齿侧间隙L，则缩小时间步长Δt然后进入步骤E；
    E、比较时间步长Δt和最小时间步长tmin，如果时间步长Δt大于或等于最小时间步长tmin，则进入步骤D；如果时间步长Δt小于最小时间步长tmin，则进入步骤F；
    F、以最小时间步长tmin为步长计算相对位移
    Figure PCTCN2016070104-appb-100001
    Figure PCTCN2016070104-appb-100002
    然后进入步骤J；
    G、以碰撞时间间隔Δti为步长来计算相对位移x(n+1)，如果相对位移x(n+1)的绝对值|x(n+1)|小于或等于L-ε，则进入步骤J；如果相对位移x(n+1)的绝对值|x(n+1)|小于或等于L+ε且相对位移的绝对值|x(n+1)|大于L-ε，则进入步骤I；如果相对位移的绝对值|x(n+1)|大于L+ε，则减小碰撞时间间隔Δti，然后进入步骤H；
    H、比较碰撞时间间隔Δti和最小时间步长tmin的大小，如果碰撞时间间隔Δti大于或等于最小时间步长，则返回步骤G；如果碰撞时间间隔Δti小于最小时间步长tmin，则取最小时间步长为步长，即t(n+1)＝t(n)+tmin，然后进入步骤 I；
    I、计算相对位移
    Figure PCTCN2016070104-appb-100003
    使步长数n增大1，再取碰撞时间常数∈Δt为步长，即t(n+1)＝t(n)+∈Δt，计算相对速度为：
    Figure PCTCN2016070104-appb-100004
    Figure PCTCN2016070104-appb-100005
    计算相对位移为：x(n+1)＝x(n)，然后进入步骤J；
    J、判断t(n+1)是否大于M×Tf，如果是，则结束计算，如果否，则步长数n增大1然后返回步骤B；
    ε为误差小量，M为周期数，Tf为激励周期，t(n)为第n个时间采样点，x(n)为第n个相对位移。
  2. 根据权利要求1所述的齿轮动态传动误差确定方法，其特征在于，步骤B中，所述碰撞时间间隔计算方式如下：
    Figure PCTCN2016070104-appb-100006
    式中，
    Figure PCTCN2016070104-appb-100007
    Figure PCTCN2016070104-appb-100008
    Δti取两个值中最小的正实数，Rbg为从动轮的基圆半径，Rbp为主动轮的基圆半径，Tm为外载荷力矩的均值与从动轮转动惯量的比值，Tp为外载荷力矩的波动部分与从动轮转动惯量的比值，ω为外载荷力矩的频率，t0为系统初始时间，k为由于油膜引起的在齿轮副接触正面和侧面的啮合刚度的和与从动轮转动惯量的比值，c为由于油膜引起的在齿轮副接触正面和侧面的啮合阻尼的和与从动轮转动惯量的比值。
  3. 根据权利要求2所述的齿轮动态传动误差确定方法，其特征在于，计算碰撞事件间隔Δti时，如果Δti无正实数，则直接判定Δti＜Δt。
  4. 根据权利要求1或2或3所述的齿轮动态传动误差确定方法，其特征在于，步骤D中，缩小时间步长的具体操作为：将Δt减小减少Δt2。
  5. 根据权利要求1或2或3所述的齿轮动态传动误差确定方法，其特征在 于，步骤D中，减小碰撞时间间隔的具体操作为：将Δti减少
    Figure PCTCN2016070104-appb-100009
PCT/CN2016/070104 2015-01-05 2016-01-05 齿轮动态传动误差确定方法 WO2016110236A1 (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US15/319,228 US10274399B2 (en) 2015-01-05 2016-01-05 Gear dynamic transmission error determination

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201510003696.8A CN104502096B (zh) 2015-01-05 2015-01-05 齿轮动态传动误差计算方法
CN201510003696.8 2015-01-05

Publications (1)

Publication Number Publication Date
WO2016110236A1 true WO2016110236A1 (zh) 2016-07-14

Family

ID=52943516

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2016/070104 WO2016110236A1 (zh) 2015-01-05 2016-01-05 齿轮动态传动误差确定方法

Country Status (3)

Country Link
US (1) US10274399B2 (zh)
CN (1) CN104502096B (zh)
WO (1) WO2016110236A1 (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112035966A (zh) * 2020-07-13 2020-12-04 西安交通大学 一种基于齿轮内源激励力的齿轮振源识别方法
CN113221289A (zh) * 2021-06-01 2021-08-06 江苏科技大学 一种基于区间模型的谐波减速器动态传动误差近似方法
CN113340592A (zh) * 2021-06-03 2021-09-03 河北工业大学 考虑外部振动激励的渐开线齿轮磨损与动力学综合实验台
CN114970100A (zh) * 2022-04-22 2022-08-30 西安交通大学 基于残余传递误差的弧齿锥齿轮故障内源激励识别方法
CN115081132A (zh) * 2022-06-15 2022-09-20 北京航空航天大学 一种面向齿轮装置回差的多层次影响因素耦合分析方法
CN115098988A (zh) * 2022-04-22 2022-09-23 西安交通大学 一种基于传递误差的弧齿锥齿轮内源激励识别方法

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104502096B (zh) 2015-01-05 2016-10-19 盐城工学院 齿轮动态传动误差计算方法
CN105486219A (zh) * 2015-11-06 2016-04-13 中南大学 一种航空高速齿轮动态传动误差测量系统
CN105527100A (zh) * 2016-02-01 2016-04-27 浙江德孚力汽车变速箱有限公司 用于圆锥齿轮传动的啮合检测机构
CN106840652B (zh) * 2017-01-18 2019-03-22 北京工业大学 一种基于齿轮整体误差的动态传动误差预测方法
CN107036810A (zh) * 2017-05-04 2017-08-11 盐城工学院 一种齿轮传动误差的测量方法
CN108470102B (zh) * 2018-03-21 2021-10-22 西安科技大学 一种面向啮合性能预控的小轮齿面设计方法及加工方法
CN108444700B (zh) * 2018-07-02 2021-10-26 南昌航空大学 圆柱齿轮啮合刚度静态测量装置及测试方法
CN109342051B (zh) * 2018-10-23 2020-08-25 重庆理工大学 一种传动间隙测量方法
CN110321655B (zh) * 2019-07-11 2024-01-23 西北工业大学 齿面动态载荷分布计算方法
CN110822073B (zh) * 2019-11-20 2023-02-10 上海新时达机器人有限公司 齿隙调节方法及调节装置
CN111859576B (zh) * 2020-07-27 2024-02-02 大连交通大学 机器人用rv减速器含间隙机构传动误差计算方法
CN112464486B (zh) * 2020-12-02 2024-03-19 华北水利水电大学 一种机器人用摆线针轮减速机的传动误差数值仿真方法
CN113405444B (zh) * 2021-06-15 2023-03-14 中国航发沈阳发动机研究所 一种测量齿轮啮合区油膜厚度的装置及测量方法
DE102022002437A1 (de) 2021-08-03 2023-02-09 Sew-Eurodrive Gmbh & Co Kg Prüfvorrichtung und Verfahren zur Prüfung eines Getriebes mit einer Prüfvorrichtung
CN113607409B (zh) * 2021-10-08 2021-12-17 山东世纪智慧农业科技有限公司 一种用于齿轮的测试系统和方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0933398A (ja) * 1995-07-17 1997-02-07 Nissan Motor Co Ltd 傘歯車の歯当たり伝達誤差シミュレーション方法
CN101246083A (zh) * 2008-03-24 2008-08-20 西安电子科技大学 直齿圆柱齿轮动态啮合刚度的测量方法
CN101271038A (zh) * 2008-05-09 2008-09-24 南京航空航天大学 齿轮副传动精度测量装置和方法
CN102519723A (zh) * 2012-01-09 2012-06-27 唐进元 一种齿轮动态传动误差测控系统
CN103698124A (zh) * 2013-12-12 2014-04-02 盐城工学院 一种测量齿轮传动误差的方法
CN104502096A (zh) * 2015-01-05 2015-04-08 盐城工学院 齿轮动态传动误差计算方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103969044B (zh) * 2014-04-30 2016-05-11 东北大学 一种齿轮传动精度测试实验台
CN104074949B (zh) * 2014-06-18 2016-08-24 盐城工学院 弹性隔层式减振齿轮和降低齿轮传动误差的方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0933398A (ja) * 1995-07-17 1997-02-07 Nissan Motor Co Ltd 傘歯車の歯当たり伝達誤差シミュレーション方法
CN101246083A (zh) * 2008-03-24 2008-08-20 西安电子科技大学 直齿圆柱齿轮动态啮合刚度的测量方法
CN101271038A (zh) * 2008-05-09 2008-09-24 南京航空航天大学 齿轮副传动精度测量装置和方法
CN102519723A (zh) * 2012-01-09 2012-06-27 唐进元 一种齿轮动态传动误差测控系统
CN103698124A (zh) * 2013-12-12 2014-04-02 盐城工学院 一种测量齿轮传动误差的方法
CN104502096A (zh) * 2015-01-05 2015-04-08 盐城工学院 齿轮动态传动误差计算方法

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112035966A (zh) * 2020-07-13 2020-12-04 西安交通大学 一种基于齿轮内源激励力的齿轮振源识别方法
CN112035966B (zh) * 2020-07-13 2024-04-19 西安交通大学 一种基于齿轮内源激励力的齿轮振源识别方法
CN113221289A (zh) * 2021-06-01 2021-08-06 江苏科技大学 一种基于区间模型的谐波减速器动态传动误差近似方法
CN113340592A (zh) * 2021-06-03 2021-09-03 河北工业大学 考虑外部振动激励的渐开线齿轮磨损与动力学综合实验台
CN114970100A (zh) * 2022-04-22 2022-08-30 西安交通大学 基于残余传递误差的弧齿锥齿轮故障内源激励识别方法
CN115098988A (zh) * 2022-04-22 2022-09-23 西安交通大学 一种基于传递误差的弧齿锥齿轮内源激励识别方法
CN114970100B (zh) * 2022-04-22 2024-04-05 西安交通大学 基于残余传递误差的弧齿锥齿轮故障内源激励识别方法
CN115098988B (zh) * 2022-04-22 2024-05-10 西安交通大学 一种基于传递误差的弧齿锥齿轮内源激励识别方法
CN115081132A (zh) * 2022-06-15 2022-09-20 北京航空航天大学 一种面向齿轮装置回差的多层次影响因素耦合分析方法

Also Published As

Publication number Publication date
CN104502096B (zh) 2016-10-19
CN104502096A (zh) 2015-04-08
US20170138819A1 (en) 2017-05-18
US10274399B2 (en) 2019-04-30

Similar Documents

Publication Publication Date Title
WO2016110236A1 (zh) 齿轮动态传动误差确定方法
JP6009397B2 (ja) 位置制御装置のイナーシャ推定方法およびイナーシャ推定装置
CN107545127B (zh) 一种考虑接触的工业机器人关节刚度建模方法
JP6262081B2 (ja) 変速機クラッチトルク推定方法
KR102585751B1 (ko) Sbw시스템의 랙포스 추정방법
CN103698124B (zh) 一种测量齿轮传动误差的方法
CN103195869B (zh) 渐开线圆柱齿轮系统传动精度可靠性确定与调控方法
JP4787172B2 (ja) ウェブ張力のリアルタイム測定および位置センサを用いた制御
CN108120365B (zh) 多齿分度台分度误差的逐齿检测方法
CN112257192B (zh) 一种基于实测齿廓偏差的一齿切向综合偏差计算方法
CN109029978A (zh) 精密减速器几何回差评定方法
CN110398359A (zh) 一种机械传动链综合误差的动态检测方法及装置
CN107314845B (zh) 齿轮动态啮合力的测试方法
Huang et al. Influence of system-inherent phase in spur gear multi-state mesh based on nonlinear dynamics
CN115081132A (zh) 一种面向齿轮装置回差的多层次影响因素耦合分析方法
CN110361967B (zh) 滑模观测器的构建方法
TW201517496A (zh) 馬達參數量測方法與馬達參數量測系統
Ratanasumawong et al. Estimation of sliding loss in a parallel-axis gear pair
JP2018173702A (ja) パラメータ同定装置、駆動システム、パラメータ同定方法及びプログラム
CN108414218B (zh) 一种90°轴交角小模数锥齿轮副双面啮合测量方法
CN116663342A (zh) 一种齿轮动态啮合刚度计算方法及装置
CN110822073B (zh) 齿隙调节方法及调节装置
JP4038052B2 (ja) 歯車設計支援装置、歯車設計支援方法及び記録媒体
CN112035966B (zh) 一种基于齿轮内源激励力的齿轮振源识别方法
JP4658531B2 (ja) 工作機械の熱変位推定方法

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 16734896

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 15319228

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 16734896

Country of ref document: EP

Kind code of ref document: A1