1 # -*- coding: utf-8 -*-
2 # Copyright (C) 2014-2020 EDF R&D
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.
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.
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
18 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
20 """Construit les arêtes débouchantes"""
24 from .geomsmesh import geompy
25 from .geomsmesh import geomPublish
26 #from .geomsmesh import geomPublishInFather
29 from .sortEdges import sortEdges
31 def construitEdgesRadialesDebouchantes(idisklim, idiskout, gptsdisks, raydisks,
32 facesPipePeau, edgeRadFacePipePeau, nbsegCercle):
33 """construction des listes d'edges radiales sur chaque extrémité débouchante"""
36 # --- listes de nappes radiales en filling à chaque extrémité débouchante
38 facesDebouchantes = [False, False]
39 idFacesDebouchantes = [-1, -1] # contiendra les indices des faces disque débouchantes (facesPipePeau)
41 for i, idisk in enumerate(idisklim):
43 logging.debug("extremité %s, indices disques interne %s, externe %s",i, idisk, numout)
45 if (idisk != 0) and (idisk != len(gptsdisks)-1): # si extrémité débouchante
46 for k in range(nbsegCercle):
48 iddeb = max(0, numout)
49 idfin = max(iddeb+3,idisk+1) # il faut 3 rayons pour faire un filling qui suive le fond de fissure
50 #logging.debug("extremité %s, indices retenus interne %s, externe %s",i, idfin, iddeb)
51 comp = geompy.MakeCompound(raydisks[k][iddeb:idfin])
52 name='compoundRay%d'%k
53 geomPublish(initLog.debug, comp, name)
55 idfin = min(len(gptsdisks), numout+1)
56 iddeb = min(idfin-3, idisk) # il faut 3 rayons pour faire un filling qui suive le fond de fissure
57 #logging.debug("extremité %s, indices retenus interne %s, externe %s",i, idfin, iddeb)
58 comp = geompy.MakeCompound(raydisks[k][iddeb:idfin])
59 name='compoundRay%d'%k
60 geomPublish(initLog.debug, comp, name)
61 nappe = geompy.MakeFilling(comp, 2, 5, 0.0001, 0.0001, 0, GEOM.FOM_Default)
64 geomPublish(initLog.debug, nappe, name)
65 facesDebouchantes[i] = True
66 listNappes.append(nappes)
68 # --- mise en correspondance avec les indices des faces disque débouchantes (facesPipePeau)
69 for i, nappes in enumerate(listNappes):
70 if facesDebouchantes[i]:
71 for k, face in enumerate(facesPipePeau):
72 #logging.debug('i, k, face, nappes[0] %s %s %s %s', i, k, face, nappes[0])
73 #geomPublish(initLog.debug, nappes[0], 'lanappe')
74 #geomPublish(initLog.debug, face, 'laface')
75 edge = geompy.MakeSection(face, nappes[0])
76 if geompy.NbShapes(edge, geompy.ShapeType["EDGE"]) > 0:
77 idFacesDebouchantes[i] = k
79 logging.debug("idFacesDebouchantes: %s", idFacesDebouchantes)
81 # --- construction des listes d'edges radiales sur chaque extrémité débouchante
83 for i, nappes in enumerate(listNappes):
84 ifd = idFacesDebouchantes[i] # indice de face débouchante (facesPipePeau)
88 face = facesPipePeau[ifd]
89 edges = [edgeRadFacePipePeau[ifd]]
90 for k, nappe in enumerate(nappes):
92 obj = geompy.MakeSection(face, nappes[k]) # normalement une edge, parfois un compound d'edges dont un tout petit
94 vs = geompy.ExtractShapes(obj, geompy.ShapeType["VERTEX"], False)
96 eds = geompy.ExtractShapes(obj, geompy.ShapeType["EDGE"], False)
97 [edsorted, _, maxl] = sortEdges(eds)
100 maxl = geompy.BasicProperties(edge)[0]
101 if maxl < 0.01: # problème MakeSection
102 logging.info("problème MakeSection recherche edge radiale %s, longueur trop faible: %s, utilisation partition", k, maxl)
103 partNappeFace = geompy.MakePartition([face, nappes[k]], [] , [], [], geompy.ShapeType["FACE"], 0, [], 0)
104 edps= geompy.ExtractShapes(partNappeFace, geompy.ShapeType["EDGE"], False)
106 for ii, ed in enumerate(edps):
108 vxs = geompy.ExtractShapes(ed, geompy.ShapeType["VERTEX"], False)
109 distx = [geompy.MinDistance(vx, face) for vx in vxs]
110 distx += [geompy.MinDistance(vx, nappes[k]) for vx in vxs]
112 lgedge = geompy.BasicProperties(ed)[0]
113 logging.debug(" dmax %s, longueur edge %s",dmax, lgedge)
114 if dmax < 0.01 and lgedge > 0.01:
116 if (len(ednouv) > 0):
117 logging.debug(" edges issues de la partition: %s", ednouv)
118 for ii, ed in enumerate(ednouv):
119 geomPublish(initLog.debug, ed, "ednouv%d"%ii)
120 [edsorted, _, maxl] = sortEdges(ednouv)
121 logging.debug(" longueur edge trouvée: %s", maxl)
124 logging.info("problème partition recherche edge radiale %s", k)
125 vxs = geompy.ExtractShapes(partNappeFace, geompy.ShapeType["VERTEX"], False)
127 for ii,vx in enumerate(vxs):
128 distx = geompy.MinDistance(vx, face)
129 distx += geompy.MinDistance(vx, nappes[k])
130 logging.debug("vertex distance: %s", distx)
133 logging.debug("nombre vertex candidats %s", len(vxnouv))
135 eds = [geompy.MakeEdge(vxnouv[j],vxnouv[(j+1)%len(vxnouv)]) for j in range(len(vxnouv))]
136 [edsorted2, _, maxl] = sortEdges(eds)
138 logging.debug("lg edge: %s", maxl)
140 logging.debug("problème recherche edge radiale %s non résolu", k)
142 name = 'edgeEndPipe%d'%k
143 geomPublish(initLog.debug, edge, name)
144 listEdges.append(edges)
146 return (listEdges, idFacesDebouchantes)