applications/solvers/solidMechanics/viscoElasticSolidFoam/updateLoadingPatchTraction.H

   1 // Loading patch
   2 {
   3     label patchID = mesh.boundaryMesh().findPatchID("loading");
   4
   5     if(patchID == -1)
   6     {
   7         FatalErrorIn(args.executable())
   8             << "Can't find patch with name loading" << abort(FatalError);
   9     }
  10
  11     if
  12     (
  13         DU.boundaryField()[patchID].type()
  14      == directionMixedFvPatchVectorField::typeName
  15     )
  16     {
  17         directionMixedFvPatchVectorField& loadingPatchU =
  18             refCast<directionMixedFvPatchVectorField>
  19             (
  20                 DU.boundaryField()[patchID]
  21             );
  22
  23         vectorField n = mesh.boundary()[patchID].nf();
  24
  25         vectorField traction = n*(n&(n&DSigma));
  26
  27         scalarField lpMu = mu.boundaryField()[patchID];
  28         scalarField lpLambda = lambda.boundaryField()[patchID];
  29
  30         tensorField lpGradDU = gradDU.boundaryField()[patchID];
  31
  32         loadingPatchU.refGrad() =
  33         (
  34             traction
  35           - (n & (lpMu*lpGradDU.T() - (lpMu + lpLambda)*lpGradDU))
  36           - n*lpLambda*tr(lpGradDU)
  37         )
  38        /(2.0*lpMu + lpLambda);
  39     }
  40 }
  41
  42 // Bottom patch
  43 {
  44     label patchID = mesh.boundaryMesh().findPatchID("bottom");
  45
  46     if(patchID == -1)
  47     {
  48         FatalErrorIn(args.executable())
  49             << "Can't find patch with name bottom" << abort(FatalError);
  50     }
  51
  52     if
  53     (
  54         DU.boundaryField()[patchID].type()
  55      == directionMixedFvPatchVectorField::typeName
  56     )
  57     {
  58         directionMixedFvPatchVectorField& bottomPatchU =
  59             refCast<directionMixedFvPatchVectorField>
  60             (
  61                 DU.boundaryField()[patchID]
  62             );
  63
  64         vectorField n = mesh.boundary()[patchID].nf();
  65
  66         vectorField traction = n*(n&(n&DSigma));
  67
  68         scalarField bpMu = mu.boundaryField()[patchID];
  69         scalarField bpLambda = lambda.boundaryField()[patchID];
  70
  71         tensorField bpGradDU = gradDU.boundaryField()[patchID];
  72
  73         bottomPatchU.refGrad() =
  74         (
  75             traction
  76           - (n & (bpMu*bpGradDU.T() - (bpMu + bpLambda)*bpGradDU))
  77           - n*bpLambda*tr(bpGradDU)
  78         )
  79        /(2.0*bpMu + bpLambda);
  80     }
  81 }