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