Research Article  Open Access
Modeling and Verification of a New Hyperelastic Model for RubberLike Materials
Abstract
The essential evaluation criterion for the hyperelastic model is its ability to describe the mechanical behavior of rubberlike 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 rubberlike 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 vibrationisolation devices [1–3]. Without considering the time effect, rubberlike 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 rubberlike materials. The essential evaluation criterion for the hyperelastic model is its ability to reproduce the mechanical behavior of rubberlike 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 rubberlike materials in all deformation modes by using only a certain type of test data. In addition, these models cannot also perfectly reproduce the stressstrain 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 threedimensional electroelastic constitutive framework that can model the stiffness gaining during the curing process undergoing finite deformations. By improving MooneyRivlin 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 sixparameter model. In contrast to the traditional model, Crespo et al. [8] proposed a new hyperelastic model by WhatYouPrescribeIsWhatYouGet 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 networkingaveraging of the tube model and by applying a closedform of the Rayleigh exact nonGaussian 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 rubberlike 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 ArrudaBoyce model in predicting the equibiaxial data, Hossain et al. [12] compared five modified versions of the ArrudaBoyce model and two modified versions of fullnetwork model with the ArrudaBoyce model. Bahreman and Darijani [13] proposed a new polynomial hyperelastic model which is a function of the principal invariants of the left CauchyGreen strain tensor. Based on the neoHookean 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 FungDemiray 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 CauchyGreen deformation tensor can more accurately reflect the mechanical behavior and some physical effects of rubberlike materials. Under the incompressible assumption, Lin [16] derived the mathematical relationship between the right CauchyGreen deformation tensor principal invariants and used the interpolation algorithm to obtain a new hyperelastic model that can predict stressstrain 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 rubberlike 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 rubberlike 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 stressstrain relationship of rubberlike materials can be derived from the strain energy function . The relationship between the first PiolaKirchhoff 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 CauchyGreen deformation tensor . In general, rubberlike 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 stressstrain 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 rubberliked 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 stressstrain relationship of rubberlike materials is significantly nonlinear over the large deformation range. The hyperelastic models which do not contain correlation term of the second invariant I_{2} 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 rubberlike 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 ValanisLandel hypothesis [21]. The relationships between the principal invariants of the right CauchyGreen deformation tensor and the principal stretches are as follows.
When and , (11) is the neoHookean model.
When and , (11) is the MooneyRivlin 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 rubberlike 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 CauchyGreen 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 stressstrain 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 stressstrain 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.

(a)
(b)
(c)
In order to evaluate the goodness of fit of each model, the coefficients of determination R^{2} 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)
(a)
(b)
(c)
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_{2} 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 stressstrain 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.

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 highorder terms of invariant I_{1}, and the orders of invariant I_{2} are different. The order of the I_{2} terms in the Carroll model is 1/2, and the highest order of the I_{2} 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 rubberlike 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 rubberlike 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.
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^{−2} MPa, = 7.9451×10^{−1} MPa, = 5.6729×10^{−2} MPa, = 1.6134×10^{−3} 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 Yaxis, respectively. The simulation results are shown in Figure 5.
(a)
(b)
(c)
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 loadingunloading process is repeated four times, and then the fifth loading curve is taken as the final results.
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)
6. Conclusion
In this paper, a hyperelastic model for an isotropic and incompressible rubberlike materials is proposed based on the Seth strain tensor invariants. The new model is a function of the invariants of the right CauchyGreen deformation tensor, satisfying the ValanisLandel 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 highorder terms of invariants I_{1} and I_{2}, while the Yeoh model does not contain the terms of I_{2} and the Carroll model does not contain the highorder terms of I_{2}. 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 nonpneumatic tire,” Key Engineering Materials, vol. 775, pp. 554–559, 2018. View at: Publisher Site  Google Scholar
 C. Cheng, S. Li, Y. Wang, and X. Jiang, “Force and displacement transmissibility of a quasizero stiffness vibration isolator with geometric nonlinear damping,” Nonlinear Dynamics, vol. 87, no. 4, pp. 2267–2279, 2017. View at: Publisher Site  Google Scholar
 L. P. Li, “Experiment analysis about mechanical properties of rubber bushing for suspension telescopic shock absorber,” Applied Mechanics and Materials, vol. 670671, pp. 1008–1011, 2014. View at: Publisher Site  Google Scholar
 P. Steinmann, M. Hossain, and G. Possart, “Hyperelastic models for rubberlike materials: Consistent tangent operators and suitability for Treloar's data,” Archive of Applied Mechanics, vol. 82, no. 9, pp. 1183–1217, 2012. View at: Publisher Site  Google Scholar
 M. Hossain and P. Steinmann, “More hyperelastic models for rubberlike materials: consistent tangent operators and comparative study,” Journal of the Mechanical Behavior of Materials, vol. 22, no. 12, pp. 27–50, 2013. View at: Publisher Site  Google Scholar
 M. Hossain, “Modelling the curing process in particlefilled electroactive polymers with a dispersion anisotropy,” Continuum Mechanics and Thermodynamics, 2019. View at: Publisher Site  Google Scholar
 K. Yaya and H. Bechir, “A new hyperelastic model for predicting multiaxial behaviour of rubberlike materials: formulation and computational aspects,” Mechanics of TimeDependent Materials, vol. 22, no. 2, pp. 167–186, 2018. View at: Publisher Site  Google Scholar
 J. Crespo, M. Latorre, and F. Montáns, “WYPIWYG hyperelasticity for isotropic, compressible materials,” Computational Mechanics, vol. 59, no. 1, pp. 73–92, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 V. N. Khiêm and M. Itskov, “Analytical networkaveraging of the tube model: rubber elasticity,” Journal of the Mechanics and Physics of Solids, vol. 95, pp. 254–269, 2016. View at: Publisher Site  Google Scholar
 F. Zhao, “Continuum constitutive modeling for isotropic hyperelastic materials,” Advances in Pure Mathematics, vol. 06, no. 09, pp. 571–582, 2016. View at: Publisher Site  Google Scholar
 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 NonLinear Mechanics, vol. 84, pp. 1–11, 2016. View at: Publisher Site  Google Scholar
 M. Hossain, A. Amin, and M. N. Kabir, “Eightchain and fullnetwork models and their modified versions for rubber hyperelasticity: a comparative study,” Journal of the Mechanical Behavior of Materials, vol. 24, no. 12, pp. 11–24, 2015. View at: Publisher Site  Google Scholar
 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. View at: Google Scholar
 H. Bechir, L. Chevalier, M. Chaouche, and K. Boufala, “Hyperelastic constitutive model for rubberlike materials based on the first Seth strain measures invariant,” European Journal of Mechanics  A/Solids, vol. 25, no. 1, pp. 110–124, 2006. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 B. Lin, “A new model for hyperelasticity,” Acta Mechanica, vol. 208, no. 12, pp. 39–53, 2009. View at: Publisher Site  Google Scholar
 H. Xiao, “Elastic potentials with best approximation to rubberlike elasticity,” Acta Mechanica, vol. 226, no. 2, pp. 331–350, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 L. Yuan, Z.X. Gu, Z.N. Yin, and H. Xiao, “New compressible hyperelastic models for rubberlike materials,” Acta Mechanica, vol. 226, no. 12, pp. 4059–4072, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 M. M. Attard and G. W. Hunt, “Hyperelastic constitutive modeling under finite strain,” International Journal of Solids and Structures, vol. 41, no. 1819, pp. 5327–5350, 2004. View at: Publisher Site  Google Scholar
 K. Yaya, H. Bechir, and F. Bremand, “Implementation of new strainenergy density function for a grade of carbon blackfilled natural rubber in finite element code,” in Proceedings of the 6th International Conference on Advances in Mechanical Engineering and Mechanics (ICAMEM2015), Hammamet, Tunisia, 2015. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 K. Farahani and H. Bahai, “Hyperelastic constitutive equations of conjugate stresses and strain tensors for the SethHill strain measures,” International Journal of Engineering Science, vol. 42, no. 1, pp. 29–41, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 T. Beda, “An approach for hyperelastic modelbuilding and parameters estimation a review of constitutive models,” European Polymer Journal, vol. 50, no. 1, pp. 97–108, 2014. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar  MathSciNet
 O. H. Yeoh, “Some forms of the strain energy function for rubber,” Rubber Chemistry and Technology, vol. 66, no. 5, pp. 754–771, 1993. View at: Publisher Site  Google Scholar
 L. R. G. Treloar, “Stressstrain data for vulcanised rubber under various types of deformation,” Transactions of the Faraday Society, vol. 40, no. 2, pp. 59–70, 1944. View at: Publisher Site  Google Scholar
 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). View at: Google Scholar
 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). View at: Google Scholar
 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. View at: Google Scholar
Copyright
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.