Travaux Dirigés 1: Elément finis 1D
Exercice 1: Eléments finis 1D
d2u
Soit l'équation différentielle suivante : +a=0, x 0,100 (1)
dx 2
1- Chercher la solution analytique de cette équation avec a=1 et les conditions aux limites
du
sont : u(0) = 100, (100)=0
dx
2- Déterminer la forme faible de l'équation (1).
3- Résoudre numériquement l'équation (1), en utilisant une approximation linéaire à une
dimension par éléments finis, dans les cas suivants:
3.1- Le domaine ne contenant qu'un seul élément.
3.2- Le domaine contenant deux éléments.
3.3- Le domaine contenant 3 éléments.
4- Calculer le vecteur élémentaire associé dans le cas a=x .
Exercice 2: Barre en traction
On considère la déformée d'une colonne de section sous son propre
d2u
poids. EA 2 +.g.A=0, x 0, L (2)
dx
Où E le module de Young du matériau, ρ sa masse volumique, g la
constante gravitationnelle et u(x) la variable d'allongement local.
1- Déterminer les conditions aux limites du problème.
2- Déterminer la forme faible de l'équation (2).
3- Résoudre numériquement l'équation (2), en utilisant une approximation linéaire à une
dimension par éléments finis, dans les cas suivants:
3.1- Le domaine ne contenant qu'un seul élément.
3.2- Le domaine contenant deux éléments. .
4- Calculer le vecteur des réactions externes.
Exercice 3: Assemblage
1/2
Dans cet exercice on propose d'étudier la technique d'assemblage de matrices et vecteurs
élémentaires pour un maillage numéroté de deux façons différentes. Les numéros des
éléments sont donnés par des numéros encadrés.
Pour chaque élément indicé par "e", on considère les formes élémentaires générales
suivantes :
k11
e
k12e
f1e
k = e
e
, f = e
e
k12 k e22 f 2
1- Déterminer la "connectivité" du maillage pour les 2 choix.
2- Dans les 2 choix, donner les formes globales de la matrice [K] et du vecteur{F}
obtenues par assemblage des formes élémentaires.
Exercice 4: Thermique 1D
Soit le modèle mathématique suivante :
d 2T(x) dT(x)
k +a -bT(x)+f(x)=0, avec f(x)=x
dx 2 dx
Les conditions aux limites de ce problème sont :
dT
T(x = 0) = T0 , k (x=L)=F0
dx
1- Déterminer la forme faible du problème.
2- Après la discrétisation du domaine par 4 éléments de longueurs égales, calculer les
matrices et vecteurs élémentaires résultant de la discrétisation de la forme faible sur
un élément quelconque.
3- Procédez à l'assemblage du système avec la prise en compte de la condition de
Dirichlet.
2/2