Quantitative Index and Abnormal Alarm Strategy Using Sensor-Dependent Vibration Data for Blade Crack Identification in Centrifugal Booster Fans

Centrifugal booster fans are important equipment used to recover blast furnace gas (BFG) for generating electricity, but blade crack faults (BCFs) in centrifugal booster fans can lead to unscheduled breakdowns and potentially serious accidents, so in this work quantitative fault identification and an abnormal alarm strategy based on acquired historical sensor-dependent vibration data is proposed for implementing condition-based maintenance for this type of equipment. Firstly, three group dependent sensors are installed to acquire running condition data. Then a discrete spectrum interpolation method and short time Fourier transform (STFT) are applied to preliminarily identify the running data in the sensor-dependent vibration data. As a result a quantitative identification and abnormal alarm strategy based on compound indexes including the largest Lyapunov exponent and relative energy ratio at the second harmonic frequency component is proposed. Then for validation the proposed blade crack quantitative identification and abnormality alarm strategy is applied to analyze acquired experimental data for centrifugal booster fans and it has successfully identified incipient blade crack faults. In addition, the related mathematical modelling work is also introduced to investigate the effects of mistuning and cracks on the vibration features of centrifugal impellers and to explore effective techniques for crack detection.


Introduction
Blast furnace gas (BFG) is a byproduct of iron-making. With the great expansion of the iron and steel industry, the production of BFG during iron-making has increased remarkably [1]. However, blast furnace gas is characterized in low calorific value, difficulty to burn and combustion instability as a power fuel [2], so how to deal effectively with blast furnace gas is a problem that puzzles iron and steel enterprises. In recent years, many enterprises have taken to use BFG to generate electricity in order to save energy and improve benefits. Centrifugal booster fans are important pieces of equipment which are used to pressurize BGF and make sure it can go into the furnace safely and combust stably. Faults occurring on the centrifugal booster fan may lead to accidents such as unstable combustion in the furnace, blow outs, downtime, and potentially huge economic losses.
Blade crack faults (BCFs) are among the typical faults in centrifugal booster fans. Different cracks arise after long running due to the resonance, decreased anti-fatigue capability because of manufacture problems, installation issues or the work conditions [3]. This may result in blades breaking off and the unit being damaged. For a rotor system, the stiffness of the shaft would display cyclical behavior The technological process of the BFG recovery unit is as follows: The original BFG gas from the the Department of Energy and Environment Protection → Venturi meter → BFG entrance shutting door → BFG booster fan entrance control damper → three BFG booster fans → BFG booster fan exit isolation damper → Heater and its bypass → BFG gas burners.
The speed of the booster fans is divided into two levels: 744 r/min and 993 r/min, whose corresponding powers are 730 kW and 1650 kW, respectively.

Condition Data Acquisition Testing Framework
For iron-making, lots of rotating machineries work in the complex process going from iron ore to steel products. The equipment is in long-term use under complex conditions. This may lead to various types of fault and cause huge economy losses. However, the key parts of rotating machinery are not stationary and therefore not easy to change, so it is crucial to carry out effective condition monitoring and fault diagnosis.
For this purpose, a testing system focusing on extracting abnormal condition information from vibration signals is designed. This system works on the data acquisition level and network database level of the testing and diagnosis system in the steel-making plant. Off-line vibration tests are conducted for iron and steel smelting mechanical equipment using portable data acquisition devices such as the CSI2320, Telesens8823, SONY-EX and so on. These devices can selectively implement hardware integration for acquired signals to save the data.
The feature information can be extracted from the vibrations by means of traditional spectrum analysis, characteristic spectrum analysis and special feature extraction modules in the testing system. The change trends of these features can be used to judge the working condition of the equipment and the appearance of incipient faults. The software interface of the testing system is implemented based on Labview 7.0. The testing system is programmed with a mixture of Labview and Visual Studio routines with consideration to execution efficiency. In addition, SQL Server is introduced as the extended interface to access the internal database. Moreover, simple tips about the main functions are available in a help module. The system also has other functions such as saving results, report generation, etc. The whole work flow diagram of the quantitative identification research framework is displayed in Figure 2. The technological process of the BFG recovery unit is as follows: The original BFG gas from the the Department of Energy and Environment Protection Ñ Venturi meter Ñ BFG entrance shutting door Ñ BFG booster fan entrance control damper Ñ three BFG booster fans Ñ BFG booster fan exit isolation damper Ñ Heater and its bypass Ñ BFG gas burners.
The speed of the booster fans is divided into two levels: 744 r/min and 993 r/min, whose corresponding powers are 730 kW and 1650 kW, respectively.

Condition Data Acquisition Testing Framework
For iron-making, lots of rotating machineries work in the complex process going from iron ore to steel products. The equipment is in long-term use under complex conditions. This may lead to various types of fault and cause huge economy losses. However, the key parts of rotating machinery are not stationary and therefore not easy to change, so it is crucial to carry out effective condition monitoring and fault diagnosis.
For this purpose, a testing system focusing on extracting abnormal condition information from vibration signals is designed. This system works on the data acquisition level and network database level of the testing and diagnosis system in the steel-making plant. Off-line vibration tests are conducted for iron and steel smelting mechanical equipment using portable data acquisition devices such as the CSI2320, Telesens8823, SONY-EX and so on. These devices can selectively implement hardware integration for acquired signals to save the data.
The feature information can be extracted from the vibrations by means of traditional spectrum analysis, characteristic spectrum analysis and special feature extraction modules in the testing system. The change trends of these features can be used to judge the working condition of the equipment and the appearance of incipient faults. The software interface of the testing system is implemented based on Labview 7.0. The testing system is programmed with a mixture of Labview and Visual Studio routines with consideration to execution efficiency. In addition, SQL Server is introduced as the extended interface to access the internal database. Moreover, simple tips about the main functions are available in a help module. The system also has other functions such as saving results, report generation, etc. The whole work flow diagram of the quantitative identification research framework is displayed in Figure 2. 1~, The work flow diagram of the research framework.

Discrete Spectrum Interpolation Method
The rotating frequency and its harmonic components can reflect the fault features concerning misalignment, rub-impact, and dynamic unbalance. These features are usually extracted from the frequency spectrum. However, the FFT and spectrum analysis would cause energy leakage owing to time domain truncation and the interference of noise [26,27]. This may lead to great errors in the frequency, amplitude and phase in the FFT and spectrum analysis. In order to improve the accuracy, a discrete spectrum interpolation method is adopted.
Let x(t) be a harmonic signal sequence with frequency f0, amplitude A0 and phase θ0. Suppose the amplitude and phase first calculated by a Discrete Fourier Transform (DFT) and then corrected by the interpolation method are: where ( ) w X k means the kth line of the harmonic signal, i.e., the maximum value of the main lobe.
W(f 1 ) expresses the frequency spectrum for a rectangular window with the value W(f 1 ) = sin(f 1 )/ (f 1 ). Rk and Ik represent the real and imaginary parts of DFT, respectively. When rotating frequency is input, the system would correct the amplitude and phase of the rotating frequency and search for the accurate amplitude and phase of its harmonic components automatically.

Short Time Fourier Transform
Blade crack fault diagnosis is a problem that troubles scholars and engineers at home and abroad. Online and offline condition monitoring systems are widely used in rotating machinery. The key problem is how to choose the method to process the signal from the industrial field and obtain the fault characteristics. The short-time Fourier transform (STFT) is one of the earliest and the most basic methods used for time-frequency analysis [28]. The STFT is one of the most widely used algorithms in signal processing and fault diagnosis based on a detailed Fourier transform centered at each time point. In STFT, the signal is compared with window functions that are concentrated in

Discrete Spectrum Interpolation Method
The rotating frequency and its harmonic components can reflect the fault features concerning misalignment, rub-impact, and dynamic unbalance. These features are usually extracted from the frequency spectrum. However, the FFT and spectrum analysis would cause energy leakage owing to time domain truncation and the interference of noise [26,27]. This may lead to great errors in the frequency, amplitude and phase in the FFT and spectrum analysis. In order to improve the accuracy, a discrete spectrum interpolation method is adopted.
Let x(t) be a harmonic signal sequence with frequency f 0 , amplitude A 0 and phase θ 0 . Suppose the amplitude and phase first calculated by a Discrete Fourier Transform (DFT) and then corrected by the interpolation method are:Â where X w pkq means the kth line of the harmonic signal, i.e., the maximum value of the main lobe. W(∇f 1 ) expresses the frequency spectrum for a rectangular window with the value W(∇f 1 ) = sin(π∇f 1 )/ (π∇f 1 ). R k and I k represent the real and imaginary parts of DFT, respectively. When rotating frequency is input, the system would correct the amplitude and phase of the rotating frequency and search for the accurate amplitude and phase of its harmonic components automatically.

Short Time Fourier Transform
Blade crack fault diagnosis is a problem that troubles scholars and engineers at home and abroad. Online and offline condition monitoring systems are widely used in rotating machinery. The key problem is how to choose the method to process the signal from the industrial field and obtain the fault characteristics. The short-time Fourier transform (STFT) is one of the earliest and the most basic methods used for time-frequency analysis [28]. The STFT is one of the most widely used algorithms in signal processing and fault diagnosis based on a detailed Fourier transform centered at each time point. In STFT, the signal is compared with window functions that are concentrated in both the time and STFT x pτ, f q " where, w(t) is the window function which has a user defined time duration; and x(t) is the waveform signal in the time domain.

Case Study via Blade Crack Fault of a Centrifugal Booster Fan
Unit 4 in a power station composed of three imported fans. As shown in Figure 3 [24,31,32], the rotor blades in No. A fan are welded on the entrance control damper. The three fans performed well since they were first used in production, and had never been overhauled before. On 20 July 2011, the rotor broke apart during the process of switching from low speed to high speed. Pieces of the blades flew out of the volute. Figures from the scene are shown in Figure 4 [24,31,32]. The bearing box in the drive end and the coupling are crushed. Besides, the main shaft is seriously deformed. By analysis of the causes that produced the accident, we find that there were blade cracks in the booster fans.
where, w(t) is the window function which has a user defined time duration; and x(t) is the waveform signal in the time domain.

Case Study via Blade Crack Fault of a Centrifugal Booster Fan
Unit 4 in a power station composed of three imported fans. As shown in Figure 3 [24,31,32], the rotor blades in No. A fan are welded on the entrance control damper. The three fans performed well since they were first used in production, and had never been overhauled before. On 20 July 2011, the rotor broke apart during the process of switching from low speed to high speed. Pieces of the blades flew out of the volute. Figures from the scene are shown in Figure 4 [24,31,32]. The bearing box in the drive end and the coupling are crushed. Besides, the main shaft is seriously deformed. By analysis of the causes that produced the accident, we find that there were blade cracks in the booster fans.  In order to extract the vibration characteristics during crack growth and gain experience for condition monitoring on the same type of unit, the testing system as described before is used to analyze the historical data. From 4 August 2010 to 6 July 2011 (nearly a year before the accident), data acquisition of the two bearings, which support the rotor, was carried out by the industry technological service company at the sampling frequency of 2560 Hz with the length of 4096, ten times, including low speed (744 r/min) and high speed (993 r/min). Three groups of sensors are used, as shown in Figure 5. both the time and frequency domains. The STFT algorithm and the window function can be mathematically represented as follows: where, w(t) is the window function which has a user defined time duration; and x(t) is the waveform signal in the time domain.

Case Study via Blade Crack Fault of a Centrifugal Booster Fan
Unit 4 in a power station composed of three imported fans. As shown in Figure 3 [24,31,32], the rotor blades in No. A fan are welded on the entrance control damper. The three fans performed well since they were first used in production, and had never been overhauled before. On 20 July 2011, the rotor broke apart during the process of switching from low speed to high speed. Pieces of the blades flew out of the volute. Figures from the scene are shown in Figure 4 [24,31,32]. The bearing box in the drive end and the coupling are crushed. Besides, the main shaft is seriously deformed. By analysis of the causes that produced the accident, we find that there were blade cracks in the booster fans.  In order to extract the vibration characteristics during crack growth and gain experience for condition monitoring on the same type of unit, the testing system as described before is used to analyze the historical data. From 4 August 2010 to 6 July 2011 (nearly a year before the accident), data acquisition of the two bearings, which support the rotor, was carried out by the industry technological service company at the sampling frequency of 2560 Hz with the length of 4096, ten times, including low speed (744 r/min) and high speed (993 r/min). Three groups of sensors are used, as shown in Figure 5. In order to extract the vibration characteristics during crack growth and gain experience for condition monitoring on the same type of unit, the testing system as described before is used to analyze the historical data. From 4 August 2010 to 6 July 2011 (nearly a year before the accident), data acquisition of the two bearings, which support the rotor, was carried out by the industry technological service company at the sampling frequency of 2560 Hz with the length of 4096, ten times, including low speed (744 r/min) and high speed (993 r/min). Three groups of sensors are used, as shown in Figure 5. Two sensors are mounted on the two bearings. The vibration data from the horizontal, vertical and radial directions are obtained in each sensor. In addition, a sensor is mounted on the motor. The data from the vertical direction is obtained there. The details about the running status data acquisition are shown in Table 1, where D means low speed, Dh means high speed and GE means the amplitude of the envelope in g. A means radial direction, H indicates horizontal direction and V is the vertical direction.  The vibration data from the drive end in the horizontal direction is chosen for analysis. We can see that there is no a distinctive trend that would allow identifying crack faults. The waterfall plot based on FFT of No. 2H sensor data from 0 to 100 Hz is displayed in Figure 6 [23,24]. We can clearly observe the rotating frequency and its harmonic components. As time goes by, the amplitude of the rotating frequency increases first and then drops. However, the amplitude of the second harmonic drops first and then increases. From the point of view of the dynamics, the cracks close and open twice in one cycle of the rotation. The amplitude of the vibration response signal changes twice due to the change of the rotor stiffness from large to small, so the characteristics of the amplitude increase in the second harmonic could be used to indicate a blade crack.
In order to obtain more accurate amplitude and phase data of the rotating frequency and the second harmonic component, a discrete spectrum interpolation method is used to analyze the signals of the No. 2 and No. 3 sensors in the horizontal direction. The result after correction is shown in Figure 7 [23,24]. As shown in Figure 7, the amplitude of the rotating frequency first increased and then decreased. The amplitude of the second harmonic component increases constantly. A breathing crack is considered to appear from January to March in 2011 according to the results. The rotor Two sensors are mounted on the two bearings. The vibration data from the horizontal, vertical and radial directions are obtained in each sensor. In addition, a sensor is mounted on the motor. The data from the vertical direction is obtained there. The details about the running status data acquisition are shown in Table 1, where D means low speed, Dh means high speed and GE means the amplitude of the envelope in g. A means radial direction, H indicates horizontal direction and V is the vertical direction.  The vibration data from the drive end in the horizontal direction is chosen for analysis. We can see that there is no a distinctive trend that would allow identifying crack faults. The waterfall plot based on FFT of No. 2H sensor data from 0 to 100 Hz is displayed in Figure 6 [23,24]. We can clearly observe the rotating frequency and its harmonic components. As time goes by, the amplitude of the rotating frequency increases first and then drops. However, the amplitude of the second harmonic drops first and then increases. From the point of view of the dynamics, the cracks close and open twice in one cycle of the rotation. The amplitude of the vibration response signal changes twice due to the change of the rotor stiffness from large to small, so the characteristics of the amplitude increase in the second harmonic could be used to indicate a blade crack.
In order to obtain more accurate amplitude and phase data of the rotating frequency and the second harmonic component, a discrete spectrum interpolation method is used to analyze the signals of the No. 2 and No. 3 sensors in the horizontal direction. The result after correction is shown in Figure 7 [23,24]. As shown in Figure 7, the amplitude of the rotating frequency first increased and then decreased. The amplitude of the second harmonic component increases constantly. A breathing crack is considered to appear from January to March in 2011 according to the results. The rotor stiffness changes twice in every rotating cycle. This leads to the increase of the amplitude of the second harmonic component. stiffness changes twice in every rotating cycle. This leads to the increase of the amplitude of the second harmonic component.  For the purpose of extracting the non-stationary characteristics of the cracked rotor, STFT is employed to process the signal of the normal and cracked rotor. The result is displayed in Figure 8. From the figure, we can find a component that is obviously equal to the rotation frequency of the fan, which is 12.5 Hz. When cracks appear, there is a discontinuous line at the frequency of the second harmonic, which fluctuates at the rotating frequency and its second harmonic, as shown in Figure 9.    For the purpose of extracting the non-stationary characteristics of the cracked rotor, STFT is employed to process the signal of the normal and cracked rotor. The result is displayed in Figure 8. From the figure, we can find a component that is obviously equal to the rotation frequency of the fan, which is 12.5 Hz. When cracks appear, there is a discontinuous line at the frequency of the second harmonic, which fluctuates at the rotating frequency and its second harmonic, as shown in Figure 9.  For the purpose of extracting the non-stationary characteristics of the cracked rotor, STFT is employed to process the signal of the normal and cracked rotor. The result is displayed in Figure 8. From the figure, we can find a component that is obviously equal to the rotation frequency of the fan, which is 12.5 Hz. When cracks appear, there is a discontinuous line at the frequency of the second harmonic, which fluctuates at the rotating frequency and its second harmonic, as shown in Figure 9.   For the purpose of extracting the non-stationary characteristics of the cracked rotor, STFT is employed to process the signal of the normal and cracked rotor. The result is displayed in Figure 8. From the figure, we can find a component that is obviously equal to the rotation frequency of the fan, which is 12.5 Hz. When cracks appear, there is a discontinuous line at the frequency of the second harmonic, which fluctuates at the rotating frequency and its second harmonic, as shown in Figure 9.  The photo of the cracked rotor on the A fan is displayed in Figure 10. In this section, two methods are applied to analyze the blade crack fault. The vibration characteristics of the crack blade are extracted. This can help to improve the qualitative diagnosis performance for blade cracks, but the detection of crack damage based on the vibration signals is still less studied, and especially aquantitative identification method for scheduling reasonable maintenance plans is lacking.

The Proposed Quantitative Identification Method
Both Sun and Chen have attempted to propose a quantitative identification index for blade crack identification and have obtained certain achievements [24,31], but a more comprehensive abnormality alarm strategy via a quantitative identification index should be proposed to indicate the unbalancedness and implement necessary condition-based maintenance, so a quantitative identification method based on compound indexes including a traditional index and new index is developed in this section.

Largest Lyapunov Exponent Algorithm
Lyapunov exponents, which measure the exponential rates of divergence or convergence of nearby trajectories in state space, are generally calculated to characterize chaotic processes. If the largest value in the spectrum of Lyapunov exponents is positive, it means that the system is chaotic. The largest value equal to zero indicates periodic or quasi-periodic dynamics. If all Lyapunov The photo of the cracked rotor on the A fan is displayed in Figure 10. In this section, two methods are applied to analyze the blade crack fault. The vibration characteristics of the crack blade are extracted. This can help to improve the qualitative diagnosis performance for blade cracks, but the detection of crack damage based on the vibration signals is still less studied, and especially aquantitative identification method for scheduling reasonable maintenance plans is lacking. The photo of the cracked rotor on the A fan is displayed in Figure 10. In this section, two methods are applied to analyze the blade crack fault. The vibration characteristics of the crack blade are extracted. This can help to improve the qualitative diagnosis performance for blade cracks, but the detection of crack damage based on the vibration signals is still less studied, and especially aquantitative identification method for scheduling reasonable maintenance plans is lacking.

The Proposed Quantitative Identification Method
Both Sun and Chen have attempted to propose a quantitative identification index for blade crack identification and have obtained certain achievements [24,31], but a more comprehensive abnormality alarm strategy via a quantitative identification index should be proposed to indicate the unbalancedness and implement necessary condition-based maintenance, so a quantitative identification method based on compound indexes including a traditional index and new index is developed in this section.

Largest Lyapunov Exponent Algorithm
Lyapunov exponents, which measure the exponential rates of divergence or convergence of nearby trajectories in state space, are generally calculated to characterize chaotic processes. If the largest value in the spectrum of Lyapunov exponents is positive, it means that the system is chaotic. The largest value equal to zero indicates periodic or quasi-periodic dynamics. If all Lyapunov

The Proposed Quantitative Identification Method
Both Sun and Chen have attempted to propose a quantitative identification index for blade crack identification and have obtained certain achievements [24,31], but a more comprehensive abnormality alarm strategy via a quantitative identification index should be proposed to indicate the unbalancedness and implement necessary condition-based maintenance, so a quantitative identification method based on compound indexes including a traditional index and new index is developed in this section.

Largest Lyapunov Exponent Algorithm
Lyapunov exponents, which measure the exponential rates of divergence or convergence of nearby trajectories in state space, are generally calculated to characterize chaotic processes. If the largest value in the spectrum of Lyapunov exponents is positive, it means that the system is chaotic. The largest Sensors 2016, 16, 632 9 of 20 value equal to zero indicates periodic or quasi-periodic dynamics. If all Lyapunov exponents are negative then the stable critical point is an attractor [29,30]. Among all the Lyapunov exponents, the Largest Lyapunov Exponent (LLE) has aroused considerable interest for its significant practical applications. The LLE has been applied to many fields for its notable capabilities. In this paper, the LLE is calculated as an indicator of the chaotic behavior of the load demand by using Wolf's algorithm, which is given as [29,30]: where L'(t k ) and L'(t k-1 ) mean the Euclidean distances computed between the nearest neighboring points on the different trajectories of the attractor at the t k and t k-1 time steps, respectively [29]. m indicates the number of replacement steps or iteration number. Details on the calculation parameter selection are given in [29,30]. The negative value of LLE indicates normal condition and a positive value of LLE indicates non-linear conditions, then the value of LLE can be used to initially identify the mechanical system state. According to the principle of the LLE algorithm, the LLE value of condition data from the No. 2H sensor is computed and displayed in Figure 11. We can observe that the non-linear condition appears after 9 December 2010. In addition, we still need a quantitative identification index to confirm the degree of crack fault.
Sensors 2016, 16, 632 9 of 20 exponents are negative then the stable critical point is an attractor [29,30]. Among all the Lyapunov exponents, the Largest Lyapunov Exponent (LLE) has aroused considerable interest for its significant practical applications. The LLE has been applied to many fields for its notable capabilities. In this paper, the LLE is calculated as an indicator of the chaotic behavior of the load demand by using Wolf's algorithm, which is given as [29,30]: where L'(tk) and L'(tk-1) mean the Euclidean distances computed between the nearest neighboring points on the different trajectories of the attractor at the tk and tk-1 time steps, respectively [29]. m indicates the number of replacement steps or iteration number. Details on the calculation parameter selection are given in [29,30]. The negative value of LLE indicates normal condition and a positive value of LLE indicates non-linear conditions, then the value of LLE can be used to initially identify the mechanical system state. According to the principle of the LLE algorithm, the LLE value of condition data from the No. 2H sensor is computed and displayed in Figure 11. We can observe that the non-linear condition appears after 9 December 2010. In addition, we still need a quantitative identification index to confirm the degree of crack fault.

Relative Energy Ratio at Second Harmonic Frequency Component
According to the previous section, we know that when a breathing crack appears, the rotor stiffness changes twice from large to small in one rotation cycle [33,34]. The amplitude of the second harmonic increases obviously and the phase also changed significantly [35,36]. Hence, a quantitative identification method for blade crack fault and a new index are proposed to describe the degree of damage accurately when cracks grow. The new index is called relative energy ratio at second harmonic frequency component and expressed as K2f: where Af means the amplitude of the rotating frequency and A2f means the amplitude of the second harmonic. Moreover, we can find that the energy of the first six order harmonic accounts for more

Relative Energy Ratio at Second Harmonic Frequency Component
According to the previous section, we know that when a breathing crack appears, the rotor stiffness changes twice from large to small in one rotation cycle [33,34]. The amplitude of the second harmonic increases obviously and the phase also changed significantly [35,36]. Hence, a quantitative identification method for blade crack fault and a new index are proposed to describe the degree of damage accurately when cracks grow. The new index is called relative energy ratio at second harmonic frequency component and expressed as K 2f : where A f means the amplitude of the rotating frequency and A 2f means the amplitude of the second harmonic. Moreover, we can find that the energy of the first six order harmonic accounts for more than 98% of the total energy, so based on the amplitude of the frequency spectrum, the parameter j is selected as j = 6. The value of K 2f can be used to judge the degree of the damage, where a and b are thresholds that need to be decided: (1) 0 ď K 2f < a : Normal condition, the rotor works well; (2) a ď K 2f ď b: Incipient fault, a crack appears, but is small; (3) K 2f > b: Serious fault, many cracks appear.
During the period of August 2010 to December 2010, the rotor worked normally and we can think the reliability of this fan remained at the normal level of 0.9 during this time. The values of K 2f calculated from the bearing data in the horizontal direction are shown in Figure 12. According to the LLE value in Figure 11, a non-linear condition appeared after 9 December 2010. Moreover, according to the value of K 2f in December 2010 to March 2011 as shown in the Figure 12, the threshold value of a is preliminarily determined as 0.02. On 20 July 2011, the rotor broke apart during the process of switching from low speed to high speed (from 744 r/min to 993 r/min) and we can think that the reliability of this fan has decreased to 0 at this time. Based on the lowest reliability requirement the steel-making plant requires for safe running and the value of K 2f in bearing 2 from June 2011 to July 2011, the value of b is determined as 0.25 for safety's sake and the reliability of this fan has decreased below 0.5 at this time by linear mapping relationship analysis, so the parameters are determined completely.
The new index K 2f is introduced into the testing system and works well in the condition monitoring of the centrifugal booster fan. The results in Figure 12 could be employed to analyze the condition of No. A fan.
Based on the proposed quantitative identification method, the abnormality alarm strategy can be obtained, and then the condition-based maintenance actions can be arranged reasonably to ensure safe and reliable operation. To sum up, the proposed quantitative identification and abnormality alarm strategy procedure of using sensor-dependent vibration data for blade crack identification in centrifugal booster fans can be summarized by the flow chart displayed in Figure 13. than 98% of the total energy, so based on the amplitude of the frequency spectrum, the parameter j is selected as j = 6.
The value of K2f can be used to judge the degree of the damage, where a and b are thresholds that need to be decided: (1) 0 ≤ K2f < a : Normal condition, the rotor works well; (2) a ≤ K2f ≤ b: Incipient fault, a crack appears, but is small; (3) K2f > b: Serious fault, many cracks appear.
During the period of August 2010 to December 2010, the rotor worked normally and we can think the reliability of this fan remained at the normal level of 0.9 during this time. The values of K2f calculated from the bearing data in the horizontal direction are shown in Figure 12. According to the LLE value in Figure 11 The new index K2f is introduced into the testing system and works well in the condition monitoring of the centrifugal booster fan. The results in Figure 12 could be employed to analyze the condition of No. A fan.
Based on the proposed quantitative identification method, the abnormality alarm strategy can be obtained, and then the condition-based maintenance actions can be arranged reasonably to ensure safe and reliable operation. To sum up, the proposed quantitative identification and abnormality alarm strategy procedure of using sensor-dependent vibration data for blade crack identification in centrifugal booster fans can be summarized by the flow chart displayed in Figure 13. Meanwhile, the process of the proposed method and strategy for the mentioned engineering tasks in the power station can be summarized as follows: (1) Collect the sensor-dependent running condition vibration data; (2) Pre-process this vibration data using the discrete spectrum interpolation method and STFT; (3) Compute the LLE and the relative energy ratio K2f.

Test and Validation
We have also examined the performance of statistical features reported in the literature as comparisons to validate the performance of the proposed method. Some of the feature parameters have been demonstrated to be ineffective in previous publications, but in different papers, different feature parameters are applied according to the experience accumulated by different researchers. In the different applications, different feature parameters give different diagnosis performance. Thus, many feature parameters are calculated in this study. In total, 21 feature values are obtained, shown in Table 2. These features are adopted to indicate the faulty condition from the acquired vibration signals. The results are displayed in Figures 14-19. According to the analyzed figures, some results can be obtained. First, obvious trends cannot be found the majority of the mentioned 21 feature parameters, except for the feature values F4, F11 and F12, so they are of no use for blade crack identification in centrifugal booster fans. Next, the feature values of F4, F11 and F12 show relatively clear trends compared to the remaining feature values, but the proposed strategy using compound feature parameters can indicate the running condition of centrifugal booster fans by hierarchical descriptions. The contrastive results demonstrate the effectiveness of the proposed strategy for the engineering task at hand. Meanwhile, the process of the proposed method and strategy for the mentioned engineering tasks in the power station can be summarized as follows: (1) Collect the sensor-dependent running condition vibration data; (2) Pre-process this vibration data using the discrete spectrum interpolation method and STFT; (3) Compute the LLE and the relative energy ratio K 2f . (4) Confirm the degree of blade crack fault(s) of the centrifugal booster fan; (5) Conduct the corresponding maintenance management activities based on the fault degree.

Test and Validation
We have also examined the performance of statistical features reported in the literature as comparisons to validate the performance of the proposed method. Some of the feature parameters have been demonstrated to be ineffective in previous publications, but in different papers, different feature parameters are applied according to the experience accumulated by different researchers. In the different applications, different feature parameters give different diagnosis performance. Thus, many feature parameters are calculated in this study. In total, 21 feature values are obtained, shown in Table 2. These features are adopted to indicate the faulty condition from the acquired vibration signals. The results are displayed in Figures 14-19. According to the analyzed figures, some results can be obtained. First, obvious trends cannot be found the majority of the mentioned 21 feature parameters, except for the feature values F 4 , F 11 and F 12 , so they are of no use for blade crack identification in centrifugal booster fans. Next, the feature values of F 4 , F 11 and F 12 show relatively clear trends compared to the remaining feature values, but the proposed strategy using compound feature parameters can indicate the running condition of centrifugal booster fans by hierarchical descriptions. The contrastive results demonstrate the effectiveness of the proposed strategy for the engineering task at hand.

Category Parameter
Time-Domain index Frequency-Domain index

Category Parameter
Time-Dom ain index is the time domain signal, n = 1,2,…,N; N is the sample point; s(k) is the spectrum of x(n), k = 1,2,…,K; K is the number of spectrum lines; fK is the frequency of k-th spectrum line.         2H sensor by F4, F5 and F6. Figure 15. The analyzed results of A fan No. 2H sensor by F 4 , F 5 and F 6 .      Moreover, the proposed blade crack quantitative identification method is applied to identifying the running condition of No. A fan using A fan No. 3H sensor data. According to the principle of the Largest Lyapunov Exponent algorithm, the LLE value of the condition data from the No. 3H sensor is computed and displayed in Figure 20. We can observe that the non-linear condition also appears after 9 December 2010. In addition, the index K2f is also used to identify the running condition of No. A fan using A fan No. 3H sensor data and the results are displayed in Figure 21. From the result, we also can obtain a clear trend to indicate the blade cracks in centrifugal booster  Moreover, the proposed blade crack quantitative identification method is applied to identifying the running condition of No. A fan using A fan No. 3H sensor data. According to the principle of the Largest Lyapunov Exponent algorithm, the LLE value of the condition data from the No. 3H sensor is computed and displayed in Figure 20. We can observe that the non-linear condition also appears after 9 December 2010. In addition, the index K2f is also used to identify the running condition of No. A fan using A fan No. 3H sensor data and the results are displayed in Figure 21. From the result, we also can obtain a clear trend to indicate the blade cracks in centrifugal booster Moreover, the proposed blade crack quantitative identification method is applied to identifying the running condition of No. A fan using A fan No. 3H sensor data. According to the principle of the Largest Lyapunov Exponent algorithm, the LLE value of the condition data from the No. 3H sensor is computed and displayed in Figure 20. We can observe that the non-linear condition also appears after 9 December 2010. In addition, the index K 2f is also used to identify the running condition of No. A fan using A fan No. 3H sensor data and the results are displayed in Figure 21. From the result, we also can obtain a clear trend to indicate the blade cracks in centrifugal booster fans and the value is above the abnormal warning stage after 8 March 2011. The result indicates that the proposed index and alarm strategy is feasible.    [24] show the results. From the Figure 22a, we can see that the crack grows from the fan entrance and propagates along the radial direction of the damper. Figure 22b shows the welding position of the blade and the damper. Stress concentration in this position led to the growth and propagation of the crack. As shown in Figure 23, a crack appears at the welding position with a length of 3-4 cm in No. C fan. After detecting these cracks, the maintenance for No. B and No. C fans was carried out in March 2012 to prevent further accidents. These above results show the effectiveness and robustness of the proposed quantitative detection   The proposed blade crack quantitative identification method is also applied to identifying the running condition of the remaining centrifugal booster fans. As mentioned previously, there are three fans in the unit 4. While No. A fan was being repaired, the quantitative identification index is also used to examine the No.  [24] show the results. From the Figure 22a, we can see that the crack grows from the fan entrance and propagates along the radial direction of the damper. Figure 22b shows the welding position of the blade and the damper. Stress concentration in this position led to the growth and propagation of the crack. As shown in Figure 23, a crack appears at the welding position with a length of 3-4 cm in No. C fan. After detecting these cracks, the maintenance for No. B and No. C fans was carried out in March 2012 to prevent further accidents. These above results show the effectiveness and robustness of the proposed quantitative detection The proposed blade crack quantitative identification method is also applied to identifying the running condition of the remaining centrifugal booster fans. As mentioned previously, there are three fans in the unit 4. While No. A fan was being repaired, the quantitative identification index is also used to examine the No. B and No. C fans. These three fans have similar structures and functions, so the results based on the proposed method could be employed to analyze the condition of the No. B and No. C fan. After calculation and comparison, the values LLE > 0 and K 2f = 0.024 show that the No. C fan is in an incipient fault state. It should be focused on and chosen for monitoring. In addition, the fault status of the No. B fan (K 2f = 0.33) is much more serious. It needs to be repaired to avoid it breaking apart like No. A fan. During the checks, cracks were found in the No. B fan and No. C fan as expected. Figures 22 and 23 [24] show the results. From the Figure 22a, we can see that the crack grows from the fan entrance and propagates along the radial direction of the damper. Figure 22b shows the welding position of the blade and the damper. Stress concentration in this position led to the growth and propagation of the crack. As shown in Figure 23, a crack appears at the welding position with a length of 3-4 cm in No. C fan. After detecting these cracks, the maintenance for No. B and No. C fans was carried out in March 2012 to prevent further accidents. These above results show the effectiveness and robustness of the proposed quantitative detection method for blade crack faults.

Mathematical Modelling for Revealing Vibration Signal Properties
An impeller consists of a cover component, a disk component and several blades. The finite element model of an impeller is depicted in Figure 24. According to many experienced engineers, cracks initiate mostly at the weld toe on the cover sides of the blade, as shown in Figure 24. In this paper, the effects of a crack located at such a position on the vibration response of the impeller are of interest, and the cracks at other position can

Mathematical Modelling for Revealing Vibration Signal Properties
An impeller consists of a cover component, a disk component and several blades. The finite element model of an impeller is depicted in Figure 24. According to many experienced engineers, cracks initiate mostly at the weld toe on the cover sides of the blade, as shown in Figure 24. In this paper, the effects of a crack located at such a position on the vibration response of the impeller are of interest, and the cracks at other position can be modeled without much modification. For details of the mathematical modelling and mathematical formulation readers may refer to [34,37].

Mathematical Modelling for Revealing Vibration Signal Properties
An impeller consists of a cover component, a disk component and several blades. The finite element model of an impeller is depicted in Figure 24.

Mathematical Modelling for Revealing Vibration Signal Properties
An impeller consists of a cover component, a disk component and several blades. The finite element model of an impeller is depicted in Figure 24. According to many experienced engineers, cracks initiate mostly at the weld toe on the cover sides of the blade, as shown in Figure 24. In this paper, the effects of a crack located at such a position on the vibration response of the impeller are of interest, and the cracks at other position can be modeled without much modification. For details of the mathematical modelling and mathematical formulation readers may refer to [34,37]. According to many experienced engineers, cracks initiate mostly at the weld toe on the cover sides of the blade, as shown in Figure 24. In this paper, the effects of a crack located at such a position on the vibration response of the impeller are of interest, and the cracks at other position can be modeled without much modification. For details of the mathematical modelling and mathematical formulation readers may refer to [34,37].
After component synthesis, the equation of motion of a complete impeller is represented by: where M, C and K are the mass, damping and stiffness matrices of the ROM; p is the vector of the DOFs retained; b denotes the external excitation acting on the impeller; and f nl (p) is the nonlinear forces caused by intermittent contact of the crack surfaces. More details of this modelling process can be found in [34]. One of the most important issues for crack detection and identification is to search for sensitive indicators. An effective indicator should possess several features, such as robustness, monotonicity and industrial testability. The resonant frequencies discussed in the previous sections are potential indicators for the quantitative detection of crack faults. Two other kinds of frequency based indicators are studied in [34].
In sum, the frequency-based indicators for crack identification of centrifugal impellers were studied and discussed. However, an effective and reliable tool with a sensitive indicator for crack identification of impellers in operation still faces a lot of challenges at present.

Conclusions
In this paper, a vibration analysis method for the purpose of detection and quantitative identification of blade crack faults based on the amplitude of the rotating frequency is proposed. Aiming at the problem of energy leakage in FFT, a discrete spectrum interpolation method is proposed to extract the amplitude and phase accurately first. Then a quantitative identification and abnormality alarm strategy based on compound indexes including the Largest Lyapunov Exponent and relative energy ratio of the second harmonic frequency component is proposed. The results show that the proposed method is feasible. In the future, more effective signal processing methods should be studied and used to extract the characteristics of blade crack faults. Moreover, dynamic modeling and analysis of cracked rotor blades is necessary and urgent in future work. More reasonable effective indexes could be constructed to indicate the crack initiation and propagation from the point of view of dynamic analysis. Furthermore, although the proposed method shows good performance, more reasonable parameter selection for the terms a and b in K 2f should be studied based on plenty of running condition sensor-dependent vibration data in the future, and urgent demands, including quantitative diagnosis and fault location techniques, still remain to be established for scheduling reasonable maintenance plans.