Rotordynamics

Last updated

Rotordynamics (or rotor dynamics) is a specialized branch of applied mechanics concerned with the behavior and diagnosis of rotating structures. It is commonly used to analyze the behavior of structures ranging from jet engines and steam turbines to auto engines and computer disk storage. At its most basic level, rotor dynamics is concerned with one or more mechanical structures (rotors) supported by bearings and influenced by internal phenomena that rotate around a single axis. The supporting structure is called a stator. As the speed of rotation increases the amplitude of vibration often passes through a maximum that is called a critical speed. This amplitude is commonly excited by imbalance of the rotating structure; everyday examples include engine balance and tire balance. If the amplitude of vibration at these critical speeds is excessive, then catastrophic failure occurs. In addition to this, turbomachinery often develop instabilities which are related to the internal makeup of turbomachinery, and which must be corrected. This is the chief concern of engineers who design large rotors.

Contents

Rotating machinery produces vibrations depending upon the structure of the mechanism involved in the process. Any faults in the machine can increase or excite the vibration signatures. Vibration behavior of the machine due to imbalance is one of the main aspects of rotating machinery which must be studied in detail and considered while designing. All objects including rotating machinery exhibit natural frequency depending on the structure of the object. The critical speed of a rotating machine occurs when the rotational speed matches its natural frequency. The lowest speed at which the natural frequency is first encountered is called the first critical speed, but as the speed increases, additional critical speeds are seen which are the multiples of the natural frequency. Hence, minimizing rotational unbalance and unnecessary external forces are very important to reducing the overall forces which initiate resonance. When the vibration is in resonance, it creates a destructive energy which should be the main concern when designing a rotating machine. The objective here should be to avoid operations that are close to the critical and pass safely through them when in acceleration or deceleration. If this aspect is ignored it might result in loss of the equipment, excessive wear and tear on the machinery, catastrophic breakage beyond repair or even human injury and loss of lives.

The real dynamics of the machine is difficult to model theoretically. The calculations are based on simplified models which resemble various structural components (lumped parameters models), equations obtained from solving models numerically (Rayleigh–Ritz method) and finally from the finite element method (FEM), which is another approach for modelling and analysis of the machine for natural frequencies. There are also some analytical methods, such as the distributed transfer function method, [1] which can generate analytical and closed-form natural frequencies, critical speeds and unbalanced mass response. On any machine prototype it is tested to confirm the precise frequencies of resonance and then redesigned to assure that resonance does not occur.

Basic principles

The equation of motion, in generalized matrix form, for an axially symmetric rotor rotating at a constant spin speed Ω is

where:

M is the symmetric Mass matrix
C is the symmetric damping matrix
G is the skew-symmetric gyroscopic matrix
K is the symmetric bearing or seal stiffness matrix
N is the gyroscopic matrix of deflection for inclusion of e.g., centrifugal elements.

in which q is the generalized coordinates of the rotor in inertial coordinates and f is a forcing function, usually including the unbalance.

The gyroscopic matrix G is proportional to spin speed Ω. The general solution to the above equation involves complex eigenvectors which are spin speed dependent. Engineering specialists in this field rely on the Campbell Diagram to explore these solutions.

An interesting feature of the rotordynamic system of equations are the off-diagonal terms of stiffness, damping, and mass. These terms are called cross-coupled stiffness, cross-coupled damping, and cross-coupled mass. When there is a positive cross-coupled stiffness, a deflection will cause a reaction force opposite the direction of deflection to react the load, and also a reaction force in the direction of positive whirl. If this force is large enough compared with the available direct damping and stiffness, the rotor will be unstable. When a rotor is unstable, it will typically require immediate shutdown of the machine to avoid catastrophic failure.

Jeffcott rotor

The Jeffcott rotor (named after Henry Homan Jeffcott), also known as the de Laval rotor in Europe, is a simplified lumped parameter model used to solve these equations. A Jeffcott rotor consists of a flexible, massless, uniform shaft mounted on two flexible bearings equidistant from a massive disk rigidly attached to the shaft. The simplest form of the rotor constrains the disk to a plane orthogonal to the axis of rotation. This limits the rotor's response to lateral vibration only. If the disk is perfectly balanced (i.e., its geometric center and center of mass are coincident), then the rotor is analogous to a single-degree-of-freedom undamped oscillator under free vibration. If there is some radial distance between the geometric center and center of mass, then the rotor is unbalanced, which produced a force proportional to the disk's mass, m, the distance between the two centers (eccentricity, ε) and the disk's spin speed, Ω. After calculating the equivalent stiffness, k, of the system, we can create the following second-order linear ordinary differential equation that describes the radial deflection of the disk from the rotor centerline. [2]

If we were to graph the radial response, we would see a sine wave with angular frequency . This lateral oscillation is called 'whirl', and in this case, is highly dependent upon spin speed. Not only does the spin speed influence the amplitude of the forcing function, it can also produce dynamic amplification near the system's natural frequency.

While the Jeffcott rotor is a useful tool for introducing rotordynamic concepts, it is important to note that it is a mathematical idealization that only loosely approximates the behavior of real-world rotors. [2]

Campbell diagram

Campbell Diagram for a Simple Rotor Campbelldiagram.png
Campbell Diagram for a Simple Rotor

The Campbell diagram, also known as "Whirl Speed Map" or a "Frequency Interference Diagram", of a simple rotor system is shown on the right. The pink and blue curves show the backward whirl (BW) and forward whirl (FW) modes, respectively, which diverge as the spin speed increases. When the BW frequency or the FW frequency equal the spin speed Ω, indicated by the intersections A and B with the synchronous spin speed line, the response of the rotor may show a peak. This is called a critical speed.

History

The history of rotordynamics is replete with the interplay of theory and practice. W. J. M. Rankine first performed an analysis of a spinning shaft in 1869, but his model was not adequate and he predicted that supercritical speeds could not be attained. In 1895, Dunkerley published an experimental paper describing supercritical speeds. Gustaf de Laval, a Swedish engineer, ran a steam turbine to supercritical speeds in 1889, and Kerr published a paper showing experimental evidence of a second critical speed in 1916.

Henry Jeffcott was commissioned by the Royal Society of London to resolve the conflict between theory and practice. He published a paper now considered classic in the Philosophical Magazine in 1919 in which he confirmed the existence of stable supercritical speeds. August Föppl published much the same conclusions in 1895, but history largely ignored his work.

Between the work of Jeffcott and the start of World War II there was much work in the area of instabilities and modeling techniques culminating in the work of Nils Otto Myklestad [3] and M. A. Prohl [4] which led to the transfer matrix method (TMM) for analyzing rotors. The most prevalent method used today for rotordynamics analysis is the finite element method.

Modern computer models have been commented on in a quote attributed to Dara Childs, "the quality of predictions from a computer code has more to do with the soundness of the basic model and the physical insight of the analyst. ... Superior algorithms or computer codes will not cure bad models or a lack of engineering judgment."

Prof. F. Nelson has written extensively on the history of rotordynamics and most of this section is based on his work.

Software

There are many software packages that are capable of solving the rotor dynamic system of equations. Rotor dynamic specific codes are more versatile for design purposes. These codes make it easy to add bearing coefficients, side loads, and many other items only a rotordynamicist would need. The non-rotor dynamic specific codes are full featured FEA solvers, and have many years of development in their solving techniques. The non-rotor dynamic specific codes can also be used to calibrate a code designed for rotor dynamics.

Rotordynamic-specific Software
NameCommercial / AcademicDescription
Dynamics R4 (Alfa-Tranzit Co. Ltd) [5] CommercialSoftware developed for design and analysis of spatial systems
AxSTREAM RotorDynamics, (SoftInWay) [6] Integrated Software platform for Rotor Dynamics, capable of lateral, torsional, and axial rotor dynamics for all widely used rotor types using the Finite Element Method on either beam or 2D-axisymmetric elements, and is capable of being automated.
Rotortest, (LAMAR - University of Campinas) [7] Finite Element Method based software, including different types of bearing solver. Developed by LAMAR (Laboratory of Rotating Machinery) - Unicamp (University of Campinas).
SAMCEF ROTOR [8] Software Platform for Rotors Simulation (LMS Samtech, A Siemens Business)
MADYN (Consulting engineers Klement) [9] CommercialCombined finite element lateral, torsional, axial and coupled solver for multiple rotors and gears, including foundation and housing.
MADYN 2000 (DELTA JS Inc.) [10] CommercialCombined finite element (3D Timoshenko beam) lateral, torsional, axial and coupled solver for multiple rotors and gears, foundations and casings (capability to import transfer functions and state space matrices from other sources), various bearings (fluid film, spring damper, magnetic, rolling element)
NASTRAN (MSC / Hexagon)CommercialFinite element code which can be used to compute Campbell diagrams through SOL 107.
iSTRDYN (DynaTech Software LLC) [11] Commercial2-D Axis-symmetric finite element solver
FEMRDYN (DynaTech Engineering, Inc.) [12] Commercial1-D Axis-symmetric finite element solver
DyRoBeS (Eigen Technologies, Inc.) [13] Commercial1-D beam element solver
RIMAP (RITEC) [14] Commercial1-D beam element solver
XLRotor (Rotating Machinery Analysis, Inc.) [15] Commercial1-D beam element solver, including magnetic bearing control systems and coupled lateral-torsional analysis. A powerful, fast and easy to use tool for rotor dynamic modeling and analysis using Excel spreadsheets. Readily automated with VBA macros, plus a plugin for 3D CAD software.
ARMD (Rotor Bearing Technology & Software, Inc.) [16] CommercialFEA-based software for rotordynamics, multi-branch torsional vibration, fluid-film bearings (hydrodynamic, hydrostatic, and hybrid) design, optimization, and performance evaluation, that is used worldwide by researchers, OEMs and end-users across all industries.
XLTRC2 (Texas A&M) [17] Academic1-D beam element solver
ComboRotor (University of Virginia) [18] Combined finite element lateral, torsional, axial solver for multiple rotors evaluating critical speeds, stability and unbalance response extensively verified by industrial use
MESWIR (Institute of Fluid-Flow Machinery, Polish Academy of Sciences) [19] AcademicComputer code package for analysis of rotor-bearing systems within the linear and non-linear range
RoDAP (D&M Technology) [20] CommercialLateral, torsional, axial and coupled solver for multiple rotors, gears and flexible disks(HDD)
ROTORINSA (ROTORINSA) [21] CommercialFinite element software developed by a French engineering school (INSA-Lyon) for analysis of steady-state dynamic behavior of rotors in bending.
COMSOL Multiphysics, Rotordynamics Module add-on (Rotordynamics Module) [22]
RAPPID - (Rotordynamics-Seal Research) [23] CommercialFinite element based software library (3D solid and beam elements) including rotordynamic coefficient solvers

See also

Related Research Articles

<span class="mw-page-title-main">Angular momentum</span> Conserved physical quantity; rotational analogue of linear momentum

In physics, angular momentum is the rotational analog of linear momentum. It is an important physical quantity because it is a conserved quantity – the total angular momentum of a closed system remains constant. Angular momentum has both a direction and a magnitude, and both are conserved. Bicycles and motorcycles, flying discs, rifled bullets, and gyroscopes owe their useful properties to conservation of angular momentum. Conservation of angular momentum is also why hurricanes form spirals and neutron stars have high rotational rates. In general, conservation limits the possible motion of a system, but it does not uniquely determine it.

<span class="mw-page-title-main">Electric motor</span> Machine that converts electrical energy into mechanical energy

An electric motor is an electrical machine that converts electrical energy into mechanical energy. Most electric motors operate through the interaction between the motor's magnetic field and electric current in a wire winding to generate force in the form of torque applied on the motor's shaft. An electric generator is mechanically identical to an electric motor, but operates in reverse, converting mechanical energy into electrical energy.

Hmethods are used in control theory to synthesize controllers to achieve stabilization with guaranteed performance. To use H methods, a control designer expresses the control problem as a mathematical optimization problem and then finds the controller that solves this optimization. H techniques have the advantage over classical control techniques in that H techniques are readily applicable to problems involving multivariate systems with cross-coupling between channels; disadvantages of H techniques include the level of mathematical understanding needed to apply them successfully and the need for a reasonably good model of the system to be controlled. It is important to keep in mind that the resulting controller is only optimal with respect to the prescribed cost function and does not necessarily represent the best controller in terms of the usual performance measures used to evaluate controllers such as settling time, energy expended, etc. Also, non-linear constraints such as saturation are generally not well-handled. These methods were introduced into control theory in the late 1970s-early 1980s by George Zames, J. William Helton , and Allen Tannenbaum.

<span class="mw-page-title-main">Centrifugal compressor</span> Sub-class of dynamic axisymmetric work-absorbing turbomachinery

Centrifugal compressors, sometimes called impeller compressors or radial compressors, are a sub-class of dynamic axisymmetric work-absorbing turbomachinery.

<span class="mw-page-title-main">Rigid body dynamics</span> Study of the effects of forces on undeformable bodies

In the physical science of dynamics, rigid-body dynamics studies the movement of systems of interconnected bodies under the action of external forces. The assumption that the bodies are rigid simplifies analysis, by reducing the parameters that describe the configuration of the system to the translation and rotation of reference frames attached to each body. This excludes bodies that display fluid, highly elastic, and plastic behavior.

<span class="mw-page-title-main">Euler's Disk</span> Scientific educational toy

Euler's Disk, invented between 1987 and 1990 by Joseph Bendik, is a trademarked scientific educational toy. It is used to illustrate and study the dynamic system of a spinning and rolling disk on a flat or curved surface. It has been the subject of several scientific papers.

Dynamical simulation, in computational physics, is the simulation of systems of objects that are free to move, usually in three dimensions according to Newton's laws of dynamics, or approximations thereof. Dynamical simulation is used in computer animation to assist animators to produce realistic motion, in industrial design, and in video games. Body movement is calculated using time integration methods.

In rotordynamics, the rigid rotor is a mechanical model of rotating systems. An arbitrary rigid rotor is a 3-dimensional rigid object, such as a top. To orient such an object in space requires three angles, known as Euler angles. A special rigid rotor is the linear rotor requiring only two angles to describe, for example of a diatomic molecule. More general molecules are 3-dimensional, such as water, ammonia, or methane.

In classical mechanics, the Newton–Euler equations describe the combined translational and rotational dynamics of a rigid body.

<span class="mw-page-title-main">Rotation around a fixed axis</span> Type of motion

Rotation around a fixed axis or axial rotation is a special case of rotational motion around an axis of rotation fixed, stationary, or static in three-dimensional space. This type of motion excludes the possibility of the instantaneous axis of rotation changing its orientation and cannot describe such phenomena as wobbling or precession. According to Euler's rotation theorem, simultaneous rotation along a number of stationary axes at the same time is impossible; if two rotations are forced at the same time, a new axis of rotation will result.

In solid mechanics, in the field of rotordynamics, the critical speed is the theoretical angular velocity that excites the natural frequency of a rotating object, such as a shaft, propeller, leadscrew, or gear. As the speed of rotation approaches the object's natural frequency, the object begins to resonate, which dramatically increases system vibration. The resulting resonance occurs regardless of orientation. When the rotational speed is equal to the natural frequency, then that speed is referred to as a critical speed.

<span class="mw-page-title-main">Bently Nevada</span> Technology and service company

Bently Nevada is an asset protection and condition monitoring hardware, software and service company for industrial plant-wide operations. Its products are used to monitor the mechanical condition of rotating equipment in a wide variety of industries including oil and gas production, hydroelectric, wind, hydrocarbon processing, electric power generation, pulp and paper, mining, water and wastewater treatment. The company was founded in 1961 by Don Bently. Bently Nevada is headquartered in Minden, Nevada, about one hour south of Reno. Don Bently was the first to manufacture a commercially successful eddy-current proximity probe which measured vibration in high-speed turbomachinery by allowing the direct observation of the rotating shaft. The company also performed research in the field of rotordynamics, furthering knowledge of machinery malfunctions such as shaft cracks and fluid-induced instabilities. Its research also helped refine the equations used to describe vibratory behavior in rotordynamic systems.

<span class="mw-page-title-main">Vibration</span> Mechanical oscillations about an equilibrium point

Vibration is a mechanical phenomenon whereby oscillations occur about an equilibrium point. Vibration may be deterministic if the oscillations can be characterised precisely, or random if the oscillations can only be analysed statistically.

<span class="mw-page-title-main">Campbell diagram</span> Plot in engineering

A Campbell diagram plot represents a system's response spectrum as a function of its oscillation regime. It is named for Wilfred Campbell, who introduced the concept. It is also called an interference diagram.

In physics, nonlinear resonance is the occurrence of resonance in a nonlinear system. In nonlinear resonance the system behaviour – resonance frequencies and modes – depends on the amplitude of the oscillations, while for linear systems this is independent of amplitude. The mixing of modes in non-linear systems is termed resonant interaction.

Rotating unbalance is the uneven distribution of mass around an axis of rotation. A rotating mass, or rotor, is said to be out of balance when its center of mass is out of alignment with the center of rotation. Unbalance causes a moment which gives the rotor a wobbling movement characteristic of vibration of rotating structures.

Transformation of three phase electrical quantities to two phase quantities is a usual practice to simplify analysis of three phase electrical circuits. Polyphase a.c machines can be represented by an equivalent two phase model provided the rotating polyphases winding in rotor and the stationary polyphase windings in stator can be expressed in a fictitious two axes coils. The process of replacing one set of variables to another related set of variable is called winding transformation or simply transformation or linear transformation. The term linear transformation means that the transformation from old to new set of variable and vice versa is governed by linear equations. The equations relating old variables and new variables are called transformation equation and the following general form:

 [new Variable] = [transformation matrix][old variable]  [old Variable] = [transformation matrix][new variable]

Dynamic Substructuring (DS) is an engineering tool used to model and analyse the dynamics of mechanical systems by means of its components or substructures. Using the dynamic substructuring approach one is able to analyse the dynamic behaviour of substructures separately and to later on calculate the assembled dynamics using coupling procedures. Dynamic substructuring has several advantages over the analysis of the fully assembled system:

<span class="mw-page-title-main">Critical embankment velocity</span>

Critical embankment velocity or critical speed, in transportation engineering, is the velocity value of the upper moving vehicle that causes the severe vibration of the embankment and the nearby ground. This concept and the prediction method was put forward by scholars in civil engineering communities before 1980 and stressed and exhaustively studied by Krylov in 1994 based on the Green function method and predicted more accurately using other methods in the following. When the vehicles such as high-speed trains or airplanes move approaching or beyond this critical velocity, the vibration magnitudes of vehicles and nearby ground increase rapidly and possibly lead to the damage to the passengers and the neighboring residents. This relevant unexpected phenomenon is called the ground vibration boom from 1997 when it was observed in Sweden for the first time.

Agnieszka Muszyńska (1935–2024) was a Polish-American mechanical engineer specializing in rotordynamics and the vibrations of rotating machinery.

References

Notes

  1. Liu, Shibing; Yang, Bingen (2017-02-22). "Vibrations of Flexible Multistage Rotor Systems Supported by Water-Lubricated Rubber Bearings". Journal of Vibration and Acoustics. 139 (2): 021016–021016–12. doi:10.1115/1.4035136. ISSN   1048-9002.
  2. 1 2 Matsushita, Osami; Tanaka, Masato; Kanki, Hiroshi; Kobayashi, Masao; Keogh, Patrick (2017). Vibrations of Rotating Machinery (1st ed.). Japan: Springer. pp. 199–206. ISBN   978-4-431-55455-4.
  3. Myklestad, Nils (April 1944). "A New Method of Calculating Natural Modes of Uncoupled Bending Vibration of Airplane Wings and Other Types of Beams". Journal of the Aeronautical Sciences (Institute of the Aeronautical Sciences). 11 (2): 153–162. doi:10.2514/8.11116.
  4. Prohl, M. A. (1945), "A General Method for Calculating Critical Speeds of Flexible Rotors", Trans ASME, 66: A–142
  5. "ROTORDYNAMICS OF TURBOMACHINERY, Software and Engineering services, Alfa-Tranzit Co". www.alfatran.com. Retrieved 2022-01-11.
  6. "AxSTREAM - Rotordynamics | Turbomachinery Design Software". www.softinway.com. 2016-07-27. Retrieved 2022-01-11.
  7. "Home". www.fem.unicamp.br. Retrieved 2022-01-11.
  8. "Rotor dynamics analysis: SAMCEF Rotors is up to the challenge". Archived from the original on 2012-09-29. Retrieved 2013-06-04.
  9. "Home". www.madyn.de. Retrieved 2022-01-11.
  10. "MADYN 2000 for Rotordynamics". www.delta-js.ch. Retrieved 2022-01-11.
  11. "Archived copy". Archived from the original on 2021-05-07. Retrieved 2022-01-11.{{cite web}}: CS1 maint: archived copy as title (link)
  12. "FEMRDYN". dynatechengr.com. Retrieved 2022-01-11.
  13. "Dyrobes – A Revolution in Rotor Dynamics Software" . Retrieved 2022-01-11.
  14. "RITEC - RIMAP (Rotordynamic Analysis Software Package)". Archived from the original on 2013-10-04. Retrieved 2013-10-01.
  15. "Home". XLRotor. Retrieved 2022-01-11.
  16. "ARMD Rotor Dynamics" . Retrieved 2022-01-01.
  17. "Turbomachinery Laboratory". turbolab.tamu.edu. Retrieved 2022-01-11.
  18. "Rotating Machinery and Controls Laboratory".
  19. "the Szewalski Institute of Fluid-flow Machinery". Archived from the original on 2008-12-15.
  20. "D&M; Technology Co., Ltd". Archived from the original on 2021-05-06.
  21. "LaMCoS". lamcos.insa-lyon.fr. Retrieved 2022-01-11.
  22. "Rotordynamics Software for Predictive Analyses of Rotating Machines". www.comsol.com. Retrieved 2022-01-11.
  23. "Rotordynamics-Seal Research - RSR". rda.guru. Retrieved 2022-01-11.