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