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

книги / Принципы и практика решения задач по общей физике. Механика. Физика макросистем

.pdf
Скачиваний:
8
Добавлен:
12.11.2023
Размер:
3.47 Mб
Скачать

может перемещаться без трения по горизонтальной плоскости, упруго закреплен с помощью вертикальной пружины с коэффициентом жесткости k (рис. 3.3). В положении

равновесия пружина не растянута и ее длина равна l0 . Для определения потенциальной

энергии пружины выразим ее удлинение через смещение груза x :

 

 

 

l = l

2 + x2 l .

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

k

 

2 2

 

 

2

 

 

U (x) =

 

kl =

 

 

( l0 + x l0 ) .

 

Рис. 3.3

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Предполагая

колебания

малыми

( x << l0 ), последнее выражение

можно представить в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U (x) =

kx

4

= αx

4

,

 

 

 

 

 

 

 

 

 

8l

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

где α =

k

.

Заметим, что дифференциальное уравнение колебаний

8l 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

(второй закон Ньютона) оказывается нелинейным. В самом деле сила, действующая на груз,

F(x) = −

dU

 

= −

kx3

.

dx

 

2l 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

И уравнение колебаний приобретет вид

 

 

m

d 2 x

+

 

k

 

 

x

3

= 0 .

 

dt2

 

2l

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

Таким образом, из предположения малости колебаний не обязательно следует линейное описание системы!

121

Для определения периода колебаний подставим U (x) в выражение (2):

A

 

 

dx

 

 

 

 

T = 4

 

 

 

 

 

 

.

 

2α

 

 

 

 

 

 

0

 

 

 

4

x

4

 

 

 

 

 

 

 

m

A

 

 

 

 

 

 

 

 

 

 

 

 

Введем замену переменных u = x / A . Тогда выражение для периода можно записать в виде

 

4

1

du

 

T =

m

0

.

A

2α

1 u4

Этот интеграл выражается через гамма-функцию Г(z) :

 

 

 

 

Г(

5

 

1

du

 

 

4)

 

0

 

=

π

 

 

1,311 .

1 u4

Г(

3

 

 

 

 

4)

 

И для периода колебаний получаем окончательное выражение:

T = 3,71 mk lA0 .

В этом выражении можно видеть проявление одного из характерных свойств нелинейных систем – их неизохронность, т.е. зависимость периода свободных колебаний от амплитуды.

Еще один пример нелинейной колебательной системы – физический и математический маятники. Колебания физических и математических маятников являются изохронными (т.е. период колебаний не зависит от амплитуды) только при достаточно малых амплитудах. При увеличении амплитуды колебаний они становятся неизохронными и, кроме того, не являются чисто гармоническими (в них появляются гармоники более высоких порядков).

Для установления зависимости периода колебаний (и соответственно частоты) от амплитуды применим энергетический подход к незатухающим колебаниям физического маятника.

122

 

1

dϕ 2

Из закона сохранения энергии следует

 

I

 

 

=

2

dt

 

 

 

 

= mga(cos ϕ−cos ϕ0 ) , где ϕ – угол отклонения маятника от положе-

ния равновесия; ϕ0

– максимальное значение угла отклонения (ам-

плитуда колебаний); I – момент инерции маятника; a

– расстояние

от точки подвеса до центра инерции маятника (рис. 3.4).

 

Введем в рассмотрение период колебаний фи-

 

зического маятника в приближении малых амплитуд

 

 

 

 

T

= 2π

I

.

 

 

 

 

 

 

 

 

 

0

 

mga

 

 

 

 

 

 

 

Тогда закон сохранения энергии можно будет

 

представить в виде

 

 

 

 

 

 

 

dϕ

=

4π

sin2 ϕ0

sin2 ϕ .

 

 

 

T0

 

 

dt

2

2

Рис. 3.4

Для определения периода колебаний разрешим

это уравнение относительно dt , проинтегрируем его по углу от ϕ = 0

до ϕ = ϕ0

и умножим результат на четыре. В итоге получаем

 

 

T

ϕ0

 

 

 

dϕ

 

 

 

 

 

 

 

 

 

T =

 

0

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

π

 

sin2

ϕ0

sin2 ϕ

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

Путем последовательной

замены

 

переменных x =

sin(ϕ/ 2)

 

 

k

и x = sin u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

полученное выражение

для

периода колебаний можно

привести к полному эллиптическому интегралу первого рода:

 

 

 

 

2T

 

π/ 2

 

du

 

 

 

 

 

 

 

 

T =

 

 

0

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

π

 

1

k

2

sin

2

 

 

 

 

 

 

 

 

 

0

 

 

u

 

 

 

где обозначено k = sin(ϕ0 / 2) .

123

Данный интеграл не выражается через элементарные функции. Однако учитывая, что k sin x 1, подынтегральное выражение мож-

но разложить в ряд по формуле бинома Ньютона и после почленного интегрирования приходим к формуле

T =T

 

 

1

sin2

ϕ

 

 

1 3

2

sin4

ϕ

 

 

1

+

 

 

0

+

 

 

 

 

0

+... .

 

 

 

 

0

 

 

4

 

2

 

2 4

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С достаточной для практики точностью можно учесть только первые два члена этого ряда:

T =T

1

+

1

sin2

ϕ0

.

(3)

 

0

 

4

 

2

 

 

 

 

 

 

 

Значения относительного отклонения в процентах периода колебаний T, рассчитанного по формуле (3), от значения T0 для различных амплитуд следующие:

ϕ0, град

10

20

30

40

50

60

70

80

90

(TT0) / T0, %

0,2

0,7

1,7

2,9

4,5

6,3

8,2

10,3

12,5

Итак, мы выяснили, что колебания математического и физического маятника изохронны только при малых амплитудах. Точно так же колебания материальной точки, движущейся по дуге окружности, обращенной выпуклостью вниз, под действием силы тяжести будут в общем случае неизохронными. В связи с этим возникает вопрос: существует ли такая гладкая кривая, при движении по которой под действием силы тяжести материальная точка будет совершать гармонические изохронные колебания при любых амплитудах? Такая задача была решена Гюйгенсом, который показал, что такой кривой является перевернутая циклоида. При доказательстве Гюйгенс опирался в значительной степени на кинематические свойства циклоиды как кривой, которую описывает фиксированная точка окружности, катящейся без скольжения по прямой. Найдем период колебаний материальной точки, движущейся по циклоиде (циклоидальный маятник), исходя из более простых физических соображений.

124

Рассмотрим циклоиду, обращенную выпуклостью вниз. В соответствии с этим примем, что окружность расположена ниже горизонтальной прямой, по которой она катится (эта прямая на рис. 3.5 изображена штриховой линией). За ось X примем параллельную ей прямую, смещенную вниз на диаметр ок-

ружности 2R . Пусть точка A на катящейся окружности, описывающая циклоиду, в исходном положении находится на оси Y в наивысшей точке. Если окружность при качении повернется на угол ϕ, то ее центр C пере-

местится вправо на расстояние

Rϕ.

Рис. 3.5

 

При этом точка A сместится относительно центра влево на расстояние Rsin ϕ и вниз на расстояние R(1 cos ϕ) . Тогда прямоугольные коор-

динаты точки A примут значения:

x = R(ϕ−sin ϕ), y = R(1+ cos ϕ) .

Эти уравнения задают циклоиду в параметрическом виде (параметр ϕ является функцией времени ϕ = ωt , где ω равна угловой

скорости вращения колеса). Пусть теперь x и y означают координа-

ты материальной точки, совершающей циклоидальные колебания под действием силы тяжести. Тогда потенциальная энергия точ-

ки U = mgy , а кинетическая – E = m(x2 + y2 )/ 2 . Найдя производные x, y , после элементарных преобразований получаем

U = 2mgR cos2 ϕ/ 2, E = 2mR2ϕ2 sin2 ϕ/ 2.

Введем новую переменную q = cos ϕ/ 2 . Тогда q = −(1/ 2)ϕsin ϕ/ 2 . Эта переменная q может быть принята за обобщенную координату,

определяющую положение колеблющейся точки, а ее производная q – за обобщенную скорость. В этих обозначениях

U = 2mgRq2 , E =8mR2 q2 .

125

Данные соотношения показывают, что потенциальная энергия является квадратичной функцией координаты q , а кинетическая энергия – квадратичная функция скорости q с постоянными коэффициентами. Отсюда сразу делаем вывод о том, что при любых амплитудах колебания циклоидального маятника будут изохронными и гармоническими с периодом T = 2π 4R / g .

При малых амплитудах колебаний движение циклоидального маятника почти не отличается от движения математического маятника по окружности радиусом, равным радиусу кривизны циклоиды в нижней точке. Найдем его значение. Точка A (см. рис. 3.5) описывает данную циклоиду независимо от того, катится ли колесо равномерно или с ускорением, важно только, чтобы оно не проскальзывало. Проще всего рассмотреть, разумеется, равномерное качение колеса. Такое качение получается в результате сложения равномерного вращения колеса вокруг оси с угловой скоростью ω и равномерного поступательного движения, линейная скорость v которого равна ωR .

Во всех инерциальных системах отсчета материальная точка имеет одно и то же ускорение. Поэтому находить его можно в любой такой системе отсчета. Ясно, что ускорение точек обода колеса связано только с его вращением вокруг оси. Поэтому ускорение любой точки обода a направлено по радиусу к центру колеса и определяет-

ся выражением

 

a = v2 / R .

(4)

Значит, и в низшей точке циклоиды ускорение элемента обода колеса равно v2 / R и направлено вверх.

Рассмотрим теперь движение этой же точки обода как движение по циклоиде. Скорость в любой точке траектории направлена по касательной к ней; значит, в низшей точке циклоиды скорость направлена горизонтально. Ускорение же, как мы выяснили, направлено верти-

126

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

a =

V 2

,

(5)

r

 

 

 

где V – скорость точки обода в ее низшем положении; r – искомый радиус кривизны циклоиды.

Скорость любой точки обода катящегося колеса равна векторной сумме скорости поступательного движения колеса и линейной скорости вращения вокруг оси. При отсутствии проскальзывания эти скорости равны по модулю, а в нижней точке и направлены одинаково. Поэтому V = 2v и, сравнивая формулы (4) и (5), находим

r = 4R .

Так как радиус кривизны в нижней точке как раз и равен длине математического маятника l , то мы попутно получили и формулу для периода колебаний математического маятника при малых амплитудах T = 2π l / g .

127

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

3.2.1.Математический маятник между наклонными стен-

ками. Между двумя наклонными стенками с малым углом полураствора α находится шарик на нити длиной l (рис. 3.6). Толчком ему

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

Понятно, что если начальная скорость шарика настолько мала, что он не ударяется о стенки, то период колебаний рассчитывается как для обычного математического маятника:

Рис. 3.6

 

T =T0 = 2π

l

.

 

 

 

 

 

 

 

g

Найдем, при какой скорости шарика vmin это возможно. Исхо-

дя из закона сохранения энергии, получаем

 

 

 

1

mvmin

2

= mgl(1 cos α) .

2

 

 

 

 

 

 

При малом угле α находим vmin = α gl . Если же начальная скорость шарика больше vmin , то он начинает упруго ударяться о стенки, при-

чем после каждого удара скорость шарика просто изменяет знак. Это означает, что после удара участок 12 зависимости угла отклонения от времени (рис. 3.7) повторяет участок 34, который был бы при свободных колебаниях. Уравнение

 

свободных колебаний имеет вид

 

ϕ= ϕ0 sin ω0t ,

 

где ϕ0 – угловая амплитуда малых

Рис. 3.7

свободных колебаний, ϕ0 =v/ g / l;

128

ω0 – частота свободных колебаний,

ω0

= 2π/T0 . Полагая ϕ = α , най-

дем время достижения шариком стенки:

 

 

 

 

α = ϕ

 

sin ω t

t =

1

 

arcsin(α/ ϕ

 

) .

 

 

0

ω

0

 

 

 

 

0 1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

Период же колебаний будет больше этого времени в 4 раза:

 

T = 4t =

4

arcsin(α/ ϕ

) =T

2

arcsin

α gl

при v ≥ α

gl.

 

 

 

1

ω

 

0

0 π

 

 

 

v

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.2.2. Математический маятник на тележке. На гладкой го-

ризонтальной поверхности находится тележка массой M c установ-

ленным на ней математическим маятником длиной

l и массой m

(рис. 3.8). Найти период колебаний.

 

Эта задача существенно отличается

 

от обычного математического маятника

 

тем, что точка подвеса сама совершает

 

колебания в горизонтальной плоскости.

 

Попробуем несколько видоизменить

 

реальную ситуацию. Так как тележка вме-

 

сте со стойкой жесткая, то саму эту тележ-

 

ку можно представить как материальную

Рис. 3.8

точку массы M , находящуюся на проти-

 

воположном от шарика конце нити (см. рис. 3.8).

 

В этой системе двух материальных точек нет

 

внешних сил, действующих по горизонтали. По-

 

этому в силу закона сохранения импульса в ней

 

существует точка C (рис. 3.9), не совершающая

 

горизонтальных колебаний – центр инерции. Ко-

 

нечно, эта точка слегка колеблется по вертикали,

 

но в силу малости самих колебаний маятника

 

вертикальными колебаниями можно пренебречь.

 

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

 

материальной точки массой m относительно

Рис. 3.9

точки C на нити длиной l.

 

129

k(x2 x1 ) = m2

Положение центра инерции можно определить по формуле

x

=

mi xi .

c

 

mi

Помещая начало отсчета оси x в точке с массой m , легко находим l:

 

l

=

M

 

l .

 

 

 

 

 

m + M

 

 

 

 

Для периода колебаний получаем

 

 

 

 

T = 2π

 

l

= 2π

 

 

l

 

 

.

 

g

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

g 1

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

3.2.3. Два шарика на пружинке. На гладкой горизонтальной поверхности находятся два шарика массами m1 и m2 , связанные пружинкой жесткости k . Найти частоту колебаний такой системы

(рис. 3.10).

Рассмотрим несколько вариантов решения данной задачи, каждый из которых иллюстрирует

Рис. 3.10 различные физические подходы. Первый вариант. Так как мы наблюдаем одновременное движе-

ние двух тел, то разумно описать движение каждого из них. Пусть x1 и x2 – отклонения от положения равновесия первого и второго шарика. Запишем для них второй закон Ньютона:

k(x

x ) = m

d 2 x

,

1

dt2

2

1

1

(1)

d 2 x2

dt2

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

130

Соседние файлы в папке книги