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
24 from .distance2 import distance2
26 from .fissError import fissError
28 # -----------------------------------------------------------------------------
29 # --- zone de defaut extraite du maillage
31 def creeZoneDefautMaillage(maillagesSains, shapeDefaut, tailleDefaut,
32 nomZones, coordsNoeudsFissure):
34 Identification de la zone à remailler, opérations sur le maillage
36 La zone à remailler est définie à partir d'un objet géométrique
37 ou à partir d'un jeu de points et d'une distance d'influence.
38 @param maillagesSains : (le maillage de l'objet initial, booleen isHexa)
39 @param shapeDefaut : objet géométrique représentant la fissure
40 (selon les cas, un point central, ou une shape plus complexe,
41 dont on ne garde que les vertices)
42 @param tailleDefaut : distance d'influence définissant la zone à remailler:
43 tous les éléments du maillage initial qui pénètrent dans cette zone
45 @param nomZones : préfixe des noms de groupes créés dans le maillage initial. S'il y a un groupe de noeuds
46 @coordsNoeudsFissure : jeu de points donné par une liste (x1,y1,z1, x2,y2,z2, ...)
47 @return (origShapes, verticesShapes, dmoyen) liste id subShapes,
48 listes noeuds de bord, longueur arête moyenne bord
52 maillageSain = maillagesSains[0]
53 isHexa = maillagesSains[1]
54 lists = maillageSain.CreateHoleSkin(tailleDefaut, shapeDefaut, nomZones, coordsNoeudsFissure)
56 logging.debug("lists={}".format(lists))
60 verticesShapes = list()
62 cumul = 0 # somme des distances carrées entre point ordonnés (taille des arêtes)
63 nb = 0 # nombre d'arêtes évaluées
67 origShapes.append(aShape)
68 logging.debug(" shapeId {}".format(aShape))
71 for inode in range(1, len(aList)):
72 xyz = maillageSain.GetNodeXYZ(aList[inode])
74 cumul += distance2(xyz, xyz0)
77 #logging.debug(" node %s %s", aList[inode], xyz)
78 vertices.append(geompy.MakeVertex(xyz[0], xyz[1], xyz[2]))
80 verticesShapes.append(vertices)
84 texte = "La zone à remailler n'est pas détectée correctement.<br>"
85 texte += "Cause possible :<ul>"
86 texte += "<li>La distance d'influence est trop petite. "
87 texte += "L'ordre de grandeur minimal correspond à la taille des mailles du maillage sain dans la zone à remailler.</li></ul>"
88 raise fissError(traceback.extract_stack(),texte)
90 dmoyen = math.sqrt(cumul/nb) # ~ taille de l'arête moyenne du maillage global
91 return origShapes, verticesShapes, dmoyen