applications/solvers/solidMechanics/elasticNonLinTLSolidFoam/correctDirectionMixedTL.H

   1 {
   2   forAll(mesh.boundary(), patchID)
   3     {
   4       if(U.boundaryField()[patchID].type()
   5          == solidDirectionMixedFvPatchVectorField::typeName
   6          )
   7         {
   8           solidDirectionMixedFvPatchVectorField& loadingPatch =
   9             refCast<solidDirectionMixedFvPatchVectorField>
  10             (
  11              U.boundaryField()[patchID]
  12              );
  13
  14           tensorField Finv = inv(I + gradU);
  15           vectorField newN = Finv & n.boundaryField()[patchID];
  16           newN /= mag(newN);
  17           loadingPatch.valueFraction() = sqr(newN);
  18
  19           //- set gradient
  20           loadingPatch.refGrad() =
  21             (
  22              //Traction
  23              ( (mu.boundaryField()[patchID] + lambda.boundaryField()[patchID]) * (n.boundaryField()[patchID] & gradU.boundaryField()[patchID]) )
  24              - ( mu.boundaryField()[patchID] * (n.boundaryField()[patchID] & gradU.boundaryField()[patchID].T()) )
  25              - ( mu.boundaryField()[patchID] * ( n.boundaryField()[patchID] & (gradU.boundaryField()[patchID] & gradU.boundaryField()[patchID].T()) ) )
  26              - ( lambda.boundaryField()[patchID] * tr(gradU.boundaryField()[patchID]) * n.boundaryField()[patchID] )
  27              - ( 0.5 * lambda.boundaryField()[patchID] * tr(gradU.boundaryField()[patchID] & gradU.boundaryField()[patchID].T()) * n.boundaryField()[patchID] )
  28              )
  29             /
  30             (2.0*mu.boundaryField()[patchID] + lambda.boundaryField()[patchID]);
  31         }
  32     }
  33 }