Next Article in Journal
Controllability of Semilinear Stochastic Generalized Systems in Hilbert Spaces by GE-Evolution Operator Method
Next Article in Special Issue
The Soliton Solutions of the Stochastic Shallow Water Wave Equations in the Sense of Beta-Derivative
Previous Article in Journal
Modern Dimensional Analysis-Based Heat Transfer Analysis: Normalized Heat Transfer Curves
Previous Article in Special Issue
Reductions and Exact Solutions of Nonlinear Wave-Type PDEs with Proportional and More Complex Delays
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Squeeze Film Effect with a High-Pressure Boundary in Aerostatic Bearings

Center for Precision Engineering, Harbin Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(3), 742; https://doi.org/10.3390/math11030742
Submission received: 20 December 2022 / Revised: 29 January 2023 / Accepted: 31 January 2023 / Published: 1 February 2023

Abstract

:
The squeeze film effect was discussed in several fields, but mostly under the same pressure boundary conditions. However, pressures at the inlet and outlet are different for aerostatic bearings. In this paper, the dynamic Reynolds equation group, with the stiffness and damping pressure written separately, is deducted and numerically solved with a high-pressure boundary for a parallel flat and circular thin film. The circular thin film considers the two results of the supply pressure boundary inside and outside. All dynamic pressure distribution and stiffness curves are given in a dimensionless form, and a comparative analysis of squeeze film characteristics with and without external pressure is conducted. From the calculation results, it can be concluded that the squeeze effect shows damping for zero-frequency and stiffness for infinite-frequency for compressible lubricants. The dynamic pressure in the static high pressure region is also high at high frequencies affected by gas compressibility. Based on these analytical results, the transfer functions of the thin film are given to further analyze the dynamic performance of aerostatic bearings, and the shape of the response curve approximates an exponential decay form, even when the amplitude increases to 10% of the gas film thickness.

1. Introduction

When two surfaces separated only by a thin air film have relative normal motion, the air film has a reaction force against surfaces to block movement tendency, which is the squeeze film effect. The squeeze film effect is reflected in many mechanical parts, with the common feature being that all of the boundary is at atmospheric pressure.
For example, the squeeze film damping effect is quite important at ball raceway contact in ball bearing; the damping of ball bearings is mainly derived from the squeeze film effect [1,2,3]. The squeezed film at ball raceway contact is an inherent characteristic of ball bearings and may also occur at the contact of the outer ring with a bearing seat. When the outer ring of the bearing is designed to slide freely and bear an additional load of the axial thermal expansion, there is a small gap in the outer ring of the bearing, and the squeeze film effect is more obvious. The air damping effect of micro-electromechanical systems (MEMSs) is also very obvious. Micro-electromechanical systems such as micro-accelerometers and micro-gyroscopes are limited by technology and packaging conditions, and their moving structures generally work in an environment with a certain gas pressure. As the characteristic size of a micro-structure shrinks, its surface area-to-volume ratio gradually increases, and the air damping effect becomes extremely significant [4]. In the capacitive micro-machined accelerometer described in [5], when the central plate moves up and down, the air caught in the gap between the plates is squeezed, then force is exerted on both plates, and the squeeze film effect can be significant, which provides effective damping in a wide-frequency band. Besides these two cases, a squeeze film before clutch contact affects its performance. A wet clutch soaks the disc in the transmission oil, adopts oil-immersed heat dissipation, and has a high heat dissipation efficiency. However, the oil medium delays the physical contact between the parallel disc clutch surfaces, and its performance is seriously affected [6,7]. The squeeze film effect also contributes to tire slipping on wet roads. When a car drives over a water-covered road, the lubricating effect of the water film significantly reduces the coefficient of the friction between a tire and road, and the tires are then prone to slipping. When the tire speed exceeds a certain limit value, the dynamic pressure generated by the water flow is sufficient to lift the tire off the ground. Once hydroplaning occurs, tires will completely lose their drive and braking forces [8]. Similarly, walking with rubber soles on wet or icy surfaces leads to similar problems [9]. Typical examples of making full use of the squeeze film effect are squeeze film dampers and vibro-tactile texture reproduction techniques. A squeeze film damper is installed between the rolling bearing and support structure, and an oil gap is formed between the outer rings of connected oil film. Squeeze oil film dampers can effectively absorb the vibration energy of rotors and ensure their smooth operation. It is one of the most important research achievements affecting high-speed rotor dynamics in recent decades. At present, squeeze film damping technology is widely used in high-speed rotating machinery such as aero-engines, turbo-compressors, and industrial gas engines [10]. Texture tactile reproduction is a specific hardware device that generates physical stimuli that act on users when in contact with a textured surface, thereby simulating the surface characteristics of objects such as roughness, hardness, and viscosity. Texture haptic technology is used to generate ultrasonic waves via a touch screen, and the squeeze film effect is used to create an air gap between the user’s finger and the display, thus reducing the friction and increasing smoothness [11]. It can be seen that the principle of this technique is identical to that of squeeze film bearings.
There is one common feature in the aforementioned occasions where the squeeze film effect exists: all of the boundary is at atmospheric pressure, or more precisely, the boundary pressure is the same. The problem with constant pressure boundary conditions has previously been studied [12,13], in particular by Langlois [14], Gross [15], Pan [16], and Blech [17]. However, for externally pressurized air bearings (EPABs), the pressure decreases from a high-pressure in the inlet to an atmosphere condition in the outlet, which possesses different characteristics compared with the same-pressure boundary. In the early days of air bearing development, discussions of stability were mainly aimed at understanding whether a bearing is affected by a pneumatic hammer or at assessing the stability margin [18,19]. There has been much interest in evaluating the damping capacity of bearings, which is quite critical because EPABs are famous for having relatively little damping. Analytical, numerical, and experimental methods were used to study the dynamic performance of EPABs. In addition to experiments, in general, there are two methods used to solve the dynamic characteristics of an EPAB: one is the time-linearization or perturbation method, and the other is the time-domain iterative method. The system is solved through the time domain response and the dynamic stiffness parameters are identified from it. The basic idea of the perturbation method is to superimpose a gap change caused by displacement or velocity in the steady-state gap as a small disturbance, and to approximate the transient Reynolds equation with a small amount, ignoring all high-order terms [20]. The dynamic characterization of EPABs includes an evaluation of stiffness and damping coefficients as a function of excitation frequency. A change in bearing capacity obtained by dynamic pressure integration divided by the disturbance value is the dynamic stiffness of the gas film. Licht first used the perturbation method to test the accuracy of lumped parameter models of single-bore thrust bearings [21]. Boffey analyzed the effect of additional rubber rings using Licht’s calculation method, noting that external damping can improve stability [22], which was experimentally verified [23]. The advent of high-performance calculation technology made it possible to build complex numerical models. Stiffler modeled the feedhole boundary between the central disk region and the exterior annular region as a line source, and the dynamic pressure distribution of rectangular and circular EPABs were investigated [24,25]. Miyatake and Yoshimoto used computational fluid dynamics (CFD) technology to numerically simulate the static and dynamic performance of a circular EPAB with small feed holes [26]. By introducing the concept of an infinitesimal radius, Bhat overcame the mathematical difficulty of the conformal transformation of disc thrust gas bearing center point by the finite element method (FEM), analyzed the dynamic characteristics of the capillary-restricted air bearing pad, and determined the dynamic stiffness with the disturbance frequency [27]. These distributed parameter models can be used to accurately analyze and solve for the dynamic characteristics of EPABs, but cannot provide a direct link between design parameters and the performance index. The influence of different parameters on the dynamic characteristics of bearings is coupled together, and it is difficult to form a systematic design method. To solve this problem, Colombo et al. proposed a simple lumped parameter (LP) model of infinite parallel plates [28], which is based on the discretization of the Reynolds equation and treats the parallel plates as a combination of gas tank and pneumatic restrictors. This simple LP model is able to correctly predict the damping capacity, unlike the analytical solution. This basic idea has been adopted to predict the dynamic performance of EPABs, which can be seen as an assembly of orifices and gas film [29,30,31]. Starting from this point, if we can determine the influence of various parameters of the bearing on its transfer function, then the dynamic performance of EPABs can be appropriately designed, according to the zero-pole configuration of the control system.
It should be noted that the boundary conditions for normal orifice-type aerostatic bearings dynamically change, since an orifice not only generates high pressure but also controls the mass flowing into the gas film surface, which is a major difference from the constant pressure boundary found in this paper. In fact, the dynamic characteristics of the gas film solved in this paper are a simplification of a special kind of orifice-type bearing known as surface-restriction aerostatic bearings, whose main feature is that the boundary pressure is different, but the external pressure is constant. There are very few studies on this type of gas bearing, and there is no research on its dynamic characteristics. In this paper, the squeeze film effect of thin films inside EPABs with a high-pressure boundary is determined, and compared with that of same-boundary conditions. For both boundary types, this paper offers an important view on the dynamic performance of air bearings: air bearings behave according to the squeeze number. For low squeeze numbers, when the exciting frequency is low or the film is thick, the viscous forces dominate, and the air is sucked in or forced out of the film when the air gap height oscillates. If there is no orifice, the air gap should be treated as a viscous damper. For high squeeze numbers with an increase in exciting frequency or with small air gaps, the compressibility effects become stronger, which makes the air gap behave like a spring. The key to analyzing the squeeze film problem is to determine the time-dependent pressure distribution inside the film caused by the relative motion. When gas compressibility is considered, the results of the squeezed film solution are complicated, and its properties cannot be intuitively seen, even for a simple disk boundary. Therefore, this paper aimed to calculate the dynamic pressure distribution for common components of EAPBs, i.e., infinite long thin films and circular films considering both inside flow and outside flow.

2. The Compressible Reynolds Equation and Its Perturbation Form

2.1. Reynolds Equation

When the bearing speed is small or the bearing is not moving, the aerodynamic effect can be neglected, and the Reynolds equation considering the ideal gas compressibility can be expressed as [32]:
· h 3 12 η P abs P abs = ( h P abs ) t
Here, is the Laplacian operator, h is the thin film thickness, η is the dynamic viscosity of the lubricants, and P abs is the absolute pressure inside a thin film.
The expansion of Equation (1) in the Cartesian coordinate system is [32]:
x h 3 12 η P abs P abs x + y h 3 12 η P abs P abs y = ( h P abs ) t
The expansion of Equation (1) in the cylindrical coordinate system is [32]:
1 r r h 3 12 η r P abs P abs r + 1 r θ h 3 12 η 1 r P abs P abs θ = ( h P abs ) t
When solving the static characteristics of EPABs, the term on the right side of Equation (1) is zero. Using the differential equation f d f = d f 2 / 2 , it can be expressed with the pressure square P abs 2 as the variable [33]:
x h 3 P abs 2 x + y h 3 P abs 2 y = 0

2.2. The Perturbation Form of the Reynolds Equation

The general problem of fluid lubricated bearings can be analyzed in terms of the three following aspects [34]: (a) steady-state conditions; (b) transient operating conditions; and (c) effects due to small-amplitude harmonic changes in geometry. The steady-state conditions are the starting point of all calculations, describing the conditions for a balanced operation of the entire system. This step requires the calculation of the bearing capacity, friction loss, mass flow rate, etc. The transient operating conditions are considered for the time-transient solution of a coupled dynamical system, and the results can be used to obtain dynamic stability or study the nonlinear behavior in the case of large offsets.
The position of an air bearing is not fixed. In the process of cutting or under micro-shocking, the bearing vibrates regularly around the equilibrium position. According to Newton’s second law, the dynamic displacement of the air film can be calculated by the following formula:
h 1 = v 0 t + W 1 ( h 1 ( t ) ) + external force M d t 2 , W 1 = P abs 1 d Ω
The direct calculation of the above formula not only includes a the space integral of dynamic pressure distribution, but also includes the time integral of the position change, which relatively complicates the solution process. The effects due to small-amplitude harmonic changes in the geometry are considered to determine the motion of the bearing in the normal direction of the air film and only determines the dynamic behavior linearized near a given operating point; that is, the response process of the bearing after a small impact is analyzed, and the perturbation method is used to handle the time term in the Reynolds equation [20]. Note that variables such as the pressure and clearance in the gas film are linearly superimposed by the steady-state item and the time-varying perturbation item, that is:
P abs = P abs 0 + P abs 1 , h = h 0 + h 1
The subscript 0 indicates the steady-state item, which is completely unrelated to time t. The subscript 1 indicates the dynamic item, and the amplitude of the dynamic item is much smaller than the steady-state item. The perturbation method linearizes the dynamic behavior of the air film near the steady-state operating point, and the superposition of the input leads to the superposition of the output. Therefore, various tools in the linear system can be used to analyze the characteristics of the system. The gas film is equivalent to the superposition of a spring and damper, and the stiffness and damping characteristics of the gas film are solved first.
By bringing the expression of the variable in Equation (6) into Equation (2) and omitting the high-order small amount, the steady-state Reynolds Equation (4) and perturbation equations can be obtained, respectively [35]:
2 · h 0 3 P abs 0 P abs 1 + 3 · h 0 2 h 1 P abs 0 2 = 24 η h 0 P abs 1 t + P abs 0 h 1 t
We can transpose and rearrange the equations as follows:
h 0 P abs 1 t + 1 12 η · h 0 3 P abs 0 P abs 1 = P abs 0 h 1 t 1 8 η · h 0 2 h 1 P abs 0 2
The above equations show that the dynamic pressure and its rate of change are coupled with the air film gap and its rate of change (the squeeze motion speed). The thrust of the air film on the air-floating surface is obtained by the pressure integral, so the spatial distribution of the corresponding air film pressure cannot be obtained by only specifying the dynamic force; that is, the characteristics of the air film can only be analyzed via the flexibility method: firstly, we specify the displacement form on the air film, and then solve the pressure change. Suppose that h 1 changes periodically in the following form:
h 1 = Re ( h ˜ 1 e j ω t ) = Re ( h 1 c + j h 1 s ) e j ω t = h 1 c cos ( ω t ) h 1 s sin ( ω t )
d h 1 d t = Re ( j ω h ˜ 1 e j ω t ) = ω [ h 1 c sin ( ω t ) + h 1 s cos ( ω t ) ]
For notation convenience, the real part is omitted, and a single-frequency vibration is directly expressed in the form of a complex exponential. The wavy line indicates a complex variable. The pressure distribution P abs 1 ( x , y , t ) produces a corresponding periodic variation:
P abs 1 = Re ( P ˜ 1 e j ω t ) = Re ( P 1 c + j P 1 s ) e j ω t = P 1 c cos ( ω t ) P 1 s sin ( ω t )
P abs 1 t = Re ( j ω P ˜ 1 e j ω t ) = ω [ P 1 c sin ( ω t ) + P 1 s cos ( ω t ) ]
By bringing the complex expressions of displacement and pressure into Equation (7), the complex form of the perturbation Reynolds equation is obtained [35]:
2 · h 0 3 P abs 0 P ˜ 1 + 3 · h 0 2 h ˜ 1 P abs 0 2 = 24 η j ω h 0 P ˜ 1 + j ω P abs 0 h ˜ 1
In fact, the above equation can be obtained by performing Laplace transformation on Equation (7) and letting s = j ω . By using the flexibility method to determine the characteristics of the gas film, the form of displacement can arbitrarily specified. Therefore, let h 1 c = 0 , and we have h 1 = Re j h 1 s e j ω t = h 1 s sin ( ω t ) , P abs 1 = Re ( P 1 c + j P 1 s ) e j ω t , Thereby,
2 · h 0 3 P abs 0 P 1 s 24 η ω h 0 P 1 c = 3 · h 0 2 h 1 s P abs 0 2
12 η ω h 0 P 1 s + · h 0 3 P abs 0 P 1 c = 12 η ω P abs 0 h 1 s
The above equations and Equation (13) indicate that the real part and imaginary part of the dynamic pressure are mutually coupled; the amplitude of dynamic pressure P abs 1 ( x , y , t ) is proportional to the forced displacement amplitude h 1 . However, the phase of dynamic pressure has a nonlinear relationship with forced displacement frequency ω .

2.3. Stiffness and Damping of the Gas Film

Regardless of the influence of the fluid–solid coupling effect, the sides of the thin film are treated as a single rigid body, and a transfer function is used to describe the dynamic characteristics of the bearing. In this kind of model, the plates composing lubricating film move as one. When the bearing is in translation, h 1 x = h 1 y = 0 ; therefore, h 1 in Equation (7) can be written outside the brackets, namely
2 · h 0 3 P abs 0 P abs 1 + 3 h 1 · h 0 2 P abs 0 2 = 24 η h 0 P abs 1 t + P abs 0 d h 1 d t
Equation (13) then changes to
2 · h 0 3 P abs 0 P ˜ 1 + 3 h ˜ 1 · h 0 2 P abs 0 2 = 24 η j ω h 0 P ˜ 1 + j ω P abs 0 h ˜ 1
and Equation (14) changes to
2 · h 0 3 P abs 0 P 1 s 24 η ω h 0 P 1 c = 3 · h 0 2 P abs 0 2 h 1 s
12 η ω h 0 P 1 s + · h 0 3 P abs 0 P 1 c = 12 η ω P abs 0 h 1 s
The dynamic stiffness of the gas film is defined as follows (note that the positive and negative of the calculation formula are related to the definition of the displacement direction, a sign of the real part and imaginary part of the dynamic pressure):
K b ( ω ) = K ( ω ) + j ω C ( ω ) = W ˜ 1 d Ω h ˜ 1 = Ω P ˜ 1 d Ω h ˜ 1 = Ω P 1 c d Ω · h 1 c + Ω P 1 s d Ω · h 1 s h 1 c 2 + h 1 s 2 + j ω Ω P 1 s d Ω · h 1 c Ω P 1 c d Ω · h 1 s ω ( h 1 c 2 + h 1 s 2 )
Obviously, the calculation of dynamic stiffness from the above formula is complicated. When forcing a periodically varying displacement is h 1 s sin ( ω t ) by making h 1 c = 0 , the formulas for calculating the dynamic stiffness and damping of the air film are as follows:
K ( ω ) = Ω P 1 s d Ω h 1 s , C ( ω ) = Ω P 1 c d Ω ω h 1 s

3. Analytical Results with Same Pressure Boundaries

3.1. A One-Dimensional Infinite Width Flat Air Film

Figure 1a shows a one-dimensional infinitely wide flat lubricating film. The gas film is surrounded by ambient pressure, which means that only two flat plates that are close to each other form a thin film in the air, and there is no flow when the upper plate does not move.
When fluid compressibility is considered, the dynamic Reynolds equation of the one-dimensional lubricating film in the Cartesian coordinate system is obtained from Equation (2):
h 3 d 2 P abs 2 d x 2 = 24 η d ( h P abs ) d t
By substituting the perturbation form in Equation (6) and the complex form in Equation (9), Equation (17) under the condition of a one-dimensional flat gas film is obtained:
d 2 ( P abs 0 P ˜ 1 ) d x 2 + 3 2 h ˜ 1 h 0 d 2 ( P abs 0 2 ) d x 2 = 12 η j ω h 0 2 P ˜ 1 + 12 η j ω P abs 0 h 0 3 h ˜ 1
When there is no high-pressure gas supply at the boundary, that is, the atmospheric pressure boundary, the steady-state pressure distribution is P abs 0 ( x ) = P a , which further simplifies Equation (23) with boundary conditions:
d 2 ( P ˜ 1 ) d x 2 = 12 η j ω P a h 0 2 P ˜ 1 + 12 η j ω h 0 3 h ˜ 1 , P ˜ 1 ( 0 ) = 0 , P ˜ 1 ( l ) = 0 , d P ˜ 1 d x l 2 = 0
The analytical solution to this differential equation is:
P ˜ 1 ( x ) = P a h ˜ 1 h 0 cosh [ σ s e j π / 4 ( x l 1 2 ) ] cosh ( σ s e j π / 4 1 2 ) 1 , σ s = 12 η ω l 2 P a h 0 2
where σ s is the squeeze number. The expanded expression is as follows:
P 1 s ( x * ) = 2 cosh k + cos k cosh k 2 cos k 2 cosh k x * cos k x * + sinh k 2 sin k 2 sinh k x * sin k x * 1 P 1 c ( x * ) = 2 cosh k + cos k sinh k 2 sin k 2 cosh k x * cos k x * cosh k 2 cos k 2 sinh k x * sin k x * and , x * = x l 1 2 , k = σ s 2
Figure 2 and Figure 3 show the distribution of the dynamic pressure of a one-dimensional flat lubricating film under different compression numbers (no external air supply; only an atmospheric pressure boundary). Figure 2 shows a schematic diagram of the distribution of the pressure P 1 s in the same phase with external excitation, which produces dynamic stiffness. Figure 3 shows a schematic diagram of the 90-degree-phase pressure distribution P 1 c , which produces damping. Since the boundary condition is a constant pressure, the dynamic pressure at the boundary is zero; thus, intuitively, the compressibility of the gas causes some of it to be trapped in the central area, and it cannot be removed and inhaled in time.
The definition of the dynamic stiffness Equation (20) can be used to obtain the dimensionless dynamic stiffness and damping expressions of the air film:
K * ( ω ) = K ( ω ) b l P a / h 0 = Re 1 2 tanh σ s e j π / 4 / 2 σ s e j π / 4 = 1 1 k sinh k + sin k cosh k + cos k
C * ( ω ) = C ( ω ) η b l 3 / h 0 3 = 12 σ s Im 1 2 tanh σ s e j π / 4 / 2 σ s e j π / 4 = 12 σ s 1 k sinh k sin k cosh k + cos k
Figure 4 shows the variation in the dimensionless dynamic stiffness and the damping coefficient with the squeeze number. It is known that the dynamic coefficients of the incompressible fluid lubricating film are all constant, the stiffness is zero, and the damping is constant. As a comparison, the dynamic coefficient of the compressible gas film varies with the excitation frequency: when the excitation frequency is ω 0 , the squeeze number σ s 0 , and the gas film behaves like the incompressible lubricating film; with the excitation frequency increasing, the squeeze number increases, and when σ s > 100 , the air film only shows rigidity, and the damping is almost negligible.
As shown in Figure 5, the dynamic stiffness amplitude varies with the squeeze number: as the squeeze number increases, the dynamic coefficient of the incompressible fluid lubricating film increases linearly, but when the squeeze number σ s > 10 , the dynamic stiffness amplitude of the air film no longer increases, and the dimensionless number approaches the limit 1. With the increase in the squeeze number, the rigid restoring force of the air film increases, but the damping force decreases. The frequency with the same amplitude is called the cutoff frequency. For a one-dimensional flat plate on the atmospheric boundary, the squeeze number corresponding to the crossover frequency is approximately 10.
The above analysis shows that whether the squeeze film exhibits damping or rigidity depends on the size of the squeeze number, and the squeeze number, i.e., σ s = 674.5 based on Table 1, is calculated according to Equation (25). It can be seen that, since the gap of the bearing is small, the stiffness effect produced by the squeeze film is relatively large. When analyzing the dynamic characteristics of the bearing, the squeeze film effect must be considered.

3.2. One-Dimensional Circular Air Film

As shown in Figure 1b, the normal movement of a gas film is superimposed on the basis of a circular plate gas film. When fluid compressibility is considered, the dynamic Reynolds equation of the one-dimensional lubricating film in the cylindrical coordinate system is obtained from Equation (3):
1 r r h 3 12 η r P abs P abs r = ( h P abs ) t
By substituting the perturbation form in Equation (6) and the complex form in Equation (9), the form of Equation (17) inside of a one-dimensional annular flat plate gas film is:
d d r d ( P abs 0 P ˜ 1 ) d r + 3 2 h ˜ 1 h 0 d P abs 0 2 d r + 1 r d ( P abs 0 P ˜ 1 ) d r + 3 2 h ˜ 1 h 0 d P abs 0 2 d r = 12 η j ω h 0 2 P ˜ 1 + 12 η j ω P abs 0 h 0 3 h ˜ 1
The solution without supply pressure (the atmospheric pressure boundary) will now be discussed. Under this condition, the steady-state pressure distribution is P abs 0 ( r ) = P a , That is, Equation (30) is further simplified, and there are boundary conditions (the inner hole radius is equal to zero):
d 2 P ˜ 1 d r 2 + 1 r d P ˜ 1 d r = 12 η P a h 0 2 j ω P ˜ 1 + 12 η h 0 3 j ω h ˜ 1 , P ˜ 1 ( R 2 ) = 0 , d P ˜ 1 d r 0 = 0
The analytical solution to this differential equation is:
P ˜ 1 ( r ) = P a h ˜ 1 h 0 J 0 ( σ s e j π / 4 r * ) J 0 ( σ s e j π / 4 ) 1 , σ s = 12 η ω R 2 2 P a h 0 2 , r * = r R 2
The expanded expression is as follows:
P 1 s ( r * ) = ber σ s ber ( σ s r * ) + bei σ s bei ( σ s r * ) ( ber σ s ) 2 + ( bei σ s ) 2 1 P 1 c ( r * ) = bei σ s ber ( σ s r * ) ber σ s bei ( σ s r * ) ( ber σ s ) 2 + ( bei σ s ) 2
In applied mathematics, the Kelvin functions ber ν ( x ) and bei ν ( x ) are the real and imaginary parts, respectively, of J ν x e 3 π i 4 , where x is real, and J ν ( z ) is the ν th order Bessel function of the first kind.
The definition of the dynamic stiffness of the air film, as given in Equation (20), can be used to obtain the dynamic stiffness and damping expressions of the air film:
K * ( ω ) = K ( ω ) π R 2 2 P a / h 0 = Re 1 2 J 1 ( σ s e j π / 4 ) σ s e j π / 4 J 0 ( σ s e j π / 4 ) = 1 2 σ s ber σ s bei 1 σ s ber 1 σ s bei σ s ber 1 σ s + bei 1 σ s ( ber σ s ) 2 + ( bei σ s ) 2
C * ( ω ) = C ( ω ) η π R 2 4 / h 0 3 = 12 σ s Im 2 J 1 ( σ s e j π / 4 ) σ s e j π / 4 J 0 ( σ s e j π / 4 ) 1 = 2 σ s ber σ s ber 1 σ s + bei 1 σ s bei σ s bei 1 σ s ber 1 σ s ( ber σ s ) 2 + ( bei σ s ) 2
Figure 6 shows the changing trend of the dimensionless dynamic stiffness and damping coefficient with the squeeze number, and Figure 7 shows the dynamic stiffness amplitude with the squeeze number. Obviously, the squeeze film characteristics of the parallel circular plate and flat lubricating film are similar, but the dimensionless damping at zero frequency becomes 1.5 (this is consistent with the damping obtained when compressibility is not considered), and then the crossover frequency shifts to the left. The squeeze number corresponding to the crossover frequency is approximately equal to σ sc 6 .
Similarly to a guide railway, by substituting Equation (32), the squeeze number of the common thrust bearing air film shown in Table 2 is estimated σ s = 168.6 . Since the clearance of the bearing is small, the stiffness effect of the squeeze film is relatively large. The result of the above formula shows that the squeeze film effect must be considered when analyzing the dynamic characteristics of the thrust bearing.

4. Numerical Results with Different Pressure Boundaries

The static and perturbation Reynolds equations are all elliptic differential equations without time, so they can be conveniently solved by the FEM. With the development of modern computing technology, the related literature has increased, so the basic theory of the FEM and the Reynolds equation solving process adopting FEM will not be expanded upon in detail here. The interested reader can refer to [36,37].
The static characteristics of the aerostatic bearing and the dynamic characteristics are evaluated considering the stiffness and damping coefficients. The flowchart for solving the Reynolds equation is shown in Figure 8. Firstly, the steady-state Reynolds equation is solved and the static performance is determined. Secondly, the coupled perturbation Reynolds equations are solved to obtain the dynamic performances of the air bearing.

4.1. A One-Dimensional Infinite Width Flat Air Film

The first example of a numerical solution to the perturbed Reynolds equation is as follows: The dynamic characteristics in the flat plate with external pressure equal gap as shown in Figure 9a are calculated. For the one-dimensional flat lubricating film shown in Figure 9a, when compressibility is considered, the steady-state pressure distribution is determined by Equation (2), and the dynamic pressure distribution is determined by Equation (23). In order to make the results more universal, the following variable substitutions are used:
P abs 0 = P abs 0 * P a , P d = P a / σ , P ˜ 1 = P ˜ 1 * P a , h ˜ 1 = h ˜ 1 * h 0 , x = x * l , σ s = 12 η ω l 2 P a h 0 2
where σ is the reciprocal of dimensionless air supply pressure. The dimensionless Equation (23) is as follows:
d 2 ( P abs 0 * P ˜ 1 * ) d x * 2 + 3 2 h ˜ 1 * d 2 ( P abs 0 * 2 ) d x * 2 = j σ s P ˜ 1 * + j σ s P abs 0 * h ˜ 1 * , P abs 0 * = P abs 0 * ( σ , x * )
The above dimensionless form shows that we only need to solve for different supply pressure ratios and squeeze numbers to fully understand the dynamic characteristics of one-dimensional flat films. The FEM was used to solve the above equations. Figure 10 and Figure 11 are the real part P 1 s and imaginary part P 1 c of the dynamic pressure distribution with different air supply pressures σ and different squeeze numbers σ s . It can be seen from those figures that, as the gas supply pressure increases, the high pressure area of the dynamic pressure is also shifted to the left, where more gas is compressed and concentrated, since the high-pressure area in the steady-state pressure distribution is near the origin of the x axis.
Integrating the real and imaginary parts of the dynamic pressure, Figure 12 shows the dimensionless dynamic stiffness and damping coefficient of the one-dimensional flat plate with a gas supply boundary, and its definition is the same as that under the atmospheric pressure gas supply boundary in Equation (27). Because the air supply pressure is greater than the atmospheric pressure, the value of the dimensionless dynamic stiffness is greater than 1. The maximum value of stiffness can be seen to increase from 1 to 6 when the supplying pressure changes from zero to 9 bar. As a comparison, the dimensionless dynamic stiffness in Figure 4 is always less than 1.
Improving the supply pressure does not improve the amplitude of the damping coefficient, but delays its decreasing frequency. It can be seen in Figure 12 that the half-damping squeeze number is increasing from 10 to 70 when the supplying pressure changes from 0 to 9 bar. It is somewhat strange that the value of stiffness is decreasing with the improvement in the supply pressure in the low-frequency region. The effects of the increasing damping and the decreasing stiffness together increase the cut-off frequency from 10 to 70 with an increase in supply pressure, as shown in Figure 13. Figure 13 clearly shows that the increase in the air supply pressure increases both the damping force and the stiffness force.

4.2. A One-Dimensional Annular Air Film Flowing from the Inside and the Outside

The second example of a numerical solution to the disturbed Reynolds equation is as follows: the one-dimensional flow model in cylindrical coordinate system shown in Figure 9b is calculated. Gas will flow from the inside to the outside if the inner diameter R 0 is of a high-pressure boundary and the outer diameter R 2 is of an atmospheric pressure boundary. In the same way, gas will flow from the outside to the inside if the outer diameter R 2 is of a supplying boundary. When considering gas compressibility, the steady-state pressure distribution of the one-dimensional annular lubricating film is determined by Equation (3), and the dynamic pressure distribution is determined by Equation (30). The following variable substitutions are used:
P abs 0 = P abs 0 * P a , P d = P a / σ , P ˜ 1 = P ˜ 1 * P a , h ˜ 1 = h ˜ 1 * h 0 , r = r * R 0 , L 2 = R 2 / R 0 , σ s = 12 η ω ( R 2 R 0 ) 2 P a h 0 2
where σ is the reciprocal of the dimensionless air supply pressure. The dimensionless Equation (30) is as follows:
d d r * + 1 r * d ( P abs 0 * P ˜ 1 * ) d r * + 3 2 h ˜ 1 * d P abs 0 * 2 d r * = j σ s P ˜ 1 * + j σ s P abs 0 * h ˜ 1 * , P abs 0 * = P abs 0 * ( σ , L 2 , r * )
The above dimensionless form of P abs 0 * indicates that we need to analyze the dynamic characteristics of a one-dimensional ring squeeze film with different external and internal diameter ratios, supply pressure ratios, and squeeze numbers. Using the FEM to solve the above equations, the real part of the dynamic pressure P 1 s and the imaginary part P 1 c with different air supply pressures σ and different squeeze numbers σ s (the ratio of outer and inner diameters L 2 = 1.5 ) are illustrated in Figure 14 and Figure 15, respectively. These figures show that more gas is compressed and concentrated close to the inner diameter with the increase in gas supply pressure, since the high-pressure area is located on the inner edge. The high-pressure area of the dynamic pressure also shifts to the inner diameter, which is similar to that of a one-dimensional flat plate.
The real part and imaginary part of the dynamic pressure are integrated to obtain the dimensionless dynamic stiffness and damping coefficient of the ring-type thin film, as shown in Figure 16, and its definition is similar to the definition under the atmospheric pressure gas supply boundary:
K * ( ω ) = K ( ω ) π ( R 2 2 R 0 2 ) P a / h 0 , C * ( ω ) = C ( ω ) η π ( R 2 2 R 0 2 ) ( R 2 R 0 ) 2 / h 0 3
Figure 16 shows the dynamic characteristic coefficient under different ratios of outer and inner diameters. It is noted that the dimensionless stiffness and damping ordinates in the figures are exactly the same. With the increase in the ratio of outer and inner diameters, the maximum stiffness slightly decreases, and the damping coefficient remains almost unchanged, which shows that the dimensionless value is basically irrelevant to the ratio of inner and outer diameters, and indicates that the aforementioned dimensionless definition is appropriate.
By carefully analyzing the definition of the dimensionless damping coefficient in Equation (40), it is clear that, due to the roughly quadratic relationship between the reference value and the ratio of the outer and inner diameters R 2 / R 0 , the radius ratio is key in the design of thrust bearings.
The solid line in Figure 16 represents the result of the flow from inside to outside, and the dotted line represents the performance curve of the flow from outside to inside. As the disturbance frequency increases, the dimensionless stiffness of the flow from the outside is larger than that from the inside. The definition in Equation (44) shows that, under the same gap, the larger the area of the high pressure zone, the larger the load value obtained by integration, and the larger the value of the maximum stiffness. Obviously, the high-pressure boundary at the outer edge R 2 is larger than the high-pressure area at the inner edge R 0 . Therefore, the flow direction has little effect on the damping coefficient, but as the radius ratio increases, the flow direction has a greater impact on the stiffness coefficient.
Because the air supply pressure is greater than the atmospheric pressure, the value of the dimensionless dynamic stiffness is greater than 1. It is clear that the maximum value of the stiffness increases from 1 to 5 when the supplying pressure changes from 0 to 6 bar. As a comparison, the dimensionless dynamic stiffness in Figure 6 is always less than 1.
Improving the supply pressure does not improve the amplitude of the damping coefficient, but it does delay its decreasing frequency. In the intermediate frequency region, the air supply pressure has a very strong influence on the damping (note that the abscissa is a logarithmic coordinate); it can be seen in Figure 16 that the half-damping squeeze number is increasing from 10 to 70 when the supplying pressure changes from 0 to 9 bar. It is slightly strange that the value of the stiffness is decreasing with the improvement in the supply pressure in the low-frequency region. The effects of the increasing damping and decreasing stiffness together increase the cut-off frequency from 10 to 60 as the supply pressure increases, as shown in Figure 16. Figure 17 clearly shows that the increase in the air supply pressure increases both the damping force and the stiffness force.
Figure 18 shows the dynamic stiffness modulus of the 1D ring with supply pressure under different radius ratios when the air flows in from outside. It can be seen that the flow direction does not cause much difference compared with Figure 16 when the air flows from inside.

5. Results and Discussion

5.1. Characteristics of the Thin Film with Exciting Frequency Approaching Zero

Although the calculation results in Equation (16) show that the real and imaginary parts of the dynamic pressure are coupled, some basic conclusions can be drawn at low and high frequencies through limit analysis.
When the exciting frequency approaches zero ω 0 , Equation (16) can be simplified:
· h 0 3 P abs 0 P 1 s + 3 2 h 1 s · h 0 2 P abs 0 2 = 0 , P 1 c 0
Obviously, the first formula of the above equation is the solution formula of static stiffness. When the gap distribution is uniform, that is to say, h 0 is a constant, it can be simplified to obtain P 1 s = 0 by combining with the Laplace Equation (4). Therefore, purely parallel plates have no static stiffness, regardless of whether the pressure boundary is uniform. It is clear that if the Laplace Equation (4) is completely unrelated to the translational gap h, the pressure distribution is also completely unrelated to that gap, so the bearing capacity is also completely unrelated to it, that is, there is no stiffness.
When there is no gradient in the static pressure distribution (no additional supply pressure), the simplification can also lead to P 1 s = 0 , the static stiffness of the parallel plate is 0, and there is only a certain squeeze film damping. However, when the frequency approaches 0, so does the disturbance velocity and the damping force.
The transposition of Equation (14) (b) leads to the following:
· h 0 3 P abs 0 P 1 c ω h 1 s = 12 η h 0 P 1 s h 1 s 12 η P abs 0
It should therefore be noted that the damping coefficient (the damping force divided by the product of the frequency and amplitude) is not zero, and the gas then behaves like an incompressible liquid.

5.2. Characteristics of the Thin Film with an Exciting Frequency Approaching Infinity

When the vibration frequency approaches infinity ω + , Equation (16) is simplified:
P 1 s = P abs 0 h 1 s h 0 , P 1 c = 0
At this time, the air film behaves like a rigid spring. The above formula shows that the ultimate stiffness coefficient of the bearing is related to the static (absolute) pressure distribution, and the damping is zero. By bringing the expression of the above formula P 1 s into the stiffness calculation Formula (21), we have
K = Ω P abs 0 h 0 d Ω = Ω P abs 0 d Ω h 0
The above formula shows that K is only related to the static pressure distribution and gap distribution of the gas film, independently of other factors.

5.3. Transfer Function of the Thin Film

The above process shows that a correct analytical solution can be obtained by directly solving the perturbation Reynolds equation. However, the squeeze film effect of a compressible air film is complicated, and the expression of stiffness and damping are nonlinear functions of excitation frequency, which makes them difficult to apply directly. Therefore, using a certain simplified analysis to describe the dynamic characteristics of the air film simply by using several parameters would be highly beneficial for further design analysis.
Based on the characteristics of thin films with exciting frequencies approaching to zero and infinity, and observing the real and imaginary parts of the stiffness in Figure 5, Figure 7, Figure 13, Figure 17a,b and Figure 18a,b, the air film stiffness can be represented by the transfer function shown below:
K b ( s ) = W 1 ( s ) h 1 ( s ) = A s τ s + 1 , K b ( j ω ) = ω 2 A τ τ 2 ω 2 + 1 + j ω A τ 2 ω 2 + 1
By substituting the values at frequencies of zero and infinity, the values of variables A and τ can be obtained:
A = C ( 0 ) , τ = C ( 0 ) K ( )
Table 3 shows the values of all the calculation results found. Since the damping at zero frequency is basically independent of the air supply pressure, and the air film stiffness is positively correlated with the air supply pressure at high frequencies, the time constant is negatively correlated with the air supply pressure.
In order to verify the accuracy of the transfer function, Figure 19 and Figure 20 illustrate the difference in dimensionless stiffness and damping coefficients between the fitted transfer function and numerically calculated Reynolds equation result. The figures show that the fitting results are relatively consistent in the low-frequency region, and the stiffness value obtained in the high frequency region is relatively higher.
The above calculations illustrate the dynamic characteristics of the flat air film in the frequency domain. By transforming them into the time domain, the accuracy of parameter fitting can also be evaluated through a step response. Assume that the boundary structure of an air film is subjected to an impact, and the upper plate produces a displacement h 1 ( t ) = h u u ( t ) in an instant, where u ( t ) represents a unit step function. From the definition of stiffness, the force generated by the air film on the upper plate can be obtained:
W 1 ( t ) = L 1 { K b ( s ) h 1 ( s ) } = L 1 { A s τ s + 1 h u s } = A h u τ u ( t ) e t τ = K h u u ( t ) e t τ
Obviously, the force on the upper plate is an exponential decay function, and its decay constant is the time constant τ .
Figure 21 shows the time-domain response calculated by the FEM and the response curve obtained by inversely transforming the fitted function. Note that this figure uses dimensionless units: the ordinate represents the ratio of the air film force to K h u , and the abscissa represents the ratio of time to the decay constant τ , so the fitting result curve is identical to that in Figure 21a,b. It can be seen that the fitting results can represent the characteristics of thin gas films well. Figure 21b also shows the calculation results of the FEM at different step amplitudes. When the amplitude of the gap variation h u accounts for a small proportion of gas film gaps (in the figure, this value is 0.1% of h 0 ), the response amplitude is highly consistent with the result of linear fitting. As the amplitude of the gap variation accounts for a large proportion, a nonlinear effect gradually appears, but the shape of the response curve also approximates an exponential decay form, even when the amplitude increases to 10% of h 0 .

6. Conclusions

In contrast to the equations representing the squeeze film effect with the same pressure boundary conditions, a group of dynamic Reynolds equations with stiffness and damping pressure were derived, and parallel flat and circular thin films were analyzed considering a high-pressure boundary. The following conclusions can be drawn:
(1) The various dynamic pressure distribution and stiffness curves of the circular thin film were obtained by calculations of the supply pressure both inside and outside the boundary. The characteristics of the squeeze films with and without external pressure were compared.
(2) The squeeze effect of compressible lubricants causes damping at a zero frequency and stiffness at an infinite-frequency. Due to the influence of gas compressibility, the dynamic pressure is still high at high frequencies. The crossover squeeze number, which indicates the frequency at which the stiffness and the damping force are the same, is equal to 10 given a one-dimensional flat plate and 6 given a parallel disc with an atmospheric boundary. The crossover squeeze number increases when the supplying pressure becomes larger.
(3) Based on our models, the fitting results of the time-domain response calculated by the FEM matches well with the response curve obtained by inversely transforming the fitted function. Both accurately represent the characteristics of a thin gas film. As the amplitude of the gap variation accounts for a large proportion, a nonlinear effect gradually appears, but the shape of the response curve approximates an exponential decay form, even when the amplitude increases to 10% of the gas film thickness.

Author Contributions

Conceptualization, Y.W. and Z.Q.; methodology, Y.W.; software, Y.W.; validation, Y.W., J.X. and W.C.; formal analysis, J.X.; investigation, Y.W.; resources, J.X.; data curation, Z.Q.; writing—original draft preparation, Y.W.; writing—review and editing, J.X.; visualization, W.C.; supervision, Z.Q. and B.W.; project administration, B.W.; funding acquisition, Z.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Open Project Program of the State Key Laboratory of Applied Optics (SKLAO2021001A05); National Natural Science Foundation of China (No. 51905130); Heilongjiang Provincial Natural Science Foundation of China (No. LH2020E039).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available on request from the authors.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
η Viscosity of lubricants (Pa·s)
ω Exciting frequency (rad/s)
σ Reciprocal of dimensionless supply air pressure
σ s Squeeze number
σ sc Cut-off frequency
Laplacian operator
Ω Thin-film domain
hFilm thickness ( μ m)
jUnit of imaginary number
lLength of gas film (m)
tTime (s)
CDamping of air bearing (N·s/m)
KStiffness of air bearing (N/m)
K b Dynamic stiffness (N/m)
K Ultimate stiffness (N/m)
R 0 Radius of inner annular film (m)
R 2 Radius of outer annular film (m)
L 1 Ratio of the outer and inner diameters of annular air film
P a Ambient pressure (bar)
P abs Absolute pressure (Pa)
P d High pressure (bar)
WLoad-carrying capacity (N)
Subscripts
0Steady-state item
1Dynamic item
sReal part
cImaginary part
Superscript
˜ Complex variables
*Dimensionless variable

References

  1. Wardle, F. Ultra-Precision Bearings; Elsevier: Cambridge, UK, 2015; pp. 227–306. [Google Scholar]
  2. Walford, T.; Stone, B. The sources of damping in rolling element bearings under oscillating conditions. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 1983, 197, 225–232. [Google Scholar] [CrossRef]
  3. Wensing, J.A. On the Dynamics of Ball Bearings. Ph.D. Thesis, University of Twente, Enschede, The Netherlands, 1998. [Google Scholar]
  4. Bao, M.; Yang, H. Squeeze film air damping in MEMS. Sens. Actuators Phys. 2007, 136, 3–27. [Google Scholar] [CrossRef]
  5. Albarbar, A.; Mekid, S.; Starr, A.; Pietruszkiewicz, R. Suitability of MEMS accelerometers for condition monitoring: An experimental study. Sensors 2008, 8, 784–799. [Google Scholar] [CrossRef]
  6. Wenbin, L.; Jianfeng, H.; Jie, F.; Liyun, C.; Chunyan, Y.; Wenjing, W. Simulation of the engagement of carbon fabric wet clutch: Analytical and experimental comparison. Tribol. Int. 2015, 90, 502–508. [Google Scholar] [CrossRef]
  7. Jang, S. Manipulating Frictional Performance of Wet Clutch Engagement through Material Properties and Operating Conditions. Lubricants 2022, 10, 225. [Google Scholar] [CrossRef]
  8. Kane, M.; Do, M.T.; Cerezo, V.; Rado, Z.; Khelifi, C. Contribution to pavement friction modelling: An introduction of the wetting effect. Int. J. Pavement Eng. 2019, 20, 965–976. [Google Scholar] [CrossRef]
  9. Gronqvist, R. Mechanisms of friction and assessment of slip resistance of new and used footwear soles on contaminated floors. Ergonomics 1995, 38, 224–241. [Google Scholar] [CrossRef]
  10. Zeidan, F.Y.; Andres, L.S.; Vance, J.M. Design Furthermore, Application of Squeeze Film Dampers in Rotating Machinery. In Proceedings of the 25th Turbomachinery Symposium, Texas A&M University, Turbomachinery Laboratories, College Station, TA, USA, 13 December 1996; pp. 1–20. [Google Scholar] [CrossRef]
  11. Biet, M.; Giraud, F.; Lemaire-Semail, B. Squeeze film effect for the design of an ultrasonic tactile plate. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2007, 54, 2678–2688. [Google Scholar] [CrossRef] [PubMed]
  12. Sadd, M.H.; Stiffler, A.K. Squeeze Film Dampers: Amplitude Effects at Low Squeeze Numbers. J. Eng. Ind. 1975, 97, 1366–1370. [Google Scholar] [CrossRef]
  13. Salbu, E.O.J. Compressible Squeeze Films and Squeeze Bearings. J. Basic Eng. 1964, 86, 355–364. [Google Scholar] [CrossRef]
  14. Langlois, W. Isothermal squeeze films. Q. Appl. Math. 1962, 20, 131–150. [Google Scholar] [CrossRef]
  15. Gross, W.A. Gas Film Lubrication; John Wi1ey and Sons, Inc.: New York, NY, USA, 1962; pp. 308–317. [Google Scholar]
  16. Pan, C.H.T. The Gaseous Squeeze-Film at Moderately Large Squeeze Numbers. J. Basic Eng. 1970, 92, 766–781. [Google Scholar] [CrossRef]
  17. Blech, J.J. On isothermal squeeze films. J. Tribol. 1983, 105, 615. [Google Scholar] [CrossRef]
  18. Talukder, H.; Stowell, T. Pneumatic hammer in an externally pressurized orifice-compensated air journal bearing. Tribol. Int. 2003, 36, 585–591. [Google Scholar] [CrossRef]
  19. Grossman, R.L. Application of Flow and Stability Theory to the Design of Externally Pressurized Spherical Gas Bearings. J. Basic Eng. 1963, 85, 495–502. [Google Scholar] [CrossRef]
  20. Holmes, M.H. Introduction to Perturbation Methods; Springer Science & Business Media: New York, NY, USA, 2012; Volume 20, pp. 57–73. [Google Scholar]
  21. Licht, L.; Elrod, H. A Study of the Stability of Externally Pressurized Gas Bearings. J. Appl. Mech. 1960, 27, 250–258. [Google Scholar] [CrossRef]
  22. Boffey, D.A. A Study of the Stability of an Externally-Pressurized Gas-Lubricated Thrust Bearing with a Flexible Damped Support. J. Lubr. Technol. 1978, 100, 364–368. [Google Scholar] [CrossRef]
  23. Boffey, D.A.; Desai, D.M. An experimental investigation into the rubber-stabilization of an externally-pressurized air-lubricated thrust bearing. J. Lubr. Technol. 1980, 102, 65–70. [Google Scholar] [CrossRef]
  24. Stiffler, A.K. Analysis of the Stiffness and Damping of an Inherently Compensated, Multiple-Inlet, Circular Thrust Bearing. J. Lubr. Technol. 1974, 96, 329–336. [Google Scholar] [CrossRef]
  25. Stiffler, A.K.; Smith, D.M. Dynamic Characteristics of an Inherently Compensated, Square, Gas Film Bearing. J. Lubr. Technol. 1975, 97, 52–62. [Google Scholar] [CrossRef]
  26. Miyatake, M.; Yoshimoto, S. Numerical investigation of static and dynamic characteristics of aerostatic thrust bearings with small feed holes. Tribol. Int. 2010, 43, 1353–1359. [Google Scholar] [CrossRef]
  27. Bhat, N.; Kumar, S.; Tan, W.; Narasimhan, R.; Low, T.C. Performance of inherently compensated flat pad aerostatic bearings subject to dynamic perturbation forces. Precis. Eng. 2012, 36, 399–407. [Google Scholar] [CrossRef]
  28. Colombo, F.; Raparelli, T.; Trivella, A.; Viktorov, V. Identification of a Lumped Parameters Numerical Model of Gas Bearings: Analysis of 1D Parallel Plates. In Proceedings of the Advances in Italian Mechanism Science, Cassino, Italy, 10 December 2019; pp. 467–473. [Google Scholar] [CrossRef]
  29. Belforte, G.; Colombo, F.; Raparelli, T.; Trivella, A.; Viktorov, V. Study of the Static and Dynamic Performance of Rectangular Air Pads by Means of Lumped Parameters Models. In Proceedings of the ASME 2012 11th Biennial Conference on Engineering Systems Design and Analysis, Nantes, France, 2–4 July 2012. [Google Scholar]
  30. Colombo, F.; Lentini, L.; Raparelli, T.; Trivella, A.; Viktorov, V. Dynamic characterisation of rectangular aerostatic pads with multiple inherent orifices. Tribol. Lett. 2018, 66, 1–13. [Google Scholar] [CrossRef]
  31. Moradi, M.; Colombo, F.; Raparelli, T.; Trivella, A.; Viktorov, V. Dynamic lumped model of externally pressurized rectangular air bearings. Precis. Eng. 2019, 56, 101–112. [Google Scholar] [CrossRef]
  32. Hamrock, B.J.; Schmid, S.R.; Jacobson, B.O. Fundamentals of Fluid Film Lubrication; Marcel Dekker, Inc.: New York, NY, USA, 2004; pp. 202–228. [Google Scholar]
  33. Wu, Y.; Qiao, Z.; Xue, J.; Liu, Y.; Wang, B. Root iterative method for static performance analysis of aerostatic thrust bearings with multiple pocketed orifice-type restrictors based on ANSYS. Ind. Lubr. Tribol. 2019, 72, 165–171. [Google Scholar] [CrossRef]
  34. Castelli, V.; Pirvics, J. Review of numerical methods in gas-bearing film analysis. Tribology 1969, 2, 76. [Google Scholar] [CrossRef]
  35. Shi, J.; Cao, H.; Jin, X. Investigation on the static and dynamic characteristics of 3-DOF aerostatic thrust bearings with orifice restrictor. Tribol. Int. 2019, 138, 435–449. [Google Scholar] [CrossRef]
  36. Bathe, K.J. Finite Element Procedures; Prentice-Hall, Inc.: Hoboken, NJ, USA, 1996; pp. 642–694. [Google Scholar]
  37. Wu, Y.; Qiao, Z.; Xue, J.; Wang, B.; Chen, W. A newly treated boundary conditions to enhance accuracy of finite element analysis for orifice-type aerostatic bearings. Adv. Eng. Softw. 2022, 173, 103277. [Google Scholar] [CrossRef]
Figure 1. Squeeze film of the one-dimensional constant gap flow model: (a) infinite width flat air film; and (b) circular air film.
Figure 1. Squeeze film of the one-dimensional constant gap flow model: (a) infinite width flat air film; and (b) circular air film.
Mathematics 11 00742 g001
Figure 2. Distribution of P 1 s inside of the 1D parallel slider.
Figure 2. Distribution of P 1 s inside of the 1D parallel slider.
Mathematics 11 00742 g002
Figure 3. Distribution of P 1 c inside the 1D parallel slider.
Figure 3. Distribution of P 1 c inside the 1D parallel slider.
Mathematics 11 00742 g003
Figure 4. Dynamic coefficients of the 1D parallel slider.
Figure 4. Dynamic coefficients of the 1D parallel slider.
Mathematics 11 00742 g004
Figure 5. Dynamic stiffness modulus of the 1D parallel slider.
Figure 5. Dynamic stiffness modulus of the 1D parallel slider.
Mathematics 11 00742 g005
Figure 6. Dynamic coefficients of the 1D parallel disc.
Figure 6. Dynamic coefficients of the 1D parallel disc.
Mathematics 11 00742 g006
Figure 7. Dynamic stiffness modulus of the 1D parallel disc.
Figure 7. Dynamic stiffness modulus of the 1D parallel disc.
Mathematics 11 00742 g007
Figure 8. Flowchart of the calculation procedure.
Figure 8. Flowchart of the calculation procedure.
Mathematics 11 00742 g008
Figure 9. Squeeze film of a one-dimensional constant gap flow model with different pressure boundaries: (a) infinite width flat air film; and (b) circular air film.
Figure 9. Squeeze film of a one-dimensional constant gap flow model with different pressure boundaries: (a) infinite width flat air film; and (b) circular air film.
Mathematics 11 00742 g009
Figure 10. Distribution of P 1 s inside of the 1D parallel slider with high-pressure in the left.
Figure 10. Distribution of P 1 s inside of the 1D parallel slider with high-pressure in the left.
Mathematics 11 00742 g010
Figure 11. Distribution of P 1 c inside of 1D parallel slider with high-pressure in the left.
Figure 11. Distribution of P 1 c inside of 1D parallel slider with high-pressure in the left.
Mathematics 11 00742 g011
Figure 12. Dynamic coefficients of the 1D parallel slider with supply pressure.
Figure 12. Dynamic coefficients of the 1D parallel slider with supply pressure.
Mathematics 11 00742 g012
Figure 13. Dynamic stiffness modulus of the 1D parallel slider with supply pressure.
Figure 13. Dynamic stiffness modulus of the 1D parallel slider with supply pressure.
Mathematics 11 00742 g013
Figure 14. Distribution of P 1 s inside of the 1D parallel ring with high-pressure inside.
Figure 14. Distribution of P 1 s inside of the 1D parallel ring with high-pressure inside.
Mathematics 11 00742 g014
Figure 15. Distribution of P 1 c inside of the 1D parallel ring with high-pressure inside.
Figure 15. Distribution of P 1 c inside of the 1D parallel ring with high-pressure inside.
Mathematics 11 00742 g015
Figure 16. Dynamic coefficients of the 1D ring with supply pressure under different radius ratios.
Figure 16. Dynamic coefficients of the 1D ring with supply pressure under different radius ratios.
Mathematics 11 00742 g016
Figure 17. Dynamic stiffness modulus of the 1D ring with supply pressure under different radius ratios (flows from inside).
Figure 17. Dynamic stiffness modulus of the 1D ring with supply pressure under different radius ratios (flows from inside).
Mathematics 11 00742 g017
Figure 18. Dynamic stiffness modulus of the 1D ring with supply pressure under different radius ratios (flows from outside).
Figure 18. Dynamic stiffness modulus of the 1D ring with supply pressure under different radius ratios (flows from outside).
Mathematics 11 00742 g018
Figure 19. Fitted transfer function of dimensionless stiffness.
Figure 19. Fitted transfer function of dimensionless stiffness.
Mathematics 11 00742 g019
Figure 20. Fitted transfer function of dimensionless damping.
Figure 20. Fitted transfer function of dimensionless damping.
Mathematics 11 00742 g020
Figure 21. The time-domain response calculated by the FEM and obtained by inversely transforming the fitted function.
Figure 21. The time-domain response calculated by the FEM and obtained by inversely transforming the fitted function.
Mathematics 11 00742 g021
Table 1. Parameters of a typical aerostatic guideway.
Table 1. Parameters of a typical aerostatic guideway.
NameViscosityDisturbance FrequencyCharacteristic LengthFilm Thickness
Value 17.9 × 10 6 Pa · s 500 Hz 0.1 m 10 × 10 6 m
Table 2. Parameters of a typical aerostatic thrust bearing.
Table 2. Parameters of a typical aerostatic thrust bearing.
NameViscosityDisturbance FrequencyCharacteristic RadiusFilm Thickness
Value 17.9 × 10 6 Pa · s 500 Hz 0.05 m 10 × 10 6 m
Table 3. Parameters of a typical aerostatic thrust bearing.
Table 3. Parameters of a typical aerostatic thrust bearing.
ModelBoundary Setting A = C ( 0 ) K ( ) τ = C ( 0 ) / K ( )
1D flat filmAmbient pressure η b l 3 / ( h 0 3 ) b l P a / h 0 η l 2 / ( h 0 2 P a )
1D circular filmAmbient pressure 3 η π R 2 4 / ( 2 h 0 3 ) π R 2 2 P a / h 0 3 η R 2 2 / ( 2 h 0 2 P a )
1D flat film σ s = 1 / 3 1.028 η b l 3 / ( h 0 3 ) 2.148 b l P a / h 0 0.479 η l 2 / ( h 0 2 P a )
1D flat film σ s = 1 / 5 1.046 η b l 3 / ( h 0 3 ) 3.393 b l P a / h 0 0.414 η l 2 / ( h 0 2 P a )
1D flat film σ s = 1 / 7 1.054 η b l 3 / ( h 0 3 ) 4.652 b l P a / h 0 0.227 η l 2 / ( h 0 2 P a )
1D flat film σ s = 1 / 9 1.058 η b l 3 / ( h 0 3 ) 5.911 b l P a / h 0 0.179 η l 2 / ( h 0 2 P a )
1D annular film σ s = 1 , L 2 = 1.25 η π ( R 2 2 R 0 2 ) ( R 2 R 0 ) 2 h 0 3 π ( R 2 2 R 0 2 ) P a h 0 η ( R 2 2 R 0 2 ) h 0 2 P a
1D annular film σ s = 1 / 3 , L 2 = 1.25 1.031 η π ( R 2 2 R 0 2 ) ( R 2 R 0 ) 2 h 0 3 2.008 π ( R 2 2 R 0 2 ) P a h 0 0.513 η ( R 2 2 R 0 2 ) h 0 2 P a
1D annular film σ s = 1 / 5 , L 2 = 1.25 1.050 η π ( R 2 2 R 0 2 ) ( R 2 R 0 ) 2 h 0 3 3.158 π ( R 2 2 R 0 2 ) P a h 0 0.336 η ( R 2 2 R 0 2 ) h 0 2 P a
1D annular film σ s = 1 / 7 , L 2 = 1.25 1.058 η π ( R 2 2 R 0 2 ) ( R 2 R 0 ) 2 h 0 3 4.324 π ( R 2 2 R 0 2 ) P a h 0 0.245 η ( R 2 2 R 0 2 ) h 0 2 P a
1D annular film σ s = 1 , L 2 = 2.0 η π ( R 2 2 R 0 2 ) ( R 2 R 0 ) 2 h 0 3 π ( R 2 2 R 0 2 ) P a h 0 η ( R 2 2 R 0 2 ) h 0 2 P a
1D annular film σ s = 1 / 3 , L 2 = 2.0 1.041 η π ( R 2 2 R 0 2 ) ( R 2 R 0 ) 2 h 0 3 1.901 π ( R 2 2 R 0 2 ) P a h 0 0.547 η ( R 2 2 R 0 2 ) h 0 2 P a
1D annular film σ s = 1 / 5 , L 2 = 2.0 1.064 η π ( R 2 2 R 0 2 ) ( R 2 R 0 ) 2 h 0 3 2.940 π ( R 2 2 R 0 2 ) P a h 0 0.362 η ( R 2 2 R 0 2 ) h 0 2 P a
1D annular film σ s = 1 / 7 , L 2 = 2.0 1.075 η π ( R 2 2 R 0 2 ) ( R 2 R 0 ) 2 h 0 3 4.003 π ( R 2 2 R 0 2 ) P a h 0 0.268 η ( R 2 2 R 0 2 ) h 0 2 P a
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wu, Y.; Xue, J.; Qiao, Z.; Chen, W.; Wang, B. The Squeeze Film Effect with a High-Pressure Boundary in Aerostatic Bearings. Mathematics 2023, 11, 742. https://doi.org/10.3390/math11030742

AMA Style

Wu Y, Xue J, Qiao Z, Chen W, Wang B. The Squeeze Film Effect with a High-Pressure Boundary in Aerostatic Bearings. Mathematics. 2023; 11(3):742. https://doi.org/10.3390/math11030742

Chicago/Turabian Style

Wu, Yangong, Jiadai Xue, Zheng Qiao, Wentao Chen, and Bo Wang. 2023. "The Squeeze Film Effect with a High-Pressure Boundary in Aerostatic Bearings" Mathematics 11, no. 3: 742. https://doi.org/10.3390/math11030742

APA Style

Wu, Y., Xue, J., Qiao, Z., Chen, W., & Wang, B. (2023). The Squeeze Film Effect with a High-Pressure Boundary in Aerostatic Bearings. Mathematics, 11(3), 742. https://doi.org/10.3390/math11030742

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