- PVSM.RU - https://www.pvsm.ru -

Задача о жуке на глобусе

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

Задача о жуке на глобусе - 1

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

На шаре сидит жук массы m. В начальный момент времени система покоится, затем жук начинает ползти по шару так, что его траектория вычерчивает на шаре окружность радиуса b,quad ble R. Придя в исходную точку, жук останавливается.

На какой угол повернулся шар к моменту остановки жука, если сила тяжести отсутствует?

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

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

Обозначим через T время от начала движения жука до его остановки.

Свяжем с шаром декартову систему координат Ox'y'z' так, чтобы жук двигался в плоскости z'=h,quad (h^2+b^2=R^2). Тогда радиус-вектор жука выражается формулой

boldsymbol r(t)=hboldsymbol e'_z+b(cosalphaboldsymbol e'_x+sinalphaboldsymbol e'_y).

Функция alpha=alpha(t) по условию монотонно возрастает на отрезке [0,T] и такова, что

alpha(0)=0,quaddotalpha(0)=dotalpha(T)=0,quad alpha(T)=2pi.

Кинетический момент системы относительно точки O, очевидно, сохраняется:

boldsymbol K_O=Jboldsymbolomega+m[boldsymbol r,boldsymbol v]=0.qquad (1)

Здесь boldsymbolomega — угловая скорость шара относительно лабораторной системы координат, а boldsymbol v – скорость жука.

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

boldsymbol v=boldsymbol v_r+boldsymbol v_e,

где

boldsymbol v_e=[boldsymbolomega, boldsymbol r]

— переносная скорость жука, а

boldsymbol v_r=bdotalpha(-sinalphaboldsymbol e'_x+cosalphaboldsymbol e'_y)

— его относительная скорость.

Подставляя эти формулы в уравнение (1) и решая получившееся линейное алгебраическое уравнение относительно вектора boldsymbolomega, находим:

boldsymbol omega=-frac{m}{J+mR^2}[boldsymbol r,boldsymbol v_r],qquad(2)

где

[boldsymbol r,boldsymbol v_r]=dotalphabig(b^2boldsymbol e'_z-hbcosalphaboldsymbol e'_x-hbsinalphaboldsymbol e'_ybig).qquad(3)

Формула (2) отражает важный факт: в тот момент, когда жук останавливается на шаре, шар тоже прекращает вращение.

Другой важный факт, следующий из формул (2) и (3), заключается в следующем: вектор угловой скорости шара, заданный в связанной с шаром системе координат, представляется в виде

boldsymbolomega=dotalphaboldsymbol w(alpha),

где

boldsymbol w=-frac{mbig(b^2boldsymbol e'_z-hbcosalphaboldsymbol e'_x-hbsinalphaboldsymbol e'_ybig)}{J+mR^2}.qquad (4)

И тут начинается самое интересное. Введем лабораторную систему координат Oxyz так, чтобы при t=0 она совпадала с системой Ox'y'z'.

Через X(t) обозначим матрицу размера 3 times 3, которая является решением следующей задачи Коши для матричного дифференциального уравнения:

dot X=dotalpha (t) XW(alpha),quad X(0)=I,qquad(5)

где

W=begin{pmatrix}0 & -w_z & w_y \w_z & 0 & -w_x \-w_y & w_x &0end{pmatrix},

а w_x, w_y, w_z — компоненты вектора boldsymbol w (см. формулу (4)) в системе Ox'y'z'.

Смысл оператора X(t):mathbb{R}^3to mathbb{R}^3 следующий. Предположим, вектор boldsymbol u жестко связан с шаром. Его координаты относительно системы Ox'y'z' не меняются, а относительно системы Oxyz изменяются, поскольку этот вектор поворачивается вместе с шаром. То есть наблюдатель, находящийся в системе Oxyz, видит зависимость координат вектора boldsymbol u от времени: boldsymbol u=boldsymbol u(t).

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

boldsymbol u(t)=X(t)boldsymbol u(0).

Таким образом, X(T) — это матрица поворота шара. В этой матрице и содержится ответ на вопрос задачи.

Сделаем замену переменной tmapsto alpha в (5):

frac{d}{dalpha} tilde X=tilde XW(alpha),quad tilde X(0)=I.qquad(6)

Решением этой системы является матрица tilde X(alpha); соответственно, искомая матрица поворота шара — этоtilde X(2pi),quad tilde X(2pi)=X(T) .

Важный вывод из этого наблюдения состоит в том, что матрица поворота шара совершенно не зависит от конкретного закона движения жука alpha=alpha(t).

Маловероятно, что систему (6) можно проинтегрировать в замкнутой форме. Этот вопрос требует отдельного исследования. Однако матрицу поворота шара всегда можно найти численно, решив задачу (6) на интервале alphain[0,2pi].

Автор: drzewo

Источник [1]


Сайт-источник PVSM.RU: https://www.pvsm.ru

Путь до страницы источника: https://www.pvsm.ru/teoreticheskaya-mehanika/451734

Ссылки в тексте:

[1] Источник: https://habr.com/ru/articles/1035972/?utm_campaign=1035972&utm_source=habrahabr&utm_medium=rss