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                     (-C.internalField()[celli].xxzz()*epsilon.internalField()[celli].xx()
   8                      - C.internalField()[celli].yyzz()*epsilon.internalField()[celli].yy())
   9                     /
  10                     C.internalField()[celli].zzzz();
  11                 }
  12               gradU.correctBoundaryConditions();
  13             }