doc/salome/gui/SMESH/input/tui_viewing_meshes.doc

   1 /*!
   2
   3 \page tui_viewing_meshes_page Viewing Meshes
   4
   5 <br>
   6 \anchor tui_viewing_mesh_infos
   7 <h2>Viewing Mesh Infos</h2>
   8
   9 \code
  10 import geompy
  11 import smesh
  12 import SMESH
  13
  14 # create a box
  15 box = geompy.MakeBox(0., 0., 0., 20., 20., 20.)
  16 geompy.addToStudy(box, "box")
  17 [Face_1,Face_2,Face_3,Face_4,Face_5,Face_5] = geompy.SubShapeAll(box, geompy.ShapeType["FACE"])
  18
  19 # create a mesh
  20 tetra = smesh.Mesh(box, "MeshBox")
  21
  22 algo1D = tetra.Segment()
  23 algo1D.NumberOfSegments(3)
  24
  25 algo2D = tetra.Triangle()
  26 algo2D.MaxElementArea(10.)
  27
  28 algo3D = tetra.Tetrahedron(smesh.NETGEN)
  29 algo3D.MaxElementVolume(900.)
  30
  31 # Creation of SubMesh
  32 Regular_1D_1_1 = tetra.Segment(geom=Face_1)
  33 Nb_Segments_1 = Regular_1D_1_1.NumberOfSegments(5)
  34 Nb_Segments_1.SetDistrType( 0 )
  35 Quadrangle_2D = tetra.Quadrangle(geom=Face_1)
  36 isDone = tetra.Compute()
  37 submesh = Regular_1D_1_1.GetSubMesh()
  38
  39 # compute the mesh
  40 tetra.Compute()
  41
  42 # Creation of group
  43 group = tetra.CreateEmptyGroup( SMESH.FACE, 'Group' )
  44 nbAdd = group.Add( [ 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76 ] )
  45
  46 # Print information about the mesh
  47 print "Information about mesh:"
  48 print "Number of nodes       : ", tetra.NbNodes()
  49 print "Number of edges       : ", tetra.NbEdges()
  50 print "Number of faces       : ", tetra.NbFaces()
  51 print "          triangles   : ", tetra.NbTriangles()
  52 print "          quadrangles : ", tetra.NbQuadrangles()
  53 print "          polygons    : ", tetra.NbPolygons()
  54 print "Number of volumes     : ", tetra.NbVolumes()
  55 print "          tetrahedrons: ", tetra.NbTetras()
  56 print "          hexahedrons : ", tetra.NbHexas()
  57 print "          prisms      : ", tetra.NbPrisms()
  58 print "          pyramids    : ", tetra.NbPyramids()
  59 print "          polyhedrons : ", tetra.NbPolyhedrons()
  60
  61 # Get Information About Mesh by GetMeshInfo
  62 print "\nInformation about mesh by GetMeshInfo:"
  63 info = smesh.GetMeshInfo(tetra)
  64 keys = info.keys(); keys.sort()
  65 for i in keys:
  66   print "  %s   :  %d" % ( i, info[i] )
  67   pass
  68
  69 # Get Information About Group by GetMeshInfo
  70 print "\nInformation about group by GetMeshInfo:"
  71 info = smesh.GetMeshInfo(group)
  72 keys = info.keys(); keys.sort()
  73 for i in keys:
  74   print "  %s  :  %d" % ( i, info[i] )
  75   pass
  76
  77 # Get Information About SubMesh by GetMeshInfo
  78 print "\nInformation about Submesh by GetMeshInfo:"
  79 info = smesh.GetMeshInfo(submesh)
  80 keys = info.keys(); keys.sort()
  81 for i in keys:
  82   print "  %s  :  %d" % ( i, info[i] )
  83   pass
  84 \endcode
  85
  86
  87
  88 <br>
  89 \anchor tui_find_element_by_point
  90 <h2>Find Element by Point</h2>
  91
  92 \code
  93 import geompy
  94 import smesh
  95 import SMESH
  96
  97 # Create a geometry to mesh
  98 box = geompy.MakeBoxDXDYDZ(100,100,100)
  99
 100 # Create a mesh
 101 mesh = Mesh(box,"Mesh")
 102 mesh.AutomaticHexahedralization()
 103 mesh.Compute()
 104
 105 # Create a point
 106 x,y,z = 0, 0, 1
 107
 108 # Find all elements (except 0D ones) located at the point
 109 all_elems_except_0D = mesh.FindElementsByPoint(x,y,z)
 110 assert( len(all_elems_except_0D) == 4)
 111
 112 # Find nodes at the point
 113 nodes = mesh.FindElementsByPoint(x,y,z, SMESH.NODE )
 114 assert( len(nodes) == 0)
 115 assert( len( mesh.FindElementsByPoint(x,y,0, SMESH.NODE)) == 1)
 116
 117 # Find an edge at the point
 118 edges = mesh.FindElementsByPoint(x,y,z, SMESH.EDGE )
 119 assert( len(edges) == 1)
 120
 121 # Find faces at the point
 122 edges = mesh.FindElementsByPoint(x,y,z, SMESH.FACE )
 123 assert( len(edges) == 2)
 124
 125 # Find a volume at the point
 126 vols = mesh.FindElementsByPoint(x,y,z, SMESH.VOLUME )
 127 assert( len(vols) == 1)
 128
 129 # Find 0D elements at the point
 130 edges = mesh.FindElementsByPoint(x,y,z, SMESH.ELEM0D )
 131 assert( len(edges) == 0)
 132
 133 \endcode
 134
 135 */