doc/examples/transforming_meshes_ex09.py

   1 # Sew 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 # make two not sewed quadranges
  13 OY0 = geom_builder.MakeVectorDXDYDZ(0, 1, 0)
  14 OY1 = geom_builder.MakeTranslation( OY0, 1, 0, 0, theName="OY1" )
  15 OY2 = geom_builder.MakeTranslation( OY0, 1.01, 0, 0, theName="OY2" )
  16 OY3 = geom_builder.MakeTranslation( OY0, 2, 0, 0 )
  17 q1  = geom_builder.MakeQuad2Edges( OY0, OY1 )
  18 q2  = geom_builder.MakeQuad2Edges( OY2, OY3 )
  19
  20 shape = geom_builder.MakeCompound( [q1,q2], theName='shape' )
  21
  22 # make a non-uniform quadrangle mesh on two faces
  23 mesh = smesh_builder.Mesh(shape, "Two faces : quadrangle mesh")
  24 mesh.Segment().Arithmetic1D( 0.1, 0.4 )
  25 mesh.Segment(q1).NumberOfSegments( 5 )
  26 mesh.Quadrangle()
  27 mesh.Compute()
  28
  29 # sew free borders
  30
  31 segs1 = mesh.GetSubMeshElementsId( OY1 ) # mesh segments generated on borders
  32 segs2 = mesh.GetSubMeshElementsId( OY2 )
  33
  34 FirstNodeID1  = mesh.GetElemNode( segs1[0], 0 )
  35 SecondNodeID1 = mesh.GetElemNode( segs1[0], 1 )
  36 LastNodeID1   = mesh.GetElemNode( segs1[-1], 1 )
  37 FirstNodeID2  = mesh.GetElemNode( segs2[0], 0 )
  38 SecondNodeID2 = mesh.GetElemNode( segs2[0], 1 )
  39 LastNodeID2   = mesh.GetElemNode( segs2[-1], 1 )
  40 CreatePolygons = True
  41 CreatePolyedrs = False
  42
  43 res = mesh.SewFreeBorders(FirstNodeID1, SecondNodeID1, LastNodeID1,
  44                           FirstNodeID2, SecondNodeID2, LastNodeID2,
  45                           CreatePolygons, CreatePolyedrs )
  46 print(res)
  47 print("nb polygons:", mesh.NbPolygons())