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 # GEOM GEOM_SWIG : binding of C++ omplementaion with Python
23 # File : GEOM_Sphere.py
24 # Author : Damien COQUERET, Open CASCADE
33 # It is an example of creating a hexahedrical mesh on a sphere.
35 # Used approach allows to avoid problems with degenerated and
36 # seam edges without special processing of geometrical shapes
38 #-----------------------------------------------------------------------
50 P0 = MakeVertex(0., 0., 0.)
51 P1 = MakeVertex(-Dist, -Dist, -Dist)
52 P2 = MakeVertex(-Dist, -Dist, Dist)
53 P3 = MakeVertex(-Dist, Dist, Dist)
54 P4 = MakeVertex(-Dist, Dist, -Dist)
56 VZ = MakeVectorDXDYDZ(0., 0., 1.)
58 #Construction Elements
65 PolyLine = MakePolyline(PointsList)
67 Face1 = MakeFace(PolyLine, 1)
68 Face2 = MakeScaleTransform(Face1, P0, Factor)
69 Face3 = MakeScaleTransform(Face1, P0, -1.)
72 Sphere = MakeSphereR(Radius)
74 Block = MakeHexa2Faces(Face1, Face2)
75 Cube = MakeHexa2Faces(Face1, Face3)
77 Common1 = MakeBoolean(Sphere, Block, 1)
78 Common2 = MakeRotation(Common1, VZ, Angle90)
80 MultiBlock1 = MakeMultiTransformation1D(Common1, 21, -1, 3)
81 MultiBlock2 = MakeMultiTransformation1D(Common2, 31, -1, 3)
83 #Reconstruct sphere from several blocks
84 ShapesList.append(Cube)
85 ShapesList.append(MultiBlock1)
86 ShapesList.append(MultiBlock2)
87 Compound = MakeCompound(ShapesList)
89 Result = MakeGlueFaces(Compound, 0.1)
92 Id_Sphere = addToStudy(Sphere, "Sphere")
93 Id_Cube = addToStudy(Cube, "Cube")
95 Id_Common1 = addToStudy(Common1, "Common1")
96 Id_Common2 = addToStudy(Common2, "Common2")
98 Id_MultiBlock1 = addToStudy(MultiBlock1, "MultiBlock1")
99 Id_MultiBlock2 = addToStudy(MultiBlock2, "MultiBlock2")
101 Id_Result = addToStudy(Result, "Result")
103 #-----------------------------------------------------------------------
105 my_hexa = smesh.Mesh(Result, "Sphere_Mesh")
106 algo = my_hexa.Segment()
107 algo.NumberOfSegments(NbSeg)