基于SVD方法的INS传递对准的可观测性能分析

在惯导系统一般的误差动态方程和速度观测方程的基础上,建立了姿态传递对准所必需的弹、舰相对姿态误差方程和观测方程.介绍了基于动态系统可观测性矩阵奇异值分解的状态变量可观测度的分析方法.用奇异值分解的方法,对同时采用速度和姿态传递的INS对准模型,分析了系统变量的可观测性和可观测度,为对准方程的可观测性结构分解和误差估计性能的改善提供了必要的基础.     

基于SVD方法的lNS传递对准的可观测性能分析

在惯导系统一般的误差动态方程和速度观测方程的基础上，建立了姿态传递对准所必需的弹、舰相对姿态误差方程和观测方程。介绍了基于动态系统可观测性矩阵奇异值分解的状态变量可观测度的分析方法。用奇异值分解的方法，对同时采用速度和姿态传递的INS对准模型，分析了系统变量的可观测性和可观测度，为对准方程的可观测性结构分解和误差估计性能的改善提供了必要的基础。     

基于谱分解的SINS传递对准可观测度分析方法

针对舰载武器捷联惯导系统的快速精确传递对准问题,对"速度+姿态"匹配进行了简化,给出"速度+航向"匹配的传递对准模型。在此基础上,提出了一种基于Hermite矩阵谱分解的可观测度分析方法,该方法仅需要对阶数与系统状态维数相同的Hermite矩阵进行谱分解,算法简单。最后,应用该方法分别对三种匹配(速度匹配、"速度+姿态"匹配和"速度+航向"匹配)的系统状态进行可观测度分析。仿真试验结果表明:基于谱分解的可观测度分析方法是正确有效的,"速度+航向"匹配的对准精度在10 s内优于4′,能够实现捷联惯导系统的快速精确传递对准。     

局部可观测理论在惯导系统快速传递对准中的应用

鉴于传统的全局可观测理论很难定量分析时变系统的可观测性,从机载导弹传递对准姿态角误差的可观测性出发,首次将局部可观测性理论应用于"速度 姿态"匹配的快速传递对准中。将条件数的概念引入局部可观测矩阵,定量地计算出在三种不同机动方式下,局部可观测矩阵的条件数,用以表征机动对准过程中系统的局部可观测度。仿真结果证明,提高姿态角误差可观测度的最佳载体机动方式为机翼摇摆运动。该研究结果为实现机载武器动基座快速精确对准技术在工程中的应用提供了依据。     

基于SVD的SINS多位置对准可观测性分析

针对捷联式惯导系统多位置对准可观测性的问题,以捷联式惯导十状态误差方程为研究对象,利用卡尔曼滤波分别得到固定位置对准与二位置对准以及二位置对准与三位置对准时的方位失准角估计误差收敛情况。在分段线性定常系统理论的基础上,利用奇异值分解的方法,分别对固定位置对准以及二位置对准时的系统各状态变量的可观测性进行分析。仿真结果表明,三位置对准和二位置对准时的方位失准角估计误差达到的稳态误差是一致的,并且多位置对准能够改善捷联惯导系统各状态变量的可观测度。该研究结果不仅为确定并提高捷联惯导系统各状态的可观测度提供了途径,而且为捷联式惯导系统的可观测性达到最佳以及捷联惯导系统对准精度的快速提高提供了理论基础。     

基于可观测性分析的光纤陀螺SINS/测速仪组合系统技术

为了简化求可观测度的过程,便于导航时利用状态变量可观测性分析结果进行姿态解算。在仔细研究了分段定常系统(PWCS)的可观测性分析方法,以及基于奇异值分解的系统状态可观测度分析方法的基础上,提出了一种改进的状态可观测性分析方法,该方法在求可观测度时,不需要量测值,利用可观测性矩阵奇异值分解后的V阵即可进行可观测性分析;利用该理论得到的结果,在SINS/测速仪组合导航解算过程中只对可观测度好的俯仰角误差和横滚角误差进行反馈校正,不反馈可观测度较差的航向角误差,仿真结果证明,这种有选择的反馈校正方法可以保证姿态角在导航开始几秒钟内就快速收敛。     

基于伪距的北斗双星/SINS组合导航系统可观测度分析

在连续的分段式定常系统模型基础上,证明了通过简化分段式定常系统的可观测矩阵来分析其离散化模型可观测性能的可行性;介绍了一种基于奇异值分解的系统状态可观测度分析方法.通过对一个分段间隔的条带化可观测矩阵进行奇异值分解来判断每一个奇异值对应的系统状态及每一个系统状态对应的可观测度;根据北斗双星定位系统的双伪距、高精度原子钟及气压高度信息,建立了北斗双星/SINS伪距组合滤波模型,并将可观测度分析方法应用于该模型中.仿真结果表明:在伪距组合模式下对简化的分段式定常系统可观测矩阵进行奇异值分解的方法能有效反映北斗双星/SINS组合导航系统的可观测程度.     

基于降维滤波器的SINS快速初始对准算法（英文）

引入惯性元件信息观测的初始对准方法能够快速提高SINS的对准速度,但同时存在滤波计算量大、系统噪声与观测噪声相关以及观测值中的高频噪声影响滤波精度的问题。针对这些问题提出了一种捷联惯导快速初始对准降维滤波器设计方法,通过剔除不可观测量和合理选取状态量以降低状态方程维数,并推导了观测方程,在采用低通滤波器对惯性器件原始信息预处理基础上应用噪声相关下的Kalman滤波进行状态估计。理论分析和试验结果表明,新方法提高了对准速度,减少了计算量,水平姿态角收敛速度提高了90%,计算量减少了83.33%,并可有效抑制高频噪声对状态估计的影响。     

基于惯性参考系基准的快速传递对准方法

研究了以惯性参考系为基准的新型传递对准方法。推导了计算惯性坐标系和计算体坐标系传递对准动态误差模型,并给出了相应的"速度+姿态"观测方程。基于惯性参考系的对准方法通过链式法则将子惯导输出的姿态矩阵描述为三个变换矩阵之积,其中两个变换矩阵通过对准时间和主惯导提供的位置信息可得到精确求解,剩余的变换矩阵(子惯导体坐标系至惯性坐标系间的变换矩阵)通过子惯导陀螺仪的输出进行解算,其误差通过传递对准估计得到的失准角进行补偿。对提出的两种对准方法进行摇摆实验验证,方位对准误差优于4’(1)。与传统基于导航坐标系的方法相比,基于惯性坐标系的方法直接将误差定位到惯性坐标系上,具有算法简便的特点。     

车载捷联惯导系统静止条件下的初始对准方法研究

根据捷联惯导系统(SINS)的误差模型，从提高初始对准中卡尔曼滤波估计的可观测性和可观测度的角度出发，提出了同时利用加速度表的输出和速度误差量这两种信息作为卡尔曼滤波的测量量，对陆地车辆捷联惯导系统在静止条件下进行初始对准。可观测性定性分析、可观测度定量计算和仿真结果都表明，充分利用外部可观测信息，可以提高系统的可观测性和可观测度，加快初始对准的速度，减小估计误差。该方法已经在实际的SINS中得到应用。     

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

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

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

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

Preface: Current research progress on mechanics of high speed rail

The construction and operation of high-speed rail（HSR）grid within the past two decades in China,in terms of the scale,may be possibly comparable to any national-wide construction in the history of China,even the Great Wall.By counting railways with commercial train service at the speed of200 km/h,China has the world’s longest HSR network with over 19 369.8 km of track in service today and this number     

ACTA MECHANICA SINICA 2014 ANNUAL SUBJECT INDEX

正~~     

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