src/Tools/blocFissure/gmu/peauInterne.py

   1 # -*- coding: utf-8 -*-
   2
   3 import logging
   4 from geomsmesh import smesh
   5 import SMESH
   6 import traceback
   7 from fissError import fissError
   8
   9 from listOfExtraFunctions import lookForCorner
  10 from fusionMaillageAttributionDefaut import fusionMaillageDefaut
  11
  12 # -----------------------------------------------------------------------------
  13 # --- peau interne du defaut dans le maillage sain
  14
  15 def peauInterne(fichierMaillage, shapeDefaut, nomZones):
  16   """
  17   Retrouve les groupes de défaut dans le maillage sain modifié par CreateHoleSkin (CreeZoneDefautMaillage)
  18   On récupère le volume et la peau de la zone de défaut, les éventuelles faces et arêtes internes de cette zone.
  19   """
  20   logging.info("start")
  21   ([maillageSain], status) = smesh.CreateMeshesFromMED(fichierMaillage)
  22
  23   groups = maillageSain.GetGroups()
  24   zoneDefaut = None
  25   for grp in groups:
  26     logging.debug("groupe %s",grp.GetName())
  27     if grp.GetName() == nomZones + "_vol":
  28       zoneDefaut = grp
  29       break
  30   zoneDefaut_skin = None
  31   for grp in groups:
  32     if grp.GetName() == nomZones + "_skin":
  33       zoneDefaut_skin = grp
  34       break
  35   zoneDefaut_internalFaces = None
  36   for grp in groups:
  37     if grp.GetName() == nomZones + "_internalFaces":
  38       zoneDefaut_internalFaces = grp
  39       break
  40   zoneDefaut_internalEdges = None
  41   for grp in groups:
  42     if grp.GetName() == nomZones + "_internalEdges":
  43       zoneDefaut_internalEdges = grp
  44       break
  45
  46   # --- Le groupe ZoneDefaut ne doit contenir que des Hexaèdres"
  47
  48   info=maillageSain.GetMeshInfo(zoneDefaut)
  49   keys = info.keys(); keys.sort()
  50   nbelem=0
  51   nbhexa=0
  52   for i in keys:
  53     #print "  %s  :  %d" % ( i, info[i] )
  54     nbelem+=info[i]
  55     if "Entity_Hexa" in str(i):
  56       nbhexa+=info[i]
  57   if (nbelem == 0) or (nbhexa < nbelem) :
  58     texte = "La zone a remailler est incorrecte.<br>"
  59     texte += "Causes possibles :<ul>"
  60     texte += "<li>Les mailles à enlever dans le maillage sain n'ont pas été détectées.</li>"
  61     texte += "<li>Certaines faces du maillage sain sont à l'envers : les normales aux faces en paroi de volume doivent être sortantes.</li>"
  62     texte += "<li>Il n'y a pas que des Hexaèdres réglés linéaires dans la zone à remailler (notamment mailles quadratiques, tetraèdres non traités)</li></ul>"
  63     raise fissError(traceback.extract_stack(),texte)
  64
  65   nbAdded, maillageSain, DefautBoundary = maillageSain.MakeBoundaryElements( SMESH.BND_2DFROM3D, 'DefBound', '', 0, [ zoneDefaut ])
  66   internal = maillageSain.GetMesh().CutListOfGroups( [ DefautBoundary ], [ zoneDefaut_skin ], 'internal' )
  67   internalBoundary = smesh.CopyMesh( internal, 'internalBoundary', 0, 0)
  68
  69   maillageDefautCible = smesh.CopyMesh(zoneDefaut_skin, 'maillageCible', 0, 0)
  70   listOfCorner = lookForCorner(maillageDefautCible)
  71   logging.debug("listOfCorner = %s", listOfCorner)
  72   if len(listOfCorner) > 0:
  73       logging.info("présence de coins à la surface externe de la zone à reconstruire")
  74       zoneDefaut_skin, internalBoundary = fusionMaillageDefaut(maillageSain, maillageDefautCible, internalBoundary, zoneDefaut_skin, shapeDefaut, listOfCorner)
  75
  76   return maillageSain, internalBoundary, zoneDefaut, zoneDefaut_skin, zoneDefaut_internalFaces, zoneDefaut_internalEdges
  77