doc/salome/examples/creating_meshes_ex04.py

   1 # Editing of a mesh
   2
   3 import salome
   4 salome.salome_init()
   5 import GEOM
   6 from salome.geom import geomBuilder
   7 geompy = geomBuilder.New()
   8
   9 import SMESH, SALOMEDS
  10 from salome.smesh import smeshBuilder
  11 smesh =  smeshBuilder.New()
  12
  13 def PrintMeshInfo(theMesh):
  14     aMesh = theMesh.GetMesh()
  15     print("Information about mesh:")
  16     print("Number of nodes       : ", aMesh.NbNodes())
  17     print("Number of edges       : ", aMesh.NbEdges())
  18     print("Number of faces       : ", aMesh.NbFaces())
  19     print("Number of volumes     : ", aMesh.NbVolumes())
  20     pass
  21
  22 # create a box
  23 box = geompy.MakeBox(0., 0., 0., 20., 20., 20.)
  24 geompy.addToStudy(box, "box")
  25
  26 # select one edge of the box for definition of a local hypothesis
  27 subShapeList = geompy.SubShapeAll(box, geompy.ShapeType["EDGE"])
  28 edge = subShapeList[0]
  29 name = geompy.SubShapeName(edge, box)
  30 geompy.addToStudyInFather(box, edge, name)
  31
  32 # create a mesh
  33 tria = smesh.Mesh(box, "Mesh 2D")
  34 algo1D = tria.Segment()
  35 hyp1 = algo1D.NumberOfSegments(3)
  36 algo2D = tria.Triangle()
  37 hyp2 = algo2D.MaxElementArea(10.)
  38
  39 # create a sub-mesh
  40 algo_local = tria.Segment(edge)
  41 hyp3 = algo_local.Arithmetic1D(1, 6)
  42 hyp4 = algo_local.Propagation()
  43
  44 # compute the mesh
  45 tria.Compute()
  46 PrintMeshInfo(tria)
  47
  48 # remove a local hypothesis
  49 tria.RemoveHypothesis(hyp4, edge)
  50
  51 # compute the mesh
  52 tria.Compute()
  53 PrintMeshInfo(tria)
  54
  55 # change the value of the 2D hypothesis
  56 hyp2.SetMaxElementArea(2.)
  57
  58 # compute the mesh
  59 tria.Compute()
  60 PrintMeshInfo(tria)