src/Tools/blocFissure/gmu/partitionneFissureParPipe.py

   1 # -*- coding: utf-8 -*-
   2 # Copyright (C) 2014-2019  CEA/DEN, EDF R&D
   3 #
   4 # This library is free software; you can redistribute it and/or
   5 # modify it under the terms of the GNU Lesser General Public
   6 # License as published by the Free Software Foundation; either
   7 # version 2.1 of the License, or (at your option) any later version.
   8 #
   9 # This library is distributed in the hope that it will be useful,
  10 # but WITHOUT ANY WARRANTY; without even the implied warranty of
  11 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  12 # Lesser General Public License for more details.
  13 #
  14 # You should have received a copy of the GNU Lesser General Public
  15 # License along with this library; if not, write to the Free Software
  16 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  17 #
  18 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  19 #
  20
  21 import math
  22 import logging
  23 from .geomsmesh import geompy
  24 from .geomsmesh import geomPublish
  25 from .geomsmesh import geomPublishInFather
  26 from . import initLog
  27 from .findWireEndVertices import findWireEndVertices
  28 from .prolongeWire import prolongeWire
  29 import traceback
  30 from .fissError import fissError
  31
  32 def partitionneFissureParPipe(shapesFissure, elementsDefaut, rayonPipe):
  33   """
  34   pipe de fond de fissure, prolongé, partition face fissure par pipe
  35   identification des edges communes pipe et face fissure
  36   """
  37   logging.info('start')
  38   shapeDefaut       = shapesFissure[0] # faces de fissure, débordant
  39   fondFiss          = shapesFissure[4] # groupe d'edges de fond de fissure
  40   centreFondFiss    = elementsDefaut[15]
  41   tgtCentre         = elementsDefaut[16]
  42
  43   if geompy.NumberOfFaces(shapeDefaut) == 1:
  44     plan = geompy.MakePlane(centreFondFiss, tgtCentre, 10000)
  45     shapeDefaut = geompy.MakePartition([shapeDefaut], [plan], [], [], geompy.ShapeType["FACE"], 0, [], 0)
  46     #fondFissCoupe = geompy.GetInPlaceByHistory(shapeDefaut, fondFiss) #= inutile
  47     geomPublish(initLog.debug, shapeDefaut, 'shapeDefaut_coupe')
  48     #geomPublishInFather(initLog.debug,shapeDefaut, fondFissCoupe, 'fondFiss_coupe')
  49
  50   extrem, norms = findWireEndVertices(fondFiss, True)
  51   logging.debug("extrem: %s, norm: %s",extrem, norms)
  52   cercle = geompy.MakeCircle(extrem[0], norms[0], rayonPipe)
  53   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
  54   geomPublish(initLog.debug, cercle, 'cercle')
  55   fondFissProlonge = prolongeWire(fondFiss, extrem, norms, 2*rayonPipe)
  56   try:
  57     pipeFiss = geompy.MakePipe(cercle, fondFissProlonge)
  58   except:
  59     texte = "génération du pipe le long de la ligne de fond de fissure prolongée impossible. "
  60     texte += "Cause possible : la ligne s'autointersecte lorsqu'on la prolonge."
  61     raise fissError(traceback.extract_stack(),texte)
  62   geomPublish(initLog.debug, pipeFiss, 'pipeFiss')
  63   partFissPipe = geompy.MakePartition([shapeDefaut, pipeFiss], [], [], [], geompy.ShapeType["FACE"], 0, [], 1)
  64   geomPublish(initLog.debug, partFissPipe, 'partFissPipe')
  65   fissPipe = geompy.GetInPlaceByHistory(partFissPipe, shapeDefaut)
  66   geomPublish(initLog.debug, fissPipe, 'fissPipe')
  67   partPipe = geompy.GetInPlaceByHistory(partFissPipe, pipeFiss)
  68   geomPublish(initLog.debug, partPipe, 'partPipe')
  69
  70   edgesPipeFiss = geompy.GetSharedShapesMulti([fissPipe, partPipe], geompy.ShapeType["EDGE"])
  71   for i, edge in enumerate(edgesPipeFiss):
  72     name = "edgePipe%d"%i
  73     geomPublishInFather(initLog.debug,fissPipe, edge, name)
  74   try:
  75     wirePipeFiss = geompy.MakeWire(edgesPipeFiss)
  76   except:
  77     wirePipeFiss = geompy.MakeCompound(edgesPipeFiss)
  78     logging.debug("wirePipeFiss construit sous forme de compound")
  79   geomPublish(initLog.debug, wirePipeFiss, "wirePipeFiss")
  80
  81   wireFondFiss = geompy.GetInPlace(partFissPipe,fondFiss)
  82   edgesFondFiss = geompy.GetSharedShapesMulti([fissPipe, wireFondFiss], geompy.ShapeType["EDGE"])
  83   for i, edge in enumerate(edgesFondFiss):
  84     name = "edgeFondFiss%d"%i
  85     geomPublishInFather(initLog.debug,fissPipe, edge, name)
  86   wireFondFiss = geompy.MakeWire(edgesFondFiss)
  87   geomPublish(initLog.debug, wireFondFiss,"wireFondFiss")
  88
  89   return (fissPipe, edgesPipeFiss, edgesFondFiss, wirePipeFiss, wireFondFiss)