Yeoh (hyperelastic model)

The Yeoh hyperelastic material modelYeoh, O. H., 1993, "Some forms of the strain energy function for rubber," Rubber Chemistry and technology, Volume 66, Issue 5, November 1993, Pages 754-771.] is a phenomenological model for the deformation of nearly incompressible, nonlinear elastic materials such as rubber. The model is based on Ronald Rivlin's observation that the elastic properties of rubber may be described using a strain energy density function which is a power series in the strain invariants $I\_1,\; I\_2,\; I\_3$. [Rivlin, R. S., 1948, "Some applications of elasticity theory to rubber engineering", in Collected Papers of R. S. Rivlin vol. 1 and 2, Springer, 1997.] The Yeoh model for incompressible rubber is a function only of $I\_1$. For compressible rubbers, an dependence on $I\_3$ is added on. Since a polynomial form of the strain energy density function is used but all the three invariants of the left Cauchy-Green deformation tensor are not, the Yeoh model is also called the reduced polynomial model.

Yeoh model for incompressible rubbers

The original model proposed by Yeoh had a cubic form with only $I\_1$ dependence and is applicable to purely incompressible materials. The strain energy density for this model is written as:$W\; =\; sum\_\{i=1\}^3\; C\_i~(I\_1-3)^i$where $C\_i$ are material constants. The quantity $2\; C\_1$ can be interpreted as the initial shear modulus.

Today a slightly more generalized version of the Yeoh model is used.Selvadurai, A. P. S., 2006, "Deflections of a rubber membrane", Journal of the Mechanics and Physics of Solids, vol. 54, no. 6, pp. 1093-1119.] This model includes $n$ terms and is written as:$W\; =\; sum\_\{i=1\}^n\; C\_i~(I\_1-3)^i\; ~.$

When $n=1$ the Yeoh model reduces to the neo-Hookean model for incompressible materials.

The Cauchy stress for the incompressible Yeoh model is given by:

Uniaxial extension

For uniaxial extension in the $mathbf\{n\}\_1$-direction, the principal stretches are $lambda\_1\; =\; lambda,~\; lambda\_2=lambda\_3$. From incompressibility $lambda\_1~lambda\_2~lambda\_3=1$. Hence $lambda\_2^2=lambda\_3^2=1/lambda$. Therefore,:$I\_1\; =\; lambda\_1^2+lambda\_2^2+lambda\_3^2\; =\; lambda^2\; +\; cfrac\{2\}\{lambda\}\; ~.$The left Cauchy-Green deformation tensor can then be expressed as:If the directions of the principal stretches are oriented with the coordinate basis vectors, we have:$sigma\_\{11\}\; =\; -p\; +\; 2~lambda^2~cfrac\{partial\; W\}\{partial\; I\_1\}\; ~;~~\; sigma\_\{22\}\; =\; -p\; +\; cfrac\{2\}\{lambda\}~cfrac\{partial\; W\}\{partial\; I\_1\}\; =\; sigma\_\{33\}\; ~.$Since $sigma\_\{22\}\; =\; sigma\_\{33\}\; =\; 0$, we have:$p\; =\; cfrac\{2\}\{lambda\}~cfrac\{partial\; W\}\{partial\; I\_1\}\; ~.$Therefore,:$sigma\_\{11\}\; =\; 2~left(lambda^2\; -\; cfrac\{1\}\{lambda\}\; ight)~cfrac\{partial\; W\}\{partial\; I\_1\}~.$The engineering strain is $lambda-1,$. The engineering stress is :$T\_\{11\}\; =\; sigma\_\{11\}/lambda\; =\; 2~left(lambda\; -\; cfrac\{1\}\{lambda^2\}\; ight)~cfrac\{partial\; W\}\{partial\; I\_1\}~.$

Equibiaxial extension

For equibiaxial extension in the $mathbf\{n\}\_1$ and $mathbf\{n\}\_2$ directions, the principal stretches are $lambda\_1\; =\; lambda\_2\; =\; lambda,$. From incompressibility $lambda\_1~lambda\_2~lambda\_3=1$. Hence $lambda\_3=1/lambda^2,$. Therefore,:$I\_1\; =\; lambda\_1^2+lambda\_2^2+lambda\_3^2\; =\; 2~lambda^2\; +\; cfrac\{1\}\{lambda^4\}\; ~.$The left Cauchy-Green deformation tensor can then be expressed as:If the directions of the principal stretches are oriented with the coordinate basis vectors, we have:$sigma\_\{11\}\; =\; -p\; +\; 2~lambda^2~cfrac\{partial\; W\}\{partial\; I\_1\}\; =\; sigma\_\{22\}\; ~;~~\; sigma\_\{33\}\; =\; -p\; +\; cfrac\{2\}\{lambda^4\}~cfrac\{partial\; W\}\{partial\; I\_1\}\; ~.$Since $sigma\_\{33\}\; =\; 0$, we have:$p\; =\; cfrac\{2\}\{lambda^4\}~cfrac\{partial\; W\}\{partial\; I\_1\}\; ~.$Therefore,:$sigma\_\{11\}\; =\; 2~left(lambda^2\; -\; cfrac\{1\}\{lambda^4\}\; ight)~cfrac\{partial\; W\}\{partial\; I\_1\}\; =\; sigma\_\{22\}\; ~.$The engineering strain is $lambda-1,$. The engineering stress is :$T\_\{11\}\; =\; cfrac\{sigma\_\{11\{lambda\}\; =\; 2~left(lambda\; -\; cfrac\{1\}\{lambda^5\}\; ight)~cfrac\{partial\; W\}\{partial\; I\_1\}\; =\; T\_\{22\}~.$

Planar extension

Planar extension tests are carried out on thin specimens which are constrained from deforming in one direction. For planar extension in the $mathbf\{n\}\_1$ directions with the $mathbf\{n\}\_3$ direction constrained, the principal stretches are $lambda\_1=lambda,\; ~lambda\_3=1$. From incompressibility $lambda\_1~lambda\_2~lambda\_3=1$. Hence $lambda\_2=1/lambda,$. Therefore,:$I\_1\; =\; lambda\_1^2+lambda\_2^2+lambda\_3^2\; =\; lambda^2\; +\; cfrac\{1\}\{lambda^2\}\; +\; 1\; ~.$The left Cauchy-Green deformation tensor can then be expressed as:If the directions of the principal stretches are oriented with the coordinate basis vectors, we have:$sigma\_\{11\}\; =\; -p\; +\; 2~lambda^2~cfrac\{partial\; W\}\{partial\; I\_1\}\; ~;~~\; sigma\_\{11\}\; =\; -p\; +\; cfrac\{2\}\{lambda^2\}~cfrac\{partial\; W\}\{partial\; I\_1\}\; ~;~~\; sigma\_\{33\}\; =\; -p\; +\; 2~cfrac\{partial\; W\}\{partial\; I\_1\}\; ~.$Since $sigma\_\{22\}\; =\; 0$, we have:$p\; =\; cfrac\{2\}\{lambda^2\}~cfrac\{partial\; W\}\{partial\; I\_1\}\; ~.$Therefore,:$sigma\_\{11\}\; =\; 2~left(lambda^2\; -\; cfrac\{1\}\{lambda^2\}\; ight)~cfrac\{partial\; W\}\{partial\; I\_1\}\; ~;~~\; sigma\_\{22\}\; =\; 0\; ~;~~\; sigma\_\{33\}\; =\; 2~left(1\; -\; cfrac\{1\}\{lambda^2\}\; ight)~cfrac\{partial\; W\}\{partial\; I\_1\}~.$The engineering strain is $lambda-1,$. The engineering stress is :$T\_\{11\}\; =\; cfrac\{sigma\_\{11\{lambda\}\; =\; 2~left(lambda\; -\; cfrac\{1\}\{lambda^3\}\; ight)~cfrac\{partial\; W\}\{partial\; I\_1\}~.$

Yeoh model for compressible rubbers

A version of the Yeoh model that includes $I\_3\; =\; J^2$ dependence is used for compressible rubbers. The strain energy density function for this model is written as:where , and $C\_\{i0\},\; C\_\{k1\}$ are material constants. The quantity $C\_\{10\}$ is interpreted as half the initial shear modulus, while $C\_\{11\}$ is interpreted as half the initial bulk modulus.

When $n=1$ the compressible Yeoh model reduces to the neo-Hookean model for compressible materials.

References