Isotropic, Hyperelastic-Plastic Material

Constitutive Law

This MPM Material is an isotropic, elastic-plastic material in large strains using a hyperelastic formulation. The elastic regime for this material is identical to a Mooney except that it only allows a Neohookean elastic regime (with G = G₁ and G₂ = 0).,

The formulation of finite strain plasticity is based on the notion of a stress free intermediate configuration and uses a multiplicative decomposition of the deformation gradient F given by:

      [math]\displaystyle{ \mathbf{F} = \mathbf{F}_{e}. \mathbf{F}_{p} }[/math]

Where F_e and F_p are the elastic and plastic deformation gradient tensors respectively, with det (F_p), that supposes the plastic flow to be isochoric. The Neo-Hookean elastic stored energy, represented by its uncoupled volumetric-deviatoric internal energy form, is consistent with the fundamental idea that the elastic-plastic deviatoric response is assumed to be uncoupled from the elastic volumetric response

In finite strain plasticity, the stored energy is based on the additive decomposition of the stored energy into elastic W_e and plastic F_p internal energies. The elastic stored energy is related to the intermediate configuration and the plastic stored energy is expressed in term of plastic state variables α.

      [math]\displaystyle{ W =W_{e} (\mathbf{B}_{e}) + W_{p} (\alpha) }[/math]

The stored Neo-Hookean stored energy, W_e, is identical to the Mooney Energy W and dependent on entered small-strain, bulk modulus (κ), small-strain, shear modulus (G = G₁), and dilation energy option (UJOption). The value of G₂ is always zero in this material.

In associative plasticity, the plastic storage energy is represented by the plastic flow condition. The plastic flow model considered here is isotropic hardening. It is handled by any hardening law available in the code (see Hardening Laws). The associative flow rate is defined by the principle of maximum plastic dissipation (Simo J C, 1988a and 1988b). It is given, in the present context, by:

      [math]\displaystyle{ L_{v} \mathbf{B}^e = \mathbf{F} {\delta\over t} (\mathbf{\bar C}{^{p-1}}) \mathbf{F^T} = - {2\over 3} {\gamma} {\rm Tr}(\mathbf{B}_{e}) )\mathbf{n} \qquad {\rm with} \qquad \mathbf{n} = {\mathbf{\tau^{d}}\over ||\mathbf{\tau^{d}}||} }[/math]

Where L_v is the Lie derivative of the deviatoric part of the elastic left Cauchy-Green strain tensor [math]\displaystyle{ \bigl(\mathbf{\bar B}_{e}\bigr) }[/math]. It represents the plastic strain rate that is a tensor normal to the yield surface in the stress space; n is a normal to the yield surface and γ is the consistency parameter also called the plastic multiplicator. In addition, a isotropic hardening law is needed. It is represented by the rate equation, as in the linear theory:

      [math]\displaystyle{ {d {\alpha}\over dt } = (2/3)^{1\over 2} {\gamma} }[/math]

Material Properties

The material properties are set using

Property	Description	Units	Default
E	Elastic modulus	MPa	none
[math]\displaystyle{ {\nu} }[/math]	Poisson's ratio	none	none
alpha	Thermal expansion coefficient	ppm/M	40

History Variables

None

Examples

These commands model polymer as an isotropic hyperelastic-plastic material with a particular linear isotropic hardening (using scripted or XML commands):

Material "polymer","polymer","HEIsotropic"
   E 3100
   nu .4
   yield 72
   Ep 1000   
   alpha 60
   rho 1.2
 Done

 
 <Material Type="24" Name="Polymer">
   <rho>1.2</rho>
   <K>5166.67</K>
   <G1>1107.14</G1>
   <yield>72</yield>
   <Ep>1000</Ep>
   <alpha>60</alpha>
 </Material>

G1 represents G in the formulation on top

References

• Simo J. C. and T. J. R. Hughes (2000), "Computational Inelasticy", Interdisciplinary Applied Mechanics, Volume 7. Springer Edition.

• Simo J. C. (1988a), "Framework for finite elastoplasticity. Part I", Computer Methods in Applied Mechanics and Engineering, 66: 199-219.

• Simo J. C. (1988b), "Framework for finite elastoplasticity based on maximum dissipated energy and the multiplicative decomposition. Part II: Computational aspects", Computer Methods in Applied Mechanics and Engineering, 68: 1-31.