METHOD AND DEVICE FOR ESTIMATING THE STATE OF WEAR OF A JOURNAL BEARING
20200096418 ยท 2020-03-26
Inventors
- Sebastian Nowoisky (Michendorf, DE)
- Noushin Mokhtari Molk Abadi (Berlin, DE)
- Mateusz GRZESZKOWSKI (Berlin, DE)
Cpc classification
F16C2233/00
MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
G01N2291/0258
PHYSICS
F16C2361/65
MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
G01N29/40
PHYSICS
F16C17/246
MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
G01N29/449
PHYSICS
International classification
G01L1/25
PHYSICS
F16C17/24
MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
G01L1/10
PHYSICS
Abstract
A method for estimating the state of wear of a plain bearing comprises: establishing a time profile of at least one friction event from a structure borne noise signal by a mathematical friction event model, determination of a measure, which characterizes at least one friction event based on a time duration of the at least one friction event, the amplitude of the at least one friction event and/or an integral measure over the at least one friction event, combination of the measure, accumulated over time, with an angle indication for the at least one friction event at the circumference of the plain bearing in order to determine a spatially resolved wear model of the plain bearing, wherein the spatially resolved wear model is obtained by an estimating filter, and outputting a signal in accordance with the wear model to characterize the state of wear.
Claims
1. A method for estimating the state of wear of a plain bearing having a shaft mounted therein, in particular rotating therein, wherein at least one time-dependent structure borne noise signal of the plain bearing is recorded by at least one structure borne noise sensor, in particular precisely one structure borne noise sensor, comprising a) establishing the time profile of at least one friction event from the structure borne noise signal by means of a mathematical friction event model, b) determination of a measure, which characterizes at least one friction event based on the time duration of the at least one friction event, the amplitude of the at least one friction event and/or an integral measure over the at least one friction event, c) combination of the measure, accumulated over time, with an angle indication for the at least one friction event at the circumference of the plain bearing in order to determine a spatially resolved wear model of the plain bearing, wherein the spatially resolved wear model (W) is obtained by means of an estimating filter, and d) outputting of a signal in accordance with the wear model to characterize the state of wear.
2. The method according to claim 1, wherein the spatially resolved wear model is obtained by means of a Kalman filter, an extended Kalman filter and/or a regression.
3. The method according to claim 1, wherein the number of friction events is counted and used to generate a signal for service life estimation.
4. The method according to claim 1, wherein the friction events comprise mixed friction and/or solid body friction events.
5. The method according to claim 1, wherein a correlation of the frequencies of the friction events and/or between the feature determined (from the structure borne noise) and the degree of wear is carried out.
6. The method according to claim 1, wherein current operating data, in particular the operating duration, a temperature, in particular an oil temperature, load data, acoustic data in frequency ranges outside the friction events and/or a measured rotational speed are included in the service life estimation, and/or operating data of structurally identical and/or structurally similar plain bearings are included in the service life estimation and/or priori information is included in the service life estimation.
7. (canceled)
8. (canceled)
9. The method according to claim 1, wherein the friction event model is an envelope curve model which is obtained by means of the following steps: a) filtering the structure borne noise signal, b) subsequent calculation of an envelope curve for the filtered structure borne noise signal, c) subsequent smoothing of the envelope curve.
10. The method according to claim 9, wherein the angle indication for the at least one friction event is obtained from the envelope curve model, in which a) combination of the data for the smoothed envelope curve with a rotation angle signal dependent on the rotation of the shaft in the plain bearing is carried out, b) calculation of at least one maximum is carried out, which correlates with the at least one friction event, from the combined data from step a) for the determination of an angle indication for the at least one friction event at the circumference of the plain bearing.
11. The method according to claim 10, wherein the rotation angle signal is determined and/or generated by pattern recognition or by means of a reference pulse by an incremental encoder, in particular a magnetic reference pulse.
12. The method according to claim 10, wherein the rotation angle signal is generated exclusively by the movement of the shaft and/or of the plain bearing, in particular by at least one magnetic element of the shaft and/or in the plain bearing and a correspondingly associated magnetic sensor.
13. The method according to claim 10, wherein the rotation angle signal is generated actively by means of at least one pulse, in particular a zero pulse or a multiplicity of pulses of the incremental encoder.
14. The method according to claim 9, wherein the filtering system has a high pass filter, in particular with a cutoff frequency between 50 and 300 kHz, in particular between 80 and 150 kHz.
15. The method according to claim 9, wherein the calculation of the envelope curve is performed by means of a Hilbert transformation or by averaging over a predetermined quantity of filtered structure borne noise data points.
16. The method according to claim 9, wherein the envelope curve is smoothed by means of a smoothing filter, in particular a Savitzky-Golay filter.
17. The method according to claim 1, wherein the plain bearing is arranged in a planetary gearbox, in particular a planetary gearbox of a wind turbine, a vehicle or an aircraft engine.
18. The method according to claim 1, wherein kinematic motion data and/or structure borne noise events of the planetary gearbox, in particular the motion data and/or structure borne noise events of the movements of the sun gear, planet carrier and/or planet gears, are filtered out.
19. A device for estimating the state of wear of a plain bearing having a shaft mounted therein, in particular rotating therein, with respect to at least one friction event, wherein at least one time-dependent structure borne noise signal can be recorded by at least one structure borne noise sensor, in particular precisely one structure borne noise sensor of the plain bearing, having a first computation means for establishing the time profile of at least one friction event from the structure borne noise signal by means of a mathematical friction event model, a second computation means for determining a measure, which characterizes the at least one friction event based on the time duration of the at least one friction event, the amplitude of the at least one friction event and/or an integral measure over the at least one friction event, a third computation means for combining the measure, accumulated over time, with an angle indication for the at least one friction event at the circumference of the plain bearing in order to determine a spatially resolved wear model of the wear model of the plain bearing (1), wherein the spatially resolved wear model is obtained by means of an estimating filter, and a signaling means for outputting a signal in accordance with the model to characterize the state of wear.
20. The device according to claim 19, wherein the at least one structure borne noise sensor is arranged on the end of a holder of the plain bearing.
21. The device according to claim 19, wherein the at least one structure borne noise sensor has a piezoelectric element for recording the structure borne noise.
22. The device according to claim 19, wherein the at least one structure borne noise sensor is arranged in the immediate vicinity of the circumference of the plain bearing, in particular in the immediate vicinity of the introduction of a force.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
[0047] The proposed solution will be discussed in connection with the exemplary embodiments illustrated in the figures.
[0048]
[0049]
[0050]
[0051]
[0052]
[0053]
[0054]
[0055]
[0056]
[0057]
[0058]
[0059]
[0060]
[0061]
[0062]
DETAILED DESCRIPTION
[0063] The monitoring of hydrodynamic plain bearings for mixed friction events is described below with reference to a number of illustrative embodiments.
[0064] Theoretically, a hydrodynamic plain bearing 1 has an infinite service life as long as the shaft 6 and the lining of the plain bearing 1 are separated from one another by a loadbearing lubricating film. As soon as these two components come into contact, there is mechanical friction (mixed friction, solid body friction), which ultimately leads to damage. As a result, the plain bearing 1 can lose its ability to function since it is no longer possible for a loadbearing lubricating film to form in the presence of increased abrasion or damage.
[0065] Embodiments relating to the estimation of the state of wear are described below, proceeding in one step on the basis that the time profile of a friction event a, b, c, d is determined from a structure borne noise signal S by means of a mathematical friction event model R. First of all, therefore
[0066] One known method for monitoring hydrodynamic plain bearings 1 uses shaft orbit plots (see J. Deckers, Entwicklung einer Low-Cost Krperschallsensorik zur berwachung des Verschleiverhaltens von wlz-oder gleitgelagerten Kreiselpumpen kleiner Leistung, [Development of a Low-Cost Structure Borne Noise Sensor System for Monitoring the Wear Behavior of Low-Power Centrifugal Pumps with Rolling or Plain Bearings] Dissertation, Gerhard Mercator-Universitat Duisburg, Duisburg, 2001).
[0067] Here, the shaft orbit within the plain bearing is detected by two position sensors mounted orthogonally on the plain bearing 1. The two phase-shifted position signals detected in this case are represented in polar coordinate form in orbit plots. These plots represent the rotation angle-dependent movement of the supported shaft 6 transversely to the axial shaft axis.
[0068] To detect the phase position, a key phasor (reference transducer) is used.
[0069] If the shaft 6 then leaves the permitted orbit, a rubbing contact (i.e. a mixed friction event) between the shaft 6 and the lining of the plain bearing 1 has taken place. This can be identified in the shaft orbit plot.
[0070] By means of this method, it is not only possible to identify a rubbing contact process but also to ascertain the intensity and position of the contact in the circumferential direction of the plain bearing lining. Mixed friction events a, b, c, d are discussed below but the technical teaching is not restricted to this type of friction analysis.
[0071] The mutual contact between roughness peaks on the two sliding partners in a mixed friction event (rubbing contact process) causes structure borne noise with a frequency of up to 2 MHz in the plain bearing 1.
[0072] As compared with other diagnostic methods, the use of structure borne noise analysis offers advantages in respect of early detection of bearing damage to the plain bearing 1 (see M. Fritz, A. Burger and A. Albers, Schadensfrherkennung an geschmierten Gleitkontakten mittels Schallemissionsanalyse, [Early Detection of Damage to Lubricated Sliding Contacts by Means of Noise Emission Analysis], Institut fr Maschinenkonstruktionslehre and Kraftfahrzeugbau, Bericht, Universitt Karlsruhe, 2001; P. Raharjo, An Investigation of Surface Vibration, Airborne Sound and Acoustic Emission Characteristics of a Journal Bearing for Early Fault Detection and Diagnosis, Dissertation, University of Huddersfield, May 2013)
[0073] By means of suitable signal processing and feature extraction algorithms, it is possible to distinguish viscous friction, which does not affect the service life, from mixed and solid body friction, which reduce the service life. However, the algorithms used for diagnosis by means of structure borne noise assess the state of friction only globally, not locally, over the circumference of the plain bearing 1, i.e. there is no angular resolution in the identification of the mixed friction events a, b, c, d.
[0074] But a knowledge of local mixed friction processes a, b, c, d is essential for characterization of the state of wear of plain bearings 1. Repeated friction at the position =20 for example (plotted in
[0075] One phenomenon which occurs with the superposition of a high-frequency carrier signal and a low-frequency useful signal is amplitude modulation.
[0076] In the case of local contact between the shaft 6 and the lining of the plain bearing 1, there is likewise amplitude modulation of the structure borne noise signal (see M. Leahy, D. Mba, P. Cooper, A. Montgomery and D. Owen, Experimental investigation into the capabilities of acoustic emission for the detection of shaft-to-seal rubbing in large power generation turbines, Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology, Vol. 220, No. 7, p. 607-615, 2006; A. Albers and M. Dickerhof, Simultaneous Monitoring of Rolling-Element and Journal Bearings Using Analysis of Structure-Born Ultrasound Acoustic Emissions, in International Mechanical Engineering Congress & Exposition, Vancouver, British Columbia, Canada, 2010.)
[0077] The mixed friction events a, b, c, d occur in dependence on the rotational speed of the shaft 6. The mixed friction events a, b, c, d themselves each generate a structure borne noise signal, which is of significantly higher frequency than the rotation frequency of the shaft 6. Overall, a structure borne noise signal S in which a low-frequency rotation frequency and a higher-frequency structure borne noise signal are superposed is recorded.
[0078] The diagrammatic signal profile of the structure borne noise signal S can be seen in
[0079] In the case of viscous friction, with which no rubbing contact processes occur between the shaft 6 and the lining of the plain bearing 1, no amplitude modulations occur.
[0080] As mentioned above, the exact circumferential position (i.e. the angle ) at which mixed friction events a, b, c, d occur represents important information, including for service life prediction. Thus, the accumulation of mixed friction events a, b, c, d at one circumferential position can be interpreted as a measure of the wear of the plain bearing lining. Service life prediction can also be improved. It is also possible, for example, to form the integral over one revolution in each case and this takes account of the intensity and duration of the friction event at the same time.
[0081] The intention is, with appropriate monitoring, to keep both the complexity of the measurement chain and the costs for the production of a product (in this case the diagnostic or prediction system) as low as possible. Reducing sensor numbers, ideally the use of just one sensor, simplifies the measurement chain and also allows a significant cost reduction.
[0082] Embodiments which use various properties of the detected structure borne noise signal S to detect mixed friction events a, b, c, d are described below.
[0083] An envelope curve, also referred to as an envelope, envelops a family of curves (e.g. that of a structure borne noise signal S as per
[0084] For the application presented here, however, it is fully sufficient to know the behavior of the amplitudes for the determination of the envelope. All that is required is to ascertain where local maxima and minima occur. The envelope curve is obtained by determining the RMS (root mean square or quadratic mean) over a predetermined quantity of data points. Alternatively, it is also possible to use a Hilbert transformation. This also provides a measure of the energy of the structure borne noise signal, wherein the peak values and curve shape are taken into account in each case.
[0085] The curve formed by an envelope has bends and/or sharp angles which need to be smoothed to allow a judgment about local maxima and minima. For this purpose, it is possible, for example, to use low-order approximation polynomials to achieve the best possible smoothing. One possibility for smoothing is to use the Savitzky-Golay filter (see A. Savitzky and M. J. E. Golay, Smoothing and Differentiation of Data by Simplified Least Squares Procedures, Anal. Chem., July 1964).
[0086] This method smooths a signal by section-by-section fitting of a polynomial function to the signal. This fitting process employs the method of least squares between the matrix X and the vector y:
y=X b
[0087] The solution for b with the aid of the least squares is
b=(X.sup.TX).sup.1X.sup.Ty.
[0088] The estimated values used for smoothing are:
=Xb=X(X.sup.TX).sup.1X.sup.Ty=Hy
[0089] First of all, embodiments of the method for mixed friction localization over the circumference of the plain bearing lining by means of structure borne noise measurement are explained in more detail below.
[0090] In
[0091] In a first step 201, the amplitude-modulated structure borne noise signal S is subjected to high-pass filtering in order as far as possible to attenuate disturbance signals from the surroundings and signals which have nothing to do with the mixed friction events a, b, c, d. In the present case, a cutoff frequency of 100 kHz is used. In other embodiments, it is also possible to use other cutoff frequencies. Frequencies in the range of 50 to 300 kHz are generally appropriate since the mixed friction events are very largely above these cutoff frequencies.
[0092] In the following step 202, the envelope of the structure borne noise signal S is formed, e.g. by means of Hilbert transformation. The method used here forms the average over a certain number B of signal points n, (e.g. 800) and stores this value in a vector.
[0093] As explained above, it is possible when using the RMS to determine the energy of the envelope, but this has sharp bends or spikes. To enable the local maxima and minima to be ascertained with greater numerical accuracy, a 3rd-order Savitzky-Golay filter can be used for example to smooth the structure borne noise signal S (step 203). In other embodiments, it is also possible to use other filters.
[0094] Using a first computation means, a time profile was thereby obtained by means of a mathematical friction event model R.
[0095] In step 204, a measure M, which characterizes the at least one mixed friction event a, b, c, d based on the time duration T.sub.a, T.sub.b, T.sub.c, T.sub.d of the mixed friction event a, b, c, d, the amplitude A.sub.a, A.sub.b, A.sub.c, A.sub.d of one friction event a, b, c, d and/or an integral measure over the friction event a, b, c, d is then determined, using a second computation means.
[0096] The measures are explained in greater detail in conjunction with
[0097] The three periodically occurring mixed friction events a, b, c each have different lengths T.sub.a, T.sub.b, T.sub.c, i.e. they are a measure of the friction contact time over the circumference of the plain bearing 1. The lengths of time are thus a measure of the magnitude of the contact angle. The smaller the contact angle, the larger is the kurtosis of the signal at the same intensity. Kurtosis describes the steepness or sharpness of the signal. As the contact angle becomes larger, kurtosis decreases. Thus, kurtosis can be used as a measure M of the duration of the mixed friction events.
[0098] In addition or as an alternative, it is also possible to use the intensity of the mixed friction events a, b, c as a measure M for the mixed friction events a, b, c. The intensity can for example be determined by way of the amplitudes A.sub.a, A.sub.b, A.sub.c (see
[0099] T.sub.a, T.sub.b, T.sub.c, with the result that a combined measure M is obtained which allows equally for the duration and the amplitude. It is also possible to determine the integral by means of a squared area (with or without time weighting), where the larger amplitudes (and, in the case of time weighting, events of longer duration) are given greater consideration.
[0100] In all cases, measures which, whether weighted or unweighted, characterize the duration and intensity of the mixed friction events are obtained.
[0101] The measure M, accumulated over time, is then combined in step 205 with an angular position (t) for the mixed friction events a, b, c, d at the circumference of the plain bearing in order to determine a spatially resolved wear model W of the wear model of the plain bearing 1.
[0102] The way in which the angular position (t) can be determined in one embodiment is described below.
[0103] In the illustrated embodiment, an incremental encoder is used to output pulses to the plain bearing 1. These pulses are rotation angle signals Z, which are dependent on the rotation of the shaft 6 in the plain bearing 1.
[0104] In this application, the zero pulse signal (Z signal) of the incremental encoder is used to identify the exact angular position (t) of the mixed friction events a, b, c, d. Precisely one revolution of the shaft 6 in the plain bearing 1 takes place between two square wave signals of an incremental encoder. Both signals, both the structure borne noise signal S and the Z signal Z are recoded simultaneously (step 205). For improved accuracy, it is also possible to use more than one pulse signal per revolution.
[0105] When the processed structure borne noise signal S and the Z signal Z (step 205) are superposed, it is possible to make an accurate association between the maximaresulting from the structure borne noise of the mixed friction eventsand the angular position (step 206), wherein a signal D is output in accordance with the model wear model W, and this is transmitted to an engine control system, for example.
[0106] The signals can be processed and evaluated by means of a computer 30. The signals of the structure borne noise sensor 3 can be transmitted to the computer 30 in a conventional manner, if appropriate via an amplifier.
[0107] Here, each maximum represents a rubbing contact between the shaft and the plain bearing lining, i.e. a mixed friction event a, b, c. Given a knowledge of the signal relating to the measured angular positions, each maximum can be associated with an angle which characterizes the point of rubbing contact (i.e. the mixed friction event a, b, c) in the lining of the plain bearing 1.
[0108]
[0109] In step 203, the individual mixed friction events a, b, c are identified. In a counter (step 207) these are counted individually, and a sum is formed. Here, it is possible to ascertain, with the aid of features (RMS, envelope curve, kurtosis) and pattern recognition, whether there is mixed friction (or solid body friction). A distinction is not necessarily drawn between location, duration and intensity. Mixed friction events which occur are then counted up until a critical value is reached. After this, a warning signal is output.
[0110] Together with the result of step 206 and/or of signal D, the remaining service life is estimated (step 208). This process can also incorporate current operating data B.
[0111] However, it is also possible for the count E to be output and displayed separately.
[0112] If the method is used to estimate the state of wear in an engine, for example, the service life estimation could also include temperatures, such as the oil temperature in a planetary gearbox, load data, acoustic data in frequency ranges outside the friction events and/or a rotational speed, for example, as current operating data. If, for example, the engine had been operated very gently throughout, e.g. with few takeoffs and landings, a relatively long service life would be displayed. If there were a change in the operating behavior, e.g. more short distances were flown, this would be reflected by a shorter estimate for the service life.
[0113] Results in which both mixed friction andfor comparisonviscous friction occurred are illustrated below. For this purpose, use was made of embodiments which are illustrated in
[0114] In
[0115] In the illustrated embodiment, just one structure borne noise sensor 3 is required and, in the illustrated embodiment, is arranged offset slightly sideways from the center, in the vicinity of the circumference of the plain bearing 1 and on the front side of the plain bearing device 10. The structure borne noise generated during the operation of the plain bearing 1 is transmitted well to the structure borne noise sensor 3 by the solid bodies. It is worthwhile here to arrange the structure borne noise sensor 3 in the vicinity of the introduction of an external force. A Physical Acoustics WD 100-900 kHz wideband sensor can be used as the structure borne noise sensor 3, for example. The structure borne noise sensor 3 can have a piezoelectric element.
[0116] In the illustrated embodiment, a force F.sub.N (i.e. a bearing load) is imposed from above (see
[0117]
[0118] It is expected that there is no modulation in the signal obtained when the plain bearing 1 is operated with viscous friction since no rubbing contact occurs between the shaft 6 and the lining of the plain bearing 1. In the case of mixed friction events, it should be possible to detect a modulation in the signal.
[0119] Under a constant load F.sub.N, a falling rotational speed ramp was run. Each rotational speed was held for three seconds. This makes it possible, at a constant load F.sub.N, to move from viscous friction, which occurs at high rotational speeds, into the range of mixed friction.
[0120] The structure borne noise signal S and the Z signal Z of the incremental encoder are illustrated by way of example in
[0121] The Z signal Z is the signal of the incremental encoder and is output once per revolution as a square wave signal. Precisely one revolution of the shaft 6 takes place between two square wave signals.
[0122]
[0123] As is apparent, mixed friction events a, b, c, d take place between the shaft 6 and the lining of the plain bearing 1. The modulation in the signal can be observed in
[0124] In
[0125] This structure borne noise signal S is subsequently processed using one embodiment of the method, as illustrated, for example, in conjunction with
[0126] The time is plotted on the x axis of
[0127]
[0128] Since the maxima are to be determined numerically, the envelope of the structure borne noise signal S is smoothed using the Savitzky-Golay filter (
[0129] The signal that is then formed can be seen in
[0130] From the above descriptions, it is clear that a plain bearing 1 can be efficiently monitored for rubbing contacts (i.e. mixed friction events a, b, c, d) by arranging a structure borne noise sensor 3 in the vicinity of the plain bearing 1. Together with the pulse generator and a computer for evaluating the data, it is possible in this way to efficiently monitor a plain bearing 1 in a motor or an aircraft engine, for example.
[0131] The monitoring of plain bearings 1 in an epicyclic planetary gearbox 20 is described below as one possible application.
[0132]
[0133] The planetary gearbox 20 can be driven via the sun gear 23, which rotates at an angular speed .sub.S. The planet gears 22 roll on the sun gear 23 and in the ring gear 21, which is assumed to be fixed here. The shafts 6 of the planet gears 22 are supported on the carrier 24 by means of plain bearings 1, with the result that the planet gears 22 rotate at an angular speed .sub.P. The carrier 24, which forms the output of the planetary gearbox 20 in the illustrated embodiment, rotates around the axis of the sun gear 23 at the angular speed .sub.C.
[0134] In the embodiment shown in
[0135] As an alternative, it is also possible for the structure borne noise sensor 3 to be arranged on the corotating carrier 24, as illustrated in
[0136] Here, the illustration of the epicyclic planetary gearbox 20 should be taken to be only illustrative. In other embodiments, it is possible to use five or more planet gears, for example. It is also possible to choose a different mechanism, that is to say that the input and output differ from the example in
LIST OF REFERENCE SIGNS
[0137] 1 Plain bearing [0138] 2 Holder of plain bearing [0139] 3 Structure borne noise sensor [0140] 4 Supporting hole for introducing force [0141] 5 Oil supply line [0142] 6 Shaft [0143] 7 First support bearing [0144] 8 Second support bearing [0145] 9 Electric motor [0146] 10 Plain bearing device [0147] 20 Planetary gearbox [0148] 21 Ring gear [0149] 22 Planet gears [0150] 23 Sun gear [0151] 24 Carrier [0152] 30 Computer [0153] a, b, c, d Rubbing contact points, friction events [0154] A.sub.a, A.sub.b, A.sub.c, A.sub.d Amplitude of a friction event [0155] B Current operating data [0156] D Output signal [0157] E Count [0158] F.sub.N Force introduced [0159] I.sub.a, I.sub.b, I.sub.c Integral measure over one friction event [0160] M Measure of friction event [0161] R Friction event model [0162] S Structure borne noise signal [0163] T.sub.a, T.sub.b, T.sub.c, T.sub.d Time duration of a friction event [0164] W Wear model [0165] Z Signal of pulse generator