doc/salome/examples/ghs3d_optimization.py

   1 # Copyright (C) 2007-2016  CEA/DEN, EDF R&D, OPEN CASCADE
   2 #
   3 # This library is free software; you can redistribute it and/or
   4 # modify it under the terms of the GNU Lesser General Public
   5 # License as published by the Free Software Foundation; either
   6 # version 2.1 of the License, or (at your option) any later version.
   7 #
   8 # This library is distributed in the hope that it will be useful,
   9 # but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 # Lesser General Public License for more details.
  12 #
  13 # You should have received a copy of the GNU Lesser General Public
  14 # License along with this library; if not, write to the Free Software
  15 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 #
  17 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 #
  19
  20 import salome
  21 salome.salome_init()
  22
  23 from salome.geom import geomBuilder
  24 geompy = geomBuilder.New(salome.myStudy)
  25
  26 from salome.smesh import smeshBuilder
  27 smesh =  smeshBuilder.New(salome.myStudy)
  28
  29 # create a disk
  30 disk = geompy.MakeDiskR(100., 1, theName="disk")
  31
  32 # triangulate the disk
  33 mesh = smesh.Mesh( disk )
  34 cadsurf = mesh.Triangle( smeshBuilder.MG_CADSurf )
  35 cadsurf.SetQuadAllowed( True )
  36 mesh.Compute()
  37
  38 # extrude the 2D mesh into a prismatic mesh
  39 mesh.ExtrusionSweepObject( mesh, [0,0,10], 7 )
  40
  41 # split prisms into tetrahedra
  42 mesh.SplitVolumesIntoTetra( mesh )
  43
  44 # copy the mesh into a new mesh, since only a mesh not based of geometry
  45 # can be optimized using MG-Tetra Optimization
  46 optMesh = smesh.CopyMesh( mesh, "optimization" )
  47
  48 # add MG-Tetra Optimization
  49 mg_opt = optMesh.Tetrahedron( smeshBuilder.MG_Tetra_Optimization )
  50 mg_opt.SetSmoothOffSlivers( True )
  51 mg_opt.SetOptimizationLevel( smeshBuilder.Strong_Optimization )
  52
  53 # run optimization
  54 optMesh.Compute()
  55
  56 print "Nb tetra before optimization", mesh.NbTetras()
  57 print "Nb tetra after  optimization", optMesh.NbTetras()
  58
  59 # End of script
  60