连续旋转的光纤陀螺全温失准角快速建模与补偿方法

在全温范围内应用的光纤陀螺,其输入轴失准角随温度的变化是影响光纤陀螺惯性系统性能的重要指标之一。特别是在大角速率或者高精度应用时,失准角的变化误差甚至超过零偏漂移误差和标度因数误差。采用温度补偿技术是一种提升光纤陀螺温度性能的有效方法,其中建立精确的温度模型是关键。提出了一种连续旋转的光纤陀螺全温失准角快速建模补偿方法。基于单轴速率转台的连续旋转,可以有效识别光纤陀螺失准角在全温范围内的变化拐点,提高建模和补偿的精度。试验结果表明,某型光纤陀螺全温输入轴失准角变化约14″,补偿后全温输入轴失准角变化小于1″,精度提高了一个数量级以上。在高精度光纤陀螺惯性系统中,该方法可用于指导光纤陀螺失准角的实时温度补偿技术研究及工程实现。     

用于改善光纤陀螺全温标度因数变化的环圈制作工艺

全温标度因数变化大小直接影响光纤陀螺的应用效果。随着中高精度光纤陀螺在复杂工况下的广泛应用,对陀螺标度因数变化的要求也在不断提高。光纤长度和光纤环平均直径在全温范围内的变化量是全温标度因数变化最大的影响因素。本文介绍了控制该变化量的常规方式,并提出通过不同排纤方式和不同光纤环窗口比的方法来进一步减小这种变化量的方法。通过仿真和试验证实:某型中等精度光纤环在不同制作工艺条件下全温标度因数变化量有将近300×10~(-6)的差值。     

光纤陀螺IMU全温三方位速率/一位置标定及分段线性插值补偿方法

针对光纤陀螺惯性测量单元全温动态环境下测量误差问题,提出一种全温三方位正反速率/一位置标定及分段线性插值补偿方法,建立了光纤陀螺惯性测量单元误差模型,在每个恒温点设计三方位正反速率/一位置标定方案。采用分段线性插值算法实时补偿系统零偏和标度因数温度误差,系统全温环境下的测量精度提高5倍左右。车载实验结果表明,采用该方法后系统4200 s纯惯性姿态测量误差小于1°。     

光纤陀螺标度因数分段标定的工程实现

基于光纤陀螺标度因数存在非线性和不对称性,从工程实用角度出发,介绍了光纤陀螺速率标定方法,给出了光纤陀螺标度因数分段标定的工程实现方案,阐述了速率标定测试数据最优分段原则,找出了使分段标定后不连续光纤陀螺模型变成连续模型的解决措施.多个光纤陀螺速率标定的研究结果显示:分段标定能较大程度上减小光纤陀螺标度因数非线性度;利用实测数据进行离线导航计算的结果表明:对光纤陀螺标度因数分段标定能减小标度因数误差对导航精度的影响,文中分段标定方法切实可行.     

大动态光纤陀螺标度因数双线性插值补偿模型

标度因数的温度特性和非线性是影响大动态光纤陀螺精度的重要因素,在航空、航天等对动态精度要求较高的应用场合需要对标度因数的温度特性和非线性误差进行补偿。通过对大动态光纤陀螺标度因数的误差分析,得出环境温度和输入角速率是影响标度因数误差的主要因素。建立了一种基于双线性插值的补偿模型,对大动态光纤陀螺的标度因数的温度特性和非线性进行综合补偿。在温度范围为-40℃~+60℃、角速率范围为0~7200(°)/s的条件下,标度因数误差由补偿前超过1.3′10~(-3)降低为小于5′10~(-6),标度因数精度提升了两个数量级,验证了补偿模型的有效性。补偿算法复杂度低,易于工程实现。     

光纤陀螺标度因数迟滞模型分析与补偿技术

高精度惯性导航系统对由温度引起的光纤陀螺标度因数变化指标提出了很高的要求。采用温度补偿技术是一种提升标度因数性能的有效方法,其中建立精确且普适的温度模型是关键。提出并分析了光纤陀螺温度与标度因数模型的迟滞现象。通过分析和试验表明,标度因数模型的迟滞现象是由光纤陀螺结构的热不均匀性造成的,采用多温度点采样来修正标度因数模型的方法可以有效避免模型的迟滞现象,提升标度因数模型的补偿效果,使光纤陀螺可以适应各种温度变化的环境。在-40℃~+60℃范围内同时对光纤环圈和光源的温度进行采集,并利用光源温度与平均波长的关系来修正标度因数模型,通过模型修正可以将光纤陀螺全温标度因数稳定性指标由常规模型下的36×10~(–6)提升到12×10~(–6)。     

光纤陀螺标度因数的测试误差分析

为了进一步提高光纤陀螺标度因数的测试精度,对光纤陀螺标度因数测试过程进行理论分析,确定了影响光纤陀螺标度因数测试误差的主要因素,并进行了计算机仿真和实验验证。结果表明:由于安装误差、北向地速分量以及转台速率精度的影响,光纤陀螺测试起始位置和采样时间的选择均会给小速率标度因数不对称性和非线性度的测试带来误差,而大速率标度因数的测试基本不受影响;通过对各输入速率点进行整圈采样,可以有效地降低小速率标度因数的测试误差,使其测试精度提高1个量级以上,实现对光纤陀螺标度因数性能更加准确的测试。     

双惯导联合旋转调制光纤陀螺标度因数误差自校正方法

旋转调制光纤陀螺航海惯导系统中,光纤陀螺标度因数误差会与地球自转角速度耦合产生等效的天向和北向陀螺漂移误差,也会与船体摇摆角速度以及惯性测量单元旋转调制角速度耦合产生短时动态误差,限制了长航时航海惯性导航精度。通过使用两套三轴旋转调制光纤陀螺航海惯导系统进行联合旋转调制,提出一种光纤陀螺标度因数误差在线估计与自校正方法。根据两套三轴旋转调制光纤陀螺航海惯导系统的水平旋转轴空间夹角关系建立观测方程,实现在线估计滤波。半实物仿真结果表明,自主导航过程中光纤陀螺标度因数误差在线估计精度优于1 ppm,利用输出校正方式在线补偿光纤陀螺标度因数误差导致的惯导定位误差,有效抑制了两套三轴旋转调制光纤陀螺航海惯导系统定位误差的增长。实际转台模拟实验中,两套三轴旋转调制光纤陀螺惯导系统300 h纯惯性导航整体定位最大误差分别减小25%和40%。算法采用地心地固坐标系,因此也适用于极区导航情况。     

光纤陀螺启动过程标度因数补偿方法

针对光纤陀螺启动过程标度因数变化大、稳定时间长的问题,提出标度因数补偿方案。分析了光纤陀螺启动过程中标度因数误差及超辐射发光二极管平均波长随温度变化误差产生的物理机制,建立了启动过程中标度因数误差的数学模型。进一步提出了一种通过测量温控电桥电路THERMIN端电压实时补偿启动过程标度因数的方案。试验结果表明,启动过程中(2 s内)光纤陀螺标度因数误差峰峰值从约25 000×10~(-6)降低到小于300×10~(-6),大幅提高了启动过程标度因数性能,满足了武器系统的快速启动需求。     

基于角速率积分法的光纤陀螺标度因数不对称度测量方法

标度因数不对称度是评价光纤陀螺的一项重要指标,对其进行精确测量在高精度导航应用中具有重要意义。传统的测试方法受转台速率控制精度限制,很难精确测量小于1×10-6的不对称度。首次提出基于角速率积分的标度因数不对称度测量方法。该方法给定转台正反方向转动的角度,由固定在转台上的被测光纤陀螺进行角速率测量,并对输出值积分,从而得到标度因数不对称度。该方法基于转台位置控制,避免了转速不稳定及正反向转速不对称等因素造成的影响。还对可能引起测量误差的因素进行了分析。最后采用角速率积分法测得高精度光纤陀螺标度因数不对称度小于1×10-6。     

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

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

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

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