Abstract

Monitoring the vibrations of high-speed rotating blades is significant to the security of turbomachineries. Blade tip timing (BTT) is considered as a promising technique for detecting blade vibrations without contact online. However, extracting blade vibration characteristics accurately from undersampled BTT signals measured at varying rotational speed (VRS) has become a big challenge. The existing two methods for this issue are restricted within the order bandwidth limitation and require prior information and precise sensor installation angles, which is often unpractical. To overcome these difficulties, a compressed sensing-based order analysis (CSOA) method was proposed. Its feasibility comes from the sparsity of BTT vibration signals in the order domain. The mathematical model for the proposed method was built, and the optimizing principles for sensor number and sensor arrangement were given. Simulated and experimental results verified the feasibility and advantages of the proposed method that it could extract order spectrum accurately from BTT vibration signals measured at VRS without the drawbacks in the existing two methods.

1. Introduction

High-speed rotating blades are critical components in turbomachinery, such as gas turbine engines and aircraft engines [1]. They often suffer from multifrequency vibrations resulted from various exciting forces, including centrifugal force, aerodynamic force, and impact force. These vibrations reduce the service life of the blades and may lead to cracks and even fractures, resulting in catastrophic accidents [2]. Thus, it is quite necessary to monitor blade conditions online to find the damages in time. The traditional approach is to install strain gauges on the blades, which is rather costly because it requires high-quality telemetry units or slip rings and complex installation process [3, 4].

BTT technique is an alternative to strain gauges and has become a hot issue due to its advantages of noncontact measurement, easy installation, and the ability to monitor all the blades simultaneously [5]. The BTT system precisely measures the time of arrival (TOA) of each blade tip at each sensor which is installed in the casing around the bladed disk [6]. A reference sensor is mounted near the spindle to provide a reference time once per revolution, from which the TOA of nonvibrating blade can be derived [7]. The vibration displacement can be obtained via the difference in the TOAs of the vibrating blade and nonvibrating blade through which blade parameters including vibrating frequency and mode can be analyzed [8]. However, the BTT vibration signal is undersampled because the overall sampling frequency, which is determined by rotational speed and limited sensor number, is usually less than twice the maximum blade vibrating frequency [911].

Until now, some methods have been proposed to extract blade vibration characteristics from undersampled BTT signals. Witos and Wisnioch [12] in the Polish Air Force have been studying the BTT technique since 1993 and conducted excellent researches. They extracted TOA signal components and diagnostic symptoms by a numerical process. Moreover, they employed a phase mapping method for analysis. Beauseroy and Lengellé [13] proposed a new analytical method for multifrequency BTT signals with a group of regularly spaced optical sensors. Salhi et al. [14] introduced a subspace modelling method to identify the natural frequency of blades using two sensors installed at the interval of 180°. Lin et al. [15] presented a reconstruction algorithm for multimode blade vibration signals based on the sparse representation theorem. Hu et al. [9] reconstructed BTT signals measured by nonuniformly placed sensors by considering it as a bandpass signal. Pan et al. [16] proposed a method to extract the spectrum of multifrequency blade vibrations and alleviate the measurement uncertainties of BTT signals. These methods are all based on time domain and constant speed assumption. However, the rotational speed may change frequently and result in inconstant sampling frequency, which invalidates these methods. Besides, the methods in [1316] require strictly precise sensor installation angles, which are difficult to realize.

Although extracting blade vibration characteristics from BTT signals measured at VRS has become an obstacle, few methods were proposed for solving this problem. Zhan et al. [17] proposed an interpolation algorithm to reconstruct blade vibration signals measured by uniformly arranged sensors, which made compensation for the effects of VRS. However, this method required strictly precise sensor installation angles and prior information, including central frequency and signal bandwidth. Chen et al. [18] improved the work in [17] and proposed two-order analysis-based reconstruction methods for BTT vibration signals at VRS by considering the signals as bandpass. Compared to [17], the first method in [18] substituted angular domain for time domain and still required uniformly installed sensors. The second method accurately captured high-engine orders of blade vibrations by using two sensors with arbitrary installation angles. However, neither of them can reconstruct vibration signals without the prior information of the central engine order and order bandwidth of the signal. Besides, the order bandwidth is quite limited by a sensor number.

Compressed sensing (CS) theory was proposed to accurately reconstruct spectrum-blind (no prior information) sparse signals from undersampled samples [1922]. Since the BTT signals measured at VRS have sparsity in the order domain, this paper proposed an order analysis method to obtain the order spectrum of BTT vibration signals at VRS based on CS theory. Compared to the existing methods for BTT vibration signals at VRS in [17, 18], the proposed method has three major advantages: firstly, it requires no prior information of the central engine order and order bandwidth of signals; secondly, it is still feasible to the signals exceeding the order bandwidth limitation as stated in [18]; and thirdly, it allows considerable installation angle errors.

This paper is organized as follows: the BTT mechanism was elaborated in Section 2. A CSOA mathematical model was built in Section 3. Then, the optimizing principles for the number and arrangement of the sensors were provided in Section 4. The existing methods in [18] were described in detail in Section 5. Simulations in Section 6 and experiments in Section 7 were conducted to verify the feasibility of the proposed method and the advantages over the existing methods. Finally, conclusions were drawn in Section 8.

2. BTT Mechanism

As shown in Figure 1, several sensors are mounted in the stationary casing around a rotating bladed disk to measure the TOAs of each blade at each sensor. Sensor r is installed close to the spindle with a reflective tape to provide a reference time for each revolution. The time pulse signals are recorded for the first time when the reflective tape passes Sensor r. The TOAs and the reference time are usually acquired by extracting the rising edge time or falling edge time from the recorded time pulse signals.

Considering Sensor r as the reference position, the angular positions of the ith sensor (Sensor i) and the bth blade (Blade b) at the time when the reflective tape passes Sensor r are represented as and , respectively. At the nth revolution, the actual TOA of Blade b at Sensor i and the duration of rotation are denoted as and , respectively. If the blades do not vibrate, the TOA of Blade b at Sensor i during the nth revolution is termed as the expected TOA, , which can be derived bywhere is the average rotational frequency from the time the nth revolution starts to the expected TOA, which can be derived by the spline interpolation method.

Practically, the blades usually vibrate under various incentives, causing a deflection when arriving at the sensors, which results in a time difference between the actual TOA and expected TOA. The deflection of Blade b when it arrives at Sensor i during the nth revolution can be obtained bywhere denotes the distance from the axis to the tip of Blade b. Thus, by matching the sequences and , a time-based vibrating displacement sequence of Blade b arriving at Sensor i can be acquired. Through analyzing all the sequences, the vibrating characteristics of blades can be obtained, thus providing information for detecting blade damages in real time.

3. CSOA Method

Consider a real continuous angular-domain blade vibration signal for Blade b, and m sensors are installed in the casing. Here, the expected rotating angles (no blade vibrations) at which the blades pass the sensors are considered as angular sampling points. Since the installation angles remain unchanged during the measurement process, the sampling intervals of a single sensor on the angular domain are identical. A sampling vibrating signal of Blade b measured by Sensor can be described as follows:

Let , then equation (3) can be rewritten as

The sampling function can be expressed as

Since is a periodic function with the sampling interval, , it has a Fourier expansion:where is the Fourier coefficient, which can be given by

Let . Its Fourier transform can be calculated aswhere E is the engine order. The discrete Fourier transform (DFT) of is given by

Equation (9) can be tied to equation (8) since they are transforms of the identical signal, which can be formulated aswhere is set to avoid the aliasing effect according to Shannon’s sampling theorem since the vibrating displacement of Blade b is sampled once per revolution by Sensor i. is an integer chosen to cover all the nonzero terms in order spectrum, which can be obtained bywhere is the maximum order of the signal. It can be seen from equation (10) that can be considered as a linear combination of shifted copies of . Thus, the CSOA mathematical model can be derived as follows by rewriting equation (10) in matrix form:where

The sensing matrix is constituted by

Since the BTT signal is sparse in order domain, the analysis of the signal can be transformed into a typical CS problem as follows:

Figure 2 provides a schematic of the proposed CS mathematical model. The abscissa and the ordinate of the graph, respectively, denote and with varying k. is set to cover the order scale of the blade vibration signal . It is considered the blade vibration signal has several orders arbitrarily distributed in the order domain. It can be seen that the order spectrum of the blade vibration signal can be extracted by finding the nonzero terms in the order scale, which can be transformed into the typical CS problem shown in equation (18) and solved by existing CS algorithms.

Since the CS algorithms can find all the nonzero terms in ignorance of their distribution, the proposed method requires no prior information of central frequency and order bandwidth of the blade vibration signal. In addition, the key to solve this CS problem is constructing the sensing matrix , which comes down to the choices of m, L, and .

4. Parameter Optimization

It can be concluded from Section 3 that the parameter optimization for the sensing matrix contributes greatly to successfully extract the order spectrum from the blade vibration signal. L can be chosen via equation (15), while the optimizations of m and require further investigations. This is an optimizing problem of the number and arrangement of the sensors.

4.1. Sensor Number Optimization

According to the CS theory [19, 20], is the unique sparsest solution of equation (18) and can be accurately extracted with a great possibility if it satisfies in the following equation:where is the smallest number of columns from which are linearly dependent [2325]. According to this definition, . Thus, the sensor number can be related to as follows:

However, equation (20) is a sufficient but not necessary condition for the unique sparsest solution of equation (18) and great possibility of recovering . Thus, sensor numbers which are smaller than may also be allowable.

4.2. Sensor Arrangement Optimization

According to [26, 27], can be accurately recovered by orthogonal matching pursuit (OMP) [28, 29] or basis pursuit (BP) [30, 31] algorithms if it satisfies the following condition [25]:where is the two-sided coherence [32] of , which is given bywhere is the ith column vector of . It can be inferred from equation (21) that the smaller the is, the larger the is allowed, which in fact allows a greater number of the blade vibration engine orders. Thus, the sensor arrangement principle is to choose proper installation angles to minimize as much as possible.

5. Detailed Description of the Existing Method

The existing two methods presented by Chen [18] considered the blade vibration signals as bandpass based on Shannon’s sampling theorem in the angular domain. The first method uses several uniformly mounted sensors to obtain uniformly sampling signals in the angular domain with a constant sampling interval . The first method can be expressed aswhere is the reconstructed signal and is the Hilbert transform of the sampled blade vibration signal. is the central order, and denotes the sampling order. However, the order bandwidth of the signal cannot exceed in equation (23). Besides, this method requires strictly precise installation angles to ensure uniform sampling in the angular domain.

The second method uses two nonuniformly mounted sensors. It is based on the second-order sampling theorem [33]. Its reconstruction formula is given bywhere and are the blade vibration signals sampled by the two sensors, respectively. is the angle between the two sensors. is the order bandwidth, and can be given bywhere . The second method does not require precise installation angles. However, it cannot reconstruct blade vibration signal whose order bandwidth exceeds one.

6. Numerical Simulations

6.1. Simulated Blade Vibration Signal

It is assumed that each blade is identical and is not coupled with others. In addition, they vibrate in a single degree of freedom under the simulated synchronous exciting force, which can be given bywhere and are the amplitude and the order of the force, respectively. is the rotating angle, which is calculated as follows:where is the rotational frequency varying with time.

A mathematical model for blade vibration can be derived aswhere , c, and k are the equivalent mass, damping coefficient, and elastic coefficient of the blade. x is the blade tip displacement.

Based on equations (26)–(28), a simulated signal acquisition system for blade tip vibration at VRS was built, as shown in Figure 3. Chirps 1∼3 generate swept frequency cosine with instantaneous frequencies equal to , , and times the rotational frequency, respectively. Chirp 4 provides swept frequency sine with an instantaneous frequency equal to the rotational frequency. The switches enable the system to output the blade tip displacements each time it satisfies the following condition:where the angular position of the blade and the ith sensor installation angle is set in hit crossings. Port 1 outputs the original blade tip displacements, Port 2∼j outputs the displacements measured by the sensors, and Port j + 1 is used for calculating the starting time sequences of each revolution.

6.2. Simulation for the Order Bandwidth Limitation

To verify that the CSOA method is still feasible to the BTT vibration signals which exceed the order bandwidth limitation of the existing methods in [18], simulated signals of different generated by the system in Figure 3 were analyzed. The simulation parameters are shown in Table 1 and have been optimized according to the theoretical analysis in Sections 3 and 4. In addition, white Gaussian noise is added to the simulated signals.

The natural frequency of the blade can be calculated by

Figure 4 shows the simulated blade tip displacements varying with the rotating frequency and the data sampled by the sensors. The blade vibrates resonantly when the excitation frequency reaches .

Figure 5 shows the order spectrum obtained via the CSOA method. It is obvious that the order spectrums of the two sets of simulated blade tip displacements are accurately obtained.

For comparison, the data were also analyzed with the existing methods. The sensor number for the first method and for the second method in [18] are set to 7 and , respectively. The results are shown in Figure 6.

The simulation results suggested that the two methods in [18] failed to acquire the exact order spectrum of the signal beyond the bandwidth limitation. In conclusion, the CSOA method is effective whether the signal is within or exceeding the order bandwidth limitation of the existing methods proposed in [18]. Moreover, L in CSOA can be enlarged to meet the needs of a larger order bandwidth.

6.3. Simulation for Installation Angle Errors

To demonstrate that the CSOA method allows considerable installation angle errors, the angles were set as , where is the deviation from the expected installation angles and randomly generated within . Other parameters are shown in Table 1. The simulation results via the CSOA method are shown in Figure 7.

In Figure 7, the order spectrum is precisely acquired despite the installation angle errors, which proves that the CSOA method is still effective under the condition of inaccurate installation angles. This characteristic is quite useful since the installation angle errors exist in practical applications.

7. Experiments

For further verification on the effects of the proposed CSOA method, an experimental rig was built, as shown in Figure 8. The rig consists of an electrical motor (installed under the test bench), a bladed disk with 32 uniformly distributed straight blades (labelled as Blades 1∼32), a spindle, a casing, six uniformly distributed permanent magnets, and four optical fiber sensors (Sensors 1∼3 were mounted on the casing in precise alignment along the circumstance, and Sensor r was the reference sensor). The operating frequency of a high-speed data acquisition card is 5 MHz. Sensor r was installed vertically, and the installation angles of Sensors 1∼3 were set as , , and , respectively. These angles were corrected into , , and using the TOAs measured by the sensors at a constant speed. The parameters of the blades are presented in Table 2.

Finite element simulation analysis of the blade was conducted according to the parameters in Table 2. Figure 9 shows that the modal shapes of the first three vibration modes were bending, bending, and torsion, respectively. And Figure 10 shows the Campbell diagram of the blade model. Point A (4397, 2199) was selected, where the blade resonated at the first mode and 30 times the rotating frequency. Thus, the rotating speed was varied around point A to obtain blade resonant vibration.

According to the modal analysis result of the blade, the rotating speed was increased continuously from 3850 rpm to 4850 rpm. Blade 27 was arbitrarily chosen for analysis. The blade tip displacement varied with the rotating speed and maximized at the resonant speed range (RSR), 4422∼4426 rpm (Figure 11(a)), which was close to the modal analysis result. Figure 11(b) shows that the blade vibrated at 30 times the rotating frequency. The consistency of the modal analysis and experimental result provided a strong evidence for the validity of the CSOA method.

The experimental data of all the blades were further analyzed with the CSOA method and the two methods in [18], as shown in Table 3. represents the vibration order obtained by the CSOA method. and are the vibration orders identified, respectively, by the existing two methods with the presupposed central order (PCO) equal to . The RSRs were between 4387 rpm and 4448 rpm, which suggested that the blades all vibrated at 30 times of the rotating frequency. varied from 29.97 to 30.06, which further validated the proposed CSOA method. When the PCO was set as 30, both existing methods identified the vibration orders precisely. However, the changes in PCO lead to incorrect results, which demonstrated that the success of these methods depended on the PCO. Thus, the CSOA method has an advantage over the existing two methods when the central order of blade vibration is unknown.

8. Conclusions

Currently, BTT is considered as a promising technique for online monitoring blade vibrations at a high-rotational speed. However, VRS brings obstacles to extracting blade vibration characteristics from undersampled BTT signals. To solve this problem, this paper proposed a CSOA method for BTT signals measured at VRS, which has good advantages over the existing methods. A CSOA mathematical model was built. Simulations and experiments were performed to demonstrate the advantages and validity of the CSOA method. The main results are summarized as follows:(1)The order spectrum of undersampled BTT vibration signals measured at VRS can be accurately obtained through the proposed CSOA method.(2)No prior information of the central order and order bandwidth of BTT signals is needed for analysis via the CSOA method.(3)The CSOA method is effective whether the BTT vibration signals are within or exceeding the order bandwidth limitation of the existing methods in [18].(4)The CSOA method allows considerable installation angle errors, which is more applicable to the actual situation.

Data Availability

The data used to support this research are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the National Basic Research Program of China (Grant no. 2015CB057404) and the Natural Science Foundation of Hunan Province, China (Grant no. 2019JJ50721).