doc/examples/quality_controls_ex04.py

   1 # Free Edges
   2
   3 from mechanic import *
   4
   5 aFilterMgr = smesh_builder.CreateFilterManager()
   6
   7 # Remove some elements to obtain free edges
   8 # Criterion : AREA > 95.
   9 area_margin = 95.
  10
  11 aFilter = smesh_builder.GetFilter(SMESH.FACE, SMESH.FT_Area, SMESH.FT_MoreThan, area_margin)
  12
  13 anIds = mesh.GetIdsFromFilter(aFilter)
  14
  15 mesh.RemoveElements(anIds)
  16
  17 # Criterion : Free Edges
  18 aBorders = mesh.GetFreeBorders()
  19
  20 # create groups
  21 aGroupF = mesh.CreateEmptyGroup(SMESH.FACE, "Faces with free edges")
  22 aGroupN = mesh.CreateEmptyGroup(SMESH.NODE, "Nodes on free edges")
  23
  24 # fill groups with elements, corresponding to the criterion
  25 print("")
  26 print("Criterion: Free edges Nb = ", len(aBorders))
  27 for i in range(len(aBorders)):
  28   aBorder = aBorders[i]
  29   print("Face # ", aBorder.myElemId, " : Edge between nodes (", end=' ')
  30   print(aBorder.myPnt1, ", ", aBorder.myPnt2, ")")
  31
  32   aGroupF.Add([aBorder.myElemId])
  33   aGroupN.Add([aBorder.myPnt1, aBorder.myPnt2])