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=ec5f5d3980ce2082a7bc23da15bdf32f6f7aadbb;hb=ccb5e3c25bf204e839c0094f50f543eacedba5e6;hpb=79b1ac2b6df9117f16f11d444b1f165d477a1813 diff --git a/src/SMESH_SWIG/ex19_sphereINcube.py b/src/SMESH_SWIG/ex19_sphereINcube.py index ec5f5d398..4528d5323 100644 --- a/src/SMESH_SWIG/ex19_sphereINcube.py +++ b/src/SMESH_SWIG/ex19_sphereINcube.py @@ -1,29 +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. +# +# 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 geompy +import SMESH, SALOMEDS +from salome.smesh import smeshBuilder +smesh = smeshBuilder.New(salome.myStudy) # Geometrie # ========= @@ -42,23 +50,23 @@ 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)) @@ -66,83 +74,74 @@ sphere_troue = MakeCut(sphere_pleine, boite) #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 = MakePlane(sphere_centre, MakeVectorDXDYDZ( 1, 0, 1), plan_trim) -f2 = MakePlane(sphere_centre, MakeVectorDXDYDZ(-1, 1, 0), plan_trim) -f3 = MakePlane(sphere_centre, MakeVectorDXDYDZ( 1, 1, 0), plan_trim) -f4 = MakePlane(sphere_centre, MakeVectorDXDYDZ( 1, 0, -1), 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_decoupee = MakePartition(solids, sphere_outils, [], [], ShapeType["SOLID"]) -solids = geompy.SubShapeAll(sphere_troue,geompy.ShapeType["SOLID"]) -sphere_decoupee = MakePartition(solids, [f1], [], [], ShapeType["SOLID"]) -solids = geompy.SubShapeAll(sphere_decoupee,geompy.ShapeType["SOLID"]) -sphere_decoupee = MakePartition(solids, [f2], [], [], ShapeType["SOLID"]) -solids = geompy.SubShapeAll(sphere_decoupee,geompy.ShapeType["SOLID"]) -sphere_decoupee = MakePartition(solids, [f3], [], [], ShapeType["SOLID"]) -solids = geompy.SubShapeAll(sphere_decoupee,geompy.ShapeType["SOLID"]) -sphere_decoupee = MakePartition(solids, [f4], [], [], ShapeType["SOLID"]) -solids = geompy.SubShapeAll(sphere_decoupee,geompy.ShapeType["SOLID"]) -sphere_partie = geompy.MakeCompound(solids) +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_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"]) -solids = geompy.SubShapeAll(cube_trou,geompy.ShapeType["SOLID"]) -cube_decoupe = MakePartition(solids, [f1], [], [], ShapeType["SOLID"]) -solids = geompy.SubShapeAll(cube_decoupe,geompy.ShapeType["SOLID"]) -cube_decoupe = MakePartition(solids, [f2], [], [], ShapeType["SOLID"]) -solids = geompy.SubShapeAll(cube_decoupe,geompy.ShapeType["SOLID"]) -cube_decoupe = MakePartition(solids, [f3], [], [], ShapeType["SOLID"]) -solids = geompy.SubShapeAll(cube_decoupe,geompy.ShapeType["SOLID"]) -cube_decoupe = MakePartition(solids, [f4], [], [], ShapeType["SOLID"]) -solids = geompy.SubShapeAll(cube_decoupe,geompy.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(solids) +cube_decoupe = geompy.MakeCompound([cube_decoupe]) -cube_partie = GetBlockNearPoint(cube_decoupe, MakeVertex(-cube_cote, 0, 0)) -cube_bloc = RemoveExtraEdges(cube_partie) +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 # ================== @@ -150,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 # ======= @@ -187,3 +186,8 @@ hexa.Compute() # -------------------- hexa_groupe = hexa.Group(groupe) + +# Update object browser +# --------------------- + +salome.sg.updateObjBrowser(True)