考虑脆性损伤和渗流的圆形水工隧洞解析解

针对脆性岩石的损伤特性引入脆性损伤模型,导出三维各向同性损伤情况下的损伤变量表达式。在考虑损伤、不同工况主应力对应关系和渗流作用影响下,推导出圆形水工隧洞洞周围岩的弹塑性应力和塑性半径的解析解。对岩石的脆性特性的研究表明,考虑脆性损伤时的洞周围岩塑性半径比不考虑时有一定程度提升;围岩脆性特征越明显,塑性半径越大。无论洞周围岩应力状态如何,考虑渗流影响后,洞周围岩塑性区范围均会有所扩展。在内水压从0逐渐增加过程中,塑性半径随内水压增加而逐渐减小,直至塑性区消失,在超出临界内压后,塑性区重新出现,塑性半径逐渐增大。     

深埋圆形巷道围岩塑性区位移及特征曲线新解和参数分析

考虑中间主应力、脆性软化、剪胀特性、塑性区弹性模量及弹性应变定义等综合影响,提出深埋圆形巷道弹塑性计算的基本假定,进而建立围岩塑性区位移及特征曲线新解,对其进行可比性分析与比较验证,并详细探讨各因素影响特性。所得围岩塑性区位移及特征曲线新解是不需要借助任何数值算法的理论解析解,具有广泛的适用性和很好的可比性。研究结果表明:巷道结构的强度理论效应即中间主应力影响显著,相应的支护可以减弱或改用轻型支护;不考虑脆性软化或剪胀特性都将低估围岩实际位移,工程设计偏危险;半径相关的塑性区弹性模量得到的围岩变形及特征曲线处于上、下限之间,能体现巷道开挖卸载受扰程度的距离渐进性变化;塑性区弹性应变应优先选用更合理和准确的第三定义,且各因素之间存在相互影响。     

考虑岩体扩容和塑性软化的软岩巷道变形解析

考虑岩体的扩容和塑性软化特性,引进扩容梯度和软化模量的概念,推导出均匀介质中软岩巷道应力和变形的理论解答。与其他理论模型进行比较,验证了理论模型的正确性。通过算例分析了岩体的扩容梯度和软化模量对围岩塑性区、破裂区半径以及围岩变形和压力的影响。分析结果表明,考虑岩体的扩容和塑性软化特性使得分析更加准确,研究成果对软岩巷道支护设计与施工具有一定指导意义。     

考虑围岩软化特性和应力释放的圆形隧道黏弹塑性解

将围岩的塑性应变软化特性引入到考虑应力释放的圆形隧道黏弹塑性解中,并且在围岩的软化和残余强度阶段考虑围岩的塑性体积膨胀特性,提出了考虑塑性软化以及塑性体积膨胀和围岩应力释放的圆形隧道弹塑性解。当软化系数k = ∞、膨胀系数h = s时,该解转化为黏弹-脆塑性解;当k = 0、h = s时,则转化为黏弹-理想塑性解,进一步令h = s = 1,则转化为不考虑塑性体积膨胀的黏弹-理想塑性解。通过具体实例计算,分析了掌子面与研究断面间距x、围岩的软化系数k、膨胀系数h和s、支护结构等对围岩塑性区、破碎区半径和变形的影响。当开挖面与研究断面间距x在(0~4)D（D为隧道直径）范围内,随着时间增加塑性圈和破碎区迅速增大;超过4D,塑性区和破碎区半径增量逐渐变小,趋于稳定值;围岩中包含塑性区和破碎区时,二者半径的比值只取决于围岩的性质,与支护结构无关,但支护结构可以限制塑性区及破碎区的范围;考虑应变软化和塑性体积膨胀时,围岩径向位移和塑性区及破碎区半径均大于不考虑应变软化和塑性体积膨胀时的结果;软化系数k增大,围岩位移、塑性区和破碎区半径增加、塑性区半径和破碎区半径之间的比值变小。得到的结果对于隧道工程设计和施工具有一定的指导性和参考价值。     

弹塑性软化模型下隧洞围岩变形与支护压力分析

为深入探究岩体非线性软化应变特性对隧洞围岩力学行为的影响,提出了一种考虑软化系数的岩体弹塑性软化模型,并推导了圆形隧洞围岩应力-应变的求解方法。然后,通过将本文方法与Brown方法、换算后的Mohr-Coulomb屈服准则方法进行对比,验证了本文方法的合理性。最后,探究了软化系数对围岩变形及支护压力影响规律。结果表明：软化系数较大时,采用弹塑性软化模型与弹脆塑性模型算得结果基本相同;对围岩变形及稳定性起主要控制作用的是松动区域范围;控制塑性区域出现时的临界支护压力在数值上与弹塑性区域交界处径向应力是相等的。     

考虑渗流及软化的海底隧道围岩弹塑性分析

在海水渗流的情况下,考虑隧道围岩屈服后强度衰减（塑性软化）的特点,用弹性-塑性软化-塑性残余三线性应力-应变模型推导出了海底隧道围岩塑性软化区半径、塑性残余区半径、洞周边位移、围岩内任一点应力及围岩压力的解析表达式,并与仅考虑渗流及渗流和软化皆不考虑的其他两种情况进行了对比,研究表明：在渗流和材料软化中,后者对隧道围岩稳定最为不利。分析了不同软化程度对结果的影响,该结果为海底隧道围岩稳定性分析等方面提供了理论依据。     

一种求解应变软化岩体中球腔问题的数值方法

均匀地应力场中球形洞室开挖问题可以简化为空间轴对称问题,为分析和评价球形洞室开挖的稳定性,采用有限差分方法对平衡方程与变形协调方程进行离散,将洞室围岩潜在塑性区划分为一系列的球壳,并考虑岩体的应变软化行为,可以得到满足线性Mohr-Coulomb（H-C）准则和非线性Hoek-Brown（H-B）屈服准则的围岩体塑性区范围及相应的应力与变形分布。计算结果表明,随着离散步的增加,数值解逐渐收敛,文中解与已有文献的结果完全吻合,且计算效率更高。对一组强度参数遵循CWFS（黏聚力弱化与摩擦角强化）规律的围岩体内塑性区范围、变形与应力的分布进行了分析,并与常规的应变软化模型的计算结果进行了对比,CWFS模型的结果更合理一些。     

基于广义SMP准则及应力路径法的柱孔扩张分析

柱孔扩张理论广泛用于旁压试验机制分析、沉桩挤土效应等岩土工程问题中。基于经典的Vesic孔扩张理论,采用能够考虑中主应力效应的广义SMP屈服准则,并结合有限应变理论对均质土体中柱孔扩张问题进行分析,根据应力路径假设提出确定塑性区平均体应变的解析步骤。通过大量变参数计算给出具有不同刚度指标 、泊松比 和内摩擦角 的土中柱孔扩张极限状态下的塑性区平均体应变、孔扩张半径比和孔扩张系数。分析结果表明, 、 、 越大,塑性区平均体应变越小,塑性区半径比越大,孔扩张系数越大;随着 增加, 的变化对孔扩张系数的影响更显著。同时,将结果与基于Mohr-Coulomb准则的孔扩张解答进行比较,探讨不同土性参数时中主应力效应对孔扩张系数的影响。文中分析方法及给出的考虑中主应力效应的柱孔扩张系数表可为原位试验分析及桩侧摩阻力估算提供参考。     

基于Hoek-Brown准则的隧洞围岩变形研究

由于在许多实际条件下,比如节理岩体中,线性的M-C准则不太适用,非线性的Hoek-Brown比较适用,因此,可以尝试使用这一非线性屈服准则对洞室变形进行研究。研究隧洞变形时,将围岩分为弹性区、应变软化区、塑性流动区。采用Hoek-Brown准则和非关联流动法则对洞室变形进行了理论推导;软化区域围岩参数随着塑性变形增加而变化,解析法难以求得应力,采用龙格－库塔方法进行数值计算,求解得到塑性软化区和流动区半径,并最终求得洞室变形。通过算例计算表明,在不考虑软化区和流动区时,方法和Carranza-Torres计算结果相差甚小;随着原岩应力的增加,膨胀角对洞室变形的影响增大。     

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

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

单轴压缩非贯通节理岩石尖端塑性区及弹塑性断裂模型研究

基于Drucker-Prager（下简称D-P）准则,建立压缩载荷作用下的非贯通节理岩石的弹塑性断裂模型。针对节理岩石小范围屈服翼裂纹尖端塑性区,推导了D-P屈服准则的纯I、纯II及I、II复合型3种翼裂纹无量纲塑性区径长函数,并与Mises准则的塑性区进行对比;结果表明,D-P准则的I型和复合型塑性区较Mises屈服准则的塑性区大,且其II型及I、II复合型塑性区在翼裂纹上下表面不连续。进一步,引入断裂软化因子以表征节理岩石裂隙断裂扩展后的断裂软化规律,考虑非贯通节理岩石复合型断裂软化,是由于节理尖端翼裂纹应变能密度超过最小应变能密度导致其成核扩展引起的,提出用应变能密度的指数函数形式表征断裂软化变量的演化;塑性屈服函数采用Borja等的应力张量3个不变量的硬化/软化函数,反映塑性内变量及应力状态对硬化函数的影响;建立节理岩石的弹塑性断裂本构关系及其数值算法,并用回映隐式积分算法编制了弹塑性断裂模型的程序。以单轴压缩下非贯通节理岩石为例,分析岩石断裂韧度、节理摩擦系数和节理倾角等参数的影响,结果表明,所提出的弹塑性断裂模型与数值和试验结果比较吻合。     

基于切向应变软化的深埋圆形隧道围岩弹塑性分析

为了研究深埋圆形隧道围岩的应变软化行为,首先通过引入应变软化度,即软化区与残余区交界面处的切向应变和弹塑性交界面处的切向应变的比值,来表征应变软化程度。其次,推导了求解软化区与残余区位移方程,并将该方法求得的位移曲线与Lee法和Cui法对比,验证了方法的正确性。最后,通过算例分析,研究不同力学模型对塑性区半径的影响,以及应变软化度对塑性区中软化区与残余区比例和位移的影响。结果表明：弹塑性交界面处的围压与软化度和支护力均无关,塑性区半径与支护力有关;软化度直接影响软化区与残余区的比例,随着软化度的增加,软化区范围增加,相应的残余区范围减小;软化区与残余区交界面的位移随着软化度的增加而逐渐增加,应变软化度的增加,支护力对其影响增加;洞壁处的位移受支护力的影响较大,随支护力的增加而减小。     

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.     

基于Drucker-Prager准则的岩石弹塑性损伤本构模型研究

大多数岩石材料软化本构模型在硬化函数中引入塑性内变量来表示材料的硬化/软化性质,但并不能反映岩石微裂隙损伤对材料力学性能的影响及单轴拉伸和压缩所表现的初始屈服强度f0与屈服极限fu的差异。基于D-P准则同时考虑塑性软化及损伤软化,建立岩石类材料的弹塑性本构关系及其数值算法。塑性屈服函数采用Borja等的应力张量的硬化/软化函数,反映塑性内变量及应力状态对硬化函数的影响;由于岩石损伤软化是微裂隙扩展所导致的体积膨胀引起的,因此,提出用体积应变表征岩石损伤变量的演化,并用回映隐式积分算法编制了岩石的弹塑性损伤本构程序。对单轴压缩及拉伸荷载作用下的岩石材料试验进行数值模拟,结果表明,所提出的岩石弹塑性损伤本构模型可以较好地符合岩石材料的力学特性。     

饱和土体小孔扩张问题的弹塑性解析解

将土体在圆孔扩张过程中的应力分布分为三个区域，基于修正剑桥模型，推导了圆孔扩张过程中土体在三个区域的超孔隙水压力的解析表达式，对比分析了固结比对土体中应力分布的影响。分析推导可为沉桩、旁压试验、静力触探等岩土工程问题提供理论依据。     

非线性硬化与非线性软化的巷、隧道围岩塑性分析

在符合岩体实际的弹性、应变非线性硬化和具拐点的非线性软化的三段式光滑连接应力-应变关系条件下,求得静水压力下圆形巷、隧道围岩塑性硬化区和软化区的应力分布表达式。证明了所求得的巷道围岩径向和切向应力曲线在硬化区半径 和软化区半径 处均光滑连接,特别是所求得的切向应力 - 曲线没有尖峰向上的应力集中,此结果与日本学者久武胜保的试验结果一致。指出基于理想弹、塑性本构模型的Kastener解有明显缺陷,其 - 曲线尖峰向上的应力集中与实际不符。     

基于能量耗散的球孔扩张理论解答

以空间准滑动面（SMP）准则为基础,推导了扩底桩扩孔压力的理论解。从能量耗散的角度分析球孔扩张的全过程,利用应力不变量推导了符合球孔扩张的屈服准则;化简微分方程得到了弹塑性区应力表达式,进而求出位移、应变表达式;分别利用体积守恒和能量守恒性推导出扩孔压力的表达式。该法考虑了塑性区弹性变形,并得到了扩孔压力p、塑性区半径R与扩孔半径a的关系。算例分析表明,该方法计算的扩孔压力与现场试验得出的结果较好地吻合,塑性区半径和扩孔压力均随扩孔半径的增加而增大,但增幅逐渐减小而趋于稳定值,剪胀角对塑性区半径和扩孔压力影响显著,随着剪胀角的增加,塑性区半径和扩孔压力明显增加。