src/Tools/blocFissure/gmu/partitionneFissureParPipe.py

   1 # -*- coding: utf-8 -*-
   2
   3 import math
   4 import logging
   5 from .geomsmesh import geompy
   6 from .geomsmesh import geomPublish
   7 from .geomsmesh import geomPublishInFather
   8 from . import initLog
   9 from .findWireEndVertices import findWireEndVertices
  10 from .prolongeWire import prolongeWire
  11 import traceback
  12 from .fissError import fissError
  13
  14 def partitionneFissureParPipe(shapesFissure, elementsDefaut, rayonPipe):
  15   """
  16   pipe de fond de fissure, prolongé, partition face fissure par pipe
  17   identification des edges communes pipe et face fissure
  18   """
  19   logging.info('start')
  20   shapeDefaut       = shapesFissure[0] # faces de fissure, débordant
  21   fondFiss          = shapesFissure[4] # groupe d'edges de fond de fissure
  22   centreFondFiss    = elementsDefaut[15]
  23   tgtCentre         = elementsDefaut[16]
  24
  25   if geompy.NumberOfFaces(shapeDefaut) == 1:
  26     plan = geompy.MakePlane(centreFondFiss, tgtCentre, 10000)
  27     shapeDefaut = geompy.MakePartition([shapeDefaut], [plan], [], [], geompy.ShapeType["FACE"], 0, [], 0)
  28     #fondFissCoupe = geompy.GetInPlaceByHistory(shapeDefaut, fondFiss) #= inutile
  29     geomPublish(initLog.debug, shapeDefaut, 'shapeDefaut_coupe')
  30     #geomPublishInFather(initLog.debug,shapeDefaut, fondFissCoupe, 'fondFiss_coupe')
  31
  32   extrem, norms = findWireEndVertices(fondFiss, True)
  33   logging.debug("extrem: %s, norm: %s",extrem, norms)
  34   cercle = geompy.MakeCircle(extrem[0], norms[0], rayonPipe)
  35   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
  36   geomPublish(initLog.debug, cercle, 'cercle')
  37   fondFissProlonge = prolongeWire(fondFiss, extrem, norms, 2*rayonPipe)
  38   try:
  39     pipeFiss = geompy.MakePipe(cercle, fondFissProlonge)
  40   except:
  41     texte = "génération du pipe le long de la ligne de fond de fissure prolongée impossible. "
  42     texte += "Cause possible : la ligne s'autointersecte lorsqu'on la prolonge."
  43     raise fissError(traceback.extract_stack(),texte)
  44   geomPublish(initLog.debug, pipeFiss, 'pipeFiss')
  45   partFissPipe = geompy.MakePartition([shapeDefaut, pipeFiss], [], [], [], geompy.ShapeType["FACE"], 0, [], 1)
  46   geomPublish(initLog.debug, partFissPipe, 'partFissPipe')
  47   fissPipe = geompy.GetInPlaceByHistory(partFissPipe, shapeDefaut)
  48   geomPublish(initLog.debug, fissPipe, 'fissPipe')
  49   partPipe = geompy.GetInPlaceByHistory(partFissPipe, pipeFiss)
  50   geomPublish(initLog.debug, partPipe, 'partPipe')
  51
  52   edgesPipeFiss = geompy.GetSharedShapesMulti([fissPipe, partPipe], geompy.ShapeType["EDGE"])
  53   for i, edge in enumerate(edgesPipeFiss):
  54     name = "edgePipe%d"%i
  55     geomPublishInFather(initLog.debug,fissPipe, edge, name)
  56   try:
  57     wirePipeFiss = geompy.MakeWire(edgesPipeFiss)
  58   except:
  59     wirePipeFiss = geompy.MakeCompound(edgesPipeFiss)
  60     logging.debug("wirePipeFiss construit sous forme de compound")
  61   geomPublish(initLog.debug, wirePipeFiss, "wirePipeFiss")
  62
  63   wireFondFiss = geompy.GetInPlace(partFissPipe,fondFiss)
  64   edgesFondFiss = geompy.GetSharedShapesMulti([fissPipe, wireFondFiss], geompy.ShapeType["EDGE"])
  65   for i, edge in enumerate(edgesFondFiss):
  66     name = "edgeFondFiss%d"%i
  67     geomPublishInFather(initLog.debug,fissPipe, edge, name)
  68   wireFondFiss = geompy.MakeWire(edgesFondFiss)
  69   geomPublish(initLog.debug, wireFondFiss,"wireFondFiss")
  70
  71   return (fissPipe, edgesPipeFiss, edgesFondFiss, wirePipeFiss, wireFondFiss)