[发明专利]一种小天体探测器着陆动力学模拟系统有效
申请号: | 201510977475.0 | 申请日: | 2015-12-23 |
公开(公告)号: | CN105631099B | 公开(公告)日: | 2019-06-28 |
发明(设计)人: | 阮晓钢;肖尧;张晓平;黄静;朱晓庆;陈志刚;林佳;柴洁;陈岩;伊朝阳;刘冰;李诚 | 申请(专利权)人: | 北京工业大学 |
主分类号: | G06F17/50 | 分类号: | G06F17/50 |
代理公司: | 北京思海天达知识产权代理有限公司 11203 | 代理人: | 沈波 |
地址: | 100124 *** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明公开了一种小天体探测器着陆动力学模拟系统,该系统包括探测器参数编辑器、小天体重力计算引擎、动力学模拟引擎、API绑定层;小天体重力计算引擎根据目标天体的三位多面体模型或者球谐系数参数,完成探测器当前位置重力加速度的计算;动力学计算引擎根据当前的参数配置生成探测器的3自由度/6自由度动力学模型,并根据API接口绑定层传入的发动机控制量、星历数据库、天体参数数据库相关数据和小天体重力计算引擎,完成探测器的运动和姿态计算,并通过API接口绑定层返回给客户端程序,最终完成了小天体探测器着陆动力学模拟计算,使客户端可以通过不同的语言接口对模拟系统进行调用,传输相关数据。 | ||
搜索关键词: | 一种 天体 探测器 着陆 动力学 模拟 系统 | ||
【主权项】:
1.一种小天体探测器着陆动力学模拟系统,其特征在于:该模拟系统主要包括:探测器参数编辑器、小天体重力计算引擎、动力学计算引擎、API绑定层;探测器参数编辑器快速完成系统各参数配置,包括探测器机体坐标系的确定、转动惯性矩阵的计算,发动机的参数配置、天体固连坐标系的确定、着陆坐标系的确定;小天体重力计算引擎根据目标天体的三维多面体模型或者球谐系数参数,完成探测器当前位置重力加速度的计算;动力学引擎根据当前的参数配置自动生成探测器的3自由度/6自由度动力学模型,并根据API接口绑定层传入的发动机控制量、星历数据库、天体参数数据库相关数据和小天体重力计算引擎,完成探测器的运动和姿态计算,并通过API接口绑定层返回给客户端程序,最终完成了小天体探测器着陆动力学模拟计算;系统包括探测器参数编辑器、小天体重力计算引擎、动力学计算引擎、API绑定层;1)参数编辑器的基本原理是借助于计算机图形学相关技术,通过小天体三维多面体模型或者探测器的三维模型快速确定探测器的各项参数和相关坐标系,其实施方案如下,(1)小天体的质心、最小和最大转动惯量方向的计算方法:首先,对导入小天体三维多面体模型进行处理,将所有不是三角形的面全部再次分割成三角形,每个三角形面和坐标原点组成了一个四棱柱;其次,根据天体平均密度,计算每个四棱柱的体积、质量、质心位置;第三,计算天体质心其中,n为四面体的总数,Mi和ri分别表示第i个四面体的质量和质心位置;(2)根据探测器的模型、各部件材质对应的密度参数,使用(1)同样的方法计算探测器的质心、最大转动惯量方向;(3)借助(1)、(2)的计算结果,在小天体的三维模型上确定天体固连坐标系和着陆坐标系原点的位置,在探测器的三维模型上确定机体坐标系,进而计算出探测器在机体坐标系下的转动惯性矩阵;(4)借助探测器的三维模型和发动机在模型中的位置,确定各发动机推力方向、作用点参数,推力大小、比冲和燃耗比需要根据实际情况手动输入确定;2)小天体重力计算引擎根据目标天体的三维多面体模型或者球谐系数参数,完成探测器当前位置重力加速度的计算;多面体重力模型如下:其中,r为计算点的位置向量,G为重力常数,ρ为小天体密度,edge和face为三维多面体模型所有棱边和三角形平面的集合,棱边e为三角形平面A和三角形平面B的共同边,其两端分别为顶点1和顶点2;re为计算点到棱边e上任意一点的向量;为平面A的外法线方向单位向量,平面A内的边e的外法线方向单位向量,的平面B相关的法向量,其定义与平面A对应的定义一致,Ee为3×3矩阵;其中r1,r2分别为计算点到棱边e两个端点的距离,e12为棱边的长度;rf为计算点到平面f任意一点的向量,为平面f的单位法向量,Ff为3×3矩阵;α=rirjrk+ri(rj·rk)+rj(rk·ri)+rk(ri·rj),β=ri·(rj×rk),ri、rj、rk分别为计算点到三角形平面f的3个顶点的向量,ri、rj、rk分别为对应的向量长度;引力势能U对r进行一阶求导,得到重力加速度,其表达式如下:多面体重力模型能够计算任意形状,但密度均匀的小天体任意点的重力加速度;球谐系数重力模型其表达式如下:其中,r为计算点距离目标天体质心的距离,λ为经度,φ为纬度,GM为目标天体重力常数,R为目标天体参考半径,为规范化系数,其具体数值由重力计算引擎从小天体球谐函数稀疏参数数据库中读取和确定;另r=[x y z]T为计算点在天体固连坐标系中的笛卡尔坐标表示,计算点的重力加速度为引力势能U对r的一阶导数:其中,分别为引力势能对r、φ、λ的一阶偏导数;其中,球谐系数重力模型计算速度快,计算结果精确,但是当计算点在最小半径圆之内,即r<R时计算结果会发散;由于小天体形状极其不规则,天体地表附近的部分位置处于最小半径圆之内,在探测器着陆过程中,当探测器位置在球谐系数重力模型最小半径圆之外,系统使用球谐系数重力模型计算重力加速度,以提高计算速度,当探测器穿过最小半径圆后,系统强制使用多面体重力模型计算重力加速度;3)动力学计算引擎具体实现方案如下:根据当前的参数配置生成探测器的3自由度/6自由度动力学模型;探测器质心运动方程为:其中,m为探测器的质量,t为时间,P为探测器推力矢量,G为当前位置目标天体的重力加速度与各大天体对探测器的引力摄动之和,目标天体的重力加速度由重力计算引擎根据当前位置实时计算得到,各大天体对探测器的引力摄动则根据探测器当前位置和时间、星历数据库、天体参数数据库,计算太阳系各大天体相对于探测器的位置,进而根据万有引力定律计算得到;ωe为目标天体的自转角速度,r为探测器质心在惯性坐标系中的位置;探测器绕质心转动的动力学模型为:其中,I为探测器的转动惯量矩阵,ωT为探测器旋转角速度,t为时间,MC为控制力矩矢量,Mk为附加哥氏力力矩矢量;上述两个方程构成了探测器的6自由度动力学模型,当设定为对探测器进行运动学仿真计算时,系统只生成质心运动方程,探测器的姿态由客户端经API接口绑定层传入;根据生成探测器的3自由度/6自由度动力学模型、小天体重力引擎计算和各大天体引力摄动计算结果,并根据API接口绑定层传入的发动机控制量,借助Runge‑Kutta算法,对探测器动力学模型进行积分,完成探测器的运动和姿态计算,并通过API接口绑定层返回给客户端程序,最终完成了小天体探测器着陆动力学模拟计算;4)API接口绑定层完成C/C++、Python、Matlab、Fortran、Java五种语言的接口绑定,使客户端能够通过不同的语言接口对模拟系统进行调用,传输相关数据,其调用顺序如下:(1)初始化,确定仿真初始时间、仿真时间间隔、探测器初始位置、姿态、速度、仿真结束条件;(2)设定控制量输入方式,通过回掉函数进行计算,或者每个时间点,由客户端计算并调用接口传入控制量;(3)启动动力学计算引擎开始进行仿真计算,计算结果由函数返回值进行返回;(4)到达仿真时间长度或者仿真结束条件,计算结束并退出。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京工业大学,未经北京工业大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201510977475.0/,转载请声明来源钻瓜专利网。