无陀螺捷联惯性导航技术

无陀螺捷联惯性测量装置与传统捷联惯导系统的主要区别是角速度的获取方式不同,角速度的解算精度是无陀螺捷联惯性导航系统的核心问题,决定了系统的性能及实际应用的可行性。本文剖析了无陀螺捷联惯性测量装置的误差来源,建立了无陀螺捷联惯性导航系统角速度解算数学模型,并重点探讨了加速度计元件误差对角速度解算精度的影响。进行了无陀螺捷联惯性测量装置试验,结果表明,尽管计算误差较大,但无陀螺捷联惯性测量装置可以反映出运动平台的角运动规律,实际应用中对加速度计精度和计算机速度要求较高,另外应寻找更好的算法尽量补偿角速度解算误差。     

权值动态分配的GFSINS多角速度加权平均融合算法

针对单一角速度算法解算精度不高,提出将最优加权平均法应用于GFSINS角速度解算中。针对实际应用中各种算法无法准确获取到角速度解算误差方差的情况,给出了一种对方差进行实时估计,再基于最优权值分配原则来动态分配权值的加权融合算法;对积分法、开方法、微分法和对数法四种角速度算法进行了有效融合,并完成了仿真。仿真结果表明,即使精度最低的积分法参与加权融合,也能进一步提高融合精度,同时还验证了加权平均法在系统角速度解算中的有效性,融合后,角速度精度比对数法提高了一倍。     

无陀螺惯性测量单元设计及仿真分析

在简述无陀螺惯性测量单元(GF-IMU)测量角速度的一般原理的基础上,提出了一种九加速度计的配置方案及其角速度解算方法.该算法利用加速度计输出的测量信息直接求得所测量的角速度的绝对值,没有积分运算引入的迭代误差,在一定程度上提高了载体角速度的解算精度.针对该九加速度计配置方案在解算载体角速度时,存在需要开方及符号误判等不足,提出一种解算角速度的微分算法,彻底消除解算过程中开方及符号判断带来的误差,提高了系统精度和实时性.仿真分析证明了方案的可行性和算法的有效性.     

基于十二加速度计的GFSINS安装误差标定及补偿

针对高速旋转弹姿态解算精度不高的问题,提出用微小型加速度计构成无陀螺微惯性测量单元GFMIMU,为旋转弹提供完整的导航信息.设计并制作了一种十二加速度计的惯组实物模型,并在此基础上分别提出了方向安装误差角及位置安装误差的标定方案.根据存在安装误差的角速度解算精确数学模型,采用一种优化算法--最小二乘优化迭代方法,可以精确地补偿加速度计的安装误差.仿真结果表明,经补偿后误差被有效抑制,在20 s内姿态角的累积误差均能保证在以内,满足高精度姿态解算要求,具有一定实际应用价值.     

单陀螺多加速度计捷联惯性导航解算方法

六加速度计无陀螺惯导系统误差随时间发散比较严重。为了有效提高导航系统精度,提出了一种单陀螺仪多加速度计(五加速度计)的捷联惯性导航解算方法。该导航解算方法通过合理配置5个加速度计和1个陀螺仪,可不经积分而直接解算角速度,完全消除了加速度计输出方程中角加速度项的影响,能使在姿态和位置解算时分别减少1次积分,从而有效抑制误差随时间发散。给出了单陀螺多加速度计捷联惯导姿态和位置解算原理的理论推导过程,并对该导航解算方法进行了仿真。在仿真时间为80 s时,与无陀螺惯导相比,该方法的姿态解算和位置解算精度均提高了60%以上。     

GFSINS姿态角速度双路组合方案设计

针对无陀螺捷联惯导系统解算载体姿态角速度精度不高的系统瓶颈,分析现有的一种九加速度计配置方案角速度信息输出特征,提出通过改进解算方法,在载体上同时获取两套相对完整而独立的姿态角速度信息从而构成姿态角速度双路组合方案的设想,并且在该组合方案中引入反馈型自适应神经网络,通过设置合理的神经网络学习周期,提供更逼近真实值的加速度计输出,结合包含加速度计误差少、没有加速度计误差积累两种优点的姿态角速度辅助算法,以获得更高精度的载体姿态角速度输出,进一步提高无陀螺捷联惯导系统导航精度.仿真结果验证了该设计思想可行性,并且证明了角速度双路组合方案相对普通角速度解算方案在精度上的优势.     

一种用于无陀螺捷联惯导系统的角速度融合算法

无陀螺捷联惯导系统(GFSINS)是用加速度计的合理空间组合解算出载体的角速度。载体角速度的解算精度是GFSINS的技术关键。在分析GFSINS九加速度计配置方案的基础上,提出一种新的角速度融合算法,消除该方案解算过程中开方计算及符号判断带来的误差。该算法还明显抵消了加速度计输出中包含的常值零点偏移误差和温度漂移误差等,具有实时性好、计算量小、通用性强的优点。仿真计算表明该算法可行,并能在一定程度上提高系统解算精度。     

无陀螺仪惯性系统构型中安装误差分析与标定

针对一种立方构型的无陀螺仪惯性系统，分析了构型安装误差对惯性系统测量精度的影响及构型误差的标定方法。首先分析了立方构型下载体姿态与加速度的解算方法，然后给出了安装误差对于系统输出的影响，最后通过分析安装误差的敏感性因子给出了一种有针对性的标定方法。理论分析与仿真计算表明，惯性系统角速度与线速度输出是加速度计输出的线性组合，构型安装误差对无陀螺仪惯性系统测量精度的影响非常明显，在给出的标定方案中，方向安装误差标定精度较位置安装误差要高，而且转台角速度对位置安装误差的标定精度影响很大，对方向安装误差影响较小；无陀螺仪惯性系统中安装误差不可忽略，必须进行标定，据此提出的一种标定方案简单、切实可行。     

一种惯性测量装置火箭橇试验误差分离方法

惯性测量装置火箭橇试验可以提供大过载环境,为标定加速度计高阶误差项系数提供必要条件。为了利用位置外测信息准确辨识加速度计高阶误差项系数,提出一种针对火箭橇试验的惯性测量装置误差模型参数辨识方法。在给出加速度计高阶误差模型后,该方法通过建立惯性测量装置位置遥测误差和误差模型参数间的线性函数关系构建位置环境函数,之后利用最小二乘法估计出误差系数数值。仿真结果表明,该方法在辨识加速度计高阶误差项参数上具有很高的精度;应用该方法对一次火箭橇试验数据进行误差分离后,通过对解算结果进行显著性分析,证明了辨识出的高阶误差项系数的有效性。     

旋转式光学陀螺捷联惯导系统的旋转方案设计

在光学陀螺捷联式惯性系统中,利用系统旋转补偿技术可对陀螺组件和加速度计组件的输出误差进行调制,从而抑制系统的误差发散,提高导舷精度.通过分析惯性测量组件的误差模型和旋转式捷联系统误差传播方程,解释了旋转误差补偿的机理.针对惯性测量组件输出误差的特性,设计单轴正反转停和双轴转位的系统旋转方案.在摇摆状态下分别对无旋转、单轴和双轴三种方案进行长时间导航仿真,对旋转补偿误差的能力进行了比较.结果表明:旋转能够抑制长期的定位误差发散,在角运动状态下旋转系统能比无旋转系统保持更好的姿态精度.     

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

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

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

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