applications/solvers/solidMechanics/elasticPlasticNonLinULSolidFoam/rotateFields.H

   1 //--------------------------------------------------//
   2 //- rotate fields
   3 //--------------------------------------------------//
   4 {
   5   Info << "Rotating fields" << endl;
   6
   7   volTensorField F = I + gradDU;
   8
   9   volTensorField Finv = inv(F);
  10
  11   volScalarField J = det(F);
  12   if(min(J.internalField()) < 0)
  13     {
  14       FatalErrorIn(args.executable())
  15         << "Negative Jacobian -> cell negative volumes!!!"
  16         << exit(FatalError);
  17     }
  18
  19   rho = rho/J;
  20   n = mesh.Sf()/mesh.magSf();
  21
  22   // rotate strain and stress fields
  23   //epsilon = symm(Finv & epsilon & Finv.T());
  24   epsilon = transform(Finv, epsilon);
  25
  26   //sigma  = 1/J * symm(F.T() & sigma & F);
  27   sigma = (1/J)*transform(F.T(), sigma);
  28
  29   // Note: for strict large strain approach, the constitutive stiffness
  30   // would be rotated here - see elasticOrthoNonLinULSolidFoam.
  31 }