Code_Aster ®
Version
8.2
Titre :
DEMO006 Optimisation du rayon de courbure d'une tuyauterie coudée
Date :
05/09/05
Auteur(s) :
J. LAVERNE, F. LEBOUVIER
Clé : V3.03.131-A
Page :
1/6
Organisme(s) : EDF-R&D/AMA, DeltaCAD
Manuel de Validation
Fascicule V3.03 : Statique linéaire des plaques et coques
Document V3.03.131
DEMO006 Optimisation du rayon de courbure
d'une tuyauterie coudée
Résumé :
Ce test de démonstration illustre l'utilisation du langage python dans Code_Aster.
L'objectif de ce cas-test est d'optimiser la valeur du rayon de courbure d'une tuyauterie via une boucle python.
Le rayon de courbure est modifié itérativement. La valeur optimale est obtenue lorsque la contrainte maximum
de Von Mises est inférieure à un seuil.
Le langage python, dans ce cas-test permet itérativement :
·
de modifier le fichier décrivant la géométrie du coude,
·
de lancer GMSH afin de mailler le coude,
·
de récupérer la contrainte maximum de Von Mises,
·
d'évaluer un critère d'arrêt des itérations,
·
de lancer GMSH pour effectuer un post-traitement interactif.
Ce cas-test reprend le problème du cas-test forma01f : tuyauterie coudée, constituée d'un matériau élastique
linéaire, soumise à une force appliquée à son extrémité, modélisée par des éléments de coques DKT.
Remarque :
Pour avoir les fonctions interactives de visualisation, mettre interactif = 1 au début du fichier de
commande.
Manuel de Validation
Fascicule V3.03 : Statique linéaire des plaques et coques
HT-66/05/005/A
Code_Aster ®
Version
8.2
Titre :
DEMO006 Optimisation du rayon de courbure d'une tuyauterie coudée
Date :
05/09/05
Auteur(s) :
J. LAVERNE, F. LEBOUVIER
Clé : V3.03.131-A
Page :
2/6
1
Problème de référence
1.1 Géométrie
L'étude concerne une tuyauterie comprenant deux tuyaux droits et un coude [Figure 1.1-a].
Les données géométriques du problème sont les suivantes :
·
la longueur LG des deux tuyaux droits est de 3 m,
·
le rayon Rc initial du coude est de 0.3 m,
·
l'angle du coude est de 90 degrés,
·
l'épaisseur des tuyaux droits et du coude est de 0.02 m,
·
et le rayon extérieur Re des tuyaux droits et du coude est de 0.2 m.
LG
D
B
section D
section B
RC
C
O
section C
Z
Y
e
L
Z
G
X
Re
X
A
section A
Figure 1.1-a
Remarque :
La géométrie du problème présente une symétrie par rapport au plan (A,X,Y).
1.2
Propriétés de matériaux
Matériau élastique linéaire isotrope. Les propriétés du matériau sont celles de l'acier A42 :
·
le module d'Young :
E = 1.8 1011 Pa,
·
le coefficient de Poisson :
= 0.3,
Manuel de Validation
Fascicule V3.03 : Statique linéaire des plaques et coques
HT-66/05/005/A
Code_Aster ®
Version
8.2
Titre :
DEMO006 Optimisation du rayon de courbure d'une tuyauterie coudée
Date :
05/09/05
Auteur(s) :
J. LAVERNE, F. LEBOUVIER
Clé : V3.03.131-A
Page :
3/6
1.3
Conditions aux limites et chargements
· Conditions aux limites : encastrement au niveau de la section A,
· Chargement : force constante FY dirigée selon l'axe Y et appliquée sur la section B.
La valeur de FY est calculée à partir :
du rayon moyen : RMOY = 0.19 ,
de la force totale appliquée :
2
FTOT = 500000 N / m ,
Son expression est la suivante :
FY = FTOT (2 RMOY ) ( 418828.8)
· La contrainte de Von Mises limite est de 2.0E+09 N/m²
2
Solution de référence
2.1
Méthode de calcul utilisée pour la solution de référence
La solution de référence est obtenue numériquement. Il s'agit donc uniquement d'un test de non
régression.
2.2
Résultats de la modélisation
Le valeur optimale du rayon de courbure, respectant le critère mises < 2.109 N.m-2 est de 1.1m,
obtenue après 5 itérations. A chacune des itérations le rayon de courbure est augmenté de 0.2m.
Rayons de
Contrainte max de
Itérations
courbure
Von Mises
1
0.3 m
3.3315 E+09 N/m²
2
0.5 m
2.7647 E+09 N/m²
3
0.7 m
2.4256 E+09 N/m²
4
0.9 m
2.1727 E+09 N/m²
5
1.1 m
1.9670 E+09 N/m²
2.3
Incertitude sur la solution
Solution de non régression.
Manuel de Validation
Fascicule V3.03 : Statique linéaire des plaques et coques
HT-66/05/005/A
Code_Aster ®
Version
8.2
Titre :
DEMO006 Optimisation du rayon de courbure d'une tuyauterie coudée
Date :
05/09/05
Auteur(s) :
J. LAVERNE, F. LEBOUVIER
Clé : V3.03.131-A
Page :
4/6
3 Modélisation
A
3.1
Caractéristiques de la modélisation
Modélisation en éléments coques (DKT).
Figure 3.1-a : géométrie et maillage du demi cylindre (RC initial)
3.2
Caractéristiques du maillage
Le maillage est régénéré à chaque itération en raison de la modification du rayon de courbure de la
géométrie. Cependant les caractéristiques topologiques du maillage sont inchangées :
·
1013 mailles (900 TRIA3, 110 SEG2, 3POI3)
·
507 noeuds.
Les groupes de mailles correspondent à :
·
GM30 <=> surface du TUYAU
·
GM28 <=> section B (effort)
·
GM31 <=> point A1 (-R, 0, 0)
·
GM27 <=> section A (encastrement)
·
GM29 <=> SYMETRIE
Manuel de Validation
Fascicule V3.03 : Statique linéaire des plaques et coques
HT-66/05/005/A
Code_Aster ®
Version
8.2
Titre :
DEMO006 Optimisation du rayon de courbure d'une tuyauterie coudée
Date :
05/09/05
Auteur(s) :
J. LAVERNE, F. LEBOUVIER
Clé : V3.03.131-A
Page :
5/6
3.3 Commandes
Aster
Ce paragraphe décrit l'algorithme utilisé dans le fichier de commandes et présente les commandes
Code_Aster utilisées.
·
Initialisation de certaines variables :
Rayon de courbure : Rc = 0.3
Critère de convergence : crit = 2.0E+09 (contrainte de Von Mises cible)
·
On entre dans la boucle python dont le critère d'arrêt porte sur la contrainte max de Von
Mises. Tant que le critère n'est pas respecté, on effectue les instructions suivantes :
Reconstruction du maillage :
on
redéfinit
Rc dans le fichier géométrique .geo de GMSH
on lance gmsh via python pour générer le fichier de maillage .msh
Lecture du maillage (PRE_GMSH) et génération du maillage (LIRE_MAILLAGE). On utilise
DEFI_GROUP pour renommer les groupes de mailles selon la correspondance :
# GM30 <=> TUYAU
# GM28 <=> EFOND
# GM31 <=> A1
# GM27 <=> ENCAST
#
GM29
<=>
SYMETRIE
Définition des éléments finis utilisés (AFFE_MODELE). Les tuyaux droits et le coude sont
modélisés par des éléments de coque (DKT).
Réorientations des normales aux éléments : on utilise MODI_MAILLAGE pour orienter tous
les éléments de la même façon, avec une normale tournée vers l'intérieur.
Définition et affectation du matériau (DEFI_MATERIAU et AFFE_MATERIAU). Les
caractéristiques mécaniques sont identiques sur toute la structure.
Affectation des caractéristiques des éléments coques (AFFE_CARA_ELEM) : épaisseur,
vecteur V définissant le repère de dépouillement (mot-clé ANGL_REP)
Définition des conditions aux limites et du chargement (AFFE_CHAR_MECA).
La tuyauterie est encastrée en sa base, sur tous les noeuds situés dans le plan Y=0.
La tuyauterie présente un plan de symétrie Z=0.
On calcule un effort réparti FY dirigé selon l'axe Y et appliquée à la section B, (l'effort
réparti est tel que la résultante 2Pi* RMOY. FY = FTOT, FTOT étant la force totale
que l'on désire appliquer). Pour appliquer l'effort sur la section B, on utilisera
FORCE_ARETE.
Résolution du problème élastique linéaire (MECA_STATIQUE).
Calcul du champ de contraintes par éléments aux noeuds pour chaque cas de charge
(option 'SIGM_ELNO_DEPL'). Les contraintes sont calculées dans le repère local défini
pour chaque élément à l'aide du vecteur V (mot clé ANG_REP précédent). Utiliser
NIVE_COUCHE pour définir le niveau de calcul dans l'épaisseur.
Calcul du champ de contraintes équivalentes par éléments aux noeuds calculé à partir du
champ de contrainte (option 'EQUI_ELNO_SIGM').
Calcul des champs précédents aux noeuds (options 'SIGM_NOEU_DEPL',
Boucle python d'optimisation du rayon de courbure
'EQUI_NOEU_SIGM')
Impression des résultats (IMPR_RESU).
Détermination d'une table contenant les calculs de moyennes du champ de contraintes
équivalentes aux noeuds.(`POST_RELEVE_T')
Extraction de la composante VMIS de la table précédente via python.
Test d'arrêt :
Si VMIS est supérieure à CRIT, alors on réitère avec Rc=Rc+0.2
Sinon on sort de la boucle python.
Manuel de Validation
Fascicule V3.03 : Statique linéaire des plaques et coques
HT-66/05/005/A
Code_Aster ®
Version
8.2
Titre :
DEMO006 Optimisation du rayon de courbure d'une tuyauterie coudée
Date :
05/09/05
Auteur(s) :
J. LAVERNE, F. LEBOUVIER
Clé : V3.03.131-A
Page :
6/6
4
Résultats de la modélisation A
4.1 Valeurs
testées
Rayons de
Contrainte max de Von
Itérations
courbure
Mises
1
0.3 m
3.3315E+09 N/m²
2
0.5 m
2.7647E+09 N/m²
3
0.7 m
2.4256E+09 N/m²
4
0.9 m
2.1727E+09 N/m²
5
1.1 m
1.9670E+09 N/m²
5
Synthèse des résultats
Le rayon de courbure optimal respectant le critère : contrainte max de Von Mises < à 2.0E+09 est de
1.1m. Il a été obtenu après 5 itérations, et la contrainte maximum de Von Mises trouvée est de
1.9670E+09 N/m².
Manuel de Validation
Fascicule V3.03 : Statique linéaire des plaques et coques
HT-66/05/005/A
Document Outline