Resume Cours p202 10nov23

Table des matières
0.1 Plan du cours . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5

0.2 Références . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1 Oscillateur harmonique 7
1.1 Mouvement périodique et mouvement oscillant . . . . . . . . . . . . . . 7
1.2 Système masse-ressort . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.3 Aspects énergétiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.4 Pendule simple . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.5 Circuit électrique oscillant . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.6 Résolution de l’équation di↵érentielle de l’OH . . . . . . . . . . . . . . 13
2 L’oscillateur harmonique amorti, oscillations forcées 16

2.1 Forces de frottement fluide . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.1.1 Mouvement amorti, temps d’amortissement . . . . . . . . . . . 16
2.2 Equation di↵érentielle de OH amorti (OHA) . . . . . . . . . . . . . . . 17
2.3 Résolution de l’équation di↵érentielle de l’OHA . . . . . . . . . . . . . 18
2.3.1 L’équation caractéristique et les trois régimes . . . . . . . . . . 18
2.3.2 Recherche de solutions réelles pour le régime pseudo-périodique 19
2.4 Solutions de l’OHA en notation complexe et les trois régimes . . . . . . 20
2.4.1 Régimes apériodique, critique, pseudo-périodique . . . . . . . . 21
2.5 Aspects énergétiques, facteur de qualité . . . . . . . . . . . . . . . . . . 22
2.6 Oscillations forcées, résonance . . . . . . . . . . . . . . . . . . . . . . . 23
2.6.1 Résonance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.6.2 Fréquence de résonance , largeur de bande . . . . . . . . . . . . 25
2.7 Impédance mécanique, impédance électrique . . . . . . . . . . . . . . . 27
2.7.1 Impédance mécanique . . . . . . . . . . . . . . . . . . . . . . . 27
2.7.2 Impédances électriques, circuit RLC . . . . . . . . . . . . . . . . 28
3 Au delà de l’OH simple 30

3.1 Oscillations anharmoniques . . . . . . . . . . . . . . . . . . . . . . . . 30
3.2 Oscillateurs couplés . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
3.3 Battements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2
3.4 Chaı̂ne infinie d’oscillateurs . . . . . . . . . . . . . . . . . . . . . . . . 33
3.4.1 Relation de dispersion . . . . . . . . . . . . . . . . . . . . . . . 34
3.4.2 Vitesse de groupe, vitesse de phase . . . . . . . . . . . . . . . . 34
4 Cordes vibrantes, équation d’onde de D’Alembert 35

4.1 Les ondes, une approche descriptive . . . . . . . . . . . . . . . . . . . . 35
4.2 Corde vibrante . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
4.3 Equation de d’Alembert et la forme générale des solutions . . . . . . . 38
4.3.1 opérateur D’Alembertien . . . . . . . . . . . . . . . . . . . . . . 39
4.4 Solutions de l’équation de d’Alembert . . . . . . . . . . . . . . . . . . . 39
4.4.1 Solution de l’équation d’onde à 1D . . . . . . . . . . . . . . . . 39
4.4.2 Ondes harmoniques . . . . . . . . . . . . . . . . . . . . . . . . . 40
4.4.3 Détermination des constantes d’intégration à partir des condi-
tions initiales . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
4.5 Impédance mécanique d’une corde . . . . . . . . . . . . . . . . . . . . . 43
5 Réflexion, ondes stationnaires, modes propres d’oscillation 45

5.1 Corde fixée à une extrémité, réflexion . . . . . . . . . . . . . . . . . . . 45
5.2 Corde fixée aux deux extrémités . . . . . . . . . . . . . . . . . . . . . . 46
5.3 Ondes stationnaires . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
5.4 Modes propres . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
5.5 Oscillations forcées, résonance sur une corde tendue . . . . . . . . . . . 50
6 Ondes longitudinales dans les solides 52

6.1 Elasticité des solides . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
6.2 Ondes longitudinales dans les milieux élastiques . . . . . . . . . . . . . 53
6.3 Impédance mécanique . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
6.4 Aspects énergétiques des ondes dans un milieu élastique . . . . . . . . . 56
6.4.1 Densité d’énergie . . . . . . . . . . . . . . . . . . . . . . . . . . 56
6.4.2 Flux d’énergie . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
7 Acoustique, propagation du son dans les fluides 59

7.1 Compressibilité des fluides . . . . . . . . . . . . . . . . . . . . . . . . . 59
7.1.1 Cas du gaz parfait . . . . . . . . . . . . . . . . . . . . . . . . . 59
7.1.2 Autres fluides . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
7.2 Ondes longitudinales dans les fluides . . . . . . . . . . . . . . . . . . . 61
7.2.1 Cas du gaz parfait . . . . . . . . . . . . . . . . . . . . . . . . . 62
7.2.2 Cas général des ondes sonores dans un fluide . . . . . . . . . . . 63
7.3 Impédance acoustique . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
3
8 Intensité sonore, les sources sonores 67
8.1 Aspects énergétiques, flux d’énergie . . . . . . . . . . . . . . . . . . . . 67
8.1.1 Densités d’énergie cinétique et d’énergie potentielle . . . . . . . 67
8.1.2 Flux d’énergie et intensité acoustique . . . . . . . . . . . . . . . 69
8.1.3 Intensité sonore moyenne . . . . . . . . . . . . . . . . . . . . . . 70
8.2 Niveau acoustique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
8.3 sources sonores . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
8.3.1 Instruments à cordes . . . . . . . . . . . . . . . . . . . . . . . . 73
8.3.2 Instruments à vent, tuyaux sonores . . . . . . . . . . . . . . . . 74
9 Quelques phénomènes ondulatoires 77

9.1 Propagation des ondes à 3 dimensions . . . . . . . . . . . . . . . . . . . 77
9.1.1 Equation des ondes à 3D . . . . . . . . . . . . . . . . . . . . . . 77
9.1.2 Ondes planes . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
9.1.3 Ondes sphériques . . . . . . . . . . . . . . . . . . . . . . . . . . 78
9.2 E↵et Doppler . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
9.2.1 Notes sur l’e↵et Doppler . . . . . . . . . . . . . . . . . . . . . . 80
9.3 Onde de choc . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
10 Réflexion et transmission 84
10.1 Onde à l’interface entre deux milieux . . . . . . . . . . . . . . . . . . . 84
10.1.1 Réflexion et transmission du flux d’énergie . . . . . . . . . . . . 86
11 Interférences, Di↵raction 88
11.1 Description ondulatoire de la lumière . . . . . . . . . . . . . . . . . . . 88
11.2 Interférences . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
11.2.1 Superposition de deux ondes, e↵et de la di↵érence de phase . . . 89
11.2.2 Expérience des trous d’Young . . . . . . . . . . . . . . . . . . . 90
11.2.3 Interféromètre de Michelson . . . . . . . . . . . . . . . . . . . . 91
11.3 Di↵raction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
4
Plan, références
Module de 46 heures, 11 cours (x 2 h = 22 h) et 12 TD (x 2 = 24 h) Plan
0.1 Plan du cours

1. Oscillateur harmonique : équation di↵érentielle, circuit LC
2. OH amorti, oscillations forcées et résonance
3. Corde vibrante . Equation d’Onde de l’Alembert et ses solutions
4. Réflexion des ondes, ondes stationnaires
5. Modes propres d’oscillations (corde fixée à une extrémité)
6. Acoustique, ondes dans les fluides
7. Aspects énergétiques pour les ondes acoustiques, impédance mécanique
8. Sensibilité de l’oreille, cordes vibrantes et tuyaux acoustiques
9. Ondes planes et ondes sphériques
10. E↵et Doppler, Interférences
Non traitées en Phys202 :
1. Au delà de l’OH simple, couplage de deux oscillateurs, Chaı̂ne d’oscillateur
2. Ondes longitudinales dans les solides, aspects énergétiques
3. Réflexion et transmission des ondes aux interfaces
4. Di↵raction des ondes
5
0.2 Références
1. Cours de Physique de Berkeley physics, Volume 1, Mécanique, C. Kittel, W.D.
Knight, M.A. Ruderman, traduction P. Lallemand, Ed. Armand Colin / Dunod
2. Cours de Physique de Berkeley physics, Volume 3, Ondes, Franck S. Crawford,
Jr traduction P. Léna, Ed. Armand Colin / Dunod
3. Toute la Thermodynamique, la Mécanique des fluides et les ondes mécaniques,
L. Bocquet, J.P. Faroux, J. Renault Ed. Dunod
4. Physique 3. Ondes, optique et physique moderne Resnick, Halliday traduction
Nguyen, Martel , Ed. du renouveau pédagogique
5. Physique 3 Ondes, Optique et physique moderne, René Lafrance Chenelière
Éducation
6. Vibrations et ondes Tamer Bécherrawy , Ed. Hermes Science / Lavoisier
7. Vibrations, ondes et optique - Ondes mécaniques (Volume 2) Tamer Bécherrawy
, Ed. Hermes Science / Lavoisier
6
Chapitre 1
Oscillateur harmonique
1.1 Mouvement périodique et mouvement oscillant

- Une fonction f (t) de R ! R est périodique s’il existe un réel T 2 R⇤ , T > 0 tel
que
f (t + T ) = f (t) 8t
Notons que si T existe, alors la relation ci-dessus est également vérifiée pour
toutes les valeurs Tk = kT avec k 2 Z. On appelle T la période de la fonction f , la
plus petite valeur T strictement positive pour laquelle la relation f (t+T ) = f (t)
est vérifiée pour tout t. Cette définition se généralise aisément au cas de fonctions
non scalaire d’une variable réelle, à valeur dans R ou pas (complexe, entier . . .).
- On appelle un mouvement périodique l’état dynamique d’un système physique
(mécanique, électrique, . . .) dont l’état, la position pour un point matériel par
exemple, est décrit par une fonction périodique. Un mouvement périodique se
répète donc identique à lui même, avec une périodicité T . L’inverse de la période
⌫ = T1 est appelée la fréquence et correspond au nombre de fois où le mouvement
se répète en une unité de temps. La période a la dimension d’une temps, avec
la seconde comme unité SI, et la fréquence a comme dimension l’inverse d’une
temps, et le Hertz comme unité SI. Une fréquence peut bien sûr être inférieure
à l’unité si T > 1.
- Les mouvements périodiques correspondent souvent à une mouvement autour
d’une position d’équilibre. Une oscillation est un mouvement de va et vient
autour de cette position d’équilibre.
- Les mouvements d’un pendule, d’une masse suspendue à un ressort, ou d’une
bille au fond d’un bol sont des exemples d’oscillations périodiques. Notons que
les mouvements des planètes autour du soleil ou de la lune autour de la terre
sont d’autres exemples de mouvement périodique.
Ce chapitre est consacré à l’étude de oscillations périodiques de système mécanique ou
électrique. Notons que les systèmes oscillants, avec un mouvement périodique ont été
7
et sont à la base de toutes les horloges, et même de la mesure de l’écoulement du temps
au sens large.
1.2 Système masse-ressort

On étudie le mouvement à une dimension d’un système composé d’une masse m
attachée à l’une des extrémités d’un ressort de raideur k, l’autre extrémité étant fixe
(attaché à un mur par exemple). On note x(t) la position de la masse en fonction
du temps et on suppose que la position d’équilibre F = 0 correspond à x = 0. Le
ressort exerce une force de rappel F = k x u~x sur la masse m, proportionnelle à son
allongement.
Notons qu’un ressort a une longueur au repos, en absence de contrainte `0 . La force
exercée par le ressort est donnée par la Loi de Hooke 1 : qui exprime la proportionnalité
entre la déformation du ressort selon son axe ` = ` `0 , allongement ` > 0 ou
compression ` < 0, et la force exercée. Il faut exercer une force de traction pour
allonger le ressort et une force de compression pour diminuer sa longueur. La loi de
Hooke s’écrit :
F= k `= k (` `0 )
où k est la raideur du ressort ou la constante de rappel. Le signe moins exprime le fait
que le ressort s’oppose à la déformation. Il faut néanmoins être attentif lors de l’écriture
des composantes du vecteur force. Noter par exemple que les forces s’exerçant aux deux
extrémités du ressort sont de sens opposé.
Revenons au système masse-ressort. On choisit l’origine de l’axe Ox afin que la
position x = 0 coincide avec le ressort au repos, x = 0 correspond à ` = `0 . On applique
le principe fondamental de la dynamique et on obtient l’équation di↵érentielle satisfaite
par la position de la masse m.
k
mẍ = kx ! ẍ + x=0
mr
k
ẍ + !02 x = 0 avec !0 =
m
2⇡ 1 !0
!0 = ⌫0 = =
T0 T0 2⇡
Il s’agit d’une équation di↵érentielle linéaire de second ordre, dont les solutions s’écrivent
sous forme de somme de fonctions sinus et cosinus (fonctions harmoniques). Il y a deux
1. Robert Hooke, scientifique anglais du XVII siècle, 1635-1703
8
constantes d’intégration qui doivent être déterminées à partir des conditions initiales.
x(t) = C cos(!0 t) + S sin(!0 t)

ẋ(t) = v(t) = !0 ( C sin(!0 t) + S cos(!0 t))
ẍ(t) = a(t) = !02 (C cos(!0 t) + S sin(!0 t)) = !02 x(t)
La solution peut s’écrire de manière équivalente, à l’aide d’une seule fonction sinus (ou
cosinus), une amplitude A et une phase .
x(t) = A cos(!0 t + ')

ẋ(t) = !0 A sin(!0 t + ')
x(0) = A cos(') ẋ(0) = A!0 sin(')
x(t) = A sin(!0 t + 's )
ẋ(t) = !0 A cos(!0 t + 's )
x(0) = A sin('s ) ẋ(0) = A!0 cos('s )
⇡ ẋ(0)
's = ' = !0 tan(') = !0 cot('s )
2 x(0)
Expression de A, ' en fonction de C, S.
A cos(!0 t + ') = A (cos ' cos(!0 t) sin ' sin(!0 t))

= C cos(!0 t) + S sin(!0 t)
C = A cos ' S = A sin '
S
A2 = C 2 + S 2 tan ' =
C
p
Lepmouvement s’e↵ectue avec la pulsation !0 = k/m, donc une période T0 =
2⇡ m/k.
r
k
Pulsation : !0 = 2⇡⌫0 =
m r
1 2⇡ m
Periode : T0 = = = 2⇡
⌫0 !0 k
Celle-ci (T) ) est indépendante de l’amplitude du mouvement A. A noter le lien
entre un mouvement circulaire uniforme dans un plan et le mouvement oscillatoire. En
utilisant les coordonnées polaires (r, ) :
r(t) = Cte = R (t) = !0 t + '0

x(t) = R cos( (t)) y(t) = R sin( (t))
x(t) = R cos(!t + '0 )
9
1.3 Aspects énergétiques
La force F (x) = kx peut s’écrire comme le gradient de l’énergie potentielle U (x) :
1 2
U (x) = kx
2
dU
F (x) = kx = F= gradU = rU
dx
L’énergie mécanique totale, égale à la somme de l’énergie cinétique et de l’énergie
potentielle est conservée au cous du mouvement. On a une conversion de l’énergie
potentielle en énergie cinétique et vice versa. On note K(t) l’énergie cinétique de la
masse m au cours du mouvement, U (t) son énergie potentielle, et E(t) = K(t) + U (t)
l’énergie mécanique totale. En notant que k = m!02 , on montre facilement que l’énergie
mécanique totale est conservée, égale à sa valeur E0 à l’instant intial E0 = E(t = 0) :
1 1
K(t) = mẋ2 = (A!0 sin(!0 t + '))2
2 2
1 2 1
U (t) = kx = k (A cos(!0 t + '))2
2 2
1 ⇥ ⇤
E(t) = K(t) + U (t) = A2 m!02 cos2 (!0 t + ) + k cos2 (!0 t + ')
2
1 2
E(t) = kA = E0 = Cte
2
On peut aussi calculer la moyenne de l’énergie cinétique et de l’énergie potentielle en
notant que la moyenne de sin2 (!0 t) ou cos2 (!0 t) sur une période est égale à 1/2 :
Z
1 T 1 1
hKi = K(t)dt = m!02 A2 = kA2
T 0 4 4
Z T
1 1
hU i = K(t)dt = kA2 = hKi
T 0 4
E0
hKi = hU i =
2
1.4 Pendule simple

Le pendule simple est constitué d’une masse ponctuelle m, fixée à un pivot par une
tige rigide de longueur ` et de masse négligeable. Le pendule peut se déplacer dans un
10
Figure 1.1 – Illustration de la conservation de l’énergie mécanique totale pour un os-
cillateur harmonique. La courbe rouge représente l’énergie potentielle U (x) en fonction
de la position de la masse, la ligne noire horizontale la valeur de l’énergie totale E0 et
la flèche verticale correspond à l’énergie cinétique K(x) = E0 U (x).
plan, la trajectoire de la masse étant un arc de centre de rayon ` et centré sur le pivot.
La masse m est soumise à l’action de la pesanteur, ou le poids (P = mg u~z ), à la force
de réaction de la tige, alignée avec celle-ci. A noter que la tige peut être remplacée par
un fil non extensible. On note ✓ la position angulaire du pendule, ✓ = 0 représentant
la position d’équilibre stable du pendule, lorsque la tige est alignée avec l’axe vertical
vecuz et la masse vers le bas. On projète les forces et l’accélération selon une direction
tangente (k, ~u✓ ) à la trajectoire (arc de cercle) et selon une direction perpendiculaire
(?, ~ur ) à la trajectoire. L’application du principe fondamental de la dynamique selon la
direction tangente, où seule intervient la projection du poids permet d’écrire l’équation
di↵érentielle satisfaite par l’angle ✓ :
d2 ✓
ak = ` = `✓¨
dt2
Pk = mg sin ✓~u✓
m`✓¨ = mg sin ✓
g
✓¨ + sin ✓ = 0
`
On obtient une équation di↵érentielle non linéaire d’ordre un, pas facile à résoudre.
Mais si on se limite à des petits écarts par rapport à la position d’équilibre ✓ ⌧ 1, et
11
qu’on remplace sin ✓ par son développement limité à l’ordre 1, on obtient :
sin ✓ ' ✓ pour ✓ ⌧ 1

g
✓¨ + sin ✓ = 0
`
g
✓¨ + ✓ = 0 pour ✓ ⌧ 1
` r
g
✓¨ + !02 ✓ = 0 !0 =
`
On retrouve une équation di↵érentielle d’oscillateur harmonique.qLes mouvements de

faible amplitude du pendule se font avec une période T = 2⇡ g` indépendante de
l’amplitude des oscillations.
1.5 Circuit électrique oscillant

On considère un circuit électrique formé d’un condensateur de capacité C et d’une
bobine avec une résistance négligeable et une auto-inductance L, connectés en série. La
di↵érence de potentielle aux bornes du condensateur s’écrit VC = Cq où q(t) est la charge
électrique stockée dans le condensateur à l’instant t. Le courant électrique i(t) = dqdt
dans le circuit correspond à la variation de la charge aux bornes du condensateur et la
di
di↵érence de potentiel aux bornes de la bobine s’écrit VL = L dt . En appliquent la loi
des maille VC + VL = 0, on obtient l’équation di↵érentielle satisfaire par la charge aux
bornes du condensateur :
1 dq
VC (t) =q(t) i(t) =
C dt
di d2 q
VL (t) = L = L 2
dt dt
d2 q 1
VL (t) + VC (t) = 0 ! 2
+ q(t) = 0
dt rLC
d2 q 1
2
+ !02 q(t) = 0 !0 =
dt LC
On retrouve de nouveau une équation di↵érentielle de type Oscillateur harmonique.

La position d’équilibre mécanique correspond à q(t) = 0 et i(t) = 0 8t. Sinon, la charge
aux bornes du condensateur, de même que le courant dans le circuit et les tensions aux
bornes du condensateur et de la bobine auront un comportement oscillatoire périodique.
On peut établir une correspondance entre les quantités mécaniques et électrique :
12
Position : x(t) ! Charge : q(t)
Vitesse : ẋ ! Courant : i(t)
1
Force : F = kx ! Potentiel electrique : VC (t) = q(t)
C
1
Raideur ressort : k ! Capacite :
C
masse : m ! Inductance : L
1 1 2
Energie potentielle : U = kx2 ! Uelectrostatique = q
2 2C
1 1
Energie cinetique : K = mẋ2 ! Emagnetique = L i2
2 2
1.6 Résolution de l’équation di↵érentielle de l’OH

L’équation di↵érentielle de l’OH est une équation di↵érentielle satisfaite par une
fonction réelle x(t) de la forme :
d2 x
+ !02 x(t) = 0 !0 2 R; !02 > 0
dt2
Il s’agit d’une équation di↵érentielle linéaire d’ordre 2 sans second membre. On rappelle
que toute combinaison linéaire de solutions d’une équation di↵érentielle linéaire est
également une solution de l’équation di↵érentielle.
En présence d’un second membre, par exemple lorsqu’une force extérieure addition-
nelle Fext (t) s’applique aussi sur la masse attachée à un ressort :
mẍ = kx + Fext (t)

2
dx Fext (t)
2
+ !02 x(t) = fext (t) avec fext (t) =
dt m
la solution générale de l’équation avec second membre s’écrit comme la somme de

l’équation générale de l’équation di↵érentielle sans second membre et d’une solution
particulière avec second membre.
On va chercher ici les solutions d’une équation di↵érentielle linéaire de second ordre
de la forme :
ẍ + ax = 0
13
l’ED de l’OH ci-dessus correspond au cas où x(t), a 2 R et (a = !02 ) > 0. On multiplie
l’équation par ẋ :
ẍẋ + axẋ = 0
d d 2
(ẋ)2 + a x = 0
dt dt
d 2
ẋ + ax2 = 0
dt
On intègre et on obtient l’expression satisfaite par ẋ en fonction d’une première constante

d’intégration notée B. On note que B > 0 pour a > 0, sauf si x(t) = 0 où B = 0.
ẋ2 + ax2 = B
p
ẋ = B ax2
Le cas B = 0 est possible si a < 0. On obtient alors la solution x(t) :

Cas B = 0 avec a < 0
La solution x(t) s’écrit sous forme de la somme de deux fonctions exponentielles, l’une
croissante et l’autre décroissante, et les deux constantes d’intégration C+ , C :
p
ẋ2 = ax2 ! ẋ = ±b ⇥ x(t) b2 = a b= a
x(t) = C+ ebt + C e bt
Cas B > 0 :
p
Pourple cas général où B > 0, on peut écrire l’équation
p ẋ = B ax2 sous la forme
ẏ ⇠ 1 ± y 2 , en posant y(t) = ↵x(t) avec ↵ = ±a/B selon que a > 0 ou a < 0
p
ẋ = B ax2 B > 0 et B ax2 > 0
r r
p a 2 a
ẋ = B 1 x On pose : y(t) = ± x(t)
B B
B p
ẏ = p 1 ⌥ y2
±a
Si on connaı̂t la solution générale z(t) de l’équation l’équation di↵érentielle ż =

H(z), alors y(t) = z( t) satisfait
p l’équation di↵érentielle ẏ = H(y). On va donc
essayer de résoudre l’ED ż = 1 ⌥ z 2 , ce qui se fait assez facilement en se rappelant
14
p
que l’intégrale de 1/ 1 ± z 2 , correspond à arcsin et arcsinh :
dz p
ż = = 1 z2 cas a > 0
dt
Z z(t) Z t
dz
p = dt
z0 1 z2 0
arcsin(z(t)) arcsin(z0 ) = t
z(t) = sin(t + arcsin(z0 ) = sin(t + t0 )
dz p
ż = = 1 + z2 cas a > 0
dt
z(t) = sinh(t + asinh(z0 )) = sinh(t + t0 )
Pour obtenir y(t), il suffit de remplacer t par t dans les expressions ci-dessus, avec
= pB±a . La solution x(t) de l’équation di↵érentielle initiale s’obtient en se rappelant
p
que y(t) = ±a/Bx(t). Or l’ED de départ étant linéaire, x(t) = Cy(t) en serait une
solution pour 8(t).
La solution générale pour le cas a < 0 est la même que celle établie ci-dessus pour
B = 0, et une somme de fonction sinus et cosinus pour a > 0 :
p p
a<0 : x(t) = C+ exp at + C exp at
p p
a>0 : x(t) = C cos at + S sin at
15
Chapitre 2
L’oscillateur harmonique amorti,

oscillations forcées
2.1 Forces de frottement fluide

Les forces de frottements sont en général associées aux phénomènes de dissipation,
où l’énergie est dégradée en chaleur. Les forces de frottement solide en mécanique sont
en général représentées par un force d’intensité constante, ayant la même direction que
le déplacement, donc la vitesse, mais de sens opposé. Les forces de frottement solide
sont dues à l’interaction entre l’objet qui se déplace et la surface sur laquelle il se
déplace. Elles ont parallèles à la surface et ont une intensité proportionnelle à la force
de réaction du support, normale à celui-ci.
Les forces de frottement fluide regroupe des forces de frottement dont l’intensité
est proportionnelle à la vitesse de déplacement, et bien sûr, avec un sens opposé au
déplacement.
d~r
Ff f /
dt
Ce type de force apparaı̂t en e↵et lors de mouvement de faible vitesse d’un objet dans
un fluide, liquide ou gaz. Dans les gaz, elles peuvent souvent être masquées par les
forces de traı̂née aérodynamique, qui sont elles proportionnelles au carré de la vitesse.
Une résistance dans un circuit électrique a un e↵et analogue à des forces de frottement
fluide, provoquant une chute du potentiel électrique VR = Ri(t) proportionnelle au
courant électrique qui la traverse.
2.1.1 Mouvement amorti, temps d’amortissement

On considère dans un premier temps une masse m soumise uniquement à une force
de frottement fluide aura un mouvement de type exponentielle amortie. On considère
un mouvement à une dimension, selon l’axe Ox et une force Ff f = ⌘ ẋ où ⌘ est le
16
coefficient de frottement fluide, ayant la dimension d’une force divisée par une vitesse,
ce qui correspond à une masse/temps (unité SI kgs 1 ).
⌘ 1 ⌘
mẍ = ⌘ ẋ ! ẍ + ẋ = 0 ⌧
=m
m
On écrit l’équation di↵érentielle satisfaite par la vitesse v(t) = ẋ en fonction de la
constante ⌧ = m ⌘
qui représente le temps caractéristique du mouvement. Il s’agit d’une
équation di↵érentielle linéaire de premier ordre qui se résout facilement, introduisant
une première constante d’intégration v0 qui correspond à la vitesse à t = 0.
1
v̇ + v = 0
⌧ ✓ ◆
t
v(t) = ẋ = v0 exp
⌧
L’expression de l’équation du mouvement x(t) s’obtient en intégrant l’expression de
v(t) et qui amène l’introduction d’une seconde constante d’intégration, x0 = x(t = 0),
la position à t = 0 On peut aussi calculer l’expression de l’énergie cinétique qui décroit
aussi exponentiellement, avec un temps caractéristique deux fois plus court ⌧ /2 :
✓ ✓ ◆◆
t
x(t) = x0 + v0 ⌧ 1 exp
⌧
✓ ◆
1 2 2t
K(t) = mv0 exp
2 ⌧
⌘
Noter que ⌧ avec ⌧1 = m a la dimension d’un temps et correspond au temps ca-
ractéristique d’amortissement de la vitesse. Au bout du temps ⌧ , la vitesse est réduite
d’un facteur e ' 2.7 par rapport à la vitesse initiale, et d’un facteur e2 ' 7.4 à t = 2⌧ .
L’énergie cinétique correspondant au carré de la vitesse, le temps d’amortissement
(relaxation) de l’énergie cinétique est donc la moitié de celle de la vitesse = ⌧2 .
2.2 Equation di↵érentielle de OH amorti (OHA)

On ajoute un dispositif exerçant une force de frottement fluide sur le système masse-
ressort étudié au chapitre 1, un appendice plongeant dans un bain d’huile par exemple.
L’équation di↵érentielle satisfaite par la position x(t) de la masse s’écrit sous la forme
q
m k
suivante en introduisant les deux constantes ⌧ = ⌘
homogène à un temps et !0 = m
,
homogène à l’inverse d’un temps :
mẍ = Fressort + Ff f = kx ⌘ ẋ
⌘ k
ẍ + ẋ + x = 0
m m
1 1 ⌘ k
ẍ + ẋ + !02 x = 0 = ; !02 =
⌧ ⌧ m m
17
L’équation di↵érentielle satisfaite par la variable dynamique d’autres systèmes de
type OH peut également se mettre sous une forme analogue. Si on considère par exemple
un circuit électrique avec un condensateur et une bobine connectée en série, auquel on
ajoute une résistance, on obtient :
di 1 dq
L + Ri + q = 0 i(t) = = q̇
dt C dt
R 1
q̈ + q̇ + q = 0
L LC
1 1 R 1
q̈ + q̇ + !02 q = 0 = ; !02 =
⌧ ⌧ L LC
Comme on va le voir dans la suite, les solutions correspondent le plus souvent à des os-
cillations amorties et cette équation di↵érentielle est appelée l’équation de l’Oscillateur
Harmonique Amorti (OHA).
Cette équation est souvent écrite sous une forme légèrement di↵érente en introdui-
sant un nouveau coefficient sans dimension noté ⇠, caractérisant l’intensité de la force
de frottement fluide. Rappelons que ⌧ et !10 ont la même dimension, celle d’un temps. ⌧
est le temps d’amortissement, et T0 = !2⇡0 est la période d’oscillation de l’OH simple. ⇠
s’écrit comme 4⇡ fois le rapport entre la période d’oscillation T0 = 2⇡/!0 , en absence
d’amortissement, et le temps caractéristique d’amortissement ⌧ , défini ci-dessus, dans
la section 2.1.1 et correspond à l’intensité des forces de frottement fluide, rapportée à
celle de la force de rappel.
1 1 T0 1
⇠= = ! ⌧=
2⌧ !0 4⇡ ⌧ 2⇠!0
1 1 ⌘
ẍ + ẋ + !02 x = 0 =
⌧ ⌧ m
1 ⌘
ẍ + 2⇠!0 ẋ + !02 x = 0 ⇠= p
2 km
2.3 Résolution de l’équation di↵érentielle de l’OHA

2.3.1 L’équation caractéristique et les trois régimes
On recherche des solutions sous la forme x(t) = Ae t . On calcule la dérivée première
et la dérivée seconde et on remplace dans l’équation di↵érentielle et on voit que doit
satisfaire une équation algébrique de second degré, appelée équation caractéristique.
x(t) = Ae t ! ẋ(t) = A e t ẍ(t) = A 2 e t
2 1
+ + !02 = 2
+ 2⇠!0 + !02 = 0
⌧
18
Le déterminant s’écrit = 4!02 (⇠ 2 1). Il y a deux racines réelles, négatives toutes
les deux, si > 0, ce qui correspond à un coefficient d’amortissement ⇠ > 1, supérieur
à l’unité. Le système a un mouvement apériodique dans ce régime avec un fort amor-
tissement. On a une seule racine réelle double dans cas limite = 0, ⇠ = 1, appelé
régime critique. Le régime pseudo-périodique correspond au cas < 0, ⇠ < 1, où il
y a deux racines complexes, avec une partie réelle négative et une partie imaginaire.
⇣ p ⌘
⇠>1! >0 : = ! ⇠ ± ⇠ 2 1
12 0
⇠=1! =0 : 1 = 2 = ⇠!0
⇣ p ⌘
⇠<1! <0 : 12 = ⇠!0 ± i!p = !0 ⇠±i 1 ⇠2
La solution de l’ED dans le régime pseudo-périodique correspond à des oscillations

amorties, produit d’une exponentielle décroissante et d’une fonction sinusoı̈dale, comme
discuté ci-dessous. !p = !0 (1 ⇠ 2 ) est la pulsation de ces oscillations amorties.
2.3.2 Recherche de solutions réelles pour le régime pseudo-

périodique
On peut aussi adopter une approche plus physique et rechercher directement des
solutions de la forme x(t) = A exp( t) sin(!p t+') qui correspondent à des oscillations
amorties. On va déterminer les conditions qui doivent être satisfaites par !p et en
calculant ẋ et ẍ et en injectant les expressions dans l’équation di↵érentielle de l’OHA.
On verra que !p et seront indépendant de la phase , et de l’amplitude qui devront
être déterminées par les conditions initiales.
✓ ◆
2
!p2 + !02 Ae t sin(!p t + ') +
⌧
⇣ !p ⌘
2!p + Ae t cos(!p t + ') = 0
⌧
Les deux expressions algébriques contenant !p et devant les termes sin(!p t) et

cos(!p t) doivent s’annuler simultanément. On obtient alors :
✓ ◆
1 1
!p 2 =0 ! =
⌧ 2⌧
✓ ◆2
1
!02 + 2
= !p2 ! !p2 = !02 = !02 1 ⇠2
⌧ 2⌧
19
Notons qu’une solution réelle pour !p n’existe que si !0 > 2⌧1 . Dans ce cas, la solution
s’écrit sous la forme suivante, ou bien à l’aide du coefficient sans dimension ⇠ avec
1
2⌧
= ⇠!0 :
✓ ◆ ✓ ◆2
t 1
x(t) = A0 exp sin(!p t + '0 ) !p2 = !02
2⌧ 2⌧
= ⇠!0 !p2 = !02 1 ⇠2
⇣ p ⌘
x(t) = A0 exp ( ⇠!0 t) sin !0 1 2
⇠ t + '0 pour ⇠ < 1
La solution que nous venons de trouver, sous forme d’oscillations amorties, n’existe
que si ⇠ < 1 (ou 2⌧1 < !0 ). C’est le régime pseudo-périodique de l’OHA, avec un
mouvement sinusoı̈dal de pulsation !p , dont l’amplitude décroit avec le temps. Notons
que l’amortissement abaisse la fréquence des oscillations par rapport à l’OH non amorti
(!p < !0 ). Dans le cas d’un amortissement faible, pour ⇠ ⌧ 1 (ou !0 ⌧ 1 et T0 ⌧ ⌧ ),
on peut faire un DL à l’ordre 1 en ⇠ et on obtient l’expression approchée de !p ' !0
et à l’ordre 2 :
✓ ◆
2 2 2 1 2
!p = !0 (1 ⇠ ) ! !p ' !0 1 ⇠ !p < !0
2
Voir ci-dessous, la section 2.4.1 pour une discussion plus détaillée des trois régimes
pour les solutions de l’OHA.
2.4 Solutions de l’OHA en notation complexe et les

trois régimes
On se place en notation complexe et on définit x(t) = Re(x̃(t)) avec x̃(t) 2 C.
On recherche donc des solutions complexes de l’ED sous la forme x̃(t) = X̃0 e t , avec
X̃0 , 2 C. Le calcul de la dérivée par rapport au temps s’e↵ectue aisément, l’application
de l’opérateur d/dt se réduit à une multiplication par :
d t d2 2 t
x̃(t) = X̃0 e x̃(t) = X̃0 e
dt dt2
On remplace les expressions de dx̃/dt et d2 x̃/dt2 dans l’équation di↵érentielle de l’OHA
et on voit que doit satisfaire une équation algébrique de second degré, appelée
l’équation caractéristique.
x(t) = Re(x̃(t)) x̃(t) = X̃0 e t

2C
2
+ 2⇠!0 + !02 X̃0 e t = 0
! 2 + 2⇠!0 + !02 = 0 (Eq. caracteristique)
20
On peut noter que cette équation algébrique de 2nd degré admet toujours deux
racines complexes et que les conditions initiales déterminent X̃0 , cad l’amplitude et la
phase.
= 4!02 ⇠ 2 1 ⇠<1! <0

⇣ p ⌘
1,2 = !0 ⇠ ± ⇠2 1
Lorsque ⇠ > 1, on a deux racines réelles, et une racine double 1 = 2 = ⇠!0 pour
le cas ⇠ = 1. Pour ces deux régimes, apériodique, respectivement critique, les racines
sont réelles et négatives, les solutions sont donc des exponentielles décroissantes :
⇣ ⇣ p ⌘ ⌘
x(t) = X̃0 exp !0 ⇠ ± ⇠ 2 1 t x(t) = X̃0 exp ( !0 t)
Finalement, dans le cas ⇠ < 1, le déterminant est négatif, et on a deux racines com-
plexes, avec une partie réelle négative et une partie imaginaire.
⇣ p ⌘
1,2 = !0 ⇠ ± i 1 ⇠2
⇣ p ⌘
x(t) = X̃0 exp( xi!0 t) exp ±i!0 1 ⇠ 2 t
h ⇣ p ⌘ ⇣ p ⌘i
2
x(t) = X̃0 exp( xi!0 t) cos !0 1 ⇠ t ± i sin !0 1 ⇠ t 2
2.4.1 Régimes apériodique, critique, pseudo-périodique

On résume ci-dessous l’expression des solutions réelles pour les trois régimes.
• Régime apériodique , ⇠ > 1 ! > 0. L’équation caractéristique admet deux
racines réelles négatives ( 1,2 < 0). La solution générale s’écrit comme la somme
de deux exponentielles décroissante, avec deux amplitudes A, B dont les valeurs
dépendent des conditions initiales.
x(t) = A exp( 1 t) + B exp( 2 t)

⇣ p ⌘
1,2 = !0 ⇠ ± ⇠2 1 1,2 2R 1,2 <0
• Le régime critique correspond au cas limite où ⇠ = 1 ! = 0 cad = 2⌧1 =

!0 . L’équation caractéristique admet une seule racine réelle , = !0 . On
vérifie aisément t exp( !0 t) est aussi solution de l’ED, en plus de exp( !0 t).
21
Figure 2.1 – Mouvement d’un oscillateur harmonique amorti en régime pseudo-
périodique pour trois valeurs du coefficient d’amortissement.
La solution générale de l’équation di↵érentielle s’écrit donc comme le produit

d’une exponentielle amortie et d’une fonction affine du temps :
x(t) = exp( !0 t) [A + Bt] (⇠ = 1 ! ⇠!0 = !0 )
Le régime critique est celui qui assure le retour le plus rapide à la position
d’équilibre, quelque soit les conditions initiales.
• Régime pseudo-périodique , ⇠ < 1 ! < 0. L’équation caractéristique admet
deux racines complexes, avec une partie réelle négative et une partie imaginaire
pure, comme nous l’avons vu ci-dessus. La partie imaginaire de 1,2 se traduit
alors par des fonctions harmoniques cos !t et sin !t :
x(t) = exp( ⇠!0 t) [C cos(!p t) + S sin(!p t)]

p
x(t) = A exp( ⇠!0 t) cos(!p t + ') !p = !0 1 ⇠2
Le pseudo-période Tp = !2⇡p est égale à l’intervalle de temps séparant deux

passages par la position d’équilibre, dans le même sens de déplacement. La figure
2.1 représente l’allure du mouvement pour un oscillateur harmonique amorti en
régime pseudo-périodique.
2.5 Aspects énergétiques, facteur de qualité

1. Aspects énergétiques, facteur de qualité. On se contente de se placer dans le cas
des amortissement faibles, où on néglige la variation d’amplitude t au cours d’une
22
période d’oscillation. On calcule ainsi la puissance moyenne dissipée P (t) au
cours de chaque cycle. Ici, par la moyenne, il faut comprendre la moyenne sur un
cycle d’oscillations. L’amplitude des oscillations A(t) décroı̂t exponentiellement
avec le temps
✓ ◆ ✓ ◆
t 2 2 t
A(t) = A0 exp ! A(t) = A0 exp
2⌧ ⌧
✓ ◆
1 1 t 1
hKi(t) = m!02 A(t)2 = kA20 exp = kA20 exp ( 2⇠!0 t)
4 4 ⌧ 4
1
hU i(t) = kA(t)2 = hKi(t)
4
1 1
hEi(t) = kA(t)2 = kA20 exp ( 2⇠!0 t)
2 2
d hEi
hP i(t) = hEi = 2⇠!0 hEi =
dt ⌧
2. Facteur de qualité On définit le facteur de qualité Q comme le rapport entre

l’énergie emmagasinée (stockée) et l’énergie disspiée durant un cycle d’oscilla-
tions (période Tp ), multipliée par 2⇡
hEi hEi
Q = 2⇡ = !p
hP iT hP i
!p
Q = !p ⌧ =
2⇠!0
1
Q ' pour ⇠ ⌧ 1 ! ! ' !0 Q 1
2⇠
Valeurs typiques (d’après cours de physique de Berkeley, mécanique, Chap 7) -
lien avec la largeur en fréquence (Q élevé, faible largeur en fréquence)
Terre (tremblement de terre) ⇠ 102 103
Corde de piano ⇠ 103
Cavité micro-ondes cuivre 104
Atome excité ⇠ 107
Noyau excité (57 Fe) ⇠ 1012
2.6 Oscillations forcées, résonance

2.6.1 Résonance
1. On considère un OHA soumis à une force extérieure harmonique F (t) = F0 cos(!t),
ou dans le cas d’un circuit RLC, où un générateur produisant une di↵érence de
23
potentiel sinusoı̈dal V (t) = V0 cos(!t) est ajouté dans le circuit. On recherche
la solution générale de l’équation di↵érentielle :
r
k 1 ⌘
mẍ + ⌘ ẋ + kx = F0 cos(!t) !02 = ; ⇠= p
m 2 km
F0
ẍ + 2⇠!0 ẋ + !02 x = f0 cos(!t) f0 =
m
Ou, en utilisant la première notation introduite avec le temps d’amortissement

⌧ :
1 1 ⌘ F0
ẍ + ẋ + !02 x = f0 cos(!t) = ; f0 =
⌧ ⌧ m m
2. Solution de l’ED, somme de la solution générale de l’ED sans second membre,
plus une solution particulière xp (t) de l’ED avec second membre. On utilise la
notation complexe pour rechercher x̃p (t) :
x̃p (t) = X̃0 ei!t

f˜0
X̃0 =
(!02 ! 2 ) + 2i⇠!0 !
Ou bien en utilisant la première notation utilisant ⌧ :
f˜0
X̃0 =
(!02 ! 2 ) + i!/⌧
X̃0 1 1
= !a =
f˜0 (!02 2
! ) + i!!a ⌧
3. Il faut bien sûr ajouter à cette solution particulière, la solution générale xg (t)
de l’OHA et dont les constantes d’intégration, l’amplitude initiale et la phase
dépendent de l’état du système à t = 0, position et vitesse. Or nous avons
vu que xg (t) revient vers zéro, quelque soit le régime, avec une exponentielle
décroissante, ⇠ exp( ⇠!0 t). Le système perd la mémoire de ses conditions ini-
tiales, au bout d’une temps t = n⌧ avec n de quelques unités. On s’intéresse
donc uniquement dans la suite à la solution particulière xp (t), qui correspond
au mouvement en régime permanent, après le régime transitoire où s’e↵ace
petit à petit l’influence des conditions initiales.
4. Le mouvement en régime permanent pour une excitation harmonique est est
un mouvement harmonique avec la même fréquence (ou la période) que celle
de l’excitation, avec une réponse du système caractérisée par le gain g(!) et le
24
déphasage (!). Le gain est le rapport entre l’intensité de l’excitation et l’am-
plitude du mouvement induit ou de la réponse, alors que le déphasage représente
la di↵érence de phase entre l’excitation et le mouvement induit. La figure 2.2
représente la variation du gain g( !!0 ) pour quelques valeurs du coefficient d’amor-
tissement ⇠.
|X̃0 | 1
Gain : g(!) = =q
|f˜0 | 2
(!02 ! 2 ) + (2⇠!0 !)2
1 1
g(!) = 2
q
!0
(1 (!/!0 )2 )2 + 2⇠(!/!0 )2
!
X̃0 2⇠!0 !
Phase : (!) = Phase tan =
˜
f0 (!02 ! 2 )
Ou bien en utilisant la première notation utilisant ⌧ :

|X̃0 | 1
Gain : g(!) = =q
|f˜0 | (!02
2
! 2 ) + (!/⌧ )2
!
X̃0
Phase : = Phase
f˜0
!/⌧
tan =
(!02 !2)
2.6.2 Fréquence de résonance , largeur de bande

1. La fréquence de résonance correspond à la fréquence pour laquelle le gain est
maximum. L’amplitude des oscillations atteint alors sa valeur maximale pour
une intensité fixée de l’excitation extérieure.
2. Pour ⇠ ⌧ 1 , le gain atteintp son maximum pour ! ' !0 où le premier terme de
la racine carré s’annule (!02 ! 2 )2 + (2⇠!0 !)2 ' 2⇠!0 !.
3. On peut calculer la valeur exacte de la pulsation de résonance (!res ) à partir
de la dérivée de l’expression du gain. La dérivée s’annule lorsque le gain est
extrémum.
dg 8!!02 ⇠ 2 4! (!02 ! 2 )
= ⇣ ⌘3/2
d! 2
2 4! 2 !02 ⇠ 2 + (!02 ! 2 )
dg
= 0 ! 2⇠ 2 !02 = !02 !res
2
d! p
!res = !0 1 2⇠ 2
25
Figure 2.2 – La courbe de gain g(!) en fonction de !!0 pour trois valeurs du coefficient
d’amortissement sans dimension ⇠, en bleu pour ⇠ = 0.1, orange pour ⇠ = 0.15 et vert
pour ⇠ = 0.2.
Il n’y a pas de résonance si ⇠ > p12 et pour ⇠ ⌧ 1 , on a !res ' !0 (1 ⇠ 2 ),

inférieure à la fréquence propre de l’oscillateur amorti !res < !p .
4. A basses fréquences, pour ! ⌧ !0 , on a le déphasage qui s’approche de 0, et
le gain g ! !12 . Le mouvement de la masse suit les variations de la force, avec
0
F0
une amplitude déterminée par la raideur du ressort X0 ' k
.
5. Pour ! = !res , on obtient le gain maximum qui s’approche du facteur de qualité
Q fois le gain à très basses fréquences (Q 1 pour ⇠ ⌧ 1). Le déphasage
⇡
s’approche de 2 - la force ou l’excitation extérieure est alors en phase avec la
vitesse de déplacement.
1 1
g(!res ) = 2
p
2⇠!0 1 ⇠ 2
1 1
g(!res ) = ' 2
= g(0) = Qg(0) pour ⇠ ⌧ 1
2⇠!0 2⇠
6. A hautes fréquences, pour ! !0 , le déphasage tend vers ! ⇡ et le gain

décroit comme l’inverse du carré de la pulsation
1
g(!) ' pour ! !0 , ⇠ ⌧ 1
!2
26
7. Largeur de bande : La réponse de l’oscillateur reste importante dans une bande
de fréquence (ou de pulsation) autour de la fréquence de résonance. Par conven-
tion, on définit la bande passante comme la bande en fréquence ou en pulsation
[!1 , !2 ] tel que la réponse du système, ou le gain ici, soit supérieur ou égal à
p1 fois le gain à la résonance, ce qui correspond à avoir une puissance absorbée
2
supérieure à la moitié ( 12 ) de la puissance maximale absorbée. On se place dans
1
le cas où ⇠ ⌧ 1. Le gain à la résonance !res ' !0 est égal à g(!res ) = 2⇠! 2.
0
Les deux pulsations !1 , !2 définissant la bande passante doivent alors vérifier la
condition :
1 1
g(!1 ) ' g(!2 ) ' p
2 2⇠!02
2
(!1,2 !02 ) = ±2⇠!02
✓ ◆
1
!1,2 = !0 (1 ± ⇠) = !0 1±
2Q
!0 !0
! = (!2 !1 ) = ! Q=
Q !
On voit que le facteur de qualité Q, défini comme le rapport entre l’énergie

stockée dans l’oscillateur et celle dissipée en une période ( au facteur 2⇡ près),
correspond aussi au rapport entre la fréquence de résonance et la largeur de
bande, pour les faibles amortissements. Plus le facteur de qualité augmente, et
meilleure est la définition (précision) de la fréquence d’oscillation de l’OH. On
va donc rechercher un système oscillant avec un grand facteur de qualité pour
construire une horloge par exemple.
2.7 Impédance mécanique, impédance électrique
2.7.1 Impédance mécanique
Par analogie au cas électrique (Z = V /I), on définit l’impédance mécanique comme

le rapport entre la force harmonique appliquée et la vitesse v(t) = ẋ. Comme dans
le cas
p électrique,2 l’impédance dépend de la fréquence ou la pulsation Z(!). (Notation
j= 1 2 C, j = 1).
En partant de l’expression de x̃p (t), on calcule celle de la vitesse et on en déduit
27
l’impédance :
1 F0 j!t 1
x̃p (t) = X̃0 ej!t = e 2⇠!0 =
(!02 !2)+ 2j⇠!0 ! m ⌧
j! F0 j!t
ṽ(t) = j! X̃0 ej!t = e
(!02 2
! ) + 2j⇠!0 ! m
1 1 m
= j (!02 ! 2 ) + 2j⇠!0 !
v(! F0 !
On obtient alors l’expression de l’impédance mécanique :
F (!) m⇣ 2 2 !⌘
Z(!) = = j !0 ! + j
v(!) ! ⌧
m
Z(!) = j !02 ! 2 + 2j⇠!0 !
!
On peut écrire l’impédance en fonction du coefficient de frottement visqueux (fluide)

⌘, coefficient de proportionnalité de le force de frottement fluide par rapport à la
vitesse : :
m m 2 m 2
Z(!) = +j ! !02 = ⌘ + j ! !02
⌧ !
✓ ◆ ✓! ◆
m k k
Z(!) = + j m! = ⌘ + j m!
⌧ ! !
m
Z(! ! !0 ) ' =⌘
⌧
L’impédance à la pulsation propre dépend uniquement du coefficient de frottement

fluide.
2.7.2 Impédances électriques, circuit RLC

On peut facilement montrer que les impédances électriques, qui correspondent à une
généralisation de la notion de résistance électrique R, et de la loi d’Ohm U = RI, se
combinent comme des résistances, lorsque des dipôles électriques sont assemblés pour
constituer un circuit. L’impédance équivalente Zeq de deux dipôles Z1 , Z2 connectés en
série ou en parallèle peut être obtenue par les expressions suivantes :
Serie : Zeq = Z1 + Z2
1 1 1
Parallele : = +
Zeq Z1 Z2
28
Une résistance électrique a une impédance indépendante de la fréquence, et celle
d’un condensateur de capacité C et d’une bobine d’auto-inductance L s’écrivent :
1
ZR (!) = R ZC (!) = ZL (!) = jL!
jc!
On peut alors calculer facilement l’impédance d’une résistance, d’une bobine et d’un
condensateur connectés en série. La résistance peut aussi correspondre simplement à
la résistance parasite de la bobine. Soumis à une di↵érence de potentiel sinusoı̈dale de
pulsation ! et d’amplitude V0 , le circuit sera parcouru par un courant sinusoı̈dale de
même pulsation et d’amplitude I0 avec V0 = |ZRLC |I0 . Le déphasage entre le courant
et le potentiel est donné par la phase de ZRLC . L’impédance ZRLC peut s’écrire en
1
fonction de !02 = LC qui est la pulsation caractéristique d’un circuit LC. L’impédance
d’une auto-inductance et d’un condensateur connectés en série s’annule pour ! = !0 :
1
ZRLC (!) = R + jL! +
jC!
2
LC! 1
ZRLC (!) = R + j
C!
! /!02 1
2
ZRLC (!) = R + j
C!
(! 2 !02 )
ZRLC (!) = R + j
C!!02
On peut réécrire l’expression ZRLC ci-dessus en introduisant une deuxième pulsation

1
caractéristique !a = RC , en notant que ⌧ = RC est le temps caractéristique de décharge
d’un condensateur à travers une résistance.
1
ZRLC (!) = RC!) + j ! 2 /!02 1
C! ✓ ✓ 2 ◆◆
1 ! !
ZRLC (!) = +j 1
C! !a !02
29
Chapitre 3
Au delà de l’OH simple
3.1 Oscillations anharmoniques

1. Cas d’un terme d’ordre supérieur dans le potentiel (écart par rapport à la loi de
Hooke F = kx (on suit le développement du cours de Berkeley , Mécanique,
Chapitre 7
1 2 1
V (x) = kx ✏kx3
2 3
Fx (x) = kx + ✏kx2
On obtient l’équation di↵érentielle non linéaire (terme ✏kx2 )
mẍ + kx ✏kx2 = 0
ẍ + !02 x ✏x2 = 0
On va chercher une solution avec les premiers termes d’un développement en

série de Fourier. On ne garde que les termes d’ordre 1 au plus en ✏, e et
x(t) = e + A (cos(!t) + q cos(2!t)) e ⌧ 1, q ⌧ 1

2
ẍ = ! (cos(!t) + 4q cos(2!t))
1
x2 (t) ' A2 cos2 (!t) + . . . = A2 (1 + cos(2!t))
2
En remplaçant dans l’équation di↵érentielle, une équation avec un terme constant,

un terme en cos(!t) et un terme en cos(2!t).
✓ ◆ ✓ ◆
2 1 1 2
A ! !02 cos(!t) + A 2
4q! + q!02 2 2
✏A!0 cos(2!t) + !0 e ✏A =0
2 2
30
Les trois termes doivent s’annuler, on obtient alors :
!2 !02 = 0 ! ! = !0
1 1
3q ✏A = 0 ! q = ✏A
2 6
1 2 1
e ✏A = 0 ! e = ✏A2
2 2
Les oscillations libres s’e↵ectuent autour d’une position moyenne hxi = e
décalée par rapport à la position d’équilibre, à la même fréquence que l’OH,
sans le terme non linéaire. Par contre, l’équation de mouvement présente une
composante à la fréquence double (2nd harmonique).

1 2 1
x(t) = ✏A + A cos(!t) ✏ cos(2!t)
2 6
2. On peut traiter par exemple le cas du pendule simple, en ajoutant le terme

3
suivant dans le développement de Taylor de sin ✓ '= ✓ ✓3! .
m ` ✓¨ + mg sin ✓ = 0
✓ 3
◆
g ✓
✓¨ + ✓ = 0
` 6
Voir s’il est possible d’écrire une solution approchée sous la forme ✓(t) = A (sin(!t) + q sin(3!t))
3. Autre exemple : génération de la fréquence double dans des systèmes non linéaires
3.2 Oscillateurs couplés

1. On considère deux oscillateurs mécaniques identiques (x1 (t), x2 (t)) couplés à
travers un ressort de raideur q :
k q
x¨1 + x1 + (x1 x2 ) = 0
m m
k q
x¨2 + x2 + (x2 x1 ) = 0
m m
Ce système se diagonalise facilement en introduisant deux variables z1 (t), z2 (t)

1 1
z1 = p (x1 + x2 ) z2 = p (x1 x2 )
2 2
1 1
x1 = p (z1 + z2 ) x2 = p (z1 z2 )
2 2
31
En calculant la somme et la di↵érence des deux équations di↵érentielles ci-
dessus, on obtient :
k
z¨1 + z1 = 0
m
k + 2q
z¨2 + z2 = 0
m
On a deux modes d’oscillations avec les deux fréquences propres !a < !b
k k + 2q
!a2 = ; !b2 =
m m
2
x1 (t) = x2 (t) = xa (t) = A cos(!a t + ) ! z2 (t) = 0 z1 (t) = p xa (t)
2
2
x1 (t) = x2 (t) = xb (t) = B cos(!b t + ) ! z1 (t) = 0 z2 (t) = p xb (t)
2
2. Le cas de deux oscillateurs mécaniques distincts (x1 (t), x2 (t)) couplés à travers
un ressort de raideur q :
k1 q
x¨1 + x1 + (x1 x2 ) = 0
m1 m1
k2 q
x¨2 + x2 + (x2 x1 ) = 0
m2 m2
On peut écrire le système des deux équations di↵érentielles couplées sous forme
matricielle :
✓ ◆
d2 x1
Ẍ = = AX
dt2 x2
!
k1 +q q
m1 m1
A = q k2 +q
m2 m2
On diagonalise la matrice A et on obtient les deux valeurs propres, correspon-

dant aux carrés des deux pulsations propres :
k1 k2 q q
= + + +
m1 m2 m1 m2
k1 k2 + q(k1 + k2 )
=
m1 m2
2 1⇣ p
2
⌘
!a,b = ± 4
2
32
3.3 Battements
Voir Ref [2] , Berkeley-Ondes, pages 28-30 de la version française -
1 = A cos(!1 t) 2 = A cos(!2 t)
= 1+ 2 = A cos(!1 t) + A cos(!2 t)
On introduit les deux fréquences !M , !mod :

1 1
!M = (!1 + !2 ) !mod = (!1 !2 )
2 2
(t) = Amod (t) cos(!M t) Amod (t) = 2A cos(!mod t)
développer cette section
3.4 Chaı̂ne infinie d’oscillateurs

Voir Ref [3] , Chapitre 22, pages 393-397
On considère des masses m, reliées entre elles par des ressorts de raideur ↵, de
longeur au repos a. L’ensemble constitue une chaı̂ne infinie d’oscillateurs, qu’on suppose
déployée le long de l’axe Ox. A l’équilibre, les masses identifiées par un numéro n 2 Z
se trouvent aux positions xn = na. Les masses peuvent osciller autour de leurs positions
d’équilibrePet on note un (t) le mouvement de la masse numéro n. On applique le PFD
2
m dtd 2 ~r = F au mouvement un (t) d ela masse numéro n qui subit deux forces des
deux ressorts qui la relient aux masses n 1 et n + 1.
mün = ↵ (un un 1 ) ↵ (un un+1 )

↵
ün + (2un un 1 un+1 ) = 0
m
On cherche des solutions de type ondes sinusoı̈dales progressives. On utilise la notation
complexe (expliquer la notation complexe pour les ondes). A noter que cette onde, se
propageant vers x > 0 ne prend des valeurs que pour des valeurs discrètes xn = na.
un (t) = A exp (i(!t kxn ))
On injecte dans l’équation algébrico-di↵érentielle ci-dessus, en notant que :
xn 1= xn a xn+1 = xn + a ! un 1 (t) = eika un (t) un+1 (t) = e ika

un (t)
↵
!2 + 2 eika + e ika = 0
m
↵
! 2 + 2 (1 cos(ka)) = 0
m
33
2
On pose !m = 4↵
m
, on utilise la relation sin2 ✓ = 12 (1 cos(2✓)) :
✓ ◆
2 2 2 ka 2 ↵
! = !max sin !max =4
2 m
✓ ◆
ka
! = !max sin
2
3.4.1 Relation de dispersion

Des ondes longitudinales se propagent sur cette chaı̂ne. Seules les valeurs de la
fonction d’onde (x, t) = A exp (i(!t kx)) aux points xn = na ont un sens physique.
Lorsqu’on ajoute p 2⇡
a
à k, un (t) = (xn , t) reste inchangé :
2⇡ 0
k0 = k + p ! eik xn = eikxn epn2⇡ = eikxn
a
On considérera donc uniquement un domaine limité pour le nombrev d’onde k, pour
lesquelles la relation de dispersion s’écrit :
⇡ ⇡
k +
a a ✓ ◆
ka
!(k) = !max sin
2
Les valeurs négatives de k < 0 représentent une onde se propageant vers x < 0, alors
que les valeurs positives k > 0 correspondent à une onde se propageant vers x > 0.
Seule des ondes avec une pulsation ! < !max peuvent se propager sur cette chaı̂ne, au
delà, on a une onde évanescente.
3.4.2 Vitesse de groupe, vitesse de phase

on peut calculer la vitesse de phase et la vitesse de groupe :
✓ ◆
!(k) !max ka
v' = = ⇥ sin
k k 2
✓ ◆
d! a ka
vgrp = = sign(k) !max ⇥ cos
dk 2 2
Milieu presque non dispersif pour |ka|ll1, la vitesse de groupe atteint zéro pour
ka = ±⇡. Dans ce dernier cas, on a deux masses consécutives qui oscillent en opposition
de phase.
un (t) = ±A cos(!t) un 1 (t) = un+1 (t) = un (t)
34
Chapitre 4
Cordes vibrantes, équation d’onde

de D’Alembert
4.1 Les ondes, une approche descriptive

1. Quelques exemples de phénomènes ondulatoires : des vagues qui se propagent
à la surface de l’eau, les cordes vibrantes, le son, les ondes électromagnétiques
dont la lumière, les ondes sismiques. A l’exception des ondes électromagnétiques,
les autres exemples cités sont des ondes mécaniques qui se propagent dans un
milieu matériel élastique.
2. Une caractérisation inspirée de la référence [4] , Tamer Bécherrawy, chapitre 4 :
Une onde est une perturbation qui se propage, propagation de l’énergie et d’autres
grandeurs physiques (impulsion par exemple) . Onde caractérisée par sa direction
et sa vitesse de propagation . . .on retrouve certains phénomènes en commun avec
les particules : propagation en ligne droite (milieux homogènes), la réflexion ou
le transport de l’énergie. En revanche, l’extension spatiale, l’interférence et la
di↵raction sont spécifiques aux ondes
3. Lorsqu’une partie d’un milieu élastique est écartée de sa position d’équilibre,
cette perturbation est communiquée de proche en proche à d’autres parties du
milieu et elle se propage jusqu’à de longues distances. Une onde transporte
l’énergie, donc aussi l’information, sur de longues distances, sans déplacement
macroscopique de la matière formant le milieu élastique. On verra que cette
vitesse de propagation augmente avec la rigidité du milieu (coefficient d’élasticité
plus grande), et à l’inverse, elle diminue pour des milieux plus dense (masse
volumique plus élevée).
4. Une chaı̂ne infinie d’oscillateurs (voir chapitre 3) peut être utilisée pour com-
prendre la notion de propagation de proche en proche, à travers l’interaction
entre deux oscillateurs successifs.
35
5. Une onde est appelée transversale lorsque la perturbation, ou le mouvement
des éléments du milieu se fait perpendiculairement à la direction de propagation,
alors que ce mouvement est parallèle à la direction de propagation pour les ondes
longitudinales. Dans certains cas, le déplacement des éléments du milieu peut
avoir deux composantes, l’une parallèle et l’autre perpendiculaire à la direction
de déplacement, comme par exemple dans le cas de certaines vagues à la surface
de l’eau.
6. La propagation peut se faire à une dimension, le long d’une corde ou d’un ressort
par exemple, à deux dimensions, comme les vagues à la surface d’un étang, ou
à trois dimensions, comme c’est le cas du son dans l’air.
7. Les ondes sont caractérisées par la vitesse de propagation c dans le milieu : une
perturbation se produisant à t = 0 en un point O se retrouve en tous point M ,
à une distance d = ct de O au bout d’un temps t. Elles sont représentées par un
champ dépendant du temps (~r, t) (ou (x, t) à une dimension).
8. Les ondes obéissent au principe de superposition, le déplacement du milieu,
résultant de deux ondes 1,2 , atteignant un point M , s’écrit comme la somme
algébrique des deux ondes, ce qui donne naissance à des phénomènes d’in-
terférences.
(M, t) = 1 (M, t) + 2 (M, t)
4.2 Corde vibrante

1. Modélisation des déformations et des contraintes sur une corde tendue, infini-
ment souple. La corde étant souple (sans raideur), la tension de la corde applique
une force localement tangente à la courbe en tout point. On note µ la masse
linéique de la corde et T sa tension, avec T(x) représentant le vecteur force due
à la tension en un point de la corde.
On étudiera les mouvements de la corde dans le plan xOy, la corde étant tendue
selon Ox, sa position d’équilibre étant confondue avec l’axe Ox. On néglige ou
on ignore l’e↵et du poids (attraction gravitationnelle) sur le mouvement de la
corde. Chaque point de la corde est repéré par son abcisse x, et son déplacement
sera noté ~s(x, t), avec deux composantes, l’une selon Ox notée u(x) et selon Oy,
noté simplement Y (x).
On isole par la pensée un élément de la corde situé en x, de longueur x, on fait
tendre ensuite x vers zéro. La résultante des forces s’exerçant sur cet élément de
corde est la di↵érence entre la tension de la corde s’exerçant des deux côtés. On
2
projète la résultante des forces selon Ox et Oy, on applique le PFD (m ddt2~r = F)
à l’élément de corde. On note ↵(x) l’angle entre la tangente de la corde en un
point donné et l’axe Ox, ou la droite confondue avec la corde non perturbée.
36
On admettra que la tension T = |T| est constante le long de la corde et ne varie
pas au cours du temps.
2. En un point M de la corde de coordonnée x0 , la force de tension s’applique
selon la tangente à la corde en ce point, caractérisée par l’angle ↵. On note
T , T+ la force exercée par la partie gauche (x < x0 )de la corde sur la partie
droite (x > x0 ) en M et respectivement, celle exercée par la partie droite sur
la partie gauche. Les projections des deux forces selon Ox et Oy (voir la figure
4.1) s’écrivent :
Tx = T cos ↵ Ty = T sin ↵
Tx+ = T cos ↵ Ty+ = T sin ↵
3. On applique le principe fondamentale de la dynamique à l’élément de la corde
pour établir l’équation di↵érentielle, plus précisément deux équations aux dérivées
partielles (EDP ) satisfaite par u(x, t) et le déplacement transverse Y (x, t). On
va donc projeter la force de tension T selon les deux directions Ox, le long de
la corde non perturbée, et la direction transverse, correspondant à Y
@Y
tan(↵(x, t)) =
@x
@2
µ x u = Tx + Tx+ = T [ cos(↵(x)) + cos(↵(x + x))]
@t2
@2
µ x 2 Y = Ty + Ty+ = T [ sin(↵(x)) + sin(↵(x + x))]
@t
On se place dans l’approximation des petits angles ↵(x, t) ⌧ 1, et en ne gardant
que les termes de premier ordre
@Y
tan(↵) ' sin(↵) ' ↵ '
@x
cos(↵) ' 0
@2 @2
µ 2Y = T 2Y
 2 @t @x
@ T @2
Y =0
@t2 µ @x2
4.
@2
µ x 2Y = T [ sin(↵(x)) + sin(↵(x + x))]
@t
@Y
↵ ⌧ 1 ! tan(↵) ' sin(↵) ' ↵ '
@x
@2 T @2
Y = Y
@t2 µ @x2
37
Figure 4.1 – Représentation schématique de la force de tension T, tangente à la corde,
et de sa projection selon les deux axes Ox et Oy.
5. Equation d’onde : on vient d’établir dans un cas particulier, celui de la corde

tendue infiniment souple, et avec les approximations ci-dessus, une équation
di↵érentielle appelée équation d’onde ou équation de D’Alembert qui décrit la
propagation des ondes. On peut réécrire cette
q équation sous une forme plus
T
générale en définissant le paramètre c = µ
qui correspond à la vitesse de
propagation des ondes sur la corde.
@2 1 @2 T
Y (x, t) Y (x, t) = 0 avec c2 =
@x2 c2 @t2 µ
4.3 Equation de d’Alembert et la forme générale

des solutions
L’équation de propagation des ondes a été établie dans le cas à une dimension
par Jean Le Rond d’Alembert (Paris 1717-1783), mathématicien, physicien, philo-
sophe français à la suite de l’étude des vibrations d’une corde de violon. On re-
trouve cette équation lors de l’étude de nombreux phénomènes physiques, les ondes
électromagnétiques en particulier. Elle s’écrit en général en fonction d’un paramètre c
qui correspond à la vitesse de propagation des ondes (vitesse de phase). La perturba-
tion (déplacement de la corde, ou de la surface d’un liquide, surpression dans le cas
des ondes acoustiques, . . .) sera représentée par le champ (~r, t) , ou (x, t) pour un
problème à une dimension. Noter que le champ représentant l’onde, noté dépend
du point de l’espace et du temps, donc au moins deux variables, (x, t), ou 4 variables,
(x, y, z, t) dans le cas plus général de propagation à 3 dimensions.
38
4.3.1 opérateur D’Alembertien
Nous avons établi l’équation d’onde à une dimension :

@2 1 @2
(x, t) = 0
@x2 c2 @t2
Il s’agit d’une équation di↵érentielle aux dérivées partielles de second ordre, satisfaite
par l’onde (x, t), fonction de la coordonnée spatiale x et du temps t.
Pour une onde se propageant dans un milieu à 3 dimensions, l’équation di↵érentielle
ci-dessus s’écrit sous la forme générale ci-dessous, à l’aide de l’opérateur Laplacien
r2 . Celui-ci correspond à l’application de l’opérateur gradient à la l’onde, suivi de
l’opérateur divergence. En coordonnées cartésiennes, l’opérateur laplacien s’écrit :
@2 @2 @2
r2 = r.r = + +
@x2 @y 2 @y 2
l’équation de D’Alembert s’écrit en coordonnées cartésiennes à 3D :

1 @2
r2 (x, y, z, t) = 0
c2 @t2

@2 @2 @2 1 @2
+ + (x, y, z, t) = 0
@x2 @y 2 @y 2 c2 @t2
1 @2
r2 = ⇤ Operateur D0 Alembertien
c2 @t2
4.4 Solutions de l’équation de d’Alembert
4.4.1 Solution de l’équation d’onde à 1D
Considérons l’équation d’onde à une dimension satisfaite par une quantité :
@2 1 @2
=0
@x2 c2 @t2
39
On e↵ectue le changement de variable (x, t) ! (u, v) :
(x, t) ! (u, v) u = x + ct , v = x ct
@ @ @
= +
@x @u @v
@2 @2 @2 @2 @2
= + + +
@x2 @u2 @v@u @u@v @v 2
@2 @ 2
@ 2
@ 2
= +2 + 2
@x2 @u
2 @u@v @v
@ @ @
= c
@t @u @v
 2
@2 @ @2 @2
= c2 2 +
@t2 @u2 @u@v @v 2
@2 @2
On substitue dans l’équation di↵érentielle de départ : les termes @u2
et @v 2
s’annulent
et l’équation di↵érentielle se réduit à :
@2
= 0 ! (u, v) = f (u) + g(v)
@u@v
Les solutions de l’équation d’onde s’écrivent donc sous la forme de la somme de
deux fonctions dérivables quelconques f (u) et g(v) :
@2 1 @2
= 0 ! (x, t) = f (x + ct) + g(x ct)
@x2 c2 @t2
g(x ct) correspond à une onde qui se propage vers les x croissants, alors que f (x + ct)
représente une propagation dans le sens opposé, vers les x décroissants.
4.4.2 Ondes harmoniques

Une forme particulière des solutions de l’équation de d’Alembert s’écrit à l’aide des
fonctions harmoniques (sinus ou cosinus) :
(x, t) = A cos (kx !t) propagation vers x > 0
(x, t) = B cos (kx + !t) = B cos ( kx !t) propagation vers x < 0
Ou en notation complexe :
(x, t) = Aei(!t kx)
+ Bei(!t+kx)
La linéarité de l’équation de d’Alembert, associée aux propriétés de la transformée de
Fourier permet alors d’écrire n’importe quelle solution de cette équation di↵érentielle
comme une combinaison linéaire d’ondes harmoniques, plus exactement d’une intégrale
avec une infinité non dénombrable de fonctions harmoniques.
40
1. On considère une onde progressive sinusoidale (harmonique) qui s’écrit (x, t) =
0 cos (kx !t + '). On peut toujours absorber la phase à l’origine ' par un
changement de l’origine des temps t ! t t0 avec !t0 = '. On prendra donc
' = 0 dans la suite. La phase totale (x, t) correspond à l’argument de la
fonction harmonique.
(x, t) = cos (kx !t) = 0 cos (x, t)

0
⇣ ⇣ ! ⌘⌘
(x, t) = 0 cos k x t
k ✓ ◆
⇣ ! ⌘ k
(x, t) = kx !t = k x t =! t+ x
k !
2. L’onde harmonique est une fonction doublement périodique, avec une périodicité
temporelle T et une périodicité spatiale , appelée longueur d’onde ; On définit
également le nombre d’onde k :
2⇡ 1
T = = periode
! ⌫
!
⌫= frequence
2⇡
2⇡ !
= = T longueur d0 onde
k k
La pulsation ! correspond à la vitesse angulaire ! = 2⇡

T
, le taux de variation de
la phase par rapport du temps @t . Le nombre d’onde k = 2⇡ est analogue à la
@
vitesse angulaire, mais dans le domaine spatial. En e↵et, l’onde harmonique est
périodique selon la coordonnée spatiale et la longueur d’onde est la période
spatiale.
(x, t) = !t ± kx
(x + , t) = (x, t + T ) = (x, t) ± 2⇡
@ 2⇡
! = =
@t T
@ 2⇡
k = =
@x
3. On appelle la relation de dispersion l’expression de la pulsation en fonction

du nombre d’onde !(k). La vitesse de phase, v = !/k correspond à la vi-
tesse de propagation des ondes de pulsation ! dans le milieu. Cette vitesse est
indépendante de ! dans le milieux linéaires. La vitesse de groupe d!/dk corres-
pond quant à elle à la vitesse de propagation de l’information ou de l’énergie.
41
! = !(k) Relation de dispersion
!(k)
v = Vitesse de phase
k
d!
vgrp = (k) Vitesse de groupe
dk
4. Milieu non dispersif : la vitesse de groupe et la vitesse de groupe sont identiques

dans un milieu non dispersif, pour !(k) = ck où c est une constante positive :
c 2 R, c > 0. La relation de dispersion, la vitesse de phase et la vitesse de groupe
pour une corde avec c2 = T /µ où T est la tension de la corde, et µ sa masse
linéique :
s s
T T
!(k) = c k = k v' = vgrp =
µ µ
4.4.3 Détermination des constantes d’intégration à partir des

conditions initiales
On peut déterminer le mouvement de la corde à partir de la solution générale de
l’équation d’onde, et les conditions initiales, la position et la vitesse de la corde en tous
points à t = 0, définies par les deux fonctions h(x) et j(x) :
@
(x, t = 0) = h(x) (x, t = 0) = j(x)
dt
On cherche les deux fonctions f (z) et g(z) tel que
(x, t) = f (x ct) + g(x + ct)
Noter qu’on a exprimé les fonctions f et g en fonction de la variable z qui ne correspond
pas ici à la coordonnée spatiale, mais simplement l’argument des ces fonctions. La
variable z a été introduite pour être distincte de la cordonnée spatiale x et du temps t.
La vitesse de déplacement en tout point de la corde s’exprime en fonctions des
dérivées f 0 (z) et g 0 (z) :
@
= c ( f 0 (x ct) + g 0 (x + ct))
@t
Les conditions initiales à t = 0 se traduisent en contraintes pour les fonctions f, g.
f (z) + g(z) = h(z)
df dg
( f 0 (z) + g 0 (z)) c = j(z) f 0 (z) = ; g 0 (z) =
dz dz
42
On intègre la seconde équation , et on combine les deux nouvelles équations pour
obtenir les expressions de f (z) et g(z). Il faut se rappeler que la primitive (l’intégrale)
d’une fonction est définie à une constante près. Soit une fonction k(z) et sa primitive
K(z) :
Z b
dK
= k(z) ! k(z)dz = K(b) K(a)
dz a
Pour obtenir une expression explicite de f (z), g(z) en fonction h(z) et j(z), nous
sommes amenés à formuler une hypothèse supplémentaire. On admettra que la corde
revient à sa position d’équilibre pour x ! ±1, ce qui revient à dire que les fonctions
f (z), g(z), ainsi que h(z) et j(z) tendent vers zéro lorsque x tend vers ±1.
Z
1 z
f (z) + g(z) = j(z 0 )dz 0
c 1
Z z
1 1
f (z) = h(z) j(z 0 )dz 0
2 2c
Z z1
1 1
g(z) = h(z) + j(z 0 )dz 0
2 2c 1
• Vitesse initiale nulle :

Si la corde est lâchée avec une vitesse initale nulle à t = 0 , cad j(x) = 0, la solution
s’écrit :
1 1
(x, t) = h(x ct) + h(x + ct)
2 2
L’onde s’écrit comme la somme de deux ondes, chacune ayant la même forme que h(x)
, mais mais d’amplitude moitié (facteur 12 ), l’une se propageant vers la droite, vers les
x croissants, et l’autre, en sens inverse, vers la gauche. Les deux déformations sont
strictement superposées à t = 0, avec des vitesses opposées en tout point, ce qui se
traduit pas une vitesse totale nulle en tout point de la corde à t = 0.
4.5 Impédance mécanique d’une corde

On rappelle que l’impédance mécanique correspond au rapport entre la force et la
vitesse. On considère une onde se propageant vers les x croissants ! +1, sur une
corde infinie :
(x, t) = f (x ct)
Soit un point x0 de la corde, la partie gauche de la corde x < x0 exerce sur la partie
droite x > x0 une force F dont la composante Fy s’écrit (noter le signe moins) :
43
(Rajouter un schéma)
d
Fy = T sin ↵ tan ↵ =
dx
d
↵ ⌧ 1 ! Fy ' T tan ↵ = T
dx
d 1d
(x, t) = f (x ct) ! =
dx c dt
Td
Fy = +
c dt
p
Or, ddt est justement la vitesse de la corde, et en se rappelant que c = T /µ, on obtient
l’expression de l’impédance mécanique de la corde :
T d
Fy
Zcorde = = cd dt
vy dt
T p
Zcorde = = Tµ
c
L’impédance peut-être a↵ectée d’un signe si on veut l’utiliser dans des expressions
algébriques. Si on écrit l’expression de l’impédance pour une onde se propageant vers
x ! 1, on a :
d 1d
(x, t) = f (x + ct) ! =+
dx c dt
Td p
Fy = ! Zcorde = Tµ
c dt
44
Chapitre 5
Réflexion, ondes stationnaires,

modes propres d’oscillation
5.1 Corde fixée à une extrémité, réflexion

On considère les ondes qui se propagent avec la vitesse c sur une corde tendue, mais
fixée rigidement en un point x = L. La corde s’étend donc dans le domaine x < L et
la fixation de la corde en x = L impose une valeur nulle de déplacement de la corde en
x = L : (x = L, t) = 0. . On recherche comment adapter la solution de l’équation des
ondes afin de répondre à cette condition aux limites (x = L, 8t) = 0. On part de la
solution générale de l’équation de D’Alembert, écrite comme une somme d’une onde se
propageant vers les x croissants, définie par la fonction f (z) et une onde se propageant
vers les x décroissants, définie par la fonction g(z). rouge
(x, t) = f (x ct) + g(x + ct)
(x = L, t) = f (L ct) + g(L + ct) = 0
Ceci fixe donc la forme de la fonction g(z) en fonction de celle de f (z).
g(L + ct) = f (L ct) 8t ! g(z) = f (2L z)
On obtient g(z) à partir de la fonction f (z) à l’aide de trois opérations, inversion du
sens de l’axe z (argument de la fonction) , z ! z, un décalage le long de l’axe z :
f ( z) ! f (2L z), et enfin inversion du signe de la fonction, ce qui correspond à une
multiplication par 1 : f (2L z) ! f (2L z). L’e↵et de ces trois opérations qui
permettent d’obtenir g(z) à partir de f (z) est illustré sur la figure 5.1.
(x, t) = i (x, t) + r (x, t)
i (x, t) = f (x ct)
r (x, t) = g(x + ct) = f (2L x ct) = f (2L (x + ct))
45
Figure 5.1 – Illustration des opérations pour obtenir la fonction g(z) représentant
l’onde réfléchie à partir de la fonction f (z) correspondant à l’onde incidente.
On peut écrire l’onde totale comme la somme d’une onde incidente i qui se propage
vers les x croissants et qui subit une réflexion en x = L, donnant naissance à une onde
réfléchie r , qui se propage vers les x décroissants.
(x, t) = i (x, t) + r (x, t)

(x, t) = f (x ct) f (2L (x + ct))
5.2 Corde fixée aux deux extrémités

Soit une corde fixée aux deux extrémités, en x = 0 et en x = L. On part de
la solution obtenue ci-dessus pour une corde fixée à une extrémité, en x = L, et on
impose la condition de déplacement nul (0, t) = 0 en x = 0, en plus de celle en x = L.
(x, t) = f (x ct) f (2L (x + ct))

(0, t) = (L, t) = 0 !
f ( ct) f (2L ct) = 0 ! f (z) = f (z + 2L) 8z
La fonction f (z) est donc une fonction périodque de période 2L avec f (z) = f (z +
2L) 8z. Ceci entraı̂ne également une périodicité temporelle ⌧ = 2L
c
f (z) = f (z + 2L) ! f (ct) = f (ct + 2L)

✓ ✓ ◆◆
2L
! f (ct) = f c (t + = f (c ((t + ⌧ ))
c
2L c
! ⌧= ⌫=
c 2L
46
La fonction f (z) peut donc s’écrire sous forme d’une série de Fourier, somme de fonc-
tions harmoniques , sinus et cosinus, avec des harmoniques multiples d’une fréquence
c
de base ⌫1 = 2L
X h ⇣ z⌘ ⇣ z ⌘i
f (z) = a0 + an cos n⇡ + bn sin n⇡
n2N,n>0
L L
5.3 Ondes stationnaires

Pour trouver la forme de (x, t), pour une oscillation à une fréquence donnée, on
peut procéder de plusieurs manières. Dans un premier temps, nous allons simplement
écrire que l’onde incidente est une onde harmonique de pulsation ! et l’onde réfléchie a
la même pulsation, mais se propage en sens inverse, avec une amplitude opposée à celle
de l’onde incidente. On suppose en outre que la corde est fixée en x = 0 et s’étend vers
les x croissants, x > 0. Il est également utile de rappeler la relation entre la pulsation
2⇡ c
et le nombre d’onde ! = ck, et la longueur d’onde = = 2⇡ .
k !
(x, t) = i (x, t) + r (x, t)
i (x, t) = A cos (!t kx + ')
r (x, t) = A cos (!t + kx + ')
On peut noter que la fonction cos étant paire cos(✓) = cos( ✓), les expressions ci-
dessous pour i et r sont équivalentes :
cos(kx !t) = cos(!t kx) cos( kx !t) = cos(!t + kx)
L’expression de l’onde correspondant à la superposition des deux ondes i et r

s’annule en x = 0. On développe les deux cos en utilisant la relation trigonométrique
cos(a + b) = cos a cos b sin a sin b
et on obtient :
(x, t) = A [cos ((!t + ') kx) cos ((!t + ') + kx)]

(x, t) = 2A sin(kx) sin(!t + ')
(x, t) = A(t) ⇥ F(t)
On voit que (x, t) se met sous la forme du produit d’une fonction A(t) = 2A sin(kx)
de la cordonnée spatiale x et d’une fonction du temps F(t) = sin(!t + '). L’allure de
47
la corde pour une onde stationnaire est représentée sur le premier schéma de la figure
5.2, en haut. Il s’agit de l’expression d’une onde stationnaire. On appelle les noeuds
de vibration les points xnoeud où l’amplitude de l’onde s’annule : sin(kxnoeud ) = 0. Les
ventres de vibration correspondent aux points xventre où l’amplitude des vibrations est
maximale : sin(kxnoeud ) = 1. Les positions des noeuds et des ventres sont indiquées sur
le schéma en bas de la figure 5.2.
Noeuds : sin (kxnoeud ) ! kxnoeud = n⇡ , n 2 N

✓ ◆
1
Noeuds : sin (kxventre ) ! kxventre = +n ⇡ , n2N
2
✓ ◆
n⇡ 1
xnoeud = =n xventre = n +
k 2 2 2
Deux noeuds ou deux ventres successifs sont distants d’une demi longueur d’onde /2
et un noeud et un ventre d’un quart de longueur d’onde /4, comme indiqué sur la
figure 5.2.
5.4 Modes propres

Si on veut que l’expression de l’onde stationnaire obtenue ci-dessus corresponde à
une onde qui se propage sur une corde fixée aux deux extrémités, il faut imposer la
condition de déplacement nul en x = L aussi, celle en x = 0 étant déjà vérifiée.
(x = L, t) = 2A sin(kL) sin(!t + ') = 0

sin(kL) = 0 ! kL = n⇡ n 2 N, n 1
Seules les ondes dont ayant des fréquences ou pulsation particulière !n peuvent
exister sur une corde fixée aux deux extrémités, ce sont les modes propres d’oscillation
d’une corde fixée aux deux extrémités.
n (x, t)
= 2A sin(kn x) sin(2⇡⌫n t + 'n )
n (x = 0, t) = n (x = L, t) = 0 8t
2⇡ ⇡ ⇡c
kn = n =n !n = n
2L L L
L’expression ci-dessus des modes propres peut aussi être établie à partir de l’ex-
pression générale de l’onde en fonction de f (z) obtenue pour une corde fixée aux deux
extrémités et le développement en série de Fourier de f (z) donnée ci-dessus. On écrit
48
Figure 5.2 – Allure de la corde à di↵érents instants t1 < t2 < t3 < t4 < t5 pour une
onde stationnaire (schéma du haut). Enveloppe des oscillations de la corde pour une
onde stationnaire, et position des noeuds et des ventres (schéma du bas).
49
l’expression de pour chacune des composantes cos et sin de la série de Fourier :
(x, t) = f (x ct) f ( (x + ct)) = f (x ct) f ( x ct)

h ⇣ n⇡ ⌘ ⇣ n⇡ ⌘i
cos
n (x, t) = a n cos (x ct) cos ( x ct)
L L
cos
(x, t) = a n [cos (!n t k n x) cos (! t + k x)]
n
h ⇣ n⇡ ⌘ ⇣n n⇡ n ⌘i
sin
n (x, t) = bn sin (x ct) sin ( x ct)
L L
sin
n (x, t) = bn [ sin (!n t kn x) + sin (!n t + kn x)]
On développe les termes cos sin, et on obtient :

cos
n (x, t) = 2an sin(kn x) sin(!n t)
sin
n (x, t) = 2bn sin(kn x) cos(!n t)
n (x, t) = A sin(kn x) cos(!n t + 'n )
Il y a une infinité dénombrable de modes propres d’oscillation de la corde, chacune

de ces modes peut-être vue comme la superposition de deux ondes se propageant en
sens inverse. L’amplitude est fixé en chaque point. Il y a des noeuds n 1 où l’ampli-
tude s’annule et la corde est immobile, et n ventres d’oscillations (ou l’amplitude est
maximum.
n (x, t)
= An (x) cos(!n t + 'n )
An (x) = A sin(kn x)
Les modes propres d’oscillation sur une corde fixée à ses deux extrémités sont à
des fréquences ⌫n = n⌫1 , multiples entières de la fréquence du mode fondamentale
⌫1 = ⌫fond
c
⌫fond = ⌫1 = Pour n2N n>0
2L
c c
⌫n = n !n = n⇡
2L L
2L 2⇡ ⇡
n = kn = =n
n n L
5.5 Oscillations forcées, résonance sur une corde

tendue
Si on excite d’une corde fixée aux deux extrémités en appliquant une force extérieure
à une fréquence ⌫e , la réponse de la corde de la corde, ca.d l’amplitude de ses vibrations,
50
sera d’autant plus importante que la fréquence de l’excitation extérieure est proche de
la fréquence de l’un des modes propres de la corde. On peut calculer la réponse dans le
cas particulier où une extrémité de la corde de longueur L est fixée et qu’on impose un
mouvement oscillatoire d’amplitude a et de pulsation !e = 2⇡⌫e à l’autre extrémité :
f (x = 0, t) = 0 (x = L, t) = a cos(!e t)
!e
f (x, t) = A sin(ke x) cos(!e t + ') ke =
c
f (x = L, t) = A sin(ke L) cos(!e t + ') = a cos(!e t)
a
A= '=0
sin ke L
a
f (x, t) = sin(ke x) cos(!e t)
sin kL
L’amplitude des oscillations est d’autant plus importantes que l’extrémité excitée
se trouverait proche d’un noeud des ondes stationnaires . . .
!e
(x = L, t) = a cos(!e t) ke =
c
a
f (x, t) = sin(ke x) cos(!e t)
sin kL
51
Chapitre 6
Ondes longitudinales dans les

solides
6.1 Elasticité des solides

Lorsqu’on soumet un solide à une force, on constate un alongement ou un raccour-
cissement (variation de la taille) dans la direction de la force, selon qu’il s’agit d’une
traction ou d’une compression. On observe également une variation inverse dans le
plan transverse, mais de plus faible magnitude. Dans le régime élastique, où le solide
revient à sa forme et dimension initiale lorsque la contrainte est supprimée, la force est
proportionnelle à l’alongement relative. Le coefficient de Poisson, noté ⌫ caractérise le
rétrécissement dans le plan transverse, pour un allongement. Le coefficient de Poisson
est défini comme
`?
⌫=
`k
Considérons une barre métallique de longueur L0 et de section S. On soumet cette
tige à une force de traction F , appliquée selon l’axe de la barre et on constate une
élongation de la tige, avec la force qui est proportionnelle à l’allongement relatif. On
constate aussi que losrqu’on fait varier la section de la barre, l’intensité de la force
est proportionnelle à la section de la barre. On appelle contrainte ( = FS ) la force
exercée par unité de surface, homogène à une de pression. Le module d’Young 1 , noté
E en général, exprime la proportionnalité entre la contrainte et l’allongement relatif.
E est homogène à une pression, unité Pa.
F L
= = E
S L
1. Thomas Young, physicien britannique 1773-1829
https://en.wikipedia.org/wiki/Thomas Young (scientist)
52
Figure 6.1 – Illustration module d’Young, zone de déformation élastique, plastique,
rupture
Remarque : figure + explication domaine élastique, déformation plastique, rup-

ture
Cette proportionnalité peut s’expliquer si on considère une forme de potentiel d’in-
teraction atomique avec un minimum autour de la distance inter-atomique moyenne
a (solide cristallin). La force exercée entre deux atomes lorsque ceux-ci sont écartés
de leur position d’équilibre peut alors être exprimée à travers la loi de Hooke (pro-
portionalité à l’allongement). L’allongement a sera proportionnel à L/L. si on tient
compte de la densité de surface du nombre de chaı̂nes d’atomes, on obtient bien une
force proportionnelle à la section S.
L L
L= a! a=a
a L
Notons aussi que la contrainte est une quantité algébrique (signée). > 0 cor-
respond à un allongement ( L > 0) et < 0 correspond à un a raccourcissement
( L < 0). Si on isole une tranche par la pensée, elle subit une force de traction par les
tranches voisines pour > 0, donc orientée selon la normale sortante des deux surfaces
délimitant la tranche, et une force de compression pour < 0, donc orientée selon la
normale entrante.
6.2 Ondes longitudinales dans les milieux élastiques

On considère un milieu élastique , de masse volumique ⇢ et caractérisé par le module
d’Young E. On considère un barreau de forme cylindrique d’axe Ox. On note u(x, t)
53
le déplacement selon ox d’une tranche d’épaisseur infinitésimal située en x0 . Analysons
donc ce qui se passe pour une tranche d’épaisseur `, ` ! 0 située entre x0 et x0 + `. On
peut écrire la contrainte se développant dans cette tranche en fonction de u(x :
` ! ` + (u(x0 + `) u(x0 ))
` u(x0 + `) u(x0 ) @u
= '
` ` @x
` @u
(x, t) = E = E
` @x
Intéressons nous maintenant au mouvement de cette tranche d’épaisseur `, ell ! 0,

de section S , donc de masse m = ⇢S`. Elle est soumise à gauche et à droite, en x0 et en
x0 + ` à des contraintes T+ , T . On applique le principe fondamental de la dynamique
à cette tranche m~a = F~ :
X
Fx = T + T = S [ (x0 + `) (x0 )]
X @ @ 2u
Fx = S `=SE 2 `
@x @x
@ 2u @ 2u
⇢S` 2 = S E 2 `
@t @x
2 2
@ u E @ u
2
=
@t ⇢ @x2
Les perturbations de position dans un milieu élastique, notées ici u(x) dans le cas à
1D satisfont aussi l’équation des ondes
p de D’Alembert. Ce sont des ondes longitudinales
qui se propagent à la vitesse cs = E/⇢, et correspond à la propagation du son dans les
solides. Les contraintes locales (x, t), équivalent aux surpressions dans dans les fluides
(voir chapitre 7 ) satisfont aussi l’équation des ondes et correspondent à la propagation
du son dans les solides.
@ 2u ⇢ @ 2u E
⇤ u(x, t) = = 0 c 2
s =
@x2 E @t2 ⇢
@2 ⇢ @2
⇤ (x, t) = =0
@x2 E @t2
6.3 Impédance mécanique

L’impédance est définie comme le rapport de la contrainte (x, t) et de la vitesse de
déplacement v(x, t) = @u/@t. Les dérivées partielles du champ de déplacement u(x, t)
54
p
⇢ (kg/m3 ) E (GPa) Coe↵. Poisson ⌫ cs = E/⇢ (m/s)
Aluminium (Al) 2700 62 ⇠ 0.3 ⇠ 4800
Laiton (alliage Cu , Zn) 8470 ⇠ 100 ⇠ 0.35 ⇠ 3400
Fer (Fe) 7860 208 ⇠ 0.25 ⇠ 5150
Acier (alliage Fe, C) 7850 210 ⇠ 0.25 ⇠ 5150
Tungstene (W) 19350 ⇠ 400 ⇠ 4550
Verre pyrex 2320 ⇠ 50-100 0.20-0.30 4600 6500
Plexiglas (verre acrylique) 1180 ⇠3 ⇠ 1600
Caoutchouc 950 ⇠ 2.2 - 2.7 ⇠ 1600
Table 6.1 – Caractéristiques mécaniques pour quelques matériaux Masse vo-

lumique ⇢ en kg/m3 , le module d’Young E, en GPa = 109 Pa, et le coef-
ficient de Poisson ⌫, (sans dimension). La dernière colonne donne
p une valeur
approchée de la vitesse du son calculée à partir de cs = E/⇢. (source :
https://www.simulationmateriaux.com/ComportementMecanique/ )
par rapport au temps et espace vérifie la relation suivante pour une onde progressive :
s
E
u(x, t) = f (x ± ct) c=
⇢
@u @u
v(x, t) = = ±c
@t @x
@u
(x, t) = E
@x
E @u E p
@x
Z = = = ± = ± ⇢E
v ±c @u
@x
c
On obtient donc l’expression de l’impédance mécanique d’une onde longitudinale

se propageant dans un solide élastique :
p
Z = ± E⇢ (± pour cos(!t ± kx))
Notez qu’on peut considérer l’impédance comme une quantité algébrique, positive pour
une onde se propageant vers les x croissants (x ! +1) et négative pour une onde se
propageant vers les x décroissants (x ! 1).
55
6.4 Aspects énergétiques des ondes dans un milieu
élastique
6.4.1 Densité d’énergie
Les ondes longitudinales dans un milieu solide élastique ou les perturbations de
pression dans un fluide, correspondant aux ondes sonores (acoustiques) sont associées
à de l’énergie mécanique. Il y a en e↵et de l’énergie cinétique associée au mouvement
des éléments du fluide ou du milieu élastique, ainsi que de l’énergie potentielle, similaire
à celle stockée dans un ressort comprimé, associée aux perturbations de pression ou les
contraintes et allongements correspondant dans un solide élastique.
Considérons un élément de solide élastique de section S et de longeur L0 de module
d’Young E, subissant un allongement L. Lorsqu’on a un allongement L, il faut
exercer une force F = SE LL0 et pour étirer la barre, il faut fournir un travail W =
SE LL0 d( L). En intégrant sur L , avec 0 < L < L, on obtient l’expression du
travail à fournir, qui correspond donc à l’énergie potentielle élastique stockée dans le
volume V = S ⇥ L0 . On note Eelast l’énergie stockée dans le volume V = SL0 par sa
compression ou son étirement et ✏elast représente la densité d’énergie élastique (énergie
par unité de volume) :
Z L
L SE 1
Eelast = SE d( L) = ( L)2
0 L 0 L 0 2
✓ ◆2
W 1 L 1 2
✏elast = = E =
SL0 2 L0 2E
On peut donc écrire l’expression de la densité d’énergie potentielle élastique ✏elast
stockée (énergie par unité de volume) en fonction de des perturbations longitudinales
de position u(x, t). L’expression de la densité d’énergie cinétique ✏cin ( énergie cinétique
par unité de volume) s’obtient directement (noter que @u @t
est la vitesse de l’élément de
solide en mouvement :
✓ ◆2
1 @u
✏elast (x, t) = E
2 @x
✓ ◆2
1 @u
✏cin (x, t) = ⇢
2 @t
" ✓ ◆ ✓ ◆2 #
2
1 @u @u
✏(x, t) = ✏cin (x, t) + ✏elast (x, t) = ⇢ +E
2 @t @x
⇤ Egalité de l’énergie élastique et de l’énergie cinétique

On remarque que les deux densité d’énergie cinétique et élastique sont égales à tout
56
instant et en toute position pour une onde progressive se propageant vers les x croissants
(x ! +1) ou les x décroissants.
@u @u
u(x, t) = f (x ± ct) ! = ±c
"@t ✓ ◆@x ✓ ◆2 #
2
E @u @u
c2 = ! ⇢ =E
⇢ @t @x
✏cin (x, t) = ✏cin (x, t)
Cette égalité stricte n’existe que pour des ondes progressives ; en e↵et, si on prend
l’exemple d’une onde stationnaire qui s’écrit comme la superposition deux ondes se
propageant en sens opposé :
u(x, t) = u+ (x, t) + u (x, t) = A [cos(kx !t) cos(kx + !t)]
u(x, t) = 2A sin(kx) sin(!t)
@u
= 2Ak cos(kx) sin(!t)
@x
@u
= 2A! sin(kx) cos(!t)
@t
on voit aisément la di↵érence de comportement des deux densités d’énergie, dont les
variations sont décalées temporellement et selon x.
6.4.2 Flux d’énergie

On considère le travail des forces dues aux contraintes dans un solide élastique (ou
des forces de pression dans un fluide) sur une tranche de fluide de section. On calcule
la puissance P traversant une section unité, en se rappelant que = E @u @x
représente
la contrainte (force / unité de surface). La puissance développée par une force dont le
point d’appui se déplace à la vitesse ~v s’exprime comme F~ .~v . La puissance P qui passe
à travers une surface de section unité s’écrit donc :
@u @u @ 2u
P(x, t) = E = E (6.1)
@x @t @x@t
@u @u
P(x, t) = (x, t)v(x, t) (x, t) = E v(x, t) = (6.2)
@x @t
Soit une petite tranche de section S et de longueur ` avec ` ! 0. La variation

temporelle de l’énergie dans cette tranche correspond à la di↵érence du flux entrant et
sortant d’énergie :
@
(✏cin (x, t) + ✏elast (x, t)) S` = P P+ = P(x0 , t) P(x0 + `, t)
@t
@ @
✏(x, t) = P
@t @x
57
Calculons le gradient du flux d’énergie, et la dérivée temporelle de la densité
d’énergie :
 2
@ @ u @u @u @ 2 u
P = E +
@x @x2 @t @x @x@t
2
@ @ u @u @u @ 2 u
✏(x, t) = ⇢ 2 +E
@t @t @t @x @t@x
@2u
Or u(x, t) satisfait l’équation de D’Alembert, on remplace donc @t2
dans la seconde
2
équation par E⇢ @@xu2 et obtient alors
@ @
✏(x, t) + P(x, t) = 0
@t @x
P(x, t) représente donc le flux de puissance (énergie par unité de temps et par unité
de surface) transporté par l’onde longitudinale qui se propage selon Ox dans le milieu
élastique. Dans le cas général dune onde se propageant selon une direction définie par
le vecteur d’onde ~k, on peut définir le vecteur P ~ = P ~k et l’équation de conservation
|~k|
de l’énergie s’écrit sous forme locale ou intégrale sous les formes suivantes :
@ ~ r, t) = 0
✏(~r, t) + div P(~
@t
~ = r.P
div P ~ = @Px + @Py + @Pz
@x @y @z
ZZZ ZZ
@ ~
~ dS
✏(~r, t)d3~r = P.
@t V S
58
Chapitre 7
Acoustique, propagation du son

dans les fluides
Acoustique : du grec akouein , entendre - Relatif à la perception des sons. Partie de

la physique qui traite des sons et des ondes sonores. L’Acoustique, selon le dictionnaire
Larousse 1 en ligne est défini comme :
Science qui étudie les propriétés des vibrations des particules d’un
milieu susceptible d’engendrer des sons, infrasons ou ultrasons, de les
propager et de les faire percevoir.
L’étude des ondes sonores, qui correspondent à la propagation des perturbations de
pression dans un fluide s’apparente à l’étude des ondes longitudinales dans un milieu
élastique (voir section 6.2 ). En e↵et, on va voir ci-dessous que les perturbations de
pression p jouent le même rôle que les contraintes dans un solide. L’elastcité du
fluide peut être caractérisé par le coefficient de compressibilité qui est équivalent,
comme on va le voir ci-dessous, à l’inverse du module d’Young (E = 1 ).
7.1 Compressibilité des fluides

7.1.1 Cas du gaz parfait
Un gaz, plus généralement un fluide, est un milieu élastique. La figure 7.1 ci-dessous
représente schématiquement l’élasticité du gaz. Une augmentation du volume du gaz
entraı̂ne une baisse de pression, alors qu’à l’inverse, une diminution de volume s’ac-
compagne d’une augmentation de pression.
On rappelle l’expression de l’équation d’état d’un gaz parfait et la relation de La-
place, qui caractérise une évolution adiabatique réversible (isentropique) du gaz. La
1. https://www.larousse.fr/dictionnaires/francais/acoustique/798
Voir Aussi : https://fr.wikipedia.org/wiki/Acoustique
59
Figure 7.1 – Représentation schématique de la force de tension T, tangente à la corde,
et de sa projection selon les deux axes Ox et Oy.
constante des gaz parfaits R = 8.31J/K/mol correspond au produit du nombre d’Avo-

gadro par la constante de Boltzmann R = NA kB et l’indice adiabatique est égal au
Cp
rapport des chaleurs spécifiques à pression et à volume constant = :
CV
pV = nRT R = 8.31J/K/mol
pV = Cte = Cp /CV
On va utiliser ces deux relations pour exprimer la variation relative de pression en

fonction de celle du volume pour des évolutions isotherme, respectivement isentropique.
dp dV
isotherme : T = Cte ! d(pV ) = pdV + V dp = 0 ! =
p V
1 dp dV
isentropique : d(pV ) = V dp + pV dV = 0 ! =
p V
En partant de l’équation d’état écrite pour une mole, n = 1, et en introduisant la

masse molaire M, on peut exprimer la masse volumique ⇢ en fonction de la température
et la pression, et décrire l’état du fluide à travers les trois variables T, p, ⇢, ce qui permet
d’éliminer la variable volume V travailler localement en tout point du fluide.
V mol RT
n = 1 ! pV mol = RT ! p =
M M
M M p RT
⇢(T, p) = mol = p ! =
V RT ⇢ M
60
On se place dans le cadre de petites perturbations du fluide autour de l’état
d’équilibre p0 , V0 , T0 .
p T ⇢
p = p0 + p ⌧1 T = T0 + T T0
⌧1 ⇢ = ⇢0 + ⇢ ⌧1
p0 ⇢0
7.1.2 Autres fluides

dans le cas général d’un fluide, l’élasticité du milieu peut être caractérisé par le
coefficient de compressibilité qui caratérise une variation relative de volume dV /V
suiet à une variation de pression.
✓ ◆
1 @V
=
V @p
La valeur du coefficient dépend bien sûr du fluide, mais également des conditions
d’évolution. Deux type d’évolution sont en général considéré, au moins pour le gaz,
une compression ou détente isotherme ( T ) ou adiabatique. L’entropie ne varie pas
dans le cas d’une évolution adiabatique réversible : S = Cte ! S . La variation de
pression p induite par une variation de volume joue le rôle de la contrainte dans un
solide suite à son allongement ou compression. L’inverse du coefficient joue donc le
rôle du module d’Young :
L 1 V
=E p=
L V
Les coefficients de compressibilité isotherme et isentropique pour un gaz parfait s’écrivent :
1 dV 1 1
dp = ! T = S =
V p p
7.2 Ondes longitudinales dans les fluides

On peut faire un raisonnement similaire au cas du solide élastique, pour analyser
la propagation des perturbations dans le fluide. On considère donc un tube de fluide,
cylindrique de section S, d’axe Ox et on note u(x, t) le déplacement d’une tranche de
fluide par rapport à sa position d’équilibre. Ce tube de fluide et une tranche d’épaisseur
` est représenté sur la figure 7.2. Pour une tranche de fluide d’épaisseur `, situé à
l’équilibre entre x0 et x0 + `, la variation de volume s’écrit :
V @u
V = S (u(x0 + `) u(x0 )) =) = pour ` ! 0
V0 @x
61
Figure 7.2 – Représentation schématique du déplacement des couches de fluide u(x, t)
pour une onde sonore se propageant selon Ox.
7.2.1 Cas du gaz parfait

On se place dans le cas de petits déplacements d’éléments de fluide, créant des
perturbations p, V adiabatiques (sans échange de chaleur), qu’on suppose réversible,
donc isentropique :
dp dV
= (isentropique)
p V
p dp V
' '
p0 p V0
@
p(x, t) = p0 u(x, t)
@x
L’élasticité du gas parfait correspond au fait qu’une variation de volume provoque

une variation de pression proportionnelle à dV V
. Je rappelle que la perturbation de
pression p joue le rôle de la contrainte dans les solides élastisques, mais avec une
convention de signe opposée.
@
isotherme : p= p0 u(x, t) E ! p0
@x
@
isentropique : p= p0 u(x, t) E ! p0
@x
La tranche de fluide est soumise aux forces de pressions, mais seul le gradient rp des
forces de pression apparaı̂t dans la résultante des forces qui s’exercent sur la tranche.
On ne considère ici que des variations selon x donc le gradient se réduit à une dérivée
@
partielle par rapport à x (r ! @x ).
On applique le PFD de la dynamique à la tranche de fluide de section S et d’épaisseur
` avant perturbation, situé à l’équilibre entre x0 et x0 + ` ; sa masse m ne varie pas et
62
est égale à S`⇢0 :
@2 @
S`⇢0 2
u(x, t) = S [p(x0 ) p(x0 + `)] = S` p(x, t)
@t @x
@2 p0 @ 2
u(x, t) = + u(x, t)
@t2 ⇢0 @x2

@2 1 @2 2 p0 RT0
u(x, t) = 0 c s = =
@x2 c2s @t2 ⇢0 M
p0 est donc l’équivalent du module d’Young E, caractérisant l’élasticité du milieu

gazeux. La di↵érence de signe vient d’une convention di↵érente, p > 0 correspond à
une compression, alors que < 0 pour la compression d’un solide.
Les perturbations longitudinales de position du fluide se comportent comme une
onde et se propage à la vitesse du son cs . celle-ci s’exprime en fonction de la température
absolue du gaz T , sa masse molaire M, l’indice adiabatique et la constante des gaz
parfaits R
r
RT Cp
cs = R = 8.31 J/K/mol =
M CV
De même, les perturbations de pression p se propagent à cette même vitesse, p
s’exprimant comme la dérivée de u(x, t). En e↵et, le son correspond à des perturbations
adiabatiques de l’état thermodynamique du gaz, qui donne naissance aux perturbations
de pression qui se propagent jusqu’à nos oreilles. Les perturbations de pression satisfont
donc également l’équation de D’Alembert :
@u
p(x, t) = p0
 @x
@2 1 @2
p(x, t) = 0
@x2 c2s @t2
Il n’y a pas de dispersion pour la propagation des ondes acoustiques dans un gaz
parfait (toutes les fréquences se propagent à la même vitesse).
r
RT0
!(k) = cs k = k
M
7.2.2 Cas général des ondes sonores dans un fluide

Voir la référence [R3] , (Boquet, Faroux, Renault), chapitre 24 pour établir l’équation
des ondes satisfaite par la pression et la vitesse dans un fluide à partir des équations
d’Euler du mouvement dans un fluide.
63
q
RT0
M (g/mol) = Cp /CV cs = M (m/s)
Hydrogène (H2 ) 2 1.41 1325
Helium (He) 4 1.66 1017
Eau (vapeur, H2 O) 18 1.33 429
Air 29 1.4 347
Oxygène (O2 ) 32 1.4 330
Argon (Ar) 40 1.67 322
(CO2 ) 44 1.3 271
Table 7.1 – Masse molaire M en grammes, rapport des chaleurs spécifiques =

Cp /CV ou exposant adiabatique à qla température T0 ⇠ 300K, et la vitesse du son
RT0
calculée pour T0 = 300K avec cs = M
On peut exprimer l’élasticité du fluide à travers le coefficient de compressibilité .

Une compression qui correspond à une diminution du volume provoque une augmen-
tation locale de la pression ( V < 0 ! p > 0) et une détente se traduit par une
diminution de pression ( V > 0 ! p < 0).
✓ ◆
1 @V
T = Isotherme
V @p T
✓ ◆
1 @V
S = Isentropique
V @p S
En notant que la masse volume est inversement proportionnel au volume, on peut

écrire :
1 d⇢ dV
⇢/ ) =
V ⇢ V
✓ ◆ ✓ ◆
1 @⇢ @p 1 1
= ) =
⇢ @p @⇢ ⇢
✓ ◆ ✓ ◆
@p 1 1 @p 1 1
= =
@⇢ S ⇢ S @⇢ S ⇢ T
Dans le cas de perturbations de position des éléments de fluide, ceux-ci subissent
des forces de rappel vers la position d’équilibre proportionnelles aux perturbations de
pression, donc / 1/ . la vitesse de propagation du son va pourvoir s’écrire en fonction
de @p/@⇢. On a vu que les équations de mouvement pour une tranche de fluide prenaient
la même forme que celles des élements d’un milieu élastique, si on remplace le module
d’Young par : ✓ ◆
1 dV 1
dp = : E!
S V S S
64
. Les perturbations de pression p(x, t) s’exprime en fonction du champ de déplacement
à travers le coefficient de compressibilité :
1 @u
p(x, t) = (x, t)
S @x
La vitesse de propagation des ondes dans un fluide s’exprime en fonction du coefficient

de compressibilité isentropique S et de la masse volumique du fluide ⇢ :
1 1 d⇢ dV
E! ⇢/! =
S V ⇢ V
✓ ◆
@p 1 E 1 1
E!⇢ = !=
@⇢ S S ⇢ ⇢ S
⇣ ⌘ 1 1
@p
c2s = @⇢
=
S ⇢ S
Les perturbations de pression p(~r, t) dans un fluide satisfont donc l’équation de propa-
gation des ondes qui s’écrit sous la forme ci-dessous :
✓ ◆ ✓ ◆
2 1 @2 @p 1
r p(~r, t) = 0 c2s = =
c2s @t2 @⇢ S ⇢ S
@2 @2 @2
r2 = + +
@x2 @y 2 @z 2
Si on applique la relation ci-dessus donnant l’expression de la vitesse de propagation

cs du son dans un fluide au cas du gaz parfait, on obtient bien évidemment le résultat
déjà établie au paragraphe 7.2.1
✓ ◆
@p 1 M
c2s = = ⇢0 = p0
@⇢ s s ⇢0 RT0
r
1 RT0
isentropique : s = cs =
p0 M
7.3 Impédance acoustique

On défnit l’impédance acoustique comme le rapport entre la supression p et la
vitesse de déplacement du fluide :
p p
Zac = =
v (@u/@t)
65
Soit une onde acoustique d’amplitude A et de pulsation ! :
u(x, t) = A cos(!t kx) ! = cs k

@u
v(x, t) = = !A sin(!t kx)
@t ✓ ◆
1 @u k
p(x, t) = = A sin(!t kx)
s @x s
On obtient l’expression de l’impédance acoustique en fonction de la masse volumique

⇢ et du coefficient de compressibilité S :
p k 1 1 1 p
Zac = = = = ⇢ S
v ! s cs S cs
r
1 ⇢
Zac = = cs ⇢ =
cs S S
Cette expression peut aussi être obtenue dans le cas général, d’une onde harmo-
nique ou pas, satisfaisant l’équation de d’Alembert. En e↵et, les dérivées spatiales et
temporelles du champ de perturbation u(x, t) sont reliées entre elles :
@u @u
v(x, t) = = ⌥cS
@t @x
1 @u
p(x, t) =
S @x
p 1
Zac = =±
v cs S
A noter le signe le changement de signe entre la première et la 3ème équation, et le signe
moins (-) pour la première équation pour une onde se propageant vers les x croissants.
Dans le cas d’un gaz parfait, l’expression de l’impédance acoustique s’obtient aisément
en se rappelant que la pression statique s’écrit : p0 = RT 0
M 0
⇢ et que 1/ S = p0 :
r r
M RT0
gaz parfait : Zac = p0 = ⇢0
RT0 M
Finalement, on a l’expression de l’impédance acoustique en fonction de la masse volu-
mique du milieu ⇢ et de la vitesse de propagation des ondes cs :
r
p 1 ⇢
Zac = = ⇢ cs = =
v S cs S
On note que l’impédance peut-être considérée comme une quantité algébrique, po-
sitive (Z > 0) pour une onde se propageant vers les x croissants, et négative (Z < 0)
pour une onde qui se propage vers les x décroissants.
66
Chapitre 8
Intensité sonore, les sources sonores
8.1 Aspects énergétiques, flux d’énergie

8.1.1 Densités d’énergie cinétique et d’énergie potentielle
En suivant un raisonnement analogue à celui exposé dans la section 6.4, on peut
établir une expression la densité d’énergie dans un fluide en présence de petits déplacements
ou de perturbations de pression, en remplaçant le module d’Young E par 1/ S (le S
en indice de S signifie à entropie constante) :
✓ ◆
1 E @p 1
E! c2s = =! =
s ⇢ @⇢ s s ⇢0
On obtient alors l’expression de la densité d’énergie potentielle associée aux pertur-

bations de pression dans le fluide :
✓◆2
1 1 @u
✏elast (x, t) =
2 S @x
✓ ◆2
1 @u
✏elast (x, t) = p0 gaz parfait
2 @x
On peut aussi établir cette expression en considérant une tranche de fluide, d’épaisseur
infinitésimale `, de volume V = S`. La surpression s’écrit 1S @u @x
. Cette tranche de
fluide se comporte comme un ressort de longueur au repos `, et de constante de raideur
kfluide , auquel il faut appliquer une force Fp pour compenser la surpression et la compri-
mer ou la dilater. On peut donc calculer l’énergie potentielle stockée Eelast de manière
analogue à celle d’un ressort. Cette quantité d’énergie s’avère être proportionnelle au
67
volume, on en déduit alors la densité d’énergie élastique ✏elast :
@u
` = ` |Fp | = kfluide `
@x ✓ ◆
1 @u S
|Fp | = ! kfluide =
S @x S`
✓ ◆2
1 1 1 @u
Eelast = kfluide ( `)2 = S`
2 2 S @x
✓ ◆2
1 1 @u
✏elast =
2 S @x
Enfin, en adoptant un point de vue encore di↵érent, partant de l’expression du
travail des forces de pression | W | = |pdV |, et le fait qu’une variation du champ de
déformation @u
@x
provoque une variation du volume de la tranche de fluide, on peu écrire :
✓ ◆ ✓ ◆
@u @u
dV = Sd `=Vd
@x @x
✓ ◆ ✓ ◆
1 @u @u
W = d
S @x @x
Il faut fournir du travail pour comprimer ou dilater la transche de fluide. En
intégrant l’expression ci-dessus entre 0 et une valeur donnée de ⌥ = @u @x
on obtient
le travail W à fournir, qui est stockée sous forme d’énergie potentielle dans la tranche
de fluide de volume V . On obtient la densité d’énergie potentielle, ou encore densité
d’énergie élastique ou de pression en divisant par le volume :
Z Z ⌥ ✓ ◆
V @u
W = W = d =
S 0 @x
✓ ◆2
1 V @u
W =
2 S @x
✓ ◆2
W 1 1 @u
✏elast (x, t) = =
V 2 S @x
On écrit la densité d’énergie associée à l’onde sonore, comme la somme de la densité
d’énergie cinétique ✏cin (~r, t) et de la densité d’énergie élastique ou de pression ✏elast (~r, t) :
✏(~r, t) = ✏cin (x, t) + ✏elast (x, t)

" ✓ ◆ ✓ ◆2 # ✓ ◆
2
1 @u 1 @u 1 @⇢
✏(~r, t) = ⇢ + S =
2 @t S @x ⇢ @p s
✓ ◆
1⇥ 2 ⇤ 1 @u
✏(~r, t) = ⇢v + s p2 Note : p =
2 S @x
68
En se rappelant que p = ±Zac @u @t
et l’expression de Zac = ⇢S (voir le paragraphe
7.3), on peut montrer l’égalité de la densité d’énergie cinétique et élastique.
✏cin (~r, t) = ✏elast (~r, t) ✏(~r, t) = s p2
8.1.2 Flux d’énergie et intensité acoustique

Comme dans le cas de la propagation des ondes dans un solide élastique, on peut
écrire l’expression du vecteur flux d’énergie en fonction de la surpression p et de la
vitesse locale du fluide v ( Voir la section 6.4.2 , l’équation 6.2 en particulier) :
~ = pv p2
Jac = P Iac = pv = Zac v 2 =
Zac
Jac représente la densité de courant d’énergie acoustique, sa direction correspond

à la direction de propagation du son, et son amplitude I, appelé intensité acoustique
représente la puissance par unité de surface (unité W/m2 ) transportée par l’onde so-
nore. q
⇢
En utilisant les expressions de cs = p⇢1 S , et de Zac = S
, on peut montrer
facilement que l’intensité acoustique correspond au produit de la densité d’énergie
acoustique par la vitesse de propagation du son :
Iac = |Jac | = ✏(~r, t) ⇥ cs

r
2 1 2 S
Iac = s p ⇥ p = p ⇥
⇢ S ⇢
1
Iac = p2
Zac
La relation écrite ci-dessous est satisfaite par la densité d’énergie et la densité de
courant d’énergie afin d’assurer la conservation d’énergie. La variation temporelle de
la densité d’énergie est compensée par div Jac :
@
✏(~r, t) + r.Jac = 0
@t
En se plaçant dans le cas à une dimension, on vérifie aisément que le flux d’énergie
acoustique Jac défini ci-dessous satisfait l’équation de conservation d’énergie. Calculons
la variation temporelle de la densité d’énergie, ainsi que le gradient du flux d’énergie. Il
faut aussi se rappeler que p satisfait l’équation de D’Alembert pour relier les dérivées
temporelles et spatiales de la surpression p(x, t) :
69
@ @ 1 @
✏(x, t) = 2 S p p = 2 S cs p p
@t @t cs @t
@ 1 @
✏(x, t) = 2 p p Onde progressive
@t Zac @x
@ 1 @
Jac = 2 p p
@x Zac @x
8.1.3 Intensité sonore moyenne
L’expression de l’intensité acoustique obtenue ci-dessous correspond à la puissance

instantanée transportée par l’onde. Or, cette puissance a une évolution temporelle
rapide, suivant celle de la surpression. Pour un son, avec une fréquence de 1 kHz par
exemple, l’intensité acoustique varie entre zéro et une valeur maximum un miller de fois
par seconde. On est en général intéressé par la valeur moyenne de l’intensité acoustique,
moyennée sur un temps correspondant à quelques fois la période des fréquences les
plus basses présentes dans le signal sonore, donc typiquement une fraction de seconde.
On va montrer que l’intensité acoustique moyenne est égale à la moitié de l’intensité
acoustique pic, pour une onde sonore harmonique.
On considère donc une onde sonore harmonique de fréquence ⌫ et de période ⌧ =
1
⌫
. On note ap l’amplitude de la surpression acoustique, et on va calculer l’intensité
acoustique moyenne sur durée d’une période, en un point donné, à l’origine. Notons
que l’intensité acoustique moyenne ne dépend pas de la position du point, pour une
onde plane, une onde se propageant selon une direction donnée, comme celle considérée
ici. De même, l’ajout d’un phase dans la fonction sinus ou cosinus ne changerait pas le
calcul de la valeur moyenne sur une période.
p(x, t) = ap cos(2⇡⌫t ± kx)

✓ ✓ ◆◆2
a2p t
Iac (x = 0, t) = cos 2⇡
Zac ⌧
Z ✓ ✓ ◆◆2
a2p 1 ⌧ t
hIac i = ⇥ cos 2⇡ dt
Zac ⌧ 0 ⌧
Il est facile de montrer que l’intégrale de sin2 ou cos2 sur une période vaut 12 . Il suffit
en e↵et d’exprimer sin2 ✓ ou cos2 ✓ en fonction de cos(2✓), puis noter que l’intégrale de
70
sin ou cos sur un nombre entier de périodes est nulle.
cos(2✓) = cos2 ✓ sin2 ✓ = 2 cos2 ✓ 1 = 1 2 sin2 ✓
Z 2⇡
cos(2✓)d✓ = 0
0
Z 2⇡ Z 2⇡ ✓ ◆
1 2 1 1 1
Ic = cos ✓ d✓ = + cos(2✓) d✓ =
2⇡ 0 2⇡ 0 2 2
Z 2⇡ Z 2⇡ ✓ ◆
1 1 1 1
Is = sin2 ✓ d✓ = cos(2✓) d✓ =
2⇡ 0 2⇡ 0 2 2
On obtient alors l’expression de l’intensité acoustique moyenne hIac i.
1 a2p
hIac i =
2 Zac
Le plus souvent, on omet le terme moyenne et on note simplement Iac , la quantité
2
qui correspond à hIac i, et qu’on écrit sous la forme Iac = 12 Zpac , où p correspond à
l’amplitude de la surpression et pas sa valeur instantanée.
8.2 Niveau acoustique

2
L’intensité sonore Iac = 12 Zpac représente le flux de puissance (énergie par unité de
temps et unité de surface, unité SI W/m2 transporté par l’onde sonore. L’intensité
sonore est en pratique exprimée sur une échelle logarithmique, par rapport à un niveau
ref
de référence Iac , appelé aussi le niveau sonore et exprimé en dB :
✓ ◆
Iac ref
IdB = 10 log 10 ref
Iac = 10 12 W/m2
Iac
ref
Le niveau de référence Iac a été choisie pour être proche du seuil d’audition humaine
à 1 KHz. Rappelons que l’oreille humaine entend les sons dans la gamme 20 Hz à 20
kHz. Les courbes isosoniques, appelées quelques fois les courbes de Fletcher-Munson
correspondent à des contours représentant un même niveau de perception d’intensité
sonore par l’oreille humaine.
( voir https://fr.wikipedia.org/wiki/Courb isosonique
ou https://en.wikipedia.org/wiki/Equal-loudness contour )
— 30-40 dB représente un environnement peu bruyant
— 55-60 dB pour une conversation
— 80-100 dB environnement bruyant, cris, musique très fort
— 110-120 dB Outillage électrique ou pneumatique , tronçonneuse, avions
— 140 Avion à réaction au décollage
(voir https://noiseawareness.org/info-center/common-noise-levels/
71
Figure 8.1 – Courbes isosoniques, représentant des contours
d’égal niveau de perception humaine d’intensité sonore (source
https://commons.wikimedia.org/wiki/File:Lindos1.svg )
8.3 sources sonores

— La voie : le passage de l’air à travers les cordes vocales les fait vibrer. Le timbre de
la voie est donné en grande partie par les modes de résonance des cavités buccale
et nasale. Le timbre de la voie change si on inhale un gaz de densité di↵érente de
celle de l’air. Si on respire de l’hélium, la vitesse du son étant inversement propor-
tionnelle à la masse molaire, la voie devient plus aigue, alors qu’elle devient plus
grave si on respire un gaz avec une masse molaire plus élevée que l’air, l’hexfluo-
rure de sou↵re par exemple. Voir : https://www.youtube.com/watch?v=CWdjOpU-p6Y
— Cordes vibrantes : de nombreux instruments de musique (piano, guitare, vio-
lon . . .) utilise les cordes vibrantes comme source primaire du son. Les cordes
étant fixées aux deux extrémités, les vibrations ne peuvent se produire que
selon les modes propres d’oscillation de la corde - à la fréquence fondamen-
tale et ses harmoniques (voir le paragraphe 5.4 ). Les instruments de musique
combinent en général les cordes vibrantes avec une autre structure résonante,
permettant d’amplifier sélectivement certaines fréquences. Un son est perçu
agréablement par l’oreille si les sons qui la composent - les notes de musiques -
ont des fréquences avec des rapports rationnels ( fa /fb = m/n avec m, n 2 N).
Dans la gamme tempérée utilisée dans la musique occidentale, un ensemble de
douze notes, les sept notes principales (Do Ré Mi Fa Sol La Si) et les cinq
notes intermédiaires (Do# Ré# Fa# Sol# La# ) constitue les douze demi tons
d’une octave. La fréquence de deux demi-tons successifs est dans le rapport
72
12
p
2 = 21/12 . Les fréquences des notes sont doublées pour deux octaves succes-
sives. La3 = 220 Hz, La4 = 440 Hz.
Voir par exemple :
http://maths-sciences-lp.ac-amiens.fr/sites/maths-sciences-lp.ac-amiens.fr/IMG/dossier musique/musique dossier.pdf
https://lecompositeur.com/theorie/frequence-des-notes-en-hertz/
— Tuyaux sonores : Les instruments à vent (orgue, flûte, trombone, saxophone
. . .) utilisent des tuyaux, souvent fermés d’un côté comme résonateur. C’est la
longueur du tuyau qui détermine la fréquence fondamentale.
— Les instruments à percussion utilisent les mode propres de vibrations d’un solide
de forme déterminée, ou celles d’une membrane élastique tendue sur un support
rigide (tambour). La physique des modes de vibrations d’une membrane tendue
s’approche de celle d’une corde tendue fixée à ses deux extrémités, mais bien
plus riche et complexe, compte tenu de la nature à 2 dimensions de modes de
propagation des ondes.
— Enfin, les dispositifs électro-acoustiques (haut-parleur) permettent de produire
un son à partir d’un signal électrique. Ils peuvent être constitués d’une mem-
brane associée à électro-aimant ou un actuateur piézo-électrique. Le déplacement
de la membrane, source du son est contrôlé par le signal électrique appliqué.
8.3.1 Instruments à cordes

Les instruments à corde sont équipés de plusieurs cordes vibrantes qui sont la source
des sons, les di↵érentes notes de musique, caractérisée chacune par sa fréquence. Soit
une corde vibrante, de masse linéique µ, tendue avec une tension T , de longueur L,
fixée aux deux extrémités, Nous avons établi au chapitre 5 l’expression des ondes sta-
tionnaires sur une corde vibrante :
(x, t) = A sin(kx) sin(2⇡⌫t + ')
On a vu qu’en imposant les conditions aux limites, (0, t) = (L, t) = 0, on obtenait
une série de valeurs discrètes pour k, qui correspondent aux fréquences de résonance
de la corde.
n⇡ 2⇡ 2L
kn L = n⇡ ! kn = k= ! n=
s L n
T ccorde ccorde
ccorde = , = ! ⌫n = n
µ ⌫ 2L
2L 1 ccorde
n = = ⌫n = n = n⌫1 n 2 N, n = 1, 2, 3 . . .
n n 2L
La corde, en oscillant, produit un son comprenant la fréquence ⌫1 du mode fondamental
et de ses harmoniques, ⌫n , multiples entiers de ⌫1 . Pour produire les di↵érentes notes
73
de musique, on peut jouer sur la longueur, la tension, mais aussi la masse linéique de
la corde. Pour une longueur et un type de corde, masse linéique fixée, on peut accorder
la corde en ajustant la tension de la corde. Dans un piano, chaque touche actionne
un marteau qui vient frapper une, voir un ensemble de deux ou trois cordes. D’une
note à l’autre, la longeur et la tension de la corde varie, et des cordes de di↵érentes
épaisseur, donc de masse linéique, sont utilisées. Notez que les fréquences des modes
propres d’oscillation d’une corde tendue, donc des sons produits sont déterminées par
la vitesse de propagation des ondes sur la corde, et pas de la vitesse du son.
8.3.2 Instruments à vent, tuyaux sonores

On considère ici un modèle très simplifié des instruments à vent, où l’instrument est
assimilé à un tuyau, d’axe Ox, fermé ou ouvert aux extrémités. C’est l’établissement
d’ondes sonores stationnaires, qui se propagent dans l’air contenu dans ces tuyaux,
qui détermine les fréquences de résonance. Pour comprendre l’e↵et d’une extrémité
fermée ou ouverte du tuyau, il faut noter l’onde sonore est caractérisée par un champ
de déplacement u(x, t), et le champ de surpression associé p(x, t). La surpression étant
@
proportionnel à la dérivée spatiale @x du champ de déplacement, elle est déphasé de
⇡/2 par rapport à ce dernier pour une onde harmonique :
u(x, t) = au sin(kx + ') sin(2⇡⌫t)

p(x, t) = ap cos(kx + ') sin(2⇡⌫t)
Dans le cas d’une onde stationnaire, les noeuds de surpression correspondent aux
ventres de déplacement et vice versa. Le déplacement s’annule sur une extrémitée
fermée du tuyau, alors que la surpression est maximum, alors qu’elle s’annule à une
extrémité ouverte, où le déplacement est maximal.
Tuyau fermé : Noeud de déplacement Ventre de surpression
Tuyau ouvert : Ventre de déplacement Noeud de surpression
Les deux extrémités fermées ou ouvertes :

Pour un tuyau de longueur L avec les deux extrémités fermées, le champ de déplacement
s’annule en x = 0 et en x = L. Les modes propres du tuyau s’apparentent à ceux d’une
corde vibrante, comme pour le cas d’un tuyau avec les deux extrémités ouvertes. C’est
la surpression qui s’annule aux deux bouts ; sa forme d’onde est la même que celle sur
la corde fixée aux extrémités. Les longueurs d’onde et les fréquences de résonance ont la
même expression que celles déterminées ci-dessus, pour la corde vibrante. Par contre,
il faut remplacer ccorde par la vitesse de propagation du son dans l’air cs .
2L 1 cs
n = = ⌫n = n = n⌫1 n 2 N, n = 1, 2, 3 . . .
n n 2L
Pour un tuyau fermé des deux côtés, c’est le champ de déplacement qui doit s’annuler
74
Figure 8.2 – Champ de surpression p(x, t) en bleu et champ de déplacement u(x, t)
(orange) pour les trois premiers modes n = 1, 2, 3. La série à gauche correspond à un
tuyau fermé aux deux extrémités, alors que la série de droite correspond à un tuyau
ouvert aux deux extrémités.
aux deux extrémités, en x = 0 et en x = L.
u(x = 0, t) = u(x = L, t) = 0 8t ! ' = 0 ; sin(kL) = 0

⇡ 2L
kn = n ; n = n 2 N, n = 1, 2, 3 . . .
L n
Pour un tuyau ouvert des deux côtés, on a
⇡
p(x = 0, t) = p(x = L, t) = 0 8t ! ' = ± ; cos(kL ± ⇡/2) = 0
2
⇡ 2L
kn = n ; n = n 2 N, n = 1, 2, 3 . . .
L n
Les champs de surpression et de déplacement pour les trois premiers modes propres
d’oscillation d’un tuyau fermé ou ouvert aux deux extrémités sont représentés sur la
figure 8.2.
Une extrémité ouverte et une extrémité fermée :
On a dans ce cas un ventre de surpression du côté fermé et un noeud du côté ouvert,
ce qui se traduit par :
u(x = 0, t) = u(x = L, t) = 0 8t ! sin(') = 0 ; sin(kL + ') = 1

⇡ ⇡
' = 0 ; kL = + m⇡ = + n⇡ m = 0, 1, 2, 3 . . . , n = m + 1
2 2
⇡ 4L
kn = (2n 1) ! n= n 2 N, n = 1, 2, 3 . . .
2L 2n 1
75
Figure 8.3 – Champ de surpression p(x, t) en bleu et champ de déplacement u(x, t)
(orange) pour les trois premiers modes n = 1, 2, 3, pour un tuyau fermé à gauche et
ouvert à droite.
La distance entre deux noeuds ou deux ventres consécutifs est d’une demi longueur
d’onde /2, alors qu’un noeud et un ventre consécutif ne sont séparés que d’un quart
de longueur d’onde. Les longueurs d’onde des modes propres du tuyau satisfont alors
la condition :
1 3 5 2n 1
1 =L 2 =L 3 =L ... n =L
4 4 4 4
4L 4L 4L
1 = 4L 2 = 3 = 4 =
3 5 7
La longueur d’onde et la fréquence des modes propres d’un tuyau sonore fermé à une
extrémité et ouvert à l’autre s’écrivent donc :
4L
n = = (2n 1) 1 1 = 4L
2n 1
cs cs
⌫n = (2n 1) = (2n 1)⌫1 ⌫1 = n 2 N, n = 1, 2, 3 . . .
4L 4L
Contrairement au cas de la corde, c’est la vitesse du son dans le gaz qui détermine les
fréquences des modes de résonance du tuyau et donc des sons produits. La figure 8.3
montre la configuration des champs de surpression et de déplacement pour un tuyau
fermé à une extrémité et ouvert à l’autre.
76
Chapitre 9
Quelques phénomènes ondulatoires
9.1 Propagation des ondes à 3 dimensions

9.1.1 Equation des ondes à 3D
L’équation de propagation des ondes (~r, t) dans l’espace à trois dimensions ~r 2 R3
s’écrit :

1 @2
⇤ = r2 (x, y, z, t) = 0
c2 @t2
@2 @2 @2
r2 = + +
@x2 @y 2 @y 2
Rappelons que r2 correspond à l’opérateur Laplacien, appliqué à une fonction scalaire
f (x, y, z) 2 R :
f = r2 f = div (grad f ) = r. (r)
9.1.2 Ondes planes

Les solutions générales peuvent s’écrire sous forme de superposition d’ondes planes.
Une onde plane est un solution de l’équation de d’Alembert à trois dimensions, similaire
à la solution 1D, et définissant une onde se propageant dans une direction définie
par un vecteur unitaire ~u = (ux , uy , uz ). Les surfaces d’onde, qui correspondent au
lieu des points où l’onde prend la même valeur à un instant donné sont des plans
perpendiculaires à ~u.
q
(~r, t) = f (~u.~r ct) |~u| = u2x + u2y + u2z = 1
(x, y, z, t) = f (xux + yuy + zuz ct)
77
Une onde plane harmonique de pulsation ! s’écrit sous la forme :
⇣ ⌘
~
(x, y, z, t) = A cos k.~r !t ~k = k~u ! = c k
9.1.3 Ondes sphériques

On peut écrire l’opérateur Laplacien r2 en coordonnées sphériques, appliquée à
une fonction f (r, ✓, ') :
 ✓ ◆
2 1 @2 1 @ @f @ 2f
r f (r, ✓, ') = (rf ) + sin ✓ +
r @r2 r2 sin ✓ @✓ @✓ @'2
On peut aisément vérifier que les solutions à symétrie sphérique (r, ✓, ') = (r) de
l’équation de D’Alembert s’écrivent sous la forme :
A B
(r, ✓, ') = (r) = f (r ct) + g(r + ct)
r r
Il s’agit donc d’ondes sphériques divergentes, se propageant vers les r croissants et
convergentes. se propageant vers les r décroissants. Les surfaces d’onde sont de forme
sphériques, centrée sur l’origine (ou la source). L’amplitude de l’onde (déplacement,
surpression, . . .) varie comme l’inverse de la distance à l’origine (la source), alors que la
densité d’énergie et le flux d’énergie J(r) transportée par l’onde (ou l’intensité) varie
comme r12 ; l’intégrale du flux sur des sphères concentriques restent constante.
1
A(r) = A0 amplitude
r
1
J(r) = J0 2 intensite
ZZ r ZZ
J(r)dS = J(R) dS = J(R)4⇡R2 = 4⇡ J0
Sphere R Sphere R
9.2 E↵et Doppler

Lorsque la source d’une onde ou bien le récepteur ou l’observateur se déplace par
rapport au milieu dans le quel se propage une onde, la fréquence de l’onde perçue par
le récepteur est di↵érente de celle de l’onde émise. Cette fréquence apparaı̂t plus élevée
lorsque l’émetteur et le récepteur se déplace l’un vers l’autre en s’approchant, alors
que la fréquence perçue est plus faible lorsque l’émetteur et le récepteur s’éloignent
l’un de l’autre. Ce phénomène porte le nom de Christian Johann Doppler 1 (1803-
1853), physicien et mathématicien autrichien, qui a émis l’hypothèse selon laquelle la
1. Johann Doppler https://en.wikipedia.org/wiki/Christian Doppler
78
Figure 9.1 – Gauche : représentation schématique des ondes planes (source : wikipedia)
- Droite : représentation schématique des ondes sphériques
fréquence ou la longueur d’onde de la lumière dépend de la vitesse relative de l’émetteur

et de l’observateur.
Pour comprendre le phénomène, imaginons que l’émetteur émet une impulsion à
intervalle régulier, à t0 = 0, t1 = ⌧, . . . tn = n⌧ . Ces impulsions se propagent à la vitesse
c dans un milieu vers un récepteur situé à une distance d de la source à t = 0. La
fréquence à l’émission est égale à ⌫e = ⌧1 .
1. Cas 1 : émetteur et récepteur immobile. La distance entre les deux est également
constante, égale à d. Le récepteur reçoit les impulsions aux instants t0n = dc + n⌧ .
La fréquence de réception des impulsions est identique à celle de leur émission,
⌫r = ⌫e .
2. Cas 2 : récepteur immobile, émetteur se déplaçant avec la vitesse ve dans la
direction du récepteur. Une vitesse positive ve > 0 correspond à l’émetteur et le
récepteur qui s’approchent l’un de l’autre, alors ve < 0 correspond où ces deux
objets. La distance émetteur-récepteur varie au cours du temps et on note dn
cette distance aux instants d’émission des impulsions et qui s’écrit dn = d nve ⌧ .
Les impulsions seront reçus aux instants t0n :
dn d ve
t0n = tn + = n⌧ + n⌧
c c c
d ⇣ ve ⌘ ⇣ ve ⌘
t0n = + n⌧ 1 ) ⌧r = ⌧ 1
c c c ⇣
⌫e c ve ⌘
⌫r = = ⌫e ⌫r ' ⌫e 1 + pour ve ⌧ c
1 ve /c c ve c
La fréquence perçue ⌫r sera donc plus élevée si ve > 0, dans le cas où l’émetteur
s’approche du récepteur et plus faible dans le cas contraire, pour ve < 0, d’un
79
facteur c cve .
3. Cas 3 : l’émetteur immobile et le récepteur se déplace avec la vitesse vr selon
la direction le reliant à l’émetteur. vr > 0 s’il se rapproche de l’émetteur, vr < 0
dans le cas contraire. Les impulsions seront reçues aux instants t0n . Au moment
où l’impulsion n est émise , le récepteur se trouve à une distance dn = d vr n⌧ ,
mais il continue à se rapprocher durant la propagation de l’impulsion.
✓ ◆
0 d vr n⌧ d vr
tn = tn + = + n⌧ 1
c + vr c + vr c + vr
✓ ◆ ✓ ◆
0 d c c
tn = + n⌧ ) ⌧r = ⌧
c + vr c + vr c + vr
⇣ vr ⌘ c + vr
⌫r = ⌫e 1 + = ⌫e
c c
4. Cas 4 : Dans le cas où l’émetteur et le récepteur sont tous les deux en mouve-
ment, on combine les deux e↵ets ci-dessus, on multiplie donc les deux facteurs
de changement de fréquence. Si on définit les deux vecteurs unitaires parallèlle
à la ligne reliant l’émetteur (E) et le récepteur (R), le vecteur ~uER orienté de
l’émetteur vers le récepteur, et ~uRE orienté du récepteur vers l’émetteur, les deux
quantités algébriques ve et vr correspondront aux vitesses projetés de l’émetteur
sur ~uER et du récepteur sur ~uRE . Les deux expressions ci-dessous pour ve et vr ,
écrites sous forme de produit scalaire de deux vecteurs, sont valables de manière
générale, pour une direction quelconque des vecteurs vitesse de l’émetteur V~E
et du récepteur V~R . L’e↵et Doppler s’annule dont pour des mouvements de
l’émetteur et du récepteur perpendiculaire à la ligne ER.
ve = V~E .~uER vr = V~R .~uRE
Ces deux expressions définissent ve et vr selon la convention de signe définie
ci-dessus, avec des vitesses projetées positives si l’émetteur et le récepteur s’ap-
prochent et négatives s’ils s’éloignent. L’expression générale de l’e↵et Doppler
s’écrit, dans le cas où l’émetteur et le récepteur sont tous les deux en mouvement
s’écrit :
✓ ◆ ✓ ◆
c + vr 1 + vr /c
⌫r = ⌫e = ⌫e
c ve 1 ve /c
⇣ vr ⌘ ⇣ ve ⌘
⌫r ' ⌫e 1 + 1+ pour |ve | ⌧ c
c c
9.2.1 Notes sur l’e↵et Doppler

1. Les ondes mécaniques se propagent dans un milieu. Lorsque l’émetteur est im-
mobile, la forme des surfaces d’onde est inchangé, et c’est le récepteur qui se
80
déplace à travers les crêtes et les creux des ondes, alors que ceux-ci se déplacent
aussi. Lorsque l’émetteur se déplacent, les surfaces d’onde se déforment.
2. Les ondes électromagnétiques ne se propagent pas dans un milieu et c’est la

relativité restreinte qui s’applique. L’e↵et Doppler pour les ondes EM dépend
uniquement de la vitesse relative de l’émetteur et de l’observateur
3. Considérons une onde plane harmonique se propageant dans un milieu élastique,
un fluide par exemple, au repos dans le référentiel R. Soit un référentiel R’ en
mouvement par rapport à R avec la vitesse v~0 . L’expression de cette onde dans
le référentiel R où le fluide est au repos s’écrit :
⇣ ⌘
(~r, t) = A exp ~k.~r !t ~k = k~u
~r = r~0 + v~0 t
⇣ ⌘
(r~0 , t) = A exp ~k.r~0 (! ~k.v~0 )t
En utilisant l’expression de damns le référentiel R’ et en identifiant R’ au

référentiel où l’émetteur, respectivement le récepteur est immobile, on obtient :
✓ ◆
1
!r = ! k~u.v~r !e = ! ! !r = !e 1 ~u.v~r
c
✓ ◆
1
!e = ! k~u.v~e ! !r = ! = !e + k~u.v~e = !e 1 + ~u.v~e
c
Figure 9.2 – E↵et Doppler. Figure de gauche et au centre : représentation schématique

des surfaces d’onde pour une source immobile et une source en mouvement. Figure de
droite : représentation de l’e↵et Doppler, pour une source (émetteur) en mouvement (
(source : wikipedia) )
81
9.3 Onde de choc
Lorsque la source se déplace à une vitesse supersonique, les formules établies dans
le paragraphe 9.2 ne s’appliquent plus. Les perturbations de pression issues du mou-
vement de la source ne peuvent plus se propager et s’éloigner de la source, celle-ci se
déplaçant à une vitesse au moins aussi grande que le son (vitesse supersonique). Les
ondes successives interfèrent de manière constructive sur un cône dont le sommet est
confondu avec la source en mouvement. Ce cône est le lieu d’une variation brutale de
la pression appelé onde de choc. L’angle au sommet du cône ✓ peut être déterminé
en écrivant la condition s’isochronisme des ondes issues des positions successives de la
source. Supposons que la source se trouve à l’origine, en x = 0 à t = 0. Elle se déplace
à une vitesse v, supérieure à la vitesse c de propagation du son v > c. A un instant
t < 0, la source se trouvait en un point M de coordonnée xM = vt et l’onde émise à
t < 0 a parcouru une distance r = c|t| depuis M , et se trouve donc sur une sphère de
rayon r = c|t| et centrée sur M . L’onde de choc est le cône constitué par la surface
tangente à l’ensemble de ces sphères pour t < 0 :
c
d = v|t| r = d sin ✓ = c|t| ! sin ✓ =
v
L’intersection de ce cône avec un plan contenant la trajectoire de la source est représentée
sur la figure 9.3.
— Ondes de choc se formant en di↵érents points du fuselage d’un avion en vol
supersonique
— La lumière Cherenkov est produite par une particule chargée qui se déplace à
une vitesse plus grande que celle de propagation de la lumière dans un milieu
transparent. C’est un phénomène électromagnétique analogue à l’onde de choc
sonore
82
Figure 9.3 – représentation schématique de l’onde de choc à un instant donné t0 .
La source se trouve en S à t0 et les cercles représentent les surfaces d’ondes émises par
la source à di↵érents instants t < t0 , et les points noirs les positions successives de la
source à ces instants. Les droites SH et SH 0 (en rouge) représentent le lieu des points
aux di↵érentes surfaces d’onde, correspondant à l’onde de choc.
83
Chapitre 10
Réflexion et transmission
10.1 Onde à l’interface entre deux milieux
On étudie la réflexion et la transmission d’une onde incidente sur la surface de

séparation entre deux milieux élastiques. On suppose que la surface de séparation est
confondue avec le plan yOz ; (1) dénote le milieu situé dans le demi-espace x < 0 et
(2) celui du demi-espace x > 0. Exemples :
— Deux cordes de masses linéiques di↵érentes µ1 , µ2 , soudées en x = 0 et tendues

selon Ox avec une même tension T
— Propagation du son dans deux fluide séparés par une membrane de masse
négligeable située dans le plan yOz.
— l’interface entre deux lignes de transmission de signal électrique - exemple câble
coaxial connecté à une paire torsadée
Soit une onde progressive plane i (x, t) qui se propage selon Ox vers x ! +1 dans
le milieu (1). Au passage de la surface de séparation, elle donne naissance à une onde
réfléchie r (x, t) dans le milieu (1), se propageant vers les x décroissants et une onde
transmise t (x, t) qui se propage dans le milieu (2) vers les x croissants.
i (x, t)= A1 exp (i(kx !t)) x0

r (x, t) = B1 exp (i( kx !t)) x0
t (x, t) = A2 exp (i(kx !t)) x>0
1 (x, t) = i (x, t) + r (x, t) 2 (x, t) = t (x, t)
84
Figure 10.1 – Représentation schématique des ondes incidente, transmise et réfléchie
à l’interface entre deux milieux, d’impédance Z1 et Z2 .
On va imposer la continuité du déplacement, mais également des contraintes appliquées

au passage de l’interface en x = 0.
@ 1 @ 2
1 (0, t) = 2 (0, t) ! (0, t) = (0, t)
@t @t
@ 1,2
v1,2 (x, t) = = i! 1,2 (x, t)
@t
@ 1 @ 2
1 (0, t) = 2 (0, t) ! (0, t) = (0, t)
@x @x
Ces conditions de continuité impose l’égalité de la pulsation ! entre l’onde incidente,

l’onde réfléchie et l’onde transmise. En se rappelant la définition de l’impédance Z = v :
1 (x, t) = Z1 [vi (x, t) vr (x, t)]
2 (x, t) = Z2 vt (x, t)
Les deux conditions de continuité énoncées ci-dessus, celle du déplacement ou de la

vitesse (a), et ensuite celle de la constrainte ou de la surpression p s’écrivent alors
sous la forme :
(a) : A1 + B1 = A2
(b) : Z1 (A1 B1 ) = Z2 A2
85
On a deux équations et trois inconnues , mais on peut déterminer les rapports d’am-
plitude en réflexion r et en transmission t :
B1 A2
r= t=
A1 A1
1+r =t Z1 (1 r) = Z2 t
2Z1 Z1 Z 2
t= r=
Z 1 + Z2 Z1 + Z 2
10.1.1 Réflexion et transmission du flux d’énergie

Le flux d’énergie, ou puissance par unité de surface pour une onde longitudinale
dans un solide ou un son dans un fluide s’écrit :
2
P(x, t) = (x, t)v(x, t) ! P = Zv 2 =
Z
p2
Iac = p v = Zac v 2 =
Zac
Le flux d’énergie transportée par les ondes incidente Pi , réfléchie Pr et transmise Pt
s’écrivent donc :
Pi = Z1 A21 ! 2
✓ ◆2
Z1 Z2
Pr = Z1 B12 ! 2 = Z1 A21 ! 2
Z1 + Z2
✓ ◆2
2Z1
Pt = Z2 A22 ! 2 = Z2 A21 ! 2
Z1 + Z 2
On peut alors définir les coefficients de réflexion R et de transmission T , correspondant
à la fraction du flux d’énergie réfléchie par l’interface vers le milieu (1) et transmise au
milieu (2). On vérifie aisément que la somme des deux coefficients de réflexion et de
transmission en flux d’énergie est égale à 1, reflet de la conservation d’énergie. Le flux
incident d’énergie se divise en deux parties, une partie transmise vers le milieu (2) et
une partie renvoyée vers le milieu 1.
✓ ◆2
Z 1 Z2
R =
Z 1 + Z2
✓ ◆2
Z2 2Z1 4Z1 Z2
T = =
Z1 Z 1 + Z2 (Z1 + Z2 )2
R+T = 1
— Cas Z1 ⌧ Z2 , par exemple passage du son de l’air dans un solide. On s’approche

du cas d’une réflexion totale, pas d’onde transmise dans le milieu (2), et un
86
régime d’onde stationnaire dans le milieu 1 : (A2 ' 0, B1 ' A1 ). Toute l’énergie
de l’onde incidente est réfléchie vers le milieu 1 (R ' 1, T ' 0) :
Z1 Z1
T '4 R'1 4
Z2 Z2
— Cas Z1 ' Z2 , deux milieux d’impédance proches, deux gaz de masses mo-
laires similaires par exemple séparés par une membranne de masse surfacique et
d’élasticité négligeable. On s’approche du cas d’une transmission totale (t = 1)
et pas de réflexion : (A2 ' A1 , B1 ' 0). Toute l’énergie de l’onde incidente est
transmise vers le milieu 2 et continue à s’y propager R ' 0, T ' 1.
1 Z1 Z2
Z= (Z1 + Z2 ) e= e⌧1
2 2Z
Z1 = Z(1 + e) Z2 = Z(1 e)
T = 1 e2 R = e2
— Cas Z1 Z2 , par exemple passage du son d’un solide dans l’air. On a un

coefficient de transmission en amplitude qui est proche de l’unité (t ' 1), tout en
ayant une amplitude de réflexion égale à l’amplitude incidente ; (A2 ' 2A1 , B1 '
A1 ). On a la presque totalité de l’énergie de l’onde incidente qui est réfléchie
vers le milieu 1, tout en ayant une fraction transmise, avec une amplitude de
déplacement de l’onde transmise égale à presque deux fois l’amplitude de l’onde
incidente (R ' 1, T ' 0) :
Z2 Z2
T '4 R'1 4
Z1 Z1
87
Chapitre 11
Interférences, Di↵raction
11.1 Description ondulatoire de la lumière

La lumière est une onde électromagnétique, une combinaison d’un champ électrique
et d’un champ magnétique oscillant, et se propageant dans l’espace (le vide) ou dans
un milieu diélectrique. L’onde électromagnétique est une onde transverse, les champs
électrique E et magnétique B se trouvant dans un plan perpendiculaire a la direction
de propagation. Ces champs étant vectoriels, l’onde électromagnétique ne peut être
totalement caractérisée par sa fréquence et sa direction de propagation et il faut en
outre préciser son état de polarisation (rectiligne, circulaire, elliptique, non polarisé).
Le comportement et l’interaction des ondes électromagnétiques varie énormément
avec leur fréquence, depuis le domaine radio, aux plus basses fréquences, suivi des
micro-ondes ( ⇠ mm), puis l’infra-rouge, le domaine visible, puis l’ultra-violet. Aux
plus hautes fréquences, on trouve les rayons X et . Le domaine visible correspond au
domaine de sensibilité de l’oeil humain et représente une toute petite partie du spectre
électromagnétique avec 2 [0.4 0.8] µm.
Les di↵érentes longueurs d’onde sont perçues en tant que couleur par le cerveau
humain, et les couleurs de l’arc en ciel vont du rouge aux plus grandes longueurs
d’onde, ⇠ 0.8µm jusqu’au bleu ⇠ 0.4µm, en passant par l’orange, le jaune, et le
vert.
La lumière se propage à la vitesse c = 299792458 m/s dans le vide. Dans un milieu
diélectrique (transparent), la lumière se propage à la vitesse c/n où n > 1 est l’indice
du milieu. Celui-ci varie avec la fréquence, ce qui explique la formation des arc-en-ciel
et la dispersion de la lumière par un prisme.
L’optique géométrique, est une approximation de l’optique géométrique dans le cas
des ondes planes, valable loin des sources lumineuses. Un faisceau lumineux correspond
à un ensemble de rayons lumineux qui se propagent en ligne droite dans un milieu
transparent. Les lois de la réflexion et de la réfraction décrivent le comportement des
88
rayons lumineux incidents sur l’interface de deux milieux d’indices n1 et n2 .
⇣ ⌘ ⇣ ⌘
î = r̂ n1 sin î1 = n2 sin î1
La théorie ondulatoire de la lumière a été développé par Christian Huygens (1629-

1695, La Haye, Pays-Bas), mathématicien, astronome et physicien néerlandais, perfec-
tionné par le physicien français Augustin Jean Fresnel (1788-1827). Chaque point d’une
surface d’onde se comporte comme une source secondaire d’ondes sphériques selon le
principe d’Huygens-Fresnel. La surface d’onde après propagation est obtenue par la
surface tangente à l’ensemble des ondes sphériques.
11.2 Interférences
11.2.1 Superposition de deux ondes, e↵et de la di↵érence de
phase
On considère la superposition de deux ondes de même amplitude et de même
fréquence, mais avec une di↵érence de phase :
✓ ◆
1 (t) = A cos(!t) = A cos !t +
2 2
✓ ◆
2 (t) = A cos(!t + ) = A cos !t + +
2 2
✓ ◆ ✓ ◆
(t) = 1 (t) + 2 (t) = 2A cos cos !t +
2 2
L’amplitude de l’onde résultante 2A cos 2 varie entre 0 , pour = ⇡(+2k⇡) à ±2A

pour = 2⇡(+2k⇡). L’intensité de l’onde I (sonore, lumineuse), qui correspond à la
densité d’énergie et le flux de puissance est proportionnelle au carré de l’amplitude,
l’intensité varie donc entre 0 et 4 fois l’intensité d’une seule onde :
2
I1,2 / 1,2 I(t) / (t)2
✓ ◆ ✓ ◆
2 2 2
I(t) / 4A cos cos !t +
2 2
✓ ◆
I(t) = I( ) cos2 !t +
2
✓ ◆
I( ) / 4A2 cos2 = 2A2 (1 + cos )
2
89
Pour deux ondes provenant d’une même source, le déphasage dépend de la di↵érence
de marche L entre les deux ondes ( di↵érence de distance de propagation ).
L
= kL = 2⇡
1 2 3 4 5 6
-1
-2
Figure 11.1 – Variation de l’amplitude (courbe bleue) et de l’intensité, carré de

l’amplitude (courbe orange) en fonction du déphasage lors de l’interférence de deux
ondes.
Note : Conditions d’interférence , cohérence entre les deux ondes
11.2.2 Expérience des trous d’Young

Une onde plane incidente (lumière) se propage selon Ox et arrive sur un écran,
parallèle à yOz et placé à l’origine O, en x = 0 perpendiculairement à Ox. L’écran est
percé de deux fentes, parallèles à Oy et situées en z = ± a2 . Un écran k yOz est aussi
placé perpendiculairement à Ox, à une distance de l’écran, en x = L avec L a. Les
deux fentes peuvent être assimilées à deux sources cohérentes d’ondes. Soit un point
M de l’écran, à une distance r de l’origine, de coordonnée xM = L, zM = z, ou (r, ✓)
en polaire. On note r1 , r2 les distances de M aux deux fentes.
p
s r = L2 + z 2 z = rssin ✓
✓ ◆2 ✓ ◆2
2
a2 2
a2
r1 = L + z r2 = L + z +
2 2
r r
1 a2 1 a2
r1 = r 1 a sin ✓ + 2 r2 = r 1 + a sin ✓ + 2
r 4r r 4r
On calcule la di↵érence de marche r = r2 r1 des ondes provenant des deux fentes
jusqu’en M , en e↵ectuant un développement limité à l’ordre 1 en (a/r). La condition
90
d’interférence constructive se traduit alors par la condition suivante sur les angles ✓n :
r = r2 r1 ' a sin ✓
rn = n ! sin ✓n = n
a
On obtient l’expression de l’amplitude de l’onde en un point de l’écran comme la
superposition (interférence) des deux ondes issues des deux fentes. On note = c/⌫ la
longueur d’onde de l’onde et on se place dans l’approximation des petits angle ✓ ⌧ 1.
⇣⇡ ⌘
A(zM ) = 2A0 cos r
⇣⇡ ⌘
A(zM ) = 2A0 cos a sin ✓
On se place dans l’approximation des petits angle ✓ ⌧ 1 :
✓ ⌧ 1 ! sin ✓ ' tan ✓ ' ✓

zM zM
tan ✓ = ! ✓'
L L
l’amplitude et l’intensité au point M de l’écran de coordonnée zM s’écrivent sous la
forme :
⇣ ⇡a z ⌘
M
A(zM ) ' 2A0 cos
✓ ✓ L ◆◆
2⇡a zM
I(zm ) ' I0 1 + cos
L
On obtient les expressions ci-dessous pour les positions znc des franges brillantes
(maximum d’intensité), et celles des franges sombres znd , correspondant à un minimum
d’intensité et à des interférences destructives.
zn
✓n ' 'n
L a
znc ' n L
✓a ◆
d 1
zn ' n+ L
2 a
11.2.3 Interféromètre de Michelson

L’interféromètre de Michelson (Albert Abraham Michelson, Etats-Unis, 1852-1931)
est constitué d’une lame semi-réfléchissante M0 (la lame séparatrice) et deux miroirs
M1,2 situés aux distances L1,2 de la séparatrice, le long des deux bras du dispositif.
Un faisceau de lumière monochromatique cohérente est séparé en deux par M0 , et
91
recombinés dans le bras de sortie après réflexion sur les deux miroirs. La di↵érence de
marche entre les deux rayons est égale au double de la di↵érence de longueur entre les
deux bras. L’amplitude et l’intensité de la lumière dans le bras de sortie s’écrivent :
2⇡
= 2 (L2 L1 )
✓ ◆
2⇡
A = A0 cos (L2 L1 )
✓ ◆
2 2⇡
I = I0 cos (L2 L1 )
✓ ✓ ◆◆
4⇡
I = I0 1 + cos (L2 L1 )
Figure 11.2 – Interféromètre de Michelson. schéma à gauche

(https://www.renishaw.fr/fr/comment-fonctionnent-les-systemes-interferometriques--38612),
image à droite ( https://www.pierron.fr/interferometre-de-michelson.html )
11.3 Di↵raction
Le phénomène de di↵raction correspond à une déformation ou distorsion des sur-
faces d’onde lorsque celle-ci rencontre un obstacle. La lumière (ou le son) semble se
répandre dans di↵érentes directions lorsqu’elle passe à travers une petite ouverture.
On observe de même des bords flous des ombres portées des objets de petite taille.
92
La di↵raction peut se comprendre à partir du principe d’Huygens-Fresnel évoqué
ci-dessus. Si on considère une onde plane arrivant sur un obstacle, par exemple un
écran percé d’un trou de petit diamètre a ou une fente rectangulaire a ⇥ b, chaque
point de l’ouverture se comporte comme une source secondaire d’onde émettant dans
toutes les directions. Ces ondes secondaire interfèrent. L’amplitude en un point donné
sera déterminée par le résultat de l’interférence de toutes les ondes provenant de tous
les points de l’ouverture.
Lorsqu’on se place loin de l’obstacle ou de l’ouverture, on est dans le cas de di↵rac-

tion en champ lointain ou l’approximation de Fraunhofer. Dans le cas contraire, pour
un point situé proche de l’ouverture, on serait dans le cas de la di↵raction de Fresnel.
En champ lointain, on peut calculer l’intensité dans une direction donnée en calcu-
lant l’amplitude résultant de l’interférence d’ondes plane issues des di↵érents points de
l’ouverture. Le point clé est alors la détermination de la di↵érence de marche, ou le
déphasage des ondes issues des di↵érents points en fonction de la direction de propa-
gation.
On se place dans des conditions similaires à l’expérience des trous d’Young. Une
onde plane de pulsation ! = 2⇡⌫ et de nombre d’onde k = 2⇡/ = !/c se propageant
dans la direction Ox arrive sur un écran confondu avec le plan yOz. L’écran est percé
d’une fente rectangulaire de grande longueur b selon Oy (b ! +1), mais de faible
largeur a selon Oz, a/2 < z < a/2. On se place loin de la fente, dans une direction
~u = (cos ✓, 0, sin ✓) dans le plan xOz, caractérisé par l’angle ✓.
Soit un point M de la fente, de coordonnée z, la di↵érence de marche de l’onde

issue de ce point par rapport à l’onde issue O s’écrit `(z) = z sin ✓. On note u la
coordonnée le long de la direction ~u à partir de l’origine O et on utilise la notation
complexe. L’onde EM (t, u) issue du point M de la fente de coordonnée z s’écrit :
`(z) = z sin ✓
EM (t, u) = A0 exp (i (!t k(u `)))
✓ ✓ ◆◆
z sin ✓
EM (t, u) = A0 exp i !t ku + 2⇡
On a un maximum d’intensité I / E 2 lorsque toutes les ondes sont en phase, donc

pour ✓ = 0 et un minimum quand une fraction notable des ondes partielles sont en
opposition de phase.
93
L’amplitude complexe de l’onde A(✓) dans la direction ✓ s’écrit :
Z +a/2
E(t, ✓) = EM (t, ✓) dz = A(✓) exp (i!t)
a/2
Z +a/2
A(✓) = A0 dz exp (ik `(z))
a/2
Z +a/2 ✓ ◆
sin ✓
A(✓) = A0 dz exp i2⇡z
a/2
a sin ✓
On pose z 0 = z/a et s = et ; on reconnaı̂t l’expression de la transformée de
Fourier d’une fonction porte :
Z +1/2
A(✓) = A0 a dz 0 exp ( 2i⇡z 0 s)
1/2
A0 sin (⇡s)
A(✓) =
a ⇡s
A0 sin ⇡ a sin ✓
A(✓) =
a ⇡ a sin ✓
L’intensité de l’onde correspond au carré de l’amplitude. L’intensité lumineuse suite

à la di↵raction de la lumière par une fente rectangulaire allongée (b ) varie donc
selon la l’angle ✓ selon l’expression :
!2
sin ⇡ a sin ✓
I(✓) = I0
⇡ a sin ✓
La figure de di↵raction par une ouverture circulaire de diamètre D a une symétrie

circulaire et l’intensité dépend de l’angle ✓. L’intensité s’exprime en fonction de l’angle
✓ à l’aide de la fonction de Bessel de premier ordre J1 (x) :
!2
J1 ⇡ D sin ✓
I(✓) = I0
⇡ D sin ✓
Le premier zéro de la fonction de Bessel se trouve en x ' 1.22⇡ : J1 (1.22⇡) ' 0. L’angle
correspondant au premier minimum d’intensité est donc donnée par l’équation :
sin ✓1 = 1.22
D
94
1.0
0.8
0.6
0.4
0.2
-30 -20 -10 10 20 30
1.0
0.8
0.6
0.4
0.2
5 10 15 20 25 30
Figure 11.3 – Di↵raction - Intensité en fonction de l’angle I(✓) avec l’angle en degrés.
En haut, fente rectangulaire atimesb, b , a = 2 en bleu, a = 5 en jaune/orange et
a = 10 en vert. En bas, ouverture circulaire de diamètre D, D = 2 en bleu, D = 5
en jaune/orange et D = 10 en vert.
95

Resume Cours p202 10nov23

Transféré par

Droits d'auteur :

Formats disponibles

Resume Cours p202 10nov23

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Resume Cours p202 10nov23

Transféré par

Droits d'auteur :

Formats disponibles

Table des matières

0.1 Plan du cours . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5

2 L’oscillateur harmonique amorti, oscillations forcées 16

3 Au delà de l’OH simple 30

4 Cordes vibrantes, équation d’onde de D’Alembert 35

5 Réflexion, ondes stationnaires, modes propres d’oscillation 45

6 Ondes longitudinales dans les solides 52

7 Acoustique, propagation du son dans les fluides 59

9 Quelques phénomènes ondulatoires 77

Module de 46 heures, 11 cours (x 2 h = 22 h) et 12 TD (x 2 = 24 h) Plan

0.1 Plan du cours

1.1 Mouvement périodique et mouvement oscillant

1.2 Système masse-ressort

1. Robert Hooke, scientifique anglais du XVII siècle, 1635-1703

x(t) = C cos(!0 t) + S sin(!0 t)

x(t) = A cos(!0 t + ')

A cos(!0 t + ') = A (cos ' cos(!0 t) sin ' sin(!0 t))

r(t) = Cte = R (t) = !0 t + '0

1.4 Pendule simple

sin ✓ ' ✓ pour ✓ ⌧ 1

On retrouve une équation di↵érentielle d’oscillateur harmonique.qLes mouvements de

1.5 Circuit électrique oscillant

On retrouve de nouveau une équation di↵érentielle de type Oscillateur harmonique.

1.6 Résolution de l’équation di↵érentielle de l’OH

mẍ = kx + Fext (t)

la solution générale de l’équation avec second membre s’écrit comme la somme de

On intègre et on obtient l’expression satisfaite par ẋ en fonction d’une première constante

Le cas B = 0 est possible si a < 0. On obtient alors la solution x(t) :

Si on connaı̂t la solution générale z(t) de l’équation l’équation di↵érentielle ż =

L’oscillateur harmonique amorti,

2.1 Forces de frottement fluide

2.1.1 Mouvement amorti, temps d’amortissement

2.2 Equation di↵érentielle de OH amorti (OHA)

2.3 Résolution de l’équation di↵érentielle de l’OHA

x(t) = Ae t ! ẋ(t) = A e t ẍ(t) = A 2 e t

La solution de l’ED dans le régime pseudo-périodique correspond à des oscillations

2.3.2 Recherche de solutions réelles pour le régime pseudo-

Les deux expressions algébriques contenant !p et devant les termes sin(!p t) et

2.4 Solutions de l’OHA en notation complexe et les

x(t) = Re(x̃(t)) x̃(t) = X̃0 e t

= 4!02 ⇠ 2 1 ⇠<1! <0

2.4.1 Régimes apériodique, critique, pseudo-périodique

x(t) = A exp( 1 t) + B exp( 2 t)

• Le régime critique correspond au cas limite où ⇠ = 1 ! = 0 cad = 2⌧1 =

La solution générale de l’équation di↵érentielle s’écrit donc comme le produit

x(t) = exp( !0 t) [A + Bt] (⇠ = 1 ! ⇠!0 = !0 )

x(t) = exp( ⇠!0 t) [C cos(!p t) + S sin(!p t)]

Le pseudo-période Tp = !2⇡p est égale à l’intervalle de temps séparant deux

2.5 Aspects énergétiques, facteur de qualité

2. Facteur de qualité On définit le facteur de qualité Q comme le rapport entre

2.6 Oscillations forcées, résonance

Ou, en utilisant la première notation introduite avec le temps d’amortissement

x̃p (t) = X̃0 ei!t

Ou bien en utilisant la première notation utilisant ⌧ :

Ou bien en utilisant la première notation utilisant ⌧ :

2.6.2 Fréquence de résonance , largeur de bande

Il n’y a pas de résonance si ⇠ > p12 et pour ⇠ ⌧ 1 , on a !res ' !0 (1 ⇠ 2 ),

6. A hautes fréquences, pour ! !0 , le déphasage tend vers ! ⇡ et le gain

On voit que le facteur de qualité Q, défini comme le rapport entre l’énergie