Method of analyzing a vibratory signal derived from rotation of at least one moving part belonging to a rotary mechanism

Patent metadata

Field	Value
Publication number	US-11480460-B2
Application number	US-201916508696-A
Country	US
Kind code	B2
Filing date	Jul 11, 2019
Priority date	Jul 12, 2018
Publication date	Oct 25, 2022
Grant date	Oct 25, 2022

Abstract

A method of analyzing a vibratory signal derived from rotation of at least one moving part belonging to a rotary mechanism forming all or part of a drive train for transmitting drive torque, the rotary mechanism being fitted to an aircraft and the method comprising at least one first measurement step including measuring vibration in at least one direction and generating a vibratory signal representative of the operation of the rotary mechanism as a function of time, the first measurement step being performed by means of at least one vibration sensor; and at least one second measurement step including measuring an angular position of the moving part, the moving part having at least one degree of freedom to move in rotation about a respective axis of rotation Z. Such an analysis method makes it possible to determine at least one usable range for the vibratory signal.

First claim

What is claimed is:

1. A method of analyzing a vibratory signal derived from rotation of at least one moving part belonging to a rotary mechanism forming all or part of a drive train for transmitting drive torque, the rotary mechanism being fitted to an aircraft and the method comprising:
at least one first measurement step including measuring vibration in at least one direction and generating a vibratory signal representative of the operation of the rotary mechanism as a function of time, the first measurement step being performed by means of at least one vibration sensor; and
at least one second measurement step including measuring an angular position of the moving part(s), the moving part(s) having at least one degree of freedom to move in rotation about a respective axis of rotation Z, the second measurement step(s) serving to count a determined number n of rotations of the moving part(s) about the respective axis of rotation Z;
wherein the method further comprises at least:
a preprocessing step for calculating a plurality of arguments of complex numbers generated from complex vibrations of the vibratory signal measured by the vibration sensor(s) over a cycle of the predetermined number n of rotations of the moving part(s) about the respective axis of rotation Z;
a first analysis step for determining an angle offset relating to the plurality of arguments, the first analysis step serving to generate a first analysis result D 1 , the first analysis result D 1 being a condition that is satisfied or not satisfied as a function of the angle offset;
a second analysis step for determining a first distortion, referred to as “low” distortion, of the plurality of arguments, the second analysis step serving to generate a second analysis result D 2 , the second analysis result D 2 being a condition that is satisfied or not satisfied as a function of the first distortion;
a third analysis step for determining local instability of the plurality of arguments, the third analysis step serving to generate a third analysis result D 3 , the third analysis result D 3 being a condition that is satisfied or not satisfied as a function of the local instability;
a fourth analysis step for determining a second distortion, referred to as a “high” distortion, of the plurality of arguments, the fourth analysis step serving to generate a fourth analysis result D 4 , the fourth analysis result D 4 being a condition that is satisfied or not satisfied as a function of the second distortion, the second distortion being distinct from the first distortion;
a validation step for determining at least one usable time range for the vibratory signal, the validation step depending simultaneously on the first analysis result D 1 , on the second analysis result D 2 , on the third analysis result D 3 , and on the fourth analysis result D 4 ;
the validation step being successful in determining the at least one usable time range for the vibratory signal as each of the first analysis result D 1 , the second analysis result D 2 , the third analysis result D 3 , and the fourth analysis result D 4 simultaneously is a respective condition that is not satisfied during part of the time that the vibratory signal is generated;
a storing step for storing the vibratory signal only in the at least one usable time range in an onboard memory of the aircraft;
a monitoring step for monitoring wear of the at least one moving part based on the vibratory signal over the at least one usable time range stored in the onboard memory; and
wherein the preprocessing step, the first analysis step, the second analysis step, the third analysis step, the fourth analysis step, the validation step, and the storing step are performed on board the aircraft while the aircraft is in flight whereby the at least one usable time range for the vibratory signal is determined directly while the aircraft is in flight with a quantity of data that is stored on the onboard memory pertaining to the vibratory signal being limited to only the vibratory signal in the at least one usable time range.

2. The method according to claim 1 , wherein the validation step is further performed on the ground, the aircraft including an onboard memory for continuously storing the vibratory signal.

3. The method according to claim 1 , wherein the method includes a data transmission step enabling data representative of the vibratory signal to be transmitted to at least one ground station.

4. The method according to claim 3 , wherein the data transmission step takes place while being simultaneously dependent on the first analysis result D 1 , on the second analysis result D 2 , on the third analysis result D 3 , and on the fourth analysis result D 4 .

5. The method according to claim 1 , wherein the preprocessing step comprises:
a calculation substep for calculating a first moving window Fourier transform from the vibratory signal; and
a second calculation substep for calculating a first matrix ANG 1→n of n arguments ANG k′ , where k varies over the range 1 to n, where n corresponds to a number of rotations of the moving part(s) about the respective axis of rotation Z.

6. The method according to claim 5 , wherein the first analysis step comprises:
a calculation substep for calculating a second matrix of angle cosines P 1→n from n cosine values P k of the first matrix ANG 1→n , and such that:
P k =cos(ANG k )

a calculation substep for calculating a mean value P moy from the n cosine values P k constituting the second matrix P 1→n ;
an identification substep for identifying a median signal P med by calculating the minimum Euclidean distance between the mean value P moy and each of the n cosine values P k constituting the second matrix P 1→n ;
a calculation substep for calculating n inter-correlation values {circumflex over (R)} k , such that dim({circumflex over (R)} k )=2Z c −1, the n inter-correlation values {circumflex over (R)} k being calculated by taking the convolution product between the median signal P med and each of the n cosine values P k constituting the second matrix P 1→n , where Z c corresponds to a number of teeth of a stationary ring co-operating with the moving part(s);
a calculation substep for calculating n time offset values T k′ , such that T k corresponds to an abscissa value for an absolute maximum of a curve representative of the n cross-correlation values {circumflex over (R)} k , k varying over the range 1 to n, and when k=P med′ T P med =argmax({circumflex over (R)} P med )=0; and
a first diagnosis substep for generating the first analysis result D 1 that is a condition that is satisfied if at least one of the n values of the time offset T k , where k varies over the range 1 to n, is greater than a first predetermined threshold value S 1 , and conversely in a condition that is not satisfied if all of the n values of the time offset T k are less than or equal to the first predetermined threshold value S 1 .

7. The method according to claim 6 , wherein the second analysis step comprises:
a resetting substep of resetting the n cosine values P k constituting the second matrix P 1→n by a number corresponding to the n values of a time offset T k in order to generate a third matrix C 1→n of n reset values C k for k varying over the range 1 to n;
a calculation substep for calculating a second moving window Fourier transform from the median signal P med ;
a calculation substep for calculating a mean argument value φ moy for a predetermined harmonic H m of the second Fourier transform, the mean argument value φ moy being such that:
φ moy =atan2(y,x)=Arg(x+iy)

where x=Re(Hm) is the real part of the predetermined harmonic H m of the second Fourier transform and
y=Im(Hm) is the imaginary part of the

Assignees

• Airbus Helicopters

Inventors

• Morel Herve
• Rodriguez Jonathan

Classifications

• G01H1/04

  of vibrations which are transverse to direction of propagation · CPC title

• G01H1/003 Primary

  of rotating machines (G01H1/10 takes precedence) · CPC title

• G01M13/028

  Acoustic or vibration analysis · CPC title

• B64C27/001

  Vibration damping devices · CPC title

• B64C27/12

  Rotor drives · CPC title