src/SMESH_SWIG/PAL_MESH_041_mesh.py

   1 import geompy
   2 import salome
   3
   4 import StdMeshers
   5
   6 #-----------------------------GEOM----------------------------------------
   7
   8 #----------Vertexes------------
   9 p1 = geompy.MakeVertex(20.0,30.0,40.0)
  10 p2 = geompy.MakeVertex(90.0,80.0,0.0)
  11 p3 = geompy.MakeVertex(30.0,80.0,200.0)
  12
  13 #----------Edges---------------
  14 e1 = geompy.MakeEdge(p1,p2)
  15 e2 = geompy.MakeEdge(p2,p3)
  16 e3 = geompy.MakeEdge(p3,p1)
  17
  18 #----------Wire----------------
  19 ListOfEdges = []
  20 ListOfEdges.append(e3)
  21 ListOfEdges.append(e2)
  22 ListOfEdges.append(e1)
  23 wire1 = geompy.MakeWire(ListOfEdges)
  24
  25
  26 #----------Face----------------
  27 WantPlanarFace = 1
  28 face1 = geompy.MakeFace(wire1,WantPlanarFace)
  29
  30 Id_face1 = geompy.addToStudy(face1,"Face1")
  31
  32
  33
  34 #-----------------------------SMESH-------------------------------------------
  35 smesh = salome.lcc.FindOrLoadComponent("FactoryServer", "SMESH")
  36
  37 # -- Init --
  38 plane_mesh = salome.IDToObject( Id_face1)
  39 smesh.SetCurrentStudy(salome.myStudy)
  40
  41 mesh = smesh.CreateMesh(plane_mesh)
  42
  43 smeshgui = salome.ImportComponentGUI("SMESH")
  44 smeshgui.Init(salome.myStudyId)
  45
  46 id_mesh = salome.ObjectToID(mesh)
  47 smeshgui.SetName( id_mesh, "Mesh_1")
  48
  49
  50 print"---------------------Hypothesis"
  51
  52
  53 #---------------- NumberOfSegments
  54 numberOfSegment = 9
  55
  56 hypNbSeg = smesh.CreateHypothesis( "NumberOfSegments", "libStdMeshersEngine.so" )
  57 hypNbSeg.SetNumberOfSegments( numberOfSegment )
  58
  59 print hypNbSeg.GetName()
  60 print hypNbSeg.GetNumberOfSegments()
  61 smeshgui.SetName(salome.ObjectToID(hypNbSeg), "Nb. Segments")
  62
  63
  64 #--------------------------Max. Element Area
  65 maxElementArea = 200
  66
  67 hypArea200 = smesh.CreateHypothesis("MaxElementArea","libStdMeshersEngine.so")
  68 hypArea200.SetMaxElementArea( maxElementArea )
  69 print hypArea200.GetName()
  70 print hypArea200.GetMaxElementArea()
  71
  72 smeshgui.SetName(salome.ObjectToID(hypArea200), "Max. Element Area")
  73
  74 print"---------------------Algorithms"
  75
  76 #----------------------------Wire discretisation
  77 algoWireDes = smesh.CreateHypothesis( "Regular_1D", "libStdMeshersEngine.so" )
  78 listHyp = algoWireDes.GetCompatibleHypothesis()
  79
  80 print algoWireDes.GetName()
  81 smeshgui.SetName(salome.ObjectToID(algoWireDes), "Ware descritisation")
  82
  83 #----------------------------Triangle (Mefisto)
  84 algoMef = smesh.CreateHypothesis( "MEFISTO_2D", "libStdMeshersEngine.so" )
  85 listHyp = algoMef.GetCompatibleHypothesis()
  86
  87 print algoMef.GetName()
  88
  89 #----------------------------Add hipothesis to the plane
  90 mesh.AddHypothesis( plane_mesh, hypNbSeg )     # nb segments
  91 mesh.AddHypothesis( plane_mesh, hypArea200 )   # max area
  92
  93 mesh.AddHypothesis( plane_mesh, algoWireDes )  # Regular 1D/wire discretisation
  94 mesh.AddHypothesis( plane_mesh, algoMef )      # MEFISTO 2D
  95
  96 smeshgui.SetName(salome.ObjectToID(algoMef), "Triangle (Mefisto)")
  97
  98 print "---------------------Compute the mesh"
  99
 100 smesh.Compute(mesh, plane_mesh)
 101
 102 salome.sg.updateObjBrowser(1)
 103