2 # Tetrahedrization of a simple box. Hypothesis and algorithms for
3 # the mesh generation are global
14 # -----------------------------------------------------------------------------
16 ##geom = salome.lcc.FindOrLoadComponent("FactoryServer", "Geometry")
17 ##myBuilder = salome.myStudy.NewBuilder()
20 myBuilder = geompy.myBuilder
28 box = geompy.MakeBox(0., 0., 0., 100., 200., 300.)
30 idbox = geompy.addToStudy(box,"box")
32 print "Analysis of the geometry box :"
33 subShellList=geompy.SubShapeAll(box,ShapeTypeShell)
34 subFaceList=geompy.SubShapeAll(box,ShapeTypeFace)
35 subEdgeList=geompy.SubShapeAll(box,ShapeTypeEdge)
37 print "number of Shells in box : ",len(subShellList)
38 print "number of Faces in box : ",len(subFaceList)
39 print "number of Edges in box : ",len(subEdgeList)
43 smeshgui = salome.ImportComponentGUI("SMESH")
44 smeshgui.Init(salome.myStudyId)
48 # ---- create Hypothesis
50 print "-------------------------- create Hypothesis"
52 print "-------------------------- NumberOfSegments"
56 hyp1=gen.CreateHypothesis("NumberOfSegments")
57 hypNbSeg=hyp1._narrow(SMESH.SMESH_NumberOfSegments)
58 hypNbSeg.SetNumberOfSegments(numberOfSegments)
59 hypNbSegID = hypNbSeg.GetId()
60 print hypNbSeg.GetName()
62 print hypNbSeg.GetNumberOfSegments()
64 idseg = smeshgui.AddNewHypothesis( salome.orb.object_to_string(hypNbSeg) )
65 smeshgui.SetName(idseg, "NumberOfSegments")
67 print "-------------------------- MaxElementArea"
71 hyp2=gen.CreateHypothesis("MaxElementArea")
72 hypArea=hyp2._narrow(SMESH.SMESH_MaxElementArea)
73 hypArea.SetMaxElementArea(maxElementArea)
74 print hypArea.GetName()
76 print hypArea.GetMaxElementArea()
78 idarea = smeshgui.AddNewHypothesis( salome.orb.object_to_string(hypArea) )
79 smeshgui.SetName(idarea, "MaxElementArea")
81 print "-------------------------- MaxElementVolume"
83 maxElementVolume = 500
85 hyp3=gen.CreateHypothesis("MaxElementVolume")
86 hypVolume=hyp3._narrow(SMESH.SMESH_MaxElementVolume)
87 hypVolume.SetMaxElementVolume(maxElementVolume)
88 print hypVolume.GetName()
89 print hypVolume.GetId()
90 print hypVolume.GetMaxElementVolume()
92 idvolume = smeshgui.AddNewHypothesis( salome.orb.object_to_string(hypVolume) )
93 smeshgui.SetName(idvolume, "MaxElementVolume")
95 # ---- create Algorithms
97 print "-------------------------- create Algorithms"
99 print "-------------------------- Regular_1D"
101 hypothesis=gen.CreateHypothesis("Regular_1D")
102 regular1D = hypothesis._narrow(SMESH.SMESH_Regular_1D)
103 regularID = smeshgui.AddNewAlgorithms( salome.orb.object_to_string(regular1D) )
104 smeshgui.SetName(regularID, "Wire Discretisation")
106 print "-------------------------- MEFISTO_2D"
108 hypothesis=gen.CreateHypothesis("MEFISTO_2D")
109 mefisto2D = hypothesis._narrow(SMESH.SMESH_MEFISTO_2D)
110 mefistoID = smeshgui.AddNewAlgorithms( salome.orb.object_to_string(mefisto2D) )
111 smeshgui.SetName(mefistoID, "MEFISTO_2D")
113 print "-------------------------- NETGEN_3D"
115 hypothesis=gen.CreateHypothesis("NETGEN_3D")
116 netgen3D = hypothesis._narrow(SMESH.SMESH_NETGEN_3D)
117 netgenID = smeshgui.AddNewAlgorithms( salome.orb.object_to_string(netgen3D) )
118 smeshgui.SetName(netgenID, "NETGEN_3D")
120 # ---- init a Mesh with the boxe
123 idmesh = smeshgui.AddNewMesh( salome.orb.object_to_string(mesh) )
124 smeshgui.SetName(idmesh, "MeshBox")
125 smeshgui.SetShape(idbox, idmesh)
127 # ---- add hypothesis to the boxe
129 print "-------------------------- add hypothesis to the boxe"
131 ret=mesh.AddHypothesis(box,regular1D)
133 ret=mesh.AddHypothesis(box,hypNbSeg)
135 ret=mesh.AddHypothesis(box,mefisto2D)
137 ret=mesh.AddHypothesis(box,hypArea)
139 ret=mesh.AddHypothesis(box,netgen3D)
141 ret=mesh.AddHypothesis(box,hypVolume)
144 smeshgui.SetAlgorithms( idmesh, regularID)
145 smeshgui.SetHypothesis( idmesh, idseg )
146 smeshgui.SetAlgorithms( idmesh, mefistoID )
147 smeshgui.SetHypothesis( idmesh, idarea )
148 smeshgui.SetAlgorithms( idmesh, netgenID )
149 smeshgui.SetHypothesis( idmesh, idvolume )
151 sg.updateObjBrowser(1)
154 print "-------------------------- compute the mesh of the boxe"
155 ret=gen.Compute(mesh,idbox)
157 log=mesh.GetLog(0) # no erase trace
162 sg.updateObjBrowser(1)