applications/solvers/stressAnalysis/solidDisplacementFoam/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         if (thermalStress)
  17         {
  18             const volScalarField& T = Tptr();
  19             sigma = sigma - I*(rho*threeKalpha*T);
  20         }
  21
  22         volScalarField sigmaEq
  23         (
  24             IOobject
  25             (
  26                 "sigmaEq",
  27                 runTime.timeName(),
  28                 mesh,
  29                 IOobject::NO_READ,
  30                 IOobject::AUTO_WRITE
  31             ),
  32             sqrt((3.0/2.0)*magSqr(dev(sigma)))
  33         );
  34
  35         Info<< "Max sigmaEq = " << max(sigmaEq).value()
  36             << endl;
  37
  38         runTime.write();
  39     }