# Phase Demodulation of Rotor Torsional Vibration Measurement under Time-Varying Speed

Authors: Moli Chen; Nan Zheng; Jun Zheng; Zhu Jin; Guihuo Luo
Journal: Shock and Vibration
 (2022)

---

## Abstract

When torsional vibration is measured with the zebra-tape method, the modulation signals can be contaminated by electrical noise, environmental noise, step noise, sensor nonlinearities, and so on. Although the amplitudes of these types of noise are very weak relative to the amplitude of the ideal pulse test signal, they may introduce large-scale random broad-spectrum noise to the signal phases. Such noise is difficult to remove using frequency domain, time-frequency-domain, or threshold sampling methods and seriously affects the accuracy of torsional vibration measurements. This paper presents a phase demodulation algorithm based on downsampling and local resampling (DSLR) to improve the accuracy of torsional vibration measurements. To verify the proposed DSLR algorithm, torsional vibration simulations and experiments are conducted under a time-varying rotation speed. The results show that the angular displacement signal obtained from the torsional vibration test using zebra tapes is the superposition of the dynamic and static angular displacements, with the latter generated because of the varying rotation speed with respect to time. The DSLR algorithm effectively reduces the phase noise and demodulates the pulse signal phases. The work in this paper provides a method for obtaining reference torsional vibration measurements of rotors based on the zebra-tape method under time-varying speeds.


---

## Body

## 1. Introduction

The speed of rotation of propellers, aero-engines, or gas turbines fluctuates greatly as they start up, accelerate, decelerate, and stop. The rapidly time-varying rotation speed causes torsional vibration in the rotor shaft. If the instantaneous speed of the rotor reaches the critical torsional speed and remains there for a certain period of time, the excessive torsional vibrations that are produced will seriously threaten the safety of the body. Torsional vibration tests are important means of revealing the variation law of rotor torsional vibration under time-varying speeds.The zebra-tape method is commonly used to measure rotor torsional vibration [1, 2], as the tapes can be installed on any part of the shaft that is physically accessible [3]. Using noncontact magnetoelectric or eddy sensor with a high-frequency response, the pulse signal is measured as the tapes pass through the sensor probe. By demodulating the pulse phases and extracting the torsional vibration signal, the torsional vibration can be measured. In this process, pulse phase demodulation is one of the most important tasks, and the accuracy of the phases extracted from the pulse signal directly affects the accuracy of the torsional vibration measurements.The bending vibration and torsional vibration of machines, nonlinearities in the sensor, and defects in the tape mean that the pulse signal observed by the sensor will not be an ideal square wave or sine wave but will instead feature some modulated frequency, amplitude, and phase components. In addition, the modulated signal will have some interference from the geometric errors inherent in toothed disk structures [4, 5], the electric noise of cable transmission, environmental noise, and so on. All of these factors directly affect the pulse arrival time of the modulated signal, resulting in pulse phase noise [6]. Thus, phase demodulation is a very difficult task in the torsional vibration measurement process.To date, various phase demodulation methods have been developed, such as the fast Fourier transform [7, 8], Hilbert–Huang transform [9], filtering [10, 11], order tracking filter [12, 13], wavelet decomposition [14], variable mode decomposition [15, 16], and harmonic signal decomposition [17]. These algorithms are mostly based on the theory of frequency domain or time-frequency domain transformation and provide effective torsional vibration phase demodulation with a clear frequency distribution, such as for harmonic components and the geometric errors of toothed disk structures.For electrical noise and environmental noise, however, although the amplitudes of such noises are very weak relative to the amplitude of the pulse signal, they will cause extreme pulse redundancy. This redundancy makes it very difficult to calculate the pulse duty cycle. Moreover, the phase noise caused by this redundancy is random broad-spectrum noise, so it is difficult to denoise using frequency domain or time-frequency domain algorithms, which seriously affects the accuracy of pulse phase demodulation.To improve the accuracy of torsional phase demodulation based on zebra tapes under a time-varying rotation speed, this paper derives the relationship between the pulse phase difference of the zebra tapes, the static angular displacement caused by the time-varying rotational speed, and the torsional angular displacement. A phase demodulation algorithm based on downsampling and local resampling (DSLR) is then proposed. The accuracy of this algorithm is verified through simulation examples of torsional vibration. The algorithm is applied to measure the torsional vibration angle displacement of a propeller shaft under the starting transition process, and the phase demodulation algorithm is further evaluated.


## 2. Testing Theory

As shown in Figure1, zebra tapes are installed on the rotating shaft sections, and two noncontact magnetoelectric sensors are used to collect the pulse signals of the tapes passing through the sensors. One of the two sensors is used to observe/measure the phases of torsional vibration, while the other records the reference phases. The signals from the sensors are amplified and adjusted by an acquisition instrument and then transmitted to a computer.Figure 1
Schematic diagram of torsion measuring devices on zebra-tape method.When torsional vibration occurs in the rotor, the time at which the tapes pass the sensors advances or recedes, and this modulates the phases of the test signals. The phase differences (i.e., the advanced or lagged time) caused by phase modulation can be written as(1)τt=tτt−t0t.tτt is the time at which the pulse extrema (i.e., maxima or minima) are generated by the tapes passing sensor1 under the condition of torsional vibration, and t0t is the time at which the pulse extrema without torsional vibration pass sensor2. We call tτt the observation phase and t0t the reference phase; the corresponding signals are called the observation signal and the reference signal.According to (1), the angular displacement response (in degrees) for the rotor can be expressed as(2)θt=360×nt60τt=θτt−θ0t,θτt=6nttτt,θ0t=6ntt0t,where n(t) is the rotation speed of the rotor (rpm). Even if the expressions of θt, θτt, andθ0t are discretized, (2) shows that θt can be calculated when the length of θτt is different from that of θ0t because of the different zebra tapes, with digital interpolation used to ensure that the data sequences have the same length.If we consider the variable static angular displacement generated by a time-varying speed, (2) needs to be modified. The modification method is described as follows.The relationship between the rotor control speedn¯t and the static torque T¯t applied to the rotor is as follows:(3)n¯t=9549PtT¯t=9549PtKn¯θ¯t,where P(t) is the input power for the rotor and Kn¯ is the slowly varying torsional stiffness of the observed shaft segment under the control speed n¯t as the gyroscope torque changes with the bending-torsion coupled vibration. (The number 9549 is the conversion coefficient when the speed is in rpm, the power is in kW, and Kn¯ is in N·m/°.)According to (3), the static torsional angle of the rotor is given by(4)θ¯t=9549PtKn¯n¯t.It can be inferred from (2) and (4) that the angular displacement of the rotor obtained from the test signals is the superposition of the dynamic and static angular displacements. Thus, (2) can be modified to(5)θt=θ∼t+θ¯t=6ntτt,where θ%t is the vibrational angular displacement, which is a high-frequency component compared with the static displacement component θ¯t induced by the time-varying speed.According to (5), the vibrational angular displacement of the rotor under the time-varying speed is given by(6)θ%t=6ntτt−θ¯t.In practice, it is not necessary to calculateθ¯t with (4) after obtaining the rotor speed nt and the phase difference τt, because θ¯t can be determined and eliminated using the trend term decomposition method. Once θt has been detrended, θ%t can be identified.


## 3. Phase Demodulation Algorithm

From (1) and (2), determining the angular response of the torsional vibration requires t0t and tτt, which are the times of the extrema in the pulse signals. Therefore, an algorithm based on DSLR is constructed to extract the maxima from the pulse signal. First, the upper envelope of the pulse signal is denoised by the downsampling method, and then, the maximum sequence is obtained by locally resampling the original pulse signal. After calculating the angular displacement of the torsional vibration using (6), the phase demodulation of the signals from the two sensors can be realized.


### 3.1. Denoising the Pulse Phases

We extract the extrema of a pulse signal using a downsampling method to decrease the redundant noise of the extrema.The interference from weak electrical noise, environmental noise, and step noise makes it difficult to identify the pulse extrema directly using the envelope algorithm or threshold algorithm. The purpose of downsampling is to reduce the occurrence of weak random noises in the pulse signal and avoid redundant extrema, so as to improve the extraction accuracy of pulse phases using more accurate extrema.Taking phase denoising based on the extraction of maxima as an example, the downsampling pulse phase denoising algorithm is outlined as follows.Let the signal of the sensor be(7)X=Xk=xtk|k=1,2,…,N,where N is the sampling length of signal X.If the sampling frequency ofX is fs, the maximal testing speed of the rotor is nmax, and the number of tapes is Z, and then, the interval Δk of the downsampling should satisfy the following formula:(8)2<Δk<15fsZnmax.The downsampling method resamples data from the original signal everyΔk points. The physical meaning of 15fs/Znmax is a quarter of the number of sampling points in one rotation at the maximal test speed (⋅ represents the rounding down operator). Even if the method decreases the signal sampling frequency, the final sampling frequency is still greater than four times the corresponding frequency of the maximum test speed. Thus, it meets the requirements of the sampling theorem for the digital acquisition of speed pulse signals, and no frequency aliasing will occur.Downsampling to the sequenceX yields(9)X∼=X∼i=xt1+iΔk|i=0,1,2,…,NΔk.The downsampling process is illustrated in Figure2. The maximum B is less than the minimum A in Figure 2(a)), so the threshold sampling method cannot be applied to determine the maxima of the signal. Though the accuracy is degraded, the downsampling method can remove weak noise from the signal by decreasing the sampling frequency (Figure 2(b)). The effect is better at larger signal-to-noise ratios or in the case of less noise.Figure 2
Downsampling sequence: (a) the original signal (X) and (b)X∼ in (X).We first look for the maxima ofX% and then search for the maxima of the original signal X near the peaks of X%. The procedure for identifying the peaks of X% is as follows.The differential sequencedX% is obtained as follows:(10)dX%=dXi%=Xi+1%−Xi%|i=0,1,2,…,NΔk−1.If theith point in sequence dX% satisfies(11)dXi%>0,dXi%dXi+1%≤0,thenXi+1% is the maximum of the downsampling sequence X%. The points that satisfy (11) can be expressed as the sequence(12)X⌢0=Xi+1%|dXi%>0anddXi%dXi+1%≤0,where X⌢0 is the result of sequence X being downsampled, the physical meaning of which is the upper envelope of the sequence X after being sampled 1/Δk times. This gives the maxima in the time period Δk/f. The result of the downsampling method is described in Figure 3, where the triangles denote the components of sequence X⌢0.Figure 3
Downsampling result.


### 3.2. Extracting the Pulse Phases

X⌢0 is not the maximum sequence of X. To improve the accuracy of the phase calculation and obtain the upper envelope of X, we need to resample X locally.Let the resampling interval be 2Δk. If X˜i+1=xt1+iΔk+Δk is an element of X⌢0, the locally resampling sequence is defined by(13)Xi+1L=xt1+iΔk,xt1+iΔk+1,xt1++iΔk+2,L,xt1++iΔk+Δk,xt1++iΔk+Δk+1,xt1++iΔk+Δk+2,xt1+iΔk+2Δk−1.The meaning ofXi+1L is illustrated in Figure 4.Figure 4
Local resampling sequence.The maximum of sequenceXi+1L (Figure 5) can be expressed as(14)X⌢i+1=maxXi+1L.Figure 5
Local resampling result: (a) local resampling sequence and (b) the extracted pulse phase point.Then, the upper envelope sequence ofX is given by(15)X⌢=X⌢i+1.Through this resampling procedure, the maximum sequenceX⌢ is extracted from X by determining the maxima within 2 Δk of the points in X⌢0. This algorithm avoids the interference of weak noise and improves the accuracy of the calculated upper envelope. The upper envelope of the signal X is clarified through (15).The pulse phases extracted fromX are the abscissa sequence t⌢ to X⌢ in (15), which can be expressed as(16)t⌢=t⌢1,t⌢2,....The physical meaning of DSLR is further illustrated in Figure6, where the crosses denote extrema that are not peaks in their cycles. These are the redundant points to be denoised.Figure 6
The physical meaning of the DSLR phase demodulation algorithm.


### 3.3. Phase Demodulation for Torsional Vibration

Using the DSLR algorithm, the pulse phase sequencest⌢π and t⌢0 can be obtained from the observation signal and the reference signal. The phase difference sequence is then calculated as(17)τ=τi=t⌢τi−t⌢oi|i=1,2,....By substituting (17) into (2), the angular displacement response of the rotor is given by(18)θ=θi=6niτi|i=1,2,...,where ni is the discrete speed, which can be measured experimentally.Note thatθ is a function of time t, but its elements are nonuniform samples after the envelope extraction operation. To obtain the torsional vibration signals at equal time intervals, an interpolation operator must be applied to the sequence θ; this agrees with Badaoui [18].According to (6), if only the dynamic vibration signal θ% is needed, the static component θ¯t can be determined using a trending algorithm (e.g., moving average method). This component can then be removed from θ in (18) to obtain a pure vibration signal.


### 3.4. Phase Demodulation Method

A flowchart of the phase demodulation method for torsional vibration is shown in Figure7. The proof about the DSLR algorithm is in appendix A.Figure 7
Flowchart of phase demodulation algorithm.


## 3.1. Denoising the Pulse Phases

We extract the extrema of a pulse signal using a downsampling method to decrease the redundant noise of the extrema.The interference from weak electrical noise, environmental noise, and step noise makes it difficult to identify the pulse extrema directly using the envelope algorithm or threshold algorithm. The purpose of downsampling is to reduce the occurrence of weak random noises in the pulse signal and avoid redundant extrema, so as to improve the extraction accuracy of pulse phases using more accurate extrema.Taking phase denoising based on the extraction of maxima as an example, the downsampling pulse phase denoising algorithm is outlined as follows.Let the signal of the sensor be(7)X=Xk=xtk|k=1,2,…,N,where N is the sampling length of signal X.If the sampling frequency ofX is fs, the maximal testing speed of the rotor is nmax, and the number of tapes is Z, and then, the interval Δk of the downsampling should satisfy the following formula:(8)2<Δk<15fsZnmax.The downsampling method resamples data from the original signal everyΔk points. The physical meaning of 15fs/Znmax is a quarter of the number of sampling points in one rotation at the maximal test speed (⋅ represents the rounding down operator). Even if the method decreases the signal sampling frequency, the final sampling frequency is still greater than four times the corresponding frequency of the maximum test speed. Thus, it meets the requirements of the sampling theorem for the digital acquisition of speed pulse signals, and no frequency aliasing will occur.Downsampling to the sequenceX yields(9)X∼=X∼i=xt1+iΔk|i=0,1,2,…,NΔk.The downsampling process is illustrated in Figure2. The maximum B is less than the minimum A in Figure 2(a)), so the threshold sampling method cannot be applied to determine the maxima of the signal. Though the accuracy is degraded, the downsampling method can remove weak noise from the signal by decreasing the sampling frequency (Figure 2(b)). The effect is better at larger signal-to-noise ratios or in the case of less noise.Figure 2
Downsampling sequence: (a) the original signal (X) and (b)X∼ in (X).We first look for the maxima ofX% and then search for the maxima of the original signal X near the peaks of X%. The procedure for identifying the peaks of X% is as follows.The differential sequencedX% is obtained as follows:(10)dX%=dXi%=Xi+1%−Xi%|i=0,1,2,…,NΔk−1.If theith point in sequence dX% satisfies(11)dXi%>0,dXi%dXi+1%≤0,thenXi+1% is the maximum of the downsampling sequence X%. The points that satisfy (11) can be expressed as the sequence(12)X⌢0=Xi+1%|dXi%>0anddXi%dXi+1%≤0,where X⌢0 is the result of sequence X being downsampled, the physical meaning of which is the upper envelope of the sequence X after being sampled 1/Δk times. This gives the maxima in the time period Δk/f. The result of the downsampling method is described in Figure 3, where the triangles denote the components of sequence X⌢0.Figure 3
Downsampling result.


## 3.2. Extracting the Pulse Phases

X⌢0 is not the maximum sequence of X. To improve the accuracy of the phase calculation and obtain the upper envelope of X, we need to resample X locally.Let the resampling interval be 2Δk. If X˜i+1=xt1+iΔk+Δk is an element of X⌢0, the locally resampling sequence is defined by(13)Xi+1L=xt1+iΔk,xt1+iΔk+1,xt1++iΔk+2,L,xt1++iΔk+Δk,xt1++iΔk+Δk+1,xt1++iΔk+Δk+2,xt1+iΔk+2Δk−1.The meaning ofXi+1L is illustrated in Figure 4.Figure 4
Local resampling sequence.The maximum of sequenceXi+1L (Figure 5) can be expressed as(14)X⌢i+1=maxXi+1L.Figure 5
Local resampling result: (a) local resampling sequence and (b) the extracted pulse phase point.Then, the upper envelope sequence ofX is given by(15)X⌢=X⌢i+1.Through this resampling procedure, the maximum sequenceX⌢ is extracted from X by determining the maxima within 2 Δk of the points in X⌢0. This algorithm avoids the interference of weak noise and improves the accuracy of the calculated upper envelope. The upper envelope of the signal X is clarified through (15).The pulse phases extracted fromX are the abscissa sequence t⌢ to X⌢ in (15), which can be expressed as(16)t⌢=t⌢1,t⌢2,....The physical meaning of DSLR is further illustrated in Figure6, where the crosses denote extrema that are not peaks in their cycles. These are the redundant points to be denoised.Figure 6
The physical meaning of the DSLR phase demodulation algorithm.


## 3.3. Phase Demodulation for Torsional Vibration

Using the DSLR algorithm, the pulse phase sequencest⌢π and t⌢0 can be obtained from the observation signal and the reference signal. The phase difference sequence is then calculated as(17)τ=τi=t⌢τi−t⌢oi|i=1,2,....By substituting (17) into (2), the angular displacement response of the rotor is given by(18)θ=θi=6niτi|i=1,2,...,where ni is the discrete speed, which can be measured experimentally.Note thatθ is a function of time t, but its elements are nonuniform samples after the envelope extraction operation. To obtain the torsional vibration signals at equal time intervals, an interpolation operator must be applied to the sequence θ; this agrees with Badaoui [18].According to (6), if only the dynamic vibration signal θ% is needed, the static component θ¯t can be determined using a trending algorithm (e.g., moving average method). This component can then be removed from θ in (18) to obtain a pure vibration signal.


## 3.4. Phase Demodulation Method

A flowchart of the phase demodulation method for torsional vibration is shown in Figure7. The proof about the DSLR algorithm is in appendix A.Figure 7
Flowchart of phase demodulation algorithm.


## 4. Simulation of the Phase Demodulation of Torsional Vibration

As the focus of this paper is the extraction of the torsional vibration signal, the static torsional angle component caused by the time-varying speed is not considered in the simulation analysis.The simulation model with a time-varying rotation speed is expressed as(19)ωt=2π60n1sin0.2πt+2+10πcos0.6πt+0.8+2π60n2.Here, the units ofω are rad/s, n1=1000 rpm, and n2=2000 rpm.Assume that the number of zebra tapesZ = 60. The sampling frequency of the pulse signal is set to fs=200kHz and the number of sampling points is N=3.5×105 to reduce the distortion of the pulse signal. A frequency of 200 kHz ensures that the number of sampling points on every tape passing the rotation signal is not less than 62 (2 × 105/3200 = 62.5, where 3200 rpm is the maximum speed). The speed with respect to time is shown in Figure 8.Figure 8
The simulation model of time-varying speed.


### 4.1. Evaluation of Reference Phases

Based on the known time-varying speed, the reference phaset0t can be obtained by the phase demodulation of the sensor signal under the condition of no torsional vibration. At the same time, in the process of estimating the reference phase, the phase denoising and extraction performance of the DSLR method will be verified.We simulate the reference signal with no torsional vibration as follows:(20)xt=sinZωt+0.02η1t,where η1t is a uniformly distributed random noise in the interval (0, 1) and t is the uniform speed-up time. The evolution of xt in Figure 9 is a harmonic signal with noise, but with no advanced or lagged phase term. Thus, we use xt to simulate the reference signal with zero torsional vibration and no phase modulation.Figure 9
The pulse signal of nontorsional vibration: (a)xt and (b) the local amplification of xt.
(a)(b)To verify the superiority of the proposed denoising and extraction algorithms, we apply both the DSLR algorithm and the conventional peak-finding algorithm to extract the pulse phases ofxt. According to (8), we set the downsampling interval Δk  = 6 in the DSLR algorithm. The results are shown in Figure 10. The pulse phases extracted by the DSLR method are consistent with the peak values of the original reference signal with no redundancy or omission, whereas those given by the conventional peak-finding method include the redundant phase points caused by the noise signal. In Figure 10, the stars that do not coincide with the circles are redundant. These redundant points are caused by sensor or environmental noise, and their number and phase are uncertain, which makes it difficult to match the demodulated phase sequence t⌢τ to the reference phase sequence t⌢0. This complicates the demodulation of the torsional vibration using (1) or (17). The accuracy of the zebra-tape method is seriously affected by these redundant points.Figure 10
Comparison of DSLR and find-peaks methods ont0t.The reference phase sequencet0t extracted by the DSLR method is shown in Figure 11. The DSLR method is very effective in denoising and extracting the pulse phases.Figure 11
Reference phaset0t.


### 4.2. Extraction of Torsional Vibration

In this section, we build pulse samples from a torsional sensor and simulate and verify the torsional vibration demodulation algorithm.Let us assume that the torsional angular displacement signal is as follows:(21)θt=2sin64πt.The advance time of the pulse signal caused by torsional vibration is given by(22)τt=πθt180ωt.Modulating the phase ofxt in (20) with τt, the pulse signal of the torsion sensor with noise is given by(23)yt^=sinZωt+0.02η2t,where t^=t+τt and η2t is uniformly distributed random noise in the interval (0, 1).Comparing (23) with (20), it can be seen that the variable t^ in (23) is the result of the torsional angular displacement modulation of the variable t in (20). Besides the different sensor noise, the right-hand terms of both equations are consistent, which ensures that yt^ is the phase modulation result of xt.Figure12 shows the phase difference between the signal pulses with and without torsional vibration. The torsional vibration described by (21) is a constant amplitude, constant frequency harmonic signal. Under the modulation of time-varying speed, the whole frequency trend of the phase difference τt decreases with increasing speed, while the amplitude increases.Figure 12
Phase differenceτt.Figure13 shows the pulse phases extracted by the DSLR method and the conventional peak-finding method in the presence of torsional vibration. The curve is similar to the calculation result in Figure 10, and the superiority of the proposed algorithm is again apparent.Figure 13
Comparison of DSLR and find-peaks methods ontτt.The modulation phase signaltτt obtained by the phase extraction of the modulation signal is shown in Figure 14. This has the same length as the reference phase t0t, so we can demodulate the phases of the torsional vibration using (17). The true and calculated torsional signals of the rotor are shown in Figure 15. Compared with the true torsional signal in Figure 15(a)), the demodulated torsional signal in Figure 15(b) is not distorted in frequency, amplitude, or phase—the largest error in the amplitude is only 1.3%, which verifies the correctness and high accuracy of the proposed DSLR algorithm for the demodulation of torsional vibration.Figure 14
Modulated phasestτt.Figure 15
Results of torsional vibration angle: (a) actual and (b) calculated.
(a)(b)


## 4.1. Evaluation of Reference Phases

Based on the known time-varying speed, the reference phaset0t can be obtained by the phase demodulation of the sensor signal under the condition of no torsional vibration. At the same time, in the process of estimating the reference phase, the phase denoising and extraction performance of the DSLR method will be verified.We simulate the reference signal with no torsional vibration as follows:(20)xt=sinZωt+0.02η1t,where η1t is a uniformly distributed random noise in the interval (0, 1) and t is the uniform speed-up time. The evolution of xt in Figure 9 is a harmonic signal with noise, but with no advanced or lagged phase term. Thus, we use xt to simulate the reference signal with zero torsional vibration and no phase modulation.Figure 9
The pulse signal of nontorsional vibration: (a)xt and (b) the local amplification of xt.
(a)(b)To verify the superiority of the proposed denoising and extraction algorithms, we apply both the DSLR algorithm and the conventional peak-finding algorithm to extract the pulse phases ofxt. According to (8), we set the downsampling interval Δk  = 6 in the DSLR algorithm. The results are shown in Figure 10. The pulse phases extracted by the DSLR method are consistent with the peak values of the original reference signal with no redundancy or omission, whereas those given by the conventional peak-finding method include the redundant phase points caused by the noise signal. In Figure 10, the stars that do not coincide with the circles are redundant. These redundant points are caused by sensor or environmental noise, and their number and phase are uncertain, which makes it difficult to match the demodulated phase sequence t⌢τ to the reference phase sequence t⌢0. This complicates the demodulation of the torsional vibration using (1) or (17). The accuracy of the zebra-tape method is seriously affected by these redundant points.Figure 10
Comparison of DSLR and find-peaks methods ont0t.The reference phase sequencet0t extracted by the DSLR method is shown in Figure 11. The DSLR method is very effective in denoising and extracting the pulse phases.Figure 11
Reference phaset0t.


## 4.2. Extraction of Torsional Vibration

In this section, we build pulse samples from a torsional sensor and simulate and verify the torsional vibration demodulation algorithm.Let us assume that the torsional angular displacement signal is as follows:(21)θt=2sin64πt.The advance time of the pulse signal caused by torsional vibration is given by(22)τt=πθt180ωt.Modulating the phase ofxt in (20) with τt, the pulse signal of the torsion sensor with noise is given by(23)yt^=sinZωt+0.02η2t,where t^=t+τt and η2t is uniformly distributed random noise in the interval (0, 1).Comparing (23) with (20), it can be seen that the variable t^ in (23) is the result of the torsional angular displacement modulation of the variable t in (20). Besides the different sensor noise, the right-hand terms of both equations are consistent, which ensures that yt^ is the phase modulation result of xt.Figure12 shows the phase difference between the signal pulses with and without torsional vibration. The torsional vibration described by (21) is a constant amplitude, constant frequency harmonic signal. Under the modulation of time-varying speed, the whole frequency trend of the phase difference τt decreases with increasing speed, while the amplitude increases.Figure 12
Phase differenceτt.Figure13 shows the pulse phases extracted by the DSLR method and the conventional peak-finding method in the presence of torsional vibration. The curve is similar to the calculation result in Figure 10, and the superiority of the proposed algorithm is again apparent.Figure 13
Comparison of DSLR and find-peaks methods ontτt.The modulation phase signaltτt obtained by the phase extraction of the modulation signal is shown in Figure 14. This has the same length as the reference phase t0t, so we can demodulate the phases of the torsional vibration using (17). The true and calculated torsional signals of the rotor are shown in Figure 15. Compared with the true torsional signal in Figure 15(a)), the demodulated torsional signal in Figure 15(b) is not distorted in frequency, amplitude, or phase—the largest error in the amplitude is only 1.3%, which verifies the correctness and high accuracy of the proposed DSLR algorithm for the demodulation of torsional vibration.Figure 14
Modulated phasestτt.Figure 15
Results of torsional vibration angle: (a) actual and (b) calculated.
(a)(b)


## 5. Torsional Vibration Test of a Propeller Shaft

The torsional vibration of a civil aircraft propeller shaft was measured under the start-up condition. Figure16 shows the test devices used to induce torsional vibration in the propeller shaft. The speed controller adjusted the speed of the engine according to the output power. A reducer with a reduction ratio of 1.69 allowed the engine to rotate the propeller. Zebra tapes were installed on the input shaft of the propeller, and two magnetoelectric sensors were used to collect the torsional vibration signals. One sensor recorded the passing pulses from the zebra tapes on the propeller shaft, the signal of which was used to determine the observation phases of torsional vibration. The other sensor recorded the passing pulses from the toothed disk on the engine rotor, the signal of which was used to determine the reference phases. The pulse signals from both sensors were output to an acquisition conditioner before being sent to a computer.Figure 16
The torsional vibration test apparatus.The zebra tapes were machined specially from 0.2-mm-thick metal material for the test. First, 60 uniform zebra tapes were glued on the input shaft of the propeller. The tapes at 0° and 180° in the circumferential direction were then torn off to provide a rotation speed flag and for rotor balance. The disk controlling the engine speed had 60 teeth, but the zebra tape at 0° was removed for the speed test. A torsional spring was installed at the connection between the toothed disk and the shaft of the engine rotor, so as to weaken the influence of shaft torsion and improve the accuracy of the engine speed measurements. The toothless disk and the torsional spring are original components of the engine rotor itself rather than added for this torsional vibration test.The time-varying speed of the propeller shaft at a signal sampling frequency offs = 216 kHz is shown in Figure 17. The speed of the propeller shaft climbed slowly from 1200 rpm with some slight fluctuations up to 28 s and then entered the fast nonuniform climbing stage. The signals of the test sensors are shown in Figures 18 and 19. Figure 18 shows that the output voltage of the propeller sensor increased continuously, indicating some bending vibration on the input shaft of the propeller. The missing tapes modulated the amplitude, frequency, and phase of the observation signal of the propeller shaft, as shown in Figure 19.Figure 17
Time-varying speed of a propeller shaft.Figure 18
Observation signal of the propeller shaft.Figure 19
Amplified observation signal of the propeller shaft.The reference signals from the engine rotor are shown in Figure20 and 21. Comparing Figure 19 with Figure 21, it can be seen that the number of missing zebra tapes on the propeller shaft is twice that of the toothed disk, which is consistent with the actual working condition. The torsional spring ensures that highly accurate reference phases can be obtained by demodulating the signal in Figure 21.Figure 20
Reference signal from the engine rotor.Figure 21
Amplified reference signal.Using the DSLR algorithm, we demodulated the pulse signals in Figures18 and 20 to obtain tτt and t0t in turn. From formula (2) and the speed signal in Figure 17, we obtain θτt and θ0t, respectively; the speed signal to calculate θ0t, especially, was 1.69 times that in Figure 17.However, the lengths of discreteθτt and discrete θ0t were unequable, because of the different speeds of the propeller and the engine rotor leading both pulse signals had different pulse numbers, we cannot use formula (2) to calculate θt directly.Given the reduction ratio is 1.69, the length of discreteθ0t is 1.69 times that of θτt; thus, an operator about linear interpolation to θ0t was done to make the lengths of discrete θτt and discrete θ0t equal, and then, formula (2) was used to calculate torsional vibration θt. The final demodulated torsional vibration signal is shown in Figure 22.Figure 22
Torsional angle with static torsion effect of the propeller shaft.The propeller shaft exhibited slowly varying static torsional displacement as the speed changes with time, and it was speculated that the reason introducing the static term is shown in formula (3). The torsional displacement displayed a downward trend as the speed increases. The possible reason was that the running power of the propeller gradually matched the power output to it by the engine, which reduced the torsional vibration.Applying the detrending operator to the demodulated signal in Figure22 not only removes the influence of static torsional displacement caused by time-varying speed but also weakens the interference noise caused by the missing zebra tapes. The torsional vibration angles of the propeller shaft with respect to time and speed are shown in Figures 23 and 24. The torsional vibration of the propeller shaft decreases with speed, and reaches a maximum value at 770 rpm. Contrasting this with Figure 17, it appears that the torsional vibration response decreases as the speed changes rapidly from 800 rpm to 900 rpm, which indicates that the torsional vibration of the rotor is suppressed to a certain extent when the instantaneous speed becomes too fast. When the propeller shaft speed exceeds 900 rpm, although the propeller speed still increases rapidly, the torsional vibration response also increases. This indicates that the relationship between the time-varying speed and the torsional vibration response is nonlinear, and the time-varying speed is not the main factor affecting the torsional vibration of the rotor. Specifically, the shapes of the curves in Figure 24 before 800 rpm are caused by the slight jitter of the rotation speed.Figure 23
vibration angle vs. time of the propeller shaft.Figure 24
Vibration angle vs. rotation speed of the propeller shaft.


## 6. Conclusion

To extract the torsional angular displacement from pulse signals, this paper has proposed the DSLR phase demodulation algorithm and described its application to a working propeller shaft. Experiments were conducted in which the static angular displacement was considered in a torsional vibration test for a rotor with a time-varying speed based on the zebra-tape method.Under the condition of a time-varying rotational speed, the angular displacement obtained by the zebra-tape method was found to be the superposition of dynamic and static torsional angular displacements. After removing the static angular displacement component, the remaining torsional vibration signal is proportional to the rotation speed with some phase difference. The calculation of the phase difference is dependent on the phase of the observation signal and the reference phase.For a rotor torsional vibration pulse signal with time-varying speed and weak noise, the secondary sampling of the pulse signal based on DSLR can improve the demodulation accuracy of the torsional vibration phases.When the instantaneous rotation speed is too high, the torsional vibration response amplitude of the rotor is inhibited, but this is not always the most important factor affecting the torsional vibration response amplitude of the rotor.



---

*Source: 9222022-2022-04-11.xml*