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, \
172 return elementsDefaut
174 # ---------------------------------------------------------------------------
175 def genereMaillageFissure(self, geometriesSaines, maillagesSains, \
176 shapesFissure, shapeFissureParams, \
177 maillageFissureParams, elementsDefaut, step, \
178 mailleur="MeshGems"):
179 maillageFissure = construitFissureGenerale(shapesFissure, shapeFissureParams, \
180 maillageFissureParams, elementsDefaut, \
181 step, mailleur, self.numeroCas)
182 return maillageFissure
184 # ---------------------------------------------------------------------------
185 def setReferencesMaillageFissure(self):
186 if self.references is not None:
187 self.referencesMaillageFissure = self.references
189 self.referencesMaillageFissure = dict( \
190 Entity_Quad_Quadrangle = 0, \
191 Entity_Quad_Hexa = 0, \
193 Entity_Quad_Edge = 0, \
194 Entity_Quad_Triangle = 0, \
195 Entity_Quad_Tetra = 0, \
196 Entity_Quad_Pyramid = 0, \
197 Entity_Quad_Penta = 0 \