applications/solvers/solidMechanics/elasticNonLinIncrTLSolidFoam/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()
  17         << endl;
  18
  19     volScalarField sigmaEq
  20       (
  21        IOobject
  22        (
  23         "sigmaEq",
  24         runTime.timeName(),
  25         mesh,
  26         IOobject::NO_READ,
  27         IOobject::AUTO_WRITE
  28         ),
  29        sqrt((3.0/2.0)*magSqr(dev(sigma)))
  30        );
  31
  32     Info<< "Max sigmaEq = " << max(sigmaEq).value()
  33         << endl;
  34
  35     //- Calculate Cauchy stress
  36     volTensorField F = I + gradU;
  37     volScalarField J = det(F);
  38
  39     //- update density
  40     rho = rho/J;
  41
  42     volSymmTensorField sigmaCauchy
  43       (
  44        IOobject
  45        (
  46         "sigmaCauchy",
  47         runTime.timeName(),
  48         mesh,
  49         IOobject::NO_READ,
  50         IOobject::AUTO_WRITE
  51         ),
  52        (1/J) * symm(F.T() & sigma & F)
  53        );
  54
  55     //- Cauchy von Mises stress
  56     volScalarField sigmaCauchyEq
  57       (
  58        IOobject
  59        (
  60         "sigmaCauchyEq",
  61         runTime.timeName(),
  62         mesh,
  63         IOobject::NO_READ,
  64         IOobject::AUTO_WRITE
  65         ),
  66        sqrt((3.0/2.0)*magSqr(dev(sigmaCauchy)))
  67        );
  68
  69     Info<< "Max sigmaCauchyEq = " << max(sigmaCauchyEq).value()
  70         << endl;
  71
  72     volTensorField Finv = inv(F);
  73     volSymmTensorField epsilonAlmansi
  74       (
  75        IOobject
  76        (
  77         "epsilonAlmansi",
  78         runTime.timeName(),
  79         mesh,
  80         IOobject::NO_READ,
  81         IOobject::AUTO_WRITE
  82         ),
  83        symm(Finv & epsilon & Finv.T())
  84        );
  85
  86     //- boundary traction
  87     volVectorField traction
  88       (
  89        IOobject
  90        (
  91         "tractionCauchy",
  92         runTime.timeName(),
  93         mesh,
  94         IOobject::NO_READ,
  95         IOobject::AUTO_WRITE
  96         ),
  97        mesh,
  98        dimensionedVector("zero", dimForce/dimArea, vector::zero)
  99        );
 100     forAll(traction.boundaryField(), patchi)
 101       {
 102         tensorField Fpatch = I + gradU.boundaryField()[patchi];
 103
 104         traction.boundaryField()[patchi] =
 105           n.boundaryField()[patchi] & (sigma.boundaryField()[patchi] & Fpatch);
 106       }
 107
 108     runTime.write();
 109   }