Command manual

Command list


Material properties
"Optional title"
did, erode, noic
$\sigma_s$, $K_c$, $t_s$, $\alpha_s$, $\beta_s$

Parameter definition

did Unique damage identification number
erode Element erosion flag
0 $\rightarrow$ failed element is not eroded
1 $\rightarrow$ failed element is eroded
2 $\rightarrow$ node splitting at failure (crack plane orthogonal to max principal strain)
3 $\rightarrow$ node splitting at failure (crack plane orthogonal max principal stress)
noic Flag to turn off cracking along interface between different materials
0 $\rightarrow$ material interface cracks are allowed
1 $\rightarrow$ material interface cracks are not allowed
$\sigma_s$ Threshold stress (maximum principal stress) for initiation of fracture.
$K_c$ Stress intensity factor for crack propagation (only used with node splitting)
default: not used
$t_s$ Time to initiate fracture at threshold stress
default: not used
$\alpha_s$ Exponent controlling time to initiate fracture
default: not used
$\beta_s$ Parameter controlling the pressure dependency
default: not used


This is a brittle fracture criterion. The material cracks once the damage parameter, $D$, has evolved from 0 to 1. The damage is defined as:

$D = \displaystyle{ \frac{1}{t_s} \int_0^{t} H(\bar{\sigma}_1 - \sigma_s) \cdot (\bar{\sigma}_1 / \sigma_s )^{\alpha_s} } \mathrm{d}t$

$\bar{\sigma}_1$ is defined as:

$\bar{\sigma}_1 = \sigma_1^{dev} - (1-\beta_s) \cdot p $

where $\sigma_1^{dev}$ is the maximum deviatoric principal stress and $p$ is the pressure. With $\beta_s = 0$, $\bar{\sigma}_1$ equals the maximum principal stress. Note that the damage only grows if $\bar{\sigma}_1 \geq \sigma_s.$

Crack propagation is controlled by a stress intensity criterion (if having node splitting activated). The stress intensity $K_I$ is estimated for the integration points surrounding the crack tip. The crack will propagate if $K_I > K_c$ (Modus I crack).