applications/solvers/solidMechanics/elasticNonLinULSolidFoam/calculateDivDSigmaLargeStrainExp.H

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