Université du Littoral – Côte d’Opale

Licence des sciences, Calais

Vibrations et ondes
Dmitrií Sadovskií
Département de physique, Université du Littoral, LPCA
e-mail: sadovski@univ-littoral.fr, web: http://pca3.univ-littoral.fr/~sadovski

Calais, automne 2019

Horaires
cours-TD en 12–14 séances de 3h. habituellement 1330 –1645 lundi, salle C104 à partir de 23/9 et jusqu’à . . .
les TP sont suivis par un autre enseignant (Christophe Przygodsky)

Contenu du cours
Ce cours traditionnel est destiné aux étudiants de la troisième année en licence sciences du tronc commun physique–
chimie à Calais. Il a été enseigné auparavant par les Professeurs Hadj Abdelhak (2006–2007) et Robin Bocquet (2008). Le
cours aura les thèmes principales suivantes.
Oscillateur harmonique libre (2 séances)
– Differents exemples (mécanique, électricité).
– Rappels sur la méchanique Newtonienne. Forces.
– Linéarité et non-linéarité.
– Notions de la méchanique Lagrangienne et Hamiltonienne.
– Diagramme de phases, l’espace des phases.
Oscillateur harmonique libre amorti (1 séance)
Oscillations forcées par une force harmonique (1–2 séances)
Oscillateur anharmonique (1 séance) Notions de la théorie des perturbations. NB : ce thème est facultatif, il est addressé
principalement aux étudiants qu’ont choisi l’option Physique et qui font le projet numérique dans le cadre du cours
«Physique numérique» (D. Sadovskií, Ch. Przygodsky)
Partiel sur les oscillations (2h). TD correction immediate (1–2h).
Systèmes linéaires à plusieurs degrés de liberté (2 séances) Solution générale par les utils matriciels. Modes propres.
Cas non-linéaire : l’idée du couplage, résonances. Ce thème sert comme une transitions aux ondes (systèmes au
nombre infinie des degrés de liberté).
Ondes nondispersives (2 séances)
– La corde. Ondes transversales et longitudinales, planes, cylindriques, sphériques.
– Propagation. Interférence et diffraction.
– Ondes sinusoïdales progressives et stationnaires. Noeuds, conditions limite. Résonateurs. Son.
– Equation d’onde (d’Alembert). Equation de corde.
– Terminaison. Réflexion.
Ondes dispersives (1 séance)

1
2 5 décembre 2019 Introduction D. Sadovskií

Modalités de contrôles des connaissances

Contrôle continu Pour pouvoir mieux interagir et démontrer leur niveaux de compréhension (ou incompréhension. . . ),
les étudiants passent rapidement au tableau en TD. En cours les étudiants seront parfois demandés de répondre par
écrit aux 2–3 questions rapides. Il n’y aura pas des notes (sauf si l’on demande exprès), et il ne faut pas avoir peur.
C’est la participation qui compte.
Partiel Il y a un partiel (DS) écrit au mi-cours. Ce partiel couvre le thème «Vibrations». Il est noté sur 20.
L’examen final porte sur tout le cours et est noté sur 20,
la note finale est définie avec le règle du sup contre la moyenne examen+DS.

Livres
– «Oscillations» de Ph. Chen & R. Guillemard, rappels des cours et exercises (existe dans la BU). Voir Chap. 1 et 2.
– «Mécanique» (physique théorique Tome 1) de L.D. Landau & E.Lifchitz, Mir 1988. Pour la partie vibrations : le
plus compacte cours de méca classique.
– «Classical Mechanics» par David Morin (Cambridge University Press, 2008)
– «Ondes mécaniques et sonores» de H. LUmbroso, 2me année , Dunod, Paris, 1999. Voir Chap. 1 et 2 pour la partie
ondes, un peu trop avancé.

Rappels sur le cours «Vibration et ondes» License 3me année. Calais

Version vibondes 2019-12-05 13:09 2019 c Dmitrií A. Sadovskií. Toute dissémination et
reproduction autres que dans le cadre de l’enseignement à l’Université du littoral sont interdites.
La version la plus recente de ce manuel est disponible en format éléctronique (pdf) sur
http://purple.univ-littoral.fr/etudes/vibondes.pdf

Vos suggestions, améliorations et corrections sont bienvenues.

Merci de les envoyer par mél à sadovski@univ-littoral.fr
D. Sadovskií Mécanique 5 décembre 2019 3

Rappels sur la mécanique

La méchanique est historiquement la première théorie dynamique. Dans une telle théorie, on peut trouver l’état du système
à n’importe quel moment du temps (ou plus généralement, pour n’importe quelle valeur de la variable indépendante) si on
connaît son état à un moment du temps donné.
L’évolution du système en temps est décrit par une équation différentielle ordinaire (é.d.o.) ou un système de tels
équations. Soit K le nombre de degrés de libertè, t temps, qk les coordonnées, q̇k = dqk /dt les vitesses, où k = 1, . . . , K.

L’hypothèse centrale de la mécanique : pour caractériser l’état du système à tout mo-

ment du temps donné il suffit de donner sa position q et sa vitesse (drivée première) q̇. Ainsi
l’evolution d’un système mécanique est décrit par des é.d.o.’s d’ordre ≤ 2, c’est à dire les
équations qui contiennent les positions q = (q1 , . . . , qK ) et ses drivées premières q̇ et se-
condes q̈. On appelle éspace des configurations l’espace de dimension K avec coordonnées
qk , et on considère la trajectoire t 7→ (q(t), q̇(t)) du système dans l’espace des phases de
dimension 2K.

Lagrange On écrit les équations de Lagrange pour la fonction de Lagrange (ou Lagran-
gien) L,
d ∂L ∂L
− = 0, k = 1, . . . , K.
dt ∂ q̇k ∂qk Joseph Louis Lagrange
NB : on appelle q est q̇ coordonnées et vitesses généralisées. Pour des systèmes simples avec des masses ponctuelles sans
constraints, on distingue énergie cinétique T (q, q̇) et énergie potentielle V (q), telles que L = T − V

Newton Pour des masses ponctuelles on écrit leurs énergie cinétique

X X
T = 21 mk q̇k2 = 12 mk vk2
k k

On assume aussi que V est une fonction réelle de q(t) seulement. On compute l’accélération

d ∂L
= mk q̈k = mk v̇k = mk ak
dt ∂ q̇k
Isaak Newton
et la force
∂L ∂V (q)
Fk = − =− , k = 1, . . . , N soit F~ = −∇q V (q)
∂qk ∂qk
NB : Newton (1643–1727) formulait sa théorie mécanique bien avant Lagrange (1736–1813) mais l’approche de ce dernier
est la plus générale. La mécanique de Newton (et Kepler) et un cas particulier de la mécanique Lagrangienne.

Hamilton On définit la quantité du movement

∂L
pk :=
∂qk
et la fonction d’Hamilton (ou Hamiltonien)
K
X
H(q, p, t) := q̇k pk − L(q, q̇, t).
k=1

L’évolution du système est donnée par sa trajectoire (q(t), p(t)) qui satisfait K é.d.o’s
d’ordre 1 William Rowan Hamilton
∂H ∂H
q̇k = , ṗk = − , k = 1, . . . , K.
∂pk ∂qk
Notons que dans le cas simple (N particules ponctuelles, K = 3N )

p~n := mn~r˙n := mn~vn et H = T + V.

NB : Cette approche est la plus utile dans le cas sans dissipation où l’énergie est conservée. On connait aussi des applica-
tions en mécanique quantique, où l’énergie est conservée (et quantifiée) et où H devient l’opérateur différentiel Ĥ.
4 5 décembre 2019 Oscillateur harmonique libre, et libre amorti D. Sadovskií

1 Oscillateur harmonique libre

1.1 Cours
– Introduction, présentation et plan du cours
– Rappels sur la méchanique Newtonienne. Forces.
– Pendule mathématique
– Linéarité et non-linéarité.
– Diagramme de phases, l’espace des phases.
– Diagramme de phases complète du pendule

1.2 TD
– Questionnaire éclair de la rentrée (20 min), correction
– oscillateur à ressort posé verticalement

2 Oscillateur harmonique libre, et libre amorti

2.1 Cours
Ressort vertical
– Questionnaire éclair : oscillations libres (15 min)
– Récapitulatif sur séance 1 : équations et leur solutions
– Solution générale matricielle des équations différentielles ordinaires linéaires sur
l’exemple simple d’un oscillateur libre

(a) hyperbolique (b) elliptique (c) focale

F IG . 1 – Portraits locaux dans un espace des phases R2 avec les coordonnées (x, ẋ)

2.2 TD oscillations libres

Exercice 2.1. Lequel des systèmes dont le portrait de phases est presenté sur la figure 1 est stable ? unstable ?

Exemple en circuits électrique LC, RC, RLC Rappels : Les paramètres du système
sont inductivité L résistance ohmique R, et capacité C ; l’état du système est décrit en
C utilisant les charges q(t) et les courants
C
R
L dq
j(t) =
dt
la variable indépendante est le temps t. L’énergie potentielle corresponde à la tension u.
Nous avons
dj q duC 1 dq
uL = −L , uC = et donc = ,
dt C dt C dt
Question : Comparez ce(s) circuit(s) à un système mécanique (ressort, pendule, . . . ). Quels sont les équivalents mécaniques
de R, L, et C ?
Exercice 2.2. Oscillation électrique libre non-amortie (circuit LC). Charge et décharge d’un condensateur. Considérez
le condensateur C branché sur la résistance R comme montré sur la figure ci-à-gauche. Dérivez l’équation différentielle
pour la tension u(t) aux bornes de C. Soit u(0) > 0 au moment t = 0. Trouvez u(t) pour t > 0.
Exercice 2.3. Amortissement (circuit RC). Charge et décharge d’un condensateur. Considérez le condensateur C branché
sur la résistance R comme montré sur la figure ci-à-gauche. Dérivez l’équation différentielle pour la tension u(t) aux
bornes de C. Soit u(0) > 0 au moment t = 0. Trouvez u(t) pour t > 0.
D. Sadovskií Oscillateur harmonique libre, et libre amorti 5 décembre 2019 5

Pendule
Exercice 2.4 (pendule plan). Considerez le pendule mathématique plan à faible amplitude, fig. 2a. En particulier, prenez
l’example du fameux pendule de Foucault, fig. 2b, une balle de 28 kg suspendue au bout d’une corde en acier de 1.5 mm
de diamètre (corde de piano) et de longueur l = 67 m à un point situé au centre de la coupole du Panthéon. Dessinez le
portrait de l’équilibre de ce système, comparez à fig. 1b.
Exercice 2.5 (pendule circulaire). Considerez le pendule mathématique plan à l’amplitude arbitraire, voir la fig. 2c. Quel
est l’espace des phases de ce système ? Combien points d’équilibre y-a-t-il ? Quelle est leur stabilité ? Dessinez le portrait
local de chaque équilibre, comparez aux fig. 1a-1b. Dessinez le portrait global du système dans l’espace des phases.

ϕ y
l
ϕ
mg
mg
x x
(a) (b) (c)

F IG . 2 – Pendule mathématique restreint dans un plan : a à faible amplitude, b pendule de Foucault au Panthéon, c circulaire

Solution On peut choisir les coordonnées cartésiennes (x, y) et polaires (r, ϕ) comme indiqué sur la fig. 2a. A noter que
ϕ = 0 sue l’axe x et ϕ > 0 pour y > 0. Le rayon r ≡ l est fixe et par conséquent, ce système possède un degré de liberté
car les deux coordonnées x et y sont liées par la relation

f (x, y) = x2 + y 2 − l2 = 0 ⇒ 1
2 df /dt = ẋx + ẏy = 0 (2.1)

Autrement, on constate que la position du système peut être définie par une seule coordonnée, l’angle ϕ, et qu’on pourrait
donc utiliser ϕ comme une coordonnée généralisée. On remarque que

x = r cos ϕ et y = r sin ϕ avec r ≡ l = const.

On calcule les dérivées (à savoir faire abseulement !)

ẋ = −l sin ϕϕ̇, ẏ = l cos ϕϕ̇, (2.2)

Les énergies cinétique et potentielle sont, respectivement,

T = 1
2 m(ẋ2 + ẏ 2 ) = 1
2 m l2 ϕ̇2 et V = −mg (x − l) = −mg l (cos ϕ − 1) ≈ 21 mg l ϕ2 (2.3)

Notez que V (x, y) peut être choisi à une constante près ; ici son zero est à l’équilibre du système avec x = l et y = 0.
Pour approximer V (ϕ), on utilise le développement limité de cos ϕ pour |ϕ| ≪ π/2 :

cos ϕ ≈ 1 − 12 ϕ2 + 1 4
24 ϕ + O(ϕ6 ) et sin ϕ ≈ ϕ − 61 ϕ3 + O(ϕ5 ).

En notant qu’en coordonnées polaires, l’énergie totale

H≈ 1
2 m l2 ϕ̇2 + 1
2 mg l ϕ2

est celle d’un oscillateur harmonique en ϕ, nous obtenons la solution. Autrement, on élimine x et ẋ de T + V en utilisant
(2.1) et on simplifie en utilisant le développement par rapport à y/l < ymax /l = ε ≪ 1. Ainsi on obtient
p
T + V = 21 mẏ 2 (1 − y 2 /x2 ) + mg(l − l2 − y 2 ) ≈ 12 mẏ 2 + 21 m(g/l) y 2 .

où on utilise y 2 /x2 ∼ O(y 2√ /l2 ) ≪ 1 et le développement 1 − y 2 /l2 ≈ 1 − 21 y 2 /l2 + O(ε4 ). Notez aussi, que
p

Tmax = Vmax et que ẏmax = gl ymax /l. Par conséquent, le terme ẏ 2 y 2 /x2 ∼ O(ε4 )gl peut être négligé.
6 5 décembre 2019 Oscillateur harmonique libre, et libre amorti D. Sadovskií

Bobbing buoy

Exercice 2.6 (la bouée sautillant). Une bouée flotte sur l’eau calme. Elle consiste d’une large masse m restant entièrement
immergée sous l’eau et d’un tuyau cylindrique vertical attachée à la masse, vidé, enfermé, et sortant d’eau. L’aire de la
section du tuyau est de S [m2 ], la densité de l’eau ρH2 O égale à 1 [kg/L]. On ignore les frottements et la masse de tuyau.
Au moment t = 0, la bouée est légèrement enfoncée dans l’eau de δ par rapport à son position d’équilibre et relâchée.
Trouver la période des oscillations et la solution complète pour le mouvement de la bouée.

Exercice 2.7 (le ballon sautillant). Remplacer la bouée dans l’exercice 2.6 par un ballon de foot de rayon R et masse m
dont la partie submergée à l’équilibre est de profondeur d < R (fig. 3a). Trouver d, la période, et la solution complète.

patm

R
y(0) r
d
d p

densité ρ V0

(a) (b)

F IG . 3 – (a) Le ballon sautillant dans l’eau au moment t = 0, voir exercice 2.7 ; (b) experience de Rüchardt.

Solution On choisit l’axe verticale y avec l’origine à la surface d’eau et la coordonnée y donnant la position de la ligne
de flottaison (waterline). Soit V0 le volume submergé sous l’eau à l’équilibre y = 0. Quand la bouée ou le ballon sortent
de l’eau, donc si y > 0, ce volume et la force d’Archimède diminuent. On somme les forces agissant sur la bouée ou le
ballon et applique la 2me loi de Newton. Pour la bouée cela donne immédiatement
r
ρgS m
r
−mg + ρg V (y) = −mg + ρgV0 − ρgSy = −ρgSy = mÿ, d’où ω0 = et T = 2π .
m ρgS

Dans le cas du ballon, on peut exprimer d en sachant que à l’équilibre la somme des forces est nul
d  
d
Z
−mg + ρgV0 = 0 avec V0 = π R2 − (R − z)2 dz = πd2 R −


,
0 3

mais cela nécessite une analyse des racines d’une équation cubique. Ainsi dans le cas R ≫ 1 on trouve d ≈ m1/2 (πρR)−1/2
et on continue la même racine. Par la suite, on cherche à exprimer le volume submergé en fonction de y
Z 0 Z 0 Z 0

2 
2
V (y) − V0 = S(y) dy = π ρ(y) dy = π R2 − R − (d − y) dy .
y y y

Strictement dit, les oscillations du ballon ne sont pas harmoniques. Dans l’approximation y ≪ d < R on obtient
r
ρg πr2 p
mÿ = ρg V (y) − V0 ≈ −ρgπd (2R − d) y et par conséquent ω0 = avec r = d (2R − d)


m
Exercice 2.8. Considerez l’experience de Rüchardt 1 permettant de mesurer la constante adiabatique κ du gaz dans la
bouteille (fig. 3b) de volume V0 ≫ 41 πd2 ∆y. Expliquez l’analogie et la difference avec l’exercice 2.7. Dans l’hypothèse
d’un procès adiabatique et sans frottements (cf. sec. 8.1.2), exprimez la force de rappel et trouvez la pulsation propre ω0 .

1 voir E. Rüchardt, Eine einfache Methode zur Bestimmung von C /C , Physikalische Zeitschrift, 30 58–59 (1929), et la méthode modifiée de W.
p v
F. Koehler, A Laboratory experiment on the determination of γ for gases by self-sustained oscillations, Am. J. Phys. 19, 113 (1951).
D. Sadovskií Oscillateur harmonique libre, et libre amorti 5 décembre 2019 7

3 Oscillateur harmonique libre amorti

3.1 Cours
– Questionnaire éclair : forme matricielle (15 min)
– Récapitulatif sur la représentation matricielle
– Notions de la méchanique Lagrangienne et Hamiltonienne
– Pertes, introduction de friction
– Equation différentielle et sa solution formelle

3.1.1 L’équation et sa solution générale

On considère l’équation de 2ème ordre de forme plus générale
ω0 ω0
ẍ + ẋ + ω02 x = 0, où Q ∈ R>0 et γ :=
Q 2Q
sont appelés, respectivement, facteur de qualité et constante de l’amortissement. On considère u = (u1 , u2 ) et on assume
que u1 := x et u2 := ẋ. On obtient
u̇1 = u2  
0 1
ω0 ⇒ u̇ = Au avec A =
u̇2 = −ω02 u1 − u2 −ω02 −ω0 /Q
Q
On cherche les valeurs propres λ1,2 de la matrice A.
√ ω0 p
det(A − λ) = λ(λ + ω0 /Q) + ω02 = 0, avec D= 1 − 4Q2
Q
et par conséquent, " #
p
−1 ± 1 − 4Q2
λ1,2 = ω0
2Q

Régime Qualité valeurs propres de A Solution x(t)

1
apériodique 0<Q< 2 λ1,2 ∈ R<0 c1 exp(λ1 t) + c2 exp(λ1 t)
critique Q = 12 λ1 = λ2 = −ω0 ∈ R<0 (c1 + c2 t) exp(−ω0 t)
oscillatoire Q > 21 λ1,2 = −γ ± iω ∈ C, ω ∈ R>0 exp(−γt)(c1 cos ωt + c2 sin ωt)
non-amorti Q=∞ iλ1,2 ∈ R c1 sin ω0 t + c2 cos ω0 t

3.1.2 Régime oscillatoire, quasi-période

Dans la solution générale on assume c1 = 0 et on cherche les maximums de la fonction
x(t) = exp(−γt) sin ωt
La condition pour les points critiques est
où
 
ẋ = exp(−γt) −γ sin ωt + ω cos ωt = 0 ⇒ ω cos ωt − γ sin ωt = 0 ⇒ cos(ωt + φ) = 0, tan φ = γ/ω.
La solution générale (pour les maximums et les minimums) est
π
ωtk + φ = + πk, k = 0, 1, 2, . . .
2
La quasi-période est donnée par la condition
2π
Tk = tk+1 − tk−1 =
ω
Notez (développement limité au voisinage γ ≪ ω0 ) que
ω ≈ ω0 − γ/2 et φ ≈ γ/ω0

3.2 TD
– Obtenez les expressions pour tk et déphasage φ. Ce dernier est appelé aussi retardement. Pourquoi ?
– autres problèmes mécaniques : disc, torsion
– two springs in sequence, in parallel
– molécule diatomique
8 5 décembre 2019 Oscillateur harmonique amorti forcé D. Sadovskií

4 Oscillateur harmonique amorti forcé

4.1 Cours
– Questionnaire éclair : forme de la solution avec differents λ
– Récapitulatif sur oscillateur harmonique libre amorti, diagramme de phase
– Équation différentielle avec la force harmonique, solution
– Résonance
Exercice 4.1 (régime sinusoïdal). Pourquoi le régime forcé avec une force harmonique ∝ sin Ωt est tellement impor-
tant ? Donnez au moins trois raisons pourquoi nous utilisons l’électricité au courant alternatif (AC). Aide : analyzer la
production, la transportation, et la conversion.

4.1.1 Force harmonique

On considère l’équation différentielle ordinaire nonhomogène de 2e ordre
ω0
ẍ + ẋ + ω02 x = ẍ + 2γ ẋ + ω02 x = a0 cos Ωt, (4.1)
Q
Exercice 4.2 (fréquences dangereuses). Une voiture de entre axe L = 5m est conçue pour passer sur un obstacle de
l’hauteur h = 0.2m avec une vitesse entre 60 et 100 km/h. Quelle fréquence externe Ω doit on anticiper pour régler ses
amortisseurs ? Quelle sera l’amplitude a0 ?
On considère u = (u1 , u2 , u3 , u4 ) et on assume que u1 := x, u2 := ẋ, u3 := cos Ωt, et u4 := sin Ωt. On obtient
u̇1 = u2  
ω0 0 1 0 0
u̇2 = −ω02 u1 − u2 + a0 u3 −ω02 −ω0 /Q a0 0 
Q ⇒ u̇ = Au avec A =  
 0 0 0 −Ω
u̇3 = −Ω u4
0 0 Ω 0
u̇4 = Ω u3
On cherche les valeurs propres λ1,2 et λ3,4 de la matrice A. L’équation séculaire est
det(A − λ·14 ) = (λ2 + Ω2 )(λ2 + λω0 /Q + ω02 ) = 0.
Notez que A = diag(A1 , A2 ) + A0 où A1 décrit notre système (libre, sans force), A2 décrite le système externe (force),
et A0 , une matrice nilpotente, donne le couplage entre ces deux systèmes. Notez aussi que
det(A − λ·1) = det(A1 − λ·12 ) det(A2 − λ·12 ).
Les solutions sont " p #
−1 ± 1 − 4Q2
λ1,2 = ω0 , λ3,4 = ±iΩ
2Q

4.1.2 Amplitude & phase en régime stationnaire

Considérez la solution particulière (générale pour t → ∞)
x(t) = a+ exp(iΩt) + a− exp(−iΩt)
correspondant au valeurs propres λ3,4 où a± ne dependent du temps. Pour trouver a± on calcule
ẋ = ia+ Ω exp(iΩt) − ia− Ω exp(−iΩt) et ẍ = −a+ Ω2 exp(iΩt) − a− Ω2 exp(−iΩt),
et les substitue dans l’équation (4.1). En supposant que les coefficients devant exp(±iΩt) doivent être nuls, on obtient
a− = ā+ = 1
2 a0 ā et − aΩ2 + 2aγΩi + ω02 a − a0 = 0,
et par conséquent
1 (ω02 − Ω2 ) − 2γΩi
a= =
(ω02 − Ω2 ) + 2γΩi (ω02 − Ω2 )2 + 4γ 2 Ω2
On trouve l’amplitude et la phase
1 2γΩ
|a| = p et ϕ = arg a = − arctan ,
(ω02 − Ω2 )2 + 4γ 2 Ω2 ω02 − Ω2
respectivement, et exprime la solution
x(t) = 12 a0 |a| exp(i(Ωt + ϕ)) + exp(−i(Ωt + ϕ)) = a0 |a| cos(Ωt + ϕ)


D. Sadovskií Oscillateur harmonique amorti forcé 5 décembre 2019 9

Exercice 4.3. Tracez |a| et ϕ en fonction de δ = Ω/ω0 pour differents valeurs de Q = ω0 /(2γ).
1 1
|a| = p 2 = 2p
(ω0 − Ω2 )2 + ω02 Ω2 /Q2 ω0 (1 − δ 2 )2 + δ 2 /Q2
Trouver Ω pour laquelle |a| devient maximale et, en simplifiant pour le cas Q >> 1, trouvez la valeur max(|a|) corres-
pondant. Pour quelle valeur de δ la valeur de |a| est diminuée à la moitié de sa valeur maximale ?

3 Q=3
Q=2
Q=1
Q=0.5
2.5 Q=0.25

2
Amplitude

1.5

0.5

0
0 0.5 1 1.5 2
Frequence ω/ω0

Solution En solvant d|a|/dδ = 0 pour δ nous obtenons

p 2Q2
δmax = 1 − 1/(2Q2 ) ≈ 1 − 1/(2Q)2 + . . . , et max(|a|)ω02 = p ≈ Q + ...
4Q2 − 1
Notez que le maximum n’existe que pour Q > 2−1/2 . Le déphasage ϕ(δ, Q) est donné par
2γΩ 1
ϕ = − arctan = arctan
ω02 − Ω2 Q(δ 2 − 1)
0
Q = 10
π Q=3
− Q=2
4
Q=1
Déphasage ϕ

π
−
2

3π √
− Q = 1/ 2
4 Q = 0.5
Q = 0.25
−π
0 0.5 1 1.5 2
Rapport des fréquences Ω/ω0
On vérifie que ϕ = 0 pour Ω = 0 (force constante) en se rappelant que dans ce cas la solution stationnaire est tout
simpliment la solution particulière constante avec x decalé par rapport à sa position d’équilibre libre
x(t) = a cos(Ωt + ϕ) = a0 ω0−2 .
10 5 décembre 2019 Oscillateur harmonique amorti forcé D. Sadovskií

4.2 TD : Cirquits électriques en régime sinusoïdal

Filtres CR et RC en régime sinusoïdal. Considérez les deux circuits ci-dessous.
filtre RC (passe bas) filtre CR (passe haut)
u1 = u u 2 = uC u2 = uR

R R
C C
u(t) u(t)
masse u1 = u

Dans chaque cas, les repères u1 dit entrée du filtre et u2 dit sortie du filtre indiquent le branchement des deux voies de
l’oscilloscope (voir les TP de 2ème année). On remarque, que dans les deux cas u1 = u(t), la tension à la sortie de la
source (par ex., générateur des fréquences). Notons aussi, que dans le cas du circuit RC on observe u2 = uC , la tension
aux bornes du condensateur C, tandis que pour CR on observe u2 = uR , la tension aux bornes de la résistance ohmique
R.
Exercice 4.4. Peut on mesurer directement (par ex., avec un oscilloscope) la tension uR (uC ) dans le cas du circuit RC
(CR) ? Pourquoi ? Expliquez le règle des masses.
On considère le cas de la tension d’entrée u(t) = U sin Ωt sinusoïdale purement alternative (U0 = 0) de pulsation Ω,
amplitude U et phase ϕ0 = 0. On utilise la loi des mailles pour les tensions uR et uC aux bornes de la résistance R et le
condensateur C branchés en série, la loi d’Ohm pour R, la relation entre la charge Q et la tension uC , et enfin la définition
du courant j pour obtenir le système des équations suivant
 
Q(t) dQ(t)
u(t) = U sin Ωt, u(t) = uR (t) + uC (t), uR (t) = Rj(t), uC (t) = , j(t) = = Q̇ . (4.2a)
C dt
Exercice 4.5. Peut on visualiser la tension uR (uC ) dans le cas du circuit RC (CR) en utilisant un oscilloscope à deux
voies avec la possibilité d’ajouter les voies ? Comment ? Expliquez la loi des mailles.
On cherche la tension uR (t) pour laquelle on obtient l’équation différentielle
uR
u̇R + − U Ω cos Ωt = 0 (4.2b)
RC
Démonstration. On dérive la première, la deuxième et l’avant dernière équation dans (4.2a)

du duC Q̇
= u̇ = U Ω cos Ωt, u̇ = u̇R + u̇C , = u̇C = .
dt dt C
Par la suite, on résoudre le système des équations linéaires pour se débarrasser de tous les fonctions, telles que uC (t),
u̇C (t), j(t), etc, qui nous n’intéressent pas.
La solution uR (t) de (4.2b) est
s
U 1 1 1
uR (t) = UR sin(Ωt + ϕR ), UR = , ∆Ω = 1+ , tan ϕR = , cos ϕR = . (4.2c)
∆Ω (ΩRC)2 ΩRC ∆Ω

Amplitudes UR2 et UC2

1

0.8
amplitude relative

0.6

0.4

0.2

0
0 0.1 0.2 0.3 0.4 0.5
π arctan(ωRC)

En utilisant cette solution et la loi de mailles dans (4.2a) on obtient

U
uC (t) = UC sin(Ωt + ϕC ), UC = , ϕC = ϕR − 21 π, tan ϕC = −ΩRC. (4.2d)
∆Ω ΩRC
D. Sadovskií Oscillateur harmonique amorti forcé 5 décembre 2019 11

Exercice 4.6. Vérifiez la loi des mailles u = uR + uC pour (4.2c) et (4.2d). Vérifiez que (4.2c) est vraiment la solution de
(4.2b). Expliquez comment on peut trouver les expressions pour ∆Ω , ϕR , et UR à partir de l’équation (4.2b) et la forme
générale UR sin(Ωt + ϕR ) de sa solution.
Exercice 4.7. En utilisant (4.2c) trouvez limΩ=0+ ∆Ω et limΩ=∞ ∆Ω . En utilisant vos résultats et (4.2d) trouvez limΩ=0+ UC (Ω)
et limΩ=∞ UC (Ω). Dessinez les fonctions de gaine UR (Ω)/U et UC (Ω)/U . Quelle est leur valeur dans Ω = (RC)−1 ?
Expliquez pourquoi RC est un filtre passe bas, et CR est un filtre passe haut.
Exercice 4.8. Comparez le déphasage ϕR et ϕC de tensions uR (t) et uC (t). Présentez les sur la diagramme des phases.
Quelle tension est en avance par rapport à u(t), quelle est en retard ?

Filtres RLC série en régime sinusoïdal.

√
Exercice 4.9. Montrez que pour un filtres RLC série en régime sinusoïdal, ω0 = 1/ LC et 2γ = R/L.
Exercice 4.10. On construit un récepteur radio. Proposez un schéma et expliquez le principe. S’agit-il d’un système
(oscillateur) libre ou forcé ? Quel régime (apériodique ou oscillatoire) doit on chercher à atteindre ? Pourquoi ? En recevant
en même temps les deux stations radio de fréquences f1 et f2 et d’intensités égales I1 = I2 , on considère que l’une est
distinguée de l’autre si le rapport des amplitudes respectives A1 et A2 des signaux de la réponse du récepteur est supérieur
à 2:1. Quelle est ∆f = |f1 − f2 | minimale pour laquelle votre récepteur permettra de distinguer ces stations ?

Filtres RLC parallel en régime sinusoïdal.

√
Exercice 4.11. Montrez que ω0 = 1/ LC et que 2γ = 1/(RC).

4.3 TD : prepa DS «Vibrations»

Exercice 4.12. Une masse ponctuelle m attachée au ressort de rigidité k obtient une vitesse initiale v0 ; la déformation ini-
tiale du ressort est x0 . Trouvez l’amplitude des vibrations. Trouvez la déformation dans t = T . Prenez valeurs numériques
m = 100g, k = 10N/m, v0 = 1m/sec, et T = 1sec.
Exercice 4.13. Une pendule de longueur l = 5 m fait des petites oscillations. Apres 100 périodes, son amplitude s’est
diminuée par deux fois. Trouvez γ et qualité Q.
Exercice 4.14. Appareillage scientifique est placé sur une platforme amortisseur des vibrations posée au plancher. Le
plancher vibre avec fréquence f = 20 Hz et amplitude A. La fréquence propre de l’appareil sur l’amortisseur est f0 =
20 Hz, l’amplitude de ses vibrations est a. Trouvez la valeur maximale du rapport de A et a.
Exercice 4.15. Soit x la coordonnée d’un oscillateur harmonique non-amorti. Trouvez l’expression pour x(t) si la force
extérieure est F (t) = 0 pour t 6∈ [0, T ], F (t) = F0 pour t ∈ [0, T ]. Donnez la condition pour retrouver l’oscillateur à son
équilibre ẋ = x = 0 pour t > T .
Exercice 4.16. Trouvez l’amplitude et la phase d’un oscillateur harmonique non-amorti étant à son équilibre à t = 0. La
force extérieure est F (t) = 0 pour t 6∈ [0, T ], F (t) = ξt pour t ∈ [0, T ].
Exercice 4.17. Considérez la solution particulière de l’équation (4.1) sous forme
x(t) = A cos(Ωt + ϕ) (4.3)
et retrouvez les expressions pour l’amplitude A et le déphasage ϕ obtenus dans la section 4.1.2.

Solution de l’exercice 4.17 On compute les dérivées de x(t)

ẋ = −AΩ sin(Ωt + ϕ) et ẍ = −Ωx,
et les substitue dans l’équation (4.1) dont les termes se décomposent en utilisant
cos(Ωt + ϕ) = cos(Ωt) cos ϕ − sin(Ωt) sin ϕ et sin(Ωt + ϕ) = cos(Ωt) sin ϕ + sin(Ωt) cos ϕ.
En rassemblant les termes devant cos Ωt et sin Ωt le résultat s’écrit sous la forme
Ccos cos Ωt + Csin sin Ωt = 0.
On conclut que la solution (pour tout t) n’est possible que dans le cas Ccos = Csin ≡ 0. Par conséquent
(ω0 − Ω) cos ϕ − 2γΩ sin ϕ = X cos(α + ϕ) = a−1 avec a = A/a0 > 0,
1/2
avec X = (ω02 − Ω2 )2 + 4γ 2 Ω2

(ω0 − Ω) sin ϕ + 2γΩ cos ϕ = X sin(α + ϕ) = 0 > 0.
La solution unique de ces équations est possible si sin(α + ϕ) = 0 et cos(α + ϕ) = 1, donc α + ϕ = 0.
12 5 décembre 2019 Oscillateur harmonique amorti forcé D. Sadovskií

4.3.1 Pendule du Pohl

Considérez le TP qu’on fait parfois en mécanique (voir la photo) et en capes.
On trouve ce système dans les montres mécaniques. Ce pendule est contitué
d’un disque dont la rotation libre autour de son centre est empêchée par un
ressort spiral qui tend à ramener le disque vers sa position d’équilibre. Un
petit rouleau de rayon r, entrainé par un moteur à une fréquence de rotation
réglable Ω est relié au ressort spiral, lui donnant un déplacement angulaire
qui force les oscillations du disque. Le disque peut être freiné pour contrôler
l’effet d’amortissement. Il s’agit d’un frein électromagnétique sans contacte
exploitant les courants de Foucault. Coefficient de freinage est λ. Deux bras
relient le ressort au rouleau. Le bras quasi-horizontal du longueur L est atta-
ché au rouleau, celui attaché au ressort reste quasi-vertical et son longueur l
Pendule de Pohl
est d’ordre L. Pour assurer que les deformations du ressort sont petites, nous
devons garder r ≪ l ∼ L.

Exercice 4.18 (Pendule du Pohl). Soit ϕ et ϑ sont les angles donnant, respectivement, les positions du disque de pendule
et du rouleau de moteur par rapport à leurs centres respectifs O et O′ . Si θ(t) = Ωt, trouvez ϕ(t). Vérifiez que ϕ ≪ π.
Pour simplifier vos calculs, considérez le cas particulier de xO − xO′ = L et yO − yO′ = l (voir le dessin).

Solution de l’exercice 4.18 Plaçons l’origine des coordonnées dans O′ , prenons Q = (X, Y ) et q = (x, y) comme
coordonnées de deux bouts du bras horizontal (positions de charnières). On trouve

Q = (xO , yO ) + l(sin ϕ, − cos ϕ) et q = r(cos ϑ, sin ϑ)

La distance entre ces points est fixe et est égale à L, on impose donc la condition

kQ − qk2 − L2 = (xO + l sin ϕ − r cos ϑ)2 + (yO − l cos ϕ − r sin ϑ)2 = 0

donnant

yO 2 + xO 2 + r2 + l2 − L2 + (−2 r sin ϑ − 2 l cos ϕ) yO + (2 l sin ϕ − 2 r cos ϑ) xO + 2 l r sin(ϑ − ϕ) = 0.

On simplifie cette condition pour le cas particulier en supposant Y0 = l, X0 = L, et r = εL

2 l2 (1 − cos ϕ) + 2 L l sin ϕ + ε L (−2 L cos ϑ + 2 l (− sin ϑ + sin(θ − ϕ))) + ε2 L2 = 0,

et apres le développement limité en ϕ au voisinage ϕ = 0 on obtient l’équation

+2 l L ϕ − 2 ε L (L − l ϕ) cos ϑ + ε2 L2 + O(ϕ2 ) = 0.

(à cet ordre il suffit de substituer cos ϕ ≈ 1 et sin ϕ ≈ ϕ) dont la solution est

L cos ϑ − ε/2 L ε 
ϕ=ε ≈ε cos ϑ + cos 2ϑ + . . . .
l 1 − ε cos ϑ l 2
On conclut qu’à son ordre principal ε, la petite déformation forcée du ressort
r
cos Ωt + O ε2 où r = εL


ϕ(t) =
l
est harmonique et proportionnelle au rapport r/l ≪ 1.
Exercice 4.19 (Pendule du Pohl). En appliquant la mécanique Newtonienne à un corps rigide en rotation autour de son
axe O, obtenez l’équation de mouvement du pendule dans l’exercice 4.18 et donnez sa solution stationnaire à t → ∞.
D. Sadovskií Vibrations nonlinéaires 5 décembre 2019 13

5 Vibrations nonlinéaires à un seul degré de liberté

Rappel : degré de liberté est le nombre de coordonnées généralisées nécessaires pour caractériser de façon non-
ambiguë tout les configurations du système.

Exercice 5.1 (degrés de liberté vibrationnels). Considérez molécules (a) diatomique, (b) linéaire à 3 atomes comme
CO2 ou CH2 CH2 , (c) nonlinéaire comme H2 O, C6 H6 , etc. Dans chaque cas, quel est le nombre des degrés de liberté
vibrationnels ? Proposez une formule générale en fonction de nombre des atomes.

Rappel : développement limité Un outil indispensable de simplification et de l’analyse local, e.g., au voisinage d’un
équilibre stable (elliptique) du pendule. Dans sa forme la plus simple, il s’agit de la formule de Newton–Leibnitz :

f ′ (x) = df (x)/dx := lim (f (x) − f (x0 )) /(x − x0 ) ⇒ f (x) ≈ f (x0 ) + f ′ (x0 )(x − x0 ) + O (x − x0 )2 . (5.1)


x→x0

Plus généralement, les approximations successives en puissances de (x − x0 ) sont données par la série de Taylor

1 df 
f (x) ≈ f (x0 ) + f ′ (x0 )(x − x0 ) + 12 f ′′ (x0 )(x − x0 )2 + . . . + (x − x0 )n + O (x − x0 )n+1 . (5.2)


n! dx x0


1
En particulier, si x0 est un point d’équilibre, au voisinage de x0 nous avons f (x) ≈ f (x0 ) + 2 f ′′ (x0 )(x − x0 )2 .

5.1 Molécule diatomique, potentiels empiriques

On caractérise les molécules diatomiques2 par leur distance interatomique d’équilibre re et leur énergie de dissociation
De , voir quelques exemples ci-dessous3

molécule état De (eV) re (Å) ν0 (cm−1 ) molécule état De (eV) re (Å) ν0 (cm−1 )
LiH Σ 2.47 1.596 1406 H2 Σg 4.5 0.741 4401
CO Σ 11.2 1.128 2170 N2 Σg 9.8 1.098 2359
HCl Σ 4.4 1.275 2991 O2 3
Σg 5.2 1.208 1580
HF Σ 5.8 0.917 4138 F2 Σg 1.6 1.412 917
NO Σ 7.0 1.151 1904 Cl2 Σg 2.5 1.988 560

L’énergie totale ou l’hamiltonien est (ici µ est la masse réduite, voir problème de deux corps en cours de mécanique)

p2 m1 m 2
H= + V (x), où x = r − re , et µ = .
2µ m1 + m2

On représente V (x) par différentes fonctions empiriques, parmi lesquelles le potentiel de Morse4

2
VM (x) = De 1 − e−α(x/re ) , (5.3)

ou le potentiel de Lennard-Jones5 dit «6–12»

  
σ 12  σ 6 a b
VLJ (r) = 4ǫ − = 12 − 6 , avec r = re + x, (5.4)
r r r r

ou encore le potentiel de Kratzer

a2 a1
−
VK (r) = + a0 , avec r = re + x. (5.5)
r2 r
sont les plus connus. Enfin, dans l’approximation harmonique, on utilise le potentiel quadratique

k 2
V0 (x) = x , k = d2 V /dx2 dans x = 0, (5.6)
2
où k connu aussi comme k0 ou kxx est un appelé constante de force harmonique.
2 Voirles cours de la mécanique classique et quantique
3 Huber, K. P. & Herzberg, G. (1979) Molecular Spectra and Molecular Structure IV. Constants of Diatomic Molecules (New York, Van Nostrand,
Reinhold) ; http://hyperphysics.phy-astr.gsu.edu/hbase/tables/diatomic.html ; Computational Chemistry Comparison and
Benchmark Database. Part XIII – Vibrations (Aug. 2005) NIST Standard Reference Database 101, Release 12, http://srdata.nist.gov/
cccbdb/vibrations.asp. Notez qu’en anglais re est appelé equilibrium separation ou bond length
4 Morse, P. M. (1929) Diatomic molecules according to the wave mechanics. II. Vibrational levels. Phys. Rev. 34, 57–64.
5 Lennard-Jones, J. E. (1931) Cohesion. Proc. Phys. Soc. 43, 461-482.
14 5 décembre 2019 Vibrations nonlinéaires D. Sadovskií

En mécanique classique, on obtiens les équations hamiltoniennes du mouvement, un système des équations différen-
tielles ordinaires d’ordre 1, pour les variables canoniques conjuguées x et p
dx ∂H p dp ∂H dV
= ẋ = = , = ṗ = − =− = −F (x).
dt ∂p µ dt ∂x dx
Pour passer en mécanique classique newtonienne, on observe que
ṗ 1
ẍ = = − F (x),
µ µ
et on considère donc F (x) comme une force agissant sur la particule de masse µ.

Unités. On utilise le système des unités atomiques (u.a.) dans lequel ~ = 1. L’u.a. d’énergie Eh , nommé Hartree, est
égale à 2625.5 kJ/mol, 627.5 kcal/mol, 27.211 eV, ou 219474.6 cm−1 . Le cm−1 est utilisé souvent en spectroscopie pour
exprimer les fréquences de vibration et/ou les énergies. Ainsi, pour obtenir la fréquence en Hz on multiplie par c en cm.
L’u.a. de longueur a0 , ou rayon de Bohr, est égale 0.529 177 210(2) Å. Les u.a. de masse et de charge sont, respectivement,
la masse me et la charge e d’électron.
Exercice 5.2 (Choix de la fonction potentielle, calcul des paramètres). Trouvez les relations entre caractéristiques mo-
léculaires re , De , et ω0 , et les paramètres des fonctions potentielles VM , VLJ , et/ou VK . Peut on reproduire correctement
les trois caractéristiques ? Faites le calcul pour une molécule de votre choix. Pour les petits vibrations autour re , donnez
une approximation polynomiale ; déterminez les constants de force harmonique kxx et cubique kxxx . Faites le dessin de
différents fonctions potentielles et des approximations quadratique et cubique. Dessinez les trajectoires du système dans
l’espace de phases (x, ẋ) comme les niveaux des valeurs constantes de H pour V (x) et ses différentes approximations.
Expliquez les différences. Que pouvez vous dire sur la vitesse d’avancement du système le long de la trajectoire (la vitesse
de phase) ? Donnez un graphe qualitative de t 7→ x(t) pour V0 (x) et V (x). Commentez la différence.

5.2 Pendule plan à grande amplitude

Exercice 5.3. Retournez au pendule dans l’exercice 2.5 et prenez m = 1 et ℓ = g. Donnez les termes nonlinéaires
corrigeant l’approximation harmonique (exercice 2.4) pour son énergie et son équation de mouvement, cf. fig. 4.

2.5

1.5
vitesse angulaire dϕ/dt ∈ R [rad/sec]

0.5

−0.5

−1

−1.5

−2

−2.5

−1 −0.75 −0.5 −0.25 0 0.25 0.5 0.75 1

angle du pendule ϕ ∈ S1 [π rad]
F IG . 4 – Les trajectoires du pendule (fig. 2c) dans son espace des phases P = S1 ×R ; la région en gris clair corresponde aux vibrations.
D. Sadovskií Vibrations à plusieurs degrés de liberté 5 décembre 2019 15

6 Vibrations à plusieurs degrés de liberté

Exercice 6.1 (une chaîne de masses couplées). Considérez la chaîne dans la fig. 5 de n points matériels de masses mi ,
i = 1, . . . , n, couplées par des ressorts de rigidités kj , j = 1, . . . , n + 1. Les points peuvent se déplacer sans frottement
le long de l’axe horizontal (axe x). A priori, les différentes masses et ressorts ne sont pas identiques. Quel est le nombre
des degrés de liberté ? Proposez les coordonnées dans l’espace des configurations. Dessinez les forces agissant sur mi .

k1 k2 k3 kn kn+1
m1 m2 m3

mn−1 mn

F IG . 5 – Une chaîne de points matériels dans R1 couplés par des ressorts.

Exercice 6.2 (molécule triatomique linéaire non-pliant). Quelles sont les vibrations de la molécule CO2 ? Afin de modéli-
ser ses vibrations d’élongation, considérez la chaîne de trois points matériels dans R1 avec m1 = m3 = mO et m2 = mC
couplés par deux ressorts de rigidité k1 = k2 =: k. Utilisez la mécanique newtonienne ou hamiltonienne pour obtenir les
équations de mouvement. Ces équations, sont elles séparables ? Donnez leur solution générale.

Solution de l’exercice 6.2 Soit x1 , x2 , et x3 les coordonnées des atomes O, C, et O de O−C−O, et p1 , p2 , et p3 les
impulsions (quantités de mouvement) correspondant (conjugués). L’énergie totale
3 3
X pi k  X pi k 2
(r1 − re )2 + (r2 − re )2 = q1 + q22 ,

H= + +
i=1
2mi 2 i=1
2mi 2

où re est la distance C−O à l’équilibre, et r1 = x2 − x1 et r2 = x3 − x2 sont les deux distances C−O instantanées, est
conservée. Nous pouvons utiliser les équations de mouvement sous forme hamiltonienne6 :
∂H p1 ∂H
ẋ1 = = , ṗ1 = − = kq1 , ⇒ mO ẍ1 = kq1 , (6.1a)
∂p1 mO ∂x1
∂H p2 ∂H
ẋ2 = = , ṗ2 = − = −kq1 + kq2 , ⇒ mC ẍ2 = −k(q1 − q2 ), (6.1b)
∂p2 mC ∂x2
∂H p3 ∂H
ẋ3 = = , ṗ3 = − = −kq2 . ⇒ mO ẍ3 = −kq2 . (6.1c)
∂p3 mO ∂x3
Notons immédiatement que, conformément au principe de Galilée,
mO ẍ1 + mC ẍ2 + mO ẍ3 = 0.
C’est à dire, le centre des masses du système donné par mO x1 + mC x2 + mO x3 n’accélère pas. Les deux autres combi-
naisons indépendantes de (x1 , x2 , x3 ) représentent les mouvements internes, donc les deux vibrations d’élongation. Ainsi,
pour les combinaisons
q+ := q1 − q2 = r1 − r2 = 2x2 − (x1 + x3 ) et q− := q1 + q2 = r1 + r2 − 2re = −(x1 − x3 ) − 2re (6.2)
nous obtenons les deux équations différentielles séparées (ou découplées)
k 2mO + mC
q̈+ + q+ = 0 et q̈− + k q− = 0. (6.3)
mO mO m C
On appelle de telles combinaisons les modes normaux, ou modes propres de vibration.
Exercice 6.3 (principe de symétrie). Quelle(s) est(sont) la(les) symétrie(s) du système dans l’exercice 6.2 (molécule
O−C−O libre non-pliant) ? Quelles sont les propriétés de symétrie de q± dans l’éq. (6.2) ? Y-a-t-il une relation avec la
séparabilité des équations (6.3) ?
Exercice 6.4. Considérez la chaîne dans la fig. 5 de n points matériels de masses mi ≡ m, i = 1, . . . , n, couplées par des
ressorts de rigidités kj ≡ k, j = 1, . . . , n + 1. Donnez les équations de mouvement et proposez leur solution.
Exercice 6.5. Pour le système dans l’exercice 6.2, remplacez le potentiel harmonique par le potentiel de Morse VM (r)
dans l’éq. (5.3) et trouvez le terme de couplage principal (cubique) entre q+ et q− . (NB : utiliser développent limité.)

6 Si non, faites le bilan des forces pour chaque atome et appliquez la 2ème loi de Newton. C’est encore plus simple.
16 5 décembre 2019 Ondes nondispersives D. Sadovskií

7 Ondes nondispersives (séance 1)

Exemples La corde. Les vagues . . .

Ondes scalaires
1. Onde plane et non-déformable. Fonction d’onde f (at ∓ bx).
2. Cas particulier important : onde sinusoïdale.
3. Types géometriques : plane, circulaire, cylindrique, sphérique
4. Ondes transversales et longitudinales.
Caractéristiques & paramètres (résumé) :
1. Période T , fréquence ν, et longeur d’onde λ.
2. Vitesse de propagation (vitesse de phase) v, nombre d’onde k ou vecteur d’onde ~k.
Equation d’onde (equation d’Alembert)

7.1 Cours
Procès ondulaire dans R : Soit s une caractéristique physique (dépression, déplacement,
distortion, etc), t temps, x coordonnée le long l’axe x. Considérez la fonction

s = f (φ) = f (t − x/v) = f (t − kx x) avec v > 0 (7.1a)

Une telle fonction décrit un procès ondulaire. Elle est appelée ainsi la fonction d’onde. Son
argument φ = t−x/v est appelé phase. Au present, considerez un observateur qui se déplace
avec une vitesse v, donc sa coordonnée est x(t) = x0 + v t. Pour un tel observateur la phase
restera constante :
s̃ = f t − x(t)/v = f (−x0 /v) = const.

Alors la fonction d’onde (7.1a) décrit une propagation en direction positive. De façone simi-
laire on obtient que
s = f (t + x/v) avec v > 0 (7.1b)
correspond à une propagation en direction negative. La constante v > 0 donne la vitesse de Christiaan Huygens
propagation. On peut ainsi introduire la forme plus générale d’une fonction d’onde

s = f (ωt − kx x) = f ω(t ∓ x/v) avec v > 0 et kx = ±ω/v, (7.1c)

où kx 6= 0 dit nombre d’onde est un nombre positif ou negatif selon la direction de propagation du procès, tandis que ω
donne l’échelle du temps.

Exercice 7.1. Considerez le procès décrit par la fonction d’onde f ω(t ∓ x/v) où

f (φ) = exp(−φ2 /2) (la «cloche» Gausienne).

Faites le graphes pour différentes valeurs de t. Que pouvez vous dire sur la direction de propagation ?

Procès ondulaire dans R3 : Les positions en espace euclidienne physique à trois dimensions sont données par le vecteur
~r ; en coordonnées cartésiennes ~r = (x, y, z)T . On obtient

s = f (ωt − ~k · ~r) = f (ωt − kx x − ky y − kz z) avec k = k~kk > 0 et v = ω/k (7.2)

où le vecteur d’onde ~k pointe la direction de propagation.

Equation d’onde (Euler 1748, d’Alambert 1747, connue comme équation de d’Alembert)

df df
Pour la dérivée complète de f on utilisera la notation suivante : f ′ (φ) = =
dφ d(t ∓ x/v)
D. Sadovskií Ondes nondispersives 5 décembre 2019 17

Dérivons s = f (t ∓ x/v) dans (7.1) deux fois sur x et sur t. Cela donne, respectivement,
∂2s 1 ∂2s
= 2 f ′′ et = f ′′
∂x2 v ∂t2
(Notez que la double dérivation nous débarrasse du signe ∓ devant x/v.) Par consequent
∂2s 2
2∂ s
= v (7.3)
∂t2 ∂x2
De façon similaire on obtient poue les ondes en éspace physique R3 avec coordonnées (x, y, z)
∂2s
 2
∂2s ∂2s

2 ∂ s
= v + + = v 2 ∆s. (7.4)
∂t2 ∂x2 ∂y 2 ∂z 2
où ∆ est l’opérateur de Laplace ou Laplacien.
Exercice 7.2. Faites en détailles le calcul pour obtenir (7.4) à partir de (7.2).

Ondes sinusoïdales dans R Dans ce cas très important, la fonction f (φ(t, x)) dans (7.1c) est une harmonique
s = A cos(ωt − kx x + φ0 ), (7.5)
avec amplitude A > 0 (constante réelle positive), pulsation ω, nombre d’onde kx , et—s’il le faut—la phase «initiale» φ0
dans le point x = 0 au moment t = 0. On utilise également la représentation complexe de (7.5) basée sur exp(ωt − kx x).
Exercice 7.3. Vérifiez que la fonction d’onde (7.5) satisfait l’équation de d’Alembert (7.3).
A chaque moment du temps, l’onde se présente dans l’espace R comme une fonction harmonique A cos(kx x) sans fins de
période spatiale
λ = 2π/kx (7.6)
dit longeur d’onde. On voit que la vitesse de propagation (ou la vitesse de phase) est
λ
v = ω/kx = = λν, (7.7)
T
où T = 2π/ω [sec] est la période temporelle des oscillations, et ν = ω/(2π) [Hz] est leur fréquence.

Ondes sinusoïdales dans R2 et R3 On denote les coordonnées cartésiennes, cylindriques, et sphériques, respectivement,
(x, y, z), (ϕ, ρ, z), et (ϕ, θ, r). Pour différents symétries agissant dans l’espace on obtient :
Ondes planes
s = A cos(ωt − ~k · ~r + φ0 ), (7.8a)
Ondes circulaires et cylindriques
A
s = √ cos(ωt − kρ ρ + φ0 ), (7.8b)
ρ
Ondes sphériques
A
s= cos(ωt − kr r + φ0 ), (7.8c)
r
Exercice 7.4. Pour chaque fonction d’onde (7.8) et pour un moment fixe du temps, prenez un point dans l’espace et
dessinez l’ensemble de tout points dans l’espace dans lequels la phase est de même valeur. Cet ensemble est appelé le
front de l’onde. Dans chaque cas, quelle est sa symétrie ?
Exercice 7.5. Pour chaque fonction d’onde (7.8), vérifiez que l’équation de d’Alembert (7.3) est satisfait. Notez que en
coordonnées cylindriques et sphériques le Laplacien s’exprime, respectivement,
1 ∂2f ∂2f ∂2f 1 ∂2f ∂2f
 
1 ∂ ∂f 1 ∂f
∆f = ρ + 2 2+ 2 = + + + , (7.9a)
ρ ∂ρ ∂ρ ρ ∂φ ∂z ∂ρ2 ρ ∂ρ ρ2 ∂φ2 ∂z 2
∂2f
   
1 ∂ ∂f 1 ∂ ∂f 1
∆f = 2 r2 + 2 sin θ + 2 2 . (7.9b)
r ∂r ∂r r sin θ ∂θ ∂θ r sin θ ∂φ2
Solution : Pour le cas d’onde sphérique, f ne depend que de r. Par conséquant, nous avons à vérifier que
∂2s v2 ∂ v2 ∂
 
2 2 2 ∂s
= v ∆s ⇒ −ω s = r = (− cos ξ + kr sin ξ) = −v 2 k 2
∂t2 r2 ∂r ∂r r2 ∂r
où ξ(t, r) := (ωt − kr + φ0 ).
18 5 décembre 2019 Ondes nondispersives D. Sadovskií

7.2 Ondes électromagnétiques

Exercice 7.6. Rappel sur l’equation d’onde (sec. 7.1). Montrez que tout un processus de propagation ondulaire d’une
grandeur s le long de l’axe x est décrit par une fonction du type s = s (t ∓ x/v). Ainsi donnez la définition de la phase,
de la vitesse de phase, et du front, et précisez la direction de propagation. En comparant les deuxièmes dérivées spatiale
et temporelle de s, obtenez l’équaion d’onde dans une dimension et la généralisez pur 3D.

7.2.1 Rappel sur la théorie de Maxwell

Coulomb 1785, Oersted 1820, Faraday 1831, Maxwell 1861–1864 . . .
Equations de Gauss (avec ρ : R3 → R la densité de charge) :
div D = ρ et div B = 0, où B = µµ0 H et D = εε0 E. (7.10a)
Equations de l’induction (avec j : R3 → R3 la densité du courant de conductivité) :
∂D ∂B
rot H = j + et rot E = − (7.10b)
∂t ∂t
Exercice 7.7. Enoncez les équations de Maxwell en forme locale (7.10) dans le vide (donc
en absence des courants et des charges électriques). Appliquez l’identité du calcul véctoriel James Clerk Maxwell
2
∇ × (∇ × A) = ∇(∇ · A) − ∇ A, avec div A := ∇ · A, et rot A := ∇ × A
√
pour en déduire que la vitesse de la lumière (dans le vide) est c = 1/ µ0 ε0 = 2.99792458 108 m/sec, avec µ0 = 4π 10−7
H/m, et ε0 = 1/(µ0 c2 ) F/m.
This velocity is so nearly that of light, that it seems we have strong reasons to conclude that light itself
(including radiant heat, and other radiations if any) is an electromagnetic disturbance in the form of waves
propagated through the electromagnetic field according to electromagnetic laws. James Clerk Maxwell, 1865

7.2.2 Equation d’onde électromagnétigue plane

Exercice 7.8. Considérez les champs B(r, t) = B0 cos(ωt + kz) ex et E(r, t) = E0 cos(ωt + kz) ey . Conforment ils aux
équations de Maxwell (7.10) ? Sous quelles conditions ? Appliquez (7.10) et précisez la relation entre B0 et E0 , ω et k.
Exercice 7.9. Considérez une onde plane dont le front est parallèle au plan {z = 0} (plan x0y).
i. Quelles conditions sur les dérivées spatiales doit-on imposer dans ce cas ?
ii. Simplifiez les équations du groupe induction (7.10b) et montrez que Bz et Ez sont constants (statiques). En
appliquant les équations du groupe Gauss (7.10a) montrez que Bz et Ez ne dependent de z. Pour la suite mettez
Bz = Ez = 0.
iii. Parmi les équations restantes du groupe induction (7.10b), choisissez une paire décrivant une onde plane polarisée
rectilignement selon l’axe x. Que impose cela sur les composantes de E ? Quelles composantes sont présentes dans
cette paire des équations ? Que pouvez vous conclure sur l’orientation des champs électriques et magnétiques dans
cet onde ?
iv. Dérivez ces équations soit par rapport au temps soit par rapport à la variable spatiale restante (laquelle ?) et
obtenez ainsi les équations d’onde pour la composante E(z, t).
v. Quelle est la vitesse de propagation ? Quelle est la forme générale de la fonction Ex (z, t) ?
vi. Obtenez aussi l’équation d’onde pour la composante du champ magnétique H(z, t) (laquelle ?) induite par la
variation de E(z, t). Obtenez la relation entre E(z, t) et H(z, t). Ces deux champs, sont-ils en phase ?
vii. Exprimez E(r, t) pour une onde plane, harmonique, polarisée rectilignement, dont la diréction de propagation
est donnée par le vecteur d’onde k. Dessinez le trièdre direct (k, E, B), déduisez la relation entre E(z, t) et H(z, t).

7.2.3 Énergie des ondes électromagnétigues. Intensité.

Exercice 7.10. En sachant que l’énergie per unité du volume dW/dτ des champs électrique et magnétique est respecti-
vement 12 εε0 E2 et 21 µµ0 H2 trouvez l’énergie dans un volume parcouru par l’onde pendant le temps dt. Ainsi trouvez le
flux d’énergie et donnez l’expression du vecteur de Poynting P. Expliquez sa relation avec l’énergie du champ électroma-
gnétique (dans le vide) et sa signification physique (comme vecteur).
Exercice 7.11. Soit E = (Ex , 0, 0), où Ex = E0 sin(ωt − kx) et k = 2π/λ, se propage dans le vide. Trouvez la grandeur
instantanée P (x, t) et l’orientation du vecteur de Poynting P. Trouvez la valeur moyenne P appelée intensité.
Exercice 7.12. Onde stationnaire. Considérez deux ondes E1 et E2 de mêmes pulsation ω, norme kkk du vecteur d’onde
et amplitude E0 qui se propagent dans les directions opposées selon l’axe x (ex et − ex ). Trouvez le champ résultant
E. Representez graphiquement E(x, t) pour t = 0, T /4, T /2, . . . (où T est la période). Utilisez une des équationas de
Maxwell pour retrouver B(x, t). Enfin trouvez P (x, t) et P .
D. Sadovskií Ondes nondispersives 5 décembre 2019 19

7.3 L’équation de Schrödinger, les ondes de de Broglie, et l’atome d’hydrogène

L’équation de Schrödinger7 , connue aussi sous le nom de l’équation ondulaire (wave equa-
tion, Wellengleichung),
∂
i~ ψ(t, x) = Ĥ ψ(t, x), (7.11)
∂t
où l’opérateur différentiel Ĥ représente l’hamiltonien (l’énergie) du système, est l’équa-
tion fondamentale de la mécanique quantique ressemblant à (7.3). Sa solution ψ(t, x) (ici
en version spatiale) décrit la probabilité de trouver la particule dans le point x au temps t.
Remarquablement, ψ(t, x) est couramment nommée fonction ondulaire (wavefunction, Wel-
lenfunktion) ou fonction d’onde. Pour une particule libre de masse m, l’hamiltonien est égal
à l’énergie cinétique
p̂2 1 ∂2
Ĥ = = −~2 ,
2m 2m ∂x2 Erwin Schrödinger
et les solutions de l’équation ondulaire (7.11)

ψ(t, x) = A exp(i ω t − i k x) (7.12)

ont la forme d’une onde plane complexe progressive dont l’amplitude A est déterminée par la normalisation de la dis-
tribution de probabilité8 . Les constantes d’intégration ω ∈ R≥0 et k ∈ R correspondent à l’énergie et à la quantité de
mouvement classiques, respectivement
~2 k 2
E = ~ω = et p = ~ k,
2m
mais la célérité ω/k de (7.12) ne donne que la moitié v/2 de la vitesse classique v. Ceci est dû au fait que (7.12) est
delocalisée sur tout R et ne corresponde pas directement à une particule classique.
Considérons le système de deux corps, le proton et l’électron, constituant l’atome
~v
d’hydrogène. La masse du proton de 1836me est assez grande pour assumer que le proton
~v + d~v reste immobile dans 0 ∈ R3 . Dans le champ central Coulombien du proton, l’électron
rn
dϕ classique reste sur des orbites circulaires (ou, plus généralement, elliptiques) de rayon r.
Si, comme suggéré par de Broglie 9 , on représente l’électron par l’onde (7.12) de longueur
λ = 2π/k, il est impératif (voir la figure ci à gauche) de satisfaire la condition
2π
2πrn = nλ = n~ où n = 1, 2, 3, . . . , (7.13)
mv
pour une onde stationnaire le long d’un cercle de rayon r = rn , car tout une autre onde
se auto-détruit (c’est à dire par l’interférence avec elle même10 ). D’autre part, un mouve-
ment classique à vitesse de valeur absolue constante v le long de ce cercle subit une accélération radiale

d~v v dϕ v2
ar = = =
dt dt rn
due à la force colombienne de l’attraction entre l’électron et le proton. Par conséquent, selon la 2me loi de Newton,

v2 e2
me ar = = 4πε0 2 (7.14)
rn rn

où me et e sont, respectivement, la masse et la charge de l’électron. On résout (7.13) et (7.14) pour arriver à

~2 n2 1 me e2
rn = 4πε0 et v = .
me2 e2 4πε0 ~ n

En particulier, r1 = 0.529Å = 5.29m−11 , nommé rayon de Bohr a0 , est le rayon de la première orbite (état atomique 1s).

7 Erwin Schrödinger, Quantisierung als Eigenwertproblem , Annalen der Physik 79 p.361 et 489 ; 80, p.437 ; 81, p.109 (1926).
8 L. D. Landau et E. M. Lifshitz, Mécanique quantique, théorie non-relativiste, 3me édition, Nauka, Moscou, 1974 ; voir la fin de la chap. III§21.
9 Louis de Broglie, Recherches sur la théorie des quanta (thèse de physique), Paris, 1924.
10 Comparer au TP «Vitesse du son» de L2 où le son est propagé dans un trombone à circonférence variable ; voir la sec. 8.1.4 et 8.1.5.
20 5 décembre 2019 Ondes nondispersives D. Sadovskií

8 Ondes nondispersives (séance 2)

8.1 Cours
On continue de couvrir des différents aspects des ondes non-dispersives. On donne surtout des exemples physiques,
notamment la corde (un prototype pour un nombre de systèmes importants) et l’acoustique.

Exemple Obtenir l’équation d’onde pour la corde.

Ondes scalaires Rappel sur ondes transversales et longitudinales.
Ondes acoustiques discussion préliminaire, TD expérimental.
Ondes vectorielles Exemple : ondes électromagnétiques. Rappels. Polarisation.
Propagation Onde progressive et stationnaire. Noeuds, conditions limite.
Superposition des ondes. Phénomène d’interférence. Exemple : deux sources.
Principe de Huygens Phénomène de diffraction.
Atténuation Terme additionnel dans l’équation d’onde, solution dans le cas sinosoïdal.
Dispersion Définition, terme additionnel dans l’équation d’onde.

8.1.1 Exemple physique : équation d’onde pour une corde

Prenez une corde, élastique et flexible, tendue initialement, en position de l’équilibre, le long de l’axe x (horizontal). La
tension T > 0 et uniforme. (Notez que ici T donne la valeur absolue des forces de tension.) La masse linéique de la corde
µ > 0 [kg/m] est constante elle aussi. On considére le mouvement strictement transversal dand tout points de la corde,
voir la figure ci-dessous, où ce mouvement est vertical.

masse µdx
≈T

α(x + dx)
✻

α(x)
≈T f (t, x + dx)
✻
f (t, x)
❄ ❄

T T

x x + dx

Le déplacement de chaque point x de la corde à chaque instant du temps t est donné par la fonction f (t, x). Considérez un
petit segment de la corde de longeur dx situé au point x. Deplacé de son équilibre, ce segment peut être étiré faiblement
(car la corde est supposée élastique), mais sa masse restera la même, dm = µdx. En plus, si les déplacements sont petits,
ces étirements sont négligeables et la force de tension restera inchangée ≈ T . La valeur absolue de cette force projetée
sur l’axe vertical sera T | sin α(t, x)| où α(t, x) est l’angle que fait la corde avec l’axe horizontal au point x au moment t.
Par définition de f (t, x) nous avons 
∂f 
α(t, x) = .
∂x t,x

En plus, parce que tout est petit, cet angle est lui aussi petit et

∂f 
sin α(t, x) ≈ α(t, x) = .
∂x t,x

Notons aussi que le mouvement de chaque point est accéléré et l’accélération

∂ 2 f 

a = a(t, x) =
∂t2 t,x
D. Sadovskií Ondes nondispersives 5 décembre 2019 21

varie avec position sur la corde x et temps t. Alors, on a enfin tout pour appliquer la deuxème loi de Newton

∂ 2 f 
  
∂f  ∂f 
(dm)a = µdx 2  ≈ −T sin α(t, x) + T sin α(t, x + dx) ≈ −T +T
∂t t,x ∂x t,x ∂x t,x+dx

On constate que la somme des forces sera nulle pas justement si le déplacement vertical f (t, x) ne varie avec x, mais
aussi si sa dérivée ∂f /∂x est constante ! En d’autres mots, la corde doit étre courbée (légèrement !). La courbure est
naturellement donnée par la deuxème dérivée
 
∂f  ∂f 
−
∂x t,x+dx ∂x t,x ∂ 2 f 

:= .
(x + dx) − x ∂x2  t,x

Alors nous obtenons l’équation d’Alembert pour la corde

∂2f  ∂ 2 f  ∂ 2 f  2 
   
2 ∂ f
p
µdx 2  = T dx ⇒ = v , avec v = T /µ. (8.1)
∂t t,x ∂x2 t,x ∂t2 t,x ∂x2 t,x

où v est la vitesse de propagation de phase des ondes transversales dans la corde tendue, voir (7.3).

8.1.2 Ondes acoustiques, notions préliminaires, surpression

Il s’agit d’ondes mécaniques, comme dans le chapitre précédent, mais le mouvement de
vibration est longitudinal, parallèle à la direction de propagation. Dans un gaz, la pression
est souvent assez élevée, de l’ordre de la pression atmosphérique, P ≈ 105 Pa. Quand une
vibration est produite et se propage, elle est véhiculée par une surpression locale et dépendant
du temps p(x, t), avec |p| ≪ P . Comme les phénomènes acoustiques sont très rapides,
l’équilibre des températures n’est pas atteint. L’évolution d’un volume V est reliée à celle de
la pression par la loi adiabatique

P V κ = cte ⇒ dP/P + κ dV /V = 0, avec κ = Cp /Cv (8.2)

où Cp et Cv sont des chaleurs spécifiques à pression constante et à volume constant. La

constante adiabatique κ est connue aussi comme coefficient de Lapalace. La forme diffé-
rentielle (8.2) au point de fonctionnement suffit, avec la valeur de κ en ce point ; la forme
intégrale impliquant que κ reste constant sur une grande plage n’est pas nécessaire.
Orgue
Exercice 8.1. Proposez une expérience permettant de mesurer la vitesse du son. Duquel
matériel aurez vous besoin ? Expliquez le principle des mesures.

Exercice 8.2. Réalisez cette expérience avec le matériel de démonstration : haut parleur (HP) d’impédance 8 ou 4 Ohm,
génerateur basses fréquences (GBF) avec sortie HP pré-amplifiée, oscilloscope à deux canaux, microphone avec un pré-
ampli et, possiblement, un filtre passe bande.

8.1.3 Équation d’ondes acoustiques

Considérons un cylindre de section constante S le long de l’axe d’abscisse x contenant un gaz de pression P , température
T et masse volumique µ. Considérons un petit volume V = S dx de masse m = µ S dx se trouvant aurepos dans x.
Imaginons une déformation (longitudinale !) de l’état d’équilibre telle que le point x se déplace vers x + u(x, t). Par
conséquent, le volume V se déforme et devient

  ∂u 
S x + dx + u(x + dx, t) − x − u(x, t) = S dx + S dx = V + dV.
∂x x,t

On obtient la variation relative dV /V = ∂u/∂x et utilise (8.2) pour trouver la surpression


∂u 
p = dP = −κP . (8.3)
∂x x,t

La somme des forces agissant sur V de gauche S P + p(x, t) et de droite −S P + p(x + dx, t) est
   

∂ 2 p 
  
    ∂p  ∂p 
F = S P + p(x, t) − S P + p(x + dx, t) = −S  dx = −V = κP V ,
∂x x,t ∂x x,t ∂x2 x,t
22 5 décembre 2019 Ondes nondispersives D. Sadovskií

où on applique (8.3). Par consequent, la 2me loi de Newton pour V s’écrit

∂ 2 u  ∂ 2 p 
 
m ẍ = m ü = µ V = F = κP V ,
∂t2 x,t ∂x2 x,t

d’où on obtient l’équation d’onde

∂2u ∂2p κP RT
2
= v 2, avec v 2 = =κ . (8.4)
∂t ∂x µ M
Ici on suppose l’équation d’état de gaz parfait qui pour une mole prend la forme P V = P (M/µ) = RT , avec la masse
molaire M et la température absolue T .
En prenant M = 29g (moyenne entre l’oxygène et l’azote selon leurs proportions respectives), T = 300K, et κ = 1.4
pour un gaz diatomique, on calcule la vitesse du son v = 347 m/s. Notez que v est largement inférieure à la vitesse c de
la lumière. Cela explique le décalage entre la perception de la lumière et celle du son d’un éclair.

Remarque La relation (8.3) rappelle la formulation F = kdx de la loi de Hook (loi des ressorts) ou encore f = A∂l/l
avec module d’Young A. On peut montrer que la vitesse du son dans les solides est
A
v2 = avec la masse linéique µ.
µ

8.1.4 Ondes stationnaires. Noeuds, conditions limite

La flûte, l’orgue, etc

A sin(ωt − kx) + A sin(ωt + kx) = 2A sin(ωt + π/2) cos(−kx − π/2)

Les notes Le fréquence νn de la note située n demi-tons au-dessus de la note de départ ν0 est donnée par la formule

νn = ν0 × an avec a = 21/12

On s’aperçoit que la fréquence croit de manière géométrique par rapport à la note.

Exercice 8.3. On joue de la guitare. La corde pour la note La3 de la 3ème octave donne la fréquence νLa3 = 440 Hz
(fréquence de référence). Le manche de la guitare est de longueur L = 0.5 m, la masse linéique de la corde est µ = 7 g/m.
1. Caractérisez le type des ondes (transversal, longitudinal, sinusoïdal, progressif, stationnaire etc).
2. Quelle est la longueur d’onde ? Quelle est sa forme ? Pourquoi ?
3. Proposez la fonction d’onde f (t, x) correspondante.
4. Quelle est (approximativement) la vitesse de propagation ?
5. Quelle est (approximativement) la tension de corde nécessaire pour jouer cette note ? Quelle masse m aurait on
besoin de suspendre pour produire une telle force ?
6. A présent, on fait vibrer la corde et on entends d’autres notes en plus de la note La désirée de la 3ème octave.
Pourquoi ? Quelles sont les fréquences ν et longueurs d’onde λ de ces notes ? Pourquoi ?
Commentaire : Assumez les ondes non-dispersives et négligez atténuation.
Exercice 8.4. Quel est la longueur du pipe de l’orgue donnant la note Re de la 2ème octave, fréquence 147 Hz ? Comment
le longeur varie-t-il pour les notes de la même octave ?
D. Sadovskií Ondes nondispersives 5 décembre 2019 23

8.1.5 Interférence
x
Expérience de Young On considère deux sources ponctuelles
(ou deux fentes parallèles à l’axe z) situées dans les points A et A
B, avec |AB| = ℓ et AB parallèle à Ox, et l’écran situé dans O
le plan xOz parallèlement à AB à la distance L ≫ ℓ (voir la fi- y
gure). L’axe Oy perpendiculaire à AB intersecte AB au milieux B
X
entre A et B. Prenez point X dans l’écran xOz à la distance x
du centre O.
Au départ de chaque source, l’onde est donnée par Experience de Young
I0 cos(ωt − k.r + φ0 ), avec φ0 = 0.
A l’arrivée au point X situé dans le plan de l’écran, nous avons11
   
|BX| + |AX| |BX| − |AX|
I(x, t) = Ir cos(ωt − k |AX|) + Ir cos(ωt − k |BX|) = 2Ir cos ωt − k cos k ,
2 2
où dans la limite de L ≫ ℓ, en prenant ε = ℓ/L et x ≈ ℓ, les distances |AX| et |BX| sont données par12
s 2 " #
 2
ℓ (ℓ/2 − x)
|AX| = L2 + −x ≈L 1+ , (8.5a)
2 2L2
s 2 " #
 2
ℓ (ℓ/2 + x)
|BX| = L + 2 +x ≈L 1+ , (8.5b)
2 2L2

et par conséquent
|BX| + |AX| ℓ2 + x2 |BX| − |AX| xℓ
≈ L(1 + ) = L + O(ε2 ), et ≈− .
2 2L2 2 2L
Alors 
ℓk  x
I(x, t) = 2Ir cos(ωt − kL) cos − x = 2Ir cos(ωt − kL) cos πε (8.6)
2L λ
√ √
Pour concretiser ce résultat, notez que pour les fentes (donc l’ondes cylindriques) Ir = I0 / r ≈ I0 / L+O(ε2 ), et qu’on
enregistre typiquement la moyenne en temps (durant une période T = 2π/ω) du carré de l’amplitude, i.e., la puissance

1 T
Z
P (x) = [I(x, t)2 ] = I(x, t)2 dt.
T 0
Exercice 8.5. Considérez deux sources ponctuelles émettant les ondes harmoniques (sinusoïdales) E1 et E2 de même
pulsation ω, de même kkk, de même phase initiale ϕ0 = 0 et de même amplitude E0 qui se propagent dans le vide.
Les sources sont situées dans les points (0, 0, 0) et (0, 0, 1). Trouvez le champ E(t) au point M = (l, 0, 0) de l’axe x.
Trouvez la grandeur instantanée PM (t) du vecteur de Poynting P et sa valeur moyenne P M . Précisez le sens physique de
P . Trouvez le (les) distances l > 0 pour laquelle(s) P M et maximale.

8.1.6 Diffraction
Le principe de Huygens : tout point du front d’onde agit comme une source secondaire.

Experience de Fraunhofer pour une seule fente d’une petite largeur : Au present, considerez la fente AB (coupée dans
l’ecran parallèlement à l’axe z). Selon le principe de Huygens, tout point de AB représente une source ponctuelle secon-
daire (des ondes cylindriques). En prenant à chaque fois deux points A′ et B ′ de AB symétriquement par rapport au centre
O, et en utilisant le résultat (8.6), nous trouvons l’intensité au point X comme une somme
Z ℓ 
|A′ B ′ | k
  
2L ℓk
I(x, t) = 2Ir cos(ωt − kL) cos − x d|A′ B ′ | = 2Ir cos(ωt − kL) sin x .
0 2L kx 2L
Par conséquent
2L
P (x, t) = P0 sync(2πx/a) avec la période spatiale a =
λ = 2ℓ/F,
ℓ
où F est parfois appelé nombre de Fresnel. On constate, que ℓ/L doit être assez petit, ou, également, F ≪ 1 pour voir
l’alternance des intensités sur l’axe x à l’échelle plus grande que ℓ.
11 rappelez que cos a + cos b = 2 cos a+b cos a−b
12 utilisez
√
2 2
le développement limité pour 1 + ε2 ≈ 1 + 21 ε2 + . . .
24 5 décembre 2019 Ondes nondispersives D. Sadovskií

8.1.7 Terminaison, impedance

Génération. Impédance Dans la sec. 8.1.1 et sa figure on implique que la corde est infinie, c’est à dire, elle va de −∞
à +∞ le long de l’axe x. Evidemment, cela est une idéalisation (car tout qui est bien ou mal arrive à sa fin). Au présent
donc, considerons une corde qui commence à x = 0 et qui va vers +∞. La fonction d’onde est
p
f (t − x/v) = f (φ) avec la vitesse de phase v = T /µ.

Quelle est la force périodique F (t) nécessaire pour initier (générer) les ondulations transversales au point x = 0 ?. Selon
le bilan des forces dans la sec. 8.1.1
  
 ∂f  1 ∂f  p ∂f 
F (t) = −T sin α(t, x) x=0 ≈ −T
 = −T = Tµ ,
∂x t,x=0 v ∂t t,x=0 ∂t t,x=0

où nous avons utilisé

∂f ∂f
f ′ (φ) = = −v .
∂t ∂x
Par l’analogie avec les oscillations forcées (sec. 4), on observe que ∂f /∂t est la vitesse ẏ(t) dans x = 0. Le coefficient
p
Z0 = T µ, avec T = Z0 v et µ = Z0 /v,

est appelé l’impédance caractéristique. Elle dépend seulement de la corde et sa tension, c’est à dire, l’impédance est la
propriété du système. Dans des autres systèmes nous utiliserons Z0 et v comme leur paramètres.

Terminaison Retournons la situation précédant. Prenons une corde qui arrive de x = −∞ et qui se termine à x = 0.
Quelle est la force Fterm (t) nécessaire pour compenser la partie «manquant» de la corde ? Par les mêmes arguments

∂f 
Fterm (t) = −Z0 ,
∂t t,x=0

c’est à dire, la même force F (t) mais avec le signe opposé. Que se passe-t-il si cette compensation n’est pas faite ?

Exercice 8.6. Proposez l’analogie électrique pour la terminaison (appelée dumper en anglais).

8.1.8 Réflexion et transmission

Terminaison correcte (Z0 ) et incorrecte (Z1 ). Solution générale comme une somme de deux ondes (voir ondes station-
naires, sec. 8.1.4), l’onde incident et l’onde réfléchie, avec les amplitudes différentes qui se propagent en directions
opposées. La fonction d’onde est
f (t, x) = fi (t − x/v) + fr (t + x/v)
Coefficient de réflexion R.
D. Sadovskií Ondes nondispersives 5 décembre 2019 25

8.1.9 Atténuation
Exemples : La corde dans une liquide visqueuse. Diapason avec un résonateur. Dans
chaque cas expliquez l’origine des pertes.

Exercice 8.7. Imaginez un diapason avec un résonateur «parfait» dont les parois n’absorbent
pas les vibrations et ne laissent pas les ondes acoustiques sortir. On frappe le diapason. Va-
t-il sonner infiniment ? Si non, quel serait l’origine des pertes ? Comment pourrait on mettre
en évidence—voir enregistrer ces pertes ?

Retournez à l’équation de la corde (sec. 8.1.1). Soit −β(∂f /∂t) dx une force de frottement Diapason
agissant sur le segment de corde dx. On obtient donc l’équation suivante

∂2f  ∂ 2 f  ∂ 2 f  2 
     
∂f  ∂f  2 ∂ f
p
µdx 2  = T 2
dx − β dx ⇒ 2
+ γ = v 2
, avec v = T /µ et γ = β/µ.
∂t t,x ∂x t,x  ∂t t,x
 ∂t t,x  ∂t t,x
 ∂x t,x


Est-elle une équation d’onde ? C’est à dire, peut on trouver une solution de type f (ωt ∓ kx) ? La réponse est oui si le
nombre d’onde k est considéré comme un nombre complexe, k ∈ C.

Exercice 8.8. Comme une solution de l’equation obtenue ci-dessus, essayez la fonction

f (t, x) = A exp(i(ωt − kx)),

et trouvez la relation entre k et (ω, Γ, v). Considerez la limite de faible atténuation. Précisez la définition quantitative de
cette limite (i.e., trouvez le petit paramètre ǫ ≪ 1). Faites le comparison avec le cas des vibrations atténuées.
Solution. Suite au calcul des dérivées on obtient
r
ω γ
−ω 2 f + i ωγf = −k 2 v 2 f ⇒ k=± 1−i ∈C
v ω

Dans la limite de faible atténuation 2ǫ = γ/ω ≪ 1. On écrit

r  γ 2
ω γ  γ
k=± 1+ exp(−iφ/2) où φ = arctan ≈
v ω ω ω
ω  γ 
≈± 1−i
v 2ω
Alors la condition Re k/ Im k = ǫ ≪ 1 signifie que dans une longueur | Im k|−1 = 2v/γ s’insèrent
plusieurs logueurs
d’onde λ = 2π/| Re k|. Faites le dessin de la fonction d’onde f (ωt − kx) = exp i ωt − i (Re k)x exp −(Im k)x .


26 5 décembre 2019 Ondes nondispersives D. Sadovskií

9 Ondes dispersives
Rappels. Dispersion en optique Indice de réfraction n, que signifie-t-il ?. Lumière mono-
chrome : la relation entre couleur, fréquence ν et la longueur d’onde λ. Dispersion (décom-
position) de la lumière blanche à travers d’un prisme.

Dispersion, corde Retournez à l’équation de la corde. Si elle est courbée, à son chaque
point elle possède un moment13 de pliage (ou torque). On a l’équivalent de la loi d’Hook
pour la deformation linéaire
Prisme
∂2f
τ (x) = −kt avec kt > 0.
∂x2
Donc le torque est proportionnel à la courbure avec coefficient négatif signifiant que une corde pliée tente de se replier.
Ce torque correspondent à une force de réaction perpendiculaire à la corde et donc dirigée presque verticalement13
∂τ ∂3f
Ft ≈ = −kt 3 .
∂x ∂x
En considérant la somme des forces à deux extrémités du ségment dx on trouve
∂Ft ∂4f
dFt = dx − kt 4
∂x ∂x
Comme dans le cas d’atténuation (sec. 8.1.9), cette force doit être ajoutée à l’équation de la corde (8.1). Nous avons ainsi
l’équation d’onde nonlinéaire
"  #
∂ 2 f  2 
∂ 4 f 
  
∂f  2 ∂ f
+γ =v + 2d avec d = − 21 kt v 2 .
∂t2 t,x ∂t t,x ∂x2 t,x ∂x4 t,x

Exercice 9.1 (relation de dispersion). Nous ne pouvons plus assurer que n’importe quelle fonction f (t ∓ x/v) peut
satisfaire cette équation, car une relation entre ω et k nous sera imposée. Essayez une fonction particulière
f (t − x/v) = A exp(i ωt − i kx)
Trovez la relation entre ω et k, appelée la relation de dispersion. Considerez le cas de faible dispersion.
Solution. Suite au calcul des dérivées
∂f ∂2f ∂2f ∂3f ∂4f
= i ωf, = −ω 2 f, = −k 2 f, = −k 2 (−i k)f, = k4 f
∂t ∂t2 ∂x2 ∂x3 ∂x4
on obtient p ω p
ω = ±vk 1 − 2dk 2 et vφ = = v 1 − 2dk 2 ,
|k|
où vφ = vφ (k) est appelée vélocité de phase. Pour d > 0 (et bien sur 2dk 2 < 1) elle est inférieure à la vitesse de
propagation v. Soit ǫ = 2dk 2 ≪ 1 (condition de faible dispersion). En développant, on obtient la forme standard de la
relation de dispersion
ω ≈ vk − dk 3

Dispersion normale et anormale Les cas de d < 0 et d > 0 correspondent à la dispersion normale et anormale.

Dispersion, ondes stationnaires Exemple : une corde de longueur L fixée à deux extrémités (guitare). Soit ν1 la fré-
quence fondamentale, accordée par exemple à 440 Hz (La 3ème). Quelles sont les fréquences des harmoniques si elles
sont dispersives ? Les conditions aux limites impliquent que les vecteurs d’onde sont
π
kn = n , avec k = 1, 2, 3, . . .
L
Pour les fréquences on obtient
v dπ 2 3
2πνn = ωn = vkn − dkn3 ⇒ νn = n− n
2L 2L3
Exercice 9.2. Combinez les effets de la dispersion et faible atténuation. Donnez la solution.

13 Rapplez que par rapport au point O le moment M ~ de la force F

~ appliquée dans le point O′ est un vecteur OO ~ ′∧F ~ . Ici on considère le torque ∂τ
par rapport au point O = x avec |OO~ ′ | = ∂x ; la direction de OO
~ ′ est presque le long de l’axe x et F
~ se rapproche à la verticale.
D. Sadovskií Questionnaire éclair QE1 15-12-2019 Q1

Questionnaire éclair de la rentrée

Répondez, svp, sur cette feuille et au verso s’il vous est nécessaire.
Prénom/Nom
D’où, de quelle formation venez vous ?
Votre moyenne en physique
Votre moyenne en math
Que voulez vous faire après L3 ?
Que attendez vous du cours «Vibrations» ? En quoi ce cours vous peut être utile ?

Répondez, svp, aux questions seulement si vous savez les réponses ! Vous avez 20 minutes.
Problème 1 Soit a et b des nombres réels, a, b ∈ R et

z = a + ib, z ∈ C,

un nombre complexe. On considère la representation de Euler

z = |z| exp(iϕ), |z| ∈ R>0 et ϕ ∈ [0, 2π),

avec |z| un nombre réel non-négatif, et angle ϕ ∈ [0, 2π).

Comment appelle-t-on |z| et ϕ ? Exprimez |z|, tan ϕ, cos ϕ, et sin ϕ en fonction de a et b.
Problème 2 On considère au présent une fonction réelle

x : R → R : t 7→ x(t)

du temps t, telle que

dx/dt := ẋ = ax, a ∈ R.
Quelle est la forme explicite de x(t) ? Quelle est son allure pour a > 0 ? Pour a < 0 ? Quelle est la valeur x(0).
Problème 3 On considère deux autres fonctions réelles

y1 : R → R : t 7→ y1 (t) et y2 : R → R : t 7→ y2 (t)

du temps t, telles que

ÿ1 = −ω 2 y1 et ÿ2 = ω 2 y2 , avec ω ∈ R>0 .
Quelle peut être la forme explicite de y1 (t) et y2 (t) ? Combien y-a-t-il des solutions indépendants dans chaque cas ?
Problème 3∗ Au present, trouvez la solution y1 (t) et y2 (t) pour le système

ÿ1 = −ω 2 y1 et ÿ2 = ω 2 y1 , avec ω ∈ R>0 .

Problème 4 Enfin on considère une fonction complexe

z : R → C : t 7→ z(t)

du temps réel t, telle que

d2 z/dt2 := z̈ = −ω 2 z, ω ∈ C \ {0}.
Donnez la forme explicite de z(t). Considérez des cas différents pour des valeurs ω ∈ R>0 , iω ∈ R \ {0}, et une valeur
générique avec Im ω 6= 0 et Re ω 6= 0. Comparez au problème 3.
Problème 5 On se rapproche enfin de notre sujet. Considerez un ressort sans poids. Sa longeur est l, et sa constante de
raideur est k. Le ressort est suspendu verticalement au plafond de l’hauteur h, et à son bout bas libre on attache un corps
ponctuel de masse m. On laisse le corps vibrer en le long de l’axe vértical z autour du point z0 . Quel est la valeur de z0 si
le plancher est situé au z = 0. Faites un dessin.
Q2 22-12-2019 Questionnaire éclair QE2 D. Sadovskií

Questionnaire éclair : oscillations harmoniques libres

Prénom/Nom

Répondez, svp, sur cette feuille et au verso s’il vous est nécessaire.
Répondez, svp, aux questions seulement si vous savez les réponses ! Vous avez 15 minutes.
Problème 1 Considerez deux systèmes (a) et (b) à 1 degré de liberté composés chaqun d’un ressort de rigidité k avec un
bout fixe, et avec une masse ponctuelle m attachée au bout libre. Dans le cas (a) le ressort et posé horizontalement sur un
plan, dans le cas (b) il est suspendu verticalement. Il n’y a pas de friction, ni d’autres pertes, mais on tient compte, bien
sur, de la force pesanteur. Quelle est la pulsation propre ω0 des oscillations libres dans chaque cas ?

Problème 2 On considère une fonction réelle

x : R → R : t 7→ x(t)

du temps réel t, telle que

d2 x/dt2 := ẍ = −ω 2 z, ω ∈ R>0 .
Au moment t = 0, les valeurs initiales sont x(0) = x0 et ẋ(0) = v0 > 0. Trouvez x(t).
Problème 3 Lequel des deux systèmes dont le portrait de phases est presenté sur la figure ci-dessous est stable ? unstable ?
D. Sadovskií Questionnaire éclair QEQ3 29-12-2019 Q3

Questionnaire éclair : representation matricielle

Prénom/Nom

Répondez, svp, sur cette feuille et au verso s’il vous est nécessaire.
Répondez, svp, aux questions seulement si vous savez les réponses ! Vous avez 15 minutes.
Problème 1 Considerez l’équation différentielle ordinaire de deuxième ordre

ẍ + ω02 x = 0

qui représent les oscillations harmoniques. Comment peut on remplacer cet équation par un système des équations de
premier ordre. Combien équations de premier ordre a-t-on besoin ? Quelle est la représentation matricielle de ce système.
Problème 2 Donnez la forme matricielle correspondant à l’équation différentielle ordinaire de deuxième ordre

ẍ + γ ẋ + ω02 x = 0, avec γ, ω0 ∈ R \ {0}

Problème 3 Trouvez les valeurs propres de la matrice suivante

 
a b
A= , avec a, b, c ∈ R \ {0}
b c

Quelle est la relation entre det A et les valeurs propres de A ?

Q4 12-12-2019 Questionnaire éclair QE4 D. Sadovskií

Questionnaire éclair : oscillations harmoniques libres amorties

Problème 1 Le mouvement d’un système à un degrè de liberté est décrit par l’équation différentielle de 2ème ordre
ω0 ω0
ẍ + ẋ + ω02 x = 0, où Q ∈ R>0 et γ :=
Q 2Q

Pour Q = 1, trouvez la solution particulière x(t) : t 7→ R telle que

x(0) = 1 et ẋ(0) = −1

En sachant que l’énergie cinétique T et potentielle V du système sont, respectivement,

1 1 2 p
mẋ2 et kx , où ω0 = k/m,
2 2
trouvez l’evolution de l’énergie totale T + V du système.
Problème 2 Considerez le cas Q = 1/4 et trouvez la solution et l’énergie du système pour les mêmes conditions initiales.
D. Sadovskií Questionnaire Q5 12-2019 Q5

Questionnaire : fonctions d’une et de plusieures variables

Problème 1. Considérez la fonction connue en physique comme potentiel de Lennard–Jones

b a
fLJ : R>0 → R : x 7→ 12
− 6, avec b, a > 0.
x x
Quelle est l’image de cette fonction dans R ? Tracez le graph de fLJ (x). Possède-t-elle un (ou plusieures) point(s) station-
naire(s) ? Si oui, trouvez le(s). Exprimez la coordonnée x0 de ce point et la valeur stationnaire fLJ (x0 ) en terme de a et b.
En supposant que x est la distance entre deux molécules, et que fLJ (x) donne l’énergie d’interaction moléculaire, trouvez
l’énergie de dissociation D de l’agrégat. En développant fLJ (x) au voisinage de x0 ,

fLJ (x) ≈ fLJ (x0 ) + 12 k(x − x0 )2 + . . .

trouvez la constante de rigidité k (constante de force) du lien entre les molécules comme fonction de a et b. Enfin,
exprimez a et b en fonction de D et k.
Problème 2. Considérez la fonction connue en physique comme potentiel de Hénon–Heiles

fHH : R2 → R : (x, y) 7→ 12 (x, y)2 + ε( 13 x3 − xy 2 ), avec 0 < ε ≪ 1.

Quelle est l’image de cette fonction dans R ? Calculez ∇fHH . La fonction fHH , possède-t-elle un (ou plusieurs) point(s)
stationnaire(s) ? Si oui, trouvez et caractérisez le(s). Dans le plan R2x,y , tracez les lignes des niveaux constantes (re-
présentation «cartographique») de fHH (x, y). Dans le(s) point(s) stationnaires (x0 , y0 ) trouvez les valeurs stationnaires
fHH (x0 , y0 ). Pour le(s) point(s) stable(s), trouvez le(s) constante(s) de rigidité kx et ky . Considérez les applications
g : R2 → R2 , dites opérations de symétrie. Existe-t-il une (ou plusieures) opérations g pour laquelle (lesquelles) fHH
reste invariante, c’est à dire f ◦ g ≡ f .
Problème 3. Dans l’espace R3 , considérez le système de trois charges électriques ponctuels placés dans r1 , r2 , r3 (pensez
d’une molécule H3 ). Le potentiel Coulombien créé par chaque charge est 1/kr − ri k. Soit r1 , r2 , r3 forment un triangle
équilatéral dans le plan (x, y, 0). Pour la fonction potentielle du système f , trouvez ∇f . Montrez que f possède un point
critique x0 au milieu du triangle. Caractérisez ce point. Pour un trinagle assez grand, développez f dans x0 jusqu’à les
termes réfléchissant la symétrie du système.