Cut hexahedra

Uniform Cut

To make a uniform cut of an hexahedra in python mode, you need the following arguments:

  • edge : an edge on the hexahedra to cut.

  • nbCuts: the number of cuts.

Use the function cutUni:

elts = doc.cutUni(edge, nbCuts)

GUI command: Uniform Cut

Custom Cut

To make a custom cut of an hexahedra you need the following arguments:

  • edge: an edge on the hexahedra to cut.

  • tl : the list of the lengths of the layers in ascendant order.

Use the funtion cut:

elts = doc.cut(edge, tl)

GUI command: Custom Cut

Operations on elts: Elements

Example

 1# -*- coding: utf-8 -*-
 2# Copyright (C) 2009-2026  CEA, EDF
 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####### Test for Cut Hexa ###############
22
23import hexablock
24
25
26doc = hexablock.addDocument("Cut Hexa Test")
27
28size_x = 2
29size_y = 1
30size_z = 1
31
32grid  = doc.makeCartesianTop (size_x,size_y,size_z)
33arete = grid.getEdgeK (0, 0, 0)
34
35doc.saveVtk ("decoupe1.vtk")
36nbCuts = 2
37grid2 = doc.cutUni (arete, nbCuts)
38doc.saveVtk ("decoupe2.vtk")
39
40#tablen = []
41#reste  = 1
42#abscisse = 0
43#for nl in range(5):
44#    reste     /= 2
45#    abscisse  += reste
46#    tablen.append (abscisse)
47tablen = [2, 3, 8, 14, 18]
48arete = grid2.getEdge(0)
49grid3 = doc.cut(arete, tablen)
50doc.saveVtk ("decoupe3.vtk")