660 likes | 859 Vues
第三章 航天器姿态运动学和动力学. 3.1 航天器的姿态运动学. 3.2 航天器的姿态动力学. 3.3 航天器的一般运动方程. 3.4 姿态干扰力矩. 第三章 天器的姿态运动学和动力学. 航天器的 姿态运动学 是从几何学的观点来研究航天器的运动,它只讨论航天器运动的几何性质,不涉及产生运动和改变运动的原因;而航天器的 姿态动力学 则是研究航天器绕其质心运动的状态和性质。所以航天器姿态的运动方程须由两部分组成,一部分为通过坐标变换关系得出的运动学方程,另一部分则是以牛顿动力学定律 ( 如动量矩定律 ) 为基础的动力学方程。
E N D
第三章 航天器姿态运动学和动力学 3.1 航天器的姿态运动学 3.2 航天器的姿态动力学 3.3 航天器的一般运动方程 3.4 姿态干扰力矩
第三章 天器的姿态运动学和动力学 航天器的姿态运动学是从几何学的观点来研究航天器的运动,它只讨论航天器运动的几何性质,不涉及产生运动和改变运动的原因;而航天器的姿态动力学则是研究航天器绕其质心运动的状态和性质。所以航天器姿态的运动方程须由两部分组成,一部分为通过坐标变换关系得出的运动学方程,另一部分则是以牛顿动力学定律(如动量矩定律)为基础的动力学方程。 本章中将航天器视作刚体。
3.1 航天器的姿态运动学 3.1.1常用参考坐标系 坐标系形式很多,每种坐标系都有其自己的特点,因此也就只适用于一定的范围,所以根据具体情况选择坐标系是必要的。一般来说,讨论航天器姿态运动常用的坐标系,主要有4种。
1.惯性坐标系 所有的运动都要参照的基本坐标系是惯性坐标系, 2.质心平动坐标系 这是一个与惯性坐标系密切相关的坐标系。原点O位于航天器质心,OX,OY,OZ轴分别与某一惯性坐标系的坐标轴保持平行。 3.质心轨道坐标系 简称轨道坐标系。这是一个以航天器质心为原点的正交坐标系,如图3.1所示。
4.本体坐标系Oxyz 又称为星体坐标系。在此坐标系中,原点0在航天器质心,Ox,Oy,Oz三轴固定在航天器本体上。若Ox,Oy,Oz三轴为航天器的惯量主轴,则该坐标系称为主轴坐标系。
3.1.2 航天器的姿态运动学方程 在坐标系确定以后,航天器上任何一点的位置就可以在固联于星体的本体坐标系Oxyz中表示;若要描述三轴稳定航天器的对地定向运动,则要借助于质心轨道坐标系 ;若要讨论自旋卫星的章动运动时,就必须运用质心平动坐标系OXYZ。而各种坐标系之间的关系可以通过一系列旋转角来表示,这些旋转角称为欧拉角。具体地说可以通过3个欧拉角 , , 来确定本体坐标系Oxyz相对于其他坐标系的位置。
以坐标系Oxyz和OXYZ为例,星体轴的位置可通过3次旋转达到OXYZ坐标轴的位置。旋转顺序具有多种形式,但不能绕一个轴连续旋转两次,因为连续两次旋转等同于绕这个轴的一次旋转。为此可以得出两类12种可能的旋转顺序如下:以坐标系Oxyz和OXYZ为例,星体轴的位置可通过3次旋转达到OXYZ坐标轴的位置。旋转顺序具有多种形式,但不能绕一个轴连续旋转两次,因为连续两次旋转等同于绕这个轴的一次旋转。为此可以得出两类12种可能的旋转顺序如下: 一类:1-2-3, l-3-2,2-3-1,2-1-3,3-1-2,3-2-1; 二类:3-1-3, 2-l-2,1-2-1,3-2-3,2-3-2,1-3-1。 显然,一类是每轴仅旋转一次,二类是某一轴不连续地旋转两次。下面详细介绍被称为经典欧拉转动顺序的“3-1-3”旋转和“1-2-3”旋转。
1.“3-1-3”旋转 (1)OXYZ一绕OZ (“3”)轴转 角 :如图3.2所示,这两个坐标系之间的变换矩阵为 (3.1)
(2) 绕 (“1”)轴转 角 :如图3.3所示,这两个坐标系之间的变换矩阵为 (3.2)
(3) 绕 (“3”)轴转 角 :如图3.4所示,这是最后一次旋转,此时已达到了航天器的本体坐标系Oxyz。两者的变换矩阵可推导为 (3.3)
综合以上变换,坐标系OXYZ与Oxyz之间的直接转换关系即为综合以上变换,坐标系OXYZ与Oxyz之间的直接转换关系即为
若令 ,则通过A可以把质心平动坐标系OXYZ中表示的矢量分量变换成为本体坐标系Oxyz中表示的分量,即 (3.4)
若坐标系Ozyz中的分量已知,需要确定坐标系OXYZ中的分量,则利用两个坐标系之间正交变换的逆矩阵就等于它的转置矩阵这一性质,即若坐标系Ozyz中的分量已知,需要确定坐标系OXYZ中的分量,则利用两个坐标系之间正交变换的逆矩阵就等于它的转置矩阵这一性质,即 得到 (3.5)
其中 (3.6) (3.7)
这样,利用经典欧拉转动,通过 3个欧拉角就将航天器的本体坐标系Oxyz和质心平动坐标系相互联系起来了。 基于欧拉转动顺序”3-1-3”,可以进一步将航天器的空间转动角速度ω在本体坐标系中的分量 用欧拉角表示,从而推导出航天器的姿态运动学方程。
如图3.5所示。将角速度 沿 和 轴分解,则 , 和 在正交坐标系 中的分量分别为: 轴为 , 轴为 , 轴为 。再将 轴和 轴分量按Ox和Oy轴分解,其结果表示如下: (3.8)
或者以逆形式表示,即 (3.9) 式(3.8)或(3.9)即为航天器的一组姿态运动学 方程。
2.“1-2-3”旋转 类似地,也可以通过欧拉“1-2-3”旋转将航天器的不同坐标系相互联系起来。例如从 出发,进行以下3次旋转: (1) 绕 (“l”)转 角 (2) 绕 (“2”)转 角 (3) 绕 (“3”)转 角 于是坐标系Oxyz和 之间的坐标变换关系即为 (3.10)
(3.11) 式中 (3.12)
同样可得按照2-3-1,3-1-2,1-3-2,2-1-3,3-2-1等不同转动顺序的变换关系。当 时,即在小角度变化情况下, 可近似为 (3.13) 其中欧拉角 分别称为俯仰角、偏航角和滚动角,而Oz,oy,Oz轴分别称为航天器的滚动轴、俯仰轴和偏航轴。
相应地,利用“l-2-3”姿态角也可以将 的分量 表示出来,得到另一组航天器的姿态运动学方程,即 (3.14) 或者以逆形式表示为 (3.15)
3.2航天器的姿态动力学 作为刚体的航天器的姿态动力学是以刚体的动量矩定理为基础的。因此在确定了描述航天器姿态运动的各种坐标系和运动学之后,了解刚体的动量矩定理就成为研究航天器姿态动力学的一个重要条件。
3.2.1动量矩定理 首先考察质点,如图3.6所示,力 对点 的矩 (3.16) 其中矢径 ,且A在力的作用线上。因此,力矩矢量 ,垂直于由 和 作用线组成的平面,并且 的指向按右手规则来确定。类似地,质点的动量 对点0的矩可表示成 (3.17) 它垂直于质点的矢径 和动量 所组成的平面,且 的指向也由右手规则确定。
静力学里曾指出,力对于通过点O的任一轴,例如Oz轴静力学里曾指出,力对于通过点O的任一轴,例如Oz轴 的矩,等于它对点O的矩在该轴上的投影 ,并且可以写成 = 该动量矩具有量纲 在国际单位制中,动量矩的常用单位是 。
设坐标系Ozyz是固定直角坐标系,以矢径r与牛顿第二定律的方程作叉乘,有设坐标系Ozyz是固定直角坐标系,以矢径r与牛顿第二定律的方程作叉乘,有 等号右端就是力F对原点O的矩 ,左端可以改造为 但 ,所以上式等号右端第二项等于零(两个平行矢量的叉积等于零),而第一项就是质点对点O的动量矩矢量 对时间的导数。于是得
(3.18) 即质点对任意固定点的动量矩对时间的导数,等于该质点所受的力对同一点的矩。这就是质点的动量矩定理。 若 =O,则 =常矢量。即若质点所受的合力对某固定点的矩恒等于零,则 质点对同一点的动量矩守恒。该结论说明了质点动量矩守恒的条件。 动量矩定理很容易由质点推广到质点系。按式(3.18)对质点系内每个质点写出动量矩方程,然后相加,得
其中末等号左端方括号中就是整个质点系对固定点O的动量矩,用Ho代表,即其中末等号左端方括号中就是整个质点系对固定点O的动量矩,用Ho代表,即 等号右端等于质点系所受合外力对点O之矩的矢量和,用Mo代表。内力成对地出现,它们对任一点之矩的矢量和恒等于零。于是有 (3.19) 可见,质点系对任一固定点的动量矩对时间的导数,等于该质点系所受全体外力对同一点之矩的矢量代数和。这就是质点系动量矩定理。 特殊情况:若 ,则Ho =常矢量。
3.2.2姿态动力学方程 设航天器在空间以角速度 旋转,其动量矩为Ho。为了方便起见,基准点选航天器本体坐标系Oxyz的原点,也即航天器质心0,M是作用在航天器相对于质心0的合外力矩,所以航天器的动量矩即为 (3.20) 式中,矢量r是刚体内相对于质心的矢径;dr/dt是质量元dm在空间相对于质心的速度矢量;m为航天器的总质量。于是在本体坐标系中,刚体的 和M可以分别表示成
(3.21) (3.22) (3.23) (3.24) 式中, 是航天器本体坐标系各轴的单位矢量,上两式右端的系数则是相应矢量沿各坐标轴的分量。将式(3.21)对时间t求取导数,求动量矩H在空间的变化率,即 (3.25) 由于刚体在空间中以的角速度进行旋转,所以与其固连的本体坐标系各轴方向也在相应变化。
以知坐标轴单位矢量的导数公式是 (3.26) 代入式(3.25),并根据动量矩定理得 (3.27) 因 所以式(3.27)在航天器本体坐标系中可以展开为
其在各轴的分量表示为 (3.29a) 或表示成矩阵矢量形式,即 (3.29b) 式(3.29a)或(3.29b)称为欧拉力矩方程式。
同理,对式(3.23)求导也可得 若刚体内各质点相对于质心的位置不变,式(3.20)描述的动量矩即为 (3.30) 利用矢量叉乘公式,有 代入(3.30),并考虑到式(3.22),则
(3.31a) 即 (3.13b) 式中,I为惯性矩阵;Ix,Iy,Iz分别为刚体绕坐标轴Ox,Oy,Oz的转动惯量; 称为惯量积。它们分别为
惯量积的数值可正可负,它们与坐标系的选取密切有关。如果在某一坐标系中, ,则该坐标系称为主轴坐标系,OX,Oy,Oz轴就是刚体的主惯量轴。 因此,如果取航天器的本体坐标系为主轴坐标系,则有 (3.32) 把它们代人欧拉力矩方程(3.29),并忽略质量变化就可以以得到 (3.33) 这就是基于本体坐标系的航天器的姿态动力学方程组,也称为欧拉动力学方程。
3.3 航天器的一般运动方程 3.3.1六自由度运动方程 设作为刚体的航天器质量为m,质心为O,坐标系是质心轨道坐标系,坐标系Oxyz是本体坐标系,且坐标轴Ox,Oy,Oz取为航天器主惯量轴,坐标系是惯性坐标系,F是所有作用在航天器上的合外力矢量,M是所有作用在航天器上相对于O点的合外力矩矢量。 根据牛顿第二定律,相对于质心O的动力学方程在惯性坐标系中的投影式为 (3.34) 式中, 为F在坐标系 各轴上的投影分量。
实际上,式(3.34)当中的由式(2.7)、(2.8)和(2.9)表示,即实际上,式(3.34)当中的由式(2.7)、(2.8)和(2.9)表示,即 (2.8) (2.7) (2.9) 而在第二章中讨论的二体轨道运动方程式(2.21)正是式(3.34)在以下特殊条件下的极坐标形式: (1)式(2.7)中n=2; (2)式(2.8)中 = 0。
又根据对质心的动量矩定理,航天器绕质心O运动的姿态动力学方程在本体坐标系Oxyz中的投影式为又根据对质心的动量矩定理,航天器绕质心O运动的姿态动力学方程在本体坐标系Oxyz中的投影式为 (3.33)
式中, 是M沿航天器主惯量轴的分量; 是航天器空间转动角速度 沿主惯量轴的分量,它们与欧拉角 的关系是 (3.15) 联立式(3.34)、(3.15)和(3.33)三组方程就得到了刚性航天器一般运动的全部运动