1 # Copyright (C) 2007-2008 CEA/DEN, EDF R&D, OPEN CASCADE
3 # Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 # CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
6 # This library is free software; you can redistribute it and/or
7 # modify it under the terms of the GNU Lesser General Public
8 # License as published by the Free Software Foundation; either
9 # version 2.1 of the License.
11 # This library is distributed in the hope that it will be useful,
12 # but WITHOUT ANY WARRANTY; without even the implied warranty of
13 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 # Lesser General Public License for more details.
16 # You should have received a copy of the GNU Lesser General Public
17 # License along with this library; if not, write to the Free Software
18 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
20 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
22 # =======================================
31 # Creer un cylindre surplombe d'une demi-sphere le tout troue par un petit cylindre.
32 # Decouper en hexahedre et mailler.
38 cylindre_hauteur = 400
41 trou_z = cylindre_rayon/2
48 cylindre_base = MakeVertex(0, 0, 0)
49 cylindre_dir = MakeVectorDXDYDZ(1, 0, 0)
50 cylindre = MakeCylinder(cylindre_base, cylindre_dir, cylindre_rayon, cylindre_hauteur)
55 dome_sphere = MakeSpherePntR(cylindre_base, cylindre_rayon)
56 dome = MakeFuse(dome_sphere, cylindre)
61 cheminee_base = MakeVertex(-cylindre_hauteur/2, 0, trou_z)
62 cheminee_trou = MakeCylinder(cheminee_base, cylindre_dir, trou_rayon, 2*cylindre_hauteur)
63 cheminee = MakeCut(dome, cheminee_trou)
65 # Decoupage et reparation
66 # -----------------------
68 blocs_plan1 = MakePlane(cheminee_base, MakeVectorDXDYDZ(0, 1, 0), plan_trim)
69 blocs_plan2 = MakePlane(cheminee_base, MakeVectorDXDYDZ(0, 0, 1), plan_trim)
71 blocs_part = MakePartition([cheminee], [blocs_plan1, blocs_plan2], [], [], ShapeType["SOLID"])
73 piece = RemoveExtraEdges(blocs_part)
78 piece_id = addToStudy(piece, "ex17_dome1")
83 # Maillage hexahedrique
84 # ---------------------
86 hexa = smesh.Mesh(piece, "ex17_dome1:hexa")
89 algo.NumberOfSegments(20)