A spatial estimation model for continuous rock mass characterization from the specific energy of a TBM

Summary  Basic principles of the theory of rock cutting with rolling disc cutters are used to appropriately reduce tunnel boring machine (TBM) logged data and compute the specific energy (SE) of rock cutting as a function of geometry of the cutterhead and operational parameters. A computational code written in Fortran 77 is used to perform Kriging predictions in a regular or irregular grid in 1D, 2D or 3D space based on sampled data referring to rock mass classification indices or TBM related parameters. This code is used here for three purposes, namely: (1) to filter raw data in order to establish a good correlation between SE and rock mass rating (RMR) (or tunnelling quality index Q) along the chainage of the tunnel, (2) to make prediction of RMR, Q or SE along the chainage of the tunnel from boreholes at the exploration phase and design stage of the tunnel, and (3) to make predictions of SE and RMR or Q ahead of the tunnel’s face during excavation of the tunnel based on SE estimations during excavation. The above tools are the basic constituents of an algorithm to continuously update the geotechnical model of the rock mass based on logged TBM data. Several cases were considered to illustrate the proposed methodology, namely: (a) data from a system of twin tunnels in Hong Kong, (b) data from three tunnels excavated in Northern Italy, and (c) data from the section Singuerlin-Esglesias of the Metro L9 tunnel in Barcelona. Correspondence: G. Exadaktylos, Department of Mineral Resources Engineering, Technical University of Crete, Chania, Greece     

Geostatistical Evaluation of the Mechanical Properties of Rock Mass for TBM Tunnelling by Seismic Reflection Method

Summary. The evaluation of the rock mass mechanical properties by the seismic reflection method and TBM driving is proposed for TBM tunnelling. The relationship between the reflection number derived from the three-dimensional seismic reflection method and the rock strength index (RSI) derived from TBM driving data is examined, and the methodology of conversion from the reflection number to the RSI is proposed. Furthermore a geostatistical prediction methodology to provide a three-dimensional geotechnical profile ahead of the tunnel face is proposed. The performance of this prediction method is verified by actual field data.     

土岩组合岩体中抗拔桩极限承载力的确定

为全面探究土岩组合岩体中抗拔桩的极限承载力,结合工程岩土参数及试验数据,运用Flac3D数值分析软件对其进行数值模拟分析,即可得到土岩组合岩体中抗拔桩的极限承载力。基于被动状态下的Kotter极限平衡方程式求解土层提供的抗拔力,根据岩石强度,基于Hoek-Brown破坏准则求解抗拔桩嵌岩端岩体的抗拉强度,从而可计算得到嵌岩端岩体的抗拔力;由静力平衡原理,叠加土层及嵌岩端岩体提供的抗拔力及破坏锥体重量,即可得到土岩组合岩体中嵌岩抗拔桩的极限承载力理论解析式。在嵌岩深度较小的情况下,该解析式的理论计算值与数值模拟分析值相接近,但随着嵌岩深度的增加,理论计算值会偏离数值计算值。故结合数值模拟试验值,对提出的极限承载力理论解析式作进一步的修正,得到修正后的极限承载力解析式能反映嵌岩端岩石风化程度、嵌岩深度、土层厚度、桩长对极限承载力的影响。运用修正后的解析式对该地质条件下不同抗拔桩的极限承载力计算表明：数值模拟结果与理论计算结果相吻合,说明所建立的抗拔桩极限承载力解析式的方法是可行的。同时,运用该方法可确定类似工程中嵌岩抗拔桩的极限承载力。     

TBM滚刀破岩过程影响因素数值模拟研究

全断面岩石掘进机（TBM）的破岩效率主要受刀盘设计和岩体特征的影响。采用颗粒流方法建立了岩石试件与滚刀的数值模型,对TBM滚刀破岩过程的影响因素进行了分析。研究表明,单刃滚刀交替作用下强度较低的岩石中易形成规则的张拉裂纹而生成较大岩碴;强度较高的岩石中滚刀的侧向挤压促使形成块度较小的片状岩碴;双刃滚刀作用下岩石表面受到强烈挤压后出现较大的拉应力,使岩石更易破碎,且仅在强度较高的岩石中才易形成透镜状岩碴。滚刀破岩过程中存在能耗最小的最佳间距,该最佳间距随着岩石强度的增大而减小。滚刀破岩过程中,结构面对裂纹扩展具有显著的控制性作用,并阻隔损伤向结构面下的岩石中渗透;随着结构面与滚刀侵入方向夹角的减小,结构面将引导裂纹向岩石深部扩展,而当夹角较大时,结构面则会引导裂纹横向扩展,易导致大块岩碴的形成     

复杂矿区三维地质可视化及数值模型构建

复杂三维地质可视化和数值模型的构建是采矿、岩土及水利水电工程中常常遇到的关键问题。针对数值模拟软件在复杂矿区地质体模型构建及单元网格划分等前处理问题上存在的不直观、工作量大等缺点,结合矿业软件在复杂三维地质体建模的强大优势,在研究线框模型和块体模型差异的基础上,提出了基于线框模型的复杂矿区三维地质可视化及数值模型构建技术。根据线框模型特点不同,以地表和采空区模型为型,分别阐述了基于Surpac的表面模型和实体模型无缝转换成MIDAS/GTS数值模型的有效方法,并将该技术成功应用于冬瓜山铜矿复杂采空区群围岩稳定性研究中。结果表明,该方法是可行、有效的,从而为矿山开采设计、安全分析及复杂地质体结构的岩土工程数值模拟分析提供了一种新的建模方法。     

Three-Dimensional Back-Analysis of an Instrumented Shallow Tunnel Excavated by a Conventional Method

The excavation of a shallow tunnel induces deformations of the soil volume in the vicinity and above the tunnel and consequently on the nearby buildings. The range of these deformations depends among other on the geological conditions, the geometry of the tunnel, and the excavation method. In this context, this research focuses on the 3D numerical modeling of a shallow tunnel instrumented during its construction, located on the Toulouse (France) subway line B for which the excavation has been carried out in a conventional manner in an over consolidated molassic geological context. The objective of this analysis is to estimate the tunnel behavior in terms of vertical and horizontal movements of the surrounding soil and the deformations of the existing buildings. The explicit finite differences numerical code FLAC^3D is used to model the various implementation phases of the work where the fluid–soil interaction is taken into account through an undrained coupled analysis. The results of this 3D model are compared to those of the in situ measurements in order to validate the geotechnical characteristics of the molasses. The latter are a useful basis for the back-analysis of the different monitoring sections implemented in areas where the tunnel excavation is made by TBM with pressurized front.     

基于透明土孔隙介质双渗流模型试验研究

对土体内部渗流可视化研究,透明土材料可以作为的重要研究手段之一,能实现的岩土体内部渗流直接观测研究。利用透明土技术,建立双重孔隙型介质渗流物理模型,利用有限差分软件FLAC3D建立双重孔隙型介质结构三维数值分析模型,进行测量数据和数值模拟计算对比分析。模型试验结果和模拟结果表明,基于透明体的双渗透实验能实现渗流的可视化,以及FLAC3D数值模拟结果与实验所测的流速曲线拟合良好。     

Midas／GTS软件在边坡三维稳定分析中的应用

Midas／GTS岩土数值模拟软件具有简洁的界面、前后处理功能强大的岩土材料模型库,能满足大部分岩土体的破坏模式。建立工程边坡的三维数值模型,并对其计算分析．结果接近真实情况,且计算结果相对安全。本文采用该软件与有限元强度折减法相结合模式,对三维边坡实例的各种工况进行数值模拟,得到边坡的最大剪应变、最大剪应力出现的部位和数值大小以及边坡整体安全系数,由此确定边坡的薄弱带和潜在滑面,能客观反映边坡稳定性。     

A Combined Three-Dimensional Geological-Geostatistical-Numerical Model of Underground Excavations in Rock

Summary. This paper exploits geological and borehole geotechnical data obtained in the exploratory phase of a tunneling project to investigate in a first place if the kriging interpolation scheme may effectively reproduce the spatial variability of rock mass quality (Rock Mass Rating, RMR) in the vicinity of tunnels. For this purpose a quick solver in Fortran has been developed that performs variography analysis of 3D spatial data, fast kriging estimations of RMR between borehole sampling locations at the centroids of the elements of the numerical model, and model validation. For the purpose of an integrated underground excavation design, a step further is made by incorporating into the 3D mechanical numerical model of the rock mass, the three-dimensional (3D) solid geological model, thus coupling the geology with the ground (geotechnical) model (i.e. each element of the numerical model is assigned a geological material). The mechanical properties of each finite difference cell (or Representative Elementary Volume) of the ground model were then prescribed according to its geological type, the spatial heterogeneity of the rock mass expressed quantitatively with the kriging model, and the upscaling calculations of the mechanical properties of the intact rocks determined in the laboratory, based on the size-effect (strength dependence on size) and Damage Theory. Furthermore, a preliminary numerical simulation of the advance of unsupported tunnels in the model of the heterogeneous rock mass was performed for illustration purposes.     

Numerical Analysis of Deformation Control of Deep Foundation Pit in Ulanqab City

The deformation of a deep foundation pit in Ulanqab city was analyzed by in situ monitoring and numerical simulation data. First, a limit equilibrium method was used to evaluate the stability of foundation pit under non-supporting and ribbed plate anchored retaining wall support. The results show that the ribbed plate anchored retaining wall support could control the soil deformation effectively in deep foundation pit. Additionally, a three-dimensional finite difference model of excavation zone and reinforcement structure was established using FLAC3D software to simulate the excavation of deep foundation pit. The simulation results were verified by monitoring field data. The simulation results show that all the bolts are in the tension state, and the maximum axial force is located at the end of the bolt, with a value of 46.4 kN. Foundation pit excavation is a continuous stress release process, making the rock and soil of side wall of foundation pit slide towards the free surface. The largest deformation is 4.0 mm, near the side wall of foundation pit, consistent with the deformation monitoring results.

     

TBM Performance Analysis in Pyroclastic Rocks: A Case History of Karaj Water Conveyance Tunnel

Karaj Water Conveyance Tunnel (KWCT) is 30-km long and has been designed for transferring 16 m³/s of water from Amir-Kabir dam to northwest of Tehran. Lot No. 1 of this long tunnel, with a length of 16 km, is under construction with a double shield TBM and currently about 8.7 km of the tunnel has been excavated/lined. This paper will offer an overview of the project, concentrating on the TBM operation and will review the results of field performance of the machine. In addition to analysis of the available data including geological and geotechnical information and machine operational parameters, actual penetration and advance rates will be compared to the estimated machine performance using prediction models, such as CSM, NTNU and Q_TBM. Also, results of analysis to correlate TBM performance parameters to rock mass characteristics will be discussed. This involves statistical analysis of the available data to develop new empirical methods. The preliminary results of this study revealed that the available prediction models need some corrections or modifications to produce a more accurate prediction in geological conditions of this particular project.     

Anchoring Parameters Optimization of Tunnel Surrounding Rock Based on Particle Swarm Optimization

In the highway tunnel project, due to the uncertainty and complexity of geotechnical parameters, it is difficult to design and optimize the anchorage parameters. In order to solve this problem, based on the Zhenfengling tunnel project as the research object, 3D tunnel model was established by ANSYS software, the model was imported into FLAC^3D software for numerical analysis and calculation for the displacement and stress of surrounding rock during construction, using orthogonal experiment methods to analyze factors affecting the stability of surrounding rock. The influence of anchor length, anchor spacing, anchor diameter, the thickness and elastic modulus of spray layer on the displacement of arch waist, arch crown and arch bottom of the tunnel was obtained by using range and variance analysis. The regression model of the relationship between arch surrounding rock displacements and anchorage parameters was determined by fitting regression. The particle swarm optimization algorithm was used to optimize the anchorage parameters in combination with the tunnel section cost formula, and the parameters were compared with those not optimized. Finally, the anchor length is designed as 3 m, the anchor spacing is designed as 1.4 m, the anchor diameter is designed as 21 mm, the thickness of spray layer is designed as 24 cm, and the elastic modulus of spray layer is designed as 24 GPa. The stability requirements can be met, and the economic efficiency is also greatly improved. The cost of support with optimized parameters is 19.4% lower than that of the original design parameters.

     

Effect of spatial variability of cross‐correlated soil properties on bearing capacity of strip footing

Geotechnical engineering problems are characterized by many sources of uncertainty. Some of these sources are connected to the uncertainties of soil properties involved in the analysis. In this paper, a numerical procedure for a probabilistic analysis that considers the spatial variability of cross‐correlated soil properties is presented and applied to study the bearing capacity of spatially random soil with different autocorrelation distances in the vertical and horizontal directions. The approach integrates a commercial finite difference method and random field theory into the framework of a probabilistic analysis. Two‐dimensional cross‐correlated non‐Gaussian random fields are generated based on a Karhunen–Loève expansion in a manner consistent with a specified marginal distribution function, an autocorrelation function, and cross‐correlation coefficients. A Monte Carlo simulation is then used to determine the statistical response based on the random fields. A series of analyses was performed to study the effects of uncertainty due to the spatial heterogeneity on the bearing capacity of a rough strip footing. The simulations provide insight into the application of uncertainty treatment to geotechnical problems and show the importance of the spatial variability of soil properties with regard to the outcome of a probabilistic assessment. Copyright © 2009 John Wiley & Sons, Ltd.     

Beam–solid contact formulation for finite element analysis of pile–soil interaction with arbitrary discretization

This paper presents an embedded beam formulation for discretization independent finite element (FE) analyses of interactions between pile foundations or rock anchors and the surrounding soil in geotechnical and tunneling engineering. Piles are represented by means of finite beam elements embedded within FEs for the soil represented by 3D solid elements. The proposed formulation allows consideration of piles and pile groups with arbitrary orientation independently from the FE discretization of the surrounding soil. The interface behavior between piles and the surrounding soil is represented numerically by means of a contact formulation considering skin friction as well as pile tip resistance. The pile–soil interaction along the pile skin is considered by means of a 3D frictional point‐to‐point contact formulation using the integration points of the beam elements and reference points arbitrarily located within the solid elements as control points. The ability of the proposed embedded pile model to represent groups of piles objected to combined axial and shear loading and their interactions with the surrounding soil is demonstrated by selected benchmark examples. The pile model is applied to the numerical simulation of shield driven tunnel construction in the vicinity of an existing building resting upon pile foundation to demonstrate the performance of the proposed model in complex simulation environments. Copyright © 2014 John Wiley & Sons, Ltd.     

基于PSO-MLSSVR的土体参数反演方法在深基坑工程中的应用

为了确保基坑工程安全,常常会采用数值模拟的方法预测支护结构的位移,其中岩土体力学参数的选取对于结果的影响最大.本文使用了一种粒子群(PSO)算法结合多输出最小二乘支持向量回归机(MLSSVR)的基坑土体参数位移反分析法,以深圳某深基坑的支护桩顶水平位移监测数据为依据,基于正交设计生成具有代表性的土体参数组合,通过有限元...     

基于MatDEM的砂土侧限压缩试验离散元模拟研究

离散元法基于非连续介质力学理论,尤其适用于砂土等离散介质体的数值模拟研究。利用岩土体离散元模拟软件MatDEM的二次开发功能,研发了砂土侧限压缩试验三维离散元模拟器。对三个不同级配砂土试样进行了侧限压缩试验,并且进行了与之相对应的数值模拟,通过分析对比试验结果与数值模拟结果验证了所开发模拟器的有效性。模拟结果表明:离散元法可以很好地反应砂土压缩过程中的配位数变化;每个模拟样品中,粒径较小的单元受到较大的平均压力,导致平均位移较大;数值计算结果的主要误差是由离散元颗粒自身的泊松比引起的。研究突破了常规土力学研究方法的局限性,为今后岩土工程离散元模拟研究提供了参考。      

Risk-based characterisation of an urban building site

The evaluation of the underground soil stratigraphy is a key aspect in geotechnical site characterisation. However, these means of site exploration are only pinholing subsoil conditions and expert knowledge is needed to understand subsoil conditions in order to build a reliable geological-geotechnical model. This contribution employs a geostatistical simulation methodology for the simulation of random fields representing geological uncertainty. This combines borehole data and expert knowledge via a mathematical framework. Moreover a risk-based site characterisation scheme is developed for urban site characterisation. This novel characterisation scheme offers additional insight into the effects of large-scale, geological spatial variability by using fragility curves to quantify these effects.     

基于Hardin曲线的土体边界面本构模型在ADINA软件中的实现

为了提出一种适合于岩土地震数值模拟的土体本构模型,基于土体动应力-应变关系的Hardin曲线及其在非等幅往返荷载下的Pyke修正,采用von Mises准则在偏应力平面上构造边界面,以反向加载点和当前应力点的连线在边界面上投影的比例作为硬化参数,推导了塑性硬化模量并给出该边界面本构的具体增量表述。在有限元软件ADINA中通过自定义材料的二次开发实现了该本构模型,并利用动三轴试验对该本构模型进行了验证。数值模拟与试验结果的对比表明,本构模型能如实反映土体的应力-应变关系。针对实际工程场地的地震反应,应用边界面本构模型在ADINA中进行了二维数值模拟,与SHAKE91的计算结果进行了对比,说明了该本构模型应用于岩土地震工程问题的合理性。     

挤压地层双护盾TBM围岩变形及应力场特征研究

为研究挤压地层双护盾隧道掘进机（TBM）作用下围岩变形及应力场特征,采用FLAC3D建立了完整模型,并详细阐述了隧道掘进机（TBM）施工过程中的模拟方法,重点分析了隧洞纵横断面内围岩位移场、应力场、塑性区特征。模拟结果表明,两腰下部范围内的围岩与TBM护盾发生接触并产生挤压,拱顶并未接触;受刀盘与护盾连接处的尺寸高差和前后护盾的锥度影响导致仰拱围岩内出现3次加卸载,仰拱内部环向应力和径向应力均大于拱顶和两腰,而且其主应力方向与径向线斜交,受扰动剧烈,但仰拱下方70°范围内的围岩基本处于弹性状态;横向断面内围岩塑性区自上而下逐渐减小,且距掌子面越远塑性区范围越大,但后盾塑性区范围变化不大。