X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;ds=sidebyside;f=src%2FTools%2FblocFissure%2Fgmu%2FcreeZoneDefautMaillage.py;h=fcef6543283fe59b0aa3238cfb87d39b5daf44fb;hb=60bb6273b99dc88e3d6f360d7ae5792a77c68d06;hp=fa6edc8498e37f297c21a2065aa09f813078cdf7;hpb=d5029840731bccaa1718e65f0abf3b19198c7293;p=modules%2Fsmesh.git
diff --git a/src/Tools/blocFissure/gmu/creeZoneDefautMaillage.py b/src/Tools/blocFissure/gmu/creeZoneDefautMaillage.py
index fa6edc849..fcef65432 100644
--- a/src/Tools/blocFissure/gmu/creeZoneDefautMaillage.py
+++ b/src/Tools/blocFissure/gmu/creeZoneDefautMaillage.py
@@ -1,73 +1,96 @@
# -*- 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
+#
+"""Zone de defaut extraite du maillage"""
import logging
-from .geomsmesh import geompy
import math
-from .distance2 import distance2
import traceback
+
+from .geomsmesh import geompy
+from .distance2 import distance2
from .fissError import fissError
# -----------------------------------------------------------------------------
-# --- zone de defaut extraite du maillage
def creeZoneDefautMaillage(maillagesSains, shapeDefaut, tailleDefaut,
nomZones, coordsNoeudsFissure):
- """
- Identification de la zone à remailler, opérations sur le maillage
- de l'objet sain.
+ """Identification de la zone à remailler, opérations sur le maillage de l'objet sain.
+
La zone à remailler est définie à partir d'un objet géométrique
ou à partir d'un jeu de points et d'une distance d'influence.
@param maillagesSains : (le maillage de l'objet initial, booleen isHexa)
@param shapeDefaut : objet géométrique représentant la fissure
- (selon les cas, un point central, ou une shape plus complexe,
- dont on ne garde que les vertices)
+ (selon les cas, un point central, ou une shape plus complexe, dont on ne garde que les vertices)
@param tailleDefaut : distance d'influence définissant la zone à remailler:
tous les éléments du maillage initial qui pénètrent dans cette zone
sont détectés
@param nomZones : préfixe des noms de groupes créés dans le maillage initial. S'il y a un groupe de noeuds
@coordsNoeudsFissure : jeu de points donné par une liste (x1,y1,z1, x2,y2,z2, ...)
- @return (origShapes, verticesShapes, dmoyen) liste id subShapes,
- listes noeuds de bord, longueur arête moyenne bord
+ @return (origShapes, verticesShapes, dmoyen) liste id subShapes, listes noeuds de bord, longueur arête moyenne bord
"""
logging.info("start")
-
+ texte = "distance d'influence (tailleDefaut) = {}".format(tailleDefaut)
+ logging.debug(texte)
+ #print (" shapeDefaut : {}".format(shapeDefaut))
+ #print (" nomZones : {}".format(nomZones))
+
maillageSain = maillagesSains[0]
- isHexa = maillagesSains[1]
+ #isHexa = maillagesSains[1]
lists = maillageSain.CreateHoleSkin(tailleDefaut, shapeDefaut, nomZones, coordsNoeudsFissure)
- logging.debug("lists=%s", lists)
+ #print("lists = {}".format(lists))
trace = True
- origShapes = []
- verticesShapes = []
+ origShapes = list()
+ verticesShapes = list()
- cumul = 0 # somme des distances carrées entre point ordonnés (taille des arêtes)
- nb = 0 # nombre d'arêtes évaluées
+ cumul = 0. # somme des distances carrées entre points ordonnés (taille des arêtes)
+ nbar = 0 # nombre d'arêtes évaluées
for aList in lists:
aShape = aList[0]
origShapes.append(aShape)
- logging.debug(" shapeId %s", aShape)
- vertices = []
- xyz0 = None
- for inode in range(1, len(aList)):
- xyz = maillageSain.GetNodeXYZ(aList[inode])
- if xyz0 is not None:
- cumul += distance2(xyz, xyz0)
- nb += 1
+ vertices = list()
+ xyz0 = maillageSain.GetNodeXYZ(aList[1])
+ ###vertices.append(geompy.MakeVertex(xyz0[0], xyz0[1], xyz0[2])) pourquoi pas le premier ?
+ #print (" node 0 {} : {}".format(aList[1],xyz0))
+ #print (xyz0)
+ for node in aList[2:]:
+ xyz = maillageSain.GetNodeXYZ(node)
+ #print (" node {} : {}".format(node,xyz))
+ cumul += distance2(xyz, xyz0)
xyz0 = xyz
- #logging.debug(" node %s %s", aList[inode], xyz)
+ #logging.debug(" node %s %s", node, xyz)
vertices.append(geompy.MakeVertex(xyz[0], xyz[1], xyz[2]))
- pass
+ nbar += len(aList) - 2
verticesShapes.append(vertices)
- pass
+ #print ("nbar = {}".format(nbar))
+ #print ("cumul = {}".format(cumul))
- if (nb == 0) :
+ if (nbar == 0):
texte = "La zone à remailler n'est pas détectée correctement.
"
texte += "Cause possible :