UVSQ - Maı̂trises de Physique et Mécanique - MP3 16 juin 2003

Examen de Méthodes Numériques

Approximation de Solutions
16 juin 2003 — Durée 1 heure 30
tous documents autorisés

Approximation d’un problème d’électrostatique 1D

On recherche une approximation du potentiel électrique V (x)
dans un milieu diélectrique isotrope sous l’hypothèse électrosta- V0 D0
tique. Les diverses symétries permettent de considérer le problème
comme unidimensionnel. La permittivité électrique du matériau ρ
est variable et notée (x). Le domaine considère d’épaisseur L
contient une distribution de charges volumiques notée ρ(x). Il est ε
x
soumis sur son bord x = 0 à un potentiel imposé nul V0 = 0 et
sur son bord x = L à un déplacement électrique non nul connu D0 0 L
(voir figure). Problème électrostatique 1D

Les équations de Maxwell pour ce problème 1D se ramène à :

 
∂ ∂V (x)
(x) + ρ(x) = 0, ∀x ∈ [0, L] (1)
∂x ∂x
V (0) = 0, en x = 0 (2)
∂V (x)
(L) + D0 = 0, en x = L (3)
∂x |L

1. Etablir la formulation intégrale forte de ce problème. On note P (x) la fonction de projection

dans le domaine et P (x) celle utilisée pour les conditions aux limites.

2. En déduire la formulation intégrale faible. Proposer des relations entre P (x) et P (x) qui permette
de simplifier au maximum cette formulation. Quel est l’intérêt de cette formulation faible.

3. On cherche une approximation du potentiel solution sous la forme :

N
X
V (x) = qi φi (x) (4)
i=0

On applique la méthode de Galerkin à la formulation faible établie dans la question précédente.

On considère que la base d’approximation choisie est telle que la condition en x = 0 est automa-
tiquement vérifiée, c’est à dire que les fonctions de base sont telles que φi (0) = 0. Donner la
forme générale des termes de la matrice de raideur et du vecteur des forces généralisées.

4. On souhaite obtenir une approximation polynomiale quadratique de la solution dans le cas

particulier où (x) = cste =  et ρ(x) = ax (où a est donné). Pour cela, on choisi :
x x
φ1 (x) = , φ2 (x) = ( )2
L L
Construire la solution approchée. Est-ce la solution exacte?
1
5. On souhaite construire les caractéristiques d’un élément fini linéaire pour ce problème. Proposer
des fonctions de base linéaires pour cet élément. Donner l’expression de la matrice de rigidité
élémentaire pour le cas particulier de la question précédente.
 UVSQ - Maı̂trises de Physique et Mécanique - MP3 Juin 2003 – Eléments de correction
Eléments de correction
1. La formulation intégrale forte du problème est : trouver le potentiel V (x) tel que
Z L  
∂ ∂V (x) ∂V
[ (x) + ρ(x)]P (x)dx + V (0)P (0) + ((L) |L + D0 )P (L) = 0.
0 ∂x ∂x ∂x
où P (x) et P (x) sont les fonctions de projection.
2. Pour obtenir la formulation faible, on réalise une intégration par partie du premier terme. La
formulation du problème devient : trouver le potentiel V (x) tel que
Z L  L
∂V ∂P ∂V
− (x) dx + (x) P (x) +
0 ∂x ∂x ∂x 0
Z L
∂V
ρ(x)P (x)dx + V (0)P (0) + ((L) |L + D0 )P (L) = 0, ∀P (x), P (x)
0 ∂x
soit
Z L
∂V ∂P ∂V ∂V
− (x) dx + (L) |L P (L) − (0) |0 P (0)+
0 ∂x ∂x ∂x ∂x
Z L
∂V
ρ(x)P (x)dx + V (0)P (0) + ((x) |L + D0 )P (L) = 0, ∀P (x), P (x)
0 ∂x
En faisant le choix particulier suivant pour les fonctions de projection :
P (L) = −P (L) et P (0) = P (0)
la formulation se simplifie en : trouver le potentiel V (x) tel que
Z L Z L  
∂V ∂P ∂V
(x) dx − ρ(x)P (x)dx + (0) |0 − V (0) P (0) + D0 P (L) = 0
0 ∂x ∂x 0 ∂x

3. Pour appliquer la méthode de Galerkin, on choisit une base de fonctions φi (x), telle que φi (0) = 0.
Les fonctions de base sont choisies comme fonctions de projection. La formulation devient donc :
trouver le potentiel approximé V (x) tel que
Z L Z L
∂V ∂φj
(x) dx − ρ(x)φj (x)dx + D0 φj (L) = 0, ∀j ∈ [0, . . . , N ]
0 ∂x ∂x 0

En cherchant des approximations dans la base, le problème devient : trouver les composantes qi
de l’approximation V (x) = N
P
i=0 i φi (x) telles que :
q
Z L Z L
∂φi ∂φj
(x)qi dx − ρ(x)φj (x)dx + D0 φj (L) = 0, ∀j ∈ [0, . . . , N ]
0 ∂x ∂x 0
Il s’agit bien d’un système linéaire sur les qi :
 
 q0
 

[K]{q} = {f }, avec q= ..
 . 
qN
 

où les termes de la matrice de rigidité [K] s’espriment sous la forme :

Z L
∂φi ∂φj
kij = (x) dx
0 ∂x ∂x
2
et les termes du vecteur forces généralisées {f } sous la forme :
Z L
fi = ρ(x)φi (x)dx − D0 φi (L) = 0
0
 UVSQ - Maı̂trises de Physique et Mécanique - MP3 Juin 2003 – Eléments de correction
x
4. Dans le cas paticulier (x) = cste =  et ρ(x) = ax, et pour les fonctions de base : φ1 (x) = L et
φ2 (x) = ( Lx )2 , on obtient :
L L
x2 aL2
Z Z
 
k11 = dx = f1 = dx − D0 .1 =
a − D0
0 L2 L 0 L 3
L Z L 3
aL2
Z
  x
k12 = 2Xdx = f2 = a 2 dx − D0 .1 = − D0
0 L3 L 0 L 4
Z L
 4
k11 = 4
4x2 dx =
0 L 3L

qui conduit au système :

aL3 D0 L

 q1 + q2

 = −
3 


3
 q1 + 4 q2 = aL − D0 L



3 4 
dont la solution est :
7 aL3 D0 L


 q 1 = −
12   7 aL2 D0

 aL 2
soit V (x) = ( − )x − x
3 12   4
 q2 = − aL



4

Il ne s’agit bien évidemment pas de la solution exacte. On peut voir en particulier que :
∂ ∂V aL
( )=− 6= −ax
∂x ∂x 2

5. Dans la cas d’une base de fonctions de type éléments finis linéaires, les fonctions de base utilisées
sur un domaine 1D de longeur L sont :
x x
φ1 (x) = 1 − et φ2 (x) =
L L
qui donnent pour les termes de la matrice de rigidité :
Z L
 
k11 = 2
dx =
0 L L
Z L
 
k12 = − 2 dx = −
0 L L
Z L
 
k22 = 2
dx =
0 L L

soit pour la matrice :

 
 1 −1
[K] =
L −1 1
3
 UVSQ - Maı̂trises de Physique et Mécanique - MP3 2 septembre 2003

Examen de Méthodes Numériques

Approximation de Solutions
2 septembre 2003 — Durée 1 heure
tous documents autorisés

Approximation de l’équation de Laplace 1D

On recherche une approximation de la fonction scalaire u(x) obéissant à l’équation de Laplace sur un
domaine [0, h] et à diverses conditions aux limites.

Les équations de ce problème 1D sont :

∂ 2 u(x)
= f (x), ∀x ∈ [0, h] (1)
∂x2
u(0) = 0, en x = 0 (2)
∂u(x)
= a, en x = h (3)
∂x |h

où f (x) est une fonction scalaire connue et a un scalaire connu.

1. Etablir la formulation intégrale forte de ce problème. On note P (x) la fonction de projection

dans le domaine et P (x) celle utilisée pour les conditions aux limites.

2. En déduire la formulation intégrale faible. Proposer des relations entre P (x) et P (x) qui per-
mettent de simplifier au maximum cette formulation. Quel est l’intérêt de cette formulation
faible.

3. On cherche une approximation de la solution sous la forme :

N
X
u(x) = qi φi (x) (4)
i=0

On applique la méthode de Galerkin à la formulation faible établie dans la question précédente.

On considère que la base d’approximation choisie est telle que la condition en x = 0 est automa-
tiquement vérifiée, c’est à dire que les fonctions de base sont telles que φi (0) = 0. Donner la
forme générale des termes de la matrice de raideur et du vecteur des forces généralisées.

4. On souhaite obtenir une approximation polynomiale quadratique de la solution dans le cas

particulier où f (x) = bx (où b est un scalaire donné). Pour cela, on choisit :
x x
φ1 (x) = , φ2 (x) = ( )2
h h
Construire la solution approchée. Est-ce la solution exacte?

4
 UVSQ - Maı̂trises de Physique et Mécanique - MP3 Sept 2003 – Eléments de correction
Eléments de correction
1. La formulation intégrale forte du problème est : trouver la fonction u(x) telle que
Z h 2
∂ u(x) ∂u
[ 2
− f (x)]P (x)dx + (u(0) − 0)P (0) + ( |h − a)P (h) = 0.
0 ∂x ∂x
où P (x) et P (x) sont les fonctions de projection.
2. Pour obtenir la formulation faible, on réalise une intégration par partie du premier terme. La
formulation du problème devient : trouver la fonction u(x) telle que
Z h  h
∂u ∂P ∂u
− dx + P (x) −
0 ∂x ∂x ∂x 0
Z h
∂u
f (x)P (x)dx + u(0)P (0) + ( |h − a)P (h) = 0, ∀P (x), P (x)
0 ∂x
soit
Z h
∂u ∂P ∂u ∂u
− dx + |h P (h) − |0 P (0)−
0 ∂x ∂x ∂x ∂x
Z h
∂u
f (x)P (x)dx + u(0)P (0) + ( |h − a)P (h) = 0, ∀P (x), P (x)
0 ∂x
En faisant le choix particulier suivant pour les fonctions de projection :
P (h) = −P (h) et P (0) = P (0)
la formulation se simplifie en : trouver la fonction u(x) telle que
Z h Z h  
∂u ∂P ∂u
dx + f (x)P (x)dx + |0 − u(0) P (0) − aP (h) = 0
0 ∂x ∂x 0 ∂x

3. Pour appliquer la méthode de Galerkin, on choisit une base de fonctions φi (x), telle que φi (0) = 0.
Les fonctions de base sont choisies comme fonctions de projection. La formulation devient donc :
trouver la fonction approximée u(x) telle que
Z h Z h
∂u ∂φj
dx + f (x)φj (x)dx − aφj (h) = 0, ∀j ∈ [0, . . . , N ]
0 ∂x ∂x 0

En cherchant des approximations dans la base, le problème devient : trouver les composantes qi
de l’approximation u(x) = N
P
i=0 φi (x) telles que :
q i
Z h Z h
∂φi ∂φj
qi dx + f (x)φj (x)dx − aφj (h) = 0, ∀j ∈ [0, . . . , N ]
0 ∂x ∂x 0
Il s’agit bien d’un système linéaire sur les qi :
 
 q0
 

[K]{q} = {f }, avec q= ..
 . 
qN
 

où les termes de la matrice de rigidité [K] s’expriment sous la forme :

Z h
∂φi ∂φj
kij = dx
0 ∂x ∂x
5
et les termes du vecteur forces généralisées {f } sous la forme :
Z h
fi = − f (x)φi (x)dx + aφi (h) = 0
0
 UVSQ - Maı̂trises de Physique et Mécanique - MP3 Sept 2003 – Eléments de correction
4. Dans le cas particulier f (x) = bx et pour les fonctions de base : φ1 (x) = hx et φ2 (x) = ( hx )2 , on
obtient :
Z h Z h 2
1 1 x bh2
k11 = 2
dx = f1 = − b dx + a.1 = − +a
0 h h 0 h 3
Z h Z h 3
1 1 x bh2
k12 = 3
2xdx = f2 = − b dx + a.1 = − +a
0 h h 0 h2 4
Z h
1 2 4
k22 = 4
4x dx =
0 h 3h

qui conduit au système :

bh3

 q1 + q2

 = ah −
3


3
 q1 + 4 q2 = ah − bh



3 4
dont la solution est :
7

 q = ah − bh3
 1

 12 7 2 bh
soit u(x) = (a − bh )x + x2
3 12 4
 q2 = bh



4

Il ne s’agit bien évidemment pas de la solution exacte. On peut voir en particulier que :

∂2u bh
= 6= bx
∂x2 2

5. Dans la cas d’une base de fonctions de type éléments finis linéaires, les fonctions de base utilisées
sur un domaine 1D de longueur h sont :
x x
φ1 (x) = 1 − et φ2 (x) =
h h
qui donnent pour les termes de la matrice de rigidité :
Z h
1 1
k11 = 2
dx =
0 h h
Z h
1 1
k12 = − 2 dx = −
0 h h
Z h
1 1
k22 = 2
dx =
0 h h

soit pour la matrice :

 
1 1 −1
[K] =
h −1 1

6
 UVSQ - Maı̂trises de Physique et Mécanique - MP3 14 juin 2004

Examen de Méthodes Numériques

Approximation de Solutions
14 juin 2004 — Durée 1 heure 30
tous documents autorisés

Approximation d’un problème de fil tendu

On cherche une approximation de la déformée w(x) d’un fil inextensible tendu (par une tension T
constante) et soumis à une densité linéique de force f (x). Ce fil de longueur L est bloqué à ses deux
extrémités (figure de gauche).

w(x)

T 0 L T T 0 L/3 2L/3 L T
x x

f(x) F

Les équations de ce problème 1D sont :

∂ 2 w(x)
T + f (x) = 0, ∀x ∈ [0, L] (1)
∂x2
w(0) = 0, en x = 0 (2)
w(x) = 0, en x = L (3)

1. Etablir la formulation intégrale forte de ce problème. On note P (x) la fonction de projection

dans le domaine et P (x) celle utilisée pour les conditions aux limites.
2. En déduire la formulation intégrale faible. Proposer des conditions sur P (x) et P (x) qui per-
mettent de simplifier au maximum cette formulation. Quel est l’intérêt de cette formulation
faible?
3. On cherche une approximation de la déformée solution sous la forme :
N
X
w(x) = qi φi (x) (4)
i=0
On applique la méthode de Galerkin à la formulation faible établie dans la question précédente.
On considère que la base d’approximation choisie est telle que les conditions en x = 0 et x = L
sont automatiquement vérifiées, c’est-à-dire que les fonctions de base sont telles que φi (0) = 0
et φi (L) = 0. Donner la forme générale des termes de la matrice de raideur et du vecteur des
forces généralisées.
4. On souhaite obtenir une approximation polynomiale quadratique (N = 2) de la solution dans le
cas particulier où f (x) = F pour L/3 < x < 2L/3 et f (x) = 0 sinon (figure de droite). Expliquer
pourquoi cela revient à ne faire l’approximation qu’avec la seule fonction :
φ2 (x) = x(x − L)
Construire la solution approchée. Est-ce la solution exacte?
7
5. On souhaite construire les caractéristiques d’un élément fini linéaire pour ce problème. Proposer
des fonctions de base linéaires pour cet élément. Donner l’expression de la matrice de rigidité
élémentaire. Proposer un exemple de maillage pour le problème particulier de la question 4.
 UVSQ - Maı̂trises de Physique et Mécanique - MP3 Juin 2004 – Eléments de correction
Eléments de correction
1. La formulation intégrale forte du problème est : trouver le déplacement w(x) tel que
Z L
∂ 2 w(x)
[T + f (x)]P (x)dx + w(0)P (0) + w(L)P (L) = 0
0 ∂x2
quelles que soient les fonctions de projection P (x) et P (x).
2. Pour obtenir la formulation faible, on réalise une intégration par partie du premier terme. La
formulation du problème devient : trouver le potentiel V (x) tel que
Z L  L
∂w ∂P ∂w
− T dx + T P (x) +
0 ∂x ∂x ∂x 0
Z L
f (x)P (x)dx + w(0)P (0) + w(L)P (L) = 0, ∀P (x), P (x)
0
soit
Z L
∂w ∂P ∂w ∂w
− T dx + T |L P (L) − T |0 P (0)+
0 ∂x ∂x ∂x ∂x
Z L
f (x)P (x)dx + V (0)P (0) + w(L)P (L) = 0, ∀P (x), P (x)
0

En faisant le choix particulier suivant pour les fonctions de projection :

P (L) = P (L) = P (0) = P (0) = 0
la formulation se simplifie en : trouver le déplacement w(x) tel que
Z L Z L
∂w ∂P
T dx − f (x)P (x)dx = 0, ∀P (x)
0 ∂x ∂x 0

3. Pour appliquer la méthode de Galerkin, on choisit une base de fonctions φi (x), telle que φi (0) = 0.
Les fonctions de base sont choisies comme fonctions de projection. La formulation devient donc :
trouver le déplacement approximé w(x) tel que
Z L Z L
∂w ∂φj
T dx − f (x)φj (x)dx = 0, ∀j ∈ [0, . . . , N ]
0 ∂x ∂x 0

En cherchant des approximations dans la base, le problème devient : trouver les composantes qi
de l’approximation w(x) = N
P
i=0 i φi (x) telles que :
q
Z L Z L
∂φi ∂φj
T qi dx − f (x)φj (x)dx = 0, ∀j ∈ [0, . . . , N ]
0 ∂x ∂x 0
Il s’agit bien d’un système linéaire sur les qi :
 
 q0
 

[K]{q} = {f }, avec q= ..
 . 
qN
 

où les termes de la matrice de rigidité [K] s’espriment sous la forme :

Z L
∂φi ∂φj
kij = T dx
0 ∂x ∂x
8
et les termes du vecteur forces généralisées {f } sous la forme :
Z L
fi = f (x)φi (x)dx
0
 UVSQ - Maı̂trises de Physique et Mécanique - MP3 Juin 2004 – Eléments de correction
4. Une approximation polynomiale de degré deux s’écrit :

w(x) = q0 + q1 x + q2 x2

La satisfaction des conditions aux limites en x = 0 et x = L impose :

q0 = 0 et q1 = −Lq2

L’approximation est donc :

w(x) = q2 x(L − x)

Dans le cas de chargement proprosée on obtient :

L 2L/3
T L3 13T L2
Z Z
k22 = T (2x − L)2 dx = f2 = F x(x − L)dx = −
0 3 L/3 162

qui conduit à l’équation :

T L3 13T L3
q2 = −
3 162
dont la solution est :
13F 13F
q2 = − soit w(x) = x(L − x)
54T 54T

Il ne s’agit bien évidemment pas de la solution exacte. On peut voir en particulier que :

∂2w 13F
T 2
=−
∂x 27T
ne vérifie pas l’équation (1).

5. Dans la cas d’une base de fonctions de type éléments finis linéaires, les fonctions de base utilisées
sur un domaine 1D de longeur h sont :
x x
φ1 (x) = 1 − et φ2 (x) =
h h
qui donnent pour les termes de la matrice de rigidité :
Z h
h T
k11 = 2
dx =
0 L h
Z h
h T
k12 = − 2 dx = −
0 L h
Z h
h T
k22 = 2
dx =
0 L h

soit pour la matrice :

 
T 1 −1
[K] =
h −1 1
9
 UVSQ - Maı̂trises de Physique et Mécanique - MP3 2 septembre 2004

Examen de Méthodes Numériques

Approximation de Solutions
2 septembre 2003
tous documents autorisés

Déplacement des spires d’un ressort pesant

On considère un ressort pesant (de masse linéı̈que µ et de

longueur L et de raideur par unité de longueur k) retient ω k, µ
M
une masse M pouvant coulisser sans frottement sur une tige. r
On cherche à approximer le déplacement u(r) de chacun des 0
points du ressort lorsque l’ensemble est entraı̂né en rotation L
à une vitesse angulaire fixe ω.
Les équations de ce problème 1D sont :
∂ 2 u(r)
k + µω 2 r = 0, ∀r ∈ [0, L] (1)
∂r2
u(0) = 0, en r = 0 (2)
∂u(r)
k = M ω 2 L, en r = L (3)
∂r |L
où k, L, µ, M et ω sont des scalaires fixés.

1. Etablir la formulation intégrale forte de ce problème. On note P (r) la fonction de projection

dans le domaine et P (r) celle utilisée pour les conditions aux limites.
2. En déduire la formulation intégrale faible. Quel est l’intérêt de cette formulation faible.
3. Pour chercher une approximation de ce problème, on souhaite utiliser la méthode de Ritz. On
l’applique au principe de minimum de l’énergie potentielle du ressort. C’est-à-àdire que le
prolème est formulé de la manière suivante : Le champ de déplacement solution est tel que
u(r = 0) = 0, et miminise l’énergie potentielle :
Z Z L
1 L ∂u(r) ∂u(r)
Ep (u) = k dr − µω 2 ru(r)dr − M ω 2 Lu(L)
2 0 ∂r ∂r 0

Et on cherche une approximation de la solution sous la forme :

N
X
u(r) = qi φi (r) (4)
i=0

On considère que la base d’approximation choisie est telle que la condition en r = 0 est automa-
tiquement vérifiée, c’est à dire que les fonctions de base sont telles que φi (0) = 0. Donner la
forme générale des termes de la matrice de raideur et du vecteur des forces généralisées. On
rappelle que le système linéaire est obtenu par minimisation de la fonctionnelle Ep par rapport
à chacune des variables qi .
4. On souhaite obtenir une approximation polynomiale quadratique de la solution. Pour cela, on
choisit N = 2 et : 10
r r
φ1 (r) = , φ2 (r) = ( )2
L L
Construire la solution approchée. Est-ce la solution exacte?
 UVSQ - Maı̂trises de Physique et Mécanique - MP3 Sept 2004 – Eléments de correction
Eléments de correction
1. La formulation intégrale forte du problème est : trouver la fonction u(r) telle que
Z L
∂ 2 u(r) ∂u
[k + µω 2 r]P (r)dr + (u(0) − 0)P (0) + (k |L − M ω 2 L)P (L) = 0.
0 ∂r2 ∂r

quelques soient les fonctions de projection P (r) et P (r).

2. Pour obtenir la formulation faible, on réalise une intégration par partie du premier terme. La
formulation du problème devient : trouver la fonction u(r) telle que
Z L  L
∂u ∂P ∂u
− k dr + k P (r) +
0 ∂r ∂r ∂r 0
Z L
∂u
µω 2 rP (r)dr + u(0)P (0) + (k |L − M ω 2 L)P (L) = 0, ∀P (r), P (r)
0 ∂r
soit
Z L
∂u ∂P ∂u ∂u
− k dx + k |L P (L) − k |0 P (0)+
0 ∂r ∂r ∂r ∂r
Z h
∂u
µω 2 rP (r)dr + u(0)P (0) + (k |L − M ω 2 L)P (L) = 0, ∀P (r), P (r)
0 ∂r

3. Pour appliquer la méthode de Ritz, on choisit une base de fonctions φi (r), telle que φi (0) = 0.
Les fonctions de base sont choisies comme fonctions de projection. La formulation devient donc :
trouver la fonction approximée u(qi ) qui minimise
Z L Z L
1 ∂φi (r) ∂φj (r)
Ep (qi ) = kqi qj dr − µω 2 rqi φi dr − M ω 2 Lqi φi (L)
2 0 ∂r ∂r 0

En exprimant le minimum de la fonctionnelle Ep par rapport à chacune des variables qi , ont

obtient :
Z L Z L
∂Ep ∂φi (r) ∂φj (r)
=( k dr)qj − µω 2 rφi dr − M ω 2 Lφi (L) = 0, ∀j ∈ [0, . . . , N ]
∂qi 0 ∂r ∂r 0

Il s’agit bien d’un système linéaire sur les qi :

 
 q0
 

[K]{q} = {f }, avec q = ..
 . 
 
qN

où les termes de la matrice de rigidité [K] s’expriment sous la forme :

Z L
∂φi ∂φj
kij = k dr
0 ∂r ∂r

et les termes du vecteur forces généralisées {f } sous la forme :

Z L
11
fi = µω 2 rφ 2
i (r)dr + M ω Lφi (L)
0
 UVSQ - Maı̂trises de Physique et Mécanique - MP3 Sept 2004 – Eléments de correction
4. Dans le cas particulier N = 2 et pour les fonctions de base : φ1 (r) = Lr et φ2 (r) = ( Lr )2 , on
obtient :
Z L Z L
k k r2 µL
k11 = 2
dr = f1 = − µω 2 dx + M ω 2 L.1 = Lω 2 ( + M)
0 L L 0 L 3
Z L Z L 3
k k 2r µL
k12 = 3
2rdr = f2 = − µω 2
dr + M ω 2 L.1 = Lω 2 ( + M)
0 L L 0 L 4
Z L
k 2 4k
k22 = 4
4r dr =
0 L 3L

qui conduit au système : 


 L2 ω 2 µL
 q1 + q2
 =
k
(
3
+ M)


 2 2
 q1 + 4 q2 = L ω ( µL + M )

3 k 4
dont la solution est :


 L2 ω 2 7µL

 1 q = (M + )
k 12 L2 ω 2 7µL µω 2 L3 2
soit u(r) = (M + )r − r

 2 3 k 12k 4k
 q2 = − µω L

4k

Il ne s’agit bien évidemment pas de la solution exacte. On peut voir en particulier que :

∂2u µω 2 L3
k = − 6= −µω 2 r
∂r2 2

12