src/Tools/blocFissure/CasTests/fissureCoude_5.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 """Fissure dans un coude"""
  21
  22 import os
  23
  24 from blocFissure.gmu.fissureCoude  import fissureCoude
  25
  26 class fissureCoude_5(fissureCoude):
  27   """problème de fissure du Coude : - adaptation maillage"""
  28
  29   nomProbleme = "fissureCoude_5"
  30   geomParams = dict()
  31   meshParams = dict()
  32   shapeFissureParams = dict()
  33   maillageFissureParams = dict()
  34   referencesMaillageFissure = dict()
  35
  36 # ---------------------------------------------------------------------------
  37   def setParamGeometrieSaine(self):
  38     """
  39     Paramètres géométriques du tuyau coudé sain:
  40     angleCoude
  41     r_cintr
  42     l_tube_p1
  43     l_tube_p2
  44     epais
  45     de
  46     """
  47     self.geomParams = dict(angleCoude = 90,
  48                            r_cintr    = 1200,
  49                            l_tube_p1  = 1600,
  50                            l_tube_p2  = 1200,
  51                            epais      = 40,
  52                            de         = 760)
  53
  54   # ---------------------------------------------------------------------------
  55   def setParamMaillageSain(self):
  56     self.meshParams = dict(n_long_p1    = 16,
  57                            n_ep         = 3,
  58                            n_long_coude = 15,
  59                            n_circ_g     = 20,
  60                            n_circ_d     = 20,
  61                            n_long_p2    = 12)
  62
  63 # ---------------------------------------------------------------------------
  64   def setParamShapeFissure(self):
  65     """
  66     paramètres de la fissure pour le tuyau coude
  67     profondeur  : 0 < profondeur <= épaisseur
  68     rayonPipe   : rayon du pipe correspondant au maillage rayonnant
  69     lenSegPipe  : longueur des mailles rayonnantes le long du fond de fissure (= rayonPipe par défaut)
  70     azimut      : entre 0 et 360°
  71     alpha       : 0 < alpha < angleCoude
  72     longueur    : <=2*profondeur ==> force une fissure elliptique (longueur/profondeur = grand axe/petit axe).
  73     orientation : 0° : longitudinale, 90° : circonférentielle, autre : uniquement fissures elliptiques
  74     lgInfluence : distance autour de la shape de fissure a remailler (si 0, pris égal à profondeur. A ajuster selon le maillage)
  75     elliptique  : True : fissure elliptique (longueur/profondeur = grand axe/petit axe); False : fissure longue (fond de fissure de profondeur constante, demi-cercles aux extrémites)
  76     pointIn_x   : optionnel coordonnées x d'un point dans le solide, pas trop loin du centre du fond de fissure (idem y,z)
  77     externe     : True : fissure face externe, False : fissure face interne
  78     """
  79     self.shapeFissureParams = dict(profondeur  = 10,
  80                                    rayonPipe   = 2.5,
  81                                    lenSegPipe  = 2.5,
  82                                    azimut      = 180,
  83                                    alpha       = 40,
  84                                    longueur    = 200,
  85                                    orientation = 0,
  86                                    lgInfluence = 50,
  87                                    elliptique  = False,
  88                                    externe     = False)
  89
  90   # ---------------------------------------------------------------------------
  91   def setParamMaillageFissure(self):
  92     """
  93     Paramètres du maillage de la fissure pour le tuyau coudé
  94     Voir également setParamShapeFissure, paramètres rayonPipe et lenSegPipe.
  95     nbSegRad = nombre de couronnes
  96     nbSegCercle = nombre de secteurs
  97     areteFaceFissure = taille cible de l'arête des triangles en face de fissure.
  98     """
  99     self.maillageFissureParams = dict(nomRep        = os.curdir,
 100                                       nomFicSain    = self.nomProbleme,
 101                                       nomFicFissure = self.nomProbleme + "_fissure",
 102                                       nbsegRad      = 5,
 103                                       nbsegCercle   = 8,
 104                                       areteFaceFissure = 5)
 105
 106   # ---------------------------------------------------------------------------
 107   def setReferencesMaillageFissure(self):
 108     from salome.smesh import smeshBuilder
 109     if smeshBuilder.NETGEN_VERSION_MAJOR < 6:
 110       self.referencesMaillageFissure = dict( \
 111                                              Entity_Quad_Quadrangle = 5366, \
 112                                              Entity_Quad_Hexa = 7120, \
 113                                              Entity_Node = 54832, \
 114                                              Entity_Quad_Edge = 722, \
 115                                              Entity_Quad_Triangle = 1570, \
 116                                              Entity_Quad_Tetra = 11747, \
 117                                              Entity_Quad_Pyramid = 532, \
 118                                              Entity_Quad_Penta = 496 \
 119       )
 120     else:
 121       self.referencesMaillageFissure = dict( \
 122                                              Entity_Quad_Quadrangle = 5366, \
 123                                              Entity_Quad_Hexa = 7120, \
 124                                              Entity_Node = 54175, \
 125                                              Entity_Quad_Edge = 721, \
 126                                              Entity_Quad_Triangle = 1482, \
 127                                              Entity_Quad_Tetra = 11324, \
 128                                              Entity_Quad_Pyramid = 532, \
 129                                              Entity_Quad_Penta = 496 \
 130       )