基于1范数正则化的模型修正方法在结构损伤识别中的应用

以基于灵敏度分析的有限元模型修正方法为基础,提出了一种基于1范数正则化过程的结构损伤识别方法。通过与以Tikhonov正则化为代表的二次型正则化过程相比较,本文的理论分析表明1范数正则化方法在迭代计算过程中能根据上一迭代步损伤识别结果自适应地调整正则化项中的损伤参数权系数,从而显著改善了Tikhonov正则化识别结果过度光滑的缺陷,更利于识别结构的局部损伤。为解决引入1范数造成的数值计算困难,文中还对基于1范数正则化的模型修正算法进行了改进。以二维框架模型为例的损伤识别数值模拟表明：1范数正则化方法与模型修正方法相结合可以有效抑制实测模态参数中噪声的影响,体现出较好的鲁棒性;在模态噪声水平达到10%的情况下,仍能有效抑制噪声干扰,凸显结构局部损伤位置,准确识别损伤程度。     

基于模态参数灵敏度的损伤方程组求解正则化方法研究

基于模态参数的结构损伤识别方法是振动损伤识别领域中应用最为广泛的方法.利用模态参数灵敏度构建结构损伤方程组,对其进行求解可以识别结构损伤位置和程度.由于实际工程中模态参数不完备性和噪声的影响,结构损伤方程易出现病态问题,直接求解可能产生错误的结果.为了解决这一问题,可以引入正则化方法进行求解.然而,各类正则化方法的基本...     

面载荷识别的TSVD正则化方法

线性系统的响应可以表示为单位脉冲响应函数与动态载荷的卷积分,经离散后得到载荷识别反问题。引入TSVD正则化方法处理反问题的不适定性,对不同噪声的实测响应进行了面载荷的反演求解,并对识别精度做出了分析。算例表明,TSVD正则化方法对响应误差的适应性强,识别精度高。     

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

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

求解非线性反问题的大范围收敛梯度正则化算法

基于同伦映射的思想，改进了求解非线性反问题的梯度正则化算法。通过路径跟踪有效地拓宽了梯度正则化算法求解的收敛范围。对于正则化参数的修正，通过引入拟Sigmoid函数，提出了一种下降速率可调的连续化参数修正方法，在保证迭代稳定的条件下，得到较好的计算效率，同时保证该算法具有很好的抵抗观测噪声能力。实际算例表明，该方法收敛范围宽，计算效率高，在存在较强观测噪声的条件下也能得到很好的反演结果。     

基于L1正则化的随机梁式结构静力损伤识别方法

提出了一种结合摄动法和L1正则化方法的随机梁式结构静力损伤识别方法。考虑初始模型误差和测量误差的影响,建立了关于随机损伤指数的控制方程,并将摄动法和L1正则化方法相结合,对随机损伤指数的控制方程进行求解,进而从概率的角度对结构的损伤进行识别。损伤试验结果表明,和传统的最小二乘求解法相比,本文方法能够更为准确地识别多处局部损伤的位置及大小,对实际结构损伤检测具有较好的参考价值。     

一种用于重力梯度动态测量的载体环境引力梯度补偿方法

重力梯度仪动态测量时,重力梯度敏感器一直稳定在地理坐标系下,载体姿态变化使载体质量分布相对敏感器的位置发生变化,形成载体环境引力梯度变化。为提高重力梯度仪动态测量精度,提出一种基于Tikhonov正则化的载体环境引力梯度补偿方法。首先,推导了载体环境引力梯度的解析模型,建立了引力梯度变化的回归方程。然后,针对回归算子病态性问题,提出了Tikhonov正则化方法,通过半物理仿真确定最优正则化参数,使补偿量的误差控制在2%以内。最后,利用该参数处理船载试验实测数据,结果表明:所提出的方法对载体环境梯度变化补偿具有明显的效果,可将两路重力梯度测量信号内符合中误差分别降低19 E和21 E,补偿后重力梯度测量精度达到10 E的精度水平。     

接触力优化设计的熵正则化方法

包含单侧接触的机械和结构系统中,接触面间的局部接触力很大时,容易导致结构局部磨损和破坏。利用结构优化的手段极小化接触面上的最大接触力是调整接触面上接触力分布的一个有效途径。然而在数学上,目标函数的不可微和约束中的接触条件导致该结构优化问题属于非光滑和非凸优化问题,因而常用优化算法的优化结果是不可靠的。为此,本文提出接触力优化的一个熵正则化模型,并利用基于正则化过程的序列二次结构优化算法求解。所提的方法在数值上避免了处理不可微问题,而且正则化参数的调整还有助于获得更小的局部最优解。以离散结构为例的数值算例验证了本文的模型和算法。     

材料物性参数识别的梯度正则化方法

本文对梯度正则化方法（Ｇｒａｄｉｅｎｔ－Ｒｅｇｕｌａｒｉｚａｔｉｏｎ Ｍｅｔｈｏｄ）作了进一步的研究,给出一种建立了梯度正则化迭代算法和选择正参数的简明实用方法。文中椭圆算子方程参数识别算例不仅说明了ＧＲ法具有广泛的适应性和一定的抗噪声能力,而且收敛速度较快,具有较大的收敛范围。     

微分方程反问题的梯度正则化法

本文在微分方程反问题研究当前发展状况的基础上,提出了一种求解微分方程反问题的比较统一的方法——梯度正则化法。该方法将反问题的求解分为两个主要部分,一是展开以求补充条件对未知参数的梯度矩阵,二是正则化以解不适定的线性方程组。梯度正则化法从反问题的共性入手,未附加任何特殊的约束,所以可以适用各种类型的反问题,且在求解时不受空间维数的限制。     

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

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

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

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