北斗三频精密姿态测量中整周模糊度的单历元可靠固定（英文）

载波相位模糊度固定是GNSS高精度姿态测量中的关键问题。针对附有基线长约束的姿态测量场景提出了一种可靠的Bei Dou三频模糊度单历元快速固定算法。首先使用无几何模型可靠地固定两个超宽巷模糊度;其次使用这两个超宽巷模糊度计算一个宽巷模糊度;然后在固定模糊度的宽巷观测值的辅助下,使用几何相关模型对第三个窄巷模糊度(该模糊度需与两个超宽巷模糊度线性独立)的浮点解进行估计,并使用C-LAMBDA方法对其进行固定;最后对三个频点上的原始模糊度进行恢复并利用其进行单历元姿态测量。多组Bei Dou三频实测数据的处理结果表明,提出的单历元模糊度固定方法具有固定成功率高和计算速度快的特点,能够有效地避免周跳的影响,非常适合高采样、动态复杂环境下的多频GNSS姿态测量数据处理应用。     

GPS快速定位中相位模糊度动态解算的一种正则化方法

针对GPS快速定位中少数历元组成的法方程存在严重病态性的问题,研究了GPS单频整周模糊度快速解算的新方法:即首先采用改进SVD分解获得系数矩阵的精确奇异值,避免了较小奇异值抖动的影响;然后基于改进SVD分解结合法矩阵病态性特点,设计了一种改进Tikhonov正则化方法,合理构造了正则化矩阵,有效抑制法矩阵的病态性。算例表明:与LS估计-LAMBDA方法和Tikhonov正则化-LAMBDA方法相比,新方法能够有效降低法矩阵的条件数约3个数量级,仅解算4个历元数据,浮点解偏离真值的方差从41.89减小到1.04,可以快速获得准确、稳定的模糊度浮点解。模糊度固定性能结果进一步表明,新方法显著地提高模糊度搜索效率和成功率,解算成功率提高100%。     

网络RTK参考站间模糊度动态解算中病态方程的一种解算方法

在网络RTK参考站间的模糊度估计中,若误差方程严重病态,将导致模糊度解与其准确值偏差较大或整周模糊度无法固定,因此提出了一种适于网络RTK模糊度动态解算的新方案：1）法方程病态性的判断;2）Tikhonov正则化解算病态方程;3）LAMBDA方法搜索固定整周模糊度。同时,深入研究了Tikhonov正则化矩阵的构造方法和正则化参数的选取准则。最后以实例验证了采用此方案解算病态方程是可行的,通过选取合适的正则化参数可以解得准确的整周模糊度;详细讨论了选择不同的正则化参数对模糊度解算结果的影响。     

三频差分GNSS/INS紧组合模型

北斗卫星导航系统已经提供三频信号的位置服务。针对目前尚无关于北斗三频信号的GNSS/INS紧组合研究的现状,将三频模糊度解算算法与紧组合结合,构建了三频差分GNSS/INS紧组合模型。给出了新模型的观测方程和模糊度解算方法,分析了新模型的优点。进行了车载组合导航实验对比新模型与经典模型,结果表明:新模型的解算时间消耗远少于经典模型,在9颗可视卫星时,新模型的总解算时间仅为经典模型的55.53%;卫星数的增加对新模型解算时间的影响不显著;新模型能够获得和经典模型相当的精度,位置差别均方根在毫米级,速度差别均方根在厘米级,姿态差别均方根在0.01°级别。新模型利用了三频模糊度解算速度快的优点,又有效地削弱了其不稳定性对定位结果的影响,为基于三频信号的GNSS/INS紧组合做了前沿性研究。     

基于BDS-3五频超宽巷/宽巷组合的中长基线单历元定位方法

为充分发挥BDS-3五频信号的定位优势,提出了一种基于BDS-3五频超宽巷/宽巷组合的中长基线单历元定位方法。首先构建了伪距/载波联合的无几何无电离层(GIF)分步解算模型,分析了各模糊度的先验解算精度,分析结果表明:五频组合中有三个线性无关的超宽巷可采用GIF模型单历元可靠固定。其次针对第四个宽巷采用GIF模型单历元固定精度低的问题,采用几何相关模型,对伪距、三个模糊度固定的超宽巷与第四个宽巷方程进行联立解算,并利用LAMBDA方法进行模糊度的搜索固定。理论分析结果表明,BDS-3固定四个超宽巷/宽巷模糊度后,站星距离精度可达0.206 m,相比三频BDS-2和GPS分别提升60.08%和57.35%。实测数据结果表明,BDS-3超宽巷/宽巷模糊度基本可实现单历元可靠固定,模糊度固定后在北/东/天三个方向上可分别实现0.161 m、0.187 m和0.457 m的定位精度。     

单频GPS/BDS/MEMS IMU紧组合模糊度固定抗差模型

针对高精度GNSS/INS组合定位中,伪距多路径误差严重影响模糊度固定效率及系统定位精度的问题,引入抗差估计函数,建立INS辅助GPS/BDS模糊度快速固定抗差模型。首先建立了单频GPS/BDS/INS紧组合动力学模型以及观测模型,继而研究了紧组合系统INS辅助下GNSS模糊度快速分解模型,分析了伪距观测值粗差对于模糊度参数估值的影响。为减弱异常观测值对于模糊度固定成功率的影响,提出了附有INS定位约束的模糊度解算抗差算法。最后利用实测导航定位数据验证了算法的有效性。结果表明:对比不考虑观测值粗差影响的模糊度固定算法,模糊度固定抗差算法显著提高了模糊度固定的可靠性,抗差算法的模糊度固定效率不受粗差值大小的影响;对于模拟的伪距多粗差情形,抗差算法对GPS/INS、BDS/INS和GPS/BDS/INS三种组合方案的模糊度固定Ratio均值分别提高了100.8%、47.7%以及19.5%;在INS定位松约束条件下,抗差算法提高了模糊度固定成功率。     

变基线约束的GNSS单频单历元姿态测量新算法

针对基于值域的GNSS姿态测量算法没有考虑工程应用中基线小动态变化的问题,推导了模糊度反约束值域搜索算法的搜索误差模型,研究了该误差模型对算法成功率的影响,分析了利用基线残差最小固定模糊度在工程应用中的不可靠性,提出了变基线约束的单频单历元姿态测量新算法:以基线长度小动态变化为约束,以惯性导航设备(INS)提供的粗略方位角、俯仰角为辅助,通过对基线长度、俯仰角、方位角三维搜索确定模糊度搜索域,以基于模糊度残余最小的无参化目标函数固定模糊度。试验结果表明,新算法能够在基线动态变化±10 cm、方位角偏差±5°、俯仰角偏差±5°情况下,实现模糊度的单频单历元固定和载体高精度姿态测量,并具有较高的成功率。     

GPS载波相位双差整周模糊度在航求解

鉴于整周模糊度的快速解算是GPS载波相位测姿的一个主要技术难题,从系统实时性、动态性要求出发,给出了一个整周模糊度的快速解算的新方法。其主要特点为:利用三角函数约束条件减少模糊度搜索范围,优选指标不受基线运动状态影响,可利用多历元动态基线数据进行平差。实验结果验证了该算法的有效性。     

一种改进的LLL模糊度规约算法

整周模糊度的高效解算是GNSS高精度数据处理中的关键,基于格论进行GNSS模糊度估计时需要通过格基规约来实现最优整周模糊度向量的快速搜索。针对高维情况下常规LLL规约算法辅助整周模糊度解算存在规约耗时较长和规约性能有限的问题,引入最小列旋转QR分解技术对基向量进行预排序,采用延后尺度规约和部分尺度规约来减少规约过程中的冗余尺度规约,以改善LLL算法的执行效果。分别通过模拟和实测数据进行实验,结果表明:改进后的LLL算法可以明显降低格基规约耗时,实测环境下其规约效率相比于传统方法提高了约10倍,且能够保证较好的规约性能,从而有效提升高维模糊度的解算效率。     

基于伪距宽巷组合的GPS模糊度固定方法

中长基线模糊度快速解算是GPS网络差分技术的核心内容。根据网络差分模糊度固定只用于基线解算而不用于定位的特点,通过对常用GPS原始观测数据组合方式的分析,提出伪距宽相组合的数据处理方法,有效消除了电离层、对流层传播误差,形成抗差性强、大气误差自由的GPS组合观测值;在此基础上引入双差伪距宽相组合进行模糊度浮点解并建立法方程,应用高度角与大气误差的关系构造权阵,对常规LAMBDA算法进行了改进,形成一种适合GPS网络差分的中长基线模糊度解算方法。对三个参考站同步观测数据的实际测试结果表明：使用该方法网络模糊度解算时间小于300s,基线长超过60km,并满足闭合性原则。     

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

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

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

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