Méthode des éléments finis simplifiés | Sciences de l'ingénieur - Structures et contraintes
Introduction à la méthode des éléments finis
Découvrez la méthode des éléments finis pour l'analyse des structures
Définition de la méthode des éléments finis
Qu'est-ce que la méthode des éléments finis ?
La méthode des éléments finis (MEF) est une technique numérique utilisée pour résoudre des problèmes d'ingénierie et de physique. Elle consiste à diviser un domaine continu en sous-domaines (éléments) reliés entre eux par des points (nœuds), permettant de résoudre des équations différentielles partielles complexes.
- 1 Discrétisation du domaine en éléments simples
- 2 Approximation des inconnues dans chaque élément
- 3 Assemblage des contributions élémentaires
- 4 Résolution du système d'équations
Types d'éléments
Classification des éléments
Les éléments 1D sont utilisés pour modéliser des structures linéiques comme les barres, les poutres et les câbles. Un élément barre est caractérisé par deux nœuds et ne supporte que des efforts de traction ou compression.
Les éléments 2D sont utilisés pour modéliser des structures planes. Les triangles sont les plus simples avec 3 nœuds, les quadrangles avec 4 nœuds offrent une meilleure précision. Les éléments peuvent être de type membrane (sollicitations dans le plan) ou plaque (flexion hors-plan).
Les éléments 3D modélisent des volumes complets. Les tétraèdres (4 nœuds) sont les plus simples, les hexaèdres (8 nœuds) offrent plus de précision. Ils sont utilisés pour des analyses complexes de contraintes dans des volumes tridimensionnels.
Maillage
Discrétisation du domaine
Le maillage consiste à subdiviser un domaine continu en éléments discrets. Le maillage doit être suffisamment fin pour capturer les variations locales des contraintes et déplacements, mais pas trop pour éviter une surcharge computationnelle.
- Forme des éléments (éviter les éléments trop allongés)
- Graduation (transition progressive entre zones fines et grossières)
- Continuité (compatibilité des déplacements entre éléments)
- Densité (zones critiques avec maillage plus fin)
Matrices de rigidité
Relation force-déplacement
La matrice de rigidité relie les forces nodales aux déplacements nodaux pour un élément donné. Pour un élément barre, la matrice de rigidité est une matrice 2×2 qui dépend de la rigidité axiale EA/L :
Les matrices de rigidité élémentaires sont assemblées pour former la matrice de rigidité globale du système. Chaque contribution élémentaire est ajoutée à la position correspondante dans la matrice globale selon la numérotation des nœuds.
Résolution du système
Équation d'équilibre
L'équation fondamentale de la méthode des éléments finis est :
Où K est la matrice de rigidité globale, U le vecteur des déplacements inconnus, et F le vecteur des forces nodales.
Avant de résoudre le système, il faut appliquer les conditions aux limites : déplacements imposés (appuis fixes) et forces appliquées. Les déplacements connus sont fixés, les forces connues sont placées dans le vecteur F.
Post-traitement
Analyse des résultats
Après résolution du système, on peut calculer les contraintes dans chaque élément à partir des déplacements nodaux. Pour un élément barre, la contrainte est donnée par :
- Déformée de la structure (déplacements)
- Répartition des contraintes
- Cartographie des déformations
- Évolution temporelle (dynamique)
Exercice d'application
Système à deux barres
Considérons un système composé de deux barres en série. La barre 1 a une longueur L₁ = 1m, une section S₁ = 10 cm² et un module d'élasticité E₁ = 200 GPa. La barre 2 a L₂ = 2m, S₂ = 5 cm² et E₂ = 70 GPa. Une force F = 100 kN est appliquée au nœud 3.
1. Calculer les matrices de rigidité élémentaires.
2. Former la matrice de rigidité globale.
3. Résoudre le système pour obtenir les déplacements.
4. Calculer les contraintes dans chaque barre.
Solution de l'exercice
Correction détaillée
Pour la barre 1 : K₁ = (E₁S₁)/L₁ = (200×10⁹ × 10×10⁻⁴)/1 = 2×10⁷ N/m
Pour la barre 2 : K₂ = (E₂S₂)/L₂ = (70×10⁹ × 5×10⁻⁴)/2 = 1.75×10⁷ N/m
Les matrices de rigidité sont :
Après assemblage et application des conditions aux limites (nœud 1 fixe), la matrice globale devient :
Soit : K = [[3.75×10⁷, -1.75×10⁷], [-1.75×10⁷, 1.75×10⁷]]
Le système à résoudre est : [[3.75×10⁷, -1.75×10⁷], [-1.75×10⁷, 1.75×10⁷]] × [u₂, u₃] = [0, 100×10³]
Résolution : u₂ = 4×10⁻³ m, u₃ = 10×10⁻³ m
Contrainte dans barre 1 : σ₁ = E₁ × (u₂ - 0)/L₁ = 200×10⁹ × 4×10⁻³ = 800 MPa
Contrainte dans barre 2 : σ₂ = E₂ × (u₃ - u₂)/L₂ = 70×10⁹ × 6×10⁻³/2 = 210 MPa
Résumé
Points clés
- Discrétisation du domaine en éléments
- Définition des matrices de rigidité élémentaires
- Assemblage de la matrice de rigidité globale
- Application des conditions aux limites
- Résolution du système d'équations
- Calcul des contraintes et déformations
- Adaptabilité à des géométries complexes
- Flexibilité pour différents types de matériaux
- Précision contrôlable par le raffinement du maillage
- Capacité à traiter des problèmes multi-physiques
Conclusion
Félicitations !
Continuez à pratiquer pour renforcer vos compétences en ingénierie