"Optional title"
did, erode, noic
$\sigma_s$, $K_c$, $t_s$, $\alpha_s$, $\beta_s$
Parameter definition
Variable | Description |
---|---|
did | Unique damage identification number |
erode | Element erosion flag |
noic | Flag to turn off cracking along interface between different materials |
$\sigma_s$ | Threshold stress (maximum principal stress) for initiation of fracture. |
$K_c$ | Stress intensity factor for crack propagation (only used with node splitting) |
$t_s$ | Time to initiate fracture at threshold stress |
$\alpha_s$ | Exponent controlling time to initiate fracture |
$\beta_s$ | Parameter controlling the pressure dependency |
Description
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).