Добавил:
nastia.sokolowa2017@yandex.ru Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
ДифУр-яЧ1-2.pdf
Скачиваний:
0
Добавлен:
19.02.2024
Размер:
449.07 Кб
Скачать

Общее решение однородного уравнения таково: yо = C1ex +C2ex , а общее решение уравнения y = yн + yо = e2x (0.3cos x +0.6sin x) +C1ex +C2ex .

7. Уравнение колебаний

Рассмотрим задачу о механических колебаниях. Пусть груз массы m покоится на упругой рессоре, закрепленной в некоторой точке A (рис. 8.1).

m

Положение равновесия

y

A A

Рис. 8.1

Отклонение груза от положения равновесия обозначим переменной y. Будем считать положительным отклонение вниз и отрицательным – отклонение вверх. В положении равновесия сила веса уравновешивается

упругостью рессоры. Предположим, что сила F1 , стремящаяся вернуть груз в положение равновесия (так называемая восстанавливающая сила), пропорциональна отклонению, т. е.

F1 = −ky (k = const, k > 0)

(1)

(k – положительная постоянная для данной рессоры, так называемая жесткость рессоры). Предположим еще, что движению груза препятствует сила сопротивления F2 , направленная в сторону, противоположную направ-

лению движения и пропорциональная скорости v движения груза, т.е.

 

 

F

= −λv = −λ

d y

(λ = const, λ > 0) .

 

(2)

 

 

 

 

 

 

2

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Составим дифференциальное уравнение движения груза на рессоре. В соответствии со вторым законом

Ньютона получаем mw = F

+ F

, где

w =

d 2 y

, или

 

 

 

 

 

 

 

 

 

1

2

 

 

 

dt2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

d 2 y

= −k y λ

d y

(k > 0, λ

> 0) .

(3)

 

 

 

 

 

dt2

dt

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение (3) является линейным дифференциальным уравнением второго порядка с постоянными коэф-

фициентами, его можно записать в виде

 

 

 

 

 

 

 

 

 

d 2 y

+ p

d y

+ q y = 0

,

(4)

 

dt2

dt

 

 

 

 

 

 

 

 

где

λ

 

 

 

k

 

 

 

 

p =

,

q =

.

 

(5)

 

m

 

 

 

 

 

 

 

m

 

 

Уравнение (4) называется уравнением свободных колебаний.

Рассмотрим случай, когда точка A прикрепления рессоры совершает вертикальное движение по закону z = ϕ(t) .

Это имеет место тогда, когда нижний конец рессоры прикреплен к катку, движущемуся вместе с рессорой и грузом по неровности (рис. 8.2).

При этом восстанавливающая сила F1 будет выражаться формулой F1 = −k(y +ϕ(t)), а сила сопротивления F2 – формулой F2 = −λ(y′+ϕ(t)).

45

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

z = ϕ(t)

 

 

 

 

 

 

 

 

t

 

 

 

 

Рис. 2

 

 

 

Дифференциальное уравнение движения груза на рессоре принимает вид

 

m

d 2 y

= −k

(y +ϕ(t))λ(y

(6)

dt2

 

+ϕ (t))

или

 

 

 

 

 

 

 

 

 

 

d 2 y

+ p

d y

+ q y =

f (t) ,

(7)

 

 

dt2

dt

 

 

 

 

 

 

 

где

= λ

 

k ,

f (t) = − kϕ(t) + λϕ(t) .

 

p

, q =

(8)

 

m

 

m

 

 

 

m

 

Уравнение (7) является неоднородным линейным дифференциальным уравнением второго порядка с по-

стоянными коэффициентами. Это уравнение называется уравнением вынужденных колебаний.

 

7.1. Свободные колебания

Рассмотрим уравнение свободных колебаний, т. е. уравнение

 

d 2 y

+ p

d y

+ q y =

0 .

 

Соответствующее

 

dt

2

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

характеристическое уравнение k2 + pk + q = 0 имеет корни k = −

p

 

 

p2

 

q,

k

2

= −

p

+

 

p2

q , в зависи-

 

 

 

 

 

 

 

1

2

 

 

 

 

 

4

 

 

 

 

 

2

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

мости от которых получаем общее решение уравнения свободных колебаний.

 

 

F2 = 0 ,

тогда

λ = 0 и p = 0

Рассмотрим сначала случай, когда отсутствует сила сопротивления, т.

е.

(см. (2) и (5)). Характеристическое уравнение в этом случае принимает вид k2 + q = 0 и имеет комплексные

корни k1 = − iβ,

k2 = iβ , где β =

q

, поэтому общее решение уравнения свободных колебаний

y′′+ q y = 0

определяется формулой

 

 

 

 

 

y = C1 cos βt +C2 sin βt .

(1.1)

 

 

 

 

 

 

 

 

 

Преобразуя правую часть этой формулы, введя новые постоянные A и ϕ0 : C1 = Asinϕ0 , C2

= Acosϕ0 ,

 

 

 

 

 

= arctg

C1

.

 

откуда A = C2

+C2 , ϕ

0

 

 

 

1

2

 

 

C2

 

 

 

 

 

 

 

 

Подставляя выражения для C1 и C2 в формулу (1.1), получаем общее решение

 

 

 

 

 

 

 

 

 

 

y = Asin(βt +ϕ0 ) ,

(1.2)

где A и ϕ0 – произвольные постоянные величины. Колебания в данном случае называются гармоническими. Интегральными линиями являются синусоиды.

Определение 1. Периодом колебаний называется промежуток времени T, за который аргумент синуса изменится на 2π .

В рассматриваемом случае T = 2βπ .

Определение 2. Частотой колебаний называется число колебаний за время 2π .

В данном случае частота равна β .

Определение 3. Амплитуда колебаний есть величина наибольшего отклонения от положения равновесия.

46

 

 

Из формулы (1.2) видно, что это наибольшее отклонение равно A .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определение 4. Начальной фазой называется величина ϕ0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Исследуем случай, когда F 0, т. е.

з 0 . Рассмотрим три возможности:

 

p2

 

< q,

p2

> q,

p2

= q .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

4

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если

p2

< q , то характеристическое уравнение k2 + pk + q = 0 имеет комплексные корни

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k =α +iβ, k

2

=α iβ , где α = −

 

p

<

0, β = q

p2

.

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

Общее решение уравнения (4) определяется формулой

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = eαt (C

cos βt +C

2

sin βt)

 

 

 

 

 

 

 

 

 

(1.3)

или

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = Aeαt sin(βt +ϕ0 ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.4)

 

 

В качестве амплитуды здесь рассматривают величину

 

Aeαt

,

зависящую от времени. Поскольку α < 0 ,

то эта величина стремится к нулю при t → ∞, и колебания называются затухающими.

 

 

 

 

 

 

Если

p2

> q , то характеристическое уравнение имеет вещественные отрицательные различные корни k

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и k2 . Общее решение уравнения определяется формулой

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = C ek1t +C

ek2t

(k

< 0, k

2

< 0) ,

 

 

 

 

 

 

 

 

 

 

(1.5)

 

 

 

 

 

 

1

 

2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

из которой видно, что y 0 при t → ∞ . Колебаний в данном случае не будет (сила сопротивления велика по сравнению с жесткостью рессоры).

Если

p2

= q , то корни характеристического уравнения равны между собой

(k

= k

2

= −

p

< 0)

, поэтому

 

 

4

 

 

 

 

 

1

 

2

 

 

общее решение уравнения (4) определяется формулой

 

 

 

 

 

 

 

 

 

 

 

 

 

y = (C

+C

t)e

p

t .

 

 

 

 

 

 

(1.6)

 

 

2

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

Здесь также y 0 при t → ∞ , однако, y стремится медленнее, чем в предыдущем случае.

7.2. Вынужденные колебания

Обратимся к уравнению вынужденных колебаний, т. е. к уравнению (7). Рассмотрим практически важный случай, когда возмущающая внешняя сила является периодической и изменяется по закону f (t) = asinωt .

В этом случае уравнение (7) запишется так:

 

d 2 y

+ p

d y

+ q y = asinωt .

 

(2.1)

 

dt2

dt

 

 

 

 

 

 

Исследуем уравнение (2.1) при p = 0 (сила сопротивления отсутствует) и p 0,

p2

< q .

4

 

 

 

 

 

 

Если сила сопротивления отсутствует, т. е. p = 0 , то уравнение вынужденных колебаний принимает вид

d 2t + q y = asinωt . dt2

Соответствующее однородное уравнение имеет общее решение

yн =C1 cos βt +C2 sin βt (β2 = q)

или

yн = Asin(βt +ϕ0 )

Предположим, что частота внешней силы не равна частоте собственных колебаний, т. е.

частное решение неоднородного уравнения (2.2) имеет вид

yн = P cosωt +Q sinωt ,

Подставляя эту функцию и вторую производную в уравнение (2.2), находим коэффициенты

P = 0, Q = q aω2 .

(2.2)

(2.3)

(2.4)

β ω , тогда

(2.5)

P и Q :

47

Следовательно,

 

 

a

 

 

 

y

н

=

sinωt .

(2.6)

q ω2

 

 

 

 

 

Общее решение уравнения (2.2) определяется формулой y = yо + yн , т. е.

 

y = Asin(βt +ϕ0 ) +

 

a

sinωt .

(2.7)

 

q ω2

 

 

 

 

 

 

 

Если частота собственных колебаний совпадает с частотой внешней силы, т. е. β = ω , то частное решение уравнения (2.2) следует искать в следующем виде:

 

 

 

yн =t (P cosωt +Q sinωt) .

(2.8)

Подставляя функцию (2.8) и ее вторую производную в уравнение (2.2), находим

 

 

a

 

~

 

a

 

P = −

 

,

Q = 0, т. е. y

= −

 

t cosωt .

 

2ω

2ω

 

Итак, общее решение уравнения (2.2) в этом случае имеет вид y = Asin(βt +ϕ0 ) 2aω t cosωt , или

y = Asin(βt +ϕ0 )

a

t cos βt ,

(2.9)

2β

 

 

 

поскольку ω = β .

Из последней формулы видно (см. второе слагаемое алгебраической суммы), что амплитуда колебания неограниченно возрастает при неограниченном возрастании времени t . Рассматриваемое явление (частота собственных колебаний совпадает с частотой внешней силы) называется резонансом.

При выполнении условий p 0 и p42 < q характеристическое уравнение k2 + pk + q = 0 имеет ком- плексно-сопряженные корни k1 =α +iβ, k2 =α iβ .

Общее решение однородного уравнения y′′+ p y′+ q y = 0 , где y = y(t) , определяется формулой

y

о

= eαt (C

cos βt +C

2

sin βt) или y

о

= Aeαt

sin(βt +ϕ

0

) .

(2.10)

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Частное решение неоднородного уравнения (2.1) ищем в виде функции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yн = P cosωt +Q sinωt .

 

 

 

 

 

 

 

(2.11)

Подставляя функцию (2.11) и ее производные yо, yо′′

в уравнение (2.1), находим коэффициенты P и Q:

 

 

 

P =

 

 

pωa

 

 

,

Q =

 

(q ω2 )a

 

 

 

.

 

 

(2.12)

 

 

(q ω

2

)

2

2

 

2

(q ω

2

)

2

+ p

2

2

 

 

 

 

 

 

 

 

 

+ p ω

 

 

 

 

 

 

 

ω

 

 

 

 

 

 

 

Вводим новые постоянные величины

 

 

A1 и ϕ1 , полагая

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P = A1 sinϕ1,

 

Q = A1 cosϕ1 ,

 

 

 

 

 

 

 

 

 

 

 

 

(2.13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

, tgϕ =

P

.

 

 

 

 

 

A = P2 +Q2

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

(q

ω2 )2 + p2ω2

 

 

 

 

 

1

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Формула (2.11) с учетом равенств (2.12) и (2.13) принимает вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yн =

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

sin(ωt +ϕ1) .

 

 

(2.14)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(q ω

2 )2

+ p2ω2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общее решение уравнения (2.1) определяется формулой y = yо + yн

или

 

 

 

y = Aeαt sin(βt +ϕ

0

) +

 

 

 

 

 

 

a

 

 

 

 

 

 

sin(ωt +ϕ ) .

 

 

(2.15)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(q ω2 )2 + p2ω2

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, отклонение y состоит из суммы двух слагаемых. Первое из них (решение соответствующего однородного уравнения) определяет затухающие колебания (α = − 2p < 0, ибо p > 0) , это слагаемое убывает

при возрастании t. Следовательно, через некоторый промежуток времени главное значениеy будет определяться вторым слагаемым, представляющим вынужденные колебания. Частотаω этих колебаний равна частоте внешней

силы: f (t) = asinωt . Амплитуда вынужденных колебаний тем больше, чемменьше p и разность q ω2 .

48