一种岩石损伤本构模型在地下强爆炸中的应用

为了描述地下爆炸波的传播,建立了一种考虑塑性硬化、剪胀和损伤软化效应的岩石本构模型,并应用于地下强爆炸自由场的数值计算,获得的速度和位移波形以及峰值速度衰减曲线等与国外地下试验测量数据和同类计算相比十分接近,从而验证了本构模型的有效性。     

水饱和岩石中爆炸应力波传播的数值模拟

基于连续介质力学和不相融混合物理论,假定组分间无相对运动,采用B-W-N-B有效应力准则,在屈服面中引入孔隙影响因子,提出了一种多孔含水介质流固耦合的本构模型,并给出了孔隙的演化方程,对水饱和凝灰岩介质中的爆炸应力波传播作了数值模拟。数值计算给出的应力波形与实测结果有良好的一致性。采用本文中所述流固耦合本构模型,可很好地预测水饱和岩石中爆炸波的演化规律。     

饱和土中爆炸波传播问题的数值模拟

应用流体弹塑性力学的框架描述了爆炸荷载作用下饱和土介质的本构模型,包括体积压缩关系和强度关系。应用LS-dyna软件分析了三相饱和土介质中的爆炸波传播及其与结构的相互作用问题,并与先前完成的试验结果进行了对比,两者吻合较好。饱和土中爆炸波传播的一些特殊现象被重现,表明本构模型的描述是可行的。     

土中爆炸波与地下结构相互作用计算方法研究综述

土中爆炸波与结构相互作用是地下工程结构抗爆计算的理论基础,理论研究中通常在宏观试验基础上,引入合理的介质动力本构模型分析波动问题。针对工程防护中的计算模型和解析方法,综述了岩土介质中爆炸波与结构相互作用的现状和研究进展,内容包括土的动力本构模型,自由场中爆炸波的传播、爆炸波与箱型、圆柱形和直墙圆拱型等地下结构的相互作用等。重点阐述了常用的非饱和土和饱和土的动力本构模型、考虑松弛效应的土的综合弹塑性模型,一维平面波在不均匀介质和分层介质中的传播及地下结构爆炸荷载计算方面的研究成果,分析了各种模型和计算方法的适用范围和存在问题。     

爆炸波作用下三相饱和土本构研究

基于对爆炸荷载作用下三相饱和土宏观试验现象的综述，并结合以往对饱和土中气体运动规律的分析，对一维爆炸波作用下三相饱和土本构关系进行了描述。     

混凝土中爆炸应力波衰减规律的数值模拟研究

基于Kong-Fang混凝土材料模型和LS-DYNA的多物质ALE算法,开展混凝土中爆炸波衰减规律的数值模拟研究。首先,基于已有实验数据对材料模型参数和数值算法的可靠性进行了验证,在此基础上分析球形装药在混凝土自由场中爆炸波衰减规律,利用量纲分析和数值模拟拟合了球形装药在混凝土自由场中近区爆炸波峰值应力计算公式并明确其适用范围;然后,分析装药埋深对混凝土中装药正下方不同距离处爆炸波峰值应力分布的影响,建立了耦合系数与装药埋深和测点距离之间的定量关系。结果表明:Kong-Fang混凝土材料模型可实现对混凝土中爆炸波传播衰减规律的高精度数值模拟;定义混凝土中装药质量系数和耦合常数,可定量描述装药埋深和测点距离对峰值应力耦合系数的影响;建立的混凝土中近区爆炸波峰值应力计算公式可较准确地快速预测不同装药埋深、不同测点距离和不同混凝土强度时爆炸波峰值应力。研究结果可为混凝土结构抗爆设计和爆炸毁伤评估提供参考。     

饱和和非饱和介质的弹塑性损伤模型

针对饱和和非饱和工程材料变形和强度特点,给出了一个具体的弹塑性损伤本构方程,该本构方程可以描述材料性质的劣化过程（即损伤和软化）以及饱和程度对材料变形强度特性的影响。此外,还讨论了该弹塑性损伤本构方程的有限元实施方法,编制了专用的有限元程序。对常规三轴试伯进行了有限元分析,得出了一些具体的结果。     

花岗岩介质中地下爆炸震源函数研究

结合地震反演方法和数值计算方法，研究了地下爆炸震源函数的特征以及介质特性对震源函数的影响。从水饱和花岗岩地下爆炸波数值模拟计算出发，重点分析了含水量对震源函数特性的影响。结果表明，含水量会提高地下爆炸的耦合强度，对于震源函数的稳态值，干岩要比水饱和岩石低，且稳态值随含水量的增大而增大；拐角频率随含水量的增大而降低；在水饱和花岗岩中，震源函数的高频衰减指数高于f-2；而过冲特性并不明显。     

形状记忆合金的热力学本构方程

利用形状记忆因子,直观地描述了形状记忆合金(SMA)实现其超弹性与形状记忆效应的热力学宏观过程.根据热力学基本原理并假设SMA为各向同性材料,建立了描述复杂应力状态下SMA热力学行为的三维本构方程.该本构方程具有直观的表述形式且所含材料常数均可以通过宏观实验测定,克服了现存SMA三维本构方程由于表述形式复杂和材料参数不...     

土体的本构模型和超重力物理模拟

数值模拟和物理模拟是分析土体沉降和稳定性的主要手段. 本构模型作为描述土体应力应变关系的数学表达式, 是数值模拟的基础. 土体具有碎散性, 这一基本物理特性导致了其具有压硬性、摩擦性和剪胀性, 这是土的力学特性区别于金属的主要特征, 在土体的本构模型中必须反映这3个基本特性. 传统土力学将土体的变形和强度分离考虑, 分别采用弹性理论和基于刚塑性模型的极限平衡理论分析, 虽然应用广泛, 但由于不能全面地反映土的基本力学特性, 计算结果的精度常常难以满足定量分析的需要. 剑桥模型作为第一个全面反映压硬性、摩擦性和剪胀性的弹塑性本构模型, 实现了变形和强度的统一, 能较好地描述饱和正常固结黏土的应力应变关系, 被视为是现代土力学的开端; 统一硬化模型通过引入一个独特的硬化参数进一步发展了剑桥模型, 将适用范围扩大到超固结黏土. 作者认为, 未来岩土体本构模型研究的挑战是: 如何考虑岩土体在受力过程中土骨架相变与多场耦合, 以解决目前本构模型尚无法定量分析的能源、交通、环境和水利相关的重大岩土工程问题. 超重力物理模拟具有缩尺效应和缩时效应, 克服了常重力物理模拟中模型的应力水平低于原型的缺点, 特别适用于大尺度、长历时问题的模拟. 相较数值模拟, 超重力物理模拟的优势在于能够检验本构模型的合理性, 揭示本构模型无法描述的未知特性. 最后, 介绍了采用数值模拟和物理模拟联合分析大直径钢管桩水平受荷特性的工程案例.      

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

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

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

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