applications/solvers/solidMechanics/icoFsiElasticNonLinULSolidFoam/createStressFields.H

   1     Info<< "Reading incremental displacement field DU\n" << endl;
   2
   3     volVectorField DU
   4     (
   5         IOobject
   6         (
   7             "DU",
   8             runTime.timeName(),
   9             stressMesh,
  10             IOobject::MUST_READ,
  11             IOobject::AUTO_WRITE
  12         ),
  13         stressMesh
  14     );
  15
  16     volTensorField gradDU = fvc::grad(DU);
  17
  18     volVectorField Usolid
  19     (
  20         IOobject
  21         (
  22          "Usolid",
  23          runTime.timeName(),
  24          stressMesh,
  25          IOobject::READ_IF_PRESENT,
  26          IOobject::AUTO_WRITE
  27          ),
  28         DU
  29      );
  30
  31
  32     Info<< "Reading incremental displacement field DV\n" << endl;
  33     volVectorField DV
  34     (
  35         IOobject
  36         (
  37             "DV",
  38             runTime.timeName(),
  39             stressMesh,
  40             IOobject::READ_IF_PRESENT,
  41             IOobject::AUTO_WRITE
  42         ),
  43         fvc::ddt(DU)
  44     );
  45
  46
  47     Info<< "Reading accumulated velocity field V\n" << endl;
  48     volVectorField Vs
  49     (
  50         IOobject
  51         (
  52             "Vs",
  53             runTime.timeName(),
  54             stressMesh,
  55             IOobject::READ_IF_PRESENT,
  56             IOobject::AUTO_WRITE
  57         ),
  58         stressMesh,
  59         dimensionedVector("zero", dimVelocity, vector::zero)
  60     );
  61
  62
  63     Info << "Reading accumulated stress field sigma\n" << endl;
  64     volSymmTensorField sigma
  65     (
  66         IOobject
  67         (
  68             "sigma",
  69             runTime.timeName(),
  70             stressMesh,
  71             IOobject::READ_IF_PRESENT,
  72             IOobject::AUTO_WRITE
  73         ),
  74         stressMesh,
  75         dimensionedSymmTensor("zero", dimForce/dimArea, symmTensor::zero)
  76     );
  77
  78
  79     Info << "Reading incremental stress field DSigma\n" << endl;
  80     volSymmTensorField DSigma
  81     (
  82         IOobject
  83         (
  84             "DSigma",
  85             runTime.timeName(),
  86             stressMesh,
  87             IOobject::READ_IF_PRESENT,
  88             IOobject::AUTO_WRITE
  89         ),
  90         stressMesh,
  91         dimensionedSymmTensor("zero", dimForce/dimArea, symmTensor::zero)
  92     );
  93
  94     constitutiveModel rheology(sigma, DU);
  95     volScalarField rho = rheology.rho();
  96     volScalarField mu = rheology.mu();
  97     volScalarField lambda = rheology.lambda();
  98
  99     volTensorField F = I + gradDU.T();
 100     volTensorField DF = F - I;
 101     volScalarField J = det(F);
 102
 103
 104     word solidDdtScheme
 105     (
 106      stressMesh.schemesDict().ddtScheme("ddt(" + DU.name() +')')
 107     );
 108
 109 //     if
 110 //     (
 111 //         solidDdtScheme != fv::EulerFixedMeshDdtScheme<vector>::typeName
 112 //      && solidDdtScheme != fv::backwardFixedMeshDdtScheme<vector>::typeName
 113 //     )
 114 //     {
 115 //         FatalErrorIn(args.executable())
 116 //             << "Selected temporal differencing scheme: " << solidDdtScheme
 117 //             << ", instead: "
 118 //             << fv::EulerFixedMeshDdtScheme<vector>::typeName
 119 //             << " of "
 120 //             << fv::backwardFixedMeshDdtScheme<vector>::typeName
 121 //             << exit(FatalError);
 122 //     }
 123
 124
 125 //     IOdictionary rheologyProp
 126 //     (
 127 //         IOobject
 128 //         (
 129 //             "rheologyProperties",
 130 //             runTime.constant(),
 131 //             stressMesh,
 132 //             IOobject::MUST_READ,
 133 //             IOobject::NO_WRITE
 134 //         )
 135 //     );
 136
 137
 138 // dimensionedVector fb(rheologyProp.lookup("bodyForce"));