Cours 8ITEC

Simulation Mécanique

MEF, maillage, conditions limites, post-traitement

Objectifs pédagogiques

Comprendre le principe de la méthode des éléments finis (MEF)
Maîtriser la création d'un maillage adapté au problème
Définir correctement les conditions aux limites et les chargements
Interpréter les résultats de simulation et valider un dimensionnement

Introduction : Tester virtuellement avant de fabriquer

Exemple concret : Crash-test virtuel automobile

Avant de détruire une voiture réelle en crash-test (coût : ~300 000€), les constructeurs réalisent des milliers de simulations numériques. La simulation MEF prédit les déformations, les zones de rupture et les efforts sur les occupants avec une précision de 95%.

Pourquoi la simulation mécanique ?

Réduction des coûts : moins de prototypes physiques
Optimisation : tester des variantes rapidement
Sécurité : valider sans risque de casse
Compréhension : visualiser ce qui est invisible (contraintes internes)

Définition : La simulation mécanique utilise des modèles numériques pour prédire le comportement d'une structure soumise à des sollicitations (forces, pressions, températures).

1La Méthode des Éléments Finis (MEF)

Principe fondamental

La MEF divise une géométrie complexe en un grand nombre de petits éléments simples (triangles, tétraèdres...). Les équations de la mécanique sont résolues sur chaque élément, puis assemblées pour obtenir le comportement global.

Les étapes d'une simulation MEF

1. Pré-traitement

Géométrie CAO → Maillage → Matériaux → Conditions limites

2. Résolution

Calcul matriciel : [K] × {u} = {F}

3. Post-traitement

Visualisation des résultats → Analyse → Validation

Logiciels de simulation courants

Logiciel	Éditeur	Spécialité
ANSYS	Ansys Inc.	Multi-physique (référence)
SolidWorks Simulation	Dassault	Intégré à la CAO
Abaqus	Dassault	Non-linéaire, crash
Fusion 360 Simulation	Autodesk	Cloud, accessible

Equation fondamentale : [K] × {u} = {F}
ou [K] = matrice de rigidite, {u} = deplacements, {F} = forces appliquees

Formule de Von Mises (critere de rupture)

σ_VM = √[((σ₁-σ₂)² + (σ₂-σ₃)² + (σ₃-σ₁)²) / 2]

σ₁, σ₂, σ₃ = contraintes principales

Exemple de verification :

Piece en acier S355 (Re = 355 MPa), simulation donne σ_{VM max} = 180 MPa

Cs = Re / σ_VM = 355 / 180 = 1,97 → OK (Cs > 1,5)

2Le Maillage : Clé de la Précision

Qu'est-ce que le maillage ?

Le maillage (mesh) est la discrétisation de la géométrie en éléments finis. Plus le maillage est fin, plus le résultat est précis, mais plus le calcul est long.

Types d'éléments

Éléments 2D

Triangle (Tri3, Tri6) : flexible, s'adapte à toutes géométries
Quadrangle (Quad4, Quad8) : plus précis, géométries régulières
Usage : coques, plaques

Éléments 3D

Tétraèdre (Tet4, Tet10) : maillage automatique
Hexaèdre (Hex8, Hex20) : plus précis, difficile à générer
Usage : pièces massives

Critères de qualité du maillage

Critère	Définition	Valeur idéale
Aspect ratio	Rapport longueur/largeur	< 3:1
Skewness	Déformation angulaire	< 0,5
Jacobien	Distorsion de l'élément	> 0,3

Raffinement local

On raffine le maillage dans les zones critiques (congés, concentrations de contraintes) tout en gardant un maillage grossier ailleurs pour optimiser le temps de calcul.

3Conditions aux Limites et Chargements

Importance des conditions aux limites

Les conditions aux limites définissent comment la pièce est maintenue (encastrement, appui...) et quelles charges elle subit. Une erreur ici invalide totalement les résultats.

Types de liaisons (blocages)

Liaison	DDL bloqués	Exemple réel
Encastrement	6 DDL (Tx, Ty, Tz, Rx, Ry, Rz)	Poutre soudée au mur
Appui simple	1 à 3 translations	Poutre posée sur support
Pivot	5 DDL (rotation libre)	Axe de roue
Glissière	5 DDL (translation libre)	Tiroir

Types de chargements

Forces et moments

Force ponctuelle : F en un point (N)
Force répartie : sur une ligne (N/m)
Moment : couple de rotation (N.m)

Pressions et autres

Pression : sur une surface (Pa = N/m²)
Gravité : poids propre
Température : dilatation thermique

Règle de Saint-Venant : Les conditions aux limites n'affectent les contraintes que localement. À distance suffisante, le champ de contraintes est uniforme.

4Post-traitement et Interprétation

Grandeurs calculées

La simulation fournit les déplacements aux noeuds, d'où sont dérivées les déformations puis les contraintes via la loi de Hooke.

Types de contraintes

Contrainte	Symbole	Usage
Von Mises	σ_VM	Critère de rupture (matériaux ductiles)
Principales (σ1, σ2, σ3)	σ_I, σ_II, σ_III	Directions principales
Cisaillement max	τ_max	Critère de Tresca

Critère de validation

Coefficient de sécurité : Cs = σ_limite / σ_{max calculée}

Cs > 1,5 : conception sécuritaire (cas général)
Cs > 2 à 3 : structures critiques (levage, aéronautique)
Cs < 1 : RUPTURE - redimensionner !

Convergence du maillage

Pour valider un résultat, on vérifie la convergence : raffiner le maillage jusqu'à ce que le résultat ne change plus significativement (<5%).

Méthode : Maillage 1 → σ = 120 MPa → Maillage 2 (2× plus fin) → σ = 125 MPa → Écart 4% → OK

Attention aux singularités : Les contraintes tendent vers l'infini aux angles vifs. Il faut ajouter des congés ou ignorer ces zones localisées.

Résumé en 5 points clés

1La MEF discrétise une géométrie en éléments simples pour résoudre [K]{u} = {F}.
2Le maillage doit être raffiné dans les zones critiques (aspect ratio <3, skewness <0,5).
3Les conditions aux limites (liaisons + charges) sont essentielles à la validité du modèle.
4La contrainte de Von Mises est le critère principal pour valider un dimensionnement.
5La convergence doit être vérifiée : résultat stable malgré le raffinement du maillage.

Mini-Quiz

Question 1 : Que représente la matrice [K] dans l'équation de la MEF ?

a) La matrice des forces appliquées

b) La matrice de rigidité de la structure

c) La matrice des déplacements

Réponse : b) [K] relie les déplacements aux forces via les propriétés du matériau et la géométrie

Question 2 : Pourquoi raffiner le maillage dans les zones de congés ?

a) Pour accélérer le calcul

b) Pour capturer les concentrations de contraintes

c) Pour réduire la taille du fichier

Réponse : b) Les gradients de contraintes sont élevés dans ces zones, un maillage fin est nécessaire

Question 3 : Un coefficient de securite Cs = 0,8 signifie que :

a) La piece est surdimensionnee

b) La piece est correctement dimensionnee

c) La piece va se rompre

Reponse : c) Cs < 1 signifie σ_max > σ_limite → rupture certaine

Question 4 : Simulation sur alu 6061-T6 (Re=276 MPa) donne σ_VM=120 MPa. Quel est le Cs ?

a) 0,43

b) 2,3

c) 156

Reponse : b) Cs = 276/120 = 2,3 → piece bien dimensionnee

Question 5 : Qu'indique un aspect ratio de 10:1 sur un element de maillage ?

a) Element de bonne qualite

b) Element trop allonge (mauvaise qualite)

c) Element trop petit

Reponse : b) Aspect ratio > 3:1 indique un element deforme qui degrade la precision

Cours précédent Cours suivant