Salome HOME
Copyright update 2022
[modules/smesh.git] / src / Tools / blocFissure / CasTests / fissureCoude_6.py
1 # -*- coding: utf-8 -*-
2 # Copyright (C) 2014-2022  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 - ASCOU08"""
21
22 import os
23
24 from blocFissure.gmu.fissureCoude  import fissureCoude
25
26 class fissureCoude_6(fissureCoude):
27   """ASCOU08"""
28
29   nomProbleme = "fissureCoude_6"
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 = 40,
48                            r_cintr    = 854,
49                            l_tube_p1  = 1700,
50                            l_tube_p2  = 1700,
51                            epais      = 62.5,
52                            de         = 912.4)
53
54   # ---------------------------------------------------------------------------
55   def setParamMaillageSain(self):
56     self.meshParams = dict(n_long_p1    = 15,
57                            n_ep         = 3,
58                            n_long_coude = 16,
59                            n_circ_g     = 30,
60                            n_circ_d     = 30,
61                            n_long_p2    = 15)
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   = 4,
81                                    lenSegPipe  = 4,
82                                    azimut      = 90,
83                                    alpha       = 20,
84                                    longueur    = 133,
85                                    orientation = 0,
86                                    lgInfluence = 30,
87                                    elliptique  = False,
88                                    externe     = True)
89
90   # ---------------------------------------------------------------------------
91
92   def setParamMaillageFissure(self):
93     """
94     Paramètres du maillage de la fissure pour le tuyau coudé
95     Voir également setParamShapeFissure, paramètres rayonPipe et lenSegPipe.
96     nbSegRad = nombre de couronnes
97     nbSegCercle = nombre de secteurs
98     areteFaceFissure = taille cible de l'arête des triangles en face de fissure.
99     """
100     self.maillageFissureParams = dict(nomRep        = os.curdir,
101                                       nomFicSain    = self.nomProbleme,
102                                       nomFicFissure = self.nomProbleme + "_fissure",
103                                       nbsegRad      = 5,
104                                       nbsegCercle   = 12,
105                                       areteFaceFissure = 5)
106
107   # ---------------------------------------------------------------------------
108   def setReferencesMaillageFissure(self):
109     from salome.smesh import smeshBuilder
110     if smeshBuilder.NETGEN_VERSION_MAJOR < 6:
111       self.referencesMaillageFissure = dict( \
112                                              Entity_Quad_Quadrangle = 7392, \
113                                              Entity_Quad_Hexa = 9888, \
114                                              Entity_Node = 61217, \
115                                              Entity_Quad_Edge = 833, \
116                                              Entity_Quad_Triangle = 1024, \
117                                              Entity_Quad_Tetra = 6725, \
118                                              Entity_Quad_Pyramid = 452, \
119                                              Entity_Quad_Penta = 408 \
120                                            )
121     else:
122       self.referencesMaillageFissure = dict( \
123                                              Entity_Quad_Quadrangle = 7392, \
124                                              Entity_Quad_Hexa = 9888, \
125                                              Entity_Node = 61307, \
126                                              Entity_Quad_Edge = 833, \
127                                              Entity_Quad_Triangle = 1014, \
128                                              Entity_Quad_Tetra = 6800, \
129                                              Entity_Quad_Pyramid = 452, \
130                                              Entity_Quad_Penta = 408 \
131       )