基于边界配点法浅埋隧道开挖引起 周围土体变形的分析

对于浅埋隧道施工引起土层位移的解析解,其边界条件的数学处理上比较困难,因此实际工程的应用中受到限制。为解决这一缺陷,基于弹性力学的Airy应力函数,将半无限平面内土体的应力及位移分布转化为解析函数。然后,采用边界配点的方法,控制地表处的应力边界条件( = 0, = 0)及隧道周边土体的位移边界条件( , )。最后,采用最小二乘法,确定土体应力及位移函数的各项系数,求得浅埋隧道在周边土体产生径向位移的作用下,地表及深层土体位移的半数值半解析解。算例分析表明,该方法可以考虑任意荷载及位移作用的边界条件,能够充分发挥解析法和数值法两者的优点,且计算结果与实测结果较吻合,在预测地表及深层土层位移中具有一定的实用价值,可为进一步采取工程措施控制地层变形提供理论依据。     

Modelling poroelastic hollow cylinder experiments with realistic boundary conditions

A general poroelastic solution for axisymmetrical plane strain problems with time dependent boundary conditions is developed in Laplace domain. Time‐domain results are obtained using numerical inversion of the Laplace transform. Previously published solutions can be considered as special cases of the proposed solution. In particular, we could reproduce numerical results for solid and hollow poroelastic cylinders with suddenly applied load/pressure (Rice and Cleary, Rev. Geophys. Space Phys. 1976; 14 :227; Schmitt, Tait and Spann, Int. J. Rock Mech. Min. Sci. 1993; 30 :1057; Cui and Abousleiman, ASCE J. Eng. Mech. 2001; 127 :391). The new solution is used to model laboratory tests on thick‐walled hollow cylinders of Berea sandstone subjected to intensive pressure drawdown. In the experiments, pressure at the inner boundary of the hollow cylinder is observed to decline exponentially with a decay constant of 3–5 1/s. It is found that solutions with idealized step‐function type inner boundary conditions overestimate the induced tensile radial stresses considerably. Although basic poroelastic phenomena can be modelled properly at long time following a stepwise change in pressure, realistic time varying boundary conditions predict actual rock behaviour better at early time. Experimentally observed axial stresses can be matched but appear to require different values for α and ν than are measured at long time. The proposed solution can be used to calculate the stress and pore pressure distributions around boreholes under infinite/finite boundary conditions. Prospective applications include investigating the effect of gradually changing pore pressure, modelling open‐hole cavity completions, and describing the phenomenon of wellbore collapse (bridging) during oil or gas blowouts. Copyright © 2004 John Wiley & Sons, Ltd.     

A new closed-form solution for circular openings modeled by the Unified Strength Theory and radius-dependent Young’s modulus

A new closed-form solution is presented for the stress and displacement distribution surrounding circular openings with finite external radii that are subject to uniform internal and external pressures under plane strain conditions. The specific solution for a deep circular tunnel in an infinite rock mass is also provided. It is assumed that the rock mass is elastic–brittle–plastic and governed by the Unified Strength Theory (UST). In the plastic zone, the radius-dependent Young’s modulus (RDM) model and a non-associated linear flow rule were adopted to establish the radial displacement solution. The new closed-form solution obtained in this paper is a series of results rather than one specific solution; hence, it is suitable for a wide range of rock masses and engineering structures. The traditional solutions, which are based on the Mohr–Coulomb failure criterion and the Generalized Twin Shear Stress yield criterion, can be categorized as special cases of this proposed solution. This new solution agrees reasonably well with the results of a borehole collapse test, a secondary development numerical simulation and an additional closed-form solution using the generalized non-linear Hoek–Brown failure criterion. Parametric studies were conducted to investigate the effects of intermediate principal stress, RDM and dilatancy on the results. It is shown herein that the effects of intermediate principal stress and dilatancy are significant; the RDM model is recommended as the optimum approach for calculating radial displacement and support pressure.     

The effect of Poisson's ratio and the far-field boundary conditions on the accuracy of finite element calculations

It is shown that a finite element calculation which approximates an ‘infinite medium’; problem by a mesh with finite boundaries will yield greater accuracy when stress boundary conditions are applied on the far-field boundary than is obtainable with displacement boundary conditions. In particular, with Poisson's ratio close to 0.5, the accuracy of the latter model is severely impaired, whereas the stress boundary condition model is unaffected for Poisson's ratio of 0.49 and a reasonable mesh. The eight-node quadratic isoparametric element displays superb accuracy for the axisymmetric thick cylinder with either type of boundary condition.     

Thermal stresses in borehole heat exchangers

The backfilling materials of borehole heat exchangers (BHE), particularly the grout material, must provide a suitable thermal contact and ensure durability to the induced thermal stresses because of the heat loading. In this paper, the thermal stresses that occurred in BHEs because of heat injection or extraction is investigated with an analytical solution of a hollow cylinder model that is adapted for time‐dependent heat loading, the geometry of a BHE, and the thermo‐mechanical properties of surrounding ground conditions. Firstly, the hollow cylinder model is solved with the considered boundary conditions in 2D plane stress. Secondly, the temperature differences at the inner and outer circles of the cylinder are evaluated with the heat line source models for continuous and discontinuous loading to observe the impact of the heat loading schedule. The developed analytical solution for thermal stress investigation is validated with numerical models. It is demonstrated that the analytical solutions agree well with numerical results for two types of BHE configurations (co‐axial and single U‐shaped pipes). Furthermore, the calculated maximum stresses are compared with the tensile strength of grout materials obtained from Brazilian tests. It is predicted that the thermal contraction of the grout, partially constrained by the surrounding rock, generates tensile stresses that may lead to cracking in the BHE. According to the results, the stiffness of rock has a primary role on the developed tensile stresses, and the relationship between the thermal conductivity of the ground and of the grout induces a proportional impact on the magnitude of thermal stresses. Copyright © 2015 John Wiley & Sons, Ltd.     

Analytical solution of the harmonic waves diffraction by a cylindrical lined cavity in poroelastic saturated medium

This paper presents a model for the analysis of plane waves diffraction at a cavity in an infinite homogeneous poroelastic saturated medium, lined by a lining composed of four equal segments. An elastic boundary layer is placed between the cavity lining and the infinite porous medium. The boundary layer is simulated by ‘elastic boundary conditions’ in which the bulk matrix stress is proportional to the relative displacement between the lining and the surrounding medium matrix boundary. In addition, fluid impermeability through the intermediate layer is assumed. For the frequencies, that differ from the pseudoresonanse frequencies, the problem was reduced to the problem of an ideal elastic medium. A closed‐form analytical solution of the problem was obtained using Fourier–Bessel series, the convergence of which was proven. It was shown that the number of series terms required to obtain a desired level of accuracy can be determined in advance. The influence of the medium porosity on the medium dynamic stress concentration was studied. Copyright © 2006 John Wiley & Sons, Ltd.     

Drained cavity expansions in soils of finite radial extent subjected to two boundary conditions

The problem of drained cavity expansion in soils of finite radial extent is investigated. Cylindrical and spherical cavities expanded from zero radius subjected to either constant stress or zero displacement at the finite boundary are considered. The generalised analytical solution procedure presented enables more advanced constitutive models to be implemented than have been possible in previous studies. Results generated for Sydney quartz sand highlight substantial differences between cavity limit pressures for the two boundary conditions and for boundaries of finite and infinite radial extent. This is of significance in accounting for boundary effects when interpreting cone penetration tests conducted in calibration chambers. Copyright © 2011 John Wiley & Sons, Ltd.     

正交各向异性岩体中非圆形水工隧洞的解析解

将围岩和衬砌分别视作均质、连续的线弹性正交各向异性和各向同性体,并充分考虑衬砌的支护滞后效应和隧洞运行时的内水压力作用,运用复变函数方法中的幂级数解法,提出了正交各向异性岩体中任意形状水工隧洞的力学解析方法。以直墙半圆拱形水工隧洞为例,所获得的解析解可精确满足衬砌内边界的应力边界条件以及围岩与衬砌接触面的应力、位移连续条件,同时还将解析结果与ANSYS数值结果对比分析,吻合良好。利用获得的解析解,讨论了围岩开挖面上不同的各向异性程度、不同的弹性对称面角度以及隧洞内不同的水压荷载对衬砌以及围岩上应力和位移分布的影响。     

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     

饱和稀疏裂隙岩体三维水流-传热过程中位移和应力的一种半解析计算方法

针对高放射核废深地质处置库近场环境,建立分布热源作用下饱和裂隙岩体三维水流-传热过程中位移和应力的一种半解析计算方法：采用Goodier热弹性位移势和Laplace变换计算由温度梯度产生的温梯位移和应力;考虑单一裂隙的情况,利用经典弹性力学的Boussinesq解和Cerruti解计算为满足边界条件的约束位移和应力,与温梯位移和应力叠加,可得总体热位移和应力;把裂隙面离散为矩形单元集合,采用极坐标系下的解析法计算包含奇点的单元积分,采用数值法计算与分布热源有关和不含奇点的单元积分。与基于裂隙面法向一维热传导假设的一种解析解对比,结果表明,半解析法与解析法的计算结果基本一致,但由于半解析法考虑岩石的三维热传导,因温度时空分布和演变的不同而导致不同的温梯应力。针对一个假想单裂隙岩体三维水流-传热过程,计算温梯位移和应力、约束位移和应力、总体位移和应力;结果表明,裂隙水流-传热可能对位移和应力的分布和演变有显著影响,距离分布热源较近的岩石因升温膨胀受到约束而出现压应力,而距离分布热源较远的岩石则可能因协调收缩受到约束而出现拉应力。     

Explanation of locking of four‐node plane element by considering it as elastic Dirichlet‐type boundary value problem

The aim of this study is to arrive at a better understanding of the phenomenon of locking of low‐order compatible displacement type of finite elements in particular for the hour‐glass mode of the plane four‐node element and dilative materials. To this end the properties of finite elements are investigated in an analytical way, where a finite element is considered as a plane boundary value problem with prescribed boundary displacement (Dirichlet problem). In this paper for the sake of simplicity the simplest possible linear comparison solid, namely isotropic linear elasticity, is applied, although recognizing fully that for a dilative material elasto‐plasticity would be more realistic. From the study described in this paper it is concluded that locking of the four‐node element is not due to any particular numerical formulation of this compatible finite element since, even the analytical solution suffers from this problem. The locking of this element is not related to incompressibility of the material either as the analytical solution shows locking to occur at a parameter set which differs significantly from the one in case of incompressibility. It is shown that locking is a consequence of the combination of the dilative material behaviour and the compatible displacement type of boundary conditions, which leads to infinite isotropic stresses in the element. These infinite isotropic stresses occur at the limit of uniqueness of the solution, which for this element is shown to occur outside the parameter range of the sufficiency of uniqueness. Copyright © 2000 John Wiley & Sons, Ltd.     

Elastic solutions for a transversely isotropic half-space subjected to a point load

We rederive and present the complete closed-form solutions of the displacements and stresses subjected to a point load in a transversely isotropic elastic half-space. The half-space is bounded by a horizontal surface, and the plane of transverse isotropy of the medium is parallel to the horizontal surface. The solutions are obtained by superposing the solutions of two infinite spaces, one acting a point load in its interior and the other being free loading. The Fourier and Hankel transforms in a cylindrical co-ordinate system are employed for deriving the analytical solutions. These solutions are identical with the Mindlin and Boussinesq solutions if the half-space is homogeneous, linear elastic, and isotropic. Also, the Lekhnitskii solution for a transversely isotropic half-space subjected to a vertical point load on its horizontal surface is one of these solutions. Furthermore, an illustrative example is given to show the effect of degree of rock anisotropy on the vertical surface displacement and vertical stress that are induced by a single vertical concentrated force acting on the surface. The results indicate that the displacement and stress accounted for rock anisotropy are quite different for the displacement and stress calculated from isotropic solutions. © 1998 John Wiley & Sons, Ltd.     

偏心受压厚壁圆筒支护结构的应力分析

厚壁圆筒由于能充分发挥材料的抗压特性而越来越多地应用到基坑支护结构工程之中,传统的设计计算主要是在均布压力作用下厚壁圆筒的拉梅解。为了解决实际工程中经常遇到的偏心受压问题,给出了一个偏心荷载作用下厚壁圆筒的应力解析解,并结合工程实例对圆筒支护结构受力进行简化处理后,将解析解应用于计算偏心受压厚壁圆筒支护结构中的拉应力,特别是讨论了地下水位不均匀变化、地面不均匀堆载和土层不均衡开挖等可能会引起较大的圆筒环向拉应力和在圆筒内、外表面产生竖向裂缝,从而对圆筒支护结构产生不利影响。有利于指导对圆筒截面的配筋和控制基坑的开挖施工。     

冻结黏土空心圆柱试样制样方法

冻土空心圆柱仪是室内研究冻土在车辆荷载等复杂应力路径下力学特性的主要仪器,试验中为实现主应力轴方向在垂直于径向平面内连续旋转,其试样多采用薄壁空心圆柱试样,但由于空心圆柱试样薄壁的特点,增加了制样的困难,成为制约冻土空心圆柱试验研究进展的主要阻碍。为此,借鉴冻土三轴试验的制样方法,基于实验室现有的制样条件,设计了一套重塑冻结黏土空心圆柱试样的制样装置,可提高制样效率,减小制样过程中人为因素的影响;随后,基于该制样装置,提出了一种冻结黏土空心圆柱试样的制样方法,给出了制备重塑冻结黏土的制样步骤。最后,通过对所制备试样的物理力学性质测定,发现采用该制样方法可以制备含水率及干密度均匀、力学性质稳定的重塑冻结黏土试样;同时对两个平行制备的冻结黏土试样的不同主应力轴方向的定向剪切试验结果进行了比较,证明了所提出的制样方法具有良好的可重复性,且所得的冻土力学特性与已有研究成果吻合较好,可用于进一步系统研究复杂条件下冻结黏土的静、动力学特性。     

Semi‐analytical elastostatic analysis of unbounded two‐dimensional domains

Unbounded plane stress and plane strain domains subjected to static loading undergo infinite displacements, even when the zero displacement boundary condition at infinity is enforced. However, the stress and strain fields are well behaved, and are of practical interest. This causes significant difficulty when analysis is attempted using displacement‐based numerical methods, such as the finite‐element method. To circumvent this difficulty problems of this nature are often changed subtly before analysis to limit the displacements to finite values. Such a process is unsatisfactory, as it distorts the solution in some way, and may lead to a stiffness matrix that is nearly singular. In this paper, the semi‐analytical scaled boundary finite‐element method is extended to permit the analysis of such problems without requiring any modification of the problem itself. This is possible because the governing differential equations are solved analytically in the radial direction. The displacement solutions so obtained include an infinite component, but relative motion between any two points in the unbounded domain can be computed accurately. No small arbitrary constants are introduced, no arbitrary truncation of the domain is performed, and no ill‐conditioned matrices are inverted. Copyright © 2002 John Wiley & Sons, Ltd.     

Strength criterion for cross-anisotropic sand under general stress conditions

By incorporating the fabric effect and Lode’s angle dependence into the Mohr–Coulomb failure criterion, a strength criterion for cross-anisotropic sand under general stress conditions was proposed. The obtained criterion has only three material parameters which can be specified by conventional triaxial tests. The formula to calculate the friction angle under any loading direction and intermediate principal stress ratio condition was deduced, and the influence of the degree of the cross-anisotropy was quantified. The friction angles of sand in triaxial, true triaxial, and hollow cylinder torsional shear tests were obtained, and a parametric analysis was used to detect the varying characteristics. The friction angle becomes smaller when the major principal stress changes from perpendicular to parallel to the bedding plane. The loading direction and intermediate principal stress ratio are unrelated in true triaxial tests, and their influences on the friction angle can be well captured by the proposed criterion. In hollow cylinder torsional shear tests with the same internal and external pressures, the loading direction and intermediate principal stress ratio are related. This property results in a lower friction angle in the hollow cylinder torsional shear test than that in the true triaxial test under the same intermediate principal stress ratio condition. By comparing the calculated friction angle with the experimental results under various loading conditions (e.g., triaxial, true triaxial, and hollow cylinder torsional shear test), the proposed criterion was verified to be able to characterize the shear strength of cross-anisotropic sand under general stress conditions.     

温度变化对深基坑钢管支撑温度效应影响分析

在基坑围护工程中,大量使用内支撑结构体系,在支撑体系中采用最多的是钢管支撑,因为其属于金属材料,温度变化对支撑结构的影响不容忽视。基于弹性热力学基本原理推导出空心圆筒热应力公式,同时考虑了围护墙的侧向水平位移,定量分析变温对基坑钢管支撑轴力影响,实测轴力值及数值解验证了计算公式的正确性。结果表明:解析解大于数值解,且解析解大于实测值,偏于保守,但差值不大可作为设计安全储备,便于工程师在今后设计过程中作为参考实例。     

Exact solution for drained spherical cavity expansion in saturated soils of finite radial extent

This study presents an exact solution for drained spherical cavity expansion subjected to constant stress and zero displacement boundary conditions in finite medium. The solution is exact as no simplification is involved in the solution procedure in the plastic region and can be solved as an initial value problem. The effect of finite radial extent on the results of cavity expanding to a certain radius is accounted for through the initial values at the elastic-plastic boundary. A critical state constitutive model with a nonassociated flow rule is adopted. The model parameters from literature are then used in generating the results for cavity expansion in infinite and finite radial extent to highlight the difference. Also, the results are presented in a way that can be used to account for the boundary effect for the interpretation of cone penetration tests conducted in calibration chambers.     

垂直界面附近椭圆形夹杂对SH波的散射与地震动研究

利用复变函数法和Green函数法给出了SH波对垂直界面附近椭圆形夹杂散射问题的解析解答。首先,将待求的半空间模型沿垂直边界分割为区域Ⅰ和区域Ⅱ两个直角域。通过保角映射的方法将区域Ⅰ内椭圆形夹杂的外域映射为单位圆外域,并利用镜像方法构造出两个区域内满足直角域边界条件的散射波场及适用的Green函数;其次,利用界面契合思想,通过在界面处添加附加力系的方法建立起满足界面处位移和应力连续条件的无穷代数方程组,并截断有限项求解;最后,给出了求解地表位移幅值的具体算例。结果表明,入射波数、入射角度、夹杂位置、垂直界面以及材料参数都对地表位移幅值的分布有影响。     

The Applicability of the Ground Response Curve to Tunnelling Problems that Violate Rotational Symmetry

The applicability limits of the closed-form solution to the problem of ground response to tunnelling are sounded out by systematically investigating the effect of deviations from some of the important assumptions underlying the closed-form solution. The ground response curve (GRC) expresses the relationship between tunnel support pressure and the radial displacement of the tunnel boundary on the basis of a rotationally symmetric model. The assumptions underlying rotational symmetry are a circular tunnel, a hydrostatic and uniform initial stress field, an isotropic and homogeneous ground and uniformly distributed support pressure. Deviations from these assumptions generally necessitate potentially time-consuming numerical analyses. The paper revisits the classical problem of tunnel excavation in a linearly elastic, perfectly plastic ground obeying the Mohr–Coulomb yield criterion, and analyses the effects of non-uniformity and anisotropy of the initial stress field and of a non-circular tunnel geometry. The results show that the GRC also provides a reasonably accurate approximation of average tunnel convergence for a wide range of ground conditions that violate rotational symmetry.