Effects of Texture Bottom Profile on Static and Dynamic Characteristics of Journal Bearings

Li, Peng; Zeng, Fang; Xiao, Sen; Zhen, Dong; Zhang, Hao; Shi, Zhanqun

doi:https://doi.org/10.1155/2021/7068744

Shock and Vibration

On this page

Abstract Introduction Results and Discussion Conclusions Appendix Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Dynamic Analysis and Control Applied in Nonlinear Rotor Systems

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 7068744 | https://doi.org/10.1155/2021/7068744

Effects of Texture Bottom Profile on Static and Dynamic Characteristics of Journal Bearings

Peng Li,¹Fang Zeng,¹Sen Xiao,¹Dong Zhen,¹Hao Zhang,¹and Zhanqun Shi¹

Academic Editor: Hulun Guo

Received02 Sept 2021

Revised05 Oct 2021

Accepted22 Oct 2021

Published02 Dec 2021

Abstract

The purpose of this paper is to numerically study the effect of texture bottom profile on static, dynamic, and stability performance parameters of hydrodynamic journal bearings. The different performance parameters of square textured journal bearings with different bottom profiles are numerically investigated and compared with those of smooth journal bearing. There are five bottom profiles of this square texture: flat, curved, isosceles triangle (T1), oblique triangle (T2), and oblique triangle (T3). The static and dynamic coefficients are calculated by solving the steady-state Reynolds equation and the perturbation equations with FDM numerical technique. The performance characteristics under different texture distribution, depth, and bottom profiles are studied, and the current numerical results show that the selection of texture parameters is crucial to improve the static, dynamic, and stability performances of hydrodynamic journal bearing. Meanwhile, it is also found that the square texture with a flat bottom profile has a higher improvement in the values of static performance parameters in comparison to those other bottom profiles. Moreover, the simulation results indicate that the dynamic and stability performances improvement of textured journal bearing is also significant, especially when the eccentricity ratio is smaller.

1. Introduction

The surface texture method originated from surface roughness but is different from the surface roughness. Importantly, as a microbearing, the surface texture can obtain local hydrodynamic pressure and increase the load-carrying capacity due to its microwedge effect [1–3]. Furthermore, the introduction of textures of reasonable size and density can also reduce the friction and increase the service life of the mechanical seal [4–7]. Moreover, the method of introducing surface texture into friction pair has been proved to be an effective mean to improve the tribological performance characteristics.

In the past few decades, some theoretical and experimental studies on textured journal bearings have been published by numerous researchers. Ausas et al. [8] solved the Elrod and Adams [9] model with the mass conservation algorithm and found that the minimum film thickness and friction torque were increased slightly for fully rectangular textured journal bearings. Lu and Khonsari [10] carried out experimental investigation to analyze the performances of journal bearings with different dimples and found that fully textured journal bearing with reasonable texture parameters can minimize the friction coefficient under mixed lubrication condition. Another, Yamada et al. [11] carried out numerical and experimental investigation to study textured journal bearings with square dimples and reported that larger eccentricity ratio and smaller attitude angle were obtained compared with the smooth journal bearing. In addition, Tala-Ighil et al. [12, 13] numerically compared the performance of spherical and cylindrical textured journal bearings with different distributions and found that partial texture region has a positive effect on the performance of journal bearings while full texture region has a negative effect. Brizmer and Kligerman [14] demonstrated that partial texture distribution can obtain greater load-carrying capacity than full texture distribution for spherical textured journal bearing. Also, Zhang et al. [15] optimized the spherical texture parameters to achieve greater load-carrying capacity of journal bearings. Kango et al. [16] found that microgrooved journal bearing can reduce the friction coefficient and average temperature better than spherical textured journal bearing, and the performance gain is more obvious at lower eccentricity ratio. Researchers [17–19] used a multiobjective optimization method to optimize the texture parameters of the microgrooved journal bearing to obtain greater load-carrying capacity and smaller friction coefficient. Lately, Filgueira Filho et al. [20] experimentally analyzed the performances of textured journal bearings with chevron, saw tooth, oblong dimple, and aligned dimple and found that the chevron textured journal bearing exhibited smaller eccentricity under reasonable texture geometrical parameters. Besides, some researchers [21–23] also analyzed the influence of protruded and dimple textures on the load-carrying capacity and friction of journal bearings and compared the results with that of smooth journal bearing.

Besides the above research studies on the static performance characteristics of textured journal bearings, some researchers have found that surface texture also affects dynamic and stability performance parameters in recent years. Importantly, stiffness and damping coefficients have significant effects on the dynamic performance of rotating systems as important criteria for designers to evaluate performance gains [24]. Jiang et al. [25] investigated the influence of spherical convex and concave textures on dynamic performance of journal bearings and reported that stability critical speed of journal bearing could be obtained with appropriate texture geometric parameters. Further, Yamada et al. [26, 27] numerically and experimentally analyzed the influence of full square texture on the dynamic performance and linear stability-threshold shaft speed of journal bearing; the results show that lower cross-coupled stiffness coefficients and higher stability-threshold speed of journal bearing were obtained. Besides, Matele and Pandey [28] studied the dynamic characteristics of journal bearings with circular, square, and densely distributed square textures and reported that the dynamic performances of journal bearing with properly texture geometric parameters have been significantly improved. Similarly, the effects of triangular-shaped and chevron-shaped textures on the static, dynamic, and stability performances of journal bearings have been studied, and the optimal performance parameters have been found by researchers [29–32]. In addition, Meng et al. [33] carried out a calculative investigation to study the effects of simple rectangular texture and compound texture on the dynamic performance of journal bearings and found that locally distributed compound texture can increase the direct stiffness coefficients, damping coefficients, and critical speed of journal bearing. Recently, Singh and Awasthi [34, 35] numerically investigated the effects of texture parameters (texture depth and its distribution) on the static, dynamic, and stability performances of journal bearings with cylindrical, spherical, triangular, and kite textures; the results show that the performance parameters were significantly improved when the texture distribution is in the pressure rising region of bearings.

In numerous above works, the microtextures of several geometric shapes introduced on the bearing surface have been studied such as triangular, chevron, and spherical. Further, it is worth noting that the texture internal structure also plays an important role for tribological performances. Nanbu et al. [36] numerically investigated the effects of texture bottom profiles on lubrication performance and found that the microwedge effect of the texture bottom profile was an important parameter. In addition, Shen and Khonsari [37] carried out numerical and experimental investigation to demonstrate that dimple with flat bottom profile produces higher load-carrying capacity than those triangular profiles. Similarly, Rahmani et al. [38, 39] numerically studied the tribological performances of flat bottom textured journal bearing and found the maximum load-carrying capacity and minimum friction coefficient than those triangular profiles. Qiu et al. [40] conducted parametric optimization analysis on the textures of six different texture shapes and found that the ellipsoidal texture with curved bottom profiles could obtain the minimum friction coefficient and the maximum stiffness by comparison. Some researchers [41–43] used optimization algorithms to optimize the bottom profile of the texture, which improved the tribological performance. However, it should be noted that the selection principle of the texture optimal design parameters is not universal because the improvement of tribological performance is closely related to the operating conditions [44].

Although researchers have done numerous studies on the influence of texture parameters on performance characteristics of hydrodynamic journal bearing, under specific texture parameters, texture bottom profile affects static performance characteristics of journal bearing, such as load-carrying capacity, and friction force. And each bottom profile has a difference in static performance gain of journal bearing. However, to the knowledge of the authors, a comparative study of the effects of texture bottom profiles (flat, curved, isosceles triangular T1, oblique triangular T2, and oblique triangular T3) on the dynamic and stability performance characteristics of journal bearings by changing texture parameters (texture distribution and depth) and eccentricity ratios has not been reported. Therefore, the purpose of this paper is to systematically study the static, dynamic, and stability performance characteristics of square textured journal bearing with different bottom profiles. Thus, the texture distribution and depth are considered, and then the performance parameters of textured journal bearings are compared with those of smooth journal bearing to obtain the optimal texture parameters under specific conditions.

2. Theoretical Analysis

2.1. Reynolds Equation

In the study, the schematic diagram of textured journal bearing with the coordinate system and the expanded surface of textured journal bearing is shown in Figures 1(a) and 1(b). Besides, the square texture with different bottom profiles such as flat, curved, isosceles triangular T1, oblique triangular T2, and oblique triangular T3 is introduced on the bearing surface in Figure 1(c).

(a)

(b)

(c)

In the fluid lubrication state, the Reynolds equation based on the assumptions of isothermal, laminar, and incompressible Newtonian fluids in Cartesian coordinate can be written aswhere represents the fluid film pressure, represents the fluid film dynamic viscosity, denotes the journal speed in the direction, and denotes the time.

Besides, in order to unify the full fluid film region and cavitated region, the film rupture and the reformation boundary can be automatically determined in the numerical solution process, and the dimensionless dependent variable and the switch function are introduced as [45].where represents the reference pressure and represents the cavitation pressure.

Moreover, the fluid density can be showed as follows [46]:where represents the local fluid density at the cavitation region.

The dimensionless parameters used can be expressed aswhere denotes the radius of bearing, represents length of bearing, denotes radial clearance, denotes angular velocity by the rotation of the journal, and and represent the circumferential and axial coordinates of journal bearing.

Combined with equations (2) to (5), the dimensionless form of the Reynolds equation with mass conservation cavitation effect can be written as follows [46]:

And then the dimensionless fluid film pressure can be calculated by

2.2. Fluid Film Thickness Equations

The dimensionless fluid film thickness for textured journal bearing can be described aswhere represents the eccentricity ratio of journal bearing, represents the circumferential axial coordinates of journal bearing, and represents the dimensionless fluid film thickness contributed by surface texture. Besides, the expression for a square texture with different bottom profiles such as flat, curved, isosceles triangular T1, oblique triangular T2, and oblique triangular T3 can be calculated as follows [36].

For smooth surface,

For flat bottom profile,

For curved bottom profile,

For isosceles triangular T1 bottom profile,

For oblique triangular T2 bottom profile,

For oblique triangular T3 bottom profile,

2.3. Fluid Film Perturbation Pressure Equations

For the calculation of the dimensionless fluid film perturbation pressure, the dimensionless fluid film pressure is obtained from the steady-state Reynolds equation (10a) without considering the time term and equation (7), and then combined with the fluid film thickness equations (8) and ((9a) to (9f)), the dimensionless form of steady-state Reynolds equation and fluid film perturbation pressure equations in coordinate systems can be obtained as follows:(1)The dimensionless form of the steady-state Reynolds equation is as follows:(2)The dimensionless form of fluid film perturbation pressure equation of is as follows:(3)The dimensionless form of fluid film perturbation pressure equation of is as follows:(4)The dimensionless form of fluid film perturbation pressure equation of is as follows:(5)The dimensionless form of fluid film perturbation pressure equation of is as follows:

2.4. Boundary Conditions

The boundary conditions of the dimensionless steady-state Reynolds equation (10a) without considering the time term for journal bearing are given as follows:where represents the environment pressure.

And the periodic boundary conditions are given as

The boundary conditions of the dimensionless fluid film perturbation pressure equations (10b)–(10e) for journal bearing are given as

Here, represents the dimensionless fluid film perturbation pressure.

2.5. Performance Characteristics of Journal Bearing

The expressions of the static, dynamic, and stability performance parameters for journal bearing are presented in this section. Besides, it should be noted that the coordinate system is used in the calculation of the dimensionless performance parameters; once the equilibrium position of journal bearing is determined, the coordinate system needs to be reduced by an angle of and then converted to the Cartesian coordinate system.

2.5.1. Static Performance Characteristics

In order to evaluate the performance characteristics of journal bearings for steady-state condition, when determining the journal center equilibrium position for a given eccentricity ratio, the dimensionless static performance characteristics can be calculated. Among them, the dimensionless static performance parameters include the load-carrying capacity, attitude angle, friction force, friction coefficient, and end leakage flow rate.

(1) The Load-Carrying Capacity. The load-carrying capacity can be determined by the method of force balance. The components of load-carrying capacity, in the circumferential and axial directions, can be obtained by integrating the hydrodynamic pressure generated in the clearance gap of journal bearing, whose dimensionless expressions are given as follows [47]:

The load-carrying capacity is defined as the total force that the journal bearing system can support, whose dimensionless expressions are expressed as follows:

The components of dimensionless load-carrying capacity in the Cartesian coordinate can be obtained as follows:

(2) The Attitude Angle. The attitude angle of the journal bearing is expressed as follows:

(3) The Friction Force and the Friction Coefficient. The dimensionless form of viscosity shear stress due to lubricant shearing behavior can be obtained as follows:

The frictional force of lubricant on the journal surface can be calculated by integrating the shear stress, whose dimensionless expression is given as follows:

The dimensionless friction coefficient of lubricant on the journal surface can be calculated as follows:

(4) The End Leakage Flow Rate. The dimensionless end leakage flow rate has two components such as from the left-hand bearing end and from the right-hand bearing end which are given as follows:

The total dimensionless end leakage flow rate can be calculated as follows:

2.5.2. Dynamic Performance Characteristics

When journal center is disturbed from the static equilibrium position, the dimensionless dynamic performance characteristics can be calculated. Among them, the dimensionless dynamic performance parameters include the fluid film stiffness coefficients and fluid film damping coefficients. Importantly, the stiffness and damping coefficients are the spring and dashpot systems, which constrain the movement of the journal during operation [24].

(1) The Fluid Film Stiffness Coefficients. After calculating the fluid film perturbation pressure, the fluid film stiffness coefficients can be calculated by integrating the fluid film perturbation pressure, whose dimensionless expressions in the coordinate systems are expressed as follows [47]:

The dimensionless fluid film stiffness coefficients in the Cartesian coordinate system are expressed as follows:

The direct stiffness coefficients and reflect the ability to resist external loads in the and directions, respectively. The larger the values of the direct stiffness coefficients are, the better the stability of the journal bearing system is. The cross-coupled stiffness coefficients and reflect the unstable energy input of the bearing. The smaller the absolute values of the cross-coupled stiffness coefficients are, the better the stability of the journal bearing system is. Therefore, the stiffness coefficients have significant influence on the overall dynamic characteristics of the journal bearing system [24].

(2) The Fluid Film Damping Coefficients. After calculating the fluid film perturbation pressure, the fluid film damping coefficients can be calculated by integrating the fluid film perturbation pressure, whose dimensionless expressions in the coordinate systems are expressed as follows [47]:

The dimensionless fluid film damping coefficients in the Cartesian coordinate system can be obtained as follows:

The damping coefficients reflecting the speed of the journal are restored to stability after being impacted by external loads. The larger the values of the damping coefficients are, the better the stability of the journal bearing system is. Similarly, the damping coefficients have significant influence on the overall dynamic characteristics of the journal bearing system [24].

2.5.3. Stability Performance Characteristics

(1) The Critical Mass. For the journal bearing system, the dimensionless form of linearized motion equation is expressed with the known eight dynamic coefficients as follows [47, 48]:where and are the displacement, , , and are velocity, and and are the accelerations, for the journal in and the directions. Importantly, it should be noted here that the system is considered stable when the journal mass is smaller than the critical mass . The dimensionless critical mass is written as follows [47, 48]:

Here,

(2) The Threshold Speed of Stability. The threshold speed of stability is an important parameter to evaluate the stability of the journal bearing system. Importantly, the value of threshold speed of stability is higher than that of the journal speed, indicating that the journal bearing system will become more stable. After calculating the dimensionless value of the critical mass, the dimensionless value of threshold speed of stability is written as follows [47, 48]:

Here, is the reaction force of fluid film .

3. Solution Procedure and Numerical Methods

In the present study, in order to calculate the influence of texture parameters on the static, dynamic, and stability performance characteristic parameters of journal bearing, the global numerical calculation process is divided into two parts: the steady-state parameters numeration (step 2 to step 8) and dynamic parameters numeration (step 9 to step 11). Moreover, the computation process is operated in MATLAB R2014a. Figure 2 illustrates the flowchart of the global solution procedure, and the details of numerical procedure and methods are shown in the Appendix.

4. Results and Discussion

In this section, the static, dynamic, and stability performance characteristic parameters of journal bearings are numerically calculated. Besides, five different bottom profiles of a square texture, flat, curved, isosceles triangle (T1), oblique triangle (T2), and oblique triangle (T3) have been introduced on the bearing surface, and the values of circumferential texture region, texture depth, and eccentricity ratio are considered. Furthermore, performance characteristic parameters of textured journal bearings are numerically compared with those of smooth journal bearing. Here, the operating parameters and texture parameters in the current study are listed in Table 1. It should be emphasized that, unless otherwise specified, the values of these parameters remain unchanged in the following theoretical analysis.

4.1. Validation

In order to verify the accuracy of current model numerical solution, model validation has been shown in the two types. In the first type, for the verification study, the static performance characteristics of smooth, partially textured, and fully textured journal bearings are calculated. Table 2 presents the results of comparison of the current model compared with Tala-Ighil et al. [12]. Besides, in order to clearly show the validation results, the dimensionless fluid film thicknesses and fluid film pressures of smooth, partially textured, and fully textured journal bearings are shown in Figure 3. In the second type, for the verification study the static and dynamic performance characteristics of partially textured journal bearing with texture region located at rising pressure region, Figure 4 presents the results of comparison of the current model compared with Singh and Awasthi [34], It is important to emphasize that because of the external force acting in the direction, the direct stiffness coefficient is the key factors to consider for dynamic performance; the direct stiffness coefficient and the load-carrying capacity can be improved in rising pressure texture region. In particular, the performance gains are superior at lower eccentricity ratios. Moreover, the geometric parameter values of journal bearing in the validation are the same as those in the reference studies, and the numerical differences obtained are all within the acceptable range.

(a)

(b)

(c)

(d)

(e)

(f)

(a)

(b)

(c)

(d)

(e)

(f)

4.2. Influence of the Circumferential Texture Region

In this subsection, the influences of the circumferential texture region on the static, dynamic, and stability performance parameters of smooth/square textured journal bearings with flat, curved, triangle T1, triangle T2, and triangle T3 bottom profiles are observed in Figures 5 and 6. Here, five circumferential texture regions are taken into consideration, rising pressure texture region (0°–128.5°), first half texture region (0°–180°), middle texture region (128.5°–257°), second half texture region (180°–360°), and full texture region (0°–360°), and the corresponding number of textures along the circumferential direction is 5, 7, 5, 7, and 14, respectively.

(a)

(b)

(c)

(d)

(e)

(f)

(a)

(b)

(c)

(d)

(e)

(f)

As seen from Figures 5(a) to 5(f), the influence of the circumferential texture region on the static performance parameters of textured journal bearings is shown, such as the load-carrying capacity (), the maximum fluid film pressure (), the friction coefficient (), the end leakage flow rate (), the attitude angle (), and the friction force (). From Figures 5(a) and 5(b), it is observed that the load-carrying capacity () and the maximum film pressure () of textured journal bearings are higher than those of smooth journal bearing in rising pressure texture region and first half texture region, while other texture distributions have negative effects. Furthermore, it is found that square textured journal bearing with flat bottom profile has the highest values of and amongst other textured journal bearings with different bottom profiles in rising pressure texture region. In Figure 5(c), the variation values of the friction coefficient () are opposite to those of the load-carrying capacity () of textured journal bearings in each texture region. Moreover, it is found that square textured journal bearing with flat bottom profile has the lowest value of amongst other textured journal bearings with different bottom profiles in rising pressure texture region. In Figure 5(d), the end leakage flow rate () of textured journal bearings is higher than that of smooth journal bearing in rising pressure texture region, first half texture region, and full texture region, and then values are lower than those of smooth bearing in other texture distributions. Furthermore, it is found that square textured journal bearing with flat bottom profile has the highest value of the end leakage flow rate () amongst other textured journal bearings with different bottom profiles in rising pressure texture region. In Figure 5(e), the attitude angle () of textured journal bearings is lower than that of smooth bearing in rising pressure texture region, first half texture region, and full texture region, and then values are higher than those of smooth bearing in other texture regions. Moreover, it is found that square textured journal bearing with flat bottom profile has the lowest value of amongst other textured journal bearings with different bottom profiles in first half texture region. In Figure 5(f), the friction force () of textured journal bearings is lower than that of smooth bearing in each texture region. Besides, it is found that square textured journal bearing with flat bottom profile has the lowest value of amongst other textured journal bearings with different bottom profiles in full texture region.

As seen from Figures 6(a) to 6(f), the influence of the circumferential texture region on the dynamic and stability performance parameters of textured journal bearings is shown, such as the direct stiffness coefficients , the direct damping coefficients , the critical mass , and the threshold speed of stability (). In Figure 6(a), the direct stiffness coefficients () of textured journal bearings are lower than those of smooth bearing in each texture regions except in first half texture region. In addition, it is found that square textured journal bearing with flat bottom profile has the highest value of amongst other textured journal bearings with different bottom profiles in first half texture region. In Figure 6(b), the direct stiffness coefficients () of textured journal bearings are higher than those of smooth bearings in rising pressure texture region and second half texture region; nevertheless, the direct stiffness coefficients () of textured journal bearings are lower than those of smooth bearings in other texture regions. Moreover, it is found that square textured journal bearing with flat bottom profile has the highest value of amongst other textured journal bearings with different bottom profiles in second half texture region. From Figures 6(c) and 6(d), it is observed that the direct damping coefficients of textured journal bearings are lower than those of smooth bearing in each texture region. Furthermore, it is found that smooth journal bearing has the highest values of and amongst textured journal bearings with different bottom profiles in each texture region. In Figure 6(e), the critical mass () of textured journal bearings is higher than that of smooth bearings in rising pressure texture region and first half texture region; nevertheless, the critical mass () of textured journal bearings is lower than that of smooth bearings in other texture regions. Besides, it is found that square textured journal bearing with flat bottom profile has the highest values of amongst other textured journal bearings with different bottom profiles in first half texture region. In Figure 6(f), it is observed that the threshold speed of stability () of textured journal bearings is lower than that of smooth bearing in rising pressure texture region, and then the threshold speed of stability () of textured journal bearings is higher than that of smooth bearing in other texture regions. Furthermore, it is found that square textured journal bearing with flat bottom profile has the lowest values of amongst other textured journal bearings with different bottom profiles in rising pressure texture region.

Based on the above analysis of performance parameters of journal bearing, the results indicate that different texture regions have different effects on the performances of textured journal bearings with different bottom profiles. It is found that the textured journal bearing can maximize the load-carrying capacity (), the maximum fluid film pressure (), and lubrication flow ratio () and minimize the friction coefficient (), compared with smooth journal bearing. It is important to emphasize that because of the external force acting in the direction, the direct stiffness coefficient () and the critical mass () are the key factors to consider for dynamic performance. The direct stiffness coefficient () and the critical mass () can also be improved simultaneously in rising pressure texture region. Moreover, the performance of textured journal bearing with flat bottom profile is better than that of textured journal bearings with other bottom profiles. Therefore, the performance of square textured journal bearings in rising pressure texture region will be analyzed in the following sections.

4.3. Influence of the Texture Depth

In this subsection, it should be emphasized that the selected texture region is the rising pressure texture region of journal bearing. The influences of the texture depth () on the static, dynamic, and stability performance parameters of smooth/square textured journal bearings with flat, curved, triangle T1, triangle T2, and triangle T3 bottom profiles are shown and discussed in Figures 7 and 8, and the texture depth () ranges from 0 to 1.2 in the study.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 7

Effect of the texture depth on static performance characteristics for smooth/textured journal bearings: (a) the load-carrying capacity (); (b) the maximum fluid film pressure (); (c) the friction coefficient (); (d) the end leakage flow rate (); (e) the attitude angle (); (f) the friction force (). Bearings are lower than that of smooth bearing in the texture depth () ranges. However, the end leakage flow rate () keeps increasing and the friction force () keeps decreasing, respectively.

(a)

(b)

(c)

(d)

(e)

(f)

As seen from Figures 7(a) to 7(f), the influence of the texture depth () on the static performance parameters of textured journal bearings is shown, such as the load-carrying capacity (), the maximum film pressure (), the friction coefficient (), the end leakage flow rate (), the attitude angle (), and the friction force (). From Figures 7(a) and 7(b), the results show that the load-carrying capacity () and the maximum film pressure () of textured journal bearings are higher than those of smooth bearing in the texture depth () ranges. Moreover, with the increase in the texture depth (), the load-carrying capacity () and the maximum film pressure () of textured journal bearings firstly increase and then decrease, respectively. It can also be found that there is an optimum value of texture depth for the load-carrying capacity () and the maximum film pressure (), and that depends on the operating conditions of the bearing. Furthermore, the optimal texture depth of textured journal bearings with different bottom profiles is different, and square textured journal bearing with flat bottom profile has the highest values of and amongst other textured journal bearings with different bottom profiles. From Figures 7(c) and 7(d), the results show that the friction coefficient () and the attitude angle () of textured journal bearings are lower than those of smooth bearing in the texture depth () ranges. Contrary to the changing trend of the load-carrying capacity (), with the increase in the texture depth (), the friction coefficient () and the attitude angle () of textured journal bearings firstly decrease and then increase, respectively. Furthermore, it is also found that square textured journal bearing with flat bottom profile has the lowest values of and amongst other textured journal bearings with different bottom profiles. From Figures 7(d) and 7(f), the results show that the end leakage flow rate () of textured journal bearings is higher than that of smooth bearing and the friction force () of textured journal.

As seen from Figures 8(a) to 8(f), the influence of the texture depth () on the dynamic and stability performance parameters of textured journal bearings is shown, such as the direct stiffness coefficients (), the direct damping coefficients (, ), the critical mass (), and the threshold speed of stability (). And from the results, the direct stiffness coefficients (), the direct damping coefficients (, ), and the threshold speed of stability () of textured journal bearings are lower than those of smooth bearing in the texture depth () ranges; moreover, with the increase in the texture depth (), the direct stiffness coefficients (), the direct damping coefficients (, ), and the threshold speed of stability () of textured journal bearings are always decreased in Figures 8(a), 8(c), 8(d), and 8(f), respectively. Furthermore, it is found that square textured journal bearing with flat bottom profile has the lowest values of , , , and amongst other textured journal bearings with different bottom profiles. From Figures 8(b) and 8(e), the direct stiffness coefficient () and the critical mass () firstly increase and then decrease, respectively. It can also be found that there is an optimum value of texture depth for the direct stiffness coefficient () and the critical mass (), and that depends on the operating conditions of the bearing. Furthermore, the optimal texture depth of textured journal bearings with different bottom profiles is different, and square textured journal bearing with flat bottom profile has the highest values of and amongst other textured journal bearings with different bottom profiles.

Based on the above analysis of performance parameters of journal bearing, the results indicate that different texture depth () has different effects on the performances of textured journal bearings with different bottom profiles, and the optimal texture depth of each bottom profiles is also different. There is an optimum value of texture depth that maximize the load-carrying capacity (),maximum fluid film pressure (), lubrication flow ratio (), direct stiffness coefficient (), and critical mass () and minimize the friction coefficient (), respectively, compared with smooth journal bearing. Moreover, for the texture depth () within the selected range, the performance of textured journal bearing with flat bottom profile is better than that of square textured journal bearings with other bottom profiles.

4.4. Influence of the Eccentricity Ratio

In this subsection, it should be emphasized that the selected texture region is the rising pressure texture region of journal bearing. The influences of the eccentricity ratio () on the static, dynamic, and stability performance parameters of smooth/square textured journal bearings with flat, curved, triangle T1, triangle T2, and triangle T3 bottom profiles are presented and discussed in Figures 9 and 10, and the eccentricity ratio () ranges from 0.2 to 0.8 in the study.

(a)

(b)

(c)

(d)

(e)

(f)

(a)

(b)

(c)

(d)

(e)

(f)

As seen from Figures 9(a) to 9(f), the influence of the eccentricity ratio () on the relative difference on static performance parameters of textured journal bearings relative to smooth bearing is presented, such as the load-carrying capacity (), the maximum film pressure (), the friction coefficient (), the end leakage flow rate (), the attitude angle (), and the friction force (). And the results show that, with the increase in the eccentricity ratio (), the load-carrying capacity (), the maximum film pressure (), the end leakage flow rate (), and the friction force () of smooth journal bearing increase, while the friction coefficient () and the attitude angle () of smooth journal bearing decrease accordingly. From Figures 9(a) to 9(c), the load-carrying capacity () and the maximum film pressure () of textured journal bearings are higher than those of smooth bearing, and the friction coefficient () of textured journal bearings is lower than that of smooth bearing for the eccentricity ratio () ranges from 0.2 to 0.6; however, the trend of change of these values is the opposite for the eccentricity ratio () exceeding 0.6. Furthermore, for the value of the eccentricity ratio 0.3 (), it is also found that square textured journal bearing with flat bottom profile has the highest values of , , and the lowest value of amongst other textured journal bearings with different bottom profiles. In Figure 9(d), regardless of the change of eccentricity ratio (), the end leakage flow rate () of textured journal bearings is always higher than that of smooth bearing. Furthermore, for the value of the eccentricity ratio 0.3 (), it is also found that square textured journal bearing with flat bottom profile has the highest values of amongst other textured journal bearings with different bottom profiles. From Figures 9(e) to 9(f), the attitude angle () and the friction force () of textured journal bearings are always lower than those of smooth bearing. Furthermore, for the value of the eccentricity ratio 0.3 (), it is also found that square textured journal bearing with flat bottom profile has the lowest values of and amongst other textured journal bearings with different bottom profiles.

From Figures 10(a) to 10(f), the influence of the eccentricity ratio () on the relative difference in dynamic and stability performance parameters of textured journal bearings relative to smooth bearing is shown, such as the direct stiffness coefficients (), the direct damping coefficients (, ), the critical mass (), and the threshold speed of stability (). And from the results, with the increase in the eccentricity ratio (), the direct stiffness coefficients (), the direct damping coefficients (, ), the critical mass (), and the threshold speed of stability () of smooth journal bearing also increase accordingly. In Figures 10(a), 10(c), 10(d), and 10(f), regardless of the change of eccentricity ratio (), the direct stiffness coefficients (), the direct damping coefficients (, , ), and the threshold speed of stability () of textured journal bearings are always lower than those of smooth bearing. Furthermore, for the value of the eccentricity ratio 0.3 (), it is also found that square textured journal bearing with flat bottom profile has the lowest values of , , , and amongst other textured journal bearings with different bottom profiles. From Figures 10(b) and 10(e), the critical mass () and the direct stiffness coefficient () of textured journal bearings are higher than those of smooth bearing for the eccentricity ratio () ranges from 0.2 to 0.45, and then it is the opposite for the eccentricity ratio greater than 0.45. Furthermore, for the value of the eccentricity ratio 0.3 (), it is also found that square textured journal bearing with flat bottom profile has the highest values of and amongst other textured journal bearings with different bottom profiles. It is important to emphasize that because of the external force acting in the direction, the direct stiffness coefficient () is the key factor to consider for dynamic performance.

Based on the above analysis of performance parameters of journal bearing, the results indicate that different eccentricity ratios () have different effects on the performances of textured journal bearings with different bottom profiles. It can be found that the larger eccentricity ratio (), minimum film thickness will be smaller, and the hydrodynamic pressure effect will be more obvious. And the eccentricity ratio () is smaller, the relative difference of performance parameters of journal bearings provided by surface texture is more obvious. However, with the increase in eccentricity ratio (), the performances improvement become gradually weakened. Moreover, for the eccentricity ratio () within the selected range, the results show that square textured journal bearing with flat bottom profile can obtain better performance improvement.

5. Conclusions

In this paper, the effects of a square texture with flat, curved, isosceles triangle (T1), oblique triangle (T2), and oblique triangle (T3) bottom profiles on the static, dynamic, and stability performance parameters of journal bearings are numerically studied at different texture distribution and depth. Firstly, the film pressure is obtained through solving the steady-state Reynolds equation, and then the perturbation pressures are also obtained through solving the perturbation equations. Next, the static, dynamic, and stability performance parameters of journal bearings can be calculated. Finally, the performance characteristic parameters of the textured journal bearings with difference bottom profiles are compared with those of the smooth journal bearing. The relevant conclusions are drawn as follows:(1)Reasonable geometrical parameters of surface texture can not only improve the static performance but also improve the dynamic performances and further enhance the stability of the journal bearing system, compared with the smooth journal bearing.(2)For the circumferential texture region, it is found that textured journal bearing can maximize the load-carrying capacity (), maximum fluid film pressure (), and lubrication flow ratio () and minimize the friction coefficient (), compared with smooth journal bearing in rising pressure texture region. It is important to emphasize that because of the external force acting in the direction, the direct stiffness coefficient () and the critical mass () are the key factors to consider for dynamic performance; the direct stiffness coefficient () and the critical mass () can also be improved simultaneously in rising pressure texture region.(3)For the texture depth (), there is an optimum value of texture depth that maximize the load-carrying capacity (), maximum fluid film pressure (), lubrication flow ratio (), direct stiffness coefficient (), and the critical mass () and minimize the friction coefficient (), respectively, compared with smooth journal bearing. Moreover, the optimal texture depth of each bottom profiles is also different.(4)For the texture bottom profile, the results indicate that different bottom profiles have different effects on the performances of textured journal bearings, such as flat, curved, isosceles triangle (T1), oblique triangle (T2), and oblique triangle (T3) bottom profiles. Moreover, for the five different bottom profiles, the performance of textured journal bearing with flat bottom profile is better than that of textured journal bearings with other bottom profiles.(5)For the eccentricity ratio, the larger eccentricity ratio (), minimum film thickness will be smaller, and the hydrodynamic pressure effect will be more obvious. And the eccentricity ratio () is smaller, the relative difference of performance parameters of textured journal bearings provided by surface texture is more obvious, compared with smooth journal bearing. However, with the increase in eccentricity ratio (), the performances improvement become gradually weakened.

Based on the study of the performance characteristics of textured journal bearings, while designing, designers should focus on those optimal values of texture region, depth, and bottom profile given, which can maximize the static and dynamic performance and enhance the stability of journal bearings.

Appendix

The details of numerical procedure and methods are shown as follows:(1)Input data values of the eccentric ratio , the length to diameter ratio , and the texture and mesh parameters.(2)Guess the initial values of the dimensionless fluid film thickness , the dependent variable , and the switch function .(3)Calculate the dimensionless fluid film thickness for the smooth and textures different bottom profiles by equations (8) and ((9a)–(9f)).(4)Here, we set the circumferential and axial mesh points to 314 and 100, respectively, to ensure that each texture occupies 15 mesh points in length and width. The finite difference method is used to discretize the steady-state Reynolds equation (10a), and the difference forms are as follows:(5)Further, the dependent variable and the switch function are obtained by solving the following difference equations with the over-relaxation iterative method. Here, , , , , and represent the coefficients of the discretized steady-state Reynolds equation (10a), and the relaxation factor is set to 1.9.(6)The expression of the convergence condition of the dependent variable and the switch function is judged by the following equation. If not, update and and then return step 4. Here, the convergence precisions for and are set to .(7)Calculate the components of dimensionless load-carrying capacity and by equations (13a) and (13b). Further, the initial values (, ) are given, and calculate the attitude angle through the twofold secant method of equation as follows:(8)The expression of the convergence condition of the attitude angle is judged by the following equation. If not, update and then return step 7. Here, the convergence precision is set to .(9)Similarly, the dimensionless fluid film perturbation equations (10b) to (10e) are also separated into difference forms, and the difference forms are as follows:(10)Next, the fluid film perturbation pressures are also obtained by solving the following difference equations with the over-relaxation iterative method. Here, the relaxation factor is set to 1.9.(11)The expression of the convergence condition of the fluid film perturbation pressures is judged by the following equation. If not, update and then repeat step 10. Here, the convergence precision is set to .(12)Calculate the static, dynamic, and stability performance parameters for journal bearing by equations ((13a)–(13c)) to (25).

Data Availability

The figure and table data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest or personal relationships that would affect the work reported in this article.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (Grant nos. 52175085 and 51875166).

References

D. B. Hamilton, J. A. Walowit, and C. M. Allen, “A theory of lubrication by microirregularities,” Journal of Basic Engineering, vol. 88, no. 1, pp. 177–185, 1966.
View at: Publisher Site | Google Scholar
J. N. Anno, J. A. Walowit, and C. M. Allen, “Microasperity lubrication,” Journal of Lubrication Technology, vol. 90, no. 2, pp. 351–355, 1968.
View at: Publisher Site | Google Scholar
J. N. Anno, J. A. Walowit, and C. M. Allen, “Load support and leakage from microasperity-lubricated face seals,” Journal of Lubrication Technology, vol. 91, no. 4, pp. 726–731, 1969.
View at: Publisher Site | Google Scholar
I. Etsion and L. Burstein, “A model for mechanical seals with regular microsurface structure,” Tribology Transactions, vol. 39, no. 3, pp. 677–683, 1996.
View at: Publisher Site | Google Scholar
I. Etsion, Y. Kligerman, and G. Halperin, “Analytical and experimental investigation of laser-textured mechanical seal faces,” Tribology Transactions, vol. 42, no. 3, pp. 511–516, 1999.
View at: Publisher Site | Google Scholar
I. Etsion, “Improving tribological performance of mechanical components by laser surface texturing,” Tribology Letters, vol. 17, no. 4, pp. 733–737, 2004.
View at: Publisher Site | Google Scholar
I. Etsion, “State of the art in laser surface texturing,” Journal of Tribology, vol. 127, no. 1, pp. 248–253, 2005.
View at: Publisher Site | Google Scholar
R. Ausas, P. Ragot, J. Leiva, M. Jai, G. Bayada, and G. C. Buscaglia, “The impact of the cavitation model in the analysis of microtextured lubricated journal bearings,” Journal of Tribology, vol. 129, no. 4, pp. 868–875, 2007.
View at: Publisher Site | Google Scholar
H. G. Elrod and M. L. Adams, “A computer program for cavitation and starvation problems,” Cavitation and Related Phenomena in Lubrication, vol. 103, pp. 37–41, 1974.
View at: Google Scholar
X. B. Lu and M. M. Khonsari, “An experimental investigation of dimple effect on the stribeck curve of journal bearings,” Tribology Letters, vol. 27, no. 2, pp. 169–176, 2007.
View at: Publisher Site | Google Scholar
H. Yamada, H. Taura, and S. Kaneko, “Static characteristics of journal bearings with square dimples,” Journal of Tribology, vol. 139, no. 5, Article ID 051703, 2017.
View at: Publisher Site | Google Scholar
N. Tala-Ighil, M. Fillon, and P. Maspeyrot, “Effect of textured area on the performance of hydrodynamic journal bearing,” Tribology International, vol. 44, no. 3, pp. 211–219, 2011.
View at: Publisher Site | Google Scholar
N. Tala-Ighil, P. Maspeyrot, M. Fillon, and A. Bounif, “Effects of surface texture on journal-bearing characteristics under steady-state operating conditions,” Proceedings of the Institution of Mechanical Engineers - Part J: Journal of Engineering Tribology, vol. 221, no. 6, pp. 623–633, 2007.
View at: Publisher Site | Google Scholar
V. Brizmer and Y. Kligerman, “A laser surface textured journal bearing,” Journal of Tribology, vol. 134, no. 3, Article ID 031702, 2012.
View at: Publisher Site | Google Scholar
H. Zhang, G. N. Dong, M. Hua, F. F. Guo, and K. S. Chin, “Parametric design of surface textures on journal bearing,” Industrial Lubrication & Tribology, vol. 67, no. 4, pp. 359–369, 2015.
View at: Publisher Site | Google Scholar
S. Kango, R. K. Sharma, and R. K. Pandey, “Comparative analysis of textured and grooved hydrodynamic journal bearing,” Proceedings of the Institution of Mechanical Engineers - Part J: Journal of Engineering Tribology, vol. 228, no. 1, pp. 82–95, 2014.
View at: Publisher Site | Google Scholar
C. P. Liu, W. Y. Li, X. Q. Lu, and B. Zhao, “Effects of double parabolic profiles with groove textures on the hydrodynamic lubrication performance of journal bearing under steady operating conditions,” Mechanics and Industry, vol. 21, no. 3, pp. 1–14, 2020.
View at: Publisher Site | Google Scholar
C. P. Liu, N. Zhong, X. Q. Lu, and B. Zhao, “A multiobjective optimization of journal bearing with double parabolic profiles and groove textures under steady operating conditions,” Mechanics and Industry, vol. 21, no. 3, pp. 1–16, 2020.
View at: Publisher Site | Google Scholar
X. Y. Zhang, C. P. Liu, and B. Zhao, “An optimization research on groove textures of a journal bearing using particle swarm optimization algorithm,” Mechanics and Industry, vol. 22, no. 1, pp. 1–11, 2021.
View at: Publisher Site | Google Scholar
I. C. M. Filgueira Filho, A. C. Bottene, E. J. Silva, and R. Nicoletti, “Static behavior of plain journal bearings with textured journal - experimental analysis,” Tribology International, vol. 159, Article ID 106970, 2021.
View at: Google Scholar
F. M. Meng, L. Zhang, Y. Liu, and T. T. Li, “Effect of compound dimple on tribological performances of journal bearing,” Tribology International, vol. 91, no. 3, pp. 99–110, 2015.
View at: Publisher Site | Google Scholar
C. Gui and F. M. Meng, “Comparative study of spherical dimple and bump effects on the tribological performances of journal bearing,” Proceedings of the Institution of Mechanical Engineers - Part J: Journal of Engineering Tribology, vol. 233, no. 1, pp. 139–157, 2018.
View at: Publisher Site | Google Scholar
S. Sharma, A. Sharma, G. Jamwal, and R. K. Awasthi, “The effect of V-shape protruded and dimple textured on the load-carrying capacity and coefficient of friction of hydrodynamic journal bearing: a comparative numerical study,” Proceedings of the Institution of Mechanical Engineers - Part J: Journal of Engineering Tribology, vol. 235, no. 5, pp. 997–1011, 2020.
View at: Publisher Site | Google Scholar
K. Maharshi, T. Mukhopadhyay, B. Roy, L. Roy, and S. Dey, “Stochastic dynamic behaviour of hydrodynamic journal bearings including the effect of surface roughness,” International Journal of Mechanical Sciences, vol. 142-143, pp. 370–383, 2018.
View at: Publisher Site | Google Scholar
Y. Y. Jiang, Y. Gao, and Q. An, “Dynamic analysis of the rotor supported by journal bearing with micro-spherical surface texturing,” Proceedings of the Institution of Mechanical Engineers - Part K: Journal of Multi-Body Dynamics, vol. 228, no. 4, pp. 355–365, 2014.
View at: Publisher Site | Google Scholar
H. Yamada, S. Kaneko, and H. Taura, “Effect of square-dimple size on dynamic characteristics of textured journal bearings,” Proceedings of the 9th IFToMM international conference on rotor dynamics. Mechanisms and machine science, vol. 21, pp. 883–893, 2015.
View at: Publisher Site | Google Scholar
H. Yamada, H. Taura, and S. Kaneko, “Numerical and experimental analyses of the dynamic characteristics of journal bearings with square dimples,” Journal of Tribology, vol. 140, no. 1, Article ID 011703, 2018.
View at: Publisher Site | Google Scholar
S. Matele and K. N. Pandey, “Effect of surface texturing on the dynamic characteristics of hydrodynamic journal bearing comprising concepts of green tribology,” Proceedings of the Institution of Mechanical Engineers - Part J: Journal of Engineering Tribology, vol. 232, no. 11, pp. 1365–1376, 2018.
View at: Publisher Site | Google Scholar
S. Sharma, G. Jamwal, and R. K. Awasthi, “Numerical study on steady state performance enhancement of partial textured hydrodynamic journal bearing,” Industrial Lubrication & Tribology, vol. 71, no. 9, pp. 1055–1063, 2019.
View at: Publisher Site | Google Scholar
S. Sharma, G. Jamwal, and R. K. Awasthi, “Dynamic and stability performance improvement of the hydrodynamic bearing by using triangular-shaped textures,” Proceedings of the Institution of Mechanical Engineers - Part J: Journal of Engineering Tribology, vol. 234, no. 9, pp. 1436–1451, 2019.
View at: Publisher Site | Google Scholar
S. Sharma, G. Jamwal, and R. K. Awasthi, “Enhancement of steady state performance of journal bearing using chevron shaped surface texture,” Proceedings of the Institution of Mechanical Engineers - Part J: Journal of Engineering Tribology, vol. 233, no. 12, pp. 1833–1843, 2019.
View at: Publisher Site | Google Scholar
G. Jamwal, S. Sharma, and R. K. Awasthi, “The dynamic performance analysis of chevron shape textured hydrodynamic bearings,” Industrial Lubrication & Tribology, vol. 72, no. 1, pp. 1–8, 2019.
View at: Publisher Site | Google Scholar
F. M. Meng, Y. F. Zhang, L. L. Su, H. Y. Yu, and Y. Zheng, “Dynamic characteristics of compound textured journal bearing,” Proceedings of the Institution of Mechanical Engineers - Part J: Journal of Engineering Tribology, vol. 235, no. 7, pp. 1312–1334, 2020.
View at: Publisher Site | Google Scholar
N. Singh and R. K. Awasthi, “Influence of dimple location and depth on the performance characteristics of the hydrodynamic journal bearing system,” Proceedings of the Institution of Mechanical Engineers - Part J: Journal of Engineering Tribology, vol. 234, no. 9, pp. 1500–1513, 2020.
View at: Publisher Site | Google Scholar
N. Singh and R. K. Awasthi, “Influence of texture geometries on the performance parameters of hydrodynamic journal bearing,” Proceedings of the Institution of Mechanical Engineers - Part J: Journal of Engineering Tribology, vol. 235, pp. 1–17, 2021.
View at: Publisher Site | Google Scholar
T. Nanbu, R. Ning, Y. Yasuda, D. Zhu, and Q. J. Wang, “Micro-textures in concentrated conformal-contact lubrication: effects of texture bottom shape and surface relative motion,” Tribology Letters, vol. 29, no. 3, pp. 241–252, 2008.
View at: Publisher Site | Google Scholar
C. Shen and M. M. Khonsari, “Effect of dimple's internal structure on hydrodynamic lubrication,” Tribology Letters, vol. 52, no. 3, pp. 415–430, 2013.
View at: Publisher Site | Google Scholar
R. Rahmani, A. Shirvani, and H. Shirvani, “Optimization of partially textured parallel thrust bearings with square-shaped micro-dimples,” Tribology Transactions, vol. 50, no. 3, pp. 401–406, 2007.
View at: Publisher Site | Google Scholar
R. Rahmani, I. Mirzaee, A. Shirvani, and H. Shirvani, “An analytical approach for analysis and optimisation of slider bearings with infinite width parallel textures,” Tribology International, vol. 43, no. 8, pp. 1551–1565, 2010.
View at: Publisher Site | Google Scholar
M. Qiu, B. R. Minson, and B. Raeymaekers, “The effect of texture shape on the friction coefficient and stiffness of gas-lubricated parallel slider bearings,” Tribology International, vol. 67, no. 4, pp. 278–288, 2013.
View at: Publisher Site | Google Scholar
W. Wang, Y. Y. He, J. Zhao, Y. Li, and J. B. Luo, “Numerical optimization of the groove texture bottom profile for thrust bearings,” Tribology International, vol. 109, pp. 69–77, 2017.
View at: Publisher Site | Google Scholar
W. Wang, Y. Y. He, J. Zhao, J. Mao, Y. Hu, and J. Luo, “Optimization of groove texture profile to improve hydrodynamic lubrication performance: theory and experiments,” Friction, vol. 8, pp. 83–94, 2018.
View at: Publisher Site | Google Scholar
N. D. Hingawe and S. P. Bhore, “Multi-objective optimization of the design parameters of texture bottom profiles in a parallel slider,” Friction, vol. 8, pp. 726–745, 2020.
View at: Publisher Site | Google Scholar
D. Gropper, L. Wang, and T. J. Harvey, “Hydrodynamic lubrication of textured surfaces: a review of modeling techniques and key findings,” Tribology International, vol. 94, pp. 509–529, 2016.
View at: Publisher Site | Google Scholar
S. Xiong and Q. J. Wang, “Steady-state hydrodynamic lubrication modeled with the payvar-salant mass conservation model,” Journal of Tribology, vol. 134, no. 3, Article ID 031703, 2012.
View at: Publisher Site | Google Scholar
V. Bruno, M. Patrick, and F. Jean, “Cavitation in dynamically loaded journal bearings using mobility method,” Wear, vol. 193, no. 2, pp. 155–162, 1996.
View at: Google Scholar
J. W. Lund, “Review of the concept of dynamic coefficients for fluid film journal bearings,” Journal of Tribology, vol. 109, no. 1, pp. 37–41, 1987.
View at: Publisher Site | Google Scholar
R. Sinhasan and K. C. Goyal, “Transient response of a circular journal bearing lubricated with non-Newtonian lubricants,” Wear, vol. 156, no. 2, pp. 385–399, 1992.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Peng Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

598

Downloads

583

Citations