Minimizing Misalignment Effects in Finite Length Journal Bearings

: This paper focuses on a method to reduce the detrimental effects that occur due to the misalignment in journal bearings by approaching it with the more complete model of a ﬁnite length bearing. Such a drawback is quite common in industrial applications, and it is generally accepted that misalignment causes a signiﬁcant thinning in the ﬁlm thickness in the area that is close to the bearing edges. Therefore, removing a certain volume of material from the inner surface of the bearing (bushing) over a distance that is at the bearing edges provides an additional clearance to compensate for the clearance reduction that is due to misalignment. A numerical solution that is used in this work is based on the ﬁnite difference method where the Reynolds boundary conditions are considered in the solution scheme, thereby, using an iterative procedure to identify the cavitation zone. A three-dimensional misalignment model is incorporated in the solution in order to provide a more realistic presentation of the deviations and errors that there are in comparison with the ideal aligned case. It has been found in the present work that the edge modiﬁcation increases the thickness of the lubricant layer considerably and reduces the pressure spikes that are associated with the presence of misalignment. The suggested design also reduces the coefﬁcient of friction in comparison with that of the misaligned case. Furthermore, this method helps in reducing the asymmetry of the hydrodynamic pressure ﬁeld that results from the misalignment. This method enables the operation of journal bearings over a wider range of misalignment levels without sacriﬁcing the load-carrying capacity of the bearing by maintaining a relatively thicker layer of lubricant at the critical positions that are not so due to the effects of misalignment.


Introduction
Journal bearings are widely used in industrial uses which represent an essential part of many applications such as generators, turbines, pumps and compressors.The shaft rotates inside the bearing where the clearance is normally in the order of tens of microns, and which is filled with a useful lubricant.Such relatively thin clearance provides limited tolerances against installation errors, manufacture errors and shaft deformation under load.Misalignment is the common term that is used when the journal bearing operates under these conditions.In general, a misalignment in journal bearing is unavoidable but it is possible, to some extent, to limit its consequences on the system's performance.Extensive works have been published to explain the misalignment effects on journal bearings characteristics.Ebrat et al. in 2004 [1] studied the dynamic characteristics of journal bearings considering their misalignment, where they emphasized that the effect of shaft misalignment is very important in assessing the rotor dynamic behavior and the interaction of the crankshaft and block of the internal combustion engine.Sun and Gui in 2004 [2] found that a misalignment due to shaft deformation has significant effects on the film thickness and the pressure distributions in the full film lubrication of a bearing.Jang and Khonsari in 2015 [3] explained that a misalignment is detrimental to the overall performance of the bearings.Jamali and Al-Hamood in 2018 [4] studied the 3D misalignment effects in journal bearings, where they found that misalignment causes severe thinning in the film thickness and significant increasing in the pressure levels.Recently, Song et al. in 2022 [5] showed that the friction increases obviously due to a misalignment in the mixed lubrication area.In the mixed lubrication regime of a lubricated contact, the surface roughness significantly affects the performance of the contact.Mixed lubrication, which is also known as the partial lubrication regime, has significant considerations in a wide range of applications such as the internal combustion engine.The operation of many components is usually classified to be under this regime of lubrication.In this case, the applied load is supported by the surface asperities and the fluid film of the lubricant which results in metal-to-metal contact as well as elastohydrodynamic lubrication regime.In this regime (mixed lubrication), the calculation of the coefficient of friction is relatively difficult as a result of the continuous changing of the instantaneous surface topography over the contact zone due to the movement of the surfaces.Therefore, the solution of the hydrodynamic lubrication regime problems is relatively simple.However, solving any contact problem that is based on the mixed lubrication regime is related to the ratio of the film thickness to the combined surface roughness [6].
The friction (as well as wear) in journal bearings has been under investigation by researchers for decades due to its significant effects on the performance of this type of bearing.Muskat and Morgan [7] showed that the friction coefficient for finite length journal bearings is greater than those of infinite length bearings for fixed Sommerfeld variables and it is in direct proportion with the increasing of the bearing's length.Dufrane et al. [8] proposed a model to study the wear effect on hydrodynamic lubrication.They found that there is an optimum value for the film thickness that is related to the progress of wear in the bearings.Unlu et al. [9] used a new approach that was based on experimental and artificial neural network approaches to determine the friction coefficients in radial bearings.Dry and lubricated conditions were examined under different values of velocity and load.The misalignment effects on the rate of wear progress in journal bearings have been examined numerically by researchers using finite difference or finite element methods.Goenka [10] used a numerical solution for the analysis of the journal bearings.His model can be used to analyze the bearings when they are under misalignment conditions without needing to employ a high-computing cost for both partial and full arc journal bearings.Nikolakopoulos and Papadopoulos [11] analyzed misaligned journal bearings based on a numerical solution for the Reynolds equation.Both of the linear and nonlinear characteristics of journal bearings were presented.Bouyer and Fillon [12] presented a more significant experimental study to evaluate the performance of journal bearings (100 mm diameter) that were under misalignment torque.Their results showed that the bearings' performance is greatly affected by this applied misalignment torque.For example, the minimum film thickness is decreased by 80% due to misalignment.Boedo and Booker [13] also investigated, numerically, the steady-state and transient behaviors of misaligned bearings.Sun et al. [14] used a special test bench to examine the journal bearings' performance considering the journal misalignment due to the shaft bending under load.They found that there were changes in the values and distributions of the pressure field and film thickness as a result of the presence of misalignment.Padelis et al. [15] developed a model to determine the relationships among misalignments, the coefficients of friction and wear depths in journal bearings.They used a numerical solution to solve Reynolds equation that was based on the finite element method.The presented functions are related to the friction coefficient with the misalignment conditions and wear depth for different values of the Sommerfeld number.
The optimization of the bearing's geometry was investigated experimentally by Nacy in 1997 [16] where the bearing was chamfered at the edges in order to control the bearing side leakage.Also, Bouyer and Fillon in 2004 [17] used defects on the geometry of the bearing in order to investigate the bearing's performance when it was under misalignment torque.Strzelecki in 2005 [18] modified the bearing over the whole length of it by creating a hyperboloidal profile where he found that such a profile is useful in carrying extreme loads under misalignment.Chasalevris and Dohnal in 2016 [19] used journal bearings of variable geometry to improve the stability of the system.More recently Ren et al. in 2021 [20] studied numerically the effects of the axial profile parameters on the bearing's performance, where their results showed that using quadratic profile significantly improves the bearing's performance.
Booker et al. [21] presented an interesting review for an engine-bearing design that was based on conformal elastohydrodynamic analyses.In their review, physical models and solution methods were examined.Allmaier and Offner [22] also reviewed the elastohydrodynamic simulation of journal bearings.They addressed a wide range of topics such as polymer coatings, the mixed lubrication regime, the simulation of the thermal behavior and a low-viscosity lubricant.This review emphasized that journal bearings are still facing new challenges which motivate the researcher to develop more accurate and new methods to address the previously mentioned topics.
Recently, Krishnkant et al. [23] explained that shaft misalignment is considered as a major issue that hampers the satisfactory operation of the bearing system, and perfect alignment is rather difficult to achieve.In practical applications of the journal bearing, there is always some amount of misalignment that exists.Song et al. [5] established a mixed lubrication model for misaligned bearings while considering the effects of turbulence and cavitation using a finite difference scheme.Their results showed that the misalignment increases the friction obviously at the mixed lubrication area.Guo et al. [24] used profile modification for journal bearings while considering the misalignment effects in order to minimize the bearing wear due to the journal defection under the external load.They concluded that bearing modification provides an engineering approach for the anti-wear design of this type of bearing.Also, there are many other research papers that have investigated the aforementioned misalignment effects, but they have done so marginally and without an in-depth study and analysis of the problem [25][26][27].
This work aims at improving the geometrical design of journal bearings that are working under misalignment in more than one plane.It has been found that, despite the extensive work that is available in the literature about journal bearings, this topic (geometrical design) still requires more investigation in order to minimize the consequences of the severe levels of misalignment.A wide range of misalignment deviations as well as bearing modification parameters in both the longitudinal and radial directions are considered in this work.Such a comprehensive study is necessary to identify the optimum position for the modification along the bearing's length as well as the depth of modifications in the radial direction.This study is achieved numerically using a finite difference method where the Reynolds boundary conditions are used in the solution scheme for the finite length bearings.The results are presented in a dimensionless form to serve the general applications for the outcome of this work.

Mathematical Models and Numerical Solution
The problem of journal bearings is governed by the Reynolds equation in additional to the film thickness equation which is essentially a geometrical relationship.For the steady-state case and incompressible flow, the Reynolds equation is given by [28]: where U: mean velocity, p: pressure, h: film thickness and η: viscosity.The geometry of the journal bearing is shown in Figure 1. Figure 1a illustrates the coordinates, attitude angle (∅), the angle θ, which is measured counter-clockwise from the positive y-axis, and the eccentricity (e), which is the distance between the bearing center and the shaft center.Figure 1b shows a 3D representation of the bearing.
where U: mean velocity, p: pressure, h: film thickness and : viscosity.The geometry of the journal bearing is shown in Figure 1. Figure 1a illustrates the coordinates, attitude angle (∅), the angle , which is measured counter-clockwise from the positive y-axis, and the eccentricity (e), which is the distance between the bearing center and the shaft center.Figure 1b shows a 3D representation of the bearing.The film thickness is calculated by [29]: where  = is the eccentricity ratio and C is the radial clearance between the bearing and the shaft.
Using dimensionless variables, these equations become: where: The film thickness is calculated by [29]: where ε r = e C is the eccentricity ratio and C is the radial clearance between the bearing and the shaft.

Misalignment Model
The model of misalignment that is used in this paper is essentially adopted from a previous work of one of the authors of this paper [4].This reference provided a detailed explanation for the misalignment model for its solution.Figure 1c shows an exaggerated schematic drawing of the 3D misalignment model.Using the 3D misalignment model in the solution scheme provides a more realistic presentation of the misalignment case where the deviations and errors that are in comparison with the ideal, aligned case can be considered in the analyses.The dimensionless deviations along the bearing's length in the horizontal as well as the vertical directions are given by: where δ is the dimensional value of the misalignment, ∆ = δ/c is the dimensionless misalignment and Z = z/L.
It is worth mentioning that the attitude angle in the case of misaligned bearing is a function of the z position which can be given by [4]: e cos ∅−δv(z) .for z > L/2

Bearing Modification
The new shape of the inner surface at the two bearing ends is schematically illustrated in the Figure 2. The resulting gap in the circumferential direction, which is a function of the position z along the bearing's length, is given by:

Coefficient of Friction
The coefficient of friction is also calculated in this work in order to investigate the effect of bush chamfering on this important contact characteristic.The friction coefficient that is calculated in this work is based on the model which is given by Lund and Thomsen [30] which is, The dimensionless parameters A and B in Equations ( 6)-( 8) are: A = a/C and B = b/L.The values of A and B are used to investigate the adopted method of removing the material from the bearing's inner surface.Scaling the design parameters to the clearance, C and the bearing length, L gives us a clearer picture about the amount of change that occurs in the bearing's geometry.
The coupling of the Equations ( 4)-( 8) provides the gap between the misaligned shaft and the modified bearing.

Coefficient of Friction
The coefficient of friction is also calculated in this work in order to investigate the effect of bush chamfering on this important contact characteristic.The friction coefficient that is calculated in this work is based on the model which is given by Lund and Thomsen [30] which is, where W is the total load, the friction force is F f = P u , u = Rω, W and the power loss (P) is: where F X and F Y are the bearing forces in the horizontal and vertical direction, respectively, which can be easily calculated by the integration of the resulting pressure field.

Numerical Solution
The analytic solution of the hydrodynamic problem in journal bearings is possible for particular cases, where the approximate solutions are obtained.These approximate solutions involve ignoring one of the pressure gradients in the Reynolds equation.The pressure gradient in the circumference direction (∂p/∂θ) is ignored in the case of short bearings, and the gradient in the longitudinal direction (∂p/∂θ) is ignored in the assumption of long bearings.These assumptions give acceptable results for a certain length-to-diameter ratio (L/D) of the bearing.However, for finite length bearings where neither of the pressure gradients can be ignored, a numerical solution is required.In this work, the governing equations are solved numerically using the finite difference method.This solution is explained in reference [4].The solution space is divided to N × M points, where N and M are the numbers of points in the circumferential and longitudinal direction of the bearing, respectively.Part of the solution space is shown in Figure 3 where point (i, j) and its surrounding four mesh points are illustrated.of the pressure gradients can be ignored, a numerical solution is required.In this work, the governing equations are solved numerically using the finite difference method.This solution is explained in reference [4].The solution space is divided to  ×  points, where  and  are the numbers of points in the circumferential and longitudinal direction of the bearing, respectively.Part of the solution space is shown in Figure 3 where point (, ) and its surrounding four mesh points are illustrated.By using the central difference to discretize the pressure gradient as well as the variation of the film thickness in the circumferential direction, it gives us, ∂P ∂θ a = The film thickness at the boundaries of the control volume that shown in Figure 3 is used as an average value which can be given by: Similarly, the gradient in the z direction is discretized as follows: By substituting Equations ( 10)- (20) into Equation ( 2), the pressure at point i, j can be given by: where ∆Z and ∆θ are the mesh steps.The oil film thickness in terms of the i, j position is, Designs 2022, 6, 85 8 of 17

Solution Procedure
Equation ( 21) is solved at all points on the solution domain except for the boundaries where the pressure values are set to zero.This step requires a first assumption of the pressure field over this solution domain.The new calculated pressure values are then used as the input values to calculate another pressure field using Equation (21).The pressure distribution in the solution space is obtained iteratively by using a successive over relaxation method to reduce the required number of iterations.The relaxation factor (RF) that is used in this case is 1 < RF < 2 with there being an optimum value of 1.93.The calculated pressure (P new ) that is based on Equation (11) is updated using the following equation: where P old is the pressure value from the previous iteration at the i, j node.
It is worth mentioning that the overall convergence of the solution requires the convergence of the pressure distribution which is given by ∑|P (i,j)new −P (i,j)old | ∑ P (i,j)old < 10 −7 in addition to the convergence of the load which is explained below.
After the convergence of the solution is achieved, the pressure values are available at all positions in the solution domain.Therefore, the bearing forces in the radial and tangential forces can be easily calculated numerically, using the following discrete forms in order to calculated the total load W : The calculated load is compared with the input load to achieve an accuracy limit that is ∓10 −5 relative to the input load.If the calculated load is not within this tolerance limit, then the eccentricity ratio is continuously updated until the load convergence is obtained.
Equations ( 5)-( 8) are also easily discretized and combined with Equation ( 10) to obtain the total gap between the surfaces at any i and j node.The Reynolds boundary conditions are used in the solution scheme using an iterative procedure to identify the cavitation zone where p = 0 and dp dθ = 0 are satisfied at the limit of this zone.It is worth mentioning that the iterations for the misaligned and modified design cases (under overrelaxed solution) are continued until the results converge to the same load in all of the cases.This load corresponds to an eccentricity ratio of 0.6 in the aligned case.A flow chart for the solution scheme is shown in Figure 4. cases.This load corresponds to an eccentricity ratio of 0.6 in the aligned case.A flow chart for the solution scheme is shown in Figure 4.

Results and Discussions
The numerical solution is first performed to examine the required mesh size in both directions that is sufficient enough to obtain independent results on the number of nodes in the solution space.The effect of misalignment on the pressure distribution is shown in Figure 5 for a finite length bearing where / = 1.25.Figure 5a shows the case when ∆ℎ = ∆ = 0, which represents a perfectly aligned bearing.It can be seen how the pressure distribution is uniform and symmetric about the line  = 0.5 (half of the bearing length).Figure 5b-d shows the results for the cases ∆ℎ = ∆ = 0.4 , ∆ℎ = ∆ = 0.55 and ∆ℎ = ∆ = 0 .58,respectively.The case of ∆ℎ = ∆ = 0.4 is chosen as the

Results and Discussions
The numerical solution is first performed to examine the required mesh size in both directions that is sufficient enough to obtain independent results on the number of nodes in the solution space.The effect of misalignment on the pressure distribution is shown in Figure 5 for a finite length bearing where L/D = 1.25.Figure 5a shows the case when ∆h o = ∆v o = 0, which represents a perfectly aligned bearing.It can be seen how the pressure distribution is uniform and symmetric about the line Z = 0.5 (half of the bearing length).Figure 5b-d  The influences of bearing edge chamfering on the pressure distributions are shown in Figure 6.The result of the misaligned case (without modification) that is presented in previous figure is repeated here for the purpose of comparison, as shown in Figure 6a, where ∆ℎ = ∆ = 0.58. Figure 6b-d shows the effect of the modification on the pressure distribution where  =  = 0.1,  =  = 0.2 and  =  = 0.25, respectively.These values are selected to be presented in this figure after implementing a wide range of modification parameters  and  in the solution scheme.The influences of bearing edge chamfering on the pressure distributions are shown in Figure 6.The result of the misaligned case (without modification) that is presented in previous figure is repeated here for the purpose of comparison, as shown in Figure 6a, where ∆h o = ∆v o = 0.58. Figure 6b-d shows the effect of the modification on the pressure distribution where A = B = 0.1, A = B = 0.2 and A = B = 0.25, respectively.These values are selected to be presented in this figure after implementing a wide range of modification parameters A and B in the solution scheme.
It can be seen from this figure that the modification returns to some degree the symmetrical pressure distribution despite the presence of a misalignment.This outcome is expected to have significant effects on the dynamic behavior of the system, particularly when it is under dynamic load such as impact load and position perturbation, which will be addressed comprehensibly in future work.
The corresponding film thickness distributions are shown in Figure 7. Figure 7a shows the case of the aligned bearing, while Figure 7b shows the results of the misaligned case where ∆h o = ∆v o = 0.58. Figure 7c illustrates the effect of modification on the film thickness of a misaligned bearing where ∆h o = ∆v o = 0.58 and A = B = 0.25.These figures explain how the misalignment disturbs the perfectly extruded shape of film thickness along the bearing's length which results in a reduction in the film thickness and a concentration of the pressure spikes that are close to the film thinning position which is related to the misalignment.It can be seen from this figure that the modification returns to some degree the symmetrical pressure distribution despite the presence of a misalignment.This outcome is expected to have significant effects on the dynamic behavior of the system, particularly when it is under dynamic load such as impact load and position perturbation, which will be addressed comprehensibly in future work.
The corresponding film thickness distributions are shown in Figure 7. Figure 7a shows the case of the aligned bearing, while Figure 7b shows the results of the misaligned case where ∆ℎ = ∆ = 0.58. Figure 7c illustrates the effect of modification on the film thickness of a misaligned bearing where ∆ℎ = ∆ = 0.58 and  =  = 0.25.These figures explain how the misalignment disturbs the perfectly extruded shape of film thickness along the bearing's length which results in a reduction in the film thickness and a concentration of the pressure spikes that are close to the film thinning position which is related to the misalignment.Figure 8 shows a more detailed comparison for the previous results at the section of maximum pressure in the circumferential direction.Figure 8a shows the results of the pressure variation, whereas Figure 8b shows the corresponding results about the film thickness for three cases of a perfectly aligned bearing, a misaligned bearing (∆h o = ∆v o = 0.58) and a modified misaligned bearing (where A = B = 0.25).It is clear from these figures that removing material from the bearing's inner surface at the edges brings back the shape of both the pressure and film thickness at a level that is very close to the normal case of an aligned bearing at the position of maximum pressure.Figure 8 shows a more detailed comparison for the previous results at the section of maximum pressure in the circumferential direction.Figure 8a shows the results of the pressure variation, whereas Figure 8b shows the corresponding results about the film thickness for three cases of a perfectly aligned bearing, a misaligned bearing (∆ℎ = ∆ = 0.58) and a modified misaligned bearing (where  =  = 0.25).It is clear from these figures that removing material from the bearing's inner surface at the edges brings back the shape of both the pressure and film thickness at a level that is very close to the normal case of an aligned bearing at the position of maximum pressure.
The values of  and  for the cases that are shown in Figures 5 and 6 are shown in Figure 9a,b, respectively.The maximum pressure increases by 51.63 %, in comparison with the aligned case when ∆ℎ = ∆ = 0.58 and the corresponding film thickness reduces by 83.5%.The bearing chamfering reduces the  by 26.38% and increases the  to be 0.254, which was 0.066 (an increase of 284.85%) when  =  = 0.25 for the misaligned bearing.This important outcome is expected to have positive consequences on the bearing life through the limitation of the occurrence of partial contact between the asperities (mixed regime) of the shaft and the journal bushing when a very thin layer of oil is interposed in the resulting wedge.More details about the effect of the modification parameters on the maximum pressure and the minimum film thickness of the misaligned bearings are shown in Table 1.In this table, two extreme misalignment cases that were under wide range of chamfer parameters are investigated.The misalignments in these cases have the following values: ∆ℎ = ∆ = 0.55 and ∆ℎ = ∆ = 0.58.On the other hand, the modification parameters are  =  = 0 (Unmodified bearings),  =  = 0.1,  =  = 0.2,  =  = 0.25 and  =  = 0.3.It can be seen that in both cases, the modification reduces the maximum pressure values and increases the minimum film thickness when the chamfer parameters are less than 0.3.When  =  = 0.3, the minimum film thickness continues to increase but, on the other hand, the maximum pressure value also increases.This can be attributed to the fact that as the parameter  increases (modification parameter in the radial direction), the modification part of the bearing will be less involved in supporting the load.This can be explained in other words; that the pressure will generate over a smaller part of the modified portion of the bearing.Despite that the pressure values are slightly increased when  =  = 0.3, higher values for the chamfer parameters are not recommended.

Mis.
Mis. and mod.More details about the effect of the modification parameters on the maximum pressure and the minimum film thickness of the misaligned bearings are shown in Table 1.In this table, two extreme misalignment cases that were under wide range of chamfer parameters are investigated.The misalignments in these cases have the following values: ∆h o = ∆v o = 0.55 and ∆h o = ∆v o = 0.58.On the other hand, the modification parameters are A = B = 0 (Unmodified bearings), A = B = 0.1, A = B = 0.2, A = B = 0.25 and A = B = 0.3.It can be seen that in both cases, the modification reduces the maximum pressure values and increases the minimum film thickness when the chamfer parameters are less than 0.3.When A = B = 0.3, the minimum film thickness continues to increase but, on the other hand, the maximum pressure value also increases.This can be attributed to the fact that as the parameter A increases (modification parameter in the radial direction), the modification part of the bearing will be less involved in supporting the load.This can be explained in other words; that the pressure will generate over a smaller part of the modified portion of the bearing.Despite that the pressure values are slightly increased when A = B = 0.3, higher values for the chamfer parameters are not recommended.the previous figures, where the film thickness is reduced significantly due to the presence of a misalignment.It can be seen that in the first case, the coefficient of friction ( f ) is 2.218 and in the second case, the misalignment increases the friction coefficient by 14.56%.On the other hand, chamfering the bush edges for the misaligned case reduces the coefficient of friction by 9.76%, in comparison to that which is achieved with the misaligned case.However, as the misalignment reduces, the thickness of the lubricant layer that is close to the bush edges increases.The surface roughness of the shaft and the bush will certainly have its contribution to the friction between the two surfaces.In such case, the coefficient of friction in the misaligned case is clearly more than the calculated value is as based on the approach that adopted in this work.Therefore, a more realistic analysis is required to consider the surface features that are at least over the region where the film thickness is relatively thin.Such a solution requires an elastohydrodynamic analysis of the contact problem in order to incorporate the deformation of the rough surfaces.Nevertheless, chamfering the bush edges increases the minimum film thickness significantly, as illustrated previously (see Figure 9).This is highly expected to minimize the surface features' effects on the characteristics of this contact problem.However, the elastohydrodynamic analysis of the contact problem of the journal bearing while considering the surface features that are under misalignment as well as bush chamfering will be studied by the authors in future works.

Conclusions
In this work, the effects of chamfering the bearing inner surface on the pressure and film thickness distributions are presented in detail.A 3D misalignment model is incorporated in the analysis scheme where a numerical solution is used to solve the problem of the finite length journal bearing.The solution procedure involves the use of the Reynolds boundary conditions, where the limit of the cavitation zone is determined iteratively.the results showed that such a design overcomes the problem of the thinning of the lubricant layer that is due to misalignment by increasing the film thickness significantly (almost three times) and helping to reduce the pressure spike that result due to misalignment by 26.38%.The coefficient of friction which is one of the most important characteristics of this type of bearing is also reduced by 9.76% due to introducing the chamfering of the bush edges in comparison with the misaligned case.Another important outcome has been found is the reduction of non-symmetry in the 3D pressure distribution which is expected to have marked influence on the dynamic response of the rotor-bearing system through its impact on the dynamic coefficients.

Figure 1 .
Figure 1.Schematic drawing of the journal-bearing.(a) Side view of aligned bearing, (b) 3D representation and (c) misalignment model used in this work.

Figure 1 .
Figure 1.Schematic drawing of the journal-bearing.(a) Side view of aligned bearing, (b) 3D representation and (c) misalignment model used in this work.

Figure 2 .
Figure 2. Bearing with modified edges.(a) 3D full bearing, (b) 3D section and (c) Schematic drawing of the bearing after removing material from the inner surface.

Figure 2 .
Figure 2. Bearing with modified edges.(a) 3D full bearing, (b) 3D section and (c) Schematic drawing of the bearing after removing material from the inner surface.

Figure 3 .Figure 3 .
Figure 3. Part of solution domain.By using the central difference to discretize the pressure gradient as well as the variation of the film thickness in the circumferential direction, it gives us, ( , ) ( , )

Figure 4 .
Figure 4. Flow chart for the solution procedure.

Figure 4 .
Figure 4. Flow chart for the solution procedure.

Figure 7 . 17 Figure 8 .
Figure 7. Dimensionless film thickness (L/D = 1.25).(a): Aligned (∆h o = ∆v o = 0), (b): misaligned (∆h o = ∆v o = 0.58) and (c): modified bearing (∆h o = ∆v o = 0.58 and A = B = 0.25 ).The values of P max and H min for the cases that are shown in Figures5 and 6are shown in Figure9a,b, respectively.The maximum pressure increases by 51.63 %, in comparison with the aligned case when ∆h o = ∆v o = 0.58 and the corresponding film thickness reduces by 83.5%.The bearing chamfering reduces the P max by 26.38% and increases the H min to be 0.254, which was 0.066 (an increase of 284.85%) when A = B = 0.25 for the misaligned bearing.This important outcome is expected to have positive consequences on the bearing life through the limitation of the occurrence of partial contact between the asperities (mixed regime) of the shaft and the journal bushing when a very thin layer of oil is interposed in the resulting wedge.

Figure 8 .
Figure 8.Effect of misalignment and profile modification on (a): pressure variation and (b): film thickness at the section of maximum pressure.

Table 2 .
The values for the coefficient of friction calculations for different cases.