Dynamics similarity design and verification of rotor system

Lijia Chen1 , Guihuo Luo2 , Jiaqi Han3 , Xi Kuang4 , Nan Zheng5 , Xuemin Liao6

1, 2, 3, 4, 5, 6College of Energy and Power Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing, 210016, China

2Corresponding author

Vibroengineering PROCEDIA, Vol. 28, 2019, p. 211-216. https://doi.org/10.21595/vp.2019.21048
Received 20 September 2019; accepted 28 September 2019; published 19 October 2019

Copyright © 2019 Lijia Chen, et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Abstract.

In order to study the dynamics similarity of original model similar to normal model, the similarity criteria and the similarity ratio of the normal and original models for the rotor system were derived by the dimension analysis method. ANSYS was used to numerically calculate the critical speeds, modal shapes and harmonic response of the original and normal models of rotor system. The analysis results show that, for the rotor system, the dynamic characteristics of the normal and original models satisfy the requirement of the similarity criteria perfectly. The dynamic characteristics of the original model can be predicted accurately by the corresponding normal model.

Dynamics similarity design and verification of rotor system

Keywords: normal model, dynamical similarity, dynamic characteristics.

1. Introduction

In terms of rotor dynamics test and dynamic similarity, many scholars have carried out a lot of useful explorations and achieved many results. Guoqiang H aimed at the simulation of the torsional vibration frequency of the rotor system, and used the physical structure similarity principle to design an analog rotor system similar to the physical structure of the 600 MW steam turbine shafting, The simulated rotor can better simulate the dynamic characteristics of the original steam turbine shaft [1]. Paul designed a rotor tester for rotor dynamics research based on similarity theory [2]. Guozhi Chen established a simulated rotor to simulate the dynamic characteristics of a gas generator rotor [3]. Wu had analyzed the dynamic similarity of the scale model rotor and the actual rotor system based on the similarity theory. Which showed that the frequency of the scale model is proportional to the actual rotor, and the scale model can be used to predict vibration characteristics to the actual rotor [4]. Zhongyong Gao applied the similarity theory to the design of the rotor tester. A simple rotor test model design method was proposed after the similarity of the model rotor and the actual rotor is discussed by the Ranlin’s formula and the matrix iteration formula [5]. Luo had put forward the similar condition of model rotor and actual rotor after researched equation analysis method and dimensional analysis method [6].

The similarity criteria of the normal model and the original model for the rotor system is obtained by the dimensional analysis method in this paper. The critical speeds of the normal model and the original model were calculated by the finite element software ANSYS, and the similarity is verified.

2. Dimensional method

According to the structural characteristics of the rotor system and its differential equation of motion, the relevant physical quantities of the rotor model include: excitation force F, excitation frequency ω0, geometry size l, time t, elastic modulus E, stress σ, strain s, Poisson ratio μ, damping ratio ξ, material density ρ, displacement response U, critical speeds f, bearing stiffness coefficient k and bearing damping coefficient c.

For the normal vibration operation of the rotor, the relationship of the physical quantities of the rotor system in the linear elastic range can be expressed as:

(1)
f ( F , l , t , E , s , e , x , m , r , w 0 , f , U , k , c ) = 0 .

According to the π- principle, π criteria can be derived as follows:

(2)
π 1 = F E l 2 ,         π 2 = t E l ρ ,           π 3 = ε ,           π 4 = ξ ,           π 5 = μ ,           π 6 = σ E , π 7 = ω 0 l ρ E ,           π 8 = f l ρ E ,           π 9 = U l ,           π 10 = k l E ,           π 11 = c l 2 E ρ .

In the similar model design of the rotor system, all physical quantities of the normal model and the original model must be proportional. The normal model and the original model are respectively represented by subscripts m and p. It is assumed that geometrical dimensions lm and lp satisfy the similarity ratio Cl, the times tm and tp satisfy the similarity ratio Ct, the densities ρm and ρp satisfy the similarity ratio Cρ, the elasticity modulus Em and Ep satisfy the similarity ratio CE, the excitation frequencies ω0m and ω0p satisfy the similarity ratio Cω0, the excitation forces Fm and Fp satisfy the similarity ratio CF, the displacement responses Um and Up satisfy the similarity ratio CU, the velocity responses Vm and Vp satisfy the similarity ratio Cv, the acceleration responses am and ap satisfy the similarity ratio Ca, and the support stiffness coefficients km and kp satisfy the similarity ratio Ck, the bearing damping coefficients cm and cp satisfy the similarity ratio Cc, the scale formula can be obtained.

According to the scaling formula, the dynamic similarity of the normal model and the original model can be written as [7]:

(3)
C t = t m t p = C l C p C E ,           C ω 0 = ω 0 m ω 0 p = C E C l C p ,           C F = F m F p = C E C l 2 ,
C U = U m U p = C l ,           C v = v m v p = C E C p ,           C a = a m a p = C v C t = C E C l C P ,
C f = f m f p = 1 C t = C E C l C ρ ,           C k = k m k p C E C l ,           C c = c m c p = C l 2 C E C p .

If the r-order mode of normal model and original model are ψrm and ψrp respectively, then ψirm/ψirp=Cl. It can be seen that when the model size is fixed, there is a direct proportion of ψrm and ψrp, and the modal shape is the complex amplitude ratio of each measurement point, substantially. If the modal shapes of two models are standardized on the elements in the same position, then:

(4)
ψ r m = ψ r p .

The modal shapes are exactly the same for the normal model and the original model in the same order.

In rotor model design, the normal model and original model are generally made of the same material, then CE=Cρ=1. If the geometric scale is Cl=λ, then the physical quantities of normal model and the original model are shown in Table 1.

Table 1. The scale of each physical quantity of normal model and original model of the rotor system

Parameter
Similarity ratio
Proportion
Feature length
C l
λ
Density
C ρ
1
Elastic modulus
C E
1
External force
C F
λ 2
Time
C T
λ
Excitation frequency
C ω 0
1/λ
Stiffness coefficient
C k
λ
Damping coefficient
C c
λ 2

3. Establishment of normal model of rotor system

As can be seen from Fig. 1. the original model consists of a rotating shaft and a disk, and two bearings with elasticity. Basic data of the original model is shown in Table 2. The finite element model of the rotor system is established according to the geometric model of the rotor system, as shown in Fig. 2. The rotating shaft is simulated by BEAM188, the disk is simulated by MASS21, and the bearing is simulated by COMBI214.

Fig. 1. Structure of rotor system of original model

 Structure of rotor system of original model

Fig. 2. Finite element model of the rotor system

 Finite element model of the rotor system

Table 2. Basic data of the original model

Parameter
Data
Quality of the disk / kg
1.401
Moment of inertia about diameter / (kg·m3)
0.002
Polar moment of inertia / (kg·m3)
0.00136
Bearing stiffness coefficients / (N/m)
4.348×107
Elasticity modulus / GPa
207.8
Poisson’s ratio
0.3
Density / (kg/m3)
7806

In order to analyze the dynamic similarity between the normal model and original model, according to the similarity criteria, the normal model is established. All the relevant parameters shown in Table 3 are strictly in accordance with the similarity relation for the similarity transformation.

Table 3. Basic data of the normal model of the rotor system

Basic parameters
Original model
Normal model
Geometric scale / λ
1
0.5
Material density / (kg/m3)
7806
7850
Modulus of elasticity / (GPa)
207.8
207.8
Stiffness coefficient /(N/m)
4.378×107
2.189×107

4. Similarity analysis of dynamic characteristics of normal model and original model

The finite element models of the normal model and the original model are established respectively, and the modal analysis is carried out by QR method. The critical speeds of the normal model and the original model can be obtained with simulation analysis that has been shown in Table 4.

Table 4 and 5 show that the critical speeds of normal model are twice as much as that of the original model and the strain energy distribution is exactly the same.

Table 4. The first three order critical speeds of normal model and original model

Frequency order
Original model / rpm
1/2 model / rpm
Error / %
1st – order
17167
34335
0.0029
2nd – order
48988
97977
0.0010
3rd – order
95653
191305
0.0005

The modal shapes of the normal model of the rotor system can be obtained by using ANSYS. At the same time, the modal displacements of each order can be extracted and normalized. The modal displacement diagram of each node is shown in Fig. 3.

Table 5. Strain energy distribution of the first 3rd order for the normal model and the original model

Frequency order
Original model
Normal model
Rotor
Bearing I
Bearing II
Rotor
Bearing I
Bearing II
1st – order
47.4
44.8
7.8
47.4
44.8
7.8
2nd – order
26.1
0.40
73.5
26.1
0.40
73.5
3rd – order
64.7
31.1
4.2
64.7
31.1
4.2

Fig. 3. Comparison and the first three critical modes of the original model and the normal model

 Comparison and the first three critical modes of the original model and the normal model

a) 1st -order critical mode before normalization

 Comparison and the first three critical modes of the original model and the normal model

b) 1st -order critical mode after normalization

 Comparison and the first three critical modes of the original model and the normal model

c) 2nd -order critical mode before normalization

 Comparison and the first three critical modes of the original model and the normal model

d) 2nd -order critical mode after normalization

 Comparison and the first three critical modes of the original model and the normal model

e) 3rd -order critical mode before normalization

 Comparison and the first three critical modes of the original model and the normal model

f) 3rd -order critical mode after normalization

With the similarity analysis of the normal model and the original model for the rotor system, it can be found that dynamic characteristics of the normal model for the rotor system is similar to the dynamic characteristics of the original model. Therefore, the modal shapes and the critical speeds of the original model for the rotor system can be accurately predicted by the normal model.

In order to compare and analyze the similarity of the dynamic response for the normal model and the original model, the harmonic response analysis is carried out. With the harmonic response analysis, the steady-state response of the complex structure under the harmonic load can be determined. The damping effect need to be considered when studying the forced vibration response characteristics of the rotor. If the original model damping is 5 N/(m/s), and then the normal model damping is 25 N/(m/s) according to the similarity ratio.

The modal superposition method is used to analysis rotor harmonic response. The exciting forces applied to the normal model and the original model are 2.5 N and 10 N, respectively. With harmonic response analysis, the steady state response of each rotor system at the disk is obtained. The frequency response curves are shown in Fig. 4 and Fig. 5.

Fig. 4. The frequency response curve of the normal model and the original model at the disk

The frequency response curve of the normal model and the original model at the disk

a) Response curve before normalization

The frequency response curve of the normal model and the original model at the disk

b) Response curve after normalization

Fig. 5. The frequency response curve of the normal model and the original model at the bearing

The frequency response curve of the normal model and the original model at the bearing

a) Response curve before normalization

The frequency response curve of the normal model and the original model at the bearing

b) Response curve after normalization

Fig. 4 and Fig. 5 show that the displacement-frequency curve of the normal model is basically similar to the curve of original model. It can be seen that the frequency dynamic response of two models are similar, so the normal model can accurately predict the frequency response characteristics of the original model subjected to external loads.

5. Conclusions

The normal similarity theory of the rotor system dynamic model was discussed in this paper. Based on the dimensional analysis method, the similarity criteria of the normal and original models for the rotor system were derived and the dynamic similarity relationship of two models were established. the dynamic characteristics and dynamic response similarity were analyzed by numerical simulation. the main conclusions show that:

1) The critical speeds of the normal and original models for the rotor system satisfy the determined relationship. in other words, if the geometrical size of the model is scaled down to 1/n times of the original model, the critical speeds of the normal model will increase to n times of the original model. In addition, the vibration modes and the dynamic characteristics of the two models are completely identical.

2) With the application of similarity criteria, the amplitudes of response curves of the normal model is consistent with the amplitudes of the original model. in addition, dynamic response of two models are similar. Therefore, the normal model can accurately predict the forced vibration response characteristics of the original model under external loads.

Therefore, the normal model can accurately predict the dynamic vibration characteristics of the original model rotor system such as modal shapes, critical speeds and harmonic response.

References

  1. Guoqiang H., Tao Y., Shuhong H. Research and development on integrated experimental platform of 600 MW turbine shafting. Research and Exploration in Laboratory, Vol. 29, Issue 2, 2010, p. 8-10. [CrossRef]
  2. Allaire Paul E., Li G. X., Lin Z. L., et al. Modeling of a high speeds rotor test rig with active magnetic bearings. Journal of Vibration and Acoustics, Vol. 128, 2006, p. 269-281. [Publisher]
  3. Guozhi C., Guofan W. Analysis and experimental study of the high speeds flexible rotor dynamic behaviors. Journal of Aerospace Power, Vol. 21, Issue 3, 2006, p. 563-568. [CrossRef]
  4. Wu J. J. Prediction of lateral vibration characteristics of a full-size rotor-bearing system by using those of its scale models. Finite Elements in Analysis and Design, Vol. 43, 2007, p. 803-816. [Publisher]
  5. Zhongyong G. Application of similarity theory in flexible rotor test. Machine Tool and Hydraulics, Vol. 5, 1987, p. 88-92. [CrossRef]
  6. Luo Z., Li J. Z., Yan Y. L., et al. Design method of dynamic similarity of the rotor-bearing system. Applied Mechanics and Materials, Vol. 300, 2013, p. 1148-1151. [Publisher]
  7. Baxi C. B., Telengator A., Razvi J. Rotor scale model tests for power conversion unit of GT-MHR. Nuclear Engineering and Design, 2011. [Publisher]