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