Corrigé détaillé (partiel) du TD : Modélisation et simulation mécanique
Fichier source : [Link]
Remarque : Ce document contient un corrigé détaillé pour la PARTIE I (éléments barre) et
quelques exercices de la PARTIE II. Si vous souhaitez le corrigé complet de tous les exercices,
dites-le; je continuerai et ajouterai les autres exercices dans un fichier mis à jour.
Partie I — Structures de type barre
Exercice 1 (élément barre 1D) - Principes et matrice de raideur
1) Approximation du champ de déplacement
Pour un élément barre unidimensionnel de longueur L, avec deux nœuds 1 et 2 et déplacements
scalaires u1 et u2, on adopte l'interpolation linéaire :
u(x) = N1(x) u1 + N2(x) u2
avec N1(x) = 1 - x/L, N2(x) = x/L (coordonnée x mesurée depuis le nœud 1).
2) Coordonnées généralisées en fonction des degrés de liberté
Les coordonnées généralisées α = [u1 u2]^T.
3) Champ de déformation
La déformation ε(x) = du/dx = dN1/dx u1 + dN2/dx u2 = (-1/L) u1 + (1/L) u2 = (1/L)(u2 - u1).
4) Matrice de rigidité élémentaire
L'énergie de déformation de l'élément : U = 1/2 ∫_0^L E A ε^2 dx = 1/2 (EA/L) (u2 - u1)^2.
D'où la matrice de raideur élémentaire (en coordonnées nodales u1,u2):
k_e = (EA/L) [[1, -1], [-1, 1]].
Application : arbre en trois tronçons
— Tronçons de sections A1, A2=A/2, A3=A1 et même longueur L.
On écrit k1 = EA/L, k2 = E(A/2)/L = (EA)/(2L), k3 = EA/L.
Assembler la matrice globale et appliquer conditions limites et charges. Les déplacements
inconnus se déterminent en résolvant K_global * u = F_reduite. Les contraintes dans chaque
tronçon: σ = E * ε = E*(Δu / L) avec Δu différence de déplacements nodaux de l'élément.
Exercice 2 (barre 2D) — Treillis de 3 membres
Données numériques (pour applications) : E = 200 GPa, L = 2 m, Section tubulaire carrée
extérieure 100 mm, ép=5 mm, P = 250 kN.
Section A et aire efficace sont calculées :
A_tube = a_ext^2 - a_int^2 où a_int = a_ext - 2*ép = 100 - 10 = 90 mm => A = 100^2 - 90^2 =
10000 - 8100 = 1900 mm^2 = 1.9e-3 m^2.
1) Matrice de raideur d'un élément barre 2D
Pour une barre faisant un angle θ avec l'horizontale, la matrice élémentaire (4x4) en coordonnées
nodales [u1 v1 u2 v2] est :
k_e = (EA/L) * [ c^2 c s -c^2 -c s
c s s^2 -c s -s^2
-c^2 -c s c^2 c s
-c s -s^2 c s s^2 ]
où c = cosθ, s = sinθ.
2) Assemblage
Identifier les nœuds 1,2,3 et assembler chaque k_e dans K_global 6x6 (2 dof par nœud).
Appliquer conditions d'appui (ex: encastré / fixe) pour réduire le système et résoudre K_red *
u_red = F_red.
3) Calcul numérique (exemple simplifié)
Pour le maillage élément 1: 1-2, 2-3, 1-3, déterminer les angles : si la géométrie est un triangle
équilatéral de côté L, alors θ_12 = 0°, θ_23 = 60°, θ_13 = 60° (exemple). Calculer k_e pour
chaque élément et assembler.
Ensuite résoudre pour V1..V4 (notation du TD) et retrouver contraintes: σ = (E/L) * (u2_local -
u1_local).
Exercice 3 (système de manutention, 5 barres) — méthode
1) Justification u1 = u3 = 0 : symétrie et liaisons (conditions données).
2) Partitions : ne garder que degrés de liberté non nuls (ex: v1,v3, etc.).
3) Établissement de K_global : assembler en tenant compte des rotations et transformations
locales->globales.
4) Résolution : K_red * u_red = F_red -> calcul de v1 et v3.
5) Vérification résistance de la barre 1-3 : calcul de la contrainte σ = (E/L) * Δl/A et comparaison
avec σ_admissible = 240 MPa.
Exercice 4 (treillis 5 éléments) — calculs détaillés (méthode)
Procédure :
1. Pour chaque élément (AB, BC, AD, DC, BD), écrire k_e en ne gardant que les dof non nuls.
2. Assembler K_global (taille réduite en tenant compte des appuis et symétries).
3. Appliquer charges (P donné) et résoudre pour déplacements.
4. Calculer déformations et contraintes dans l'élément BD via projection sur l'axe de la barre.
Exemple numérique fourni dans le TD : E=210 GPa, S=100 mm^2, L=1 m, P=100 kN. On
calcule k=EA/L = 210e3*100e-6/1 = 21 kN/mm? (attention unités). Veiller à l'homogénéité des
unités (m->mm).
Partie II — Structures de type poutre (extraits)
Exercice 5 (élément poutre en flexion)
Interpolation de la flèche v(x) pour un élément 2-nœuds (Hermite cubique) :
v(x) = N1(x) v1 + N2(x) θ1 + N3(x) v2 + N4(x) θ2
où N1..N4 sont les fonctions de forme Hermite (polynômes cubiques) définies sur [0,L].
Matrice de rigidité élémentaire en flexion (Euler-Bernoulli) :
k_e = (EI/L^3) * [[12, 6L, -12, 6L], [6L, 4L^2, -6L, 2L^2], [-12, -6L, 12, -6L], [6L, 2L^2, -6L,
4L^2]].
Énergie de déformation : U = 1/2 {v}^T k_e {v}.
Applications numériques (exemple encastrement - charge concentrée):
Résoudre le système assemblé et extraire déplacements et rotations nodales. Puis calculer
moments et contraintes flectrices σ = M*y/I.
---
Fin du corrigé partiel fourni ici.
Si vous voulez :
- le corrigé complet (tous les exercices) au format Word, je peux le produire mais il est
volumineux — précisez si vous souhaitez tout en un seul fichier ou plusieurs fichiers par partie.
- des étapes numériques détaillées (calculs manuels chiffrés) pour un exercice précis, indiquez
lequel et je l'ajoute.
