月面大范围空间机动初步优化_任宽.pdf
《月面大范围空间机动初步优化_任宽.pdf》由会员分享,可在线阅读,更多相关《月面大范围空间机动初步优化_任宽.pdf(7页珍藏版)》请在咨信网上搜索。
1、2023年第1期 导 弹 与 航 天 运 载 技 术(中英文)No.1 2023 总第392期 MISSILES AND SPACE VEHICLES Sum No.392 收稿日期:2020-09-18;修回日期:2021-01-23 文章编号:2097-1974(2023)01-0052-07 DOI:10.7654/j.issn.2097-1974.20230111 月面大范围空间机动初步优化 任 宽,申 麟,李 扬,吴胜宝,胡冬生(中国运载火箭技术研究院研究发展部,北京,100076)摘要:月面大范围机动技术将是未来探索以及开发月球的关键技术之一,月面机动轨迹优化可以显著减少机动过程的
2、燃料消耗。主要针对月面大范围机动轨迹进行初步优化,将整个机动轨迹分为上升段、自由漂浮段、第1次动力减速段、垂直动力着陆段,前面3个飞行阶段是对整个飞行任务的初制导,在达到飞行射程的同时使飞行终端满足垂直着陆段的降落要求,通过垂直着陆段修正初制导误差,使飞行器满足目标点精确着陆要求。首先,建立各个飞行阶段的运动学模型,将前3个飞行阶段作为一个轨迹优化问题,采用遗传优化算法对优化问题进行求解。另外,将垂直动力着陆段轨迹优化问题通过无损凸化和离散的方式转化为一个有限维的二阶锥凸问题,通过凸优化求解器求解以实现在线轨迹优化。仿真结果表明,采用的轨迹优化方法具有较高的任务适应性和鲁棒性,可以适应不同机动
3、任务需求实现高精度垂直动力着陆。关键词:月面大范围机动技术;机动轨迹优化;垂直动力着陆;凸优化 中图分类号:V476.3 文献标识码:A Initial Optimization of Lunar Large-scale Space Maneuver Ren Kuan,Shen Lin,Li Yang,Wu Sheng-bao,Hu Dong-sheng(R&D Department,China Academy of Launch Vehicle Technology,100076,Beijing)Abstract:The large-scale lunar maneuvering techn
4、ology will be one of the key technologies for future exploration and development of the moon.The optimization of lunar maneuvering trajectory can significantly reduce fuel consumption during the maneuvering process.The initial optimization of the lunar large-scale maneuvering trajectory is focused o
5、n,and the entire maneuvering trajectory is divided into ascending section,free-floating section,first dynamic deceleration section,vertical dynamic landing section.The first three flight sections are the initial guidance for the flight mission.While reaching the flight range,the flight terminal meet
6、s the landing requirements of the vertical power landing section,and the vertical power landing section is to correct the initial guidance error so that the aircraft can accurately land at the target point.Firstly,the kinetic model of each flight section,the first three flight sections is establishe
7、d are considered as a trajectory optimization problem,and genetic optimization algorithm is used to solve the optimization problem.In addition,the trajectory optimization problem of the vertical dynamic landing section is transformed into a finite-dimensional SOCP(second-order cone program)problem t
8、hrough lossless convexification,and solved by the convex optimization problem solver to achieve online trajectory optimization.The simulation results show that the trajectory optimization method adopted has high task adaptability and robustness,and can achieve high-precision vertical power landing a
9、dapting to different maneuver task requirements.Key words:large-scale lunar maneuvering technology;optimization of lunar maneuvering trajectory;vertical power landing;convex optimization 0 引 言 月球一直是载人探索最热门目标星球之一,美国在2017年正式宣布将重启登月计划并建立永久性月球基地,各国也纷纷将载人登月与月球开发作为未来航天发展的重点1,2。从人类发展角度来说,探索与开发月球是人类走出地球摇篮的必
10、经之路,因此,研究月面大范围机动技术对未来开发月球、建立月球基地、探索月面空间运输及地月空间运输技术具有重要意义。月面空间运输与地球空间运输有较大差别,月球引力只有地球的六分之一,同时,月球没有大气,不用考虑气动、气压对空间运输的影响。月面大范围机动轨迹设计与传统弹道设计有一定区别,需要考虑动力减速的过程,飞行任务终端需要满足较为严格的终 任 宽等 月面大范围空间机动初步优化 53第1期 端约束,包括位置和速度约束,同时,通常需要采用在线制导的方式以实现精确动力软着陆。由于整个月面空间机动将会经历多个飞行阶段,所以月面机动轨迹优化应该是一个整体优化的过程,各个飞行阶段相互制约,相互影响,并且需
11、要适应不同发射任务。本文主要研究中远射程(L500 km)的月面机动轨迹优化,将整个机动轨迹分为 4 个阶段,分别是上升段、自由漂浮段、第 1 次动力减速段、垂直动力着陆段。根据月面机动各飞行阶段特性建立其动力学模型,将整个机动过程轨迹优化问题分为 2 个轨迹优化问题,分别是月面机动任务轨迹优化问题和终端垂直动力段在线轨迹优化问题。随后,针对几种飞行任务的机动轨迹优化问题进行了数值分析,同时考虑垂直着陆段初始状态偏差进行蒙德卡洛仿真分析,验证了轨迹优化算法的有效性及多任务适应性。1 月面机动动力学模型 整个月面大范围机动轨迹如图 1 所示。图1 月面机动轨迹示意 Fig.1 Schematic
12、 Diagram of Lunar Maneuvering Trajectory 1.1 月面机动参考坐标系 为便于各个飞行阶段的动力学建模,本文采用图 2所示的 3 组坐标系。平春分点平赤道面平赤经平赤纬ZMYMOMXMyOxzOTXTYTZT 图2 月面机动坐标系示意 Fig.2 Schematic Diagram of Lunar Maneuvering Coordinate System a)月心惯性坐标系MMMMO X Y Z:原点MO位于月心,参考面与月球平赤道面重合,MMO X轴指向平春分点,MMO Z轴指向平天极,MMO Y轴由右手法则确定。b)发射坐标系Oxyz:原点O位于发
13、射点,Ox轴平行于地面指向目标方向,Oy轴垂直于地面向上,Oz轴与x、y轴构成右手坐标系。c)目标点坐标系TTTTO X Y Z:原点TO位于目标降落点,TTO X垂直于目标点水平面指向上方,TTO Z位于水平面内指向月球北极,TTO Y轴由右手法则确定。考虑到月球自转周期较长,忽略月球自转的影响,假设月球为匀质圆球,不考虑月球摄动以及地球引力的影响,假设飞行器在飞行过程中处于瞬时平衡状态。1.2 上升段 在月面大范围机动轨迹设计过程中,上升段是飞行器最初动力加速的阶段,上升段的终端参数决定了整个机动轨迹的射程。控制飞行器在发射平面内飞行,忽略侧向运动,建立发射坐标系下月面机动的动力学方程为
14、3m3specos()sin()()xyxyxvyvPvtxmRPvtyRmRPmI g=-=-+=?(1)式中 x,y,xv,yv分别为飞行器的位置和速度分量;R为月心到飞行器质心的距离;月球引力常数;mR为月球半径;P为发动机推力;()t为相对于发射坐标系的俯仰角;spI为发动机海平面比冲;eg为地球重力加速度。对于上升段轨迹,需设计的是俯仰程序角,一般将上升段轨迹分为垂直上升段、程序转弯段、瞄准段。设计俯仰程序角表达式为 111ss100sk ,2()(),(/22 ,tttttttt ttttt=-=+-?(2)式中 0为终端俯仰角;?为俯仰角速率,可以设计为一个常值。上升段需要优化的
15、控制变量是终端俯仰角0和上升段工作时间kt。1.3 自由漂浮段 自由漂浮段是继上升段结束、无动力漂浮的飞行段,这段占整个飞行轨迹最大,耗时最久。在第 1.2节的基础上,对上升段积分可得到上升段终点参数kkk,vxy,若已知自由漂浮段的终端高度h,可以计算自由漂浮段的射程和终端参数。a)自由漂浮段射程计算。导 弹 与 航 天 运 载 技 术(中英文)2023年 54 在假设月球为匀质圆球的条件下,飞行器自由漂浮段轨迹应在上升段终点的飞行参数决定的轨道平面内,该平面过月心。自由漂浮段的角射程c为在发射平面内,自由漂浮段起点(上升段终点)到自由漂浮段终点所对应的地心角3,角射程c可根据式(3)计算得
16、到。2c4tan22+-=BBACA (3)2ekkekk ekkek22kkke2kkkkkkkkk2(1tan)()2tan()()/arctan=+-+=-=+=+=+=+ArrrBrCrrrxyRrRhvrxRy (4)式中 er为自由漂浮段终端与月心的距离;kr为上升段终端与月心的距离;k为能量参数,表示轨道上一点的动能的两倍与势能之比;k为上升段终端与当地水平面的倾角;k为上升段的射程角;k为上升段终端在发射坐标系下的速度倾角。因此,自由漂浮段的射程为 kcc=LR (5)b)自由漂浮段段终端参数计算。在发射坐标系下,自由漂浮段段终端参数为 eekceekc2kekeeekccek
17、kesin()cos()22()1()(sinsin)=+=+-=-+=-+=-+-xryrRvvrrTEEeEEn (6)式中 k为上升段的射程角;eeee,x y v分别为在发射坐标系下自由漂浮段终端位置、速度和速度倾角;cT为自由漂浮段飞行时间。其中的轨道参数为 k2kk2kkk3222kkke22ee21(2)coscosarccos1cos(1)=-=+-=-=-rar venar vr vrEea (7)式中 a为轨道半长轴;e为轨道偏心率;n为飞行器在椭圆轨道上飞行的平均角速度;e自由漂浮段终端与当地水平面的倾角;E为偏近点角。1.4 第 1 次动力减速段 第 1 次动力减速段是
18、在自由漂浮段之后,开始动力减速,主要作用是在降低飞行器的飞行速度的同时对准目标点,使其终端的速度,速度倾角满足实现垂直动力软着陆的要求。第 1 次动力减速段在发射坐标系下的运动学方程与式(1)一样。但推力在发动机关机后为零,即:eddp(),()0,=P ttttP tttt (8)式中 et为第 1 次动力减速段开始时间;dt为发动机关机时间;pt动力减速段到达指定高度的结束时间。对于第 1 次动力减速段,同样需要设计俯仰程序角,可将第 1 次动力减速段分为固定程序角减速段、程序转弯段、垂直下降段,如图 3 所示。设计俯仰程序角表达式为 1222nn21mmnp,()(),(90)/290,
19、tttttttt ttttt=-=+-+-?(9)式中 1为初始固定俯仰角;?为俯仰角速率,可以设计为常值;m为目标射程角,90-m为发射点坐标系下垂直于目标点水平面的俯仰角。对于第 1 次动力减速段,需要优化的控制变量包括初始高度h,固定俯仰角1,固定程序角减速段工作时间12edttt=-,垂直动力减速时间2dndttt=-。1.5 垂直动力着陆段 前 3 个飞行阶段设计采用跟踪标准轨迹的制导方法(即标称制导法),本文主要对月面机动任务标称轨迹进行优化,而最后垂直动力着陆段要求将飞行器精确地降落到目标点且终端速度接近于零。为了实现垂直动力着陆段在线精确制导,本文采用基于凸优化的在线制导方式,
20、切换高度设置在 2000 m,主要研究基于凸优化的在线轨迹优化方法。考虑在 2000 m 高度范围内,重力加速度mg变化在 0.2%左右,所以在垂直动力着陆段可以忽略月球引力的变化。另外,虽然前 3 个阶段标称轨迹的设计中未考虑侧向运动,但是在跟踪制导过程中会产生较大的终端位置与速度误差,垂直动力着陆段需要修正这些制导误差,因此,需要考虑飞行器的侧向运动,所以在目标点坐标系TTTTO X Y Z下,考虑三自由度的垂直动力着陆段的动力学方程可以表示为 任 宽等 月面大范围空间机动初步优化 55第1期 mspe()()()()/()()()ttttm ttm tI grvvgPP=+=-?(10)
21、式中 ()tr,()tv分别为飞行器在目标点坐标系下的位置和速度矢量。2 月面机动任务轨迹优化 月面机动任务轨迹优化主要针对从上升段到第 1 次动力减速段结束这部分轨迹,分析月面机动任务,上升段终端状态(位置和速度)将决定整个机动轨迹的射程,另外,动力减速的开始高度、减速方式和工作时间对于整个机动任务燃料消耗也有显著影响。根据1.5节的月面机动动力学方程,选择上升段关机时间2t,上升段终端俯仰角0,第 1 次动力减速段开机高度h,初始固定俯仰角1,固定程序角减速时间1dt,垂直动力减速时间2dt作为优化变量。由此可得:a)优化变量:20112,d,d thtt=x (11)b)目标函数:min
22、+上升段第1次动力减速段Jtt=x (12)设置发动机固定推力工作,可以以发动机工作时间作为优化变量,上升段t为上升段发动机工作时间,第1次动力减速段t为第 1 次动力减速段发动机工作时间。c)终端约束:终端射程角约束:pm=(13)终端速度倾角约束:pm905+(14)终端速度约束:pmaxvv (15)终端高度约束:22ppm()xyRRh+-=(16)式中 p为第 1 次动力减速段终端的射程角;m为目标射程角;p为在发射坐标系下第 1 次动力减速段终端速度倾角;pp,xy,在发射坐标系下的终端位置。采用经典的四阶龙格库塔法对运动学方程积分以得到整个机动轨迹的终端参数,采用遗传算法对整个优
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 月面大 范围 空间 机动 初步 优化 任宽
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。