Ansys Calcul

Computational approach of systems design in engineering
Calcul Scientifique pour la Conception de Systèmes, Partie 2 (CSCS2)
D. Lemosse1,2
1
Département de Mécanique, INSA Rouen Normandie
2
Laboratoire de Mécanique de Normandie, EA 3828
Mai 
D. Lemosse (INSA de Rouen) Computational approach of systems design in engineering Mai  1 / 134
Outline
1 Généralités
2 Ansys
3 Modèles éléments finis en mécanique des solides
4 Dimensionnement mécanique
5 Modèles de résistance des matériaux
6 Génération de maillages
7 Problèmes non linéaires
9 Analyse modale
10 Flambement
11 Résolution par méthodes de puissance
Généralités
1 Généralités
2 Ansys
9 Analyse modale
10 Flambement
Généralités Place de la simulation dans la chaine de conception mécanique
Place de la simulation dans la chaine de conception mécanique
Processus « classique » de conception d’un produit
Coût induit par une modification selon les phases
Objectifs de la simulation en conception mécanique
Évaluation des efforts dans une structure

Les systèmes technologiques sont souvent hyperstatiques avec des configurations
géométriques et mécaniques complexes.
Les méthodes analytiques ne permettent pas l’évaluation des efforts internes.
Vérifier la tenue d’un prototype

Vérifier que la modélisation de la conception de la structure réelle vérifie le cahiers de
charges mécanique, à savoir des contraintes inférieures aux limites dimensionnante :
élastique ou rupture.
Optimiser un système (géométrique, matériau)

On cherche à améliorer une conception existante et fonctionnelle généralement pour
en réduire le coût, soit par la diminution de la matière utilisée, soit par l’amélioration
des performances.
Défauts/qualités des méthodes numériques
Les qualités
peu coûteuses (investissement matériel),
tout est maîtrisé (étude de sensibilité, paramétrique),
répétabilité et quasi-infinité des tests,
accès à des informations physiquement inaccessibles.
Les défauts
phase de modélisation très importante (geom, cpt ...),
multiple approximations (maillage, interpolation ...),
le contrôle des paramètres peut influencer le résultat,
expérience très importante (choix du maillage ...).
Equations d’équilibre des milieux continus déformables
La résolution d’un problème de mécanique des milieux déformables est souvent

difficile.
On peut quelques fois obtenir une solution à partir d’hypothèses permettant

l’introduction d’approximations a priori de la solution.
Mais dans le cas général, cette stratégie ne permet pas d’obtenir une solution vérifiant
l’ensemble des équations du système.
Afin d’obtenir une solution malgré tout, on met en place des stratégies de résolution
permettant d’obtenir une solution approchée acceptable (par exemple au sens des
moindres carrés).
Généralités Principe de la méthode des éléments finis
Mécanique des milieux continus
L’équation d’état est écrit en temps et en espace.

~ (σ) + ~fv = ρ~γ
div
à laquelle on ajoute les conditions aux limites (CL) et les conditions initiales (CI).
On approxime les équations différentielles d’abords en espace pour réduire l’EDP en

EDO, puis en temps.
f (U (x, t) , t) = 0 −→ g (Ui , t) = 0 −→ h (Ui , tj ) = 0
(U (x, t)) fonction continue en temps et en espace, i ∈ [1, Ni ] discrétisation en espace,

j ∈ [1, Nj ] discrétisation en temps.
Approximation des équations
On utilise plusieurs stratégies de résolutions numériques des équations différentielles
Résolution des EDO

Variable unique (ou dissociées)
Schéma dédié : différences finies
Résolution des EDP

Plusieurs variables (associées)
Schémas dédiés :
Eléments finis en structure
Volumes finis en fluide
Eléments de frontière en acoustique
Principe de la méthode des éléments finis
Méthode des éléments finis

Méthode de résolution d’équations aux dérivées partielles s’appuyant sur trois outils
mathématiques :
formulation faible (principe des travaux virtuels)
discrétisation
interpolation
On parle classiquement de formulation en déplacement (déplacements introduits
explicitement en inconnues)
Formulation de Galerkin
Champs réel et virtuel ayant même approximation
Formulation isoparamétrique
Géométrie et cinématique ayant la même approximation
Formulation en déplacements d’une barre
On étudiera le cas statique d’une barre à section constante de matériau homogène en

élasticité linéaire.
~ (σ) + ~fv = ~0
div
dσx (x) dx (x) d 2 u (x)
+ fx (x) = E · + fx (x) = E · + fx (x) = 0
dx dx dx 2
d 2 u (x)
Z
δU · E · + fx (x) dV = 0
V dx 2
Formulation faible de barre
δWint + δWext = 0
Z
dδu (x) du (x)
δWint = − ·E · dV
V dx dx
Z Z
δWext = δU · fx (x) dV + δU · Fx (x) dS
V ∂V
Interpolation cinématique sur une barre
Un élément de barre linéaire à deux noeuds induit une interpolation linéaire des
déplacements à partir des valeurs nodales u1 et u2 .
x x
u (x) = 1 − · u1 + · u2 = N1 (x) · u1 + N2 (x) · u2
L L
La formulation de Galerkin induit la même approximation pour le champ virtuel.
δu (x) = N1 (x) · δu1 + N2 (x) · δu2
Approximation des champs de déformation par éléments
du (x) dN1 (x) dN2 (x) h i u

dN1 (x) dN2 (x) 1
= · δu1 + · δu2 = = [B] {Ue }
dx dx dx dx dx u2
dδu (x)
= [B] {δUe }
dx
Matrice de raideur d’une barre
Ecrire le travail virtuel des efforts internes sous forme condensée matricielle fait
apparaitre la matrice élémentaire de raideur [Ke ]
δWint = − {δUe }T [Ke ] {Ue }
Z
[Ke ] = [B]T ES [B] dx
x
Matrice de raideur de l’élément de barre linéaire à deux noeuds

ES 1 −1
[Ke ] =
L −1 1
Généralités Boucle de modélisation
Boucle de modélisation
Modèle de simulation en 4 phases :
Choix de modélisation analytique.
Choix et réalisation du modèle éléments finis.
Résolution numérique.
Analyse de la validité et des résultats de la modélisation.
Boucle de modélisation - Analyse de l’Objet de l’Etude
Etape 1 : Choix de modélisation analytique
Boucle de modélisation - Analyse de l’Objet de l’Etude (suite)
Choix de la pièce (ensemble de pièces) à étudier.

Définir la géométrie qui modélisera la pièce.
Remplacer l’extérieur par des modèles mécaniques (torseurs, limites).
Faire une collection des données quantitatives (matériaux, expériences).
Introduire les hypothèses mécaniques de modélisation : topologie, linéarités ou
non linéarités.
Représenter la géométrie par un modèle numérique (CAO)
Introduire les simplifications géométriques (courbures, sym)
Découpage de la géométrie en zones simples (assemblages, partitions)
Boucle de modélisation - Modélisation par Eléments Finis
Etape 2 : Choix et réalisation du modèle éléments finis
Boucle de modélisation - Modélisation par Eléments Finis (suite)
Réaliser un maillage de la géométrie

placer des noeuds aux points géométriques importants
mailler tout l’espace
raffiner dans les zones sensibles (mesures, formes complexes)
chercher des éléments de bonne qualité (Q4 ou degré 2)
Traduire les conditions aux limites sur la discrétisation
déplacements imposés aux noeuds (restrictions, contraintes)
efforts ponctuels et répartis transformés en efforts généralisés
Les simplifications de type symétrie sont introduites en déplacements imposés
Il faut éliminer les mouvements de corps rigides
La modélisation analytique et géométrique conditionnent la topologie du modèle
éléments finis :
1D – barre, poutre, coque axisymétrique
2D – plan, axisymétrie, coque
3D – modèle massif
A topologie donnée, il existe plusieurs possibilités de modèles en formulation DKT,
Q4G, Wilson ...
support d’approximation (noeuds, degré)
Il existe des modèles sans support phénoménologique physique ressort, rigide, masse
et raideur
Boucle de modélisation - Résolution Numérique
Etape 3 : Résolution numérique
Boucle de modélisation - Analyse des Résultats
Etape 4 : Analyse de la validité et des résultats de la modélisation
Boucle de modélisation - Analyse des Résultats
Motivations du post-traitement
Le post traitement c’est
L’évaluation de la qualité du modèle (convergence et corrélation).
L’analyse des quantités mécaniques traduisant la qualité de la pièce étudiées
(vérification ou dimensionnement).
Résultats évalués
La MEF évalue les déplacements (degré de liberté aux noeuds) et éventuellement
variables internes (dans l’élément).
Il faut préciser la sauvegarde des données internes aux éléments (contraintes,
déformations, énergies de déformations ...)
Validité de la modélisation
Validité de la modélisation (suite)
Types de dimensionnements
Cinématique (rigidité)
assurer une position (antenne)
assurer un déplacement (débattement d’un capteur)
Statique (résistance)
état de contrainte inférieur à un état limite
non plastification (limite élastique pour les matériaux ductiles) - fonctionnement
normal
non rupture (limite à la rupture) - accident
Fatigue
chargement cyclique avec effort modéré (ou non)
tenue ou nombre de cycles
Flambement
structures élancées en état de compression
phénomène de bifurcation - instabilité
Ansys
1 Généralités
2 Ansys
9 Analyse modale
10 Flambement
Ansys Présentation
Présentation de ANSYS
Ansys est un éditeur de logiciel spécialisé dans la simulation numérique. ANSYS

développe, promeut et assure le support de ses logiciels de simulation servant à
prédire le comportement d’un produit dans son environnement.
F IGURE – Par Tmilnthorp sur Wikipédia
L’entreprise a son siège principal à Canonsburg en Pennsylvanie aux USA. Elle

possède plusieurs sièges et filiales en Europe, Asie et Amérique du nord.
Ansys Présentation
Présentation de Ansys
ANSYS commercialise de nombreux logiciels dans une suite ANSYS Softwares.

Les logiciels phares de cette suite sont :
ANSYS Structural
Calcul statique ;
Calcul dynamique (modale, harmonique, temporel, spectral, vibration aléatoire) ;
Calcul transitoire ;
Flambement linéaire et non linéaire ;
Comportements non linéaires (matériau et géométrique).
ANSYS Mechanical
Il possède les capacités d’ANSYS Structural, en y ajoutant ;
Analyse thermique (conduction, rayonnement, radiation, convection)
Couplage thermo-mécanique ;
Calcul acoustique et couplage acoustique-mécanique ;
Couplage électro-mécanique et électro-thermo-mécanique.
Ansys Présentation
Présentation de Ansys
ANSYS commercialise de nombreux logiciels dans une suite ANSYS Softwares.
Les logiciels phares de cette suite sont :
ANSYS CFX et Fluent
Ces logiciels permettent de réaliser des simulations de mécanique de fluides.
ANSYS AutoDyn et LS-DYNA

Ces logiciels permettent de résoudre les problèmes de mécanique avec de grandes
déformations et de grandes vitesses de déformations, telles que pour les crashs ou les
explosions. Ces logiciels utilisent une formulation explicite.
ANSYS Electromagnetics et Ansoft

Ces logiciels permettent de réaliser des simulations de phénomènes
électromagnétiques.
ANSYS multiphysique
Ce logiciel intègre tous les solveurs de calcul implicite avec un couplage
électromagnétique-thermo-mécanique..
Ansys Présentation
Historique
F IGURE – Portrait de J.A. Swanson - Pittsburgh university
En 1963, le Dr. John A. Swanson travaille sur l’analyse des contraintes des
réacteurs nucléaires à Westinghouse Astronuclear Labs à Pittsburgh, c’est le
début de la vulgarisation de la méthode des éléments finis.
Ansys Présentation
Historique
En 1963, le Dr. John A. Swanson travaille sur l’analyse des contraintes des
réacteurs nucléaires à Westinghouse Astronuclear Labs à Pittsburgh, c’est le
début de la vulgarisation de la méthode des éléments finis.
En 1970, John A. Swanson quitta Westinghouse Astronuclear Labs et développe
le logiciel nommé Swanson Analysis Systems, ce logiciel permettait de simuler le
comportement transitoire des composants nucléaires sous chargement thermique
et de pression.
En 1980, ANSYS fonctionne sur un VAX 11-780 mini-ordinateur.
En 1994, Swanson Analysis System devient ANSYS.
En 1995, ANSYS fonctionne sous Windows 95.
En 2005, ANSYS acquiert AutoDyn.
En 2006, ANSYS achète Fluent.
En 2010, ANSYS achète ANSOFT et MAXWELL.
Ansys Présentation
Environnements logiciels
ANSYS possède deux environnements pour exécuter ces solveurs :
ANSYS Classique ou APDL

C’est la première solution logicielle développée. Elle est destinée à la construction de
modèles éléments finis à géométrie simple, facilement réalisable à l’aide d’opérations
géométriques basiques. A partir de cet environnement, l’utilisateur construit
directement un modèle élément finis en utilisant le langage scripté ANSYS Parametric
Design Language (APDL).
ANSYS Workbench
C’est une plateforme qui propose une approche différente dans la construction des
modèles, en réutilisant le code ANSYS initial. Elle est particulièrement adaptée au
traitement de cas à géométrie complexe (nombreux corps de pièces).
Ansys Présentation
Interface de ANSYS APDL

Pour exécuter le logiciel, il faut sélectionner :
Démarrer > Tous les programmes > ANSYS 2019 R3 > Mechanical APDL 2019 R3
Ansys Présentation
Les fichiers de ANSYS

ANSYS écrit plusieurs fichiers lui permettant de stocker les données ou les résultats.
Le nom de ces fichiers est toujours le même, avec une extension qui change. Par
défaut, le nom du fichier est file.xxx. Pour changer le nom du fichier de travail
(jobname), il faut le faire dés le début de l’étude : File > Change Jobname >
"nomDuFichier"
File type File name File format

Log file Jobname.log ASCII
Error file Jobname.err ASCII
Output file Jobname.out ASCII
Database file Jobname.db Binary
Results files :
Structural or coupled Jobname.rst Binary
Thermal Jobname.rth Binary
Magnetic Jobname.rmg Binary
FLOTRAN Jobname.rfl Binary
Load step file Jobname.Sn ASCII
Graphics file Jobname.grph ASCII
Element matrices Jobname.emat Binary
Ansys Présentation
Choix du type d’analyse
Il faut ensuite sélectionner, dans le menu "Preferences", le type d’étude que l’on va
effectuer.
Ansys Présentation
Organisation du menu principal

Le menu de ANSYS est organisé d’une manière chronologique selon la construction
du modèle de simulation. Ainsi, pour réaliser une analyse par éléments finis, nous
devons :
Préparation de la modélisation Preprocessor /prep7
¯
Choix de l’élément fini à utiliser en fonction de l’étude à réaliser.
Affectation des propriétés physiques (section, épaisseurs ...).
Création des matériaux.
Création du modèle géométrique.
Maillage et création des éléments.
Résolution Solution /sol

Application des conditions aux limites et des chargements.
Résolution
Post traitement General Postproc /post1 et TimeHist Postpro /post26

Traitement et affichage des résultats
Ansys Réalisation d’un modélisation élément finis
Création des modèles mécaniques
Une particularité de Ansys est que le modèle mécanique (milieu et matériau) doit être
défini afin de réaliser la discrétisation (le maillage).
Modèle mécanique
Dans le menu Preprocessor
Le modèle de comportement mécanique est défini par Element Type ;
Les propriétés comme l’aire d’une barre sont définis par Real Constants ;
Les sections et épaisseurs sont définis par Sections ;
Les lois matériaux et leurs valeurs sont définis par Material Props ;
Géométrie
Même s’il est possible de créé un maillage ex-nihilo noeud par noeud, et élément par
élément, on a le plus souvent besoin de définir une géométrie support du maillage.
Modèle géométrique : menu Modeling

Ce menu regroupe tous les outils classiques de création d’objets de CAO. Les outils de
création des objets élémentaires (points, lignes, surfaces, volumes). Des outils de
modifications et de transformations (extrusion, duplication, déplacements ...).
Maillage
La discrétisation dans ANSYS ce fait selon les étapes suivantes :
Maillage : menu Meshing

Sélection de la physique du modèle dans Mesh Attributes ;
Définition de la taille des éléments du maillage dans Size Cntrls ;
Opération de discrétisation dans Mesh.
Ansys Cas d’étude et résolution
Cas d’étude
Bien que cela puisse être défini dans le meun Preprocessor, on définit de préférence
le cas d’étude dans le menu Solution.
Cette sélection, réalisée dans le sous menu Analysis Type, donne accès aux
conditions aux limites et aux options de résolutions dédiées.
Conditions aux limites
Define Loads
Dans ce menu, on peut définir les conditions aux limites du problème : déplacements
imposés, forces et moments ponctuels, pressions, températures, inerties et gravité, etc.
Load Step Opts

Dans ce menu, on peut définir plusieurs séries de conditions aux limites. Ces
conditions sont appliquées successivement comme un chargement temporel (mais
quasi statique dans les cas ne prenant pas en compte la dynamique).
Résolution
Avant d’exécuter la simulation, on peut préciser les options de résolution
Solution Controls
Ce menu permet de choisir
Le type de calcul (Linéaire, non linéaire) ;
Le mode de pilotage de la résolution et le nombre de pas de calcul ;
La fréquence de sauvegarde des résultats ;
Les options de résolution (direct ou itératif).
Ansys Analyse des résultats
Affichage des résultats
Une fois que la résolution est terminée, on peut choisir les jeux de résultats à utiliser
pour représenter la déformation de la pièce ou représenter les champs de contraintes
et de déformations.
Ansys Un langage scripté
Un langage scripté
Ce qu’il faut bien comprendre, c’est que les opérations réalisées successivement par
l’interface logicielle peuvent être réalisées directement en ligne de commande dans un
langage de programmation spécifique au logiciel.
!*————————————* !
!* Creation du type element * !
!*————————————* !
ET,1,SHELL181
!*——————————* !
!* Creation du materiau * !
!*——————————* !
MPTEMP„„„„
MPTEMP,1,0
MPDATA,EX,1„2e11
MPDATA,PRXY,1„0.3
Modèles éléments finis en mécanique des solides
1 Généralités
2 Ansys
9 Analyse modale
10 Flambement
Modélisation des solides
Familles de modèles éléments finis
Mécanique des solides déformables

Éléments massifs
Éléments d’élasticité plane (contrainte ou déformation plane)
Éléments axisymétriques
Éléments de barre
Résistance des matériaux

Éléments plaques et coques
Éléments de poutre
Modèles rhéologiques
Eléments ressort
Eléments rigides
Eléments de masses
Mécanique des solides déformables
Directement obtenus à partir des équations de la mécanique des milieux continus.

Eventuellement avec des hypothèses de comportement global (élasticité plane par
exemple).
Les degrés de liberté sont les déplacements liés à des points matériels.
Les efforts sont traduits par un tenseur des contraintes complet ou partiel.
Modèles de mécanique des solides déformables

3D - Eléments massifs
2D - Eléments d’élasticité plane (contraintes ou déformations plane)
2D - Eléments axisymétriques
1D - Eléments de barre
Modèles éléments finis en mécanique des solides Modélisation des pièces massives
Eléments massifs
Topologie 3D
Cinématique déplacements (U, V , W ) aux noeuds.
Statique Tenseur des contraintes dans l’élément (σx , σy , σz , τxy , τxz , τyz ).
Propriétés physiques -
Matériau Relation de comportement et données matériaux.
Supports géométriques tétraèdre, hexaèdre, prisme.
Eléments massifs - exemple
Maillages structurés et partition

Il faut toujours préférer les éléments héxaédriques. L’obtention de maillage structuré.
Ces maillages sont obtenus par transformation trans-finies ou par extrusions. Sur des
géométrie complexe, ce type de maillage nécessite le partitionnement de la géométrie
complète en éléments plus simples.
Modèles éléments finis en mécanique des solides Modélisation particulières
Eléments plans
Topologie 2D
Cinématique déplacements (U, V ).
Statique Tenseur des contraintes dans l’élément (σx , σy , τxy ).
Propriétés physiques Epaisseurs.
Matériau Relation de comportement (hypothèse σp ou p ) et données matériaux.
Supports géométriques triangle, quadrangle (dans le plan O, X ~,Y ~ ).
Eléments plans - exemple
Eléments axisymétriques
Topologie 2D
Cinématique déplacements (Ur , Uz ) aux noeuds.
Statique Tenseur des contraintes dans l’élément (σr , σz , σθ , τrz ).
Matériau Relation de comportement axisymétrique et données matériaux.

Supports géométriques triangle, quadrangle (dans le plan O, ~er , ~ez ).
Eléments de barre
Topologie 1D
Cinématique déplacements (U) aux noeuds.
Statique Tenseur des contraintes dans l’élément (σx ).
Propriétés physiques Aire de la section.
Supports géométriques linéique.
Dimensionnement mécanique
1 Généralités
2 Ansys
9 Analyse modale
10 Flambement
Dimensionnement mécanique
Cinématique (rigidité)
assurer une position (antenne)
assurer un déplacement (débattement d’un capteur)
Statique (résistance)
état de contrainte inférieur à un état limite
non plastification (limite élastique pour les matériaux ductiles) - fonctionnement
normal
non rupture (limite à la rupture) - accident
Fatigue
chargement cyclique avec effort modéré (ou non)
tenue ou nombre de cycles
Flambement
structures élancées en état de compression
phénomène de bifurcation - instabilité
Dimensionnement mécanique Etat limite - contrainte maximale
Contrainte maximale
Limite élastique σ < σlimite

Critère de Von Mises [σ] → σVM < σlimite
1
q
σVM = √ (σ1 − σ2 )2 + (σ2 − σ3 )2 + (σ3 − σ1 )2
2
Matériaux ductiles
Insensible en hydrostatique.
Matériaux métalliques
Dimensionnement mécanique Etat limite - contrainte maximale
Facteurs de sécurité
Prendre en compte les variations des efforts, en particulier la valeur maximale.

Prendre en compte aussi les écarts des caractéristiques matériaux.
σref
Fdim = γF · Fref → [σ] σdim =
γσ
On va donc chercher à vérifier
σ (Fdim ) < σdim
Les coefficients sont classiquement entre 1.5 et 4, mais peuvent aller jusqu’à 10 sur
les engins de levages pouvant mettre en jeu la vie d’usagers. La valeur du coefficient
dépends de la maitrise de son évaluation.
Dimensionnement mécanique Tenue à la fatigue
Fatigue
Cas classique : sollicitations inférieures

aux limites élastique mais avec un nombre
important de charge-décharge (>1000).
Causes : singularités géométriques

(échelle µm) qui cause des concentrations
de contraintes locales qui génèrent des
fissures qui se propagent.
Plusieurs modélisations du phénomène

existent, empiriques et analytiques suivant
le type de sollicitation (en particulier du
ratio contrainte/cycles).
Courbe de Wöhler
Résultats empiriques : courbe de Wöhler

Conditions expérimentales : sollicitations sinusoidale, contrainte homogène sur
éprouvette, matériau homogène, conditions normalisée température état de surface.
Contrainte limite d’endurance
On définit la contrainte σD au dessous de laquelle la durée de vie d’un matériau (en

nombre de cycles) est infini.
Alliages ferreux : σD existe.
Alliages non ferreux : σD n’existe pas, par convention c’est la valeur pour une
durée de vie de 5 · 108 cycles.
Matériau σD
Fontes et aciers moulés 0.4Rm
Aciers laminés 0.5Rm
Aciers haute résistance [0.3Rm , 0.4Rm ]
Alliages légers coulés 0.3Rm
Alliages légers 0.4Rm
Alliages de titane [0.3Rm , 0.6Rm ]
Lois analytiques empiriques de fatigue
Loi de Basquin N · S m = C → log (N) = log (C) − m · log (S)


S

log
 0.9Rm 
Fer log (N) = 3 
1 + 
σD 
log
0.9Rm

S

log
 0.9Rm 
1 + 1.9
Non ferreux log (N) = 3  
σD 
log
0.9Rm
Histoire du chargement
Pas un sinus parfait, niveaux de contrainte différents au cours du temps
Sur chaque zone on définit :

σmin σmin + σmax
R= et σ= =
σmax 2
Modification de la courbe de Wöhler
En fonction de l’historique de sollicitation, on modifie la courbe de durée de vie

potentielle.
Des lois empiriques de durée de vie sont Etat de la structure

proposées :
σD (R) σm R σm Etat
Soderberg + =1
σD (−1) Re >0 >0 Traction
σD (R) σm >0 <0 Compression
Goodman + =1 <0 tous Alternée
σD (−1) Rm
Vérification de la tenue en fatigue
Chaque cycle d’amplitude σai de Ni cycles est associé à une valeur Ni0 de cycles
admissibles.
Ni
On calcul l’endommagement par fatigue de chaque cycle : ei = 0
Ni
L’endommagement sur la durée de vie complète est : e = Σei
On conclue alors :
e < 1 Pas d’endommagement par fatigue.
e ≥ 1 Fissuration pouvant induire la rupture par fatigue.
Modèles de résistance des matériaux
1 Généralités
2 Ansys
9 Analyse modale
10 Flambement
Obtenus à partir des équations de la résistance des matériaux, ces dernières dérivées
des équation de la mmc à l’aide d’hypothèses cinétiques et cinématiques.
Les degrés de liberté sont les déplacements du support géométrique (ligne moyenne
par exemple) et les rotations de la matière transverse (fibre ou section transverse),
assemblés dans un torseur cinématique.
Les efforts sont traduits par le torseur de cohésion.

2D - Eléments plaques et coques
1D - Eléments de poutres
Modèles de résistance des matériaux Hypothèses mécaniques de plaques
Géométrie des plaques
Éléments structuraux minces suivants une direction.

La géométrie est remplacée par un support à topologie surfacique −→ maillage
surfacique topologie.
Rem : la surface moyenne est classiquement utilisée, mais il est quelque fois plus
simple d’utiliser la surface inférieure ou supérieure (pour les composites multicouches
par exemple).
Il existe deux modélisations :
Plaque : la surface de description est plane et l’épaisseur est constante.
Coque : la surface de description n’est pas plane et/ou l’épaisseur est variable.
Comportement de la fibre épaisseur - plaque épaisse de Mindlin
Fibre indéformable
Une fibre normale à la surface moyenne a une cinématique de solide rigide :
Pas de raccourcissement ;
Pas de flexion ;
Pas de gauchissement.
Son mouvement est donc la superposition d’une translation et d’une rotation.
Fibre normale - plaque mince de Kirchhoff

Dans le cas de plaque de très faible épaisseur (ep << L/10), le cisaillement transverse
est négligeable. Cela a pour incidence qu’une fibre orthogonale à la surface moyenne
avant déformation, le reste au cours de la transformation.
Cela conduit à :
Une simplification du champ de contrainte car τt = 0.
On peut identifier les rotations de la fibre épaisseur avec les gradients du champ
de déplacement transverse w(x, y ).
Cinématique des plaques
Le champ de déplacements prend la forme

~M = U
U ~M + R
~ M~0 M
0
En introduisant les rotations θx et θy autour des axes x et y , et sous les hypothèses

des petites perturbations, le champ cinématique s’écrit :
 
u0 (x, y ) + z · βx (x, y )
~ M =  v0 (x, y ) + z · βy (x, y ) 
U
w0 (x, y )
Modèles de résistance des matériaux Modélisation aux éléments finis de plaques
Eléments plaques et coques
Topologie 2D
Cinématique déplacements (U, V , W ), et rotations de flexion (θx , θy ) aux noeuds.
Statique Torseur de cohésion aux noeuds (Nx , Ny , Txy , Txz , Tyz , Mfx , Mfy ).
Propriétés physiques Epaisseurs.
Supports géométriques triangle, quadrangle.
Modèles de résistance des matériaux Modélisation aux éléments finis de plaques
Eléments plaques et coques - exemples
Génération de maillages
1 Généralités
2 Ansys
9 Analyse modale
10 Flambement
Génération de maillages
Représentation d’un maillage
Structure des données de maillage

Du point de vue de la structuration des données, un maillage est constitué d’une liste
de coordonnées de noeuds (table des coordonnées) et de la liste des éléments décrit
par une collection de noeuds (table des connectivités).
Liste des noeuds Liste des éléments

2411 2412
1 1 1 11 180 115 1 1 7 8
7.8658D+01 -3.5735D+01 0.0000D+00 213 214 240 239 369 370 396 395
2 1 1 11 181 115 1 1 7 8
7.4658D+01 -3.5735D+01 0.0000D+00 214 215 241 240 370 371 397 396
... ...
Ces tables peuvent être générées manuellement (ou semi manuellement en utilisant
des scripts) ou par construction géométrique (par exemple par extrusion).
Génération de maillages Maillage libre
Les maillages libres
Algorithmes de génération automatique

Des algorithmes génèrent de façon automatique une discrétisation de la géométrie. Le
triangle et le tétraèdre sont les seules formes qui assurent la fermeture géométrique.
Il existe trois type de stratégies :
frontale :
on remplit l’espace géométrique depuis l’extérieur vers l’intérieur, cela produit des
mailles déformées au centre.
arbre :
on génère des boites englobantes de la géométrie que l’on découpent
successivement jusqu’à avoir une taille adéquate.
Voronoï :
on remplit l’espace intérieur d’un nuage de points que l’on associes ensuite avec
des cellules de Voronoï. On génère ensuite une maille entre chaque noeuds
contigus
Un lissage du maillage doit être réalisé, par une technique de triangulation de

Delaunay ou par une relaxation élastique.
Méthodes frontales
La méthode frontale consiste à mailler le domaine, de la façon la plus régulière

possible en respectant une taille imposée par l’opérateur, depuis la frontière (ligne ou
surface) vers l’intérieur jusqu’à remplissage complet du domaine.
sens du maillage
L’algorithme adapte les tailles et fait disparaître les noeuds surnuméraires, mais n’est
pas exempt de défaut en particulier sur la fin du processus.
Méthodes frontales

Avance du front
Méthodes frontales

zone vide
Arbre de construction
Le domaine est d’abord enfermé dans une boite englobante. Cette boite est découpée
suffisamment de fois pour atteindre une taille imposée par l’opérateur et pour capter
chaque point géométrique.
Arbre de construction
Les boites sont ensuite maillées à l’aide de motifs prédéfinis. Ces motifs imposent de
plus un niveau de découpage de façon à assurer les transitions entre domaines fins et
grossiers.
Cellule de Voronoï
Le domaine est rempli d’un nuage de points disposés aléatoirement.
Cellule de Voronoï
Des cellules de Voronoï sont construites autour de chaque point afin de définir le
voisinage de chaque point.
Cellule de Voronoï
Un élément est alors construit entre les points des cellules de Voronoï formant trois
voisinages concomitants.
Génération de maillages Maillage structuré
Projection transfini
On parle de maillage réglé ou transfini. Il consiste à associer une forme de référence

facile à mailler à la forme réelle.
C’est la seule façon de réaliser des maillages hexaédrique et carré contrôlés.
Maillage régulier
Il faut faire attention à la manière d’associer les domaines réels aux domaines de
référence.
Certaines transformations induisent des déformations localisées importantes.
Sur l’exemple de l’arche, si l’on choisit les points inférieurs comme base de référence
pour définir un carré de projection, le maillage subit un déformation d’ensemble, ce qui
conduit à un maillage relativement régulier.
Maillage déformé
Si l’on choisit les points extrêmes, de façon à englober la pièce, cela induit une
déformation localisée sur un coté du maillage. Certaine mailles n’assureront pas
l’unicité de la transformation (i.e. des éléments avec des Jacobiens changeant de
signe).
Génération de maillages Partition
Découpage des pièces
Les partitions sont des découpages virtuels qui permettent de décomposer la pièce à
mailler complexe en une combinaison de pièces simples : en particulier pouvant être
maillées de façon réglée.
Les frontières sont de type ligne ou surface. Ces frontières

assurent (imposent) la continuité des maillages.
Compatibilité des maillages
Il faut assurer la compatibilité des maillages.
Maillages incompatibles
Il existe des stratégies de collage de maillages non compatibles, non actuellement

utilisées de façon classique dans les codes industriels, mathématiquement simples
mais délicates à mettre en oeuvre techniquement.
Génération de maillages Degré des éléments
Degré des éléments
Les éléments finis s’appuient sur des monômes de Pascal.
1
x y
x2 xy y2
x3 x 2y xy 2 y3
···
Théoriquement les degrés sont infinis, mais dans la pratique seuls les degrés 1 et 2
sont utilisés, rarement degré 3.
Génération de maillages Degré des éléments
Eléments Serendip
Serendipity
mot inventé en 1754 par le philosophe anglais Sir Horatio WALPOLE, qualifie la faculté
à trouver la bonne information par hasard. Le mot provient du roman "The three
princes of Serendib" à qui le hasard apportait la solution à des situations fâcheuses.
Les éléments à base carrée sont de type Serendip : on fait disparaître un noeud
interne par rapport aux éléments de Hermite, ce qui implique une certaine forme de
base polynomiale.
Les éléments de degrés supérieurs à 1 sont très sensibles aux déformations de

maillages.
Génération de maillages Qualité géométrique des éléments
Critères géométriques
Il existe des critères géométriques afin d’évaluer la qualité des éléments du maillage
réalisé. Ces critères quantifient l’écart entre la valeur théorique, basée sur un élément
géométriquement parfait, et la valeur de l’élément du maillage.
Seuls les isomorphismes respectent cela :

translation, rotation et allongements
proportionnels.
Problèmes non linéaires
1 Généralités
2 Ansys
9 Analyse modale
10 Flambement
Problèmes non linéaires Types de non linéarités
Non linéarités géométriques

Les déplacements/rotations sont grands −→ modification du domaine.
Les déformations sont grandes −→ conservation des termes du second ordre
Non linéarités matériaux

Comportements élastique non linéaire.
Comportement plastiques.
Matériaux hyper-élastiques.
Non linéarités de contact

Les déplacements/rotations sont « controllés » −→ Modification des conditions
aux limites.
Les efforts de contact dépendent du problème −→ Modification du second
membre.
Formulations et mesures suivant le type de non linéarité traitée
Type Hypothèses Formulation Mesures

Matériau hpp Material Non σ
rdc non lin. Linear Only
Grands dép. [∇U] négl. Total Lagrange σ et A
Petites déf. rdc lin. ou non Update Lagrange S et E
Grands dép. hors hpp Total Lagrange σ et log
Grandes déf. rdc lin. ou non Update Lagrange S et E
Mesures de déformations
Mesures de Contraintes : linéarisée
σ Cauchy A Almansi
S Piola-Kirchhoff E Green-Lagrange
log Hencky-Mises
Problèmes non linéaires Comportement non linéaire géométrique
Flexion d’une poutre, modélisation linéaire
L0 ·F
Contrainte de flexion : σ = y · Igz
F
Contrainte de cisaillement : τ = S
2
Contrainte équivalente de Von Mises : σVM = σ2 + 3 · τ 2
Flexion d’une poutre, modélisation non linéaire géométrique
F ·sin(α) L·F F ·sin(α) L0 ·F ·cos(α)

Contrainte de traction et flexion : σ = S
+y · Igz
' S
+y · Igz
F ·cos(α)
Contrainte de cisaillement : τ = S
Non linéarité géométrique
Particularité
Les déplacements et les rotations sont grands. −→ Modification du domaine.
Les déformations sont grandes. −→ Conservation des termes du second ordre.
L’équation d’équilibre du milieu continu

~ ([σ]) + f~v = ~0
div
et affaiblie par une formulation variationnelle

Z Z Z
− tr ∇U ~ ? · σ dV + ~ ? · ~fV dV +
U U~ ? · ~fS dS = 0
Ω Ω ∂ΩF
~ = U~d et U
~ ? = ~0 sur ∂ΩU avec ∇U
~? = ∂U~?
On impose U ∂x t+∆t
Cinématiques des solides déformables
On a
∂x t+∆t ∂x t+∆t ∂x t
=
∂x 0 ∂x t ∂x 0
donc h i h ih i
F0t+∆t = Ftt+∆t F0t
Changement de référentiel
La relation de changement de référentiel
∂u ? ∂u ? ∂x t h i−1
t+∆t
= t t+∆t
= [F ? ] · Ftt+∆t
∂x ∂x ∂x
permet d’exprimer l’équilibre sur la configuration actuelle à partir d’un état précédent
Z Z
Wint = − tr ∇U ~ ? · σ dV = − tr E ? · Stt+∆t dV
Ωt+∆t Ωt
On doit pour cela introduire la contrainte de Piola-Kirchhoff de seconde espèce Stt+∆t

et le champ des déformations virtuelles
h i
1 iT
T
h
E ? = · tr Ftt+∆t · [F ? ] + [F ? ] · Ftt+∆t
2
Forces extérieures
Force conservative
Indépendante de l’évolution du domaine.
Vecteur force constant.
Force suiveuse
Modifié par l’évolution du domaine.
Vecteur force variable.
Mesures des contraintes
La mesure de contrainte doit être adaptées à l’étude traitée :

Cauchy σ t+∆t
τtt+∆t = σ t+∆t · det Ftt+∆t

Kirchhoff
−1
Ttt+∆t = Ftt+∆t · τ t+∆t

Piola-Kirchhoff I
t+∆t t+∆t
t+∆tt −T
Piola-Kirchhoff II St =T · Ft
t T
σR = Rtt+∆t · σ t+∆t · Rtt+∆t

Co-rotationnelle
Mesures des déformations
On cherche des lois de comportement objectives permettant de décrire le matériau :

σ t+∆t = H t+∆t
ou
S t+∆t = H Ett+∆t
Ces lois nécessitent l’introduction de mesures de déformations adaptées :

Déformation ingénieur
t+∆t
Green-Lagrange h
1 iT h i
Ett+∆t = tr Ftt+∆t · Ftt+∆t
− [I]
2
Formulation éléments finis
On écrit l’équilibre sur une configuration Ωt connue

Z
Wint = − tr E ? · Stt+∆t dV = − hU ? i · {R}
Ωt
On doit alors exprimer le tenseur des déformations de Green-Lagrange

n o h n oi n o
Ett+∆t = [BL ] + BNL Utt+∆t · Utt+∆t
et le tenseur des déformations virtuelles

n o h n oi
{E ? } = δEtt+∆t = [BL? ] + BNL
?
Utt+∆t · {U ? }
Ce qui donne
Z n o
T
{R} = ([BL? ] + [BNL
?
]) · [H] · ([BL ] + [BNL ]) · Utt+∆t dV
Ωt
Problèmes non linéaires Résolution de systèmes non linéaires
Résolution non linéaire
On construit un processus itératif U i dont la limite {U ∞ } correspond au champ de

t+∆t
déplacement Ut .
Un développement au premier ordre nous permet d’écrire :
∂R U i
i+1
R U = R U + ∆U ' R U i +
i
· ∆U
∂U
On réalise les opérations suivantes
Calcul du résidu {R} = R U i .

1
2 Test : si kRk < rr , fin du processus.

∂R (U i )
3 Calcul de la matrice tangente [KT ] = ∂U
.
4 Résolution du système linéaire [KT ] · {∆U} = − {R}.
Mise à jour U i+1 = U i + {∆U}.

5
6 Retour en 1.
Formulation Lagrangienne Totale
Formulation Lagrangienne Actualisée
1 Généralités
2 Ansys
9 Analyse modale
10 Flambement
Modèles de résistance des matériaux Elements de poutre
Eléments de poutre
Topologie 1D
Cinématique déplacements (U, V , W ), , rotations de flexion (θy , θz ) et rotations de
torsion (θx ) aux noeuds.
Statique Torseur de cohésion aux noeuds (Nx , Ty , Tz , Mtx , Mfy , Mfz ).
Propriétés physiques Sections (aires et moments quadratiques).
Supports géométriques linéique.
Modèles de résistance des matériaux Elements de poutre
Eléments de poutre - exemple
Modèles de résistance des matériaux Modèles rhéologiques
Eléments utilisés classiquement et régulièrement.

Dont les spécifications sont rhéologiques (comportementales) avec des supports
géométriques n’étant pas forcément définis.
Eléments ressorts
Eléments rigides
Eléments de masses
Eléments ressort
Topologie 0D (ponctuelle)
Cinématique déplacements (U, V , W ), et rotations (θx , θy , θz ) aux noeuds.
Statique Torseur aux noeuds (Fx , Fy , Fz , Mx , My , Mz ).
Propriétés physiques Raideurs axiales et raideurs de rotations.
Matériau -
Supports géométriques -
Topologie 1D (linéaire)
Cinématique Déplacements (u) aux noeuds.
Statique Force axiale (Fx ) aux noeuds.
Propriétés physiques Raideur axiale.
Matériau -
Supports géométriques linéique à deux noeuds.
Eléments rigides
Topologie -
Cinématique Déplacements (U, V , W ), et rotations (θx , θy , θz ) au noeud support
(maître).
Statique Torseur statique au noeud support (Fx , Fy , Fz , Mx , My , Mz ).
Matériau -
Supports géométriques un noeud support (maître) et autant d’autre noeuds que
nécessaire (esclaves).
Eléments de masses
Topologie 0D
Cinématique déplacements (U, V , W ), et rotations (θx , θy , θz ) au noeud.
Statique -
Propriétés physiques Masse et moments d’inerties.
Matériau -
Supports géométriques -
Modèles rhéologiques - exemples

Modélisation du telescope TAROT (CNRS - La Silla - Chilli)
Analyse modale
1 Généralités
2 Ansys
9 Analyse modale
10 Flambement
Analyse modale Equilibre dynamique
Equilibre dynamique
On se place dans le cadre de l’analyse linéaire.
hypothèse des petites perturbations (HPP)

petits déplacements et petites déformations
élasticité linéaire
forces conservatrices
On a alors :
Wint + Wext = Wacc avec Wacc = δU T · M · Ü

∂ 2 U(x,t)
La matrice de masse M vient de l’interpolation de Ü (x, t) = ∂t 2
sur le même
espace que U (x, t).
M est une matrice symétrique, définie et positive.
Analyse modale Equilibre dynamique
Matrice de masse de l’élément de barre
Si l’on prend l’exemple de l’élément linéaire de traction compression à deux noeuds, la

partie dynamique est donnée par l’expression
Z x=L
∂ 2 u (x, t)
Wacc = δu (x) · ρS · dx
x=0 ∂t 2
avec x x
ü (x, t) = 1 − · ü1 (t) + · ü2 (t)
L L
aussi

ü1
Wacc = δu1 δu2 [M]
ü2
et
ρSL 2 1
[M] =
6 1 2
Analyse modale Modes propres
Modes propres
On cherche la solution homogène de cette équation différentielle ordinaire.
K · U (t) + M · Ü (t) = 0
La solution générale est recherchée sous forme harmonique : U (t) = V · ejωt

V est le mode de vibration, la forme géométrique ;
ω
f = 2π
est la fréquence de vibration.
On cherche à résoudre :

K − ω2 · M · V = 0 avec V 6= 0
On va donc obtenir (ωi , Vi ), les fréquences (pulsations) et modes propres. On a N

modes propres (éventuellement multiples), avec N la taille des matrices K et M.
On utilise fréquemment une procédure d’itérations inverses.
méthode des puissances (depuis les hautes fréquences) ;
méthode inverse (depuis les basses fréquences).
Les modes sont normés, donc sans amplitude. Les sollicitations extérieures induisent
l’amplitude des vibrations.
Analyse modale Exemples
Analyse modale des bâtiments

L’analyse modale des bâtiments peut avoir pour motivation
l’évaluation de la réponse du bâtiment à une sollicitation sismique, vérification en
dimensionnement critique (ELU)
l’évaluation des vibrations du bâtiment soumis à une charge dynamique
permanente (vent ou matériel mécanique), vérification en fonctionnement (ELS)
Analyse modale d’un vilebrequin

Obtenir les modes propres d’un vilebrequin peut avoir pour motivation
evaluation des modes de vibration
dimensionnement d’usage
Analyse modale d’une turbine

Chercher les modes propres d’une turbine permet :
l’évaluation des modes de vibration en fonctionnement ;
le dimensionnement d’usage ;
les charges sur la plateforme.
Flambement
1 Généralités
2 Ansys
9 Analyse modale
10 Flambement
Flambement Modèles analytiques
Bifurcation mécanique
Le flambement est un mode de déformation critique qui s’applique à des structures

élancées soumises à de la compression. Il s’agit du passage brusque d’un état de
compression à un état de flexion.
Cette déformation est due en majeure partie au phénomène de rotation du tenseur des
contraintes lors de la déformation.
Le flambement est un mode de défaillance qui correspond à un état singulier de

l’équilibre mécanique − → plusieurs possibilités de réponses peuvent être suivies par la
structure : bifurcation.
Flambement Rotule - Rotule
Mfz = −y (x) · F
E · I · y 00 (x) = Mfz
F
y 00 + ω 2 · y = 0 avec ω2 =
E ·I
y (x) = A · cos (ω · x) + B · sin (ω · x)
y (0) = 0 implique A=0
y (L) = 0 implique B · sin (ω · L) = 0
ω·L=k ·π
π2 · E · I
F = k2 ·
L2
Flambement Encastrée - Libre
Mfz = (y (L) − y (x)) · F

y 00 + ω 2 · y = ω 2 · y (L)
y (x) = A · cos (ω · x) + B · sin (ω · x) + ω 2 · y (L)
y (0) = 0 implique A = −ω 2 · y (L)
y 0 (0) = 0 implique B=0
F · y (L)
y (x) = (1 − cos (ω · x))
E ·I

L 2 E ·I
cos · F =1−
E 2 · I2 F
Flambement Modélisation EF du flambement linéaire
Matrice de raideur géométrique
La matrice des raideurs tangentes d’une structure traduit une linéarisation du

comportement autour d’un point d’équilibre.
Cette matrice peut se décomposer sous la forme d’une somme de trois composantes :
Matrice de raideur tangente
KNL (U) = KL + KG (U) + Kσ (U)
KL est la partie linéaire du comportement mécanique ;

KG est la partie associée aux grandes déformations ;
Kσ traduit les grands déplacements.
C’est la modification de la matrice Kσ lors de la déformation qui produit le bifurcation.
Flambement linéaire
On s’intéresse généralement au flambement linéaire : transition depuis l’état linéaire à

l’état non linéaire
En effet, après une déformation par flambement, il y a presque systématiquement

défaillance de la structure.
Cette hypothèse permet ainsi d’affirmer qu’avant la bifurcation, le comportement est
linéaire et vérifie
K .U0 = F0
On cherche la transition de l’état linéaire à l’état non linéaire en augmentant l’effort de
référence :
Fext = λ · F0
λ est le facteur de charge.
Facteur de charge critique
Au moment du flambement, le système doit vérifier l’expression non linéaire de

l’équilibre :
[KNL (U)] · U = Fext
Mais comme on est à la limite de l’état linéaire, on vérifie aussi
U = λ · U0
Les grandes déformations sont négligées, aussi obtient-on le système :
[KL ] · λ · U0 + [Kσ (λ · U0 )] · λ · U0 = λ · F0
ce qui s’écrit encore

[KL + λ · Kσ (U0 )] · U0 = F0
La notion de bifurcation se traduit mathématiquement par une matrice singulière. On

cherche donc les valeurs critiques λc pour lesquelles on a
det ([KL ] + λc · [Kσ (U0 )]) = 0
Flambement Exemples de flambement
Exemples de flambement
Flambement global de poutres encastrées - encastrées
Flambement local du voile de la poutre
Flambement global de l’ame de la poutre
Flambement global par plissement
Flambement global par déversement
Flambement local par cloquage (effondrement) d’une poutre creuse
Résolution par méthodes de puissance
1 Généralités
2 Ansys
9 Analyse modale
10 Flambement
Position du problème
Trouver les valeurs propres de la matrice carrée [A] de dimension N.
Trouver les N couples (λi , Vi ) tels que
[A] · Vi = λi · Vi
Les Vi forment une base vectorielle, que l’on peut normer et classer par λi strictement
croissants (les modes multiples sont traités différement).
On peut donc écrire un vecteur X quelconque sous la forme
X = Σi=N
i=1 (αi · Vi )
Les αi sont des scalaires qui représentent les coordonnées du vecteur X dans la base
(Vi ).
Méthodes de puissance
Si l’on prend un vecteur quelconque écrit dans la base des vecteurs propres
X = Σi=N
i=1 (αi · Vi )
L’application de la matrice [A] à ce vecteur implique le résultat suivant

i=N
[A] · X = [A] · Σi=1 (αi · Vi ) = Σi=N i=N
i=1 (αi · [A] · Vi ) = Σi=1 (αi · λi · Vi )
Si on applique m fois cette multiplication et que l’on prend en compte l’organisation

croissante des λi , on obtient le résultat suivant
m i=N
λm

 [A] · X = Σi=1 (αi ·
 i ·V
i ) m
= λN · Σi=1 αi · λλi
m i=N
· Vi
 N
' αN · λm N · VN

On obtient ainsi une évaluation de la valeur propre VN .
Méthodes de puissance (suite)
On peut ensuite reprendre le même processus avec le vecteur
Y = X − (X · VN ) · VN = X − αN · VN
L’application de la matrice [A] à ce vecteur donne
[A] · Y = Σi=N−1
i=1 (αi · λi · Vi )
Ainsi, en appliquant à nouveau m fois cette multiplication, on obtient
[A]m · Y ' αN−1 · λm

N−1 · VN−1
On reproduit N fois ce processus de façon à identifier depuis VN jusqu’à V1 .
Méthodes de puissance inverse
On préfère généralement obtenir d’abord les plus petites valeurs propres (qui
correspondent aux basses fréquences ou premiers modes critiques de flambement).
Pour cela, on inverse la stratégie précédente

h i 1
A−1 · X = Σi=N
i=1 αi · · Vi
λi
Réitérer m fois l’application mène à

h im 1
A−1 · X ' α1 · m · V1
λ1
On peut ensuite remonter de V1 jusqu’à VN par projections successives.

Ansys Calcul

Transféré par

Droits d'auteur :

Formats disponibles

Ansys Calcul

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Ansys Calcul

Transféré par

Droits d'auteur :

Formats disponibles

Computational approach of systems design in engineering

Calcul Scientifique pour la Conception de Systèmes, Partie 2 (CSCS2)

3 Modèles éléments finis en mécanique des solides

5 Modèles de résistance des matériaux

7 Problèmes non linéaires

8 Modèles de résistance des matériaux

11 Résolution par méthodes de puissance

Place de la simulation dans la chaine de conception mécanique

Processus « classique » de conception d’un produit

Coût induit par une modification selon les phases

Objectifs de la simulation en conception mécanique

Évaluation des efforts dans une structure

Vérifier la tenue d’un prototype

Optimiser un système (géométrique, matériau)

Défauts/qualités des méthodes numériques

Equations d’équilibre des milieux continus déformables

La résolution d’un problème de mécanique des milieux déformables est souvent

Equations d’équilibre des milieux continus déformables

On peut quelques fois obtenir une solution à partir d’hypothèses permettant

Equations d’équilibre des milieux continus déformables

Mécanique des milieux continus

L’équation d’état est écrit en temps et en espace.

On approxime les équations différentielles d’abords en espace pour réduire l’EDP en

f (U (x, t) , t) = 0 −→ g (Ui , t) = 0 −→ h (Ui , tj ) = 0

(U (x, t)) fonction continue en temps et en espace, i ∈ [1, Ni ] discrétisation en espace,

Approximation des équations

On utilise plusieurs stratégies de résolutions numériques des équations différentielles

Résolution des EDO

Résolution des EDP

Principe de la méthode des éléments finis

Méthode des éléments finis

Formulation en déplacements d’une barre

On étudiera le cas statique d’une barre à section constante de matériau homogène en

Formulation faible de barre

Interpolation cinématique sur une barre

δu (x) = N1 (x) · δu1 + N2 (x) · δu2

Approximation des champs de déformation par éléments

du (x) dN1 (x) dN2 (x) h i u 

Matrice de raideur d’une barre

δWint = − {δUe }T [Ke ] {Ue }

Matrice de raideur de l’élément de barre linéaire à deux noeuds

Boucle de modélisation - Analyse de l’Objet de l’Etude

Etape 1 : Choix de modélisation analytique

Boucle de modélisation - Analyse de l’Objet de l’Etude (suite)

Choix de la pièce (ensemble de pièces) à étudier.

Boucle de modélisation - Modélisation par Eléments Finis

Etape 2 : Choix et réalisation du modèle éléments finis

Boucle de modélisation - Modélisation par Eléments Finis (suite)

Réaliser un maillage de la géométrie

Boucle de modélisation - Résolution Numérique

Etape 3 : Résolution numérique

Boucle de modélisation - Analyse des Résultats

Etape 4 : Analyse de la validité et des résultats de la modélisation

Boucle de modélisation - Analyse des Résultats

Validité de la modélisation (suite)

3 Modèles éléments finis en mécanique des solides

5 Modèles de résistance des matériaux

7 Problèmes non linéaires

8 Modèles de résistance des matériaux

du (x) dN1 (x) dN2 (x) h i u