Research Article  Open Access
Numerical Simulation and Experimental Study on the GasSolid Coupling of the Aerostatic Thrust Bearing with Elastic Equalizing Pressure Groove
Abstract
Aiming at the problem of low stiffness of aerostatic bearing, according to the principle of gassolid coupling, this paper designs a kind of aerostatic thrust bearing with elastic equalizing pressure groove (EEPG) and investigates the effect of elastic equalizing pressure groove (EEPG) on the stiffness of aerostatic bearing. According to the physical model of the bearing, one deduces the deformation control equation of the elastic equalizing pressure groove and the control equation of gas lubrication, using finite difference method to derive the control equations and coupling calculation. The bearing capacity and stiffness of aerostatic bearing with EEPG in different gas film clearance are obtained. The calculation results show that the stiffness increased by 59%. The results of numerical calculation and experimental results have good consistency, proving the gassolid coupling method can improve the bearing stiffness.
1. Introduction
Aerostatic bearings are widely used in the field of high speed and high accuracy due to their good performance [1, 2]. However, because the air has obvious compressibility, aerostatic bearing is difficult to obtain high stiffness, so that its application field is confined to the condition of light load. Therefore, it is very important to improve the stiffness of the aerostatic bearing and expand its application field.
At present, there are two methods to improve the stiffness of aerostatic bearings: one way is the change of geometric structure of aerostatic bearing, increasing the supply pressure, increasing the throttle hole quantity, decreasing throttle hole diameter, and setting the equalizing pressure groove and a throttling cavity with different shapes on the aerostatic bearing surface to obtain a high stiffness [3–5]. Nishio et al. [6] studied the effect of orifice diameter less than 0.05 mm on the performance of aerostatic thrust bearings. They found that the smaller the diameter of the orifice, the greater the stiffness of the bearing. Hao and Jia [7] studied the effect of specialshaped throttling cavity on the performance of aerostatic bearing. Charki et al. [8] studied the performance of aerostatic bearings with multiple orifices. The results show that the performance of the bearing is not only related to the flow rate but also closely related to the number and diameter of the orifices. Belforte et al. [9] studied the influence of the bearing groove on the bearing pressure distribution, gas consumption, and stiffness. This work was also done in some studies [10–12]. Another way is through the electric or pneumatic precision electronic components to control the throttle effect of aerostatic bearings, in order to obtain high stiffness [13]. Aguirre et al. study [14] is driven by an electromagnetic drive mechanism to control the gas film clearance changes, so as to improve the bearing stiffness. Ghodsiyeh et al. [15] adopt pneumatic control as a semiactive compensation and improve the aerodynamic stiffness of the air bearing, and the main parameters affecting the static performance of the system are discussed. Compared with the above two methods, with only changing the geometry of aerostatic bearings, the structure is simple, but it is difficult to obtain a larger stiffness. Compared with the first kind of research, the second research method, through active control and dynamic compensation, can get larger stiffness, but this kind of aerostatic bearing has complex structure, with high degree of mechanical and electrical integration.
Through the above analysis, it can be seen that the gas film clearance is an important factor that affects the stiffness of aerostatic bearing. In this paper, the gassolid coupling method is proposed to enhance the stiffness of aerostatic bearing with EEPG, and the gas film clearance can be changed dynamically through the deformation of EEPG. The mathematical model of the gassolid coupling of EEPG is derived. The numerical simulation and experimental results show that the EEPG can improve the stiffness of aerostatic bearing.
2. Physical Model and Calculation Method
2.1. Working Principle
This section mainly introduces the twodimensional physical model and its working principle. This new type of bearing consists of two parts: the bearing body and the elastic thin plate. The elastic thin plate is connected with the bearing body through vacuum brazing, and the EEPG is mainly obtained by the deflection deformation of elastic thin plate in this paper. The new bearing structure is shown in Figure 1. The radius is R, there is an annular airway in the bearings of the internal part, and the elastic thin plate is covered on the annular airway. There are four orifices that are evenly distributed on the elastic thin plate, and the diameter is . The minimum radius of the annular airway is , and the maximum radius is . The width value of the annular airway is the difference between the maximum radius and the minimum radius, that is, .
Figure 2(a) is the bearing of the EEPG processing schematic. Firstly, high pressure gas is injected into the bearing to make the elastic deformation of the elastic thin plate, as shown in Figure 2(a)(). And then grind the convex part of the elastic thin plate, as shown in Figures 2(a)() and 2(a)(). Then stop gas supply, the deformation is restored, elastic plate will form a groove in an annular passage, and the concave shallow groove is the initial EEPG.
(a) Manufacture process
(b) Working state
The bearing is working, the elastic thin plate is affected by the pressure difference between the supply pressure and the gas film pressure, and thin elastic deformation is generated, thus forming an EEPG. The gas film clearance is small, with gas film pressure distribution near the supply pressure, and the pressure difference of thin plate is small, so the basic elastic thin plate is not deformed. Elastic equalizing pressure groove depth of aerostatic bearing is substantially equal to the initial processing produce. At this time, the working state is equal to the working state of the aerostatic bearing of the REPG (rigid equalizing pressure groove). The gas film clearance increases to a certain range, the gas film pressure is less than the supply pressure. The elastic thin plate pressure difference increases, resulting in deformation of thin plate, but less than the depth of the initial equalizing pressure groove. At the moment, the depth of the EEPG is equal to initial elastic equalizing pressure groove minus the deformation of the thin plate. The gas film clearance increases to limit clearance, with distribution of film pressure close to standard atmospheric pressure. At the same time, the elastic thin plate is only affected by the supply pressure, the EEPG disappears, and the equalizing pressure effect disappears. Figure 2(b) shows the different working states of bearings.
It is worth noting that the gas supply pressure at work should be greater than processing pressure (). In this work, the EEPG of the new aerostatic bearing will change dynamically. With the increase of the gas film clearance, the bearing capacity of the bearing decreases rapidly, and the stiffness formula , and the faster the decrease of bearing capacity, the bigger the bearing stiffness. Compared with the aerostatic bearing with REPG, the aerostatic bearing with EEPG has better adaptability and higher stiffness.
2.2. The Governing Equation
2.2.1. Elastic Equalizing Pressure Groove (EEPG) Deformation Control Equation
According to the assumption of small deflection bending theory of thin plate [16], the basic equation of strain component of bending deformation can be obtained, as shown in and herewhere , and are components of strain and is Bending deflection.
The stress component is obtained by (3):
That is,where is Young’s modulus and is Poisson ratio.
Taking a parallelepiped element from the elastic thin plate, each edge length is . On the cross section which perpendicular to the axis. The stress component is synthesized as a bending moment; the stress component is synthesized as torque; the stress component is synthesized as a transverse shear force, as shown in and herewhere is Bending stiffness and is elastic thin plate thickness.
By the same argument, on the cross section which is perpendicular to the axis, the stress component is synthesized as a bending moment; the stress component is synthesized as torque; the stress component is synthesized as a transverse shear force, as shown in
The equilibrium equation of thin plate is
From (5), (7), and (8) one can get the deformation control equation, as shown in where is uniform pressure.
The calculation model of this paper is a circle, so formula (9) is transformed into polar coordinates:
2.2.2. Gas Lubrication Control Equation
The general form of the Reynolds equation is shown in where is gas film clearance, is density, and is the pressure.
The ideal gas state equation is shown in
is standard atmospheric pressure, and is standard gas density.
In low speed state, there is negligible relative sliding speed. Because the static gas bearing is calculated in this paper, the pressure does not change with time, and the timeterm is equal to 0. Therefore, the right side of (13) is equal to 0. Because of the existence of EEPG, the gas film clearance can not be eliminated from the Reynolds equation, and (13) is simplified and transformed it into polar coordinates:
The mass flow rate of the bearing is equal to the mass flow of the outlet bearing.
Here,where is mass flow in, is mass flow out, is acceleration of gravity, is dynamic viscosity, is outlet pressure, is supply pressure, is temperature, is gas constant, is adiabatic index, is flow coefficient is orifice diameter, and is flow velocity coefficient.
Given the boundary conditions and the simultaneous equations (14) and (15) to obtain the numerical solution, we can get the pressure distribution of the gas film surface and calculate the bearing capacity .
The calculation formula of bearing stiffness is as follows:where is stiffness, is the change volume of gas film clearance, and is the change volume of bearing capacity.
The calculation formula of stiffness is obtained:
2.3. The Boundary Conditions
The EEPG generates deflection deformation in the range of , when or , and no elastic deformation, .
The pressure at the orifice of the gas lubrication governing equation, the pressure at the edge of the bearing .
2.4. Difference Scheme
The governing equations (10) and the governing equations (14) are discretized separately. This paper uses 2order accuracy central difference and the control equation (10) of the 3order and 4order difference schemes such as formula (20). The governing equation (14) uses 1order and 2order difference scheme, which is omitted in this paper. The 4order finite difference scheme in polar coordinates is shown in Figure 3.where are grid number in the direction of , respectively; are grid distance in the direction of , respectively.
2.5. GasSolid Coupling Calculation
The calculation process of this paper is the coupling process of the deformation control equation and the gas lubrication governing equation of EEPG. The flow chart of gassolid coupling is shown in Figure 4. Firstly, the calculation parameters are initialized, and the boundary conditions and coupling precision are set, according to the EEPG deformation control equation to calculate the initial deformation of EEPG. Then the pressure distribution of the gas film surface is solved by the gas lubrication control equation. In the coupling process, the coupling condition is the difference between the pressure distribution of the gas film surface and processing pressure, to solve the deformation control equation of EEPG, get the deformation of elastic thin plate after the completion of the first coupling. Then, solve the gas lubrication control equation, and the numerical values of the pressure distribution of the film surface is obtained. The EEPG is also obtained, marked as couplingh1. Set the results of the first coupling output as the initial conditions, and then solve elastic pressure groove deformation control equation and gas lubrication control equation again, the numerical values of the pressure distribution of the film surface and EEPG are obtained. The value of EEPG is couplingh2; the difference between the value of the EEPG of the two coupling outputs is used to judge the convergence condition . If , coupling ends; if , the calculation returns to the beginning of the first coupling. Through several iterations, until the error is less than the coupling accuracy, the coupling process is completed, and the numerical results of the EEPG after the gassolid coupling are given. Finally, the gas lubrication control equation is solved, and the pressure distribution of gas film surface is obtained; calculate bearing capacity and stiffness.
3. Test Rig
Build the aerostatic bearing stiffness test rig, as shown in Figure 5. The experimental bench is loaded by the jack, and the load is transmitted to the tested bearing through the force sensor and the force transmission rod, and the working gas film clearance of the bearing is measured by the microdisplacement sensor. The maximum output pressure of the test bench is 5000 N, and the measurement accuracy of the microdisplacement sensor is 0.1 μm.
4. Results and Discussion
4.1. Theoretical Result
According to the calculation method of the previous section, the calculation parameters of gas bearing, as shown in Table 1, are selected to calculate the gassolid coupling. The deformation curves of the EEPG and the film pressure distribution of the aerostatic bearing are obtained under different gas film clearance, as shown in Figure 6, where (i) is the elastic pressure groove coupling deformation graph and (ii) is the gas film pressure distribution. Suppose that this calculation is a computational model 1.

(a) μm
(b) μm
(c) μm
(d) μm
In order to compare the performance of the aerostatic bearing with EEPG, this paper calculated the performance of aerostatic bearing with REPG in the same performance parameters. Calculation parameters are also shown in Table 1, and the depth value of the pressure equalizing groove does not change with the pressure. At this time, we can obtain the film pressure distribution of the aerostatic bearing with REPG by solving (12), as shown in Figure 7. Suppose that this calculation is a computational model 2.
(a) μm
(b) μm
(c) μm
(d) μm
According to the calculation results of model 1 and model 2, the numerical performance curves of the 2 kinds of aerostatic bearings can be obtained, as shown in Figure 8, where model 1 is the aerostatic bearing with EEPG and model 2 is the aerostatic bearing with REPG. Figure 8(a) is the bearing capacity comparison curve, Figure 8(b) is the stiffness contrast curve, and Figure 8(c) is a gas flow rate curve.
(a) Bearing capacity curve
(b) Stiffness curve
(c) Mass flow curve
It can be seen from Figure 8(a) that the bearing capacity decreases with the increase of gas film clearance. The bearing capacity of model 1 overall compared to model 2 bearing capacity is small, but when the gas film clearance is less than or equal to 5 μm ( μm), model 1 and model 2 have the same bearing capacity. This is because when the gas film clearance is small, the pressure on the two sides of the elastic thin plate of model 1 are nearly equal, and the pressure difference is approximately equal to 0. Therefore, the elastic pressure groove has almost no deformation, as shown in Figure 6(a)(i), and the gas film pressure distribution is shown in Figure 6(a)(ii).
When the gas film clearance is from 6 μm to 15 μm (15 μm > > 6 μm), the bearing capacity of model 1 is less than that of model 2. Because the air supply remained unchanged, with the gas film clearance increases, the bearing capacity is decreased. Thus, the pressure difference between the upper and the lower surfaces of the elastic thin plate is increased, and the EEPG is forced to deform outwards, the sealing ability of the EEPG decreases, and the bearing capacity decreases rapidly. The deformation of the EEPG is shown in Figures 6(b)(i) and 6(c)(i). The deformation of the EEPG is uneven, and the maximum deformation occurs between the 2 orifices. The gas film pressure distribution is shown in Figures 6(b)(ii) and 6(c)(ii).
When the gas film clearance is more than 16 μm ( μm), model 1 bearing capacity is much smaller than model 2. Because the pressure difference of the elastic thin plate is greater than the initial working pressure, the deformation of EEPG is greater than the initial EEPG deformation, so the EEPG sealing effect disappeared, and bearing force decreased rapidly. The deformation of the EEPG is shown in Figure 6(d)(i), and the gas film pressure distribution is shown in Figure 6(d)(ii) ( μm).
As shown in Figure 8(b), the stiffness of model 1 and model 2 increases first and then decreases with the increase of gas film clearance. The maximum stiffness of model 1 is in the film gap μm, and the maximum stiffness of model 1 is 139. The maximum stiffness of model 2 is in the film gap μm, and the maximum stiffness of model 1 is 87. Model 1 stiffness increases compared to model 2 by 59.7%, but when the gas film clearance is greater than 9 μm, the stiffness of model 2 is greater than that of model 1.
It can be seen from Figure 8(c) that the mass flow rate of the bearing increases with the increase of the gas film clearance. When the gas film clearance is μm, the mass flow of model 1 is much less than that of model 2. This is the result of the gassolid coupling of EEPG.
According to the above analysis, it is found that the performance of aerostatic thrust bearing with EEPG is better than that of the aerostatic bearing with REPG.
4.2. Experimental Result
Figure 9 is the experimental test bearing, processing parameters as shown in Table 1, the bearing diameter is equal to 60 mm, and orifice diameter is equal to 0.3 mm. The experimental bearings were treated by grinding and polishing, and the surface roughness was less than 0.2 μm, but the region of EEPG was not included.
The experimental results are shown in Figure 10. The results of numerical calculation and experimental test results are basically the same. The gas film clearance increases, and the bearing capacity decreased rapidly. The experimental results show that the maximum stiffness of the bearing with EEPG is obtained at μm and is consistent with the numerical calculation and well validated the numerical gassolid coupling calculation is correct. However, the bearing capacity of the test bearing is less than the bearing capacity of the numerical calculation, because of the influence of the surface roughness of the test bearing, the position accuracy of the orifice, and the grinding precision of the EEPG.
(a) Bearing capacity curve
(b) Stiffness curve
5. Conclusion
In this paper, the bearing characteristics of aerostatic bearing with EEPG are studied by numerical calculation and experiment. Through the research, we can draw the following conclusions.(1)In this paper, the use of gassolid coupling adaptive mode of aerostatic bearing is feasible in theory and practice.(2)The stiffness of aerostatic thrust bearing with EEPG is increased by 40% than that of the aerostatic bearing with REPG.(3)The structure and location of the EEPG have a great influence on the performance of the bearing, so how to optimize the design warrants further research.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors would like to acknowledge that the project was supported by the National Natural Science Foundation of China (51375384), China Postdoctoral Fund (2016M602937XB), and Xi’an Technological University gas lubrication research and innovation team construction project.
References
 J. Du, G. Zhang, T. Liu, and S. To, “Improvement on load performance of externally pressurized gas journal bearings by opening pressureequalizing grooves,” Tribology International, vol. 73, no. 5, pp. 156–166, 2014. View at: Publisher Site  Google Scholar
 S. Gao, K. Cheng, S. Chen, H. Ding, and H. Fu, “CFD based investigation on influence of orifice chamber shapes for the design of aerostatic thrust bearings at ultrahigh speed spindles,” Tribology International, vol. 92, pp. 211–221, 2015. View at: Publisher Site  Google Scholar
 Y. F. Li, Y. H. Yin, H. Yang, X. N. Liu, J. Mo, and H. L. Cui, “Modeling for optimization of circular flat pad aerostatic bearing with a single central orificetype restrictor based on CFD simulation,” Tribology International, vol. 109, pp. 206–216, 2017. View at: Google Scholar
 Y. R. Jeng and S. H. Chang, “Comparison between the effects of singlepad and doublepad aerostatic bearings with pocketed orifices on bearing stiffness,” Tribology International, vol. 66, no. 7, pp. 12–18, 2013. View at: Publisher Site  Google Scholar
 J.A. Zhang, H. Dong, F.J. Ma, and B. Liu, “Performance analysis of the flotation cushion with elastic pressure equalizing groove for aerostatic slideway,” Applied Mechanics and Materials, vol. 152154, pp. 743–748, 2012. View at: Publisher Site  Google Scholar
 U. Nishio, K. Somaya, and S. Yoshimoto, “Numerical calculation and experimental verification of static and dynamic characteristics of aerostatic thrust bearings with small feedholes,” Tribology International, vol. 44, no. 12, pp. 1790–1795, 2011. View at: Publisher Site  Google Scholar
 A. P. Hao and Y. H. Jia, “Research on the characteristics of bearing behavior of the profiled throttling cavity air bearing,” in Proceedings of the 7th AsianPacific Conference on Aerospace Technology and Science, APCATS 2013, pp. 367–377, May 2013. View at: Publisher Site  Google Scholar
 A. Charki, K. Diop, S. Champmartin, and A. Ambari, “Numerical simulation and experimental study of thrust air bearings with multiple orifices,” International Journal of Mechanical Sciences, vol. 72, no. 1, pp. 28–38, 2013. View at: Publisher Site  Google Scholar
 G. Belforte, F. Colombo, T. Raparelli, A. Trivella, and V. Viktorov, “Comparison between grooved and plane aerostatic thrust bearings: Static performance,” Meccanica, vol. 46, no. 3, pp. 547–555, 2011. View at: Publisher Site  Google Scholar
 T.J. Zi, J.A. Zhang, Z.W. Lu, and B. Liu, “Numerical analysis of sewing air floating needle bar based on slit type guide rail,” Applied Mechanics and Materials, vol. 670671, pp. 866–872, 2014. View at: Publisher Site  Google Scholar
 H. Dong, X.L. Zhao, and J.A. Zhang, “Static characteristic analysis and experimental research of aerostatic thrust bearing with annular elastic uniform pressure plate,” Advances in Mechanical Engineering, vol. 7, no. 3, pp. 1–13, 2015. View at: Publisher Site  Google Scholar
 S. H. Chang, C. W. Chan, and Y. R. Jeng, “Numerical analysis of discharge coefficients in aerostatic bearings with orificetype restrictors,” Tribology International, vol. 90, pp. 157–163, 2015. View at: Publisher Site  Google Scholar
 T. Raparelli, V. Viktorov, F. Colombo, and L. Lentini, “Aerostatic thrust bearings active compensation: Critical review,” Precision Engineering, vol. 44, pp. 1–12, 2016. View at: Publisher Site  Google Scholar
 G. Aguirre, F. AlBender, and H. Van Brussel, “A multiphysics model for optimizing the design of active aerostatic thrust bearings,” Precision Engineering, vol. 34, no. 3, pp. 507–515, 2010. View at: Publisher Site  Google Scholar
 D. Ghodsiyeh, F. Colombo, T. Raparelli, A. Trivella, and V. Viktorov, “Diaphragm valvecontrolled air thrust bearing,” Tribology International, vol. 109, pp. 328–335, 2017. View at: Publisher Site  Google Scholar
 H. W. Liu, Theory of plate and shells, Zhejiang University Press, 1987.
 G. Belforte, T. Raparelli, V. Viktorov, and A. Trivella, “Discharge coefficients of orificetype restrictor for aerostatic bearings,” Tribology International, vol. 40, no. 3, pp. 512–521, 2007. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Zhao Xiaolong 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.