applications/solvers/stressAnalysis/solidEquilibriumDisplacementFoam/createFields.H

   1 Info<< "Reading displacement field D\n" << endl;
   2 volVectorField D
   3 (
   4     IOobject
   5     (
   6         "D",
   7         runTime.timeName(),
   8         mesh,
   9         IOobject::MUST_READ,
  10         IOobject::AUTO_WRITE
  11     ),
  12     mesh
  13 );
  14
  15 Info<< "Reading/setting displacement correction field Dcorr\n" << endl;
  16 volVectorField Dcorr
  17 (
  18     IOobject
  19     (
  20         "Dcorr",
  21         runTime.timeName(),
  22         mesh
  23     ),
  24     D
  25 );
  26 Dcorr *= 0.0;
  27
  28 Info<< "Calculating stress field sigmaD\n" << endl;
  29 volSymmTensorField sigmaD
  30 (
  31     IOobject
  32     (
  33         "sigmaD",
  34         runTime.timeName(),
  35         mesh,
  36         IOobject::NO_READ,
  37         IOobject::NO_WRITE
  38     ),
  39     mu*twoSymm(fvc::grad(D)) + (lambda*I)*tr(fvc::grad(D))
  40 );
  41
  42 volTensorField sigmaExp
  43 (
  44     IOobject
  45     (
  46         "sigmaExp",
  47         runTime.timeName(),
  48         mesh,
  49         IOobject::NO_READ,
  50         IOobject::NO_WRITE
  51     ),
  52     (lambda - mu)*fvc::grad(Dcorr) + mu*fvc::grad(Dcorr)().T()
  53   + (lambda*I)*tr(fvc::grad(Dcorr))
  54 );