1 # -*- coding: utf-8 -*-
2 # Copyright (C) 2014-2020 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
20 """construction de la fissure générale - maillage"""
27 from salome.smesh import smeshBuilder
31 from .geomsmesh import geompy
32 from .geomsmesh import smesh
34 from .putName import putName
35 from .enleveDefaut import enleveDefaut
36 from .shapeSurFissure import shapeSurFissure
37 from .regroupeSainEtDefaut import RegroupeSainEtDefaut
38 from .triedreBase import triedreBase
40 # -----------------------------------------------------------------------------
42 def construitFissureGenerale_c(maillageSain, meshBoiteDefaut, \
43 zoneDefaut, zoneDefaut_skin, zoneDefaut_internalFaces, zoneDefaut_internalEdges, \
45 maillageFissureParams, \
47 """construction de la fissure générale - maillage"""
49 logging.info("Usage du mailleur %s", mailleur)
51 nomRep = maillageFissureParams['nomRep']
52 nomFicSain = maillageFissureParams['nomFicSain']
53 nomFicFissure = maillageFissureParams['nomFicFissure']
55 O, _, _, _ = triedreBase()
57 # pour aider l'algo hexa-tetra à ne pas mettre de pyramides à l'exterieur des volumes repliés sur eux-mêmes
58 # on désigne les faces de peau en quadrangles par le groupe "skinFaces"
59 group_faceFissOutPipe = None
60 group_faceFissInPipe = None
61 groups = meshBoiteDefaut.GetGroups()
63 if grp.GetType() == SMESH.FACE:
64 if grp.GetName() == "fisOutPi":
65 group_faceFissOutPipe = grp
66 elif grp.GetName() == "fisInPi":
67 group_faceFissInPipe = grp
69 # le maillage NETGEN ne passe pas toujours ==> on force l'usage de MG_Tetra
71 logging.info("Maillage avec %s", mailleur)
72 if ( mailleur == "MeshGems"):
73 algo3d = meshBoiteDefaut.Tetrahedron(algo=smeshBuilder.MG_Tetra)
75 algo3d = meshBoiteDefaut.Tetrahedron(algo=smeshBuilder.NETGEN)
76 hypo3d = algo3d.MaxElementVolume(1000.0)
77 hypo3d.SetVerboseLevel( 0 )
78 hypo3d.SetStandardOutputLog( 0 )
79 hypo3d.SetRemoveLogOnSuccess( 1 )
80 putName(algo3d.GetSubMesh(), "boiteDefaut")
81 putName(algo3d, "algo3d_boiteDefaut")
82 putName(meshBoiteDefaut, "boiteDefaut")
84 is_done = meshBoiteDefaut.Compute()
85 text = "meshBoiteDefaut.Compute"
87 logging.info(text+" OK")
89 text = "Erreur au calcul du maillage.\n" + text
93 _ = meshBoiteDefaut.GetMesh().UnionListOfGroups( [ group_faceFissOutPipe, group_faceFissInPipe ], \
95 maillageSain = enleveDefaut(maillageSain, zoneDefaut, zoneDefaut_skin,
96 zoneDefaut_internalFaces, zoneDefaut_internalEdges)
97 putName(maillageSain, nomFicSain+"_coupe")
98 _, normfiss = shapeSurFissure(facesPortFissure)
99 maillageComplet = RegroupeSainEtDefaut(maillageSain, meshBoiteDefaut, \
100 None, None, 'COMPLET', normfiss)
102 logging.info("conversion quadratique")
103 maillageComplet.ConvertToQuadratic( 1 )
104 logging.info("groupes")
105 groups = maillageComplet.GetGroups()
106 grps = [ grp for grp in groups if grp.GetName() == 'FONDFISS']
107 _ = maillageComplet.GetMesh().CreateDimGroup( grps, SMESH.NODE, 'FONDFISS' )
109 logging.info("réorientation face de fissure FACE1")
110 grps = [ grp for grp in groups if grp.GetName() == 'FACE1']
111 _ = maillageComplet.Reorient2D( grps[0], normfiss, grps[0].GetID(1))
113 logging.info("réorientation face de fissure FACE2")
114 plansim = geompy.MakePlane(O, normfiss, 10000)
115 fissnorm = geompy.MakeMirrorByPlane(normfiss, plansim)
116 grps = [ grp for grp in groups if grp.GetName() == 'FACE2']
117 _ = maillageComplet.Reorient2D( grps[0], fissnorm, grps[0].GetID(1))
118 _ = maillageComplet.GetMesh().CreateDimGroup( grps, SMESH.NODE, 'FACE2' )
120 logging.info("export maillage fini")
121 fichierMaillageFissure = os.path.join (nomRep , '{}.med'.format(nomFicFissure))
122 maillageComplet.ExportMED(fichierMaillageFissure)
123 putName(maillageComplet, nomFicFissure)
124 logging.info("fichier maillage fissure %s", fichierMaillageFissure)
126 return maillageComplet