1 # -*- coding: utf-8 -*-
4 from .geomsmesh import geompy
6 from .distance2 import distance2
8 from .fissError import fissError
10 # -----------------------------------------------------------------------------
11 # --- zone de defaut extraite du maillage
13 def creeZoneDefautMaillage(maillagesSains, shapeDefaut, tailleDefaut,
14 nomZones, coordsNoeudsFissure):
16 Identification de la zone à remailler, opérations sur le maillage
18 La zone à remailler est définie à partir d'un objet géométrique
19 ou à partir d'un jeu de points et d'une distance d'influence.
20 @param maillagesSains : (le maillage de l'objet initial, booleen isHexa)
21 @param shapeDefaut : objet géométrique représentant la fissure
22 (selon les cas, un point central, ou une shape plus complexe,
23 dont on ne garde que les vertices)
24 @param tailleDefaut : distance d'influence définissant la zone à remailler:
25 tous les éléments du maillage initial qui pénètrent dans cette zone
27 @param nomZones : préfixe des noms de groupes créés dans le maillage initial. S'il y a un groupe de noeuds
28 @coordsNoeudsFissure : jeu de points donné par une liste (x1,y1,z1, x2,y2,z2, ...)
29 @return (origShapes, verticesShapes, dmoyen) liste id subShapes,
30 listes noeuds de bord, longueur arête moyenne bord
34 maillageSain = maillagesSains[0]
35 isHexa = maillagesSains[1]
36 lists = maillageSain.CreateHoleSkin(tailleDefaut, shapeDefaut, nomZones, coordsNoeudsFissure)
38 logging.debug("lists=%s", lists)
44 cumul = 0 # somme des distances carrées entre point ordonnés (taille des arêtes)
45 nb = 0 # nombre d'arêtes évaluées
49 origShapes.append(aShape)
50 logging.debug(" shapeId %s", aShape)
53 for inode in range(1, len(aList)):
54 xyz = maillageSain.GetNodeXYZ(aList[inode])
56 cumul += distance2(xyz, xyz0)
59 #logging.debug(" node %s %s", aList[inode], xyz)
60 vertices.append(geompy.MakeVertex(xyz[0], xyz[1], xyz[2]))
62 verticesShapes.append(vertices)
66 texte = "La zone à remailler n'est pas détectée correctement.<br>"
67 texte += "Cause possible :<ul>"
68 texte += "<li>La distance d'influence est trop petite. "
69 texte += "L'ordre de grandeur minimal correspond à la taille des mailles du maillage sain dans la zone à remailler.</li></ul>"
70 raise fissError(traceback.extract_stack(),texte)
72 dmoyen = math.sqrt(cumul/nb) # ~ taille de l'arête moyenne du maillage global
73 return origShapes, verticesShapes, dmoyen