doc/salome/examples/filters_ex39.py

   1 # Elements of a domain
   2
   3 import salome, SMESH
   4 salome.salome_init()
   5 from salome.geom import geomBuilder
   6 geompy = geomBuilder.New()
   7 from salome.smesh import smeshBuilder
   8 smesh =  smeshBuilder.New()
   9
  10 # create two boxes to have two domains in the mesh
  11
  12 box1 = geompy.MakeBoxDXDYDZ( 100,100,100 )
  13 box2 = geompy.MakeTranslation( box1, 200, 0, 0 )
  14 boxes = geompy.MakeCompound( [box1, box2] )
  15 box1, box2 = geompy.SubShapeAll( boxes, geompy.ShapeType["SHAPE"], "box")
  16
  17 vertex = geompy.SubShape( box1, geompy.ShapeType["VERTEX"], [1] )
  18
  19 # create a mesh
  20
  21 mesh = smesh.Mesh( boxes )
  22 mesh.Segment(box1).NumberOfSegments( 5 )  # to have different nb of elements on the boxes
  23 mesh.Segment(box2).NumberOfSegments( 10 )
  24 mesh.Quadrangle()
  25 mesh.Hexahedron()
  26 mesh.Compute()
  27
  28 # Create filters with FT_ConnectedElements criterion by pointing a domain in different ways:
  29
  30 # using point coordinates in box_1
  31 nodeFilter = smesh.GetFilter( SMESH.NODE, SMESH.FT_ConnectedElements, "=", "1.,2,10", mesh=mesh )
  32 print("Nb. nodes in box_1:", len( nodeFilter.GetIDs()))
  33
  34 # using point coordinates in box_2
  35 edgeFilter = smesh.GetFilter( SMESH.EDGE, SMESH.FT_ConnectedElements, "=", [202,1,1 ], mesh=mesh )
  36 print("Nb. segments in box_2:", len( edgeFilter.GetIDs()))
  37
  38 # using a geom vertex of box_1
  39 faceFilter = smesh.GetFilter( SMESH.FACE, SMESH.FT_ConnectedElements, "=", vertex, mesh=mesh )
  40 print("Nb. faces in box_1:", len( edgeFilter.GetIDs()))
  41
  42 # using node ID in box_2
  43 voluFilter = smesh.GetFilter( SMESH.VOLUME, SMESH.FT_ConnectedElements, "=", 10, mesh=mesh )
  44 print("Nb. volumes in box_2:", len( voluFilter.GetIDs()))
  45