一种大角加速度下的星跟踪算法

在飞行器大角加速度机动情况下,星敏感器星跟踪算法预测坐标的偏差变大导致星点扫描范围变大,星点识别时间变长,数据更新率降低,为此提出了一种采用联合星敏感器和陀螺信息预测星像坐标的星跟踪算法.该方法首先由组合导航输出的上一时刻的姿态数据和一个周期内修正后的陀螺数据积分推算下一时刻的飞行器姿态,然后通过坐标转换计算出已知恒星下一时刻在星敏感器坐标系上的星图坐标.理论分析和仿真结果表明,该方法飞行器最大角加速度30(°)/s2,最大角速度30(°)/s的情况下,预测星敏感器视场内恒星的星像坐标的偏差在10个像素以内,可使星图扫描和星点匹配快速完成,提高星敏感器的工作效率.     

卫星惯性/星光组合自主定姿方法研究

利用星敏感器测量所得的恒星星光方向信息,采用六状态变量推广卡尔曼滤波,得到卫星姿态误差和陀螺漂移误差信息并进行相应的修正。仿真结果表明,这种方法修正效果良好,采用中等精度陀螺组件就可以实现高精度定姿;仿真结果还验证了陀螺随机误差、星敏感器测量噪声和滤波周期等因素对定姿精度的影响。     

一种适用于星敏感器的大角速度下星跟踪算法

传统的星敏感器预报方法是利用当前四元数和陀螺输出角速度来预测下帧姿态四元数,再通过预测的四元数计算视场内恒星的理想星像坐标,最后在这些理想坐标为中心的范围内来获取实际的星像坐标。但是,由于陀螺漂移的存在,随着时间的增加,预测的理想坐标与实际坐标偏差越来越大,获取实际星像坐标扫描的区域越来越大,星敏感器的数据更新率也随之降低。为此,本文提出了一种适用于星敏感器的大角速度下星跟踪算法。首先根据连续两帧的姿态偏差计算星敏感器的角速度,利用角速度来计算星敏感器当前时刻与最近上一帧的姿态偏差。然后以该姿态偏差来预报当前时刻的星敏感器姿态输出,这样不仅提高了数据更新率,而且满足飞行器在大角速度机动时单纯采用星敏感器进行姿态测控的要求。最后,采用2010年外场试验实际观星数据对该算法进行了仿真验证。     

多探头星敏感器星图融合姿态确定精度

针对多探头星敏感器星图融合姿态确定精度问题做了研究。首先提出多探头星敏感器星图融合方法;其次在此基础上以三个探头为例,推导了三探头星敏感器星图融合姿态确定误差协方差矩阵公式;然后利用协方差椭球对该协方差矩阵的性质做了详细分析,明确了探头间构型、星图融合参考系选择以及各探头星点选取等因素对星图融合姿态确定精度的影响关系;仿真结果表明,最优构型的多探头星敏感器通过星图融合可输出同精度的三轴姿态,并且光轴精度较单探头星敏感器提升了30%。研究成果对多探头星敏感器的工程研制具有理论指导意义。     

一种双星敏感器联合确定船体姿态方法

单星敏感器横滚测量精度偏低会影响航向测量精度,采用双星敏感器联合的方法可以有效提高船姿精度,其中关键是蒙气差修正和联合姿态确定方法。目前蒙气差修正一般采用基于恒星观测矢量的修正方法,联合姿态确定方式一般采用基于双视轴指向联合的方法。以上在蒙气差中需要引入初始船体姿态,而且修正迭代解算过程复杂,再者单星敏感器需要识别4颗以上恒星才能获得视轴指向,如果其中一个无法获得视轴指向,双星敏感器无法解算船体姿态。针对上述情况,提出了一种新的双星敏感器联合确定船体姿态的方法,包括基于恒星参考矢量的蒙气差修正方法和基于恒星观测矢量的联合姿态确定方法。基于该方法蒙气差修正不再需要初始船姿,而且在理论上只需双星敏感器各自至少有1颗可识别恒星,即可完成船体姿态角的解算。试验结果表明该方法降低了双星敏感器获取船体姿态角所需的可识别恒星数量限制,提高了双星敏感器获取船体姿态的能力。     

基于姿态关联帧叠加的星图信噪比增强方法

星敏感器在近地面使用时,由于星光信号容易被强烈的天空背景噪声所淹没,导致其无法在白天正常工作。为了提高星敏感器白天探星的能力,提出了一种星敏感器动态条件下基于姿态关联叠加提高星图信噪比的方法。该方法利用陀螺仪测量出与参考帧相关的每个星图帧之间的姿态变化量,实现不同帧星图之间的平移和旋转变换(即相关变换)。最后,将相关帧星图叠加到一起,使星图信噪比得到■倍增长(n为关联星图的帧数),从而提高全天时星敏感器在白昼条件下的工作能力。仿真结果表明:该方法能够有效地提高星图信噪比,且受惯性器件误差和安装角误差的影响较小,具有较强的可行性。     

位置未知条件下惯性/星敏组合系统海上牵引启动方法

针对地理位置未知条件下惯导系统海上启动对准实现问题,提出了基于锚泊条件的惯性/星敏组合导航系统位置未知条件下的海上牵引启动方法。建立了惯性/星敏信息的迭代滤波模型,研究了惯性?/星敏牵引启动方法、分析了自主启动精度、提出了自主启动实现方案。通过惯性基准和星敏感器联合使用,惯性基准辅助星敏感器进行观星定位,观星定位数据对惯性基准平台反馈修正,基于惯性/星敏位置和姿态信息在牵引启动下的相互校准方式,组合系统经过多次迭代修正,构造高精度水平基准,完成惯性/星敏组合系统海上位置未知条件下纬度、经度自估计,实现船舶对准后的正常航行。实验结果表明,采用所提出的惯性/星敏组合系统牵引启动方法获得的经纬度自估计精度与星敏感器的测量精度相当。     

在飞航导弹中用星敏感器修正捷联惯导陀螺漂移

星敏感器是一种高精度的姿态测量装置。研究了星敏感器和陀螺的特点,对星敏感器工作原理和修正陀螺漂移技术进行了原理分析。在不利用外界提供的姿态和位置信息的情况下,采用卡尔曼滤波的信息融合算法,建立组合导航系统的状态方程和量测方程,利用星敏感器输出的载体相对于惯性空间的姿态信息来修正捷联惯导的陀螺漂移。设计飞航导弹的典型飞行轨迹,通过数学仿真,对上述算法的有效性进行了验证,结果表明星敏感器能够有效地补偿捷联惯导由于陀螺漂移带来的误差,明显提高了导航定位精度。     

一种小卫星高精度自主定轨/定姿一体化新方法

基于传统小卫星对轨道和姿态参数确定采用分别计算的复杂模式,提出了一种利用地磁场和天文信息同时确定卫星轨道和姿态参数的新方法。首先通过分析小卫星轨道动力学J2模型和卫星姿态动力学模型,建立系统状态方程。其次将三轴磁强计与地磁场模型参考值的矢量作差,分析微分差值与状态变量的数学关系,建立定位/定姿观测方程。利用星敏感器提供的高精度姿态信息,建立定姿观测方程,同时利用星敏感器间接敏感地平观测折射恒星,建立定位观测方程。最后提出基于信息融合的先进滤波算法,并通过对多种导航模式进行数值仿真及结果分析,论证所设计一体化方法提高了系统定轨/定姿的精度和可靠性。     

基于双微机电惯性测量单元的船用自主水平姿态估计算法

针对船载自主水平姿态基准系统需求,提出了一种基于双微机电惯性测量单元的方案。其中一套固联船体,另一套为转位方式,以两套微机电惯性测量单元间相对速度、位置、姿态为观测约束,基于状态变换卡尔曼滤波实现不依靠外界辅助信息的自主式水平姿态测量算法,减弱了陀螺、加速度计零偏误差对水平姿态精度的影响。通过采用三轴角运动模拟转台进行了海况模拟实验验证。无舰船主惯导提供舰船运动信息情况下,虽然船的质心水平机动会影响水平姿态的绝对精度,但不会影响船上不同位置点之间的相对姿态测量精度;而当舰船上有高精度主惯导系统可提供速度参考时,即使有常值速度偏差和舒勒周期速度误差,仍可实现绝对姿态测量精度优于0.02°(1σ)。所提出的姿态测量方案可以在海况恶劣且无外界参考信息的情况下建立高性价比的全船统一姿态参考基准。     

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

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

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


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).