applications/solvers/solidMechanics/elasticNonLinULSolidFoam/calculateDivDSigmaLargeStrainExp.H

   1 //----------------------------------------------------//
   2 //- sigma explicit non linear 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       + 0.5*lambda*tr(gradDU & gradDU.T())*I
  11       + ((sigma + DSigma) & gradDU),
  12         "div(sigma)"
  13     );
  14 }
  15 else if(divDSigmaLargeStrainExpMethod == "surface")
  16 {
  17     divDSigmaLargeStrainExp = fvc::div
  18     (
  19         muf * (mesh.Sf() & fvc::interpolate(gradDU & gradDU.T()))
  20       + 0.5*lambdaf * (mesh.Sf() & (fvc::interpolate(gradDU && gradDU)*I))
  21       + (mesh.Sf() & fvc::interpolate( (sigma + DSigma) & gradDU ))
  22     );
  23 }
  24 else
  25 {
  26     FatalError
  27         << "divDSigmaLargeStrainMethod not found!"
  28         << exit(FatalError);
  29 }