X-Git-Url: http://git.salome-platform.org/gitweb/?p=modules%2Fsmesh.git;a=blobdiff_plain;f=src%2FSMESH_SWIG%2Fex19_sphereINcube.py;h=d59b048a56b6e148114efdf180765d0c762db363;hp=c303c7072def0901d50c4eb8486cf9bb8c8ebb25;hb=1821a9c35f90be93d3d55b1e1db9114902eaf6f8;hpb=bd8f1aee7c78f7d2eb82bd4fec5e08c9e3d280ce diff --git a/src/SMESH_SWIG/ex19_sphereINcube.py b/src/SMESH_SWIG/ex19_sphereINcube.py index c303c7072..d59b048a5 100644 --- a/src/SMESH_SWIG/ex19_sphereINcube.py +++ b/src/SMESH_SWIG/ex19_sphereINcube.py @@ -1,5 +1,5 @@ # -*- coding: iso-8859-1 -*- -# Copyright (C) 2007-2013 CEA/DEN, EDF R&D, OPEN CASCADE +# Copyright (C) 2007-2014 CEA/DEN, EDF R&D, OPEN CASCADE # # Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN, # CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS @@ -7,7 +7,7 @@ # 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. +# version 2.1 of the License, or (at your option) any later version. # # This library is distributed in the hope that it will be useful, # but WITHOUT ANY WARRANTY; without even the implied warranty of @@ -23,10 +23,15 @@ # ================================== # -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 # ========= @@ -45,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)) @@ -69,74 +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"]) -sphere_decoupee = MakePartition([sphere_troue], [f1], [], [], ShapeType["SOLID"]) -sphere_decoupee = MakePartition([sphere_decoupee], [f2], [], [], ShapeType["SOLID"]) -sphere_decoupee = MakePartition([sphere_decoupee], [f3], [], [], ShapeType["SOLID"]) -sphere_decoupee = MakePartition([sphere_decoupee], [f4], [], [], 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 = geompy.MakeCompound([sphere_decoupee]) -sphere_partie = GetBlockNearPoint(sphere_decoupee, MakeVertex(-sphere_rayon, 0, 0)) -sphere_bloc = RemoveExtraEdges(sphere_partie) +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"]) -cube_decoupe = MakePartition([cube_trou], [f1], [], [], ShapeType["SOLID"]) -cube_decoupe = MakePartition([cube_decoupe], [f2], [], [], ShapeType["SOLID"]) -cube_decoupe = MakePartition([cube_decoupe], [f3], [], [], ShapeType["SOLID"]) -cube_decoupe = MakePartition([cube_decoupe], [f4], [], [], 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 = 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 # ================== @@ -144,24 +149,22 @@ 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 # ======= -smesh.SetCurrentStudy(salome.myStudy) - # Create a hexahedral mesh # ------------------------