doc/salome/examples/modifying_meshes_ex18.py

   1 # Uniting a Set of Triangles
   2
   3 import salome
   4 import smesh
   5
   6 # create an empty mesh structure
   7 mesh = smesh.Mesh()
   8
   9 # create the following mesh:
  10 # .----.----.----.
  11 # |   /|   /|   /|
  12 # |  / |  / |  / |
  13 # | /  | /  | /  |
  14 # |/   |/   |/   |
  15 # .----.----.----.
  16
  17 bb = [0, 0, 0, 0]
  18 tt = [0, 0, 0, 0]
  19 ff = [0, 0, 0, 0, 0, 0]
  20
  21 bb[0] = mesh.AddNode( 0., 0., 0.)
  22 bb[1] = mesh.AddNode(10., 0., 0.)
  23 bb[2] = mesh.AddNode(20., 0., 0.)
  24 bb[3] = mesh.AddNode(30., 0., 0.)
  25
  26 tt[0] = mesh.AddNode( 0., 15., 0.)
  27 tt[1] = mesh.AddNode(10., 15., 0.)
  28 tt[2] = mesh.AddNode(20., 15., 0.)
  29 tt[3] = mesh.AddNode(30., 15., 0.)
  30
  31 ff[0] = mesh.AddFace([bb[0], bb[1], tt[1]])
  32 ff[1] = mesh.AddFace([bb[0], tt[1], tt[0]])
  33 ff[2] = mesh.AddFace([bb[1], bb[2], tt[2]])
  34 ff[3] = mesh.AddFace([bb[1], tt[2], tt[1]])
  35 ff[4] = mesh.AddFace([bb[2], bb[3], tt[3]])
  36 ff[5] = mesh.AddFace([bb[2], tt[3], tt[2]])
  37
  38 # unite a set of triangles
  39 print "\nUnite a set of triangles ... ",
  40 res = mesh.TriToQuad([ff[2], ff[3], ff[4], ff[5]], smesh.FT_MinimumAngle, 60.)
  41 if not res: print "failed!"
  42 else:       print "done."
  43
  44 salome.sg.updateObjBrowser(1)