基于鞍点估计及其改进法的可靠性灵敏度分析

鞍点估计可以直接逼近非正态变量空间中线性功能函数概率分布, 进而得出 功能函数的失效概率. 在此基础上进行了基于鞍点估计的可靠性灵敏度分析. 对于非线 性功能函数, 尤其是强非线性功能函数, 基于鞍点估计进行可靠性及灵敏度分析时存在较大的误 差, 为此建立了基于鞍点估计的改进方法------鞍点线抽样方法的可靠性灵敏度分析. 在标 准化的变量空间中利用线抽样方法的样本点将系统失效概率转化为一系列线性响应功能函数 失效概率的平均值, 从而可靠性灵敏度转化为一系列线性响应功能函数的失效概率对随机变 量分布参数偏导数的平均值, 再采用鞍点概率估计方法直接估计非正态变量标准化空间中这 一系列线性响应功能函数的失效概率及可靠性灵敏度. 通过比较两种方法的基本思想、实现 过程和算例结果可以发现: (1) 第1种方法只适用于线性程度较好的功能函数的情 况, 其误差主要来源于非线性极限状态函数的线性化; (2) 改进方法给出的是失效 概率及失效概率对随机变量分布参数偏导数的估计值, 这些估计值随样本点数的增 加而趋于真值, 并且该方法可以考虑功能函数的非线性对失效概率的影响, 因 此具有广泛的适用范围.     

非正态变量可靠性分析的鞍点线抽样方法

结合鞍点概率分布估计和传统线抽样方法的优点,提出了非正态变量可靠性分析的鞍点线抽样方法.传统的线抽样方法对非正态变量可靠性问题进行分析时需将非正态变量等价转换为标准正态变量,这种非线性转换将增加响应功能函数的非线性程度,进而加大了转换后响应函数失效概率估计的难度.所提鞍点线抽样方法则无需将非正态变量转化为标准正态变量,它利用鞍点概率分布估计方法可以直接估计非正态变量空间中线性响应函数概率分布的特点,并利用线抽样方法可以将非线性功能函数的失效概率转化为一系列线性功能函数失效概率平均值进行估计的优点,实现了非正态变量空间非线性功能函数失效概率的高精度估计.鞍点线抽样方法使用前需将变量进行标准化变换,这种变换是线性的,通过对变量的标准化变换可以消除变量的量纲,从而使得标准化变量空间概率分布更具规律性.理论推导可以证明:鞍点线抽样方法在基本变量服从正态分布时将退化为传统的线抽样方法.算例验证结果表明:针对非线性功能函数的可靠性问题,鞍点线抽样方法比传统的直接鞍点估计具有更高的精度,比直接Monte Carlo模拟有更高的效率.     

正态变量相关情况下可靠性灵敏度分析的新方法

基于独立正态变量情况下可靠性灵敏度分析的线抽样法,提出了一种求解正态相关变量情况下可靠性灵敏度的新方法。在所提方法中,首先将正态相关变量等效变换为正态独立变量,然后利用线抽样方法独立完成等效独立变量情况下失效概率对独立变量的所有分布参数的灵敏度分析,最后依据等效变换前后变量分布参数之间的解析关系和复合函数求导公式,求得...     

基本变量对失效概率重要性测度分析的新方法

在研究矩独立基本变量对响应分布影响的重要性测度的基础上,定义了基本变量对失效概率的重要性测度.基本变量对失效概率的重要性测度可以直接全面地给出基本变量对结构安全影响的重要程度,因而该重要性测度将更具有工程指导作用.文中分析了所定义基本变量对失效概率影响重要性测度的基本性质,并基于鞍点线抽样在求解可靠度时不受随机变量分布形式的限制及其效率和精度较高的优点,提出了求解该重要性测度的鞍点线抽样方法.     

基于子集模拟和重要抽样的可靠性灵敏度分析方法

针对工程实际中大量存在的小失效概率问题，提出了基于子集模拟和重要抽样的可靠性灵敏 度分析方法. 在子集模拟重要抽样可靠性分析方法中，通过引入合理的中间失效事件，将小 的失效概率表达为一系列较大的条件失效概率的乘积，而较大的条件失效概率则可通过构造 中间失效事件的重要抽样密度函数来高效求解. 基于子集模拟重要抽样可靠性分析的思想， 论文将可靠性灵敏度转化为条件失效概率对基本变量分布参数的偏导数形式，推导了基于子 集模拟和重要抽样的可靠性灵敏度估计值及估计值方差的计算公式，并采用算例对所提方法 进行了验证. 算例结果表明所提方法具有较高的计算精度和效率，并且适用单个和多个失效 模式系统.     

????????Ч?????????????????·???

在研究矩独立基本变量对响应分布影响的重要性测度的基础上，定义了基本变量对失效概率 的重要性测度. 基本变量对失效概率的重要性测度可以直接全面地给出基本变量对结构安全 影响的重要程度，因而该重要性测度将更具有工程指导作用. 文中分析了所定义基本变量对 失效概率影响重要性测度的基本性质，并基于鞍点线抽样在求解可靠度时不受随机变量分布 形式的限制及其效率和精度较高的优点，提出了求解该重要性测度的鞍点线抽样方法.     

可靠性灵敏度函数及其特征指标的条件概率模拟求解方法

可靠性分析中基本变量分布参数为区间均匀变量时,失效概率为分布参数的函数。基于条件概率马尔科夫链模拟,提出了一种可靠性灵敏度函数的求解方法,并提出了一种新的可靠性灵敏度度量指标,它为参数可靠性灵敏度函数在参数空间上的期望。文中推导了线性极限状态正态变量下全局灵敏度函数及新指标的计算式.并提出了高效的基于条件概率马尔科夫链...     

结构可靠性灵敏度分析的方向（重要）抽样法

方向抽样法是在标准正态空间极坐标系下,通过对矢径的方向进行随机抽样来分析结构可靠度的.但是当极限状态面接近平面时,方向抽样法就没有优势了.为了提高方向抽样法的效率,提出了三种基于方向(重要)抽样法的可靠性灵敏度分析方法.根据独立标准正态空间中基本变量的x<'2>分布特性及矢径与随机变量分布参数的关系,推导失效概率对基本随机变量分布参数的可靠性灵敏度分析的计算公式.该文所提的可靠性及灵敏度计算方法有较高的计算效率和精度,对于高度非线性极限状态方程问题亦有很强的适应性.     

含模糊随机变量的多模式系统广义失效概率计算的内插迭代线抽样方法

通过正则化基本变量的度量空间,定义了单个基本变量同时具有模糊和随机双重不确定性时的广义失效概率.在广义失效概率的计算中,模糊随机变量被等价变换为随机变量.从而使得广义失效概率的计算变换为随机失效概率的计算.当模糊随机变量的密度函数和隶属函数均为正态型时,推导了其等价概率密度函数的形式和参数.采用自适应线抽样方法对基本变量同时具有模糊和随机不确定性时的多模式广义失效概率进行了计算,并采用数值算例对自适应线抽样广义失效概率计算方法的效率和精度进行了验证.算例分析表明该方法的计算结果是合理的,并且由于自适应线抽样法具有较高的效率和精度,因而所提方法具有一定的工程意义.     

对称抛物型模糊变量下的可靠性及灵敏度分析

当模糊变量的隶属函数为正态型时,可以将模糊随机可靠性及可靠性灵敏度转化为随机可靠性及可靠性灵敏度,并利用复合函数求导法则求解模糊随机失效概率对正态型隶属函数分布参数的灵敏度.对于对称抛物型隶属函数,文中提出了"改进最大最小"法和"改进等面积"法两种近似等价正态化方法,从而将模糊随机可靠性问题转换为随机可靠性问题,并利用线抽样方法分析之.算例结果表明,由于"改进最大最小"法所得的等价正态型隶属函数能在函数图形尾部更好地近似对称抛物型隶属函数,因而"改进最大最小"法更适用于模糊变量的隶属函数为对称抛物型分布时模糊随机可靠性及可靠性灵敏度的近似计算.     

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

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

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

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

Preface: Current research progress on mechanics of high speed rail

The construction and operation of high-speed rail（HSR）grid within the past two decades in China,in terms of the scale,may be possibly comparable to any national-wide construction in the history of China,even the Great Wall.By counting railways with commercial train service at the speed of200 km/h,China has the world’s longest HSR network with over 19 369.8 km of track in service today and this number     

ACTA MECHANICA SINICA 2014 ANNUAL SUBJECT INDEX

正~~     

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