applications/solvers/solidMechanics/elasticOrthoNonLinULSolidFoam/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 - a cell volume has become negative!"
  22         << exit(FatalError);
  23     }
  24
  25   rho = rho/J;
  26
  27   n = mesh.Sf()/mesh.magSf();
  28
  29   //- rotate strain
  30   //epsilon = symm(Finv & epsilon & Finv.T());
  31   epsilon = transform(Finv, epsilon);
  32
  33   //- rotate stress
  34   //sigma  = 1/J * symm(F.T() & sigma & F);
  35   sigma = (1/J) * transform(F.T(), sigma);
  36
  37   //- rotate elastic constitutive tensor
  38   C = transform(F.T(), C);
  39
  40   // - update implicit stiffness tensor
  41   forAll(K, celli)
  42     {
  43       K[celli].xx() = C[celli].xxxx();
  44       K[celli].yy() = C[celli].yyyy();
  45       K[celli].zz() = C[celli].zzzz();
  46     }
  47   K.correctBoundaryConditions();
  48 }