applications/solvers/solidMechanics/elasticPlasticNonLinTLSolidFoam/calculateDivDSigmaNonLinExp.H

   1 if(divDSigmaNonLinExpMethod == "standard")
   2   {
   3     divDSigmaNonLinExp = fvc::div
   4       (
   5        ( mu * (
   6                (gradDU & gradU.T())
   7                + (gradU & gradDU.T())
   8                + (gradDU & gradDU.T())
   9                ) )
  10        + ( lambda * 0.5 * tr( (gradDU & gradU.T())
  11                               + (gradU & gradDU.T())
  12                               + (gradDU & gradDU.T())
  13                               ) * I )
  14        + ( DSigma & gradU )
  15        + ( (sigma + DSigma) & gradDU ),
  16        "div(sigma)"
  17        );
  18   }
  19  else if(divDSigmaNonLinExpMethod == "surface")
  20    {
  21      divDSigmaNonLinExp =
  22        fvc::div(
  23                 mesh.magSf()
  24                 *(
  25                   ( muf * (n & fvc::interpolate( (gradU & gradDU.T()) + (gradDU & gradU.T()) + (gradDU & gradDU.T()) )) )
  26                   + ( 0.5*lambdaf * (n * tr(fvc::interpolate( (gradU & gradDU.T()) + (gradDU & gradU.T()) + (gradDU & gradDU.T()) ))) )
  27                   + (n & fvc::interpolate( (DSigma & gradU) + ((sigma + DSigma) & gradDU) ))
  28                   )
  29                 );
  30
  31      // divDSigmaNonLinExp =
  32      //   fvc::div
  33      //   (
  34      //         ( muf * mesh.Sf()
  35      //           & fvc::interpolate(
  36      //                              (gradDU & gradU.T())
  37      //                              + (gradU & gradDU.T())
  38      //                              + (gradDU & gradDU.T())
  39      //                              ) )
  40      //         + ( lambdaf * 0.5* tr(
  41      //                               fvc::interpolate(
  42      //                                                (gradDU & gradU.T())
  43      //                                                + (gradU & gradDU.T())
  44      //                                                + (gradDU & gradDU.T())
  45      //                                                )
  46      //                               ) * mesh.Sf() )
  47      //         + (mesh.Sf() & fvc::interpolate( (DSigma & gradU )
  48      //                                          + ( (sigma + DSigma) & gradDU )
  49      //                                          ) )
  50      //         );
  51    }
  52  else
  53    {
  54      FatalError << "divDSigmaExp method " << divDSigmaExpMethod << " not found!" << endl;
  55    }