damage_elasticity#
Description#
This is a simple behavior for elastic energy based damage. It is
duplicated by the gen_evp **damage option with type
*elastic. This behavior is given as an example of a simple user
model programmed without ZebFront. The source is available in the
developer handbook.
(221)#\[\ten \sigma = (1-D)\tenf D_{el}:(\ten{\varepsilon}-\ten{\varepsilon}_{th})\]
(222)#\[\bar{Y} = \frac{1}{2}\ten \varepsilon_{el}:\tenf D_{el}: \ten \varepsilon_{el}\]
(223)#\[D = \alpha \left[ \mbox{max}\sqrt{\bar{Y}} - \sqrt{Y_0}\right]\]
Syntax#
***behavior damage_elasticity
\(~\,\) **elasticity <ELASTICITY>
\(~\,\) **Y0 COEFFICIENT
\(~\,\) **alpha COEFFICIENT
Stored Variables
The grad variable is the gradient of temperature, and the flux is the heat flux.
prefix |
size |
description |
default |
|---|---|---|---|
|
T-2 |
total (small deformation) strain |
yes |
|
T-2 |
Cauchy stress |
yes |
|
S |
\(\bar{Y}\) |
no |
|
S |
damage \(D\) |
yes |
Example#
***behavior damage_elasticity
**elasticity isotropic
young 200000.
poisson 0.0
**Y0 20.0
**alpha 0.001
***return