src/Tools/blocFissure/gmu/construitEdgesRadialesDebouchantes.py

   1 # -*- coding: utf-8 -*-
   2
   3 import logging
   4
   5 from geomsmesh import geompy
   6 import GEOM
   7 from sortEdges import sortEdges
   8
   9 def construitEdgesRadialesDebouchantes(idisklim, idiskout, gptsdisks, raydisks,
  10                                        facesPipePeau, edgeRadFacePipePeau, nbsegCercle):
  11   """
  12   construction des listes d'edges radiales sur chaque extrémité débouchante
  13   """
  14   logging.info('start')
  15
  16   # --- listes de nappes radiales en filling à chaque extrémité débouchante
  17
  18   facesDebouchantes = [False, False]
  19   idFacesDebouchantes = [-1, -1] # contiendra les indices des faces disque débouchantes (facesPipePeau)
  20   listNappes =[]
  21   for i, idisk in enumerate(idisklim):
  22     numout = idiskout[i]
  23     logging.debug("extremité %s, indices disques interne %s, externe %s",i, idisk, numout)
  24     nappes = []
  25     if  (idisk != 0) and (idisk != len(gptsdisks)-1): # si extrémité débouchante
  26       for k in range(nbsegCercle):
  27         if i == 0:
  28           iddeb = max(0, numout)
  29           idfin = max(iddeb+3,idisk+1) # il faut 3 rayons pour faire un filling qui suive le fond de fissure
  30           #logging.debug("extremité %s, indices retenus interne %s, externe %s",i, idfin, iddeb)
  31           comp = geompy.MakeCompound(raydisks[k][iddeb:idfin])
  32           name='compoundRay%d'%k
  33           geompy.addToStudy(comp, name)
  34         else:
  35           idfin = min(len(gptsdisks), numout+1)
  36           iddeb = min(idfin-3, idisk) # il faut 3 rayons pour faire un filling qui suive le fond de fissure
  37           #logging.debug("extremité %s, indices retenus interne %s, externe %s",i, idfin, iddeb)
  38           comp = geompy.MakeCompound(raydisks[k][iddeb:idfin])
  39           name='compoundRay%d'%k
  40           geompy.addToStudy(comp, name)
  41         nappe = geompy.MakeFilling(comp, 2, 5, 0.0001, 0.0001, 0, GEOM.FOM_Default)
  42         nappes.append(nappe)
  43         name='nappe%d'%k
  44         geompy.addToStudy(nappe, name)
  45         facesDebouchantes[i] = True
  46     listNappes.append(nappes)
  47
  48   # --- mise en correspondance avec les indices des faces disque débouchantes (facesPipePeau)
  49   for i, nappes in enumerate(listNappes):
  50     if facesDebouchantes[i]:
  51       for k, face in enumerate(facesPipePeau):
  52         #logging.debug('i, k, face, nappes[0] %s %s %s %s', i, k, face, nappes[0])
  53         #geompy.addToStudy(nappes[0], 'lanappe')
  54         #geompy.addToStudy(face, 'laface')
  55         edge = geompy.MakeSection(face, nappes[0])
  56         if geompy.NbShapes(edge, geompy.ShapeType["EDGE"]) > 0:
  57           idFacesDebouchantes[i] = k
  58           break
  59   logging.debug("idFacesDebouchantes: %s", idFacesDebouchantes)
  60
  61   # --- construction des listes d'edges radiales sur chaque extrémité débouchante
  62   listEdges = []
  63   for i, nappes in enumerate(listNappes):
  64     ifd = idFacesDebouchantes[i] # indice de face débouchante (facesPipePeau)
  65     if ifd < 0:
  66       listEdges.append([])
  67     else:
  68       face = facesPipePeau[ifd]
  69       edges = [edgeRadFacePipePeau[ifd]]
  70       for k, nappe in enumerate(nappes):
  71         if k > 0:
  72           obj = geompy.MakeSection(face, nappes[k]) # normalement une edge, parfois un compound d'edges dont un tout petit
  73           edge = obj
  74           vs = geompy.ExtractShapes(obj, geompy.ShapeType["VERTEX"], False)
  75           if len(vs) > 2:
  76             eds = geompy.ExtractShapes(obj, geompy.ShapeType["EDGE"], False)
  77             [edsorted, minl,maxl] = sortEdges(eds)
  78             edge = edsorted[-1]
  79           else:
  80             maxl = geompy.BasicProperties(edge)[0]
  81           if maxl < 0.01: # problème MakeSection
  82             logging.debug("problème MakeSection recherche edge radiale %s, longueur trop faible: %s, utilisation partition", k, maxl)
  83             partNappeFace = geompy.MakePartition([face, nappes[k]], [] , [], [], geompy.ShapeType["FACE"], 0, [], 0)
  84             edps= geompy.ExtractShapes(partNappeFace, geompy.ShapeType["EDGE"], False)
  85             ednouv = []
  86             for ii, ed in enumerate(edps):
  87               vxs = geompy.ExtractShapes(ed, geompy.ShapeType["VERTEX"], False)
  88               distx = [geompy.MinDistance(vx, face) for vx in vxs]
  89               distx += [geompy.MinDistance(vx, nappes[k]) for vx in vxs]
  90               dmax = max(distx)
  91               logging.debug("  dmax %s",dmax)
  92               if dmax < 0.01:
  93                 ednouv.append(ed)
  94             logging.debug("  edges issues de la partition: %s", ednouv)
  95             for ii, ed in enumerate(ednouv):
  96               geompy.addToStudy(ed, "ednouv%d"%ii)
  97             [edsorted, minl,maxl] = sortEdges(ednouv)
  98             logging.debug("  longueur edge trouvée: %s", maxl)
  99             edge = edsorted[-1]
 100           edges.append(edge)
 101           name = 'edgeEndPipe%d'%k
 102           geompy.addToStudy(edge, name)
 103       listEdges.append(edges)
 104
 105   return (listEdges, idFacesDebouchantes)