Modélisation C ============== Dans le cas de la modélisation en éléments coques, le maillage consiste en la discrétisation de la **surface moyenne** de la tuyauterie. La géométrie étant symétrique par rapport au plan :math:`(A,X,Y)` , on ne maille qu’une demi-surface. Conditions aux limites et chargement: encastrement aux deux extrémités de la tuyauterie, et pression à la surface interne. .. image:: images/100000000000033C000001D7B34283E122499990.png :width: 5.5654in :height: 3.2319in .. _RefImage_100000000000033C000001D7B34283E122499990.png: A1 P Géométrie --------- On peut créer cette géométrie en définissant les points :math:`A1` , :math:`P` et :math:`\mathrm{A2}` , puis l’arc de cercle :math:`\mathit{Base}` (Menu New Entity → Basic → Point / Arc).Il suffit ensuite de créer le premier tuyau droit :math:`\mathrm{AC}` à partir de l’arc de cercle :math:`\mathrm{Base}` par le menu New Entity → Generation → Extrusion: Vector = OY, Height = 3. Pour créer le coude, il faut récupérer l’extrémité du tuyau :math:`\mathrm{AC}` en appliquant le menuNewEntity→Explode (Edge), puis créer un vecteur parallèle à l’axe Z et passant le point O :math:`(\mathrm{0.6,}\mathrm{3.0,}0.)` . Ensuite générer la géométrie du coude par le menu New Entity → Generation → Revolution. Enfin, appliquer la même démarche pour le tuyau :math:`\mathrm{DB}` (Explode puis Extrusion). Créer un "compound" (Menu New Entity → Build → Compound) en sélectionnant les trois parties de la tuyauterie. On créera ensuite les groupes desmailles où on désire mettre des conditions limites: Base, Symetrie, Efond et surface du tuyau (Menu New Entity → Group → Create Group). On créera aussi le groupe point :math:`\mathrm{A1}` . Maillage -------- Lancer le module Mesh de la plate-forme Salome-Meca. Le maillage est défini par le menu Mesh → Create Mesh. Sélectionner la géométrie à mailler, puis l’algorithme et l'hypothèse de discrétisation par dimension: * 2D Quadrangle: Mapping. * 1D Wire Discretisation avec l’hypothèse basique Number of Segment(15 segments par edge. Puis calculer le maillage(Menu Mesh → Compute). Pour permettre un raffinement différent selon les edges, on créera un sous-maillage (Menu Mesh → Create Sub-mesh) définissant l’hypothèse basique Number of Segment sur la circonférence, par exemple 10 segments sur :math:`\mathit{base}` et l’hypothèse supplémentaire «Propagation of 1D hypothese on Opposite Edges». Puis calculer le maillage(Menu Mesh → Compute). Créer les groupes demailles correspondants aux groupe géométriques(Menu Mesh→Create Groups from Geometry). Exporter lemaillage au format MED. Création et lancement du cas de calcul (via AsterStudy) ------------------------------------------------------- Lancer le module AsterStudyde la plate-forme Salome-Meca. Puis en colonne gauche, cliquer surl'onglet Case View. On définit le fichier de commandes du cas de calcul. Nota:ajouter des commandes par le menu Commands→ Show All. Les principales étapes de ce calcul mécanique pour la création et le lancement du cas de calcul sont les suivantes: * Lire le maillage au format MED: Commande LIRE_MAILLAGE. * Définir les éléments finis utilisés: Commande AFFE_MODELE. La tuyauterie seramodélisée par des éléments de coque (DKT). * Orienter des normales aux éléments: Commande MODI_MAILLAGE / ORIE_NORM_COQUE pour orienter tous les éléments de la même façon, avec une normale tournée vers l’intérieur du tuyau (étant donné la convention de signe sur la pression) afin de donner une valeur positive à la pression (utiliser le groupe surface). * Définir le matériau: Commande DEFI_MATERIAU: E, NU, ALPHA * Affecter le matériau: Commande AFFE_MATERIAU. Les caractéristiques mécaniques sont identiques sur toute la structure. * Affecterles caractéristiques des éléments coques: Commande AFFE_CARA_ELEM / COQUEpour définir l’épaisseur. * Affecterles conditions aux limites mécaniques et les chargements: Commande AFFE_CHAR_MECA: * Il y a un encastrement sur le groupe de mailles :math:`\mathrm{Base}` et :math:`\mathrm{Efond}` , et des conditions de symétrie (déplacement normal :math:`\mathrm{DZ}` nul et rotations :math:`\mathrm{DRX}` et :math:`\mathrm{DRY}` nulles) sur le groupe de mailles :math:`\mathrm{Symetrie}` : DDL_IMPO. * Une pression interne :math:`P` : PRES_REP. Il faut convertir la pression P à la surface interne en la pression à la surface moyenne. * Résoudre le problème élastique : Commande MECA_STATIQUE: CARA_ELEM, CHAM_MATER, EXCIT, MODELE. * Imprimer les déplacements et les contraintesau format MED: Commande IMPR_RESU.. * Pour lancerle cas de calcul, en colonne gauche, cliquer sur l'onglet History View.