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