X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FTools%2FblocFissure%2Fgmu%2FgetCentreFondFiss.py;h=709343082cce284eff42816b84464986f1e82aaa;hb=93f22a734d078eeadea90256715452905b3a69a6;hp=6c04d36275da6ca5fc7b6f9ded5f49cc69b07787;hpb=5d68554076bbca0e1e95fb0db215a6c2b84b6c54;p=modules%2Fsmesh.git diff --git a/src/Tools/blocFissure/gmu/getCentreFondFiss.py b/src/Tools/blocFissure/gmu/getCentreFondFiss.py index 6c04d3627..709343082 100644 --- a/src/Tools/blocFissure/gmu/getCentreFondFiss.py +++ b/src/Tools/blocFissure/gmu/getCentreFondFiss.py @@ -1,30 +1,54 @@ # -*- coding: utf-8 -*- +# Copyright (C) 2014-2021 EDF R&D +# +# This library is free software; you can redistribute it and/or +# modify it under the terms of the GNU Lesser General Public +# License as published by the Free Software Foundation; either +# version 2.1 of the License, or (at your option) any later version. +# +# This library is distributed in the hope that it will be useful, +# but WITHOUT ANY WARRANTY; without even the implied warranty of +# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +# Lesser General Public License for more details. +# +# You should have received a copy of the GNU Lesser General Public +# License along with this library; if not, write to the Free Software +# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA +# +# See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com +# +"""Identification du centre de fond de fissure,""" import logging -from geomsmesh import geompy import bisect +from . import initLog + +from .geomsmesh import geompy +from .geomsmesh import geomPublish +from .geomsmesh import geomPublishInFather + publie = False def getCentreFondFiss(shapesFissure): - """ - identification du centre de fond de fissure, - transformation fond de fissure en edge unique (seulement pour la procédure insereFissureGenerale). + """Identification du centre de fond de fissure, + + transformation fond de fissure en edge unique (seulement pour la procédure construitFissureGenerale). On distingue le cas d'utilisation de la procédure insereFissureLongue par le nombre d'éléments de shapesFissure. """ global publie logging.debug("start") - + fondFiss = shapesFissure[4] # groupe d'edges de fond de fissure - if len(shapesFissure) == 6: # procédure insereFissureGenerale, et edge fond de fissure fournie explicitement + if len(shapesFissure) == 6: # procédure construitFissureGenerale, et edge fond de fissure fournie explicitement edgeFondExt = shapesFissure[5] else: edgeFondExt = None - + if len(shapesFissure) > 6: # procédure insereFissureLongue (fissure plane, plusieurs edges sur le fond de fissure) centreFondFiss = shapesFissure[1] tgtCentre = None - else: # procédure insereFissureGenerale, détermination edge unique et milieu de l'edge + else: # procédure construitFissureGenerale, détermination edge unique et milieu de l'edge if geompy.NumberOfEdges(fondFiss) > 1: if geompy.NbShapes(fondFiss, geompy.ShapeType["WIRE"]) > 0: # wire aWire = fondFiss @@ -35,8 +59,8 @@ def getCentreFondFiss(shapesFissure): edgeFondExt = fondFiss aWire = geompy.MakeWire([fondFiss], 1e-07) if not publie: - geompy.addToStudy(aWire, "wireFondFissExt") - + geomPublish(initLog.debug, aWire, "wireFondFissExt") + lgWire = geompy.BasicProperties(aWire)[0] edges = geompy.ExtractShapes(aWire, geompy.ShapeType["EDGE"], True) lgEdges = [geompy.BasicProperties(ed)[0] for ed in edges] @@ -50,7 +74,7 @@ def getCentreFondFiss(shapesFissure): logging.debug("lgsumEdges %s", lgSumEd) logging.debug("id edge: %s, lgOnEdge: %s, lgEdge: %s",iedr, lgOnEdge, lgEdges[iedr]) if iedr > 0: # il y a une edge avant celle du milieu - if geompy.MinDistance(edges[iedr-1], geompy.MakeVertexOnCurve(edges[iedr], 0.0 )) < 1.e-3: # edge orientée croissante + if geompy.MinDistance(edges[iedr-1], geompy.MakeVertexOnCurve(edges[iedr], 0.0 )) < 1.e-3: # edge orientée croissante centreFondFiss = geompy.MakeVertexOnCurve(edges[iedr], lgOnEdge/lgEdges[iedr]) else: centreFondFiss = geompy.MakeVertexOnCurve(edges[iedr], 1.0 - lgOnEdge/lgEdges[iedr]) @@ -61,9 +85,9 @@ def getCentreFondFiss(shapesFissure): centreFondFiss = geompy.MakeVertexOnCurve(edges[iedr], 1.0 - lgOnEdge/lgEdges[iedr]) else: # on ne sait pas comment est orientée l'edge unique, mais ça n'a pas d'importance centreFondFiss = geompy.MakeVertexOnCurve(edges[iedr], lgOnEdge/lgEdges[iedr]) - geompy.addToStudyInFather(aWire, centreFondFiss, "centreFondFiss") + geomPublishInFather(initLog.debug,aWire, centreFondFiss, "centreFondFiss") tgtCentre = geompy.MakeTangentOnCurve(edges[iedr], lgOnEdge/ lgEdges[iedr]) - + if edgeFondExt is None: # fond de fissure non fourni explicitement sous forme d'edge try: edgeFondExt = geompy.MakeEdgeWire(aWire, 0.0005, 1e-07) @@ -71,7 +95,8 @@ def getCentreFondFiss(shapesFissure): logging.debug("erreur MakeEdgeWire sur fond de fissure, on fait sans") edgeFondExt = None if not publie and edgeFondExt is not None: - geompy.addToStudy(edgeFondExt, "edgeFondExt") - + geomPublish(initLog.debug, edgeFondExt, "edgeFondExt") + publie = True + return edgeFondExt, centreFondFiss, tgtCentre