applications/solvers/solidMechanics/elasticPlasticSolidFoam/calculatePlasticResidual.H

   1 {
   2   scalar magDEp = gMax(mag(DEpsilonP.internalField()));
   3
   4   if (magDEp < SMALL)
   5     {
   6       magDEp = SMALL;
   7     }
   8
   9   plasticResidual =
  10     gMax
  11     (
  12      mag
  13      (
  14       DEpsilonP.internalField()
  15       - DEpsilonP.prevIter().internalField()
  16       )
  17      /magDEp
  18      );
  19 }