小波多尺度分析在重力仪数据处理中的应用研究

为了有效地消除重力异常畸变对重力仪测量数据的影响，得到更高精度的重力异常测量值，在分析小波多尺度分析原理的基础上，首次将其应用到重力仪数据处理中，并与自适应卡尔曼滤波进行了对比分析，以实际数据与处理后数据偏差平方和的平均数作为衡量两种数据处理方法性能的指标。理论分析和仿真实验表明：小波多尺度分析和自适应卡尔曼滤波都能在一定程度上抑制重力异常畸变对重力仪测量数据的影响，但小波多尺度分析的性能优于自适应卡尔曼滤波。     

基于自寻优小波降噪算法的海洋重力数据滤波

为抑制船测重力数据中的噪声，提升海洋重力场数据精度，提出了可自动修正阈值参数的自寻优小波降噪算法。首先比较了IIR滤波器、FIR滤波器、Kalman滤波方法、传统小波降噪方法的滤波效果，分析了几种方法的特点及适用性。其次，根据海洋重力信号大多属于缓变的特点，借鉴遗传算法思想，设计了自寻优小波降噪算法。最后，利用湖上船载原子重力仪测量数据进行了实验验证，实验结果表明，自寻优小波降噪算法解决了传统小波降噪方法泛化能力弱、对复杂噪声分量滤除不彻底的问题，与IIR滤波器、FIR滤波器、Kalman滤波方法和传统小波降噪方法等滤波算法相比，滤波结果更接近有用信号，不需要对数据进行截短，滤波后信号精度提升了17%以上。     

分形滤波在高精度海洋重力仪数据处理中的应用研究

为有效地消除各种外界干扰噪声对高精度海洋重力仪测量值的影响,提高重力异常测量值的精度,在分析基于分形理论的滤波算法基础上,首次将其应用到高精度海洋重力仪系统数据处理中,并与自适应卡尔曼滤波进行对比分析,以实际信号数据与处理后信号数据的均方差作为衡量两种数据处理方法好坏的依据。理论分析和仿真实验表明:分形滤波方法和自适应卡尔曼滤波都能在一定程度上消除干扰噪声对重力异常信号的影响,但在相同背景条件下,分形滤波的性能优于自适应卡尔曼滤波。     

基于多尺度几何变换的导航雷达图象去噪仿真

导航雷达在采集、传输和显示过程中,由于多种因素的影响导致最终形成的图像中含有大量的噪声,影响了使用者对导航信息的分析和应用.传统的雷达图像去噪算法大多采用小波变换,但这种方法存在边缘模糊等问题.为了去除导航图像的噪声并解决小波变换中存在的边缘模糊问题,本文提出用基于多尺度几何变换的图像去噪方法对导航雷达图像进行处理,并利用基于多尺度几何变换的方法(包括基于Curvelet系数维纳滤波去噪方法和基于Contourlet域去噪方法)和基于小波变换的BayesShrink方法分别对含有模拟杂波和噪声的导航雷达图像进行仿真实验.实验结果表明:与基于小波变换的图像去噪方法相比,基于多尺度几何变换去噪方法能够更加有效去除雷达杂波和噪声.     

MEMS陀螺随机噪声的多尺度时间序列建模

基于小波多尺度分析方法,使用db4小波,将MEMS陀螺仪输出数据进行深度为4的多尺度分解,在各尺度上进行信号重建,对重建后的各尺度信号进行时间序列建模,各尺度的时间序列模型的输出和作为陀螺仪的随机噪声估计.多尺度时间序列建模方法提高了建模精度,将预测误差的方差降低了一个数量级,该方法已经成功应用到某型微小型飞行器的微导航系统中.     

流体饱和多孔隙介质波动方程多尺度反演

基于多尺度的思想,将小波多分辨分析和多尺度方法相结合,应用于流体饱和多孔隙介质孔隙率的反演。利用小波变换,将原始反问题分解为不同尺度上的一系列子反问题,并按照尺度从粗到细的顺序依次求解。在每一个尺度上,都采用稳定、收敛快的正则化高斯牛顿法求解,次一级尺度上求出的“全局最优解”作为上一级的初始解,依此类推,直到求出原始问题的真正的全局最优解。通过与传统的正则化高斯牛顿法相比较,显示了小波多尺度法是一个大范围收敛、能够有效节省计算量的方法,数值模拟的结果也表明了方法的有效性。     

海洋重力测量误差补偿技术

海洋重力场信息在勘探矿源和导航定位等方面都具有重要意义.进行海洋重力实时测量时,重力仪会受到各种外界扰动力的影响,再加上重力敏感器本身稳定性和惯性平台系统性能影响,重力敏感器的输出需要进行一系列数据处理和补偿后才能得到当地重力异常值,研究了重力敏感器安装角误差标定、零位漂移估计和格值修正等重力数据预处理方法.分析了海洋重力测量数据处理流程,主要包括零点漂移补偿、水平加速度误差补偿、厄特弗斯效应修正、高度修正和噪声滤波处理等.对每个数据处理过程都提出了具体补偿算法,并分析了补偿后的重力测量误差,将以上重力数据处理方法应用到实际重力测量,结果表明重力仪能够准确测量出当地重力值,其精度为1 mGa1.     

区间五次Hermite样条多小波Euler-Bernoulli梁单元

首先采用区间五次Hermite样条函数,分别构造了三节点梁的边界和中间节点的多小波尺度函数,然后,基于小波多辨分析思想,构建了梁单元位移多尺度近似空间的基函数系;最后,采用最小势能原理,得到弯曲梁的平衡方程,从而构造了区间五次Hermite样条多小波Euler-Bernoulli梁单元。算例结果表明,该小波单元可通过改变尺度来重新划分网格,从而可自由调节单个小波单元的计算精度,其计算精度与在相同网格划分下采用传统三节点Hermite梁单元计算的完全一致;与其它小波单元相比较,该小波单元具有计算简单明了,物理意义明确,易于理解的特点。     

自适应Kalman和零相移滤波算法在重力信号处理中的对比

为了有效消除海洋重力仪测量信号的噪声,提高重力数据的获取精度,根据随机过程理论,借助基于二阶高斯一马尔可夫异常位模型的重力异常协方差函数,得到海洋重力测量中重力异常信号的状态方程.对sage-husa滤波算法和零相移算法进行了理论分析,为了抑制零相移滤波器的首尾数据畸变,求解了滤波器的初始状态.根据实测重力数据,进行了去噪仿真试验.理论分析和仿真结果表明,sage-husa与零相移滤波算法均能较好地抑制采样重力数据中的噪声干扰,但零相移滤波算法的性能优于sage-husa滤波器.     

关于陀螺信号处理中小波基选取的研究

在建立陀螺漂移的数学模型的基础上，使用不同的小波基对陀螺输出信号进行多尺度小波分解，去除由噪声产生的小波系数，从而达到消除噪声的目的。从不同小波基的滤波效果来看，双正交样条小波Biorl.5具有较好的去噪效果，从而也验证了小波对称性在陀螺信号去噪中的重要性。     

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

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

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

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