4 ### This file shows how to get the information that are displayed when using Mesh Information for elements/nodes
11 import salome_notebook
12 notebook = salome_notebook.NoteBook()
13 sys.path.insert(0, r'/local00/home/B61570/work_in_progress/mesh_info')
20 from salome.geom import geomBuilder
25 geompy = geomBuilder.New()
27 Box_1 = geompy.MakeBoxDXDYDZ(200, 200, 200)
28 geompy.addToStudy( Box_1, 'Box_1' )
29 bottom = geompy.CreateGroup(Box_1, geompy.ShapeType["FACE"])
30 geompy.UnionIDs(bottom, [31])
36 import SMESH, SALOMEDS
37 from salome.smesh import smeshBuilder
39 smesh = smeshBuilder.New()
43 NETGEN_3D_Parameters_1 = smesh.CreateHypothesisByAverageLength( 'NETGEN_Parameters', 'NETGENEngine', 25, 0 )
44 Mesh_tetra = smesh.Mesh(Box_1,'Mesh_tetra')
45 status = Mesh_tetra.AddHypothesis( Box_1, NETGEN_3D_Parameters_1 )
46 NETGEN_1D_2D_3D = Mesh_tetra.Tetrahedron(algo=smeshBuilder.NETGEN_1D2D3D)
47 isDone = Mesh_tetra.Compute()
49 raise ("Could not compute mesh: "+Mesh_tetra.GetName())
52 Mesh_quadratic = smesh.Mesh(Box_1,'Mesh_quadratic')
53 NETGEN_1D_2D_3D_1 = Mesh_quadratic.Tetrahedron(algo=smeshBuilder.NETGEN_1D2D3D)
54 NETGEN_3D_Parameters_2 = NETGEN_1D_2D_3D_1.Parameters()
55 NETGEN_3D_Parameters_2.SetMaxSize( 34.641 )
56 NETGEN_3D_Parameters_2.SetMinSize( 0.34641 )
57 NETGEN_3D_Parameters_2.SetSecondOrder( 1 )
58 isDone = Mesh_quadratic.Compute()
60 raise ("Could not compute mesh: "+Mesh_quadratic.GetName())
64 Mesh_hexa = smesh.Mesh(Box_1,'Mesh_hexa')
65 Regular_1D = Mesh_hexa.Segment()
66 Number_of_Segments_1 = Regular_1D.NumberOfSegments(15)
67 Quadrangle_2D = Mesh_hexa.Quadrangle(algo=smeshBuilder.QUADRANGLE)
68 Hexa_3D = Mesh_hexa.Hexahedron(algo=smeshBuilder.Hexa)
69 isDone = Mesh_hexa.Compute()
71 raise ("Could not compute mesh: "+Mesh_hexa.GetName())
74 Mesh_poly = smesh.CreateDualMesh(Mesh_tetra, 'dual_Mesh_1', True)
78 Mesh_prism = smesh.Mesh(Box_1,'Mesh_prism')
79 Regular_1D_1 = Mesh_prism.Segment()
80 Number_of_Segments_2 = Regular_1D_1.NumberOfSegments(15)
81 NETGEN_1D_2D = Mesh_prism.Triangle(algo=smeshBuilder.NETGEN_1D2D,geom=bottom)
82 NETGEN_2D_Parameters_1 = NETGEN_1D_2D.Parameters()
83 NETGEN_2D_Parameters_1.SetMaxSize(35)
84 NETGEN_2D_Parameters_1.SetMinSize(0.3)
85 Prism_3D = Mesh_prism.Prism()
86 isDone = Mesh_prism.Compute()
88 raise ("Could not compute mesh: "+Mesh_prism.GetName())
91 Mesh_pyramids = smesh.Mesh(Box_1,'Mesh_pyramids')
92 Regular_1D_2 = Mesh_pyramids.Segment()
93 Number_of_Segments_3 = Regular_1D_2.NumberOfSegments(15)
94 Quadrangle_2D_1 = Mesh_pyramids.Quadrangle(algo=smeshBuilder.QUADRANGLE)
95 NETGEN_3D = Mesh_pyramids.Tetrahedron()
96 bottom_1 = Mesh_pyramids.GroupOnGeom(bottom,'bottom',SMESH.FACE)
97 isDone = Mesh_pyramids.Compute()
99 raise ("Could not compute mesh: "+Mesh_pyramids.GetName())
102 ## Set names of Mesh objects
103 smesh.SetName(NETGEN_1D_2D_3D.GetAlgorithm(), 'NETGEN 1D-2D-3D')
104 smesh.SetName(NETGEN_3D_Parameters_1, 'NETGEN 3D Parameters_1')
105 smesh.SetName(Mesh_tetra.GetMesh(), 'Mesh_tetra')
106 smesh.SetName(Mesh_hexa.GetMesh(), 'Mesh_hexa')
108 if salome.sg.hasDesktop():
109 salome.sg.updateObjBrowser()
111 # Look in SMESH_GUI/SMESHGUI_MeshInfo.cxx +1666 for list of what is
113 def face_info(mesh, elem_id):
115 Print equivalent of Mesh Information for a face
117 elem_type = mesh.GetElementGeomType(elem_id)
119 conn = mesh.GetElemNodes(elem_id)
123 position = mesh.GetElementPosition(elem_id)
124 pos = f"{position.shapeType} #{position.shapeID}"
126 grav_center = mesh.BaryCenter(elem_id)
128 normal = mesh.GetFaceNormal(elem_id, normalized=True)
130 aspect_ratio = mesh.GetAspectRatio(elem_id)
131 #aspect_ratio = mesh.FunctorValue(SMESH.FT_AspectRatio, elem_id, isElem=True)
133 warping = mesh.GetWarping(elem_id)
134 #warping = mesh.FunctorValue(SMESH.FT_Warping, elem_id, isElem=True)
136 min_angle = mesh.GetMinimumAngle(elem_id)
137 #min_angle = mesh.FunctorValue(SMESH.FT_MinimumAngle, elem_id, isElem=True)
139 taper = mesh.GetTaper(elem_id)
140 #taper = mesh.FunctorValue(SMESH.FT_Taper, elem_id, isElem=True)
142 skew = mesh.GetSkew(elem_id)
143 #skew = mesh.FunctorValue(SMESH.FT_Skew, elem_id, isElem=True)
145 area = mesh.GetArea(elem_id)
146 #area = mesh.FunctorValue(SMESH.FT_Area, elem_id, isElem=True)
148 diameter = mesh.GetMaxElementLength(elem_id)
149 #diameter = mesh.FunctorValue(SMESH.FT_MaxElementLength2D, elem_id, isElem=True)
151 min_length = mesh.FunctorValue(SMESH.FT_Length2D, elem_id, isElem=True)
160 - X: {grav_center[0]}
161 - Y: {grav_center[1]}
162 - Z: {grav_center[2]}
168 - Aspect Ratio: {aspect_ratio}
170 - Minimum Angle: {min_angle}
174 - Element Diameter 2D: {diameter}
175 - Minimum Edge Length: {min_length}
179 def volume_info(mesh, elem_id):
181 Print equivalent of Mesh Information for a volume
183 elem_type = mesh.GetElementGeomType(elem_id)
185 if elem_type in [SMESH.Entity_Polyhedra, SMESH.Entity_Quad_Polyhedra]:
189 while face_conn != []:
190 face_conn = mesh.GetElemFaceNodes(elem_id, iface)
192 conn.append(face_conn)
193 nb_nodes = len(mesh.GetElemNodes(elem_id))
195 conn = mesh.GetElemNodes(elem_id)
199 position = mesh.GetElementPosition(elem_id)
200 pos = f"{position.shapeType} #{position.shapeID}"
202 grav_center = mesh.BaryCenter(elem_id)
204 aspect_ratio = mesh.GetAspectRatio(elem_id)
205 #aspect_ratio = mesh.FunctorValue(SMESH.FT_AspectRatio3D, elem_id, isElem=True)
207 volume = mesh.GetVolume(elem_id)
208 #volume = mesh.FunctorValue(SMESH.FT_Volume3D, elem_id, isElem=True)
210 jacob = mesh.GetScaledJacobian(elem_id)
211 #jacob = mesh.FunctorValue(SMESH.FT_ScaledJacobian, elem_id, isElem=True)
213 diameter = mesh.GetMaxElementLength(elem_id)
214 #diameter = mesh.FunctorValue(SMESH.FT_MaxElementLength3D, elem_id, isElem=True)
216 min_length = mesh.FunctorValue(SMESH.FT_Length3D, elem_id, isElem=True)
225 - X: {grav_center[0]}
226 - Y: {grav_center[1]}
227 - Z: {grav_center[2]}
229 - Aspect Ratio 3D: {aspect_ratio}
231 - Scaled Jacobian: {jacob}
232 - Element Diameter 3D: {diameter}
233 - Minimum Edge Length: {min_length}
237 def node_info(mesh, node_id):
239 coord = mesh.GetNodeXYZ(node_id)
241 conn_edge = mesh.GetNodeInverseElements(node_id, SMESH.EDGE)
242 conn_face = mesh.GetNodeInverseElements(node_id, SMESH.FACE)
243 conn_vol = mesh.GetNodeInverseElements(node_id, SMESH.VOLUME)
245 position = mesh.GetNodePosition(node_id)
246 pos = f"{position.shapeType} #{position.shapeID}"
249 vec[0:len(position.params)] = position.params
260 - Volumes: {conn_vol}
272 volume_info(Mesh_tetra, 3000)
274 volume_info(Mesh_hexa, 3000)
276 volume_info(Mesh_poly, 3000)
278 volume_info(Mesh_prism, 1400)
280 volume_info(Mesh_pyramids, 8176)
282 volume_info(Mesh_quadratic, 1180)
289 face_info(Mesh_tetra, 147)
291 face_info(Mesh_hexa, 1464)
293 face_info(Mesh_poly, 771)
295 face_info(Mesh_quadratic, 138)
302 node_info(Mesh_tetra, 152)
304 node_info(Mesh_tetra, 32)
306 node_info(Mesh_tetra, 2)