有限厚度板穿透裂纹前缘附近三维弹性应力场分析

通过三维有限元计算来研究有限宽度、有限厚度含有穿透裂纹板的裂纹前缘应力场，从中找出应力强度因子与板的厚度、裂纹长度之间的关系，同时还分析了裂尖的三维约束程度和三维约束区的大小。分析结果表明：应力强度因子沿厚度的分布是不均匀的，应力强度因子的最大值及其位置与厚度有关；有限厚度板中面应力强度因子(KI)m-p及最大应力强度因子(KI)max均大于平面应力或平面应变的应力强度因子。对有限厚度裂纹问题，按平面应力或平面应变来考虑是不安全的；板中面的应力强度因子(KI)m-p及最大应力强度因子(KI)max是厚度B／a的函数；板的中面离面约束系数Tx最大，自由面(z=B)Tx=0。沿厚度方向裂尖附近的离面约束系数Tx也是z／B和B／a的函数，随着厚度的增加离面约束系数Tx增大，离中面越近离面约束系数Tx越大。Tx随着x的增大急剧减小，三维约束影响区域大小大约为板厚的一半，且裂纹长度a／W对应力强度因子沿厚度变化规律及Tx影响区域大小影响较小。     

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

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

复合材料断裂力学参量的实验测定

本文采用三维云纹干涉仪同步测量含裂纹复合材料弯曲板的面内位移场及离面位移场，并对离面位移场云纹图进行计算机自动采集和处理，编制了云纹数字图像处理软件。文中推导了正交异性复合材料弯曲问题的离面位移场表达式，得到了弯曲板裂纹尖端位移w与应力场强度因子K1的关系式。通过实测位移，结合裂尖位移场推导出应力强度因子K1值。     

双轴载荷作用下源于椭圆孔的分支裂纹的一种边界元分析

利用一种边界元方法来研究双轴载荷作用下无限大板中源于椭圆孔的分支裂纹．该边界元方法由Crouch与Starfied建立的常位移不连续单元和笔者提出的裂尖位移不连续单元构成．在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界,文中算例说明本数值方法对计算平面弹性裂纹的应力强度因子是非常有效的。该文对双轴载荷作用下无限大板中源于椭圆孔的分支裂纹的数值结果进一步证实本数值方法对计算复杂裂纹的应力强度因子的有效性,同时该数值结果可以揭示双轴载荷及裂纹体几何对应力强度因子的影响。     

I型裂纹尖端约束应力区模型及其解析解

考虑了I型裂纹尖端损伤区域内三种不同的约束应力分布形式,即右三角分布形式(情况A)、均匀分布形式(情况B)、左三角分布形式(情况C),并采用复变函数方法求得了应力强度因子与裂纹张开位移的解析解;在此基础上,通过数值计算得到了应力强度因子和裂纹张开位移随约束应力区长度、约束应力大小以及分布形式的变化规律。研究结果表明:随裂尖材料损伤程度的增加,裂尖损伤区内约束应力减小,应力强度因子和裂纹张开位移增大;约束应力的分布形式对应力强度因子和裂纹张开位移有显著影响;相对于其他区域,约束应力对裂纹尖端区域裂纹张开位移的影响较大。然而,对于裂尖损伤区域的形成与作用荷载、材料性质、构件几何尺寸之间的关系,还需要进行更为深入的研究。     

三维约束下幂硬化双材料界面裂尖场分析

本文用Ｊ２形变理论分析了三轴应力约束下，由具有不同硬化指数和硬化系数的幂硬化材料构成的界面裂纹问题．通过在本构方程中引入三轴应力约束因子Ｔｚ，我们考虑了三维约束效应对裂尖场的影响．计算结果表明，三轴应力约束对裂尖场的分布，尤其是对径向应力σｒ和有效应力σｅ有重要影响．本文还探讨了三轴约束和材料性能对裂尖应力奇异性的影响     

基于有限元方法的航空发动机叶片应力强度因子计算

为研究叶片裂纹尖端的应力奇异性,以某型航空发动机压气机叶片为例,利用有限元方法研究了叶片裂纹尖端应力强度因子的计算方法,并研究了旋转叶片振动状态下裂尖应力强度因子随裂纹长度的变化规律。建立计算模型时,在裂纹尖端划分了三维奇异单元,在裂尖外围划分了过渡单元。计算结果表明：研究旋转叶片振动状态下的裂尖应力奇异性,仅利用I型应力强度因子就具有足够的精度;对于同一裂纹,绝大多数情况下叶盆面应力强度因子大于叶背面应力强度因子,故研究叶片应力强度因子时只需研究叶盆应力强度因子即可;随着裂纹扩展,叶盆面I型应力强度因子不断增大。本文的研究方法及结论为进一步研究叶片的裂纹扩展规律及损伤容限奠定了基础。     

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

进行了Ⅰ/Ⅱ复合型载荷作用下航空结构铝合金CTS试样线弹性的全场三维有限元计算,分析了复合型离面应力约束因子Tz和面内约束T应力的分布特性,研究了厚度和载荷条件对应力各分量及应力三轴性水平Rσ、应变能密度U0的影响以及这些量在实验中观察到的裂纹起裂方向上的特点。结果表明：（1）平面内约束T应力在此种试样形式下为零;（2...     

带表面裂纹LY12板材弹塑性抗裂特性的研究

本文对不同的热处理状态(CZ、R)、厚度(B=6、8、10和14mm)和裂纹深度与板厚比(α_√B=0.45～0.55和0.20～0.30)的带表面裂纹LY12铝合金有限宽板进行了塑性抗裂特性的试验研究。提出了等效裂纹扩展假设,测得了起裂时表面裂纹咀的张开角口α1,它是常数;导出了裂尖张开角αT的近似计算公式,并测得裂纹扩展时αT的变化规律;测得了稳定扩展和失稳扩展过程中裂纹塑性张开角(dδ_(MP)/da)_S和(dδ_(MP)/da)_U,它们是不同的常数;还测得了不同的热处理状态、厚度和α_0/B带表面裂纹LY12板材从起裂到失稳的裂纹等效扩展量,它约为板厚的4%。     

圆柱形界面相裂纹在扭转载荷作用下的动态断裂

本文研究了位于界面相中的圆柱形界面裂纹的扭转冲击问题.采用Laplace、Fourier变换和位错密度函数将混合边值问题转化为求解Cauchy核奇异积分方程,利用Laplace数值反演技术计算了动态应力强度因子.讨论了材料特性和结构的几何尺寸对动态应力强度因子的影响.结果表明,随着界面相厚度的增加,无量纲化的动态应力强度因子减小.当裂纹靠近剪切弹性模量大的材料时,无量纲化的动态应力强度因子增大,反之减小.界面相两侧不同的材料组合对裂尖动态应力强度因子的影响是随着剪切弹性模量和质量密度的比值的增加而减小.界面相中裂纹长度对裂尖动态应力强度因子的影响比其他因素的影响大.     

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

正~~     

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

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

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

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

人造地震动合成中的位移误差分析

三角函数级数法是合成人工地震动常用算法之一,但是通过对加速度积分求取位移时,却存在与零线漂移相类似的位移漂移现象。