doc/salome/examples/grouping_elements_ex05.py

   1 # Union of groups
   2
   3 import SMESH_mechanic
   4 import SMESH
   5
   6 smesh  = SMESH_mechanic.smesh
   7 mesh   = SMESH_mechanic.mesh
   8 salome = SMESH_mechanic.salome
   9
  10 # Criterion : AREA > 20
  11 aFilter = smesh.GetFilter(SMESH.FACE, SMESH.FT_Area, SMESH.FT_MoreThan, 20.)
  12
  13 anIds = mesh.GetIdsFromFilter(aFilter)
  14
  15 print "Criterion: Area > 20, Nb = ", len( anIds )
  16
  17 # create a group by adding elements with area > 20
  18 aGroup1 = mesh.CreateEmptyGroup(SMESH.FACE, "Area > 20")
  19 aGroup1.Add(anIds)
  20
  21 # Criterion : AREA = 20
  22 aFilter = smesh.GetFilter(SMESH.FACE, SMESH.FT_Area, SMESH.FT_EqualTo, 20.)
  23
  24 anIds = mesh.GetIdsFromFilter(aFilter)
  25
  26 print "Criterion: Area = 20, Nb = ", len( anIds )
  27
  28 # create a group by adding elements with area = 20
  29 aGroup2 = mesh.CreateEmptyGroup( SMESH.FACE, "Area = 20" )
  30
  31 aGroup2.Add(anIds)
  32
  33 # create union group : area >= 20
  34 aGroup3 = mesh.UnionListOfGroups([aGroup1, aGroup2], "Area >= 20")
  35 print "Criterion: Area >= 20, Nb = ", len(aGroup3.GetListOfID())
  36 # Please note that also there is UnionGroups() method which works with two groups only
  37
  38 # Criterion : AREA < 20
  39 aFilter = smesh.GetFilter(SMESH.FACE, SMESH.FT_Area, SMESH.FT_LessThan, 20.)
  40
  41 anIds = mesh.GetIdsFromFilter(aFilter)
  42
  43 print "Criterion: Area < 20, Nb = ", len(anIds)
  44
  45 # create a group by adding elements with area < 20
  46 aGroup4 = mesh.CreateEmptyGroup(SMESH.FACE, "Area < 20")
  47 aGroup4.Add(anIds)
  48
  49 # create union group : area >= 20 and area < 20
  50 aGroup5 = mesh.UnionListOfGroups([aGroup3, aGroup4], "Any Area")
  51 print "Criterion: Any Area, Nb = ", len(aGroup5.GetListOfID())
  52
  53 salome.sg.updateObjBrowser(True)