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