applications/solvers/solidMechanics/elasticOrthoSolidFoam/createFields.H

   1     Info<< "Reading field U\n" << endl;
   2     volVectorField U
   3     (
   4         IOobject
   5         (
   6             "U",
   7             runTime.timeName(),
   8             mesh,
   9             IOobject::MUST_READ,
  10             IOobject::AUTO_WRITE
  11         ),
  12         mesh
  13     );
  14
  15     volTensorField gradU //= fvc::grad(U);
  16     (
  17         IOobject
  18         (
  19             "grad(U)",
  20             runTime.timeName(),
  21             mesh,
  22             IOobject::NO_READ,
  23             IOobject::NO_WRITE
  24         ),
  25         mesh,
  26         dimensionedTensor("zero", dimless, tensor::zero)
  27     );
  28
  29     volSymmTensorField epsilon
  30     (
  31         IOobject
  32         (
  33             "epsilon",
  34             runTime.timeName(),
  35             mesh,
  36             IOobject::READ_IF_PRESENT,
  37             IOobject::AUTO_WRITE
  38         ),
  39         mesh,
  40         dimensionedSymmTensor("zero", dimless, symmTensor::zero)
  41     );
  42
  43     volSymmTensorField sigma
  44     (
  45         IOobject
  46         (
  47             "sigma",
  48             runTime.timeName(),
  49             mesh,
  50             IOobject::READ_IF_PRESENT,
  51             IOobject::AUTO_WRITE
  52         ),
  53         mesh,
  54         dimensionedSymmTensor("zero", dimForce/dimArea, symmTensor::zero)
  55     );
  56
  57     volVectorField divSigmaExp
  58     (
  59         IOobject
  60         (
  61             "divSigmaExp",
  62             runTime.timeName(),
  63             mesh,
  64             IOobject::NO_READ,
  65             IOobject::NO_WRITE
  66         ),
  67         mesh,
  68         dimensionedVector("zero", dimForce/dimVolume, vector::zero)
  69     );
  70
  71
  72 //- rheology
  73 constitutiveModel rheology(sigma, U);
  74 volSymmTensor4thOrderField C = rheology.C();
  75 surfaceSymmTensor4thOrderField Cf = fvc::interpolate(C, "C");
  76 volDiagTensorField K = rheology.K();
  77 surfaceDiagTensorField Kf = fvc::interpolate(K, "K");
  78 volScalarField rho = rheology.rho();
  79
  80 surfaceVectorField n = mesh.Sf()/mesh.magSf();