first push

NonLinearSolver/materialLaw/mlawHyperelastic.cpp

@@ -859,7 +859,7 @@ void mlawPowerYieldHyper::predictorCorrector_nonAssociatedFlow(const STensor3& F
       q1->getRefToDissipationActive() = true;
        // plasticity
       int ite = 0;
-      int maxite = 100; // maximal number of iters
+      int maxite = 1000; // maximal number of iters
 
 
       //Msg::Error("plasticity occurs f = %e",f);
@@ -919,9 +919,14 @@ void mlawPowerYieldHyper::predictorCorrector_nonAssociatedFlow(const STensor3& F
 
         if(ite > maxite){
           Msg::Error("No convergence for plastic correction in mlawPowerYieldHyper nonAssociatedFlow Maxwell iter = %d, f = %e!!",ite,f);
+	  if(fabs(f)<100.*_tol)
+	    break;
+	  else
+	  {
             P(0,0) = P(1,1) = P(2,2) = sqrt(-1.);
             return;
 	  }
+        }
       };
 
       q1->_DgammaDt = Dgamma/this->getTimeStep();