applications/solvers/solidMechanics/viscoElasticSolidFoam/updateTractionFreePatchTractionIncrement.H

   1
   2 {
   3     label patchID = mesh.boundaryMesh().findPatchID("freeSurface");
   4
   5     if(patchID == -1)
   6     {
   7         FatalErrorIn(args.executable())
   8             << "Can't find patch with name freeSurface" << abort(FatalError);
   9     }
  10
  11     if
  12     (
  13         DU.boundaryField()[patchID].type()
  14      != tractionDisplacementIncrementFvPatchVectorField::typeName
  15     )
  16     {
  17         FatalErrorIn(args.executable())
  18             << "Bounary condition on " << DU.name()
  19                 <<  " is "
  20                 << DU.boundaryField()[patchID].type()
  21                 << "for patch: " << mesh.boundaryMesh()[patchID].name()
  22                 <<  ", instead "
  23                 << tractionDisplacementIncrementFvPatchVectorField::typeName
  24                 << abort(FatalError);
  25     }
  26
  27     tractionDisplacementIncrementFvPatchVectorField& tractionPatch =
  28         refCast<tractionDisplacementIncrementFvPatchVectorField>
  29         (
  30             DU.boundaryField()[patchID]
  31         );
  32
  33     vectorField n = mesh.boundary()[patchID].nf();
  34
  35     vector traction = vector::zero;
  36     vector DTraction = vector::zero;
  37
  38     vectorField relaxedTraction =
  39         (
  40             n
  41           & (
  42                 sigma.boundaryField()[patchID]
  43               + DSigmaCorr.boundaryField()[patchID]
  44             )
  45         );
  46
  47     tractionPatch.DTraction() = DTraction + (traction - relaxedTraction);
  48
  49     tractionPatch.DPressure() = 0.0;
  50 }
  51
  52