smesh/src/Tools/blocFissure/CasTests/fissureCoude_1.py

118 lines
5.4 KiB
Python
Raw Normal View History

2014-01-09 14:20:44 +00:00
# -*- coding: utf-8 -*-
2021-03-23 17:44:27 +03:00
# Copyright (C) 2014-2021 EDF R&D
2019-02-14 14:55:47 +03:00
#
# This library is free software; you can redistribute it and/or
# modify it under the terms of the GNU Lesser General Public
# License as published by the Free Software Foundation; either
# version 2.1 of the License, or (at your option) any later version.
#
# This library is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
# Lesser General Public License for more details.
#
# You should have received a copy of the GNU Lesser General Public
# License along with this library; if not, write to the Free Software
# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
#
# See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
#
2021-03-17 14:13:23 +01:00
"""Fissure dans un coude"""
2014-01-09 14:20:44 +00:00
2021-01-26 15:12:35 +01:00
import os
2014-01-09 14:20:44 +00:00
from blocFissure.gmu.fissureCoude import fissureCoude
class fissureCoude_1(fissureCoude):
2021-01-21 16:09:37 +01:00
"""problème de fissure du Coude : adaptation maillage"""
nomProbleme = "fissureCoude_1"
2021-03-17 14:13:23 +01:00
geomParams = dict()
meshParams = dict()
shapeFissureParams = dict()
maillageFissureParams = dict()
referencesMaillageFissure = dict()
2014-01-09 14:20:44 +00:00
# ---------------------------------------------------------------------------
def setParamGeometrieSaine(self):
"""
Paramètres géométriques du tuyau coudé sain:
angleCoude
r_cintr
l_tube_p1
l_tube_p2
epais
de
"""
self.geomParams = dict(angleCoude = 60,
r_cintr = 1200,
l_tube_p1 = 1600,
l_tube_p2 = 1200,
epais = 40,
de = 760)
# ---------------------------------------------------------------------------
def setParamMaillageSain(self):
self.meshParams = dict(n_long_p1 = 16,
n_ep = 3,
n_long_coude = 15,
n_circ_g = 20,
n_circ_d = 20,
n_long_p2 = 12)
# ---------------------------------------------------------------------------
def setParamShapeFissure(self):
"""
paramètres de la fissure pour le tuyau coude
profondeur : 0 < profondeur <= épaisseur
rayonPipe : rayon du pipe correspondant au maillage rayonnant
lenSegPipe : longueur des mailles rayonnantes le long du fond de fissure (= rayonPipe par défaut)
azimut : entre 0 et 360°
alpha : 0 < alpha < angleCoude
longueur : <=2*profondeur ==> force une fissure elliptique (longueur/profondeur = grand axe/petit axe).
orientation : 0° : longitudinale, 90° : circonférentielle, autre : uniquement fissures elliptiques
lgInfluence : distance autour de la shape de fissure a remailler (si 0, pris égal à profondeur. A ajuster selon le maillage)
elliptique : True : fissure elliptique (longueur/profondeur = grand axe/petit axe); False : fissure longue (fond de fissure de profondeur constante, demi-cercles aux extrémites)
pointIn_x : optionnel coordonnées x d'un point dans le solide, pas trop loin du centre du fond de fissure (idem y,z)
externe : True : fissure face externe, False : fissure face interne
"""
self.shapeFissureParams = dict(profondeur = 10,
rayonPipe = 2.5,
lenSegPipe = 2.5,
azimut = 160,
alpha = 20,
longueur = 400,
orientation = 90,
lgInfluence = 50,
elliptique = False,
externe = True)
# ---------------------------------------------------------------------------
def setParamMaillageFissure(self):
"""
Paramètres du maillage de la fissure pour le tuyau coudé
Voir également setParamShapeFissure, paramètres rayonPipe et lenSegPipe.
nbSegRad = nombre de couronnes
nbSegCercle = nombre de secteurs
areteFaceFissure = taille cible de l'arête des triangles en face de fissure.
"""
2021-01-26 15:12:35 +01:00
self.maillageFissureParams = dict(nomRep = os.curdir,
2014-01-09 14:20:44 +00:00
nomFicSain = self.nomCas,
2021-03-03 18:47:59 +01:00
nomFicFissure = self.nomCas + "_fissure",
2014-01-09 14:20:44 +00:00
nbsegRad = 5,
nbsegCercle = 6,
areteFaceFissure = 5)
# ---------------------------------------------------------------------------
def setReferencesMaillageFissure(self):
2021-01-21 10:07:24 +01:00
self.referencesMaillageFissure = dict( \
Entity_Quad_Quadrangle = 6042, \
Entity_Quad_Hexa = 7794, \
Entity_Node = 65337, \
Entity_Quad_Edge = 874, \
Entity_Quad_Triangle = 2058, \
Entity_Quad_Tetra = 16037, \
Entity_Quad_Pyramid = 738, \
Entity_Quad_Penta = 672 \
)