参数辨识双位置对准改进算法

对捷联惯导参数辨识双位置对准算法进行了研究。指出现有的绕方位轴转90°得到第二位置的方案不是最优方案,绕方位轴转180°作为第二位置是最优转位方案。为在参数辨识对准中使用上述最优转位方案,提出了回溯导航算法,给出了回溯参数辨识初始对准算法。根据回溯参数辨识对准算法对参数辨识双位置对准算法进行了改进。最后采用光纤捷联惯导对改进前后的对准算法分别进行了半实物仿真试验。试验结果表明,改进算法有效地提高了对准精度和测漂精度。     

等效转动矢量在捷联惯导系统二位置对准中的应用

捷联惯导系统的初始对准精度受加速度计零偏和陀螺漂移的限制。为减小对准误差,引入等效转动矢量对二位置对准方法进行了研究。通过对捷联惯导系统的静基座误差方程作Lyapunov变换,得到转动过程中加速度计零偏和陀螺漂移误差的传播方程。当系统由第一位置转到第二位置时,利用等效转动矢量的方法得到方程的状态转移矩阵,并由第一位置对准的约束条件推导出加速度计零偏和陀螺漂移的解析解。当二位置对准采用绕方位轴转动时,分析解的稳定性得出最优转动角度是180°。     

一种绕任意轴旋转的捷联惯导系统多位置初始对准方法（英文）

基于单轴旋转SINS的多位置对准方法通常要求旋转轴和方位轴重合,轴重合误差会影响多位置对准精度。针对这一问题,提出了一种可绕任意轴旋转的多位置初始对准方法,几何分析指出,该旋转轴可沿任意方向,不需要与方位轴重合,它只需绕某任意旋转轴转过两个已知角度即可解算出IMU的零偏,提高对准精度。仿真实验显示,提出的绕任意轴旋转的多位置对准方法相比单个位置下的解析粗对准方法的对准精度显著提高,航向角误差由5'降低到0.2'以下。     

捷联惯导系统双位置快速抗干扰对准方法

经典捷联惯导系统双位置参数辨识精对准采用水平速度作为量测值,对部分参数辨识时间较长,且收敛速度受递推计算初值影响。提出了一种利用单位时间的水平比力作为量测值的双位置快速抗干扰对准方法,并忽略误差构成的高阶项,可以将三次曲线转换为线性模型进行估计,提高了初始对准的收敛速度。双位置对准试验结果表明,双位置快速抗干扰对准方法算法计算量小,算法简单,对准速度快,抗干扰性强,对准精度高且不受递推计算初值的影响,具有较强的工程实用价值。     

基于双速度模式的多级惯性系动基座粗对准算法

由于传统载体系速度(V~b)辅助惯性系对准算法(IA)中的近似忽略项,不可避免的存在模型误差问题将会影响捷联惯性导航系统对准精度,且会随载体速度增大而更加明显。对传统多级IA算法进行改进,提出基于双速度模式的惯性系初始对准算法来解决其模型误差问题且不牺牲其对准精度的稳定性,同时利用回溯计算来逐步提高对准精度,保证对准快速性。最后,通过仿真和车载半实物仿真试验进行验证。车载半实物仿真下,传统V~b辅助IA算法、传统多级IA算法以及改进多级IA算法在120 s的航向对准精度分别为-0.779°、1.747°和-0.314°,5次半实物仿真航向对准精度均方误差分别为0.7082°、2.6482°和0.3392°。试验结果表明采用所提出的双速度模式的多级惯性系粗对准算法具有优越的性能。     

基于参数识别视速度的双矢量粗对准方法（英文）

针对传统基于视速度双矢量粗对准中,由于传感器随机噪声的影响,存在对准精度差,收敛速度慢的缺点,提出了一种新型自适应Kalman滤波的参数识别粗对准方法。该方法通过对视速度运动进行建模,设计采用自适应Kalman滤波对模型参数进行参数识别,从而有效地消除视运动中的随机噪声,提高粗对准的精度和收敛速度。由于自适应滤波的特点,新方法不需要对传感器误差进行统计,使其在实际系统中具有更加广泛的应用价值。针对双矢量粗对准的计算特点,设计了一种矢量重构算法,从而尽可能地规避双矢量共线性问题,加快了粗对准的收敛过程。仿真与转台实验表明,与传统方法对比,新方法在相同的对准时间内具有更高的对准精度,在相同的对准精度下,具有更高的收敛速度。转台实验的最终对准精度为-0.1391°,标准差为0.012°。     

平台惯性导航系统初始对准技术研究及调试

本文研究平台惯性导航系统初始自对准技术，实际上是采用任意双位置平台方位锁定。文中首先介绍自对准原理，然后进行详细的理论分析，最后以某型号平台惯性导航系统为应用，介绍数字精调平回路、平台方位锁定回路的设计及初始自对准调试情况。对应于系统随机漂移为０．０５°／ｈｒ的指标，平台方位自对准精度达到１７．４２角分（１σ），且重复性较好，自对准调试结果得到了光学传递对准的检验。随后进行了车载导航实验，导航精度达到４．５５ｎｍ／ｈｒ（１σ）。     

抗扰动的捷联惯导系统回溯参数辨识对准法

针对捷联惯导系统参数辨识对准法航向角估计时间长且收敛速度受到计算初值影响大的问题,提出了将对准用的陀螺和加速度计采样数据存储后按时间进行正向和逆向排列,形成正向序列和逆向序列。由正向导航的解算公式推导出逆向导航的解算公式,利用公式的一致性提出了回溯误差参数辨识对准法。利用抗扰动的惯性系对准算法为回溯参数辨识法提供对准初值,应用回溯参数辨识法对误差参数反复进行逆向和正向辨识后对对准结果进行修正,能够有效提高对准精度。实验结果表明,利用相同时间的采样数据,该算法航向角的对准误差仅为改进前算法的25%。     

激光陀螺速率偏频系统初始对准方法研究

针对速率偏频激光陀螺转动的工作特性，建立了新的系统误差方程，并在静基座条件下，采用卡尔曼滤波方法实现了系统的初始对准。仿真结果表明，这种初始对准方法对准速度较快，水平对准精度比抖动偏频系统提高了一个数量级，方位对准精度也有较大提高，并且对两个水平加速度计随机常值偏置有很好的估计效果。     

一种改良Kalman滤波参数辨识粗对准方法

针对传统基于g信息的粗对准的捷联惯导系统中,受传感器噪声的影响,存在效视运动无法提取和双向量共线的缺点,提出了一种基于改良Kalman滤波的参数辨识粗对准方法。该方法通过构建视在重力在初始载体系中的映射模型,利用改良Kalman滤波进行模型参数辨识,然后通过识别参数重新构建视在重力在初始载体系中的映射,解决了由于传感器噪声导致有效视运动无法正常提取的缺点。利用识别参数具有随估计次数增多得到优化的特点,构造初始时刻和最终时刻向量,避免双向量共线问题。利用改良Kalman滤波算法的自适应特点,优化参数识别精度与速度。转台实验表明,采用改良Kalman滤波方法航向对准精度为-0.0414°,标准差为0.041°,而传统RLS方法得到的航向精度为-0.0738°,标准差为0.128°。由此可知,本文提出的方法性能更优。     

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

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

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

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