XNAV/UVNAV/SINS组合导航在航天器轨道机动中的应用

针对X射线脉冲星导航在航天器轨道机动过程中精度不高甚至发散的问题,提出一种将X射线脉冲星导航结合惯性导航和紫外敏感器的组合导航方法。以航天器在惯性系中的位置、速度、姿态四元数和惯性导航设备误差作为系统状态变量,用X射线探测器测量X射线脉冲到达时间,用紫外敏感器测量中心天体质心相对于航天器的方向矢量和距离以及航天器在惯性系中的姿态四元数,用扩展卡尔曼滤波器估计组合导航系统状态。仿真结果验证了该组合导航方法的可行性,能够解决轨道机动中X射线脉冲星单独导航的误差过大(位置误差达107m)问题,且该组合导航具有较高的导航精度,在轨道机动前、机动中和机动后导航位置误差均在100 m以内。     

基于脉冲星和小行星的组合导航在深空巡航段的应用

为提高航天器在深空探测巡航段的自主导航能力,提出一种利用脉冲星和小行星信息的自主导航系统.利用单X射线探测器测量脉冲星辐射的光子到达时间,利用导航相机获取小行星图像,经信号处理,获取脉冲星的脉冲到达时间和小行星中心点对应的像元像线,利用集中式扩展卡尔曼滤波算法进行状态估计,并修正星载时钟钟差.以美国“深空一号”任务深空巡航段轨道为例进行了仿真实验,仿真结果表明了该方法的有效性,能够完成深空探测器在巡航段的自主导航任务,同时,组合导航解决了基于小行星光学信息单独导航时的星载时钟钟差发散问题,且同光学导航相比,位置估计精度提高了57.18％.     

基于单探测器的X射线脉冲星深空导航算法

提出一种基于单探测器的X射线脉冲星深空导航算法.该算法利用单个X射线探测器在不同时刻的观测数据来构建观测模型,因此航天器不必搭载多个X射线探测器,大幅降低了X射线脉冲星深空导航系统的质量和功耗.针对不同的观测方案,利用线性化可观性分析方法对导航系统进行了可观性分析.分析结果表明,航天器应当分时段观测不同的脉冲星来提高导航性能.最后,给出了相应的仿真算例.仿真结果表明利用该算法能获得与传统算法相近的导航精度.      

基于点云矩形面特征的故障航天器位姿测量

针对非合作航天器的相对导航问题,提出了一种利用点云矩形面特征测量非合作航天器位姿的方法。首先从点云数据中提取矩形面;然后根据矩形面点云数据计算出点云分布矩阵,通过特征值分解求出相对位置和姿态,并解决了因矩形面对称而产生的多解问题;最后设计了卡尔曼滤波器,确定目标星与追踪星的相对位置参数以及目标星的姿态、角速度。仿真结果表明:相对位置的估计精度优于0.005m,目标星姿态精度优于0.1°,验证了该方法的有效性。     

基于脉冲星导航系统的X射线探测器

X射线脉冲星导航中脉冲到达时间（TOA）是脉冲星导航系统中基本的测量量,其测量精度直接影响导航定位的精度。国内外对用于X射线脉冲星导航的探测器研究还处于探索阶段。通过比较不同类型探测器的特性,根据导航定位精度对探测器各项性能指标的要求,提出一种X射线探测器系统设计方案,为用于脉冲星导航的X射线探测器的研制提供一种参考。本方案用CCD作为探测器主体,整个探测器由光学结构、CCD及电子学系统组成。     

惯性/卫星/里程计多信息融合方法及在铁路测绘中的应用

针对基于惯性技术对铁路基础设施进行精确测绘的需求提出一种多信息融合惯性基准方案,为测量测绘提供高精度位置和姿态参考。对载体运动特点和车载状态下惯性/里程组合导航航向角误差可观性进行分析,认为天向陀螺漂移和航向误差是导致测量精度下降的主要因素,针对该问题设计了基于双向滤波、双向平滑的多信息融合方案,针对缺乏绝对位置基准的应用情况,引入"正矢"概念和相对定位精度的评判方法。仿真及试验结果表明,在陀螺常值漂移0.2(°)/h条件下,该方案相对定位精度优于0.3 mm(300 m弦正矢),显著提高了车载铁路线路测绘位置、姿态基准精度,降低了对惯性器件的要求,利用中、低精度器件实现了高精度测量定位。     

一种基于加速度计和天文的卫星自主导航新方法

直接敏感地平是一种典型自主天文导航方法,该方法简单可靠,易于实现,但是由于常用卫星轨道动力学J2模型精度有限,地球敏感器精度较低,因此导航精度不高。加速度计是测量运载体线加速度的常用惯性导航设备,当航天器在轨运行时,星载加速度计能够测量航天器所受发散力。结合上述两种方法的特点,提出一种将加速度计和天文相结合的自主天文导航新方法。在常用卫星轨道动力学模型基础上,引入大气阻力和太阳光压系数模型作为自主导航系统状态方程的一部分,并建立近地空间环境下星载加速度计的测量模型,将其与直接敏感地平均作为导航系统观测方程。设计基于信息融合的自主导航滤波方法,通过对多种导航模式进行数值仿真及结果分析,结果表明所设计方法提高了系统定位精度62.8%和速度精度63.9%,增强了系统可靠性。     

柱面坐标系下航天器仅测角相对导航算法

针对采用星载无源光学相机实现对空间非合作目标的相对导航时轨道状态不可观测/弱可观测的问题,研究了基于柱面坐标系相对运动动力学模型的仅测角相对导航算法。首先,在柱面坐标系下建立了三维的相对轨道运动动力学模型和视线角测量模型;然后,基于该模型进行了仅测角相对导航的状态可观测性分析,定性地获得了相对轨道状态可观测性更佳的几何条件;接着,为了解决强非线性模型的估计问题,建立了基于平方根容积卡尔曼滤波的仅测角相对导航估计算法及长航时导航滤波方案;最后,对所提算法进行了数值仿真实验验证。仿真结果表明柱面系下动力学模型的轨道曲率捕获能力强,能够为仅测角相对导航提供足够可观测性,相对轨道估计误差可收敛至公里级,能够满足中远程空间态势感知、交会的精度要求。     

一种航天飞行器的INS/CNS自主导航方案

传统的天文导航方法由于受水平基准精度的制约,难以进一步提高定位精度。针对INS/CNS组合导航方法对水平基准依赖的问题,提出了一种新颖的基于星光折射技术的INS/CNS自主组合导航方案。将惯性系下的非线性惯导误差传播方程作为系统状态方程,将星敏感器测得姿态和星光折射信息作为量测,采用UKF滤波算法,构成全面最优的INS/CNS组合导航。仿真结果表明,星敏感器精度为3″时,导航系统的定位精度优于200 m,姿态精度优于3″,导航定位精度随所使用的折射星数目增多明显提高,且方案在系统大角度误差条件下仍然适用。     

基于模糊迭代均方根容积卡尔曼滤波的天基非合作目标跟踪

针对非合作航天器相对导航中测量噪声不确定的问题,提出了一种模糊迭代均方根容积卡尔曼滤波算法,实现对非合作目标相对状态的测量。该算法利用容积点均方根迭代策略和模糊推理系统实时调整改进容积卡尔曼滤波的量测噪声协方差阵权值,修正量测噪声协方差阵,使其接近真实噪声值,从而提高目标跟踪算法的自适应能力,提高了滤波精度。通过建立数学仿真模型,分别采用扩展卡尔曼滤波、容积卡尔曼滤波以及模糊迭代均方根容积卡尔曼滤波进行跟踪仿真,仿真结果表明,与标准容积卡尔曼滤波相比,该改进算法能够提高13.17%的跟踪精度。     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

Effet de l'anisotropie élastique cristalline sur la distribution des facteurs de Schmid à la surface des polycristauxInfluence of crystalline elasticity anisotropy on Schmid factor distribution at the free surface of polycrystals

Experimental studies of the plasticity mechanisms of polycrystals are usually based on the Schmid factor distribution supposing crystalline elasticity isotropy. A numerical evaluation of the effect of crystalline elasticity anisotropy on the apparent Schmid factor distribution at the free surface of polycrystals is presented. Cubic elasticity is considered. Order II stresses (averaged on all grains with the same crystallographic orientation) as well as variations between averages computed on grains with the same crystallographic orientation but with different neighbour grains are computed. The Finite Element Method is used. Commonly studied metals presenting an increasing anisotropy degree are considered (aluminium, nickel, austenite, copper). Concerning order II stresses in strongly anisotropic metals, the apparent Schmid factor distribution is drifted towards small Schmid factor values (the maximum Schmid factor is equal to 0.43 instead of 0.5) and the slip activation order between characteristic orientations of the crystallographic standard triangle is modified. The computed square deviations of the stresses averaged on grains with the same crystallographic orientation but with different neighbour grains are a bit higher than the second order ones (inter-orientation scatter). Our numerical evaluations agree quantitatively with several observations and measures of the literature concerning stress and strain distribution in copper and austenite polycrystals submitted to low amplitude loadings. Hopefully, the given apparent Schmid factor distributions could help to better understand the observations of the plasticity mechanisms taking place at the free surface of polycrystals. To cite this article: M. Sauzay, C. R. Mecanique 334 (2006).     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).     

Dynamics and Breakup of Viscoelastic Liquids (A Review)

The phenomena of hydrodynamic breakup of liquid jets, drops, films, bridges, and filaments are reviewed for liquids with viscoelastic properties. The reasons for breakup are capillary instabilities, collisions with rigid obstacles, and other forms of dynamic action. The relationship between the properties of the liquids and the features of the breakup process is discussed.     

Flow stability analysis and excitation using pulsating jets

Classical flow stability applied to transition from laminar to turbulent flow may also describe the behavior of vorticity fluctuations created by a pulsating jet placed along a solid boundary. A numerical laminar flow experiment involving a pulsating jet placed along the surface of a duct with flow separation downstream, resulted in eliminating most part of the separated flow region. Applying the same approach to a turbulent flow, it was possible to develop a turbulent stability flow formulation and apply successfully turbulent pulsating jet flow separation control. To cite this article: D. Skamnakis, K. Papailiou, C. R. Mecanique 333 (2005).