Long blade vibration model for turbine-generator shafts torsional vibration analysis

Similar documents
PRACTICAL CONSIDERATIONS IN HUMAN-INDUCED VIBRATION

Chapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II

Functions of Random Variables

2C09 Design for seismic and climate changes

A New Method for Solving Fuzzy Linear. Programming by Solving Linear Programming

Ahmed Elgamal. MDOF Systems & Modal Analysis

A Method for Damping Estimation Based On Least Square Fit

Dynamic Analysis of Axially Beam on Visco - Elastic Foundation with Elastic Supports under Moving Load

Stationary states of atoms and molecules

FREQUENCY ANALYSIS OF A DOUBLE-WALLED NANOTUBES SYSTEM

CH E 374 Computational Methods in Engineering Fall 2007

DYNAMICS. Systems of Particles VECTOR MECHANICS FOR ENGINEERS: Seventh Edition CHAPTER. Ferdinand P. Beer E. Russell Johnston, Jr.

best estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best

ON THE MOTION OF PLANAR BARS SYSTEMS WITH CLEARANCES IN JOINTS

7.0 Equality Contraints: Lagrange Multipliers

Lecture 07: Poles and Zeros

MOLECULAR VIBRATIONS

d dt d d dt dt Also recall that by Taylor series, / 2 (enables use of sin instead of cos-see p.27 of A&F) dsin

Some Different Perspectives on Linear Least Squares

2/20/2013. Topics. Power Flow Part 1 Text: Power Transmission. Power Transmission. Power Transmission. Power Transmission

Simple Linear Regression

ENGI 4421 Propagation of Error Page 8-01

KURODA S METHOD FOR CONSTRUCTING CONSISTENT INPUT-OUTPUT DATA SETS. Peter J. Wilcoxen. Impact Research Centre, University of Melbourne.

Non-degenerate Perturbation Theory

The Finite Element Analysis of Magnetic Field of the Hybrid Excitation Linear Synchronous Motor with Permanent Magnet and Electric Excitation

Analysis of Variance with Weibull Data

Evaluation study on training of 100m sprint athletes based on improved D-S evidence theory

RECURSIVE FORMULATION FOR MULTIBODY DYNAMICS

MEASURES OF DISPERSION

[ L] υ = (3) [ L] n. Q: What are the units of K in Eq. (3)? (Why is units placed in quotations.) What is the relationship to K in Eq. (1)?

A New Family of Transformations for Lifetime Data

Estimation of Stress- Strength Reliability model using finite mixture of exponential distributions

ANALYSIS ON THE NATURE OF THE BASIC EQUATIONS IN SYNERGETIC INTER-REPRESENTATION NETWORK

Transforming Numerical Methods Education for the STEM Undergraduate Torque (N-m)

Comparison of Analytical and Numerical Results in Modal Analysis of Multispan Continuous Beams with LS-DYNA

Analytical Study of Fractal Dimension Types in the Context of SPC Technical Paper. Noa Ruschin Rimini, Irad Ben-Gal and Oded Maimon

CS 2750 Machine Learning. Lecture 7. Linear regression. CS 2750 Machine Learning. Linear regression. is a linear combination of input components x

STRUCTURAL SYSTEM IDENTIFICATION BASED ON SUBSYSTEM ANALYSES

A Conventional Approach for the Solution of the Fifth Order Boundary Value Problems Using Sixth Degree Spline Functions

{ }{ ( )} (, ) = ( ) ( ) ( ) Chapter 14 Exercises in Sampling Theory. Exercise 1 (Simple random sampling): Solution:

The Geometric Least Squares Fitting Of Ellipses

Engineering Vibration 1. Introduction

Standard Deviation for PDG Mass Data

( t) ( t) ( t) ρ ψ ψ. (9.1)

Derivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations

L5 Polynomial / Spline Curves

UNIT 7 RANK CORRELATION

Physics 114 Exam 2 Fall Name:

Journal of Chemical and Pharmaceutical Research, 2014, 6(7): Research Article

Beam Warming Second-Order Upwind Method

Lecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem

Research Article An Improved Fourier Series Method for the Free Vibration Analysis of the Three-Dimensional Coupled Beams

EECE 301 Signals & Systems

Chapter 8. Inferences about More Than Two Population Central Values

Introduction to local (nonparametric) density estimation. methods

Lecture 8 IEEE DCF Performance

3.1 Introduction to Multinomial Logit and Probit

THE MODELING AND SIMULATION OF A ROBOTIC ARM

PHYS Look over. examples 2, 3, 4, 6, 7, 8,9, 10 and 11. How To Make Physics Pay PHYS Look over. Examples: 1, 4, 5, 6, 7, 8, 9, 10,

FEMA 451B Topic 4 Handouts Structural Dynamics of MDOF 1. Symbol Styles Used in this Topic M U. Matrix or vector (column matrix) m u

A New Mathematical Approach for Solving the Equations of Harmonic Elimination PWM

BERNSTEIN COLLOCATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS. Aysegul Akyuz Dascioglu and Nese Isler

The internal structure of natural numbers, one method for the definition of large prime numbers, and a factorization test

Can we take the Mysticism Out of the Pearson Coefficient of Linear Correlation?

MODELING AND OPTIMIZATION OF TWO-LAYER POROUS ASPHALT ROADS ABSTRACT 1- INTRODUCTION

Ordinary Least Squares Regression. Simple Regression. Algebra and Assumptions.

44 Chapter 3. Find the 13 term and the sum of the first 9 terms of the geometric sequence 48, 24, 12, 6, 3, 3 2 Solution 2

G S Power Flow Solution

Capacitated Plant Location Problem:

STIFFENED MULTI-BAY COUPLED SHEAR WALLS

3D Reconstruction from Image Pairs. Reconstruction from Multiple Views. Computing Scene Point from Two Matching Image Points

Chapter Business Statistics: A First Course Fifth Edition. Learning Objectives. Correlation vs. Regression. In this chapter, you learn:

Multi Objective Fuzzy Inventory Model with. Demand Dependent Unit Cost and Lead Time. Constraints A Karush Kuhn Tucker Conditions.

A Robust Total Least Mean Square Algorithm For Nonlinear Adaptive Filter

The Study on Direct Adaptive Fuzzy Controllers

Remote sensing image segmentation based on ant colony optimized fuzzy C-means clustering

Lecture 16: Backpropogation Algorithm Neural Networks with smooth activation functions

RELIABILITY OF THE IEC (1995) STANDARD FOR DETERMINING BREAKDOWN VOLTAGE OF ELECTROINSULATING OILS

DEMONSTRATION AND METHOD FOR CALCULATING THE EFFICIENCY OF DIFFERENTIAL MECHANISM

Solving Constrained Flow-Shop Scheduling. Problems with Three Machines

A NEW FINITE ELEMENT CONSIDERING SHEAR LAG

DISPLACEMENT-BASED SEISMIC DESIGN OF MDOF CONTINUOUS GIRDER BRIDGES

SUBCLASS OF HARMONIC UNIVALENT FUNCTIONS ASSOCIATED WITH SALAGEAN DERIVATIVE. Sayali S. Joshi

On Submanifolds of an Almost r-paracontact Riemannian Manifold Endowed with a Quarter Symmetric Metric Connection

A Variable Structure Model Reference Adaptive Control For MIMO Systems

A Penalty Function Algorithm with Objective Parameters and Constraint Penalty Parameter for Multi-Objective Programming

Basic Concepts in Numerical Analysis November 6, 2017

An Innovative Algorithmic Approach for Solving Profit Maximization Problems

12.2 Estimating Model parameters Assumptions: ox and y are related according to the simple linear regression model

D. L. Bricker, 2002 Dept of Mechanical & Industrial Engineering The University of Iowa. CPL/XD 12/10/2003 page 1

Comparison of Dual to Ratio-Cum-Product Estimators of Population Mean

Application of Calibration Approach for Regression Coefficient Estimation under Two-stage Sampling Design

Symmetry of the Solution of Semidefinite Program by Using Primal-Dual Interior-Point Method

Quantization in Dynamic Smarandache Multi-Space

Log1 Contest Round 2 Theta Complex Numbers. 4 points each. 5 points each

Bounds on the expected entropy and KL-divergence of sampled multinomial distributions. Brandon C. Roy

Analysis of Lagrange Interpolation Formula

residual. (Note that usually in descriptions of regression analysis, upper-case

The theoretical background of

Transcription:

Avalable ole www.ocpr.co Joural of Checal ad Pharaceutcal Research, 05, 7(3):39-333 Research Artcle ISSN : 0975-7384 CODEN(USA) : JCPRC5 Log blade vbrato odel for turbe-geerator shafts torsoal vbrato aalyss ezheg J ad Yuog Gu North Cha Electrc Power Uversty, Chagpg Dstrct, Beg, Cha ABSRAC he ubrella type vbrato of the last stage blade dsk of the low pressure turbe ca be duced by the exteral exctato whch caused by the torsoal vbrato of the turbe-geerator shafts. However, there s o applcable odel for aalyzg the blade dsk vbrato due to shft torsoal vbrato. Accordg to the characterstc of ubrella type blade dsc vbrato, a blade odel for aalyzg the vbrato of the log blades durg the torsoal vbrato of the shafts s gve ths paper. A paraeter adustg ethod based o sestvty aalyss s preseted so as to ake the heret characterstc of the blade dsc odel close to that of the real blade dsc. Keywords: Blade dsc vbrato odel; Paraeter Adustet; urbe-geerator; orsoal vbrato. INRODUCION At preset, splfed odel for shaft-blade cobed vbrato s to cosder the shafts of the stea turbe geerator ut as a ass-sprg odel. I ths case, the turbe blades are cosdered as the brach structures of the shafts odel ad ther radal oveets are gored. o accout the teracto of the blades durg vbrato, the fluece of the shroud has bee aalyzed [-5], ad the shroud s cosdered as sprgs [6]; or cosdered as dapg [7-8]. Durg shaft lateral vbrato, blades dfferet posto have dfferet vbrato states, therefore, these kds of odel ca effectvely sulate the teractos aog dfferet blades durg lateral vbrato. However, durg torsoal vbrato, blades dfferet posto receve the sae exteral torque fro the turbe shft. herefore, the characterstc of the last stage blade dsk vbrato due to shaft torsoal vbrato s ubrella type vbrato wthout radus. here s o relatve oveet aog the blades these two odels uder ubrella type vbrato. Hece, these two odels caot stulate the pact caused by the shroud o the blade durg shaft vbrato. herefore, t s ecessary to establsh a ew vbrato odel for sulatg the vbrato of the log blades durg the torsoal vbrato of turbe-geerator shafts. Log Blade Vbrato Model o aalyze the pact caused by the blade shroud o the log blades durg torsoal vbrato aalyss, a ass-sprg odel ca be establshed accordg to the followg ethods: durg shaft torsoal vbrato, the blade vbrato ca be cosdered as lateral vbrato relatve to the blade root. he, each blade ca be cosdered as a ass-sprg odel wth oe ed fxed o the shaft, whch has oly x drecto vbrato. he ass-sprg odel ca be show the followg Fg.. Fg. Blade vbrato equvalet odel 39

ezheg J ad Yuog Gu J. Che. Phar. Res., 05, 7(3):39-333 Paraeter Adustet. urbe blades ca be odeled to luped ass torsoal vbrato odel wth degrees of freedo, the u-daped free torsoal vbrato dfferetal equato s as follows: I θ & + Kθ 0 & () Where I ad K are the oet of erta atrx ad stffess atrx, respectvely; θ ad && θ are ass dsplaceet atrx ad accelerato atrx, respectvely. Substtutg θφs (ωt+ψ) to Eq. (), the atural frequecy ω ad vbrato odes φ ca be obtaed fro: ( K I ) φ 0 ω () Based o vbrato echacs, the relatoshps betwee atural frequecy ad vbrato odes are as follows: K ω I φ (3) ( ) 0 [ I ] φ φ (4) Where ω ad φ are the atural frequecy ad vbrato ode after regularzato of the tch order, respectvely. Multply Eq. (3) by φ fro left ad dfferetatg t, followg relatos are obtaed: ϕ I ω ϕ K I + I + K I 0 (5) x x x x x ( ω ) ϕ ϕ ( ω ω ) ϕ ϕ ( ω ) Isertg Eqs. (3) Ad Eq. (4) to Eq. (5), the atural frequecy sestvty of structure paraeter s obtaed by: φ ω x x I ω x ω φ (6) Where x stads for the elastc stffess K. he sestvty of ω to the th sprg s gve by: φ φ ω K ω [( φ ) ( φ ) ] ω + (7) All the sestvtes of each order atural frequeces ca be calculated based o Eq.(7). he value of the sestvty reflects the rate of chage of atural frequecy wth the torsoal rgdty varato. he atural frequecy of torsoal vbrato ω ca be expressed by the oet of erta I ad torsoal stffess K as follows: ( I I K L K ) ω f (8) L, Wth aylor seres expaso Eq. (8) ca be wrtte whe gorg the secodary ad ts above odfed ters: ω ω ω ωo ω I + K (9) I Eq. (9) descrbe the relatoshp betwee the varato of torsoal vbrato atural frequecy ω ad structure paraeters I K. he relatoshp s cosdered to be lear due to the the varato of structure paraeters are 330

ezheg J ad Yuog Gu J. Che. Phar. Res., 05, 7(3):39-333 relatvely sall. So the secodary ad secodary above odfed ters the aylor expaso were gored. he torsoal oet of erta has bee olded accurately eough ad the dapg coeffcet has lttle fluece o the torsoal vbrato heret characterstcs whe adustg the torsoal vbrato odel. herefore, oly cosderg the torsoal rgdty varato Eq. (9) ca be splfed as: ω K ω (0) orsoal vbrato atural frequecy devato ω ca be acheved based o the actual value ω by otorg ad aalyss ad the calculato result ω o by the torsoal vbrato odel. Ad the the torsoal rgdty of adaptve adustet quatty K ca be solved by a equato set based o Eq. (0), whch ca ake the orgal vbrato odel becoe ore accurate. he equato set ca be wrtte as atrx expresso [A] [X] [B]: ( ) ω L ω B () ( ) X () K L K ω ω L A M L M (3) ω ω L he soluto of ths equato set whch s torsoal rgdty adustet quatty K desred to kow ca be obtaed fro atrx trasforato equatos, that s XA/B. Fro the above, the paraeter adustg procedure s show Fg.. Fg. he logc schee of torsoal vbrato odel ole adaptve adustet Aalyss ad Applcato ake a last stage blade the low pressure turbe of a 000MW turbe-geerator ut as a research obect. Assue the actg stress at the blade tp posed by the shroud were proportoal to the dsplaceet of blade tp relatve to ts root, show Fg.3. Based o fte eleet aalyss results of the ature characterstcs of blade vbrato, the ature characterstcs of ass-sprg odel ca be show Fg.4 ad table. Fg.3 Mass-sprg equvalet odel 33

ezheg J ad Yuog Gu J. Che. Phar. Res., 05, 7(3):39-333 0.5 x 0-0.5 - -.5 ode(orgal) ode(orgal) ode(strg-mass Model) ode(strg-mass Model)..4.6.8 R() Fg.4 Mass-sprg odel vbrato odel Let kg able Mass-sprg vbrato odel ature characterstcs ω ω / 07.743H z 6 0 / 403.34H z ω ω ( ) X ( X) ( X) ( X) 3.7436 0.398 -.667 k N calculate ass-sprg odel paraeters through paraeter odfcato ethod show the fg. Ad the calculato value s show the followg table. k 3 able ass-sprg odel calculato value ω ω ( ) k / X ( X) 3 0.5 3.73 0.3660 -.3660 0.93-0.36-0.040.93.8674 0.540 4.75 0.3775 -.936 0.8 0.949-0.093.404.063 0.5533 4.586 0.356 -.8 0.944 0.393-0.07.5958.306 0.5660 3.94 0.399 -.678 0.0996 0.393-0.0005.6954.4409 0.5665 3.7365 0.39 -.63 It ca be see fro the table, after 4 tes of calculato, the odel ature frequecy rato odel s alost the sae wth blade vbrato characterstcs. Ad the, t s oly requred to adust ad k to ake the odel rotatoal erta relatve to shaft ceterle be the sae wth that of actual blade. Ad the st phase ature frequecy of the odel s the sae wth that of blade. So, the ass paraeter ad rgdty paraeter of sgle blade odel ca be calculated. I-addto, oe crcle has 85 blades. he paraeters of the whole crcle blade ass-sprg odel ca be calculated. he ass ad rgdty values of sgle blade ad oe crcle blades s show the followg table 3. able 3 Mass-sprg odel paraeter k k k3 3.35kg 84.75k g 9.0kg 766.70k g.343 07N/.46 09N/ CONCLUSION 3.783 07N/.7866 09N/ 7.6087 06N/ 3 6.4674 08N/ It has bee proved that the heret characterstcs of the log blade ubrella type vbrato due to torsoal vbrato of turbe-geerator shafts ca be accurately sulated usg the blade odelg ethod ad the paraeter adustg ethod. Whch ca sure the accuracy of the blade vbrato respose aalyss. REFERENCES []D. Cha, A.Sha. Joural of Egeerg for Gas urbes ad Power, 003, 5 (4): p.04-049. [] K. H. Koh, J.H.Grff, S.Flpp, et al. Gas urbes Power, 005, 7 (4): p.856-86. 33

ezheg J ad Yuog Gu J. Che. Phar. Res., 05, 7(3):39-333 [3] Eder Cgeroglu, H.Nevzat Ozguve. Joural of Soud ad Vbrato, 006, 95 (3-5): p.08 043. [4] Aldo A.Ferr, Waye E.Whtea. Joural of Soud ad Vbrato, 007, 306 (3): p. 400 48. [5]Eder Cgeroglu, NgA, Cha-Hsag, et al. Nolear Dyac, 007, 50. [6]J.H.Grff,.M.Hoosac. Joural of Vbrato, Acoustcs, Stress ad Relablty Desg. 984, 06(): p. 04~0. [7]S..We, C.Perre. Joural of Vbrato, Acoustcs, Stress ad Relablty Desg.988, 0(4): p.49~438. [8]S..We, C.Perre. Joural of Vbrato, Acoustcs, Stress ad Relablty Desg.988, 0(4): p.439~449. 333