src/SMESH_SWIG/ex10_grid4geometry.py

   1 #  -*- coding: iso-8859-1 -*-
   2 # Copyright (C) 2007-2021  CEA/DEN, EDF R&D, OPEN CASCADE
   3 #
   4 # Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   5 # CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   6 #
   7 # This library is free software; you can redistribute it and/or
   8 # modify it under the terms of the GNU Lesser General Public
   9 # License as published by the Free Software Foundation; either
  10 # version 2.1 of the License, or (at your option) any later version.
  11 #
  12 # This library is distributed in the hope that it will be useful,
  13 # but WITHOUT ANY WARRANTY; without even the implied warranty of
  14 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  15 # Lesser General Public License for more details.
  16 #
  17 # You should have received a copy of the GNU Lesser General Public
  18 # License along with this library; if not, write to the Free Software
  19 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  20 #
  21 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  22 #
  23
  24 # =======================================
  25 #
  26 import salome
  27 salome.salome_init()
  28 import GEOM
  29 from salome.geom import geomBuilder
  30 geompy = geomBuilder.New()
  31
  32 import SMESH, SALOMEDS
  33 from salome.smesh import smeshBuilder
  34 smesh =  smeshBuilder.New()
  35
  36 # Geometry
  37 # ========
  38
  39 # Element of a grid compound by a square with a cylinder on each vertex build by using partition
  40
  41 # Values
  42 # ------
  43
  44 ox = 0
  45 oy = 0
  46 oz = 0
  47
  48 arete   =  50
  49 hauteur = 100
  50 rayon   =  10
  51
  52 g_trim = 1000
  53
  54 # Box
  55 # ---
  56
  57 piecePoint = geompy.MakeVertex(ox, oy, oz)
  58
  59 pieceBox = geompy.MakeBoxTwoPnt(piecePoint, geompy.MakeVertex(ox+arete, oy+hauteur, oz+arete))
  60
  61 # Cut by cylinders
  62 # ----------------
  63
  64 dirUp = geompy.MakeVectorDXDYDZ(0, 1, 0)
  65
  66 pieceCut1 = geompy.MakeCut(pieceBox , geompy.MakeCylinder(piecePoint                        , dirUp, rayon, hauteur))
  67 pieceCut2 = geompy.MakeCut(pieceCut1, geompy.MakeCylinder(geompy.MakeVertex(ox+arete, oy, oz      ), dirUp, rayon, hauteur))
  68 pieceCut3 = geompy.MakeCut(pieceCut2, geompy.MakeCylinder(geompy.MakeVertex(ox      , oy, oz+arete), dirUp, rayon, hauteur))
  69 pieceCut4 = geompy.MakeCut(pieceCut3, geompy.MakeCylinder(geompy.MakeVertex(ox+arete, oy, oz+arete), dirUp, rayon, hauteur))
  70
  71 # Compound by make a partition of a solid
  72 # ---------------------------------------
  73
  74 dir = geompy.MakeVectorDXDYDZ(-1, 0, 1)
  75
  76 tools = []
  77 tools.append(geompy.MakePlane(geompy.MakeVertex(ox+rayon, oy, oz      ), dir, g_trim))
  78 tools.append(geompy.MakePlane(geompy.MakeVertex(ox      , oy, oz+rayon), dir, g_trim))
  79
  80 piece = geompy.MakePartition([pieceCut4], tools, [], [], geompy.ShapeType["SOLID"])
  81
  82 # Add in study
  83 # ------------
  84
  85 piece_id = geompy.addToStudy(piece, "ex10_grid4geometry")
  86
  87 # Meshing
  88 # =======
  89
  90 # Create a hexahedral mesh
  91 # ------------------------
  92
  93 hexa = smesh.Mesh(piece, "ex10_grid4geometry:hexa")
  94
  95 algo = hexa.Segment()
  96 algo.NumberOfSegments(10)
  97
  98 hexa.Quadrangle()
  99
 100 hexa.Hexahedron()
 101
 102 # Mesh calculus
 103 # -------------
 104
 105 hexa.Compute()
 106
 107 # Update object browser
 108 # ---------------------
 109
 110 salome.sg.updateObjBrowser()