ROTOR DYNAMIC ANALYSIS OF STEAM TURBINE ROTOR
Int. J. Mech. Eng. & Rob. Res. 2014Nagaraju Tenali and Srinivas Kadivendi, 2014ISSN 2278 ñ 0149 www.ijmerr.comVol. 3, No. 1, January 2014 2014 IJMERR. All Rights ReservedResearch PaperROTOR DYNAMIC ANALYSIS OF STEAMTURBINE ROTOR USING ANSYSNagaraju Tenali1* and Srinivas Kadivendi1*Corresponding Author: Nagaraju Tenali, tenali.n1830@gmail.comRotor dynamics is a field under mechanics. Mainly deals with the vibration of rotating structures.In recent days, the study about rotordynamics has gained more importance within steamturbine industries. The main reason is steam turbine consists of many rotating parts constitutesa complex dynamic system. While designing rotors of high speed turbo machineries, it is ofprime importance consider rotor dynamics characteristics in to account. And the world we areliving in today is pushing the technology harder and harder. The products need to get better andtoday they also need to be friendlier to the environment. To get better products we need betteranalysis tools to optimize them and to get closer to the limit what the material can withstand.The modeling features for rotor and bearing support flexibility are described in this thesis. Byintegrating these characteristic rotor dynamics features into the standard FEA- modal, harmonicand transient analysis procedures found in ANSYS we can analyze and determine the designintegrity of rotating equipment. Some ideas are presented to deal with critical speeds calculationusing ANSYS. This Thesis shows how elements BEAM188 and COMBI214 are used to modelthe shaft and bearings, respectively.The purpose of a standard rotor dynamics analysis of Steamturbine rotor is to enable an engineer to characterize the lateral dynamics design characteristicsof a given design. With the Campbell plots, we can determine critical speeds and system stability.These techniques, along with a same modeling and results are also calculated from TMS-050to verify ANSYS result with testing result for unbalance response.Keywords: Ansys, Critical speed, Rotor, Rotor dynamics, Steam turbine, TMS-050, VibrationsINTRODUCTIONmachines are accompanied by higherrequirements for their reliability. To increaseoperational life of turbo machines is also oneof the main tasks of quality improvement. Inthis connection at present, when developingand mastering the steam turbines, moderncomputational and experimental methods areused to determine strength and reliabilitySteam turbine plant is an integral part ofthermal power station. Thereforedevelopment, construction and improvementof steam turbine are an important field ofdevelopment of power industry. Growth inpower and more complicated design of turbo1Department of Mechanical Engineering, D V R & D H S MIC College of Technology, Kanchikacherla, Krishna District - 521 180,A.P., India338
Int. J. Mech. Eng. & Rob. Res. 2014Nagaraju Tenali and Srinivas Kadivendi, 2014is to enable an engineer to characterizethe lateral dynamics design characteristicsof a given design. While analysis of somerotating equipment may require analysisspecific to the unit, a general method hasemerged for performing the standard lateralanalysis.characteristics. Rotor dynamics is the branchof engineering that studies the lateral andtorsional vibrations of rotating shafts, with theobjective of predicting the rotor vibrations andcontaining the vibration level under anacceptable limit. The principal componentsof a rotor-dynamic system are the shaft orrotor with disk, the bearings, and the seals.The shaft or rotor is the rotating componentof the system. Basically there are three formsof vibrations associated with the motion ofthe rotor: torsional, axial and lateral. Torsionalvibration is the dynamics of the shaft in theangular/rotational direction. Normally, this islittle influenced by the bearings that supportthe rotor. Axial vibration is the dynamics ofthe rotor in the axial direction and is generallynot a major problem. Lateral vibration, theprimary concern, is the vibration of the rotorin lateral directions.Fundamental EquationThe general form of equation of motion forall vibration problems is given by,.(1.1)Where,[M] symmetric mass matrix[C] symmetric damping matrix[K] symmetric stiffness matrix[f] external force vector[u] generalized coordinate vector Inrotordynamics, this equation of motion canbe expressed in the following general form[3],The bearings play a huge part indetermining the lateral vibrations of the rotor.In this thesis, we will study the basic conceptsof the lateral rotor dynamics of turbomachinery. With ever increase in demand forlarger size and velocity in modern machines,Rotor Dynamics became more and more animportant subject in the mechanicalengineering design. It is well know thattorsional vibration in rotating machines,reciprocating machines installation andgeared system, whirling of rotating shaft, theeffect of flexible bearing, instabilities due toasymmetriccross-sectionshafts,hydrodynamics bearings, hysteresis,balancing of rigid and flexible rotor can beunderstood only on the basis of rotordynamics studies. Rotor dynamics is anextremely important branch of the disciplineof dynamics that pertains to the operation andbehavior of a huge assortment of rotatingmachines. The purpose of a standardrotor dynamics analysis and design audit.(1.2)The above mentioned equation (1.2)describes the motion of an axially symmetricrotor, which is rotating at constant spin speedabout its spin axis. This equation is justsimilar to the general dynamic equationexcept it is accompanied with skewsymmetric gyroscopic matrix, [C gyro] andskew-symmetric circulatory matrix [H]. Thegyroscopic and circulatory matrices [C gyro]and [H] are greatly influenced by rotationalvelocity . When the rotational velocity ,tends to zero, the skew-symmetric termspresent in equation (1.2) vanish andrepresent an ordinary stand still structure.The gyroscopic matrix [C gyro] containsinertial terms and that are derived from kineticenergy due to gyroscopic moments acting on339
Int. J. Mech. Eng. & Rob. Res. 2014Nagaraju Tenali and Srinivas Kadivendi, 2014the rotating parts of the machine. If thisequation is described in rotating referenceframe, this gyroscopic matrix [C gyro] alsocontains the terms associated with Cariolesacceleration. The circulatory matrix, [H] iscontributed mainly from internal damping ofrotating elements (XU Yang et al., 2004).When the rotor is rotating at constantrotational speed, the equation of motionfor the mass center can be derived fromNewton’s law of motion and it is expressedin the following form.(1.3)TheoryThe concept of rotor dynamics can be easilydemonstrated with the help of generalizedLaval-Jeffcott rotor modal as shown in Figure1.(1.4)The above equations can be re-written as,.(1.5)Figure 1: Generalized Laval- JeffcottRotor Model.(1.6)is the phase angle of the massWhere,unbalance. The above equations of motionsshow that the motions in X and Y directionare both dynamically and statically decoupledin this model. Therefore, they can be solvedseparately.Determination of naturalfrequenciesFor this simple rotor model, the undampednatural frequency, damping ration and thedamped natural frequency of the rotor modelfor X and Y direction can be calculated fromThe generalized Laval-Jeffcott rotor consistsof long, flexible mass less shaft with flexiblebearings on both the ends. The bearings havesupport stiffness of KX and KY associatedwith damping CX and CY in x and y directionrespectively. There is a massive disk of mass,m located at the center of shaft. The centerof gravity of disk is offset from the shaft geometric center by an eccentricity of e. Themotion of the disk center is described by twotranslational displacements (x, y) as shownin Figure 2.(1.7)Steady state response tounbalanceFor single unbalance force, as present in thiscase, the can be set to zero. Therefore theequations (1.5) and (1.6) becomes,Figure 2: End View of Laval- Jeffcott Rotor.(1.8).(1.9)340
Int. J. Mech. Eng. & Rob. Res. 2014Nagaraju Tenali and Srinivas Kadivendi, 2014Then the solution for the response is,Table 1: Sectional D)(mm)Temp 109.007560211.0071603116.506360418.0013660MODELING AND DESIGNDATA INPUT554.006360618.0011060733.5010060The modeling features for rotor and bearingsupport flexibility are described in this thesis,and shows how elements BEAM188,COMBI214 are used to model the shaft andbearings. And MASS21 used to model theadditional 4831.0010560491.0010560.(1.11)RotorBEAM188 Element Description: BEAM188is suitable for analyzing slender to moderatelyStubby/thick beam structures. BEAM188 isa linear (2-node) or a quadratic beam elementin 3-D. BEAM188 has six or seven degreesof freedom at each node, with the numberof degrees of freedom depending on thevalue of KEYOPT(1). When KEYOPT (1) 0(the default), six degrees of freedom occurat each node. These include translations inthe x, y, and z directions and rotations aboutthe x, y, and z directions. When KEYOPT (1) 1, a seventh degree of freedom (warpingmagnitude) is also considered. This elementis well-suited for linear, large rotation, and/or large strain nonlinear applications.Figure 3: Beam Geometry341
Int. J. Mech. Eng. & Rob. Res. 2014Nagaraju Tenali and Srinivas Kadivendi, 2014Where, L Length of the each section (mm)Two bearings used in this thesis, one at frontside (at section 08) and other at rear side (atsection 43). And “COMBI-214” Element usedfor Modeling of the Bearing in ANSYS.COMBI214 Element Description: 2-D SpringDamper Bearing. COMBI214 has longitudinalas well as cross-coupling capability in 2-Dapplications. It is a tension compressionelement with up to two degrees of freedomat each node: translations in any two nodaldirections (x, y, or z). COMBI214 has twonodes plus one optional orientation node. Nobending or torsion is considered. The springdamper element has no mass.D Diameter of each section (mm)Blades (Additional Disks Masses)MASS21 Element Description: MASS21 is apoint element having up to six degrees offreedom: translations in the nodal x, y, and zdirections and rotations about the nodal x, y,and z axes. A different mass and rotary inertiamay be assigned to each coordinatedirection.Figure 4: Mass 21 GeometryTable 3: Bearing DetailsLength Diameter Section(mm)(mm) LocationType ofbearingFront Bearing401008Tilting padRear Bearing4010043Tilting padTable 4: Properties for 1st BearingSpeed(RPM)KXX 103KXYKYXKYY 103CXX 103CXYCYXCYY 001594941020038.25342
Int. J. Mech. Eng. & Rob. Res. 2014Nagaraju Tenali and Srinivas Kadivendi, 2014Bearing DetailsTable 2: Desk Input DataDiskNo.Mass ofDisk 36.549Figure 5: Combi 214 GeometryTable 5: Properties for 2nd BearingSpeed(RPM)KXX 103KXYKYXKYY 103CXX 103CXYCYXCYY 51.0213000131057001544951020038.25Rotor Material PropertiesTable 6: Rotor with MeshingTable 6: Material PropertyYoung Modulus ‘E2.1 1011 N/m2Poisson Ratio (µ )0.25Density ‘ñ’7800 kg/m3343
Int. J. Mech. Eng. & Rob. Res. 2014Nagaraju Tenali and Srinivas Kadivendi, 2014Figure 7: Finite Element Model of RotorFigure 9: Mode Shape 2RESULTS AND DISCUSSIONThe different analysis carried out to the rotordynamic integrity of the steam turbine rotorunder given loads.Figure 10: Mode Shape 3Analysis TypesA. Modal analysisB. Harmonic analysisC. Transient analysisA. Modal AnalysisFigure 8, Figure 9, Figure 10 and Figure11shows fist four mode shape and dampednatural frequency at the operating speed(11800 rpm) by Ansys.Figure 11: Mode Shape 4Figure 8: Mode Shape 1344
Int. J. Mech. Eng. & Rob. Res. 2014Nagaraju Tenali and Srinivas Kadivendi, 2014Table no 8 Shows comparisons of the twodifferent tools for undamped naturalfrequency (Hz) at the operating speed (11800rpm).Table No.9 Show comparison of thedamped natural frequency (Hz) of twodifferent tools at the operating speed (11800rpm).Table No.8 Undamped natural frequency(Hz) at the operating speed.Critical speed and Campbell diagramanalysis.In this analysis, a number of Eigenfrequency analyses are performed on thesteam turbine rotor model for the speed rangestarting from 0 rpm to 12000 rpm with anincrement of 150 rpm using multiple loadsteps.Table 8: Undamped Natural Frequency(Hz) at the Operating SpeedMode 13591361Figure 13: Damped Critical Speed by AnsysFigure 12. Shows Critical speeds of thesystem throughout the full range speed. Wecan find out the natural frequency of thesystem by interpolating.Figure 12: Critical Speed Map by TMS-050Figure 14: Damped Critical Speed by TMS-050Table 9: Damped Natural Frequency(Hz) at the Operating SpeedMode 4.9451.9Figure 13 show damping critical speed(Campbell diagram) of the system fromAnsys. Figure 14 Show critical map methodto find out damped critical speed of thesystem by TMS-050 software.Often, rotor critical speeds correspond tonatural frequencies of the system. Steamturbine rotor is supported by two tilting padbearings. Typically, stiffness and damping345
Int. J. Mech. Eng. & Rob. Res. 2014Nagaraju Tenali and Srinivas Kadivendi, 2014coefficients of the bearing are varied withrotating speed, and in this case, naturalfrequencies of the system are varied. Whena natural frequency equals to the rotatingspeed, the rotating speed is called criticalspeed. TMS-050 series software gives onlynumerical damped natural frequencycorresponding to the speed. In order to findout damped critical speed it is necessary toconvert this numerical data into graphicalrepresentation.Figure 16: Unbalance Response at 2nd BearingTable 10: Damped Critical SpeedCritical SpeedAnsys rpmTMS-050 rpm174517780281458045Figure 17: Unbalance Response at 1st BearingTable No.10 shows damped critical speed oftwo different tools Ansys and TMS-050B. Harmonic AnalysisIn this section, it will show unbalanceresponse of the system at the bearinglocation by apply unbalance force at thecenter position to find out displacementswhich is very sinusoidal at the same knownfrequency and comparison of all the resultwith shop test result.Figure 18: Unbalance Response at 2nd BearingFigure 15: Unbalance Response at 1st Bearing346
Int. J. Mech. Eng. & Rob. Res. 2014Nagaraju Tenali and Srinivas Kadivendi, 2014Figure 19: Unbalance Response at 1st BearingTable 12: Unbalance Response at 2nd BearingBearingLocationAnsysTMS-050Shop Test120.152220.5212.5139.8C. Transient AnalysisIn this section, it will show the response of astructure to arbitrary time-varying loads at thebearing location to find out stability of thesystem at the different operating speed withthe seal effect. If the amplitude of the systemis decrease with time, that means system isstable otherwise system is unstable. Also itwill calculate log-decrement (ld) value of thesystem at different speed and comparison ofthe ld value.Figure 20: Unbalance Response at 1st BearingFigure 21: Response at 1st BearingFigure 15 and Figure 16 shows unbalanceresponse at bearing 1 and 2 respectively fromAnsys , Figure 17 and Figure 18 showsunbalance response at bearing 1 and 2respectively from TMS-050 and Figure 19and Figure 20 shows experimental unbalanceresponse at bearing 1 and 2 respectively.Figure 22: Response at 2nd BearingTable 11: Unbalance Response at 1st BearingBearingLocationAnsysTMS-050Shop Test1101210.5220.52120.12347
Int. J. Mech. Eng. & Rob. Res. 2014Nagaraju Tenali and Srinivas Kadivendi, 2014From Eigen frequency analysis of steamturbine rotor model, Eigen frequencies of thesteam turbine rotor for different rotationalspeeds are calculated. The Eigenfrequencies obtained from Ansys and TMS050 are closed to each other for most of themodes. The number of critical speedscalculated from Ansys model and TMS-050is fair. The Campbell diagram generated fromAnsys is very similar to critical speed diagramof TMS-050.Figure 21 and Figure 22 shows the transientresponse at bearing 1 and 2 respectively.Table 13: Comparison of ld 3182.127118001.9972.447From harmonic analysis, the maximumdisplacement of the rotor and bearing loadfor the applied unbalance loading aredetermined. Peak values of the responsecurves obtained from Ansys and TMS-050are relatively close to each other. And alsothe results obtained from two tools are closedto experimental results. From Transientanalysis, the amplitude of response isdecreases with increase the time, whichmeans system is stable. And also logdecrement values calculated from Ansys andTMS-050 are good agreement.Table 14: Comparison of ld 5631.213118001.1320.918Finally the results obtained from variousanalyses are under acceptable limits. So thesystem is safe for working with given bearingvalves and rotor loads. Apart from these,Ansys software could be an effective tool forrotor dynamics calculation in many aspects.It has got some extra additional features thanTMS-050. Ansys has the capability to handlemore complex geometry.Table 13 and Table 14 Show thecomparison of the ld value at differentoperating speed for the 1st dampedfrequency and 2nd damped frequencyrespectively.CONCLUSIONThus the main objective of the thesis work tobuild and to perform Rotor dynamics analysisof steam turbine rotor model using Ansys isaccomplished. It has been shown throughsimulations and comparisons, the resultsobtained from Ansys model and TMS-050 arein good agreement with each other. The Rotormade of multiple steps is modeled andanalyzed for different boundary conditions inAnsys and TMS-050. The analysis summaryis as follows.REFERENCES1.348Abdul Ghaffar Abdul Rahman, Ong ZhiChao and Zubaidah Ismail (2010),“Effectiveness of Impact-SynchronousTime Averaging in determination ofdynamic characteristics of a rotordynamic system”, Technical Paper,Mechanical Engineering Department,Faculty of Engineering, University ofMalaya, 50603 Kuala Lumpur, Malaysia.
Int. J. Mech. Eng. & Rob. Res. 20142.Nagaraju Tenali and Srinivas Kadivendi, 20149.Carlos Alberto Bavastri, Euda MaradaSilva Ferreira, José João de Espíndolaand Eduardo Márcio Oliveira Lopes(2005), “Modeling of Dynamic Rotorswith Flexible Bearings Using ViscoelasticMa
ROTOR DYNAMIC ANALYSIS OF STEAM TURBINE ROTOR USING ANSYS Nagaraju Tenali 1* and Srinivas Kadivendi *Corresponding Author: Nagaraju Tenali, tenali.n1830@gmail.com Rotor dynamics is a field under mechanics. Mainly deals with the vibration of rotating structures. In recent days, the study abo
17 Table 3. Compressed Water and Superheated Steam (continued) 0.01 MPa (ts 45.806 C) 0.02 MPa (t s 60.058 C) 0.03 MPa (t s 69.095 C) v ρh s t, C v h s t, C v ρ h s 26 446. 0.037 814 3076.7 9.2827 300 13 220. 0.075 645 3076.5 8.9625 300 8811.0 0.113 49 File Size: 630KBPage Count: 60Explore furtherCalculator: Superheated Steam Table TLV - A Steam .www.tlv.comSuperheated Steam Tables - Gilson Enggilsoneng.comCalculator: Superheated Steam Table TLV - A Steam .www.tlv.comCalculator: Superheated Steam Table TLV - A Steam .www.tlv.comSteam Table Calculator Superheated Steam Region Spirax .www.spiraxsarco.comRecommended to you b
flapping trimmed to near zero. The new RTA rotor balance provides increased accuracy in measuring rotor hub loads. This five- component rotor balance measures rotor lift, drag and side forces, together with the rotor pitching and rolling moments. Also incorporated was an instrumented flex coupling to measure rotor torque.
Salient pole rotor - the individual rotor poles protrude from the center of the rotor, characterized by concentrated windings, non-uniform air gap, larger rotor diameters, used in applications requiring low machine speed and a large number of machine poles (example - hydroelectric generation). 2. Cylindrical rotor - the individual rotor poles .
Page 2 June 2015 Large Steam Power Plants Siemens Steam Turbines for coal-fired Steam Power Plants Power output 120 MW to 700 MW Max. steam parameters Main steam / Hot reheat steam 177 bar / 600 ºC / 620 ºC 2,570 psi / 1,110 ºF / 1,150 ºF SST-5000 series for coal-fired Steam Power Plants
2. Use the steam control to select the amount of steam you want ( m low, l medium, h high). 3.2. Squeeze the steam button to produce steam, release it to stop. NOTE: When you first start your steam station and pull the trigger for steam ironing, there will be a delay as your steam station pumps water from the reservoir to prime the system.
Superheated steam. Properties of steam Steam tables To determine various steam properties . Guide C, Section 4, or IOP Guide Total enthalpy, h g h f (sensible) h fg (latent) Extract from superheated steam tables Extract from the saturated steam tables Do you know how to use these steam tables? (Source: www.spiraxsarco.com .
525.9835Hz, then using Ansys software to carry out modal analysis under different constraint conditions on the screw rotor, the modal that the screw rotor, the finite element . frequency using the dunkery calculation method of screw rotor. The screw rotor shaft is simplified as uniform, equivalent diameter d by type (1) to calculate 36.6mm, d i
Annual Women’s Day Sunday, August 24 Congratulations on a blessed Youth Day!! Enjoy your break during the month of August. Women’s Day Choir Rehearsals July 31, August 7, 14, 19, 21 . Beginners Andrew Ash Chaz Holder Primary Deion Holder Nia Belton Junior William Ash Deondrea Belton Intermediate RaShaune Finch Jaylin Finch Advanced Rayanna Bibbs Tavin Brinkley Adult #2 Jeffry Martin Joseph .
