圆及椭圆轨道释放绳系子星的方法研究

研究绳系卫星系统质心沿椭圆轨道运动时，子星释放的最佳初始条件，运用变量替代简化动力学控制方程；对圆轨道运动的母星，采用轨道转移方法释放子星，计算弹射变轨参数及过渡轨道终了真近点角。本文取消母星质量远大于子星质量的假设，考虑子星释放过程中绳系卫星系统质心随子星运动的变化，所得结果更具有普遍意义。数值模拟结果表明，所述方法的有效性。     

绳系卫星横向振动的控制方法

研究绳系卫星绳索横向振动问题。首先建立绳索振动方程,由Pontrgasin方法求得抑制绳索横向振动的纵向优化控制力。实现方式为开启母星体内的电机控制绳长使星体产生加速度,星体惯性力导致绳轴向附加张力。      

绳系卫星部署阶段的动力学分析

绳系卫星系统部署阶段末时刻状态决定着其编队飞行的初始状态，为了研究不同因素对绳系卫星系统部署阶段运动的影响，本文建立了一种简单的平面哑铃模型，将系绳视作有阻尼的弹簧，两个卫星视作刚体，考虑重力梯度力矩。通过仿真，发现系绳释放的速度和副星推力是影响绳系卫星系统部署阶段稳定性的主要因素。     

二体绳系卫星系统冲击效应动力学与试验研究

为了研究冲击效应对二体绳系卫星系统动力学特性的影响，开展了基于绝对节点坐标法的柔性绳索模型及绳系卫星系统运动特性分析研究。首先介绍了绝对节点坐标法柔索模型的建立过程，并进行了动力学仿真实验，在仿真时充分考虑了柔索的轴向刚度矩阵、弯曲刚度矩阵以及质量矩阵，可在不增加计算量的同时保留系绳真实特性。最后，搭建了三自由度的二体绳系卫星地面模拟试验系统，在规格为20 m×30 m的大型气浮平台上开展了试验验证，对末尾释放阶段和初始回收阶段系绳的张力以及子星的运动状态进行了分析。验证结果表明：地面模拟试验与仿真实验相比，系绳张力误差均在5%以内，验证了该建模方法的有效性和所搭建的试验系统的可靠性，同时为绳系卫星的未来发展和地面试验的进一步展开奠定了基础。     

绳系系统轨道机动过程中的面内摆角抑制

摘要： 针对绳系卫星系统轨道机动过程中面内摆角抑制问题进行控制器设计.仅考虑面内摆角子系统下,针对绳系系统的动力学模型,提出了分层迭代非线性滑模控制方法.证明控制方法设计的控制律使系统所有状态渐近稳定.通过数学仿真验证方法的有效性和优势.     

基于非线性模型预测的绳系系统系绳摆振控制

针对绳系系统离轨稳定控制问题,开展了系绳和绳端卫星构成的绳系系统在拖拽离轨过程中系绳摆动稳定控制方法研究。在考虑拖拽离轨过程中约束下,首先建立了包含绳端卫星的姿态运动的模型,并建立了绳系系统的离轨动力学方程和便于控制器设计的简化动力学方程。其次根据模型预测原理以最优化方法设计参考轨迹,最后以模型预测控制方法为基础设计了稳定系绳摆动的非线性模型预测控制方法。使用MATLAB软件平台仿真,验证了所设计的参考轨迹能完成目标和模型预测控制器有好的跟踪能力。      

绳系卫星系统耦合振动频率的确定

将N(N≥2)体绳系卫星系统的各节系绳分为足够多的小段,如共有M-1段,每小段上的系绳张力可认为是常量。这样由N 体问题转化来的M 点问题即可逼真系绳张力沿系绳的变化。利用中间连接点的边界条件,M 点边界值问题可简化为等效的两点边界值问题,从而平面内耦合振动频率的求解归化为四阶行列式的求值。由绳系卫星系统运动中心概念推导出的每小段系绳张力的平均值,被用于系统振动频率的求解。并给出一套适合于编制计算程序的求解轨道平面内耦合振动频率的数学表达式,以及一组模拟计算结果。     

空间微重力组合体及其控制

空间微重力组合体是绳系卫星系统的升级版,有更强大的功能和应用前景。可是在无控状态下,新旧两款系统的留位（系绳伸出长度不变）和回收（系绳回卷收缩）运动都不稳定,必须置入有效的控制算法才能使其稳定运行。从绳系卫星的简化动力学模型入手,抓住微重力和动态运动中的惯性力相互关系,推导出基础型和扩展型的距离速率控制算法,简明有效地抑制这些系统的不稳定运动和改进空间组合体在轨交会对接控制技术。该算法也可应用于有分布质量和弹性变形的连续系绳的动力学模型以及其他更复杂的绳系卫星（家族）系统的控制中。     

基于动力学递推算法的绳系卫星系统刚柔耦合多体模型

利用动力学递推算法建立了空间绳系卫星系统(TSS)的铰接柔性杆动力学模型.所研究的绳系系统包含两颗刚体绳端卫星和一段连接两颗卫星的柔绳,柔绳离散为一系列球铰连接的弹性杆,考虑杆不均匀的纵向一维形变,并利用球铰的运动模拟柔绳的弯曲和扭转;然后基于递推算法推导得出了该刚柔混合离散模型的动力学方程.数值仿真结果表明:所建立的绳系卫星系统模型能够很好地模拟空间系绳的各向摆振和扭转,此外,递推算法的引入能够有效地减少离散模型动力学方程的维数,从而大幅减少计算量.      

空间绳系拖拽系统摆动特性与平稳控制

考虑了任务星与废星的姿态运动以及系统组合体的面内外姿态运动,建立了绳系拖拽离轨系统动力学与控制模型,以切向常值推力下绳系拖拽轨道转移为任务过程,分析了任务星在喷气和零动量轮的限制姿态反馈控制条件下飞行时,废星姿态摆动、系统组合体面内外摆动和任务星姿态运动的规律及相互影响关系。采用留位和阻尼控制相结合的系绳张力复合控制方法,并结合任务星姿态控制,确保绳系拖拽转移安全平稳进行。仿真结果表明:常值推力下绳系拖拽轨道转移时,牵挂点偏置诱发的废星姿态周期性摆动会激发绳系组合体的面内外同频率高阶摆动,星体姿态运动是任务星姿态扰动力矩产生的主要因素;采用张力复合控制可有效消除废星姿态摆动并保持星间相对距离,结合任务星姿态控制,可实现离轨过程的平稳与安全,大幅减少任务星的姿控能耗。      

一种基于火箭视位置与视速度的航天器初轨确定方法

为了充分利用航天器初始轨道确定的信息源,不断完善地面测控系统定轨手段和方法,提出了一种基于火箭视位置与视速度的航天器初始轨道确定方法.首先对火箭视位置与视速度测量弹道的原理进行了分析;然后利用火箭秒节点计算机字(视速度)以及某一点的初始弹道(初始位置矢量、速度矢量),通过对火箭的运动方程进行逐秒积分,以获得每秒的位置矢量和速度矢量,再转换为弹道数据,并利用此数据直接确定航天器初始轨道.最后通过某太阳同步轨道卫星实测数据仿真计算,证明了方法的正确性和可靠性.      

轨道机动过程中推力加速度的在线最小方差估计

定位和跟踪空间机动目标时，对目标运动建模受发动机推力的不确定性影响，通过统计处理离散的雷达观测数据实时估计发动机的推力，进而定位和跟踪机动目标便是本文所要研究和解决的问题，本文在地心惯性系建立了常推力轨道机动过程中连续变质量运动模型和离散雷达量测模型，机动过程中质量秒耗量和排气速度作为表征轨控发动机推力的两个近似常量，应用扩展卡尔曼滤波对离散雷达测量数据进行序贯统计处理得到发动机推力的最小方差估计；文中详细地给出了线性化量测模型的变分方程和观测矩阵；仿真结果表明该算法能快速、准确地在线估计轨控发动机的等效推力。     

一种基于星载GPS的月球探测器导航算法

为验证GNSS旁瓣信号对中高轨航天器在轨飞行的导航支持能力,"嫦娥5号"飞行试验器搭载了GPS导航接收机。由于目前使用的星载导航解平稳性不够,为此,根据月球探测器不同飞行段落轨道动力学特性,建立了一种基于星载GPS数据的月球探测器实时导航算法。该算法充分考虑了月球探测器不同段落的动力学特性,基于UKF滤波设计,并引入了自适应策略,以降低滤波参数选取的难度。在分析星载GPS数据特点基础上,月地返回段实测数据的处理结果验证了算法的可行性和有效性。     

同步三星目标运动状态快速检测方法

针对星载无源定位系统难以从定位结果及时准确地检测目标运动状态的问题,利用同步三星时差定位体制的定位误差特性,根据静止假设下同步双星时/频差与三星定位结果应相匹配的原理,提出了基于定位一致性的目标运动状态快速检测算法.先通过定位结果的误差分析掌握目标实际位置所在区域,然后在静止假设下统计该区域的频差范围,并将定位时刻的频差观测值与之比较,实现目标运动状态检测,并且推导了检测概率的理论表达式.仿真和试验结果验证了算法的有效性,时/频差估计误差标准差为10 μs/100 mHz时,对速度大于40 m/s的目标,检测概率高于99%.实测数据表明该算法显著提高了系统目标运动状态检测的准确性和时效性.     

On the possibility of detecting some terms of the diurnal polar motion by the study of satellite orbits

The high precision of satellite laser ranging data obtained by the laser stations and the improvement of the models of forces acting on the satellite permit us to estimate the pole position with a few milliarsecond precision. In addition to the well-known precession and nutation in a space-fixed system the earth's rotation axis is affected by a diurnal motion which is superimposed to the polar motion in the terrestrial system. However, the way the pole position is usually estimated by determining a constant value and assuming a linear drift over a few days, the diurnal polar motion is averaged to zero.After some study of the geometry of the system and the relation between the diurnal polar motion and the satellite orbit, we have performed a series of simulations to see if this variation can be separated from variations in the orbit. The main result is that a specific quasi-diurnal term (corresponding to the fortnightly term in space) can be estimated with a milliarcsecond accuracy. Due to the effect of the tides at the same frequency, no significant value has been determined using real data.     

基于IMM-UKF的航天器相对位姿确定算法

在航天器相对导航过程中,相对距离测量信息容易受到干扰,测量误差有较大的不确定性,通常基于单一模型的滤波算法无法对噪声进行辨识,很难获得精确的导航结果。针对应用Clohessy-Wiltshire(C-W)方程受到圆轨道假设的限制问题,研究了建立在惯性坐标系下的近距离相对运动方程(Lawden方程),建立了基于这两个方程的模型集。根据导航系统测量敏感器的特点,设计基于Rodrigues参数及无迹卡尔曼滤波(UKF)的交互式多模型(IMM)视觉相对位姿动态估计算法(IMM-UKF),在保证计算效率的前提下,确保相对轨道姿态确定的稳定性和精确性。数值仿真验证了算法的有效性和先进性。     

基于视线测量和轨道预报高轨非合作目标相对导航方法

对于非合作目标,由于中远距离星上相对测量手段有限,大多情况仅能获得视线角信息.仅视线测量相对导航方法在GEO轨道条件下滤波精度低、可观测性差.提出一种基于星间视线方位测量和轨道预报信息结合的非合作目标相对导航方法.建立基于星间相对运动模型的状态方程和基于星间视线测量和轨道预报信息的观测方程,分别选取了扩展卡尔曼滤波和无迹卡尔曼滤波两种方法,仿真分析了轨道预报信息精度和滤波方法对导航精度的影响.     

月地转移轨道快速设计方法

月地转移轨道设计一般分为初步轨道设计和精确轨道设计.其中,初步轨道设计的准确性是确保后续精确轨道设计收敛的关键.提出了一种基于Lambert算法的月地转移轨道快速设计方法.以出月球影响球的时刻、位置和速度为中间变量,将轨道分为地心段和月心段分别进行计算.将探测器飞出月球影响球至指定再入点的地心段轨道简化为一个Lambert问题进行求解,提出了通过牛顿迭代法求解月地转移轨道Lambert问题的方法,避免了Lambert问题求解时大量的超几何函数和级数计算,提高了计算效率.在月心段轨道的快速计算中,提出了根据探测器出影响球速度矢量、月球停泊轨道倾角和近月点高度计算月心双曲线轨道根数的新方法.通过迭代计算,使得两段轨道在月球影响球处的位置和速度连续,从而获得一条完整的满足两端约束的双二体月地转移轨道.该方法计算速度快,精度相对较高.计算结果可以作为后续精确轨道设计的初值.      

Attitude stability and periodic attitudes of rigid spacecrafts on the stationary orbits around asteroid 216 Kleopatra

In this work, equilibrium attitude configurations, attitude stability and periodic attitude families are investigated for rigid spacecrafts moving on stationary orbits around asteroid 216 Kleopatra. The polyhedral approach is adopted to formulate the equations of rotational motion. In this dynamical model, six equilibrium attitude configurations with non-zero Euler angles are identified for a spacecraft moving on each stationary orbit. Then the linearized equations of attitude motion at equilibrium attitudes are derived. Based on the linear system, the necessary conditions of stability of equilibrium attitudes are provided, and stability domains on the spacecraft’s characteristic plane are obtained. It is found that the stability domains are distributed in the first and third quadrants of the characteristic plane and the stability domain in the third quadrant is separated into two regions by an unstable belt. Subsequently, we present the linear solution around a stable equilibrium attitude point, indicating that there are three types of elemental periodic attitudes. By means of numerical approaches, three fundamental families of periodic solutions are determined in the full attitude model.