doc/salome/examples/transforming_meshes_ex10.py

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