Next Article in Journal
VRChem: A Virtual Reality Molecular Builder
Next Article in Special Issue
Experimental and Numerical Dynamic Identification of an Aerostatic Electro-Spindle
Previous Article in Journal
Screening Support System Based on Patient Survey Data—Case Study on Classification of Initial, Locally Collected COVID-19 Data
Previous Article in Special Issue
Porous Gas Journal Bearings: An Exact Solution Revisited and Force Coefficients for Stable Rotordynamic Performance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parameter Sensitivity Analysis on Dynamic Coefficients of Partial Arc Annular-Thrust Aerostatic Porous Journal Bearings

1
Department of Automotive Engineering, Yeungnam University, Gyeongsan 38541, Korea
2
Melentiev Energy Systems Institute of Siberian Branch of Russian Academy of Sciences, 664033 Irkutsk, Russia
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2021, 11(22), 10791; https://doi.org/10.3390/app112210791
Submission received: 30 September 2021 / Revised: 10 November 2021 / Accepted: 11 November 2021 / Published: 15 November 2021
(This article belongs to the Special Issue Gas Bearings: Modelling, Design and Applications)

Abstract

:
Aerostatic bearings are widely used in high-precision devices. Partial arc annular-thrust aerostatic porous journal bearings are a prominent type of aerostatic bearings, which carry both radial and axial loads and provide high load-carrying capacity, low air consumption, and relatively low cost. Spindle shaft tilting is a resource-demanding challenge in numerical modeling because it involves a 3D air flow. In this study, the air flow problem was solved using a COMSOL software, and the dynamic coefficients for tilting degrees of freedom were obtained using finite differences. The obtained results exhibit significant coupling between the tilting motion in the x-and y-directions: cross-coupled coefficients can achieve 20% of the direct coefficient for stiffness and 50% for damping. In addition, a nonlinear behavior can be expected, because the tilting motion within 3°, tilting velocities within 0.0012°/s, and relative eccentricity of 0.2 have effects as large as 20% for direct stiffness and 100% for cross-coupled stiffness and damping. All dynamic coefficients were fitted with a polynomial of eccentricity, tilting, and tilting velocities in two directions, with a total of six parameters. The resulting fitting coefficient tables can be employed for the fast dynamic simulation of the rotor shaft carried on the proposed bearing type.

1. Introduction

Porous aerostatic bearings are widely used in machinery, owing to their relatively high load-carrying capacity and low air consumption [1,2,3]. Similar to other types of aerostatic bearings, owing to their high accuracy and zero contamination, they are suitable for high-precision devices. In addition, the friction drag of a porous aerostatic bearing is low [4,5], its motion accuracy is higher than that of a conventional orifice-type aerostatic bearing [5,6], and its load thresholds at high rotation speeds are the highest among other types of aerostatic bearings [5,7]. The production of porous materials with desirable permeabilities is relatively complicated. However, the amount of porous material required can be reduced via design optimization [8].
The optimal governing equations required for porous air bearing analysis remain under discussion. Zhong et al. [9] experimentally determined Ergun’s equation coefficients for the pressure drop in sintered metal porous media for air bearings. Belforte et al. [10] experimentally verified Forchheimer’s law for porous air bearing applications. Zhong et al. [11] verified the accuracy of Forchheimer’s law for air flow through sintered metal porous media and Darcy’s law, under slight pressure drops. Recently, it has been reported that a deep learning-based approach can be employed for high-precision and low-cost analysis of the 3D flow state inside a porous material [12].
Various numerical methods and engineering packages have been adopted for air pressure calculations in porous bearings. Huang et al. [13] analyzed the pressure in a porous conveyor air bearing using the FLUENT software. Cui et al. [14] used FLUENT software to numerically analyze the effect of manufacturing errors on the load-carrying capacity and stiffness of an annular-thrust porous aerostatic bearing. In addition, the deformation of porous thrust bearings was numerically investigated based on the fluid–solid interaction method by Wang et al. [15]. Van Ostayen et al. [16] used the COMSOL software for active aerostatic bearing analysis. Hwang and Khan [17] conducted 3D finite element analysis, while Plante et al. [18] proposed the shed 1D generalized flow theory. Otsu et al. [19] and Hosokawa et al. [8] used the FDM.
Several papers [19,20,21,22] deal with the numerical calculation of stiffness and damping coefficients corresponding to eccentricity or parallel shaft displacement (Figure 1a). However, insufficient data are available on stiffness and damping for shaft tilting (Figure 1b) in porous air bushing. Cui et al. [14,23] analyzed the non-flatness effect of the air-bearing surface on the stiffness and damping of an aerostatic porous bearing, which is computationally similar to modeling the tilting.
The rotor in ultraprecision devices is usually supported by journal and thrust bearings, which share a common air supply system but separate porous pads [14,23]. Annular-thrust bearings share the same porous pad carrying loads in both the annular and thrust directions (Figure 2). This dual role decreases the manufacturing cost while triggering flow coupling [23,24,25,26,27]. A partial arc bearing has a higher load-carrying capacity but lower stiffness than a full-arc bearing [27].
The numerical modeling of the rotor motion, including the effect of shaft tilting, as well as the modeling of a partial arc annular-thrust porous aerostatic bearing is substantially time-consuming, because it requires a 3D air flow solution. Dynamic modeling can be performed in a less resource-consuming manner if the full bearing model is replaced with stiffness and damping coefficients. The nonlinear dynamic motion can be modeled, provided the influences of eccentricity, tilting, and velocities on the dynamic coefficients are included. In this study, the air flow in the porous pad and air lubrication film between the porous pad and solid surface of the shaft were numerically modeled using the COMSOL Multiphysics® software. The pressure distribution and total forces and moments were calculated for low eccentricity, low tilting angles, and tilting speeds. The stiffness and damping coefficients were determined using finite differences.

2. Materials and Methods

The porous air bearing was modeled using the COMSOL Multiphysics software, “Darcy’s law”, and “Thin-Film Flow, Shell” modules. Here, the static solution is sufficient in determining the response of the displacement and velocity of the spindle. In addition, the porous bushing is a hollow cylinder parallel to the z-axis with a cut sector α. Figure 3 illustrates the different boundary conditions.
The pressure in the porous pad is modeled with the “Darcy’s law” module:
( ρ κ μ p   ) = 0
where p, ρ, κ , and μ denote the air pressure, air density, permeability coefficient, and dynamic viscosity of air, respectively. ∇ is the gradient operator in three-dimensional space. The pressure at the air-bearing surface (see Figure 3a) is modeled with the “thin-film flow, shell” module:
t ( ρ h ( 1 2 v t h 3 12 μ t p f ) ) = ρ t h + Q
where p f , h, v t , and Q represent the air film pressure, air film thickness, tangential spindle surface velocity, and air flux, respectively. t is the gradient operator in the two-dimensional space of the air film surface and ∇ is the partial derivative operator. At the interface between the porous pad and the air-bearing surface, the pressure and air flux must be continuous. In the “Darcy’s law” module, the “Pressure” boundary condition is adopted at the surfaces that open to the air lubrication film, which is expressed as:
p = p f .
In the “Thin-Film Flow, Shell”, the continuity of the air flux is expressed by admittance Y in the “Perforation” interface, which is expressed as:
Y = Q · n ρ   p f
where n = ( n x , n y , n z ) is the air-bearing surface normal vector. The air film thickness in the “Thin-Film Flow, Shell” module depends on the spindle eccentricity ( ε x ,   ε y ), tilting ( t x ,   t y ), and tilting velocities ( t ˙ x , t ˙ y ) (refer to Figure 1):
h = c + u · n u · t h
t h = v · n v · t h
u = ( c ε x + t x ( z L / 2 ) ; c ε y + t y ( z L / 2 ) ; t x x t y y )
v = ( y ω + t ˙ x ( z L / 2 ) ; x ω + t ˙ y ( z L / 2 ) ; t ˙ x x t ˙ y y )
where u is the spindle surface displacement, and v is the spindle surface velocity.
The isothermal compressibility of the air was modeled using the user-defined density in both modules. Here, the density is related to the air pressure in the “Variables” block as follows:
ρ = ρ 0   p p 0 ,
where the subscript 0 refers to the atmospheric conditions, and p0 = 1 bar. The air supply is modeled by another “Pressure” interface in the “Darcy’s law” module with a constant pressure of ps (see Figure 3b). The outflow at the air-bearing edges is given by the “Border” interface in the “Thin-Film Flow, Shell” module with a fixed pressure of p0 (see Figure 3c). The other boundaries were sealed.
After the static problem was solved and the air pressure was calculated, the moments in the x and y tilt directions were calculated, using the “Surface Integration” in the “Derived Values” block, according to the following expressions:
M x =   p ( n x ( z L / 2 ) n z x ) d x d y
M y =   p ( n y ( z L / 2 ) n z y ) d x d y
where x′ and y′ are the local coordinates of the air film surface.
Stiffness and damping for tilting degrees of freedom were calculated by using finite differences:
K t x , t x = ( M x ( ε x , ε y , t x + Δ t x , t y , t ˙ x , t ˙ y ) M x ( ε x , ε y , t x , t y , t ˙ x , t ˙ y ) ) / Δ t x
K t x , t y = ( M x ( ε x , ε y , t x , t y + Δ t y , t ˙ x , t ˙ y ) M x ( ε x , ε y , t x , t y , t ˙ x , t ˙ y ) ) / Δ t y
K t y , t x = ( M y ( ε x , ε y , t x + Δ t x , t y , t ˙ x , t ˙ y ) M y ( ε x , ε y , t x , t y , t ˙ x , t ˙ y ) ) / Δ t x
K t y , t y = ( M y ( ε x , ε y , t x , t y + Δ t y , t ˙ x , t ˙ y ) M y ( ε x , ε y , t x , t y , t ˙ x , t ˙ y ) ) / Δ t y
D t x , t x = ( M x ( ε x , ε y , t x , t y , t ˙ x + Δ t ˙ x , t ˙ y ) M x ( ε x , ε y , t x , t y , t ˙ x , t ˙ y ) ) / Δ t ˙ x
D t x , t y = ( M x ( ε x , ε y , t x , t y , t ˙ x , t ˙ y + Δ t ˙ y ) M x ( ε x , ε y , t x , t y , t ˙ x , t ˙ y ) ) / Δ t ˙ y
D t y , t x = ( M y ( ε x , ε y , t x , t y , t ˙ x + Δ t ˙ x , t ˙ y ) M y ( ε x , ε y , t x , t y , t ˙ x , t ˙ y ) ) / Δ t ˙ x
D t y , t y = ( M y ( ε x , ε y , t x , t y , t ˙ x , t ˙ y + Δ t ˙ y ) M y ( ε x , ε y , t x , t y , t ˙ x , t ˙ y ) ) / Δ t ˙ y

3. Results

The airbearing moments, tilting stiffness, and damping were calculated for ε x = 0.1 ,   0.2 ,   0.3 , ε y = 0.0 ,   0.1 , t x = 0.003 ° ,   0.00275 ° ,   ,   0.003 ° , and t y = 0.0 ° ,   0.00025 ° , , 0.003 ° , respectively; with p s = 8   bar , bearing length L = 0.11 m, shaft diameter D = 0.1 m, porous cylinder thickness b = 0.01 m, clearance c = 10 μm, rotation speed = 1000 rpm, bearing arc α = 240°, porosity = 0.4, and permeability = 8.33 × 10 17 m2. The film variation due to tilting in the given range does not exceed 57% for simultaneous change of t x and t y . Together with 30% variation due to eccentricity, it gives a minimal film thickness of 13% of clearance, or 1.3 μm. The air-bearing pressure, forces, and moments were calculated for all combinations of parameters t x , t y , ε x , and ε y , and their combinations with the following pairs of tilting velocity values: ( t ˙ x , t ˙ y ) { ( 0 , 0 ) , ( Δ t ˙ x , 0 ) , ( 2 Δ t ˙ x , 0 ) , ( 0 , Δ t ˙ y ) , ( 0 , 2 Δ t ˙ y ) } . The total number of cases was 13000. The steps in the differences (12–19) were Δ t x = Δ t y = 0.00025 ° , and Δ t ˙ x = Δ t ˙ y = 0.00573 ° / s . The maximum and average absolute errors of the finite differences are presented in Table 1.

3.1. Preliminary Analysis

As illustrated in Figure 4, the direct stiffness coefficient Ktx,tx has values in the range of 80–140 kN·m/rad. Within the range under consideration, Ktx,tx depends linearly on ε x and nonlinearly on t x , t y , t ˙ x , t ˙ y . The parameters t x , t ˙ x , and ε x exhibit significant effects (more than 10%). Ktx,tx is almost insensitive to t x at low negative tilting angles; however, it decreases linearly with t x at positive and high negative tilting angles (see Figure 4a). The effect of the x-tilting velocity t x ˙ is the diminishing of Ktx,tx, which is more significant at negative tilting angles. Ktx,tx increases with decreasing t y , increasing t ˙ y and ε x , and is insensible to ε y (see Figure 4b,c).
The cross-coupled stiffness coefficient Ktx,ty is approximately 10 times smaller than the direct stiffness coefficient Ktx,tx, and lies in the range of 6–15 kN·m/rad (see Figure 5). It is highly sensitive to all parameters, except for t ˙ x . Ktx,ty increases with decreasing t x and t ˙ y , and increasing t y (see Figure 5a,b). Eccentricity ε x can increase or decrease Ktx,ty, depending on ε y (see Figure 5c).
Kty,tx, which is another cross-coupled stiffness coefficient, is of the same order of magnitude; however, it is negative, and falls in the range of −18–0 kN·m/rad (refer to Figure 6). It is not sensitive to t ˙ x , and the t ˙ y effect on it is also negligible (see Figure 6a,b). The magnitude of Kty,tx is higher at lower t x and ε y and at higher t y   and ε x .
The direct stiffness coefficient Kty,ty is higher than that of Ktx,tx. Kty,ty falls in the range of 110–140 kN·m/rad (refer to Figure 7). This is triggered by the cut arc, which increases force and decreases the stiffness in the x-direction [27]. The y-tilting angle t y exerts a strong nonlinear effect on Kty,ty. The stiffness coefficient Kty,ty is almost constant for a low positive t y , and then increases with increasing t y . In contrast to Ktx,tx, which is primarily sensitive to “its own” parameter t x   and significantly less sensitive to t y . Kty,ty is equally sensitive to both t x   and t y . In addition, it is equally sensitive to both higher ε x and ε y , with a negligible t ˙ y effect, and no t ˙ x effect. Furthermore, it is higher for lower t x and higher t y ,   ε x , and ε y .
Because the simultaneous variation of tilting velocities was not studied, the effect of t ˙ x was solely studied for Dtx,tx and Dty,tx, while the effect of t ˙ y was solely studied for Dtx,ty and Dty,ty. Instead of the missing parameter, the coupling between t x   and t y is comprehensively presented in Figure 8, Figure 9, Figure 10 and Figure 11.
The direct damping coefficient Dtx,tx is positive and falls in the range 0.2–2.2 kN·m s/rad (refer to Figure 8). It is higher for lower t x , and higher t y and higher t ˙ x (see Figure 8a,b). The effect of eccentricity is also substantial, nonlinear, and strongly coupled between ε x and ε y (see Figure 8c).
The cross-coupled damping coefficient Dtx,ty does not exceed 1 kN m s/rad in magnitude, being mostly negative (see Figure 9). Dtx,ty rapidly decreases with decreasing t x or increasing t y in the negative t x   range, while in the positive t x   range, it is much less sensible both to t x and t y (see Figure 9a). Dtx,ty decreases with increasing t y and is almost insensible to t ˙ y (see Figure 9b). The dependence of Dtx,ty on eccentricity parameters is highly nonlinear. Dtx,ty can take relatively high positive values for zero ε x . and positive ε y , but in most of the range, the sensitivity to ε x and ε y is low (see Figure 9c).
The cross-coupled damping coefficient Dty,tx falls in the range −0.6–0.13 kN m s/rad (see Figure 10). It is negative for positive t y and negative t x , otherwise it is positive. Similar to Dtx,ty, Dty,tx rapidly decreases with decreasing t x or increasing t y in the negative t x   range, while in the positive t x   range, it is much less sensible both to t x and t y . However, for t x > 0.002 °, Dty,tx increases with increasing t y (see Figure 10b), while Dtx,ty decreases with increasing t y for any t x (see Figure 9a). For zero t y ,  Dty,tx decreases with increasing t x and decreasing t ˙ x . Similar to all previously considered damping coefficients, the dependence of Dty,tx on the eccentricity parameters is highly nonlinear.
In contrast with all other damping coefficients, the direct damping coefficient Dty,ty does not significantly depend on ε x and ε y (refer to Figure 11). It falls in the range 0.6–2.0 kN·m s/rad. It is higher for lower t x , and higher t y , similar to Dtx,tx. Dty,ty is almost insensible to the y-tilt velocity t ˙ y . It slightly increases with increasing ε x and increasing ε y .
In summary, the direct stiffness and damping coefficients were positive. The cross-coupled stiffness coefficients are approximately 10 times smaller than the direct stiffness coefficients, and the cross-coupled damping coefficients are approximately 2–3 times smaller than the direct damping coefficients. Hence, a stable tilting motion was expected. However, information on the rotor inertia and the bearing at the other end is required to make precise conclusions on the dynamic behavior. Moreover, the results above exhibit a strong coupling between the eccentricity and tilting motion, including the nonlinearity of the moments relative to the position parameters. This can make the rotor motion severely complicated.

3.2. Regression

Although it is possible to obtain the direct solution of the time-dependent problem in the model described above, the process is significantly time consuming because the air flow problem must be solved at every step. Conversely, the air-bearing forces and moments can be calculated using the interpolated stiffness and damping coefficients. In this study, each dynamic coefficient for each combination of ε x , ε y , t ˙ x , and t ˙ y was first fitted with a quadratic polynomial of t x and t y , and then each coefficient of the polynomial was fitted with a linear function of   ε x , ε y , t ˙ x , and t ˙ y . The accuracy of the fitting is expressed in the form of determination coefficients, as presented in Table 2. The accuracy is relatively good, although it can be improved for K ˜ t x , t y ,   D ˜ t x , t x , D ˜ t x , t y , and   D ˜ t y , t x by providing more points and increasing the fitting order in ε x and ε y degrees of freedom.
The fitting results are presented in Table 3, Table 4, Table 5, Table 6, Table 7, Table 8, Table 9 and Table 10. An example of the table interpretation is given below for K t x , t x :
K ˜ t x , t x = ( 1.09533 + 1.49129 ε x + 0.17719 ε y 0.02161 T ˙ x ) + ( 0.30745 3.94666 ε x + 0.10915 T ˙ x + 0.00478 T ˙ y ) T x + ( 0.15003 + 0.09324 ε x 1.16823 ε y + 0.01582 T ˙ x + 0.0114 T ˙ y ) T y + ( 0.3232 + 2.45267 ε x + 2.42515 ε y + 0.16265 T ˙ x 0.28434 T ˙ y ) T x T y + ( 7.99445 ε x 0.26854 ε y 0.35236 T ˙ x + 0.14364 T ˙ y ) T x 2 + ( 0.34256 + 0.41188 ε x + 1.09915 ε y 0.13117 T ˙ x + 0.17535 T ˙ y ) T y 2
Here, the dimensionless variables are adopted:
K i , j K ˜ i , j = p 0 D L 3 16   c = 83 , 187   N   m rad
D i , j D ˜ i , j = p 0 D L 2 2   ( 0.0001 rad s ) = 1.51 × 10 5   N   m   s rad
t x T x = t y T y = c L = 0.000182   rad = 0.0104 °
t ˙ x T ˙ x = t ˙ y T ˙ y = 0.0001 rad s = 0.00573 ° 1 s

4. Discussion

The tilting stiffness and damping analyzed in this study could be significant for the spindle shaft tilting motion, provided the shaft length is not significantly longer than the bearing length. Otherwise, the shaft tilting will be controlled by the air-bearing force applied at a long distance from the shaft center, while the tilting stiffness and damping will influence the shaft bending. In addition, bending can also significantly influence the air-bearing forces, moments, and dynamic coefficients, because of the convexity and concavity of the bearing surface [14,23]. Based on the obtained results, the tilting motion is independently stable because the direct stiffness and damping coefficients are positive and significantly higher in amplitude than the cross-coupled coefficients. However, the tilting dynamic coefficients were substantially influenced by shaft eccentricity. Hence, the stability of the shaft parallel and tilting motions must be analyzed together. The increased rotor speed would result in a more pronounced hydrodynamic component in the bearing forces, moments, stiffness, and damping components, as well as in more significant x-y coupling, that is typical for the journal bearings.

5. Conclusions

The problem of tilting stiffness and damping analysis is hardly ever addressed in the literature since the fundamental books on rotodynamic and tribology where the theoretical analysis was presented for idealized types of bearings. This ignorance can be explained by the minor role of individual bearing moments, compared to the total shaft-bearing systems moments in most applications. However, the tilting moments exist, and they become more and more important in high-precision machinery. Besides the bearing moments, the tilting motion greatly affects the bearing forces, as well as other types of air film distortions, which are widely considered in recent publications.
The purpose of this study was to predict the behavioral characteristics according to the design and operating parameters of porous aerostatic bearings that can be utilized in applications requiring high-precision positioning, such as semiconductor manufacturing equipment. Most importantly, the highly efficient and reliable numerical analysis approach proposed in this study can contribute to increasing system stability and reducing costs and risks at the product level prior to the actual manufacturing process.

Author Contributions

Conceptualization, P.H. and P.K.; methodology, P.K.; software, P.K.; validation, P.H. and P.K.; formal analysis, P.K.; investigation, P.H. and P.K.; resources, P.H. and P.K.; data curation, P.H.; writing—original draft preparation. P.K.; writing—review and editing. P.H. and S.-W.K.; visualization, P.K.; supervision, P.H.; project administration, P.H.; funding acquisition, P.H. and S.-W.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partly supported by Korea Institute for Advancement of Technology (KIAT) grant funded by the Korea Government (MOTIE) (P0002092, The Competency Development Program for Industry Specialist) and Korea Institute of Energy Technology Evaluation and Planning (KETEP) grant funded by the Korea Government (MOTIE) (20214000000010, Gyeongbuk Regional Wind Energy Cluster Human Resources Development Project) and State Assignment Project (No. FWEU-2021-0005) of the Fundamental Research Program of the Russian Federation 2021–2030.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available by request.

Acknowledgments

The research was performed with use of the equipment of High-Temperature Circuit Multi-Access Research Center (Project No. 13.ЦKП.21.0038).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Akilian, M.; Forest, C.R.; Slocum, A.H.; Trumper, D.L.; Schattenburg, M.L. Thin optic constraint. Precis. Eng. 2007, 31, 130–138. [Google Scholar] [CrossRef]
  2. New Way Air Bearings. Air Bearing Application and Design Guide; IBS Precision Engineering: Eindhoven, The Netherlands, 2006. [Google Scholar]
  3. Khim, G.; Park, C.H. A rotary stage in a vacuum using an air bearing. Vacuum 2014, 105, 39–45. [Google Scholar] [CrossRef]
  4. Andres, L.S.; Cable, T.A.; Zheng, Y.; De Santiago, O.; Devitt, D. Assessment of Porous Type Gas Bearings: Measurements of Bearing Performance and Rotor Vibrations. In Proceedings of the ASME Turbo Expo 2016: Turbomachinery Technical Conference and Exposition, Seoul, Korea, 13–17 June 2016. No. V07BT31A031. [Google Scholar]
  5. Gao, Q.; Chen, W.; Lu, L.; Huo, D.; Cheng, K. Aerostatic bearings design and analysis with the application to precision engineering: State-of-the-art and future perspectives. Tribol. Int. 2019, 135, 1–17. [Google Scholar] [CrossRef]
  6. Cappa, S.; Reynaerts, D.; Al-Bender, F. Reducing the Radial Error Motion of an aerostatic Journal Bearing to a Nanometer Level: Theoretical Modelling. Tribol. Let. 2014, 53, 27–41. [Google Scholar] [CrossRef]
  7. Lie, K.N.; Su, J.C. A numerical analysis of optimum air journal bearings. J. Aeronaut. Astronaut. Aviation. Ser. A 2014, 46, 49–60. [Google Scholar]
  8. Hosokawa, T.; Somaya, K.; Miyatake, M.; Yoshimoto, S. Static characteristics of aerostatic thrust bearings with multiple porous inlet ports. J. Tribol. 2015, 137, 1–8. [Google Scholar] [CrossRef]
  9. Zhong, W.; Li, X.; Tao, G.L.; Kagawa, T. Measurement and Determination of Friction Characteristic of Air Flow through Porous Media. Metals 2015, 5, 336–349. [Google Scholar] [CrossRef]
  10. Belforte, G.; Raparelli, T.; Viktorov, V.; Trivella, A. Permeability and Inertial Coefficients of Porous Media for Air Bearing Feeding Systems. ASME J. Tribol. 2007, 129, 705–711. [Google Scholar] [CrossRef]
  11. Zhong, W.; Li, X.; Liu, F.H.; Tao, G.L.; Lu, B.; Kagawa, T. Measurement and Correlation of Pressure Drop Characteristics for Air Flow Through Sintered Metal Porous Media. Transp. Porous Media 2014, 101, 53–67. [Google Scholar] [CrossRef]
  12. Zeng, C.; Wang, W.; Cheng, X.; Zhao, R.; Cui, H. Three-dimensional flow state analysis of microstructures of porous graphite restrictor in aerostatic bearings. Tribol. Int. 2021, 159, 106955. [Google Scholar] [CrossRef]
  13. Huang, T.Y.; Shen, S.C.; Lin, S.C.; Hsu, S.Y. Pressure Distribution in the Air Film and the Porous Conveyor Air Bearing. Appl. Mech. Mater. 2014, 479, 380–384. [Google Scholar] [CrossRef]
  14. Cui, H.; Wang, Y.; Yue, X.; Huang, M.; Wang, W. Effects of manufacturing errors on the static characteristics of aerostatic journal bearings with porous restrictor. Tribol. Int. 2017, 115, 246–260. [Google Scholar] [CrossRef]
  15. Wang, W.; Cheng, X.; Zhang, M.; Gong, W.; Cui, H. Effect of the deformation of porous materials on the performance of aerostatic bearings by fluid-solid interaction method. Tribol. Int. 2020, 150, 106391. [Google Scholar] [CrossRef]
  16. Van Ostayen, R.A.J.; Van Beek, A.; Munnig-Schmidt, R. Design and Optimization of an Active Aerostatic Thrust Bearing. In Proceedings of the ASPE, Dallas, TX, USA, 14–19 October 2007. [Google Scholar]
  17. Hwang, P.; Khan, P.V. Non-Iterative Finite Element Scheme for Combined Annular-Thrust Porous Aerostatic Bearings Analysis. In Proceedings of the ASME/STLE International Joint Tribology Conference, Denver, CO, USA, 7–10 October 2017. [Google Scholar]
  18. Plante, J.S.; Vogan, J.; El-Aguizy, T.; Slocum, A.H. A design model for circular porous air bearings using the 1D generalized flow method. Prec. Eng. ISPEN J. 2005, 29, 336–346. [Google Scholar] [CrossRef]
  19. Otsu, Y.; Miyatake, M.; Yoshimoto, S. Dynamic Characteristics of Aerostatic Porous Journal Bearings with a Surface Restricted Layer. ASME J. Tribol. 2011, 133, 011701-1. [Google Scholar] [CrossRef]
  20. Zhang, X.; Lin, B. Theoretical research on deformation of porous material in air bearing. Appl. Mech. Mater. 2012, 215–216, 779–784. [Google Scholar] [CrossRef]
  21. Yoshimoto, S.; Kohno, K. Static and Dynamic Characteristics of Aerostatic Circular Porous Thrust Bearings (Effect of the Shape of the Air Supply Area). ASME J. Tribol. 2001, 123, 501–508. [Google Scholar] [CrossRef]
  22. Su, J.C.T.; Lie, K.N. Rotor dynamic instability analysis on hybrid air journal bearings. Tribol. Int. 2006, 39, 238–248. [Google Scholar] [CrossRef]
  23. Cui, H.; Wang, Y.; Yue, X.; Li, Y.; Jiang, Z. Numerical analysis of the dynamic performance of aerostatic thrust bearings with different restrictors. Proc. Inst. Mech. Eng. Part J. Eng. Tribol. 2019, 233, 406–423. [Google Scholar] [CrossRef]
  24. Zang, Y.; Hatch, M.R. Analysis of Coupled Journal and Thrust Hydrodynamic Bearing using Finite Volume Method. ASME Adv. Inf. Storage Process. Syst. 1995, 1, 71–79. [Google Scholar]
  25. Rahman, M.; Leuthold, H. Computer simulation of a coupled journal and thrust hydrodynamic bearing using a Finite-Element Method. ASME Adv. Inf. Storage Process. Syst. 1996, 2, 103. [Google Scholar]
  26. Wang, Y.; Wang, Q.J.; Lin, C. Mixed Lubrication of Coupled Journal-Thrust-Bearing Systems Including Mass Conserving Cavitation. ASME J. Tribol. 2003, 125, 747–755. [Google Scholar] [CrossRef]
  27. Hwang, P.; Khan, P.V. Static Study of the Arc Effect on Annular-Thrust Aerostatic Porous Bearings. In Proceedings of the 5th World Tribology Congress (WTC 2013), Torino, Italy, 8–13 September 2013. [Google Scholar]
Figure 1. Two spindle motion modes: (a) eccentricity; (b) tilting. sb and sr denote the bearing and rotor axes, respectively.
Figure 1. Two spindle motion modes: (a) eccentricity; (b) tilting. sb and sr denote the bearing and rotor axes, respectively.
Applsci 11 10791 g001
Figure 2. Annular-thrust porous aerostatic bearing.
Figure 2. Annular-thrust porous aerostatic bearing.
Applsci 11 10791 g002
Figure 3. Porous pad boundaries: (a) airbearing surface carrying both radial and axial loads; (b) air supply surface; (c) outflow to atmospheric pressure; (d) outflow to atmospheric pressure in xy-plane with axes and bearing partial arc.
Figure 3. Porous pad boundaries: (a) airbearing surface carrying both radial and axial loads; (b) air supply surface; (c) outflow to atmospheric pressure; (d) outflow to atmospheric pressure in xy-plane with axes and bearing partial arc.
Applsci 11 10791 g003
Figure 4. Parameter sensitivity of the direct stiffness coefficient Ktx,tx: (a) Effects of the tilting angle and tilting velocity in the x-direction; (b) Effects of the tilting angle and tilting velocity in the y-direction; (c) Effects of relative eccentricities. All unreferenced parameters have zero values.
Figure 4. Parameter sensitivity of the direct stiffness coefficient Ktx,tx: (a) Effects of the tilting angle and tilting velocity in the x-direction; (b) Effects of the tilting angle and tilting velocity in the y-direction; (c) Effects of relative eccentricities. All unreferenced parameters have zero values.
Applsci 11 10791 g004
Figure 5. Parameter sensitivity of the cross-coupled stiffness coefficient Ktx,ty: (a) Effects of tilting angle and tilting velocity in the x-direction; (b) Effects of the tilting angle and tilting velocity in the y-direction; (c) Effects of relative eccentricities. All unreferenced parameters have zero values.
Figure 5. Parameter sensitivity of the cross-coupled stiffness coefficient Ktx,ty: (a) Effects of tilting angle and tilting velocity in the x-direction; (b) Effects of the tilting angle and tilting velocity in the y-direction; (c) Effects of relative eccentricities. All unreferenced parameters have zero values.
Applsci 11 10791 g005
Figure 6. Parameter sensitivity of the cross-coupled stiffness coefficient Kty,tx: (a) Effects of the tilting angle and tilting velocity in the x-direction; (b) Effects of the tilting angle and tilting velocity in the y-direction; (c) Effects of relative eccentricities. All unreferenced parameters have zero values.
Figure 6. Parameter sensitivity of the cross-coupled stiffness coefficient Kty,tx: (a) Effects of the tilting angle and tilting velocity in the x-direction; (b) Effects of the tilting angle and tilting velocity in the y-direction; (c) Effects of relative eccentricities. All unreferenced parameters have zero values.
Applsci 11 10791 g006
Figure 7. Parameter sensitivity of the direct stiffness coefficient Kty,ty: (a) Effects of the tilting angle and tilting velocity in the x-direction; (b) Effects of the tilting angle and tilting velocity in the y-direction; (c) Effects of relative eccentricities. All unreferenced parameters have zero values.
Figure 7. Parameter sensitivity of the direct stiffness coefficient Kty,ty: (a) Effects of the tilting angle and tilting velocity in the x-direction; (b) Effects of the tilting angle and tilting velocity in the y-direction; (c) Effects of relative eccentricities. All unreferenced parameters have zero values.
Applsci 11 10791 g007
Figure 8. Parameter sensitivity of the direct damping coefficient Dtx,tx: (a) Effects of the tilting angle and tilting velocity in the x-direction; (b) Effects of the x and y tilting angles; (c) Effects of the relative eccentricities. All unreferenced parameters have zero values.
Figure 8. Parameter sensitivity of the direct damping coefficient Dtx,tx: (a) Effects of the tilting angle and tilting velocity in the x-direction; (b) Effects of the x and y tilting angles; (c) Effects of the relative eccentricities. All unreferenced parameters have zero values.
Applsci 11 10791 g008
Figure 9. Parameter sensitivity of the cross-coupled damping coefficient Dtx,ty: (a) Effects of x and y tilting angles; (b) Effects of the tilting angle and tilting velocity in the y-direction; (c) Effects of relative eccentricities. All unreferenced parameters have zero values.
Figure 9. Parameter sensitivity of the cross-coupled damping coefficient Dtx,ty: (a) Effects of x and y tilting angles; (b) Effects of the tilting angle and tilting velocity in the y-direction; (c) Effects of relative eccentricities. All unreferenced parameters have zero values.
Applsci 11 10791 g009
Figure 10. Parameter sensitivity of the cross-coupled damping coefficient Dty,tx: (a) Effects of the tilting angle and tilting velocity in the x-direction; (b) Effects of x and y tilting angles; (c) Effects of relative eccentricities. All unreferenced parameters have zero values.
Figure 10. Parameter sensitivity of the cross-coupled damping coefficient Dty,tx: (a) Effects of the tilting angle and tilting velocity in the x-direction; (b) Effects of x and y tilting angles; (c) Effects of relative eccentricities. All unreferenced parameters have zero values.
Applsci 11 10791 g010
Figure 11. Parameter sensitivity of the direct damping coefficient Dty,ty: (a) Effects of x and y tilting angles; (b) Effects of the tilting angle and tilting velocity in the y-direction; (c) Effects of relative eccentricities. All unreferenced parameters have zero values.
Figure 11. Parameter sensitivity of the direct damping coefficient Dty,ty: (a) Effects of x and y tilting angles; (b) Effects of the tilting angle and tilting velocity in the y-direction; (c) Effects of relative eccentricities. All unreferenced parameters have zero values.
Applsci 11 10791 g011
Table 1. Finite difference errors.
Table 1. Finite difference errors.
Stiffness (kN·m/rad)Damping (kN·m·s/rad)
Ktx,txKtx,tyKty,txKty,tyDtx,txDtx,tyDty,txDty,ty
max0.1240.0660.0890.1200.0110.0080.0030.003
average0.0220.0060.0100.0160.0030.0020.00040.0003
Table 2. Coefficients of determination for dimensionless dynamic coefficient fitting.
Table 2. Coefficients of determination for dimensionless dynamic coefficient fitting.
Stiffness (kN·m/rad)Damping (kN·m·s/rad)
K ˜ t x , t x K ˜ t x , t y K ˜ t y , t x K ˜ t y , t y D ˜ t x , t x D ˜ t x , t y D ˜ t y , t x D ˜ t y , t y
R20.990.770.970.960.890.830.910.96
Table 3. Fitting results for K ˜ t x , t x .
Table 3. Fitting results for K ˜ t x , t x .
1 ε x ε y T ˙ x T ˙ y
11.095331.491290.17719−0.021610
T x −0.30745−3.9466600.109150.00478
T y −0.150030.09324−1.168230.015820.0114
T x T y −0.32322.452672.425150.16265−0.28434
T x 2 07.99445−0.26854−0.352360.14364
T y 2 0.342560.411881.09915−0.131170.17535
Table 4. Fitting results for K ˜ t x , t y .
Table 4. Fitting results for K ˜ t x , t y .
1 ε x ε y T ˙ x T ˙ y
10.1097300.51941−0.00802−0.02162
T x −0.13722−0.006−0.807450−0.02419
T y 0−1.68577−0.422860.080250.193
T x T y 0.229043.62763−2.1016400.74463
T x 2 −0.12131.020971.429660.06678−0.15351
T y 2 0.500933.62374−2.31180−0.34029
Table 5. Fitting results for K ˜ t y , t x .
Table 5. Fitting results for K ˜ t y , t x .
1 ε x ε y T ˙ x T ˙ y
1−0.07381−0.096730.78029−0.00428−0.00686
T x 0.110420.41995−1.260020.014490.00965
T y −0.24779−1.110430.334550.020620.02685
T x T y 1.066052.902280−0.208870.24493
T x 2 −0.72174−0.864341.501620.07247−0.11828
T y 2 −0.61270.695470.933970.10771−0.20109
Table 6. Fitting results for K ˜ t y , t y .
Table 6. Fitting results for K ˜ t y , t y .
1 ε x ε y T ˙ x T ˙ y
11.425690.479050.2985100.0155
T x −0.30423−1.1056500.06510.12299
T y −0.642440.25697−3.214670−0.10432
T x T y −1.0113105.084340.17214−0.52047
T x 2 0.432631.335880−0.033280.27025
T y 2 3.536341.579832.17612−0.132590.69417
Table 7. Fitting results for D ˜ t x , t x .
Table 7. Fitting results for D ˜ t x , t x .
1 ε x ε y T ˙ x T ˙ y
10.00970.0035800.001560
T x −0.02028−0.050930−0.001170
T y −0.00168−0.00288−0.016830.000230
T x T y −0.01595−0.043940.092860.001580
T x 2 0.034040.194380−0.002250
t y 2 0.017080.04876−0.01605−0.00150
Table 8. Fitting results for D ˜ t x , t y .
Table 8. Fitting results for D ˜ t x , t y .
1 ε x ε y T ˙ x T ˙ y
10.0024−0.009620.012650−0.00171
T x −0.002990.01579−0.0265400.00227
T y −0.00896−0.01581000.00029
T x T y 0.036540.10701−0.0510300.00167
T x 2 −0.00931−0.052370.054140−0.00165
T y 2 −0.01297−0.031910.062420−0.0018
Table 9. Fitting results for D ˜ t y , t x .
Table 9. Fitting results for D ˜ t y , t x .
1 ε x ε y T ˙ x T ˙ y
10.000720.001140.0067400
T x −0.00151−0.00441−0.0180200
T y −0.008330−0.01345−0.002080
T x T y 0.033250.08433−0.0282800
T x 2 −0.00695−0.015930.0464100
T y 2 −0.00795−0.015520.04750.000930
Table 10. Fitting results for D ˜ t y , t y .
Table 10. Fitting results for D ˜ t y , t y .
1 ε x ε y T ˙ x T ˙ y
10.009870.004260.002890−0.00014
T x −0.00619−0.01127−0.001700.00026
T y −0.009430.01182−0.0492700.00208
T x T y −0.01952−0.050930.103360−0.00248
T x 2 0.015140.04059−0.0194100.0007
T y 2 0.069670.05133−0.0767600.00259
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hwang, P.; Khan, P.; Kang, S.-W. Parameter Sensitivity Analysis on Dynamic Coefficients of Partial Arc Annular-Thrust Aerostatic Porous Journal Bearings. Appl. Sci. 2021, 11, 10791. https://doi.org/10.3390/app112210791

AMA Style

Hwang P, Khan P, Kang S-W. Parameter Sensitivity Analysis on Dynamic Coefficients of Partial Arc Annular-Thrust Aerostatic Porous Journal Bearings. Applied Sciences. 2021; 11(22):10791. https://doi.org/10.3390/app112210791

Chicago/Turabian Style

Hwang, Pyung, Polina Khan, and Seok-Won Kang. 2021. "Parameter Sensitivity Analysis on Dynamic Coefficients of Partial Arc Annular-Thrust Aerostatic Porous Journal Bearings" Applied Sciences 11, no. 22: 10791. https://doi.org/10.3390/app112210791

APA Style

Hwang, P., Khan, P., & Kang, S. -W. (2021). Parameter Sensitivity Analysis on Dynamic Coefficients of Partial Arc Annular-Thrust Aerostatic Porous Journal Bearings. Applied Sciences, 11(22), 10791. https://doi.org/10.3390/app112210791

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

Article Metrics

Back to TopTop