Non-Regular Precession of a Gyrostat in Three Uniform Fields

封面

如何引用文章

全文:

详细

This article presents a solution to the problem of the conditions of gyrostat non-regular precession in three homogeneous fields, in which the ratio of precession and proper rotation velocities is constant. The case of a gyrostat with axial dynamic symmetry, the axis of its proper rotation coinciding with the axis of symmetry of the gyrostat, is highlighted. It is shown that the precession is possible only at a precession rate twice the rate of its proper rotation, and the gyrostatic moment deflected from the axis of symmetry by some angle ε. An expression for each of the rates is obtained through elementary functions of time. At 0 < ε < ε*, the motion is periodic, at ε ≥ ε*, the velocity tends to zero and the solid makes no more than one revolution around the axis of its proper rotation, the angle ε* is expressed through the constant nutation angle θ. A relationship has been found between the nutation angle and the ratio of the axial and equatorial moments of inertia, under spherical symmetry cosθ = 1/4. The set of permissible positions of the centers of force at arbitrary given angles between the lines of force of homogeneous fields and for the special case of orthogonal fields is indicated.

全文:

  1. Введение. Классическая задача движения тяжелого твердого тела с неподвижной точкой имеет много обобщений для различных силовых полей. Наиболее исследован случай, когда поле одно или действуют несколько полей с общей осью симметрии. Получены [1] решения для тяжелого тела в магнитном поле, для гиростата под действием потенциальных и гироскопических сил [2, 3].

Важным случаем движения является прецессия и случай симметричного тяжелого тела хорошо известен. Гриоли была найдены [4] условия регулярной прецессии несимметричного тела вокруг оси, отклоненной от вертикали. В наших работах [5–8] показано, что прецессия тела с полостью, заполненной жидкостью, также возможна при отсутствии динамической симметрии. Обзор прецессий твердого тела и гиростата под действием сил различной природы приведен в [9, 10].

Случай, когда направления полей заданы двумя или тремя векторами в инерциальном пространстве, изучен в значительно меньшей степени и исследования в этой области активно проводятся в настоящее время. Первые примеры регулярной прецессии несимметричного твердого тела и гиростата в двух [11] и трех [12] однородных полях были построены Х. Яхья. В этих решениях оси прецессии и собственного вращения ортогональны, а скорости прецессии и собственного вращения совпадают; эти решения можно считать аналогами прецессии Гриоли для двух и трех полей. В наших работах описаны возможные случаи прецессии твердого тела и гиростата в двух [13] и трех [14] однородных полях; найден новый случай [14], когда скорость прецессии вдвое больше скорости собственного вращения, угол между осями прецессии и собственного вращения задан равенством cosθ = 1/6.

Была рассмотрена [15] регулярная прецессия гиростата в трех полях, одно из которых – осесимметричное, и для частного случая, когда скорости прецессии и собственного вращения равны, поля ортогональны и ось прецессии совпадает с осью симметрии неоднородного поля, получены условия, связывающие параметры системы. В нашей работе [16] выполнено исследование всех возможных случаев регулярной прецессии в данной суперпозиции трех полей, найдены конфигурационные условия и центры приведения сил. Показано [16], что прецессия возможна при скорости прецессии равной, вдвое большей или вдвое меньшей скорости собственного вращения. Для известного случая [15] с равными скоростями прецессии и собственного вращения указаны новые решения с осью прецессии, отклоненной от оси симметрии неоднородного поля. Найдены [16] новые случаи регулярной прецессии, когда отношение скоростей прецессии и собственного вращения равно двум либо одной второй. Показано, что в частном случае гиростата, гиростатический момент которого направлен по оси собственного вращения и в случае твердого тела скорость прецессии может быть вдвое меньше скорости собственного вращения, только если угол нутации задан равенством sinθ = 4/5.

Г.В. Горром была рассмотрена задача о нерегулярной прецессии вокруг вертикали динамически симметричного тела в трех однородных ортогональных полях, при которой отношение скоростей прецессии и собственного вращения постоянно [17–19]. В нашей работе [20] проанализированы возможные случаи нерегулярной прецессии динамически симметричного тела в трех однородных полях с постоянным отношением скоростей при произвольных углах между силовыми линиями полей и с произвольным направлением оси прецессии в инерциальном пространстве. В частном случае сферической симметрии тела при скорости прецессии вдвое меньшей или вдвое большей скорости собственного вращения угол нутации определен равенством cosθ = 1/4, при равных скоростях cosθ = 1/2.

В настоящей статье приведено решение задачи об условиях нерегулярной прецессии гиростата в трех однородных полях, при которой отношение скоростей прецессии и собственного вращения постоянно. Условия получены для случая гиростата с осевой динамической симметрией, ось собственного вращения которого совпадает с осью симметрии гиростата. Показано, что нерегулярная прецессия возможна только при скорости прецессии вдвое большей скорости собственного вращения и гиростатическом моменте, отклоненном от оси симметрии. При угле отклонения ε, меньшем ε*, движение периодическое, при ε ≥ ε* скорость стремится к нулю и тело совершает не больше одного оборота вокруг оси собственного вращения; угол ε* выражен через постоянный угол нутации θ. Решение найдено в элементарных функциях. Найдена связь между отношением осевого и экваториального моментов инерции гиростата и углом нутации, при сферической симметрии cosθ = 1/4. Указано множество допустимых положений центров приведения сил при произвольных заданных углах между силовыми линиями однородных полей и для частного случая ортогональных полей.

  1. Постановка задачи. Для описания движения гиростата вокруг неподвижной точки под действием трех полей используем уравнения [12]

Iω˙+ω×Iω+σ=M=α1×u1+α2×u2+α3×u3 (2.1)

α˙i+ω×αi=0; i=1, 2, 3 (2.2)

Здесь ( ) – производная по времени в системе отсчета, связанной с телом; векторы ui постоянны в этой системе, ui =  piOCi, Ci – центры приведения сил, I – оператор инерции тела в неподвижной точке, ω – угловая скорость тела, единичные векторы αi задают направления сил каждого из полей, σ – гиростатический момент, M – момент действующих на тело сил относительно неподвижной точки O.

Прецессия тела задается равенством

ω=ωrm+ωpρ (2.3)

Единичные векторы m и ρ постоянны, соответственно, в подвижной и инерциальной системах. Скалярные функции ωr (t) и ωp (t) – это величины скоростей собственного вращения и прецессии. Прецессия называется регулярной, если обе скорости ωr и ωp постоянны, и нерегулярной, если хотя бы одна из скоростей непостоянна [10].

Рассмотрим прецессии гиростата, для которых, как и в работах [17–20] для твердого тела, отношение скоростей постоянно

ωp/ωr=κ=const (2.4)

Ниже решается следующая задача: при заданных в инерциальной системе отсчета направлениях полей αi и оси прецессии ρ определить при каких ограничениях на оператор I, гиростатический момент σ, векторы ui и отношение скоростей κ гиростат может совершать прецессию и найти зависимость скоростей прецессии и собственного вращения от времени.

Векторная функция ρ (t) удовлетворяет уравнению (2.2), которое, при учете равенства (2.3), становится линейным

ρ˙+ωrtm×ρ=0 (2.5)

Пусть (l1, l2, l3) – некоторый связанный с телом ортонормированный правый базис такой, что l3 = m. Решение уравнения (2.5) имеет вид

ρ=sinθsinτl1+cosτl2+cosθl3; dτ=Ωdt (2.6)

Здесь Ω = ωr (t) – скорость собственного вращения, произвольный параметр θ – это постоянный угол между осями собственного вращения и прецессии (угол нутации), cosθ = (m, ρ). При Ω = const в рассматриваемом случае, когда отношение скоростей  постоянно, прецессия является регулярной.

Векторные функции ω (t), αi (t), как и ранее [14, 16, 20], задаются в связанном с телом ортонормированном базисе (l1, l2, l3) равенствами:

ω=Ωω~, ω~=κsinθsinτl1+cosτl2+1+κcosθl3 (2.7)

αi=Rsi; i=1, 2, 3, ρ=Rl3 (2.8)

Элементы матрицы оператора поворота R в базисе (li) следующие:

r11=cos2θ2cosκ+1τsin2θ2cosκ1τ, r31=sinθsinκτr12=cos2θ2sinκ+1τ+sin2θ2sinκ1τ, r13=sinθsinτr21=cos2θ2sinκ+1τsin2θ2sinκ1τ, r32=sinθcosκτr22=cos2θ2cosκ+1τsin2θ2cosκ1τ, r23=sinθcosτ, r33=cosθ   (2.9)

Функции αi (t), заданные равенствами (2.8), являются решениями линейных (при заданной формулой (2.7) функции ω (t)) уравнений (2.2) при произвольных постоянных (в связанной с телом системе отсчета) векторах si. При известных функциях αi (t) момент внешних сил M задан в подвижной системе отсчета.

Задача, решаемая в работе, состоит в нахождении условий обращения в тождество равенства (2.1) при функциях ω, αi (t), M, заданных равенствами (2.7)–(2.9). Из формул (2.7) получим

ω˙=Ω˙ω~+Ω2κsinθcosτl1sinτl2 (2.10)

Уравнение (2.1) записывается в виде

ΩdΩdτf+Ω2g+Ωh=M (2.11)

f=Iω~, g=ω~×Iω~+κsinθIcosτl1sinτl2, h=ω~×σ (2.12)

Зависимость τ = τ (t) найдем из равенства

t=dτΩτ (2.13)

В общем случае, когда поля не ортогональны, удобно, как ранее в наших работах [14, 16, 20] преобразовать формулу для момента сил следующим образом. Зададим векторы ni и оператор G равенствами

n1=s2×s3s1,s2,s31  2  3, ui=Gni; i=1, 2, 3 (2.14)

Здесь ⟨a, b, c⟩ = (a, b × c), (1 2 3) – знак циклической перестановки.

Всюду в работе рассматриваем случай неприводимых полей и считаем векторы  некомпланарными, тогда, в силу равенств (2.8), ⟨s1, s2, s3⟩ ≠ 0.

Имеем формулу [16] для суммы моментов внешних сил

M=Gl1×Rl1+Gl2×Rl2+Gl3×Rl3 (2.15)

Ниже, в разд. 3, приведена система трех тождеств, выполнение которых необходимо и достаточно для существования искомого решения, описывающего прецессию гиростата с осевой динамической симметрией, для которой отношение скоростей прецессии и собственного вращения постоянно. В разд. 4 доказано, что для совместности названной системы тождеств в случае трех полей, неприводимых к двум или одному полю, необходимо, чтобы отношение скоростей было равно одному, двум или одной второй. Дальнейший анализ показывает, что случаи κ = 1 и κ = 1/2 возможны только при Ω = const, то есть при регулярной прецессии. Кроме того, показано, что вектор гиростатического момента должен быть отклонен от оси динамической симметрии гиростата. Разд. 5 и 6 содержат основные результаты статьи. В разд. 5 показано, что при κ = 2 система трех тождеств совместна и система уравнений (2.1), (2.2) имеет в случае трех неприводимых полей решение, описывающее нерегулярную прецессию с постоянным отношением скоростей. Получено необходимое условие, выражающее угол нутации через отношение осевого и экваториального моментов инерции гиростата. Дано описание допустимых положений центров приведения сил. Найдено выражение Ω (τ) через элементарные функции. В разд. 6 выполнен анализ возможных движений гиростата. Записана в элементарных функциях зависимость от времени компонент угловой скорости гиростата. Показано, что при угле ε отклонения гиростатического момента от оси симметрии гиростата, меньшем ε*, движение периодическое, при ε ≥ ε* скорость стремится к нулю и тело совершает не больше одного оборота вокруг оси собственного вращения; угол ε* выражен через постоянный угол нутации θ.

  1. Прецессия гиростата с осевой симметрией. Определяющие тождества. Рассмотрим прецессию гиростата, имеющего осевую динамическую симметрию, ось симметрии которого совпадает с осью собственного вращения. В этом случае Iij = 0, ij, I1 = I2. Оси l1, l2 выберем так, что σ2 = 0. Запишем заданные формулой (2.12) параметры f, g, h

f=κsinθI1sinτl1+cosτl2+1+κcosθ I3l3 (3.1)

g=λ1cosτl1sinτl2, λ1=defκsinθ1+κcosθ I3κcosθ I1 (3.2)

h=κsinθσ3cosτl1sinτl2+1+κcosθ σ1l2κsinθ σ1cosτl3 (3.3)

Отметим следующее. Если гиростатический момент коллинеарен оси симметрии, то σ1 = 0 и векторы g и h коллинеарны, что упрощает анализ уравнения (2.11). Проведенная проверка показала, что в этом случае искомое решение, описывающее прецессию с постоянным отношением скоростей осесимметричного гиростата в трех неприводимых однородных полях, существует, только если Ω = const. Прецессия является регулярной, все возможные случаи такой прецессии описаны в нашей работе [16]. Всюду ниже при построении условий нерегулярной прецессии с постоянным отношением скоростей (то есть при условии (2.4)) считаем σ1 ≠ 0, гиростатический момент  при этом отклонен от оси симметрии.

Запишем уравнение (2.11) в проекциях на оси

κsinθI1sinτΩdΩdτ+λ1cosτΩ2+κsinθσ3cosτΩ=M1 (3.4)

κsinθI1cosτΩdΩdτλ1sinτΩ2+1+κcosθσ1κsinθσ3sinτΩ=M2 (3.5)

1+κcosθI3ΩdΩdτκsinθσ1cosτΩ=M3 (3.6)

Уравнения (3.4), (3.5) эквивалентны системе

κsinθI1ΩdΩdτ+1+κcosθσ1cosτΩ=M1sinτ+M2cosτ=defM~1 (3.7)

λ2=κsinθ2I1+1+κcosθ2I3 (3.8)

Из уравнений (3.6) и (3.7) находим

λ2σ1cosτΩ=1+κcosθI3M~1κsinθI1M3 (3.9)

λ2ΩdΩdτ=κsinθM~1+1+κcosθM3 (3.10)

Здесь

λ2=κsinθ2I1+1+κcosθ2I3 (3.11)

Таким образом, исходное векторное уравнение (2.11) эквивалентно системе трех уравнений (3.8)–(3.10).

Из формул (2.9) и (2.15) получим компоненты Mi момента в базисе (li) и затем запишем величины M~i, заданные в формулах (3.7) и (3.8)

M~1=cosθG13cosτG23sinτ+G32sinκτG31cosκτ++sinθ2G12+G21cosκ+1τ+G11G22sinκ+1τ++G12G21cosκ1τ+G11+G22sinκ1τ (3.12)

M~2=cosθG23cosτ+G13sinτ+G32cosκτ+G31sinκτsinθ2(2G33+G11G22cosκ+1τG12+G21sinκ+1τG11+G22cosκ1τ+G12G21sinκ1τ) (3.13)

M3=sinθG13cosτ+G23sinτ++sin2θ2G21G12cosκ1τG11+G22sinκ1τ++cos2θ2G12+G21cosκ+1τ+G11G22sinκ+1τ (3.14)

Правые части уравнений (3.9), (3.10) записываются в виде

κsinθM~1+1+κcosθM3==sinθG13cosτ+G23sinτ+κG32sinκτG31cosκτ++1+cosθ2κ+1G12+G21cosκ+1τ+G11G22sinκ+1τ++1cosθ2κ1G12G21cosκ1τ+G11+G22sinκ1τ (3.15)

1+κcosθI3M~1κsinθI1M3=1+κcosθI3G32sinκτG31cosκτ++λ3G13cosτG23sinτ++12λ4G12+G21cosκ+1τ+G11G22sinκ+1τ++12λ5G12G21cosκ1τ+G11+G22sinκ1τ (3.16)

λ3=cosθ1+κcosθI3+κsinθ2I1λ4=sinθ1+κcosθ I3κ1+cosθ I1λ5=sinθ1+κcosθI3+κ1cosθI1 (3.17)

  1. Предварительный анализ. Необходимо определить условия, при которых равенства (3.8)–(3.10) тождественно выполнены для некоторой функции Ω (τ) ≠ const и найти эту функцию. Из формулы (3.10) получим

λ22Ω2=sinθG13sinτG23cosτG31sinκτG32cosκτ++1+cosθ2(G12+G21sinκ+1τG11G22cosκ+1τ)++1cosθ2(G12G21sinκ1τG11+G22cosκ1τ)+C (4.1)

Если κ ≠ 1, то необходимы условия

G11=G22, G12=G21 (4.2)

Формулы (3.9) и (4.1) при условиях (4.2) записываются в виде

λ2σ1cosτΩ=1+κcosθI3G32sinκτG31cosκτ++λ3G13cosτG23sinτ+λ5G12cosκ1τ+G11sinκ1τ (4.3)

λ22Ω2=sinθG13sinτG23cosτG31sinκτG32cosκτ++1cosθG12sinκ1τG11cosκ1τ+C (4.4)

Если функцию Ω (τ) из равенства (4.3) подставить в равенство (4.4), то получим тождество, одним из условий выполнения которого при κ ≠ 1, κ ≠ 2, κ ≠ 1/2 является G11 = G12 = 0. Учитывая условие (4.2), получим, что матрица  – вырожденная и имеем приводимый случай. Таким образом, искомое решение, описывающее нерегулярную прецессию гиростата в трех полях с постоянным отношением скорости прецессии к скорости собственного вращения, может существовать только в одном из указанных случаев. В нашей работе [20] показано, что нерегулярная прецессия твердого тела в трех однородных полях возможна, когда отношение скорости прецессии к скорости собственного вращения равно одному, двум или одной второй. Ниже показано, что аналогичное решение для гиростата существует, только если скорость прецессии вдвое больше скорости собственного вращения.

Рассмотрим случай κ = 1. Условия (3.9), (3.10) записываются в виде

λ2σ1cosτΩ=1+cosθI3G32sinτG31cosτ+λ3G13cosτG23sinτ++12λ4G12+G21cos2τ+G11G22sin2τ+12λ5G12G21 (4.5)

λ2Ω˙=λ2ΩdΩdτ=sinθG13+G31cosτ+G23+G32sinτ++1+cosθG12+G21cos2τ+G11G22sin2τ (4.6)

В соответствии с формулой (3.11) λ2 ≠ 0 и, так как σ1 ≠ 0, то правая часть в формуле (4.5) должна делиться на cosτ, для этого необходимо выполнение условий

1+cosθI3G32λ3G23=0, λ4G12+G21=λ5G12G21 (4.7)

При выполнении этих условий формула (4.5) упрощается:

Ω=λ4(G12+G21cosτ+G11G22sinτ)+λ3G131+cosθI3G31λ2σ1 (4.8)

Если из формулы (4.8) подставить Ω в равенство (4.6) и сравнить коэффициенты при cos2τ, sin2τ, то получим необходимые условия

2xy=y, x2y2=x, x=G11G22ν, y=G12+G21ν, ν=21+cosθλ2σ12λ42 (4.9)

Если y ≠ 0, то x = 1/2 и второе уравнение (4.9) не имеет действительных решений. Если y = 0, то x = 0 или x = 1. В первом случае из формулы (4.8) следует Ω = const. Остается только случай y = 0, x = 1 и получаем условия

G12=G21, G11G22=ν (4.10)

Формула (4.8) и равенство (4.6) принимают вид

Ω=1λ2σ1λ4νsinτ+λ3G131+cosθ  I3G31 (4.11)

λ2ΩdΩdτ=sinθG23+G32sinτG13+G31cosτ+1+cosθνsin2τ (4.12)

Равенства (4.11) и (4.12) совместны при условиях

G23+G32=0, λ4νλ2σ12λ3G131+cosθI3G31=sinθG13+G31 (4.13)

Так как λ3+1+cosθI3=(1+cosθ)2I3+(sinθ)2I10, то из первых условий (4.7) и (4.13) следует G23 = G32 = 0.

Осталось рассмотреть возможность выполнения тождества (3.8), которое при полученных условиях записывается в виде

1+cosθσ1sinτsinθ  σ3Ωλ1Ω2==sinθ  2νcos2τ+2G33G11G22+cosθG13+G31sinτ

Сравнение здесь слагаемых с cos2τ приводит либо к условию ν = 0, но тогда Ω = const, либо к условию 1+cosθ2λ1λ4+sinθλ2=0. Данное условие приводится к виду (1+cosθ)2I3+(sinθ)2I1=0 и не может быть выполнено. Таким образом, при κ = 1 искомого решения нет.

Покажем, что этот случай при κ = 1/2 также невозможен. Из формулы (3.9) при учете условия (4.2) получим

λ2σ1cosτΩ=acosτ+bsinτ+ccosτ2+dsinτ2 (4.14)

Здесь a=λ3G13, b=λ3G23, c=λ5G121+cosθ2I3G31, d=1+cosθ2I3G32λ5G11.

Для делимости правой части в формуле (4.14) на cosτ необходимо b = c = d = 0. Формула принимает вид λ2σ1Ω = const. Заданный формулой (3.11) параметр λ2 больше нуля, следовательно σ1Ω = const. При σ1 ≠ 0 отсюда следует Ω = const и прецессия будет регулярной. В разд. 3 показано, что рассмотрение случая σ1 = 0 также приводит к условию Ω = const. Таким образом, при κ = 1/2 искомая нерегулярная прецессия невозможна.

  1. Построение условий прецессии в случае κ = 2. Формулы (4.3), (4.4) при условиях (4.2) дают

λ2σ1cosτΩ=1+2cosθI3G32sin2τG31cos2τ++λ3G13+λ5G12cosτ+λ5G11λ3G23sinτ (5.1)

λ22Ω2=sinθG31sin2τ+G32cos2τsinθG23+1cosθG11cosτ++sinθG13+1cosθG12sinτ+C (5.2)

Для делимости на cosτ в формуле (5.1) необходимы условия

G31=0, λ5G11λ3G23=0 (5.3)

При этом

λ2σ1Ω=21+2cosθI3G32sinτ+λ5G12+λ3G13 (5.4)

Отсюда следует G32 ≠ 0, иначе Ω = const, и получаем

Ω2=1λ2σ12(1+2cosθI3G32221cos2τ++41+2cosθI3G32λ5G12+λ3G13sinτ+λ5G12+λ3G132) (5.5)

Сравнение с формулой (5.2) приводит к условию

sinθG23+1cosθG11=0 (5.6)

Формула (5.2) при условиях (5.3), (5.6) дает

Ω2=2λ2sinθG32cos2τ+sinθG13+1cosθG12sinτ+C (5.7)

Сравнивая формулы (5.5) и (5.7), приходим к условиям

1λ2σ12(1+2cosθ)I3G322=1λ2sinθG322λ2σ121+2cosθI3G32λ1G12+λ2G13=1λ2sinθG13+1cosθG12

Из первого условия находим

G32=sinθλ2σ12(1+2cosθI3)2 (5.8)

Из второго условия при учете формулы (5.8) следует связь

sinθG13+1cosθG12=2sinθ1+2cosθI3λ5G12+λ3G13 (5.9)

По формуле (3.13) получим

M~2=cosθG32cos2τ+sinθG11+cosθG23cosτ++cosθG13sinθG12sinτsinθG33 (5.10)

В левой части равенства (3.8) при учете формулы (5.4) нет слагаемых с cosτ, поэтому в M~2 такие члены должны отсутствовать и из формулы (5.10) получаем условие

sinθG11+cosθG23=0 (5.11)

Определитель системы (5.6), (5.11) не равен нулю, следовательно

G11=G22=G23=0 (5.12)

Формула (5.10) принимает вид

M~2=cosθG32cos2τ+cosθG13sinθG12sinτsinθG33 (5.13)

Добьемся теперь выполнения равенства (3.8), используя формулы (5.4) и (5.13). Сравнивая коэффициенты при cos2τ, получим связь

2λ1λ2σ12(1+2cosθ I3)2G3221λ2σ11+2cosθ2I3G32σ1=cosθG32

Подставив сюда G32 из формулы (5.8), приходим к условию

2λ1sinθ1+2cosθ2I3=μcosθ

Данное условие при использовании формул (3.2), (3.11) для λ1, λ2 приводится к виду

4cosθ1cosθI1=14cos2θI3 (5.14)

Отметим, что в случае сферической симметрии I1 = I2 = I3 из условия (5.14) получаем

cosθ=14 (5.15)

Такое же значение отмечено в работах [17] и [20] для частного случая твердого тела со сферической симметрией.

При имеющемся ограничении I3 < 2I1 допустимые значения угла нутации определяются условием

0<cosθ<1231 (5.16)

Сравнивая в тождестве (3.8) коэффициенты при sinτ и свободные члены, получим связи

14λ1λ2σ12I3G32λ5G12+λ3G13+λ2sinθG12cosθG131+2cosθ  =4σ3σ1sinθI3G32 (5.17)

G33=2σ3λ5G12+λ3G13λ2σ1+λ121+2cosθI3)2G322+λ5G12+λ3G132sinθλ2σ121+2cosθ2I3G32λ2sinθ (5.18)

Равенства (5.9) и (5.17) позволяют найти величины G12 и G13, из равенства (5.18) затем найдем .

Заданные формулами (3.17) параметры λ3 и λ5 при связи (5.14) записываются в виде

λ5=sinθ1+2cosθ2cosθI3, λ3=1cosθ1+2cosθ2cosθI3 (5.19)

Условие (5.9) при учете формул (5.19) преобразуется к виду

1cosθG12+sinθG13=0 (5.20)

Полученная связь позволяет записать G12, G13 в виде

G12=sinθg, G13=1cosθg (5.21)

С учетом формул (5.19), (5.20) получим

λ5G12+λ3G13=sinθG12cosθG13I3=1cosθ1+2cosθI3g (5.22)

Параметр g определим из равенства (5.17)

g=4σ3cosθσ1sinθ12cosθG32

Подставляя сюда G32 из формулы (5.8), получим

g=4σ1σ314cos2θI3 (5.23)

Формулы (5.21) принимают вид

G12=4σ1σ3sinθ14cos2θI3, G13=4σ1σ31cosθ14cos2θI3 (5.24)

Из формулы (5.18) при учете формул (5.21) и (5.23) находим

G33=σ12I31+2cosθ+2σ32I112cosθ (5.25)

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

I3I1=4cosθ1cosθ14cos2θ (5.26)

и элементы матрицы G, определяющей положения центров приведения, задавались формулами

G12=G21=4σ1σ3sinθI314cos2θ, G13=4σ1σ31cosθI314cos2θ , G32=σ12I3sinθcosθ  1+2cosθG33=σ12I31+2cosθ+2σ32I112cosθ, G11=G22=G23=G31=0 (5.27)

Здесь σ1 и σ2 – проекции гиростатического момента на экваториальную плоскость и на ось симметрии эллипсоида инерции.

Из формулы (5.4) получаем выражение для скорости собственного вращения

Ω=2σ1I3sinθ1+2cosθ  sinτ+σ3I1 (5.28)

Допустимые положения центров приведения сил определяются при заданной матрице G из формул (2.14).

В случае ортогональных полей единичные векторы α1, α2, α3 образуют правую ортогональную тройку, тогда, в соответствии с формулами (2.8), единичные векторы s1, s2, s3 образуют левую ортогональную тройку. Из формул (2.14) получим ni = si, i = 1, 2, 3. Если ось прецессии совпадает с осью симметрии неоднородного поля, то ρ = α3 и, в соответствии с формулой (2.8), s3 = l3.

  1. Анализ движения гиростата. Получим явную зависимость переменной τ от времени и зависимость от времени компонент угловой скорости гиростата. Запишем формулу (5.28) в виде

Ω=a+bsinτ, a=σ3I1, b=2σ1I3sinθ1+2cosθ (6.1)

Интеграл (2.13) выражается через элементарные функции. При a2 > b2 движение – периодическое, при b2a2 скорость вращения стремится к нулю.

Из формул (5.26), (6.1) получим

a2b2=ctg2ε41cosθcos2θ1+cosθ12cosθ2,  ab=ctgεtgθ22cosθ12cosθ (6.2)

Здесь ε – угол между вектором гиростатического момента σ и осью симметрии l3, σ1 = σsinε, σ3 = σcosε. Движение периодическое при |ε| < ε*, иначе – затухающее

tg2ε*=41cosθcos2θ1+cosθ12cosθ2,  tgε*=2cosθ12cosθtgθ2 (6.3)

В случае сферической динамической симметрии гиростата cosθ = 1/4,

tg2ε*=35 ; ε*3746'

Случай 1: a = ±b

t+C=1atgτ2π4, sinτ=±1a2(t+C)21+a2(t+C)2, cosτ=2at+C1+a2(t+C)2ωr=12ωp=Ω=2a1+a2(t+C)2

Носитель гиростата приближается к состоянию покоя, Ω > 0 при t > ∞.

Случай 2: a2 > b2

tg12τπ2=a2b2abtgt~2; t~=a2b2t+const

Отсюда находим

cosτπ2=acost~babcost~, sinτπ2=a2b2sint~abcost~, (6.4)

ωr=12ωp=Ω=a2b2abcost~ (6.5)

Скорость собственного вращения Ω и вдвое большая скорость прецессии являются периодическими функциями времени с периодом T=2π/a2b2.

Знак Ω совпадает со знаком параметра a и, в силу формулы (6.1), со знаком σ3. Учитывая связь (2.6) переменных τ и t, получим, что при увеличении времени t переменная τ монотонно возрастает, если σ3 > 0 (гиростатический момент образует острый угол с осью собственного вращения) и монотонно убывает, если σ3 < 0.

Наибольшее и наименьшее значения Ω задаются равенствами

maxΩ=a+b, minΩ=ab

Обозначим

ft~,χ=1χ2  sint~  1χcost~  ; χ=ba

При |χ| ≈ 1 вращение тела происходит очень неравномерно, так как, в соответствии с формулой (6.5) почти при всех t~ скорость Ω близка к нулю и только в малой окрестности t~max скорость возрастает до величины, близкой к 2a. На интервале (0, 2π) функция ft~, χ имеет максимум (f = 1) и минимум (f = –1) в точках t~1,2, заданных равенством cost~1,2=χ. При |χ| ≈ 1 величина ft~, χ почти всюду близка к нулю, кроме малых окрестностей точек t~10, t~22π. Обозначив δ = (1 − χ2)1/4, при малых δ получим оценку

ft~,χ<δ1+Oδ2 при t~2δ,2π2δ (6.6)

Пусть a, b > 0. Из формул (6.4), (6.6) следует, что на большей части периода параметр τ находится вблизи значения τ = π / 2. При прохождении почти полного периода, соответствующему изменению t~ на величину 2π – 4δ, происходит поворот вокруг оси собственного вращения на малый угол Δτ ≈ 2δ (и поворот вокруг оси прецессии на угол, вдвое больший). Затем, за малую часть периода Δt~=4δ происходит поворот на угол Δτ ≈ 2π – 2δ. Отметим, что аналогичное исследование неравномерности вращения тела с полостью, наполненной жидкостью, выполнено в нашей работе [8].

Случай 3: b2 > a2. После интегрирования получаем

t~=lnΦ; Φ=atgτ2+bb2a2atgτ2+b+b2a2 (6.7)

Отсюда следует t~ при τ → τ1, 2, где

tgτ12=b2a2ba, tgτ22=b2a2+bacosτ1,2=±b2a2a, sinτ1,2=ab (6.8)

Пусть τ1, τ2 ∈ (−π, π). В зависимости от знаков величин a и b скорость Ω, заданная равенством (6.1), может быть положительна либо отрицательна в любой момент времени. Обозначим

Δτ=τ2τ1, γ=arcsinab

Рассмотрим возможные случаи, используя формулы (6.8)

a>0, b>0: τ1=γ, τ2=γπ, Δτ=2γπ<0, Ω<0, Φ<0a<0, b>0: τ1=γ, τ2=πγ, Δτ=π2γ>0, Ω>0, Φ<0a>0, b<0: τ1=πγ, τ2=γ, Δτ=2γπ<0, Ω<0, Φ<0a<0, b<0: τ1=γπ, τ2=γ, Δτ=π2γ>0, Ω>0, Φ<0

Так как 0 < γ < π / 2, то во всех случаях максимально возможный поворот вокруг собственной оси равен |Δτ| = π − 2γ < π.

Учитывая, что во всех случаях Φ < 0, из формулы (6.7) получим

tgτ2=b2a2bb2a2+bexpt~aexpt~+1cosτ=b2a2b2a2cht~+bsht~b2cht~+bb2a2sht~+a2, sinτ=a  bcht~+b2a2sht~+bb2cht~+bb2a2sht~+a2ωr=12ωp=Ω=aa2b2b2cht~+bb2a2sht~+a2

При t > ±∞ получим

cosτ>b2a2b, sinτ>ab

Носитель гиростата приближается к состоянию покоя, Ω → 0. Угол собственного вращения изменяется на величину, меньшую π, угол прецессии – меньше, чем на 2π. При χ → 1 получаем |τ2 − τ1| = π. Ось собственного вращения не может совершить вокруг оси прецессии больше одного оборота.

Заключение. Известные точные решения задачи о вращении твердого тела в суперпозиции трех однородных полей описывают регулярную прецессию [12, 14], либо нерегулярную прецессию с постоянным отношением скоростей прецессии и собственного вращения [17–20]. Для гиростата в трех однородных полях также получены [12, 18] условия регулярной прецессии. В настоящей работе получены условия нерегулярной прецессии с постоянным отношением скоростей гиростата с осевой симметрией в трех неприводимых однородных полях. Показано, что возможен только случай, когда скорость прецессии вдвое больше скорости собственного вращения и гиростатический момент отклонен от оси собственного вращения. Приведено выражение каждой из скоростей через элементарные функции времени. Выделены случаи периодического и затухающего движений.

×

作者简介

V. Olshansky

Institute of Precision Mechanics and Control of the RAS

编辑信件的主要联系方式.
Email: olshanskiy_vlad@mail.ru
俄罗斯联邦, Saratov

参考

  1. Bogoyavlensky O.I. Euler equations on finite dimensional Lie algebras arising in physical problems // Math. Phys. Commun., 1984, vol. 95, pp. 307–315.
  2. Yehia H.M. On the motion of a rigid body acted upon by potential and gyroscopic forces I. The equations of motion and their transformation // J. Theor. & Appl. Mech., 1986, vol. 5, no. 5, pp. 747–754.
  3. Yehia H.M., El-kenani H.N. Effect of the gravity and magnetic field to find regular precessions of a satellite-gyrostat with principal axes on a circular orbit // J. Appl. Comput. Mech., 2021, vol. 7(4), pp. 2120 – 2128.
  4. Grioli G. Esistenza e determinazione delle prezessioni regolari dinamicamente possibili per un solido pesante asimmetrico // Ann. Mat. Pura e Appl., 1947, vol. 26, iss. 3–4, pp. 271–281.
  5. Ol’shanskii V.Yu. On the regular precession of an asymmetric liquid-filled rigid body // Mech. of Solids, 2018, vol. 53 (suppl. 2), pp. 95–106.
  6. Ol’shanskii V.Yu. New cases of regular precession of an asymmetric liquid-filled rigid body // Celest. Mech. Dyn. Astron., 2019, vol. 131, iss. 12, art. no. 57.
  7. Ol’shanskii V.Yu. Analysis of regular precession conditions for asymmetrical liquid-filled rigid bodies // Celest. Mech. Dyn. Astron., 2020, vol. 132, iss. 9, art. no. 46.
  8. Ol’shanskii V.Yu. Semi-regular precession of an asymmetrical rigid body filled with a liquid // Mech. of Solids, 2021, vol. 56, iss. 8, pp. 1500–1513.
  9. Gorr G.V., Kovalev A.M. Motion of a Gyrostat. Kyev: Nauk. Dumka, 2013. 408 p.
  10. Gorr G.V., Maznev A.V., Shchetinina E.K. Precession Motions in Rigid Body Dynamics and Dynamics of Linked Rigid Bodies Systems. Donetsk: Donetsk National Univ. 2009. 222 p. (in Russian)
  11. 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, vol. 2, iss.3, pp. 200–205.
  12. 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, vol. 25, iss. 2, pp. 216–219.
  13. Ol’shanskii V.Yu. Regular precession of a rigid body in two uniform fields // Mech. Res. Commun., 2023, vol.127, art. no. 104041.
  14. Ol’shanskii V.Yu. Regular precession of a gyrostat in three uniform fields // Mech. of Solids, 2022, vol.57, iss. 8, pp. 1873–1884.
  15. Hussein A.M. Precessional motion of a rigid body acted upon by three irreducible fields // Rus. J. Nonlin. Dyn., 2019, vol. 15, iss. 3, pp. 285–292.
  16. Ol’shanskii V.Yu. Regular precession of a gyrostat in three force fields // Mech. of Solids, 2023, vol. 58, iss. 7, pp. 2515–2530.
  17. Gorr G.V. One class of resonance precession motions of a rigid body under the action of three homogeneous force fields // JAMM, 2023, vol. 87, iss. 1, pp. 3–18.
  18. Gorr G.V. Statement of the problem on precessions of a rigid body with a fixed point in three homogeneous force fields. Precession-isoconic motions of a rigid body // Izv. RAS. Mech. of Solids, 2023, no.3, pp. 123–134. (in Russian).
  19. Gorr G.V. On a class of precessions of a rigid body with a fixed point under the action of forces of three homogeneous force field // Rus. J. Nonlin. Dyn., 2023, vol. 19, iss. 2, pp. 249–264.
  20. Ol’shanskii V.Yu. Nonregular precession of a rigid body in three uniform fields // JAMM, 2024, vol. 88, iss. 1, pp. 17–33.

补充文件

附件文件
动作
1. JATS XML

版权所有 © Russian Academy of Sciences, 2024

Согласие на обработку персональных данных с помощью сервиса «Яндекс.Метрика»

1. Я (далее – «Пользователь» или «Субъект персональных данных»), осуществляя использование сайта https://journals.rcsi.science/ (далее – «Сайт»), подтверждая свою полную дееспособность даю согласие на обработку персональных данных с использованием средств автоматизации Оператору - федеральному государственному бюджетному учреждению «Российский центр научной информации» (РЦНИ), далее – «Оператор», расположенному по адресу: 119991, г. Москва, Ленинский просп., д.32А, со следующими условиями.

2. Категории обрабатываемых данных: файлы «cookies» (куки-файлы). Файлы «cookie» – это небольшой текстовый файл, который веб-сервер может хранить в браузере Пользователя. Данные файлы веб-сервер загружает на устройство Пользователя при посещении им Сайта. При каждом следующем посещении Пользователем Сайта «cookie» файлы отправляются на Сайт Оператора. Данные файлы позволяют Сайту распознавать устройство Пользователя. Содержимое такого файла может как относиться, так и не относиться к персональным данным, в зависимости от того, содержит ли такой файл персональные данные или содержит обезличенные технические данные.

3. Цель обработки персональных данных: анализ пользовательской активности с помощью сервиса «Яндекс.Метрика».

4. Категории субъектов персональных данных: все Пользователи Сайта, которые дали согласие на обработку файлов «cookie».

5. Способы обработки: сбор, запись, систематизация, накопление, хранение, уточнение (обновление, изменение), извлечение, использование, передача (доступ, предоставление), блокирование, удаление, уничтожение персональных данных.

6. Срок обработки и хранения: до получения от Субъекта персональных данных требования о прекращении обработки/отзыва согласия.

7. Способ отзыва: заявление об отзыве в письменном виде путём его направления на адрес электронной почты Оператора: info@rcsi.science или путем письменного обращения по юридическому адресу: 119991, г. Москва, Ленинский просп., д.32А

8. Субъект персональных данных вправе запретить своему оборудованию прием этих данных или ограничить прием этих данных. При отказе от получения таких данных или при ограничении приема данных некоторые функции Сайта могут работать некорректно. Субъект персональных данных обязуется сам настроить свое оборудование таким способом, чтобы оно обеспечивало адекватный его желаниям режим работы и уровень защиты данных файлов «cookie», Оператор не предоставляет технологических и правовых консультаций на темы подобного характера.

9. Порядок уничтожения персональных данных при достижении цели их обработки или при наступлении иных законных оснований определяется Оператором в соответствии с законодательством Российской Федерации.

10. Я согласен/согласна квалифицировать в качестве своей простой электронной подписи под настоящим Согласием и под Политикой обработки персональных данных выполнение мною следующего действия на сайте: https://journals.rcsi.science/ нажатие мною на интерфейсе с текстом: «Сайт использует сервис «Яндекс.Метрика» (который использует файлы «cookie») на элемент с текстом «Принять и продолжить».