高等动力学课件：lecture16-3D

1、 Lecture 163D刚体动力学内容回顾Lagrange情形三个首次积分章动角控制方程内容回顾三个根互不相等的情况（图b）21内容回顾存在三重根的条件存在两种情况：u1=u2=u3=1,如图(b)所示，这时只对应一种状态，即刚体只能绕竖直轴转动。Case1 : u2=u3时,当u2=u3,必有u2=u3=1f(u)代数表达式中也可能存在两个重根的情况。Case 2: u1=u20dotdot地球的进动与章动现象LLL=023.5春分点秋分点夏至冬至地球由于形状的非完全球形，与太阳和月亮的万有引力将会引起附件力矩，进而产生进动现象。Legendre多项式Legendre函数是指如下微分方程的

2、解勒让德方程的解在|x| 1 时可写成标准级数收敛的幂级数形式，当n = 0, 1, 2,. 时，随n 值变化方程的解相应变化，构成一组由正交多项式组成的多项式序列，这组多项式称为勒让德多项式Legendre 多项式具有如下级数的形式并且有多项式展开偏心力矩的表示：假设：地球的形状是一个标准的椭圆形状。太阳或月球认为是一个质点M问题：如何求出偏心力矩？如何分析地球在该偏心力矩作用下的运动？椭球体与一固定质点的万有引力势函数M,Crrimirii如图所示的质点mi与M质点之间的万有引力势可表示为偏心力矩的表示：对上式采用Legendre多项式展开后，表示为取Legendre 多项式的前三项近似对

3、第0项，即对第1次多项式M,Crrimirii对第二次多项式其中trJ是椭球体的惯量张量的迹，Jr是绕r轴的转动惯量势能函数表示因此，二阶近似的势能函数是假设地球绕质心的惯量椭球的主惯量为J1, J2, J3. 且J1=J2, , , 分别为r相对地球惯量椭球的方向余弦。势能函数表示xyzOrarccosOxy为地球运动的轨道平面，r为地球质心到太阳质心的连线的方向。O是地球的自转轴。O为地球的质心。三个方向余弦存在如下的关系则由二阶项引起的势能为假设地球绕太阳运动的轨道近似为圆轨道（r不变），为得到地球的平均进动角速度，二次项引起的势能在地球绕太阳一周的平均势能为地球的进动与章动现象地球的运

4、动等效为受到偏心万有引力作用下的刚体的定点运动情况。在平均意义上，地球做规则进动。则有地球的能量为考虑地球自转的角速度dot远远大于进动角速度，并且进动角速度通过平均化近似为常数，将能量对时间微分进一步简化为地球的进动与章动现象以上方程与Lagrange情形下出现重根的情况相同。令将考虑椭球性质的万有引力势的二阶近似函数代入根据Keppler第三定律，并设r为地球椭圆轨道的半长轴，椭圆轨道的圆频率为则进动角速度与椭圆轨道运动的角速度的比为地球自转轴的角度为在太阳万有引力作用下，进动完成一周所需要的时间大约为81,000年。同样的分析，由于月球对地球的万有引力作用也会引起进动，两者共同作用将导致

5、进动完成一周所需要的时间大约为26,000年。当然，在进动的同时，会有天章动现象。刚体一般运动的动力学方程两个坐标系: OxyzO,i,j,k惯性参考系，C,e1,e2,e3过质心的固连系 质心运动定理 对质心动量矩定理 与欧拉角的关系利用以上关系，可以得到利用Euler角和质心位置表示的刚体一般运动的6个二阶微分方程。在给定12个初始条件后，可以利用积分得到微分方程的解。均匀球体在重力场中的运动质心的运动：只受重力的作用。质心为一抛物线轨迹。姿态动力学分析：相对质心的动量矩守恒（Euler情形）。对球体来说，J1=J2=J3角速度矢量在空间中将是一固定矢量。与给定的初始条件相关。由于球体的对

6、称性，在任一方向的角速度均保持为常量。炮弹的稳定性轴对称炮弹在运动过程中，其对称轴将与质心运动方向发生一个偏转角。假设空气阻力的合力R作用在对称轴的A点上，并且方向与质心运动方向相同。问炮弹的自转角速度3应为多大，才能保证炮弹的弹头总是指向前方。炮弹的稳定性炮弹质心的运动：受空气阻尼力，并在重力场作用下质点的运动（见例4.12）,假设R是与速度平方的阻尼力，则轨迹方程为炮弹的稳定性对质心的动量矩定理:（Lagrange情形）轴对称刚体，可利用莱沙尔坐标系建立动力学方程炮弹的发射过程类似一个定轴转动的刚体的释放过程。初始章动角为零。要使章动角在外部阻力作用下，摄动为小量。应满足如下的条件。在慢进

7、动的情况下，可忽略p的高阶项要产生以上慢进动的条件为即要求自转的角速度必须足够高圆盘在粗糙面上的纯滚动Comments:只考虑圆盘纯滚动的情况，显然，由于接触点处的摩擦作用可能会改变圆盘接触的性质。其运动情况将会复杂的多。Euler将这一问题归结为力学中的12个难题之一。圆盘在粗糙面上的纯滚动圆盘的基本参数：半径 a; 质量m, 对质心处的转动惯量分别为：J1=J2, J3运动约束情况：接触点A处无相对滑动。坐标系统:惯性坐标系统： I, j, k质心平动坐标系统：C, x, y, z圆盘固连坐标系统：C , , 莱沙尔坐标系统: C n, s, 受力情况：重力mg, 接触点的约束力法向支撑力

8、Fz切向摩擦力Fn, Fm假设条件：接触点无相对滑动。不考虑滚动摩阻圆盘在粗糙面上的纯滚动各个坐标系之间的变换关系为将运动约束方程映射到O. n, m, k中矢量CA在盘面上，且指向接触点运动约束方程则定义圆盘在粗糙面上的纯滚动接触点处的运动约束方程在O,n,m,k中的表示为以上是定义在速度水平上的约束方程，且不可积，称为非完整约束方程。质心运动定理 ：并写为坐标系O, n, m, k中的分量形式圆盘在粗糙面上的纯滚动在莱沙尔坐标系中应用质心动量矩定理莱沙尔参考系的角速度外力矢量在莱沙尔坐标系下的分量形式质心动量矩方程在莱沙尔坐标系中的分量形式为圆盘在粗糙面上的纯滚动消去约束反力则则质心动量矩方程的第三式变为圆盘在粗糙面上的纯滚动根据约束方程则则质心动量矩方程的第一式变为圆盘稳态规则进动定倾角的运动：=0圆盘的动力学方程为由于则圆盘稳态规则进动三个常数之间的关系为能够成立的条件是根据运动约束方程，质心的运动规律为以上表明：质心的轨迹做匀速圆周运动，轨迹的半径为当0=/2时，则根据三常数条件3=0, 圆盘直立自转不滚p=0,圆盘在一个竖直平面内做纯滚动。摄动分析圆盘沿直线纯滚动的标称运动条件受到摄动后假设纯滚动条件仍然保持，略去高于2阶的小量。圆盘