2 # Tetrahedrization of a simple box. Hypothesis and algorithms for
3 # the mesh generation are global
12 geom = salome.lcc.FindOrLoadComponent("FactoryServer", "GEOM")
13 smesh = salome.lcc.FindOrLoadComponent("FactoryServer", "SMESH")
15 smeshgui = salome.ImportComponentGUI("SMESH")
16 smeshgui.Init(salome.myStudyId);
18 # -----------------------------------------------------------------------------
25 box = geompy.MakeBox(0., 0., 0., 100., 200., 300.)
27 idbox = geompy.addToStudy(box,"box")
29 print "Analysis of the geometry box :"
30 subShellList=geompy.SubShapeAll(box,ShapeTypeShell)
31 subFaceList=geompy.SubShapeAll(box,ShapeTypeFace)
32 subEdgeList=geompy.SubShapeAll(box,ShapeTypeEdge)
34 print "number of Shells in box : ",len(subShellList)
35 print "number of Faces in box : ",len(subFaceList)
36 print "number of Edges in box : ",len(subEdgeList)
39 ### ---------------------------- SMESH --------------------------------------
41 # ---- create Hypothesis
43 print "-------------------------- create Hypothesis"
45 print "-------------------------- NumberOfSegments"
49 hypNbSeg=smesh.CreateHypothesis("NumberOfSegments", "libStdMeshersEngine.so")
50 hypNbSeg.SetNumberOfSegments(numberOfSegments)
52 print hypNbSeg.GetName()
53 print hypNbSeg.GetId()
54 print hypNbSeg.GetNumberOfSegments()
56 smeshgui.SetName(salome.ObjectToID(hypNbSeg), "NumberOfSegments_10")
58 print "-------------------------- MaxElementArea"
62 hypArea=smesh.CreateHypothesis("MaxElementArea", "libStdMeshersEngine.so")
63 hypArea.SetMaxElementArea(maxElementArea)
65 print hypArea.GetName()
67 print hypArea.GetMaxElementArea()
69 smeshgui.SetName(salome.ObjectToID(hypArea), "MaxElementArea_500")
71 print "-------------------------- MaxElementVolume"
73 maxElementVolume = 500
75 hypVolume=smesh.CreateHypothesis("MaxElementVolume", "libStdMeshersEngine.so")
76 hypVolume.SetMaxElementVolume(maxElementVolume)
78 print hypVolume.GetName()
79 print hypVolume.GetId()
80 print hypVolume.GetMaxElementVolume()
82 smeshgui.SetName(salome.ObjectToID(hypVolume), "MaxElementVolume_500")
84 # ---- create Algorithms
86 print "-------------------------- create Algorithms"
88 print "-------------------------- Regular_1D"
90 regular1D=smesh.CreateHypothesis("Regular_1D", "libStdMeshersEngine.so")
91 smeshgui.SetName(salome.ObjectToID(regular1D), "Wire Discretisation")
93 print "-------------------------- MEFISTO_2D"
95 mefisto2D=smesh.CreateHypothesis("MEFISTO_2D", "libStdMeshersEngine.so")
96 smeshgui.SetName(salome.ObjectToID(mefisto2D), "MEFISTO_2D")
98 print "-------------------------- NETGEN_3D"
100 netgen3D=smesh.CreateHypothesis("NETGEN_3D", "libNETGENEngine.so")
101 smeshgui.SetName(salome.ObjectToID(netgen3D), "NETGEN_3D")
103 # ---- init a Mesh with the boxe
105 mesh = smesh.CreateMesh(box)
106 smeshgui.SetName(salome.ObjectToID(mesh), "MeshBox")
108 # ---- add hypothesis to the boxe
110 print "-------------------------- add hypothesis to the box"
112 mesh.AddHypothesis(box,regular1D)
113 mesh.AddHypothesis(box,hypNbSeg)
115 mesh.AddHypothesis(box,mefisto2D)
116 mesh.AddHypothesis(box,hypArea)
118 mesh.AddHypothesis(box,netgen3D)
119 mesh.AddHypothesis(box,hypVolume)
121 salome.sg.updateObjBrowser(1)
123 print "-------------------------- compute the mesh of the boxe"
124 ret=smesh.Compute(mesh,box)
127 log=mesh.GetLog(0) # no erase trace
130 print "Information about the MeshBox:"
131 print "Number of nodes : ", mesh.NbNodes()
132 print "Number of edges : ", mesh.NbEdges()
133 print "Number of faces : ", mesh.NbFaces()
134 print "Number of triangles : ", mesh.NbTriangles()
135 print "Number of volumes : ", mesh.NbVolumes()
136 print "Number of tetrahedrons: ", mesh.NbTetras()
138 print "probleme when computing the mesh"