src/Tools/blocFissure/gmu/casStandard.py

   1 # -*- coding: utf-8 -*-
   2 # Copyright (C) 2014-2022  EDF R&D
   3 #
   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.
   8 #
   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.
  13 #
  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
  17 #
  18 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  19 #
  20 """Cas standard"""
  21
  22 import os
  23 import logging
  24
  25 from .geomsmesh import geompy, smesh
  26 from .geomsmesh import geomPublish
  27 from .geomsmesh import geomPublishInFather
  28
  29 from . import initLog
  30
  31 from .fissureGenerique import fissureGenerique
  32
  33 from .initEtude import initEtude
  34 from .triedreBase import triedreBase
  35 from .genereMeshCalculZoneDefaut import genereMeshCalculZoneDefaut
  36 from .creeZoneDefautDansObjetSain import creeZoneDefautDansObjetSain
  37 from .construitFissureGenerale import construitFissureGenerale
  38 from .putName import putName
  39
  40 O, OX, OY, OZ = triedreBase()
  41
  42 class casStandard(fissureGenerique):
  43   """problème de fissure standard, défini par :
  44
  45   - un maillage sain (hexaèdres),
  46   - une face géométrique de fissure, qui doit légèrement dépasser hors du volume maillé
  47   - les noms des groupes d'arêtes ou leurs numéros d'arêtes (edges au sens de GEOM) correspondant au fond de fissure
  48   - les paramètres de maillage de la fissure
  49   """
  50   referencesMaillageFissure = None
  51
  52   # ---------------------------------------------------------------------------
  53   def __init__ (self, dicoParams, references = None, numeroCas = None):
  54     initEtude()
  55     self.references = references
  56     self.dicoParams = dicoParams
  57     self.numeroCas = numeroCas
  58
  59     if 'nomProbleme' in self.dicoParams:
  60       self.nomProbleme = self.dicoParams['nomProbleme']
  61
  62     if 'nomCas' in self.dicoParams:
  63       self.nomCas = self.dicoParams['nomCas']
  64     elif 'nomres' in self.dicoParams:
  65       self.nomCas = os.path.splitext(os.path.split(self.dicoParams['nomres'])[1])[0]
  66     elif ( self.numeroCas is not None ):
  67       self.nomCas = self.nomProbleme +"_%d"%(self.numeroCas)
  68     else:
  69       self.nomCas = 'casStandard'
  70
  71     if 'reptrav' in self.dicoParams:
  72       self.reptrav = self.dicoParams['reptrav']
  73     else:
  74       self.reptrav = os.curdir
  75
  76     if 'lenSegPipe' in self.dicoParams:
  77       self.lenSegPipe = self.dicoParams['lenSegPipe']
  78     else:
  79       self.lenSegPipe = self.dicoParams['rayonPipe']
  80
  81     if 'step' in self.dicoParams:
  82       step = self.dicoParams['step']
  83     else:
  84       step = -1 # exécuter toutes les étapes
  85
  86     if 'aretesVives' not in self.dicoParams:
  87       self.dicoParams['aretesVives'] = 0
  88
  89     # valeur par défaut : exécution immédiate, sinon execution différée dans le cas d'une liste de problèmes
  90     if ( self.numeroCas is None ):
  91       self.executeProbleme(step)
  92
  93   # ---------------------------------------------------------------------------
  94   def genereMaillageSain(self, geometriesSaines, meshParams):
  95     logging.info("genereMaillageSain %s", self.nomCas)
  96
  97     ([objetSain], _) = smesh.CreateMeshesFromMED(self.dicoParams['maillageSain'])
  98     putName(objetSain.GetMesh(), objetSain.GetName(), i_pref=self.numeroCas)
  99
 100     return [objetSain, True] # True : maillage hexa
 101
 102   # ---------------------------------------------------------------------------
 103   def setParamShapeFissure(self):
 104     """
 105     paramètres de la fissure pour méthode insereFissureGenerale
 106     lgInfluence : distance autour de la shape de fissure a remailler (A ajuster selon le maillage)
 107     rayonPipe   : le rayon du pile maillé en hexa autour du fond de fissure
 108     convexe     : True : la face est convexe (vue de l'exterieur) sert si on ne donne pas de point interne
 109     pointIn_x   : optionnel : coordonnée x d'un point dans le solide sain (pour orienter la face - idem avec y,z)
 110     """
 111     logging.info("setParamShapeFissure %s", self.nomCas)
 112     if 'pointInterieur' in self.dicoParams:
 113       self.shapeFissureParams = dict(lgInfluence = self.dicoParams['lgInfluence'],
 114                                      rayonPipe   = self.dicoParams['rayonPipe'],
 115                                      lenSegPipe  = self.lenSegPipe,
 116                                      pointIn_x   = self.dicoParams['pointInterieur'][0],
 117                                      pointIn_y   = self.dicoParams['pointInterieur'][1],
 118                                      pointIn_z   = self.dicoParams['pointInterieur'][2])
 119     else:
 120       self.shapeFissureParams = dict(lgInfluence = self.dicoParams['lgInfluence'],
 121                                      rayonPipe   = self.dicoParams['rayonPipe'],
 122                                      lenSegPipe  = self.lenSegPipe)
 123
 124   # ---------------------------------------------------------------------------
 125   def genereShapeFissure( self, geometriesSaines, geomParams, shapeFissureParams, \
 126                                 mailleur="MeshGems"):
 127
 128     lgInfluence = shapeFissureParams['lgInfluence']
 129
 130 #   Contrôle de 'brepFaceFissure' pour les anciennes versions
 131     if ( 'brepFaceFissure' in self.dicoParams ):
 132       self.dicoParams['CAOFaceFissure'] = self.dicoParams['brepFaceFissure']
 133     cao_file = self.dicoParams['CAOFaceFissure']
 134     suffix = os.path.basename(cao_file).split(".")[-1]
 135     if ( suffix.upper() == "BREP" ):
 136       shellFiss = geompy.ImportBREP(cao_file)
 137     elif ( suffix.upper() == "XAO" ):
 138       (_, shellFiss, _, l_groups, _) = geompy.ImportXAO(cao_file)
 139     fondFiss = geompy.CreateGroup(shellFiss, geompy.ShapeType["EDGE"])
 140 #   Contrôle de 'edgeFissIds' pour les anciennes versions
 141     if ( 'edgeFissIds' in self.dicoParams ):
 142       self.dicoParams['edgeFiss'] = self.dicoParams['edgeFissIds']
 143     if isinstance(self.dicoParams['edgeFiss'][0],int):
 144       geompy.UnionIDs(fondFiss, self.dicoParams['edgeFiss'] )
 145     else:
 146       l_groups = geompy.GetGroups(shellFiss)
 147       l_aux = list()
 148       for group in l_groups:
 149         if ( group.GetName() in self.dicoParams['edgeFiss'] ):
 150           l_aux.append(group)
 151       geompy.UnionList(fondFiss, l_aux )
 152     geomPublish(initLog.debug, shellFiss, 'shellFiss' )
 153     geomPublishInFather(initLog.debug, shellFiss, fondFiss, 'fondFiss' )
 154
 155
 156     coordsNoeudsFissure = genereMeshCalculZoneDefaut(shellFiss, self.dicoParams['meshBrep'][0] ,self.dicoParams['meshBrep'][1], \
 157                                                      mailleur, self.numeroCas)
 158
 159     centre = None
 160     return [shellFiss, centre, lgInfluence, coordsNoeudsFissure, fondFiss]
 161
 162   # ---------------------------------------------------------------------------
 163   def setParamMaillageFissure(self):
 164     self.maillageFissureParams = dict(nomRep           = self.reptrav,
 165                                       nomFicSain       = self.nomCas,
 166                                       nomFicFissure    = self.nomProbleme + "_fissure",
 167                                       nbsegRad         = self.dicoParams['nbSegRad'],
 168                                       nbsegCercle      = self.dicoParams['nbSegCercle'],
 169                                       areteFaceFissure = self.dicoParams['areteFaceFissure'],
 170                                       aretesVives      = self.dicoParams['aretesVives'])
 171
 172   # ---------------------------------------------------------------------------
 173   def genereZoneDefaut(self, geometriesSaines, maillagesSains, shapesFissure, shapeFissureParams, maillageFissureParams):
 174     elementsDefaut = creeZoneDefautDansObjetSain(geometriesSaines, maillagesSains, shapesFissure, shapeFissureParams, maillageFissureParams, \
 175                           self.numeroCas)
 176     return elementsDefaut
 177
 178   # ---------------------------------------------------------------------------
 179   def genereMaillageFissure(self, geometriesSaines, maillagesSains, \
 180                                   shapesFissure, shapeFissureParams, \
 181                                   maillageFissureParams, elementsDefaut, step, \
 182                                   mailleur="MeshGems"):
 183     maillageFissure = construitFissureGenerale(shapesFissure, shapeFissureParams, \
 184                                                maillageFissureParams, elementsDefaut, \
 185                                                mailleur, self.numeroCas)
 186     return maillageFissure
 187
 188   # ---------------------------------------------------------------------------
 189   def setReferencesMaillageFissure(self):
 190     if self.references is not None:
 191       self.referencesMaillageFissure = self.references
 192     else:
 193       self.referencesMaillageFissure = dict( \
 194                                             Entity_Quad_Quadrangle = 0, \
 195                                             Entity_Quad_Hexa = 0, \
 196                                             Entity_Node = 0, \
 197                                             Entity_Quad_Edge = 0, \
 198                                             Entity_Quad_Triangle = 0, \
 199                                             Entity_Quad_Tetra = 0, \
 200                                             Entity_Quad_Pyramid = 0, \
 201                                             Entity_Quad_Penta = 0 \
 202                                            )