Isotropic Kelvin Material

As proposed by Kelvin [13] material strain $D\hackscore{ij}=v\hackscore{i,j}+v\hackscore{j,i}$ can be decomposed into an elastic part $D\hackscore{ij}^{el}$ and visco-plastic part $D\hackscore{ij}^{vp}$:

$$D\hackscore{ij}=D\hackscore{ij}^{el}+D\hackscore{ij}^{vp}$$ (136)

with the elastic strain given as

$$D\hackscore{ij}^{el'}=\frac{1}{2 \mu} \dot{\sigma}\hackscore{ij}'$$ (137)

where $\sigma'\hackscore{ij}$ is the deviatoric stress (Notice that $\sigma'\hackscore{ii}=0$). If the material is composed by materials $q$ the visco-plastic strain can be decomposed as

$$D\hackscore{ij}^{vp'}=\sum\hackscore{q} D\hackscore{ij}^{q'}$$ (138)

where $D\hackscore{ij}^{q}$ is the strain in material $q$ given as

$$D\hackscore{ij}^{q'}=\frac{1}{2 \eta^{q}} \sigma'\hackscore{ij}$$ (139)

where $\eta^{q}$ is the viscosity of material $q$. We assume the following betwee the the strain in material $q$

$$\eta^{q}=\eta^{q}\hackscore{N} \left(\frac{\tau}{\tau\hackscore{t}^q}\right)^{\frac{1}{n^{q}}-1} \tau=\sqrt{\frac{1}{2}\sigma'\hackscore{ij} \sigma'\hackscore{ij}}$$ (140)

for a given power law coefficients $n^{q}\ge1$ and transition stresses $\tau\hackscore{t}^q$, see [13]. Notice that $n^{q}=1$ gives a constant viscosity. After inserting equation 7.52 into equation 7.51 one gets:

$$D\hackscore{ij}'^{vp}=\frac{1}{2 \eta^{vp}} \sigma'\hackscore{ij} \frac{1}{\eta^{vp}} = \sum\hackscore{q} \frac{1}{\eta^{q}} \;.$$ (141)

and finally with 7.49

$$D\hackscore{ij}'=\frac{1}{2 \eta^{vp}} \sigma'\hackscore{ij}+\frac{1}{2 \mu} \dot{\sigma}\hackscore{ij}'$$ (142)

The total stress $\tau$ needs to fullfill the yield condition

$$\tau \le \tau\hackscore{Y} + \beta \; p$$ (143)

with the Drucker-Prager cohesion factor $\tau\hackscore{Y}$, Drucker-Prager friction $\beta$ and total pressure $p$. The deviatoric stress needs to fullfill the equilibrion equation

$$-\sigma'\hackscore{ij,j}+p\hackscore{,i}=F\hackscore{i}$$ (144)

where $F\hackscore{j}$ is a given external fource. We assume an incompressible media:

$$-v\hackscore{i,i}=0$$ (145)

Natural boundary conditions are taken in the form

$$\sigma'\hackscore{ij}n\hackscore{j}-n\hackscore{i}p=f$$ (146)

which can be overwritten by a constraint

$$v\hackscore{i}(x)=0$$ (147)

where the index $i$ may depend on the location $x$ on the bondary.