$$E(\theta)=E_x+E_y=\frac{J_x\omega_x^2}{2}+\frac{J_y\omega_y^2}{2}=\frac{L_x^2}{2J_x}+\frac{L_y^2}{2J_y}=\frac{L^2}{2 J_y}+\frac{L^2}{2}\left(\frac{1}{J_x}-\frac{1}{J_y}\right) \cos ^2 \theta$$
Кинетическая энергия и момент импульса сохраняются, и $\cos^2\theta$ может быть получен из $\text{A2}$. Следовательно, набор значений, которые может принимать $\theta$, дискретен, и $\theta$ не может изменить свое начальное значение. Поэтому ось симметрии движется вокруг $\vec{L}$, описывая конус с углом раствора $2\theta_0$. Следовательно, $$\psi=2\theta_0$$
Ось симметрии покоится относительно вращающейся системы отсчета, потому что $\theta(t)=\theta_0$, и ось симметрии всегда принадлежит плоскости $Cx_1y_1$. Следовательно, $\vec{\omega}_s$ должен быть коллинеарен оси симметрии в любой момент времени. Таким образом, $$\gamma_s(t)=0.$$ Проецируя сумму $\vec{\Omega}+\vec{\omega}_s$ на оси $Cx$ и $Cy$ для любого момента времени $t$ $$\omega_s+\Omega\cos\theta=\omega_x=\frac{L\cos\theta}{J_x}$$$$\Omega\sin\theta=\omega_y=\frac{L\sin\theta}{J_y}.$$ Откуда $$\Omega(t)=\frac{L}{J_y}.$$ Учитывая, что $\theta(t)=\theta_0$ $$\omega_s(t)=\left(\frac{1}{J_x}-\frac{1}{J_y}\right)L\cos\theta_0.$$ $\omega_s$ также не зависит от времени.
Взаимодействие остатка топлива с баком может быть рассмотрено как внутренняя сила. Следовательно, как и раньше, внешние моменты и силы не действуют на систему, и момент импульса сохраняется. Так как $J_x < J_y$, Из $\text{A2}$ можно показать, что $E(\cos\theta)$ достигает минимума при $\theta=\pi/2$. Таким образом, $$\theta_2=\frac\pi2$$
$\text{B1}$ показывает, что после переходного процесса ступень будет вращаться вокруг оси, перпендикулярной оси симметрии. Конечная угловая скорость $$\omega_2=\frac{L}{J_y}=\frac{\sqrt{J_x^2\cos^2\gamma_1+J_y^2\sin^2\gamma_1}}{J_y}\omega_1 \approx \frac{5}{9} ~рад /с \approx 0.556 ~рад /с$$
Разрежем сферу на кольца так, чтобы $\dot{\vec{B}}$ был перпендикулярен их плоскостям и введем угол $\varphi$, как показано на рисунке.
По закону Фарадея, ЭДС, вызванная в таком кольце изменением поля, равна $$\mathcal{E}=d\Phi/dt=S\dot B=\pi R^2\sin^2\varphi\dot B.$$ Сопротивление кольца равно $$dr=\frac{2\pi\rho R\sin\varphi}{DRd\varphi}.$$ Ток в кольце $$dI=\mathcal{E}/dr=\frac{1}{2\rho}DR^2\dot B\sin\varphi d\varphi.$$ И, наконец, магнитный момент $$d\mu=SdI=\frac{\pi}{2\rho}DR^4\dot B\sin^3\varphi d\varphi=\frac{1}{4\rho}\dot BdJ,$$ где $dJ$ – момент инерции кольца единичной плотности относительно центральной оси, параллельной $y$. таким образом, $$\mu=\frac{1}{4\rho}J\dot B=\frac{2\pi}{3\rho}DR^4\dot B,$$ где $J$ – момент инерции сферы единичной плотности относительно оси, проходящей через его центр. С учетом направления, $$\mu_x=0$$ $$\mu_y=-\frac{2\pi}{3\rho}DR^4\dot B$$ $$\mu_z=0$$
Момент сил задается уравнением $\vec{M}=[\vec{\mu},\vec{B}]$. Он направлен вдоль оси $z$ и равен $$M_z=\mu B\sin\alpha=\frac{2\pi}{3\rho}DR^4B\dot B\sin\alpha$$
Пусть $R_O$ – радиус орбиты. Формула для поля диполя в точке $\vec{r}=(R_O\cos u,R_O\sin u,0)$ и магнитным моментом $\vec{\mu}=(0,-\vec{\mu}_E,0)$ дает $$B_X=-\frac{3}{2}\frac{\mu_0\mu_E}{4\pi R_O^3}\sin2u,$$ $$B_Y=(1-3\sin^2u)\frac{\mu_0\mu_E}{4\pi R_O^3},$$ $$B_Z=0.$$ В точке, где орбита пересекает плоскость экватора ($u=0$), магнитное поле равно $$B_X=0,$$ $$B_Y=\frac{\mu_0\mu_E}{4\pi R_O^3},$$ $$B_Z=0.$$ Отсюда $B_0=\frac{\mu_0\mu_E}{4\pi R_O^3}.$ Наконец, магнитное поле Земли равно $$B_{EX}(u)=-\frac32B_0\sin2u$$ $$B_{EY}(u)=\frac12(3\cos2u-1)B_0$$ $$B_{EZ}(u)=0$$
Использование пункта $\text{C2}$ требует поиска производной магнитного поля в СО ступени. Введем СО $xyz$, у которой $z$ коллинеарно $Z$, и плоскость $xy$ повернута на угол $\beta$ относительно $XY$. Магнитное поле в этой СО равно $$B_x=B_X\cos\beta+B_Y\sin\beta,$$ $$B_y=-B_X\sin\beta+B_Y\cos\beta,$$$$B_z=B_Z=0.$$ Тогда производная магнитного поля $$\dot B_x=\dot B_X\cos\beta+\dot B_Y\sin\beta+(-B_X\sin\beta+B_Y\cos\beta)\dot\beta=$$$$=(B_X'(u)\cos\beta+B_Y'(u)\sin\beta)\dot u+(-B_X\cos\beta-B_Y\sin\beta)\dot\beta,$$$$\dot B_y=-\dot B_X\sin\beta+\dot B_Y\cos\beta+(-B_X\cos\beta-B_Y\sin\beta)\dot\beta=$$$$=(- B_X'(u)\sin\beta+B_Y'(u)\cos\beta)\dot u+(-B_X\cos\beta-B_Y\sin\beta)\dot\beta,$$ $$\dot B_z=0.$$ Подставляя $\dot u=2\pi/T$ и $\dot\beta=\omega$ мы получаем, что момент направлен вдоль $z$ и равен $$M_Z=\frac{2\pi}{3\rho}DB^2_0R^4\left(\frac{3\pi}{T}(3-\cos2u)-\frac\omega2(5-3\cos2u)\right)$$
Усредним $M_z$ по $u$ и используем полученное выражение вместо $M_Z$. Это поможет избавиться от членов, содержащих $\cos2u$ :$$\langle M_z\rangle=\frac{2\pi}{3\rho}DB_0^2R^4\left(\frac{9\pi}{T}-\frac{5\omega}{2}\right).$$ Так как момент направлен вдоль оси вращения, он не меняет ее направления. Это означает, что полученная формула всегда описывает динамику вращательного движения ракеты. Мы предполагаем, что переходный процесс завершен. Тогда из $\text{B1}$ следует, что ступень ракеты вращается вокруг оси, которая перпендикулярна оси симметрии. Таким образом, момент импульса ступени равен $$L_Z=J_y\omega.$$ Так как $\dot L_Z=M_Z$, уравнение на угловую скорость имеет вид $$\dot\omega=\frac{2\pi}{3J_y\rho}DB_0^2R^4\left(\frac{9\pi}{T}-\frac{5\omega}{2}\right).$$ Его решение это $$\omega(t)=\frac{18\pi}{5T}+\left(\omega_2-\frac{18\pi}{5T}\right)e^{-\delta t},$$ где $\delta=\frac{5\pi}{3J_y\rho}DB_0^2R^4$.
Из $\text{C5}$ следует, что угловая скорость стремится к $18\pi/2T$. Таким образом отношение двух периодов $$\frac{T}{T_s(\infty)}=\frac{T\omega(\infty)}{2\pi}=9/5=1.8$$