Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

telnov-mechanika-and-TO

.pdf
Скачиваний:
28
Добавлен:
28.03.2016
Размер:
2.42 Mб
Скачать

 

x +2gx + w2x = 0,

где 2g =

b

.

 

 

(56.3)

 

 

 

 

 

 

0

 

 

 

 

m

 

 

 

Как обсуждалось ранее, ищем решения в виде

 

 

 

 

 

 

 

 

 

 

 

x(t) = ReZ ,

где

Z = Aeiwt ,

 

 

 

 

 

(56.4)

w и A – пока неопределенные величины. Подставляя (56.4) в (56.3)

имеем

Aeiwt (-w2

 

 

 

 

 

 

 

 

 

 

+iw2g + w2) = 0 .

 

 

(56.5)

 

 

 

 

 

0

 

 

 

 

 

 

Поскольку левый сомножитель ненулевой, то получаем характери-

стическое уравнение

 

 

 

 

 

 

 

 

 

 

 

 

w2 -iw2g -w2= 0 .

 

 

 

 

 

(56.6)

 

 

 

 

0

 

 

 

 

 

 

 

Решение этого уравнения квадратного уравнения

 

 

 

 

 

 

 

 

w = ig -g2 + w2 .

 

 

 

 

 

(56.7)

 

 

 

 

 

0

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

Z = A1e

-gt+i w2-g2 t

+A2e

-gt-i w2-g2

t

.

(56.8)

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

Если w0

> g (трение мало)

 

 

 

 

 

 

 

 

 

Z =e-gt (A cos Wt +iA sin Wt +A cos Wt -iA sin Wt),

(56.9)

 

1

 

1

2

 

 

 

2

 

 

где

W = w2- g2 .

 

 

 

 

 

 

 

(56.10)

 

 

0

 

 

 

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

x(t) = ReZ =e-gt (c

cos Wt +c sin Wt) ,

(56.11)

 

 

 

1

 

2

 

 

 

 

 

 

где константы c1,c2 – являются вещественными (действительными)

числами, значения которых находятся из начальных условий. Действительно, из (56.9) следует, что c1 = Re(A1 +A2), c2 = Im(A2 -A1) ).

Таким образом, при малом затухании решение похоже на свободные колебания, только амплитуда затухает как e-gt , и частота несколько сдвинута: W = w20- g2 .

Время жизни колебаний (при g w0 )

141

t =

1

=

2p

´

w

0

 

»T

w

0

ºT

Q

,

(56.12)

g

w

0

2pg

 

2pg

p

 

 

 

 

 

 

 

 

 

где добротность

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q =

w0

.

 

 

 

 

 

(56.13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2g

 

 

 

 

 

 

Найдем решение задачи с учетом начальных условий x(0) = x0 , x ¢(0) = v0 . Полагая в (56.11) t = 0 сразу получаем c1 = x0 . Вторую константу находим, дифференцируя (56.11) и приравнивая v0 при

t= 0 :

x ¢(t) = -ge-gt (c1 cos Wt +c2 sin Wt) +e-gt (-c1W sin Wt +c2Wcos Wt) ,

¢

= -gc1 +c2W, откуда c2

=

v0 + gx

0

.

 

 

 

 

 

 

x (0) º v0

 

 

W

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В результате получаем искомое решение при g < w0

 

 

 

 

 

x(t) =e-gt (x

 

cos Wt +

v0

+ gx0

sin Wt).

 

(56.14)

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

W

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если w0 < g , то в (56.8) выражение под корнем становится отрица-

тельным и решение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z = A e

-gt- g2-w2

t

+A e

-gt+ g2-w2

t

,

 

(56.15)

 

 

 

 

0

 

 

 

 

 

0

 

 

 

1

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

- g2-w2

t

 

 

 

 

g2-w2

t

ö

 

 

x = ReZ =e

-gt ç

 

 

0

 

+c e

 

0

 

÷

(56.16)

 

 

 

c e

 

 

 

 

 

 

 

 

÷.

 

 

 

 

 

ç 1

 

 

 

 

 

 

 

2

 

 

 

 

÷

 

 

 

 

 

 

ç

 

 

 

 

 

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

è

 

 

 

 

 

 

 

 

 

 

 

 

 

В этом случае решение апериодическое, при затухании совершается менее одного колебания. Константы c1, 2 легко находятся из начальных

условий, аналогично предыдущему случаю (найдите сами). Наконец, рассмотрим случай g = w0 . Если подставитьg = w0 в

(56.16), то получим x = c e-gt , это неверно! Одна константа не может

удовлетворить двум начальным условиям (координата и скорость). В теории дифференциальных уравнений показывается, что решением

уравнения x +2gx + w2

= 0, при g = w

0

является

 

0

 

 

 

 

 

 

x(t) = (c

+c

t)e-gt .

(56.17)

 

1

2

 

 

 

 

142

Это можно показать следующим образом. Положим

g2 -w02 = e 0 . С учетом малости e и того, что ex » 1 +x при малом x , формулу (56.16) можно записать в виде

x(t) =e-gt (c (1 -et) +c (1 + et) =e-gt [(c

+c ) + e(c

-c )t]. (56.18)

1

2

1

2

2

1

Делая переобозначая c1 +c2 c1 ,

e(c2 -c1) c2

получаем (56.17).

Начальные условия учитываются аналогично (56.14). В результате получается

x(t) =e-gt[x

0

+(x

g +v

)t] .

(56.19)

 

0

0

 

 

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

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

Рассмотрим движение осциллятора с затуханием под действием внешней гармонической силы F(t) = F0 cos wt (для удобства, выбором

начала отсчета времени, делаем начальную фазу равной нулю). Уравнение движения

x +2gx + w2x =

F0

cos wt .

(57.1)

 

0

m

 

 

 

 

 

Как обсуждалось в §55, общим решением данного уравнения является сумма решения однородного уравнения (без правой части) и частного решения неоднородного (с правой частью) уравнения. Найдем это ча-

стное решение. Записываем силу в виде F = F0eiwt и ищем решение в виде Z = Aeiwt . Подставляя в (57.1) получаем

Aeiwt (-w2 +iw2g + w2) =

F0

eiwt ,

(57.2)

 

 

 

 

 

0

 

m

 

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A =

 

F0

 

 

.

 

(57.3)

 

m(w2-w2 +iw2g)

 

 

0

 

 

 

 

 

 

Решение неоднородного уравнения

 

 

 

 

 

 

x(t) = Re(Aeiwt ) =

F0 cos(wt + d)

= a cos(wt + d),

(57.4)

m (w2-w2 )2

+ 4g2w2

 

 

 

 

 

 

 

0

 

 

 

 

 

 

где относительная фаза между координатой и силой

143

*

Рис. 42

tg d = sin d

=

-2gw

.

(57.5)

 

cos d

 

w2-w2

 

 

 

0

 

 

Легко найти, что амплитула достигает максимума при

w* = w2

-2g2 .

(57.6)

0

 

 

Зависимость амплитуды от частоты схематически изображена на рисунке 42. Амплитуда колебаний равна:

w = 0

 

 

a =

 

F0

 

=

F0

,

 

 

 

 

(57.7)

 

mw02

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

w = w*

 

a =

 

 

 

 

 

F0

 

 

 

 

 

,

(57.8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2mg w2

- g2

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

F w2

 

w w

 

 

a =

 

0

 

=

 

0

 

0

.

(57.9)

0

 

 

mw2

 

 

 

 

 

 

 

 

 

 

 

 

k w2

 

В случае малого трения, w0

g

 

 

 

 

 

 

 

 

 

 

w* » w

 

 

a =

F0

 

=

F0

 

w0

=

F0

Q,

(57.10)

0

2mgw

0

mw2

 

2g

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

т.е. амплитуда в резонансе в Q раз (Q—добротность) больше, чем при статическом воздействии.

При g w0

a µ

1

 

 

 

 

 

 

»

1

(57.11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(w0 + w)2(w0 -w)2 + 4g2w2

2w0 Dw2 + g2

 

и при Dw2 = g2 амплитуда падает в

 

2 раз. Относительная ширина

резонансной кривой

 

 

 

 

 

 

 

 

 

 

 

 

2Dw

=

2g

=

1

.

 

(57.12)

 

 

w

 

 

 

 

 

 

0

 

w

0

Q

 

 

Исследуем теперь поведение относительной фаза d, даваемой (57.5). Эта формула записана так, что числитель равен синусу, а знаменатель косинусу этой фазы. Отсюда находим относительный сдвиг фазы координаты и силы

144

w = 0 d = 0,

w = w

0

d = -p .

(57.13)

 

2

 

w w0

 

d = -p

 

Интересно посмотреть сдвиг фазы для скорости. Если x = cos wt , то

æ

 

ö

 

ç

p÷

, то есть скорость опережает координату

x ¢ = -sin wt = cosçwt +

 

÷

ç

2

÷

 

è

ø

 

по фазе на p2 , а, значит, сдвиг фазы скорости относительно силы будет

w = 0

d = p,

 

 

 

2

,

(57.14)

w = w0

d = 0

w w0

d = -p

,

 

 

2

 

 

т.е. в резонансе скорость и сила находятся в фазе, при этом мощность, закачиваемая в осциллятор максимальна.

Общее решение уравнения колебаний под действием внешней силы есть сумма решения однородного уравнения и найденного частного решения, (57.15)

x(t) =

F0

cos(w t + d)

 

 

+(56.11), (56.16) или (56.17)

(57.16)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m (w

2-w2)2 + 4g2w2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для w0 > g , w0

< g и w0

= g , соответственно.

 

 

 

 

 

 

 

 

Рассмотрим еще отдельно случай g = 0 , при этом

 

 

 

 

 

x(t) =

F0 cos w t

 

 

+c

cos w

t +c

 

sin w

t.

(57.17)

 

 

m(w2-w2)

 

 

 

 

 

1

0

 

 

 

 

2

 

0

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть x(0) = 0,

x(0) = 0 , тогда получается c2

 

= 0 и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

 

 

 

 

ö

æ

 

 

ö

 

 

 

 

 

 

 

 

 

 

 

 

çw + w0

 

÷

çw -w0

÷

 

 

 

 

 

 

 

 

 

2F sinç

 

 

 

 

t÷sinç

 

 

t÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F (cos w t -cos w

 

t)

 

 

0

ç

 

2

 

 

÷

ç

 

2

÷

 

0

 

 

 

ç

 

 

 

÷

ç

 

÷

x(t) =

0

 

 

 

 

= -

 

 

 

è

 

 

 

 

 

ø

è

 

 

ø

.(57.18)

m(w2-w2)

 

 

 

 

 

m(w

+ w)(w -w)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

0

 

 

 

 

Левая половина формулы описывает быстрые колебания на суммарной

145

частоте, а правая медленные биения на разностной частоте, т.е. амплитуда то нарастает, то убывает. При w = w0

x µ t sin w0t .

(57.19)

амплитуда растет линейно со временем.

Потери энергии осциллятора при малом затухании

При малом затухании, w0 g

dáEñ = f

v = -bv2 = -2g mv2

= -4g mv2

= -2g E , (57.20)

dt

тр

 

2

 

 

 

 

отсюда

 

 

 

 

 

 

dáEñ = -2gáEñ

áEñ = E e-2gt .

(57.21)

 

 

dt

0

 

 

 

 

 

Работа сил при вынужденных колебаниях

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

P = -dáEñ = 2gáEñ = 2gm x2 = 2gmw2

x2

 

dt

0

. (57.22)

 

 

 

2

2

 

=

gw0F0

 

 

m [(w2-w2 )2

+ 4g2w2 ]

 

0

 

 

 

Например, при облучении атомов светом происходит рассеяние света, мощность дается формулой (57.22), где F0 =eE0 .

§ 58. Параметрический резонанс

При вынужденных колебаниях к осциллятору прикладывается внешняя сила. Но внешнее воздействие может также осуществляться путем изменения во времени параметров системы. При некоторых частотах изменения параметров могут возникать параметрический резонанс.

Рассмотрим такой резонанс на частном примере, обычных качелях.

При раскачке качелей мы внизу встаем, а в верхней точке приседаем. В верхней точке качели покоятся, и приседание на их движение не влияет. В нижней точке встаем на dl , при длине качелей l . Из сохранения момента импульса

146

 

vl

æ

 

ö

 

 

ç

 

dl ÷

 

v + dv =

 

» v ç1

+

÷.

(58.1)

 

 

l -dl

ç

 

÷

 

 

è

 

l ø

 

Найдем изменение энергии

E = mv2

,

dE = mvdv = mv2 dl

= 2E dl .

(58.2)

2

 

l

l

 

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

 

dN =

2dt

,

 

 

 

 

 

 

 

(58.3)

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

отсюда

dE

=

4dl dt

= 4

 

w

0

dl

dt ,

(58.4)

E

l

T

 

2p l

 

 

 

 

 

 

 

 

E = E0 exp(

2w0

dl

t) .

 

(58.5)

 

p

 

l

 

 

 

 

 

 

 

 

 

 

 

Поскольку есть затухание E = E e-2gt , то для роста амплитуды нуж-

 

 

 

 

0

 

 

 

 

 

 

 

 

 

но

 

 

 

 

 

 

 

 

 

 

 

 

 

2w

0

dl

> 2g

или

dl

>

gp

=

p

.

(58.6)

 

p

 

l

l

w

0

2Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

§ 59. Адиабатический инвариант

Рассмотрим следующую задачу. Пусть маленький шарик летает между неподвижной стенкой и тяжелой пластиной удаляющейся от стенки со скоростью u , упруго отражаясь от них. При этом пластина

 

X

движится очень медленно, так что за

 

один период движения шарика отно-

v

 

 

сительное изменении расстояния

 

 

 

между стенкой и пластиной мало.

 

 

 

На фазовой плоскости шарик опи-

 

 

 

сывает траекторию, показанную на

 

 

X

рис. 43 сплошной линией. Через

 

 

x

много периодов движения фазовая

 

Рис. 43

траектория изменится, как схема-

147

тично показано пунктирой линией. Как меняется площадь S внутри фазовой кривой?

Поскольку S = xv, то ее изменение

 

DS = vDx +xDv .

(59.1)

Изменение модуля скорости шарика при упругом отскакивании от

стенки удаляющеяся со скоростью u равно

 

Dv = -2u .

(59.2)

Это следует из того, что в системе стенки шарик до столкновения имеет скорость v -u , которая меняется на противоположную при отскакивании. При переходе в неподвижную систему от этой скорости нужно отнять скорость стенки, в результате, скорость шарика после отскока будет v -2u . Изменения расстояния пластины от стенки за один период движения шарика

Dx = uT = u

2x

.

(59.3)

 

 

v

 

Подставляя (59.2),(59.3) в (59.1), получаем

 

DS = v u 2x -x 2u = 0 .

(59.4)

v

 

Итак, при медленном, адиабатическом, изменении параметров системы (расстоянии между стенкой и пластиной) фазовая кривая остается (почти) замкнутой, ее форма меняется, но площадь остается постоян-

ной, является адиабатическим инвариантом.

Найдем теперь сохраняющуюся величину при медленном изменении параметров (жесткости) гармонического осциллятора.

Имеем E = mx2 2 + kx22 , где k медленно меняется. Продифференциру-

ем и произведем усреднение по большому времени (значков усреднения не пишем, но подразумеваем)

 

 

2

 

2xx

 

2mxx

 

2

 

 

 

 

 

 

 

 

 

 

kx

 

 

 

kx

+x (kx +mx) =

k E

.

(59.5)

E =

 

 

 

+k

 

+

 

=

 

k 2

2

 

2

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dE =

1 dk

ln E =

1 ln k+const

 

E

=

E

 

= const ,

(59.6)

 

 

 

 

w0

E

2 k

 

 

 

2

 

 

 

 

k

 

 

 

 

148

где w

 

=

k

, m считается постоянной. В случае математического

0

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

маятника (шарик на нитке), w2

= g

, значит инвариантом является ве-

 

 

 

 

 

0

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

личина I = E

l

. Итак, при медленном изменении параметров ос-

 

 

 

 

 

 

 

g

 

 

 

 

 

 

 

 

 

 

циллятора инвариантом является величина

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I =

 

E

.

 

 

 

 

(59.7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку E =

 

mw2x2

 

 

 

 

 

 

 

 

 

 

 

 

o 0

, а максимальная скорость v

 

= w

x

 

, то вели-

 

 

 

0

0

 

 

 

 

 

2

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

чина произведение амплитуды на скорость, т.е. площадь внутри фазовой кривой, является адиабатическим инвариантом

x0v0 = const .

(59.8)

149

Г Л А В А VII

ВОЛНЫ В УПРУГОЙ СРЕДЕ

§ 60. Упругие среды

Деформации тел под действием сил называются упругими, если при снятии напряжения их размеры восстанавливаются. При небольших относительных деформациях относительное удлинение пропорционально приложенной силе, деленной на поперечное сечение

 

T = -P = F

= E Dl

º Ee º kDl

(60.1)

 

 

 

 

 

 

 

S

 

 

 

l

 

 

S

 

F

 

 

 

F

где E

 

 

 

 

Dl

,

 

 

 

 

 

 

 

 

 

 

 

 

 

– модуль Юнга, e = l

 

 

 

Z

 

 

 

 

 

F

 

 

ES

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l . Данное соотношение на-

 

Рис. 43

 

 

 

 

k

º Dl =

 

 

 

 

 

 

 

 

 

 

 

 

 

зывается законом Гука.

 

 

 

 

 

 

 

 

При продольном растяжении (вдоль Z) поперечные размеры стерж-

ня уменьшаются

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dx = Dy

= -m Dz

º -m Dl

(60.2),

 

 

 

x

 

 

y

 

 

 

z

 

 

l

 

где m– коэффициент Пуассона.

 

 

 

 

 

 

Изменение объема при растяжении в одном направлении

 

 

DV

 

 

Dx

 

Dy

Dz

 

Dz

(1 -2m),

(60.3)

 

V

 

=

 

x

+

 

y +

 

z

 

= z

так, что при m ~ 1

объем сохраняется. Для резины m ~ 0.5 , металлов

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m ~ 0.3 , алмаза 0.07.

 

 

 

 

 

 

 

 

 

 

 

 

При равномерном обжатии

 

 

 

 

 

 

 

 

 

 

 

e

 

= e = e = -P

(1 -2m),

(60.4)

 

 

 

x

 

y

 

z

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

отсюда

150

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]