applications/solvers/solidMechanics/deprecatedSolvers/newContactStressFoam/createFields.H

   1     mesh.weights();
   2
   3     Info<< "Reading field U\n" << endl;
   4     volVectorField U
   5     (
   6         IOobject
   7         (
   8             "U",
   9             runTime.timeName(),
  10             mesh,
  11             IOobject::MUST_READ,
  12             IOobject::AUTO_WRITE
  13         ),
  14         mesh
  15     );
  16
  17     // Create gradient field
  18     volTensorField gradU = fvc::grad(U);
  19
  20     volSymmTensorField sigma
  21     (
  22         IOobject
  23         (
  24             "sigma",
  25             runTime.timeName(),
  26             mesh,
  27             IOobject::READ_IF_PRESENT,
  28             IOobject::AUTO_WRITE
  29         ),
  30         mesh,
  31         dimensionedSymmTensor("zero", dimForce/dimArea, symmTensor::zero)
  32     );
  33
  34     // Create rheology model
  35     rheologyModel rheology(sigma);
  36
  37     // Create contact problem
  38     contactProblem contact(U, gradU);