脆性断裂的非局部力学理论

本文提出一种脆性材料断裂的非局部力学理论,内容包括:Ⅰ、Ⅱ、Ⅲ型Griffith裂纹的非局部弹性应力场,裂纹尖端邻域非局部弹性应力场的渐近形式,脆性开裂的最大拉应力准则。文中给出了这种理论应用于三种基本型裂纹和Ⅰ-Ⅱ、Ⅰ-Ⅲ复合型裂纹临界开裂条件的计算结果,并把它们与一些试验资料和最小应变能密度因子理论进行了对比。     

扭力轴三维裂纹扩展寿命仿真研究

对疲劳载荷谱作用下三维表面裂纹,采用双重边界元理论求解裂纹前沿的应力应变场,运用Forman理论、最小应变能密度法和Elber模型,计算裂纹前沿各点的扩展长度、扩展方向和应力强度因子等特征量.根据增量步下裂纹几何形状的改变,对裂纹面进行网格重划分和迭代计算,模拟三维裂纹的扩展,预测裂纹扩展寿命.扭力轴表面裂纹扩展的仿真结果表明该方法合理可行.     

一种新的确定应力强度因子的数值方法

应力强度因子是一个非常重要的参数,可以用来估算裂纹和切口的断裂.这篇论文提供了一种基于包含应力集中区域一定体积上的平均应变能密度,来确定应力强度因子的数值方法.对于I型或是II型裂纹的单一加载方式,应力强度因子都可以直接从一定体积上的平均应变能密度的表达式求得其解,但是对于I-II复合型裂纹,情况相对复杂.因此,作者们提出了利用围绕切口尖端一定体积上几组不同关于裂纹切口平分线对称区域上的平均应变能密度,来拟合复合加载下I型和II型应力强度因子的数值方法.为了验证,计算了I-II复合型裂纹的半圆形三点弯曲试样应力强度因子,并与文献中给出的应力强度因子进行了比较.结果表明,提出的数值方法可靠,为平均应变能密度准则的工程应用提供了一种新的思路.     

复合型裂纹扩展的应变能准则

现有的复合型裂纹扩展准则主要可分为两大类,一是以裂纹扩展前的应力-应变场为基础,二是基于能量平衡观点的能量释放率理论.例如,最大切向应力准则,最大主切向应力准则,最大切向应变准则和应变能密度准则都属于前者;Griffith-Irwin能量准则则属于后者.所有的这些准则几乎都假设裂纹扩展将沿着以裂纹尖端为中心的射线方向,并且在此方向上的某个基本函数存在着一个极值,当此极值达到临界值时,裂纹开始失稳扩展.这些准则的着眼点都是只考察裂纹尖端附近某一点处的某个力学参量(如应力、应变、应变能密度或能量释放率等)的变化.本文所考察的则是裂纹尖端周围某一特定区域内总应变能的变化.     

考虑夹杂相互作用的复合陶瓷夹杂界面的断裂分析

复合材料中夹杂含量较高时,夹杂间的相互作用能显著改变材料细观应力应变场分布,基体和夹杂中的平均应力应变水平也会发生较大变化,导致复合材料强度等力学性能发生显著变化. 为修正单一夹杂模型运用在实际材料中的误差,基于相互作用直推估计法,建立一种考虑含夹杂相互作用的夹杂界面裂纹开裂模型. 首先根据相互作用直推估计法,得到残余应力和外载应力共同作用下夹杂中的平均应力,再计算无限大基体中相同的夹杂达到相同应力场时的等效加载应力,将此加载应力作为含界面裂纹夹杂的等效应力边界条件,在此边界条件下求得界面裂纹尖端的应力强度因子,进而得到界面裂纹开裂的极限加载条件,并分析了夹杂弹性性能、含量、热残余应力、夹杂尺寸等因素对界面裂纹开裂条件的影响. 结果表明,方法能够有效修正单夹杂模型运用在实际材料中的误差,较大的残余应力对界面裂纹开裂有重要的影响,夹杂刚度的影响并非单调且比较复杂;在残余应力较小时,降低柔性夹杂刚度或者增大刚性夹杂刚度都有利于提高材料强度;扩大夹杂尺寸将导致裂纹开裂极限应力显著降低,从而降低材料强度.      

MY准则解析复合型裂尖塑性区

首次以MY(平均屈服)准则对I-II复合型裂纹在小范围屈服下的裂尖塑性区进行了分析,分别获得了平面应力和平面应变状态下塑性区尺寸的解析解。这两解表明,塑性区尺寸是材料屈服强度、应力强度因子、极角θ的函数。与Tresca准则、TSS屈服准则获得的解以及Mises解比较表明:Tresca准则预测塑性区上限,TSS屈服准则预测塑性区下限,MY准则预测的塑性区居于Tresca与TSS塑性区之间,逼近Mises解。另外,文中讨论了平面应力和平面应变状态下裂纹尖端的开裂问题,结果表明:当裂纹角β=π4时,平面应力状态下裂纹沿0-θ=0.2952π方向开裂;平面应变状态下裂纹沿0-θ=0.3188π方向开裂。     

紧凑拉剪试样中Ⅰ/Ⅱ复合型穿透裂纹的三维特性分析

进行了Ⅰ/Ⅱ复合型载荷作用下航空结构铝合金CTS试样线弹性的全场三维有限元计算,分析了复合型离面应力约束因子Tz和面内约束T应力的分布特性,研究了厚度和载荷条件对应力各分量及应力三轴性水平Rσ、应变能密度U0的影响以及这些量在实验中观察到的裂纹起裂方向上的特点。结果表明：（1）平面内约束T应力在此种试样形式下为零;（2）复合型裂纹的三维效应区与厚度成正比,为0.4-0.5B,与载荷条件基本无关,但是厚度效应随着载荷角的减小逐渐变小,到II型载荷时基本消失;（3）离面约束因子Tz随着径向和厚度尺寸的增加逐渐减小,但在周向基本不发生变化;（4）最小应变能密度U0min的方位能够表征此种材料三维复合型断裂的裂纹起裂方向。研究结果为建立三维复合型断裂准则提供了理论基础。     

复合V形切口脆断的形状改变能密度因子准则

沿袭裂纹尖端附近的形状改变能密度因子概念思路,定义V形切口尖端附近的形状改变能密度因子.基于V形切口尖端附近的奇异应力场,建立了Ⅰ-Ⅱ复合V形切口脆性断裂的形状改变能密度因子准则.应用该准则对复合载荷下双边切口试件进行了理论预测,并把预测结果与实验结果、应变能密度因子准则的预测结果进行了比较.结果表明,形状改变能密度因子准则预测的起裂角和断裂曲线与实验值更吻合.     

黏弹性界面断裂分析的增量“加料”有限元法

提出了一种适用于黏弹性界面裂纹问题的增量“加料” 有限元方法. 利用弹性界面裂纹尖端位移场的解答,通过对应原理和拉普拉斯逆变换近似方法,得到了黏弹性界面裂纹的尖端位移场. 用该位移场构造了黏弹性界面裂纹“加料” 单元和过渡单元位移模式,推导了增量“加料” 有限元方程,求解有限元方程可获得应力强度因子和应变能释放率等断裂参量. 建立了典型黏弹性界面裂纹平面问题“加料” 有限元模型,计算结果表明,对于弹性/黏弹性界面裂纹和黏弹性/黏弹性界面裂纹,该方法都能得到相当精确地断裂参量,并能很好地反映蠕变和松弛特性,可推广应用于黏弹性界面断裂问题的计算分析.      

黏弹性功能梯度材料裂纹问题的有限元方法

针对组分材料体积含量任意分布的黏弹性功能梯度材料裂纹问题建立有限元分析途径. 通过Laplace变换,将黏弹性问题转化到象空间中求解,基于反映材料非均匀的梯度单元和裂纹尖端奇异特性的奇异单元计算象空间中的位移、应力和应变场,应用虚拟裂纹闭合方法得到应变能释放率,分别由应力和应变能释放率确定应力强度因子. 给出这些断裂参量在物理空间和象空间之间的对应关系,由数值逆变换求出其在物理空间的相应值. 文中分析两端均匀受拉的黏弹性边裂纹板条,首先针对松弛模量表示为空间函数和时间函数乘积的特殊梯度材料进行计算,结合对应原理验证方法的有效性. 然后分析组分材料体积含量具有任意梯度分布的情形,由Mori-Tanaka方法预测象空间中的等效松弛模量. 计算结果表明,蠕变加载条件下,应变能释放率随时间增加,其增大程度与黏弹性组分材料体积含量相关. 由于梯度材料的非均匀黏弹性性质,产生应力重新分布,导致应力强度因子随时间变化,其变化范围与组分材料的体积含量分布方式有关.     

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

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

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

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