EHL Military Paper
EHL Military Paper
EHL Military Paper
AD NUMBER
AD877768
FROM
Distribution authorized to U.S. Gov't.
agencies and their contractors;
Administrative/Operational Use; NOV 1970.
Other requests shall be referred to Air
Force Aero Propulsion Lab.,
Wright-Patterson AFB, OH 45433.
AUTHORITY
ELASTROHYDRODYNAMIC
0LUBRICATION
J. M. McGrew
A. Gu
H.S. Cheng
CD S. F. Murray
November 1970
U~~uuLL
NOTICE
When Government drawings, specifications, or other data are used for any purpo.c
the United States Governmen,: thereby incurs no responsibility nor any obligati'on
whatsoever; and the fact that the government may have formulated, furnished, or i,.
any way supplied the said drawings, specifications, or other data, is not to be
foreign governments of foreign nationals may be made only with prior approval
of the Air Force Aero; Propulsion Laboratory (APFL), Wright-Patterson Air Force
IN IVT_ ..
..... .
log iF SECTIO 1
IT
Copies of this report should not be returned unless return is required by security
DISCLAIMER NOTICE
-iHcLq-c,_;"7
ib)A wcj
BEST
AVAILABLE COPY
AFAPL-TR-70-27
ELASTROHYDRODYNAM IC
LUBRICATION
J. M. McGrew
A. Gu
H.S. Cheng
S.F. Murray
November 1970
968 Albany-Shaker Road, Latham, New York under USAF Contract No.
3048 Task No. 304806. The work was administered under the di-
February 1970.
ii
r
ABSTRACT
topography.
iv
TABLE OF CONTENTS
Page
-f
TABLE OF CONTENTS (Continued)
Page
a. Approximate Theorieg 77
b. Numerical Solutions to the Energy Equation
c. Temperature Measurements 81
5." FRICTION 85
a. Rolling and Sliding Friction ........................... 85
b. Spinning Friction' . 106
6. STRESS 111
a. Calculation of Contact Stresses Without Traction 113
(I) Geometric Definitions-- - - - 113
(2) Theory for Line Contact---------------------------- 115
(3) Theory for Elliptical Contict 119
(4) Controlling Stress for Fatigue 119
b. Calculation of Contact Stresses with Traction 121
7. SURFACE TOPOGRAPHY 6
vi
TABLE OF
CONTENTS
(Continued)
A
b. Ball Bearing,')Dynamics 134
(1) The Concept of Raceway Control 134
(2) Ball Bearing Solutions .......----
........- 137
9. ROLLING ELEMENT BEARING FAILURE MODES 142
*a. Classification of Failure Modes for Rolling
Element Bearings 142
b. Contact Failures 144
(1) Wear Failures 144
(a) Mild Wear Failure Mode 144
(b) Severe Wear Failure Mode' 151
vii
TABLE OF CONTENTS (Continued)'
Page
b. Results 230
V. CONCLUSIONS AND RECOMMENDATIONS 234
1. CONCLUSIONS 234
2. RECOMMENDATIONS 237
VI. REFERENCES /
APPRFEENCES
APPENDICES..
-----------------------------------------------------------
. . . . . . . . . . ._ o
,
viii
!
I
LIST OF FIGURES
Page
- ix
KT
LIST OF FIGURES (Continued)
Page
x
LIST Or FIGURES (Continued)
Page
xi
LIST OF FIGURES (Continued)
Page
91. Real Area of Contact vs. Specific Film Thickness ---------------- 187
92. Coefficient of Friction vs. Shear Rate Parameter for Heating
Parameter = 5 x I0 7- - - - - - - - - - - - - - - - - - - - -- - - -- - -- -- - - - - - - - - - - - - 189
93. Coefficient of Friction vs. Shear Rate Parameter for Heating
Parameter = 10-6 ----------------------------------------------- 190
94. Coefficient of Friction vs. Shear Rate Parameter for .=eating
Parameter = 5 x 106------------------------------------------- 191
95. Coefficient of Friction vs. Shear Rate Parameter for Heating
Parameter 10----------------------------------------------- 192
96. Coefficient of Friction vs. Shear Rate Parameter for Heating
Parameter = 5 x 10-------------------------------------------- 193
97. Coefficient of Friction vs. Shear Rate Parameter for Heating
Parameter = 5 x 10------------------------------------------- 194
98. Coefficient of Friction vs. Shear Rate Parameter for Heating
Parameter = 10----------------------------------------------- 195
99a. Velocities with Respect to Bearing Centerline ------------------- 202
99b. Relative Velocities with Respect to Rolling Element ------------- 202
100. l00 vs. o for Lubricants Tested in the ASHE Report ------------- 203
101. Computer Input, Example I -------------------------------------- 206
102. Computer Output, Example -------------------------------------- 207
103. Computer Input, Example 2 -------------------------------------- 216
104. Computer Output, Example 2 ------------------------------------- 217
105. Computer Input, Example 3 -------------------------------------- 222
106. Computer Output, Example 3 ------------------------------------- 223
107. A Contact Ellipse---------------------------------------------- 228
108. Frictional Torque as a Function of Spinning Velocity --------- 231
109. Effects of Spinning Velocity on Frictional Forces and Torque --- 232
xii
LIST OF TABLES
Page
IV. Typical Values for Isothermal Balk Modulus at 25C and 0 psig 37
XXI. Sample Bearing Number 3, Low Speed Angular Contact Bearing .......- 221
xiii
NOMENCLATURE
A2( 1 + -. - in-1
x Rx2
B + 1
2Ry I 2Ry 2
E E
e Protrusion width , in
e
e b
f Coefficient of friction
2
G Shear modulus , lb/in
xiv
NOMENCLATURE (Continued)
h
h*
h Film thickness , in
h* Film thickness near the exit where pressure gradient vanishes, iii
k Wear coefficient
p Pressure , lb/in 2
, lb/in 2
PHZ Maximum Hertzian pressure
0f
o 1 2-
2
16KR
,02oU2
24.
Kof
.. Qm K fT
R Radius of ball , in
xv
NOMENCLATURE (Continued)
T Temperature, OR
T Inlet temperature, OR
UIU
2 Rolling speed at contact of body number 1 and number 2, in/sec
U
E' R'
x
ur Rolling velocity of race, in/kec
us U2 - U1 I n/sac
vV + V2 in/sec
2
VSpeed perpendicular to the direction of rolling, in/sec
E'R'
x
xvi
I!
T
- NOMENCLATURE (Continued)
b
x Contact zone dxit coordinate, in
*. Temperature-viscosity coefficient, OF
1211 U
AE'R'
2 x
w
xvii
NO4ENCLATURE (Continued)
"I
s Spinning angular velocity of ball, rad/sec
xvi1
xviii
.4.1
SECTION I
INTRODUCTION AND SUIMMARX
Recent experimental daca show that highly loade'd angular contact bearings
using high temperature synthetic lubricants have much shorter fatigue lives
than classical bearing theory (Lundberg-palmgren) predicts. Th~se failures
appear to be due to the inadequate lubricant film in the contact zone. This in
turn permits metal to metal contact of local surface asperities, progressive
surface damage, and ultimate failure. Thus successful analysis of such applica-'
tions must include meaningful'prediction of Lhe separating film thickness. This,
in turn requires a thorough knowledge of elastohydrodynamic lubrication and the ,,,
parameters which govern it,. This report is addressed to such predictions.
This problem has been under study for some eighteen years-both here and abroad.
Gradually considerable insight has been achieved into elastohydrodynamic be-
havior. Unfortunately, perhaps because of the diverse aspects of this re-
search, there has been little work done in compiling and interrelating the
various results obtained and published. A critical study of the current itate-
of-the-art is both important and timely. It permits optimum utilization of
available knowledge and improved planning for continued research. The first
result is particularly important. si'Ice there is real immediate need for im-
proved predictions of safe operational limits of load, speed, temperature and
"* lubricant for high performance rolling contact bearings. Longer range, as
additional information is accumulated, still better techniques for analyzing
these variables can be evolved.
This report covers the work conducted during the first phase of a threerphase
program in elastohydrodynamic research as follows:
I.I
Phase II - Based on the experimental and analytical needs determined as
the conclusions of Phase I, experimental programr will be carried out
determining which factors are significant in defining behavior in actual
contact zones, critical lubricant properties, and the interaction becween
the two. In addition, any necessary extensions of analytical descrip-
tions of contact zone hydrodynamics and the influence of contact zone
dynamics on overall ball-bearing dynamics will be carried out.
With these input data the program calculates the following quantities success-
ively by means of a series of subroutines:
1. The nomil film thickness based on isothermal theory (the nominal
--ft ickness is the film thickness at the entrance region of the
elattohydrodynamic contact).
2.- The reduction of film thickness due to thermal and side leakage
effects.
3. The exit protrusion width and depth.
/ 4. Percentage of area of contact.
5. Contact friction.
6. Pressure distribution.
7. Mid-film and surface temperatures.
8. Subsurface and surface stress distributions.
2
I
!j
I.T The use of sulroutines provides a convenient way to replace any old methods
of caiculation with any new techniques which may be developed in the future
-- for any one of these variables. Using this approach, one does not have to
reurite the main cowputer program as the state-of-the-art of elastohydrody-
namic lubrication advances.
The computer program is primarily intended for use on rolling element bearing
contact problems. However, it has been made general enough so that it is
applicable to many other heavily loaded concentrated contact problems. Thus,
it should also find use in the design of gears, cams, and pivots.
t3
.. f"/
SECTION II
STATE-OF-THE-ART IN ELASTOHYDR DYNAMIC UBRICATION
The theory of heavily loaded concentrat.ed contacts has been studied extensive-
ly in recent years. It is an interdiszip'inary theory involving the physics,
chemistry, mechanics, and metallurgy of the interacting surfaces. As such,
it is not an easy problem to grasp.and reduce to practical engineering pro-
cedures.
laid for a preliminary design procedure which is discussed in Section III and
attention is focused on those remaining problems which have not yet yielded
to the continuing research effort in this field.
This chapter has two main objectives. The first is to present engineers with
a simplified picture of the 'hemical and physical behavior of a variety of
petroleum and synthetic lubricants, and the second is to summarize the lubri-
cant properties which are 'known to be important for elastohydrodynamic lubri-
cation.
a. Development of SyntheticOils
4
44
jI
ties of the oil such as viscosity, density, etc. However, they do
have a strong influence on oxidative stability since many of these
compounds are naturally-occurring oxidation inltibitors. They also
have a marked beneficial effect on the boundary lubrication charac-
teristicc of the oil.
The average molecular weight of zhe compounds in the oil will deter-
mine the viscosity grade, while the predominant molecular structures
will govern its physical characteristics.
The need for high performance lubricating oils is a fairly recent de-
velopment. Up to the time of World War II, conventional petroleum
oils were adequate for most machinery. Although synthetic hydrocar-
bons and polyalkylene glycols were commercially available before 1940,
their use was restricted to very speciqlized applications. Even the
petroleum oils that were being used at that time were relatively
crude forerunners of the highly compound motor oils that are availa-
ble today. Many of these oils contain ten to fifteen percent of syn-
thetic compounds, added to improve viscosity characteristics and
chemical stability.
During World War II., the synthesis of new fluids as substitutes for
petroleum oils was given considerable attention, particularly in
Germany. Some work was also done in the United States, mainly on
instrument lubricants, but this effort was limited because there was
no pressing need for replacing petroleum oils. 'The development of
the aircraft gas turbine engine for military use marked the turning
point where the favorable properties of synthetic oils began to over-
shadow other factors, such as cost and a lack of background experience.
The first aircraft turbine engines were lubricated with grade 1010
petroleum oil. Problems with meeting the -65F low temperature start-
ing requirement caused a shift to a grade 1005 petroleum, but this
oil was excessively volatile and had very low viscosity at normal
engine operating temperatures. It was this low temperature problem,
and the need for fluids with lower volatility characteristics,
which really sparked the interest in synthetic oils.
5
cant, specification MIL-L-7808. Later modifications of this oil
(now MIL-L-7808G) include esters other than the diesters. This is
classed as a Type I oil. Problems were anticipated with this oil in
many turbo-prop and heLicopter transmissions because the viscosity
at the normal operating temperature was not thought to be high enough
to provide adequate lubrication for heavily-loaded gears. A com-
promise lubricant, he MIL-L-23699 oil was developed bM relaxing the
low temperature starting requirement from -65 0 F to -40 F. This per-
mitted the use of a higher viscosity ester-base fluid which was
thought to bea more effective gear lubricant. This is classed as a
Type II oil.
All of the oils described above are based on synthetic esters. How-
ever, there are also many other classes of synthetic oils which have
also been studied extensively. Some of these are listed in Table I.
It is difficult to rank these fluids on the basis of physical proper-
ties because small changes in structure, or the use of additives, can
change the ranking appreciably.
Most of these synthetics are now being used for specialized applica-
tions, but each has certain disadvantages which have prevented them
from gaining wide acceptance. Of these fluids, the silicones and the
polyphenyl ethers are probably of the greatest current interest for
high temperature bearing usage.
In the last few years, there has been a growing awareness,-of the im-
portance of elastohydrodynamic (EHD) effects on the lubrication of
rolling contact bearings and gears. Under the extremely high unit
pressures which are generated in the contacyareas, the lubricant
viscosity can increase many-fold, and can thus support appreciable
loads. However, the oil will also be subjected to high temperatures
and high shear stresses, both of which tend to negate the increases
in viscosity with pressure. The combined effect of temperature, pres-
sure and shear on lubricant viscosity is still imperfectly understood.
In addition, there are other physical characteristics of the lubri-
cant such as: density, compressibility, dissolved gases, etc., that
may play an important role in elastohydrodynamic lubrication.
6
TABLE I
CANDIDATE SYNTHETIC LUBRICANTS
VISCOSITY
LUBRICANT rEkPERATURE THERMAL OXIDATION BOUNDARY
-TYPE CHARACTERISTICS STABILITY RATE LUBRICATION
Dimethyl and
methylphenyl A B B D
silicones
Halogenated BEA B B B
silicones
Fluorinated
and
h nd D C-A B-A B-A
chlorinated
hydrocarbons
Synthetic
hydrocarbons
and C-B B C* B-A*
super-refined
mineral oils
Polyphenyl D-C A C C
ethers
7
b. Basic Relationships Between Structure and Physical Properties
The forces that hold these molecules together to form a liquid are
basically electrical in nature. Even though each individual molecule
is electrically neutral, shifts in electron density can and do occur
within the molecule, and this results in stray fields of forces which
induce unbalanced forces in neighboring molecules. This is analogous
to the situation where two magnets of equalpole strength are arranged
as shown below.
S
N
8
of electrical force can set up large attractive forces for neighboring
molecules. This effect is much more important if the molecules do not
have a symmetrical chemical structure. For example, if this molecule
with the match stick shape is symmetrical except for a different group
"* of atoms at one end (polar group), there is a built-in electron unbal-
ance (dipole moment) in the molecul-, which is much more potent than
the induced fields of force (induced dipole moment) described above.
" Because these electrical forces fall off very rapidly with dis-
tance, liquids have only short range order. Therefore, liquids
can respond to any type of shear stress.
9
/
valence bonds. The simplest member of the series is a gas, methane,
which consists of one carbon atom surrounded by four. hydrogen atoms.
H-C -H
The next member of the series, ethane, has two carbon atoms:
H H
H A ____- H
H
I T H
SHt
Some of the physical properties of this series are listed in Table II.
At room temperature and atmospheric pressure, the short chain members
of the series (up to 4 carbon atoms in the chain) are gases. This is
because the volume and surface area of the molecules is too small to
exert attractive forces large enough to keep the molecules together
unless they are forced together under pressure. If the temperature is
lowered so that the kinetic energy of the molecules is reduced, then
the gases will condense to the liquid state.
10
V
LTABLE II
SOME NORMAL PARAFFIN HYDROCARBONS
4olecular
Formula Name Boiling Point Melting Point Density at 20C
oc oC 2/cc
(at 15 mm)
11i
VI
molecules which contain more carbon atoms in the chain are more vis-
cous and have a higher resistance to flow. As the temperature is
raised, the kinetic energy of the molecules increases and this ther-
mal agitation reduces the forces of attraction, thus decreasing the
viscosity of the liquid.
-C - C - 0 - C - C - C - C - C-
tether linkage
H H H
HC CH HC NC " H
HC( , .CH 4ZHC,
H H H
Benzene Naphthalene
12
These aromatic compounds have higher attractive forces between mole-
cules than the paraffinic-base oils and, as will be noted later, this
influences many of their physical properties.
-OH alcohol"group
-COOH acid group
-NH2 amine group
-CONH2 amide group
These groups could be atta 'ied to a carbon atom at any point on the
hydrocarbon chain (by substituting the group for a hydrogen atom).
However, the.greatest benefit in boundary lubrication is achieved by
having the group attached at one end. For example, stearic acid has
the following struc ure:
H HI H
H-C - C -d - C - COOH
'14
portion stands up away from the surface. The adherence of the layer
depends on two factors, the strength of the dipole moment which is
pulling the molecule down against the surface and the strength of the
lateral forces of attraction which are holding the hydrocarbon por-
tions of the molecules together. Both of these forces are important
and it has been shown that good boundary lubrication is only obtained
13
when the chain length of the hydrocarbon portion of the molecule is
greater than 14 carbon atoms.
he bulk modulus does not enter into the film thickness calcula-
tions, but it is believed to be a signific~nt factor which in-
fluences both temperature and pressure peaks in the filml.
In addition 'to the above, there are other properties which may have an
important influence on the results obtained. These are: ,
Fluid density
Gas solubility
Wetting effects and surface films
Thermal properties of the oil
Effect of 'polar compounds
bottcm-most layer will not move at all, the uppermost will move
with a velocity, U, and each intermediate layer will move with
a velocity directly proportional to its distance from the sta-
tionary bottom plane. This orderly type of movement in parallel
layers is known as streamline, laminar or.viscous flow. The
F required to impart motion to the layers
force per unit area (A)
is called the shear stress, while the movement of one layer of
oil relative to another is the shear strain rate. The rite of
14
shear (S) of a particular layer, sometimes called the velocity
gradient, is defined as the ratio of its velocity to its perpen-
dicular distance from the stationary surface, and is constant
U UI U2
=
for each layer: S =h h etc. Newton correctly deduced
F
Shear Stre;s A
= Rate of Shear = U (1)
h
FF U U
h2 h
b 1 h2
Fig. I
15
TABLE III
CONVERSION CONSTANTS
Multiply By To Obtain
Multiply By To Obtain
Multiply By To Obtain
Multiply By To Obtain
16
NEWTONIAN LIQUID
A B
17
PLASTIC THIXOTROPIC
-* >- PSEUDO-PLASTIC
1-4 1-4E-4
0
E U
0 W
0
L..
U3,
0
SHEAR STRESS
18 .,I-9674
treats the subject from a practical viewpoint. Murphy and Zisman
(Ref. I) discus- the subject as ft relates to structural chemis-
try.
19
of the attempts to understand the combined effects and correla-
ti:ns, and a br:ef presentation of some less well understood
rheol3g:cal considerations will be given.
20
where: V _ kinematic viscosity, Cs
k = constant between 0.6 and 0.8
T = absolu:te temperature in degrees Rankine
nC = are constants for the given fluid
The'Walther eqsation has been employed by the ASTM to gener-
ate Athvscosity-temperature charts, (Ref. 16) which are
21
large variation of viscosity with temperature and the other
is said to have a VI of 100 because it has a relatively
small change of viscosity with temperature. One fluid is
selected from each of the standard sets so that all three
fluids have the same viscosity at 210F and then the viscos-
ity index of the unknown fluid is defined as the following
v+A = BT (6)
22
I-.-= w4O
_ _
-3 .0
MILLIO C3
j0
uju
toJ '-j0
L ----- 00-
LL-_ R 0
00
0 0
>~
";om 1 "4
8. 0
If 0 04
0
OF
~mn-j
Ilz )ISN3 AISSA Lb3N)
1111
II ko 111IY 023
B/ T () !
SAe
V AT + BS (8) !
where particular
A, B, and C are constants which must be determined for
~each case.,.
24
2II
100,000-
10,000
I tOO I
425O
Hflv57
p~raffinic fluid molecules. In general, it is to be ex-
pected-that the increase of viscosity with pressure with-
in a given class of fluids will be greater for cyclic
molecules than for the linear molecules.
0 (9)
.l=- - (10)
Lp T 6p I
rn (L + C) '-(11)
.7.
N 26
pfessions should have t replaced with Zn-l)C. Fer all
fluids investigated to date, the plot of isobars on the
temperature-viscosity charts (Fig. 4) are straight lines.
27
influenced the results but the pressure-viscosity behavior
seemed to show significant correlation.
(c) Shear'Effects
The various types of shear dependence of viscosity have been
described in the earlier section dealing with the definition
of viscosity. The-pseudoplastic and thixotropic types of
non-Newtonian behavior are the most common types of non-
Newtonian behavior observed in fluids that are used as lub-
ricants. When the apparent viscosity decreases with incrip-
sing shear rate or shear stress, this is defined a'3 pseudo-
plastic behavior, and w1 ll be the subject' of this section.
Thixotropic behavior is the decrease of viscosicy with time
atf a constart shear rate. This will9be disc ssed in n e
next section. 1.
28
has the shear rate, the shear stress is also defined. On
the other hand, it can be argued that the shear stress is
the more logical variable to use because that is a measure
of the force applied to the fluid. The shear rate is a
*"measure of its response to that force.
Wright and Crouse (Ref. 32. have also studied polymer thick-
ened fluids and attempted to correlate the results using the
product of the shear rate a-d the shear stress as a paramet-
er. This is the energy inpir rate per unit volume. Philipp-
off (Ref. 33) has employed a reduced variable approach to
try to correlate non-Newtonian data.
29.
stresses due to the forces of gravity.
30
T
.. I~0,000,,
TYPE 37 BEARING
0 DIMETHYL SILICONE 1000 CS. AT 250C
X SILICONE BLEND 1000 CS. AT 250C
8000 --" R-2 SIZE BEARING
A DIMFTHYL SILICONE1000 CS. AT 25*C
0 SILICONE BLEND 1000CS. AT 25 C
% 0 CNTISTOKE
E--- DIMETHYL SILICONE
94000
2000 -- _-
O0 2 4 6 8 10 12
COASTDON TIME-MINUTES
r/
13
i ?
31 ~F1-97-
importance to lubrication is the thixotropic behavior asso-
ciated with the mechanical and thermal shear degradation
of hEgh molecular weight materials at the high shear flows
associated with lubrication. Elastohydrodynamic contacts
subject the fluid to very high shear rates and shear stress-
es, consequently dissipating large amounts of energy in the
fluid. This can caase the long chain molecules to be frag-
.4 Vented into smaller molecules, thus resulting in a-lower
molecular weight material which exhibits different viscous
behavior. As would be expected, the higher molecular weight
materials will be more isceptible to this shear degrada-
tior.
32
//
unit volume per unit time in the fluid. In elastohydrody-
namic-contacts, the energy input rate is very high'and,
even though the volume of fluid involved is small, this
degradation phenomena probably plays an important role in
elastohydrodynamics particularly when polymeric type materi-
als are employed.
33
in the case of silicones; however, there are others who
have argued to the contrary. For example, Forster (48) has
tried toexplain that viscoelasticity was responsible for
the increased load carrying capacity of spur gears which
was observed by Borsoff (49) with increasing speed in the
high speed range. Forster claimed that the loading times
[in the high speed gears was of the order which would result
in the viscoelastic behavior of the lubricant..
34
rr
whereas the latter wegt t& higher pressures (80,,000 psi)
but was limited to 10 dynes/cm shear stress (shear.rates
of only about 104 or 105 sec-lat high pressure). Qualita-
tively, nothing unusual occurs when the independent varia-
bles are c-ombined, but the quantitative prediction of the
behavior of fluids under these combined conditions is yet
to be satisfactorily accomplished.
+ C] log(l +-p-+
2000' logD (12)
Barlow and Lamb (27, 47, 52) and Phi.ippoff '33, have suc-
cessfully employed reduced variables (basically a shift i
frequency) to graphically correlate pressure, temperature
and shear rate.data from the oscillating crystal'type visco-
meter. The data used was limited to pressures below 15,000
psi. They have not attempted to put the data into an expres-
sion. -
35,
0r
molecular structureand free volume.
(3) uk Modulus
36
I.I
Table'IV (abstracted from Ref. 25) shows some typical values for
various types of petroleum and synthetic oils. The dimethyl
silicone oils have.the low-est bulk-moduli (68), i.e. they are
the most compresgible oils. - /
t-ABLE IV i !
* TYPICAL VALUES FOR ISOTHERMAL BULK
" ' MODULUS AT'25C AND q/,PSIG
37
() Lubricant Density
The densiLy ranges for some typical petroleum and synthetic oils
are given in Table V.
TABLE V
TYPICAL DENSITIES OF VARIOUS OILS
38
I
.1
Actually, only those oils which contain halogens (chlorine,
fluorine, bromine, and ioiine) ir the molecule are markedly
different in density. Like the silicones, these halogenated
oils have weak electrical forces holding the molecules toge-
ther, and the molecular spacing is far apart. However, the
halogen atoms replace hydrogen atoms in the molecale and their
atomic weight is sc much higher than hydrcgen that they have a
large influence on density.
TABLE VI
COEFFICIENTS OF !HERMkL EXPANSION
FOR VARIOUS TYPES OF OILS
Coefficient of Expansion,
0
Type of Oil Fractional Volume Change/
39
r.I
data obtained. Henry's Law, which states that the solubility
of a gas in a liquid is directly proportional to the pressure
of the gas above the liquid at a definite temperature, seems
to hold fnr most of the gas-oil systems although the solubility
does decrease at higher pressures, above 75 psia. The solubil-
ity of gases in oils generally increases with increasing tempera-
ture except for carbon dioxide.
TABLE VI!
SOLUBILITY OF VtiRIOUS GASES
IN PETROLEUM OIL OR DIESTERS
Helium 2.0
Hydrogen 5.1
Nitrogen 10.0
Air 15.0
Argon 24.5
Carbon Dioxide 94.0
40
I
the bulk of the liquid is much greater than the number in the
vapor above. Because of this inward attraction, the surface tries
to contract to the smallest possible area with the result that the
surface behaves as if it were in a state of tension. This force,
which is the Surface Tension, C, has a characteristic value for
each liquid. it is the same at every point and in all directions
along the surface of the fluid. It may be defined as the forces
in dynes acting at right angles to any line of one centimeter
length on the surface. The existence of a surface implies an in-
terface between two media, e.g., a liquid and air or a liquid and
solid, and the surface tension depends on the nature of both media.
CGL
0 can be any
gas, G Liquid, L value between
90 and 180 0
CGS z CLS
solid, S
where CGS' CLS and CGL are the surface tensions of the gas-solid,
liquid-solid and gas-liquid surfaces respectively. The contact
angle 0 depends on the three interfacial tensions, but whether it
is greater or less than 900 depends on the relative magnitude of
for and C If the gas-solid tension (CGS) is greater than that
for the liqid-solid interface (CL5), then the cos 9 must be
positive and 9 is less than 90 . if the reverse is Lrae, then 9
0
must be between 900 and 180 .
In the former case, e.g., water
on clean glass, the water is said to "wet" the glass. [This is
a matter of degree since complete wetting implies a contact angle
of 00]. In the latter case, e.g., mercury on glass, the liquid
is said not to "wet" the glass. These examples can be illus-
trated as follows:
/I/
Wetting, Water on Glass Non-Wetting, Mercury on Class
9 less than 900 9 greater than 900
41
The primary reason for introduzing this subject into a discussion
of elastolydrodynamic lubrication is because the wetting charac-
teristics of commericai lubricants and bearing surfaces are not
well-establi43 " or controlled. in both hydrodynamic and elasco-
hydrodynamic calculations, che assumption is made that the layer
of oil i-nediacely adjacent t: the bearing surface is firmly ad-
herent and immobile. It is also assumed that the properties of
-the lubricant at the solid-lhquid interface are typical of the
bulk properties of the fluid. If, for example, a steel bearing
surface is being .sed with a pure lubricant that wets the steel
readily, then these assumptions are reasonable. However, if the
steel surfaze is covered with somei form of a reaction film, or if
the oil'contains a surface-active compound which is preferentially
adsorbed at the .brsnt-metal. interface, then the wetting char-
acteristics of the b9Je-til snd the steel surface can be drastic-
ally changed. For conventions! hydrodynamic bearings, this
would have no signif:cant effect. The magnitude of the film
thicknesses frvolved would mask such thin film effects. However,
in elastonydrodyr.amic lubrication, particularly under marginal
conditions, the characteristics of the surface films might be im-
portant. If the fluid did not wer the beariog surfaces, or if
the physical properties of the fluid were significantly different
in the boundary layer (because of preferential adsorption or con-
centration of surface-active compounds on the surface), then thisI
could affect the film thickness calculations.
Two studies can be cited which show that wetting or surface film
formation could be important. In the first study (71) a number
of lubricated, instrument size, ball bearings, which appeared to
I
be in satisfactory condition before being put on test, showed sur-
face distress and erratic torque results after short periods of
running time. Examination disclosed that the bearings were con-
I
taminated with some unknown films and that the bearings were poor-
ly wetted by the lubricant. Once these bearings were cleaned so
that gocd wetting was obtained, no further difficulties were
experienced.
Another example of this type was observed in the NASA high tempera-
ture ball bearing tests described in Ref. 72. A polyphenyl
ether was found to be satisfactory when running in air at 600 F,
but failed at the same temperature when run under nitrogen.
Johnson (73) has hypothesized that these failures were due to the
formation of non-wetting surface films which prevented adequate
heat transfer through the oil. It has also been suggested that
some form of a contaminant film might inhibit reactions between
the lubricant and the bearing surface, thus affecting boundaL
lubrication. Regardless of the mechanisms involved, these
studies show that sirface films may play an important role in
elastohydrodynamic lubrication.
42
Under heavily loaded, high speed conditions, ir has been found
that bearings lubricated with silicone oils may show an excessive
temperat-re ri-se with no tendency for stabilization. This could
also be the res--lt of non-wetting or reactio&-inhibiting films.
143
TABLE VIII
2
where K is in BTU/hr.ft/ oFin.
p6 0 = density in grams/cc at 60 F
44
I
\p
ments such as chlorine, ,:ulfjr, phosphorous, etc.) are also used
as additives to prevent metal seizure. It is beyond Tie scope
of this discussion to go into any detail on the mechanisms of
lubrication in this regime. The subject has been covered very
well in many references, e.g. Refs. 3, 4, 6. and 7. However,
the fact that the th.icknesses of elastohydrodynamic films
are
often on the same order of magnitude as the surface roughness
values indicates that metal to metal contact probably occurs
frequently at least to some degree. The influence of these
sliding contacts, and subsequent surface reations, on elasto-
hydrodynamic lubrication has never been resolved. Fein (66 and
75) has pointed out that this coild be an important consideration.
Among the candidate oils which are being onsidered in this program,
the following lubricants are of particular interest-
1) Ester base oils meeting the }'L-L-7808 specification
2) Ester base oils meeting The MIL-L-23699 specification
3) Mineral oils
4) Polyphenyl ethers
5) The dimethyl, methylphenyl, and halogenated silicone oils
-5
45
2. FILM THICKNESS
Theories and experiments have shown that the fjim shape of a line con-
tact is reasonably flat throughout the conjunction, with a constricted
section at the exit. The minimum film thickness at the exit is ap-
proximately 70-75 percent of the nominal film in the flat region.
Since the nominal film thickness determines the level of the separa-
tion in the conjunction, the calculation of this quantity has been
recognized as the most important problem in elastohydrodynamic lub-
rication and has received considerable attention. In 1949, Grubin
and Vinogradova (82) first successfully obtained the solution of the
Reynolds equation at the inlet section by assuming the inlet deforma-
tion profile to be identical to the Hertzian profile. This approach
circumvents the necessity of solving the coupled elasticity nd Rey-
nolds equation and offers great simplicity in determining the separa-
tion' between the line contacts. Grubin's results can be summarized
in the following expression in dimensionless form:
8
h -- 11
00R 1.95 (ELI) 1 (15)
x l
where
E = E'
UoU
E'R'
x
F w
E'R'
x
46
I1
I
2 2
El 2 E+E
1 2
press-ire-viscosity expdnent
U = rolling velocity 1 - U )/2
U0 inlet vi'scosity
F = modulu& of elasticity
= Poisson's ratio
47
Lt-
11 E0. 60.7
_2 1.6 013 (16)
x
In the high speed region the only experimental data are those mea-
sured at Battelle Memorial Institute in the past five years (18-?3).
These data show a drastic reduction of film thickness when speed
increases, as indicated in Fig. 9. Furthermore, a strong load-
. dependance can be seen in Table IX and in the circumferential film
profiles shown in Fig. 10.
48
Icn
~ -
wn cn
wl L&JMl
Llg U) cn k
on LL
0 cr-
;2! 0
WQ C LJ~
0R 3
dzLLJO JJ
N*C D ) 0 <~ D"i~L
MommD
00 to C).Z :) Lc
If 1 -j - c C)- <X)
CMC)
L'jCILLj3
5- L0
oa-
49-X 0 OL-
REINO EE ORDINATE M, -
50 ABSCISSA
N, , R EC ! 0- 0 -
O-IZ,--2WSON
AND HIGGINSON 2
I A-4 ------
5 10 0
C-Os 051
50 00 504100
2 2-961
-4---- -
CD A N
CDE
l a-0
5 :l --
C:)
N: :
C:I
N r-N
0-0
1 LOAD __
+ F-.
I,.~ ~
-J-HJLl ~ c
- -L.
wL
za:
U Z > CD CD -1
0:: 0 '- vILn
a- wiQ LA I-. L
a. U
- ~a. a: 2
-~~ 0 -V
CD Q0 C Q U-
90099- w~ OM
LLCN -C: o
-r
00 CD < 0: 0nC
I r- C-1 C) -
.::
- -
z 11 C'),
a: a.
x <.~
LLL
51
-400
3M
Maximum calculated
_
{
Hertz pressure
80.000 psi
100,000 psi..
140,000 psi
U,
c
E 000-
LL
U;
C
xx
Direction Of 5.
rolling
-15 -10 -5 0 5 10 15
Circumferential Distance, -' inch A-46219
52 1471-9686
r
TABLE IX
MEASURED FILM THICKNESS BETWEEN ROLLING DISKS
LUBRICATED WITH A HIGHLY REFINED MINERAL OIL (21)
T 01-
-
l b-sec in/I PH~ l h o~l6
1
T
R 0.2in 2
in2/Ib psi in/see Qxh11
m h D
10
1 17.5 37.1 0.47
14 21 54.0 0.388
17.5 i5.5 51.0 0.30s
" 20 13 49.3 0.264
Qm KfT (17)
0 U 0.7
h DOW (OE') 0 )
R' 1.6 0.13 (18)
x w 7
h
(19)
ShDOW
53
Fig. 9 plots two groups of cT curves, the ratio of the actual film
thickness to that predicted by the isothermal theory, calculated from
the present theory. The solid curves are for a lubricant with a mod-
erate temperature viscosity exponent, while the dotted curves are for
a lubricant with strong viscosity temperature variance Conditions
used in the experiments are listed in Table IX. In the experimental
data, the value or is assumed to be the ratio of the measured film to
that calculated by the Dowson and lHigginson formula. These data are
also presented in two groups. The open dots are for a super refined
mineral oil and the black dots are for polyphenyl ether. It is seen
that the data for polyphenyl ether blends smoothly with those obtained
for the mineral oil. In general, the trend of the theoretical curves
appears to be quite in harmony with the experimental data. However,
quantitatively the measured film thicknesses still lie considerably
below the calculated value particularly for highly loaded cases. Dis-
crepancies as much as 2.5 to I still exist at higher loads. There
are several reasons to which this discrepancy may be attributed. First,
the fact that the discrepancy becomes greater for thinner film thick-
nesses suggests racher strongly that the measured film could be some-
what lower than the actual film in the contact due to any slight axial
misalignments between the two disks. Secondly, the assumption of a
flooded inlet in the theory may be unrealistic. A delayed pressure
generation at the inlet due to a restricted lubricant supply or due to
early formation of vortices in the inlet region would certainly cause
drastic reduction in the calculated film thickness.
To sum up, Lhe film separation between steel rolling and sliding con-
tacts can he accurately predicted by film thickness formulae proposed
by Grubin (82) or by Dowson and Higginson (91) for rolling speed up
to 250 in/sec. For high speed contacts, the reduction of film thick-
ness due to thermal effects may be determined by the data provided by
Cheng (93). However, there is still a disagreement between these data
and the BM measured data with regard to the effect of load. It is
uncertain whether the discrepancy is due to certain factors neglected
in the analysis or due to the inherent inaccuracy of the x-ray tech-
nique at these small film thicknesses.
54
r
b. Film Thickness and Film Shape of Elliptical Contacts
The film thickness at the inlet of a ball contact was first analyzed
by Archard and Cowking (107). They successfully obtained an analyti-
cal solution to the three-dimensional inlet problem by assuming a
Hertzian deformation profile. Their results show that the reduction
of film thickness due to side leakage is less than 10 percent. How-
ever, experiments conducted by Archard and Kirk (86) show that the
reduction of film thickness due to side leakage in a ball contact zone
is about 40 percent, considerably greater than that predicted by the
theory. From the experimental data they found that film thickness
follows a simple power law dependence on velocity, load, and material
parameters, similar to that of Gr-bin's.
-i V
0
O = tan
U
U ~U 1+2 U2
55
Elastohydrodynamic Contactb
4
(a) (a)
(b)
(b) (b)
56
rI -9692
B CONSTANT FILM
THICKNESS LINES
A A
DIRECTION
OF ROLLING
MINIMUM
B FILM
hmin
0 hmin
57
1-957 5
.17
VI 4. V2
V
2
R
x
R'
0
o
inlet viscosi:y
CY pressure-jiscoiitv exponent
3. PRESSURE r)lSTRIBUTION
58
HIGH A
MEDiUM A
HIGH A
59
has a strong influence upon the pressure peak. increasing the pressure-vis-
cosity exponent has the same influence as increase A, it gives rise to a high-
er peak and moves the peaK toward the entrance region.
In 1961 Archard, Gair and Hirst (92) employWd the same inverse proced-
ure and .found extensive solutions for the heavily londed cases. They
separatedthe contact zone into four regions, two in the inlet and
two in the outlet region. For each of these regions, a different itera-
tive procedure was employed in order to bring convergence. The
60
Lxep
After Hertz
.0.6
-0.4 -0.5
0
'~~~ ~
-(C-I-L-
61 %Trl-9687
_.-G.4 -02 0.0
,o
0.5
-1 6
62 MTI-9705 u.
5
10
9 4
V -, 7
0
O 6
2 2
0
-2 -0 1
X, pounds
63
.1T1-9697
111
P KlO
-4-3-2-I0 X, pounds I
U "- [0] 0 foi cow---] (1 10', [2] 101 [3] I0"r [4] 10', [5] 1-0! [5] I0- [6] I0"*
W-3xIO's G=5000
64
MTI-9696
agreement between the pressure distributions by Dowson and Higginson
(91) and by Archard, Gair and Hirst (92) is good except for the in-
tensity of the pressure spike.
d. Pressure Measurement
Subsequent attempts were made by Orcitt (i18) and Kannel (119) in mea-
suring the pressure'distribution in heavily loaded contacts. The
pressure sensing in a narrow contact region was made possible by us-
ing a vapor deposited manganin strip having approximately .00]" in
width and a few microLnches in thickness. Orcitt's experiment was
conducted with a glass disc which limits the maximum contact pressure
to-about 50,000 psi. At this pressure, be showed that the pressure
distribution deviate considerably from Herizan and indicates some
65
09 -I
f~ ~ ____U7926 K SEC
(25% SLIP) -
07 -(ISOTHERMAL
,til 05
-- 15 - o5 0 0 ,0- 14
-9571
66
MT1I-95 71 1
I
_U_7926 N /S
0_ t _ *1 1623 355~iii
07- . I 1160 1
05 15 LB/IK f
04 -SLIP RATIO 025
03 '
67
MrI-9570
23- PERCENTAGE SLIP
22-3
21-
20-
17-
15-
68
MrI-9615
marked disturbances at the exit region. However, no evidence of
pressure spikes was detected.
4. TEMPERATURE DISIRIBUTION
69
2
160 1-PRESSURE-103 Ib/ifl
140.,
I MEASURED PRESSURE PROIES
P, -126 000 Iblifi2
12 P 1600l/n
P 1 1 104 000 Ib/in2
REGION 40
20-
30 20 10 0 10 30
CIRCU.MFERENTIAL DISTANCE-103 in
T = 115TF.
RS = 1400 ft/min.
Lubricant: polyphenyl-ether.
(Courtesy of! Proc. Instni. Mech. Engrg., Vol. 180, pt. 3B,
1965-66).
70
MTI-9613
T 15 F
Rligspe~ed 4,500 fpii 0
T - Lubricant: polyphenyl cther
-20~~ 0 . 00 2
LOO&n TmsOf
Itt ROII'ng~eed
---- 2900fpm
71
4k
CL0
rc
00
4-1
uu
E 0
%1 4
C.44
44 W
to 00
01
wW-
724
- Rolling speed =4,500 fpm
Lubricant: polyphenyl ether
NO.Ompo
Twv~okaOD
115F
.0 -0 0 D2
T 115 F
Rolling speed =2,900 fpm
LubiCNFIRtf
-2DM~ -9007002
or 4
00o "4
0'O-
-
cc
---
0 $u4
* C.)
%C a
00 0JlS~
isd0 0 M
14-
054
-4
-44d
N 0~
0W
00u
-4 -0 W
V-4cc.
-474
250
-Highly refined mineral oil
Polyphenyl ether
200.
150..
ISOI
I')I
0-
25,A sec
- - -
0k
- -" I 'I I
Time-
75
MTI-9694
250
Highly refined mineral oil
-- Polyphenyl ether
200
in 150
0
x
.100
50
r/
12.5A eC
Time-
76
TII'-9693
!
T tion of the lubricant. A brief summary will be given to various theories
available in determining the temperature in the lubricant film.
a. Approximate Theories
The assumption of a mean viscosity across the film was later relaxed
by Cheng (93), and he developed a rigorous two dimensional finite
difference method in determining the temperature distributions across
as well as along the film. Typical three dimensional temperature
surface within the-film for a pair of steel rollers are shown in
Fig. 32. In this figure
AT
9 T (22)
0
77
C!'
a)~f
CD,
D0
+ CI
E OF 6
78 NT-964
*1 20
100
ATMEAN
W-715
LWIN.
70
60
40-
30 - 2L
AAT
0- ~
q".5 0 05 10 14
79
HTI-9569
0
-a')f
10 0
800
AT = temperature rise, OR
: x =exit coordinate
the Hertzian whichin.is approximately equal to half of
width,
h film thickness
s = (U2 - I)/U2
The temperature plot is divided into two parts at the pressure peak
in order to show the clear details before and after the peak. Since
this is a high speed case, the surface temperature rise at y = 0 is
practically insignificant comparing to the mid-film temperature rise.
The effects of rolling and. sliding speed on the surface and mid-film
temperature rise are shown in Figs. 33 and 34. For a constant slip
ratio, the mid-film temperature always increases with the rolling
speed on the surface temperature would depend upon the slip ratio.
At high slip ratio, the surface temperature decreases with the roll-
ing speed.
c. Temperature Measurements
81
250
2002
0.1
10
50-
W=715 LB/IN.
82
MT1-1529)
*40
0-
1-0.05
0-
K 50 100 200
UIN./SEC
500 1000
83
Krl- 1531
/ 2
100-20 \10020
01818
80-16 80 16
4 Nx0 g x10"
Hx106 10
70- 14 70 14
!60-? 6012
408
50- 10
' i 40- 8
30 030 6
0403
0J "~ 24681012
0 3-, C)35 03~ 3
- -1 0
-2 -1 0 0
-2 x (b) U- 1o-- 05x1O--
- -s
(a) U 0-IO V=-O'2xIO 1 W=3"19xIO
110 22
100-20
90 18s
016 ..
1 4
Nxl I10 PERCENTAGE
SLIP
70 14 (x. xOO)
14-
TEMPERATURE
60 12 60 12RISE
(C) 2 0
50 10 10-
40 -e
5
30 6 6-
2 115 20
-3 0 5 1 34567 8 2 1
00
-2 -' 1 x
x 0 1-2 -1 x . 1
1
(c) U- 10- 1 V-035x10 W-3.19X1O (d) U-10-11 W-'=-3"19xO-8
84
?11-9612
ture profiles for a pure rolling case are shown in Fig. 36 in conjunc-
tion with measured pressure and film profile. The variation of mea-
sured surface temperature against the slip ratio (U2 -Ul)/U 2 is shown
in this figure and, they are compared with the theoretical value cal-
culated by the thermal elastohydrodynamic theory. It is seen that
the agreemdnt is fair at high slip ratio and rather unsatisfactory
at low sliding speeds.
Recently, Kannel et al (102) have made a series of surface temperature
measurements under pure rolling conditions and for loads higher than
those used by Orcutt. A titanium wire deposited on a silica insula-
tor was used as the sensing element. Figs. 37 and 38 show the effects
of load and speed'on the surface temperature. In general, the temp-
erature level is considerably higher than those measured by Orcutt.
5. FRICTION
Early measurement of the friction between gear teeth and between cir-
cular discs have placed emphasis mainly upon finding empirical corre-
85
40
*20
U1 -250 in/sec.
I U=251 in/sec.
10 U 150 in/sec.
U = 10 in/sec.
0 01020.3 04
U2 -U1
Slip Ratio,
86 kfrl-970b
60
Temperatures corrected
-Measured Temperatures
E240 Max Hlertz 04: 00 4.0
S pressures (psi) 128,.120
20
0. 20 40 60 s0 100 120
Circumferential Distance, 10' inch -46
T =L15F
amb
Rolling Speed =1400 fpm
Lubricant: polyphenyl ether
Fig. 37 Variation in Disk Temperature with Loading
(United States Air Force Tech. Report No.
ASD-TDR-61-643, pt. VI, 1966).
Measured temperatures
90fpm
Rolling speed l4Ofpm
40- 20e4p60r0atures
T b '=115 F
Lubricant: polyphenyl ether
Fig. 38 Variation in Disk Temperature with Rolling Speed
(United States Air Force Tech. Report No. ASD-TDR-
61-643, pt. VI, 1966).
87
H? 1-9701
lations with the design parameters. Typical examples of contributions
are by Cameron and Newman (123), Hisharin (124), and Benedict and
Kelley (125). But none of these papers were designed to promote the
basic understandings of the traction between elastohydrodynamic con-
tacts.
Crook (126) used two kinds of rolling disc machines in measuring the
friction in the line contact as a function of sliding speed. In the
small sliding speed region he used a four disc machine which consists
of a center disc surrounded by three equally spaced outer discs (Fig.
39a). Since the center disc is free-floating, the measured torque
will not contain any extraneous torque from the supporting bearings.
For this reason the four-disc machine gives very accurate friction
torque measurement at very small sliding speeds. However, in the
region of high slip, the four-disc machine is not suitable because
it cannot maintain a stable sliding speed. For high sliding speeds,
Crook used the conventional two-disc machine as shown in Fig. 40. In
order to investigate how the friction decreases with the sliding speed
in the high slip region, the rotation of both discs must be controlled
by variable speed motors, one of which acts as a motor and the other
acts as a generator. In this manner, Crook was able to measure the
friction characteristics throughout the entire region of sliding speed
using the four-disc machine in the low slip region and the two-disc
machine in the high speed region.
Using the four disc machine, Crook found a profound influence of roll-
ing speed upon the frictional torque in the low slip region. In this
region the slope of the traction vs. slip curve is equal to a so-
called effective viscosity divided by the oil film thickness. By mea-
suring the slope of the traction curve and by calculating the oil film
thickness from the existing elastohydrodynamic theory the effective
viscosity can be readily evaluated. If the thermal effects and the
non-Newtonian effects of the lubricant are both absent in this region,
then the effective viscosity will stay constnat as a function of roll-
ing speed. However, this condition was not found in Crook's experi-
mental results. On the contrary, he found a marked influence of the
rolling speed on the effective viscosity of the lubricant. As shown
in Fig. 41, the ratio of the effective viscosity to the supply viscos-
ity can be reduced by almost two orders of magnitude as the rolling
speed increases from 0 to 1200 cm/sec. Such a drastic reduction does
not appear to be due to thermal effects only. Crook speculated that
it was the visco-elastic effect of the lubricant which prevented it
from reaching the static viscosity in the short time interval as it
passes through the contact zone. The visco-elastic effect has been
more thoroughly investigated analytically by Dyson in a recent paper
(128).
By using the two-disc machine, Crook was able to extend the friction
data in the high slip region for loads ranging from 7.5 to 20 x 107
dynes/cm and rolling speed from 400 cm/sec to 1200 cm/sec. All the
friction curves showed the same basic trend which is characterized
by an ascending portion at a small sliding speed and a descending
88
F
gerCan
Lad
Fig.
39 Th For DriskMchn,
*89
.9TI-9699
Fig. 40 The Two Disk Machine
(Courtesy of Philosophical
Trans., Series A255, 281,
1961).
90
MTI-9695
I
friction at high sliding speeds. Fig. 42 shows that an increase in
load does not change the basic characteristics of the friction curve
much, but does increase the level of the friction force. Similarly,
in Fig. 43 it was shown that a change in rolling speed also effects
the level of the friction force. When the rolling speed increases,
the friction level decreases.
Using this simplified thermal theory, Crook was able to calculate the
coefficient of friction or the iffective viscosity as a function of
sliding speed. Fig. 44 shows the comparison between the calculated
effective viscosity and the measured effective viscosity as a function
of sliding speed for a load of 1.2 x 108 dynes/cm and a rolling speed
of 400 cm/sec. There are three calculated curves representing three
different values of the ratio of the thermal conductivity of the lub-
ricant to the temperature viscosity coefficient. It is 5een that even
if this ratio is increased by almost one order of magnitude, it is
still not possible to predict the sharp reduction of the effective
viscosity at small sliding speeds. Based on this evidence, Crook
concluded that the friction force atosmall sliding speed cannot be
accurately predicted by considering the thermal effects only.
91
10~
io2 b
0 600 1200
U cms
-~ 0
40,.031 r-4
* 4
4 0
0
C0
V.4
"MO0.019 0
0 200' 400
(U2-U1 ) (cms)
141
5.2 10
4-4
00
L)
0200 400
211
10000
F 5000 %*
%
I %%
%% %%b%%
%.
C%.
%
5'00
0 10 20 30
211
U2 36.2 IN./SEC
0.01-
94
NrI-97 17
iI
drastic reductions of friction can exist if the lubricant viscosity
is shear rate dependent, according to Ree Eyring. However, in all
T their calculated data the friction force was found to be depeident
upon 1/h which means if the tolling speed increases the friction
force will fall according to 1/h. On the other hand, all the experi-
mental data gathered thus far have shown the proportionality is far
greater than 1/h and in most instances are more nearly proportional
to 1/h2 . Thus, one may conclude that the inclusion of the Ree Eyring
model in the friction analysis will not be able to predict a reduc-
"tion of friction force at low rolling speeds.. This is really not
surprising, since the Ree Eyri*ng model is primarily a shear dependent
Newtonian modeland is not likely to predict the loss of pressure-
viscosity dependance during short transients.
95
05
44
z
19
0
01
1/4 1 4- 16 .64 256
SLIDING SPECU
ii
-
0.08 1
Y1\p,243x 05Ibfiin 2
0.7-%
0.0
I176k
of Inst. ,
Meh r. Vol 182 30,1676)
0-04 2 2596
0 151 5 10
SLDN PEWId-r/
U,=cntn 7i/;O t 0C
_ 0-07 -
I
0.06
- 0.05
0~U~4V
Aor-
30 50 9
9O* U. ins
0"02 0o 0 7"8 "
S- 16
L AL A 34
x x 47
S001 + 76
V 153
** 260
-0 10 20 .30 40 50 60
1 4 _ _ _ __ _ _ _
011-2
97
KTi-9619
simplified thermal theory. These curves are for rolling speeds of
16 in/sec and 260 in/sec and are based upon apparent viscosities de-
duced from the low-slip friction data. It is quite evident that the
calculated traction coefficients can be as much as one or more than
one order of magnitude higher than the experimental friction data.
Discrepancies become larger when the maximum Hertzian pressure becomes
greater.
Dowson and Holmes (132) modified the Crook's four-disc machine and
investigated the effect of surface quality upon traction character-
istics of rolling and sliding contacts. By plotting friction against
the surface roughness, they showed (Fig. 55) that the friction
93
O04 -
~II,- *
-J - ii U16M/5
-noJo0o2" -.-.--x- -- - -
01
- .. eT-, - ,, ,
00 , 10 Is 20 25 30
SLIING SPEED,
WU-U -i/s
.... P~z-- 70001b/
Po - 243 000 lbf/'.
--
T- 30C, K -6-,0 x 10- ca /,n s deSC.
- -- ,176 000bffm
T .3WCK.6.0xIlO-,I.,sdeXC.
0
3
X10
1007 ib41 1
7I
9c-
StIEARPLANETEMPEIRATURE,
99
[-70
008-
A
006 r+
Rog
0 0
;7 003- Smooth
I discsS
0-02 ci-
0011- ~ ~
0 -- h __ -___ U
A C
r
Velocty Pofil
35
Slip 0200
25
_ _ 450I
___
0 20 -
256 Ibf/m
xx
0,11 iIbfl ++
101
9618
.TI
.4
initially decreases with surface roughness, reaches a minimum, and
then increases steadily with surface roughness. Unlike Johnson and
Jefferies' conclusion on the effect of surface roughness on friction,
Downson and Holmes found that the influence of surface quality is
quite pronounced. However, these two results may not necessarily
in direct contradiction because the load used by Dowson and Holmes was
much smaller than thatused in Johnson and Jefferies' experiments.
In the non-linear region, the strain becomes. large and the use of the
viscoelastic properties as determined in oscillatory shear with small
strains becomes questionable. Dyson was able to establish a relation
between the shear stre.ss in continuous shear and the properties
measured with oscillating crystals under small strains. This re-
lation was obtained for'the Barlow-Lanib jnodel as wellas for 4 the
Maxwell model. He showed that the characteristics of the ascending
friction curve can be qualitatively explained better with the Barlow-
Lamb model than with the Maxwell model.
In the thermal region, close agreement was found by Dyson between the
measured mean shear stress-and the limiting shear stress based on
the viscoelastic properties from the oscillatory shear. However,
he hadto postulAte that the'shear modulus must be dependent upon
the contact pressure and mid-film temperature.
102
! /
LINEAR I NON-
REGION ILINEAR'
HTHERMAL REGION
A-A
103 MI -% 0
true viscosity
________ ____11_ % Nsm-2 0t
E_____ 0___ 9_____ zero rolling speed
N' 08a 3
6__0_7___ 4 -_10
006
ei4 mentoT +
"6n-a6
~0 2
0.4.
00 104
0.0-979
true viscosity
,N-2Lat
0
* 7 zero rolling speed
_ _
T1_
-. 1.0 __ _ _ _ __0_ _ _ _ _ _ _ _ _ _ _
E 069
00
0 7 xperim_
0 .
02
0.3
0-1-
0- 00
V. 004
105
* ~ MTI-9719
generation due to sliding only, the calculated tractions for the
following four viscosity models:
Model 1 - Newtonian
uexp T (23)
0
where ', , and y' are constants. Fig. 59 shows the friction
calculate4,by their simplified theory agress very well with
those obtained from the full elastohydrodynamic theory (12).
Model 2 - Newtonian
tI +exp (T-T)
(Yp) (24)
5
0
By adjusting the values for y' and 5: in this model, they showed
(Fig. 60), that the calculated friction could be made to agree
very closely to those measured by Crook (127).
where S is- the shear rate, G is the sheer modulus, and r is the
shear stress. By using this model they showed, Fig. 61, that a
value of the shear modulus G of 5000 psi is required to produce
significantly different traction than that calculated by the
Newtonian model.
With this model they found the effect of shear modulus is con-
siderably less than that produced by the Maxwell model.
b. Spinning-Friction
106
I.[
S ~ N
I;%
S 0
It
iU 7 15 S.
SW oI
UUz
i! ,00
F 198
i10
I.,
g.lO~/ / -
-4
.23 p
U-t57 is/sec
Z0 model Eq (93
,---- crook do%
20 40 0 ' 120 40 o6
108
MTI-9691
40
46
20
40 60 SO 1-0 20 84 S
109
Mn-9690
Driver
Dr-%e ~ D n~e Deftpne
-- Uper bearing
- - ar n - - S ha ft
Wan.ifolc -
-Dead-%%eight
/load assembly'
ar
Supprt
.e~be .. ~'Collar-/ -Fulcrum
Shaft-/ Weight
tip- specimen
Support
Lower test-
housing
assembly,
L--j CO-7933
110 HrI-9709
frictional torque varies markedly with contact groove geometry, con-
tact stress, lubricant type, and speed.
1. Exponential Law
I' = (10)
2. Power Law
e p (27a)
t = Io l forpfp
[
eeI= P* + &2 (p - P*)] for p >p* (27b)
When two elastic bodies are in contact with each other inder a heavy load,
the stresses near the contact region can often reach unusually high values due
to large localized pressures in the contact. Materials subjected to these
.1 repeated high stresses will gradually develop a crack and eventually chip away
from the contacting body. This phenomenon is commonly known as surface pitting
which is one of the major failure modes in rolling element bearings. In order
to predict the conditions under which pitting occurs, it is necessary first
to make accurate calculations of the stress field in the contact. This sec-
tion deals with the analysis and calculation of these stresses in the contact-
ing bodies.
The first published classical work dealing with stresses in elastic contacts
is by H. Hertz of Germany in 1881 (138). He obtained the pressure distribution
[ 1ii
.06- - rfoREtiCAL PERCENT
S---E.'rZRNTAL COWOPMIn'Y
IN.-LB-60
I
112
Mri-9688 J
in the contact for two elastic bodies having counterformal surface. When the
surface geometries can be approximated by two paraboloids, the pressure dis-
tribution in the contact zone is found to be an ellipsoid. He also solved
for the stresses in the body in terms of a potential function. Since Hertz,
numerous works have appeared dealing with analytical solutions of stress in
the contact. A detailed account of these solutions can be found in a paper by
Smith and Liu (139) and in a review by Lubkin (140).
However, these theories are all based on the assumption that the normal or
tangential loads in the contact are symmetrical with respect to the center of
the contact. This assumption is valid for a dry contact, but it is invalid
for a lubricated contact. Recent elastohydrodynamic theories (Refs 93 and 141)
discussed'earlier have shown that the pressure distribution including the
effects of a lubricant becomes greatly deviated from Hertzian - particularly
in the exit of the contact where drastic changes of the pressure profile are
found. Horeover high local frictional forces are shown to exist if a small
sliding velocity is introduced in the contact. Consequently, there is a
need of a method to calculate stresses in a solid considering an arbitrarily
distributed normal and tangential load. We first review the calculation of
contact stresses without traction.
113
[\
P
Y?
zj
a.-
xI
b z
Xb z
Fig.
6PonCoc n otatElpe
114
n?-92
cause the surface of the disks to be deformed elastically
over a region surrounding the initial point of contact,
f thereby bringing the two bodies into contact over a small
area"as-show-h'in Fig. 64b.
For two bodies in, line contact, the solution of the contact stress
problem is relatively simple. The area of contact when a load of w
i .lb per unit length is applied is a long narrow rectangle of width 2b
in the x direction and length t in the y direction. The stresses at
ponts on the z axis below the surface are given by:
2
ii )I_1 + (z/b) -1zb) b
ax - I + z b)
21 (28)
0
y= -2v[ 1+( )- (.]) 29
az = - F (z/b)
I+ 2
bb
I(30) A
where
b : 2w/Tr (31)
I- 2
lI-V2 R'
22
A =4, + "'R ,1
SI ,)- E2 (32)
'xl " 2 x22
I0
are given by
115
1/2
HZ [ b]HZ
IT:I
Fig. 65 Line Contact Geometry.
116 HTI-9659
r
(b) Maximum Shear Stress
01 = -0.1856 A (34)
02 = -0.7861 (35)
03 = 0.9718 vi (36)
and
T = 0.3- (37)
max A
B /R + I/R
B 1 1 I + iRx 22 > 50 (39)
2
yl 'y
then the maximum stresses can be calculated as for the line con-
tact previously discussed. This means that when the above in-
117
MAGNITUDE OF THE STRESSES z
ON AND BELOW SURFACE
x _ X
Irz
118 PfrI-9650
equality is satisfied, the area of contact between the rolling
element and race is very nearly a long narrow rectangle. For
many cases, this situation holds.
3
b n V 4(-A+ -E (41)
L -1 1
B +- i___ I/in
X, xl
(42)
- 1 1
A = - + 1/in
2Ryl 2Ry 2
119
+0.30
40.25 --
0.20 -20
0.15
-0.05 - ,
. 0-o,
\ /
-0.15 -
-0.20-- - -
-0.25
-0.3 ---
-2.5 -1.5 -0.5 05 1.5 2.5
x/b
.5) E/
io.2 "~maxQ
Cm.
120 MTI-9667
the z axis.
This conclusion is borne out by the work of Greenert (146) who ob-
tained excellent correlation of the orthogonal shear stress range
with the fatigue life of AISI 52100 rollers.
If the contacting angle, 9 (see Fig. 65) is less than 5 degrees, the
stress and deformation in the solid can be accurately calculated by
neglecting the curvature of the boundary. Using the elasticity
solution for a unit load on the boundary of an elastic half-space
one obtains:
121
0.8
0.7-
0.6 \
E b/b
b/o
122
MTI-9624
F
= 6 xyt (x,y,g) = -2 i -.)2 ]
Saxn (x,y, TT R4
-3
~
8 ay n 6a(x'y, =
_ --~ 2 7--.4
R (43)
6 2-
(X'y,) =
-- 14--
4
R
where
8a n, 6a__nrn
1 = dimensionless normal and shear stress
components due to a unit normal load at
x=
1: x = X/
"" = y/x*
= /x
123
z (JJF)
Lzll
KF
p(
124
24rI-9645
0.30
-o #=30 #=60 4zD
0.10
X
0.00 -
[I w/
H -0.10
-0.20
-01.301
-3.00 -2-OO -1.00 0.00 1.00 2-00
DISTANCE - X
125 T-
' ,.
yp() yn + () 6ay t dc
= g xyn + 4yt4)
where
/ (E)= 2 -U x )(U
x Q) x bh( C)
To illustrate the use of the program, a sample case was run wdith no
tangential load and a normal-loadir- given by
w -- 2
p = I- (45)
7. SURFACE TOPOGRAPHY
126
r
begin to touch at the tips of the higher asperities. The study of the aver-
age distance through which the one surface moves into tha other surface and
the total real area of contact resulting from thds movement provides import-
ant information about the friction and wear which occurs in the contact zone,
IThe theories of friction and wear are based on the fact that real
surfaces are rough so that the area of contact is limited to dis-
crete spots on the two surfaces. It has long been agreed that the
most convenient way of eYxplaining Amonton's laws -- that the fric-
tion is proportional to the load and independent of the apparent
area -- is to take the real area of contact to be proportional to
the load.
To find the total load and the total area, we substitute the excess
height(z-h ) for C and sum over all load-bearing asperities which
gives,
127
P f(z-ho) NO(z)dz A g(z-h) NO/(z)dz (47)
Thus, when we know (a) the distribution of asperity heights; (b) the
deformation law for a single asperity; we can imzediately derive the
total load and the total area in terms of the separation and from it
the load-area relation, itself.
Greenwood and Williamson (151) showed that the real area o' contact
for the elastic deformation was proportional to the load. Greenwood
(149) extended the analysis to show that this was a consequence of
the statistics of the surface and not the particular mode of defor- -.
mation.
C
128
mathematics. It is important to the realism of the theoretical model
identify the characteristics of typical surfaces and to determine
which of these can accurately be represented by a Gaussian distribu-
tion. Williamson et al (152) have divided surface preparation pro-
cesses into two classes, cumulative and extreme value.
Extreme value processes are ones for,which the final surface re-
flects only t1'e extreme of the events which occurred at each
point. Extreme value processes nearly always involve the rempv-
al of material such as in grinding or peening. Turning, shaping,
and milling are processes.id which the events are neither cumu-
lative nor random and create non-Gaussian surfaces.,'Grinding
on the other hand, can yield a height distribution which is
approximately Gathsian. In grinding, the grit is located ,/ith
a Gaussian height distribution, and creates a topography which
is the result of many interacting Gaussian probabilities. Its
heightdistribution appears as a gentle curve on probability
papers. Fig.-73 which is a distribution presented by Williamson
et al (144) supports this argument. For most zork there is usu-
ally little error in treating ground surfaces as Gaussian. This
is fortinate for lubrication engineers, as most rolling element
t(bearing surfaces are 'finished by grinding.
(3) Pure Topographies and Transitional Topographies
129
0 0
I-3
UJU
130
II
1 99.9-00
I I- L / :
195
f80. I
! So-
U,
50-
~2S
2,-
1 1
T-
I I I I I I I I
I
! .o
I 131
1(11-9617
et al (152) termed a transitional topography.
Ao A 1 (2 2
-2.n/ 2 2)
c I
~h o
Ac = A 0(h o) (50)
// /
132
i z
II 00
a=so
50-
I A
133
.fTI-9616
mitted through film lubricated contacts between rolling elements
and races with negligible slip. Film thickness predictions confirm
the presence of a hydrodynamic contact mechanism in rolling element
bearings. No detectable cage or rolling element slip can therefore
be expected in highly loaded bearings. The actual motion of a roll-
ing element in a given bearing will be located between the limiting
cases computed for rigid and elastic components.
Harris (154) using the Dowson and Higginson (153) formulation cal-
culated the cage speed in a high speed cylindrical roller bearing..
As illustrated in Fig. 75, he was able to predict quite reasonably
the cage speed versus load for a high speed roller bearing. How-
ever, discrepancies between theory and experiment occurred at low
loads corresponding to the intermediate lubrication regime of rigid
bearing surfaces lubricated by an oil having pressure dependent
fluid properties.
Boness (155) showed experimentally that oil supply cam have a signif-
icant effect on cage and roller motion. Using the Dowson and
Higginson model (153), Boness calculated the'normal and tractive
forces at the contact for restricted oil supply to the lubricated
contact by varying the location of the point in the oil film where
the pressure commences to build up. The theoretical results imply
that reducing the oil supply to a minimum value required to main-
tain full hydrodynamic conditions reduces the cage slip by up to 75
percent of its fully flooded value.
134
430-0
000T
aOPm
0
0 . . .. . .
- o 0 o i~o~ 0 20
[LI
35 0
.'TI-9 0
0 o
w 0
-.
LL- - DO
4
0
0
41
4J
40
P4
136
will roll at the outer contact and spin and roll at the
inner race. The ball motion about its own axes and about
the bearing axes is said to be controlled by the raceway
at which rolling occurs, the outer raceway in this case.
Fig. 77 shows the forces and moments acting on the ball
under the influence of raceway control. Eq. 51 from Ref.
216 establishes the criterion for raceway control.
PI al (a2 C1 ) > P2 a2 e-
2 (51)
137
I:I
OUTER RACEWAY CONTROL X=0
INNER RACEWAY CONTROL X I
I
I
I
I.P2 % I
2 (-FiLs w
II
I
I
:1
13B xrl-gs I
L0
IY
0- NORMAL FORCE
F-FRICTION FORCE
-M MOMENT
XYZ DIRECTION IN CONTACT
-
PC ZONE-
FC- g- GYROSCOPIC
C- CAGE
FV- VISCOUS DRAG
x Z11 7
.1 FV FY2
a1
CF
II-. - C.
FC
139 MTI-9649
052
3 BALLS 1.125 DIAMETER, 9000 RPM
o-SHEVCHENKO AND BOLAN [157] TEST DATA
0.50 ---- RACEWAY CONTROL THEORY
- HARRIS ANALYSIS [156
0
.0.48-
S0.46-
0.0
0.44
Cl)
-0.44 - "--
0.42
140 HTI-961,8
this discrepancy
141
/
a In a perfect geometry bearing under pure thrust
load there are no force patterns tending to cause
equilibrium in a few revolutions. The meagerness
of the lubricant supply complicates the situation
because the balls impact the separator causing
sbarp increases in force and high accelerations.
More lubricant tends to soften the impact.
142
TABLE X
FAILURE MODES OF ROLLIMl-ELEMENT BFARINGS
I. CONTACT FAILURES
1. Wear Failures
a. Mild Wear
b. Severe Wear
2. F railures
a. Cold Flow
-- b. Overheat Softening
1. Bulk Failures
a. Overload Cracking
b. Overheat Cracking
c. Bulk Fatigue
d. Fretting of Fit Surface
e. Permanent Dimensional Changes
2. Cage Failures
3. Lubricant Failures
143
I
various faults and defects. Fatigue failures of rolling-element bear-
ings have been systematically classified by Wren and Moyer (166) with
illustrative examples of each failure mode. The Russian, Edigaryan
(167) claims that the life of rolling-element bearings is dependent on
only three basic failure mechanisms. A recent classification of roll-
ing-element failure modes is that by Tallian (168). Other good reviews
are given by Cheng and Orcutt (169), Smalley et al (170), O'Connor and
Boyd (171), and Moyar and Morrow (172).
* Degree of Lubrication
In bearings which are\eitherunlubricated or poorly lub-
ricated such as with jet, engine fuel or cryogenic fluids,
wear is the predominate failure mode. Wear in these
fluids ie so great that a sufficient number of cycles
never accumulates to form the familiar subsurface spald
of contact fatigue.
144
e Amount of Lubricant Contaminatio1
Amount of Sliding
A
[1
' a rolling element and its track. These slip mechanisms arise
from applied loading and the geometric conformity of the system.
Creep measurements arising from such slip are presented and com-
pared with theoretical predictions. The wear resulting from
this slip was determined using the radioactive tracer measure-
ments.
145
N, satisfactory explanation of this velocity dependence was
fFund.
V =-kPI (52)
where T
P= Load
I = Distance slid
Ac = k' P (53)
146
Substituting Eq. (53) into Eq. (52) one obtains the wear rate
k
V/t=(-) A, (54)
This equation states that'the wear rate is, for a given set of
conditions, lubricant and surface geometry, proportional to the
total area of asperity contacts existing simultaneously.
A /A F(h )
where
I A
c
=
*
area occupied by film thickness less than h
0
A
-k' total Hertzian contact area
147
/'!
amounts to removal of a layer of approximately 400 microinches
thickness from the rolling track during the time the balls
travel one million inches on each other. This is quite discern-
ible wear for a rolling contact and will, in long run deetroy
the accuracy of the configuration. Thus, a contact operating
under the very high Hertz pressure pmax = 680 kpsi cannot run
.long with the thin elastohydrodynamic film corresponding to
= . For higher values of go 3, the wear rate is about
l0 times lower and, at these levels, it is negligible as a
failure mode.
The wear law cited states that the wear constant k is (as a
first approximation) independent of the elastohydrodynamic con-
dition, and characterizes wear properties (boundary lubrication
properties) of the lubricant/material combination of the rolling 2
Bamberger et. al. (80) have evaluated the wear preventing action
in rolling contact of three advanced lubricants, a synthetic
paraffinic oil, a fluorocarbon, and a 5P4E polyphenyl ether.
Their results show that if partial elastohydrodynamic conditions
prevailed, bearing failure was mainly from wear. They demon-
strated that, by reducing bearing load for the paraffinic and
fluorocarbon oils to insure a full elastohydrodynamic film,
failure was primarily due to fatigue. With the 5P4E polyphenyl
ether, wear and surface distress occurred, indicating partial
elastohydrodynamic lubrication for all operating conditions.
Their results validate the use of the elastohydrodynamic film
parameter as criterion for evaluating the effectiveness of an
elastohydrodynamic film.
148
0.5 inches long. The larger roller was ground with a 50 cham-
fer on each end to limit the contact length, L, to 0.180
inches -nd reduce end effects. Most of Landen's data was taken
at a load.w, of 13,000 lb/in of contact length and at this
load the contact width, 2b, was 0.056 inches.
Both Landen (179) and Queener (180) have shown Cia: the total
wear can be expressed as the sum of the constant wear ml = kPI
and the transient wear m = mo (I-E-nt). The total wear is
then
-
m = kP +m 2 (l-e ) (57)
149
?1
F TABLE XI
.6x10- 6 0
0
0 0
E
.4
o 0 OContinuous
u 0 Colored Surface
0 0 Film (5)(6)
0
00
I .2
~- ,Not Rippled
(1)(2)
l 0_
0 0.2 0.4 0.5 0.8 i
150
Krl-681.2
materials and slip ratios of interest in gas turbine bearing
practice. Determination of these coefficients should form a
part of any continuing program in the bearing failure area.
151
resulting from elastohydrodynamic action will tend to reduce the
chance of smearing.
152
/L
TABLE XII
SURFACE DISTRESS DATA OF BODENSIEK (184)
-1.
Fr
153
many investigators (187-196) in the well accepted critical tempera-
ture concept, which states that smearing occurs when a critical
value of peak surface temperature (temperature at inlet plus peak
transient temperature in-contact) is reached, regardless of the
hydrodynamic film thickness. The broad concept of a critical temp-
erature i's consistent with modern boundary, lubrication theory which
says that good lubricants form adsorbed or themically'reacted films
on the surfaces. These films protect the underlying surfaces as
they slide over one another so.that friction is low and surface
damage is negligible; however, as the temperature of the surfaces is
increased, a point is reached beyond which the films melt, are dis-
solved by the lubricant, or are desorbed from the surfaces. When
this happens, friction goes up, still more heat is generated, and
the surfaces eventually weld and tear with catastrophic results if
the sliding speeds are high enough.
In the previous discussion, the specific film thickness and the maxi-
mum surface temperature have been identified as the major parameters
governing smearing. A graphical representation of thp region of
smearing or severe wear is shown in Fig. k. This plot uses the in-
verse of the maximum surface temperature as the abscissa with the
specific film thickness as the ordinate.
154
.
Smearing possi-
ble in the pre- No wear, no smearing
sence of dirt
and foreigh
particles.
0 I/T*
00-
(Critical Temperature)
155
____
fi
A principal barrier to wide-spread practical application of elasto-
hydrodynamic principles, such as the film thickness formula, is the
fact that the surface temperature at the inlet regions must be
known with soire accuracy in order to establish the inlet viscosity.
Blok (193), as early as 1937, offered formulae for the flash tempera-
ture rise occurring in a loded contact which, if added to the inlet
surface temperature oi blank temperature, gives the maximum surface
temperature in the contact. In spite of the more recent observa-
tions by Fein (192) to the effect that a transition temperature
criterion is oversimplified, it has gzeat usefulness in judging the
likelihood of severe wear failure. Archard (194) has sumnarized
flash temperature calculations for easy practical use. He shows
that an upper bound of the attainable flash temperature in a con-
tact can be calculated by assuming clat metal to metal contact
occurs over a single contiguous area and by assuming that plastic
flow takes place over the entire extent of this area.
termined.
156
conditions of operation will be determined primarily by boundary
lubrication effects.
Leach and Kelley (188) have shown that there is a critical tempera-
*" ture which defines the transition between effective boundary lubri-
cation and severe wear (smearing). This transition point is the
temperature at which polar constituents in the oil are either
melted or de-adsorbed from the contacting surfaces and effective
boundary lubrication no longer exists. If an elastohydrodynamic
film separates the surfaces, then effective lubrication can persist
above this temperature transition point; however, any metal to
metal contacts through the film will result in surface damage.
157
Surface Lubricating Films
Each of these films has its own unique properties and thus, differ
in the manner in which it provides surface protection. It would
be expected that each has its own unique method of failure. In
general, the failures of these films may be categorized as follows: ""
When sliding commences between two surfaces, they are usually mis-
aligned and their surface asperities are interlocked. Any films
present on the surface at this time will be broken away. However,
as sliding proceeds, the surface temperature increases and the sur-
faces "run in". Surface films are then more easily formed and less -
likely to be broken away. If the rate of formation is fast enough,
effective sliding will result. The rate of formation of films will
depend upon the type of film formed. Generally speaking, some sort
of a reaction is necessary so the rate of formation is strongly
depending upon the temperature and the concentration of reactants
in the contact area.
Destruction of Films
158
Tlum
i K. (58)
n api m are constants, n < 1 depending upon the load m < 1 depend-
ing upon the velocity and the geometry of the contact. These high
temperatures can cause the destruction of the surface filmd for a
variety of reasons:
Melting
Dissolving
Decomposition
Evaporation
The destruction of the film allows metal contact and the possibility
of seizure.
Deformation
Material Seizure
Once the fluid film thickness approaches and the protective lubri-
cant films have been broken through, a seizure process can result,
under certain conditions. These conditions are a function of the
material properties, the friction or adhesion, and the geometry of
the contact.
159
Fig. 82 The Two-Dimensional Plasticity Model
2 +9 2 = 9-2 (59)
m
If r4 is high the process will continue and necking will occur away
from the interface. This will ultimately lead to fracture as shown
by the curves dotted lines of Fig. 82. That this process occurs
has been experimentally verified by the work of Greenwood (202),
Cocks (203) and Antler (204). A precise failure criterion has not
been developed; however, it has been suggested that it is associa-
ted with the decrease in area in the necking region, the work
hardening whicu strengthens the interface region, or the ductility
of the material. Whatever the explicit criterion, it is apparent
that it is associated with the ability of the material to deform.
This also conforms to common experience - that brittle materials
are muc better in sliding than a ductile one.
160
I
the more resistant the material must be to deformation to prevent
surface damage.
It has been shown by Johnson (205) that the load at which a metal
first yields in rolling contact is given by
a = 4-r (61)
max s
161
This represents an increase of 66 percent in load carrying capa-
city above the load for first yield given by Equation (60). Tan-
gential forces introduced by sliding or microslip cause an appre-
ciable reduction in the shakedown limit as shown in Fig. 83.
162
I.
Continuously
1.5 Plastic
FF
0F
0
z
0.5
163 MTI-965 1
Sha t-down Limit rI,
-4)6
NC
rl0 10-
-.
C
0-5 0 4 10 3 24
164 -o
fatigue failure, and that the stress field is the result of a
Hertzian normal pressure distribution. The constant proportion-
ality and exponents need to bo determined experimentally, and,
therefore the theory provides a failure criterion rather than
absolute values. As such it enjoys popularity and is the basis
of conventional life calculations.
165
tude of, the maximum shear and octahedral stress and shift them
nearer to the surface. These stress alterations increase the
likelihood of surface fatigue. The degree of worsening depends
on the operating conditions and the pressure-viscosity exponent
of the lubricant. Increase in speed or viscosity in the full
film elastohydrodynamic regime tends to increase the magnitude
of the secondary pressure peak and hence the maximum stresses.
The effect of surface traction is most severe when oils of a
high pressure-viscosity exponent are used under small tr moder-.
ate sliding to roll ratios.
166
16
fatigue in this regime. Pioneering work, in-this area has been
contributed by Christensen (211, 212). He found that che maxi-
mum pressure between two spheres approaching each other under
7a
| constant load can exceed the Hertzian prebsure if the speed
and the pressure-viscosity exponent are large.
167
44
0 g.-,
t -, w
-Jma
- - -
(in)
U..j
(l) CL 0
Go V
j 0
oxo
ui w
LUx
-4
La.. 00
I 168
plicit in this discussion, however, is the very significant
additional fact that the influence of film thickness cannot be
considered apart from other variables. For example, the load
effect occurs primarily through its influence on lubricant film
thickness in the contact zone, and it follows that other varia-
bles that also influence that film thickness will shift the
failure mode. Thus, increases in speed and in lubricant viscos-
ity, decreases in temperature, and decreases in surface rough-
ness can be expected to shift the transition zones from one
failure mode to another.
169
/I
SECTION III
ELASTOHYDRODYNAMIC PERFORMANCE CALOULATION
The design tool which has been the major goal of this study consists of a
computer program which calculates all major variables in an elastohydrody-
namic contact. The calculations are based either on the most reliable
theory of elastohydrodynamic lubrication or on the empirical formula'deriv-
ed from the most recent experimental data. As outlined in the previous
section, the input to the computer program consists of:
With these input data the program calculates the following quantities success-
ively by means of a series of subroutines:
The use of subroutines provides a convenient way to replace any old methods
of calculation with any new techniques which may be developed in the future
for any one of these variables. Using this approach, one does not have to
rewrite the main computer program as the state-of-the-art of elastohydrody-
namic lubrization advances.
170
/
/
/'
The computer program is primarily intended for use on rolling element bear-
ing contact problems. However, it has been made general enough so that it
is applicable to many other heavily loaded concentrated contact problems.
Thus, it should also find use in the design of gears, cams, and pivots.
1. FLOW DIAGRAM
2. ANALYSIS
The analysis for the variables listed in the flow diagram are described
in the following sections.
The semi-major and -minor axes of the contact ellipse for a given
total load P, geometry and material properties are calculated from
Equations 40 and 4 1 on page 119.
w = 3P (63)
4a
w
E'R' (4
(64)
x
U
=
0
E'R' (65)
x
where
171
Read and Write Input Data
ED
D , Qm' Cr' '
m
B', ', D-
EDR
, E
J
Subroutine FTFS
Calculate:
I. the thermal reduction -factor I0
T
ir
Calcul te:
1. the isothermal minimum film
thickness (hmin) s
2. the minimum film thickness hin
Subroutine EBHM
Calculate the protrusion width e
172 1/ HI-9675
I
I
Subroutine ABAR
Calculate:
,1. the specific film thickness h /a
o n
2. the real''area of contac, A
Subroutine FRCTN
Calculate the coefficient of friction f
Subroutine PRESS
I Generate. the Pressure Distribution P(x)
1
Subroutine TEMPT
-Calculate:
Subroutine STRESS
Calculate the subsurface stress distribution
Figure 86 (Continued)
173 z--9T7
a 00 inlet viscosity lb-sec/in 2
PIHz 3 P
E' 2 yabE' for the elliptical contact (67)
and
PIIZ I -"
" R for the line contact (68)
x
c. I.ubricant Paramezers
C, a a rE'
C!
art = 2" P Z
1* ~.puz(69)
' - 5/T
0
P
2 T 0{
d. Thermal Parameters
2
40 (UI + U 2)
QM " 2KfT
; cfUa3
Qc 0 2
fx (70)
K2 1/2 RE /
TICK 4wE3 ,'
i 174 .
J e. Thermal Reduction Factor
T
0
"" PH
H- dimensionless Hertz pressure
E-U U2 - U 1
I
s slip ratio U
U2
In Ref. 105, the curves for c>t indicate that the most influential
parameters is the viscous heating parameter Qm" The other four para-
E4
and s all play secondary roles in governing the
thermal reduction factor. In particular, the influence of the slip
ratio is vpry mild. A review of these curves shows that within the
range of the parameters investigated, the thermal reduction factor
can be represented by
PUZ
T =
fl(l - 0.1 s)(l + f3 -j) (72)
175
II
TABLE XIII
VARIATIONS OF f WITH R9SPECT TO
,~, and Q
5 1
T. 105
* 1 0,B/' _
o2
_ -T .
0.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0
0.01 .99 .98 '.98 .99 .97 .97 .98 .98 .98 .
0.02 .98 .96 .96 .98 .95 .94 .96 .95 .94
0.05 .96 .94 .93 .95 .92 .90 .93 .90 .89 -
0.2 .89 .86 .81 .86 .81 .76 .84 .78 .74
10.0 .36 .27 .22 .34 .22 .17 .34 .21 .15
176
VARIATION OF WITH RESPECT TO
t103
~
15.71~ ~57 ~~~~n
157 35 114m14 14
M.
f. Side Leakage Reduction Factor
Rx -( C -- -n
'
(73)
where C, nl, and n 2 are numerical functions of a/b. These values are
shown in Table XV.
TABLE XV
, I=a/b C nI n2
- 5 '- C
S- [ -- U ] - (nl)'=5 n2 - (n2)Y=5 (73)
In calculating the minimum film thickness for the isothermal case, the
formulas proposed by Dowson and Higginson (153) based on their numeri-
cal solution is used.
178
multiplying this by the thermal and side leakage reduction factors,
one obtains the actual film tffickness
h. =_ (h minisso (7,5)tM
h. Protrusion Width
0.4
- ' 02[2 oU 0 0-:76 (6
2411
2411 )O" E'R1J (76]
I X]
where K is a proportional constant and
I
I =-( (- 3)d (77)
0
e x -(b -e
e)
1.0
179
TABLE XVI " (
ANALYSIS OF e BASED ON CHENG'S ISOTHERMAL THEORY (REF. 93)
0o w h
s = 0.0 E'R'x E'R'x aE' e H=h
e h
xo
0.0 11.75 x 10 - 5
2.828 x 10 5650 .890 .875
180
TABLE XVII
ANALYSIS OF e BASED
ON CHENG'S THERMAL THEORY
(L F. 93)
U2 -U 1
---- - :u
- = u0U
oq_w-
U2 E'R' 0R'E'R'
evE'
EtR e =e/b H =- eE
x x
e h0
0.005 2.87 x 10-11 2.828 x 10 - 5 4473
0.005 .49 .854
6.01 x 10 -11 2.828 x 10- 5 4473 .65
0.005 .859
11.662 x 10 -11 2.828 x 10 - 5 44?3 .81?
0.005 .859
18.00 x 10-11 2.828 x 10- 5 4473 1.03 .849
0.020 2.715 x 10- 2.828 x 10- 4473 .465 .854
0.020 5.360 x 10 11 -
2.828 x 10 -5
4473 .65
0.020 .854
9.91 x 10 - I1 2.828 x 10 - 5 4473 .83
0.020 .81-9
16.38 x 10- 11 2.828 x 10- 4473 1.04 .838
0.050 2.50 x 10 11 -
2.828 x 10 -5
4473 .49
0.05C .838
4.87 -c 1011 2.828 x 105 4473 .65
0.050 .838
8.34 x i0 -11 2.828 x 10 - 5 4473 .835
0.050 .838
14.55 x 10- 11 2.828 x 10- 5 4473 1.02 .832
0.100 2.238 x 10- 1 1 2.828 x 10- 5 4473
0.100 1.434 x 10- 1 1 .50 .827
2.828 x 10- 5 4473
0.100 7.611 x 10-11 .667 .838
2.828 x 10 -5 4473 .85 .832
0.100 12.61 x 10-11 2.828 x 10 - 5 4473 1.04 .827
0.250 1.76 x 10 1 -I
2.828 x 10 -5
4473 .50 .795
0.250 3.62 x i0 - 11
0.250 2.828 x 10- 5 4473
5.85 x 10- 1 1 2.828 x 10 - 5 .667 .805
4473 .833
0.250 -I .805
5.52 x 10 1 2.828 x 10- 5 4473 1.00 .805
181
It is noted that if the normalized protrusion shape ItI()does not
vary much with the operating conditions, the I, will not be a func-
0
t:on of OE', E. , and W. Consequently, the above exponential re-
x
lations should give the correct variation of the e with the speed,
load -Ipressure-viscosity parameter.
e = C ~
a ) . (78)
TABLE XVIII
VARIATION OF e AND C WITH SLIP
U2 - U
I
e C'
U2
These values are used in the computer program to predict the protru-
sion width.
182
0
-4
x
C4 cn
C-4
C
Ul
I I CNC
$4 0
In~ Cn JJ
ClC
s-I C183
41
-4
jG)j
0
4.'
01
CJ184
JM.4
C\ J E
' 4
-44
cu)
18i
for various slip ratios. It is seen that .1. is close to 0.8, which
is the value used in the computer program regardless of the operat-
ing parameters.
Area of Contact
Whenever the ratio of the nominal film thickness to the effective
root mean square value of the surface roughness of an elastohydro-
dynamic contact is smaller than 4, there will be occasional surface
contact. The extent of the contact area is dependent upon the ratio
of the elastohydrodynamic film thicknezz to surface roughness, tne
shape of the elastohydrodynamic film distribution, and the surface
topography. Tallian, et a] (177) have developed a combined graphi-
cal-numerical approach to predict the extent of the surface contact
for surface topography either determined experimentally or generated
by a random process. Their method gives the relation between the
contact duration, that is the fraction of time during which there is
at least on asperity in contact, and the specific film thickness.
While this method has an advantage of being able to accommodate any
arbitrary surface topography, it lacks simplicity and is somewhat
inconvenient for extended computations. There is a need for a simple
analysis to study the extent of surface contacts in a partial elasto-
hydrodynamic contact. The analysis described here is aimed at achiev-
ing this objective.
This simple analysis is based upon the following two major assumptions:
186
It
99.99
,99.9
99.8 2
98
98 - h
0
he
0 20 h 0.
* ~o 0
U
* 5.0
$4
h he
2.0 0, -==0.
1.0 eh.
0.5
0.2 h
1.e
0.1 h
0
0.05 e =
h
0
[Fig.
CF
[ ~187 1-9
upon rt-e area of contact. These curves also indicate that the pro-
tLusion depth has a stronger effect than the protrusion widthA
J. Fr~ctional Coefficient
8 0-s
lo_ (79)
2 = Kf
G3 o PHZ
188
ri~~~_____________-
144)
I I
o~~~~~ _-*__ -_
-C. __ ~~ f------
-J~E 0- -4
U).
-I--- 0 CU
0 PI
.4--4
IAI
I~ II i
"~\
~I'l4_
" I
' ___ j (0 '4
~ .~Io
a).-
ci 6~ 'i C; 0
NOID18 00I311130
~~189
Iin
cc'
, oo
It 4)
ItI N
_ _ _ _o__ _ _ I I 1 _ __' _ _ _ _ )
4o
, ; \... S
I0 440c
IL u
1 "1 Ni.'
-'-i------ i Io o o4 )
I l
it_1_ 0( -A J-
4N_1_1
0 0
_____ I I. '
\-
cS\-A 'NP
0 _; C 0~
0 001NOI13 30 IN31IUJ330 -I
190
01
0 )
-4J
C3
I0
i~iI i -4
U 41
0A4
- 0c 4J
UU)W
0 0
_li __ (o_ U,
fcII
0 1I0
r' 0
I z4
1 l. Il----- r-- 0
W1311
.. N11A
191-
[
__ __ 1 _
_ _2_ _
wt
I ______M
Icc
Win
> r-
00
4. 4.1)
4f7
520N
L- 41
-- Colt ~4
% 0 0.
(L . - % 0
11 Ln11
a j
"a0 "a'
V\ Ic%
_____ ____ 0C
6~~ 0 0 0;
M010183 .40 INi3IJ3303 -1
192
II?
mI
I'IT
oc
INN
'o >.
W4J
0 0-
( ~ I
s-
t4
CCD
\
q.
WCD
19
-IT
n 5.1
\ \aL
L'' rL
II V 1
il I I.-
% 4 r=
*" , 44
, \%
1141
05 0
" S S 0
N 2.
\S ,
% %
- 0 Q C0 0 i
NOILO183J 30 IN3131JA303 - i
194
ICI
to
zw
10 c
Cu *0
'\
"
', ", 1 , (n
>
%0
____ % ' - % ) a
",,"f" ~" 0o
0'iIt 0 .0
"" ' x 4- 40
\\ 1 ",,
'S4 _ .PI
'. C. 0.J..-4
%'%S0
'
9 (V ' ","" 1-,-1
\ % 41
%%
0 0 c; C) 0 0 C) 'S S
W113I43
-4 IN331 IA03 -
195'
film.
k. Pressure Distribution
Almost all elastohydrodynamic theories thus far have shown that the
pressure profile for heavily loaded contacts does not deviate much
from the Hertzian elliptical pressure profile. The deviation occurs
mainly in the inlet and the exit regions. Due to the elastohydrody-
namic action'there is a gradual build-up in pressure ahead of the
Hertzian region which is absent in the case of a dry contact. In
the contact zone the pressure profile is essentially Hertzian. In
the exit region, the elastohydrodynamic theories show that there
exists a pressure spike right at the point the protrusion begins.
After the pressure spike, the pressure decays extremely rapidly to
the ambient pressare. This has been previously illustrated in Fig. 13.
*The pressure spike occurs exactly at the point where the pro-
trusion begins.
OPrcssure increase before the spike and the pressure decay after
the spike are assumed to be parabolic.
196
integration of the pressure profile between the center line of
the contact and the exit is equal to half of the load.
1. Temperature Distribution.
The thermal analysis for this method is based upon the assump-
tion that the viscosity across the film is constant and corres-
ponds to the mean temperature across the film. Using this
assumption, the temperature across the film can be integrated
t o analytically. With this analytical integration, the finite
difference approximation is only necessary along the contact
region but not across the bearing surfaces. This permits one
to solve for the surface temperature and the mid-film tempera-
ture distribution along tne contact region, without the calcu-
lation of the detailed temperature across the lubricant film.
197
f:im is removed by employing a finite d:fferernce approximation
of the temperature derivatives across the film. With the
friLte difference formulation, the viscosity variation across
-he film can be incorporated. This program has been usel ,n
solving the thermal elastohydrodynamic problem in Ref. 93,
whi:ch also gives a detailed accouit of the numerical analysis.
198
I
f is critically dependent upon the value of pressure-viscosity coeffi-
cient O at high pressure and temperatures defined in Equation (10)
on page 26. It has been shown by Crook (126-and Johnson ,(129)
that if the value of a is independent of pressure, both the calcu-
lated friction and temperature will-,be too high. This discrepancy
becomes increasingly worse as the contacting pressure increases.
Recently, Allen, Zaretsky, and Townseno (137) have shown that by
using a redced pressure-viscosity coefficient at higher pressures,
good agreement between the calculated and measured friction can be
otLained. This composite pressure-viscosity model is adopted in
this preliminary design program 'o calculate temperature pr-files
in the contact.
for p < p
Ln ) = VS + Ka' 1 (p - p*) + [ + iV(p - p*) , (81)
To
0 0
for p > p
n. Stress Distributions
199
Il
SECTION IV
APPLICATIONS
Three examples are given in this section to illustrate the use of the
Elastohydrodynamic Performance Program. Typical bearings from high,
moderate, and low speed applications have been selected to demonstrate
features of the program and its various input-output options. Detailed
input and output instructions along with the program listing are given
in Appendix VII. First, however, some supplementary input instructions
specifically f6r a roller-race or a ball-race contact are listed below.
The input should be the linear velocities of the two contacting bodies
at the point of contact. In bearing applications,it is often nece-
ssary to perform a coordinate transformation so that the center of the
ball/roller is fixed in the transformed coordinates in order to compute
200
[
the linear velocities of the ball/roller and of the contacting race.
Al This is illustrated in Fig. 99.
In Fig. 99, Q. and C are the angular velocities of the outer and
inner races with respect to the bearing center. o R is the orbital
velocity of the ball/roller with respect to the bearing center and
is the angular velocity of the >all/roller about its own center.
in a reference frame in which the ball/roller center, 0, is station-
ary the angular velocities of the ball/roller and of the races are
shown in.Fig. 99b.
In general, these constants are not available from the lubricant ven-
dor. In order to determine these constants from the standard viscos-
ity data which are generally available from the vendor, the pressure-
temperature-viscosity data in the 1953 ASME pressure-viscosity report
(23) have been analyzed by a least square fit type of regression pro-
gram based on the above viscosity model.
Using the results from this regression analysis, the following proce-
dure can be used to determine o, 0, and y from the standard lubricant
data given by the vendor which consists of
201
OUTER RACE
BALL / ROLLER R R 0
INNER RACE
ro ri ""'
OR
202 14 -9664
I I'
CL9 0 - -- - -
03b
('4J
_ _ _ _ ____
0 __ _ __
0 $4
00
0*
- - - - - - - - - _ _ __ __ 0
00
00
91 NI'01I
* - - - - - - _ . - - - - - _ _ __203
r
polate between these two curves.
(2) Calculate v by
60
100 + -Y (85)
0
(4) Calculate 1 by
For the first example we consider a deep groove, split inner race ball
bearing running at 15, 100 rpm and lubricated with MIL 7808 oil at
220 F. The specific bearing dimensions are summarized in Table XIX. -.
TABLE XIX
SAMPLE BEARING NUMBER 1
DEEP GROOVE BALL BEARING
With the above geometrical data, the necessary speed input information
for the ball outer race contact can be calculated by the method des-
cribed on page 200. The results are summarized below
204
I Orbital angular velocity of ball w_= 8400 rpm. 2
Speed of ball at contact (Body No. 1) Ul = wo " R 60
I= 1995 in/sec.
Speed of outer race at contact (Body No. 2) U2 (0 - w R)
/ ro0 6R0
= 2099 in/sec.
The individual ball loads were determined from the bearing internal
load distribution for a thrust load of 1560 pounds. The computer
input for the elastohydrodynamic performance program is shown in
Fig. 101,
The output sheet is presented in Fig. 102. Some important results are:
x/b = 0.
*Maximum shear st ess in the ball at (y/b = 0.01) = 0.08003 x
205
100
-~rz
-44
C0 0
-J
4D 4 Ma
OX 0
m N 'As *m
-4
(4) 4 D PJ,-I
000 00
'D C, -4 a xn
.0 o0 0 .0
206
PERFORMANCE OF AN CLASTOHYIRODYNAMIC CONTACT
GEOMETRIAL DATA
MATERIAL DATA
OPERATING CONDITIONS
LUBRICANT DATA
207
y1-9579
00
Ina
ct- z
* 0
o co 0
(UA1
- 4p.
ujS
0z
LAJ w
00
U LL.
is. 0
uj 3a.
n 0
495
W b- ul - IA -
j A
0 0 0
S
w B S S 0X
J i h i hi h i h
00aa a208 a
Ii 1111a 0
A 6
m 04%0 q O
I -x
Nx
o a.
6- CL %
4 4V4
4 15. Z IA 0l -
a m - 1
U) m 0 - A..
) 4 IZ 0
Iii a 40 M4
I a~ -z -0 -0
4 C 4 6... 3 Ia
u U.I
tr, z
IA D
3~ ~~ 3 0 0 0 M %
Z~ *n * T II wUI
4 153
W 15. W w *JjJ 4 1.3 n.J2095.
FwICTIONAL COFF. 3.52367E-02
ABW NUOMLL) FILM PwIESSU6RE
-4.ODOOOFeJU 9*~t)2ObE#01 ,.Do
-3.5000JF#0U 6.vd734E*01 4,42364E-05
-3.00000F*00 4.A335A#01 i,32W6E-O4
-2.50000E400 3.0s,64A*0.1 419O0iE-O4
-Z.OOUOOF>OO i.hat'tjE*U1 1.OO374L-0J3
-j.755J9F#90 1.09O7bLO1l kob19b1EU03
-i.5107tiF#00 b.4514dE*OO 3*9847$IE>O3
-1.38847F#00 '.hZ7O2E#OL b.39044E-03
-1.26b11EOO0 3.()971/t#00 1.16232E-02
-1.2O0,1EOO0 2.4t)311E*0O 1.66364E-02
-!.1438bE#00 1.498'1LOO0 2.54753E-02
-j.O8me7lEUU 1.4451l4LOO0 4,Z918EU02
-1.02155E#00 1.1174nL*OO o.3OI9SEO02
-1.010W#00O 1.O7',1lt+OO '.6337EO02
210
Mrl-9581
[I
ITRDUz 1 OU= 6.00000E-01" ERFRx 9.58697E-02
211
ARRAY OF Y-COORD. FOk S1RE,.S CALCu.
2
1.U000E-04 .00uOo0-0/
Yu 1.000OOE-02
212
HTI-9588
i
1
I XST SA SY PSI PSy
sAY PS2 PAy Os
-.61153 -. 6*dd -.. 0105 -.01330 -.46051 -. 5054? -. 32166 .02245 .07822
SHEAR STRESS AT E~e .3).00000 DEGEES
I -. 01330 l .01
uJu -- 6 -
A
1ST Sx - SI SAY PSI PSY PSz PIY OS
.00000 -. 62-363 -. b345)d .0219,. -. 60773 -.6b53 -.41967 .02240 .10089
SHEAR STRESS AT EVERY 4..uOOOO DE6NEES
-- -. 02195 .U044b .(i!1'
213
XST Sx SY SAY PSX PSY PSZ PTXY OS
75140 -.45101 -.'2ie -.Olt7b -.41449 -.45780 -.29047 o02166 #07091
SHEAR STRESS AT EvFY 4b.U0000 DEGkEES
-.01575 -.01486 .0115
XST SA SY
TAY PSX PSY PSZ PIXY OS
.84462 -.46406 -.486d3 -.00639 -.46131 -.48959 -o31665 oOI414 *07574
SHEAR STRESS AT EVERY 45.00000 DEGREES
-.00839 .01138 .00839
102 (continued)
214
MTI-9586
n
3l The thermal reduction factor is significant, .81, which means that
the actual film thickness is 19 percent less than isothermal theory
V
TABLE XX
SAMPLE BEARING NUMBER 2
.ROLLER BEARING
215
r Ok
z
-
X N
w
Aj 49 w co
In0 0 0:
Uta.
00
41
~ ~~.
Aj~ JwU 41u ij
w i - w wtd
Lw
F- 0N0 I 0 q
zn:- 0.PCl
- ~ ~ .-. 4E
N -4
216
J PENFO.RMANCE OF AN CLASIOh4YDHOOYNAMIC CONTACT
GEOMETRIAL DATA
II MAfERIAL DA
-~ OPERATING CONDITIONS
217
MI- 'S2
rr
m~~- c 14 z4
L) 0
ASS IS
V 41
'A 0
'A'
00
w u.
41.3>
IA N x~- -
w V) M W M a a0
i *A x , w w z ' '
14 or3 x1 CJ 4 w4
o a218
a a
II0 . .
0l
00 0 0 A
to) isu 0 80 .5
uj
3 *A
Ai 0 0 4
j ui
8 A ~ 0 i
0l 'A
'0 -t C'
4 z
-
-D
T -i
0 - U
I- ZJ w
go 1.-& -
3 WNa 3
219
0.94815
*Maximum Temperature Rise in Lubricant = 0.0253456 X (220 + 460)
= 17.3 0 F, occurring at
x/b = 0.94815
OMaximum Normal Compressive Stress in the Roller at (y/b = 0.01)
= 0.6712 x y x 2.333 x 105 = 2.46 x 105 psi, occurring at x/b
=0
iT
*Maximum Shear tress in the Roller at (y/z = 0.01) = 0.08629 x
x 2.313 x 10 = 3.16 x 104 psi, occurring at x/b = 0.94815
N = c9a (88)
'-
.. TABLE XXI
SAMPLE BEARIN
LOW SPEED ANGULAR CONTACT BEARING
The input and output data for the sampel calculation are shown in Figs.
i5C and 106. Several conclusions can be drawn from this data, First,
internally generated thermal effects are unimportant since the thermal
and isothermal film thicknesses are almost identical. Additional per-
tinent information is summarized below.
[*MaximumHertz Pressure
0 Minimum Film Thickness
= 47,222 psi
= 9.94372 x 10" 0 in
*Percentage of Area Contact = 4.6%
* Specific Film - 1.75782
*Friction in Rolling Direction = 0
*Since the Normalized Protrusion Width 5 = 0.553517, the
Pressure Calculations Breaks Down. A Message was Printed Out.
Due to the low loading, the calculated fatigue life of, this bearing is
extremely long. Because it is operating close to the partial elasto-
hydrodynamic region, mild wear is a possible failure mode. To inves-
tigate this failure mode, the elastohydrodynamic program was repeated
for different loads and oil temperatures. The results are summarized
.in Table XXII.
221
Z A
x
m 10
44
4t4
* 0 4.1 .0 4
0
4 .n
~- 0
I~ -Z r
L)% 4
* A* L 'A.:. J W a
1 10
22
wi. PER4FORMANCE Of AN ELASTOHYOI4ODYNANIC CONTACT
GEOr4ETRIAL i)ATA
MATERIAL DATA
OPER~ATING CONDITIONS
223
MTfI-9589
w i& S
:0 I-I
0 0o
-J
0 0
0 (n
JI.-
II
U I
Ai ILi I.-.
0330 5 3
A-
0 4 I 0
* -
* r
u- cr D -
0P - L .i
4i 3 4 c *.
0, CL) 0- V
U S224
0 a 0
m000 400
0I I
*- in i A j n a
(VO a aD 'x O
* 0 0 3c 0 0
*- ~ . O x'b-
*i -CD-
w : .
N 0L
4t
4 9 j L 0 l
x 0
m 4
Ln 0
z~ + '.) L)< -
-A) . A 0 N _ * 0 rl r4
D Aj
4Z L
w i
L:
IL X)
Li tLi2.
-c * 0 0L
I.- Z 4
w 0 0 (
w w - w 0. 4 L LU I-C
4 fl 225
TABLE XXII
SAMPLE BEARING NUMBER 3
LOW SPEED ANGULAR CONTACT BEARING FILM
THICKNESS AND AREA RATIO AS A FUNCTION OF OIL TEMPERATURE
h 0 sin) ToA A
This table shows that at the upper end of the temperature range, the
film thickness is small compared to the surface roughness, resulting
in a high degree of asperity contact. At the low end of the tempera- .
ture range, the film thickness is large, implying that wear will not
occur. Such a range of alternatives suggests a great need for care-
ful temperature control if this bearing is to operate successfully.
The table suggests one other important point, the surface finish is a
major parameter and must be carefully measured and controlled. A
change from a 4 jin to a 2 gin finish would probably reduce the wear
rate, initially at least by about 2/3.
226
formance of a ball bearing, the frictional forces between individual ball
and race must be calculated, assuming ball motion about all axes
The commonly used methods (215, 216) assume Coulomb friction at contact.
A recent paper (137) employs elastohydrodynamic lubrication theory for
line contacts to calculate the frictional torque on a spinning ball in
a non-conforming groove. In that paper Crook's simple friction analysis
(120) was used to calculate friction but frierional stresses in thedir-
ection of the major axis of the contact ellipse, outside a circle of
radius of the semi-minor axis of the ellipse, were neglected, and the
analytical results were made to match the measurements by bending appro-
priately the pressure viscosity curve at an arbitrarily chosen high
pressure.
a. Method of Analysis
u u- sD
y - u (89)
1.5P2
Tab (90)
227
V y
A u
Ur
I 1I
u0s
ffi6 a
T I d
I I
ig10 o atEIp
I I-66
228
U 1.'7 "
.[t
h -R' 1. 6 (91)
w(1- (93)
0b a2
The total frictional force acting on the ball in the x-direction, Fx,
is simply the sum of all dFx over the whole ellipse, i.e.,
yra b
F = dF f f w dy (94)
y-a -B
(95)
y 2 lu - 03sy- uJ
a A
F = S ydx dy (96)
-a -A
Slu + fy Ur dy (97)
-a
229
Finally, the frictional torque M about the normal to the ellipse is
z
calculated as follows:
y=a a A
=- y dF + T y x dx dy (98)
yz-a -a -A
a f w ( s2 b (1- Y
r--
b. Results
A few cases of pure spinning were studied to compare the present method
of analysis with the experimenta data by Allen, et al (137). The
input data are: E1 = E 2 = 3 x 10 psi, vI = V2 = 0.3, a= 0.92 x 10-4/
-5 2
psi, 01 =
0.028/F, Kf = 0.1 BTU/F-hr-ft, o g = 6 x 10 lb-secin
at 30 0C, To = 300 C, u = v = ur = 0, R = 0.25 in., r = c, and RG =
0.255, 0.275, 0.3 in. (corresponding to percent conformity = 51, 55;
60 respectively). The calculated spinning torque, Mz, at spinning
velocity a) 'of 1050 rpm as a function of maximum Hertz stress toge-
ther with the xperimental data taken from Ref. 137 is shown in F*
108. It is seen that the agreement is qtite good.
230
JrC% I C4
- a.
OD b
h..4
'44
0-M0
IL, N-4%b
CL 0
-j-
H-'m
ij ..
ql--ui A,
40 *
(0 1
CL
44
oM -
o 0
0 0
C5 C5
232
increases with spinning speed, but the rate of increase decreases at
high value of sptnntng speed, The frictional force in the rolling
direction, F , decreases with spinning speed. This is because that
the effect of sliding speed decreases as spinning speed increases,
The existence of a small hump for the later frictional force, F , is
doe to the fact that at small spinning speed Fy increases as thl re-
sult of the slight increase of the frictional coefficients. At very
large spinning speed the effect of lateral sliding speed vanishes.
I 23
233
SECTION V
CONCLUSIONS AND RECOMMENDATEON3
C,
Although technical work on elastohydrodynamic lubrication has been going on
for some eighteen years both here and abroad, there has been littlelopportun-
ity to compile and assess the significance and the interrelationshiks of the
diverse research results available today. This report represents qIch an
assessment and, in it, all of the available data on heavily loaded rolling and
sliding contact lubrication has been surveyed, compiled and summari ed. Out
of this assessment, an immediately useful design procedure has bee generated
and a plan for continued research evolved as outlined below.
1. CONCLUSIONS
a. Film Thickness
234
changes from one Gaussian distr.bution to another, as the
bearing wears.
c. Friction
d. Pressure Distribution
235
the pressure decays extremely rapidly to the ambient
pressure.
e. Temperatute Distribution
f. Stress Distributon
236
traction. A numerical solution based on the elasticity
solution for an elastic half-space was used in the calcula-
tion procedure to determine the three principal stresses,
the maximum shear stress as well as the octahedral stress
for any given point in the contact zone,
2. RECOMENDATIONS
a. Friction
/JI
I,'
of th:s report. The empirical correlation described ir.Sectic.-
III is based entirely on the limited data published by Johnson
and Cameron (129). Further experimental friction measarements
are required to verify this correlation over wider ranges of
pressure, temperature, shear rate, and lubricant properties.
Such empirical friction data is needed to serve as a-stop gap
design measure until a fuller understanding of lubricanr rheclo-
gy is achieved.
b. Film Thickness
d. Lubricant Rheology
238
high temperature, pressure, and shear rates. This then is the
S "bi problem, and unfortunately the most difficult. This lack
of knowledge about lubricant rheology under severe conditions
cannot simply be overcome by extending the range of measurements
using present experimental methods, as these methods are not
suitable under such operating conditions.
Fo
239
REFERENCES
1. Murphy. C.1. and Zisman, W.A., -Structural Guides for Synthetic Lubricant
Development," Ind. and Eng. Chem., Vol. 42, No. 12, Dec. 1950, pp. 2415-
2420.
7. Bowden, F.P. and Tabor, D., The Friction and Lubrication of Solids, Vol.2,
Oxford Press, 1964.
9. Wilcock, D.F. and Booser, E.R., Bearing Design and Application, McGraw-
Hill Book Company, Inc., New York, 1957.
10. Georgi, C.W., Motor Oils and Engine Lubrication, Reinhold Publishing Corp.,
New York, 1950.
11. "Viscosity-I and Viscosity-IlI" Lubrication, Anon., Vol. 47, Nos.. 1 and 2,
1961, Published by Texaco, Inc., New York.
240
/ .
13. Ibid. Estimation of Molecular Weight of Petroleum Oils from Viscosity
Measurements (D2502), Vol. 17.
17. Ibid? Calculating Viscosity Index from Kinematic Viscosity (D2270), Vol.
17.
18. Gerrard, J.E. and Philippoff, W., "Viscous Heating and Capillary Flow,"
4th International Congress of Rheology, 1963, Paper No. 51.
19. Gerrard, J.E., Steldler, F.E., and Appeldoorn, J.K., "Viscous Heating in
Capillaries: The 'diabatic Case," ACS Pecroleum Division Meeting, Chicago,
Illinois, Sept. 1964.
20. Gerrard, J.E., Steidler, F.E., and Appeldoorn, J.K., "Viscous Heating in
Capillaries: The Isothermal Wall Case," ACS Petroleum Division Meeting,
Atlantic City, N.J., Sept. 1965.
21. Novak, J.D. and Winer, W.O., "Some Measurements of High Pressure Lubri-
cant Rheology," Trans. ASME J. Lubr. Tech., Vol. 90F, pp. 580-591, 1968.
22. Bridgeman, P.W., Collected Works, Harvard University Press.
23. "Viscosity and Density of Over 40 Lubricating Fluids of Known Composition,
at Pressures to 150,000 psi and Temperatures to 425 F," ASME, New York,
1966.
24. Cameron, A., Principles of Lubrication, Wiley & Sons, New York, 1966.
25. Ed. by Ling, F.F., Klaus, E.E. and Fein, R.S., Boundary Lubricaticn - Al
Appraisal of World Literature, ASME, New York, 1959, Chapter 8 "Physical
Properties of Lubricants," by Appeldoorn, J.K'.
27. Barlow, A.J. and Lamb, J., "The Visco-Elastic Behavior of Lubricating
Oils Under Cyclic Shearing Stress," Proc. Royal Society, Vol. 253A, pp.
52-69, 1959.
28. Anderson, W.J. and Carter, T.L., "Effect of Lubricant Viscosity and Type
on Ball Fatigue Life," ASLE Trans., Vol. 1, No. 2, 1958, pp. 266-272.
241
30. Selby, T.W., "The Non-Newtonian Characteristics of Lubricating Oils,"
ASLE Trans.. Vol. 1, 1958, pp. 68-81.
32. Wright, W.A. and Crouse, W.W., Jr., "A New Concept in Generalizing Non-
Newtonian Fluid Flow Data," ASLE Trans., Vol. 8, No. 2, pp. 184-190, 1965.
34. Currie, C.C. and Smith, B.F., "Flow Characteristics or Organo Polysilo-
xane Fluids and Greases," Industrial and Engineering Chemistry, Vol. 42,
pp. 2457-2562, 1950.
37.. Lewis, P; and Murray, S.F., "Effect of Lubricant 'Propert.ies on the Coast-
down Characteristics of Instrument Bearings," Lubr. Engrg., Vol. 18, No.
.1, 1962.
39. Porter, R.S., Cantow, M.J.R., and Johnson, J.F., "Polymer Degradation VI-
Distribution Changes on Polyisobutene Degradation in Laminar Flow," Poly-
mer, Vol. 8, pp. 87-96, 1967.
42. Selby, T.W., "A Low Temperature Study of the Relationship Between Engine
Cranking Speed and Lubricant Viscosity," ASLE Trans., Vol. 2, No. 2, 1959,
pp. 208-216.
242
/
J 46. Dyson, A. and Wilson, A.R., "Film Thicknesses in Elastohydrodynamic Lub-
rication by Silicones Fluids," Proc. of I.M.E., Vol. 180, pt. 3K, pp.
97-112, 1965-66.
47. Barlow, A.J., Harrison, G., and Lamb, J., "Visco-elastic Relaxation of
Polydimethyl Siloxane Liquids," Proc. Royal Society, Vol. 282A, pp. 228-
251, 1964.
49. Borsoff, V.N., "On the Mechanism of Gear Lubrication," Trans. ASME, J.
of Basic Engrg., 1959.
50. Smith, F.W., "Some Aspects of Nonlinear Behavior in Lubricants Under Ex-
treme Stress," Trans. ASME. J. Lubr. Tech,, Vol. 90F, pp. 549-554, 1968.
52. Lamb, J., "Physical Properties of Fluid Lubricants: Rheological and Vis-
coelastic Behavior," Proc. I.M.E., Vol. 182, pt. 3A, pp. 293-310.,.
See alo, Roelands, C.J.A., Vlugter, JtE., and Waterman, H.I., "The Vis-
cosity-Temperature-Pressure Relationship of Lubricating Oils and its
Correlation withChemical Constitution," ASME Paper No. 62-WA-178.
56. Griest, E.M., Webb, W., and Schiessler, R.W., "Effect of Pressure on Vis-
cosity of Higher Hydrocarbons and Their Mixtures," J. Chem. Phys, Vol.
29, No. 4, pp. 711-720, 1958.
57. Hogenboon, D.L., Webb, W., and Dixon, J.A., "Viscosity of Several Liquid
Hydrocarbons as a Function of Temperature, Pressure and Free Volume," J.
Chem. Phys., Vol. 46, No. 7, pp. 2586-2598, 1967.
58. Cheng; H.S. and Sternlicht, B,, "A Numerical Solution for the Pressure-
* Temperature and Film Thickness Between Two Infinitely Long, Lubricated
Rolling and Sliding Cylinders Under High Loads." Trans. ASME, J. Basic
Engrg, Vol. 87D, pp. 695-707.
59. Dowson, D. and Whitaker, A.V., "A Numerical Procedure for the Solution
of the Elastohydrodynamics Problems of Rolling and Sliding Contacts Lub-
ricated by a Newtonian Fluid," Proc. I.M.E., Vol. 180, pt. 3B, pp. 57-71,
1965-66.
243
60. Foord, C.;... Hammann, W.C., and Cameron, A., "Evaluation of Lubricants
Using Opti:al Elastohydrodynamics," ASLE Trans., Vol. 11, pp. 31-43, 1968.
62. Deryagin, V.V., Abrikosove, I.I., and Lifshitz, E.M., 0Iarterly Rev., Vol.
10, p. 295, 1956.
63. Cameron, A. and Gohar, R., "Theoretical and Experimental Studies of the
Oil Film in Lubricating Pcint Contact," Proc. Royal Soc., Vol. 291A, pp.
520-536, 1966.
64. Needs, S.J., "Boundary Film Investigation," Trans. ASME, Vol. 62, pp.
331-345, 1940.
66. Fein, R.S. and Kreuz, K.E., "Discussion on 'Boundary Lubrication"', Pro-
ceeding of the NASA Symposium on Interdisciplinary Approach to Friction
and Wear, 1967, pp. 6.2.1-6.2.25.
69. Hopkins, V., Wilson, D.R., and Bolze, C., "Isothermal Bulk Modulus of
Selected Fluids to 700 F and 10,000 psig," Trans. ASME. J. of Basic Engrg.
1964.
70. Beerbower, A. and Green, D.F., "The Behavior of Lubricating Oils in Inert
Gas Atmosphere," ASI ITrans., Vol. 4, No. 1, 1966, p. 87.
71. Greeman, P.A., "Gyro Ball Bearings - Technology Today," Presented at 6th
AGARD Guidance and Control Meetings, "Inertial Navigation Components,"
Germany, May 7-9, 19686'
4.
72. Rhoads, W.L. arnd Sibley, L.B., "Supersonic Transport Lubrication System
Investigation," NASA CR-54313, Third Semiannual Report, May 20, 1966.
74. Murray, F.F. and Lewis, P., "Vaduum Operation of Oil-Lubricated Bearings,"
MTI 68TR65. October 1968, Work done for NASA Goodard Space Flight Center
under Contract NAS5-10433.
244
I 75. Fein, R.S. and Krautz, K.L., "Chemistry of Boundary Lubrication of Steel
by Hydrocarbons," AILE Trans., Vol. 8, No. 1, 1965, p. 29.
76. Klaus, E.E. and!ersh, R.E., "Paraffinic Resin - A High Quality Lubricant
and Newtonian .. Improver for Hydraulic Fluids," SAE Paper 630673,
Milwaukee, Wis onsin , September 1965.
77. "Super Refined Hydraulic Oils," Anon. Bulletin from Humble Oil and Refin-
ing Company, Houston, Texas.
3 78. Parker, R.J., Zaretsky, E.V., and Anderson, W.J., "Rolling Contact Lubri-
cation Studied with Polyphenyl Ethers at Reduced Pressures," NASA TN D-
3130.
80. Parker, R.J., Bamberger, E.N., and Zarets~y, E.V., "Evaluation of Lubri-
jIPresented
cants for High Temperature Ball Bearing Applicatiors," ASME Paper 67LUB21,
at ASLE/ASME Lubr. Conf., Chicago, October 17-19, 1967.
81. Godfrey, D., "Chemical Changes in Steel Surfaces During Extreme Pressure
82.Lubrication," ASLE Trans., Vol. 5, No. 1, 1962, p. 57.
82. Grubin, A.N. and Vinogradova, I.E., Central Scientific Research Institute
for Technology and Mechanical Engineering, Book No. 30 (Moscow); D.S.I.R.,
Trans. No. 337, 1949.
83. Crook, A.W., "The Lubrication of Rollers - I," Phil. Trans. Series A 250,
I 387, 1958.
84. Crook, A.W., "The Lubrication of Rollers - II. Film Thickness with Rela-
8tion to Viscosity and Speed," Phil. Trans., Series A 254, 234, 1961.
245
91. Dowson, D. and Higginson, G.R., "A Numerical Solution to the Elastohy-
drodynamic Problem," F. Mech. Engrg. Sci. 1 (1), 6, 1959.
92. Archard, u.D., Gair, F.C., and Hirst, W., "The Elastohydrodynamic Lubri-
cation of Rollers," Proc. R. Soc. Series A 262, 51, 1961.
97. Martin, H.M., "Lubrication of Gear Teeth," Engineering, London, Vol. 102,
199, 1961.
98. Weber, C. and Saalfeld, K., "Schmierfilm bei Walzen mit Verformuny,"
Zeits. ang. Mech., Vol. 34, 64, 1954.
106. Archard, J.F. and Kirk, MJ., "Film Thickness for a Range of Lubricants
Under Severe Stress," J. Ach. En'rg. Sci. 6, 101, 1964.
246
i 107. Archard, J.F. and Cowk1i. E.W., "A Simplified Treatment of Elastohy-
109. Hingley, C.J., Schmidt, W.E., and Tataiah, K., "A Study of the Gecmerry
of Elastohydrodynamic Films in Point Contact," SKF Report AL66L035.
112. Osterle, J.F. and Stephenson, R.R., "A Direct Solution of the Elastohy-
drodynamic Lub:ication Problem," Trans. Amer. Soc. Lubr., Engrs,. 5, (2)',
365, 1962.
113. Wernick, R.J., "Some Computer Results in the Direct Iteration Soluction
of the Elastohydrodynamic Equations," MTI.Technical Report 62TR38, 1963.
II
114. Dowson, D. and Whitaker, A.V., "The Isothermal Lubricatfon of Cylinders"
Presented at the ASME-ASLE International !,brication Conference (Wash-
K rington, D.C.) Oct. 1964, ASLE Preprint 64-LC-22.
115. Sternlicht, B., Lewis, P., and Flynn, P., "Theory of Lubrication and
247
121. Bell, J.C., Kannel, JW., and Allen, C.M., "The Rheological Behavior of
the Lubricant in the Contact Zone of a Rolling Contact System," J.
Basic Engrg. Trans. ASME 86, Series D, No. 3, 423, 1964. 1
123. Cameron, A. and Newman, A.D., Proc. Conference of Steam Turbine Research
and Development, London: Institute of Mechanical Engineers.
129. Johnson, K.L. and Cameron, R., "Shear Behavior of Elastohydrodynamic Oil
F Ilm at High Rolling Contact Pressures," Proc. of Instn. of Mech. Engrs.
1967-1968, Vol. 182, p. 307.
130. Jefferis, J.A. and Johnson, K.L., "Sliding Friction Between Lubricated
Rollers," Proc. of Instn. of Mech. Engrs., 1967-1968, Vol. 182.
132. Dowson, P. and Wholmes, T.L., "Effect of Surface Quality Upon the
Traction Characteristics of Lubricated Cylindrical Contacts," Proc. of
Instn. Mech. Engrs., Vol. 182.
133. Barlow, A.J. and Lamb, J., "The Visco-elastic Behavior of Lubricating
Oils Under Cyclic Shearing Stress," Proc. R. Soc., A, 253, p. 52, 1959.
134. Kannel, J.W. and Walowit, J.A., "Simplified Analyses for Tractions-Be-
tween Rolling-Sliding EHD Contacts," to be published in ASME Trans.
J. Lubr. Tech.
135. Dietrich, M.W., Parker, R.J., Zaretsky. E.V., and Andefson, W.J., "Con-
tact Conformity Effects on Spinning Torque and Friction," ASME Paper
No. 68-Lub-l0.
248
I\
1 136. Parker, R.J., Zaretsky, E.V., and Anderson, W.J., "Spinning Friction
Coeffictents with Three Lubricants," J. of Lubr. Tech. Trans. of &SME,
Series F, Vol. 90, No. 1, January 1968, p. 330.
137. Allen, C., Townsend, D., and Zaretsky, E.V., "Elastohydrodynamic Lubri-
cation of a Spinning Ball in a Non-Conforming Groove," ASME Paper No.
169-LubS-ll.
138. Hertz, H., "Uber die Beruehrung fester elastischer Koerper," 1881;
"Ueber die Beruehrung fester elastischer Koerper und ueber die Haerte,"
1882: reprinted in "Gesammelte Werke von Heinrich Hertz, Vol. 1, pp.
155-173, 174-196, Leipzig, 1895, English translations in H. Hertz,
"Miscellaneous Papers" translated by D.E. Jones and G.A. Schott, pp.
146-162, 163-183, MacMillian, New York, 1896.
139. Smith, J.O. and Liu, C.K., "Stresses Due to Tangential and Normal Loads
i on an Elastic Solid with Application to Some Contact Stress Problems,"
I Trans. of ASME J. of Applied Mechanics, June, 1953, pp. 157-166.
142. Thomas, H.R. and Hoersch, V.A., "Stresses Due to the Pressure of One
Elastic Solid Upon Another," University of Illinois Bulletin 212, July
15, 1930.
146. Greenert, W.J., "The Toroid Contact Roller Test as Applied to-the Study
of Bearing Materials," J. of Basic Engrg., Tr~ns. ASME, March 1962.
* 148. Cheng, H.S. and Wernick, R.J., "Stresses in Lubricated Rollers Consid-
ering Arbitrarily Distributed Normal and Tangential Loads," MTI 65TR61,
* December 196f.
f 149. Greenwood, J.A., "The Area of Contact Between Rough Surfaces and Flats"
r J. of Lubr. Tech., Trans. ASME, Vol. 89, Series F, 1967, p. 81.
249
*t
150. Tallian, T.E. et al, "Lubricant Films in Rolling Contact of Rough Sur-
faces," ASLE Trans., Vol,. 7, 1964. 'k
151.- Greenwood, J.A. and Williamsor, J.B.P.,, "The Contact of Nominally Flat
Surfaces," Proc. of R. Soc., London, Vol. 295, Seiies A, 1966, p. 300.
152. Williamson, J.B.P., Pullen, J. ., and Hunt, R.T., 1969 Surface Mechanics
New York, ASME.
155. Boness, R.J., "The Effect of Oil Supply on Cage and Roller Motion in a
Lubricated Roller Bearing," J. of Lubr. Tech., Trans. ASME, Series F,
Vol. 92, No. 1, 1970, pp. 39-53.
157. Shevchenko, R. and Bolan, P., "Visual Study of Ball Motion in a High-
Speed Thrust Bearing," SAE Paper No. 37, January 14-18, 1957.
158. Poplawski, J.V. and Mauriello, J.A., "Skidding in Lightly Loaded High-
Speed Ball Thrust Bearings," ASME Preprint 69-LubS-20.
161. Allan, R.K., "Rolling Bearings," Sir Isaac Pitman & Sons, Ltd., 1946.
162. Palmgren, A., "Ball and Roller Bearing Engineering," S.H. Burbank & Co.,
Third Edition, 1959.
163. Kaufman, H.N. and Walp, H.P., "Interpreting Service Damage in Rolling
Type Bearings - A Manual on Ball and Roller Bearing Design," Amer. Soc.
of Lubr. Eng., 1953.
164. Bisson, E.E. and A,,derson, W.J., "Advanced Beariing Technology," Nat'l
Aeronautics and Space Admin., Wash., D.C., NASA SP-38, 1964.
166. Wren, J.J., and Moyer, C.A., "Modes of Fatigue Failures in Rolling
Ele.ent Bearings," Proc. Instn. Mech. Engrs., 179, Pt. 3D, 1964,
-. 236.
250
L
I
167. Edigaryan, F.S., "Basic Types of Wear in Anti-Friction Bearings,"
Russian Engr. Jour., XLVI, 4, 1966.
168. Tallian, T.E., "Special Research Report on Rolling Contact Failure Con-
trol Through Lubrication, SKF Report AL66Q028, Sep 1966.
171. O"Connor, J.J. and Boyd, J., "Standard Handboolo of Lubrication Engineer-
ing," McGraw Hill Book Co., New York, 1968.
172. Moyar, G.J., and Morrow, A.V., "Surface Failure of Bearings and Other
Rolling Elements," Univ. of Illinois Engr. Experiment Station, Bulletin
469, 1964.
174. Eschmann, Paul, "Rolling Bearing Wear Life," ASME Paper 67-WA/Lub-2, 67.
175. Brothers, B.G., and Hailing, J., "Effect of Geometric Conformity Between
I Rolling Bodies on the Slip and Wear in the Contact Region," Proc.
Instn. Mech. Engrs., 3rd Cony. Lub and Wear, 179, Pt. 3, 1964, p. 134.
176. Hailing, J., and Brothers, B.G., "Wear Due to the Microslip Between a
S Rolling Body and its Track," ASME Paper 64-Lub-30, 1964.
177. Tallian, T.E., McCool, J.I., and Sibley, L.B., "Partial Elastohydro-
dynamic Lubrication in Rolling Contact," Elastohydrodynamic Lubrication
Symposium in Leeds, Paper No. 14, Inst. of Mech Engrs, London, 1965.
3 178. Tallian, T.E., Brady, E.F., McCool, JoI., and Sibley, L.B., "Lubricant
Film Thickness and Wear in Rolling Contact, ASLE Trans, 8, 4, 1965,
p. 411.
179. Landen, E.W., Slow Speed Wear of Steel Surfaces Lubricated by Thin
Oil Films," ASLE Trans. 11 6-18 (1968).
180. Queener, C.A., Smith, T.C. and Mitchell, W.L., "Transient Wear of
Machine Parts," Wear 8, 391 (1965).
181. Devine, M.J., Lamson, E.R., Cerim, J.P., and Carrol, R.J., Session 4,
Lubrication and Wear Conference, London, 196 .
182. Cocks, M.F., "The Formation of Wedges of Displace Metal Between Sliding
Metal Surfaces," Wear 8, 1965, p. 85.
251
183. Cocks, M.F., "Shearing of Junctions Between Metal Surfaces," Wear, Jul/
Au;, 9, 4. 1966, p. 320.
185. ibrahim, M., Cameron, A., "Oil Film Thickness and the Mechanism of
Sco-ing in Gear Teeth," Paper 20, Proc. Conf. Lub. and Wear, Bournemouth,
Inst. Mech Engrs 1963.
186. Peacock, L.A., and Rhoads, W.H., "Extreme Temperature Aerospace Bearing
Lubrication Systems," NASA CR-72446, 1968.
187. DeGruch, V.J., "Development of an Edge Type Disk Machine and Preliminary
Studies of Vai-ous Gear Materials Lubricant Combinations," Proc. Conf.
Lubrication and Wear, Bournemo.th, I:st. Mech Engrs, 1964, p. 160.
188. Leach, E.E., Kelley, B.W., "Temperature - The Key to Lubricant Capacity"
ASLE Trans, 8, 1965, p. 271.
189. Matveesky, R.M., "The Critical Temperature of Oil with Point and Line
Contact Machines," Jour Basic En:., Trans ASM 87, 1965, p. 754.
191. Fein, R.S., Rowe, C.N., Krenz, K.L., "Transition Temperatures in Slid-
ing Systems, ASLE Trans 2, 1, 1959, p. 50.
192. Fein, R.S., "Transition Temperatures with the Four Ball Machine," ASLE
Trans 3, 1, 1960, p. 34.
252
I
199. Courtney-Pratt, J.S., and Eisner, E., "The Effect of Tangential Force
on the Contact of Metallic Bodies," Proc. Roy Soc., Ser A, 238, 1957,
p. 529.
200. Parker, R.C. and Hatch, D., "The Static Coefficient of Friction and
the Area of Contact," Proc. Phy. Soc. 63, 1950, p. 185.
201. Bowden, F.P. and Tabor, D., "The Friction and Lubrication of Solids,
Part II, Oxford (London) 1964.
203. Cocks, M., "Role of Displaced Metal in the Sliding of Flat Metal
Surfaces," J. Appl. Phy. 35, 6, 1964, p. 1807.
205. Johnson, K.L., "A Review of the Theory of Rolling Contact Stresses,"
Wear 9, (1966) 4-19.
206. McKelvey, R.D. and Moyer, C.A. "The Relation Between Critical Maximum
Compressive Stress and Fatigue Life Under Rolling Contact," Paper No. 1
Proceedings of the Symposium on Fatigue in Rolling Contact, Institution
of Mechanical Engineers, London, 1963.
207. Tsukizoe, T. and Hisakdo, T., "On the Mechanism of Contact Between Metal
Surfaces - The Penetrating Depth and the Average Clearance," Journal of
Basic Engineering, Trans. ASME, Series D, Vol. 87, No. 3, Sept. 1965,
pp. 666-674.
209. Dawson, P.H. "The Effect of Metallic Contact and Sliding on the Shape
of the S-N Curve for Pitting Fatigue," Paper No. 4 Proc. Symp. Fatigue
.1 in Rolling Contact, I. Mech. E. 1962.
253
213. Cheng, H.S., "A Numerical Solution of the Elastohydrodynamic Film
Thickness in an Elliptical Contact," Journal of Lubrication Technology,
Trans. ASME, Vol. 92, Series F, No. 1, pp. 155-162.
214. Dowson, D. and 'Whitaker,A.V.. "A Numerical Procedure for the Solution
. of the E1zszohydrodynamic Problem of Rolling and Sliding Contacts Lubri-
cated by a Newtonian Fluid," Lubrication and Wear Group Symposium on
Elastohydrodynamic Lubrication, Paper 4, Institution Mech. Engrs.,
London.
215. Jones, A., "Ball Motion and Sliding Friction in Ball Bearings", ASME
Transactions, Journal of Basic Engineering, March 1959.
I5
APPENDIX I
TYPICAL PROPERTIES OF VARIOUS CLASSES OF LUBRICANTS
Mineral oil base lubricants are available in a range of bulk viscositiEs. Some
typical grades are listed below:
Automotive:
SAE 1OW 41 6.0 0.870 410 -15
30 114 11.3 0.891 460 - 5
Industrial gear:
80 69 7.9 0.934 365 -25
90 287 20.4 0.930 450 -10
140 725 34 0.937 500 0
Turbine:
Grade 1010 10 2.5 0.864 280 -70
Light 32 5.4 0.872 410 0
Medium 65 8.2 0.877 455 10
Heavy 99 10.8 0.885 470 10
Hydraulic:
Light 32 4.8 0.887 370 -45
Medium 67 ' 7.3 0.895 405 -15
Heavy 196 14.0 0.901 495 10
Extra low t'mp. i4 5.2 0.884 230 -80
Wide temp. 56 I10.5 0.871 310 -45
Contidued
255
TABLE I-I TYPICAL VISCOSITIES OF STANDARD PETROLEUM OILS - CONTINUED
It is important to note that these are nominal physical property data. Varia-
tions from batch to batch are normally on the order of 5 to 10 percent. Pure
hydrocarbonswith known structure and physical properties are available, but
these are too low in viscosity to be used as lubricants.
Special oil blends can also be obtained from some of the major oil companies
which have the same bulk viscosity at some specified temperature, say 100 F,
but widely different viscosity-temDerature and viscosity-pressure character-
istics. These oils are prepared by blending to obtain predominantly aromatic
or paraffinic structures.
Synthetic hydrocarbons are also available from some of the oil companies.
These offer the advantages of having known chemical structures and higher
purities. Table I-III shows some physical property data on typical oils from
one source.
Although the mineral oils are of major importance from a commercial stand-
point, the uncertainties regarding purity, chemical structure, etc., make
them poor candidates for laboratory studies when the goal is to gain e better
understanding of some phenomenon.
256
TABLE I-I TYPICAL PHYSICAL PROPERTIES OF SUPER-REFINED MINEFAL
OIL (FR& REF. 77) CONFcRMING TO .LO 7277
257
"A'A
TABLE I-III TYPICAL PROPERTIES FOR MOBIL SYNTHETIC HYDROCARBON FLUIDS*
PROPERTY A B C
I
I
258
I
1-2. POLYPHENYL ETHERS
The polyphenyl ethers were the result of extensive effoxts by many investiga-
tors to develop high temperature, radiation-resistant lubricants. Chemically,
the two go hand in hand, i.e., those organic materials with the best resist-
ance to oxidation or thermal degradation are also the most resistant to
radiation damage. A typical structure for the polyphenyl ethers is as follows:
.41C - 0 , C - 0 - C- 0 - -
SC- 0 C C
The density and bulk modulus of the polyphenyl ethers is high, indicating
close packing and low compressibility for these oils. The 5P4E and the 4P3E have
high pour points ( + 1OF to +35F) and relatively poor viscosity-temperature
behavior. In many respects, the polyphenyl ethers are similar in structure
to the aromatic base petroleum oils (except for the ether linkages, and the
fact that they are pure compounds). For this reason, the physical behavior
of oils in each of thexQ two classes is similar.
The results of several investigations (Refs. 78-80) have indicated that the
polyphenyl ethers are poor lubricants in rolling contact bearings. This was
ascribed to the fact that these oils do not form elastohydrodynamic films at
high unit stresses. However, the reasons for this behavior are not clear at
this time.
ypical physical property data on the four ring and five ring polyphenyl
ethers are given in Table A-IV. No data on viscosity-shear behavior were
found.
1-3. ESTERS
259
TABLE I-IV TYPICAL PHYSICAL PROPERTIES OF POLYPHENYL ETHERS
4-Ring 5-Ring
Viscosity-Centistokes
280F 13,000 cs
1000F 63.2 355
2100F 6.0 13.1
4000F 1.31 1.25 (5000F)
Density - g/nil
Specific Heat-BTU/lb/F
0
100OF 0.36 ,0.40
Thermal Conductivity
BTU/hr/ft/ F
100O -- 0.078
0 0.077
400 F --
260
:i
!
Most of the synthetic ester lubricants are being made from organic acids
having two aci groups (dibasic acids). Three of the most important members
of this class are sebacic acid, adipic acid and azelaic acid. These have
the following structure:
HOOC C - C C - COOH
"" H H
n
By reacting these di-basic acids with suitable alcohols, diesters can be pre-
pared with very satisfactory physical properties. One alcohol, 2-ethyl
hexanol, has been used widely for this purpose. It has the structure:
H H H H H H
H C C - C - C - C - C - OH
H H H H HHCH H
HCH
H
Using the letter R to represent the hydrocarbon portion of the dibasic acid,
and R' to represent the hydrocarbon chain of the alcohol, the reaction be-
tween the dibasic acid and two molecules of alcohol is:
If the dibasic acid is sebacic and the alcohol is 2-ethylhexanol, the- result-
ing diester is di (2-ethylhexyl) sebacate. This diester is typical of those
MIL-L-7808 formulations which contain diester base fluids. By compounding
the oil with a suitable additive package, a satisfactory lubricant can be
obtained for aircraft gas turbine engines.
In general, the esters have physical properties which are very similar to the
petroleum oils. This is not unexpected, since they are predominantly straight
chain hydrocarbons with the same atomic make-up. They have lower pour points
and better viscosity-temperature properties than conventional mineral oils.
Some typical datp are given in Table I-V.
261
TABLE I-V TYPICAL PHYSICAL PROPERTY DATA ON ESTER'BASE LUBRICANTS
Viscosity-Centistokes
Specific Gravity
+10F.936 0.987
+210%F .892 0.939
+300%F 0.854 0.904
Thermal Conductivity
BTU/hr/ft 10r
(a) 0-67-22
V (b) 0-70-4
262
1I
1-4. SILICONES
The silicone oils are linear polymers with a structural shape that is similar
to the straight chain hydrocarbons. Instead of a chain of - C - C-carbon
atoms, the silicones have a chain of repeating - -0-Si-O-Si- units for a
backbone. Some typical structures are shown below:
CH OH
13 1 .
H3C - CH
S-0- Si - Si - CH3 where n can be any integer
3 Ch 3 Ck3
C n
These molecules have weaker intermolecular forces than the hydrocarbons, and
the molecular chains are much more flexible. Because of their wide molecular
spacing, the silicones are easily compressed and have low bulk moduli. A
generalized relationship for the bulk modulus of silicones s a function of
pressure, temperature and density is presented in Ref. 68. !The most outstand-
ing characteristic of the silicones is their small change i viscosity with
temperature. The dimethyl silicones are particularly good prom this respect.
*. Pour points of most of the dimethyl silicones are on the order of -60 F. By
synthesizing molecules wit. branched chains, the pour point can be reduced to
below -100 F.
The -C-Si bond in the dimethyl silicones is the weak link as far as oxidative
or therial degradation is concerned. To improve the resistance of the molecule
to temperature effects, phenyl groups can be substituted for methyl groups.
This results in some loss in favorable properties, such as pour point and
viscosity-temperature characteristics.
Aside from their relatively high cost, the most outstanding disadvantage of
'the &ilicones is their poor boundary lubricating characteristics for steel
sliding on steel. Considerable effort has been expended on this problem
with little success. The effect is specific since it has been found that
steel vs. bronze combinations are effectively lubricated by the silicone oils.
263
It has been hypothesized (Ref. 81) that the silicone oil decomposes at the
metal surface to form inorganic silicides which inhibit conventional anti-
wear or extreme-pressure additives from reacting and providing beneficial sur-
face filns.
I
!
264
1-5 EXTRACTS FROM SPECIFICATION MIL-L-9236B, 4 MARCH 1960
Intended Use: Lubricating oil for specific models of aircraft gas turbine
engines.
Requirements
r Evaporation loss: Not more than 15% by weight when tested at 400 F
fcr 6 1/2 hours.
* Gear Tests:
e Bearing stabilization temperature: The oil shall not cause the bee-
ing temperature to stabilize over 500F when tested in the ERDCO
bearing test rig. Time required to reach stabilization temperature
shall be less than 60 minutes. Bulk oil temperature 4000F.
* Used Oil control: Increase in oil viscosity at 100 0 F after the engine
test shall not exceed 75%. Must be as good or better than the refer-
ence oil.
265
V
o Additional requirements include:
266
1-6 EXTRACTS FROM4 PROPOSED SPECIFICATION MIL-L-27502. SUPERSEDING MIL-L-
9236B, 4 MARCH 1960
Lubricating Oil, Aircraft Turbine Engine, Ester Base
Intended Use: For specific engines and power equipment requiring an e.ter
base oil with an approximate temperature range capability of
-40 to +428 0 F.
Requirements
Elastomer compatibility
Corrosion characteristics
Deposition characteristics
Compatibility, turbidity
267
APPENDIX II
ANALYSIS OF PROTRUSION WIDTH
d 6
= uo(Ul + U2) (h -h*) (11-I)
3
dx h
e- p
where q = 1-
d h*2 H3
Referring to the diagram below, the limits in the exit region are from
I-e to = 1.0.
h*h
e e
Figure II-i
268
Z'i
1.0
H - I ( 1-4)
e )'
1- = 6p, ( Ul + U2 b
0 -1
6a o(U I
2
+ U2)b
J i.0 H I
3--- dx (11-5)
0~x
=H -
-- .( -- ) d d 'Equation(II- becomes
*InhE'*6ducing x' =- dx
=
6c Lo(U1 + U2)b H - I d' - (11-6)
2 3
h
h* 0o H
--
LetI1
Let I I
= -J dH
C3
It follows that
.2
__ 2 12 (11-7)
2 Rx U+
6,y uu0(U U
: I .x e R, R' 22
u UU)\Ji I
\ E' I.
2 R'
"*
+
H* .2 0(U I U2) 1.4 .26
269
- 1 H* 1 (11-10)
241 (1U + u2)
1 X ,) (0 E'R' -
x x
Substituting (11-9) into (II-10), one obtains the final expression for the
protrusion width,
- K 2 L0(U1 + U2)
e - 241 ( 2 E'R' ) .
,x
2.70
-U. APPENDIX III
ANALYSIS OF AREA OF CONTACT
III-1 ANALYSIS
Smooth-Film Profiles
m m 0
The isothermal film profiles obtained by Dowson and Higginson (153) may
be used to estimate the value of 6 and H for given load and speed para-
meter..
For a ground surface, Williamson (152) has shown that the height distribution
or any arbitrary cross-section of the surface profile is very close to
Gaussian distribution. This evidence lends support to the assumption used
earlier by Tsukizae and Hisakado (207) to investigate the penetrating depth
and the average clearance between two contacting surfaces.
If the height 2.ta !n Fig. III-2a is divided into a number of equal segments
Au, then the number of crosses made by a straight line at a distance u from
the median line within a sample length L is represented by the distribution
histogram Fig. III-2c. The probability density of the surface height distri-
bution can be written as,
271
of Contact
2b -_ _ _ _ _ _ _ _ _ _ _
2.72
J f (u)
= 1 exp [ (f )2] 111
where
ar = standard deviatiGn
From the height distribution, one can readily determine the area of intersec-
tion between an ideal, rigid, plane surface and a rough surface at f = f as
follows:
f(u) du
A = L . f, u (111-2)
+t f (u) du
where Acis the area of intersection per unit width. If one ignores the change
of the surface profile due to the plastic deformation, then the area of inter-
section is simply the contacting area. Fig. III-2b shows the variation of the
contacting area curve as a function of u. This is generally known as the
Abbott's bearing area curve.
where
"[f (u) du
S..F (h) = h
J"f (p) du
It follows that the total area of contact per unit width becomes
Ac = F (h) dx (111-4)
a
273
Median Xu
T-11 du
Width Direction
Length Directio
/ R Median Plane -
274
Introducing
c 2a'
X
a
h
1
Ac F/ (h) dx- (111-5)
c 2
-(H((111-7)
Ji"f (t) dt
! F () =
ff (t) dt
275
!f~t) dt -of (t)dt
2f~f( t) dt
2 24(t)II-
2 !,
where
t2
1 2
J 0 ,,2 e dt (IIl-10)
h
= specific filmo thickness, _C0
Gb
h
Hm = exit film ratio, -
0
276
-. of Contact
b b
VN
Iie
277
APPENDIX IV
SIMPLIFIED ANALYSIS OF FRICTION
f PZ h e et ( + 1) d(2E (IV-l)
TTpHZ h0 +1)l
LIo 81 u s 2
where 4 8K ) e
= e) p - o(TTo)
a = pressure-viscosity coefficient
In the derivation of Eq. (IV-l) the effect of the temperature rise of the
contacting surfaces h, been neglected. Crook (120) shows the surface temper-
ature rises are not very imp'jrtant to the calculation of frictional coefficient
From Eq. (IV-l) one sees that the coefficient of friction, f, is governed
by the following three paraseters:
218
ru
LL0 US
G1= P h
Rz 0
2
G8Kf ' (IV- 2)
G3 = *PHZ
It should be noted that all the frictional coefficient graphs are for a given
inlet lubricant temperature 30C. Change of inlet temperature affects the
frictional coefficient due to the dependence of the limiting shear stress
upon the inlet temperature.effect.
279
APPENDIX V
CALCULATION OF APPROXIMATE ELASTOTI ORODYNAMIC PRESSURE PROFILE
The entire conLact region is divided into three regions: the inlet,
the central, and the exit region as shown in Fig. V-i. The inlet
region is bounded by -4.0<x < _, where 7 is the point where the
pressure profile begins to matc# the Hettaian pressure profile. It is
assumed that the transition pressure between the inlet and the central
region occurs when
SPA
e = 150 (-)-
(V)
or
'PA = 5
- 5 b
_ 5- (V-2)
It follows that
Xa J- (-VA-
__ (V-3)
At the far side of the inlet region (-4.0 < x < -. 2.0), form equal divi-
sions are used. It is followed by non-uniform grids as shown in
Fig. V-1. The central region is bounded by 3 < < (1 - 27) with 12
non-uniform spacings before the center line and 4 uniform spacings
after the center line. The exit region is divided into 8 uniform
spacings of e/4.
280
Ix
__ _ 17
0 aI bz - 0
0- JO '-a .44
C14,
00
0a
41 'AT[ *bZ
0 C,, J
4)0
K: JO
IE ba1 '-a
9/1 0
k "4
jo tqa*b
Zc
m
~*o;0 z-
I,,
1- - 6s2 (h*) 2+
2-4.U2)b fH 0 H3-1 -
1 - e-u" "
It follows that the ratio becomes
i- e- A
1 e S (V-5)
-e - PA SA
where
SH- f -C Hl H3 dx and
SH'-dx (V-7)
)p = - i_2 (V-9)
IT PZ 18
282
!-
To determine the peak pressure pB' it is further assumed that the gain-
ing in load represented by the area A
bolic pressure can be written as
p(x) =
5.
D
x)
Accordingly, the rising para-
V- )+
XD
n +
D (V-lO)
where
= D and -x -2
PD T D
"* where
2 -\l2
PBB = e-(I-
)
S - It follows that the gaining in load WG before the spike can be written
as
2- (V-12)
283
of contact
284
I ., ~ and the loss in load WL after the spike is
IL I ':
" .~ l- x p~Bd (V-t3)
WL-D Dx1 dx
I e
For extremely lightly loaded cases, e may exceed 0.5. In which case,
the method described here will break down and there will be an error
message printed out in computer program and both the pressure, the tem-
perature, and the stress calculations will be bypassed in the program.
285
APPENDIX VI
INLET TEMPERATURE CALCULATION
In the method the temperature profile of the lubricant in the inlet region is
assumed to be a parabolic in terms of surface temperatures and center film
temperature, i.e.,
2
2
(VI-I)
T(x y) = 1 2
(2TI +2T Ch
2 - 4T)2
- +(-3T "
T2 + 4T,c ) h + TI,
T,
h h h
cff (pu PUK
+pvT)dy -y f 2dy f Y dx
l
h
u L
p oT
dy (VI-2)
In the above equation, c, o, K and i are the specific heat, density, heat
conductivity and viscosity of lubricant respectively, u and v are the velocity
components in the x and y direction respectively.
Since the surface Lemperatures do not rise appreciably in the inlet region,
T and T ere assumed equal to the constant inlet temperature T . Lubri-
cant denity is set to be constant. The film thickness is given by tC
Hertzian displacement profile. In non-dimensionalized form the profile is
286
as follows:
h = l+B x x -1)-
22
n ( ++- x -1)
] (VI-3)
where
h
h
0
h = film thickness
- 22U UI 2
3 (l
1 U1 (i
u 3(l- )y + [-- - y)
y+u1+u2 (I-41
-N
* where
u = U
u =flow velocity
y h
Under the assumptions mentioned above and with the temperature profile Eq. VI-
1, the non-dimensionalized integrated energy equation takes the following
-287
(I
form:
1 .2
where
PoCf(U 1 + U )b3
= 16 Kf R 2
2
0 (U1 + U2 )
Qm - 2KJo -
LI M viscosity L at T0
T
LI
IL
LO
This equation is integrated, using Eqs. :(Vj.,3).:and (VI-4) for h and u respec-
tively, to qbtain the center film temperature Qc" The viscosity function
exp 8 is employed in the integration.
0
The integration is started where h - 6. The' temperature there 'is assumed equal
to the inlet temperature T o. The integration stops near the entrance to the
contact zone where there is no longer reverse flow. This position is calcu-
lated from the velocity profile Eq. (VI-4).
The above calculation is contained in the subroutine THILT in the main elasto-
hydrodynamic computer program.
288
,,, I
3 F .
APPENDIX VII
EIASTOHYDRODYNAMIC PERFORMANCE COMPUTER PROGRAM
The Fortran statements for the input routine are listed below.
Input Routine
289
Card 3 FORMAT (8E10.3)
CONS1 = thermal conductivity of body number 1, BTU/F-hr-ft.
CONS2 = thermal conductivity of body number 2, BTU/F-hr-ft.
CPSI = specific heat of body number 1, BTU/lb-F.
CPS2 = specific heat of body number 2,3 BTU/lb-F.
RHOSI = density of body number 1, lb/in .
RHOS2 = density of body number 2, lb/in 3 .
Card 4 FORMAT (8E10.3)
EPST = calculating the convergence factor for temperature
at each grid point, recommended value is 0.0001.
EPSRD = this is a convergence factor for the iterations to
determine the secondary pressure viscosity coeffi-
cient, recommended value is 0.0005.
290
K NSTRS control for stress calculations
NSTRS = 0, calculating stresses
NSTRS = 1, bypass stress calculations
Card 6 FORMAT.(1715)
THIS CARD IS ONLY REQUIRED FOR NSTRS --0.
.=Jew + ( + yp)11
291
CPF = specific heat of the lubricant, BTU/lb-F.
RHOF = density of the lubricant, lb/in 3 .
b. Thermal Parameters:
PO (U2
QM 2KfT
) 3
CONV PCf (U2 + U1
16Kf (R') 2
292
Dl =
D2C
= ( Kl)" (:w' 3
UP1: /
c. Contact Dimensions:
_ (U2 + U1
UBAR = U = "-
2E'R'
X
G a
AFA = 1PH
GHA
BTA = 2
1 TPHZ
T
0
""T
WBTA = "--
o
Equivalent line contact load ; for the elliptical contact.
PHZ/ED =.
E'R' 4w
"X
"" ~PHZ/ED =PZ
El
4w
2=
Slide to Roll Ratio U2 + U
I
293
e. Reduction Factors:
Self explanatory.
Protrusion width e
b
h - hm
Protrusion Depth = h
h
0
h
Specific Film
XBAR =.
b
h
Normalized Film
h 00
Pressure =
11
2, PHZ
Note: For extremely lightly loaded cases such that the protrusion
width t exceeds 0.5, the pressure calculation in the program breaks
down (see Appendix V). A message will be printed out for such
cases and the pressure, the temperature, and the stress calcula-
tions will be bypassed.
1
294
!
7f
k. Stress Calculations:
XST x
b
= inlet viscosity.
T = inlet temperature.
295
L%
A "J A
I - S
1 -T21
I(j-I-- I
L 1
e = protrusion width -
eb
296
PROGRAM PN45(INPUTOUTPUTTAPE5=INPUTTAPE6SOUTPUTTAPE22)
C PERFORMANCE CALCULATION OF ELASTOHYOROOYNAMIC
I C
C
CONTACT
297
C R~EAD INPUT DATA
DIMENSION AMDA(24)oFMA(24),AOBA(24)
DIMENSION SLIPA(9),CCA(9),EEA(9)
COMMON K1K*~KtOKt~t~K*RN*I
COMMON RX1,RX2,RYlHY2,RMS1,RMS2,E1,E2,vI9V2
COMMON CONS1,CONS2.CPS1,CPS2,RHOS1,RHOS2,WTTEMIU2,UI
COMMON VIS.ALF,8ETGAMCONFCPF ,RHOFPSPSTAR
COMMON AFA,8TA*GMA*A98vAM9SLIP9 RXP.RYPEP.OM.CONVO(2),COMP,AHZ,
lPHZ#WBARvUHARPHiZEGG.EOOWBOA 9Ww
COMMON FCOEEBARFEETFEE5,HMINtSPFILMA8ARSHM
COMMON EPSTEPSRDDTF.TWF1IrIRDUITINRITE
COMMON Xi(503,hE(50),P(SO),T(50),TI(2,50),F(2,50),SNS(50)
COMMON XSTTXDEXtYSTT.YDEYXST(31),Y(31) ,JJFSoKKFSNTHENOPL
NRz5
NWu6
DATA(AMDA(1),I1292)/ 0.0,O.05,O.1,0.15,0O.225O.3,O.35,0.4,0.45
1'0.5,0.55'0.6.O.659O.790.75,0.8,0.85,0.9,0.92.0.9490.96,0.98e1.O/
DATA(FMA(I).I=1,24)/1.0,1.03591.072.1.112,1559,19591.24,1.299
11.35,1.4191.48591.b6.1.65.1 .7791.89592.06592.29,2.60.3.O956.3o3973 -
293*8275.4,5319*5*94696*/
DATA(AOBA(I).191924)/1.0,O.935,O.875,0.82,0.765,0.710Op6,tO.615,
1 0.S7,0.525,0.480.0O.4490.4,0.358.0.3j2,0.282.0.24,O.I95,0.14875,
20el2897909107S9O9O83ld4,0.055l89o0.0/-
DATA(SLIPA(I) ,Iulg9)/0090.O5,0,).0.jO1.590.25~3o0.35gO.4/
DATA(CCA(I),Iul,9)/0.5,0.71,0.92,11lo23,1.3,1.35,1.38,l.4/
DATA(EEA(I) .131,9)/0.4,0.41,0.42,O.4255,0.4285,0.43,0.4390.43,0.43
I/ 1-) I,
120 READ (NR93) RXRX2sRY1,RY2,RMS1,RNS2
READ (NR,3) EloE2tVl.V2
READ (NR931 CONSICONS2,CPS1,CPS2,ROSRsOS2
READ(NR93) EPST9EPSRDsDTF9RDU
READ(NR*2) ITTt ITRDUoNRITENPTSNSTRS
IF(tdSTRS *EQ. 1) Go To 90
READ(NR92) JJFS.NIHE
NOPL,1
READ(NR93) ( YCJ) JuliJJFS)-
90 WRITE(NW9500)
I
500 FORMAT41hl)T
WRITE (NW96)
WRITE(NW94)
READ(NRol)
READ(NR93) WTvTEMI,1J29U1
READ(NR,3)VISALF.BETGAMCONFCPFRHOF
READ (NR93) PSTAR*GFC
READ(NR92)LINCONgNEXTtNGF
WRITE (NW91)
WRITE(NW*4)
WRITE (NuelO)
WRITE (NW911) RX1wRX2
IF(LINCOt4.EQ.1) 60 70 600
WRITE (NU912) RY19RY2
600 WRITE (NW913) RMS1,RMS2
WRITE (NW930)
WRITE (NW*14) E1,E2
WRITE (NW915) VI9V2
WRITE (NW,16) CONS1,CONS2
WRITE (NW.17) CPS19CPS2 -
WRITE (NW*31)
298
IF(LINCON .EQ- 1 ) WRI7E (NW,57) WT
7T IF(LINCON .EQ. 0 ) WRIIE (NW,19) WT
WRITE (NW,21) TEMI
WRITE (NW920) UlU2
WRITE (NW,32)
WRITE (NW,221 VIS
3
WRITE (NW92 ) ALF
WRITE (NW924) 83EI
WRITE (Nk*25) GAN
WRITE (NW926) CONF
WRITE CNW,27) CPF
WRITE (NW*28) RHOF
PS=200000.
TWF11I.
RXP= j.0/RX1+1.0/RX2
RYP= 1.0/RY1*1.0/RY2
0
TEMI=TEMI*46 .
AA=0.5*RXP
DDO0.5RYP
RXP = 1.0/RXP
RYP = .0/RYP
AMD=(AADD0)/(AADD0) 4
CALL TLU (AIDFMAMDA9FMA92 ) 24
CALL TLU (AMDAOBAMDAA0BA9 )
P1=3*1416
EP = Pi /((i.o-Vl**2)/El1-(1 O-V2**2)/E2)
33 3333 3
SB=FM*(0.75*P1*WT/ (AADD)/EP)**0.
wW 0,75*WT/S3
IF(LINCON .EQ. I ) WW=WT
EDOW= EP*2*0/PI
WB ARWW/EDOW/RXP
UBARVYIS (UlU2)*0e5/EDOW/RXP
IF(LINCON.EQeI) 60 70 100
IF(AMD-0.9) 101,1019100
100 SA= 2.0*(WW*RXP/EP)*0hS5
PHZE= (WW/EP/RXP) #*O*S/PI
PHZ=EPPZE
GO TO 105
101 SA= SB*AOB
PHZE= 1 5*WT/ (PI*SA*St5*EP)
105 CONTINUE
A1.0=S
SLIP=(U2-Ul)/U2
SRRAT=(U2-U1 )/(U2*U )*2.0 2
AM46o*VIS*(Ul*U2) *EP*HXP/WW**
PHZzEP*PHZE
IV(PHZ.LT * 1 ES) PSI AR=PHZ/2*
QMzVIS (Ul.U2)**2/ (CONF*TEM1 )*2.315
Flu (WW/EP) O*3/RXP
FZ=SQRT (FI)
Fla(FI/4*)**O.ZS
GI=PI*RHOSI*CPSI*CONSI
62xPI *RHQS2*CPS2*CO4S2
CONVIRHOF*CPF*(U14J2) .VZ/CONF*21600o
D(2)sCONF0.5SSQI.' t./(62*J2) )/Fl/208.
AFA sALF*PHZ*Pl/2oO
BTAZI3EI/1E41
GMAxGAMOPHZP/2o* /lEMI
GG a ALF*EDOW
299
BOAxI .0/AO8
IF(LINCON.EQ*1) BOAzlOe
CALL FTFS (FEET9FEES)
IF(LINCONi.EQ. 1) FEESx1.0
SHMITzRXP* GG*-0.6*U8AR'*0.7/WBAR*40*I3 *1.6
SHMz SHMIT *FEET*FEES
B02oO*WW/(EP*SHMI
CALL 'lLU(SLIPqCC*SLIPA9CCA99)
CALL TLU(SLIPtEE*SLIPA*EEA,9)
EBAR=CC*GO* .
EBARxEBAROIJBAR**EE
ESARxEBAR~wBAR*O (-0.74)
HMIlNxO&S
SPFILMS~HM/SQRT (RNS1**RMS2**2)p1000000.
PT~z1 ."HMIN
C WRITE OUTPUT DATA
WRITE .NN,45)
WRITE '(NW946)OM9CONV
WRITE INW94?)D(1)9D(2)
WRITE (NW935)
WRITE (NW936)
IF(LINCON*E~o.) GO TO 601
WRITE (tdb.37)SASB
WRITE CNW938)BOAtFM
GO TO 602
601 WRITE(NW960) SA
602 WRITE (P4W939)
WRITE (NW942)UBAR
WRITE- (tdW43)GGAFA
WRITE (NW944)BTA9GMA
WRITE (NW940)WWoWBAR
WRITE (NW*41)PHZoPhZE
WRI-TE(NW*55)8*SRRAT
* WRITE(NW956)AM
WRITE (NW949)
WRITE (NW948) FEETvFEES
WRITE (NW950) S#IM179SHiM
WRITE(NW,5l)
WRITE (NW,52)E@ANPTD
CALL AREA (EBARqHMItISPFILMvABAR)
ABAR.100O.*ABAR
WRITE (NW9S3) ABARoSPFILN
IF(NGF.EQa1) GO TO 110
CALL FRICT
110 IF(NGF.EQ.1) FCOEzGFC
WRITE (NW954) FCOE
IF(NPTS.EQ.1) GO 10 200
IF(EBAReGE.0.5) GO 70 603
CALL PRESS
CALL TEMPT(RDU)
IF(NSTRS *EQv 1) 60 TO 200
KKFSx27
DO 190 KzlKKFS
KKKzK,14
190 XST(K)xXH(KKK)
XSTT=-l .0
h
XDEX=O.25
YSTT=-0 .5
YDEYO .25
CALL STRESS
GO TO 200
300
*603 WRITE(NW*604)
200 IF(NEXT oEQ*'1) GO TO 120
CALL EXIT
1 FORMAT(72H
I
I 2 FORMAT(1715)
3 FORMAT(8E10.3)
4 FORMAT (/
5 FORMAT( 40X 45H4 PERFORMANCE OF AN ELASTOHYDRODYNAMIC CONTACT
6 FORMAT( 40X 5***********#********e* )
10 FORMAT(// 40X 1711 GEOMETkIAL DATA t//)
11 FORMAT(55H RAD, OF BODY 192 IN THE DIRECTION OF ROLLING,(INP
19 E12.5,2XE12e5/)
18 FORMAT(55H DENSITY OF BODY 1,2,(LB/IN**3)
It E12.592XE12.5/)
19 FORMAT(55H TOTAL LOAD, (LB)
1* E12.5/)
20 FORMAT(55H SURFACE VELOCITY OF BODY 192,(IN/SEC)
30 FORMAT(/// 40X 14H1 MATERIAL DATA 9///)
31 FORMAT(/// 40X 21H1 OPERATING'CONDITIONS 9///)
32 FORMAT(/// ISH LUBRICANT DATA 9///P
21 FORMAT(55H INLET TEMPERATURE OF THE LUBRICANI,(DEG. F)
1. E12*5/)
22 FORMAT(55H INLET VISCOSITY,(LB-SEC/IN**2)
It E12*5/)
23 FORMAT(SSH PRESSURE-VISCOITY COEF., (IN**2/LB)
24 FORMAT(55H TEMPERATURE-VISCOSITY COEF.,(DEG. R)
1. E12.5/)
25 FORMAT(55H PRESSURE-rEMPERATURE-VISCOSITY COEF. (IN*iDEG.R/LB)
I, -. 19 E12,5/)
26 FORMAT(55H THERMAL CONDUCTIVITY OF LUB.,(9/OEG.F-HR-FT)
19 E1295/)
27 FORMAT(55H SPECIFIC HEAT OF LU89, B/LB-DEG*F)
1t E12*5/)
28 FORMAT(551 DENSITY OF LU8.,(LB/IN**3)
35 FORMAT(/// 40X 221.CALCULATED PARAMETERS I
36 FORMAT W/ 40X 1811 CONTACT DIMENSION II
37 FORMAT(35H "ALF HERTZ WIDTH(ROLLING DIRe)9IN. 9X E12*598X
1 38HHALF HERTZ WIDTH(PERP. TO ROLLING)9IN* 2X E1295)
38 FORMAT(41 A/B 40X E12oS 98X IN1K 39X E12.5/)
39 FORMAT( /// 40X 3711 LOAD9SPEEO.AND LUBRICANT PARAMETERS /1//)
40 FORMAT(33H EQUIV# LINE CONTACT LOAD(LBIN) lIX E12*5 98X 4HWBAR
:1 136X E1295 I
301
41 FORMATC20H HERTZ PRESSURE(PSI) 24X E12,598X 6HPHZ/ED 34X E12.5/)
42 FORMAT(511 U8AR 39X E12,5 /)
43 FORMAT(3t 66 41X (12.59 8X 31IAFA 37X (12.5/)
44 FORMAT(41 BTA 40X E12.59 8X 3hGMA 37X E12.5/)
45 FORMAT(/// 40X 2011 THERMAL PARAMETERS II
46FORMAT( 3H ON 41X EIZ959 8X 4hCONV 36X E1295/)
47 FORMAT( 3H1 DI 41X (12.5. OX 2HD2 38X (12.5/)
48 FORMAT(25k THERMAL REDUCTI"ON FACTOR 19X E12.59 OX 19HSIDE tEAKAGE
IFACTOR 21X E12.5/)
49 FORMAT(/// 40X 1811 REDUCTION FACTORS II
50 FORMAT(21H IS0TH. MIN* FIINCIN) 23X (12.5 '.8K I3HMIN. FILM(IN)
127X E12.5 /)
51 FORMAT(/// 40X 3611 DIMENSIONLESS PROTRUSION PARAMETERS //I)
52 FORMAT0iTH PROTRUSION WIDTH V7K El2o5t 8X 16HPROTRUSION DEPTH
124X E12.5/)
53 FORMAT (27N PERCENTAGE OF AREA CONTACT 17X E12*5.8X l4hSPECIFIC
VILM 26X E12.5 /)
54 FORMAT(17H FRIC7IONAL COEF. 27X E12.5/)-
55 FORMAT(16H 8=290*Ww/EP/SHM 28X (12.5, 8X 36HSLIDE TO ROLL RAT!
10=2WU2-UI)/(U2*Ul) 4X (12.5)
56 FORMAT(31 AN=6.*V1S*(U1.U2)*EP*RXP/WW**2 13X (12.5)
57 FORMAT(55H LINE CONTACT LOAD*(LB/IN)
19,E12*5/)
60 FORMAT(35H HALF HERTZ WIDINh(OLLING OIR*)9IN* 9X (12.5)
604 FORMAT(/// 41H1 THE PROTRUSION WIDTH IS GREATER THAN 0.5 58H PR
1ESSUREtTEMPERATUREgANO STRESS CALCULATIONS ARE BYPASSED
END
SUBROUTINE TEMPT (RDU)
C START%
DIMENSIOk AD(2),ADC(2).SUMT(2),TT(2),AF(2,50).FUC2,50).PP(50),
1 DCC?)
COMMON KIKCKGKAKOKDKB.KEKFKRNRNW
COMMON 'RK,RX2-RY1,t2RMS1,RMS2,E1.E2,ViV2
COMMON CONS1,CONS2,CPS1,CPS2.NHOS1,RHOS2,WTTEMI9u2,UI
COMMON VISALFBEY.GAMCONFCPFRHOFPSPSTAR
COMMON AFA9BTA9GMA9A,89AMqSLIP* RXP.RYPEPOMCONvO(2) .COMP.AHZ,
1PHZ*WBARtUBARvPHZE*GGoEDO~W90A *WW
COMMON FCOE.EBARFEETFEESHMINSPFILM.ABAR.SHM
COMMON EPSTEPSRDOTFTWFIITIWUTTNRITE
COMMON Kt(50)HE(50),P(50),T(50),TI(2,50),F(2.50).SHS(50)
FUNA(KY)zAFA*K. (BTA.GM4A*X)*(1./Y-1.)
FUNB(KUtZ) z(U*(0.33333*U.Oo5*COMP*K).Z/6.O)
FUNC(KU.Z):u*(3.0*(U.COMP*K) *L)
FUND(X.YZ)zRDU*AFA*(X-Z).(81A.HDU*GMA*(K-Z))*(1./Y-1.)
ROUT zRDU
WRITE (NW9802)
U1BARU/(Ul*U2)
SLP=(SLIP/(2*0-SLIP))
SLPG=SLP*SLP
IF(GM-0*5 )77597759776 -
775 XINIn-SQRT(390/8*190)
TINI=1*0
G3=FUNB(P(Kl) .0.75.SLPO)
OELTTSEXP(FUNA(P(KI).1.0))*3*UM
GOTO 777
776 CALL THILT(Q)4,I8TACONVBU1BARKNITINb;NR1TE)
IF(TINI.EO.(1.)) 60 TO 775
DELTT=TINI-1*0
302
777 00 4 K*KI*KF
T(Kizl.0
TI (19K)=1.0
4 TJ(2,K)=1.0
Fj DO 755 KxIKA
TEMPzXINI-XHtK)
IF(TEMP)75097559755
750 KT=K
TEMP=DELIT / (4.O.X11(KT))
T (K)=TINI
T(K.1)=TINI(X-dKT.1I)-XH(KT) )*TEMP
A T(K.2)=TINI.(XH(KT.2)-XH(KT))*TEMP
T(K,3)=TINI, (XH(KJ.3)-XH(KT) )'IEMP
I755 GO TO 401
CONTINUE
401 FABBA=AM*BB/SORTF( A)
500,NKF =KF
I Fl =B/A**O.25
FI=SQR7F( (2.0-SLIP)/(I.0-SLIP))
F12=SQRTF (2.0-SLIP)
AD(I)=D(1)*Fl
AD (2)=0 (2)'Fl
CABB=CONV*SQRTF (A)/B/B
FF=EXPF(FUNA(P(KT)9,T(KT)))
FH=I .0-1 .0/HE (1T)
F10 =FUNC(P(KT)oFISLPO)
F(lKT)=AD(1)*(FIO-2.O*SLIPFI)FF/HE(KT)*QM
F(2,KT)=AD(2)*(FlO.2oO*SLIP'Fh)'FF/HE(KT) '014
KKT=KT* 1
IF (NRITE)88987988
87 WRITE(NW92)
88 CONTINUE
ITRDU~l
RDU=RDUT*0.95
6 Do 40 K=KKTgNKF
ii HlHE(K)
A21 KKK =K-KT.1
C CHECK KKK LESS OR GREATER THAN 3
IF (KKK-3) 10910.25
10 GO TO (11,11,20),KKK
11 SUMT(1)=T(KI)
*1 FCDET=CAeB*H1*FDEI /60
F2O=CABB*FDET
F23=(FH*3.0'(FH.COMP*P(K) ).SLPQ)*FF
F31=F23-2 .0*FF*SLP*FH
303
F32zF23*2o0*FF*SLP*FH
F(2.K)zAD(2)*(TI(1,K)-I1(2,K).F32*QM-F20)/HI
DO 50 Jzl*2
50 TT(J)xSUMT(J)SQRTF XH(K)-XII(K-1))0(1.3333333*F(J.K)*
1 0,66666666*F(JtK-1))
F7a(-4CTT(1)TT(2))*0*5*(QM *5-FCOET))/T(K))
IF(F7)760,760,190
760 T(K)=(T(K)-1.0)*0*95.1*0
TM2*T (K)
GO TO 763
190 FC2zkOGF(F7)
IF -IT-1) 20092009205
200 DTmDTF*(T(K)-T(K-)-)
TMI zT(K)
T (K)aT (K) .DT
GO TO 220
205 DFCx(FC2-FCI)/(TM2-TMI)
DTK=FC2/OFC*TWF 1
T141 xT(K)
T(K)xT(K)-DTK
IF (ABSF(OTK)-EPST) 210921092r5
215 IF (IT-ITT) 2209720.8tb
220 FC~xFC2
TM2xT(K)
TI(19K)alTT(1)
TI (2,I%)xT(2)
763 IF(NRIJE)76697649766
764 IF(ITROU-1)7669765*766
765 WRITE(t4W994) FI1FCOET*F(29K)9F(19K)9G59F7
WRITE(NW,3)XH(K),T(K),IT71(1,K),TI(2,K),FC2
766 IF(IT-ITT).770.770985
770 IT aIT+I
GO TO 12
210 IF(NRITE ) 81980.81
80 WRITE(NW93) XI(K)9T(K)9IT9TI(IK)9TI(2#K)
81 GO TO 40
85 WRITE(NW*86)
NERRal
G0O7010Q*
20 Fl 1wSORTF'(XH(K)-XH(K-M)
F2 xSQRTF(XH(K)-X1(K-2))
DO 21 14.1,2
21 SUMT(M)s*T(KI)*0*5*(F(M9K-1) /F1.F(M#K-2)I/F0
1 (XH(K-1)-XI1(K-2))
KEOx2
GO TOJ 15
25 6O TO (35930)tKEO
30 KKFuK-3
KEOxI
GO TO 33
35 KKFzK-4'
33 KKzK-1
DO 31 M=1*2
DO 31 J=KltKK
31 FO(MJ)xF (MqJ)/SQRTF(XH(K)-XH(J))
DO 32 M1012
SUMT(14)xTIKI)
0O 32 JzKltKKF92
304
32-SUMT(M)=SUMT(M)+(FD(M.J) *4.0*FD(MtJ*1) I+FD(MtJ.2)
I 1U )*(XH(J,1)-X1(J)) /3.0
GO TO (609150)oKEO
150 DO 61 M=192
fT61 SUMT(M)=SUT(M)05*FDMIK'1) +FD(MK-2) )*(XH(
1 K-1)-XH(K-2))
60 T (K)=(T (K-1)-T(K-2) )/(XH(K-1)-XH(K-2) )*(XH(K)-XH(K-1) ) T (K-I)
DO 160 M=192
160 TJ(MtK)=rI(MgK-1)
GO TO 15
40 CONTINUE
F3=(VIS*(U2eUl)/SIM)/(WW/AHZ)*SLP
DO 720 K=KIKF
* - IF(P(K)-P5,TB )11097109711
710 F6=EXP(FUNA(P(K),1 (K))
6O TO 712
711 F6=EXP(AFA*PSTO.FUND.(P(K),T(K)'PSTB))
712 SHS(K)=F6/iE(K)*F3
720 CONTINUE
FRC=0.0
KKF=KF-2
DO 725 K=KI,'KKF*2
725 FRC=FRC,(SHS(K),4.*SHS(K,1,,SHS(K,2))*(XH(K,1)-XH(K))/3.0
FRC2=FHC-FCOE,
IF(ITROU-1)660*660,661
660 RDUl=RDU
ROU=RDUT*1*05
6040O 730
661 DFRC=(FRC2-FRC )/ (RDU2-RDU1)
IF(ABSF(FRC2)-EPSRD) 100,100,726
726 DRDU=FRC2/DFRC
IDUI=RDU
-a RDU=RDU-DRDU
IF(ITRDU-ITT) 73097309740
730 FRCI=FRC2
RDU2=RDU
ERFR=FI1C2
WRITE(NW995) ITRDU-PRDUoERF14
SITROU=ITRDU41
-. GO TO 6
740 WRITE(NW,96)
GO TO 100
100 CONTINUE
600 CONTINUE
WRITE (NW9801)
DO 800 K=KIKF
800 WRITE(NW,94) XH(K),P(KhoT(K),TI(1,K),TI(2,K),SHS(K)
RETURN
1 FORMAT (1HO,6E12*5)
2 FORMAT (53H1 XII 7 IT TI T2
3 FORMAT (IH 9F12979E12*591493E12*5)
8 FORMAT W6 12.1)
72 FORMAT (20H1 EXP ARG TOO BIG 9I49E12*5)
86 FORMAT(20H DIVERGE TEMP)
-*90 FORMAT(F12*7914)
93 FORMAT (lIH 91714)
94 FORMAT (IH 96F12&7)
*195 FORMAT(/7H ITROU=9 15,6X,5HRDU= ,E12*5,6X, 5HERFR=tE12.5/)
96 FORMAT(20H DIVERGE RDU)
801 FORMAT(7X 4HXBAR 6X 81IPRESSURE 5X 6HT MEAN 8A 2HT1 1OX 2HT2 4X liii
305
ISHEAR FORCE) --
802 FORMAT(///)
END
SUBROUTINE THILT(QM#BETAOCB 9uIAXXTT9NRITE)
DIMENSION TC(200),A(200),VT(9),VIST(9 )9VS0(9),HLZOO)'
Fl (ARG)=-ARG*SORT (ARG*2-1.)-ALOG(-ARG*SQRT(ARG*02-1.))
SPil .-2.*IUl
/
DXz0.05
HST=6.
TIzI.
IF(XS.GE.(-I.)) XS=-l.-o00I
11 XS=XS-((IlSP-1.)/B-FZ(XS))/(2.*SQRT(XS**2-l.))
IF(XS.GE.(-1.)) XS=-l.-0.001,
MP=V,*B*F1 (XS)
IF (ABS(HP-HSP).GT*(O.1)) GO TO 11
XSsXS.0.S10./
IF(XSP.GE.(1I.)) XSP=-1*05
- A~Xz-S0RT((t1ST-I* /H,1.0)
IF(XO*GE*(-I*)) XOz-1.-O.001
10 XOX-(S-o/-IX)/2*QTk*2l)
IF(XOGE.(-l.)) XOx-I.-0*01
HSxl.*B*FI (XO)
IF 'A8S(IS-HST)oGT.tO.1)) 60 TO 10
X~sAINT (X0*10*)/10.
H(1)sl*#B*FI (XO)
IF(NRITE-.EQ.l) GO TO 90
WRITE46980). XOH(1)
90 IF(XO*GE.(-1.l )IGO TO 65-
TC(l)zrl
X(1)=XO
MzO
N=O
K=0
J=1.
H(J)=l(J)l
J PH=12 -1/H(J)
Al:-2.*PI,/5...CP/3.*2.*Ul/3.
A2='-Pt/IO.,SP/6..IJ/3.
IF(M.E(J.11 GO TO 40
IF(IK.EQ* 1) GO TO 60
22 AASxAA -
25 YxO*
SUMSO.
DO 30 K=199
VIST(K)=EXP(BETA*vT(,)
VSG(K)36.*(PH*Y)*2,l2.*PHIY*(S-3.*PH),(SP-3.*PI)*
SUM=SUM.VIST lK)*VS U(K)
30 Y=Y.0.125
TGLz:(SUM-VIST(1)*tVSO(I)/2.-VIST(91*VSO(9)/2.)*0.l2S,
11=B*e2.OM/QC.IGL/H(J)/2.
IF(N.EQ.l) G0 TO 50
AA=(8J-82)*6X/2.*AA
.X(J)=X(J)tDX/2,,
60O TO 2G
306
[1 40 TC(J)=(AA/H(J)-A2)/Al
N= 1
GO TO 25
50 N=O:-
-AA=(B1-2)*DX.AAS
6O TO 70
65XX-1.J
TT=1.
~ :"; 70 CONTINUE
6IF(NRITE.EQ.1)l GO
WRITE(6980) XXTT
TO 95
80 FORMAT(4EI4.4)
95 CONTINUE
RETYJ'N
END RC
*DIMEJSIO G(3tG()A37iFT(13?7tZ7*?tZ-.
COWO KIKC ,KGiftAKO.KKBKIEKF.KIt'R.NW
COMMON RX1,RX2,RY1,RY2,RMS1,RMS2,E1,Eaev1,V2
COMMON CON4S1lCONSZ-,CPS1,CPS2RiOS1,.RHQS2,WT,1EM19u29,u.
~'COMMON VISALftBET.GAMCONFCPFIHOFPSVPSTAR.
.4 COMMON AFA98rA'9iiMAA,~,AM9SLIP# RXPRYPEPQMCONVD(2Y~,COM.,AHZ,
.4 1PHZqWBAR.UBARPHZE.GG*EDOWB0A 9WW
COMMON FCOE.EBARFEETFEES.HMINSPFILMABARSHM
COMMON EPSTEPSRDDT- ,.TWFI1TRDUITT.NRITE
COMMON XH(50),HE(50)-,i'50),T(50),TI(250)F250)SiS(56)
>1 1 0.01359
f)ATAFTN()91=113)/*000290*0.79 OVItv.0*00239 0.006?9
0.021, 0;9379 0.046, 0.0529 0.0615, 0o0685q'0.;0775/
DiATA(FTN(1),I=1 4,26 /0.9000,79'0,0029 0.00399 0.00659 0.01500
307-
DATA(FrN(I)91=14491S6)/ 0.0129 0.02259 0.029. 0.034b, 0.043,
1 0.0489 0.0529 0.0599 0.0639 0.0669 0.070. 090739 0.077/
- - I 368
DATA(FTP4(I),I=5349546)/0*00719 0.0159 0.0239 0.03, 0.049
I 0.045. 0.048, 0,051. 0.0529 0.053. 0.0549 0.055. 0.055/
DATA(FTN(1).1=5479559)/o00003,.00009,.00018,.0003..0009.
1 0.00199 0.0039 0.009. 0.016. 0.023. 0.037. 0,045. 0.051/
- DATA(FTN(1) .I=560.572)/,00019.00026..00058..0009..00299
1 0o0057. 0.00929 0.021. 0.0299 0.0349 0.041. 0.041. 0.051/
DATA(FTN(I).1=5739585)/.000259.00074.0.0015. 0,00249 0,00129
1 0.013. 0o018. 0.028. 0,035. 0.039. 0*045. 0.048. 0,051/
DATA,(FTNII)tI=5869598)/0.000549 0.0016. 0o00339 0,0052. 0.014,
1 0*0229 0.02659 0.0359 .04. 0.042. 0.04659 0.049. 0.051/
DATA(FTN(I)91=599q611)/0o0013v 0.00369 0.0069. 0.011. 0.023.
1 0.029. 0.03259 0.03959 0.0449 0.0469 0.0485. 0.050. 0.051/
4 DATA(FTN(I),I=612.624)/0.0015. 0.041. 0.008, 0.0129 0.02359
1 0.0294 Bl-011. n~A 4-1h 0.04659 0.0499 0.05059 0.051/
DATA(FTN(I).1=625q637)/0o00189 0o00529 0.00949 0.014. 0.02559
B -- 1 0.031. 0.0359 0.042. 0.0459 0.047. 0.05. 0.051. 0.051/
DATA(AG1(I),I=1913)/1.E-89 3oE-8. 6.E-89 I.E-7. 3.E-79 6.E-7.
I 1.E-69 3-E-69 6.E-69 I.E-S. 3-E-59 1.E-49 1.E-3/
DATA(A62(1).I=197)/5.E-79 1.E-6. 5oE-69 I.E-5. 5.E-59 SoE-49 1.E-3
f
WRITE (NW920)
j TEPP=TEMI/546.-1.
VIS3=VIS*EXP (BTA*TEMP)
ALF3=ALF+GA#4*TE1MP/TEM i
OM3=QM/V IS*VIS3*T EM /546.
FEETT=FEET
FEEST=FEES
QMT=QM
QM=SM3
BTAT=BTA
BTA=BTA*TEM 1/546.
f AFA=AFA/ALF*ALF3
CALL FTFS(FEETgFEES)
IF(NRITE 5.0.51 O 0
50WRITE(NW931) 61. 62, 3, VIS3 ETFE~VI39M
30 FORMAT (8E411G.) .. 4 2,1..1 3E2561VSZ 25
501 IF(G1.EQ.0.S3/VIO)TOo*AF
245*Oo/ETTFE
1IF(62.EO..) GO TO 245
62=10VI*U-l*2BT(4.EM)8*CN*7/30
IGAGI(3))40.0.
IF(NRTE)
595309
60 UZ(I.J)=FTN(13*ItJ)
GO TO 100
70 UZCI.j)=EXP(ALOG(FTN(1.1.J) )*(ALUG(FTN(2,IJ) )-ALOG(FTN(1.I.J) ))/(
IALOG(AGI(2))-ALOG(AGI(1)))*(ALOG(Gl)-ALOG(AGI(l))))
100 CONTINUE
DO 200 K=197
IF((;3-AG3(7)) 1409160#160
140 IF(G3.LT.AG3(l)) GO TO 170
CALL TLu(G3,VZ(K)9AG3,LeZ(jvK)v7)
GO TO 200
(160 VZ(K)=UZ(7*K)
GO TO 200
170 VZ(K)*UZ(1.K)-(UZ(2.K)-UZ(1.K) )/(AG3(2)-AG3(l) )*(AG3(l)-G3)
IF (VZ(Ku.LEsO.0) VZcI~z1.E-5
200 CONTINUE
IF(G2-AG2(f)) 2049206.206
20'. JF(62*LT.AG2(l)) GO TO 207
CALL TLUG(62#Fik3tAG2*VZ*7)
GO TO 210
206 FR3=EXP(ALOG(VZ(b)).(ALOG(VZ(7))-ALOG(VZ(6)))/(ALOG(A62(7))-ALOG(A
162(6 ))(ALOG(62)-ALOG(AG2(6))))
GO TO 210
207 FR3=IXPALOGVZ(2))(ALOG(VZ())-ALOG(VZ(2)))1v(ALOG(AG2(l))-ALOG(A
IG2(I))*ALOGG2)AL~bAG2()*)
210 IF(NRITE) 20b.205924.0
205 wRITE(Nw*2) FP3
2 FORP4AT(5t' FR3=oE12.5)
5 FORMATC7F9S.)
WRITE(NW95) UZ
WRITE(NW95) VZ
240 FCOE=ALO(G(FR3)-0.149*(ALOG(THEIA)-ALOG(30.))
FCOE =E xi A COE)
GO TO 250
245 FCOF.=0.
250 CONTINUE
RETURN
END
SUBRWOUTINE TLU(Ao89C.IJ.N)
C LINEAR INTEkPOLATIUN ROUTINE
C A= INDEPENDENT VAkIAHLE
C d= DEPENGENT VAPIAHLE (ANSWER)
C C= INflEPENULENT TAHLf
C U= DEPENDUENT TABLE
C N= NUMBEW OF POINTS IN TAtBLE
c INDEPENT TAB3LE MUST 8E SORTED. EITHER ASCENDING OR L)ESCLNDING
DIM4ENSION~ C(iJ.DC1)
IFIN-1) 1.293
1 8=0.
GO TO 100
2 H=D(l)
Go TO 100
3 1L1,
MU=N
8 IF(MU-MLI) /I5*i59
9 M=(MU*ML)/2
IF (C(l1)-C(2))I11.2.10
10 IF(C(M)-A)13912914
11 IF(A-C(PMfl13912,14
12 B=D(M)
GO TO 100
13 MU=M
600TO8
310
14 MLzM
60 TO 8
15 BzD(ML).(D(NU)-D(ML))*( (A-C(ML) )/(C(MU)-C(ML)))
100 RETURN
END
SUBROUTINE TLUG(A98,C9DN)
C I LOGRITHMIC INTERPOLATION
C A= INDEPENDENT VARIABLE
C B= DEPENDENT VARIABLE (ANSWER)
C C= INDEPENDENT TABLE
C D= DEPENDENT TABLE
C N= NUMBER OF POINTS IN TABLE
C INDEPENT TABLE MUST BE SORTED9, EITHER ASCENDING OR DESCENDING
DIMENSION C(1)90(1)
IF(N-1)19293
1 8=0.
GO TO 100
2 B=D(1)
GO TO 100
3 ML=l
MU=N
8 IF(MU-ML-1) 1591599
9 M=(MU.ML)/2
IF(C(l)-C(2) )1192910
10 IF(C(M)-A)13912914
11 IF(A-C(M))13912914
12 8=D(M)
GO TO 100
13 MU=M
GO TO 8
14 ML=M 4
GO TO 8
15 B=EXP(ALOG(D(ML)),(ALOG(D(MU))-ALOG(D(ML)))*((ALOG(A)-ALOG(C(NL)))
1 /(ALOG(C(MU))-ALOG(C(ML)))))
I 100 RETURN
END
SUBROUTINE PRESS
COMMON KIKCKGKAKOKDKBKEKFKRNRNW
COMMON RX1,RX2,RY1,'RY2,RMS1,RMS2,ElE2,V1,V2
COMMON CONS1,CONS2,CPS1,CPS2,RHOS1,RHOS2,WTvIEMlU2,UI
COMMON VISALFBETGAMCONFCPFRHOFPSPSTAR
COMMON AFABTA9GMA*A9H9AMSLIP9 RXPRYPEPOMCONV,D0(2) .COMPAHZ,
IPHZtWBARUBARtPHZE,669EOOWBOA ,WW
COMMON FCOEEBARFEETFEES,-HMINSPFILMABARSHM
COMMON EPSTEPSROD7DFTWFI1TRDUITTNRITE
COMMON XH(50),HE(50),P(50)iTI50),TI(2950),F(2,50ISHS(50)
DIMENSION SUMA150)
C ESTABLISH X-H(K),
KC=9._
KG=15
KA=17
KO=29
KD=33
KB=37
KF=41
KR=41
XH(1)=-4.0
DX=0.5
KKI'=Kl,1
U 100 5b N=KK1,KK4,
311
50 XH(K)sxh(K-1) *UX
P(KA)=5.O/AFA
XH(KA)=-SQRT (1.0-(3.1416*P(KA)*0.5)**2)
DX=-XH(KA)/4.0
XH(6)=XH(53 .DX
XH(7)=XH(6) .DX
XH(8)=XP.'7) .DX*0.5
XH 19) zXh (8) DXO .5
00 55 K=10913
55 XH(K)zXXH(K-1J.OXOO.25
DXz(1.0.XH(KA) 300.S
KKAzKA-3
DO 60 K=KKAKA
60 XH(K)=XH(K-1)*DX
DXs-XH (KA),8.
KKA:KA.1
KA4 KA.4
D0 65'K=KKAtKA4
65 XH(K)BXH(K-1)*DX*0.25
XH(KA.5)zXHIKA*4) .DX*0.5
XH(KA*6)xXH(KA+5) .DX*0.5
KA7=KA. 7
DO 66 KSKA7,KO
66 XH(K)*XH(K-1)+OX
DXsO.25*(1.0-2oOOEBAR)
KKO=KO. I
.D0 70-RzKKOKD
70 XH(K)*XHCK-1)+0X'
DXu0.25*EBAR
KKD=KD. 1
KkKF=KF-I1-
DO 75 K=KKDKKKF
?5 XH(K)nXH(K-l)+DX
XII KF) .0.9999999999
IF(NRITE)779?697?
76 WRITE(tNw,2)(AH(K)#KuI.KF) -
C ESTABLISH HEMK
77 KKAuKA-l
DEXHUXH (KA-2) -XHI(KA)
DO 80 Kxl.KKA
XHKK=Xi4(K)*DEXH
TEMPuSQRT (XHKK 0*2-1.0)
80.HE(K)*I.0.B* (ABS(XHKK *1EMP)-ALOG(A8S(ABS(XHKl )+TEMP)))
DO 85 K=KAKB
85 HE(K)=1.0
KK92KB+ I
* ' DO 90 K=KKB*KF
90 HE(K)=1.0-(1.0-H4IN)*SIN(3.1416*(XH(K.-XHCKB) )/(l.0-X&1(KB)))
IF (NRITE)92991992
91 WRITE(NW92) (HE(K)9K=1,KF)
C ESTABLISH P(K)
* .92 SUMA(1)=O0A
DO 100 Kx29KA
FK1x(HE (K-I )-1.0)/hE (K-I) **3
FK2=(HE(K)-.)/4E(K)**3
100 SUMA(K)=SUMA(K-1).(FKI.VK2J*0,5*(XH(K)-XH(K-1))
P(1)=0.
DO 105 K=29K.
P(K)=AL06l(.0-SUA(K)/SUMA(KA)*(.O-EXP(-AFA*P(KA))))/AFA
-- 312
105 P (K)=-P (K)
P12=0*63662
KKA=KA. 1
DO0 110 K=KKAKD
ii110 P(K)=SQRT(1.0-XH(K)**2)*PI2
P(KD)=PI2*SGRT(1.0-(IoO-2.0*EBAR)**2)
PBB=PI2*SQRT(1*0-(I*O-EBAR)**2)
KKD=KD'1
KKB=KB.1
DO 115 K=KKBoKF
-Fl=-PB8lv((1.0XH(K-1))/EBAR)**2,PI2*SQRT(1.0-XH(K1I)**2)
F2=-PBB*((1.0-XI1(K))/EBAR)**2,PI2*SQRT(1.0-XH(K)**2)
115 WS=W5.(Fl.F2)*05*(XH(K)-X1(,K-1))
SUM4=0.0
SUM5=0*0
DO 120 K=KKDKB
Fl=( (XH(K-1)-XI4(KD) )/EBAR)**2
F2=( (XH(K)-XH(KD) )/EBAR)**2
F3=P(KD)*(1.0-Fl)-PI2*SQRT(1.0-XH(K-1)**2)
F4=P(KD)*(1.O-F2)-PI2*SQRT (1.0-XH(K)**2)
SUM4--SUM4,(F1,F2)o.*5*(XH(K)-XH(K-1))
120 SUM5-SUM54(F3,F4)0.5*(XH(K)XI(K-1))
l(KB)=(W5-SUM5)/SUM4
Do 125 K=KKDtKB
125 P(K)=(P(KR)-P(KD))*((XH(K)-XH(KD))/EBAR)**2.P(KD)
F ~DO 126 K=KKBKF
126 P(K)=PBB*((1.0-XI1(K))/EBAR)**2
WRITE (NW,1)
14 1 1 FORMAT(I0X 4HXBAR 6X 12HNORMLZD FILM 5X 8HPRESSURE
DO 130 K=KIKF
130 WRITE(NW,2)XH(K),HE(K),P(K)
-*2 FORMAT(3X 6(3X E12.5))
RETURN
END
SUBROUTINE STRESS
C STRESS CALCULATION FOR EHD (MODIFIED AND COMPILED ON SCI-TEK)
C REVISED FUR TIME DEPENDENT RESULTS, 2-14-69 9PN452 9 REEL NOo//13
C BASED ON ORIGINAL VERSION~ OF 4-9-64
C NSW DIAGNOSTIC OUTPUT CONTROL, NSW= 1 DIAGNOSTIC OUTPUT
DIMENSION SX(31),SY(31).TXY(31),PSX(31),PSY(31),PTXY(31)
DIMENSION AP(31),AQ(31),BP(31),BO(31),CP(31),CO(31)
DIMENSION ZU2(31)gZUO(31),OS(31)tPSZ(31)9 RTXY(50931)
DIMENSION I8UF(1000)9XARRY(31),YARRY(3I) 9Q0M)
COMMON KIKCKGKAKO.KDKBKEKFKRNRNW
COMMON RX1,RX2,RY1.RY2,RM4S1,RMS2,E1,E2,VlV2
COMMON CONS1,CONS2,CPS1,CPS2,RHOS1,HHOS2,WTTEMl.U2,U1
COMMON VISALFBETGAMCONFCPFRHOF9PSPSTAR
COMMON AFA9BTA*GMA9AB9AMqSLIPv RXPRYPE ,QM.CONVD(2),COMPAH'Zt
1PHZWBARUBAR9PI$ZEGGEDOW*BOA ,WW
COMMOh! FCOEEBARFEETFEESHMINSPFILMABARSHM
COMMON EPSTEPSRDDTFTWFIITRDUITTNRITE
COMMON Xh(50),HE(50),P(50).T(50),TI(2,50),F(250.,SHS(50)
COMMON XSTTXDEXYSTTYOEYXST(31),Y(31) ,JJFSKKFS#NTHEtNOPL
A WRITE (6,700)
JJF=JJFS
K KF=KKF S
~1 ANU=V1-
DO 20 K=KIKF
20 Q(K)=SHS(K)
NSW=2
520 FORMAT (75H-STARTING VALUE OF X INTERVAL OF X -STARTING VALUE 0
IF Y INTERVAL OF Y )
313
522 FORMAT (4XE12.5.,7XE12.58XE12.5,7XE12.5)
152 ANTZNTHE-1 -
DELAz3. 14159265/ANT
DE~i9oo/ANT
KFPI=KKF.1
KFP2zKFP 1.1
WRITE (6.512)
WRITE(692) (XST(K)tK=1,KKF)
WRITE (6.701)
WRITE (69514)
WRITE '-(692) (Y(J)q Jul, JJF)
WRITE (69701)
00 300 Lxl*JJF
YJsY(L) -
Y2uYJ*2
WRITE (698) YJ
XH7=XST (J)
UKIX7 -XH(KI)
UKF=Xk7 -XH(KF)
YKI2uY2*UKI**2 0064
YKF2xY2+UKF**2 0065
YP.KlaY2/YK12 0066-
Y2KFxY2/YKF2 0067
YUKI*YJ*L5KI/YK12 0068
YUKF=YJ*UKF/YKF2 0069
SS~P(I *UK-P(F)YUF-QKl*Y2I*(KF*YKF0070 -
SXYxP(Kl)@Y2KI-P(KF)*Y2KFQ(Kl)*YUK1-Q(KF4)*YUKF 0071
SSY.-SSX 0072
GO TO (58959)tNSW
58 WRITE(692) SSXoSXYtSSY
59 N2=KF-2005
DO 350 IaKA9N2,2 07
IIw(I1/ *007
IF(L-1) 3569355,356 --
Z~sZU(11)0082
ZOauO(II) 0083
Z13XH7 -XH(1*1)
Po=P(I) -~0085
*PlwP(I+.1 0086
P2xP(1#2) 0087
50 WRITE(692);ovZI9Z2vP0*PI9P2,HSQ
5NPOu1 0090
357 QQAz(Z2*(PO*Z1-2.64P1ltZ0).P2*Z0*Zl)*HSQ. 0091
00OQxz(PO(Z1,Z2)-4.0*P1*Z1,P2*(ZOZ1))*HSO 0092
QQC2(P0-C~o0'P1.P2)*HSQ 0093
GO TO (52953)ONSW
52 WRITE1692)QQA*QQBPQQC
53 GO TO (3589359)9 NPO 0096
358 AP(II)00QA 0097
OP(II)SQQB 0098
CP(II)xQgc 0099
PozO(I) 0100
PlzO(-I.1 0101
P2zQ(1.2) 0102 -
NPOx2 0103
Go To 357 0104
314
359 AG(11)=QQA 0105
'1 0QD=)QQ0 0106
CO(I1)=QQC 0107
356 AI=AP(II) 0108,
CI=CP(JI) 0109
CY2=CI*Y2 0110
AP3Y=A1-3*0*CY2 0112
AP2Y=AI-2*0*CY2 01
APY=AI.CY2 0113
AI=A(11)0114
CI=CQ(II) 01
CQY2=CI*Y2 0116
A03Y=AI-3.0*COY2 01
A02Y=AI-2.0*COY2 0118
AQY=AI+CQVw2 0119
Z2=ZU2(11) 0121
7O=ZUO(11) 02
I FLOG=ALOG( (Y2,Z2**2)/(Y2,ZO**2))
FTANATANi (Z2/YJ)-AIAN (ZO/YJ)
CP1=CP(JI) 0124
CQI=Co(I1) 0125
BPI=BP(I1) 0126
801=80(11) 0127
UY= (Z2-ZO) *YJ 0128
u3Y=3 * 0uY 0129
UU2=2.O (Z2-ZO) 0130
U2Y2=Z2**2-L0**2 0131
YAN=YJ*F TAN 0132
YV=JFO 0133
aOG=YJ*YLOG 0134
GO TO (55956)oNSW
55 WRITE (6,2)AP3Y.AP2YAPYAG3YA,2YAOYU3YUU2,UZY2,YTANYLOGY2LOG 01
56 SSX=SSX-CPI*U3Y-AP3Y*FTAN-PI*YLOG(HQIUU2CIU2Y2-3.08QIYTAN 0137
1,AQ2Y.*FLOG) 0138
SX=X-P*TNCIYLGCIUYA3*TN81YO 0139
350 SS'=SSYCP1*UY-APY*FTAN-B0I#YTAN-CI*Y2.OG -0140
SX(J)=-SSX*0.3175 0141
SY(J)=-SSY*0.3l75 0142
TXY(J)=-SXY*0.3175 0143
SHEA(0.5iSX(J)-Y(J)))2+TXY(J)2 0144
315
fill PTXY(J)uSQRT (SHEA)
53x(SX(J)*SY(J )Qe5 -0146
PSA(J)sPTXY(J)*S3' 0147
PSY(J)aS3-PTXY(j) 0148
PSZ(J)*(PSX(J)*PSY(." u*ANU
OS(J)=SQRT ((PSX(j) -ISY(J)) 0**(PSX(J)-PSZ(J))**2.(PSY(J)-PSZ(J)) 0150
10*2)/3*0 0151
WRITE (6910)
WRITE (6,12) XST (J),SX(J) SY (J),IXY (J).PSX (J),PSY(J) ,PSZ(U) ,PTXY (J),
105(J)
ANG2=Oo
D0 410~ MMUINTHE
S199SIN(ANG2)
* CO9wCOS(ANG2)
RTXY(JMM)s0.5*(SX(J)-SYCJ) )*SI4.TXY4J)*CO9
10 ANG2*ANG2*DELA
* ~ RITE(69518)DEG
WRITE' (69.2) (RTXY(J9MM),MMzINTHE)
WRITE (6914) -
400 CONTINUE
ANGL*O.
GO TO (3009602)*NOPL
602 NKTz1
603 D0 600 T1.KKF
IF(NKT-1) 604.6049600
604 XARRY(I)xXST(I)
600 YARRY(I)*RTXY(I9NKT)
XARRY(KFP1)aXSTT
XARRY (KFP2) aXDEX
YARRY(KFP1)*YSTT
YARRY (KFP2) UYOEY
IF(NKToNE*I). GO TO 530
CALL'PLO7S(IUuF*100.2)
CALL PLOT.(0*090959-3)
CALL AXIS(090,090,IOHUISTANCE-X,-10,5.0,0.0,XARRY(KFP1) ,AARRY(KFP2
2 FORMAT( 6E14*5)
4 FORe4AT(15IS)
8 FORMAT (5X2HY=*IXE12*5//)
10 FORMAT (5A3HXST.8XZHSX.8X2HSYSX3HSXY,7X3HPSX,7X3HPSY, fl,3HPSZ,6X4HP
1TXY#7X2H0S)
316
12 FORMAT(12(IXF9.5))
14 FORMAT (IN//)
500 FORMAT(14H POISSON RATIO)
Z
502 FORMAT (3511 KF JJF NC KKP NTHE NGPL
504 FORMAT(34H DIMENSIONLESS NORMAL LOADING QK)
506 FORMAT38H DIMENSIONLESS TANGENTIAL LOADING P(K))
508 FORMAT(37H DIMENSIONLESS WIDTH XH(K)
512 FORMAT(37H ARRAY OF X-COORDs FOR STRESS CALCU.
514 FORMAT(37H ARRAY OF Y-COORD* FOR STRESS CALCU.
518 FORMAT(23H SHEAR STRESS AT EVERY9F1O.5,8H DEGREES)
700 FORMAT(///)
'71FORMAT(/).
\SUBROUTINE
END
AREA(EBARHMINSPFILMA8AR)
DIMENSION HH(5I)vFH(51)
II KS=9
,Cl=l./SQRT(2o)
IF(SPFILM) 10911911
10 FHO=0.5.0.5*ERROR(ABS(SPFILM)*Cl)
GO TO 12
11 FHO=0.5-0.5*ERROR(ABS(SPFILM)*Cl)
12 ABAR=0.50 (2.'-EHAR) *FHO
FKS=KS
DO 20 K=1,KS
TEMPO0.5*ERROR(ABS(HH(K) )*Cl)
IF(HH(K) )14415915
14 FI(K)=0*5*TEMP
GO TO 20
15 FH(K)=095-TEMP
20 -CONTINUE
DX3=(EBAR/(FKS-1*))/3.0
KKS=KS-2
DO 3.0 K=19KKS92
30 ABAR=ABARDX3*0.5*(FH(K),4.0*FH(K,1),FHiK.2))
RETURN
END
FUNCTION ERROR (X)
E=1 .0/(1.040*3275911*X)
S=(((((0.940646070*E)-1.287822453)*E.1.25969513)*E-0.2521'28668)
1*E+0*225836846)*E
317.
LIATA /(ALPHA(I).I=1,3)/15.71923.564,312419/
DATA CFAB(ljI=1911) /1.0..99..98..96..939.899.8,.129.629.479
1 .36/
DATA (FAB(I)9I=12922) /1.0..989.969.949,.986..78..65..539.379
2 .27/
DATA (FAS(I),i=23933) /1...99..98,.939.89..81..66,.56,.44,.3,
3 .221
DATA (FAS(I).=34944) /l.,.99,.98.95,.9..86..77,.679.57,.449
4 .34/
DATA (FAB(1)9I=45955) /1...979.95..929.87..81..7,.6..479.329.22/
DATA (FAB(1)9I=56966) /I...979.94,99.849.76..62,5,369.24,.17/
Zx DATA (FAB(I)oI=67977) /1...98..96..93,9..84..759.66..569.439
5.934/
DATA (FA8(I)tI=78,88) /1.,.98.9S95 ,.985,.78,.66,569.459.319
6 .21/
DATA (FAB(I)91=89999) /1...98,.949.89,.82,.74..59,.46,.35,.23,
7 .15/
DATA(82(1)9I=195) /0.,0o,0.,47*5*107o8/
DATA(82(I)91=69l0) /0.,3*6912.8,48.8,75o6/
DATA(B2(1) .1=11915) /0.93.95912.4940e79.5I.4/
DATA(QM8(I),1=1.5)/0.0.0.5.1.0.5.0.10./
DO 100 I113
DO 100 J=193
IF(GN.LEoAQM(11)) GO T0 20
U(I.J)=EXP(ALOG(FAB(11,ItJ) )-(ALOG(FAB(I0,1,J) )-ALOG(FAS(I1.1.J)))
1/(ALOG(AQM(11)?-ALOG(AOM(10)))*(ALOG(QM)-ALOG(AQM(Il))))
6O TO 100
20, CALL TLU(ON.U(I.J) ,AQM.FAB(1,I.J),I1)
100 CONTINUE
AVl~sAFA*100000./P1Z
ALBTzBTA/AFm
IF(NRITEoNE.0) GO To 101
50 FORMAT(8E12*4)
WRITE(NW950) ON.AFI4.ALBT
101 DO 200 KxI,3
CALL TlU(ALBTV(K).BETA.U(I1.K)93)
200 CONTINUE
CALL TLU(AFMXFABoALPHAoV,3)
Do 300 11=193
DO 300 JJ=193
IF(ON.LE.014B(5)) GO TO 30
W(1IJJ)=EXP(ALOG(2(5!IJJ))(ALOG(B(5Ii.JJ))ALOG(32(4,Ii.J
IJ)) )/(ALOG(OMB(5) )-ALOG(QMB(4) ))*(ALOG(OM)-ALOG(QM815))))
G0 TO 300
30 CALL TLU(ONW(IIJJ).QMB,82(1,II.JJ)95)
300 CONTINUE
DO &00 KKki#3
CALL TLU(ALBT*AX(IKK),IETAgW(1.KK)93)
318
400 CONTINUE
CALL TLIJ(AFMX12*ALPHA9XA,3)
CALL TLU(I3OA9XCvGAA9C,4)
CALL TLU(HOAgXN19GAAtN1.4)
CALL TLIBOAgXN2.GAA9N294)
IF(BOA*GJ.5.) PHIS=J.
IF(BOA.GT.5,) GO T0 450
P4IS=(XC/C(1) )*( (UBAR*GG)**(XNI-N1 (1)) )*C(PHZE*(3,i41l,9/2.) )**(
IXN2-N2(I)))
450 PHIT=XFAB*(1..(-.1*SLIP))*(.-X82P1ZE)
1010 FORMAT(//,1A,3HF =qE12o5)
1011 FORMAT(1X94HO2 =gE12'.5)
1012 FORNAT(1X,3HC =9E12.5)
10123 FORMAT(1A,4HN1 =9E12.5)
1014 FORiMATC1Xv4HN2 .=9E12.5)
RETURN
END
319
APPENDIX VIII
All the calculations mentioned in Section IV-B are included in the computer
program BALFTN. In addition to calculating the frictional forces and torque
at a ball-race contact, it provioes the option to compute the derivatives of
forces and torque with respect to speeds and load. With the aid of Fig. VIII-
1, where a ball-race contact is shown, the instruction for input cards and
the description of output data are given below.
I/ 320 /
r/
Card-3 Read R, RG, RR, Ul, U2, VI, WS
FORMAT (8E10.3)
If NGLG = 2, this card should be omitted.
R, RG, RR = have been explained in the instruction for
Card 2.
Ul = u in Fig. VIII-1, velocity of ball in the rolling
direction x, in in/sec.
U2 = ur in Fig. VIII-I, velocity of race in x-direction,
in in/sec.
Vl = v in Fig. VIII-l, lateral sliding velocity of ball,
in in/sec.
WS = Ws in Fig. VIII-l, spinning velocity of ball about
z-axis, in rpm.
-/
/ 321
/R R
/+
Al/
C*/
Fig. /III
ABl-aeCnat
32/
MTI9/4
II
VIII-2. OUTPUT EXPLANATION
The first section contains the write-out of the input data. Symbols
and units used have been explained in the instructions for input.
The still next section represents the calculated forces and moments
acting on the ball. In each line the second number is a normalized
quantity. Contact load. P is isei as the normalizing force and the
product of the load and ball radius. PR is employed as the ,ormaliz-
ing moment.
The last section (if NDV = 1, this section is omitted) contains the
derivatives of forces and moment with respect to speeds and load.
Let Fx be the force in the x-direction acting on the ball, F the
force in the y-direction, and Mz the moment acting on the ba~l about
the z-axis:
6F
DFX/DU - 6 , *in lb-sec/in.
6F
DFX/DVS - x .in Ib-sec/r.
6F
DFX/DWS - W xs in lb/rpm.
6F
DFXiDW = P , in Ib/ib.
6F
DFY/DU = 6u , in lb-sec/in.
S6F
DFY/DV = -"V in lb-sec/in.
6F
DFY/DWS = Y in lb/rpm.
6Ws
6F
DFY/DW = 6P- in lb/lb.
323
DHZ/DU =-
Cu
in lb-sec.
6m
DAMLIDV = u ' i lb-sec.
E*1Z/DWS = 6Z , in in-lb/rpa.
6z
D&/iDW =-;- in in.
324
rC C
C
C
PROORAN BALFTN(INPbT.O)UTPUT.TAPE5ZINPUT.TAPE6ZOUTPUT)
FRICTION CALCULATION OF ANGUJLARC CONTACT BALL BEAR~ING
INPUT DESCRIPTION---
RzRAOIuS OF BALL.Ih.
PG=W.OOVE ukADJUS 04- iACE91%.
C RR=RADILUS OF RACEvIh.
L C
C
C
TJ=INLET TEMPER~ATURHE* OEG.F.
wRzINOLLING ANG. vEL. OF tIALL. mPe4.
VL=LATEkCAL ROLLING vEL. OP BAlL. RPM.
wS=SPINNING VEL. OF 6ALL. k1PM.
C
C W2=ROLLING VEL. OF iWAC~v kPM.
C UlvVl=ROLLINiG AND LA~TI'AL VEL. (IF SALL. IN/SEC.
C W=TOTAL LOA0. LBF.
C U2=.RjLLING VEL. OF oiACh. IN/Sr.L.
C VIS=INLET vISCOSIIY. LdF-StCIN**2
C ALP11A=PASSU(RE-V15COSITY LtJEF.s IN**2/LdF.
C BETATEPENATUtRE-VISCOSITY COLI-.. /OEG.F.
C FK=rTlEWHfAL CONSUCTIVIlY Of LUd1CICANT. 8fU/OEG.F-tlk-FT
DIMENSION ADA(24)qFMA(24)vAOdA(24)
DIMENSION U(3).V(3).gSS(3).wT1(i).SW()SA3)S1(3).P1Z(3).Y(60).
* IiATIO(6)FFXX(39396)jFFYY(39396).TM4(3.3.6)
COMMON t,1.62vG3.UCOE. THETA
DATA(AMGA(l)#I=I#24)/ .0.9.190.15820.76bO.U.1l0.0.59.40.4
190.570.,5.6o.tO..440.75.O.t$.8.3tO.90.O240. 9b49.90.489./
GO TO 30
*20 READ(S.2) W99Gv.ksU1.U2#Vi~wS
30 RFAD(592) ALP11AvbLtA.VI~oFKvWT I
IF(OV.i(.3 60 TO 31
RFAU(S.2) fOU.DOV'.S.w7
U(2)U,Cj
U (3)=U1-OU
V(3)=V1-DV
-. WSS(2)=dS.DwS
4SS (3 =dS-LowS
MTT (i')=w7wT~
4TT (J)=wl-DwT
2 FORHAT('8EIO.3)
3 FOIRMAT(// 3H 9= E12.4.bA 4H1 l~b= E12.4#5A 411 Nkz E12.4)
31 WRITE(6.3) Ndot3.km~
4 FORM4AT(/ 411 UI= Ele.445), 4,H U2= EIZ.495X 4H1 Vi= L12.495A 4H1 WS= El
12.4)
ORITE(6941 ux9U2.VI*W5
dRITt(6s,5) ALPHA9t3ETA.'VIS
5 FORMAT(/ 711 ALPH1A= le1.4#5A 6H1 6LTA= E12.4.hX 1111 VISCOSITY= E12.4
325
o FORMAT(/ 6ml LOAU= c.12.495A 311 I:. E12.4)
DY:2?./ 10 Y
1)0 -* 1=1-3
SBII)=FP4IO.71>3.1..1bluTTII)/CAA.0O[)/EP)**O.33333
HNm=.5ALP1A.f3VI
-*O.7EO*#.O3*WAt'0.43
GIT=VIb3OIALPrA/AL)*O.hI(VIS/VIS3)**O.I
G2T=VI53*e-ETA/ I8.*tK*O.2Ib) .
1=1
J: I
K: I
L~l
N: I
48 FFY=O.
FF x0.
Tm2i:O.
JFiSS.EO.O.) rO 10 4J
YS: ItJI) -Us?) /SS
GO TO 44'
'.3 Y(I)I)..OY/2.
GO TO 4.2
4.4 NN=1
326
IF4y(ft).L!.(1.-0T/2.)) G0 T0 41
y (PW) =yl;-Oy
RAT 10("04)=1..- Y M) -*2
42 yIN"* 1 =y om-JY
RAT' 1 P?#)=.- A#1 *
OWVM)LT(1) NMZP4M-
IF 1 IV.EO.2) GO TO 4.9
1#9 DO 50 M=2*W,4
FFX=FFX*FX
VFY:FFW*FY
TFX=-FX'*CP)*St(L)
TN1= Tml' Fx
TFY=-FCGfuSS*SA (Li ***AT10 (MI /3./UtpL)Y
TM2=TM2*TFY
IF (NDV.EO.2) G~o TO 5v
WRITE cb9l0J YN,,I'Lt.1L.3CEFVT~.r
10 FO.RMATC12EI1.3)
50 CONTINUJF
TM4:TMl#7M2
FFXX (I .J#NI=FFA
FFYY(19J*NU=FTFY
IF(Msv.EU.I) $0 TO 80
IF(C.$E3).R.(.Gifl~I-. ~.i1.1.ORcL.T.I1 6 10 :
Go TO 48
55 I1
IF(J.E.3).R..6.I.OkCL-Ub1i CGO Toj 56
J=J.1
$0 10 4.8
56 irl
!F((K.GE.31.OP.CL.G7.1); GO TO 57
N=N. 1
$0 TO 48
57 K1I
IF(L.CGE.3) $0 TO 58
327
80 FFFASFFXX(19191)/WT
FFFY:FF IV0*1,91)/wI
WRITE(6911) FFAA(1921).IFFrX
11 FORMAT(// 44$ FRICTION IN ThE 01I*CTION OF ROLLING (Lb.)s E120#954.
sr30.5)
12 FORMAT(/ 41$ FRI1CTION IN ItiE LATER~AL UIIfECTIOe 11.8.)z EIZ.49bx FIS
1.5)
WRITE(6912) FFYY(1.11)eFFFY
13 FORMAT(/ 26$ SP~INNING IOR~uE (IN-LB.)s E1ZS.45X FI@.b)
WRITE(6913) TMN(1.1.IhtFTN4
IF(NOV.EO.1) 60 T0 120
wRJTE(6*81. FxOU
dl FORMAT(// Z0$ OFA/Ou L(/IN/SEL)s E12.4)
WRITE(6982) FXOV
82 FORMAT(/ 20$ Oix/OV (L8/It./sEC)= C12.4)
WRITEc6t83) FXLIWw)
d3 FORMAT(/ jdh OFA/DwS (LB/5RPM)z E12.4)
WRIJE(6.84) FXDW
84 FORMAT(/ 16m OFX/Ow (LB/Lb)x E12.4)
WR!TE(6905) FYDU
8S FORMAT(/ 20H DFY/OU (LH/IN/SEC)= E12.4)
WRITE(6.86) FYOV
86 FORMAT(/ 20m DFY/Ov (Lb/lh/SEC)z E12.4)
WRITE(6987) FYOWS
87 FORMATI/ Idn DF'r/OwS (L9/hPq= E12.4)
WRITE(6988) FYOW
88 FORMAT(/ 16h DFY/LUw (Li6/Lb)= E12.4)
WRITE (6989) TODUj
89 FORMAT ( 23$ DMZ/U, (.I-Lb/1N/!SEC)=E12.4)
WRITE(6990) TMOV
90 FORMAT(/ 23'l DMZ/OV (:N-LII/IN/SEC)= E12.4)
WRITE(69911 THOwS
91 FORMAT(/ 21$ DOZOwS (lN-Lu/R'M)z E12.4)
WRITE(6992) TM~w
92 FORMAT(/ 19$ Ln4Z/Ow (IN-LB/LB)= E12.4)
120 STOP
END
SUBROUTINE FNCTN
DIMENSION AGI(131-.AG !(?).AU3().FTN(1J.1,7).UI(17).VZ(7)
COMMON G19629639FCOE*THiETA
DATA(FIN(1).I=1,13)/0.0002b.O.0007, Oo0O14. 0.0023. 0.0067t
1 0.00359 0.0219 0.031. 0.0469 0.052. 0.0615v._065v 0.0775/
OATA'(FTN(I)9I=i49Z6) /0.00079 0.002. 0.00399 0.00659.0*0150r
1 0.02 0.029. 0.042. U.04959 0.0549 0.0625. 0o069s'0.077tj/
DATA(FTN(1),I=27*J9)/ 0.0026. 0.0078, 0.014. 0.018. 0.028.
1 0.034- 0.0385. 0.04859 0.U549 0.058. 0.0649 0.0695, 0.0775/
DArA(FTN91=40952) /0.0077. 0.017.0.0239 0.027. 0.035.
328
DATA(FJNt(I)vl63v6S)/0.019.v 0.028b. 0.0340 0.037. 0.044. 0.046.
1. 0.05 * 0.056. 0.06. 0.06359 0.069. 0o0125. 0.0775/
DATAfFTft(j),13669?b1 0*O02159 0.032. 0;0379 0.0416. 0*041S9
I 0.051. 0.054. 0.06. 0.0635. 0.066. 0.01. 0.013. 0.0775/
DATA(FTtd(I)91=79991) / 0.0245.0.0349 0.04, 0.0435. 0.06,
1 0.054. 0.0565. 0.062. 0.0655. 0.0675. 0.071. 0.0749 0.0i75i
DATA(FTN(I).1z929104)/0.0002,0.0006. OoOO129 0o.U29 0.006. 0.012 9
1 0.019. 0.036. 0.047. 0.064. Oo06459 0.072. o*PnIZ
DATA(FThf11,I=106,117)/0.0006.0.002, 0.004. 0.0665. 0.016. 0.0249
1 0.028. 0.0439 0.052, 0.0S56 u.066S. 0.012. 0.07
DATA(VTId(1)91z118.130)/ 0.00279 0.0081s 0.014. 0.01 9 0.0285,
10.035a 0.0395. 0.Oblo P1.0579 0.061. 0.068. 0.0725. 0.077/
DATA(FTtII).1=131*143)/0.00669 0.0kb. 0.021. 0.0259--t.0349 0.04 9
1 0.0445. 0.054. 0.06. 0.06.. 0.069. 0.073. 0 077/
DATA(FTMII).1=144*156)/ 0.012. 0.0225. 0. 0;49 0.0345b. 0.044.
1 0.048. 0.052. 0.0599 0.0639 .06. 0.070. 0.073. 0.077/
oATA(FIN(I);I=1579169)/ 0.019. 0.032. 0.039. 0.04159 0.05,0.05..
1 0.056, 0.061. 0.0649 0.061. 0.07159,0.0759 0.077/
DATA(VTNII)9j1170.101)/ U.021. 0.0335s 0.041. 0.04559. 0.063.
1 0.056. 0.059. 0.0649 0.068. 0.01. 0.073. 0.016. 0.017/'
DATA(VTN(1) .I=183.19b1/.00016.OOO5b0.00119.0.019.0.001!8.0.011.
1 0.0181 0.033. 0.0429 0.04d. 0.06. 0.0689 0.Or/- .kb .03
DATA(FTNf1).I=196.208J/.00079.U021. 0.00439 0.0079 .01o00
1 0.027. 0.0399 0.046. 0.062. 0.062. 0.07. 0.015/
DATA(V1N(1)vJ=109.*221)/ 0.Ob1dq 0.00,3.^ 0.01. 0.016. 0.02b551
1 0.0349 0.0379 0.048. 0.054. 19.0509 0.066. 0.074. 0.015/
DATAiF7N(I).I=22223.k/ 0.00469 0.011. 0.016. 0.024. 0.0379
U 0.039. 0.0449 0.05259 0o067, 0.061. 0.068. 0.074o 0.075f
DATA(FTNfIJ.I=235v247)/ 0.0082. 0.017. 0.02S9 0.03. 0.04Z,
IL
329
DATA(FT~d(DI.I439-'5,/0.0I2v 0.0229 0*02959 0o0359 0.04S5.
1 0.0519 0.OiSe 0.061. 0.0639 0*0659 0.067. 0.066. 0.068/
OAYAC(FIN(I 191=564o)/.00005..Q0014..000270.00045..0014
1 0.60279 0.0046. 0.015. 0.023. 0.0299 0.0429 0.052. 0.Ob5/
DATA(FTNt!).1=69941).000129.00042.00085.0014..00449
1 0.00840 0.0130 0.02'. . 0.0329 0.037. 0.04559 i.653, 0.0t5/
0ATA(ru.(1).I=4832..ovd/.0004.O0l2 0.0024. 0.004. 0.011.
b~ 0.017. 0.Oe3o 0.03459 0.039. 0.04.259 0.048. 0.0535. 0.065/
- ATA(FIN(I)9!'.95#5d7)/0.00IS9 U.00429 0.001s 0.012. 0.023.
1 0.029. 0.0339 0.039b. 0.04.39 0.045. 0.0soo,*0.05JS. 0.066/
DATA;FlN(I)*I=5oSd.20)/Q.00349 0.0084.. 0.014. 0.019S9 0.o0s.
1 0.0369'*0.U. 0*04.69 0.048. 0.049. 0.051. 0.05'. 0.065/
OATA(FT9N(DI.I133)/0.00569 0.013. 0*0Z9 0.0e6t 0.Q3759
1 0.04.d. 0.045. 0.0499 0.056. 0.O5et 0*0655 0.051#. 0.0!5/
DATAIFTN(IJI:9=*3.i46,)/0.0071. 0.015.e 0.023. 0.03. 0.04.
1 0.045. 0.04d. 0.051- U.0529 0.0539 0,0549 0.055. 0.055/
DATA(FIN(1).I=,'.7.65'9)/.00003..00009..00018,.000J..0009.
1 0.0019 0.00J9 0.009. 0.016. 0.0e3. 0.0s1. 0.0459 0.051/
O)ATA(F7N(I).1='0.61e/.00019.O026..00058,.000V.9.00299
1 0.0057. 0.009?. 0.021. 0.0299 0.034. 0.041. 0.0479 0.051/
DATA(FrN(I).I= 73.6)dS)/.000e~5..0007490.0015. 0.00249 0.00729
1 0.013. 0.0169 .0e~ 0.03 , 0.039. 0.04S. 0*0'.d. 0.051/
DATAIFTN(1U,1=tb.9,49S/0.000t)49 0.0016. 0.00a39. 0.005Z9 0.014.
I 0.6k2. .U.026i. 0.03:3, 0.04. 0.042. 0.04659 0.049. 0.051/
OATAtFI.N(1).I=,9.o11l)/0.0013. 0.003b. 0.0069. 0.0119 0.023.
1 0.0/9. 0.0356 0.03Yiv 0.044. 0.046. 0.04.859 0.050. 0.051/
DATA(Fl1NC1).I~bl .624)/0.0015, (J.U410 0.008. 0.012. 0.02359
I 0.0295. 0.013. 0.04. 0.044. 0.0465. 0.0494 0.0!)05. 0.05L/
DATA(FIN(I)=.n~r26,c3I)/0.001139 0.005e, OeV0949 0.014. 0.02S59
1 0.031. 0.('3to 0.042. 0.04 . 0.0479 0.05. 0.0510 0.051/
60 U7C!.J)=IPTNI.,.1.j)
60 It' luO
70 I*J)=t.PALOG(Pi NtlI~l-J)*ALOG(FTN(2,I9J))-ALOG(ITN(1.19J)l))/-
1ALOG(AGI(2fl-ALOG(AC,1(1fli'(ALOG(bl)-ALOGCAGI(l))))
100 CONTINUE
00 200 A=197
IF(63-AGj(7)) 14U,160.1c'U
140 IF(G3.LT.AG3(I)) OQ Tu 170
CALL TLUfC6i.V/(K).AG3.uZ(1.*K),7)
60 TO e00
160 VZ(K)=uiC7.K) .
GO To 200
170 VZ(K)=01(.K)-(LL I/K)-UL(1.K))/(AG3(/1-AG3(1) )*(AG3(1)-G3)
IF (VZ(FK).LE.0.0) VZ(K)=1.E-tp
200 CONTINUE
330 *
X CALL TLIJ6(629FR39AG29VZ97j
6010O2101
206 FR3sEXP(ALOG(VL(6)).(ALOG(VZ(1))-ALOG(VL(6)))/(ALOG(AGZ(7))-ALOG(A
162(6)) )*(ALOG(62)-ALO4,(AG2(6):))
162(2)) )O(ALO6(G2)-ALO6(AG2(2)))
210 FCOE=ALOG(FR3)-0.1'9*(ALOG(T"EIA-ALO&(30.))
FCOEsEXP (FCOE)
GO TO 250
240 FCOE=O.
250 CONTINUE
RETUIkN
END
SUBROUTINE TLU(A~bvC909N)
DIMENSION C(1)#D(l)
IF (N-i) 1.2.3
IO T81O100
2 8=0(1)
Go To 100
3 ML1I
8IF(MU-ML-1) 15.1b.Y
9 M=(MU*ML)/2
10 IF(C(M)-A)13.l,14
11 IF(A-C('4))I3,le*14
12 8=n (m)
60 TO 100
13 MU=H
- GO TO
14ML=M
15 G=(-OTU I(M)-((A-C(ML))/(CiS4U)-C(ML)))
100 RETUriN
ENr)
SUBROU7INE TLUt,(A969C909N)
DIMENSION C(1)90(1)
IF (N-i) 1,2.3
L1 4=0.
Go To 100
2 8=D(1)
GO To 100
3 ML1l
MU=N
8 IF(MU-ML-1) 15915,9i
9 M=(MUML)/Z
10 IF(C(14)-A)13912914
11 IF(A-C(M))13*Ie914
12 O=D(M)
GO TO 100
14 ML=14
GO TO H
15 B=EXPLILOG(D(ML))*(AL~b(V(MU))-ALOGAUt(ML)))'((ALOG(A)-ALUG(C(ML)))
1 /(ALUG(C(MU))-ALOG(C(mLWf))
100 RFTUHN
END
F 331
"
11CT A Tr- 'TT
Secunty Classification
DOCUMENT CONTROL DATA -R & D
(Security classification of title, body of abstract and indexing annotation must be entered when the overall Aport is classifed)
I ORIGINATING AC TIVITY (Corporate author) 2ZA.REPORT SECURITY CLASSIFICATION
N/A
3 REPORT TITLE
Contract F33615-69-C-1305
b. PROJECT NO. 3048 None
Task 304806
c. sb. OTHER REPORT NO(S) (Any OM*rnumberat
,at may be asaigned
th e
Tport)
,, None
10. DISTRIBUTION STATEMENT
This document ir subject to special export controls and each transmittal to foreign
governments or foreign nationals may be made only with prior approval of the Air Force,
Aero Propulsion Laboratory, Wright-Patterson Air Force Base, Ohio 45433
I1. SUPPLEMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY
DD INOV01473 UNCLASSIFIED
Security Cl~ssification
Security Classification
14. LINK A LINK 9 LINK C
It
KEY WORDS
ROLE WT ROLE WT ROLE WY
Security Classification