Cours de mcanique

M12-Chute libre avec frottements

Introduction

Nous avons modlis au chapitre prcdent le corps qui chute dans le champ de pesanteur en
considrant que les frottements de lair taient ngligeables. Cette supposition nayant quune
utilit thorique, nous complexifions ici notre modle en tenant compte de ces frottements :
comment ceux-ci vont modifier la trajectoire du corps qui chute ?
Ce sera loccasion de voir que ces frottements peuvent tre de deux types, linaires ou quadratiques, nous avons alors rencontr deux types dquations diffrentielles : la rsolution de la
premire ne nous posera pas de problme ; mais la rsolution de la seconde est moins aise :
nous en profiterons pour voir une mthode numrique itrative permettant dapprocher la forme
de la solution : la mthode dEuler.
Enfin parmi les deux modles de forces de frottement, lequel se rvle le plus juste pour tudier
le parachutisme ? Nous tenterons une rponse laide de la mcanique des fluides.

Problme 3

Un parachutiste de masse 80 kg ralise un saut depuis un hlicoptre. La premire partie du

saut, celle qui nous intresse ici, est ralise sans parachute. Quelles sont les caractristiques de
celle-ci sachant que les frottements de lair ne sont pas ngligeables ?

Systme
Le systme tudi est le sauteur considr ponctuel.

Rfrentiel et base

On tudie son mouvement dans un rfrentiel terrestre li au sol ( son point de chute), ce
rfrentiel est considr galilen pendant la dure de la chute. On utilisera une base cartsienne
une dimension pour suivre lvolution du sauteur : un axe Oz vertical ascendant avec origine
au point de chute constituera le repre dtude.
On considre en effet que le mouvement du parachutiste est strictement vertical.

Forces

5.1

Bilan des forces

Le sauteur est soumis son poids not P , force distance exerce par la Terre sur lui.
Il est soumis aux forces de frottements de lair, modlises par une force de contact note

f . Cette force peut aussi tre nomme rsistance de lair.

Mcanique

5.2

M12-Chute libre avec frottements

5.2

Deux types de forces de frottements

Deux types de forces de frottements

La force de frottements de lair peut prendre deux formes :

Frottements linaires
Dans le cas dune vitesse faible, la force de frottement est proportionnelle la vitesse :

f = k
v

(1)

On parle de frottements linaires. k est une constante qui dpend de la nature du fluide et
des caractristiques de lobjet.
Par exemple pour une sphre de rayon r, on a k = 6 r o est la viscosit du fluide.
Frottements quadratiques
Dans le cas dune vitesse importante, la force de frottement est proportionnelle au carr de
la vitesse :

f = k 0 v
v
(2)
On parle de frottements quadratiques. k 0 est aussi une constante qui dpend du fluide et
des caractristiques de lobjet mais elle prend une autre forme que k :
Son expression est du type k 0 = 12 Cx S avec la viscosit du fluide, S la surface frontale de
lobjet et Cx le coefficient de traine (appel dans le langage courant coefficient de pntration
dans lair) qui dpend de la gomtrie du corps. Par exemple, voici trois gomtries et trois
valeurs de Cx :

Cx = 1.32

Cx = 0.45

Cx = 0.04

Ce coefficient de traine peut se calculer pour une sphre lisse (sans rugosit) dans le cas
dcoulement faible vitesse ( faible nombre de Reynolds), il dpend alors du nombre de
Reynolds.
Pour des coulements turbulents ( grand nombre de Reynolds > 103 ), on mesure le Cx en
soufflerie. En sachant quil est constant pour un corps donn.

Utilisation de la 2me loi de Newton

F = m
a P + f = m
a

(3)

Rsolution du problme dans le cas de frottements linaires

7.1

Equation diffrentielle

Le PFD donne : m
g k
v = m
a.

Mcanique

M12-Chute libre avec frottements

7.2

Solution de lquation diffrentielle

On projette maintenant cette relation sur laxe Oz vertical ascendant.

dvz
+ k vz = m g
dt
k
dvz
+ vz = g

dt
m

m g k vz = m a m

(4)
(5)

On obtient donc une quation diffrentielle en vz , linaire du premier ordre coefficients

constants. On sait rsoudre cette quation mathmatiquement.
Une fois lexpression de la vitesse vz obtenue, on en dduira la position par intgration.
Une notation particulire
Souvent ce type dquation sera crite ainsi :
dvz
vz
+
= g
dt

avec

m
k

La notation fait rfrence un temps. En effet, la grandeur =

de la fonction v = f (t), comme nous allons le voir par la suite.

7.2

m
est un temps caractristique
k

Solution de lquation diffrentielle

7.2.1

Principe

Une quation diffrentielle linaire avec second membre se rsout en deux temps :
on cherche dabord la solution sh de lquation homogne, cest dire lquation
sans second membre ;
on cherche une solution particulire sp , cest dire une solution qui a mme forme
que le second membre (si le second membre est constant, la solution particulire
recherche sera une constante).
La solution de lquation diffrentielle avec second membre est la somme de la solution
homogne et de la solution particulire : s = sh + sp .
Attention, dans la solution de lquation homogne apparaissent souvent des constantes
(une si lquation est du premier ordre, deux si elle est du deuxime ordre). La
dtermination de ces constantes laide des conditions initiales doit tre mene en tenant
compte de la solution particulire.

7.2.2

Pour notre problme

On recherche la solution de lquation complte (5), qui est une vitesse.

Solution de lquation homogne
dv v
Equation homogne :
+ =0
dt

t
Solution : vh = A exp

avec A une constante



On peut vrifier en drivant une fois vh que cette solution vrifie lquation homogne.

Mcanique

M12-Chute libre avec frottements

7.3

Courbe |vz | = f (t)

Solution particulire
Le second membre tant constant (gal g), on cherche une solution particulire vp = cste.
dvp
Alors
= 0 et on obtient vp = g .
dt
Solution globale
On a donc :

t
vz = A exp

+ g

(6)

On peut maintenant dterminer A laide des conditions initiales :

A t = 0 : v(t = 0) = 0 = A g A = g

(7)

Et finalement :


vz = g exp

(8)

Attention, rappelons que cette vitesse est ngative puisque le corps qui chute se dirige
suivant laxe Oz descendant.

7.3
7.3.1

Courbe |vz | = f (t) et caractristiques

Courbe

On souhaite visualiser la norme de la vitesse en fonction du temps. Son expression est donc :


|vz | = g 1 exp



(9)

Voici la courbe obtenue :

vlim

69.5

|vz |(m.s1 )

60

40
Cas des frottements linaires

20

10

20
t(s)

30

40

Figure 1 Vitesse du parachutiste dans le cas de frottements linaires

La vitesse augmente dabord fortement, puis de plus en plus faiblement pour atteindre une
valeur limite.

Mcanique

7.3.2

M12-Chute libre avec frottements

7.3

Courbe |vz | = f (t)

Vitesse limite

La valeur de la vitesse limite peut tre obtenue en calculant la limite de |vz (t)| quand le
temps tend vers linfini :
t
lim |vz (t)| = lim g 1 exp
t
t



= g

(10)

On peut galement la trouver partir de lquation diffrentielle (5). En effet, la vitesse

limite est constante, on a ainsi :
dvz lim vz lim
vz lim
+
= g 0 +
= g
dt

vz lim = g
|vz lim | = g

(11)
(12)
(13)

7.3.3

Temps caractristique
m
La grandeur =
est caractristique de lvolution de la vitesse dans le temps. Dans ce
k
type dvolution, on parle de rgime transitoire et de rgime permanent :
le rgime est transitoire tant que la vitesse volue ;
le rgime est permanent lorsque la vitesse limite est atteinte.
Le temps est un bon indicateur pour savoir quand on passe dun rgime lautre : on
considre quau bout de 5 , le rgime permanent est atteint.
Dtermination de
On peut obtenir la valeur de graphiquement : on cherche labscisse du point dintersection
entre la tangente la courbe en t = 0 et lasymptote quand t de la courbe |vz | = f (t).
On obtient ainsi la limite entre rgime transitoire et rgime permanent :
vlim

69.5

|vz |(m.s1 )

60

40
Rgime transitoire

Rgime permanent

20

10

20

30 5 40
t(s)

50

Figure 2 Diffrents rgimes lors de la chute avec frottements

Mcanique

7.4

M12-Chute libre avec frottements

7.4

Obtention de la position

(14)

Obtention de la position

La fonction z = f (t) sobtient en intgrant la fonction vz = f (t) :



vz = g exp

1 = g exp

=z(t) = g 2 exp

g t + cste

(15)

La constante est obtenue laide de la condition initiale de position :

A t = 0, z = h donc g 2 + cste = h cste = h + g 2
On a finalement :


z(t) = g 2 1 exp



g t+h

(16)

On prend comme altitude de dpart h = 4000 m.

0.5

4,000

0.4

z(m)

3,000
2,000
1,000
0

1 exp

5,000

0.3
0.2
0.1

10

20
t(s)

30

40

10

20
t(s)

30

40

Figure 3 Position du parachutiste dans le cas de frottements linaires

La forme de la courbe z = f (t) montre que la position varie quasi linairement par rapport
au temps, cest--dire quon a pratiquement z = a t +b avec
 a la pente ngative.
t
En effet, nous pouvons voir que lexpression 1 exp
est trs petite et que lexpression

de z(t) crite lquation (16) tend vers :

z(t) = g t + h

(17)

Il sagit bien dune droite de pente g ngative.

Rsolution dans le cas de frottements quadratiques

8.1

Equation diffrentielle

Le PFD donne : m
g k0 v
v = m
a.
Dans loptique dutiliser la mthode dEuler, nous allons utiliser un axe Oz vertical descendant afin de travailler avec une vitesse positive.
dvz
+ k 0 vz2 = m g
dt
dvz
k0

+ vz2 = g
dt
m

m g k vz2 = m a m

(18)
(19)

Mcanique

M12-Chute libre avec frottements

8.2

Vitesse limite

m
On pourra galement introduire le temps caractristique 0 = 0 . Cette quation diffrentielle
k
nest pas linaire, nous ne pouvons pas la rsoudre facilement.

8.2

Vitesse limite

Par contre, nous pouvons dores et dj connatre la vitesse limite :

dvz
Lorsque
= 0 alors vz lim =
dt

gm 0
= g .
k0

Cette quation diffrentielle, complexe rsoudre, va tre loccasion dutiliser une mthode
de rsolution numrique itrative : la mthode dEuler.

8.3

Rsolution par la mthode dEuler

8.3.1

Ce quest la mthode dEuler

La mthode dEuler est une mthode numrique itrative qui permet dobtenir une solution
approche dune quation diffrentielle partir des conditions initiales.
Rappels mathmatiques

Drive = cfficient directeur de la tangente la courbe.

Calcul dune drive en un point aise :
60

v(m.s1 )

v = vB vA
40
A
t = tB tA
20

10

15
t(s)

20

25

Figure 4 Calcul de la drive dune courbe en un point

Daprs la dfinition mathmatique de la drive :


dv
dt

=
t=10 s

Si on ralise un zoom sur la courbe :

v
t

(20)

Mcanique

M12-Chute libre avec frottements

60

8.3

Rsolution par la mthode dEuler

v(m.s1 )

v = vB vA
40
A

t = tB tA
20

v
t

10

15
t(s)

20

25

Figure 5 Drive et temps infinitsimal

On peut alors crire, en considrant un intervalle de temps t suffisamment petit :


dv
dt

=
t

v
t

(21)

On peut alors exprimer la petite variation de vitesse v qui se produit pendant le petit intervalle
de temps t grce lquation diffrentielle :
Si

dv
= A v2 + B
dt

alors

v = (A v 2 + B) t

lorsque t 0

(22)

Mise en uvre
On part de la condition initiale, la valeur de v(t = 0) = v0 ;
on choisit le pas de calcul, soit la valeur de t ;
on calcule :
v1 = v0 + v = v0 + (A v02 + B) t

(23)

et ainsi de suite :
vi+1 = vi + (A vi2 + B) t

(24)

un tableur viendra nous assister dans la rptition des calculs.

le choix du pas de calcul t doit tre judicieux : il faut prendre un intervalle suffisamment
petit pour que lapproximation soit valable, mais pas trop petit afin que les calculs ne
soient pas trop longs.
8.3.2

Utilisation de cette mthode dans notre cas

Obtention de la vitesse en fonction du temps

Pour utiliser cette mthode, il nous faut la valeur des cfficients A et B qui apparaissent
dans lquation diffrentielle :
dvz
k0
+ vz2 = g
dt
m
dvz
k0

= vz2 + g = A vz2 + B
dt
m

(25)
avec A =

k0
et B = g
m

(26)

La valeur de B est donc connu, on peut valuer la valeur de A partir de la vitesse limite
atteinte par un parachutiste. Celle-ci, qui dpend de la position (group ou droit comme un i)
du sauteur lors de la chute, est de 69,5 m.s1 . Donc :
r

vzlim =

gm
=
k0

g
g
9,81
= A =
=
= 2,0 102 SI
A
|vzlim |2
69,52
8

(27)

Mcanique

M12-Chute libre avec frottements

8.3

Rsolution par la mthode dEuler

On connat galement la vitesse initiale : vz (t = 0) = 0. On peut donc appliquer la mthode

en choisissant un pas t judicieux. On prendra par exemple t = 0,3 s. Alors :
v1 = v0 + (A v02 + B) t = B t = 2,94 m.s1
v2 = v1 + (A v12 + B) t = 5,87 m.s1
...
A laide dun tableur, on rpte les calculs jusquau temps voulu. On peut ensuite tracer la
courbe vz = f (t).
Ci-dessous, on a trac les courbes pour des pas de calculs diffrents. On remarque quil ny
a pas de diffrences entre nos trois tests.

Mcanique

M12-Chute libre avec frottements

9. Quel type de frottements choisir ?

vlim

69.5

vz (m.s1 )

60
Cas des frottements

40

quadratiques

t = 0.3
t = 0.7
t = 0.1

20

10

20
t(s)

30

40

Figure 6 Mthode dEuler applique aux cas des frottements quadratiques

Quen est-il de la position en fonction du temps ?
Pour obtenir la courbe de position en fonction du temps, on part de la donne de vitesse et
d
on calcule la distance parcourue par la formule classique v = . On utilise cette formule pour
t
chaque ligne du tableur dans lequel on a exploit la mthode dEuler.
On a ensuite chang lorigine pour prendre le point de dpart du parachutiste 4000 m.
5,000

z + 4000 ()m)

4,000
3,000
2,000
1,000
0

10

20
t(s)

30

40

Figure 7 Position en fonction du temps dans le cas de frottements quadratiques

Quel type de frottements est le plus appropri pour ltude

du mouvement dun parachutiste ?

Il faut faire appel la mcanique des fluides pour rpondre cette question. En effet, on
peut changer de point de vue, et plutt que de considrer la chute du parachutiste dans lair,
on tudie lcoulement de lair autour du parachutiste fixe.
Cest coulement est souvent complexe, il nest pas seulement caractris par la vitesse relative
v du fluide, mais par un nombre sans dimension appel nombre de Reynolds :

10

Mcanique

M12-Chute libre avec frottements

Re =

vd

(28)

9. Quel type de frottements choisir ?

Re : nombre de Reynolds sans dimension.

v : vitesse relative de fluide en m.s .

d : taille caractristique de lcoulement en m.

: masse volumique du fluide en kg.m3 .

: viscosit du fluide en Pa.s.

On distingue alors plusieurs types dcoulement :

si Re < 1, lcoulement est dit laminaire. Dans ce cas, la force de frottements fluides est

proportionnel la vitesse : frottements linaires, f = k

v ;
si Re > 103 , lcoulement est dit turbulent. Alors la force de frottements fluides est

quadratique : f = k v
v.
Dans le cas de la chute du parachutiste, le fluide est lair, sa viscosit est denviron
= 1,7 105 Pa.s ; le nombre de Reynolds a de grande chance dtre suprieur 103 :
lcoulement est turbulent et les frottements quadratiques.

Rfrences
"Physique Tout-en-un MPSI PCSI PTSI" - Marie-Nolle Sanz / Anne-Emmanuelle Badel
/ Franois Clausset - Editions Dunod 2008 ;
Le site Culture Sciences Physiques de lENS de Lyon.
: http ://owl-ge.ch/IMG/pdf/frottement.pdf

11