Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Основы математических знаний для изучения физики..pdf
Скачиваний:
21
Добавлен:
15.11.2022
Размер:
4 Mб
Скачать

п л2~ п

(3.25)

/=1 аг /=1

Как видно, в результирующий вектор сил F входит лишь сумма внешних сил. Сумма всех внутренних сил оказывается равной нулю в силу третьего закона Ньютона: f jk = -/* , . Записывая левую часть уравнения (3.25) с учетом

определения центра масс, получаем искомое уравнение

d2rc _

d^.

F

(3.26)

т ^ 2

 

dt

Отсюда следует, что при отсутствии внешних сил или при равенстве нулю результирующей внешних сил скорость центра масс постоянна: vc = const.

3.3. Колебания пружинного маятника

Различного рода колебания широко представлены в нашем мире. Достаточно вспомнить, что без звуковых колебаний мы остались бы глухими, а без электромагнитных - слепыми. Широко распространены колебания и в механике. Неслучайно поэтому физики, механики и математики с давних пор уделяли изучению колебаний много внимания. Книг по теории колебаний множество. До сих пор пользуется популярностью книга академика Л.И. Мандельштама [8], который является одним из создателей нелинейной теории колебаний.

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

Уравнение движения (второй закон Ньютона) имеет вид

тх = -кх.

(3.27)

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

Рис. 12. Пружинный маятник Ось х соответствует направлению колебаний, начало координат выбрано в

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

Уравнение (3.27) обычно записывают в виде

* + <nJjr = 0, cogs—.

(3.28)

т

 

Общее решение этого уравнения может быть записано в двух

эквивалентных формах

 

x(t) = с, C O S ( CD/ ) + с2sin(co0/) = A cos(co0/ + ф0).

(3.29)

Вычислите вторую производную от (3.29), подставьте ее в (3.28) и убедитесь, что (3.29) есть действительно решение уравнения (3.27) при любых начальных условиях. Между параметрами первой и второй форм существует взаимно однозначное соответствие

А =ф1 + с1, tg(<p0) = -^ -.

(3.30)

ci

 

Вторая форма записи удобна тем, что в ней сразу видна амплитуда колебаний. Так как аргумент во второй форме определяет фазу колебания, величину ср0 называют начальной фазой. Первая форма удобнее для

удовлетворения начальных условий.

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

*(0) = лг0, х(0) = хо,

(3.31)

где х0, х0 - заданные значения соответственно координаты и скорости в

начальный момент времени.

Согласно первой форме решения в начальный момент времени должно быть выполнено соотношение

*(0) = с, cos(0) + с2sin(0) = с} =х0.

(3.32)

Взяв производную от функции х(/) и положив / = 0, получим второе соотношение

С учетом (3.32), (3.33) конкретное решение может быть записано в виде

x(t) = х0cos(co0 t) + — sin(co0 t) .

(3.34)

«о

 

Видно, что, если в начальный момент времени

скорость равна нулю,

врешении останется только первое слагаемое. И наоборот, если отклонение

вначальный момент времени равно нулю, в уравнении останется только второе слагаемое. В общем же случае присутствуют оба слагаемых и в соответствии

с формулой (3.30) амплитуда колебаний А = yjxl + (x0/ со0)2

Для полного понимания свойств полученного решения полезно вычислить потенциальную и кинетическую энергию колеблющегося тела:

£ ,(0 = f * 2(0. Ek= ^ x \ t ) ,

(3.35)

и убедиться в том, что сумма потенциальной и кинетической энергии от времени не зависит:

Е = Ер +Ек = ^ х 20+ j x t = const.

(3.36)

Опыт говорит о том, что любые колебания без подпитки энергии затухают. Сохранение энергии при колебаниях говорит о том, что в использованной модели не учтены диссипативные процессы. Ближе к реальной ситуации модель с учетом силы вязкого сопротивления пропорциональной скорости. Учет этой силы приводит к необходимости решать следующее дифференциальное уравнение:

тх - -кх - Ъх,

(3.37)

где Ь - коэффициент сопротивления >0).

Поделив все слагаемые на массу, запишем это уравнение в «стандартном»

виде:

 

х + 2р* + ©о* = 0, р = 0,5Ытг O)Q= к/т

(3.38)

Коэффициент Р называют коэффициентом затухания. С точки зрения

математики (3.38) - дифференциальное уравнение с постоянными коэф­ фициентами. Со времен Эйлера на первом этапе нахождения решения таких уравнений отыскивают корни так называемого характеристического уравнения. Характеристическое уравнение получается при подстановке в дифференциальное уравнение замены х = ехр(^0- В нашем случае эта замена дает (после сокращения на ехр(7.0) квадратное уравнения для X:

X,2 + 2{ЗХ.+ C0Q= 0.

(3.39)

Корни квадратного уравнения таковы:

V 2 = -P ± VPJ-0)0 •

(3.40)

Видно, что корни характеристического уравнения различны. Из теории известно, что в этом случае решение дифференциального уравнения надо разыскивать в виде линейной комбинации двух решений

х = а1• exp (A.,f) + а2• exp (k2t) .

(3.41)

Так мы и поступим. Но вначале нужно разобраться - какие корни уравнения являются вещественными, а какие комплексными. Допустим вначале, что корни вещественные. Тогда должно выполняться неравенство

р2>Шо, Р>®0-

(3.42)

Это неравенство предполагает, что коэффициент трения велик. Из формул (3.36) следует, что при выполнении (3.42) оба корня Х.р Х2 отрицательны и, следовательно, согласно (3.41) решение является экспоненциально затухающим. Этот случай соответствует так называемым

ангармоническим колебаниям.

Наиболее интересен случай обратного неравенства

Р2<С02, Р<ю0.

(3.43)

Это неравенство предполагает, что коэффициент трения мал. При выполнении этого неравенства возможна замена:

VP2- CO2=V-(®O"P 2) = ф э о - Р 2) >

с помощью введения мнимой единицы / = V—Т - Таким образом, при малой величине трения корни характеристического уравнения комплексные:

Решение (3.41) для этих корней может быть записано в виде

x = ar e(_p+'“)' + аг • е(' мь>' = e’13' (a, • e+'“' + a2 • e ^ ').

(3.45)

А теперь можно воспользоваться знаменитой формулой Эйлера, связывающей экспоненциальную зависимость с тригонометрическими функциями:

е/г = cosг +/sinz .

Но можно обойтись и без формулы Эйлера, заметив, что в (3.45) выделена экспоненциальная зависимость от времени. Это обстоятельство подсказывает искать решение в виде произведения двух функций:

*(0 = w(/)-e~p'

(3.46)

Подстановка (3.46) в дифференциальное уравнение (3.38) дает знакомое нам дифференциальное уравнение для функции w:

W + (CDQ- p2)w= w +a)2w= 0, (b = yj(o20 - p 2

(3.47)

Сравнение этого уравнения с уравнением для незатухающих колебаний (3.28) позволяет утверждать, что функция w(f) будет иметь решение вида (3.29) . Следовательно, можно положить, что

w(t) = с, cos(a>0+ с2sin(Ш) = В0cos(cbt + ф0).

(3.48)

Окончательное решение в соответствии с (3.46) может быть записано в

виде

 

х(() = В0• е"р/ • cos(cof + q>0) .

(3.49)

Из этого решения видны два отличия по сравнению со случаем без учета трения. Первое отличие заключается в том, что амплитуда колебаний уменьшается со временем по экспоненте. Второе отличие заключается в том,

что колебания происходят с уменьшенной циклической частотой со = yj(o20- $ 2

Оба эти отличия понятны с физической точки зрения - вязкое трение всегда уменьшает энергию системы и уменьшает скорость движения. Но точные зависимости этих эффектов получены математикой.

Обсудим характеристики затухания. Время, за которое начальная

амплитуда уменьшается в е раз (е

- основание натуральных логарифмов)

называется временем релаксации. Время релаксации т = 1/р.

За время релаксации совершится число колебаний

т _

&

_ \j(£>о-ft2

f

2тср

2лр

Если представить собственную частоту в виде со0 = ф (напомним, что р в рассматриваемом случае больше единицы), формула (3.50) может быть записана в виде п = у[р2-1 /(2п). Например, при р =100, и « 16.

Характеристиками затухания являются также логарифмический коэффициент затухания

8= 1п(- A(t)

n

>

A(t + f )

 

и добротность колебательной системы

 

 

E{t)

Q =2n

E (t)-E {t +f )

(3.51)

(3.52)

Из определения добротности следует, что это есть произведение 2л на отношение энергии в произвольный момент времени к убыли этой энергии за промежуток времени от t до (t +f ). Так как энергия пропорциональна амплитуде колебаний А(г), для добротности может быть получена формула

A \t)

(3.53)

б = 2л

l - e - 2pf

й Т 5-'

A2(t) - A 2(t +f )

 

При малых значениях логарифмического коэффициента затухания (5<с1)знаменатель в формуле (3.53) «26 (см. раздел 1.5) и добротность может

быть оценена как Q ^ n /b =nn .

3.4. Вращательное движение

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

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

О
Рис. 13. Вращение тела вокруг вертикальной оси

ускорение. Для описания вращательного движения используют вектор элементарного поворота dcp. По модулю этот вектор равен углу поворота тела вокруг оси за интервал времени dt, а направлен он вдоль оси вращения. В правой системе координат его однозначное направление вдоль оси вращения определяется по правилу правого винта (из конца вектора бф поворот тела происходит против часовой стрелки).

Кинематической характеристикой вращения является угловая скорость

_ dw J .

(3.54)

со = ——=ф

dt

Вращение называется равномерным, если угловая скорость постоянна: со = const. В этом случае общим решением уравнения (3.54) является линейная зависимость угла поворота от времени:

ср(/) = со/ + ср(0).

(3.55)

При вращении линейная скорость произвольной точки тела А (рис. 13)

v=-^- = [fflxF] = [raxp].

(3.56)

Радиус-вектор произвольной точки г отсчитывается от произ­ вольной точки на оси вращения, принятой за начало координат, а век­ тор р проведен от оси вращения по кратчайшему расстоянию к рассмат­ риваемой точке тела (см. рис. 13).

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

Время одного полного оборота тела при равномерном вращении дает значение периода вращения Т

Период вращения обратно пропорционален частоте вращения п\

т_ \ _ 2 к

пп

Угловое ускорение определяется как производная от угловой скорости:

_ _ dco __d2cp

(3.58)

Зная угловую скорость и угловое ускорение, можно вычислить ускорение произвольной точки N вращающегося тела (рис. 14):

а =— = —[©xr] = [sxp] + [ajxv].

(3.59)

/ '

 

 

 

I\

 

С учетом (3.56) второе слагаемое в (3.59)

V

 

Чччв

может быть представлено иначе:

 

 

 

а = [coxp] + [^cox[coxp]J = [сохр]-со2 - р . (3.60)

 

 

Первое слагаемое в (3.60) называют

 

Рис. 14. Угловые скорость

тангенциальным (касательным) ускорением

 

и ускорение

ах, а второе нормальным ускорением ап,

 

 

a, = [ i x p] = [sxF],

ап= - а 2 ■р

(3.61)

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

Пока описана лишь кинематика вращательного движения. Для описания динамики вращательного движения необходимо ввести еще два понятия - момента силы и момента импульса.

Моментом силы относительно какой-либо точки О называют векторную величину

М = х / J, A/ = /r -r-sin(a) = /r -/.

(3.62)

Величину / называют плечом силы.

Момент импульса материальной точки относительно фиксированной

точки О вычисляется с помощью векторного произведения:

 

Li=mi[rl *vl\ = [ri x p ] .

(3.63)

Соответственно, моментом импульса механической системы, состоящей из п материальных точек, является сумма

 

£ = Ё А =!Ё [ ^ ха ]-

 

(3.64)

Продифференцировав это выражение по времени, получим

 

dL Л d г—. _ л

г— - I

- dPi

 

¥ = § ^ г' * й ] - 5

г ><й ,+ Гй* dt

=Z к1x dt

(3.65)

При вычислении производной от момента импульса учтено, что векторное произведение скорости на импульс равно нулю в силу параллельности соответствующих векторов. Для получения связи момента импульса с моментом силы следует воспользоваться вторым законом Ньютона, дающим связь производной от импульса по времени с силой:

Ф

= у м 1= м

(3.66)

dt

Фактически уравнение динамики вращательного движения уже получено. Оно утверждает, что производная от момента импульса системы определяется суммарным моментом сил

= М

(3.67)

dt

 

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

Lx =lx<ox, Ly = ly% , Lz =/2со2.

(3.68)

Для вычисления момента инерции тела относительно выбранной оси рассмотрим вначале, как обычно, момент инерции материальной точки. Для упрощения вывода будем предполагать, что вращение происходит вокруг одной фиксированной оси. При этом нам можно опустить индекс, соответствующий выбранной оси. Тогда момент инерции материальной точки с номером / и массой mi вычисляется по формуле (рис. 15)

L = m;r 2

Момент инерции системы мате­ риальных точек выражается соответст­ вующей суммой моментов инерции, а в случае сплошного тела - соответствующим интегралом по объему тела:

1= jpdVr2 = jpr2dV,

(3.70)

V V

(3.69)

пnti

W

i О

Рис. 15. Параметры момента инерции материальной точки

где р - плотность; г - расстояние до оси вращения.

Сложность вычисления интеграла (3.70) зависит от его формы. В случае простейших форм типа стержня, цилиндра, шара, куба вычисление интеграла не вызывает затруднений. Кроме того, вычисление интегралов облегчается при использовании теоремы Штейнера (Гюйгенса-Штейнера). Эта теорема позволяет по значению момента инерции относительно центра масс /с вычислить момент инерции относительно любой оси по формуле

I =Jc + m d \

(3.71)

где d - расстояние выбранной оси до центра масс.

Теорема существенно упрощает вычисления моментов инерции, так как требует лишь знания момента инерции тела относительно центра масс /с.

Момент инерции тела относительно центра масс /с особенно легко вычисляется для тел с осевой симметрией. Примеры вычисления моментов инерции рассмотрим позднее, а пока вернемся к закону вращательного движения. С учетом формул (3.67), (3.68) имеем для описания вращения около одной оси уравнение

/ — = / ^ = Л/

(3.72)

d/ dr

Как видно, закон вращательного движения с математической точки зрения эквивалентен закону движения материальной точки. При этом полезно знать «соответствие» величин поступательного и вращательного движения. Укажем это соответствие, сравнивая одномерное движение по оси х с вра­ щением около одной оси, описываемым углом ср (таблица).

Соответствие характеристик движения - поступательного и вращательного

Характеристика

Поступательное

Вращательное

 

движение

движение

Инерционность

Масса т

Момент инерции I

Координата

X

 

Угол ср

Скорость

V -X

Угловая скорость

 

 

 

со = ф

Ускорение

a = v =x

Угловое ускорение

 

 

 

8 = СО= ф

Сила

F

Момент силы М

 

~ В

силу

соответствия

 

математических

уравнений

для

 

вращения и поступательного движения

 

решение

уравнения (3.72)

находится

 

точно таким же способом, как для

 

поступательного движения.

 

 

 

Приведем

пример вычисления

 

момента инерции для цилиндрического

 

тела (рис.

16) с тремя параметрами -

 

Я (высота

цилиндра),

R{ (внутренний

 

радиус цилиндра),

(внешний радиус

 

цилиндра).

Плотность

цилиндра

р

Рис. 16. Параметры инерции

полагаем постоянной. Момент инерции

цилиндрического тела

зависит

от

выбранной

оси

вращения.

 

Будем

считать,

что

ось

вращения

соответствует оси цилиндра. Из соображений симметрии ясно, что ось вращения проходит через центр масс. Выделим в цилиндре бесконечно тонкий цилиндрический слой толщиной d r. В соответствии с формулой (3.69) момент инерции этого слоя

d/ = dm г2= (р • Я • 2пг • dr)r2

Полный момент инерции равен интегралу:

Rl \ D4 _ R4

/ = jd/ = р-2л*Я jV3dr = р - п Н —

-----L

я,

я,

2

Эту формулу можно упростить, используя значение полной массы тела