1@Parser Implicit;
2@Behaviour Norton;
3@Algorithm NewtonRaphson_NumericalJacobian;
4@RequireStiffnessTensor;
5@MaterialProperty real A;
6@MaterialProperty real m;
7@StateVariable real    p;
8@ComputeStress{  sig = D*eel;}
9@Integrator{
10  real seq = sigmaeq(sig);
11  Stensor n = Stensor(0.);
12  if(seq > 1.e-15){
13    n   = 1.5*deviator(sig)/seq;
14  }
15  feel = deel + dp*n-deto;
16  fp   = dp - dt*A*pow(seq,m);
17}
18@TangentOperator{
19  Stensor4 Je;
20  getPartialJacobianInvert(Je);
21  Dt = D*Je;
22}
23
24