src/Tools/blocFissure/gmu/partitionneFissureParPipe.py

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