# An Optimization Method for the Main Influence Factors on the Performance of Aerostatic Spindle

Authors: Dongju Chen; You Zhao; Ri Pan; Jinwei Fan
Journal: Shock and Vibration
 (2019)

---

## Abstract

This paper proposes a novel optimization method to improve the performance of aerostatic spindle. The optimization process is based on the main factors affecting the load and stiffness of aerostatic spindle, which are determined by cross-correlation analysis. The maximum values of load capacity and stiffness are obtained from the relationship between the optimization function and the main influence factors. After optimization, the corresponding natural frequency has been improved, which indicates that the optimization method proposed in this paper is more efficient and reliable in the field of precision machining. The results will enforce the industrial application of aerostatic spindle system.


---

## Body

## 1. Introduction

The aerostatic spindle is widely used in precision and ultraprecision machining by the advantages of low friction and high accuracy [1]. When the aerostatic bearing is working, the air supply is provided by an external air compressor. Then, the gas enters the air cavity through the orifice to form a layer of gas film which can support the workpiece to machine and achieve frictionless lubrication. The bearing capacity of the aerostatic spindle mainly comes from the gas film between the aerostatic bearing and the spindle. For ultraprecision machining, the stability of aerostatic spindle should be guaranteed. In fact, the performance of the aerostatic spindle is completely different due to some factors, such as structural parameters of air cavity, orifice diameter, gas film thickness, and air supply pressure. And the performance of the aerostatic spindle is a key factor affecting the quality of the machine [2]. Besides, it will be also directly reflected to the surface of the workpiece and affect the quality of the workpiece. Therefore, it is important to study the performance of the aerostatic spindle to improve the precision of ultraprecision machining.There are two major aspects for the optimization of aerostatic spindle: one focuses on the static characteristics analysis of aerostatic spindle to optimize basic parameters [3]. Aguirre et al. [4] proposed an active compensation strategy to achieve nanometer position control, which can improve the performance of the thrust bearing. Chang and Jeng [5] researched the performance of aerostatic bearing by optimizing the bearing diameter, orifice diameter, orifice structure, gas supply pressure, and so on. Li et al. [6] proposed an optimization model by the maximum Mach number to improve dynamic stability. However, this is only one aspect of the optimization of aerostatic spindle.Another aspect of research is that the final performance is determined by optimization methods. Shie and Shih [7] obtained a numerical solution of the pressure distribution between the surface of the aerostatic bearing and the worktable using the finite difference method. Li et al. [8, 9] proposed a coupling optimization method to optimize the spindle structure by magnetic circuit calculation and magnetic field calculation. Wang and Chang [10] used Pareto ranking method and genetic algorithm to optimize the stiffness and air mass flow of porous air bearing. Besides, they used the hypercube-diving method to research the multiobjective optimization of porous air bearing [11]. These optimization methods have different objective function for the different structure of gas cavity and the aerostatic bearing with orifice has a high load capacity. The performance of gas film in microscale is different from the macroscale, and the corresponding optimization in microscale of structure of gas film is seldom.In this paper, an optimization method is proposed to optimize the factors affecting the performance of the aerostatic spindle. These factors include gas film thickness, diameter of orifice, and the diameter of stomatic distribution circle. The relationship between load capacity, stiffness, and the main influence factors determined by correlation analysis is deduced. The results are verified from static and dynamic characteristics. And the corresponding simulation and experiment of the bearing capacity and stiffness of the aerostatic spindle are carried out to achieve enhanced performance in widely application.


## 2. Pressure Distribution of Aerostatic Spindle

### 2.1. Structure of Aerostatic Spindle

The circular flat pad aerostatic bearing with eight orifice-type restrictor is shown in Figure1. According to the structure characteristics and actual working situation of the aerostatic bearing, the main design parameters contain bearing clearance h, orifice diameter d3, the diameter of orifice distribution circle d2, the rotate speed of the spindle, and so on. The bearing clearance h, orifice diameter d3, and the diameter of orifice distribution circle d2 are considered as design variables because they have more influence on the aerostatic bearing performance, and the ranges are 5–20 μm for h, 2 mm–5 mm for d3, and 75–85 mm for d2, respectively. The air supply pressure ps is given between 0.3 MPa and 0.6 MPa.Figure 1
Structure of bearing.
(a)
(b)


### 2.2. Flow Field Analysis in Microscale

Microfactors will lead to decrease of the gas film density and increase of the gas compressibility. Therefore, the traditional Reynolds equation is no longer applicable to describe the microscale flow. The flow factorQ embodied microfactor whose equation is Q = C1 + C2/PH is introduced into the Reynolds equation, where C1 and C2 can be obtained from the literature [12]. So, the Reynolds equation after introducing the flow factor Q is as follows:(1)∂∂xph3Q∂p∂x+∂∂zph3Q∂p∂z=6μU∂ph∂x+12μ∂ph∂t,where p is the film pressure, h is the bearing film thickness, μ is the dynamic viscosity of the lubricating gas, U is the circumferential speed of the bearing surface, and U=ωR0. ω is the angular velocity of spindle rotation, R0 is the bearing radius, and t is the time. Equation (1) is converted into the polar coordinate for the convenience of calculation.(2)∂∂rph3Q∂p∂r+∂r2∂θph3Q∂p∂θ=6μU∂ph∂r+12μ∂ph∂t.Equation (2) is dimensionless for the convenience of calculation. The characteristic pressure of the film is p0 and the characteristic thickness of film is bearing radius clearance h0. Let p=pP0, p0=6μR02/h02, r=RR0, and h=hH0. The characteristic length of the distribution circle radius of the orifice is R. P is nondimensional pressure, and H is nondimensional bearing film thickness. The time term is ignored in solving the static characteristic of bearing. So, the dimensionless form of Reynolds equation with flow factor Q is as follows:(3)∂∂RPH3Q∂P∂R+∂R2∂θPH3Q∂P∂θ=ω∂PH∂R.Two key steps for solving Reynolds equation by finite difference method are (1) dividing the solution domain into nodal grids and (2) transforming partial differential Reynolds equation into difference equation. Finally, the pressure distribution ofpa can be obtained by substituting boundary conditions, as shown in Figure 2.Figure 2
Performance of aerostatic spindle in axial direction. (a) Gas pressure distribution. (b) Relation between bearing stiffness and the gas film clearance.
(a)
(b)Then, use the disturbance method to establish unsteady Reynolds equation. The basic idea is to assume that when the input disturbance is a sinusoidal or cosine function, the output is the superposition of sinusoidal and cosine functions when the object of study is a first-order linear system. So, the integral of the pressure disturbance term in phase with the input disturbance can be defined as dynamic stiffness, and the integral of the pressure disturbance term with phase lag can be defined as dynamic damping coefficient.(4)h=h0+h0sinω1t,p=pa+pbsinω1t+pccosω1t,where pa is the film pressure in the steady state, h0 is steady-state film thickness, ω1 is the perturbation frequency, pb is the pressure perturbation term in phase with the input perturbation, and pc is the pressure perturbation term with phase lag. Substituting equation (4) in equation (1), let pA=pa⋅pb, pB=pa⋅pc, the disturbance quantities in quadratic and above are removed. The perturbed Reynolds equation is as follows:(5)∂2pA∂r2+1r∂pA∂r+1r2∂2pA∂θ2=−12μω1PBpah02,∂2pB∂r2+1r∂pB∂r+1r2∂2pB∂θ2=−12μω1pa2+pApah02.The dynamic performance is solved by finite difference method. First, the perturbation Reynolds equation is discretized and then meshed. Finally, the boundary conditions of steady state and perturbation are introduced into the equation to solve the problem. By substituting the steady pressurepa and the steady film thickness h0 obtained from the steady-state solution into the perturbation equation, the dynamic pressure and the dynamic film gap can be obtained. pb and pc can be obtained by solving equation (5).


### 2.3. Performance Analysis of Aerostatic Spindle

According to the pressure distributions, the load capacityW is calculated by integrating the pressure along the lower surface of bearing.(6)W=8∫0π/4∫d1/2d4/2prdrdθ.At the same time, the stiffness of bearing can be obtained by calculating the increment ratio of the bearing capacity and bearing clearance.(7)K=dWdh.The relationship between stiffness and bearing clearance is shown in Figure2(b). The bearing capacity and stiffness of the aerostatic spindle are affected by many factors, including structural and working parameters. Next, we will evaluate the effect of the main factors on the spindle performance.


## 2.1. Structure of Aerostatic Spindle

The circular flat pad aerostatic bearing with eight orifice-type restrictor is shown in Figure1. According to the structure characteristics and actual working situation of the aerostatic bearing, the main design parameters contain bearing clearance h, orifice diameter d3, the diameter of orifice distribution circle d2, the rotate speed of the spindle, and so on. The bearing clearance h, orifice diameter d3, and the diameter of orifice distribution circle d2 are considered as design variables because they have more influence on the aerostatic bearing performance, and the ranges are 5–20 μm for h, 2 mm–5 mm for d3, and 75–85 mm for d2, respectively. The air supply pressure ps is given between 0.3 MPa and 0.6 MPa.Figure 1
Structure of bearing.
(a)
(b)


## 2.2. Flow Field Analysis in Microscale

Microfactors will lead to decrease of the gas film density and increase of the gas compressibility. Therefore, the traditional Reynolds equation is no longer applicable to describe the microscale flow. The flow factorQ embodied microfactor whose equation is Q = C1 + C2/PH is introduced into the Reynolds equation, where C1 and C2 can be obtained from the literature [12]. So, the Reynolds equation after introducing the flow factor Q is as follows:(1)∂∂xph3Q∂p∂x+∂∂zph3Q∂p∂z=6μU∂ph∂x+12μ∂ph∂t,where p is the film pressure, h is the bearing film thickness, μ is the dynamic viscosity of the lubricating gas, U is the circumferential speed of the bearing surface, and U=ωR0. ω is the angular velocity of spindle rotation, R0 is the bearing radius, and t is the time. Equation (1) is converted into the polar coordinate for the convenience of calculation.(2)∂∂rph3Q∂p∂r+∂r2∂θph3Q∂p∂θ=6μU∂ph∂r+12μ∂ph∂t.Equation (2) is dimensionless for the convenience of calculation. The characteristic pressure of the film is p0 and the characteristic thickness of film is bearing radius clearance h0. Let p=pP0, p0=6μR02/h02, r=RR0, and h=hH0. The characteristic length of the distribution circle radius of the orifice is R. P is nondimensional pressure, and H is nondimensional bearing film thickness. The time term is ignored in solving the static characteristic of bearing. So, the dimensionless form of Reynolds equation with flow factor Q is as follows:(3)∂∂RPH3Q∂P∂R+∂R2∂θPH3Q∂P∂θ=ω∂PH∂R.Two key steps for solving Reynolds equation by finite difference method are (1) dividing the solution domain into nodal grids and (2) transforming partial differential Reynolds equation into difference equation. Finally, the pressure distribution ofpa can be obtained by substituting boundary conditions, as shown in Figure 2.Figure 2
Performance of aerostatic spindle in axial direction. (a) Gas pressure distribution. (b) Relation between bearing stiffness and the gas film clearance.
(a)
(b)Then, use the disturbance method to establish unsteady Reynolds equation. The basic idea is to assume that when the input disturbance is a sinusoidal or cosine function, the output is the superposition of sinusoidal and cosine functions when the object of study is a first-order linear system. So, the integral of the pressure disturbance term in phase with the input disturbance can be defined as dynamic stiffness, and the integral of the pressure disturbance term with phase lag can be defined as dynamic damping coefficient.(4)h=h0+h0sinω1t,p=pa+pbsinω1t+pccosω1t,where pa is the film pressure in the steady state, h0 is steady-state film thickness, ω1 is the perturbation frequency, pb is the pressure perturbation term in phase with the input perturbation, and pc is the pressure perturbation term with phase lag. Substituting equation (4) in equation (1), let pA=pa⋅pb, pB=pa⋅pc, the disturbance quantities in quadratic and above are removed. The perturbed Reynolds equation is as follows:(5)∂2pA∂r2+1r∂pA∂r+1r2∂2pA∂θ2=−12μω1PBpah02,∂2pB∂r2+1r∂pB∂r+1r2∂2pB∂θ2=−12μω1pa2+pApah02.The dynamic performance is solved by finite difference method. First, the perturbation Reynolds equation is discretized and then meshed. Finally, the boundary conditions of steady state and perturbation are introduced into the equation to solve the problem. By substituting the steady pressurepa and the steady film thickness h0 obtained from the steady-state solution into the perturbation equation, the dynamic pressure and the dynamic film gap can be obtained. pb and pc can be obtained by solving equation (5).


## 2.3. Performance Analysis of Aerostatic Spindle

According to the pressure distributions, the load capacityW is calculated by integrating the pressure along the lower surface of bearing.(6)W=8∫0π/4∫d1/2d4/2prdrdθ.At the same time, the stiffness of bearing can be obtained by calculating the increment ratio of the bearing capacity and bearing clearance.(7)K=dWdh.The relationship between stiffness and bearing clearance is shown in Figure2(b). The bearing capacity and stiffness of the aerostatic spindle are affected by many factors, including structural and working parameters. Next, we will evaluate the effect of the main factors on the spindle performance.


## 3. Degree Analysis of the Influence Factors of the Aerostatic Bearing

### 3.1. Influence Factors on the Performance of Aerostatic Spindle

Theoretical calculations show that the factors affecting performance of aerostatic spindle include the bearing clearanceh, orifice diameter d3 and the diameter of orifice distribution circle d2, the depth of the air cavity d and the pressure supply ps. However, the influence extent of these factors is different. For example, some factors improve the performance of the aerostatic bearing, and others reduce the performance. Therefore, the first step is to assess the influence extent of these factors.


### 3.2. Analysis of the Impact Extent of Influence Factors

In order to analyze the influence of the single factor more intuitively on the final performance of aerostatic bearing, the cross-correlation analysis method is proposed, which uses similarity of two series as a function of the lag of one relative to the other to evaluate correlation extent of the single parameter on the performance of aerostatic spindle. This is also known as a sliding dot product or sliding inner-product. To characterize the correlation between two random variablesx and y1,2, which produce yield data sets x=x1,…,xn and y1,2=y1,…,yn, the correlation coefficient is usually used to express the degree of linear correlation between the two variables [13].(8)Rxyk=∑i=1nxi−x¯yki−yk¯∑i=1nxi−x¯2∑i=1nyki−yk¯2,k=1,2,where x is the influence factors including the bearing clearance h, orifice diameter d3, the diameter of stomatic distribution circle d2, the depth of the air cavity d, and the pressure supply ps. The ranges are 5 μm–20 μm for h, 2 mm–5 mm for d3, 75 mm–85 mm for d2, 0.75 mm–1.0 mm for d, and 0.3 MPa–0.6 MPa for ps. i is the signal number of correlation analysis and y1 and y2 are the load capacity and stiffness of aerostatic spindle, respectively.After analysis the above, the influence extent of single factor on the stiffness of the aerostatic bearing is obtained, as shown in Figure3. From this, we can get that the correlation coefficients of the influence factors h, d3, and d2 are all more than 0.7; however, the correlation values of influence factor d and ps are both less than 0.6, which indicates the main influence factors are the bearing clearance h, orifice diameter d3, and the diameter of stomatic distribution circle d2. Therefore, next, we will optimize the three parameters.Figure 3
Correlation coefficient of the stiffness and influence factors.


## 3.1. Influence Factors on the Performance of Aerostatic Spindle

Theoretical calculations show that the factors affecting performance of aerostatic spindle include the bearing clearanceh, orifice diameter d3 and the diameter of orifice distribution circle d2, the depth of the air cavity d and the pressure supply ps. However, the influence extent of these factors is different. For example, some factors improve the performance of the aerostatic bearing, and others reduce the performance. Therefore, the first step is to assess the influence extent of these factors.


## 3.2. Analysis of the Impact Extent of Influence Factors

In order to analyze the influence of the single factor more intuitively on the final performance of aerostatic bearing, the cross-correlation analysis method is proposed, which uses similarity of two series as a function of the lag of one relative to the other to evaluate correlation extent of the single parameter on the performance of aerostatic spindle. This is also known as a sliding dot product or sliding inner-product. To characterize the correlation between two random variablesx and y1,2, which produce yield data sets x=x1,…,xn and y1,2=y1,…,yn, the correlation coefficient is usually used to express the degree of linear correlation between the two variables [13].(8)Rxyk=∑i=1nxi−x¯yki−yk¯∑i=1nxi−x¯2∑i=1nyki−yk¯2,k=1,2,where x is the influence factors including the bearing clearance h, orifice diameter d3, the diameter of stomatic distribution circle d2, the depth of the air cavity d, and the pressure supply ps. The ranges are 5 μm–20 μm for h, 2 mm–5 mm for d3, 75 mm–85 mm for d2, 0.75 mm–1.0 mm for d, and 0.3 MPa–0.6 MPa for ps. i is the signal number of correlation analysis and y1 and y2 are the load capacity and stiffness of aerostatic spindle, respectively.After analysis the above, the influence extent of single factor on the stiffness of the aerostatic bearing is obtained, as shown in Figure3. From this, we can get that the correlation coefficients of the influence factors h, d3, and d2 are all more than 0.7; however, the correlation values of influence factor d and ps are both less than 0.6, which indicates the main influence factors are the bearing clearance h, orifice diameter d3, and the diameter of stomatic distribution circle d2. Therefore, next, we will optimize the three parameters.Figure 3
Correlation coefficient of the stiffness and influence factors.


## 4. Numerical Model

### 4.1. Relationship between Performance and Influence Factors

In order to find the relationship between the performance of aerostatic spindle and the influence factors. An arrayX with these design variables which include the bearing clearance h, orifice diameter d3, and the diameter of stomatic distribution circle d2 is determined, X=x1,x2,x3=h,d3,d2. The industrial applications require high bearing capacity and stiffness, but their values vary with design variables. Therefore, the multiobjective optimization containing the bearing capacity and stiffness is proposed, and the multiobjective optimization functions Fx including the maximum bearing capacity function f1x and the maximum stiffness function f2x are established.(9)Fx=N1f1x+N2f2x,where N1 and N2 represent weighting factors which reflect the specific gravity of bearing capacity and stiffness in the total objective function, respectively. And, they are calculated by the equation:(10)Ni=1fiX∗,where fiX∗ represents the optimal solution of single constraint problem for the i-th subobjective as the objective function.


### 4.2. Function Determination of the Performance of Aerostatic Spindle

According to the assumption of uniform radial flow of gas, the motion equation in cylindrical coordinates can be simplified as follows:(11)∂p∂r=μ∂2vr∂z2,∂p∂θ=0,∂p∂z=0.By integrating equation (11) with z twice and combining the equivalent mass flow equation and the gas state equation, we can get the expression of pressure p as follows:(12)p12=pa2−pa2−ps2lnr/d4lnd2/d41/2,d4<r<d2,(13)p22=p02−p02−pa2lnr/d2lnd1/d21/2,d2<r<d1.The expression of bearing capacity can be obtained by substituting equations (12) and (13) in equation (6):(14)W=pa−psπd22−d422lnd2/d4+p0−paπd12−d222lnd1/d2−πd12−d42p0+πd32pa.Then, we can getf1x.(15)f1x=pa−psπx32−d422lnx3/d4+p0−paπd12−x322lnd1/x3−πd12−d42p0+πx22pa.And then, we define the gauge pressure ratio by using Powell’s method.(16)kg=pa−p0ps−p0.So, we write the bearing capacity as follows:(17)W=kgpa−psπd22−d422lnd2/d4+p0−paπd12−d222lnd1/d2−πd12−d42p0+πd32pa.Then, we introduce the slit factorG.(18)G=pa/ps1−pa/ps1/21+pa/ps⋅24μ2RT1/2pa⋅CDd32lnd1/d28h3,(19)kg=21+1+4/G21/2.According to the equation (7), we can get the objective function f2x of the stiffness as follows:(20)f2x=−0.98x1ps−paπ2x2x32−d42lnx3/d4,where x1 represents the gas film thickness h, x2 is the orifice diameter d3, and x3 is the diameter of distribution circle d2. d1 represents the external diameter of bearing, d4 represents the inner diameter of bearing, and ps represents the air supply pressure. p0 represents the atmospheric pressure and the discharge coefficient CD = 0.8.


## 4.1. Relationship between Performance and Influence Factors

In order to find the relationship between the performance of aerostatic spindle and the influence factors. An arrayX with these design variables which include the bearing clearance h, orifice diameter d3, and the diameter of stomatic distribution circle d2 is determined, X=x1,x2,x3=h,d3,d2. The industrial applications require high bearing capacity and stiffness, but their values vary with design variables. Therefore, the multiobjective optimization containing the bearing capacity and stiffness is proposed, and the multiobjective optimization functions Fx including the maximum bearing capacity function f1x and the maximum stiffness function f2x are established.(9)Fx=N1f1x+N2f2x,where N1 and N2 represent weighting factors which reflect the specific gravity of bearing capacity and stiffness in the total objective function, respectively. And, they are calculated by the equation:(10)Ni=1fiX∗,where fiX∗ represents the optimal solution of single constraint problem for the i-th subobjective as the objective function.


## 4.2. Function Determination of the Performance of Aerostatic Spindle

According to the assumption of uniform radial flow of gas, the motion equation in cylindrical coordinates can be simplified as follows:(11)∂p∂r=μ∂2vr∂z2,∂p∂θ=0,∂p∂z=0.By integrating equation (11) with z twice and combining the equivalent mass flow equation and the gas state equation, we can get the expression of pressure p as follows:(12)p12=pa2−pa2−ps2lnr/d4lnd2/d41/2,d4<r<d2,(13)p22=p02−p02−pa2lnr/d2lnd1/d21/2,d2<r<d1.The expression of bearing capacity can be obtained by substituting equations (12) and (13) in equation (6):(14)W=pa−psπd22−d422lnd2/d4+p0−paπd12−d222lnd1/d2−πd12−d42p0+πd32pa.Then, we can getf1x.(15)f1x=pa−psπx32−d422lnx3/d4+p0−paπd12−x322lnd1/x3−πd12−d42p0+πx22pa.And then, we define the gauge pressure ratio by using Powell’s method.(16)kg=pa−p0ps−p0.So, we write the bearing capacity as follows:(17)W=kgpa−psπd22−d422lnd2/d4+p0−paπd12−d222lnd1/d2−πd12−d42p0+πd32pa.Then, we introduce the slit factorG.(18)G=pa/ps1−pa/ps1/21+pa/ps⋅24μ2RT1/2pa⋅CDd32lnd1/d28h3,(19)kg=21+1+4/G21/2.According to the equation (7), we can get the objective function f2x of the stiffness as follows:(20)f2x=−0.98x1ps−paπ2x2x32−d42lnx3/d4,where x1 represents the gas film thickness h, x2 is the orifice diameter d3, and x3 is the diameter of distribution circle d2. d1 represents the external diameter of bearing, d4 represents the inner diameter of bearing, and ps represents the air supply pressure. p0 represents the atmospheric pressure and the discharge coefficient CD = 0.8.


## 5. Optimization Process

### 5.1. Optimization of Single Objective Function

Using equation (15) as the fitness function to write the M file program for optimizing the bearing capacity function. The constrains are also transformed to the program of M file, that is, the ranges of h, d3 and d2 are 5–20 μm, 2 mm–5 mm and 75–85 mm respectively. The M files reflecting the bearing capacity and constrains are imported to the genetic algorithm toolbox of MATLAB software. The process of solution is shown in Figure 4. The result of optimization is shown in Figure 5. When the design variable X = [8.301 3.213 82.531], that is, bearing clearance is 8.301 μm, orifice diameter is 3.213 mm, and the diameter of stomatic distribution circle is 82.531 mm, the optimal result of single objective function f1x is obtained. The optimal value of the bearing capacity is max f1x = 1254.2.Figure 4
Import of goal function and constrains.Figure 5
Optimization results of capacity.Use equation (20) as another fitness function to write the M file program for optimizing the stiffness function. The constrains are also transformed to the program of M files. The M files reflecting the stiffness and constrains are imported to the genetic algorithm toolbox. The process of solution is shown in Figure 6. The result of optimization is shown in Figure 7. When the design variable X = [12.342 2.912 82.125], that is, bearing clearance, stomatal diameter, and the diameter of stomatic distribution circle are 12.342 μm, 2.912 mm, and 82.125 mm respectively, the optimal result of single objective function f2x is obtained. The optimal value of the stiffness is max f2x = 135.4.Figure 6
Import of goal function and constrains.Figure 7
Optimization results of stiffness.The bearing capacity and stiffness of aerostatic spindle are optimized by using genetic algorithm toolbox in simulation system software that combines the optimization theory and genetic algorithm. First, the performance of aerostatic bearing is optimized by the single objective function of the capacity and the stiffness. The optimization results are shown in Table1.Table 1
Optimal results of single objective function.
Optimization values
Optimization parameters
h (μm)
d
3 (mm)
d
2 (mm)
W = 1254.2 N
8.301
3.213
82.531
K = 135.4 N/μm
12.342
2.912
82.125It can be seen from Table1 that the contradiction is present in the optimal result of the stiffness and load capacity function, respectively. Thus multiobjective optimization function is used to optimize the performance of the aerostatic spindle in order to get ideal consistent results.


### 5.2. Optimization of Multiobjective Function

The optimal values of bearing capacity and stiffness calculated above are brought into equation (10) to obtain the weighting factors N1 and N2. Substituting formulas (15) and (19) in (9), the multiobjective function including bearing capacity, stiffness, and influence parameters of aerostatic bearing is determined.(21)FX=0.00079f1X+0.00741f2X.The fitness function is shown in formula (12). The multiobjective function and constrains are imported to the genetic algorithm toolbox in the simulation. The optimization process of the solution is shown as Figure 8.Figure 8
Optimization process.The result of optimization is shown in Figure9. When the design variable X = [10.521 3.501 82.214], that is, bearing clearance is 10.521 μm, orifice diameter is 3.501 mm, and the diameter of stomatic distribution circle is 82.214 mm, the optimal result of multiobjective function FX is max FX = 1.8682.Figure 9
Optimization results of multiple objective.


## 5.1. Optimization of Single Objective Function

Using equation (15) as the fitness function to write the M file program for optimizing the bearing capacity function. The constrains are also transformed to the program of M file, that is, the ranges of h, d3 and d2 are 5–20 μm, 2 mm–5 mm and 75–85 mm respectively. The M files reflecting the bearing capacity and constrains are imported to the genetic algorithm toolbox of MATLAB software. The process of solution is shown in Figure 4. The result of optimization is shown in Figure 5. When the design variable X = [8.301 3.213 82.531], that is, bearing clearance is 8.301 μm, orifice diameter is 3.213 mm, and the diameter of stomatic distribution circle is 82.531 mm, the optimal result of single objective function f1x is obtained. The optimal value of the bearing capacity is max f1x = 1254.2.Figure 4
Import of goal function and constrains.Figure 5
Optimization results of capacity.Use equation (20) as another fitness function to write the M file program for optimizing the stiffness function. The constrains are also transformed to the program of M files. The M files reflecting the stiffness and constrains are imported to the genetic algorithm toolbox. The process of solution is shown in Figure 6. The result of optimization is shown in Figure 7. When the design variable X = [12.342 2.912 82.125], that is, bearing clearance, stomatal diameter, and the diameter of stomatic distribution circle are 12.342 μm, 2.912 mm, and 82.125 mm respectively, the optimal result of single objective function f2x is obtained. The optimal value of the stiffness is max f2x = 135.4.Figure 6
Import of goal function and constrains.Figure 7
Optimization results of stiffness.The bearing capacity and stiffness of aerostatic spindle are optimized by using genetic algorithm toolbox in simulation system software that combines the optimization theory and genetic algorithm. First, the performance of aerostatic bearing is optimized by the single objective function of the capacity and the stiffness. The optimization results are shown in Table1.Table 1
Optimal results of single objective function.
Optimization values
Optimization parameters
h (μm)
d
3 (mm)
d
2 (mm)
W = 1254.2 N
8.301
3.213
82.531
K = 135.4 N/μm
12.342
2.912
82.125It can be seen from Table1 that the contradiction is present in the optimal result of the stiffness and load capacity function, respectively. Thus multiobjective optimization function is used to optimize the performance of the aerostatic spindle in order to get ideal consistent results.


## 5.2. Optimization of Multiobjective Function

The optimal values of bearing capacity and stiffness calculated above are brought into equation (10) to obtain the weighting factors N1 and N2. Substituting formulas (15) and (19) in (9), the multiobjective function including bearing capacity, stiffness, and influence parameters of aerostatic bearing is determined.(21)FX=0.00079f1X+0.00741f2X.The fitness function is shown in formula (12). The multiobjective function and constrains are imported to the genetic algorithm toolbox in the simulation. The optimization process of the solution is shown as Figure 8.Figure 8
Optimization process.The result of optimization is shown in Figure9. When the design variable X = [10.521 3.501 82.214], that is, bearing clearance is 10.521 μm, orifice diameter is 3.501 mm, and the diameter of stomatic distribution circle is 82.214 mm, the optimal result of multiobjective function FX is max FX = 1.8682.Figure 9
Optimization results of multiple objective.


## 6. Discussion

### 6.1. Comparison of Static Performance by Optimization Parameters

The comparison of optimization parameters obtained by multiobjective functionFX and original parameters are shown in Table 2. Using optimization parameters and original parameters, the bearing capacity and stiffness of the different bearing clearance are simulated. The bearing performance before and after optimization is shown in Figures 10 and 11. Compared with original parameters, the performance of bearing with optimized parameters is improved. The bearing capacity and the stiffness are increased 7.98% and 9.17%, respectively.Table 2
Comparison of optimization parameters and original parameters.
h (μm)
d
3 (m)
d
2 (mm)
W (N)
K (N/μm)
Original parameters
12
3.5
82.5
1052
120
Optimization parameters
10.521
3.501
82.214
1136
131Figure 10
Comparison of bearing capacity.Figure 11
Comparison of bearing stiffness.The aerostatic spindle system and the static stiffness measurement system of aerostatic bearing are, respectively, established by the optimization parameters of the bearing, as shown in Figure12. An inductive micrometer is used to measure the variation of the displacement of the spindle table, and it is necessary to perform zero adjustment first. The specific operation is to determine the height of the induction microprobe by adjusting the height of the measuring gantry. The contact between the probe and the table is within ±2 μm to ensure the accuracy of the test.Figure 12
Measurement system of aerostatic bearing static stiffness.Air compressor provides gas to the aerostatic spindle system, and the supply pressure is 0.5 MPa. When the mass block as load is applied to worktable of aerostatic spindle, the workbench will have vertical displacement change. The change is sensed by the probe of the inductive micrometer and transmitted to the inductive micrometer which shows the change in display. The weight of mass block is 50 N, and the method of stepwise loading is used to perform multiple measurements and average. According to loaded load of the aerostatic bearing and the displacement change detected by the inductance micrometer, the static stiffness value is calculated by the stiffness calculation formula. The values of stiffness test are shown in Figure 13 which shows that the stiffness after optimization is more close to the values of stiffness test. So, the correctness and the feasibility of optimization method in this paper are verified by stiffness test.Figure 13
Comparison of bearing stiffness.


### 6.2. Comparison of Dynamic Response by Optimization Parameters

By optimizing the parameters, the static performance of the aerostatic spindle is improved; however, the dynamic performance of the spindle can be affected. This is because the corresponding frequency varies with the change of bearing capacity and stiffness. Here, the COMBIN14 spring-damper unit is used in modal analysis of the spindle model, in which the thrust bearings can be seen as the longitudinal spring-damper. The node of each orifice in the gas cavity is only stretched or compressed and not affected by bending and torsion. The radial bearings can be seen as the longitudinal and torsional spring-damper, and the node of each orifice in the gas cavity has the rotary element with three degrees of freedom in addition to the tension-compression in the normal direction. The model of aerostatic spindle is shown in Figure14.Figure 14
Model of aerostatic spindle.And the natural frequencies of the first six orders before and after optimization are calculated and shown in Table3, in which we can get that the corresponding frequency has been increased 5% after optimization.Table 3
First 6-order frequency of spindle in two cases.
Orders
1
2
3
4
5
6
Original parameters (Hz)
494
496
1439
4857
4876
5591
Optimization (Hz)
516
518
1503
5075
5094
5842The dynamic characteristics of the aerostatic spindle are important in actual work, and the performance in axial direction is the key factor affecting the machining accuracy of the machine tool. In order to better verify the optimization effect of aerostatic bearing, the axial dynamic response is analyzed. The motion equation which has natural frequency and damping coefficient is proposed.(22)x¨+2ωnξx˙+ωn2x=0,(23)ξ=C2mωn.where ωn is the natural frequency, m is the mass of the spindle, which is 3.7 kg, and C is the damping coefficient. The solution of above equation is as shown.(24)x0t=1−e−ξωnt1−ξ2sinωdt+arctan1−ξ2ξ.And the response curve is shown in Figure15 which indicates that the response speed after optimization has been improved, but it still in the same order of magnitude with the values before optimization.Figure 15
Dynamic response of aerostatic spindle.


## 6.1. Comparison of Static Performance by Optimization Parameters

The comparison of optimization parameters obtained by multiobjective functionFX and original parameters are shown in Table 2. Using optimization parameters and original parameters, the bearing capacity and stiffness of the different bearing clearance are simulated. The bearing performance before and after optimization is shown in Figures 10 and 11. Compared with original parameters, the performance of bearing with optimized parameters is improved. The bearing capacity and the stiffness are increased 7.98% and 9.17%, respectively.Table 2
Comparison of optimization parameters and original parameters.
h (μm)
d
3 (m)
d
2 (mm)
W (N)
K (N/μm)
Original parameters
12
3.5
82.5
1052
120
Optimization parameters
10.521
3.501
82.214
1136
131Figure 10
Comparison of bearing capacity.Figure 11
Comparison of bearing stiffness.The aerostatic spindle system and the static stiffness measurement system of aerostatic bearing are, respectively, established by the optimization parameters of the bearing, as shown in Figure12. An inductive micrometer is used to measure the variation of the displacement of the spindle table, and it is necessary to perform zero adjustment first. The specific operation is to determine the height of the induction microprobe by adjusting the height of the measuring gantry. The contact between the probe and the table is within ±2 μm to ensure the accuracy of the test.Figure 12
Measurement system of aerostatic bearing static stiffness.Air compressor provides gas to the aerostatic spindle system, and the supply pressure is 0.5 MPa. When the mass block as load is applied to worktable of aerostatic spindle, the workbench will have vertical displacement change. The change is sensed by the probe of the inductive micrometer and transmitted to the inductive micrometer which shows the change in display. The weight of mass block is 50 N, and the method of stepwise loading is used to perform multiple measurements and average. According to loaded load of the aerostatic bearing and the displacement change detected by the inductance micrometer, the static stiffness value is calculated by the stiffness calculation formula. The values of stiffness test are shown in Figure 13 which shows that the stiffness after optimization is more close to the values of stiffness test. So, the correctness and the feasibility of optimization method in this paper are verified by stiffness test.Figure 13
Comparison of bearing stiffness.


## 6.2. Comparison of Dynamic Response by Optimization Parameters

By optimizing the parameters, the static performance of the aerostatic spindle is improved; however, the dynamic performance of the spindle can be affected. This is because the corresponding frequency varies with the change of bearing capacity and stiffness. Here, the COMBIN14 spring-damper unit is used in modal analysis of the spindle model, in which the thrust bearings can be seen as the longitudinal spring-damper. The node of each orifice in the gas cavity is only stretched or compressed and not affected by bending and torsion. The radial bearings can be seen as the longitudinal and torsional spring-damper, and the node of each orifice in the gas cavity has the rotary element with three degrees of freedom in addition to the tension-compression in the normal direction. The model of aerostatic spindle is shown in Figure14.Figure 14
Model of aerostatic spindle.And the natural frequencies of the first six orders before and after optimization are calculated and shown in Table3, in which we can get that the corresponding frequency has been increased 5% after optimization.Table 3
First 6-order frequency of spindle in two cases.
Orders
1
2
3
4
5
6
Original parameters (Hz)
494
496
1439
4857
4876
5591
Optimization (Hz)
516
518
1503
5075
5094
5842The dynamic characteristics of the aerostatic spindle are important in actual work, and the performance in axial direction is the key factor affecting the machining accuracy of the machine tool. In order to better verify the optimization effect of aerostatic bearing, the axial dynamic response is analyzed. The motion equation which has natural frequency and damping coefficient is proposed.(22)x¨+2ωnξx˙+ωn2x=0,(23)ξ=C2mωn.where ωn is the natural frequency, m is the mass of the spindle, which is 3.7 kg, and C is the damping coefficient. The solution of above equation is as shown.(24)x0t=1−e−ξωnt1−ξ2sinωdt+arctan1−ξ2ξ.And the response curve is shown in Figure15 which indicates that the response speed after optimization has been improved, but it still in the same order of magnitude with the values before optimization.Figure 15
Dynamic response of aerostatic spindle.


## 7. Conclusions

In this paper, the performance of aerostatic spindle has been improved by a novel optimization approach, and the results of before and after optimization are compared in static and dynamic aspect. In the optimization process, the main factors affecting the bearing capacity and stiffness of aerostatic bearings are determined by cross-correlation analysis, and the corresponding actual experimental platform of the aerostatic spindle system is established for analysis. The following major conclusions are obtained:(1)
A new optimization method for the performance of aerostatic spindle considering the main influence factors was proposed and proved experimentally. Taking the load capacity and stiffness of the spindle as the optimize function in the multiobjective function. The relationship between the optimize function and the main influence factors is proposed. Then, the maximum values are obtained by genetic algorithm.(2)
The influence factors affecting the performance of the aerostatic spindle were evaluated, whose influence degree were obtained by cross-correlation analysis. And the correlation coefficients of more than 0.7 were determined, which are imported into the optimization process as the independent variables.(3)
The dynamic characteristics after optimization were analyzed, which show that the optimization method proposed in the paper provided effective parameters for the aerostatic spindle. The corresponding natural frequency and the response speed have been increased when the response amplitude is in the same order magnitude with the value before optimization. The optimization of aerostatic bearing is proceeded combining the theory with microscale characteristics. The performance of aerostatic bearing is improved.



---

*Source: 8124029-2019-04-15.xml*