doc/examples/filters_ex18.py

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