Известия РАН. Механика твердого тела, 2019, № 3, стр. 69-75
О ПЕРЕХОДНЫХ ПРОЦЕССАХ В ДИНАМИКЕ ЭЛЛИПСОИДА ВРАЩЕНИЯ НА ПЛОСКОСТИ C ТРЕНИЕМ
М. А. Муницына *
Московский физико-технический институт (государственный университет)
Долгопрудный, Россия
* E-mail: munitsyna@gmail.com
Поступила в редакцию 31.07.2017
После доработки 31.07.2017
Принята к публикации 02.11.2017
Аннотация
Рассматривается задача о движении эллипсоида вращения, центр масс которого совпадает с геометрическим центром, а ось вращения является осью динамической симметрии, на горизонтальной плоскости с трением. При некотором классе начальных условий дается качественный анализ динамики. Приводится геометрическая интерпретация результатов и результаты численных экспериментов.
1. Постановка задачи. Рассмотрим задачу о движении эллипсоида вращения по горизонтальной плоскости с трением. Будем считать, что центр масс эллипсоида совпадает с его геометрическим центром, а ось вращения является осью динамической симметрии. Введем подвижную систему координат $C{{{\xi }}_{{\text{1}}}}{{{\xi }}_{{\text{2}}}}{{{\xi }}_{3}}$, оси которой совпадают с главными центральными осями инерции эллипсоида. Пусть $a$ и $c$ – экваториальный и осевой радиусы эллипсоида соответственно, и уравнение его поверхности имеет вид $f({\xi }) = ({\xi }_{1}^{2} + {\xi }_{2}^{2}){\text{/}}{{a}^{2}} + {\xi }_{3}^{2}{\text{/}}{{c}^{2}}$ = 0. Поскольку вектор восходящей вертикали γ определяется равенством
(1.1)
${\mathbf{\gamma }} = - \frac{{{\text{grad}}f({\mathbf{r}})}}{{\left\| {{\text{grad}}f({\mathbf{r}})} \right\|}}$(1.2)
${{r}_{1}} = - \frac{{{{a}^{2}}{{{\gamma }}_{1}}}}{{\rho }},\quad {{r}_{2}} = - \frac{{{{a}^{2}}{{{\gamma }}_{2}}}}{{\rho }},\quad {{r}_{3}} = - \frac{{{{c}^{2}}{{{\gamma }}_{3}}}}{{\rho }}$Теоремы об изменении количества движения и кинетического момента эллипсоида, условие постоянства вектора ${\kern 1pt} {\gamma }$ и условие безотрывности движения имеют соответственно вид
(1.3)
$m\mathop {\mathbf{v}}\limits^ \bullet \; + [{\mathbf{\omega }},m{\mathbf{v}}] = - mg{\mathbf{\gamma }} + N{\mathbf{\gamma }} + {\mathbf{F}}$(1.4)
$J\mathop {\mathbf{\omega }}\limits^ \bullet + [{\mathbf{\omega }},J{\mathbf{\omega }}] = [{\mathbf{r}},N{\mathbf{\gamma }} + {\mathbf{F}}] + {\mathbf{M}}$Здесь m – масса эллипсоида, ${\kern 1pt} {\mathbf{v}}{\kern 1pt} $ и ${\kern 1pt} {\mathbf{\omega }}{\kern 1pt} $ – векторы скорости центра масс и угловой скорости соответственно, ${\mathbf{J}} = {\text{diag}}(A,A,C)$ – центральный тензор инерции эллипсоида, $g$ – ускорение свободного падения, $N \geqslant 0$ – величина нормальной составляющей реакции опорной плоскости, ${\kern 1pt} {\mathbf{F}}{\kern 1pt} $ и ${\kern 1pt} {\mathbf{M}}{\kern 1pt} $ – сила и момент трения соответственно. При заданной модели трения
(1.7)
${\mathbf{F}} = {\mathbf{F}}({\mathbf{v,\omega ,\gamma ,}}N),\quad {\mathbf{M}} = {\mathbf{M}}({\mathbf{v,\omega ,\gamma }},N)$2. Абсолютно гладкая плоскость. Приведем результаты исследований динамики тела на абсолютно гладкой плоскости [1–3] в случае эллипсоида вращения.
При F = 0, M = 0 уравнения (1.3)–(1.6) допускают интеграл энергии $H = h = {\text{const}}$,
(2.1)
$H = \frac{m}{2}{{{\mathbf{v}}}^{2}} + \frac{1}{2}(J{\mathbf{\omega ,\omega }}) - mg({\mathbf{r,\gamma }})$(2.2)
$K = A({{{\omega }}_{1}}{{{\gamma }}_{1}} + {{{\omega }}_{2}}{{{\gamma }}_{2}}) + C{{{\omega }}_{3}}{{{\gamma }}_{3}}$При $k = \pm C{\omega }$ система (1.3)–(1.6) имеет решения вида
(2.4)
$\begin{array}{*{20}{c}} {{{v}_{1}} = ucos{\omega }t,\quad {{v}_{2}} = - usin{\omega }t,\quad {{v}_{3}} = 0,\quad \quad {{{\omega }}_{1}} = {{{\omega }}_{2}} = 0} \\ {{{{\omega }}_{3}} = {\omega },\quad \quad {{{\gamma }}_{1}} = 0,\quad {{{\gamma }}_{2}} = 0,\quad {{{\gamma }}_{3}} = \pm 1,\quad N = mg} \end{array}$При $k \ne \pm C{\omega }$ (${{{\gamma }}_{3}} \ne \pm 1$) указанные выше интегралы позволяют провести редукцию к одномерному движению и представить (2.1) в виде
(2.5)
$\frac{{\alpha }}{2}{{{\dot {\theta }}}^{2}} + V({\theta }) = h - \frac{1}{2}(C{{{\omega }}^{2}} + m{{u}^{2}})$Движения эллипсоида с постоянным углом нутации ${\theta }$ соответствуют критическим точкам функции $V({\theta })$. В частности, при $k = 0$ и
в точках ${\theta } = \pm {\pi /}2$ достигается минимум $W({\theta })$, и соответствующие решения системы (1.3)–(1.6)(2.6)
$\begin{array}{*{20}{c}} {N = mg,\quad {{v}_{1}} = usin{\varphi }cos{\omega }t,\quad {{v}_{2}} = - usin{\varphi }sin{\omega }t,\quad {{v}_{3}} = ucos{\varphi }} \\ {{{{\omega }}_{1}} = 0,\quad {{{\omega }}_{2}} = 0,\quad {{{\omega }}_{3}} = {\omega },\quad {{{\gamma }}_{1}} = sin{\omega }t,\quad {{{\gamma }}_{2}} = cos{\omega }t,\quad {{{\gamma }}_{3}} = 0} \end{array}$Кроме того, при
(2.7)
$\begin{array}{*{20}{c}} {{{v}_{1}} = {{u}_{1}}(cos{\psi }cos{\varphi } - sin{\psi }cos{{{\theta }}_{0}}sin{\varphi }) + {{u}_{2}}(sin{\psi }cos{\varphi } + cos{\psi }cos{{{\theta }}_{0}}sin{\varphi }) + O({\varepsilon })} \\ {{{v}_{2}} = {{u}_{1}}( - cos{\psi }sin{\varphi } - sin{\psi }cos{{{\theta }}_{0}}cos{\varphi }) + {{u}_{2}}( - sin{\psi }sin{\varphi } + cos{\psi }cos{{{\theta }}_{0}}cos{\varphi }) + O({\varepsilon })} \\ {{{v}_{3}} = {{u}_{1}}sin{\psi }sin{{{\theta }}_{0}} - {{u}_{2}}cos{\psi }sin{{{\theta }}_{0}} + O({\varepsilon })} \\ {{{{\omega }}_{1}} = O({\varepsilon }),\quad {{{\omega }}_{2}} = O({\varepsilon }),\quad {{{\omega }}_{3}} = {\omega }} \\ {{{{\gamma }}_{1}} = cos{\varphi }sin{{{\theta }}_{0}} + O({\varepsilon }),\quad {{{\gamma }}_{2}} = sin{\varphi }sin{{{\theta }}_{0}} + O({\varepsilon })} \\ {{{{\gamma }}_{3}} = cos{{{\theta }}_{0}} + O({\varepsilon }),\quad N = mg + O({\varepsilon })} \end{array}$3. Плоскость с трением скольжения. Рассмотрим случай плоскости с трением скольжения (${\mathbf{M}} = 0$). Будем считать, что сила трения обращается в нуль только одновременно со скоростью проскальзывания. Тогда из рассмотренных выше движений реализуются те, на которых скорость точки контакта эллипсоида с опорной плоскостью равна нулю, т.е. решения вида (2.4) при $u = 0$ и решения (2.6) при $u = a{\omega }$.
Для исследования устойчивости указанных решений необходимо конкретизировать модель трения. Например, при классическом вязком трении ${\mathbf{F}} = - {\kappa }({\mathbf{v}} + [{\mathbf{\omega ,r}}])$ этот вопрос исследован в [4] для однородного эллипсоида. При произвольных моментах инерции A и C равномерные вращения (2.4) с неподвижным центром масс ($u = 0$) вытянутого вдоль оси симметрии эллипсоида (c > a) устойчивы при ${{{\omega }}^{2}} > {\omega }_{I}^{2}$, а для сжатого – при выполнении обратного неравенства. Равномерные качения (2.6) без проскальзывания ($u = a{\omega }$) вытянутого (сжатого) эллипсоида, параметры которого удовлетворяют неравенству $(c - a)(C - A) < 0$, устойчивы при ${{{\omega }}^{2}} > {\omega }_{{II}}^{2}$ (${{{\omega }}^{2}} < {\omega }_{{II}}^{2}$). Здесь
Из этих результатов в частности следует, что если вытянутый вдоль эллипсоид расположить на плоскости так, что ось его симметрии горизонтальна, и сообщить центру масс достаточно большую горизонтальную скорость, то в силу неустойчивости соответствующих качений ось симметрии не будет оставаться в горизонтальном положении. При этом она может занять вертикальное положение, соответствующее устойчивым решениям (2.4).
Для исследования указанных переходных процессов заметим, что проекция уравнения (1.4) на ось динамической симметрии эллипсоида имеет вид
(3.1)
$C{{{\dot {\omega }}}_{3}} = ({{f}_{1}}{{{\gamma }}_{2}} - {{f}_{2}}{{{\gamma }}_{1}})\frac{{{{a}^{2}}}}{{\rho }}$(3.2)
$\dot {K} = - \frac{{{{c}^{2}} - {{a}^{2}}}}{{\rho }}({{f}_{1}}{{{\gamma }}_{2}} - {{f}_{2}}{{{\gamma }}_{1}}){{{\gamma }}_{3}}$(3.3)
$\dot {U} = \frac{{([{\mathbf{v}},{\mathbf{\gamma }}],[{\mathbf{F}},{\mathbf{\gamma }}])}}{{mU}}$Выберем в качестве указанного периодического решения (2.7), откинув слагаемые малого порядка. Тогда при классическом вязком трении из (3.1), (3.2), (3.3) получим
(3.4)
$\begin{array}{*{20}{c}} {{{{{\dot {\omega }}}}_{3}} = - {\kappa }\frac{{{{a}^{4}}\mathop {sin}\nolimits^2 {\theta }}}{{{\text{C}}{{{\rho }}^{2}}}}\left( {1 + \frac{{{\rho }({{u}_{1}}cos{\psi } + {{u}_{2}}sin{\psi })}}{{{{a}^{2}}{{{\omega }}_{3}}sin{\theta }}}} \right){{{\omega }}_{3}}} \\ {{\dot {\theta }} = - {\kappa }\frac{{{{a}^{2}}{{c}^{2}}sin{\theta }cos{\theta }}}{{C{{{\rho }}^{2}}}}\left( {1 + \frac{{{\rho }({{u}_{1}}cos{\psi } + {{u}_{2}}sin{\psi })}}{{{{a}^{2}}{{{\omega }}_{3}}sin{\theta }}}} \right)} \\ {\dot {U} = - \frac{{\kappa }}{m}U\left( {1 + \frac{{{{a}^{2}}{{{\omega }}_{3}}({{u}_{1}}cos{\psi } + {{u}_{2}}sin{\psi })sin{\theta }}}{{{\rho }{{U}^{2}}}}} \right)} \end{array}$Здесь ${{u}_{1}}$, ${{u}_{2}}$ ($u_{1}^{2} + u_{2}^{2} = {{U}^{2}}$) и ${\psi }$ – неизвестные функции времени.
Считая малой скорость центра масс (${\text{|}}U{\text{|}} = a{{{\omega }}_{3}}$) из (3.4) получим
(3.5)
${{{\dot {\omega }}}_{3}} = - {\kappa }\frac{{{{a}^{4}}\mathop {sin}\nolimits^2 {\theta }}}{{{{{\rho }}^{2}}C}}{{{\omega }}_{3}},\quad {\dot {\theta }} = - {\kappa }\frac{{{{a}^{2}}{{c}^{2}}}}{{2C{{{\rho }}^{2}}}}sin2{\theta }$Уравнения (3.4) в совокупности с равенствами (2.7) приближенно описывают динамику эллипсоида на плоскости с вязким трением при малой скорости центра масс и начальных условиях, близких к (2.7). Заметим, что при выполнении неравенства ${\text{|}}U{\text{|}} > a{{{\omega }}_{3}}$ правая часть последнего равенства (3.4) отрицательна и величина проекции скорости центра масс на опорную плоскость убывает, причем при ${\text{|}}U{\text{|}} \gg a{{{\omega }}_{3}}$ по закону $U\, = \,U(0){\text{exp}}( - {\kappa }t{\text{/}}m)$.
На фиг. 1 тонкими кривыми представлены результаты численного интегрирования точных уравнений движения (1.3)–(1.6) эллипсоида c параметрами $m = 1$ кг, $a = 0.05$ м, $c = 0.15$ м, на плоскости с коэффициентом вязкого трения ${\kappa } = 0.1$ кг/с при начальных условиях ${{{\omega }}_{1}}(0) = {{{\omega }}_{2}}(0) = 0$, ${{{\omega }}_{3}}(0) = 400$ c–1, ${{{\gamma }}_{2}}(0) = 0$, ${{{\gamma }}_{3}}(0) = 0.05$ (величина ${{{\gamma }}_{1}}(0)$ определялась геометрическим интегралом), т.е. в начальный момент времени эллипсоид имеет достаточно большую угловую скорость, направленную вдоль его оси симметрии, и располагался на плоскости так, что эта ось почти горизонтальна. В экспериментах, обозначенных цифрами $(1)$, $(2)$ и $(3)$ принимались значения ${{v}_{2}}(0) = 0$, ${{v}_{2}}(0)$ = = 20 м/с и ${{v}_{2}}(0) = - 20$ м/с соответственно, и ${{v}_{1}}(0) = {{v}_{3}}(0) = 0$, т.е. в случае (1) начальная скорость центра масс отсутствует, а в случаях (2) и (3) – направлена вдоль плоскости, причем в случае (2) в начальный момент отсутствует проскальзывание в точке контакта с плоскостью.
Во всех трех экспериментах с течением времени косинус угла нутации эллипсоида увеличивается до единицы, угловая скорость уменьшается, но остается практически сонаправленной с осью симметрии (зависимости ${{{\omega }}_{1}}(t)$ и ${{{\omega }}_{2}}(t)$ мало отличаются от нуля и не приводятся на рисунке), скорость центра масс уменьшается до нуля в случаях (2), (3) и мало отличается от нуля в случае (1). Таким образом, финальным движением эллипсоида являются устойчивые вращения (2.4) (u = 0).
Жирной кривой на фиг. 1 представлены результаты интегрирования приближенных уравнений (3.5) при соответствующих начальных условиях. Они хорошо совпадают с результатами эксперимента без начальной скорости центра масс (1), с осредненными значениями соответствующих переменных при начальном движении без проскальзывания (2) и с их финальными значениями в случае (3). Разность финальных значений угловой скорости в интегрировании точных и приближенных уравнений объясняется потерями энергии на переходных процессах.
На фиг. 2 представлена зависимость отношения величины реакции опорной плоскости к весу эллипсоида в эксперименте (2). В остальных экспериментах соответствующие результаты аналогичны. Отрыва эллипсоида от плоскости не наблюдается.
4. Геометрическая интерпретация результатов. Заметим, что при любой скорости центра масс из (3.4) следует уравнение
(4.1)
$\frac{{{{{\gamma }}_{3}}{\kern 1pt} {\omega }_{3}^{{{{c}^{2}}/{{a}^{2}}}}}}{{1 + \sqrt {1 - {\gamma }_{3}^{2}} }} = {\text{const}}$Таким образом при движениях эллипсоида с начальными условиями вида (2.7) при достаточно большой угловой скорости на плоскости параметров $({\text{|}}{{{\omega }}_{3}}{\text{|}},{\text{|}}{{{\gamma }}_{3}}{\text{|}})$ движение происходит в малой окрестности соответствующих кривых (4.1), заканчиваяcь на прямой ${\text{|}}{{{\gamma }}_{3}}{\text{|}} = 1$. Причем при малой скорости центра масс это движение монотонно.
На фиг. 3 пунктиром обозначены соответствующие равенствам (4.1) кривые, а сплошными кривыми – результаты соответствующих численных экспериментов без начальной скорости центра масс. При большой начальной угловой скорости эллипсоида указанные кривые практически совпадают. При ее уменьшении нарушается условие ${\varepsilon } \ll 1$ и наблюдается существенное различие результатов, а финальными движениями эллипсоида являются регулярные прецессии [4] или устойчивые качения (2.6).
5. Плоскость с поликомпонентным трением. При наличии трения верчения или качения (${\mathbf{M}} \ne 0$) система (1.3)–(1.6) имеет решения вида (2.6), (2.4), в которых следует положить $u = 0$, ${\omega } = {\omega }(t),$ где зависимость ${\omega }(t)$ определяется моделью трения. Уравнения (3.1), (3.2) в этом случае принимают вид
(5.1)
$C{{{\dot {\omega }}}_{3}} = ({{f}_{1}}{{{\gamma }}_{2}} - {{f}_{2}}{{{\gamma }}_{1}})\frac{{{{a}^{2}}}}{{\rho }} + {{m}_{3}},\quad \dot {K} = - \frac{{{{c}^{2}} - {{a}^{2}}}}{{\rho }}({{f}_{1}}{{{\gamma }}_{2}} - {{f}_{2}}{{{\gamma }}_{1}}){{{\gamma }}_{3}} + ({\mathbf{\gamma }},{\mathbf{M}})$Заметим, что если мощность трения верчения и качения мала по сравнению с мощностью трения скольжения (например, в моделях трения [7], [8]), то порядок последних слагаемых правых частей (5.1) мал по сравнению с первым, и приближенные уравнения (3.4) также будут достаточно хорошо описывать соответствующие переходные процессы. Таким образом, движение на плоскости $({\text{|}}{{{\omega }}_{3}}{\text{|}},{\text{|}}{{{\gamma }}_{3}}{\text{|}})$ также будет происходить в окрестности кривых (4.1), с той только разницей, что по достижении прямой ${\text{|}}{{{\gamma }}_{3}}{\text{|}} = 1$ изображающая точка не остановится, а будет двигаться вдоль этой прямой влево, а при пересечении с соответствующей прецессиям кривой будет двигаться вдоль нее в сторону уменьшения ${\text{|}}{{{\gamma }}_{3}}{\text{|}}$ [9].
На фиг. 4 представлены численные результаты эксперимента с начальными условиями (1) в уже рассмотренном случае вязкого трения (жирная черная кривая), при классическом сухом трении (жирная серая кривая), поликомпонентном сухом трении [10] (тонкая серая кривая) и поликомпонентном вязком трении [8] (тонкая черная кривая), а также соответствующие кривые на плоскости $({\text{|}}{{{\omega }}_{3}}{\text{|}},{\text{|}}{{{\gamma }}_{3}}{\text{|}})$.
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (16-01-00338, 18-01-00335, 19-01-00140) и Программы № 29 Президиума РАН.
Список литературы
Маркеев А.П. Динамика тела, соприкасающегося с твердой поверхностью. М.; Ижевск: Ин-т компьют. исслед., 2014. 496 с.
Буров А.А., Карапетян А.В. О несуществовании дополнительного интеграла в задаче о движении тяжелого твердого эллипсоида по гладкой плоскости // ПММ. 1985. Т. 49. Вып. 3. С. 501.
Карапетян А.В., Рубановский В.Н. О бифуркации и устойчивости перманентных вращений тяжелого трехосного эллипсоида на гладкой плоскости // ПММ. 1987. Т. 51. Вып. 2. С. 260.
Муницына М.А. Движения сфероида на горизонтальной плоскости с вязким трением // ПММ. Т. 76. Вып. 2. 2012. С. 214–223.
Моргунов Б.И., Волосов В.М. Метод осреднения в теории нелинейных колебательных систем. М.: МГУ, 1971. 508 с.
Блехман И.И. Вибрационная механика. М.: Физматлит, 1994. 400 с.
Журавлев В.Ф. О модели сухого трения в задаче качения твердых тел // ПММ. 1998. Т. 62. Вып. 5. С. 762–767.
Карапетян А.В. Двухпараметрическая модель трения и ее свойства // ПММ. 2009. Т. 73. В. 5. С. 515–519.
Муницына М.А. Динамика вращения сфероида на горизонтальной плоскости с трением // ПММ. 2018. Т. 82. Вып. 1. С. 1–9.
Муницына М.А. Модель трения в случае плоского эллиптического контакта тела с опорной плоскость // Нелинейная динамика. Т. 8. Вып. 4. 2012. С. 705–712.
Муницына М.А. Двухкомпонентная модель вязкого трения // ХI всероссийский съезд по фундаментальным проблемам теоретической и прикладной механики: сборник докладов (Казань, 20–24 августа 2015 г.). С. 2668–2670.
Дополнительные материалы отсутствуют.
Инструменты
Известия РАН. Механика твердого тела