Stress History Based Degradation With Custom Rules
This model implements the universal damage model assuming the degradation factor depends on the stress history.
material CustomStressDegradation (1) (2) (3) [4]
# (1) int, unique material tag
# (2) int, host intact material tag
# (3) int, tension degradation expression tag
# [4] int, compression degradation expression tag, default: (3)
It is possible to assign different degradation evolution rules for tension and compression.
The degradation expression shall take one scalar input and return a vector of two scalars, of which the first is the degradation factor while the second is the derivative of the degradation factor with respect to the input.
The degradations for positive/negative stresses are defined separately.
$$ D_{pos}=D_{pos}(\bar{\sigma}{max}),\quad D{neg}=D_{neg}(\bar{\sigma}_{min}). $$
In which $$\bar{\sigma}{max}$$ is the maximum stress (positive) of the whole loading history, and $$\bar{\sigma}{min}$$ is the minimum stress (negative) of the whole loading history.
To use this model, expressions of degradation evolution rules must be defined first. Assuming the degradation is an exponential function, and will be activated when the effective stress is larger than 20,
the corresponding derivative is then
One can then define the expression in a plain text file as follows.
// file: exp_t
if(x>20){
y[0]:=exp((-x+20)/20);
y[1]:=-y[0]/20;
}else{
y[0]:=1.;
y[1]:=0.;
}
To use it, one shall use the SimpleVector
expression.
expression SimpleVector 1 x y|2 exp_t
Here we use
We use an elastic host material for demonstration.
material Elastic1D 1 1000
expression SimpleVector 1 x y|2 exp_t
material CustomStressDegradation 2 1 1
materialTest1D 2 1E-3 10 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 105 110 115 120 60
We can define a different rule for compression.
// file: exp_c
if(x<-10){
y[0]:=exp((x+10)/30);
y[1]:=y[0]/30;
}else{
y[0]:=1.;
y[1]:=0.;
}
material Elastic1D 1 1000
expression SimpleVector 1 x y|2 exp_t
expression SimpleVector 2 x y|2 exp_c
material CustomStressDegradation 2 1 1 2
materialTest1D 2 1E-3 10 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 105 110 115 120 60