基于岩石蠕变及D-P准则的深部巷道围岩弹塑性分析

为深入分析岩石蠕变及中间主应力对深部高应力软岩巷道围岩稳定性及分区变形特性的影响机理, 通过剖析长期荷载作用下深部巷道围岩的变形特征, 揭示了基于围岩蠕变特性的深部巷道围岩四分区变形机理, 并依据Druck-Prager准则及关联性流动法则考虑了中间主应力、岩体扩容及应变软化的影响, 推导出巷道围岩各分区应力、变形和半径解析解. 结合工程算例, 通过对现场监测数据以及不同力学模型的计算结果进行对比分析, 论证了本文理论的科学性, 并揭示了不同强度准则和围岩参数对深部巷道围岩各分区形态的影响规律. 结果表明, 忽视深部巷道围岩的蠕变特性将导致初始黏聚力的取值大于实际值, 从而造成围岩的理论承载能力大于实际承载能力; 在[0, 0.7]区间内提高中间主应力系数可有效控制围岩塑性区及破碎区范围的扩展; 当初始黏聚力及内摩擦角减小时, 中间主应力系数对围岩塑性区、破碎区范围及巷道变形的影响力显著提高. 研究成果可为地下工程支护设计及围岩稳定性评估提供理论参考.      

沙曲矿区地应力及采场地压的综合监测与分析

作为典型复合型顶板的矿山,沙曲矿区的工程地质条件受到多条构造带影响,采煤工作面顶底板均为不稳定性软岩,并含少量裂隙水。为了调查该矿区的地应力条件及地压规律,本文在沙曲矿区的主要开拓巷道中开展了6个地点的地应力测试,结果表明该矿区的地应力场比较均匀,属于准静水应力场。采用现场应力套孔解除和室内围压率测定的方法得到最大水平主应力数值在20.11～22.87(MPa)之间,侧压系数在2.25～2.85范围内,属较高应力; 针对24101等工作面的大采高综采巷道进行了煤柱支承压力的现场监测,结果表明24101工作面的巷道煤柱应力分布特点可以分为塑性破坏区、峰值应力区和原岩应力区; 同时,对南翼14204工作面顶板进行了离层监测和锚杆测力计载荷观测,结果进一步揭示了随着时间的推移,巷道顶板稳定性的衰减主要是由于离层数量的增加所造成的规律。该地应力及采场地压的综合监测及分析结果对沙曲矿区的支护设计和安全生产具有重要指导和实用价值。     

巷道围岩塑性状态判定分析方法

本文依照莫尔-库仑屈服条件和非关联流动准则, 从地应力场及地下巷道开挖工程特点出发, 提出对围岩塑性屈服状态进行判定的方法。对无支护巷道和水泥锚杆支护巷道围岩进行有限元分析表明, 巷道开挖以后及时支设锚杆, 对出现塑性屈服的围岩起到了等效提高其刚度和强度的作用, 从而增加了围岩稳定性。与现场实测结果对比说明, 非关联流动准则比关联流动准则更符合实际。     

基于Hoek-Brown准则的破裂围岩应力分析

深部煤矿巷道在高围压的作用下,围岩普遍出现破裂,而且围岩破裂范围在扩大的同时,往往出现继续破坏的现象,并形成不同的分区.论文针对长的圆形巷道,将围岩分成破裂、塑性和弹性三个区域,采用Hoek-Brown准则,进行非关联弹塑性分析,获得了其应力和变形的封闭解析解,并说明该解析解是不唯一的.通过利用在弹塑交界处应力连续的条件、以及在破裂和塑性交界处径向应变连续的条件,获得了确定围岩破裂区和塑性区半径的解析算式.最后,通过算例并分析了破裂区和塑性区应力和应变的分布特点、以及破裂区范围的影响.利用所获得的结果,可以为巷道的稳定性分析以及支护设计提供理论依据.     

锚杆支护方式对深埋洞室隔冲击效果影响的数值模拟研究

与浅埋洞室相比,随着埋深增加而产生的高地应力现象使深埋洞室安全性评价更加复杂。本研究在考虑地应力平衡的前提下,采用有限元与无限元相结合的方法,对冲击荷载作用下深埋洞室的锚杆支护方式进行了设计与优化。计算结果表明：在平衡初始地应力之后,Mises应力会随着锚杆长度、疏密程度的增加略有增长。当锚杆长度和疏密程度增至某一特定值时,洞室洞壁的位移和塑性应变将达到最小值,洞室围岩将达到相对安全状态。锚杆间隔布置下,长锚杆数量的增加将提高洞室围岩的安全性。当对洞室拱脚进行重点支护时,整个洞壁的位移均有所减小。     

矩形巷道围岩弹性变形能积聚特征分析

煤矿冲击地压主要发生在巷道中,其主要原因之一是巷道围岩积聚了大量的弹性能。为得出矩形巷道围岩弹性变形能积聚特征,降低巷道支护成本,推导了巷道冲击破坏失稳能量准则,并建立了矩形巷道围岩能量积聚计算模型,理论分析了采深、巷道断面尺寸和煤层厚度对矩形巷道围岩能量积聚影响规律,得出:矩形巷道积聚的弹性能随采深的增加而增大,采深越深,巷道积聚的弹性能增长速率越快。巷道围岩积聚能量随巷道断面尺寸增加而增大。当煤层厚度小于巷道影响范围时,巷道积聚能量随煤层厚度增加而增大。在实际工程中,尽可能减小巷道断面尺寸,尽可能沿顶、底板布置巷道。研究结果为冲击地压巷道布置和降低巷道支护成本提供了理论依据。     

ANALYSIS OF ACCUMULATION CHARACTERISTICS OF ELASTIC DEFORMATION ENERGY OF SURROUNDING ROCK IN RECTANGULAR ROADWAY1)

煤矿冲击地压主要发生在巷道中,其主要原因之一是巷道围岩积聚了大量的弹性能。为得出矩形巷道围岩弹性变形能积聚特征,降低巷道支护成本,推导了巷道冲击破坏失稳能量准则,并建立了矩形巷道围岩能量积聚计算模型,理论分析了采深、巷道断面尺寸和煤层厚度对矩形巷道围岩能量积聚影响规律,得出：矩形巷道积聚的弹性能随采深的增加而增大,采深越深,巷道积聚的弹性能增长速率越快。巷道围岩积聚能量随巷道断面尺寸增加而增大。当煤层厚度小于巷道影响范围时,巷道积聚能量随煤层厚度增加而增大。在实际工程中,尽可能减小巷道断面尺寸,尽可能沿顶、底板布置巷道。研究结果为冲击地压巷道布置和降低巷道支护成本提供了理论依据。     

超千米深井巷道围岩变形破坏时空演化规律实验研究

针对超千米深井巷道围岩大变形破坏的难题,通过相似模拟实验,分析了巷道周围的应力应变演化规律和倾角对巷道稳定性的影响,揭示了超千米深井巷道变形破坏时空演化规律。实验结果表明:由于岩层倾角的存在,巷道底板两角出现应力集中现象。随外界载荷的增大,巷道拱部左上角与底板右下角首先发生破坏,随后向围岩深部扩展,且出现大范围变形破坏,巷道左帮和右帮应力集中现象加剧;巷道围岩各位置处产生的位移量差别较大,巷道拱部与底板产生位移量较大,两帮位移量相对较小,巷道变形量最大的位置达到30mm以上;超千米深井巷道围岩变形破坏特征分为3个阶段:缓慢变形阶段、剧烈变形阶段和严重变形破坏阶段,巷道围岩产生大变形非对称性破坏,且破坏范围较大,顶底板最大破坏深度达到34mm,巷道围岩失去承载能力,巷道空间完全封闭。将上述实验结果应用于工程实践,提出华丰煤矿超千米深井巷道全断面高强预应力锚网支护,特殊部位加强支护的方案,达到了理想的效果。     

深部回采巷道围岩零应变交界圈现象的模型试验研究

以淮南矿区口孜东矿深部回采巷道为工程背景,采用真三轴模型试验装置模拟了从浅埋静水压力、深埋静水压力、初掘采动应力(K=1.5为1.5倍900m埋深静水压力,后同)到回采动压(K=2~3)等不同应力条件下,大、小两种尺寸的直墙拱形回采巷道围岩应变状态,获得了浅埋条件下巷道围岩浅部拉应变、深部零应变,以及深埋静水压力及初掘采动应力下,巷道围岩出现"零应变交界圈"现象及深埋回采高应力下非线性大应变的围岩应变特征。应变监测结果反映了深部回采巷道围岩变形破坏的动态演化过程,为巷道的稳定性控制提供了一定的思考。     

流变岩体中让压支护作用下隧道力学行为研究

高地应力深埋软岩隧道大变形问题已成为隧道工程建设领域的突出难题. 根据高地应力深埋软岩隧道的变形特征, 基于"围岩能量吸收、变形释放"的让压支护是解决软岩隧道大变形问题的有效方法. 针对流变岩体中深埋圆形隧道在让压支护作用下的力学响应问题, 通过引入分数阶微积分理论, 采用Abel黏壶元件建立了改进的分数阶Burgers蠕变模型来表征围岩的时效变形. 此外, 通过在让压支护不同变形阶段引入刚度修正系数, 克服了传统支护未能考虑围岩变形释放的问题. 据此, 本文推导了在考虑支护延迟安装影响下, 不同变形阶段围岩与让压支护相互作用的解析解. 为了验证理论研究的正确性, 对一算例进行了不同解答及工程结果的比对, 吻合较好. 最后, 参数研究结果表明: 围岩与让压支护间的相互作用受蠕变本构模型分数阶阶数影响较大. 隧道的位移或支护压力与让压位移、支护刚度修正系数间存在线性比例关系, 但由于刚度修正系数仅保持在较小的变化范围内, 隧道的位移或支护压力变化并不显著.      

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

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

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

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