Vibrations Et Ondes

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 94

REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE

Ministère de l’Enseignement Supérieur et de la Recherche Scientifique


Université TAHRI Mohammed Béchar
Faculté des Sciences Exactes
Département des Sciences de la Matière

N° d’ordre : UTMB/FSE/PP/

Filière : Physique

Spécialité : 2eme Année Licence SM-ST

Module : Vibrations et ondes

Vibrations et Ondes
Cours et Exercices
Bennaceur Said

Expertisé par :
Dr . Kadri Siham
Dr . Atouani Toufik

Année Universitaire :2020/2021


Préambule

Le présent polycopié, intitulé : « Vibrations et Ondes Cours et


Exercice » est élaboré et présenté en conformité au canevas
relatif à la formation Licence LMD-S3 dans le domaine
Science de la matière (SM) et Science et Technologie (ST).Ce
cours est structuré en deux parties :
La première, répartie en Trois chapitres, traite le problème
des vibrations. Le premier chapitre porte les Généralités sur
les vibrations. Qui décrit la Définition d’un mouvement
vibratoire, exemples de systèmes vibratoires et mouvements
périodiques. Dans le deuxième chapitre est destiné en
premier partie pour les systèmes linéaires à un degré de
liberté ou il traite les oscillations libres, l’oscillateur
harmonique, pulsation propre d’un oscillateur harmonique et
l’énergie d’un oscillateur harmonique. La deuxième partie de
ce chapitre est concernée pour traiter les oscillations libres
amorties. Concernant la troisième partie de ce chapitre traite
les oscillations libres forcées qui prend en compte le cas
d’une excitation sinusoïdale (résonance, déphasage). Et enfin
la quatrième partie de ce chapitre est consacrée pour les
oscillations amorties forcées où elle décrit les deux modes de
vibration (transitoire et permanente), avec une comparaison
entre systèmes oscillants mécaniques et électriques. Le
troisième chapitre et le dernier de cette partie cite les
vibrations aux plusieurs degrés de liberté, car il commence
par une étude sur les systèmes à deux degrés de liberté et il
généralise avec les systèmes a N degrés de liberté. La
deuxième partie qui est intitulée par « les ondes
mécaniques », consacrée au traitement des phénomènes de
propagation des ondes. Car elle est divisée en quatre
chapitres, le premier donne des généralités sur les ondes
mécaniques initialisant par leurs classifications puis son
équation générale finalisant par les caractérisations de ces
ondes mécaniques. Le deuxième chapitre traite les ondes
transversales sur une corde qui produisent, avec les ondes
incidentes et réfléchies les ondes stationnaires. Le troisième
chapitre est consacré aux ondes longitudinales dans les
fluides finalisant par l’effet Doppler.
Le quatrième chapitre est consacré aux Ondes élastiques
dans les solides.
Les cours présentés avec un enchaînement logique, chaque
nouveau concept défini est clarifié par des exemples simples
et utiles, une série d’exercices dans chaque chapitre venant
enrichir le cours, le tout a été réalisé avec l’esprit de
permettre une meilleure assimilation par l’étudiant.

Dr. Bennaceur Said


SOMMAIRE

Préambule
PARTIE I : VIBRATION
CHAPITRE 01 : Généralités sur les vibrations
1. Définition d’un mouvement périodique 02
2. Définition d’une oscillation 02
3. Définition d’un mouvement sinusoïdal 03
4. Nombre de liberté 04
5. Représentation complexe d’un mouvement vibratoire 04
6. Définition des séries de Fourier 04
Exercices 05
CHAPITRE 02 : Systèmes linéaires à un degré de liberté
2.1 Les oscillations libres 06
2.1.1 Définition d’un oscillateur libre 06
2.1.2 Définition d’un oscillateur harmonique 06
2.1.3 Equation du mouvement 06
2.1.4 Energie d’un oscillateur harmonique (E) 08
2.1.5 Equation de Lagrange 10
2.1.6 Conditions d’équilibre et d’oscillation 11
Exercices 13
2.2 Oscillateur amortie à un degré de liberté 14
2.2.1 Définition d’un oscillateur amortie 14
2.2.2 Equation de Lagrange 14
2.2.3 Equation du mouvement 14
2.2.3.1 Résolution de l’équation du mouvement 15
2.2.3.2 Décrément logarithmique 17
Exercices 18
2.3 Oscillations amorties forcées 20
2.3.1 Définition d’un oscillateur amortie forcé 20
2.3.2 Equation de Lagrange 20
2.3.3 Equation du mouvement 20
2.3.3.1 Résolution de l’équation du mouvement 21
2.3.3.2 Résonnance 23
2.3.3.3 Résonance de phase 24
2.3.3.4 Bande passante 24
2.3.3.5 Analogie électromécanique 26
Exercices 27
CHAPITRE 03 : Systèmes linéaires à plusieurs degrés de liberté
3.1 Systèmes linéaire de plusieurs degrés de liberté 29
3.1.1 Les type de couplages 29
3.2 Systèmes linéaire deux degrés de liberté cas libre 32
3.2.1 Equation du mouvement 32
3.2.2 Les modes d’oscillations 33
3.3 Système force à deux degrés de libertés 34
3.4 Oscillations forcées d’un système non amorti à deux degrés de liberté 35
3.5 Généralisation aux systèmes à n degrés de liberté 35
Exercices 39
PARTIE II : LES ONDES MECANIQUES
CHAPITRE 04 : Généralités sur les ondes mécaniques
4.1 Classification des ondes 41
4.2 Intégrale général de l’équation générale d’ondes progressives 42
4.2.1 Solution de l'équation de propagation 42
4.3 Vitesse de phase, vitesse de groupe 43
4.4 Front d’onde 44
4.4.1 Exemple des ondes planes, ondes sphériques 44
4.5 Réflexion et transmission des ondes 46
Exercices 48
CHAPITRE 05 : Ondes transversales sur une corde
5.1 Introduction 50
5.2 Equation de propagation 50
5.3 Impédance caractéristique 51
5.4 Ondes stationnaires 52
5.4.1 Définition 52
5.4.2 Solution stationnaire de l’équation de l’Alembert 52
5.4.3 Réflexion et transmission 53
5.4.4 Types de condition aux limites pour une corde vibrante 53
5.4.4.1 Deux extrémités fixes 53
5.4.4.2 Une extrémité fixe et l’autre libre 55
5.4.4.3 Deux extrémités libres de la corde 57
5.4.4.4 Autres types de conditions aux limites 58
5.5 Partie expérimentale 59
5.5.1 Expérience générale 59
5.5.2 Energie emmagasinée sur une corde vibrante 59
5.5.3 Énergie cinétique 59
5.5.4 Energie potentielle 60
Exercices 61
CHAPITRE 06 : Ondes longitudinales dans les fluides
6.1 Introduction 65
6.2 Ondes planes dans un tuyau cylindrique 65
6.2.1 Equation d’ondes dans un gaz 67
6.2.2 Ondes planes dans un tuyau cylindrique 67
6.2.2.1 Equation d’ondes dans un fluide 69
6.2.2.2 Vitesse du son 69
6.2.2.3 Impédance acoustique et impédance caractéristique 70
6.2.2.4 Energie transportée par une onde 71
6.2.2.5 Intensité 74
6.2.2.6 Niveau sonore 75
6.2.2.7 Coefficients de réflexion et de transmission d’ondes 75
(conditions aux limites)
6.2 Effet Doppler 77
6.3.1 Historique 77
6.3.2 Définition 78
Exercices 80
CHAPITRE 07 : Ondes élastiques dans les solides
7.1 Introduction 85
7.2 Définition et propriétés élastiques d’un solide 85
7.3 Equation d’onde dans barreau solide 86
7.4 Résolution de l’équation de d’Alembert 87
7.4.1 Un barreau infini 88
7.4.2 Impédance mécanique 88
7.4.3 Aspect énergétique 88
Exercices 90
Partie I : VIBRATIONS

Dr. Bennaceur Said Page 1


CHAPITRE 01 : Généralités sur les vibrations

CHAPITRE 01 : Généralité sur les vibrations


1. Définition d’un mouvement périodique
Un mouvement est dit périodique s’il répète identique à lui-même pendant des
intervalles de temps égaux.
Exemples :

 Le mouvement de révolution de la Lune: La lune effectue un cycle complet de


révolution autour de la terre en environ 29 jours.

Figure 1-1 : Mouvement de révolution de la lune

 Les battements du cœur: un battement de cœur est une succession de


contractions et de relâchement des muscles cardiaques qui actionnent des valves
et provoquent la circulation du sang dans le corps.

Figure 1-2 : Électrocardiogramme

2. Définition d’une oscillation


On appelle oscillation, un mouvement qu’il s’effectue autour d’une position
d’équilibre.

Dr. Bennaceur Said Page 2


CHAPITRE 01 : Généralités sur les vibrations

Exemples :
 Masse –ressort :
Soit un ressort de longueur à vide x0. On
modifie sa longueur en exerçant une force de
tension à son extrémité libre, en le comprimant
ou l'étirant.
Figure 1-3 : Système masse-ressort

 Circuit électrique :
Soit un circuit contenant une bobine d'inductance( L )et
d'un condensateur de capacité( C).La charge et la
décharge du condensateur s'accompagnent d'oscillations
électriques.
Figure 1-4 : Système bobine -condensateur

3. Définition d’un mouvement sinusoïdal


On appelle un mouvement sinusoïdal si sa forme :
𝑔(𝑡) = 𝐴𝑠𝑖𝑛(𝜔𝑡 + 𝜑)

[ 𝜈 Où     (1-1)
2
𝑔(𝑡) = 𝐴𝑐𝑜𝑠(𝜔𝑡 + 𝜙)
Avec g(t)
A : amplitude A
φ : phase initiale (rad)
∅ : phase initiale (rad)
ω : pulsation
T T (s)
t : temps

Figure 1-5 : Forme du mouvement sinusoïdal

 Période (T) en seconde (s) : c’est l’intervalle du temps.


 Fréquence (f) en (Hertz): le nombre des répétitions par seconde.
1
𝑓=𝑇 (1-2)
Avec :
T : est la période
 Pulsation (ω) en (rad/s) :
2
  2f (1-3)
T

Dr. Bennaceur Said Page 3


CHAPITRE 01 : Généralités sur les vibrations

Avec :
T : est la période
f : est la fréquence
 Pour la périodicité en math s’exprime : g (t  T )  g (t )
4. Nombre de liberté
On définit le nombre de liberté (dll) par la relation suivante :
dll  n  m (1-4)

Remarque : le nombre de liberté est le nombre d’équation à étudier.


n : Nombre des coordonnées généralisées indépendantes ou liées
m : Nombre des liaisons entre les coordonnées
5. Représentation complexe d’un mouvement vibratoire
 Représentation complexe : pour faciliter les calculs nous transformons les
grandeurs sinusoïdales en forme exponentielles grâce à la forme d’Euler :
e j  cos   j sin  (1-5)

6. Définition des séries de Fourier


La grandeur périodique peut être exprimée par les sommes des fonctions de sinus et
cosinus afin de la manipuler physiquement et mathématiquement.
Cette somme est appelée série de Fourier.
La série de fourier d’une fonction f(t) périodique de période (T) , est définie par :

f (t )  a 0   a n cos nt    bn sin nt  (1-6)
n 1 n 1

 (a0 , an et bn) les coefficients de Fourier


T T
T
2
1 , an   f (t ) cos(nt )dt , bn  2  f (t ) sin(nt )dt
a0 
T 
0
f (t )dt
T 0 T 0

 La pulsation (ω) est appelée pulsation fondamentale

2
 (1-7)
T

 Les pulsations (nω) multiples de ω sont appelées les harmoniques

Dr. Bennaceur Said Page 4


CHAPITRE 01 : Généralités sur les vibrations

Exercices
Exercice 01 :
Soit la grandeur sinusoïdale g(t) représentée ci-
contre.
1. Calculer l’amplitude, pulsation et la phase
initiale.
2. Déduire l’équation g(t).

Exercice 02 :
Un mouvement vibratoire est caractérisé par le déplacement suivant :

x(t )  4 cos(25t  )
2
Ou x en centimètres, t en secondes et la phase en radians.
1. Déterminer l’amplitude maximale
2. Donner la pulsation propre, la fréquence et la période du mouvement.
3. Exprimer la phase initiale (déphasage à l’origine).
4. Calculer le déplacement, la vitesse et l’accélération aux instants t=0s et t=0.5s.
Exercice 03 :
Un mobile M décrit un mouvement rectiligne sinusoïdal si son abscisse x est une fonction
sinusoïdale du temps: x(t )  2 cos(4t  6)
Trouver à l’aide de la représentation complexe
.
1. la vitesse du mouvement rectiligne sinusoïdal x(t )
..
2. l’accélération du mouvement rectiligne sinusoïdal x(t )

Exercice 04 :
Un mouvement harmonique est décrit par x(t)=X cos (100t+). Les conditions initiales
sont x(0)=4m et x (0) =1m/s.
a) Calculer X et .
b) Exprimer x(t) sous la forme x = A cos (ωt) + B sin (ωt) et en déduire les valeurs de A.

Dr. Bennaceur Said Page 5


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

CHAPITRE 02 : Systèmes linéaires à un degré de liberté


2.1Les oscillations libres
2.1.1 Définition d’un oscillateur libre
Un oscillateur libre est un système oscillant en absence de toute force d’excitation.
2.1.2 Définition d’un oscillateur harmonique
Un oscillateur harmonique est un oscillateur qui est ramené à sa position d'équilibre
lorsqu'il est déplacé d'une certaine distance en raison d'une force de rappel opposée au
mouvement.
F  c  x (2-1)
Avec
c : une constante positive.
x : la distance parcourue par l'oscillateur

Exemple : Système Masse –ressort


La masse est écartée légèrement de sa position d'équilibre et relâchée sans
vitesse initiale, l'expérimentateur constate que celle-ci se met à osciller autour de
cette position d'équilibre. La tension est la force de rappel qui ramène à sa
position d’équilibre.

Figure 2-1 : Système masse –ressort horizontal

T  k  x (2-2)

2.1.3 Equation du mouvement


L’équation du mouvement d’un oscillateur harmonique est de la forme
𝑞̈ + 𝜔02 𝑞 = 0 (2-3)

q sont les coordonnées généralisées, pour la mécanique 𝑞(𝑥, 𝑦, 𝑧, 𝜃, 𝜑 … … ) et


l’électricité 𝑞(𝑖, 𝑢, 𝑄, … … )
La solution de l’équation (2-2) s’écrit comme :
𝑞(𝑡) = 𝐴𝑠𝑖𝑛(𝜔02 𝑡 + 𝜑)
{ (2-4)
𝑞(𝑡) = 𝐴𝑐𝑜𝑠(𝜔02 𝑡 + ∅)
𝜋
∅=𝜑+2

Dr. Bennaceur Said Page 6


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

Avec :
A : amplitude
φ : phase initiale
∅ : phase initiale
ω : pulsation
t : temps

Exemple : Étude du mouvement d’un oscillateur harmonique, ressort (k) lié avec une
masse (m). Soit une masse accrochée à l’extrémité d’un ressort verticale sans masse.
Cette masse se déplace sans frottement sur le plan vertical. A t=0, on écarte ce point de
sa position d’équilibre d’une grandeur x puis on le lâche sans vitesse initiale.

Figure 2-2 : Système masse –ressort vertical

En appliquant la méthode de Newton

a. Equilibre
    
F  0 P  T1  0 P  T1  0

mg  kx0  0 (2-5)
b. Mouvement
    
F  ma P  T2  m a P  T2  ma
..
mg  k ( x0  x)  ma mg  kx0  kx  m x
 
0
.. ..
 kx  m x m x kx  0

L’équation de mouvement
.. k
x x0 (2-6)
m

En comparant cette équation par l’équation (2-3) on déduit que :

Dr. Bennaceur Said Page 7


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

𝑘 𝑘
𝜔02 = 𝑚 𝜔0 = √𝑚 (2-7)

𝜔0 est la pulsation libre.


2.1.4 Energie d’un oscillateur harmonique (E)
L’énergie d’un oscillateur harmonique est la somme de l’énergie cinétique (T) et
potentielle (U).

E  T U (2-8)
 Energie cinétique (T)
.
L'énergie cinétique d'un objet de masse m, ayant vitesse ( q ) (donc liée au
mouvement) vaut :
. . . .
1 .2
 Translation T  mq q  (v, x, y, z ) (2-9)
2
1 .2
 Rotation T  I  (2-10)
2
.
IΔ moment d’inertie  vitesse angulaire
 Energie potentielle (U)
Il y a plusieurs formes qui présentent l’énergie potentielle
 Energie potentielle de la force pesanteur :
Une masse ( m) se trouve dans un champs gravitationnel constant (g), condense
une énergie potentielle pesanteur est sous forme :

U P  mgh (2-11)

 Energie potentielle élastique


- pour Un ressort d’une constante de raideur (k) et d’une déformation (x). L’énergie
Potentielle est donnée par :
1 2
Uk  kx (2-12)
2
x

Dr. Bennaceur Said Page 8


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

- Un ressort de torsion d’un constant de raideur (k) , d’un déformation (θ)

1
U k  k 2 (2-13)
2

 Conservation de l’énergie totale


L'énergie ne peut être ni créée ni détruite, mais transférée seulement d'un système à
un autre et transformée d'une forme à une autre. Donc si l’énergie totale ( cinétique et
potentielle ) d’un système sera invariable durant son mouvement nous conclurons
que :
dE
0 (2-14)
dt
Cette équation de conservation va donner l’équation du mouvement.

Exemple
Soit une masse m accrochée à l’extrémité d’un ressort horizontal sans masse de raideur
k. Cette masse se déplace sans frottement sur le plan horizontal. A t=0, on écarte ce point
de sa position d’équilibre d’une grandeur x puis on le lâche sans vitesse initiale.

Figure 2-3 : oscillateur mécanique : système masse –ressort horizontal

 Energie totale
E=T+U (2-15)
 Energie cinétique
1 .2
T mx (2-16)
2
 Energie potentielle
1
U  U p  U k  U k  kx2 (2-17)
 2
0

1 .2 1 2
E  m x  kx (2-18)
2 2

Dr. Bennaceur Said Page 9


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

 .

1 . 2 1 2
d  m x  kx 
dE 2 2 
Le système est conservatif 0 dE
  
dt dt dt

dE . .. . .. .
 m x x  kx x  0 m x  kx  0 (2-19)
dt
.. k .
Equation de mouvement x x 0 (2-20)
m
2.1.5 Equation de Lagrange
L’équation de Lagrange permet de déterminer l’équation du mouvement des systèmes
mécaniques. Elle est décrite par l’équation suivante :

 
d  L  L
  q  Fi
(2-21)
dt   q.
 i  i

Ou ( L) est le Lagrangien qui est une fonction explicite des coordonnées généralisées et
des vitesses généralisées
L=T-U (2-22)

T : est l’énergie cinétique totale du système


U: est l’énergie potentielle totale du système
.
qi et qi sont les coordonnées et les vitesses généralisées

Fi sont les forces généralisée associées à qi .


Pour le cas du système conservatif à un degré de liberté, l’équation (2-21) se réduit à :
 
d  L  L
  q  0
(2-23)
dt   q. 
 
Exemple
Soit le système précèdent (figure 2-3), on utilisant la méthode de Lagrange écrire
l’équation de mouvement et déduire la pulsation propre.
 Lagrangien
L=T-U (2-24)

 Energie cinétique
1 .2
T  mx (2-25)
2

Dr. Bennaceur Said Page 10


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

1
U  U p  U k  U k  kx2 (2-26)
 2
0

1 .2 1 2
L
m x  kx (2-27)
2 2
 Equation de Lagrange

 
d  L  L
  x  0
(2-28)
dt   x. 
 

Après les dérivations nous recevons cette équation


..
m x kx  0 (2-29)
Si nous devisons cette équation sur (m) nous obtenons l’équation de mouvement
.. k
x x0 (2-30)
m
Par la transposition avec l’équation générale 𝑥̈ + 𝜔02 𝑥 = 0

Donc nous tirons la pulsation

𝑘
𝜔0 = √𝑚 (2-28)

2.1.6 Conditions d’équilibre et d’oscillation


 Condition d’équilibre
Suivant les notions précédentes, la force est définie comme dérivante d’une énergie potentielle
U (2-29)
F 
q

Notre système est en équilibre donc F=0


Donc la condition d’équilibre sera comme suit
𝜕𝑈
=0 (2-30)
𝜕𝑞

 Condition d’oscillation
Un équilibre est stable si la constante de l’oscillateur harmonique est supérieure à
zéro(0).
F
C 0 (2-31)
q

U
Et comme on a décrit apparemment F  (2-32)
q

Dr. Bennaceur Said Page 11


CHAPITRE 02 : Systèmes linéaires à un degré de liberté


On appelle aussi condition d’oscillation  U 
2
0 (2-33)
q 2  q 0

 Si l’équilibre est instable don on n’aura pas des oscillations

 2U 
0 (2-34)
q 2  q 0

Dr. Bennaceur Said Page 12


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

Exercices
Exercice 01
Le système ci-contre peut tourner librement autour du point O. La boule est supposée
ponctuelle et la tige sans masse. sin     et cos   1
1. Trouver l’énergie cinétique T et l’énergie potentielle U.
2. Trouver l’équation du mouvement.
3. Trouver la pulsation propre  0  sachant que m=1kg, L=2m,
k=2N/m, g=10m/s2
Exercice 02
On considère le circuit électrique constitué d’un condensateur de capacité C branché aux
bornes d’une bobine d’inductance L. Initialement le condensateur est chargé, puis on
ferme l’interrupteur S. En écrivant la loi des mailles, donner l’équation différentielle du
système en fonction de la charge q(t) du condensateur ; déduire la pulsation des
oscillations du circuit.

Exercice 03
Une tige de longueur totale L+l et de masse négligeable, porte à son extrémité supérieure
une masse ponctuelle m. L’autre bout de la tige est relié à un ressort de raideur k. Celui-ci
n’était pas déformé à l’équilibre et supposé rester horizontal lors des petits mouvements.
La tige peut tourner librement autour du point O. À l’équilibre la tige était verticale.
1. Trouver l’énergie potentielle U et l’énergie cinétique T du système
2. Trouver l’équation du mouvement et la pulsation propre 0
3. Trouver la condition d’oscillation du système.

Dr. Bennaceur Said Page 13


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

2.2 Oscillateur amortie à un degré de liberté


2.2.1 Définition d’un oscillateur amortie
Un système soumis à un frottement est dit amortie, les frottements visqueux sont de
la forme : f  v (2-35)
α : coefficient de frottement
Est schématisé par l’amortisseur

2.2.2 Equation de Lagrange


En introduisant l’amortisseur, l’équation de Lagrange devienne :
 
d  L  L .
   (2-36)
dt   q.  q q

 

Si on introduisant la fonction de dissipation


. . .
 1 .2 Translation ( x, y, z )
 D  q
1 . 2
D   q2  (2-37)
. 2
2  
 D  v 2    l  
1 1 Rotation
 2 2  
D .
Avec une petite comparaison nous remarquons que .
  q2 (2-38)
q
Donc l’équation de Lagrange de système amortie s’écrit :
 
d  L  L D
  q   .
(2-39)
dt   q.  q
 

2.2.3 Equation de mouvement


L’équation de mouvement est sous la forme :
𝑞̈ + 2𝜆𝑞̇ +𝜔02 𝑞 = 0 (2-40)
(λ) coefficient d’amortissent
Le facteur de qualité est exprimé comme suit :
0
 (2-41)
2
Exemple
Soit une masse (m) est fixée à un ressort horizontal de raideur .(k) et un amortisseur de
coefficient de frottement (α). Trouver l’équation du mouvement

Dr. Bennaceur Said Page 14


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

Figure 2-4 : Système amorti [masse, ressort] horizontal


 Lagrangien
L=T-U (2-42)
 Energie cinétique
1 .2
T  mx (2-43)
2
 Energie potentielle
1
U  U p  U k  U k  kx2 (2-44)
 2
0

1 .2 1 2
L  m x  kx (2-45)
2 2
 Equation de dissipation
1 .2
D x (2-46)
2
 Equation de Lagrange
 
d  L  L D
dt   x.   x   .
 x
 
Après les dérivations nous recevons cette équation
.. .
m x kx   x (2-47)

Si nous devisons cette équation sur(m) nous acquérons l’équation de mouvement

..  . k
x x x0 (2-48)
m m
Si nous superposions cette équation et l’équation générale du mouvement :

𝑥̈ + 2𝜆𝑥̇ +𝜔02 𝑞 = 0 (2-49)

𝑘 
Nous obtenons 𝜔0 = √𝑚 et 
2m
2.2.3.1 Résolution de l’équation du mouvement
La résolution de l’équation du mouvement ( 2-40) sera supposée comme suit :
. ..
q(t )  Aert q(t )  Are rt q(t )  Ar 2 e rt (2-50)

Dr. Bennaceur Said Page 15


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

En remplaçant dans l’équation ( 2-26), on aboutit à :

Ar 2 e rt  2Are rt  0 Aert  0
2

(r 2  2r  0 ) Aert  0 r 2  2r  0  0


2 2
(2-51)

Donc il y a trois solutions suivant le sens de déterminant

  (b 2  ac )  2  0
2

  2  0  0 le mouvement apériodique (  0.5) :


2

Nous remarquons un grand amortissement imposé sur le système .En effet, ce système
une fois lâché de sa position d’équilibre ne fait que revenir à sa position d’équilibre sans
faire d’oscillation, tellement l’amortissement appliqué est fort.
La solution de notre équation précédente sera
comme suit :

q(t )  A1e r t  A2 e r t
1 2

𝑟1 = −𝜆 − √𝜆2 − 𝜔02

𝑟1 = −𝜆 + √𝜆2 − 𝜔02

−(𝜆+√𝜆2 −𝜔02 )𝑡 −(𝜆−√𝜆2 −𝜔02 )𝑡


𝑞(𝑡) = 𝐴1 𝑒 + 𝐴2 𝑒 (2-52)

  2  0  0 le mouvement critique (  0.5) :


2

Le régime est dit critique car il correspond à un amortissement critique pour lequel on
bascule d’un régime ou il n'y a plus d'oscillations vers pseudopériodique.

q(t )  A1e r t  A2te r t


1 2

r1  r2  r  

q(t )  ( A1  A2t )e t (2-53)

  2  0  0 le mouvement pseudo-périodique (  0.5) :


2

L'amplitude décroissante au cours du temps en raison de frottement.

Dr. Bennaceur Said Page 16


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

q(t )  A1e r t  A2 e r t
1 2

r1    j 0  2
2

r2    j 0  2
2

−(𝜆+𝑗√𝜔02 −𝜆2 )𝑡 −(𝜆−𝑗√𝜔02 −𝜆2 )𝑡


𝑞(𝑡) = 𝐴1 𝑒 + 𝐴2 𝑒 (2-54)

le calcul montre que q(t ) s'écrit sous la forme équivalente :

q(t )  A1e t cos(t   ) (2-55)

Pseudo-pulsation :   0 2  2 (2-56)

2𝜋 2𝜋
Pseudo –périodique : 𝑇= = (2-57)
𝜔
√𝜔02 −𝜆2

2.2.3.2 Décrément logarithmique


On définit le décrément logarithmique ( ) par :   ln
q(t n ) (2-58)
q(t n  T )

Où q(t n ) et q(t n  T ) représentent les amplitudes des oscillations aux instants

(t n ) et (t n  T ) ; généralement ces deux instants correspondant à deux extrema


successifs de même signe. Cette quantité mesure la décroissance des amplitudes.

Ae t cos(t n   )
n
Ae t n

  ln  (t T )  ln  (t T ) (2-59)
Ae cos( (t n  T )
n
Ae n

Ae t n

  ln  (t T )  ln e T
Ae n

  T (2-60)

Dr. Bennaceur Said Page 17


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

Exercices
Exercice 01 :
Soit le système mécanique vibratoire représenté sur la figure ci-contre. Si G est le centre de
gravité de la barre de masse M et de longueur L.
• Trouver l’équation différentielle du mouvement. Déduire  0 et 

• Ecrire l’équation du mouvement dans le cas   0

Exercice 02
Le système est constitué de 2 masses m1 et m2, d’une tige de masse négligeable de
longueur L et d’un ressort k et d’un amortisseur de
coefficient de frottement visqueux (α).
 Ecrire l’équation différentielle du du mouvement,
sachant que le système effectue des oscillations de
faible amplitude.
 Déterminer la pulsation propre du système.
 Trouver l’équation du mouvement, sachant
.
que θ(0)=θ0 et  (0)  0 . On donne :
m1 k
 m2  m ,  k et    0
4 4
Exercice 03
Dans son mouvement va-vient, le disque ci-dessous suspendu au ressort reste en contact
avec le mur qui le fait tourner sur lui-même. L’ensemble des frottements est symbolisé
par la force f disque  v appliquée au centre du disque. A l’équilibre le ressort était

allongé de y0.
1. Trouver l’énergie potentielle U de système en fonction
de y.
2. Simplifier U à l’aide de la condition d’équilibre.
3. Trouver l’énergie cinétique T et la fonction de
dissipation D.
4. Trouvez Lagrangien puis l’équation du mouvement.

Dr. Bennaceur Said Page 18


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

5. Sachant que k= 13.5N/m. M=1kg trouver la valeur maximale que le coefficient α


ne doit pas atteindre pour le système oscille.
6. Pour α = 9N.s/m, trouve la nature du mouvement ainsi que l’équation horaire y(t) (
.
y(0)=1cm, y(0)  0 .)

7. Si α = 3N.s/m, trouver le temps τ nécessaire pour que l’amplitude diminué à 1 de


5
sa valeur, calculer le décrément logarithmique δ du mouvement.

Exercice 04
La figure ci-contre représente un circuit série. Initialement,
S
le condensateur est chargé, puis on ferme
l’interrupteur S. On donne L=1H et C= 0.01 μF . R
C
on désigne par q(t) la charge du condensateur à L
l’instant t.
1. Ecrire l’équation différentielle qui décrit le circuit en fonction de la variable q. On
précisera les valeurs de λ et ω0.
2. Dans quel cas le système oscille ? quelle est la valeur de résistance critique ?
3. Tracer approximativement q(t) pour R=100Ω et R=500Ω.

Dr. Bennaceur Said Page 19


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

2.3 Oscillations amorties forcées


2.3.1 Définition d’un oscillateur amortie forcé
On appelle oscillateur amortie forcé lorsque qu’il est soumis à un frottement et une
force extérieure appelée force d’excitation.
2.3.2 Equation de Lagrange
Il y a deux formes de l’équation de Lagrange suivant le mode de déplacement :

 
d  L  L D
dt   q.   q   .  F (t ) Translation
(2-61)
 q
 

 
d  L  L D
dt   q.   q   .   (t ) Rotation
(2-62)
 q
 

Ou F(t) force d’excitation


μ(t) moment de la force d’excitation
2.3.3 Equation de mouvement
L’équation du mouvement des systèmes forcés est de la forme :
𝐹(𝑡)
𝑞̈ + 2𝜆𝑞̇ + 𝜔02 𝑞 = (a) Constante (2-63)
𝑎

Exemple
On appliquant dans le système décrit dans la section 2.2.3 une force extérieure
sinusoïdale:
F (t )  f 0 cos(t ) (2-64)

Figure 2-5 : Système amorti forcée (masse, ressort) horizontal

Dr. Bennaceur Said Page 20


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

 Lagrangien
L=T-U (2-65)
 Energie cinétique
1 .2
T  mx (2-66)
2
 Energie potentielle
1
U  U p  U k  U k  kx2 (2-67)
 2
0

1 .2 1 2
L  m x  kx (2-68)
2 2
 Equation de dissipation
1 .
D   x2 (2-69)
2
 Equation de Lagrange
 
d  L  L D (2-70)
dt   q.   q   .  F (t )
 q
 

Après les dérivations nous recevons cette équation


.. .
m x kx   x f 0 cos(t ) (2-71)
..  . k f cos(t )
x x x  0
m m m (2-72)
La constante a = m
Pour la pulsation propre et le coefficient d’amortissement
𝑘 
𝜔0 = √𝑚  (2-73)
2m
2.3.3.1 Résolution de l’équation du mouvement
D’après la figure (2-1) est présentée l’enregistrement du mouvement d’un oscillateur
amorite forcée nous remarquons deux régimes : le premier est court s’appelle transitoire
et le deuxième est long s’appelle permanent. Donc la solution de l’équation du
mouvement se compose de deux solutions
q(t)

Figure 2-6 : l’enregistrement du mouvement d’un oscillateur amorite forcé.

Dr. Bennaceur Said Page 21


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

q(t )  qT (t )  q P (t ) (2-74)
Ou : qT (t) solution transitoire
qP (t) solution permanente
Puisque : qT (t )  q P (t ) q (t )  q P t   A cos(t   ) (2-75)

(Ω) Pulsation de force excitatrice


(φ) Déphase entre le système et la force excitatrice
A l’aide la représentation complexe nous déduisons ces deux paramètres
Il faut noter q(t) est la partie réelle de 𝐴𝑒 𝑗(Ω𝑡+𝜑) : 𝑞(𝑡) = 𝑅𝑒[𝐴𝑒 𝑗(Ω𝑡+𝜑) ]

q (t )  q P t   A cos(t   )  A j ( t  )
(2-76)
F (t )  F0 cos(t )  F0 e j ( t )

𝑞̇ (𝑡) = 𝐴Ω𝑒 𝑗(Ω𝑡+𝜑) (2-77)

𝑞̈ (𝑡) = −𝐴Ω2 𝑒 𝑗(Ω𝑡+𝜑) (2-78)

j ( t   ) F0 ( jt )
( 2  2j   0 ) Ae 
2
e
a
F0
A( 2  2j   0 )  Be  j
2
On pose B=
a
A(( 2   0 )  2j )  Be  j
2

Nous utilisons le conjugué

A(( 2   0 )  2j )  Be  j
2

A(( 2   0 )  2j )  Be j
2

L’amplitude du mouvement est donc :


𝐹0
𝐴= 𝑎
(2-79)
2
√(𝜔02 −Ω2 ) +4𝜆2 Ω2

D’après l’équation

A(( 2   0 )  2j )  Be  j  cos   j sin 


2

cos   ( 2   0 )
2

sin   2

La phase 𝜑 du mouvement (déphasage entre 𝑞 (𝑡) et 𝐹 (𝑡) ) est donner par :

 2
tan   (2-80)
(  2   0 )
2

Dr. Bennaceur Said Page 22


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

2.3.3.2 Résonnance
La pulsation d’excitation  pour laquelle l’amplitude 𝐴 atteint son maximum est appelée
A
pulsation de résonance (d’amplitude)  R . 𝐴 est maximale lorsque 0


A
F0
 
4  0   2  8 2 
2
 
0 A a 0

  
3
2
2 0   2  4 2  2
2 2

A
0 −4(𝜔02 − Ω2𝑅 ) + 8𝜆2 = 0


 R   0  22
2
(2-81)

Si nous introduisons l’équation (2-81) dans l’équation ( 2-79), on obtient :


F0
Amax  (2-82)
a 4  02  4 4
2

0
Nous introduisons le facteur de qualité   dans l’équation (2-81)
2
1
Ω𝑅 = 𝜔0 √1 − 2𝜙 (2-83)

A cette pulsation, l’amplitude est :


F0 
Amax  (2-84)
a 2
1
0
1
4 2
1
Donc par la résonance il faut que :  0 2  22 0 1 − 2𝜙 > 0
1
∅>
√2
1
Le facteur de qualité doit donc être supérieur à ( 2) l’amortissement doit être faible

Figure 2-7 : évolution de d’amplitude en fonction la pulsation

Dr. Bennaceur Said Page 23


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

2.3.3.3 Résonance de phase



Lorsque Ω𝑅 = 𝜔0 tan    
2
lim tan( )  0 
 F(t) et q(t) en phase
0 2

lim tan( )  
𝜑=0 F(t) et q(t) en quadrature phase
  R

lim tan( )  0
   F(t) et q(t) en opposition phase


Figure 2-8 : évolution de la phase en fonction la pulsation


2.3.3.4 Bande passante
La bande passante est par définition l’ensemble des pulsations pour lesquelles l’amplitude
Amax
A
2
Nous définissions la largeur de la bande passante    2  1 ou  1 et  2
Amax
correspondent à A (1 ) =A(  2 )=
2
Si on pose  R   0 dans l’équilibre (2-82), on trouve
F0
Amzx  a (2-85)
2 0

Dr. Bennaceur Said Page 24


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

Figure 2-9 : évolution de d’amplitude en fonction la pulsation


𝐹0

42  2   02   2   82  02
𝐴𝑚𝑎𝑥 𝐹
= = 2𝑎𝜆𝜔0 √2
𝑎 2
√2 √(𝜔2 2 2 2 2 0
0 −Ω ) +4𝜆 Ω

Plus  est faible, plus le pic de résonnance est étroit (Ω petit) , on peut alors considéré
dans ΔΩ que Ω proche de ωo
D’où la bande passante de l’oscillateur : B     2  1  2 (2-86)
Donc le facteur de qualité sera comme suit :
0 0
  (2-87)
2 B
2.3.3.5 Analogie électromécanique
Nous pouvons faire une analogie entre le système mécanique (translation et rotation) et
le système électrique RLC est regroupé dans le tableau suivant :

Dr. Bennaceur Said Page 25


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

Tableau 2-1 : Analogie des paramètres électromécanique

Dr. Bennaceur Said Page 26


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

Exercices
Exercice 1
Une masse m, suspendue par un ressort de raideur k et un amortisseur de coefficient de
frottement, oscille verticalement sous l’effet d’une excitation F de la forme F (t)=F0cosΩ
1. Trouver l’énergie cinétique T, l’énergie potentielle U, et la fonction de dissipation
D.
2. Trouver le Lagrangien puis l’équation du mouvement.
3. Trouver, à l’aide de la représentation complexe, la solution
permanente de l’équation du mouvement. (Préciser son amplitude A
et sa phase)
4. Donner la condition de résonance et la pulsation de résonance ΩR:
5. Donner la bande passante B pour un amortissement faible:( λ˂˂ ω0)
Exercice 2
Dans le système ci-contre, la boule est ponctuelle et la tige est de longueur totale 3l et de
masse négligeable. Avec F (t)=F0cosΩ.
1. Trouver l’énergie cinétique T, l’énergie potentielle U, et la fonction de dissipation D.
2. Trouver le Lagrangien puis l’équation du mouvement.
3. Trouver, à l’aide de la représentation complexe, la
solution permanente de l’équation du mouvement.
(Préciser son amplitude A et sa phase)
4. Déduire la pulsation de résonance ΩR:
5. Donner les pulsations de coupure Ω1; Ω2 et la bande
passante B pour un amortissement faible ( λ˂˂ ω0).
6. Calculer ΩR , B, et le facteur de qualité si m=1kg,
k =15N/m, l=0,5m, α=0,5N.s/m, g=10m.s-2.
Exercice 3
Dans le système ci-contre, un fil autour du disque (de masse négligeable) est inextensible
et non glissant, est lié avec deux masses m1 , m2 et un ressort avec
un raideur.
1. Trouver l’énergie cinétique T, potentielle U, et la fonction de
dissipation D.
2. Trouver le Lagrangien puis l’équation du mouvement.
(F=F0cosΩt.)

Dr. Bennaceur Said Page 27


CHAPITRE 02 : Systèmes linéaires à un degré de liberté

3. Trouver, à l’aide de la représentation complexe, la solution permanente


de l’équation du mouvement. (Préciser son amplitude A et sa phase φ)
4. Déduire la pulsation de résonance ΩR.
5. Donner les pulsations de coupure Ωc1; Ωc2 et la bande passante B
Pour un amortissement faible ( λ˂˂ ω0).
6. Calculer ΩR; B; et le facteur de qualité si m1=2kg, m2=1kg, k=10N/m, α=0,1N.s/m.

Dr. Bennaceur Said Page 28


CHAPITRE 03 : Systèmes linéaires à plusieurs degrés de liberté

CHAPITRE 03 : Systèmes linéaires à plusieurs degrés de liberté

3.1 Systèmes linéaire de plusieurs degrés de liberté


Pour l’étude des systèmes à plusieurs degrés de liberté, il est nécessaire d’écrire plusieurs
équations différentielles du mouvement que l’on peut obtenir à partir des équations de
Lagrange :

(3-1)

3.1.1 Les types de couplages


a) Couplage Elastique : Le couplage dans les systèmes mécaniques sera par élasticité.
Pour les systèmes électriques, les circuits couplés sont présentés par capacité, ce qui est
équivalent au couplage par élasticité.

Figure 3-1 : Couplage systèmes élastique et équivalant électrique.

b) Couplage Visqueux : Le couplage dans les systèmes mécaniques est évident par
amortisseur. Dans les systèmes électriques, on aperçoit les circuits couplés par résistance,
équivalents au couplage par amortisseur

Figure 3-2 : couplage systèmes visqueux et équivalant électrique.

Dr. Bennaceur Said Page 31


CHAPITRE 03 : Systèmes linéaires à plusieurs degrés de liberté

c) Couplage Inertiel : Le couplage dans les systèmes mécaniques est assuré par inertie.
Dans les systèmes électriques, on trouve les circuits couplés par inductance, équivalents
au couplage par inertie

Figure 3-3 : couplage systèmes inertiels et équivalant électrique

3.2 Systèmes linéaire deux degrés de liberté cas libre


Un système à 2 degrés de liberté possède 02 coordonnées généralisées, 02 équations
différentielles et 02 pulsations propres (ω1, 𝜔2).

Figure 3-4 : Mouvement oscillatoire d’un système couplé à deux degré de liberté.
3.2.1 Equation du mouvement
Exemple
Soit le système libre décrit précédemment. Les deux variables indépendantes sont 𝑥1 et
𝑥2. 𝑘0 est appelé élément de couplage.
Calculer les énergies cinétique et potentielle du système, en utilisant la formule de
Lagrange établir les équations différentielles du mouvement. En déduire les pulsations
propres du système.
 Lagrangien
L=T-U (3-2)
 Energie cinétique
. .
1 1
T m x12  m x 2 2 (3-3)
2 2

Dr. Bennaceur Said Page 32


CHAPITRE 03 : Systèmes linéaires à plusieurs degrés de liberté

 Energie potentielle
1 2 1 1
U kx1  k 0 ( x1  x 2 ) 2  kx22 (3-4)
2 2 2
. .
1 1 1 1 1
L m x12  m x 2 2  kx12  k0 ( x1  x2 )2  kx22 (3-5)
2 2 2 2 2

 Equations de Lagrange

d  L  L ..
 0 m1 x1  (k1  k 0 ) x1  k 0 x 2  0
dt   x.  x1
 1 (3-6)
d  L  L ..
m 2 x 2  (k 2  k 0 ) x 2  k 0 x1  0
 0
dt   x.  x2
 2

j ( t 1 )
x1 (t )  A1 cos(t  1 )  A1
j (t  2 )
x 2 (t )  A2 cos(t   2 )  A2 (3-7)

Nous introduisons
𝑘1 +𝑘0 𝑘
(−𝜔2 + ) 𝐴1 − 𝑚0 𝐴2 = 0 (−𝜔2 + 𝑎)𝐴1 − 𝑏𝐴2 = 0
𝑚1 1
{ 𝑘 ⟹{ (3-8)
𝑘 +𝑘
− 𝑚0 𝐴1 + (−𝜔2 + 2𝑚 0 ) 𝐴2 = 0 −𝑐𝐴1 + (−𝜔2 + 𝑑)𝐴2 = 0
2 2

Ces équations admettent des solutions si le déterminant sera nul

(3-9)
3.2.2 Les modes d’oscillations
Définition du mode d’oscillation : le mode d’oscillation est l’état dans lequel les
éléments dynamiques du système effectuent une oscillation harmonique avec la même
pulsation qui correspond à une de ses deux pulsations.
Les deux solutions réelles et positives 𝜔1 et 𝜔2 de cette équation sont appelées pulsations
propres ou normales. La plus petite est appelée la fondamentale, l’autre est appelée
l’harmonique.
Pour notre cas nous présentons les deux modes comme suit :

Premier mode propre : Pour 𝜔=𝜔1 : La vibration est


en phase

Dr. Bennaceur Said Page 33


CHAPITRE 03 : Systèmes linéaires à plusieurs degrés de liberté

(3-10)

Deuxième mode propre : Pour 𝜔=𝜔2 : La vibration est en


opposition de phase

(3-11)

3.3 Système forces et amorti a deux degrés de libertés

Figure 3-5 : Mouvement oscillatoire d’un système amorti et forcé


couplé à deux degré de liberté.
Exemple : en utilisant l’équation de Lagrange déduire les équations du mouvement du
système de la figure 3-5 qui se compose de deux masses m1 et m2 reliées par 03 ressorts
ces constants de raideurs k1, k2 et k0 glissent sur un plan horizontal. Sont soumises à deux
forces de frottements représentées par la présence d’amortisseur des coefficients
d’amortissement identiques (α) et une force externe dépondant du temps
 Lagrangien

L=T-U (3-12)

 Energie cinétique et Energie potentielle


. .
1 1
T m1 x12  m2 x 22 (3-13)
2 2
1 1 1
U  kx12  k 0 ( x1  x 2 ) 2  kx22 (3-14)
2 2 2

. .
1 1 1 1 1
L m1 x12  m2 x 2 2  k1 x12  k 0 ( x1  x 2 ) 2  k 2 x 22 (3-15)
2 2 2 2 2
. .
1 1 F (t )  F0 cos(t )
D   1 x1   2 x 2 et
2 2

2 2
 
d  L  L D
m1 x1  (k1  k 0 ) x1   1 x1  k 0 x 2  F0 cos t 
.. .
   .  F1 t 
dt   x.  x1  x1
 1 (3-16)
  .. .
d  L  L D m2 x 2  (k 2  k 0 ) x 2   2 x 2  k 0 x1  0
    .
dt   x  x2
.
 x2
 2

Dr. Bennaceur Said Page 34


CHAPITRE 03 : Systèmes linéaires à plusieurs degrés de liberté

3.4 Oscillations forcées d’un système non amorti à deux degrés de liberté
Nous considérons notre système mécanique symétrique décrit dans la figure 3-5.
(𝑘1 = 𝑘2 = 𝑘 𝑒𝑡 𝑚1 = 𝑚2 = 𝑚) ou nous appliquons une force extérieure de forme
sinusoïdale au premier sous système qui s’exprime comme suit :

m1 x1  (k1  k 0 ) x1  k 0 x 2  F0 cos t 


..

(3-17)
D=0 et F≠0
..
m2 x 2  (k 2  k 0 ) x 2  k 0 x1  0
Régime permanant
j ( t 1 )
x1 (t )  A1 cos(t  1 )  A1
j (t )
j (t  2 ) x1 (t )  A1
x2 (t )  A2 cos(t   2 )  A2 (3-18)
j (t )
j (t ) x 2 (t )  A 2
F (t )  F0 cos(t )  F0 e

(3-19)

(3-20)

Figure 3-6 : Evolution de l’amplitude en fonction de la pulsation

Dr. Bennaceur Said Page 35


CHAPITRE 03 : Systèmes linéaires à plusieurs degrés de liberté

3.5 Généralisation aux systèmes à n degrés de liberté


Energie cinétique généralisée
Un système à n degrés de libertés possède n variables : x1, x2, x3…………… xn.
Énergie cinétique généralisée

1 .
1 .
1 . 2
T m1 x12  m2 x 22  ................  mn x n
2 2 2
Si on suppose : m1  m2  .......  mn  m

1 . . .
1 n .2
T m( x12  x 22  ...........  x n2 )  m x i (3-21)
2 2 i 1
𝒙𝟏
𝒙
On pose vecteur vertical (colonne) |𝒙 > = ( 𝟐.. ) ..
.
𝒙𝒏
On pose vecteur horizontal (ligne) < 𝒙| = (𝒙𝟏 , 𝒙𝟏 , … … . 𝒙𝒏 )

𝒙̇ 𝟏
𝒙𝟏 𝒏 𝒙̇ 𝟐 𝒏
𝒙
⟨𝒙|𝒙⟩ = (𝒙𝟏 , 𝒙𝟐 , … … . 𝒙𝒏 ) ( 𝟐.. ) = ∑ 𝒙𝒊 ↔ ⟨𝒙̇ |𝒙̇ ⟩ = (𝒙̇ 𝟏 , 𝒙̇ 𝟐 … . . 𝒙̇ 𝒏 ) . ↔ ∑ 𝒙̇ 𝟐𝒊
𝟐
.
𝒙𝒏 𝒊=𝟏 . 𝒊=𝟏
𝒙̇
( 𝒏)
𝒎 𝟎 𝟎…… 𝟎 𝒙̇ 𝟏
𝟎 𝒎 𝟎…… 𝟎 𝒙̇ 𝟐
𝟏 𝟏 𝟏
𝑻 = 𝟐 𝒎⟨𝒙̇ |𝒙̇ ⟩= 𝟐 (𝒙̇ 𝟏 , 𝒙̇ 𝟐 … . . 𝒙̇ 𝒏 ) 𝟎 𝟎 𝒎
. … . . 𝟎 . 𝑻 = 𝟐 ⟨𝒙̇ |𝓣|𝒙̇ ⟩
. .
( 0 0 0 ……..m ) (𝒙̇ 𝒏 ) (3-22)

𝓣 : Matrice carrée (nxn), Opérateur associé à T. Les éléments de 𝝉 sont déduits des
dérivées

Énergie potentielle généralisée


𝟏 𝟏 𝟏 𝟏
𝑼 = 𝟐 𝒌𝟏 𝒙𝟐𝟏 + 𝟐 𝒌𝟐 𝒙𝟐𝟐 + 𝟐 𝒌𝟑 𝒙𝟐𝟑 + … … … … + 𝟐 𝒌𝒏 𝒙𝟐𝒏 (3-23)

Si k1= k2= k1=……………= kn.=k


𝟏 𝟏
𝑼 = 𝟐 𝒌(𝒙𝟐𝟏 + 𝒙𝟐𝟐 + 𝒙𝟐𝟑 + ⋯ 𝒙𝟐𝒏 ) = 𝟐 𝒌 ∑𝒏𝒊=𝟏 𝒙𝟐𝒊 (3-24)

𝒌 𝟎 𝟎…… 𝟎
𝟎 𝒌 𝟎…… 𝟎
𝟏 𝟏
𝑼 = 𝟐 𝒌⟨𝒙|𝒙⟩ = 𝑼 = 𝟐 ⟨𝒙|𝝁|𝒙⟩ 𝒕𝒆𝒍 𝒒𝒖𝒆 µ = 𝟎 𝟎 𝒌. … . . 𝟎 (3-25)
0 0 0. …….k
( )

Dr. Bennaceur Said Page 36


CHAPITRE 03 : Systèmes linéaires à plusieurs degrés de liberté

Le Lagrangien (la fonction de Lagrange) :


1 1
L=T-U= 2 ⟨𝑥̇ |𝒯|𝑥̇ ⟩ − 2 ⟨𝑥|𝜇|𝑥⟩ (3-26)
𝑑 𝜕𝐿 𝜕𝐿
Equation de Lagrange : ( ) − (𝜕𝑥 ) = 0
𝑑𝑡 𝜕𝑥̇ 𝑖 𝑖

𝜕𝐿 1 𝜕
𝜕𝑥̇ 𝑖
𝒯
= 2 𝜕𝑥̇ ⟨𝑥̇ | |𝑥̇ ⟩
𝑖
{ 𝜕𝐿 1 𝜕
(3-27)
𝜕𝑥𝑖
= − 2 𝜕𝑥 ⟨𝑥|𝜇|𝑥⟩
𝑖

𝜕 𝜕
On a : 𝜕𝑥̇ ⟨𝑥̇ |𝒯|𝑥̇ ⟩ = 2⟨𝐼𝑖 |𝒯|𝑥̇ ⟩ 𝑒𝑡 ⟨𝑥|𝜇|𝑥⟩ = 2⟨𝐼𝑖 |µ|𝑥̇ ⟩ avec Ii vecteur unité.
𝑖 𝜕𝑥𝑖

⟨𝐼1 ⁄= (1,0,0,0 … … 0) , ⟨𝐼2 ⁄= (1,0,0,0 … … 0) , … … … … . ⟨𝐼𝑛 ⁄= (1,0,0,0 … … 0).


𝜕
⟨𝐼 |
𝜕𝑥̇ 𝑖 𝑖
𝒯|𝑥̇ ⟩ ⟹ 𝑑𝑡𝑑 (𝜕𝑥
𝜕𝐿
̇
) = ⟨𝐼𝑖 |𝒯|𝑥̇ ⟩ (3-28)
𝑖

𝜕𝐿
= −⟨𝐼𝑖 |µ|𝑥⟩ ⟹ ⟨𝐼𝑖 |𝒯|𝑥̈ ⟩ + ⟨𝐼𝑖 |𝜇|𝑥⟩ = 0, c ′ est l′ équation de Lagrange
𝜕𝑥𝑖
𝑑2
⟹ 𝒯/𝑥̈ ⟩ + 𝒰/𝑥⟩ = 0 ⟹ /𝑥̈ ⟩ + 𝒯 −1 = 0 ⟹ 𝑑𝑡 2 /𝑥⟩ + ℒ/𝑥⟩ = 0 (3-29)

ℒ = 𝒯 −1 . 𝜇: Opérateur associé au Lagrange ℒ, 𝒯 −1 : Matrice inverse de 𝒯

Equation aux valeurs propres:


𝑑2
La solution de l’équation /𝑥⟩ + ℒ/𝑥⟩ = 0 peut être sous la forme complexe :𝑪𝒆𝒋𝝎𝒕
𝑑𝑡 2

𝑑2
/𝑥⟩ = −𝜔2 /𝑥⟩ ⟹ −𝜔2 /𝑥⟩ + ℒ/𝑥⟩ = 0 ⟹ (ℒ − 𝜔2 )/𝑥⟩ = 0 , (3-30)
𝑑𝑡 2

(ℒ − 𝜔2 )/𝑥⟩ = 0 ⟹ 𝐷𝑒𝑡 [ℒ − 𝜔2 𝐼] = 0 c’est l’équation aux valeurs propres. 𝜔𝑖 :


valeur propre.

 Valeur propres :
Exemple : ℒ = ( 𝐿𝐿11 𝐿12
𝐿22
),
21

2
L’équation aux valeurs propres : ℒ − 𝜔2 𝐼 = ( 𝐿𝐿11 𝐿12
𝐿22
) − 𝜔2 ( 01 0
1
) = ( 𝐿𝐿11 −𝜔𝐿 𝐿12
2 )
21 21 22 −𝜔

𝐿 −𝜔2 𝐿12
| 𝐿 11 𝐿22 −𝜔2
| = 0 ⟹ (𝐿11 − 𝜔2 )( 𝐿22 − 𝜔2 ) − 𝐿21 . 𝐿21 = 0 (3-31)
21

⃗⃗𝑖 , 𝜔1 ↦ ⃗⃗⃗
A chaque valeur propre 𝜔𝑖 correspond un vecteur𝑉 𝑉1 (𝑥𝑥11 ) , 𝜔2 ↦ ⃗⃗⃗
𝑉2 (𝑥𝑥21 )
12 22

⃗⃗⃗1 = 0 ⟹ | 11 𝐿 −𝜔12 𝐿12


 Pour : 𝜔 = 𝜔1 : (ℒ − 𝜔12 )𝑉 𝐿
| (𝑥11 )
𝐿22 −𝜔12 𝑥12
= (00) (3-32)
21

Dr. Bennaceur Said Page 37


CHAPITRE 03 : Systèmes linéaires à plusieurs degrés de liberté

(𝐿11 − 𝜔21 )𝑥11 + 𝑥12 𝐿12 = 0


⟹{ ⟹ 𝑥11 𝑒𝑡 𝑥12 ? (3-33)
𝐿21 𝑥11 − (𝐿22 − 𝜔21 )𝑥12 = 0

⃗⃗⃗2 = 0 ⟹ | 11 𝐿 −𝜔22 𝐿12


 Pour : 𝜔 = 𝜔2 : (ℒ − 𝜔22 )𝑉 𝐿
| (𝑥12 )
𝐿22 −𝜔22 𝑥22
= (00) (3-34)
21
⟹ 𝑥12 𝑒𝑡 𝑥22

Solutions des équations différentielles :

𝑥1(𝑡) 𝑥 𝑥
( ) = 𝐴𝑉 ⃗⃗⃗2 𝑒 𝑗(𝜔2 𝑡+𝜑2) = 𝐴 ( 11 ) 𝑒 𝑗(𝜔1 𝑡+𝜑1 ) + 𝐵 ( 21 ) 𝑒 𝑗(𝜔2 𝑡+𝜑2 )
⃗⃗⃗1 𝑒 𝑗(𝜔1 𝑡+𝜑1) + 𝐵𝑉
𝑥2 (𝑡) 𝑥12 𝑥22

𝑥1(𝑡) = 𝐴𝑥11 𝑒𝑗(𝜔1𝑡+𝜑1) + 𝐵𝑥21 𝑒𝑗(𝜔2𝑡+𝜑2)


{ (3-35)
𝑥2(𝑡) = 𝐴𝑥12 𝑒𝑗(𝜔1𝑡+𝜑1) + 𝐵𝑥22 𝑒𝑗(𝜔2𝑡+𝜑2)

Dr. Bennaceur Said Page 38


CHAPITRE 03 : Systèmes linéaires à plusieurs degrés de liberté

Exercices
Exercice 01 :
On considère les oscillations libres du système à deux degrés de liberté
1. Calculer les énergies cinétique et potentielle du système ;
𝑚
2. Pour k1 = k2 = k et m1= m2= , et en utilisant la formule de Lagrange établir les
2

équations différentielles du mouvement. En déduire les pulsations propres du


système.

Exercice 02 :

On considère deux circuits électriques ( R , Lind ,Cap ) couplés représenté par la figure
ci-dessous :

 Quel est le nombre de degré de liberté ?


 Déterminer le Lagrangien du système.
 Donner les équations du mouvement
On néglige les résistances des deux circuits. On prend les nouvelles grandeurs
1
Physiques tel que: L1ind = L2ind = Lind et C1ap= C2ap = Cap et 𝜔02 =
𝐿𝑖𝑛𝑑 𝐶𝑎𝑝

 Etablir les nouvelles équations différentielles du mouvement.


 En déduire les pulsations propres du système en fonction de 𝜔0
 Donner les solutions générales.
 Quel est le modèle mécanique équivalent ?

Dr. Bennaceur Said Page 39


CHAPITRE 03 : Systèmes linéaires à plusieurs degrés de liberté

Exercice 03 :

Sur un arbre OO’ horizontal et fixe, de masse négligeable, encastré à ses


extrémités O et O’, sont fixés trois disques (D1), (D2) et (D3) de centres respectifs
O1, O2 et O3 et de même moment d’inertie J par rapport à leur axe commun OO’.
On désignera 𝜃1(t) 𝜃2(t) et 𝜃3(t), les angles angulaires de rotation de chacun des
trois disques par rapport à leur position de repos

de liberté

Les quatre partis OO1, O1O2, O2O3 et O3O’de l’arbre ont la même constante de
torsion C.
𝐶
On posera 𝜔02 =
𝐽
Régime libre :
 Déterminer le Lagrangien de ce système.
 Etablir les équations différentielles du second ordre vérifiées par les angles
𝜃1(t), 𝜃2(t) et 𝜃3(t).
 En déduire les trois pulsations propres 𝜔1p, 𝜔2p et 𝜔3p de ce système en
fonction de 𝜔0.
 Déterminer pour chaque des trois modes propres, les amplitudes angulaires
des disques D2 et D3 si l’amplitude angulaire du disque D1 est A= 1 radian.
 Calculer l’énergie mécanique totale ET de cette chaîne de trois disques, pour
chacun des modes propres, en fonction de C et de l’amplitude angulaire 𝜃1
du disque D1.
Régime forcé :
On applique au seul disque (D1) un couple moteur sinusoïdal de moment
Γ(𝑡) = Γ0 cos(𝜔𝑡) La pulsation est réglable et d’amplitude Γ0 .
𝜔 2
 Etablir en fonction du paramètre 𝑋 = (𝜔 ) les amplitudes angulaires A1, A2 et A3
0

de chacun des disques en régime forcé.


 Pour quelles valeurs de X ce système est-il en résonance ?

Dr. Bennaceur Said Page 40


Partie II
LES ONDES
MECANIQUES

Dr. Bennaceur Said Page 40


CHAPITRE 04 : Généralités sur les ondes mécaniques

CHAPITRE 04 : Généralités sur les ondes mécaniques


4.1 Classification des ondes
Une onde est un phénomène vibratoire qui se propage. Elle peut être classée en deux
familles: ondes mécaniques et ondes électromagnétiques.
Ces ondes ont deux principaux modes de propagation.
a- ondes mécaniques (élastiques)
• onde acoustique (son)
• onde sismique
• onde à la surface d’un liquide
• onde dans une corde
b- ondes électromagnétiques
• onde lumineuse (lumière)
• onde radio
Et avec deux principaux modes de propagation :
a- mode longitudinale (onde L)
La vibration est parallèle au sens de propagation (onde de “compression”)
Exemple : propagation de la déformation d’un ressort

Figure 4-1 : propagation dans le ressort


b- mode transversale (onde T)
La vibration est perpendiculaire au sens de propagation (onde de “cisaillement”).
• Exemple : propagation de la déformation d’une corde

Figure 4-2 : propagation dans une corde


Dr. Bennaceur Said Page 41
CHAPITRE 04 : Généralités sur les ondes mécaniques

4.2 Intégrale générale de l’équation générale d’ondes progressives.


Dans les phénomènes vibratoires traités dans les chapitres précédents, nous nous sommes
intéressés à des phénomènes ou des grandeurs physiques qui dépendaient d'une seule
variable, le temps t. Nous allons maintenant examiner toute une série de phénomènes qui
sont décrits par une fonction qui dépend à la fois du temps t et d'une variable d'espace, x
par exemple. Ces phénomènes sont régis par une équation aux dérivées partielles, appelée
équation d'Alembert ou équation d'onde ou encore équation de propagation à une
dimension de la forme :
2A 1 2A
 0 (4-1)
x 2 V 2 t 2
A : la fonction d’onde qui se propage
x : la direction de propagation
t : temps
Dans laquelle V est une grandeur physique qui a les dimensions d'une vitesse et sera
appelée dans la suite vitesse de propagation.
4.2.1 Solution de l'équation de propagation
L’équation d’onde simple (4-1) peut se décomposer sous la forme :
 A 1 A  A 1 A 
    0 (4-2)
 x V t  x V t 
  1    1  
Ou     A  0 (4-3)
 x V t  x V t 

La dernière formulation de l’équation d’onde montre qu’il existe deux vitesses dans le
système : +V et −V, associées à deux sens de propagation opposés. Pour chercher des
solutions à l’équation d’onde, on introduit deux nouvelles variables :
X (x , t )  x Vt
 (4-4)
Y (x , t )  x Vt
qui sont des variables d’espace associées à deux référentiels mobiles aux vitesses +V et –
V respectivement. Dans l’écriture en (x, t), le référentiel a une position fixe, tandis que
les variables X et Y sont associés à deux référentiels mobiles.
Les opérateurs de dérivation partielle de l’équation d’onde (4-1) doivent être transformés
et on utilise pour cela :

Dr. Bennaceur Said Page 42


CHAPITRE 04 : Généralités sur les ondes mécaniques

  X  Y  
    (4-5)
x X x Y x X Y
De même, la transformation de l’opérateur ∂/∂t est :
  X  Y    
  V    (4.6)
t X t Y t  Y X 
Dans l’équation (4-1), nous avons besoin de l’opérateur de dérivée seconde en x, donc
∂2/∂x2 :

2       
2 2

     ,
x 2
 x   X Y 
et l’opérateur ∂2/∂t2 est

2       
2 2

   V 2   
t 2
 t   Y X 
Par ce changement de variables, l’équation (4.1) devient une équation intégrable :
2A   A    A 
    0 (4.7)
X Y X  Y  Y  X 
forme 1 forme 2 forme 3

La deuxième forme de cette équation montre que la dérivée de A par rapport à Y ne


dépend pas de X, donc la fonction A ne dépend que de Y . Mais la troisième forme
montre que la dérivée de A par rapport à X ne dépend pas de Y , donc la fonction A ne
dépend que de X . La forme générale de la solution A(X , Y ) est donc la somme de deux
fonctions indépendantes qui ne dépendent chacune que d’une variable, X ou Y :
A  x , t   F  X   G Y   F  x Vt   G  x Vt  (4-8)

Les deux fonctions F et G sont des fonctions d’onde arbitraires et ne sont déterminées que
par les conditions aux limites et initiales du système. L’équation d’onde simple n’impose
aucune forme particulière pour ces fonctions, du moment qu’elles sont dérivables.
4.3 Vitesse de phase, vitesse de groupe
La vitesse de phase caractérise la vitesse de déplacement d'une onde plane. Elle dépend
éventuellement de la fréquence si le milieu est dispersif. En ce sens cette vitesse n'a pas
vraiment de sens physique car l'onde plane n'est qu'une représentation mathématique utile
du paquet d'ondes qui lui est un vrai signal physique. La vitesse de groupe caractérise le
déplacement de l'ensemble du paquet, c'est si on veut une sorte de vitesse moyenne du
groupe d'ondes planes qui définit le signal.

Dr. Bennaceur Said Page 43


CHAPITRE 04 : Généralités sur les ondes mécaniques

Vitesse de phase Par simple géométrie, il est clair que la crête de l’onde (son extremum)
avance d’une longueur λ à chaque période τ et donc la vitesse de la crête, aussi appelée
vitesse de phase, est :
υp = λ τ = λν (4-9)
La vitesse de la crête de l’onde, ou vitesse de phase.
On peut écrire : v = υp / λ ou λ = υp / v (4-10)
On peut dire de υp = λ ν que la crête avance d’une distance λ , v fois par unité de temps.
Cette relation pour (υp) est valable pour toute onde, parce qu’elle est purement
géométrique. Selon le mécanisme physique responsable de la génération de cette onde, la
vitesse (υp) peut être ou ne pas être une fonction de la longueur d’onde/fréquence ou peut
être constante et indépendante de la longueur d’onde/fréquence. Dans ce dernier cas,
l’onde est dite harmonique. Nous y reviendrons.
4.4 Front d’onde
Nous pouvons aussi appeler surface d'onde ou nous définissons comme
une surface d'égale phase d'une onde, c'est-à-dire que ces points ont mis le même temps
de parcours depuis la source.

Front d’onde

Figure 4-3 : Surface d’onde (front d’onde)


4.4.1 Exemple des ondes planes, ondes sphériques
 Onde plane
Une onde est plane quand les fronts d’ondes sont des surfaces planes. Cela correspond à
⃗ orienté suivant une seule direction de l’espace. Les fronts d’ondes
un vecteur d’onde 𝑘
sont des plans parallèles entre eux séparés d’une distance 𝜆 = 2𝜋/𝑘 et tous
⃗ . Si le vecteur d’onde est orienté suivant la direction x , le
perpendiculaires à la direction𝑘
⃗ = 𝑒𝑥 et l’onde plane monochromatique est notée :
vecteur propagation est 𝑘
𝐴(𝑥 , 𝑡) = 𝐴0 𝑐𝑜𝑠(𝑘𝑥 − 𝑣𝑡) (4-11)

Dr. Bennaceur Said Page 44


CHAPITRE 04 : Généralités sur les ondes mécaniques

Figure 4-4 : Illustrations des fronts d’onde plane (a) et de l’amplitude d’une onde se
propageant selon l’axe x (b).

 Onde sphérique
Quand la propagation admet une symétrie ponctuelle dans un espace isotrope à trois
dimensions, l’équation d’onde admet comme solutions les combinaisons linéaires des
ondes monochromatiques sphériques. La direction de propagation est donnée par :
⃗ = 𝑘𝑟 𝑒𝑟 ,
𝑘 𝑘𝜃 = 0, 𝑘𝜑 = 0
Et l’opérateur laplacien en coordonnées sphériques est tel que :
1 𝜕2
Δ𝐴 = 𝑟 𝜕𝑟 2 (𝑟𝐴) (4-12)

L’équation d’onde (4-12) s’écrit donc


1 𝜕2 1 𝜕2 𝐴
(𝑟𝐴) − =0 (4-13)
𝑟 𝜕𝑟 2 𝑐 2 𝜕𝑡 2

On retrouve ici une équation d’onde à une dimension d’espace dont la solution r A se
décompose sous la forme d’une somme F (r − ct) + G(r + ct). La fonction d’onde
sphérique est donc de la forme
1 1
A(r, t)= 𝐴(𝑟, 𝑡) = 𝑟 𝐹(𝑟 − 𝑐𝑡) + 𝑟 𝐺(𝑟 + 𝑐𝑡) (4-14)

Et l’onde monochromatique sphérique s’écrit sous la forme


𝐴0
𝐴(𝑟, 𝑡) = 𝑟0 ‖ ± 𝜔𝑡)
𝑐𝑜𝑠(𝑘‖𝑟 − ⃗⃗⃗ (4-15)
𝑟

Avec 𝑟0 = (𝑟0 , 𝜃0 , 𝜑0 ) la position de la source ou du point d’absorption et 𝑟 = (𝑟,


𝜃, 𝜑) les coordonnées sphériques du point considéré.

Dr. Bennaceur Said Page 45


CHAPITRE 04 : Généralités sur les ondes mécaniques

Figure 4-5 : Illustrations des fronts d’onde sphérique (a) et de


l’amplitude de l’onde sphérique (b).
4.5 Réflexion et transmission des ondes
Le dioptre est le plan xOz. Le demi-espace correspondant à y > 0 est rempli d’un milieu 1
d’indice optique n1 le demi espace correspondant à y < 0 est empli d’un milieu d’indice
optique n2.Nous considérons deux demi-espaces séparés par un dioptre plan.

Figure 4-5 : Réflexion et transmission


⃗ i se dirigeant vers le
Nous considérons une onde plane progressive de vecteur d’onde 𝑘
⃗ i fait un angle θi (angle d’incidence) avec l’axe 0y qui est
dioptre. Le vecteur d’onde 𝑘
normal à l’inferface. A l’interface entre les deux milieux, on constate deux phénomènes :
– Après la traversée du dioptre, l’onde est déviée. C’est le phénomène de réfraction.
⃗ t fait un
L’onde transmise est toujours une onde plane progressive. Son vecteur d’onde 𝑘
angle θt avec l’axe Oy.
– Une onde plane est renvoyée par le dioptre vers le milieu 1. C’est le phénomène de
réflexion. Cette onde réfléchie est elle aussi une onde plane progressive. Son vecteur
⃗ r fait un angle θr avec l’axe Oy.
d’onde est 𝑘
Notons immédiatement que les pulsations de ces trois ondes sont identiques.

Dr. Bennaceur Said Page 46


CHAPITRE 04 : Généralités sur les ondes mécaniques

Exemple : Le champ électromagnétique dans les deux demi-espaces.


Dans le demi-espace y > 0 le champ électrique 𝜀 (x, y > 0, z, t) est la superposition des
champs des ondes incidente et réfléchie :

𝜀𝑖0 𝑒 𝑖(𝑘⃗⃗⃗𝑖 .𝑟−𝜔𝑡) + ⃗⃗⃗⃗⃗


𝜀 (x, y > 0, z, t) = ⃗⃗⃗⃗ 𝜀𝑟0 𝑒 𝑖(𝑘⃗⃗⃗⃗𝑟 .𝑟−𝜔𝑡) . (4-16)
Tandis que dans le demi espace y < 0, une seule onde est présente. Le champ électrique
𝜀 (x, y< 0, z, t) est donc celui de l’onde transmise :

𝜀𝑖0 𝑒 𝑖(𝑘⃗⃗⃗𝑖 .𝑟−𝜔𝑡) .


𝜀 (x, y < 0, z, t) = ⃗⃗⃗⃗ (4-17)
Comme chacune des trois ondes (incidente, réfléchie et transmise) est une onde plane, le
champ magnétique s’exprime simplement en fonction des vecteurs d’onde et du champ
électrique. En notant α l’indice qui correspond à chacune des ondes (α = i, r ou t) :
⃗⃗⃗⃗⃗
𝑘 𝑛 ⃗⃗⃗⃗⃗
𝑘
𝛽𝑖0 𝑒 𝑖(𝑘⃗⃗⃗⃗⃗𝛼.𝑟−𝜔𝑡) = 𝜔𝛼 ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝛽𝛼 = ⃗⃗⃗⃗⃗ × 𝜀𝑖0 𝑒 𝑖(𝑘⃗⃗⃗𝑖 .𝑟−𝜔𝑡) = 𝑐 ‖𝑘⃗⃗⃗⃗⃗𝛼 ‖ × ⃗⃗⃗⃗⃗⃗
𝜀𝛼0 𝑒 𝑖(𝑘⃗⃗⃗⃗⃗𝛼.𝑟−𝜔𝑡) .
𝛼

Dr. Bennaceur Said Page 47


CHAPITRE 04 : Généralités sur les ondes mécaniques

Exercices
Exercice 1 :
Vérifier que les fonctions suivantes sont des solutions de l’équation D’Alembert
x
x i  (t  )
a) u (x , t )  A sin( (t  )) , b) u (x , t )  Ae V
, c) u (x , t )  A cos(k (x Vt ) , d)
V
x
i  (t  )
u (x , t )  Ae V
e) u (x , t )   (x Vt )2
On suppose que x, t et V sont respectivement la position, le temps et la vitesse de
propagation
Exercice 02 :
Soit des ondes sinusoïdales définies par les deux équations suivantes :
A1 (x , t )  2cos(6.28t 15x ) et A2 (x , t )  5cos(9.42t  22x )
Quelle est la vitesse de phase υp pour chaque onde
1  2 1  2
Sachant que cos(1 )  cos(2 )  2cos( ) cos( ) trouver la vitesse de phase υp et
2 2
de groupe VG de l’onde résultante
Exercice 03 :
La figure ci-dessous représente l’aspect d'un ébranlement se propageant le long d'une corde
indéfiniment longue. L'aspect de la corde est représenté aux instants t =2 s et t =5 s.

1- Calculer la vitesse de propagation de cet ébranlement. Quelle est la position de


l'ébranlement à t = 0?
2- A quel instant t1 le point B aura-t- il une amplitude maximale ?
3- A quel instant t2 le point B revient-il à sa position initiale ?
4- A quel instant t3 le point M retrouvera-t-il pour la première fois une élongation nulle ?
5- Tracer sur un même graphe les amplitudes yA(t) et yB(t) des points A et B de la corde.
6- Représenter la courbe de vitesse 𝑦̇ 𝐵 (𝑡) du point B.

Dr. Bennaceur Said Page 48


CHAPITRE 04 : Généralités sur les ondes mécaniques

Exercice 04 :
Une onde transversale sinusoïdale, se propage le long d’une corde dans le sens des x
négatifs à la vitesse V = 10m/s. La figure ci-dessous illustre le déplacement des particules
sur la corde en fonction de la position à un instant t.

Calculer :
a) l’amplitude de cette onde, b) sa longueur d’onde, c) sa vitesse de phase, d) sa période,
e) la vitesse maximale d’une particule sur la corde. f) Quelle est l’équation de cette onde ?

Exercice 05 :
1- Deux ébranlements de même amplitude mais de signes opposés se dirigent l’un vers
l’autre sur la corde. La figure ci-contre les représente à t=0 . Si les deux ébranlements sont
distants de d = 60cm, dessinez la configuration de la corde aux instants 10, 15 et 30 ms.
2- Répondre aux mêmes questions si l’ébranlement de droite est du même signe que celui de
gauche.

Dr. Bennaceur Said Page 49


CHAPITRE 05 : Ondes transversales sur une corde

CHAPITRE 05 : Ondes transversales sur une corde


5.1 Introduction
Considérons une corde tendue, rectiligne selon la coordonnée x, et de longueur infinie.
Nous allons étudier la propagation d’un faible ébranlement le long de la corde.
Supposons que cet ébranlement se produise suivant l’axe 0y.
Etudions l’équation du mouvement de cette corde. Nous dénoterons par T la tension à
laquelle est soumise la corde. On considère en un point d’abscisse x un segment très court
de cette corde, de longueur Δx. La masse Δm du segment est donnée par :

Δm = μ.Δx (5-1)
Où μ est la densité linéique de masse de la corde, c’est-à-dire la masse par unité de
longueur qui s’exprime en (kg/m).
5.2 Equation de propagation

Figure 5-1 : Onde sur une corde vibrante

Dans une situation hors équilibre, le segment n’est plus droit, il présente une courbure.
Nous considérons des mouvements d’oscillation de la corde de petite amplitude

(5-2)
⃗ (𝑥, 𝑡) = 𝑢(𝑥, 𝑡)𝑒𝑥
𝑢
Après l’approximation
𝜕𝑢
sin(𝜃)/𝑥 = tan(𝜃)/𝑥 = 𝜕𝑥 /𝑥
𝜕𝑢
sin(𝜃)/𝑥+∆𝑥 = tan(𝜃)/𝑥+∆𝑥 = 𝜕𝑥 /𝑥+∆𝑥

Dr. Bennaceur Said Page 50


CHAPITRE 05 : Ondes transversales sur une corde

L’allongement du segment est négligé aussi par l’approximation


La tension T= Constante
La force appliquée sur le segment dans la direction y est la résultante de la force
appliquée au point x qui force appliquée vers le bas est égale en module à
𝜕𝑢
𝐹(𝑥, 𝑡) = 𝑇𝑠𝑖𝑛(𝜃)/𝑥 ≈ 𝑇 tan(𝜃)/𝑥 = 𝜕𝑥 /𝑥 (5-3)

Pour la force au point x+𝛥𝑥 qui vers le haut


𝜕𝑢
𝐹(𝑥 + 𝛥𝑥, 𝑡) = 𝑇𝑠𝑖𝑛(𝜃)/𝑥+𝛥𝑥 ≈ 𝑇 tan(𝜃)/𝑥+𝛥𝑥 = / (5-4)
𝜕𝑥 𝑥+𝛥𝑥
La somme de ces forces dans la direction y est donc :
𝜕𝑢 𝜕𝑢 𝜕2 𝑢 (5-5)
R= [𝜕𝑥 /𝑥+∆𝑥 − 𝜕𝑥 /𝑥 ] = 𝑇 𝜕𝑡 2 ∆𝑥

Nous pouvons appliquer maintenant la deuxième loi de Newton :

𝜕2 𝑢
𝑅 = ∆𝑚 𝜕𝑡 2
(5-6)
𝜕 2𝑢 𝜕 2𝑢
𝜇 2 =𝑇 2
𝜕𝑡 𝜕𝑥

𝑇
Si on définit 𝑉 = √𝜇 son dimension d’une vitesse

𝜕 2𝑢 1 𝜕 2𝑢 (5-7)
− =0
𝜕𝑥 2 𝑉 2 𝜕𝑡 2
On appelle cette équation l’équation d’onde de la corde, Cette équation est l’équation de
d’Alembert à une dimension.
5.3 Impédance caractéristique

Si nous supposons notre onde une onde progressive harmonique nous pouvons définir par
𝑢(𝑥, 𝑡) = 𝑈0 cos(𝜔𝑡 − 𝑘𝑥) (5-8)
Avec représentation complexe
𝑢(𝑥, 𝑡) = 𝑈0 𝑒 𝑗(𝜔𝑡−𝑘𝑥) (5-9)
𝜔 2𝜋
Ou 𝑘 = = module de vecteur d’onde, 𝜆 la longueur d’onde
𝑉 𝜆

Nous définissons la force en un point la projection selon (Oy) de la force exercée en ce


point par la partie gauche de la corde sur la partie droite :
𝜕𝑢
𝐹 = −𝑇 𝜕𝑥 (5-10)

𝐹(𝑥, 𝑡) = 𝑗𝑘𝑇𝑈0 𝑒 𝑗(𝜔𝑡−𝑘𝑥) (5-11)


La vitesse de particules s’écrit :

Dr. Bennaceur Said Page 51


CHAPITRE 05 : Ondes transversales sur une corde

𝜕𝑢
𝑢̇ (𝑥, 𝑡) = = 𝑗𝜔𝑈0 𝑒 𝑗(𝜔𝑡−𝑘𝑥)
𝜕𝑥
Nous remarquons que la vitesse
𝑢̇ (𝑥, 𝑡) en phase avec la force F
On appelle impédance en un point le rapport de la l’amplitude complexe de la force à
l’amplitude complexe de la vitesse de particule
𝐹𝑦
𝑍(𝑥) =
𝑢𝑦̇
Dans le cas d’une onde progressive, on obtient :
𝑍(𝑥) = 𝜇𝑉 = √𝜇𝑇

La quantité √𝜇𝑇 définit l’impédance caractéristique de la corde

𝑍𝐶 = 𝜇𝑉 = √𝜇𝑇 (5-12)
Nous obtenons une propriété de l’onde progressive plane
𝑍(𝑥) = 𝑍𝐶 ∀𝑥
5.4 Ondes stationnaires
5.4.1 Définition
Une onde stationnaire est représentée, en notation réelle, sous la forme du produit d’une
fonction de l’espace et d’une fonction du temps :
𝑦(𝑥, 𝑡) = 𝐹(𝑥) × 𝐺(𝑡) (5-13)

5.4.2 Solution stationnaire de l’équation de l’Alembert


Soit la solution de l’équation d’Alembert :
𝑦(𝑥, 𝑡) = 𝐹(𝑥) × 𝐺(𝑡)
𝜕 2 𝑦(𝑥, 𝑡) 1 𝜕 2 𝑦(𝑥, 𝑡) " (𝑥)
1
∀𝑥, ∀𝑡, 0 = − = 𝐺(𝑡)𝐹 − 𝐹(𝑥)𝐺 " (𝑡)
𝜕𝑥 2 𝑉 2 𝜕𝑡 2 𝑐2

La solution stationnaire de l’équation de d’Alembert s’écrit :


𝜔
𝑦(𝑥, 𝑡) = (𝐶𝑠𝑖𝑛𝑘𝑥 + 𝐷𝑐𝑜𝑠𝑘𝑥)𝑐𝑜𝑠(𝜔𝑡 + 𝜑) Avec 𝑘 = (5-14)
𝑉

√𝐶 2 + 𝐷2 : Amplitude d’onde
𝜔 : Pulsation d’onde ( rad/s)
f : fréquence
2𝜋
k : nombre d’onde ( k= 𝜆 )

𝜆 : longueur d’onde (m)

Dr. Bennaceur Said Page 52


CHAPITRE 05 : Ondes transversales sur une corde

V : vitesse de phase ( m/s)


2𝜋
T : période temporelle (w= 𝑇 )

Nous pouvons considérée l’onde stationnaire comme superposition de deux ondes


progressives de mémé amplitudes de sens opposé.
Donc les constantes C et D sont définies par les conditions initiales et les conditions aux
limites.
𝑻
Avec 𝒇. 𝝀 = √ 𝝁𝟎 (5-15)
T0 : tension
µ : densité linéique

5.4.3 Réflexion et transmission


Lors du passage de l’onde mécanique incidente du milieu 1 vers le milieu 2 il existe une
onde transmise vers le milieu 2 et une onde réfléchie vers le milieu 1.

Figure 5-2: Modes de transmission et de réflexion d’une onde mécanique


Au point de discontinuité, on applique les conditions de continuité telle que :

𝑦𝑖 (0, 𝑡 ) + 𝑦𝑟 (0, 𝑡 ) = 𝑦𝑡 (0, 𝑡 )


𝜕𝑦𝑖 (𝑥,𝑡) 𝜕𝑦𝑟 (𝑥,𝑡) 𝜕𝑦𝑡 (𝑥,𝑡)
[ ] +[ ] =[ ] (5-16)
𝜕𝑥 𝑥=0 𝜕𝑥 𝑥=0 𝜕𝑥 𝑥=0

5.4.4 Types de condition aux limites pour une corde vibrante

5.4.4.1 Deux extrémités fixes

Lorsque la corde est tendue par les deux extrémités, on aura l’interférence des
ondes incidentes et les ondes réfléchies, d’où l’apparition des ondes stationnaires.
Dans ce cas-là, on applique les conditions aux bords, comme suit :

Dr. Bennaceur Said Page 53


CHAPITRE 05 : Ondes transversales sur une corde

𝑦(0, 𝑡 ) = 𝑦(𝑎, 𝑡 ) = 0 (5-17)

Figure 5.3 : Modèle de la corde tendue


On obtient les modes propres :
Appliquons ces deux conditions à la solution 𝑦(𝑥, 𝑡 )
𝑦(0, 𝑡 ) = 𝐷𝑐𝑜𝑠(𝜔𝑡 + 𝜑) = 0 → 𝐷 = 0
𝑦(𝐿, 𝑡 ) = 𝐶𝑠𝑖𝑛𝑘𝐿 𝑐𝑜𝑠(𝜔𝑡 + 𝜑) = 0 → 𝑠𝑖𝑛𝑘𝐿 = 0
𝑘𝐿 = 𝑃𝜋 → avec p=0,±1, ±2, … … ….
𝜋
𝑘𝑝 = 𝑃
𝐿
𝜋 𝑇
𝜔𝑝 = 𝑣𝑘𝑝 = 𝑝 √ (5-18)
𝐿 𝜇
𝜔𝑝 est la pulsation avec laquelle oscillent toutes les parties de la corde dans un mode P
donné (𝑘𝑝 )

Pour P=0 la corde à l’équilibre

𝜋 𝜋
Mode 1 : 𝑘1 = 𝜔1 = 𝑣 𝐿
𝐿

2𝜋 𝜋
Mode 2 : 𝑘2 = 𝜔2 = 2𝑣 𝐿 = 2 𝜔1
𝐿

3𝜋 𝜋
Mode 3 : 𝑘3 = 𝜔3 = 3𝑣 𝐿 = 2 𝜔1
𝐿

Ainsi de suite ………..

𝑦𝑝 (𝑥, 𝑡) = 𝐶𝑝 sin 𝑘𝑝 𝑥𝑐𝑜𝑠(𝜔𝑝 𝑡 + 𝜑𝑝 )

Dr. Bennaceur Said Page 54


CHAPITRE 05 : Ondes transversales sur une corde

𝑦(𝑥, 𝑡) = ∑∞
𝑝=1 𝑦𝑝 (𝑥, 𝑡) (5-19)

5.4.4.2 Une extrémité fixe et l’autre libre

Supposons qu’à x=0, il existe un excitateur vertical qui applique une force sur
l’extrémité de la corde.

y
Libre

L
x=0 Fixe
A un instant donné, la force verticale applique par la corde à travers son extrémité
sur l’excitateur est : y T
ϴ
𝐹( 𝑐𝑜𝑟𝑑𝑒 𝑠𝑢𝑟 𝑒𝑥𝑐𝑖𝑡𝑎𝑡𝑒𝑢𝑟) = 𝑇𝑠𝑖𝑛𝜃
𝜕𝑦
𝐹( 𝑐𝑜𝑟𝑑𝑒 𝑠𝑢𝑟 𝑒𝑥𝑐𝑖𝑡𝑎𝑡𝑒𝑢𝑟) = 𝑇𝑐𝑜𝑠𝜃𝑡𝑎𝑛𝜃 ≈ 𝑇 𝜕𝑥 /𝑥=0
L
x=0 Fixe

Maintenant, s’il n’existe aucun excitateur en x=0, l’extrémité gauche sur libre, la
force de réaction 𝐹( 𝑐𝑜𝑟𝑑𝑒 𝑠𝑢𝑟 𝑒𝑥𝑐𝑖𝑡𝑎𝑡𝑒𝑢𝑟) (composante verticale) doit être nulle :
𝜕𝑦
𝑇 / = 0
𝜕𝑥 𝑥=0

𝜕𝑦(0,𝑡)
=0
𝜕𝑥

Ceci est la condition pour une extrémité libre, on applique cette condition à
équation d’onde :
𝑦(𝑥, 𝑡) = (𝐶𝑠𝑖𝑛𝑘𝑥 + 𝐷𝑐𝑜𝑠𝑘𝑥)𝑐𝑜𝑠(𝜔𝑡 + 𝜑) (5-20)

Libre
𝜕𝑦
𝑦(0, 𝑡) = 0, / = 0
𝜕𝑥 𝑥=𝐿

Fixe
𝑦(0, 𝑡) = 0 D=0 L

Dr. Bennaceur Said Page 55


CHAPITRE 05 : Ondes transversales sur une corde
𝜕𝑦
/ = 0 C k cosk L = 0
𝜕𝑥 𝑥=𝐿

𝜋
cosk L = 0 k L= p p=1,3,5,7…………
2

𝑦𝑝 (𝑥, 𝑡) = 𝐶𝑠𝑖𝑛𝑘𝑝 𝑥𝑐𝑜𝑠(𝜔𝑝 𝑡 + 𝜑𝑝 )

𝜋 𝜋𝑃
kp = p , 𝜔𝑝 = 𝑣𝑘𝑝 = 𝑣
2𝐿 2𝐿

Mode 1: p=1
𝜋 𝜋
k1 = 𝜔1 = 𝑣
2𝐿 2𝐿

λ1= 4L 0
2L
𝜋
𝑦1 (𝑥, 𝑡) = 𝐶1 sin( 𝑥)𝑐𝑜𝑠(𝜔1 𝑡 + 𝜑1 )
2𝐿

Mode 2: p=2

3𝜋 𝜋
K2 = 𝜔2 = 3𝑣
2𝐿 2𝐿

4𝐿
λ2 = 0
3 L

3𝜋
𝑦2 (𝑥, 𝑡) = 𝐶2 sin( 𝑥)𝑐𝑜𝑠(𝜔2 𝑡 + 𝜑2 )
4𝐿

𝜆
2
Mode 3: p=3

5𝜋 𝜋
K3 = 𝜔3 = 5𝑣
2𝐿 2𝐿

4𝐿
λ3 = 0
5 L

5𝜋
𝑦3 (𝑥, 𝑡) = 𝐶3 sin( 𝑥)𝑐𝑜𝑠(𝜔3 𝑡 + 𝜑3 )
2𝐿

𝑃𝜋
𝑦(𝑥, 𝑡) = ∑ 𝐶𝑝 sin ( 𝑥) 𝑐𝑜𝑠(𝜔𝑝 𝑡 + 𝜑𝑝 )
2𝑙 (5-21)
𝑝=1,3,5..

Dr. Bennaceur Said Page 56


CHAPITRE 05 : Ondes transversales sur une corde

5.4.4.3 Deux extrémités libres de la corde

Nous appliquons simplement la condition


𝜕𝑦
=0 Pour x=0 et x=L
𝜕𝑥

𝑦(𝑥, 𝑡) = (𝐶𝑠𝑖𝑛𝑘𝑥 + 𝐷𝑐𝑜𝑠𝑘𝑥)𝑐𝑜𝑠(𝜔𝑡 + 𝜑)


𝜕𝑦
/ = 𝐶𝑐𝑜𝑠(𝜔𝑡 + 𝜑) = 0 C=0
𝜕𝑥 𝑥=0

𝜕𝑦
/ = 𝐷𝑘𝑠𝑖𝑛𝑘𝐿 = 0 kL=p 𝜋 𝑝 = 0, ±1, ±2, ±3 … … ..
𝜕𝑥 𝑥=𝐿

𝑝𝜋
𝑘𝑝 = p=1,2,3,……..
𝐿

Mode 1: p=1
𝜋 𝜋
k1 = 𝜔1 = 𝑣
𝐿 𝐿

λ1= 2L
0
𝜋 L
𝑦1 (𝑥, 𝑡) = 𝐷1 cos( 𝑥)𝑐𝑜𝑠(𝜔1 𝑡 + 𝜑1 )
𝐿
𝜆
Mode 2: p=2 2

2𝜋 𝜋
k2 = 𝜔2 = 2𝑣 𝐿
𝐿

λ2 = 𝐿 0
L
2𝜋
𝑦2 (𝑥, 𝑡) = 𝐷2 cos( 𝑥)𝑐𝑜𝑠(𝜔2 𝑡 + 𝜑2 )
𝐿
𝜆
Mode 3: p=3 2

3𝜋 𝜋
k3 = 𝜔3 = 3𝑣
𝐿 2𝐿

2𝐿
λ3 = 0
3 L

3𝜋
𝑦3 (𝑥, 𝑡) = 𝐷3 cos( 𝑥)𝑐𝑜𝑠(𝜔3 𝑡 + 𝜑3 )
𝐿

Dr. Bennaceur Said Page 57


CHAPITRE 05 : Ondes transversales sur une corde


𝑃𝜋
𝑦(𝑥, 𝑡) = ∑ 𝐷𝑝 cos ( 𝑥) 𝑐𝑜𝑠(𝜔𝑝 𝑡 + 𝜑𝑝 )
𝐿 (5-22)
𝑝=1

5.4.4.4 Autres types de conditions aux limites


Il existe d’autres types de conditions aux limites qu’on peut rencontrer dans les
dispositifs :
La corde est attachée à une masse qui roule verticalement dans une fente

y
F
L’équation en x=0 s’écrit :
𝜕𝑦(0, 𝑡)
𝑚𝑦̈ (0, 𝑡 ) = 𝑇
𝜕𝑥
x

La corde est attachée à un ressort de raideur k

Dans ce cas y
𝜕𝑦(0, 𝑡) F
𝑘𝑦(0, 𝑡 ) = 𝑇
𝜕𝑥

x=0
La corde est attachée à un amortisseur de coefficient β

Dans ce cas y
F
𝜕𝑦(0, 𝑡)
𝛽𝑦̇ (0, 𝑡 ) = 𝑇
𝜕𝑥

β
5.5 Partie expérimentale
Le premier dispositif expérimental qui illustre l’expérience de la production des ondes
stationnaires est celui de « MELDE » représenté dans la figure 5-3

Dr. Bennaceur Said Page 58


CHAPITRE 05 : Ondes transversales sur une corde

Figure 5-4 : Dispositif expérimental de « MELDE »

5.5.1 Expérience générale


Une corde de longueur L est placée horizontalement sur deux tiges en fer, à l’une des
extrémités de la corde nous attachons un poids de masse M. Un vibreur fait vibrer la
corde selon une intensité variable. Nous constatons qu’il se produit un mouvement
particulier de la corde qui présente des ventres et des nœuds en différents points.

Figure 5-4: Expérience de la corde tendue

5.5.2 Energie emmagasinée sur une corde vibrante


Il faut noter pour des ondes stationnaires. Il n’y a pas de phénomène de propagation, ni de
la matière ni de l’énergie :
Néanmoins, ces calculs nous serons très utiles lors de l’étude des phénomènes de
propagation.
5.5.3 Énergie cinétique
L’énergie cinétique d’un élément infinitésimal de la corde est donnée par :
1 𝜕𝑦 2
𝑑𝐸𝐶 = 2 𝜇𝑑𝑠 ( 𝜕𝑡 ) (5-23)
𝜕𝑦 2
𝑑𝑠 = √𝑑𝑥 2 + 𝑑𝑦 2 = 𝑑𝑥√1 + (𝜕𝑥 )

Dr. Bennaceur Said Page 59


CHAPITRE 05 : Ondes transversales sur une corde

Pour de faibles déviations de la corde de sa position d’équilibre, on peut négliger le terme


𝜕𝑦 2
(𝜕𝑥 ) devant le 1 et écrire 𝑑𝑠 ≈ 𝑑𝑥.

1 𝜕𝑦 2
𝑑𝐸𝐶 = 𝜇𝑑𝑥 ( )
2 𝜕𝑡 (5-24)
On définit la densité d’énergie cinétique
𝑑𝐸𝐶 1 𝜕𝑦 2
𝑒𝑐 = = 𝜇( ) (5-25)
𝑑𝑥 2 𝜕𝑡
5.5.4 Energie potentielle

L’énergie potentielle est égale au travail nécessaire pour étirer une partie
infinitésimale de la corde de longueur dx , pour quelle devenue d’une longueur ds.

F
𝑑𝑊 = 𝑇(𝑑𝑠 + 𝑑𝑦)

2 ds
𝜕𝑦
𝑑𝑠 = √𝑑𝑥 2 + 𝑑𝑦 2 = 𝑑𝑥√1 + ( )
𝜕𝑥

En faisant l’approximation :

2 2 F
√1 + (𝜕𝑦) ≈ 1 + 1 (𝜕𝑦)
𝜕𝑥 2 𝜕𝑥 dx
On se trouve donc avec

1 𝜕𝑦 2
𝑑𝐸𝑝 = 2 𝑇𝑑𝑥 (𝜕𝑥 ) (5-26)

On définit la densité d’énergie potentielle :


𝑑𝐸𝑝 1 𝜕𝑦 2 (5-27)
𝑒𝑝 = = 𝑇( )
𝑑𝑥 2 𝜕𝑥
Donc la densité d’énergie totale :

Pour obtenir l’énergie totale emmagasinée dans une corde vibrante de longueur L, on
calcule
1 𝜕𝑦 2 𝜕𝑦 2
𝑒(𝑥, 𝑡) = 𝑒𝑐 (𝑥, 𝑡) + 𝑒𝑝 (𝑥, 𝑡) = [𝜇 ( ) + 𝑇 ( ) ]
2 𝜕𝑡 𝜕𝑥
(5-28)
𝐿
𝐸𝑇 = ∫ 𝑒(𝑥, 𝑡)
0

Dr. Bennaceur Said Page 60


CHAPITRE 05 : Ondes transversales sur une corde

Exercices
Exercice 01 :
Une corde AB horizontale, de masse m et de longueur L=2m , est tendue avec une
tension T réglable. La corde est reliée à son extrémité A d’abscisse x=0 par un ressort de
raideur K dont l’extrémité C est reliée à un vibreur. Celui-ci communique au ressort un
mouvement sinusoïdal transversal d’amplitude s0=1cm et de pulsation ω=100π rad/s :
S(t)=s0 e jωt
L’autre extrémité B d’abscisse x=L de la corde est reliée à un amortisseur de coefficient
de frottement visqueux α =0.2 N.s.m-1 . Des guides parfaitement glissants n’autorisent
que les mouvements transverses des points A, B et C. Un point M d’abscisse x est repéré
à l’instant t par sa petite élongation y(x, t) . L’onde se propageant de A vers B est
progressive.
1- Pour quelle tension T, exprimée en fonction de α et , λ l’onde qui se propage le long de
la corde est progressive. En déduire la vitesse de phase V de l’onde et sa longueur d’onde
λ.
2- Montrer que les extrémités A et B de la corde vibrent en phase.
3- Exprimer les forces qui s’exercent au point A et écrire l’équation différentielle du
mouvement de ce point.
4- Etablir les expressions de l’amplitude y0 et de la phase φ de l’élongation, y(0,t) du
point A. En déduire l’élongation , y(x,t) d’un point quelconque de la corde. Dans quel cas
la phase φ est nulle ? Que vaut dans ce cas l’amplitude de vibration de chaque point de la
corde ?

Dr. Bennaceur Said Page 61


CHAPITRE 05 : Ondes transversales sur une corde

Exercice 02 :
Deux cordes 1 et 2, de masses linéiques µ1 et µ2 , sont soudées à la jonction O et sont
tendues horizontalement suivant l’axe Ox avec une tension T. On choisit l’abscisse x=0 à
la jonction O des deux cordes. Soit une onde incidente de la forme :
𝑦𝑖 (𝑥, 𝑡) = 𝑎𝑖 𝑒 𝑗(𝜔𝑡−𝑘1 𝑥)
Cette onde donne naissance, à la jonction O, une onde réfléchie 𝑦𝑟 (𝑥, 𝑡)dans la région
x<0 et une transmise 𝑦𝑡 (𝑥, 𝑡) dans la région des x>0.

1- Ecrire les expressions des 𝑦1 (𝑥, 𝑡)pour la région des x<0 et y2(x, t) pour la région des
x>0
2- Exprimer la tension instantanée en tout point transversale instantanée :
6 en tout point M d’abscisse x de la corde 1(région x<0).
3- Ecrire les deux équations de continuités au niveau de la jonction O pour le
déplacement et la tension transversale.
4- En déduire les coefficients de réflexion r et de transmission t en déplacement en
fonction de µ1 et µ2.
5- Application numérique : calculer r et t lorsqu’un fil d’acier (corde 1) de diamètre
d1=2 mm est relié à la jonction O à un autre fil d’acier (corde 2) de diamètre d2=1 mm.
6- Application numérique : mêmes questions, l’onde incidente venant de la droite pour
aller vers la jonction O.
Exercice 03 :
Une corde de longueur infinie, de masse linéique µ , est tendue sous la tension T. La
direction de la corde au repos est parallèle à l’axe x’Ox. Une masse m est fixée à l’origine
x=0 (figure a). On négligera le poids de la corde et celui de la masse m devant les autres
forces.On considère une onde incidente de déplacement de la forme :

𝑦𝑖 (𝑥, 𝑡) = 𝑎𝑖 𝑒 𝑗(𝜔𝑡−𝑘1 𝑥)

Dr. Bennaceur Said Page 62


CHAPITRE 05 : Ondes transversales sur une corde

Venant de la gauche (région x<0) et se propageant dans le sens des x croissants. Elle
donne naissance en x=0 à une onde réfléchie 𝑦𝑟 (𝑥, 𝑡) et à une onde transmise 𝑦𝑡 (𝑥, 𝑡).On
notera r et t les coefficients complexes de réflexion et de transmission relatifs aux
amplitudes des déplacements au point O.
1- Donner l’expression complexes de l’onde transversale, 𝑦1 (𝑥, 𝑡) pour la région des 0
x<0 et de l’onde transversale, 𝑦2 (𝑥, 𝑡)pour la région des x>0.
2- Montrer que la force transversale qui s’exerce sur la masse m s’écrit :

𝜕𝑦2 𝜕𝑦1
𝐹𝑦 (0, 𝑡) = 𝑇 [( ) −( ) ]
𝜕𝑥 𝑥=0 𝜕𝑥 𝑥=0

3- Montrer que, compte tenu de la condition de continuité du déplacement en x=0 et de


l’équation du mouvement de la masse m , les coefficients de réflexion r et de
transmission t s’écrivent :
𝑗𝑚𝜔 2𝑍𝐶
𝑟=− 𝑒𝑡 𝑡= 𝑎𝑣𝑒𝑐 𝑍𝐶 += √𝑇𝜇
𝑗𝑚𝜔 + 2𝑍𝐶 𝑗𝑚𝜔 + 2𝑍𝐶

Quelles sont les limites de ces deux coefficients lorsque𝑚 → 0𝑒𝑡 𝑚 → ∞?


4- Montrer que le système de la figure a) est équivalent au système de la figure b)
constitué d’une corde semi-infinie occupant la région des x<0 et fermée en O par une
impédance Z0 dont on calculera l’expression en fonction de Zc , m et ω . Donner le
coefficient de réflexion r en fonction de l’impédance réduite𝑍0′ = 𝑍0 ⁄𝑍𝐶 .
5- Calculer l’impédance mécanique ramenée Z(x) en tout point d’abscisse x ( x< 0) , en
fonction de Zc , 𝑍0′ , k et x. Que vaut cette impédance pour les positions particulières
𝜆 𝜆
suivantes : 𝑎) 𝑥𝑛 = 𝑛 2 𝑏) 𝑛𝑛′ = (2𝑛 + 1) 4

Ou λ est la longueur d’onde et n=0,-1,-2………

Dr. Bennaceur Said Page 63


CHAPITRE 05 : Ondes transversales sur une corde

6- Ecrire l’onde résultante de déplacement y1(x, t) sous la forme : 𝑦1 (𝑥, 𝑡) = 𝐴(𝑥) ×


𝑦𝑖 (𝑥, 𝑡)
7- Déterminer les positions et l’amplitude Amin des nœuds de déplacement ainsi que les
positions et l’amplitude Amax des ventres de déplacement. En déduire le taux d’ondes
𝐴𝑚𝑎𝑥
stationnaire défini par le rapport : 𝜌 = 𝐴𝑚𝑖𝑛

7- A quelle distance se trouve le premier maximum le plus proche de la terminaison x= 0.

Dr. Bennaceur Said Page 64


CHAPITRE 06 : Ondes longitudes dans les fluides

CHAPITRE 6 : Ondes longitudinales dans les fluides


6.1 Introduction
L’acoustique est la science des ondes sonores dont nous pouvons percevoir une partie
avec nos oreilles. La physique des ondes sonores ne se limite pas à la propagation du son
dans l’air car on peut, avec les mêmes équations, décrire les ondes sonores dans tous les
fluides, liquides ou gaz. Malgré la complexité de la description des fluides en
mouvement, l’équation d’onde acoustique est la plus simple qui soit, tant que les
variations de pression sont petites par rapport à la pression atmosphérique. La vitesse de
propagation des fluctuations de pression, la vitesse du son, peut être calculée à partir de
deux paramètres du fluide : sa masse volumique et sa compressibilité.
La matière est classiquement divisée en trois états : l’état solide, l’état liquide et l’état
gazeux. Les deux derniers états sont regroupés sous le nom de fluide. À la différence
d’un solide, un fluide ne peut conserver une forme propre : par exemple un volume d’eau
va adopter la forme du récipient et un gaz va occuper tout le volume disponible. Les
composants des fluides, atomes ou molécules, n’ont pas de positions fixes et changent
constamment de voisins.
6.2 Ondes planes dans un tuyau cylindrique
6.2.1 Equation d’ondes dans un gaz.
On définit alors, les Ondes élastiques dans les fluides comme des ondes mécaniques qui
se propagent dans les gaz ou dans les liquides.
On suppose que les fluides soient parfaits, il n’y aura pas d’absorption.
L’onde élastique dans l’air est due à la propagation de la variation de pression
Dans la suite, nous utiliserons les symboles suivants pour étudier l’onde acoustique qui se
propage suivant l’axe des x :
x : coordonnée à l’équilibre d’une particule du milieu.
𝑈𝑥 : composante suivant l’axe des x du déplacement de particule par rapport à la position
d’équilibre.
ρ0 : masse volumique du fluide à l’équilibre
P : pression instantanée en un point quelconque
P0 : pression à l’équilibre
p = P − P0 : surpression ou pression acoustique
c : vitesse de propagation de l’onde

Dr. Bennaceur Said Page 67


CHAPITRE 06 : Ondes longitudes dans les fluides

Figure 6-1 : propagation d’une onde élastique sinusoïdale plane dans l’air
Dans un gaz, la pression est souvent de l’ordre de la pression atmosphérique, P= 105 Pa
Une vibration mécanique est véhiculée par une surpression locale: p( x,t ) telle que :
|𝑝| ≪ 𝑃
Sachant que les mouvements aux fréquences sonores sont trop rapides pour qu’il y ait des
échanges thermiques, les variations de pression sont donc adiabatiques et l’évolution d’un
volume V est reliée à celle de la pression comme suit :
𝑑𝑃 𝑑𝑉 (6-1)
+𝛾 =0
𝑃 𝑉
𝐶𝑃
Ou 𝛾 = est le rapport des chaleurs spécifiques à pression constante et à volume
𝐶𝑉

constant.
On obtient la relation suivante : PVγ=cte
On entend par particule, un élément de volume contenant des millions de molécules de
telle sorte qu’il puisse être considéré comme continu, mais toutefois suffisamment petit
pour que les grandeurs acoustiques comme la pression, la masse volumique et la vitesse
de particule puissent être considérées comme constantes dans cet élément de volume.
Dans ce qui suit, nous négligerons les effets de la gravitation de telle sorte que P0 et ρ0
sont uniformes dans tout le milieu. On suppose d’autre part que le milieu est homogène,
isotrope et parfaitement élastique, c’est-à-dire non dissipatif.
6.2.2 Ondes planes dans un tuyau cylindrique
6.2.2.1 Equation d’ondes dans un fluide
Beaucoup de problèmes de bruit industriel sont liés à la propagation d’ondes sonores dans
des conduites et des tuyaux. Cette propagation est conditionnée par la longueur d’onde
acoustique λ telle que :

Dr. Bennaceur Said Page 68


CHAPITRE 06 : Ondes longitudes dans les fluides
𝑐
𝜆= ≫𝐷 (6-2)
𝑓

D : Diamètre de la conduite sonore


c : Célérité du son dans l’air qui est égale à 344m/s à 20°c
f : Fréquence de l’onde
Soit une membrane vibrante émise une onde plane progressive dans un fluide uniforme,
de masse volumique ρ et de pression P0 à l’équilibre. Cette onde se propage dans la
direction x positif.
Le phénomène de propagation des ondes sonores dans le fluide est dû principalement par
une succession de compressions et de détentes des tranches de fluide voisines.
On considère une tranche de fluide d’épaisseur dx située aux abscisses x et dx .Soient les
pressions P(x) et P(x+dx) agissant sur les plans x et dx respectivement qui génèrent le
mouvement de la tranche comme le montre la figure 6-2.

Soient U( x ) et U( x+ dx ) les déplacements à l’instant t des plans d’abscisse x et x +dx

respectivement

Figure 6-2: Propagation d’onde acoustique dans un fluide

Dr. Bennaceur Said Page 69


CHAPITRE 06 : Ondes longitudes dans les fluides

En appliquant la loi de la dynamique de Newton :


𝜕2 𝑈
𝑑𝑚 𝜕𝑡 2 = 𝐹(𝑥 + 𝑑𝑥) − 𝐹(𝑥) (6-3)

F( x ) et F( x+ dx ) sont des forces d’actions appliquée aux plans d’abscisses x et x+dx


respectivement.
La résultante des forces s’écrit :
𝐹(𝑥 + 𝑑𝑥) − 𝐹(𝑥) = −𝑆(𝑃(𝑥 + 𝑑𝑥) − 𝑃(𝑥))
Avec :
𝜕𝑝
(𝑃(𝑥 + 𝑑𝑥) − 𝑃(𝑥)) = 𝑑𝑝 = 𝑑𝑥
𝜕𝑥
D’où :
𝜕 2𝑈 𝜕𝑃
𝑑𝑚 2
= −𝑆 𝑑𝑥
𝜕𝑡 𝜕𝑥
On définit la surpression p d’un fluide compressible comme suit :
1 𝜕𝑈
𝑝=− (6-4)
𝜒 𝜕𝑥

Ou χ est appelé le coefficient de compressibilité.


En injectant la valeur de surpression à l’équation du mouvement ; on obtient :
𝜕 2𝑈 𝜕 1 𝜕𝑈
𝜌𝑆𝑑𝑥 2
= −𝑆 (− )𝑑𝑥
𝜕𝑡 𝜕𝑥 𝜒 𝜕𝑥
D’où l’équation de la propagation des ondes sonores dans le fluide s’écrit :
𝜕2 𝑈 1 𝜕2 𝑈 𝜕2 𝑈 𝜕2 𝑈
= 𝜒𝜌 𝜕𝑥 2 ⟺ = 𝑉 2 𝜕𝑥 2 (6-5)
𝜕𝑡 2 𝜕𝑡 2

1
Avec 𝑉 = √𝜌𝜒 est appelé la célérité de l’onde sonore dans le fluide.

6.2.2.2 Vitesse du son

Le phénomène de propagation étant un processus adiabatique, on définit la vitesse du son


dans un fluide :

𝛾𝑃0
𝑉=√ (6-6)
𝜌0

Exemple : Dans l’air, dans les conditions normales T = 20◦C et P0 = 105N·m−2, γ = 1.4
et ρ0 = 1.29 kg·m−3, on en déduit 𝑉 ≈ 330𝑚. 𝑠 −1 .
 La valeur de la pression à l’équilibre dépend fortement de la température. Pour
une mole de gaz parfait, on a :

Dr. Bennaceur Said Page 70


CHAPITRE 06 : Ondes longitudes dans les fluides

𝛾𝑅𝑇
𝑉=√ (6-7)
𝑀

M représente la masse molaire du gaz.


 Dans un gaz parfait, la vitesse de propagation du son est proportionnelle à la
racine carrée de la température mesurée en ◦K.
6.2.2.3 Impédance acoustique et impédance caractéristique.

Une onde acoustique sinusoïdale, s’écrit :


𝑥
𝑝(𝑥, 𝑡) = 𝑝0 𝑐𝑜𝑠 [𝜔(𝑡 − 𝑉)] (6-8)

On définit le module du vecteur d’onde k par


𝜔
𝑘=
𝑉
D’où
𝑝(𝑥, 𝑡) = 𝑝0 𝑐𝑜𝑠(𝜔𝑡 − 𝑘𝑥) (6-9)
En notation complexe, l’onde progressive sinusoïdale s’écrit

𝑝(𝑥, 𝑡) = 𝑝0 𝑒 𝑗(𝜔𝑡−𝑘𝑥) (6-10)


La relation liant la pression acoustique et la compressibilité, avoir
𝜕𝑈
𝑝(𝑥, 𝑡) = −𝜅 𝜕𝑥 (6-11)

Permet d’écriture
1
𝑈(𝑥, 𝑡) = − ∫ 𝑝(𝑥, 𝑡)𝑑𝑥
𝜅
1
𝑈(𝑥, 𝑡) = − ∫ 𝑝0 𝑒 𝑗(𝜔𝑡−𝑘𝑥) 𝑑𝑥
𝜅
𝑝
0
𝑈(𝑥, 𝑡) = 𝑗𝜅𝑘 𝑒 𝑗(𝜔𝑡−𝑘𝑥) (6-12)

𝑝0
𝑈(𝑥, 𝑡) = 𝑒 𝑗(𝜔𝑡−𝑘𝑥) (6-13)
𝑗𝜔𝜌0 𝑉
La variation de cette dernière expression par rapport au temps permit d’écrire la vitesse
𝜕𝑈(𝑥,𝑡) 𝑝
𝑈̇(𝑥, 𝑡) = 𝜕𝑥 = 𝜌 0𝑉 𝑒 𝑗(𝜔𝑡−𝑘𝑥) (6-14)
0

On constate que pour une onde progressive la vitesse de particule est en phase avec la
pression acoustique.
On appelle impédance acoustique en un point le rapport de l’amplitude complexe de la
pression à l’amplitude complexe de la vitesse de particule :

Dr. Bennaceur Said Page 71


CHAPITRE 06 : Ondes longitudes dans les fluides
𝑝
𝑍(𝑥) = 𝑈̇ (6-15)
Dans le cas d’une onde progressive, on obtient :
Z (x) = ρ0V (6-16)
Le produit ρ0V définit l’impédance acoustique caractéristique du fluide

Zc = ρ0V (6-17)

On obtient une propriété de l’onde plane progressive :


Z (x) = Zc ∀ x (6-18)

6.2.2.4 Energie transportée par une onde.

 Energie cinétique :

Le déplacement d’une position du fluide d’une position du fluide comprise entre les deux
plans x et x+dx est donné par U(x)

A
x

Δx
x x+Δx
A : section du plan

L’énergie cinétique de cette position lors de son déplacement est donc égale à :
1 𝜕𝑈(𝑥, 𝑡) 2 (6-19)
𝑑𝐸𝐶 = (𝜌0 . 𝐴. Δ𝑥) ( )
2 𝜕𝑡
Dans le cas ou Δ𝑥 est très petit, on peut définir l’énergie cinétique volumique ou la
densité volumique d’énergie cinétique comme

𝑑𝐸𝐶 1 𝜕𝑈(𝑥, 𝑡) 2 (6-20)


𝑒𝑐 = = 𝜌0 ( )
Ou 𝑑𝑉 = 𝐴. Δ𝑥 𝑑𝑉 2 𝜕𝑡

 Energie potentielle :
Le déplacement du fluide s’accompagne d’une petite variation du volume de la position
considérée ci-dessus sous l’effet de la surpression. C’est comme si on vient de
comprimer un ressort en fournissant un travail, qui va être emmagasinée dans cette
position puis dégage une fois que l’effet de compression ait cessé.
L’énergie potentielle est donc égale au travail nécessaire pour changer le volume de la
position, donc :
𝑑𝑊 = −𝑃𝑑𝑉 > 0 (6-21)

Dr. Bennaceur Said Page 72


CHAPITRE 06 : Ondes longitudes dans les fluides

Le signe moins est présent du fait que dV est négatif (compressibilité)


𝜕𝑉
Mais on a 𝑑𝑉 = 𝜕𝑃 𝑑𝑃
𝜕𝑉
Ce qui donne 𝑑𝑊 = −𝑃 𝜕𝑃 𝑑𝑃
𝜕𝑃 1
Rappelons que le « Bulk modulus » 𝐵 = −𝑉 𝜕𝑉 ou 𝐵=𝜒

(𝜒)𝑐𝑜𝑝𝑟𝑒𝑠𝑠𝑖𝑏𝑖𝑙𝑖𝑡é
Puisque le volume V de la position du fluide est égal à A.dx, on réécrit l’expression du
travail sous forme :
𝑝
𝑑𝑊 2 = 𝐴. 𝑑𝑥 𝐵 𝑑𝑃 (6-22)

Le travail accumulé par le fluide lorsque la surpression passe de 0 à ΔP à partir de la


pression d’équilibre peut s’écrire sous la forme :
1 𝑃0 +Δ𝑃
𝑑𝑊 = 𝐴. 𝑑𝑥 ∫ 𝑃 𝑑𝑃
𝐵 𝑃0

𝑃0 +Δ𝑃
1 𝑃2
𝑑𝑊 = 𝐴. 𝑑𝑥 [ ]
𝐵 2 𝑃0
11
𝑑𝑊 = 𝐴. 𝑑𝑥 [(𝑃 + Δ𝑃)2 − 𝑃0 2 ]
𝐵2 0
1
𝑑𝑊 = 𝐴. 𝑑𝑥 [(Δ𝑃)2 + 2𝑃0 Δ𝑃]
2𝐵
On sait que :
𝜕𝑈
Δ𝑃 = −𝐵 𝜕𝑥 (6-23)

Ce qui, en remplaçant Δ𝑃 dans l’équation ci-dessous, donne


1 𝜕𝑈
𝑑𝑊 = 𝐴. 𝑑𝑥 [(𝐵 𝜕𝑥 )2 + 2𝑃0 Δ𝑃] (6-24)
2𝐵

Remarquons que Δ𝑃, la surpression, peut prendre des valeurs positives comme sa valeur
moyenne (dans le temps et dans l’espace est nulle, et donc le terme 𝑃0 Δ𝑃 ne nous
renseigne pas sur le transfert d’énergie.
Dans le cas d’une onde sinusoïdale :
Δ𝑃 = Δ𝑃𝑚𝑎𝑥 cos(𝜔𝑡 − 𝑘𝑥)
1 𝑇
〈Δ𝑃〉 = ∫ Δ𝑃 cos(𝜔𝑡 − 𝑘𝑥) 𝑑𝑡 = 0
𝑇 0 𝑚𝑎𝑥
Ou encore
1 𝜆
〈Δ𝑃〉 𝑇 = ∫ Δ𝑃 cos(𝜔𝑡 − 𝑘𝑥) 𝑑𝑥 = 0
𝜆 0 𝑚𝑎𝑥

Dr. Bennaceur Said Page 73


CHAPITRE 06 : Ondes longitudes dans les fluides

L’expression du travail devient donc


𝐵𝑆 𝜕𝑈 2
𝑑𝑊 = ( 𝜕𝑥 ) 𝐴𝑑𝑥 (6-25)
2

Ce qui est l’énergie potentielle emmagasinée dans la position de volume dV= A.dx
𝐵𝑆 𝜕𝑈 2
𝑑𝐸𝑃 = ( ) 𝐴𝑑𝑥
2 𝜕𝑥
La densité volumique de l’énergie potentielle de la position est :

𝑑𝐸𝑃 𝐵𝑆 𝜕𝑈 2 (6-26)
𝑒𝑝 = = ( )
𝑑𝑉 2 𝜕𝑥

La densité volumique de l’énergie forme :

1 𝜕𝑈(𝑥, 𝑡) 2 𝐵𝑆 𝜕𝑈 2 (6-27)
𝑒𝑡 = 𝑒𝑐 +𝑒𝑝 = 𝜌0 ( ) + ( )
2 𝜕𝑡 2 𝜕𝑥
L’énergie totale sonore d’un tube de longueur L obtenue, en faisant l’intégrale :
𝐿
𝐸𝑡 = ∫0 𝑒𝑡 . 𝐴. 𝑑𝑥 (6-28)

1 𝐿
𝜕𝑈(𝑥, 𝑡) 2 𝐿
𝜕𝑈 2
𝐸𝑡 = [𝐴. 𝜌0 ∫ ( ) 𝑑𝑥 + 𝐴𝐵𝑆 ∫ ( ) 𝑑𝑥 ] (6-29)
2 0 𝜕𝑡 0 𝜕𝑥

Dans le cas d’une onde progressive sinusoïdale monochromatique :

𝑈(𝑥, 𝑡) = 𝑈𝑚𝑎𝑥 cos(𝜔𝑡 ± 𝑘𝑥)


𝜕𝑈
On aura : = −𝜔𝑈𝑚𝑎𝑥 sin(𝜔𝑡 ± 𝑘𝑥)
𝜕𝑡

𝜕𝑈
Et = ±𝑘𝑈𝑚𝑎𝑥 sin(𝜔𝑡 ± 𝑘𝑥)
𝜕𝑥

𝐴 𝐿 𝐿
𝐸𝑡 = 2 [𝜌0 ∫0 𝜔2 U𝑚𝑎𝑥
2
sin2 (𝜔𝑡 ± 𝑘𝑥)𝑑𝑥 + 𝐵𝑆 ∫0 𝑘 2 U𝑚𝑎𝑥
2
sin2 (𝜔𝑡 ± 𝑘𝑥)𝑑𝑥] (6-30)

2𝜋
Si L est un multiple de la longueur d’onde : L= n λ = 𝑛 𝑘

𝐿 𝐿
∫0 sin2 (𝜔𝑡 ± 𝑘𝑥)𝑑𝑥 Serait égale à 2
d’où

𝐴 𝜔2 U𝑚𝑎𝑥
2
𝐵𝑆 𝑘 2 U𝑚𝑎𝑥
2
𝐸𝑡 = [𝜌0 𝐿 +𝐿 ]
2 2 2
(6-31)
𝐴𝐿 2
𝐸𝑡 = U [𝜌 𝜔2 + 𝐵𝑆 𝑘 2 ]
4 𝑚𝑎𝑥 0

Dr. Bennaceur Said Page 74


CHAPITRE 06 : Ondes longitudes dans les fluides

La relation de dispersion d’écrit comme :


𝜔 = 𝜗𝑠𝑜𝑛 𝑘 (6-32)
Ce qui permet d’écrire :
𝑉
2 𝐵 (6-33)
Où V=A.L 𝐸𝑡 = 4 U𝑚𝑎𝑥 [𝜌0 + 𝜗2 𝑆 ] 𝜔2
𝑠𝑜𝑛

𝐵𝑆
On sait aussi que : 𝜗 2 𝑠𝑜𝑛 = 𝜌0
𝐸𝑡 1 2
D’où = 4 𝜌0 U𝑚𝑎𝑥 𝜔2 (6-34)
𝑉

6.2.2.5 Intensité de l’onde sonore


L’intensité acoustique est la puissance transportée par les ondes sonores, par unité de
surface, mesurée perpendiculairement à la direction de ce transfert. Pour calculer cette
’intensité il faut calculer l’énergie lorsque l’onde traverse une surface S perpendiculaire à
la direction de propagation pendant un intervalle de temps.

Figure 6-3 : Intensité de l’onde acoustique

Cette énergie dE est égale à l’énergie contenue dans un volume S V dt et elle égale à

𝑑𝐸 = 𝜀 𝑆 𝑉𝑑𝑡 (6-35)
𝑑𝐸
D’où la puissance 𝑝= =𝜀𝑆𝑉 (6-36)
𝑑𝑡

On déduit l’expression de l’intensité de l’onde acoustique


1
𝐼(𝑡) = 𝑆 𝑝 (6-37)

𝐼(𝑡) = 𝜀 𝑉 (6-38)
𝑝02
𝐼(𝑡) = 𝑐𝑜𝑠 2 (𝜔𝑡 − 𝑘𝑥) (6-39)
𝜌0 𝑉

On appelle intensité de l’onde acoustique la valeur moyenne

Dr. Bennaceur Said Page 75


CHAPITRE 06 : Ondes longitudes dans les fluides
1 𝑡+𝑇
𝐼 = 𝑇 ∫𝑡 𝐼(𝑡)𝑑𝑡 (6-40)

𝑝2 𝑝2
𝐼 = 2𝜌0𝑉 = 2𝑍0 (6-41)
0 𝐶

6.2.2.6 Niveau sonore


Le domaine de fréquences accessibles à une oreille humaine s’étend environ de 20 Hz
à 20kHz. La gamme des intensités sonores accessible est très large. Le seuil
d’audition correspond à une intensité sonore (pour une oreille moyenne à environ
1500 Hz ) 𝐼0 = 10−12 𝑊𝑎𝑡𝑡/𝑚2 . Le seuil de douleur se situe à 1𝑊𝑎𝑡𝑡/𝑚2 . Il est
donc intéressant d’utiliser une échelle logarithmique pour repérer les intensités
sonores.
Le niveau sonore est défini en décibels (dB) par
𝐼
𝐿 = 10 𝑙𝑜𝑔 𝐼 (6-42)
0

𝐼0 = 10−12 𝑊𝑎𝑡𝑡/𝑚2 .

Exemples : niveau d’un réveil électronique est 20 dB , une pièce ( salle) silencieuse
=30dB

Une rue animée = 75 dB .

En titre d’expérience calculons l’amplitude de déplacement d’une onde sonore d’un


niveau sonore très élevé :L=120dB , de fréquence f=100Hz.
On sait que
𝐼
𝐿 = 10 log 𝐼 𝐼 = 1𝑊𝑎𝑡𝑡/𝑚2
0

1 2 2
D’autre part on a 𝐼0 = 2 𝜌0 𝑈𝑚 𝜔 𝜐𝑠𝑜𝑛 𝜌0 = 1.29𝑘𝑔/𝑚3 𝜐𝑠𝑜𝑛 = 337𝑚/𝑠

2𝐼 2×1
𝑈 = √𝜌 2
= √1.29×(2×𝜋×1000)2 ×337 = 10.79𝓃𝑚
0 𝜔 𝜐𝑠𝑜𝑛

6.2.2.7 Coefficients de réflexion et de transmission d’ondes (conditions aux limites)


Soit deux milieux fluides semi-infinis séparés par une surface plane. Choisissons un
repère orthonormé de telle sorte que le plan yOz coïncide avec la surface de séparation.
Lorsque une onde acoustique provenant de −∞, se propageant dans le premier dans la
direction de l’axe des x arrive à la surface de séparation, elle donne naissance à deux
ondes.
 une onde réfléchie qui se propage dans le premier milieu dans le sens des x
décroissant.

Dr. Bennaceur Said Page 76


CHAPITRE 06 : Ondes longitudes dans les fluides

 une onde transmise qui se propage dans le second milieu dans le sens des x
croissant.

Figure 6-4 : Réflexion à une interface fluide-fluide

L’onde résultante dans le premier milieu (x ≤ 0) est caractérisée par :

𝑝1 (𝑥, 𝑡) = 𝑝𝑖 (𝑥, 𝑡) + 𝑝𝑅 (𝑥, 𝑡) (6-43)


𝑗(𝜔𝑡−𝑘1 𝑥) 𝑗(𝜔𝑡−𝑘1 𝑥)
𝑝1 (𝑥, 𝑡) = 𝑝𝑖 𝑒 + 𝑝𝑅 𝑒
1
𝑈̇1 (𝑥, 𝑡) = 𝑍 [𝑝𝑖 𝑒 𝑗(𝜔𝑡−𝑘1 𝑥) − 𝑝𝑅 𝑒 𝑗(𝜔𝑡−𝑘1 𝑥) ] (6-44)
1

Dans le deuxième milieu :


𝑝2 (𝑥, 𝑡) = 𝑝𝑇 (𝑥, 𝑡)
𝑝2 (𝑥, 𝑡) = 𝑝𝑇 𝑒 𝑗(𝜔𝑡−𝑘1 𝑥)
1
𝑈̇2 (𝑥, 𝑡) = 𝑍 𝑝𝑇 𝑒 𝑗(𝜔𝑡−𝑘𝑥) (6-45)
2

On déduit
𝑝𝑖 + 𝑝𝑅 = 𝑃𝑇
{ 1 (𝑝 − 𝑝 ) = 1 𝑃
𝑧 𝑖 𝑅 𝑍 𝑇
1 2 (6-46)
𝑝𝑅 𝑝𝑇
1+ 𝑝 = 𝑝
𝑖 𝑖
Ou encore { 𝑝𝑅 𝑍1 𝑝𝑇
1− = (6-47)
𝑝𝑖 𝑍2 𝑝𝑖

Nous déduisons
 le coefficient de réflexion pour la pression.
𝑝𝑅
𝑅𝑃 = (6-48)
𝑝𝑖

 le coefficient de transmission pour la pression.

𝑝𝑇
𝑇𝑃 = (6-49)
𝑝𝑖
Les deux relations de continuité s’écrivent comme suit :

Dr. Bennaceur Said Page 77


CHAPITRE 06 : Ondes longitudes dans les fluides

1 + 𝑅𝑃 = 𝑇𝑃
𝑍
1 − 𝑅𝑃 = 𝑍1 𝑇𝑃
2 (6-50)

On en déduit les coefficients de réflexion et de transmission


𝑍 −𝑍
𝑅𝑃 = 𝑍1 +𝑍2
1 2
(6-51)
2𝑍1
𝑇𝑃 = 𝑍
1 +𝑍2

En tenant compte des relations pi = Z1 𝑈̇𝑖 , pR = −Z1 𝑈̇𝑅 et pT = Z2 𝑈̇𝑇 , on peut calculer
les coefficients de réflexion et de transmission pour la vitesse de particules et pour le
déplacement de particules :
𝑍 −𝑍
𝑅𝑈̇ = 𝑅𝑈 = 𝑍1 +𝑍2
1 2

2𝑍1
𝑇𝑈̇ = 𝑇𝑈 = 𝑍
1 +𝑍2

En tenant compte des relations𝐼𝑖 = 𝑝𝑖 /2𝑍1 , 𝐼𝑅 = 𝑝𝑅2 /2𝑍1 et , 𝐼𝑇 = 𝑝𝑇2 /2𝑍2 , on peut
calculer les coefficients de transmission et réflexion pour intensité acoustique

𝐼𝑅 𝑍1 −𝑍2 2
𝛼𝑅 = =[ ]
𝐼𝑖 𝑍1 +𝑍2
(6-52)
𝐼𝑇 4𝑍1 𝑍2
𝛼𝑇 = =
𝐼𝑖 [𝑍1 + 𝑍2 ]2

6.3 Effet Doppler


6.3.1 Historique
L’effet Doppler a été découvert en 1842 par Johann Christian Doppler (1803-1853
/Autrichien). Puis il a aussi été redécouvert 6 ans plus tard par Hippolyte Fizeau (1819-
1896/Français), d’où cet effet porte aussi le nom d’effet Doppler-Fizeau.

Figure 6-5 : Photos Doppler et Fizeau

Dr. Bennaceur Said Page 78


CHAPITRE 06 : Ondes longitudes dans les fluides

6.3.2 Définition
La longueur d’onde des ondes émises par une source en mouvement dépend de l’endroit
d’où elles sont observées. L’effet Doppler est présent dans notre quotidien : radar, sirène
d’un camion en mouvement…Cet effet se manifeste aussi bien par les ondes sonores que
les ondes lumineuses. On suppose une source S émettant une onde dans l’air, de longueur
d’onde λs de vitesse vs. Simultanément, la source se déplace avec une vitesse v

Figure 6-6 : Le schéma simplifie de l’effet doppler

Lorsque la source se rapproche de l’observateur B, les ondes émises par la source ont une
longueur d’onde λ inférieure à λs. Réciproquement, l’observateur A reçoit des ondes dont
la longueur λ est supérieure à λs.
𝑉𝑆 𝑉𝑆
On a ainsi, entre la fréquence de l’émission 𝑓𝑠 = et la fréquence𝑓 = de l’onde
𝜆𝑆 𝜆
reçue par les observateurs, les relations

 Pour l’observateur A : on a alors 𝑓 < 𝑓𝑠


𝑉𝑆 (6-53)
𝑓 = 𝑓𝑠 × (1 − )
𝑉
 Pour l’observateur B : on a alors 𝑓 > 𝑓𝑠
𝑉𝑆 (6-54)
𝑓 = 𝑓𝑠 × (1 + )
𝑉
 Pour que ces relations soient cohérentes, l’air sera supposé immobile par rapport
au sol.
 La vitesse de la source (v) est très inférieure à la vitesse de l’onde (vs).

Dr. Bennaceur Said Page 79


CHAPITRE 06 : Ondes longitudes dans les fluides

 Pour l’observateur B, le son perçu sera plus aigu, au contraire de l’observateur A


qui aura une perception du son plus grave.
Exemple :
Soit c = 340 m/s ; vitesse du véhicule Ve = 25 m/s (90 km/h) ; fréquence de la sirène
400Hz.Un observateur immobile voyant approcher un camion de pompier entendra la
sirène a une fréquence
340
𝑓= × 400 = 431.75 ≈ 432𝐻𝑧
(340 − 25)
Et lorsque le camion s'éloigne ?
Le schéma ci-contre simplifie permet de répondre a la
question : Lorsque l'émetteur aura dépassé le récepteur,
les fronts d'onde ont de plus en plus de distance a
parcourir pour atteindre le récepteur. L'intervalle de
temps entre deux fronts d'onde à la réception est
maintenant supérieur à l'intervalle réel lors de
l'émission. Chaque front d'onde est espace du précèdent
de :
𝑑 = 𝑑𝑓 + 𝑑𝑒 = 𝑐. 𝑇𝑒 + 𝑉𝑒 . 𝑇𝑒 = (𝑐 + 𝑉𝑒 ). 𝑇𝑒
Le son perçu par le récepteur a donc une fréquence apparente d’expression :
𝑐
𝑓= .𝑓
𝑐 + 𝑉𝑒 𝑒
𝑐
< 1 donc 𝑓 < 𝑓𝑒 . L’observateur perçoit un son plus grave.
𝑐+𝑉𝑒

Reprenons notre exemple :


c = 340 m/s ; Ve = 25 m/s (90 km/h) ; fe =400 Hz.
L'observateur immobile voyant s'éloigner le camion, entendra la sirène a une fréquence
340
𝑓= × 400 = 472.6 ≈ 373𝐻𝑧
(340 + 25)
A noter que l’écart n’est pas « symétrique » : gain de 32Hz à l’approche, perte de
27 Hz à l’éloignement.

Dr. Bennaceur Said Page 80


CHAPITRE 06 : Ondes longitudes dans les fluides

Exercices
Exercice 01 :
Soit I l’intensité acoustique d’une onde sonore et l’intensité de référence prise égale à
10-12 W/m2 On appelle IdB le niveau d’intensité acoustique exprimé en décibel (dB) :
1
𝐼𝑑𝐵 = 10 log 10

1° Calculer en décibels le seuil de douleur pour l’oreille humaine sachant que son
intensité est 10 W/m2
2° Un des haut-parleurs d’une chaîne haute-fidélité (Hi-Fi) fonctionne et le niveau
d’intensité acoustique où se trouve l’auditeur vaut 65 dB. Quelle est l’intensité du son
perçu ?
3° Les deux haut-parleurs fonctionnent. Le seuil de douleur est-il atteint ? Quel est le
niveau acoustique pour l’auditeur.
Exercice 02 :

1° Chez l’homme, l’amplitude de pression maximale tolérable par l’oreille (seuil de


douleur) est 𝑝𝑚 = 28𝑃𝑎 . Quelle est l’amplitude de déplacement 𝑢𝑚 d’une telle onde
sonore dans l’air à la fréquence = 1000 𝐻𝑧 . Calculer son intensité acoustique.
2° Mêmes questions si l’amplitude de pression minimale tolérable par l’oreille humaine
(seuil d’audibilité) est𝑝𝑚 = 2.8 . 10−5 𝑃𝑎 à la même fréquence.
𝐾𝑔
On donne pour l’air : 𝜌 = 1.2 𝑚3 𝑒𝑡 𝑉 = 340 𝑚/𝑠 .

Exercice 03 :
On définit l’impédance acoustique spécifique d’un milieu par le rapport, en notation
complexe, de la pression acoustique sur la vitesse des particules :
𝑃
𝑍= (𝑒𝑥𝑝𝑟𝑖𝑚é𝑒 𝑒𝑛 𝑅𝑎𝑦𝑙𝑠)
𝑢̇
1° Calculer Z pour une onde progressive plane sinusoïdale se propageant vers les x
𝐾𝑔
croissants. Application numérique : on donne pour l’air 𝜌 = 1.21 𝑚3 𝑒𝑡 𝑉 = 343 𝑚/𝑠

2° Une onde incidente , se propageant dans 𝑢(𝑥, 𝑡) = 𝑢0 𝑒 𝑗(𝜔𝑡−𝑘𝑥) , se propageant dans


l’air, arrive en incidence normale sur la surface d’un solide dont l’impédance acoustique
𝑃(0,𝑡)
spécifique en x = 0 vaut 𝑍0 = .
𝑢̇ (0,𝑡)
𝑃(𝑥,𝑡)
a) Calculer, en un point d’abscisse x(𝑥 ≤ 0) , l’impédance ramenée 𝑍(𝑥) =
𝑢̇ (𝑥,𝑡)

en fonction de 𝜌, 𝑉, 𝑍0 𝑒𝑡 tan(𝑘𝑥) .

Dr. Bennaceur Said Page 81


CHAPITRE 06 : Ondes longitudes dans les fluides

b) Quelle est la valeur de l’impédance 𝑍(−𝜆⁄4) ramenée au point 𝑥 = −𝜆⁄4


Etudier les cas particuliers 𝑍0 → ∞ ; 𝑍0 = 𝜌𝑉 ; 𝑍0 = 𝑗𝑋 (𝑜𝑢 𝑋 𝑒𝑠𝑡 𝑢𝑛 𝑟é𝑒𝑙 𝑝𝑜𝑠𝑖𝑡𝑖𝑓) et
donner pour chacune de ces impédances 𝑝(−𝜆⁄4, 𝑡) et 𝑢̇ (−𝜆⁄4 , 𝑡)
c) Mêmes questions pour 𝑥 = −𝜆⁄2
Exercice 04 : (effet Doppler) principe du Radar
Considérons une source sonore, immobile, de fréquence f. Considérons maintenant un
objet s'approchant de la source sonore à la vitesse v. Au même endroit que la source
sonore, un récepteur sonore capte le son de la source sonore réfléchi par l’objet en
mouvement. On constate que la fréquence f ' du son reçu a une fréquence différente du
son émis f.
1. De quel phénomène s'agit-il ?
2. On remarque également que f ' est supérieure à f. Pourquoi ?
Nous allons maintenant établir la relation mathématique liant f ' et f :
On note T la période de la source sonore et T’ la période du son reçu.
Considérons le son émis par la source à l’instant t = 0.
Le son se propage et atteint l’objet : on note d la distance qui sépare la source sonore et
l'objet à ce moment.
3. Quelle durée (τ) met le son émis à l'instant t = 0 pour aller jusqu'à l'objet et revenir ?
4. Quelle durée (τ’) met le son émis à l'instant t = T pour aller jusqu'à l'objet et revenir ?
5. Quelle est la relation entre τ, τ’, T et T’ ?
6. En déduire que :
v : vitesse de l’objet
1 𝑐
c : vitesse du son dans l’air 𝑓′ = 𝑓 2𝑉 = 𝑓 𝑐−2𝑉
1−
𝑐

Dr. Bennaceur Said Page 82


CHAPITRE 07 : Ondes élastiques dans les solides

CHAPITRE 07 : Ondes élastiques dans les solides


7.1 Introduction
On définit une propagation d’onde dans un milieu matériel comme étant une perturbation
évolutive du milieu sous l’action d’une excitation. Puisque le milieu est constitué de
plusieurs particules en interaction entre elles, les forces internes sont responsables du
déplacement de la perturbation. Cette propagation dépend des propriétés physiques du
milieu où l’onde se propage.
7.2 Définition et propriétés élastiques d’un solide
Les solides sont des milieux matériels organisés et ordonnés .Les mouvements de chaque
élément d’un solide sont fortement corrélés à ces voisins.
Les milieux solides ont toujours une élasticité qui rend possible la propagation des
déformations produites par une force extérieure. Les caractéristiques des ondes observées
dépendent du milieu ; ces ondes peuvent être longitudinales ou transverses, la célérité
dépend de la densité et des coefficients d’élasticité.
Considérons le cas de la compression et des ondes longitudinales. Les déformations dans
les milieux solides dues aux vibrations et aux ondes qui en résultent peuvent être
analysées en partant de leurs propriétés d’élasticité, décrites en particulier par la loi de
Hooke qui s’écrit :
(7-1)
𝜎 =𝐸×𝜀
𝜎 : La contrainte appliquée (la force appliquée par unité de surface N/m2( ou pa)
𝜀 : La déformation élastique (l’allongement relatif suivant la direction d’application de la
force).
E : le module d’Young (Pa)

Figure 7-1 : Courbe Contrainte-déformation dans un solide

Dr. Bennaceur Said Page 85


CHAPITRE 07 : Ondes élastiques dans les solides

7.3 Equation d’onde dans un barreau solide


Soit un barreau solide continu, isotrope et homogène de masse volumique 𝜌 de
faible section S avec un module d’élasticité E. le barreau est soumis à une force de
traction qui exerce normalement à sa section droite sa direction est confondue avec
l’axe (ox). Ce barreau est le siège d’une onde longitudinale qui est décrite par la
fonction suivante :
⃗ (𝑥, 𝑡) = 𝑢(𝑥, 𝑡)𝑒⃗⃗⃗𝑥
𝑢 (7-2)

Les mouvements sont de faibles amplitudes, nous donnent des déformations élastiques

Figure 7-2 : Propagation d’un déplacement longitudinal dans un barreau solide

Parmi les hypothèses, on considère un élément du barreau de longueur dx à l’équilibre, de


masse dm=ρSdx de longueur dx’ en étirement.
En applique la relation fondamentale de la dynamique sur cet élément, les forces exercées
sont :
 La force exercée sur la section droite 𝐹(𝑥 + 𝑑𝑥)
 La force exercée sur la section gauche −𝐹(𝑥)
∑ 𝐹 = 𝑑𝑚 × ⃗⃗⃗⃗
𝑎𝑐 (7-3)

Pour le mouvement longitudinal :

𝜕2 𝑢
𝑎 = (𝜌𝑆𝑑𝑥) 𝜕𝑡 2 𝑒𝑥 (7-4)

Dr. Bennaceur Said Page 86


CHAPITRE 07 : Ondes élastiques dans les solides

Par la projection sur l’axe ( ox) :


𝜕2 𝑢
𝐹(𝑥 + 𝑑𝑥) − 𝐹(𝑥, 𝑡) = 𝜌𝑆𝑑𝑥 𝜕𝑡 2 (7-5)

𝑑𝑥 et infinitésimal, on peut ainsi utiliser le développement de Taylor au premier ordre :

𝜕𝐹(𝑥) 𝜕2 𝑢
𝑑𝑥 = 𝜌𝑆𝑑𝑥 𝜕𝑡 2 (7-6)
𝜕𝑥

La force de traction appliquée sur l’élément du barreau fait apparaitre une déformation
élastique 𝜀 qui vérifie la loi de Hooke :
𝐹 (7-7)
=𝐸×𝜀
𝑆
𝜕(𝑆𝐸𝜀) 𝜕 2𝑢
𝑑𝑥 = 𝜌𝑆𝑑𝑥 2 (7-8)
𝜕𝑥 𝜕𝑡
Un modèle simplifié permet d’exprimer la déformation élastique en fonction du
déplacement longitudinal.

Lors de l’étirement, le petit élément dx se trouve à la fois déplacé et déformé :

𝑑𝑥 ′ = 𝑑𝑥 + [𝑢(𝑥 + 𝑑𝑥) − 𝑢(𝑥)] (7-9)

Le développement de Taylor au premier ordre permet d’écrire :


𝜕𝑢(𝑥)
𝑢(𝑥 + 𝑑𝑥) = 𝑢(𝑥) + 𝑑𝑥 (7-10)
𝜕𝑥

𝜕𝑢(𝑥) 𝑑𝑥 ′ −𝑑𝑥
= =𝜀 (7-11)
𝜕𝑥 𝑑𝑥

𝜀 Est allongement relatif de segment dx

En remplaçant (7-5) dans (7-8), on obtient l’équation de propagation d’une onde


longitudinale dans un barreau solide :
(7-12)
𝜕𝑢2 (𝑥,𝑡) 𝐸 𝜕𝑢2 (𝑥,𝑡)
− =0
𝜕𝑡 2 𝜌 𝜕𝑥 2

Cette équation est l’équation de d’Alembert à une dimension d’une onde élastique dans
un solide.
𝐸
𝑐 = √𝜌
(7-13)
Est la vitesse de propagation de cette onde, elle dépend des caractéristiques
macroscopiques du barreau.
7.4 Résolution de l’équation de d’Alembert
La solution générale de l’équation de d’Alembert unidimensionnelle s’écrit sous la
forme :

Dr. Bennaceur Said Page 87


CHAPITRE 07 : Ondes élastiques dans les solides
𝑥 𝑥
𝑢̇ (𝑥, 𝑡) = 𝑓 (𝑡 − 𝑐 ) + 𝑔 (𝑡 + 𝑐 ) (7-14)

𝑓 𝑒𝑡 𝑔 doivent néanmoins satisfaire les conditions auxquelles est soumis le milieu de


propagation :
 Les conditions aux limites
 Les conditions initiales
 Les conditions de continuités
7.4.1 Un barreau infini
La solution de l’équation de d’Alembert est une onde plane progressive qui se propage
vers les x croissants (ou vers les x décroissants) , lorsque le barreau est considéré infini de
longueur 𝐿 ≫ 𝑆 , dans le cas des vibrations sinusoïdales, la solution est de la forme :

𝑢(𝑥, 𝑡) = 𝑢𝑚 𝑐𝑜𝑠(𝜔𝑡 ± 𝑘𝑥 + 𝜑) (7-15)


En notation
(7-16)
𝑢̅(𝑥, 𝑡) = 𝑢̅𝑚 𝑒 𝑗(𝜔𝑡±𝑘𝑥)
7.4.2 Impédance mécanique
L’impédance mécanique en point est le rapport de l’amplitude complexe de la force à
l’amplitude complexe de la vitesse de particule :
̅̅̅̅
𝜕𝑢
𝐹 (𝑥) −𝑆𝐸 ̅̅̅̅
𝑍(𝑥) = ̅
= ̅̅̅̅
𝜕𝑢
𝜕𝑥
(7-17)
𝑢
𝜕𝑡

Alors dans le cas d’une Onde plane progressive harmonique (OPPH), on obtient :
𝑆𝐸𝐾
𝑍(𝑥) =
𝜔
𝜔
La relation de dispersion :𝑘 =
𝑐
𝐸
La vitesse de propagation : 𝑐 2 =
𝜌

Alors : 𝑍(𝑥) = 𝑆√𝜌𝐸 = 𝑆𝑍𝐶 [𝐾𝑔/𝑠]


(7-18)
𝑍𝑐 = √𝜌𝐸 = 𝜌𝑐 [𝑁/𝑚2 ]
𝑍𝑐 Définit l’impédance caractéristique du barreau
7.4.3 Aspect énergétique
 Énergie cinétique :

L’énergie cinétique dEC d’une tanche élémentaire dV=Sdx dans un barreau solide lors
de passage de l’onde de est donné par :

Dr. Bennaceur Said Page 88


CHAPITRE 07 : Ondes élastiques dans les solides

1 𝜕𝑢 2 1 𝜕𝑢 2
𝑑𝐸𝐶 = 2 𝑑𝑚 ( 𝜕𝑡 ) = 2 𝜌𝑆𝑑𝑥 ( 𝜕𝑡 ) (7-19)
 L’énergie potentielle calculée à partir des forces qui interviennent lors du passage
de l’onde
1 1 𝜕𝑢 2 (7-20)
𝑑𝐸𝑝 = 𝐸𝑆𝑑𝑥(𝜀)2 = 𝐸𝑆𝑑𝑥 ( )
2 2 𝜕𝑥
Alors, dans le cas d’une OPPH qui se propage vers les x croissants, les densités
d’énergie cinétique et potentielle s’écrivent :

𝑑𝐸𝐶 1 𝜌𝑆𝑑𝑥 𝜕𝑢 2 1 𝑥
𝜀𝑐 = =2 ( 𝜕𝑡 ) = 2 𝜌𝜔2 𝑢𝑚
2
𝑠𝑖𝑛2 𝜔(𝑡 − 𝑐 )
𝑑𝜏 𝑆𝑑𝑥
(7-21)
𝑑𝐸𝑝 1 𝐸𝑆𝑑𝑥 𝜕𝑢 2 1 𝜔2 2 𝑥
𝜀𝑝 = = ( ) = 𝐸 2 𝑢𝑚 𝑠𝑖𝑛2 𝜔(𝑡 − )
𝑑𝜏 2 𝑆𝑑𝑥 𝜕𝑥 2 𝐶 𝑐
La densité d’énergie totale :
1 𝑥
𝜀𝑇 = 𝜀𝑐 + 𝜀𝑝 = 2 𝜌𝜔2 𝑢𝑚
2
(1 − 𝑐𝑜𝑠2𝜔(𝑡 − 𝑐 ) (7-22)

La densité d’énergie moyenne :


1 𝑇1 𝑥
〈𝜀𝑇 〉 = ∫0 𝜌𝜔2 𝑢𝑚
2
(1 − 𝑐𝑜𝑠2𝜔(𝑡 − 𝑐 )𝑑𝑡
𝑇 2 (7-23)

1
〈𝜀𝑇 〉 = 𝜌𝜔2 𝑢𝑚
2
2

La puissance instantanée reçue au plan x :

𝑝(𝑥, 𝑡) = 𝐹. ⃗ = −𝑆𝐸 (𝜕𝑢) (𝜕𝑢)


⃗⃗⃗ 𝑉 (7-24)
𝜕𝑥 𝜕𝑡
𝜔2 2 𝑥 𝑥
𝑝(𝑥, 𝑡) = 𝑆𝐸 𝑢𝑚 𝑠𝑖𝑛2 𝜔 (𝑡 − ) = 𝜌𝜔2 𝑐 𝑢𝑚
2
𝑠𝑖𝑛2 𝜔(𝑡 − )
𝐶 𝑐 𝑐
La puissance moyenne de l’ordre est :
1 𝑇 1
〈𝑝〉 = ∫ 𝑝(𝑥, 𝑡)𝑑𝑡 = 𝜌𝜔2 𝑐 𝑢𝑚
2
𝑆 (7-25)
𝑇 0 2
∅ = 〈𝑝〉 = 〈𝜀𝑇 〉 𝑆𝑐
∅ est le flux d’énergie qui traverse la surface S perpendiculaire à la direction de
propagation
L’intensité de l’onde est définit par la puissance moyenne par unité de surface :

〈𝑝〉 (7-26)
𝐼= == 〈𝜀𝑇 〉. 𝑉 [𝑊𝑎𝑡𝑡/𝑚2 ]
𝑆

Dr. Bennaceur Said Page 89


CHAPITRE 07 : Ondes élastiques dans les solides

Exercices
Exercice 01 :
Soit U une onde mécanique longitudinale se propageant suivant l’axe Ox dans un barreau
cylindrique homogène indéformable de masse volumique ρ, de module de
Young E de longueur l et de section droite S.
1. Ecrire l’équation de propagation d’Alembert.
2. Déterminer la solution U(x, t) de l’onde.
3. Déterminer l’impédance mécanique Z(x) à la position x. En déduire l’impédance
caractéristique du milieu Zc.
Exercice 02 :
Un barreau cylindrique d’aluminium, de masse volumique 𝜌 de longueur l et de section
droite∑, subit un allongement relatif :
△𝑙
𝐸∑
Sous l’effet d’une force F d’étirement dans le sens de l’axe 0x du barreau ; la constante E
est le module d’Young du métal. On négligera les variations de la section du barreau.
Lors du passage d’une onde acoustique, l’élément de barreau compris entre les plans de

section voisins d’abscisses x,x dx se déplacent respectivement de s (x,t) et s(x+dx,t) à

l’instant t par rapport à leur position d’équilibre ;


 Montrer que l’élongation s(x,t) obéit à une équation de propagation d’ondes.
 Calculer la célérité V de ces ondes dans le barreau d’aluminium pour lequel
𝜌=2700kg/m3 et E= 7.1010 Pa
Soit une onde plane progressive acoustique d’amplitude a0 :
𝑥
𝑠(𝑥, 𝑡) = 𝑎0 𝑐𝑜𝑠 (𝑡 − )
𝑉
 Calculer dans ce cas la variation relative maximale 𝛿𝑚𝑎𝑥 de volume de l’élément

de barreau ( x, x +dx ) entre les instants 0 et t ;

 La tension T(x, t) du barreau au niveau de sa section droite d’abscisse x,


l’instant t.
Le barreau de longueur l, fixé à une de ses extrémités O, est libre à l’autre extrémité A.
On cherche la solution de l’équation de propagation de l’onde sous la forme :
𝑠(𝑥, 𝑡) = 𝑔(𝑥)𝑠𝑖𝑛𝜔𝑡

Dr. Bennaceur Said Page 90


CHAPITRE 07 : Ondes élastiques dans les solides

 Déterminer la fonction g(x), on notera a l’amplitude de cette fonction spatiale. On


prend comme conditions aux limites :
𝜕𝑠
𝑠(𝑥 = 0, 𝑡) = 0 𝑒𝑡 (𝑥 = 𝑙, 𝑡) = 0
𝜕𝑥

Déterminer les fréquences propres du barreau en fonction de E, l, ρ et d’un entier N.


Dans les conditions où le son le plus grave a une fréquence f0=2KHz,
Déterminer :
 La longueur l= OA du barreau cylindrique.
 L’énergie cinétique moyenne 𝐸̅𝑐 (𝑡) de ce barreau en fonction de sa masse M et
de la vitesse maximale Vmax d’un élément du barreau.
Exercice 03
On se propose d’étudier la propagation d’une onde transversale à la surface S d’une
membrane tendue. On considère une membrane rectangulaire dans l’espace plan Oxyz et
dont les axes sont Ox,Oy,Oz . Soit un élément ds dont les cotes sont soumises à des
tensions linéaires τ , comme le montre la figure … comme suit :
 Etablir l’équation de propagation de l’onde sachant que la membrane a une masse
surfacique σ.
 Trouver les solutions de l’équation différentielle par la méthode des séparations
des variables.
 En déduire la forme de la solution générale de l’équation de propagation.

Dr. Bennaceur Said Page 91


Références

[1] H.DJELOUAH, Vibrations et ondes mécaniques, OPU , (2011).


[2] T.BECHERRAWY, Vibrations et Ondes, Tomes 1-4, Ed .Hermes –Lavoisier, (2010).
[3] J.BRUNEAUX , Vibrations et Ondes, Ed. Marketing,( 2010).
[4] P. DENEVE, « Mécanique », Edition ELLIPSES, ISBN 2-7298-8751-2, (1987).
[5] M. TAMINE, O. LAMROUS, « Vibrations et Ondes »,Edition OPU, ISBN 1-02-
3698, (1993).
[6] R. GABILLARD, « Vibrations et Phénomène de propagation », Edition Dunod,
(1972).
[7] M. BALKANSKI, C. SEBENE, « Ondes et phénomènes vibratoires », Edition
Dunod, (1973).
[8] C. KING, GEORGE “ Vibrations and waves” ,John Wiley & Sons, (2013).
[9] KALISKI, SYLWESTER , and L. SOLARZ . “Vibrations and waves”. Elsevier
Science, (2013).
[10] Crowell CROWELL, BENJAMIN . “Vibrations and waves”. Vol. 3. Light and
Matter, 2(000).
[11] T . BECHERRAWY, . "Vibrations, ondes et optique. Vol. 3, ondes
électromagnétiques." (2008).
[12] L. RONAN. « Ondes et vibrations: fondamentaux et applications à l'acoustique et à
la diffusion de la chaleur ». Edition Dunod, (2017).
[13] D.P.TCHOTO, and P. Mathey. “Electromagnétisme: Ondes et propagation guidée ».
Edition Dunod, (2017).
[14] FOUAD BOUKLI HACENE. «Vibrations et Ondes Mécaniques ». Université
Hassiba BENBOUALI de CHLEF. Année Universitaire : (2012 /2013).
[15] JEAN MARC RICHARD, « Ondes et Vibrations », Laboratoire de Physique
Subatomique et Cosmologie, http //lpsc.in2p3.fr/théorie/Richard, (2009).
[16] T. FREEGARDE., “Introduction to the Physics of Waves”, Cambridge University
Press (2013).

Dr. Bennaceur Said Page 92

Vous aimerez peut-être aussi