Next Article in Journal
An Arbitrary Lagrangian–Eulerian Formulation for Modelling Cavitation in the Elastohydrodynamic Lubrication of Line Contacts
Next Article in Special Issue
Experimental and Numerical Analysis on the Seizure of a Carbon-Filled PTFE Central Groove Journal Bearing during Start-Up Period
Previous Article in Journal
Influence of Laser Pulse Number on the Ablation of Cemented Tungsten Carbides (WC-CoNi) with Different Grain Size
Previous Article in Special Issue
Comparison of Perturbed Reynolds Equation and CFD Models for the Prediction of Dynamic Coefficients of Sliding Bearings
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

EHD Effects in Lubricated Journal Bearing

by
Pavel Pavlovich Usov
Department Cathedra of high mathematics 2, National Research University of Electronic Technology, 124498 Zelenograd, Moscow, Russia
Lubricants 2018, 6(1), 12; https://doi.org/10.3390/lubricants6010012
Submission received: 30 October 2017 / Revised: 16 January 2018 / Accepted: 16 January 2018 / Published: 23 January 2018
(This article belongs to the Special Issue Friction and Lubrication of Sliding Bearings)

Abstract

:
This paper presents a numerical analysis of the influence of deformation of infinite parallel cylindrical solids in partial journal bearing on the oil film characteristics. The stationary elastohydrodynamic EHD problems for three design models of bearings are considered: (1) The bearing in which the basic contribution to the elastic displacement of the surface brings the thin elastic liner; (2) The elastic cylinder and the elastic bushing which is modeled by the elastic space with a cylindrical cut; (3) The elastic cylinder and the elastic bushing in the presence of the thin elastic liner with a small module of elasticity. It is shown that when the minimum film thickness is fixed and deformations of the elastic solids increase, then the load capacity increases, reaches a maximum, and then decreases. The deformation of solids can raise load capacity many times over. When the deformation of solids increases from zero, the pressure distribution changes from the distribution of pressure in the case of rigid bodies to the distribution of pressure which takes place with the dry contact of elastic bodies.

Graphical Abstract

1. Introduction

Nowadays, most hydrodynamic journal bearings are required to operate with decreasing film thicknesses. When a bearing is highly loaded, the elastic displacements can exceed the minimum film thickness and can essentially influence the bearing characteristics.
The deformation effects on the characteristics of a plain journal bearing have been studied in a number of works. Higginson [1] and Angra et al. [2] presented an analysis of the effects of the elastic deformation of the bearing liner on the performance of a journal bearing. The calculations were made at high values of the minimum film thickness. The results of the calculations show that at a fixed minimum film thickness the load grows as the liner becomes more flexible, but this growth is not great.
An analysis of the EHD behavior of a bearing operating under conditions which tend to be severe was presented by Bendaoudl et al. [3]. The influence of mechanical deformations were presented at one load. It was shown that the minimum film thickness was 18% less when deformations were not taken into account in the calculations.
The significance of the deformation effects on a plane journal bearing subjected to severe operating conditions was studied by Bouyer et al. [4]. It was shown that mechanical deformations significantly decrease the maximum pressure, significantly modify the profile of the lubricant gap, but do not have a great influence on the minimum film thickness magnitude. These conclusions were made on a basis of the analysis of the results of calculations for two loads.
Chetti [5] and Osman [6] studied the load dependence on eccentricity for rigid bodies and deformable bodies. The obtained results for elastic bearings was lower than those for rigid bearings for the whole range of eccentricity.
In all previous work the solutions were proposed only for several loads and for large minimum film thicknesses. They do not give a representation of how mechanical deformations influence the characteristics of a lubricant layer in the transition from small deformations to big deformations.
In the present work, a numerical analysis of the influence of deformations of infinite parallel cylindrical solids in partial journal bearings on the oil film characteristics is presented. It is shown that when the minimum film thickness is fixed and the deformation of the elastic solid increases, then the load capacity increases, reaches a maximum, and then decreases. The deformation of solids can raise the load capacity many times over.

2. Design Model One

Let us consider the conformal contact of a hard shaft and thin elastic liner that is inserted into a rigid bush (Figure 1). The inner radius of the liner is a little different from the radius of the shaft R , which is Δ / R < < 1 , where Δ = R 1 R , R 1 is the inner radius of the liner.
Let the line of center (Figure 1) of the bearing and shaft form with the line of action of the load angle θ , which is defined when solving the problem. The angular coordinate φ is measured from the line of action of the load in the direction of movement of the surface of the shaft.
Consider the case of a partial journal bearing. In this case, the inlet coordinate of a lubricant layer is determined by the position of the input oil-groove, from which the oil arrives under slight pressure into the gap. Then the inlet coordinate of a lubricant layer is equal to ϑ . The coordinate of the output point of a lubricant layer φ c is unknown and is determined when solving the problem. This coordinate should satisfy the condition φ c ϑ c . Here ϑ c is the coordinate of the inlet boundary of the output oil-groove in which the lubricant arrives after passing the zone of high pressure.
If condition (1) is satisfied, then the deformation of the shaft and bush can be neglected when determining the film thickness, and only the deformation of the liner must be taken into account.
ε R > > E E 1 , ε R > > E E 0
Here ε is the thickness of the liner, E ,    E 0 and E 1 are the elastic modules of liner, shaft and bush respectively.
The elastic displacement of the shaft surface w 0 is proportional of the shaft radius R and inversely proportional to the modulus of elasticity of the shaft E 0 , that is w 0 = C 1 R / E 0 . The elastic displacement of the bush surface is w 1 = C 2 R / E 1 and the elastic displacement of the liner surface is w = C 3 ε / E . Here C 1 , C 2 and C 3 are some constants. The conditions (1) follow from these expressions and from the conditions w 0 < < w and w 1 < < w .
When the additional condition ε R < < 1 is correct, then the displacement of the liner surface w can be defined as [7,8]
w = ε ( 1 + ν ) ( 1 2 ν ) ( 1 ν ) p ˜ E
Here p ˜ is pressure which acts from the lubricant layer to the surface of the liner, and ν is the Poisson’s ratio of liner. According to (2), the film thickness h ˜ can be defined as:
h ˜ ( φ ) = Δ e cos ( φ θ ) + ε ( 1 + ν ) ( 1 2 ν ) 1 ν p ˜ E
Here e is the eccentricity. The pressure we define from the Reynolds equation
1 R 2 d d φ ( ρ ˜ h ˜ 3 12 μ ˜ d p ˜ d φ ) u 2 R d ( ρ h ˜ ) d φ = 0
Here ρ ˜ and μ ˜ are the oil density and viscosity respectively, and u is the velocity. The oil density and viscosity depend on the pressure as follows
ρ ˜ = ρ 0 ( 1 + C 1 p ˜ 1 + C 2 p ˜ ) ,   μ ˜ = μ 0 exp ( α p p ˜ )
Here ρ 0 and μ 0 are the oil density and viscosity at atmospheric pressure respectively, С 1 and С 2 are constant and α p is pressure coefficient of viscosity.
The boundary conditions for the pressure are of the form
p ˜ ( ϑ ) = p ˜ ( φ c ) = 0 d p ˜ d φ | φ = φ c = 0 ,   if   φ c < ϑ c
Condition (5) means that if φ c > ϑ c is obtained when using the condition d p ˜ d φ | φ = φ c = 0 then it is necessary to accept φ c = ϑ c and to solve the problem without using the condition d p ˜ d φ | φ = φ c = 0 In this paper we will suppose that the angles ϑ and ϑ c are such that the condition d p ˜ d φ | φ = φ c = 0 is satisfied. Integrating Equation (3) in the domain ( φ , φ c ) and taking into account boundary condition (5), the following equation can be obtained
1 R 2 ( ρ ˜ h ˜ 3 12 μ ˜ d p ˜ d φ ) + u 2 R ρ ˜ h ˜ = u 2 R ρ ˜ 0 h ˜ c
The conditions of the balance of forces look like
R ϑ φ c ( p ˜ cos φ τ ˜ sin φ ) d φ = W ,   R ϑ φ c ( p ˜ sin φ + τ ˜ cos φ ) d φ = 0
Here W is the force, and τ ˜ = μ ˜ u h ˜ + h ˜ 2 1 R d p ˜ d φ is tangent stress acting from a lubricant layer on a shaft.
The following dimensionless variables are defined
h = h ˜ Δ ,   p = p ˜ p L ,   p L = Δ ( 1 ν ) ε ( 1 + ν ) ( 1 2 ν ) E ,   ρ = ρ ˜ ρ 0 ,   μ = μ ˜ μ 0
Here p L is characteristic elastic contact pressure. When p ˜ = p L the displacement of the liner surface w equals the radial clearance Δ .
In dimensionless variables, the equations and conditions can be represented as follows
h ( φ ) = 1 χ cos ( φ θ ) + p
ρ h 3 exp ( G p ) d p d φ = V ( ρ h h c ) ,   h c = 1 χ cos ( φ c θ )
ρ = 1 + K 1 p 1 + K 2 p
ϑ φ c ( p sin φ + τ cos φ ) d φ = 0 ,   τ = ψ ( V 6 exp ( G p ) h + h 2 d p d φ )
θ φ c ( p cos φ τ sin φ ) d φ = P
p ( ϑ ) = p ( φ c ) = 0
Here ψ = Δ R , χ = e Δ is relative eccentricity
V = 6 u R μ 0 Δ 2 p L = p * p L = 6 μ 0 u R ε ( 1 + ν ) ( 1 2 ν ) ( 1 ν ) Δ 3 E
p * = 6 u R μ 0 Δ 2 is the characteristic hydrodynamic pressure, K 1 = p L C 1 , K 2 = p L C 2 , G = α p p L , and P is dimensionless force which is equal to
P = W R p L
Dimensionless parameter V equals the ratio of the characteristic hydrodynamic pressure to the characteristic elastic contact pressure. When this parameter is small the pressure p ( φ ) is also small compared to characteristic elastic contact pressure, and deformations of the liner have little effect on the solution of the problem. It means that at small values for parameter V the solution of the problem differs slightly from the solution for rigid bodies. When parameter V increases, the influence of the deformations of the liner on the pressure and film thickness also increases.

2.1. Asymptotic Properties of the Solution

Let us consider the solution of the problem at large values of parameter V and when the influence of pressure on lubricant density and viscosity is not taken into account, that is at ρ = 1 , μ = 1 . It follows from (7) that in this case
h h c = h 3 V d p d φ
Since the pressure and pressure derivative are limited then
h h c   at   V
It follows from (6) and (14) that p p at V where
p = h c 1 + χ cos ( φ θ )
Let us consider the case when ϑ c = ϑ and φ c = ϑ . The angle θ is proportional to ψ when function p ( φ ) looks like (15). Parameter ψ is usually small ( ψ < 10 3 ) and we can accept θ = 0 . It is not difficult to obtain in this case
lim V χ = 1 h c cos ϑ ,   and   p = ( 1 h c ) ( cos φ cos ϑ 1 )
lim V P = ( 1 h c ) ( ϑ cos ϑ sin ϑ ) ,   lim V p max = ( 1 h c ) ( 1 cos ϑ 1 )

2.2. Numerical Method

The input parameters of the problem are K 1 ,    K 2 ,    G ,    V and χ . After setting these parameters, the functions h ( φ ) ,    p ( φ ) ,    ρ ( φ ) , μ ( φ ) ,    τ ( φ ) , and parameters φ c and θ can be found by solving the system of Equations (6)–(8) under conditions (9) and (11). The dimensionless load P can be found using the expression (10) after that.
The angles θ ,    φ c are defined by iteration. The initial approach of these parameters is assigned at first. Thus, the initial values of θ ,    φ c are specified before the calculations.
The mesh of N nodes { φ k } , k = 1 , 2 , , N , φ 1 = ϑ , φ N = φ c is entered. It follows from Equations (6)–(8), that
h N = h c = 1 χ cos ( φ c θ )
h i = h ( φ i ) = L i + y i ,    i = 1 , 2 , , N 1
L i = 1 χ cos ( φ i θ ) + p i + 1
y i = p i p i + 1 ,   p i = y i + p i + 1
ρ i + 1 ( y i + L i + h i + 1 ) 3 [ 1 exp ( G y i ) ] = 4 G V ( Δ φ ) exp ( G p i + 1 ) [ ρ i + 1 ( y i + L i + h i + 1 ) 2 h c ] i = 1 , 2 , , N 1
ρ i = 1 + K 1 p i 1 + K 2 p i
Equations (9) and (10) can be presented as (24) and (25) after replacing the integrals by finite sum.
k = 1 N a k ( p k sin φ k + τ k cos φ k ) = 0
Δ φ k = 1 N a k ( p k cos φ k τ k sin φ k ) = P
τ k = ψ ( V 6 exp ( G p k ) h k + h k 2 p k + 1 p k φ 2 φ 1 ) ,    k = 1 , 2 , , N 1
τ N = ψ V 6 h c ;    a 1 = a N = 1 2 ;    a k = 1 ,    k = 2 , , N 1
As θ and φ c are known, we can calculate h N = h c on Formula (18) and L N 1 on Formula (20). The condition p N = 0 must be taken into account. As h N , h c , L N 1 , p N are known, the value of y N 1 can be found using the Equation (22) and expressions (23). Really, after substitution of expression (23) into the Equation (22) and the assignment of i = N 1 , we receive a non-linear equation with respect to y N 1 . After solving this equation, the values of h N 1 and p N 1 can be defined under the Formulas (19) and (21).
Further, the value of L N 2 can be calculated and so on. As a result, the values of h i and p i can be calculated consistently at i = N 1 , i = N 2 , , i = 1 . Thus, the values of h i    ( i = 1 , 2 , , N ) and p i    ( i = 1 , 2 , , N 1 ) can be calculated if the values of θ ,    φ c are assigned. After that τ k ( k = 1 , 2 , , N ) can be computed on Formulas (26), (27) and parameter T can be determined.
T = Δ φ k = 1 N ( φ k φ k 1 ) ( p k sin φ k + τ k cos φ k )
The validity of the conditions (28) indicates that the solution of the problem is found.
| p 1 | < ε p ,   | T | < ε T
Here ε p and ε T are the parameters which determine the precision of the solution. They are set before the calculations.
The flowchart for determining θ and φ c is shown in Figure 2. Equation (22) is solved with respect to y using Newton iterations.
The number of nodes N determines the accuracy of the solution. The influence of number N on the accuracy of the solution is illustrated by the results of the calculation presented in Table 1. The calculation was made at the next values of parameter: ϑ = 1.2 , V = 0.01 , χ = 1.3 , G = 0 , K 1 = K 2 = 0 .
The results in Table 1 indicate that the results of the calculation have high precision when N 20 .

2.3. Discussion

At first, we will consider the results when the influence of pressure on lubricant density and viscosity is not taken into account, that is, the decision of a problem at G = 0 , K 1 = 0 , K 2 = 0 . Figure 3, Figure 4, Figure 5 and Figure 6 present the results at ϑ = 1.4 radians.
The ratio (29) is proportional to the dimension load W and does not depend on the elastic properties of the liner. The parameter V is inversely proportional to the modulus of elasticity of the liner. Therefore, the dependence of ratio P / V on parameter V characterizes the influence of the elasticity of the liner on the load.
P V = W Δ 2 6 μ 0 u R 2 = P p * R
Chetti [5] and Osman [6] presented the dependence of load from eccentricity for rigid bodies and deformable bodies. They came to the conclusion that at given eccentricities the load for rigid bearing is greater than that for elastic bearing.
Figure 3 shows the dependences of ratio P / V and minimum film thickness h min on V at χ = 0.9 . It follows from these results that the dimension load decreases when the elastic deformation of the liner increases as in [5,6]. However, the presented results also demonstrate that at a given eccentricity the minimum film thickness essentially increases when the elastic deformation increases. This means that the influence of deformations on the film thickness may be shown more demonstrably by the dependences of the ratio P / V on V at h min = c o n s t .
Figure 4 shows the dependences of P V on V at different values of h min . When V 0 then the solution of the problem approaches the solution for the rigid bodies. Therefore
lim V 0 P V = W h d ( h min ) p L R
Here W h d is the load in a pure hydrodynamic regime which depends on the minimum film thickness, that is W h d = W h d ( h min ) .
The value of coordinate φ c does not exceed the value of ϑ = 1.4 in all the presented results of the calculation. The pressure distribution approaches the pressure distribution for the dry contact of a rigid cylinder and an elastic liner at a large value of V . It follows from (17) that at ϑ = 1.4 and a large value of V there will be
P V 7.25 V ( 1 h min )   or   lg ( P V ) = lg ( 7.25 ( 1 h min ) ) lg V
Relationship (30) is presented in Figure 4 as a linear dependence.
It follows from the results presented in Figure 4 that the ratio P V as a function of V has a maximum. When the deformations of the liner are small, they lead to an increase in the load capacity of the lubricant layer. However, when the deformations of the liner are large, they can decrease this load capacity.
The elasticity of the liner can raise the load capacity of the lubricant layer many times over. For example, at h min = 0.02 this increase can be 18 times. At h min = 0.2 this increase may not be more than two times.
Figure 5 shows the dependences of h ( φ ) and p ( φ ) / p max at various values of parameter V and at h min = 0.02 . It follows from these results that the point of minimum function h ( φ ) moves to the left when V increases. The pressure distribution gradually passes from specific hydrodynamic distribution to pressure distribution in dry contact.
Higher values of parameter V correspond to a more flexible liner. The shaft penetrates into the liner more deeply if the liner is more flexible. This leads to an increase of the size of the contact zone when parameter V increases. The point of minimum film thickness was located near the output point of a lubricating layer. As a result, the point of minimum film thickness shifts towards right as parameter V increases.
Figure 6 shows the dependences of the maximum of pressure p max on parameter V at different values of h min . It follows from these results that the dependence of lg p max on lg V is close to linear dependence up to the limiting value
p max = ( 1 cos 1.4 1 ) ( 1 h min ) = 4.88 ( 1 h min )
The results presented above are received without taking into account the influence of oil density and viscosity variation. This influence is illustrated by the curves in Figure 7.
Figure 7 shows the dependences of ratio P / V on V at h min = 0.02 . Curve 1 is obtained when the influence of oil density and viscosity variation are not taken into account. Curve 4 presents results when only the influence of the viscosity variation is taken into account. Curve 2 presents results when only the influence of density variation is taken into account. Curve 3 presents results when the influence of density and viscosity variations is taken into account.
The presented results indicate that the density variation with the variation of pressure leads to a small variation in the load capacity of the lubricant layer. Basically, this influence becomes appreciable only in a small vicinity of the value of V at which the function P V ( V ) has a maximum.
The influence of the viscosity variation on load capacity also is small and it also is maximal in a small vicinity of the value of V at which the function P V ( V ) has a maximum. However, the influence of viscosity variation on load capacity is more important than the influence of density variation, because curves 3 and 4 are quite close to each other compared to curves 1 and 2.

3. Design Model Two

Now we will determine the deformation of the elastic shaft and elastic bushing using the solution of the plane problem of theory elasticity for the elastic circle and for the elastic plane with circular cut.
A general solution to the problem of elasticity for circular rings is given in a book [9]. It follows from this solution that the elastic displacement of the surface of elastic circles can be determined under the formula
w 0 ( φ ) = 2 ( 1 ν 0 2 ) R π E 0 ϑ φ c p ˜ ( ξ ) [ 1 2 + cos ( φ ξ ) ln | sin φ ξ 2 | ] d ξ + ( 1 2 ν 0 ) ( 1 + ν 0 ) R 2 π E 0 ϑ φ c p ˜ ( ξ ) [ ( φ ξ ) 2 π θ ( φ ξ ) ] d ξ + C ˜ 1 cos φ + C ˜ 2 sin φ
Here C ˜ 1 and C ˜ 2 are constants that must be specified when solving a specific problem,
θ ( φ ξ ) = { 1 , i f    ϑ ξ φ 0 ,    i f    φ < ξ φ c
The elastic displacement of the surface of the circular cut in elastic plane can be determined under the formula
w 1 = 2 ( 1 ν 1 2 ) R 1 π E 1 ϑ φ c p ˜ ( ξ ) cos ( φ ξ ) ln | sin φ ξ 2 | d ξ + ( 1 2 ν 1 ) ( 1 + ν 1 ) R 1 2 π E 1 ϑ φ c p ˜ ( ξ ) [ ( φ ξ ) 2 π θ ( φ ξ ) ] d ξ + C 3 cos φ + C 4 sin φ
Let us consider the case of E 0 = E 1 = E and ν 0 = ν 1 = ν . Taking into account the expression (31) and (32), the film thickness can be determined under the formula
h ˜ ( φ ) = Δ δ cos ( φ θ ) 4 R ( 1 ν 2 ) π E ϑ φ c p ˜ ( ξ ) [ cos ( φ ξ ) ln | 2 sin φ ξ 2 | 1 4 ]
The following dimensionless variables are defined
h = h ˜ Δ ,   p = p ˜ p S ,   p S = 1 8 Δ R E ( 1 ν 2 ) ,   ρ = ρ ˜ ρ 0 ,   μ = μ ˜ μ 0 ,   χ = δ Δ
In dimensionless form the basic equations look like
h ( φ ) = 1 χ cos ( φ θ ) 1 2 π ϑ φ c p ( ξ ) [ cos ( φ ξ ) ln | sin φ ξ 2 | 1 4 ] d ξ
ρ h 3 exp ( G p ) d p d φ = V ( ρ h h 0 )
h 0 = h ( φ c ) ,   ρ = 1 + K 1 p 1 + K 2 p
τ = ψ ( V 6 μ h + h 2 d p d φ )
ϑ φ c ( p sin φ + τ cos φ ) d φ = 0
ϑ φ c ( p cos φ τ sin φ ) d φ = P
Here
V = 6 u R μ 0 Δ 2 p S = p * p S = 48 μ 0 u R 2 ( 1 ν 2 ) Δ 3 E
K 1 = p S C 1 ,   K 2 = p S C 2 ,   G = α p p S ,   P = W R p S

3.1. Numerical Method

When calculating the value of h at the mesh point φ = φ i , the function p ( φ ) is approximated by a piecewise linear function (43)
p ( x ) = k = 2 N 1 p k S k ( φ )
where S k ( φ ) = { φ φ k 1 φ k φ k 1 ,     φ k 1 φ φ k φ k + 1 φ φ k + 1 φ k ,    φ k φ φ k + 1 0 ,    φ [ φ k 1 ,    φ k + 1 ] .
Substituting (43) into (35) and using the integrals (44), (45) the expression (46) can be obtained
cos t ln | sin t 2 | d t = I 1 ( t ) + c ,    I 1 ( t ) = sin t ln | sin t 2 | 1 2 ( t + sin t )
t cos t ln | sin t 2 | d t = I 2 ( t ) + c 1 ,    I 2 ( t ) = ( t sin t + cos t 1 ) ln | sin t 2 | 1 2 ( t 2 2 + t sin t + 2 cos t )
h i = 1 χ cos ( φ i θ ) 1 2 π k = 2 N 1 ( I k , i 1 + I k , i 2 ) p k ,   i = 1 , 2 , , N
I k , i 1 = φ i φ k 1 φ k φ k 1 [ I 1 ( φ k φ i ) I 1 ( φ k 1 φ i ) 1 4 ( φ k φ k 1 ) ] + 1 φ k + 1 φ k 1 [ I 2 ( φ k φ i ) I 2 ( φ k 1 φ i ) 1 8 ( φ k φ i ) 2 + 1 8 ( φ k 1 φ i ) 2 ]
I k , i 2 = φ k + 1 φ i φ k + 1 φ k [ I 1 ( φ k + 1 φ i ) I 1 ( φ k φ i ) 1 4 ( φ k + 1 φ k ) ] + 1 φ k + 1 φ k [ I 2 ( φ k + 1 φ i ) I 2 ( φ k φ i ) 1 8 ( φ k + 1 φ i ) 2 + 1 8 ( φ k φ i ) 2 ]
The Reynolds’s Equation (36) can be written in the finite difference form as
ρ i + 1 ( h i + h i + 1 ) 3 [ exp ( G p i + 1 ) exp ( G p i ) ] = 4 G V ( Δ φ ) { [ ρ i + 1 ( h i + h i + 1 ) 2 h c ] + 2 h c 3 d p d φ | φ = c } i = 1 , 2 , , N 1
The conditions of the balance of forces looks like in the case of mathematical model 1.
The variable h i , h i + 1 and ρ i + 1 can be excluded from Equation (47) by using expressions (23) and (46). As a result, these equations will represent the system of N 1 equations with respect to N 2 variables p 2 , p 3 , , p N 1 . This system of equations together with the equation of balance of forces defines the variables p 2 , p 3 , , p N 1 , θ and the variable φ c . The specified system of equations is solved using the Newton iterations method.

3.2. Discussion

Figure 8 presents results that are analogous to the results presented in previous section. The calculation where made at G = 0 , K 1 = 0 , K 2 = 0 and ϑ = 1.2 . The solid curves present the results of the calculation of the discussed design model two. The dashed curves present the results of calculation of design model three, which will be discussed later.
The presented results indicate that the solution of this problem has the same properties as the solution of the problem for bearings with a thin liner. In particular, the function Θ ( V ) = P ( V ) | h min = c o n s t V ( V ) has a maximum which shifts to the right with the growth of h min . The growth of the load capacity of the lubricant layer caused by the deformation of bodies can be more than in 3.5 times when h min = 0.05 . This growth can reach values of 1.5 when h min = 0.2 .
The calculations indicate that as in the case of model one, variation of the viscosity and density of the lubricant slightly influences film thickness and pressure.

4. Design Model Three

It follows from the results presented above that deformations can essentially raise the load capacity of a lubricant layer. The maximum load capacity is provided at a certain value of parameter V which looks like (12) in a bearing with a thin liner when the shaft and bush deformations can be neglected, and looks like (41) in the case of a homogeneous elastic shaft and a homogeneous elastic bush.
If the bearing is projected on a certain regime it is desirable to pick up the parameters so that the deformations as much as possible raise the loading capacity of the bearing. However, if materials with the specified module of elasticity are used and the specified lubricant is used, then there are a few possibilities to obtain a maximal load capacity by changing the value of parameter V . However, there is a possibility of increasing the load at the expense of the deformation of bodies by the introduction of the thin elastic liner with a small module of elasticity.
The given circumstance is illustrated by the results of the calculations for the scheme presented in Figure 9. The results of the calculations for this bearing scheme are shown in Figure 8.
The scheme of the bearing consists of an elastic shaft, elastic bush end elastic thin liner, the modules elasticity of which is small. In this case the film thickness can be defined as
h ˜ ( φ ) = Δ δ cos ( φ θ ) 4 R ( 1 ν 2 ) π E ϑ φ c p ˜ ( ξ ) [ cos ( φ ξ ) ln | 2 sin φ ξ 2 | 1 4 ] + ε ( 1 + ν l ) ( 1 2 ν l ) 1 ν l p ˜ E l
Here ν l ,    E l are the Poisson’s coefficient and module elasticity of liner, respectively. In dimensionless variables defined as (34), the Equation (48) looks like
h ( φ ) = 1 χ cos ( φ θ ) 1 2 π ϑ φ c p ( ξ ) [ cos ( φ ξ ) ln | sin φ ξ 2 | 1 4 ] d ξ + β p
where β = 1 8 π ( 1 + ν l ) ( 1 2 ν l ) ( 1 ν 2 ) ( 1 ν l ) ε R E E l .
The value of β = 0 corresponds to the case ε = 0 , that is, the case when the elastic liner is absent. The increase of parameter β corresponds to the increase in the liner thickness.
The dashed curves in Figure 8 show the dependences of the function Θ ( V ) = P ( V ) | h min = c o n s t V ( V ) at different values of h min and at β = 0.2 . The results presented in Figure 8 indicate that the thin liner leads to an increase in the maximum bearing load capacity and leads to a decrease in the value of parameter V at which the function Θ ( V ) has a maximum.

5. Example of Bearing Calculation

Determine the minimum film thickness of the bearing mill roll, the construction of which corresponds to design model two. The parameters of the bearing are: ϑ = 1.2   radian , R = 0.3   m , W = 1.02 10 7   N / m , ψ = Δ / R = 8 10 4 , u = 10   m / s , μ 0 = 0.01   Pa s , α p = 1.7 10 8   P a 1 , ν = 0.3 , E = 2 10 11   Pa .
According to the Formulas (34), (41) and (42) we find p s = 2.2 10 7   Pa , V = 0.142 , P = 1.545 and P / V = 10.88 . These values of V and P / V correspond to points in Figure 8, which is between the curves h min = 0.07 and h min = 0.1 . The calculations indicate that at G = 0 ,    K 1 = 0 these values of V and P / V correspond to h min = 0.079 and h ˜ min = R ψ h min = 19   μ m .
It follows from the results of the calculation that when V 0 and P / V = c o n s t = 10.88 then h min 0.0367 and h ˜ min 8.8   μ m . This means that the deformation increased the load capacity more than two times.
Consider the influence of the viscosity variation on the minimum film thickness. Using (42) we have G = 0.374 . It follows from the solution of system of equation that h min = 0.09 and h ˜ min = 21.6   μ m at this value of G and at K 1 = 0 . Thus the increase in lubricant viscosity leads to an increase in minimum film thickness of 13.7 percent.
Consider the influence of density variation on minimum film thickness. Using (42), we have K 1 = 0.0132 ,    K 2 = 0.0374 . The calculations give h min = 0.0906 and h ˜ min = 21.74   μ m at these values of parameters K 1 , K 2 and at G = 0.374 . Thus the increase in lubricant density leads to an increase in minimum film thickness of 0.6 percent.

6. Conclusions

The stationary EHD problems for three design models of bearings were considered:
(1)
The bearing in which the basic contribution to the deformation brings a thin elastic liner;
(2)
The elastic cylinder and the elastic bushing, which is modeled by an elastic space with a cylindrical cut;
(3)
The elastic cylinder and the elastic bushing in the presence of a thin elastic liner with a small module of elasticity.
It was shown that when the minimum film thickness is fixed and the deformations of the elastic solids increase, then the load capacity increases, reaches a maximum, and then decreases. Deformations of solids can raise load capacity many times over.
The maximum load capacity takes place at a certain value of dimensionless parameters which is proportional to the speed, the lubricant viscosity and inversely proportional to the elasticity module. It was shown that the presence of the thin elastic liner in the third settlement scheme of the bearing led to a reduction of the value of the dimensionless parameters at which the load capacity has a maximum. This means that it is possible to essentially raise the load capacity of bearings at operating conditions by the introduction of a liner with a small module of elasticity.
When the deformations of solids increase from zero, the pressure distribution changes from the distribution of pressure in the case of rigid bodies to the distribution of pressure that takes place in the dry contact of elastic bodies.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Higginson, G.R. The Theoretical Effects of Elastic Deformation of the Bearing Liner on Journal Bearings Performance. In Proceedings of the Symposium on Elastohydrodynamic Lubrication; Institution of Mechanical Engineers: Westminster, UK, 1965; Volume 180, pp. 31–38. [Google Scholar]
  2. Angra, S.; Mehta, N.P.; Rattan, S.S. Effects of elastic deformation of bearing liner on performance of finite offset-halves pressure dam bearings. Tribol. Lett. 1996, 2, 273–285. [Google Scholar] [CrossRef]
  3. Bendaoud, N.; Mehala, K.; Youcefi, A.; Fillon, M. An experimental and numerical investigation in elastohydrodynamic behaviour of a plain cylindrical journal bearing heavily loaded. Proc. Inst. Mech. Eng. J. 2012, 226, 809–818. [Google Scholar] [CrossRef]
  4. Bouyer, J.; Fillon, M. On the Significance of Thermal and Deformation Effects on a Plain Journal Bearing Subjected to Severe Operating Condition. J. Tribol. 2004, 126, 819–822. [Google Scholar] [CrossRef]
  5. Chetti, B. Analysis of a circular journal bearing lubricated with micropolar fluids including EHD effect. Ind. Lubr. Tribol. 2014, 66, 168–173. [Google Scholar] [CrossRef]
  6. Osman, T.A. Effect of lubricant non-Newtonian behavior and elastic deformation on the dynamic performance of finite journal plastic bearing. Tribol. Lett. 2004, 17, 31–40. [Google Scholar] [CrossRef]
  7. Rades, M. Dynamic Analysis of an Inertial Foundation Model. Int. J. Solids Struct. 1972, 8, 1353–1372. [Google Scholar] [CrossRef]
  8. Lahmar, M.; Haddad, A.; Nicolas, D. Elastohydrodynamic analysis of one-layered journal bearings. Proc. Inst. Mech. Eng. 1998, 212, 193–205. [Google Scholar] [CrossRef]
  9. Galachov, M.A.; Usov, P.P. Differential and Integral Equations of Mathematical Theory of Lubrication; Nauka: Moscow, Russia, 1990; p. 290. [Google Scholar]
Figure 1. Scheme of journal bearing. Rigid cylinder, rigid bush and elastic liner.
Figure 1. Scheme of journal bearing. Rigid cylinder, rigid bush and elastic liner.
Lubricants 06 00012 g001
Figure 2. Flowchart for determining θ and φ c .
Figure 2. Flowchart for determining θ and φ c .
Lubricants 06 00012 g002
Figure 3. Dependences of ratio P / V and minimum film thickness h min on V at χ = 0.9 .
Figure 3. Dependences of ratio P / V and minimum film thickness h min on V at χ = 0.9 .
Lubricants 06 00012 g003
Figure 4. Ratio P V as a function of V at G = 0 , K 1 = 0 , K 2 = 0 different values of h min .
Figure 4. Ratio P V as a function of V at G = 0 , K 1 = 0 , K 2 = 0 different values of h min .
Lubricants 06 00012 g004
Figure 5. Dimensionless film thickness and ratio p p max as a function of angle coordinate φ at G = 0 , K 1 = 0 , K 2 = 0 ,    h min = 0.02 and at different values of parameter V .
Figure 5. Dimensionless film thickness and ratio p p max as a function of angle coordinate φ at G = 0 , K 1 = 0 , K 2 = 0 ,    h min = 0.02 and at different values of parameter V .
Lubricants 06 00012 g005
Figure 6. Dimensionless maximum pressure as a function of V at G = 0 , K 1 = 0 , K 2 = 0 and at different values of h min .
Figure 6. Dimensionless maximum pressure as a function of V at G = 0 , K 1 = 0 , K 2 = 0 and at different values of h min .
Lubricants 06 00012 g006
Figure 7. Dependences of ratio P / V on V at h min = 0.02 .
Figure 7. Dependences of ratio P / V on V at h min = 0.02 .
Lubricants 06 00012 g007
Figure 8. Ratio P V as a function of V at G = 0 , K 1 = 0 , K 2 = 0 different values of h min .
Figure 8. Ratio P V as a function of V at G = 0 , K 1 = 0 , K 2 = 0 different values of h min .
Lubricants 06 00012 g008
Figure 9. Scheme of journal bearing. Elastic cylinder, elastic bush and elastic liner.
Figure 9. Scheme of journal bearing. Elastic cylinder, elastic bush and elastic liner.
Lubricants 06 00012 g009
Table 1. The values of parameters P and h min at various values of N .
Table 1. The values of parameters P and h min at various values of N .
N 2050100200500
P 0.330560.335850.336460.336600.33660
h min 0.042950.043120.043140.043150.04315

Share and Cite

MDPI and ACS Style

Usov, P.P. EHD Effects in Lubricated Journal Bearing. Lubricants 2018, 6, 12. https://doi.org/10.3390/lubricants6010012

AMA Style

Usov PP. EHD Effects in Lubricated Journal Bearing. Lubricants. 2018; 6(1):12. https://doi.org/10.3390/lubricants6010012

Chicago/Turabian Style

Usov, Pavel Pavlovich. 2018. "EHD Effects in Lubricated Journal Bearing" Lubricants 6, no. 1: 12. https://doi.org/10.3390/lubricants6010012

APA Style

Usov, P. P. (2018). EHD Effects in Lubricated Journal Bearing. Lubricants, 6(1), 12. https://doi.org/10.3390/lubricants6010012

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop