[论文解读] General solution to the Euler-Poisson equations of a free Lagrange top directly for the rotation matrix

Alexei A. Deriglazov|arXiv (Cornell University)|Mar 4, 2023

Aerospace Engineering and Control Systems被引用 5

一句话总结

本文提出了一种自由拉格朗日陀螺的欧拉-泊松方程的一般解析解，直接以旋转矩阵 $ R_{ij}(t) $ 的形式表示，无需使用欧拉角等参数化方法。解通过三角函数表达，并明确将旋转矩阵与庞索的运动几何图像联系起来，揭示了其可分解为关于角动量轴和实验室z轴的旋转。

ABSTRACT

The Euler-Poisson equations para determinar the rotation matrix of a rigid body can be solved without using of particular parameterization like the Euler angles. For the free Lagrange top, we obtain and discuss a general analytic solution, and compare it with the Poinsot's picture of motion.

研究动机与目标

• 推导自由刚体的欧拉-泊松方程的一般解析解，不依赖于欧拉角等参数化方法。
• 直接从运动方程表达旋转矩阵 $ R_{ij}(t) $ 的时间演化，保持完整的正交结构。
• 将解与庞索的运动几何图像联系起来，特别是惯性椭球在不变平面上的滚动。
• 证明旋转矩阵可从仅涉及角速度和一个基向量的约化系统解中重构。
• 提供一种全局有效、非局部的刚体转动描述，避免基于角度参数化固有的奇点。

提出的方法

• 通过拉格朗日作用量将刚体动力学表述为约束系统，导出关于 $ R_{ij}(t) $ 和 $ \boldsymbol{\Omega}(t) $ 的欧拉-泊松方程。
• 利用系统的运动积分——转动能量（3）、角动量（4）和正交性约束（5）——来限制解空间。
• 将完整系统约化为最小集合：求解欧拉方程（8）和一个基向量 $ \boldsymbol{\gamma}(t) $ 的矢量演化方程（9）。
• 通过求解三个初始条件为 $ (1,0,0), (0,1,0), (0,0,1) $ 的标准正交基向量 $ \mathbf{G}_1(t), \mathbf{G}_2(t), \mathbf{G}_3(t) $ 来构建完整的旋转矩阵 $ R(t) $。
• 以三角函数形式表达解，显式依赖于初始角动量分量和转动惯量。
• 将旋转矩阵分解为 $ R(t) = R_{\mathbf{m}}(t) \times R_{OZ}(t) $，表示关于角动量轴和实验室z轴的连续旋转。

实验结果

研究问题

• RQ1 能否在不引入欧拉角等参数化方法的情况下，直接求解自由刚体的欧拉-泊松方程以获得旋转矩阵？
• RQ2 旋转矩阵 $ R_{ij}(t) $ 的一般解如何与庞索对刚体运动的几何描述相关联？
• RQ3 拉格朗日陀螺的旋转矩阵在三角函数和物理常数下的显式形式是什么？
• RQ4 如何从仅涉及角速度和一个基向量的约化系统解中重构完整的旋转矩阵？
• RQ5 旋转矩阵分解 $ R(t) = R_{\mathbf{m}}(t) \times R_{OZ}(t) $ 在锥面运动和惯性椭球滚动方面的几何解释是什么？

主要发现

• 旋转矩阵 $ R_{ij}(t) $ 的一般解以显式三角函数形式导出，对所有时间有效且无奇点。
• 由于初始条件和 $ R_{ki}R_{kj} = \delta_{ij} $ 积分的守恒，解自动满足正交性约束 $ R^T R = I $。
• 旋转矩阵可分解为 $ R(t) = R_{\mathbf{m}}(t) \times R_{OZ}(t) $，其中 $ R_{\mathbf{m}}(t) $ 绕角动量矢量旋转，$ R_{OZ}(t) $ 绕实验室z轴以频率 $ \phi = \frac{(I_1 - I_3)m_3}{I_1 I_3} t $ 旋转。
• 瞬时角速度 $ \boldsymbol{\omega}(t) $ 被显式导出，其分量依赖于 $ m_2, m_3, I_1, I_3 $ 和一个与时间相关的余弦项。
• 角速度大小满足 $ \boldsymbol{\omega}^2 = k^2 + \frac{I_1 + I_3}{I_1 I_3} m_3 \phi $，其中 $ k = \frac{|\mathbf{m}|}{I_1} $，证实了复合旋转频率的存在。
• 该解与庞索的图像完全一致：体固定轴的端点在圆锥面上描绘出圆周轨迹，且惯性椭球在不变平面上无滑动地滚动。