doc/salome/examples/filters_ex18.py

   1 # Borders at multi-connection
   2
   3 import salome
   4 salome.salome_init_without_session()
   5 from salome.geom import geomBuilder
   6 geompy = geomBuilder.New()
   7 import SMESH
   8 from salome.smesh import smeshBuilder
   9 smesh =  smeshBuilder.New()
  10
  11 # make a mesh on a box
  12 box = geompy.MakeBoxDXDYDZ(100,100,100)
  13 mesh = smesh.Mesh( box, "Box" )
  14 mesh.Segment().NumberOfSegments(10)
  15 mesh.Quadrangle()
  16 mesh.Hexahedron()
  17 mesh.Compute()
  18 # copy all elements with translation and merge nodes
  19 mesh.TranslateObject( mesh, [10,0,0], Copy=True )
  20 mesh.MergeNodes( mesh.FindCoincidentNodes( 1e-5 ))
  21
  22 # get mesh edges with number of connected elements (faces and volumes) == 3
  23 filter = smesh.GetFilter(SMESH.EDGE, SMESH.FT_MultiConnection, 3)
  24 ids = mesh.GetIdsFromFilter(filter)
  25 print("Number of border edges with 3 faces connected:", len(ids))