Ⅰ-Ⅱ复合型尖V形切口脆断准则

基于最大周向应力和Seweryn-Novozhilov准则,用线弹性奇异应力场,给出Ⅰ-Ⅱ复合型尖V形切口的脆性断裂准则;并把裂纹作为切口张角为零的V形切口,这样就把V形切口问题与裂纹问题的断裂准则统一起来.为了验证该准则,本文采用有机玻璃板材加工多种V形切口试件进行了拉伸实验;同时,应用该准则对文献[16]中复合载荷下双边切口试件进行了起裂方向和临界载荷预测,并与其实验结果相比较.结果表明,本文所给出的Ⅰ-Ⅱ复合型尖V形切口脆性断裂准则与实验是相吻合的,适合于实际工程应用。     

与界面垂直相交V形切口的边界元分析方法

建立了边界元法计算各向同性结合材料中与界面垂直相交V形切口奇异应力场的分析方法。首先将V形切口尖端附近区域的位移场和应力场用Williams渐近展开式表达,将其代入弹性力学基本方程中,由插值矩阵法获得应力奇异性指数及其对应的位移函数;然后在V形切口尖端区域挖取一个小扇形域,将该扇形区域的位移场表示为有限项奇性指数和特征角函数的线性组合,并对挖去小扇形域后的剩余结构建立边界积分方程;最后将扇形区域位移场表达式和边界积分方程联合求出其切口尖端位移场的组合系数,从而获得各向同性结合材料中与界面垂直相交V形切口尖端的应力强度因子。本文的计算结果与现有结果对比吻合良好,表明了本文方法的有效性。     

插值矩阵法分析双材料平面V形切口奇异阶

对二维V形切口问题提出奇异阶分析的一个新方法.首先,以V形切口尖端附近位移场沿其径向渐近展开为基础,将其线弹性理论控制方程转换成切口尖端附近关于周向变量的常微分方程组特征值问题,然后将数值求解两点边值问题的插值矩阵法进一步拓展为求解一般常微分方程组特征值问题,插值矩阵法是在离散节点上采用微分方程中待求函数的最高阶导数作为基本未知量.由此,V形切口的应力奇性阶问题通过插值矩阵法获得,同时相应的切口附近位移场和应力场特征向量一并求出.     

非奇异应力对中央含V形切口试样断裂性能的影响

采用Williams渐近展开式表达V形切口尖端附近区域的位移场和应力场,将其代入弹性力学基本方程中,应力奇异性指数及其对应的位移和应力角函数由求解常微分方程组获得。由于在远离切口尖端的区域无应力奇异性,将切口尖端应力奇异性区域移出后,应用边界元法分析无应力奇异性的剩余结构;将Williams渐近展开式与弹性力学边界积分方程结合,解出切口尖端附近应力奇异性区域的各应力场渐近展开项系数,从而获得切口尖端附近区域的完整应力场;基于此,研究了非奇异应力项对中央含V形切口试样的表观断裂韧度和临界荷载预测值的影响。结果表明:考虑非奇异应力项时,脆性断裂的表观断裂韧度和临界荷载的预测值要比忽略非奇异应力项时的预测值更接近实验值。     

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

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

压电切口张开角和深度对其尖端力电损伤场的影响

基于三维各向异性压电损伤本构理论,导出了广义平面应力问题的损伤本构方程,并据此分析了压电薄板板边V形切口尖端附近的力电损伤,研究了切口张开角和深度对切口尖端力电损伤的影响规律．结果发现：和张开角对切口尖端损伤的影响相比,深度的影响更为明显;在张开角对切口尖端力损伤的影响规律方面,压电材料与一般弹塑性材料存在明显差异,原因在于压电切口尖端力电载荷比会随着深度的改变发生很大变化;不同深度下张开角与切口尖端力、电损伤关系曲线随着张开角的增大由发散逐渐会聚,不同张开角下深度与切口尖端力、电损伤关系曲线随着切口加深由会聚逐渐发散,并且电损伤曲线表现得更为明显．     

V形切口应力强度因子的一种边界元分析方法

将V形切口结构分成围绕切口尖端的小扇形和剩余结构两部分. 尖端处扇形域应力场表示成关于尖端距离$\rho$的渐近级数展开式,从线弹性理论方程推导出了一组分析平面V形切口奇异性的常微分方程特征值问题,通过求解特征方程,得到前若干个奇性指数和相应的特征向量. 再将切口尖端的位移和应力表示为有限个奇性阶和特征向量的组合. 然后用边界元法分析挖去小扇形后的剩余结构. 将位移和应力的线性组合与边界积分方程联立,求解获得切口根部区域的应力场、应力幅值系数和整体结构的位移和应力. 从而准确计算出平面V形切口的奇异应力场和应力强度因子.      

基于应力比值法的V形切口应力强度因子分析

本文基于有限元分析技术建立了一种应力比值方法,用于计算V形切口的应力强度因子。该方法不需要在V形切口尖端采用反映应力奇异性的奇异单元。求解时,首先给定参考问题的广义应力强度因子,然后利用待求问题的应力值与参考问题的应力值之间的比值来求解待求问题的广义应力强度因子。算例采用切口尖端应力方法分析了平板的V形切口问题。计算结果表明,该方法计算精度较高,能够方便地用于求解相关的工程问题。     

Ⅰ-Ⅱ复合型裂纹等ε_θ线体积应变能断裂准则

利用裂纹尖端附近能量的变化去建立Ⅰ-Ⅱ型裂纹断裂准则是有效的方法之一。本文利用Ⅰ-Ⅱ型裂纹尖端周向应变ε_θ所围成的区域内体积应变能的变化,建立了基于等ε_θ线内体积应变能Ⅰ-Ⅱ复合型裂纹断裂准则;讨论了Ⅰ-Ⅱ复合型裂纹断裂中的两种应力状态下的临界荷载比值;并将准则的预测值与最大正应力理论和应变能密度因子准则预测值以及已有部分材料的实验结果进行了对比分析。结果表明,本文的理论预测值与铬锰钢、混凝土材料以及PMMA材料的实验值相比基本吻合。说明该准则对于工程材料的复合型断裂问题分析具有一定的理论意义和参考价值。     

U形切口的有限断裂准则

有限断裂力学准则综合了应力和能量参数,假设裂纹或切口端部有限裂纹长度的增长．特别地,该有限裂纹的长度不是材料的基本常数,而是与构件的结构有关．基于U形切口两种形式：点方式和线方式有限断裂准则,对文献中的铝合金U形切口三点弯曲断裂实验进行了分析验证．一方面基于材料的断裂韧度和抗拉强度,预测切口件断裂载荷;另一方面根据几组不同的切口根部半径及其对应的临界切口应力强度因子,同时估算材料的断裂参数：断裂韧度和抗拉强度．将点方式和线方式两种不同形式有限断裂准则的预测结果,与平均周向应力准则、最大周向应力准则以及文献中相关结果进行了比较得出：无论是预测断裂载荷还是估算材料断裂参数,线方式有限断裂准则,与文献中相关结果比较吻合,尤其是估算的断裂韧度精度较高．     

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

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

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

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