ANALYSIS OF ROCKBOLT REINFORCEMENT USING BEAM-COLUMN THEORY

A simple analytical procedure that applies classical beam-column theory for evaluating passive rockbolt roof reinforcement is presented in this paper. The analytical model is derived from first principles and is capable of modelling any number of reinforcing bolts. Each rockbolt is modelled as a linear spring and the model allows for non-uniform bolt spacing. In this study the rock beam is assumed to be isotropic and linearly elastic for the sake of simplicity. However, the analytical model can be extended to include anisotropic rockmass as well as inelastic material behaviour. The solution to the coupled set of governing equations is obtained by using a simple numerical solution procedure. The results from the analytical model indicate that the critical buckling load of a rock beam is strongly influenced by the ambient rock modulus. For salt-rock excavations the rock modulus typically declines with time due to various phenomena, and a diminished modulus could seriously compromise roof stability. The other main conclusion of this study is that rockbolts loose their effectiveness in restraining a roof beam once its critical buckling load is approached. In such a situation, increasing bolt stiffness does not improve its reinforcing action on a roof beam but it enhances the possibility of bolt failure due to anchor pull-out. © 1997 by John Wiley & Sons, Ltd. Int. J. Numer. Anal. Meth. Geomech., vol. 21, 241–253 (1997)     

Analytical and Numerical Study of the Mechanics of Rockbolt Reinforcement around Tunnels in Rock Masses

Summary  This paper addresses the problem of quantifying the mechanical contribution of rockbolts installed systematically around tunnels excavated in rock masses. The mechanical contribution referred to here is that of increased stress confinement and decreased tunnel convergences as compared with corresponding stresses and displacements obtained for non-reinforced tunnels. The problem is treated analytically first by presenting a closed-form solution for stress and displacement distributions around a circular tunnel excavated in elastic material and reinforced by grouted or anchored rockbolts. The analytical solution assumes that rockbolts are regularly spaced around the tunnel and that axi-symmetry conditions of geometry and loading apply. The results obtained with the closed-form solution are shown to be equivalent to the results of the same problem solved with traditional numerical methods. Based on the analytical and numerical results and by introducing dimensionless ratios that allow to quantify the increase of radial stresses and the decrease of radial displacements in the reinforced region of the tunnel, the paper shows that reinforcement can have a significant mechanical effect (i.e., increasing the confinement and decreasing the convergences) in tunnels excavated in rock masses of poor to very poor quality. The paper analyzes then the mechanical contribution of rockbolt reinforcement when the rock mass is assumed to behave elasto-plastically. For this case, it is shown that rockbolt reinforcement can also have a critical effect in controlling the extent of the plastic failure zone and the convergences of the tunnel. Correspondence: C. Carranza-Torres, Department of Civil Engineering, University of Minnesota, Duluth Campus, 1305 Ordean Court, Duluth, USA     

软土虹吸排水完整井非稳定流模型及解析解

软土地基的处理一直是工程界的难题,大多数处理方式采用排出软土中水的措施。虹吸排水免动力,且操作简单,可以降低地下水。根据该技术的这些优点,可以将其应用于软土地基处理领域。为了得到虹吸排水作用下土体内部的水位及虹吸管内流量的变化,针对软土中地下水运动进行理论分析。基于Theis理论,针对潜水含水层完整井的轴对称井流微分方程,采用Boltzmann变换,推求得到虹吸过程中软土地基中的水位和流量的显式解析解,并进行室内物理模拟试验验证。结果表明,解析解与试验数据较吻合。与前人经典解析解进行对比,表明相对误差绝对值范围为0～15%,误差较小。推导得到的解析解相比经典解析解而言,计算更简单,更有利于工程应用。     

Layout design of rockbolts for natural ground reinforcement

The present study addresses a layout design of rockbolts for reinforcing natural ground structures applying a special optimization method, called multiphase layout optimization. Rockbolts are used to tighten loosed natural ground, and the layout of rockbolts are determined without sufficient information about the physical properties of the ground materials. Because of this uncertainty, unexpected deformation often occurs at the excavation surface of natural ground. In that case, it is requested to determine an effective layout of the additional rockbolts promptly with respect to the actual deformation at the construction site. However, it is not easy to determine the optimal layout because of its complexity, and consequently, it has no choice but to determine the layout in an empirical way. This study introduces a numerical approach to determine an optimal layout of rockbolts with respect to arbitrarily possible deformation of natural ground. The objective is to maximize the stiffness of the overall ground structure reinforced with rockbolts. For optimization, a gradient‐based optimization scheme is applied because of its numerical efficiency. It was verified from a series of numerical examples that this method has great potential to improve the stiffness of the overall ground structure and shows a certain applicability to a practical design. Copyright © 2013 John Wiley & Sons, Ltd.     

Elastic Solution for Deep Tunnels. Application to Excavation Damage Zone and Rockbolt Support

Summary  A new formulation is presented for deep circular tunnels in rock with cylindrical anisotropy. The formulation is an exact solution since it satisfies equilibrium, strain compatibility, and the anisotropic constitutive model. Complete solutions have been found for two scenarios: tunnel with excavation damage zone, and tunnel with rockbolt support. The solution is based on the assumption of a deep, circular tunnel in a medium with two homogeneous zones: an inner zone surrounding the tunnel, which is either isotropic or anisotropic, and an outer zone, for the remainder of the medium, which is isotropic. Plane strain conditions, elastic response of rock, rockbolts and support, and simultaneous excavation and support installation are also assumed. For tunnels surrounded by an excavation damage zone with reduced rock properties, the tangential stresses and the radial deformations at the tunnel wall are very sensitive to both the magnitude of stiffness reduction of the damaged rock and the size of the damaged zone. The effect of the rockbolts on the rock is approximated by treating the rockbolt-rock composite as a material with cylindrical anisotropy with stiffnesses related to the properties of the rock and rockbolts, and spacing of the rockbolts. Comparisons between the analytical solution and a numerical method show small differences and provide confidence in the approach suggested.     

三种盾构管片配筋设计方法对比及现场实测分析

通过对比修正惯用法的解析解、ANSYS数值模拟解和ANSYS梁-弹簧模型法数值模拟解,结合东莞至惠州城际轨道交通工程DK97+898.238~DK100+897段管片内力现场实测结果,深入探讨各计算方法的合理性和准确性,并得出以下结论:修正惯用法解析解管片内力及配筋最大,修正惯用法数值解次之,梁-弹簧法数值解与现场实测反算值最小。修正惯用法解析解偏保守,容易造成设计配筋过大和浪费;建模较为复杂的梁-弹簧法及现场实测可反馈设计,其结果表明,修正惯用法数值模拟结果比解析解更接近工程实际,本工程内比解析法可节省配筋22.6%左右,且仍有富余量。     

Y型桩桩端阻力产生附加应力的分析计算

Y型沉管灌注桩截面形式复杂,截面几何特性存在4个独立控制变量、8个积分区间。Y型桩桩端均布荷载在地基内部任意点产生附加应力系数的解析表达式求解困难,目前尚无得出简单的解析表达式。借助数学分析软件Mathematica的NIntegrate数值积分功能,以及Geddes推导应力解的思想,得出了Y型桩桩端均布荷载作用在地基内部任意点竖向附加应力系数的数值计算方法,经校核有很高精度。分析了Y型截面4个独立变量外包圆半径R、模板弧度?、开弧间距s、夹角度数? 对附加应力系数的影响及采用Geddes应力解计算产生的误差。分析了随荷载作用面积变化、计算深度变化Y型、三角形、矩形、方形、圆形、圆环、H形、T形截面桩端均布荷载作用在地基内部时截面中心点下附加应力系数的差异及变化规律。将考虑Y型桩截面形状的附加应力计算方法应用于Y型桩单桩的沉降计算,计算结果与静载荷试验数据吻合较好。     

Consolidation of double‐layered ground with vertical drains

Vertical drains are usually installed in subsoil consisting of several layers. Due to the complex nature of the problem, over the past decades, the consolidation properties of multi‐layered ground with vertical drains have been analysed mainly by numerical methods. An analytical solution for consolidation of double‐layered ground with vertical drains under quasi‐equal strain condition is presented in this paper. The main steps for the computation procedure are listed. The convergence of the series solution is discussed. The comparisons between the results obtained by the present analytical method and the existing numerical solutions are described by figures. The orthogonal relation for the system of double‐layered ground with vertical drains is proven. Finally, some consolidation properties of double‐layered ground with vertical drains are analysed. Copyright © 2001 John Wiley & Sons, Ltd.     

Coupling Rheological Behavior of Rock Mass Supported with Rockbolts Based on the Point Load Model

Rockbolts can be regarded as a permanent support structure for underground engineering. However, the coupled rheological mechanism between rockbolts and rock mass remain poorly understood. In this research, the rock mass supported with discretely mechanically or frictionally coupled rockbolts with the point load model was investigated. At first, the elastic solutions of the coupling model were solved. Subsequently, the viscoelastic analytical solutions were acquired by applying the Laplace inverse transforms. Finally, the effect of viscosity coefficients and supporting parameters on the coupling model rheological behavior were explored. The results indicated that the angle of rockbolts support has a greater influence on the radial stress and tangential stress of the rock mass in the elastic state, but has little influence on the rheological state. Moreover, the roof, the interface between the left wall and right wall should be supported strengthen; the viscosity coefficient has little influence on the radial stress and tangential stress in the rheological state and has a negative correlation with the radial displacement. In order to limit the displacement in the rheological state of the rock mass, we choose the rockbolts with a large viscosity coefficient.

     

边坡负压排水非稳定流研究

自启动负压排水方法是一种新兴的间歇性波动排水方法,为对该方法的实现机理进行理论研究,建立了可供解析研究的负压排水模型。先引入狄拉克函数表述负压排水模型中点汇的数理特性,并对此模型推求出潜水剖面二维流运动方程及潜水面边界条件;再使用二维傅里叶正逆变换求解负压排水模型的定解问题,得到负压排水渗流场解析解。设置两种工况分别计算了自启动负压排水方法下数值模拟和理论解得到的边坡浸润线曲线,并进行了对比分析。结果表明：在不同排水孔位置工况的3组条件中,排水孔位置越深,理论解误差越大,其中排水孔位置为（24,10）时,靠近排水孔区域理论解与数值模拟结果误差最大;在排水孔位置固定、不同初始水头高度工况的3组条件中,左、右水头越高,理论解误差越大,其中左水头为28 m、右水头为8 m时,理论解与数值模拟结果误差最大。该自启动负压排水渗流场解析解在简化条件下的浸润线理论解与数值模拟结果吻合度较高。     

An experimental,analytical and numerical study on the thermal efficiency of energy piles in unsaturated soils

Energy piles are bi-functional foundation elements used as structural support as well as ground heat exchangers for shallow geothermal energy systems. Because they are relatively short, energy piles may be partially embedded in unsaturated soils. Saturation conditions influence the thermal properties of the ground and therefore the heat exchange rate, which in turn affects the efficiency of energy piles. This paper combines analytical, experimental and numerical investigations to evaluate the heat exchange rate of energy piles partially or fully embedded in unsaturated soils. The proposed analytical solution is based on the cylindrical heat source theory that treats the soil as a semi-infinite, homogeneous, and isotropic medium. The solution from this theory is multiplied by a function developed analytically in this paper and the outcome is the heat exchange rate for energy piles in unsaturated soils. The proposed function depends on soil saturation, soil and pile thermal properties, and pile geometry. The analytical solution was compared against a finite element solution; which was in turn validated against results from laboratory experiments. Very satisfactory agreements between the analytical, numerical and experimental outputs were observed. The proposed method can be used for a quick and simple evaluation of the efficiency of energy piles in unsaturated conditions. The proposed analytical solution can also be a useful tool for the verification of numerical codes developed for the design of energy piles in unsaturated soils.     

Analytical study of ground responses induced by the excavation of quasirectangular tunnels at shallow depths

The construction of quasirectangular tunnels at shallow depths is becoming increasingly common in urban areas to efficiently utilize underground space and reduce the need for backfilling. To clarify the mechanical mechanism of the stresses and displacements around the tunnels, this study proposes analytical solutions that precisely account for quasirectangular tunnel shapes, the ground surface, the tunnel depth, and the ground's elastic/viscoelastic properties. The Schwarz alternating method combined with complex variable theory is employed to derive the elastic solution, and convergent and highly accurate solutions are obtained by superposing the solutions in the alternating iterations. Based on the solution and the extended corresponding principle for the viscoelastic problem, the time-dependent analytical solutions for the displacement are obtained for the ground assuming any viscoelastic model. The analytical solutions agree well with the finite element method (FEM) numerical results for models that are completely consistent, and qualitatively agree with field data. Furthermore, based on the stress solution combined with the Mohr-Coulomb failure criterion, the predicted initial plastic zone and propagation directions around the tunnels are qualitatively consistent with those determined by the limit analysis. A parametric study is performed to investigate the influences of the rectangular/quasirectangular tunnel shape, burial depth, and supporting pressure on the ground stresses and displacements.     

双重介质中的一类核素迁移数学模型及其反演

探讨了孔隙与单裂隙双重介质中的一类核素迁移数学模型及其反演问题。该核素迁移模型是一个耦合的抛物型方程组定解问题。若已知排污点的核素浓度变化规律,利用Laplace变换及其逆变换方法,求得了核素迁移模型正问题的解析解;反之,由下游裂隙中某个点的实测核素浓度,利用偏微分方程的叠加原理和反问题的拟解法,反求出核素迁移模型反问题的解,即排污点的核素状态。最后,给出核素迁移模型的正问题和反问题的数值模拟。数值结果表明,正问题的解析解能够刻画核素的迁移规律,也显示出所提反问题方法能有效地反演核素污染源。     

Dynamic sliding of tetrahedral wedge: The role of interface friction

The role of interface friction is studied by slow direct shear tests and rapid shaking table experiments in the context of dynamic slope stability analysis in three dimensions. We propose an analytical solution for dynamic, single and double face sliding and use it to validate 3D‐DDA. Single face results are compared with Newmark's solution and double face results are compared with shaking table experiments performed on a concrete tetrahedral wedge model, the interface friction of which is determined by constant velocity and velocity stepping, direct shear tests. A very good agreement between Newmark's method on one hand and our 3D analytical solution and 3D‐DDA on the other is observed for single plane sliding with 3D‐DDA exhibiting high sensitivity to the choice of numerical penalty value. The results of constant and variable velocity direct shear tests reveal that the tested concrete interface exhibits velocity weakening. This is confirmed by shaking table experiments where friction degradation upon multiple cycles of shaking culminated in wedge run out. The measured shaking table results are fitted with our 3D analytical solution to obtain a remarkable linear logarithmic relationship between friction coefficient and sliding velocity that remains valid for five orders of magnitude of sliding velocity. We conclude that the velocity‐dependent friction across rock discontinuities should be integrated into dynamic rock slope analysis to obtain realistic results when strong ground motions are considered. Copyright © 2011 John Wiley & Sons, Ltd.     

An improved numerical approach in surrounding rock incorporating rockbolt effectiveness and seepage force

The solutions of stress and displacement of a circular opening excavated in brittle and strain-softening rock mass incorporating rockbolts effectiveness and seepage force are presented in this study. The evolution equation is reconstructed for the strength parameters that incorporate these factors. Based on the evolution equation, an improved numerical method and stepwise procedure are presented which are compatible with the Mohr–Coulomb (M–C) and the generalized Hoek–Brown (H–B) failure criteria, respectively. Then given three interaction mechanisms between rockbolts and surrounding rock, solutions for stress and displacement are proposed in line with the improved numerical method and numerical stepwise procedure. The proposed approach can be reduced to Fahimifar and Soroush’s (Tunn Undergr Space Technol 20:333–343, 2005) solutions for special cases. The proposed method was validated by field monitoring data and FLAC results of Yanzidong tunnel. Examples under the M–C and generalized H–B failure criteria for rock mass are generated through MATLAB programming. Moreover, parametric studies are conducted to highlight the influence of rockbolts effectiveness in combination with seepage force on the stress and displacement of very good, average, and very poor surrounding rock. Results show that in this case, stress confinement is higher and tunnel convergences are lower than the corresponding stresses and displacements obtained in non-reinforced tunnels. Displacement and plastic radius are also higher than those without considering seepage force.     

Isogeometric analysis for unsaturated flow problems

Unsaturated flow problems in porous media often described by Richards’ equation are of great importance in many engineering applications. In this contribution, we propose a new numerical flow approach based on isogeometric analysis (IGA) for modeling the unsaturated flow problems. The non-uniform rational B-spline (NURBS) basis is utilized for spatial discretization whereas the stable implicit backward Euler method for time discretization. The nonlinear Richards’ equation is iteratively solved with the aid of the Newton–Raphson scheme. Owing to some desirable features of an efficient numerical flow approach, major advantages of the present formulation involve: (a) numerical oscillation at the wetting front can be avoided or facilitated, simply by using either an h-refinement or a lumped mass matrix technique; (b) higher-order exactness can be obtained due to the nature of the IGA features; (c) the approach is straightforward to implement and it does not need any transformation, e.g., Kirchhoff transformation or filter algorithm; and (d) in contrast to the Picard iteration scheme, which forms linear convergences, the proposed approach can however yield quadratic convergences by using the Newton–Raphson method for solving resultant nonlinear equations. Numerical model validation is analyzed by solving a three-dimensional unsaturated flow problem in soil, and its derived results are verified against analytical solutions. Numerical applications are then studied by considering three extensive examples with simple and complex configurations to further show the accuracy and applicability of the present IGA.     

层状地基中群桩水平振动

在考虑地基土分层的基础上,采用动力Winkler地基模型模拟桩土相互作用并运用传递矩阵,求解层状地基中的单桩和群桩的阻抗函数.在计算动力相互作用因子时考虑了被动桩与土的相互作用.最后将相互作用因子和群桩阻抗的本文解与精确解进行对比,验证了本文方法的有效性.     

Two-dimensional meshfree modelling of contaminant transport through saturated porous media using RPIM

This paper presents a new numerical tool to model the two-dimensional contaminant transport through saturated porous media using a meshfree method, called radial point interpolation method (RPIM) with polynomial reproduction. In RPIM, an approximate solution is constructed entirely in terms of a set of nodes and no characterisation of the interrelationship of the nodes is needed. The advection–dispersion equation with sorption is considered to illustrate the applicability of the RPIM. The Galerkin weak form of the governing equation is formulated using 2D meshfree shape functions constructed using thin plate spline radial basis functions. MATLAB code is developed to obtain the numerical solution. Three numerical examples are presented and the results are compared with those obtained from the finite element method and analytical solutions. In order to test the practical applicability and performance of the RPIM, two case studies of contaminant transport through landfill liners are presented. A good agreement is obtained between the results of the RPIM and the field investigation data.     

A finite volume / discontinuous Galerkin method for the advective Cahn–Hilliard equation with degenerate mobility on porous domains stemming from micro-CT imaging

A numerical method is formulated for the solution of the advective Cahn–Hilliard (CH) equation with constant and degenerate mobility in three-dimensional porous media with non-vanishing velocity on the exterior boundary. The CH equation describes phase separation of an immiscible binary mixture at constant temperature in the presence of a conservation constraint and dissipation of free energy. Porous media / pore-scale problems specifically entail images of rocks in which the solid matrix and pore spaces are fully resolved. The interior penalty discontinuous Galerkin method is used for the spatial discretization of the CH equation in mixed form, while a semi-implicit convex–concave splitting is utilized for temporal discretization. The spatial approximation order is arbitrary, while it reduces to a finite volume scheme for the choice of element-wise constants. The resulting nonlinear systems of equations are reduced using the Schur complement and solved via inexact Newton’s method. The numerical scheme is first validated using numerical convergence tests and then applied to a number of fundamental problems for validation and numerical experimentation purposes including the case of degenerate mobility. First-order physical applicability and robustness of the numerical method are shown in a breakthrough scenario on a voxel set obtained from a micro-CT scan of a real sandstone rock sample.