Command manual

Command list


Material properties
"Optional title"
mid, $\rho$, $E$, $\nu$
cid${}_1$, cid${}_2$, $\xi$, $\eta$, $\sigma_{dev}^{max}$, $B_0$, $B_1$

Parameter definition

mid Unique material identification number
$\rho$ Density
$E$ Young's modulus
$\nu$ Poisson's ratio
cid${}_1$ ID of a CURVE or FUNCTION defining pressure versus inelastic volumetric compaction
cid${}_2$ ID of a CURVE or FUNCTION defining deviatoric yield stress due to grain adhesion
default: no adhesion
$\xi$ Parameter controlling the shape of the yield surface
$\eta$ Parameter controlling the shape of the yield surface
$\sigma_{dev}^{max}$ Upper cap to deviatoric yield strength
default: no upper cap
$B_0$ Failure parameter
default: no failure
$B_1$ Failure parameter
default: 0


This is a material model for granular media, where inelastic deviatoric flow and volumetric flow are coupled. The shape of the yield surface is shown in the figure below. Note that the yield surface is split into two regions (I) and (II). Volumetric flow (compaction) is assumed to only occur in region II (associated flow). Plastic flow in region I is purely deviatoric (J2). In region II the surface is elliptic:

$\Phi = \displaystyle{\sqrt{ \left( \frac{\sigma_{dev}^{eff}}{\xi p_c}\right)^2 + \left( \frac{p - \xi p_c}{p_c-\xi p_c}\right)^2 }} - 1$

cid1 describes the compaction pressure as function of the effective volumetric strain:

$p_c = \hat p_c(\epsilon_{vol}^{eff})$.

cid2 describes the grain adhesion due to compaction (see yield surface in figure):

$\sigma_a = \hat\sigma_a(\epsilon_{vol}^{eff})$.

Damage $D$ is assumed to grow at inelastic deviatoric deformations in region I:

$\displaystyle{ \dot D = \frac{ \dot\epsilon_{dev}^{eff}}{B_0} \left( 1 - \frac{p}{p_c} \right)^{B_1}}$

Damage falls back to zero whenever the material undergoes further compaction (region II). As $D$ reaches 1 the adhesive components of the yield stress are set to zero and pressure is not allowed to drop below 0.

$D$ is set to 1 (instant failure) if meeting the pressure cut-off criterion, i.e. if $p = -\sigma_a$.

Yield surface and flow law
Yield surface and flow law