{'fraction_non_alphanumeric': 0.034037742958520426, 'fraction_numerical': 0.012151851037774358, 'mean_word_length': 4.508128675198893, 'pattern_counts': {'":': 1, '<': 0, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 7, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
A New Method to Calculate Water Film Stiffness and Damping for Water Lubricated Bearing with Multiple Axial Grooves

Water lubricated guide bearings for hydro turbines and pumps are conventionally designed with multiple axial grooves to provide effectively cooling and flushing away abrasives. Due to the variety of groove configuration in terms of number and size, a predication of their performance is difficult. This paper deals with an analytical procedure to investigate groove effect on load capacity, stiffness and damping for this type of bearing where it is considered as an assembly of many inclined slide bearings. The result can be applied to bearings made of hard materials combined with low bearing pressure.


Introduction
As we all known, water-lubricated guide bearings for hydro turbines and pumps are conventionally designed with multi-axial grooves. These grooves are provided for purpose of effectively cooling the bearing and flushing away abrasives. However, due to the variety of groove design in terms of its number and size, a prediction of bearing performance in terms of load capacity, stiffness and damping characteristics is very difficult. The author of this paper [1] introduced an analytical method to investigate groove effect on the Sommerfeld Number and coefficients of stiffness and damping based on inclined slide bearing solutions for bearings with rigid surface. However, the quality and accuracy of the solution depends on how close the geometry of the inclined slide bearing to represent the actual wedge shape at individual pads of the grooved bearing, especially the wedge shape of those pads which are loaded most. This paper examined three different geometric shapes of inclined slide bearings and provided a solution that would have satisfactory accuracy.
A brief review of available literature is useful to understand the development of this method. For rigid surface plain bearings with no grooves, in terms of steady operation, besides the classic solution of long bearing theory by Sommerfeld [2] and short bearing theory by DuBois and Ocvirk [3], there are several excellent analytical solutions for finite length bearings. The finite length bearing theory by Childs et al. [4,5] is one of the excellent solutions. Another good analytical solution is proposed by Capone et al. [6]. Numerical solutions with using finite difference and finite element methods are abundant. The evaluation of them is not the focus of this paper. For bearings designed with multi-axial grooves, Pai et al. [7] published a number of works on steady performance and dynamic stability of simple rotor. Ren [8] published a paper on calculation of water film thickness of water lubricated bearing with multi-axial grooves for steady state operation. On stiffness and damping coefficients of non-grooved plain bearings, classical short bearing solution is the most popular one. The solution by Childs et al. [4,5] is strongly recommended for finite length bearings. It needs to mention that above works are related to rigid surface bearings. For deformable surface bearings, the effect of surface deformation is considered. Lahmar et al. [9] provides a procedure to simultaneously evaluate both static and dynamic performance with small perturbation method. Recent development was focused on CFD and FSI (fluid structure interaction) [10][11][12][13][14][15][16][17]. The effect of turbulence on stiffness and damping was investigated in Refs. [18,19]. In review on available information, methods to determine stiffness and damping are mainly relying on numerical simulation for bearings with multiple axial grooves. The objective of this paper was to provide a semi-analytical method to investigate the groove effect on load capacity, stiffness and damping based on infinite length and rigid surface which is approximately valid for bearings made of hard polymers under relatively low bearing pressure [10].
For the calculation results to be useful, certain conditions for bearings have to be applied. First of all, the ratio of bearing length to the width of bearing pad must be greater than 3.0 or higher. Secondly, the bearing pressure shall be relatively low so that the surface deformation effect doesn't overwhelmingly change the result. In spite of this paper doesn't include the effect of elastic deformation of bearing surface in terms of elastohydrodynamic lubrication, the results are valid approximately for polymer bearings with higher hardness and lower pressure that most pump and turbine guide bearings are the case. The result is not suitable to water lubricated bearings with rubber staves that needs special treatment either through experiment [20][21][22][23][24][25][26][27][28][29] or numerical analysis. Experimental study on bearings with multiple axial grooves demonstrates that a relatively rigid surface of bearing pad more easily forms hydrodynamic pressure than soft surface [22]. This has been demonstrated in the elastohydrodynamic study on sliding bearings [30]. A practice engineering application of stiffness and damping was shown in Ref. [31].

Stiffness and Damping of Inclined Slide Bearings with Different Geometries
The idea to evaluate the load capacity (Sommerfeld Number), the stiffness and damping coefficients of a circular bearing with multiple axial grooves is that the circular journal bearing can be considered as an assembly of many simple inclined slide bearings (Figure 1), so that the analytical results for an inclined slide bearing ( Figure 2) can be used as building blocks to form a calculation method. Without loss of generality, Figure 1 demonstrates pressure created by bearing pad only. In reality, water pressure in grooves varies from negligible to significant. Since this procedure works with non-dimensional functions, the water pressure in grooves can be easily added back to dimensional pressure. The implementation of this idea starts from evaluating the dynamic characteristics of the sliding bearing shown in Figure 2. To avoid of confusion, an inclined slider is also called inclined slide bearing. A sliding pad or bearing pads are always referring to the bearing surface between two neighboring grooves. As indicated in the introduction, the load capacity of the entire bearing depends on the load capacity of the individual inclined sliding pads. The accuracy of the solution is a direct function of how the geometry of the inclined sliding bearing to represent the wedge shape of the individual pads of the circular bearing. In following sections, three useful geometries of inclined slide bearing are examined, the proximity to the wedge shape of the main bearing is compared.

Linear Inclined Slide Bearing
The linear inclined slider was used in previous work [1]. The water film is presented with a linear function as follows: The main functions of the solution are: (2) Grooved bearing as an assembly of sliding pads

Damping function (non-dimensional)
It is noticed that the right side of Eqs.
(2)-(5) is a function of the ratio of film thickness at leading edge to the film thickness at trailing edge only. For the purpose of evaluating stiffness and damping, the film thickness at leading and trailing edges is considered as a function of time. This is shown in Figure 3.
The non-dimensional dynamic load is expressed by stiffness and damping function: According to Ref.
[1], the final load on the linear inclined slide bearing is expressed with Eq. (7) is the force of the individual slider under dynamic motion. It is the fundamental relationship between the slider force and slider displacement and squeeze velocity. As long as the four functions expressed from Eqs. (2)-(5) are known, the load capacity of the inclined slide bearing is fully defined. Therefore, the whole subject turns to finding a set of functions expressed in Eqs.

Exponential Gap Slide Bearing
The shape of exponential gap slide bearing is expressed with following equation: The four functions of result are as follows: 1. Non-dimensional load capacity 2. Location of static load center (non-dimensional distance) 3. Stiffness function (non-dimensional) 4. Damping function (non-dimensional)

Parabolic Gap Slide Bearing
The author of this paper derived the four functions for a parabolic gap slider (see Appendix). The film thickness expression is presented by: The resulting functions are as follows.
1. Non-dimensional load capacity  Figure 4 shows the basic functions for these three inclined slide bearings. Since the non-dimensional load capacity function is exactly the half of the stiffness function for all type of geometries, it was not repeated in Figure 4.
It is to notice that the stiffness function of the parabolic gap and exponential gap is higher than the linear gap for all leading to trailing edge film ratios. This implied if the parabolic function is in good agreement with the actual bearing clearance (pads), a circular bearing simulated with the parabolic function will have a higher load capacity. The dynamic load center is slightly different from static load center Figure 4(d). Appendix provides the definition of load center ratio R(η). This paper used static load center for Sommerfeld number evaluation and dynamic load center for stiffness and damping evaluation for all three types of sliding bearings.

Assembly Procedure
To evaluate which type of slider geometry best suitable for building the circular bearing, the assembly procedure must be presented first. The first step of the assembly procedure is to define the location angles of each pad relative to a rotating co-ordinate frame r-φ [8]. Figure 5 shows a circular journal bearing with multiple axial grooves under a steady operational condition. By given load and shaft speed, the shaft center is offset from bearing center with an eccentricity of "e". The connecting line between bearing center and shaft center is in-line with r-axis of the rotating coordinate system r-φ. Assuming the bearing is fixed in position and the load is vertical as shown on the figure, the r-φ coordinate system has an attitude angle "Φ" with respect to the loading direction, namely, the y-axis. The attitude angle changes depending on load, shaft speed and groove numbers.
A second co-ordinate frame x-y is defined in line with load direction. In this system, y-axis is in load direction and x-axis is perpendicular to the load direction.
In Figure 5, the r-axis divides the entire bearing into two equal halves. All pads underneath r-axis (in the sense of Figure 5) have convergent angles in shaft rotating direction and therefore are able to create hydrodynamic lifting forces. All pads above r-axis have divergent angles in shaft rotating direction and are therefore not able to create hydrodynamic lifting forces. In theory, the divergent bearing half could create a vacuum and therefore a negative pressure. However, since in practice, almost in all the cases, outside source of lubricant will be supplied to the bearing grooves, the divergent half will not create a negative pressure, but keep the same level of pressure as supplied lubricant as in grooves. It is therefore acceptable to assume the pressure on that half of bearing as zero.
This assumption is corresponding to half Sommerfeld or Gumbel boundary condition.
Since the bearing is assumed to be fixed, all angles ( α Li , α Ti = 1, 2, 3, …, N/2) defining the positions of grooves will change with attitude angle which is an unknown parameter. One set of groove location angles only defines a particular equilibrium of steady operation. For calculation purpose, a set of "floating numbers" are assigned to the pads underneath the r-axis. As a rule, no matter how the attitude angle to change, it is always the first convergent pad underneath r-axis at the minimum water film location is assigned number "1". Other pads are enumerated clockwise with number 2, 3, 4,… in sequence.
After having defined the pad location angles, the film thickness ratio of leading to trailing edge under steady operation condition is The film thickness at any leading and trailing edge is The second step of the assembly procedure is to calculate the force contribution of each pad to support the entire bearing load. Figure 6 illustrates the supporting force from one pad with location angles α Li and α Ti . Considering Eqs. (19) and (20), the force by each individual pad can be calculated with using Eq. (7) obtained from previous section. All terms in Eq. (21) are referred to pad number "i": Figure 6, the pad load is considered to be in direction pointing to bearing center. The projection of bearing load to r-φ co-ordinate system is then (18)  In Eqs. (22), (23), the function K (η) is one of the function K L (η) , K E (η) or K P (η) depending on which one is chosen. The same applied to function C(η) . For the entire circular bearing, the dynamic part of pad force is only caused by a very small change of bearing eccentricity Δe and attitude angle ΔΦ, the film thickness at trailing edge can be correlated to these small eccentricity and attitude angle change as follows The same can be applied to the pad velocity. However, the velocity must refer to the entire pad, not just the (24) �h Ti = cos α Ti · �e + sin α Ti · e · �Φ.
trailing edge. This implies that the pad has no rotation about its load center at any instant of dynamic motion. Therefore the velocity becomes: The function A(η) is one of the A L (η), A E (η) or A P (η) depending on which type of sliding bearing chosen and A d (η) is one of the A Ld (η), A Ed (η) or A Pd (η) depending on which type of sliding bearing chosen. Inserting Eqs. (24) and (25) into Eqs. (22) and (23), their matrix form can be expressed as The coefficients of stiffness and damping for the pad with index "i" are:

Quality Comparison of Slide Bearing Geometry
Above section proposed an idea of using an array of inclined slide bearings to build a circular journal bearing with multiple axial grooves. However, the quality of this approach depends on how closely the individual slider represents the shape of individual pad at any given location and eccentricity ratio.
In true sense of grooved bearing, the non-dimensional form of water film thickness at any given location angle "α" (Figure 7) is When evaluating the individual inclined slide bearing, the non-dimensional water film thickness is expressed as a function of the film thickness ratio of leading edge to trailing edge as well as the local coordinate "s" (Figure 7). Therefore, the non-dimensional water film thickness for the pad number "i" in terms of above mentioned water film thickness ratio and local coordinate can be presented in following form: Where s = s/B is non-dimensional local coordinate. The true bearing may have grooves with chamfer or round fillet. Since any chamfer and fillet will be too large for water film formation. Therefore, chamfer and fillet must be considered as part of grooves, not surface of pads. All other expression of water film thickness for bearing pad "i" can also expressed with local coordinate and the water film thickness ratio at leading and trailing edge.

For linear inclined slide bearing
For exponential inclined slide bearing For parabolic inclined slide bearing Eqs. (29)-(31) presents the water film thickness that are intended to be used to replace Eq. (28). The purpose of doing so is to simplify the problem by solving Reynold's Equation at the inclined slide bearing level rather than at the full bearing level. It is to notice that all these equations return the film thickness at trailing edge of pad "i" which is (1 + ε · cos α Ti ) when s = 0 . By the same token, they return the film thickness at leading edge which is (1 + ε · cos α Ti ) · η i when s = −1 . To evaluate which function from Eqs. (29)-(31) is the best approach to Eq. (28), one set of square root errors was defined. These are: Eqs. (32)-(34) were derived from non-dimensional water film thickness and are function of eccentricity ratio, number and distribution of grooves. They are valid for any size bearings with any number of grooves. In following, a 12-groove bearing is examined for the square root errors. This bearing will have six bearing pads taking load. Assuming the location of the minimum film thickness falls into the very center of a groove, so that the first pad named with one will have the entire pad being loaded. Figure 8 shows the result for the bearing with 12 grooves. It showed that the parabolic inclined slider has the least error for the first pad for eccentricity ratio from 0.9 to 0.999 which is the range of most interest for water lubricated guide bearings. The linear slider seems to be best suitable for second pad at high eccentricity ratio and rest of other pads. The exponential slider seems to be suitable for pads except for first one at lower eccentricity ratio. However, this is only observations from a bearing with 12 grooves. Investigation on different number of grooves showed that the error for exponential slider changes rapidly with increasing eccentricity ratio, meaning that for small eccentricity ratio, errors are small, for large eccentricity ratio, errors are big. This is especially true for first and second pad. The linear slide bearing is insensitive to eccentricity ratio. Therefore, in following evaluation, a scheme that parabolic slider to first pad and (32) linear slider to rest of pads is applied. This paper investigated only three types of sliding bearings. It is certainly there must be other types of sliding bearings that would fit for the purpose. As demonstrated in Ref. [8], the first pad takes the most load of entire bearing. Even though the parabolic slider is only able to simulate the first loaded pad, it is still a significant improvement of load capacity in comparison to the approach with all linear sliders.

Effect of Groove Number
Ref. [8] provides a procedure to calculate the load capacity and water film under steady operational condition. The non-dimensional load capacity of bearing was determined as follows: where The total non-dimensional supporting force contributed by all pads underneath r-axis is then the sum of all components above: where λ is a number less than 1.0. In Figure 5, if the position of minimum film thickness is located within the pad number 1, only a part of this pad will take load. The number "λ" gives the percentage of the pad that takes load. Its value is unknown at beginning of a calculation. It depends on attitude angle "Φ", width of bearing pads "B", width of grooves as well as the relationship between bearing loading direction and the pad position to which the load pointing to. For vertical bearings, such as hydro turbine guide bearing and vertical pump bearing, loading direction is undefined. In this case, practical calculation can be done by assuming λ = 0.5 which presents a condition with least load capacity. For horizontal bearings, load direction can be easily defined. Parameter λ can be determined by an iteration. At beginning, first to assume an initial value, for example 0.5, then calculate the attitude angle according to Eq. (40) below and subsequently calculate the new λ-value. With new value run calculation again, get another attitude angle. Repeating the same procedure until a satisfactory result obtained. In evaluation on Figures 9, 10, 11 and 12, λ = 1.0 was used. A full analysis of effect of λ-value on Sommerfeld Number is worth of a full separate paper to discuss.
The resultant force in dimensional form is then The attitude angle is: According to conventional definition of Sommerfeld Number for circular bearings, it is The Sommerfeld Number for grooved bearings can be derived with using Eq. (39): The Sommerfeld Number, defined by Eq. (42) is a function of number of grooves and eccentricity ratio. Its reciprocal defines the load capacity of a bearing while the reciprocal of Eq. (41) is the actual load to the bearing. It is evident that for low eccentricity ratio (less than 0.9), there is no significant difference between the modeling with all linear sliders and the one with mixed sliders, namely, the first one using parabolic slider and rest of them linear sliders. However, for high eccentricity ratio (ε > 0.9), this reflects the case of high loading, the difference can be significant, Figure 9. The Sommerfeld Number simulated with all linear pads can be up to 5 times higher than the Sommerfeld Number with mixed pad geometry for the example investigated. By definition, higher Sommerfeld Number means lower load capacity. This is reflected in Figure 10.
The Sommerfeld Number ratio shown in Figure 10 is the Sommerfeld Number with all linear sliders divided by the Sommerfeld Number with mixed pads in which the first pad is parabolic and rest of them linear. Water lubricated guide bearings can be subject to eccentricity ratio as high as 0.999. In this case, an all linear modeling definitely under estimate the bearing loading capacity. Based on the quality comparison of sliding bearing geometries (39) in previous section, the parabolic gap is more closed to true shape of bearing clearance of the first pad. Therefore, the mixed scheme must more closely present the true bearing performance. It is an improvement of all linear slider modeling, especially for large eccentricity ratio. Ren et al. [1] in their previous paper quantitatively demonstrated that the load capacity of grooved bearing is lower than that of non-grooved bearings. The Sommerfeld Number of grooved bearing modeled with all linear pads was compared with the bench mark Sommerfeld Number, namely the formulation from Refs. [4,5]. A similar comparison is made here for the grooved bearing with mixed type of inclined sliders to the non-grooved bearings. The Sommerfeld Number of the solution by Childs is again used as bench mark for the comparison.
The massive curve on Figure 11 is the Sommerfeld Number of non-grooved bearing. Other curves are for grooved bearing with different number of grooves. Definitely, the grooves reduce the load capacity of a bearing. The more grooves, the bigger is the load capacity reduction.
The ratio of Sommerfeld Number for grooved to nongrooved bearing is shown in Figure 12. This provides a better visualization as how much the load capacity reduction can be expected. From Figure 12 is to see reducing the number of grooves is an effective way to increase load capacity of a grooved bearing.

Effect of Groove Size
From the defining equation of Sommerfeld Number of grooved bearing Eq. (42), it is a function of ratio d/B. For a fixed number of grooves, the size (width) of groove will take away a part of bearing surface which results in a narrower bearing pad. This increases the d/B ratio. Therefore, for a real load capacity of a practice design, groove effect must be taken into consideration, especially grooves with round or fillet corners.

Stiffness and Damping Coefficients
Following the similar procedure in Refs. [1,8], the nondimensional stiffness and damping coefficients for a circular bearing with multi-axial grooves were obtained by summarizing the coefficients of stiffness and damping over all supporting pads and are expressed as (1 + ε · cos α T 1 ) 3 · K P (η 1 ) · cos(π − ϑ 1 ), (1 + ε · cos α T 1 ) 3 · K P (η 1 ) · cos(π − ϑ 1 ), (1 + ε · cos α T 1 ) 3 · K P (η 1 ) · sin(π − ϑ 1 ),  43) and (44), it can be seen that they only depend on location angles and the number of grooves. This means they are changing with different groove configurations. The same is applied to the damping coefficients. For purpose to make comparison with other available methods, here a new group of non-dimensional coefficients of stiffness and damping is defined as follows: (43) (1 + ε · cos α T 1 ) 3 · K P (η 1 ) · sin(π − ϑ 1 ), (1 + ε · cos α Ti ) 3 · cos 2 (π − ϑ i ) (1 + ε · cos α T 1 ) 3 · cos 2 (π − ϑ 1 ), The Sommerfeld number in Eqs. (47) and (48) is the one defined by Eq. (42). The final dimensional coefficients of stiffness and damping are as follows: The non-dimensional stiffness group in Eq. (47) and non-dimensional damping group in Eq. (48) is directly comparable with existing circular bearing results such as long and short bearing theory and others. In this paper, the stiffness and damping coefficients from Childs and Moes [4,5] are of particular interest. They are considered to be accurate for non-grooved plain bearings and used for verifying the correctness of Eqs. (47) and (48). Same as for non-dimensional coefficients of stiffness, the coefficients of damping have also been compared at same groove number and L/D ratio. In Figures 13 and  14, the stiffness and damping coefficients by Childs and Moes was based on L/D = 2.0 and for grooved bearings was based on groove number = 8.
It is understandable that stiffness coefficient K yy for grooved bearing is slightly greater than that for nongrooved bearing. This is because the pressure is more concentrated on the area around of loading center due to grooves. The same reason may explain why the coefficient of stiffness K xx is lower than that of non-grooved bearing. The cross coefficient of stiffness K yx shows a different behavior from non-grooved bearing. Another noticeable characteristic is that the turning point of the cross-stiffness coefficient K xy into negative is shifted to lower eccentricity ratio.
The damping coefficient C yy is almost identical for both non-grooved bearing and grooved bearing in this particular geometrical condition. The coefficient C xx of grooved bearing is lower than that of non-grooved bearing. The cross-damping coefficient C xy = C yx has larger difference for low and high eccentricity ratio and small difference for intermediate eccentricity ratio.

Influence of the Number of Grooves
As stated in previous section, the coefficients of stiffness and damping are not only a function of eccentricity ratio, but also the number of grooves. Figure 15 is a comparison between two bearings with 8 grooves and 12 grooves, respectively. The effect on the coefficients of stiffness is different. An increased groove number has an insignificant effect on K yy while it reduces K yx . For crosscoefficient of stiffness K xy , the increase of groove number shifts the turning point to negative to lower eccentricity ratio. Figure 16 presents a comparison between the coefficients of damping for two bearings with 8 grooves and 12 grooves respectively. Again, the groove number has less effect on C yy while affecting other coefficients significantly.

Conclusions
This paper provides a new method to calculate the load capacity, the coefficients of stiffness and damping for water lubricated guide bearings with multi-axial grooves. The focus is on the effect of grooves and groove number. The paper doesn't include the effect of surface deformation. The result is an approximation and can be applied to water lubricated bearings made from hard polymers combined with lower pressure or other materials, such as Lignum Vitae wood and ceramics. The paper uses a socalled mixed scheme which means using parabolic slider for the first pad only, rest of the pads uses linear slider. The stiffness and damping of the grooved bearing was investigated considering groove effect. The coefficients of stiffness and damping demonstrated different characteristics from those with no grooves. Since the coefficients of stiffness and damping are function of eccentricity ratio and number of grooves, the effect of number of grooves was studied in great depth. It showed that the number of grooves has less effect on the coefficient K yy and C yy while it has a larger effect on other coefficients of stiffness and damping. Further research in considering surface deformation with using similar modeling could be an interesting subject. Small perturbation method means to find a solution of Eq. (A4) not far from the steady state solution with a linearization approach. This implies to find a solution, such as where p * 0 is non-dimensional pressure under steady operation, p * 1 is perturbation amplitude of a dynamic pressure on top of the pressure under steady operation. In true

Abbreviations
sense p * 1 is a coefficient of the non-dimensional dynamic pressure. δ is a small perturbation, which is a small number much less than 1.0. Its physical meaning is the ratio of amplitude change of film thickness to the minimum film thickness under steady operation. Insert Eqs. (A5) and (A6) into Eq. (A4), and equating the coefficients of zero order of "δ" on left and right side of Eq. (A4), it is resulted in an equation for pressure p * o By the same token, by equating the coefficients of first order of "δ" on left and right side of Eq. (A4), the coefficient of dynamic pressure p * 1 will fulfill following equation In this procedure, all other terms with orders equal to and higher than δ 2 are neglected.
The boundary conditions for the non-dimensional pressurep * are To fulfill these conditions, the non-dimensional pressure on steady operation p * 0 as well as the real and (A9) p * = 0 for x * = 0 and x * = −1.
imaginary part of non-dimensional dynamic pressure all need to be zero on the boundaries. This is expressed as First is to find the solution of Eq. (A7). By integrating twice of Eq. (A7), the non-dimensional pressure on steady operation is expressed in following form (A10) The boundary condition for p * 0 requires C 2 = 0, and Insert C 1 into Eq. (A11), the final non-dimensional pressure on steady operation takes form as below The load capacity function is the integration of the non-dimensional pressure (Eq. (A13)) The final result after implementation of the integration is  It is interesting to notice that there is a similarity of right side of Eq. (A7) and the first term on right side of Eq. (A8). Since the solution of Eq. (A7) creates the load capacity function Eq. (A15), the first real term on the right side of Eq. (A8) must generate the stiffness function. This concludes that the stiffness function is just equal to two times of the load capacity function by amount, therefore K P (η) = 2 (η − 2) · tan −1 √ η − 1 + √ η − 1 η 2 tan −1 √ η − 1 + (η + 2 3 ) √ η − 1 . Corresponding real part of non-dimensional dynamic pressure coefficient will be The next task is to find imaginary part of the nondimensional dynamic pressure coefficient p * 1 which needs to fulfill following equation: Following similar procedure to solve Eq. (A7), after integration twice of Eq. (A18), the imaginary part of nondimensional dynamic pressure coefficient is expressed with Utilizing boundary condition to above equation, the two constants are Inserting them into Eq. (A19) and integrating it over from x* = −1 to x* = 0, the damping function is as following: The total pad pressure appears as complex function which is The location of load center under steady operation and dynamic vibration is slightly different. The location of load center for steady operation is calculated with (A17) p * 1,r = −4 η 2 tan −1 √ η − 1 + (η + 2/3) (A20) C 3 = 8 4η 2 − 2(η + 1) 3η 2 √ η − 1 · tan −1 √ η − 1 + (η − 1)(3η + 2) , C 4 = 4 η − 1 .
(A21) C P (η) = − 0 −1 p * 1,i (x * , η)·dx * = 2(2η + 1) (A23) And the location of load center for dynamic load only is calculated with Since p * 1r is two times of static pressure p * o and has dominate amount in comparison to p * 1i , the value of Eq. (A24) is not very much different from the value from Eq.
Similarly this ratio is also defined for exponential and linear slider (see Figure 4d). This paper used static load center for Sommerfeld Number evaluation and dynamic load center for stiffness and damping evaluation for all three types of sliding bearings. The notion A Ed and A Ld presents the dynamic load center of exponential and linear slider respectively. x · p * 2 1r + p * 2 1i dx 0 −1 p * 2 1r + p * 2 1i dx .