applications/solvers/solidMechanics/elasticIncrSolidFoam/createFields.H

   1     Info<< "Reading field DU\n" << endl;
   2     volVectorField DU
   3     (
   4         IOobject
   5         (
   6             "DU",
   7             runTime.timeName(),
   8             mesh,
   9             IOobject::MUST_READ,
  10             IOobject::AUTO_WRITE
  11         ),
  12         mesh
  13     );
  14     volTensorField gradDU = fvc::grad(DU);
  15
  16     volVectorField U
  17     (
  18         IOobject
  19         (
  20             "U",
  21             runTime.timeName(),
  22             mesh,
  23             IOobject::READ_IF_PRESENT,
  24             IOobject::AUTO_WRITE
  25         ),
  26         mesh,
  27         dimensionedVector("zero", dimLength, vector::zero)
  28     );
  29
  30     volSymmTensorField DEpsilon
  31     (
  32         IOobject
  33         (
  34             "DEpsilon",
  35             runTime.timeName(),
  36             mesh,
  37             IOobject::READ_IF_PRESENT,
  38             IOobject::AUTO_WRITE
  39         ),
  40         mesh,
  41         dimensionedSymmTensor("zero", dimless, symmTensor::zero)
  42     );
  43
  44     volSymmTensorField epsilon
  45     (
  46         IOobject
  47         (
  48             "epsilon",
  49             runTime.timeName(),
  50             mesh,
  51             IOobject::READ_IF_PRESENT,
  52             IOobject::AUTO_WRITE
  53         ),
  54         mesh,
  55         dimensionedSymmTensor("zero", dimless, symmTensor::zero)
  56     );
  57
  58     volSymmTensorField DSigma
  59     (
  60         IOobject
  61         (
  62             "DSigma",
  63             runTime.timeName(),
  64             mesh,
  65             IOobject::READ_IF_PRESENT,
  66             IOobject::AUTO_WRITE
  67         ),
  68         mesh,
  69         dimensionedSymmTensor("zero", dimForce/dimArea, symmTensor::zero)
  70     );
  71
  72     volSymmTensorField sigma
  73     (
  74         IOobject
  75         (
  76             "sigma",
  77             runTime.timeName(),
  78             mesh,
  79             IOobject::READ_IF_PRESENT,
  80             IOobject::AUTO_WRITE
  81         ),
  82         mesh,
  83         dimensionedSymmTensor("zero", dimForce/dimArea, symmTensor::zero)
  84     );
  85
  86     volVectorField divDSigmaExp
  87     (
  88         IOobject
  89         (
  90             "divDSigmaExp",
  91             runTime.timeName(),
  92             mesh,
  93             IOobject::NO_READ,
  94             IOobject::NO_WRITE
  95         ),
  96         mesh,
  97         dimensionedVector("zero", dimensionSet(1,-2,-2,0,0,0,0), vector::zero)
  98     );
  99
 100     constitutiveModel rheology(sigma, DU);
 101
 102     volScalarField rho = rheology.rho();
 103
 104     volScalarField mu = rheology.mu();
 105     volScalarField lambda = rheology.lambda();
 106     surfaceScalarField muf = fvc::interpolate(mu, "mu");
 107     surfaceScalarField lambdaf = fvc::interpolate(lambda, "lambda");
 108
 109     surfaceVectorField n = mesh.Sf()/mesh.magSf();