applications/solvers/solidMechanics/elasticOrthoNonLinULSolidFoam/checkPlaneStress.H

   1 if(rheology.planeStress())
   2   {
   3     //- add higher order terms
   4     volScalarField higherTerms = -0.5*volTensorField(gradDU&gradDU.T()).component(tensor::ZZ);
   5
   6     forAll(gradDU.internalField(), celli)
   7       {
   8         gradDU.internalField()[celli][tensor::ZZ] =
   9           ((-C.internalField()[celli][symmTensor4thOrder::XXZZ]*DEpsilon.internalField()[celli][symmTensor::XX]
  10             - C.internalField()[celli][symmTensor4thOrder::YYZZ]*DEpsilon.internalField()[celli][symmTensor::YY])
  11                     /
  12            C.internalField()[celli][symmTensor4thOrder::ZZZZ])
  13           -higherTerms.internalField()[celli];
  14       }
  15     forAll(gradDU.boundaryField(), patchi)
  16       {
  17         forAll(gradDU.boundaryField()[patchi], facei)
  18           {
  19             gradDU.boundaryField()[patchi][facei][tensor::ZZ] =
  20               ((-C.boundaryField()[patchi][facei][symmTensor4thOrder::XXZZ]*DEpsilon.boundaryField()[patchi][facei][symmTensor::XX]
  21                 - C.boundaryField()[patchi][facei][symmTensor4thOrder::YYZZ]*DEpsilon.boundaryField()[patchi][facei][symmTensor::YY])
  22                         /
  23                C.boundaryField()[patchi][facei][symmTensor4thOrder::ZZZZ])
  24               - higherTerms.boundaryField()[patchi][facei];
  25           }
  26       }
  27   }