# -*- coding: utf-8 -*-
import logging
-from geomsmesh import smesh
+from .geomsmesh import smesh
import SMESH
+import traceback
+from .fissError import fissError
-from listOfExtraFunctions import lookForCorner
-from fusionMaillageAttributionDefaut import fusionMaillageDefaut
+from .listOfExtraFunctions import lookForCorner
+from .fusionMaillageAttributionDefaut import fusionMaillageDefaut
# -----------------------------------------------------------------------------
# --- peau interne du defaut dans le maillage sain
groups = maillageSain.GetGroups()
zoneDefaut = None
for grp in groups:
- #print " ",grp.GetName()
+ logging.debug("groupe %s",grp.GetName())
if grp.GetName() == nomZones + "_vol":
zoneDefaut = grp
break
if grp.GetName() == nomZones + "_internalEdges":
zoneDefaut_internalEdges = grp
break
+
+ # --- Le groupe ZoneDefaut ne doit contenir que des Hexaèdres"
+
+ info=maillageSain.GetMeshInfo(zoneDefaut)
+ keys = list(info.keys()); keys.sort()
+ nbelem=0
+ nbhexa=0
+ for i in keys:
+ #print " %s : %d" % ( i, info[i] )
+ nbelem+=info[i]
+ if "Entity_Hexa" in str(i):
+ nbhexa+=info[i]
+ if (nbelem == 0) or (nbhexa < nbelem) :
+ texte = "La zone a remailler est incorrecte.<br>"
+ texte += "Causes possibles :<ul>"
+ texte += "<li>Les mailles à enlever dans le maillage sain n'ont pas été détectées.</li>"
+ texte += "<li>Certaines faces du maillage sain sont à l'envers : les normales aux faces en paroi de volume doivent être sortantes.</li>"
+ 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>"
+ raise fissError(traceback.extract_stack(),texte)
nbAdded, maillageSain, DefautBoundary = maillageSain.MakeBoundaryElements( SMESH.BND_2DFROM3D, 'DefBound', '', 0, [ zoneDefaut ])
internal = maillageSain.GetMesh().CutListOfGroups( [ DefautBoundary ], [ zoneDefaut_skin ], 'internal' )
maillageDefautCible = smesh.CopyMesh(zoneDefaut_skin, 'maillageCible', 0, 0)
listOfCorner = lookForCorner(maillageDefautCible)
- print "listOfCorner = ", listOfCorner
+ logging.debug("listOfCorner = %s", listOfCorner)
if len(listOfCorner) > 0:
- print " /!\ SUITE DU SCRIPT EN CONSTRUCTION /!\\"
+ logging.info("présence de coins à la surface externe de la zone à reconstruire")
zoneDefaut_skin, internalBoundary = fusionMaillageDefaut(maillageSain, maillageDefautCible, internalBoundary, zoneDefaut_skin, shapeDefaut, listOfCorner)
return maillageSain, internalBoundary, zoneDefaut, zoneDefaut_skin, zoneDefaut_internalFaces, zoneDefaut_internalEdges