applications/solvers/solidMechanics/elasticPlasticNonLinULSolidFoam/calculateDivDSigmaNonLinExp.H

   1 if(divDSigmaNonLinExpMethod == "standard")
   2 {
   3     divDSigmaNonLinExp = fvc::div
   4     (
   5         ( mu * (gradDU & gradDU.T()) )
   6       + ( 0.5 * lambda * tr(gradDU & gradDU.T()) * I )
   7       + ( (sigma + DSigma) & gradDU ),
   8         "div(sigma)"
   9     );
  10 }
  11 else if(divDSigmaNonLinExpMethod == "surface")
  12 {
  13     divDSigmaNonLinExp = fvc::div
  14     (
  15         mesh.magSf()*
  16         (
  17             ( muf * (n & fvc::interpolate( gradDU & gradDU.T() )) )
  18           + ( 0.5*lambdaf * (n * tr(fvc::interpolate( gradDU & gradDU.T() ))) )
  19           + (n & fvc::interpolate( (sigma + DSigma) & gradDU ))
  20         )
  21     );
  22 }
  23 else
  24 {
  25     FatalError << "divDSigmaExp method " << divDSigmaExpMethod << " not found!" << endl;
  26 }