src/Tools/blocFissure/gmu/construitFissureGenerale.py

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