src/Tools/blocFissure/gmu/creeZoneDefautMaillage.py

   1 # -*- coding: utf-8 -*-
   2
   3 import logging
   4 from .geomsmesh import geompy
   5 import math
   6 from .distance2 import distance2
   7 import traceback
   8 from .fissError import fissError
   9
  10 # -----------------------------------------------------------------------------
  11 # --- zone de defaut extraite du maillage
  12
  13 def creeZoneDefautMaillage(maillagesSains, shapeDefaut, tailleDefaut,
  14                            nomZones, coordsNoeudsFissure):
  15   """
  16   Identification de la zone à remailler, opérations sur le maillage
  17   de l'objet sain.
  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
  26   sont détectés
  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
  31   """
  32   logging.info("start")
  33
  34   maillageSain = maillagesSains[0]
  35   isHexa = maillagesSains[1]
  36   lists = maillageSain.CreateHoleSkin(tailleDefaut, shapeDefaut, nomZones, coordsNoeudsFissure)
  37
  38   logging.debug("lists=%s", lists)
  39
  40   trace = True
  41   origShapes = []
  42   verticesShapes = []
  43
  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
  46
  47   for aList in lists:
  48     aShape = aList[0]
  49     origShapes.append(aShape)
  50     logging.debug("  shapeId %s", aShape)
  51     vertices = []
  52     xyz0 = None
  53     for inode in range(1, len(aList)):
  54       xyz = maillageSain.GetNodeXYZ(aList[inode])
  55       if xyz0 is not None:
  56         cumul += distance2(xyz, xyz0)
  57         nb += 1
  58       xyz0 = xyz
  59       #logging.debug("    node %s %s", aList[inode], xyz)
  60       vertices.append(geompy.MakeVertex(xyz[0], xyz[1], xyz[2]))
  61       pass
  62     verticesShapes.append(vertices)
  63     pass
  64
  65   if (nb == 0) :
  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)
  71
  72   dmoyen = math.sqrt(cumul/nb) # ~ taille de l'arête moyenne du maillage global
  73   return origShapes, verticesShapes, dmoyen