]> SALOME platform Git repositories - modules/smesh.git/blob - src/SMESH_SWIG/ex30_tepal.py
Salome HOME
Nommage des maillags intermédiaires
[modules/smesh.git] / src / SMESH_SWIG / ex30_tepal.py
1 #  -*- coding: iso-8859-1 -*-
2 # Copyright (C) 2007-2021  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()