考虑损伤界面的多相复合材料总体平均力学性能的预测

在广义自治有限元迭代法的基础上研究了损伤面对多相复合材料的总体平均横向力学性能的影响。给出了对于不同体分比和损伤界面性质的典型复合材料的数值结果。数值分析表明损伤界面的厚度及其损伤程度对复合材料总体平均横向平面应变体积模量和剪切模量有显著影响。本文还计算了损伤界面区域的局部应力场。     

内聚力界面单元与复合材料的界面损伤分析

推导了一种基于内聚力模型无厚的界面单元,用来模拟复合材料纤维与基体之间的界面层．研究了纤维周期分布的复合材料受横向荷载时,在界面不同的强韧性条件下其界面损伤演化的规律和对复合材料整体性质的影响     

纤维复合材料损伤过程的数值模拟

利用界面断裂力学和有限元法数值模拟纤维增强复合材料的细观损伤过程，研究各种主要破坏模式之间的相互转变和影响，指出以断裂能和混合度表示的界面性能是控制复合材料损伤过程的主要细观参数。分析了界面韧度对破坏性能的影响，探讨了基于破坏模式控制的复合材料韧度设计的新途径。     

基于构型力内变量的界面损伤问题研究

本文基于材料构型力的基本理论和损伤力学中含有内变量的热力学框架,提出了新的损伤变量定义方式,为研究界面损伤问题提供了一种新思路。首先,基于双相弹性体的能量分析,给出界面材料构型力表达式,通过构型力的离散化方法,实现了其在有限元中的数值计算。其次,定义构型力为界面损伤内变量,进而提出一种新的损伤演化模型,并采用刚度劣化的方法,对该界面损伤模型进行数值实现。最后,通过对复合材料界面损伤问题（有裂纹或无裂纹）进行数值模拟,分析了其界面损伤发展趋势,探讨了此模型的合理性和优越性。基于构型力内变量的界面损伤模型,可为研究复合材料的界面损伤失效问题提供一种普适性的方法。     

材料强度对层板抗弹性能和损伤特性的影响

根据纤维增强复合材料宏细观结构特征,基于ABAQUS软件平台,建立了层合板高速冲击损伤三维有限元分析模型。该模型在复合材料层间引入界面单元模拟层间分层,采用三维粘弹性本构,结合Hashin失效准则模拟单层板面内损伤.利用该模型,深入研究了复合材料层板的抗弹性能和损伤特性,数值分析结果与实验结果吻合良好,证明了该方法的合理有效性。通过数值分析,详细探讨了材料强度参数对层板抗弹性能和损伤特性的影响规律,获得了一些有价值的结论。     

界面对多相介质平均弹性性能和应力场的影响

在建立多相介质夹杂问题的自治模型基础上，应用自治有限元法研究了多相介质的平均弹性性能问题，通过数例讨论了在不同体分比下，界面对多相介质总体平均弹性性能的影响，结果表明，界面厚度对多相介质横向平均剪切模量和横向平均体积模量有明显的影响，而对纵向平均拉伸模量影响不大。为了进一步研究界面破坏特征，本文还分析了多相介质界面区域的应力场。     

复合材料层板分层伴以横向裂纹扩展的三维有限元分析

本文采用三维有限元模型,对典型铺设的[0_2/±45_2/90_2]_s碳/环氧复合材料层板中分层伴以横向裂纹的产生和扩展导致的层间应力分布进行了分析。计算结果表明层间裂纹首先在90°层中部出现并开裂至相邻界面处而产生横向裂纹,横向裂纹的出现引起局部分层按三角形状扩展;并指出分层损伤过程是一个主导性的稳定扩展过程,这是导致刚度下降的主要因素。最后,数值计算结果与实验结果比较,两者是吻合的。     

金属基复合材料的强度与损伤分析

用细观计算力学的方法分析了金属基复合材料(MMC)多重损伤与强度的关系,采用唯象的内聚力模型模拟纤维/基体界面的脱粘和采用G-T模型描述韧性基体的损伤.并用上述模型分析了长纤维增强MMC在横向荷载作用下损伤演化的规律,讨论了不同界面性质与材料强度及损伤、破坏模式之间的关系.     

金属基复合材料和强度与损伤分析

用观察计算力学的方法分析了金属基复合材料（ＭＭＣ）多重损伤与强度的关系,采用唯象的内聚力模型模拟纤维／基体界面的脱粘和采用Ｇ－Ｔ模型描述韧性基体的损伤。并用上述模型分析了长纤维增强ＭＭＣ在横向荷载作用下损伤演化的规律,讨论了不同界面性质与材料强度及损伤、破坏模式之间的关系。     

复合材料刚度性能表征的协同损伤力学模型

针对复合材料层合板的弥散型损伤,提出一个刚度性能表征的协同损伤力学模型. 该模型兼顾了微观物理损伤响应和宏观材料刚度性能表征. 从微观角度,建立细观RVE 模型求解裂纹表面张开位移和滑开位移,以此定义损伤张量,并在宏观上通过对材料应变和损伤表面位移进行均匀化处理,建立单向板或层合板的损伤刚度矩阵和损伤张量之间的联系. 以基体裂纹为例,详细分析并建立了横向裂纹和纵向裂纹的损伤本构. 计算了[±θ/90₄]_S 铺层层合板中基体横向裂纹对刚度性能的影响,结果表明该方法能够准确地预测复合材料层合板由损伤导致的刚度性能衰减.      

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

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

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

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