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
22 from .geomsmesh import geompy
23 from .geomsmesh import geomPublish
24 from .geomsmesh import geomPublishInFather
26 from .extractionOrientee import extractionOrientee
27 from .getSubshapeIds import getSubshapeIds
29 # -----------------------------------------------------------------------------
31 # --- faces toriques et volumes du tore
33 def facesVolumesToriques(tore, plan, facesDefaut):
35 Extraction des deux faces et volumes du tore partitionné, qui suivent la génératrice elliptique.
36 @param tore : le tore partitionné et coupé.
37 @param plan : le plan de coupe
38 @return (facetore1,facetore2) les 2 faces selon la génératrice
42 centre = geompy.MakeVertexOnSurface(plan, 0.5, 0.5)
43 normal = geompy.GetNormal(plan, centre)
44 reference = geompy.MakeTranslationVector(centre, normal)
46 [facesInPlan, facesOutPlan, facesOnPlan] = extractionOrientee(plan, tore, reference, "FACE", 1.e-2, "faceTorePlan_")
47 [facesInSide, facesOutSide, facesOnSide] = extractionOrientee(facesDefaut, tore, reference, "FACE", 1.e-2, "faceTorePeau_")
48 facesIdInPlan = getSubshapeIds(tore, facesInPlan)
49 facesIdOutPlan = getSubshapeIds(tore, facesOutPlan)
50 facesIdOnSide = getSubshapeIds(tore, facesOnSide)
51 facesIdInSide = getSubshapeIds(tore, facesInSide)
52 facesIdOutSide = getSubshapeIds(tore, facesOutSide)
53 #facesIdInOutSide = facesIdInSide + facesIdOutSide
56 for i, faceId in enumerate(facesIdInPlan):
57 if faceId not in facesIdOnSide:
58 facetore1 = facesInPlan[i]
60 for i, faceId in enumerate(facesIdOutPlan):
61 if faceId not in facesIdOnSide:
62 facetore2 = facesOutPlan[i]
65 #[facetore1,facetore2] = geompy.GetShapesOnShape(pipe0, tore, geompy.ShapeType["FACE"], GEOM.ST_ON)
67 geomPublishInFather(initLog.debug, tore, facetore1, 'facetore1' )
68 geomPublishInFather(initLog.debug, tore, facetore2, 'facetore2' )
70 [volumeTore1, volumeTore2] = geompy.ExtractShapes(tore, geompy.ShapeType["SOLID"], True)
71 geomPublishInFather(initLog.debug, tore, volumeTore1, 'volumeTore1' )
72 geomPublishInFather(initLog.debug, tore, volumeTore2, 'volumeTore2' )
74 return facetore1, facetore2, volumeTore1, volumeTore2