applications/solvers/solidMechanics/elasticPlasticSolidFoam/writeFields.H

   1 if (runTime.outputTime())
   2 {
   3     volScalarField epsilonEq
   4     (
   5         IOobject
   6         (
   7             "epsilonEq",
   8             runTime.timeName(),
   9             mesh,
  10             IOobject::NO_READ,
  11             IOobject::AUTO_WRITE
  12         ),
  13         sqrt((2.0/3.0)*magSqr(dev(epsilon)))
  14     );
  15
  16     Info<< "Max epsilonEq = " << max(epsilonEq).value() << endl;
  17
  18     volScalarField epsilonPEq
  19     (
  20         IOobject
  21         (
  22             "epsilonPEq",
  23             runTime.timeName(),
  24             mesh,
  25             IOobject::NO_READ,
  26             IOobject::AUTO_WRITE
  27         ),
  28         sqrt((2.0/3.0)*magSqr(dev(epsilonP)))
  29     );
  30
  31     Info<< "Max epsilonPEq = " << max(epsilonPEq).value() << endl;
  32
  33     volScalarField sigmaEq
  34     (
  35         IOobject
  36         (
  37             "sigmaEq",
  38             runTime.timeName(),
  39             mesh,
  40             IOobject::NO_READ,
  41             IOobject::AUTO_WRITE
  42         ),
  43         sqrt((3.0/2.0)*magSqr(dev(sigma)))
  44     );
  45
  46     Info<< "Max sigmaEq = " << max(sigmaEq).value() << endl;
  47
  48     volScalarField sigmaHyd
  49     (
  50         IOobject
  51         (
  52             "sigmaHyd",
  53             runTime.timeName(),
  54             mesh,
  55             IOobject::NO_READ,
  56             IOobject::AUTO_WRITE
  57         ),
  58         (
  59             sigma.component(symmTensor::XX)
  60           + sigma.component(symmTensor::YY)
  61           + sigma.component(symmTensor::ZZ)
  62         )/3
  63     );
  64     Info<< "Max sigmaHyd = " << max(sigmaHyd).value() << endl;
  65
  66     runTime.write();
  67 }