Pour un élément linéaire 1D, généralement utilisé dans la méthode des éléments finis pour des
problèmes unidimensionnels (par exemple, pour des barres ou des poutres), les fonctions de
forme sont des fonctions d'interpolation qui sont linéaires en fonction de la coordonnée
spatiale.
Contexte
Dans le cas d'un élément linéaire 1D, nous avons deux noeuds : un noeud à x1x_1x1 et un
noeud à x2x_2x2. Les fonctions de forme doivent être définies de manière à ce qu'elles soient
égales à 1 au noeud correspondant et égales à 0 à l'autre noeud.
1. Forme des fonctions de forme
Les fonctions de forme N1(x)N_1(x)N1(x) et N2(x)N_2(x)N2(x) doivent être linéaires, c’est-
à-dire de la forme :
N1(x)=a1+b1xN_1(x) = a_1 + b_1 xN1(x)=a1+b1x N2(x)=a2+b2xN_2(x) = a_2 + b_2 xN2
(x)=a2+b2x
où a1,b1,a2,b2a_1, b_1, a_2, b_2a1,b1,a2,b2 sont des coefficients à déterminer. Ces fonctions
doivent satisfaire les conditions suivantes :
N1(x1)=1N_1(x_1) = 1N1(x1)=1, N1(x2)=0N_1(x_2) = 0N1(x2)=0
N2(x1)=0N_2(x_1) = 0N2(x1)=0, N2(x2)=1N_2(x_2) = 1N2(x2)=1
Cela signifie que N1(x)N_1(x)N1(x) vaut 1 au noeud 1 (situé en x1x_1x1) et 0 au noeud 2
(situé en x2x_2x2), tandis que N2(x)N_2(x)N2(x) vaut 1 au noeud 2 et 0 au noeud 1.
2. Détermination des coefficients des fonctions de forme
Nous devons déterminer les coefficients a1,b1,a2,b2a_1, b_1, a_2, b_2a1,b1,a2,b2 pour
chaque fonction de forme en imposant les conditions aux noeuds.
Pour N1(x)N_1(x)N1(x) :
1. N1(x1)=1N_1(x_1) = 1N1(x1)=1, donc : a1+b1x1=1a_1 + b_1 x_1 = 1a1+b1x1=1
2. N1(x2)=0N_1(x_2) = 0N1(x2)=0, donc : a1+b1x2=0a_1 + b_1 x_2 = 0a1+b1x2=0
En résolvant ce système d'équations, nous obtenons les valeurs de a1a_1a1 et b1b_1b1.
Pour N2(x)N_2(x)N2(x) :
1. N2(x1)=0N_2(x_1) = 0N2(x1)=0, donc : a2+b2x1=0a_2 + b_2 x_1 = 0a2+b2x1=0
2. N2(x2)=1N_2(x_2) = 1N2(x2)=1, donc : a2+b2x2=1a_2 + b_2 x_2 = 1a2+b2x2=1
En résolvant ce système d'équations, nous obtenons les valeurs de a2a_2a2 et b2b_2b2.
3. Résolution du système pour N1(x)N_1(x)N1(x) et N2(x)N_2(x)N2(x)
Pour N1(x)N_1(x)N1(x) :
À partir des équations :
a1+b1x1=1a_1 + b_1 x_1 = 1a1+b1x1=1 a1+b1x2=0a_1 + b_1 x_2 = 0a1+b1x2=0
Soustrayons la deuxième équation de la première pour éliminer a1a_1a1 :
(b1x1−b1x2)=1(b_1 x_1 - b_1 x_2) = 1(b1x1−b1x2)=1 b1(x1−x2)=1b_1 (x_1 - x_2) = 1b1
(x1−x2)=1 b1=1x1−x2b_1 = \frac{1}{x_1 - x_2}b1=x1−x21
Ensuite, substituons cette valeur de b1b_1b1 dans l'une des équations pour déterminer
a1a_1a1 :
a1+b1x1=1a_1 + b_1 x_1 = 1a1+b1x1=1 a1+1x1−x2x1=1a_1 + \frac{1}{x_1 - x_2} x_1 =
1a1+x1−x21x1=1 a1=1−x1x1−x2a_1 = 1 - \frac{x_1}{x_1 - x_2}a1=1−x1−x2x1
a1=x2−x1x1−x2=−x1−x2x1−x2a_1 = \frac{x_2 - x_1}{x_1 - x_2} = -\frac{x_1 - x_2}{x_1 -
x_2}a1=x1−x2x2−x1=−x1−x2x1−x2
Pour N2(x)N_2(x)N2(x) :
À partir des équations :
a2+b2x1=0a_2 + b_2 x_1 = 0a2+b2x1=0 a2+b2x2=1a_2 + b_2 x_2 = 1a2+b2x2=1
Soustrayons la première équation de la deuxième pour éliminer a2a_2a2 :
b2(x2−x1)=1b_2 (x_2 - x_1) = 1b2(x2−x1)=1 b2=1x2−x1b_2 = \frac{1}{x_2 - x_1}b2=x2
−x11
Ensuite, substituons cette valeur de b2b_2b2 dans l'une des équations pour déterminer
a2a_2a2 :
a2+b2x1=0a_2 + b_2 x_1 = 0a2+b2x1=0 a2+1x2−x1x1=0a_2 + \frac{1}{x_2 - x_1} x_1 =
0a2+x2−x11x1=0 a2=−x1x2−x1a_2 = -\frac{x_1}{x_2 - x_1}a2=−x2−x1x1
4. Expression finale des fonctions de forme
En résolvant les équations, nous obtenons les fonctions de forme :
N1(x)=x2−xx2−x1N_1(x) = \frac{x_2 - x}{x_2 - x_1}N1(x)=x2−x1x2−x
N2(x)=x−x1x2−x1N_2(x) = \frac{x - x_1}{x_2 - x_1}N2(x)=x2−x1x−x1
Conclusion
Les fonctions de forme pour un élément linéaire 1D sont :
N1(x)=x2−xx2−x1N_1(x) = \frac{x_2 - x}{x_2 - x_1}N1(x)=x2−x1x2−x
N2(x)=x−x1x2−x1N_2(x) = \frac{x - x_1}{x_2 - x_1}N2(x)=x2−x1x−x1
Ces fonctions sont linéaires et satisfont les conditions d'interpolation : elles sont égales à 1 au
noeud correspondant et égales à 0 à l'autre noeud.
Les fonctions de forme N1(x) et N2(x) sont linéaires, et ont la forme générale suivante :
N1(x)=a1+b1x
N2(x)=a2+b2x
Nous devons maintenant déterminer les coefficients a1,b1,a2,b2 pour chaque fonction de
forme en imposant les conditions aux noeuds.
3. Conditions aux nœuds : Les conditions aux noeuds sont les suivantes :
Ces conditions assurent que chaque fonction de forme est égale à 1 au noeud auquel elle est
associée et égale à 0 au noeud opposé.
4. Résolution du système pour N1(x) et N2(x)
Nous allons maintenant résoudre le système d'équations pour déterminer les coefficients.
5. Expression des fonctions d'interpolation
Les fonctions d'interpolation linéaires pour un élément 1D sont donc :