考虑机抖激光陀螺信号滤波特性的圆锥算法修正

捷联惯导系统一般用圆锥算法来补偿其圆锥运动漂移,标准的圆锥算法是以理想的陀螺信号为输入设计的,而机抖激光陀螺常用数字低通滤波器来滤除抖动偏频信号,由于滤波器的非理想性,滤波后信号的幅频特性发生畸变,引入了较大的圆锥算法误差。参考标准圆锥算法误差公式,用相对圆锥误差分析方法比较研究了算法的误差变化特性,证明误差大小与滤波器通带特性有关。基于经典圆锥运动,推导了数据滤波后修正的圆锥算法公式。修正算法考虑了滤波器幅频特性的影响,补偿了滤波器引入的圆锥算法误差。仿真表明:修正算法在保证滤波器较小延时的同时,能明显减小算法精度损失,提高姿态算法整体精度水平。     

激光捷联系统中的优化圆锥误差补偿算法

激光捷联系统中采用低通滤波器消除激光陀螺角增量输出中机械抖动引入的干扰信号,同时也对陀螺敏感的外界惯性输入角速度信号进行了频率整形,产生了视在圆锥误差,此时传统圆锥补偿算法未考虑滤波器影响补偿精度严重降低.针对本系统采用了31阶低通滤波器对陀螺的角增量输出整形,分析了其引入的视在圆锥误差,基于滤波器的频率特性,采用五子样圆锥误差补偿算法,即在旋转矢量更新周期内有五个陀螺采样信号,可以构成四种不同时间间隔的陀螺输出角增量信号的叉积,利用这些叉积的线性组合更新旋转矢量.仿真结果表明,对经过滤波器整形的陀螺输出角增量进行补偿,优化的圆锥补偿算法的补偿精度明显优于传统圆锥补偿算法,使系统姿态角的精度提高了两个数量级.     

多凹点滤波器在LINS中的应用

分析了激光捷联惯性导航系统中的激光陀螺频谱特性,指出激光捷联惯性导航系统中三个激光陀螺的抖动频率会互相干扰。提出一种多凹点有限冲击响应低通滤波器频域设计方法,通过指定三个抖动频率点幅值为小量,在幅频特性曲线上产生三个凹陷点。该滤波器采用标准滤波算法,同时实现了对激光捷联惯性导航系统三个激光陀螺的滤波和机抖信号的陷波。给出了一个24阶多凹点滤波器设计实例,幅频特性分析说明,由于三个陷波点的集中,该滤波器还可以在机抖频率可能存在的频段上,提供额外的80 d B以上的衰减,可防止激光陀螺抖动频率漂移后陷波性能下降。试验结果表明,多凹点滤波器输出数据波动幅度比原激光捷联惯性导航系统使用的滤波器结果减小5~8倍,三轴摇摆试验的导航精度也有一定提高。     

四频差动激光陀螺三种误差模型对寻北精度的影响

分析了四频差动激光陀螺漂移信号的特性,将陀螺输出的漂移误差信号分为常值漂移误差和随时间变化的一次项、二次项漂移误差,并据此建立陀螺漂移误差模型,分别对陀螺漂移进行零次拟合、一次拟合及二次拟合.针对这些模型结合寻北推导了误差的补偿算法,并通过寻北实验精度比较,验证了不同误差模型的补偿效果.实验结果表明,就本文实验所用陀螺,含二次项误差的模型寻北精度较高,使寻北精度从零次拟合模型的1密位降低到0.5密位.     

考虑幅相频畸变的激光惯导圆锥补偿算法设计

激光惯导系统中角速度测量元件激光陀螺通过抖动轮安装在台体上,从频域角度看,其输出信号与台体真实的运动在幅值和相位存在微小畸变,这种畸变会引起圆锥补偿算法误差。为减弱该误差项,提出一种考虑激光陀螺幅相频特性的圆锥补偿算法的设计方法。在传统圆锥补偿算法的基础上,改进的算法在设计过程中针对激光陀螺的谐振频率,结合惯导系统减振器的设计带宽,采用加权最小二乘的FIR滤波器设计思路,优化圆锥补偿算法的相关系数。最终,仿真试验结果表明在标准圆锥运动环境下采用改进的圆锥补偿算法,相对算法误差由31%减小到0.04%,为高精度激光陀螺惯导系统圆锥补偿算法的优化设计提供了理论支撑。     

激光陀螺POS惯性数据滤波及时延补偿

激光陀螺作为位置姿态测量系统(POS)的核心传感器,其精度直接决定激光陀螺POS系统精度,围绕机抖激光陀螺信号去噪的需求,基于FIR滤波器建模的方法,设计了FIR数字低通滤波器;针对滤波器导致的信号延迟问题,根据FIR数字滤波器群延迟特性建模,提出了激光陀螺POS数据时延补偿方法。静态实验结果表明,设计的FIR数字低通滤波器降低了激光陀螺抖动噪声功率达80dB。进一步通过飞行实验表明,提出的方法降低了激光陀螺POS系统姿态角误差,与POS/AV510相比航向角误差由0.017°降低到0.01°,俯仰角误差由0.007°降低到0.005°,横滚角误差由0.016°降低到0.005°,满足了机载InSAR对激光陀螺POS精度要求。     

基于状态约束的MIMU/磁强计组合姿态估计滤波算法

构建了基于MEMS技术的陀螺、加速度计、磁强计及空速计组合的微小型飞行器姿态测量系统.研究了基于四元数的扩展卡尔曼滤波算法.取姿态误差四元数和陀螺随机漂移构建状态向量,通过误差四元数微分方程和陀螺随机误差模型建立卡尔曼滤波状态方程,采用速度信息实时补偿加速度计输出值得到重力矢量,利用重力矢量估计水平姿态,通过滤波补偿姿态误差,降低了对陀螺的精度要求.将状态向量之间的约束方程作为伪量测方程引入到量测模型中,解决了由于状态向量相互约束导致的滤波发散和奇异.动态飞行滤波噪声的自适应调整增强了系统性能.仿真和实验表明,该滤波算法能够有效避免系统的漂移,提高系统测量精度和稳定性.     

恒速偏频激光陀螺系统静基座初始对准中等效东向陀螺漂移估计

以单轴恒速偏频激光陀螺系统为研究对象,在分析IMU传感器误差的基础上,建立了合理有效的静基座初始对准滤波器模型。针对系统连续旋转运行的特性,提出了简洁适用的滤波器估计误差检验方法,利用自主设计的原理样机验证了恒速偏频技术的实际可行性,对滤波算法进行了实验测试。实验结果表明,初始对准滤波算法能够稳定有效地估计IMU传感器误差,且等效东陀螺漂移估计精度优于0.0004(°)/h,该系统具有很高的工程应用潜力。     

基于前向线性预测算法的光纤陀螺零漂的神经网络建模

在详细分析光纤陀螺零漂的基础上,提出了先用滤波算法对光纤陀螺信号进行预处理,然后采用RBF神经网络对滤波后的信号进行建模的方法。针对光纤陀螺信号特点分别采用FLP算法、小波滤波算法、解相关变步长LMS自适应滤波算法对其进行了预处理,比较三种滤波方法,小波滤波算法效果优于其它两种预处理方法,但针对基于预处理后的陀螺信号采用RBF神经网络进行建模时,小波滤波预处理后的信号在建模精度上却是最差的,而对FLP算法滤波后的信号进行RBF建模,建模精度提高了两个数量级。结果表明:基于FLP算法的RBF神经网络在光纤陀螺中的建模是有效的,可大大提高建模的精度。     

基于自适应联邦滤波的传递对准算法

为了提高传递对准的对准速度和对准精度,研究了基于自适应联邦滤波的"速度 姿态角"传递对准的算法,在状态变量中加入了安装误差角和挠曲变形角,在算法中进行了估计并补偿。为了保证滤波的实时性,采用联邦滤波的方法,分别建立了速度匹配子滤波器、姿态角匹配子滤波器和主滤波器的模型,给出了状态方程和量测方程。在子滤波器中利用模糊控制器对噪声特性进行了自适应调整以解决系统噪声和量测噪声是未知情况下滤波发散或者精度不高的缺点从而增强系统的鲁棒性。最后在载体匀速直线加三轴摇摆的运动模型下进行了仿真,结果表明该方法能够有效地估计安装误差角和部分挠曲变形角,并且能够以一定精度完成初始对准。     

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

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

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

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