失效概率矩独立全局灵敏度指标的高效算法

基于失效概率的矩独立全局灵敏度指标能够有效地分析输入变量的不确定性对结构系统失效概率的影响程度. 然而,目前以抽样方式来计算该灵敏度指标的方法都不能最大程度地利用样本. 因此,研究了在准确计算该指标的基础上如何提高样本的利用率. 基于所证明的连续区间上的全方差公式,提出了基于空间分割及重要抽样法来高效计算该指标的方法,其仅需一组样本,且计算量与输入变量的维数无关. 该方法首先通过重要抽样密度抽取一组样本,使得抽取到的样本以较大的概率落入失效域从而加快计算的收敛速度,其次,通过重复利用这一组样本来计算出各个输入变量的基于失效概率的矩独立全局灵敏度指标,大大提高了样本的利用率和计算效率. 验证算例的计算结果,说明了所提方法在计算效率、计算精度、收敛性及稳健性方面都较已有同类方法高,具有更好的工程适用性.      

含模糊变量和区间变量的结构可靠性优化设计方法

提出了基于能度可靠性理论的含模糊变量和区间变量的结构优化设计方法.运用能度可靠性理论建立了含模糊变量和区间变量的结构混合可靠性模型,采用结构的最大失效可能度作为混合模型的可靠性度量,给出最大失效可能度的计算方法,并以混合可靠性模型的最大失效可能度为约束条件建立了混合结构优化设计模型.实例计算表明了本文方法的有效性.     

含弧形裂纹复合陶瓷强度的尺度效应

认为含弧形裂纹复合陶瓷由随机方向的三相胞元与有效介质构成,用细观力学的方法研究了复合陶瓷的损伤失效和强度。首先确定三相胞元的外载应变,再依据复合陶瓷在损伤过程中的细观应力场和广义热力学力,计算出三相胞元内基体和颗粒的损伤等效应力,当基体和颗粒的损伤等效应力分别等于两者的极限应力时,得到基体和颗粒的破坏应力。然后,根据混合型应力强度因子计算弧形裂纹扩展时的能量释放率,进而得到界面的破坏应力。最后综合考虑基体、颗粒和和界面损伤影响,获得含弧形裂纹复合陶瓷的宏观强度及其尺度效应。     

具确定弱区域正交各向异性薄板的弹塑性屈曲分析

基于弹塑性力学和损伤理论,建立了一个与应力球张量有关的具损伤正交各向异性材料的混合硬化屈服准则,该准则无量纲化后与各向同性材料的Mises准则同构,在此基础上,建立了正交各向异性材料的增量型和全量型弹塑性损伤本构方程,并以具确定弱区域正交各向异性矩形薄板为例,根据屈曲时的能量准则和全量理论,以等效塑性应变为内变量,对其弹塑性屈曲问题进行了分析,讨论了几何参数和弱区域对正交各向异性薄板弹塑性屈曲临界应力的影响.     

采用重要抽样法的结构动力可靠度计算

首次对比分析了结构动力可靠度计算的三种重要抽样法,并对部分方法进行了补充修正.单元失效域法补充了依据随机教决定抽样区间的产生方法,根据单元失效域的条件概率和权重系数给出重要抽样密度函教.方差放大系数法直接通过激励过程的特性给出重要抽样密度函数的具体表达式.功率谱法的重要抽样密度函数仅为激励幅值的函数,根据结构反应的功率谱密度增大激励幅值的方差,建议幅值样本值的联合概率密度函数可表示为幅值样本值分量的概率密度函数的连乘形式.结果表明:对于线性体系三种方法的计算效率均比Monte-Carlo法有显著提高,而单元失效域法的计算效率又比另两种方法高.     

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

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

动力损伤后的脆性岩石静力蠕变断裂模型研究

应力波动力扰动下脆性岩石的静力蠕变特性,对深部地下工程围岩变形的评价有重要的实践意义.动力载荷作用导致的局部细观裂纹损伤严重影响脆性岩石蠕变力学行为.基于细观裂纹扩展与应力关系模型、动力扰动损伤演化函数、静动力载荷演化路径函数与黏弹性本构模型,提出一种应力波动力扰动下脆性岩石蠕变断裂特性的宏细观力学模型.其中动力损伤通过控制岩石内部细观裂纹数量变化实现.模型描述了应力波动力扰动下岩石的应变时间演化曲线,解释了岩石动力扰动下蠕变失效特性.研究了不同应力波幅值及周期影响下的脆性岩石应变-时间关系曲线,并通过试验结果验证了模型的合理性.讨论了动力损伤变化形式,突变发生时刻,突变量的大小对岩石蠕变失效特性的影响.分析了应力波幅值、周期对岩石动态动力损伤效应以及蠕变失效特性的影响.主要研究结果:动力损伤的变化值越大,岩石蠕变失效发生时间越短.冲击载荷扰动期间,动力损伤发生的时刻及增加的形式,对动力扰动后的岩石应变及蠕变破坏时间影响很小.动力损伤变化量随应力波幅值增加、周期减小而加速增大.应力波幅值越大、周期越小,岩石发生蠕变失效时间越短.     

一种基于凸集-概率混合模型的结构可靠性分析法

本文建立一种凸集-概率混合模型下的结构可靠性分析方法。考虑椭球模型和区间模型两种不同类型的凸集模型,在标准空间内通过拉丁超立方抽样生成样本点,通过矩阵变换将其转换到凸集空间内,得到凸集变量的样本点;将凸集变量样本点代入极限状态方程,从而混合可靠性模型转变为纯概率可靠性模型;利用Laplace渐进积分法对每个极限状态方程进行失效概率求解,统计结果的最大值和最小值,得到失效概率的上、下界,研究了凸集变量的个数对结果的影响,通过失效概率的变异系数描述计算结果的稳定性;通过三个算例验证了计算结果的精度,并采用混合模型的蒙特卡洛法进行对比计算。研究表明：本文所提方法计算精度和效率均较高;凸集模型的类型会对结果产生较大影响;为使结果趋于稳定,椭球模型所需的凸集样本点个数多于区间模型;若凸集样本点数目相同,椭球模型的失效概率计算结果变异系数较小,稳定性高于区间模型。     

一种基于凸集-概率混合模型的结构可靠性分析法

本文建立一种凸集-概率混合模型下的结构可靠性分析方法。考虑椭球模型和区间模型两种不同类型的凸集模型，在标准空间内通过拉丁超立方抽样生成样本点，通过矩阵变换将其转换到凸集空间内，得到凸集变量的样本点；将凸集变量样本点代入极限状态方程，从而混合可靠性模型转变为纯概率可靠性模型；利用Laplace渐进积分法对每个极限状态方程进行失效概率求解，统计结果的最大值和最小值，得到失效概率的上、下界，研究了凸集变量的个数对结果的影响，通过失效概率的变异系数描述计算结果的稳定性；通过三个算例验证了计算结果的精度，并采用混合模型的蒙特卡洛法进行对比计算。研究表明：本文所提方法计算精度和效率均较高；凸集模型的类型会对结果产生较大影响；为使结果趋于稳定，椭球模型所需的凸集样本点个数多于区间模型；若凸集样本点数目相同，椭球模型的失效概率计算结果变异系数较小，稳定性高于区间模型。     

航空涡轮盘多轴蠕变-疲劳寿命预测及对比研究

针对几种经典和新发展的蠕变-疲劳寿命模型开展综述介绍,并建立预测航空涡轮盘在循环热-机蠕变-疲劳载荷谱下蠕变-疲劳行为的数值流程,对某型航空涡轮盘的蠕变-疲劳损伤和寿命进行预测和对比．结果表明：等效应变法与临界平面法得出的疲劳损伤差距较小,等效应变法由于数值计算简单,工程适用性更强．寿命-时间分数(TF)法由于无法考虑应力松弛效应,给出了最为保守的蠕变损伤预测,其对盘体应力三轴度引起的损伤不敏感;延性耗竭法(DE)法仅以蠕变应变率作为损伤因素,虽考虑多轴蠕变因子的影响,但是给出的蠕变损伤过小;修正应变能密度耗竭(MSEDE)法综合考虑蠕变应变与应力松弛,并且考虑多轴蠕变因子与弹性跟随效应的影响,结合疲劳损伤模型可以给出合理的蠕变、疲劳损伤比例,其预测结果更加合理．     

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

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

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

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