doc/salome/examples/transforming_meshes_ex07.py

   1 # Sew Meshes Border to Side
   2
   3 import geompy
   4 import smesh
   5
   6 # create two faces of a box
   7 box1 = geompy.MakeBox(0., 0., -10., 30., 20., 25.)
   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(9)
  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 border to side
  35 # FirstNodeIDOnFreeBorder, SecondNodeIDOnFreeBorder, LastNodeIDOnFreeBorder,
  36 # FirstNodeIDOnSide, LastNodeIDOnSide,
  37 # CreatePolygons, CreatePolyedrs
  38 mesh.SewBorderToSide(5, 45, 6, 113, 109, 0, 0)