第一章 二体问题

1、第一章 航天器运动主讲教师：杏建军2021年7月16日2授课内容v 绪论v 二体相对运动方程v 二体相对运动方程的求解31.1 课程的主要研究内容课程的主要研究内容主要研究内容主要研究内容：人造物体（航天器）在空间（距离地面100 km以上）自然力和人为控制力作用下运动的一门学科。课程名称：航天器动力学基础与应用课程名称：航天器动力学基础与应用主要包括天体引力和大气阻力航天器的运动航天器的运动：包括质心运动和姿态运动，相应的课程为航天器质心动力学和航天器姿态动力学。航天器质心运动和姿态运动是解耦的，因此可以分开研航天器质心运动和姿态运动是解耦的，因此可以分开研究和控制究和控制。为什么？41.1

2、 课程的主要研究内容课程的主要研究内容课程主要研究内容：课程主要研究内容：1、航天器运动数学建模（、航天器运动数学建模（物理规律的研究物理规律的研究）2、数学模型的求解（、数学模型的求解（微分方程的求解微分方程的求解）3、数学模型解算条件的提供（、数学模型解算条件的提供（微分方程解算条件的提供微分方程解算条件的提供）4、航天器运动规律的应用（、航天器运动规律的应用（具体的工程实践具体的工程实践）物理规律的研究：物理规律的研究：牛顿定理和万有引力定理牛顿定理和万有引力定理非线性微分方程的求解非线性微分方程的求解：数学模型解算条件的提供数学模型解算条件的提供：初始轨道的测量和最优控制：初始轨道的测

3、量和最优控制航天器运动规律的应用航天器运动规律的应用：各种航天任务的轨道设计：各种航天任务的轨道设计( , , , )ftrr r u51.2 课程的作用课程的作用航天器质心运动规律的研究是一切航天任务的开始的基础航天器质心运动规律的研究是一切航天任务的开始的基础任何航天任务的设计都要满足其基本的运动规律（例如：任何航天任务的设计都要满足其基本的运动规律（例如：在万有引力作用下，航天器作圆锥曲线运动在万有引力作用下，航天器作圆锥曲线运动）讨论问题：有没有可能设计任意飞行的航天器？讨论问题：有没有可能设计任意飞行的航天器？61.2 课程的作用课程的作用航天器相对运动实例航天器相对运动实例71.2

4、 课程的作用课程的作用对于人造地球卫星来说，航天器轨道高度是任务设计的对于人造地球卫星来说，航天器轨道高度是任务设计的关关键参数键参数。为什么？轨道高度决定了：轨道高度决定了：发射成本发射成本效载荷规模。如雷达、光学相机、通信卫星发射机效载荷规模。如雷达、光学相机、通信卫星发射机对地球的覆盖范围。对地球的覆盖范围。对热点地区的覆盖特性对热点地区的覆盖特性决定了一些有特殊用途的轨道决定了一些有特殊用途的轨道81.2 课程的作用课程的作用对于深空探测器，航天器轨道设计决定了整个任务过程对于深空探测器，航天器轨道设计决定了整个任务过程卡西尼号是17国参与的土星探测任务，历时过6年8个月、32亿千米。

5、卡西尼号的诡异飞行轨迹卡西尼号的诡异飞行轨迹为什么这样设计？直接飞行只需12.5亿千米行星助力飞行，节省燃料。如果直接飞往土星，需要70吨推进剂，卡西尼号总重6.4吨91.2 课程的作用课程的作用亚星三号卫星亚星三号卫星 亚星三号是美国休斯公司为香港亚洲卫星公司制作的通信卫星，于1997年12月由质子火箭发射进入地球静止轨道，但由于火箭故障，进入了轨道倾角为51的无用轨道，发射失败，香港卫星公司向保险公司索赔2亿美金。一个航天史上的故事由此开始一个航天史上的故事由此开始101.3 学科发展史学科发展史从严格意义上来说，航天器动力学开始于上世纪从严格意义上来说，航天器动力学开始于上世纪50年代前

6、苏年代前苏联发射第一颗人造地球卫星，但它的起源非常久远。联发射第一颗人造地球卫星，但它的起源非常久远。古典天文学古典天文学中国是世界上天文学起步最早，发展最快的国家之一。中国是世界上天文学起步最早，发展最快的国家之一。早在尧舜时代就设置了天文官。早在尧舜时代就设置了天文官。1. 古希腊也是古典天文学什么发达的国家，取得了辉煌的古希腊也是古典天文学什么发达的国家，取得了辉煌的成就，并对欧洲各国的文化影响很大。主要成果包括：成就，并对欧洲各国的文化影响很大。主要成果包括：确定地球的形状和大小；日月的远近和大小；日心说等。确定地球的形状和大小；日月的远近和大小；日心说等。111.3 学科发展史学科发

7、展史古典天文学的社会需求古典天文学的社会需求制定历法的需要，知道农业生产。如我们常说的制定历法的需要，知道农业生产。如我们常说的24节气节气预测天灾人祸，旦夕祸福。（在古代，占星术和天文学预测天灾人祸，旦夕祸福。（在古代，占星术和天文学是没有明显的区别的）是没有明显的区别的）古典天文学研究方法古典天文学研究方法没有理论指导，没有先进的观测手段没有理论指导，没有先进的观测手段兴趣，长期不懈的观测，积极思考兴趣，长期不懈的观测，积极思考121.3 学科发展史学科发展史天体力学：应用天体力学：应用力学规律力学规律研究天体的运动和形状研究天体的运动和形状 天体力学以天体力学以数学数学为主要研究手段（微

8、积分），以牛顿为主要研究手段（微积分），以牛顿万万有引力定律有引力定律为基础。为基础。天体力学的发展天体力学的发展奠基期奠基期（从古典天文学到十九世纪后期），标志性成果：（从古典天文学到十九世纪后期），标志性成果：开普勒提出三大定理；牛顿创立微积分，发现万有引力开普勒提出三大定理；牛顿创立微积分，发现万有引力定理；欧拉、达朗贝尔、拉格朗日、拉普拉斯创立分析定理；欧拉、达朗贝尔、拉格朗日、拉普拉斯创立分析力学，建立了天体力学的力学基础，提出了摄动理论的力学，建立了天体力学的力学基础，提出了摄动理论的分析方法；海王星的发现（理论的实际应用）。分析方法；海王星的发现（理论的实际应用）。131.3 学

9、科发展史学科发展史天体力学的发展天体力学的发展发展期发展期（从十九世纪后期到二十世纪（从十九世纪后期到二十世纪50年代）。研究对年代）。研究对象新增加了太阳系的小天体。研究方法新增了定性方法象新增加了太阳系的小天体。研究方法新增了定性方法和数值方法，定性方法由庞加莱和李雅普洛夫创建，数和数值方法，定性方法由庞加莱和李雅普洛夫创建，数值方法最早追溯到高斯（最小二乘定轨）。值方法最早追溯到高斯（最小二乘定轨）。新时期新时期（二十世纪（二十世纪50年代以后，航天器动力学出现）。年代以后，航天器动力学出现）。研究对象新增了人造物体，物体运动的预报精度与观测研究对象新增了人造物体，物体运动的预报精度与观

10、测精度大大提高，相应的摄动分析方法、定性方法和数值精度大大提高，相应的摄动分析方法、定性方法和数值方法也有了相应的发展。方法也有了相应的发展。141.3 学科发展史学科发展史航天器动力学航天器动力学二十世纪二十世纪50年代以后，随着人造天体的发射，航天器动年代以后，随着人造天体的发射，航天器动力学出现。力学出现。与天体力学相比，研究对象发生了变化（人造物体）。与天体力学相比，研究对象发生了变化（人造物体）。与自然天体相比，人造物体的受力物体增加了人为控制力，与自然天体相比，人造物体的受力物体增加了人为控制力，运动形式更为复杂；物体的预报精度与观测精度大大提高。运动形式更为复杂；物体的预报精度与

11、观测精度大大提高。研究的基本方法没有变化：以微积分为数学基础，以摄研究的基本方法没有变化：以微积分为数学基础，以摄动分析方法、定性分析方法和数值方法为手段。动分析方法、定性分析方法和数值方法为手段。151.4 教程和参考书教程和参考书1、航天器轨道动力学，、航天器轨道动力学，赵钧赵钧编著，哈工大出版社，编著，哈工大出版社，20112、航天器轨道动力学与控制，、航天器轨道动力学与控制，杨嘉摨杨嘉摨主编，宇航出版社，主编，宇航出版社，1995（注：国内航天器领域经典专著）（注：国内航天器领域经典专著）3、Fundamentals of Astrodynamics and Applications(

12、Second Edition)，Vallado，D.V. Microcosm Press, 2001 （注：国外（注：国外经典教材）经典教材）4、An introduction to the mathematics and methods of astrodynamics, Richard H. B. AIAA, 1999 （注：（注：MIT教材）教材）161.5 考核方式和成绩评定考核方式和成绩评定考核方式考核方式考核内容考核内容成绩比例（成绩比例（%）平时到课率、课堂回平时到课率、课堂回答问题及研讨答问题及研讨基础知识，学习主动性基础知识，学习主动性2020课后作业课后作业综合应用知识解决

13、具体综合应用知识解决具体工程问题的能力工程问题的能力1010文献阅读与专题报告文献阅读与专题报告自主学习，分析问题和自主学习，分析问题和主动交流的能力主动交流的能力2020期末闭卷理论考试期末闭卷理论考试学生掌握基本概念及基学生掌握基本概念及基本理论的程度本理论的程度505017授课内容v 绪论v 二体相对运动方程v 二体相对运动方程的求解18牛顿第二定理牛顿第二定理：Force = Mass Acceleration 21211122Gm mdmrrdt2rrr2.1 万有引力定律和牛顿第二定律万有引力定律和牛顿第二定律21212222Gm mdmrrdt2rrr万有引力定理万有引力定理：任

14、意两个质点有通过连线方向上的力相互吸引。该引力大小与它们质量的乘积成正比，与它们距离的平方成反比。m1m2r1r21921211122Gm mdmrrdt2rrr2.2 二体系统线动量守恒定理二体系统线动量守恒定理21212222Gm mdmrrdt2rrr21 12 220dmmdtrr1 12 212mmmmcmrrr220ddtcmr内力不改变系统的质心内力不改变系统的质心2021211122Gm mdmrrdt2rrr2.3 二体相对运动方程二体相对运动方程21212222Gm mdmrrdt2rrr2121222() ()()G mmdrrdt21rrrr212223()=G mmd

15、dtrrr rrr212.3 二体相对运动方程二体相对运动方程121 1221220Gmmmam ammr12aa另一种推导思路另一种推导思路质点质点m1与质点与质点m2所受的相互引力大小相等，方向相反所受的相互引力大小相等，方向相反212231()mG mmmrrr 2122raaaar2222222ddddtdtdt121rrraa222.4 讨论讨论1：m1相对系统质心的相对运动方程相对系统质心的相对运动方程推导质点推导质点m1相对系统质心的相对系统质心的相对运动方程相对运动方程212223()=G mmddtrrr rrr质点质点m2相对相对m1的相对运动方程的相对运动方程m1m2r1

16、r2rcmcm32222312()Gmddtmmr rr231.4 讨论问题讨论问题2如果在地球附近发射一颗人造地球卫星，其相对地心的如果在地球附近发射一颗人造地球卫星，其相对地心的相对运动方程是什么？相对运动方程是什么？如果向月球发射一个月球探测器，其相对地月质心的相如果向月球发射一个月球探测器，其相对地月质心的相对运动方程是什么？对运动方程是什么？如果向月球发射一个月球探测器，其相对地球质心的相如果向月球发射一个月球探测器，其相对地球质心的相对运动方程是什么？对运动方程是什么？242.5 二体相对运动方程的标量表示二体相对运动方程的标量表示矢量形式的非线性微分方程矢量形式的非线性微分方程2

17、12223()=G mmddtrrr rrr为了在实际中应用，需要转换为标量形式为了在实际中应用，需要转换为标量形式yzxyzxriii223223223=d xxdtrd yydtrd zzdtr 直角坐标的基矢量不直角坐标的基矢量不随时间变化随时间变化252.6 讨论讨论3： 二体相对运动方程的极坐标形式二体相对运动方程的极坐标形式极坐标的基矢量是随时间变化极坐标的基矢量是随时间变化xiyirii222220drrdtrddrrdtdt 262.7 课后作业课后作业通过复习或者自学理论力学中拉格朗日方程的有关内容，通过复习或者自学理论力学中拉格朗日方程的有关内容，应用拉格朗日方程推导讨论问

18、题应用拉格朗日方程推导讨论问题1、3。下节课随机抽点。下节课随机抽点学生上台上述，并请所有学生提交推导过程作业。学生上台上述，并请所有学生提交推导过程作业。作业关键点：作业关键点：拉格朗日函数拉格朗日函数 = 系统动能系统动能-系统的势能；系统的势能；拉格朗日方程；拉格朗日方程；1. 广义坐标。广义坐标。阅读：阅读：An introduction to the mathematics and methods of astrodynamicsp97-p99；p110。27授课内容v 绪论v 二体相对运动方程v 二体相对运动方程的求解283.1 近似方法近似方法非线性微分方程非线性微分方程223=

19、ddtrrr微分方程解决的基本方法：解析法和数值法微分方程解决的基本方法：解析法和数值法解析法：求出非线性微分方程的解析解。解析法：求出非线性微分方程的解析解。hg 算例：自由落体运动算例：自由落体运动293.1 近似的方法近似的方法基本原理：泰勒级数展开基本原理：泰勒级数展开数值积分法：通过数值积分的方法求解数值积分法：通过数值积分的方法求解(,)yfxy1(1)( )y iy iKhh步长步长1( ),( )(,)xiyiKfxy经典方法：龙格库塔方法经典方法：龙格库塔方法303.1 近似的方法近似的方法龙哥库塔方法龙哥库塔方法Matlab函数：函数：ode45t,y=ode45(orbi

20、t,tspan,y0,options)function dy = orbit(t,y,flag,mu)%函数说明输入输出函数说明输入输出dy =zeros(6,1);r = sqrt(y(1)2 + y(3)2 + y(5)2);dy(1) = y(2);dy(2) = -mu/r3*y(1);dy(3) = y(4);dy(4) = -mu/r3*y(3);dy(5) = y(6);dy(6) = -mu/r3*y(5);end223223223=d xxdtrd yydtrd zzdtr 313.1 近似的方法近似的方法function robit_computermu = 3.98600

21、4418e+14; %地球引力常数tspan = 0:60:86400;options = odeset(AbsTol,1e-15,RelTol,1e-12,NormControl,on);x0 = -5292392.072;-4862.201380;3111662.355; -4136.781314; 3101114.660;-4147.028008;T,Y = ode45(orbit,tspan,x0,options,mu);plot3(Y(:,1),Y(:,3),Y(:,5)32-1-0.500.51x 107-505x 106-505x 106333.1 近似的方法近似的方法近似级数展

22、开方法近似级数展开方法220002()( )( )( )()()2 !ttd r tdr tr tr tttd td t An introduction to the mathematics and methods of astrodynamicsp110-p114有二体问题的近似级数展开详细推导和高效迭代算法，请大家自学并完成问题3-5。343.2 二体问题的解析解二体问题的解析解二体运动方程二体运动方程223=ddtrrr角动量矢量常数角动量矢量常数223=ddtr rrrr()=0dddtdtvrvrconst hrv353.2 二体问题的解析解二体问题的解析解开普勒第二定理开普勒第二定

23、理二体问题角动量是常数二体问题角动量是常数角动量在极坐标下的表示角动量在极坐标下的表示dddrrrrdtdtdtrrrirviii2drdt zhrvi面积化率面积化率212dAdrdtdt363.2 二体问题的解析解二体问题的解析解偏心率矢量常数偏心率矢量常数223=ddtrrhrh2zhr riidddtdt rii拉格朗日矢量拉格朗日矢量()dddtdtrivhconstrvhree偏心率矢量偏心率矢量e373.2 二体问题的解析解二体问题的解析解开普勒第一定理开普勒第一定理re rvh rr r()abc矢量运算公式矢量运算公式2hrre rvh rrv hr rcosrefe r其中

24、其中f是是e与与r两个矢量的夹角，定义为两个矢量的夹角，定义为真近点角真近点角383.2 二体问题的解析解二体问题的解析解2cosrefhr2/1cos1coshprefef圆锥曲线方程，当圆锥曲线方程，当e在（在（0，1）之间时为椭圆，即开）之间时为椭圆，即开普勒第一定理普勒第一定理 f= 0为近地点，为近地点，f = 180为远地点为远地点11papprree22(1)parrapae393.2 二体问题的解析解二体问题的解析解开普勒第三定理开普勒第三定理2(1)1cos1cospaerefef椭圆的面积椭圆的面积22201122dAr dAr d403.2 二体问题的解析解二体问题的解析

25、解行星（人造地球卫星）扫过一圈的时间定义为行星行星（人造地球卫星）扫过一圈的时间定义为行星的运动周期的运动周期T222(1)1222paeaehT行星的运动周期行星的运动周期32aT定义椭圆上的平均运动定义椭圆上的平均运动n32nTa开普勒第三定理开普勒第三定理32aconstT这个定理是正确的吗？413.2 二体问题的解析解二体问题的解析解能力积分常数能力积分常数11revhr 222222121121121rrh vhre evhvhrvhhvhvrv h12222221hvvearr423.2 二体问题的解析解二体问题的解析解能量积分常数能量积分常数122222vvarra 动能势能动能动能+势能势能=常数常数 在中心引力作用系统能力守恒在中心引力作用系统能力守恒此公式又成为此公式又成为活力公式活力公式（vis-viva equations）,是航天是航天动力学里面非常重要的一个恒等式。动力学里面非常重要的一个恒等式。433.2 二体问题的解析解二体问题的解析解活力公式的应用，航天器速度的计算活力公式的应用，航天器速度的计算222vra讨论讨论