applications/solvers/solidMechanics/elasticOrthoSolidFoam/setPlaneStressGradU.H

   1 //- set gradU.zz() for plane stress
   2 if(rheology.planeStress())
   3 {
   4     forAll(gradU.internalField(), celli)
   5     {
   6         gradU.internalField()[celli].zz() =
   7             (
   8               - C.internalField()[celli].xxzz()*epsilon.internalField()[celli].xx()
   9               - C.internalField()[celli].yyzz()*epsilon.internalField()[celli].yy()
  10             )
  11             /
  12             C.internalField()[celli].zzzz();
  13     }
  14     gradU.correctBoundaryConditions();
  15 }