考虑锚杆作用的深埋软岩隧道黏弹塑性力学响应解析

深埋软岩隧道围岩表现出显著的塑性软化与剪胀特性,而当下的理论分析很少同时考虑这两点,导致预测结果与隧道实际变形行为存在一定误差.为解决该问题,本文基于Kelvin-Voigt流变模型和Mohr-Coulomb强度准则,考虑了塑性阶段时围岩软化与剪胀特征,并引入了掌子面空间约束效应,建立了深埋软岩隧道黏弹-塑性计算分析模...     

不同排水条件下非饱和土中柱孔扩张问题的解析分析

现有的圆柱孔扩张理论已可为诸如石油工程中井筒稳定性鉴定、 及旁压和圆锥贯入实验分析等提供理论依据, 但在非饱和地基压力注浆, 复合地基处理等实际工程问题中却鲜有应用. 基于弹塑性理论和非饱和土力学原理, 采用统一强度理论, 对非饱和土中柱形小孔扩张问题进行了解析研究. 首先将柱孔周围土体分为弹性区和塑性区, 并考虑在弹性区遵循小应变理论, 在塑性区遵循大应变理论, 同时考虑了中间主应力及粒间吸力对非饱和土体强度的影响. 其次应用有效应力表示的统一强度准则, 在本构关系、几何方程、动量平衡方程等基本方程的基础上, 结合相应的边界条件, 最终获得了不同排水条件下柱孔扩张时周围弹塑性区域内的应力场、应变场、位移场及极限扩孔压力的解析表达式. 通过数值算例和参数分析, 在与现有的饱和及非饱和土中柱孔扩张理论进行退化验证的同时, 分析了吸力、剪胀参数、中主应力效应参数及初始径向有效应力等对弹塑性区域内的应力场、应变场及位移场的影响规律, 验证了本文理论的正确性及有效性, 以期为实际工程问题提供合理的理论依据.      

不同拉压特性的双层厚壁圆筒极限承载力解答

采用统一强度理论并考虑材料拉伸与压缩弹性模量的差异性，建立均匀内压作用下双层厚壁圆筒的应力表达式，获得了其内压相应的弹性极限解答、塑性极限解答，并分析拉压强度比、拉压模量系数、统一强度理论参数、半径比及分层半径对弹性、塑性极限内压的影响规律．研究结果表明：弹性、塑性极限内压随拉压强度比的增加而减小，但随统一强度理论参数、半径比的增加而增大；弹性极限内压随分层半径的增加呈现先增大后减小变化，随拉压模量系数的增加而一直减小；塑性极限内压与拉压模量系数、分层半径无关．应用于实际工程时，可根据所得结果选择合理的壁厚及分层半径，再根据材料特性确定其他参数，以便更加准确地计算结构的受力状况．     

被动桩中土拱效应特征与影响参数研究

采用平面有限元方法,对被动桩在粘性、无粘性土体条件下的土拱形成机理进行了分析。从土拱效应的4个主要方面:土拱形状、桩周土体塑性(拉裂)区的分布、竖向位移等值线及桩后土体残余荷载分担比,探讨了改变桩、土参数,如粘聚力、内摩擦角、弹性模量、泊松比、剪胀角、桩土接触面特性时,土拱效应的变化规律。研究表明,土性不同时,土拱形成过程中土体的受力和破坏形式不同;泊松比、剪胀角、桩土接触特征对土拱效应的影响最明显;强度高、剪胀角大、桩土接触面粗糙或泊松比小的土体,易形成土拱。根据研究结果,指出Tom io Ito排桩塑性绕流土压力理论计算公式尚存不足。     

循环加卸载损伤大理岩的动力学特性

利用MTS 815电液伺服岩石实验系统进行上限应力为80%、85%、90%、95%单轴抗压强度的大理岩单轴压缩循环加卸载实验，每种上限应力条件分别设置20、40、60、80次循环。再利用分离式Hopkinson压杆对损伤岩样进行动力学实验。分析了循环加卸载上限应力及循环次数对大理岩塑性应变的影响，揭示了大理岩动态力学参数和破碎吸收能随损伤变量的演化规律。实验结果表明:塑性应变与循环次数呈正相关，且上限应力越大，塑性应变趋于稳定所需的循环次数也会增大；动态单轴抗压强度、动态弹性模量随损伤变量增加呈指数衰减；破碎吸能占比以损伤变量D=0.343为临界点分为两个阶段，D＜0.343时，破碎吸能占比稳定在10%左右，数值约为13 J，当D＞0.343时破碎吸能占比随损伤变量增加不断增大。研究结果可为岩体工程的设计、施工及支护参数的选取提供参考。     

???????????????????????

对具有不同拉压模量的厚壁球壳,采用双剪统一强度理论推导了其扩张问题的应力及位移的 统一解. 分析了不同模量、不同模型控制参数对厚壁球壳扩张时的扩张压力和应力场的影响. 结果表明：厚壁球壳弹性极限压力、应力场、位移场等均随着模量控制参数、模型参数的变 化而变化,在$\alpha<1$的情况下(即$E^ + < E^ - $),可以明显提高球壳的弹 性极限压力$p_e $; 厚壁球壳塑性极限压力与材料的拉压模量无关,与模型参数$\eta $有关,且随$\eta$的增加,先增大后减小. 因此若采用经典的弹性理论和单一的 模型参数对厚壁球壳进行设计计算,会带来较大的误差.     

基于统一强度理论的非静水压圆形隧道塑性区半径的脆塑性摄动解

基于统一强度理论和弹?脆?塑性模型, 综合考虑围岩强度的中间主应力效应和脆性软化, 采用摄动法建立了非静水压圆形隧道塑性区半径的解析解, 继而探讨所得解析解的适用范围, 并与文献复变函数法、摄动法、数值模拟和总荷载不变法进行对比, 最后分析各因素对隧道塑性区形状和大小的影响规律. 研究结果表明: 所建立的圆形隧道塑性区半径摄动解析解为反映中间主应力效应不同程度的系列解答, 可退化为理想弹塑性模型解答且得到文献中4种方法的正确性和合理性验证, 适用于塑性区完全包围隧道的情况, 对应的隧道塑性区为双轴对称的类椭圆, 具有广泛的理论意义和工程应用价值; 摄动参数对隧道类椭圆形塑性区的大小和长/短轴变化都有明显影响; 隧道塑性区范围随中间主应力效应、围岩峰后强度参数的增加均显著减小, 说明不考虑中间主应力效应的Mohr?Coulomb强度准则解答偏保守, 弹?脆?塑性模型相比理想弹塑性模型更适合隧道塑性分析.      

考虑轴向力和渗透力时软化围岩隧道解析

为了研究轴向应力和渗透力共同作用下软化围岩的应力与位移的变化及分布规律. 基于摩尔-库伦屈服准则及应力-应变软化模型并考虑轴向应力和渗透力的共同作用,将整个塑性区分为有限个同心圆环,以弹塑性交界面处的应力、应变为初始值,并采用微小径向应力增量逐步求出各个圆环上的应力应变及塑性区半径,据此重构了考虑渗透力和轴向力共同作用下软化围岩应力应变特性的逐步求解方法. 利用该方法,推导出软化围岩应力应变的解. 计算结果表明:在考虑轴向应力作用下,塑性区半径和隧道围岩位移都随着渗透力的增加而有所减小;当轴向应力为最小主应力时,渗透力的影响更为显著. 这说明渗透力的存在对于隧道围岩的应力应变分布以及塑性半径和围岩的位移有不可忽略的影响.      

SOLUTION AND ANALYSIS OF CIRCULAR TUNNEL FOR THE STRAIN-SOFTENING ROCK MASSES CONSIDERING THE AXIAL IN SITU STRESS AND SEEPAGE FORCE

为了研究轴向应力和渗透力共同作用下软化围岩的应力与位移的变化及分布规律. 基于摩尔-库伦屈服准则及应力-应变软化模型并考虑轴向应力和渗透力的共同作用,将整个塑性区分为有限个同心圆环,以弹塑性交界面处的应力、应变为初始值,并采用微小径向应力增量逐步求出各个圆环上的应力应变及塑性区半径,据此重构了考虑渗透力和轴向力共同作用下软化围岩应力应变特性的逐步求解方法. 利用该方法,推导出软化围岩应力应变的解. 计算结果表明：在考虑轴向应力作用下,塑性区半径和隧道围岩位移都随着渗透力的增加而有所减小;当轴向应力为最小主应力时,渗透力的影响更为显著. 这说明渗透力的存在对于隧道围岩的应力应变分布以及塑性半径和围岩的位移有不可忽略的影响.     

冻融循环和不均匀冻胀下寒区隧道的塑性解答

季冻区长期冻融循环造成围岩强度和变形性能劣化,使得寒区隧道易进入塑性状态,且寒区隧道围岩呈现以径向冻胀为主的不均匀冻胀.合理考虑寒区隧道冻融循环对围岩性能的劣化以及围岩的不均匀冻胀属性,基于Mohr-Coulomb准则推导了寒区隧道冻胀力、应力与位移的塑性解答,同时给出相应的弹性解答和冻结围岩弹-塑性状态的判定方法,对所得解答进行讨论和对比验证,最后探讨了冻融循环、不均匀冻胀与体积冻胀率对寒区隧道应力分布、塑性区半径、洞壁位移和冻胀力的影响规律.研究表明:本文解答具有广泛的适用性和良好的可比性,并得到文献塑性解答的退化验证;冻胀力、洞壁位移与塑性区半径随冻融次数增加分别增大20.3%、8.44倍、2.16倍,以量化长期冻融循环造成围岩性能的劣化效应;冻结围岩由均匀冻胀转变为不均匀冻胀时冻胀力增大42.8%,但塑性区半径几乎无变化;4种体积冻胀率参数均显著影响冻胀力,尤其是水热迁移系数可使冻胀力增大123.6%.本文结果可为季冻区隧道设计与冻害问题解决提供一定的理论依据.     

Micromechanical study of elastic moduli of loose granular materials

In micromechanics of the elastic behaviour of granular materials, the macro-scale continuum elastic moduli are expressed in terms of micro-scale parameters, such as coordination number (the average number of contacts per particle) and interparticle contact stiffnesses in normal and tangential directions. It is well-known that mean-field theory gives inaccurate micromechanical predictions of the elastic moduli, especially for loose systems with low coordination number. Improved predictions of the moduli are obtained here for loose two-dimensional, isotropic assemblies. This is achieved by determining approximate displacement and rotation fields from the force and moment equilibrium conditions for small sub-assemblies of various sizes. It is assumed that the outer particles of these sub-assemblies move according to the mean field. From the particle displacement and rotation fields thus obtained, approximate elastic moduli are determined. The resulting predictions are compared with the true moduli, as determined from the discrete element method simulations for low coordination numbers and for various values of the tangential stiffness (at fixed value of the normal stiffness). Using this approach, accurate predictions of the moduli are obtained, especially when larger sub-assemblies are considered. As a step towards an analytical formulation of the present approach, it is investigated whether it is possible to replace the local contact stiffness matrices by a suitable average stiffness matrix. It is found that this generally leads to a deterioration of the accuracy of the predictions. Many micromechanical studies predict that the macroscopic bulk modulus is hardly influenced by the value of the tangential stiffness. It is shown here from the discrete element method simulations of hydrostatic compression that for loose systems, the bulk modulus strongly depends on the stiffness ratio for small stiffness ratios.     

STUDY ON DAMAGE BIFURCATION AND INSTABILITY OF ROCK-LIKE MATERIALS

The critical bifurcation orientation and its corresponding hardening modulus for rock-like geomaterials are derived by considering the effect of stiffness degradation and volumetric dilatancy under the assumption of isotropic damage. The dependency of the localized orientation on the degree of damage and initial Poisson's ratio of rock is examined and the bifurcation behavior of the uniaxial compression sample under the plane-stress condition is compared with that under plane-strain condition. It is shown that the localization orientation angle intimately depends on both the initial Poisson's ratio and degree of damage for the rock sample under the uniaxial compression condition. As the initial Poisson's ratio or degree of damage increases, the orientation angle of the plane on which localization tends to be initiated gets to decrease. At the same time, the localization orientation angle of a rock sample under the plane-stress condition is larger than that under the plane-strain condition.     

冻融循环和不均匀冻胀下寒区隧道的塑性解答

季冻区长期冻融循环造成围岩强度和变形性能劣化,使得寒区隧道易进入塑性状态,且寒区隧道围岩呈现以径向冻胀为主的不均匀冻胀.合理考虑寒区隧道冻融循环对围岩性能的劣化以及围岩的不均匀冻胀属性,基于Mohr-Coulomb准则推导了寒区隧道冻胀力、应力与位移的塑性解答,同时给出相应的弹性解答和冻结围岩弹-塑性状态的判定方法,对所得解答进行讨论和对比验证,最后探讨了冻融循环、不均匀冻胀与体积冻胀率对寒区隧道应力分布、塑性区半径、洞壁位移和冻胀力的影响规律.研究表明:本文解答具有广泛的适用性和良好的可比性,并得到文献塑性解答的退化验证;冻胀力、洞壁位移与塑性区半径随冻融次数增加分别增大20.3%、8.44倍、2.16倍,以量化长期冻融循环造成围岩性能的劣化效应;冻结围岩由均匀冻胀转变为不均匀冻胀时冻胀力增大42.8%,但塑性区半径几乎无变化;4种体积冻胀率参数均显著影响冻胀力,尤其是水热迁移系数可使冻胀力增大123.6%.本文结果可为季冻区隧道设计与冻害问题解决提供一定的理论依据.     

2021-07期目录

在陈宗基院士关于地下硐室长期稳定性力学问题研究的基础上,采用Sadovsky院士关于复杂地质岩体的等级构造学说,围绕深部岩体非均匀构造与封闭应力固有的统计力学属性,研究了岩体非均匀变形与封闭应力特性,以及深部硐室围岩的长期稳定性等两个力学问题。给出了岩体非均匀构造与封闭应力的数学表征;根据质量守恒定律得到了计算深部硐室围岩长期变形的一般公式;得到了围岩变形中劈裂扩容变形占主要部分的结论,并且阐明了深部围岩卸荷时更易出现劈裂破坏的原因。给出了劈裂破坏形态的演进序列与扩容位移的计算方法。将围岩松动圈范围、破裂区位置和边壁位移的计算结果与锦屏一级电站厂房现有监测数据进行了对比,两者相当吻合。     

爆破地震波作用下城市隧道结构动力响应的敏感性

选取纵波作为研究对象,采用波函数展开法,推导了无限弹性介质中复合衬砌结构隧道在爆破地震波作用下动力响应问题的解析解,结合南京红山南路隧道群开挖工程开展了隧道结构动应力集中因数的敏感性分析,并通过曲线拟合分别得到了隧道围岩及内衬环向动应力集中因数随隧道结构动力学及几何参数变化的回归方程。分析结果表明:隧道围岩和二次衬砌层的弹性模量及围岩泊松比对考察点处环向动应力的影响较大,而初期支护层弹性模量的影响几乎可以忽略不计;衬砌层厚度变化对围岩环向动应力的影响要远大于内衬,在隧道施工设计时可以通过增加衬砌层的厚度平衡隧道结构的受力状态,但增加初期支护层的效果不如二次衬砌层。     

Tailored heterogeneity increases overall stability regime of composites having a negative-stiffness inclusion

Recent theoretical and experimental results have shown the possibility of enormous increases in composite material overall elastic stiffness, damping, thermal expansion, piezoelectricity, etc., when the composite contains a tuned non-positive-definite (i.e., negative stiffness) constituent. For such composite materials to have practical utility, they must be stable. Recent research has shown they can be, for a limited range of constituent negative stiffness. This research has treated linear elastic composite materials with homogeneous phases, via the energy method and full dynamic stability analyses.In the present work, we first show how to analyze the composites previously treated by the comprehensive but simpler static stability approach, obtaining closed-form results. We then employ this approach to show that permitting heterogeneity of the positive-definite phase can substantially increase the range of constituent negative stiffness while maintaining overall composite stability. We first treat the positive-definite phase heterogeneity as piecewise homogeneous, and then treat it as continuously-varying. In the continuously-varying heterogeneity case, we seek the radially optimal distribution of the elastic moduli in the coatings, under constant coating average moduli constraint, to permit the most negative possible inclusion stiffness while maintaining overall composite stability. This is accomplished for three coating cases: constant bulk modulus but arbitrarily radially-varying shear modulus; constant shear modulus but arbitrarily radially-varying bulk modulus; and both moduli arbitrarily radially varying. We find the optimal coatings to be: a heterogeneous one with shear modulus being a specific continuously decreasing function of radius for the first case; a homogeneous one for the second case; and a heterogeneous one with both moduli being either Dirac-delta or Heaviside-step decreasing functions of radius for the last case (if the coating moduli are unrestricted in magnitude or have upper limits, respectively). The results show a substantial increase in the permissible inclusion negative stiffness range is provided by coating heterogeneity, while maintaining overall composite stability. Such an increased range of constituent negative stiffness provides an enlarged tuning parameter range for the development of novel, high-performance composite materials.     

Analysis of main shock of thrust fault earthquake by catastrophe theory

The relationship between work and energy increment of a thrust fault system with quasi-static deformation can be decomposed into two parts: volume strain energy and deviation stress energy. The relationship between work and energy increment of the deviation stress of a simplified thrust fault system is analyzed based on the catastrophe theory. The research indicates that the characteristics displayed by the fold catastrophe model can appropriately describe the condition of earthquake generation, the evolvement process of main shock of thrust fault earthquake, and some important aftershock properties. The bigger the surrounding press of surrounding rock is, the bigger the maximum principal stress is, the smaller the incidences of the potential thrust fault surface are, and the smaller the ratio between the tangential stiffness of surrounding rock and the slope is, which is at the inflexion point on the softened zone of the fault shearing strength curve. Thus, when earthquake occurrs, the larger the elastic energy releasing amount of surrounding rock is, the bigger the earthquake magnitude is, the larger the half distance of fault dislocation is, and the larger the displacement amplitude of end face of surrounding rock is. Fracturing and expanding the fault rock body and releasing the volume strain energy of surrounding rock during the earthquake can enhance the foregoing effects together.     

Loading rate and confining pressure effect on dilatancy,acoustic emission,and failure characteristics of fissured rock with two pre-existing flaws

Investigating the dilatancy, acoustic emission and failure characteristics of fissured rock are significant to ensure their geotechnical stability. In this paper, the uniaxial and triaxial compression experiments with AE monitoring under different loading rates were carried out on fissured rock specimens with the same geometrical distribution of two pre-existing flaws. The dilatancy and AE activity of these specimens were discussed, and the effects of the confining pressure and loading rate on the mechanical parameters and failure characteristics were analyzed. The results show that the exponential strength criterion is more suitable than the Mohr–Coulomb strength criterion to characterize the strength characteristics of fissured rock. The crack evolution and failure characteristics of fissured rock specimens are more complicated than those of intact rock specimens. The failure characteristics of the fissured rock follow the tensile shear coalescence model, crack branching occurs with increasing the loading rate, and the multi-section coalescence model is verified with increasing the confining pressure. The phenomena of stress drop and yield platform usually occur after the dilatancy onset, the specimen does not fail instantaneously, and the propagation and coalescence of cracks cause a sharp increase in the AE signals, circumferential strain, and volumetric strain.     

A micromechanics-based elastoplastic damage model for granular materials at low confining pressure

This paper is devoted to the formulation of a micromechanics-based constitutive model for granular materials under relatively low confining pressure. The constitutive formulation is performed within the general framework of homogenization for granular materials. However, new rigorous stress localization laws are proposed. Some local constitutive relations are established under the consideration of irreversible thermodynamics. Macroscopic plastic deformation is obtained by considering local plastic sliding in a limit number of families of contact planes. The plastic sliding at each contact plane is described by a non-associated plastic flow rule, taking into account pressure sensitivity and normal dilatancy. Nonlinear elastic deformation related to progressive compaction of contacts is also taken into account. Material softening is described by involving damage process related to degradation of microstructure fabric. The proposed model is applied to some typical granular materials (sands). The numerical predictions are compared with experimental data.     

侧向受拉脆性岩石压缩本构及蠕变失效研究

轴向压缩作用下,脆性岩石侧向应力严重影响岩石力学特性。侧向压应力影响下的轴向压缩岩石力学行为已经得到广泛研究,然而侧向拉应力对轴向压缩岩石力学行为影响研究很少。本文基于脆性岩石翼型裂纹扩展模型中,初始裂纹面法向应力与剪切应力的正负方向为判断依据,研究了侧向拉应力对轴向压缩力学行为的影响。发现恒定的侧向拉应力作用下,轴向压缩应力渐进变化过程中,脆性岩石内部细观初始裂纹面的法向应力只能为压缩应力,不存在拉应力情况。分析了从侧向压应力到拉应力转化过程中,脆性岩石轴向压应力与细观裂纹扩展长度关系、轴向压应力与轴向应变关系、岩石峰值强度、裂纹启裂应力及初始弹性模量的变化规律。并分析了侧向拉应力对岩石蠕变裂纹长度、裂纹速率、轴向应变及应变率演化曲线,以及对蠕变失效时间及稳态蠕变应变率的影响。讨论了侧向拉压应力突变转化以及侧向拉应力分级增大对轴向压缩岩石蠕变演化行为影响。该研究为深部地下工程围岩稳定性评价提供了一定理论依据。