doc/salome/examples/grouping_elements_ex04.py

   1 # Edit a Group
   2
   3 import SMESH_mechanic
   4
   5 smesh  = SMESH_mechanic.smesh
   6 mesh   = SMESH_mechanic.mesh
   7 salome = SMESH_mechanic.salome
   8
   9 # Get ids of all faces with area > 35
  10 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_Area, smesh.FT_MoreThan, 35.)
  11
  12 anIds = mesh.GetIdsFromFilter(aFilter)
  13
  14 print "Criterion: Area > 35, Nb = ", len(anIds)
  15
  16 # create a group by adding elements with area > 35
  17 aGroup = mesh.CreateEmptyGroup(smesh.FACE, "Area > 35")
  18 aGroup.Add(anIds)
  19
  20 # Get ids of all faces with area > 40
  21 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_Area, smesh.FT_MoreThan, 40.)
  22
  23 anIds = mesh.GetIdsFromFilter(aFilter)
  24
  25 print "Criterion: Area > 40, Nb = ", len(anIds)
  26
  27 # create a group of elements with area [35; 40] by removing elements with area > 40 from group aGroup
  28 aGroup.Remove(anIds)
  29
  30 # print the result
  31 aGroupElemIDs = aGroup.GetListOfID()
  32
  33 print "Criterion: 35 < Area < 40, Nb = ", len(aGroupElemIDs)
  34
  35 j = 1
  36 for i in range(len(aGroupElemIDs)):
  37   if j > 20: j = 1; print ""
  38   print aGroupElemIDs[i],
  39   j = j + 1
  40   pass
  41 print ""
  42
  43 salome.sg.updateObjBrowser(1)