doc/salome/examples/modifying_meshes_ex03.py

   1 # Add 0D Element on Element Nodes
   2
   3 import smesh, SMESH, geompy
   4
   5 # create a geometry
   6 box = geompy.MakeBoxDXDYDZ( 10, 10, 10 )
   7 face = geompy.SubShapeAll( box, geompy.ShapeType["FACE"])[0]
   8
   9 # make 3D mesh
  10 mesh = smesh.Mesh( box )
  11 mesh.AutomaticHexahedralization(0)
  12
  13 # create 0D elements on all nodes of the mesh
  14 res = mesh.Add0DElementsToAllNodes( mesh )
  15
  16 # find 0D elements on all nodes of the mesh, all found nodes are added to a new group
  17 groupName = "0Dmesh"
  18 res = mesh.Add0DElementsToAllNodes( mesh, groupName )
  19 mesh.RemoveGroupWithContents( res ) # remove all found 0D elements
  20
  21 # create 0D elements on all nodes of a sub-mesh, with group creation
  22 groupName = "0Dsubmesh"
  23 submesh = mesh.GetSubMesh( face, "faceSM")
  24 res = mesh.Add0DElementsToAllNodes( submesh, groupName )
  25
  26 # create 0D elements on all nodes of a group
  27 group = mesh.Group( face, "faceGroup" )
  28 res = mesh.Add0DElementsToAllNodes( group )
  29
  30 # remove all 0D elements
  31 mesh.RemoveElements( mesh.GetIdsFromFilter( smesh.GetFilter( SMESH.ELEM0D,
  32                                                              SMESH.FT_ElemGeomType,
  33                                                              "=",SMESH.Geom_POINT )))
  34
  35 # create 0D elements on all nodes of some elements
  36 res = mesh.Add0DElementsToAllNodes( mesh.GetElementsId() )
  37
  38 mesh.RemoveElements( mesh.GetElementsByType( SMESH.ELEM0D ))
  39
  40 # create 0D elements on some nodes
  41 nodes = range(1,10)
  42 res = mesh.Add0DElementsToAllNodes( mesh.GetIDSource( nodes, SMESH.NODE ))