圆柱孔收缩的线弹性-完全塑性解及其应用

针对隧道开挖引起的地层损失现象,将开挖过程简化为圆柱孔收缩问题,借鉴文献[1]关于圆柱孔收缩的计算方法引入临界状态土力学理论,给出了圆柱孔收缩问题的线弹性-完全塑性解析解,并结合隧道开挖的工程特点给出了地层损失比以及地表土体最大沉降与隧道支护压力之间的关系式。将该方法的计算结果与Grant的离心机试验实测结果进行了对比,并通过国内外的工程实例对该方法的计算结果进行了验证,证明计算结果准确合理,可为隧道开挖引起的土体变形分析提供借鉴。     

移动荷载作用下土岩组合基坑吊脚桩变形分析

以青岛地区特有的土岩组合地质条件为背景,采用现场监测和Plaxis有限元模拟相结合的方法,研究了土岩基坑中吊脚桩在龙门吊移动荷载作用下的变形规律及动力响应。结果表明：桩身变形模拟结果与实测值吻合较好,基坑的变形主要发生在基坑上部软弱土层,吊脚桩嵌岩处产生应力集中;在龙门吊移动荷载作用下,桩顶水平位移较大,但其动力响应最小,而嵌岩处水平位移较小,但其动力响应最大;嵌岩处、桩身最大正弯矩处及最大负弯矩处的土压力动力响应较大,且移动荷载刚经过时刻影响最大。研究成果可为类似土岩结合地区深基坑支护设计提供参考。     

Hybrid finite-discrete element modelling of dynamic fracture and resultant fragment casting and muck-piling by rock blast

A state-of-the-art review is conducted to highlight the fracture mechanism in rock blast and advantages and limitations of various methods in modelling it. A hybrid finite-discrete element method (FEM-DEM) is implemented to simulate rock fracture and resultant fragment muck-piling in various blasting scenarios. The modelled crushed, cracked and long radial crack zones are compared with those in literatures to calibrate the hybrid FEM-DEM. Moreover, the hybrid modelling reproduces the rock fragmentation process during blasting. It is concluded that the hybrid FEM-DEM is superior to continuous and discontinuous methods in terms of modelling dynamic fracture of rock under blast-induced impact load.     

“Transverse Action” – A model benchmark exercise for numerical analysis of the Callovo-Oxfordian claystone hydromechanical response to excavation operations

The Callovo-Oxfordian claystone (COx) is considered as a potential geological formation to host an industrial radioactive waste repository in France. A detailed understanding of the thermo-hydro-mechanical (THM) behavior of the COx is a key issue for design of different repository structures and the safety calculations of the project. More particularly, numerical modeling of induced fracture networks around drifts excavated at the main level of the Andra’s Meuse/Haut-Marne Underground Research Laboratory (M/HM URL) and short and long term behavior of the COx around these drifts are of great interest. Several constitutive models have been developed/used in the framework of the R&D and simulation programs of Andra. A model benchmark exercise has been launched since 2012 to provide an overall view of the developed models regarding the in situ observations. In this view, two series of test cases, respectively at material point scale and at drift excavation scale are defined. Different kinds of constitutive models based on the elasto-visco-plasticity concept, continuum damage mechanics, the rigid block spring method and two-scale computational homogenized model (CHM) are used within this exercise. The obtained results show that accounting for material anisotropy and strain localization treatment techniques can improve the obtained results when elasto-visco-plastic models are used. Damage mechanics based approaches and methods accounting for discontinuities through discrete elements provide also interesting insights especially when fracturing processes must be modeled. However, more efforts are necessary to improve the robustness of these kinds of approaches in the complex context of COx response to excavation works.     

Numerical 3D modeling of the effects of strain rate and confining pressure on the compressive behavior of Kuru granite

This paper deals with numerical modeling of the compressive behavior of granite rock under high strain rate dynamic loading and wide range of confining pressure. For this end, a constitutive model based on damage mechanics and viscoplasticity for rock is formulated and implemented in explicit dynamics FEM. The viscoplastic part is based on a simple power law type yield criterion that incorporates the rate-dependency with a linear viscosity term. Moreover, a Rankine type of tensile cut-off is employed. The damage part of the model is formulated with separate scalar damage variables in tension and compression. The model is calibrated for Kuru granite and validated with the experimental data from dynamic compression tests at the strain rate of 600 1/s up to 225 MPa of confining pressure. The numerical simulations demonstrate that, despite the underlying continuum modeling approach, the model captures the correct experimental failure modes, including the transition from single-to-multiple fragmentation, as well as the dynamic compressive strengths at different confining pressures.     

Numerical investigation of the opening effect on the mechanical behaviours in rocks under uniaxial loading using hybrid continuum-discrete element method

Openings including their size, shape and distribution in rock play a significant role in the performance of rock related structures. The well-established knowledge in this area can contribute to the engineering practices, for example, underground space design, planning and optimisation in Civil and Mining Engineering and wellbore stability in Drilling Engineering, among others. Thus, understanding the failure mechanism of rock with openings is theoretically and practically meaningful. Laboratory testing on rock or rock-like materials with openings have been studied extensively in the literature, which, however, primarily focuses on the cracks/fractures. In this paper, a comprehensive numerical study on the effect of non-banded openings, i.e., circular, rectangular, and triangular opening, on the rock mechanical behaviour is performed using a hybrid continuum-discrete element method. It is revealed that the proposed simulation method can reproduce reasonably the crack initiation and propagation, and predict well the change of the mechanical behaviour due to the openings. In addition, the influence of the opening shape and opening ratio (=area of opening/specimen area) on the mechanical behaviour is also investigated.     

川中高石梯-磨溪地区震旦系灯影组碳酸盐岩岩石类型及分布规律

四川盆地川中地区震旦系气藏储量规模大,储层分布面积广。在Folk(1959)碳酸盐岩分类方案与张荫本等(1996)对粘结岩分类方案研究的基础之上,结合四川盆地震旦系灯影组的具体情况,提出新的划分方案。研究表明:(1)高石梯-磨溪地区震旦系灯影组碳酸盐岩岩石类型可分为3个亚类,10个微类。该方案不同之处在于把凝块(石)云岩归位颗粒岩,同时根据藻的形态把藻云岩进一步细分;(2)沉积岩类的发育和分布受沉积层序的控制。海侵体系域(TST)沉积阶段,以泥-粉晶云岩、微生物岩等正化学岩沉积为主。高位体系域(HST)沉积阶段,以凝块云岩、砂屑云岩等异化颗粒岩沉积为主。纵向上,岩性具有从暗色泥晶云岩→藻云岩→颗粒云岩→浅色泥晶云岩演化的特征,即泻湖亚相→藻丘亚相→颗粒滩亚相→台坪亚相的规律;(3)凝块云岩、砂屑云岩等颗粒岩是主要的储集岩类;(4)灯影组碳酸盐岩不同岩类的发育和分布对储层的发育起到控制作用。对灯影组的油气勘探,除了对岩溶作用的研究之外,还需要重视沉积古地貌、岩相展布特征的分析,有助于有效的预测灯影组储层"甜点区"。     

华北盆地强震孕育的动力学机制研究

强震活跃区通常出现在地壳变形剧烈的地区.但是,华北盆地地壳变形速率小,而破坏性强震却十分频繁.文中利用有限元方法对华北盆地发震断层下方可能存在的速度异常体进行了力学分析,计算结果显示：无论这个异常区是低速异常还是高速异常,在水平构造力的作用下（无论是挤压还是引张力作用）,应力及应变能密度都在空间某些区域积聚.根据本文的工作及前人的研究结果,认为华北盆地强震孕育的物理机制如下：华北克拉通在减薄、破坏过程中,浸入的岩浆或上涌的地幔物质经过漫长的地质时期演化,成为低速异常体.当断裂带下面存在这种低速异常体时,水平向构造力的作用,会使异常体上方应力集中、应变能积聚.当异常体上方的应力水平达到某种限度,介质发生破坏.在破坏过程中,异常体释放的能量,使断裂带上的应变能集聚速度加快;另一方面,异常体的破坏使断层的闭锁面积减小.随着外部构造应力的持续加载,当发震断层上的应力及能量积累达到一定极限,断层产生错动,发生强震.同时,文中结果还显示,断层在从闭锁、积累能量到解锁、发生强震的过程中,其下方的低速异常体起着重要的推进作用.