applications/solvers/solidMechanics/icoFsiElasticNonLinULSolidFoam/rotateSolidFields.H

   1 //--------------------------------------------------//
   2 //- rotate fields
   3 //--------------------------------------------------//
   4 {
   5   Info << "Rotating fields" << endl;
   6
   7   F = I + DF;
   8
   9   //U += DU;
  10   Usolid += DU;
  11
  12   //epsilon += DEpsilon;
  13
  14   sigma += DSigma;
  15
  16   //volTensorField Finv = inv(F);
  17
  18   J = det(F);
  19
  20   rho = rho/J;
  21
  22   //epsilon = symm(Finv.T() & epsilon & Finv);
  23
  24   sigma  = 1/J * symm(F & sigma & F.T());
  25 }