考虑闭合效应和三维应力约束的表面裂纹扩展模拟

假定承受Ⅰ型常幅载荷下的表面裂纹在扩展中的形状保持为半椭圆，利用Newman半椭圆表面裂纹应力强度因子公式计算应力强度因子。提出了等效厚度的概念，利用穿透直裂纹的研究结果，考虑表面裂纹扩展中塑性致闭和三维应力约束效应。基于Elber模型建立了三维表面裂纹扩展模型。数值模拟了表面裂纹扩展过程，研究了裂纹形状变化及规律，计算了裂纹扩展寿命。将计算结果与有关试验结果进行了对比，吻合较好。     

表面裂纹疲劳扩展的数值模拟

建立了一种无形状约束的模拟表面裂纹在线弹性断裂力学条件下疲劳扩展的数值方法,并研究了表面疲劳裂纹形状演化和裂纹尖端应力强度因子(SIF)的分布特征。该方法以三维有限单元技术和Paris疲劳裂纹扩展规律为基础,并在裂纹扩展增量计算中考虑了裂纹闭合影响。本文第一部分主要介绍模拟三维疲劳裂纹扩展的数值方法的理论背景和相关的技术细节。着重分析和讨论基于三维有限单元法计算裂纹SIF所涉及的几个主要问题:裂纹尖端单元网格密度对估算精度的影响;自由表面的影响及其修正方法;裂纹尖端非正交单元网格的影响及修正方法。     

疲劳裂纹扩展过程中△K理论值的修正

通过聚碳酸酯紧凑拉伸试件的疲劳裂纹扩展速率试验和同时进行的光弹性实验,作者发现随着疲劳裂纹的扩展在裂纹尖端附近产生逐渐加强和扩大的残余应力场,揭示出了裂纹尖端附近的闭合效应的存在,它将使应力强度因子降低.因此,疲劳裂纹扩展过程中应力强度因子幅度△K 的理论值应该修正。本文提出了一种简便有效的修正方案.     

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

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

用三维光弹法确定轴表面横向半椭圆裂纹应力强度因子及其修正系数

用三维光弹法测定了轴表面横向半椭圆裂纹的应力强度因子，并分析计算了应力强度因子修正系数，结果表明，应力强度因子的分布型式对加载条件和裂纹几何不敏感；修正系数对裂纹深度呈现对称分布趋势，仅取决于裂纹几何参数。     

基于裂纹闭合效应的高载迟滞预测模型修正

裂纹闭合效应通常是导致I 型裂纹扩展在高载作用下发生迟滞效应的主要因素之一．本文采用汽车薄板QSTE340TM 材料，针对不同应力比，高载比条件下疲劳裂纹扩展行为进行了实验研究．论文通过断面分析，针对各参数对裂纹闭合效应的具体影响进行了分析讨论，认为裂纹作用区域随裂纹扩展而动态变化，从而提出了一种对有效应力强度因子幅的修正方法．通过在原有模型中引入幂函数形式的动态变量α，表征裂纹闭合效应的作用比例随裂纹长度的动态变化，取得了较好的效果．     

估算裂纹应力强度因子的新方法

根据裂纹形状与裂纹尖端应力强度因子分布之间的固有关系，在线弹性断 裂力学条件下，提出了一种按已知I型裂纹应力强度因子分布规律求裂纹形状及相应应力强 度因子的无梯度迭代法. 通过有限厚度、有限宽度板穿透裂纹和表面裂纹的数值模拟实例验 证了所提出方法的有效性和实用性，并对不同应力强度因子分布规律对裂纹形状以及相应的 应力强度因子大小的影响进行了分析和讨论. 所提出的方法有助于提高实际扩展裂纹应 力强度因子的估算精度以及更合理地预测疲劳裂纹形状演化.     

滚动接触下的轨道表面三维裂纹研究方法

针对已有的求解三维裂纹的应力强度因子方法工作量大、难以实现复杂的三维模型建立的问题,提出了一种求解广义移动荷载下三维裂纹尖端应力强度因子的普遍方法。采用包络式三维实体建模方法,有效避免了三维裂纹有限元模型网格划分时因受限于拓扑误差容忍而出现网格划分失败的现象;通过MATLAB三维插值施加节点荷载,解决了不规则三维裂纹有限元模型施加复杂空间荷载的难题。计算结果表明:荷载中心距离裂纹中心约4mm时裂纹尖端应力强度因子达到最大值7.41×107Pa·m1/2;轮轨接触疲劳以张开型为主。此外,给出了最大应力强度因子条件下裂纹尖端塑性区的分布及裂纹扩展情况,可为轮轨滚动接触疲劳的研究及设计提供一定的参考。     

三维斜置半圆形表面裂纹扩展数值模拟

以裂纹前沿局部坐标中主法平面内的最大周向拉应力作为控制参数,对第一主应力拉应力裂纹扩展准则进行了补充和修正.修正后的拉应力裂纹扩展准则称为最大周向拉应力裂纹扩展准则.它不仅可以有效地确定裂纹扩展方位、扩展方向和扩展条件,而且可以有效地确定裂纹扩展步长.根据修正的裂纹扩展准则,给出了三维裂纹曲面的数学描述方式,可以有效地模拟扩展曲面的弯折和扭曲.     

基于子模型法的带有表面裂纹钢丝应力强度因子研究

钢丝裂纹应力强度因子是进行钢丝疲劳断裂寿命评估、疲劳裂纹扩展分析和钢丝断裂强度评估等工作的重要参数。本文首先介绍了裂纹扩展分析软件FRANC3D,然后基于子模型法模拟研究了拉伸荷载作用下带有表面裂纹钢丝的应力强度因子,裂纹种类包括直线形裂纹和半圆形裂纹,最后拟合得到拉伸荷载作用下带表面裂纹钢丝的应力强度因子形状修正系数表达式,分析了利用该公式进行承载力评估时产生误差的原因。研究结果表明,利用子模型模拟分析拉伸荷载作用下带有表面裂纹的钢丝应力强度因子时计算精度高,计算速度快,对计算机硬件要求低;利用该方法得到的钢丝裂纹应力强度因子,在进行索承式桥梁吊索安全性能评估时,评估结果更精确。     

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

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

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

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