Next Article in Journal
Characterization of the Mechanical Properties, Water Sorption, and Solubility of Antibacterial Copolymers of Quaternary Ammonium Urethane-Dimethacrylates and Triethylene Glycol Dimethacrylate
Next Article in Special Issue
Effect of Impact Velocity and Angle on Impact Wear Behavior of Zr-4 Alloy Cladding Tube
Previous Article in Journal
Effects of Solution Temperature on Tensile Properties of a High γ′ Volume Fraction P/M Superalloy
Previous Article in Special Issue
Three-Body Abrasive Wear Behavior of WC-10Cr3C2-12Ni Coating for Ball Mill Liner Application
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Sealing Effect Improvement Prediction of Flat Rubber Ring in Roller Bit Based on Yeoh_Revised Model

1
School of Petroleum Engineering, China University of Petroleum (East China), Qingdao 266580, China
2
Mechatronic Engineering College, Southwest Petroleum University, Chengdu 610500, China
*
Authors to whom correspondence should be addressed.
Materials 2022, 15(16), 5529; https://doi.org/10.3390/ma15165529
Submission received: 13 June 2022 / Revised: 3 August 2022 / Accepted: 9 August 2022 / Published: 11 August 2022

Abstract

:
In a roller bit, the flat rubber ring (FRR) often needs to apply a certain amount of compression to ensure that its rotation and static sealing surfaces can be stably sealed. For the predicted Mises stress, values smaller than the actual Mises stress due to soft single-axis compression (SAC) stress are predicted by the Yeoh (N = 3) model. To more reasonably predict stress under the static compression of the FRR in the roller bit, the sealing effect of the FRR based on the SAC contact stress and the calculated Mises stress was evaluated by the Yeoh_revised model. Based on the assumption that hydrogenated nitrile-butadiene rubber (HNBR) is isotropic and incompressible, first, we derived the fitting formulas for three types of constitutive models and the Jacobi matrix of the Yeoh_revised model and developed hyperelastic constitutive subroutines. Simultaneously, the accuracy of three models (Yeoh, Yeoh_revised and Ogden) was evaluated by the goodness of fit (R2) to data from three kinds of tensile experiment tests. The highest R2 is 0.9771 with the Yeoh_revised model, which merges the advantages of the other two fitting models and effectively improves the Yeoh model’s soft property of SAC contact stress. Additionally, by measuring on-site FRR wear, the maximum Mises stress on the sealing surface calculated based on the Yeoh_revised model is about twice that of the Yeoh model, and the maximum Mises stress on the rotation contact sealing surface is higher than that on the outside (static sealing) surface, which makes the aging of the rotation surface more severe. Thus, it was demonstrated that, on the premise of ensuring FRR sealing contact stress, the Yeoh_revised model can more reasonably predict the sealing effect of the FRR to more precisely calculate Mises stress than the Yeoh model. This also contributes to FRR structure optimization to prolong the service life of the FRR in the roller bit.

1. Introduction

Generally speaking, the rubber of the FRR used in a roller bit is an initially isotropic and incompressible material. The study of its mechanical constitutive model is mainly divided into two categories: phenomenological theory and thermodynamic statistical theory based on molecular chain networks. Based on the former theory, Mooney, M. [1] and Rivlin, R. S. [2] proposed the classic Mooney–Rivlin (M-R) hyperelastic constitutive model, which represents the strain energy density function with deformation tensor invariant Ii (i = 1, 2, 3). The second-order term of this model is the most widely used hyperelastic constitutive model in current small deformation analyses, with significant errors based on the Neo-Hookean model when fitting experimental data with a tensile ratio greater than 1.4. On the basis of the Mooney–Rivlin model, O.H. Yeoh [3] only considered the effect on the strain energy density function from the first strain invariant (I1) while retaining (I1-3) at less than or equal to the third-order terms to accurately predict the uniaxial and plane tensile test stress–strain relationship, but the “soft” phenomenon was observed compared to experimental data. Ogden, R. W. [4] represented the strain energy density function with the elongation ratio λi (i = 1, 2, 3), which can accurately predict the stress–strain data of uniaxial stretching, plane stretching and equibiaxial stretching experiments simultaneously, but it is not suitable for predicting single-stretching experimental data.
Based on the latter theory, Guth et al. [5] proposed the Neo-Hookean model, called the classical Gaussian statistical mechanics model, which is only suitable for small deformation situations. Lopez-Pamies, O. [6] proposed an evolving porous Neo-Hookean model to characterize the macroscopic response of 2D isotropic porous Neo-Hookean solids with random and particulate microstructures, providing ample motivation to carry out further analyses for more general 3D material systems. Treloar, L. R. G., et al. [7] and Treloar, L. R. G. and G. Riding [8] proposed that the Neo-Hookean model no longer applies when the end of a Gaussian chain exceeds the vector length by 1/3 of the chain length. Guth, E. et al. [5] and Wang, M. C. and E. Guth [9] modified the ideal Gaussian theory under the distribution of the Langevin inverse function, which more accurately simulates the nonlinear elastic characteristics of rubber materials under large deformation, and further proposed a simplified three-chain model. ARRUDA, E. M. and M. C. Boyce [10] proposed an eight-chain molecular model based on the three-chain non-Gaussian network model and showed that the Arruda–Boyce model may be preferable when only uniaxial tensile test data are available. On the previous basis, P. D. Wu and E. van der Giessen [11] established a more accurate full network model, which can fit between the three-chain and eight-chain molecular models. Yang, L. and L. Yang [12] proposed that under certain assumptions, the Gent model [13,14] and the Arruda–Boyce model are consistent. Lopez-Pamies, O. [15] proposed a two-term model, which can be regarded as a generalized Neo-Hookean model, as well as a variant of the Arruda–Boyce (eight-chain) model, and can degenerate into the Yeoh configuration, so it also fits equibiaxial stretching experimental data with the soft phenomenon. Based on this model, Huang, Z.-P. [16] proposed a Gent-Gent hyperelastic constitutive model, which improved the Gent fit. Zhou, L. et al. [17] provided further evidence of the good performance of the equibiaxial deformation prediction of the Gent-Gent model by using it to study the inflation of a circular plane membrane.
Based on the experimental data of L. R. G. TRELOAR [18,19], Xiao-ling, H. et al. [20] fitted five common prediction models for stress–strain test data (Neo-Hookean, Mooney–Rivlin, Yeoh, Ogden and Arruda–Boyce) and employed a model selection strategy for three basic test datasets. The results show that only Ogden (N = 3) could completely fit the experimental data when the three types of experiments were sufficient, avoiding the soft phenomenon of equibiaxial tensile stress fitting. Xue-bing, L. and W. Yin-tao [21] proposed a revised Yeoh hyperelastic material constitutive model that effectively overcomes the “partial soft” properties of the Yeoh model in predicting equibiaxial stretching curves. In a large strain range, the stress–strain relationship of uniaxial, plane and equibiaxial tension compression could be accurately predicted simultaneously, but it was not compared with the Ogden (N = 3) model and needs to be verified by the finite element method (FEM).
Among the constitutive models used in FRR seal analysis, the M-R model is the most extensively used; for instance, Niu, S [22] and Zhang, J. and J. Xie. [23] investigated the sealing performance of an O-ring based on FEM with the M-R model, where the effects of pressure and pre-compression, fluid pressure, the friction coefficient, etc., on the sealing performance analysis were studied. In addition, Liao, B. et al. [24] obtained the time-varying reliability of the O-ring, and Hu, Y. et al. [25] studied the fretting wear of the O-ring, as well as the life prediction, safety and service life of the seal. Yi Zhou and L. Wang [26] proposed that O-ring rubber circles are superior to other cross-sectional rubber rings, while Zhang, H. and J. Zhang [27] reported that the O-ring can be effectively replaced by a D-ring in a static seal. Liang, B. et al. [28] proposed that the material and geometric parameters had a greater influence on the reliability of the rubber O-ring, so Zhang, L. and X. Wei. [29] proposed a new structure of a butterfly rubber ring by analyzing the factors affecting the sealing performance of the rubber ring on the groove side. The Yeoh model is also widely used. Zhou, C. et al. [30,31] studied the sealing performance of a combined sealing structure composed of a rubber D-ring, O-ring and wedge ring used in a high-pressure hydrogen storage vessel by elucidating the swelling mechanisms due to dissolved hydrogen.
The Yeoh_revised constitutive model is the most accurate in fitting experimental data (R2 = 0.9771) compared to Yeoh and Ogden models and is especially suitable for predicting Mises and contact stresses under SAC deformation. In this study, we evaluated the sealing effect of SAC contact stress and Mises stress based on the Yeoh_revised model for a novel FRR in a roller bit under SAC deformation. The results clearly show that the Yeoh_revised model can predict the Mises stress more accurately than the Yeoh model and provides a more reasonable prediction of the aging of the FRR, which further ensures a more stable seal. It also provides specifications for its size optimization in the future.

2. Problem Description and Analysis Method

2.1. FRR Static Compression Analysis in a Roller Bit

FRRs are often used for sealing in roller bits. The sealing effect is related to the influence of the working environment, such as the moving force in the roller bit, down-hole temperature, ambient pressure, etc. [32,33]. During the sealing process, the contact stress of the sealing surface is related to the amount of FRR compression, the fluid pressure difference on both sides of the FRR, and FRR hardness. Generally speaking, higher ambient temperature and internal Mises stress tend to accelerate the aging of rubber [24,25], resulting in rubber hardness increases that make the wear of the sealing surface inconsistent [34]. The higher the hardness of the FRR surface, the slower the wear. Due to inconsistent hardness, a wear gradient of the FRR in the roller bit will occur, and the sealing effect of the FRR will finally deteriorate.
To reduce the difference in wear, this study calculated the contact stress and Mises stress of an FRR in the static sealing process based on the Yeoh_revised model and compared it to the calculation results of the Yeoh model. The aim is to achieve a more accurate prediction of the maximum Mises stress on the FRR under the premise of ensuring a stable seal. Meanwhile, it provides a reference for reasonably designing the size of the FRR structure, prolonging the service life and ensuring a more reliable seal.
The novel structure scheme of the FRR is shown in Figure 1a, in which the center symmetry form and structural parameters are as follows: seal length l = 6.25 mm, height h = 3.0 mm, oblique length b = 2.3 mm, arc radius r = 3.84 mm, oblique horizontal angle e = 11°, arc corresponding to central angle β = 25°, oblique chamfer length α = 0.6 mm, etc.
Figure 1b is the schematic diagram of the FRR force in the roller drill, where the outside surface and rotation surface of the FRR are pre-compressed by a sealed groove and axle journal, respectively. During the actual sealing process, the two ends of the rotation surface profile of the FRR are subjected to fluid pressure (e.g., 20 MPa). Meanwhile, to simplify the calculation, the ambient temperature is assumed to be 120 °C, the friction factor of the sealing contact surface is 0.08, the sealed groove and axle journal are rigid bodies only considering FRR deformation, and the contact, load and constraint of the whole model have axial symmetry. Cloud diagrams of FEM calculation results for contact stress and Mises stress under FRR static compression are shown in Section 3.2. By setting the contact path (Figure 1b), the contact stress and Mises stress curves can be obtained as shown in Section 3.2 too.

2.2. Hyperelastic Experiment and Fitted Model

To determine the hyperelastic constitutive model of HNBR materials, uniaxial tensile (UT), planar tensile (PT) and equivalent tensile (ET) tests [35,36,37] at 120 °C were performed based on the good resistance of HNBR at high temperatures, as shown in Figure 2. In the basic dense elastomer hyperelastic tensile experiment executed by Axel Products, Inc., the parameters were set as follows: strain speed: 0.01 s−1, five times per level; strain levels: 11%, 24%, 30% and 39%. Three UT specimens, three PT specimens and three ET specimens were cut from the provided slabs (slab measurements: 150 mm in length by 150 mm in width, 1.0 to 2.0 mm thick). In the test, the specimens were loaded slowly with between zero force and a user-defined stretch level for 5 loadings and unloadings at up to 4 maximum strain levels so as to examine the initial stress–strain behavior and the “stabilized” stress–strain behavior in each of the maximum strain conditions. In addition, a schematic diagram of the tensile sample force is given in Figure 3. Considering the HNBR as an isotropic incompressible material, the main tensile ratio λi (i = 1,2,3) and the Cauchy stress in related directions (Figure 3), Equation (1) is obtained.
S T : λ 1 = λ , λ 2 = λ 0.5 , λ 3 = λ 0.5 , σ 2 = σ 3 = 0 E T : λ 1 = λ , λ 2 = λ , λ 3 = λ 2 , σ 3 = 0 P T : λ 1 = λ , λ 2 = 1 , λ 3 = λ 1 , σ 3 = 0

2.2.1. Mathematical Stress Formula

HNBR materials are usually regarded as incompressible and isotropic hyper-elastomers, and their mechanical properties can be described by the strain energy density equation. Under isothermal conditions, the strain energy density equation of isotropic super-elastic materials can be represented as in Equation (2) as a function of the three strain invariants of the left/right Cauchy–Green deformation tensor b/C, as defined in Equation (3).
ϕ = ϕ ( I C , I I C , I I I C ) , ϕ = ϕ ( I b , I I b , I I I b )
in which:
I b = I C = t r C = λ 1 2 + λ 2 2 + λ 3 2 I I b = I I C = 1 2 [ ( t r C ) 2 t r C 2 ] = λ 1 2 λ 2 2 + λ 1 2 λ 3 2 + λ 2 2 λ 3 2 I I I b = I I I C = det C = λ 1 2 λ 2 2 λ 3 2
By definition, the relationship between type II Piola–Kirchhoff stress S and strain energy density ϕ (strain energy per unit volume) is expressed as Equation (4) [13].
S = 2 ϕ C
Combining Equation (2) with Equation (3) and deriving Equation (4) by the chain rule of the composite function [21], we derived Equation (5).
S = 2 ϕ C = 2 [ ( ϕ I C + I C ϕ I I C ) I ϕ I I C C + I I I C ϕ I I I C C 1 ]
We consider the relationship between type II Piola–Kirchhoff stress S and Cauchy stress σ as in Equation (6) and combine C = F T F , b = F F T , I b = I C , I I b = I I C , I I I b = I I I C and J = ( I I I b ) 0.5 . Using the Cayley–Hamilton theorem [38] in Equation (7), Equation (6) can be transformed into Equation (8).
σ = 1 J F S F T
b 2 = I b b + I I b I I I I b b 1
σ = 2 ( I I I b ) 1 2 [ ( ϕ I b + Ι C ϕ I I b ) b ϕ I I b b 2 + I I I b ϕ I I I b I ] = 2 ( I I I b ) 1 2 [ ϕ I b b ϕ I I b I I I b b 1 + ( I I b ϕ I I b + I I I b ϕ I I I b ) I ] = P I + ϕ I b b ϕ I I b b 1
Depending on the incompressibility of the rubber material, Ι Ι Ι b = Ι Ι Ι C = 1 , and P is hydrostatic pressure. For the isotropic materials, σ, b and b−1 are coaxial [39]; therefore, in the spindle coordinate system, the principal Cauchy stress in Equation (9a) and the nominal principal stress in Equation (9b) are obtained.
σ i = P + 2 ϕ I b λ i 2 2 ϕ I I b λ i 2
P i = σ i λ i = P λ i + 2 ϕ I b λ i 2 ϕ I I b λ i 3
Eliminating uncertain hydrostatic pressure terms, Equation (10) is as follows.
σ i σ j = 2 ϕ I b ( λ i 2 λ j 2 ) 2 ϕ I I b ( λ i 2 λ j 2 )

2.2.2. The Three Models’ (Yeoh (N = 3), Yeoh_Revised (N = 3) and Ogden (N = 3)) Fitted Stress

Yeoh (N = 3) uses the higher-order term of I1 to correct the strain energy function of the Neo-Hookean model. The modified model is called the Reduced Polynomial model, whose strain energy function is expressed in Equation (11). Based on Drucker stability [40], the material constraint inequalities are defined in Equation (12).
ϕ = C 10 ( I 1 3 ) + C 20 ( I 1 3 ) 2 + C 30 ( I 1 3 ) 3
C 10 > 0 , C 20 < 0 , C 30 > 0
Combining Equations (9b) and (10), the relationship between the nominal stress Pi and the stretching ratio λi of the Yeoh (N = 3) model under three basic deformation modes is obtained, as shown in Equation (13) [38].
P y e o h S T = i = 1 3 2 i C i 0 ( λ 2 + 2 λ 1 3 ) i 1 ( λ λ 2 ) P y e o h E T = i = 1 3 2 i C i 0 ( 2 λ 2 + λ 4 3 ) i 1 ( λ λ 5 ) P y e o h P T = i = 1 3 2 i C i 0 ( λ 2 + λ 2 2 ) i 1 ( λ λ 3 )
Since the IIC term in Equation (2) is completely discarded, the Yeoh model predicts the ET stress with the “soft” phenomenon. Based on the above considerations, by introducing the first part of IIC [20], the following strain energy equation is proposed (see Equation (14)), and based on Drucker stability [40], the material constraint inequality is defined in Equation (15).
ϕ = C 10 ( I 1 3 ) + C 20 ( I 1 3 ) 2 + C 30 ( I 1 3 ) 3 + C 01 ( I 2 3 )
C 10 > 0 , C 20 < 0 , C 30 > 0 , C 01 > 0
Combining Equations (9b) and (10), the relationship between the nominal stress Pi and the stretching ratio λi of the Yeoh_revised (N = 3) model under three basic deformation modes is obtained, as shown in Equation (16).
P y e o h _ r e S T = i = 1 3 2 i C i 0 ( λ 2 + 2 λ 1 3 ) i 1 ( λ λ 2 ) + 2 C 01 ( 1 λ 3 ) P y e o h _ r e E T = i = 1 3 2 i C i 0 ( 2 λ 2 + λ 4 3 ) i 1 ( λ λ 5 ) + 2 C 01 ( λ 3 λ 3 ) P y e o h _ r e P T = i = 1 3 2 i C i 0 ( λ 2 + λ 2 2 ) i 1 ( λ λ 3 ) + 2 C 01 ( λ λ 3 )
The strain energy function of the Ogden model is shown in Equation (17) [4,38], where μi and αi are arbitrary constants (which can be non-integers). Based on Drucker stability [40], material constraint inequalities are shown in Equation (18).
ϕ = i = 1 N μ i α i ( λ 1 α i + λ 1 α i + λ 1 α i 3 )
i = 1 N μ i α i > 0
Combining Equations (9b) and (10), the relationship between the nominal stress Pi and the stretching ratio λi of the Ogden (N = 3) model under three basic deformation modes is obtained as shown in Equation (19).
P o g d e n S T = 2 i = 1 3 μ i α i ( λ α i 1 λ 1 2 α i 1 ) P o g d e n E T = 2 i = 1 3 μ i α i ( λ α i 1 λ 2 α i 1 ) P o g d e n P T = 2 i = 1 3 μ i α i ( λ α i 1 λ α i 1 )

2.2.3. Model Goodness of Fit

Using the hyperelastic constitutive fitted model obtained in Section 2.2.2, least squares fitting was performed on experimental data [18,32]. The model goodness of fit R2 was introduced to evaluate the fitting quality, and the deviation square sum (SSdev) and the residual square sum (SSres) were calculated.
R 2 = 1 S S r e s S S d e v
S S r e s = i = 1 N ( P i P i ) 2
S S d e v = i = 1 N ( P i P ¯ i ) 2
where Pi is the experimental value; P i ¯ is the average of Pi; P i ^ is the model fit value; and N is the number of test data points involved in fitting. A larger goodness-of-fit value, R2, indicates the higher overall goodness of the model fit. The Yeoh fitting formula in Equation (13), Yeoh_revised fitting formula in Equation (16) and Ogden fitting formula in Equation (19) were fitted to three sets of tensile experimental data. The fitting results and residual analysis are shown in Figures 6 and 7, respectively.

2.3. Validation of the Fitted Constitutive Parameters in FEM

In engineering, the vast majority of dense HNBR materials are often considered mechanically incompressible during the deformation process [41,42], which still satisfies the accuracy requirements of the solution. In general, the hyperelastic Jacobi matrix of incompressible HNBR materials is often solved in the form of the stress–strain rate; S , τ is not objective [39], so C tan cannot be used for the hyperelastic Jacobi matrix of rubber materials [39]. Since τ o and τ are objective, the Jacobi matrix that can be used to solve the constitutive model according to its objectivity is shown in Equations (21a) and (21b). Finally, the Jacobi matrix Cijkl of the material is derived from the variation in the Jaumann-Zaremba rate of the Kirchhoff stress tensor [39].
S = C t a n : E τ = τ l τ τ l T F 1 τ F T = C t a n : F T D F τ = ( F ¯ F : C t a n : F T ¯ F T ) : D = L : D
τ = τ + D τ + τ D = L : D D τ + τ D = 2 H : D L = L + 2 H C i j k l = L ^ i j k l / J = ( L i j k l + 2 H i j k l ) / J
where C tan is known as the material tangent elasticity tensor, τ is the derivative of the Kirchhoff stress tensor with time, τ is the Oldroyd rate of the Kirchhoff stress tensor, τ o is the Jaumann-Zaremba rate of the Kirchhoff stress tensor, L or L ^ is the spatial tangent elasticity tensor in the current configuration, H is the symmetric part of L ^ , l is the velocity deformation gradient tensor, and D is the symmetric part of l.

2.3.1. Calculation of Yeoh_Revised Jacobi Matrix (Incompressible)

In order to verify the correctness of the HNBR hyperelastic configuration model, this section first establishes the Jacobi matrix of the Neo-Hookean model and compares the results with the Neo-Hookean Jacobi matrix given by the help document in ABAQUS [43]. Helmholtz free energy per unit reference volume (ϕ) of the Neo-Hookean model is shown in Equation (22) [44,45].
ϕ = ϕ ˜ + ϕ v o l = C 10 ( I C ˜ 3 ) + 1 D 1 ( J 1 ) 2
where ϕ ˜ , ϕ v o l represent the shear and dilatation parts of Helmholtz free energy, respectively, and Di characterizes the compressibility of HNBR.
According to the definitions in Equations (4) and (5), the second type of Piola–Kirchhoff stress (S) is shown in Equation (23a), where S ˜ , S v o l represent the shear and dilatation parts of Piola–Kirchhoff stress (S), S ¯ ˜ represents the stress (S) defined in the pure dilatation configuration (Figure 4), and P is the fourth-order projection tensor.
S = 2 ϕ ( C ) C = S ˜ + S v o l
S v o l = J ϕ v o l ( J ) J C 1 = 2 D 1 J ( J 1 ) C 1 S ˜ = 2 ϕ ˜ ( C ˜ ) C = J 2 3 [ C ˜ C : 2 ϕ ˜ ( C ˜ ) C ˜ ] = J 2 3 P : S ¯ ˜ = J 2 3 S ¯ ˜ d e v
S ¯ ˜ = 2 ϕ ˜ ( C ˜ ) ( C ˜ ) ,   P = C ˜ C = I I 1 3 C 1 C
Substituting Equations (22), (23b) and (23c) into Equation (23a), the specific expression of S is obtained in Equation (24).
S = S ˜ + S v o l = 2 C 10 J 2 3 [ I t r C 3 C 1 ] + 2 D 1 J ( J 1 ) C 1
To determine the Neo-Hookean Jacobi matrix, Equation (24) is substituted into Equation (21a) to obtain the material tangent elasticity tensor ( C tan ), as shown in Equation (25), and the spatial tangent elasticity tensor ( L ) is shown in Equation (26).
C tan = 2 S C = 4 3 C 10 J 2 3 [ C 1 I + I C 1 ] + [ 4 9 C 10 J 2 3 t r C + 2 D 1 ( 2 J 1 ) J ] C 1 C 1 + [ 2 3 C 10 J 2 3 t r C 2 D 1 ( J 2 J ) ] [ C 1 ¯ C 1 + C 1 ¯ C 1 ]
L = F ¯ F : C tan : F T ¯ F T = 4 3 C 10 J 2 3 [ I b + b I ] + [ 4 9 C 10 t r b J 2 3 + 2 D 1 ( 2 J 1 ) J ] I I + [ 2 3 C 10 t r b J 2 3 2 D 1 ( J 2 J ) ] [ I ¯ I + I ¯ I ]
In fact, the material Jacobi matrix is obtained from the variation in the Jaumann-Zaremba rate of the Kirchhoff stress tensor ( τ o ). According to the relationship between τ o and τ (see Equation (21b)), Hijkl and L ^ i j k l component expressions are shown in Equations (27) and (28), respectively. The detailed derivation process of Hijkl is shown in Appendix A.
2 H i j k l = C 10 ( δ i k b ˜ j l + δ i l b ˜ j k + b ˜ i k δ j l + b ˜ i l δ j k ) [ 2 3 C 10 t r b ˜ 3 2 D 1 ( J 2 J ) ] ( δ i k δ j l + δ i l δ j k )
L i j k l = L i j k l + 2 H i j k l = 4 3 C 10 [ δ i j b ˜ k l + b ˜ i j δ k l ] + [ 4 9 C 10 t r b ˜ + 2 D 1 ( 2 J 1 ) J ] δ i j δ k l + C 10 ( δ i k b ˜ j l + δ i l b ˜ j k + b ˜ i k δ j l + b ˜ i l δ j k )
Finally, the Neo-Hookean Jacobi matrix component Cijkl of rubber material is shown in Equation (29). The result is consistent with the help document given in the ABAQUS help file [46].
C i j k l = L i j k l / J = 4 3 J C 10 [ δ i j b ˜ k l + b ˜ i j δ k l ] + [ 4 9 J C 10 t r b ˜ + 2 D 1 ( 2 J 1 ) ] δ i j δ k l + C 10 J ( δ i k b ˜ j l + δ i l b ˜ j k + b ˜ i k δ j l + b ˜ i l δ j k )
Likewise, based on Equation (14), the strain energy density equation for the Yeoh_revised (incompressible) model is obtained [38] (see Equation (30)).
ϕ = ϕ ˜ + ϕ v o l = i = 1 3 C i 0 ( I C 3 ) i + C 01 ( I I C 3 ) + γ ( J 1 )
where C ˜ = C , and γ is the introduced Lagrange multiplier and satisfies γ = P = σ k k / 3 .
Equation (30) is substituted into Equation (23a), and after the calculation of Equation (31), Equation (32) is obtained.
S v o l = J ϕ v o l ( J ) J C 1 = J γ C 1 S ˜ = 2 ϕ ˜ ( C ˜ ) C = 2 ( ϕ ˜ ( C ˜ ) I C + ϕ ˜ ( C ˜ ) I I C I C ) I ϕ ˜ ( C ˜ ) I I C C = i = 1 3 2 i C i 0 ( I C 3 ) i 1 I + 2 C 01 ( I C I C )
S = S ˜ + S v o l = i = 1 3 2 i C i 0 ( I C 3 ) i 1 I + 2 C 01 ( I C I C ) + J γ C 1
Similar to the calculation of Equations (25)–(28), the Yeoh_revised model (incompressible) Jacobi matrix component Cijkl of HNBR material is shown in Equation (33), and the detailed derivation process is shown in Appendix B.
C i j k l = L i j k l / J = 1 J [ 8 C 20 + 24 C 30 t r b 72 C 30 4 C 01 ] b i j b k l 4 C 01 J b i k b j l + 1 J [ J γ ] δ i j δ k l + 1 J [ C 10 + 2 C 20 ( t r b 3 ) + 3 C 30 ( t r b 3 ) 2 + C 01 t r b ] ( δ i k b j l + δ i l b j k + b i k δ j l + b i l δ j k ) 1 J C 01 ( δ i k b j p b p l + δ i l b j p b p k + b i p b p k δ j l + b i p b p l δ j k )

2.3.2. Notes on Developing Subroutines

In this study, a subroutine program based on Equation (33) was used to verify the superiority of the Yeoh_revised model (incompressible) in predicting SAC by FEM, and it was compared with the Yeoh and Ogden models. Since the Jacobi matrix expressed in Equation (33) is derived based on the premise of incompressibility, UMAT or UHYPER subprograms can be directly developed [47].

2.3.3. Verification of the Equivalent CAE Simulation Design for the ET (SAC) Test

On the premise that HNBR material volume is incompressible, combined with experimental data obtained at constant temperature (such as 120 °C), different λi values under ET boundary conditions were converted to different compressive axial displacements (uzi) based on Equation (1) (see Table 1). The contact stress between the picked mesh element and pressure plate was extracted, as shown in Figure 5, by using an axisymmetric model for calculations; the finite element mesh type of rubber material is CAX4H, and the number of CAX4H is 300. The pressure plate is designated as an analytical rigid body, its mesh type is RAX2, and the pressure plate and rubber are set to be in frictionless contact (equivalent to ET). The contact stresses under different λi are plotted together in Figure 8, which is based on the three constitutive models obtained using Equations (13), (16) and (19), respectively. In addition, the Yeoh_revised Jacobi matrix was derived on the basis of Equation (33), and the other two Jacobi matrix models usually exist in popular FEM software.

3. Results and Discussion

3.1. Fitting Results Analysis

In Figure 6 and Figure 7, based on the three types (UT, ET and PT) of tensile experimental data of HNBR at 120 °C, the Yeoh_revised (N = 3) constitutive model’s overall goodness of fit (R2 = 0.9771) is the most accurate compared with the Yeoh (N = 3) constitutive and Ogden (N = 3) models, and the maximum negative fitting residuals are: −0.3782, −0.155, −0.2141, respectively; the Yeoh_revised model is improved compared to the Yeoh constitutive model, which predicts ET (or SAC) stress with the softness phenomenon. This is because the ET experimental fitting formula holds one more item 2 C 01 ( λ 3 λ 3 ) , which gradually increases as λ increases, than the Yeoh model in Equation (16). The obtained curve in Equation (16) shows a more obvious “upturning” phenomenon than that in Equation (13), which indicates better agreement with ET experimental data. Compared with the Ogden (N = 3) model, the Yeoh_revised (N = 3) model achieves a more accurate fit with fewer parameters (i.e., four), so it proves to be the most suitable for stress prediction in SAC deformation. Meanwhile, the Yeoh_revised (N = 3) model still has the smallest fluctuation in UT and PT tensile values, as shown in Figure 7b,c. Thus, in general, the Yeoh_revised (N = 3) model possesses the best fitting effect and merges the advantages of the other two fitting models. The fitting coefficient values of the three models, the deviation square sum (SSdev) and the goodness of fit (R2) are shown in Table 2.

3.2. Equivalent FEM Verification Results

In Figure 8, when the tensile ratio (λ) increases from 1 to 1.10 at an experimental temperature of 120 °C, the experimental values are consistently higher than the CAE-calculated values of the Yeoh_revised, Yeoh and Ogden models. Nevertheless, the three constitutive models produce results with different size degrees; in general, Yeoh_revised is the smallest, Yeoh is the largest, and Ogden is in the middle. For example, when λ is 1.10, the Yeoh_revised CAE calculation model is 10.17% smaller, while Yeoh_revised and Ogden CAE calculation models are 10.17% and 14.16% smaller, respectively, which again indicates that the Yeoh_revised model will significantly improve the soft SAC stress based on the Yeoh model. For other stretch ratios, the comparison between the CAE-calculated values of the three constitutive models and experiments is shown in Table 3.
Meanwhile, in Figure 8, it is clear that as λ increases from 1 to 1.10, the numerical values from the experiment and the CAE-calculated values of Yeoh_revised and Ogden are all linearly distributed, while the Yeoh CAE-calculated values are nonlinearly distributed, and their deviation from experimental values rises as λ increases. Therefore, the Yeoh model is seldom suitable for predicting SAC stress under large deformation.
In order to compare the difference between the experimentally fitted values (i.e., Equations (13), (16) and (19)) and the CAE-calculated values (i.e., Equation (33)) for the same λ, the fitted values of experimental SAC (ET) data were collated together with the CAE-calculated values of the three constitutive models in Table 4. The results demonstrate that the maximum deviation between the fitted values of each constitutive model and the corresponding CAE-calculated values never exceeds ±0.5% at 120 °C, which, importantly, meets the requirement for engineering accuracy, so the accuracy of the fitting parameters of each constitutive model obtained by the least squares method in this paper is proved.

3.3. FEM Results of FRR under Static Compression

The cloud diagram of contact stress FEM calculation results under FRR static compression is shown in Figure 9. For a more intuitive comparison, the contact stresses calculated based on two constitutive models in the contact path (Figure 1b) is extracted and plotted in Figure 10, where the static pre-compression contact stresses are basically the same distribution, and the difference between the contact stresses after being squeezed by fluid is also subtle. In comparison, based on the Yeoh_revised model, the contact stress value on the rotation surface is slightly larger than the value calculated by the Yeoh model between points A and B in the contact path, while the contact stress value on the outside surface between points D and E shows the opposite pattern. After being squeezed by fluid, the contact stress of the two contact surfaces increases above the fluid pressure (i.e., 20 MPa) in the cavity. Meanwhile, the contact stress after fluid extrusion based on the two constitutive models is basically the same, which indicates that the FRR structure with this structural parameter set would ensure that it is sealed. The Mises stress expression is shown in Equation (34).
σ = 1 2 [ ( σ 1 σ 2 ) 2 + ( σ 2 σ 3 ) 2 + ( σ 3 σ 1 ) 2 ]
Under static compression, the maximum Mises stress calculated based on the two constitutive models is located in the F, G and C areas of the FRR section, as shown in Figure 11a,c, and only the maximum Mises stress values are different. According to Equations (20a) and (34), the Yeoh_revised model predicts principal stress σi (i = 1, 2, 3) values more accurately in each direction, so it also desired the value of Mises stress more precise. For instance, the maximum Mises stress value obtained based on the Yeoh_revised model is 1.437 MPa, which is greater than the Yeoh model value of 1.413 MPa.
After being squeezed by fluid, the maximum Mises stress also occurs in areas where the two sealing surfaces are in contact with the fluid, i.e., A, B, D and E. It is also found that the cloud graph in Figure 11d is not consistent with Figure 11b; based on the Yeoh_revised model, the maximum Mises stress after fluid extrusion is mainly distributed in the contact areas between the two sealing contact surfaces and the liquid cavity (i.e., A, B, D and E), while the Yeoh model hardly reflects this distribution law. Extracting the Mises stress values on the contact path (Figure 1b), as shown in Figure 12, indicates that Mises stress values at the contact areas (i.e., A, B) are basically the same. Owing to the symmetrical deformation of the FRR, the contact areas are at the junction of the rotation sealing surface and fluid of both ends, and the two contact areas (D, E) on the outer sealing surface are also identical. Meanwhile, based on the two constitutive models, Mises stress concentration occurs in the four regions (A, B, D and E) after being squeezed by fluid (i.e., the FRR sealing process), and the Mises stress obtained by the Yeoh_revised model is twice the value obtained by the Yeoh model; the larger the value, the easier it is to age and harden in the four regions. Simultaneously, Mises stress in the contact areas (i.e., A and B) at the junction of the rotating surface and fluid is the maximum on the contact path, which indicates that this area is the most fragile part of the entire FRR structure and is likely to cause the sealing effect of the FRR to deteriorate [48].

4. FRR Field Application

Taking the well named MaHW6415 as an example, the lithology is gray mudstone, the estimated ambient temperature is 90 °C, and the pump pressure is 20 MPa (fluid pressure difference: ±0.5 MPa). After the total number of rotations of the roller bit reached about 980,000 cycles, the FRR was extracted, its appearance was checked, and the wear condition was analyzed, as shown in Figure 13.
There are worn lines on the rotation surface and outside surface in Figure 13b,d. The larger the distance between the center line of the stress concentration area (i.e., thick dotted line) and the wear gradient boundary line (i.e., thin dotted line), the larger the hardened area of the rubber contact surface, the more severe the aging of the rubber contact surface, and the worse the sealing effect. Meanwhile, the distance between the thick dotted line and the thin dotted line of the cross-sectional stress concentration area in Figure 13d is higher than that in Figure 13b, which demonstrates that the maximum Mises stress on the rotation sealing surface is higher than that on the outside sealing surface. Consequently, the spacing distance will be reduced by optimizing the FRR structural parameters to reduce the maximum Mises stress value. The schematic diagrams and test results of rotation and outside surfaces are shown in Figure 14 and Table 5, respectively.
In Table 5, the inner-diameter wear value of the rotation surface before and after use is 0.52 mm; the maximum wear amount of the section length is about 0.30 mm at 90° or 180°, respectively, so eccentric wear occurs, which shows that the FRR is not completely antisymmetrically deformed inside the roller drill. The hardness of the FRR rotation surface is hardened by 11HA before and after use, while the outside surface is only hardened by 4HA, which also proves that the maximum Mises stress on the rotation sealing surface is higher than that on the outside sealing surface in actual operation.

5. Conclusions

This paper mainly focuses on the SAC contact stress and Mises stress of the Yeoh_revised model compared to the Yeoh model, which was also applied to evaluate the sealing effect of an FRR in a roller bit under SAC deformation. The results clearly show that the Yeoh_revised model can predict the Mises stress more accurately than the Yeoh model and provides theoretical guidance for mitigating the aging of FRRs and ensuring a more stable seal. It also further provides specifications for its size optimization. The following important conclusions are drawn from this research:
  • Comparing the fitted values with the FEM-calculated data of three constitutive models, it is demonstrated that the maximum deviation between the fitted value of each constitutive model and the corresponding CAE-calculated value never exceeds ±0.5% at 120 °C, which proves the accuracy of the fitting values of the parameters of each constitutive model obtained through the least squares method in this paper.
  • Compared with the Yeoh model, the Yeoh_revised and Ogden models both address the soft phenomenon encountered by the Yeoh constitutive model when predicting stress in ET (SAC) tensile tests. Moreover, the Yeoh_revised model shows the greatest improvement, with an R2 up to 0.9771, in fitting the experimental values, and its maximum underestimation is reduced to half of that of the Yeoh model.
  • The Yeoh_revised constitutive model is the most accurate in fitting the experimental data. Compared with Ogden, it achieves more accurate fitting with fewer parameters (i.e., four), while the number of fitted parameters (i.e., six) is higher in the Ogden model; therefore, it is more suitable for Mises stress analysis of an FRR in a roller bit under SAC deformation.
  • On the premise of ensuring the stability of sealing FRR contact stress, the maximum Mises stress obtained with the Yeoh_revised model is 1.437 MPa greater than the Yeoh model’s value of 1.413 MPa before FRR extrusion by fluid. The Yeoh_revised model is more accurate in predicting Mises stress. In the sealing process (i.e., after FRR extrusion by fluid), the Mises stress obtained with the Yeoh_revised model is twice the value obtained with the Yeoh model, which provides a more reasonable prediction for reducing FRR aging and further ensuring a more stable seal. It also provides specifications for its size optimization in the future.

Author Contributions

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

Funding

This work was supported by the Sichuan Regional Innovation Cooperation Project (NO. 2021YFQ0041).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no conflict of interest.

Appendix A

The derivation process of Equation (27) is as follows:
I I i j k l s y m = 0.5 ( I I i j k l + I I ¯ i j k l ) = 0.5 ( δ i k δ j l + δ i l δ j k )
where I I = I ¯ I , I I ¯ = I ¯ I is a minor symmetric fourth-order unit tensor [38].
According to the symmetry of D:
D i p = I I i p k l s y m D k l = 0.5 ( δ i k δ p l + δ i l δ p k ) D k l D p j = I I p j k l s y m D k l = 0.5 ( δ p k δ j l + δ p l δ j k ) D k l
Substitute in:
D i p b ˜ p j = 0.5 ( δ i k b ˜ j l + δ i l b ˜ j k ) D k l b ˜ i p D p j = 0.5 ( b ˜ i k δ j l + b ˜ i l δ j k ) D k l
D i p I p j = 0.5 ( δ i k δ j l + δ i l δ j k ) D k l I i p D p j = 0.5 ( δ i k δ j l + δ i l δ j k ) D k l
D i p b ˜ p q b ˜ q j = 0.5 ( δ i k b ˜ j p b ˜ p l + δ i l b ˜ j p b ˜ p k ) D k l b ˜ i p b ˜ p q D q j = 0.5 ( b ˜ i p b ˜ p k δ j l + b ˜ i p b ˜ p l δ j k ) D k l
According to Equation (21b), Equations (A6) and (A7) can be obtained:
τ = J σ = F S F T = 2 C 10 J 2 3 [ b t r b 3 I ] + 2 D 1 J ( J 1 ) I
2 H : D = 2 C 10 ( D b ˜ + b ˜ D ) [ 2 C 10 t r b ˜ 3 2 D 1 ( J 2 J ) ] ( D I + I D )
Substituting Equations (A3) and (A4) into Equation (A7), Equation (27) can be obtained.

Appendix B

Equation (32) can be redefined as Equation (A8):
S = S ˜ + S v o l = i = 1 3 2 i C i 0 ( I C 3 ) i 1 I + 2 C 01 ( I C I C ) + J γ C 1 = S C 10 + S C 20 + S C 30 + S C 01 + S v o l
where:
S C 10 = 2 C 10 I
S C 20 = 4 C 20 ( I C 3 ) I
S C 30 = 6 C 30 ( I C 3 ) 2 I
S C 01 = 2 C 01 ( I C I C )
S v o l = J γ C 1
According to the definition in Equation (25), Equation (A14) is obtained from Equation (A9):
C C 10 tan = 2 S C 10 C = 0
C C 20 tan = 2 S C 20 C = 8 C 20 I I
C C 30 tan = 2 S C 30 C = ( 24 C 30 I C 72 C 30 ) I I
C C 01 tan = 2 S C 01 C = 4 C 01 I I 4 C 01 I ¯ I
C v o l tan = 2 S v o l C = γ J C 1 C 1 γ J ( C 1 ¯ C 1 + C 1 ¯ C 1 )
Combining Equations (A14)–(A18), Equation (A19) is obtained:
C tan = ( 8 C 20 + 24 C 30 I C 72 C 30 + 4 C 01 ) I I 4 C 01 I ¯ I + γ J C 1 C 1 γ J ( C 1 ¯ C 1 + C 1 ¯ C 1 )
According to Equation (21a), Equation (A19) becomes Equation (A20), which is the spatial tangent elasticity tensor (L) of the Yeoh_revised model (incompressible).
L = ( 8 C 20 + 24 C 30 I b 72 C 30 + 4 C 01 ) b b 4 C 01 b ¯ b + γ J I I γ J ( I ¯ I + I ¯ I )
Substituting Equations (A3)–(A5) into Equation (A20), the fractional formula of L of the Yeoh_revised model (incompressible) can be obtained, as shown in Equation (A21).
L i j k l = ( 8 C 20 + 24 C 30 I b 72 C 30 + 4 C 01 ) b i j b k l 4 C 01 b i k b j l + γ J δ i j δ k l γ J ( δ i k δ j l + δ i l δ j k )
Similarly, on the basis of Equation (32), according to Equation (21b), Equations (A22) and (A23) can be obtained.
τ = J σ = F S F T = i = 1 3 2 i C i 0 ( I b 3 ) i 1 + 2 C 01 I b b 2 C 01 b b + J γ I
2 H : D = i = 1 3 2 i C i 0 ( I b 3 ) i 1 + 2 C 01 I b ( D b + b D ) 2 C 01 ( D b b + b b D ) + J γ ( D I + I D )
Substituting Equations (A3)–(A5) into Equation (A23), Equation (A24) can be obtained.
2 H i j k l = i = 1 3 i C i 0 ( I b 3 ) i 1 + C 01 I b ( δ i k b j l + δ i l b j k + b i k δ j l + b i l δ j k ) C 01 ( δ i k b j p b p l + δ i l b j p b p k + b i p b p k δ j l + b i p b p l δ j k ) + J γ ( δ i k δ j l + δ i l δ j k )
Combining Equations (A21) and (A24), the fractional formula of L ^ of the Yeoh_revised model (incompressible) can be obtained, as shown in Equation (A25).
L ^ i j k l = L i j k l + 2 H i j k l = ( 8 C 20 + 24 C 30 I b 72 C 30 + 4 C 01 ) b i j b k l C 01 ( δ i k b j p b p l + δ i l b j p b p k + b i p b p k δ j l + b i p b p l δ j k ) 4 C 01 b i k b j l + γ J δ i j δ k l + i = 1 3 i C i 0 ( I b 3 ) i 1 + C 01 I b ( δ i k b j l + δ i l b j k + b i k δ j l + b i l δ j k )
Finally, on the basis of Equation (A25), the Yeoh_revised model (incompressible) Jacobi matrix component Cijkl is obtained, as shown in Equation (33).

References

  1. Mooney, M. A Theory of Large Elastic Deformation. J. Appl. Phys. 1940, 11, 582–592. [Google Scholar] [CrossRef]
  2. Rivlin, R.S. Large Elastic Deformations of Isotropic Materials. IV. Further Developments of the General Theory. Philos. Trans. R. Soc. London Ser. A Math. Phys. Sci. 1948, 241, 379–397. [Google Scholar]
  3. Yeoh, O.H. Some forms of the strain energy function for rubber. Rubber Chem. Technol. 1993, 66, 754–771. [Google Scholar] [CrossRef]
  4. Ogden, R.W. Large deformation isotropic elasticity—On the correlation of theory and experiment for incompressible rubberlike solids. Proc. R. Soc. Lond. A. 1972, 326, 565–584. [Google Scholar] [CrossRef]
  5. Guth, E.; Wack, P.E.; Anthony, R.L. Significance of the Equation of State for Rubber. J. Appl. Phys. 1946, 17, 347–351. [Google Scholar] [CrossRef]
  6. Lopez-Pamies, O.; Idiart, M.I. An Exact Result for the Macroscopic Response of Porous Neo-hookean Solids. J. Elast. 2009, 95, 99–105. [Google Scholar] [CrossRef]
  7. Treloar, L.R.G.; Hopkins, H.G.; Rivlin, R.S.; Ball, J.M. The Mechanics of Rubber Elasticity. Proc. R. Soc. A Math. Phys. Eng. Sci. 1976, 351, 301–330. [Google Scholar] [CrossRef]
  8. Treloar, L.R.G.; Riding, G. A non-gaussian theory for rubber in biaxial strain. Proc. R. Soc. A Math. Phys. Eng. Sci. 1979, 369, 261–280. [Google Scholar]
  9. Wang, M.C.; Guth, E. Statistical Theory of Networks of Non-Gaussian Flexible Chains. J. Chem. Phys. 1952, 20, 1144–1157. [Google Scholar] [CrossRef]
  10. Arruda, E.M.; Boyce, M.C. A three-dimensional constitutive model for the large stretch behavior of rubber elastic materials. J. Mech. Phys. Solids. 1993, 41, 389–412. [Google Scholar] [CrossRef]
  11. Wu, P.; van der Giessen, E. On Improved 3-D Non-Gaussian Network Models for Rubber Elasticity. Mech. Res. Commun. 1992, 19, 427–433. [Google Scholar] [CrossRef]
  12. Yang, L.; Yang, L. Note on Gent’s Hyperelastic Model. Rubber Chem. Technol. 2018, 91, 296–301. [Google Scholar] [CrossRef]
  13. Gent, A.N. A New Constitutive Relation for Rubber. Rubber Chem. Technol. 1996, 69, 59–61. [Google Scholar] [CrossRef]
  14. Gent, A.N. Elastic instabilities in rubber. Rubber Chem. Technol. 2005, 40, 165–175. [Google Scholar] [CrossRef]
  15. Lopez-Pamies, O. A new-based hyperelastic model for rubber elastic materials. Comptes Rendus Mécanique 2010, 338, 3–11. [Google Scholar] [CrossRef]
  16. Huang, Z.-P. A Novel Constitutive Formulation for Rubberlike Materials in Thermoelasticity. J. Appl. Mech. 2014, 81, 041013. [Google Scholar] [CrossRef]
  17. Zhou, L.; Wang, S.; Li, L.; Fu, Y. An evaluation of the Gent and Gent-Gent material models using inflation of a plane membrane. Int. J. Mech. Sci. 2018, 146, 39–48. [Google Scholar] [CrossRef]
  18. Hu, X.L.; Liu, X.; Li, M.; Luo, W.B. Selection strategies of hyperelastic constitutive models for carbon black filled rubber. Eng. Mech. 2014, 31, 34–42. (In Chinese) [Google Scholar]
  19. Treloar, L.R.G. Stressstrain data for vulcanised rubber under varioustypes of deformation. J. Chem. Soc. 1943, 40, 59–70. [Google Scholar]
  20. Treloar, L.R.G. The Elasticity of a Network of Long-Chain Molecules.3. J. Chem. Phys. 1945, 42, 84–94. [Google Scholar]
  21. Li, X.; Wei, Y. An improved Yeoh constitutive model for hyperelastic material. Eng. Mech. 2016, 33, 38–43. (In Chinese) [Google Scholar] [CrossRef]
  22. Niu, S. Sealing Performance Analysis of Rubber O-ring in Static Seal Based on FEM. Int. J. Eng. Adv. Res. Technol. 2015, 1, 32–34. [Google Scholar]
  23. Zhang, J.; Xie, J. Investigation of Static and Dynamic Seal Performances of a Rubber O-Ring. J. Tribol. 2018, 140, 042202. [Google Scholar] [CrossRef]
  24. Liao, B.; Sun, B.; Yan, M.; Ren, Y.; Zhang, W.; Zhou, K. Time-Variant Reliability Analysis for Rubber O-Ring Seal Considering Both Material Degradation and Random Load. Materials 2017, 10, 1211. [Google Scholar] [CrossRef]
  25. Hu, Y.; Han, C.; Zhang, J.; Luo, Z. Fretting Wear of Rubber Sealing Ring Caused by Fluid Pressure Fluctuation. Mechanics 2021, 27, 321–326. [Google Scholar] [CrossRef]
  26. Zhou, Y.; Wang, L. Study on Optimum Structural Design of Roller Bit Bearing Double Rubber Ring Seal. Int. J. Sci. Res. 2016, 5, 298–303. [Google Scholar]
  27. 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]
  28. Liang, B.; Yang, X.; Wang, Z.; Su, X.; Liao, B.; Ren, Y.; Sun, B. Influence of Randomness in Rubber Materials Parameters on the Reliability of Rubber O-Ring Seal. Materials 2019, 12, 1566. [Google Scholar] [CrossRef]
  29. Zhang, L.; Wei, X. A Novel Structure of Rubber Ring for Hydraulic Buffer Seal Based on Numerical Simulation. Appl. Sci. 2021, 11, 2036. [Google Scholar] [CrossRef]
  30. 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]
  31. Zhou, C.; Zheng, J.; Gu, C.; Zhao, Y.; Liu, P. Sealing performance analysis of rubber O-ring in high-pressure gaseous hydrogen based on finite element method. Int. J. Hydrog. Energy 2017, 42, 11996–12004. [Google Scholar] [CrossRef]
  32. Alacqua, S.; Capretti, G.; Perosino, A.; Veca, A. Sealing Ring Interposed between the Block and the Cylinder Head of an Internal Combustion Engine having a Composite Structure. U.S. Patent 8100410B2, 2 April 2013. [Google Scholar]
  33. Szczypinski-Sala, W.; Lubas, J. Tribological Characteristic of a Ring Seal with Graphite Filler. Materials 2020, 13, 311. [Google Scholar] [CrossRef]
  34. Archard, J.F. Contact and Rubbing of Flat Surfaces. J. Appl. Phys. 1953, 24, 981–988. [Google Scholar] [CrossRef]
  35. Sasso, M.; Palmieri, G.; Chiappini, G.; Amodio, D. Characterization of hyperelastic rubber-like materials by biaxial and uniaxial stretching tests based on optical methods. Polym. Test. 2008, 27, 995–1004. [Google Scholar] [CrossRef]
  36. Ghoreishy, M.H.R. Determination of the parameters of the Prony series in hyper-viscoelastic material models using the finite element method. Mater. Des. 2012, 35, 791–797. [Google Scholar] [CrossRef]
  37. Lev, Y.; Faye, A.; Volokh, K.Y. Thermoelastic deformation and failure of rubberlike materials. J. Mech. Phys. Solids 2019, 122, 538–554. [Google Scholar] [CrossRef]
  38. Abaqus 6.14 Theory Guide. Fitting of Hyperelastic and Hyperfoam Constants. 2014. Available online: http://130.149.89.49:2080/v6.14/books/stm/default.htm (accessed on 10 September 2020).
  39. Chaves, E.W.V. Notes on Continuum Mechanics. In International Center for Numerical Methods in Engineering (CIMNE); Springer: Barcelona, Spain, 2013. [Google Scholar]
  40. Abaqus 6.14 Theory Guide. Hyperelastic Material Behavior. 2014. Available online: http://130.149.89.49:2080/v6.14/books/stm/default.htm (accessed on 10 September 2020).
  41. Chester, S.A.; Anand, L. A coupled theory of fluid permeation and large deformations for elastomeric materials. J. Mech. Phys. Solids 2010, 58, 1879–1906. [Google Scholar] [CrossRef]
  42. Chester, S.A.; Anand, L. A thermo-mechanically coupled theory for fluid permeation in elastomeric materials: Application to thermally responsive gels. J. Mech. Phys. Solids. 2011, 59, 1978–2006. [Google Scholar] [CrossRef]
  43. Chadwick, P. Thermo-Mechanics of Rubberlike Materials. Philos. Trans. R. Society. A Math. Phys. Eng. 1974, 276, 371–403. [Google Scholar]
  44. Zeng, N.; Haslach, H.W. Thermoelastic Generalization of Isothermal Elastic Constitutive Models for Rubber. Rubber Chem. Technol. 1996, 69, 313–324. [Google Scholar] [CrossRef]
  45. Heinrich, G.; Kaliske, M.; Klüppel, M.; Mark, J.E.; Straube, E.; Vilgis, T.A. The Thermoelasticity of Rubberlike Materials and Related Constitutive Laws. J. Macromol. Sci. Part A Pure Appl. Chem. 2003, 40, 87–93. [Google Scholar] [CrossRef]
  46. Writing User Subroutines with ABAQUS. 2020. Available online: https://www.ymcn.org/2542987.html (accessed on 17 November 2020).
  47. Abaqus 6.14 Theory Guide. User Subroutine to Define a Material’s Mechanical Behavior. 2014. Available online: http://130.149.89.49:2080v6.14/books/sub/default.htm (accessed on 10 September 2020).
  48. Xiong, S.; Salant, R.F. A Numerical Model of a Rock Bit Bearing Seal. Tribol. Trans. 2000, 43, 542–548. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the cross-sectional structure (a) and contact deformation of FRR (b) in a roller bit. In the contact path (green dotted line in Figure 1b), the critical point on the rotation surface contact with the outer and inner fluid zone is A and B, point C is the middle point in the contact path and the critical point on the outside surface contact with the outer and inner fluid zone is D and E.
Figure 1. Schematic diagram of the cross-sectional structure (a) and contact deformation of FRR (b) in a roller bit. In the contact path (green dotted line in Figure 1b), the critical point on the rotation surface contact with the outer and inner fluid zone is A and B, point C is the middle point in the contact path and the critical point on the outside surface contact with the outer and inner fluid zone is D and E.
Materials 15 05529 g001
Figure 2. Three tensile test devices (a1c1), test section magnification (a2c2) and test samples’ extrinsic features (a3c3) after tests.
Figure 2. Three tensile test devices (a1c1), test section magnification (a2c2) and test samples’ extrinsic features (a3c3) after tests.
Materials 15 05529 g002
Figure 3. Schematic representation of uniaxial tensile (a), equivalent tensile (b), planar tensile (c) of stretched samples and deformation directions (d).
Figure 3. Schematic representation of uniaxial tensile (a), equivalent tensile (b), planar tensile (c) of stretched samples and deformation directions (d).
Materials 15 05529 g003
Figure 4. Multiplicative decomposition of deformation.
Figure 4. Multiplicative decomposition of deformation.
Materials 15 05529 g004
Figure 5. The FEM schematic diagram for rubber SAC.
Figure 5. The FEM schematic diagram for rubber SAC.
Materials 15 05529 g005
Figure 6. The results of three models fitted to HNBR stretching experiments at 120 °C: Yeoh model (a), Yeoh_revised model (b) and Ogden model (c).
Figure 6. The results of three models fitted to HNBR stretching experiments at 120 °C: Yeoh model (a), Yeoh_revised model (b) and Ogden model (c).
Materials 15 05529 g006
Figure 7. Residual analysis of three models for fitting results of HNBR tensile experiment at 120 °C: UT tensile (a), ET tensile (b) and PT tensile (c) results.
Figure 7. Residual analysis of three models for fitting results of HNBR tensile experiment at 120 °C: UT tensile (a), ET tensile (b) and PT tensile (c) results.
Materials 15 05529 g007
Figure 8. The SAC stress compared with FEM calculations at 120 °C.
Figure 8. The SAC stress compared with FEM calculations at 120 °C.
Materials 15 05529 g008
Figure 9. Contact stress cloud diagram of FRR in a roller bit: (a) pre-compression contact stress based on Yeoh model; (b) contact stress after fluid extrusion based on Yeoh model; (c) pre-compression contact stress based on Yeoh_revised model; (d) contact stress after fluid extrusion based on Yeoh_revised model.
Figure 9. Contact stress cloud diagram of FRR in a roller bit: (a) pre-compression contact stress based on Yeoh model; (b) contact stress after fluid extrusion based on Yeoh model; (c) pre-compression contact stress based on Yeoh_revised model; (d) contact stress after fluid extrusion based on Yeoh_revised model.
Materials 15 05529 g009
Figure 10. Contact stress on the contact path of FRR in a roller bit. In the contact path (green dotted line in Figure 1b), the critical point on the rotation surface contact with the outer and inner fluid zone is A and B, point C is the middle point in the contact path and the critical point on the outside surface contact with the outer and inner fluid zone is D and E.
Figure 10. Contact stress on the contact path of FRR in a roller bit. In the contact path (green dotted line in Figure 1b), the critical point on the rotation surface contact with the outer and inner fluid zone is A and B, point C is the middle point in the contact path and the critical point on the outside surface contact with the outer and inner fluid zone is D and E.
Materials 15 05529 g010
Figure 11. Mises stress cloud diagram of FRR in a roller bit: (a) pre-compression Mises stress based on Yeoh model; (b) Mises stress after fluid extrusion based on Yeoh model; (c) pre-compression Mises stress based on Yeoh_revised model; (d) Mises stress after fluid extrusion based on Yeoh_revised model. In the contact path (green dotted line in Figure 1b), the critical point on the rotation surface contact with the outer and inner fluid zone is A and B, point C is the middle point in the contact path and the critical point on the outside surface contact with the outer and inner fluid zone is D and E.
Figure 11. Mises stress cloud diagram of FRR in a roller bit: (a) pre-compression Mises stress based on Yeoh model; (b) Mises stress after fluid extrusion based on Yeoh model; (c) pre-compression Mises stress based on Yeoh_revised model; (d) Mises stress after fluid extrusion based on Yeoh_revised model. In the contact path (green dotted line in Figure 1b), the critical point on the rotation surface contact with the outer and inner fluid zone is A and B, point C is the middle point in the contact path and the critical point on the outside surface contact with the outer and inner fluid zone is D and E.
Materials 15 05529 g011
Figure 12. Mises stress on the contact path of FRR in roller bit. In the contact path (green dotted line in Figure 1b), the critical point on the rotation surface contact with the outer and inner fluid zone is A and B, point C is the middle point in the contact path and the critical point on the outside surface contact with the outer and inner fluid zone is D and E.
Figure 12. Mises stress on the contact path of FRR in roller bit. In the contact path (green dotted line in Figure 1b), the critical point on the rotation surface contact with the outer and inner fluid zone is A and B, point C is the middle point in the contact path and the critical point on the outside surface contact with the outer and inner fluid zone is D and E.
Materials 15 05529 g012
Figure 13. The on-site wear effect on FRR in a roller bit. The wear effect on the outside and rotation surfaces of FRR before and after use is shown in (a) and (c), respectively. The magnification of the outside and rotation surfaces’ wear effect is shown in (b) and (d), respectively. In the contact path (green dotted line in Figure 1b), the critical point on the rotation surface contact with the outer and inner fluid zone is A and B, the critical point on the outside surface contact with the outer and inner fluid zone is D and E.
Figure 13. The on-site wear effect on FRR in a roller bit. The wear effect on the outside and rotation surfaces of FRR before and after use is shown in (a) and (c), respectively. The magnification of the outside and rotation surfaces’ wear effect is shown in (b) and (d), respectively. In the contact path (green dotted line in Figure 1b), the critical point on the rotation surface contact with the outer and inner fluid zone is A and B, the critical point on the outside surface contact with the outer and inner fluid zone is D and E.
Materials 15 05529 g013
Figure 14. Schematic diagram of wear measurement of FRR in a roller bit.
Figure 14. Schematic diagram of wear measurement of FRR in a roller bit.
Materials 15 05529 g014
Table 1. SAC displacements corresponding to different tensile ratios λi.
Table 1. SAC displacements corresponding to different tensile ratios λi.
λi11.021.041.061.081.10
uzi06.16 mm7.26 mm8.30 mm9.28 mm10.21 mm
Table 2. Comparison of the fitted results of tensile stress based on three models.
Table 2. Comparison of the fitted results of tensile stress based on three models.
YeohYeoh_RevisedOgden
T120 °CT120 °CT120 °C
C101.3745C100.36u10.8152
C20−1.4273C20−1.3323α1−0.0003
C302.3639C303.0578u20.8152
/ / C010.9965α20.0005
/ / / / u30.8152
/ / / / α3−0.0002
SSdev0.7124SSdev0.2659SSdev0.6465
R20.9385R20.9771R20.9442
Table 3. Comparison between CAE-calculated and experimental values.
Table 3. Comparison between CAE-calculated and experimental values.
λTest_DataYeoh_CaeDeviationYeoh_Re_CaeDeviationOgden_CaeDeviation
T = 120 °C1.000.00000.000000.000000.00000
1.020.41370.3159−23.65%0.3213−22.34%0.2892−30.10%
1.040.73330.601−18.04%0.6315−13.88%0.5753−21.55%
1.061.02500.8451−17.55%0.9219−10.06%0.8559−16.50%
1.081.33461.048−21.48%1.195−10.46%1.1300−15.33%
1.101.63091.219−25.26%1.465−10.17%1.4000−14.16%
Notes: deviation = (xx_cae − test_date)/test_date × 100%.
Table 4. Comparison between CAE-calculated and fitted values.
Table 4. Comparison between CAE-calculated and fitted values.
T = 120 °C
λYeoh_CaeYeoh_FitDeviationλOgden_CaeOgden_FitDeviation
1.000.00000.000001.000.00000.00000
1.020.31590.3173−0.43%1.020.28920.2905−0.46%
1.040.60100.60080.04%1.040.57530.57520.01%
1.060.84510.84380.15%1.060.85590.85460.15%
1.081.04801.04660.13%1.081.13001.12900.09%
1.101.21901.21570.27%1.101.40001.39590.29%
λYeoh_re_caeYeoh_re_fitDeviation
1.000.00000.00000
1.020.32130.3227−0.43%
1.040.63150.63150.01%
1.060.92190.92070.13%
1.081.19501.19340.13%
1.101.46501.46080.29%
Notes: deviation = (xx_cae − xx_fit)/xx_cae × 100%.
Table 5. Flat sealing ring size and hardness measurements before and after use.
Table 5. Flat sealing ring size and hardness measurements before and after use.
ObjectBore Size/mmLength of Different Phase-Angle
Cross-Sections (l)
Height of Different Phase-Angle
Cross-Sections (h)
Hardness/HA
90°180°270°90°180°270°Outside SurfaceRotation Surface
Before use54.706.256.256.256.253.903.903.903.909090
After use55.226.035.965.956.053.923.923.883.8994101
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, W.; Wang, C.; Fan, P.; Kuang, Y.; Dong, Z. The Sealing Effect Improvement Prediction of Flat Rubber Ring in Roller Bit Based on Yeoh_Revised Model. Materials 2022, 15, 5529. https://doi.org/10.3390/ma15165529

AMA Style

Zhou W, Wang C, Fan P, Kuang Y, Dong Z. The Sealing Effect Improvement Prediction of Flat Rubber Ring in Roller Bit Based on Yeoh_Revised Model. Materials. 2022; 15(16):5529. https://doi.org/10.3390/ma15165529

Chicago/Turabian Style

Zhou, Wei, Chengwen Wang, Peng Fan, Yuchun Kuang, and Zongzheng Dong. 2022. "The Sealing Effect Improvement Prediction of Flat Rubber Ring in Roller Bit Based on Yeoh_Revised Model" Materials 15, no. 16: 5529. https://doi.org/10.3390/ma15165529

APA Style

Zhou, W., Wang, C., Fan, P., Kuang, Y., & Dong, Z. (2022). The Sealing Effect Improvement Prediction of Flat Rubber Ring in Roller Bit Based on Yeoh_Revised Model. Materials, 15(16), 5529. https://doi.org/10.3390/ma15165529

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