X-Git-Url: http://git.salome-platform.org/gitweb/?p=modules%2Fsmesh.git;a=blobdiff_plain;f=src%2FTools%2FblocFissure%2Fgmu%2FpartitionneFissureParPipe.py;h=c3279b0a39643992f256361682e0f2639b66d858;hp=7d6289797d00843eac503d3ba6d55c3c77646098;hb=93f22a734d078eeadea90256715452905b3a69a6;hpb=5482b99d07dd144fd5be299e722f39a81de3b5be diff --git a/src/Tools/blocFissure/gmu/partitionneFissureParPipe.py b/src/Tools/blocFissure/gmu/partitionneFissureParPipe.py index 7d6289797..c3279b0a3 100644 --- a/src/Tools/blocFissure/gmu/partitionneFissureParPipe.py +++ b/src/Tools/blocFissure/gmu/partitionneFissureParPipe.py @@ -1,15 +1,39 @@ # -*- coding: utf-8 -*- +# Copyright (C) 2014-2021 EDF R&D +# +# 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. +# +# 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 +# +"""Pipe de fond de fissure, prolongé, partition face fissure par pipe""" -import math import logging -from geomsmesh import geompy -from geomsmesh import geomPublish -from geomsmesh import geomPublishInFather -import initLog -from findWireEndVertices import findWireEndVertices -from prolongeWire import prolongeWire +import math +import traceback + +from .geomsmesh import geompy +from .geomsmesh import geomPublish +from .geomsmesh import geomPublishInFather -def partitionneFissureParPipe(shapesFissure, elementsDefaut, rayonPipe): +from . import initLog +from .findWireEndVertices import findWireEndVertices +from .prolongeWire import prolongeWire +from .fissError import fissError + +def partitionneFissureParPipe(shapesFissure, elementsDefaut, rayonPipe, \ + nro_cas=None): """ pipe de fond de fissure, prolongé, partition face fissure par pipe identification des edges communes pipe et face fissure @@ -24,41 +48,46 @@ def partitionneFissureParPipe(shapesFissure, elementsDefaut, rayonPipe): plan = geompy.MakePlane(centreFondFiss, tgtCentre, 10000) shapeDefaut = geompy.MakePartition([shapeDefaut], [plan], [], [], geompy.ShapeType["FACE"], 0, [], 0) #fondFissCoupe = geompy.GetInPlaceByHistory(shapeDefaut, fondFiss) #= inutile - geomPublish(initLog.debug, shapeDefaut, 'shapeDefaut_coupe') - #geomPublishInFather(initLog.debug,shapeDefaut, fondFissCoupe, 'fondFiss_coupe') - + geomPublish(initLog.debug, shapeDefaut, 'shapeDefaut_coupe', nro_cas) + #geomPublishInFather(initLog.debug,shapeDefaut, fondFissCoupe, 'fondFiss_coupe', nro_cas) + extrem, norms = findWireEndVertices(fondFiss, True) logging.debug("extrem: %s, norm: %s",extrem, norms) cercle = geompy.MakeCircle(extrem[0], norms[0], rayonPipe) cercle = geompy.MakeRotation(cercle, norms[0], math.pi/3.0 ) # éviter d'avoir l'arête de couture du pipe presque confondue avec la face fissure - geomPublish(initLog.debug, cercle, 'cercle') + geomPublish(initLog.debug, cercle, 'cercle', nro_cas) fondFissProlonge = prolongeWire(fondFiss, extrem, norms, 2*rayonPipe) - pipeFiss = geompy.MakePipe(cercle, fondFissProlonge) - geomPublish(initLog.debug, pipeFiss, 'pipeFiss') + try: + pipeFiss = geompy.MakePipe(cercle, fondFissProlonge) + except: + texte = "génération du pipe le long de la ligne de fond de fissure prolongée impossible. " + texte += "Cause possible : la ligne s'autointersecte lorsqu'on la prolonge." + raise fissError(traceback.extract_stack(),texte) + geomPublish(initLog.debug, pipeFiss, 'pipeFiss', nro_cas) partFissPipe = geompy.MakePartition([shapeDefaut, pipeFiss], [], [], [], geompy.ShapeType["FACE"], 0, [], 1) - geomPublish(initLog.debug, partFissPipe, 'partFissPipe') + geomPublish(initLog.debug, partFissPipe, 'partFissPipe', nro_cas) fissPipe = geompy.GetInPlaceByHistory(partFissPipe, shapeDefaut) - geomPublish(initLog.debug, fissPipe, 'fissPipe') + geomPublish(initLog.debug, fissPipe, 'fissPipe', nro_cas) partPipe = geompy.GetInPlaceByHistory(partFissPipe, pipeFiss) - geomPublish(initLog.debug, partPipe, 'partPipe') - + geomPublish(initLog.debug, partPipe, 'partPipe', nro_cas) + edgesPipeFiss = geompy.GetSharedShapesMulti([fissPipe, partPipe], geompy.ShapeType["EDGE"]) - for i, edge in enumerate(edgesPipeFiss): - name = "edgePipe%d"%i - geomPublishInFather(initLog.debug,fissPipe, edge, name) + for i_aux, edge in enumerate(edgesPipeFiss): + name = "edgePipe{}".format(i_aux) + geomPublishInFather(initLog.debug,fissPipe, edge, name, nro_cas) try: wirePipeFiss = geompy.MakeWire(edgesPipeFiss) except: wirePipeFiss = geompy.MakeCompound(edgesPipeFiss) logging.debug("wirePipeFiss construit sous forme de compound") - geomPublish(initLog.debug, wirePipeFiss, "wirePipeFiss") - + geomPublish(initLog.always, wirePipeFiss, "wirePipeFiss", nro_cas) + wireFondFiss = geompy.GetInPlace(partFissPipe,fondFiss) edgesFondFiss = geompy.GetSharedShapesMulti([fissPipe, wireFondFiss], geompy.ShapeType["EDGE"]) - for i, edge in enumerate(edgesFondFiss): - name = "edgeFondFiss%d"%i + for i_aux, edge in enumerate(edgesFondFiss): + name = "edgeFondFiss{}".format(i_aux) geomPublishInFather(initLog.debug,fissPipe, edge, name) wireFondFiss = geompy.MakeWire(edgesFondFiss) - geomPublish(initLog.debug, wireFondFiss,"wireFondFiss") + geomPublish(initLog.always, wireFondFiss, "wireFondFiss", nro_cas) - return (fissPipe, edgesPipeFiss, edgesFondFiss, wirePipeFiss, wireFondFiss) \ No newline at end of file + return (fissPipe, edgesPipeFiss, edgesFondFiss, wirePipeFiss, wireFondFiss)