applications/solvers/solidMechanics/elasticPlasticSolidFoam/createHistory.H

   1 OFstream * forceFilePtr(NULL);
   2 OFstream * stressFilePtr(NULL);
   3 word historyPatchName(mesh.solutionDict().subDict("solidMechanics").lookup("historyPatch"));
   4 label historyPatchID = mesh.boundaryMesh().findPatchID(historyPatchName);
   5 if(historyPatchID == -1)
   6   {
   7     Warning << "history patch " << historyPatchName
   8             << " not found. Force-displacement will not be written"
   9             << endl;
  10   }
  11  else if(Pstream::master())
  12   {
  13     fileName historyDir = "history";
  14     mkDir(historyDir);
  15
  16     {
  17       fileName forceFileName(historyDir/"forceDisp_"+historyPatchName+".dat");
  18       Info << "\nForce-displacement for patch " << historyPatchName
  19            << " will be written to " << forceFileName
  20            << endl;
  21       forceFilePtr = new OFstream(forceFileName);
  22       OFstream& forceDispFile = *forceFilePtr;
  23       forceDispFile << "#Disp(mm)\tForce(N)" << endl;
  24     }
  25
  26     {
  27       fileName stressFileName(historyDir/"stressStrain_"+historyPatchName+".dat");
  28       Info << "\nStress(Engineering Small Stress)-strain(Engineering Small Strain) for patch "
  29            << historyPatchName
  30            << " will be written to " << stressFileName
  31            << endl;
  32       stressFilePtr = new OFstream(stressFileName);
  33       OFstream& stressStrainFile = *stressFilePtr;
  34       stressStrainFile << "#Strain(-)\tStress(Pa)" << endl;
  35     }
  36   }