3 \page tui_viewing_meshes_page Viewing Meshes
6 \anchor tui_viewing_mesh_infos
7 <h2>Viewing Mesh Infos</h2>
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"])
20 tetra = smesh.Mesh(box, "MeshBox")
22 algo1D = tetra.Segment()
23 algo1D.NumberOfSegments(3)
25 algo2D = tetra.Triangle()
26 algo2D.MaxElementArea(10.)
28 algo3D = tetra.Tetrahedron(smesh.NETGEN)
29 algo3D.MaxElementVolume(900.)
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()
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 ] )
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()
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()
66 print " %s : %d" % ( i, info[i] )
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()
74 print " %s : %d" % ( i, info[i] )
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()
82 print " %s : %d" % ( i, info[i] )
89 \anchor tui_find_element_by_point
90 <h2>Find Element by Point</h2>
97 # Create a geometry to mesh
98 box = geompy.MakeBoxDXDYDZ(100,100,100)
101 mesh = Mesh(box,"Mesh")
102 mesh.AutomaticHexahedralization()
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)
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)
117 # Find an edge at the point
118 edges = mesh.FindElementsByPoint(x,y,z, SMESH.EDGE )
119 assert( len(edges) == 1)
121 # Find faces at the point
122 edges = mesh.FindElementsByPoint(x,y,z, SMESH.FACE )
123 assert( len(edges) == 2)
125 # Find a volume at the point
126 vols = mesh.FindElementsByPoint(x,y,z, SMESH.VOLUME )
127 assert( len(vols) == 1)
129 # Find 0D elements at the point
130 edges = mesh.FindElementsByPoint(x,y,z, SMESH.ELEM0D )
131 assert( len(edges) == 0)