Table des matières

1 Motivation 2
1.1 Qu’est ce qu’une EDP . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2 Éléments finis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3

2 Éléments Finis de degré 1 pour le problème de Dirichlet homogène en di-

mension n = 1 4

3 Éléments Finis de degré 2 pour le problème de Dirichlet homogène en di-

mension n = 1 9
Chapitre 1

Motivation

1.1 Qu’est ce qu’une EDP

Une EDP est une équation dont l’inconnue est une fonction et qui relie la fonction à ses
dérivées partielles.
L’exemple le plus simple est celui du mouvement d’un corps à un point sur la droite. La
variable x correspond alors au temps et le mouvement est d’écrit par l’équation

u00 (x) = f (u(x))

C’est le célèbre principe fondamentale de la dynamique. Une EDP est alors une relation entre
les variables et les dérivées partielles de u. On se limite pour les énoncés au cas de fonctions
d’une seule variable, mais les notions qui suivent se généralisent facilement aux fonctions de n
variables réelles.

Figure 1.1 – EDP en 1D

L’EDP elle-même ne décrit qu’une contrainte sur trois que notre fonction doit satisfaire,
elle doit également satisfaire certaines conditions aux limites et une certaine condition initiale

2
Mini Projet 1 Master MASI

1.2 Éléments finis

La méthode des éléments-finis (MEF) est une méthode d’approximation numérique de so-
lutions de problèmes aux limites statiques ou dynamiques tels que : diffusion thermique -
mécanique des milieux continus (solides et fluides) - électromagnétisme mais en fait, absolu-
ment tous les problèmes d’équations aux dérivées partielles (EDP) aux limites. II s’agit, comme
dans toutes les méthodes numériques, de trouver une approximation discrète. Pour faire bref,
d’un problème différentiel aux limites linéaire, on trouve une formulation variationnel associées
équivalente, dont on calcule une approximation de la solution en projetant sur un espace de
dimension finie, ce qui revient à résoudre au final un système linéaire. L’appellation éléments
finis vient de la décomposition du domaine d’étude en éléments : ils sont souvent représentés
par un maillage.

3 FS Agadir (2020-2021)
Chapitre 2

Éléments Finis de degré 1 pour le

problème de Dirichlet homogène en
dimension n = 1

Ce qui suit a été extrait du ”MÉMOIRE DE MASTÈRE : ELEMENTS

FINIS : DU CLASSIQUE AU ISOGEOMETRIQUE” avec de nombreuses
corrections.

Description mathématique
Considérons le problème de Dirichlet homogène en dimension n = 1 :
−αu00 (x) + βu(x) = f (x) sur ] a, b[


u(a) = u(b) = 0
où α, β sont des constantes et f est une fonction continue.
La formulation variationnelle de ce problème s’écrit : Trouver la fonction u
appartenant à H01 [a, b] telle que pour tout v ∈ H01 [a, b] on ait :
Z b Z b Z b
0 0
αu (x)v (x)dx + βu(x)v(x)dx = f (x)v(x)dx
a a a

L’existence et l’unicité d’une solution pour ce problème se démontrent en utili-

sant le théorème de Lax-Milgram.
Introduisons une discrétisation de l’intervalle [a, b] en N sous-intervalles ou
éléments Ti = [xi−1 , xi ]. Les éléments Ti n’ont pas forcément même longueur.
V0,h est alors l’espace des fonctions continues affines par morceaux (affines sur
les segments Ti ) et nulles aux extrémités a et b.
Rappelons que chaque fonction vh ∈ V0,h est déterminée de manière unique
par la donnée de ses valeurs aux points xi pour i = 2, . . . , N. L’espace V0,h est

4
Mini Projet 1 Master MASI

de dimension N − 1 et il est engendré par la base de Lagrange qui est formée

des N − 1 fonctions wi ∈ V0,h définies par :

1 si i = j,
wi (xj ) = δij = ∀i = 2, . . . , N et ∀j = 2, . . . , N
0 sinon
D’où :  x−xi−1

 xi −xi−1 si x ∈ [xi−1 , xi ]
xi+1 −x
wi (x) = xi+1 −xi si x ∈ [xi , xi+1 ]

 0 si x ∈
/ [xi−1 , xi+1 ]
Écriture du problème approché
Écrivons le problème approché dans V0,h : Trouver la fonction uh appartenant
à V0,h telle que pour tout vh ∈ V0,h on ait :
Z b Z b Z b
αu0h (x)vh0 (x)dx + βuh (x)vh (x)dx = f (x)vh (x)dx (2)
a a a

Le problème étant linéaire, l’égalité est vraie pour tout vh si et seulement si

elle est vraie pour une base de l’espace vectoriel V0.h , c’est-à-dire
(2) est vérifiée ∀vh ∈ V0,h ⇔ (2) est vérifiée ∀wi pour i = 2, . . . , N .
D’autre part, si uh est solution du problème approché dans V0,h , on peut
l’exprimer dans la base des wi comme suit :
N
X
uh (x) = uj wj (x)
j=2

avec uj = uh (xj ) est le valeur approchée de la solution exacte au point xj . Ainsi

le problème approché s’écrit : Chercher u2 , u3 , . . . , uN tels que
N Z
X b Z b  Z b
αwj0 (x)wi0 (x)dx + βwj (x)wi (x)dx uj = f (x)wi (x)dx
j=2 a a a

Posons Rb
Fi := a f (x)wi (x)dx
Rb Rb
Aij := a αwj0 (x)wi0 (x)dx + a βwj (x)wi (x)dx.
On remarque que le problème approché prend la forme d’un système linéaire de
N − 1 équations à N − 1 inconnues, qui peut s’écrire sous la forme matricielle
suivante :
AU = F
Le calcule des coefficients de la matrice de raideur, de masse et du second
membre plus de l’assemblage se trouve dans le Script.

5 FS Agadir (2020-2021)
 Mini Projet 1 Master MASI

On a alors le programme suivants en MATLAB :

1 function [ ] = EFdegre1 ( )
2 %% C l e a r t h e workspace
3 clc ; clear all ;
4 % Equation p a r a m e t e r s
5 alpha = 1 ; beta = 1 ;
6 N = 1 6 ; h = 1/N;
7 a = 0; b = 1;
8 % C a l c u l e d e s M a t r i c e s é l é m e n t a i r e s :
9 K = ( a l p h a /h ) ∗ [ 1 , −1; −1 , 1 ] ;
10 M = ( b e t a /3∗h ) ∗ [ 1 , 0 . 5 ; 0 . 5 , 1 ] ;
11 % Maillage
12 I 1 = l i n s p a c e ( a , b ,N+1) ;
13 f u n c t i o n y=f ( x )
14 y=x ˆ 4 ;
15 end
16 % Approximated s o l u t i o n
17 Uapp = s o l u t i o n d e g r e 1 (N, h , I 1 , @f ,K,M) ;
18 % Exacte s o l u t i o n
19 c2 = (37 −24∗ exp ( 1 ) ) / ( exp ( 1 )−exp ( −1) ) ;
20 c1 = −24−c2 ;
21 Uex = c1 ∗ exp ( I 1 ) + c2 ∗ exp(− I 1 )+ I 1 . ˆ 4 + 12∗ I 1 . ˆ 2 + 24 ;
22 % Error
23 e r r = Uapp−Uex ’ ;
24 e r r o r = max( abs ( e r r ) ) ;
25 disp ( ’ erreur ’ )
26 disp ( error )
27 % Ploting
28 I 1 = l i n s p a c e ( a , b , 200 ) ;
29 I 2 = l i n s p a c e ( a , b ,N+1) ;
30 s = c1 ∗ exp ( I 1 ) + c2 ∗ exp(− I 1 )+ I 1 . ˆ 4 + 12∗ I 1 . ˆ 2 + 24 ;
31 h o l d on
32 royalBlue = [0 ,0.4470 , 0 . 7 4 1 0 ] ;
33 r oyalGray = 1 / 2 5 5 ∗ [ 2 0 0 , 2 0 0 , 2 0 0 ] ;
34 p l o t ( I 1 , s , ’ c o l o r ’ , r o y a l B l u e , ’ LineWidth ’ , 2 ) ;
35 p l o t ( I 2 , Uapp , ’ s ’ , ’ MarkerEdge ’ , ’ k ’ , ’ MarkerFace ’ , royalGray , ’ M a r k e r S i z e
’ ,4) ;
36 hold o f f
37 end
38 % Assemblage de l a m a t r i c e g l o b a l e B
39 f u n c t i o n [ B ] = m a t r i x A d e g r e 1 (N, K, M)
40 B = z e r o s (N+1,N+1) ;
41 f o r i = 1 :N
42 for j = 1:2
43 m = i+j − 1 ;
44 for k = 1:2
45 n = i+k − 1 ;
46 B( m , n ) = B( m , n )+K( j , k )+M( j , k ) ;
47 end
48 end
49 end

6 FS Agadir (2020-2021)
 Mini Projet 1 Master MASI

50 end
51 % Formule de t r a p e z e
52 f u n c t i o n [ ap ] = a p r o x i m a t i o n t r a p e z e s ( f )
53 ap = 1 / 2 ∗ ( f ( 0 )+f ( 1 ) ) ;
54 end
55 % Les f o n c t i o n s de forme Phi1 e t Phi2
56 f u n c t i o n y = phi1 ( x )
57 y = 1−x ;
58 end
59

60 f u n c t i o n y = phi2 ( x )
61 y = x ;
62 end
63 % Le v e c t e u r s e c o n d membre
64 f u n c t i o n [ l ] = member2 s ( e l F i n i s , h , s , f )
65 x1 = e l F i n i s ( s ) ;
66 l=z e r o s ( 2 , 1 ) ;
67 f u n c t i o n y = r1 ( x )
68 y = f ( x1+h∗x ) ∗ p h i 1 ( x ) ;
69 end
70 f u n c t i o n y = r2 ( x )
71 y = f ( x1+h∗x ) ∗ p h i 2 ( x ) ;
72 end
73 l ( 1 ) = a p r o x i m a t i o n t r a p e z e s ( @r1 ) ;
74 l ( 2 ) = a p r o x i m a t i o n t r a p e z e s ( @r2 ) ;
75 end
76

77 f u n c t i o n [ F ] = vectMember2 degre 1 ( e l F i n i s , h ,N, f )

78 F = z e r o s (N+1 , 1 ) ;
79 f o r s = 1 :N
80 l = member2 s ( e l F i n i s , h , s , f ) ;
81 for i = 1:2
82 i g = s+i −1;
83 F( i g ) = F( i g )+h∗ l ( i ) ;
84 end
85 end
86 end
87 % C a l c u l e du s o l u t i o n a p p r o c hé e
88 f u n c t i o n [ U ] = s o l u t i o n d e g r e 1 (N, h , I , f ,K,M)
89 A = m a t r i x A d e g r e 1 (N, K,M) ;
90 A( 1 , : ) = 0 ; A( 1 , 1 ) =1;
91 A(N+1 , : ) = 0 ; A(N+1,N+1)=1;
92 F = vectMember2 degre 1 ( I , h ,N, f ) ;
93 F(1) = 0 ;
94 F(N+1) = 0 ;
95 U = A\F ;
96 end

7 FS Agadir (2020-2021)
Mini Projet 1 Master MASI

L’exécution de ce Script nous donne la figure suivante

Figure 2.1 – La solution exacte + approchée

L’erreur commise vaut

1.352671032083912e − 04

8 FS Agadir (2020-2021)
 Chapitre 3

Éléments Finis de degré 2 pour le

problème de Dirichlet homogène en
dimension n = 1

Commentaires
• On considère toujours le même problème, cette fois on utilise un polynôme
de degré deux, ce qui nécessite 3 fonctions de forme chacune de degré deux. Et
les matrices seront de forme 3 × 3.
• Cette fois le maillage nous donne la table de connectivité avec l’ajout d’un
point intermédiaire dans chaque élément Ti .
• En utilisant les polynômes de degré deux, on se trouve dans l’obligation
de travailler avec trois points ce qui veut dire trois fonctions de formes.

Le calcule des coefficients de la matrice de raideur, de masse et du second

membre plus de l’assemblage se trouve dans le Script suivant.

1 f u n c t i o n [ ] = EFdegre2 ( )
2 %% C l e a r t h e workspace
3 clc ; clear all ;
4 % Equation p a r a m e t e r s
5 a l p h a =1; b e t a =1;
6 a =0; b=1;
7 N=16; h=1/N;
8 % M a t r i c e é l e m e n t a i r e
9 K=(2∗ a l p h a /h ) ∗ [ 7 / 6 , −4/3 , 1 / 6 ; −4/3 , 8 / 3 , −4/3;1/6 , −4/3 , 7 / 6 ] ;
10 M=( b e t a /2∗h ) ∗ [ 4 / 1 5 , 2 / 1 5 , −1/15;2/15 , 1 6 / 1 5 , 2 / 1 5 ; −1/15 , 2 / 1 5 , 4 / 1 5 ]
;
11 % Maillage
12 I 1=l i n s p a c e ( a , b , 2∗N+1) ;
13 I 2=l i n s p a c e ( a , b , 2 ∗N+1) ;
14

9
 Mini Projet 1 Master MASI

15 f u n c t i o n y=f ( x )
16 y=x ˆ 4 ;
17 end
18

19 % Approximated s o l u t i o n
20 Uapp=s o l u t i o n d e g r e 2 (N, h , I 2 , @f ,K,M) ;
21 % Exacte s o l u t i o n
22 c2 = (37 −24∗ exp ( 1 ) ) / ( exp ( 1 )−exp ( −1) ) ;
23 c1 =−24−c2 ;
24 Uexact=c1 ∗ exp ( I 1 ) + c2 ∗ exp(− I 1 )+ I 1 . ˆ 4 + 12∗ I 1 . ˆ 2 + 24 ;
25 % Error
26 e r r=Uapp−Uexact ’ ;
27 disp ( ’ erreur ’ )
28 d i s p (max( abs ( e r r ) ) )
29 % Plot
30 I 1=l i n s p a c e ( a , b , 200 ) ;
31 s=c1 ∗ exp ( I 1 ) + c2 ∗ exp(− I 1 )+ I 1 . ˆ 4 + 12∗ I 1 . ˆ 2 + 24 ;
32 h o l d on
33 royalBlue = [0 ,0.4470 , 0 . 7 4 1 0 ] ;
34 r oyalGray = 1 / 2 5 5 ∗ [ 2 0 0 , 2 0 0 , 2 0 0 ] ;
35

36 p l o t ( I 1 , s , ’ c o l o r ’ , r o y a l B l u e , ’ LineWidth ’ , 2 ) ;
37 p l o t ( I 2 , Uapp , ’ s ’ , ’ MarkerEdge ’ , ’ k ’ , ’ MarkerFace ’ , royalGray , ’ M a r k e r S i z e
’ ,4) ;
38 hold o f f
39 end
40 % Assemblage de l a m a t r i c e g l o b a l e B
41 f u n c t i o n [ ms]= m a t r i x A d e g r e 2 (N, K,M)
42 ms=z e r o s ( 2 ∗N+1 ,2∗N+1) ;
43 f o r i e =1:N
44 f o r i =1:3
45 i g =2∗ i e+i −2;
46 f o r j =1:3
47 j g =2∗ i e+j −2;
48 ms( i g , j g )=ms( i g , j g )+K( i , j )+M( i , j ) ;
49 end
50 end
51 end
52 end
53 % Approximation de l ’ i n t e g r a l e de f
54 f u n c t i o n [ ap]= a p r o x i m a t i o n 2 p o i n t s ( g )
55 w=[ 1 , 1 ] ;
56 xpq=[− s q r t ( 1 / 3 ) , s q r t ( 1 / 3 ) ] ;
57 gpq =[g ( xpq ( 1 ) ) , g ( xpq ( 2 ) ) ] ;
58 ap =0;
59 f o r i =1:2
60 ap=ap+w( i ) ∗ gpq ( i ) ;
61 end
62 end
63 % Les f o n c t i o n s de forme Phi1 , Phi2 e t p h i 3
64 f u n c t i o n y=p h i 1 d 2 ( x )
65 y=1/2∗x ∗ ( x−1) ;
66 end

10 FS Agadir (2020-2021)
 Mini Projet 1 Master MASI

67

68 f u n c t i o n y=p h i 2 d 2 ( x )
69 y=(1−x ) ∗(1+x ) ;
70 end
71

72 f u n c t i o n y=p h i 3 d 2 ( x )
73 y=1/2∗x ∗ ( x+1) ;
74 end
75 % Le v e c t e u r s e c o n d membre
76 f u n c t i o n [ l ]= m e m b e r 2 d e g r e 2 s ( e l F i n i s , h , s , f )
77 x2=e l F i n i s ( 2 ∗ s ) ;
78 l=z e r o s ( 3 , 1 ) ;
79 f u n c t i o n y=r 1 ( x )
80 y=f ( x2+h/2∗ x ) ∗ p h i 1 d 2 ( x ) ;
81 end
82 f u n c t i o n y=r 2 ( x )
83 y=f ( x2+h/2∗ x ) ∗ p h i 2 d 2 ( x ) ;
84 end
85 f u n c t i o n y=r 3 ( x )
86 y=f ( x2+h/2∗ x ) ∗ p h i 3 d 2 ( x ) ;
87 end
88 l ( 1 )=a p r o x i m a t i o n 2 p o i n t s ( @r1 ) ;
89 l ( 2 )=a p r o x i m a t i o n 2 p o i n t s ( @r2 ) ;
90 l ( 3 )=a p r o x i m a t i o n 2 p o i n t s ( @r3 ) ;
91 end
92

93 f u n c t i o n [ F]= vectMember2 degre 2 ( e l F i n i s , h ,N, f )

94 F=z e r o s ( 2 ∗N+1 , 1 ) ;
95 f o r s =1:N
96 l=m e m b e r 2 d e g r e 2 s ( e l F i n i s , h , s , f ) ;
97 f o r i =1:3
98 i g =2∗ s+i −2;
99 F( i g )=F( i g )+h/2∗ l ( i ) ;
100 end
101 end
102 end
103 % C a l c u l e du s o l u t i o n a p p r o c hé e
104 f u n c t i o n [U]= s o l u t i o n d e g r e 2 (N, h , e l F i n i s , f ,K,M)
105 A=m a t r i x A d e g r e 2 (N, K,M) ;
106 A( 1 , : )=0 ;A( 1 , 1 ) =1;
107 A( 2 ∗N+1 , : ) =0;A( 2 ∗N+1 ,2∗N+1)=1;
108 F=vectMember2 degre 2 ( e l F i n i s , h ,N, f ) ;
109 F ( 1 ) =0;
110 F( 2 ∗N+1)=0;
111 U=A\F ;
112 end

11 FS Agadir (2020-2021)
Mini Projet 1 Master MASI

L’exécution de ce Script nous donne la figure suivante

Figure 3.1 – La solution exacte + approchée

L’erreur commise vaut

4.543847223264716e − 07

Conclusion

Conclusion Nous avons pu, dans ce travail, simuler numériquement le

problème de Dirichlet homogène à l’aide d’une approche dite Éléments
finis. La variation du maillage, ou du second membre permet aussi d’ob-
tenir différentes résultats. Le prolongement de ce code serait d’étendre
son champ d’application à des des dimensions supérieurs. Il serait, par
exemple intéressant d’introduire une modélisation de configurations de
contrôle non destructif. Ce type de modélisation représenterait plus
fidèlement la réalité.

12 FS Agadir (2020-2021)