初始应力各向异性状态下圆柱孔扩张机制分析

采用静止土压力系数K0描述土体初始应力各向异性程度,建立了压力控制边界条件的柱孔扩张数值模型,分析压力为边界条件的圆柱孔扩张特性。计算结果表明,初始应力条件各向异性时,柱孔周围的土体的径向位移是不相等的,扩张后的柱孔呈椭圆形,初始应力各向异性是柱孔非对称扩张的内在原因;塑性区的分布具有明显的方向性,塑性区最大半径位于柱孔周围初始大主应力方向上;在相同的扩张压力下,随不排水抗剪强度增加,土体的塑性区减小,且土体初始应力各向异性造成的土体塑性区分布各向异性程度也随之减小;在相同的扩张压力下,柱孔初始应力越大,其塑性区范围越大,且土体初始应力各向异性造成的土体塑性区分布的各向异性程度也越大。     

压力控制的圆孔扩张数值模拟分析

受沉积历史的影响,实际工程中土体的初始应力往往呈各向异性,此时传统圆孔扩张理论的假定条件不再成立,故其适用性也受到限制。借助于FLAC有限差分数值软件,建立了以压力为圆孔扩张边界的二维圆孔扩张模型,从圆孔的变形、孔周土屈服范围、圆孔扩张产生的超静孔水压力分布等方面进行分析,获得在初始应力各向异性的条件下压力控制圆孔扩张过程土体响应规律。计算分析结果表明,初始应力各向异性时,压力控制的圆孔扩张孔口径向位移、塑性区分布、超静孔隙水压力影响范围各个方向不相等;塑性区的分布具有明显的方向性,塑性区最大半径位于孔周土体初始大主应力方向上,并且其值比在相同的扩张压力作用下各向同性初始应力条件下的塑性区半径大,因此传统的初始等应力条件下位移控制的圆孔扩张理论用于分析各向异性初始应力的工程是偏于不安全的。     

K0固结饱和土柱孔扩张问题弹塑性分析

采用K0固结各向异性土体本构模型,将柱孔扩张后周围土体分为弹性区和塑性区,根据柱孔扩张理论和边界条件,推导出K0固结状态下饱和天然土体柱孔扩张问题弹塑性区的应力、塑性区半径以及超孔隙水压力的理论解答。同时,通过算例与修正剑桥模型解答进行对比分析,结果表明,土的不同初始固结状态对柱孔扩张后孔周围的应力和超孔隙水压力产生很大影响,采用考虑K0固结诱发各向异性土体本构模型所得到的应力和超孔隙水压力解答大于修正剑桥模型的解答,但塑性区影响半径却明显小于后者。     

考虑各向异性、剪胀和渗流的柱孔扩张问题统一解

基于统一强度理论,考虑渗流、各向异性、剪胀特性以及中主应力等因素的综合影响,推导了柱孔扩张问题弹塑性解答。结果表明,富水条件下的扩孔问题中考虑渗流体积力的影响是非常有必要的,不考虑渗流影响的扩孔压力比实际工程偏小;渗流体积力越大,塑性区半径越小。另外,传统的扩孔解答假定土体为各向同性初始应力状态高估了极限扩孔压力,也低估了扩孔过程中圆孔周边的塑性区范围,造成工程安全等级偏低。该结果为海底圆形隧道以及渗流条件下地下水平构筑物施工弹塑性分析提供了理论依据,对工程设计具有一定参考价值。     

考虑边界效应的非饱和土圆柱孔扩张问题解析

圆柱孔扩张理论是分析沉桩挤土及旁压试验等岩土工程问题的有效方法,然而已有研究中多假定土体为无限连续介质,当圆柱孔周围存在约束边界时并不适用。采用临界状态模型描述非饱和土的应力-应变关系,在此基础上引入反映土体比体积-吸力关系的液相本构方程,推导了非饱和土中存在固定边界时,圆柱孔扩张问题在不同排水条件下的半解析解。通过计算不同工况,对比分析边界效应、非饱和吸力效应和排水条件对圆柱孔扩张响应的影响。结果表明,受边界效应影响,扩孔压力随着孔径增大而持续增大,这与无边界时扩孔压力趋于某一稳定值的情况有所不同;当圆柱孔扩张相同孔径时,周围固定边界范围越小,所需扩孔压力越大。土体初始吸力越大,扩孔压力也越大,相同应力增量引起的体积压缩量减小,说明土体出现吸力硬化效应。在孔周相同位置处,不排水比排水条件下的土体饱和度高、吸力小,对应的有效应力分量也小;距孔壁越近,该差异越明显。     

考虑中主应力时土体劈裂灌浆力学机制的大变形分析

在塑性力学和大变形理论的基础上分析土体在劈裂灌浆初始阶段的力学机制,将劈裂灌浆的初始阶段视为无限土体中的圆孔扩张问题,并将圆孔周围土体中的应力分布分为两个区域。在弹性区中土体服从小变形假设,在塑性区中服从大变形假设。假设初始劈裂灌浆压力是圆孔扩张到极限时极限扩孔压力;当土体在灌浆压力作用下,大小主应力正好换位时,则出现劈裂时的灌浆压力就为二次劈裂灌浆压力值,推导出塑性区半径、劈裂灌浆的竖向劈裂灌浆和水平向劈裂灌浆压力的理论解答,同时也获得了弹塑性区中的应力场分布规律。该理论的计算结果与工程实测结果比较接近,初步证实了该理论的可靠性。     

大应变静力触探数值模拟及锥形因子影响因素分析

以模拟圆锥在不排水黏土中的静力贯入和分析其影响因素为目标,假设土体为均质弹性—完全塑性材料且服从Mises屈服准则,采用任意拉格朗日—欧拉(ALE)网格划分技术确保锥尖土体在大应变条件下的网格质量,进行大应变有限元数值模拟,并分析了稳定状态下土体刚度指数、原位应力状态和锥尖粗糙程度对塑性区半径与锥形因子的影响,获得了锥形因子表达式。模拟结果表明:塑性区随着刚度指数的增大而增大,锥尖周围塑性区的径向扩张处于柱形孔扩张和球形孔扩张之间,更接近于球形孔扩张;锥形因子随土体刚度指数、锥尖粗糙程度的增大而增大,随土体原位应力状态参数的增大而减小;得到的锥形因子表达式可以量化土体刚度指数、原位应力状态和锥尖粗糙程度的影响,具有较高的精确度。     

复杂应力条件下土体各向异性及其建模思路

由于加荷方式不同,土体在复杂应力状态下在各主应力方向上应力-应变关系表现出显著应力各向异性,在常规三轴试验基础上,采用经典弹塑性理论各向同性土体模型对此不能合理描述。通过真三轴试验,总结应力各向异性柔度矩阵规律,结合试验规律进行相应理论研究,用非线性各向异性弹性矩阵代替弹塑性模型的弹性矩阵,用具有各向异性屈服准则的弹塑性模型描述塑性部分,建立非线性各向异性弹性-塑性模型,可以改善柔度矩阵矩阵形态,反映复杂应力状态下土体应力各向异性特征。     

饱和土中劈裂灌浆压力研究

在塑性力学和大变形理论的基础上,分析土体在劈裂灌浆初始阶段的力学机理。将劈裂灌浆的初始阶段视为无限土体中的圆孔扩张问题,并将圆孔周围土体中的应力分布分为两个区域。在弹性区中土体服从小变形假设,在塑性区中服从大变形假设。利用水力压裂理论和圆孔扩张理论,假设孔隙水压力大于或等于初始应力与扩张应力增量之和时,土体将出现拉应力,若拉应力超过土的抗拉强度,土体将产生开裂,此时的灌浆压力就为初始劈裂灌浆压力。当土体在灌浆压力作用下,大小主应力正好换位时,则出现劈裂时的灌浆压力就为二次劈裂灌浆压力值,推导出饱和土体劈裂灌浆压力的理论公式。该理论的计算结果与工程实测结果比较接近,初步证实了该理论的可靠性。     

基于SMP准则柱孔扩张问题相似解

以柱孔扩张理论为基础,将柱孔扩张后塑性区半径b作为时间刻度,引入扩孔速度概念,将塑性区内任意一点的应力、应变描述成塑性区半径b和柱孔扩张后扩张体单元距孔心距离r的函数,采用平面应变条件下无黏性土的SMP准则和不相关联流动法则,考虑到无黏性土的剪胀性,推导出柱孔扩张后极限扩张力 的相似解。通过参数分析表明,土体强度和剪胀性对塑性区半径和极限扩张力产生很大的影响,即随着土体强度和剪胀性的增加,极限扩张力明显增加,相应塑性区半径却随之减小。通过与基于摩尔-库仑准则解答进行对比分析表明,考虑中间主应力影响的SMP准则所得到的极限扩张力明显高于基于摩尔-库仑准则的解答。该解答为无黏性土中进行静力触探（CPT）试验和静压桩等岩土工程问题提供必要的理论依据。     

Undrained expansion of a cylindrical cavity in clays with fabric anisotropy: theoretical solution

This paper presents a novel, exact, semi-analytical solution for the quasi-static undrained expansion of a cylindrical cavity in soft soils with fabric anisotropy. This is the first theoretical solution of the undrained expansion of a cylindrical cavity under plane strain conditions for soft soils with anisotropic behaviour of plastic nature. The solution is rigorously developed in detail, introducing a new stress invariant to deal with the soil fabric. The semi-analytical solution requires numerical evaluation of a system of six first-order ordinary differential equations. The results agree with finite element analyses and show the influence of anisotropic plastic behaviour. The effective stresses at critical state are constant, and they may be analytically related to the undrained shear strength. The initial vertical cross-anisotropy caused by soil deposition changes towards a radial cross-anisotropy after cavity expansion. The analysis of the stress paths shows that proper modelling of anisotropic plastic behaviour involves modelling not only the initial fabric anisotropy but also its evolution with plastic straining.     

Analytical solution of undrained cylindrical cavity expansion in saturated soil under anisotropic initial stress

This paper presents an analytical solution for undrained elasto-plastic cylindrical cavity expansion in saturated soil under anisotropic initial stress. The problem is formulated by assuming small-strain deformation in the elastic zone and large-strain deformation in the plastic zone. Plastic yielding is determined by the Tresca failure criterion and an associated flow rule. Two stress functions are used to describe the stress state in the two zones around the cavity. The elasto-plastic boundary can subsequently be determined by solving the two stress functions under the stress boundary conditions. Additionally, the cavity pressure-expansion relationship, the total stress and the excess pore pressure around the cavity wall under anisotropic initial stress can be easily obtained by application of a numerical integration. The results show that the cavity pressure and excess pore pressure under the isotropic initial total stress (K = 1) are larger than those under the anisotropic initial stress (K > 1 and K < 1), which is true at all states of the expansion. The higher value of K develops lower stress and pore pressure around the cavity wall at the ultimate states. However, the stress and excess pore pressure are not sensitive to the value of K. The present solution may be used for analyzing the uplift capacity of plate anchors in soils and Horizontal Directional Drilling (HDD) problems such as the tunneling, and pipeline installation.     

Undrained cavity expansion in anisotropic soils with isotropic and frictional destructuration

This paper proposes a semi-analytical solution of undrained cylindrical cavity expansion in anisotropic soils with both isotropic and frictional destructuration. The rigorous derivation based on the general form of the SANICLAY model with destructuration is provided following a standardized solving procedure, and the features of anisotropy and structuration are then invoked in the cavity expansion solution by adopting the non-associated hierarchical model. Cavity expansion tests in both structured and unstructured clays with various overconsolidation ratio are conducted to investigate the evolutions of effective stresses, excess pore pressure, anisotropic parameters and structuration factors during cylindrical expansion. The results show that the effective stresses at the cavity wall are lower after expansion and the cavity excess pore pressure is oppositely higher in structured clays with slightly smaller plastic regions. The evolutions of anisotropy for structured clays appear to follow similar patterns to unstructured cases, whereas the degree of anisotropy is further developed with gradual loss of inter-particle bonds. Finally, the proposed solution is applied to predict the limit pressure of pressuremeter tests in Bothkennar clay, showing its ability for interpretation of in situ testing data in natural structured clays.

     

SH波作用下各向异性弹性半空间中浅埋圆孔附近的动应力集中研究

研究了稳态水平剪切（SH）波作用下含有圆柱形孔洞的各向异性弹性半空间的动力响应问题,探讨介质的各向异性对孔附近动应力的影响规律。运用复变函数和多极坐标方法在变换空间中构造出各向异性弹性半空间中圆孔的散射波函数,所构造的散射波函数能够自动满足半空间水平表面上应力自由的边界条件和无穷远处的Sommerfeld辐射条件。根据圆孔表面应力自由的边界条件,采用最小二乘法求解所构造的散射波函数的系数。分析了SH波以不同角度自弹性半空间内入射时,不同介质特性、入射波数和埋深对圆孔附近动应力集中系数的影响。数值结果表明,各向异性弹性半空间中孔附近的动应力比各向同性情况有明显增高。     

A rigorous semi‐analytical solution for undrained cylindrical cavity expansion in critical state soils

This paper presents a generalized, rigorous and simple large strain solution for the undrained expansion of a vertical cylindrical cavity in critical state soils using a rate‐based plasticity formulation: the initial stress field is taken as anisotropic, that is with horizontal stresses that differ from the vertical stress, and the soil is assumed to satisfy any two‐invariant constitutive model from the critical state (Cam‐clay) family; no simplifying assumption is made during the mathematical derivation; calculating the effective stresses around the cavity requires the solution of a nonlinear equation by means of the Newton–Raphson method in combination with quadrature. Cavity expansion curves and stress distributions in the soil are then presented for different critical state models (including the modified Cam‐clay model). The solution derived can be useful for estimating the instantaneous response of saturated low‐permeability soils around piles and self‐boring pressuremeters and can serve as trustworthy benchmark for numerical analysis codes. Copyright © 2016 John Wiley & Sons, Ltd.     

Two-dimensional elastoplastic analysis of cylindrical cavity problems in Tresca materials

This paper presents analytical elastic-plastic solutions for static stress loading analysis and quasi-static expansion analysis of a cylindrical cavity in Tresca materials, considering biaxial far-field stresses and shear stresses along the inner cavity wall. The two-dimensional static stress solution is obtained by assuming that the plastic zone is statically determinate and using the complex variable theory in the elastic analysis. A rigorous conformal mapping function is constructed, which predicts that the elastic-plastic boundary is in an elliptic shape under biaxial in situ stresses, and the range of the plastic zone extends with increasing internal shear stresses. The major axis of the elliptical elastic-plastic boundary coincides with the direction of the maximum far-field compression stress. Furthermore, considering the internal shear stresses, an analytical large-strain displacement solution is derived for continuous cavity expansion analysis in a hydrostatic initial stress filed. Based on the derived analytical stress and displacement solutions, the influence of the internal shear stresses on the quasi-static cavity expansion process is studied. It is shown that additional shear stresses could reduce the required normal expansion pressure to a certain degree, which partly explains the great reduction of the axial soil resistance due to rotations in rotating cone penetration tests. In addition, through additionally considering the potential influences of biaxial in situ stresses and shear stresses generated around the borehole during drillings, an improved cavity expansion approach for estimating the maximum allowable mud pressure of horizontal directional drillings (HDDs) in undrained clays is proposed and validated.