src/Tools/blocFissure/gmu/construitFissureGenerale.py

   1 # -*- coding: utf-8 -*-
   2 # Copyright (C) 2014-2020  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 """procédure complète de construction d'une fissure générale"""
  21
  22 import logging
  23 import salome
  24 from .geomsmesh import geompy
  25 from .geomsmesh import geomPublishInFather
  26 from . import initLog
  27 import GEOM
  28 from .geomsmesh import smesh
  29 from salome.smesh import smeshBuilder
  30 import SMESH
  31
  32 from .orderEdgesFromWire import orderEdgesFromWire
  33 from .restreintFaceFissure import restreintFaceFissure
  34 from .partitionneFissureParPipe import partitionneFissureParPipe
  35 from .construitPartitionsPeauFissure import construitPartitionsPeauFissure
  36 from .compoundFromList import compoundFromList
  37 from .identifieFacesEdgesFissureExterne import identifieFacesEdgesFissureExterne
  38 from .calculePointsAxiauxPipe import calculePointsAxiauxPipe
  39 from .elimineExtremitesPipe import elimineExtremitesPipe
  40 from .construitEdgesRadialesDebouchantes import construitEdgesRadialesDebouchantes
  41 from .creePointsPipePeau import creePointsPipePeau
  42 from .ajustePointsEdgePipeFissure import ajustePointsEdgePipeFissure
  43 from .construitMaillagePipe import construitMaillagePipe
  44 from .mailleAretesEtJonction import mailleAretesEtJonction
  45 from .mailleFacesFissure import mailleFacesFissure
  46 from .mailleFacesPeau import mailleFacesPeau
  47
  48 from .construitFissureGenerale_a import construitFissureGenerale_a
  49 from .construitFissureGenerale_b import construitFissureGenerale_b
  50 from .construitFissureGenerale_c import construitFissureGenerale_c
  51
  52 # -----------------------------------------------------------------------------
  53
  54 def construitFissureGenerale(shapesFissure, shapeFissureParams, \
  55                              maillageFissureParams, elementsDefaut, \
  56                              step=-1, mailleur="MeshGems"):
  57   """procédure complète fissure générale"""
  58   logging.info('start')
  59   logging.info("Usage du mailleur %s", mailleur)
  60
  61   shapeDefaut       = shapesFissure[0] # faces de fissure, débordant
  62
  63   rayonPipe = shapeFissureParams['rayonPipe']
  64   if 'lenSegPipe' in shapeFissureParams:
  65     lenSegPipe = shapeFissureParams['lenSegPipe']
  66   else:
  67     lenSegPipe = rayonPipe
  68
  69   nbsegRad          = maillageFissureParams['nbsegRad']      # nombre de couches selon un rayon du pipe
  70   nbsegCercle       = maillageFissureParams['nbsegCercle']   # nombre de secteur dans un cercle du pipe
  71   areteFaceFissure  = maillageFissureParams['areteFaceFissure']
  72   lgAretesVives     = 0
  73   if 'aretesVives' in maillageFissureParams:
  74     lgAretesVives   = maillageFissureParams['aretesVives']
  75
  76 # Point interne
  77   pointInterne = construitFissureGenerale_a (geompy, shapeFissureParams)
  78
  79   # fillings des faces en peau
  80   facesDefaut              = elementsDefaut[0]
  81   #centresDefaut            = elementsDefaut[1]
  82   #normalsDefaut            = elementsDefaut[2]
  83   #extrusionsDefaut         = elementsDefaut[3]
  84   dmoyen                   = elementsDefaut[4]
  85   bordsPartages            = elementsDefaut[5]
  86   #fillconts                = elementsDefaut[6]
  87   #idFilToCont              = elementsDefaut[7]
  88   maillageSain             = elementsDefaut[8]
  89   internalBoundary         = elementsDefaut[9]
  90   zoneDefaut               = elementsDefaut[10]
  91   zoneDefaut_skin          = elementsDefaut[11]
  92   zoneDefaut_internalFaces = elementsDefaut[12]
  93   zoneDefaut_internalEdges = elementsDefaut[13]
  94   #edgeFondExt              = elementsDefaut[14]
  95   centreFondFiss           = elementsDefaut[15]
  96   #tgtCentre                = elementsDefaut[16]
  97   if lgAretesVives == 0:
  98     lgAretesVives = dmoyen
  99
 100   # --- restriction de la face de fissure au domaine solide :
 101   #     partition face fissure étendue par fillings, on garde la face interne
 102
 103   facesPortFissure = restreintFaceFissure(shapeDefaut, facesDefaut, pointInterne)
 104
 105   # --- pipe de fond de fissure, prolongé, partition face fissure par pipe
 106   #     identification des edges communes pipe et face fissure
 107
 108   (fissPipe, edgesPipeFiss, edgesFondFiss, wirePipeFiss, wireFondFiss) = \
 109                           partitionneFissureParPipe(shapesFissure, elementsDefaut, rayonPipe)
 110   edgesFondFiss, edgesIdByOrientation = orderEdgesFromWire(wireFondFiss)
 111
 112   for i,edge in enumerate(edgesFondFiss):
 113     geomPublishInFather(initLog.debug, wireFondFiss, edge, "edgeFondFiss%d"%i)
 114
 115   # --- peau et face de fissure
 116   #
 117   # --- partition peau défaut - face de fissure prolongée - wire de fond de fissure prolongée
 118   #     il peut y avoir plusieurs faces externes, dont certaines sont découpées par la fissure
 119   #     liste de faces externes : facesDefaut
 120   #     liste de partitions face externe - fissure : partitionPeauFissFond (None quand pas d'intersection)
 121
 122   partitionsPeauFissFond = construitPartitionsPeauFissure(facesDefaut, fissPipe)
 123
 124   # --- arêtes vives détectées (dans quadranglesToShapeNoCorner
 125   #                             et quadranglesToShapeWithCorner)
 126
 127   aretesVivesC = compoundFromList(bordsPartages, "areteVive")
 128   aretesVivesCoupees = list()  # ensemble des arêtes vives identifiées sur les faces de peau dans l'itération sur partitionsPeauFissFond
 129
 130   # --- inventaire des faces de peau coupées par la fissure
 131
 132   ptEdgeFond, fsFissuExt, edFisExtPe, edFisExtPi, \
 133     facesPeaux, edCircPeau, ptCircPeau, gpedgeBord, gpedgeVifs, edFissPeau, ptFisExtPi, \
 134     edgeRadFacePipePeau, facesPipePeau = \
 135     construitFissureGenerale_b( partitionsPeauFissFond, \
 136                                 edgesPipeFiss, edgesFondFiss, wireFondFiss, aretesVivesC, \
 137                                 facesDefaut, centreFondFiss, rayonPipe, aretesVivesCoupees )
 138
 139   # --- identification des faces et edges de fissure externe pour maillage
 140
 141   (faceFissureExterne, edgesPipeFissureExterneC, wirePipeFissureExterne, edgesPeauFissureExterneC) = \
 142       identifieFacesEdgesFissureExterne(fsFissuExt, edFisExtPe, edFisExtPi, edgesPipeFiss)
 143
 144   # --- preparation maillage du pipe :
 145   #     - détections des points a respecter : jonction des edges/faces constituant la face de fissure externe au pipe
 146   #     - points sur les edges de fond de fissure et edges pipe/face fissure,
 147   #     - vecteurs tangents au fond de fissure (normal au disque maillé)
 148
 149   (centres, gptsdisks, raydisks) = calculePointsAxiauxPipe(edgesFondFiss, edgesIdByOrientation, facesDefaut, \
 150                                                            centreFondFiss, wireFondFiss, wirePipeFiss, \
 151                                                            lenSegPipe, rayonPipe, nbsegCercle, nbsegRad)
 152
 153   # --- recherche des points en trop (externes au volume à remailler)
 154   #     - on associe chaque extrémité du pipe à une face filling
 155   #     - on part des disques aux extrémités du pipe
 156   #     - pour chaque disque, on prend les vertices de géométrie,
 157   #       on marque leur position relative à la face.
 158   #     - on s'arrete quand tous les noeuds sont dedans
 159
 160   (idFillingFromBout, idisklim, idiskout) = elimineExtremitesPipe(ptEdgeFond, facesDefaut, centres, gptsdisks, nbsegCercle)
 161
 162   # --- construction des listes d'edges radiales sur chaque extrémité débouchante
 163
 164   (listEdges, idFacesDebouchantes) = construitEdgesRadialesDebouchantes(idisklim, idiskout, gptsdisks, raydisks, \
 165                                                                         facesPipePeau, edgeRadFacePipePeau, nbsegCercle)
 166
 167   # --- création des points du maillage du pipe sur la face de peau
 168
 169   (gptsdisks, idisklim) = creePointsPipePeau(listEdges, idFacesDebouchantes, idFillingFromBout, \
 170                                              ptEdgeFond, ptFisExtPi, edCircPeau, gptsdisks, idisklim, nbsegRad)
 171
 172   # --- ajustement precis des points sur edgesPipeFissureExterneC
 173
 174   gptsdisks = ajustePointsEdgePipeFissure(edgesPipeFissureExterneC, wirePipeFissureExterne, gptsdisks, idisklim)
 175
 176    # --- maillage effectif du pipe
 177
 178   (meshPipe, meshPipeGroups, edgesCircPipeGroup) = construitMaillagePipe(gptsdisks, idisklim, nbsegCercle, nbsegRad)
 179
 180   # --- edges de bord, faces défaut à respecter
 181
 182   (internalBoundary, bordsLibres, grpAretesVives) = mailleAretesEtJonction(internalBoundary, aretesVivesCoupees, lgAretesVives)
 183
 184   # --- maillage faces de fissure
 185
 186   (meshFaceFiss, _, grpEdgesPeauFissureExterne, _) = \
 187       mailleFacesFissure(faceFissureExterne, \
 188                          edgesPipeFissureExterneC, edgesPeauFissureExterneC, \
 189                          meshPipeGroups, areteFaceFissure, rayonPipe, nbsegRad, \
 190                          mailleur)
 191
 192   # --- maillage faces de peau
 193
 194   meshesFacesPeau = mailleFacesPeau(partitionsPeauFissFond, idFillingFromBout, facesDefaut, \
 195                                     facesPeaux, edCircPeau, ptCircPeau, gpedgeBord, gpedgeVifs, edFissPeau, \
 196                                     bordsLibres, grpEdgesPeauFissureExterne, grpAretesVives, \
 197                                     edgesCircPipeGroup, dmoyen, rayonPipe, nbsegRad, \
 198                                     mailleur)
 199
 200   # --- regroupement des maillages du défaut
 201
 202   listMeshes = [internalBoundary.GetMesh(), meshPipe.GetMesh(), meshFaceFiss.GetMesh()]
 203   for mfpeau in meshesFacesPeau:
 204     listMeshes.append(mfpeau.GetMesh())
 205
 206   meshBoiteDefaut = smesh.Concatenate(listMeshes, 1, 1, 1e-05,False)
 207
 208 # Maillage complet
 209   maillageComplet = construitFissureGenerale_c( maillageSain, meshBoiteDefaut, \
 210                                                 zoneDefaut, zoneDefaut_skin, zoneDefaut_internalFaces, zoneDefaut_internalEdges, \
 211                                                 facesPortFissure, \
 212                                                 maillageFissureParams, \
 213                                                 mailleur )
 214
 215   if salome.sg.hasDesktop():
 216     salome.sg.updateObjBrowser()
 217
 218   logging.info("maillage fissure fini")
 219
 220   return maillageComplet