applications/solvers/solidMechanics/elasticPlasticNonLinULSolidFoam/calculateDivDSigmaLargeStrainExp.H

   1 //----------------------------------------------------//
   2 //- sigma explicit large strain explicit terms
   3 //----------------------------------------------------//
   4 if(divDSigmaLargeStrainExpMethod == "standard")
   5 {
   6     divDSigmaLargeStrainExp = fvc::div
   7     (
   8         mu*(gradDU & gradDU.T())
   9       + 0.5*lambda*(gradDU && gradDU)*I //- equivalent to 0.5*lambda*(I*tr(gradDU & gradDU.T()))
  10       + ((sigma + DSigma) & DF.T()),
  11         "div(sigma)"
  12     );
  13 }
  14 else if(divDSigmaLargeStrainExpMethod == "surface")
  15 {
  16     divDSigmaLargeStrainExp = fvc::div
  17     (
  18         muf * (mesh.Sf() & fvc::interpolate(gradDU & gradDU.T()))
  19       + 0.5*lambdaf * (mesh.Sf() & (fvc::interpolate(gradDU && gradDU)*I))
  20       + (mesh.Sf() & fvc::interpolate( sigma & DF.T() ))
  21       + (mesh.Sf() & fvc::interpolate(DSigma & DF.T() ))
  22     );
  23 }
  24 else
  25 {
  26     FatalError
  27         << "divDSigmaLargeStrainMethod not found!"
  28         << exit(FatalError);
  29 }
  30
  31 //- relax
  32 divDSigmaLargeStrainExp.relax();