含共线刚性夹杂各向异性的平面问题

应用昨变函数方法，给出了含共是性线夹杂各向异性体平面问题的一般解；对于一个或二个夹二个夹杂的情形，给出了封闭形式的应力奇异性系数解；结果表明，应力奇异性系数与材料常数和εｘ＾∞有关，这里εｘ＾∞为无限元处ｘ方向的线应变。     

含共线刚性线夹杂压电介质的广义二维问题

基于Stroh理论，研究了含N个共线刚性线夹杂压电介质的广义二维问题，给出了介质内的复势函数、夹杂内的电场和夹杂尖端附近场奇异性系数的解析表达式。结果表明场奇异性系数取决于无限远处的应变。     

周期分布多边形夹杂奇异性应力干涉问题的研究

采用一种新型的杂交元模型和一种单胞模型来解决周期分布多边形夹杂角部的奇异性应力相互干涉的问题。新型杂交元模型是基于广义Hellinger-Reissner变分原理建立的,其中奇异性应力场分量和位移场分量是采用有限元特征分析法的数值特征解得到的。使用当前的新型杂交元模型,只需要在夹杂角部邻域的周界上划分一维单元,避免了像传统有限元模型那样需要划分高密度二维单元。文中给出了代表奇异性应力场强度的夹杂角部广义应力强度因子数值解,并考虑材料属性、夹杂尺寸和夹杂位置关系的影响。算例中,考虑了夹杂和基体完全接合的情况,并给出了考核例。结果表明:当前模型能得到高精度数值解,且收敛性好;与传统有限元法和积分方程方法相比,该模型更具有通用性,为非均质材料的细观力学分析打下了基础。     

电致伸缩材料内电极附近应力奇异性研究

应用复变函数的方法,研究电致伸缩材料内置电极附近的应力奇异性.基于精确的电边界条件,采用Hilbert理论以及复变函数中的Cauchy积分与留数定理,首先分别给出了柔性电极和刚性电极的复势函数解,然后就这两种极限情况,讨论了电极刚度对应力场奇异性的影响.研究结果表明:无论对柔性电极还是刚性电极,Max-well应力的应力场均呈现r-1阶的奇异性,但对于前者总应力奇异性系数为零,而对于后者总应力奇异性系数与基体的材料常数有关.     

线夹杂和线夹杂的相互作用

以短纤维复合材料为工程背景，本文利用线夹杂的工程计算模型以及无限平面中单夹杂的基本解，导出了线夹杂和线夹杂相互作用的平面问题的奇异积分方程。给出了夹杂端点的应力强度因子和夹杂界面应力的表达式，并作了具体的数值计算。     

QFP结构中引脚焊缝界面端的奇异性热应力问题研究

由于引脚、印制电路板和焊接剂的热-机材料属性不同,在受到热载荷或机械载荷时,引脚焊接界面端会产生奇异性应力,有可能产生界面开裂.为了基于界面端奇异场来评价QFP结构引脚界面端力学行为,本文拟采用数值方法求解引脚焊缝任意角度尖劈界面端的应力强度系数.具体步骤为:首先,基于高次内插有限元特征分析法确定两相任意角度尖劈界面端的奇异性指数和应力角分布函数,并引入常数热应力项,获得热-机耦合奇异性应力场表达式;采用有限元分析技术和最小二乘拟合法来获得应力强度系数的数值解.文中考察了热-机材料属性对热载荷下焊接剂/印制电路板界面端应力强度系数的影响,并给出改善界面端热应力状态的建议.     

任意形状热夹杂位移场的三角形单元离散算法

平面夹杂模型在纤维增强型复合材料中有广泛应用.复合材料内部通常含有不规则形状夹杂,而夹杂物的存在能严重影响材料的机械力学性能,往往导致应力集中及裂纹萌生等失效先兆.先前关于多边形夹杂的研究大多数关注受均匀本征应变下的应力/应变解,而对位移的分析较少. 基于格林函数方法和围道积分,本文给出了平面热夹杂边界线单元的封闭解析解,可方便应用于受任意分布本征应变的任意形状平面热夹杂位移场的数值计算.当夹杂受均匀本征应变时, 只需将该夹杂边界进行一维离散,因而本文方法可直接得出受均匀分布热本征应变的任意多边形夹杂位移场的封闭解析解.当夹杂区域存在非均匀分布本征应变时,可将该区域划分为足够小的三角形单元进行数值计算. 众所周知,应力应变场在多边形夹杂顶点处具有奇异性,容易导致数值计算上的处理困难及相应的数值稳定性问题; 然而本文工作表明,在多边形顶点处位移场是连续有界的, 因而数值稳定性较好.本文算法可以便捷高效地通过计算机编程实现. 文中给出的验证算例,均体现了本文离散方法的高精度、以及计算编程的鲁棒性.      

动态裂纹尖端弹粘塑性场的研究

本文利用弹一粘一塑性材料力学模型，对动态扩展裂纹尖端的指数奇异性和对数奇异性进行了渐近分析。文中假定，弹性阶段的粘性效应可以略去，仅在塑性应变中粘性才起作用，对于这种模型，推导出了其率敏感型的本构关系。以Ⅱ型裂纹为例，进一步推导了两种奇异性下裂纹尖端场的渐近微分控制方程，并进行了数值仿真分析。同时讨论了粘性系数α、马赫数M^2对裂纹尖端应力应变场的影响，即，弹粘塑性材料扩展裂纹的奇异性取决于其粘性系数和马赫数，粘性系数较大时，裂纹尖端场具有对数奇异性；粘性系数较小时，裂纹尖端场具有指数奇异性。修正了文献中对数奇异性区域的大小；解释了文献中过渡区的成因；给出了过渡区尖端应力场解的形式，从而建立了裂纹尖端场的统一解。     

任意形状热夹杂位移场的三角形单元离散算法

平面夹杂模型在纤维增强型复合材料中有广泛应用.复合材料内部通常含有不规则形状夹杂,而夹杂物的存在能严重影响材料的机械力学性能,往往导致应力集中及裂纹萌生等失效先兆.先前关于多边形夹杂的研究大多数关注受均匀本征应变下的应力/应变解,而对位移的分析较少.基于格林函数方法和围道积分,本文给出了平面热夹杂边界线单元的封闭解析解,可方便应用于受任意分布本征应变的任意形状平面热夹杂位移场的数值计算.当夹杂受均匀本征应变时,只需将该夹杂边界进行一维离散,因而本文方法可直接得出受均匀分布热本征应变的任意多边形夹杂位移场的封闭解析解.当夹杂区域存在非均匀分布本征应变时,可将该区域划分为足够小的三角形单元进行数值计算.众所周知,应力应变场在多边形夹杂顶点处具有奇异性,容易导致数值计算上的处理困难及相应的数值稳定性问题;然而本文工作表明,在多边形顶点处位移场是连续有界的,因而数值稳定性较好.本文算法可以便捷高效地通过计算机编程实现.文中给出的验证算例,均体现了本文离散方法的高精度、以及计算编程的鲁棒性.     

压电螺型位错和共线界面刚性线夹杂的干涉效应

研究了压电材料中压电螺型位错和共线界面导电刚性线夹杂的电弹干涉效应．运用复变函数解析延拓技术与奇性主部分析方法,获得了该问题的一般解答．作为算例,求出了界面含一条刚性线夹杂时两种压电介质区域广义应力函数的封闭形式解．导出了作用在位错上的像力和刚性线夹杂表面剪应力和电位移的解析表达式．讨论了界面刚性线长度,两种材料的剪切模量比和压电系数比对位错力和刚性线表面剪应力的影响规律．为进一步研究该类问题提供了一个基本解。     

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

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

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

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