Cylinders and Pipes Construction

Cylinder

Simple Cylinder

To make a simple cylinder grid in python mode, the following data are required:

  • nbR : number of hexahedra on radial.

  • nbA : number of hexahedra along the perimeter of the cylinder.

  • nbH : number of hexahedra along the axis of the cylinder.

Use the function makeCyinderTop:

elts = doc.makeCylinderTop(nbR, nbA, nbH)

GUI command: Simple Cylinder

Uniform Cylinder

The following data are required:

  • origin : origin’s coordinates of the cylinder (select a vertex).

  • vx : the base of the cylinder (select a vector).

  • vz : the axis of the cylinder (select a vector).

  • rint : the radius of the hole in the cylinder .

  • rext : the radius of the cylinder.

  • angle : angle of the cylinder around the Z axis.

  • hauteur: the height of the cylinder.

  • nbR : number of hexahedra on radial.

  • nbA : number of hexahedra along the perimeter of the cylinder.

  • nbH : number of hexahedra along the axis of the cylinder.

Use the function makeCylinderUni to make a uniform cylinder:

elts = doc.makeCylinderUni(origin, vx, vz, rint, rext, angle, hauteur, nbR, nbA, nbH)

GUI command: Uniform Cylinder

Custom Cylinder

The following data are required:

  • origin : origin’s coordinates of the cylinder (select a vertex).

  • vx : the base of the cylinder (select a vector).

  • vz : the axis of the cylinder (select a vector).

  • tr : a list of radiuses in ascendant order.

  • ta : a list of angles in ascendant order.

  • th : a list of heights in ascendant order.

Use the function makeCylinder to make a custom cylinder:

elts = doc.makeCylinder(origin, vx, vz, tr, ta, th)

GUI command: Custom Cylinder

Operations on elts: Elements

Example

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
# -*- coding: utf-8 -*-
# Copyright (C) 2009-2023  CEA, EDF
#
# This library is free software; you can redistribute it and/or
# modify it under the terms of the GNU Lesser General Public
# License as published by the Free Software Foundation; either
# version 2.1 of the License, or (at your option) any later version.
#
# This library is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
# Lesser General Public License for more details.
#
# You should have received a copy of the GNU Lesser General Public
# License along with this library; if not, write to the Free Software
# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
#
# See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
#

####### Make Cylinder Test ########## 

import hexablock


# Simple Cylinder -----

doc  = hexablock.addDocument ("Make Cylinder Test")

nbR = 8
nbA = 10
nbH = 8

cylinder0 = doc.makeCylinderTop (nbR, nbA, nbH)
cylinder0.saveVtk ("makeCylinderTop.vtk")
   
# Uniform Cylinder -----

origin = doc.addVertex (0, 0, 5)
vx     = doc.addVector (1, 0, 0)
vz     = doc.addVector (0, 0, 1)
rint   = 2
rext   = 4
angle  = 300
hauteur = 1


cylinder1 = doc.makeCylinderUni (origin, vx, vz, rint, rext, angle, hauteur, nbR, nbA, nbH)
cylinder1.saveVtk ("makeCylinderUni.vtk")

# Custom Cylinder

origin = doc.addVertex (5, 0, 0)
tr = [10, 20, 30, 40]
ta = [45, 90, 135, 180, 225]
th = [5, 30, 40, 60]

cylinder2 = doc.makeCylinder(origin, vx, vz, tr, ta, th)
cylinder2.saveVtk("makeCylinder.vtk")

Make Cylinders

To make two cylinders in T shape the following data are required for each cylinder:

  • orig : the origin of the cylinder (select a vertex).

  • vz : the axis of the cylinder (select a vector).

  • rext : the radius of the cylinder.

  • h : the height of the cylinder.

One of the two cylinders must be bigger than the other.

Use the function makeCylinders:

elts = doc.makeCylinders (orig1, vz1, rext1, h1, orig2, vz2, rext2, h2)

GUI command: Cylinders

Operations on elts: Elements

Example

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
# -*- coding: utf-8 -*-
# Copyright (C) 2009-2023  CEA, EDF
#
# This library is free software; you can redistribute it and/or
# modify it under the terms of the GNU Lesser General Public
# License as published by the Free Software Foundation; either
# version 2.1 of the License, or (at your option) any later version.
#
# This library is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
# Lesser General Public License for more details.
#
# You should have received a copy of the GNU Lesser General Public
# License along with this library; if not, write to the Free Software
# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
#
# See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
#

####### Make Cylinders Test ########## 

import hexablock
import math

    
doc  = hexablock.addDocument ("Make Cylinders Test")

rext1 =  2*math.sqrt (2.0)
rext2 =  3*math.sqrt (2.0)
h1   =  16
h2   =  16
xl1  =  -8

orig1  = doc.addVertex ( 0, 0,  xl1)
orig2  = doc.addVertex (-8, 0,  0)
vz1   = doc.addVector ( 0, 0,  1)
vz2   = doc.addVector ( 1, 0,  0)
   
cylinders = doc.makeCylinders (orig1, vz1, rext1, h1, orig2, vz2, rext2, h2)
cylinders.saveVtk("makeCylinders.vtk")

Pipe

Simple Pipe

To make a simple pipe grid in python mode, the following data are required:

  • nbR : number of hexahedra on radial.

  • nbA : number of hexahedra along the perimeter of the pipe.

  • nbH : number of hexahedra along the axis of the pipe.

Use the function makePipeTop:

elts = doc.makePipeTop(nbR, nbA, nbH)

GUI command: Simple Pipe

Uniform Pipe

The following data are required:

  • origin : the origin’s coordinates of the pipe (select a vertex).

  • vx : the base of the pipe (select a vector).

  • vz : the axis of the pipe (select a vector).

  • rint : the radius of the hole in the pipe .

  • rext : the radius of the pipe.

  • angle : angle of the pipe around the Z axis.

  • hauteur: the height of the pipe.

  • nbR : number of hexahedra on radial.

  • nbA : number of hexahedra along the perimeter of the pipe.

  • nbH : number of hexahedra along the axis of the pipe.

Use the function makePipeUni to make a uniform pipe:

elts = doc.makePipeUni(origin, vx, vz, rint, rext, angle, hauteur, nbR, nbA, nbH)

GUI command: Uniform Pipe

Custom Pipe

The following data are required:

  • origin : origin’s coordinates of the pipe (select a vertex).

  • vx : the base of the pipe (select a vector).

  • vz : the axis of the pipe (select a vector).

  • tr : a list of radiuses in ascendant order.

  • ta : a list of angles in ascendant order.

  • th : a list of heights in ascendant order.

Use the function makePipe to make a custom pipe:

elts = doc.makePipe(origin, vx, vz, tr, ta, th)

GUI command: Custom Pipe

Operations on elts: Elements

Example

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
# -*- coding: utf-8 -*-
# Copyright (C) 2009-2023  CEA, EDF
#
# This library is free software; you can redistribute it and/or
# modify it under the terms of the GNU Lesser General Public
# License as published by the Free Software Foundation; either
# version 2.1 of the License, or (at your option) any later version.
#
# This library is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
# Lesser General Public License for more details.
#
# You should have received a copy of the GNU Lesser General Public
# License along with this library; if not, write to the Free Software
# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
#
# See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
#

####### Make Pipe Test ########## 

import hexablock
        
doc  = hexablock.addDocument ("Make Pipe Test")

# Simple Pipe -----

nbR = 8
nbA = 10
nbH = 8

pipe0 = doc.makePipeTop (nbR, nbA, nbH)
pipe0.saveVtk ("makePipeTop.vtk")

# Uniform Pipe -----

orig = doc.addVertex (0, 0, 0)
vx   = doc.addVector (1, 0, 0)
vz   = doc.addVector (0, 0, 1)
rint = 1
rext = 3
angle = 360
hauteur = 2

pipe1 = doc.makePipeUni (orig, vx, vz, rint, rext, angle, hauteur, nbR, nbA, nbH)
pipe1.saveVtk ("makePipeUni.vtk");
   

# Custom Pipe -----

origin = doc.addVertex (0, 5, 0)
tr = [10, 20, 30, 40]
ta = [45, 90, 135, 180, 225]
th = [5, 30, 40, 60]

pipe2 = doc.makePipe(origin, vx, vz, tr, ta, th)
pipe2.saveVtk("makePipe.vtk")

Make Pipes

To make two pipes in T shape the following data are required for each pipe:

  • orig : the origin of the pipe (select a vertex).

  • vz : the axis of the pipe (select a vector).

  • rint : the internal radius of the pipe.

  • rext : the radius of the pipe.

  • h : the height of the pipe.

One of the two pipes must be bigger than the other.

Use the function makepipes:

elts = doc.makePipes(orig1, vz1, rint1, rext1, h1, orig2, vz2, rint2, rext2, h2)

GUI command: Pipes

Operations on elts: Elements

Example

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
# -*- coding: utf-8 -*-
# Copyright (C) 2009-2023  CEA, EDF
#
# This library is free software; you can redistribute it and/or
# modify it under the terms of the GNU Lesser General Public
# License as published by the Free Software Foundation; either
# version 2.1 of the License, or (at your option) any later version.
#
# This library is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
# Lesser General Public License for more details.
#
# You should have received a copy of the GNU Lesser General Public
# License along with this library; if not, write to the Free Software
# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
#
# See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
#

####### Make Pipes Test ########## 

import hexablock
import math

    
doc  = hexablock.addDocument ("Make Pipes Test")

rext1 =  2*math.sqrt (2.0)
rext2 =  3*math.sqrt (2.0)
rint1 =  rext1/2
rint2 =  rext1
h1   =  16
h2   =  16
xl1  =  -8

orig1  = doc.addVertex ( 0, 0,  xl1)
orig2  = doc.addVertex (-8, 0,  0)
vz1   = doc.addVector ( 0, 0,  1)
vz2   = doc.addVector ( 1, 0,  0)
   
pipes = doc.makePipes (orig1, vz1, rint1, rext1, h1, orig2, vz2, rint2, rext2, h2);
pipes.saveVtk("makePipes.vtk")