applications/solvers/solidMechanics/elasticOrthoAcpSolidFoam/updateCrackPropagationPatchesInternalFaces.H

   1 {
   2     labelHashSet selectedFaces;
   3
   4     const labelListList& edgeFaces = mesh.edgeFaces();
   5
   6     const scalarField& cohesiveZoneI = cohesiveZone.internalField();
   7
   8     forAll (crackPropagationPatches, patchI)
   9     {
  10         label curPatchID = crackPropagationPatches[patchI];
  11
  12         const labelList& patchEdges =
  13             mesh.boundaryMesh()[curPatchID].meshEdges();
  14
  15         forAll(patchEdges, edgeI)
  16         {
  17             label curEdge = patchEdges[edgeI];
  18
  19             const labelList& curFaces = edgeFaces[curEdge];
  20
  21             forAll(curFaces, faceI)
  22             {
  23                 if (curFaces[faceI] < mesh.nInternalFaces())
  24                 {
  25                     if (cohesiveZoneI[curFaces[faceI]] > (1.0 - SMALL))
  26                     {
  27                         if (!selectedFaces.found(curFaces[faceI]))
  28                         {
  29                             selectedFaces.insert(curFaces[faceI]);
  30                         }
  31                     }
  32                 }
  33             }
  34         }
  35     }
  36
  37
  38 //     if (notchPatchID != -1)
  39 //     {
  40 //         const labelList& notchPatchEdges =
  41 //             mesh.boundaryMesh()[notchPatchID].meshEdges();
  42
  43 //         forAll(notchPatchEdges, edgeI)
  44 //         {
  45 //             label curEdge = notchPatchEdges[edgeI];
  46
  47 //             const labelList& curFaces = edgeFaces[curEdge];
  48
  49 //             forAll(curFaces, faceI)
  50 //             {
  51 //                 if (curFaces[faceI] < mesh.nInternalFaces())
  52 //                 {
  53 //                     if (cohesiveZoneI[curFaces[faceI]] > (1.0 - SMALL))
  54 //                     {
  55 //                         selectedFaces.insert(curFaces[faceI]);
  56 //                     }
  57 //                 }
  58 //             }
  59 //         }
  60 //     }
  61
  62 //     if (mesh.boundary()[cohesivePatchID].size())
  63 //     {
  64 //         const labelList& cohesivePatchEdges =
  65 //             mesh.boundaryMesh()[cohesivePatchID].meshEdges();
  66
  67 //         forAll(cohesivePatchEdges, edgeI)
  68 //         {
  69 //             label curEdge = cohesivePatchEdges[edgeI];
  70
  71 //             const labelList& curFaces = edgeFaces[curEdge];
  72
  73 //             forAll(curFaces, faceI)
  74 //             {
  75 //                 if (curFaces[faceI] < mesh.nInternalFaces())
  76 //                 {
  77 //                     if (cohesiveZoneI[curFaces[faceI]] > (1.0 - SMALL))
  78 //                     {
  79 //                         selectedFaces.insert(curFaces[faceI]);
  80 //                     }
  81 //                 }
  82 //             }
  83 //         }
  84 //     }
  85
  86     crackPropagationPatchesInternalFaces = labelList(selectedFaces.toc());
  87 }