Journal-Bearing Databook (PDFDrive)
Journal-Bearing Databook (PDFDrive)
Journal-Bearing
Databook
With Contributions by
T. Someya, J. Mitsui, J. Esaki, S. Saito,
Y Kanemitsu, T. Iwatsubo, M. Tanaka,
S. Hisa, T. Fujikawa, H. Kanki
This work is subject to copyright. All rights are reserved. whether the whole or part of the material is
concerned, specifically the rights of translation, reprinting. reuse of illustrations, recitation. broadcast-
ing, reproduction on microfilms or in other ways. and storage in data banks. Duplication o f this
publication or parts thereof is only permitted under the provisions o f the German Copyright Law of
September 9. 1965, in its version of June 24, 1985, and a copyright fee must always be paid. Violations
fall under the prosecution act of the German Copyright Law.
© Springer-Verlag Berlin Heidelberg 1989
Originally published by Springer-Verlag Berlin Heidelberg New York in 1989
The use of registered names. trademarks. etc. in this publication docs not imply. even in the absence
of a specific statement. that such names are exempt from the relevant protective laws and regula tions
and therefore free for general use.
Typesetting : Daten- und Lichtsatz-Servicc, Wurzburg
v
Companies Supplying Data (in alphabetic order)
Ebara Corporation
Hitachi, Ltd.
Hitachi Zosen Corporation
Ishikawajima-Harima Heavy Industries Co., Ltd.
Kawasaki Heavy Industries, Ltd.
Kobe Steel, Ltd.
Mitsubishi Electric Corporation
Mitsubishi Heavy Industries, Ltd.
Mitsui Engineering & Shipbuilding Co., Ltd.
Nippon Kokan K. K.
Nippon Seiko K. K.
Toshiba Corporation
Toyoda Machine Works, Ltd.
Ebara Corporation
Hitachi, Ltd.
Hitachi Zosen Corporation
lshikawajima-Harima Heavy Industries Co., Ltd.
Mitsubishi Electric Corporation
Mitsubishi Heavy Industries, Ltd.
Mitsui Engineering & Shipbuilding Co., Ltd.
Toshiba Corporation
VI
Preface
Journal bearings, which are used in all kinds of rotating machinery, do not only
support static loads, such as the weight of rotors and load caused by transmitted
torque of reduction gears, but are, in addition almost the only machine element that
is able to suppress various exciting forces acting on the rotating shaft.
As rotating machines have become large and multi-staged, while compactness,
high speed, and high output have also been realized in recent years, not only has
the bearing load increased, but also the magnitude and variety of exciting forces.
Therefore, the role and importance of journal bearings have increased tremendous-
ly.
In particular, for the design of rotating machines with low vibration levels and
high reliability, knowledge of the exact characteristic data of bearings, and especial-
ly of the stiffness or spring coefficients and the damping coefficients of oil films in
bearings, is essential.
However, the amount of reliable data now applicable to practical design is
limited. Through the activity of the Research Subcommittee on Dynamic Charac-
teristics of Journal Bearings and Their Applications (designated as PSC 28), estab-
lished and organized in June 1979 through May 1982 within the Japan Society of
Mechanical Engineers (JSME), these coefficients, together with static characteris-
tics, have been calculated and also measured on a number of new test rigs. These
data are the principal contents of the present databook, in which almost all types
of bearing actually in use are included. The book can be said to be the fullest
collection of information currently available. In addition it contains explanations
on the characteristics of oil films, and also examples of practical applications, so
that it will be useful to either newcomers to the field or experts, as well as to research
workers and designers.
As mentioned above, the characteristic data in this databook are based on the
hard work of the JSME Research Subcommittee on Dynamic Characteristics of
Journal Bearings and Their Applications. The names and affiliations of the mem-
bers are listed before, and their help is gratefully acknowledged. Thanks are also
due to the companies listed for allowing us to publish the data.
This volume is an English version, with some revisions, of the Japanese original,
published under a licence agreement between JSME and Springer-Verlag. The
translation was supplied by the authors of this databook, and their selfless dedica-
tion is much appreciated. Special thanks are due to Prof. M. Tanaka and Dr. S.
Saito, who have devotedly acted as secretaries throughout the progress of the
project since 1979. Timely help from Mr. M. Takahashi, the manager of JSME, and
Mr. A. von Hagen, acquisitions editor of Springer-Verlag, in publishing the book
is acknowledged with thanks.
VII
Contents
Nomenclature . . . . . . . . XIII
2. Calculated Data 13
IX
Contents
X
Contents
4. Explanations . . . . . . . . . . . . . . . .
. . . . . . 231
4.1 Method of Calculation for Bearing Characteristics (By J. Mitsui) . 231
4.1.1 Generalized Reynolds Equation for Laminar Flow 231
4.1.2 Reynolds Equation for Turbulent Flow . . . 232
4.1.3 Method of Calculation for Reynolds Equation 233
4.1.4 Boundary Condition for Oil Film Pressure 234
4.1.5 Calculation of Static Bearing Characteristics . 236
4.1.6 Method of Calculation for Dimensionless Stiffness and Damping
Coefficients of Oil Film . . . . . . . . . . . . . . . . 237 . .
4.2 Methods of Measurement for Journal Bearing Characteristics
(By J. Esaki) . . . . . . . . . . . . . . . . . . . . . . . . .
240
4.2.1 Outline of Performance Characteristic Tests of Journal Bearing 240
4.2.2 Methods of Measurem'e nt for Oil Film Coefficients . . 242
4.2.3 Example of a Test Rig for Dynamic Characteristics of
Journal Bearings . . . . . . . . . . . . . . . . . 243
4.2.4 Experimental Results Obtained with Circular Bearings with
Two Axial Grooves . . . . . . . . . . . . . . 245
4.3 Estimation of Effective Oil Film Temperature (By S. Saito) 248
4.3.1 Estimation from Calculated Data 248
4.3.2 Estimation from Measured Data . . . . . . . . . 249
XI
Contents
XII
N otnenclature
C = R - r Radial clearance
Cb Assembled clearance
CP Machined clearance
C;j = Cpwc;)W Dimensionless damping coefficient
D Bearing diameter
F Oil film force
H Heat generated in the film per second
K;j = C P kij/ W Dimensionless stiffness or spring coefficient
L Bearing width
N Rotational speed (rps)
Qe = qe/ (RwCPL) Dimensionless circumferential oil flow at the trailing edge of
pad or at the boundary of film rupture
Q1 = q1/ (RwCPL) Dimensionless forced oil flow
Q, = q,/ (RwCPL) Dimensionless side flow of oil
R = D/2 Bearing radius
Re = CPU fv Reynolds number
S = 11NLD (R /Cp) 2 / W Sommerfeld number
T Temperature (0 C)
U Journal surface velocity
W Bearing load
cij Damping cofficient of oil film (i,j = x, y)
d Journal diameter
e Journal eccentricity
f Oil film force
h Frictional coefficient on journal
h = C + e cos () Oil film thickness
hmin Minimum oil film thickness
kij Stiffness or spring coefficient of oil film (i,j = x, y)
mP = 1 - Cb/CP Preload factor
p Oil film pressure
p1 Feed pressure of oil
Pm = W / (LD) Mean bearing pressure
r = d/ 2 Journal radius
t Time
x, y, z Vertical, horizontal and-axial coordinates
oc Angle of pivot position from leading edge
jJ Pad angle
e = e/Cb Journal eccentricity ratio
() Circumferential angle
Jl Viscosity of lubricant
v Kinematic viscosity of lubricant
Q Density of lubricant
<jJ Attitude angle
If/= CP/R Clearance ratio
w = 2 nN Angular velocity of journal
Bold letters denote vectors or matrices
XIII
1. Introduction
This databook is concerned with the static and dy- Chapter 1 deals with general and fundamental
namic characteristics of oil films in the journal bear- items relating to the databook as a whole and how
ings most frequently used in high-speed rotating to use it; Chaps. 2 and 3 contain the characteristic
machinery. These characteristics, which exert a fun- data of journal bearings; and Chap. 4 explains the
damental influence upon the vibration and stability procedures applied to obtain the characteristics, the
of rotors, have been obtained both theoretically influences of various parameters on them, and their
and experimentally and are presented in tables and applications to vibration analysis of rotors.
graphs. The book is intended to contribute to the Numerical values are given wherever possible in
proper and vibration-free design of rotating ma- the form of dimensionless numbers, and otherwise
chinery. with SI units in most cases.
g) Others
(I) Hydrodynamic
bearing
With 3 or 4 pads
( i i ) Tilting pad - {
beanng Foil bearing
LOP or LBP
Floating bu sh bearing
Journal iii)Others .. J
~ Pourous metal bearing
0
bearing
(2) Hydrostatic
bea nng
--1r - - - - - - - - - With 4 or 5 pockets
00
0
,___ _ _ _ _ _ _ _ _ With axial grooves
between pockets
If/ =
cR: clearance ratio, (2)
c
m = 1 - C b : preload factor. (3)
p
a Machining b Assembly
Geometrical parameters have to be defined fur-
Fig. 1.2a, b. Main dimensions of journal bearings. L , length ther as needed, e.g., number of segments (pads),
angular breadth of segments, offset of segments,
load direction relative to segment, positition, shape
and size of oil grooves.
Depending upon the main dimensions just de-
designated as C:
scribed, limiting space is obtained in which the cen-
C = CP = Cb = R- r, (1 ') tre of the journal can move freely. This space is
L: bearing length (m). called the clearance "circle", and is in fact a circle
with radius C in the case of a 360° cylindrical bear-
From these main dimensions the following di- ing; it has a different form in other cases, however,
mensionless geometrical parameters, which govern as shown in Fig. 1.3.
2
1.2 Static Characteristics of Journal Bearings
Clearance circle
This section explains how to determine the charac- static characteristics are also required for determi-
teristics of oil film in a journal bearing in static nation of the dynamic characteristics.
equilibrium with the aid of this databook. These
3
1 Introduction
1.2.1 Equilibrium Position where J..l is the absolute viscosity of the lubricant
of Journal Centre (Pa · s). For f..l, the appropriate mean value has to be
taken when, because of heat dissipation, variations
in temperature and hence in viscosity are not negli-
A journal supported by journal bearings and rotat- gible. The units of variables W, N, D, and L are N,
ing with a constant rotational speed N is assumed 1/s, m, and m, respectively. For various bearing
to be loaded vertically and downwards by an exter- geometries these functional relationships
nal force of constant magnitude W, as shown in
Fig. 1.5. e = e(S)
(7)
Then, the centre 0 1 of the loaded journal shifts 4J=4J(S)
from the centre OB of the bearing, so that a conver-.
gent clearance space is formed between the journal are presented m tables and graphs in Chaps. 2
and the bearing surface, into which the oil is forced and 3.
by its viscosity and adhesion, so that hydrodynamic After determination of the values e and ¢, the
pressure,p, is generated in the oil film. This pressure minimum film thickness, hmin ' and its circumferen-
is called the wedge film pressure. tial location are obtained for the given geometry of
In the remaining divergent clearance space the the aforementioned clearance circle. For example,
lubricant flow is usually cavitated by aeration (ex- in the case of 360°cylindrical bearings the following
pulsion of air or gas), and the pressure in this space equation is obtained:
is close to the atmospheric pressure under normal hmin = C(1 - e) (m) (8)
operating conditions.
The resultant force P of pressure p is balanced and the location of hmin is given by ¢.
with the external load W. Because of the asymmetri-
cal distribution of p against the line OB 0 Jo the cen-
tre 0 1 of the journal does not move vertically or in
the direction of the applied load, but shifts, as 1.2.2 Frictional Coefficient
shown in Fig. 1.5, along a sort of semicircle when of Journal Bearings
the load or rotational speed is changed. For a given
set of geometrical parameters, including L : D and
mP as defined in Sect. 1.1.2, the eccentricity ratio e !he tangential stress of the oil film acting on the
and attitude angle 4J are functions of a dimension- JOUrn~ results in frictional torque T., (N m). If we set
less parameter called the Sommerfeld number, de- T., = F1 r, the coefficient of friction of the journal
fined as follows: bearing, f1o is defined as follows:
(6)
fJ=~= TJ. (9)
W rW
4
1.2 Static Characteristics of Journal Bearings
and therefore
(13)
0.15
tiE-
0.10
/"
/
J~':1::-
0.08
<>:- /
0.06 /
'=>'oo
II
eli' 0.04 /
/
/
0.02 0. 0 1 '---'---'----'----'-----''---'---'----'----'---'
0 0.2 0.4 0.6 0.8 d/L 1.0 0 0.2 0.4 0.6 0.8 a/L 1.0
a Circular hole b Axial oil groove
5
1 Introduction
X
and axial groove.
In case of 360° cylindrical bearing with central/
circumferential grooves the following equation ap-
plies:
3
HHJ
nC · ( 1 + 3 e2) · p 1 (m 3 /s)
=~
4
qf 2 (19) Fig. 1.8. Heat balance in journal bearing
3f.i,- H: heat dissipated in oil film
D H 1 : net heat flowing out with oil
H 2 : heat transfer from bearing housing to surroundings
where I.: is the effective length of bearing or L minus H 3 : heat conduction from bearing housing to foundation
the axial width of the groove. H 4 : heat transfer to coolant
Therefore, the total oil flow is obtained from H 5 : heat transfer from journal
Eqs. (11)-(19). H 6 : heat transfer from machine parts to bearing
8: mean temperature of oil film
q = q(J) + qf = q. + (1 - x) qe + qf 8A: surface temperature of bearing housing
80 : temperature of surroundings
or W: external load
q = RwCPL{Q. + (1- x) Qe} + q1 (m 3 /s) (20) 1}: frictional torque
6
1.2 Static Characteristics of Journal Bearings
Since the heat H dissipated in the oil film per Eq. (24) becomes
second is given by the product of frictional torque
H 2 = &A(.9- 9 0 ) (W). (28)
'F.J and angular speed w, Eq. (9) gives:
Therefore, Eqs. (21)', (22), (23), and (28) yield, for the
H = T1 w =f1 rWw (W). (22)
thermal equilibrium of a bearing, the following
It should be pointed out that there are also cases, equation:
such as thrust bearings of the pad type, in which
heat dissipation through churning is high [6] and f 1 r Ww = C(Jq(92 - 91 ) + &A(9- 90 ) (W). (29)
has to be added to the right-hand side of Eq. (22). No reliable relation is yet established between the
The net heat H 1 flowing out with the oil is given mean temperature 9 and the oil inlet and oulet
by temperatures, 91 and 92 • But, as long as the temper-
H1 = Ceq(.9 2 - .9 1 ) (W) (23) ature difference .9 2 - 9 1 is not too wide, a weighted
mean can be used for 9, so that:
where
9 + (1 -A) .92 ( C) (30)
c: specific heat of oil (J/kg oq, = ).91 0
7
1 Introduction
( Start )
f
CD Bearing type
D,L, r, C,, Main parameters to be
t - - given at the beginning
W, N, pf, o'ioo . . I> of design
c,Q, A, a, x, ..!
i
CD LID, R=~
C"=R-r ... Eq. (1)
Calculation of
S_ · · · Eq. (2) t - - geometrical parameters
1/J= R
m =I-S
Cp · · · Eq. (3)
(J
i
G) Assume{} {}: Mean temperature of oil film
1
<D p=fi ({}) .. For example,
Eq. (35)'
fl. is calculated according
to temperature-viscosity relationship
~
CD S= 1!:!!2!:__
2
ljJW
· · · Eq. (6) Sommerfeld number S, is calculated
i
c= c(S) }
· · · Eq. (7)
¢=¢(S) r - - - According to data in Chapter 2 or 3
for example,
hmin
· · · Eq. (8)
i
CD !J= 1/Jfunc(S) · · · Eq. (10)
1 - - According to data in Chapter 2
Qs= Qs(S)
} Eq. (16)
Q, =Q/S)
~
® q= Eq. (20) Oil flow q is calculated
® {}
oow
•
= Right hand side
of Eq. (32) {}"'w: New value for {}
@)
1{}-.,?>ncw l <0
• (For
example
Yes @ {}={}new @ •
Dynamic characteristics are
calculated by interpolation
1 oq {}2 = {} -,\,?>I of data given in Chapter 2
· · · Eq. (33)
1-..l or 3 using the values of S
tNo determined at @
@ {}=<&new {} <&--<&-
=T=T ··· Eq. (34)
l
2 -{}1
'
@ ~r {}obtained at@, calculations
4 - @are repeated once again
i
( Stop
@ Change values
of appropriate
parameters
Yes @
excessively high? 'iNa
Is{}, {}2 --&-1 or frictional loss
given at G) @
Yes
Is hmin too small"!
Yes @ .No
Is thermal expansion excessively large?
'No
size of bearing and operational parameters. Start- After passing through Q3) "' @ one can proceed
ing from step CD in Fig. 1.9, the condition of heat to @, where dynamic characteristics are obtained
balance is satisfied at (j)'. from the data given in Chaps. 2 or 3 using the value
If the mean temperature 9, temperature increase for S determined above at @ '.
9 2 - 9 1 or frictional loss is excessively high at Q3), As will be explained later, in Sect. 4.3.1, it is also
cooling oil flow should be increased and/or clear- possible to estimate the temperature increase by
ance and other sizes of bearing corrected, and the assuming the ratio !;_ of the heat H 1 to the heat
calculation should be repeated from step aJ. If at (jJ) dissipation H: fi = H 1 / H, with approximately
the minimum oil film thickness hmin is too small to fi = 0.6 in most cases. Then Eq. (32) has to be re-
maintain full film lubrication, the parameters in CD placed by the following equation
should be changed. Further, if at @ the change of
bearing clearance by thermal expansion is not neg-
ligible, the change of parameters C, If/, m given at aJ 9 = (1- A)f1 rWwl;_ + 91 . (32)'
should be also taken into account. c (!q
9
1 Introduction
10
1.4 How to Use the Characteristic Data of Bearings
11
1 Introduction
journal rotation changes, but they have to be calcu- Sects. 1.2.4 and 4.3. This databook is intended to
lated afresh if the surfaces are asymmetrical. contribute to the prediction of vibrational be-
It goes without saying that the coordinates x, y of haviour, including that of the stability of a rotating
the journal centre and the positive directions of shaft supported by a journal bearing. How much
components Fx, FY of the oil film force acting on the the dynamic characteristics of oil film influence the
journal (not on the bearing) are those defined in vibrational behaviour depends upon the shape of
Sect. 1.1.3. In cases of other coordinates systems the the shaft system. When this influence seems to be
conversion given in Sect. 4.1 is needed. large, as in the case of an overhanging shaft, the
accuracy of calculated vibrational behaviour
should be checked carefully.
1.4.4 Some Notes on Application In general, calculated and experimental data
have to be compared against each other, checked,
of the Data and then applied. Accurate prediction of other
forces than that of the bearing, such as the forces of
For effective use of the present databook, it is most labyrinths, discs, seals etc., is also important, as is
important to give the appropriate mean tempera- the correct evaluation of the influences of structures
ture of the oil film, as already mentioned. Methods that support the bearings.
of predicting the temperature are described in
References
1. Glienicke J (1970) VDI-Z, Fortschritt- Berichte 1: 22 5. Vogelpohl G (1967) Betriebssichere Gleitlager. Springer,
2. Birembaut MY (1977) Raideur et Amortissement des Paliers Berlin Heidelberg New York
Hydrodynamiques. Mem Techn CETIM (Paris) 32 6. Neal PB (1982) Heat transfer in pad thrust bearings. PIME
3. JSME (1982) Dynamic characteristics of journal bearing and 196: 217 - 228
vibration of rotating machineries - Application to design 7. Schaffrath G (1969) Forschung im lngenieur-Wesen 35 (6)
based on the newest data (in Japanese). JSME, Tokyo (Text- 8. Someya T (1976) Transaction of JSME 42 : 360, 2599 (in
book no. 548 for the JSME lecture course) Japanese)
4. Hirano F (1957) Transactions of JSME 23: 134, 717; Junkat- 9. Someya T (1976) Transactions of JSME 42 : 362, 3303;
su (JSLE) (1962) 7: 5, 225 (in Japanese) Junkatsu (JSLE) (1975) 20 : 3, 191 (in Japanese)
12
2. Calculated Data
In this chapter, calculated data for 59 kinds of bearing geometries is collected. The format of the data
collections, except those for calculations 56-59, is as follows.
First page : Geometry of the bearing whose data is calculated, outline of the calculation method, and
tabulation of the calculated results
Second page: Figures for eccentricity, attitude angle, and oil flow versus Sommerfeld number, and for
friction coefficient versus Sommerfeld number
Third page: Figures for spring coefficients versus Sommerfeld number and for damping coefficients
versus Sommerfeld number
The main symbols used are as follows (see Nomenclature).
S: Sommerfeld number=
J.1NLD
W (R)
Cp
2
e: Eccentricity ratio = ~b
</>: Attitude angle (degrees)
Q: Nondimensional flow= (Rm~PL)
f 1 If/- 1 : Friction coefficient for journal
Kii: Nondimensional spring coefficient = C~ij (i,j = x, y)
13
Calculation No. 1
5.96 . 0750 82.6 .0663 l . 00 111.5 1.77 l3 . 6 -13. 1 2 .72 27. 2 2.06 2.06 14.9
4.43 .100 80.4 • 0880 .999 83.0 1. 7 5 10.3 -9.66 2.70 20.6 2 ,08 2,08 11.4
2.07 .200 71.6 .170 .999 39.5 1. 8 8 5.63 -4.41 2.56 11.2 2.22 2.22 6.4 3
l. 24 .300 64.1 .243 • 997 24.4 2.07 4.27 -2.56 2.34 8.50 2.32 2.32 4.73
• 798 .400 57.5 .309 • 986 16.4 2.39 3 .75 -1. 57 2.17 7.32 2 .24 2.24 3,50
.517 • 5 00 51.1 .366 . 967 11. 3 2.89 3.57 -. 924 2.03 6.81 2.10 2.10 2.60
.3?3 .600 44.7 .416 • 939 7.75 3. 65 3. 62 -.427 1.92 6.81 2.08 2.08 2.06
.187 .700 38.2 .459 .900 5.13 4. 92 3.88 .0235 l. 83 7.32 2.16 2.16 l. 70
.135 .750 34.9 .478 .875 4.07 5. 90 4.11 • 2 58 1.80 7,65 2 .10 2.10 l. 46
.0315 .900 22.7 .524 .758 l. 56 l3 .8 5.76 l. 24 l. 72 10.6 2.03 2.03 ,846
.00499 .975 12.5 .543 .663 .54 43. 3 9.61 4. 2 6 l. 99 22.4 2.93 2. 93 .637
14
Calculation No. 1
0.8
0.6
~ cu Qs -e.
0.4 0.4
30°
1s•
0 0
0.005 0.02 0.05 0.1 0. 2 0.5 2 5
s
lOOr----------------------------------------------------------------,
15
Calculation No. 1
s
lOOr---- ----------- ----------- ----------- ----------- ----------- ---,
50
20
·~
0.02 0.05
s
16
Calculation No. 2
Static Characteristi cs
s Spring Coefficients Damping Coefficients
c ¢ Qs Q, !J<P-1 Kxx Kxy Kyx Kyy Cxx Cxy Cyx Cyy
2.95 .0500 83.0 . 0351 .999 55.1 20.1 -19.8 2.90 ~0.7 1.65 1.65 1~.5
2.38
1. 96 . 075 80.8 . 0523 .999 36.6 2.18 13.5 -13 .o 2,88 27.4 1.66 1.66 9.81
1. 45 .100 78.7 .0692 • 998 27.3 2.04 10.2 -9.62 2. 77 20.8 1. 67 1. 67 7.55
.689 .200 69 . 9 .132 • 998 13.2 1. 99 5.45 -4.38 2.56 11.3 1. 78 1. 78 4.39
.424 .300 61.0 .181 .998 8.41 2.20 3.96 -2.52 2.31 8.67 1. 99 1. 99 3.50
.284 .400 54.6 .220 . 991 5.93 2.37 3.34 -1.60 2.06 7.57 2.10 2.09 2 . 94
.194 .500 49.1 .251 • 972 4.37 2.67 3.09 -1. 00 1. 89 6.99 2.00 2.00 2.31
. 0469 .BOO 32.1 .302 .855 1. 69 6.10 3.82 .420 1. 74 8.08 2.15 2.15 1. 34
.0195 .900 24.0 .300 .770 1. 04 11.9 5.28 1. 32 1. 80 10.1 2.06 2.06 .916
.00866 . 950 18 .0 .297 . 710 .684 24.0 7.81 2.79 2.01 14.2 2 .44 2.44 • 791
17
Calculation No. 2
goo
0.8
75°
0.6 60°
45°
0.4
30°
0.2 0·2
15°
18
Calculation No. 2
....
~
20
....
0
19
Calculation No. 3
16.8 .0500 87.5 .0232 2.00 294 1.90 20 . 0 -19.0 1.53 40.1 1.43 1.50 36.4
8.25 .100 85.3 .0467 1.99 145 1.88 10 . 1 - 9.53 1.52 20.2 1.52 1.51 18.8
3.90 .200 80.7 .0942 1.98 69 . 9 2 .03 5.34 - 4 . 30 1.46 10 . 3 1.46 1.46 8.87
2.35 .300 75.6 .143 1.96 43.3 2.32 3.84 - 2.16 1.31 7.21 1.27 1.26 4.90
1.49 .400 69.6 .196 1.93 28.7 2. 63 3, 21 - 1.25 1.22 5.76 1.19 1.20 3.21
1.16 .450 65.9 .224 1.90 23.4 2.92 3. 02 - • 737 1.10 5 . 33 1.02 1.02 2,27
.904 .500 60.6 .255 1.67 16.5 3.37 2.89 - .200 .939 4.99 ,742 .744 1.33
.658 .550 53.4 .287 1.62 14.0 3.65 2.96 .0660 .940 5.05 . 704 .709 .954
.469 .600 45.6 .315 1. 76 10.5 4.39 3.21 .166 1.04 5 .49 . 844 .850 ,867
.336 .650 39.3 .337 1.71 7.94 5.02 3.51 .267 1.14 6 . 10 1 .01 1 . 02 .832
.239 .700 34.2 .357 1.66 6.05 5.78 3. 86 .319 1.24 6.89 1.22 1.23 .835
.111 .800 26.7 .392 1.57 3.42 6.04 4.66 . 361 1.46 9.28 1.86 1.87 .940
.0353 .900 19.7 .424 1.47 1.61 14.9 5.88 .840 1.39 ll.5 1. 74 1. 76 .627
.0129 .950 15.0 .438 1.40 .907 26.6 7.80 1.44 1.50 15.6 1.99 1.97 . 527
20
Calculation No. 3
2.0
1.8
1.6
1.4
1.2
~ 1.0 w -e..
0.8
0-6
0.4
0·2
21
Calculation No. 3
50.--------------------------------------------------------------,
·~
::.::
2
Kyy
0.5
50.------------------------------------------------------------,
C xy
22
Calculation No.4
8.57 .0500 87.9 .0130 2.00 151 2.24 20.0 -17.7 1.26 40.0 1.29 1.25 35.2
4.24 .100 85.9 .0260 2.00 74 .a 2.00 10.1 -8.89 1.26 20.1 1.25 1.24 17.7
2.01 .200 81.8 .0523 1.99 36.2 2.10 5.26 -3.88 1.20 10.3 1.18 1.17 8.04
1.22 .300 77.1 .0797 1.98 22.7 2.26 3.73 -2.28 1.16 7.10 1.14 1.15 5.01
• 794 .400 71.8 .109 1.95 15.4 2.60 3.03 -1.11 1.02 5.56 .916 .920 2.80
.636 .450 68.4 .125 1.93 12.7 2.84 2.82 -.622 .919 5.09 .781 .781 1.97
.498 .500 63.8 .142 1.90 10.3 3.13 2.69 -.290 .856 4.80 .692 .697 1.42
.366 .550 56.3 .161 1.86 7.91 3.58 2. 71 .0346 .824 4. 77 .608 .613 .921
.260 .600 47.6 .178 1.80 5.92 4.11 2.92 .196 .903 5.15 .695 .701 • 767
.186 .650 40.1 .190 1. 74 4.51 4. 71 3.23 .291 1.01 5.79 .853 .860 .723
.134 .700 34.2 .201 1.68 3.50 5.41 3.58 .349 1.11 6.62 1.05 1.06 .719
.0661 .BOO 26.1 .218 1.59 2.12 7.38 4.52 .368 1.34 9.10 1.67 1.68 .824
.0242 .900 19 .a .235 1.49 1.16 13.6 5.37 .855 1.30 10.9 1.55 1.56 .563
.00980 .950 15.2 .243 1.42 .728 24.2 7.37 1.41 1.47 15.3 1.97 1.94 .528
23
Calculation No.4
24
Calculation No. 4
SOr-------------------------------------------------------1
...,
~
...,
0
25
Calculation No. 5
1.72 0.500 6 2.5 0.5 0 4 0. 9 73 27.7 4.5 9 8.55 -2.7 7 5.3 0 1 4.0 6. 1 5 5.9 0 8.44
1.0 9 0.5 50 57.8 0. 5 3 7 0.9 7 0 1 8.3 4.6 7 6.8 0 - 0.8 3 2 4.1 9 1 0.8 4.3 8 4.3 7 5.18
0.7 4 2 0.6 0 0 5 2.7 0.5 6 4 0.9 6 9 1 2.9 4.7 2 5.7 9 -0.2 7 3 3.3 4 9.2 4 3.4 6 3.4 6 3.60
0.517 0.6 50 4 7. 9 0.5 9 4 0.9 56 9. 52 4.8 5 5.2 2 -0.0993 2.8 5 8.08 2.9 5 2.9 5 2.7 2
0.36 2 0.7 0 0 4 3.2 0.611 0.9 52 7. 0 9 5.4 3 5.1 1 0.231 2.5 0 8.30 2.8 4 2.84 2.11
0.249 0.7 50 38.6 0.6 3 1 0. 9 3 9 5.2 7 6.2 0 4.8 5 0.414 2.21 8.0 1 2.6 3 2.64 1.70
0.165 0.8 0 0 3 4.0 0.6 4 6 0.9 21 3.8 5 7.5 7 5.1 5 0.697 2.1 7 8.58 2.4 2 2.4 3 1.52
0.10 2 0.8 50 2 9.5 0. 6 4 6 0.9 2 0 2.7 1 9. 6 9 5.33 0.985 2.0 5 8.95 2.3 9 2.4 0 1.21
0.05 58 0. 9 0 0 2 4.0 0.6 4 4 0.9 0 1 1. 7 9 14. 1 6.28 1.38 1. 9 9 10.8 2.3 3 2.3 4 1.0 I
0.0219 0.9 50 17.1 0.6 31 0.8 6 7 0.98 5 30.0 7.86 2.91 1.82 1 3.0 1.77 1.78 0.5 55
26
Calculation No. 5
50r-------------------------------------- --,
27
Calculation No. 5
50r----------------------------------------,
....,
:<
20~--------------------------------------------~
....,
0
28
Calculation No.6
2.55 .2 63.2 .737 2.59 18.4 2.22 5.94 -1.16 1.20 11.8 1.22 1.23 2.44
1.49 .3 54.6 .997 2.27 11.8 2.42 4.62 - .843 1.27 9.12 1.33 1.34 1. 89
.931 .4 48.1 1. 20 1.97 8.35 2.70 4.12 - .695 1.38 8. 01 1.45 1.46 1.64
.588 .5 43.2 1.35 1. 66 6.11 3.14 3.93 - .564 1.46 7.52 1.57 1.59 1.49
.359 .6 38.5 1.46 1.36 4.94 3.98 2.82 - .285 1.38 7.43 1.65 1. 68 1.34
. 203 .7 33.5 1.50 1.05 3.24 5.24 3.95 - .0348 1.36 7.61 1.61 1.66 1.10
.145 .75 30.9 1.50 .892 2.70 6.23 4.05 .133 1.34 7.87 1.60 1.64 .982
.0984 .8 28.0 1.47 .774 2.19 7.66 4.29 .324 1.34 8.34 1.60 1.65 .879
.0612 .85 24.8 1.41 .574 1.71 9.89 4.73 .560 1.37 9.13 1.62 1.69 .778
.0328 .9 21.0 1.32 .409 1.25 14.0 5.64 .897 1.47 10.6 1. 68 1. 76 .675
.0124 .95 15.9 1.16 .241 . 767 26.6 7.79 1.86 1.60 13.9 1. 79 1.92 .546
29
Calculation No.6
1.0,---------------------------------------- -----------------,
50,-----------------------------------------------------~
30
Calculation No.6
50.---------------------------------------------------------.
50,----------------------------------------------------.
31
Calculation No. 7
4 .72 .1 77.6 .500 2.86 39.7 1.65 10.5 -3.97 1.71 20.7 1.00 1.58 7.26
2.20 .2 67.6 .941 2.56 19 . 6 1. 78 5.85 -2.14 1. 75 11.5 1.59 1.71 4.16
1.31 .3 60.1 1.31 2.26 12.7 1.94 4.54 -1.57 1.82 8.70 1. 75 1.81 3.14
.838 .4 53.9 1.62 1.97 9.04 2.31 3.94 -1.04 1. 73 7.54 1.85 1. 89 2.59
• 539 •5 48.1 1.87 1.67 6.68 2.85 3. 71 - .661 1.66 7.04 1.85 1.91 2.21
.335 .6 42.4 2.05 1.34 4.95 3. 68 3.70 - . 324 1.59 6.90 1. 78 1.82 1.66
.192 •7 36.4 2.16 1.06 3.56 4.97 3.83 .0236 1.51 7.16 1. 74 1. 78 1.31
.139 .75 33.3 2.18 .899 2.96 5.93 4.00 .197 1.51 7.48 1. 74 1.77 1.17
.0947 .8 30.1 2.17 • 740 2.39 7.23 4.37 .3 57 1.55 7.98 1. 73 1. 78 1.03
.0594 .85 26.4 2.13 .578 1.85 9.42 4.86 .615 1.57 8.81 1.75 1.80 .895
.0322 .9 22.2 2.04 .413 1.33 13.7 5.65 1.06 1.57 10.3 1.81 1.88 .768
.0123 .95 16.5 1.86 .242 .799 25.9 7.98 1.97 1. 73 13.8 1.90 2.01 .596
32
Calculation No.7
1.0.---------------------------------------- -----------------.
50.--------------------------------------- ---------------,
33
Calculation No. 7
50~----------------------------------------------------~
5
s
50~------------------------------------------------------~
·~
5
s
34
Calculation No. 8
2.19 .1 71.9 . 271 2.94 14.9 2 . 04 10.7 -1 . 40 .920 21.4 .924 .923 2.82
1.03 •2 59.0 .477 2.65 7.57 2.19 6.07 - .776 .966 12.3 1.08 1.01 1.69
.615 .3 50.4 .623 2.35 5.03 2.38 4.69 - .574 1.04 9.63 1.17 1.13 1.36
.394 .4 44.0 .723 2.05 3.67 2.61 4.17 - .488 1.12 8.51 1.26 1.26 1.21
.256 .5 39.4 . 790 1. 74 2.80 2.92 3. 89 - . 438 1.21 7.98 1.37 1.39 1.14
.165 .6 35.7 .832 1.44 2.18 3.47 3.75 - .342 1.25 7.82 1.51 1. 53 1.10
.100 .7 31.8 .837 1.13 1. 70 4.52 3.68 - .127 1.22 7.94 1.57 1.61 1.00
. 0750 .75 29.7 .824 .964 1.49 5.34 ), 76 .0096 1.22 8.10 1.55 1.60 . 914
.0539 .8 27.4 .798 .802 1. 29 6.49 3.93 ,170 1.24 8.43 1.55 1.61 .833
.0361 .85 24.7 . 755 .636 1.08 8 . 33 4 . 32 .372 1.30 9.08 1.58 1.64 .754
.0213 .9 21.3 .691 .468 .867 12.2 4.97 .829 1.34 10.3 1.65 1. 72 .672
.00922 .95 16.5 .590 • 296 .607 23.8 6.91 1.92 1.49 13 . 5 1.81 1.92 .568
35
Calculation No. 8
QJ w
50,------------------------------------- ---------------,
5
s
36
Calculation No. 8
50~------------------------------------------------------~
50,---------------------------------------------------------,
Cxx "-._
·~
5
s
37
Calculation No. 9
1.65 .1 75.2 .373 2.91 14.0 1.67 10.5 - 2 .59 1.36 21.0 1.06 1.31 4.97
.786 .2 63.8 .677 2.63 7.11 1. 72 5.87 -1.43 1.40 11 . 8 1.42 1.44 2.93
.479 .3 55.8 .907 2.34 4.78 1.89 4.52 -1.08 1.46 9 . 11 1.54 1. 56 2.92
.315 .4 50.0 1.09 2.05 3.54 2.06 3.93 - .904 1.51 7.90 1.64 1. 66 1.98
.212 •5 45.4 1.22 1. 76 2.76 2.49 3.55 - .615 1.45 7.32 1. 72 1. 75 1.77
.141 .6 40.7 1.31 1.45 2.20 3.10 3.40 - .374 1.41 7 .09 1.71 1. 75 1.51
.0883 •7 35.8 1.33 1.14 1. 75 4.15 3.42 - .0613 1.36 7.17 1.65 1. 70 1.23
.0674 .75 33.2 1.32 .982 1.54 4.93 3.54 .0980 1.36 7.41 1.66 1. 70 1.11
.0493 .8 30.3 1.29 .819 1.34 6.00 3.78 .2 62 1.39 7.81 1.67 1.71 1.00
.0336 .85 27.0 1.24 .654 1.13 7.91 4.17 .565 1.41 8 .53 1.71 1. 76 .897
.0202 .9 23.0 1.15 .487 .902 11.4 5.02 .952 1.51 9.87 1.77 1.84 . 781
.00898 .95 17.5 1.02 .315 .627 23.1 7.03 2.06 1.63 1).1 1.91 2.03 .640
38
Calculation No. 9
1.0r----------------------------------------------------,
50r-----------------------------------------------------~
39
Calculation No. 9
50r---------------------------------------------------------,
5
s
50r------------------------------------------------------.
...,
0
40
Calculation No. 10
4.79 .0556 87.2 . 292 .530 104 61.5 33 . 1 -43.3 2. 06 123 -43.3 -42.5 36.6
3.06 .0868 87.2 .293 .524 65.9 39.4 21.7 -28.0 1. 43 80.0 -27.7 -27.7 23.4
1. 93 .137 86 . 9 .196 . 520 41 .6 25.1 13.7 -17.7 1.16 50.2 -17.9 -17.5 15.0
.993 .262 85.9 .312 .502 21.4 13.5 7.50 -9.20 .897 26.5 -8.47 -8.34 8.28
.496 . 438 84 . 4 . 343 .469 10.8 7.88 4.52 -4.60 .831 14.0 -3.42 -3.36 4. 71
.306 .608 80.8 . 370 .441 6.99 6.06 3. 7 3 -2.82 . 940 10.5 -1.13 -1.11 3.45
.209 . 718 77.2 .402 . 404 4. 91 5.66 3.60 -1.67 1. 16 8.74 -.0503 - . 0240 2.72
.0507 .970 56.8 .436 • 353 1. 61 9.47 5. 46 . 945 1. 94 9.03 2.29 2.31 1. 38
.0303 .995 48.6 . 424 . 359 1. 12 13.2 6.45 1. 72 2.12 10 . 4 2.58 2.62 1.16
.0207 1.01 42.5 .407 . 370 . 891 17.2 7. 81 2.33 2.24 11.8 2. 70 2.73 1.02
.0100 1.01 34.1 . 372 .400 . 590 29 . 6 9.68 3. 28 2.39 15 . 3 2.88 2.93 .875
41
0.2
0 0~--~~--~~~~~~~~--~~~~~~~~----~----~~~~ 0°
0.005 0.02 0. 01 0.05 0.1 0.2 0.5 1 2 5 10
100~--------------------------------------------------------,
50
20
10
5
-;-
"S.
~
0.5
0 - 2~--~~~~~~~~~~~--~~~~~~~----~~--~~~~
0.005 0. 01 0.02 0.05 0. 1 0.2 0.5 1 2 5 10
s
42
Calculation No. 10
....,
~
50
20
10
5
....,
G
-Cxy • - C y x
1
0.5
0.2~~~~--~~~~~~~~~~--~~~~~----~--~~~~~
0.005 0.01 0.02 0.05 0.1 0.5 1 2 5 10
s
43
Calculation No. 11
5.25 .0313 79.0 .374 . 362 135 224 55.7 -120 23.0 314 101 -101 66 . 2
3.00 . 0547 79,0 • 375 . 362 76.7 129 30.2 -68.6 13.2 176 -57.8 -57.7 38.0
2.02 .0828 79.0 .375 .362 50 .7 87.4 19.3 -46 . 1 8.97 112 -38.7 -38.5 25 .5
1. 01 . 167 78.8 . 377 .360 25.8 43.9 9. 85 -23 . 1 4.61 61.9 -19.0 -18.9 13. 0
.515 .317 78.3 . 381 .354 13.2 24 .0 5.29 -11.8 2.60 30.2 -9.42 -9.38 6.86
• 304 .505 77.3 .394 . 341 8.01 15.5 4.07 -6. 77 1. 89 19. 7 - 5.27 -5.2 4 4.42
.201 .662 76.0 .402 .331 5 . 52 12 .3 3. 74 -4.61 1. 56 15.1 - 2.89 -2 .82 3.13
.0985 . 921 73.0 . 430 .296 3.01 10.3 4.0 1 - 1.86 1. 50 10.8 - .272 -.245 1. 95
.0524 1.06 67.6 .445 . 280 1. 79 11.6 4.79 - .23 1 1. 62 9. 72 .6 19 .632 1.28
.0100 1.13 47.1 .416 . 300 . 602 28 .0 9.23 2.89 2. 14 15. 6 2. 71 2.75 . 841
44
Calculation No. 11
1.2
1.0 1.0
goo
0.8 0. 8
75°
0.6 0. 6 60°
C1 w -e.
45°
0.4 0.4 Qs
0.2 0.2
0 0 oo
0.005 10
s
100~-----------------------------------------------------------,
50
20
10
5
-;
-s.
~
0.5
0.2~~~~--~~~~~~~~~--~~._~~~~~----~~~._~~~.
0.005 0.01 0.02 0.05 0.1 0.2 0.5 1 2 5 10
s
45
Calculation No. 11
100r---------------------------------------------~----------~
50
20
10
0.5
1 2 5 10
s
46
Calculation No. 12
5.06 .0227 76 . 9 .412 .282 154 474 72.0 -215 49.8 548 157 -157 92.5
3.01 .0383 76 . 9 . 412 .282 91.7 282 40 . 6 -128 29.7 319 -94.3 -94.1 54.4
2.04 .0571 76.9 . 412 .281 62.1 186 26.5 -86.7 18.2 215 -63.8 -63 . 6 36.9
.983 .108 76.8 .412 .281 30 . 0 93.7 12.7 -43.2 9. 79 106 -30.8 -30.7 17.9
.477 .254 76.8 .414 .279 14.7 46.6 6.83 -20.4 4. 79 50.9 -15.0 -15.0 8.77
.286 .395 76.4 . 417 .275 8.94 29.5 4 . 91 -12.3 3.24 33 . 0 -8 . 45 -8.42 5.45
.197 .544 75.8 .422 .270 6.27 22.4 3.88 -8.44 2. 35 24.7 -5.76 - 5.75 3.93
.0999 . 862 74.7 .434 .256 3.41 15.6 3.82 -3.75 1. 61 16.1 -2.13 -2.09 2. 24
.0504 1. 11 71.6 .456 . 231 1. 92 14.7 4.36 -1.11 1.40 12.3 .0850 .0901 1. 32
.0306 1. 22 68.0 .457 .228 1. 29 16.9 5.21 .0379 1. 49 11.2 .582 .586 . 978
.0199 1.28 63.9 .458 .226 . 958 20.1 6.32 1.08 1. 62 12.1 1. 24 1. 25 . 866
.00976 1. 30 57.0 .452 .230 .618 29.8 9.04 2.59 1. 93 15.8 2.63 2.66 .840
47
Calculation No. 12
1.4.-----------------------------------------------------------~
1.2
1.0 1.0
0.8 0.8
0-6
C1 w
0.4
0.2 0.2
0 0
0.005 0.01 0.02 0.05 0.1 0.2
s
100~--------------------------------------------------------~
50
20
10
-' 5
~
~
0.5
0 . 2~~~~--~--~~~~~----~--~~~~~----~--~._~~~
0.005 0.01 0.02 0.05 0.1 0.2 0.5 1 2 5 10
s
48
Calculation No. 12
20
10
:::;.
~
0.5
50
20
10
0.5
0·2~~~~~~--~~~~~~--~~~~~~~----~~~~~~~
0.005 0.01 0.02 0.1 0.2 0.5 1 2 5 10
s
49
Calculation No. 13
4.99 .Q272 91.5 .177 .541 99.0 140 77.4 -85.9 -2. 15 282 -106 -107 73.3
2.93 .0459 91.5 .177 .541 58.5 82.7 45.1 -50.0 -1.34 167 -62.6 -62.7 43.2
2.08 .0647 91.5 .177 .540 40.3 58.7 31.2 -35.4 -.823 117 -44.3 -44.5 30.7
• 961 .140 91.1 .180 .537 19.2 27 .3 15. 4 -16.4 -.101 55.2 -20.0 -20.0 14.8
.500 .265 89.9 • 190 .524 10.0 14.6 7.81 -8.61 .424 28.0 -9.58 -9.59 8.00
.305 . 393 87.8 .201 .507 6.42 9.53 5.36 -5.23 . 798 18.1 -5.10 -5.10 5.25
.207 .514 85.8 .218 .476 4.45 7.01 4. 27 -3.59 . 977 13.3 -2.55 -2.51 3.82
.104 • 712 78 . 6 .248 .423 2.39 5.03 3.40 -1.43 1. 22 8. 82 .0753 .0761 2.33
.0496 .885 67.4 .280 .364 1. 35 5.78 3.78 .0881 1. 41 7.60 1.12 1.14 1. 42
.0291 . 968 57.6 .267 .353 . 931 8.13 4.51 .775 1.60 8.44 1.84 1.87 1. 11
.0199 • 990 50.5 .260 . 357 .733 11.0 5.29 1. 52 1. 76 9.48 2.20 2.23 • 978
.00993 1.01 39.0 .229 .383 .509 19.1 7.42 3.20 2.09 11.9 2.63 2.65 .856
50
Calculation No. 13
1.2
1.0 1.0
goo
0.8 0.8
75°
0. 6 0.6 60°
w -e.
CJ
45°
0.4 0.4
30°
0. 2 0.2
15°
0 0 oo
0.005 0.01 0.02 0.05 0.1 0.2 0.5 1 2 10
s
100~------------------------------------------------------~----~
50
20
10
0 . 2~~~~--~~--~~~~--~~~~~~~~----~_.~~~~~
0.005 0. 01 0.02 0.05 0. 1 0.2 0.5 1 2 5 10
s
51
Calculation No. 13
....
::<
s
100~----------------------------------------------------------~
10
....
0
2
- Cxy•- Cyx
0.2~~~----~--~~~~~L--L~~~_._L~~~----~~--~~~_u
0.005 0.01 0.02 0.05 0.1 0.2 0.5 1 2 5 10
s
52
Calculation No. 14
4.90 .0172 81.2 .218 .379 133 391 117 -198 30.7 609 -197 -197 120
2.92 .0289 81.2 .218 .379 78.8 234 67.2 -118 18.2 361 -118 -117 70.5
1. 97 .0430 81.2 .218 .379 53.2 158 44.6 -79.6 12.3 241 -79.2 -79.2 47.3
• 956 .0898 81.2 . 219 . 378 25.5 77.3 20.9 -39.0 6.07 113 -39.4 -39.2 23.8
.521 . 165 81.0 .219 . 377 14.1 42. 3 12.5 -21.1 3.42 63.5 -20.7 -20.7 12.9
• 305 .277 80.7 .222 .372 8.26 25.6 7.40 -12.3 2.23 38.1 -12.2 -11.8 7.80
.192 .427 80.0 .228 .363 5.33 16.8 5.30 -7.80 1. 71 24.7 -7.00 -6.98 5.15
.100 .655 78.3 .239 .344 2. 93 10.7 3.82 -3 . 97 1. 34 15.0 - 2. 81 -2.79 3.01
.0513 .925 74.0 .256 . 311 1. 66 8.92 3.55 -1.61 1. 24 11. 1 - . 490 -.481 1.80
.0313 1.07 69.1 .263 .292 1.15 9.61 4.02 -.401 1. 36 10.0 . 33C . 345 1. 29
.0202 1.15 64.0 .265 .281 .865 11.7 4.80 .570 1.56 10.3 . 87 . 873 . 995
.0104 1. 15 53.9 .247 . 296 .559 20.2 6.60 2.24 1. 90 11.5 1.55 1.58 • 731
53
Calculation No. 14
1.2 €
1.0 1.0
0.8 0.8
0.6
0.4
0.2
100~------------------------------------------------~------~
10
5
-'
-s.
~
0.5
54
Calculation No. 14
...,
~
5
...,
0
-Cxy'- Cyx
55
Calculation No. 15
3.01 .0221 78.1 .234 .2 96 99.5 483 83.3 -213 44.9 610 -184 -184 101
1. 91 .0346 78.1 .234 .296 63 . 5 306 54.9 -135 28.6 391 -117 -116 64.1
. 996 . 0659 78.1 .234 .296 33.2 159 29.5 -70.5 14.9 206 -60.6 60.5 33.4
.503 .128 78.0 .235 .295 16.8 80.9 15 . 5 -35.6 7.61 106 -30.5 -30.4 16.9
. 313 .203 78.0 . 235 .294 10.6 50.9 10.1 -22.1 4.82 65.6 -18.9 -18.9 10.6
.208 . 303 77.8 .236 .292 7. 00 35.1 6. 60 -14.7 3.36 42.5 -12.1 -12.0 7.14
.0973 .584 76.7 .243 . 281 3.41 18.8 4.18 -6.85 1. 95 22.3 -5.30 -5.29 3.63
.0480 • 903 74.7 .252 .263 1.85 13.5 3.68 -3.01 1. 46 14.0 -1.54 -1.52 2.08
.0307 1.11 72.3 .263 .245 1.30 12.8 4.14 -1.53 1.38 12.6 -.299 -.265 1. 52
.0208 1. 22 69.2 .262 .239 . 976 13.8 4.72 -.559 1.45 11.8 .529 .532 1. 21
.00987 1. 34 62.1 .261 .233 .587 21.7 6.88 .209 1. 92 13.2 1.61 1.63 .836
56
Calculation No. 15
1.4
1.2
1.0 1.0
90"
0.8 0.8
0.4 0.4
0·2 0· 2
100~----------------------------------------------------------~
50
20
10
0 . 2~~~--~~~~~~~~~--~~~~~~~----~----~~. .~
0.005 0.01 0.02 0.05 0.1 0.2 0.5 1 2 5 10
s
57
Calculation No. 15
....
~
50
20
10
0 5
0.5
58
Calculation No. 16
4.70 .0500 55·7 0.126 0.760 150.5 24.8 31.6 -33.4 22.0 69.8 -0.572 -0.462 65.8
2.33 .100 57·3 0.127 0.760 75.0 13.9 15.8 -17.5 10.1 36.0 -0.844 -1. 36 31.9
1.115 .200 58.5 0.130 0.755 36.5 8.44 7.84 -8.72 4.94 19.2 -0.852 -0.816 15.4
0.695 .300 58.7 0.136 0.750 23.2 6.84 5·70 -5·76 2.30 14.3 -0.220 -0.364 9.26
0.475 .400 59.0 0.143 0.740 16.4 6.28 4.56 -4.03 2.76 11.6 0.109 0.125 6.58
0.33 .500 57·8 0.151 0.7)0 11. 9 6.28 4.24 -3.16 2.04 10.4 0.608 0 .232 4.22
0.233 .600 57.0 0.161 0.715 8.8 6.29 3·92 -2.32 2.10 9.34 0.728 0.740 4.08
0.159 .?00 55.2 0.171 0.705 6.45 7·67 3.88 -1.49 1.88 8.70 0.954 0.982 2.96
0.100 .Boo 52.2 0.185 0.685 4.45 8.76 4.18 -0.604 1.66 8.74 1.24 1. 24 1.99
0.0505 .900 46.6 0.209 0.655 2.65 11.92 5.30 0 .348 1.60 10.I,Q 1.82 1.83 1.23
0.0267 .950 39·7 0.227 o.63o 1.66 17.92 6.10 1.28 1.50 11.20 1.64 1.66 0.690
59
Calculation No. 16
1.0
90°
0.8
75°
0·6
OJ (I)
0.4
200~-------------------------------------------------.
100
60
Calculation No. 16
50~------------------------------------------------------~
100.---------------------------------------------------~
....
0
61
Calculation No. 17
ttJ
Boundary condition: Reynolds condition
Number of grid points: Axial direction: 10
(uniform spacing) Cir-
cumferential direction:
60
2.16 .0500 40.0 0.155 0.520 96.7 49.5 J6.0 -J9.J 44.7 84.0 -1.J8 -1.25 78.6
1.06 .100 42.0 0.155 0.520 46.7 26.J 17.7 -20.7 21.2 4).8 -1.66 -1.60 J8.7
0.510 .200 44.5 0.156 0.520 2).1 15.J 8.16 -10.5 10.0 22.8 -1.42 -1.40 17.4
0.)22 .)00 45.9 0.158 0.517 14.9 12.4 5.16 -7-59 6.)) 15.9 -1.57 -1.55 11.6
0.224 .400 47.5 0.161 0.513 10.6 10.9 4.59 -5.64 4.89 14.5 -0.9)0 -0.912 7.86
0.161 .500 48.1 0.16) 0.510 7-90 10.8 ).81 -4.74 ).81 12.) -0.921 -0.906 6.84
0.117 .600 48.5 0.167 0.50) 6.0J 10.9 ).60 -).75 J.21 11.J -O.J60 -0.525 4.47
o.o84o -700 48.5 0.171 0.500 4.60 11.7 ).42 -).06 2.84 10.6 -0.265 -0.253 4.05
0.0592 .Boo 48.5 0.175 o.49J 3.50 1) .1 ).54 -2.)0 2.5) 10.2 -0.004 -0.015~ ).27
0.0)79 .900 47.) 0.179 0.48) 2.50 15-7 4.08 -1.)) 2.16 10.5 0.)9) 0.405 2.29
0.0288 -950 4605 0.182 o.4Bo 2.04 17-9 4.59 -0.786 1.97 11.0 0.618 0.6)0 2.00
62
Calculation No. 17
1.0~------------------------------------------------~
2~--~---L~~~~~----~--~~~~~~----~~--~
0.01 0.2 0.5 2 5
s
63
50r-------------------------------------------..-----.
5
s
64
Calculation No. 18
1.27 .0500 31.0 0.169 o.t.oo 71.3 73-2 )6.1, -lt2.lt 61t.B BB.It -J.It5 -J.Jt 85.6
0.625 .too JJ-7 0.170 o.t.oo J5.J JB.It 18.2 -2J.It )0.9 t.B.o -2.52 -2.ltlt ltlt.O
0.299 .200 )6.6 0.170 0.)98 17.1 21.9 8.52 -12.5 llt.2 26.2 -2.16 -2.12 20.5
0.187 .)00 )8.0 0.171 0.)98 10.8 17.1 5-56 -8.56 8.92 19.2 -1.Bo -1.77 12.2
0.131 .~too t.o.G 0.172 0.)95 7-77 15.) lt.28 -6.88 6.1,8 16.0 -1 . 67 -1.67 9.0it
0.0950 .500 lt2.1 0;173 0.)95 5-85 11t.B ).It? -5.68 5.0lt 1J .Bit -1.52 -1.50 6.92
0.0700 .600 lt).2 0.17lt O.J9J lt.5J 15.8 J.OJ -lt.72 lt.2lt 12.5 -1.17 -1.29 5-20
0.0516 .700 t.t..o 0.176 0.)90 J.5J 16.0 2.9J -lt.Oit J-55 11.8 -0.9ltlt -0.972 lt.28
O.OJ71t .Boo 1,1,.7 0.178 O.J85 2-73 17.8 J.OJ -J.J2 ).10 11.lt -0.688 -0.668 ).56
0.0256 .900 1,1,.8 0.179 O.JBJ 2.06 20.8 J,l,t, -2.1,1 2.70 11.6 -0.298 -0.288 2 .58
o.o2o1t -950 1,1,.7 o.1Bo O.JBo 1.75 2).2 ).86 -1.92 2.1,1,1, 12.0 -0.128 -0.112 2.28
65
Calculation No. 18
66
Calculation No. 18
100~----------------------------------------~
....,
0
2
s
67
Calculation No. 19
0.856 .0500 27.0 0.178 0.)26 57-2 94.5 )9.6 -46.8 84.0 99.0 -4.16 -4.01 94.0
0 .42 .100 29.1 0.178 0.)26 28.2 50.0 18.4 -24.) 39 .4 50. 5 -).98 - ).90 45-7
0.202 .200 )2.0 0.178 0.)26 13-7 28.4 8.00 -1).9 18.1 27.0 -2. 79 -2-75 22.5
0.128 .)00 )4.4 0.179 0.)26 8.82 21.9 4.75 -9.60 11.7 19-5 -1.91 -1.88 1).2
0.0892 .400 )6.2 0.179 0 .)24 6.)0 19.4 4.05 -7.80 8.20 17-7 -1.96 -1.9) 9-55
0.0652 -500 )8.) 0.180 0.)24 4.76 18. 5 ).64 - 6 .60 6.)5 16.5 -1.72 -1.70 7-35
0.0484 .600 40.0 0.18o O.J22 J.70 18. 7 ).2) -5.65 5-15 15.1 -1. 50 -1.48 5-80
O.OJ57 -700 41.0 0.181 O.J20 2.88 19.8 ).08 -4.76 4.)8 11•• 4 -1.11 -1.09 4.27
0.0262 .Boo 42.1 0.182 0.)18 2.26 21.6 ).58 -4.16 J.76 15.1 -0.790 -0.775 ).64
0.0182 .900 42.6 0.182 0.)18 1.72 25.0 4.14 -J.J9 ).16 15-7 -0.)88 -O.J7J 2.89
0.0151 -950 4).2 0.18) 0.)16 1.50 27.) 4-57 -2.98 2.97 16.2 -0.175 -0.160 2.67
68
Calculation No. 19
69
Calculation No. 19
lOOr---------------------------------~~
70
Calculation No. 20
71
Calculation No. 20
72
Calculation No. 20
Kxx
Kxy
.....
0
73
Calculation No. 21
74
Calculation No. 21
75
Calculation No. 21
Kxx
Kxy
5
:;::.
::.:::
Kyy
2
1
-Kyx
.....
0
76
Calculation No. 22
77
Calculation No. 22
75°
60°
45°
30°
78
Calculation No. 22
5
:;:.
:.::
Kyy
2
....
c3
79
Calculation No. 23
80
Calculation No. 23
81
Calculation No. 23
82
Calculation No. 24
83
Calculation No. 24
1.4r-----------------------------------------------------,14
Of
50..---------------------------------,
5
s
84
Calculation No. 24
sor-----------------------------------------------------~
....
~
5
s
5
s
85
Calculation No. 25
86
Calculation No. 25
1.4r---------------------------------------------------~14
5or-----------------------------------------------------~
87
Calculation No. 25
50.------------------------------------------------------.
5
s
50.------------------------------------------------------.
88
Calculation No. 26
89
Calculation No. 26
I.4r-----------------------------------------------------~14
0.2
5
s
90
Calculation No. 26
50~------------------------------------------------------~
5
s
sor-------------------------------------------------------~
"-cxx=Cyy
...,
0
0.02 0.05 0. 1 5
s
91
Calculation No. 27
92
Calculation No. 27
1.4r-------------------------------------------------------~14
0.8
QJ
0.6
0.4
0.2
5
s
93
Calculation No. 27
sor---------------------------------------------------------,
94
Calculation No. 28
95
Calculation No. 28
t.4r------------------------------,I4
sor---------------------------~
96
Calculation No. 28
sor-----------------------------------------------------~
5
s
....
G
5
s
97
Calculation No. 29
98
Calculation No. 29
1.4 r - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - , 1 4
5
s
99
Calculation No. 29
sor---------------------------------------------------------.
sor-------------------------------------------------------~
100
Calculation No. 30
101
Calculation No. 30
sor-------------------------------------------------------
5
s
102
Calculation No. 30
sor-------------------------------------------------------~
5
s
50
....,
G
5
s
103
Calculation No. 31
104
Calculation No. 31
1.4~------------------------------------------------------~14
Qt
5or-------------------------------------------------------~
-
I
-s. 5
~
0.5
0.005 0.05 0. 1 0.2 0.5 1 2 5
s
105
Calculation No. 31
50~--------------------------------------------------------,
5
s
sor----------------------------------------- ----------------,
5
s
106
Calculation No. 32
107
Calculation No. 32
1.4.---------------------------------------------------------.14
5
s
108
Calculation No. 32
5Qr-------------------------------------- ----------------.
:::-
:::<
5
Kxx
5
s
sor---------------------------------------------------------,
.....
c3
5
s
109
Calculation No. 33
110
Calculation No. 33
1.0.-------------------------------------------------------~
¢= 0
t
g_l:-oo~s:-~-.~.-:!-l::-=-----=~--L---L..,...i,~u..Jo.J..-1___jo.L2---lL.-.J.....J-L....J....J....J.L____:r::::r:::::t::±s::::~:::t::t::.J...ll
s
111
Calculation No. 33
....
0
0.5
112
Calculation No. 34
113
Calculation No. 34
1.0
goo
75°
0.6
cy -e.
(<) Q,
~= 0
114
Calculation No. 34
10
s
115
Calculation No. 35
116
Calculation No. 35
1.0
90°
75°
0.6
o- -e.
w 45°
Q.
~= 0
~. LOO...I.5-L...I.~-:----:...l...:--'---..JW-LL.i..OU.-l--OL.2-.1..-~0L.5::i:I:J::I:t::==:t:2=::1--L.J.5...L.L..L.Jlt
s
0.2~~~~~----~--~~~~~-----'----L-~~~~----~--~~~~~
0. 005 0.01 0.02 0.05 0.1 0.2 1 10
s
117
Calculation No. 35
...,
:<
10
s
...,
G
0.005 0.01 0. 05 10
s
118
Calculation No. 36
Static Char act eris tics Spring Coefficients Damping Coeffici ents
s JJ.p-1
€ ¢ Qs Q, K xx K xy K yx K yy Cxx C xy Cyx C yy
119
Calculation No. 36
1.4~----------------------------------------------------~14
50
5
s
120
Calculation No. 36
sor---------------------------------------------------------,
....,
:<
5
s
sor---------------------------------------------------------,
....,
G
5
s
121
Calculation No. 37
122
Calculation No. 37
1.0
go•
0.8
75°
0.6
a -e.
(<)
0.4
~.~0~05~~~~~~-J~~UU~--~O.L2~~~~~~~==~z==~~±5~~
s
123
Calculation No. 37
10
s
124
Calculation No. 38
125
Calculation No. 38
l.Or--------------------------------------------
-----------------~
10
0.4
0. 2
126
Calculation No. 38
10
s
127
Calculation No. 39
128
Calculation No. 39
1.0~--------------------------------------------------------------,
0.8
0.6
w
0.4
0.2
129
Calculation No. 39
50
130
Calculation No. 40
131
Calculation No. 40
goo
75°
60°
-e.
45°
':E 0.2
s
0.2 2 5 10
s
132
Calculation No. 40
0.02
s
100-------------------------------------------------------------,
50
20
133
Calculation No. 41
134
Calculation No. 41
2.0 1.0
1.8 go•
1.6 0.8
75°
1.4
1.2 0.6 60°
(::§ 1.0 (<) -e.
45°
0.8 0. 4
0.6 30°
0.4 0.2
15°
0.2
0 0 o·
0.005 0.01 0.02 0.05 0.1 0.2 0.5 1 2 5 10
s
-
I
~ 5
~
0.5
0.2 0.5 1 2 5 10
s
135
Calculation No. 41
50
20
10
...,
~ 5
0.5
...,
c.> 5
Cyy
1 Cxy= 0
Cyx= 0
0.5
136
Calculation No. 42
137
Calculation No. 42
138
Calculation No. 42
100~--------------------------------------------------~-----------,
50
20
10
0.5
100~--------------------------------------------------~---------,
10
Cxx
::;:.
Cj
5
Cyy
I
Cxy= 0
Cyx= 0
139
Calculation No. 43
140
Calculation No. 43
1.0 1.0
go•
0.8 0.8
75°
0 0 o•
0.005 0.01 0.02 0.05 0. 1 0.2 0.5 1 2 5 10
s
141
Calculation No. 43
500
200
100
50
~
10
500
200
100
0
·~
50
Cxy= 0
20 Cyx= 0
10
Cxx
5
Cyy
2
0.005 0.01 0.02 0. 05 0.1 0.2 0.5 1 2 5 10
s
142
Calculation No. 44
143
Calculation No. 44
7.5
1.2
7.0
0.4
6.5
0.2
o1 0
0.2 2 5
so.-------------------------------------------------------~--~
20
10
5
-
I
~
~
0. 5
144
Calculation No. 44
50~------------------------------------------------------------~
20
10
5
...,
:<
0.5
0.2~~~~~~~----~--~~~~~U-----~--L-~~~~L-----L-~~~
0.002 0.005 0. 01 0.02 0. 05 0. 1 0. 2 0. 5 1 2 5
s
50r-----------------------------------------------------------~
20
Cxx
0.5
145
Calculation No. 45
146
Calculation No. 45
1.2.---------------------------------------------~
1.0
4.5 0.8
w 0.6
0.4
0.2
oL OL-~~~L---~~~-LLLLU----L--L-L~_u_u
0.005 0 . 02 0. 05 0. 2 0. 5 1
s
50,_----------------------------------------~
10
5
-I
-s.
~
0.5
147
Calculation No. 45
20
10
0.5
50~----------------------------------------------~
20
10
....
G 5
0. 5
148
Calculation No. 46
149
Calculation No. 46
1.0r-----~----------------------------------------,
3.0 0.8
z
0.6
C1 1<.1
0.4
2.5
0.2
ol 0
0.005 0.01 0. 02 0.05 0. 1 0.2 0. 5 1
s
50r---------------------------------------;
20
10
0.5
150
Calculation No. 46
50r-----------------------------------------~
50r-----------------------------------------~
20
10
·~
5
0
151
Calculation No. 47
152
Calculation No. 47
1.0
2.5 0.8
0.6
!':
cy w
0. 4
2 0.2
50.---------------------------------------
20
5
'7
-s.
~
0.5
153
Calculation No. 47
50r---------------------------------------- -.
0.5
50
' 20
10
5
::::;-
u
0.5
154
Calculation No. 48
155
Calculation No. 48
1.0 go•
0.8
so•
0.6
w '9..
0.4 Jo•
0.2
0.02 0. 2 0.5
s
0.005 0.05 1 2 5 10
156
Calculation No. 48
50 ~-------------------------------------------------------------------
20
10
5
...,
:.2
...,
G
157
Calculation No. 49
158
Calculation No. 49
1.0 r---------------------------------------------------------------,90°
oo
50r-----------------------------------------~-----------------,
20
10
-
I
5
~
~
0.5
159
Calculation No. 49
50 ~-------------------------------------------------------------,
20
10
0.5
5
:::o
\..)
0.5
160
Calculation No. 50
161
Calculation No. 50
l.Or---------------------------------------------------------------,900
60°
sor-----------------------------------------------~------------------~
0.5
162
Calculation No. 50
50~----------------------------------------------------------------~
.....
0
163
Calculation No. 51
164
Calculation No. 51
l.Or---------------------------------------------------------------~90"
60°
30°
o·
0.005 0. 01 0.02 0.2 1
165
Calculation No. 51
50~------------------------------------------------------------,
·~
166
Calculation No. 52
167
Calculation No. 52
8.0
1.0
7.5
0.8
w
0.6
7.0 0.4
100~----------------------------------------------------------~
50
20
0.5
168
Calculation No. 52
10~---------------------------------------------------------,
0.5
::;:.
::.::
0.2
0.1
0.05
5
s
..., 0.5
0
0.2
0.02 0.2 2 5
169
Calculation No. 53
170
Calculation No. 53
1.0
0.8
4.5
0.6
~
~ w
0.4
4.0 0.2
o'L
50~---------------------------------------
-----~
20
171
Calculation No. 53
SOr---------------------------------------------~
50~-----------------------------------------,
·~
172
Calculation No. 54
173
Calculation No. 54
1.0
0.8
3.0
0.6
~
a w
0.4
2.5 0.2
0.02 0. 2
s
.50r---------------------------------------------,
20
174
Calculation No. 54
50~------------------------------------~
20
0.02 0. 2 0.5
s
50r-------------------------------------,
20
·~
175
Calculation No. 55
176
Calculation No. 55
1.0
2.5 0.8
0.6
~
0 w
0. 4
2.0
0. 2
oL 0
0.02 0· 2 0.5
50~-----------------------------------
0.5
177
Calculation No. 55
50~------------------------------------~
20
~ 5
SOr------------------- -------------------,
20
178
Calculation No. 56
Note: Symbols not defined here are defined in the " Nomen-
HI cl a ture " section.
179
Calculation No. 56
l.Or-------------------------------~----,100
180
Calculation No. 56
Kxy =c 0
Kyx=c 0
0 .1.____---l~-L--L--L-L..ILWL...l..---..l------'-L__IL..L...J......l....U
0.01 0.02 0.05 0.1 0. 2 0.5 1
w
lOOr---------------------------~----~
Cxy=c 0
Cyx=c 0
·~
G Cxx
181
Calculation No. 57
5. Bearing geometry
s
Kxx Kxy Cxx Cxy
ez f·lj!-1
Kyy -Kyx Kzz Cyy -Cyx Czz
182
Calculation No. 57
Ir-----------~--------------------------------------~
-1
-2
1~----~~---L--~-L~-L~~------~~--~--~~~~~w
0.1 0.2 1 2 5 10
s
183
Calculation No. 57
-TI_______-~3~-------------.
lOr-----------~o~.rg_o,.~s--To_____
...,
G
184
Calculation No. 58
5. Bearing geometry
e, : eccentricity ratio
= 1 - h,/h0
cc groove angle
h0 : groove depth
y: land/groove width ratio
k: number of grooves
h,: bearing clearance
R: radius at the side of large
diameter
z p: half cone angle
I: axial length of cone
/ 9 : axial length of groove
185
Calculation No. 58
1~---------------------------------------------------,
0.5
-0.5
10
s
0. 9 0.5 0 -0 . 5 -1
50~~~----------~~------------~--------~r--------ft
186
Calculation No. 58
....,
1
::<
0.2
sr-~o.rg----------~~--------------~o__________-_o~-~5_______
-~I
....,
G
187
Calculation No. 59
5. Bearing geometry
(J)
t X
z
az
r=-
a,
e.: eccentricity ratio = 1 - h,/ h0 h,: bearing clearance
IX: groove angle Ri: inner radius
h0 : groove depth R.: radius of groove end
y: land/groove width ratio R 0 : outer radius
k: number of grooves
Spring Damping
Static Characteristics Coefficients
Coefficients
s
E:z ft/J-' Kzz Czz
.134 . 9 2. 6 2 1l . 3 340
7 . 05 - 5 12.3 . 1 55 .136
1 7 . 2 -1 23 . l .0515 .0286
188
Calculation No. 59
lr----------------------------------------------------,
0.5
.: 0
-0.5
0.1 20
s
189
Calculation No. 59
190
3. Experitnental Data
This chapter contains the following experimentally determined data for eleven sets of test bearings:
Circular bearing with two oil grooves 5 sets
Two-lobe bearing 3 sets
Tilting pad bearing 3 sets
Total 11 sets
The data was determined by means of eight different rigs. Test bearing sets 1 and 9 were tested on the same
rig; bearing sets 2, 7 and 10 were also tested on one rig. A different test rig was used for each of the other
bearing sets.
The data is listed in the order of test bearing no. as follows :
Table 3.o.1 Experimental results for oil film coefficients
Table 3.o.2 Test conditions
Fig. 3.o.1 Experimental results for oil film coefficients
Fig. 3.o.2 Cross section of bearing test rig
where o stands for the test bearing no.
In these figures, the experimentally determined values of the dimensionless oil film coefficients are shown
as functions of the Sommerfeld number, and the curves theoretically determined have been plotted and
drawn in as solid lines for comparison. Each experimental value can be dimensionalized by reference to its
test condition.
All formulations of bearing characteristics, static or dynamic, require evaluation of the pressure distribu-
tion of the oil film. Ideally, characteristics should be formulated in terms of parameters definable at the
design stage and easily measurable during service. Then, for a given pressure profile and, hence, a given
bearing load, the thickness and viscosity distributions of the oil film and the journal speed need to be known.
Details on these points are discussed in Sect. 4.1 with reference to calculation method and in Sect. 4.2 with
reference to measurement method.
191
Test Bearing No. 1
Test Sommer- Static characteristics Dimensionless stiffness coefficients Dimensionless damping coefficients
no. feld
number Eccen- Attitude Oil flow Temper-
tricity angle rate ature
ratio Kxx Kxy Kyx Kn cxx cxy cyx c yy
s e ¢ deg. q cm 3 js roc
1 0.025 83 45.7 7.522 1.553 0.755 1.074 - 0.459 -0.470 -0.360 0.291
2 0.035 83 47.0 7.974 2.218 1.334 1.212 0.222 -0.391 -0.479 0.356
3 0.037 83 40.9 8.191 2.260 1.185 1.231 0.322 -0.345 -0.369 0.391
4 0.052 83 40.9 8.639 3.083 1.927 1.530 1.507 -0.0020 -0.220 0.606
5 0.070 133 55.7 6.728 3.512 1.727 1.825 -0.167 -2.749 -0.837 -0.630
6 0.084 117 49.8 5.482 3.456 1.441 1.844 -0.085 -2.474 -0.882 -0.659
7 0.088 100 42.5 5.911 2.876 1.273 1.409 1.066 -0.610 -0.681 0.361
8 0.126 191 60.5 2.792 2.913 0.0064 1.741 3.795 0.288 1.322 1.261
9 0.133 83 42.3 4.050 3.237 0.431 1.733 3.958 1.042 1.081 1.827
10 0.158 167 54.3 2.761 3.160 0.131 1.798 3.314 -0.141 0.979 1.064
11 0.195 183 59.0 2.057 3.373 -0.295 1.611 4.854 0.887 1.466 2.317
12 0.195 167 53.2 2.672 3.274 -0.143 1.870 4.893 0.871 1.766 1.838
13 0.225 100 42.3 2.483 3.188 -0.235 0.755 4.529 0.829 0.374 2.397
14 0.284 167 54.8 2.972 4.056 -0.537 1.669 5.847 -0.885 1.654 2.877
15 0.336 150 58.0 1.753 3.584 -1 .059 1.392 7.020 1.857 1.605 4.015
16 0.442 133 51.3 2.934 4.624 -0.745 1.694 6.783 -1.244 0.951 3.179
17 0.498 167 55.4 7.971 9.659 -0.330 1.872 15.77 -1.377 -0.554 0.132
18 0.715 133 57.1 3.185 8.338 -1.701 2.551 16.90 -1.374 2.179 4.011
19 0.913 100 45.3 2.615 7.425 -2.075 2.339 17.63 2.599 0.794 5.058
20 1.458 133 56.8 2.615 14.60 -7.075 2.952 33.44 4.130 3.053 14.09
21 1.947 117 43.4 4.841 14.24 -7.170 2.711 36.40 4.740 -0.767 15.25
Notes:
• At a given steady load-speed combination, test vig temperatures were stabilized and data were separately obtained, first for the
determination of eccentricity ratio and attitude angle, and secondly for the determination of oil film coefficients. In order to examine
the effects of the thermal drift of proximity sensors and the bearing distortion, data for a complete check on the bearing centre was
also recorded after the dynamic tests.
• Determination of effective film temperature. The weighted average value (7;. + 2 I;,",)/3 is assumed to be the effective temperature,
where 7;. is inlet oil temperature and I;,", is outlet temperature of the lower oil pocket.
• Effective viscosity for Sommerfeld number is calculated from the effective film temperature.
192
Test Bearing No. 1
Notes:
• Properties of the lubricant
# 90 turbine oil
Specific gravity @ 15/4 °C 0.866
Kinetic viscosity mm 2 /s @ 37.8 °C 28.66
@ 50.0 °C 19.90
@ 98.9 °C 4.906
• Determination of amplitude ratio and dynamic load. The ratio of vibration amplitude to bearing clearance is adjusted within the
range 1/6 ~ 1/7 by changing the dynamic load 2 ~ 4 kN (peak).
193
Test Bearing No. 1
1.0 100°
0: £ }
Experiment L
~;.:.p
.8 80° I
!
i
{
.6 60°
!
"' ·=t=· ~
~c--.
.4 40° .... II
I
jjl
.2 200 rI
·0 .02 .05 .1 .2 .5
oo
2
Bearing material : White metal WJH
Backing metal S45C
a S=Y..!Y. ( R )'
P• "C; Bearing geometry
r-------------------------- --------0
35 f-.
20 f-.
0
0
0 Cxx
Uis
II
15 ,_
-Sf-
I I I I
.02 .05 .1 .2 .5 2 .02 .2 .5 2
0 0
s=Y..!Y.
Pm Cp
(..!i.)' c
s=Y..IY.(..!i
Pm Cp
. )'
b
Fig. 3.1.1 a-c. Experimental results for oil film coefficients. a Eccentricity ratios and attitude angles. b Dimensionless stiffness
coefficients. c Dimensionless damping coefficients
194
Test Bearing No. 1
195
Test Bearing No. 2
Test Sommer- Static characteristics Dimensionless stiffness coefficients Dimensionless damping coefficients
no. feld
number Eccen- Attitude Oil flow Temper-
tricity angle rate ature
ratio Kxx Kxy Kyx K,, cxx cxy c,x en
s e </> deg. q cm 3/s r oc
1 0.0393 0.95 24 70.0 55.0 12.5 2.61 2.52 1.92 2.23 0.099 0.101 0.290
2 0.0486 0.91 28 75.0 55.0 9.58 2.07 2.16 2.09 1.00 0.187 0.010 0.395
3 0.0588 0.80 29 103.3 64.2 5.86 1.48 1.13 1.48 1.18 0.355 0.085 0.419
4 0.0655 0.85 33 75.0 54.7 9.99 2.54 2.75 2.59 1.08 0.185 0.058 0.572
5 0.0720 0.73 31 130.0 70.3 4.95 1.60 0.843 1.46 1.68 0.544 0.176 0.476
6 0.0752 0.75 34 100.0 63.0 6.33 1.90 1.34 1.66 1.59 0.337 0.208 0.545
7 0.0779 0.64 27 146.7 78.0 4.76 1.81 0.933 1.49 i 2.11 0.740 0.226 0.497
8 0.0785 0.73 30 175.0 86.8 4.11 1.93 0.609 1.43 2.33 0.730 0.285 0.584
9 0.0925 0.68 36 130.0 69.4 5.20 2.03 0.890 1.70 2.31 0.528 0.392 0.663
10 0.1003 0.62 33 148.3 76.8 4.30 2.12 0.513 1.70 2.61 0.633 0.603 0.804
11 0.1023 0.66 37 183.3 85.2 4.13 2.25 0.524 1.60 2.93 0.763 0.585 0.828
12 0.1034 0.76 40 70.0 53.3 11.8 4.34 3.98 3.95 1.94 0.031 0.366 1.04
13 0.1066 0.68 40 96.7 61.5 6.99 2.65 1.65 2.09 2.43 0.267 0.481 0.825
14 0.1260 0.63 46 125.0 68.8 5.55 2.69 0.960 2.03 3.32 0.420 0.730 0.983
15 0.1354 0.53 44 146.7 76.4 5.00 2.86 0.745 2.06 3.77 0.531 0.700 1.01
16 0.1377 0.59 44 175.0 84.6 3.91 2.67 0.305 1.79 4.04 0.434 0.932 1.34
17 0.1611 0.61 58 91.7 60.6 7.42 3.54 1.72 2.54 3.57 0.031 0.880 1.43
18 0.1899 0.56 58 121,7 68.5 5.23 3.26 0.796 2.49 4.81 0.016 1.21 1.62
19 0.2032 0.45 57 141.7 75.9 5.11 3.00 1.01 2.06 4.85 0.319 0.990 1.27
20 0.2096 0.50 56 163.3 84.0 3.58 3.10 0.053 2.00 5.37 0.070 1.33 1.98
21 0.2131 0.57 53 65.0 51.7 13.3 7.07 5.02 6.73 4.02 0.481 1.12 2.28
22 0.2977 0.47 62 80.0 60.8 6.99 4.53 1.37 3.54 6.03 0.610 1.78 2.74
23 0.3552 0.44 76 96.7 70.0 5.32 4.17 0.233 3.19 7.90 0.627 2.12 3.73
24 0.3924 0.31 78 120.0 77.1 3.70 4.47 1.10 3.14 7.73 0.227 2.39 3.98
25 0.3950 0.37 73 130.0 85.0 3.21 3.83 1.79 2.69 7.82 0.606 1.93 5.07
26 0.4229 0.38 70 63.3 51.4 6.99 4.70 0.126 3.73 7.20 0.822 1.46 3.71
27 0.6143 0.33 75 70.0 62.7 8.73 5.86 0.344 4.84 10.5 1.04 2.33 5.70
28 0.6642 0.30 88 80.0 72.8 7.53 5.50 1.40 5.17 13.0 0.139 2.53 8.56
29 0.7358 0.17 90 91.7 79.6 6.07 7.46 4.00 4.77 11.2 0.366 2.09 9.07
30 0.7800 0.19 85 105.0 87.5 7.91 6.29 5.88 2.80 14.2 2.98 0.388 12.8
196
Test Bearing No. 2
Notes:
• Properties of lubricant
# 90 turbine oil
Specific gravity @ 20/4 °C 0.865
Kinetic viscosity mm 2 fs @ 40 °C 31.64
@ 50 °C 21.00
@ 10o oc 5.267
• Effective oil film temperature and equivalent viscosity for calculation of Sommerfeld number. The average bearing metal temper-
ature over 20 points is assumed to be the effective temperature of oil film, and the equivalent viscosity is calculated on the basis
of the effective temperature.
• Determination of exciting force. The ratio of vibration amplitude to bearing clearance is adjusted within the range from 1/6 to 1/20
by changing the exciting force (1 ~ 3.5 kN, peak).
197
100°
80°
60°
40° .....
20°
00
2
Bearing material: White metal ASTM B23
S=i!N(_B_)' Backing metal SS41
a Pm Cp
Bearing geometry
---Theory ---Theory
·~
<>
a
... ::;:.. c.ij:S:
II
Cxx
c.ijs ·~
0
II 0
0 G 0
0
:<
·~
0
0 0 10
Kxy cP
0
0
5 0 Cyy
Cxy, Cyx
0 0
0 Kyx
"'
-5 -5
b
s=lili (_B_ )'
P• Cp c s=iili(li_
pm Cp
r
Fig. 3.2.1 a-c. Experimental results for oil film coefficients. a Eccentricity ratios and attitude angles. b Dimensionless stiffness
coefficients. c Dimensionless damping coefficients
198
Test Bearing No. 2
199
Test Bearing No.3
1 0.0549 0.835 26.7 56.9 11.3 4.58 2.02 1.98 7.01 1.10 2.01 0.889
2 0.0854 0.793 32.4 55.8 7.80 4.67 1.12 2.11 8.35 2.46 3.08 1.82
3 0.114 0.755 35.6 55.6 5.47 4.10 1.795 2.23 8.18 2.82 3.32 2.24
4 0.123 0.749 29.6 64.4 5.97 4.24 0.493 1.95 7.29 1.72 2.54 1.54
5 0.182 0.685 32.7 62.3 3.98 3.44 0.0932 1.84 6.63 2.05 2.63 1.95
6 0.244 0.627 36.0 60.8 3.12 2.90 0.0989 1.69 5.94 1.87 2.18 1.80
7 0.318 0.572 39.5 60.1 2.53 2.65 0.125 1.48 5.66 1.79 1.71 1.64
8 0.408 0.521 43.7 59.3 1.99 2.62 0.194 1.37 5.77 2.08 1.40 1.61
9 0.501 0.467 49.4 58.5 1.60 2.45 0.203 1.33 5.66 2.08 1.22 1.55
10 0.564 0.464 62.1 62.0 2.43 2.04 -0.354 1.22 5.94 2.08 1.36 2.19
11 0.700 0.418 68.3 61.8 2.39 1.90 -0.534 1.10 6.01 1.66 1.34 2.51
12 0.871 0.372 75.8 61.0 2.48 1.83 -0.777 0.967 5.91 0.953 1.35 2.87
13 1.06 0.328 82.8 60.7 2.54 2.16 -1.01 0.811 6.03 0.640 1.34 3.45
14 1.31 0.287 90.4 60.8 3.12 2.40 -1 .45 0.763 4.81 -1.59 1.38 4.57
15 1.60 0.233 98.5 60.9 3.70 3.52 -1.80 0.671 7.92 - 2.77 0.988 5.76
16 2.21 0.157 108 61.3 4.17 6.03 -2.94 0.162 17.4 -3.38 -1.87 7.83
Test no. Sommerfeld Revolutional Frequency of Static load Bearing mean Supply oil Lubricant inlet Lubricant outlet
number speed dynamic load pressure pressure temperature temperature
s N rps wD/2n Hz WkN Pm MPa p 1 MPa T;n °C 7;,., oc
1 0.0549 26.7 45 4.50 0.901 0.1 50 54
2 0.0854 26.7 45 3.00 0.600 0.1 50 53
3 0.114 26.7 45 2.25 0.450 0.1 50 53
4 0.123 66.7 45 4.02 0.804 0.1 50 58
5 0.182 66.7 45 2.88 0.577 0.1 50 58
Notes:
• Properties of the lubricant
Turbine oil ISO VG32
Specific gravity @ 15/4 oc 0.862
Kinetic viscosity @ 37.8 oc 34.6
@ 50.0 oc 21 .3
@ 98.9 oc 5.50
• Determination of effective film temperature. The effective film temperature is assumed to be the average of four temperatures
measured at the locations shown in Fig. 3.3.2 by means of copper-constantan (CC) thermocouples implanted within 3 mm of the
bearing metal surfaces.
• Effective viscosity for Sommerfield number. Effective viscosity is calculated from the effective film temperature.
• Determination of amplitude and frequency of the dynamic load. Horizontal and vertical dynamic loads are set up so that the
amplitude of the vibration response in each direction reaches approximately one-tenth of the minimum bearing clearance it is or
taken as less than 300 N.
200
Test Bearing No. 3
1.0 100°
{>
.8
{> w
0 : € l 0
{>
80° 100
I
6 : ¢ J Experiment
.6 60°
IU Theory
"' C>
C> ~
.4 40°..,
"" ~~'*"~~-+ ~
~
~
.2 €
20°
Holes for insertion
.0
ao L=50 of thermocouples
.02 .05 .1 .2 .5 2
Bearing material : White metal WJ2
S=yjj_ (lL)' Backing metal SS41
a P• Cp
Bearing geometry
25
20
0 :
A
A
v
: Kxy
•
o : Kyy
Kxx
. Ky:r
l Experiment
20,------------------------------------.
15 --Theory
..,-
·~
~::: 10
II
·~
·::::=--- Cn
A
<>
A
A
- 5 '-----'---'--'-'--'-'-'-'-..____ __,___.__._._-'-'-'-'-'----'---'---'--' I I I I I I I
.01 .02 .05 .1 .2 .5 2 . 02 .05 .1 .2 .5 2
b
S= yJj_
P• Cp
(.Jl_)' c s"' yJj_ ( .Jl_ )'
pm Cp
Fig. 3.3.1 a-c. Experimental results for oil film coefficients. a Eccentricity ratios and attitude angles. b Dimensionless stiffness
coefficients. c Dimensionless damping coefficients
201
Test Bearing No. 3
202
Test Bearing No. 4
Notes:
• Properties of the lubricant
# 90 turbine oil
Specific gravity @ 15/4 °C 0.866
Kinetic viscosity mm 2 /s @ 37.8 oc 32.5
@ 50.0 oc 20.0
@ 98.9 oc 5.33
• The effective film temperature is assumed to the temperature at 3 mm below bearing metal surface.
• The effective viscosity for the Sommerfeld number is calculated from the effective film temperature.
• The exciting force range is set at 0.981 kN to 1.96 kN, and the ratios of vibration amplitude to the bearing clearance vary from 1/100
to 1/4.
203
Test Bearing No. 4
w
!
1.0 100°
0 0
.8 80°
50
o:e
.6 .:. : ; 60°
--Theory ....
.4 40°
.2 20°
00
0
.02 .05 .1 .2 .5 2 70
Bearing material : White metal WJ2
a
s=u.E. (.li. )'
pm Cp Backing metal SS41
Bearing geometry
---Theory ---Theory
....·~
.:,-
a
Js
II Js
·~
II
:.:: ·~
-5 -5
b c
s=u.E. (.li.)
p .. Cp
2
Fig. 3.4.1 a-c. Experimental results for oil film coefficients at N = 12.5 rps. a Eccentricity ratios and attitude angles. b Dimensionless
stiffness coefficients. c Dimensionless damping coefficients
204
Test Bearing No. 4
205
Test Bearing No. 4
2-1 0.0373 0.941 16.8 - 42.5 5.09 1.98 0.150 1.19 0.601 -0.106 0.120 0.235
2-2 0.0519 0.881 19.1 - 40.7 6.71 2.65 0.605 1.38 1.34 0.0313 0.440 0.384
2-3 0.0643 0.857 20.2 - - 6.13 2.26 0.215 1.30 1.63 0.156 0.399 0.475
2-4 0.0844 0.810 23.2 - 39.0 6.12 2.78 0.146 1.29 1.98 0.136 0.404 0.526
2-5 0.126 0.745 27.8 - 38.2 5.73 2.73 0.108 1.19 2.20 0.0846 0.390 0.626
2-6 0.164 0.690 32.0 - - 5.90 3.11 0.0879 1.24 2.74 0.241 0.657 0.865
2-7 0.239 0.603 41.3 - - 5.29 2.98 0.0137 1.16 3.04 0.458 0.796 1.14
2-8 0.308 0.548 45.5 - - 4.31 2.80 -0.0105 1.07 3.34 0.657 0.830 1.30
2-9 0.432 0.452 53.0 - 34.8 3.57 2.73 - 0.331 1.02 3.56 0.805 0.777 1.71
2-10 0.734 0.333 59.5 - - 3.42 3.20 -0.863 0.778 4.41 1.20 0.632 2.81
2-11 1.14 0.130 90.0 - 34.1 4.01 2.72 -1.68 0.859 5.69 -0.612 -0.0517 4.42
206
Test Bearing No. 4
1.0 100°
0
tJ.
.8
80°
0
0: €
.6 } Experiment
tJ.:; 60°
...
--Theory
.4 40°.....
0
.2 20°
0
0
.02 .05 .1 .2
oo
.5 2
0 : Kxx 6
. Kxy } Experiment
: 0 : Cxx A :
. C
Cxy }Experiment
0 : Kyx o . K>Y o·: Cyx 0 . yy
---Theory - - - Theory
.;,-
·~
a
U\s
·~
II
·~
~ G
.02 .2 . 02 2
b
S=Jdf_ (.ll.)'
P• Cp c
Fig. 3.4.3a-c. Experimental results for oil film coefficients at N = 25 rps. a Eccentricity ratios and attitude angles. b Dimensionless
stiffness coefficients. c Dimensionless damping coefficients
207
Test Bearing No. 4
Table 3.4.5. Experimental results for oil film coefficients at N = 50 rps in test bearing no. 4
Test Sommer- Static characteristics Dimensionless stiffness coefficients Dimensionless damping coefficients
no. feld
number Eccen- Attitude Oil flow Temper-
tricity angle rate ature
K xx Kxy K yx KYY cxx cxy cyx cyy
ratio
s e rjJ deg. q cm 3/s r oc
3-1 0.0511 0.822 21.0 - 52.0 4.98 2.18 1.19 1.66 0.714 -0.103 0.0571 0.162
3-2 0.0711 0.786 24.0 - 48.5 5.16 2.54 1.22 1.84 1.31 0.0686 0.351 0.392
3-3 0.0880 0.762 26.0 - - 5.54 2.40 0.289 1.09 1.53 -0.0294 0.296 0.474
3-4 0.116 0.714 29.5 - 44.8 4.63 2.23 0.519 0.972 1.44 -0.0495 0.0628 0.434
3-5 0.172 0.654 33.3 - - 4.46 2.44 0.644 0.897 1.98 0.0795 0.268 0.641
3-6 0.224 0.603 36.8 - 41.1 3.58 2.47 0.0896 0.698 2.58 0.384 0.455 0.921
3-7 0.327 0.532 41.4 - - 2.75 2.16 -0.134 0.435 3.14 0.610 0.639 1.26
3-8 0.421 0.476 49.2 - 39.5 2.33 2.28 -0.388 0.276 3.76 0.880 0.761 1.61
3-9 0.591 0.397 55.6 - - 1.75 2.31 -0.726 -0.173 3.95 1.04 0.583 2.15
3-10 1.01 0.278 69.0 - 39.8 1.67 2.75 -1.48 -0.912 5.28 0.867 0.357 3.33
208
Test Bearing No.4
1.0 ,-----------------------------------------,100°
.8
.6
--Theory
.4
.2
0
.02 .05 .1 .2 .5
S= 0Y_ ( _!i_ )'
a Pm Cp
.;,-
~~ a
U\s U\s
II
II
u
·~
:.;:
·~
Kxy
-5
.02
Fig. 3.4.4a-c. Experimental results for oil film coefficients at N = 50 rps. a Eccentricity ratios and attitude angles. b Dimensionless
stiffness coefficients. c Dimensionless damping coefficients
209
Test Bearing No. 5
Notes:
• Properties of the lubricant
* 90 turbine oil
Specific gravity @ 15/4 °C 0.871
Kinetic viscosity mm 2 js @ 30 °C 34.6
@ 50.0 oc 15.4
@ sooc 6.50
• Measurement method for inlet oil temperature ('I;") and outlet oil temperature (T,.,). 7;" and T0 . , are measured by using CC
thermocouples mounted in the inlet oil hole and the outlet oil hole.
• Effective viscosity for Sommerfeld number. The effective film viscosity is assumed to correspond to the temperature measured within
3 mm of the bearing metal surface in the direction of load.
210
Test Bearing No. 5
1.0 100°
.8 80°
0 :
e~t
£
} Experi m 60°
.6 "' :~ '/, '//#-
....
-
- - Theory
.4 40°
Hf~a-+~~a- q
Inlet hole for
supply lubricant
.2 20°
1:'1'1I J1IV1/.'n'
I
00
.02 .2 .5 2
Bearing material: White metal WJ2
a
S= y]1_ (-lL)'
P• C p Bearing geometry
Backing metal SS41
.:;
-~
a
Uj=::
II
J
-5 -5
b
S=MN(_B__)'
pm C p c
Fig. 3.5.1 a-c. Experimental results for oil film coefficients. a Eccentricity ratios and attitude angles. b Dimensionless stiffness
coefficients. c Dimensionless damping coefficients
211
Test Bearing No. 5
® ® ® ®
212
Test Bearing No. 6
Test Sommer- Static characteristics Dimensionless stiffness coefficients Dimensionless damping coefficients
no. feld
number Eccen- Attitude Oil flow Temper-
tricity angle rate ature
Kxx Kxy Kyx Kyy cxx cxy cyx Cy y
ratio
s e ¢ deg. q cm 3js roc
1 0.033 0.973 60.9 430 67.1 5.90 2.04 0.375 1.162 2.72 1.329 0.160 1.159
2 0.037 1.00 68.0 440 64.3 5.30 1.860 0.350 1.208 2.50 1.310 0.192 1.131
3 0.045 0.969 69.3 440 65.2 4.48 1.451 0.267 1.351 2.24 1.313 0.233 1.042
4 0.054 0.939 76.8 458 65.8 3.87 1.259 0.208 1.502 2.06 1.319 0.268 1.882
5 0.068 0.931 83.6 442 67.0 3.60 1.200 0.173 1.672 2.21 1.328 0.431 0.855
6 0.080 0.885 84.3 450 67.0 3.61 1.197 0.814 1.883 2.39 1.323 0.723 0.966
7 0.086 0.858 86.2 455 66.4 3.66 1.248 0.0388 2.06 2.47 1.330 0.851 0.953
8 0.096 0.862 88.2 435 63.0 3.73 1.324 0 2.27 2.65 1.353 1.00 0.971
9 0.110 0.836 90.3 450 59.5 3.84 1.436 - 0.0518 2.47 2.78 1.371 1.308 1.099
10 0.133 0.781 92.6 440 66.0 4.03 1.592 -0.201 2.81 2.85 1.401 1.596 1.124
Test no. Sommerfeld Revolutional Frequency of Static load Bearing mean Supply oil Lubricant inlet Lubricant outlet
number speed dynamic load pressure pressure temperature temperature
s N rps wnf2n Hz WkN Pm MPa p1 MPa T;. oc T,ut oC
1 0.033 60 36.40.45.51.4 388.5 2.51 -0.051 49.0 67.1
2 0.037 60 36.40.45.51.4 383.9 2.48 - 0.051 47.0 64.3
3 0.045 60 36.40.45.51.4 311.1 2.01 -0.054 47.8 65.2
4 0.054 60 36.40.45.51.4 251.9 1.63 -0.053 47.9 65.8
5 0.068 60 36.40.45.51.4 191.3 1.23 -0.065 47.9 67.0
6 0.080 60 36.40.45.51.4 163.1 1.06 -0.069 47.7 67.0
7 0.086 60 36.40.45.51.4 154.8 1.00 -0.071 47.7 66.4
8 0.096 60 36.40.45.51.4 158.3 1.02 -0.069 41.6 63.0
9 0.110 60 36.40.45.51.4 153.3 0.99 -0.065 40.3 59.5
10 0.133 60 36.40.45.51.4 97.1 0.63 -0.071 47.5 66.0
Notes:
• Properties of lubricant
Turbine oil IS032
Specific gravity @ 15/4 oc 0.868
Kinetic viscosity mm 2 js @ 37.8 °C 34.6
@ 50.00C 21.3
@ 98.9 °C 5.50
• Measurement method for inlet oil temperature (7;.) and outlet oil temperature (T,u1). 7;. and T,u, are measured by chromel-alumel
(CA) thermocouples attached to the inside of bearing casing and each pipe.
• Equivalent viscosity to calculate Sommerfeld number. T,u, is assumed as an effective temperature of oil film. Therefore, measured
oil pressure distribution corresponds to the calculated oil pressure distribution at T,u•·
• Condition of amplitude ratio and exciting force. Exciting force is 35.5 - 58.8 kN, so that the ratio of vibration amplitude to bearing
clearance will be 1/15 - 1/6.
213
Test Bearing No. 6
1.0 100°
0 ; E
A:{> } Experiment
.8 goo
--Theory
... .6 A 60°
....
.4 40°
.2 20°
0
.01 .02 .05 .1 .2 .5
Bearing material: White metal WJH
a
S=y_N (
pm Cp
Ji..)' Backing metal SS41
Bearing geometry
20 20,-----------------------------------~
10 f-
::;:-
~K.,
~
·~
"a
~:;:: ~:;::
II 5>-- Kxy I~ 5
~ 0
D D
G ~
Kyy
0
~ Kyx
-5 f-1 I -5
I I I
.01 .02 .05 .1 .2 .5 .08 .2 .5
b S=~.N( Cp
R)'
c
S=y_N
Pm
(.R.)'
Cp
Fig. 3.6.1 a-c. Experimental results for oil film coefficients. a Eccentricity ratios and attitude angles. b Dimensionless stiffness
coefficients. c Dimensionless damping coefficients
214
Test Bearing No. 6
215
Test Bearing No. 7
216
Test Bearing No. 7
Notes:
• Properties of the lubricant
# 90 turbine oil
Specific gravity @ 20/4 oC 0.865
Kinetic viscosity mm2 f s @ 40°C 31.64
@ 50°C 21.00
@ 100 °C 5.267
• Effective oil film temperature and equivalent viscosity for the calculation of Sommerfeld number. The average bearing metal
temperature over 20 points is assumed to be the effective temperature of the oil film, and the equivalent viscosity is calculated on
the basis of the effective temperature.
• Determination of exciting force. The ratio of vibration amplitude to bearing clearance is adjusted within the range of 1/6 to 1/20
by changing the exciting force (1 ~ 3.5 kN, peak).
217
Test Bearing No. 7
0
1.0 0 0 0
0 0 0 0 0 ~0 100°
0
00 cfJ
8CID
.8 ~ 80°
0 : <
} Experiment
b. : ~
.6 --Theory 60°
....
.4 40°
I>
.2 20°
0 00
.002 . 005 .02 . 05 .1 .2
Bearing material : White metal ASTM B23
S=y}l ( Jl.)' Backing metal SS41
a pm Cp
Bearing geometry
0 : Kxx b. : Kxy }
Experiment
<> : Kyx o :K,,
0 : Cxx
---Theory b. : Cxy } Experiment
() : Cyx 0 . c,,
---Theory
0
0 0
Cxx
0 ·~
0 a"
Uis
0 0
~~ 0
Uis 0
0
co 0
Kxx II ·~
0 u
:J
0
b. oS 0
Kyx
<> <>
e
-5
<>
.005 .02 .1 .2 .002 .01 .02 .1 .2
b
s=Mli(
Pm Cp
1l.)' c
S=y}l ( Jl.)'
P• Cp
Fig. 3.7.1 a-c. Experimental results for oil film coefficients. a Eccentricity ratios and attitude angles. b Dimensionless stiffness
coefficients. c Dimensionless damping coefficients
218
Test Bearing No. 8
Test Sommer- Static characteristics Dimensionless stiffness coefficients Dimensionless damping coefficients
no. feld
number Eccen- Attitude Oil flow Temper-
tricity angle rate ature
ratio Kxx Kxy Kyx Kyy cxx cxy c,x cyy
s e ¢ deg. q cm 3 /s T °C
Test no. Sommerfeld Revolutional Frequency of Static load Bearing mean Supply oil Lubricant inlet Lubricant outlet
number speed dynamic load pressure pressure temperature temperature
s N rps wv/2n Hz WkN Pm MPa p 1 MPa T;. oc 7;,., oc
1 0.005 60 80 18.865 4.911 0.102
2 0.008 59.5 80 13.720 3.571 0.102
3 0.013 58.5 80 10.290 2.679 0.102
4 0.029 58 80 5.145 1.339 0.102
5 0.046 59 80 3.430 0.893 0.102
6 O.Q78 59 80 2.058 0.536 0.102
7 0.148 60.5 80 1.029 0.268 0.102
Notes:
• Properties of the lubricant
Turbine oil ISO 32
Specific gravity @ 20/4 oc 0.87
Kinetic viscosity mm 2 /s @ 30 °C 43.7
@ 50 °C 19.2
@ 100 °C 6.09
• Equivalent viscosity to calculate Sommerfeld number. Viscosity at mean temperature of bearing surface was used.
• Eccentricity ratio and attitude angle are measured independently of oil film dynamic characteristics.
219
Test Bearing No. 8
1.6r-----------------------------------~
1.4
1.2
1.0 100°
0 : £ }
"' A : (> Experiment
80°
.8
-Theory
.6 60°
....
.4 40°
.2 20°
0 0
.002 .2
S=MN(Jl.)'
a P• Cp
25 25r-----------------------------------~
0 Kxz
0
0: Kxx 6: Kxy}
() : Kyz a .. K yy Experiment
- - - Theory ·~ 0 0
u
a
rj:s:
II
0: Cxx 6 : Cxy}Experiment
1:. ·~ () : Cyx D: Cyy
G
--Theory
Kxy Cyy
1:.
Kyy
()
D
0
Kyx -5
.002 .02 .2
S=P-N (JL)'
P. Cp
b c
Fig. 3.8.1a-c. Experimental results for oil film coefficients. a Eccentricity ratios and attitude angles. b Dimensionless stiffness
coefficients. c Dimensionless damping coefficients
220
Test Bearing No. 8
221
Test Bearing No. 9
Test Sommer- Static characteristics Dimensionless stiffness coefficients Dimensionless damping coefficients
no. feld
number Eccen- Attitude Oil flow Temper-
tricity angle rate ature
Kxx Kxy Kyx Kyy cxx Cxy Cyx Cyy
ratio
s e cjJ deg. q cm 3 js roc
1 0.033 800 59.3 2.598 -0.048 -0.223 2.149 0.306 0.095 -0.117 0.507
2 0.057 775.0 55.3 2.653 -0.0040 -0.182 2.211 0.382 0.168 -0.232 0.766
3 0.087 766.7 52.8 2.675 0.054 -0.211 2.245 0.362 0.078 -0.300 1.009
4 0.214 750.0 47.8 2.800 0.135 -0.202 2.378 1.061 -0.017 -0.403 1.672
5 0.305 716.7 46.3 2.717 0.189 -0.180 2.344 1.457 -0.143 -0.369 1.989
6 0.442 750.0 45.3 2.805 0.198 -0.277 2.391 1.915 -0.230 -0.288 2.442
7 0.620 750.0 46.3 2.933 0.159 -0.544 2.544 2.650 -0.346 -0.563 3.280
8 1.444 733.3 37.5 3.625 0.118 -0.957 3.135 3.980 - 0.423 -0.890 4.248
9 3.652 716.7 32.5 4.849 0.052 -1.383 4.189 5.982 -0.756 -1.288 5.795
Notes:
• Properties of the lubricant
# 90 turbine oil
Specific gravity @ 15/4 °C 0.866
Kinetic viscosity mm 2 js @ 37.8 °C 28.66
@ 50.0 °C 19.90
@ 98.9 °C 4.906
• Equivalent viscosity for Sommerfeld number. The effective film temperature is assumed to be the average of two temperatures
measured at points each 25° from vertical (x-axis) by means of CA (chromel-alumel) thermocouples 3 mm from the bearing metal
surface. Effective viscosity is calculated from the effective film temperature.
• Determination of the exciting force. The ratio of vibration amplitude to bearing clearance is adjusted within the range of 1/6 - 1j7
by changing the exciting force.
222
Test Bearing No. 9
r---{ IT·--
35 1- 35
20 f- 20
.:;
·~
... ·~ a
~~
~~ 15 f- II
II ·~
15
~
·~
0
10 - 10
51-
El
Kxx= Kyy 0 5
0
8 8 8~ 0 8 8
e e~-~
0 0
~
A
0 ~
0
0 ~ ij ~ ¢ ¢
~
0
-5 -
I I l_ I
.02 .05 .1 .2 .5 2 .02 .05 .1 .2 2
b
s=0:!..(ll.)'
P• Cp c s=0:!.. (.fl.)'
P• Cp
Fig. 3.9.1 b, c. Experimental results for oil film coefficients. b Dimensionless stiffness coefficients. cDimensionless damping coefficients
223
Test Bearing No. 10
Test Bearing No. 10: Tilting Bad Bearing with Five Pads
D = 100 mm, L =50 mm, CP = 0.216 mm, mP = 0.51, p= 60°,
r:t/p = 0.5, LOP
Test Sommer- Static characteristics Dimensionless stiffness coefficients Dimensionless damping coefficients
no. feld
number Eccen- Attitude Oil flow Temper-
tricity angle rate ature
ratio Kxx Kxy Kyx Kn cxx Cxy cyx cyy
s e cP deg. q cm 3 js T °C
1 0.029 0.85 533 53 .6 10.35 -0.17 -0.45 1.23 1.43 0.07 0.17 0.929
2 0.040 0.79 553 57.7 9.44 -0.22 -0.57 1.72 1.54 0.04 0.16 1.26
3 0.046 0.77 517 53 .0 11.16 -0.23 -0.84 1.70 2.25 0.12 0.30 1.34
4 0.048 0.71 558 59.3 9.17 -0.15 -0.59 2.00 1.27 0.03 0.05 1.41
5 0.052 0.68 558 62.4 9.09 -0.34 -0.46 2.12 1.55 -0.07 -0.41 1.83
6 0.052 0.61 553 70.4 9.03 -0.18 -0.54 2.91 1.32 1.06 -0.08 2.10
7 0.054 0.65 558 66.4 9.37 -0.40 -0.64 2.63 1.48 0.20 -0.32 2.05
8 0.065 0.71 517 51.0 9.95 -0.17 -0.88 2.17 2.54 0.16 -0.01 1.72
9 0.065 0.68 538 56.7 10.16 -0.31 -0.64 2.32 2.18 0.24 0.13 1.74
10 0.073 0.63 555 58.9 9.56 -0.39 -0.97 2.94 1.48 0.88 0.08 2.07
11 0.083 0.61 563 60.9 10.54 -0.69 -0.58 3.16 2.54 0.10 -0.58 2.51
12 0.088 0.62 533 54.4 10.00 -0.37 -0.82 2.99 2.64 0.26 0.02 2.48
13 0.089 0.58 557 63.4 10.24 -0.52 -0.58 3.66 2.56 0.74 -0.57 2.99
14 0.090 0.55 547 67.0 10.42 -1.41 -0.47 3.56 3.89 0.59 -0.45 3.09
15 0.098 0.62 517 51.7 10.12 -0.20 -0.79 2.94 3.56 0.40 0.15 2.39
16 0.106 0.58 553 56.9 9.52 -0.54 -0.73 3.28 2.25 0.27 -0.24 2.79
17 0.117 0.54 558 59.1 10.40 -0.40 -0.93 4.27 2.19 1.59 -0.28 3.35
18 0.126 0.52 555 61.5 12.00 -0.89 -1.29 5.20 3.37 1.87 -1.15 3.62
19 0.136 0.47 542 64.2 11.28 -0.95 -0.60 5.52 4.37 2.20 -1.49 4.12
20 0.142 0.53 600 54.2 9.88 -0.04 -0.10 4.20 3.92 0.62 -0.14 3.82
21 0.160 0.48 558 54.3 10.32 -0.87 -0.79 4.84 3.79 1.00 0.0 4.69
22 0.191 0.46 555 56.5 12.16 0.13 -1.23 5.76 2.96 1.39 -0.58 4.20
23 0.198 0.48 533 50.4 9.92 -0.15 -0.78 5.54 5.74 0.77 -0.02 4.27
24 0.207 0.40 525 61.6 13.36 -0.25 -0.72 7.20 6.77 2.43 -2.34 5.97
25 0.210 0.42 555 59.1 16.6 1.91 0.44 7.56 4.56 1.77 -2.81 4.87
26 0.298 0.39 508 52.2 12.00 -0.45 -0.48 7.97 6.50 1.58 -0.05 6.13
27 0.374 0.37 558 53.5 12.56 -0.76 -1.59 9.60 7.09 1.87 0.32 8.46
28 0.422 0.33 558 53.8 15.04 2.11 -2.73 11.92 7.11 2.99 - 0.70 8.71
29 0.463 0.31 558 56.1 18.40 1.68 -2.83 12.72 9.61 6.87 -2.21 9.43
30 0.494 0.29 537 57.6 17.76 5.41 -3.13 14.24 10.55 7.39 -4.09 8.89
224
Test Bearing No. 10
Test no. Sommerfeld Revolutional Frequency of Static load Bearing mean Supply oil Lubricant inlet Lubricant outlet
number speed dynamic load pressure pressure temperature temperature
s N rps Wv/2n Hz WkN Pm MPa p 1 MPa T;. oc T.ut oC
1 0.029 83.3 190 14.7 2.94 0.140 40.0 45.4
2 0.040 133 190 14.7 2.94 0.139 39.7 48.3
. 3 0.046 83.3 190 9.80 1.96 0.136 40.5 45.8
4 0.048 167 190 14.7 2.94 0.135 39.0 49.2
5 0.052 200 190 14.7 2.94 0.135 40.0 52.9
6 0.052 267 190 14.7 2.94 0.132 40.3 63.6
7 0.054 233 190 14.7 2.94 0.129 40.3 57.9
8 0.065 83.3 190 7.35 1.47 0.139 40.0 43.3
9 0.065 133 190 9.80 1.96 0.137 39.5 48.8
10 0.073 167 190 9.80 1.96 0.132 40.5 50.5
11 0.083 200 190 9.80 1.96 0.131 39.5 53.9
12 0.088 133 190 7.35 1.47 0.136 39.5 47.1
13 0.089 233 190 9.80 1.96 0.131 40.0 57.1
14 0.090 267 190 9.80 1.96 0.135 40.0 62.9
15 0.098 83.3 190 4.90 0.980 0.140 40.5 45.7
16 0.106 167 190 7.35 1.47 0.133 38.5 49.3
17 0.117 200 190 7.35 1.47 0.133 40.0 52.7
18 0.126 233 190 7.35 1.47 0.133 40.0 56.7
19 0.136 267 190 7.35 1.47 0.137 40.0 62.0
20 0.142 133 190 4.90 0.980 0.137 39.0 47.9
21 0.160 167 190 4.90 0.980 0.133 38.0 48.2
22 0.191 200 190 4.90 0.980 0.135 39.0 52.7
23 0.198 83.3 190 2.45 0.490 0.140 40.5 45.4
24 0.207 267 190 4.90 0.980 0.134 40.0 62.2
25 0.210 233 190 4.90 . 0.980 0.135 40.0 56.1
26 0.298 133 190 2.45 0.490 0.137 39.0 47.7
27 0.374 167 190 2.45 0.490 0.135 39.0 49.1
28 0.422 200 190 2.45 0.490 0.135 39.0 51.3
29 0.463 233 190 2.45 0.490 0.135 40.5 55.7
30 0.494 267 190 2.45 0.490 0.139 39.5 62.2
Notes:
• Properties of the lubricant
# 90 turbine oil
Specific gravity @ 20/4 °C 0.865
Kinetic viscosity mm 2 /s @ 40 °C 31.64
@ 50 °C 21.00
@ 100 °C 5.267
• Effective oil film temperature and equivalent viscosity for calculation of Sommerfeld number. The average temperature over five
pads, measured at the pivot position for each pad, is assumed to be the effective temperature of oil film, and the equivalent viscosity
is calculated on the basis of the effective temperature.
• Determination of exciting force. The ratio of vibration amplitude to bearing clearance is adjusted within the range of 1/6 to 1/20
by changing the exciting force.
225
Test Bearing No. 10
1.0 r - - - - - - - - - - - - - - - - - - - - ,
0 o Experiment
Theory
.4 ·+-
.2
...:::-
0
0 ..,.
·~
Uj:;::
0 a
II 0 Uj:;::
c II
:<
·~
0
oc
·~
0 0 c G
0
10
Kxx
.o..~~.
0
+.o.
<> <>
<>
Ooo <>
-5 <> <>
<>
.02 .05 .1 .2 .5 2 . 02 .2 .5 2
Fig. 3.10.1 a-c. Experimental results for oil film coefficients. a Eccentricity ratios and attitude angles. b Dimensionless stiffness
coefficients. c Dimensionless damping coefficients
226
Test Bearing No. 11
Test Bearing No. 11: Tilting Pad Bearing with Five Pads
D = 76.4 mm, L = 31.6 mm, CP = 0.086 mm, mP = 0, fJ = 58°,
a/fJ = 0.5, LOP
Test Sommer- Static characteristics Dimensionless stiffness coefficients Dimensionless damping coefficients
no. feld
number Eccen- Attitude Oil flow Temper-
tricity angle rate ature
ratio Kxx Kxy Kyx Kyy cxx cxy cyx cyy
s e ¢ deg. q cm 3 /s PC
1 0.106 0.88 0.6 168 42.0 3.03 -0.093 -0.353 0.395 0.264 -0.080 -0.173 0.344
2 0.144 0.85 1.0 169 41.5 3.38 -0.075 -0.320 0.414 0.234 -0.059 0.020 0.436
3 0.169 0.81 1.2 166 46.8 2.98 0.056 -0.025 0.300 0.286 -0.011 -0.118 0.592
4 0.224 0.74 0.0 167 47.0 3.34 0.058 0.074 0.284 0.769 -0.021 -0.052 0.716
5 0.306 0.69 -1.4 167 46.0 4.04 0.073 0.033 0.275 1.40 0.039 -0.197 1.14
6 0.496 0.65 1.5 167 45.8 4.28 -0.146 -0.226 0.177 2.89 0.007 0.446 1.61
7 0.708 0.54 2.1 167 46.0 3.72 0.494 0.336 0.553 2.95 0.095 0.048 1.90
8 0.786 0.57 1.6 166 44.5 3.75 -0.324 -0.269 0.242 3.66 -0.008 0.632 2.13
9 1.48 0.44 -2.2 166 45.8 2.79 -0.204 -0.088 0.361 3.73 -0.574 0.133 3.42
10 1.76 0.46 -5.5 164 41.5 2.72 -0.551 0.181 0.833 4.13 -0.772 -0.346 3.77
Test no. Sommerfeld Revolutional Frequency of Static load Bearing mean Supply oil Lubricant inlet Lubricant outlet
number speed dynamic load pressure pressure temperature temperature
s N rps Wn/2n Hz WkN Pm MPa p 1 MPa oc r..
7;,., oc
1 0.106 49.7 44.9 9.48 3.92 0.1 40.3 42.0
2 0.144 49.9 44.9 7.14 2.95 0.1 39.5 41.5
3 0.169 99.2 99.9 9.49 3.92 0.1 41.0 46.8
4 0.224 99.8 99.8 7.09 2.93 0.1 41.0 47.0
5 0.306 102 100 4.74 1.96 0.1 41.3 46.0
6 0.496 104 99.7 3.56 1.47 0.1 41.0 45.8
7 0.708 100 99.6 2.36 0.979 0.1 41.0 45.8
8 0.786 103 99.3 2.36 0.977 0.1 40.3 44.5
9 1.48 101 99.9 1.19 0.490 0.1 40.8 45 .8
10 1.76 102 100 1.20 0.495 0.1 38.5 41.5
Not es:
• Properties of the lubricant
# 140 turbine oil
Specific gravity @ 15/4 °C 0.874
Kinetic viscosity mm 2 /s @ 37.8 °C 61.8
@ 50.0 °C 35.4
@ 98.9 °C 7.78
• Measurement method for inlet oil temperature (T;.) and outlet oil temperature (7;,.,). T;. and 7;,., are measured by means of
thermocouples.
• Effective viscosity for Sommerfeld number. The effective film viscosity is assumed to have the value corresponding to the measured
outlet temperature.
• Determination of amplitude ratio and dynamic load. The dynamic loads are varied from 0.423 kN to 2.05 kN with an increase in
the static load. The ratio of vibration amplitude to bearing clearance is adjusted within the range 1 : 11 to 1 : 13.
227
Test Bearing No. 11
l.Or----------------------------------------
0.8c~ ~w
0 Experiment ' 0
0.61- ~ 0
- Theory o~
00
0.4t- ~
0.2
0 I I I 'I
0.02 0.05 0.1 0. 2 0.5 1 2
Bearing material : White metal WJ2
s=Mli (lL)' Backing metal S20C
a P• Cp
Bearing geometry
.:,-
·~
s
dl~
II
·~
Cxx
01-------
- 5 - 5
c
S= Mli
P• Cp
(lL) 2
Fig. 3.11.1 a-c. Experimental results for oil film coefficients. a Eccentricity ratios. b Dimensionless stiffness coefficients. c Dimension-
less damping coefficients
228
Test Bearing No. 11
229
4. Explanations
(/l
With the following dimensionless variables
ap = ~ au) ap- a (flaw) (1) f) = x jr, 1J = yjh, ( = zjr
ax ay ay ' az - ay ay
u = uj(rw), v = vj(rw), w= wj(rw)
0
ap
u= -
(yJ-Y dy - J0h-flY dy J0 -dyfl j J0h -dy)
y
where
ax 0 fl fl
y dylh dy
+ Uh· J- J -, (3) (7)
ofl ofl
w = ap
az
(s l:
0 fl
dy -1 ydy
0 fl
J dy/1 dy)
0 fl 0 fl
0
231
4 Explanations
Eq. (10):
0
a ( 3 aA) a ( 3 oA) aH
aB H aB + ac H ac = 6 ao · (11)
where rx = f1
au
ay- QU V'
I I
(12)
e = ejC: eccentricity ratio
</>: attitude angle
rz = f1
aw
ay - (! v w'
I I
Substituting Eq. (8) into Eq. (5), with some omis- the turbulent flow and Q is the density of the lubri-
sion we get the following equation (If/ ~ 1 is taken cant. The second term on the right-hand side of
into consideration): Eq. (12) is called the Reynolds stress. Shearing
stress for the laminar flow is obtained by setting
3 3
Reynolds stress is given by the following equa-
tion based on the mixing length model of Prandtl
= :() (H G) + e</> sin 0 + el cos 0.
1 (9) [4]. Flow in the sliding direction is assumed to be
predominant
The Reynolds equation for the isoviscous case
is obtained by setting E = 1/12 and G = 1/2 in - eu v I I = e[2 . 'au' au
oy . oy,
(13)
Eq. (9):
I
ay . aw
- ev w = e[2 . 'awl ay,
I
232
4.1 Method of Calculation for Bearing Characteristics
where Gx and Gz are called the turbulent flow coef- Applying the transform by Eq. (19) into Eq. (10),
ficients and are the correctional terms of viscosity the following equation is obtained :
caused by the turbulent diffusion.
The turbulent flow coefficient is given by the fol- az A a2 A
-2+ - + C1 -
aA + C -aA + C A- =C .(20)
lowing form by using the local Reynolds number, ae a(
2 ae 2 o( 3 4
Re1:
1 Generally, coefficients C 1 and C 2 are not equal to
zero. cl and c2 are zero only in the isoviscous case.
Gx = 12(1 + rxxRefx}'
(17) C 3 and C4 are given by Eq. (21) for the isoviscous
1 case:
Gz = 12(1 + rxzRef'} ·
233
4 Explanations
x.
Azj
"''
(n)
j M-1 N-1
L L
No
1- +'~) < J c • (26)
j- 1
"'' i=2 j= 2
(n
Aii
i- 1 i+l
t= 0 The allowable error, J 0 , must by very small, in the
~--------------------~---8 order of 10- 6 "' 10 -!2, to ensure a good level of
8= 0 8= 2 1r
Fig. 4.1.2. Finite-difference mesh accuracy in the calculation.
234
4.1 Method of Calculation for Bearing Characteristics
p~----------------~h
L! D=O. 7
e=O.S
:
+
:z::
II
I'<:
"....
~
"'0.
....
0
8
e -0 . 1 Circumferential
"'"' angle 8
a Sommerfeld boundary condition "'
"2
0
-0.2
·;:;; -0 .3
p~----------------,h c
8"' -0.4
6
-0.5
-0.6 Isoviscous theory
-0.7
:o ( H3 ~~) + :( ( H3 ~~)
= - 6e(1 - 2¢') sinO- 12e' cosO. (29)
c Reynolds boundary condition The first term on the right-hand side of Eq. (29)
represents pressure generation by the wedge film
Fig. 4.1.3a-c. Pressure boundary conditions for fully circular
journal bearings effect and the second term, pressure generation by
the squeeze film effect.
For A 1 , which is the solution of Eq. (29) with
yields the following conditions: e' = 0 and 1 - 2¢' = 1, the Reynolds boundary
condition [Eq. (27)] is used. The following boundary
A =0 at 8=0 condition is used for A 2 , the solution of Eq. (29)
oA withe'= 1 and 1 - 2¢' = 0:
A= 0, ao = 0 at 8 = 0* (27)
A =0 at ' = 0, (* oA2 = o at 0 = 0,2n
ao
where (* is the dimensionless bearing width and is
defined as (* = Ljr. (30)
The film pressure obtained using the Reynolds The film pressure for the case that the two terms
condition is positive in the region ofO < 0 < 0* and of the right-hand side in Eq. (29) exist simulta-
zero in 0* < 0 < 2 n. 0* is not the prescribed value, neously can be expressed as follows:
but is determined during the calculation. 0* is a
function of ( . If the axes of the journal and the (31)
bearing are parallel to each other, 0* yields a con-
Some conditions are proposed for Eq. (31):
vex and symmetrical curve in the axial direction.
Figure 4.1.4 shows the process used to determine (i) A 2 = 0 at the location of A 1 = 0 (0 = 0, 0*)
0*. (Fig. 4.1.5 a).
Next, the boundary condition in the unsteady (ii) A 2 exists in the region of 0 < 0 < 0* and
state is investigated. In this case, two boundary A 2 = 0 at 0 = O* (Fig. 4.1.5b).
235
4 Explanations
0/J
d Condition (iv)
a) Bearing Load-Carrying Capacity The friction coefficient is defined as the ratio of the
and Attitude Angle friction force to the bearing load. It is given by the
following equation:
The bearing load-carrying capacity is defined as the
reaction force of the lubricating film. Therefore, it h If/ - 1 = n: SS(_!_ auI ) dO dC ,
can be obtained by integrating the film pressure ' A H a17 ~=1 (36)
calculated from the static Reynolds equation
[Eq. (11 )]. The dimensionless film reaction force, F, JB If/- 1 = n: SS(_!_ au I ) dOd' '
and its components with respect to the eccentricity ' A H a17 ~=0
direction and its perpendicular direction, .F:, Fq, , are where h and h are the friction coefficients on the
236
4.1 Method of Calculation for Bearing Characteristics
journal and bearing surfaces, respectively. oii/81'/ film reaction force, F, acting on the journal centre,
can be calculated as follows : 0 1 , can be expressed as a linear function of the
au 2 aA (
a11 = H ae 17 - 2 + 1 ·
1) (37)
displacement and its velocity. Figure 4.1.7 shows
the Cartesian coordinate system, OB- xy, with the
origin of the bearing centre, OB. When F is ex-
panded in Taylor's series and the terms of second
The friction coefficient calculated in this data-
and higher powers are disregarded, the following
book is that on the journal surface. It is calculated
equation is obtained:
from the first part of Eq. (36).
II
equations: F0 x, F0 Y: components of F0 in x and y direc-
tions at the static equilibrium
Qs = (L/D) !(Hw){=od1'Jdei, (38) point 0 10
{* 1 ·= dfdt: differentiation by time t
Qe = (L/D) J0 0J(Hii)o=o~d1JdC (39) iJ Fx/ ox lo etc.: partial differential coefficients at
= = =
where (* L fr, ii uf(rw), w w/(rw). (}~ repre-
sents the location of Reynolds boundary at ( = 0.
static equilibrium point
The second and third terms on the right-hand
side of Eq. (40) are called the stiffness and damping
matrices k, c of the oil film, respectively. Then, k and
c can be written as follows:
4.1.6 Method of Calculation for iJFx
Dimensionless Stiffness ay
and Damping Coefficients of Oil Film k=
iJFY iJFY
ax ay o (41)
c
The stiffness and damping coefficients of the lubri-
cating film were defined in Sect. 1.3.1. In this sec- iJFx
-
tion, the theoretical procedures for calculation of ax ay
c=
these coefficients are explained. iJFY iJFY
If the journal vibration about its static equilibri- ox ay o
um point, 0 1 o, is very small, the nonlinear dynamic
The dimensionless stiffness and damping coeffi-
cients are commonly defined by the following equa-
=(K
tion:
K xx Kxy) = CP k
K yx Kyy W ' (42)
C (Cxx
= Cxy) = Cvw c.
Cyx Cyy W
When these coefficients are applied to Eq. (40), the
following dimensionless equation is obtained:
237
4 Explanations
+ {(ar.l
o<f> + F_q,) sm</>
. (aF"' o<f> I F_) sm</>
. }L1</>
as described in Sect. 4.1 .5, it is convenient to express
t_!le left-hand side of Eq. (43) by using£., Fq,, F0 , and 0
0 0-
0
- 0, 0
Foq,· Then the left-hand side is rewritten using these
+ (ar.j . . ~. aF"'I
oel cos ..~.) L1e
four terms, which are proportional to the small
increments of displacement from OJo• that is, oel sm 0 '1'0 - 0 '1'0
I
Substituting Eqs. (44), (45) and (46) into the left- X (oF"''
oe 0- 1 oF:'
eo o</> 0 ) -eoFFo'
,
0
238
4.1 Method of Calculation for Bearing Characteristics
1 1
0.5,L-~~~~~~~~~~~~~--~=-~~~~·
Fig. 4.1.8. Partial derivatives of the film reaction
0.05 0.5
force at the static equilibrium point Perturbation Je, J¢. Je,' J¢'
239
4 Explanations
where F, and F, are the components of F in the References for Sect. 4.1
OB- uv system. The film spring and damping coef-
ficients of the OB- uv system can be calculated [1] Dowson, D., Int. J. Mech. Sci., Vol. 4 (1962), 159.
from those of the OB - x y system by applying the [2] Pinkus, 0., Sternlicht, B., Theory of Hydrodynamic Lubri-
following conversion: cation, (1961), McGraw-Hill.
[3] Hinze, J.O., Turbulence, (1975), McGraw-Hill.
r
cos 2 fJ sinfJcosfJ sinfJcosfJ sin 2 fJ ~ 348.
*- -sinfJcosfJ cos 2 B -sin2 B sinBcosB [8] Wada, S. and Hashimoto, H., Journal of JSLE, Vol. 44,
C - - sinBcosB - sin 2 B cos 2 B sinO cosO No. 382, (1978), 2140.
sin 2 8 -sinO cos() - sinBcosB cos 2 8 [9] Elrod, H.G., Ng, C.W., Trans. ASME, Ser. F, Vol. 89,
No. 3, (1967), 346.
(54) [10] Vogelpohl, G., VDI-Forschungsheft, 386, (1937).
[11] Akasaka, T., Numerical Calculation, (1977), Korona-sha.
Theoretical procedures that can be used to ob- [12] Glienicke, J., Dissertation, Karlsruhe (1966).
tain the bearing characteristics have been briefly (13] Someya, T., Dissertation of the University of Tokyo,
described. These procedures or methods can be ap- (1972).
plied for all journal bearings. Concerning the film [14] Nakanishi, Y., Chikira, S., et al., Hitachi Zosen Technical
thickness function H, although only one case of the Review, Vol. 37, No.1, (1976), 49.
[15] Dandage, S.R., Dissertation, Univ. of Wisconsin, (1974).
circular journal bearing is shown here, H takes var- [16] Pafelias, T.A., Dissertation, Rensselaer Polytechnic Institu-
ious function forms according to the bearing types. te, (1974).
In addition, in the case of pivoted pad bearings, the
balance of the moment of each pad about the pivot
point should be also taken into consideration in
calculation of the film reaction force.
240
4.2 Methods of Measurement for Journal Bearing Characteristics
Dynamic loads act on the bearing to greater or centre and the reaction forces of the bearing, i.e., the
lesser degrees under real operating conditions. The dynamic characteristics of the bearing. Devices that
operating conditions can be classified in two broad can be used to apply the dynamic load include elec-
groups: those in which the journal centre moves all trodynamic, electrohydraulic, unbalance weight-
around the bearing clearance, as in the case of bear- type and cam-type exciters. The test apparatus
ings in reciprocating engines; and those in which shown in Fig. 4.2.2 has two electrohydraulic ex-
the journal centre moves with a small displacement citers which act on the test bearing casing with
around its static position of equilibrium. small movements. Details of the apparatus are ex-
The dynamic load is required to measure the plained in Sect. 4.2.3.
relation between those movements of the journal
Test bearing
b) Measurements of Friction Torque
Support bearing
In a case where a test bearing rests freely on its
shaft, the friction torque of the bearing can be ob-
tained by measuring the reaction force of an arm
overhung from the bearing. In other cases, the
torque can be measured directly by a rotary torque
Test shaft
meter or a torsion shaft coupled into the drive train
between the driving and the driven shafts. Other-
Static load
wise, in the case of bearings in actual machines, the
frictional loss is usually estimated from the oil feed
Hydraulic cylinder rate and the temperature rise.
for static loading The frictional coefficient of a well-lubricated
bearing is normally in the order of 0.001, and the
friction force is much smaller than the static load.
/
Consequently, to obtain better accuracy of the fric-
Fig. 4.2.1. An example of a static loading mechanism in a bear- tion torque, losses other than those in the test bear-
ing test rig ing should be measured in advance.
(D Test bearing casing (J) Tension rod Shaft diameter Max. 280 mm
0 Bellows cylinder ® Support ball bearing regular size 140 mm
(for static loading) I00 ~ 10000 rpm
® Flexible diaphragm coupling Speed range
0 Displacement sensor
® Rotational speed sensor Static load Max. 50kN
@ Accelerometer from three directions
® Test shaft
Dynamic load Max. 10 kN, peak sinusoidal
(i) Load cell
@ Test bearing DC ~soo Hz
@ Hydraulic exciter head @ Flexible bar
(for dynamic loading) (for prevention of side inclination)
Fig. 4.2.2. An example of a test rig for dynamic characteristics of journal bearings - cross section of the main parts
241
4 Explanations
c) Measurements Recorded and Data Processing Hagg and Sankey [1] derived the four coefficients
which disregarded the effect of cross-coupling be-
With static characteristics the following variables tween the major and minor axes of the journal lo-
are measured in addition to the friction torque: for cus. However, to eliminate the four cross-coupling
control purposes the rotational speed, steady state terms and to reduce the original eight coefficients to
load, oil feed pressure, oil feed rate and oil inlet four equivalent coefficients with the nature of the
temperature, for example; and for experimental bearing dynamic characteristics, each vibration
data, such variables as the steady-state locus of the phenomenon has to be analysed individually, as
shaft, the oil film pressure and temperature distri- discussed elsewhere [4].
bution in the bearing. In dynamic testing, such dy- In the method adopted in [3] and [7], after the
namic high-frequency magnitudes as the exciting four stiffness coefficients had been determined from
forces and the resultant displacements, velocities the static equilibrium locus of the journal, the resid-
and accelerations must be measured in addition. ual four damping coefficients were obtained from
These numerous data are processed by computer. the dynamic response of the journal to unbalance
The practical example is described in Sect. 4.2.3 g. excitation.
242
4.2 Methods of Measurement for Journal Bearing Characteristics
To obtain the eight oil film coefficients, two or more Fig. 4.2.3. Photograph of test rig for dynamic characteristics of
separate and different loadings must be applied in journal bearings
the above methods, for which it is essential to
ensure that the duty parameter for the bearing
remains constant throughout the test sequence.
The methods in group c can be modified to ob-
tain the eight coefficients at once, by simultaneous test bearing ©, which is set in the test bearing
application of sinusoidal forces in two directions casing Cl), floats on its shaft. Three sets of com-
with different frequencies and measurement of the pressed-air bellows (blare attached to the test bear-
resulting motions in two directions. This is the two- ing casing horizontally and vertically. The bellows
directional compound sinusoidal wave excitation are used for static loading of the test bearing and
method [10, 13, 16], which was adopted for test are so flexible that the test bearing casing can orbit
bearings 1, 5 and 9 in Chap. 3. freely around the shaft. The concentric setting and
Details of the analytic procedure for extraction of operation of the test bearing is secured by means of
the eight coefficients by the methods in groups c an isotropic arrangement of eight flexible bars @.
and d can be found elsewhere [10]. In addition, two hydraulic exciters are located
left and right of the test bearing casing at 45° from
vertical.
4.2.3 Example of a Test Rig for Dynamic Each exciter head @ is connected to the test
bearing casing by means of the tension rod (/), and
Characteristics of Journal Bearings the crevice joints at both attachment points, to the
test bearing casing and the rig bed. These exciter
The test rig developed by Glienicke [5] for the mea- connectors are able to impose negligible constraint
surement of oil film coefficients is frequently quoted on the test bearing casing.
as a typical one, which has been followed by many Since each of the static loads provided by three
other apparatuses, refined mainly in connection bellows and each of the bias loads from two exciters
with the loading mechanism and the processing of is independently adjustable, the static load vector
the measurements. can rotate through 360°. This function is adequate
Eight test rigs are shown in Chap. 3, with the for bearings with a noncircular cylindrical bore,
experimental results obtained with each. In the which the equilibrium curve itself is dependent on
paragraphs below a representative survey is made the angular position of the static load.
of the rig shown in Fig. 3.1 .2, in which the com- The two exciters make it possible to apply dy-
pound sinusoidal wave excitation method can be namic sinusoidal forces together with any relative
applied. Figure 4.2.2 shows the configuration of phase and magnitude. The ratio between the two
the principal mechanical components of the rig. exciting frequencies can also be adjusted to any
Figure 4.2.3 is a photograph of the rig. ratio with an integer up to 20 (e.g. 1 : 2, 3 : 4).
The test shaft (j]) is supported at both ends in Provision is made for the following measure-
comparatively rigid-support ball bearings ®. The ments:
243
4 Explanations
a) Static Displacement and Vibration Response meter of pulse generating type through an FV (fre-
quency to voltage) converter.
The relative displacements between test bearing
and journal are measured by four pairs of contact-
f) Measurement from Inside Rotating Shaft
less eddy current proximity sensors Q), located on
each side of the bearing in the horizontal and verti- The thickness and pressure distribution of the oil
cal directions. film around the journal surface provide useful infor-
Eddy current proximity sensors offer the advan- mation on the operating conditions of the bearing.
tages of noncontacting, compact and lightweight In order to measure these variables by minimum
design and ensure broad ranges of frequency re- transducers, eddy current proximity sensors,
sponse and amplitude linearity. Some of them, how- piezoresistive pressure transducers and CA thermo-
ever, are likely to cause zero drift with temperature couples are incorporated in the rotating shaft.
changes, which must be compensated for static dis- These transducers are very small and light. Their
placement measurements. outputs are directed through a slip ring.
The inertia of the test bearing and its associated
structure bypasses a proportion of the forcing, so
g) Data Processing
that the true force across the oil film differs from the
exciter output force. To compensate this effect by Figure 4.2.4 is a schematic of the data processing
means of mass cancellation, the absolute accelera- system with a minicomputer. For determination of
tions of the test bearing casing are measured by four test conditions or for control purposes, low-speed
pairs of piezoelectric accelerometers @, each locat- data such as the rotational speed, static load, oil
ed at the position of a proximity sensor. feed pressure, oil feed rate, oil inlet temperature,
and also certain significant temperatures in the
b) Static and Dynamic Loads bearing are fed into the minicomputer through an
analog-to-digital converter, digital voltmeter and
The static loading applied by the compressed-air scanner.
bellows is measured by Bourdon-type pressure The exciting forces and the resultant displace-
gauges. ments, velocities and accelerations from which the
The static bias load and dynamic load applied by oil film coefficients are determined are high-speed
the hydraulic exciters are measured by the load cells data. These high-speed data are also fed into the
G) attached to the exciter heads. The total static minicomputer simultaneously with the low-speed
force vector is the sum of the loads applied by three data via a high speed A/D converter and a sample
bellows, two exciters and the gravitational load of and hold amplifier.
test bearing assembly itself. In this mode, sampling periods should be
matched with those of all high-speed data, to allow
c) Rotational Speed of Test Shaft
exact FFT calculations. In this system, the sam-
pling interval is synchronized with the lower of the
The magnetic sensor @ indicates shaft orientation exciting frequencies.
and provides a pulse for an accurate rotational The data processing program provides the fol-
speed indicator. This pulse signal allows the exciters lowing functions:
to be run in synchrony with the test shaft, and the
resulting sinusoidal forces can also be combined at
the same amplitude to provide a rotating force of
constant magnitude, as required. Low speed data 100 ch
0
I YODAC-8
d) Inlet Oil Temperature and Bearing Metal 99
Temperature
Chromel-alumel (CA) thermocouples are used for
all temperature measurements. The white metal
temperature within 3 mm of the surface is measured
at each of the test bearing locations to obtain the
mean oil film temperature.
0
I
AID
e) Supply Oil Pressure and Oil Flow Rate 15
CONVERTER
244
4.2 Methods of Measurement for Journal Bearing Characteristics
245
4 Explanations
----- Theory
90 L!D=0.5
"0.8
0 .... 75
-~ 0.6 " 60
01
;:. §
·o
•t::
c: 0.4 "
-o 45
§ --Theory -~
30
~ 0.2 L!D=0.5 ~
15
0~~~~~~~~~~~~~--~~~ 0
0.01 0.02 0.05 0.1 0.2
1 0.5 2 5 0.01 0. 02 0.05 0.1 0.2 0.5 2 5
Sommerfeld number S Sommerfeld number S
Fig. 4.2.7. Eccentricity ratio of circular bearings with two oil Fig. 4.2.8. Attitude angle of circular bearings with two oil
grooves grooves
15 15
:.:' 0
~
c 10 0 10
"
,•
0
·u
e 0
"0
0
0 0 0
9 (J 0
•
(.)
0 0 0
9Q
"'"'
0
~
"
0
:.:
Q Q
@ e
";::J 5 Q e e 5
"'
"'
""'
i3 • •
.9
rg
" 0
.§ ----- Theory 0
(:! L!D=0.5
0.01 0.02 0.05 0.1 0.2 0.5 1 2 5 0.01 0.02 0.05 0.1 0.2 0.5 2 5
a Sommerfeld number S b s
15,---- -- -- -- - - - - - - -- -- -- ------, 15
10 10
0
~
~ 5 0
~
5 ooo
0 0
0 0 e
..
0.01 0.02 0.05 0.1 0.2 0.5 5 0. 01 .02 0.05 0.1 0.2 0.5 2 5
c s d s
Fig. 4.2.9. Dimensionless stiffness coefficients of oil film in circular bearings with two oil grooves
246
4.2 Methods of Measurement for Journal Bearing Characteristics
15.-----------------------------. 15.-----------------------------.
---Theory 0
LID=0.5 0
.. 10
Q:> •
.,.o
•
5 ..
.. •
• 0
•
0.01 0.02 0.05 0.1 0.2 0.5 2 5 0.01 0.02 0.05 0.1 0. 2 0.5 2 5
a Sommerfeld number S b s
15.-----------------------------. 15r-----------------------------~
10 10
0
0
..
(.) " J
5 5 .. .. 0
•
. .. • .. 8 •
• • •tJ(J;t{"
..a.
0 J> .. (JC> ..
~_.oo
. ..
• e
J ~~el'~i,~ •
0 0 , .....
0 89
.... e
.. 0 ,.a;."\~o'"~
0.01 0.02 0.05 0.1 0.2 0. 5 2
0
This relationship of modification is obtained by ratios, e > 0.7) and light loadings (i.e., low eccentric-
interpolation of calculated results as discussed in ity ratios, e < 0.3). This certainly cannot be ex-
Sect. 4.5. plained in terms of errors in the nominal clearance
Calculated values are based on constant viscosity assumed or in the measurement of movement in the
assumptions using the finite difference technique. journal or bearing. An allowance for variable vis-
Details of the calculation method are described in cosity certainly seems to improve the agreement
Sect. 4.1. between theory and experiment, but there is still a
Although these results are derived from a combi- significant difference. One possible explanation is
nation of five different experimental data and are nonlinearity of the oil film.
somewhat scattered, it can be said that the experi-
mental data shows fairly good agreement with the
theoretical curves. The experimental data itself b) Theoretical Assumption of Boundary Conditions
should be still examined, mainly in the individual for the Oil Film Extent and Pressures
test condition, as previously mentioned. The pres- at the Film Boundaries
ent data, however, shows the following points of
special interest. Stiffness coefficients show better agreement be-
tween theory and experiment than do damping co-
efficients. It is necessary at this point to examine
a) Nonlinearity of Oil Film boundary conditions for the oil film extent and
pressures at the film boundaries.
The discrepancies between isoviscous theory and Since stiffness coefficients are ratios of the reac-
experimental results become comparatively signifi- tion force change of the oil film to the statical dis-
cant at both high loadings (i.e., high eccentricity placement of the journal centre, only the wedge film
247
4 Explanations
pressure is reflected in these values; damping coeffi- References for Sect. 4.2
cients, in contrast, depend on both wedge film pres-
sure and squeeze film pressure, because the damp- [1] Hagg, A.C., Sankey, G .O., J. Appl. Mech., 23-2 (1 956),
ing force of the oil film arises from both pressures. 302.
The boundary conditions should still be examined [2] Hagg, A.C., Sankey, G.O., J. Appl. Mech., 25-1 (1958),
141.
when the mean bearing pressure and the oil inlet [3] Nakagawa, E., Proc. JSME, 104 (1963), 31.
condition seem to be of particular concern. [4] Kanki, H., Kawakami, T., Esaki, J., Proc. JSME, 810-16
(1981), 57.
[5] Glienicke, J., Proc. IME, 181-3B (1966- 67), 116; Disserta-
tion T.H. Karlsruhe (1966).
c) Reductions in Damping Coefficients [6] Mitchell, J.R., Holmes, R., Ballengooyen, H.V., Proc.
at High Eccentricity Ratios IME, 180-3K (1965-66), 90.
[7] Woodcock, J.S., Holmes, R., Proc. IME, 184-3L (1969- 70),
111.
Among four damping coefficients, Cxx (where x is [8] Morton, P.G., J. Lub. Techno!., Trans. ASME, Ser. F, 93-1
the direction of static load) in particular shows the (1971), 143.
[9] Morton, P.G., IUTAM Symp. Lyngby, Dynamics of Ro-
widest discrepancy between theory and experiment, tors, (1974), 350.
the experimental results being lower. At high eccen- [10] Someya, T., Dissertation University of Tokyo (1972).
tricity ratios the damping coefficients derived ex- [11] Nakanishi, T., Chikira, K., Simoda, H., Hitachi Zosen
perimentally also differ qualitatively from the theo- Technical Review, 37-1 (1976), 49.
[12] Nakanishi, T., Chikira, K., Simoda, H., Hitachi Zosen
retical values. This tendency for a drop in damping Technical Review, 38-4 (1977), 224.
coefficients at high eccentricity ratios is observed [13] Esaki, J., Furukawa, T., Matsumoto, I., Mitsubishi Juko
with almost all test bearings and seems to reflect the Giho, 16-1 (1979), 61.
actual state of the coefficients. [14] Hisa, S., Matsuura, T., Someya, T., International Confer-
The method of measurement for the dynamic ence on Vibrations in Rotating Machinery, Churchill Col-
lege, Cambridge (1980), 223.
characteristics of journal bearings is discussed in [15] Masuda, T., Suganami, T., Oishi, N., Taguma, Y., Proc.
this section. The rigs, transducers and data process- JSME, 813-5 (1981), 91.
ing systems are still subject to improvement. [16] Mitsui, J., Trans. JSME, C, 48-428 (1982), 565.
The dynamic characteristics of oil film depend on in Sect. 1.2.4, H is newly defined as follows, disre-
the temperature of oil film, because oil viscosity is garding H 2 :
heavily dependent on its temperature. If an error is (1)
made in the estimation of oil temperature, the dy-
namic characteristics calculated may be incorrect. where I(is the ratio of heat used for the temperature
Therefore, it is of the utmost importance that the rise of oil film to the total heat generated; the value
effective oil film temperature be estimated correctly 0.6 is commonly used. This practice means that
if the best possible use is to be made of this data- 60% of the heat generated is consumed to raise the
book. Because the temperature is different in the temperature of the oil film and the remainder is
circumferential and axial directions in an actual oil removed from the bearing by heat transfer.
film, however, analysis must take the change of vis- The next equation is obtained for the heat bal-
cosity into account, as shown in Sect. 4.7, for cor- ance of journal bearing using Eqs. (22) and (23) in
rect estimation. In this section a comparatively sim- Sect. 1.2.4, and Eq. (1):
ple method of estimating the effective oil film
temperature is explained, which allows the best pos- (2)
sible use of this databook; it is calculated under the
assumption of constant temperature (viscosity). e
The effective temperature of oil film is finally
given by the following equation, which is derived
from Eq. (2) above and Eq. (30) in Sect. 1.2.4:
248
4.4 Boundary Condition for Oil Film Pressure Distribution
4.5
viscosity is estimated. Next, the Sommerfeld num- a Tilting pad
5 pad
ber S is calculated on the basis of the effective vis- 4.0 LID=0.5
cosity. The friction coefficient, jj,, and the oil flow, L.O.P.
Qs and Qe, are obtained by interpolating data 3.5
shown in this databook. The new effective tempera- mp
•...
0.3: !':.
~
....• 1.5
4.3.2 Estimation from Measured Data "q
1.0
For calculation of the pressure distribution of jour- The following are well known to be the condi-
nal bearings, the following must be known: (1) the tions in which circumferential oil film rupture oc-
pressure at the oil feed position, (2) the pressure at curs:
the bearing end, and (3) the circumferential region
within which the oil film exists. In this databook the a) Sommerfeld condition (no rupture)
pressures at the oil feed position and at the bearing b) Gumbel condition
end are set at zero. This means that the ambient c) Reynolds condition
pressure and the oil supply pressure are equal to the d) Prandtl-Hopkins condition
pressure of atmosphere. The third item (3) is dis- e) Coyne-Elrod condition
cussed in this section. f) Floberg condition
249
4 Explanations
100.,....-----------------------,
These conditions can be applied to the wedge - Boundary condition A
action, but there are not yet any established theo- J 50 -- Boundary condition B
ries that would apply to the squeeze action. Though
-;::
conditions (e) and (f) are considered to fit approxi- <!)
'(]
mately the actual situation, these require compli- s <!)
The two are different in the way the oil film force is ·~
5
calculated, the static and dynamic force compo.-~ E
nents being separated for one and not the other. Iii Q 0. 5
the latter method, Reynolds condition is commonly
applied to the total pressure. In this section, the
following two conditions are examined fbr the fox:,- Sommerfeld number S
mer method. ·
Fig. 4.4.2. Damping coefficient
A. The pressure distribution by squeeze action is
calculated on the assumption that there is no film both these conditions is shown in Fig. 4.4.2. It can-
rupture independent of that caused by wedge ac- not be judged from the example whether either is
tion. The pressure is then integrated within the re- more realistic than the other, though the values
gion of positive oil film pressure by wedge action, obtained with condition B are smaller than those
which is calculated under the Reynolds condition obtained for condition A within the region where
(Fig. 4.4.1 a). the Sommerfeld number is small.
B. The pressure distribution by squeeze action is Next, the stability of the Jeffcott rotor is dis-
calculated on the assumption that the pressure by cussed with reference to these two values. The sta-
squeeze action also becomes zero at the boundary bility chart is shown in Fig. 4.4.3. Here, A. and v are
of film rupture by wedge action calculated on the the modified Sommerfeld number (S/v) and the on-
assumption that the Reynolds condition obtains set speed (w jfg/C), respectively. Further, S is the
(Fig. 4.4.1 b). Sommerfeld number, w is the rotational angular
speed, g is the acceleration of gravity, and C is the
radial bearing clearance. Though the onset speeds
with A and B are equal with large A., they are very
different from each other with small A., and it should
Squeeze action be noted that the onset speed in A becomes ex-
tremely low. Small A. means, physically, that the
...... bearing mean pressure and the bearing radial clear-
ance are large. It is well known that self-excited
a
vibration does not occur in such a case. Therefore,
- Stability limit
---- ( C,-y+ C yx) !( C xx+ Cn )
I
Wedge action 0.8 8
~ Squeeze action
I
I
l
J -----....... '
+~ 0.6
b
C,..) '\ unstable
'\
....... ~~
Fig. 4.4.1 a, b. Pressure boundary condition
B.C. (B) ''\\
\
calculation of the damping coefficients, because 0 . 001 0.01 0.05 0.1 0.5 1 5 10
). ~ Siv
these depend on the squeeze action. An example of
the calculation of damping coefficients for a cylin- Fig. 4.4.3. Sta bility chart. B.C.(A), boundary condition A ;
drical journal bearing with two axial grooves with B.C.(B), boundary condition B
250
4.5 Influences of Length : Diameter Ratio on Journal Bearing Characteristics
it may be concluded that condition B is more real- Here, as the term (kxy- kyx) is the same for con-
istic than condition A. ditions A and B, the difference in stability is caused
Next, the subject is discussed from the viewpoint by the first, second and third terms on the right-
of energy. The infinitely small vibration around the hand side. Now, the assumption ax= aY is made for
equilibrium point of journal is expressed as follows : the sake of simplicity. The energy dissipation, that
is, the stability of the rotor, can then be discussed on
Ax= ax cos(wt + <f>x), the basis of the ratio (cxy + cyx)l(cxx + Cyy), shown
L1y = aY sin (w t + </>y} . in Fig. 4.4.3. With small A values the ratio is smaller
The x andy components of oil film reaction force with condition B than with A. This is consistent
are given by the equations below : with the idea that the energy dissipation is larger
with B than with A and the stability with B is
L1Fx = + kxyL1y + CxxL1x + CxyL1y,
kx xL1x superior to that with A.
L1FY = kyxL1x + kyyL1y + CyxL1x + cYYL1y . From the discussions above it may be concluded
that condition B is more realistic than condition A.
Consequently, the energy dissipated during one Therefore, the data in this book has been calculated
cycle of vibration is obtained as follows: on the basis of condition B whenever the static and
2x dynamic oil film forces are computed separately.
E= J (L1Fx· x + L1Fy· y) dt
0
=nw 2 axcxx+nw
2 2 2
aycyy
- nw axay(cxy + Cyx) sin(<f>x- </>y)
2
102 102
10 1
-L!D=0.5
1.0 -·-LID=0.1
0 L.......................1o.__,...
, ..............................1.._o--~....._...........~~ 10- ' 1o-l 10-' 10-l 1
a s<ff b s<ff c s<ff
Fig. 4.5.1 a-c. Effective Sommerfeld no. for L/D ratio. 5 - tilting pad bearing, mP = 0.5, P= 60°, rx/ P= 0.5, LOP
251
4 Explanations
mally short bearing theory (3]. These simple meth- The value of a decreases with decreasing angle of
ods imply that bearing characteristics for different arc p, and is nearly equal to 1.0 for multilobe bear-
L : D ratios can be approximat ed by introducing an mgs.
effective Sommerfeld number defined by The revised data obtained from Eq. (1) is shown
in Fig. 4.5.1 for calculation no. LBD-3, where
(L/D) }a (L/D)0 is 0.5. It is found that for different L : D ratios
Seff = { (L /D)o S. (1)
the eccentricity ratio e, stiffness coefficient K;j, or
Here, the denominato r (L/ D)0 is the L : D value damping coefficient C;j can be represented by a
that has been selected as standard. single curve. Though the correction factor depends
The value for exponent a is obtained as follows on the bearing type, it is useful in practice to express
from calculated data in this databook: data as a function of the effective Sommerfeld num-
ber for the estimation of data for different L: D
a) for partial arc bearings: values.
log a= 0.75 logp- 1.488 (2)
where p = angle of partial arc (degrees), provid-
ed that References for Sect. 4.5
p= 180° if p ~ 180°
[1] Ocvirk, F.W., NACA TN 2808, (1952).
b) for two-lobe bearings and five-pad tilting pad [2] Warner, P.C., J. Basic Eng., Trans. ASME, Ser. D, 85-2
bearings: (1963), 247.
[3] Nakagawa, E., and Aoki, H., Trans. JSME, 31-229 (1965),
a~ 1.0 (3) 1398.
252
Calculation No. LBD-1
L!D=0.35
Static characteristics Stiffness coefficients Damping coefficients
s rp /Jtp-1
€ Q. Q. Kxx Kxy Kyx K.vy Cxx Cxy Cyx Cyy
11.5 . 0750 83.6 .0700 1. 00 215. 1.52 13.6 -13.1 2.68 27.1 2.23 2.23 17.3
8.55 .100 81.2 .0929 1. 00 160. 1. 59 10.3 -9.68 2.69 20.6 2.25 2.25 13.2
3.99 .200 72.4 .180 1. 00 76.0 1. 81 5.69 -4.42 2.56 11.2 2.40 2.40 7.42
2.38 .300 65.0 .260 .995 46 . 5 2.03 4.36 -2.57 2.37 8,43 2.38 2.38 5.09
1. 51 .400 58.2 .333 .983 30.8 2.40 3.86 -1.57 2.21 7.23 2.19 2.19 3.52
.961 ,500 51.5 .400 • 965 20.8 2.96 3.73 -.894 2.08 6.86 2.34 2.24 2.82
.587 ;600 45.0 .462 . 937 13.9 3.83 3.82 -.375 1. 98 6.80 2.10 2.10 2.09
• 327 .700 38.1 .518 . 894 8.82 5.23 4.10 .0633 1.86 7.19 2.02 2.02 1.59
.153 .800 30.7 .569 .834 5.03 7.85 4,67 .513 1. 76 8,22 1. 99 1. 99 1.18
.0466 .900 21.9 .616 .745 2.22 14.7 5.90 1. 07 1. 63 0.4 1. 80 1.80 .724
,0160 .950 16.1 .643 • 701 1.12 30.1 8.66 2.60 l. 84 15.9 2.43 2.43 .693
253
Calculation No. LBD-1
L!D=O.S
Static characteristics Stiffness coefficients Damping coefficients
s /J¢·1
e ¢ Qs Qe Kxx Kxy Kyx Kyy Cxx Cxy Cyx Cyy
5.96 .0750 82.6 .0663 1. 00 111.5 1. 77 13.6 -13.1 2. 72 27.2 2.06 2.06 14.9
4.43 .100 80.4 . 0880 .999 83.0 l. 75 10.3 -9.66 2.70 20.6 2.08 2.08 11.4
2.07 .200 71.6 .170 . 999 39.5 1. 88 5. 63 -4.41 2.56 11.2 2.22 2.22 6.43
1. 24 . 300 64.1 .243 . 997 24.4 2.07 4 .2 7 -2 . 56 2.34 8.50 2.32 2.32 4.73
. 7 98 . 4 00 57.5 .309 .986 16.4 2.39 3.75 -1.57 2.17 7.32 2.24 2.24 3.50
.517 • 500 51.1 .366 . 967 11.3 2.89 3.57 -.924 2.03 6.81 2.10 2.10 2.60
. 323 .600 44.7 .416 . 939 7.75 3.65 3.62 -.427 1.92 6.81 2.08 2.08 2.06
.187 . 700 38.2 .459 • 900 5.13 4. 92 3.88 .023~ 1. 83 7.32 2.16 2.16 1. 70
.135 . 750 34.9 .478 .875 4.07 5.90 4.11 .258 1.80 7.65 2.10 2.10 1. 46
.0582 .850 27.2 .510 .806 2.30 9.56 4.92 .805 1. 73 9.12 2.06 2.06 1. 06
.0315 .900 22.7 .524 .758 l . 56 13. 8 5.76 1. 24 l. 72 10.6 2.03 2.03 .846
.00499 .975 12.5 .543 .663 .54 43.3 9.61 4.26 1. 99 22.4 2. 93 2. 93 .637
L! D=0.75
Static characteristics Stiffness coefficients Damping coefficients
s ht/J-1
e ¢ Qs Q. Kxx Kxy K yx Kyy Cxx Cxy Cyx Cyy
2.99 .0750 82.3 . 0596 .999 56.0 1. 83 13.5 -13.1 2.71 27.2 1. 82 1. 82 11.8
2.22 .100 80.1 .0791 .999 41.7 1. 79 10.3 -9.66 2.67 20.6 1.84 1.84 9.08
1. 05 .200 70. 8 .152 . 999 20.0 1. 93 5.54 -4.39 2.55 11.2 1. 96 1. 96 5.19
.637 .300 62.2 .210 . 998 12.6 2.17 4.10 -2.54 2.31 8.64 2.16 2.16 4.04
.418 .400 55.9 .261 .988 8.66 2.39 3.52 -1.60 2.09 7.44 2.11 2.11 3.11
.278 .500 50.3 .308 • 972 6.16 2.76 3. 32 -.960 1. 97 6.91 2.11 2.11 2.53
.110 .700 38.3 .367 .905 3.12 4.46 3.52 -.0564 1. 77 7.14 2.03 2,03 1. 61
.0590 .BOO 31.8 .384 .849 2,08 6.54 4, OS .437 1. 74 8.01 2.04 2.03 1.26
.00963 • 950 17.6 .396 .707 .741 24.9 7,88 2.65 1. 94 14,5 2.41 2.41 .760
254
Calculation No. LBD-1
LID=!
Static characteristics Stiffness coefficients Damping coefficients
s Qs Q, JJ,p-1
€ ~ K xx Kxy K yx Kyy Cxx Cxy Cyx C yy
2.95 .0500 83.0 • 0351 .999 55.1 2.38 20.1 -19.8 2.90 ~0.7 1.65 1.65 1~.5
1. 96 . 075 80.8 . 0523 .999 36.6 2 . 18 13.5 -13.0 2.88 27.4 1. 66 J.. 66 9.81
.284 . 400 54.6 .220 . 991 5.93 2.37 3.34 -1.60 2.06 7.57 2.10 2.09 2.94
.194 .500 49.1 .251 .972 4.37 2 .6 7 3.09 -1.0 0 1. 89 6,99 2. 00 2.00 2.31
.00866 . 950 18.0 .297 .710 .684 24.0 7.81 2. 7 9 2.01 14.2 2 .4 4 2.44 . 7 91
L!D=1.5
Static characteristics Stiffness coefficients Damping coefficients
s Q, Ijt/J-1
€ ~ Qs Kxx Kxy Kyx Kyy C xx C xy Cyx C yy
1. 82 .0500 82.7 .0268 .999 34 .1 2.~ 7 20.0 -19. 7 2.95 40.9 1. 47 1. 47 11.3
.902 .100 78 .4 .0528 .997 16. 9 2.09 10. 1 -9. 62 2 .7 8 20. 8 1. 4 9 1. 4 9 5.9 2
.434 .200 68.3 .0985 .996 8.34 2 . 10 5.28 -4.34 2.59 11.4 l . 59 l . 59 3.53
.376 .800 32.4 .198 .859 1. 40 5.61 3 . 55 .375 1.71 7.81 2 .03 2 .03 1. 28
.00786 .950 18.4 .188 .712 .637 22.9 7.64 2.81 2.05 14.1 2.51 2.51 .83
255
Calculation No. LBD-1
L!D=0.35
--- - - L!D=l.S
---- L! D=0.35
2
s
256
Calculation No. LBD-1
- - - - - - - LID=0.35
------- LID=0.75
- - - - - - LID=1.5
20
,,
2
~
/
/
I
I
I
I
I
s
Ioor--------------------------------------------------------~
- - - - - - - L/ D= 0. 35
------- LID=0.75
----- -------
257
Calculation No. LBD-2
L/D=0.5
Static characteristics Stiffness coefficients Damping coefficients
s Q, /Jt/J-1
€ ~ Qs Kxx Kxy Kyx Kyy Cxx Cxy Cyx Cyy
3.52 .050 78.8 .351 139. 145. 28 . 9 -75.1 14.6 196. -68.0 -67.7 43 .1
1. 75 .100 78.7 .351 69 .3 72.9 14.4 -37.5 7. 31 97 . 6 -33.7 -33. 6 21.6
.858 .200 78.5 .353 36.7 37.2 7. 71 -18.9 3.92 48.6 -16.3 -16.2 10.7
.554 .300 78.2 .357 22.3 25.1 5.54 -12.3 2.75 32.5 -10.3 -10.2 7.12
.397 .400 77.9 .362 16.2 19.1 4.53 -8.86 2.19 24.5 -7.13 -7.11 5.30
.299 .500 77.4 .368 12.5 15.5 4.02 -6.71 1.89 19.7 -5.04 -5.02 4.19
.230 .600 76.7 .375 9.89 13.3 3.60 -5.11 1.65 16.1 -3.63 -3.61 3.39
.178 .700 75.7 .383 7.97 11.8 3.54 -3.87 1.53 13.9 -2.39 -2.37 2.81
.136 .800 74.5 .392 6.28 10.8 3.51 -2.80 1.43 12.1 -1.43 -1.42 2.31
.101 .900 72.7 .402 5.39 10.3 3.74 -1.82 1.40 11.2 -.441 -.432 1.91
.0841 .950 71.4 .407 4.99 10.3 3. 77 -1.28 1.35 10.6 -.045 -.036 1.68
258
Calculation No. LBD-2
L!D=l
Static characteristics Stiffness coefficients Damping coefficients
s
€ ¢ Qs Q, lj<P-' Kxx Kxy Kyx Kyy Cxx Cxy Cyx Cyy
1. 73 .050 81.4 . 196 66.4 132 . 33.7 -65.4 11.6 210. -69.4 -69.1 41.5
.863 .100 81.3 .197 33.1 65.9 16.8 -35.6 5.!)9 104. -34.4 -34.3 20.7
.425 .200 81.1 .199 16.4 33.7 9.07 -16 . 5 3.08 52.1 -16.7 -16.6 10.4
.276 .300 80.7 .201 10.7 22.7 6.33 -10.8 2.16 34.5 -10.6 -10.5 6.92
.200 .400 80.3 .205 7.87 17.1 5.03 -7.87 1. 74 25.8 -7.37 -7.33 5.19
.152 .500 79.7 .209 6.12 13.8 4.28 -6.01 1. 54 20.6 -5.24 -5.21 4.15
.119 .600 78.9 .214 4.92 11.7 3.75 -4.65 1.41 16.8 -3.77 -3.75 3.41
.0942 .700 78.0 .219 4.05 10.2 3.49 -3.59 1.35 13.9 -2.73 -2.71 2.83
.0742 .BOO 76.7 .225 3.40 9.25 3.34 -2.69 1.30 11.9 -1.85 -1.84 2.35
.0575 .900 74.9 .231 2.98 8.73 3.27 -1.83 1.25 10.3 -1.08 -1.07 1.92
.0497 .950 73.8 .234 2.92 8.42 3.26 -1.37 1.23 9. 77 -. 717 -.710 1.72
LID=1.5
Static characteristics Stiffness coefficients Damping coefficients
s ¢ Qs Q, fJcp-1 Kyy Cyx Cyy
€ Kxx Kxy Kyx Cxx Cxy
1.38 .050 82.5 .126 51.9 124. 36.0 -61.0 10.9 218. -69.7 -69.4 40.6
.682 .100 82.4 .126 25.7 62.0 17.2 -30.2 5.33 106. -34.3 -34.2 20.2
.336 .200 82.2 .127 12.8 31.5 9.49 -15.3 2.78 52.9 -16.7 -16.6 10.1
.219 .300 81.8 .129 8.39 21.2 6.69 -10.1 1.96 35.2 -10.6 -10.5 6.78
.159 .400 81.3 .131 6.16 16.0 5.28 -7.39 1.60 26.3 -7.31 -7.30 5.11
.122 .500 80.7 .134 4.81 13.0 4.38 -5.66 1.39 20.7 -5.33 -5.31 4.08
.0959 .600 79.9 .138 3.89 10.9 3.88 -4.42 1.30 16.6 -3.99 -3.98 3.35
.0765 .700 78.9 .141 3.23 9.53 3.47 -3.42 1.25 13.7 -2.93 -2.92 2.80
.0610 .800 77.6 .145 2.74 8.55 3.29 -2.61 1.23 11.5 -2.09 -2 .08 2.33
.0479 .900 75.8 .149 2.45 7.97 3.22 -1.84 1.22 9.99 -1.31 -1.30 1.92
.0419 .950 74.6 .151 2.45 7.69 3.12 -1.40 1.20 9.40 -.924 -.918 1.72
259
Calculation No. LBD-2
L!D=0.5
------ LI D= I
- - - L!D=l.5
1.0~------------------------------------------------------~
L!D= O.S
LID= I
L!D=l.5
260
Calculation No. LBD-2
Ll D=0 .5
------ LID=!
- - - LID=l.5
- - - L ! D= 0.5
------LID=!
---LID=1.5
...,
G
261
Calculation No. LBD-3
L/ D=0.3
Static characteristics Stiffness coefficients Damping coefficients
s
e QIN ljl/1-' Kxx K yy Cxx Cyy
262
Calculation No. LBD-3
L!D=O.S
- 7
LID-o
Static characteristics Stiffness coefficients Damping coefficients
s
€ QIN ht/J-1 Kxx Kyy Cxx Cyy
263
Calculation No. LBD-3
5.0
---------
:; 4.5
01
0.6
---- --- ---
0.4--
4.0 ---L/D=0.7
or
s
- - - L!D=0.3
------ L/D=0 . 5
- - - L!D=0.7
264
Calculation No. LBD-3
---L!D=0.3
LID=0.5
20 - - - L!D=0.7
....
~
50 ---L!D=0.3
L / D=0.5
20 - - - L/D=0 .7
265
4 Explanations
\
0.8
- · - Rc=aooo
ity. In the superlaminar region the load-carrying
capacity and the friction coefficient increase with 0.6
increasing Reynolds number and the discrepancies
between the actual bearing characteristics and the 0.4
data calculated by the laminar theory widen. In this
section the influences of superlaminar flow are dis- ~.
cussed with reference to the data listed below 0.2
~
(pp. 267- 277). ~
0
0. 1
a s,,,
Calculation no. Bearing type Specifications
TBL-1 Two-lobe L:D = 1, mP = a.67, p = 1 7ao,
bearing 2 oil grooves of 10° each
TBL-2 Four-pad L: D = 1, mP = a, p= sao, 10 -Re=O
tilting pad rxfp = a.5, LBP ----- Re=3000
bearing
266
Calculation No. TBL-1
Re= 0
Static characteristics Stiffness coefficients Damping coefficients
s
€ ¢ Qs Qe h"'-1 Kxx Kxy K yx Kyy Cxx Cxy Cyx Cyy
4.90 .0172 81.2 .218 .379 133 391 117 -198 30.7 609 -197 -197 120
2.92 . 0289 81.2 .218 .379 78.8 234 67.2 -118 18.2 361 -118 -117 70.5
1.97 .0430 81.2 .218 .379 53.2 158 44.6 -79.6 12.3 241 -79.2 -79.2 47.3
• 956 .0898 81.2 .219 .378 25.5 77.3 20.9 -39.0 6.07 113 -39.4 -39.2 23.8
.521 .165 81.0 . 219 .377 14.1 42.3 12.5 -21 . 1 3.42 63.5 -20.7 -20.7 12.9
• 305 . 277 80 . 7 . 222 .372 8.26 25.6 7.40 -12.3 2.23 38.1 -12.2 -11.8 7. 80
.192 .427 80.0 .228 .363 5.33 16.8 5.30 -7.80 1. 71 24.7 -7 . 00 -6.98 5.15
.100 .655 78.3 .239 .344 2.93 10.7 3.82 -3.97 1. 34 15.0 -2.81 -2.79 3.01
.0513 .925 74 . 0 .256 .311 1. 66 8.92 3.55 -1.61 1. 24 11.1 -.49C -.481 1.80
.0313 1 . 07 69.1 .263 . 292 1.15 9. 61 4.02 -.401 1. 36 10.0 • 33( .345 1.29
.0202 1.15 64.0 .265 .281 • 865 11.7 4.80 . 570 1.56 10.3 • 87 .873 • 995
.0104 1. 15 53.9 .247 . 296 .559 20.2 6.60 2.24 1. 90 11.5 1.55 1.58 • 731
267
Calculation No. TBL-1
Re =3000
Static characteristics Stiffness coefficients Damping coefficients
s
e ¢ Qs Q, ljtV' Kxx Kxy Kyx Kyy Cxx Cxy Cyx Cyy
4. 96 .0119 85.2 .250 .393 366 567 209 -299 24.9 908 -318 -307 203
3.11 .0189 85.2 .250 .393 230 356 131 -188 15.6 570 -199 -193 127
2.09 .0283 85.2 .250 • 393 153 237 86.0 -126 10.4 377 -134 -129 84.9
• 974 .0629 85.2 .250 .393 71.5 111 39.8 -58.6 5 . 05 176 -62.3 -61.2 39.8
.517 .119 85.0 . 251 .392 37.8 59 . 5 21 . 3 -31 . 3 2.80 93.6 -32.8 -31.9 21.1
. 315 .195 84.6 .253 . 387 23.2 36.6 13.7 -19.3 1. 71 58.3 -20.0 -19.2 13.1
. 196 .310 84.3 .258 .380 14 . 4 23.4 8.52 -12.0 1.32 36.3 -12.1 -11.7 8.42
.0992 .535 83.0 .274 .353 7.37 13.0 5.21 -6.19 1.11 20.0 -5.64 -5.44 4.77
.0514 .812 81.1 .303 . 299 3.85 8.58 3.62 -3.01 1.02 11.9 -2.24 -2. 14 2.92
.0306 1. 07 77.5 . 342 . 212 2.43 7. 55 3.91 -1.48 1.04 9.37 -.422 -.340 2.12
.0202 1. 18 73.3 . 353 .169 1.68 8.05 4.48 -.431 1. 38 8.99 .678 .743 1.72
.0105 1. 26 63.8 . 352 .139 .986 12.4 6.13 1. 41 1. 88 10.2 1. 78 1.86 1.28
Re=5000
4.88 .0095 85.9 .258 . 394 467 682 262 -365 26.5 1110 -394 -377 255
3.02 .0154 85.9 .258 .394 289 422 163 -226 16.3 685 -244 -233 157
2.07 .0225 85.9 .258 . 394 198 290 110 -155 11.2 466 -167 -160 107
1.02 .0459 85.9 .259 .393 96.9 142 53.0 -76.4 5.52 227 -82.4 79.1 52.9
.514 .0927 85.9 . 259 .393 48.5 72.6 26.1 -38.5 2.84 110 -41.4 -39.7 26.8
.301 .158 85 . 8 .260 .390 28.6 42.5 16.2 -22.6 1. 75 68.0 -24.0 -23.2 15.9
.202 .253 85.4 .263 .388 19.2 29.0 11.0 -15.2 1.25 45.8 -15.9 -15.5 10.9
.0997 .458 84.6 .275 .380 9.51 15.2 5. 81 -7 . 60 . 943 22.8 -7.45 -7.12 5.74
.0486 .750 82.8 • 302 • 315 4.69 9.11 3.95 -3.65 .874 12.9 -3.06 -2.71 3.17
.0299 1.01 79.1 • 338 .251 3.00 7.24 3.62 -2.01 • 980 9. 92 -1.09 -.955 2.36
.0199 1.16 75.4 .354 .203 2.08 7.46 4.00 -.904 1.16 8.72 .0796 .165 1.85
.0103 1. 30 67.4 .372 .129 1.14 12.9 5.56 1. 74 1.64 8.42 1.28 1.40 1.27
268
Calculation No. TBL-1
1.4
€
- - · - - Re= 0
1.2 "",... -------- Re=3000
,•,
~ Re =5000
1.0 1.0
''
\\ '
----- ---- ----- ----- -
\
0.8 0.8 ---------------------~------
0.6 0 .6
c::y w
Qs
45
0.4 0.4
0.2 0.2 - 30
15
0 0
0. 005
s
50
20
10
~
'
~
2
Re= 0
1 Re =3000
R e =5000
0.5
0.2~~~--~-~~~~~~-~-~~~~~---L-~~~
0.005 0.01 0.02 0.05 ~~
0.1 0.2 0.5 2 5 10
s
269
Calculation N o. TBL-1
Re= 0
Re =3000
Re=5000
5 10
s
-- Re= 0
I
I
Re=5000
I
0. 05 0.1 0.2 0.5 1 2 10
s
270
Calculation No. TBL-1
1. 0
9
8
z5
0.1
9
8
l .....;;
.......... !'-...
""'
5 ........
4
I'-
3
2
I'..
I'-..
I' Gx
0.0 1
9
8
"""'
7
6
Gz
4
3
0.00 1
3 4 5 6 7 8910' 3 4 8
5 6 7 910' 345678910'
10 2
271
Calculation No. TBL-2
4. Special items
mP: 0
P: 80°
rxfP: 0.5 y
Calculation method: FDM
Boundary conditions: Reynolds condition
Number of grid points: Axial direction: 18
(per pad) uniform spacing)
Circumferential di-
rection: 18 (uniform
spacing) X
Re= 0
Static characteristics Stiffness coefficients Damping coefficients
s
€ rp Q. Qe JJ.p-l Kxx Kxy Kyx Kyy Cxx Cxy Cyx Cyy
272
Calculation No. TBL-2
Re=lOOO
Static characteristics Stiffness coefficients Damping coefficients
s
€ ¢ Qs Q, JJ.p-1 Kxx Kxy Kyx Kyy Cxx Cxy Cyx Cyy
Re =3000
Static characteristics Stiffness coefficients Damping coefficients
s
€ ¢ Q. Qe JJ.p-1 Kxx Kxy Kyx Kyy Cxx Cxy Cyx Cyy
273
Calculation No. TBL-2
Re =5000
Static characteristics Stiffness coefficients Damping coefficients
s
€ ¢ Qs Qe JJ,p-l Kxx Kxy Kyx Kyy Cxx Cu Cy x Cyy
R e=BOOO
Static characteristics Stiffness coefficients Damping coefficients
s
€ ¢ Qs Qe JJ,p-• Kxx Kxy Kyx Kyy Cxx Cxy Cyx Cyy
274
Calculation No. TBL-2
Re= 12000
Static characteristics Stiffness coefficients Damping coefficients
s
€ ¢ Q., Q,. !J<P- ' K.rx Kx.v K.vx Kyy Cxx Cxy Cyx Cyy
275
Calculation No. TBL-2
1.8 0.6 60 .
~ (<) -e..
0.2 0.4 40 .
0 0 0
10
200~----------------------------------------------------------~
10
s
276
Calculation No. TBL-2
50r---------------------------------------------------------------.
10
50~----------------------------------------------------------~
....,
G
10
s
277
4 Explanations
In fluid film bearings, heat is generated in the film When the dimensionless film temperature
due to journal rotation and is dissipated by heat T=lyKifi 2 (T- T.n)f(pincv) is introduced, the di-
transfer. The hydrodynamic lubrication theory mensionless energy equation can be expressed as
dealing with the effect of heat generation and dissi- follows:
pation is called the thermohydrodynamic theory _, ar _,ar _, ar _ o2r
ao + v a;; + w a[ = P., a,2 + <P
1
(THD theory). This theory can give the precise per- u (3)
formances of fluid film bearings. The isoviscous the-
ory (IV theory) in contrast, cannot give the bearing where
temperatures, by definition, because it is based on
the assumption of constant temperature and there-
fore also of constant viscosity throughout the film.
The effective film temperature has been usually as-
sumed in the IV theory to obtain the bearing char-
v' = lf/- 1 H(v- ii'llfl ~;) (4)
!G:
acteristics. Although this method is very simple, the
prediction of bearing performance is inaccurate, be-
cause no generally accepted method of estimating v =- If/ H + ~~ - ; ~: ~~) d'l
the effective temperature has yet been proposed.
Recent investigation [1, 2] has shown that an effec-
tive viscosity applicable to the calculation for all
bearing characteristics does not exist and that ap-
plication of the THD theory is necessary for precise and j1 = fll flin (dimensionless viscosity coefficient),
prediction of the bearing characteristics. An outline P., = yKC 2 cvj2 (Peclet number), and !lin is the vis-
of the THD theory is described in this section. cosity coefficient of the supplied oil.
The steady state Reynolds equation for THD For the THD case the Reynolds equation, the
theory can be obtained by settings' = 0 and¢' = 0 energy equation of the lubricating film and the heat
in Eq. (10) of Sect. 4.1.1: conduction equations in the bearing and the jour-
nal should be solved simultaneously. In the case of
aea (H 3
E
aA)
ao + a(a (H 3
E
aA)
a( = aea (H G) (1) unsteady state condition, Reynolds [Eq. (9) in
Sect. 4.1.1] is simplified by using the two-dimen-
sional viscosity distribution flm(O, () which can be
where E and G are parameters expressing the effect calculated by averaging jl(O, 17, ()in the film thick-
of the film viscosity distribution and can be calcu- ness direction and given as follows [3]
lated from Eqs. (6) and (7) in Sect. 4.1.2. When E
1
and G are given, the film temperature is calculated
by solving the energy equation of the film. The
flm(O, () = Jfl(O, 17, () d17. (5)
0
energy equation for the fluid film bearings is given
by the following equation: E and G are expressed as follows:
ay + (ow)
ooy2T + fl {(ou)
2 2
= , 2
JA ay } (2) Substituting these values forE and G into Eq. (9) in
Sect. 4.1.1, the following Reynolds equation for the
where unsteady state case can be obtained:
(H (H
y: specific gravity of the lubricant
K: specific heat of the lubricant a 3
aA) a 3
aA)
J: heat equivalent of work ae flm . ae + o( flm . o(
2: heat conductivity of the lubricant
T: film temperature
ft: viscosity coefficient of the lubricant
= 6 ~; + 12 (s</J' sine- 8 1
cos 0). (7)
278
4.7 Analysis of Bearing Characteristics with Due Consideration for Viscosity Variation
4. 7.2 Examples of Temperature ing obtained by THD theory. In all figures, the
Distribution in the Bearing Bush Sommerfeld number S* = fl;nNLD /(W'f/ 2 ) is plot-
ted along the abscissa. In Fig. 4.7.2, fl;n is the viscos-
ity coefficient of the supplied oil. The eccentricity
Figure 4.7.1 shows the typical temperature distribu- ratio and the attitude angle are shown against S*
tion at the central plane of the bearing bush. The with Peclet number Pe as a parameter. The friction
experimental isotherm (i) is found to be in good coefficient fJ '1'- 1 and the dimensionless film stiff-
agreement with the theoretical one (ii) by the THD ness and damping coefficients K, C are shown in
theory [8). On the metal surface, the point of maxi- Figs. 4.7.3, 4.7.4 and 4.7.5, respectively. The Peclet
mum temperature is located near the position of the number is one of the most important parameters
minimum film thickness and the minimum temper- which affect the temperature gradient in the film
ature, near the oil groove at the crown. The heat thickness direction.
which flows to the normal direction of the Figure 4.7.6 shows the experimental and theoret-
isotherms is conducted in the bearing bush and also ical results [4] for dimensionless load-carrying ca-
dissipated to the surroundings. The temperature on pacity - the eccentricity ratio curve for the circular
the outer surface of the bearing bush varies consid- journal bearing with one oil groove at the crown of
erably in the circumferential direction. the upper metal part. The experimental condition is
as follows:
a Experimental b Theoretical
Fig. 4.7.1 a, b. Isotherms in the bearing bush. Eccentricity ratio = 0.7, Attidue angle = 43o . _ .
N = 2502 rpm; W = 5.61 N; P;n = 98 kPa; T;n = 40.0 °C; I;, = 29.2 oc; Clearance ratio = 0.00157; Lubncant = 90 turbme otl
0.8 75'
279
4 Explanations
100r-------------------,
50
20
10
Pe =7 .05
0.5
100
50
Kxx
,, , .
, ;
,, , ,
20 ;
, ,, ,
10 Kxy ,,, ,,
,
:::;-
,,
~
,
, ;",
,
5 ,' ;
- -·
2
Kyy
0.5
----- Isoviscous theory
Fig. 4.7.4. Dimensionless stiffness coefficients of oil film
0.2 according to the THD theory and the isoviscous theory
0.005 0.01 0.02 0.05 0.5 1 2 5 10
S*=p.;nNLD/(W¢ 2 )
100r-------------------~
50
20
10
·~
1:...) 5
- - THO theory
Pe=3.52
0.5
----- Isoviscous theory
0. 2~~-:-:-....,...::-::-'---';:-~~'-:--::'-..,..-'......L::i+.I...LL~---,!,--"0......0..~-LUJ
0.005 0.01 0.02 0.05 0.1 0.2 0.5 2 5 10
Fig. 4.7.5. Dimensionless damping coefficients of oil
S*=p;nNLD/(W¢' 2 ) film according to the THD theory and the isoviscous
theory
280
4.7 Analysis of Bearing Characteristics with Due Consideration for Viscosity Variation
I.Or---------------------------------------,
• 500rpm
+ 1750rpm
THD theory ~.,
(2\~
0.8
0 • 2750rpm
·~
.A. 3500rpm
.... 0 6
a ·
:g
1l 0.4 THD theory
4~~~
./_,.""'
8 (17 SOrpm) ,.:-.
,. ""- Isoviscous theory
w ~z
0.2
0~~~~~~~~--~~~--~----~~--~
Fig. 4.7.6. 0.01 0.5 2 5 100.1 0.2
Experimental and theoretical eccentricity ratios
P~ = W¢' /(p.;n (J)LD)
2
X I0-3
The experimental results are obtained by varying
25 the bearing load under the constant revolutional
0 1000rpm speed. The THD theory agrees well with the exper-
<> 1500rpm iment. With the increase of the speed, the eccentric-
20
•... 1750rpm
2500rpm 0
ity ratio against the same Po* is found to increase
because of the decrease in the film effective viscosi-
0 2750rpm ty.
6 3000rpm Figure 4.7.7 shows the experimental and theoret-
ical frictional coefficients. The frictional coefficient
..:::: decreases with the increase in speed because of the
..... 15
c:
OJ
·u same effect as described above.
e
OJ Little data on the stiffness and damping coeffi-
0
u
cients of the film for the THD case has been provid-
:f<
c:
0
ed so far [5, 6]. Recent research work showed that
:~ 10
~ the theoretical film coefficients deriving from the
THD theory were in good agreement with the ex-
perimental ones [7].
, TI!Dthomy
(2500rpm)
5
,!3
~~ References for Sect. 4. 7
0~~----~----~------~----~ [1] Seireg, A., Ezzat, H., Trans. ASME, Ser. F, Vol. 95, No. 2
0.05 0.1 0.2 0.5 1.0 (1973), 187.
[2] Yu, T.S., Szeri, A.Z., Trans. ASME, Ser. F, Vol. 95, No. 3
s• =p.;nNLD/(W<p2) (1973), 298.
(3] Mitsui, J., Yamada, T., Bull. of the JSME, Vol. 22, No. 172
Fig. 4.7.7. Experimental and theoretical friction coefficients
(1979), 1491.
[4] Mitsui, J., Trans. of JSME, Ser. C, Vol. 49, No. 438 (1983),
199.
(5] Craighead, I.A., Dowson, D., et al., I. Mech. Engrs., Paper
II (i) (1980), 47.
[6] Boncompain, R., Frene, J., I. Mech. Engrs., Paper I (iii), 33.
[7] Mitsui, J., Trans. of JSME, Ser. C, Vol. 48, No. 428 (1982),
565.
281
4 Explanations
4.8.1 Introduction becomes the centre of gravity, which means that the
phase lag becomes about 180°. This phenomenon
The vibration of a rotor supported by journal bear- becomes clearer at (6). At the previous rotating
spee~ it is mainly unbalance vibration that appears,
ings is different from that of one supported by ball
bearings: the rotor becomes unstable in the high- but m modes (7) and (8) the unstable vibration
speed range and its vibration is nonsynchronous mode whose frequency is the same as the natural
and different from the unbalance vibration. This frequency is beginning to be superposed on the un-
section is devoted mainly to synchronous vibration , balance vibration mode. In modes (9) and (10) the
. unstable vibration mode becomes more prominent
I.e. unbalance vibration.
Figure 4.8.1 shows the difference between unbal- and the unbalance vibration mode is eliminated.
ance vibration and unstable vibration. Fig- This unstable vibration is called "oil whip".
ure 4.8.1 a shows the phase angle, which is defined In this section the method of obtaining the unbal-
as the angle between the unbalance location and the ance response is shown for a simple model, and
maximum amplitude. Figure 4.8.1 b shows the un- resonance amplitude, vibration mode and whirl or-
t~:~:~:
balance response curve, unstable vibration curve (1)
and measured vibration curve. Figure 4.8.1 c shows
the frequency for the two vibration modes. Fig-
r"7t"7t"71l
ure 4.8.2 shows vibrational modes for the point
number on the response curve shown in Fig. 4.8.1 b. (2)
A pulse mark on the abscissa means an unbalance
point fixed in the rotor, i.e. the pulse is marked once
per revolution.
n~~~
The vibrational mode (1) shows the case of low (4)
rotating speed and in this case the amplitude is
small and the phase angle¢ is almost zero. In mode
(2) the amplitude becomes large and the deflection
b:::~
of the rotor can not instantly follow the unbalance
r-\C: ~~=:
(6)
force. In mode (3) the rotating speed coincides with
the natural frequency of the rotor system (resonan-
ce frequency), so that the phase angle is 90°. In
mode (4) the phase angle becomes more than 90° (7) (8) A (\ 0 (\ 0 (\
and in mode (5) the rotor deflection cannot follow
the unbalance force and a centre of spin of the rotor
(9)0!~
i'"'"l v:
a
j 9: :_======:?=====~:==============- 10
y-
Fig. 4.8.2. Vibration wave
Unbalance vibration
4
/Measured vibration
8
5 67 •
b
m : Mass of disc
k : Shaft stiffness
l : Span length
o· : Static equilibrium point
0 J : 1ournal centre
WI
rpm 0 • : Center of disc
G : Mass centre
Fig. 4.8.1 a-c. Relations among phase angle, amplitude and fre-
quency Fig. 4.8.3. Model rotor and coordinate system
282
4.8 Unbalance Response of Rotors in Journal Bearings
bit are theoretically represented. From the view- 4.8.3 Unbalance Response
point of rotor-bearing system design the method of
selecting the bearing type and bearing configura-
In this analysis, the force of gravity on the rotor is
tion so as to reduce the amplitude of the critical
neglected. Then unbalance response is obtained as
speed is discussed.
a particular solution of Eq. (3). The unbalance force
in Eq. (3) is represented as
4.8.2 Equation of Motion
f(t) = {0, 0, eRw 2 coswt, 0, eRw 2 sinwt, Of
For the sake of simplicity, the Jeffcot rotor model
= {0, 0, eRw 2, 0, 0, O}r coswt
shown in Fig. 4.8.3 is considered, and it is assumed
that the rotor is located at the centre of the shaft + {0,0,0,0,eRw 2,0}T sinwt
and that bearing characteristics on both sides are =fc coswt + fs sinwt (4)
equal. The equilibrium point of the journal centre 0'
is selected as the origin of the rotor coordinate thus the particular solution of Eqs. (3) and (4) is
system. The coordinate of mass centre of the rotor
is selected as (x, y) and the journal centre is selected
x = ~ coswt + ~ sinwt
as (x 0 , y 0 ). The equation of motion is represented in = {~ 1 , ~ 2 , ... , ~ 6 f coswt
ref. [2] as + {F. 1 , F. 2 , ... , F. 6 } T sin w t . (5)
mx + k(x- x 0 ) = meRw 2 coswt + mg,
Substituting Eq. (5) into Eq. (3) we get
my+ k(y- y 0 ) = meRw 2 sinwt, (1)
wAFS + BPc =fc,
k(x- x 0 ) = J>,;, k(y- Yo)=~. (6)
- wAPc +BPs= fs·
where m is rotor mass, k shaft stiffness, w rotating
speed of rotor, eR eccentricity of the rotor and P,; This equation is rewritten in matrix form as
and ~ the journal bearing forces which are repre-
sented as
[ B wAJ
- wA B
[Pc] {/c}.
Ps
=
fs
+ kxy Yo+ cxxXo + CxyYo + mg,
P,; = kxxXo
P k (2)
y = yxXo + k yyYo + CyxXo + CyyYo,
0 0
w:]-1{~:}.
ing coefficients of the journal bearing, respectively.
These values are listed in the tables and figures in
this book in the dimensionless form as
[~] = [-~A (7)
{ x} =
y
{PC2
cos w t
Pcs coswt
+ IS 2 s~n w
+ Pss smwt
t} . (8)
where CP is a radius clearance of the journal bear-
ing, W is the static load of the bearing and w, the The coordinates of the shaft at the bearing are rep-
rotating speed of the rotor. Equation (1) and (2) can resented as
be rewritten in matrix form as
Ai+Bx=f(t) (3) {x0 } = {Pc
3 coswt + f's 3 s~nwt}.
(9)
Yo Pc 6 coswt + f's 6 smwt
where
x1 x 0 0 Cxx 0 0 Cxy In the numerical calculation the bearing coeffi-
x2 X 0 0 Cyx 0 0 Cyy cients must be calculated for each rotating speed,
x3 Xo 1 0 0 0 0 0 because the bearing coefficients are functions of the
{x} = (A)=
x4 y 0 1 0 0 0 0 rotating speed.
Xs y 0 0 0 1 0 0
x6 Yo 0 0 0 0 1 0
283
4 Explanations
Rotating direction
4.8.5 Numerical Calculation
x,Re x, Re of Unbalance Response
Fig. 4.8.4. Calculation of elliptic whirling amplitude
In this section unbalance response is calculated for
a simple rotor system to illustrate the effect of bear-
The motion of a disc centre in the horizontal and ing types and some bearing parameters. The data of
vertical directions is represented as the model rotor is shown in Table 4.8.1. Figure 4.8.5
x(t) = Pe 2 cosOJt + f's 2 sinOJt, shows the unbalance response and whirling orbit
(10) for a cylindrical bearing. The axis of the abscissa
y(t) =Pes cosOJt + Pss sinOJt. being a rotating speed nondimensionalized by the
If the motion of the centre of the journal is to be natural frequency of a simply supported rotor OJ 8 ,
represented, ~ 2 and P. 2 are replaced with ~ 3 and and the axis of ordinate, an amplitude nondimen-
Jl. 3 , and ~sand P.s with ~ 6 and }J. 6 • Equation (10) sionalized by a mass eccentricity eR.
is then written as In this case, eigenvalues of the rotor/bearing sys-
tem (A.= o: + ift) are o: 1 =- 0.0965, P1 = 0.707 OJ8 ,
x(t) = JPr!2 + N 2 sin(OJt + qJx), (11) o: 2 = - 0.0514 and P2 = 0.966 OJ so that peak val-
8,
y(t) = JPr!s + Ns sin(OJt + qJy), ues for the response are near the first and second
natural frequencies. As I o: 1 1 is larger than Io: 21, the
where
qJ =tan- 1 -Pc2 Table 4.8.1. Basic numerical model
x Ps2'
Rotor mass m (kgf) 70
The whirling orbit of the rotor supported by jour- Shaft stiffness k (kgfjmm) 1.428
nal bearings is usually elliptical, and it is represent-
Stiffness ratio of shaft-bearing y 0.904
ed by the semimajor axis A, the semiminor axis B,
the direction of the principal axis and the whirling L:D 0.5
direction. A and Bare obtained from Eq. (10) as
r(t) = x(t) + jy(t) U=F)
= Pe 2 cos OJ t + Ib sin OJ t 100
+ j (Pes cosOJt + Pss sinOJt). (12)
Applying the relation e;o = cos 0 + i sinO to Eq. (12),
we get 10
284
4.8 Unbalance Response of Rotors in Journal Bearings
285
4 Explanations
100 100
~10
.. 10
--
-.::
---- r=5.0
.1 - - (1) Circular bearing =2.0
............ (IO)mp=l/2 =r.o
------ (10mP=2/3 =o.5
- - - (J2)mp=3/4 s; =1.62 =0.2
0 .5 1.5 0 .5 1.5
a Dimensionless rotating speed w/ ws Dimensionless rotating speed w/ w,
Two-lobe bearing Fig. 4.8.8. Unbalance response as a function of y
E =0.6
<I)
-o s; = o.1535
.c..e 1
E
-<
- - (1) Circular beanng
2
.1 ············ (40)mp =0
---- - (40mp = 1/2
- - - (42)mp =2/3
---- (43)mp=3/ 4
0 .5 1.5
b Dimensionless rotating speed w/ ws
Tilting pad bearing (TP 5, LBPJ r : Relative flexibility of shaft
Fig. 4.8.7. Influence of preload factor on unbalance response Fig. 4.8.9. Logarithmic decrement J
286
4.8 Unbalance Response of Rotors in Journal Bearings
...
- __Ys.,.---.....~~ ,_
(!.)
--
u
(!.)
0.5
'~-
comes unstable. The threshold speed for insta-
"0
u bility is large as the preload factor become large.
'§
0
:E 0.1 0.5 w
m, (l)
s
.3
-0.5 r : 0.5 References for Sect. 4.8
s~ : o.3
- 1 : : Dimensionless rotating speed [1] Gasch, R . and PI t; ncr, H ., Rotordynamik, Springer-Ver-
s
lag, 1975, Cha p. 7.
b Logarithmic decrement oof three-lobe bearing ( L: D= 1) [2] lwatsubo, T., Tra ns. of JSME (C), 46-405, May 1980, p . 493.
[3] Kuro hashi, M . e t a!, Tra ns. of JSME (C), 47-422, Oct. 1981,
Fig. 4.8.10. Rotating speed vs. logarithmic decrement p. 1277.
287
4 Explanations
quite different features from that of any kind of Rotational speed (J)
forced vibration, such as unbalance vibration. Fig. 4.9.2. Inertia effect of oil whip
288
4.9 Stability of Rotors Supported by Plain Bearings
Frequency OJK
~ Oil whip
2 (J)K
(J)K
Fig. 4.9.3.
Rotational speed OJ Oil whirl and oil whip
289
4 Explanations
0 0 Co C2 c4 c6 llv
0 10 15
vc = (A 1 A 3 AV{(Ai + A 2 A;- A 1 A 4 As) WK
· (As+ cxA 1 ))) 112 (6) \ 12WK
where
\
1oil
\whip
''A
j,' '
Oil whirl
.5 \ , '- Stable
Vc
wcf(g /Cp) 112
=
(g: acceleration o~ gravity, CP: mean radial
\ I ....... A, L'
290
4.9 Stability of Rotors Supported by Plain Bearings
291
4 Explanations
292
4.9 Stability of Rotors Supported by Plain Bearings
\J:;::/.
6
5
"'
;.
1/ 4
~1
0 0
.02 . 05 .1 .2 .5 .02 . 05 .1 .2 .5 2
A A
10
9
a= 0
12
(4) T hree-lobe bearing
v
(3) Two-lobe bearing
Oil groove 10' x 2 Oil groove 10' x 3
8 L!D=O . S
LID=0.5 ~ 1/16 10
1/lp =0.5
mp=0.75 ~
6 8
1/ 16
~1/4
;. 6
4 1/ 4
2 ~1 2
o~~~~~~~~~~~~~~
. 01 . 02 . OS .1 . 2 .5 1 2 5 10
A
Ja=O
(5) Three-lobe bearing
Oil groove 10' x 3 Oil groove 10' x 4
12 L!D=O.S L!D=O . S
mp=O. 75 mp= O.S
10" 10
8 ~1/16 8
6 6
4 ~1/4 1/4
OL-~~~~~~--~~~~~--~
.01 .02 .05 .1 .2 .5
,\
..\
Fig. 4.9.8. Sta bility charts (1-6)
293
4 Explanations
4 ~ 1/4
where kxy= - kyx
= Q is often assumed. A consid-
..__._ erable part off is still unknown and further inves-
2 • • • •I tigation is needed .
The stability of a rotating system which has such
0 a destabilizing force can be checked by means of the
.005 .01 .2 .05 .2
above-mentioned eigenvalue analysis of the equa-
A
tion of motion derived with due consideration of
Fig. 4.9.8. Stability chart (7) the force. The sign of the real part a; of each eigen-
value s; shows whether or not the system is stable.
Furthermore, the stability margin can be evaluated
by calculating the damping ratio (; or logarithmic
4.9.7 Flow-Induced Vibration decrement <5; as follows:
(; = - a;/(a? + b?) 112 , (18)
Oil whip or oil whirl is caused by the action of the <5; = - 2na;/b;. (19)
oil film in a journal bearing. In addition, modern
high-performance turbomachines often experience As it is usually difficult to estimate the magnitude
another type of self-excited vibration of the rotor, of the destabilizing force fin actual machines, engi-
possibly due to the action of fluid flow in the ma- neers are advised to design machines so as to keep
chine. (; or <5; positive and well above zero at the rated
While the operating pressure of fluid in the ma- load and speed.
chine has increased rapidly in recent times, rotors The following measures are recommended for
have become more lightweight and the operating suppression of the subsynchronous whirl [4, 5].
speed higher. Therefore, the ratio of power to rotor 1. Raising first critical speed by increasing shaft
mass is increasing, and such rotors readily become diameter or decreasing bearing span.
subject to various vibration problems. As will be 2. Enhancing the asymmetry of rigidity of both oil
seen later in Sects. 4.10- 12, steam turbines, pumps film and bearing support in the horizontal and
and compressors have different causes of such self- vertical directions or increasing oil film damp-
excited vibrations. The most effective countermea- mg.
sure against the vibrations in all types of machinery 3. Increasing the system damping by adding
is to improve the damping capability of the rotating squeeze film dampers.
system by designing better journal bearings.
According to recent reports [4, 5], the self-excited
vibration experienced in high-speed turbocompres-
sors are characterized as follows:
1) It depends on load. When load is low it does not
occur at all or the vibration amplitude is negligi-
References for Sect. 4.9
ble. As load increases to a critical point the vi-
bration suddenly starts or the vibration ampli- (1] B.L. Newkirk and H.D. Taylor, General Electric Review,
Vol. 28, No. 8 (Oct. 1925), pp. 559-568.
tude suddenly increases to a violent level. A (2] Y. Hori, Trans ASME, Ser. E, Vol. 26, No.2 (June 1959),
further increase of the load would endanger the pp. 189-198.
machine. [3] E.C. Peste! and F.A. Leckie, Matrix Method in Elastome-
2) The vibration often occurs when the shaft speed chanics, McGraw-Hill, 1963.
[4] J.C. Wachel, NASA CP-2250 (1982), pp. 1- 19.
is higher than the first critical speed and the [5] M. Kurohashi eta!., Proceedings of 3rd Int'l Conf. on Vi-
frequency of vibration is lower than the rotating brations in Rotating Machinery, I Mech E, (Oct. 1984),
speed and near the first natural frequency of the pp. 51-58.
294
4.10 Application to Steam Turbine Generator
4.10.1 Introduction
3000
_i:r-
J.4 poles
-t &7
2 poles,'/ / I
2000 ~I
A significant requirement of a large steam turbine
generator is that it should remain reliable for a long 1000
2 poles r-iV
period to ensure the public electric power supply.
For safe operation of the generator unit over a long
~VJIV lj;V
period the rotor vibration is important. Recently, ~ ~v
the problem of vibration has become more impor- r~v~v 0 3000rpm EL
tant, since turbine generators are larger and their ~ /' C. 3000rpm PAD
)( 3600rpm EL
output higher, in addition to which frequent start
X
-
0
100
0 3600rpm PAD
up and shutdown are necessary in accordance with • 1500rpm EL
the fluctuation in power demand. In addition, im- ~
mediate regulation and stabilization of rotor vibra- X
JV x --EL (2 poles)
;>;: ---·PAD(2 poles)
tion levels is required to shorten the test period, ---EL (4 poles)
which implies high fuel costs. Although there are 10 I I I I I I I I
various kinds of reliable analysis and evaluation 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
methods for rotor vibration, it is necessary to un- Journal diameter in
derstand the rotor supporting conditions and the
dynamic characteristics of the bearing before effec- Fig. 4.10.2. Relationship between unit output and bearing size
tive utilization of these methods is possible. In this
section the application of the dynamic characteris- Table 4.10.1 Main geometry of bearings
tics of bearings in large steam turbine generators Type
will be introduced, with examples of unbalance vi- Item
Pad bearing Elliptical bearing
bration and stability in the rotor-bearing system.
LjD 0.5- 1.0 0.5- 1.0
CjD 1.3/1000- 2/1000 0.9/1000- 1.5/1000
m o- o.s 0.5
4.10.2 Bearing Types and Operating m, preload; C, diametral clearance; D, diameter; L, length
Conditions
generator rotors. Furthermore, the bearing size has
been inevitably increased in proportion to the in-
a) Bearing Types crement in output, as shown in Fig. 4.10.2. The
A large steam turbine generator consists of a num- main dimensions and an outline structure of typical
ber of rotors, as shown in Fig. 4.10.1. The bearing bearings are shown in Table 4.1 0.1 and Fig. 4.1 0.3.
types for this multi-span system are selected with
due account taken of the stability of each rotor and b) Operating Conditions
the load capacity of each bearing. A pad bearing is
generally applied to a high/intermediate pressure The position of a journal centre expressed by an
rotor and an elliptical bearing, to low-pressure and eccentricity ratio and an attitude angle has a pro-
700MW G
295
4 Explanations
nounced influence on the dynamic characteristics of (1) Between pad bearing (2) Load on pad bearing
a bearing, but the position of this centre changes
naturally with the operating conditions, which vary
during the process from start up to attainment of
the rated output of the generator. For instance, the
alignment of bearings varies significantly with the
operating conditions, causing variations in bearing
load, as shown in Fig. 4.10.4. Therefore, each bear-
ing level must be set with prior consideration of the
variations in alignment. Also the rotor-bearing sys-
tem must be designed in such a way as to prevent Oil outlet Oil inlet
Oil inlet
too much variation in a bearing load. However,
some variations will generally take place, amount- <D Carrier ring CD Outer ring
0.8
Relative level changes ( # 4- # 5 )
0. 7
0. 6
"i 0.5
G) 0.4
>
~
"'>
·~
0.3
G) 0.2
~
0.1
"'
E al:;::
o::;:
::>""2- •-L
g=jg_~ 500
>-cn- 440MW
760 4000 400 730mmHg
6 8 10 12
(hr) Fig. 4.10.4. Alignment changes during operation
296
4.10 Application to Steam Turbine Generator
a) Calculation Model
The rotor-bearing system of a large steam turbine
generator has a complicated configuration, consist-
ing of many rotors, discs and blades, operating as a
rotating system. The transfer matrix method is use-
ful for analysis of the vibration of such a rotating
system. The configurations of a low pressure rotor,
a support system, and a rotor-bearing system, are
shown in Fig. 4.10.7 a-c, respectively. The rotating
system is roughly classified into a rotor body and a
bearing system. The rotor body can be replaced
with an equivalent lumped mass and a weightless
0 : Measured values } elastic shaft. The bearing system is divided into a
e : Calculated values 20-in-dia tilting pad bearing bearing section and a bearing support section, each
0 : Measured values } of which consists of a spring and a damping ele-
+: Calculated values 32-in-dia elliptical bearing
ment.
Fig. 4.10.6. a Mean bearing pressure vs. journal centre.
b Rotor speed vs. journal centre b) Support System
The support system is one of the important factors
for determination of the vibration characteristics of
stand the characteristics based on practical condi- a rotating machine. A bearing system can be divid-
tions. ed into a bearing section and a support section, as
illustrated in Fig. 4.1 0.7 b. The support section is
either installed directly on the concrete foundation
or attached to the casing. One example of the verti-
4.10.3 Unbalance Response cal stiffness ratio for each of these cases is shown in
Table 4.1 0.2 with regard to the bearing support,
The admissible maximum amplitude of rotor vibra- rotor, and bearing oil film, and implies that stiffness
tion is in a range from 4/100 to 6/100 mmp-p, in is at almost the same level at the LP section, in each
operating conditions. case. This means that the stiffness of the support
Thorough investigation and quality control are section is also an important factor for determina-
required for rotor vibration at each stage of design, tion of the characteristics of vibration.
manufacture, and assembly, so that the vibration
data of the rotor system extended by 20- 60 m in Table 4.10.2 Stiffness ratio of support elements
total length, weighing 1000-8000 kN, and consist- LP section HP/IP section
ing of 3-6 rotor units may be below the admissible
amplitude regardless of the above-mentioned varia- Rotor 1.00 2.00
Bearing 1.29 2.10
tions in the supporting conditions. Fortunately, ro- Pedestal 1.40 17.9
tor-balancing technique has been greatly improved
297
4 Explanations
ca : Damping
,, ' '"""'~ }'''""''
a b
} Bearing
Equivalent mass
} Pedestal
c
298
4.10 Application to Steam Turbine Generator
lOOr---r--------.-----T--------. Table 4.10.3 Influence coefficient (10- 3 11m/g- mm, degrees)
---- Measured before improvement
- - Calculated
- - - Calculated after improvement
Bearing no. *1 #2 3* #4
Measured 0.04, 3° 0.31, 130° 1.11, 31 ° 1.08, 213°
Calculated 0.17, 3W 0.29, 91 o 0.99, 60° 1.07, 242°
i: position of measurement
j: corrected position
k: balancing speed
(a) i · k = j; [W] = - [a] - l [A]
(b) i · k > j; [aV [a] [W] + [aV [A] = o
W; corrected weight
Balancing on an actual machine can be facilitat-
ed if the influence coefficient has been obtained in
advance by calculating the unbalance response.
Table 4.10.3 shows the comparison between the Fig. 4.10.11. Change of vibration vector
299
4 Explanations
X max. o·
(C/Cc <1)
WI / wn w lwn
,,
Jso·
5 Fig. 4.10.14. Change of vibration vector due to load
4
v
v
• • On the other hand, Fig. 4.10.14 shows the change
in vibration vector due to load for the same rotor-
• bearing system as in Table 4.1 0.4. Table 4.1 0.4
shows that the Q factor is decreased to a large ex-
• • •
tent and the measured value is close to the calculat-
ed value. Furthermore, Fig. 4.10.14 shows that the
change in the vibration vector due to load is also
,, ', '
• small. These results show the effect of a lowered Q
factor. Meanwhile, this unit has been put into oper-
0 ----
- - 5 10 15 20
ation without field balancing.
Q-factor
4.10.4 Self-Excited Vibration
Fig. 4.10.13. Q-factor vs, number of balancing shots
Self-excited vibration of a steam turbine generator
The approximate relation between the number of is of two kinds: oil whip caused by the characteris-
balance shots applied for an actual machine and the tics of bearing oil film and steam whirl caused by
Q factor is shown in Fig. 4.10.13, in which the num- steam force. In addition, unstable vibration due to
ber of shots tends to increase in proportion to the shortage of oil supply to a bearing with a large bore
increase in the Q-factor. This shows the significance can be regarded as an unusual example. The dy-
of the Q factor. namic characteristics of bearing oil film relate to
Various methods of decreasing the Q factor have both cases.
been tested by calculation for the rotor-bearing sys- Some typical cases are introduced.
tem having the characteristics of vibration shown
by the broken line in Fig. 4.1 0.9. a) Oil Whip
As an example of applying the results of the test
to an actual machine, the effect of increasing the Oil whip is a well-known phenomenon; its causes
pedestal stiffness is introduced. The chain line in have b een clarified and preventive measures estab-
Fig. 4.1 0.9 shows the characteristics of vibration af- lished. Accordingly, it now rarely occurs in an actu-
ter the increase in pedestal stiffness, providing an al machine. With sufficient prior investigation into
example of comparison between a Q factor based stability attempts are made to apply a pad bearing
on the measured data and the calculated value, a s with excellent stability to a high-intermediate-pres-
shown in Table 4.1 0.4. sure rotor subjected to a comparatively low bearing
load. Figure 4.10.15 shows the results of subjecting
those bearings used in a large steam turbine-gener-
Table 4.10.4 Improvement in Q factor ator to their stability limit values. A certain design
Before improvement After improvement margin is maintained on the stability limit value for
each bearing, while both operational conditions
Measured 14.5 4,9
and change of the bearing load due to ageing effect
Calculated 13.6 5.6
are taken into account.
300
4.10 Application to Steam Turbine Generator
I
,Y ~(~EL
,
0 GEL
D PAD
a
_, / I. .I
...__ GEL
a 2.0
I
, joo~ ' D
0
'-.._,I PAD
0+---~~--~----~----~----~--~~
0 0.4 0.8 1.2
So=;%(%)'
Fig. 4.10.16. Model of flexible rotor
Fig. 4.10.15. Stability limit of bearings in operating condition
301
4 Explanations
mmp-p Record of vibration meter nomenon is called steam whirl, and is generated by
steam force exerted on the rotor. It is similar to the
= 0.09 oil whip caused by the properties of a bearing and
= 0.06 is characterized by its whirl with the natural fre-
= 0.03
= :::::o Time
(Normal) I==
quency of the rotor-bearing system.
The force which generates self-excited vibration
is caused by the action of a rotor which rotates
I
302
4.10 Application to Steam Turbine Generator
w;
___!_ A 2 +1 k e /k s - 1 0 Xs
0 -1 Ys
=0 (8)
A So A
0 - Pxx
w
+ Yxx +-
'7
W Pxy + Yx y .xb
A A S
0 WPyx+ Yyx - pyy
w
+ Yyy + -I]0 Yb
where substitutions w;
= ksjm, 1J = qjC, and defi- of several high-pressure rotors, supported by six-
nitions q = static deflection of rotor, C = bearing pad bearings, is shown Fig. 4.1 0.25. It can be seen in
clearance, S0 = (pmff.lW) (C/R) 2 , Yii = 2kii · S0 · Cf this figure that stability tends to decrease as kefks
(mg), pii = 2 cii · S0 · C · wj(mg) are applied. increases. That is, the destabilizing force increases
303
4 Explanations
.,
.
"0
.~
4/4 "i5.
E a
"0 ~
Load '"
"'0 !roo
...l 2/4 Amplitude
(( P-P)
7
'
X.o __) $
0 mm
4.11.1 Introduction pots, and fertilizer plants. They are, however, the
focus of debate from the viewpoint of vibration
Turbo-compressors, which include centrifugal and technology, since they are inherently susceptible to
axial flow types, are used in such diverse fields as the vibration of vanes, shafts, and the frame sup-
petrochemical and refining plants, natural gas de- porting the main body. Generally, a rotor-bearing
304
4.11 Application to Turbo-Compressors
system readily causes forced vibration due to the a) For Critical Speed
unbalance force or a self-excited vibration due to
the oil film force in bearing parts. Designers are, The actual critical speed shall not encroach upon
therefore, warned to give adequate consideration to the specified operating speed ranges, the margin of
these properties of the turbo-machine. which is at least 15% below or 20% above the
The lateral vibration system is made up of shaft, continuous speed. The amplification factor shall
impellers, couplings, bearings, etc. The weight of the not exceed 8 while going through criticals. (Values
impeller and coupling, and the rigidity of the shaft less than 5 are preferable.) Previously, undamped
and bearings are the determinants for the critical critical speeds were used to judge whether critical
speed of the system, whereas the internal friction of speed did not encroach upon the specified operat-
the system, the fluid dynamic exciting force induced ing speed range. However, damped natural frequen-
at sealing points and impeller and the reaction force cies and amplification factors obtained by consider-
of the oil film provided for journal bearings, etc., ing the oil film effect of bearings were introduced
have a marked influence on the resonance level and after revision.
threshold of self-excited vibration as a damping or
destabilizing effect. Therefore, control of damping
influence is indispensable for reduction of the am-
plitude in the resonance frequency range, protec- b) For Allowable Vibration Level
tion against unstable vibrations, and avoidance of
resonance. However, quantitative determination of The double amplitude of vibration in any plane
the internal friction and the fluid dynamic exciting measured on the shaft near the bearings shall not
force is almost impossible at present. For this rea- exceed 25.4 x j12000fN flm, where N is the maxi-
son, it must be pointed out that the dynamic char- mum continuous speed in revolutions per minute.
acteristics of the bearing oil film play an important
role.
Particularly in the case of a centrifugal compres-
sor, which is intended for higher speeds and pres- c) For Balancing
sure, the rotor is operated at more than the critical
speed in many cases. In such machines the vibration The rotating element shall be dynamically balanced
amplitude of rotor should be suppressed each time in multiple planes during assembly. The maximum
it passes through the critical speed at startup or admissible unbalance force on any bearing at the
shutdown. In addition, the impeller and the maximum continuous speed shall not exceed 10%
labyrinth seal induce a destabilizing force in these of the static loading of that bearing.
machines. Accordingly, it is necessary to positively
These requirements have the following implica-
suppress the vibration resulting from such a desta-
tions for the design and manufacturing process.
bilizing force by means of a damping action provid-
ed by an oil film. (i) In the design stage the rotor system, including
In this section we discuss the design procedure shaft, impeller and bearings, is modelled to suit
for rotor-bearing systems and the related applica- the vibration system, and damped natural fre-
tion of bearing data, which is mainly relevant to quencies and amplification factors are calculat-
centrifugal compressors. Briefly, "Rotor Dynam- ed by the complex eigenvalue analysis tech-
ics" as described under API Standard [1], which nique. These parameters are then evaluated.
was widely made available as the design and fabri- If these values are not desirable, design
cation standard for turbo-machines, are outlined parameters such as shaft diameter or bearing
first, followed by the analytical method which satis- dimensions are modified.
fies the API Standard. (ii) At the manufacturing and assembly stage, ro-
tating parts are balanced by using balancing
machines. Before release from the workshop, a
test is carried out with the machine running to
check that the vibration amplitude is lower
than the admissible values.
4.11.2 Dynamic Design of Rotor Systems (iii) In test operations at the users' factories, unex-
According to API Standard pected vibration may occur due to the different
load or stiffness of base conditions, or other
causes. Therefore, vibration values are moni-
The American Petroleum Institute Standard, tored and checked again. If unexpected vibra-
API 617, for turbo-compressors, which was revised tion occurs an investigation is promptly insti-
in 1979, specifies the following in the section "Dy- tuted; vibration analysis can be effectively used
namics" : to decide on corrective action.
305
4 Explanations
a) A Simplified Analysis
(
When Eq. (1) is used, the following parameters are
calculated
2m
Qs 542 radjs, w = 403 radjs,
=
X = 1.237, 17 = 0.452.
Fig. 4.11.1. Jelfcott rotor system
Then, substituting these values into Eq. (1),
Q jQs = 0.755, Q = 409 radjs, J = 0.478
and AF = 6.6
are obtained. Alternatively, in the case of Fig.
4.11.3,
k,
kbj(2ks) = x/2 = 0.619, cbwjkb = 11/X = 0.365
give the values QjQs = 0.75, AF = 7. Either may be
k, k, = k, + i we, used. In this example the effect of the oil films low-
ered the natural frequency by 25% compared with
k,k.
-K=---= the simply supported condition, and the amplifica-
ks + k, tion factor is 7. In addition, the optimum damping
b
a value that minimizes AF is given as 11 = (1 + x)/2
Fig. 4.11.2. Equivalent vibration system from Eq. (1).
306
4.11 Application to Turbo-Compressors
307
4 Explanations
tion, while (Q)Qs) and JY are those in y-directional The eigenvalues s and eigen vectors V can be nu-
vibration, which can be obtained from Eq. (1). merically obtained by the QR method from Eq. (6).
As Q increases, J corresponding to forward whirl, For this purpose, the dimensions of shaft, weight
which has a negative sign, decreases, then becomes and moment of inertia of the impeller, dynamic
zero at Qmax' and eventually becomes negative. characteristics of bearings, etc. must be fed in, so
Since a negative J means the occurrence of whirl, that M, C and K can be automatically arranged in
the greater Qmax is desirable. the computer program and that the pair of s and V
An example of the model in Fig. 4.11.1, having can be calculated through the QR method. When s
is a complex value expressed as rx + fJ i the following
m = 148 kg, = 43.4 MNjm
ks
equation becomes applicable:
kxx = 69.8 MNjm, kyy = 5.41 MNjm
x =Vest= Ve" 1 (cosfJt + i sin[Jt). (7)
c xx = 84.2 KNs jm, cy Y = 24.4 KNs jm
The above equation implies a free vibration, as seen
is shown as follows. For x- andy-directional vibra- in Fig. 4.11.4, whose amplitude varies with e"1 and
tion, whose angular frequency is equal to fJ. When rx > 0,
Qx/Qs = 0.806, Jx = 0.531 unstable vibration occurs.
fly /Qs = 0.343, JY = 2.12 In order to express the vibration characteristics,
the logarithmic decrement J = - 2 n rxj [J, amplifica-
are obtained by applying Eq. (1).
Substituting these values into Eq. (4) provides
tion factor A F = J
rx 2 + fJ 2 / ( - 2 a) and natural fre-
quency f = fJ/(2 n) are used instead of rx and fJ.
Bvr = 0.266, Bmi = 0.0946, Whirl directions of vibration and distribution of
amplitude can be judged from the vector V repre-
BDi = 0.0152, qmax = 0.278,
senting vibration mode. An example of the above
and then Qmax = qmax X ks = 12.1. process is shown in Fig. 4.11.5 b.
It is also possible to determine the amplitude
caused by unbalance forces by solving the following
equations
b) Finite Element Method Analysis
Mx+Cx+Kx= f . (8)
When gyroscopic effect, difference in characteristics
between right and left bearings, and asymmetrical Equation (8) can be reduced to the following form:
distribution of mas or stiffness along with the axial [K- M w 2 + iwC] X= U w 2 (9)
direction must be taken into account, or modes of
a higher order must be determined, the simplified where
method is not considered adequate; FEM can then X: amplitude
effectively be used. w: angular velocity
Where a rotor system is modelled as a multi-mass U: vector representing unbalance
point system, the free vibration equation can be
expressed as follows: At the design stage, it will become possible to check
whether the limiting value of amplitude specified in
M X+ c X+ KX = 0 (5) the API Standard is complied with by obtaining X.
where
x: displacement vector
M: mass matrix
C: damping matrix
K: stiffness matrix
all of which are determined through the finite ele- time
ment procedure. M comprises the mass of each
point and the moment of inertia of the impeller. C
involves the gyroscopic effect and the damping co-
efficients at the bearings. K involves the spring coef-
ficients of bearings and the fluid dynamic cross-cou-
pling stiffness, in addition to the stiffness of the
shaft. When X = vest is introduced in Eq. (5), it can
be converted to
(6)
Fig. 4.11.4. Wave form in free vibration
308
4.11 Application to Turbo-Compressors
~·
Rotating speed
bearing characterized by a higher stability is pre-
ferred in the high-speed centrifugal compressor.
Also, bearing diameter D = 76.2 mm, viscosity
10400rpm f.1 = 10.8 MPa · s and clearance ratio of 2 CPI D
Log. decrement
= 2.5/1000 were used.
-0.01 For the purpose of vibration evaluation, an ei-
Nat ural frequency genvalue analysis with reference to the FEM model
59.95Hz
shown in Fig. 4.11.7 was conducted to assess the
(b) Results
rigidity of the shaft. The critical speed under simply
supported condition Q s was found to be 542 radjs
Fig. 4.11.5. Analysis model and calculated results (5126 rpm).
The mark 0 represents the node where mass and
moment of inertia of impeller is taken into consider-
In this chapter only the finite element method is ation. The resulting Qs shows that the rotor is flex-
applied to the analysis, but the transfer matrix ible, since it is found to be lower than the service
method is also available. speed.
Next, the damped critical speed and damping
ratio, which take into account the oil film charac-
4.11.4 Application to a Centrifugal teristics of bearings, are discussed using the simpli-
fied method. The bearing data given in dimension-
Compressor less form has been converted to dimensional form,
as shown in Table 4.11.1.
A centrifugal compressor shown in Fig. 4.11.6 is Assuming that the equivalent mass m is
used as an example of a situation in which the 148 kg = 1/2 the total mass, the equivalent spring
bearing data and the analytical procedure are ap- constant ks can be calculated as ks = m, Q; =
plied. Usually, from the given specifications, such as 148 x 542 2 = 43.4 MNjm. The natural frequency
those involving pressure and flow rate, the machine
type is determined first and then the number of
stages, dimensions of the impellers and speed are
fixed in the performance design. Further, the di- 1960
1510
mension of impellers and consideration of the seal 0
....,_
0 N
;::;
....
N
and arrangement of bearings make the design con- Q I--... t~ .n.
cept of the rotor system more concrete. As a result, t--riM" - ~ ~
it can be said that the major vibration design ~ ' - ...J FJ:I 1.I
parameters are type of bearing, configuration of
sealing parts, and shaft dimension.
The example discussed here assumes the rotor
i wf~ ' --;:?p14
' -------------------- W=1450N
system design of a centrifugal compressor with ser- Q, =542 rad/s
k xx kyy cx x Cyy
speed
13200 53.2 7.23 48.7 22.2
7620 69.9 5.41 84.2 24.4
Suction Discharge Power 3880 95.0 3.54 169 26.9
3000kW 2586 114 2.65 257 28.1
1602 137 1.84 414 29.6
Fig. 4.11.6. Centrifugal turbo-compressor
309
4 Explanations
2nd 1.0
Vertical 2.4
200 1.7 l.S 2.0
0.5
Ci
II
1 -so
_2"N
"'
E
'-. ...u"'
>.
p"'
~ 100
.;,
"c:r 1 st 0
AF=(6.6)
Cross-coupling stiffness
Q MN/m
-0.5
5000 10000
Rotor speed N" rpm
Fig. 4.11.8. Natural frequency, amplification factor vs. rotor Fig. 4.11.9. Decrease in log decrement due to cross-coupling
speed stiffness
and the amplification factor relevant to each speed curve for Q - () in Fig. 4.11.9 falls at the top right of
are calculated on the basis of m, ks and the bearing the graph.
data using the procedure specified in paragraph A verification through FEM should be made in
4.11.3 a.(i), and they are plotted with broken lines, the final design stage, in addition to the simplified
as shown in Fig. 4.11.8. Numerical values in this method discussed above. Results obtained from
figure represent the amplification factor. The solid FEM are plotted in Figs. 4.11.8 and 9. The results
diagonal line indicates the frequency equivalent to of FEM analysis shown in Fig. 4.11.8 confirm that
the rotor speed. Point <D, at which this line crosses in the first mode, the natural frequency and ampli-
the broken line, indicates the critical speed. The fication factor predicted by the simplified method
requirements of the API Standard will be adequate- agree well with those yielded by FEM, while in the
ly satisfied, since the critical speed of 4500 rpm is second mode the natural frequency is equal to
equal to 59% of the service speed and the amplifica- 153% of the service speed. Thus, the provisions in
tion factor is estimated to be 5.8. the API Standard are adequately complied with.
Meanwhile, in the case of a high-pressure cen- The logarithmic decrement under the influence of a
trifugal compressor, the possibility of unstable vi- destabilizing force is shown in Fig. 4.11.9. In this
bration during operation under load must be taken case, although Qmax is a little smaller than that ob-
into account. Whenever a disturbance occurs due tained through simplified analytical procedure, it is
to the fluid dynamic cross-coupling stiffness Q at well proved that Qmax has a sufficient margin, in
the impeller or labyrinth,() is reduced and vibration contrast with Q obtained through the Lund equa-
readily occurs. However, where Qmax specified in tion. Altogether, it can be confirmed that the rotor
paragraph 4.11.3 a.( iii) exceeds Q, the condition of and the bearing system of the present turbo-com-
() > 0 can be maintained, and stability is always pressor have been designed to complete satisfac-
retained. In the present example Qmax is equal to tion.
12.1 MN/m, as shown in specimen calculation (iii). If the compressor in Fig. 4.11. 7 had been
As for Q, however, Lund [6] introduced Q = p T/ equipped with circular bearings or two-lobe bear-
(2 r h) using the so-called Alford force, where T, r, h ings instead of the tilting pad bearings, the system
and pare torque, effective radius of impeller, outlet stability would have been as shown in Fig. 4.11.10.
width of impeller and factor, respectively. Assuming Here the former is designed in accordance with data
that T= 3.68 kN · m, r = 150 mm, h = 10 mm, and no. 2(i.e. 10° x 2, L: D = 1.0 and 2 CP/D = 2/1000)
P= 3, then Q is 3.68 MN/m. As Q < Qmax> the and the latter with data no. 13 (i.e. 10° x 2,
design therefore proves to have an adequate mar- L: D = 1.0, mp = 0.5, p = 170°, 2Cp jD = 4/1000
gm. (2Cb/D = 2/1000)). D, Wand 11 of the tilting pad
However, the destabilizing force is quantitatively bearing are used in these cases without change. As
not sufficiently understood. For instance, there is the oil whip attributable to the oil film characteris-
some possibility of the existence of a negative tics must be considered in these cases, () is plotted
damping force, in addition to the cross-coupling against speed in Fig. 4.11.1 0, referring to the first
stiffness force. Consequently, not only the condition forward mode. [Point ® in the diagram corre-
of Qmax > Q but also an additional margin obtained sponds to that illustrated in the calculation exam-
by elevating() near to Q = 0 is desirable, so that the ple of paragraph 4.11.3 a.(ii).] It appears that an oil
310
4.11 Application to Turbo-Compressors
2.0
----- Simplified method
- - FEM analysis
Self-fexcited
1.0 -::1.. whirl
~ 61Hz
"0
::l
-~ {7
[
'<> P..
;: E
"'8 0
..:X: A Pd =4.9
c
"'... I
Unbalance
""' ¢
'J
Q response
.0 173Hz
0
-1 Pd=5.5 ¢
-1.0
& A
d 0 100 200
Frequency Hz
bearing
Fig. 4.11.11. Relation between discharge pressure and vibration
-2.o ~o.. _ _ _ __.__ _ _ _ ___
0
~...- .....~
311
4 Explanations
312
4.12 Application to Pumps
Q
types: horizontal and vertical installations. Typical
Resin
pump structures are shown in Fig. 4.12.1. Cylindrical bore All pumps
Plastic
0
Q)
.D Lemon bore
As the applications and types of pump vary over an Horizontal and
·a
u
(multilobe)
vertical pumps
extensive range, various types of bearings are ap- "'
G
0 bearing
>,
plied. Further, bearings are made by many manu- Ci Axial grooved
facturers and almost all types of bearings are used. Grooved bearing Circumferential Vertical pumps
Typical bearing types installed in pumps are shown grooved
in Table 4.12.1. The types most frequently used in Spiral grooved
0
medium-sized and small pumps are ball and roller Tilting pad Number of Horizontal
bearings, while journal bearings are mostly used in bearing pads varies pumps
Hydrostatic bearing
0 Pocket type Sodium pumps
Radial bearing
Discharge pipe
313
4 Explanations
large pumps. However, journal bearings are also tional conditions. Therefore, estimation of the force
applied in pumps used for liquids precluding the is difficult at the design stage. Another feature of
use of ball and roller bearings. pumps is that water-lubricated bearings are used.
Various configurations of bearings are selected There are three major differences between water-
according to the desired applications from the view- lubricated and oil-lubricated bearings.
points of performance, cost and maintenance. Hy- 1. Turbulent flow can easily be generated in water-
drostatic bearings are used in special cases. In the lubricated bearings since water has low viscosity.
selection of any bearings the emphasis is on operat- 2. Environmental pressure can be higher than at-
ing experiences.
mospheric pressure when used internally in a
In addition, liner ring seals used in pumps can
pump.
also be considered as bearings from the aspect of
3. The virtual mass of water can have a negligible
vibration characteristics [2 -4]. influence in oil-lubricated bearings, in relation to
point 1 above [7].
All these points can affect both the static and the
4.12.3 Differences Between Pumps dynamic characteristics of the bearings. They will
and Other Rotating Machines be discussed in later sections.
314
4 .12 Application to Pumps
{
F min F max Horizontal shaft - reaction force due to rotor mass
Rigid coupling - misalignment force
NO
hydraulic radial force is usually given by the equa- Imaich's research [5] and can be regarded as typical
tion below. The dimensionless radial force used in data of single-volute casing pump. In the figure, the
the equation is obtained experimentally value shown by Stepanoff is frequently used as a
conservative design evaluation. As shown in
F,. = Kr · y · H · D · B'. (1)
Fig. 4.12.3, Kr becomes 0.36 near the shutoff point
where in a pump of this type. Therefore, to reduce the
radial load, balanced casings are used. Typical bal-
F,.: Hydraulic radial force anced casings are double-volute and multivane dif-
Kr: Dimensionless radial force coefficient fuser casings. As shown in Fig. 4.12.4 [9], the radial
y: Specific gravity of pumping fluid forces in these casings are significantly smaller than
H: Head of impeller under scrutiny in the single-volute casings.
D: Impeller outlet diameter Moreover, the direction of the radial force is sim-
B': Impeller outlet width (including shrouds) ple in single-volute casings, as illustrated in
Fig. 4.12.5, while in double-volute and multivaned
Some research workers have investigated hy- diffuser casings the direction exhibits a complicated
draulic radial force. Figure 4.12.3 derives from variance, as shown in Figs. 4.12.6 and 4.12.7.
315
4 Explanations
0. 40 ..----r---.----.,..- -..,--""T""---r---- ,
-+-Tongue I
(Static component)
--...-Tongue III
0.351--=...d--+- (Static component)
--o-- Tongue I
(static and dynamic) max.
0. 30 1-..I..L:.ol-~-fr- - - - Stepanoff
Volute
Tongue I III
Fig. 4.12.3. Radial force coefficients in single volute casing pump [5]. "Tongue" refers to a tongue-shaped insert in the casing of the
pump in this experiment
0. 40
0. 35
o Double-volute casing
• Multi vane diffuser casing
0. 30 90"
0. 25
:..: 0. 20
0 . 15 0.792
0. 898
0.982
0. 10 1.036
1.095
1.171
0.05
0
~
p= ~ !'....
0 0. 2 0. 4 0. 6 0.8 1.0 1.2 1.4
Q / Qn
270"
Fig. 4.12.4. Radial force coefficients in double-volute casing Fig. 4.12.5. Radial force vectors in single-volute casing pump
pump and multivanes diffuser casing pump
The greater part of the hydraulic excitation force, volute and multivane diffuser casings [9]. Typical
which plays an important role in the excitation results are shown in Fig. 4.12.8. For single-volute-
forces acting on pumps, has not yet been clarified, type casings, the data obtained by Imaichi et al. is
and few research studies have been conducted. presented in Fig. 4.12.11. It should be noted that
Kanki et al. [9] have obtained the overall values of different models and ordinates are used in
the hydraulic excitation force, which varies as the Figs. 4.12.9, 4.12.10 and 4.12.11. In these figures, Q,
flow rate changes, from model pumps with double- represents the rated flow rate. They can only be
316
4.12 Application to Pumps
270"
270"
Values noted in the figure show flow rate Values noted in the figure show flow rate
Fig. 4.12.6. Radial force vector in double-volute casing pump. Fig. 4.12.7. Radial force vector in multivane diffusor casing
Values noted in the figure show flow rate pump. Values noted in the figure show flow ra te
"'
0
"
DiffUSER CASING
i
NUMBER Of IMPELLER VANES : l : 7
u
,;:
: CENTER
,.
; ECCENTRICITY 0.20
~
i5
~~~
<
"'
1 MINUIT l
"'
"'
0.05 ~
z:
0
;;;
z:
0
0
1171
20
UJI
40
ISOJ
60
1671 1831
80
11001
100
[1171
120 "i5
"'0
:I:
i
DOUBLE VOLUTE CASING
~
~ u
0 0.1$ ~
g< 0
u
..,..
u 1 MINUIT{ ~
V>
:ii :::>
z: "'
~
0 0.10
~
50 <
i5
<
"'
Vl
0.05 "'
w
~
z:
Q
"'
~
00 100 120
0 "i5
20 40 60
317
4 Explanations
IIl
0
LOW CYCLE COMPONENT =£
40r---. .----------------------~ "
-<>DOUBLE VOLUTE+ 7 c.';
IMPELLER Z= 6 O.O II
+DOUBLE VOLUTE+
"- IMPELLER Z= 7 0.06 ~
6 w
~ 301-------\-----,-6-DIFFUSER CASING+ u
~ IMPELLER Z= 6 r;:
w 0.05 I'.
0 _..DIFFUSER CASING+ w
::::> 0
E-< IMPELLER Z= 7 u
:::1
p..
::a 20 1-----+-----\1----.------.-----r-----i
< 0 0.2 0. 4 0.6 0.8 1.2 1.4
0 Q! Q.
<
0
....l Fig. 4.12.11. Examples of dynamic hydraulic radial force mea-
sured in single-volute casing pump
Table 4.12.3
NZ COMPONENT
40 -o-DOUBLE VOLUTE+ IMPELLER ~
z= 6 0.07 k;
II
2n n
Negative pressure part
Water-lubricated bearing
Ex. Normal bearing
Oil film seal (compressor, etc.)
318
4.12 Application to Pumps
0 so
~~--
0.2 .,
Non
v
Cavi;/ ... 360" Brg.
\
"' 1/ , ~·
v ,'
0 0. 4 LID=l.O
·c:;
·.:
c
"
(.)
(.)
0. 6 ~' Cavi. ~
t:Ll ~.,
'
·~ r.;;: Cxx -:.,.;:;7 7
"'
0.8 /_ Kxy
~, ...
' v
c'
f...- ~
~
."-. ~Cyy /
........... ~"--
1.0
10 - 3 ro-' ro-• '' ~
s Ky~ ~ a, ~
Sommerfeld Number p.N(R /C)'
P,.
..)
~
.""'-,
Vb~ ~
lj.
Fig. 4.12.12. Eccentricity versus Sommerfeld number for 360°, Kxxlli.._L
L : D = 1 bearings \
Cxy ~Kyx
1
\ I
I
0.5 II
0 0.2 0. 4 0.6 0.8 1.0
Eccentricity <
100
a bearing centre
b rotation
.\ 360" Brg.
LID=l.O
c direction of eccentricity
319
4 Explanations
Black obtained the following equation for dy- flow (0 for zero rotation and 1 for high-
namic coefficients: speed rotation)
C =Radial clearance of annular seal
A
nRP FY
[Fx] L = Axial length of annular seal
Ra = Axial flow Reynolds number in annular
seal
= -[flo-tfl 2W2 T 2
- 2flt w T
t~ 1 WT
flo - 4fl2 w T
2 2] [X]
Y
R, = Circumferential flow Reynolds number in
annular seal
T= Average fluid passing time
_ [ flt T fl 2w T 2] [xy:] w = Angular velocity of shaft rotation
- fl2W T 2 - flt T
J [x]
R = Seal radius
_ [fl 20T
2 P = Seal differential pressure
fl2 0T2 ji . (2)
flo, flp fl 2 are given in Fig. 4.12.16 for inlet loss of
In Eq. (2), A is the friction coefficient for axial c; = 0.5, and P= 0, 0.5, 1.0.
flow, and therefore A varies with liner ring configu- The analysis mentioned above is performed on
ration and Reynolds number. As for annular seals, the basis of short bearing theory. The following
a semi-empirical equation proposed by Yamada approximations are therefore given to take account
et al. [16], in which the Reynolds number for cir- of seal length, which has not been sufficiently con-
cumferential flow is taken into account, is used to sidered hitherto [2]:
determine the coefficient:
A = 0.079 R;; 114 { 1 + (7 /8 R,/ RY} 318 (3) flo
(RL) = 1
flo
+ 0.28 (L /R) 2 '
a= A(C/L) = Parameter of resistance for axial
flow flt
(RL) = flt
1 + 0.23 (L/R) 2 ' (4)
P= (7 /8 R,/Ra) 2 / {1 + (7 /8 R,/RY} = Parameter
of the ratio of axial and circumferential
fl 2
(L)R=
fl2
1 + 0.06 (L/R) 2
-
0
0.8
~ I--
1.0
0. 7
0.6
v
v ..........
/3 t-- 4.12.5 Guidelines for Evaluation
0. 5 v 1-- 1--
0. 5
1-'o
0.4 / ./
v ..,_.. 1.0
Recently, various analysis programmes have been
I ~v -- - -- developed for shaft system vibration, and vibration
0.3
0. 2 ~, ,-
,
-- -- constant A
analysis can be performed easily. However, in
0.1 Ill pumps, few shaft systems exhibit high sensitivity to
vibration except some high-energy pumps, and
0
II' highly accurate vibration analysis need not neces-
2 3 4 5 6 7 8 10
9 11
sarily be applied. Moreover, in most cases vibration
I
on the top part of the motor is attributable to the
2.0 constant A vibration characteristics of the structure system. On
--p~
l 0 the other hand, in the development of a new high-
~ I-- - ~/3
1.5
,_,, f-- energy machine, vibration analysis of the shaft sys-
- ~ 1---
0.5
1.0 1.0 tem is indispensable and the method of study de-
~ ~ :::::: ~
F::::: 1--
0.5 ~
~
scribed in the previous section must be followed in
~ F" detail [10, 12].
0 2 3 4 5 6 7 8 9 10 11 Examples of calculation models and evaluation
of the results are described. In practical terms
1. 0 r: t--1.0
0.5~/3
judgement entails consideration not only of the cal-
culated result but of the actual operational data,
Ill +-- 0
!-',
0 .6 ~ '::"""- applied data.
~...
320
4.12 Application to Pumps
Fig. 4.12.18. Critical speed separation margin with due consid- Stability can be a problem in a high-speed, high-
eration of damping effect pressure pump. The annular seals, the impeller, and
• "' .........bc=I.O%
2C
0 6 ····· ····D,=4 .6%
v ----- ,
?- 100
c. ......,, Diffuser
s 80 I
3 I
150
60
<::
s" 40
I
"u
"'
c.
V)
cs
20
0
l u
0~ 100
0 20 40 60 80 100 120 3
a Flow rate Q% ~
"sc.
E
·g.
"0
50
;.:::l
0.5
FULL 0 4
JOURNAL 0.3
co
BEARING 0.2
0.1
..,. 0.05
.g
8
~ 0.025
·o.
6
Q 0
\
Unstable Stator mode \
-0.025 in y direction \
\
- 0. 05 ';;-----;-f;;;:;;;-----;:;f;;-;;----;=:-----'''-;;t;;;;--~~
0 1000 2000 3000 4000 5000
SHAFT SPEED(RPM)
a
0.1 10
Sommerfeld number ( S )
n, : First natural frequency of the system
Fig. 4.12.20. Stability chart of various bearings
(MOTOR)
1 1000
2" 500
the oil film bearing can all be sources of destabiliz-
~ 750
ing forces acting on the pump. Figure 4.12.19 shows 6
an example of a self-excited vibration that occurred 3
f-o
t'
in a vertical pump. This is a complex phenomenon z
~ 500
contributed to by the temperature of the pumping ::E
~
fluid and impeller radial clearance. The stability of u
<
the pump is generally evaluated from a stability ....l
ll.. 250
chart such as that shown in Fig. 4.12.20 for the Ul
5
bearings and by conducting complex eigenvalue
analysis. Furthermore, Fig. 4.12.21 shows examples o~lt-~-----~~---~~-_j
0 500 1000 1500
of measurement and analysis with consideration for (n,) ( 2 nJ
the oil whip developed at natural frequencies of the b SHAFT SPEED (RPM)
casing [15).
Fig. 4.12.21. a Complex eigenvalues calculated for vertical
pump. b Vibration change as a function of speed motor solo
operation
Important points relating to the design and analysis References for Sect. 4.12
of a pump shaft system with journal bearings have
been described. These can be summarized as fol- [1) Makay, E. et al., Survey of Feed Pump Outages, EPPI
lows: Report, FP-754, (1978).
[2) Black, H.F. et al., Effects of High Pressure Ring Seals on
1) Bearings for hydraulic machines such as pumps Pump Rotor Vibrations, ASME Paper, No. 71-WA/FE-38.
must be designed with larger margins than for [3] Manabe, N. et al., A study on the Dynamic Characteristics
machines treating gas and electric machines. of Seals in Pump, Misubishi Juko Giho, Vol. 16, No.6, (in
Therefore, it is necessary to give sufficient con- Japanese) (1979).
[4] Takagi, M. et al., Analysis and Design of Centrifugal Pump
sideration to the actual data of similar machines, considering Rotor Dynamics, Proceedings of I. Mech. E,
since hydraulic excitation force is difficult to pre- (1980).
dict accurately. [5] Imaichi, K. et al., Experimental Study on the Radial Thrust
2) The following three items are important for fu- in Centrifugal Pump, Transaction of JSME Vol. 27, No.
ture works. 294 (1971) (in Japanese).
[6] Stepanoff, A.J., Centrifugal and Axial Flow Pumps. (1957)
• Hydraulic radial force John Wiley & Sons.
• Hydraulic excitation force [7) Agostinelli, A. et al., Trans. ASME, Ser. A, 82-1, (1960-4),
• Characteristics and effects of annular seals 120.
322
4.12 Application to Pumps
[8] Biheller, J.J., Trans. ASME, Ser. A, 87-3, (1965-7), 319. [13] Shiraki, K . eta!., A New Vibration Criteria for High Speed
[9] Kanki, H. eta!. , Experimental Research on the Hydraulic Large Capacity Turbomachinery, Proceedings of the
Excitation Force on the Pump Shaft. ASME Paper, 81- Eighth Turbomachinery Symposium, (1979).
DET-71 , (1981). [14] Manabe, N., Vibration of High Speed Pump in Low Flow
[10] Matthew, C.E .K, Solution of the Subsunchronous Whirl Region Operation, Turbomachinery, Vol. 8, No. 10 (1980)
Problem in High Pressure Hydrogen Turbomachinery of (in Japanese).
Space Shuttle Main Engine. AIAA/SEA 14th Joint Propul- [15] Tomisawa, M. eta!., Vibration Analysis of Vertical Electric
sion Conference, (1978-7). Motor, Mitsubishi Electric Giho, Vol. 55, No.3 (1981) (in
[11] Kanki, H. eta!., Experimental study on the Dynamic Char- Japanese).
acteristics of Pump Annular Seals, Proceedings of Rotor [16] Yamada, Y. , Flow Resistance of Rotating Cylinder, Trans-
dynamic conference I. MechE. C297/84, (1984). action JSME, Vol. 27, No. 180 (1961) (in Japanese).
[12] Manabe, N. eta!., Development of Super High Pressure [17] Kaneko, S. eta!., Static and Dynamic Characteristics of
Centrifugal Pump for Water Jet Bedrock-Crushing, Mitsu- Annular Plane Seals, Proceedings of Rotor dynamics con-
bishi Heavy Industries, Ltd. Technical Review 1981. ference I. Mech E. C278/84, 1984.
323
R. E. Roberson, R. Schwertassek
J.P. Gupta
Approaches to
Modeling of
Friction and
Wear
Proceedings of the Workshop on the Use of Surface
Deformation Models to Predict Tribology Behavior,
Columbia University in the City of New York,
December17-19, 1986
1988. 59 figures. XIII, 173 pages. ISBN 3-540-96656-0