Command manual

Command list


Material properties
Attention: This command is in the beta stage and the format may change over time.
"Optional title"
mid, $\rho$, $G$
$K_0$, $K_L$, $p_0$, $p_L$, $\epsilon_L$, $n$, $f_t$, $f_c$
$r_s$, $r_p$, $\epsilon_t$, $\epsilon_c$, $c$, $c_{dec}$, $\xi$, bulk

Parameter definition

mid Unique material identification number
$\rho$ Density
$G$ Shear modulus
$K_0$ Initial bulk modulus
$K_L$ Bulk modulus at full compaction
$p_0$ Hydrostatic pressure at crush limit
$p_L$ Hydrostatic pressure at full compaction
$\epsilon_L$ Volumetric strain at full compaction
$n$ Exponent controlling the shape of the compaction curve
default: 1
$f_t$ Uniaxial tensile strength
$f_c$ Uniaxial compressive strength
$r_s$ Strain rate shear hardening parameter
$r_p$ Strain rate pressure hardening parameter
$\epsilon_t$ Uni-axial tensile failure strain
$\epsilon_c$ Volumetric crushing failure strain
$c$ Viscosity parameter
$c_{dec}$ Viscous decay coefficient
default: not used
$\xi$ Dimensionless parameter controlling the shape of the yield surface cap
default: $\xi = f_c/2p_0$
bulk Flag to activate bulking
0 $\rightarrow$ no bulking
1 $\rightarrow$ bulking activated


This is a concrete model with different failure mechanisms in compression and tension. The material is assumed to have a pressure dependent shear resistance. Inelastic deformation is a combination of shearing and dilatation. Inelastic dilatation is interpreted as crushing that gradually reduces the shear resistance of the material. Deviatoric inelastic strains eventuelly lead to the formation of macroscopic cracks. Node splitting is used for the representation of such cracks. Note that node splitting can not occur in MERGE or REFINE interfaces.

The total stress $\sigma$ is the sum of an elastic component $\sigma^e$ and a viscous component $\sigma^v$.

$\sigma = \sigma^e + \sigma^v$


$\sigma^e = 2G \epsilon_{dev}^e - p \mathbf{I}$

$\epsilon_{dev}^e$ is the deviatoric elastic strain and $p$ is the pressure. The viscous stress component is defined as:

$\sigma^v = c \dot\epsilon$

where $\dot\epsilon$ is the total strain rate. Note that the flow criteria (in both tension and compression) are evaluated using the elastic stresses. In hydrostatic loading the elastic bulk modulus $K$ and the compaction pressure $p_c$ are assumed to grow with the ineastic compaction strain $\epsilon_v^p$:

$K = \left( 1 - \frac{\epsilon_v^p}{\epsilon_L} \right) K_0 + \frac{\epsilon_v^p}{\epsilon_L} K_L$
$\displaystyle{ p_c = p_0 + ( p_L - p_0 ) \left( \frac{\epsilon_v^p}{\epsilon_L} \right)^n + r_p \dot\epsilon_v^p }$

The volumetric compaction strain $\epsilon_v^p$ can not exceed $\epsilon_L$. The deviatoric yield stress $\sigma_y$ is defined as (se figure below):

$\sigma_y = \left\{ \begin{array}{ccc} 0 & : & p \leq p_s \\ \eta (p-p_s) & : & p_s \lt p \lt \xi p_c \\ \eta (p_c-p_s) \displaystyle{\sqrt{1 - \left( \frac{1 - p/p_c}{1 - \xi} \right)^2}} & : & \xi p_c \leq p \leq p_c \end{array} \right. $

where $\eta$ and $p_s$ (hydrostatic pressure cut-off) are derived internally to correctly match the given (quasistatic) uniaxial tensile and the compressive strengths $(f_t, f_c)$.

$\eta = \displaystyle{\frac{3}{2} \cdot \frac{f_c - 2f_t}{f_t + f_c} \cdot \frac{1}{(1-D_c) \cdot \mathrm{cos}(\theta)+D_c}}$
$p_s = -\displaystyle{\left( \frac{f_t f_c}{f_c - 2f_t} \cdot (1-D_t) + \frac{r_s \dot\epsilon_{eff}^p}{\eta} \right) \cdot (1 - D_c) \cdot (1 - D_0)}$

$\theta$ is the Lode angle, as define in the figure below. $D_t$ is the tensile damage, $D_c$ is the crushing damage and $D_0$ is the initial defect level. Note that $p_s=0$ for a fully damaged material. Also note the shear rate hardening, controlled by the parameters $r_s$ and $\dot\epsilon_{s0}$. Initial defects are optional and can be defined with the command INITIAL_DAMAGE_RANDOM.

The volumetric crushing damage $D_c$ is defined as:

$\displaystyle{D_c = \mathrm{min}(1, \frac{\epsilon_v^p}{\epsilon_c})}$

Tensile/shear cracking is controlled by the damage variable $D_t$. Node splitting is activated when $D_t$ reaches 1.

$\displaystyle{D_t = \int_0^{\epsilon_{eff}^p} \frac{1}{\epsilon_f}\mathrm{d}\epsilon_{eff}^p}$

$\epsilon_f$ is a pressure dependent failure strain:

$\displaystyle{\epsilon_f = \epsilon_t + \frac{\mathrm{max}(0, p)}{p_c}(\epsilon_c-\epsilon_t)}$

If $G_I$ is defined, a post fracture unloading is applied. The unloading (crack softening) process dissipates the energy $G_I$ per unit area of cracking.

Pressure and Lode angle dependent yield surface with cap
Pressure and Lode angle dependent yield surface with cap