EP3732457B1 - Procédé et dispositif de surveillance d'un roulement équipant un dispositif tournant - Google Patents

Procédé et dispositif de surveillance d'un roulement équipant un dispositif tournant Download PDF

Info

Publication number
EP3732457B1
EP3732457B1 EP18842779.3A EP18842779A EP3732457B1 EP 3732457 B1 EP3732457 B1 EP 3732457B1 EP 18842779 A EP18842779 A EP 18842779A EP 3732457 B1 EP3732457 B1 EP 3732457B1
Authority
EP
European Patent Office
Prior art keywords
cyclic
bearing
frequency
defect
vibration signal
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.)
Active
Application number
EP18842779.3A
Other languages
German (de)
English (en)
Other versions
EP3732457A1 (fr
Inventor
Dany ABBOUD
Mohamed EL BADAOUI
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.)
Safran SA
Original Assignee
Safran SA
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 Safran SA filed Critical Safran SA
Publication of EP3732457A1 publication Critical patent/EP3732457A1/fr
Application granted granted Critical
Publication of EP3732457B1 publication Critical patent/EP3732457B1/fr
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/04Bearings
    • G01M13/045Acoustic or vibration analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H1/00Measuring characteristics of vibrations in solids by using direct conduction to the detector
    • G01H1/003Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M15/00Testing of engines
    • G01M15/04Testing internal-combustion engines
    • G01M15/12Testing internal-combustion engines by monitoring vibrations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M15/00Testing of engines
    • G01M15/14Testing gas-turbine engines or jet-propulsion engines

Definitions

  • the invention relates to the general field of machines or rotating devices, and relates more particularly to bearings fitted to such machines or devices, such as, for example, ball or roller bearings.
  • Bearings and in particular ball or roller bearings, are rotating mechanical devices widely used in rotating machines and whose role is to drive transmission shafts into rotation.
  • Such a bearing is generally composed of two coaxial rings (an inner ring and an outer ring) between which are placed rolling elements (balls or rollers), held in a cage.
  • the rotary movement of the shaft is provided by the movement of the rolling elements.
  • bearings Due to their elementary role in mechanical transmission chains, bearings are permanently stressed and therefore weakened. They may suffer from deterioration in the surface condition or in the shape of the inner or outer rings, as well as of the rolling elements. Damage to a bearing can cause the machine in which it is used to stop unexpectedly.
  • the principle on which the vibration analysis is based is as follows: when a rolling element (ball or roller) of the bearing comes into contact with a defective surface of the bearing (for example due to a spalling or a crack on the surface of one of the bearing rings), periodic shocks occur which produce a vibration at a specific frequency, characteristic of the defect affecting the bearing, and which change the statistical structure of the vibratory signal acquired by the accelerometer.
  • a theoretical value of the specific frequency of the fault can be easily obtained from the geometric characteristics of the bearing and its kinematics.
  • the defect affecting the bearing is then detected by identifying, in the vibratory signal acquired by the accelerometer, a component at the specific frequency of this defect, for example by carrying out a spectral analysis or an analysis of the envelope of the vibratory signal.
  • the document EP 2,693,176 describes a method for detecting defects in a bearing by vibratory analysis attempting to solve the aforementioned problems.
  • This method is based on a pre-processing by means of an autoregressive filter of a vibratory signal acquired by means of an accelerometer placed on the housing of the system comprising the bearings to be monitored, followed by an analysis of the signal envelope. vibration resulting from the pretreatment.
  • this detection method presents a relatively high complexity in terms of implementation, in particular because of the preprocessing carried out.
  • the efficiency of the pre-treatment carried out largely depends on the determination of several parameters (eg order of the autoregressive filter, etc.) which can prove to be complicated in practice.
  • Another example of a method is described by P.
  • the invention can be applied to various types of bearing, and in particular to a ball bearing or to a roller bearing.
  • the invention therefore proposes a technique for monitoring a bearing fitted to a rotating device based on a cyclostationary analysis of the vibratory signal acquired by an accelerometric sensor.
  • the accelerometric sensor is placed in such a way that the vibratory signal that it acquires contains a vibratory signature of the bearing.
  • the accelerometric sensor can be placed near the bearing or on the housing of the system which includes the bearing to be monitored, but it can also be placed at a point remote from the bearing. provided that a sufficient signal-to-noise ratio is ensured (typically of the order of 5%).
  • the cyclostationary analysis proposed by the invention makes it possible to extract from the vibratory signal a diagnostic indicator of a fault likely to affect the bearing (eg fault affecting an internal ring of the bearing, fault affecting an outer ring of the bearing, etc. .) which, when compared to a threshold determined for this fault, makes it possible to detect whether the bearing is affected or not by said fault.
  • This diagnostic indicator is advantageously calculated after having applied a (pre-) treatment to the vibratory signal consisting in ridding it of its deterministic component; this deterministic component is typically generated by rotating elements other than the monitored bearing (which itself contributes to the random component of the vibratory signal), and capable of interfering with the vibratory signal acquired by the sensor.
  • the derived indicator in accordance with the invention is determined for a given fault liable to affect the bearing since it depends on an actual fault frequency estimated from the characteristic theoretical frequency of the fault and from a deviation maximum determined around this theoretical frequency, each fault resulting in a different theoretical frequency and a deviation maximum possible of the real frequency of the fault compared to this theoretical frequency.
  • an indicator associated with a fault on the internal ring of the bearing an indicator associated with a fault on the outer ring of the bearing, an associated indicator to a defect on the bearing cage, an indicator associated with a defect on a rolling element of the bearing, etc.
  • the cyclic coherence of the processed vibratory signal makes it possible to bring out a signature of the defect, including when it is "drowned” in a significant background noise (linked for example to a plurality of defective rotating elements located nearby. of the monitored bearing).
  • a significant background noise linked for example to a plurality of defective rotating elements located nearby. of the monitored bearing.
  • the inventors had the judicious idea of exploiting the cyclostationary characteristic of the vibration signals generated by the faults liable to affect the bearings fitted to rotating machines. This cyclostationarity is reflected in the presence in the vibration signals of a pattern of periodic appearance, exploited by the invention.
  • the fault diagnosis indicators proposed by the invention are also calculated taking into account an estimate of the “real” frequency of the fault, which may be different from the theoretical characteristic frequency of the fault which can be easily obtained from the geometric characteristics of the bearing and its kinematics. This avoids having indicators biased by imprecise knowledge of the frequency of the fault.
  • the diagnostic indicators proposed by the invention provide a reliable and precise estimate of the severity of the faults which may affect the bearings being monitored.
  • the indicators proposed by the invention are also easier to calculate so that the monitoring carried out by virtue of the invention can be carried out in real time.
  • the invention therefore allows reliable and early detection of a defect affecting, if applicable, a bearing, and the location of this defect, with low implementation complexity.
  • the invention therefore makes it possible to minimize the number of sensors required on board a system including several bearings for the monitoring of these bearings, which results in a substantial saving in terms of size compared to certain monitoring techniques. prior art.
  • the invention therefore has a preferred but non-limiting application in the field of aeronautics.
  • the rotating device one bearing of which is monitored, can thus be installed in an aircraft engine, the monitored bearing being, for example, a ball or roller bearing conventionally fitted to such an engine.
  • the invention is also aimed at an aircraft engine comprising at least one bearing fitted to a rotating device of the aircraft engine, at least one accelerometric sensor capable of acquiring a vibratory signal containing a vibratory signature of the bearing, and a monitoring device. of the bearing according to the invention.
  • the vibratory signal has been acquired by the accelerometer sensor during a stationary speed of the rotating device (in other words when its operating parameters, such as for example its average speed of rotation, its pressure, its temperature, its charge, etc., are constant or almost constant).
  • cyclic frequencies are considered “normalized” with respect to the mean frequency of mechanical rotation of the bearing; such standard frequencies are also more commonly referred to as “orders” and make it possible to overcome the variability in the speed of rotation of the rotating device.
  • the processing step comprises a step of spectral whitening of the vibratory signal.
  • Such a whitening of the vibratory signal spectrum is particularly simple and quick to implement. It consists in dividing the Fourier transform of the vibratory signal by its modulus while conserving its phase, and can therefore be carried out blindly, that is to say without knowing the characteristic frequencies of the deterministic component of the signal. Although the random component of the vibration signal, which contains the signature of the defect affecting the bearing, if any, can be weakened by this whitening, this weakening has little impact on the detection of the defect due to the indicator.
  • diagnostic method proposed by the invention which is calculated from the integrated cyclic coherence of the processed vibratory signal. The calculation of the integrated cyclic coherence of a signal in fact implicitly comprises a bleaching operation so that the bleaching operation carried out during the treatment step has no or little effect on the diagnostic indicator. calculated.
  • This way of calculating the integrated cyclic coherence makes it possible to keep the information contained in the vibratory signal acquired by the sensor which relates to the fault, if any, affecting the bearing (and which is described by the cyclic frequency), while maximizing the signal-to-signal ratio. over-noise thanks to the average carried out over a predetermined spectral frequency band.
  • This frequency band is for example taken equal to [k1. ⁇ f; k2. ⁇ f], where ⁇ f denotes the spectral resolution considered to estimate the cyclic coherence, and k1 and k2 two real or integer numbers.
  • ⁇ f denotes the spectral resolution considered to estimate the cyclic coherence
  • k1 and k2 two real or integer numbers.
  • the cyclic correlation of the processed vibratory signal is estimated by means of a Welch estimator.
  • This estimator has the advantage of having remarkable statistical properties and more specifically of providing an estimate with a low quadratic error.
  • Such an estimator is moreover very effective in a context where high rotational speeds are considered (eg when the rotating device is on board an aircraft such as an airplane or a helicopter).
  • estimators can be used to estimate the cyclic correlation of the processed vibratory signal, such as for example an estimator by cyclic periodogram, an estimator by smooth cyclic periodogram, an estimator by cyclic modulation spectrum, etc.
  • the step of estimating the real frequency of the fault comprises the calculation of the integrated cyclic coherence for a plurality of cyclic frequencies included in a defined interval between the theoretical frequency characteristic of the fault minus the defined maximum deviation for this theoretical frequency and the theoretical characteristic frequency of the fault plus the maximum deviation defined for this theoretical frequency, the actual frequency of the fault corresponding to the cyclic frequency among said plurality of cyclic frequencies for which the calculated integrated cyclic coherence is maximum.
  • This embodiment makes it possible, by calculating the integrated cyclic coherence for a plurality of cyclic frequencies included in a defined interval as proposed, to compensate for the lack of knowledge of the exact frequency of the fault affecting the bearing, if applicable.
  • This difficulty is resolved in this embodiment by considering that the real frequency of the fault is that which maximizes the integrated cyclic coherence over the considered interval defined around the theoretical frequency of the fault. It should be noted that in the case of a healthy bearing (i.e.
  • the value obtained for the actual frequency of the fault is not strictly speaking linked to a fault; this is not a problem because, in this case, the indicator calculated as proposed in the invention will have a value which will not exceed the threshold considered during the comparison step (in fact, the probability that another rotating element coincides with the defined interval is zero or almost).
  • two consecutive cyclic frequencies of said plurality of cyclic frequencies are chosen such that the ratio ⁇ / d ⁇ is an integer.
  • the integer M is between 6 and 10.
  • the number of harmonics present in the vibratory signal depends on the impulsivity of the vibratory signal generated by this fault.
  • the impulsivity of a signal reflects the presence of high amplitude peaks in this signal.
  • this impulsivity itself depends, in a known manner, on the severity of the defect affecting the bearing. Therefore, considering a large value of M to determine the diagnostic indicator offers the possibility of taking into account more information on the condition of the bearing and the severity of the fault affecting this bearing. However, a large value of M may increase the estimation error. affecting the calculated indicator and reducing its effectiveness in terms of early detection of the fault.
  • the inventors have observed in practice that a number M between 6 and 10 presents a good compromise making it possible to ensure both early detection of a defect affecting the bearing while ensuring that the diagnostic indicator reflects sufficiently relevant the severity of this defect.
  • the various steps of the monitoring method are determined by computer program instructions.
  • the invention is also aimed at a computer program on an information medium, this program being capable of being implemented in a monitoring device or more generally in a computer, this program comprising instructions adapted to the implementation of the steps of a monitoring method as described above.
  • This program can use any programming language, and be in the form of source code, object code, or intermediate code between source code and object code, such as in a partially compiled form, or in any other. desirable form.
  • the invention also relates to an information medium readable by a computer, and comprising instructions of a computer program as mentioned above.
  • the information medium can be any entity or device capable of storing the program.
  • the support can include a storage means, such as a ROM, for example a CD ROM or a microelectronic circuit ROM, or else a magnetic recording means, for example a floppy disk or a hard disk.
  • the information medium can be a transmissible medium such as an electrical or optical signal, which can be conveyed via an electrical or optical cable, by radio or by other means.
  • the program according to the invention can in particular be downloaded from an Internet type network.
  • the information medium can be an integrated circuit in which the program is incorporated, the circuit being adapted to execute or to be used in the execution of the method in question.
  • the figure 1 shows, in its environment, a monitoring device 1 of a bearing 2 fitted to a rotating mechanical device 3, according to the invention in a particular embodiment.
  • the bearing 2 is a ball bearing and the rotating device 3 equipped with this ball bearing is for example a compressor of a turbojet 4 conforming to invention.
  • the bearing 2 is assumed here to allow the drive in rotation of a shaft of the compressor 3 of the turbojet 4 according to a predefined axis of rotation.
  • the bearing 2 can be a roller bearing
  • the rotating device 3 can be any rotating mechanical device of which an element is driven in rotation by means of the bearing 2
  • other aircraft engines than a turbojet can be considered.
  • the invention can also be applied in a context other than the aeronautical context.
  • the monitoring device 1 is able to carry out monitoring of the bearing 2, from a vibration signal acquired by a sensor 5 fitted to the turbojet 4 over at least a predefined period of time of duration T.
  • the sensor 5 is an accelerometer placed here on the housing of the compressor 3 so as to capture the vibrations emitted by the bearing 2 and more particularly by its elements (vibratory signature of the bearing within the meaning of the invention).
  • a ball bearing is composed of various elements, and more particularly of two coaxial rings (a so-called inner or inner ring and a so-called outer or outer ring) between which are placed balls kept spaced apart by a cage. In this way, the balls can roll between the inner ring and the outer ring.
  • the invention is intended to allow the detection of a defect affecting at least one of these elements.
  • the accelerometer 5 to enable it to acquire a vibratory signal containing the vibratory signature of the bearing 2 does not pose any difficulty for the person skilled in the art, and depends on the rotating mechanical device considered and the context of use of the latter. . It is not described in detail here.
  • the accelerometer 5 is for example placed on the housing of the compressor 3.
  • the accelerometer 5 can be placed in a location further away from the bearing 2 strictly speaking since it makes it possible to acquire a signal containing a vibratory signature of the bearing 2 with a signal-to-noise ratio preferably greater than 5%.
  • the temporal vibratory signal noted X (t) acquired by the accelerometer 5 over the period T when the rotating mechanical device 3 is in rotation, is processed by an acquisition system integrated here in the accelerometer 5 comprising a conditioner, a filter analog anti-aliasing, a sampler block and an analog-to-digital converter.
  • an acquisition system integrated here in the accelerometer 5 comprising a conditioner, a filter analog anti-aliasing, a sampler block and an analog-to-digital converter.
  • Such an acquisition system delivers a digital signal sampled at a predefined sampling frequency Fs and derived from the vibratory signal X (t) acquired by the accelerometer.
  • the sampling frequency Fs is chosen which is sufficiently large to retain the kinematic and dynamic information of the bearing 2. In an aeronautical context such as that envisaged in figure 1 , Fs is taken equal for example to 50kHz.
  • L denotes in an equivalent manner the digital length of the vibratory signal (ie the number of samples Xb (n) derived by the acquisition system from the vibratory signal X (t) acquired by the accelerometer 5).
  • the sampled vibratory signal (vibratory signal within the meaning of the invention), the samples of which are referenced in the remainder of the description by X b (n), where n denotes an integer greater than 1, is then transmitted to the monitoring device 1 for analysis with a view to detecting the presence of a possible defect on the bearing 2.
  • this analysis is intended to be carried out in real time, and the monitoring device 1 is on board the vehicle. aircraft propelled by the turbojet 4 (for example in a computer of this turbojet 4).
  • the monitoring device 1 can be located in a remote device, for example in the example considered here, in a device on the ground capable of communicating via a telecommunications network with the accelerometer 5 or with a computer of the turbojet 4 able to recover the vibratory signal X (t) acquired by the accelerometer 5.
  • the monitoring device 1 has the hardware architecture of a computer, as illustrated on figure 2 .
  • It comprises in particular a processor 6, a RAM 7, a ROM 8, a non-volatile flash memory 9 as well as communication means 10 allowing in particular the monitoring device 1 to communicate with the accelerometer 5 in order to obtain the vibratory signal s (t) generated by the bearing 2 and acquired by the latter.
  • These communication means comprise for example a digital data bus or any other communication interface, in particular a communication interface on a telecommunications network, when the monitoring device 1 is not on board the aircraft propelled by the turbojet. 4.
  • the read only memory 8 of the monitoring device 1 constitutes a recording medium in accordance with the invention, readable by the processor 6 and on which is recorded a computer program PROG in accordance with the invention.
  • the monitoring device 1 also has a notification module 1H, capable of notifying a user or a remote device of the existence of a fault on the bearing 2, if applicable.
  • This notification module 1H can rely in particular on the means of communication 10 of the monitoring device 1 or on input / output means thereof, such as for example a screen or a microphone capable of signaling the detection of a fault on bearing 2 to a user installed near the monitoring device 1.
  • the figure 3 illustrates, in the form of a flowchart, the main steps of a monitoring method according to the invention in a particular embodiment in which it is implemented by the monitoring device 1 to monitor the bearing 2 fitted to the compressor 3 of the turbojet 4.
  • the accelerometer 5 is configured here to acquire the vibration signal X (t) when the compressor 3 is operating in a stationary or quasi-stationary regime: in the example considered here, this means that the speed of rotation of the compressor and / or its load (ie its torque) are constant or quasi-constant (variation of less than 5% of the speed of rotation of the compressor and of its load) .
  • the vibratory signal X (t) acquired by the accelerometer 5 is sampled by the acquisition system of the accelerometer 5, at the sampling frequency Fs.
  • the sampled vibratory signal Xb (n), n 1, ..., L, is transmitted by the acquisition system of the accelerometer 5 to the monitoring device 1, and more particularly to its first obtaining module 1A (step E10).
  • the sampled vibratory signal Xb (n) is referred to more simply in the remainder of the description as a vibratory signal.
  • the elimination of the deterministic component of the vibratory signal is carried out by the processing module 1B via a spectral whitening operation.
  • the processing module 1B first of all calculates the discrete Fourier transform of the vibratory signal Xb (n).
  • the processing module 1B divides the discrete Fourier transform obtained by the modulus of the discrete Fourier transform, and calculates the discrete Fourier transform inverse of the result obtained.
  • This laundering operation has the advantage of being simple to implement.
  • the presence of a defect on bearing 2 is manifested by a component in the cyclostationary signal which results in a periodic auto-covariance function.
  • This property of cyclostationarity is mainly due, on the one hand, to the recurrence of the fault on a periodic basis (linked to the rotation of the rotating device driven by the bearing), and on the other hand, to the presence of a fluctuation between arrival time of impacts caused by sliding of the bearing balls 2.
  • Each defect d likely to affect a bearing is characterized by a specific frequency. No limitation is attached to the nature of the defect, which may relate to one or more components of the bearing. Thus, for example, the defect d may be a defect in the outer ring of the bearing 2, a defect in the inner ring of the bearing 2, a defect in the cage of the bearing 2 or a ball of the bearing 2. Each of these defects is characterized by its own characteristic frequency.
  • the second obtaining module 1C also obtains for each defect d considered, a maximum deviation ⁇ d determined around the theoretical frequency ⁇ d th characteristic of the fault d.
  • ⁇ d A ⁇ d th with A denoting a real constant independent of the fault d.
  • the cyclic coherence of a signal is a statistical measurement which makes it possible to measure for each so-called cyclic frequency the degree of correlation between the signal and this same signal shifted in frequency.
  • a coherence close to one for a cyclic cyclic frequency ⁇ indicates a strong correlation between the components of the signal considered at the frequencies f and f- ⁇ .
  • the cyclic frequency ⁇ is a frequency expressed in hertz.
  • estimators can be used as a variant to estimate the cyclic correlation of the processed vibratory signal, such as for example a cyclic periodogram estimator, a smooth cyclic periodogram estimator, a cyclic modulation spectrum estimator, etc.
  • the cyclic correlation and the cyclic coherence thus calculated are functions of the spectral frequency k ⁇ f and are indexed by the cyclic frequency ⁇ (they are therefore also functions of this cyclic frequency), the latter being expressed in hertz.
  • the spectral frequency brings out the dynamic characteristics of the system considered, while the cyclic frequency is attached to the modulations (mainly linked to the fault affecting, if applicable, bearing 2).
  • the inventors therefore wished to keep the information cyclic carried by the cyclic frequency, but averaged the information carried by the spectral frequency over a frequency band chosen so as to maximize the signal-to-noise ratio.
  • the first 1D calculation module calculates a so-called integrated cyclic coherence by evaluating the mean of the square of the amplitude of the cyclic coherence with respect to the cyclic frequency over a frequency band chosen so as to maximize the signal-to-signal ratio. over-noise (step E50).
  • This frequency band is defined here by the interval [k1 ⁇ f ; k2 ⁇ f ], where k1 and k2 denote two integers.
  • Fs 50kHz
  • the mean calculated to obtain the integrated cyclic correlation has been considered as the square of the amplitude of the cyclic coherence.
  • the integrated cyclic coherence calculated by the first calculation module 1D is then used by the monitoring device 1 to derive a diagnostic indicator of the fault d.
  • the monitoring device via its estimation module 1E, first of all estimates the “real” (or even exact) frequency of the fault d (step E60).
  • this real frequency characteristic of the fault is included in the interval [ ⁇ d th - ⁇ d ; ⁇ d th + ⁇ d ].
  • the estimation module 1E calculates the integrated cyclic coherence for a plurality of cyclic frequencies ⁇ included in this interval.
  • a resolution denoted respectively ⁇ between two consecutive cyclic frequencies considered by the estimation module 1E chosen such that the ratio ⁇ / ⁇ is an integer is assumed. This criterion makes it possible to optimize the cost of calculating the integrated cyclic coherence at the various cyclic frequencies considered.
  • the monitoring device calculates a diagnostic indicator of the fault d (step E70).
  • This indicator advantageously measures the cyclostationarity in the signal X, and consequently the contribution of the defect d associated with the frequency ⁇ d act .
  • M has a predetermined value chosen between 6 and 10.
  • M designates the number of harmonics of the fault d considered.
  • the number of harmonics present in the signal generated by the fault depends on the impulsivity of the signal, which in turn is related to the severity of the fault.
  • M provides more information about the state of the fault.
  • M makes it possible to improve the earliness of the detection of the defect.
  • the inventors have observed that in the case of a distributed fault, few harmonics are present in the integrated cyclic coherence, therefore the sum of a large number of harmonics leads to an increase in the error of estimation, and therefore reduces the efficiency of the calculated indicator in terms of early detection.
  • the estimation module 1F can estimate different diagnostic indicators by considering different values of the integer M.
  • the diagnostic indicator ⁇ X d M thus calculated is supplied to the comparison module 1G of the monitoring device 1.
  • the latter compares the indicator ⁇ X d M with a predetermined alert threshold for the fault d, denoted THR (d) (step E80).
  • This threshold is chosen so as to allow detection of the fault d. It can be determined empirically beforehand by observing healthy bearings and bearings with the defect d, or by a statistical calculation under the assumption of stationarity in the healthy case.
  • the value of the THR (d) threshold must not be chosen too high so as to ensure early detection while not being too low either to avoid multiplying false alarms.
  • the comparison module 1G determines that the diagnostic indicator exceeds the THR (d) threshold (ie is greater than this one) (yes response to test step E90), then it detects that the bearing 2 is affected by the default d (step E100).
  • the monitoring device 1 is able to carry out a differential diagnosis and identify which The fault affects bearing 2, i.e. locating the origin or origins of the damage to bearing 2 (inner ring, outer ring, balls, cage).
  • the vibratory signal X (t) acquired by the accelerometer 5 was acquired during a stationary speed of the compressor 3.
  • the signal originating from a single accelerometer has been considered.
  • the invention can however be applied to several accelerometers. Likewise, it allows several bearings to be monitored simultaneously.
  • the invention thus proposes a robust monitoring technique for bearings integrated in a rotating device.
  • the inventors have been able to observe, during various experiments, that this technique allows early detection of defects affecting bearings, earlier in particular than the technique of the prior art proposed in the document.
  • EP 2,693,176 and this in various contexts (accelerometer placed on the defective bearing or far from it, in the presence or not of electromagnetic interference).
  • the invention has been described with reference to a vibratory signal acquired by means of an accelerometer.
  • the monitoring method proposed by the invention can however also be applied to an acoustic signal, acquired for example by means of a microphone or any other acoustic sensor, and containing an acoustic signature of the bearing that is to be monitored.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Combustion & Propulsion (AREA)
  • Acoustics & Sound (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Description

    Arrière-plan de l'invention
  • L'invention se rapporte au domaine général des machines ou des dispositifs tournants, et concerne plus particulièrement les roulements équipant de telles machines ou dispositifs, comme par exemple les roulements à billes ou à rouleaux.
  • Elle a une application privilégiée mais non limitative dans le domaine de l'aéronautique dans lequel de nombreuses machines tournantes sont utilisées.
  • Les roulements, et en particulier les roulements à billes ou à rouleaux, sont des dispositifs mécaniques tournants largement utilisés dans les machines tournantes et dont le rôle est d'entraîner en rotation des arbres de transmission. Un tel roulement est généralement composé de deux bagues coaxiales (une bague intérieure et une bague extérieure) entre lesquelles sont placées des éléments roulants (des billes ou des rouleaux), maintenus dans une cage. Le mouvement rotatif de l'arbre est assuré par le mouvement des éléments roulants.
  • En raison de leur rôle élémentaire dans les chaînes de transmission mécanique, les roulements sont sollicités en permanence et donc fragilisés. Ils peuvent subir des dégradations de l'état de surface ou de la forme des bagues intérieure ou extérieure, ainsi que des éléments roulants. L'endommagement d'un roulement peut entraîner un arrêt intempestif de la machine dans lequel il est utilisé.
  • Afin d'éviter une telle situation, il est connu de surveiller l'état des roulements équipant une machine tournante pour pouvoir détecter de façon précoce la présence de défauts sur ces roulements. Une des techniques les plus utilisées pour procéder à cette surveillance s'appuie sur une analyse vibratoire de signaux acquis au moyen de capteurs tels que des accéléromètres placés à proximité des roulements.
  • Le principe sur lequel repose l'analyse vibratoire est le suivant : lorsqu'un élément roulant (bille ou rouleau) du roulement entre en contact avec une surface défectueuse du roulement (en raison par exemple d'un écaillage ou d'une fissure à la surface d'une des bagues du roulement), il se produit des chocs périodiques qui produisent une vibration selon une fréquence spécifique, caractéristique du défaut affectant le roulement, et qui viennent changer la structure statistique du signal vibratoire acquis par l'accéléromètre. Une valeur théorique de la fréquence spécifique du défaut peut être obtenue aisément à partir des caractéristiques géométriques du roulement et de sa cinématique. La détection du défaut affectant le roulement se fait alors en identifiant dans le signal vibratoire acquis par l'accéléromètre, une composante à la fréquence spécifique de ce défaut, par exemple en réalisant une analyse spectrale ou une analyse de l'enveloppe du signal vibratoire.
  • Toutefois, les systèmes de surveillance actuels mettant en œuvre une analyse vibratoire se heurtent à plusieurs difficultés pratiques.
  • En effet, il existe aujourd'hui un besoin de surveiller des équipements de plus en plus complexes, notamment dans le domaine de l'aéronautique, qui peuvent comprendre plusieurs éléments tournants en plus des roulements (ex. des compresseurs, des engrenages, des soufflantes, etc.). Cette surveillance, pour des raisons évidentes d'encombrement, doit en outre être réalisée avec un nombre restreint d'accéléromètres. Il en résulte que la composante représentative du défaut affectant le roulement peut être masquée par des bruits provenant d'autres sources de vibrations parasites (i.e. d'autres éléments tournants de l'équipement considéré), de sorte que les méthodes utilisées dans les systèmes de surveillance actuels (par exemple l'analyse d'enveloppe) ne sont pas toujours en mesure de détecter ce défaut. Ceci peut entraîner une détection tardive du défaut voire une incapacité à détecter le défaut.
  • Pour pallier cet inconvénient, il est possible d'avoir recours à des techniques de séparation de sources. Toutefois, ces techniques ont un coût de calcul élevé, et sont souvent impossibles à mettre en œuvre en temps réel.
  • En outre, une difficulté particulière qui apparaît lorsqu'on surveille des roulements par analyse vibratoire est que les fréquences réelles des défauts qui affectent ces roulements peuvent différer des valeurs théoriques calculées pour les besoins de l'analyse vibratoire à partir des caractéristiques géométriques du roulement et de sa cinématique. Ceci est dû au glissement des roulements d'une part, et d'autre part, au fait qu'un endommagement du roulement s'accompagne généralement d'un phénomène de frottement qui a tendance à ralentir la rotation du roulement endommagé et à modifier sa cinématique.
  • Le document EP 2 693 176 décrit un procédé de détection de défauts d'un roulement par analyse vibratoire tentant de résoudre les problèmes précités. Ce procédé s'appuie sur un prétraitement au moyen d'un filtre autorégressif d'un signal vibratoire acquis au moyen d'un accéléromètre placé sur le boîtier du système comprenant les roulements à surveiller, suivi d'une analyse de l'enveloppe du signal vibratoire résultant du prétraitement. Toutefois ce procédé de détection présente une complexité en termes d'implémentation relativement élevée, en raison notamment du prétraitement effectué. En outre, l'efficacité du prétraitement réalisé dépend en grande partie de la détermination de plusieurs paramètres (ex. ordre du filtre autorégressif, etc.) qui peut s'avérer compliquée en pratique. Un autre exemple de procédé est décrit par P. BORGHESANI ET AL: "Application of cepstrum pre-whitening for the diagnosis of bearing faults under variable speed conditions",MECHANICAL SYSTEMS AND SIGNAL PROCESSING, vol. 36, no. 2, 1 avril 2013, pages 370-384.
  • Objet et résumé de l'invention
  • La présente invention a pour but notamment de remédier aux insuffisances précitées de l'état de la technique en proposant un procédé de surveillance d'un roulement équipant un dispositif tournant, comprenant :
    • une étape d'obtention d'un signal vibratoire acquis par un capteur accélérométrique, ledit signal vibratoire contenant une signature vibratoire du roulement ;
    • une étape de traitement du signal vibratoire comprenant l'élimination d'une composante déterministe du signal vibratoire ;
    • une étape d'obtention, pour un défaut déterminé susceptible d'affecter le roulement, d'une fréquence théorique caractéristique de ce défaut et d'une déviation maximale déterminée autour de cette fréquence théorique ;
    • une étape de calcul, en fonction d'une fréquence cyclique, d'une cohérence cyclique dite intégrée du signal vibratoire traité moyennée sur une bande de fréquences spectrales prédéterminée ;
    • une étape d'estimation d'une fréquence réelle du défaut à partir de la cohérence cyclique intégrée, de la fréquence théorique caractéristique du défaut et de la déviation maximale déterminée autour de cette fréquence théorique ;
    • une étape de calcul d'un indicateur de diagnostic du défaut en sommant un nombre entier M de cohérences cycliques intégrées du signal vibratoire évaluées en M fréquences cycliques égales respectivement à M harmoniques de la fréquence réelle estimée du défaut ;
    • une étape de comparaison de l'indicateur de diagnostic du défaut avec un seuil prédéterminé pour ce défaut ; et
    • en cas de dépassement du seuil par l'indicateur de diagnostic, une étape de détection du défaut sur le roulement.
  • Corrélativement, l'invention vise aussi un dispositif de surveillance d'un roulement équipant un dispositif tournant, ce dispositif de surveillance comprenant :
    • un premier module d'obtention, configuré pour obtenir un signal vibratoire acquis par un capteur accélérométrique, ledit signal vibratoire contenant une signature vibratoire du roulement ;
    • un module de traitement du signal vibratoire, configuré pour éliminer une composante déterministe du signal vibratoire ;
    • un deuxième module d'obtention, configuré pour obtenir pour un défaut déterminé susceptible d'affecter le roulement, une fréquence théorique caractéristique de ce défaut et une déviation maximale déterminée autour de cette fréquence théorique ;
    • un premier module de calcul, configuré pour calculer, en fonction d'une fréquence cyclique, une cohérence cyclique dite intégrée du signal vibratoire traité moyennée sur une bande de fréquences spectrales prédéterminée ;
    • un module d'estimation, configuré pour estimer une fréquence réelle du défaut à partir de la cohérence cyclique intégrée, de la fréquence théorique caractéristique du défaut et de la déviation maximale déterminée autour de cette fréquence théorique ;
    • un deuxième module de calcul d'un indicateur de diagnostic du défaut, configuré pour calculer ledit indicateur en sommant un nombre entier M de cohérences cycliques intégrées du signal vibratoire évaluées en M fréquences cycliques égales respectivement à M harmoniques de la fréquence réelle estimée du défaut ;
    • un module de comparaison, configuré pour comparer l'indicateur de diagnostic du défaut avec un seuil prédéterminé pour ce défaut et pour détecter le défaut sur le roulement en cas de dépassement du seuil par l'indicateur de diagnostic.
  • L'invention peut s'appliquer à divers types de roulement, et notamment à un roulement à billes ou à un roulement à rouleaux.
  • L'invention propose donc une technique de surveillance d'un roulement équipant un dispositif tournant basée sur une analyse cyclostationnaire du signal vibratoire acquis par un capteur accélérométrique. Le capteur accélérométrique est placé de telle façon que le signal vibratoire qu'il acquiert contient une signature vibratoire du roulement. Il convient de noter qu'à cet effet, le capteur accélérométrique peut être placé à proximité du roulement ou sur le boîtier du système qui comprend le roulement que l'on cherche à surveiller, mais il peut également être placé en un point éloigné du roulement dès lors qu'un rapport signal-sur-bruit suffisant est assuré (typiquement de l'ordre de 5%).
  • L'analyse cyclostationnaire proposée par l'invention permet d'extraire du signal vibratoire un indicateur de diagnostic d'un défaut susceptible d'affecter le roulement (ex. défaut affectant une bague interne du roulement, défaut affectant une bague externe du roulement, etc.) qui, lorsqu'il est comparé à un seuil déterminé pour ce défaut, permet de détecter si le roulement est affecté ou non dudit défaut. Cet indicateur de diagnostic est avantageusement calculé après avoir appliqué un (pré-)traitement au signal vibratoire consistant à le débarrasser de sa composante déterministe ; cette composante déterministe est typiquement générée par d'autres éléments tournants que le roulement surveillé (qui contribue quant à lui à la composante aléatoire du signal vibratoire), et susceptibles de parasiter le signal vibratoire acquis par le capteur.
  • On note que l'indicateur dérivé conformément à l'invention est déterminé pour un défaut donné susceptible d'affecter le roulement puisqu'il dépend d'une fréquence réelle du défaut estimée à partir de la fréquence théorique caractéristique du défaut et d'une déviation maximale déterminée autour de cette fréquence théorique, chaque défaut se traduisant par une fréquence théorique différente et par une déviation maximale possible de la fréquence réelle du défaut par rapport à cette fréquence théorique. Ainsi, en dérivant plusieurs indicateurs relatifs à une pluralité de défauts susceptibles d'affecter le roulement (ex. un indicateur associé à un défaut sur la bague interne du roulement, un indicateur associé à un défaut sur la bague externe du roulement, un indicateur associé à un défaut sur la cage du roulement, un indicateur associé à un défaut sur un élément roulant du roulement, etc.), il est possible d'identifier aisément le ou les éléments du roulement qui présentent un défaut. Ceci permet d'opérer une maintenance ciblée sur le roulement.
  • De façon avantageuse, la cohérence cyclique du signal vibratoire traité permet de faire ressortir une signature du défaut y compris lorsque celle-ci est « noyée » dans un bruit de fond important (lié par exemple à une pluralité d'éléments tournants défectueux situés à proximité du roulement surveillé). En effet, les inventeurs ont eu l'idée judicieuse d'exploiter la caractéristique cyclostationnaire des signaux de vibrations générés par les défauts susceptibles d'affecter les roulements équipant les machines tournantes. Cette cyclostationnarité se traduit par la présence dans les signaux de vibrations d'un motif d'allure périodique, exploité par l'invention.
  • Il en découle une robustesse des indicateurs de diagnostic de défaut proposés par l'invention. Ces indicateurs sont en outre calculés en tenant compte d'une estimation de la fréquence « réelle » du défaut, qui peut être différente de la fréquence théorique caractéristique du défaut qui peut être obtenue aisément à partir des caractéristiques géométriques du roulement et de sa cinématique. On évite ainsi d'avoir des indicateurs biaisés par une connaissance imprécise de la fréquence du défaut. Au contraire, les indicateurs de diagnostic proposés par l'invention fournissent une estimation fiable et précise de la sévérité des défauts affectant le cas échéant les roulements surveillés.
  • Les indicateurs proposés par l'invention sont de plus faciles à calculer de sorte que la surveillance opérée grâce à l'invention peut être réalisée en temps réel.
  • L'invention permet donc une détection fiable et précoce d'un défaut affectant le cas échéant un roulement, et la localisation de ce défaut, moyennant une faible complexité d'implémentation.
  • En outre, grâce à l'invention, plusieurs roulements peuvent être surveillés en parallèle au moyen d'un seul capteur (par exemple un seul capteur accélérométrique). L'invention permet donc de minimiser le nombre de capteurs requis à bord d'un système incluant plusieurs roulements pour la surveillance de ces roulements, ce qui se traduit par un gain substantiel en termes d'encombrement par rapport à certaines techniques de surveillance de l'art antérieur.
  • L'invention a donc une application privilégiée mais non limitative dans le domaine de l'aéronautique. Le dispositif tournant dont un roulement est surveillé peut ainsi être embarqué dans un moteur d'aéronef, le roulement surveillé étant par exemple un roulement à billes ou à rouleaux équipant classiquement un tel moteur.
  • Ainsi, l'invention vise également un moteur d'aéronef comprenant au moins un roulement équipant un dispositif tournant du moteur d'aéronef, au moins un capteur accélérométrique apte à acquérir un signal vibratoire contenant une signature vibratoire du roulement, et un dispositif de surveillance du roulement selon l'invention.
  • On note toutefois que l'invention peut être appliquée dans de nombreux autres domaines, comme par exemple à des éoliennes, à des moteurs de voiture, à des engrenages, etc.
  • Dans un mode particulier de l'invention, le signal vibratoire a été acquis par le capteur accéléromètrique lors d'un régime stationnaire du dispositif tournant (autrement dit lorsque ses paramètres de fonctionnement, tels que par exemple sa vitesse moyenne de rotation, sa pression, sa température, sa charge, etc., sont constantes ou quasi-constantes).
  • Toutefois cette hypothèse n'est pas limitative et l'invention peut également s'appliquer lorsque le signal vibratoire a été acquis lors d'un régime variable du dispositif tournant, préférentiellement lorsque les variations du régime restent de l'ordre de 15-20%. Dans ce cas, on considère dans l'invention des fréquences cycliques « normalisées » par rapport à la fréquence moyenne de rotation mécanique du roulement ; de telles fréquences normalisées sont aussi désignées plus communément sous l'appellation d' « ordres » et permettent de s'affranchir de la variabilité du régime de rotation du dispositif tournant.
  • Dans un mode particulier de réalisation, l'étape de traitement comprend une étape de blanchiment spectral du signal vibratoire.
  • Un tel blanchiment du spectre du signal vibratoire est particulièrement simple et rapide à mettre en œuvre. Il consiste à diviser la transformée de Fourier du signal vibratoire par son module tout en conservant sa phase, et peut donc être réalisé de manière aveugle, c'est-à-dire sans connaître les fréquences caractéristiques de la composante déterministe du signal. Bien que la composante aléatoire du signal vibratoire, qui contient la signature du défaut affectant le cas échéant le roulement, puisse être affaiblie par ce blanchiment, cet affaiblissement n'a que peu d'impact sur la détection du défaut en raison de l'indicateur de diagnostic proposé par l'invention qui est calculé à partir de la cohérence cyclique intégrée du signal vibratoire traité. Le calcul de la cohérence cyclique intégrée d'un signal comprend en effet implicitement une opération de blanchiment de sorte que l'opération de blanchiment réalisée lors de l'étape de traitement n'a pas ou peu d'effet sur l'indicateur de diagnostic calculé.
  • Dans un mode particulier de réalisation, l'étape de calcul de la cohérence cyclique intégrée comprend :
    • l'estimation, pour une fréquence cyclique donnée, de la corrélation cyclique du signal vibratoire traité en fonction d'une fréquence spectrale ;
    • le calcul, à partir de la corrélation cyclique estimée, de la cohérence cyclique du signal vibratoire traité pour ladite fréquence cyclique donnée en fonction de la fréquence spectrale ; et
    • la moyenne, sur ladite bande de fréquences spectrales prédéterminée, de l'amplitude au carré de la cohérence cyclique du signal vibratoire traité calculée pour ladite fréquence cyclique donnée, le résultat de ladite moyenne fournissant la cohérence cyclique intégrée pour ladite fréquence cyclique donnée.
  • Cette manière de calculer la cohérence cyclique intégrée permet de conserver les informations contenues dans le signal vibratoire acquis par le capteur qui concernent le défaut affectant le cas échéant le roulement (et qui sont décrites par la fréquence cyclique), tout en maximisant le rapport signal-sur-bruit grâce à la moyenne effectuée sur une bande de fréquences spectrales prédéterminée.
  • Cette bande de fréquences est par exemple prise égale à [k1.Δf ;k2.Δf], où Δf désigne la résolution spectrale considérée pour estimer la cohérence cyclique, et k1 et k2 deux nombres réels ou entiers. Une telle plage de fréquences permet, en ajustant de façon adéquate les nombres k1 et k2, d'une part de minimiser la contribution de bruit stationnaire en basse fréquence (i.e. en dessous de k1.Δf) tout en éliminant le repliement cyclique à haute fréquence (i.e. au-dessus de k2.Δf).
  • On note qu'en variante, il est possible dans la moyenne calculée, d'élever l'amplitude de la cohérence cyclique à un ordre différent de l'ordre 2, comme par exemple à l'ordre 1 ou à l'ordre 4. Toutefois, les inventeurs ont constaté que plus l'ordre considéré est élevé, moins la détection du défaut est précoce.
  • Dans un mode particulier de réalisation, la corrélation cyclique du signal vibratoire traité est estimée au moyen d'un estimateur de Welch.
  • Cet estimateur présente l'avantage de disposer de propriétés statistiques remarquables et plus spécifiquement de fournir une estimation ayant une erreur quadratique faible.
  • En outre, il est particulièrement simple et facile à implémenter, et requiert un coût de calcul relativement faible.
  • Un tel estimateur est par ailleurs très efficace dans un contexte où des hautes vitesses de rotation sont considérées (ex. lorsque le dispositif tournant est embarqué à bord d'un aéronef tel qu'un avion ou un hélicoptère).
  • Bien entendu, d'autres estimateurs peuvent être utilisés pour estimer la corrélation cyclique du signal vibratoire traité, comme par exemple un estimateur par périodogramme cyclique, un estimateur par périodogramme cyclique lisse, un estimateur par spectre de modulation cyclique, etc.
  • Dans un mode particulier de réalisation, l'étape d'estimation de la fréquence réelle du défaut comprend le calcul de la cohérence cyclique intégrée pour une pluralité de fréquences cycliques comprises dans un intervalle défini entre la fréquence théorique caractéristique du défaut moins la déviation maximale définie pour cette fréquence théorique et la fréquence théorique caractéristique du défaut plus la déviation maximale définie pour cette fréquence théorique, la fréquence réelle du défaut correspondant à la fréquence cyclique parmi ladite pluralité de fréquences cycliques pour laquelle la cohérence cyclique intégrée calculée est maximale.
  • Ce mode de réalisation permet, en calculant la cohérence cyclique intégrée pour une pluralité de fréquences cycliques comprises dans un intervalle défini tel que proposé, de compenser la non connaissance de la fréquence exacte du défaut affectant le cas échéant le roulement. Cette difficulté est résolue dans ce mode de réalisation en considérant que la fréquence réelle du défaut est celle qui maximise la cohérence cyclique intégrée sur l'intervalle considéré défini autour de la fréquence théorique du défaut. Il convient de noter que dans le cas d'un roulement sain (i.e. ne présentant pas de défaut), la valeur obtenue pour la fréquence réelle du défaut n'est pas liée à proprement parler à un défaut ; ceci n'est pas gênant car, dans ce cas, l'indicateur calculé comme proposé dans l'invention aura une valeur qui ne dépassera par le seuil considéré lors de l'étape de comparaison (en effet, la probabilité pour qu'un autre élément tournant coïncide avec l'intervalle défini est nulle ou quasiment).
  • Dans un mode particulier de réalisation, deux fréquences cycliques consécutives de ladite pluralité de fréquences cycliques, notées respectivement α et α + Δα, sont choisies telles que le rapport α/dα est un nombre entier.
  • L'application de ce critère pour choisir les fréquences cycliques auxquelles la cohérence cyclique intégrée est calculée permet d'optimiser le coût de calcul.
  • Dans un mode particulier de réalisation, le nombre entier M est compris entre 6 et 10.
  • Le nombre d'harmoniques présentes dans le signal vibratoire dépend de l'impulsivité du signal vibratoire généré par ce défaut. De façon connue, l'impulsivité d'un signal traduit la présence de pics d'amplitudes élevées dans ce signal. Or cette impulsivité dépend elle-même de façon connue, de la sévérité du défaut affectant le roulement. Par conséquent, considérer pour déterminer l'indicateur de diagnostic une grande valeur de M offre la possibilité de tenir compte de davantage d'informations sur l'état du roulement et la sévérité du défaut affectant ce roulement. Toutefois, une grande valeur de M peut augmenter l'erreur d'estimation entachant l'indicateur calculé et réduire son efficacité en termes de précocité de la détection du défaut.
  • Au contraire, considérer une faible valeur de M permet d'améliorer la précocité de la détection du défaut.
  • Les inventeurs ont constaté en pratique qu'un nombre M compris entre 6 et 10 présente un bon compromis permettant d'assurer à la fois une détection précoce d'un défaut affectant le roulement tout en s'assurant que l'indicateur de diagnostic reflète de façon suffisamment pertinente la sévérité de ce défaut.
  • Dans un mode particulier de réalisation, le procédé de surveillance comprend en outre une étape de notification du défaut comprenant au moins une information parmi :
    • une identification du roulement défectueux ;
    • une indication d'un élément défectueux sur ledit roulement ; et/ou
    • une indication d'une sévérité du défaut détecté sur le roulement.
  • Ceci permet d'avoir une maintenance adaptée sur le roulement surveillé.
  • Dans un mode particulier de réalisation, les différentes étapes du procédé surveillance sont déterminées par des instructions de programmes d'ordinateurs.
  • En conséquence, l'invention vise aussi un programme d'ordinateur sur un support d'informations, ce programme étant susceptible d'être mis en œuvre dans un dispositif de surveillance ou plus généralement dans un ordinateur, ce programme comportant des instructions adaptées à la mise en œuvre des étapes d'un procédé de surveillance tel que décrit ci-dessus.
  • Ce programme peut utiliser n'importe quel langage de programmation, et être sous la forme de code source, code objet, ou de code intermédiaire entre code source et code objet, tel que dans une forme partiellement compilée, ou dans n'importe quelle autre forme souhaitable.
  • L'invention vise aussi un support d'informations lisible par un ordinateur, et comportant des instructions d'un programme d'ordinateur tel que mentionné ci-dessus.
  • Le support d'informations peut être n'importe quelle entité ou dispositif capable de stocker le programme. Par exemple, le support peut comporter un moyen de stockage, tel qu'une ROM, par exemple un CD ROM ou une ROM de circuit microélectronique, ou encore un moyen d'enregistrement magnétique, par exemple une disquette (floppy disc) ou un disque dur.
  • D'autre part, le support d'informations peut être un support transmissible tel qu'un signal électrique ou optique, qui peut être acheminé via un câble électrique ou optique, par radio ou par d'autres moyens. Le programme selon l'invention peut être en particulier téléchargé sur un réseau de type Internet.
  • Alternativement, le support d'informations peut être un circuit intégré dans lequel le programme est incorporé, le circuit étant adapté pour exécuter ou pour être utilisé dans l'exécution du procédé en question.
  • Brève description des dessins
  • D'autres caractéristiques et avantages de la présente invention ressortiront de la description faite ci-dessous, en référence aux dessins annexés qui en illustrent un exemple de réalisation dépourvu de tout caractère limitatif. Sur les figures :
    • la figure 1 représente de façon schématique un dispositif de surveillance d'un roulement, conforme à l'invention dans un mode particulier de réalisation ;
    • la figure 2 représente l'architecture matérielle du dispositif de surveillance du roulement de la figure 1 ; et
    • la figure 3 illustre, sous forme d'ordinogramme, les principales étapes du procédé de surveillance selon l'invention dans un mode particulier de réalisation dans lequel il est mis en œuvre par le dispositif de surveillance de la figure 1.
    Description détaillée d'un mode de réalisation
  • La figure 1 représente, dans son environnement, un dispositif de surveillance 1 d'un roulement 2 équipant un dispositif mécanique tournant 3, conforme à l'invention dans un mode particulier de réalisation.
  • Dans l'exemple envisagé à la figure 1, le roulement 2 est un roulement à billes et le dispositif tournant 3 équipé de ce roulement à billes est par exemple un compresseur d'un turboréacteur 4 conforme à l'invention. Le roulement 2 est supposé ici permettre l'entraînement en rotation d'un arbre du compresseur 3 du turboréacteur 4 selon un axe de rotation prédéfini.
  • Cet exemple n'est toutefois donné qu'à titre illustratif, et l'invention peut s'appliquer également dans d'autres contextes. Ainsi par exemple, le roulement 2 peut être un roulement à rouleaux, le dispositif tournant 3 peut être n'importe quel dispositif mécanique tournant dont un élément est entraîné en rotation au moyen du roulement 2, et d'autres moteurs d'aéronef qu'un turboréacteur peuvent être envisagés. L'invention peut également s'appliquer dans un autre contexte que le contexte aéronautique.
  • Conformément à l'invention, le dispositif de surveillance 1 est apte à opérer une surveillance du roulement 2, à partir d'un signal vibratoire acquis par un capteur 5 équipant le turboréacteur 4 sur au moins une période de temps prédéfinie de durée T. Le capteur 5 est un accéléromètre placé ici sur le boîtier du compresseur 3 de sorte à capter les vibrations émises par le roulement 2 et plus particulièrement par ses éléments (signature vibratoire du roulement au sens de l'invention).
  • De façon connue en soi, un roulement à billes est composé de divers éléments, et plus particulièrement de deux bagues coaxiales (une bague dite interne ou intérieure et une bague dite externe ou extérieure) entre lesquelles sont placées des billes maintenues espacées par une cage. De la sorte, les billes peuvent rouler entre la bague interne et la bague externe. L'invention est destinée à permettre la détection d'un défaut affectant l'un au moins de ces éléments.
  • Le placement de l'accéléromètre 5 pour lui permettre d'acquérir un signal vibratoire contenant la signature vibratoire du roulement 2 ne pose aucune difficulté pour l'homme du métier, et dépend du dispositif mécanique tournant considéré et du contexte d'utilisation de ce dernier. Il n'est pas décrit en détail ici. Dans le cas d'un arbre d'un compresseur d'un turboréacteur tel qu'envisagé à la figure 1, l'accéléromètre 5 est par exemple placé sur le boîtier du compresseur 3.
  • En variante, l'accéléromètre 5 peut être placé en un endroit plus éloigné du roulement 2 à proprement parler dès lors qu'il permet d'acquérir un signal contenant une signature vibratoire du roulement 2 avec un rapport signal-sur-bruit préférentiellement supérieur à 5%.
  • Le signal vibratoire temporel noté X(t) acquis par l'accéléromètre 5 sur la période T lorsque le dispositif mécanique tournant 3 est en rotation, est traité par un système d'acquisition intégré ici dans l'accéléromètre 5 comprenant un conditionneur, un filtre analogique d'anti-repliement, un échantillonneur bloquer et un convertisseur analogique-numérique. Un tel système d'acquisition est connu en soi, et n'est pas décrit en détail ici. Il délivre un signal numérique échantillonné à une fréquence d'échantillonnage prédéfinie Fs et dérivé du signal vibratoire X(t) acquis par l'accéléromètre. On choisit la fréquence d'échantillonnage Fs suffisamment grande pour conserver les informations cinématiques et dynamiques du roulement 2. Dans un contexte aéronautique tel que celui envisagé à la figure 1, Fs est prise égale par exemple à 50kHz. On suppose ici que la période T sur laquelle est acquis le signal vibratoire X est un multiple L de la période d'échantillonnage soit : T=L/Fs. L désigne de façon équivalente la longueur numérique du signal vibratoire (i.e. le nombre d'échantillons Xb(n) dérivés par le système d'acquisition à partir du signal vibratoire X(t) acquis par l'accéléromètre 5).
  • Le signal vibratoire échantillonné (signal vibratoire au sens de l'invention), dont les échantillons sont référencés dans la suite de la description par Xb(n), où n désigne un entier supérieur à 1, est ensuite transmis au dispositif de surveillance 1 pour analyse en vue de détecter la présence d'un éventuel défaut sur le roulement 2. Dans le mode de réalisation décrit ici, cette analyse est destinée à être réalisée en temps réel, et le dispositif de surveillance 1 est embarqué à bord de l'aéronef propulsé par le turboréacteur 4 (par exemple dans un calculateur de ce turboréacteur 4).
  • En variante, le dispositif de surveillance 1 peut être localisé dans un dispositif distant, par exemple dans l'exemple envisagé ici, dans un dispositif au sol apte à communiquer via un réseau de télécommunications avec l'accéléromètre 5 ou avec un calculateur du turboréacteur 4 apte à récupérer le signal vibratoire X(t) acquis par l'accéléromètre 5.
  • Dans le mode de réalisation décrit ici, le dispositif de surveillance 1 a l'architecture matérielle d'un ordinateur, telle qu'illustrée à la figure 2 .
  • Il comprend notamment un processeur 6, une mémoire vive 7, une mémoire morte 8, une mémoire flash non volatile 9 ainsi que des moyens de communication 10 permettant notamment au dispositif de surveillance 1 de communiquer avec l'accéléromètre 5 pour obtenir le signal vibratoire s(t) généré par le roulement 2 et acquis par ce dernier. Ces moyens de communication comprennent par exemple un bus de données numériques ou toute autre interface de communication, notamment une interface de communication sur un réseau de télécommunications, lorsque le dispositif de surveillance 1 ne se trouve pas à bord de l'aéronef propulsé par le turboréacteur 4.
  • La mémoire morte 8 du dispositif de surveillance 1 constitue un support d'enregistrement conforme à l'invention, lisible par le processeur 6 et sur lequel est enregistré un programme d'ordinateur PROG conforme à l'invention.
  • Le programme d'ordinateur PROG définit des modules fonctionnels et logiciels ici, configurés pour mettre en œuvre un procédé de surveillance du roulement 2 conforme à l'invention. Ces modules fonctionnels s'appuient sur et/ou commandent les éléments matériels 6-10 du dispositif de surveillance 1 cités précédemment. Ils comportent notamment ici, comme illustré sur la figure 1 :
    • un premier module d'obtention 1A, configuré pour obtenir le signal vibratoire s(t) acquis par l'accéléromètre 5 ;
    • un module de traitement 1B du signal vibratoire s(t), configuré pour éliminer une composante déterministe du signal vibratoire ;
    • un deuxième module d'obtention 1C, configuré pour obtenir pour un défaut déterminé noté d susceptible d'affecter le roulement 2, une fréquence théorique caractéristique de ce défaut notée β d th
      Figure imgb0001
      et une déviation maximale notée δβd déterminée autour de cette fréquence théorique. Cette déviation maximale caractérise l'écart maximale possible que l'on s'accorde entre la fréquence théorique caractéristique du défaut et la fréquence réelle ou exacte de ce défaut ;
    • un premier module de calcul 1D, configuré pour calculer, en fonction de la fréquence cyclique notée ici α, une cohérence cyclique dite intégrée du signal vibratoire traité, moyennée sur une bande de fréquences spectrales prédéterminée choisie ici de sorte à maximiser le rapport signal-sur-bruit ;
    • un module d'estimation 1E, configuré pour estimer la fréquence réelle du défaut notée β d act
      Figure imgb0002
      à partir de la cohérence cyclique intégrée calculée par le premier module de calcul 1D, de la fréquence théorique β d th
      Figure imgb0003
      caractéristique du défaut et de la déviation maximale δβd déterminée autour de cette fréquence théorique ;
    • un deuxième module de calcul 1F d'un indicateur de diagnostic du défaut, configuré pour calculer cet indicateur en sommant un nombre entier M de cohérences cycliques intégrées du signal vibratoire traité évaluées en M fréquences cycliques égales respectivement à M harmoniques de la fréquence réelle estimée du défaut ;
    • un module de comparaison 1G, configuré pour comparer l'indicateur de diagnostic du défaut calculé par le deuxième module de calcul 1F avec un seuil prédéterminé pour ce défaut et pour détecter le défaut sur le roulement en cas de dépassement du seuil par l'indicateur de diagnostic.
  • Dans le mode de réalisation décrit ici, le dispositif de surveillance 1 dispose en outre d'un module de notification 1H, apte à notifier un utilisateur ou un équipement distant de l'existence d'un défaut sur le roulement 2 le cas échéant. Ce module de notification 1H peut s'appuyer notamment sur les moyens de communication 10 du dispositif de surveillance 1 ou sur des moyens d'entrée/sortie de celui-ci, comme par exemple un écran ou un microphone apte à signaler la détection d'un défaut sur le roulement 2 à un utilisateur installé à proximité du dispositif de surveillance 1.
  • Les fonctions de ces différents modules sont décrites plus en détail maintenant en référence aux étapes du procédé de surveillance selon l'invention.
  • La figure 3 illustre, sous forme d'ordinogramme, les principales étapes d'un procédé de surveillance selon l'invention dans un mode particulier de réalisation dans lequel il est mis en œuvre par le dispositif de surveillance 1 pour surveiller le roulement 2 équipant le compresseur 3 du turboréacteur 4.
  • On suppose ici que l'accéléromètre 5 est configuré de sorte à acquérir sur une période de temps de durée T=L/Fs, un signal vibratoire X(t) d'accélération généré par le roulement 2 alors que le compresseur 3 est en rotation entraînée par le roulement 2. L'accéléromètre 5 est configuré ici pour acquérir le signal vibratoire X(t) lorsque le compresseur 3 fonctionne dans un régime stationnaire ou quasi-stationnaire : dans l'exemple envisagé ici, cela signifie que la vitesse de rotation du compresseur et/ou sa charge (i.e. son couple) sont constantes ou quasi-constantes (variation de moins de 5% de la vitesse de rotation du compresseur et de sa charge). On note que la durée T ne doit pas être choisie trop courte de sorte à assurer la précision des indicateurs de diagnostic dérivés par l'invention. Pour une fréquence d'échantillonnage Fs=50kHz, une durée T=2s représente un bon compromis.
  • Le signal vibratoire X(t) acquis par l'accéléromètre 5 est échantillonné par le système d'acquisition de l'accéléromètre 5, à la fréquence d'échantillonnage Fs. Le signal vibratoire échantillonné Xb(n), n=1,...,L, est transmis par le système d'acquisition de l'accéléromètre 5 au dispositif de surveillance 1, et plus particulièrement à son premier module d'obtention 1A (étape E10). Le signal vibratoire échantillonné Xb(n) est désigné plus simplement dans la suite de la description par signal vibratoire.
  • Le signal vibratoire Xb(n), n=1,...,L est fourni par le module d'obtention 1A au module de traitement 1B du dispositif de surveillance 1. De façon connue en soi, le signal vibratoire Xb(n), n=1,...,L comprend une composante déterministe et une composante aléatoire ; c'est dans cette composante aléatoire que se manifestent le cas échéant, les vibrations liées au(x) défaut(s) affectant le roulement 2.
  • Pour mieux faire ressortir les vibrations liées au(x) défaut(s) affectant le cas échéant le roulement 2, le module de traitement 1B traite le signal vibratoire Xb(n), n=1,...,L afin d'éliminer sa composante déterministe (étape E20).
  • Dans le mode de réalisation décrit ici, l'élimination de la composante déterministe du signal vibratoire est réalisée par le module de traitement 1B via une opération de blanchiment spectral. Pour réaliser cette opération, le module de traitement 1B calcule tout d'abord la transformée de Fourier discrète du signal vibratoire Xb(n). On note DFTn→m {Xb(n)} cette transformée de Fourier discrète. Elle est définie de façon connue par : DFT n k W X n = n = 0 W 1 X n e j 2 π nΔt kΔf
    Figure imgb0004
    où : Δf = 1 Δt . W = F S W
    Figure imgb0005
    désigne la résolution fréquentielle, W représentant la taille de la fenêtre sur laquelle est calculée la transformée de Fourier discrète, et Δt = 1 F S
    Figure imgb0006
    désigne la résolution temporelle.
  • Puis le module de traitement 1B divise la transformée de Fourier discrète obtenue par le module de la transformée de Fourier discrète, et calcule la transformée de Fourier discrète inverse du résultat obtenu. Autrement dit, le module de traitement 1B calcule le signal X(n) suivant, pour n=1,...,L (signal vibratoire traité au sens de l'invention) : X n = IDFT m n L DFT n m Xb n DFT n m { Xb n
    Figure imgb0007
    IDFT m n L
    Figure imgb0008
    désigne la transformée de Fourier discrète inverse définie comme suit IDFT n k W X n = 1 W n = 0 W 1 X n e + j 2 π nΔt kΔf
    Figure imgb0009
    en considérant les mêmes notations utilisées lors de la définition de la transformée de Fourier discrète.
  • De la sorte, la phase du signal vibratoire Xb(n) est conservée.
  • Cette opération de blanchiment présente l'avantage d'être simple à mettre en œuvre.
  • En variante, d'autres techniques peuvent être appliquées par le module de traitement 1B pour éliminer la composante déterministe du signal Xb(n). Une telle technique est décrite par exemple dans le document de N. Sawalhi et al. Intitulé «Signal pre-whitening using cepstrum editing (liftering) to enhance fault detection in rolling element bearings », Proceedings of the 24th International Congress on Condition Monitoring and Diagnostics Engineering Management (COMADEM), 2011. Elle consiste à annuler tout le cepstre complexe du signal vibratoire en conservant uniquement la composante du signal relative à la fréquence zéro.
  • Le signal vibratoire traité X(n), n=1,...,L (et blanchi) est transmis au premier module de calcul 1D du dispositif de surveillance 1 pour qu'il analyse la cyclostationnarité du signal. En effet, la présence d'un défaut sur le roulement 2 se manifeste par une composante dans le signal cyclostationnaire qui se traduit par une fonction d'auto-covariance périodique. Cette propriété de cyclostationnarité est due principalement d'une part, à la récurrence du défaut sur une base périodique (liée à la rotation du dispositif tournant entraîné par le roulement), et d'autre part, à la présence d'une fluctuation entre les temps d'arrivée des impacts causés par le glissement des billes du roulement 2.
  • Chaque défaut d susceptible d'affecter un roulement est caractérisé par une fréquence spécifique. Aucune limitation n'est attachée à la nature du défaut d, celui-ci pouvant porter sur un ou plusieurs éléments du roulement. Ainsi par exemple, le défaut d peut être un défaut de la bague externe du roulement 2, un défaut de la bague interne du roulement 2, un défaut de la cage du roulement 2 ou d'une bille du roulement 2. Chacun de ces défauts est caractérisée par une fréquence caractéristique qui lui est propre.
  • Cette fréquence caractéristique du défaut d considéré peut être estimée (obtenue) aisément de façon théorique par le deuxième module d'obtention 1C du dispositif de surveillance 1 à partir de la connaissance des caractéristiques géométriques du roulement 2 et de sa cinématique (étape E30). On note β d th
    Figure imgb0010
    la fréquence théorique caractéristique du défaut d. Les caractéristiques géométriques du roulement 2 peuvent être obtenues aisément à partir de la fiche technique du roulement 2. La cinématique du roulement (i.e. sa vitesse de rotation) peut être obtenue via un capteur placé de façon appropriée au sein du compresseur 3, de façon connue en soi. Ainsi, à titre illustratif :
    • la fréquence théorique caractéristique fd=bext d'un défaut d de la bague externe du roulement est donnée par : f d = bext = fr 2 N B 1 D B cosψ D p
      Figure imgb0011
    • la fréquence théorique caractéristique fd=bint d'un défaut d de la bague interne du roulement est donnée par : f d = bint = f r 2 N B 1 + D B cosψ D p
      Figure imgb0012
    • la fréquence théorique caractéristique fd=c d'un défaut d de la cage du roulement est donnée par : f d = c = f r 2 1 D B cosψ D p
      Figure imgb0013
    • la fréquence théorique caractéristique fd=bille d'un défaut d d'une bille du roulement est donnée par : f d = bille = f r 2 N B 1 D B cosψ D p 2
      Figure imgb0014
      où :
      • fr désigne la fréquence de rotation mécanique du roulement 2 ;
      • NB désigne le nombre de billes du roulement 2 ;
      • DB désigne le diamètre des billes du roulement 2 ;
      • DP désigne le diamètre moyen du roulement 2 ; et
      • ψ désigne l'angle de contact du roulement 2.
  • Le deuxième module d'obtention 1C obtient également pour chaque défaut d considéré, une déviation maximale δβd déterminée autour de la fréquence théorique β d th
    Figure imgb0015
    caractéristique du défaut d. Dans l'exemple envisagé ici, on suppose par souci de simplification que cette déviation maximale est la même en pourcentage quel que soit le défaut d envisagé, soit δβ d = d th
    Figure imgb0016
    avec A désignant une constante réelle indépendante du défaut d. En pratique, les inventeurs ont constaté qu'il suffit de choisir la déviation maximale égale à quelques pourcents de la fréquence théorique caractéristique du défaut ; par exemple A = 3%.
  • Puis, dans le mode de réalisation décrit ici, le premier module de calcul 1D du dispositif de surveillance 1 calcule la cohérence cyclique du signal vibratoire traité X(n), n=1,...,L (étape E40).
  • De façon connue, la cohérence cyclique d'un signal est une mesure statistique qui permet de mesurer pour chaque fréquence dite cyclique le degré de corrélation entre le signal et ce même signal décalé en fréquence. Une cohérence proche de un pour une fréquence cyclique cyclique α indique une forte corrélation entre les composantes du signal considéré aux fréquences f et f- α.
  • Pour calculer la cohérence cyclique du signal X(n) à la fréquence cyclique α, n=1,..,L, le premier module de calcul 1D utilise ici un estimateur de Welch de la corrélation cyclique du signal, qui s'écrit de la façon suivante, à un facteur près : S 2 X α ^ kΔf = s S 1 DFT n k Nw w s n X n e jπαnΔt * DFT n k Nw w s n X n e jπαnΔt
    Figure imgb0017
    où :
    • * désigne l'opérateur de conjugaison ;
    • w(n) désigne une fenêtre glissante (ex. fenêtre de Hanning, demi-sinus, etc) et ws (n) = ws (n - sR)est la version décalée de la fenêtre glissante avec 1<R<Nw, Nw désignant un nombre entier (ex. une puissance de 2) représentant la taille de la fenêtre de Welch, Nw-R désignant le recouvrement entre les fenêtres ;
    • S est le plus grand entier inférieur ou égal à (L-Nw)/R+1 ;
    • Δf désigne la résolution fréquentielle ou spectrale, qui est égale à 1/(Δt. Nw) . On choisira préférentiellement la taille de la fenêtre de Welch Nw de sorte à avoir une résolution spectrale de l'ordre de quelques centaines de hertz dans le cas d'application envisagé ici ;
    • k est un nombre entier, désignant le canal spectral considéré.
  • On considère ici comme résolution de la fréquence cyclique Δα = 1/T pour calculer la corrélation cyclique.
  • On note que dans l'estimateur de Welch défini ci-dessus, la fréquence cyclique α est une fréquence exprimée en hertz.
  • Bien entendu, d'autres estimateurs peuvent être utilisés en variante pour estimer la corrélation cyclique du signal vibratoire traité, comme par exemple un estimateur par périodogramme cyclique, un estimateur par périodogramme cyclique lisse, un estimateur par spectre de modulation cyclique, etc.
  • Puis le module de calcul 1D dérive à partir de l'estimateur de Welch ainsi calculé la cohérence cyclique du signal X(n), n=1,...,L de la façon suivante : γ 2 X α ^ kΔf = S 2 x α ^ kΔf S 2 X 0 ^ kΔf α / 2 S 2 x 0 ^ kΔf + α / 2 1 / 2
    Figure imgb0018
  • On note que la corrélation cyclique et la cohérence cyclique ainsi calculées sont des fonctions de la fréquence spectrale kΔf et sont indexées par la fréquence cyclique α (elles sont donc également fonctions de cette fréquence cyclique), celle-ci étant exprimée en hertz. La fréquence spectrale fait ressortir des caractéristiques dynamiques du système considéré, tandis que la fréquence cyclique s'attache aux modulations (principalement liées au défaut affectant le cas échéant le roulement 2). Les inventeurs ont donc souhaité garder l'information cyclique portée par la fréquence cyclique, mais moyenné l'information portée par la fréquence spectrale sur une bande de fréquence choisie de sorte à maximiser le rapport signal-sur-bruit.
  • A cet effet, le premier module de calcul 1D calcule une cohérence cyclique dite intégrée en évaluant la moyenne du carré de l'amplitude de la cohérence cyclique par rapport à la fréquence cyclique sur une bande de fréquences choisie de sorte à maximiser le rapport signal-sur-bruit (étape E50). Cette bande de fréquences est définie ici par l'intervalle [k1Δf ;k2Δf], où k1 et k2 désignent deux entiers. La cohérence cyclique intégrée calculée par le premier module de calcul 1D est ainsi donnée par l'expression suivante : I 2 X α ^ = 1 k 1 k 1 + 1 k 1 k 2 γ 2 X α ^ kΔf 2
    Figure imgb0019
    Dans l'exemple envisagé ici, pour une fréquence d'échantillonnage Fs=50kHz, on choisit k1Δf∼1kHz et k2Δf∼20kHz. Ces valeurs permettent de minimiser la contribution du bruit stationnaire en basse fréquence (en dessous de k1Δf) et d'éliminer le repliement cyclique à haute fréquence (au dessus de k2Δf).
  • Bien entendu ces valeurs ne sont données qu'à titre illustratif et d'autres valeurs peuvent être considérées. De manière générale, on s'assurera préférentiellement que k2Δf < Fs - αmax αmax désigne la valeur cyclique maximale recherchée et donc dépend de la fréquence caractéristique théorique du défaut et de son ordre d'harmoniques (on tient compte ici du principe que l'information utile pour diagnostiquer un défaut du roulement ne se trouve généralement pas aux extrémités).
  • On note par ailleurs que dans le mode de réalisation décrit ici, on a considéré dans la moyenne calculée pour obtenir la corrélation cyclique intégrée le carré de l'amplitude de la cohérence cyclique. En variante, on peut envisager d'élever la cohérence cyclique à un autre ordre que l'ordre 2, par exemple à l'ordre 1 ou à l'ordre 4.
  • La cohérence cyclique intégrée calculée par le premier module de calcul 1D est alors utilisée par le dispositif de surveillance 1 pour dériver un indicateur de diagnostic du défaut d.
  • Plus spécifiquement, à cet effet, le dispositif de surveillance 1, via son module d'estimation 1E, estime tout d'abord la fréquence « réelle » (ou encore exacte) du défaut d (étape E60). Comme mentionné précédemment et par définition, cette fréquence réelle caractéristique du défaut est comprise dans l'intervalle [ β d th δβ d
    Figure imgb0020
    ; β d th + δβ d
    Figure imgb0021
    ]. Pour estimer cette fréquence réelle, le module d'estimation 1E calcule la cohérence cyclique intégrée pour une pluralité de fréquences cycliques α comprises dans cet intervalle. Dans le mode de réalisation décrit ici, on suppose une résolution notées respectivement Δα entre deux fréquences cycliques consécutives considérées par le module d'estimation 1E choisie telle que le rapport α/Δα soit un nombre entier. Ce critère permet d'optimiser le coût de calcul de la cohérence cyclique intégrée aux différentes fréquences cycliques considérées.
  • Puis le module d'estimation 1E estime la fréquence réelle du défaut, notée β d act
    Figure imgb0022
    , comme étant la valeur de la fréquence cyclique correspondant à la valeur maximale des cohérences cycliques intégrées calculées, soit : β d act = argmax β d th δβ d α β d th + δβ d I 2 X α ^
    Figure imgb0023
    On note que dans le cas où le roulement 2 ne présente pas le défaut d, la valeur β d act
    Figure imgb0024
    n'a pas de signification à proprement parler et ne présente pas la fréquence d'un défaut du roulement. Ceci ne constitue pas un problème en soi ; en effet, il est fortement improbable que le roulement 2 soit affecté d'un défaut distinct du défaut d ayant une fréquence réelle comprise dans l'intervalle [ β d th δβ d
    Figure imgb0025
    ; β d th + δβ d
    Figure imgb0026
    ].
  • Puis le dispositif de surveillance 1, par l'intermédiaire de son deuxième module de calcul 1F, calcule un indicateur de diagnostic du défaut d (étape E70). Cet indicateur, noté μ X d
    Figure imgb0027
    , est la somme des cohérences cycliques intégrées évaluées en M fréquences cycliques égales respectivement aux M harmoniques de la fréquence réelle β d act
    Figure imgb0028
    estimée du défaut, M désignant un nombre entier, soit : μ X d M = m = 1 M I 2 X d act ^
    Figure imgb0029
  • Cet indicateur mesure avantageusement la cyclostationnarité dans le signal X, et par conséquent la contribution du défaut d associé à la fréquence β d act
    Figure imgb0030
    .
  • Dans le mode de réalisation décrit ici, M a une valeur prédéterminée choisie entre 6 et 10.
  • On note que M désigne le nombre d'harmoniques du défaut d considérées. Le nombre d'harmoniques présent dans le signal généré par le défaut dépend de l'impulsivité du signal, qui à son tour est liée à la sévérité du défaut. Ainsi, prendre une grande valeur de M fournit plus d'informations sur l'état du défaut. Toutefois, prendre une plus petite valeur de M permet d'améliorer la précocité de la détection du défaut. En effet, les inventeurs ont constaté que dans le cas d'un défaut distribué, peu d'harmoniques sont présentes dans la cohérence cyclique intégrée, donc la somme d'un grand nombre d'harmoniques conduit à une augmentation de l'erreur d'estimation, et donc réduit l'efficacité de l'indicateur calculé en termes de précocité de détection.
  • Dans une variante de réalisation, le module d'estimation 1F peut estimer différents indicateurs de diagnostic en considérant différentes valeurs de l'entier M.
  • L'indicateur de diagnostic μ X d M
    Figure imgb0031
    ainsi calculé est fourni au module de comparaison 1G du dispositif de surveillance 1. Celui-ci compare alors l'indicateur μ X d M
    Figure imgb0032
    avec un seuil d'alerte prédéterminé pour le défaut d, noté THR(d) (étape E80). Ce seuil est choisi de sorte à permettre la détection du défaut d. Il peut être déterminé préalablement de façon empirique en observant des roulements sains et des roulements présentant le défaut d, ou par un calcul statistique sous l'hypothèse de stationnarité dans le cas sain. La valeur du seuil THR(d) ne doit pas être choisie trop élevée de sorte à assurer une détection précoce tout en n'étant pas trop basse non plus pour éviter de multiplier les fausses alarmes.
  • Si le module de comparaison 1G détermine que l'indicateur de diagnostic dépasse le seuil THR(d) (i.e. est supérieur à celui-ci) (réponse oui à l'étape test E90), alors il détecte que le roulement 2 est affecté du défaut d (étape E100).
  • Dans le mode de réalisation décrit ici, cette détection est alors notifiée par le module de notification 1H via un message d'alarme, comprenant tout ou partie des informations suivantes :
    • identification du roulement défectueux (à savoir ici le roulement 2) ;
    • indication de l'élément défectueux sur le roulement (dépend du défaut d détecté) ; et
    • indication de la sévérité du défaut détecté sur le roulement 2 (donnée par la valeur de l'indicateur μ X d M
      Figure imgb0033
      ).
  • Puis la surveillance du roulement 2 est reprise (répétition des étapes E10 à E100).
  • On note que ce qui vient d'être décrit pour le défaut d peut être réalisé pour différents défauts susceptibles d'affecter le roulement 2. De la sorte, le dispositif de surveillance 1 est en mesure de réaliser un diagnostic différentiel et d'identifier quel défaut affecte le roulement 2, c'est-à-dire de localiser l'origine ou les origines des dégradations subies par le roulement 2 (bague interne, bague externe, billes, cage).
  • Dans le mode de réalisation décrit ici, le signal vibratoire X(t) acquis par l'accéléromètre 5 a été acquis lors d'un régime stationnaire du compresseur 3.
  • Dans un autre mode de réalisation de l'invention, on peut considérer un signal vibratoire acquis lorsque le régime de rotation du compresseur 3 connait une certaine variabilité, typiquement jusqu'à 15-20%.
  • Dans un tel contexte, on considère préférentiellement une fréquence cyclique normalisée notée α̃ par rapport à la fréquence de rotation fr moyenne en hertz du roulement 2 sur la période d'acquisition T considérée, soit : α ˜ = α f r
    Figure imgb0034
    Une telle fréquence normalisée est également appelée plus communément « ordre ».
  • Dans ce mode de réalisation, lors de l'étape E40, le premier module de calcul 1D utilise alors l'estimateur de Welch de la corrélation cyclique du signal à phase corrigée, qui s'écrit de la façon suivante, à un facteur près : S 2 X α ˜ ^ kΔf = s S 1 DFT n k Nw w s n X n e α ˜ θ n / 2 * DFT nΔt kΔf Nw w s n X n e α ˜ θ n / 2
    Figure imgb0035
    où θ(n) désigne la mesure angulaire d'un arbre de référence (arbre de rotation du compresseur 3 ici) qui permet d'évaluer la fréquence de rotation fr du roulement.
  • Le premier module de calcul 1D dérive ensuite à partir de cet estimateur de Welch la cohérence cyclique à phase corrigée du signal X(n), n=1,...,L de la façon suivante : γ 2 X α ˜ ^ kΔf = S 2 X α ˜ ^ kΔf S 2 X 0 ^ kΔf α ˜ . f r / 2 S 2 X 0 ^ kΔf + α ˜ . f r / 2 1 / 2
    Figure imgb0036
  • Les autres étapes du procédé restent inchangées dès lors que l'on considère des fréquences normalisées au lieu de fréquences exprimées en hertz. Ainsi, les fréquences théoriques des défauts sont les mêmes que celles introduites précédemment au facteur fr près (i.e. on divise ces fréquences théoriques par fr).
  • Par ailleurs, dans les modes de réalisation décrits ici, on a considéré par souci de simplification, le signal issu d'un seul accéléromètre. L'invention peut toutefois s'appliquer à plusieurs accéléromètres. De même elle permet de surveiller simultanément plusieurs roulements.
  • L'invention propose ainsi une technique de surveillance robuste des roulements intégrés dans un dispositif tournant. Les inventeurs ont pu constater, lors de diverses expériences, que cette technique permet une détection précoce des défauts affectant les roulements, plus précoce notamment que la technique de l'art antérieur proposée dans le document EP 2 693 176 , et ce dans divers contextes (accéléromètre placé sur le roulement défectueux ou éloigné de celui-ci, en présence ou non d'interférences électromagnétiques).
  • On note que l'invention a été décrite en référence à un signal vibratoire acquis au moyen d'un accéléromètre. La méthode de surveillance proposée par l'invention peut toutefois également s'appliquer à un signal acoustique, acquis par exemple au moyen d'un microphone ou de tout autre capteur acoustique, et contenant une signature acoustique du roulement que l'on souhaite surveiller.

Claims (15)

  1. Procédé de surveillance d'un roulement (2) équipant un dispositif tournant (3), comprenant :
    - une étape d'obtention (E10) d'un signal vibratoire acquis par un capteur accélérométrique, ledit signal vibratoire comprenant une signature vibratoire du roulement ;
    - une étape de traitement (E20) du signal vibratoire comprenant l'élimination d'une composante déterministe du signal vibratoire ;
    - une étape d'obtention (E30), pour un défaut déterminé susceptible d'affecter le roulement, d'une fréquence théorique caractéristique de ce défaut et d'une déviation maximale déterminée autour de cette fréquence théorique ;
    - une étape de calcul (E40-E50), en fonction d'une fréquence cyclique, d'une cohérence cyclique dite intégrée du signal vibratoire traité moyennée sur une bande de fréquences spectrales prédéterminée ;
    - une étape d'estimation (E60) d'une fréquence dite réelle du défaut à partir de la cohérence cyclique intégrée, de la fréquence théorique caractéristique du défaut et de la déviation maximale déterminée autour de cette fréquence théorique ;
    - une étape de calcul (E70) d'un indicateur de diagnostic du défaut en sommant un nombre entier M de cohérences cycliques intégrées du signal vibratoire évaluées en M fréquences cycliques égales respectivement à M harmoniques de la fréquence réelle estimée du défaut ;
    - une étape de comparaison (E80-E90) de l'indicateur de diagnostic du défaut avec un seuil prédéterminé pour ce défaut ; et
    - en cas de dépassement du seuil par l'indicateur de diagnostic, une étape de détection (E100) du défaut sur le roulement.
  2. Procédé de surveillance selon la revendication 1 dans lequel l'étape de traitement comprend une étape de blanchiment spectral du signal vibratoire.
  3. Procédé de surveillance selon la revendication 1 ou 2 dans lequel l'étape de calcul de la cohérence cyclique intégrée comprend :
    - l'estimation (E40), pour une fréquence cyclique donnée, de la corrélation cyclique du signal vibratoire traité en fonction d'une fréquence spectrale ;
    - le calcul (E40), à partir de la corrélation cyclique estimée, de la cohérence cyclique du signal vibratoire traité pour ladite fréquence cyclique donnée en fonction de la fréquence spectrale ; et
    - la moyenne (E50), sur ladite bande de fréquences spectrales prédéterminée, de l'amplitude au carré de la cohérence cyclique du signal vibratoire traité calculée pour ladite fréquence cyclique donnée, le résultat de ladite moyenne fournissant la cohérence cyclique intégrée pour ladite fréquence cyclique donnée.
  4. Procédé de surveillance selon la revendication 3 dans lequel la corrélation cyclique du signal vibratoire traité est estimée au moyen d'un estimateur de Welch.
  5. Procédé de surveillance selon l'une quelconque des revendications 1 à 4 dans lequel l'étape d'estimation (E60) de la fréquence réelle du défaut comprend le calcul de la cohérence cyclique intégrée pour une pluralité de fréquences cycliques comprises dans un intervalle défini entre la fréquence théorique caractéristique du défaut moins la déviation maximale définie pour cette fréquence théorique et la fréquence théorique caractéristique du défaut plus la déviation maximale définie pour cette fréquence théorique, la fréquence réelle du défaut correspondant à la fréquence cyclique parmi ladite pluralité de fréquences cycliques pour laquelle la cohérence cyclique intégrée calculée est maximale.
  6. Procédé de surveillance selon la revendication 5 dans lequel deux fréquences cycliques consécutives de ladite pluralité de fréquences cycliques, notées respectivement α et α + Δα, sont choisies telles que le rapport α/dα est un nombre entier.
  7. Procédé de surveillance selon l'une quelconque des revendications 1 à 6 dans lequel le nombre entier M est compris entre 6 et 10.
  8. Procédé de surveillance selon l'une quelconque des revendications 1 à 7 comprenant en outre une étape de notification dudit défaut comprenant au moins une information parmi :
    - une identification du roulement défectueux ;
    - une indication d'un élément défectueux sur ledit roulement ; et/ou
    - une indication d'une sévérité du défaut détecté sur le roulement.
  9. Procédé de surveillance selon l'une quelconque des revendications 1 à 8 dans lequel le signal vibratoire a été acquis par le capteur accélérométrique lors d'un régime stationnaire du dispositif tournant.
  10. Procédé de surveillance selon l'une quelconque des revendications 1 à 9 dans lequel la fréquence cyclique est normalisée par rapport à la fréquence de rotation du roulement.
  11. Programme d'ordinateur (PROG) comportant des instructions pour l'exécution des étapes du procédé de surveillance selon l'une quelconque des revendications 1 à 10 lorsque ledit programme est exécuté par un ordinateur.
  12. Support d'enregistrement lisible par un ordinateur sur lequel est enregistré un programme d'ordinateur selon la revendication 11.
  13. Dispositif de surveillance (1) d'un roulement (2) équipant un dispositif tournant (3), ledit dispositif de surveillance comprenant :
    - un premier module d'obtention (1A), configuré pour obtenir un signal vibratoire acquis par un capteur accélérométrique, ledit signal vibratoire comprenant une signature vibratoire du roulement ;
    - un module de traitement (1B) du signal vibratoire, configuré pour éliminer une composante déterministe du signal vibratoire ;
    - un deuxième module d'obtention (1C), configuré pour obtenir pour un défaut déterminé susceptible d'affecter le roulement, une fréquence théorique caractéristique de ce défaut et une déviation maximale déterminée autour de cette fréquence théorique ;
    - un premier module de calcul (1D), configuré pour calculer, en fonction d'une fréquence cyclique, une cohérence cyclique dite intégrée du signal vibratoire traité moyennée sur une bande de fréquences spectrales prédéterminée ;
    - un module d'estimation (1E), configuré pour estimer une fréquence réelle du défaut à partir de la cohérence cyclique intégrée, de la fréquence théorique caractéristique du défaut et de la déviation maximale déterminée autour de cette fréquence théorique ;
    - un deuxième module de calcul (1F) d'un indicateur de diagnostic du défaut, configuré pour calculer ledit indicateur en sommant un nombre entier M de cohérences cycliques intégrées du signal vibratoire évaluées en M fréquences cycliques égales respectivement à M harmoniques de la fréquence réelle estimée du défaut ;
    - un module de comparaison (1G), configuré pour comparer l'indicateur de diagnostic du défaut avec un seuil prédéterminé pour ce défaut et pour détecter le défaut sur le roulement en cas de dépassement du seuil par l'indicateur de diagnostic.
  14. Dispositif de surveillance selon la revendication 13 dans lequel le roulement (2) est un roulement à billes ou à rouleaux et le dispositif tournant est embarqué dans un aéronef.
  15. Moteur d'aéronef (4) comprenant au moins un roulement (2) équipant un dispositif tournant (3) du moteur d'aéronef, au moins un capteur accéléromètrique (5) apte à acquérir un signal vibratoire comprenant une signature vibratoire dudit roulement, et un dispositif de surveillance (1) du roulement selon la revendication 13 ou 14.
EP18842779.3A 2017-12-28 2018-12-19 Procédé et dispositif de surveillance d'un roulement équipant un dispositif tournant Active EP3732457B1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR1763312A FR3076348B1 (fr) 2017-12-28 2017-12-28 Procede et dispositif de surveillance d'un roulement equipant un dispositif tournant
PCT/FR2018/053421 WO2019129956A1 (fr) 2017-12-28 2018-12-19 Procédé et dispositif de surveillance d'un roulement équipant un dispositif tournant

Publications (2)

Publication Number Publication Date
EP3732457A1 EP3732457A1 (fr) 2020-11-04
EP3732457B1 true EP3732457B1 (fr) 2021-09-29

Family

ID=61599430

Family Applications (1)

Application Number Title Priority Date Filing Date
EP18842779.3A Active EP3732457B1 (fr) 2017-12-28 2018-12-19 Procédé et dispositif de surveillance d'un roulement équipant un dispositif tournant

Country Status (4)

Country Link
US (1) US11333576B2 (fr)
EP (1) EP3732457B1 (fr)
FR (1) FR3076348B1 (fr)
WO (1) WO2019129956A1 (fr)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR3079976B1 (fr) * 2018-04-09 2020-04-17 Safran Module d'acquisition pour un systeme de surveillance d'une machine tournante, systeme et procede de surveillance
CN110261116A (zh) * 2019-07-08 2019-09-20 华南理工大学 一种轴承故障检测方法及装置
CN111189640A (zh) * 2020-01-09 2020-05-22 珠海格力电器股份有限公司 一种轴承故障监测方法及采用其的监测装置、洗衣机
FR3119234B1 (fr) * 2021-01-25 2023-10-20 Safran Procédé de surveillance de l’état de composants mécaniques sur une ligne d’arbre, dispositif et système de surveillance associés
CN113029569B (zh) * 2021-03-11 2022-04-29 北京交通大学 一种基于循环强度指数的列车轴承自主故障识别方法
CN113776834B (zh) * 2021-10-11 2022-09-02 山东大学 基于离散余弦循环谱相干的滚动轴承故障诊断方法及系统
WO2022152336A1 (fr) * 2021-10-26 2022-07-21 Západočeská Univerzita V Plzni Procédé de surveillance de vibration d'aube de turbine
CN114636554A (zh) * 2022-03-15 2022-06-17 联合汽车电子有限公司 电驱动系统轴承故障监测方法和装置
FR3134892B1 (fr) 2022-04-22 2024-04-26 Safran Méthode de surveillance d’une machine tournante pour la détection d’un défaut d’un roulement d’aéronefs

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1528377A3 (fr) * 1998-05-20 2006-12-06 DLI Engineering Corporation Appareil de diagnostic avec un enregistreur et un analyseur de données de vibration
FR2803036B1 (fr) * 1999-12-23 2002-10-11 Snecma Detection de l'endommagement de pieces d'un moteur
GB0308467D0 (en) * 2003-04-11 2003-05-21 Rolls Royce Plc Method and system for analysing tachometer and vibration data from an apparatus having one or more rotary components
US8380447B2 (en) * 2010-04-14 2013-02-19 Simmonds Precision Products, Inc. Quantification of condition indicators in the presence of synchronous noise
US20130096848A1 (en) * 2011-10-13 2013-04-18 Charles Terrance Hatch Methods and systems for automatic rolling-element bearing fault detection
FR2994261B1 (fr) * 2012-07-31 2014-07-18 Eurocopter France Procede de detection de defauts d'un roulement par analyse vibratoire

Also Published As

Publication number Publication date
FR3076348A1 (fr) 2019-07-05
US20210063276A1 (en) 2021-03-04
FR3076348B1 (fr) 2020-01-17
WO2019129956A1 (fr) 2019-07-04
EP3732457A1 (fr) 2020-11-04
US11333576B2 (en) 2022-05-17

Similar Documents

Publication Publication Date Title
EP3732457B1 (fr) Procédé et dispositif de surveillance d&#39;un roulement équipant un dispositif tournant
EP2693176B1 (fr) Procédé de détection de défauts d&#39;un roulement par analyse vibratoire
EP1111364B1 (fr) Détection de l&#39;endommagement de pièces d&#39;un moteur
EP2344728B1 (fr) Procede de surveillance d&#39;un turboreacteur
EP3250974B1 (fr) Procédé, système et programme d&#39;ordinateur pour phase d&#39;apprentissage d&#39;une analyse acoustique ou vibratoire d&#39;une machine
EP3658881B1 (fr) Procede et dispositif de recherche d&#39;un defaut susceptible d&#39;affecter un dispositif mecanique tournant de transmission de puissance
EP2507598B1 (fr) Procédé et dispositif de surveillance de vibrations en torsion d&#39;un arbre rotatif d&#39;une turbomachine
FR3079977A1 (fr) Procede et dispositif de surveillance d&#39;une machine tournante
EP3785006B1 (fr) Procédé et dispositif de surveillance d&#39;un système à engrenages
EP3531098B1 (fr) Procede de surveillance et de detection de la formation d&#39;une degradation dans au moins une piece mobile d&#39;un mecanisme tournant et systeme associe
EP4384792A1 (fr) Procédé de détection d&#39;un défaut de roulement d&#39;un système rotatif et système de surveillance mettant en oeuvre ce procédé
EP3136111B1 (fr) Procédé pour mesurer la vitesse de rotation d&#39;une roue de vehicule
WO2021089936A1 (fr) Procédé, dispositif et programme d&#39;ordinateur de surveillance d&#39;une machine tournante d&#39;un aéronef
WO2023203304A1 (fr) Méthode de surveillance d&#39;une machine tournante pour la détection d&#39;un défaut d&#39;un roulement d&#39;aéronefs
FR3087265A1 (fr) Procede et dispositif de surveillance d&#39;une piece mecanique en rotation
WO2023187285A1 (fr) Méthode de surveillance du régime moteur d&#39;un aéronef
WO2022157453A1 (fr) Procede de surveillance de l&#39;etat de composants mecaniques sur une ligne d&#39;arbre, dispositif et systeme de surveillance associes
FR3132766A1 (fr) Estimation d’une amplitude de flottement d’une soufflante d’une turbomachine
WO2023073314A1 (fr) Procede de surveillance d&#39;un systeme mecanique rotatif d&#39;un aeronef
FR3140946A1 (fr) Procédé de détection d’un mode de fonctionnement d’une machine tournante, notamment pour un aéronef en cours de vol

Legal Events

Date Code Title Description
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: UNKNOWN

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE

PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20200616

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

DAV Request for validation of the european patent (deleted)
DAX Request for extension of the european patent (deleted)
GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: GRANT OF PATENT IS INTENDED

INTG Intention to grant announced

Effective date: 20210521

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE PATENT HAS BEEN GRANTED

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

Free format text: NOT ENGLISH

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

Ref country code: AT

Ref legal event code: REF

Ref document number: 1434602

Country of ref document: AT

Kind code of ref document: T

Effective date: 20211015

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602018024430

Country of ref document: DE

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

Free format text: LANGUAGE OF EP DOCUMENT: FRENCH

REG Reference to a national code

Ref country code: LT

Ref legal event code: MG9D

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: FI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

Ref country code: RS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

Ref country code: SE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

Ref country code: NO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: HR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211229

Ref country code: LT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

REG Reference to a national code

Ref country code: NL

Ref legal event code: MP

Effective date: 20210929

REG Reference to a national code

Ref country code: AT

Ref legal event code: MK05

Ref document number: 1434602

Country of ref document: AT

Kind code of ref document: T

Effective date: 20210929

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LV

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20211230

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: AT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220129

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

Ref country code: RO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220131

Ref country code: PL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

Ref country code: NL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

Ref country code: ES

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

Ref country code: CZ

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

Ref country code: AL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

REG Reference to a national code

Ref country code: DE

Ref legal event code: R097

Ref document number: 602018024430

Country of ref document: DE

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MC

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

Ref country code: DK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

26N No opposition filed

Effective date: 20220630

REG Reference to a national code

Ref country code: BE

Ref legal event code: MM

Effective date: 20211231

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LU

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20211219

Ref country code: IE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20211219

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

Ref country code: BE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20211231

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LI

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20211231

Ref country code: CH

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20211231

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CY

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SM

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929

Ref country code: HU

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT; INVALID AB INITIO

Effective date: 20181219

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 20231124

Year of fee payment: 6

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: FR

Payment date: 20231122

Year of fee payment: 6

Ref country code: DE

Payment date: 20231121

Year of fee payment: 6

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20210929