Известия РАН. Механика твердого тела, 2023, № 3, стр. 99-105

ОБ ОТНОСИТЕЛЬНЫХ РАВНОВЕСИЯХ ТЯЖЕЛОЙ БУСИНКИ НА РАВНОМЕРНО ВРАЩАЮЩЕЙСЯ ШЕРОХОВАТОЙ СПИЦЕ

А. А. Буров a*, В. И. Никонов a**

a Федеральный исследовательский центр “Информатика и управление” Российской академии наук
Москва, Россия

* E-mail: jtm@narod.ru
** E-mail: nikon_v@list.ru

Поступила в редакцию 21.04.2022
После доработки 29.06.2022
Принята к публикации 30.06.2022

Полный текст (PDF)

Аннотация

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

Ключевые слова: относительные равновесия, сухое трение, равномерно вращающаяся система отсчета, шероховатая спица, тяжелая бусинка

1. Постановка задачи. Пусть шероховатая прямолинейная спица $\ell $ совершает равномерное вращение с постоянной угловой скоростью $\omega > 0$ вокруг вертикальной оси, отстоящей от спицы на расстоянии a. Предположим, что на спицу нанизана тяжелая бусинка P. Требуется найти положения равновесия бусинки относительно спицы в предположении о том, что коэффициент сухого трения равен $\mu = {\text{tg}}\alpha $, где $0 \leqslant \alpha < \pi {\text{/}}2$ – угол трения. Ускорение силы тяжести считается равным g. Предполагается, что спица образует с восходящей вертикалью угол $0 \leqslant \theta \leqslant \pi {\text{/}}2$, см. рис. 1.

Рис. 1.

Механическая система: вращающаяся шероховатая спица $\ell $ с нанизанной на ней бусинкой.

2. Общие условия существования относительных равновесий. В равномерно вращающейся вместе со спицей системе отсчета бусинка будет находиться в равновесии под действием (рис. 1) силы тяжести $m{\mathbf{g}}$, центробежной силы ${{{\mathbf{F}}}_{c}} = m{{{\mathbf{g}}}_{c}}$, а также реакции связи, состоящей из нормальной компоненты ${\mathbf{N}}$ (не изображена на рисунке) и силы трения ${\mathbf{T}} = T{\mathbf{e}}$, где ${\mathbf{e}}$ – единичный вектор, направленный вверх вдоль спицы. Пусть масса бусинки принята за единицу. Тогда, если ${\mathbf{f}} = {{{\mathbf{g}}}_{c}} + {\mathbf{g}}$, то условие равновесия запишется в виде

(2.1)
${\mathbf{f}} + {\mathbf{N}} + {\mathbf{T}} = 0$

Согласно закону Кулона–Амонтона величины нормальной реакции и силы трения стеснены неравенством

${\text{|}}{\mathbf{T}}{\text{|}} \leqslant \mu {\text{|}}{\mathbf{N}}{\text{|}} \Leftrightarrow {\text{|}}{\mathbf{T}}{\text{|}} - \mu {\text{|}}{\mathbf{N}}{\text{|}} \leqslant 0$

Домножая левую и правую части последнего неравенства на неотрицательную11 величину ${\text{|}}{\mathbf{T}}{\text{|}} + \mu {\text{|}}{\mathbf{N}}{\text{|}}$, получим

${\text{|}}{\mathbf{T}}{{{\text{|}}}^{2}} - {{\mu }^{2}}{\text{|}}{\mathbf{N}}{{{\text{|}}}^{2}} \leqslant 0 \Leftrightarrow ({\mathbf{T}},{\mathbf{T}}) - {{\mu }^{2}}({\mathbf{N}},{\mathbf{N}}) \leqslant 0$

Здесь и далее (·, ·) – скалярное произведение. Принимая во внимание ортогональность векторов ${\mathbf{T}}$ и ${\mathbf{N}}$, последнее из неравенств запишем в виде

(2.2)
$({\mathbf{T}} + {\mathbf{N}},{\mathbf{T}} + {\mathbf{N}}) - (1 + {{\mu }^{2}})({\mathbf{N}},{\mathbf{N}}) \leqslant 0$

Из уравнения (2.1) находим

${\mathbf{T}} + {\mathbf{N}} = - {\mathbf{f}}$
и неравенство (2.2) может быть записано как

(2.3)
$({\mathbf{f}},{\mathbf{f}}) - (1 + {{\mu }^{2}})({\mathbf{N}},{\mathbf{N}}) \leqslant 0$

Далее, путем скалярного умножения левой и правой частей уравнения (2.1) на ${\mathbf{e}}$ и дальнейших преобразований, находим величину $T$:

$({\mathbf{e}},{\mathbf{f}}) + ({\mathbf{e}},{\mathbf{T}}) = ({\mathbf{e}},{\mathbf{f}}) + T = 0 \Rightarrow T = - ({\mathbf{e}},{\mathbf{f}})$

Так как

${\mathbf{N}} = - {\mathbf{T}} - {\mathbf{f}} = ({\mathbf{e}},{\mathbf{f}}){\mathbf{e}} - {\mathbf{f}}$
то
$({\mathbf{N}},{\mathbf{N}}) = ({\mathbf{f}},{\mathbf{f}}) - {{({\mathbf{f}},{\mathbf{e}})}^{2}}$
и неравенство (2.3) приобретает вид

(2.4)
$({\mathbf{f}},{\mathbf{f}}) - (1 + {{\mu }^{2}})(({\mathbf{f}},{\mathbf{f}}) - {{({\mathbf{f}},{\mathbf{e}})}^{2}}) \leqslant 0 \Leftrightarrow (1 + {{\mu }^{2}})({\mathbf{f}},{\mathbf{e}}{{)}^{2}} - {{\mu }^{2}}({\mathbf{f}},{\mathbf{f}}) \leqslant 0$

Это неравенство составляет основной предмет дальнейших исследований.

3. Зависимость относительных равновесий от параметров. Пусть $O{{P}_{0}}$ – общий перпендикуляр к оси вращения и спице, причем точка O принадлежит оси вращения, а точка ${{P}_{0}}$ – спице. Введем вращающуюся вместе со спицей систему отсчета $O{{x}_{1}}{{x}_{2}}{{x}_{3}}$, ось $O{{x}_{3}}$ которой направлена вдоль оси вращения, ось $O{{x}_{1}}$ – вдоль общего перпендикуляра, а ось $O{{x}_{2}}$ дополняет их до правой тройки. При этом ${{\overleftarrow {OP} }_{0}} = (a{{,0,0)}^{T}}$, e = (0, sinθ, ${\text{cos}}\theta {{)}^{T}}$, ${\mathbf{g}} = (0,0, - g{{)}^{T}}$. Вводя параметр $p \in {\mathbf{R}}$, отвечающий за удаление бусинки P от точки ${{P}_{0}}$ на прямой, зададим ее положение вектором

$\overleftarrow {OP} = ({{x}_{1}},{{x}_{2}},{{x}_{3}}{{)}^{T}} = {{\overleftarrow {OP} }_{0}} + p{\mathbf{e}} = (a,p\sin \theta ,p\cos \theta {{)}^{T}}$

Тогда ${{{\mathbf{g}}}_{c}} = {{\omega }^{2}}{{(a,p\sin \theta ,0)}^{T}},$ ${\mathbf{f}} = ({{\omega }^{2}}a,{{\omega }^{2}}p\sin \theta , - g{{)}^{T}},$

$({\mathbf{f}},{\mathbf{e}}) = {{\omega }^{2}}p{{\sin }^{2}}\theta - g\cos \theta ,\quad ({\mathbf{f}},{\mathbf{f}}) = {{\omega }^{4}}({{a}^{2}} + {{p}^{2}}{{\sin }^{2}}\theta ) + {{g}^{2}}$
и неравенство (2.4) сводится к неравенству
$(1 + {{\mu }^{2}}){{({{\omega }^{2}}p{{\sin }^{2}}\theta - g\cos \theta )}^{2}} - {{\mu }^{2}}({{\omega }^{4}}({{a}^{2}} + {{p}^{2}}{{\sin }^{2}}\theta ) + {{g}^{2}}) \leqslant 0$
приводимому к виду

$\begin{gathered} {{\omega }^{4}}{{\sin }^{2}}\theta ({{\cos }^{2}}\alpha - {{\cos }^{2}}\theta ){{p}^{2}} - 2{{\omega }^{2}}g\cos \theta {{\sin }^{2}}\theta p + \\ \, + {{g}^{2}}({{\cos }^{2}}\alpha - {{\sin }^{2}}\theta ) - {{a}^{2}}{{\omega }^{4}}{{\sin }^{2}}\alpha \leqslant 0 \\ \end{gathered} $

Введение безразмерной координаты $p{\kern 1pt} '{\text{:}}$ $p = ap{\kern 1pt} '$ и дальнейшее отбрасывание штриха позволяют представить это неравенство в виде

$A{{p}^{2}} + Bp + C \leqslant 0$
(3.1)
$A = {{{{\unicode{230} }}}^{2}}{{\sin }^{2}}\theta ({{\cos }^{2}}\alpha - {{\cos }^{2}}\theta ),\quad B = - 2{{\unicode{230} }}\cos \theta {{\sin }^{2}}\theta $
$C = {{\cos }^{2}}\theta - (1 + {{{{\unicode{230} }}}^{2}}){{\sin }^{2}}\alpha ,\quad {{\unicode{230} }} = \frac{{a{{\omega }^{2}}}}{g} \geqslant 0$

Решения неравенства (3.1) относительно p определяют на спице $\ell $ области, заполненные относительными равновесиями (ОЗОР). Зависимость этих областей от безразмерного параметра ${{\unicode{230} }}$ и от углов α и θ составляет предмет дальнейших исследований.

а) Если $A = 0$ при ${{\unicode{230} }} = 0$, т.е. спица не вращается или пересекается с осью вращения, то неравенство (3.1) записывается как

(3.2)
${{\cos }^{2}}\theta - {{\sin }^{2}}\alpha \leqslant 0$

Неравенство выполнено при $\pi {\text{/}}2 - \theta \leqslant \alpha $ и не выполнено в противном случае. Это классический результат, иллюстрирующий смысл угла трения α [1] (см. также [2]). б) Если $A = 0$ при $\theta = 0$, т.е. спица вертикальна, то неравенство (3.1) приводится к виду

(3.3)
${{\cos }^{2}}\alpha - {{{{\unicode{230} }}}^{2}}{{\sin }^{2}}\alpha \leqslant 0$
которое выполнено при $0\; \leqslant {\text{tg}}\alpha \leqslant {{{{\unicode{230} }}}^{{ - 1}}}$.

в) Если $A = 0$ при ${{\cos }^{2}}\alpha - {{\cos }^{2}}\theta = 0$, то неравенство (3.1) приводится к виду

(3.4)
$ - 2{{\unicode{230} }}\cos \theta {{\sin }^{2}}\theta p + {{\cos }^{2}}\theta - (1 + {{{{\unicode{230} }}}^{2}}){{\sin }^{2}}\theta \leqslant 0$

Тогда, так как $0 \leqslant \alpha < \pi {\text{/}}2$, то и $0 \leqslant \theta < \pi {\text{/}}2$, и $\cos \theta \ne 0$. Случай θ = 0 разобран выше, поэтому решение неравенства (3.4) имеет вид

$p \geqslant {{p}_{ \star }},\quad {{p}_{ \star }} = \frac{{{{s}^{2}} - {{{\sin }}^{2}}\theta }}{{2{{\unicode{230} }}{{s}^{2}}{{{\sin }}^{2}}\theta \cos \theta }},\quad s = \frac{1}{{\sqrt {2 + {{{{\unicode{230} }}}^{2}}} }} \leqslant \frac{{\sqrt 2 }}{2}$

Граничная точка ${{P}_{ \star }}$, определяемая значением $p = {{p}_{ \star }}$,

• совпадает с точкой P0 при $\sin \theta = s$,

• располагается выше точки P0 при $\sin \theta \in (0;s)$,

• располагается ниже точки P0 при $\sin \theta \in \left( {s,1} \right]$.

г) Если $A < 0$, т.е. ${{\cos }^{2}}\alpha - {{\cos }^{2}}\theta < 0$, то в левой части неравенства (3.1) стоит квадратный трехчлен относительно p, дискриминант которого имеет вид

$D = {{\cos }^{2}}\alpha + {{{{\unicode{230} }}}^{2}}({{\cos }^{2}}\alpha - {{\cos }^{2}}\theta )$

Неравенство (3.1) при этом приводится к виду

$(p - {{p}_{ - }})(p - {{p}_{ + }}) \geqslant 0,\quad {{p}_{ \pm }} = \frac{{\sin \theta \cos \theta \pm \sin \alpha \sqrt D }}{{{{\unicode{230} }}\sin \theta ({{{\cos }}^{2}}\alpha - {{{\cos }}^{2}}\theta )}}$

Величины ${{p}_{ \pm }}$ вещественны если

$D \geqslant 0 \Leftrightarrow {{\cos }^{2}}\alpha < {{\cos }^{2}}\theta \leqslant {{\cos }^{2}}\alpha \left( {\frac{1}{{{{{{\unicode{230} }}}^{2}}}} + 1} \right) = K$

Если K ≥ 1, то угол θ < α, если $K < 1$, то используя обозначение $\beta = \arccos K$, имеем

${{\cos }^{2}}\alpha < {{\cos }^{2}}\theta \leqslant {{\cos }^{2}}\beta \Leftrightarrow \beta \leqslant \theta < \alpha $

В этом случае относительные равновесия занимают область, отвечающую значениям параметра

$p{\kern 1pt} :\; - \infty < p \leqslant {{p}_{ - }} \cup {{p}_{ + }} \leqslant p < + \infty $
т.е. располагаются вне интервала $({{p}_{ - }};{{p}_{ + }})$.

При $D < 0$ вещественных корней нет, а с учетом $A < 0$, неравенство (3.1) выполнено при любых значениях $p$, т.е. ОЗОР – все точки спицы.

д) При выполнении условия $A > 0$, т.е. ${{\cos }^{2}}\alpha - {{\cos }^{2}}\theta > 0$, неравенство (3.1) имеет вид

$(p - {{p}_{ - }})(p - {{p}_{ + }}) \leqslant 0$

Относительные равновесия занимают область, отвечающую значениям параметра

$p{\kern 1pt} :\;{{p}_{ - }} \leqslant p \leqslant {{p}_{ + }}$

Поскольку ${{\cos }^{2}}\alpha - {{\cos }^{2}}\theta > 0$, то $D \geqslant 0$ и величины ${{p}_{ \pm }}$ вещественны.

4. Графическое представление зависимости относительных равновесий от параметров. Для иллюстрации различных случаев расположения ОЗОР на плоскости (θ, α) изобразим три кривые. Кривая ${{\Gamma }_{1}}$: $D = 0$ отделяющую область, в которой изучаемый квадратный трехчлен имеет два вещественных корня, от области, в которой вещественных корней нет. ${{\Gamma }_{2}}$: $C = 0$, разделяет области, где эти корни одного или разных знаков и, наконец, ${{\Gamma }_{3}}$: $A = 0$, разделяет области, где ветви изображающей трехчлен параболы направлены вниз или вверх. Кривые ${{\Gamma }_{1}}$, ${{\Gamma }_{2}}$ и ${{\Gamma }_{3}}$ разбивают квадрат S = {0 ≤ θ ≤ π/2, $0 \leqslant \alpha \leqslant \pi {\text{/}}2\} $ на пять областей $I.$$V.$ При этом

${{\alpha }_{ \star }} = {{\Theta }_{ \star }} = \arccos \sqrt {\frac{{1 + {{{{\unicode{230} }}}^{2}}}}{{2 + {{{{\unicode{230} }}}^{2}}}}} ,\quad {{\alpha }_{{ \star \star }}} = \arctan \frac{1}{{{\unicode{230} }}}$
${{\left. {{{\alpha }_{ \star }}} \right|}_{{{{\unicode{230} }} = 0}}} = \pi {\text{/}}4,\quad {{\left. {{{\alpha }_{ \star }}} \right|}_{{{{\unicode{230} }} \to + \infty }}} = 0$

Для ${{\unicode{230} }} = 1{\text{/}}2$ эти области изображены на рис. 2.

Рис. 2.

Области $I.$$V.$ на плоскости (θ, α), отвечающие различным расположениям ОЗОР, при ${{\unicode{230} }} = 1{\text{/}}2$.

При неограниченном увеличении параметра ${{\unicode{230} }}$, области $I.$, $III.$ и $V.$ (см. рис. 2) уменьшаются, исчезая в пределе. При этом кривая ${{\Gamma }_{2}}$ сливается с граничным отрезком α = 0, а кривая ${{\Gamma }_{1}}$ приближается к диагонали α = θ. В предельном случае, оставшиеся области $II.$ и $IV.$ разбивают квадрат S на равные треугольники с общим основанием на диагонали квадрата S, соответствующей кривой ${{\Gamma }_{3}}$.

При неограниченном уменьшении параметра ${{\unicode{230} }}$ до нуля область $II.$ уменьшается и сливается с граничным отрезком $\alpha = \pi {\text{/}}2$. При этом кривая ${{\Gamma }_{2}}$ приближается к диагонали $\alpha = \pi {\text{/}}2 - \theta $. В пределе, при ${{\unicode{230} }} = 0$, т.е. в случае а), области $I.$, $III.$, $IV.$ и $V.$ разбивают квадрат S на равные треугольники с общей вершиной и основаниями, противоположными этой вершине, на сторонах квадрата S. При этих условиях неравенство (3.2) выполнено в точках диагонали $\alpha = \pi {\text{/}}2 - \theta $.

В случае б) неравенство (3.3) выполнено в точках отрезка θ = 0, $\alpha \in [0,{{\alpha }_{{ \star \star }}}]$.

В случае в) неравенство (3.4) выполнено в точках кривой ${{\Gamma }_{3}}$.

Случаю г) отвечают области $I.$$III.$

На рис. 3 серым закрашены области, где существуют ОЗОР. Области построены для случае ${{\unicode{230} }} = 0.5$, на рисунках схематично изображена спица $\ell $, что позволяет проиллюстрировать области ОЗОР. Пунктирной линией отмечены точки границы при ${{p}_{ - }}({{\unicode{230} }})$ = 0, штриховой линией отмечены точки границы при ${{p}_{ + }}({{\unicode{230} }}) = 0$. Точки $({{{{\unicode{230} }}}_{i}},{{p}_{i}})$, $i = I,II$, III отвечают случаям, когда ${{p}_{ - }}({{\unicode{230} }}) = {{p}_{ + }}({{\unicode{230} }})$. Для области $I.$ ОЗОР представляется в виде объединения двух лучей, направленных вверх и вниз по спице, при этом начало верхнего луча лежит выше точки ${{P}_{0}}$, начало нижнего луча лежит ниже точки ${{P}_{0}}$. Для области $II.$ вся спица заполнена ОЗОР. Для области $III.$ ОЗОР опять представляется в виде объединения двух лучей, направленных вверх и вниз по спице, но при этом начала обоих лучей лежат ниже точки ${{P}_{0}}$, т.е. сама точка ${{P}_{0}}$ принадлежит ОЗОР.

Случаю д) отвечают области $IV.$ и $V.$, см. рис. 3. Для области $IV.$ ОЗОР представляется в виде отрезка, концы которого лежат по разные стороны от точки P0, и точка P0 принадлежит ОЗОР. Для области V., ОЗОР также является отрезком, но целиком лежащим выше точки P0.

5. Выводы. Рассмотренная в учебнике А. Зоммерфельда ([1], стр. 73) задача о тяжелой бусинке на вращающейся спице представляет собой замечательный пример системы, в которой при прохождении параметра, например, коэффициента трения, через определенное значение, множество неизолированных относительных равновесий, заполняющее целую прямую, исчезает22. Как показало настоящее исследование, выполненное в продолжение ряда исследований по существованию, устойчивости и бифуркациям равновесий в системах с трением (см., например, [413]), в рамках несложного обобщения задачи А. Зоммерфельда обнаруживается достаточно богатая (с точки зрения динамических эффектов) картина бифуркаций относительных равновесий.

Рис. 3.

Серым закрашена область, где существуют ОЗОР. $\ell $ – спица. Параметры (θ, α) приняты следующими: $I.$: $(\pi {\text{/}}8,\pi {\text{/}}4)$, $II.$: $(\pi {\text{/}}8,4\pi {\text{/}}10)$, $III.$: $(\pi {\text{/}}4,\pi {\text{/}}3)$, $IV.$: $(3\pi {\text{/}}8,\pi {\text{/}}4)$, $IV.$: $(\pi {\text{/}}4,\pi {\text{/}}8)$. Точки сращивания корней на рисунках $I.$$III.$ ${{{{\unicode{230} }}}_{{I.}}} = 1.189207115$, ${{p}_{{I.}}} = - 2.197368226$; ${{{{\unicode{230} }}}_{{II.}}} = 0.3549196403$, pII. = $ - 3.433844707$; ${{{{\unicode{230} }}}_{{III.}}} = 1$, ${{p}_{{III.}}} = - 2.828427124$.

Список литературы

  1. Sommerfeld A. Vorlesungen über Theoretische Physik, Band 1: Mechanik. Leipzig: Akademische Verlagsgesellschaft. 1943.

  2. Розенблат Г.М. Сухое трение и односторонние связи в механике твердого тела. М.: URSS. 2011. 208 с.

  3. Painlevé P. Leçons sur le frottement. Paris: Hermann. 1895. 111 p.

  4. Leine R.I., van Campen D.H. Bifurcation phenomena in non-smooth dynamical systems // Europ. J. Mechanics – A/Solids. 2006. V. 25. P. 595–616. https://doi.org/10.1016/j.euromechsol.2006.04.004

  5. Leine R.I. Bifurcations of equilibria in non-smooth continuous systems // Physica D: Nonlinear Phenomena. 2006. V. 223. № 1. P. 121–137. https://doi.org/10.1016/j.physd.2006.08.021

  6. Иванов А.П. Об устойчивости равновесия в системах с трением // Прикладная математика и механика. 2007. Т. 71. № 3. С. 427–438.

  7. Иванов А.П. Бифуркации в системах с трением: основные модели и методы // Нелинейная динамика. 2009. Т. 5. № 4. С. 479–498.

  8. Burov A.A. On bifurcations of relative equilibria of a heavy bead sliding with dry friction on a rotating circle // Acta Mechanica. 2010. V. 212. № 3–4. P. 349–354. https://doi.org/10.1007/s00707-009-0265-1

  9. Иванов А.П. Основы теории систем с трением. Ижевск: РХД. 2011. 302 с.

  10. Иванов А.П. О равновесии систем с сухим трением // Прикладная математика и механика. 2015. Т. 79. № 3. С. 317–333.

  11. Шалимова Е.С. О движении тяжелой точки по сфере, вращающейся вокруг не проходящей через ее центр вертикальной оси, при наличии сухого трения // Нелинейная динамика. 2016. Т. 12. № 3. С. 369–383. https://doi.org/10.20537/nd1603006

  12. Буров А.А., Косенко И.И., Шалимова Е.С. Об относительных равновесиях массивной точки на равномерно вращающемся астероиде // Доклады Академии наук. 2017. Т. 475. № 3. С. 269–272. https://doi.org/10.7868/S0869565217210071

  13. Буров А.А., Никонов В.И., Шалимова Е.С. Движение массивной точки по поверхности однородного шара со сферической полостью // Прикладная математика и механика. 2021. Т. 85. № 4. С. 528–543. https://doi.org/10.31857/S0032823521040032

Дополнительные материалы отсутствуют.