# Nonlinear finite elements/Homework 11/Solutions/Problem 1/Part 15

## Problem 1: Part 15: Finding the plastic flow parameter

The discretized form of the Kuhn-Tucker conditions in conjunction with the consistency condition gives us

${\displaystyle f({\boldsymbol {\sigma }}_{n+1},\alpha _{n+1},T_{n+1})=0~.}$

Use this condition and the relations you have derived in the previous sections to arrive at a nonlinear equation in ${\displaystyle \Delta \gamma }$  that can be solved using Newton iterations.

The yield function is

${\displaystyle f={\sqrt {\cfrac {3}{2}}}~{\sqrt {\mathbf {s} :\mathbf {s} }}-\left[\sigma _{0}+B\alpha ^{n}\right]\left[1-\left({\cfrac {T-T_{0}}{T_{m}-T_{0}}}\right)\right]}$

Therefore the discretized form of Kuhn-Tucker + consistency is

${\displaystyle {\sqrt {\cfrac {3}{2}}}~{\sqrt {\mathbf {s} _{n+1}:\mathbf {s} _{n+1}}}-\left[\sigma _{0}+B\alpha _{n+1}^{n}\right]\left[1-\left({\cfrac {T_{n+1}-T_{0}}{T_{m}-T_{0}}}\right)\right]=0}$

or,

${\displaystyle \mathbf {s} _{n+1}:\mathbf {s} _{n+1}={\cfrac {2}{3}}\left[\sigma _{0}+B\alpha _{n+1}^{n}\right]^{2}\left[1-\left({\cfrac {T_{n+1}-T_{0}}{T_{m}-T_{0}}}\right)\right]^{2}}$

Now,

${\displaystyle \mathbf {s} _{n+1}=\mathbf {s} _{n+1}^{\text{trial}}-2~\mu ~{\sqrt {\cfrac {3}{2}}}~\Delta \gamma ~\mathbf {n} _{n}}$

Therefore,

{\displaystyle {\begin{aligned}\mathbf {s} _{n+1}:\mathbf {s} _{n+1}&=\mathbf {s} _{n+1}^{\text{trial}}:\mathbf {s} _{n+1}^{\text{trial}}-4~\mu ~{\sqrt {\cfrac {3}{2}}}~\Delta \gamma ~\mathbf {s} _{n+1}^{\text{trial}}:\mathbf {n} _{n}+4~\mu ^{2}~{\cfrac {3}{2}}~(\Delta \gamma )^{2}~\mathbf {n} _{n}:\mathbf {n} _{n}\\&=\mathbf {s} _{n+1}^{\text{trial}}:\mathbf {s} _{n+1}^{\text{trial}}-4~\mu ~{\sqrt {\cfrac {3}{2}}}~\Delta \gamma ~\mathbf {s} _{n+1}^{\text{trial}}:\mathbf {n} _{n}+6~\mu ^{2}~(\Delta \gamma )^{2}\end{aligned}}}

Plugging into the discretized yield condition, we have

${\displaystyle {\cfrac {2}{3}}\left[\sigma _{0}+B\alpha _{n+1}^{n}\right]^{2}\left[1-\left({\cfrac {T_{n+1}-T_{0}}{T_{m}-T_{0}}}\right)\right]^{2}=\mathbf {s} _{n+1}^{\text{trial}}:\mathbf {s} _{n+1}^{\text{trial}}-4~\mu ~{\sqrt {\cfrac {3}{2}}}~\Delta \gamma ~\mathbf {s} _{n+1}^{\text{trial}}:\mathbf {n} _{n}+6~\mu ^{2}~(\Delta \gamma )^{2}}$

or,

${\displaystyle \left[\sigma _{0}+B\alpha _{n+1}^{n}\right]^{2}\left[1-\left({\cfrac {T_{n+1}-T_{0}}{T_{m}-T_{0}}}\right)\right]^{2}={\cfrac {3}{2}}~\mathbf {s} _{n+1}^{\text{trial}}:\mathbf {s} _{n+1}^{\text{trial}}-6~\mu ~{\sqrt {\cfrac {3}{2}}}~\Delta \gamma ~\mathbf {s} _{n+1}^{\text{trial}}:\mathbf {n} _{n}+9~\mu ^{2}~(\Delta \gamma )^{2}}$

Also

{\displaystyle {\begin{aligned}\alpha _{n+1}&=\alpha _{n}+\Delta \gamma ~{\cfrac {{\boldsymbol {\varepsilon }}_{n}^{p}:\mathbf {n} _{n}}{\lVert {\boldsymbol {\varepsilon }}_{n}^{p}\rVert _{}}}\\T_{n+1}&=T_{n}+{\sqrt {\cfrac {3}{2}}}~{\cfrac {\chi ~\Delta \gamma }{\rho _{n}~C_{p}}}~\lVert \mathbf {s} _{n}\rVert _{}\end{aligned}}}

Therefore,

${\displaystyle \left[\sigma _{0}+B\left\{\alpha _{n}+\Delta \gamma ~{\cfrac {{\boldsymbol {\varepsilon }}_{n}^{p}:\mathbf {n} _{n}}{\lVert {\boldsymbol {\varepsilon }}_{n}^{p}\rVert _{}}}\right\}^{n}\right]^{2}\left[1-\left({\cfrac {T_{n}+{\sqrt {\cfrac {3}{2}}}~{\cfrac {\chi ~\Delta \gamma }{\rho _{n}~C_{p}}}~\lVert \mathbf {s} _{n}\rVert _{}-T_{0}}{T_{m}-T_{0}}}\right)\right]^{2}={\cfrac {3}{2}}~\mathbf {s} _{n+1}^{\text{trial}}:\mathbf {s} _{n+1}^{\text{trial}}-6~\mu ~{\sqrt {\cfrac {3}{2}}}~\Delta \gamma ~\mathbf {s} _{n+1}^{\text{trial}}:\mathbf {n} _{n}+9~\mu ^{2}~(\Delta \gamma )^{2}}$

The nonlinear equation in ${\displaystyle \Delta \gamma }$  is

{\displaystyle {\begin{aligned}g(\Delta \gamma )&=0\\&=9~\mu ^{2}~(\Delta \gamma )^{2}-6~\mu ~{\sqrt {\cfrac {3}{2}}}~\Delta \gamma ~\mathbf {s} _{n+1}^{\text{trial}}:\mathbf {n} _{n}\\&\qquad -\left[\sigma _{0}+B\left\{\alpha _{n}+\Delta \gamma ~{\cfrac {{\boldsymbol {\varepsilon }}_{n}^{p}:\mathbf {n} _{n}}{\lVert {\boldsymbol {\varepsilon }}_{n}^{p}\rVert _{}}}\right\}^{n}\right]^{2}\left[1-\left({\cfrac {T_{n}+{\sqrt {\cfrac {3}{2}}}~{\cfrac {\chi ~\Delta \gamma }{\rho _{n}~C_{p}}}\lVert \mathbf {s} _{n}\rVert _{}-T_{0}}{T_{m}-T_{0}}}\right)\right]^{2}\\&\qquad +{\cfrac {3}{2}}~\mathbf {s} _{n+1}^{\text{trial}}:\mathbf {s} _{n+1}^{\text{trial}}\end{aligned}}}