Vol. 3  No. 8 (2011) , Article ID: 6824 , 5 pages DOI:10.4236/eng.2011.38098

Finite Deformation and Viscoelasticity Modeling and Test

Tibi Beda1,2, Yvon Chevalier3, Kokou-Esso Atcholi4, Essole Padayodi4, Jean-Claude Sagot4

1Industrial and Mechanical Engineering Department, Ecole Nationale Supérieure Polytechnique, Yaounde, Cameroon

2Département de Physique, Faculté des Sciences, Université de Ngaoundéré, Ngaoundéré, Cameroon

3Laboratoire d’Ingénierie des Structures Mécanique et des Matériaux, Saint Ouen, France

4 Laboratoire Systèmes & Transport (SeT EA 3317), Université de Technologie de Belfort-Montbéliard, Belfort Cedex,France


{kokou-esso.atcholi, essole.padayodi, jean-claude.sagot}

Received January 19, 2011; revised July 10, 2011; accepted July 20, 2011

Keywords: Elastomeric materials, Nonlinear behavior, Material testing, Compressibility modeling, BKZ model, Mechanical properties


A model is considered as a representation of compressive and incompressive elastomeric materials in nonlinear behavior. Applications are done on one hand by the characterisation of polyurethane 60 - 65 shore A (a compressive material), and on the other hand by the characterisation of polyurethane 95 shore A and fluorosilicone, both incompressive materials. The Rivlin energy expression is used for incompressive materials. Linear vibrations superposed on static large deformation, which is most often the real using state of elastomeric materials, are studied. Relative experimental and numerical results presen-ted show good predictions.

1. Introduction

Elastomeric materials (Rubber, polyurethane, fluorosilicone) are widely used in industry over the past decades; it shows evidence of progress in this field of material researches. Rubberlike materials, compact or porous, have very complex mechanical behavior: for a small loading, this kind of material presents large strains (many percents); stress depends on strain velocity, incompressibility and temperature. Reason why, the thermomechanical behavior of these materials, take into account nonlinear viscoelastic effects: geometrical nonlinearity and nonlinear constitutive law.

The nonlinear behavior of materials can be described by a lot of models based on undeformed configuration (Lagrangien description) in which, behavior of hyperelastic materials needs relationship between Piola-Kirchhoff stress tensor S and Green-Lagrange strain tensor E. In deformed configuration (Eulerian description) mechanical behavior of materials is described by relationship between the Cauchy stress tensor [1] and the strain tensor of Almansi A. The first standpoint is useful in prediction and computational approach of mechanical systems whereas the second standpoint is used for experimental analysis of materials.

We tried to present mechanical behavior of rubber materials in the frame of non-linear viscoelasticity under small compressibility and large deformations. The final model is a combination of hyperelasticity and viscoelasticity. We pointed out two directions based on tests when using this kind of model: relaxation and vibration test under various loading.

2. Modeling Nonlinear Elasticity

A lot of models can be used to describe non-linear viscoelastic behavior of materials. The common model is a mathematical model based on multiple convolutions. The stress tensor of Piola-Kirchhoff S is related to strain Green-Lagrange tensor E by Volterra-Frechet series [2-5]. This is a very complex approach which requires special tests on materials to determine kernel of each rank. Lai and Findley [2] propose parabolic kernel (in time) for rank 1, 2 and 3. We take up a more pragmatic way to describe viscoelastic behavior of material [6,7] based on linear approach. In Schapery model [7], nonlinearities are expressed in term of variable multiplying factors (stress, strain, times …). The model presented is close to O’Down and Knauss [8]; and Valanis and Landel model [9] using single convolution to describe memory effect of material while elasticity is expressed in term of hyperelasticity. The Piola-Kirchhoff stress tensor S(t) = is related to Green-Lagrange strain tensor =, see Beda and Chevalier works [10], by the following equations for positive values of time :


where w is the strain energy density and

is the fourth rank relaxation tensor (independent functions for isotropic materials). The constitutive relations (1) are consistent with previous model if the following compatibility relations between strain energy density and relaxation functions are satisfied [10].


The nonlinear viscoelasticity described by relation (1) is a general formulation and can be widely applied. The model can be update by two classical viscoelastic tests: transient test and vibration test under large deformations.

3. Transient Relaxation Test

To take into account the compressibility of isotropic rubberlike materials the strain energy density is based on B.K.Z (Bernstein, Kearsley, Zapas) model:




are the three invariants of the right Cauchy strain tensor, where F is the deformation gradient operator, Ln a parameter.

allows us to bring back to incompressible material [11], is the penalty function and is a coupling factor between compressive and incompressible materials. When analyzing polyurethane materials, the following modeling (Rivlin hyperelasticity law) is held:


with M ≤ 3 and Aij material characteristics related to the distortional response.

The penalty function depends on bulk modulus K of material; it means that the relative variation of pressure is proportional to relative volume variation.


The coupling function Q is given in terms of series:


where coefficients an are determined by tensile test and mathematical regression [12,13], N ≤ 2.

Table 1 relates Rivlin constants for polyurethane given by using uniaxial tensile tests [14].

For compressible material (polyurethane 60 - 65 shore A) the coupling function is given by curves, Figure 1, with bulk modulus K ≈ 23 MPa.

Table 1. First part of energy deformation density of polyurethane. Result of Uniaxial tensile test under elongation λ (1 ≤ λ ≤ 1.2).

Figure 1. Coupling function Q versus invariant of right Cauchy Green tensor.

For incompressible materials (polyurethane 95 shore A) function Q and G can be neglected in relation (3).

Concerning the memory effect of material, we carried out tensile relaxation tests and retained for relaxation function usual Prony’s series model in which parameters depend on Green-Lagrange tensor. In the case of tensile test, the relaxation tensor R(l, t) is given by the following expression:


where λ is the stretch.

Tensile relaxation tests are carried out under four different stretches (λ = 1.05, λ = 1.10, λ = 1.15, λ = 1.20).

One remark that for the polyurethane 95 shore A, an incompressive material, the second order model describes quite well the material behavior prediction.

Figure 2 relates the evolution of magnitude factors and time constants versus stretch λ. It shows that the stretch λ weakens stiffness of material (decreases with increasing λ) and reinforces the elasticity (t increases with increasing λ).

Let us examine the consistency of the model and the compatibility, relations (2).

Table 2 shows that assessment of the compatibility, relation (2) is satisfactory (discrepancy ≈ 5%). Finally we try a comparison between the actual Cauchy stress tensor s and the Cauchy stress tensor predicted by the model at two strain velocities.

Figure 3 shows that the model gives a better prediction at low strain velocity which is close to static behavior.

Figure 2. Variation of magnitude factor (ri) and time constant (ti), (i = 1, 2, 3 relation 8), versus elongation λ for polyurethane 95 shore A.

Table 2. Comparison between static and relaxation behavior of polyurethane 95 shore A. Validity of compatibility.

(a) (b)

Figure 3. Comparison between model and test at various strain velocities.

4. Vibrations Tests

In industrial context, the great majority of rubber materials are working under static load (large strain) on where small vibrations ε(t) are superposed.

See for example rubber connection in vibrating structures (rotating engine, car engine, etc. …).

The total strain is then given by the following relations:


in which.

Constitutive relations (1) are expanded till first order:

and the Piola-Kirchhoff stress tensor is a sum of two stress tensors S0(t) and S1(t) given by relations (10) and (11):



S0(t) is static stress, constant for large value of time while S1(t) relates linear viscoelastic behavior around large strain.

Experimental setups involving longitudinal vibration of beams [15-17] are presented in Figure 4.

Figure 5 shows experimental curves of Young modulus of fluorosilicone rubber material. Analytic formulas for constitutive equations of viscoelastic materials are provided to predict behavior of material under large range of frequencies and different stretches. To describe the behavior of realistic viscoelastic material, a fractional derivative model [18-20] is used. The complex Young modulus can be written in the following form:

Figure 4. Experimental setups and data recording system used to evalu-ate complex Young modulus of rubber material.

Figure 5. Amplitude of complex Young modulus of a fluorosilicone rubber material.


where E0, f0, f1 and α are material parameters depending of the stretch l. Numerical values of parameters are drawn out of tensile vibration test, Figure 6, and are exhibited in Table 3.

Results presented by forward-looking updating technique (tests and fractional derivative model) are in keeping with viscoelastic properties of materials [21]: the Young modulus gain increases with increasing loading while the damping factor decreases (the frequency of maximum damping increases with increasing loading).

5. Concluding Remarks

The whole interest of this nonlinear viscoelastic model is in engineering sciences: model can predict the effect of static loading of viscoelastic behavior of materials. Experimental complex Young modulus and complex Coulomb modulus can be obtained over a large frequencies range, to include the static and high frequencies moduli under various loading. Furthermore, the effect of loading is quantized in relaxation behavior of material.

Figure 6. Young modulus of fluorosilicone rubber material versus frequen-cy at various modeling. Curves are drawn from relation (12) and updated (Table 3).

Table 3. numerical values of fractional derivative model of fluorosilicone rubber material under static load.

Some problems are open-ended. Is it allowed to use only lagrangian configuration? And vibrations tests? Finally is the model thermodynamically correct (by considering negative energetic constants)? These are among others perspectives to be looking at for sharpen the results awaited.


  1. J. Salencon, “Mécanique des Milieux Continus,” Editions Ellispe, Paris, 1988.
  2. J. S. Lai and N. W. Findley, “Stress Relaxation of Nonlinear Viscoelastic Material under Uniaxial Strain,” Transactions of the Society of Rheolog, Vol. 12, No. 2, 1968, pp. 259-280. doi:10.1122/1.549108
  3. F. J. Locket, “Nonlinear Viscoelastic Solids,” Academic Press, San Diego, 1972.
  4. A. Molinari, “Sur la Relaxation Entre Fluage et Relaxation en Viscoélasticité non Linéaire,” Comptes Rendus, Académie des Sciences, Tome 277, Série A, 1973, pp. 621-623.
  5. C. Huet, “Relations between Creep and Relaxation Function in Nonlinear Viscoelasticity with or without Aging,” Journal of Rheology, Vol. 29, No. 3, 1985, pp. 245-257. doi:10.1122/1.549789
  6. B. D. Coleman and W. Noll, “Foundation of Linear Viscoelasticity,” Reviews of Modern Physics, Vol. 33, No. 2, 1961, pp. 239-249. doi:10.1103/RevModPhys.33.239
  7. R. A. Schapery, “On the Characterization of Nonlinear Viscoelastic Materials,” Polymer Engineering Science, Vol. 9, No. 4, 1969, pp. 295-310. doi:10.1002/pen.760090410
  8. N. P. O’Dowd and W. G. Knauss, “Time Dependent Large Principal Deformation of Polymers,” Journal of the Mechanics and Physics of Solids, Vol. 43, No. 5, 1996, pp. 771-792.
  9. K. C. Valanis and R. F. Landel, “Large Axial Deformation Behavior of Filled Rubber,” Transactions of the Society of Rheology, Vol. 11, 1967, pp. 213-256. doi:10.1122/1.549080
  10. T. Beda and Y. Chevalier, “Sur le Comportement Statique et Dynamique des Élastomères en Grandes Déformations,” Mécanique Industrielle et Matériaux, vol. 50, No. 5, 1997, pp. 228-231.
  11. R. W. Ogden, “Nonlinear Elastic Deformation,” Ellis Horwood Edition, Robson, 1984.
  12. T. Beda, Y. Chevalier, H. Gacem and P. Mbarga, “Domain of Validity and Fit of Gent-Thomas and Flory-Erman Rubber Models to Data,” Express Polymer Letters, vol. 2, No. 9, 2008, pp. 615-622. doi:10.3144/expresspolymlett.2008.74
  13. T. Beda, “Combining Approach in Stages with Least Squares for Fits of Data in Hyperelasticity,” Comptes Rendus Mecanique, Académie des Sciences, Elsevier, Paris, Vol. 334, No. 10, 2006, pp. 628-633.
  14. H. Bechir, “Comportement Viscoélastique des Élastomères de Polyuréthane en Grandes Déformations— Modélisation—Validation,” Ph.D. Thesis in Mechanics, CNAM, Paris, 1996.
  15. M. Soula, T. Ving, Y. Chevalier, T. Beda and C. Esteoule, “Measurements of Isothermal Complex Moduli of Viscoelastic Materials over a Large Range of Frequencies,” Journal of Sound and Vibration, Vol. 205, No. 2, 1997, pp. 167-184. doi:10.1006/jsvi.1997.0978
  16. M. Soula, T. Ving and Y. Chevalier, “Transient Responses of Polymers and Elastomers Deduced from Harmonic Responses,” Journal of Sound and Vibration, Vol. 205, No. 2, 1997, pp. 185-203. doi:10.1006/jsvi.1997.0979
  17. Y. Chevalier, T. Beda and J. Merdrignac, “Détermination de Modules d’Young Complexes par Identification,” Mécanique Matériaux, Electricité, GAMI, Vol. 32, No. 431, 1989, pp. 55-60.
  18. T. Beda and Y. Chevalier, “Identification of Viscoelastic Fractional Complex Modulus,” American Institute of Aeronautics and Astronautics Journal, AIAA, Vol. 42, No. 7, 2004, pp. 1450-1456.
  19. Y. Chevalier and T. Beda, “Fractional Derivative in Materials and Structures Vibratory Behaviour: Identification Methods,” Journal Européen des Systèmes Automatisés, Hermès Lavoisier, Vol. 42, no. 6-8, 2008, pp. 863-877.
  20. M. Soula and Y. Chevalier, “La Dérivée Fractionnaire en Rhéologie des Polymères—Application aux Comportements ÉLastiques et Viscoélastiques Linéaires et non Linéaires des Élastomères,” ESAIM Proceeding, Vol. 5, 1998, pp. 193-204. doi:10.1051/proc:1998007
  21. A. D. Nashif, D. I. J. Jones and J. P. Henderson, “Vibration Damping,” Wiley, New York, 1984.