applications/solvers/stressAnalysis/solidEquilibriumDisplacementFoam/calculateStress.H

   1     if (runTime.outputTime())
   2     {
   3         volSymmTensorField sigma
   4         (
   5             IOobject
   6             (
   7                 "sigma",
   8                 runTime.timeName(),
   9                 mesh,
  10                 IOobject::NO_READ,
  11                 IOobject::AUTO_WRITE
  12             ),
  13             rho*sigmaD
  14         );
  15
  16         volScalarField sigmaEq
  17         (
  18             IOobject
  19             (
  20                 "sigmaEq",
  21                 runTime.timeName(),
  22                 mesh,
  23                 IOobject::NO_READ,
  24                 IOobject::AUTO_WRITE
  25             ),
  26             sqrt((3.0/2.0)*magSqr(dev(sigma)))
  27         );
  28
  29         Info<< "Max sigmaEq = " << max(sigmaEq).value()
  30             << endl;
  31
  32         runTime.write();
  33     }