Mathematical Problems in Engineering
Volume 2019 (2019), Article ID 2832059, 10 pages
https://doi.org/10.1155/2019/2832059

Research Article

Modeling and Verification of a New Hyperelastic Model for Rubber-Like Materials

Zihan Zhao,1,2 Xihui Mu,¹ and Fengpo Du¹

¹Shijiazhuang Campus, Army Engineering University, Shijiazhuang 050003, China
²Special Services Institute, Shijiazhuang 050003, China

Correspondence should be addressed to Xihui Mu; 774589304@qq.com

Received 28 March 2019; Revised 17 April 2019; Accepted 18 April 2019; Published 2 May 2019

Academic Editor: Jose Merodio

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

Abstract.

The essential evaluation criterion for the hyperelastic model is its ability to describe the mechanical behavior of rubber-like materials under different deformation modes over a large deformation range accurately. Based on the Seth strain tensor invariant, a new hyperelastic model for isotropic and incompressible rubber-like materials is proposed. In order to investigate the prediction ability of the new model, the parameters of the new model, the Yeoh model, and the Carroll model are identified by test data of 8% vulcanized rubber and two different types of carbon black filled rubber, respectively. To this end, the data of uniaxial tension and equibiaxial tension are used simultaneously. Then, the same set of model parameters is used for prediction of pure shear (plane tension) deformation. The results show that the new model not only can predict the test data of pure shear (or plane tension) accurately, but also can be reliable to describe the response of various rubber materials over a large deformation range. Finally, the finite element simulation and experiment on static stiffness of rubber bushing are carried out based on the new model. By comparison of the experimental data with the simulation data, the new model can accurately reflect the mechanical behavior of rubber bushing. The new model can be used for performance analysis of rubber products and has better application value.

1. Introduction

Rubber materials are used in engineering fields widely, such as tires, rubber tracks, vehicle seals, and vibration-isolation devices [1–3]. Without considering the time effect, rubber-like materials are generally considered to be isotropic and incompressible hyperelastic materials. It is of great significance to study the performance of rubber products to establish a reliable hyperelastic model which can accurately reflect mechanical behavior of rubber-like materials. The essential evaluation criterion for the hyperelastic model is its ability to reproduce the mechanical behavior of rubber-like materials in different deformation states over a large range accurately. The researches of Hossain et al. [4, 5] on representative hyperelastic models show that their model cannot reliably predict mechanical behavior of rubber-like materials in all deformation modes by using only a certain type of test data. In addition, these models cannot also perfectly reproduce the stress-strain relationship over a large range

In recent years, many scholars have analyzed the disadvantages of existing hyperelastic models and developed many new hyperelastic models based on the work of predecessors. Hossain [6] proposed a three-dimensional electro-elastic constitutive framework that can model the stiffness gaining during the curing process undergoing finite deformations. By improving Mooney-Rivlin model, Yaya and Bechir [7] proposed a new compressible hyperelastic model with four parameters, and the response quality of this model is equivalent to that of the Ogden six-parameter model. In contrast to the traditional model, Crespo et al. [8] proposed a new hyperelastic model by What-You-Prescribe-Is-What-You-Get formulations. The form of Crespo model is not assumed beforehand and there are no material parameters. Vu et al. [9] proposed a micromechanical model for rubber elasticity on the basis of analytical networking-averaging of the tube model and by applying a closed-form of the Rayleigh exact non-Gaussian chains. Based on Lie group methods, Zhao [10] postulated a partial differential equation for isotropic hyperelastic constitutive models. Nkenfack et al. [11] proposed a new approach named Hybrid Integral Approach to model the incompressible isotropic hyperelastic behavior of rubber-like materials. Based on the Arruda and Boyce model, this model includes an original part made of an integral density and an interleaving constraint part represented by a logarithmic function. However, this model contains six material parameters and the application is limited. Aiming at the deficiency of the Arruda-Boyce model in predicting the equibiaxial data, Hossain et al. [12] compared five modified versions of the Arruda-Boyce model and two modified versions of full-network model with the Arruda-Boyce model. Bahreman and Darijani [13] proposed a new polynomial hyperelastic model which is a function of the principal invariants of the left Cauchy-Green strain tensor. Based on the neo-Hookean model, Bechir et al. [14] proposed a new strain energy function, which can predict the test data in multiple deformation modes by only using uniaxial tensile test data. However, this model only can be reliable to describe the response of rubbers in the small deformation range. Horgan et al. [15] compared the Fung-Demiray model and the Vito model to the stress response quality in various deformation modes. They believe that the hyperelastic model with the second invariant of the right Cauchy-Green deformation tensor can more accurately reflect the mechanical behavior and some physical effects of rubber-like materials. Under the incompressible assumption, Lin [16] derived the mathematical relationship between the right Cauchy-Green deformation tensor principal invariants and used the interpolation algorithm to obtain a new hyperelastic model that can predict stress-strain relationships by only using data of uniaxial tension and equibiaxial tension. Based on the logarithmic strain tensor invariants, Xiao et al. [17] proposed an explicit method for constructing multiaxial elastic potential only through test data in uniaxial deformation mode. Based on the research of Xiao, Yu et al. [18] proposed a new hyperelastic model, which is suitable for general compressible deformation of rubber-like materials, and the parameters have direct physical meaning.

Compared with the classical model, the hyperelastic model proposed in recent years can generally better reflect the mechanical behavior of rubber-like materials in various deformation modes; while the constitutive equations are more complicated, parameters are more and derivation methods are more cumbersome. These defects limit the application of the models. In view of the above problems, a new hyperelastic model is proposed based on the Seth strain tensor invariant. The basic test data of 8% vulcanized rubber and NR55 carbon black filled rubber are used for parameters identification. It seems that the data of uniaxial tension and equibiaxial tension are sufficient to obtain prefect model parameters. Finally, the finite element model of rubber bushing is established based on the new model by ABAQUS software. By comparing experimental data with simulation data, the static stiffness characteristics of bushing are analyzed to examine the applicability of the new model furthermore.

2. The General Forms of the Hyperelastic Constitutive Model

Based on the hyperelastic theory [19], the stress-strain relationship of rubber-like materials can be derived from the strain energy function . The relationship between the first Piola-Kirchhoff stress tensor and the deformation gradient tensor is as follows.

Under the assumption of isotropic and isothermal materials, (1) can be written as

where and are the principal invariants of the right Cauchy-Green deformation tensor . In general, rubber-like materials are considered to be incompressible; then (2) can be written as

where is the Lagrangian multiplier which can be eliminated by the boundary conditions in different deformation modes. After eliminating , the corresponding stress-strain relationship can be determined by (3).

Nowadays, three kinds of tests, namely, uniaxial tension, equibiaxial tension, and pure shear (or plane tension), are usually used to determine the mechanical behavior of rubber-liked materials. Under the assumption that the base vector of the Cartesian coordinate system is , the deformation gradient is as follows for uniaxial tension, equibiaxial tension, and pure shear (or plane tension), respectively:

where is the stretch ratio. After eliminating by the boundary conditions under different deformation modes, the relationship between nominal stress and stretch ratio can be written as follows.

3. Hyperelastic Constitutive Model Based on Seth Strain Tensor

The stress-strain relationship of rubber-like materials is significantly nonlinear over the large deformation range. The hyperelastic models which do not contain correlation term of the second invariant I₂ cannot accurately reproduce the test data in multiaxial deformation state [15, 20]. Based on the invariants of Seth strain tensor, a highly nonlinear hyperelastic model is established in this paper, which can accurately reflect the mechanical behavior of rubber-like materials in various deformation modes over the large range.

When , the expression of the Seth strain tensor is defined as follows.

According to the polar decomposition theorem, the deformation gradient tensor can be written as

where and are right and left stretch tensors, respectively, and is the rotation tensor. The eigenvalues of and are the principal stretches . According to spectral decomposition, can be written as

where are orthogonal eigenvectors of . Therefore, (6) can be expressed as follows.

The first invariant of Seth strain tensor is as follows.

Following the approach of Bechir et al. [14], the first invariant of the Seth strain tensor is chosen to generalize the strain energy function. Expanding in powers of , the following function can be obtained:

where takes an arbitrary integer that is not 0. Obviously, (11) satisfies the Valanis-Landel hypothesis [21]. The relationships between the principal invariants of the right Cauchy-Green deformation tensor and the principal stretches are as follows.

When and , (11) is the neo-Hookean model.

When and , (11) is the Mooney-Rivlin model.

When takes positive and negative values simultaneously, the strain energy function can more accurately reflect the influence of each invariant on the mechanical behavior of rubber-like materials [22]. Higher order terms help to improve the goodness of fit under large deformation conditions [23]. Considering the above factors, the strain energy function is defined as follows.

Based on (12), the following expressions can be obtained.

Inserting (16) into (15), the strain energy function can be written as a function of the invariants of the right Cauchy-Green strain tensor.

4. Model Parameters Identification

The parameters of the hyperelastic model are usually determined by basic tests such as uniaxial tension, equibiaxial tension, and pure shear (or plane tension). When the parameters are identified by data of uniaxial tension and equibiaxial tension, the model can more accurately reflect the stress-strain relationship for all deformations [24]. Considering that the new strain energy is the function of invariants, in order to examine its performance, Yeoh model [25] as (18) and Carroll model [24] as (19) are selected as the comparison object.

Firstly, the test data of 8% vulcanized rubber reported by Treloar [26] are used for parameters identification. It should be emphasized that only the data of uniaxial tension and the equibiaxial tension are used simultaneously in the process of parameters identification. Then, the same values of model parameters are used to predict the stress-strain relationship in pure shear state. Curve fitting is performed by Universal Global Optimization in the 1stOpt software. The parameters values of each model are shown in Table 1. The fitting and prediction results are shown in Figure 1.

Table 1: The parameters values of each model for different rubber materials.

Model type
8% vulcanized rubber
the first carbon black filled rubber
NR55 carbon black filled rubber

New model
;
;
;

;
;
;

;
;
;

.
.
.

Yeoh model
;
;
;

;
;
;

.
.
.

Carroll model
;
;
;

;
;
;

.
.
.

(a)

(b)

(c)

Figure 1: Comparison of test data of 8% vulcanized rubber with results of each model: (a) uniaxial tension data and fitting results; (b) equibiaxial tension data and fitting results; (c) pure shear data and predicting results.

In order to evaluate the goodness of fit of each model, the coefficients of determination R² are calculated [27]:

where , , is the test values, is the model fit values, is the average values of , and N is the number of test data. In order to save space, the goodness of fit of each model is shown in Figures 1–3.

(a)

(b)

(c)

Figure 2: Comparison of test data of the first carbon black filled rubber with results of each model: (a) uniaxial tension data and fitting results; (b) equibiaxial tension data and fitting results; (c) pure shear data and predicting results.

(a)

(b)

(c)

Figure 3: Comparison of test data of NR55 carbon black filled rubber with results of each model: (a) uniaxial tension data and fitting results; (b) equibiaxial tension data and fitting results; (c) pure shear data and predicting results.

It is obvious that both of the new model and Carroll model are better than the Yeoh model for reproducing and predicting the test data of 8% vulcanized rubber. In particular, the Yeoh model exhibits a “softer” property when fitting the equibiaxial tension data. This is maybe due to the lack of the second invariant I₂ in the Yeoh model [28]. According to analysis of Carroll [24], Treloar’s data for uniaxial tension and equibiaxial tension are suitable for developing strain energy function, because the stretch values are large, λ=7.6 in uniaxial tension and λ=4.45 in equibiaxial tension. However, not all hyperelastic materials have this characteristic, for example, some carbon black filled rubber materials.

Based on the above analysis, test data of two different carbon black filled rubber materials are selected to further investigate the application scope of the new model. The formula of the first carbon black filled rubber is shown in Table 2 and the test data are shown in Figure 2 [29]. The other is NR55 carbon black filled rubber and the test data are shown in Figure 3 [30]. Likewise, the data of uniaxial tension and equibiaxial tension are used simultaneously in the process of parameters identification. Then, the same values of model parameters are used to predict the stress-strain relationship in pure shear state. The parameters values of each model are shown in Table 1. The fitting and prediction results are shown in Figures 2 and 3.

Table 2: Material formula of the first carbon black filled rubber.

Composition
Content/phr
Composition
Content/phr

Natural rubber (RSS3)
100
Microcrystalline wax
2

Zinc oxide
5
Solid benzofuran
2

Stearic acid
2
Sulfur(200 mesh)
2.5

Antioxidant
5
Accelerator
1.4

Carbon black (N550)
20

It is obvious that the stretch values of two carbon black filled rubber materials are less than that of 8% vulcanized rubber. Although the Yeoh model can keep the curves “S” shaped in all deformation states, its accuracy is worse than that of the Carroll model and the new model. Comparing the Carroll model and the new model, one finds that the curves shapes of both models are similar. By analyzing the two models, it can be found that both models contain high-order terms of invariant I₁, and the orders of invariant I₂ are different. The order of the I₂ terms in the Carroll model is 1/2, and the highest order of the I₂ terms in the new model is 2. Although results of the Carroll model and the new model are all satisfactory, the new model can be reliable to reproduce and predict the behavior of two carbon black filled rubber materials with a better approximation. The reason may be that the Carroll model is developed based on Treloar’s data, so it is more suitable for rubber materials with large stretch values. Therefore, the new model can accurately reflect the mechanical behavior of rubber-like materials for large and small strain and has a wider range of applications.

5. Model Verification by Static Stiffness Test of Rubber Bushing

In order to investigate the prediction power for mechanical behavior of rubber-like materials, the finite element model of rubber bushing is established based on the new model by the subroutine UHYPER of ABAQUS software. The static stiffness of rubber bushing is studied along with comparison between simulation and experimental results.

The object is the McPherson suspension comfort bushing, in which the rubber material is NR55 carbon black filled rubber. The bushing consists of an inner steel ring, an outer steel ring, and rubber which is bonded to steel rings by vulcanization. The finite element model of the rubber bushing is shown in Figure 4.

Figure 4: The finite element model of the rubber bushing.

The hyperelastic material is defined by the subroutine UHYPER of ABAQUS based on the new model whose parameters are the same as in Table 1 ( = 2.4658×10⁻² MPa, = 7.9451×10⁻¹ MPa, = -5.6729×10⁻² MPa, = 1.6134×10⁻³ MPa). The inner and outer steel rings are set as discrete rigid, using R3D4 units. The rubber is set as deformable body, using C3D8H hybrid units. In order to simulate the actual working condition, the contact surfaces between steel rings and rubber are set to tie constraint. The model is applied with an axial load of 1000N, a radial load of 3000N, and a rotational load of 16° around Y-axis, respectively. The simulation results are shown in Figure 5.

(a)

(b)

(c)

Figure 5: (a) The simulation result of radial stiffness; (b) the simulation result of axial stiffness; (c) the simulation result of torsional stiffness.

The static stiffness of rubber bushing is tested by MST 831 stiffness testing instrument with a loading rate of 0.01mm/s and 0.1°/s, respectively, as shown in Figure 6. In order to avoid the influence of the Mullins effect, the loading-unloading process is repeated four times, and then the fifth loading curve is taken as the final results.

Figure 6: The radial stiffness and torsional stiffness test.

It can be seen from Figure 7 that the radial stiffness curve and the axial stiffness curve are in good agreement with the experimental data. The maximum error is 7.10% and 4.75%, respectively. The torsional stiffness curve is slightly larger than the experimental data, and the maximum error is 19.40%. The reason for the errors is mainly that the vulcanization viscosity effect between the steel rings and rubber is neglected in order to reduce simulation difficulty. Conclusively, the new model can be applied to the static stiffness analysis of rubber bushing and is reliable to reflect the mechanical behavior of bushing.

(a)

(b)

(c)

Figure 7: (a) Comparison of experimental data with simulation data of radial stiffness; (b) comparison of experimental data with simulation data of axial stiffness; (c) comparison of experimental data with simulation data of torsional stiffness.

6. Conclusion

In this paper, a hyperelastic model for an isotropic and incompressible rubber-like materials is proposed based on the Seth strain tensor invariants. The new model is a function of the invariants of the right Cauchy-Green deformation tensor, satisfying the Valanis-Landel hypothesis. In order to investigate the prediction ability and applicability, the Yeoh model and the Carroll model which are based on invariants are selected as comparison objects, and the test data of 8% vulcanized rubber and two different types of carbon black filled rubber are used for parameters identification. The new model has a better prediction ability for different rubber materials than the Yeoh model and the Carroll model. The reason may be that the new model contains high-order terms of invariants I₁ and I₂, while the Yeoh model does not contain the terms of I₂ and the Carroll model does not contain the high-order terms of I₂. In addition, one advantage of the new model is that the satisfactory parameters can be identified by using uniaxial tension and equibiaxial tension data simultaneously.

Furthermore, the finite element simulation and experiment for the static stiffness of rubber bushing are studied based on the new model. The maximum error of radial stiffness, axial stiffness, and torsional stiffness is 7.10%, 4.75%, and 19.40%, respectively. The efficiency of the new model in finite element simulations is verified. The new model can be applied to the static stiffness analysis of rubber bushing and is reliable to reflect the mechanical behavior of bushing.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The present work was funded by the National Natural Science Foundation of China (61471385).

References

R. Rugsaj and C. Suvanjumrat, “Finite element analysis of hyperelastic material model for non-pneumatic tire,” Key Engineering Materials, vol. 775, pp. 554–559, 2018.

C. Cheng, S. Li, Y. Wang, and X. Jiang, “Force and displacement transmissibility of a quasi-zero stiffness vibration isolator with geometric nonlinear damping,” Nonlinear Dynamics, vol. 87, no. 4, pp. 2267–2279, 2017.

L. P. Li, “Experiment analysis about mechanical properties of rubber bushing for suspension telescopic shock absorber,” Applied Mechanics and Materials, vol. 670-671, pp. 1008–1011, 2014.

P. Steinmann, M. Hossain, and G. Possart, “Hyperelastic models for rubber-like materials: Consistent tangent operators and suitability for Treloar's data,” Archive of Applied Mechanics, vol. 82, no. 9, pp. 1183–1217, 2012.

M. Hossain and P. Steinmann, “More hyperelastic models for rubber-like materials: consistent tangent operators and comparative study,” Journal of the Mechanical Behavior of Materials, vol. 22, no. 1-2, pp. 27–50, 2013.

M. Hossain, “Modelling the curing process in particle-filled electro-active polymers with a dispersion anisotropy,” Continuum Mechanics and Thermodynamics, 2019.

K. Yaya and H. Bechir, “A new hyper-elastic model for predicting multi-axial behaviour of rubber-like materials: formulation and computational aspects,” Mechanics of Time-Dependent Materials, vol. 22, no. 2, pp. 167–186, 2018.

J. Crespo, M. Latorre, and F. Montáns, “WYPIWYG hyperelasticity for isotropic, compressible materials,” Computational Mechanics, vol. 59, no. 1, pp. 73–92, 2017.

V. N. Khiêm and M. Itskov, “Analytical network-averaging of the tube model: rubber elasticity,” Journal of the Mechanics and Physics of Solids, vol. 95, pp. 254–269, 2016.

F. Zhao, “Continuum constitutive modeling for isotropic hyperelastic materials,” Advances in Pure Mathematics, vol. 06, no. 09, pp. 571–582, 2016.

A. Nguessong Nkenfack, T. Beda, Z.-Q. Feng, and F. Peyraut, “HIA: a hybrid integral approach to model incompressible isotropic hyperelastic materials - part 1: theory,” International Journal of Non-Linear Mechanics, vol. 84, pp. 1–11, 2016.

M. Hossain, A. Amin, and M. N. Kabir, “Eight-chain and full-network models and their modified versions for rubber hyperelasticity: a comparative study,” Journal of the Mechanical Behavior of Materials, vol. 24, no. 1-2, pp. 11–24, 2015.

M. Bahreman and H. Darijani, “New polynomial strain energy function; application to rubbery circular cylinders under finite extension and torsion,” Journal of Applied Polymer Science, vol. 132, no. 13, 2015.

H. Bechir, L. Chevalier, M. Chaouche, and K. Boufala, “Hyperelastic constitutive model for rubber-like materials based on the first Seth strain measures invariant,” European Journal of Mechanics - A/Solids, vol. 25, no. 1, pp. 110–124, 2006.

C. O. Horgan and M. G. Smayda, “The importance of the second strain invariant in the constitutive modeling of elastomers and soft biomaterials,” Mechanics of Materials, vol. 51, pp. 43–52, 2012.

B. Lin, “A new model for hyperelasticity,” Acta Mechanica, vol. 208, no. 1-2, pp. 39–53, 2009.

H. Xiao, “Elastic potentials with best approximation to rubberlike elasticity,” Acta Mechanica, vol. 226, no. 2, pp. 331–350, 2015.

L. Yuan, Z.-X. Gu, Z.-N. Yin, and H. Xiao, “New compressible hyper-elastic models for rubberlike materials,” Acta Mechanica, vol. 226, no. 12, pp. 4059–4072, 2015.

M. M. Attard and G. W. Hunt, “Hyperelastic constitutive modeling under finite strain,” International Journal of Solids and Structures, vol. 41, no. 18-19, pp. 5327–5350, 2004.

K. Yaya, H. Bechir, and F. Bremand, “Implementation of new strain-energy density function for a grade of carbon black-filled natural rubber in finite element code,” in Proceedings of the 6th International Conference on Advances in Mechanical Engineering and Mechanics (ICAMEM2015), Hammamet, Tunisia, 2015.

K. C. Valanis and R. F. Landel, “The strain‐energy function of a hyperelastic material in terms of the extension ratios,” Journal of Applied Physics, vol. 38, no. 7, pp. 2997–3002, 1967.

K. Farahani and H. Bahai, “Hyper-elastic constitutive equations of conjugate stresses and strain tensors for the Seth-Hill strain measures,” International Journal of Engineering Science, vol. 42, no. 1, pp. 29–41, 2004.

T. Beda, “An approach for hyperelastic model-building and parameters estimation a review of constitutive models,” European Polymer Journal, vol. 50, no. 1, pp. 97–108, 2014.

M. M. Carroll, “A strain energy function for vulcanized rubbers,” Journal of Elasticity: The Physical and Mathematical Science of Solids, vol. 103, no. 2, pp. 173–187, 2011.

O. H. Yeoh, “Some forms of the strain energy function for rubber,” Rubber Chemistry and Technology, vol. 66, no. 5, pp. 754–771, 1993.

L. R. G. Treloar, “Stress-strain data for vulcanised rubber under various types of deformation,” Transactions of the Faraday Society, vol. 40, no. 2, pp. 59–70, 1944.

X.-B. Li and Y.-T. Wei, “An improved Yeoh constitutive model for hyperelastic material,” Engineering Mechanics, vol. 33, no. 12, pp. 38–43, 2016 (Chinese).

X. Hu, X. Liu, M. Li et al., “Selection strategies of hyperelastic constitutive models for carbon black filled rubber,” Engineering Mechanics, vol. 31, no. 05, pp. 34–42, 2014 (Chinese).

B. Chen, Research on Static/Dynamic Characteristics of Suspension Rubber Bushings and Its Application, Southwest Jiaotong University, Chengdu, China, 2014.

H. Lee S, J. Shin K, S. Msolli et al., “Prediction of the dynamic equivalent stiffness for a rubber bushing using the finite element method and empirical modeling,” International Journal of Mechanics and Materials in Design, vol. 15, no. 1, pp. 77–91, 2019.