Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
A Multi-Layer Feature Fusion Method for Few-Shot Image Classification
Previous Article in Journal
High-Throughput Effect-Directed Monitoring Platform for Specific Toxicity Quantification of Unknown Waters: Lead-Caused Cell Damage as a Model Using a DNA Hybrid Chain-Reaction-Induced AuNPs@aptamer Self-Assembly Assay
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A NARX Model-Based Condition Monitoring Method for Rotor Systems

1
School of Mechanical Engineering and Automation, Northeastern University, Shenyang 110819, China
2
School of Engineering and Material Science, Queen Mary University of London, London E1 4NS, UK
*
Authors to whom correspondence should be addressed.
Sensors 2023, 23(15), 6878; https://doi.org/10.3390/s23156878
Submission received: 8 July 2023 / Revised: 29 July 2023 / Accepted: 31 July 2023 / Published: 3 August 2023
(This article belongs to the Section Physical Sensors)

Abstract

:
In this study, we developed a data-driven frequency domain analysis method for rotor systems using the NARX (Nonlinear Auto-Regressive with eXternal input) model established by system vibration signals. We propose a model-based index of fault features calculated in a multi-frequency range to facilitate condition monitoring of rotor systems. Four steps are included in the proposed method. Firstly, displacement vibration signals are collected at multiple monitored rotating speeds. Secondly, the collected signals are processed as output data and the corresponding input data is generated. Then, NARX models are developed with input and output data to characterize the rotor system. Finally, the NRSF (Nonlinear Response Spectrum Function)-based nonlinear fault index is calculated and compared to the healthy condition. An experimental application to the misaligned rotor system is also demonstrated to verify its effectiveness. Our results indicate that the value of the index directly reflects the severity of the misaligned fault.

1. Introduction

Rotating machinery, such as engines, fans, water pumps, etc., plays an important role in aerospace, energy, transportation, petrochemical, and other industrial fields [1,2,3]. With the increasing use of rotating machinery in mechanical equipment, researchers have become more concerned about its performance requirements. They strive to improve the machines’ efficiency, reliability, and overall performance. However, the operating system of rotating machinery is often affected by various types of fault. Severe accidents can happen, resulting in economic losses. Therefore, it is necessary to monitor the status of the rotor system and detect early warning signs to avoid fatal damage.
Dynamic behavior can reflect the inherent condition of the structure system. Dynamic response changes can often be used as the point of structural damage detection. Previous studies showed that faults often introduce nonlinear dynamic characteristics to rotor systems; for example, the friction fault [4,5,6], misalignment fault [7,8,9], and crack fault [10,11] cause a multi-harmonic system output response other than the fundamental harmonic component. For this reason, the fault diagnosis method for rotor systems is often accompanied by a second harmonic vibration response. The more serious the misalignment fault is, the larger the proportion of the second harmonic [8,9]. Volterra series-based methods, such as GFRF (Generalized Frequency Response Function) and NOFRF (Nonlinear Output Frequency Response Function) [12], have also been used for the fault diagnosis of rotor systems [4,13,14]. This fault diagnosis method differs from observing the second harmonic, which only considers the output signal. NOFRFs or GFRFs are identified using the system outputs and inputs to ensure that the fault features are derived from the damage to the system rather than input excitation [4,15]. In addition, we propose a modified harmonic balance-alternating frequency/time domain method to analyze the nonlinear dynamics of a rotor system [16].
NOFRFs have been proposed as an alternative frequency domain representation of GFRFs, which circles GFRFs’ multi-dimensionality problem. They are often difficult to measure and display in practice. Current approaches to Volterra series-based fault diagnosis of rotor systems also have some problems. NOFRFs require multiple tests with different amplitudes at the same speed [13]. Furthermore, data-driven-based NOFRFs evaluated from the NARX model, created using input and output data from the system, are only applicable to conditions where the input and output data are random signals [15]. It possesses the abundant frequency and amplitude characteristics of a NARX model, which can test the different properties of a system except for rotor systems. On the other hand, the truncation order of NOFRFs depends on their presentation ability and the extent to which NOFRFs can distinguish the structural properties of different damage conditions. In other words, although the NOFRFs up to the fourth order were considered in the present research, errors can still be caused if it is not appropriate a priori [14].
To solve the problems of a NOFRF fault diagnosis method, we developed a NARX model for rotor systems from the multi-harmonic signal generated by a speed-up process. We also used NRSF (Nonlinear Response Spectrum Function) [17], a system frequency analysis tool based on the notion that the system frequency response function consists of a first-order frequency response function and other orders combined. Then, we developed a condition monitoring method through an index calculated by NRSF to quantify the severity of the nonlinear fault. The proposed method is more efficient than observing the second harmonic method in the following experiment.
In summary, there are two innovations in this paper with respect to state of the art technology: (1) using Time Synchronous Averaging (TSA) technology to improve the quality of input and output signals; (2) proposing NFIs (Nonlinear Fault Indicators) to quantify the severity of the fault.
This paper is organized as follows: Section 2 describes the theoretical background of the NARX model and the NRSFs method. The condition monitoring method for rotor systems is proposed in Section 3. Section 4 verifies the validity and feasibility of this method through an experiment. Finally, our conclusions are summarized in Section 5.

2. The NARX Model and NRSFs of a Nonlinear System

2.1. The NARX Model of a Nonlinear System

Compared to various data-driven models, NARX models are suited to condition-monitoring rotor systems due to their physically interpretable nature based on being mapped into a frequency function [18]. The NARX model represents a large class of input–output nonlinear dynamic systems [19].
y ( k ) = F [ y ( k 1 ) , , y ( k n y ) , u ( k 1 ) , , u ( k n u ) ] + e ( k )
where F [ ] is a nonlinear function that must be identified from the input u ( k ) and the output y ( k ) ; k represents the discrete time; n u and n y are the maximum input and output lags; and e ( k ) is the model prediction error assumed to be white noise.
Considering u ( k ) as a polynomial function, the NARX model (1) can be written as:
y ( k ) = θ 0 + i 1 = 1 n θ i 1 x i 1 ( k ) + i 1 = 1 n i 2 = i 1 n θ i 1 i 2 x i 1 ( k ) x i 2 ( k ) + + i 1 = 1 n i l = i l 1 n θ i 1 i l x i 1 ( k ) x i l ( k ) + e ( k )
where θ i 1 i l are the model coefficients, l being the degree of polynomial nonlinearity, n = n u + n y , and
x m ( k ) = { y ( k m ) 1 m n y u ( k m + n y ) n y + 1 m n y + n u
The polynomial NARX model (2) can also be represented as a linear-in-the-parameters model:
y ( k ) = i = 1 n M θ i p i ( k ) + e ( k )
where n M is the number of NARX model terms with full structure and n M = ( n + L ) ! / ( n ! L ! ) ; p i ( k ) , i = 1 , , n M are the monomials of the NARX model. Formula (2) is formulated by x 1 ( k ) , , x n ( k ) ; θ i , i = 1 , , n M . n M of the NARX model candidate is usually a high number. Identifying the NARX model helps select significant model terms from a large set of candidates based on input/output data so that a parsimony model structure can represent system dynamics.
For example, a system with a nonlinear differential equation model:
m y ¨ ( t ) + c y ˙ ( t ) + k 1 y ( t ) + k 3 y 3 ( t ) = u ( t )
where u ( t ) and y ( t ) are the input and output; t is continuous time; and m , k 1 , c , and k 3 are the mass, linear stiffness, damping, and nonlinear stiffness, respectively. Formula (4) can be discretized into a difference equation through a forward–backward difference approach, where y ˙ ( t ) and y ¨ ( t ) can be discretized as follows [20]:
y ˙ ( t ) = y ( t ) y ( t Δ t ) Δ t y ¨ ( t ) = y ( t + Δ t ) 2 y ( t ) + y ( t Δ t ) Δ t 2
where Δ t = 1 / f s is the time step and f s is the sampling frequency; therefore, (4) can be written as:
y ( k ) = 2 m c Δ t k Δ t 2 m y ( k 1 ) m c Δ t m y ( k 2 ) + Δ t 2 m u ( k 1 ) k 1 Δ t 2 m y 3 ( k 1 )
Let m = 1 kg , c = 12 Ns / m , k 1 = 150 N / m , and k 3 = 8 × 10 8 N / m 3 . Formula (6) becomes (7) by setting f s = 1024 Hz ,
y ( k ) = 1.9881 y ( k 1 ) 0.9883 y ( k 2 ) + 9.5367 × 10 7 u ( k 1 ) 762.9395 y 3 ( k 1 )
which can be considered as an identified NARX model of the Duffing system (5).
In practice, the NARX model identification typically uses the FROLS (Forward Regression Orthogonal Least Squares) algorithm by computing the ERR (Error Reduction Ratio) of each of the candidate terms [21], which will be introduced below.
The NARX model (4) can be written into a matrix form as:
Y = P Θ + e
where,
Y = [ y ( 1 ) , , y ( N ) ] T ; Θ = [ θ 1 , , θ n M ] T ; e = [ e ( 1 ) , , e ( N ) ] T ; P = [ p 1 , , p n M ] = [ p 1 ( 1 ) p n M ( 1 ) p 1 ( N ) p n M ( N ) ]
with p i = [ p i ( 1 ) , , p i ( N ) ] T , i = 1 , , n M . A matrix letting can be orthogonally decomposed as:
P = W A
where A is an N × n M unit upper triangular matrix W T W = diag [ d 1 , , d n M ] , with d i = w i T w i , i = 1 , , n M . Therefore,
Y = ( P A 1 ) ( A Θ ) + e = W g + e
where W = P A 1 , g = [ g 1 , , g n M ] T is an auxiliary parameter vector.
Then the ERR criterion is defined as
[ e r r ] i = ( w i T w i ) g i 2 Y T Y , i = 1 , , n M
The significant model term with the highest ERR value is selected by comparing the ERR values of each column vector of the regression matrix P . The process is repeated to determine the n ¯ M significant model terms until
1 i = 1 n ¯ M [ e r r ] i < ρ
where ρ is the threshold.
In addition, an identified NARX model can easily transform into the frequency response function, revealing core invariant behavior and facilitating the analysis, design, and fault diagnosis for the system [15].

2.2. The NRSFs Representation of Nonlinear System

When a discrete-time NARX model is identified, the N R S F 1 ( j ω ) [17], also known as the linear transfer function H 1 ( j ω ) of the NARX model, can be directly evaluated through harmonic probing method [22] in the discrete-time domain.
N R S F 1 ( j ω ) = H 1 ( j ω ) = N u = 1 n u θ N u exp ( N u j ω Δ t ) 1 N y = 1 n y θ N y exp ( N y j ω Δ t )
where N u and N y represent lags in the linear model terms.
The Volterra functional series is shown in [18]
y ( t ) = n = 1 + y n ( t )
where y n ( t ) is the nth order output of the system. Fourier transform Y ( j ω ) of y ( t ) can be expressed as:
Y ( j ω ) = n = 1 + Y n ( j ω )
where ω is the frequency and Y n ( j ω ) is the nth-order output frequency spectra obtained using Fourier transforms y n ( t ) .
Therefore,
Y ( j ω ) = Y 1 ( j ω ) + n = 2 + Y n ( j ω ) = H 1 ( j ω ) U ( j ω ) + Y ( j ω ) Y 1 ( j ω ) U ( j ω ) U ( j ω )
where,
{ N R S F 1 ( j ω ) = H 1 ( j ω ) N R S F ( j ω ) = Y ( j ω ) Y 1 ( j ω ) U ( j ω ) = Y ( j ω ) U ( j ω ) N R S F 1 ( j ω )
are defined as the NRSFs of nonlinear systems and N R S F ( j ω ) represents the system of nonlinear dynamics [17]. Based on the NARX model (8), N R S F 1 ( j ω ) can be evaluated as:
N R S F 1 ( j ω ) = 9.5367 × 10 7 exp ( j ω / 1024 ) 1 1.9881 exp ( j ω / 1024 ) + 0.9883 exp ( 2 j ω / 1024 )
Considering ω = [ 0 , 20 ]   r a d / s , Figure 1 shows the N R S F 1 ( j ω ) and N R S F ( j ω ) of the NARX model under normalized amplitude excitation in the frequency range.
System damage usually introduces nonlinear dynamic characteristics in the structure, causing nonlinear behaviors in a system that originally presented linear characteristics [23,24]. For example, when rub-impact happens to rotor systems, the vibration waveform in the time domain is distorted because of the superposition of higher frequency components on the fundamental harmonic [4,5,7]. The vibration response of misalignment faults mainly consists of the fundamental harmonic and the second harmonic. The more severe the misalignment, the larger the proportion of the second harmonic [7,9,10]. The crack fault in the rotor system can also cause fundamental, second, and third harmonic components [10,11].
When an alternative system frequency domain expression differs from GFRF and NOFRF, NRSF divides the system’s frequency response into linear and nonlinear parts without considering the GFRF dimension or the truncation order of NOFRF. The calculation results can be accurate and reliable with a higher calculation efficiency. Previous studies verified the sensitivity of NRSFs to changes in system parameters and system faults [17], whose results show that NRSFs can reflect the changes in system characteristics and provide effective feature extraction for system analysis, parameter estimation, and fault diagnosis.

3. Condition Monitoring Method for Rotor Systems

3.1. Problem Statement

The condition monitoring method in this paper can be summarized into four steps: signal collecting, data processing, data-driven modeling, and fault feature extraction. The NARX model identification and NRSF method mentioned in Section 2 can be used to represent the rotor system and feature extraction of the model. Figure 2 shows the basic idea behind the proposed method. By extracting nonlinear characteristics from the NARX model established by the FROLS algorithm based on the vibration signals, which represent changes in physical characteristics caused by failure or damage in the rotor system, and comparing the fault index with the healthy condition, we can determine system faults and their severity. However, there are still some related problems that must be solved before they can be used for condition monitoring of rotor systems in practical engineering.
Due to noise, the quality of the input and output signals significantly affects the NARX model’s ability to characterize the system. Therefore, the measured signals should be processed before being used to establish a model.
The basic idea of most fault diagnosis methods is to extract the fault features and compare them with the healthy state to judge the current system condition [15]. Therefore, a fault indicator can help conduct the above process more conveniently.

3.2. Solution of Existing Problems

The random input that a NARX model often applies does not apply to rotor systems, which often operate at several rotating speeds [1,2,3]. The NARX model established using single harmonic signals can only be investigated under the same frequency [25]. Many researchers have analyzed rotor-bearing systems based on the speed-up process [21] and the frequency sweep system identification approach [26]. In this method, the NARX model of rotor systems is established using vibration signals under different rotating speeds, which reflect the frequency characteristics of the rotor system within the speed range.
If the rotating speeds of rotor systems in engineering practice are limited, then the range of monitoring speed [ ω min , ω max ] should be selected when the collecting is convenient. Collect the displacement vibration signal x 1 ( t ) under multiple monitoring speeds within [ ω min , ω max ] at the same observation point, where i represents the i th monitoring rotating speed. Keep the observation point, collecting time t , and sampling frequency f s unchanged for each monitoring speed.
Many data processing methods can address the noise in NARX modeling and enhance the NARX model’s generalization ability, including regularization [27], AIC (Akaike Information Criterion), BIC (Bayesian Information Criterion), etc. [28]. Time Synchronous Averaging (TSA) technology [29] efficiently alleviates the influence of noise and measurement error on signals since the rotor system’s vibration response signals are periodic. Before signal division, the vibration response signals x 1 ( t ) can be resampled with a sampling frequency higher than the original sampling frequency, making it more accurate to divide the signal into segments with more complete periods. The data processing method using TSA technology is presented as follows.
Determine the number of periods N for each segment of data used for stacking and divide the signal into [ M / N ] segments, where M represents the total number of periods. Choose zero points at the start and end of each segment. Determine the sum of each segment of data and calculate the average as follows:
y i ( t ) = 1 [ M / N ] n = 1 [ M / N ] x i , n ( t )
where y i ( t ) is the processed data of the i th monitoring speed, x i , n ( t ) represents the n th segment of data for the i th monitoring speed.
Considering the duffing system (5), the output response under u ( t ) = 2 sin ( ω 1 t ) , where ω 1 = 15 rad / s , demonstrates the signal process method. To highlight improvements in the processed signal, a noise e ( t ) is mixed with the system output x 1 ( t ) by setting e ( t ) as a uniformly distributed random sequence in [ 4 , 4 ] × 10 4 . We collected 300s of the system output with a sampling frequency of f s = 1024 and obtained M = 714 periods to conduct the TSA method. Figure 3 shows the number of periods for each segment is N = 5 where an interval of two stars denotes a period. The signal is divided into [ M / N ] = 142 segments. After summation averaging (15), Figure 4 shows the output y 1 ( t ) of using TSA technology.
The rotor system is affected by a centrifugal force from an unbalanced mass. However, the centrifugal mass and distance are sometimes unknown or difficult to measure. Therefore, the input u i ( t ) for the NARX model can be generated in a harmonic form with a normalized amplitude.
u i ( t ) = sin ( ω i t )
NARX models can adapt to different system parameters by changing the coefficients of model terms so that the amplitude of the input will not affect its ability to represent the rotor system. At each rotating speed, the sampling frequency, vibration frequency, and data length of the input signal are consistent with the output signal, and the starting point is also set to zero. The sampling frequency can be appropriately reduced after data processing to create NARX models more efficiently.
In this case, several data sets under different monitoring speeds must be handled simultaneously and one common NARX model is required. This can be achieved by the data splicing method, where the inputs u i ( t ) and outputs y i ( t ) under different speeds are spliced into U ( t ) and Y ( t ) , respectively, in the same order. Then, the FROLS algorithm can be applied to establish NARX models.
To quantify the intensity of nonlinear faults and facilitate judgment of the rotor system condition, NFIs (Nonlinear Fault Indicator) are proposed in this method, which are defined as follows:
N F I = ω min ω max | N R S F ( j ω ) | d ω ω min ω max | N R S F ( j ω ) | d ω + ω min ω max | N R S F 1 ( j ω ) | d ω
where ω min and ω max are the minimum and maximum values of the frequency range, respectively. It is worth noting that NFI is a special case of NOFRF-based indices [30].
Under the investigation frequency ω n selected in [ ω min , ω max ] , we applied a single harmonic excitation with the normalized amplitude sin ( ω n t ) to the NARX model to calculate NRSF.
Compared to the fault features calculated at a single frequency, the index proposed in this method contains more information [15], resulting in a more reliable monitoring result.

3.3. Condition Monitoring Procedure

The condition monitoring method for rotor systems using NARX model and NRSF-based analysis can be summarized as follows:
Step 1: Collect the displacement vibration signals x i ( t ) under multiple monitoring rotating speeds ω i , where i represents the i th monitoring rotating speed;
Step 2: Process the collected signals using the TSA method and generate the corresponding input signal. Then, splice u i ( t ) and y i ( t ) into U ( t ) and Y ( t ) , respectively;
Step 3: Establish a NARX model to characterize the rotor system based on input U ( t ) and output Y ( t ) data using the FROLS algorithm;
Step 4: Calculate the nonlinear fault index (17) based on the identified NARX model and compare the results to healthy conditions.

4. Experimental Study

4.1. The Misaligned Rotor System Experiment

We conducted this experiment to demonstrate the condition monitoring method using an experimental rotor system with various levels of misalignment (Figure 5).
By setting different offset distances between the edge of the bearing seat and the pedestal, we achieved different levels of misalignment fault. Each increase of 1 mm offset distance was set as a misalignment level, and a total of four levels were set in the experiment, namely level 1, level 2, level 3, and level 4, respectively. The eddy current displacement sensor was aligned with the horizontal direction of the rotating shaft and connected to the LMS mobile front end to collect the experimental rotor system’s response.
Using the proposed method in Section 3.3, we performed condition monitoring of the rotor system under the level 1 misalignment fault.
Step 1: Collect the displacement vibration signals.
The displacement vibration signals were collected under rotor system speeds of ω 1 = 1550   r / m i n , ω 2 = 1650   r / m i n , and ω 3 = 1700   r / m i n at the same observation point, denoted as x 1 ( t ) , x 2 ( t ) , and x 3 ( t ) , respectively, with a sampling frequency f s 1 = 2048   Hz and collecting time of t = 33   s ;
Step 2: Process the collected signals using the TSA and data splicing methods.
The vibration response signal under ω 3 = 1700   r / m i n at the level 1 misalignment demonstrates data processing in this method. The collected vibration response signal x 3 ( t ) was interpolated by a cubic spline and sampled again with the new sampling frequency f s 2 = 20480   Hz , which was L = 10 times the original sampling frequency f s 1 . The number of periods for each segment was determined to be N = 20 and the whole M = 895 periods were divided into [ M / N ] = 44 segments. Each segment of data was summed together and an average was calculated. After the system outputs were processed, the inputs for NARX modeling were generated using (16).
The data splicing method was applied to handle the three data sets. Under the level 1 misalignment of the rotor system, the input and output data were sequentially spliced, and the data after splicing was resampled with the sampling frequency f s 3 = 1024   Hz . Figure 6 and Figure 7 show the splicing results of input U ( t ) and output Y ( t ) , respectively.
Step 3: NARX modeling using the FROLS algorithm.
By using the FROLS algorithm, we initialized the maximum time delays for the system input n u = 3 , output n y = 3 , and the degree of polynomial nonlinearity l = 3 . Based on the input U ( t ) and output Y ( t ) data after splicing, the NARX model for level 1 of the misaligned rotor system is:
y ( k ) = 1.051 y ( k 1 ) 1.896 × 10 1 y ( k 2 ) 2.708 × 10 1 y ( k 3 ) + 3.385 × 10 1 u ( k 1 ) 6.248 × 10 1 u ( k 2 ) + 5.447 × 10 3 u ( k 3 ) 2 + 3.146 × 10 2 u ( k 1 ) 2 + 7.014 × 10 3 u ( k 3 ) 3 + 3.047 × 10 1 u ( k 3 ) 4.513 × 10 2 u ( k 2 ) 2 + 1.834 × 10 1 u ( k 2 ) y ( k 1 ) + 4.819 y ( k 3 ) 2 8.429 × 10 1 u ( k 3 ) 2 y ( k 1 ) + 31.137 u ( k 1 ) y ( k 3 ) 2 20.021 u ( k 3 ) y ( k 2 ) y ( k 3 ) 5.318 y ( k 1 ) 2
Figure 8 shows the MPO (Model Predicted Output) method’s predicted output from the NARX model (24) and the real output of level 1 of the misaligned rotor system. They are almost identical, which validates the NARX model (24).
Step 4: Calculate the NRSF-based nonlinear fault index.
The frequency domain characteristics of the NARX model were investigated in the range of monitoring rotating speeds [ 1550 , 1700 ]   r / min . The investigation frequencies ω n used to calculate NRSF were set at every other 0.5   rad / s within the range [ 1550 , 1700 ]   r / min . Under the investigation frequency, a single harmonic excitation was applied to NARX models to calculate NRSF. Figure 9 shows NFI (17) calculations for the NARX model of the rotor system in the first misalignment level.
Using the same method as steps 1–4, we obtained the NFI of the rotor system in levels 2–4 of misalignment fault (Figure 8). As shown in Figure 9, the NFI values at four misalignment levels are 0.076, 0.21, 0.25, and 0.45, respectively. These values have an increasing monotonic trend with the severity of the misaligned fault, meaning that the index NFI’s value can directly reflect the severity of the misaligned fault. Therefore, the method proposed in this paper accurately identifies rotor system faults.

4.2. Comparison to the Traditional Method

Misalignment faults are often accompanied by second harmonic vibration responses. The more severe the misalignment fault is, the larger the proportion of the second harmonic [8,9]. Therefore, the amplitude of the second harmonic was always used to characterize misalignment faults in rotor systems. We compared methods to verify the advantages of the method proposed in this paper.
The vibration response spectrums of each misaligned rotor system level were calculated using the experimental data investigated beforehand. Figure 10 shows the amplitude of the second harmonic. In the traditional method, there is no obvious upward monotonic trend with the severity of the misaligned fault. For example, the amplitude of the second harmonic is gradually increasing with the severity of the misaligned fault in levels 1–3 but decreases at level 4.
The diagnosis method based on output signals is affected by system input and noise and cannot reflect the essential characteristics of the system. By contrast, we propose a model-based fault feature extraction method with an index calculated in a multi-frequency range, which is more reliable for condition-monitoring rotor systems.

5. Conclusions

In this study, we propose a NARX model-based NRSR method to facilitate condition monitoring for rotor systems. We used TSA technology to improve the quality of training data for NARX modeling. We also used an NFI indicator to quantify fault severity in rotor systems. Our specified conclusions are as follows:
(1)
NRSF-based methods are more suitable for situations where the prior conditions are insufficient and the fault characteristics are unknown since they do not have to consider the truncation order of NOFRF;
(2)
The NFI values at four misalignment levels were 0.076, 0.21, 0.25, and 0.45, which have a monotonic trend with the severity of the misaligned fault validating the proposed method’s efficiency over the second harmonic method.

Author Contributions

Conceptualization, Y.G. and Y.-P.Z.; methodology, Y.-P.Z.; software, Y.G. and C.Y; validation, Y.G. and C.Y.; resources, Z.L.; data curation, Y.G.; writing—original draft preparation, Y.G.; writing—review and editing, C.Y.; project administration, Z.L.; funding acquisition, Y.-P.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of Liaoning Province, grant number 2023-MS-031.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, X.; Wang, X.; Liu, R.; Wang, B. Influence of temperature on nonlinear dynamic characteristics of spiral-grooved gas-lubricated thrust bearing-rotor systems for microengine. Tribol. Int. 2013, 61, 138–143. [Google Scholar] [CrossRef]
  2. Wang, Y.; Zheng, Q.; Fu, D.; Zhang, H. Study on adaptive torsional vibration suppression methods for helicopter/turboshaft engine system with variable rotor speed. Asian J. Control 2021, 23, 1490–1502. [Google Scholar] [CrossRef]
  3. Misté, G.A.; Garavello, A.; Benini, E.; Gonzalez-Alcoy, M. A Methodology for Determining the Optimal Rotational Speed of a Variable RPM Main Rotor/Turboshaft Engine System. J. Am. Helicopter Soc. 2015, 60, 1–11. [Google Scholar] [CrossRef]
  4. Yao, H.; Han, Q.; Li, L.; Wen, B. Detection of rubbing location in rotor system by super-harmonic responses. J. Mech. Sci. Technol. 2012, 26, 2431–2437. [Google Scholar] [CrossRef]
  5. Zhang, W.; Meng, G.; Chen, D.; Zhou, J.; Chen, J. Nonlinear dynamics of a rub-impact micro-rotor system with scale-dependent friction model. J. Sound Vib. 2008, 309, 756–777. [Google Scholar] [CrossRef]
  6. Hua, C.; Rao, Z.; Ta, N.; Zhu, Z. Nonlinear dynamics of rub-impact on a rotor-rubber bearing system with the Stribeck friction model. J. Mech. Sci. Technol. 2015, 29, 3109–3119. [Google Scholar] [CrossRef]
  7. Li, H.; Chen, Y.; Hou, L.; Zhang, Z. Periodic response analysis of a misaligned rotor system by harmonic balance method with alternating frequency/time domain technique. Sci. China Technol. Sci. 2016, 59, 1717–1729. [Google Scholar] [CrossRef]
  8. Patel, T.H.; Darpe, A.K. Experimental investigations on vibration response of misaligned rotors. Mech. Syst. Signal Process. 2009, 23, 2236–2252. [Google Scholar] [CrossRef]
  9. Sekhar, A.S.; Prabhu, B.S. Effects of coupling misalignment on vibrations of rotating machinery. J. Sound Vib. 1995, 185, 655–671. [Google Scholar] [CrossRef]
  10. Peng, Z.K.; Lang, Z.Q.; Billings, S.A. Crack detection using nonlinear output frequency response functions. J. Sound Vib. 2007, 301, 777–788. [Google Scholar] [CrossRef]
  11. Xiang, L.; Gao, X.; Hu, A. Nonlinear dynamics of an asymmetric rotor-bearing system with coupling faults of crack and rub-impact under oil-film forces. Nonlinear Dyn. 2016, 86, 1057–1067. [Google Scholar] [CrossRef]
  12. Cheng, C.M.; Peng, Z.K.; Zhang, W.M.; Meng, G. Volterra-series-based nonlinear system modeling and its engineering applications: A state-of-the-art review. Mech. Syst. Signal Process. 2017, 87, 340–364. [Google Scholar]
  13. Liu, Y.; Zhao, Y.L.; Li, J.T.; Lu, H.H.; Ma, H. Feature extraction method based on NOFRFs and its application in faulty rotor system with slight misalignment. Nonlinear Dyn. 2020, 99, 1763–1777. [Google Scholar] [CrossRef]
  14. Lu, H.; Feng, K.; Liang, H.; Liu, Y.; Li, J.; Zhao, Y. An improved NOFRFs-based fault feature extraction method and its application to quantitative diagnosis in rotor rub-impact. J. Sound Vib. 2021, 513, 116406. [Google Scholar] [CrossRef]
  15. Peng, Z.K.; Lang, Z.Q.; Wolters, C.; Billings, S.A.; Worden, K. Feasibility study of structural damage detection using NARMAX modelling and Nonlinear Output Frequency Response Function based analysis. Mech. Syst. Signal Process. 2011, 25, 1045–1061. [Google Scholar] [CrossRef]
  16. Chen, Y.; Hou, L.; Song, H.; Lin, R.; Jin, Y.; Chen, Y. Nonlinear dynamics analysis of a dual-rotor-bearing-casing system based on a modified HB-AFT method. Mech. Syst. Signal Process. 2023, 185, 109805. [Google Scholar] [CrossRef]
  17. Ma, Y.; Luo, Z.; Liu, H.P.; Zhu, Y.P. The NRSF-SVM based method for nonlinear rotor bearing fault diagnosis. In Proceedings of the 30th Chinese Control and Decision Conference, Shenyang, China, 9–11 June 2018. [Google Scholar]
  18. Billings, S.A. Nonlinear System Identification: NARMAX Methods Nonlinear System Dentification: NARMAX Methods in the Time, Frequency, and Spatio-Temporal Domains; John Wiley and Sons: Chichester, UK, 2013; pp. 337–367. [Google Scholar]
  19. Leonyaritis, I.J.; Billings, S.A. Input-output parametric models for non-linear systems Part I: Deterministic non-linear systems. Int. J. Control 1985, 41, 303–328. [Google Scholar] [CrossRef]
  20. Ge, X.; Luo, Z.; Ma, Y.; Liu, H.P.; Zhu, Y.P. A novel data-driven model based parameter estimation of nonlinear systems. J. Sound Vib. 2019, 453, 188–200. [Google Scholar] [CrossRef]
  21. Ma, Y.; Liu, H.; Zhu, Y.; Wang, F.; Luo, Z. The NARX Model-Based System Identification on Nonlinear, Rotor-Bearing Systems. Appl. Sci. 2017, 7, 911. [Google Scholar] [CrossRef]
  22. Worden, K.; Manson, G.; Tomlinson, G.R. A Harmonic Probing Algorithm for The Multi-input Volterra Series. J. Sound Vib. 1997, 201, 67–84. [Google Scholar] [CrossRef]
  23. Peng, Z.k.; Lang, Z.Q. The Nonlinear Output Frequency Response Functions of One-Dimensional Chain Type Structures. J. Appl. Mech. 2010, 77, 1007. [Google Scholar] [CrossRef]
  24. Peng, Z.K.; Lang, Z.Q.; Billings, S.A. Analysis of Locally Nonlinear MDOF Systems Using Nonlinear Output Frequency Response Functions. J. Vib. Acoust. 2009, 131, 051007. [Google Scholar] [CrossRef]
  25. Zhu, Y.; Zhao, Y.; Lang, Z.Q.; Liu, Z.; Liu, Y. On-line Rotor Systems Condition Monitoring Using Nonlinear Output Frequency Response Functions under Harmonic Excitations. IEEE Trans. Ind. Inform. 2022, 18, 6798–6808. [Google Scholar] [CrossRef]
  26. Li, Y.Q.; Luo, Z.; He, F.X.; Zhu, Y.P.; Ge, X.B. Modeling of rotating machinery: A novel frequency sweep system identification approach. J. Sound Vib. 2021, 494, 115882. [Google Scholar] [CrossRef]
  27. Chen, S.; Chng, E.S.; Alkadhimi, K. Regularized orthogonal least squares algorithm for constructing radial basis function networks. Int. J. Control 1996, 64, 829–837. [Google Scholar] [CrossRef]
  28. Gu, Y.; Wei, H.L.; Balikhin, M.M. Nonlinear predictive model selection and model averaging using information criteria. Syst. Sci. Control Eng. 2018, 6, 319–328. [Google Scholar] [CrossRef] [Green Version]
  29. Combet, F.; Gelman, L. An automated methodology for performing time synchronous averaging of a gearbox signal without speed sensor. Mech. Syst. Signal Process. 2007, 21, 2590–2606. [Google Scholar] [CrossRef] [Green Version]
  30. Yu, C.S.; Zhu, Y.P.; Luo, H.T.; Luo, Z. Structural damage detection of spacecraft equipment based on data driven modelling and nonlinear output frequency response function. In Proceedings of the 2022 4th International Conference on Industrial Artificial Intelligence (IAI), Shenyang, China, 24–27 August 2022. [Google Scholar]
Figure 1. The NRSFs of the NARX model (8) (a) N R S F 1 ( j ω ) ; (b) N R S F ( j ω ) .
Figure 1. The NRSFs of the NARX model (8) (a) N R S F 1 ( j ω ) ; (b) N R S F ( j ω ) .
Sensors 23 06878 g001
Figure 2. The basic idea of the proposed method.
Figure 2. The basic idea of the proposed method.
Sensors 23 06878 g002
Figure 3. The system output x 1 ( t ) with noise e ( t ) .
Figure 3. The system output x 1 ( t ) with noise e ( t ) .
Sensors 23 06878 g003
Figure 4. The TSA processed result of y 1 ( t ) .
Figure 4. The TSA processed result of y 1 ( t ) .
Sensors 23 06878 g004
Figure 5. Experimental setup.
Figure 5. Experimental setup.
Sensors 23 06878 g005
Figure 6. Splicing results U ( t ) of level 1 of the misaligned rotor system.
Figure 6. Splicing results U ( t ) of level 1 of the misaligned rotor system.
Sensors 23 06878 g006
Figure 7. Splicing results Y ( t ) of level 1 of the misaligned rotor system.
Figure 7. Splicing results Y ( t ) of level 1 of the misaligned rotor system.
Sensors 23 06878 g007
Figure 8. Fitting result of the first level of misaligned rotor system.
Figure 8. Fitting result of the first level of misaligned rotor system.
Sensors 23 06878 g008
Figure 9. Fitting results of the first level of misaligned rotor systems.
Figure 9. Fitting results of the first level of misaligned rotor systems.
Sensors 23 06878 g009
Figure 10. Results of the traditional method.
Figure 10. Results of the traditional method.
Sensors 23 06878 g010
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gao, Y.; Yu, C.; Zhu, Y.-P.; Luo, Z. A NARX Model-Based Condition Monitoring Method for Rotor Systems. Sensors 2023, 23, 6878. https://doi.org/10.3390/s23156878

AMA Style

Gao Y, Yu C, Zhu Y-P, Luo Z. A NARX Model-Based Condition Monitoring Method for Rotor Systems. Sensors. 2023; 23(15):6878. https://doi.org/10.3390/s23156878

Chicago/Turabian Style

Gao, Yi, Changshuai Yu, Yun-Peng Zhu, and Zhong Luo. 2023. "A NARX Model-Based Condition Monitoring Method for Rotor Systems" Sensors 23, no. 15: 6878. https://doi.org/10.3390/s23156878

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