Maintenance Management Based on Signal Processing

1. Wind Turbines

Most of the wind turbines are three-blade units (Figure 1.) [55]. Once the wind drives the blades, the energy is transmitted via the main shaft through the gearbox (supported by the bearings) to the generator. The generator speed must be as near as possible to the optimal for the generation of electricity. At the top of the tower, assembled on a base or foundation, the housing or nacelle is mounted and the alignment with the direction of the wind is controlled by a yaw system. There is also a pitch system in each blade. This mechanism controls the wind power and sometimes is employed as an aerodynamic brake. The wind turbine features a hydraulic brake to stop itself when it is needed. Finally, there is a meteorological unit that provides information about the wind (speed and direction) to the control system.

1.1. Maintenance in Wind Turbines

Maintenance is a key tool to ensure the operation of all components of a set. One of the objectives is to use available resources efficiently. The classical theory of maintenance was focused on the corrective and preventive maintenance [9] but alternatives to corrective and preventive maintenance have appeared in recent years. One of them is Condition Based Maintenance, which ensures the continuous monitoring and inspection of the wind turbine detecting emerging faults and organizing maintenance tasks that anticipate the failure [59]. Condition Based Maintenance implies acquisition, processing, analysis and interpretation of data and the selection of proper maintenance actions. This is achieved using condition monitoring systems [27, 28]. Thereby, CBM is presented as a useful technique to improve not only the maintenance but the safety of the equipments. Byon and Ding [14] or McMillan and Ault [50] have demonstrated its successful application in wind turbines, making the CBM one of the most employed strategies in this industry. Another example of the maintenance evolution is the Reliability Centred Maintenance. It is defined as a process to determine what must be done to ensure that any physical asset works in its operating context [71]. Nowadays it is the most common type of maintenance for many industrial fields [25, 26] and it involves maintenance system functions or identifying failure modes among others maintenance tasks [52].

Figure 1.

Main parts of a turbine: (1) blades, (2) rotor, (3) gearbox, (4) generator, (5) bearings, (6) yaw system and (7) tower [36].

1.2. Condition Monitoring applied to Wind Turbines

Condition Monitoring systems operate from different types of sensors and signal processing equipments. They are capable of monitoring components ranging from blades, gearboxes, generators to bearings or towers. Monitoring can be processed in real time or in packages of time intervals. The procurement of data will be critical to determine the occurrence of a problem and determine a solution to apply. Therefore, the success of a Condition Monitoring system will be supported by the number and type of sensors used and the signal collection and processing.

Any element that performs a rotation is susceptible of being analysed by vibration. In the case of the wind turbines, vibration analysis is mainly specialized in the study of gearboxes [48, 49] and bearings [81] [85]. Different types of sensors will be required depending on the operating frequency: position transducers, velocity sensors, accelerometers or spectral energy emitted sensors.

Acoustic emissions (AE) describe the sound waves produced when a material undergoes stress as a result of an external force [35]. They can detect the occurrence of cracks in bearings [84] and blades [91] in earlier stages.

Ultrasonic tests evaluate the structural surface of towers and blades in wind turbines [22] [24]. Consistent with some other techniques, it is capable of locating faults safely.

Oil analysis may determine the occurrence of problems in early stages of deterioration. It is usually a clear indicator of the wearing of certain components. The technique is widely used in the field of maintenance, being important for gearboxes in wind turbines [47].

Thermographic technique is established for monitoring mainly electrical components [72]; although its use is extended to the search of abnormal temperatures on the surfaces of the blades [64]. Using thermography, hot spots can be found due to bad contacts or a system failure. It is common the introduction of online monitoring systems based on the infrared spectrum.

There are techniques that not being so extended, are also used in the maintenance of wind turbines. In many cases, their performance is heavily influenced by the costs or their excessive specialization, making them not always feasible. Some examples are strain measurements in blades [68]; voltage and current analysis in engines, generators and accumulators [67]; shock pulse methods detecting mechanical shocks for bearings [13] or radiographic inspections to observe the structural conditions of the [61].

1.3. Signal processing methods

Fast Fourier Transform (FFT)

The FFT converts a signal from the time domain to the frequency domain. The use of FFT also allows its spectral representation [56]. Each frequency range is framed into a particular failure state. It is very useful when periodic patterns are searched [5]. Vibration analysis also provides information about a particular reason of the fault origin and/or its severity [43]. There is extensive literature demonstrating the development of the method for rolling elements. The FFT of a function f(x) is defined as [12]:

$${\displaystyle {\int }_{-\infty }^{\infty }f(x)e^{-2i\pi xs}dx}$$ E1

This integral, which is a function of s, may be written as F(s). Transforming F(s) by the same formula, equation (2) where F(s) is the Fourier transform of f(x) is obtained.

$${\displaystyle {\int }_{-\infty }^{\infty }f(s)e^{-2i\pi \omega s}ds}$$ E2

There are a considerable number of publications regarding the diagnosis of faults for rolling machinery that justifies the models and patterns based on the Fast Fourier Transform. Misalignment is one of the most commonly observed faults in rotating machines, being the second most common malfunction after unbalance. It may be present because of improper machine assembly, thermal distortion and asymmetry in the applied load. Misalignment causes reaction forces in couplings that are the major cause of machinery vibration. Some authors evaluated numerically the effect of coupling misalignment and suggested the occurrence of strong vibrations at twice the natural frequency [70] [95], although rotating machinery can excite vibration harmonics from twice to ten harmonics depending on the signal pickup locations and directions [53].

Faults do not have a unique nature and most of the time, problems on a smaller scale are linked, e.g. in the case of misalignment, when an angular misalignment is studied, parallel misalignment (minor fault) needs to be take into account. Al-Hussain and Redmond reported vibrations for parallel misalignment at the natural frequency from experimental investigations [4].

To facilitate the diagnosis in rolling elements, some companies and researchers tabulate the most common failure modes in the frequency domain, so that the analysis can be carried out easier. Thus, the appearance of different frequency peaks determines the existence of developing problems such as gaps, unbalances or misalignments among other circumstances [31].The great advantage of these tables is that the value of the frequency peak is not a particular value and may be adapted to any situation where the natural frequency (or the rotational speed) is known.

Wavelet transform is a time-frequency technique similar to Short Time Fourier Transform although it is more effective when the signal is not stationary. Wavelet transform decompose an input signal into a set of levels at different frequencies [77]. Wavelet transforms have been applied to the fault detection and diagnosis in various wind turbine parts.

A hidden Markov model is a statistical model in which the system being modelled is assumed to be a Markov process with hidden states. A hidden Markov model can be considered as the simplest dynamic Bayesian network [8]. Ocak and Loparo presented the application for the bearing fault detection [57].

They are used when a statistical study is required. In these cases, common statistical, i.e. the root mean square or peak amplitude; to diagnose faults are employed. Other parameters can be maximum or minimum values, means, standard deviations to energy ratios or kurtosis. Moreover, trend analysis refers to the collection of information in order to find a trend.

There are many methods that, as happened with the techniques available for CM, are very specific and therefore they are used for very specific situations. Filtering methods, for example, are designed to remove any redundant information, eliminating unnecessary overloads in the process. Analysis in time domain will be a way of monitoring wind turbine faults as inductive imbalances o turn-to-turn faults. Other methodology, the power cepstrum, defined as the inverse Fourier Transform of the logarithmic power spectrum [92], reports the occurrence of deterioration through the study of the sidebands. Time synchronous averaging, amplitude demodulation and order analysis are other signal processing methodologies used in wind turbines.

2. Wavelet transform

The wavelet transform is a method of analysis capable of identifying the local characteristics of a signal in the time and frequency domain. It is suitable for large time intervals, where great accuracy is requested at low frequencies and vice versa, e.g. small regions where precision details for a deeper processing are required at higher frequencies [23]. The wavelet transform can be defined as a signal on a temporal base that is filtered successive times and whose average value is zero. These wavelets are irregular and asymmetrical [51]. The transform has many applications in control process and detection of anomalies. It enables to analyse the signal structures that depend on time and scale, being a useful method to characterize and identify signals with spectral features, unusual temporary files and other properties related to the lack of stationary. When the frequency range corresponding to each signal is known, the data can be studied in terms of time, frequency and amplitude. Therefore it is possible to see which frequencies are in each time interval, and may even reverse the wavelet transform when it is necessary. Previously to the wavelet transform, the FFT was able to work with this type of signals in the frequency domain but without great resolution in the time domain [38].

The wavelet transform of a function f(t) is the decomposition of f(t) in a set of functions and ψ _s,τ (t), forming a base. It is defined as [88] [66]:

$$W_f(s,\tau )={\displaystyle \int f(t){\psi }_{s,\tau }^{*}(t)dt}$$ E3

Wavelets transforms are generated from the translation and scale change from a same wavelet function ψ(t), called mother wavelet, which is given by equation (4):

$${\psi }_{s,\tau }(t)=\frac{1}{\sqrt{s}}\psi \left(\frac{t-\tau }{s}\right)$$ E4

where s is the scale factor, and τ is the translational factor.

The wavelets ψ _s,τ (t) generated from the same mother wavelet function ψ(t) have different scale s and location τ, but the same shape. Scale factors are always s>0. The wavelets are dilated when the scale s>1 and contracted when s<1. Thus, the changing of the value s can cover different ranges of frequencies. Large values for the parameter s correspond to lower frequencies ranges or a large scale for ψ _s,τ (t). Small values of s correspond to lower frequencies ranges or very small scales.

The wavelet transform can be continuous or discrete. The difference between them is that the continuous transform provides more detailed information but consuming more computation time while the discrete signal is efficient with fewer parameters and less computation time [17]. The Discrete Wavelet Transform coefficients are a group of discrete intervals of time and scales. These coefficients are used to formalize a set of features that characterize different types of signals. Any signal can be divided into low frequency approximations (A) and high frequency details (D). The sum of A and D is always equal to the original signal. The division is done using filters (Figure 2).

Figure 2.

Decomposition diagram.

To reduce the computational and mathematical costs due to duplication of data, a sub-sampling is usually performed, containing the half of the collected information from A and D but without losing information. It is common to accompany this information with a graphical representation where the original signal is divided in low pass filters and high pass filters [15]. When the signals are complex, the decomposition must be to further levels and it is not sufficient with two frequency bands. From this need, multilevel filters appear. Multilevel filters repeat the filtering process iteratively with the output signals from the previous level. This leads to the so called wavelet decomposition trees (Figure 3.) [2]. By decomposing a signal in more frequency bands, additional information is obtained. A suitable branch to each signal is highly recommended as more decompositions do not always mean higher quality results.

Figure 3.

Wavelet decomposition tree.

The calculation of the Continuous Wavelet Transform starts for an initial time and a scale value. The result of multiplying the two signals is integrated into the whole space of time. Subsequently, this integral is multiplied by the inverse of the square root scale value, obtaining a transformed function with a normalized energy. This process is iterative until the end of the original signal is reached and must be repeated for all the values of scale that sweep the frequency range to be studied.

2.1. Wavelet families

The concept of wavelet has emerged and evolved during the last decades. Though new families of wavelet transforms are rapidly increasing, there are a number of them that have been established with more strength over time. In most situations, the use of a particular family is set by the application.

Daubechies wavelets are the most used wavelets, representing the foundations of wavelets signal processing and founding application in Discrete Wavelet Transform. They are defined as a family of orthogonal and smooth basis wavelets characterized by a maximum number of vanishing moments. The degree of smoothness increases as long as the order is higher. Daubechies wavelets lead to more accurate results in comparison to others wavelet types and also handle with boundary problems for finite length signals in an easier way [58] [29] [60] [94]. Wavelets have not an explicit expression except for order 1, which is the Haar wavelet. The inability to present a wavelet equation by a particular formula will be the general trend for almost all types of wavelet families [76].

As above mentioned, Haar wavelets are Daubechies wavelets when the order is 1. They are the simplest orthonormal wavelets. The main drawback for Haar wavelets is their discontinuity as a consequence of not solving breaking points problems for its derivates. The Haar transform is one of the earliest examples of a wavelet transform and it is supported by a function is an odd rectangular pulse pair [33]. Haar functions are widely used for applications as image coding, edge extraction and binary logic design and are defined as [46] [41] [34] [30]:

$$H(t)=\{\begin{array}{ccc}1& 0\le t\frac{1}{2}& \\ -1& \frac{1}{2}\le t1& \\ 0& elsewhere& \end{array}$$ E5

The main advantages of the Haar wavelet are its accuracy and fast implementation compared with others methods, its simplicity and small computational costs, and its capacity for solving boundaries problems [87].

Symlet wavelet transform is an orthogonal wavelet defined by a scaling filter (a low-pass finite impulse response filter of length 2N and sum 1). Symlet wavelet transform is sometimes called SymletN, where N is the order. Symlet wavelets are near symmetric. Furthermore, they have highest number of vanishing moments for a given width [7].

Coiflet wavelets are a family of wavelets whose main characteristics are similar to the Symlet ones: a high number of vanishing moments and symmetry. Coiflet family is also compactly supported, orthogonal and capable to give a good accuracy when the original signal has a distortion. The Coiflet wavelets are defined for 5 orders [18].

Biorthogonal wavelets have become very popular because of its versatility, being capable of supporting symmetric or antisymmetric signals. They perform very well under certain boundaries conditions [97]. Moreover the Biorthogonal wavelet transform is an invertible transform. They have two sets of lowpass filters for reconstruction, and highpass filters for decomposition [32].

Along with the Haar wavelets, the Meyer family is one of the exceptions that can be represented by an equation. The Meyer wavelets have numerous applications in the theory of functions, solving differential equations, signal processing, etc. [39]. Meyer family has not compact support being this one of its drawbacks. It is defined by equation (6) [44]:

$${\lambda (\omega )|}_{[0,\infty ]}=\{\begin{array}{ll}\frac{\pi }{4}+\theta (\omega -\pi ),\hfill & \omega \in \left[\frac{2\pi }{3},\frac{4\pi }{3}\right],\hfill \\ \frac{\pi }{4}+\theta (\frac{\omega }{2}-\pi ),\hfill & \omega \in \left[\frac{4\pi }{3},\frac{8\pi }{3}\right],\hfill \\ 0,\hfill & \omega \in \left[0,\frac{2\pi }{3}\right]\cup \left[\frac{8\pi }{3},+\infty \right],\hfill \end{array}$$ E6

where $$\theta (\omega )$$ is a continuously and differentiable function equal to $$\frac{\pi }{4}$$ for $$\omega \ge \frac{\pi }{3}$$ .

2.2. Wavelet transform applications

The use of the wavelet transform has been developed over the past two decades focused on the process diagnosis and instrumentation. In 1990, Leducq introduces them in the analysis of hydraulic noise for a centrifugal pump [45]. Later other authors demonstrates its usefulness for the detection of mechanical failures and the health monitoring control in gears [74] [11] [90] [21] [82] [80]. Cracks in rotors [1], structures [73] [63] [89] [10] or composite plates [75] has been another exploitation source for wavelet transforms. In 1994, Newland researches on their properties and applications, and coins the term harmonic wavelet. Harmonic wavelets are used for ridge and phase identification in signals [54]. The results showed that the cracks found reduced the rotor speed. The effectiveness of wavelets has also been compared with the envelope detection methodology in the diagnosis of faults in the bearings, obtaining results in shorter time analysis [85].

Due to its good analytical skills in time regarding the frequency, wavelet transform is a guarantee of success in the study of transient processes. Chancey and Flowers [16] managed to discover a relation between vibration patterns and the coefficients of a wavelet. Kang and Birtwhistle [40] or Subramanian, Badrilal and Henry [78] developed techniques to find problems in power transformers. Yacamini [96] proposed a method to detect torsional vibrations in engines and generators from the stator currents.

At present, the development of techniques associated to the scopes mentioned previously are still being implemented but others wavelet transforms purposes are emerging, such as classification of linear frequency modulation signals for radar emitter recognition [83] or applications to damages caused by corrosion in chemical process installations [86]. As follow there is an explanation for some of the most examined in the scientific literature.

The application of wavelets transforms in wind turbines focuses on the implementation of adaptive controllers for wind energy conversion systems. Wavelet transform is capable of providing a good and quick approximation. The drivers studied under different noise levels achieved higher performances [69]. Other works study the monitoring and diagnosis of faults in induced generators with satisfactory results. In these cases a combination of DWTs, accompanied by statistical data and energy is proposed. The use of decomposed signals spectral components is other highly interesting technique of study. Its harmonic content has suitable characteristics to be employed in fault diagnosis as an alternative to conventional methods [3].

Rolling bearing plays an important role in rotating machines. The choice of a particular wavelet family is crucial for the maintenance and fault diagnosis. The location of peaks on the vibration spectrum can identify a particular fault. Wavelet decomposition trees are a useful tool for this identification. The mean square error extracted from the terminal nodes of a tree reports the failure and its size [17]. There are also studies focused on determining what type of wavelet is suitable for bearing maintenance [79].

The wavelet transform is a good signal analysis method when a variation of time but not of space exists. The analysis provides information about the frequency of the signal, being a solution for the engine failure detection. There are detection algorithms that identify the presence of a fault in working condition and are ahead of the shutdown of the system, reducing costs and downtimes [19] [20]. These algorithms are independent of the type of engine used. Other studies in this field, present methods to detect imbalances in the stator voltage of a three phase induction motor. The wavelet transform of the stator current is analysed. Computationally, these methods are less expensive than other existing and can detect faults in an early stage. In the same vein, monitoring fatigue damage has been studied [65].

3. Condition Monitoring for engine-generator mechanism

A novel approach for Condition Monitoring based on wavelet transforms is introduced. A system for a mechanism based on an engine and a generator will be shown. It has been designed to represent any similar mechanism located in a wind turbine, generally in the nacelle. These mechanisms are used in cooling devices (generators, gearboxes), electric motors for service crane, yaw motors, pitch motors (depending on the configuration) or pumps (oil, water) according to the sub systems configurations, ventilators, etc (Figure 4).

A set of faults are induced in different experiments: ski-slope faults, misalignment faults, angular misalignment faults, parallel misalignment faults, rotating looseness faults and external noise faults. Pattern recognition is obtained from the extraction of vibration and acoustic signals. A Fault Detection and Diagnosis method is developed from the patterns of these signals. In order to recognize the patterns, three basic steps have been followed [37]:

1. The data acquisition on the testing bench (Figure 5).

2. The extraction of the features of the experiment using specific algorithms.

3. A decision-making.

A classification has been done to obtain the optimal pattern recognitions employing the data from Fast Fourier Transform and wavelet transforms applied to the vibrations and sounds signals respectively.

Figure 4.

Different locations of a wind turbine where the CM can be used: (1) fans, (2) gear oil pump, (3) oil pump for brake and (4) water cooling pump.

3.1. Case study

The experiments were made on a mechanism consisting of an engine and a generator linked by an elastic coupling joint. The sensors employed were a current sensor, an ambient temperature sensor, another temperature sensor located in strategic points of the mechanism, a vibration sensor; and a sound sensor (microphone). The data obtained by these sensors are stored in a data acquisition board, except for the vibration which is collected directly with a vibrometer. The software employed was LabView and specific software for vibration provided by the manufacturer Kionix. The speed of the engine and its associated frequency were set by a frequency variator, and the energy is dispelled using a resistive element.

The allocation of the vibration measurements were: two points for the engine and two for the generator. Points of selection were located at the end of each machine and as close as possible to the axis which is the main rotational element of the mechanism (Figure 6).

Figure 5.

Experimental mechanism.

Figure 6.

Measuring points.

The experiments were completed for an average time of 10 seconds each one, and every experiment was repeated 3 times. Therefore, for each experiment 12 measurements of temperatures, currents, sound, velocities and vibrations were taken (Figure 7). In the case of vibration, the vibrometer is capable of storing samples for the ‘x’, ‘y’ and ‘z’ axis, in addition to a total measurement for the point studied (Figure 8).

The experiments were carried out in order to identify couplings and misalignments in different degrees. The engine has 4 rubber clamping (silemblocks), while the generator has 3 rubbers clamping. The silemblocks were located at the ends, having two on the right side of the engine and two on the left side. The generator has them placed in a triangle, two in the area closest to the coupling and one at the end. The first experiment recorded under free fault conditions, and the rest of experiments were performed when the silemblocks were removed from the engine and the generator in order to create the different degrees of decoupling (Figure 9).

Figure 7.

Data collection in LabView.

Figure 8.

Data collection with Kionix software (vibration).

The rotational speed is 1500 rpm, i.e. 25 Hz. In order to do an analysis above the natural frequency, the number of samples was increased from 25 Hz to 125 Hz, being 25 Hz the default samples. This guarantees a range 5 times bigger than the natural frequency of the engine.

Experiment	Type of experiment	Data set
1	Free fault conditions	From 1 to 12
2	Misalignment removing silemblocks from the right side of the engine	From 13 to 24
3	Misalignment removing silemblocks from the right side and the front left one of the engine	From 25 to 36
4	Generation of resistance in the coupling	From 37 to 48
5	Misalignment removing the silemblock from the right side of the generator	From 49 to 60
6	Misalignment removing 2 silemblocks near to the coupling in the generator	From 61 to 72
7	Misalignment removing the silemblock from the right side of the generator and one from the left side of the engine	From 73 to 84
8	Use of a rigid coupling	From 85 to 96

Table 1.

Experiments (1500 rpm).

The FFT of each signal has been developed in Matlab. An algorithm that allows the comparison of two signals for a given frequency was created. The main purpose is to compare pattern conditions with the signals of the rest of experiments that represent a fault and to analyse the peaks found in the natural frequency and its multiples. In some cases it is important to analyse the area located below the natural frequency. Another advantage of the program is that it is possible to obtain the amplitude values for a certain frequency range (Figure 10). With a click on a particular peak, the program provides the data.

Figure 9.

Misalignments induced removing silemblocks from the engine and the generator and experimentation with a rigid coupling.

Values for 25 Hz (natural frequency or 1X), 50 Hz (2X), 75 Hz (3X) and 100 Hz (4X) have been taken into account. Frequencies above these values have been discarded.

Figure 10.

FFT of a vibration signal.

3.2. Vibration diagnosis and results

The most common spectrums for engine-generator mechanisms are presented. Examples based on the experiments held are shown.

Ski-slope fault

A ski-slope fault appears when the spectrum begins at a high level and then it goes down slowly (Figure 11). A ski-slope shows a problem with the quality of the sensor. Sometimes it happens because the sensor has experienced a transient during the measurement process. The transient may be mechanical, thermal or electrical.

Misalignment faults

Misalignment fault appears when the centrelines of coupled shafts do not coincide. If the misaligned shaft centrelines are parallel but not coincident, then the misalignment is a parallel misalignment. If the misaligned shafts meet at a point but they are not parallel, the misalignment is angular. Most of the cases are a combination of them. The diagnosis is based on dominant vibration from the natural frequency (1X) at twice the rotational rate (2X), with increased rotational rate levels (3X, 4X, etc.) acting in the axial, vertical or horizontal directions.

Angular misalignment fault

Angular misalignment fault produces a bending moment on both shafts and this generates a strong vibration at 1X, and some others at 2X and 3X for the axial direction. There will also be strong radial components for vertical and horizontal directions (Figure 11).

Parallel misalignment fault

Parallel misalignment fault produces a shear force and a bending moment on the coupled end of each shaft. High vibration levels at 2X as well as 1X are produced in the radial direction. Most often the 2X component is higher than 1X. Depending on the coupling, there can be 3X or 4X, even reaching 8X when the misalignment is severe (Figure 11).

Rotating looseness fault

Rotating looseness fault will create harmonics or sub-harmonics every 0.5X. Even 1/3 order harmonics are possible (Figure 11).

External noise fault

It is very common to find a peak in a spectrum that is difficult to analyse. This happens because of the vibration from another machine or process. The peak will typically be at a non-synchronous frequency (Figure 11). External noise can be verified stopping the machine (or varying the speed) and seeing if the vibration is still present or checking local machines for the same frequency source.

Figure 11.

a) Angular misalignment fault (red) and pattern condition (blue), (b) parallel misalignment fault (red) and pattern condition (blue), (c) ski-slope fault (blue) and pattern condition (red) and (d) rotating looseness (blue); and external noise (red).

3.3. Vibration results

As a rule, the natural frequency (1X) has been kept as the reference. Following the same nomenclature, the peaks at 50 Hz, 75 Hz and 100 Hz have been named 2X, 3X and 4X.

Figure 12.

Vibration for point 1.

Figure 13.

Vibration for point 2.

Vibration patterns are different for the four operating points. It has been detected that the natural frequency, regardless of its amplitude, tends to predominate in the experiments associated with the end points of the set (Figures 12 and 15). Additionally, the generator’s closest point to the coupling also has a similar pattern (Figure 14). The second point differs from the rest, yielding most predominant peaks from the frequency at 50 Hz (Figure 13). To make the vibration analysis, it must be taken into consideration not only the appearance of peaks, but also the amplitude. The same diagnosis for two experiments can vary its amplitude depending on the severity of the faults found. The main symptoms appear when peaks at 0.5X, 1X, 2X and 3X, sidebands and noise sources are detected. When a failure is studied at an advanced stage, peaks at 4X are noticeable (case of rigid coupling).

Figure 14.

Vibration for point 3.

Figure 15.

Vibration for point 4.

The diagnosis of the experiments reveals that the mechanism has a minor looseness which causes the appearance of a high peak at the natural frequency in some cases, even under free fault conditions. This looseness appears because the engine and the generator are not anchored directly to the test bench. The assembly was done on a surface that has facilitated the removal of the silemblocks when the experiments required it, e.g. to create different degrees of misalignment. On the other hand, this action expands the vibration intentionally because in this way it is closer to the actual behavior of the nacelle. These frequency peaks change their trend in 1X as long as the study advances from the end of the engine to the generator. From point 2, the peak at frequencies as 2X and 3X becomes more significant and sometimes exceed the amplitude of the natural frequency.

The results for experiment 8 are also remarkable. The rigid coupling added causes a severe looseness and vibration. The growth of a frequency at 4X and a constant noise over the spectrum is observed. Although it is usual to find sidebands, peaks below 1X and high frequency peaks for all this type of experiments, this feature is unique to this last experiment. Initially, a similar diagnosis for cases 1, 4 and 8 was expected, but the behavior has been slightly different for this reason.

3.4. Wavelet transform processing approach and results

Wavelet transforms were employed to analyse the sound signals. As for the Fast Fourier Transform, an algorithm has been written with Matlab. This program plots and compares two signals. Data has been transformed in 5 decompositions named a ₄ , d ₄ , d ₃ , d ₂ and d ₁ , where each of them has an energy rate associated from the original signal (Figure 16). The algorithm also returns a percentage value per decomposition. These values of energy, the decomposition levels attached and the peak amplitudes are examined in order to look for patterns.

Functions in the time domain can be represented as a linear combination of all frequency components present in a signal, where the coefficients are the amount of energy provided by each frequency component to the original signal. The main decomposition is associated with a ₄ (main or mother wavelet) that usually has the highest energy, though it is not always necessarily the case. It has a similar pattern to the original signal. The first (d ₄ ), second (d ₃ ), third (d ₂ ) and fourth (d ₁ ) transformed signals have decreasing energy rates, being s the original signal. Usually a ₄ is the low frequency component of the original signal while d _i is the high frequency component, having d ₁ the biggest value.

It is necessary to verify that the experiments performed at 1500 rpm can be extrapolated to other speeds. In the case of wind turbines, most of the engines rotate at speeds close to 3000 rpm. A certain number of tests were done varying from 500 to 3000 rpm (at intervals of 500 rpm) in order to ensure the existence of the proportional pattern.

The results showed that regardless of the speeds or the points of study, all the graphical representations for the different decompositions of energy had the same patterns. Figure 17 indicates the existence of a similar behavior where only changes the numerical value. The biggest ones will correspond to the main signals, while the results for decompositions d ₁ and d ₂ are similar.

Data can be studied according to the evolution of a single point along the different experiments or analysing the evolution of the set points for all the experiment. Each row in Figure 18 contains two graphics, one with the amplitude peaks (left) and the other one with the energy distribution of the sound signal (right). The first two graphics correspond to the engine end (point 1). The following two graphics are the closest to the coupling (point 2). The third row belongs to the points of the generator next to the coupling (point 3), and finally, the last two graphics are for the end of the generator (point 4).

Figure 16.

Wavelet decompositions.

Figure 17.

Energies at different rotational speeds.

Figure 18.

Evolution of the frequency peaks and wavelet energy decompositions for each point in experiment 2.

Based on the distribution of the energy among the 5 different decompositions of every signal, the energy distribution for point 1, end of the engine-generator set is ruled by an almost similar pattern where each experiment has a maximum of energy in the main signal and a minimum for decomposition d ₁ or d ₂ . It means that by performing a decomposition of the signal, the energy has a closest resemblance to the original value, often exceeding 85% of the total energy, remaining a residual percentage for d ₁ or d ₂ . When the experiments are closer to the generator (points 2, 3 and 4), the energy is distributed among the 5 decompositions and not concentrated in the mother wavelet, as it is for point 1.

All the decompositions have been registered with their energy maximum and minimum values and their patterns distribution. An example for 2 experiments is shown in Table 2.

Experiment	Main	d4	d3	d2	d1	Energy
A	17,19%	9,10%	22,12%	24,95%	26,63%	167,9
B	82,32%	9,90%	5,17%	1,77%	0,84%	311,8

Table 2.

Energy distribution for experiments A and B.

Figure 19.

Energy values for point 1.

Figure 20.

Energy values for point 2.

Experiment A is associated to point 2, belonging to the engine and situated close to the coupling. Experiment B, however, is related to point 1, left end of the assembly. Experiment A has the maximum percentage of energy in d ₁ and the minimum in d _4. Furthermore, the experiment B has its maximum in the main signal and the minimum located in d ₁ . The maximum-minimum patterns are d ₁ -d ₄ and main-d ₁ respectively. Numerically, the most compensated distribution of energy is close to the coupling (experiment A – point 2) above mentioned. The patterns main-d ₁ and main-d ₂ appear for all the cases in point 1. However, the same maximum-minimum distribution is smaller for the points 2, 3 and 4. Unlike in point 1, there are different patterns for the 8 experiments in these points. Figures 19, 20, 21 and 22 represent the numerical values of the energy per point and experiment. It must be noted that the numerical values are higher or lower, depending on the type of experiment.

Figure 21.

Energy values for point 3.

Figure 22.

Energy values for point 4.

4. Conclusions

Wind turbines are complex systems that require a high level of reliability, availability, maintainability and safety. This chapter is focused on to guarantee these correct levels for mechanisms used in cooling devices for generators and gearboxes, electric motors for service crane, yaw motors, pitch motors, pumps, ventilators, etc.

The mechanism brake of the engine has been simulated linking a generator by a coupling joint. The signals collected have been:

• Vibration.

• Sound.

• Current.

• Temperature.

• Velocity.

The experiments have been done in working conditions for different points of the mechanism and considering the following failures:

• Misalignment removing silemblocks from the right side of the engine.

• Misalignment removing silemblocks from the right side and the front left one of the engine.

• Induction of resistance in the coupling.

• Misalignment removing the silemblock from the right side of the generator.

• Misalignment removing 2 silemblocks near to the coupling in the generator.

• Misalignment removing the silemblock from the right side of the generator and one from the left side of the engine.

• Using a rigid coupling.

A fault detection and diagnosis model based on the Fast Fourier Transform applied to the vibration signals; together with the wavelet transform applied to sound signals has been developed. The model detects and diagnoses correctly 100% of the failures considered.

It has been observed that for the outer ends of the engine and the generator, the appearance of a pronounced peak amplitude at the natural frequency or 2X (vibration) was associated to the maximum energy values for the main signal, the most suitable with the original, and minimum values for decomposed signals d ₁ and d ₂ (sound). In contrast, the results obtained close to the coupling did not follow a clear trend as the results were conditioned by the type of experiment. The numerical values of each peak were also taken into account in the establishment of the pattern recognitions, being different for each experiment. The same conclusion was reached for the energy values. Different models and results were expected because the objective was not to find similar patterns between different experiments, and the tests were never performed under identical conditions. The objective was to have different vibration patterns and their associated sound models in order to create a catalogue of possible scenarios for predictive maintenance in the mechanisms. Thus, it is possible to extend the range of possibilities to relate the result of an acoustic signal with the frequency domain using the Fast Fourier Transform.