TP-Cast3M

TP-treillis
L'objectif de ce TP est de comprendre la structure et le langage d'une suite d'instructions GIBIANE (procédure) pour vous permettre de débuter
rapidement avec Cast3M.
Documents associés
Cours et exercices de l'ECN [ref 1] : [Link]
Méthode des éléments finis - MEF - Cours ECN - H. OUDIN
Fichier Cast3M de l'exemple traité dans le chapitre 2 du cours : [Link]
Fichier Cast3M de la colonne traitée dans le chapitre 2 du cours : [Link]
Documentation Cast3M ([Link]
[ref 2] : Maillage - CEA - Utiliser CAst3M - F. DI PAOLA
[ref 3] : Post-traitement - CEA - Utiliser CAst3M - F. DI PAOLA

Treillis étudié c'est l'exemple traité dans le chapitre de cours sur l'étude des treillis par la MEF [ref 1]
F Nous cherchons la réponse statique du treillis plan représenté par la figure ci-contre. Ce treillis est un
assemblage (par rotules) d’éléments travaillant en traction – compression.
a Pour les calculs numériques nous utiliserons les données suivantes : a = 1 m et F = 30 KN
Barres sont en acier de section circulaire pleine identique, diamètre φ = 2 cm
a 2 Pour l'acier E = 210 GPa , Re = 240 MPa
yo
Pour le dimensionnement la contrainte maximale ne doit pas dépasser Re / 2 et la structure doit rester
xo stable (étude au flambement).

Première étape
Ouvrez le fichier [Link] dans votre éditeur et identifier les directives et opérateurs, relatifs aux 3 grandes phases de l'analyse par la méthode des
éléments que vous avez vu en cours :
• Analyse du problème Définition du modèle mathématique (type d'analyse, type d'éléments)
Discrétisation géométrique du domaine étudié
TP-Cast3M

Propriétés Mécaniques (matériau, section)

• Résolution du problème Calcul des matrices élémentaires, assemblage de la structure complète
Application des chargements externes et des conditions aux limites
Résolution du système d'équations
• post-traitement Déformée de la structure et efforts aux appuis
Quantités locales : contraintes, déformations
Pour les directives OPTI , TRAC, MESS, LIST faites le lien entre la documentation dans Notices du site du CEA et leur utilisation pour ce calcul.
De même pour les opérateurs MODE, MATE, CARA, RIGI, RESO
Ce travail est une première étape pour comprendre la philosophie du langage objet utilisé et la façon dont elle est présentée dans la documentation. Vous
venez d'ouvrir une "boite de Pandore" une fois ces premiers calculs effectués vous retournerez très rapidement sur le site du CEA pour approfondir vos
connaissances et aller plus loin avec ce code, si vous le souhaitez.
Exercice 1
Exécutez le calcul et observez ce qui se passe dans les deux fenêtres "Dos" et "graphique"
Pour terminer le calcul vous devez fermer la fenêtre Dos
Supprimer l'astérisque de la dernière ligne la directive fin devient active, relancer le calcul
La fenêtre Dos se ferme automatiquement, vous ne pouvez plus lire les résultats du calcul.
Passez la variable listing à O pour oui, relancer le calcul et ouvrer le fichier [Link] qui doit être dans le répertoire ou vous travaillez.
Vous pouvez revenir à la configuration initiale du fichier.
Supprimer l'astérisque de la ligne OPTI ECHO 0; qui devient active et relancer le calcul.
Observez la fenêtre "Dos"
Ces petits exercices vous permettent de choisir la configuration graphique qui vous convient le mieux en fonction de ce que vous avez à faire.
Bien entendu vous avez pris le temps de comparer les résultats obtenu lors de ces calculs avec ceux présentés dans le polycopié de cours, et d'en tirer les
conclusions qui s'imposent.
Exercice 2
Modifier le maillage en utilisant 4 éléments pour discrétiser la barre 1-3
Regardez et expliquez ce qui c'est passé lors de ce calcul, revenez au maillage initial.
Prenez le temps de vous renseigner sur les autres opérateurs et directives utilisés dans ce fichier [Link]. Concernant
le maillage DROI, profitez en pour regarder la documentation [Link] [ref 2],
les conditions aux limites BLOQ et FORC,

2
TP-Cast3M

le post-traitement DEFO, REAC, SIGMA,

et l'extraction d'informations utiles pouvant être réutilisées EXTR, EXCO.
Vous venez de faire vos premiers pas avec CAst3M. Vous savez exécuter des instructions, et vous déplacer dans les 3 fenêtres : Dos, Graphique, Éditeur.
Nous apprendrons un peu plus tard à manipuler la fenêtre de tracé en allant un peu plus loin sur l'exploitation et la mise en forme des résultats d'un calcul,
mais si cela vous intéresse vous pouvez diagonaliser le document [Link] [ref 3]

Deuxième étape
Vous allez maintenant modifier plus sérieusement le fichier pour programmer un calcul de dimensionnement de la structure. Ce travail vous amènera à
utiliser par vous même certains opérateurs du code.
L'objectif est de fournir un fichier dgibi autonome qui vérifie que la contrainte maximale ne dépasse pas Re/2 et que la structure reste stable dans le
domaine élastique. Vous mettrez le diamètre de la section circulaire en paramètre et les résultats suivants devront être affichés :
déplacement du point de chargement,
contrainte maximale pour la charge appliquée
charge maximale autorisée pour rester dans le domaine élastique
y-a-t-il ruine ou est-on en sécurité par rapport à la limite élastique ?
Exercice 3
Programmez le calcul de la contrainte maximale en utilisant les résultats du calcul, comparez à la valeur théorique.
Que pensez-vous de cette contrainte maximale.
Programmez le calcul de la charge limite élastique.
La structure est-elle stable pour une charge respectant le critère de dimensionnement élastique ?
Rappel la charge critique d'Euler est donnée par π 2 EI / L2c
Quelle est la charge maximale que peut supporter cette structure ?
Comment sont modifiés les résultats pour une section de diamètre φ = 4 cm ?
Fournir un fichier "dgibi" fonctionnant automatiquement pour ce treillis.
Vous testerez votre programme avec un diamètre de 2 puis de 4 cm

Troisième étape
Ouvrez le fichier [Link] dans votre éditeur, vous pouvez lancer le calcul.

3
TP-Cast3M

Que c'est-il passé et quelles sont les nouveautés de ce fichier ?

L'intérêt ici est de pouvoir comparer les résultats du calcul à l'étude théorique présenté dans le cours [ref 1].
Pour cela on a utilisé les opérateurs PROG, DIME, EVOL, MANU pour créer les courbes de la solution analytique que l'on peut tracer avec DESS
Regardez la notice de ces instructions, vous pouvez aussi lire le chapitre 3 de la [ref 3].
Regardez et comprenez la façon dont le chargement du poids propre a été introduit
Utilisation de l’opérateur MASS, pour calculer les matrices Masse élémentaires de la structure
Puis création manuelle, opérateur MANU d’un CHPOINT contenant l’accélération g appliquée aux nœuds dans la direction x

Pour retrouver les courbes analytiques à partir des résultats de calcul il faut savoir à quel type de champ on a à faire, et créer si nécessaire un champ de
valeurs par point (CHPOINT) sur la ligne où l'on compare les résultats
Pour cela on a utilisé l'opérateur CHAN qui permet de transformer le champ des contraintes (MCHAML) par élément en un CHPOINT aux nœuds.
Attention cet opérateur a calculé par défaut la valeur moyenne des contraintes aux nœuds, ce qui explique que ces valeurs correspondent à la
solution analytique alors que la valeur de la contrainte sur les éléments est différente (effort normal par élément)

Nous ne proposons pas d'autre TP sur l'étude des treillis, les problèmes de stabilité (Flambement) et de vibrations seront abordés avec l’étude des
portiques.

Dernière étape
Exercice 4 Vous allez utiliser et transformer un de vos fichiers pour calculer la structure de l'exercice MEF 5 [ref 1]
3m 2m
Pour le dimensionnement la contrainte maximale ne doit pas dépasser Re et la structure doit rester stable,
pour le calcul de la charge critique vous considérerez que les sections sont circulaires.
S1 = 25cm 2 Les barres sont en acier : E = 210 GPa , Re = 350 MPa
1.5 m

S 2 = 49cm 2 S 2 = 49cm2 Fournir un fichier "dgibi" présentant les résultats du calcul et le dimensionnement de ce treillis.
yo
Vous savez maintenant calculer la réponse statique d'un treillis plan, vous pourrez passer à des treillis
xo F=950KN tridimensionnels facilement.