Data processing method for constitutive relation of the rubber core material in spherical blowout preventer

Abstract

A data processing method for constitutive relation of large deformation rubber core in spherical blowout preventer is put forward in this article. Based on the Yeoh constitutive model and the large-deformation theory of rubber material, this article studied the constitutive relation of core material in blowout preventer. The constitutive relation experiments, the uniaxial tensile experiment, the uniaxial compression experiment, the plane tensile experiments, and the volume compression experiment were conducted and the corresponding experiment data were fitted, respectively, on basis of which the uniaxial compression data were selected for further study. For comparison, the data were, respectively, processed with a novel method and the method in GB/7757, and accordingly, two different constitutive relations were obtained and introduced into ABAQUS to simulate the experiment process. Research results show that the relative error of maximum axial deformation between the simulation using GB/7757 and the actual experiment is 34.7%, while that between the novel method and the actual experiment is 20.4%, having improved the accuracy by 14.3%, proving that the data processing method in this article can effectively improve the accuracy of finite element simulation for the core material and will contribute to the research on the performance of spherical blowout preventer.

Keywords

Spherical blowout preventer rubber large deformation constitutive relation accuracy

Introduction

In the process of drilling, there are many uncertain factors in the well, and the potential blowout risk is more likely. Once an accident occurs, it will cause extremely serious consequences, so the well control equipment becomes an indispensable important equipment in the wellsite. The oil spill in the gulf of Mexico¹ is a serious blowout accident, causing heavy losses. Spherical blowout preventer (BOP) is an important part of well control equipment used to control the wellhead pressure to avoid blowout accident during drilling, well repairing, oil testing, and so on. As an indispensable device, spherical BOP is widely used in land drilling and deep-water drilling to ensure stable operation of drilling equipment and provide safety guarantee for the wellsite.²

The rubber core is a key part of the spherical BOP and its performance is significantly affected by the core material. The rubber core is a vulnerable part affected by many factors, such as cracking (Figure 1), severe spalling on the inner surface (Figure 2), and other phenomena that often occur in the core. All these factors will degrade the sealing performance and shorten the service life of the rubber core. In order to improve the reliability of the core, finite element simulation is utilized to analyze its performance. Since accurate constitutive relation is required in the simulation of the core material, while simulations with existing constitutive relations contain unacceptable huge error, a better constitutive relation is an urgent issue that needs to be solved. Nevertheless, as a typical rubber product, the rubber core has the characteristics of large deformation, viscoelasticity, and will be easily affected by environmental factors, such as temperature, pressure, humidity, which make it quite difficult to achieve an accurate constitutive relation for researchers.

Figure 1.

Cracks in the rubber core.

Figure 2.

Severe spalling on inner surface of the rubber core.

Nevertheless, there are a variety of constitutive relation models describing the rubber material, among which the representatives are the models of Mooney and Rivlin,^3–5Yeoh,^6,7Gent,⁸ which are relatively mature and reliable. Specifically, the Yeoh⁹ model for large deformation contact/impact issues has been proposed in recent years, which has properly improved the simulation accuracy of hyperelastic materials. ZQ Yu et al.¹⁰ proposed a large deformation constitutive model for thermoplastic olefin, describing the deformation of olefinic materials in a better way. I Caylak and R Mahnken¹¹ studied the stability of mixed tetrahedrons with large deformation and made some contributions to the simulation of superelastic materials. DC Moreira and LCS Nunes¹² reported that simple shear cannot be considered as pure shear combined with a rotation when large deformation is assumed. Y-D Kwon et al.¹³ analyzed the effect of Poisson’s ratio on the simulation results of a hyperelastic material. D Forcellini and JM Kelly¹⁴ studied the large deformation behavior of elastic bearings. A-S Lectez et al.¹⁵ introduced a method for the constitutive relationship of superelastic materials with multi-axis pull-torsion tests. TG Ritto and LCS Nunes¹⁶ studied the Bayesian model of superelastic materials. J Aboudi and KY Volokh¹⁷ predicted the failure of the composite with large deformation. WK Shi et al.¹⁸ discussed the effects of the bulk compressibility on rubber isolator’s compressive behaviors. B Crawford et al.¹⁹ proposed a new constitutive relation for silicones, and the constitutive model could predict the first normal stress and extensional behavior at high to medium shear rates. JA Rodríguez-Martínez et al.²⁰ studied the mechanical response of hyperelastic spherical membranes subjected to dynamic inflation. C Gokulnath et al.²¹ put forward a methodology for obtaining implicit constitutive representations involving the Cauchy stress and the Hencky strain. Many researchers have conducted extensive studies on the mechanical properties of hyperelastic materials. However, as a kind of high molecular material, hyperelastic materials often exhibit different mechanical properties if they have different molecular arrangements at the micro level. As a result, few reports about research progress on the constitutive relation of spherical rubber core material could be found.

Based on these research achievements of rubber materials, the author of this article established a novel constitutive model for the rubber core of spherical BOP. However, since each rubber material has its own constitutive relation, the existing constitutive models of rubber materials cannot meet the requirements of spherical BOP simulation. In this article, spherical core material is taken as the research object, and four constitutive relation experiments on the material were conducted. By fitting the experimental data, reasonable experimental data were found out, and then the selected data were both processed with conventional methods and the novel methods proposed by the author. By comparing the data processing results, the most suitable and reliable data processing method was found to describe the constitutive relation of the spherical core material.

Experimental study on constitutive relation of the rubber core in spherical BOP

Introduction to experiment apparatus

In the experiment, as shown in Figure 3, a 2-kN tensile tester is utilized, of which the main parameters are that the maximum experimental load is 2000 N, the maximum experimental speed is 500 mm/min, the experimental effective width is 150 mm, the effective stretching space distance is 700 mm, and the displacement resolution is 0.1 μm where the accuracy level is <0.5.

Figure 3.

Tensile tester used in the experiment.

Uniaxial tensile experiment

The experiment is conducted according to related rules in GB/T 528-2009 Rubber, Vulcanized or Thermoplastic Determination of Tensile Stress–Strain Properties, and structural sizes²² of the test sample are shown in Figure 4.

Figure 4.

Structural sizes of the tensile test sample.

In the experiment, six standard dumbbell-like samples were tested at room temperature; specifically, the samples were stretched at a constant speed of 500 ± 50 mm/min on the tensile tester. During the experiment, the length and force of the samples were continuously monitored until the samples were torn apart. The tested samples are illustrated in Figure 5, and a set of typical experimental data are given in Table 1.

Table 1.

A part of experimental data in the uniaxial tensile experiment.

Numbering	Amount of stretch	Strain ε	Strain ε (%)	Load P (N)	Stress σ (MPa)
1	0	0	0	0	0
2	7.474964	0.3737482	7.37482	8.79889	1.0998613
3	14.75586	0.7377930	73.77930	12.06704	1.5083800
4	22.38474	1.1192370	111.92370	16.34079	2.0425988
5	29.9032	1.4951600	149.51600	22.87711	2.8596388
6	37.52874	1.8764370	187.64370	32.68158	4.0851975
7	50.02605	2.5013025	250.13025	53.54751	6.6934388
8	65.7489	3.2874450	328.74450	81.45255	10.1815688
9	84.617	4.2308500	423.08500	114.38550	14.2981875
10	103.4115	5.1705750	517.05750	141.03360	17.6292000

Figure 5.

Tested rubber samples.

ABAQUS was used to analyze and evaluate the corresponding six sets of experimental data after the experiment. One of the ideal fitting curves is shown in Figure 6, wherein the obtained Yeoh model coefficients are that C₁₀ = 0.506136047 MPa, C₂₀ = 2.300575354 × 10⁻⁰² MPa, C₃₀ = −1.769116947 × 10⁻⁰⁴ MPa. In this figure, the blue curve is the fitting curve in Yeoh model, while the red one is the experimental data curve. It should be noted that the stress of the two curves changes alternately at different abscissas, and that when the rubber deformation rate exceeds 200%, the error of the Yeoh model increases gradually. The conclusion is that the measured stress value is obviously higher than the theoretical one with the same strain.

Figure 6.

Fitting curve of tensile experimental data.

Uniaxial compression experiment

According to GB/T 7757-2009 Rubber, Vulcanized or Thermoplastic Determination of Compression Stress–Strain Properties, the sample is a cylindrical one with a diameter of φ = 29.0 ± 0.5 mm and a height of h = 12.5 ± 0.5 mm. Five samples were put in the test and were compressed at a rate of 10 mm/min at room temperature until the strain reached 44%, and then the samples were rebounded at the same speed and the same temperature, thus repeatedly and continuously compressing and loosening the sample for four cycles.²³ In the four cycles, the first three can be regarded as mechanical adjustment. Results of four cycles were recorded and result of the last cycle was processed, wherein one of the typical experimental data is shown in Table 2.

Table 2.

A part experimental data in the uniaxial compression experiment.

Numbering	Original compressive load (N)	Original compression deformation (mm)	Compressive load after treatment (N)	Compression deformation after treatment (mm)
1	44.33786	0.5425627	0	0
2	179.2249	0.943983	134.88704	0.4014203
3	321.6056	1.320596	277.26774	0.7780333
4	451.4968	1.664199	407.15894	1.1216363
5	578.8901	1.991997	534.55224	1.4494343
6	735.0093	2.375412	690.67144	1.8328493
7	894.8754	2.759528	850.53754	2.2169653
8	1071.602	3.16615	1027.26414	2.6235873
9	1259.57	3.575372	1215.23214	3.0328093
10	1464.398	3.978894	1420.06014	3.4363313
11	1706.071	4.395918	1661.73314	3.8533553
12	1994.975	4.826447	1950.63714	4.2838843
13	2327.471	5.217365	2283.13314	4.6748023
14	2759.088	5.611483	2714.75014	5.0689203
15	3311.156	5.9997	3266.81814	5.4571373

ABAQUS was used to analyze and evaluate the corresponding four sets of experimental data after the experiment. One of the ideal fitting curves is shown in Figure 7, wherein the obtained Yeoh model coefficients are that C₁₀ = 1.15619829 MPa, C₂₀ = 1.18908883MPa,C₃₀ = 0.340784687 MPa. It could be observed in Figure 7 that when the rubber compression deformation rate exceeds 40%, the error of the Yeoh model will increase, but the difference is quite small, which indicates that the overall trend of the theoretical data is basically the same as the experimental data. Since the compression deformation is the main deformation mode of the rubber core in spherical BOP, the constituent coefficient of the core material can be obtained by uniaxial compression experimental data.

Figure 7.

Fitting curve of uniaxial compression experiment data.

Plane tensile experiment

The plane tensile model is that the dimension in the width direction is much larger than that in the direction of the tensile force. Similarly, the plane compression model is that the dimension in the length/width direction is much larger than that in the direction of the tensile force (thickness direction), as shown in Figure 8.

Figure 8.

Illustration of plane experiment.

In the experiment, each of the samples, of which the sizes are 85 mm × 0.85 mm × 44 mm, was fixed with a flat fixture and tested on the tensile tester at a constant speed of 40 mm/min. With the mentioned facilities and requirements, three sets of experiments are conducted, of which a set of typical data are listed in Table 3.

Table 3.

Partial experimental data in the plane tensile experiment.

Numbering	Tensile deformation (mm)	Tensile load (N)	Tensile stress (MPa)	Stretch strain
1	0	0	0	0
2	1.0435	49.799	0.65096732	0.2087
3	1.9822	86.25	1.12745098	0.39644
4	2.9208	107.37	1.403529412	0.58416
5	3.8586	121.7	1.590849673	0.77172
6	4.9652	134.77	1.761699346	0.99304
7	9.0827	172.23	2.251372549	1.81654
8	12.9053	198.37	2.593071895	2.58106
9	19.7177	233.57	3.053202614	3.94354
10	29.8243	282.84	3.697254902	5.96486

ABAQUS was used to analyze and evaluate the corresponding three sets of experimental data after the experiment. One of the typical fitting curves is shown in Figure 9, wherein the obtained Yeoh model coefficients are that C₁₀ = 0.485940296 MPa, C₂₀ = −4.306033840 × 10⁻⁰³ MPa, C₃₀ = 2.429941023 × 10⁻⁰⁵ MPa. Figure 9 shows a large deviation between the blue curve fitted with the Yeoh model and the red curve fitted with the experimental data; it can be noted that the fitting results are not so ideal but the general trend is the correct. When the rubber deformation rate exceeds 100%, the error is obviously increased, meaning that the constitutive relationship data of plane stretching model must be used carefully.

Figure 9.

Fitting curve of plane tensile experimental data.

Volume compression experiment

The volume compression experiment is mainly used to determine Poisson’s ratio of the material, but the constitutive model of the material cannot be determined from the volume compression experiment. Nevertheless, more accurate constitutive model of the core material can be obtained by combining the volume compression experiment and one of the above three experiments.

In this experiment, cylinder samples, of which the diameter φ = 29.0 mm and the height h = 12.5 mm, were chosen to be compressed on the tester at 10 mm/min at room temperature. Related experimental data are shown in Table 4.

Table 4.

A part of experimental data in the volume compression experiment.

Numbering	The amount of compression	Strain ε	Load P (N)	Stress σ (MPa)	Volume ratio
1	0	0	0	0	1
2	0.0615185	0.004921476	94.2960	0.1427604	1.004921476
3	0.1229369	0.009834952	304.7447	0.4613714	1.009834952
4	0.1748524	0.013988192	584.5104	0.8849255	1.013988192
5	0.2501750	0.020014000	1110.9440	1.6819251	1.020014000
6	0.3315995	0.026527960	1779.7590	2.6944845	1.026527960
7	0.4107232	0.032857856	2459.2140	3.7231524	1.032857856
8	0.4666400	0.037331200	2962.3500	4.4848803	1.037331200
9	0.5210563	0.041684504	3457.9570	5.2352096	1.041684504
10	0.5709713	0.045677704	3914.6680	5.9266520	1.045677704

The experimental data were fitted by ABAQUS software, and the Yeoh model coefficients were derived by introducing the experimental data from the uniaxial compression, of which the fitting results are shown in Figure 10. The obtained Yeoh model coefficients are that C₁₀ = 1.10942387 MPa, C₂₀ = 0.109411896 MPa, C₃₀ = 0.333224613 MPa, D₁ =−5.100869789 × 10⁻⁰², D₂ =−4.054079459 × 10⁻⁰⁵, D3 = 2.245549623 × 10⁻⁰⁷. The blue curve is fitted with the Yeoh model, and the red curve is fitted with the experimental data.

Figure 10.

Fitting curve of the data combined uniaxial and volume compression experiment.

According to the fitting degree of the curves in Figure 10, the Yeoh model can basically reflect the volume compression constitutive relation of rubber when the volume compressive deformation rate is smaller than 5%, while the error will increase sharply above this value. Since Poisson’s ratio of rubber material is quite large, the volume of rubber tends to be very small, which can be concluded from the experiment. Given this, a simple method, assuming that the core rubber is incompressible material and ignoring the effect of volume deformation, the constitutive data obtained are also closer to the experimental results.

Experimental data processing

By analyzing the data of the above four sets of experiments, it can be found that the deviation between the curve fitted with Yeoh model (uniaxial and plane tensile model) and the one fitted with experimental data is relatively large. If the uniaxial or the plane tensile data are used as the constitutive data in the finite element analysis, the result will be very distorted. Besides, during the simulation, ABAQUS software indicated that the two sets of experimental data are unstable hyperelastic materials. Thus, it can be concluded that neither uniaxial tensile experimental data nor plane tensile experimental data could be used for further analysis.

On the other hand, the experimental data of uniaxial compression can be well-fitted with Yeoh model, which can objectively reflect the change of core material in spherical BOP, so the uniaxial compression experimental data were selected as the constitutive data for finite element analysis. Nevertheless, it is found that the error of the results from finite element analysis is still large compared with the experimental data, meaning that the data have to be processed more scientifically in order to improve the simulation accuracy. The original data for uniaxial compression are shown in Table 5.

Table 5.

A part of the original data in the uniaxial compression experiment.

Numbering	Original compressive load (N)	Original compression deformation (mm)	Compressive load after treatment (N)	Compression deformation after treatment (mm)
1	44.33786	0.5425627	0	0
2	179.2249	0.943983	134.88704	0.4014203
3	321.6056	1.320596	277.26774	0.7780333
4	451.4968	1.664199	407.15894	1.1216363
5	578.8901	1.991997	534.55224	1.4494343
6	735.0093	2.375412	690.67144	1.8328493
7	894.8754	2.759528	850.53754	2.2169653
8	1071.602	3.16615	1027.26414	2.6235873
9	1259.57	3.575372	1215.23214	3.0328093
10	1464.398	3.978894	1420.06014	3.4363313
11	1706.071	4.395918	1661.73314	3.8533553
12	1994.975	4.826447	1950.63714	4.2838843
13	2327.471	5.217365	2283.13314	4.6748023
14	2759.088	5.611483	2714.75014	5.0689203
15	3311.156	5.9997	3266.81814	5.4571373

Routine experimental data processing

According to related rules in GB/T 528-2009 Rubber, Vulcanized or Thermoplastic Determination of Tensile Stress–Strain Properties, the compressive stress σ is given as follows

σ = \frac{P}{S_{0}}

(1)

And, the compressive strain ε is given as follows

ε = \frac{Δ L}{L_{0}}

(2)

where P is the force applied in the stress direction, S₀ is the original cross-sectional area perpendicular to the direction of the applied force in the sample, ΔL is the deformation of the sample along the stress direction, and L₀ is the initial size of the sample in the stress direction. The data in Table 5 were processed according to GB/T 7757-2009, based on which the stress–strain data, as shown in Table 6, were obtained.

Table 6.

Stress–strain data according to GB/T 7757-2009.

Numbering	Compressive load (N)	Compressive strain ε	Compressive stress σ (MPa)
1	0	0	0
2	134.88704	0.032113624	0.204213452
3	277.26774	0.062242664	0.419771999
4	407.15894	0.089730904	0.616421954
5	534.55224	0.115954744	0.809290191
6	690.67144	0.146627944	1.045648264
7	850.53754	0.177357224	1.28767899
8	1027.26414	0.209886984	1.555235822
9	1215.23214	0.242624744	1.839811673
10	1420.06014	0.274906504	2.149912873
11	1661.73314	0.308268424	2.515795894
12	1950.63714	0.342710744	2.953184714
13	2283.13314	0.373984184	3.456570035
14	2714.75014	0.405513624	4.11002049
15	3266.81814	0.436570984	4.945828824

The Yeoh model coefficients were obtained by fitting the data in Table 6 with ABAQUS software, wherein C₁₀=1.15619829 MPa, C₂₀=1.18908883 MPa, C₃₀=0.340784687 MPa, Poisson’s ratio μ = 0.5. Accordingly, the finite element model for uniaxial compression was established in ABAQUS, as illustrated in Figure 11.

Figure 11.

Finite element model for uniaxial compression.

According to the actual situation of the experiment, the load, which is 3268.82 N and evenly distributed, from the upper iron is applied vertically downward on the sample, while the lower iron is fixed. The rubber sample is disposed between the upper and lower iron without friction. Thus, the stress distribution of uniaxial compression is obtained, as shown in Figure 12. Accordingly, the axial compression displacement of the rubber sample is achieved and shown in Figure 13.

Figure 12.

Simulation analysis results with the stress-strain data according to GB /T 7757-2009.

Figure 13.

Displacement cloud in the rubber sample with the stress-strain data according to GB /T 7757-2009.

According to the simulation results, it can be seen that when σ_max the rubber sample reaches 3.54 MPa, the simulated axial load is equivalent to the actual 3268.82 N. The axial compressed displacement obtained from the finite element software is 3.566 mm, while the measured axial displacement is 5.457 mm. The error between the simulated and the measured values reaches 34.7%, which is quite large and unacceptable. Further analysis revealed the cause of the error which is that the National Standard GB/T 7757-2009 data processing theories are based on the small deformation theory, while actually the rubber material is of large deformation category. In order to improve the accuracy of the simulation results, axial compression data are reasonably processed according to the large deformation theory.

Data processing method based on large deformation and rubber incompressibility

The experimental data of uniaxial compression were processed with the large deformation theory of incompressible rubber. Assuming that the volume of the sample remains constant before and after compression, the following equations can be obtained

V_{0} = V_{1}

(3)

V_{0} = \frac{π d_{0}^{2}}{4} L_{0}

(4)

V_{1} = \frac{π d_{1}^{2}}{4} (L_{0} - Δ L)

(5)

where V₀ is the volume of the cylindrical sample before the compression load is applied, V₁ is the volume after, d₀ is the original cross-sectional diameter of the sample before, d₁ is the cross-sectional diameter after, L₀ is the height before, and ΔL is the height increment after the compression load applied.

According to the assumption of the volume being constant before and after compression, equation (6) can be derived from equations (4) and (5)

d_{1} = \sqrt{\frac{d_{0}^{2} L_{0}}{L_{0} - L_{1}}}

(6)

Therefore, the actual compressive stress σ of the sample is as follows

σ = \frac{P}{S_{1}} = \frac{4 P}{π d_{1}^{2}}

(7)

The actual compressive strain ε of the sample is as follows

ε = \frac{Δ L}{L_{0}}

(8)

Based on the above formula, the stress and strain data, as shown in Table 7, are obtained.

Table 7.

The stress–strain partial data based on the theory of volume incompressibility.

Numbering	Compressive strain ε	Compressive stress σ (MPa)	Radius after compression (mm)
1	0	0	14.5
2	0.032113624	0.197655418	14.73858579
3	0.062242664	0.393644271	14.9734807
4	0.089730904	0.561109855	15.19788322
5	0.115954744	0.715449154	15.42164683
6	0.146627944	0.892327009	15.69635431
7	0.177357224	1.059299819	15.98683005
8	0.209886984	1.228812066	16.31260778
9	0.242624744	1.393427837	16.66143665
10	0.274906504	1.558887841	17.02828842
11	0.308268424	1.740255458	17.43408672
12	0.342710744	1.941096583	17.88503258
13	0.373984184	2.163867511	18.32632354
14	0.405513624	2.443351187	18.80602673
15	0.436570984	2.786623468	19.31738777

Through processing the data in Table 7 with ABAQUS software, the Yeoh model coefficients can be obtained, which are C₁₀=1.10942387 MPa, C₂₀=0.10 9411896 MPa, C₃₀ = 0.333224613 MPa, and Poisson’s ratio μ = 0.5. According to the actual situation of the experiment, the load, which is 3268.82 N and evenly distributed, from the upper iron is applied vertically downward on the sample, while the lower iron is fixed. The rubber sample is disposed between the upper and lower iron without friction. Thus, obtaining the stress distribution of uniaxial compression in Figure 14, accordingly, the axial compression displacement of the rubber sample is achieved and shown in Figure 15.

Figure 14.

Simulation analysis results with the stress-strain data based on the theory of volume incompressibility.

Figure 15.

Displacement cloud in the rubber sample with the stress-strain data based on the theory of volume incompressibility.

According to the simulation results, it can be seen that when σ_max the rubber sample reaches 3.54 MPa, the simulated axial load is equivalent to the actual 3268.82 N. The axial compressed displacement obtained from the finite element software is 4.343 mm, while the measured axial displacement is 5.457 mm. The error between the simulated and the measured values is 20.4%, which is decreased by14.3% compared with the result in “Routine experimental data processing”. In order to improve the accuracy of finite element analysis, the experimental data of axial compression are processed according to the large deformation theory of incompressible rubber, which is obviously better than the data processing method in the GB/T 7757-2009.

Conclusion

This article is intended to study the data processing method for constitutive relation of rubber material in spherical BOP. First, original constitutive relation data of the rubber material were obtained by uniaxial tensile experiment, uniaxial compression experiment, plane tensile experiment, and volume compression experiment. Furthermore, proper data, that is, the uniaxial compression experiment data, were screened out by utilizing ABAQUS. And then, these original data were processed with two different processing methods so that the optimal constitutive relation data are obtained.

When the data are processed with GB/T 7757 methods, the rubber constitutive relation is of low accuracy. In view of this problem, this article proposes a data processing method based on the large deformation theory of incompressible rubber. Through selection, the uniaxial compression experiment data are selected to be furtherly processed with two different processing methods, so that the optimal constitutive relation data are obtained and imported into ABAQUS software for simulation. Comparisons of the simulation results show that the calculation error using the non-compressible method is 20.4%, which is 14.3% lower than 34.7% obtained using the method of GB/T 7757. So it can be concluded that the non-compressible method significantly improves the data accuracy of the constitutive relation. This study has effectively improved the constitutive relationship accuracy of the rubber core material, laid the foundation for the simulation calculation of the spherical BOP, and will also effectively improve the performance and service life of the rubber core.

Footnotes

Handling Editor: Teen-Hang Meen

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) received no financial support for the research, authorship, and/or publication of this article.

ORCID iD

Daping Xu

References

Ramírez-León

. The oil spill in the Gulf of Mexico. In: Klapp

Cros

Fuentes

et al . (eds) Experimental and theoretical advances in fluid dynamics. Berlin: Springer, 2011, pp.103–134.

Yeo

Macgregor

Pinkstone

et al . Offshore deepwater managed pressure drilling and riser gas handling equipment design. In: SPE/IADC Managed Pressure Drilling and Underbalanced Operations Conference & Exhibition, 2015.

Mooney

. A theory of large elastic deformation. J Appl Phys 1940; 11: 582–592.

Rivlin

. Large elastic deformations of isotropic materials IV: further developments of the general theory. Philos T Roy Soc A 1948; 241: 379–397.

Rivlin

Saunders

. Large elastic deformations of isotropic materials VII: experiments on the deformation of rubber. Philos T Roy Soc A 1950; 243: 251–288.

Yeoh

. Characterization of elastic properties of carbon-black-filled rubber vulcanizates. Rubber Chem Technol 1990; 63: 792–805.

Yeoh

. Some forms of the strain energy function for rubber. Rubber Chem Technol 1993; 66: 754–771.

Gent

. A new constitutive relation for rubber. Rubber Chem Technol 1996; 69: 59–61.

Renaud

Cros

J-M

Feng

Z-Q

et al . The Yeoh model applied to the modeling of large deformation contact/impact problems. Int J Impact Eng 2009; 36: 659–666.

10.

et al . Constitutive modeling for large deformation behavior of thermoplastic olefin. Mater Design 2010; 31: 1881–1886.

11.

Caylak

Mahnken

. Stabilization of mixed tetrahedral elements at large deformations. Int J Numer Meth Eng 2012; 90: 218–242.

12.

Moreira

Nunes

LCS

. Comparison of simple and pure shear for an incompressible isotropic hyperelastic material under large deformation. Polym Test 2013; 32: 240–248.

13.

Kwon

Y-D

Kwon

S-B

et al . A finite element procedure with Poisson iteration method adopting pattern approach technique for near-incompressible rubber problems. Adv Mech Eng 2014; 6: 1–12.

14.

Forcellini

Kelly

. Analysis of the large deformation stability of elastomeric bearings. J Eng Mech 2014; 140: 1–15.

15.

Lectez

A-S

Verron

Huneau

. How to identify a hyperelastic constitutive equation for rubber-like materials with multiaxial tension–torsion experiments. Int J Nonlin Mech 2014; 65: 260–270.

16.

Ritto

Nunes

LCS

. Bayesian model selection of hyperelastic models for simple and pure shear at large deformations. Comput Struct 2015; 156: 101–109.

17.

Aboudi

Volokh

. Needleman symposium: failure prediction of unidirectional composites undergoing large deformations. J Appl Mech2015; 82: 1–15.

18.

Shi

Liu

Chen

. Effects of the bulk compressibility on rubber isolator’s compressive behaviors. Adv Mech Eng 2017; 9: 1–12.

19.

Crawford

Doherty

Spedding

et al . Constitutive relations for the extrusion of siloxane gum and silicone rubber. Asia-Pac J Chem Eng 2010; 5: 895–901.

20.

Rodríguez-Martínez

Fernández-Sáez

Zaera

. The role of constitutive relation in the stability of hyper-elastic spherical membranes subjected to dynamic inflation. Int J Eng Sci 2015; 93: 31–45.

21.

Gokulnath

Saravanan

Rajagopal

. Representations for implicit constitutive relations describing non-dissipative response of isotropic materials. Z Angew Math Phys 2017; 68: 129.

22.

GB/T 528-2009. Rubber, vulcanized or thermoplastic: determination of tensile stress-strain properties.

23.

GB/T 7757-2009. Rubber, vulcanized or thermoplastic: determination of compression stress-strain properties.