doc/salome/examples/transforming_meshes_ex09.py

   1 # Sew Free Borders
   2
   3 import geompy
   4 import smesh
   5
   6 # create two faces of the box
   7 box1 = geompy.MakeBox(0., 0., 0., 20., 20., 15.)
   8 facesList1 = geompy.SubShapeAll(box1, geompy.ShapeType["FACE"])
   9 face1 = facesList1[2]
  10
  11 box2 = geompy.MakeBox(0., 5., 0., 20., 20., 15.)
  12 facesList2 = geompy.SubShapeAll(box2, geompy.ShapeType["FACE"])
  13 face2 = facesList2[1]
  14
  15 edgesList = geompy.SubShapeAll(face2, geompy.ShapeType["EDGE"])
  16 edge1 = edgesList[2]
  17
  18 aComp = geompy.MakeCompound([face1, face2])
  19 geompy.addToStudy(aComp, "Two faces")
  20
  21 # create a mesh on two faces
  22 mesh = smesh.Mesh(aComp, "Two faces : quadrangle mesh")
  23
  24 algo1D = mesh.Segment()
  25 algo1D.NumberOfSegments(4)
  26 algo2D = mesh.Quadrangle()
  27
  28 algo_local = mesh.Segment(edge1)
  29 algo_local.Arithmetic1D(1, 4)
  30 algo_local.Propagation()
  31
  32 mesh.Compute()
  33
  34 # sew free borders
  35 # FirstNodeID1, SecondNodeID1, LastNodeID1,
  36 # FirstNodeID2, SecondNodeID2, LastNodeID2, CreatePolygons, CreatePolyedrs
  37 mesh.SewFreeBorders(6, 21, 5, 1, 12, 3, 0, 0)