src/SMESH_SWIG/ex30_tepal.py

   1 #  -*- coding: iso-8859-1 -*-
   2 # Copyright (C) 2007-2022  CEA/DEN, EDF R&D, OPEN CASCADE
   3 #
   4 # This library is free software; you can redistribute it and/or
   5 # modify it under the terms of the GNU Lesser General Public
   6 # License as published by the Free Software Foundation; either
   7 # version 2.1 of the License, or (at your option) any later version.
   8 #
   9 # This library is distributed in the hope that it will be useful,
  10 # but WITHOUT ANY WARRANTY; without even the implied warranty of
  11 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  12 # Lesser General Public License for more details.
  13 #
  14 # You should have received a copy of the GNU Lesser General Public
  15 # License along with this library; if not, write to the Free Software
  16 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  17 #
  18 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  19 #
  20
  21 # ====================================================================
  22 #
  23 import os
  24
  25 import salome
  26 salome.salome_init()
  27 import GEOM
  28 from salome.geom import geomBuilder
  29 geompy = geomBuilder.New()
  30
  31 import SMESH, SALOMEDS
  32 from salome.smesh import smeshBuilder
  33 smesh =  smeshBuilder.New()
  34
  35 # Parameters
  36 # ----------
  37
  38 results = "/tmp/ZZ"
  39
  40 radius =  50
  41 height = 200
  42
  43 # Build a cylinder
  44 # ----------------
  45
  46 base = geompy.MakeVertex(0, 0, 0)
  47 direction = geompy.MakeVectorDXDYDZ(0, 0, 1)
  48
  49 cylinder = geompy.MakeCylinder(base, direction, radius, height)
  50
  51 geompy.addToStudy(cylinder, "Cylinder")
  52
  53 # Define a mesh on a geometry
  54 # ---------------------------
  55
  56 smesh.UpdateStudy()
  57
  58 m = smesh.Mesh(cylinder)
  59
  60 # 2D mesh with BLSURF
  61 # -------------------
  62
  63 algo2d = m.Triangle(smeshBuilder.BLSURF)
  64
  65 algo2d.SetPhysicalMesh(1)
  66 algo2d.SetPhySize(5)
  67
  68 algo2d.SetGeometricMesh(0)
  69
  70 # 3D mesh with tepal
  71 # ------------------
  72
  73 algo3d = m.Tetrahedron(smeshBuilder.GHS3DPRL)
  74
  75 algo3d.SetMEDName(results)
  76 algo3d.SetNbPart(4)
  77 algo3d.SetBackground(False)
  78 algo3d.SetKeepFiles(False)
  79
  80 # Launch meshers
  81 # --------------
  82
  83 status = m.Compute()
  84
  85 # Test if ok
  86 # ----------
  87
  88 if os.access(results+".xml", os.F_OK):
  89     print("Ok: tepal")
  90 else:
  91     print("KO: tepal")
  92
  93 # Update object browser
  94 # ---------------------
  95
  96 salome.sg.updateObjBrowser()