applications/solvers/solidMechanics/elasticNonLinULSolidFoam/rotateFields.H

   1 //--------------------------------------------------//
   2 //- rotate fields
   3 //--------------------------------------------------//
   4 {
   5   Info << "Rotating fields" << endl;
   6
   7   volTensorField F = I + gradDU;
   8
   9   U += DU;
  10
  11   epsilon += DEpsilon;
  12
  13   sigma += DSigma;
  14
  15   volTensorField Finv = inv(F);
  16
  17   volScalarField J = det(F);
  18   if(min(J.internalField()) < 0)
  19     {
  20       FatalErrorIn(args.executable())
  21         << "Negative Jacobian -> cell negative volumes!!!"
  22         << exit(FatalError);
  23     }
  24
  25   rho = rho/J;
  26
  27   n = mesh.Sf()/mesh.magSf();
  28
  29   //epsilon = symm(Finv & epsilon & Finv.T());
  30   epsilon = transform(Finv, epsilon);
  31
  32   //sigma  = 1/J * symm(F.T() & sigma & F);
  33   sigma = (1/J)*transform(F.T(), sigma);
  34 }