In continuum mechanics, an Arruda–Boyce model[1] is a hyperelastic constitutive model used to describe the mechanical behavior of rubber and other polymeric substances. This model is based on the statistical mechanics of a material with a cubic representative volume element containing eight chains along the diagonal directions. The material is assumed to be incompressible. The model is named after Ellen Arruda and Mary Cunningham Boyce, who published it in 1993.[1]
The strain energy density function for the incompressible Arruda–Boyce model is given by[2]
![{\displaystyle W=Nk_{B}\theta {\sqrt {n}}\left[\beta \lambda _{\text{chain}}-{\sqrt {n}}\ln \left({\cfrac {\sinh \beta }{\beta }}\right)\right],}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aef39e94e62c561743e1881cb133ea837448d664)
where
is the number of chain segments,
is the Boltzmann constant,
is the temperature in kelvins,
is the number of chains in the network of a cross-linked polymer,
![{\displaystyle \lambda _{\mathrm {chain} }={\sqrt {\tfrac {I_{1}}{3}}},\quad \beta ={\mathcal {L}}^{-1}\left({\cfrac {\lambda _{\text{chain}}}{\sqrt {n}}}\right),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/74812c32359b8f62d32f4caf83db31119ddc2580)
where
is the first invariant of the left Cauchy–Green deformation tensor, and
is the inverse Langevin function which can be approximated by
![{\displaystyle {\mathcal {L}}^{-1}(x)={\begin{cases}1.31\tan(1.59x)+0.91x&{\text{for}}\ |x|<0.841,\\{\tfrac {1}{\operatorname {sgn}(x)-x}}&{\text{for}}\ 0.841\leq |x|<1.\end{cases}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/82d56373f76fecb2776d65779b88b44d961c49a1)
For small deformations the Arruda–Boyce model reduces to the Gaussian network based neo-Hookean solid model. It can be shown[3] that the Gent model is a simple and accurate approximation of the Arruda–Boyce model.
Alternative expressions for the Arruda–Boyce model[]
An alternative form of the Arruda–Boyce model, using the first five terms of the inverse Langevin function, is[4]
![W = C_1\left[\tfrac{1}{2}(I_1-3) + \tfrac{1}{20N}(I_1^2 -9) + \tfrac{11}{1050N^2}(I_1^3-27) + \tfrac{19}{7000N^3}(I_1^4-81) + \tfrac{519}{673750N^4}(I_1^5-243)\right]](https://wikimedia.org/api/rest_v1/media/math/render/svg/d2125f7e62d9a724faae733789a489b78c9cfe9d)
where
is a material constant. The quantity
can also be interpreted as a measure of the limiting network stretch.
If
is the stretch at which the polymer chain network becomes locked, we can express the Arruda–Boyce strain energy density as
![W = C_1\left[\tfrac{1}{2}(I_1-3) + \tfrac{1}{20\lambda_m^2}(I_1^2 -9) + \tfrac{11}{1050\lambda_m^4}(I_1^3-27) + \tfrac{19}{7000\lambda_m^6}(I_1^4-81) + \tfrac{519}{673750\lambda_m^8}(I_1^5-243)\right]](https://wikimedia.org/api/rest_v1/media/math/render/svg/ff58a4cfd63162b02ce61a6ff2d9ec9fbd205b30)
We may alternatively express the Arruda–Boyce model in the form
![{\displaystyle W=C_{1}~\sum _{i=1}^{5}\alpha _{i}~\beta ^{2i-2}~(I_{1}^{i}-3^{i})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/95598773dfc6f9002b03d04981e9b2c763c306dd)
where
and
If the rubber is compressible, a dependence on
can be introduced into the strain energy density;
being the deformation gradient. Several possibilities exist, among which the Kaliske–Rothert[5] extension has been found to be reasonably accurate. With that extension, the Arruda-Boyce strain energy density function can be expressed as
![W = D_1\left(\tfrac{J^2-1}{2} - \ln J\right) + C_1~\sum_{i=1}^5 \alpha_i~\beta^{i-1}~(\overline{I}_1^i-3^i)](https://wikimedia.org/api/rest_v1/media/math/render/svg/1877d1f641f966c84b786c743269e7b2c619e1cd)
where
is a material constant and
. For consistency with linear elasticity, we must have
where
is the bulk modulus.
Consistency condition[]
For the incompressible Arruda–Boyce model to be consistent with linear elasticity, with
as the shear modulus of the material, the following condition has to be satisfied:
![\cfrac{\partial W}{\partial I_1}\biggr|_{I_1=3} = \frac{\mu}{2} \,.](https://wikimedia.org/api/rest_v1/media/math/render/svg/e7455d028046b95bc1bfb49efab7c9beeb0fe89d)
From the Arruda–Boyce strain energy density function, we have,
![\cfrac{\partial W}{\partial I_1} = C_1~\sum_{i=1}^5 i~\alpha_i~\beta^{i-1}~I_1^{i-1} \,.](https://wikimedia.org/api/rest_v1/media/math/render/svg/9f84e04e2d0b0e9a053c6ff8082f4f9515febd00)
Therefore, at
,
![\mu = 2C_1~\sum_{i=1}^5 i\,\alpha_i~\beta^{i-1}~I_1^{i-1} \,.](https://wikimedia.org/api/rest_v1/media/math/render/svg/b256b49fe84f2eef995be763e454aebc988e31e6)
Substituting in the values of
leads to the consistency condition
![\mu = C_1\left(1 + \tfrac{3}{5\lambda_m^2} + \tfrac{99}{175\lambda_m^4} + \tfrac{513}{875\lambda_m^6} + \tfrac{42039}{67375\lambda_m^8}\right) \,.](https://wikimedia.org/api/rest_v1/media/math/render/svg/c33a127b6e45b8db983458a5c373a920f5b6a8a9)
Stress-deformation relations[]
The Cauchy stress for the incompressible Arruda–Boyce model is given by
![\boldsymbol{\sigma} = -p~\boldsymbol{\mathit{1}} +
2~\cfrac{\partial W}{\partial I_1}~\boldsymbol{B}
= -p~\boldsymbol{\mathit{1}} + 2C_1~\left[\sum_{i=1}^5 i~\alpha_i~\beta^{i-1}~I_1^{i-1}\right]\boldsymbol{B}](https://wikimedia.org/api/rest_v1/media/math/render/svg/81a3c389469fc026045d37766422646d1102044a)
Uniaxial extension[]
Stress-strain curves under uniaxial extension for Arruda–Boyce model compared with various hyperelastic material models.
For uniaxial extension in the
-direction, the principal stretches are
. From incompressibility
. Hence
.
Therefore,
![I_1 = \lambda_1^2+\lambda_2^2+\lambda_3^2 = \lambda^2 + \cfrac{2}{\lambda} ~.](https://wikimedia.org/api/rest_v1/media/math/render/svg/4b87bfcd69908379ae27cff7cd9dc61d3ea46051)
The left Cauchy–Green deformation tensor can then be expressed as
![\boldsymbol{B} = \lambda^2~\mathbf{n}_1\otimes\mathbf{n}_1 + \cfrac{1}{\lambda}~(\mathbf{n}_2\otimes\mathbf{n}_2+\mathbf{n}_3\otimes\mathbf{n}_3) ~.](https://wikimedia.org/api/rest_v1/media/math/render/svg/eb5b7a04b2de4cd2ba759594bb5b41df2b6cdd27)
If the directions of the principal stretches are oriented with the coordinate basis vectors, we have
![\begin{align}
\sigma_{11} & = -p + 2C_1\lambda^2\left[\sum_{i=1}^5 i~\alpha_i~\beta^{i-1}~I_1^{i-1}\right] \\
\sigma_{22} & = -p + \cfrac{2C_1}{\lambda}\left[\sum_{i=1}^5 i~\alpha_i~\beta^{i-1}~I_1^{i-1}\right] = \sigma_{33} ~.
\end{align}](https://wikimedia.org/api/rest_v1/media/math/render/svg/77ac3dbb74bf11083fdb2a0ac8d8dd06ee6db556)
If
, we have
![p = \cfrac{2C_1}{\lambda}\left[\sum_{i=1}^5 i~\alpha_i~\beta^{i-1}~I_1^{i-1}\right]~.](https://wikimedia.org/api/rest_v1/media/math/render/svg/86ad1c5c872b44f7d361684b586dce991f30aadc)
Therefore,
![\sigma_{11} = 2C_1\left(\lambda^2 - \cfrac{1}{\lambda}\right)\left[\sum_{i=1}^5 i~\alpha_i~\beta^{i-1}~I_1^{i-1}\right]~.](https://wikimedia.org/api/rest_v1/media/math/render/svg/04c4319b65bb9600950159f3d70435705f32d0f0)
The engineering strain is
. The engineering stress is
![T_{11} = \sigma_{11}/\lambda =
2C_1\left(\lambda - \cfrac{1}{\lambda^2}\right)\left[\sum_{i=1}^5 i~\alpha_i~\beta^{i-1}~I_1^{i-1}\right]~.](https://wikimedia.org/api/rest_v1/media/math/render/svg/bfac3b0a1d62886788be468f825931987e357b16)
Equibiaxial extension[]
For equibiaxial extension in the
and
directions, the principal stretches are
. From incompressibility
. Hence
.
Therefore,
![I_1 = \lambda_1^2+\lambda_2^2+\lambda_3^2 = 2~\lambda^2 + \cfrac{1}{\lambda^4} ~.](https://wikimedia.org/api/rest_v1/media/math/render/svg/6c83b5d6d6838e31ccc5feceaf70b7a69d341b38)
The left Cauchy–Green deformation tensor can then be expressed as
![\boldsymbol{B} = \lambda^2~\mathbf{n}_1\otimes\mathbf{n}_1 + \lambda^2~\mathbf{n}_2\otimes\mathbf{n}_2+ \cfrac{1}{\lambda^4}~\mathbf{n}_3\otimes\mathbf{n}_3 ~.](https://wikimedia.org/api/rest_v1/media/math/render/svg/e604a880183bac06582e027580d961d338573016)
If the directions of the principal stretches are oriented with the coordinate basis vectors, we have
![\sigma_{11} = 2C_1\left(\lambda^2 - \cfrac{1}{\lambda^4}\right)\left[\sum_{i=1}^5 i~\alpha_i~\beta^{i-1}~I_1^{i-1}\right] = \sigma_{22} ~.](https://wikimedia.org/api/rest_v1/media/math/render/svg/59dffd2edd4b477b7670b70e4a78afa491a088f4)
The engineering strain is
. The engineering stress is
![T_{11} = \cfrac{\sigma_{11}}{\lambda} =
2C_1\left(\lambda - \cfrac{1}{\lambda^5}\right)\left[\sum_{i=1}^5 i~\alpha_i~\beta^{i-1}~I_1^{i-1}\right] = T_{22}~.](https://wikimedia.org/api/rest_v1/media/math/render/svg/b83d49645584c6cf53bfcfeac65548a85814a1c2)
Planar extension[]
Planar extension tests are carried out on thin specimens which are constrained from deforming in one direction. For planar extension in the
directions with the
direction constrained, the principal stretches are
. From incompressibility
. Hence
.
Therefore,
![I_1 = \lambda_1^2+\lambda_2^2+\lambda_3^2 = \lambda^2 + \cfrac{1}{\lambda^2} + 1 ~.](https://wikimedia.org/api/rest_v1/media/math/render/svg/474f46d99397c3e029d63d87c7e5d9b77193ebf8)
The left Cauchy–Green deformation tensor can then be expressed as
![\boldsymbol{B} = \lambda^2~\mathbf{n}_1\otimes\mathbf{n}_1 + \cfrac{1}{\lambda^2}~\mathbf{n}_2\otimes\mathbf{n}_2+ \mathbf{n}_3\otimes\mathbf{n}_3 ~.](https://wikimedia.org/api/rest_v1/media/math/render/svg/6f2925db738eebc7cd11100bb5aee02cb1294be9)
If the directions of the principal stretches are oriented with the coordinate basis vectors, we have
![\sigma_{11} = 2C_1\left(\lambda^2 - \cfrac{1}{\lambda^2}\right)\left[\sum_{i=1}^5 i~\alpha_i~\beta^{i-1}~I_1^{i-1}\right] ~;~~ \sigma_{22} = 0 ~;~~ \sigma_{33} = 2C_1\left(1 - \cfrac{1}{\lambda^2}\right)\left[\sum_{i=1}^5 i~\alpha_i~\beta^{i-1}~I_1^{i-1}\right]~.](https://wikimedia.org/api/rest_v1/media/math/render/svg/d2c74598ecc266aad14a5e46f2e95c0b78a42437)
The engineering strain is
. The engineering stress is
![T_{11} = \cfrac{\sigma_{11}}{\lambda} =
2C_1\left(\lambda - \cfrac{1}{\lambda^3}\right)\left[\sum_{i=1}^5 i~\alpha_i~\beta^{i-1}~I_1^{i-1}\right]~.](https://wikimedia.org/api/rest_v1/media/math/render/svg/45b5df0a5f6efcbdfb033f33ef33226837031e3a)
Simple shear[]
The deformation gradient for a simple shear deformation has the form[6]
![\boldsymbol{F} = \boldsymbol{1} + \gamma~\mathbf{e}_1\otimes\mathbf{e}_2](https://wikimedia.org/api/rest_v1/media/math/render/svg/dbfda261d2e46796ad6df5c553fef8cdf6ff2c96)
where
are reference orthonormal basis vectors in the plane of deformation and the shear deformation is given by
![\gamma = \lambda - \cfrac{1}{\lambda} ~;~~ \lambda_1 = \lambda ~;~~ \lambda_2 = \cfrac{1}{\lambda} ~;~~ \lambda_3 = 1](https://wikimedia.org/api/rest_v1/media/math/render/svg/527c35f38d8ba0455e176c0982e5ddef58cdfabb)
In matrix form, the deformation gradient and the left Cauchy–Green deformation tensor may then be expressed as
![\boldsymbol{F} = \begin{bmatrix} 1 & \gamma & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{bmatrix} ~;~~
\boldsymbol{B} = \boldsymbol{F}\cdot\boldsymbol{F}^T = \begin{bmatrix} 1+\gamma^2 & \gamma & 0 \\ \gamma & 1 & 0 \\ 0 & 0 & 1 \end{bmatrix}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5ed7a851026e0a6e9d9122ed7548721b6f2aecb5)
Therefore,
![I_1 = \mathrm{tr}(\boldsymbol{B}) = 3 + \gamma^2](https://wikimedia.org/api/rest_v1/media/math/render/svg/6c4b57d8446643ef4c59ec9079abeb7b7579edf2)
and the Cauchy stress is given by
![\boldsymbol{\sigma} = -p~\boldsymbol{\mathit{1}} + 2C_1\left[\sum_{i=1}^5 i~\alpha_i~\beta^{i-1}~(3+\gamma^2)^{i-1}\right]~\boldsymbol{B}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b09423828cffc80c66de470d08d3bd90609ad570)
Statistical mechanics of polymer deformation[]
The Arruda–Boyce model is based on the statistical mechanics of polymer chains. In this approach, each macromolecule is described as a chain of
segments, each of length
. If we assume that the initial configuration of a chain can be described by a random walk, then the initial chain length is
![r_0 = l\sqrt{N}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ff6b78e0cad396268c2f6ee2775c38cee71ba680)
If we assume that one end of the chain is at the origin, then the probability that a block of size
around the origin will contain the other end of the chain,
, assuming a Gaussian probability density function, is
![p(x_1,x_2,x_3) = \cfrac{b^3}{\pi^{3/2}}~\exp[-b^2(x_1^2 + x_2^2 + x_3^2)] ~;~~ b := \sqrt{\cfrac{3}{2Nl^2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b77b842670b13db40bcbe6aa79f2279774613b2d)
The configurational entropy of a single chain from is
![s = c -k_B b^2 r^2](https://wikimedia.org/api/rest_v1/media/math/render/svg/ce9589bb79571605832723eb7f13895982a38a8d)
where
is a constant. The total entropy in a network of
chains is therefore
![\Delta S = -\tfrac{1}{2} n k_B (\lambda_1^2 + \lambda_2^2 + \lambda_3^2 - 3) = -\tfrac{1}{2} n k_B (I_1-3)](https://wikimedia.org/api/rest_v1/media/math/render/svg/c6f9d6d5d9117c66fc4649fa70c84212b97c7009)
where an affine deformation has been assumed. Therefore the strain energy of the deformed network is
![W = -\theta \, dS = \tfrac{1}{2} n k_B \theta (I_1-3)](https://wikimedia.org/api/rest_v1/media/math/render/svg/ff9fd844e9e488d4c9ea98a0a6c92448234d258e)
where
is the temperature.
Notes and references[]
- ^ a b Arruda, E. M. and Boyce, M. C., 1993, A three-dimensional model for the large stretch behavior of rubber elastic materials,, J. Mech. Phys. Solids, 41(2), pp. 389–412.
- ^ Bergstrom, J. S. and Boyce, M. C., 2001, Deformation of Elastomeric Networks: Relation between Molecular Level Deformation and Classical Statistical Mechanics Models of Rubber Elasticity, Macromolecules, 34 (3), pp 614–626, doi:10.1021/ma0007942.
- ^ Horgan, C. O. and Saccomandi, G., 2002, A molecular-statistical basis for the Gent constitutive model of rubber elasticity, Journal of Elasticity, 68(1), pp. 167–176.
- ^ Hiermaier, S. J., 2008, Structures under Crash and Impact, Springer.
- ^ Kaliske, M. and Rothert, H., 1997, On the finite element implementation of rubber-like materials at finite strains, Engineering Computations, 14(2), pp. 216–232.
- ^ Ogden, R. W., 1984, Non-linear elastic deformations, Dover.
See also[]