Next Article in Journal
IMU-Based Movement Trajectory Heatmaps for Human Activity Recognition
Previous Article in Journal
Integration of a Thermoelectric Heating Unit with Ionic Wind-Induced Droplet Centrifugation Chip to Develop Miniaturized Concentration Device for Rapid Determination of Salmonella on Food Samples Using Antibody-Functionalized SERS Tags
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Parameter Measurement of Rotors Using the Doppler Effect of Frequency-Swept Interferometry

1
Key Laboratory of Optoelectronic Technology and Systems, Ministry of Education, Chongqing University, Chongqing 400044, China
2
Artificial Intelligence Center, Peng Cheng Laboratory, Shenzhen 518000, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(24), 7178; https://doi.org/10.3390/s20247178
Submission received: 29 October 2020 / Revised: 9 December 2020 / Accepted: 14 December 2020 / Published: 15 December 2020
(This article belongs to the Section Optical Sensors)

Abstract

:
The Doppler effect of frequency-swept interferometry (FSI) is often seen as an obstacle to the dynamic ranging accuracy. However, the potential of this obstacle is rarely noticed and used. In this paper, by combining the periodical characteristics of the rotational Doppler effect, an FSI-based multi-parameter measurement method for the rotor is proposed. Through the establishment of the rotational Doppler formula of FSI, it is found that the frequency, direct component, and amplitude of the dynamic distance given by FSI can be utilized to estimate the angular velocity, axial clearance, and tilt angle of the rotor, respectively. A rotor platform and a fiber-optic FSI system were constructed, and a series of experiments were carried out to verify the proposed method. The experimental results showed that the relative errors of the measured axial clearance, angular velocity, and tilt angle were less than 3.5%. This work provides a new perspective on the multi-parameter measurement of the rotor and makes it possible to directly perform multi-parameter measurement inside the space-confined rotating machinery as only a single small-size fiber-optic probe is needed.

1. Introduction

Rotors are the core component of various rotating machines, such as motors, generators, air compressors, wind turbines, and aero engines. Although these machines are designed for different purposes, their running efficiency is mainly determined by the rotor parameters [1,2]. Therefore, the measurement of rotor parameters is of great importance for performance evaluation [3,4], fault diagnostics [5,6], health monitoring [7,8,9,10], fatigue life prediction [11,12], and dynamic control [13,14] of rotating machines.
The ideal running state of the rotor is shown in Figure 1a, where there is no contact-rubbing between the rotor and the stationary part. However, the real situation, as shown in Figure 1b, is that the rotor will deviate from the optimal operating state due to the manufacturing error, assembly tolerance, imbalance, and other complex mechanisms. Such deviation may result in different damages ranging from cosmetic damage to full destruction of the machine. To know and master the actual working state of the rotor, the synchronous measurement of the rotational speed Ω, tilt angle of the rotor α, and clearances L and D between the rotor and stationary part is necessary. To realize the measurement of these parameters, the conventional way is to integrate a number of individual subsystems, and each subsystem is dedicated to measuring a single particular parameter. Because each subsystem has its own sensor and corresponding demodulation equipment, the integrated system has a high cost and low overall reliability. Furthermore, the multi-sensor configuration requires sufficient space to allow the access and installation of multiple sensors, in practice, which is a large difficulty for space-confined rotating machines. Besides, as a side effect, the load or stress imposed by multiple sensors may, in turn, affect the running status of machines. Thus, to overcome the above challenges and better serve various applications, it is highly desired to develop a multi-parameter measurement method with fewer sensors. At present, a wide range of techniques have been reported on for the measurement of radial clearance D [15,16,17,18,19,20,21,22]. Among these techniques, the fiber-optic-based sensing approaches [17,18,19,20,21,22] are especially favored because of the intrinsic capabilities of optic fiber, such as its light weight, small size, and adaptability for embedding into structures. However, due to the narrow internal space and harsh measurement environment of the rotating machines, the measurement of the tilt angle and axial clearance under high-speed rotation is still a difficult problem that remains unsolved.
As a famous absolute distance measurement technique, frequency-swept interferometry (FSI) is widely used in various ranging applications. Since the dynamic accuracy of FSI is severely restricted by the well-known Doppler effect [23,24,25], the previously reported works mainly focused on the elimination of the Doppler error of FSI [23,24,25,26,27,28,29,30,31,32]. Because the goal of these methods was to measure the absolute dynamic distance, they were not applied for the measurement of the rotor parameters. In our previous work [32], we proposed an FSI-based method which is suitable for dynamic clearance measurement, but the core idea is still to remove the Doppler error in FSI. For absolute distance measurement, the Doppler effect is an undesirable property that needs to be eliminated, but we found, it is exactly this property that makes the multi-parameter measurement of the rotor possible. If the probe of the FSI is pointed to the rotor disk, the parameters L, Ω, and α become three factors that modulate the frequency, direct component, and amplitude of the dynamic distance given by FSI. Thus, by reversely using the characteristics of the dynamic distance, the rotor parameters L, Ω, and α can be simultaneously obtained according to the modulation relationship. Different from the thought of Doppler effect elimination, here we take the Doppler effect as a tool and utilize it to address the unresolved issues in rotor measurement. To realize the multi-parameter measurement, we first derived the modulation relationship between the Doppler effect and the parameters L, Ω, and α. Based on the modulation relationship, the estimation method of L, Ω, and α was given. In addition, experiments were conducted to prove the effectiveness of the estimation method. By using this method, the L, Ω, and α of rotors could be determined with a simple and low-cost FSI system, and as this system requires a single fiber-optic sensing probe, it has a great advantage for the multi-parameter measurement inside the space-confined rotating machines.
The remaining sections of this paper are organized as follows: In Section 2.1, the problem of the Doppler effect of FSI is reviewed. In Section 2.2, we model the Doppler effect of the rotor and propose a method to estimate the rotor parameters L, Ω, and α. Also, we analyze the impacts of the deflection angles caused by the non-parallel installation of the probe on the multi-parameter measurement error and give guidance for the probe installation location. Section 3.1 describes the experimental setup, and Section 3.2 presents the experimental procedure to validate the multi-parameter measurement method. Finally, concluding remarks are given in Section 4. In addition, Appendix A demonstrates the derivation process of the upper bound of the relative error in the estimation of α.

2. Theory

2.1. The Basics of FSI

A fiber-optic FSI system is shown in Figure 2; the light wave of the frequency-swept laser (FSL) passes through the fiber-optic circulator (FOC) and reaches the end face of the fiber probe. After the reflections of fiber probe and target, the back-propagating waves (i.e., the reference wave with frequency fr(t), and the measurement wave with frequency fm(t)) interfere and generate the FSI signal in the single-mode fiber (SMF), as shown in the enlarged figure. Then the FSI signal is routed to the photodetector (PD) by the FOC and finally sampled by the data acquisition system (DAQ). For a static target at distance L0, the FSI signal s(t) can be expressed as [30]
s ( t ) = cos { 2 π 0 t [ f r ( t ) f m ( t ) ] d t } = cos { 2 π 0 t [ f r ( 0 ) + B T t f r ( 0 ) B T ( t 2 L 0 c ) ] d t } = cos [ 2 π ( B T 2 L 0 c ) t ] = cos [ ψ ( t ) ] ,
where fr(t) = fr(0) + (B/T)t, fm(t) = fr(0) + (B/T)(t − 2L0/c), B is the sweep bandwidth of FSL, T is the sweep period of FSL, and c is the speed of light in air. Thus, using the phase ψ(t) of the measured s(t), the L0 can be determined as ((t)/dtTc/(4πB).
However, if the target moves during the frequency sweep, due to the Doppler effect, the phase ψ(t) in Equation (1) becomes [24] (note that, Jia et al. [24] gives the instantaneous frequency of the FSI signal, which equals (t)/(2πdt))
ψ ( t ) = 4 π c o t [ B T L ( t ) + f r ( t ) v ( t ) ] d t ,
where L(t) and v(t) are the real-time distance and velocity (v(t) is negative when the target approaches the fiber probe). Thus, using the (t)/dt, the measured dynamic distance Lm(t) can be obtained as [32]
L m ( t ) = T c 2 B d ψ ( t ) 2 π d t = L ( t ) + f r ( t ) B v ( t ) T L ( t ) + f a v g B v ( t ) T ,
where favg is the average of fr(t). In Equation (3), the term favgv(t)T/B is the Doppler error, which is sensitive to the target velocity v(t). This error severely restricts the dynamic ranging accuracy and needs to be eliminated. Whereas from the opposite perspective, the Lm(t) may be utilized since it contains both the information of L(t) and v(t).

2.2. Multi-Parameter Measurement Method for the Rotor

The general case of the multi-parameter measurement is shown in Figure 3. We assumed that the center of the rotor locates at the origin o (0, 0, 0) of the Cartesian coordinate system xyz, and that the rotor rotates around the z-axis with angular velocity Ω. Due to the tilt angle α, the normal vector of the rotor on = [ sin α cos ( Ω t ) ,   sin α sin ( Ω t ) ,   cos α ] will rotate around the z-axis with angular velocity Ω and angle α. We assumed the fiber probe is installed at an off-axis position o’ (r, 0, L), and the direction of the laser beam is o a   =   [ sin β cos φ ,   sin β sin φ ,   cos β ] , where β and φ are the deflection angles between the laser beam and the coordinate system xyz’ (parallel to xyz). Let the coordinates of point A be (xA, yA, zA), then we have oA   =   [ x A ,   y A ,   z A ] and o A   =   [ x A r ,   y A ,   z A L ] . Thus, according to the relations of on oA ( on · oA   =   0 ) and o a | | o A ( sin β cos φ / ( x A     r )   =   sin β sin φ / y A   =   cos β / ( z A L ) ), the following equations involving the angular velocity Ω, axial clearance L, and tilt angle α can be established.
[ sin α cos ( Ω t ) sin α sin ( Ω t ) cos α 1 0 tan β cos φ 0 1 tan β sin φ ] [ x A y A z A ] = [ 0 r L tan β cos φ L tan β sin φ ] .
Note that, the first line of Equation (4) is the result of on · oA   =   0 , the second and third lines are the equivalent transformations of sin β cos φ / ( x A     r )   =   sin β sin φ / y A   =   cos β / ( z A     L ).
According to Equation (4), zA can be calculated as:
z A = r sin α cos ( Ω t ) L sin α tan β cos ( Ω t φ ) cos α sin α tan β cos ( Ω t φ ) .
Using the angle β between the o a and the o z , the real dynamic clearance L(t) can be calculated as:
L ( t ) = | o A | = sec β ( L z A ) .
Then, substituting L(t) and v(t) = ∂L(t)/∂t into Equation (3), the dynamic clearance given by the FSI can be expressed as
L R o ( t ) = L ( t ) + f a v g T B L ( t ) t = L ( t ) f a v g T B ( z A sec β ) t = L + r tan α 1 + ( f a v g T Ω B ) 2 s i n ( Ω t + ξ ) f a v g T Ω B ( L z A ) tan α tan β sin ( Ω t φ ) cos β [ 1 tan α tan β sin ( Ω t φ ) ] ,
where subscript Ro denotes rotor, and:
ξ = arctan ( B f a v g T Ω ) .
From Equation (7), it can be found that the probe installation parameters r, β and φ, and the rotor parameters L, α, and Ω, are involved in LRo(t). For a parallel laser beam, i.e., β ≈ 0, Equation (7) can be reduced as:
L R o ( t ) = L + r tan α 1 + ( f a v g T Ω B ) 2 A R o s i n ( Ω t + ξ ) .
Equation (9) shows that the LRo(t) varies sinusoidally with angular velocity Ω and amplitude ARo. Thus, by using the frequency, periodicity, and amplitude characteristics of LRo(t), the rotating speed Ω, the axial clearance L, and the tilt angle α of the rotor can be estimated by the following equations:
Ω ^ = F T [ L R o ( t ) ] ,
L ^ = Ω ^ 2 π m t π m / Ω ^ t + π m / Ω ^ L R o ( t ) d t ,
α ^ = arctan [ A R o r ( 1 + ( f a v g B T Ω ^ ) 2 ) 1 ] ,
where the symbol ^ denotes the estimated value of the parameter, FT represents solving the angular frequency of LRo(t) by using frequency analysis method such as Fourier Transform, m is a positive integer (in this paper we set m = 1 for following calculations), and r in Equation (12) is a pre-determined value during the probe installation. Note that the LRo(t) can be obtained by ((t)/dtTc/(4πB), where ψ(t) is the phase of the measured FSI signal s(t) and which can be given by the Hilbert transform [26].
In practice, the perfect parallel installation of the probe is difficult, i.e., β is small but nonzero. To analyze the influence of β on the multi-parameter measurement, we expand Equation (7) by using the first-order approximation at the point of tanβ = 0
L R o 1 ( tan β ) = n = 0 , 1 ( tan β ) n n ! [ n L R o ( tan β ) n ] | tan β = 0 = L + r tan α 1 + ( f a v g T Ω B ) 2 sin ( Ω t + ξ ) + K tan β ,
where:
K = tan α sin ( Ω t φ ) { [ L r tan α cos ( Ω t ) ] ( 1 f a v g T Ω B ) + f a v g T Ω B r tan α sin ( Ω t ) } .
It is easy to prove that the K has a period of 2π/Ω and a mean value of 0, which indicates that the angle β does not affect the values of Ω ^ and L ^ . Comparing with Equation (9), the term Ktanβ in Equation (13) can be regarded as a perturbation on the amplitude ARo, which will lead to an error on α ^ . We calculate the relative error e α ^ , which has an upper bound as follow (see Appendix A for details):
e α ^ tan β r ( L + 2 r tan α ) .
Equation (15) indicates that the relative error is related to the parameters β, r, L, α. To achieve a certain level of estimation accuracy, it is better to specify suitable probe installation parameters before measurement according to Equation (15). For instance, e α ^ < 5% requires (L/r + 2tan α ) < 1 for β = 3°, which means even for a large tilt angle α = 5°, the estimation accuracy can be easily guaranteed by installing the probe at a location where L/r < 0.8. Note that Equation (15) holds for arbitrary φ ∈ (0, 2π).

3. Experiments and Results

3.1. Experimental Setup

To verify the proposed method, an experimental system was built, as shown in Figure 4. The FSL (Arcadia Optronix, GC-760001c, frequency-swept range 196,250 GHz-191,386 GHz, sweeping period 0.486 ms, coherent length < 95 m.) was adopted as the laser source, and the FSI signal detected by the PD was acquired by the DAQ with a sampling rate of 5 MSa/s. The probe was a cleaved single-mode fiber (output laser beam divergence was about 8°). The position of the probe can be adjusted by the motorized stage (adjusting accuracy was better than 1 μm). The rotational motion of the aluminum rotor was provided by the motor. The driving voltage of the motor was generated by the signal generator (DG4102, RIGOL) and amplified by the voltage amplifier (HAS 4011, NF Corporation). The rotational speed was monitored by the Hall sensor installed at the end of the motor.

3.2. Results and Analysis

To verify the effectiveness of the proposed method, it was necessary to know the true values of the rotor parameters (i.e., the true values of L, α, and Ω) in the experimental system. For this purpose, before loading the driving voltage on the motor, we manually rotated the rotor three cycles and measured the actual clearance L(θ) = L + rtanαcos(θ) under the condition of Ω = 0 rad/s by using the fiber probe (r = 2 mm, β < 2°). The measured L(θ) is shown in Figure 5, and the true values of L and rtanα can be obtained as L = 1987.0 μm and rtanα = 47.1 μm. As r = 2 mm, the true value of the tilt angle α can be determined as arctan (rtanα/r) = 1.35°. Thus, L = 1987.0 μm, α = 1.35°, and as for Ω, it was determined by the Hall sensor when the rotor rotates.
According to Equation (9), after loading the driving voltage, there will be a sinusoidal modulation on the LRo(t) when the rotor rotates. To verify Equation (9), the angular velocity Ω was set to 60π rad/s, and the experimental and theoretical values are plotted in Figure 6. The measured LRo_meas(t) was obtained from the phase slope (t)/dt of the raw FSI signal, and the theoretical LRo_theo(t) was calculated according to Equation (9). It was found that the LRo_meas(t) matches well with the LRo_theo(t). Moreover, comparing Figure 5 and Figure 6, it was also found, due to the rotation-induced Doppler effect, the amplitude of LRo_meas(t) increases to 180.4 μm, which was 3.83 times that of L(θ), and this value agrees well with the theoretical amplitude modulation coefficient (1 + (favgTΩ/B)2)1/2 = 3.79 given by Equation (9).
According to the dynamic FSI model of rotors in Equation (9), the direct component (DC component), frequency, and amplitude of the LRo(t) are related to the axial distance L, angular velocity Ω, and tilt angle α. Therefore, when the LRo_meas(t) was obtained, the multiple parameters could be calculated according to Equations (10)–(12). For example, using the LRo_meas(t) in Figure 6, the L and Ω were respectively estimated to be 1987.0 μm and 60π rad/s. Besides, using ARo = 180.4 μm, the α was estimated to be 1.36°, according to Equation (12), which was consistent with the true value of 1.35°. The relative measurement error of α was 1%, less than the theoretical upper bound e α ^ = tan β ( L / r + 2 r tan α ) = 3.7 % given by Equation (15).
To further verify the adaptability of the proposed method at different Ω, L, and r, a series of experiments were carried out. First, we kept the position of the fiber probe unchanged and set Ω from 60π rad/s to 160π rad/s. The multi-parameter measurement results are shown in Figure 7. In Figure 7a, the theoretical Ω was given by the Hall sensor, and the estimated Ω ^ is calculated by the proposed method. It can be seen that the angular velocity was accurately recovered, with a relative error of about 1.5%. As expected, the estimated L ^ and α ^ in Figure 7b remained constant, and their relative errors were less than 0.04% and 2.53%, respectively.
Then, we reset the Ω to the fixed value of 60π rad/s and changed the L by using the motorized stage. The L was changed from 1987 μm to 2087 μm with a step of 20 μm. The measurement results are shown in Figure 8. In Figure 8a, the theoretical L is the sum of the displacement and the initial clearance 1987 μm, and the L ^ is obtained according to Equation (11). It was evident that the estimated L ^ was consistent with the theoretical L, and the maximum relative measurement error was less than 0.06%. It can be seen from Figure 8b that the estimated Ω ^ and α ^ were independent of the variation of L ^ . The relative measurement errors of Ω ^ and α ^ were respectively less than 0.91% and 1.02%.
Next, we fixed the L and Ω respectively to 1987 μm and 60π rad/s, and changed the fiber probe position r from 2 mm to 4 mm with a step of 0.5 mm. Corresponding measurement results under different positions are illustrated in Figure 9. The maximum estimation relative errors of α ^ , L ^ and Ω ^ , shown in Figure 9 were less than 3.42%, 0.04%, and 2.11%, respectively. The proposed method still performed well in this case.

4. Conclusions

In this paper, we presented a method for simultaneously measuring the axial clearance, rotational speed, and tilt angle of the rotor by using the Doppler effect of FSI. The measurement steps are as follows: first, obtained the instantaneous phase ψ(t) of FSI signal by using the Hilbert transform and calculate the dynamic distance LRo(t) = ((t)/dtTc/(4πB) according to Equation (3); then, estimated the angular velocity Ω, axial clearance L, and tilt angle α by using the frequency, direct component and amplitude of LRo(t) according to Equations (10)–(12). Considering the practical non-parallel installation of the probe, we analyzed the influence of the non-parallel installation on the measurement error and found that the measurement error could be eliminated by installing the probe at a proper location (L/r < 0.8). The multi-parameter measurement experiments were carried out to verify the proposed method, and the relative errors of the axial clearance, rotational speed, and tilt angle were respectively lower than 0.1%, 2.2%, and 3.5%. The main difference between the previously reported FSI-based works and this work is that we used the Doppler error as a tool rather than an obstacle. This idea reduced the bulkiness of conventional, integrated multiple parameter measurement systems. Besides, the optic-fiber probe of our system is small in size, light-weight, and bendable, which is very suitable for the space-confined rotatory machines, such as generators, motors, and gear-sets. Moreover, combining our method with the reported fiber-optic radial clearance measurement method would enable a more comprehensive multi-parameter measurement of rotors.

Author Contributions

Conceptualization, B.S., W.Z., W.C., and P.Z.; Formal analysis, B.S. and W.Z.; Investigation, B.S. and W.Z.; Data analysis, B.S.; Writing—original draft, B.S.; Writing—review and editing, W.Z., W.C., and P.Z.; Supervision, W.C., and W.Z.; Funding acquisition, W.Z. All authors have read and agreed to the published version of the manuscript.

Funding

National Natural Science Foundation of China (51805054); Technology Innovation Platform Project of Aero Engine Corporation of China (AECC) (SHYS-GXDPL-18); China Postdoctoral Science Foundation (2018M643405).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Upper Bound of the Relative Error in the Estimation of α

According to Equation (7), we can consider LRo as a function of tanβ:
L R o ( tan β ) = L + r tan α 1 + ( f a v g T Ω B ) 2 s i n ( Ω t + ξ ) f a v g T Ω B ( L z A ) tan α tan β sin ( Ω t φ ) [ 1 tan α tan β sin ( Ω t φ ) ] / 1 + tan 2 β ,
where:
z A = r sin α cos ( Ω t ) L sin α tan β cos ( Ω t φ ) cos α sin α tan β cos ( Ω t φ ) = r tan α cos ( Ω t ) L tan α tan β cos ( Ω t φ ) 1 tan α tan β cos ( Ω t φ ) .
For small angle β, LRo(tanβ) can be approximated by the first-order Taylor series expansion at tanβ = 0
L R o 1 ( tan β ) = n = 0 , 1 ( tan β ) n n ! [ n L R o ( tan β ) n ] | tan β = 0 = L + r tan α 1 + ( f a v g T Ω B ) 2 sin ( Ω t + ξ ) + L R 0 tan β | tan β = 0 tan β ,
we can obtain:
L R 0 tan β | tan β = 0 = tan α sin ( Ω t φ ) [ L ( 1 f a v g T Ω B ) + r tan α 1 + ( f a v g T Ω B ) 2 sin ( Ω t + ξ ) + f a v g T Ω B r tan α cos ( Ω t ) ] = tan α sin ( Ω t φ ) [ L ( 1 f a v g T Ω B ) r tan α cos ( Ω t ) + f a v g T Ω B r tan α sin ( Ω t ) + f a v g T Ω B r tan α cos ( Ω t ) ] = tan α sin ( Ω t φ ) { [ L r tan α cos ( Ω t ) ] ( 1 f a v g T Ω B ) + f a v g T Ω B r tan α sin ( Ω t ) } .
Thus, we have
L R o 1 ( tan β ) = L + r tan α 1 + ( f a v g T Ω B ) 2 sin ( Ω t + ξ ) L R o + tan α sin ( Ω t φ ) { [ L r tan α cos ( Ω t ) ] ( 1 f a v g T Ω B ) + f a v g T Ω B r tan α sin ( Ω t ) } K tan β ,
where it can be found that the first term is the distance LRo under the condition of β = 0. Thus, the term Ktanβ induced by a small β can be regarded as a time-varying perturbation on LRo, which would lead to a change in the amplitude ARo of LRo. Thus, the change of ARo is within the range from −|Ktanβ|max to |Ktanβ|max. Here, we define e as
e = tan α ^ tan α ^ A R o + | K tan β | max r ( 1 + ( f a v g B T Ω ^ ) 2 ) 1 A R o r ( 1 + ( f a v g B T Ω ^ ) 2 ) 1 = 1 + | K tan β | max A R o ,
where α ^ and α ^ are respectively the estimation of α under the conditions of β = 0 and β 0 .
Using |a + b| ≤ |a| + |b|, |ab| = |a||b|, | sin ( Ω t ) | 1 and | cos ( Ω t ) | 1 , we can obtain:
| K tan β | max A R o | tan α sin ( Ω t φ ) { [ L r tan α cos ( Ω t ) ] ( 1 f a v g T Ω B ) + f a v g T Ω B r tan α sin ( Ω t ) } tan β | r tan α 1 + ( f a v g T Ω B ) 2 tan β r { | [ L r tan α cos ( Ω t ) ] [ 1 f a v g T Ω B ] | + | f a v g T Ω B r tan α | | sin ( Ω t ) | } | sin ( Ω t φ ) | 1 + ( f a v g T Ω B ) 2 tan β r { | [ L r tan α cos ( Ω t ) ] | | 1 f a v g T Ω B | 1 + ( f a v g T Ω B ) 2 + f a v g T Ω B 1 + ( f a v g T Ω B ) 2 r tan α } tan β r { | L r tan α cos ( Ω t ) | + r tan α } tan β r { | L | + | r tan α | | cos ( Ω t ) | + r tan α } tan β r ( L + 2 r tan α ) .
It can be seen from Equation (A7), the sin(Ωtφ) term that contains the deflection angle φ is amplified and taken as 1, thus Equation (A7) holds for arbitrary φ ∈ [0, 2π). Therefore, whatever the φ is, the relative estimation error of α ^ caused by a small β satisfies:
e α ^ = α ^ α ^ α ^ = arctan ( e tan α ^ ) α ^ α ^ arctan { [ 1 + tan β r ( L + 2 r tan α ) ] 1 tan α ^ 0 } α ^ α ^
Mathematically, for the following function y(p, q),
y ( p , q ) = p q arctan ( p tan q ) ,
we have
y q = p p sec 2 ( q ) p 2 tan 2 ( q ) + 1 = p 3 tan 2 ( q ) p tan 2 ( q ) + p + [ p tan 2 ( q ) p sec 2 ( q ) ] p p 2 tan 2 ( q ) + 1 = p ( p 2 1 ) tan 2 ( q ) p 2 tan 2 ( q ) + 1 ,
obviously, when p ≥ 1 and q ≥ 0, we have ∂y/∂q ≥ 0. So, we have
y ( p , q ) = p q arctan ( p tan q ) y ( p , 0 ) = 0 .
Therefore, in Equation (A8), if we define
{ p = [ 1 + tan β r ( L + 2 r tan α ) ] , q = α ^ ,
according to Equation (A11), Equation (A8) can be simplified as:
e α ^ arctan { [ 1 + tan β r ( L + 2 r tan α ) ] tan α ^ } α ^ α ^ α ^ [ 1 + tan β r ( L + 2 r tan α ) ] α ^ α ^ = tan β r ( L + 2 r tan α ) .
It needs to point out that, from Equation (A6) to Equation (A13), we consider the case of ARo increased by |Ktanβ|max. For the case of ARo decreased by −|Ktanβ|max, it is easy to get e α ^ tan β ( L + 2 r tan α ) / r . So the e α ^ in Equation (A13) can be regarded as the absolute value of the relative error.

References

  1. Hu, G.; Gao, W.; Chen, Y.; Zhang, D.; Tian, Y.; Qi, X.; Zhang, H. An experimental study on the rotational accuracy of variable preload spindle-bearing system. Adv. Mech. Eng. 2018, 10, 1–14. [Google Scholar] [CrossRef]
  2. Agnes, M. Vibrational diagnostics of rotating machinery malfunctions. Int. J. Rotating Mach. 1995, 1, 237–266. [Google Scholar]
  3. Vafaei, S.; Rahnejat, H. Indicated repeatable runout with wavelet decomposition (IRR-WD) for effective determination of bearing-induced vibration. J. Sound Vib. 2003, 260, 67–82. [Google Scholar] [CrossRef]
  4. Wefky, M.; Khalik, A.S.A.; Ahmed, S.; Arabawy, I.F.E. Performance evaluation of five-phase outer-rotor Permanent magnet vernier machines. In Proceedings of the Conference of the IEEE Industrial Electronics Society, Florence, Italy, 23–26 October 2016. [Google Scholar]
  5. Glowacz, A. Diagnostics of rotor damages of three-phase induction motors using acoustic signals and SMOFS-20-EXPANDED. Arch. Acoust. 2016, 41, 507–515. [Google Scholar] [CrossRef]
  6. Qu, L.; Lin, J.; Liao, Y.; Zhao, M. Changes in rotor response characteristics based diagnostic method and its application to identification of misalignment. Measurement 2019, 138, 91–105. [Google Scholar] [CrossRef]
  7. Abdul-Aziz, A.; Woike, M. Turbine rotor disk health monitoring assessment based on sensor technology and spin tests data. Sci. World J. 2013, 2013, 413587. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Abdul-Aziz, A.; Woike, M.R.; Oza, N.C.; Matthews, B.L.; Lekki, J.D. Rotor health monitoring combining spin tests and data-driven anomaly detection methods. Struct. Health Monit. 2011, 11, 3–12. [Google Scholar] [CrossRef]
  9. Gyekenyesi, A.L.; Baaklini, G.Y.; Baaklini, G.Y.; Boltz, E.S.; Shepard, S.M.; Shull, P.J. Rotor health monitoring and damage detection utilizing a disk spin simulation system. In Proceedings of the 6th Annual International Symposium on NDE for Health Monitoring and Diagnostics, Newport Beach, CA, USA, 4–8 March 2001. [Google Scholar]
  10. Matikas, T.E.; Abdul-Aziz, A.; Woike, M.; Baaklini, G.; Bodis, J.R. Turbine engine disk rotor health monitoring assessment using spin tests data. In Proceedings of the International Society for Optical Engineering, San Diego, CA, USA, 12–14 March 2012. [Google Scholar]
  11. Jain, J.R.; Kundra, T.K. Model based online diagnosis of unbalance and transverse fatigue crack in rotor systems. Mech. Res. Commun. 2004, 31, 557–568. [Google Scholar] [CrossRef]
  12. Eggers, A.J.; Digumarthi, R.; Chaney, K. Wind shear and turbulence effects on rotor fatigue and loads control. J. Sol. Energy Eng. 2003, 125, 402–409. [Google Scholar] [CrossRef]
  13. Huang, L.; Wu, Z.; Wang, K. Indirect measurement of rotor dynamic imbalance for control moment gyroscopes via gimbal disturbance observer. Sensors 2018, 18, 1873. [Google Scholar] [CrossRef] [Green Version]
  14. Qiu, J.; Tani, J.; Kwon, T. Control of self-excited vibration of a rotor system with active gas bearings. J. Vib. Acoust. 2003, 125, 328–334. [Google Scholar] [CrossRef]
  15. Jia, B.; Lei, H.; Yong, F.; Zhai, J.; Gang, Q. The development of aero-engine tip-clearance measurement technology: A simple review. In Proceedings of the 2017 IEEE 13th International Conference on Electronic Measurement & Instruments, Yangzhou, China, 20–23 October 2017. [Google Scholar]
  16. Yu, B.; Ke, H.; Shen, E.; Zhang, T. A review of blade tip clearance–measuring technologies for gas turbine engines. Meas. Control. 2020, 53, 339–357. [Google Scholar] [CrossRef] [Green Version]
  17. Garcia, I.; Zubia, J.; Beloki, J.; Arrue, J.; Durana, G.; Aldabaldetreku, G. Optical tip clearance measurements as a tool for rotating disk characterization. Sensors 2017, 17, 165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Garcia, I.; Przysowa, R.; Amorebieta, J.; Zubia, J. Tip-clearance measurement in the first stage of the compressor of an aircraft engine. Sensors 2016, 16, 1897. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Garcia, I.; Zubia, J.; Durana, G.; Aldabaldetreku, G.; Illarramendi, M.A.; Villatoro, J. Optical fiber sensors for aircraft structural health monitoring. Sensors 2015, 15, 15494–15519. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Guo, H.; Duan, F.; Wu, G.; Zhang, J. Blade tip clearance measurement of the turbine engines based on a multi-mode fiber coupled laser ranging system. Rev. Sci. Instrum. 2014, 85, 115105. [Google Scholar] [CrossRef] [PubMed]
  21. Garcia, I.; Beloki, J.; Zubia, J.; Aldabaldetreku, G.; Illarramendi, M.A.; Jimenez, F. An optical fiber bundle sensor for tip clearance and tip timing measurements in a turbine rig. Sensors 2013, 13, 7385–7398. [Google Scholar] [CrossRef] [PubMed]
  22. Jia, B.; Zhang, X. An optical fiber blade tip clearance sensor for active clearance control applications. Procedia Eng. 2011, 15, 984–988. [Google Scholar]
  23. Zhang, F.M.; Li, Y.T.; Pan, H.; Shi, C.Z.; Qu, X.H. Vibration compensation of the frequency-scanning-interferometry-based absolute ranging system. Appl. Sci. 2019, 9, 147. [Google Scholar] [CrossRef] [Green Version]
  24. Jia, X.Y.; Liu, Z.G.; Deng, Z.W.; Deng, W.; Wang, Z.Y.; Zhen, Z.L. Dynamic absolute distance measurement by frequency sweeping interferometry based Doppler beat frequency tracking model. Opt. Commun. 2019, 430, 163–169. [Google Scholar] [CrossRef]
  25. Lu, C.; Xiang, Y.; Gan, Y.; Liu, B.; Chen, F.; Liu, X.; Liu, G. FSI-based non-cooperative target absolute distance measurement method using PLL correction for the influence of a nonlinear clock. Opt. Lett. 2018, 43, 2098–2101. [Google Scholar] [CrossRef] [PubMed]
  26. Yu, W.; Pfeiffer, P.; Morsali, A.; Yang, J.; Fontaine, J. Comb-calibrated frequency sweeping interferometry for absolute distance and vibration measurement. Opt. Lett. 2019, 44, 5069–5072. [Google Scholar] [CrossRef] [PubMed]
  27. Jiang, S.; Liu, B.; Wang, H.; Zhao, B. Absolute distance measurement using frequency-scanning interferometry based on hilbert phase subdivision. Sensors 2019, 19, 5132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Tao, L.; Liu, Z.; Zhang, W.; Liu, Z.; Hong, J. Real-time drift error compensation in a self-reference frequency-scanning fiber interferometer. Opt. Commun. 2017, 382, 99–104. [Google Scholar] [CrossRef]
  29. Tao, L.; Liu, Z.; Zhang, W.; Zhou, Y. Frequency-scanning interferometry for dynamic absolute distance measurement using Kalman filter. Opt. Lett. 2014, 39, 6997–7000. [Google Scholar] [CrossRef] [PubMed]
  30. Shafir, E.; Berkovic, G. Compact fibre optic probe for simultaneous distance and velocity determination. Meas. Sci. Technol. 2001, 12, 943–947. [Google Scholar]
  31. Schneider, R.; Thurmel, P.; Stockmann, M. Distance measurement of moving objects by frequency modulated laser radar. Opt. Eng. 2001, 40, 33–37. [Google Scholar] [CrossRef]
  32. Shao, B.; Zhang, W.; Zhang, P.; Chen, W. Dynamic clearance measurement using fiber-optic frequency-swept and frequency-fixed interferometry. IEEE Photonics Technol. Lett. 2020, 32, 1331–1334. [Google Scholar] [CrossRef]
Figure 1. (a) The ideal case, and (b) the real situation of the running rotor. (L: axial clearance; D: radial clearance; α: tilt angle; Ω: angular velocity.).
Figure 1. (a) The ideal case, and (b) the real situation of the running rotor. (L: axial clearance; D: radial clearance; α: tilt angle; Ω: angular velocity.).
Sensors 20 07178 g001
Figure 2. Schematic of the fiber-optic FSI system. (FSL: frequency-swept laser; FOC: fiber-optic circulator; SMF: single-mode fiber; PD: photodetector; DAQ: data acquisition system.).
Figure 2. Schematic of the fiber-optic FSI system. (FSL: frequency-swept laser; FOC: fiber-optic circulator; SMF: single-mode fiber; PD: photodetector; DAQ: data acquisition system.).
Sensors 20 07178 g002
Figure 3. Schematic of the multi-parameter measurement of the rotor.
Figure 3. Schematic of the multi-parameter measurement of the rotor.
Sensors 20 07178 g003
Figure 4. (a) Schematic of the experimental system, and (b) photograph of the rotational motion setup. (FSL: frequency-swept laser; FOC: fiber-optic circulator; SMF: single-mode fiber; PD: photodetector; DAQ: data acquisition system; MSC: motorized stage controller; SG: signal generator; VA: voltage amplifier.).
Figure 4. (a) Schematic of the experimental system, and (b) photograph of the rotational motion setup. (FSL: frequency-swept laser; FOC: fiber-optic circulator; SMF: single-mode fiber; PD: photodetector; DAQ: data acquisition system; MSC: motorized stage controller; SG: signal generator; VA: voltage amplifier.).
Sensors 20 07178 g004
Figure 5. Variation of the actual distance L(θ) under the condition of Ω = 0 rad/s.
Figure 5. Variation of the actual distance L(θ) under the condition of Ω = 0 rad/s.
Sensors 20 07178 g005
Figure 6. Measured LRo_meas(t) and theoretical LRo_theo(t) under the condition of Ω = 60π rad/s.
Figure 6. Measured LRo_meas(t) and theoretical LRo_theo(t) under the condition of Ω = 60π rad/s.
Sensors 20 07178 g006
Figure 7. Measurement results of (a) angular velocity, (b) axial clearance, and tilt angle under different angular velocities.
Figure 7. Measurement results of (a) angular velocity, (b) axial clearance, and tilt angle under different angular velocities.
Sensors 20 07178 g007
Figure 8. Measurement results of (a) axial clearance, (b) angular velocity, and tilt angle under different axial clearances.
Figure 8. Measurement results of (a) axial clearance, (b) angular velocity, and tilt angle under different axial clearances.
Sensors 20 07178 g008
Figure 9. Measurement results of (a) tilt angle, (b) axial clearance, and angular velocity under different radii.
Figure 9. Measurement results of (a) tilt angle, (b) axial clearance, and angular velocity under different radii.
Sensors 20 07178 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shao, B.; Zhang, W.; Zhang, P.; Chen, W. Multi-Parameter Measurement of Rotors Using the Doppler Effect of Frequency-Swept Interferometry. Sensors 2020, 20, 7178. https://doi.org/10.3390/s20247178

AMA Style

Shao B, Zhang W, Zhang P, Chen W. Multi-Parameter Measurement of Rotors Using the Doppler Effect of Frequency-Swept Interferometry. Sensors. 2020; 20(24):7178. https://doi.org/10.3390/s20247178

Chicago/Turabian Style

Shao, Bin, Wei Zhang, Peng Zhang, and Weimin Chen. 2020. "Multi-Parameter Measurement of Rotors Using the Doppler Effect of Frequency-Swept Interferometry" Sensors 20, no. 24: 7178. https://doi.org/10.3390/s20247178

APA Style

Shao, B., Zhang, W., Zhang, P., & Chen, W. (2020). Multi-Parameter Measurement of Rotors Using the Doppler Effect of Frequency-Swept Interferometry. Sensors, 20(24), 7178. https://doi.org/10.3390/s20247178

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop