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