Gear Transmission Error Measurements based on the Phase Demodulation

Similar documents
PHASE DEMODULATION OF IMPULSE SIGNALS IN MACHINE SHAFT ANGULAR VIBRATION MEASUREMENTS

SIMPLE GEAR SET DYNAMIC TRANSMISSION ERROR MEASUREMENTS

Lecture on Angular Vibration Measurements Based on Phase Demodulation

NOISE AND VIBRATION DIAGNOSTICS IN ROTATING MACHINERY

Theory and praxis of synchronised averaging in the time domain

VOLD-KALMAN ORDER TRACKING FILTERING IN ROTATING MACHINERY

Compensating for speed variation by order tracking with and without a tacho signal

ANGULAR VIBRATION MEASUREMENTS OF THE POWER DRIWING SYSTEMS

Instrumentation and signal processing methods used for machine diagnostics

NOISE REDUCTION IN SCREW COMPRESSORS BY THE CONTROL OF ROTOR TRANSMISSION ERROR

SEPARATING GEAR AND BEARING SIGNALS FOR BEARING FAULT DETECTION. Wenyi Wang

Rotating Machinery Fault Diagnosis Techniques Envelope and Cepstrum Analyses

Fault Diagnosis of Wind Turbine Gearboxes Using Enhanced Tacholess Order Tracking

Real-time Math Function of DL850 ScopeCorder

Circuit Analysis-II. Circuit Analysis-II Lecture # 2 Wednesday 28 th Mar, 18

APPLICATION NOTE 3560/7702. Introduction

9LEUDWLRQ 0HDVXUHPHQW DQG $QDO\VLV

THE SINUSOIDAL WAVEFORM

Application Note. Monitoring strategy Diagnosing gearbox damage

A study of Vibration Analysis for Gearbox Casing Using Finite Element Analysis

Fundamentals of Vibration Measurement and Analysis Explained

Current-Based Diagnosis for Gear Tooth Breaks in Wind Turbine Gearboxes

NON-SELLABLE PRODUCT DATA. Order Analysis Type 7702 for PULSE, the Multi-analyzer System. Uses and Features

GEAR TRANSMISSION ERROR MEASUREMENT ACCURACY USING LOW-COST DIGITAL ENCODERS

Copyright 2017 by Turbomachinery Laboratory, Texas A&M Engineering Experiment Station

Prognostic Health Monitoring for Wind Turbines

MODEL MODIFICATION OF WIRA CENTER MEMBER BAR

Vibration Analysis on Rotating Shaft using MATLAB

Phasor. Phasor Diagram of a Sinusoidal Waveform

Detection of Wind Turbine Gear Tooth Defects Using Sideband Energy Ratio

Congress on Technical Diagnostics 1996

A Multi-Probe Setup for the Measurement of Angular Vibrations in a Rotating Shaft

IET (2014) IET.,

Cepstral Removal of Periodic Spectral Components from Time Signals

CHAPTER 6 INTRODUCTION TO SYSTEM IDENTIFICATION

CONTACTLESS MEASURING METHOD OF BLADE VIBRATION DURING TURBINE SPEED-UP

Condition based monitoring: an overview

Shaft encoders are digital transducers that are used for measuring angular displacements and angular velocities.

BLADE AND SHAFT CRACK DETECTION USING TORSIONAL VIBRATION MEASUREMENTS PART 2: RESAMPLING TO IMPROVE EFFECTIVE DYNAMIC RANGE

Gear Noise Prediction in Automotive Transmissions

Chapter 4 REVIEW OF VIBRATION ANALYSIS TECHNIQUES

Exercise 1. Consider the following figure. The shaded portion of the circle is called the sector of the circle corresponding to the angle θ.

Signal Analysis Techniques to Identify Axle Bearing Defects

A Comparison of Different Techniques for Induction Motor Rotor Fault Diagnosis

Active Vibration Isolation of an Unbalanced Machine Tool Spindle

Impact of Indexing Errors on Spur Gear Dynamics

Data Sheet. AEDT-9140 Series High Temperature 115 C Three Channel Optical Incremental Encoder Modules 100 CPR to 1000 CPR. Description.

ME scope Application Note 01 The FFT, Leakage, and Windowing

PeakVue Analysis for Antifriction Bearing Fault Detection

Extraction of tacho information from a vibration signal for improved synchronous averaging

Fault Diagnosis of Gearbox Using Various Condition Monitoring Indicators for Non-Stationary Speed Conditions: A Comparative Analysis

Linear Time-Invariant Systems

Research Article Vibration Sideband Modulations and Harmonics Separation of a Planetary Helicopter Gearbox with Two Different Configurations

Machine Diagnostics in Observer 9 Private Rules

Investigating the Electromechanical Coupling in Piezoelectric Actuator Drive Motor Under Heavy Load

Fault Detection of Double Stage Helical Gearbox using Vibration Analysis Techniques

Precision power measurements for megawatt heating controls

Fault diagnosis of massey ferguson gearbox using power spectral density

Practical performances of high-speed measurement of gear Transmission Error or torsional vibrations with optical encoders

Review on Fault Identification and Diagnosis of Gear Pair by Experimental Vibration Analysis

APPLICATION NOTE. Detecting Faulty Rolling Element Bearings. Faulty rolling-element bearings can be detected before breakdown.

Effect of parameters setting on performance of discrete component removal (DCR) methods for bearing faults detection

Sinusoids and Phasors (Chapter 9 - Lecture #1) Dr. Shahrel A. Suandi Room 2.20, PPKEE

ACOUSTIC NOISE AND VIBRATIONS OF ELECTRIC POWERTRAINS

Synchronization of Matsuoka Neural Oscillator to Meshing Frequency for Estimation of Angular Position of Gears

Optical Encoder Applications for Vibration Analysis

Fault diagnosis of Spur gear using vibration analysis. Ebrahim Ebrahimi

Combining synchronous averaging with a Gaussian mixture model novelty detection scheme for vibration-based condition monitoring of a gearbox

Theory of Telecommunications Networks

Study of Improper Chamfering and Pitting Defects of Spur Gear Faults Using Frequency Domain Technique

Fault detection of a spur gear using vibration signal with multivariable statistical parameters

Pulse Code Modulation (PCM)

Chapter 6: Alternating Current. An alternating current is an current that reverses its direction at regular intervals.

Gearbox Vibration Source Separation by Integration of Time Synchronous Averaged Signals

Appearance of wear particles. Time. Figure 1 Lead times to failure offered by various conventional CM techniques.

Frequency Demodulation Analysis of Mine Reducer Vibration Signal

Bakiss Hiyana binti Abu Bakar JKE, POLISAS BHAB

As before, the speed resolution is given by the change in speed corresponding to a unity change in the count. Hence, for the pulse-counting method

Response spectrum Time history Power Spectral Density, PSD

Novel Spectral Kurtosis Technology for Adaptive Vibration Condition Monitoring of Multi Stage Gearboxes

Applying digital signal processing techniques to improve the signal to noise ratio in vibrational signals

Chapter 3 Simulation studies

Influence of the gear geometry and the machine on the power-skiving cutter design

UNIT Explain the radiation from two-wire. Ans: Radiation from Two wire

6.1 - Introduction to Periodic Functions

Frequency Domain Analysis

CONTROLLING THE OSCILLATIONS OF A SWINGING BELL BY USING THE DRIVING INDUCTION MOTOR AS A SENSOR

Practical Machinery Vibration Analysis and Predictive Maintenance

Investigating the Electromechanical Coupling in Piezoelectric Actuator Drive Motor Under Heavy Load

Multiparameter vibration analysis of various defective stages of mechanical components

CHAPTER 14 ALTERNATING VOLTAGES AND CURRENTS

Topic 2. Signal Processing Review. (Some slides are adapted from Bryan Pardo s course slides on Machine Perception of Music)

Structure of Speech. Physical acoustics Time-domain representation Frequency domain representation Sound shaping

DIAGNOSIS OF ROLLING ELEMENT BEARING FAULT IN BEARING-GEARBOX UNION SYSTEM USING WAVELET PACKET CORRELATION ANALYSIS

END-OF-YEAR EXAMINATIONS ELEC321 Communication Systems (D2) Tuesday, 22 November 2005, 9:20 a.m. Three hours plus 10 minutes reading time.

Enayet B. Halim, Sirish L. Shah and M.A.A. Shoukat Choudhury. Department of Chemical and Materials Engineering University of Alberta

Laboratory Assignment 5 Amplitude Modulation

Degradation of BER by Group Delay in Digital Phase Modulation

Data Acquisition Systems. Signal DAQ System The Answer?

CHAPTER 9. Sinusoidal Steady-State Analysis

Transcription:

Gear Transmission Error Measurements based on the Phase Demodulation JIRI TUMA Abstract. The paper deals with a simple gear set transmission error (TE) measurements at gearbox operational conditions that means under load and during rotation. The analysis method is focused on the processing of pulse signals generated by encoders attached to both the gears in mesh. The analysis technique benefits from demodulation of a phase-modulated signal. The theory is illustrated by experiments with a car gearbox and measurement errors are discussed. Keywords: Gears, gearboxes, transmission error, encoders, phase demodulation, synchronised averaging, average tooth mesh, measurement accuracy, angular vibration. 1. Introduction Noise and vibration problems in gearing are mainly concerned with the smoothness of the drive. The parameter that is employed to measure smoothness is Transmission Error (T. E.). This parameter can be expressed as a linear displacement at a base circle radius defined by the difference of the output gear s position from where it would be if the gear teeth were perfect and infinitely stiff. Many references have attested to the fact that a major goal in reducing gear noise is to reduce the transmission error of a gear set. The basic equation for T.E. of a simple gear set is given as TE n = 2 r n (1) 1 2 ( m) Θ Θ 1 2 where n 1, n2 - teeth numbers of pinion and wheel respectively, Θ 1, Θ 2 are angles of rotation of the mentioned gears and r 2 is a wheel radius. T.E. results not only from manufacturing inaccuracies such as profile errors, tooth pitch errors and run-out, but from a bad design. The pure tooth involute deflects under load due to the finite mesh stiffness caused by tooth deflection. A gearcase and shaft system deflects due to load as well. To complicate meters, a gear axis system is a very complicate dynamic system. Department of Control Systems and Instrumentation, VŠB - Technical University of Ostrava, 17.listopadu 15, Ostrava-Poruba, 78 33, Czech Republic, tel.: 596993482 jiri.tuma@vsb.cz

While running under load one of very important parameters, tooth contact stiffness, is varying what excites the parametric vibration and consequently noise. There are many possible approaches to measuring T.E., but, as Derek Smith points out [1], in practice, measurements based on the use of encoders dominates. A possible alternative is the use of tangentially mounted accelerometers to measure the torsional accelerations of the two gears, which causes the problem of signal transfer from the rotational parts and its double integration. As it will shown later, the common encoders generating a string of pulses in hundreds per rotation are satisfying the sever requirements for measurement accuracy for the mentioned quantity in microns. The National Engineering Laboratory at East Kilbride developed a first workable system for laboratory use while Dr. R. G. Munro introduced the redesigned system for industrial use in the 196s. The sketch of the gear set consisting of the 21- and 44-tooth gears under test and attached incremental rotary encoders, designated by E1 and E2 is shown in figure 1. Both the encoders are of Heidenhain origin, the ERN 46-5 type. A perfectly uniform rotation of gear produces an encoder signal having in its frequency spectrum a single component at the frequency that is a multiple of the gear rotational frequency. As both the encoders generate 5 impulses per encoder rotation, the frequency of the single components in orders (a multiple of the encoder rotational frequency) is equal to the same number as the number of the impulses. A gearbox is a machine running in cyclic fashion, which is the reason to prefer frequency spectra of signals in term of dimensionless orders rather than frequency in Hz. The order spectra are evaluated using time records that are measured in dimensionless revolutions rather than seconds and the corresponding FFT spectra are measured in dimensionless orders rather than frequency. This technique is called order analysis or tracking analysis, as the rotation frequency is being tracked and used for analysis. The resolution of the order spectrum is equal to the reciprocal value of the revolution number per record corresponding to input data for the Fast Fourier Transform (FFT). As the measurement is focused on one complete gear revolution, the spectrum resolution is one order. Pulse signals from encoders are recorded by PULSE, the n 1 =21T Brüel & Kjær signal analyser. To simplify the phase pinion demodulation an Order Analyser instrument was employed Θ 1 which resulted in time history records corresponding to one E 1 complete gear revolution. A method of synchronised averaging in the time domain was employed for reducing random noise in Θ 2 wheel the measured data. As it is known, the order analysis is based on E 2 data resampling in such a way, that sampling frequency follows the mean frequency of shaft rotational speed during one complete shaft revolution. The mean rotational frequency is n 2 =44T evaluated by means a train of pulses generated once per a shaft revolution. The time interval between two consecutive pulses is determined by interpolation some 5 times more accurate than indicated by the actual sampling interval. Therefore the pulses Fig. 1. Measurement arrangement distribution inside this time record gives information about the instantaneous rotation angle of each of the gears under test. As a consequence of Shannon s sampling theorem a few pulses must be recorded during each mesh cycle. It means, that the number of pulses produced per encoder revolution must be a multiple of the tooth number. If five harmonics of toothmeshing frequency are required then the number of pulses per gear revolution must be at least ten times

higher than the number of teeth. The encoder generating 5 pulses per revolution seems to be an optimum. The length of resampled time record equals to 248 samples per gear revolution. The sample number is a power of two, which is required by FFT and in corresponding order spectrum, ranging to 8 orders, there is a space for ± 3 sideband components around the carrying component of 5 orders in a frequency spectrum. The frequency range of the mentioned Order Analyser in the described conditions limits the gear rotational speed to the value of 19 RPM. Gear loading has not any influence on the discussed sampling problem. As the pulse signal is not a sinusoid, its frequency spectrum consists of several harmonics of the basic pulse frequency. Variation of angular velocity during rotation causes the phase modulation of the carrying signal that can be any of the harmonic components. The phase modulation gives rise to sidebands around the carrying frequency in the frequency spectrum of the modulated harmonic signal. The phase modulation signal can be derived from the phase of the analytical signal that is evaluated using the Hilbert Transform technique, which is the topic of the following paper section. 2. Phase demodulation using the Hilbert transform An analytic time signal is a useful tool for demodulation. The analytic signal, z () t, is generally combined from real and imaginary parts z () t x() t + j y( t) = z( t) ( jϕ( t) ) = exp, (2) where y () t - the Hilbert transform of x ( t). The relationship between the FFT of the y ( t) and x( t) i ( N i) X i of the length, N, is given by Y = jsign 2. (3) As the angle of the complex values ranges from π to +π () t = arctan( y( t) x( t) ) ϕ, (4) the true angle of the analytical signal as the time function with jumps at π or +π must be obtained by unwrapping that is based on the fact that the absolute value of the difference between two consecutive angles is less than π. The principle of phase unwrapping for a harmonic signal modulated by another harmonic signal is shown in figure 2. The relationship between the phase of the analytical ϕ t, is as follows signal and the phase modulation signal, ( ) M () t = ω t + ϕ ( t) ϕ, (5) M where ω - an angular frequency of the carrier component. The phase modulation signal is a fluctuation of the phase angle around the linear term, ω t.

2π +π π 4 2 Unit -2-4,1,2,3,4,5,6,7,8,9 1 Discontinuities removing (2f f samp ϕ π) ϕ < π ϕ + 2π ϕ, ϕ > +π ϕ 2π ϕ 7 6 rad 5 4 3 2 1,2,4,6,8 1,15,1,5 -,5 -,1 -,15 rad,2,4,6,8 1 Fig. 2. Principle of phase unwrapping 3. Encoder accuracy To evaluate errors in pulses distribution against the angle of rotation, both the encoders were mounted on a shaft what ensured the same rotational speed of them. As the running was not perfectly uniform from the point of the measurement method sensitivity, both the pulse signals were under influence of phase modulation. Using the analytical method described above, the difference between modulation signals gives the error in pulse distribution. RMS of Error in radians + 168 RPM x 635 RPM o 139 RPM Circle part E 2 E 1 1/order order Fig. 3. Encoder accuracy

The frequency spectrum of the resulting error is shown in figure 3. The frequency axis is in orders. The quantity order determines a part of a circle related to the error level in the spectrum. The error level at the distance corresponding to the tooth pitch of the adjacent teeth determines the final accuracy of the T.E. measurement. As it is evident the magnitude of an error at 21 and 44 order is less than 1-5 radians, i.e. approximately 2 angular seconds. 4. T.E. measurements A perfectly uniform rotation of gear produces an encoder signal having in its frequency spectrum a single component at the frequency that is a multiple of the gear rotational frequency. As both the encoders generate 5 impulses per a shaft rotation, the frequency of the single components in order is equal to the same number as it is the number of the impulses, i. e. 5 orders. The gear speed variation as an effect of loaded teeth deflection results in the phase modulation of the impulse signal base frequency. The phase-modulated signal contains sideband components around the carrying component. The distance of the dominating sideband components from the carrying components equals to the integer multiple of the tooth number as it is shown in figures 4 and 5. The frequency axis of both the frequency spectra is in order; it means the multiples of the gear rotational frequency. The frequency of the carrying component is equal to 5 orders while the sideband component associated with the corresponding gear is at the distance of ± 21k or ± 44k (where k = 1, 2,... ) order units from the mentioned carrying component frequency. Take notice of the fact that the dominating components in both the sidebands exceed the background noise level at least 1 times or even more. Both the spectra were evaluated from time signals that are a result of synchronised averaging of 1 revolutions of gears under test. Phase demodulation, which is based V 1 1,1,1 1 1 V,1,1 1 1 Enhanced Spectrum - 21-tooth gear 395 416 437 458 479 5 521 542 563 584 65 Orders Fig. 4. Frequency spectrum of phase modulated signal generated by the E1 encoder 1 1 Enhanced Spectrum - 44-tooth gear 28 324 368 412 456 5 544 588 632 676 72 Orders on the theory of the Hilbert transform and phase unwrapping, enable the evaluation of the angular vibration of both the gears in Fig. 5. Frequency spectrum of phase modulated signal generated by the E2 encoder

mesh individually. The unwrapped phase of the frequency modulated signal that is produced by the encoder E1 is shown in figure 6. The diagram forms almost a straight-line function. After subtracting a linear term from unwrapped phase, the dependence of the phase variation on time is obtained as it is shown in figure 7. Take note of the diagram scale. The relationship between the unwrapped phase and phase modulation signal results from the formulae (5). Unwrap Angular Vibration in 36 27 18 9,,3,5,8 1, Time t /T Angular Vibration in,3,25,2,15,1,5 -,5 -,1 -,15 -,2,,2,4,6,8 1, Fig. 6. Unwrapped phase of E1 encoder signal versus time Fig. 7. Phase modulation signal versus time As it is shown in figure 8 the frequency spectrum of phase modulation signal contains a family of harmonics components to the basic 21-order component excited by toothmeshing of the 21-tooth gear. The low frequency components in the frequency spectrum result from nonuniform driving torque and overlap the toothmeshing response. The phase modulation signal in rees during the pinion revolution is shown in figure 9. The enhanced signal contains five harmonics of the toothmeshing frequency, each of them with 3 pairs of sidebands that cause the amplitude modulation of angle variation. When all these sidebands are removed a purely periodic signal is obtained. The filtration in the frequency domain can be considered as an averaging of the second stage. Therefore, one of these periods corresponding to the gear tooth pitch rotation can be taken as a representative to characterize angular vibration in average. The result of mentioned averaging is called the average toothsmesh. The term averaged toothmesh was introduced to associate vibration and noise measurement with a gear design [2]. 1,1,1 1 1 Spectrum of Ang. Vibration 21 42 63 84 15 Orders Fig. 8. Frequency spectrum of phase modulation signal

Toothmeshing frequency harmonics with 3 sideband components Enhanced Time Signal,2,1 -,1 -,2,,2,4,6,8 1, One out of 21 periods Enhanced Time Signal Average without sidebands Av erage Toothmesh,2,2,1,1 de g -,1 -,1 -,2 -,2,5 1,,2,4,6,8 1, Tooth Pitch Rotation Gear revolution Fig. 9. The second stage of angular vibration averaging for the 21-tooth gear The same average toothmesh in angular variation can be evaluated for the 44-tooth gear as it is shown in figure 1. Angle variation can be easily transformed into the arc length variation. The difference between both the angle variations gives the transmission error. The only problem consists in the true phase delay between these periodic signals because the signals from the encoders are recorded separately. Solving of this problem is based on the similarity of responses both the gears to dynamic forces acting between mating teeth, for instance in acceleration some point on the gearcase. Both the encoder pulse signals are sampled together with the acceleration signal. Two-stage averaging of the twice-measured acceleration signal gives average toothmesh responses that are delayed. The lag for the maximum correlation gives the relative delay.,2,1 -,1 Enhanced Time Signal -,2,,2,4,6,8 1, Average Toothmesh,1,5 de g -,5 -,1,5 1,5 Tooth Pitch Rotation Gear revolution Fig. 1 The second stage of angular vibration averaging for the 44-tooth gear

T.E. is given as the difference between the angular vibration signals in the arc length produced by the mating gears. The result is shown in Figure 11. All the experimental data was taken from a car gearbox. The results correspond to the rotational speed of 5 RPM at the input shaft and almost full load. The measurement method was tested at the maximum rotational speed of 125 RPM. µm 5 4 3 2 1-1 - - 2 - - 3 - - 4 - - 5,2,4,6,8 1 Tooth pitch rotation 21 21 -- tooth tooth gear gear 44-tooth 44-tooth gear gear T.E. T.E. Fig. 11. Transmission error against rotation angle in range of the tooth pitch Conclusion The paper is focused on the problem of the simple gear set transmission error measurement (T.E.). Variation of T.E. is the cause of angular vibration of both the mating gears and consequently gearcase vibration and noise. This paper deals with only one measurement method that is based on the use of encoders generating a string of 5 pulses per gear revolution and reviews the phase demodulation for evaluating of the gear angular vibration. The theory is illustrated by experimental data. Acknowledgements This research has been conducted at the Technical University of Ostrava as a part of the research project No. CEZ 3212 and has been supported by the Czech Ministry of Education. The author benefits from the research work done for the SKODA Auto Company. References [1] Derek Smith J. Gear Noise and Vibration, 1 st ed. New York Basel : Marcel Dekker Inc., 1999. ISBN: -8247-65- [2] Tůma, Jiří. Analysis of Periodic and Quasi-Periodic Signals in Time Domain, In: Proceedings of the Noise 93, St. Petersburg (Russia). Auburn (USA) : Auburn University, 1993, Volume 6, pp. 245-25.