基于OTTOSEN准则的混凝土粘塑性力学模型

研究了混凝土的弹．粘塑性动态本构关系,提出一个应变率相关的非线性混凝土模型．模型改进Bicanic的混凝土塑性间断面变化规律和Ottosen的四参数屈服准则,建立冲击荷载下的混凝土本构方程,可以应用于地震和爆炸作用下混凝土材料响应的研究．根据混凝土动态实验结果对应变率和材料强度的关系提出合理假设,考虑以往被忽略的材料峰值应变和泊松比随应变率的变化．模型包含静水压力参数和动态塑性损伤因子,可以有效地反映混凝土的动态力学行为,为其动态问题的研究提供有益的思路和有效的工具．     

基于广义非线性统一强度理论的隧道塌落拱分析

基于广义非线性统一强度的剪切破坏理论,在平面应变状态下,根据极限上限定理和隧道顶部的围岩塌落机制,在耗散能中引入了包含曲线型破裂面方程的目标函数,建立了围岩中任意断面的隧道顶部围岩的塌落机制。基于塑性位势理论,考虑塑性应变增量与塑性势函数的应力梯度成正比,由此得出了在速度间断线上任意点内能的耗散率;再运用虚功原理,建立内能耗散率和外力做功相等的关系式,从而通过变分原理得出了隧道顶部塌落面的解析表达式,由该表达式计算出隧道顶部塌落面的形状。以圆形断面隧道为例,对影响塌落拱形状的隧道半径和中间主剪应力系数等相关参数进行了讨论分析。由此得出:在平面应变状态下,隧道顶部围岩塌落体的宽度和高度随着中间主剪应力系数的不断增大而减小;塌落面的高度随着广义非线性统一强度参数的不断增大而增大,而塌落面的宽度则不断减小;塌落宽度和塌落高度随隧道半径的增大而增大。     

一种新的岩石非线性黏弹塑性蠕变模型研究

为了克服传统元件组合模型不能描述岩石蠕变过程中非线性特征的缺陷，首先根据加速蠕变阶段的应变和应变率随蠕变时间急剧增大的特点，建立黏塑性应变与蠕变时间的指数函数关系并提出非线性黏塑性体．将该非线性黏塑性体与广义Burgers蠕变模型串联，建立可以描述岩石全蠕变过程的非线性黏弹塑性蠕变模型，根据叠加原理得到一维应力状态下的轴向蠕变方程．然后基于塑性力学理论指出岩石三维蠕变本构方程建立过程中的不足之处，并给出非线性黏弹塑性蠕变模型合理的三维蠕变方程．最后采用不同应力水平下砂岩轴向蠕变试验对模型合理性进行验证，结果表明：拟合曲线与试验曲线吻合度较高，所建蠕变模型能够很好地描述砂岩在不同应力水平下的蠕变变形规律，尤其对加速蠕变阶段的非线性特征描述效果很好，验证了模型的合理性．     

压缩载荷作用下岩石的细观损伤和塑性研究

建立岩石微裂纹扩展的细观力学模型,研究了岩石的细观损伤和塑性性质.压缩载荷下微裂纹尖端翼裂纹稳定扩展表征岩石的细观损伤,采用应变能密度准则求解复合型断裂的翼裂纹扩展长度,微裂隙统计的二参数Weibull函数模型反映绝对体积应变对微裂纹分布数目影响,进而用翼裂纹扩展所表征的应力释放体积和微裂纹数目来表示含有微裂隙的岩石损伤演化变量;宏观塑性屈服函数采用Voyiadjis等的等效塑性应变的硬化函数,反映了塑性内变量对硬化函数的影响;建立岩石模型的本构关系及其数值算法,并用回映隐式积分算法编制了模型的本构程序.分析弹塑性损伤模型的围压对岩石损伤的影响,并从围压和短微裂隙长度等因素分析模型的岩石的损伤和宏观塑性特性.     

功能梯度材料纯弯曲梁的弹塑性分析

假设功能梯度材料为一理想弹塑性材料,其弹性模量和屈服强度沿梁高度方向按照幂函数变化,在小变形及平截面假设下研究功能梯度材料纯弯曲梁的弹塑性性能.根据Mises屈服准则导出了纯弯曲梁的弹性极限弯矩的解析表达式,建立了梁在弹塑性状态时截面弯矩与截面弹、塑性区分布之间的关系式,给出了梁进入塑性极限状态时中性轴的位置以及塑性极限弯矩的解析计算公式.数值算例的结果表明,功能梯度材料梁的弹塑性性能与均匀材料梁不同,其屈服不一定首先产生于截面最大应力点,而可能有多种不同的屈服模态及相应的塑性扩展.弹性模量及屈服强度的梯度变化对功能梯度材料纯弯曲梁的中性轴位置、截面弹塑性应力分布以及塑性极限弯矩均有较大影响.研究结果可为功能梯度材料梁的弹塑性分析提供一定的参考.     

压力敏感性弹塑性材料平面应力条件下起始扩展裂纹尖端的渐近场

采用一个能描述材料压力敏感性和SD效应的抛物型屈服准则,运用正交流动法则给出了材料的本构方程.在平面应力条件下,得出了裂纹尖端弹塑性应力场的基本解.在不同受力状态下(表现在采用不同的塑性混合因子),给出裂纹尖端应力场解的构造.讨论了材料参数对裂纹尖端应力场分布的影响.在自相似假设下讨论了应变的奇异性.结果表明,对于纯Ⅰ...     

粉砂岩峰后破坏区应力脆性跌落的试验和本构方程研究

对粉砂岩试样进行常规三轴压缩试验,得到了其在不同围压下的应力-应变全程曲线.结果表明,岩样应力-应变曲线峰前(破坏前区)部分的力学性质较为稳定,其力学行为可用经典强度理论进行描述;应力-应变曲线峰后(破坏后区)部分,岩样处于非稳定状态,其力学行为难以用经典强度理论来描述,而脆塑性模型恰好可以描述试验岩样的应力非连续变化特征.岩石的应力脆性跌落是有条件发生的,应力脆性跌落系数是围压的函数,并给出应力脆性跌落系数与围压的关系.对完整岩样而言,使用弹性-线性软化-残余塑性三线性非理想脆塑性模型可以描述岩样的本构特征,而对含有节理岩样而言,采用双线性弹性-线性软化-残余塑性四线型非理想脆塑性模型与试验更为相符.     

模型复合材料弹塑性界面应力分析

由纤维增强弹塑性基体所产生的界面具有弹塑性力学行为。考虑到一般材料的塑性变形都遵循幂硬化规律,对模型复合材料的界面进行弹性和应变硬化状态下的变形规律及其应力分析。以纤维拔出试验为研究模型,将界面分成弹性区和塑性区。利用界面应力剪滞理论,分别建立弹性区和塑性区的界面力学基本方程。选择适当的位移函数满足基本方程及埋入纤维的边界条件,再按位移函数求出弹性区和塑性区的界面剪应力。推导出平均界面剪应力与纤维     

基于塑性体积应变的岩石损伤变形特性实验研究

为分析岩石塑性变形与损伤的关系,在定义岩石的初始损伤和临界损伤,提出塑性体积应变分析方法,从而以塑性体积应变为损伤变量,采用归一化方法建立岩石的损伤本构模型。采用递增循环加载实验确定岩石损伤本构模型中的弹性卸载模量和弹性应变比例系数两个参数。通过实验和理论分析得出:当荷载较小时,普通单轴压缩状态下岩石损伤随荷载的增加具有减小趋势,荷载超过一定数值后,岩石损伤才开始增加;单轴递增循环压缩状态下当循环荷载大于约35%峰值强度后,卸载后岩石的损伤具有增加的趋势,小于该荷载之前具有减小的趋势。整个加载过程的理论应力-应变曲线能很好地与实验结果相吻合,在循环加载区间理论结果还能体现出岩石实验结果中的回滞环。     

岩石材料单轴压缩损伤变量与纯剪切损伤变量间的关系

本文以岩石材料损伤应变能释放率为基础，利用不同的岩石屈服准则，探讨了岩石材料单轴压缩损伤变量与纯剪切损伤变量间的关系。借助于岩石材料单轴压缩或纯剪切中任一试验，即可求得两种受力状态下的损伤变量。     

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

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

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

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