X-Git-Url: http://git.salome-platform.org/gitweb/?p=modules%2Fsmesh.git;a=blobdiff_plain;f=src%2FSMESH_SWIG%2Fex19_sphereINcube.py;h=4528d532324e46d8dde196d169ce97fcb1791e20;hp=19e9f84717e81b7dc15f65f70b049c110c6ee2be;hb=ccb5e3c25bf204e839c0094f50f543eacedba5e6;hpb=57b43b4d010e2d0a1529d3c131bbb9d416e63258 diff --git a/src/SMESH_SWIG/ex19_sphereINcube.py b/src/SMESH_SWIG/ex19_sphereINcube.py index 19e9f8471..4528d5323 100644 --- a/src/SMESH_SWIG/ex19_sphereINcube.py +++ b/src/SMESH_SWIG/ex19_sphereINcube.py @@ -1,28 +1,37 @@ -# Copyright (C) 2005 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN, -# CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS +# -*- coding: iso-8859-1 -*- +# Copyright (C) 2007-2016 CEA/DEN, EDF R&D, OPEN CASCADE # -# 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. +# Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN, +# CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS # -# 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. +# 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. # -# 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 +# 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. # -# See http://www.salome-platform.org/ +# 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 # -# CEA/LGLS 2005, Francis KLOSS (OCC) -# ================================== -from geompy import * +# ================================== +# +import salome +salome.salome_init() +import GEOM +from salome.geom import geomBuilder +geompy = geomBuilder.New(salome.myStudy) -import smesh +import SMESH, SALOMEDS +from salome.smesh import smeshBuilder +smesh = smeshBuilder.New(salome.myStudy) # Geometrie # ========= @@ -41,76 +50,98 @@ plan_trim = 1000 # Sphere # ------ -sphere_centre = MakeVertex(0, 0, 0) +sphere_centre = geompy.MakeVertex(0, 0, 0) -sphere_pleine = MakeSpherePntR(sphere_centre, sphere_rayon) +sphere_pleine = geompy.MakeSpherePntR(sphere_centre, sphere_rayon) # Cube interieur # -------------- boite_cote = sphere_rayon / 2 -boite = MakeBox(-boite_cote, -boite_cote, -boite_cote, +boite_cote, +boite_cote, +boite_cote) +boite = geompy.MakeBox(-boite_cote, -boite_cote, -boite_cote, +boite_cote, +boite_cote, +boite_cote) blocs = [boite] # Decoupage sphere # ---------------- -sphere_troue = MakeCut(sphere_pleine, boite) +sphere_troue = geompy.MakeCut(sphere_pleine, boite) + +#sphere_outils = [] +#sphere_outils.append(MakePlane(sphere_centre, MakeVectorDXDYDZ( 1, 0, 1), plan_trim)) +#sphere_outils.append(MakePlane(sphere_centre, MakeVectorDXDYDZ( 1, 0, -1), plan_trim)) +#sphere_outils.append(MakePlane(sphere_centre, MakeVectorDXDYDZ( 1, 1, 0), plan_trim)) +#sphere_outils.append(MakePlane(sphere_centre, MakeVectorDXDYDZ(-1, 1, 0), plan_trim)) + +f1 = geompy.MakePlane(sphere_centre, geompy.MakeVectorDXDYDZ( 1, 0, 1), plan_trim) +f2 = geompy.MakePlane(sphere_centre, geompy.MakeVectorDXDYDZ(-1, 1, 0), plan_trim) +f3 = geompy.MakePlane(sphere_centre, geompy.MakeVectorDXDYDZ( 1, 1, 0), plan_trim) +f4 = geompy.MakePlane(sphere_centre, geompy.MakeVectorDXDYDZ( 1, 0, -1), plan_trim) + -sphere_outils = [] -sphere_outils.append(MakePlane(sphere_centre, MakeVectorDXDYDZ( 1, 0, 1), plan_trim)) -sphere_outils.append(MakePlane(sphere_centre, MakeVectorDXDYDZ( 1, 0, -1), plan_trim)) -sphere_outils.append(MakePlane(sphere_centre, MakeVectorDXDYDZ( 1, 1, 0), plan_trim)) -sphere_outils.append(MakePlane(sphere_centre, MakeVectorDXDYDZ(-1, 1, 0), plan_trim)) +#sphere_decoupee = MakePartition(solids, sphere_outils, [], [], ShapeType["SOLID"]) -sphere_decoupee = MakePartition([sphere_troue], sphere_outils, [], [], ShapeType["SOLID"]) +sphere_decoupee = geompy.MakePartition([sphere_troue], [f1], [], [], geompy.ShapeType["SOLID"]) +sphere_decoupee = geompy.MakePartition([sphere_decoupee], [f2], [], [], geompy.ShapeType["SOLID"]) +sphere_decoupee = geompy.MakePartition([sphere_decoupee], [f3], [], [], geompy.ShapeType["SOLID"]) +sphere_decoupee = geompy.MakePartition([sphere_decoupee], [f4], [], [], geompy.ShapeType["SOLID"]) -sphere_partie = GetBlockNearPoint(sphere_decoupee, MakeVertex(-sphere_rayon, 0, 0)) -sphere_bloc = RemoveExtraEdges(sphere_partie) +sphere_partie = geompy.MakeCompound([sphere_decoupee]) + +sphere_partie = geompy.GetBlockNearPoint(sphere_decoupee, geompy.MakeVertex(-sphere_rayon, 0, 0)) +sphere_bloc = geompy.RemoveExtraEdges(sphere_partie) blocs.append(sphere_bloc) pi2 = 3.141592653/2 -sphere_dir1 = MakeVectorDXDYDZ(0, 1, 0) -sphere_dir2 = MakeVectorDXDYDZ(0, 0, 1) +sphere_dir1 = geompy.MakeVectorDXDYDZ(0, 1, 0) +sphere_dir2 = geompy.MakeVectorDXDYDZ(0, 0, 1) -blocs.append(MakeRotation(sphere_bloc, sphere_dir1, +pi2)) -blocs.append(MakeRotation(sphere_bloc, sphere_dir1, -pi2)) +blocs.append(geompy.MakeRotation(sphere_bloc, sphere_dir1, +pi2)) +blocs.append(geompy.MakeRotation(sphere_bloc, sphere_dir1, -pi2)) -blocs.append(MakeRotation(sphere_bloc, sphere_dir2, +pi2)) -blocs.append(MakeRotation(sphere_bloc, sphere_dir2, -pi2)) +blocs.append(geompy.MakeRotation(sphere_bloc, sphere_dir2, +pi2)) +blocs.append(geompy.MakeRotation(sphere_bloc, sphere_dir2, -pi2)) -blocs.append(MakeMirrorByPoint(sphere_bloc, sphere_centre)) +blocs.append(geompy.MakeMirrorByPoint(sphere_bloc, sphere_centre)) # Cube exterieur # -------------- -cube_plein = MakeBox(-cube_cote, -cube_cote, -cube_cote, +cube_cote, +cube_cote, +cube_cote) -cube_trou = MakeCut(cube_plein, sphere_pleine) -cube_decoupe = MakePartition([cube_trou], sphere_outils, [], [], ShapeType["SOLID"]) -cube_partie = GetBlockNearPoint(cube_decoupe, MakeVertex(-cube_cote, 0, 0)) -cube_bloc = RemoveExtraEdges(cube_partie) +cube_plein = geompy.MakeBox(-cube_cote, -cube_cote, -cube_cote, +cube_cote, +cube_cote, +cube_cote) +cube_trou = geompy.MakeCut(cube_plein, sphere_pleine) +#cube_decoupe = MakePartition([cube_trou], sphere_outils, [], [], ShapeType["SOLID"]) + +cube_decoupe = geompy.MakePartition([cube_trou], [f1], [], [], geompy.ShapeType["SOLID"]) +cube_decoupe = geompy.MakePartition([cube_decoupe], [f2], [], [], geompy.ShapeType["SOLID"]) +cube_decoupe = geompy.MakePartition([cube_decoupe], [f3], [], [], geompy.ShapeType["SOLID"]) +cube_decoupe = geompy.MakePartition([cube_decoupe], [f4], [], [], geompy.ShapeType["SOLID"]) + +cube_decoupe = geompy.MakeCompound([cube_decoupe]) + + +cube_partie = geompy.GetBlockNearPoint(cube_decoupe, geompy.MakeVertex(-cube_cote, 0, 0)) +cube_bloc = geompy.RemoveExtraEdges(cube_partie) blocs.append(cube_bloc) -blocs.append(MakeRotation(cube_bloc, sphere_dir1, +pi2)) -blocs.append(MakeRotation(cube_bloc, sphere_dir1, -pi2)) +blocs.append(geompy.MakeRotation(cube_bloc, sphere_dir1, +pi2)) +blocs.append(geompy.MakeRotation(cube_bloc, sphere_dir1, -pi2)) -blocs.append(MakeRotation(cube_bloc, sphere_dir2, +pi2)) -blocs.append(MakeRotation(cube_bloc, sphere_dir2, -pi2)) +blocs.append(geompy.MakeRotation(cube_bloc, sphere_dir2, +pi2)) +blocs.append(geompy.MakeRotation(cube_bloc, sphere_dir2, -pi2)) -blocs.append(MakeMirrorByPoint(cube_bloc, sphere_centre)) +blocs.append(geompy.MakeMirrorByPoint(cube_bloc, sphere_centre)) # Piece # ----- -piece_cpd = MakeCompound(blocs) -piece = MakeGlueFaces(piece_cpd, 1.e-3) +piece_cpd = geompy.MakeCompound(blocs) +piece = geompy.MakeGlueFaces(piece_cpd, 1.e-3) -piece_id = addToStudy(piece, "ex19_sphereINcube") +piece_id = geompy.addToStudy(piece, "ex19_sphereINcube") # Groupe geometrique # ================== @@ -118,18 +149,18 @@ piece_id = addToStudy(piece, "ex19_sphereINcube") # Definition du groupe # -------------------- -groupe = CreateGroup(piece, ShapeType["SOLID"]) +groupe = geompy.CreateGroup(piece, geompy.ShapeType["SOLID"]) groupe_nom = "ex19_sphereINcube_interieur" -addToStudy(groupe, groupe_nom) +geompy.addToStudy(groupe, groupe_nom) groupe.SetName(groupe_nom) # Contenu du groupe # ----------------- -groupe_sphere = GetShapesOnSphere(piece, ShapeType["SOLID"], sphere_centre, sphere_rayon, GEOM.ST_ONIN) +groupe_sphere = geompy.GetShapesOnSphere(piece, geompy.ShapeType["SOLID"], sphere_centre, sphere_rayon, GEOM.ST_ONIN) -UnionList(groupe, groupe_sphere) +geompy.UnionList(groupe, groupe_sphere) # Meshing # ======= @@ -155,3 +186,8 @@ hexa.Compute() # -------------------- hexa_groupe = hexa.Group(groupe) + +# Update object browser +# --------------------- + +salome.sg.updateObjBrowser(True)