Next Article in Journal
Accumulative Roll Bonding—A Review
Previous Article in Journal
Selected Scattering on Quasi-Ordered Hexagonal Close-Packed Al Nanodents for Tunable Output of White LEDs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Dynamic Characteristic of a Faulted Rotor System with Multi-Objective Optimization Designed SFD

1
College of Aeronautical Engineering, Civil Aviation University of China, Tianjin 300300, China
2
State Key Laboratory of Engines, Tianjin University, Tianjin 300072, China
3
Renai College, Tianjin University, Tianjin 301636, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(17), 3628; https://doi.org/10.3390/app9173628
Submission received: 31 July 2019 / Revised: 19 August 2019 / Accepted: 19 August 2019 / Published: 3 September 2019
(This article belongs to the Section Mechanical Engineering)

Abstract

:
As an important part of the turbomachinery, the rotor–bearing system has been upgraded to provide a high rotating speed in order to meet the demand of high power production. With increasing demand for stability, the squeeze film damper (SFD) has been widely used in industrial machinery because it can reduce the vibration amplitude and suppress the external force. Usually, it shows inadaptability under the different working conditions where the SFD parameters didn’t change appropriately. Therefore, the reasonable choice of operational parameters of SFD is the key solution that can provide viscous damping effectively and restrain the nonlinear vibration generated by faults. In this paper, the mathematical model of a rotor-ball bearing-SFD system considering the misalignment fault and misalignment-rubbing coupling fault is built first. Then the dynamic characteristics under typical working conditions (ω = 1000 rad/s) of the faulted rotor are discussed. The vibration attenuation effects of the SFD parameters selected by using the multi-objective optimization method on the dynamic responses are analyzed. The results show that when the rotor system operates under different states, the value and the sensitivity of optimization parameters are altered. With no fault, the amplitude of fundamental frequency decrease 23%. With the misalignment fault, the amplitude of the fundamental frequency decreases by 43.4%, the amplitude of 2× fundamental frequency decreases by 27.5%, and the amplitude of 3× fundamental frequency decreases by 66.7%. With the misalignment-rubbing coupling fault, the amplitude of fundamental frequency reduces by 7.4%, the amplitude of 2× fundamental frequency drops by 51.5%, and the amplitude of 3× fundamental frequency drops by 16.8%. Overall, the feasibility of the optimization method of the variable-structured SFD operational parameters for the faulted rotor system is verified. These parametric analyses are very helpful in the development of a high-speed rotor system and provide a theoretical reference for the vibration control and optimal design of rotating machinery.

1. Introduction

The rotor bearing system of the turbomachinery, such as turbines, pumps and compressors, might undergo nonlinear supported forces and fault excitations. With increasing demand for vibration and stability, the squeeze film damper has been widely used in industrial machinery [1].
In 1977, based on Reynolds equation of short bearing assumption, Humes [2] studied the motion law of rigid rotor system with SFD, and the theoretical model of the hydrodynamic oilfilm force was proposed in 1978. Meanwhile, the obtained results were in agreement with the experimental results; when Coorper [3] presented the experimental study on SFD, and applied it to the Conway engine successfully, the vibration amplitude of the system was reduced by 60%. Considering the effect of the cavitation, Ruggiero et al. [4] found the “oil whirl” produced by the non-linear fluid forces. Wang [5] studied a controllable SFD for a rotor system and revealed its effect mechanism. Chang [6,7,8] analyzed the dynamic behavior of a hybrid SFD-mounted gear–bearing system based on the typical unbalance approximation and discussed the rotating speed and damping effects. Inayat-Hussain [9,10] presented the dynamic analysis of a flexible rotor mounted in SFD with and without retaining springs independently, and the results showed that the onset of bifurcation increased with the stiffness of the retainer spring, which provided insights into the effect of the design parameters of SFD on the rotor response. Qin [11] derived the dynamic equations of the single-disk rotor -SFD system under the elastic support and studied the nonlinear characteristics of the vibration responses. Consequently, there is no doubt that the application of SFD in rotor system can enhance the operational stability to a large extent. However, the improvement not only depends on the advanced structure design but also relates to the correlation of SFD parameters and the real condition of the rotor system. In other words, it is important to find out a suitable parameter set of SFD to prevent system instability due to fault occurrence as well as to facilitate product design and vibration control.
In this paper, a dynamic system model of the fault rotor system with SFD is established first. Then the dynamic characteristics under typical working conditions of the faulted rotor are discussed. And the vibration attenuation effects of the oil film gap (CS), elastic support stiffness (Ka) and the lubricating oil viscosity (μs) of SFD selected by using the multi-objective optimization method on the dynamic responses are analyzed.

2. The Basic Theory

2.1. Rotor System

The simplified model of a SFD-rotor-bearing system supported by deep groove ball bearing is shown in Figure 1. Rotor equivalent lumped mass points are Ms1, Md1, Ms2, Ms3, Ms4, Ms5, Md2, Ms6 respectively, Mo1 and Mo2 are lumped mass points of rolling bearing outer ring. K is stiffness of the flexible shaft; Ka is the supporting stiffness of squirrel cage, Cz and Cp are damping coefficients of rotor in bearing and disk. δ is the gap between rotor and stator; e1 and e2 are the eccentricities of the disk; α1 and α2 are the eccentric angles of the disk; ks is the stiffness coefficient.
According to the Lagrange equation, the differential equation of motion of rotor system can be obtained as follow:
[ M ] { q ¨ } + [ C ] { q ˙ } + [ K ] { q } = { Q } ,
In the formula, [M], [C], [K] are the mass matrix, damping matrix and stiffness matrix, respectively. Their details have been further described in Appendix A.
The system non-dimensional transformation is introduced: τ = ω t , { q ¯ } = { q } μ , { q ¯ . } = { q . } μ ω , { q ¯ ¨ } = { q ¨ } μ ω 2 , where ω is rotational angular velocity of the shaft; μ is the clearance of rolling bearing. Formula (1) transforms into non-dimensional equation:
{ q ¯ ¨ } = [ M ] 1 [ C ] ω { q ¯ ˙ } [ M ] 1 [ K ] ω 2 { q ¯ } + [ M ] 1 { Q ¯ } μ ω 2 .

2.2. Rub-Impact Model

Assuming that the rotor is in a state of rubbing, the friction between the rotor and the stator is in accordance with the Coulomb’s law [12].
Set e = x p 2 + y p 2 as radial relative displacement between the rotor disk and the casing, δ is the gap between the rotor and stator at rest. When e < δ , there is no rub-impact phenomenon. When e δ , rub-impact occurs, normal force F N and Tangential friction force F T can be expressed as:
{ F T = f F N F N = k s ( e δ ) ( e δ ) .
In the formula: k s is impact stiffness coefficient, f is friction coefficient.
The rub-impact forces are decomposed in the x y coordinate system:
{ F p x F p y } = F N e [ 1 f f 1 ] { X p Y p } .

2.3. Gear Coupling Model

Figure 2 shows a schematic diagram of the gear coupling. When the gear teeth are meshing tightly to transmit a large torque, deformation of gear teeth is inevitable and the following meshing force is related to the deformation amount, meshing distance and meshing rigidity [13]. The torque-induced meshing force is given by
F T i = ϕ L i K L i ,
ϕ = T / i = 1 z [ L i K L i ( R + L i ) ] ,
where Ф is the deformation angle of the gear tooth, T is torque, Li and KLi are the equivalent meshing distance and stiffness of each gear tooth respectively, R is the radius of the root circle.
When the torque is getting transmitted by gear teeth, a relative movement between the two meshing teeth takes place, which is induced by the system vibration at the same time. As a result, the displacement of meshing nodes causes a deformation of each tooth so that another meshing force is produced, which is expressed as:
F D i = ( e sin φ i ) K L i ,
φ i = 2 π ( i 1 ) / z ,
e = x 2 + y 2 .
In the formula, φi is the position angle of the gear tooth, and e′ is the misalignment of gear coupling.
From the above analysis, it can be concluded that the meshing force of coupling misalignment can be obtained as:
F L i = ( ϕ L i + e sin φ i ) K L i .
Then, the reactions of meshing force in the X and Y directions of the rotor system are as follows:
F L X = i = 1 z ( ϕ L i + e sin φ i ) K L i cos θ i ,
F L Y = i = 1 z ( ϕ L i + e sin φ i ) K L i sin θ i .

2.4. Squeeze Film Damper

The oil film force of SFD can be calculated according to the theory of hydrodynamic lubrication. The instantaneous pressure distribution can be computed from the incompressible Reynolds equation [14]:
1 R s 2 θ ( h 3 p θ ) + h 3 2 p z 2 = 12 μ s Ω h θ + 12 μ h t .
where, p is the pressure distribution of oil film, Rs is the radius of the damper, µs is the oil viscosity, h is the oil film thickness, z is the axial displacement in the z direction, Ω is precession angular velocity of the axis, and θ is the angular displacement from the maximum oil thickness position.
The oil film bearing forces can be obtained by using short bearing theory and the Sommerfeld transform, expressed as:
{ F sx = μ s R s L 3 c 2 x 2 + y 2 [ x ( ε ˙ I 1 + ε Ω I 2 ) y ( ε ˙ I 2 + ε Ω I 3 ) ] F sy = μ s R s L 3 c 2 x 2 + y 2 [ y ( ε ˙ I 1 + ε Ω I 2 ) + x ( ε ˙ I 2 + ε Ω I 3 ) ] .
and:
I 1 = θ 1 θ 1 + π cos 2 θ ( 1 + ε cos θ ) 3 d θ ,   I 2 = θ 1 θ 1 + π sin θ cos θ ( 1 + ε cos θ ) 3 d θ , I 3 = θ 1 θ 1 + π sin 2 θ ( 1 + ε cos θ ) 3 d θ
where ε = e B c is the journal eccentricity ratio, e B is the eccentricity of the journal and ε ˙ is the first derivative of journal eccentricity ratios with respect to time.

2.5. Rolling Bearing Force Model

In the process of rotation, the varying compliance vibration is caused by rolling bearing stiffness changing periodically. Based on the Hertzian contact theory [15], the supporting forces are given as:
{ F b x = i = 1 n C b [ ( x cos θ i + y sin θ i ) cos α μ ] 3 / 2      H [ ( x cos θ i + y sin θ i ) cos α μ ] cos θ i F b y = i = 1 n C b [ ( x cos θ i + y sin θ i ) cos α μ ] 3 / 2      H [ ( x cos θ i + y sin θ i ) cos α μ ] sin θ i ,
θ i = 2 π n ( i 1 ) + ω c a g e t ,   ω c a g e = ( R i R i + R o ) ω
In the formula, Fbx, Fby are supported forces of rolling bearing in the X and Y directions respectively, Cb is the Hertz contact stiffness, n is the number of rolling element, α is contact angle of rolling bearing, θi is the position angle of the ith ball, ωcage is the angular velocity of the cage, Ri, Ro are the radius of the inner and outer ring raceway respectively, μ is the bearing clearance, H is the heaviside function.

2.6. The Optimal Latin Hyercube Design

The Latin hyercube design is used to divide each dimension coordinate interval [ x k min , x k max ] , k [ 1 , n ] in the n-dimensional space into m intervals, each interval is x k i 1 , x k i , i [1,m]. Randomly selected m points, and ensured that each level of a factor has been studied only once, leading to an n-dimensional space, the sample number m of Latin hypercube design, marked as mxn LHD [16]. The optimal Latin hypercube design is improved at the foundation of random Latin hypercube design, so all the test points are distributed as evenly as possible in the design space, so as to further enhance its space filling and equalization. Figure 3 shows test points generated by random Latin hypercube design and test points generated by optimal Latin hypercube design, respectively.

2.7. Non-Dominated Sorting Genetic Algorithm with Elite Strategy (NSGA-II)

At present, the frequently used multi-objective algorithm is mainly divided into two kinds: one method is based on non-Pareto, the other is based on Pareto. For the method based on non-Pareto, such as a genetic algorithm VEGA based on vector evaluation, because of its high probability of local optimization, the performance still needs to be improved. While based on Pareto, NSGA-II is the most effective, it has improved at the foundation of the first-generation of NSGA [17]. The mian process is shown in Figure 4.

3. Multi Objective Optimization of SFD Parameters

3.1. Initial Scheme and Experimental Design of SFD

SFD can effectively reduces the vibration of the rotor system and restrain the nonlinearity caused by the misalignment coupling fault, in this paper, a variable-structured SFD is designed, which is shown in Figure 5, the SFD oil film gap CS and elastic support stiffness Ka are changed by adjusting the relative position of moving ring and length of elastic rod. The oil viscosity μs also can be changed with the help of Magnetorheological Fluid. According to the structural characteristics of SFD, to realize the optimization of SFD structure parameters, and taking the SFD oil film gap CS, elastic support stiffness Ka and oil viscosity μs as variables, the changing range of various parameters are shown in Table 1.
First, according to the changing range of the above mentioned variables and the optimal Latin hypercube design, the design matrix table of each parameter of SFD is obtained, and a total of 80 groups is obtained, while the design matrix of the SFD parameters is shown in Table 2.
Second, according to the optimized parameters of SFD in the design matrix of above experiment, a dynamic model of the SFD rotor rolling bearing system under three states is established: normal, misalignment fault, and the misalignment-rubbing coupling fault, while the model equations of motion are solved using the fourth-order Runge–Kutta method. The model parameters of the rotor shaft, SFD, ball bearing, rub-impact, and spline coupling are listed in Appendix A. Finally, the design parameters of SFD in Table 1 are used as inputs, while the calculated frequency amplitude is used as output, introducing into RBF neural network model, so as to establish an approximate model of SFD design parameters on the system frequency amplitude. Meanwhile, 10 sets of design parameters are selected and calculated by using the optimal Latin hypercube method, and the calculated results and prediction results of an approximate model are compared.

3.2. Single Objective Parameter Optimization Analysis of the Normal Rotor System

Figure 6 shows the system response characteristic of the system in the state of no fault, it can be found that spectrum diagram mainly manifested as the fundamental frequency vibration caused by the unbalance eccentricity. The motion of rotor is an ellipse, which means the system completes a single-period movement.
Because of the vibration energy of normal system is mainly concentrated in the fundamental frequency, the optimization objective function is aiming at the minimum of its amplitude.
Figure 7 shows the calculation results of the 10 sets of normal system design parameters selected by the optimal Latin hypercube method and the comparison of the predicted results with the RBF approximation model. Because the fundamental frequency which caused by the rotating speed is the mainly response, which can get the average error value of the amplitude of one time frequency is 0.9%, the maximum error value is 2.7%. The error of the approximate model output variables can be controlled within 10%, which shows that the approximate model has high precision, and can be used to optimize the SFD parameters of the normal rotor system.
Figure 8 shows the optimization parameters of oil film gap CS, elastic support stiffness Ka, and the influence of the oil viscosity of μs on the amplitude of the fundamental frequency. The effect of three parameters to the system is strongly nonlinear, the CS has the most obvious effect.
Sequential quadratic programming method (NLPQL) is introduced to optimize SFD parameters of normal system, the objective function is the amplitude of the fundamental frequency. The optimal solution for the design parameters of SFD can be solved as: CS = 4 × 10−5; Ka = 8,008,000; μs = 0.038.
According to the design parameters, the numerical simulation of the normal SFD- rotor system is made, and the dynamic response spectrum of the system is obtained. After comparing Figure 9 and Figure 6, it can be seen that, after optimization, the amplitude of the fundamental frequency of the system dynamics response decreases from 0.8913 to 0.6858, while the amplitude is reduced by nearly 23%.

3.3. Analysis Results of Optimization Parameters of Rotor System with the Misalignment Fault

In this section, the SFD optimization parameters of the rotor system with no fault are used as the initial calculated parameters of the SFD of rotor system under misalignment fault, to analyze the dynamic characteristics of the rotor system with misalignment, in which the misalignment value is e0 = 3 × 10−4 m.
Figure 10 shows the response characteristic of the system under the misalignment fault, from Figure 10a, it can be found that the fault makes the system produce 2× fundamental frequency, 3× fundamental frequency and other high frequency components, the VC frequency of rolling bearing has growth, while the frequency spectrum is mainly fundamental frequency, 2× fundamental frequency, and 3× fundamental frequency. The calculated results are in agreement with those in reference [12]. From Figure 10b, the phase diagram of the system is shown as a net 8-shaped figure, which means the system is in quasi periodic motion.
The energy of misalignment fault system is mainly concentrated in fundamental frequency, 2× fundamental frequency, 3× fundamental frequency components, the amplitude of fundamental frequency, 2× fundamental frequency, 3× fundamental frequency of the optimization objective function for the system response is minimal.
Figure 11 is a comparison of the calculated results of 10 sets of design parameters which were selected by using the optimal Latin hypercube method under misalignment fault and the predictive results of approximate model which is built by using the radial basis function (RBF). In this figure, the average error of fundamental frequency is 1.2%, and the maximum error is 6.8%; the average error of 2× fundamental frequency is 1.8%, and the maximum error is 8.4%; the average error of 3× fundamental frequency is 0.6%, and the maximum error is 2.1%. The error of the approximate model output variables can be controlled within 10%, which shows that the approximate model has high accuracy and that it can be used to optimize the SFD parameters of the rotor system with misalignment-rubbing coupling faults.
Figure 12 shows the optimization parameters of the oil film gap CS, elastic support stiffness Ka and the impact of the viscosity of the lubricating oil μs on the amplitude of fundamental frequency, 2× fundamental frequency and 3× fundamental frequency. From Figure 12, it can be seen that the effect of the three parameters on the system is nonlinear. As a result, fundamental frequency, 3× fundamental frequency mainly influences the oil film gap CS; and the 2× fundamental frequency is mainly influenced by the support stiffness Ka.
The SFD parameters of rotor system under the misalignment fault is optimized by using NSGA-II, the parameter configuration of the algorithm is as follows: the population size is 40, algebra is 200, the crossover rate is 0.9, cross distribution index is 10, variation distribution index is 20. On these grounds, the optimal solution of SFD is CS = 2 × 10−4; Ka = 6 × 105; μs = 5 × 10−3.
According to the design parameters, the numerical simulation of the SFD- rotor system is carried out again, and the dynamic response spectrum of the system is obtained, as shown in Figure 13, and compared with the spectrum diagram in Figure 10a before optimization. It can be concluded that after optimization, the amplitude of fundamental frequency of system dynamics response decreases from 0.6916 to 0.3913, the amplitude decreases by about 43.4%; the amplitude of 2× fundamental frequency decreases from 0.0766 to 0.0555, the amplitude decreases by about 27.5%; the amplitude of 3× fundamental frequency decreases from 0.2591 to 0.0863, and the amplitude decreases by about 66.7%.

3.4. Analysis Results of Optimization Parameters of Rotor System with the Misalignment-Rubbing Coupling Fault

In this section, the SFD optimization parameters of the rotor system with the fault of misalignment referred in front section are used as the initial calculated parameters of the rotor system SFD, and the dynamic characteristics of the rotor system with misalignment-rubbing coupling faults are analyzed.
Figure 14 shows the system response characteristic of the system with misalignment-rubbing coupling faults. The fundamental frequency, 2× fundamental frequency, and 3× fundamental frequency are mainly shown in spectrum. Compared with a spectrum diagram with a misalignment fault, the rubbing coupling faults made the amplitude of fundamental frequency of the system increase a lot, and made the amplitude of the 2× fundamental frequency, 3× fundamental frequency change a small amount, and brought about an emergency of the continuous spectrum. However, the impact is modest, and the energy is low, meaning the calculated result obtained in this paper is in agreement with the experimental result of reference [18]. As seen in Figure 14a,b the phase diagram is a net 8-shaped figure, which means the system is in quasi periodic motion.
The energy of the misalignment-rubbing coupling faults is mainly concentrated in fundamental frequency, 2× fundamental frequency, and 3× fundamental frequency components, thus the amplitude of the fundamental frequency, 2× fundamental frequency, and 3× fundamental frequency of optimization objective function of the system responses is minimal.
Figure 15 is a comparison of the calculated results of 10 sets of design parameters which selected by using optimal Latin hypercube method under misalignment fault and the predictive results of approximate model which is built by using the radial basis function (RBF). Here, the average error of fundamental frequency is 0.56%, and the maximum error is 1.44%; the average error of 2× fundamental frequency is 1.47%, and the maximum error is 4.10%, the average error of 3× fundamental frequency is 0.47%, and the maximum error is 1.2%. The error of the approximate model output variables can be controlled within 10%, which shows that the approximate model has high accuracy and can be used to optimize the SFD parameters of the rotor system with misalignment-rubbing coupling faults.
Figure 16 shows the optimization parameters of the oil film gap CS, elastic support stiffness Ka and the impact of the viscosity of the lubricating oil μs on the amplitude of the fundamental, 2× fundamental and 3× fundamental frequency. From Figure 16, it can be seen that the impact of the three parameters on the system is nonlinear. Here, CS has a great influence on 2× fundamental frequency, and basically it increases linearly in the early stage. For elastic support stiffness Ka, when its value increases, the amplitude of fundamental and 3× fundamental frequency increases gradually, while the amplitude of 2× fundamental frequency is shown as the opposite trend. For the viscosity of the lubricating oil μs, when its value increases, the amplitude of fundamental frequency and 3× fundamental frequency decrease gradually, while the trend of the amplitude of 2× fundamental frequency showed decreases at first, and increases at the end.
The SFD optimization parameter of rotor system is carried out by using NSGA-II, and the parameter configuration of NSGA-II algorithm is same as the parameter configuration with misalignment faults. The optimal solution of SFD is obtained as CS = 8 × 10−5; Ka = 6 × 105; μs = 5 × 10−2.
According to the design parameters, the numerical simulation of the SFD-rotor system is carried out under the state of the misalignment-rubbing coupling fault, and the dynamic response spectrum of the system is obtained, compared with spectrum in Figure 14a and Figure 17, it can be seen that after being optimized, the amplitude of fundamental frequency of system dynamics response decreases from 1.116 to 1.033, the amplitude decreases by about 7.4%; the amplitude of 2× fundamental frequency decreases from 0.099 to 0.048, the amplitude reduces by about 51.5%, the amplitude of 3× fundamental frequency decreases from 0.214 to 0.178, and the amplitude reduces by about 16.8%.

3.5. Comparison Results of SFD Parameters of the Rotor System under Three Conditions

Table 3 is SFD optimization parameters with normal, misalignment faults, and misalignment-rubbing coupling faults. It can be found that SFD optimization parameters are different under various conditions. The optimal value of the oil film gap is at the minimum level under the normal condition, and the optimal value is increased when the misalignment fault is considered. While the optimal value is relatively reduced, the misalignment-rubbing coupling faults are considered. For the effect of elastic support stiffness, the optimal value is the largest under the normal state, and the optimal value is the same under the state of the misalignment fault and the misalignment-rubbing coupling fault, while the optimal value reduces under the normal state. For the viscosity of the oil lubrication, the optimal value is the largest under the state of misalignment-rubbing coupling fault, while the optimal value reduces under the state of misalignment fault, and the optimal value relatively increases under the normal state. The reason that SFD optimization parameters are different under various conditions may be that various conditions have different influences on the rotor system, resulting in the change of dynamic characteristics of the rotor system.

4. Conclusions

In this paper, the vibration attenuation effects of the oil film gap (CS), elastic support stiffness (Ka) and the lubricating oil viscosity (μs) of SFD selected by using the multi-objective optimization method on the dynamic responses of the faulted system were analyzed. Some conclusions drawn from the study can be summarized as follows:
(1)
For the vibration responses caused by the system with misalignment fault, three variable parameters are strongly nonlinear, in which Ka has obvious impact on fundamental, 2× fundamental and 3× fundamental frequency, CS and μs are weak., For the vibration responses caused by a misalignment-rubbing coupling faults system, in which CS has a remarkable effect on fundamental frequency, μs and Ka are sensitive to 2× fundamental and 3× fundamental frequency, respectively.
(2)
Optimizing the parameters of SFD have an effect on vibration control and stability of the rotor system with varied faults by using NLPQL and NSGA-II algorithm. The vibration responses of the system with the optimization parameters are shown as: With no fault, the amplitude of fundamental frequency decreased by 23%. With the misalignment fault, the amplitude of the fundamental frequency decreased by 43.4%, the amplitude of 2× fundamental frequency decreased by 27.5%, while the amplitude of 3× fundamental frequency decreased by 66.7%. With the misalignment-rubbing coupling fault, the amplitude of fundamental frequency decreased by 7.4%, the amplitude of 2× fundamental frequency decreased by 51.5%, and the amplitude of 3× fundamental frequency decreased by 16.8%. Overall, the feasibility of the optimization method of the variable-structured SFD operational parameters for the faulted rotor system was verified.
(3)
The research results presented in the paper can provide helpful guidance for the design of the SFD and will be validated by experiments in the next planned step.

Author Contributions

L.M. developed the model and conducted the analyses; J.W. wrote and revised the paper; G.Z. contributed analysis tools.

Funding

This work is supported by the National Natural Science Foundation of China (No. 51705355).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

M = diag (Ms1, Ms1, Md1, Md1, Ms2, Ms2, Ms3, Ms3, Ms4, Ms4, Ms5, Ms5, Md2, Md2, Ms6, Ms6, Mo1, Mo1, Mo2, Mo2)
C = diag (Cz, Cz, Cp, Cp, Cz, Cz, 0, 0, 0, 0, Cz, Cz, Cp, Cp, Cz, Cz, 0, 0, 0, 0)
K = [ k 0 k 0 0 0 0 0 0 0 0 k 0 k 0 0 0 0 0 0 0 2 k 0 k 0 0 0 0 0 0 2 k 0 k 0 0 0 0 k 0 k 0 0 0 k 0 k 0 0 k a 0 0 0 k a 0 0 k a 0 s y m k a ]
Table A1. The model parameters of rotor shaft.
Table A1. The model parameters of rotor shaft.
ParametersValues
Ms1 = Ms2 = Ms3 = Ms4 = Ms5 = Ms6 (kg)1
Md1, Md2 (kg)15, 12
K (N/m)2.5 × 107
e1, e2 (m)1 × 10−5, 5 × 10−6
Cz, Cp (N·s/m)2100, 1500
Table A2. The model parameters of SFD.
Table A2. The model parameters of SFD.
ParametersValues
Rs, L (m)0.03, 0.011
Table A3. The model parameters of ball bearing.
Table A3. The model parameters of ball bearing.
ParametersValues
n7
Cb (N/m1.5)5 × 109
μ (m)3 × 10-6
Table A4. The model parameters of rubbing.
Table A4. The model parameters of rubbing.
ParametersValues
ks (N/m)5 × 107
δ (m)1 × 10−5
f0.1
Table A5. The model parameters of spline coupling.
Table A5. The model parameters of spline coupling.
ParametersValues
z14
Li, R (m)3.07 × 10−3, 6.95 × 10−2
T (N m)1000

References

  1. Olove, B. Non-linear behavior of a flexible shaft. SFD 1997, 206, 255–260. [Google Scholar]
  2. Humes, B. The Non-Linear Performance of Squeeze Film Damper; University of Sussex: Brighton, UK, 1977. [Google Scholar]
  3. Copper, S. Preliminary investigation of oil film for the control of vibration. Mech. Eng. 1963, 5, 305–315. [Google Scholar]
  4. Ruggiero, A.; D’Amato, R.; Magliano, E.; Kozak, D. Dynamical simulations of a flexible rotor in cylindrical uncavitated and cavitated lubricated journal bearings. Lubricants 2018, 6, 40. [Google Scholar] [CrossRef]
  5. Wang, J.; Zhang, J.; Ma, L.; Yu, Y. Study on the Effect Mechanism of MRD to Rotor System. Appl. Sci. 2019, 9, 2247. [Google Scholar] [CrossRef]
  6. Chang-Jian, C.W. Non-linear dynamic analysis of a HSFD mounted gear-bearing system. Nonlinear Dyn. 2010, 62, 333–347. [Google Scholar] [CrossRef]
  7. Chang-Jian, C.W.; Kuo, J.K. Bifurcation and chaos for porous squeeze film damper mounted rotor–bearing system lubricated with micropolar fluid. Nonlinear Dyn. 2009, 58, 697–714. [Google Scholar] [CrossRef]
  8. Chang-Jian, C.W.; Yau, H.T.; Chen, J.L. Nonlinear dynamic analysis of a hybrid squeeze-film damper-mounted rigid rotor lubricated with couple stress fluid and active control. Appl. Math. Model. 2010, 34, 2493–2507. [Google Scholar] [CrossRef]
  9. Inayat-Hussain, J.I. Bifurcations of a flexible rotor response in squeeze-film dampers without centering springs. Solitons Fractals 2005, 24, 583–596. [Google Scholar] [CrossRef]
  10. Inayat-Hussain, J.I. Bifurcations in the response of a flexible rotor in squeeze-film dampers with retainer springs. Solitons Fractals 2009, 39, 519–532. [Google Scholar] [CrossRef]
  11. Qin, W.; Zhang, J.; Wang, H. Response and bifurcation of rotor with squeeze film damper supported on elastic foundation. J. Northwestern Polytech. Univ. 2006, 24, 245–248. [Google Scholar]
  12. Yuan, H.; Wang, Z.; Wen, B. Analysis on stability of dual-discs rub-impact rotor with consideration of casing elasticity. J. Vib. Shock. 2010, 29, 52–54. [Google Scholar]
  13. Zhao, G.; Liu, Z.; Ye, J. Dynamic behavior of a rotor-misaligned spline coupling system. J. Vib. Shock 2009, 28, 78–82. [Google Scholar]
  14. Wen, B.; Gu, J.; Xia, S. The Advanced Rotor Dynamics; China Machine Press: Beijing, China, 2000. [Google Scholar]
  15. Cheng, M.; Meng, G. Dynamic analysis of a rotor-ball bearing nonlinear system with Alford force. China Mech. Eng. 2011, 22, 2806–2812. [Google Scholar]
  16. Park, J.S. Optimal latin-hypercube designs for computer experiments. J. Stat. Plan. Inference 1994, 39, 95–111. [Google Scholar] [CrossRef]
  17. Lei, D.; Yan, X. Multi-Objective Intelligent Optimization Algorithm and Its Application; Science Press: Beijing, China, 2009. [Google Scholar]
  18. Li, X. Dynamic Analysis and Identification for Rotor Misalignment-Rubbing Coupling Faults; Nanjing University of Aeronautics and Astronautics: Nanjing, China, 2009. [Google Scholar]
Figure 1. The simplified model of SFD-rotor-bearing system.
Figure 1. The simplified model of SFD-rotor-bearing system.
Applsci 09 03628 g001
Figure 2. Schematic diagram of the gear coupling.
Figure 2. Schematic diagram of the gear coupling.
Applsci 09 03628 g002
Figure 3. Latin Hypercube Design. (a) Random Latin hypercube and (b) optimal Latin hypercube.
Figure 3. Latin Hypercube Design. (a) Random Latin hypercube and (b) optimal Latin hypercube.
Applsci 09 03628 g003
Figure 4. The main process of NSGA-II.
Figure 4. The main process of NSGA-II.
Applsci 09 03628 g004
Figure 5. Schematic diagram of the variable-structured SFD.
Figure 5. Schematic diagram of the variable-structured SFD.
Applsci 09 03628 g005
Figure 6. Dynamic behavior of the normal rotor system. (a) Spectrogram; (b) Phase Diagram.
Figure 6. Dynamic behavior of the normal rotor system. (a) Spectrogram; (b) Phase Diagram.
Applsci 09 03628 g006
Figure 7. Comparison between calculated results and approximate model with no fault.
Figure 7. Comparison between calculated results and approximate model with no fault.
Applsci 09 03628 g007
Figure 8. Effects of SFD design parameters on the amplitude of the fundamental frequency of the normal system.
Figure 8. Effects of SFD design parameters on the amplitude of the fundamental frequency of the normal system.
Applsci 09 03628 g008
Figure 9. Spectrum diagram of optimized system with no fault.
Figure 9. Spectrum diagram of optimized system with no fault.
Applsci 09 03628 g009
Figure 10. Dynamic behavior of system with misalignment fault. (a) Spectrum diagram; (b) Phase Diagram.
Figure 10. Dynamic behavior of system with misalignment fault. (a) Spectrum diagram; (b) Phase Diagram.
Applsci 09 03628 g010
Figure 11. Comparison between calculated results and approximate model with misalignment fault. (a) Fundamental frequency; (b) 2× fundamental frequency; (c) 3× fundamental frequency.
Figure 11. Comparison between calculated results and approximate model with misalignment fault. (a) Fundamental frequency; (b) 2× fundamental frequency; (c) 3× fundamental frequency.
Applsci 09 03628 g011
Figure 12. Effect of SFD design parameters on the amplitude of the fundamental, 2× fundamental and 3× fundamental frequency of the system with misalignment fault.
Figure 12. Effect of SFD design parameters on the amplitude of the fundamental, 2× fundamental and 3× fundamental frequency of the system with misalignment fault.
Applsci 09 03628 g012
Figure 13. Spectrum diagram of SFD optimized system with misalignment fault.
Figure 13. Spectrum diagram of SFD optimized system with misalignment fault.
Applsci 09 03628 g013
Figure 14. Dynamic behavior of system with the misalignment-rubbing coupling fault. (a) Spectrum diagram; (b) Phase diagram.
Figure 14. Dynamic behavior of system with the misalignment-rubbing coupling fault. (a) Spectrum diagram; (b) Phase diagram.
Applsci 09 03628 g014
Figure 15. Comparison between calculated results and approximate model with the misalignment-rubbing coupling fault. (a) fundamental frequency; (b) 2× fundamental frequency; (c) 3× fundamental frequency.
Figure 15. Comparison between calculated results and approximate model with the misalignment-rubbing coupling fault. (a) fundamental frequency; (b) 2× fundamental frequency; (c) 3× fundamental frequency.
Applsci 09 03628 g015
Figure 16. Effect of SFD design parameters on the amplitude of the fundamental, 2× fundamental and 3× fundamental frequency of the system with the misalignment-rubbing coupling fault.
Figure 16. Effect of SFD design parameters on the amplitude of the fundamental, 2× fundamental and 3× fundamental frequency of the system with the misalignment-rubbing coupling fault.
Applsci 09 03628 g016
Figure 17. Spectrum diagram of SFD optimized system with the misalignment-rubbing coupling fault.
Figure 17. Spectrum diagram of SFD optimized system with the misalignment-rubbing coupling fault.
Applsci 09 03628 g017
Table 1. SFD Scheme structure variable.
Table 1. SFD Scheme structure variable.
Design ParametersInitial ValueThe Changing Range/mm
CS (mm)8 × 10−5(4 × 10−5, 2 × 10−4)
Ka (N·m)3 × 106(6 × 105, 8 × 107)
μs (Pa·s)0.02(0.005, 0.01)
Table 2. The design matrix of the SFD parameters.
Table 2. The design matrix of the SFD parameters.
Design Parameter12380
CS (mm)1.4 × 10−41.1 × 10−41.5 × 10−41.8 × 10−4
Ka (N·m)6 × 1068.7 × 1063.9 × 1064.8 × 106
μs (Pa·s)4.8 × 10−22.3 × 10−22.7 × 10−24.6 × 10−2
Table 3. SFD parameters of the rotor system under three conditions.
Table 3. SFD parameters of the rotor system under three conditions.
System StatusCS (mm)Ka (N·m)μs (Pa·s)
Before optimization8 × 10−53 × 1060.02
No fault
(after optimization)
4 × 10−58 × 1060.038
Misalignment faults
(after optimization)
2 × 10−46 × 1055 × 10−3
Misalignment-rubbing coupling faults
(after optimization)
8 × 10−56 × 1055 × 10−2

Share and Cite

MDPI and ACS Style

Ma, L.; Wang, J.; Zhang, G. The Dynamic Characteristic of a Faulted Rotor System with Multi-Objective Optimization Designed SFD. Appl. Sci. 2019, 9, 3628. https://doi.org/10.3390/app9173628

AMA Style

Ma L, Wang J, Zhang G. The Dynamic Characteristic of a Faulted Rotor System with Multi-Objective Optimization Designed SFD. Applied Sciences. 2019; 9(17):3628. https://doi.org/10.3390/app9173628

Chicago/Turabian Style

Ma, Liang, Jun Wang, and Guichang Zhang. 2019. "The Dynamic Characteristic of a Faulted Rotor System with Multi-Objective Optimization Designed SFD" Applied Sciences 9, no. 17: 3628. https://doi.org/10.3390/app9173628

APA Style

Ma, L., Wang, J., & Zhang, G. (2019). The Dynamic Characteristic of a Faulted Rotor System with Multi-Objective Optimization Designed SFD. Applied Sciences, 9(17), 3628. https://doi.org/10.3390/app9173628

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