applications/solvers/solidMechanics/elasticPlasticNonLinULSolidFoam/createSolidInterfaceNonLin.H

   1 Switch solidInterfaceCorr(false);
   2
   3 solidInterfaceNonLin* solidInterfacePtr(NULL);
   4
   5 {
   6     const dictionary& stressControl =
   7         mesh.solutionDict().subDict("solidMechanics");
   8
   9     solidInterfaceCorr = Switch(stressControl.lookup("solidInterface"));
  10
  11     if(solidInterfaceCorr)
  12     {
  13         Info << "Creating solid interface nonlinear correction" << endl;
  14         solidInterfacePtr = new solidInterfaceNonLin(mesh, rheology);
  15         solidInterfacePtr->modifyProperties(muf, lambdaf);
  16
  17         //- solidInterface needs muf and lambdaf to be used for divDSigmaExp
  18         if(divDSigmaExpMethod != "surface" && divDSigmaExpMethod != "decompose")
  19         {
  20             FatalError << "divDSigmaExp must be decompose or surface when solidInterface is on"
  21                  << exit(FatalError);
  22         }
  23         if(divDSigmaLargeStrainExpMethod != "surface")
  24         {
  25             FatalError << "divDSigmaLargeStrainExp must be surface when solidInterface is on"
  26                 << exit(FatalError);
  27         }
  28     }
  29 }