Shock and Vibration
Volume 2018 (2018), Article ID 8209290, 16 pages
https://doi.org/10.1155/2018/8209290

Research Article

An Analytical Model for a Pneumatic Vibration Isolator with the Stiffness Effect of the Elastomeric Diaphragm

Dengfeng Xu,1,2 Qiang Yu,³ Fei Shen,³ Yu Zhu,1,2,3 and Gaofeng Guan³

¹State Key Laboratory of Tribology, Department of Mechanical Engineering, Tsinghua University, Beijing 100084, China
²Beijing Lab of Precision/Ultra-Precision Manufacture Equipment and Control, Tsinghua University, Beijing 100084, China
³School of Mechanical and Electrical Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China

Correspondence should be addressed to Dengfeng Xu; xudf02@mails.tsinghua.edu.cn

Received 4 April 2018; Accepted 20 May 2018; Published 5 July 2018

Academic Editor: Pedro Museros

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

Abstract.

The elastomeric diaphragm is widely used in pneumatic vibration isolators, and the relevant model is often ignored due to its complexity. Considering that the ignored model of the elastomeric diaphragm in pneumatic vibration isolators plays an important role in the discrepancy between the predicted and practical behavior, this paper develops an analytical model for the elastomeric diaphragm using the Mooney-Rivlin modeling method and elastomeric theory. Specifically, the elastomeric diaphragm consists of several segments in the axial section. After considering the structural restriction, each segment can be simplified as uniaxial stretching, and the force-strain equation can be established for each segment. By combining the equations of all segments, an analytical model of the elastomeric diaphragm can be built and solved via numerical methods. The developed model is added to the standard model of pneumatic cylinders to supply a complete analytical model for pneumatic vibration isolators. The experimental results demonstrate that the analytically predicted behavior is similar to the practical behavior. The proposed analytical model can be used as a guide for the parameter design of pneumatic isolators in practice.

1. Introduction

Pneumatic vibration isolators are widely applied in optics, semiconductors, and other fields [1, 2]. Many efforts have been devoted to the design of active control systems, which require an accurate mathematical model of the pneumatic vibration isolator [3, 4]. Harris [5] assumed that the airflow in the damping hole was laminar and first offered a linear model for the pneumatic vibration isolators. In 1984, Debra and Bryan [6] established a simpler linear model of the volume ratio of the upper chamber to the lower chamber and the inherent frequency. Erin and Wilson [7] proposed an analytical modeling method for the pneumatic vibration isolator based on the ideal gas state equation. In all previous studies, the discrepancies between the model-predicted behavior and the practically observed behavior of pneumatic vibration isolators exceeded 50% because the effect of the elastomeric diaphragm was ignored.

To obtain a more accurate model of the pneumatic isolator, Erin and Wilson proposed an improved model for pneumatic isolators by incorporating the stiffness and damping effect of the elastomeric diaphragm based on experimental investigation. Lee and Kim [8, 9] reported a nonlinear complex stiffness model that reflected the frequency and excitation amplitude, which was validated by comparison with experimental measurements. Pu et al. [10] obtained the stiffness and damping characteristics of the elastomeric diaphragm indirectly by comparing the experimental and theoretical results. Additionally, Pu established the stiffness and damping model using a similar approach in other research. Since the 1990s, analog simulation has become an attractive method for examining the pneumatic spring, and the stiffness characteristic is often studied using the finite element method (FEM) [11]. However, experimental modification or the FEM is inconvenient for the parameter optimization design of a pneumatic isolator. Significant error typically still exists between the actual behavior and the presented model after the experimental modification or FEM study. Therefore, a direct theoretical derivation method must be established for the development of an analytical model for a pneumatic vibration isolator.

Zhu et al. [12] and Li [13] proposed an analytic model of an elastomeric belt in a magnetic-levitation vibration isolator using the Mooney-Rivlin model, and the model was validated experimentally. In that research, the elastomeric belt was chosen as the research object, and an accurate model was obtained according to direct theoretical derivation. Inspired by this work, we use the same method and theory to study and analyze the elastomeric diaphragm in pneumatic vibration absorbers and develop an analytical model of the elastomeric diaphragm in this paper. Specifically, because the elastomeric diaphragm is composed of several segments in the axial section, each segment can be simplified as uniaxial stretching after considering the structural restriction, and a force-strain equation can be established for each segment. The analytical model of the elastomeric diaphragm can be built by combining all segments, and the model can be resolved by numerical solution. Experimental tests were conducted, and the results demonstrate that the analytically predicted behavior is similar to the practically observed behavior. The error in this method can be reduced to less than 10% after the elastomeric diaphragm model is added to the standard pneumatic cylinder model. Although certain scholars have proposed several improved models that also imitate the practical behavior, this paper supplies a direct theoretical derivation method with a notably small error. The proposed analytical model facilitates the practical parameter optimization design of pneumatic isolators.

2. Mathematic Model of the Elastomeric Diaphragm

2.1. Configuration of the Pneumatic Vibration Isolator

The configuration of a typical pneumatic vibration isolator is shown in Figure 1(a), which includes an air-filled cylinder, elastomeric diaphragm, and piston. The cylinder consists of a top chamber and bottom chamber that are both filled with compressed air. The elastomeric diaphragm is located between the piston and cylinder, which is the seal component of the cylinder. The payload is connected to the piston, and the flow restrictor supplies damping for the vibration isolator.

(a) Schematic of a pneumatic vibration isolator

(b) Simplification of a pneumatic vibration isolator

Figure 1: Pneumatic vibration isolator.

Figure 1(b) illustrates that the payload is supported by two components: the pressurized air in the pneumatic cylinder and the elastomeric diaphragm, the model of which is often ignored. This paper only analyzes the stiffness effect of the elastomeric diaphragm without considering the damping effect because the vibration damping effect of the elastomeric diaphragm is a highly complex phenomenon. Therefore, the total stiffness of the pneumatic spring isolators can be obtained by adding the vertical stiffness of the elastomeric diaphragm (the model typically ignored in previous analysis and research) to the vertical stiffness of the pressurized air. The total stiffness of the pneumatic spring can be expressed as follows [14]:

where represents the total stiffness of the pneumatic spring isolators, represents the vertical stiffness of the compressed air, and represents the vertical stiffness of the elastomeric diaphragm. has typically been neglected in previous analyses of the pneumatic spring, which causes a certain error between the theory and experiments. This error in the total vertical stiffness can result in a series of problems. For example, the actual inherent frequency of a pneumatic vibration absorber cannot be easily estimated with high accuracy. After considering the effect of elastomeric diaphragm, a complete analytical model for a pneumatic vibration isolator is developed, and the developed analytical model can be used in the parameter optimization design of a pneumatic vibration isolator when the total stiffness and inherent frequency of the pneumatic vibration absorber are both known.

Figure 2 shows a cross-sectional view of an elastomeric diaphragm, and the parameters of the elastomeric diaphragm are listed in Table 1. After the elastomeric diaphragm is divided into four segments, the vertical stiffness of the segments, expressed as , , , and , is deduced such that the total vertical stiffness of the elastomeric diaphragm can be determined.

Table 1: Parameters of the elastomeric diaphragm.

Symbol
Description

Radius of the piston (radius of the inner circle of the elastomeric diaphragm)

Radius of the arc of the elastomeric diaphragm

Height of the inner circle of the elastomeric diaphragm

Height of the outer circle

Thickness of the elastomeric diaphragm

“1”
Inner circle of the elastomeric diaphragm

“2”
Inner arc of the elastomeric diaphragm

“3”
Outer circle of the elastomeric diaphragm

“4”
Outer arc of the elastomeric diaphragm

(a) Size definition of the elastomeric diaphragm

(b) Parameters of the elastomeric diaphragm

Figure 2: Elastomeric diaphragm.

2.2. Force Analysis of the Elastomeric Diaphragm

The forces of the piston and elastomeric diaphragm in the vertical direction are shown in Figure 3; then,

where is the pulling force of part “2” acting on part “1” and is the pulling force of part “4” acting on part “3”. In addition, and are the reaction forces of and , respectively, and Fa, and can be rewritten as

where is the vertical force of the air in the chamber acting on the effective piston area, is the vertical force of the air outside the chamber acting on the effective piston area, is the gravity of the load, is the atmospheric pressure, is the air pressure in the chamber, is the payload mass, is the gravitational acceleration, and is the corresponding stretch ratio of a unit of part “1”, which is explained later. The stretch ratios in this paper are defined as the ratios of the stretched length to the unstretched length or original length of a small-volume element, which is a cube when it is unstretched or in the original state.

(a) Force analysis of the piston

(b) Force analysis of the elastomeric diaphragm

Figure 3: Force analysis of the piston and elastomeric diaphragm.

As shown in Figure 3(b), is the force acting on the arc of the elastomeric diaphragm in the vertical direction and is expressed as

where is the projection area of the arc of the elastomeric diaphragm, as shown in Figure 4. Then, can be expressed as

where is the corresponding stretch ratio of a unit of part “3”, which is explained later.

Figure 4: Cross-section of the dual-chamber pneumatic spring.

Substituting (8) into (7) leads to

By substituting (9) into (3) and combining the result with (2), (4), (5), and (6), the air pressure in the chamber can be obtained as

By substituting (10) into (9) and noting (3), one can obtain

2.3. Mooney-Rivlin Model

Mooney-Rivlin [15, 16] assumed that the elastic behavior was isotropic and that simple stretching satisfied Hooke’s law; then, they obtained the rubber energy density function under the premise of not changing the rubber state and compressing its volume; i.e.,

and

where , , and represent the stretch ratios in three different directions and are defined as the ratios of the stretched length to the unstretched length of a small-volume element, which is a cube when it is unstretched. These three directions, which are related to , , and , are expressed under generalized XYZ coordinates. If is defined as a corresponding stretched length, the stretch ratio can be expressed as =1+/1. Herein, the rubber satisfies the incompressibility condition, i.e., =1, which means that =0. The corresponding deformation energy density only depends on and . Under the condition of sufficiently small deformation, the form of deformation energy function can be expressed as

where and represent different material parameters.

2.4. Stiffness Analysis of the Elastomeric Diaphragm

2.4.1. Part “1” Stiffness of the Elastomeric Diaphragm

As shown in Figure 5(a), when an elastomeric diaphragm is stretched, the tension acting on part “1” can be expressed as

where is the sectional area of part “1”, which is shown in Figure 4, and can be expressed as

where is the corresponding stretch ratio of a unit of part “1”. In this paper, the letters of the subscript contain two numbers i and j, where i(1,2,3) represents the corresponding principal stress surface and j(1,2,3,4) represents the corresponding component of the elastomeric diaphragm, i.e., the inner circle of the elastomeric diaphragm, the inner arc of the elastomeric diaphragm, the outer circle of the elastomeric diaphragm, and the outer arc of the elastomeric diaphragm, respectively.

(a) Schematic diagram and force analysis of part “1”

(b) Parameters and force analysis of “1” unit of part “1”

Figure 5: Force analysis of part “1” of the elastomeric diaphragm.

Figure 5(b) shows the unit of part “1”, and the unit stress in three directions can be expressed as

where is the corresponding stretch ratio and is the unknown pressure. The stretch ratio can be expressed as =(+)/, where is the length variable in the vertical direction of part “1”.

The following equation can be obtained by combining (18) with (20):

Substituting (17) and (21) into (16) yields

The force on the unit in the cross-section is shown in Figure 6. The force of the piston acting on the unit of part “1” can be expressed as

where is the third principal stress surface of part “1” subjected to the force of the gas chamber, F₂₁ is the second principal stress surface of part “1” subjected to the force, and is the angle of the unit in the cross-section of the dual-chamber pneumatic spring.

Figure 6: Cross-section of one unit of part “1”.

and are defined as

where is the area of the third principal stress surface of part “1”, is the stress of the second principal stress surface of part “1”, and is the area of the second principal stress surface of part “1”, as shown in Figure 5(b).

Substituting (24) and (25) into (23) yields

which indicates that . Therefore, part “1” of the elastomeric diaphragm remains in contact with the piston wall, and the size of part “1” of the elastomeric diaphragm does not change when it is stretched in the circumferential direction; i.e., =1. The incompressibility of rubber leads to =1, and the following expression can be obtained:

Substituting (10) and (27) into (22) yields

To obtain a function for the stiffness of part “1” with respect to variable , the following equation can be obtained by differentiating (28) with respect to :

The complete equation is extremely long and is provided in Appendix B.

2.4.2. Part “2” Stiffness of the Elastomeric Diaphragm

As shown in Figure 7, the tension acting on part “2” can be expressed as

where originates from the pressure difference of the top chamber and ambient environment in the vertical direction.

(a) Schematic diagram of part “2”

(b) Parameters of “2” unit of part “2”

(c) Force analysis of “2” unit of part “2”

(d) Force analysis of part “2”

Figure 7: Force analysis of part “2” of the elastomeric diaphragm.

For the projected area of the unit, can be written as

where is the length from the analysis point to the center axis, is the side length of the unit in part “2”, and is the thickness of the unit in part “2”. is shown in Figure 7(a).

The vertical stress of the unit of part “2” can be expressed as

where and can be written as

where and are the stresses of the corresponding principal stress surface and is the unknown pressure. The stretch ratio can be expressed as =(+)/, where is the length variable in the vertical direction. The size of part “2” does not change when it is stretched in the circumferential direction; i.e., =1. The incompressibility of rubber leads to =1. Thus, the following equation can be obtained:

The following equation can be obtained by transferring (32) to (35):

Substituting (10), (30), (31), and (36) into (30) yields

To obtain a function of the stiffness of part “2” with respect to variable , the following equation can be obtained by differentiating (37) with respect to :

The complete equation is extremely long and is provided in Appendix B.

2.4.3. Part “3” Stiffness of the Elastomeric Diaphragm

As shown in Figure 8, when an elastomeric diaphragm is stretched, the tension acting on part “3” can be expressed as

(a) Schematic diagram and force analysis of part “3”

(b) Parameters and force analysis of “3” unit of part “3”

Figure 8: Force analysis of part “3” of the elastomeric diaphragm.

The analysis of part “3” is similar to that of part “1”. Therefore, based on the Mooney-Rivlin model, the vertical stress can be expressed as

where the stretch ratio can be expressed as =(+)/ and is the variable length in the vertical stretch of part “3”.

Additionally, is the cross-sectional area of part “3”, as shown in Figure 4. Herein, can be expressed as

Substituting (40) and (41) into (39) yields

The force on the unit in the cross-section is shown in Figure 9. The force of the piston acting on the unit of part “3” can be expressed as

where is the third principal stress surface of part “3” subjected to the force of the gas chamber, is the second principal stress surface of part “3” subjected to the force, and is the angle of the unit in cross-section of the dual-chamber pneumatic spring.

Figure 9: Cross-section of one unit of part “3”.

and are written as

where is the area of the third principal stress surface of part “3”, is the stress of the second principal stress surface of part “3”, and is the area of the second principal stress surface of part “3”, as shown in Figure 8(b).

Substituting (44) and (45) into (43) yields

where and is the arc length of the unit of part “3”, which is smaller than . Thus, the first polynomial is larger than the second polynomial of (46); i.e., . Therefore, part “3” remains in contact with the fixed ring, and the size of part “3” does not change when it is stretched in the circumferential direction; i.e., =1. However, the incompressibility of rubber leads to =1. Thus, the following equation can be obtained:

Substituting (10) and (47) into (42) yields

To obtain a function of the stiffness of part “3” with respect to variable , the following equation can be obtained by differentiating (48) with respect to :

The complete equation is extremely long and is provided in Appendix B.

2.4.4. Part “4” Stiffness of the Elastomeric Diaphragm

As shown in Figure 10, the tension acting on part “4” can be expressed as

where originates from the pressure difference of the top chamber and the ambient environment in the vertical direction.

(a) Schematic diagram of part “4”

(b) Parameters of the “4” unit of part “4”

(c) Force analysis of the “4” unit of part “4”

(d) Force analysis of part “4”

Figure 10: Force analysis of part “4” of the elastomeric diaphragm.

For the projected area of the unit, can be written as

where is the length from the analysis point to the center axis, is the side length of the unit in part “4”, is the thickness of the unit in part “4”, and is shown in Figure 10(a).

The vertical stress δ₄ of the partial unit in part “4” can be expressed as

and δ₁₄ and δ₃₄ can be rewritten as

where and are the stresses corresponding to the principal stress surface and is the unknown pressure. The stretch ratio can be expressed as =(+)/, where is the variable length in the vertical stretch of part “4”. The size of part “4” does not change when it is stretched in the circumferential direction; i.e., =1. The incompressibility of rubber leads to =1. Thus, the following equation can be obtained:

The following equation can be obtained by transferring (52) to (55):

Substituting (10), (51), (55), and (56) into (50) yields

To obtain a function of the stiffness of part “4” with respect to variable , the following equation can be obtained by differentiating (57) with respect to :

The entire equation is extremely long and is provided in Appendix B.

2.4.5. Total Stiffness of the Elastomeric Diaphragm

As shown in Figure 11, part “1” and part “2” are connected in series such that the equivalent stiffness can be expressed as

Figure 11: Total stiffness of the elastomeric diaphragm.

Part “3” and part “4” are connected in series such that the equivalent stiffness can be expressed as

The total stiffness can be expressed as

2.5. Numerical Solution of the Vertical Stiffness for the Elastomeric Diaphragm

When the load of the dual-chamber pneumatic spring is , each part of the elastomeric diaphragm is stretched to a certain length. The stretched lengths of each corresponding part are defined as , , , and . By analyzing (3) and (11), each component of the force can be expressed as

The analysis process to obtain the stretched length of each corresponding part with respect to the other variables is as follows. Substituting = into (28) and noting (62), we can obtain (B.5) of the stretched length . (B.6), (B.7), and (B.8) can also be obtained in a similar manner. Equations (B.5) to (B.8) are rather long and are provided in Appendix B.

By solving (B.5) to (B.8) in the numerical root of , , , and and substituting =, =, =, and = into , , , and , we obtain the numerical solutions of stiffness K₁, K₂, K₃, and K₄, respectively. As a result, the total stiffness can be obtained from (61).

3. Experimental Validation

To verify the theoretical analysis, we test the performance and calculate the stiffness of the specific dual-chamber pneumatic spring.

The following text describes the method for obtaining 2(C₁+C₂), where C₁ and C₂ are the Mooney-Rivlin material parameters. Classical rubber theory yields the relationship between the stiffness values of different types of rubbers and their Young’s elastic moduli [15]. The elastomeric diaphragms commonly used in design of the pneumatic spring structure consist of chloroprene rubber, nitrile butadiene rubber, butyl rubber, natural rubber, and other materials. A moderate stiffness is typically used to obtain a better overall performance. Therefore, the widely applied nitrile butadiene rubber with hardness (55±5) is used in the next two experiments. The Young’s elastic modulus changes from 2.465 MPa to 3.619 MPa via the relationship between the hardness and Young’s elastic modulus. According to the incompressibility of the rubber material, we can obtain . Hence, ranges from 0.82 MPa to 1.20 MPa, where is the shear modulus and =2(C₁+C₂). Finally, 2(C₁+C₂) also ranges from 0.82 MPa to 1.20 MPa.

The vertical total stiffness of the dual-chamber pneumatic isolator can be calculated according to the parameters in Table 2, where the stiffness of the gas is obtained based on Erin’s study. Two experiments are shown in Figures 12 and 13. The main components of the isolation platforms are composed of four dual-chamber pneumatic isolators. The microvibration of the ground can be attenuated by the dual-chamber pneumatic isolators. Thus, we can obtain the transmissibility curve using the testing instrument HANDYSAM [17]. The input and output sensors used in this study are the ultralow frequency vibration pickup sensor 991B produced by the Institute of Engineering Mechanics, China Earthquake Administration. The sensitivity of the sensors is 4×10⁻⁸ m/s. The load in experiment 1 is 50 kg, and the load in experiment 2 is 200 kg. Ground excitation used in the two experiments is just ambient vibration, whose signal in vertical direction measured by the input sensors is shown in Figure 14.

Table 2: Parameters of the dual-chamber pneumatic in the two experiments.

Symbol
Description
Value 1
Value 2
Unit

and are the Mooney-Rivlin material parameters
0.82-1.20
0.82-1.20
MPa

Part “1” vertical length
0.009
0.015
m

Part “2” or “4” projection length
0.0034
0.0050
m

Part “3” vertical length
0.014
0.020
m

Diaphragm thickness
0.0006
0.002
m

Atmospheric pressure
101325
101325
MPa

Gravitational acceleration
9.81
9.81
m/s²

Payload mass
20
375
kg

Polytropic exponent
1.4
1.4
-

Flow restrictor length
4e-4
8e-3
m

Flow restrictor diameter
1e-3
1e-3
m

Effective piston area
0.0015
0.016
m²

Air viscosity
1.824e-5
1.824e-5
Pa⋅s

Flow restriction constant
3.3640e-06
1.6019e-07
m³/(s⋅Pa)

Lower chamber volume
3.0000e-04
1.8170e-04
m³

Upper chamber volume
2.7298e-05
2.3000e-03
m³

Figure 12: Photograph of experiment 1.

Figure 13: Photograph of experiment 2.

Figure 14: Ground vibration velocity curve in vertical direction versus time.

The theoretical expression of the frequency response functions has been inferred without considering the elastomeric diaphragm [10]. In this paper, the theoretical transmissibility curves are obtained as shown in Figures 15 and 16 by substituting the actual parameters of the experimental pneumatic vibration absorber shown in Table 2 into Erin’s analytical model. The first curve is Erin’s result without consideration of the elastomeric diaphragm. The second and third curves are our results with consideration of the elastomeric diaphragm. The last curve is the experimental result with the parameters shown in Table 2. Our model yields two curves, namely, the upper and lower bounds for the inherent frequencies, which correspond to the maximum and minimum values of 2(C₁+C₂).

Figure 15: Transmissibility curves of experiment 1.

Figure 16: Transmissibility curves of experiment 2.

The results of experiment 1 are shown in Figure 15, where =1.89 Hz, =6.03-7.70 Hz, and =6.75 Hz are the inherent frequencies of the double-chamber pneumatic springs based on Erin’s study and the proposed theory in this paper that considers the elastomeric diaphragm stiffness and experiment results, respectively. The results of experiment 2 are shown in Figure 16, where =1.31 Hz, =3.38-4.30 Hz, and =3.75 Hz are the inherent frequencies based on Erin’s study and the proposed theoretical and experiment results, respectively. The value of the damping in our two experiments is related to both the length and diameter of the flow restrictor. In experiment 1, the flow restrictor length is 0.0004, and the flow restrictor diameter is 0.001; thus, the value of the damping ratio is 0.005. In experiment 2, the flow restrictor length is 0.008, and the flow restrictor diameter is 0.016; thus, the value of the damping ratio is 0.1.

Figures 15 and 16 illustrate that and are within the ranges of and , respectively, demonstrating the effectiveness of the vertical stiffness model for the elastomeric diaphragm.

The stiffness values of the experiments can be obtained from the relationship of the inherent frequency, the stiffness, and the payload; i.e.,

The calculation results of the two structures with different parameters of the dual-chamber pneumatic spring are shown in Table 3. The vertical stiffness, KE1 and KE2, in the experimental tests is within the corresponding ranges of total stiffness Ktotal obtained via calculations, demonstrating the effectiveness of the vertical stiffness model for the elastomeric diaphragm.

Table 3: Stiffnesses from the calculations and experiments.

Symbol
Description
Value 1
Value 2
Unit

Part “1” length stretched
15.0e-04 to 9.19e-04
22.0e-04 to 14.0e-04
m

Part “2” length stretched
8.53e-04 to 5.16e-04
9.90e-04 to 5.85e-04
m

Part “3” length stretched
13.0e-04 to 8.55e-04
18.0e-04 to 11.0e-04
m

Part “4” length stretched
6.26e-04 to 3.88e-04
8.23e-04 to 4.94e-04
m

Part “1” stiffness
1.87e+04 to 3.24e+04
9.97e+04 to 1.75e+05
N/m

Part “2” stiffness
3.54e+04 to 6.32e+04
2.34e+05 to 4.12e+05
N/m

Part “3” stiffness
1.90e+04 to 3.10e+04
1.04e+05 to 1.73e+05
N/m

Part “4” stiffness
4.84e+04 to 8.42e+04
2.52e+05 to 4.54e+05
N/m

Part “1” and “2” equivalent stiffness
1.23e+04 to 2.14e+04
6.99e+04 to 1.23e+05
N/m

Part “3” and “4” equivalent stiffness
1.36e+04 to 2.27e+04
7.35e+04 to 1.25e+05
N/m

Rubber diaphragm stiffness
2.59e+04 to 4.41e+04
1.43e+05 to 2.48e+05
N/m

Air stiffness
0.2850e+04
2.5406e+04
N/m

Total stiffness
2.87e+04 to 4.90e+04
1.69e+05 to 2.74e+05
N/m

Experimental stiffness
3.5975e+04
2.0819e+05
N/m

Calculated inherent frequency
6.0332 to 7.7096
3.3769 to 4.2931
Hz

Experimental inherent frequency
6.625
3.5
Hz

4. Conclusions

This work developed an analytical model of the elastomeric diaphragm using the Mooney-Rivlin modeling method and elastomeric theory. The stiffness of the elastomeric diaphragm was ignored in the previous model of pneumatic vibration isolators. Specifically, each segmented elastomeric diaphragm can be simplified as uniaxial stretching, and the Mooney-Rivlin model can be adopted to establish the analytical expression of the vertical stiffness. The analytical model of the elastomeric diaphragm is built by combining all segments and can be resolved by numerical solution. The experimental results demonstrated that the analytically predicted behavior is highly similar to the practically observed behavior in terms of the resonant frequency of the isolator after the elastomeric diaphragm model is added to the standard pneumatic cylinder model. Differing from the prior converse correction and experience design, this paper supplies a direct design method that is helpful for optimizing and designing the pneumatic spring parameters in practice.

Appendix

A.

See Table 4.

Table 4: All parameter descriptions in this paper.

Sym
Description
Sym
Description

Total stiffness of the pneumatic spring

Vertical stress of the unit of part “2”

Vertical stiffness of the air

Vertical stress of the unit of part “4”

Vertical stiffness of the elastomeric diaphragm

Corresponding stretch ratio representing the corresponding principal stress surface

Radius of the piston

Corresponding stress of the surface

Radius of the arc of the elastomeric diaphragm

Corresponding principal stress surface (i = 1, 2, 3)

Height of the inner circle of the elastomeric diaphragm

Corresponding part of the elastomeric diaphragm (j = 1, 2, 3, 4)

Height of the outer circle

Sectional area of part “1”

Thickness of the elastomeric diaphragm

Unknown pressure

“1”
Inner circle of the elastomeric diaphragm

Force of piston acting on the unit of “1”

“2”
Inner arc of the elastomeric diaphragm

Third principal stress surface of “1” subjected to the force of the gas chamber

“3”
Outer circle of the elastomeric diaphragm

Second principal stress surface of “1” subjected to the force of the gas chamber

“4”
Outer arc of the elastomeric diaphragm

Angle of the unit of “1” in the cross-section of the dual-chamber pneumatic spring

Pulling force of “2” acting on “1”

Tension acting on part“2”

Pulling force of “4” acting on “3”

Sectional area of the unit of “1”

Reaction forces of

Length from the unit of “2” to the center axis

Reaction forces of

Side length of the unit of “2”

Gravity of the load

Angle

Atmospheric pressure

Angle

Air pressure in the chamber

Unknown pressure

Payload mass

Sectional area of part “3”

Gravitational acceleration

Force of piston acting on the unit of “3”

Force acting on the arc of the elastomeric diaphragm

Third principal stress surface of “3” subjected to the force of the gas chamber

Projection area of the arc of the elastomeric diaphragm

Second principal stress surface of “3” subjected to the force of the gas chamber

Rubber energy density function

Angle of the unit of “1” in the cross-section of the dual-chamber pneumatic spring

First strain invariant

Tension acting on part“4”

Second strain invariant

Sectional area of the unit of “4”

Third strain invariant

Length from the unit of “4” to the center axis

Stretch ratios in the X direction

Side length of the unit of “4”

Stretch ratios in the Y direction

Angle

Stretch ratios in the Z direction

Unknown pressure

Material parameter

Stiffness of part “1”

Material parameter

Stiffness of part “2”

Length variable of part “1”

Stiffness of part “3”

Length variable of part “2”

Stiffness of part “4”

Length variable of part “3”

Series stiffness of part “1” and “2”

Length variable of part “4”

Series stiffness of part “3” and “4”

B. Selected Complex Equations in This Paper

Data Availability

The “Experiment 1.rar” and “Experiment 2.rar” compressed files used to support the findings of this study have been deposited in the Figshare repository (https://figshare.com/s/86ee3484e7f439547b71). Open the compressed files and then you can find Erin result, experimental result, experimental data, lower bound for the inherent frequency, and upper bound for the inherent frequency, respectively. You can use MATLAB software to open “.m” and “.fig” type files and use Microsoft Excel software to open “.xlsx” type files.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grant 51475262 and the National Science and Technology Major Project 2015ZX02104003.

References

J.-W. Liang, H.-Y. Chen, and Q.-W. Wu, “Active suppression of pneumatic vibration isolators using adaptive sliding controller with self-tuning fuzzy compensation,” Journal of Vibration and Control, vol. 21, no. 2, pp. 246–259, 2013.

L.-E. Ooi and Z. M. Ripin, “Optimization of an engine mounting system with consideration of frequency-dependent stiffness and loss factor,” Journal of Vibration and Control, vol. 22, no. 10, pp. 2406–2419, 2014.

F. C. Nelson, “Vibration isolation: a review, I. Sinusoidal and random excitations,” Shock and Vibration, vol. 1, no. 5, pp. 485–493, 1994.

F. C. Nelson, “Vibration Isolation Review: IL Shock Excitation,” Shock and Vibration, vol. 3, no. 6, pp. 451–459, 1996.

C. M. Harris, Harris’ Shock and Vibration Handbook, McGraw-Hill, New York, NY, USA, 2002.

D. B. DeBra and J. Bryan, “Design of laminar flow restrictors for damping pneumatic vibration isolators,” CIRP Annals - Manufacturing Technology, vol. 33, no. 1, pp. 351–356, 1984.

C. Erin, B. Wilson, and J. Zapfe, “An improved model of a pneumatic vibration isolator: theory and experiment,” Journal of Sound and Vibration, vol. 218, no. 1, pp. 81–101, 1998.

J.-H. Lee and K.-J. Kim, “Modelling of nonlinear complex stiffness of dual-chamber pneumatic spring for precision vibration isolations,” Journal of Sound and Vibration, vol. 301, no. 3–5, pp. 909–926, 2007.

J.-H. Lee and K.-J. Kim, “A method of transmissibility design for dual-chamber pneumatic vibration isolator,” Journal of Sound and Vibration, vol. 323, no. 1-2, pp. 67–92, 2009.

H. Pu, X. Luo, and X. Chen, “Modeling and analysis of dual-chamber pneumatic spring with adjustable damping for precision vibration isolation,” Journal of Sound and Vibration, vol. 330, no. 15, pp. 3578–3590, 2011.

T. Mankovits and T. Szabó, “Finite element analysis of rubber bumper used in air-springs,” Procedia Engineering, vol. 48, pp. 388–395, 2012.

Y. Zhu, Q. Li, D. Xu, C. Hu, and M. Zhang, “Modeling and analysis of a negative stiffness magnetic suspension vibration isolator with experimental investigations,” Review of Scientific Instruments, vol. 83, no. 9, pp. 95–108, 2012.

Q. Li, Research on an innovative vibration isolator with magnetic suspension positive stiffness combined with rubber negative stiffness system [Ph.D. thesis], Tsinghua University, Beijing, China, 2013.

S. S. Rao, Mechanical Vibrations, Tsinghua University Press, Beijing, China, 2009.

L. R. G. Treloar, The Physics of Rubber Elasticity, Oxford University Press, Oxford, UK, 1975.

A. N. Gent, Engineering with Rubber-How to Design Rubber Components, Hanser Publishers, Munich, Germany, 2001.

http://www.mfopt.com/products_detail/productId=2688.html.