applications/solvers/solidMechanics/elasticNonLinULSolidFoam/createSolidInterface.H

   1 Switch solidInterfaceCorr(false);
   2
   3 solidInterface* 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 correction" << endl;
  14         solidInterfacePtr = new solidInterface(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 }