doc/salome/examples/transforming_meshes_ex10.py

   1 # Sew Side Elements
   2
   3 import geompy
   4 import smesh
   5
   6 # create two boxes
   7 box1 = geompy.MakeBox(0.,  0., 0., 10., 10., 10.)
   8 box2 = geompy.MakeBox(0., 15., 0., 20., 25., 10.)
   9
  10 EdgesList = geompy.SubShapeAll(box2, geompy.ShapeType["EDGE"])
  11
  12 aComp = geompy.MakeCompound([box1, box2])
  13 geompy.addToStudy(aComp, "Two boxes")
  14
  15 # create a mesh on two boxes
  16 mesh = smesh.Mesh(aComp, "Two faces : quadrangle mesh")
  17
  18 algo1D = mesh.Segment()
  19 algo1D.NumberOfSegments(2)
  20 algo2D = mesh.Quadrangle()
  21
  22 algo_local = mesh.Segment(EdgesList[8])
  23 algo_local.NumberOfSegments(4)
  24 algo_local.Propagation()
  25
  26 mesh.Compute()
  27
  28 # sew side elements
  29 # IDsOfSide1Elements, IDsOfSide2Elements,
  30 # NodeID1OfSide1ToMerge, NodeID1OfSide2ToMerge, NodeID2OfSide1ToMerge, NodeID2OfSide2ToMerge
  31 mesh.SewSideElements([69, 70, 71, 72], [91, 92, 89, 90], 8, 38, 23, 58)