doc/salome/examples/creating_meshes_ex03.py

   1 # Change priority of submeshes in Mesh
   2
   3 import salome
   4 salome.salome_init()
   5 import GEOM
   6 from salome.geom import geomBuilder
   7 geompy = geomBuilder.New(salome.myStudy)
   8
   9 import SMESH, SALOMEDS
  10 from salome.smesh import smeshBuilder
  11 smesh =  smeshBuilder.New(salome.myStudy)
  12
  13 Box_1 = geompy.MakeBoxDXDYDZ(200, 200, 200)
  14 [Face_1,Face_2,Face_3,Face_4,Face_5,Face_6] = geompy.SubShapeAllSorted(Box_1, geompy.ShapeType["FACE"])
  15
  16 # create Mesh object on Box shape
  17 Mesh_1 = smesh.Mesh(Box_1)
  18
  19 # assign mesh algorithms
  20 Regular_1D = Mesh_1.Segment()
  21 Nb_Segments_1 = Regular_1D.NumberOfSegments(20)
  22 Nb_Segments_1.SetDistrType( 0 )
  23 MEFISTO_2D = Mesh_1.Triangle()
  24 Max_Element_Area_1 = MEFISTO_2D.MaxElementArea(1200)
  25 Tetrahedron = Mesh_1.Tetrahedron()
  26 Max_Element_Volume_1 = Tetrahedron.MaxElementVolume(40000)
  27
  28 # create submesh and assign algorithms on Face_1
  29 Regular_1D_1 = Mesh_1.Segment(geom=Face_1)
  30 Nb_Segments_2 = Regular_1D_1.NumberOfSegments(4)
  31 Nb_Segments_2.SetDistrType( 0 )
  32 MEFISTO_2D_1 = Mesh_1.Triangle(algo=smeshBuilder.MEFISTO,geom=Face_1)
  33 Length_From_Edges_2D = MEFISTO_2D_1.LengthFromEdges()
  34 SubMesh_1 = MEFISTO_2D_1.GetSubMesh()
  35
  36 # create submesh and assign algorithms on Face_2
  37 Regular_1D_2 = Mesh_1.Segment(geom=Face_2)
  38 Nb_Segments_3 = Regular_1D_2.NumberOfSegments(8)
  39 Nb_Segments_3.SetDistrType( 0 )
  40 MEFISTO_2D_2 = Mesh_1.Triangle(algo=smeshBuilder.MEFISTO,geom=Face_2)
  41 Length_From_Edges_2D_1 = MEFISTO_2D_2.LengthFromEdges()
  42 SubMesh_2 = MEFISTO_2D_2.GetSubMesh()
  43
  44 # create submesh and assign algorithms on Face_3
  45 Regular_1D_3 = Mesh_1.Segment(geom=Face_3)
  46 Nb_Segments_4 = Regular_1D_3.NumberOfSegments(12)
  47 Nb_Segments_4.SetDistrType( 0 )
  48 MEFISTO_2D_3 = Mesh_1.Triangle(algo=smeshBuilder.MEFISTO,geom=Face_3)
  49 Length_From_Edges_2D_2 = MEFISTO_2D_3.LengthFromEdges()
  50 SubMesh_3 = MEFISTO_2D_3.GetSubMesh()
  51
  52 # check exisiting submesh priority order
  53 [ [ SubMesh_1, SubMesh_3, SubMesh_2 ] ] = Mesh_1.GetMeshOrder()
  54 # set new submesh order
  55 isDone = Mesh_1.SetMeshOrder( [ [ SubMesh_1, SubMesh_2, SubMesh_3 ] ])
  56 # compute mesh
  57 isDone = Mesh_1.Compute()
  58
  59 # clear mesh result and compute with other submesh order
  60 Mesh_1.Clear()
  61 isDone = Mesh_1.SetMeshOrder( [ [ SubMesh_2, SubMesh_1, SubMesh_3 ] ])
  62 isDone = Mesh_1.Compute()