TS - TP Physique n°11 Eric DAINI – Lycée Paul Cézanne – Aix en Provence - http://labotp.

org

RESOLUTION D'UNE EQUATION PAR LA METHODE D'EULER (Correction)

Objectifs • Mettre en oeuvre la méthode d'Euler pour déterminer une solution numérique d'une équation différentielle.
• Trouver le modèle le plus pertinent pour l'expression de la force de frottement fluide.

I. EXPERIENCE

• Une bille d'acier tombe verticalement dans de vlim

l'huile (voir photo). Une chronophotographie de
l'expérience et un pointage Aviméca ont permis
d'obtenir le graphe v(t) ci-contre.
• La vitesse limite vlim et le temps
caractéristique τ de la chute sont rappelés sur
le graphe.
vlim(exp) = 0,94 m.s-1
• Les caractéristiques de la bille d'acier sont : -1
- Masse : m = 4,08 g; Diamètre : d = 10 mm; τexp==0,085
vlim(exp) 0,94 m.s
s
τexp = 0,085 s
Volume V.
- Masse volumique de la bille : ρ = 7 796 kg.m τ
-3

- Masse volumique de l'huile : ρh = 920 kg.m .

-3
τ
-2
- Intensité de la pesanteur: g = 9,80 m.s

II MODELISATION DE LA CHUTE

• Selon le type de force de frottement que l'on considère, l'équation différentielle du mouvement sur la vitesse s'écrit :

dv dv
Cas n°1: f = k1.v + A1 v = B (1) Cas n°2 : f = k2.v2 + A2v2 = B (2)
dt dt

 ρ 
avec A = k et B = g 1 − h 
m  ρ 

 ρ   920 
1) Valeur de la constante B avec 3 chiffres significatifs: B = g 1 − h  = 9,80 × 1 − -2
 = 8,64 m.s .
 ρ   7796 
-1 dv
2) Si v = vlim = 0,94 m.s = cte alors =0:
dt
B 8,64
dans l'équation (1): A1.vlim = B ⇔ A1 = = = 9,19 s-1
vlim 0,94

B 8,64
dans l'équation (2): A2.v²lim = B ⇔ A 2 = = = 9,78 m-1
v² lim (0,94)²

3) Expressions numériques des deux équations différentielles:

dv dv
+ 9,19.v = 8,64 (1) + 9,78.v 2 = 8,64 (2)
dt dt

III RESOLUTION NUMERIQUE D'UNE EQUATION DIFFERENTIELLE SELON LA METHODE D'EULER

1) La méthode d'Euler
dv ∆v v(t + ∆t) − v(t)
• Principe : pour des intervalles de temps ∆t "petits" on peut faire l'approximation suivante: a(t) = ≈ =
dt ∆t ∆t
soit : v(t+∆t) = v(t) + a(t) × ∆t.
∆t est appelé pas d’itération de la méthode d’Euler. L’expression de a(t) dépend du modèle utilisé.
Dans la suite on choisit : ∆t = 0,010 s.
TS - TP Physique n°11 Eric DAINI – Lycée Paul Cézanne – Aix en Provence - http://labotp.org
2) Calculs "à la main"

• Cas n°1 : modélisation en f = k 1.v: Cas n°2 : modélisation en f = k 2.v²

La bille est lâchée sans vitesse initiale à la date t0 = 0,000 s .
dv
dv a) a2(t) = = B - A1.v = 8,64 – 9,78 × v(t)
a) a1(t) = = B - A1.v = 8,64 – 9,19 × v(t) dt
dt b)
b)
-1 -2 -1 -2
t(s) v(t) en (m.s ) a1(t) en (m.s ) t(s) v(t) en (m.s ) a2(t) en (m.s )
t0 = 0,000 8,640 0,000 t0 = 0,000 0,000 0,000
t1 = to+ ∆t = 0,010 7,846 0,086 t1 = to+ ∆t = 0,010 0,086 0,086
t2 = t1+ ∆t = 0,020 7,125 0,165 t2 = t1+ ∆t = 0,020 0,172 0,172
t3 =t2 + ∆t = 0,030 6,470 0,236 t3 =t2 + ∆t = 0,030 0,256 0,256
t4 =t3 + ∆t = 0,040 5,876 0,301 t4 =t3 + ∆t = 0,040 0,336 0,336

3) Résolution avec un tableur

• Les calculs fastidieux "à la main" peuvent être automatisés avec un ordinateur et un tableur, tel qu'Excel .

Cas n°1 : modélisation de f en f = k1.v

a) Formule de calcul pour calculer t1 dans la cellule A10: = A9+B$3

formule de calcul pour calculer v(t1) dans la cellule B10: = B9+(B$5-B$4*B9)*B$3

b) Graphe de v(t) obtenu par la méthode d'Euler et les points expérimentaux:

c) Le modèle en f = k1.v n'est pas satisfaisant car le graphe v(t) obtenu avec la méthode d'Euler est situé au–dessous des points
expérimentaux.
TS - TP Physique n°11 Eric DAINI – Lycée Paul Cézanne – Aix en Provence - http://labotp.org
Cas n°2 : modélisation de f en f = k2.v2

a) formule de calcul pour calculer t1 dans la cellule A10: = A9+B$3

formule de calcul pour calculer v(t1) dans la cellule B10: =B9+(B$5-B$4*B9^2)*B$3

b) graphe de v(t) obtenu par la méthode d'Euler et les points expérimentaux.

c) Le modèle f = k2.v2 n'est pas satisfaisant car le graphe v(t) obtenu avec la méthode d'Euler passe au-dessus de certains points
expérimentaux.

Cas n°3 : modélisation de f en f = k3 . vn

• Les deux modélisations précédentes n'étant pas satisfaisantes, on envisage une troisième possibilité : f = k3.vn. L'équation différentielle
dv
correspondante s'écrit : Cas n°3 : f = k3.vn + A3vn = B
dt

a) Pour le modèle f = k1.v, les points expérimentaux sont au-dessus du graphe v(t) tandis que pour le modèle f = k2.v² certains points
expérimentaux sont au-dessous du graphe v(t).
On peut donc conclure que n est compris entre 1 et 2 dans le modèle: f = k3.vn

dv
b) Pour v = vlim = 0,94 m.s-1 = cte alors =0:
dt
B
équation (3): A3.vnlim = B ⇔ A3 = n
v lim
c) formule de calcul pour A3 dans la cellule B4: = B$5/0,94^B6
formule de calcul pour calculer t1 dans la cellule A11: = A10+B$3
formule de calcul pour calculer v(t1) dans la cellule B11: = B10+(B$5-B$4*B10^B$6)*B$3

d) La valeur de n avec le curseur qui donne la modélisation la plus satisfaisante est: n = 1,5

e) Graphe de v(t) obtenu par la méthode d'Euler et les points expérimentaux avec n = 1,5.
TS - TP Physique n°11 Eric DAINI – Lycée Paul Cézanne – Aix en Provence - http://labotp.org
f) La modélisation de f la plus satisfaisante est: f = k3.v1,5
avec:
k3 = m.A3 = 4,08.10-3 × 9,48 = 3,87.10-2 (S.I)
finalement:
f = 3,87.10-2 × v1,5

g) Lorsque le régime permanent est atteint, v = vlim = 0,94 m.s-1 donc: f = 3,87.10-2 × 0,94 = 3,5.10-2 N

h) Poids de la bille: P = m.g = 4,08.10-3 × 9,80 = 4,0.10-2 N

f + PA en régime permanent
Poussé d'Archimède: PA = ρh × V × g = 920 × (4π × (5,0.10-3)3 / 3) × 9,80 = 4,7.10-3 N
f + PA = 3,5.10-2 + 4,7.10-3 = 3,97.10-2 N ≈ 4,0.10-2 N = P

Le résultat pouvait être prévu car en régime permanent le vecteur vitesse est constant donc: P + f + PA = 0
En projection suivant un axe vertical orienté vers le bas: P – f – PA = 0 ⇔ P = f + PA .

4) Influence du pas ∆t sur la résolution de l'équation différentielle par la méthode d'Euler

a) Graphes de v(t) pour ∆t = 0,010 s et pour ∆t = 0,100 s:

b) Pour ∆t > 0,100 s, le graphe v(t) s'éloigne nettement des points expérimentaux. Pour être applicable, la méthode d'Euler doit utiliser un
pas d'intégration ∆t "pas trop grand".

c) Pour une bonne modélisation ∆t = 0,010 s.

La constante de temps est : τ = 0,085 s

On a donc: τ / ∆t = 8,5 ⇔ ∆t = τ / 8,5 ≈ τ / 10.

Une condition d'utilisation de la méthode d'Euler sur ∆t est donc: ∆t ≤ τ / 10.
TS - TP Physique n°11 Eric DAINI – Lycée Paul Cézanne – Aix en Provence - http://labotp.org
IV APPLICATION DE LA METHODE D'EULER A LA CHARGE D'UN CONDENSATEUR

• On considère un circuit (RC) série alimenté par un générateur de tension constante E.

1) Schéma du circuit (RC) et équation différentielle sur la tension uc(t) aux bornes du condensateur

i
Loi d'additivité des tensions: E = R.i(t) + uc(t)
+ Or i(t) = C.duc / dt
R
Donc: E = R.C. duc / dt + uc(t)
E En posant τ = RC et en divisant tous les termes par τ il vient:
-
C uc du c u
+ c=
E
dt τ τ

2) Exprimons uc(t+∆t):
du c ∆u c
Pour des intervalles de temps ∆t "petits" on peut faire l'approximation suivante: ≈
dt ∆t
du c ∆u c E u  E − uc 
On a: ≈ = - c or ∆uc = uc(t+∆t) – uc(t) donc il vient: uc(t + ∆t) = uc(t) +   × ∆t
dt ∆t τ τ  τ 

• Dans la cellule B10 on tape: =B9+H$5*(H$3-B9)/H$4

• Sachant qu'une solution de l'équation différentielle est uc(t) = E.(1 – exp(-t/τ)), on tape dans la cellule C10 de la colonne uc(t) théo:
=H$3*(1-EXP(-A10/H$4))

3) Allure des deux graphes:

5,0
uc(V)
4,5

4,0

3,5

3,0

uc(t) Euler
2,5 uc(t) théo

2,0

1,5

1,0

0,5

t(s)
0,0
0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9

On constate que les deux graphes sont pratiquement superposés.

4) Vérifions que la condition ∆t << τ est valable: ∆t = 0,010 s et τ = R.C = 0,10 s donc ∆t << τ avec ∆t = τ / 10.
si ∆t ≈ τ ou ∆t > τ les graphes ne se superposent plus.