(International Centre For Mechanical Sciences 297) Neville F. Rieger (Eds.) - Rotordynamics 2 - Problems in Turbomachinery-Springer-Verlag Wien (1988) PDF
(International Centre For Mechanical Sciences 297) Neville F. Rieger (Eds.) - Rotordynamics 2 - Problems in Turbomachinery-Springer-Verlag Wien (1988) PDF
ROTORDYNA MICS 2
PROBLEMS IN
TURBO MACHINERY
EDITED BY
N. F. RIEGER
STRESS TECHNOLOGY INCORPORATED
Neville F. Rieger
Rochester, New York
CONTENTS
Page
Preface
Chapter I
Introduction
by N.F. Rieger ................................................................ . I
Chapter 2
Bearing Properties
by E. Kramer ................................................................. 17
Chapter 3
Analysis of Rotors in Bearings
by E. Kramer ................................................................. 41
Chapter 4.1
Principles of Balancing and of Balancing Machines
by N.F. Rieger ................................................................ 67
Chapter 4.2
Flexible Rotor Balancing
by N.F. Rieger ................................................................ 95
Chapter 4.3
Case Histories in Balancing of High Speed Rotors
by N.F. Rieger ............................................................... 129
Chapter 5
Seal Properties
by R. Nordmann ............................................................. 153
Chapter 6
Stability of Rotors
by R. Nordmann ............................................................. 175
Chapter 7
Computer Analysis of Rotor Bearings- P.A.L.L.A.: A Package to Analyze the
Dynamic Behavior of a Rotor-Supporting Structure System
by G. Diana, A. Curani, B. Pizzigoni .......................................... .. 191
Chapter 8.1
Sensor Technology
by J. Tonnesen ............................................................•. . 261
Chapter 8.2
Experimental Techniques for Rotordynamics Analysis
by J. Tonnesen ............................................................... 269
Chapter 9
Interaction between a Rotor System and its Foundation
by L. Gaul . .........................................................•..•..... 283
Chapter 10.1
Problems of Turbine Generator Shaft Dynamics
by D. W. King, N.F. Rieger ...••...•..•......•.••.......•...................... 307
Chapter 10.2
Torsional Systems: Vibration Response by Means of Modal Analysis
by P. Schwibinger, R. Nordmann ............................. .................. 331
Chapter 10.3
Torsional Dynamics of Power Transmission Systems
by N.F. Rieger ............................. ............................. ..... 359
Chapter 11.1
Free and Forced Vibrations of Turbine Blades
by H. lrretier . ............................. ............................. ...... 397
Chapter 11.2
Flow Path Excitation Mechanisms for Turbomachine Blades
by N.F. Rieger ............................. ............................. ..... 423
Chapter 11.3
The Diagnosis and Correction of Steam Turbine Blade Problems
by N.F. Rieger ............................. ............................. ..... 453
Chapter 11.4
An Improved Procedure for Component Life Estimation with Applications
by N.F. Rieger ............................. ............................. ..... 485
Chapter 11.5
Damping Properties of Steam Turbine Blades
by N.F. Rieger ............................. ............................. .... . 515
Chapter 12.1
Magnetic Bearings
by G. Schweitzer ............................. ............................. ... 543
Chapter 12.2
Vibrations in Variable Speed Machines
by H. lrretier . ............................. ............................. ...... 571
CHAPTER I
INTRODUCTION
N.F. Rieger
Stress Technology Incorporated, Rochester, New York, USA
ABSTRACT
1.2 Definitions
1.6 References
1.1 Dimentberg Per gammon Classical shaft dynamics for turbine 1 1961
generator applications.
Whirl
Motion
Figure 1 Circular Orbit
New Orbit
Original
Orbit
Figure 5 Transient Growth to New Stable Orbit
Impulse Orbi't
Growth
BEARING PROPERTIES
E. Krimer
Technische Hochschule Darmstadt, Darmstadt, FRG
ABSTRACT
Symbols
F Journal load
F Static journal load
s
L Bearing length
R Bore radius
s Duty parameter
So Sommerfeld number
dik Uamping coefficient
e lccentricity
h Film thickness
h0 Film thickness at centric journal
18 E. Kramer
/ ' /
~L__j
(2.1)
z
w
and (2.4)
(2.5)
(2.6)
1
u = -2n -(y
ap 2 u
ax -hy) +-hy
(2. 7)
1ap 2-hy)
w =- -(y
2n az
respectively
aQ~ aQ~ ah
-+-+-
ax az at = 0 (2.8)
h -t~ 3 a 1
Q1 = Judy = ---- ~+-Uh
x 0 12n ax 2
(2.9)
h -h3 a
Q• = Jwdy = ~ iP
z 0 !~n az
a h3 ~ + a h3 ap = l}.ah + E!!
ax(n ax) az(n a-z) 6 ax 12 at (2.10)
X
Fig. 2.4 Film forces dF, F and
its components
24 E. Kramer
So= pi
nn
(2.18)
(2.21)
and ( 2. 27)
where the dimensionless terms yik' Sik depend on the bearing para-
meters and are functions of eccentricity ratio £. For the short
circular bearing they are derived in chapter 2.4.3.
2 -
dp 2= ~[e( r.-2·() sin {\p-y) -2e cos (~y )1 (2.29)
dz h
by integration to
(2.30)
Integration within the limits -L/2~z;+L/2 gives with (2.28) the
film load per unit length
. .
-L3_., e:(1-2*) sin {\p-y)-2Ticos{\p-y) - 3
q ('P· t) =- Tl I' _ ___;•;;....•- - - - - - , - - - =-~ q'('P,t). (2.31}
~ [ 1-e:cos {\p-y) 1 26
= - 30
nl :::sin(y)--:)
q('P} - 252 [1-e:cos{\p-y)] 3 ( 2' 32 )
This theoretical load (2.32) is plotted in figure 2.7, showing
regions with positive and negative loads. In reality negative
loads are impossible and therefore they will 'be assumed as zero.
28 E. Kramer
Uith this assumption and with (2.12), (2.13}, (2.14) and (2.31)
the components of the resulting force are
'( y
F1 = -K : q'(4),t)costp14), F2 = -K f q'(4),t}sirA4kf'P (2.33)
--:-+y -!+'(
- 3
with
-;-r
K - ~L R~ (2.34)
are the components with respect to the X', Y' system. The latter
may be splitted in a static part which depends on the eccentrici-
ty ratio E and in a dynamic part which is a function of the velo-
cities ~.~ and also of E:
Bearing Properties 29
0.1
0 goo
0.2
E
0,4
j 0,6
0.8
100
lgo· 10
ex So
60" pr1
30" 0.1
o· 0.01
F1 = K(Ficosy- Fbsiny)
} (2.46)
F2 = K(Fisiny+ Fbcosy) ,
(2.47)
.2£
ax = -kos
o Y
1
(2.48)
ay = -s1ny
-1 . ay _ 1
-ax
1 e , ax- ecosy
2
leads finally to
82 Fs Fs
kik = 'SOfik(£) 6 = Y;kT
with
82
Yi k = So f i k (£ ) (2.49)
32 E. Kramer
and
2s
--2""'2 [1 +a ( s) 1
( 1- s )
(2.50)
1 . 4 2 2,
f 21 (s) =--z 2 S/2"1+[(-::-) a(s)+21E;
(1-s )
2E , 4 2
f 22 (s) = - -,..--,.;
2 3 ·_1+[(-=) a(c.)+2]t. 2}
(1-E )
where
1
a ( s) = --4---=-2-~2
1+[(;:) -1]t:
Figure 2.11 shows the dimensionless stiffness coefficients yik as
a function of So;e 2 or s.
t~o~~----------~----
l;k
~+---~~~----~~~~
aF .
_1 =
aF. _oe:_
_l
~• aF. ~ •
+ _ l 2L
axk a£ axk ay axk
which leads in similar way as before to the damping coefficients
with
(2.51)
and
TI 4 2 2
2 3/2{1-[2(n) -3]e: a(e:)}
(1-e: )
4e: 4 2 2
g12(e:) = g21 (e:) =- 2 2{1-[2(:rr) -3]e: a(e:)} (2.52)
(1-e: )
TI 4 2 2+e: 2 2
g22(e:) = 2 3/2[ 1+(1T) :-ze: a(e:)]
(1-e: ) 1-e:
t~o+-~~--~---+--~~
~ik
IDf---~-~~~~--J--
- ~12 = -f3z, ~,
0.1+---~--~---+-----1-
x1 = hll FS 1 + h12 FS 2
} ( 2. 53 I
x2 = h21 FS1+h22 FS2
-1
with [hik 1 = [kik 1
X X
tions, the broken lines show the linearized one. The arrows indi-
cate the direction of Fa. It may be stated that up to F~0,5 Fs
the linear solution is rather good.
Ca X
Limit Circle
Co
The figure shows that the orbit is nearly circular for small
eccentricities and shaped like a banana for large eccentricities.
For small eccentricities the bearing behaves nearly isotropic,
but the displacements are nearly vertical to the additional force.
At high eccentricities the bearing is stiff in radial direction
and rather elastic in tangential direction. Roughly the stiffness
is in the order of Fs/o.
The theory given in chapters 2.1 to 2.3 is not only valid for
circular bore but also for other cross-sections. Furthermore par-
tial fil~ areas are possible. Special bearing types with different
geometry and also with flexible pads are used with respect to
36 E. Kramer
a ~
References
E. Krimer
Technische Hochschule Darmstadt, Darmstadt, FRG
ABSTRACT
3. 1. 1 Jeffcott Rotor
z1 = y1 + e cos~
z2 =y 2 +esin'9
mz 1 = -ky 1 - dy 1
mi 2 =- ky 2 - dy 2 - mg (3.2}
my 1 + dy 1 + ky 1 = men 2cosnt
my2 + dy 2 + ky 2 = men 2sinnt- mg
(3.6}
44 E. Kramer
:. d quantity of damping
= 2m
constants, determined by the
initial conditions
y1 = yCOS (rlt-E)
y2 = ysin (rlt-E)-y 0
where y = en 2£(1-n 2 ) 2 + (2Dn) 2 ( 112 = eV
tan E = ~ (3.10)
1-n
static displacement
Analysis of Rotors in Bearings 45
n =-
n
wk
(3.11)
D = _d_ damping factor.
2{km'
According to (3.9) for constant angular velocity r. the shaft
center describes a circular orbit with center (0, -y 0 ) and with
constant angle E (figure 3.4). The relative amplitude y/e is given
by the wellknown amplification function V [(3.10), figure 3.5)].
It shows that the amplitude y is maximum at 0-::::wk. The correspon-
ding speed
n
cr = ncr /2rr::::wk/2rr (3.12)
t
Y;
Fig. 3.4 Orbit of shaft
center in case
of unbalance excitation
y, ....
-H+-H--- 0: -----
-- 0
-0.1
,0,15
0.2
- 0.25 ---------
0.3
Fig. 3.5
Amplification function
o1-Ld::::::::=±=~r
0 2 3
..,
Fig. 3.6
Position of Wand S for
different angular velo-
cities at t=O
Yz Yz Yz
y,
w Fig. 3.7
Undamped rotor.
Qcwk Q>wk Q:oo Possible configurations
Analysis of Rotors in Bearings 47
and by (3.10)
therefore ::::J_.r£
n
cr 2n v-y: (3.13)
1,13. ~
n
cr
:::
2n v-y: (3.14)
The geometry of real rotors lies between these two models. There-
fore their first critical speed is limited by
U - T (3.16)
max max
Assuming the·first natural mode to be approximately equal to
the static mode, given by displacements yi, then the maximum
potential energy
where w(t) =io.t, which means that the exciting frequency grows
linearly. The figure shows that the envelope of amplitudes is
slightly undulated.
1 o=o.o1; t=o.o161
t 10
.!.
Jr"'
xs / 1' 1 1\1"~
4f(
5 ,'
v\11.01v
0 7'
v 0.5' 1. 2.0
t
~
\i T-
-5
-10
-e I D= o.nl
I
10 , .. ·- j
1.5 2
120
y
e
1 o=o.o2 1
10
0 0,5 15 2
T---
I w-'
~ 10- 3 Fig. 3.11
5·10- 4 Relative acceleration as
function of starting time
and critical speed
where
- = en.[(1-n.)
y. 2 2 2 + (2D.n.) 21-1/2
1 1 1 1 1
2D.n.
1 1
tan £i = -:----2
1-n·1
n D.= Nd , i = 1, 2 .
ni = W.
1 1 2 k .m
1
Analysis of Rotors in Bearings 53
2
Fig. 3.13 Orbit of shaft center
in case of unbalance
excitation.
Rotor with unequal
bearing stiffnesses
l"kwdl
forward whirl {
n=wdl (~ first critical speed)
threshold speed
backward whirl
threshold speed
r _f)_
1~
Q•wd2 ,. second critical speed)
forward whirl
'+'- n"wd2
f y. f(t)
Fig. 3.15 Rotormodel with
bearing damping
(3.25)
Analysis of Rotors in Bearings
,\
where X --
WR
WR =fi
dB
dB - ykRm,
kB
K
= kR stiffness ratio.
1
0.5 0.75
0.5
0.' de
0
(3.27)
o·
~ 05 ~
/ , ·' 075 ko
0.5-
0~~~~======~~
0 0.5 1.5 2 s 10 -
ds-
(3.28)
[K(1-n 2 )-0,5n 212+[nd- 8 (1-n 2)1 2
1000
e
10
''"'"I
0,01 ''"""0,1 '"""I1
de
100
c
10
0,1
0,1 10
fig. 3.19 Factor C [Eq. (3.29)] for small related bearing damping dB
ie 3
i
e
,___
·:
:
. t
vertoco }"
honzonto t
1. Rotor
e ·
Fig. 3.20 Unbalance res-
1 ponse of Jeff-
cott rotor with journal
bearings {161)
1 3
Iw.l . . . -
Analysis of Rotors in Bearings 59
Theoretically the rotor has two critical speeds, because the oil
film behaves anisotropically. But usually the critical speeds differ
so little that only one resonant peak appears. The example in figure
3.20 shows two peaks only for the horizontal rotor vibration. This
problem is treated rather comprehensively by J. Glienicke in /14/.
Knowledge of unbalance behaviour of Jeffcott rotor is doubtless
very important. For real rotors however more complicated models are
necessary, which are characterized mainly by:
different cross-sections of rotor along its axis,
more than two bearings,
consideration of additional effects, above all the gyroscopic
effect,
realistic models for the pedestals and the foundation.
// / / /, •· 1 .' , I I I I
5000
cpm
,,2
f 4000 ,,,
f.
3000 ,'()
,9
ta
2000
"'
f6 t
t~, 5
t3
1000
f ,2
1
0
0 2000 4000 rpm
n-
Fig. 3.22 Natural frequencies of a five-bearing rotor /20/
Order
9 ·~ ~ I I Ll1
Fig. 3.23 Natural modes of a
5,6 rPI
3,4 1-L¥l
1·2 ~-I+-1-+--oofC-A----Jt
I
y"
100
Jlm
@
6 zs zs ~
5
50 zs 1-!
0
0 1000 2000 3000
-- rpm
n-
1
A
40
1\
I 1 11,.
11 o or,r1g1
I I I~ bearings
R t
. "d
y pm I VI
I I
I I
20 I ~
0
0 1000 2000 3000 rpm
n-
References
121 Foppl, A.: Ober den ruhigen Gang von. schnell umlaufenden Hange-
spindeln, Der Civiling. 1895, S. 335.
64 E. Kramer
/13/ Greb, J.: Berechnung des Hochlaufs und Ablaufs von Rotoren.
Studienarbeit am FG Maschinendynamik, TH Darmstadt, 1985.
.
Analysis of Rotors in Bearings 65
N.F. Rieger
Stress Technology Incorporated, Rochester, New York, USA
ABSTRACT
4.1..1 Introduction
Multi-Plane Balancing:
Types of Balancers:
Calibration Types:
Rotor Types:
Rotor Supports:
a) Journal supports.
b) Journal pedestals.
c) Pedestal supports.
d) Foundation base.
Rotor Drive:
Foundations:
Electronics:
c) Operational amplifiers.
Plane Separation:
b-a
u~
0: -F 1 + U - C2 = 0:
c-b
l M2 0: U(c-b) - C1 (c-a) ::: 0: c1 = u-
c-a
b-a
lF = 0: -c1 + u - F2 = 0: F2 = u~
Resonant Machines:
f) Transporter,
g) Control room,
h) Overhead crane,
4.1.7 References
Resonant Machines
Amplitude
= 0.0
= 0.15
Hard Bearing
Machines
,....,r-- Soft Bearing H chines
Rotor Speed
Bearing carriage
Inductance Pedes ta 1
transducer Location
ODD
Vertical
1110ti_ons of
carrtage
~re negligible
Flexible members
for horizontal
1110vement
Rigid foundation
Support
Flexibility
E,I
2n
l El
_a
sin wt sin wt d(wt) El!i
wt
wt
21(
L--+--+-+--r-- +-+- fEr sinwt sin Zwt d(wt) = o
wt 0
E,I
zr
f El sinwt sin(2wt+e) d(wt) = 0
L......<-~~+f--+-f-·1---\--'-----<l
wt
I
.I.
b_±_j
Figure 10. Rigid Rotor With Two Unbalance Forces
N.F. Rieger
Stress Technology Incorporated, Rochester, New York, USA
ABSTRACT
>0 i • j
· y • j · dz • 0 {
j
• 0 i • j
"Jere c::., c.•. are natural frequencies for the ith and jth
rctc•r n1o1des tespcctively, andy. andy. are the modal ampli-
tudes. It follows that a flexfble ro~or may be bal'anced by
~li~i~ating the effect of residual unbalance mode-by-mode. As
the d~flection of the rotor at any speed may be represented
as the su~:~ of nl'ious n,odal deflections multiplied by speed-
dependant constunts, a rc·tor which has been balanced for
sruoott operation at each of its critical speeds in accordance
wit}• the abo\·e condition will remain balanced at any other
SJ;et"d.
+ . . .
e (z) + +
+ . . .
U(z)
+ ... >...Y.(z).
1 1
{"0
£ 1 ,. j
S m(z) Y1 (z)Yu(z)dz:
0 • M
1 1 - j
~0 EI(z) \
II
i -
i j
L
Wcent • S m(z)ro2 [y(z) + e(z)](oQ 1)Y 1 (z)dz
0
l
~ m(z)Y 1 (z)y(z)dz - (j i
c" 2 (z)dz
j m(z)Y 1
0 0
lienee:
v .. jo
d2
1
,.t 21 EI(z) [ ::t dz • jo
,.t 2
1
EI(z)
~
LQ i Y~ (z)J
, 2
dz
i•l
in keeping with the second orthogonality relation given
above. The virtual work of the elastic forces may be
expressed as:
Thus,
1
or
( w/ w . ) 2
1n
In o r d e r t o b a 1 a n c e t he r o t. o r td t h d i s c r e t e ma s s e s at
sr·ecifjc locations it is convenient to deterJtine the effec-
tive unbalance at these locations. The deflection from a
concentrated force F or frorr, a uniformly distributed load q
acting over a small length ~ distant c from z = 0 may be
expressed as a Dodal series, vix.,
c+~l2
= ~ q m(z)Y 1 (z)dz
c-~12
P m(c)Y (c)
i
or
-
Thus,
if
U(z) ~ m(z)e(z).
102 N.F. Rieger
•.u
-
or
w
p
a
pn
where 111 = 111 (cos 111t + i sin 111t) = 111 + i 111 as shown in
Figure 3. Now insert a trial weight T fh corre!tion plane 1.
and re-run the rotor at 111 1 • Again read w and f which
correspond to:
r-~~ ·~~
1
1
aA1 •A2 • • • u1 + T
1 1
I WB1 aB1 aB2 •Bn u2
,_
Now
r VAl-
I VAO
1
1
aAl
1
.A2
1
•An
1
u 1 + T- u 1
-u2
··o
1 1 1 1
aB1 .B2 •Bn u2
r~l-
-
1
•p1
1
•p2
-ull
108 N.F. Rieger
Thus
1 1
1 WA1 WAO
·~u • T
-
1
1 WB1 wBO
.,1 T
-
1
wp 1
T
1 < i < p
1 < q < n
Flexible Rotor Balancing 109
1 up
a
pn
82
An
a2
Bn
- \
\
wpn un
110 N.F. Rieger
wem • wem
r + i wi
em • wmo + am1u1 + . . . + amNUN
CD
- wmo +
n•
~ a U
mn n
112 N.F. Rieger
wr • w r + a r ur - a 1 U 1+
em mo m1 1 m1 1
w 1 - w i + a 1 u 1 + a m1 1
r u1r + . . • - wmo
em mo m1 1
(a 1 u 1 + arur).
mnn mnn
Let:
M
s
-I +
-~- ..li. - 0
~w r ~w 1
n n
Flexible Rotor Balancing 113
\m
L Jl a mnr L- wonr + 'n') (• mnr urn - 8 1u1)]
mn n
+ a i L- w
mn on
1
+I (•aitU~ + •r!'nU~J• 0
n
Modal Balancing:
Multi-Mode Balancing:
S u c c e s s f u 1 a p p 1 i c a t i on s o f t he I n f1 u en c e c o e f f i c i en t
technique appear to include a wider range of rotating
machinery types than those which have been reported for other
flexible rotor balancing methods. Known applications
include: a high speed (24000 rpm) pump simulator, a long
helicopter drive shaft (20 ft. unsupported span, six critical
speeds), certain ultracentrifuge applications, several small
steam turbines, and small aircraft gas turbine applications.
Notable features reported in these applications were:
o Convenience of Application
Other Methods:
4.2.10 Comment
4.2.11 References
Inventor:
rrani! c linn
by~~
H.. Attorney
BEARING
COEfFICIENTS
PEDESTAL
COEFFICIENTS
BALANCE
PLANE
r.........
I .........
BALANCE I ""l
PLANE I }
~
PROXIMITY
PROBE
100
ORIGINAL UNBALANCE =
I \
10 I \
_, '\.
L
1.0
..,..,....... / '
/
"""' !'-...
~THREE-PLANE BALANCE.
~ READINGS AT OVERHUNG DISK._+- READINGS TAKEN AT
BEARING AMPLITUDE 'RESPONSE. BEARINGS. BEARING -A-
THREE-PLANE BAL.-
=
I I I ~ :. AMPLITUDE RESPONSE.
I -THREE-PLANE BALANCE. Iii I
=READINGS AT OVERHUNG DISK. - f h ·
---
~OVERHUNG DISK RESPONSE.
---
--- --- -- ---
."'-... _\
... \
""' K ==---
~
1\.. •·•
r- FOUR-PLANE BALANCE.
~READINGS AT OVERHUNG DISK. ' -.... ·. .. .•·
.·
r-- OVERHUNG DISK RESPONSE.
I I I '·'\. \ ..... /
', .. / /
10-3
4000 5000 6000 7000 8000 9000
.10,000
\ \ //
'
11,000
SPEED, RPM
A. I B I c I D r
•'' •
400 HP on special foundations. 0.63 2t
Pumps and compressors. 0.63 15
Small turbines. 0.63
A I a I c 0
Ill Large electric motors. 0.63 5
Turbines and generators on rigid and heavy
foundations. 0.63 5 20
A I a I c I D
Large electric motors, turbines, and generators
IV on lightweight foundations. 0.63 3 10
Small jet engines. 0.63
A I B c D
v
Jet engines larger than category IV. 0.63 2 ID
CHAPTER 4.3
N.F. Rieger
Stress Technology Incorporated, Rochester, New York, USA
ABSTRACT
4.3.1 Introduction
d) .Misalignment:
Hot vs. cold alignment
Coupling forces
Thrust bearing force
Foundation distortion
Multi-bearing support alignment
Positive Indications:
4.3.3 Conclusions
First assume that the rotor has been balanced in its lowest
mode around 2300 rpm. At some higher speed, well below its
second critical speed, it now is deflecting in its second and
third modes simultaneously. The mode shapes are those shown
in Figure 6. The most suitable locations for balancing planes
are at 25 inches and 70 inches from bearing A for the second
mode, and 10 inches, 44 inches, and 90 inch~s from bearing A
for the third mode. Also assume that all correction masses
are to be added at the same radius in these planes.
When adding corrections mostly for the second mode, the total
mass added at any location must be divided between the two
planes. The proportions into which the mass is best divided
were decided by the reasoning given above. The overall
effect of these correction masses on the vibration of the
bearing pedestals is determined by summing their individual
effects in both modes. These sums are, in turn, deduced from
the corresponding magnitude of the response curves given in
Figure 7. The values involved ·are given in Table 1. A
further computation is then required to determine the overall
effect of corrections added mostly for the third mode: (see
Table 2).
4.3.6 Comments
4.3.7 References
Details of Problem:
Investigation:
:b) Impeller bear ins wear rates had been high. The
impeller bearing had been replaced three times
during the previous 12 months.
Case Histories of Balancing of High Speed Rotors 141
Calculations:
4.3.9 Conclusions
4.3.10 Recommendations
90
110
70
30 :.:
Q..
20 f
c:(
10
100
90
110
70
6(j
f.0. ~;!a
50
•o
IZOHe
30
20 1
~
10 '"'!)'
0
"""""-
40 60 80 100 20
0 20
resentation of
Figure 4 · Finite Element
Stiffene Rep End Plate
d Generator
146 N.F. Rieger
I___( r
0
_.. _ - .
10
. __ . __
20
.
_ JO 'II 50 10 lO 10 to 100
-
-----\----
~ _a...~-~-1
A
r ~-
ot,.;,;;;; incht>' tro"' M••lnt A
--, 8
~ ' '
I / co..ecllon NQul.-4 to \
I ~~~
modal delect
\
I \
I '\
I' '·
·' (a) Unbalance Response in First Mode
I'
I r
\ Correcliort requl.-4 to
\ nuftlfy tha aeconcl
\ modal defect
I
-------
_Jo --:-__-_:_ .
A , ,.
~-
''
_J_ ~-_...!!.-_!f.._.....!!...-
D•IIA"Ce "' inches from be•nftg A
I \
I \
'
''
R...li¥e deReolon
of rotor
.,1~
(b) Unbalance Response in Second Mode
R...,ll•• cleftaolon
of rotor
10 to 100
~-_.._,_11
~jt
.o
V•bt,ltO•• ot
.be111ng ,edt51.11 A
Vtbtlllhon of
M11n"o peclesllll 8
"
Figure 8. Explanations for Readings of Vibration
Vectors OA and OB:
a) Equal out-of-phase component, equal
in-phase component.
b) Equal out-of-phase component •. unequal
in-phase component.
c) Unequal out-of-phase component, equal
in-phase component.
d) Unequal out-of-phase component, unequal
in-phase component.
Case Histories of Balancing of High Speed Rotors 149
A A A
c A
lcl
1~~. c
----7ot-----=cO
Be (0!)
0 0 A
lac
\
'\ I '\ \
\ I
\
\ v '\
'\ \\I
'
a\ 8
8 8
Running Speed
2300 rpm
0.015
1 Mode 2
Model
Mode 1
I
I
0.010
Vibr•tion
Mils p-p
0.005
0
0 1,000 2,000 3,000
Speed, rpm
'Jt.~.r..J".t61,.~
,,...,... ...... fU!t·· ·~·-· \1!""~'!....1"-~~
.
=~:-~
"~""'"'" ~.
..._
. .-. ....... .. ~ ....... -.uc-o.•r......,_,...... .,..
~-- ~· ·-
t~
...,
)JS
l.ao..t
, ••
....
.....
·~
~
'~•.rc .~ ···~-~
J
... ,,,. .... .--• . .......................
)S ,..,~ ·~~ t
_.. ..::.··
M...
· '·" .... ._...
• ...... .._ _ _ _ ,...._._. ...,,-._
..... ......,~... --:','- ~ -~
t.®l...,,.,,__,.._~-roo'·~......,....&-,-.,....,r:..,....--.,
.... ·- ...
~- ) ~·
-~. 1r ._.,__..
1\Ar•••. """"" ~·· ·• ._ .(\o.t•{._·~...t-.--1-"'
.
'1-f" ........ --a.
AA'- l, l'.. ~·l ..t••,_, t ... ....._, .....................
.
~ ~
~
/
. , ./ ......,
.1#• •
.
"-:6,.~ A _ _. b .
.......... At• ... _ ...... .#d .._ ...__ • .._,. h ~ 0
~. I\'-··· ~ ~-··..,.; ..
.··
J.- ~tt'u - • • · · • · ........... ,-~..._ · • - ..,.._,... ,....,"•
/ ~-
.... / ..... T.t.,...~"~• ........ ,.._ ..,._,....... o-'.ftb"•· .. ~~_l
• ................ .,.,.,.Yflf
-
.:::~::: .. t:r.vu~ t.:'..... ·~U'"'-"
·~·"'
Drive
Coupling
Bearing
--
Un111 o( unl»..ncc Unu" .,( ua.MlitKc
.!:fht oclde<l, Di,uncC' in inches ift IC'C'OR4 mc-IC'
~ftect dl'ccl in "'iui~~MMtc
np 11 • multiple
o( dw alobncin' wciKht
uoe<l ,....,.oon fill.l. r~:"'ortK~:t:l :c::h, Bcorin~ A lkaunc 8 8C'ann1 A lcAtlftt' £
II B
'l"h.:rcforc, ac•ual \"ti'tnlion C'l' IM'ann~ •"'Uid be -+1·2l;f unus at A and - i·)x
,..,,
o(.:::.
••
.Y
s-o
f.9"
44
(arposil< oi.S.: or ....or) +0·9xf.i•
s-o -O·Sl< !:!!,
5·9
s-o
-ti·ISI<f.ir -,..,.;:;,
y 90 -2·4y +1·3Sy ~1·15y -tz+
-t-Z·lly -1·32,- ~4·S2.'" HIT
ft'TAL PM; IACH MC'tOC
Thereto«, actual wibration of bc'armcs wauW lac • 6·13.o• unico 11 A on4 6·21yocll
CHAPTER 5
SEAL PROPERTIES
R. NordiiWID
University of Kallenlautem, Kailenlautem, FRG
ABSTRACT
TURBOPUMP TURBINE
SEAL LOCATIONS
_------r
__....- PRES SURE DROP ll p = P,- p2
PRESSURE p 2
PRESSURE P,
To derive the pressure around the shaft and then the force-
motion-relationships for the vibrating rotor we are using a bulk
flow model, which was originally derived by Hirs /1/ . Childs /2/
introduced this bulk flow theory for seal elements. The first
basic idea of this theory is, that the fluid velocity distribu-
tions in radial direction are substituted by average velocities.
For a fluid element between the rotor and stator surface (Fig. 3),
located at the axial coordinate Z and the circumferential coor-
dinate 0, the average axial velocity is U (Z,0,t) and the aver~ge
z
circumferential velocity is u0 (z,0,t) . The corresponding pressure
for this locations is p(Z,0,t) and the seal radial clearance
H(0,t). The shaft circumferential velocity is U' = R O.
156 N. Nordmann
U-RQ
PRESSURE
Formulas (1) and'(2) correlate the shear stresses TR, Ts with the
Reynolds numbers, which are defined in parantheses.
't;Rd9dZ
pHRd9
equilibrium equations (3) and (4), which are shown in Fig. 5 to-
gether with the continuity equation (5)
1+m 1+m
0 0
2 }
(3)
+ (4)
Continuity Equation
auz 1 a 1 oH
H F + i ae (H 0e) + Rn at = 0 (5)
158 N. Nordmann
R = HU = .Q...J:!__Q
e v n
H ( e, t) = Ho + e:: H1
p (z,e,t) = Po + e:: p1
( 6)
uz(z,e,t) = 0 zo + e:: 0 z1
u 6(z,e,t) = 0 eo + e:: 0 e 1
6p = p 1-p 2 = {1 + ~ + 2 o} Q
2
y2 (7)
Seal Properties 159
PRESSURE p 2
for the unknowns UZ 1 ' u81 , p 1 /2/, /3/. From the pressure field
solution the reaction forces acting on the rotor due to the cir-
cular shaft motion have to be determined by integration of the
pressure along the seal and in circumferential direction. Finally
the force motion relationship is established. The dynamic system
of a seal can be modeled by a linear system with stiffness, damp-
ing and inertia terms, if small movements about the seal centre
are assumed (Fig. 7).
F m c c y k k y
y yy yy yz yy yz
= + + (8)
F m c k k z
z zz zz zy zz
The main diagonal elements in each of the matrices are equal and
the cross-coupled terms are opposite in sign. The coefficients
are mainly dependent on the pressure drop, the average axial ve-
locity V, the rotational speed n of the shaft the seal geometry
(seal length L and Radius R) and on some quantities characteri-
zing the friction in a seal. It is important to note, that the
cross coupled stiffnesses kyz = - k zy are strongly infllh:uced by
Seal Properties 161
the rotational speed n and the fluid entry swirl, which is the
circumferential velocity of the fluid at the seal entrance. This
effect may cause serious instability problems in high speed rota-
ting machinery, when the cross coupled stiffness terms become
dominant.
PRESSURE DROP 4p
AVERAGE AXIAL VELOCITY v K
AVERAGE FLOW TIME T: L/V
ROTATIONAL SPEED g
SYSTEM : SEAL
INPUT OUTPUT
FORCES DISPLACEMENTS
IDENTIFICATION PROBLEM .
TO FIND THE PARAMETERS OF THE SEAL
FRCJ.t INPUT- OUTPUT- RELATIONS
HODEL MEASUREMENT
F,
ANALYTICAL
~;=F=R=E=OU=E=N=~==R=E=S~==N=SE========~~~~~======~,
1!1:NTlFCATION OF [J'(NAMIC SEAL COEFFICIENTS
RESULTS
RIGID SHAFT IN
SEALS WITH
RIGID BEARINGS
k,)<,,,k;,.ku
c,.c,.c",c,,
m., . mlf
~ = (kyy-w2 (~
2 myy )+iwcyy )(k zz -w2 (~
2 mzz )+iwc zz )
I
-w2m
1
K (w) I K (w) k (=+m )+iwc I k +iwc
yy I yz yy 2 yy yyl yz yz
I I
----1-----
1
=2 -----------.---------- I
• I
K (w) I K {w) k
zy
+1WC
zy
1k -w2(m
-+m ) +1WC
•
zy 1 zz 1 zz 2 zz zz
INPUT OUTPUT
F~ER ANALYZER
"~
.,, ,, ,, / FOURIER RESPONSE TAPE
DISPLACEMENT TRANS FUNCTIONC
FORMATION
PICK UP
-
mobility matrix H•and the complex stiffness matrix K~should be
the unity matrix ~· By combining the measured matrix ~·with the
-
analytical matrix ~·the result will be ~ plus an additional error
matrix §, caused by measurement noise. Fig. 15 points out this
fact and shows either the complex equation or the two real equa-
IPARAMETER - ESTIMATION
rn r-r-.,
I
I s I
I
L_
I
1 --'
I
f]
E +
IT·;~o.l~ If = w~t!;-_-:__
H1 - ...ZM H'
..::--~----
+ wC H'
·-· __
=
~•
.
S'
_____.::._
A· X • E' • S'
..
A
X
CONTAINS THE MEASURED FREQUENCY
CONTAINS THE UNKNOWN PARAMETERS
RESPONSE DATA
M, D, K
-
Hlwl
CRITERION s'-MINIMAL
NORMAL EQUATIONS
The range for the rotational speed is from 0 to 6000 rpm and
for the axial fluid velocity from 0 to 14 m/sec. With this we
achieve Reynolds numbers up to 15000, when the fluid temperature
. 0 .
~s about 30 Cels~us.
.,
~
" Ill
Jl , ••
r ,.•q u•nt (Ht] Sl '"
2000 -PREDICTED
k, ,k ..
(N/m)
PREDICTED C, ,Cu
(Nslm)
-----·
0 0
-----o -----
0 n 0 R n n
5 0 0 0 0 0""' 10 v-""'-
2·10 1000
MEASURED
0 0
I
0 0 2000 4000
ROTATIONAL SPEED n I RPMJ ROTATIONAL SPEED n (RPM)
/
c,. ,c.,
,,
k,. ik., PREDICTED
/
IN/m)
(Nslm)
2·10
5
," 500
,"
,"' MEASURED
l,.._...,_o
~/ " ""'v u I" - u
0 0~--r--+---.---r-~--,
0 2000 4obo 60)0 0 2000 -4000
MEASURED
m,. m,.
(kg 1
10 -1----+---'-----jr-----1
0-6--r----+-...--+--.----1
0 2000 6000
ROTATIONAL SPEED n (RPM)
small and the cross coupled terms much more than this. Meanwhile
it is known, that the reason for the small measured cross coupled
terms is the low entry swirl in the test rig, which has a large
influence on k and k
yz zy
References
STABILITY OF ROTORS
R.NordmaDD
Univenity of Kailenlautem, Kailenlautem FRG
ABSTRACT
GAS COMPRESSOR
TURBINE
POWER 40 MW
OPERATING SPEED 4485 RPM
SHAFT MASS 16000 KG
JOURNAL CLEARANCE JOURNAL SHAFT LENGTH 4.4 m
BEARING EXCITATION BEARING JOURNAL BEARINGS: TWO LOBE d/0"' 0,8
The rigid rotor- oil film bearing model has only a few de-
grees of freedom and is relatively ~asy to handle. However, it
is a good approximation for many practical cases. It is also
useful to study the dynamic behavior of different bearing types.
Mechanical Model
200
3
>-
LJ
s
:z
....,
::J
8100
a::
Q I t I]OJ
u..
s
-' EIGENVECTOR
~
1000 ?000 ~000 400'1 5000
ROTATIONAL SP£EO RPH
Fig. 2 First Natural Frequen-
10
1 cy and Damping Constant
·r
~ o~---,----~---.---.~r-r- of the Gasturbine Com-
sooo
pressor Shaft
OPERATIONAL
SPEED
\:> -20 4485 RPM
:z
a:
~ ·30
STIFFNESS COEFFICIENTS
CENTER
OF GRAVITY
COOROINA TES
RIGID ROTOR WITH
MASS DISPLACEMENTS
Q ROTATIONAL SPEED
PRESSURE P
SI •
~ 12
801)11- SI I t, a, t...a I
llr•R-r CLEARANCE
s, = = s,(£, a,£, a)
( 1)
£ =0 a =o , £ 10 , a 10
s, = 0
= So (2)
ML
a
Fig. 6 Static Equilibrium Posi-
E I tions for Journal
E
Fig. 7 Dynamic Forces due to
Journal Vibrations
aFiL
kik = - - - Stiffness coefficient of oil film
aukL (i,k = 1 ,2)
{4)
- aFiL
cik = • damping coefficient of oil film
aukL (i,k = 1 ,2)
6r 6r It
=k -- c --- ( 5)
ik Fstat ik Fstat
For a given bearing type yik' Sik are dependent only on the Som-
merfeld number or on the static equilibrium position of the jour-
nal, respectively. Fig. 8 points out the dependence of this di-
mensionless coefficients on the Sommerfeld number for a circular
cross section bearing. The values are measured quantities, deter-
mined by Glienicke /1/. We recognize the anisotropic and nonsym-
metric behavior.
182 R. Nordmann
r:t · k F-~-F-1-:=P-jb~
+., .. '.
. .~ ~~:~~
1-'1 - .- ,,., . .
L. - :~~r
100
100 -= - ~ ~£==.:.:.~!!•_l_j
10
10
Equations of Motion
- I
knlu 1 •a'i121 • k12 1Uz-a'i111 •
u
t 11 l u 1 •a~ 2 1 • c12 l 2-aljr,1
k 22 1uz•a'll,l•kz,lu 1 -a'll 2 1•
Czz lil 2 •alil,l • c 2,1il 1 -alil 21 k 22 1u 2 - a'll11 • k21 lu1 •a'll21 •
czzl ilz -alii,) • c2, I u, •aq,21
(6)
(7)
Stability Analysis
EIGENVALUES AND
EIGENVECTORS
\ Ml:i • cu • Ku = 0
U(tl =
K l u = 0
EIGENVALUES
Fig. 12 Solution Procedure for
Eigenvalues and
EIGENVECTORS : Un= 5 +
0 I tn Eigenvectors
a t
u (t) = B e n (S sin(w t+y ) + t cos(w t+y ) (8)
n n n n n n n n
fined by w2
' 0
= g/6r (g acceleration due to gravity). So
0
is the
Sommerfeld number defined with n = wo' which is constant in the
presentation of Figure 14.
I
NATURAL VIBRATION I
D a..t ( . )
u.ltl = u,e s.srniW.t+y.J + t. coslw.t + Y.l
1,0
0
1,0 2,0 w=O lw., 3,0
~ DAMPING
Wo CONSTANT
1,0 3,0
0.
WCiR" 2,66
STABLE UNSTABLE
-1,0
>-
...
u
z
~
a
~
-2,0 ..!!-
-'
<(
a:
~
<(
z
V1
V1
1,0 -~
0
~
TRANSLATION
...~
::t:
a
Fig. 15 Eigenvalues presented
- 0,8 - 0,6 - 0,4 - 0,2 0 0,2 in the Complex Plane
DJI1ENSJONLESS DAMPING CONSTANT aJw.
The stability threshold speed WGR = 2,66 was found for the
translatory motion of the rigid rotor, characterized by a special
system parameter So • If the type of bearing is fixed, So 0 is the
0
only parameter defining the stability onset speed. For other pa-
rameters So 0 other limit speeds WGR have to be calculated. If
another bearing type is selected, the threshold speeds will also
change. Fig. 16 shows a stability diagram for the rigid rotor with
oil film bearings, pointing out the rotational speed limit WGR in
dependence of the Sommerfeld number So o • The different curves be-
long to different bearing types. Such diagrams are very useful in
the design process for rotating machinery.
i ~ :' ! ; l
~~ 1 :
1 ! I l l'
~.; -- · -- . s _ P·or2
.. Oo - ~ Wo
Wo = ~
_! :!''
.:;,·
stab le
-: J
! :
.,
0. 1 10
tor and two cylindrical journal bearings with circular cross sec-
tions. A d.c. electric motor with speed control drives the shaft
in a speed range up to 6000 rpm.
~
JOURNAL
BEARING
Fig. 17 Test Rig with Rigid Rotor and Oil Film Bearings
!
a,
i
w.
0.5 2.5
- - MEASUREMENTS
- - - CALCU.ATION
References
ABSTRACT
A ::.acKage of programs which can be used to study the static
and dynamic behaviour of a multi-supported shaft line is
illustrated. The analysis is carried out by Keeping account
of ~he interaction with a supporting structure of same
(foundatlon, casings, etc.>. Together with the structure of
the pacKage, this paper is aimed at illustrating the
~athema.tical mocel wsec ~o~ static and dynamic analysis.
7. 1. 1 Foreword
"P.A.L.L.A." is the name of a pacl<age of programs developed by a research
group worl<ing at the Sezione di Meccanica dei Sistemi of the Dipartimento
di Meccanica of the Politecnico of Milan. This package is the result of a
series of research worK, carried out over the last ten years. in
collaboration with ANSALDO, ENEL-CRTN, ENEL- DPT.SASt TECNOMASIO
ITALIANO BROWN BOVERI and 'WORTHINGTON NORD ITALIA. The pacKage
can be used to study the different problems that arise in rotor dynamics:
the analysis is carried out by l<eeping account of th• iteraction between
the shaft line and the supporting structure of same (foundation, c&se,
etc.>.
192 G. Diana et al.
<A.1>
where
-F* = {.fRJ
F* =
-F
<A.2>
f = {~:} <A.3>
In eq. <A.4> the forces .BF transmitted by the foundation are unl<nown
functions. In order to define these functions, as has been said, a modal
approach wu used. The vector g of the modal variables of the foundation
can be correlated to the displacements of the connection nodes .!F of the
foundation itself by means of the following equation:
~F=£~l9, <A.5)
194 G. Diana et al.
<A.6>
where cm,J, CrFJ and [I<FJ are diagonal matrices containing the modal
parameters of the supporting structure. If as many modes as degrees of
freedom ASsociated to the connection nodes are assumed, CcpJ is a square
matrix and can be inverted thus allowing the transformation
-1
~ = [cf l ~F <A.7>
<A.S>
where prime indicates the transposed matrix. Equation tA.S> t•> correlates
the displacements of the connection nodes ~F to the forces .BF transmitted
to the rotor through these nodes. By now substituting the expression of
the forces BF of (A.S> in equation tA.4> and by moving the terms containing
j,, iF and .!F to the left hand side, the fundamental equation <A.t> can be
rewritten as:
•• •
- -
[Hl X + [Rl X + [Kl X = F <A.9>
where CMJ, CRJ and CKJ repres!!nt the mass, damping and stiffness
matrices of the overall system that keep account not only of the
t•> If the eigenvectors are orthogonal in the rutrict sense, eq. A.7> can
I
also be written as g • ctJ !F·
P.A.L.L.A. :A Package to Analyze ... 195
inertia., damping a.nd stiffness of the rotor but also of the bearing/seals
fluid films and the effects of the carrying structure. Vector ..E in eq. <A.9)
therefore represents the vector of the single external forces applied to
the nodes of the rotor and to the connection nodes. Equation <A.9>
therefore permits, once E is Known a.nd different from zero, any problem of
forced motion to be solved. Always with eq. <A.9> on the other hand, by
assuming E= 0 it is possible to study the free motion a.nd to analyse
instability condition by solving the associated homogenous problem.
Methods of solution
As will be e:<pla.ined further on, the different problems that are
encountered when studying rotor dynamics <forced vibrations and
instability> are cast, by the P.A.L.L.A. pacKage programs, into the form of
the evaluation of the frequency response to a harmonic excitation force of
the entire rotor-fluid films-foundation system, thus adopting a more
stable algo_rithm of solution.
As has been said, all the problems are studied by analysing the frequency
response to a. harmonic e>:cita.tion force of the type:
Both here and further on, we will indicate the excitation frequency by Oe:
and the rotating speed of the rotor by .DR' both expressed in r/s <radians
per second>. With this excitation force, the forces at the connection nodes
are also of the type :
<A.12)
<A .13)
In the case of the harmonic excitation force considered, by, using the modal
approach it is easy to define the matrix of the mechanical impedances
CH.Ile;>J. In fact if we assume that a. unitary harmonic excitation force is
196 G. Diana et al.
(A.15)
(-A-2[M*l
Mt:
+ i''-tR*J
--E
+ lK*l> -o
X =-o
F * <A .16>
...
, .. ·: ·. Fo{)ndation
... ..
~
7 .. 1. 3 Package Configuration
The P.A.L.L.A. pacKage is basically interactive both as regards the data.
input as well as the choice of the type of processing to be used.
Calculation can be carried out either using the "demand" and/or "batch"
modes. The handling' of the set of the different programs is carried out by
means of a "PacKage General Menu" which carries out the various functions
by transferring control to different main programs which in turn recall
specific subroutines. Data. transfer taKes place by means of
mterconnection logical units which support the different input and/or
output files of the various programs. As specified further on, some of
these files are created interactively by the user, others automatically a.s
the output of subprograms. The preparation and use of the above
mentioned files allows the calculation to be carried out even in the "batch"
mode. The logic chart of the pacKage that evidences the different functions
is shown in fig. B.t. with the numbers of the relative "Functional BlocKs".
BlocK 0 is the "PacKage General Menu" for the handling of the various
programs: this guides the user in the choice of the different functions
carried out in the other bloc~. BlocKs 1,2,3 and 4 are recalled for the
introduction, in colloquial form, of the data relative to the shaft line,
bearings, seals and the carrying structure, which can be schematized by
means of the mechanical impedances or by means of uncoupled
mass-spring-damping systems. These same blocl<s supply the plots of the
shaft, bearings and carrying structure in order to checl< the accuracy of the
data introduced: eventual correction are carried out interactively. Once the
accuracy of the data introduced has been verified, data. files corresponding
to the different components previously mentioned are created
automatically from BlocKs 1,2, 3 and 4. BlocKs S to 13 carry out
intermediate calculation functions and are activated by means of
sequences piloted by the "General Menu" in function of the final
calculation which the user intends to carry out. BlocK s, which calculates
the overall mass and stiffness matrices of the shaft line <without fluid
films>, is carried out automatically, in that it is an obligatory step,
whatever the final type of calculation expected. The same goes for BlocK 6
which carries out the static analysis of the shaft line, analysis that
always have to be done since it supplies the loads on supports once
200 G. Diana et al.
o·
1--
i: I s~~a-r ~-11 nt
~
t
H
I
I: I ;fa!":~;!
- )i~• fi 1 t! :rti~=~n
I,I :n:: H
- ~cm~c~•~~ ;~~~~~= c·sc·a~
[ on;~-~~::-:~~:~~~~~~~-~~~~~;,---]
~ull ana.' or s:< 1e
~eou~!d l
+
~hysr~a.~ I"''Odt~ 1
---------------1-------------
• !
[~~~~::~~~:~~~~:;~~~~:~~~~~~~~~~~~~~~~~~]---~-[~~~~~~~~~~~~]
~ ShaH-1 i nt mass and st: dn!H
atr1cts
6 St~ttc anal/Sts:
al without ilu•d f1lms
OJ w1th "lu1d -films
Gr a~n r c
!3 Non~synchronous txc•tat•on
- Journai ovei ::i.ttons
- Crac-< tf~t':t!
alignment has been assigned, or the alignment, once the loads and the
pc•sition of some of the supports have been assigned. In this same BlocK it
is also possible to Kee~ ;::our;t of the - -~--.;;c;;, :; "':he lubricating film to
define the loads on s:. · · : ·"':s once the alignment has be~r- =::": :, ' :· :·r
vice-versa. For this reason it is necessary to Know the fluid film stiffness
of bearings or seals, characteristics which are calculated at a given
angular speed of the shaft line O.R in BlocKs 7 and 8, by assigning the first
attempt loads defined in BlocK 6 as input loads. To improve the calculation
an iterative procedure between BlocK 6 and BlocKs 7 and E: is followed.
BlocK 6, therefore, calculates the alignment conditions while the shaft is
static•nary <in the absence, therefore, of an oil film, or by Keeping account
of the oil film stiffness corresponding to an assigned angular speed). Once
the alignment conditions have been defined, BlocK 14, on request, gives the
stress calculation along the entire shaft line due both tO' the alignment
conditions as well as, for the horizontal shafts, to its own weight. BlocK 9
1s activated should the user wish to Keep account of the carrying structure
through the mechanical impedances. In this BlocK, the dynamic behaviour of
the carrying structure itself whose transfer functions are defined in
correspondence to the connection nodes with the shaft is simulated
through a Finite E:lement schematization. These transfer functions can
also be assigned by the user, as in the case where they are defined
starting from the e:<perimental data. BlocK 10 utilizes the transfer
functions derived from the calculation <BlocK 9), or assigned by the user to
define the modal parameters of the carrying structure. The matrix of the
mechanical impedances is also calculated from the modal parameters for an
assigned frequency .Oe;. The data processed in BlocKs 5, 6 1 7, 8, and 10 are
used in BlocK 11 to create the elastodynamic rna tri:< CE: <Ae;) J of the whole
system, matri>: defined by relation:
( B ,1)
This matri:{, besides the mass and stiffness matri:< of the shaft, contains
the rnatri:< of the mechanical impedances of the foundation and the
stiffness and damping matrices of the oil film which are a function of the
angular speed of the rotor ~ which coincides with the frequency of the
202 G. Diana et al.
- "Crank" effect
- transverse crack
BlocK 16, using matrix [E<flg>J defined in BlocK 11 and the stiffness and
damping defined in Blocl<s 7 and S carries out, at the user's request, the
stability analysis of the shaft line. Graphic outputs are also available for
Blocl<s 14 and 15. Bloc!< 3 relative to the introduction and creation of data
files for the seals is still in the implementation stage. So is Bloc!< 8
relative to the calculation of the equivalent stiffness and damping
coefficients of same.
The part of the program described in this BlocK deals with the introduction
of the geometric data of the bearings of the shaft line being examined. The
thr-ee types of bearings presently considered by the program and which can
therefore be chosen by the user are:
- multi-lobe
204 G. Diana et al.
Multi-lobe bearings
With reference to the generic i-th lobe <see fig. 2.1), the geometr1c
parameters to be assigned are :
- xci, Yci: . coordinates of the center of the lobe i-nth in a reference
system with origin in the theonc center of the !Jea.rmg ana
oriented as in fig. 2.1;
- eoi : ratio between the clearance of the lobe i-nth <assumed
circular) and the reference clearance;
- 8'oi.: start angle, measured a.s in fig. 2.1, of the 1-th lobe;
-17'-~~ : angular extension of the i-th lobe;
- JD<i> : number of cylindrical parts into wich the i-th lobe 1s a>:ially
divided <max. 31.
Assuming that the JOurnal axis is displaced by xc1,vc1 with respect to the
center of the i-th lobe. the film thicKness is defined by a function h!~ <see
[5 J) which is dependent on the bearing characteristics. When the journal
center assumes the coordinates x ,v in the reference system of fig. 2.1, the
film thicKness htm assumes the value given by:
In the case of the circular lobe considered uo till now we have equatlon
P.A.L.L.A. :A Package to Analyze ... 205
<2.3)
which specif1es the meaning of carameter eoi• For the cases in which the
bearing profile is not circular, function hEat assumes the expression:
<2.4)
Therefore, only the widths of the single cylindrical portion of the lower
lobe <b 1> and that of the upper one <b2> <fig. 2.4) have to be assigned. The
cylindrical bearing without pockets or grooves shDuld be treated as a
one-lobe bearing with a 360 deg. extension.
206 G. Diana et al.
<2.5)
in which r 0 and cp 0 are the polar coordinates of the hinge and 1\.jJ is the
generic small rotation of the pad from the position assumed to evaluate
h0 i<-fh. The quantities Xci and Yci determine what in literature is normally
indicated by preload, which in the program is assumed to be a. radial
displacement of the tilting pad center with respect to the theoretical
center of the bearing (generally the center of the circle passing through
pad hinges>. The geometric definition of a. tilting pad bearing besides the
quantities r 0 and cp 0 which appear in eq. <2.5> and the radial clearance
already defined at the beginning, (it is important to note that different
radial clearances are not expected for the var1ous pads and is therefore
always e 0 i = 1> requires, for each of the tilting pads, the following data :
- the preload of the tilting pad, defined as has already been said.
Non-dimensional ratio preload/radial clearance is used;
- the moment of inertia J 0 of the tilting pad with respect to the hinge,
assigned by means of an equivalent mass me= Jp/r 0 2,
- -\lctstart angle of the i-th tilting pad measured as in fig. 2.5;
-~~angular amplitude of the i-th tilting pad.
Moreover, the axial length b of the tilting pads must be assigned.
Roller bearings
For these type of bearings the stiffness matrix, i.e. the coefficients K>:>:t
Kxy' Kyx' Kyy in the previously defined reference system, must be directly
assigned. The introduction of the geometric data. of the bearings taKes
place bearing by bearing, with the possibility of controlllng the data
introduced by means of graphic outputs.
Mathematical Model
The carrying structur·e is schematized via F .E .M. with 12 d.o.f. "beam"
elements ([6J, [7J,l. These elements (fig. 4.1l represent homogenous beams
w1th a constant section, undergoing a>:ial action, shear, twisting moment
and bending moment along the two a>:is of inertia of the transverse
section. In Tab. 4.1 the corresponding stiffness [~J and mass [mtJ
matrices evaluated in the local reference system are shown.
The corresponding structural matrices EKe J and CMt'J calculated in the
absolute reference assumed for the overall system are calculated by means
of the well Known relation:
.
where CIGJ is the matri>: of the polar inertia moments of the rigid body
wit:1 respect to the three principal axis:
0 1 OJ0
IyG (4.3)
0 2G
208 G. Diana et al.
Once all the structural matrices of the elements that comoose the complete
model of the carrying structure have been defmed, it is coss1ble to
assemble the total stiffness matri>: [}{sJ and the mass matri:{ [MsJ using
the usual techniques. Having assigned the constraint conditions tc the
structure, after partitioning the equations of motion relative to the
constrained and unconstrained ([6J, [7J> d.o.f., the equations of motion of
the entire carrying structure therefore becomes:
(4.6)
. 4. 7.'
whose solution:
X =X eiwt
-s -s
enables us to obtain the natural oulsations of the svsterr. and relative
principal vibrating modes<*). On the other hand by assuming
r~,..
-- -so l''.D.t
t •=F ~4.9)
the supports alignment on starting up the machine. For the rotors with
horizontal axis an alignment is normally imposed so that the loads on
the bearings correspond to the so called "single beam reactions", i.e.
reactions exerted by the uncoupled rotors. Each rotor is normallv
placed on two bearings: the constraint reactions determined in this wav
are also maintained when, by means of the end joints, the
aforementioned rotor is constrained, flexurallv and torsionally, to the
other rotors composing the shaft line.
b) determination of the loads once the alignment has been assi~ned. This
method is useful when we wish to evaluate the effects of the support
displacements: e.g. thermal distortions of the foundation structure can
cause, on one or mar• supports, load variations with respect to the
nominal with consequent phenomena such as oil film instabilitv.
.
As previously indicated, the shaft line is schematized with 8 d.o.f. "beam"
finite elements. To obtain the stiffness and mass matrices of the rotor-oil
film system it is necessary to suitably assemble the stiffness matrices
cln <alrudy described in Blocl< 5> of the single elements and the
equivalent oil film stiffness matrices <see Bloc!< 7>. The d.o.f. number of
the shaft nodes is gen~rally very large so that condensation becomes
necessary (fig. 6.1). In the static analysis carried out whether the
lubricating film is l<ept account of or not, the condensation adopted, of the
static type, proves to be rigorous. The use of this approach decomposes
the generic shaft line into a set of substructures <called "superelements">
each composed of several beam elements with connection nodes between
superelement and superelement as extremitY end nodes <called
"supernodes">. As a substructure the program assumes the portions of the
rotor connecting two supernodes. Extremity nodes and bearing nodes are
always assumed to be supernodes <see fig. 6.2>. For the generic
superelement the equation of static equilibrium can be expressed as:
( 6 .1)
~e -_{~ee}
. ~e 1
(6.2)
thus obtaining:
(6.3)
or rather:
(Keel~ee+[Ke;l~e;=fee a)
(6.4)
£K,el;:$ee+[K; ;l!e;=fe, b)
(6,5)
(6.6)
(6.7)
and ·with
(6.8)
[KcJX
-ee =P
-c
In this way it is possible, to solve the static equation m rigorous form. b'!
using as degrees of freedom those relative to the single suoernodes.
Having calculated the matrices of the single superelements and the
matrices of the equivalent stiffness of the 011 film <BlocKs 7 and :?.) in th1s
way, it is possible to proceed with the assembly of the complete rotor + oil
film system as shown in figure 6.4. As regards the only shaft line, the
equations system for the static situation considered is therefore given bv:
[Kl~=£ (6.10)
From system 6.11> we can formally obtain the disolacements of the free
nodes !Land the constraint reactions .fv according to the system:
~L=[Klll-lfL-[Klll-l[Klvl~v (6.12)
fu=CKuLl~L+[Kuvl~v (6.13)
by solving system 6.12> first and by substituting the values found thus for
!L in eq. 6.13>. Coordinates !v are always, as is obvious, bearing node
coordinates. They are as has been said, in a minimum number of two per
plane, horizontal and vertical, and in a maximum number equal to twice the
number of bearings <at the bearing nodes rotations are never constrained>.
The problem, discussed in point b) of the previous paragraph, of
P.A.L.L.A. :A Package to Analyze ... 213
Characteristics
The part of the program relative to this BlocK defines the force field that
is set up between the journal and the bearing due to the lubricating film,
as a function of the relative journal-bearing position. This field is
linearized in the neighbourhood of the static equilibrium position, thus
defining the equivalent stiffness anp damping coefficients. The calculation
of the equilibrium position of the journal inside the relative bearing and
the corresponding values of the stiffness and damping constants of the oil
film is obtained by integrating Reynolds' equation for both a laminar and
turbulent regime, Keeping account of the temperature variability along the
film. Reynolds' equation is written, in cylindrical coordinates (fig. 7.1l in
the form:
7 .1)
having indicated as
Fx= Fx~x.y,x,y,.!lR)
(7.2)
Fy= FyO:x ,Y,x ,Y, .Q.R)
(7.3)
m which U represents the peripheral speed of the journal <.a.tR> and V the
Kinematic viscositv. From equation (7 .4>, having assigned a temperature T0
at the beginning of the meatust solution T~ and the local viscosity value
)J e> from eq. (7 .3> are obtained. In the c:ase of laminar flow regime the
visc:ositv is defined as:
( 7 .5)
216 G. Diana et al.
On the other hand, in the case of turbulent flow reg1me. tne orograr;.
extends the validity of Revnolds' ecuation bv mtroducing the e>:cress1ons
of a turbulent viscosity derived 1n [12J on the oas1s of ~acers [13J ana
[14J. According to the findings snown. 11: is assumed tnat there are t·,..,c,
different e>:pressions of the turbulent viscosity, i.e.:
0.6.1
in the axid direction. In [14J the values of C'O'a.nd C'1 a.s a. function of the
local Reynolds' number Rh a.re not directly given. For use mside the
program interpolating functions with expressions
were determined which g1ve results that are fa1rly close to those of the
paper cuoted above for values of Rh ranging from 900 up to 200000. Tne
experimental results obtained in a test campa1gn, earned out on elhpt1.:
bearings normally used real turbosets showed the necesslt'! to mod1fy
1n
these expressions of the turbuler.t viscos1ty. For this curpos: ·.:ee ~15J)
other modifications were also introduced to a.na.lvt1callv reproduce the
distribution of the oil film temperature measured e:·:oerimental!v. As
shown in [16J,[17J,[18J an e}:cellent agreement between the e:<perimenta.l
locus of journal center and the a.na.lvt1ca.l one was obta.med bv usmg,
instead of eq. 7 .8>, expressions:
dh
( de-
)2 dh
d.. =a + b + c (7.11)
h
Q'6' = u (7 .12)
2 12f R
Wd = 2~ n2u 2R
.KE 0 XX
s i n <D
r (7.13)
where Rxx represents the direct dissipation term of the oil film and ~ is a
corrective term. By dividing this expression by the carrying surface of
bearing rtR b, we obtain an averaged power per unit area.. Moreover, this
power wa.s "weighed" with the local pressure p<-lt> compared with the
a.vera.ge pressure Pmt thus obtaining expression:
218 G. Diana et al.
= b R rc
-- =
Pm
7 .14>
R
7 .15)
Q~
where Q,. should be evaluated through equation <7.12) a.nd Rxx can be given
the value rxx calculated by the program.
Fy<x,y,O,O,Qi = Qy
7 .16)
Fx<x,y,o,o,~ = Qx
are satisfied and where Qx and Gy are the load components on the journal,
according to the vertical x and horizontal y directions, evaluated as
previously described in Blocl< 9.
where
the dAmping constAnts oof the lubricAting -film. The numericAl values of the
derivAtives in expression 7 .18) and 7 .19) are calculated by assuming thAt
the functions F x and F y have a. parabolic trend in the neighbourhood oof the
equilibrium position and are defined by separately giving increments 8x,
8y, 8x and ~y.
------ --------- - -
·-Y
'
y
JIJ(2}=3
~
......
- --
~
In ~..J
...
2 3 4 5 6 7 8 9 10 11 12
1/3
13/35+6J1/~1 2
.13/35+6JII'/~J2 SYtt1ETRIC
J..l~
-111/210-Jr/10AI JCZ/105 +2Jv/15'A <t h11tn tr not shDIIIII art ztro)
111/210+ Jz/10AJ )2/105 +2Jz/15'A
(t1t]•f'l· l/6 -<12_~_2) <1,1>
9fi0-6Jz/5'AJl! -(11,3) (2,2>
9/70-6Jy/5'A12 (3.3)
· - - - - -f- ··- -
J~/6A _J_~~l.
13J/420-Jv/10AJ -)1 /140-Jlt'/30A
-- r - - - - - - - -
-(5.3) (5.5) -,
-131/420+Jz/lOA1 -1 2/14D-Jz/30A -(6,2) I (,,,,
TAB. 4.1 - 12 d.o.f. beam tlemtnt stiffntss and mass matricts t-)
t-)
.....
222 G. Diana et al.
····a...j
!
0 0 0 Rz -Ry
0 0 -Rz 0 Rx
0 0 Ry -Rx 0
m*
0 MRz Ry T 1l ( R/+ R/ )/m* T 121m•- RxRy T 13/m•- RxRz
I I
I 12EJz112 I I
1----------1---------- I
I I 12EJy112 I SYtt1ETRIC I
1----------1----------1-------- I
I I -6EJy/12 I 4Eoly11 I <eltmtnts not shown art ztro> I
1----------1----------1--------1--------- I
I 6EJz112 I I I 4EJz/1 I I
1----------1----------1--------1---------1-------- I
I -<1,1> I I I -<8,1> I <1,1> I I
1----------1----------1--------1---------1--------1-------- I
I I -<2,2> I -<7 1 2) I I I <2,2> I I
1----------1----------1--------1---------1--------1--------1------- I
I I (3,2) I 2EJy/1 I I I -<3,2> I <3,3> I I
l----------l----------l--------l---------1--------l--------l-------l-------l
I <4,1> I I I 2EJz11 I -<B,I> I I I <4,4) I
I I I I I I 1___ 1_ _ _ 1
I I
I 1:56 I I
1------1------ SYit1ETRI C I
I 221 I 412 I I
1------1------1------ (t1tmtnts not shown art zero> I
I I I 1:56 I I
1------1------1------1------ I
I I 1-221 I 412 I I
[ml= 1------1------1------1------1------
I :54 I 131 I I I 1:56 I
I tltl<D,2-o.2>
I( --------!___ + 11
) I
t420
1------1------1------1------1------1------ I 4
I -131 I -312 I I 1-221 I 412 I I
1------1------1------1------1------1------1------ I
I I I :54 I -131 I I I 1:56 I I
l------l------l------1------l------l------l------l------l
I I I 131 I -312 I I I 221 I 412 I
1_ _ 1_ _ 1_ _ 1_ _ 1_ _ 1_ _ 1_ _ 1_ _ 1
~ t :0.. ~ Supornades
------.---
-n-~ ~~ ~n } , { i)
Supernodes
4 IS
10
II)
~
n
Fig. 6.2 - Subdivision into superelements "
~
P.A.L.L.A. :A Package to Analyze ... 225
4 t--
element
~
2
.3 .I matrix
4 ~~
s ~ ~-
6
7 f234S 618910
8 l
9
10
[l<e] .__ w
FiQ. 6.8 - SupeNlement matrix
4 2 3 4 5 b 7 8 9 10 11 12 13
Fig. 6.4 - Overall stiffness matrix of the system <on rigid .foundation)
226 G. Diana et al.
~X
I
tx
I
<9.1)
where EMsJ, ERsJ and EKsJ are respectively the matrices of mass, damping
and stiffness of the foundation. In equation <9.1> .!s represents the nod~t
displacements of the same supporting structure, while.! 5 is the vector of
the applied excitation forces.
P.A.L.L.A. : A Package to Analyze ... 227
F =F ei~t (9.2)
-s -SO
(9.3)
In the frequency domain eq. (9.1), Keeping account of eqs. (9.2> and (9.3), can
be rewritten as:
(9.4)
(9.5)
~soc=[H<Ae>lfsoc (9.6)
fhis matri>: is evaluated for a generic pulsation ~' starting with the
modal parameters of the foundation, parameters that can be evaluated if
only one column of CH<Ae;lJ is Known, evaluated for a certain number of
pulsations De;.
The generic 1-th column of the flexibility matrix CH<O..e;>J
represents the <comple>:) vibration amplitudes of the d.o.f. of all the
connection nodes for a unitary harmonic e>:citation force applied to the
1-th d.o.f., therefore corresponding to the index of the column considered.
The generic term hKl of the column 1-th is, as has been said, complex and is
called a transfer function between the K-th d.o.f. and the 1-th d.o.f. of the
rotor-foundation ccmnection nodes. To obtain the 1-th column of the matrix
CH<Oe;>J it is enough to therefore impose a unitary harmonic excitation
228 G. Diana et al.
(10.1>
where [Ml CRl And [Kl respectively represent mass, damping and
stiffness matrices of the system, ! is the vector relative to "n" d.o.f. and
l is the vector of the excitation forces. The modal matrix r+J, generally
rectangulAr n x p <where pis the number of the vibrAtion modes considered)
is formed by p eigenvectors f (j), ordered in columns Cthese values are
obtained as eigenvectors of the matrix [M ]-1 [KJ and for a stable
conservative system are real), i.tt.:
230 G. Diana et al.
(10.3)
X (j)
"
By means of the coordinate transformation:
(10.4)
(10.5)
where r-m-J, t"-r...J and L-k-J are the mass, damping and stiffness matrices
in principal coordinates (square and diagonal matrices pxp> and g are the
principal coordinates themselves Cthe matrix t-r-J is diagonal if [RJ is of
the type: [RJ=a[MJ+b[KJ). The generic diagonal element mj of eq. U0.5)
represents the generalized mass of the j-th vibration mode r j is the
associated damping and kj the relative stiffness. If the vector of the
excitation forces! of C10.1> is harmonic:
(10.6)
(10.7>
(10.8)
(10.9)
or rather
(10.10)
Once the modal parameters mJ r j' and l<j and the corresponding vibration
modes have been defined1 with the algorithms that will be described
hereinafter and having applied the transformation <10.4) we obtain:
(10.11)
i.e. :
(10.12)
(10.13)
As can be seen, once the modal parameters are Known and once the
pulsation.O.e: of the excitation force has been assigned, the matrix CH<Ag:lJ
is numerically determined and therefore, from its inversion it is possible
to obtain <see first chapter) the matrix of the mechanical" impedances
CI<Ot!:lJ. The generic term h1<1<Ae:l of the matrix CH<Og:>J is also named as,
as has already been mentioned, transfer function and represents the
complex frequency response of the system evaluated at point K, having
applied a unitary harmonic excitation force of pulsation~ at the point 1.
From <10.13) the generic transfer function can therefore be expressed in
function of the only modal parameters as:
(10.14)
where
w-2
J = Kjm·J represents the j-th natural frequency of the system
(j)
XI< represents the generic 1<-th component of the j-th
vibration mode.
232 G. Diana et al.
r · r · r ·
6"· = ~·J· = •· -L = .,. ~ = :...L. <10.15)
J J J J r . J 2m·t~t· 2m·
CJ J J J
Method n.1
If the coupling between the vibration modes is negligible, the transfer
0
function phase presents a variation near to 180 in a frequency range near
to a resonance. A method which enables us to evaluate the damping factors
is based on the analysis of the derivative of the phase evaluated in
correspondence to the natural frequencies of the structure C19J. In fact,
phase "(( of the transfer function in a system with one d.o.f. can be
expressed as follows:
(10.16)
dl\!' 2
- = (10.17)
1 1
2'!2 = w.S
(10.18)
(10.19)
234 G. Diana et al.
Method n.2
This algorithm is based on the calculation of the frequencies wd and wt,
near to each resonance of the structure, in correspondence of which the
amplitudes of the transfer function is 11V2 of the amplitude of the
resonance peaK. The calculation of the frequencies wd and wb is performed
by determining the solutions of equation [2'0J:
(10.20)
assuming St<r Jrc.} << 1 and neglecting higher order terms we obtain:
= 1+2E!:·
- 5J (10.21)
(10.22)
or
(10.23)
The dampin!l factor associated to the j-th vibration mode can be defined by
the imaginary part which from ratio Ay<j)Kl/Ah(j)Kl [2J is:
P.A.L.L.A. : A Package to Analyze ..• 235
AYI< 1 ( j)
bhl<l(j)
= .,d·+i6'·
J J
<10 .25)
where with wdj the j-th pulsation of the damping system is indicated.
(10.26)
where with U tiVj the residue associated Xt<(j) to the j-th vibration mode is
indicated. The terms Re<D.g> and Im <De:> represent respectively the real
part and the imaginary part of the transfer function hl<l<Oa:> evaluated at
frequency .Og. Equation <10.26) is the equation of a circumference whose
center has the coordinates given by:
Xcj=-U j/26j
(10.27)
y c j=-V j/2 6j
(10.28)
from which the residue value relative to the generic j-th modtt is
proportional to the diameter of the circumference <eq. <10.26». The
accuracy of the results obtained with this method depends directly on the
solving method used to determine the transfer function and on the number
of points used to calculate the interpolating circumference.
where x1<j> represents the amplitude of the j-th vibration mode evaluated
at the point 1, Flo is the 1-th term of vector f 0 and qjo the j-th term of
vector So of eq. <10.9). By assuming resonance condition <Qg=CAJ} we will
have:
jn-r·q·
'"'"!: J JO
=i~·r·q·
J J JO
=X1(j)F1 0 (10.30)
2mj~j
a)
2mj(A)j
<10.31)
2imj~j~j2qjo = Xl(j)Flo b)
from which:
xl<j>Flo
(10.32)
2i~jf.J/Qjo
P.A.L.L.A. : A Package to Analyze ... 237
In eq. <10.32> Flo is Known, ~ j' x1(j) can be determined by the methods
wj,
described in the previous paragraphs. For the evaluation of qjo the
relation:
(10.33)
(10.34)
or
n
q JO ~ XK ( j >xk o
. --~ (10.35)
K=1
where XKo is the generic K-th element of the vector 1o' evaluated
experimentally by applying to the system an excitation force
F = F 0 exp<:iDe;t> or
by analytically calculating the latter in the same way
<see BlocKS>. By Keeping account of equation <10.34) eq. <10.32> becomes:
(10 .36)
(10.37)
(10.38)
Using this hypothesis, the error between the given transfer function
<10 .39)
The sum of the quadratic errors evaluated for the transfer function
determined in N measurement points is given by:
n 1 )
Et<6tj > =~ E2 ('>j, -. <10 .40)
k =1 . mJ
It is possible to minimize this error by imposing the following condition:
'c) ( 1/m j)
=0 (10.41)
(10.42)
Method n.2
The second method adopted tries, for each 1<-th d.o.f., to simultaneously
calculate the generalized masses mj of all the vibration modes: the user is
thus able to accept, mediate or discard the generalized masses calculated
for the generic j-th vibration mode in correspondence to the different
degrees of freedom. To identify the generalized masses the method
minimizes. for all the Wj considered in correspondence to the generic 1<-th
d.o.f., the function hl<l<'-»} (j=1,2, ..... p) defined by equation <10.14), with the
same known functions <determined experimentally or analytically> hl<l*<w;.
Analytically this can be expressed as:
10.43)
P.A.L.L.A. : A Package to Analyze ... 239
or
In the case in which the number of the sampling frequencies is higher than
the gen&ralized masses to be determined, it is possible to define the
following error function:
1
g = h*-h = h*-£Tl 10.45)
m
To minimize the quadratic error:
1 T 1
gTg = [H-[Tl - ] [ tf-[Tl - ] 10.46)
m m
it is necessary to impose th& condition:
10.47)
thus reaching:
(10.49)
The more accurate the evaluation of the single modal parameters, the more
the function hK1<tJ} approximates the starting function hKl*<Wj • The
comparison between the starting transfer function and the one obtained by
means of the synthesis operation therefore enables us to evaluate the
reliability of the modal parameters calculated. The modal parameters w.i' G' j
and mj of the carrying structure (natural frequencies, damping factors and
generalized masses) are memorized, together with the vibration modes (j),f
on a file to subsequently be used to calculate the matrix of the mechanical
impedances. As far as the vibration modes are concerned, only the mode
components corresponding to the connection points <i.e. between the
carrying structure and the shaft line> are memorized. In fact, these
components are the only ones that are actively involved in the definition
of the matrix of the mechanical impedances of the foundation. Nonetheless,
as has been mentioned, to better define the vibration modes of the
carrying structure it is convenient to also calculate the transfer function
<BlocK 9) in correspondence to foundation points different from the
connection nodes. It is important to underline that, in order to define the
matri>: of the mechanical impedances CI <.Oe;>J it is necessary, for each
single modal analysis (i.e. for each direction along which the e>:citation
force acts) to identify a number of natural frequencies at least equal to
the number of supports. Moreover, the natural vibration modes identified
must 'have significant components in correspondence to the connection
nodes between the carrying structure and the shaft line. Therefore, all the
vibration modes in which a considerable number of connection nodes
coincide with a vibration "node" (i.e. a non vibrating point) have to be
discarded. This fact might be verified e.g. for torsional vibration modes of
a foundation slab. (fig. 10.1>. The modal parameters are determined by
analysing the transfer function hK1<0e;> (with K=1 12 1... n> evaluated in n
points of the carrying structure whith a. unitary harmonic excitation force
applied at point 1. As already been said, the transfer functions can be
expressed as a function of the modal parameters only (see equation (10.14):
P.A.L.L.A. :A Package to Analyze ... 241
undeformed slab
_.,,
- ,..,. ' ..
'
'
'
- -·
;""";.,
''
,--.
''
-·· - ·-·
I
: -- .. :
~-.
: ' ' .. I ,.
___ :. --
___ - ' .,..-· -""'--
'' ' I _ ......
'
supports axis line ' ' ',, - deformation shape
The d.o.f. nu.mber associated with the shaft nodes, as already stated for
that which concerns the static analysis <BlocK 6>, is generally large and a
dynamic condensation becomes necessary. In the dynamic analysis carried
out in the frequency domain used to evaluate the response of the shaft
line to a harmonic e:<citation, this condensation is rigorous. Hereinafter
we will firstly describe the method adopted to reduce the shaft d.o.f. by
means of dynamic condensation and subsequently will show the assembly of
the condensed matrices of the shaft with the matrices of the mechanical
impedances of the foundation as well as the stiffness and damping ones of
the oil film and/or seals.
Dynamic condensation
We will now analyse, for the saKe of convenience, the shaft line shown in
fig. 11.1. Let us consider a set of beam elements corresponding to a
portion of the rotor which will be considered as a superelement. By
indicating with ~e the relative d.o.f. vector, the equations of motion of
this system can be written as:
(11.1)
where [Me J, CRe J and CKe J are respectively the mass, damping and
stiffness matrix of the generic superelement anal)tser. For the example
proposed, these matrices are asser.:tle.~ ,:..s shown in fig. 11.1. In eq. <11.1>,
.f e
represents the vector of the external excitation forces applied to the
nodes. In the P.A.L.L.A. pacKage the steady state response of the rotor to
harmonic excitation forces is analysed, so that we have:
F
-e = !-l?O
F eif:lEt (11.2)
X
~
= -eo
X ei.O.Et (11.3)
(11.4)
P.A.L.L.A. :A Package to Analyze ... 243
rSuperel emel)t I
7/ ~8
fL.= £r £~
f 2 3 4 5 6
:
E-• s~s
Itechdn/ca l impecld/Jces
2 3 4 5 6 7 8 9 S1 S2 S3 S4
1 A
2 A ITJ
3
4
5
6 A [i]
7 LBl = [E<.QE)l
8 A I.....L
9
Sl []] ';'k :: t·;.: ··: [I <.D ) ]·'''
~l}~~·Kttf:;~i
S2 [iJ_
S3 I B
S4 B I
n
Kxx <IIJl> + i.D£Rxx<.DR> 0 Kxy<~> + i.DeRxy<~>
[Al =[ 0 0 0
Kyx<~> + iDERyx<~> 0 Kyy<.OR> + iOERyy<llR>
0 0 0
CBl = -tAl
Having called Xee the vector of the "external" displacements, i.e. the
vector containing the d.o.f. of extremity nodes of the portion analysed
(nodes 1 and 10 of fig. 6.2> and kithe vector relative to the other nodes
("internal" nodes>. it is possible to renumber the equations a.s:
(11.5)
With this new order of the equations it is possible to rewrite eq. <11.4) in
oa.rti tioned form a.s:
(11.6)
[Aee<QE)J=-~E 2 [Meel+iD£[Reel+[Keel
[Ae i <QE) J=-~E 2 [Me i l+ i.QE[Re i l+[Ke i l
<11.8)
[Aie<OE>l=-ne 2 rMiel+iQE[Riel+[Kiel
[A·II· <OE> l=-n-2rM
'¥t · · l + i.QE[ RII
II · · l+ [K II
· ·1
[Aee<QE)l~eeo+[Aei<QE)l~eio=feeo a)
( 11 • 9)
[Aie<OE>l~eeo+[Ai i<QE)l~eio=feio b)
(11.10)
(11.11)
By defining
<11.12)
(11.13)
(11.14)
e(z)=az2+bz+c (12.1)
lying in an assigned orientation plane (see fig. 12.2>. For the definition of
12.1 the eccentricities e 1 <z 1 >t e2 <z 2 > and e 3 <za> must be assigned in
correspondence to three distinct points z 1' z 2 and z3 ofthe element axis.
The origin of z a>:is is assumed in correspondence to the extreme left of
the element itself.
<13.2)
248 G. Diana et al.
in which Xc and Ycare the journal displacements and X5 and Ys those of the
bearing. For small relative journal bearing displacements, in the
neighbourhood of the static equilibrium position, the components of the
force due :to the lubricating film are given by:
<13.3)
where Xr and Yr are the relative speeds. If second order effects are
neglected in a congruent manner with the linearization carried out by
<13.3>, one can, as has already been said, attribute the ovalization to a 2
per rev frequency motion of a journal with a circular section < fig. 13.2>.
Therefore, by substituting <13.2) and its first derivatives with respect to
time in <13.3), we obtain for F x and F y expressions that are both a function
of the independent variables Xc' Yc' Xs and Ys and the ovalization, which
should be considered assigned: the first are already included in the
stiffness and damping matrices of the system while the second are put in
the right hand side member. The following expressions are found with the
above mentioned substitution:
(13.4)
(13.5)
<14.1)
The relation between !ei and Zee' as has already been seen in Blod< 6 (eq.
6.5) is given by:
(14.2>
F ·=[K J·lX
-J
·-P ·
-J -J
<14.3)
D·
~fj = _l. (14.4)
2
where Jj represents the moment of inertia of the shaft section and Dj the
external diameter.
(14.5)
Known, from the solution of the forced problem, the response leeo of the
rotor at the supernodes for different frequencies <Bloc!< 1Sl, from eq.U4.5)
it is possible to obtain the response in all the points of the shaft line. To
obtain the maximum bending stresses we used the shape function of the
finite "beam" element which permits, once the displacements and the
rotations of the end nodes are Known, to obtain the bending moment in the
vertical Mx and horizontal My plane of the generic section <see fig. 14.1l,
by means of:
(14.6)
The program recognizes the bending neutral axis and calculates the
maximum value assumed by the bending moment in each section of the rotor
successively calculating the nominal stress using eq. (14.4). These
stresses will give rise to fatigue phenomena should the rotation speed of
the rotor differ from the pulsation of the dynamic forces acting on same or
should the orbit of the shaft not be circular.
X "average" circle
ovalization
bearing node
+ z
Once the elasto-dynamic matrix of the entire system ~ rotor+ oil film +
foundation CE<De;>J~ has been calculated~ as described in BlocK 11~ and once
the dynamic forces F = F 0 exp<i0gtl applied to the rotor <see BlocKs 12 and
13l have been evaluated~ it is possible to calculate <see eq.< A. 18)) the
response of the complete system by solving equation:
(15.1)
The displacements in terms of module and phase, both of the d.o.f. of the
rotors of the line as well as of the connection nodes of the foundation, can
be printed or visulized by means of graphic output that show the
displacements along the generic d.o.f. as a function of the frequency or
displacements of all the nodes of the line <spatial deformation>, for a
certain pulsa1:ion of the excitation force.
commonly Known as the shaft line. This paragraph describes the method
adopted to evaluate~ by means of the frequency response of the rotor-oil
film-foundation system to a harmonic excitation force, the overall
stability of the entire line and the of single bearings. The method used in
the P.A.L.L.A. pacKage, called "forced method" is described in paper C23J.
This method is essentially based on the determination of the energy
P.A.L.L.A. :A Package to Analyze ... 253
(16.1)
(16.2)
X = Xe~t (16.3)
where
(16.4)
If r<O, coefficient of eq. <16.4) becomes positive and solution x(t) grows
e>:ponentially in time: the system is unstable and the instability threshold
is defined by the sign of r. This sign can be determined by also analysing
the steady state response of the system to a harmonic: excitation force. In
fact, if the system of fig. 16.1 is forced by an excitation of the type
F = F 0 e>:p<i0e;tl the response will be of the type x = Xaexp<iJle:t>: the
instability conditions can be verified by the relationship between the
energy e:ic: introduced by the external force in a c:yc:le and the total energy
e:Kma>: of the system itself. For De; = Wn <excitation force frequency
coincident with the natural frequency of the system> the value of r/rc: is
given by:
r Eic
-=- ( 16 .5)
rc 41t EKmax
Since the relationship between the energy e:dc dissipated by the damping r
in a cycle and the maximum kinetic energy is:
rc = 2rnc.On ( 16. 7)
and since e:dc=Eic' for Jle: = Cl)n we have equation (16.5>. The energy Eic
introduced by the excitation force is:
( 16 .8)
where Cf represents the phase lag between the excitation force vector F 0
and the displacement vector X0 :
2ar/r c ) (16.9)
1-a2
conditions. we have positive values of sin f· To define r/rc: with the forced
method it is therefore possible to use equation <16.5) which in explicit
form is:
r F0 tr I X0 I s i n cp
=-- <16.10)
ll2mDe21Xo12
being. as we remember • ..Ce: =CA>n• From a pratic:al point of view the system
is e>:cited by an excitation force with a variable pulsation until we obtain
the maximum of the amplitudes: in these conditions <resonance conditions)
the ratio r/rc: is determined from equation (16.10> and from this the
characteristic: parameter that determines system stability ol/Wd which. as
is Known. is:
(16.11)
(16.12)
where CM J. [R J and rKJ are the square matrices of order n and.! the vector
of n independent coordinates. The solution of eq. <16.12> is given by:
(16.13)
(16.15)
(16.16)
(16.17)
where, as has been said, CM J, CRJ and CKJ, are the mass, damping and
stiffness matrices of the rotor + oil film + foundation system already
defined in BlocK 5 and I= ,! 0 e>:p<iDEt> is the harmonic: excitation force. By
substituting solution:
<-D£ 2 [MJ+iQE[RJ+[KJ>~ 0 =£ 0
(16.19)
[E(.QE) l~o=fo
The advantages of the method used lie in the fact that it is possible in
this way to Keep account of the foundation and of the carrying structure in
general, through its mechanical impedances: apropos of this it seems
opportune to underline the difficulties of reproducing the
rotor-foundation behaviour due to the difficulties of schematizing both the
foundation and the casings of the machine. This can be obviated by
experimentally defining, in correspondence to the supports, the response
of the case and the foundation to a harmonic excitation force and from this
to define the mechanical impedance matrices which can be used, if
introduced in <16.19>t both in the calculation of the frequency response as
well as the stability analysis of the shaft line considered. For further
details see ref. [22 J.
258 G. Diana et al.
7.1.20 References
E10J Diana G., Bachschmid N.- "Influenza della Struttura Portante sulle
Veloc:ita' Critic:he Flessionali di A:beri Rotanti" - L'Energia
E:lettric:a, n.9, val. LV, 197S
[17J Eiraghi E., Falco M., Pascolo P., Solari A. - "Cuscinetto lubrificato
in Regime laminare e Turbolento : Parte 2 - Lubrificazione Mista
Idrostatica-Idrodinamica" - l'E: NE:RGIA E:LETTRICA, vol. l VII, 1980.
[18J Falco M., Macchi A., Vallarino G.- "Cuscinetto Lubrificato in Regime
laminare e Turbolento : Parte 3 - Analisi Dinamica"- l'E:NE:RGIA
ELETTRICA, vel. LVII, 1980,
[22 J Cheli F., Curami A., Diana G., Vania A. - "On the Use of Modal
Analysis to Define the Mechanical Impedances of a Foundation" -
260 G. Diana et al.
SENSOR TECHNOLOGY
.J.Toanesen
The Tedmlcal Unlvenlty of DIIIIIW'k, Lyaaby, DIIUIW'k
8.1.1 Introduction
Sensors of some type must be used to measure vibrations when deal-
ing with rotordynamic problems of turbomachinery. It is important
to select the proper instrument transducer or even different types
of transducers in order to cover a particular situation. One
should consider the following questions: "What is to be measured?,
'What are your looking for?, How is the machine built?, What access
area(s} is available?" and also the follow-up questions like:
What kind of display is wanted of the phenomenon?, and what is the
time element?". That kind of question will almost eliminate the
so called universal sensor, because there is no one sensor for all
jobs.
8.1.6 References:
J. Tonnesen
The Technical Univenity of Denmark, Lyngby, Denmark
8.2.1 Introduction
Experiments will always be needed, first at all in the laboratory
or at the manufacturer's facilities, in order to validate the theo-
retical model or to provide further physical insight and just as
often to give vital feed-back for new and improved theoretical
models.
The eddy current sensors are installed in the test shaft in two
planes each located 18.5 mm from the bearing midplane . There are
two sensors in each plane and they are arranged to give both in-
line and 180° apart readings. Figure 2 shows some details of the
installation and in addition to the displacement sensors is also
seen an oil pressure sensor and a thermocouple sensor, all mounted
flush with the journal surface. The signals from the sensors are
transmitted to the stationary parts via a mercury-cell slip-ring
unit and thus simultaneous information is available about the ab-
solute eccentricity, attitude angle, oil pressure profile and
journal midplane surface temperature.
l.~-~~H~""""~-----------------------;
5000 N
1510
J
mPa·s Do
m 4.25 12.1""""
l( 12.9 5000N
6 27.8 I
27.........
5000 N
95~------------------------------~
--CALCULATED
----------------__--_
90 - - - - MEASURED
', ~OOR~
as •
'-.....
~
_.
6S00 RPM -
_
'----------------.sooo RPH
sOOO RP~-
;t 10
z
a:
::I '---~----
•
.____ __, __..... -...
,._ ~
~65
60
ss ...--- .
a:------~--
_.. - --- J!j!!i_ ~'i
~~~2=00~0._~~~~0~0~~~~~0~0~_.~~~00~~~~~00
LOAD, N
8Sr-------------------------------~
SPEED• 5000 RPM
--CALC. LOAO•S600 N
80 - ...-MEAS.
60
results for the maximum bearing bush wall temperature which occurs
in the loaded portion of the bearing a little past the minimum oil
film thickness. The agreement between calculation and measurement
is good for the region where full oil film exists, however in the
ruptured film zone there are large discrepancies. A possible ex-
planation is that there may be a backflow of cooler supply oil
from the qrooves in the unfilled space in the striated film.
8.2.4 ~amic bearing properties
The dynamic bearing parameters are obtained by applying a force
function to the system. If only the stability threshold is wanted
this is aaxnp~ relatively simple by reducing the load or for a
given load increasing the speed until whirl occurs. Under these
condition is the whirl-frequency between 40 to 50 per cent of the
running frequency. The force function used in order to determine
the stiffness and damping coefficients of the bearing is either a
sinusoidal, step or impulse force each of which have their own
advantages.
8.2.5 Rotor-Bearing Tests
Once the bearing parameters have been identified in a test rig
where the rotor dynamics influence is minimized as much as pos-
sible the next step is to investigate the combined rotor-bearing
system. As the stiffness and damping coefficients of the bearings
vary with speed, the natural frequencies of the system will also
change and a speed may be reached where the damping becomes nega-
tive. At that speed, the threshold speed, the rotor becomes un-
stable in a self-exited whirl. Hence, from this point of view, it
is important to be able to determine the damped natural frequencies
of the system at the design stage. The information is, further-
more, of considerable value for monitoring purposes and in diag-
nosting vibration measurements, especially in connection with a
FFT-analyzer. The present trend in experiments is towards measur-
ing the onset of instability, methods of stabilization and the in-
fluence of mass unbalance as well as to determine which critical
speeds must be avoided, and which of them may be safely ignored.
The shaft whirl orbit is measured at ten locations along the shaft,
namely at each of the six disks, at the two bearing pedestals on
the outboard side and in the middle of the two bearings. At each
location are capacitance displacement sensors for measuring in the
vertical and horizontal directions. The sensors in the bearings
are flooded by the supply oil, and to insure constant dielectric
properties in the gap, they are additional flooded through separate
supply holes around the tip of the sensor. From separate calibra-
tion tests the dielectric constant is found to be a factor of 2.3
less than that of air (at 50°6
with a variation of ±3 per cent in
the temperature range 20 to 80 C. It is found to be insensitive to
pressures up to 6 MPa.
I
2
~
2 3f
l.
~,If ~\ ~ 5
II II
ii!
"
IIU lJ IIIW
~~~QUENC~~Hz
:311' I
f 3
! f ROTO ~ SPEEDol7160
1
(UNS BLEI
R~
ii:
6Kt 2
2f
5
J ~ \} I Ill II
Fig.9 Test Rotor Frequency Response,
with no External Damping.
8.2.7 Stability Test Data
Some typical test results from a stability test run are shown in
figure 9. The test rotor is well balanced and should theoreti-
cally be stable up to the max. speed of 20.000 RPM,when the rotor-
bearing system consists of a rotor without the overhung disk and
with rigid bearing supports (no squeeze film damper). However the
calculations did also indicate that above 12.000 RPM the damping·
is (log.decrement) small, typically below 0.08. As may be seen
from figure 9 the rotor became· unstable at 17.160 RPM. There-
sidual unbalance in the rotor results in a bearing whirl orbit
which occupies approximately 20 per cent of the clearance just
before the instability sets in, say at 16.800 RPM. At the thres-
hold speed this whirl orbit grows rapidly, occupying about 50 per
cent of the clearance circle when the speed has been raised 50-
100 RPM.
10~
0 RPTOR SP ED= 200 PO RPM
0.. (~TABLE) "-..t
w I
a:
~J(t
1/)
w
a: 8
0..
..J 2 7
0 10'
)
u 6
3
~ s
""
J
z
~ ~4 .ftJlfi \l1500
uf0 1>0 200 300 400
FREQUENCY, Hz
• THEORY
• EXPERIMENT
8.2.9 References:
L.Gaul
Unlvenlty of tbe Federal German Armed Forces, Hambura, FRG
9.1 Introduction
The prediction of machine vibrations by theoretical approaches
as well as the modification of response after construction often
require taking the interaction between machine, foundation-struc-
ture and subsoil into account. Three examples are given. Fig. 1
shows a discretized model of a drilling machine with long founda-
tion slab on soil. The impact of static soil-structure interaction
was calculated and measured by Thurat [22]. A base for the dynamic
analysis is given by the soil model in the present paper.
A multi body model of a forging .hamner (Fig. 2) is covpled with
a viscoelastic truncated cone model of soil (Knobloch and Gaul
[17]). Thurat [22] calculated and measured the transient response.
Novak [19] treats a hammer foundation as a system of two masses
on a viscoelastic halfspace including embedment effects.
The global response of turbomachinery frame foundations e.g. the
low-tuned steel foundation with concrete raft (Dietz [7]) of Fig.
3, or the response of block foundations are calculated and studied
experimentally by small scale models in the present paper. Dynamic
284 L. Gaul
folXldation
Oil,
Piston
Pad
Joint {
Pad
I {U}
[K(iw)] I
I 1=
t{uc} .i
I
{P}
I
-{Fc }
II
[K(iw)]
l
r
{Uc)
II
{U}
II II
-{Fc }
= II
{P}
(1)
and
I I I
s {Uc} - { vc
} {F c}
II = II
[ K( iw}) I I (2)
{Uc} {Vc} {F c}
t Kr Js {Uc} {0}
[ld
. K ( 3)
l II = + [K]l II
{U} {0} {Vc}
II II
{Uc} {P}
f ~~~~~~~:-~~~~~~~~:-~~~-~~~~~~~:-~~~~]
(-<~l[M'AP]T+ [hJ:-w
[Kw>]) J
{~~~} = {~~~}
2 [Mtl>) + [Ktl>) {q>} fT}
(5)
Interaction Between a Rotor System .. . 289
Soil elements
are decomposed by
(7)
in two wave equations for the scalar and the vector potential ~
and ~k respectively
290 L. Gaul
(8)
-
CJ ••
lJ
( 12)
with A, Bk determined by the boundary conditions. The stress bound-
ary value problem is solved by superposition of harmonic vertical
displacement waves at the halfspace surface z = 0
-s '(S,y,w)
w( x,y,O,t) = H w -
p(S,y) exp.[i(Sx+yy+wt}] (13)
generated by the stress wave
-o 2 z(x,y,O,t) = p(x,y,t) = p(S,y) exp[i(Sx+yy~t)](14)
1
able measured data, re-
Shear Dilatation sults are presented for
G*liW) the constant hysteretic
~·-· Eoliwl solid and Kelvin-Voigt-
solid leadin9 t~ damping
head factors ns ,A (w) = ns ,A
and nS,A(w) = a0 ~S,A
respectively, where
Fig. 8 Clay sandwich for alternating ao = wa/vs is a frequency
dilatation and shear tests. parameter with character-
istic lenth a.
The obtained harmonic solutions are now superimposed by integrating
with respect to the wave numbers B, y in Fourier's integral theorem
1 oo -s w -
w(x,y,O,t) = 2n ff H ' (B,y,w)p(B,y) x exp [i(Bx+yy+wt)]dB dy (15)
-oo
with
- - 2pj(sin saj sin ybj\
p(S,y) - -n-\ By )
being the two-dimensional Fourier transform of the exciting stress
field (Fig. 5) at one surface element of area 4aj b·. The elastic
halfspace leads to improper integrals due to poles ~f the compli-
ance in Eq. (15) at the shear wave and the Rayleigh wave speed.
Solutions can be obtained by choosing Cauchy's principal value and
performing a contour integration in the complex plane. Here a diffe-
rent technique is used. Because the viscoelastic halfspace yields
finite resonance amplifications instead of poles (Fig. 7), the in-
tegral of Eq. (15) is no longer improper with respect to the inte-
grand and can be integrated directly.
292 L. Gaul
HALFSPACE
H H H
{Fs} =[Kitus}
=-
8
H H
H H
{Fe).(ue}
EXCAVATED HALFSPACE
H H H
{Fe!= [KeJiue}
E E H E -1 E }
{
-[Kee]-[KesH [K]-[KssD [KseJ
Base
plate
(16)
A
-for the base plate with mass m, inertia tensor J 1J
..
Interaction Between a Rotor System ... 295
A q>.
··
J .. + m e. "k X. uk
c ·· = S ( o. a: s a:s)· ( B o
L e. "k X. Pk +f~. - M. +e. "k X. Fk
B)
1J J 1J J o.= 1 1J J 1 1 1J J
mf·· c .. ) S o.p s F B ( 17)
"\ui-eijk xj ~k = - o.~l i - i '
-for d
the rotating disk with mass md , inertia tensor 0ij' angu-
lar velocity w
d .. d . . d
0ij o.j + 0 11 w(oi1 o.1- 0i3 o.3) =- Ti
..
md wi =- Kid + md e w2 ( coswt o1i-s1nwt
•
o3i ) • ( 18)
K.d w£
1 = [KH] d (20)
M.o <P£
1
K.o v£
·1
In the frequency domain Eqs. (16) to (20) lead to a substructure
equation like Eq. (1)
[ K( iw)] (21)
296 L. Gaul
a 0 = 1.0 v =0.4
lis = 0 .2 11/lls=l
mr.n I p(a 1) 3 =10
JI.n
IC•,YY
/p(ai)L20
-
0,54
:a
the response of the rotor
on a rigid foundation with
the response corresponding
..I
2
to a flexible foundation .
.
i Analogous to the frame
! foundation (Fig. 13) the
a
a.!I first three resonant ampli-
c fications are predominately
due to the foundation,
while the last two are
governed by the rotor. The
foundation influence splits
one resonant frequency of
the rotor in two of the
combined structure with
lower amplitudes.
1 2
Frequency parameter a,
Fig. 13 Response of frame foundation
on soil.
Deflection
scale
111t•11t2 '151J.m
--Rotor on flexible
----Rotor on rigid
w, foundation
/e
2 3
Shaft~ D~
.60 .0153rnm
iii
~ . ~8
....
N
.36
...
QJ
"0
::1
-a.
- - ROTOR ON EL ASTIC FRAME
AND RIGID FOUNDATION
E . 2~
< - - - - ROTOR ON ELASTIC FRAME
c AND ELAST IC FOUNDATION
~
..... .12
QJ
>
00
"' -
100 200 300 900 1000
w rod/sec
--- -- - -- --- --
E/E
r-- -- r-- --- --~· - - -
~
1\
v IW ?
I
A._ -.,..,_
0
...,..,....v -
..........
00 LIN X HZ 400
X. 107 HZ TF 1.44E-1 EU/EU AVG N 1203
00 LIN X HZ 400
Xo 107 HZ + 79.6 DEG AVG N 1203
[11] Gaul, L., 1980. Zur Dynamik der Wechselwirkung von Strukturen
mit dem Baugrund. Habil. Univ. Hannover, June.
[12] Gaul, L., 1980a. Dynamics of frame foundations interacting
with soil. Journal of Mechanical Design, Vol. 102,
pp. 303-310.
[13] Gaul, L. & 0. Mahrenholtz, 1981, 1984. Dynamische Wechsel-
wirkung zwischen Maschine, Fundament und Baugrund. Ar-
beitsbericht zum DFG-Schwerpunktprogramm Betriebsverhalten
dynamisch belasteter Maschinen (not published).
[14] Gaul, L. &M. Plenge, 1983. Ein Baugrundmodell fUr geschich-
tete BaugrUnde mit viskoelastischem Stoffverhalten~
ZAMM 63, pp. T50 - T53.
[15] Holzlohner, U., 1969. Schwingungen des elastischen Halbraumes
bei Erregung auf einer Rechteckflache. Ing.-Archiv, 18,
pp. 370-379.
[16] Huh, Y., Schmid, G. &M. Ottenstreuer, 1983. Evaluation of
kinematic interaction of soil-foundation systems by
boundary element method. SMIRT 7 Chicago, August, Paper
K814.
[17] Knobloch, W. &Gaul, 1975. Dampfungs- und Federverhalten
elastischer und viskoelastischer GrUndungen bei harmoni-
scher Erregung. Kolloq. Viskoelastische Systeme, VDI-GKE,
TU Berlin, pp. 412-451.
[18] Kramer, E., 1984. Maschinendynamik. Springer Verlag, Berlin.
[19] Novak, M., 1982. Response of hammer foundations. Proc. soil
dynamics and earthquake engineering conference Southampton,
July, pp. 783-797.
[20] Ottenstreuer, M., 1982. Frequency dependent dynamic response
of footings. Proc. soil dynamics and earthquake engineer-
ing conference, Southampton, July, pp. 799-809.
[21] Sarfeld, W. &C. Frohlich, 1980. Dynamische Wechselwirkung
von Gebauden und Fundamenten auf dem elastisch-isotropen
Halbraum, Bauingenieur 55, pp. 419-426.
[22] Thurat, B., 1978. Machine-Fundament-Baugrund.
Diss. RWTH-Aachen.
[23] Waas, G., 1972. Linear two-dimensional analysis of soil
dynamics problems in semi-infinite layered media.
Ph.D. Thesis, University of California, Berkeley.
[24] Dominguez, J., 1978. Dynamic stiffness of rectangular foun-
dations. Dep. of Civ. Eng. MIT R78-20.
CHAPTER 10.1
ABSTRACT
10.1.1 Introduction
e) Out-of-phase synchronization.
10.1.2 Notation
c Damping Coefficient
Hz Frequency, Cycles/Second
I Inertia
K Spring Constant
T(t) Time Varying Torque
e• Angle of Twist
9 Rate of Angle of Twist Angular Velocity
~ Angular Acceleration
(I) shaft Angular Velocity
., ~
Two such CMS units currently in use are the Torsional Stress
Analyzer (TSA) developed by Kraftwerk Union AG, as described
in reference [7], and the Torsional Vibration Monitor System
(TVMS) developed at General Electric Large Steam Turbine
312 D.W. King- N.F. Rieger
Sample results from the TSA are shown in Figure 8 for a TG set
with four couplings, where the shaft motions are monitored.
A code, shown in Figure 9, is used to describe the severity
of the incident and recommended act ion, if any. The second
column shows a figure for the shaft life expended in the
incident, based on a 10()4K, shaft life expenditure at first
crack initiation. The third column shows the cumulative
fatigue incurred since installation of the TSA. The accuracy
of this model is evaluated in Figure 10 for a 970 Mw, 3000
rpm, six-flow TG undergoing a sudden three-phase short circuit
(A), and for partial load rejection (B). The shaft torque
between the last low pressure turbine and the generator is
compared using the figures from the TSA and figures obtained
throughout the application of strain gages to the subject
shaft area. The results appear to endorse the TSA fatigue
life model.
The GE TVMS was first installed in early 1981. The TVMS uses
both air gap torque and angular velocity for its inputs and
produces figures for fatigue life expended in much the same
way as the TSA. The novel concept of the TVMS is the use of
Data Acquisition Systems (DAS) which sends in information
from the TG unit to a central computer to be analyzed. This
permits centralized monitoring of many remote TG units, and
tends to minimize the need for specialized test personnel at
the TG locations.
Problems of Turbine Generator Shaft Dynamics 313
typically the first four modes. These four are therefore the
primary modes which must be considered in such calculations.
Current problems associated with the use of CMS units are the
lack of a means for retrofitting the past fatigue history of
the shaft into the cumulative fatigue calculations, plus a
means for accounting for the effects of modal damping. A
means for specifying an appropriate threshold strain value
that will provide optimum use of these diagnostic tools is
also needed.
10.1. 7 References
...
t' ...
.•. L
............, r.ntt
•r•.... ..._"' unollf
Current
Transformer
. !
b
I • Rogowski Coil
Bushing
~~~-~
r-" : ______:---,
I : : :=-:-:-rr- ---=---, 1 I
I
I
I
I
I
.,.~,., .
Et.c••onc Modl'l of 1tw 51\ah System
I
I
I
I
I
I
I
I
O .. «t ~· M.f'nOiy l'roc•M Coonpu1C"
I
s..... f•ltQwl'
Pl'f lnc:•dent o.p .. of o.s·- "'
I
F,
A«t.tmula1ecl
t'.owph,. S.t•tor W•nd•na
I
I
Dtforrn•tiO,. S1rftt'"l
SNh F•t•Qvr
L_ ].F
__j
,-~~
L:::J
Figure 4 Block Diagram of Torsional
Stress Analyzer [7]
320 D.W. King- N.F. Rieger
.
w,,('
n
Wrs
l
HPT
JR
w
~ Ill
m GENERATOR
m
w
EXC
IA,EA
Ia. Ea
lc· Ec
- TORSIONAL
VIBRATION
- ELECTRICAL
FAULT
MONITOR MONITOR
J
DA.S.
CTVMDASI J D.A.S.
r-1
CEFMOASI
FAULT
I
l TORQUE
TRANSDUCER
SENSOR
l
rI I I I CON~RTER I
TRIGGER
~ CON.~TER J TRIGGER
MICRO COMPUTER
WITH
MEMORY
J MICROCOMPUT ER
WITH
MEMORY
I
I DATA
COMMUNICATIONS
I
l SYSTEM I
~ ~
_j TELEPHONE
R LINE
REMOTE COMPUTER
DATA ANALYSIS
DATA STORAGE
r --
1
I
!g{:(J·txHXJ:q,.§
.
I'------_,....::;...-~~-::::.:.-_:-_-_-_-.,:-_-_-_-;.--:_-:- ___;
."
..
~
C"
0
'i
.c.
"'
,;
.e p.u. (B)
•rr o.J
1-
"
0
~'--------M-e
•
I
Q.
~
c I
p.u.
O.!;'
,;
.e
.,•
~
!·:!,..,. . . . . . .""~
~ p.u.
Timtt
1.0 9.5 HZ
9.3 HZ
-1.0
--18.5 HZ
t.O - 18.3 HZ
-c
II) 0
~
.
Ell'F"::r:::::fu- '(...U
,~
....
·JIE5~-·· ···~
E ',
...,
II)
0
0. -1.0
..
0
~
1 ·0 ·-- 21.4 HZ
- 20.8HZ
iii
a;
..
a:
·;,
~
"ii
...0E -1.0
z 1.0
--· 40.4 HZ
-39. HZ
•• ...
...
II':' LPI
Ju1
. . =-
44.1.:.,
.....
~
I liP
,...,~
.... ,. .
~~t
-l!!!!JII"f __.;.~
-~=~~jili"r",
1 . .1 T . . - •..._...._. eM••
....-.-r••c-I
---·-1 _ ..··-··--i. . . -,
1 •
I
.... u
-~-f~~-f=-----~------------4~··~
d ,•
,..............,......... ,.,..,
•IMit I
-~~~t• ·~,,uc
''"
--·...·vNJW·fri~~J!!;:
,1 ~-· ··"····
~l'' 1 i'~'AI"l"'''f'
:
..
~ ~ ·lilll!~A!\~: ~
1
"'.........
!"In\"'"• .. " • ·• ...
4P
•Hu·•·•~····-,., ..
J J.. su·
1
.,.," .... " • -...,~.
~.... •iltmn~ii;,.~ ~~~ Four Modes
...
.. ~ __ .. 'r ........
ntn,Utlutnt IIIUU•t•tttUtff'f''''''' tt SLC'
IFI ttl t I tttl I I fIt to
I :
....
-cr{Xt{XJ--{):J-{J~
11_
0~~-----~--~~~------------
/
c=::=:::::..__ Mode 1 15 Hr
.. -1
~e 2
t_:r
1
~ ["""---~ 24.7 Hr
~
~ ~~/\ ~e3
35.6 Hz
~-~~
J -~pi 1 :
Mode 4 60.1 Hz
.
~ 0~<=:::>
;;;
i ~I /
Tors10nal v•brallon modes ol J 970 MW. 50 Hz nuclear-lueled
turbme-generator.
1 . - -
....
L •
---~.rv'....:v :: ;_ __ ...:~
, i I '
I
oo'- -- --· +- - - ·. - - ... --- - -
0 - - - - r - . - - --
";'tOO -· · - ~- --- ·- --·-·' : :
~
-
["
;;I()'!.. -·- ··-
• ~-
0
:I
~.,.. f-- - . ...
-
: o.:....
- 1·
' --'-------
. :t-
J .
w
;;oo
2
)
1- --
tI
Four Modes ~o:_ . -.,~~-.'.'.!:..:.2
.e 1 -
cr --
I ~ oo·-I
0~ l ~
; "'... 1-·
,,
1·
0· · ·· - - - -- - .... -- 4 h
123• \2l4 1234 1234 12'34.
0 ,.... - · - · - - - · --- -- 5 loiOCS.
~
;;
~
U)
Figure 14
Problems of Turbine Generator Shaft Dynamics 327
(A)
2.01
- 1.f!l
-
~
41
-~ 1.0
•.....
0 1 2 3 4 5 6 7 8 9 10
Time (sl
0.15 i
!
~ C.1~
1
r I
t I
ioos~
i r--,-______
I Dm3
i 0
500
0 100
Time (c,cles)
0 60 eo 100 120
"X.OF FW. LOAD
1
~ - ~:-· -:-· .
1 ~--~~- ~~--
: 'J.."'T<-"•' • -=-:;
I
0 .. , ....., ~
=·i
t •t ,
! !
y
. -=1
=~
:
I
!'
-
I ( ·- , ,
'I
ic:::b
-- - - --~--
• j .......... ... , ....
I ... .....,.""1 t-::;::,:ll"'t'U, ... ·-
*
iF:
:.-- ·.,,1 ....,.. !1-_'G-• ·.~)
IJn, · ~•C.cv••
,..,,""'....,.."' - - ·····
. ... ....,.. .,~ ._..:., .
s...,.f'S ... -.t....~
\JIOI· ..,....,. . c.or.,., rt":..::..;;,.~ :;,.,-s; ~r:.- t.~~
;** '- -
-
Catc.ulaled Fr<equency Meeauted Fr~uenc.Y
Uo\4 H'""t..t ViMtion (Hz) (Hz)
Rol;ng of Number or M•~o5ea tn by Double- by Frequency
I"'-Wl fo-e. a tOt System Fr•quenc, An•t,'liS from
I E•c•••••on
f'olu SuddenSho<t -
n <10 n >200
(S. C .Tut) ~rcu•t Test O.tt..
.
630 2 1
I 17.0 16 . 5 - 16 . 7
660
• 1 7. 8 8 .2 - 11 .1
1200
r--- • 1 6 .5 6 .5 - 6 .3
1170 2 2 24 . 3 24 . 7 24 .7 24 . 5
3 34 . 9 34 .6 Ji. 3 34 . 3
CHAPTD 11.2
P. Sclnrlllllla•, R. Nord....
UaiYenltJ of K.-..ue_.., K......._...., FRG
ABSTRACT
10.2.1 Introduction
r-·-·-·-·-·,
I I
NETWORK ---~j TrQnsmossion
I Iones
,L....--..-1---1
I
I TrQn sformer
GENERATOR I
Peril" s EquQtions L·--·- ·--·_j
Turbine VoltQge
govenor regulQtor
r· .,
i
I
i
i
!
HP IP LP1
! Gener Qtor -Rotor , j
L._._.-.- · - . - . - . - . _.1.....,:-=.:·.-=-··=-···--' -· - · j
TORSIONAL VIBRATION SYSTEM
Finite Element or Doscrete 11Qss Model
TURBOGENERATOR
SET
Ip
HP
FINITE ELEMENT
MODEL ( 250 DOF)
TORSIONAL
ELEMENT
K GJ [1· 1 -11 ]
- T
-e- T
Hnss-Hntrix Stiffnus-Hntri x
D =a M + 8! (2)
~. .
~~ ·.~[I
I •I
·I]
I
I· · I i•l··· I
CoiiKtdtnct
Tronstoraotion I: •lA ftH I
I;'· !1'·1:
r:' ..l'·I:
fiNITE ELE"ENT
!!1m!.
To solve (3) we'll use the modal analysis method; but before we
study the corresponding homogeneous problem <r(t) = Q).
(4)
336 P. Schwibinger- R. Nordmann
).<:.
"' ( ...... )
~\ = tile ( 5)
-; -\-
(K - Y. :.:) cp
-
=0 with I = - ). 2 ( 6)
f.=-
l
1
2':1'
£ l
(7)
f
-m
e
= [-f m1
= Me .. e
!l + -
Ke
.9.
e ( 8)
fm2
Analysis of Torsional Systems 337
i- · - -~ : ~
57.22 Hz
I l
=ljt ,._,__.\,....,.
L.,-··11~\l : . •.
,_,,...,
7:
rt til /
.,c:.'"="""=----- - _fm4
'LwiiJ"
. ~; .-""ill'i'"'-
~s ~-- _ ....:;:.cJ
118,72Hz ' -
+ ------~~~~;,; {;·;,::l fms
f~, q~
-
19e(~),(jie(~)
® --
q: f~2
L~e
.IS;e .Me
e e
q. = cp •
-1 - 1
resp. (9)
..e 2 e
q. = -w. ~ .
-1 1-1
338 P. Schwibinger- R. Nordmann
( 10)
( 11 )
des , .•
-1.
N
g_(t) = i=l
r ,.n.(t)
-1. l.
=1 n<t) ( 13)
( 14)
(15)
m.n.
l. l.
+ d.n
l.
+ k.n.
l. l.
= g.(t)
l.
(17)
with
T
m.l = cp . M cp . - Modal Mass
-l --1
'I
d. = cp. D cp •
l - l --1
- Modal Damping
( H!)
T
k. =cp. Kcp.
l -l- -l
- Modal Stiffness
g. ( t)
l
=-cp T•l -F( t) - Modal Force
n
s(t) = E cp. n.(t) ( 19)
i=1 - l l
t -D.w. (t-t)
n.(t) =- I g. ( T ) e l 10 .
SlOW. ( t-t ) d
T (20)
l rn.w. l l
l l 0
d.
D. (22)
1 = 2m.w.
1
1 10
w.
10
=~
1 1
w.1 = w.10 (23)
and
342 P. Schwibinger- R. Nordmann
k
s.
1
= + (25)
n.•1k+1 a21 a22
•
n.1k b21 b22 g.
1
k+1
n
qk(t) = E , 'k n. (t) (26)
i=1 1 1
lO
.. , rn r ad 0
Alil\ Aool'\ D
@ 'lz ol ,,.ul ('{ , l'G'~' ',, .I F' '• · I ,,}u·, ~H P, \/\..1 -: · 0 0
1n 0 q 2• l.ll/ "" "loa W "" - 'zoo®
T:z · l .. :,,.,~, '·,• ..: · I
~CJ '•I'I'• i, ltt.i),
I - lo .. 30
I
F I
0
~
•n~
Y]
I
II @ 0 qlk 0
T. C> = ·~ ,... = ..... 200 Q)
c:::::::::::::> c==T">
- )0 - 30 •
10 @)
mrad
0
I & - -
q. r "' YiiJII 'i!roo v v .... '{zoo
-10
'!'l".e ar.-.;:E:~,;.:Oe ~! -:.he thir:i ::aocie at. the cout::ling ic ~~~.al1 9 i:c
:::::r::·.:.:br. :c :l".e :.:;;bce:::en~ 1k ~he1·efore indgni ficant.
For :l".e :ransiera: ·r:t.ra:ion ..-e can 1evelop the c.:la::tic for-
c~s at an ele=en~ e in a si~~liar ~ay as w~ have done it abov~
- -.
Torsional
Moments DisRlacements
,:~,,
Lr.\. !L
~ '5' ~·
K'.t!
••
ltft nodt
c
At a torsional element act the distrubuterJ moment r A9 the: lli~-
e
crete moments rEi• the inertia moments and the internal momentG
at the nodes. Assuming the distributed moment as constant in an
element. we get for the internal moment&
[
1
- ~i ]' (27)
(.
1
Analysis of Torsional Systems 345
-f
e
e s1
f (t) [ (28)
-s
fs2
e n n
f ( t ) ~ E Ke cp ~ n . ( t ) = E f ~ n . ( t ) (29)
-s i=1 - -1 1 i=1 1111 1
-M.
e - -·- M11 •Mn•Mn•M14
...
z
i:
18 600
10 000
~ .,
e0
::1:: 0
8
~ -10 000
.:.:; -14 300
~00
I Time t/ms
COUPLING MOMENT
-
E
z
..ll::
~
E
0
l:
Cl
.!:
Mode c.
:l
0
u
E
.~
-4
r
I ~ ~ I
Mo 0·1
.t!lk
3
Jhv- Ah~o" ~yo..~oo ®
v
~
"8
I
I
I
I
-
l:
~-:~
I
~ I
I
I
v·ULP 1o't <uP' - - ""' 2oo (]) 15 I
:.;:
il E I
II
-~ 1
1
·;:: E I v
.&_k ~ § I I
I II
u
E I II
Mo -----------
Or·- - - - - - - - - -100 200 ® E I v
1000 E I v
E I I
~
-4
@ I
E ~~ I
I
~~
~ ~~ I ~~
t ~·
IE
@ ~ ~~ II
•
1:!JL3 II I
Mo 0 lmrm._,..\nnf"\nruonnn1"' ~~ II ~~
~~ 2
0 123412341234 1234 1-34
Mode
M0 = 2182 kNm
"'
T F·h(t) ( 31 )
g.(t) =
~
<p.
-~ -
where h(t) is an arbitrary 'time function'.
tn. T "
F t -D.w. ( t-T)
:!~ -
n. (t) =---
~ m.
f
w. 0
h( T) e ~ ~ 0 sin w.(t--r) dT
~
~ ~
(g)
n. (t) ( 32)
~
or
T F" (g)
'P.
n. (t) =--- n. (t)
-~ -
~ m. ~
~
(g)
where n.(t) is the modal response of the i-th SDOF to the exci-
~
T
- <il· F
w =--- -].
(33)
i m.
l.
Ill
"·-
J.
'1;-
~D•const
""~ I
I
I
1 t -D.w. (t-T)
n. (g) = max{-- I g(~) e l. 10 sin w.(t-~)d~} Sd(w.l.O ,D.) (34)
1.max w.l. 0
l. 1
350 P. Schwibinger- R. Nordmann
~m,
Fig. 13 Response-Spectrum
An example
n
I~~max., i~1 I ~i Imax (38)
But because the maxima in the different modes don't occur at the
same time this approximation overestimates the solution in ge-
neral. The experience shows that a mean sguare approximation
yields much better results. We use
IS Imax * I~.
1=1
I9i 1max
2 ( 39)
w. -w.
JO 10>0, 1 (1 S i < j ~ n) •
w.10
I9 Imax tv /1 9k 12max + 2 r
i~j I~i Imax I9j Imax (40)
w.,)0 -w.10
~ 0,1 (41)
w.10
is valid.
Analysis of Tonional Systems 353
5 100
...•
N
'o
4 80 o-o.oos • O•.O.OS
0·0•0,01 • 0•0,1
• 0•0,02S
3
-z 60
E
!40
E
!: 2
No
'3 ....
• '30
i 1 •
A
-
"' 20 0•2u60Hz
J. ,tl~ i"' ~ Wo
0 0
~o I 40 60 80 100
I
120 140
I 160
w,.wzo wlo '"'t.o '"'So w6o
,...... W./0
0 1 2 3 4
• Time-History-Method
• 10% Response-Spectrum-Method
• to check - the sum of maximum response of all modes
h(tl
'a
...E
Cl
.....
go 60
--Time-History-Method (THHJ
------ 100fo-Response-Spectrum-Hethod IRSHI
-·-·-·- RSH: Sum of max. response of all modes
yields quite good results in all shaft regions except the gene-
rator. Here the approximation yields to small displacements com-
pared to the TH-Method. Probably in the response the maxima of
two or several modes supperpose here. It is important to note
that near the couplings, where the maximum twisting and the maxi-
mum stresses occur the RSM apprcximation is fairly good.
EIGENVALUEPROBLEM EIGENVALUEPROBLEM
-!l!;,w;,mi 16°/o 16°/o -!;,wi,mi
100°/o
-17°/o
Fig. 16 Comparision of computer time used by 'Time-History'-
and 'Response-Spectrum-Method'
10.2.6 Conclusions
10.2.7 References
N.F. Rieger
Stress Technology Incorporated, Rochester, New York, USA
ABSTRACT
(1)
In matrix form:
+
[~ =
4 2 =
I1I2(1) - q(I1+12)(1) 0 (3)
(1)1 = 0
(4)
(1)2 =
f•H1+I2l] 1/2
~ I1I2
q - I (1)2
( 1 ) • 9 (5)
q 1
o Disk inertia:
= 16.00 lb.
362 N.F. Rieger
(36)
= ----
g 8(386.4)
= 0.186 lb. sec. in. 2
o Shaft stiffness:
o Natural frequency:
= 40.52 X 10 3 rad./sec.
o Angular relation:
I 2 I
'
llw I
92 = :t q i 91
!
,_ .i
r (0.186)w 2 .i
It- 91
!
L 152.7 X to 6 .
Assuming the shafts have zero inertia (or that this has been
included in the disks), and neglecting the inertia of the
gears (close to node in mode 2), the equations of equilibrium
at the gears are:
1/2
L
2n ] hz
and:
and:
F =
o.iKiK.4I4R2
91 = 2 2 2 2
11q1(q4- 14(1)2)R3 + 14q4(q1- 1 1(1)2)R2
o.2K.lK.4I4R2
94 = 2 2 2 2
11q1(q4- 1 4(1)2)R3 + 1 4q4(q1 - 1 1(1)2)R2
1
z;D =
J2(2+V) (l+V)
where:
c q
z;D = and 11)0 = (::1)
21 2w0 11
Torsional Dynamics of Power Transmission Systems 367
c
z; =
Amplitude ratio:
Torsional Dynamics of Power Transmission Systems 369
\) =
Q =
or
T
0
a = -------------2; b =
qq - (11+I2)(1)a q
Problem Details:
Machine Specifications:
Inspection Details:
Diagnostic Procedure:
Proposed Remedy:
Effectiveness of Fix:
Problem Details:
System Parameters:
Inspection Details:
Possible Causes:
Diagnostic Procedure:
The primary cause (or causes) of the timing gear slippage was
sought. Incidental possibilities were attended to as part of
the correction.
Primary Remedy:
Secondary Remedies:
10.3.8 References
~----------------------~1 (a)
91
Mode 1 2
w1 = 0
Bz
/1 f - - - - 1-
lU Shaft 2
Il
I
Shaft 1
'
I
~System ~latural
• Frequency f 1·
n
l/
y
........
s aft 1 ........._
I
~ I
- -.,I
f Shaf~ 2 I
1--- I
I
.I
I
_Speed fi RPi:
Natural frequencies w 1 • w 2
Torsion" I
Amolit!Jde
ew
(a) Undamped single degree of freedom torsional system
r, Iz
ql
e,
D-{J 62
Figure 7
Torsional Dynamics of Power Transmission Systems 383
screw
'>?~9='9='91:9!~~~7'""'~~i--N o t u ra I
rubbczr
joint ring
Plug
Natural
rubber
0~~~~~E;;;~lj oint ring
screw
c
(a) Simpl e torsi onal system with tuned damper
of damping
(b) Respo nse of simpl e system for vario us amounts
Figur e 9
Torsional Dynamics of Power Transmission Systems 385
v "'-../
1.0 ~
"" ~ t--
V- r--
y I'-. JI -
.L.O
~--
ELASlOMER ORIGINAL
COUPLING STiffNESS COUPLING Sllff NESS
FRE DUENCY•
CPf.
o~------~~----~~-----L--~----
105 4· .10 5 106 4 • 106 K
COUPLING STifFNESS LB IN./RAO.
~ 4000~----~----r-----t---~~-r--r---~
u
>-
COMPRESSOR
u
~
::I
3000~---4-----+----~----t--r--~~~
0
w
a:
~ 2000i======~======~=====-+-~....,-;--;---,
Q
~
2S IOOO,f===J.~=s;~~==F====9f==i===r====1
~
Sf COIID I DLU
57 T
Figure 17 z
;,
:::0
;·
~
..,
Torsional Dynamics of Power Transmission Systems 391
10
1-
1- I I
I
1-
1-
r-----+-~--;------+--i
~l I
I
II'
' ~J,,~~
1-
1- I
lsi CRITIC t\L
I
I 2nd CRITI AL
b
I1
,!
J
~ I
j
z
<t
0 1.0 ~
I
<t
r-
\
a:
\
~
....J
....J r-
:E
\
-
I.U
0
:::;)
1-
\I
1\
....J
a..
:::!:
\
1\
\
<t
~
\
g
t-
0.1
~TORSIONAL AMPLITUDE
I \
1- RESPONSE AT FIRST
r- COMPRESSOR.
rEXCITATION: 2-PER-REV
r- SECOND COMPRESSOR \
GEARS.
r- ORIGJt:AL MACHINE WITH-
OUT QUILl SHAFT.
I
'
L
Turbine ~~peed, rpm J
Table 2 Analogy Between Rectil inear and Torsio nal Vibrat ions
-
Symbol Unit SJD1bol Unit
Time t see t
rae!
;'
in
Displacement "' rad/sec
Velocity
"'· in/sec
in/see' ; rad/seel
Acceleration "' K in-lb/rad
Spring constant k lb/in
in-lb-see/rad
Damping coefticient e
r
lb-see/iD
dimensionleas
'r dimensionless
Damping factor
Maas ... lb-see'/in I lb-in-sec2/rad
1
~ Wn• ff
~ ll,)~ + r t.)~ :!:1 (w,.._I 1.)~)1__ +~u~ J
~
c..> D [ I
2
1 2
w~ :a k.fr ; w:- • KLJr,_ i r· lr.•I,Vr,
3
H 11 12
w,•O; Wz =u., = k(J, •Iz) 1 I,Iz.
1
4
~r:f:J 2
w,•O i (..)~ • 1\l k,k,~ (I,- n" I1)
I r.rLCK,- n•J<.&)
5
~~ I, w4 [ I,I 4 (Iz +n..I3 ~- ua.~f,Cla. +Iln.. •I.n1.)
11
!!~f.! • k',f4 (!,•I1 • n'f.,l + t,k'.. (I,•I.~ •
n•IJ •I.. n•) =-0
iii 13 1!5
11 12 13
6
~ w"(T,T2!3 )- tu'l[K,rl (rz •I 5)+K,!J(I,•I~,U
+ k,l<z. (I, +I1 •I~) =0
k,-I,w'l.
~
- 1(, 0 0
-k, l(, .. k,-r, ...... -I< ... 0
7
0 -k.. ka +k:s- rlldz =0
11 I,
-1<?»
12 13
0 0 -k3 K3 -I1 w"
'(}
Table 3 Torsional Natural Frequency (l
(cont') Equations for Simple Systems
FREQUENCY EQUATION
[
CASE
~
~
K,-X.w" -k, 0 0 s.
~r n
en
-'<. IC,.,.:~•n:t, ~ea.n: ~n,·
8 0 n: lea. 0 =0 0
"'n:-r....• .....
0 n;l<. 0
04'
~
J,
Pr.--0 r. "·~..-lj ...~
I~ r, i
T, ~ t.·I,w' -Jr. 0 0 0 ~
9 -~ lr,•~li:a•Alr. ~k& 0 "1kt 8
0 <ke. ~(ka+(s•lJ) •ltl\~ 0 :0
0 0 - t, ~ n:cr,-r,..•) 0 m:
0
[ft~ ::s
0 n4tf 0 0 ~(11-I.c.') en
r .. '<
en
g
k,-r,w• "*• 0 0 0 ~ ~
10 [}1}-IF{]r, •ta IC,+k&•Jiw' •b 0 0
0 ·ta n;t, ~te. +ta Kanl t.tlt 0
n.~r., 0 0 t<,nJ n;kJ-q;I3 w1 0
r, I.z . 0 0 "'~· 0 "'~.Z.,.,'c.•
J;~ r.
r~ ~ l<z
lei·liw" -lc· 0 0 0
11 I.e~ ilu~O -k· lee•ta.• ts·l'all - ka 0 -IC's
I<~ Ks 0 -t. ka•tcs·.r,..• -t~ 0 =0
0 0 -k• ta•a:..-:r....• -1:"4
1<4 IS~ 0 -ks 0 •k1 t'4 •t',-JS..:.~
Ei' Is
t!, ~
\0
1.11
CHAPTER 11.1
H.lrretier
University of Kusel, Kusel, FRG
ABSTRACT
1. INTRODUCTION
One of the important components of steam and qas turbines are the
blades. They re attached to the rotor surface to conv~rL
the fluid energy of the gas or the steam into rotational energy
of the rotor. A typical bladed rotor is shown in Fig. 1
In the following paper the vibrational behaviour of turbine
blades is considered whereas it is assumed that the blade
vibration is uncoupled from the disc and the shaft. This assump-
tion is valid in the most real cases. Otherwise, the coupled
vibration of the blade and the disc or the complete system of the
turbine rotor must be considered.
398 H. Irretier
a2 a2w a2w
- ~ (Ely ~) pA (3.1b)
dX~ axL ~
of the classical Bernoulli-Euler beam theory. However, these
402 H. Irretier
a a a * a 2 a2
ax [ KGA(a~ + Yz)] + ax(nx a~)+pAA (vcosa-wsina)cosa= ~
at
2 2
a a ayz ay a Yz ay
KGA(~xv + Yz) a~x(Eiz ~X + EI yz _l)=-p(I z~- Iyz~)
a a ax at~ at~
2
! = ~I + n .!$a (3.6)
+
2
2 22 I
ly = ~ (n + r,: ) dA - (3 .11)
-·C = n-
.!l K, (3 .21)
whPrP
406 H. Irretier
If the refined effects for the beam-like blade models are taken
into account, only approximate solutions are available. The
classical procedures are the Rayleigh-Ritz and Galerkin-method.
They are described in /4/ and /15/. In the last twenty years, how-
ever, computer-aided numerical methods are preferred to calculate
the natural frequency and mode shapes of beam-like blades. One
classical procedure is the method of transfer matrices while more
modern techniques are the method of finite elements and a direct
numerical integration of the differential equations by a Runge-
Kutta-method and a shooting technique. A large amount of refer-
ences to the first method is given in /3/ while the second
technique was used in /5/ and /21/.
The suitable method for the calculation of the free vibration of
shell-like and three-dimensional blades is today the method of
finite elements which was described in the chapters 3.2 and 3.3.
For all three types of blades the technique of modal transforma-
tion is the most applied method to calculate forced response of
blades due to the various types of excitation. Details on this
have been described jn the chapters 3.1.4, 3.2 and 3.3. For
details of the application and numerical problems is referred to
the numerous literature in this field.
5. RESULTS AND PARAt·1ETRIC STUDIES
General statements concerning the dependence of the vibrational
behaviour of turbine blades on the constructive parameters are
only possible for the eigenfrequencies of beam-like blades. For
their forced vibrational behaviour and in particular for shell-
like and three-dimensional blades the connections are so compli-
cated that only some typical results can be discussed here.
5.1 BEAM-LIKE BLADES
Besides the mode shapes of the bending and torsional vibrations
of blades, the influence of the constructive parameters on their
eigenfrequencies is considered in particular. Moreover, two re-
cent results concerning the forced vibration response of the beam-
like blades subjected to nozzle excitation and partial admission,
respectively, are presented.
5.1.1 FREE BENDING VIBRATIONS
5.1.1.1 MODE SHAPES
As pointed out in the chapter 3.1.1 the most simple model for the
free bending vibrations of a turbine blade with constant cross-
section, without twisting and rotation and neglecting the shear
influence and the rotary inertia effect are the differential
equations (3.1a) and (3.1b) for an Euler-Bernoulli-beam. The
corresponding mode shapes are those of a cantilever and consist of
harmonic and hyperbolic functions /10/.
If the refined effects are included in the consideration, the
blade mode shapes change more or less considerably. The influence
410 H. Irretier
x, ,' x2 2.
75.0
).I
l.,
• -=~:?
50.0 tn E•.B,
L l.•4;:-Vf
l?n Vf'
~
.!:-.so n•• "T,, T
-
1,0
r1' .J:•1
25.11- l.z I
(~-·-
l.,
~·o"
•lief
,:o2
,., -
.o
.o 1,0 2.0 3.0 4.0 5.0
.9
.8
----0--
fixed at the root and free at the end, numerical results show
that this effect can be often neglected for the first eigenfre-
quency but is of importance for the higher ones.
>.,
1.0 .----------====::;·====-=---,
>.~-
>.,~-.~ // ~ ·---
.~
/
0
.8
.7
.6
---X
twisting is completely neg-
lectable for the first
~X eigenfrequency while the
iyo =0.2 second one decreases and
'zo the third one increases
20r-------r-------+-----~ considerably. This result
is a consequence of the in-
creasing coupling between
the two bending vibrations
in y- and z-direction with
10 increasing twisting angle.
l
~-------~-------~-------~
01 I I I
Fig. 11: Influence of twisting on
the eigenfrequencies of the bending
vibrations of beam-like blades
5.1.2 FREE TORSIONAL VIBRATIONS
5.1.2.1 MODE SHAPES
The mode shapes of a cantilever with constant cross-section per-
forming free torsional vibrations are described by harmonic func-
tions which are the general solutions of equation (3.9). For
turbomachinery blades the torsional mode shapes are affected by
tapering, twisting and the other refined influences in a compli-
cated manner. A typical mode is shown by the holo~ram in Fig. 12
which 1s taken from
/25/. We notice that
the lines of equal
displacements are no
more parallel to the
edge of the blade.
Fig. 12: Mode shape of the first torsiona
vibration of a twisted turbine blade /25/ If, in addition, the
torsional vibrations
are coupled to the bending vibrations as described in chapter
3.1.3, the mode shapes are even more complicated and no general
description of the influence of the blade parameters on its mode
shapes is possible.
5.1.2.2 EIGENFREQUENCIES
As for the bending vibrations also the eigenfrequencies of the
torsional vibrations of turbine blades are considerably influ-
enced by the effects of twisting, tapering and rotation. As an
example, Fi~. 13 shows the influence of the rotation on the
fundamental eigenfrequency of the torsional vibration of a beam-
like blade. We notice again a parabolic increase of the eigen-
frequency. However, in comparison to the bending vibrations the
increase of the torsional vibration eigenfrequencies is much less
414 H. Irretier
2.0
A.,
tQ !yo =0,25
ltJ
1.9 lzo
E:ft
L
.!!2
lpa
=0.243
R L
1. 8
1. 7
1.&
1,5
.o .1 .2 .3 .4
6000
sooo
t:O x,
\ Po
\
0.9L-
blade force
Fig. 16: Rotating blade and load due to two arcs of partial
admission /26/
quency of f = 10 Hz up to a rotational frequency of 60 Hz as
shown in Ca~pbell-diagram in Fig. 17 where the resonance points
during this run-up are indicated for the first two eigenfrequen-
cies of the blade. The stress response of the blade at its root
f0 [Hz]
60 -----------=""r------
I
0 + 10
-10.0 -5.0
I I
I I i . z-118 1~
600.0 -:r-~.J..+t.....:........:~~-7----~iF-------:::::=o"'
6
f 1[Hz)
400.0
300.0
2oo.o l-A~~~--
2
1oo.o~~~;;;;~~:::;~:i::::~::::::~
.0
.0 10. 25.0 50.0 60. 75.0 fo[Hz) 100.0
~1 = .01 ~2 =.01
t, 2 t0 f2
6 f0
10 f 0
14 t 0
\
l~ ~
~ f' I:FEIOI4Hz EX 1054Hz
~
2; FE 1713 Hz EX 1717 Hz
~
\ ~
lOO
~so
fiHzl
200
--
240 3n
:!20
:oo
180
luO 2n
140
1~0
100
80 n
bO
20
Q
0 2000 lOOO n I rpml sooo
7. REFERENCES
I 1I BLOHf·1 & VOSS-Prospekt 11 Schaufeln fUr Dampfturbinen .. ,
Blohm &Vo6, Hamburg, 1976
I 21 Rao, J.S.: Turbine Blade Excitation and Vibration.
Shock and Vibration Digest 9 (1979) 15 - 23
I 31 Rao, J.S.: Turbomachine Blade Vibration.
Shock and Vibration Digest 12 (1980) 19 - 26
I 4/ Bohm, J.: Theoretische und experimentelle Parameterstudien
an schwingenden Turbinenschaufeln im Fliehkraftfeld.
Fortschr.-Ber. VDI-Zeitschrift 11, 29 (1979) 1 - 136
I 51 Montoya, J.G.: Gekoppelte Biege- und Torsionsschwingungen
einer stark verwundenen rotierenden Schaufel. BBC Brown
Boveri Mitteilungen 53 (1966) 216 - 230
I 6/ Rao, J.S.; Carnegie, W.: Solution of the Equations of
Coupled Bending-Bending-Torsion Vibrations of Turbine Blades
by the Method of Ritz-Galerkin. Int. J. of Mech. Sci. 12
(1970) 875 - 882
I 71 Carnegie, W.: Vibrations of Pre-twisted Cantilever Blading.
Proc. of the Instn. Mech. Engrs. 173, 12 (1959) 343 - 374
I 81 Carnegie, W.: Vibrations of Pre-twisted Cantilever Blading
Allowing for Rotary Inertia and Shear Deflection. J. of
Mech. Engng. Sciences 6, 2 (1964) 105 - 109
I 9/ Fu, C.C.: Computer Analysis of a Rotating Axial Turbo-
machine Blade in Coupled Bending-Bending-Torsion Vibrations.
Int. J. Num. Meth. Engng. 8 (1974)
/10/ r~eirovitch, L.: Analytical Methods in Vibrations. Collier-
Macmillan Ltd.; London, 1967
/11/ Bishop, R.E.C.; Gladwell, G.M.L.; Michaelson, S.: The Matrix
Analysis of Vibration. Cambridge University Press, Cambridge/
London/New York/Melbourne, 1965
/121 Hurty, W.C.; Rubinstein, M.F.: Dynamics of Structures.
Prentice-Hall Inc., Englewood Cliffs N.J. (USA) 1964
/1l' Lehsa, A.W.: Vibrations of Turbine Engine Blades by Shell
Analysis. Shock and Vibration Digest 12, 11 (1980) 3 - 10
/14/ Vogt, H.J.: Zur Berechnung der Eigenfrequenzen und Eigen-
formen von Schaufeln thermischer Turbomaschinen mit Hilfe
eines gekrUmmten finiten Schalenelementes. Dissertation,
TU Hannover (1974) 1 - 87
Free and Forced Vibrations 421
N.F. Rieger
Stress Technology Incorporated, Rochester, New York, USA
ABSTRACI'
11.2.1 Introduction
11.2.2 Notation
tan 8
sonic velocity in gas ("yg RT)
constant in reference [19] see below
constant in Whitehead's analysis [13]
axial distance between blade rows
426 N.F. Rieger
2
DM
dimensionless force (defined in te:~:t)
dimensionless force (defined in te:~:t)
constant in reference [19] see below
F constant in reference [19] see below
Fr Froude Number
G Glauert e:~:pansion (defined in te:~:t)
hm constant in reference [17]. Also water depth
i
J(a)
fl
Bessel function of first kind
K(a) modified Bessel function of second kind
k constant in reference [8]
1 disturbance in wavelength 2nVs /1
L non-steady a:~:ial force
Llt non-steady tangential force
Ly steady lift
Lo unsteady lift
M Mach number
N rpm
p pressure
p constant in reference [19] see below
q blade velocity due to vibration
r blade mean radius
a• gas constant
8(0') Sears' function (defined in te:~:t)
s blade pitch
t time
T temperature
T' (a) Holmes function (defined in te:~:t)
T(a) Horlock function (defined in te:~:t)
u non-steady a:~:ial velocity of gas
u mean a:~:ial velocity of gas
v non-steady tangential velocity of gas
v mean tangential velocity of gas
v. blade tangential velocity
general non-steady gas velocity
w
w mean gas velocity relative to blade
:X(a) function in reference [11]
lt a:~:ial direction
y tangential direction
Y(a) function in reference [11]
a angle of attack
IS mean flow angle relative to axial direction
r specific heat ratio (Cp/Cv)
8 constant in Whitehead's analysis [13]
e stagger angle
eLE lead edge blade surface angle relative to a:~:ial
direction
Flow Path Excitation Mechanisms 427
Constants:
D = 4A2 -"fi 2
F = [ (2A2+ )bA-D) + i (2A3-'¥1]
Bs = (l+b s ) I (1-b s )
a = arc tan U/Vs
h1 = '2 - 81
428 N.F. Rieger
2 2 2 h3
g1 h1 g1 2g1 h1 g1h1
h5- -+ -+ -+ + - - + ..!
2 3 2 3 3 15
2 h3
g1h1 g1h1 1
h6- - - + - - +
2 2 8
h3
2 1
h7 - '1 + g1h1 +
3
2 2
h8 - 12 - 11
Information Sources;
Vortex Theories:
where the 1 (1) are Beaael fuactioaa of the firat kiad. aa4
C(a) is thentbeodoraea fuactioa. The Theodora•• fuaotloa la
aiven by:
The results shown in Fiaure 6 are for the first two harmonics
(a = 1.2) of the lift ratio (unsteady lift/steady lift)
plotted aaainst the axial spacing ratio (b/c) and the pitch
ratio (Sr/Ss>• The dashed curves correspond to an elliptical
steady load distribution on both the rotor and stator blades.
while the solid curves are for an elliptical distribution on
the stator blades. and a flat-plate distribution on the rotor
blades. The authors concluded that 'the non-steady part of
the lift may be as large as 18.. of the steady lift and
therefore may be of practical importance.' This low value is
a feature of the simple model chosen. It is by no means an
upper liait for practical blade conditions.
CD
X
i l1(t - - )
u = u0 e U (7)
where J (a) are Bessel functions of the first kind and K are
modifiel Bessel functions of the second kind. The resulfs of
this analysis may be combined with the results of the first
Kemp-Sears analysis to obtain the total circulation-induced
lift fluctuation on a blade due to a aeneral periodic
disturbance of the free stream flow. This disturbance w is
resolved into components perpendicular to the chord wt and
. wc as shown in Fiaure 7. In aeneral,
parallel to the chord
the velocity of the disturbance can be expressed as a Fourier
series:
L=2x p Usee' w
on
ei~t[S(a) sin,-aT(a) cos' ]
a (11)
where U is the axial velocity, ' is the aean flow anal•• " is
excitins frequency, and a is reduced frequency. 'Dle abo••
expression for fluctuatina lift applies to flat plate
airfoils. To account for camber Horlock incla6ed the Bolaes
function, T' (a), which was derived for an airfoil of para-
432 N.F. Rieger
~- 0 0
2
p U 1 ei1 tv S(a)+u [aT(a)+( Ymax)T'(a)]
1
(12)
L
_L_
pcv2
... t[2 c-
• v
2uoUs
- -- h
1
sin w] sin1t
,v2c 1
(15)
4u U u h4 }
+{[ y2° h 3 (1-cos M) + - 0- - sin w] cos~t
.., sec' Vsec'
A a imilar ezpreaaion "141' obtained for the unsteady azial
component of lift, Ly /pc~.
L
(16)
.. neD
s
(17)
B2acose
+ i( 2 -F)}
v~ tane
This result applies for the case of a seneral stasser ansle
e. where:
water table test showed the desree to which tho water table
load transient resembled the load transient obtained from
actual bladins tests. Tho staso parameters and typical
results are shown in Fisure 4. The followins similarities
can be observed:
11.2.5 References
29) 1ohnson. R. H., 'The Hydraulic Analogy and its Use with
Time Varying Flows,' GE Research Laboratory Report, 64-
RL-(3755 C), Schenectady, NY, August 1964.
FReouENCY <Hz>
Inlet Nozzles
~p;pppj~
-- (( (( (( (( (( (( (( (( --- Moving
Blades
Tftt Conditions
Type of staae Impulse
Number of IBOVIng b'-des 1Zl
Number of .,-rtlal admission arcs 3
Number of nozzlft In arcs 10,10,12 {corresponds to 7& nozzle~
In complete drcle}
Stage pressure ratio 1.58
Stage velocity ratio o.u
Number of Instrumented blades 1
lnst. Force Direction of force measurement Tangential
Avg. Force
lnst. Force
20 Enter Leave Avg. FOrce NPF
Ripple
Blade Pitch
I I I I I I I I I I I
0 2 • & I 10 12 1• 0 2 • & I 10 12 1•
. 20r---,----T----~---,
m=1
o L...!.1:::::::5:~:Lmt:Z.J
0 0.2 0.11 0.6 0.8
b'/c r
Ratio of Unsteady Lift to Steady Lift vs. Spacing Ratio (Stator
Pitch Equals Rotor Pitch).
Figure 6. Results From Kemp-Sears Analysis
446 N.F. Rieger
~"'LA(_ ''-LJ_;/'~ x.
VELOCITY PROFILE
w • w e1v(t - &>
0
Upuru• Oowni\reaa
~ y • weiwt
ue 1 "''
~
U•
y • veiwt
U• ue iwt
At\UI\Or Dht
f1oVJNI 11.\DII
JIG
--
cr•J.OC'"'ltehea412ll 0.:1
0.5 O.J
..,...- --
o.a
o.y /"""' cr•O
/
/
O.l
O . J , _ Neuar-·llodocki:IJI
I
cr•O /
-o.4 -0.) -o.J -0.1 0.1
-. . . _ >
-.
...
'
. . _....L>- .
....... 0=0.11
-.. _.._
c:
u
...
..
c:
1. 0
-- O•O.l
.... ....
-- .. --
--- ..... -------
---
--,..
c3
::J
0.1
-----P- ---------. -
0•0.1
-
i
••
c:
0.6
Oc0.01
-
...
~
--
---
Semiactuator Disk Analysis
Ref. II
-
0
-•
•
'9
cDl
0.111 f- Stagger Angle II!'
Flat-Plate Airfoils at Zero
Incidence
::E 0.2 -
.
0.6
I
o•• 1.0
.
0 0.2
Space-Chord Ratio S /C
L a "'ICOW c elvt
0 Ill · - -
SC•t.•W
Incidence \ ;::
'' 0
o. 3 ~
IC
SC=- \ ~ It
1 I I I o. 2
' I I I
SC=O·~I'
\ I I
I
''· '~ :: o. 1
...
', ~·
\'I
~------L-------~------~---~------~------_.------..JO
-o.6 -o. s -o. 11 -o. 3 -o. 2 -0.1 0
IL/Lokeal
~ l
....
:.
1: •
.
0
c.
5
Per-rev
Harmonics
Non-Steaay N.P.F.
Blade Force
Frequency, Hz.
,--·--·--···--···-··
i 1$1 "'""·
···--...................
,CRt( IIAlJII l
I:DW·!.ll~1
:
I vs. !
! Yti.C:J n RATJII j
i
,! nu ~
j D11lt1JOII: h"C••u•l
j rcusuat tiiTJOS: ::::~!
! NOtl\. IDIItlh.SIU, AI.OJ ll.UI a:ui
i HOIIII. Y(l.RATJD: .)000 I
: !
~ ;_""""-"'"-"""'"--'-""""""""'-l
N.F. Rieaer
Stress Technoloay Incorporated, Rochester, New York, USA
ABSTRACT
11.3.1 Introdpction
Fatisue:
Corrosion:
Corrosion Fatigue:
Erosion:
at = am + a a cos wt
a
mo = aco +abo
= P/A + .MC/1
da A [AK]n
-=
dN U-R> 0 • 5
a ao
Aa = 2aao
a
A.a = 2Ktaao
11.3.7 Conclusions
11.3.8 References
da
dN
INCH PER
CYCLE
° 10
1
-6~~+-~~~~~~
2 4 6 c 10?
.6.K KSiffN'
llil1
First Tangential
\CCC&~
UUt
First Axial
Redesign
Vane
section
Original
Design Crack Overhang
origin removed
in redesign
(a)
(b)
Figure 8 Second Type Tangential Group Mode (a)
Second Type Tangential with Tiewire (b)
6
znd AXIAL
2
34 NOZZLES
:J
5
N
Q
-4
K
.2
~
rlt"W (
SL A ~Jl
I
.r .·,
...
I _L_:
rLLL.£ '////, rr//L rLLLL r..tlliLI.J
I I
/:'
.....
I
Goo()
'
-
, I
i
I '
..; 44100 I 1
~-
:r: --:
>, li'. 'fJ. r ~
u
A'u'
I _I
c:
Q)
::> I !
aQ)
'LL.L
.... 2ooo
--1
iL I JAX'
LL
1/ A!
0 T1.
0 1000 ZJXIO -
Speed RPM. "* &1100
N.F. Rieger
Stress Technology Incorporated, Rochester, New York, USA
ABSTRACT
11.4.1 Introduction
11.4.2 Notation
D Damage
Aa Stress range
Aa Strain range
am Mean stress
aT Total stress
am Mean strain
aT Total strain
11.4.3 Definitions
Engineering Stress
True Stress
a = P/A = a0 (1 + a0 )
Engineering Strain
a0 = a0 IE = A1/1
True Strain
Procedure for Component Life Estimation 489
Ae (Aa/2K') 1 /n'
p
Ae e = (Aa/2E)
R-ratio
R = amin I amax
Notch Sensitivity
q = (K f - 1)/(K t - 1)
490 N.F.Rieger
Fatigue Factor
a = constant, r = radius
Neuber's Rule
2
(ae = aEa)
1/n'
.As .Aa .Aa
= + (ii,) (1)
2 2E
(2)
a
n
As a ,
-= f (2Nf)b + , (2N )c (3)
2 E 8 f f
As
(3a)
2
a
As af'2 (2N >2b + a ' s ' (2N )b+c (4)
max = E f f f f
2
(S)
(6)
The local stresses and strains are first computed from the
loading history (keyway elastic stress o vs. time) using
equations 1 and 3. Typical results are sum:arized in Table 2
as local stress vs. strain data. For instance, a peak stress
o • of 200 ksi (point b) during the manufacturer's overload
t:st results in local plastic deformation at the keyway, and
a local stress o of 123.7 ksi. This value of o is computed
from Neuber's rule, equations 1 and 5, using Newton-Raphson
iteration, along with the local strain 8. These computations
may be summarized as follows:
200 2
08 = = 1.356
29,500
08 = o(-) +
0
0 (~,) 1/n'
E
o = 123.7 ksi
08 = 1.356
8 = L.lli = L.lli-
123 • 7 - 0.011 in./in.
0
Both the loading history and the local stress vs. local-
strain diagram show both low and high cycle loads. Low cycle
loading is represented by the hysteresis loop, points c
through g of Figure 3, whereas the high cycle loading
consists of fluctuating stress levels (between points e and
f). Although not detailed here, the total strain amplitude
aT and the mean stress om are respectively computed for the
low cycle loading to be 0.004520 in./in. and 18.3 ksi, and
for the high cycle loading to be 0.00539 in./in. and 71.1
ksi.
Procedure for Component Life Estimation 499
The fatigue damage due to both high and low cycle loading is
determined by solving equation 3 numerically. The results of
this computation are listed in Table 3, for fatigue damage
data. Table 3 also lists the damage per year which then can
be summed to determine the total fatigue damage which occurs
each year due to both types of loading. Thus:
1 n1
-= ---+ ~
N.r Nf(1) Nf(2)
= 0.1090
N.r = 9.17 yr
The calculated results showed that the R1 mode and T22 mode
were responsible for the highest dynamic stress in the upper
hook-notch region. The calculated rms total dynamic stress
Procedure for Component Life Estimation 501
was 3335 psi. The resonant stresses in mode R1 and mode T22
were 5652 psi and 5387 psi respectively.
The stress due to this overspeed was 66.1 ksi. Duriug steady
load operation the stress-strain condition was at point 4.
The steady stress due to normal speed operation was 55.9
ksi. When the vibratory stress reached the positive peak the
total stress we~t to point 5. When the stress reached a
minimum value, the stress dropped to point 6. The stress-
strain condition went back to point 2 when the unit was
brought to zero rpm.
11.4.9 Conclusions
11.4.10 References
,•..
li
..=' ···:r
j
§
.
<il ...
···l
Factory
a Overs peed Safety Safety Safety
Trip Trip Trip
Stress
Process
Load
Time
100
'iii
~ 50 Hysteresis
b loop
"'"'
~
iii
<ii
0
0
...J
Residual
- 50
~ 1- One Block
f :::~·····~
(o) Incremento! Spectrum
Tom•
rOksi
0 .004
Typical Vane
F Gas Pressure Passing
Dynamic Components Frequency
Force
Frequency
Figure 5 Excitation Spectrum
a'
•
t - - - - - - - - r l"
,)'::
I I
a'R ~ I I
Ill
Ill ~-----~\, I I
C1l
s..
+"
I\ j \
V')
/ \ / \ - tot• I dyNIRic ' ' ' " " "
u I \ /
.,_
E / v· ClyN•'c ruponse fOf" -ode 1
rtl
c: \
>, \
0
\
''
' ....... .............. _
,, ·Frequency f
€ , Strain
-i 1 ;L 2
31
3
~ 5
]4
6
7
1 18
9f ~ 10
ltzcycle
1cycle
112 cycle
1 cycle
I cycle
112 cycle
::200
.,.,
e 1so
iii
u
.
·; ; 100
'ii
c
Time
10~-----------------------------------------,
RESIDUAL STRESS
a AND STRAIN
1
NOMINAL
STRESS
....
0. 0010 l:---.,..--5tic-$t...lo...__;~.,.-~::::::;,._~-cl----l
/TIItudel
Plutlc Str•ln
Amplit<.de
a o.o o.o
b 123.74 0,01096
c,g -69,18 0.00393
d 105.6 0.01297
e 61.1 0.01229
£ 81.1 0.01161
....
MAXIMUM STIIESSES IN ILADE AND DISK ATTACHMENTS
-
lower c.... nuu pol 1SIUO pol
'-k c ... 2 103770 pol unu pol
512 N.F. Rieger
A1 863 759
R1 990 888
6 Tl 1052 919
BLADES u 161l4 ---
s 3119 ---
A1 850 763
R1 901 816
12 Tl 1002 922
BLADES u 1041l 931
s 1320 121l9
0'
....
(")
0
Table 8 Fatigue Damage Evaluation for Six-Blade Group 3
"'0
0
::s
11>
::s
....
Mode Al R1 Tl A3 T22 I
t""'
;:.
Natural Frequencies (Hz) 863 990 1052 3849 4481 tTl
"'....§'
Stress Amplitude 3.45 5.65 3.70 2.37 5.39
~
( ksi) ....
Load
o·
::s
Hfgh Cycle Fatigue 1.21 X 10•J 1. 93 X 10· 1 2. 97 X 10-J 1.26 X 10· 4 5,45 X 10-l
Case
I Damage 1/Y f
U1
....
~
CHAPTER 11.5
N.F. Rieger
Stress Technology Incorporated, Rochester, New York, USA
ABSTRACT
ll.S.l Iptroduction
Details of the test rig used, and of the test conditions and
results achieved are given in the following sections.
30,000 psi. The results showed that the material log. dec.
values were highest with high dynamic stresses. A maximum
test log. dec. value of 0.033 was obtained. Damping
dec rea sed with increase in centrifugal load, and increased
with frequency of excitation.
Material Damping:
u.s .3 Apparatus
The blade pairs were mounted in the rig in the usual manner
and durina testing were pre-loaded usina hydraulic load
cells. All blade and attachment components were contained
within the environmental chamber durina the test proaram.
Results for 6 inch blades with axial entry fir tree roots are
shown in Figure 11, and these again show an increase in log.
dec. with blade tip displacement. Data is given for tangen-
tial vlbrations and for axial vibrations, in air. The
damping relationship is seen to be almost linear in these
instances, though the vibration test amplitudes were
relatively small. Later tests established that damping
continues to increase with higher vibration amplitudes.
Increased centrifugal load causes the log. dec. values to
decrease, as was found by Lazan [17] for damping which is
related to material hysteresis.
the first time). and old roots (which had been used in
previous tests). Generally speaking there is little
difference between the values obtained under vacuum and in
air. Damping data with used roots in air gave somewhat lower
values than air tests with new roots. These tests were
conducted using a 25 inch blade loaded lengthwise with a
force of 13,400 lbs. which is smaller than typical operating
attachment loads carried by such blades.
11.5.6 Copclusions
11.5.7 References
..·
-
··~-.....~· .~-!r•• ----
•:f.,· .. _
- :;i-_'·
466l Digital
Axlol Transducer Interactive
Transducer Power S..pply Ploller
• - hu
Filter
,.,. R•O T,_ r-D
SFMCtrulll An.lyzerl .
crCJ---CJ-o
MP 100
Tangential Transducer O.t. Ceneral
Transducer Power Supply Computer
c
".
E
¥
0
·eu
z
..
·;:
o.oo
0
0
0
0
"'
N
0
0
0
0"'
..
0
"'
0
0
0
0
.:
0
Blade Tip Displacement (Mils)
Figu re 7 Loga rithm ic Decr emen t vs. Tip Disp lacem ents
and Cent rifug al Load . Tang entia l Mode 2"
sure Blad e
Strad dle Mount. Tang entia l Entry High Pres
o.so
[7"u.
~E Centrif ugal l.oad
.. 0. 30
~
0
'/
u
e .o.2
=
.
·;:
j ~30001b.
o. 1ov.....,.....,..
_soo o lb.
t:- - roo lb.
o.o0
0
0
0
0
0
"'
N
0 ..
0
0
"'
.."'
0
0
0
0
Figu re 8 Loga rithm ic Decr emen t vs. Tip Disp lacem ents
and Cent rifug al Load (Axi al). 2" Strad dle
Mount. Tang entia l Entry High Pres sure Blad e
Damping Properties of Steam Turbine Blades 535
AXIAL
0.~0
o.oo
c.
E TANGENTIAL
~ 0 . 029
~
c
"
.E
.8'
~ 0.021
\:
..J
O.OIQ
0. 000 L---...L...---!----:!::-----:!
c
0
::!
.,;
c 0
:::
0
.,;
0
c
..,c
.;
.
0
.,;
0.060 I I 1- '•
c.. o.tooJ-,..;.-,--r--r:::;=;r~
E
I 0.0, ..
~
Q
u
..v
p e
£
:i }
~ O. Ol
. 6 I 10 12
'
Bl•de Tip Displacement (In x to•')
Blede Tip Dlspl•cement (In x to•'l
AXIAL
o.ou
c 0.0)0
I...
~ \. 22, 00 lb.
... 0.020
'E \
:; 1l.-r lb.
'l:
j 0.010
'-- !--"' t.oo:o lb.
0.000
0 0.25
t
0.50 0.75 1.00
Blede Tip Dlaplac-t (Mile)
0.0)0 ---.,....---r----,~---,
c..
..
E
0.020
Q
~
.!l
E
=
.
·;:
g'
..J 0.01U
0 .ooo __________...____
0 0.003 0.006 0.009 0.012
..... Axlol
••••
I ..-- !!5'".....
~-
! ,.....-:
t.tzt
E;:;.y
Ii:
J
t.tlt
.....
I
.,;
i".: i.
- · - T i p D l . , . _ t (IN,)
0 • 03 or-----~------~--M-o-d-if-ie-d~M-id_d_l_e--,
Hook
I
c
II
0.020~----~----t--.rl~~---t-------1
E
f
~
0
..
u
"j§
.t:
·5en 0 • 010 ~-------1.1-
• •-.-.-:-1~34~0:-::0~I:-b-::C:-e-n7
tr-:;i7fu~g::a:;l1L-::oa::-::~rl
.3 ----- 9000 lb Centrifugal Load
- 5000 lb Centrifugal Load
o. 100 ,...,..--,~-.,.--,--....--..,........
..i 0. 080
1 Axial Deflection at 0. 0007"
.] 1 Standard Deviation
.
J!
E
.t: 0.040
} 0.020
o.ooo
0 3000 5000 7000 9000 11000 13000
Centrifuqal Load. Lbs.
-
....
0.070
---
0.060
I
• •
,.. ••.... ·~·
..~·
Logarithmic 0.050 I'
Decrement
6 o.on
~.
, ••
0.030 • • •
0.020
0.010
o.ooo
250 750 1250 1750 2250
0
Stress (PSI)
LP 25• Curved axial o. 02o-o. 025 0. 02o-0. 030 Adjacent root study
entry
MAGNETIC BEARINGS
G. Sebweltzer
Institute for Meebanies, ETH Zurieb, Switzerland
ABSTRACT
Magnetic bearings have some distinct advantages. They do not
generate wear and they do not need 1ubri cation. These features
make them at tractive for vacuum applications. And their dynamic
behaviour can be adjusted in a wide range, which allows active
vibration damping and control.
This chapter presents the state of the art for the design of an
electromagnetic bearing system. It introduces first the main
elements and then discusses control and system aspects. Models
for describing an elastic rotor and its active vibration control
are included. The characteristics and the losses of such a
suspension system are detailed.
Several applications are demonstrated, and future trends are
indicated.
CONTENTS
1. Introduction
2. Functional principle
3. Design goals
4. Elements of the magnetic bearing system
4.1 Model for the rotor
4.2 Sensor, controller, amplifier
4.3 Magnetic actuator
5. System aspects
5.1 Control of rigid and elastic rotors
5.2 Characteristics of the magnetic bearing system
5.3 Losses
6. Applications
7. Conclusions
8. References
544 G. Schweitzer
1. INTRODUCTION
Magnetic bearings can support a rotor in such a way that it
levitates freely without any contact. Furthermore the dynamics of
this suspension can be easily adjusted in a wide range for
various applications. These two main properties already make the
magnetic bearing a very attractive device for solving the
classical bearing problem. On the other side the magnetic bearing
is complex, expensive, usually not readily available from the
shelf, and up to now only used for some advanced machinery. In
the following the state of the art is presented so that future
trends can be derived.
Magnetic forces are generated either by permanent magnets,
electrodynamically or electromagnetically. In the constant field
of permanent magnets, however, a ferromagnetic body cannot hover
in a stable way /BR 39/, and electrodynamic forces are usually
too small or still too difficult to generate to be of actual
technical interest. They are used where small forces are
sufficient, for example in space applications for the support of
flywheels or for a micro-g-platform in a near zero-g-environment.
Or they are used where the high currents necessary for large
forces are generated by means of cryogenics as in prototypes for
an electrodynamically levitated high-speed vehicle.
It is the electromagnetic force that is used most efficiently.
For the l~vitation of guided vehicles a technology of its own has
developed /GO 84/, that basically of course has some connection
to magnetic bearings for rotors, too. Rotors have been supported
magnetically at first for physical experiments. Spectacular
2.10E7 rpm have been reached while testing the strength of small
steel balls under a centrifugal field of l.OE8 g /BY 46/. Since
then the electromagnetic rotor bearing has been applied to solve
a number of different technical problems, and therefore the
construction and the properties of the bearings differ
remarkably. Numerous patents in this field have been issued or
are pending. And a few types of bearings are already available
commercially for use in machine tools and in turbomachinery /HB
85/. More details on applications will be given in chapter 6. Now
at first the functional principle, some design goals, and the
elements of the bearing system will be discussed.
2. FUNCTIONAL PRINCIPLE
Fig. shows a hovering shaft S. Any deviation from this
reference position is measured by a sensor, the measured signal
is transformed into a control signal according to a suitable
control law. The control signal is amplified and fed as a control
current to the coils of the e 1ectromagnet, the actuator within
this control loop. When the shaft for example starts to fall down
it produces a measuring signal which leads to an increase in the
Magnetic Bearings 545
3. DESIGN GOALS
Primary design goals have been to support a rigid rotor, because
the rigidity of the rotor facilitates the control design
essentially. Magnetic bearings are quite capable of supporting
even a heavy rotor. The considerable freedom, however, in
assigning dynamic characteristics to the actuator can be used not
Magnetic Bearings 547
;- 1 _
· Sensor plane "d"
+--G_~_- _!_a--I-c .. X
Bearing pla{le "a"
cb2
I
Fig. 4: Rigid rotor model /BL 84/
-uf = [ f ax , fb x , f ay , fb y] T
(4.2)
The harmonic disturbances by unbalance are introduced by the
vector
I
stiffness are
0 0: 1 0
p = I ~ _o- ~ ~- ~ -~- ' s = 0 (4.4)
Z -1 0 0 0I
0 0: 0 0
I
Magnetic Bearings 549
a b :0 0 -I I
1 1 I0 0 Oyz mexz
~f --- -'----
I y_f = ~2 (4.5)
0 0 1a b I xz I
0 0 II 1 1 me Oyz
-x =[ -s i T] T
z T -s ( 4. 7)
l [ l [ l
E I- 1
, B=
-
-
I !1.-0 1 ~f
V=
-
-
I !1.-0 1 Y..f
(4.8)
where g is a n x 1 vector of generalized displacements, £ is the
vector of generalized forces, and the A. are structural matrices,
characterizing inertia and elasticity -Js well as gyroscopic and
nonconservative properties. Equ. (4.8) may be transformed into
the complex frequency domain and expressed by
From the so 1uti on of the ei genva 1ue prob 1em we obtain a set of
e~genvalues w1,._.. w0 and the c?rresponding real normalized
e1genvectors ~k w1th the modal matr1x ~ =[~ 1 •••• ~k····~n]·
Then the modal expansion of g
with
M= ~T Me ~ = diag (mk)
! = ~T !a ~ = diag (kk)
Magnetic Bearings 551
T
2 + K) -1 UT = I:n ~k ~k
H(s)
- = --
U(M s - - k=1 k 2 2
(4. 14}
m (s + wk )
T
! = IA £ • 1 0 = I 0 s. 1v = Iv s. 1 T = [ 1o·
T 1vT] (4. 15)
=~
lo !!m( s 2 !:\n
+ ~) -1 !!mT !A
.Y_(s) !!(s) ~ =[ 2 -1 T (4.18)
s !v !!m<s !:\n + ~) !!m !A
Here the trans f-er matrix H is a generalized dynami ca 1
flexibility matrix, since it rifates a force excitation (actuator
forces) to nodal displacements and velocities (sensor signals).
By introducing the state vector
T
~ =[ ~.
T • T]
~
an equivalent state space representation is obtained
= -e
C
l [l
V X
-e -e (4. 19)
A=- -
E ,B= .0- , c =[ !o0!!m
-e
0
r.
-M-1 K 0
-m -m -
-e U1 T
-111 -A
-e
a b
(4.20)
2Af e 1o
8 =H lfe + B (4.21)
Al ~o
Magnetic Bearings 555
FORCE UN
300.
200.
i0~.
21111.. HZ 1.4- K
FORCE f/N
300.
-- :r-t------~----.jl
.s ..... , '' t
IIIRGIIP lotmm
l
Fig. 8: Calculated (----) and Fig. 9: Force/current factor de-
measured (---,o) bea- pending on the excitation
ring force of a radial frequency
bearing (diameter = 80 mm,
a = 40 mm, n = 720)
556 G. Schweitzer
pmax =A h T (4.24)
where Pmax is the cooling power for the temperature-difference
T, A is the surface of the magnetic bearing and h is a cooling
factor, depending on the cooling medium. For a given application
one can of course use much more sophisticated models. In most
cases, the eddy-current 1asses and the hysteres i s-1 asses in the
stator are negligible compared to the copper losses. Therefore we
can calculate the maximum ampere-turns for one magnet by
P
max Aw f a
e max = n ; max = -4- PT (4.25)
w
FORCE f'/N
500.0
400. 0
300.0
POLEWIDTH b
200. 0
100. 0
2."
"· " •5 1. " 1. 5
AIRGAP lo/"'"'
F(x,i) = kX x + k. i
1
(4.26)
where the force/displacement factor k and the force/current
factor k. depend on the premagnetization~ the nominal air-gap and
the geom~try of the actuator. For small displacements and small
control currents linear characteristics as in fig. 11 are
obtained. A disadvantage of this linearization is that the
premagnetizing current causes higher copper losses: the advantage
is that this bearing is especially well applicable to control
purposes.
558 G. Schweitzer
t.O
Bearing Force ~
0
0.8
-o. a 0.8
5. SYSTEM ASPECTS
5.1 Control of rigid and elastic rotors
Any control requires information about the motion of the rotor.
Here the rotor position is measured by distance sensors. The
signals for the displacement speed can be produced in three
different ways: by electronically differentiating the
Magnetic Bearings 559
00
J = f ( 5. 1)
0
where~ and ~fare motion and control variables (eq. 4.7). When
all the state variables x, displacements z and their time
i are used, we have a complete -tor central) state
derivatives -s
feedback. Then the control ~f depends linearly on the state
(5.2)
where K is the gain matrix. In our case (fig. 12a), with 8 state
variables and 4 control variables for the two radial bearings,
the gain matrix has 32 coefficients. Obviously this large number
of coefficients makes any realization of the controller somewhat
cumbersome.
a) b)
( 5. 3)
~0 = -T0 -m
u u -r
z + -T0 -r
-m z , ~v = -Tv -m i
u -m i
u -r
+ -Tv -r
Magnetic Bearings 561
T T
M
-m-m z K z = -mJ:.
+ -m-m U n =-Um w
TA- -
(5.4)
~0m = -T0 -m-
u z , ~v m = -Tv -m-m
u z
IM
M, z, + K, Zr
cs OS
Control
y
Furthermore the control vector w does not only act on the modeled
part but on the real full-order system itself which obviously
contains the unmodeled part as well. This influence is called
"contro 1 spi 11 over". These spi 11 over terms can change and
deteriorate the behaviour of the real system and even destabilize
it.
The objective of the control design now is to derive the control
on the basis of the reduced-order system without exact knowledge
about the residual, unmodelled part.
It is advisable to use a direct output feedback, i.e. the control
forces depend linearly on the displacement and velocity
measurements
(5.5)
with the yet unkown gain matrices G.• Introducing this control
law and the measurement equations -~4. 15) into the full-order
system equations (4. 13) leads to the full-order closed-loop
system
0
!1 ••
~ + !!T IA §v Iv !! • (
! + ~ + !!
T
IA §o Io !! ) ! =- ( 5. 6)
a)
b)
imi'IH
imaH
0.1
0.01
f/Hz 1/ms
: fg
10 100 1000 0 2 4 6 8
The other rotor losses, the air losses, depend on the geometry of
the rotor and its surrounding parts. For an example the run-down
curve of a magnetically suspended rotor has been measured in air
and in vacuum (fig. 16). The braking torques caused by the various
effects are shown in fig. 17.
The stator losses are dominated by the copper losses due to the
resistance in the coils. According to equations (4.24, 4.25) they
are proportional to the square of the current i in the coils,
plc = kc * i2 (5.9)
6. APPLICATIONS
Magnetic bearings have been used for rotating machines in vacuum-
techniques because there wi 11 be no contamination by wear or
lubricants. One of the first industrial applications was for a
turbomolecular pump. An extremly low vibration level ( o.o5 mm/s
at the pump), an operational speed of 5oo Hz, and an ultrahigh
vacuum of 10 E-10 mbar free of any hydrocarbon contamination are
outstanding specifications. Now several companies are developing
and building such pumps.
566 G. Schweitzer
U/min
51!11!10..1!1
41!11!10..1!1
31!11!10..1!1
21!11!11!1.1!1
ll!ll!li!J.I!l
.l?l225
• l?l2
.l?l175
.l?llS
---
• IZil
__,....-
2 ...-....-
.IZIIZI75
....-
~....-
.IZIIZIS
-·-------------
3
.IZIIZI25
121. 121
Drehfrequenz/Hz
l?l. l?l 2IZI. 1Z1 4IZI. 1Z1 60. 0 80. 0 10IZI. l?l
Uocuum Housing
Rotor
Motor
Rodiol Beorlng
A•dol Beorlng
Rodiol Sensor
Frome
AKiol Sensor
Ouen
Furnoce
a
u.t
c
:::>
1-
...J
ll.
~
<
b
FREQUENCY
Fig. 21: Frequency response Fig. 22: Trajectory of rotor
a - undamped, shaft center for a
b - with active damper rotor touching down
on the bearing
Magnetic Bearings 569
7. CONCLUSIONS
References
/BY 46/ Beams, J.W., Young, J.l., and J.W. Moore: The production
of high centrifugal fields. J. Appl. Phys., 1946,
pp. 886-890.
H. Inetier
University of Kassel, Kassel, FRG
ABSTRACT
1• INTRODUCTION
An important field in practical mechanical engineering is the
construction of dynamically loaded machines e.g. compressors,
combustion engines, gas and steam turbines and other machines
in which rotating or accelerated components are used /1/, /2/,
/3/. These components are loaded by the forces of the machine
process and/or their own inertia forces. Thus, high stress levels
may occur in particular when the frequency spectrum of the load
excites the component or the machine in resonance.
In a lot of practical cases, these exciting frequency spectra are
time-independent. This is valid for instance for those machines
which run with constant angular velocity and with no change in
the machine process. However, important cases exist where the
exciting frequencies change with time /4/. This non-stationary
case occurs for instance in turbines which run up or run down, in
compressors with external load increase or decrease or in electro-
572 H. lrretier
0° 360°
L_ {7 {7 {7 __j
=.I~---------~~--...!_ stator blades
~- rotor blade
P rotation
Js
Ks SPEED
.__,:-+--=-+~~1--......;--t INCREASER
JlO J,
MOTOR
l~RPM
20000 HP
COMPRESSOR
5500 RPM
2
al
....1
~
w· 0
:::>
0
a:
0
1-
·2
FIRST TORSIONAL
FREQUENCY EXCITED BY
TWICE SLIP FREQUENCY
COMPONENT
·4
0 2 4 6 8 10 12
TIME tSECI
zc
1 ~=o.o31
10
la=0,02 w~l
,,
5
-5
'' ~/
!: = 0,02
15
11
_ j 0 (a>O)
0 l100 (a<O)
10
72 a.=O
IFF(-dl ~ = 0,02
.2
muU<P
c:o
10 (a>O)
=
(a<O)
2 <P' = n0 -+a:r
Fig. 9: Force at the foundation of a spring-mass-system excited
by unbalance during a passage through resonance by linear
increasing and decreasing exciting frequency /5/
4.2 RUN-UP AND SHUT-DOWN OF A FLEXIBLE ROTOR
To continue our consideration of an unbalanced flexible rotor
in chapter 2.3, now its response is discussed for the case of a
passage through resonance due to a constant angular acceleration.
The describing, linearized equations of motion (2.6) are of
similar type as for a force excited spring-mass-system. Thus, the
solution is possible in a way comparable to that one described in
the previous chapter. Fig. 10 shows the amplitude of a Laval-
rotor running through its critical speed during a run-up (a. > 0)
and a shut-down (a.< 0). The damping ratio is~ = 0,02 and curves
are plotted for various values of the constant angular accelera-
tion a. Again, the typical behaviour of a linear system running
through resonance is visible. The maximum of response occurs
after the critical speed for stationary conditions is reached,
and this maximum is less than in the case of stationary rotation
of the rotor. For acceleration lal > 0,1 no typical resonance re-
582 H. Irretier
1: ; 0,02
'""
t:O
blade force
Fig. 12: Rotating blade and load due to two arcs of partial
admission /11/
: -resonance points
I
600.0
f 1[Hz]
400.0
300.0
200.0 2
100.0
.0
.o 10 25.0 50.0 60 75.0 fofHzl 100.0
30
dma. (I) 7 - i ---- 3
2 -1- 2
20
10
00
- 10
-2 0
-3 0 U1=const) - - 9.41
00 25 5.0 7.5 t[s] 10.0
6. REFERENCES
/1/ !MechE: Proc. 3rd. Int. Conf. on Vibrations in Rotating
Machinery, York (England), Sept. 1984
i21 IFToMM: Proc. Int. Conf. on Rotordynamic Problems in Power
Plants, Rom (Italien), Sept./Oct. 1982
/3/ Holzweissig, R.; Dresig, H.: Lehrbuch der 11aschinendynamik.
Springer-Verlag Wien/~ew York, 1979
/4/ Goloskokow, E. G.; Filippow, A. P.: Instationare Schwin-
gungen mechanischer Systeme. Akademie-Verlag, Berlin, 1971
(translated from Russian)
/5/ i~arkert,R.; PfUtzner, H.: An- und Auslaufvorgange einfacher
Schwinger. Forsch. Ing.-Wes. 47,4 (1981),117-125
/5/ Gasch, R.; PfUtzner, H.: Rotordynamik. Springer-Verlag
Berlin/Keidelberg/New York, 1975
171 ~arkert, R.: Biegeschwingungsverhalten unwuchtiger elasti-
scher Rotoren bei der Resonanzdurchfahrt. VDI-Bericht 381
(1980) '155-160
/8/ Traupel, W.: Thermische Turbomaschinen. Springer-Verlag
Berl in/Heidelberg/rlew York, 1968
/9/ Rao, J.S.: Turbine Blade Vibration. Shock and Vibration
Digest 12 (1980), 19-26