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

ПОСТАНОВКА ЗАДАЧИ О ПРЕЦЕССИЯХ ТВЕРДОГО ТЕЛА С НЕПОДВИЖНОЙ ТОЧКОЙ В ТРЕХ ОДНОРОДНЫХ СИЛОВЫХ ПОЛЯХ. ПРЕЦЕССИОННО-ИЗОКОНИЧЕСКИЕ ДВИЖЕНИЯ ТЕЛА

Г. В. Горр a*

a Государственное бюджетное учреждение “Институт прикладной математики и механики”
Донецк, Россия

* E-mail: gvgorr@gmail.com

Поступила в редакцию 11.08.2022
После доработки 30.08.2022
Принята к публикации 31.08.2022

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

Аннотация

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

Ключевые слова: твердое тело, суперпозиция трех однородных силовых полей, прецессии, прецессионно-изоконические движения

1. Введение. Прецессионные движения твердого тела, имеющего неподвижную точку, в частности, регулярные прецессии, применяются в моделировании технических объектов [1]. В динамике тяжелого твердого тела и ее обобщениях на случай сложных силовых полей они изучались Д. Гриоли [2], Ф. Кляйном и А. Зоммерфельдом [3], А. Брессаном [4], автором данной статьи [5] и другими [68]. Прецессионные движения характеризуются свойством постоянства угла между двумя осями ${{l}_{1}},\;{{l}_{2}}$, проходящими через неподвижную точку, одна из которых (${{l}_{1}}$) связана с телом-носителем, а другая (${{l}_{2}}$) неподвижна в пространстве. В случае, когда одна из осей подвижной системы координат содержит ось ${{l}_{1}}$, то целесообразно такую систему координат называть прецессионной системой координат. Согласно [24], прецессионные движения подразделяются на классы: если скорости прецессии и собственного вращения постоянны, то прецессия называется регулярной; если только скорость прецессии постоянна, то прецессия называется полурегулярной прецессией первого типа; если только скорость собственного вращения постоянна, то прецессию называют полурегулярной прецессией второго типа; в других случаях прецессия называется прецессией общего типа. Наибольшее количество прецессионных движений гиростата установлено для классов регулярных и полурегулярных прецессий первого типа. Следует отметить, что уникальными случаями прецессий тяжелого твердого тела, описываемых уравнениями Эйлера–Пуассона, являются случай регулярной прецессии Д. Гриоли [2], относительно наклонной оси, случай А.И. Докшевича [9], для которого постоянно произведение скоростей прецессии и собственного вращения, прецессия А. Брессана [4] в решении В. Гесса. Полурегулярные прецессии второго типа твердого тела и гиростата относятся к меньшему числу найденных прецессий. Например, в [4] доказано, что в классической задаче полурегулярные прецессии второго типа динамически невозможны. Несмотря на это, в задаче о движении гиростата с переменным гиростатическим моментом получены некоторые решения, имеющие такое свойство [8].

Отметим, что прецессии несимметричного твердого тела с жидким заполнением изучались В.Ю. Ольшанским [1013], а регулярные прецессии тяжелого твердого тела в жидкости – В.Н. Рубановским [14]. В задаче о движении системы связанных твердых тел также установлен нетривиальный класс прецессий [15].

Анализ литературы, посвященной изучению задачи о движении гиростата в двух однородных силовых полях, показывает, что большинство статей посвящено качественным проблемам интегрируемости уравнений движения [1619] с использованием аналогии [20]. Однако, существуют и некоторые исследования, в которых рассматриваются регулярные прецессии гиростата [2123] в указанной задаче.

В данной статье дана постановка задачи о прецессиях твердого тела с неподвижной точкой в трех однородных силовых полях. Записаны уравнения движения, описывающие прецессионные движения твердого тела, в качестве подвижной системы координат принята прецессионная система координат [24] (одна из осей подвижной системы координат направлена по оси ${{l}_{1}}$). Редуцированные уравнения движения получены в независимом базисе ${\mathbf{a}} \in {{l}_{1}}$, ${\mathbf{\gamma }} \in {{l}_{2}}$ и ${\mathbf{a}} \times {\mathbf{\gamma }}$. Их интегрирование рассмотрено в случае, когда совпадают скорости прецессии и собственного вращения. Данное свойство характеризует изоконическое движение твердого тела (подвижный и неподвижный годографы вектора угловой скорости симметричны друг другу относительно касательной плоскости к соответствующим аксоидам). Таким образом, по терминологии [6] в данной статье рассматриваются условия существования прецессионно-изоконических движений твердого тела в трех однородных силовых полях. Следует отметить, что аналогами этих движений являются движение гироскопа Лагранжа с равными скоростями прецессии и собственного вращения, а также движение гироскопа Гриоли. Для случая, когда эллипсоид инерции тела является сферой, в статье показано, что угол нутации равен ${\pi \mathord{\left/ {\vphantom {\pi 3}} \right. \kern-0em} 3}$.

2. Постановка задачи. Рассмотрим движение твердого тела, имеющего неподвижную точку, в силовом поле, которое является суперпозицией трех однородных и постоянных силовых полей. Обозначим через ${\mathbf{\gamma }}{\text{,}}{{{\mathbf{\gamma }}}^{{(1)}}},{{{\mathbf{\gamma }}}^{{(2)}}}$ единичные векторы, характеризующие направления сил ${\mathbf{P}},{\mathbf{P}}_{1}^{{}}{\text{,}}{\mathbf{P}}_{2}^{{}}$ каждого из полей; $C,C_{1}^{{}},C_{2}^{{}}$ – центры приведения сил; ${\mathbf{s}} = P{\mathbf{OC}}$, ${\mathbf{r}} = {{{\mathbf{P}}}_{1}}{\mathbf{O}}{{{\mathbf{C}}}_{1}}$, ${\mathbf{p}} = P_{2}^{{}}{\mathbf{OC}}_{2}^{{}}$; $Oxyz$ – подвижная система координат, $O$ – неподвижная точка. Пусть тензор инерции тела в системе $Oxyz$ имеет значение $A = {\text{(}}A_{{ij}}^{{}}{\text{)}}$ $(i,j = \overline {1,3} )$. Тело вращается вокруг точки $O$ с угловой скоростью ω = (ω1i1 + $\omega _{2}^{{}}{\mathbf{i}}_{2}^{{}} + \omega _{3}^{{}}{\mathbf{i}}_{3}^{{}}{\text{)}}$ (${\mathbf{i}}_{1}^{{}},{\mathbf{i}}_{2}^{{}},{\mathbf{i}}_{3}^{{}}$ – единичные векторы системы $Oxyz$). Для векторов s, r, p запишем соотношения

(2.1)
${\mathbf{s}} = s_{1}^{{}}{\mathbf{i}}_{1}^{{}} + s_{2}^{{}}{\mathbf{i}}_{2}^{{}} + s_{3}^{{}}{\mathbf{i}}_{3}^{{}}{\text{,}}\quad {\mathbf{r}}{\text{ }} = r_{1}^{{}}{\mathbf{i}}_{1}^{{}} + r_{2}^{{}}{\mathbf{i}}_{2}^{{}} + r_{3}^{{}}{\mathbf{i}}_{3}^{{}}{\text{,}}\quad {\mathbf{p}} = p_{1}^{{}}{\mathbf{i}}_{1}^{{}} + p_{2}^{{}}{\mathbf{i}}_{2}^{{}} + p_{3}^{{}}{\mathbf{i}}_{3}^{{}}$

Тогда уравнения движения тела представим в виде [21]

(2.2)
$A{\mathbf{\dot {\omega }}} = A{\mathbf{\omega }} \times {\mathbf{\omega }} + {\mathbf{s}} \times {\mathbf{\gamma }} + {\mathbf{r}} \times {\mathbf{\gamma }}_{{}}^{{(1)}} + {\mathbf{p}} \times {\mathbf{\gamma }}_{{}}^{{(2)}}$
(2.3)
${\mathbf{\dot {\gamma }}} = {\mathbf{\gamma }} \times {\mathbf{\omega }},\quad {\mathbf{\dot {\gamma }}}_{{}}^{{(1)}} = {\mathbf{\gamma }}_{{}}^{{(1)}} \times {\mathbf{\omega }},\quad {\mathbf{\dot {\gamma }}}_{{}}^{{(2)}} = {\mathbf{\gamma }}_{{}}^{{(2)}} \times {\mathbf{\omega }}$
где точка над переменными ${\mathbf{\omega }},\;{\mathbf{\gamma }},\;{\mathbf{\gamma }}_{{}}^{{(1)}},\;{\mathbf{\gamma }}_{{}}^{{(2)}}$ обозначает дифференцирование по времени $t$. В формулах (2.2), (2.3) полагаем
(2.4)
${\mathbf{\gamma }} \cdot {\mathbf{\gamma }}_{{}}^{{(1)}} = 0,\quad {\mathbf{\gamma }}_{{}}^{{(2)}} = {\mathbf{\gamma }} \times {\mathbf{\gamma }}_{{}}^{{(1)}},\quad {\text{|}}{\mathbf{\gamma }}{\text{|}} = 1,\quad {\text{|}}{\mathbf{\gamma }}_{{}}^{{(1)}}{\text{|}} = 1$
то есть направления силовых полей будут характеризоваться тройкой ${\mathbf{\gamma }},\;{\mathbf{\gamma }}_{{}}^{{(i)}}$ $(i = 1,2)$. Тогда очевидны равенства ${\mathbf{P}} = P\;{\mathbf{\gamma }}$, ${\mathbf{P}}_{i}^{{}} = P_{i}^{{}}{\mathbf{\gamma }}_{{}}^{{(i)}}$ $(i = 1,2)$.

Рассмотрим прецессии тела относительно вектора ${\mathbf{\gamma }}$. Они характеризуются инвариантным соотношением (ИС) [6, 7]

(2.5)
${\mathbf{a}} \cdot {\mathbf{\gamma }} = a_{{\text{0}}}^{{}}\quad (a_{{\text{0}}}^{{}} = \cos \theta _{{\text{0}}}^{{}})$
где $\theta _{{\text{0}}}^{{}}$ – угол между векторами ${\mathbf{a}}$ и ${\mathbf{\gamma }}$ (${\mathbf{\dot {a}}} = {\mathbf{0}}$, ${\text{|}}{\mathbf{a}}{\text{|}} = 1$). Вектор угловой скорости тела на ИС (2.5) представим так [7]:

(2.6)
${\mathbf{\omega }} = \dot {\varphi }{\mathbf{a}} + \dot {\psi }{\mathbf{\gamma }}$

Переменные $\varphi {\text{,}}\;{{\psi }}$ и постоянную $\theta _{{\text{0}}}^{{}}$ можно трактовать, как углы Эйлера. Используя метод [7], запишем значение вектора ${{{\mathbf{\gamma }}}^{{(1)}}}$:

(2.7)
${{{\mathbf{\gamma }}}^{{(1)}}} = b_{0}^{{}}[a_{{\text{0}}}^{{}}{\mathbf{\gamma }}\sin (\psi + \psi _{0}^{{}}) - {\mathbf{a}}\sin (\psi + \psi _{0}^{{}}) + ({\mathbf{a}} \times {\mathbf{\gamma }})\cos (\psi + \psi _{0}^{{}})]$
где $b_{0}^{{}} = 1{\text{/}}a_{0}^{'}$ $(a_{0}^{'} = \sin \theta _{{\text{0}}}^{{}})$, $\psi _{0}^{{}}$ – постоянная.

Значение вектора ${{{\mathbf{\gamma }}}^{{(2)}}}$ найдем по второй формуле системы (2.4)

(2.8)
${\mathbf{\gamma }}_{{}}^{{(2)}} = b_{0}^{{}}[{\mathbf{a}}\cos (\psi + \psi _{0}^{{}}) - a_{{\text{0}}}^{{}}{\mathbf{\gamma }}\cos (\psi + \psi _{0}^{{}}) + ({\mathbf{a}} \times {\mathbf{\gamma }})\sin (\psi + \psi _{0}^{{}})]$

Таким образом, при получении (2.7), (2.8) полагалось, что ${\mathbf{a}} \times {\mathbf{\gamma }} \ne {\mathbf{0}}$, то есть случай равномерных вращений тела исключаем из рассмотрения. Подвижную систему координат выберем следующим образом: направим вектор i3 по вектору a. Тогда на основании ИС (2.5), первого уравнения из (2.3) имеем [6, 7]

(2.9)
${\mathbf{\gamma }} = a_{0}^{'}{\text{sin}}\varphi \cdot {\mathbf{i}}_{1}^{{}} + a_{0}^{'}\cos \varphi \cdot {\mathbf{i}}_{2}^{{}} + a_{{\text{0}}}^{{}}{\mathbf{i}}_{3}^{{}}\quad ({\mathbf{i}}_{3}^{{}} = {\mathbf{a}})$

На основании (2.6), (2.9) запишем компоненты $\omega _{1}^{{}}{\text{,}}\;\omega _{2}^{{}},\;\omega _{3}^{{}}$ вектора ${\mathbf{\omega }}$:

(2.10)
$\omega _{1}^{{}} = a_{0}^{'}\dot {\psi }{\text{sin}}\varphi ,\quad \omega _{2}^{{}} = a_{0}^{'}\dot {\psi }{\text{cos}}\varphi ,\quad \omega _{3}^{{}} = \dot {\varphi } + a_{{\text{0}}}^{{}}{{\dot {\psi }}}$

На рис. 1 приведена геометрическая трактовка прецессий тела относительно вектора ${\mathbf{\gamma }}$ ($O\xi \eta \zeta $ – неподвижная система координат).

Рис. 1.

Геометрическая трактовка прецессий твердого тела.

Замечание 1. При описании кинематических свойств в виде соотношений (2.5)–(2.10) использован метод [7], который отличается от методов, применяемых в [1013, 2123].

Замечание 2. Уравнения (2.2), (2.3) имеют интеграл энергии

(2.11)
$A{\mathbf{\omega }} \cdot {\mathbf{\omega }} - 2({\mathbf{s}} \cdot {\mathbf{\gamma }} + {\mathbf{r}} \cdot {\mathbf{\gamma }}_{{}}^{{(1)}} + {\mathbf{p}} \cdot {\mathbf{\gamma }}_{{}}^{{(2)}}) = {\text{2}}E$
где $E$ – постоянная. Как показано в [6, 7], нахождение условий существования прецессий в задачах динамики твердого тела на основании (2.11) значительно упрощается.

3. Преобразование уравнения (2.2) на ИС (2.5). Внесем в уравнение (2.2) значение ${\mathbf{\omega }}$ из (2.6) и рассмотрим полученное уравнение в базисе ${\mathbf{a}},{\mathbf{\gamma }}{\text{,}}{\mathbf{a}} \times {\mathbf{\gamma }}$ с учетом (2.7), (2.8):

(3.1)
$\begin{gathered} \ddot {\varphi }(A{\mathbf{a}} \cdot {\mathbf{a}}) + \ddot {\psi }(A{\mathbf{a}} \cdot {\mathbf{\gamma }}) - \dot {\psi }_{{}}^{2}[{\mathbf{a}} \cdot (A{\mathbf{\gamma }} \times {\mathbf{\gamma }}){\text{]}} - {\text{[}}{\mathbf{a}} \cdot ({\mathbf{s}} \times {\mathbf{\gamma }})] - \\ \, - b_{0}^{{}}\sin (\psi + \psi _{0}^{{}}) \cdot \{ a_{{\text{0}}}^{{}}[{\mathbf{a}} \cdot ({\mathbf{r}} \times {\mathbf{\gamma }}) - {\mathbf{a}} \cdot {\mathbf{p}}] + {\mathbf{p}} \cdot {\mathbf{\gamma }}\} - \\ \, - {{b}_{0}}\cos (\psi + {{\psi }_{0}}) \cdot \{ {\mathbf{r}} \cdot {\mathbf{\gamma }} - {{a}_{0}}[({\mathbf{r}} \cdot {\mathbf{a}}) + {\mathbf{a}} \cdot ({\mathbf{p}} \times {\mathbf{\gamma }})]\} = 0 \\ \end{gathered} $
(3.2)
$\begin{gathered} \ddot {\varphi }(A{\mathbf{a}} \cdot {\mathbf{\gamma }}) + \ddot {\psi }(A{\mathbf{\gamma }} \cdot {\mathbf{\gamma }}) + 2\dot {\varphi }\dot {\psi }[{\mathbf{a}} \cdot (A{\mathbf{\gamma }} \times {\mathbf{\gamma }}){\text{]}} + \dot {\varphi }_{{}}^{2}{\text{[}}{\mathbf{\gamma }} \cdot ({\mathbf{a}} \times A{\mathbf{a}})] - \\ \, - b_{0}^{{}}\sin (\psi + \psi _{0}^{{}}) \cdot \{ a_{{\text{0}}}^{{}}({\mathbf{p}} \cdot {\mathbf{\gamma }}) + [{\mathbf{a}} \cdot ({\mathbf{r}} \times {\mathbf{\gamma }}) - ({\mathbf{a}} \cdot {\mathbf{p}})]\} - \\ \, - b_{0}^{{}}\cos (\psi + \psi _{0}^{{}}) \cdot \{ a_{{\text{0}}}^{{}}({\mathbf{r}} \cdot {\mathbf{\gamma }}) - ({\mathbf{a}} \cdot {\mathbf{r}}) + [{\mathbf{a}} \cdot ({\mathbf{\gamma }} \times {\mathbf{p}})]\} = 0 \\ \end{gathered} $
(3.3)
$\begin{gathered} \ddot {\varphi }[A{\mathbf{a}} \cdot ({\mathbf{\gamma }} \times {\mathbf{a}})] + \ddot {\psi }[A{\mathbf{\gamma }} \cdot ({\mathbf{a}} \times {\mathbf{\gamma }})] + \dot {\varphi }\dot {\psi }[2(A{\mathbf{\gamma }} \cdot {\mathbf{\gamma }}) - a_{0}^{{'2}}{\text{Sp}}(A) - 2a_{{\text{0}}}^{{}}(A{\mathbf{a}} \cdot {\mathbf{\gamma }})] + \\ \, + \dot {\varphi }_{{}}^{2}[(A{\mathbf{a}} \cdot {\mathbf{\gamma }}) - a_{{\text{0}}}^{{}}(A{\mathbf{a}} \cdot {\mathbf{a}})] + \dot {\psi }_{{}}^{2}[a_{{\text{0}}}^{{}}{\text{(}}A{\mathbf{\gamma }} \cdot {\mathbf{\gamma }}) - (A{\mathbf{a}} \cdot {\mathbf{\gamma }})] - ({\mathbf{a}} \cdot {\mathbf{s}}) + a_{{\text{0}}}^{{}}({\mathbf{s}} \cdot {\mathbf{\gamma }}) - \\ \, - a_{0}^{'}[({\mathbf{p}} \cdot {\mathbf{\gamma }})\cos (\psi + \psi _{0}^{{}}) - ({\mathbf{r}} \cdot {\mathbf{\gamma }})\sin (\psi + \psi _{0}^{{}})] = 0 \\ \end{gathered} $
где ${\text{Sp}}(A) = A_{{11}}^{{}} + A_{{22}}^{{}} + A_{{33}}^{{}}$ – след матрицы A.

По аналогии с (3.1)–(3.3) распишем интеграл (2.11) на ИС (2.5), (2.6):

(3.4)
$\begin{gathered} (A{\mathbf{a}} \cdot {\mathbf{a}})\dot {\varphi }_{{}}^{2} + 2(A{\mathbf{a}} \cdot {\mathbf{\gamma }})\dot {\varphi }\dot {\psi } + (A{\mathbf{\gamma }} \cdot {\mathbf{\gamma }})\dot {\psi }_{{}}^{2} - 2{\text{\{ (}}{\mathbf{s}} \cdot {\mathbf{\gamma }}) + \\ \, + b_{{\text{0}}}^{{}}[\sin (\psi + \psi _{0}^{{}}) \cdot (a_{{\text{0}}}^{{}}{\text{(}}{\mathbf{r}} \cdot {\mathbf{\gamma }}) - ({\mathbf{r}} \cdot {\mathbf{a}}) - {\mathbf{p}} \cdot ({\mathbf{\gamma }} \times {\mathbf{a}})) + \\ \, + \cos (\psi + \psi _{0}^{{}}) \cdot ({\mathbf{a}} \cdot {\mathbf{p}} - a_{{\text{0}}}^{{}}{\text{(}}{\mathbf{p}} \cdot {\mathbf{\gamma }}) + {\mathbf{r}} \cdot ({\mathbf{a}} \times {\mathbf{\gamma }}))]\} = 2E \\ \end{gathered} $
Введем обозначения

$f_{0}^{{}}(\varphi ) = a_{0}^{'}(s_{1}^{{}}\sin \varphi + s_{2}^{{}}\cos \varphi ) + a_{{\text{0}}}^{{}}s_{3}^{{}}$
$\tilde {f}_{0}^{{}}(\varphi ) = a_{0}^{'}(s_{2}^{{}}\sin \varphi - s_{1}^{{}}\cos \varphi )$
$f_{1}^{{}}(\varphi ) = a_{0}^{'}{\text{[}}(a_{{\text{0}}}^{{}}r_{1}^{{}} + p_{2}^{{}})\sin \varphi + (a_{{\text{0}}}^{{}}r_{2}^{{}} - p_{1}^{{}})\cos \varphi - a_{0}^{'}r_{3}^{{}}]$
$f_{2}^{{}}(\varphi ) = a_{0}^{'}{\text{[}}(r_{2}^{{}} - a_{{\text{0}}}^{{}}p_{1}^{{}})\sin \varphi - (a_{{\text{0}}}^{{}}p_{2}^{{}} + r_{1}^{{}})\cos \varphi + a_{0}^{'}p_{3}^{{}}]$
(3.5)
$f_{3}^{{}}(\varphi ) = a_{0}^{'}{\text{[}}(p_{1}^{{}} - a_{{\text{0}}}^{{}}r_{2}^{{}})\sin \varphi + (p_{2}^{{}} + a_{{\text{0}}}^{{}}r_{1}^{{}})\cos \varphi ]$
$f_{4}^{{}}(\varphi ) = a_{0}^{'}{\text{[}}(r_{1}^{{}} + a_{{\text{0}}}^{{}}p_{2}^{{}})\sin \varphi + (r_{2}^{{}} - a_{{\text{0}}}^{{}}p_{1}^{{}})\cos \varphi ]$
$f_{5}^{{}}(\varphi ) = a_{0}^{'}{\text{[}}a_{{\text{0}}}^{{}}(s_{1}^{{}}\sin \varphi + s_{2}^{{}}\cos \varphi ) - a_{0}^{'}s_{3}^{{}}]$
$f_{6}^{{}}(\varphi ) = - a_{0}^{'}{\text{[}}a_{0}^{'}(r_{1}^{{}}\sin \varphi + r_{2}^{{}}\cos \varphi ) + a_{{\text{0}}}^{{}}r_{3}^{{}}]$
$f_{7}^{{}}(\varphi ) = a_{0}^{'}{\text{[}}a_{0}^{'}(p_{1}^{{}}\sin \varphi + p_{2}^{{}}\cos \varphi ) + a_{{\text{0}}}^{{}}p_{3}^{{}}]$

Сначала запишем интеграл (3.4), в силу (3.5):

(3.6)
$\begin{gathered} (A{\mathbf{a}} \cdot {\mathbf{a}})\dot {\varphi }_{{}}^{2} + 2(A{\mathbf{a}} \cdot {\mathbf{\gamma }})\dot {\varphi }\dot {\psi } + (A{\mathbf{\gamma }} \cdot {\mathbf{\gamma }})\dot {\psi }_{{}}^{2} - \\ \, - 2{\text{[}}f_{0}^{{}}(\varphi ) + b_{0}^{{}}(f_{1}^{{}}(\varphi )\sin (\psi + \psi _{0}^{{}}) + f_{2}^{{}}(\varphi )\cos (\psi + \psi _{0}^{{}}))] = 2E \\ \end{gathered} $

Затем обратимся к уравнениям (3.1)(3.3). На основании (3.5) имеем

(3.7)
$\begin{gathered} \ddot {\varphi }(A{\mathbf{a}} \cdot {\mathbf{a}}) + \ddot {\psi }(A{\mathbf{a}} \cdot {\mathbf{\gamma }}) - \dot {\psi }_{{}}^{2}[{\mathbf{a}} \cdot (A{\mathbf{\gamma }} \times {\mathbf{\gamma }}){\text{]}} + \tilde {f}_{0}^{{}}(\varphi ) - \\ \, - b_{0}^{{}}(f_{3}^{{}}(\varphi )\sin (\psi + \psi _{0}^{{}}) + f_{4}^{{}}(\varphi )\cos (\psi + \psi _{0}^{{}})) = 0 \\ \end{gathered} $
(3.8)
$\begin{gathered} \ddot {\varphi }(A{\mathbf{a}} \cdot {\mathbf{\gamma }}) + \ddot {\psi }(A{\mathbf{\gamma }} \cdot {\mathbf{\gamma }}) - 2\dot {\varphi }\dot {\psi }{\text{[}}{\mathbf{a}} \cdot ({\mathbf{\gamma }} \times A{\mathbf{\gamma }})] - \dot {\varphi }_{{}}^{2}[{\mathbf{a}} \cdot ({\mathbf{\gamma }} \times A{\mathbf{a}})] - \\ \, - b_{0}^{{}}[f_{1}^{{}}(\varphi )\cos (\psi + \psi _{0}^{{}}) - f_{2}^{{}}(\varphi )\sin (\psi + \psi _{0}^{{}}{\text{)]}} = 0 \\ \end{gathered} $
(3.9)

Представление соотношений (3.1)–(3.4) в виде (3.6)–(3.9) связано с решением задачи о замене одного из уравнений (3.7)–(3.9) интегралом (3.6). Вычислим производную по времени от левой части уравнения (3.6), используя соотношения

${\mathbf{\dot {\gamma }}} = \dot {\varphi }({\mathbf{\gamma }} \times {\mathbf{a}}){\text{,}}\quad \dot {f}_{0}^{{}} = - \varphi \tilde {f}_{0}^{{}}(\varphi ),\quad \dot {f}_{1}^{{}}(\varphi ) = \dot {\varphi }f_{3}^{{}}(\varphi ),\quad \dot {f}_{2}^{{}}(\varphi ) = \dot {\varphi }f_{4}^{{}}(\varphi )$

Если учесть в полученном уравнении соотношения (3.7), (3.8), то имеем тождество. То есть интеграл энергии можно рассматривать только вместо одного из уравнений (3.7), (3.8). Отметим, что данное свойство выполняется, если ни одно из уравнений (3.7), (3.8) не становится тождеством.

4. О регулярных прецессиях тела. Изучим случай регулярной прецессии тела: $\dot {\varphi } = n$, $\dot {\psi } = m$, где $\dot {\varphi } = n$ и m – постоянные. Он рассмотрен частично в [2123], но методом, который отличается от предложенного в данной статье. Запишем уравнения (3.6), (3.8) в виде

(4.1)
$\begin{gathered} f_{1}^{{}}(\varphi )\sin (\psi + \psi _{0}^{{}}) + f_{2}^{{}}(\varphi )\cos (\psi + \psi _{0}^{{}}) = a_{0}^{{'2}}m_{{}}^{2}H_{{\text{2}}}^{{}}(\varphi ) \\ f_{2}^{{}}(\varphi )\sin (\psi + \psi _{0}^{{}}) - f_{1}^{{}}(\varphi )\cos (\psi + \psi _{0}^{{}}) = 2nma_{0}^{'}G_{{\text{2}}}^{{}}(\varphi ) \\ \end{gathered} $
(4.2)
$\begin{gathered} H_{{\text{2}}}^{{}}(\varphi ) = \frac{1}{2}(A_{{11}}^{{}} - A_{{22}}^{{}})\cos 2\varphi + A_{{12}}^{{}}\sin 2\varphi + ..., \\ G_{{\text{2}}}^{{}}(\varphi ) = \frac{1}{2}(A_{{11}}^{{}} - A_{{22}}^{{}})\sin 2\varphi - A_{{12}}^{{}}\cos 2\varphi + ... \\ \end{gathered} $

В (4.1) не произведена замена $\psi $ на mt, поскольку данный прием может быть использован и в других случаях прецессий. Очевидно, что в силу (3.5) $f_{1}^{2}(\varphi ) + f_{2}^{2}(\varphi ) \ne 0$. Найдем из (4.1) $\sin (\psi + \psi _{0}^{{}})$, $\cos (\psi + \psi _{0}^{{}})$ и подставим их в равенство sin2(ψ + ψ0) + + $\cos _{{}}^{2}(\psi + \psi _{0}^{{}})$ = 1. Требуя, чтобы полученное уравнение было тождеством по $\varphi $, получим

(4.3)
$(m_{{}}^{2} - 4n_{{}}^{2})[(A_{{11}}^{{}} - A_{{22}}^{{}})_{{}}^{2} + 4A_{{12}}^{2}] = 0,\quad (m_{{}}^{2} - 4n_{{}}^{2})A_{{12}}^{{}}(A_{{11}}^{{}} - A_{{22}}^{{}}) = 0$

Далее будем изучать только случай m = 2n, так как другие варианты исследованы в [21, 22]. Обратимся к уравнению (3.7), полагая $\psi = 2nt$, $\varphi = nt$. Функция $(A{\mathbf{\gamma }} \cdot {\mathbf{\gamma }})$ при значении вектора ${\mathbf{\gamma }}$ из (2.7) имеет в общем случае вид

(4.4)
$\begin{gathered} (A{\mathbf{\gamma }} \cdot {\mathbf{\gamma }}) = \frac{1}{2}a_{0}^{{'2}}(A_{{22}}^{{}} - A_{{11}}^{{}})\cos 2\varphi + a_{0}^{{'2}}A_{{12}}^{{}}\sin 2\varphi + \\ \, + 2a_{0}^{'}a_{{\text{0}}}^{{}}(A_{{13}}^{{}}\sin \varphi + A_{{23}}^{{}}\cos \varphi ) + \frac{1}{2}{\text{[}}a_{0}^{{'2}}(A_{{11}}^{{}} + A_{{22}}^{{}}) + 2A_{{33}}^{{}}a_{{\text{0}}}^{2}{\text{]}} \\ \end{gathered} $

То есть функция (4.4) содержит тригонометрические функции $\cos 2\varphi $, $\sin 2\varphi $, $\cos \varphi $, $\sin \varphi $. Очевидно, что функция $f_{1}^{{}}(\varphi )\sin (2\varphi + \psi _{0}^{{}}) + f_{2}^{{}}(\varphi )\cos (2\varphi + \psi _{0}^{{}})$ не должна содержать слагаемые, содержащие $\sin 3\varphi $, $\cos 3\varphi $. Для того, чтобы указанная функция не содержала $\sin 3\varphi $, $\cos 3\varphi $, должны, в силу (3.5), выполняться условия

(4.5)
$p_{1}^{{}} = r_{2}^{{}},\quad p_{2}^{{}} = - {{r}_{1}}$

Запишем функции $f_{1}^{{}}(\varphi )$, $f_{2}^{{}}(\varphi )$ из (3.5) при условиях (4.5)

(4.6)
$\begin{gathered} f_{1}^{{}}(\varphi ) = (1 - a_{{\text{0}}}^{{}})[ - a_{0}^{'}(r_{1}^{{}}\sin \varphi + r_{2}^{{}}\cos \varphi ) - (1 + a_{{\text{0}}}^{{}})r_{3}^{{}}] \\ f_{2}^{{}}(\varphi ) = (1 - a_{{\text{0}}}^{{}})[a_{0}^{'}(r_{2}^{{}}\sin \varphi - r_{1}^{{}}\cos \varphi ) + (1 + a_{{\text{0}}}^{{}})p_{3}^{{}}] \\ \end{gathered} $

Подставим функции (4.4), (4.6) в уравнения (3.6) и учтем равенства $\psi = 2\varphi $, (4.6). Для того, чтобы полученное уравнение было тождеством по $\varphi $, необходимо положить

(4.7)
$A_{{12}}^{{}} = 0,\quad A_{{22}}^{{}} = A_{{11}}^{{}},\quad p_{3}^{{}} = 0,\quad r_{3}^{{}} = 0$

Таким образом, в силу (4.5), (4.7) для векторов ${\mathbf{p}}$ и ${\mathbf{r}}$ имеем векторные соотношения

(4.8)
${\mathbf{p}} = (r_{2}^{{}}{\text{,}} - r_{1}^{{}},0),\quad {\mathbf{r}} = (r_{1}^{{}}{\text{,}}r_{2}^{{}},0)$
то есть в силу (4.8) ${\mathbf{p}} \bot {\mathbf{r}}$, ${\mathbf{p}} \bot {\mathbf{a}}$, ${\mathbf{r}} \bot {\mathbf{a}}$ – ортогональные вектора p и r лежат в плоскости Oxy. Запишем дополнительные условия, при выполнении которых уравнение (3.6) становится тождеством по φ:

(4.9)
$\begin{gathered} 2a_{0}^{'}n_{{}}^{2}A_{{13}}^{{}}(1 + 2a_{{\text{0}}}^{{}}) - a_{0}^{'}s_{1}^{{}} - (1 - a_{{\text{0}}}^{{}})(r_{{\text{1}}}^{{}}\sin \psi _{{\text{0}}}^{{}} - r_{{\text{2}}}^{{}}\cos \psi _{{\text{0}}}^{{}}) = 0 \\ 2a_{0}^{'}n_{{}}^{2}A_{{23}}^{{}}(1 + 2a_{{\text{0}}}^{{}}) - a_{0}^{'}s_{2}^{{}} + (1 - a_{{\text{0}}}^{{}})(r_{{\text{1}}}^{{}}\cos \psi _{{\text{0}}}^{{}} + r_{{\text{2}}}^{{}}\sin \psi _{{\text{0}}}^{{}}) = 0 \\ \end{gathered} $

На основании условий (4.5)–(4.7) уравнение (3.8) будет тождеством по $t$, если имеют место равенства

(4.10)
$\begin{gathered} a_{0}^{'}n_{{}}^{2}A_{{13}}^{{}}(4a_{{\text{0}}}^{{}} + 1) - (1 - a_{{\text{0}}}^{{}})(r_{{\text{1}}}^{{}}\sin \psi _{{\text{0}}}^{{}} - r_{{\text{2}}}^{{}}\cos \psi _{{\text{0}}}^{{}}) = 0 \\ a_{0}^{'}n_{{}}^{2}A_{{23}}^{{}}(4a_{{\text{0}}}^{{}} + 1) + (1 - a_{{\text{0}}}^{{}})(r_{{\text{1}}}^{{}}{\text{cos}}\psi _{{\text{0}}}^{{}} + r_{{\text{2}}}^{{}}\sin \psi _{{\text{0}}}^{{}}) = 0 \\ \end{gathered} $

Учитывая в уравнениях (4.9) равенства (4.10), получим

(4.11)
$s_{1}^{{}} = n_{{}}^{2}A_{{13}}^{{}},\quad s_{2}^{{}} = n_{{}}^{2}A_{{23}}^{{}}$

Рассмотрим уравнение (3.9). Потребуем, чтобы оно было тождеством по $t$ при условиях (4.5), (4.7), (4.10), (4.11). Тогда получим

(4.12)
$a_{{\text{0}}}^{{}} = \frac{1}{6},\quad s_{3}^{{}} = \frac{{2n_{{}}^{2}}}{3}{\text{[}}3A_{{22}}^{{}} + 5(A_{{33}}^{{}} - A_{{11}}^{{}})]$

Итак, условием существования регулярных прецессий ($\psi = 2nt$, $\varphi = nt$) тела относительно вектора ${\mathbf{\gamma }}$ являются равенства (4.5), (4.7), (4.10)–(4.12). Первые два условия из (4.7) означают, что вектор a перпендикулярен круговому сечению эллипсоида инерции тела. Уравнения (4.10) служат для определения постоянной $\psi _{{\text{0}}}^{{}}$:

(4.13)
$\psi _{{\text{0}}}^{{}} = {\text{arcsin}}\frac{{a_{0}^{'}n_{{}}^{2}(4a_{{\text{0}}}^{{}} + 1)(r_{{\text{1}}}^{{}}A_{{13}}^{{}} - r_{{\text{2}}}^{{}}A_{{23}}^{{}})}}{{(1 - a_{{\text{0}}}^{{}})(r_{{\text{1}}}^{2} + r_{{\text{2}}}^{2})}}$
а также постоянной $n$:

(4.14)
${{n}^{4}} = \frac{{(1 - {{a}_{{\text{0}}}})(r_{{\text{1}}}^{2} + r_{{\text{2}}}^{2})}}{{(1 + {{a}_{{\text{0}}}}){{{(1 + 4{{a}_{{\text{0}}}})}}^{2}}(A_{{13}}^{2} + A_{{23}}^{2})}}$

В соотношениях (4.13), (4.14) необходимо учесть значение a0 = 1/6.

5. Прецессионно-изоконические движения. Данные движения можно описать инвариантным соотношением [6]

(5.1)
$\psi = \varphi $

В силу (5.1) запишем подвижный и неподвижный годографы в виде

(5.2)
${\mathbf{\omega }}_{{\text{M}}}^{{}} = \dot {\varphi }{\text{[}}a_{0}^{'}\sin \varphi \cdot {\mathbf{i}}_{1}^{{}} + a_{0}^{'}\cos \varphi \cdot {\mathbf{i}}_{2}^{{}} + (1 + a_{{\text{0}}}^{{}}) \cdot {\mathbf{i}}_{3}^{{}}{\text{]}}$
(5.3)
${\mathbf{\omega }}_{{\text{F}}}^{{}} = \dot {\varphi }{\text{[}}a_{0}^{'}\sin (\varphi + \psi _{{\text{0}}}^{{}}) \cdot {\mathbf{\gamma }}_{1}^{{}} - a_{0}^{'}\cos (\varphi + \psi _{{\text{0}}}^{{}}) \cdot {\mathbf{\gamma }}_{2}^{{}} + (1 + a_{{\text{0}}}^{{}}) \cdot {\mathbf{\gamma }}{\text{]}}$

Из векторных равенств (5.2), (5.3) следует, что годографы (5.2), (5.3) симметричны друг другу относительно касательной плоскости к соответствующим аксоидам. То есть движение тела характеризуется не только прецессией относительно ${\mathbf{\gamma }}$, но и изоконичностью годографов.

Исследование прецессионно-изоконических движений проведен в случае, когда эллипсоид инерции является сферой:

(5.4)
$A = {\text{diag(}}{{A}_{1}},{{A}_{1}},{{A}_{1}}{\text{)}}$

На основании условия (5.4) систему уравнений (3.6), (3.7), (3.9) представим так

(5.5)
$A_{1}^{{}}(1 + a_{{\text{0}}}^{{}})\dot {\varphi }_{{}}^{2} = E + f_{{\text{0}}}^{{}}(\varphi ) + b_{{\text{0}}}^{{}}[{{f}_{1}}(\varphi )\sin (\varphi + \psi _{{\text{0}}}^{{}}) + f_{{\text{2}}}^{{}}(\varphi )\cos (\varphi + \psi _{{\text{0}}}^{{}}){\text{]}}$
(5.6)
$A_{1}^{{}}(1 + a_{{\text{0}}}^{{}})\ddot {\varphi } = - \tilde {f}_{{\text{0}}}^{{}}(\varphi ) + b_{{\text{0}}}^{{}}[f_{{\text{3}}}^{{}}(\varphi )\sin (\varphi + \psi _{{\text{0}}}^{{}}) + f_{{\text{4}}}^{{}}(\varphi )\cos (\varphi + \psi _{{\text{0}}}^{{}}){\text{]}}$
(5.7)
$A_{1}^{{}}a_{0}^{{'2}}\dot {\varphi }_{{}}^{2} + f_{{\text{5}}}^{{}}(\varphi ) + f_{{\text{6}}}^{{}}(\varphi )\sin (\varphi + \psi _{{\text{0}}}^{{}}) + f_{{\text{7}}}^{{}}(\varphi )\cos (\varphi + \psi _{{\text{0}}}^{{}}) = {\text{0}}$

Продифференцируем левую и правую части уравнения (5.5) по $t$ и подставим в полученное уравнение $\ddot {\varphi }$ из (5.6), $f_{{\text{1}}}^{{}}(\varphi )$ и $f_{{\text{2}}}^{{}}(\varphi )$ из (3.5). Так как это уравнение должно быть тождеством по $\varphi $, то находим следующие условия на параметры задачи:

(5.8)
$s_{1}^{{}} = - r_{3}^{{}}\cos \psi _{{\text{0}}}^{{}} - p_{3}^{{}}\sin \psi _{{\text{0}}}^{{}},\quad s_{2}^{{}} = - r_{3}^{{}}\sin \psi _{{\text{0}}}^{{}} + p_{3}^{{}}\cos \psi _{{\text{0}}}^{{}}$
(5.9)
$p_{1}^{{}} = r_{2}^{{}},\quad p_{2}^{{}} = - r_{1}^{{}},\quad r_{2}^{{}}\cos \psi _{{\text{0}}}^{{}} - r_{1}^{{}}\sin \psi _{{\text{0}}}^{{}} = 0$

Исключим из уравнений (5.5), (5.7) величину $\dot {\varphi }_{{}}^{2}$. В результате получим уравнение, из которого следуют равенства

(5.10)
$\begin{gathered} (2a_{{\text{0}}}^{{}} - 1)(s_{1}^{{}} - r_{3}^{{}}\cos \psi _{{\text{0}}}^{{}} - p_{3}^{{}}\sin \psi _{{\text{0}}}^{{}}) = 0{\text{,}} \\ (2a_{{\text{0}}}^{{}} - 1)(s_{2}^{{}} - r_{3}^{{}}\sin \psi _{{\text{0}}}^{{}} + p_{3}^{{}}\cos \psi _{{\text{0}}}^{{}}) = 0 \\ \end{gathered} $
(5.11)
$(1 + a_{{\text{0}}}^{{}})E + (1 - a_{{\text{0}}}^{{}})(1 + 2a_{{\text{0}}}^{{}})s_{3}^{{}} - 2a_{{\text{0}}}^{{}}(1 - a_{{\text{0}}}^{{}})(r_{1}^{{}}\cos \psi _{{\text{0}}}^{{}} + r_{2}^{{}}\sin \psi _{{\text{0}}}^{{}}) = 0$

Если в уравнениях (5.10) положить $2a_{{\text{0}}}^{{}} - 1 \ne 0$, то получим два равенства

(5.12)
$s_{1}^{{}} = r_{3}^{{}}\cos \psi _{{\text{0}}}^{{}} + p_{3}^{{}}\sin \psi _{{\text{0}}}^{{}},\quad s_{2}^{{}} = r_{3}^{{}}\sin \psi _{{\text{0}}}^{{}} - p_{3}^{{}}\cos \psi _{{\text{0}}}^{{}}$

Приравняем величины s1 и s2; используя (5.8), (5.12), имеем нулевые значения для p3, r3, а также для s1, s2. На основании уравнения (5.5) и функций (4.6) находим, что $\dot {\varphi } = {\text{const}}$. Данный вариант, в силу (5.1), характеризует регулярную прецессию, что исключено из рассмотрения. Поэтому в (5.10) необходимо положить

(5.13)
$a_{{\text{0}}}^{{}} = \frac{1}{2},\quad \left( {\theta _{{\text{0}}}^{{}} = \frac{\pi }{3}} \right)$

Рассмотрим геометрическую интерпретацию условий (5.9). Введем вспомогательные векторы

(5.14)
${\mathbf{p}}_{*}^{{}} = r_{2}^{{}}{\mathbf{i}}_{1}^{{}} - r_{1}^{{}}{\mathbf{i}}_{2}^{{}},\quad {\mathbf{r}}_{*}^{{}} = r_{1}^{{}}{\mathbf{i}}_{1}^{{}} + r_{2}^{{}}{\mathbf{i}}_{2}^{{}}$

Тогда из (5.9), (5.14) следует, что имеют место соотношения

(5.15)
${\text{tg}}\psi _{{\text{0}}}^{{}} = \frac{{r_{2}^{{}}}}{{r_{1}^{{}}}},\quad {\mathbf{p}}_{*}^{{}} \cdot {\mathbf{r}}_{ + }^{{}} = 0$

Из равенств (5.8) в силу третьей формулы системы (5.9) получим

(5.16)
$s_{1}^{{}} = - \frac{{r_{1}^{{}}p_{3}^{{}} + r_{{\text{1}}}^{{}}r_{{\text{3}}}^{{}}}}{{\sqrt {r_{{\text{1}}}^{2} + r_{{\text{2}}}^{2}} }},\quad s_{2}^{{}} = - \frac{{r_{1}^{{}}p_{3}^{{}} - r_{{\text{1}}}^{{}}r_{{\text{3}}}^{{}}}}{{\sqrt {r_{{\text{1}}}^{2} + r_{{\text{2}}}^{2}} }}$

Равенства (5.16) позволяют связать параметры s1, s2 и $p_{3}^{{}}{\text{,}}r_{3}^{{}}$:

(5.17)
$s_{1}^{2} + s_{2}^{2} = r_{3}^{2} + p_{3}^{2}$

Преобразуем равенство (5.5) с учетом условий (5.8), (5.9), (5.16), (5.17) и значения a0 = 1/2:

(5.18)
$\dot {\varphi } = \sqrt {h_{0}^{{}} + h_{1}^{{}}\cos (\varphi + \beta _{{\text{0}}}^{{}})} $
(5.19)
$h_{0}^{{}} = - \frac{1}{{9A_{1}^{{}}}}(s_{3}^{{}} + 2\sqrt {r_{{\text{1}}}^{2} + r_{{\text{2}}}^{2}} ),\quad h_{1}^{{}} = \frac{{2\sqrt 3 }}{{3A_{1}^{{}}}}\sqrt {r_{{\text{3}}}^{2} + p_{{\text{3}}}^{2}} $
(5.20)
$\sin \beta _{{\text{0}}}^{{}} = - \frac{{s_{1}^{{}}}}{{\sqrt {p_{{\text{3}}}^{2} + r_{{\text{3}}}^{2}} }},\quad \cos \beta _{{\text{0}}}^{{}} = - \frac{{s_{2}^{{}}}}{{\sqrt {p_{{\text{3}}}^{2} + r_{{\text{3}}}^{2}} }}$

Выполним в (5.18) замену $\varphi + \beta _{{\text{0}}}^{{}} = \tilde {\varphi }$. Тогда из (5.18) имеем

(5.21)
$\dot {\tilde {\varphi }} = \sqrt {h_{0}^{{}} + h_{1}^{{}}} \sqrt {1 - k_{{}}^{2}{{{\sin }}^{2}}\frac{{\tilde {\varphi }}}{2}} \quad \left( {k_{{}}^{2} = \frac{{2h_{1}^{{}}}}{{h_{0}^{{}} + h_{1}^{{}}}}} \right)$

Положим в (5.21) $\tilde {\varphi } = 2\varepsilon $:

$\dot {\varepsilon } = \frac{{\sqrt {h_{0}^{{}} + h_{1}^{{}}} }}{2}\sqrt {1 - k_{{}}^{2}{{{\sin }}^{2}}\varepsilon } $
или

(5.22)
$\int\limits_0^{\varepsilon _{1}^{{}}} {\frac{{d\,\varepsilon }}{{\sqrt {1 - k_{{}}^{2}{{{\sin }}^{2}}\varepsilon } }}} = \frac{{\sqrt {h_{0}^{{}} + h_{1}^{{}}} }}{2}t = \tau $

В силу выполненных преобразований параметры задачи должны удовлетворять условию $h_{0}^{{}} > {\text{0}}$ ($h_{1}^{{}} > {\text{0}}$ для всех значений ${{r}_{3}}{\text{,}}p_{3}^{{}}$, так как $A_{1}^{{}} > {\text{0}}$), где h0 имеет вид, указанный в (5.19). Параметр $\beta _{{\text{0}}}^{{}}$ определим из (5.20):

(5.23)
$\beta _{{\text{0}}}^{{}} = {\text{arctg}}\frac{{s_{1}^{{}}}}{{s_{2}^{{}}}}$

Из интегрального соотношения (5.22) имеем

(5.24)
$\varepsilon (\tau ) = {\text{am}}\tau ,\quad \sin \varepsilon (\tau ) = {\text{sn}}\tau $
где ${\text{am}}\tau $, ${\text{sn}}\tau $ – эллиптические функции от τ. На основании указанных выше вычислений найдем

(5.25)
$\varphi (\tau ) = 2\varepsilon (\tau ) - \beta _{{\text{0}}}^{{}},\quad \sin \varphi (\tau ) = {\text{sin}}2\varepsilon (\tau )\cos \beta _{{\text{0}}}^{{}} - \cos 2\varepsilon (\tau ){\text{sin}}\beta _{{\text{0}}}^{{}}$

Учтем во второй формуле из (5.25) значения (5.23), (5.24). Тогда

(5.26)
$\sin \varphi (\tau ) = \frac{1}{{s_{1}^{2} + s_{2}^{2}}}{\text{[}}s_{1}^{{}}({\text{cn}}_{{}}^{{\text{2}}}\tau - {\text{sn}}_{{}}^{{\text{2}}}\tau ) - 2s_{2}^{{}}{\text{sn}}\tau {\text{cn}}\tau ]$
где в силу (5.22) $\tau = \frac{{h_{2}^{{}}}}{2}t$ $(h_{2}^{{}} = \sqrt {h_{0}^{{}} + h_{1}^{{}}} )$. Значение $\dot {\varphi }$ определим из (5.21), (5.22)
(5.27)
$\dot {\varphi } = h_{2}^{{}}{\text{dn}}\frac{{h_{2}^{{}}}}{2}t$
а для получения cosφ(τ) необходимо использовать формулу (5.26), подставив в нее $\tau = \frac{{h_{2}^{{}}}}{2}t$. Тогда функция $\cos \varphi (t) = \sqrt {1 - \sin _{{}}^{2}\varphi (t)} $. Таким образом, решение уравнений движения сферического тела для случая прецессионно-изоконических движений построено: оно выражается формулами (2.7)(2.10), в которых необходимо учесть условия (5.8), (5.9), (5.13), (5.27). Оно выражается через эллиптические функции ${\text{sn(}}\tau )$, ${\text{cn(}}\tau )$, ${\text{dn(}}\tau )$, где $\tau = \frac{{h_{2}^{{}}}}{2}t$. Представляет определенный интерес значение ${{\theta }_{0}} = 60^\circ $, а также условие, которое выражается ортогональностью векторов ${\mathbf{r}}{\kern 1pt} *{\text{,}}\,\,{\mathbf{p}}{\kern 1pt} *$. Отметим, что в случае регулярных прецессионных движений (см. п. 4) ${{\theta }_{0}} = \arccos \left( {{1 \mathord{\left/ {\vphantom {1 6}} \right. \kern-0em} 6}} \right)$.

6. Случай двух однородных силовых полей. Рассмотрим движение тела под действием двух однородных силовых полей: положим ${\mathbf{p}}{\kern 1pt} * = {\mathbf{0}}$ ($p_{{\text{i}}}^{{}} = 0$, $i = \overline {1,3} $). Тогда из условий (5.8), (5.9) найдем

(6.1)
$r_{1}^{{}} = r_{2}^{{}} = 0,\quad s_{1}^{{}} = - r_{3}^{{}}\cos \psi _{{\text{0}}}^{{}},\quad s_{2}^{{}} = - r_{3}^{{}}\sin \psi _{{\text{0}}}^{{}}$

Значения ${{a}_{0}} = {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}$ из (5.13) не изменяется. В силу (6.1) ${\mathbf{r}} = r_{3}^{{}}{\mathbf{a}}$, то есть вектор r коллинеарен вектору a (центр приведения сил второго поля находится на оси, проходящей через a). Из последних равенств (6.1) следует, что $s_{1}^{2} + s_{2}^{2} = r_{3}^{2}$ (формулы (5.16) не имеют места, так как они получены из последней формулы системы (5.9)). Выражение для h0 из (5.19) упрощается: $h_{0}^{{}} = - \frac{{s_{3}^{{}}}}{{9A_{1}^{{}}}}$, а значения h1 из (5.19) и ${\text{sin}}\beta _{{\text{0}}}^{{}}$, $\cos {{\beta }_{0}}$, в силу (6.1), принимают вид

(6.2)
$h_{1}^{{}} = \frac{{2r_{3}^{{}}\sqrt 3 }}{{3A_{1}^{{}}}},\quad {\text{sin}}\beta _{{\text{0}}}^{{}} = \cos \psi _{{\text{0}}}^{{}},\quad \cos \beta _{{\text{0}}}^{{}} = {\text{sin}}\psi _{{\text{0}}}^{{}}$

То есть и в случае (6.1), (6.2) $\varphi ({\text{t}})$, $\sin \varphi (\tau )$, $\cos \varphi (\tau )$ являются эллиптическими функциями времени. Движение тела обладает свойством периодичности (см. формулы (5.2), (5.3), (5.25), (5.26)).

Замечание 3. Из полученных выше результатов (см. пп. 4, 6) в случае, когда на тело действует сила только одного однородного поля (уравнения движения становятся уравнениями Эйлера–Пуассона), прецессионно-изоконические движения тела имеют место только в решении Лагранжа для сферического гироскопа. Данный вывод согласуется с результатами [57]. Отметим, что в решении Лагранжа изоконические движения (без дополнительного свойства прецессионности) возможны только в случае несферического распределения масс.

7. Прецессионно-изоконические движения динамически симметричного твердого тела. Рассмотрим случай, когда твердое тело обладает свойством динамической симметрии относительно вектора ${\mathbf{a}}$. Запишем вначале уравнения (3.6), (3.7) при $A = {\text{diag}}$(A1, A1, A3):

(7.1)
$h_{{\text{0}}}^{{}}(1 + a_{{\text{0}}}^{{}})\dot {\varphi }_{{}}^{2} = 2[E + f_{{\text{0}}}^{{}}(\varphi ) + b_{{\text{0}}}^{{}}(f_{{\text{1}}}^{{}}(\varphi )\sin (\varphi + \psi _{{\text{0}}}^{{}}) + f_{{\text{2}}}^{{}}(\varphi )\cos (\varphi + \psi _{{\text{0}}}^{{}})){\text{]}}$
(7.2)
$A_{3}^{{}}(1 + a_{{\text{0}}}^{{}})\ddot {\varphi } = - \tilde {f}_{{\text{0}}}^{{}}(\varphi ) + b_{{\text{0}}}^{{}}[f_{{\text{3}}}^{{}}(\varphi )\sin (\varphi + \psi _{{\text{0}}}^{{}}) + f_{{\text{4}}}^{{}}(\varphi )\cos (\varphi + \psi _{{\text{0}}}^{{}}){\text{]}}$
(7.3)
$h_{0}^{{}} = (A_{1}^{{}} + A_{3}^{{}}) + a_{0}^{{}}(A_{3}^{{}} - A_{1}^{{}})$

Метод изучения уравнений (7.1), (7.2) сохраним в прежнем виде (см. п. 5). Тогда получим следующие условия на параметры задачи: равенства (5.9) не изменяются, а вместо равенств (5.8) следует рассматривать уравнения

(7.4)
$\mu _{{\text{0}}}^{{}}s_{{\text{1}}}^{{}} = p_{{\text{3}}}^{{}}\sin \psi _{{\text{0}}}^{{}} + r_{{\text{3}}}^{{}}\cos \psi _{{\text{0}}}^{{}},\quad \mu _{{\text{0}}}^{{}}s_{{\text{2}}}^{{}} = p_{{\text{3}}}^{{}}\cos \psi _{{\text{0}}}^{{}} - r_{{\text{3}}}^{{}}\sin \psi _{{\text{0}}}^{{}}$

Здесь $\mu _{{\text{0}}}^{{}} = (A_{{\text{3}}}^{{}} - h_{{\text{0}}}^{{}}){\text{/}}A_{{\text{3}}}^{{}}$, или, в силу (7.3), имеем следующее значение $\mu _{{\text{0}}}^{{}}$ через параметры задачи:

(7.5)
$\mu _{{\text{0}}}^{{}} = - \frac{1}{{A_{{\text{3}}}^{{}}}}[(1 - a_{{\text{0}}}^{{}})A_{{\text{1}}}^{{}} + a_{{\text{0}}}^{{}}{{A}_{3}}]$

Из (7.5) следует, что $\mu _{{\text{0}}}^{{}} \ne 1$, так как в противном случае $h_{0}^{{}} = 0$, что исключаем из рассмотрения на основании значения $a_{{\text{0}}}^{{}}$: ${\text{|}}a_{{\text{0}}}^{{}}{\text{|}} > 1$, которое получаем из равенства $h_{0}^{{}}\, = \,0$.

Запишем значения $\dot {\varphi }_{{}}^{2}$ из (7.2), использовав обозначения (3.5) и условия (5.9), (7.4):

(7.6)
$\begin{gathered} \dot {\varphi }_{{}}^{2} = \frac{2}{{h_{{\text{0}}}^{{}}(1 + a_{{\text{0}}}^{{}})}}{\text{[}}a_{0}^{'}{\text{(1}} - \mu _{{\text{0}}}^{{}}{\text{)(}}s_{{\text{1}}}^{{}}\sin \varphi + s_{{\text{2}}}^{{}}\cos \varphi ) + \\ \, + E + a_{{\text{0}}}^{{}}s_{{\text{3}}}^{{}} - (1 - a_{{\text{0}}}^{{}})(r_{{\text{1}}}^{{}}\cos \psi _{{\text{0}}}^{{}} + r_{{\text{2}}}^{{}}\sin \psi _{{\text{0}}}^{{}}{\text{)]}} \\ \end{gathered} $
а также из уравнения (3.9) (полагаем $a_{0}^{{}}(A_{1}^{{}} - A_{3}^{{}}) - A_{1}^{{}} \ne 0$, так как в противном случае $s_{1}^{{}} = 0$, $s_{2}^{{}} = 0$):

(7.7)
$\begin{gathered} \dot {\varphi }_{{}}^{2} = \frac{1}{{a_{0}^{'}[a_{{\text{0}}}^{{}}(A_{1}^{{}} - A_{3}^{{}}) - A_{3}^{{}}]}}{\text{[}}a_{{\text{0}}}^{{}}{\text{(}}\mu _{{\text{0}}}^{{}} - {\text{1)(}}s_{{\text{1}}}^{{}}\sin \varphi + s_{{\text{2}}}^{{}}\cos \varphi ) + \\ \, + a_{0}^{'}(s_{{\text{3}}}^{{}} + r_{{\text{1}}}^{{}}\cos \psi _{{\text{0}}}^{{}} + r_{{\text{2}}}^{{}}\sin \psi _{{\text{0}}}^{{}}{\text{)]}} \\ \end{gathered} $

Приравнивая $\dot {\varphi }_{{}}^{2}$ из уравнений (7.6), (7.7), найдем уравнение, которое должно быть тождеством по φ. Тогда получим

(7.8)
$F(a) = a_{{\text{0}}}^{{\text{2}}}(A_{3}^{{}} - A_{1}^{{}}) + a_{{\text{0}}}^{{}}(3A_{1}^{{}} + A_{3}^{{}}) - 2A_{3}^{{}} = 0$

Если $A_{3}^{{}} = A_{1}^{{}}$, то из (7.8) следует ${{a}_{0}} = {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}$, что соответствует уже изученному варианту. При $A_{3}^{{}} \ne A_{1}^{{}}$ параметр a0 имеет значение

(7.9)
$a_{{\text{0}}}^{{}} = \frac{1}{{A_{1}^{{}} - A_{3}^{{}}}}{\text{[3}}A_{1}^{{}} + A_{3}^{{}} \pm \sqrt {8(A_{1}^{2} + A_{3}^{2}) + (A_{1}^{{}} - A_{2}^{{}})_{{}}^{2}} ]$

Для исследования промежутка изменения параметра a0 целесообразно не применять формулу (7.9), а использовать графо-аналитический метод исследования функции (7.8), обозначив через $y = {{{{A}_{1}}} \mathord{\left/ {\vphantom {{{{A}_{1}}} {{{A}_{3}}}}} \right. \kern-0em} {{{A}_{3}}}}$, а a0 – через x. Из неравенств треугольника на главные моменты инерции тела ${{A}_{i}}$ $(i = \overline {1,3} )$ следует, что $y > {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}$; из геометрического условия на параметр a0 имеем неравенство ${\text{|}}x{\text{|}} < 1$. Таким образом, на основании (7.8) для функции $y$ получим уравнение

$y = \frac{{x_{{}}^{2} + x - 2}}{{x(x - 3)}}$

Нули данной функции определены равенствами $x = 1$, $x = - 2$; производная $y{\kern 1pt} ' < 0$ при $x \ne 0$, $x \ne 3$. Следовательно, функция $y > {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}$, если $x \in (0,x_{*}^{{}})$, где

$x_{*}^{{}} = \frac{1}{2}( - 5 + \sqrt {41} )$

Из полученных результатов следует, что $a_{{\text{0}}}^{{}} \in (0{\text{,}}x_{*}^{{}})$.

Значение параметра E, который находится из метода рассмотрения уравнений (7.6), (7.7), зависит от параметров задачи, так как E – постоянная первого интеграла. Следовательно, дано обобщение решения п. 5 на случай динамически симметричного гироскопа, которое также описывается эллиптическими функциями.

Заключение. В статье дана постановка задачи о прецессионных движениях тела под действием трех однородных силовых полей. На основе подхода, указанного в [6, 7], получены три дифференциальных уравнения на функции $\varphi (t)$, $\psi (t)$, где $\varphi (t)$ – угол собственного вращения, $\psi (t)$ – угол прецессии тела. Изучена проблема замены одного из уравнений интегралом энергии. Рассмотрены случаи регулярных прецессий ($\dot {\varphi } = n$, $\dot {\psi } = m$, n и m – постоянные, удовлетворяющие условию $m = 2n$) и прецессионно-изоконических движений ($\dot {\varphi } = \dot {\psi }$). Показано, что в первом варианте угол нутации равен ${\text{arccos}}$(1/6), а во втором он имеет значения 60° для сферических гироскопов и $a_{{\text{0}}}^{{}}\, \in \,(0,x_{*}^{{}})$ для динамически симметричного гироскопа относительно вектора a. Рассмотрен вариант, когда движение тела происходит под действием сил двух однородных силовых полей. Получено сведение решения к эллиптическим функциям времени. В отличие от ранее использованных подходов, в данной статье применен векторный метод в описании прецессионных движений тела, что позволило установить решение уравнений движения в кратком виде. Редуцированные уравнения движения дают возможность изучить прецессии общего вида. Например, можно анонсировать результат о том, что исследование прецессионных движений тела в случае, когда скорость собственного вращения тела в два раза меньше его скорости прецессии, угол нутации равен ${\text{arccos(1/4)}}$. В связи с этим представляет интерес сравнение указанных в данной статье результатов с аналогичными результатами в классической задаче, которая описывается уравнениями Эйлера–Пуассона, и в задаче о движении гиростата под действием потенциальных и гироскопических сил, которая описывается уравнениями класса Кирхгофа–Пуассона. В первой задаче известно решение Д. Гриоли [2], характеризующееся прецессией тела относительно наклонной оси. Для него выполняются свойства $\dot {\psi } = \dot {\varphi }$, $\theta _{{\text{0}}}^{{}} = 90^\circ $. Во второй задаче имеет место решение А.В. Мазнева [25], отвечающее прецессии гиростата относительно вертикали и обладающее свойствами $\dot {\psi } = \dot {\varphi }$, $\theta _{{\text{0}}}^{{}} = 120^\circ $.

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

  1. Ишлинский А.Ю. Ориентация, гироскопы и инерциальная навигация. М.: Наука, 1976. 672 с.

  2. Grioli G. Esistenza e determinazione delle precessioni regolari dinamicamente possibili per un solido pesante asimmetrico // Ann. Mat. Pura Appl. 1947. S. 4. V. 26. F. 3–4. P. 271–281.

  3. Klein F., Sommerfeld A. Über die Theorie des Kreisels. N.Y.: Johnson reprint corp., 1965. 966 p.

  4. Bressan A. Sulle precessioni d’un corpo rigido costituenti moti di Hess // Rend. Semin. Mat. Univ. Padova. 1957. V. 27. P. 276–283.

  5. Горр Г.В. Прецессионные движения в динамике твердого тела и динамике систем связанных твердых тел // Прикл. мат. мех. 2003. Т. 67. Вып. 4. С. 573–587.

  6. Горр Г.В., Мазнев А.В., Щетинина Е.К. Прецессионные движения в динамике твердого тела и в динамике систем связанных твердых тел. Донецк: ДонНУ, 2009. 222 с.

  7. Горр Г.В., Мазнев А.В. Динамика гиростата, имеющего неподвижную точку. Донецк: ДонНУ, 2010. 364 с.

  8. Горр Г.В., Мазнев А.В., Котов Г.А. Движение гиростата с переменным гиростатическим моментом. Донецк: ГУ “ИПММ”, 2018. 250 с

  9. Докшевич А.И. Решения в конечном виде уравнений Эйлера−Пуассона. Киев: Наук. думка, 1992. 168 с.

  10. Ольшанский В.Ю. О регулярных прецессиях несимметричного твердого тела с жидким наполнением // Прикл. мат. мех. 2018. Т. 82. Вып. 5. С. 559–571.

  11. Ol’shanskii V.Yu. New cases of regular precession of an asymmetric liquid-filled rigid body // Celest. Mech. Dyn. Astron. 2019. V. 131. № 12. P. 57. doi.org/https://doi.org/10.1007/s10569-019-9929-x

  12. Ol’shanskii V.Yu. Analysis of regular precession conditions for asymmetrical liquid-filled rigid body // Celest. Mech. Dyn Astron. 2020. V. 132. № 9. P. 46. doi.org/https://doi.org/10.1007/s10569-020-09984-2

  13. Ольшанский В.Ю. Полурегулярная прецессия несимметричного твердого тела с жидким наполнением // Прикл. мат. мех. 2021. Т. 85. Вып. 5. С. 547–564. doi.org/https://doi.org/10.31857/S0032823521040111

  14. Рубановский В.Н. Об одном новом частном решении уравнений движения тяжелого твердого тела в жидкости // Прикл. мат. мех. 1985. Т. 49. Вып. 2. С. 212–219.

  15. Горр Г.В., Рубановский В.Н. Об одном новом классе движений системы тяжелых шарнирно связанных твердых тел // Прикл. мат. мех. 1988. Т. 50. Вып. 5. С. 707–712.

  16. Харламов М.П. Критическое множество и бифуркационная диаграмма задачи о движении волчка Ковалевской в двойном поле // Механика твердого тела. 2004. Вып. 34. С. 47–58.

  17. Харламов М.П. Особые периодические движения гиростата Ковалевской в двойном поле // Механика твердого тела. 2007. Вып. 37. С. 85–96.

  18. Рябов П.Е. Алгебраические кривые и бифуркационные диаграммы двух интегрируемых задач // Механика твердого тела. 2007. Вып. 37. С. 97–111.

  19. Харламов М.П., Яхья Х.М. Разделение переменных в одном частном случае движения гиростата в двойном поле // Механика твердого тела. 2014. Вып. 44. С. 7–15.

  20. Yehia H.M. Equivalent problems in rigid body dynamics –II // Celest. Mech. 1988. V. 41. P. 289–295. doi.org/ 01238765https://doi.org/10.1007/BF

  21. Yehia H.M. On the regular precession of an asymmetric rigid body acted upon by uniform gravity and magnetic fields // Egypt. J. Bas. Appl. Sci. 2015. V. 2. № 3. P. 200–205. doi.org/https://doi.org/10.1016/j.ejbas.2015.03.002

  22. Yehia H.M. Regular precession of a rigid body (gyrostat) acted upon by an irreducible combination of three classical fields // J. Egypt. Math. Soc. 2017. V. 25. № 2. P. 216–219. doi.org/https://doi.org/10.1016/j.joems.2016.08.001

  23. Hessein A.M. Precessional motion of a rigid body acted upon by three irreducible fields // Russ.J. Nonlin. Dyn. 2019. V. 15. № 3. P. 285–292. doi.org/https://doi.org/10.20537/nd190307

  24. Горр Г.В., Балаклицкая Т.В. О движении главных осей твердого тела, имеющего неподвижную точку, в случае прецессий относительно вертикали // Механика твердого тела. 2019. Вып. 49. С. 57−66.

  25. Мазнев А.В. Прецессионно-изоконические движения в одном решении уравнений Кирхгофа // Вісн. Донецьк. ун-ту. Сер. А: Природничі науки. 2001. Вып. 2. С. 12–16.

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