基于三维离散元堆积碎石土细-宏观力学参数相关性研究

以室内三轴压缩试验为基础,运用PFC3D离散元软件生成含有不同形状碎石块的1∶1堆积碎石土数值模型,标定得到相应的细观参数(相对误差控制在10%以内)。改变细观参数,定量探讨其与宏观力学参数(粘聚力和内摩擦角)关系,得到细观参数与粘聚力、内摩擦角的细-宏观参数关系式,并对其可靠性进行了验证。研究发现,(1)粘聚力与接触粘结力、颗粒摩擦系数呈线性正相关;同时受粘结力比值M=σc/τc的影响,当M1时,粘聚力随着M的增加而线性增加;M≥4时,粘聚力趋于稳定。(2)内摩擦角与颗粒摩擦系数呈对数关系,与接触粘结力呈抛物线关系。(3)细-宏观参数关系式是在含水率为7%室内试验结果基础上研究得出,将关系式运用到含水率为9%和11%的室内试验结果进行标定验证:依据含水率为9%和11%室内试验得到的c和φ值,运用细-宏观参数关系式计算得到细观参数,经验证其数值-室内试验结果相对误差在10%以内,说明通过细-宏观力学关系式得到的细观参数符合精度要求,关系式可靠。为后续土石混合体数值试验研究者进行细观参数标定工作提供了参考依据。     

黏性土宏细观参数相关性研究

运用颗粒流程序中的bond原理结合Mohr-Coulomb破坏准则对黏性土双轴试验进行了数值模拟,探讨了黏性土宏观强度指标与其对应的细观参数之间的关系.分析结果表明,黏性土内摩擦角与颗粒摩擦角呈线性关系,黏性土黏聚力-颗粒接触强度两者拟合曲线的斜率与颗粒间摩擦系数呈幂函数关系,且两者的相关性都很好.通过该宏细观关系,得出上海第②层褐黄色粉质黏土细观参数,对其进行了双轴试验的颗粒流模拟,计算得到的应力应变曲线与实测结果一致.研究成果对于黏性土细观参数的取值及颗粒流模拟具有指导意义.     

碎石料直剪实验的组合颗粒单元数值模拟

通过构造三维组合颗粒单元来描述颗粒间的互锁效应,对非规则颗粒材料的力学行为进行了离散元数值模拟,并通过碎石料的直剪实验进行了验证.该组合颗粒的质量与碎石块具有相同的概率分布特性,其几何形态则由不同数目、镶嵌尺寸、组合方位和粒径的球形颗粒进行随机构造.组合颗粒单元在局部与整体坐标之间的转动、力矩和方位关系通过四元素方法进行确定;颗粒之间的作用力采用具有Mohr-Coulomb摩擦定侓的Hertz-Mindlin 非线性接触模型,并考虑了非线性法向粘滞力的影响.在不同的法向应力下,对碎石料在直剪实验中的剪切应力和剪胀现象进行了离散元模拟,计算结果与实测结果相吻合;此外,在不同的法向应力和接触摩擦系数下,对碎石料的有效摩擦系数进行了计算和讨论.本文工作验证了组合颗粒单元在非规则颗粒材料的离散元模拟中的可行性.     

香港大屿山残坡积土的残余强度试验研究

粘性土的残余强度是边坡稳定性评价、桩基与土的相互作用机理研究及填土边坡设计中的重要参数。本文在综述大量文献的基础上，结合香港大屿山火山岩风化残坡积土的残余强度试验研究，分析了残余强度的测试方法和影响残余强度的因素。研究结果表明，残余强度与有效法向应力间具有明显的非线性关系；与单剪测试结果相比，多级剪测试结果明显偏高。     

颗粒摩擦对颗粒材料剪切行为影响的试验研究

通过对一种类似于土的颗粒材料--玻璃珠开展一系列室内直剪试验,研究颗粒间摩擦对颗粒材料剪切行为的影响. 试验一共考虑了4种不同的摩擦情况:干燥状态、用水浸润状态、完全淹没在水中和用油浸润状态. 分析试验结果发现,与干燥状态试样相比,用油浸润能明显降低试样的剪胀性和抗剪强度,而用水浸润和淹没在水中的方法没有产生显著的影响. 此外,通过在剪胀关系式中引入可变剪胀系数来考虑颗粒摩擦对颗粒材料剪胀性的影响,并从颗粒滑动与滚动的细观机理上初步解释了颗粒滑动摩擦角对临界状态摩擦角的影响规律.     

含随机参数的非线性黏弹性有限元计算

进行了粗粒土与结构接触面单调和循环加载试验，基于宏细观测量结果, 扩展了 损伤概念以 描述该类接触面在受载过程中的物态演化, 及由于物态演化导致的力学特性从初始状态到最终 稳定状态的连续变化过程. 揭示了接触面损伤的细观物理基础主要是接触面内土的颗粒破碎 和剪切压密这两种物态演化；指出接触面的剪胀体应变可以划分为可逆性和不可逆性剪胀体 应变两部分，其中不可逆性剪胀体应变可作为接触面损伤发展的宏观量度，因此其归一化 形式可作为一种损伤因子的定义；提出了建立粗粒土与结构接触面一种损伤本构关系的基本思路.     

非饱和固化淤泥的力学及体变特性试验研究

由于降雨、水位升降等自然因素,固化淤泥作为填土材料使用时,其饱和度(基质吸力)常处于变动状态,使固化淤泥力学及体变特性发生改变,对土工构筑物的安全带来潜在隐患。鉴于此,本文通过对不同基质吸力、净围压下的三轴固结排水剪切试验,探讨了基质吸力和净围压对非饱和固化淤泥力学及体变特性的影响。结果表明:随着净围压增大,基质吸力的强度贡献减小,100kPa净围压时基质吸力的强度贡献十分显著,300kPa净围压时基质吸力强度贡献则不明显。当基质吸力小于残余值时,固化淤泥试样的抗剪强度、无侧限抗压强度、粘聚力均随基质吸力增大呈线性增长,但当基质吸力大于残余值时强度增幅均下降,说明基质吸力大于残余值时强度贡献减小。剪切过程中试样均表现为体缩且随基质吸力增大体缩变形减小。吸力摩擦角在土-水特征曲线边界效应区和过渡区之间时几乎不变,进入残余区后开始下降。     

基于离散元模型的土石混合体直剪试验分析

土石混合体是由高强度块石和低强度土体组成的一类特殊工程地质材料,其力学性质可通过直剪试验进行确定。本文针对土石混合体的细观材料特性,分别采用球形颗粒单元和非规则组合颗粒单元模拟土体和块石材料,对其在不同含石量和颗粒粘结强度下的直剪试验过程进行离散元分析。计算结果表明,土石混合体的抗剪强度随着含石量和粘结强度的增加而增加;通过不同法向应力下直剪试验的离散元分析,确定了不同含石量下土石混合体材料的内摩擦角和粘聚力。本文研究结果有助于进一步揭示土石混合体的抗剪强度特性。     

胶结砂砾石层面直剪试验及破坏行为细观研究

本文采用宏观试验和细观模拟相结合的方法研究胶结砂砾石层面在剪切过程中的破坏行为．首先进行了不同法向应力作用下的胶结砂砾石层面直剪试验,获得了不同的剪切面破坏特征：随着法向应力的不断增加,剪切破坏面凹凸起伏程度、骨料脱落现象越发明显．其次为了深入探究层面破坏现象,按照室内试验采用的骨料级配粒径建立细观颗粒数值模型,结合物理试验对模型进行参数标定,并进行数值模型的层面直剪模拟和细观分析．结果表明,数值模型可以再现宏观直剪试验层面破坏特征;法向应力越大,层面区域颗粒发生错动和翻转的数量越多;层面破坏方式为颗粒间的张拉和剪切混合破坏,裂隙均集中在层面位置,随着法向应力的增加,裂隙的集中区域逐渐由“面”向“带”转变．     

土石混合体直剪离散元数值试验研究

工程中普遍存在的土石混合体以其复杂的工程地质力学性质越来越受研究者的关注。由于"块石"的存在使得土石混合体在细观层次上呈现明显的结构性特征,这种结构性将影响甚至控制着其变形破坏机制及宏观力学特征。本文借助三维颗粒离散元分析软件YADE分别从含石量、试样尺寸、强度特性等角度对土石混合体的力学性质开展了一系列的数值直剪试验研究工作,并取得了一些有意义的研究成果:在含石量及粒度组成相同的情况下土石混合体的宏观抗剪强度及剪胀性随着试样尺寸的增加而呈现降低趋势,而在相同试样尺寸下将随着含石量的增加而增加;内部块石的存在影响着其细观应力状态,从而影响其宏观力学特性。     

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

正~~     

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

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

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

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

人造地震动合成中的位移误差分析

三角函数级数法是合成人工地震动常用算法之一,但是通过对加速度积分求取位移时,却存在与零线漂移相类似的位移漂移现象。