doc/salome/examples/extrusion_penta_biquad.py

   1 # -*- coding: utf-8 -*-
   2
   3 import sys
   4 import salome
   5
   6 salome.salome_init()
   7 theStudy = salome.myStudy
   8
   9 import GEOM
  10 from salome.geom import geomBuilder
  11 import math
  12 import SALOMEDS
  13
  14 geompy = geomBuilder.New(theStudy)
  15
  16 O = geompy.MakeVertex(0, 0, 0)
  17 OX = geompy.MakeVectorDXDYDZ(1, 0, 0)
  18 OY = geompy.MakeVectorDXDYDZ(0, 1, 0)
  19 OZ = geompy.MakeVectorDXDYDZ(0, 0, 1)
  20 Divided_Disk_1 = geompy.MakeDividedDisk(100, 1, GEOM.SQUARE)
  21 geompy.addToStudy( O, 'O' )
  22 geompy.addToStudy( OX, 'OX' )
  23 geompy.addToStudy( OY, 'OY' )
  24 geompy.addToStudy( OZ, 'OZ' )
  25 geompy.addToStudy( Divided_Disk_1, 'Divided Disk_1' )
  26
  27 import  SMESH, SALOMEDS
  28 from salome.smesh import smeshBuilder
  29
  30 smesh = smeshBuilder.New(theStudy)
  31 aFilterManager = smesh.CreateFilterManager()
  32 Mesh_1 = smesh.Mesh(Divided_Disk_1)
  33 Regular_1D = Mesh_1.Segment()
  34 Number_of_Segments_1 = Regular_1D.NumberOfSegments(6)
  35 Quadrangle_2D = Mesh_1.Quadrangle(algo=smeshBuilder.QUADRANGLE)
  36 isDone = Mesh_1.Compute()
  37 aMaxElementLength2D0x5d7fdf0 = aFilterManager.CreateMaxElementLength2D()
  38 isDone = Mesh_1.QuadToTriObject( Mesh_1,  )
  39 Mesh_1.ExtrusionSweepObjects( [ Mesh_1 ], [ Mesh_1 ], [ Mesh_1 ], [ 0, 0, 50 ], 3, 1 )
  40 Mesh_1.ConvertToQuadratic(0, Mesh_1,True)
  41
  42 infos = Mesh_1.GetMeshInfo()
  43 print "Number of biquadratic pentahedrons:", infos[SMESH.Entity_BiQuad_Penta]
  44 if (infos[SMESH.Entity_BiQuad_Penta] != 1080):
  45   raise RuntimeError("Bad number of biquadratic pentahedrons: should be 1080")
  46
  47 ## Set names of Mesh objects
  48 smesh.SetName(Regular_1D.GetAlgorithm(), 'Regular_1D')
  49 smesh.SetName(Quadrangle_2D.GetAlgorithm(), 'Quadrangle_2D')
  50 smesh.SetName(Number_of_Segments_1, 'Number of Segments_1')
  51 smesh.SetName(Mesh_1.GetMesh(), 'Mesh_1')
  52
  53
  54 if salome.sg.hasDesktop():
  55   salome.sg.updateObjBrowser(True)