applications/solvers/solidMechanics/elasticPlasticNonLinTLSolidFoam/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(2nd Piola-Kirchoff)-strain(Green) 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 }