弹塑性材料裂纹小量稳态扩展测试方法的讨论

弹塑性材料动态断裂试验结果表明,裂纹起裂后存在一个稳态扩展过程,为了获得材料的扩展阻力,该文利用自行研制的间接杆-杆型冲击托伸试验装置,通过控制入射脉冲和消除波二次反射加载影响,成功地在冲击拉伸试验装置上实施周边切口短圆柱小试件的冲击拉伸加卸载试验,即动态起裂止裂试验,依据动态起裂止裂试验和二次加载试验的结果,采用标定曲线法获得裂纹稳态扩展的轨迹a（t）,该文还推出了裂纹稳态扩展过程的J积分计算公式,根据P-δ曲线和裂纹稳态扩展的轨迹a=a(t),由此得到了材料小量稳态扩展阻力曲线JR-△a,结果表明,JR-△a曲线不随试件尺寸的变化而变化,可以作为材料的特征曲线,从而形成了弹塑性材料裂纹稳态扩展阻力特性的J表征和测试方法。     

周边切口短圆柱试件的杆杆型冲击拉伸试验系统的弹塑性有限元分析

基于一维试验原理提出了用带有周边切口的短金属圆柱试件进行平面应变型弹塑性动态断裂韧度的测试方法；对该复杂的动力学系统进行了轴对称的弹塑性有限元分析，并计算了动态围道^J积分；根据对试件功能转换关系的分析和Rice公式的物理意义，提出了用试件两端平功载荷－两端相对位移曲线（P^-△）推广Rice公式计算试件的远场J积分，由此得到的P^-△曲线基本上消除了与裂纹运动无关的质心运动动能的影响。论证了J积分作为裂端的表征参量，且当切口深度比大于70％时，Rice公式有较高的计算精度，为平面应变型弹塑性动态断裂韧度的表征与测试提供了依据。     

一种弹塑性材料动态起裂韧度的J表征和测试方法

在数值计算的基础上提出了一个弹塑性材料动态起裂韧度的J表征和测试方法．这一方法利用自行研制的间接杆-杆型冲击拉伸试验装置，对周边切口的短圆柱试件实施基于一维试验原理的弹塑性材料动态断裂试验；利用试件两端的平均载荷—相对位移曲线(P-δ)或平均载荷-裂纹张开位移曲线(P-△)，推广Bice远场公式获得动态J积分；采用柔度变化率法确定裂纹起裂点，从而得到动态起裂韧度JJD．这一方法的优点在于利用P-δ或P-△曲线将外力对试件所做的与裂纹运动无关的质心惯性运动动能近似地从总能量中分离出来，且平均载荷P在起裂前以至失稳扩展前是单调增的，同时试验获得的P-δ和P-△曲线光滑，这使得用J积分作为裂纹尖端的表征参量以及用柔度变化率法确定起裂点和失稳点具有坚实的物理基础．     

含中心裂纹铝合金薄板裂纹稳态扩展的实验研究

一、引言对于韧性材料的含裂纹构件,通常是在大范围屈服情况下断裂的,所以必须建立弹塑性断裂理论来进行研究。用J积分判据J_(1c)或临界裂纹顶端张开位移δ_(cr)来衡量韧性材料的断裂韧性,在工程上具有实用意义。但是J_(1c)和δ_(cr)都是用来确定裂纹的初始起裂,而起裂后的裂纹稳态扩展现象很重要,特别对于硬化材料的金属薄壁构件更为明显,在裂纹缓慢稳态扩展的过程中,必须继续增加载荷,直到裂纹失稳扩展,因此要合理地确定含裂纹薄壁结构的承载能力,就需要研究裂纹稳态扩展过程。Feddersen在研究平面应力断裂问题中,用铝合金中心裂纹板试件做了大量实验,对工程设计提供了有用的分析方法,但是对于裂纹的稳态扩展过程只做了定性描述,没     

LY12CZ薄板铝材裂纹起裂、稳态扩展和失稳扩展时的抗裂特性

本文以实验为基础,研究了在弹塑性状态下,LY12CZ薄板铝材裂纹起裂,稳态扩展和失稳扩展时的行态。采用特制的涡流裂纹自动跟踪仪,较精确地得到裂纹扩展量△a_p的连续记录。在此基础上,测得真实裂纹端的张开位移COS和张开角CTOA,从而得到用它们表征的、完整的裂纹扩展阻力曲线,即COS-△a_p和CTOA-△a_p曲线。实验表明:LY12CZ薄板铝材裂纹起裂时的张开位移δ_c是一稳定的常值,在起裂后的稳态扩展阶段,张开角CTOA是一稳定的常值。本文结合理论分析后指出,对一般韧性材料,在dδ_p/da》σ_s/E时,其张开角确应是材料的抗裂特性。在实验中还发现,在裂纹扩展一定值后,实测的张开角突然下降,此现象对应于裂纹的失稳扩展。     

金属试件变形断裂过程的计算模拟分析——组合功密度模型的应用

本文作者已建立了一个既适用于裂纹体亦适用于无裂纹体的统一损伤断裂模型——组合功密度模型。在本文中,作者应用该模型并结合大应变弹塑性有限元方法对40Cr制的圆棒光滑拉伸试件、圆棒切口拉伸试件和三点弯曲裂纹试件的变形和断裂过程进行了计算模拟。模拟结果很好地再现了两种圆棒拉伸试件的实验过程;而对于三点弯曲试件,模拟得到的载荷——加载位移关系曲线、裂纹扩展量——加载位移关系曲线、J_R阻力曲线和断裂韧性J_(IC)值等均与实测结果相当吻合。证实了该模型既能用于模拟预测无裂纹体中的裂纹萌生和扩展,也能模拟预测裂纹体中的裂纹起裂和稳定扩展过程。     

Ⅰ型裂纹中低速冲击荷载下起裂韧度测试新方法

为了探寻更加合理的构型试件来研究纯Ⅰ型裂纹在冲击荷载下的起裂及扩展行为, 提出一种新构型试件, 即双倾斜底边中心裂纹试件(double inclined bottom central cracked, DIBCC)。借助于中低速落锤式冲击实验装置进行冲击实验, 通过应力波来使试件内预制裂纹起裂并扩展, 同时利用应变片测试系统监测裂纹起裂时刻, 并采用AUTODYN有限差分软件对实验过程进行数值模拟, 最后计算裂纹的动态应力强度因子, 利用实验测得的起裂时刻, 确定试件的起裂韧度。结果表明:(1)在反射拉伸波作用下, 预制裂纹两侧会产生垂直于裂纹面向外的位移, 使预制裂纹扩张, 从而使裂纹起裂。(2)数值模拟结果与实验结果在裂纹扩展路径上具有一致性, 说明本文中提出的DIBCC构型试件有效, 可以用来测试裂纹在冲击载荷下的断裂韧度。     

基于圆弧底试件的动态裂纹扩展及止裂规律研究

为了研究脆性材料的动态裂纹扩展及止裂规律，设计了一种带圆弧形底边的梯形开口边裂纹（trapezoidal opening crack with arc bottom，TOCAB）构型的试件。在落锤冲击设备加载下，对圆心角为0°、60°、90°和120°的TOCAB试件进行了冲击实验，并采用裂纹扩展计(crack propagation gauge，CPG)监测裂纹起裂和扩展时间，从而获得裂纹扩展速度。采用有限差分软件AUTODYN对落锤冲击设备和试件进行数值模拟，研究了裂纹的动态扩展过程及止裂规律。还基于实验和数值方法，计算了裂纹的临界动态应力强度因子。实验和数值结果均表明:3种弧度的TOCAB试件都可以实现运动裂纹止裂，该构型可用于研究动态裂纹止裂问题；数值计算的裂纹扩展路径与实验结果基本一致，验证了数值模型的有效性；裂纹起裂和止裂时刻的临界动态应力强度因子大于裂纹动态扩展过程中的临界动态应力强度因子。     

双边切口薄板小试件动态弹塑性有限元分析

对［１］中采用的双边切口薄板小试件进行了动态弹塑性有限元分析，计算了动态积分，研究论证了积分作为该试验系统试件裂端表征参量的可行性；对深裂纹Ｒｉｃｅ公式计算动态积分的有效性进行了验证。为文［１］提出的平面应力型动态弹塑性起裂韧度的表征与测试方法提供了进一步的论证。     

皮质骨断裂力学行为的实验研究

采用稳态裂纹扩展和疲劳裂纹扩展的实验, 研究了牛皮质骨横向和纵向裂纹扩展的断裂力学行为. 沿着两个方向制备了紧凑拉伸(CT)试件. 由于试件尺寸的限制,采用数值计算方法确定了裂纹尖端应力强度因子与裂纹长度的关系. 在实验中, 采用数字图像相关法精确测定裂纹尖端的位置. 由于裂纹沿横向扩展时有较大的偏斜, 将采用$J$积分测量其断裂韧性. 实验结果表明, 在裂纹扩展的一定范围内, 皮质骨的断裂韧度随着裂纹不断扩展而增大, 即表现出上升的阻力曲线(R-curve)性质.而皮质骨的横向裂纹扩展的断裂韧度要远远大于纵向裂纹扩展的断裂韧度, 表现出各向异性的阻力曲线行为. 在疲劳裂纹扩展中, 纵向疲劳裂纹扩展率要大于横向疲劳裂纹扩展率, 这说明皮质骨具有各向异性的疲劳裂纹扩展性质.      

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

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

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

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