Решение задачи по механике: вариативность методов и теоретическое обобщение

                       Классическая задача теоретической механики

Проанализируем известную задачу из сборника И. В. Мещерского (задача № 31.22), описывающую динамику материальной точки на криволинейной поверхности.

Условие: Малое тело массой m помещено в высшую точку идеально гладкой полусферы радиуса R. В условиях земного тяготения (ускорение g) объекту сообщают начальный горизонтальный импульс со скоростью v0. Необходимо вычислить критический угол φ, при котором происходит разрыв кинематической связи и тело покидает поверхность сферы.

Геометрическими параметрами тела при расчетах пренебрегаем.

Решение: Представим графическую интерпретацию физической модели.

 Рисунок 1 Условие задачи
 Рисунок 1. Иллюстрация к постановке задачи

  Рассмотрим несколько аналитических подходов к решению:

1. Базовый метод (уровень углубленной школьной программы): Сформулируем уравнение второго закона Ньютона в проекции на нормальную ось, направленную к центру кривизны O:

ma_ц = mg\cos\phi - N

, где

a_ц = \frac{v^2}{R}

— центростремительное ускорение, v — мгновенная скорость, N — сила нормальной реакции. В точке отрыва взаимодействие с поверхностью прекращается (N = 0). После сокращения на массу m получаем:

 \frac{v^2}{R}=gcos\phi

Воспользуемся законом сохранения механической энергии:

\frac{mv_0^2}{2} + mgR = \frac{mv^2}{2} + mgh

что эквивалентно:

\frac{mv_0^2}{2} + mgR = \frac{mv^2}{2} + mgR(1-cos\phi)

Интегрируя полученные выражения, находим:

{v_0^2} + 2gR = 3gRcos\phi

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

\phi = \arccos\left(\frac{v_0^2+2gR}{3gR}\right)

Заметим, что если начальная скорость

v = \frac{ds}{dt} = R\frac{d\phi}{dt}

откуда

\frac{d\phi}{dt} = \frac{v}{R}

Запишем уравнение движения в проекции на касательный вектор (тангенциальное ускорение):

\frac{dv}{dt} = g\sin\phi

Исключая время через деление уравнений и разделяя переменные, приходим к виду:

vdv = gR\sin\phi d\phi

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

3. Обобщение для произвольного выпуклого профиля: Рассмотрим движение тела по любой гладкой кривой y(x), обладающей отрицательной кривизной (y’’(x) < 0).

Рисунок 2 Обобщение задачи на произвольную выпуклую вниз поверхность.
Рисунок 2. Расширение задачи на случай произвольной выпуклой кривой.

Введем декартову систему координат с началом в вершине (0, H), где y(0) = H и y’(0) = 0 (условие экстремума). В произвольной точке отрыва (x, y) уравнение динамики принимает вид:

\frac{v^2}{R} = gcos\phi

Здесь R — радиус кривизны траектории, а φ — угол между вертикалью и нормалью.

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

{v^2} = gh=g(y-yO’)

, где h — проекция расстояния от мгновенного центра вращения O’ до точки нахождения тела на вертикальную ось.

Согласно энергетическому балансу:

\frac{mv_0^2}{2} + mgH = \frac{mv^2}{2} + mgy

Объединяя выражения, получаем дифференциальное уравнение:

\left(v_0^2 + 2g(H-y)\right)y'' + g(1 + (y')^2) = 0

Решение задачи Коши при заданных граничных условиях приводит к параболической траектории:

y=H-\frac{gx^2}{2v_0^2}

Точка отрыва соответствует пересечению заданной кривой y(x) с этой параболой. Если пересечение происходит только в вершине, отрыв является немедленным.

4.1. Анализ системы с диссипативными силами: Усложним модель, добавив силу трения скольжения, определяемую коэффициентом μ (Fтр = μN).

Рисунок 3 Добавление силы трения скольжения
Рисунок 3. Моделирование движения с учетом трения.

Спроецируем силы на оси динамического базиса:

ma_ц = mg\cos\phi - N m\frac{dv}{dt} = -mgsin\phi + \mu N

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

2v\frac{dv}{d\phi} = g(\mu\cos\phi - \sin\phi) - 2\mu v^2

Используя метод Бернулли для решения при начальном условии v(0) = v0, находим функциональную зависимость скорости от угла:

v^2(\phi)=\frac{2gR}{\mu^2+1}((2\mu^2+1)cos\phi-\mu\sin\phi)+(v_0^2-2gR\frac{1+2\mu^2}{1+4\mu^2})e^{-2\mu\phi}

Условие отрыва (N = 0) требует выполнения равенства

v^2 = gR \cos\varphi

. Подставляя это в решение, получаем трансцендентное уравнение:

gRcos\phi=\frac{2gR}{\mu^2+1}((2\mu^2+1)cos\phi-\mu\sin\phi)+(v_0^2-2gR\frac{1+2\mu^2}{1+4\mu^2})e^{-2\mu\phi}

Численное решение этого уравнения в интервале (0; π/2) позволяет найти искомый угол отрыва.

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

F_c=kv^2

, не меняет принципиальную разрешимость модели. Полученное уравнение:

\frac{dv}{d\phi} = g(\mu\cos\phi - \sin\phi) - \frac{1}{2}\left(2\mu v^2 + kv^2\right)

сохраняет линейность относительно v² и решается аналогичными методами.

                                                Библиографический список:

1. Мещерский И. В. Сборник задач по теоретической механике, 1980.

2. Батяев Е. А. Теоретическая механика.

3. Яблонский А. А. Курс теоретической механики.

4. Тарг С. М. Краткий курс теоретической механики.

 

Источник

Читайте также