Тут r – стала, яку ми будемо називати коефіцієнтом тертя. Знак мінус обумовлений тим,
що сила F й швидкість υ спрямовані у протилежні сторони, внаслідок чого їх проекції на вісь X мають різні знаки.
Рівняння другого закону Ньютона за наявності сили тертя має вигляд |
(31.2) |
|||||||
|
|
|
mx = −kx − rx . |
|||||
|
|
|
&& |
|
& |
|
|
|
Використаємо позначення |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
β = r /(2m) |
, |
ω02 = k / m |
, |
(31.3) |
|||
і напишемо рівняння (31.2) у вигляді |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
&& |
& |
|
2 |
. |
(31.4) |
|
|
|
x |
+ 2βx + ω0 x = 0 |
|||||
Відзначимо, що величину β в (31.4) називають коефіцієнтом загасання; ω0 |
є власною |
|||||||
частотою коливальної системи, тобто та частота, з якої коливалася б система за умови відсутності тертя. Таким чином, отримали диференціальне рівняння (31.4), що визначає поведінку коливальної величини x за наявності сили тертя. Це рівняння називають
диференціальним рівнянням загасаючих коливань.
2 Розглянемо загасаючі електричні коливання. Нехай у коливальному контурі, крім ємності C й індуктивності L , є активний опір R (рис. 31.1). Застосуємо закон Ома для
ділянки кола 1-3-2 (див. рис. 31.1): |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
IR = ϕ1 − ϕ2 + Es . |
|
|
|
|
|
(31.5) |
|||||||||||
Різницю потенціалів на конденсаторі визначимо зі співвідношення |
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
ϕ1 − ϕ2 = q1 / C = (− q)/ C . |
|
|
|
|
|
(31.6) |
|||||||||||||
Тут використали, що заряд пластини конденсатора 1 q1 = −q (див. рис. 31.1). |
|
|
|
|
|
|||||||||||||||||
Сила струму |
I є додатною, |
|
коли напрям струму |
|
+ q |
C |
|
|
|
|
|
|||||||||||
збігається з напрямом обходу ділянки кола 1-3-2, тобто за |
|
|
|
− q |
||||||||||||||||||
|
|
|
||||||||||||||||||||
годинниковою стрілкою. |
У цьому разі заряд на пластині |
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
2 |
|
|
1 |
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||
конденсатора q2 = q пов’язаний із силою струму в ділянці |
|
|
|
|
|
|
R |
|||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||
кола таким співвідношенням: |
|
|
|
|
|
|
|
|
|
I |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
I = +dq / dt = +q . |
|
|
(31.7) |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
& |
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|||
Знак «+» обумовлений тим, що, коли струм I |
|
є додатним, |
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
L |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||
заряд q2 = q збільшується ( q > 0 ). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
& |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Підставимо |
в |
(31.5) |
|
закон |
самоіндукції |
|
|
Рисунок 31.1 |
||||||||||||||
Es = −L dI / dt , співвідношення (31.6) й (31.7) й отримаємо |
|
|
|
|
|
(31.8) |
||||||||||||||||
|
|
|
|
|
Rq = −q / C − L q . |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
& |
|
&& |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Далі вводимо позначення |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
β = R /(2L) |
, |
|
ω02 = 1/(LC) |
, |
|
|
|
|
|
|
(31.9) |
||||||||
і перетворюємо рівняння (31.8) до такого вигляду: |
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
&& |
& |
2 |
. |
|
|
|
|
|
|
(31.10) |
|||||||
|
|
|
|
|
q |
+ 2βq |
+ ω0q = 0 |
|
|
|
|
|
|
|||||||||
Отримане рівняння подібне до (31.4), яке описує механічні коливання, і має таку саму назву:
диференціальне рівняння загасаючих коливань. Величину β в (31.10), як і у випадку механічної системи, називають коефіцієнтом загасання; ω0 є власною частотою
контуру, тобто та частота, з якої відбувалися б коливання за умови відсутності активного опору. Таким чином, отримали диференціальне рівняння (31.10), що визначає поведінку коливальної величини q за наявності активного опору в коливальному контурі.
66
Із порівняння формул (31.3) і (31.9) випливає, що опір R відіграє роль коефіцієнта тертя r , індуктивність L – роль маси, величина, що зворотна ємності C , – роль коефіцієнта квазіпружної сили k .
§ 32 Розв’язання диференціального рівняння загасаючих коливань. Коефіцієнт загасання, декремент загасання, логарифмічний декремент загасання, добротність [5]
1 Знайдемо розв’язок диференціального рівняння гармонічних коливань: |
|
||
&& |
& |
2 |
(32.1) |
x |
+ 2βx + ω0 x = 0 . |
||
У цьому рівнянні β – коефіцієнт загасання; ω0 |
– власна частота коливальної системи (тобто |
||
та частота, з якою коливалася б система за умови відсутності загасання). Величина x може бути механічним зміщенням частинки, електричним зарядом на конденсаторі й т.д.
Коефіцієнти β та ω0 |
визначаються параметрами коливальної системи. |
|
|
|
|
|
||||||||||||||||||||||||
Будемо шукати розв’язок рівняння (32.1) у вигляді |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
x = u exp(−βt), |
|
|
|
|
|
|
|
|
|
|
|
(32.2) |
|||||
де u – деяка, поки що |
невідома, |
функція від |
|
t . Диференціювання функції x (32.2) за |
||||||||||||||||||||||||||
змінною t |
дає |
|
|
|
x = u exp(−βt)− uβexp(− βt)= (u −βu)exp(−βt), |
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
& |
|
& |
|
|
|
|
|
|
|
|
|
& |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
u)exp(− βt). |
|
|
|||||
|
x = (u −βu)exp(−βt)− (u − βu)βexp(− βt)= (u − 2βu + β |
|
|
|||||||||||||||||||||||||||
|
&& |
&& |
|
& |
|
|
|
|
x |
& |
|
|
|
|
&& |
|
& |
|
|
|
|
|
|
|
|
|
||||
Після підстановки виразів для x , |
і x |
у рівняння (32.1) і скорочення на відмінний від нуля |
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
& |
|
&& |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
множник exp(− βt) отримаємо диференціальне рівняння для u : |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
&& |
2 |
−β |
2 |
)u = 0 . |
|
|
|
|
|
|
|
|
|
|
|
(32.3) |
||
|
|
|
|
|
|
|
|
|
|
|
|
u |
+ (ω0 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Розв’язання рівняння |
(32.3) |
залежить від |
знака |
|
|
X |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
коефіцієнта, що стоїть біля u . Розглянемо |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
A0 exp(−βt) |
|||||||||||||||||||||
випадок, коли цей коефіцієнт є додатним (тобто |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
β < ω0 ). Введемо позначення |
|
|
|
|
|
|
|
|
|
A0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A′ |
|
|
′′ |
|
A′′′ |
|
|
|
|
ω = |
ω0 − β |
, |
|
|
|
(32.4) |
|
0 |
|
|
|
|
|
|
|
A |
|
t |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
і отримаємо рівняння |
&& |
|
2 |
u |
= 0 . |
|
|
|
|
|
(32.5) |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
u |
+ ω |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Рівняння |
(32.5) |
є |
диференціальним |
рівняння |
|
|
|
|
|
|
|
T |
|
|
|
|
|
|||||||||||||
гармонічних коливань і тому його розв’язок |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
можемо записати у вигляді |
|
|
|
|
|
|
|
|
|
|
Рисунок 32.1 – Графік |
загасаючого |
||||||||||||||||||
|
|
u = A0 cos(ωt + α), |
|
|
|
|
|
коливання. Верхня штрихова крива – |
||||||||||||||||||||||
де ω є частотою загасаючих коливань (див. |
графік зміни амплітуди з часом |
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
також (32.4)). Далі підставляємо отриманий вираз для u |
в (32.2) і знаходимо у випадку |
|||||||||||||||||||||||||||||
малого тертя (β < ω0 ) розв’язок диференціального рівняння загасаючих коливань (32.1): |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
x = A e−βt cos(ωt + α) |
. |
|
|
|
|
|
|
|
|
|
|
(32.6) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Тут A0 і |
α – |
сталі, |
значення |
яких залежать |
від початкових умов, ω – |
величина, що |
||||||||||||||||||||||||
визначається формулою (32.4). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
67
2 Проведемо дослідження отриманого результату (32.6), з’ясуємо характеристики загасаючих коливань. Графік функції (32.6) наведений на рис. 32.1. Штриховими лініями показані межі, у яких знаходиться зміщення змінної величини x .
Відповідно до виду функції (32.6) величину x можна розглядати як гармонічне коливання частоти ω з амплітудою, що змінюється за законом
A(t) = A e−βt |
. |
(32.7) |
0 |
Величину (32.7) називають амплітудою загасаючих коливань. Верхня зі штрихових кривих на рис. 32.1 дає графік цієї функції, причому A0 є амплітудою в початковий момент часу.
Розглянемо послідовні найбільші відхилення величини x (вони відбуваються через період загасаючих коливань T ), наприклад, A′, A′′ , A′′′ і т.д. на рис. 32.1. Неважко з’ясувати, що відношення двох послідовних найбільших відхилень мають одне й те саме значення. Дійсно, коли A′ = A0 exp(-bt), то
|
A′ |
= |
|
A0 exp[-b(t)] |
|
= exp(bT ), |
||
|
¢¢ |
|
A0 exp[-b(t +T )] |
|||||
|
A |
|
|
|
|
|
||
|
A′′ |
= |
|
A0 exp[-b(t +T )] |
|
= exp(bT ) |
||
¢¢¢ |
|
A0 exp[-b(t + 2T )] |
||||||
A |
|
|
|
|||||
і т.д. Це означає, що таке відношення може бути характеристикою загасаючого коливання.
Таким чином, відношення значень амплітуди, що відповідають моментам часу, що відрізняється на період, дорівнює
|
|
|
A(t) |
|
|
= exp(bT ). |
|
||||||
|
|
|
A(t +T ) |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
||||
Це |
відношення називають декрементом |
загасання, а його |
логарифм – логарифмічним |
||||||||||
декрементом загасання: |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
æ |
|
A(t) |
ö |
|
|
|
|
||||
|
|
l = lnç |
|
|
|
|
|
÷ = bT |
. |
(32.8) |
|||
|
|
|
|
|
|
||||||||
|
|
|
èç A(t +T )ø÷ |
|
|
||||||||
|
Відповідно до формули (32.4) період загасаючих коливань відрізняється від періоду |
||||||||||||
вільних коливань: |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
||||||||
|
|
|
T = 2p/ w02 -b2 |
. |
(32.9) |
||||||||
При |
незначному терті ( b2 << w2 ) період |
коливань практично |
дорівнює періоду вільних |
||||||||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
коливань T0 = 2p/ w0 . Зі збільшенням коефіцієнта загасання період коливань зростає. Для характеристики коливальної системи використовується також величина
|
|
|
Q = π / λ |
, |
(32.10) |
яка називається добротністю коливальної системи.
З'ясуємо фізичний зміст коефіцієнта загасання β . Знайдемо час τ , за який амплітуда
зменшується в e раз. З визначення величини τ випливає, що e−βτ = e−1 , звідки βτ =1. Отже,
коефіцієнт загасання β є оберненим до проміжку часу τ , за який амплітуда зменшується в e разів (β = 1/ τ ). У цьому й полягає фізичний зміст коефіцієнта загасання.
З’ясуємо фізичний зміст логарифмічного декременту загасання λ . Виразивши відповідно до (32.8) β через λ і T , можна закон зменшення амплітуди з часом написати у
вигляді
A(t) = A0e−(λ /T )t .
68
Час τ , за який амплітуда зменшується в e раз, система встигає виконати Ne = τ /T коливань.
З умови exp(− λτ /T )= exp(−1) маємо, що λτ /T = λNe =1. Отже, логарифмічний декремент загасання є оберненим до числа коливань, за час яких амплітуда зменшується в e разів ( λ =1/ Ne ). У цьому полягає фізичний зміст логарифмічного декремента загасання.
Фізичний зміст добротності полягає у тому, що вона прямо пропорційна числу коливань, за час яких амплітуда зменшується в e разів (Q = π / λ = πNe ).
§ 33 Диференціальне рівняння вимушених коливань та його розв’язання [5]
1 Розглянемо механічну коливальну систему із загасанням, яка знаходиться під дією зовнішньої сили, що змінюється з часом за гармонічним законом:
Fx = F0 cosΩt . (33.1)
Під дією зовнішньої періодичної сили в системі виникають вимушені коливання. Знайдемо диференціальне рівняння, яке описує вимушені коливання. Для цього застосуємо другий закон Ньютона:
|
|
&& |
|
& |
|
|
|
|
|
||
|
|
|
mx |
= −kx − rx + F0 cosΩt . |
|
||||||
Увівши позначення |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
β = r /(2m) |
, |
ω02 = k / m |
, |
|
|
|||
перетворимо рівняння до такого вигляду: |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
&& |
|
|
& |
2 |
|
(F0 / m)cosΩt |
. |
(33.2) |
||
|
x |
+ 2βx + ω0 x = |
|||||||||
Тут β – коефіцієнт загасання; |
ω0 |
– власна |
|
частота коливальної системи; Ω |
– частота |
||||||
зовнішньої періодичної сили. Рівняння (33.2) описує вимушені коливання й називається
диференціальним рівнянням вимушених коливань.
2 Розглянемо вимушені електричні коливання у коливальному контурі з активним опором. Підключимо до коливального контуру з ємністю C , індуктивністю L й активним
опором R зовнішнє джерело змінної напруги: |
|
U = Um cosΩt |
(33.3) |
(див. рис. 33.1). Під дією зовнішньої змінної напруги у контурі виникають вимушені коливання. Отримаємо диференціальне рівняння, яке описує процеси у контурі. Для цього застосуємо закон Ома для ділянки кола 1-3-2 (див. рис. 33.1):
|
IR = ϕ1 − ϕ2 + E . |
|
|
|
(33.4) |
||||||||
Слід зазначити, що змінну напругу зовнішнього джерела U тут потрібно враховувати разом |
|||||||||||||
з ЕРС самоіндукції. Тобто загальна ЕРС, яка діє в контурі, дорівнює |
|
|
|
|
|
|
|
|
|
||||
|
E = Es +U = −L dI / dt +U . |
|
|
(33.5) |
|||||||||
Різницю потенціалів на конденсаторі визначимо зі |
I |
+ q |
|
C |
|||||||||
співвідношення |
|
|
|
|
|
|
− q |
||||||
|
|
|
|
|
|
||||||||
ϕ1 − ϕ2 = q1 / C = (− q)/ C . |
(33.6) |
|
|
|
|
|
|
|
|
|
|
||
|
2 |
|
|
|
1 |
|
|
|
|
||||
|
|
|
|
|
|
|
|||||||
Тут використали, що |
заряд пластини |
конденсатора |
1 |
U |
|
|
|
|
|
|
|
|
R |
|
|
|
|
|
|
|
|
||||||
q1 = −q (див. рис. 33.1). |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
||
Сила струму I |
є додатною, коли напрям струму |
|
|
3 |
|
|
|
|
|
||||
збігається з напрямом обходу ділянки кола 1-3-2, тобто за |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
||||
годинниковою стрілкою. У цьому разі заряд на пластині |
|
|
L |
||||||||||
конденсатора q2 = q |
пов’язаний із силою струму |
в |
|
Рисунок 33.1 |
|||||||||
|
|
69 |
|
|
|
|
|
|
|
|
|
|
|
ділянці кола таким співвідношенням:
I = +dq / dt = +q . |
(33.7) |
& |
|
Знак «+» обумовлений тим, що, коли струм I є додатним, заряд q2 = q збільшується ( q& > 0 ).
Підставимо у (33.4) співвідношення (33.5)-(33.6) з урахуванням (33.7) та (33.3) й отримуємо
|
Rq = −q / C − L q +Um cosΩt . |
(33.8) |
||||||
|
|
& |
|
|
&& |
|
|
|
Далі вводимо позначення |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
β = R /(2L) |
, |
ω02 = 1/(LC) |
, |
|
(33.9) |
|
і перетворюємо рівняння (33.8) до такого вигляду: |
|
|||||||
|
|
|
|
|
|
|
|
|
|
&& |
& |
2 |
|
|
|
. |
(33.10) |
|
q |
+ 2βq + ω0q = (Um / L)cosΩt |
||||||
Тут β – коефіцієнт загасання; ω0 – власна частота коливального контуру; Ω – частота
коливань зовнішнього джерела. Рівняння (33.10) описує вимушені коливання й називається
диференціальним рівнянням вимушених коливань.
Порівнявши диференціальне рівняння вимушених коливань для механічної системи (33.2) та для електричного коливального контуру (33.10), можемо зробити висновок, що вони є з математичної точки зору однаковими.
3 Знайдемо розв’язок диференціального рівняння вимушених коливань (33.10) (такий самий розв'язок буде й для рівняння (33.2)).
Рівняння типу (33.10) називають неоднорідними диференціальними рівняннями з сталими коефіцієнтами. З теорії лінійних диференціальних рівнянь із сталими коефіцієнтами відомо, що загальний розв’язок неоднорідного рівняння (тобто рівняння, у правій частині якого стоїть функція від t , яка не дорівнює тотожно нулю) дорівнює сумі загального розв’язку відповідного однорідного рівняння (тобто того ж рівняння, у якому права частина дорівнює тотожно нулю) і частинного розв’язку неоднорідного рівняння. Загальний розв’язок однорідного рівняння ми вже знаємо (розв’язок диференціального рівняння загасаючих коливань). Воно має вигляд
q = q0e−βt cos(ωt + α), |
(33.11) |
|
|
|
|
де ω = ω2 |
− β2 |
– частота загасаючих коливань. |
|
0 |
|
|
|
Залишається тепер знайти частинний (який не має довільних сталих) розв’язок рівняння (33.10). Будемо шукати цей розв’язок у вигляді
q = Acos(Ωt − ϕ) , |
(33.12) |
де ϕ – зсув фаз між зовнішньою напругою і викликаними нею коливаннями в контурі. Спробуємо з'ясувати, чи не існує таких значень A і ϕ , при яких функція (33.12) задовольняє рівняння (33.10). Для цього підставимо у рівняння (33.10) вираз (33.12) і його похідні:
q = −AΩsin(Ωt − ϕ), |
(33.13) |
|||
& |
|
|
|
|
&& |
= −AΩ |
2 |
cos(Ωt − ϕ), |
(33.14) |
q |
|
|||
розвертаючи одночасно sin(Ωt − ϕ) |
й cos(Ωt − ϕ) за формулами |
для синуса й косинуса |
||
різниці: |
|
|
|
|
− AΩ2 [cosϕcosΩt + sin ϕsin Ωt]− 2βAΩ[cosϕsin Ωt − sin ϕcosΩt]+ + ω20 A[cosϕcosΩt + sin ϕsin Ωt]= (Um / L)cosΩt .
Згрупувавши відповідним чином члени рівняння, отримаємо
70