Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Feature-Weighted Sampling for Proper Evaluation of Classification Models
Previous Article in Journal
Hole Concealment Algorithm Using Camera Parameters in Stereo 360 Virtual Reality System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Structure of Rubber Ring for Hydraulic Buffer Seal Based on Numerical Simulation

Key Laboratory of Fundamental Science for National Defense-Advanced Design Technology of Flight Vehicle, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(5), 2036; https://doi.org/10.3390/app11052036
Submission received: 19 January 2021 / Revised: 13 February 2021 / Accepted: 14 February 2021 / Published: 25 February 2021

Abstract

:

Featured Application

1. A novel structure of a rubber ring with butterfly geometry is delivered based on numerical simulation, and proved to be better than current used O-ring. 2. The effects of pre-compression rate and hydraulic pressure on the static sealing performance of a rubber ring are obtained. 3. The simulation method of installation boundary for the sealing analysis of a rubber ring is discussed, and a simulation method suitable for rubber ring seal analysis is proposed.

Abstract

Landing gear is a key load-bearing structure of aircraft during ground operation, and the landing capacity of landing gear is determined by the performance of buffer. To solve the problem of buffer failure caused by insufficient static sealing of a rubber ring at groove side, a new structure of a butterfly rubber ring is proposed by analyzing the factors affecting sealing performance of the rubber ring. First, the constitutive equation of rubber material is derived based on the theory of hyper-elasticity, and the material parameters are obtained by fitting the experimental data. Then, by analyzing the simulation method of installation mode and installation stroke, the simulation method suitable for calculating the sealing performance of the rubber ring is established. The linear fitting formulas with Pearson coefficient greater than 0.92 are used to discuss the influence of pre-compression rate and hydraulic pressure on the sealing performance of the rubber ring. Compared with O-ring, the contact pressure of butterfly-ring is increased by 30% in assembly state and 14% in working state. The results show that the butterfly rubber ring has excellent static sealing performance. It is concluded that improving the configuration of the sealing ring can solve the insufficient unilateral sealing of the hydraulic buffer.

1. Introduction

Landing buffer capability of landing gear buffer is the key to evaluate landing gear performance, especially for carrier-based fighter and hypersonic aircraft [1]. After buffering the ground load by the buffer, the damage caused by the ground load on the landing gear can be reduced, and the service life of the landing gear can be improved. The basic requirement of buffer design is good sealing performance [2]. Rubber material has hyperelastic performance, which is widely used in the sealing structure of buffer. As the main seal of hydraulic buffer, the sealing performance of O-ring is directly related to the working performance of buffer [3]. The research on the factors affecting the sealing performance of the rubber ring is helpful to calculate the sealing performance of the rubber ring and improve the accuracy of buffer design.
Due to the nonlinear behavior of material, geometry and contact, the finite element analysis (FEA) is widely used as a research method of rubber ring sealing performance [4]. As a kind of hyperelastic material, the constitutive equations commonly used in numerical calculation are Arruda Boyce, van der Waals, Ogden, Mooney–Rivlin (M-R), etc. Ali reviewed the widely used constitutive models of rubber materials in FEA, and pointed out that M-R model is suitable for small deformation and medium deformation [5]. Huang pointed out that the seal structure is axisymmetric, so the two-dimensional axisymmetric model can be used for the numerical calculation of the seal ring. Also, the two-dimensional axisymmetric model of the seal structure was established by using ANSYS software, and the validity of the two-dimensional axisymmetric numerical model is verified by the FEA of the sealing surface and the back support structure of the roller bit bearing seal [6,7]. Zhou established five contact pairs and analyzed the state of contact surface with the form of pseudo element. By comparing the numerical results with the experimental data, Zhou pointed out that it is feasible to simulate the contact pressure distribution of bearing seal structure through contact pairs [8]. The FEA method is an effective method to study the sealing performance of the rubber ring [9].
The most important factor in the design of hydraulic buffer seal structure is the configuration of the rubber ring. The configuration of the rubber ring determines the contact peak pressure and contact length, both of which are used as the evaluation index of seal performance. Angus Jean analyzed the influence of geometric shape on the deformation of rubber structure, and pointed out that different structural shapes have a great influence on the sealing performance of the rubber ring [10]. A lot of work has been done in the structural design of rubber sealing ring, and some new rubber sealing rings have been used in buffer seal, such as D-ring, T-ring, U-ring, X-ring, Y-ring, etc. The geometry of D-ring, which is designed to increase the contact length between groove and D-ring, could be seen as a half of O-ring. Zhou compared D-ring with O-ring in the sealing performance, and concluded the sealing performance of D-ring is better than O-ring from the perspective of contact peak stress. Also pointed the disadvantage of extrusion through buffer gap for D-ring [11]. Mose analyzed the extrusion of D-ring with hydraulic pressure 5.89 MPa, and concluded the extrusion is likely to cause sealing failure. The sealing disadvantage of D-ring is laid on stress concentration at corner under static sealing and extrusion through gap under dynamic sealing [12]. The geometric shape of X-ring is benefit for dynamic condition by decreasing friction. Shin showed the higher contact stress of X-ring than O-ring under certain research condition, but the ultimate pressure of extrusion is only 3.92 MPa for X-ring [13]. The disadvantage of X-ring is the contact length under low pressure, which is harmful to the sealing of the rubber ring [14]. Furthermore, the X-ring is not suitable for landing gear buffer, due to the limitation of hydraulic pressure. The sealing characteristic of Y-ring relied on the surface of lips to coupling, as a lip-type rubber ring. Cui analyzed the effect of hydraulic pressure on Y-ring sealing performance, and found out the location of the maximum stress and largest deformation area according to seal failure criterion [15]. The disadvantage of thresh tendency under high pressure is a weakness for Y-ring. The type of O-ring is widely used in hydraulic buffer of landing gear, due to the stability of contact length and pressure [16]. According to test data on 21 joint seals, Jahangir concluded that the properties of joint seals are differed by geometry shape and manufacture method [17]. As the seal of hydraulic buffer for landing gear, the disadvantage of lower contact pressure at groove side than piston rod side still needs to be modified for O-ring during assemble process [18].
To improve the sealing pressure at groove side of articulated landing gear buffer, the numerical model of rubber ring sealing analysis is discussed, the sealing characteristics of the rubber ring are analyzed, and the sealing performance of proposed butterfly-ring and original O-ring is compared. First, the stress–strain formula of the rubber ring is derived based on the theory of hyper-elasticity, and the material parameters of M-R are obtained by using compression test data from reference. Secondly, using nonlinear numerical technique, three typical simulation methods of installation mode and installation stroke are discussed. On this basis, the numerical model of rubber sealing ring is established, and 16 groups of 154 simulations are carried out. Based on the numerical test of O-ring, the relationship between the sealing characteristics and load of O-ring in installation mode and working mode is discussed, including pre-compression rate and hydraulic pressure. And then, by comparing with O-ring, the sealing performance of butterfly-ring is discussed.

2. Materials and Structure

2.1. Material

The material used in this buffer seal is Nitrile rubber(NBR), which is widely used in landing gear buffer, considered the pressure resistance of NBR. Based on the assumption of hyper-elasticity theory, the deformation and stress of rubber-like material which is employed by the strain–energy density function, is written as a function of the strain invariant [19].
W = f I 1 , I 2 , I 3
where W is the strain–energy function, λ represents for the invariant of the Cauchy–Green strain tensor [20]. Furthermore, the invariants of the Cauchy–Green strain tensor could be expressed by the principal stretches ratios λ 1 , λ 2 , λ 3 .
I 1 = λ 1 2 + λ 2 2 + λ 3 2 I 2 = λ 1 λ 2 2 + λ 2 λ 3 2 + λ 3 λ 1 2 I 3 = λ 1 λ 2 λ 3 2
In addition, the formula of the Cauchy stress tensor σ i j is expressed in terms of the principal stretches,
σ 1 = 2 λ 1 W I 1 + λ 2 2 + λ 3 2 W I 2 + λ 2 2 λ 3 2 W I 3 σ 2 = 2 λ 2 W I 1 + λ 1 2 + λ 3 2 W I 2 + λ 1 2 λ 3 2 W I 3 σ 3 = 2 λ 3 W I 1 + λ 1 2 + λ 2 2 W I 2 + λ 1 2 λ 2 2 W I 3
The two parameter M-R model for approximately incompressible material is adopted,
W = C 10 I 1 3 + C 01 I 2 3 + C 10 2 + C 01 1 I 3 2 1 + C 10 ( 5 ν 2 ) + C 01 ( 11 ν 5 ) 2 ( 1 2 ν ) I 3 1 2
where C 10 and C 01 are the material parameters. As the rubber material is assumed to be approximately incompressible [21], i.e., the Poisson ratio ν = 0.499 0.5 , I 3 = 1 . Therefore, the two parameter M-R model is written as,
W = C 10 I 1 3 + C 01 I 2 3
In addition, the uniaxial compression tests,
λ 1 = 1 λ 2 2 = 1 λ 3 2 = λ σ 1 = 2 λ 1 W I 1 + λ 2 2 + λ 3 2 W I 2 σ 2 = σ 3 = 0
where λ is the strain value. By Equations (5) and (6), the Cauchy stress is developed as:
σ = 2 λ 4 λ 3 1 C 10 + 2 λ 3 λ 3 1 C 01
According to the data of rubber compression experiment quoted from Reference [8], the material parameters of M-R are fitted as C 10 = 1.84 , C 01 = 0.47 with Equation (7).

2.2. Structure

The hydraulic buffer is the key device of articulated landing gear which is used to bear and transmit the vertical load from ground to fuselage. As presented in Figure 1, the rubber ring is located at the bottom of the hydraulic buffer, which plays a sealing role to prevent the leakage of the hydraulic system. The structure of O-ring is adopted as the type of inner seal and outer seal for hydraulic buffer.

3. Modeling

First, the nonlinear finite element model of seal structure is established. Then, three typical installation modes are simulated to establish the suitable installation mode, and hydraulic pressure is applied to the deformed rubber ring. On this basis, the sealing performance of the original structure and the proposed structure are analyzed and compared.

3.1. Finite Element Model

According to the rubber ring of hydraulic buffer used in articulated landing gear, considered the axisymmetric characteristic of geometry and load condition [22,23], the axisymmetric model of rubber O-ring is established with section diameter 7 mm as shown in Figure 2a. Specific geometry size of groove and piston rod are set according to the part 1 of ISO 3601 standard.
HyperMesh v.14.0 program is used to mesh the structure, because the mesh size has a direct impact on the numerical results. Therefore, by reducing the mesh size from 0.4 to 0.05, the mesh convergence is studied in 0.05 steps. When the mesh size is less than 0.2, the numerical results have no obvious change. Therefore, the mesh size is 0.2, the number of meshes for O-ring sealing structures are 4985, and the number of meshes for Butterfly-ring sealing structures are 7529. The APDL statement of ANSYS v.16.2 program is used to complete the parameter definition and solution. In the 2D axisymmetric analysis, plane 182 element is used. The material model of seal ring is the Mooney–Rivlin Hyperelastic with two parameters, and the material parameters are 1.84 and 0.47, which have been discussed in Section 2.1. The material of piston rod and groove is STS 430 steel (Elastic modulus 210 GPa, Poisson ratio 0.3). Loads applied and boundary conditions are the variable in this study, in which the pre-compression ratio is in the range of 8% to 16%, and the hydraulic pressure is in the range of 0 to 20 MPa. The contact pressure of groove side is the major cause of hydraulic leakage during assemble process, in order to enhance the contact peak pressure, the butterfly-ring is designed to elevate the contact pressure of groove side as Figure 2b shown. The maximum section diameter of butterfly-ring is the same as that of O-ring, which is 7 mm. Specific parameters are illustrated in the Appendix A, Appendix B and Appendix C.
The contact pair, which is adopted to simulate the contact nonlinear behavior [24,25,26], is defined at the squeeze area during the installation mode and working mode, as shown in Figure 3. And the element type of TARGE169 and CONTA171 are adopted to simulate the contact behavior. The rubber ring is set as contact surface [27], due to its stiffness being softer than the material of piston rod and groove (Elastic modulus 210 GPa). Augmented Lagrange method, which could lead to better conditioning than the pure penalty method, is adopted to solve the contact problem of stiffness hardening for large displacement and deformation [28,29]. The material of the groove and piston rod is STS 430 steel, the surface roughness of the piston rod is 0.05 to 0.1 μ m, and the surface roughness of the groove is 0.2 to 0.3 μ m. Considering that the surface roughness of piston rod is less than that of groove surface, the friction coefficient between piston rod and sealing ring is set as 0.1, and the friction coefficient between groove and sealing ring is set as 0.2.

3.2. Installation Mode

The boundary condition of FEA for the rubber ring is classified into two stage by installation mode and working mode. In order to study the deformation-stress of the rubber ring caused by boundary condition of assemble status, three kind of installation modes are assigned according to the Reference [30]. As schematic Figure 4 shown, piston rod axial displacement (RA), piston rod lateral displacement (RL), and groove lateral displacement (GL) are taken into account [15]. The conventional installation modes are presented in Figure 4b,c, e.g., piston rod is set with a rigid displacement in lateral direction and groove is set with a fixed constraint considered installation mode of RL. The RL and GL are considered nearly the same based on the theory of relative displacement. The detailed simulation parameters of installation mode are shown in Table A1 of the Appendix A. Additionally, piston rod is set with a rigid displacement in axial direction and groove is set with a fixed constraint considered installation mode of RA. The detailed simulation parameters are shown in Table A2 and Table A3 of the Appendix A and Appendix B.

3.3. Working Mode

The boundary condition of working mode is set as static state with the constraint of piston rod and groove, hydraulic pressure on the rubber ring, shown in Figure 5.
It should be noted that the applied area of hydraulic pressure depends on the deformation of the rubber ring after the installation mode. Therefore, after the completion of installation simulation, the hydraulic pressure should be applied according to the hydraulic pressure action region in Figure 5, and then the simulation of working mode should be executed [31,32,33]. Additionally, the hydraulic pressure applied on the rubber ring is differed in different experiment within the range of 0 to 20 MPa, as illustrated in Table A4 of the Appendix C.

4. Results

The influence of installation mode, installation stroke, pre-compression rate and hydraulic pressure on the sealing performance of the rubber ring are studied, 16 projects totally 154 experiments are conducted with numerical method [34]. Specific parameters are illustrated in Appendix. The value and corresponding code of variables are given in Table 1.

4.1. Influence of Installation Mode

As Figure 6 shown, the contact length under RL mode is compatible to that under GL mode at piston rod side. When the pre-compression rate is 8%, 12% and 16%, the contact length under RA mode is 12%, 16% and 16% longer than that under RL mode, respectively. And the regular of contact length at groove side is compatible with the regular of contact length curve at piston rod side. Furthermore, the difference of contact length at groove side is smaller than that at piston rod side between RA and RL. Additionally, the difference of contact length between RA mode and RL mode is smaller at the groove side than at the piston rod side.
As Figure 7 shown, the regular of contact peak pressure at groove side is much compatible with that at piston rod side. Also, the regular of contact peak pressure at piston rod side and groove side is just like that of contact length. The contact peak pressure under RA mode is higher than that under RL mode and GL mode at piston rod side in Figure 7a. The smallest increase rate is 20% compared with RL under pre-compression rate 8%.
Only a half of the rubber ring is compressed with installation stroke 3 mm, thus the experiment 5EM3 is disqualified. As shown in Figure 8a, with the increase of installation stroke, the increase of contact length at piston rod side, groove side and bottom sides become slower. When the installation stroke is increased to 9 mm, the limit value is reached. As shown in Figure 8b, with the increase of installation stroke, the contact peak pressure decreases to the limit value. When the installation stroke is increased to 9 mm, the limit value is reached. Additionally, the limit values of the contact peak pressure and the contact length are both obtained with the installation stroke of 9 mm.

4.2. Influence of Pre-Compression Rate

As shown in Figure 9a, the contact peak pressure at piston rod side increases linearly with the pre-compression rate, which can be fitted by linear fitting,
y = 0.45 + 0.31 x
where y is contact peak pressure, x is the corresponding pre-compression rate. The linear formula is fitted with Pearson’s 0.92. And the variation range of contact peak pressure is reflected by the bandwidth, which is in the range of 0.65 MPa to 1.69 MPa.
As shown in Figure 9b, the contact peak pressure at groove side is concentrated in a band for experiment 1M to 5M, also fitted by a linear formula,
y = 0.59 + 0.28 x
The Pearson’s is 0.97 for the fitted linear formula of contact peak pressure at groove side. The band is narrowed within the range of 0.29 MPa to 0.67 MPa.

4.3. Influence of Hydraulic Pressure

According to the criterion that contact peak pressure should be higher than the hydraulic pressure, the failure zone of seal is defined. As shown in Figure 10a, the contact peak pressures of experiment from 5A to 5E are in the narrow band range of 2.10 MPa to 3.46 MPa. And in the range of hydraulic pressure from 0 to 20 MPa, the contact peak pressures are always distributed outside the failure zone. The relationship between contact peak pressure and hydraulic pressure is fitted as follows:
y = 4.71 + 0.99 x
where y is the contact peak pressure, x is the hydraulic pressure. Also, the similar relation is fitted for grove side,
y = 3.90 + 1.20 x
For the piston rod side and groove side, the relationship between the contact peak pressure and the hydraulic pressure are both conform to Pearson’s 0.99.

4.4. Sealing Performance of Butterfly-Ring

To compare the sealing performance of butterfly-ring and O-ring, the sealing simulations are carried out under assembly mode (hydraulic pressure 0 MPa) and working mode (hydraulic pressure 5 MPa). The maximum section diameter of the two types of seals is 7 mm, and the RA installation mode is adopted with the installation stroke of 9 mm. Considering that the pre-compression rate is the main factor affecting the sealing performance of the rubber ring, the pre-compression rates from 8% to 16%. Therefore, the experiment of butterfly-ring performance is composed of 18 simulation experiments of 2 groups.
The rubber rings at the piston rod side and at groove side are required to have good sealing performance, especially at the groove side which is used to be found leakage in hydraulic buffer of articulated landing gear. The sealing performance is compared under installation mode and working mode as Figure 11 shown.
It can be found that the contact peak pressure of butterfly-ring is greater than that of O-ring at the maximum compression rate. When the pre-compression rate is 16%, the contact peak pressure at the groove side of butterfly-ring increases by 14%. And this law holds for both piston rod side and groove side. Moreover, it is worth noting that the proposed butterfly-ring improves the contact peak pressure at groove side and keeps the contact pressure at piston rod side close to that of the O-ring. Specifically, when the pre-compression rate is 8%, the contact peak pressure of butterfly-ring is 0.15 MPa lower than that of O-ring. When the pre-compression rate is 16%, the contact peak pressure of butterfly-ring is 0.45 MPa higher than that of O-ring. When the pre-compression rate is 14%, the contact peak pressure of butterfly-ring and O-ring is basically the same. As a whole, the contact peak pressure of butterfly-ring is close to that of O-ring at piston rod side, as shown in Figure 11a.
At the groove side, the contact peak pressure of butterfly-ring is always higher than that of O-ring, as shown in Figure 11b. Compared with O-ring, the contact peak pressure of butterfly-ring at groove side is increased about 30% under working status.

5. Discussion

5.1. Evaluation of Installation Mode

It can be drawn that the contact length and contact peak pressure caused by installation mode of RA are higher than that caused by installation mode of RL and GL. Considered the function between piston rod and rubber ring [35], the different interaction force for the RA, RL and GL installation mode are given in Figure 12.
The deformation of the rubber ring under RA mode is caused by shear force and squeeze force together. The deformation of the rubber ring under RL and GL is nearly the same, due to the deformation are both caused by squeeze force as illustrated in Figure 12b [36]. The deformation of the rubber ring is affected by the stiffness of rubber. For rubber material, the stiffness caused by squeeze force are different with that caused by shear force. P. B. Lindley’s research shows that the compression stiffness is more sensitive to displacement than the shear stiffness, and the relationship between the compression stiffness and displacement is of the fifth power, while the relationship between the shear stiffness and displacement is of the first power.
Thus, the pressure caused by compression stiffness would be much higher than the pressure caused by shear stiffness. It is consistent with the previous results that the contact peak pressure under RA is slightly higher than the contact peak pressure under RL and GL. Due to the contact peak pressure and contact length caused by RA, RL and GL are different, considered the three-dimensional hydraulic buffer as shown in Figure 1, the RA installation mode is approximately to the actual assemble [37]. Therefore, RA installation mode is adopted as the installation mode for FEA of rubber ring.
In addition, the installation stroke should be considered in the simulation of the rubber ring. As shown in Figure 8, the contact peak pressure arrives at the limit value with installation stroke 9 mm. Additionally, the limit value of contact length for groove side, piston rod side, and bottom side are also got stable at the installation stroke 9 mm. From the perspective of load transmission, the deformation and stress of the rubber ring is mainly triggered by installation stroke via the squeeze state. The squeeze state is simulated with contact pair defined among piston rod, rubber ring, groove as Figure 3 presented. Then, the squeeze state is calculated via the contact pair with Augmented Lagrange. In this method, the stiffness of the contact pair is rebalanced after each iteration, so the rubber ring needs enough installation stroke to reach a stable state. Therefore, the installation stroke should be 9mm in this study [38]. For different structure of the rubber rings, the installation stroke should be obtained through trial calculation.

5.2. Evaluation of Pre-Compression Rate and Hydraulic Pressure

The results of pre-compression rate experiment show that simply increasing the section diameter of the rubber ring cannot effectively increase the contact peak pressure, and properly increasing the pre-compression rate is effective for increasing the contact peak pressure. P. B. Lindley’s research shows that the relationship between compression stiffness and pre-compression rate is approximately a fifth power, so the compression stiffness of different section diameters are approximately equal under the same pre-compression rate. As shown in Figure 9, the contact peak pressure of simulation from 1M to 5M is concentrated in a narrow band. And the first-degree polynomial is used to fit the function of pre-compression rate to contact peak pressure, and the Pearson’s value is above 0.92. Obviously, the contact pressure increased by the pre-compression rate is higher than that by the section diameter.
The influence of hydraulic pressure on the contact peak pressure of the rubber ring is caused by the influence of compression stiffness. In the working mode, the compression stiffness of the rubber ring is affected by the pre-compression rate and the hydraulic pressure, so the compression stiffness of the rubber ring under the hydraulic pressure is larger than that under the pre-compression rate. Accordingly, the contact peak pressure increases with the increase of hydraulic pressure. This point is also reflected in the fitting formula at groove side. The increase rate of contact peak pressure by hydraulic pressure (1.2) is greater than that by pre-compression rate (0.28). Correspondingly, 0.99 to 0.31 at piston rod side.
Although increasing the hydraulic pressure can improve the sealing pressure, but in engineering applications, the amplitude of hydraulic pressure is determined by the working load, and it cannot always be maintained at high value. Therefore, increasing the pre-compression rate is still an effective way to improve the sealing pressure.

5.3. Evaluation of Butterfly-Ring on Sealing Performance

Compared with O-ring, the proposed butterfly-ring improves the contact peak pressure at groove side. Because the stiffness of the rubber ring varies with its geometry, and the pre-deformed butterfly-ring improves the compression stiffness and is prone to lateral deformation, so the contact peak pressure of the butterfly-ring is higher than that of the O-ring at groove side. This is also confirmed by the calculation results, as shown in Figure 11. Under installation mode, the contact peak pressure of the butterfly-ring is increased by 30% at groove side, and the contact peak pressure increases by 14–16% under working mode.
In addition, the increase of contact peak pressure of butterfly seal to load is lower than that of O-ring. It shows that the sensitivity of butterfly seal is less than that of O-ring.

6. Conclusions

To solve the insufficient sealing at the groove side of articulated landing gear buffer, the sealing performance of the proposed butterfly-ring is discussed. It is found that the pre-deformation structure of butterfly-ring improves the contact peak pressure at groove side, which is beneficial to improve the sealing performance of buffer. Also, it is pointed out that the installation mode of RA is suitable as the boundary condition for the numerical analysis of rubber ring sealing, and the installation stroke should be obtained from the extreme value of contact length and seal pressure curve. Furthermore, the linear law of contact peak pressure with pre-compression rate, and contact peak pressure with hydraulic pressure are fitted, which can be used to calculate the pre-compression rate under required conditions.

Author Contributions

L.Z. is responsible for the theoretical modeling, numerical analysis, and academic writing. X.W. provides the academic and experimental guidance. ‘A novel structure of the rubber ring for hydraulic buffer seal based on numerical simulation’. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 51805249).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
FEAFinite element analysis
M-RMooney–Rivlin
NBRNitrile rubber
RAPiston rod axial displacement
RLPiston rod lateral displacement
GLGroove lateral displacement
WStrain–energy function
IInvariant of the Cauchy–Green strain tensor
I 1 , I 2 , I 3 Three invariants of the Green deformation tensor
λ 1 , λ 2 , λ 3 Three principal stretches ratios
σ i j Cauchy stress tensor
σ 1 , σ 2 , σ 3 Three principal Cauchy stresses
ν Poisson ratio
c 01 First material parameter of M-R model
c 10 Second material parameter of M-R model

Appendix A. Parameter and Design of Installation Mode Simulation

Table A1. Parameters of installation mode experiment.
Table A1. Parameters of installation mode experiment.
Installation ModeDirectionPre-Compression Rate (%)
RANegative of Y
RLNegative of X8910111213141516
GLPositive of X
The difference of the three installation modes are the different objects to which the displacement boundary are applied, and the influence of different installation modes on hyperelastic material need to be further discussed. In order to set a precise boundary condition for finite element analysis, three typical installation modes, including RA, RL and GL, are compared and analyzed. In the simulation of installation mode, three kinds of pre-compression rates of the rubber ring with section diameter of 7 mm are selected according to the ISO 3601 standard. The installation simulation is carried out under the parameters shown in Table A1.
Table A2. Parameters and labels of installation stroke experiment.
Table A2. Parameters and labels of installation stroke experiment.
ProjectLabelBasic ParameterInstallation Stroke (mm)
5EM5EM3 Section diameter: 7 mm3
5EM55
5EM66
5EM77
5EM8Pre-compression: 16%8
5EM99
5EM1010
5EM1111
5EM1212
Since the master displacement of RA mode is the axial direction displacement of piston rod, the installation stroke should be long enough to ensure that the rubber ring is in a stable compression state. In order to determine the appropriate value of installation stroke, the section diameter and pre-compression ratio are designed according to the maximum allowable value of the buffer, and the parameters and labels are shown in Table A2.

Appendix B. Parameter and Design of Pre-Compression Rate Simulation

To analyze the influence of pre-compression rate on the sealing performance of the rubber ring, a series simulation are designed according to ISO 3601 standard. Considering that the influence of pre-compression rate on sealing performance of the rubber ring varies with the change of section diameter, and the range of section diameter is set from 1.8 mm to 7 mm. The maximum pre-compression rate is set as 16%, which is the maximum pre-compression rate of the rubber ring with section diameter from 1.8 mm to 7 mm. A total of 45 simulations have been carried out. As the regular of parameters and labels design are the same, only the parameters and labels of 1M simulation are given, as shown in Table A3.
Table A3. Parameters and labels of pre-compression rate experiment.
Table A3. Parameters and labels of pre-compression rate experiment.
ProjectLabelSection DiameterPre-Compression Rate (%)
1M1AM1.8 mm8
1BM9
1CM10
1DM11
1EM12
1FM13
1GM14
1HM15
1KM16

Appendix C. Parameter and Design of Hydraulic Pressure Simulation

Table A4. Parameters and Labels of Hydraulic Pressure Experiment.
Table A4. Parameters and Labels of Hydraulic Pressure Experiment.
ProjectLabelPre-Compression RateHydraulic Pressure (MPa)
5A5AM8%0
5AN2
5AO4
5AP6
5AQ8
5AR10
5AS12
5AT14
5AU16
5AV18
5AW20
To analyze the influence of hydraulic pressure on static sealing performance of the rubber ring, 55 simulations of 5 groups are carried out under working mode. According to the static sealing requirement of radial seal in ISO 3601, the limit hydraulic pressure is set to 20 MPa. Since the setting of parameters and labels is regularly, only the parameters and labels of 5A simulation are given in Table A4.

References

  1. Nikas, G.K. Eighty years of research on hydraulic reciprocating seals: Review of tribological studies and related topics since the 1930s. Proc. Inst. Mech. Eng. Part J J. Eng. Tribol. 2010, 224, 1–23. [Google Scholar] [CrossRef]
  2. Zhou, C.; Chen, G.; Xiao, S.; Hua, Z.; Gu, C. Study on fretting behavior of rubber O-ring seal in high-pressure gaseous hydrogen. Int. J. Hydrogen Energy 2019, 44, 22569–22575. [Google Scholar] [CrossRef]
  3. Zeng, D.; He, Q.; Li, T.; Hu, J.; Shi, T.; Zhang, Z.; Yu, Z.; Liu, R. Corrosion mechanism of hydrogenated nitrile butadiene rubber O-ring under simulated wellbore conditions. Corros. Sci. 2016, 107, 145–154. [Google Scholar] [CrossRef]
  4. Kömmling, A.; Jaunich, M.; Pourmand, P.; Wolff, D.; Hedenqvist, M. Analysis of O-Ring Seal Failure under Static Conditions and Determination of End-of-Lifetime Criterion. Polymers 2019, 11, 1251. [Google Scholar] [CrossRef] [Green Version]
  5. Ali, A.; Hosseini, M.; Sahari, B.B. A Review of Constitutive Models for Rubber-Like Materials. Am. J. Eng. Appl. Sci. 2010, 3, 232–239. [Google Scholar] [CrossRef] [Green Version]
  6. Huang, Z.; Li, G. Optimization of cone bit bearing seal based on failure analysis. Adv. Mech. Eng. 2018, 10, 168781401876748. [Google Scholar] [CrossRef]
  7. Hao, S.; Yang, J. A Novel Quality Requirement Design Method for the Quality Characteristic of Rubber Products Based on the Reliability Constraint. IEEE Access 2018, 6, 17887–17895. [Google Scholar] [CrossRef]
  8. Zhou, Y.; Huang, Z.; Tan, L.; Ma, Y.; Qiu, C.; Zhang, F.; Yuan, Y.; Sun, C.; Guo, L. Cone bit bearing seal failure analysis based on the finite element analysis. Eng. Fail. Anal. 2014, 45, 292–299. [Google Scholar] [CrossRef]
  9. Yildiz, Y.; Duzgun, M. Stress analysis of ventilated brake discs using the finite element method. Int. J. Automot. Technol. 2010, 11, 133–138. [Google Scholar] [CrossRef]
  10. Jeang, A.; Chung, C.P.; Ko, C.P. Parameters decision for an assembly under mechanical structure deformation via probabilistic optimization. Adv. Mech. Eng. 2019, 11, 168781401987924. [Google Scholar] [CrossRef] [Green Version]
  11. Zhou, C.; Chen, G.; Liu, P. Finite Element Analysis of Sealing Performance of Rubber D-Ring Seal in High-Pressure Hydrogen Storage Vessel. J. Fail. Anal. Prev. 2018, 18, 846–855. [Google Scholar] [CrossRef]
  12. Mose, B.R.; Hawong, J.S.; Nam, J.H. Stress analysis of a stepped rounded D-ring with a ratio of H1/H2 = 3.0 under uniform squeeze rate and internal pressure by photoelastic experimental hybrid method. Mater. Und Werkst. 2013, 44, 861–878. [Google Scholar] [CrossRef]
  13. Shin, D.C.; Hawong, J.S.; Lee, S.W.; Bernard, A.O.; Lim, H.S. Contact behavior analysis of X-ring under internal pressure and uniform squeeze rate using photoelastic experimental hybrid method. J. Mech. Sci. Technol. 2014, 28, 4063–4073. [Google Scholar] [CrossRef]
  14. Bernard, A.O.; Hawong, J.S.; Shin, D.C.; Dong, B. Contact behavior analysis of elastomeric x-ring under uniform squeeze rate and internal pressure before and after forcing-out using the photoelastic experimental hybrid method. J. Mech. Sci. Technol. 2015, 29, 2157–2168. [Google Scholar] [CrossRef]
  15. Cui, K.B.; Qin, J.Q.; Di, C.C.; Yang, Y.F. Finite Element Analysis and Simulation of the Sealing Performance of Y-Ring Rubber Seal. Appl. Mech. Mater. 2013, 444–445, 1379–1383. [Google Scholar] [CrossRef]
  16. Zhang, S.X.; Wang, J.J.; Gao, G.; Parke, K.; Ma, J.X. Pigment epithelium-derived factor downregulates vascular endothelial growth factor (VEGF) expression and inhibits VEGF-VEGF receptor 2 binding in diabetic retinopathy. J. Mol. Endocrinol. 2006, 37, 1–12. [Google Scholar] [CrossRef] [Green Version]
  17. Mirza, J. Joint seals for hydraulic structure in sever climates. J. Civ. Eng. Manag. 2014, 20, 38–46. [Google Scholar] [CrossRef]
  18. Yamabe, J.; Nishimura, S. Tensile Properties and Swelling Behavior of Sealing Rubber Materials Exposed to High-Pressure Hydrogen Gas. J. Solid Mech. Mater. Eng. 2012, 6, 466–477. [Google Scholar] [CrossRef] [Green Version]
  19. Orzechowski, G.; Frączek, J. Nearly incompressible nonlinear material models in the large deformation analysis of beams using ANCF. Nonlinear Dyn. 2015, 82, 451–464. [Google Scholar] [CrossRef] [Green Version]
  20. Xu, H.; Sin, F.; Zhu, Y.; Barbič, J. Nonlinear material design using principal stretches. ACM Trans. Graph. 2015, 34, 1–11. [Google Scholar] [CrossRef]
  21. Yang, H.; Yao, X.F.; Ke, Y.C.; Ma, Y.J.; Liu, Y.H. Constitutive behaviors and mechanical characterizations of fabric reinforced rubber composites. Compos. Struct. 2016, 152, 117–123. [Google Scholar] [CrossRef]
  22. Maureira, N.; de La Llera, J.; Oyarzo, C.; Miranda, S. A nonlinear model for multilayered rubber isolators based on a co-rotational formulation. Eng. Struct. 2017, 131, 1–13. [Google Scholar] [CrossRef]
  23. Yaya, K.; Bechir, H. A new hyper-elastic model for predicting multi-axial behaviour of rubber-like materials: Formulation and computational aspects. Mech. Time-Depend. Mater. 2018, 22, 167–186. [Google Scholar] [CrossRef]
  24. Tiwari, A.; Dorogin, L.; Tahir, M.; Stöckelhuber, K.W.; Heinrich, G.; Espallargas, N.; Persson, B.N.J. Rubber contact mechanics: Adhesion, friction and leakage of seals. Soft Matter 2017, 13, 9103–9121. [Google Scholar] [CrossRef] [PubMed]
  25. Sanjaya, D.P.; Fidkowski, K.J. Improving High-Order Finite Element Approximation Through Geometrical Warping. AIAA J. 2016, 54, 3994–4010. [Google Scholar] [CrossRef] [Green Version]
  26. Lee, C.R.; Jeong, H.Y. Development of Headform Impactor Finite Element Model Considering the Hyperelastic and Viscoelastic Responses of Rubber. Int. J. Automot. Technol. 2018, 19, 523–534. [Google Scholar] [CrossRef]
  27. Adams, G.G.; Nosonovsky, M. Contact modeling—Forces. Tribol. Int. 2000, 33, 431–442. [Google Scholar] [CrossRef]
  28. Belhocine, A.; Duzgun, M.; Ghazaly, N.M.; Abdullah, O.I. RETRACTED ARTICLE: Structural and contact analysis of a three-dimensional disc-pad model with and without thermal effects. Friction 2015. [Google Scholar] [CrossRef] [Green Version]
  29. Kawasaki, S.; Tada, T.; Persson, B.N.J. Adhesion and friction between glass and rubber in the dry state and in water: Role of contact hydrophobicity. Soft Matter 2018, 14, 5428–5441. [Google Scholar] [CrossRef]
  30. Shimanovsky, A.O.; Abdulkader, M.H.; Kuzniatsova, M.G. Finite Element Modelling of Contact Interaction between Spherical Indenter and Elastic-Plastic Body. Appl. Mech. Mater. 2015, 797, 307–313. [Google Scholar] [CrossRef]
  31. Zhang, H.; Zhang, J. Static and Dynamic Sealing Performance Analysis of Rubber D-Ring Based on FEM. J. Fail. Anal. Prev. 2016, 16, 165–172. [Google Scholar] [CrossRef]
  32. Yamabe, J.; Nishimura, S. Influence of fillers on hydrogen penetration properties and blister fracture of rubber composites for O-ring exposed to high-pressure hydrogen gas. Int. J. Hydrogen Energy 2009, 34, 1977–1989. [Google Scholar] [CrossRef]
  33. Yamabe, J.; Koga, A.; Nishimura, S. Failure behavior of rubber O-ring under cyclic exposure to high-pressure hydrogen gas. Eng. Fail. Anal. 2013, 35, 193–205. [Google Scholar] [CrossRef]
  34. Jana, T.; Mitra, A.; Sahoo, P. Dynamic analysis of elastically and plastically graded spherical and cylindrical contact. Proc. Inst. Mech. Eng. Part J J. Eng. Tribol. 2019, 233, 1712–1728. [Google Scholar] [CrossRef]
  35. Yeoh, O.H.; Pinter, G.A.; Banks, H.T. Compression of Bonded Rubber Blocks. Rubber Chem. Technol. 2002, 75, 549–562. [Google Scholar] [CrossRef]
  36. Brizmer, V.; Kligerman, Y.; Etsion, I. The effect of contact conditions and material properties on the elasticity terminus of a spherical contact. Int. J. Solids Struct. 2006, 43, 5736–5749. [Google Scholar] [CrossRef] [Green Version]
  37. Belhocine, A.; Omar, W.Z.W. Three-dimensional finite element modeling and analysis of the mechanical behavior of dry contact slipping between the disc and the brake pads. Int. J. Adv. Manuf. Technol. 2017, 88, 1035–1051. [Google Scholar] [CrossRef]
  38. Kožar, I.; Rukavina, T. The effect of material density on load rate sensitivity in nonlinear viscoelastic material models. Arch. Appl. Mech. 2019, 89, 873–883. [Google Scholar] [CrossRef]
Figure 1. Schematic of seals for articulated landing gear. The hydraulic buffer is a key structure of load-bearing in the articulated landing gear. And the hydraulic buffer of articulated landing gear is composed of buffer tube, groove, seals, and piston rod. To prevent the leakage of hydraulic fluid, the rubber O-rings are used as the inner seal and outer seal.
Figure 1. Schematic of seals for articulated landing gear. The hydraulic buffer is a key structure of load-bearing in the articulated landing gear. And the hydraulic buffer of articulated landing gear is composed of buffer tube, groove, seals, and piston rod. To prevent the leakage of hydraulic fluid, the rubber O-rings are used as the inner seal and outer seal.
Applsci 11 02036 g001
Figure 2. Finite element model. The original structure of seal is O-ring in (a), and the proposed structure of seal is butterfly-ring in (b), with the same maximum diameter as the O-ring. Both are assembled with groove in green mesh and piston rod in yellow mesh.
Figure 2. Finite element model. The original structure of seal is O-ring in (a), and the proposed structure of seal is butterfly-ring in (b), with the same maximum diameter as the O-ring. Both are assembled with groove in green mesh and piston rod in yellow mesh.
Applsci 11 02036 g002
Figure 3. Modeling of contact pair. Three contact pairs are modeled, contact pair 1 is between piston rod and rubber ring, contact pair 2 is between groove side and rubber ring, contact pair 3 is between groove bottom and rubber ring.
Figure 3. Modeling of contact pair. Three contact pairs are modeled, contact pair 1 is between piston rod and rubber ring, contact pair 2 is between groove side and rubber ring, contact pair 3 is between groove bottom and rubber ring.
Applsci 11 02036 g003
Figure 4. Modeling of installation mode. Three installation modes are used to simulate the process of pre-compression in public literature. RA mode means that the piston rod is applied with axial displacement and the groove is constraint with fixed displacement in (a). And RL mode means that the piston rod is applied with lateral displacement and the groove is fixed still as (b). The mode of GL means that the piston rod is fixed, and the groove is applied with lateral displacement in (c).
Figure 4. Modeling of installation mode. Three installation modes are used to simulate the process of pre-compression in public literature. RA mode means that the piston rod is applied with axial displacement and the groove is constraint with fixed displacement in (a). And RL mode means that the piston rod is applied with lateral displacement and the groove is fixed still as (b). The mode of GL means that the piston rod is fixed, and the groove is applied with lateral displacement in (c).
Applsci 11 02036 g004
Figure 5. Modeling of working mode. The pressure of hydraulic fluid, which is represented with number 2, is applied to the deformed surface of the rubber ring by pre-compression function. And the piston rod numbered 1 and the groove numbered 3 are applied with fixed constraint.
Figure 5. Modeling of working mode. The pressure of hydraulic fluid, which is represented with number 2, is applied to the deformed surface of the rubber ring by pre-compression function. And the piston rod numbered 1 and the groove numbered 3 are applied with fixed constraint.
Applsci 11 02036 g005
Figure 6. Comparison of installation mode on contact length. The contact lengths of the three installation modes are shown in (a) for piston rod side and in (b) for groove side, respectively. The three installation modes, including piston rod with axial displacement (RA mode), piston rod with lateral displacement (RL mode), and groove with lateral displacement (GL mode). The range of pre-compression rate is set between the maximum and minimum allowable value.
Figure 6. Comparison of installation mode on contact length. The contact lengths of the three installation modes are shown in (a) for piston rod side and in (b) for groove side, respectively. The three installation modes, including piston rod with axial displacement (RA mode), piston rod with lateral displacement (RL mode), and groove with lateral displacement (GL mode). The range of pre-compression rate is set between the maximum and minimum allowable value.
Applsci 11 02036 g006
Figure 7. Comparison of installation mode on contact peak pressure. The contact peak pressure of the three installation modes are shown in (a) for piston rod side and in (b) for groove side, respectively.
Figure 7. Comparison of installation mode on contact peak pressure. The contact peak pressure of the three installation modes are shown in (a) for piston rod side and in (b) for groove side, respectively.
Applsci 11 02036 g007
Figure 8. Comparison of installation stroke. The contact length in (a) and contact peak pressure in (b) under eight different installation strokes are compared at three positions, including groove side, piston rod side and bottom side.
Figure 8. Comparison of installation stroke. The contact length in (a) and contact peak pressure in (b) under eight different installation strokes are compared at three positions, including groove side, piston rod side and bottom side.
Applsci 11 02036 g008
Figure 9. Comparison of pre-compression rate. The contact peak pressure under nine different pre-compression rates are shown in (a) for piston rod side and in (b) for groove side. Five groups simulation are compared, with five sealing rings of different section diameters coded from 1 to 5. And M represents for zero hydraulic pressure.
Figure 9. Comparison of pre-compression rate. The contact peak pressure under nine different pre-compression rates are shown in (a) for piston rod side and in (b) for groove side. Five groups simulation are compared, with five sealing rings of different section diameters coded from 1 to 5. And M represents for zero hydraulic pressure.
Applsci 11 02036 g009
Figure 10. Comparison of hydraulic pressure. The contact peak pressure under different hydraulic pressures are shown in (a) for piston rod side and in (b) for groove side. Five groups simulation are compared, with five different pre-compression rates coded from A to E, and code 5 represents for the section diameter of the rubber ring (value is 7 mm). And failure zone is plotted to evaluate the sealing performance.
Figure 10. Comparison of hydraulic pressure. The contact peak pressure under different hydraulic pressures are shown in (a) for piston rod side and in (b) for groove side. Five groups simulation are compared, with five different pre-compression rates coded from A to E, and code 5 represents for the section diameter of the rubber ring (value is 7 mm). And failure zone is plotted to evaluate the sealing performance.
Applsci 11 02036 g010
Figure 11. Verification of the proposed butterfly-ring. The comparisons of contact peak pressure between O-ring and butterfly-ring are made at piston rod side in (a) and at groove side in (b), with nine different pre-compression rates. Moreover, the contact peak pressure of O-ring and butterfly-ring under assemble status and working status are compared, respectively.
Figure 11. Verification of the proposed butterfly-ring. The comparisons of contact peak pressure between O-ring and butterfly-ring are made at piston rod side in (a) and at groove side in (b), with nine different pre-compression rates. Moreover, the contact peak pressure of O-ring and butterfly-ring under assemble status and working status are compared, respectively.
Applsci 11 02036 g011
Figure 12. The force between piston rod and rubber ring in the installation mode RA is shown in (a). In the installation mode RL and GL, the interaction force is the same, and is shown in (b).
Figure 12. The force between piston rod and rubber ring in the installation mode RA is shown in (a). In the installation mode RL and GL, the interaction force is the same, and is shown in (b).
Applsci 11 02036 g012
Table 1. Value and code of variable.
Table 1. Value and code of variable.
VariableNumberSeries
Section diameterCode12345
Value (mm)1.82.73.65.37
Pre-compressionCodeABCDEFGHK
Value (%)8910111213141516
Hydraulic pressureCodeMNOPQRSTUVW
Value (MPa)02468101214161820
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, L.; Wei, X. A Novel Structure of Rubber Ring for Hydraulic Buffer Seal Based on Numerical Simulation. Appl. Sci. 2021, 11, 2036. https://doi.org/10.3390/app11052036

AMA Style

Zhang L, Wei X. A Novel Structure of Rubber Ring for Hydraulic Buffer Seal Based on Numerical Simulation. Applied Sciences. 2021; 11(5):2036. https://doi.org/10.3390/app11052036

Chicago/Turabian Style

Zhang, Lin, and Xiaohui Wei. 2021. "A Novel Structure of Rubber Ring for Hydraulic Buffer Seal Based on Numerical Simulation" Applied Sciences 11, no. 5: 2036. https://doi.org/10.3390/app11052036

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