1 # -*- coding: utf-8 -*-
2 # Copyright (C) 2014-2021 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 """Identification précise des edges et disques des faces de peau selon index extremité fissure"""
26 from .geomsmesh import geompy
27 from .geomsmesh import geomPublish
28 from .geomsmesh import geomPublishInFather
30 def identifieEdgesPeau_a(edgesFissExtPipe, facePeau, facesPeauSorted, edgesPeauFondIn, \
31 endsEdgeFond, facesPipePeau, edgeRadFacePipePeau, edgesListees):
32 """Identification précise des edges et disques des faces de peau selon index extremité fissure"""
35 for face in facesPeauSorted[:-1]: # la ou les faces débouchantes, pas la grande face de peau
36 logging.debug("examen face debouchante circulaire")
37 for i_aux,efep in enumerate(edgesFissExtPipe):
38 dist = geompy.MinDistance(face, efep)
39 logging.debug(" distance face circulaire edge %s", dist)
41 for j_aux, edpfi in enumerate(edgesPeauFondIn):
42 if geompy.MinDistance(face, edpfi) < 1e-3:
45 sharedVertices = geompy.GetSharedShapesMulti([face, edgesPeauFondIn[j_aux_0]], geompy.ShapeType["VERTEX"])
46 nameFace = "facePipePeau{}".format(i_aux)
47 nameVert = "endEdgeFond{}".format(i_aux)
48 nameEdge = "edgeRadFacePipePeau{}".format(i_aux)
49 facesPipePeau[i_aux] = face
50 endsEdgeFond[i_aux] = sharedVertices[0]
51 geomPublish(initLog.debug, face, nameFace)
52 geomPublish(initLog.debug, sharedVertices[0], nameVert)
53 edgesFace = geompy.ExtractShapes(face, geompy.ShapeType["EDGE"], True)
54 for edge in edgesFace:
55 if geompy.MinDistance(edge, sharedVertices[0]) < 1e-3:
56 edgeRadFacePipePeau[i_aux] = edge
57 geomPublish(initLog.debug, edge, nameEdge)
60 # --- edges circulaires de la face de peau et points de jonction de la face externe de fissure
61 logging.debug("facesPipePeau: %s", facesPipePeau)
62 edgesCircPeau = [None for _ in range(len(facesPipePeau))]
63 verticesCircPeau = [None for _ in range(len(facesPipePeau))]
64 for i_aux,fcirc in enumerate(facesPipePeau):
65 edges = geompy.GetSharedShapesMulti([facePeau, fcirc], geompy.ShapeType["EDGE"])
66 grpEdgesCirc = geompy.CreateGroup(facePeau, geompy.ShapeType["EDGE"])
67 geompy.UnionList(grpEdgesCirc, edges)
68 edgesCircPeau[i_aux] = grpEdgesCirc
69 name = "edgeCirc{}".format(i_aux)
70 geomPublishInFather(initLog.debug, facePeau, grpEdgesCirc, name)
71 edgesListees = edgesListees + edges
72 vertices = geompy.GetSharedShapesMulti([facePeau, fcirc], geompy.ShapeType["VERTEX"])
73 grpVertCircPeau = geompy.CreateGroup(facePeau, geompy.ShapeType["VERTEX"])
74 geompy.UnionList(grpVertCircPeau, vertices)
75 verticesCircPeau[i_aux] = grpVertCircPeau
76 name = "pointEdgeCirc{}".format(i_aux)
77 geomPublishInFather(initLog.debug, facePeau, grpVertCircPeau, name)
79 return edgesCircPeau, verticesCircPeau