模拟裂纹扩展的一种有限元局部动态子划分方法

提出了一种有限元子划分结合子结构的方法来模拟裂纹扩展问题。提出的方法中,将单元分为三类:被裂纹贯穿的单元,包含裂尖的单元和常规单元。对前两类单元进行子划分,每个单元的归类随裂纹的扩展而动态变化。覆盖一条裂纹的前两类单元子划分后构成一个子结构,子结构也是动态的,跟随裂纹的扩展而逐步扩大。本文的方法可以使裂纹沿任意路径扩展而不受初始网格的限制,裂纹扩展后无需对结构整体的网格重划分,结构整体分析的总自由度也不变。用该方法计算无限大平面中心裂纹的应力强度因子,模拟三点弯梁跨中裂纹的扩展,验证了计算精度,并进一步用该方法模拟了非均质材料中裂纹的扩展,考核了对复杂裂纹扩展问题的适用性。     

裂纹扩展独立于网格的一种有限元单元子划分方法

提出了一种有限元模拟裂纹扩展的单元子划分结合子结构的方法。本方法中,裂纹可以进入或穿过一个单元,或沿单元的边界扩展,因此裂纹可以沿任意路径扩展而不受初始网格的限制。对上述几类包含裂纹的单元按照裂纹的路径进行子划分,覆盖一条裂纹的所有子划分单元就组成了一个子结构,子结构规模随裂纹的扩展而增大。子结构中因单元子划分而新增的结点自由度,通过自由度的凝聚用初始网格结点的自由度表示,因此结构整体分析的总自由度不变。以上述方法为基础建立了裂纹萌生和扩展的准则。用本文的方法分析了单（双）材料无限大平面中心（界面）裂纹的裂尖场,验证了本文方法的精度,并模拟了颗粒复合材料中微裂纹在颗粒、基体和界面中逐步扩展的过程,考核了本文方法对复杂裂纹扩展问题模拟的适用性。     

光滑边域连续-非连续细胞自动机在裂纹开裂扩展中的应用

针对岩石裂纹开裂扩展问题，将应变光滑技术与连续-非连续细胞自动机方法相结合，构建了非连续裂纹贯穿单元和裂纹尖端单元光滑应变场，提出快速自适应光滑边域连续-非连续细胞自动机方法.构建了裂纹面位移非连续精细表征的非连续增强形函数，建立了裂纹贯穿单元和裂纹尖端单元的光滑应变矩阵求解方法，利用高斯散度定理将单元的区域积分转换为光滑域边界线积分，推导给出了光滑边域连续-非连续细胞自动机应变矩阵计算表达式，并提出了快速自适应更新方法，建立了加速因子随元胞更新而同步更新的自适应加速算法，基于此，最佳加速因子随更新自动获得，收敛速度较传统细胞自动机方法得到极大提高.利用C++编制了分析计算程序，针对多裂纹开裂扩展过程进行了模拟，并与扩展有限元进行了比较.研究发现，光滑边域连续-非连续细胞自动机方法在解的精确性、稳定性和收敛性上较扩展有限元有显著优势.     

断裂过程的有限元模拟

讨论了材料断裂过程的有限元模拟技术。基于自适应有限元的一般原理，并针对多相材料的裂纹扩展的特点，提出了一种简化的高精度和高效率有限元网格的动态重新划分策略。裂纹被假设沿着单元之间的路径连续扩展，利用节点力释放技术生成新的裂纹自由表面，发展了一种可随裂尖连续移动的网格动态加密和释放方法。这种方法已在各种裂纹问题中得以实现与应用。     

裂纹尖端弹塑性变形的实验测量与有限元计算

用显微网格数字图像处理方法，测量了Ｆｅ－Ｓｉ材料平面应力条件下Ⅰ型单边裂纹尖端附近的应变场，并且用弹塑性有限元程序进行了数值模拟。对实测结果和有限元计算模拟结果进行了对比分析。结果表明：逐级更新拉格朗日坐标体系下的弹塑性有限元计算在裂纹尖端不出现严重钝化和损伤的条件下，可以正确反映裂尖区域的弹塑性变形。在裂纹接近临界扩展时，有限元计算与实测结果之间有较大差别，其原因有待进一步研究。     

弹塑性力学问题的虚单元法

具有有限差分法特征的虚单元法，可视为是有限元法向任意多边形单元的扩展。在材料细观力学性能表征、非均质材料力学分析等非线性问题方面，传统的弹塑性有限元法具有网格数目多、效率低下等不足之处，而虚单元法使网格划分更加灵活，为材料的弹塑性力学分析等非线性问题提供了新的思路。基于增量法弹塑性力学原理和双线性投影算子，建立了弹塑性力学问题的虚单元法求解技术，提出了弹塑性力学问题虚单元法的应力更新方案，研究了弹性力学问题虚单元法的精度和收敛性，讨论了虚单元法求解弹塑性力学问题的网格依赖性。同时，开展了任意多边形和凹多边形单元的数值试验研究，结果表明，虚单元法无须分割多边形，仅需节点自由度便可求得单元刚度矩阵和应力等效荷载，程序实现简单，计算精度高，改善了传统有限元的网格依赖性和塑性区的网格奇异性。     

用可动节点有限元模型模拟裂纹扩展过程

本文针对裂纹扩展有限元模拟模型存在的问题,提出一种可动节点模型.使裂尖在单元尺度内任意移动,实现了扩展模拟连续性.加进松弛扩展,突出其所占比重,使模拟更加逼真.用弹塑性有限变形杂交元自编模拟裂纹扩展的通用FORTRAN 程序,对含中心裂纹铝板进行模拟,假定材料各向同性并遵守等向强化规律,采用米赛斯屈服准则和Drucker 塑性势流动理论做增量弹塑性静力分析,得到了与实验符合较好的结果.     

矩形平面应力扩展有限元及其程序实施

扩展有限元法(XFEM)是一种新的求解不连续问题的数值方法,由于该方法建立在传统有限元法框架内,而且避免了不连续界面扩展时的网格重新划分,所以,近年来扩展有限元法是工程力学领域里的一个研究热点.本文基于Zi和Belytschko新近提出的一种不连续位移场假设,推导了求解粘着裂纹扩展问题的矩形扩展有限元公式和离散方程.提出了基于位移加载的直接迭代的算法来求解非线性有限元方程,并给出了详细的过程.编写了Fortran程序,模拟了三点弯曲梁开裂问题,得到的结果与已知文献结果吻合,验证了方法、公式和所编写的程序的正确性.     

数学网格和物理网格分离的有限单元法(II):粘聚裂纹扩展问题中的应用

强化有限单元法将物理网格与数学网格分离开来,可以方便地描述非连续变形;粘聚区域模型是模拟断裂过程区作用最简单有效的方法,且可以避免裂纹尖端的应力奇异性.本文以平面问题为例,将强化有限单元法与粘聚区域模型相结合,利用富集数学节点描述任意粘聚裂纹扩展过程中的非连续变形问题,提出了裂纹扩展过程中数学节点富集和数学单元定义的方法.本文还导出了与平面4～8节点平面等参单元对应的8～16节点粘聚裂纹单元列武.最后,通过三点弯梁的裂纹扩展过程模拟验证了本文提出的粘聚裂纹扩展模拟方法的有效性.     

状态向量的扩展有限元方法研究

利用哈密顿正则方程的半解析法计算单元位移场和应力场,可以得到精度比较高的解.但此半解析法在计算应力尖峰区域时,该区域要细化网格.当裂纹扩展时,又要重新生成刚度矩阵进行求解,导致求解效率降低.利用扩展有限元处理裂纹的不连续性,当裂纹扩展时可以避免网格的重构.为充分利用状态向量方程和扩展有限元的优势,该文将两者结合起来分析材料的断裂问题:计算应力强度因子和模拟裂纹扩展.最后通过算例分析,验证了该文提出方案的可行性.     

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

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

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

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

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

正~~     

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