doc/examples/transforming_meshes_ex08.py

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