Dynamic Parameter Estimation of A Bend Rotor Bearing System

DOI : 10.17577/IJERTCONV6IS16015

Download Full-Text PDF Cite this Publication

Text Only Version

Dynamic Parameter Estimation of A Bend Rotor Bearing System

Dynamic Parameter Estimation of A Bend Rotor Bearing System

Sampath Kuma Kuppa

Industrial Design Department, National Institute of Technology Rourkela, Rourkela, 769008,

Odisha, India

Mohit Lal

Industrial Design Department, National Institute of Technology Rourkela, Rourkela, 769008,

Odisha, India

AbstractThe vibrational characteristics due to shaft bow effect in the rotor are different from the vibrational characteristics due to unbalance. From designing point of view it is important to know the cause of vibration is unbalance or rotor bow. Present study is mainly focused on the estimation of rotor bow parameters of a simple Jeffcott rotor model. In this model, a flexible shaft is mounted on the two bearings and a rigid disc is mounted on the shaft at midspan. Equation of motion is developed with the help of finite element method to determine the frequency response of the system. An identification algorithm is developed in the frequency domain for simultaneous estimation of the dynamic bearing parameters, residual unbalances and shaft bow. This identification algorithm is based on the least square technique. The present study pays attention to the simultaneous estimation of bearing dynamic parameters, residual unbalances and residual shaft bow parameters. The model is tested against different measurement noise and found to be robust.

KeywordsShaft Bow; Residual Unbalances; Finite Element Method (FEM); Least Squares Fit;

  1. INTRODUCTION

    Rotor dynamics is the branch of engineering that studies the lateral and torsional vibrations of rotating shafts for predicting the rotor vibrations and containing the vibration level under an acceptable limit. The principal components of a rotor-dynamic system are the shaft or rotor with disk, the bearings, coupling and seals. Flexible rotors are used in many industries is based on the design criteria. To maximizing the space available for other components such as impeller and seal, the shaft is relatively long shaft and thin geometry. Due to flexibility of rotor, bow effect arises in rotor. The bow effect is one kind rotor fault, which results in thermal defects and gravity sag. Vibration characteristics of a rotor bearing system having bowed shaft is different from vibration due to unbalance. For balancing point of view, it is important to know whether the cause of vibration is unbalance or shaft bow. The shaft bow in the large rotating machinery can lead to serious accident. Detail analysis of rotor response with bowed shaft has drawn attention of several researchers. In article [1] authors proposed a theoretical analysis of a single mass rotor with bowed shaft. The shaft is analyzed both theoretical and experimentally for balancing of the system. Reference [3] distinguished about the different types of rotor bows in rotating machines namely elastic bow, temporary elastic bow and permanent bow. Authors of [5] developed an

    the residual unbalance and bearing dynamic parameters. In article [7] authors developed a transient response analysis method for rotor transverse surface crack with a bow to estimate the effect of residual bow. During the study, the effect of residual bow is varying with the stiffness characteristics of rotating cracked shaft and dynamics of the cracked rotor. Reference [9] investigated the dynamic behavior of a geared-rotor system with viscoelastic supports under the effects of eccentricity, the transmission error of the gear mesh and the residual shaft bow is investigated under consideration of gear eccentricity, excitation of the gear transmission error and the residual shaft bow. Authors of [11- 15 & 17] proposed an algorithm to quantitatively estimate misalignment both numerically and experimentally. Ref. [16] used active magnetic bearing to suppress the excessive vibration generate due to misalignment. Authors of [18] developed algorithm to estimate speed dependent bearing and coupling misalignment parameters. The present study is focused on the estimation of rotor bow parameters. Developed an equation of motion incorporated with bow effect for of a simple rotor bearing system. A solution of this equation of motion gives the response of the given system. An identification algorithm prepared for estimating the parameters like bearing dynamic parameters, residual unbalances and rotor bow. The obtained parameters from the developed identification algorithm are comparing with assumed parameters.

  2. SYSTEM MODEL

    A simple rotor-bearing model consists of a rigid disc attached at the mid span and it is supported by two bearings at the end. The rotor shaft is considered as flexible with mass ms and diametral mass moment of inertia Is. The rotor-bearing model shown in Fig. 1 respectively, and in this study; the bearing is modeled as equivalent to spring and dampers system. The rotor is modeled as the Timoshenko beam element [2] for the brevity of this study; we are not paying attention to the gyroscopic effect. Here, B1 and B2 are the bearing. C is center of rotation the shaft, G is the center of gravity of the shaft and e is the eccentric of the rotation to the center of gravity of the shaft. Four displacements coordinates, two orthogonal translational and rotational coordinate is associated with each node (x, y, x, y). The disc has mass md and diametral mass moment of inertia is Id, the residual unbalance Ui. The bearings are modeled by eight linearized

    identification method for getting the response characteristics

    of flexible rotor model with shaft bow. Reference [6]

    c

    ij

    damping Bn

    ij

    and stiffness k Bn

    coefficients (all bearings has

    developed an identification algorithm for multi-degree-of- freedom rotor-bearing system for simultaneous estimation of

    different direct as well as cross-coupled terms, i.e. i, j =x, y

    and n=1,2 for two bearings). Since the main aim of this

    present study is to analyzed the effect of shaft bow on the rotor system and under the general shaft bow condition (i.e. shaft bow rises from

    1. Shaft Modelling

      The equations of motion of the shaft can be given as (12 degrees of freedom).

      gravity sagging) the actual shaft centerline and the theoretical

      M

      C

      K

      f t

      (1)

      shaft centerline will not be coincide which otherwise happens in shaft bow conditions. Hence, due to the self-weight of the shaft takes place which leads to bending in the shaft to the bending plane carrying a rotor bow effect on the rotor and the bearing.

      s s s s s s s

      Where,{}s is the nodal displacement vector, {f(t)}s is the force vector and mass matrix[M]s, stiffness matrix [K]s and damping matrix [C]s of the shaft are given in [10].

    2. Rigid Disc Modelling

    Fig. 1. Simple rotor bearing model

    unbalance URes =me. When rotor rotate in counter clockwise

    The disc for analysis is assumed axisymmetric and rigid with rotational and translational motion. In the mass matrix of the disc, the mass of the disc and mass moment of inertia of

    direction the relation between unbalance force of the rotor in y

    direction and x direction is given by the following formula that can be written as,

    the disc in terms are used to corresponding nodes. The equation of motion of the disc is discussed in [12], which is given below,

    E. Shaft Bow Force

    y x

    Funb jFunb

    (6)

    M

    t

    f t

    (2)

    The bow force of bent rotor bearing system i given by the

    d d d

    Where vectors {}d and {f(t)}d are the disc nodal displacement and force vectors respectively. Details of the

    formula, which is discussed in [1]. The bow force depends on

    the stiffness of the shaft, the amount of bow and phase angle between the residual unbalance.

    [M]d elemental disc matrix are given in [10].

    jt

    Fb Kr B e

    (7)

    C. Bearing Modelling

    Modelling of the bearing is given by the eight linearized damping and stiffness coefficient (all bearings have distant direct as well as cross-coupled coefficients). Equation of motion for bearing is discussed in [12], which is given below as,

    Where Kr is the rotor stiffness matrix, B is the rotor bow vector, is the location of the bow with respect to eccentricity, is the angular speed of the shaft.

    F. Assembled Equation of Motion of the System

    Combining the Eq. (1)-(7), the assembled equation of motion of the rotor bearing system is represented as

    CB

    KB

    f t

    (3)

    M C K F t

    F t

    (8)

    B B B

    Where {}B and {f}B are the bearing nodal displacement and force vectors. Details of elemental bearing are given in [10].

    D. Development of Unbalance Force

    Unbalance force in the rotor arises due to eccentricity between the center of gravity of the disc and center of rotation of the shaft. When the shaft rotates then due to the eccentricity a centrifugal forces is generated which is the cause of unbalance force, which is given by the following formula.

    unb bow

    Where M, C and K are the global mass, damping and stiffness matrices respectively. {F(t)}unb, {F(t)}bow and {}denote the vectors of generalized forces and displacements. The global rotor mass matrix M, rotor stiffness Kr and rotor displacement Cr. The damping C and stiffness K matrices combine the effects of the rotor and the bearing supports, i.e. K and C matrices for the system can be represented as,

    K Kr KB and C Cr CB

    jt

    funb t

    Funb e

    (4)

    where KB and CB are the bearing stiffness and damping

    Funb

    me2 ej

    (5)

    matrices.

    unb 2

    jt

    Where {Funb} is the residual unbalance force vector, which has

    F t

    URes e

    (9)

    phase and amplitude information an in generally the element

    Bow

    j t

    of the unbalance force vector are complex, mass of disc m and

    F t Kr B e

    (10)

    eccentricity of the rotor from the center of rotation of the shaft to the geometric center of the shaft e and the residual

    Where {URes} and {B} are the residual unbalance and residual shaft bow vector. These vectors contain the magnitude and

    phase information of residual unbalance and residual shaft bow, i.e. the element of these vectors is complex.

  3. METHODOLOGY

    This study is based on the parametric estimation of rotor bow of rotor in rotor bearing system. This has followed the steps

    Step 1: A system model is developing which consist a flexible rotor and this rotor is supported by two bearings at the end. A rigid disc is mounted on the midpoint of the shaft. Degree of freedom at each node is four and total degree of freedom of this system is 12.

    Step 2: An equation of motion has been developed for this system. Solution of this equation of motion gives information of the displacements at each node of the shaft. Response of the system is given the solution of the equation of motion.

    Step 3: For balancing of this system an identification algorithm has been developed which is based on the least square fit techniques. The regression equation is use to prepare the identification algorithm. A regression matrix is prepared which contained the responses and a regression vector, which contained the forces of the equation of motion.

    Step 4: After the estimation of the parameters. These parameters are used for developing the response and a draw a plot of response with assumed parameters and estimated parameters.

    Step 5: In this study, parameters estimated, when rotor rotates in counter clockwise direction.

    TABLE 1, represents the physical property of the bent rotor bearing system considered. These properties are used to develop the response, and further used in identification

    algorithm. The equation of motion for the bend rotor bearing system is taken from [1], after the getting equation of motion we have to calculate the solution of this equation, the solution is comes in the vector from because instead of all the constant here used only the matrices. The solution has the contained the complex numbers, the modulus of the complex numbers and direction of the complex numbers plotted with angular speed of the shaft is shown in Fig. 2. In Fig. 2, first two represents the horizontal and vertical displacement variation with respect to angular speed. The amplitude of the horizontal and vertical displacement are consider in these plots and last two represents the angular displacement in horizontal and vertical direction with respect to angular speed variation.

    Parameters

    Values (Units)

    Mass of the disc

    md = 5 kg

    Diametral

    =4/64 )

    mass

    moment

    of

    inertia

    (

    Id = 0.014 kg-m2

    Length of the shaft

    L =2.5m

    Diameter of the shaft

    d = 0.01m

    Modulus of Elasticity of shaft material (i.e. Steel)

    E = 210 GPa

    Modulus of rigidity of shaft material (i.e. Steel)

    G = 80.8 GPa

    Density of the shaft material

    = 7850 kg/m3

    The Residual Unbalance URes= md×e ,

    UXRes= 0.006 kg-

    (Uyres = -JUxres)

    m.

    where e is the eccentricity of the

    unbalance

    TABLE 1. PROPERTIES OF THE FLEXIBLE ROTOR SYSTEM

    Fig. 2. Semi-log representation of horizontal and vertical displacement with and without bow and Phase variation with the speed of the rotor

  4. DEVELOPMENT OF IDENTIFICATION ALGORITHM

    The present section deals with the response of the system to develop the identification algorithm. The generalized form of equation of motion of the system is discussed in [1], which is given below

  5. RESULTS AND DISCUSSION

    When rotor is runs in counter clockwise direction, the estimated bearing stiffness, bearing damping, shaft bow and residual unbalance parameters are shown in

    TABLE 2. From

    TABLE 2, the values of all the parameters without noise

    M C K F t

    F t

    (11)

    and with 1%, 2% and 5% noises is estimated and compared

    unb bow

    Assume the initial solution of the equation of motion is

    {}={}ejwt where {} contained the magnitude and phase information. Substitute the values of {}, {F(t)}unb and

    {F(t)}bow in the Eq.(11).

    2 M j C K

    (12)

    with the assumed values. From Fig. 3-Fig. 5, shows the comparison of all the parameters in with 1%, and 5% with the assumed parameter values. From Fig. 5, the shaft bow parameters w.r.t, 5% noise the maximum error is found in imaginary part of the residual shaft bow in x direction is58% error and rest of the parameters are below this limit. The

    2 U

    Re s r

    e j K Be j

    estimation of error in assumed and estimated parameters by

    The development of identification algorithm can be arranged suh that the all known quantity (i.e. rotor model parameters and measurable response) on the right hand side and all of the unknown quantities (i.e. bearing dynamic parameters, residual unbalance and residual shaft bow parameters) are placed on the left hand side. After rearranging the regression equation can be written as,

    Ai X B i (13)

    adding of the some percentage of noise in response ,

    therefore, the variation of assumed parameters to the estimated parameters with 1% and 5% of noise adding in the response. The variations of responses in all estimated parameters are represented in Fig. 6, which shows the variation of error in estimated parameters to the assumed parameters, and plot between percentage errors corresponding to the parameters.

    TABLE 2. ESTIMATED PARAMETERS WITHOUT NOISE AND WITH 1%, AND 5% NOISE ADDED IN RESPONSE

    Parameters

    Assumed values

    Estimated parameters

    Without Noise

    With 1%noise

    With 5% noise

    Kxx1 (N/m)

    2.50×105

    2.46×105

    2.48×105

    2.56×105

    Kxy1 (N/m)

    1.20×105

    1.22×105

    1.23×105

    1.27×105

    Kyx1 (N/m)

    1.35×105

    1.31×105

    1.29×105

    1.21×105

    Kyy1 (N/m)

    2.75×105

    2.77×105

    2.76×105

    2.72×105

    Kxx2 (N/m)

    2.75×105

    2.71×105

    2.69×105

    2.61×105

    Kxy2 (N/m)

    1.46×105

    1.48×105

    1.47×105

    1.43×105

    Kyx2 (N/m)

    1.39×105

    1.37×105

    1.37×105

    1.42×105

    Kyy2 (N/m)

    2.82×105

    2.84×105

    2.85×105

    2.87×105

    Cxx1 (Ns/m)

    300

    304.79

    303.40

    299.06

    Cxy1 (Ns/m)

    20

    23.75

    17.72

    11.95

    Cyx1 (Ns/m)

    50

    54.79

    51.40

    38.92

    Cyy1 (Ns/m)

    399

    402.75

    399.41

    386.46

    Cxx2 (Ns/m)

    315

    319.99

    317.53

    307.39

    Cxy2 (Ns/m)

    59

    62.79

    59.38

    48.98

    Cyx2 (Ns/m)

    79

    83.51

    82.84

    77.11

    Cyy2 (Ns/m)

    300

    303.19

    297.42

    274.67

    r (kg-m)

    0.0052

    0.0054

    0.0057

    0.0075

    i (kg-m)

    0.0030

    0.0031

    0.0035

    0.0046

    Bx (mm)

    0.0050

    0.0045

    0.00473

    0.002188

    By (mm)

    0.0050

    0.0055

    0.00295

    0.001254

    Fig. 3. Bearing stiffness parameters with 1% and 5% noise in simulated response

    Fig. 4. Bearing damping parameters with 1% and 5% noise in simulated response

    Fig. 5. Rotor bow and unbalance parameters with 1% and 5% noise in simulated response

    80

    0% Noise B

    y

    U

    i B

    Error Percentage

    60 1% Noise x x

    c

    1

    x

    b 5% Noise Ur xy

    40

    c

    b

    1

    yx cb2

    xy

    20

    0

    0 2 4 6 8 10 12 14 16 18 20

    Estimated Parameters

    Fig. 6. Variation of error in estimated parameters

  6. CONCLUSIONS

The main objective of the study is to estimate the parameters, which are responsible to deteriorate the efficiency of rotating machinery. There are various parameters in the rotating machines like bearing dynamic parameters, residual unbalances, shaft misalignment and residual bow, which are responsible for decreasing the efficiency. The present study is focused on the estimation of the bow effect of the rotor. Therefore, this study reveals the information related to parameters, which are obtained from the response with noise error. After that the concluding word is that the estimated parameters are from the presented algorithm is best because it helps to avoiding of the critical speeds consequently the main concern is fulfilled by this presented algorithm. The conclusion of this study is an estimation of parameters with help of presented identification algorithm give most suitable results. This algorithm gives the information of the appropriate parameters, which is using for the designing of the rotor bearing system.

REFERENCES

[1] J.C. Nicholas, E.J. Gunter, and P.E. Allaire, Effect of residual shaft bow on unbalance response and balancing of a single mass flexible rotorPart I: Unbalance response Journal of Engineering for Power, vol. 98(2), pp. 171-181, 1976.

[2] H.D. Nelson, A finite rotating shaft element using Timoshenko beam theory Journal of mechanical design, vol. 102(4), pp. 793-803, 1980.

[3] R.W. Clough, and J. Penzien, Dynamics of Structures" McGraw Hill, Inc., New York, USA, 1993.

[4] E.B. Donald, C.T. Hatch, and B. Grissom, Fundamentals of Rotating Machinery Diagnostics Bently Pressurized Bearing Company, 2002.

[5] S. Jana, and V. Arunkumar, Studies on the Effect of Shaft Bow for a Multi Mass Flexible Rotor Model in National Conference on Industrial Problems on Machines and Mechanisms, Kharagpur, 2005.

[6] R. Tiwari, Conditioning of regression matrices for simultaneous estimation of the residual unbalance and bearing dynamic parameters Mechanical Systems and Signal Processing, vol. 19(5), pp. 10821095, 2005.

[7] A. Darpe, K. Gupta, and A. Chawla, Dynamics of a bowed rotor with a transverse surface crack Journal of Sound and Vibration, vol. 296, pp. 888-907, 2006.

[8] G. Song, Z. Yang, C. Ji, and F. Wang, Theoreticalexperimental study on a rotor with a residual shaft bow Mechanism and Machine Theory, vol. 63, pp. 50-58, 2013.

[9] C. Kang, W. Hsu, E. Lee, and T. Shiau, Dynamic analysis of gear- rotor system with viscoelastic supports under residual shaft bow effect Mechanism and Machine Theory, vol. 46, pp. 264275, 2011.

[10] J.S. Rao, History of rotating machinery dynamics Springer Science & Business Media, vol. 20, 2011.

[11] M. Lal and R. Tiwari, Multifault identification in simple rotor bearingcoupling systems based on forced response measurements Mechanism and Machine Theory, vol. 51, pp. 87109, 2012.

[12] M. Lal and R. Tiwari, Identification of multiple faults with Incomplete response measurements in rotor-bearing-coupling systems In ASME Gas Turbine India Conference, pp. 613620, 2012.

[13] M. Lal, and R. Tiwari, Quantification of multiple fault parameters in flexible turbo-generator systems with incomplete rundown vibration data, Mechanical Systems and Signal Processing, vol. 41(1-2), pp.546-563, 2013.

[14] M. Lal and R. Tiwari, Identification of multiple fault parameters in a rigid-rotor and flexible-bearing-coupling system: An experimental investigation In ASME Gas Turbine India Conference, pp. 18, 2013.

[15] M. Lal and R. Tiwari, Experimental estimation of misalignment effects in rotor-bearing-coupling systems In Proceedings of the 9th IFToMM International Conference on Rotor Dynamics, Springer, pp. 779-789, 2015.

[16] S.K. Kuppa and M. Lal, Characteristic parameter estimation of AMB supported coupled rotor system In ASME Gas Turbine India Conference, pp. 1-9, 2017

[17] M. Lal and R. Tiwari, Experimental identification of shaft misalignment in a turbo-generator system Sadhana, vol. 43:80, pp. 1- 20, 2018.

[18] M. Lal and M. Satapathy, Estimation of speed dependent fault parametersin a coupled rotor-bearing system Archive of Mechanical Engineering, vol. 65 (3) pp. 327-247, 2018.

Leave a Reply