共查询到19条相似文献，搜索用时 187 毫秒
1.
岩土工程中常用的屈服准则多以压缩剪切为其破坏机制,然而硬脆性岩体的脆性破坏包括拉伸破坏、张拉剪切破坏和压缩剪切破坏3类,且随着岩体工程向深部发展,张拉剪切破坏成为了洞壁围岩的主要破坏机制。针对此问题,开展了硬脆性大理岩的室内拉剪试验,分析了大理岩拉剪破坏特征,并结合压剪试验结果,建立了考虑张拉剪切破坏机制和应力状态影响的Mohr-Coulomb准则。研究结果表明,硬脆性大理岩破裂面在拉剪应力状态和低正应力压剪应力状态下均具有张拉剪切破坏特征,高正应力压剪应力状态下则只具有压缩剪切滑移特征;拉剪应力状态下,大理岩破裂面张拉破坏特征明显,无明显剪切痕迹,剪切力固定时,剪切位移随着轴向拉力增加而增加;凝聚力和内摩擦角受应力状态影响,凝聚力随正应力增大先减小后增大,内摩擦角则随正应力的增大而减小;凝聚力、内摩擦角随正应力的变化趋势可分为4段,拉剪段、低压应力段、中压应力段和高压应力段,每段的凝聚力、内摩擦角与正应力皆可认为是线性关系,靠近抗拉强度处,内摩擦角趋近90°,凝聚力趋于无穷大;考虑张拉剪切破坏机制和应力状态影响的Mohr-Coulomb准则曲线分为两部分,可采用二次抛物线进行拟合的拉剪段和考虑凝聚力、内摩擦角随正应力演化的压剪段,由此建立的Mohr-Coulomb准则更全面、精度也更高。 相似文献
2.
建立含有非贯通层面和正交次级节理的逆层岩质边坡FLAC/PFC2D耦合计算模型,进行地震动力破坏过程模拟试验,研究了逆层岩质边坡地震动力破坏机理。试验结果证明,在地震动力破坏过程中,边坡内部层面主要产生剪切破坏,少量张拉破坏集中于逆层边坡顶部位置并且总是发生在坡体已经产生动力失稳之后,因此层面的抗拉强度并不影响逆层边坡的地震动力稳定性。坡顶正交次级节理只能产生张拉破坏,形成宏观的岩层倾倒趋势,而坡底的正交次级节理既会产生张拉破坏,也会产生剪切破坏,破坏面滑动趋势明显。动力响应坡顶放大效应和破坏面发育位置深度导致坡顶岩体的张拉倾倒早于坡底岩体的剪切滑动,与逆层边坡静力倾倒破坏顺序相反。 相似文献
3.
建立含有非贯通层面和正交次级节理的逆层岩质边坡FLAC/PFC2D耦合计算模型,进行地震动力破坏过程模拟试验,研究了逆层岩质边坡地震动力破坏机理。试验结果证明,在地震动力破坏过程中,边坡内部层面主要产生剪切破坏,少量张拉破坏集中于逆层边坡顶部位置并且总是发生在坡体已经产生动力失稳之后,因此层面的抗拉强度并不影响逆层边坡的地震动力稳定性。坡顶正交次级节理只能产生张拉破坏,形成宏观的岩层倾倒趋势,而坡底的正交次级节理既会产生张拉破坏,也会产生剪切破坏,破坏面滑动趋势明显。动力响应坡顶放大效应和破坏面发育位置深度导致坡顶岩体的张拉倾倒早于坡底岩体的剪切滑动,与逆层边坡静力倾倒破坏顺序相反。 相似文献
4.
基于抛物线型D-P准则的岩质边坡稳定性分析 总被引：3，自引：0，他引：3
岩质边坡内通常会存在部分拉剪屈服区,因此,在岩质边坡的稳定性分析中采用能够同时考虑拉剪屈服和压剪屈服的H-B准则相较M-C准则来说更为合适。推导了基于H-B准则的抛物线型D-P准则,克服了H-B准则在数值计算中的困难。针对基于抛物线型D-P准则的有限元强度折减法,证明了折减抛物线型D-P准则材料参数的合理性。结合上述研究成果,分析了茨哈峡水电站右岸泄水边坡在天然状况下和泄水雾化状况下的稳定安全度。计算结果表明,该边坡在天然状况下是稳定的,但在泄水雾化状况下将会发生失稳。因此,需对该边坡采取工程处理措施,提高其稳定安全度,以防泄洪雾化失稳。 相似文献
5.
Drucker-Prager系列屈服准则在稳定分析中的应用研究 总被引：3，自引：0，他引：3
Drucker-Prager（下简称D-P）系列屈服准则作为Mohr-Coulomb（下简称M-C）准则的修正模型在岩土工程中得到了广泛的使用。然而采用不同的D-P系列屈服准则可能会得到差距较大的结果,因此,选用合适的D-P系列屈服准则十分重要。在研究了D-P系列屈服准则与M-C准则之间的对应关系之后,指出D-P系列屈服准则对应的Lode角取值较为有限,难以反映出材料在丰富受力状态下的M-C准则屈服强度,将其取值范围扩大至－30°～30°,并以某边坡和向家坝泄12坝段抗滑稳定分析为例,对如何选用合适的D-P系列屈服准则这一问题进行了研究,提出了具体的选择方法。研究结果表明,只要选用的D-P系列屈服准则对应的Lode角能反映影响边坡和坝基稳定的关键部位的受力状态,就能够得到较好的结果,将D-P系列屈服准则对应的Lode角取值范围扩大至－30°～30°也是必要的。 相似文献
6.
Hoek-Brown准则在岩质边坡稳定分析中的优越性 总被引：1，自引：0，他引：1
介绍了最新的用于估计完整岩石或节理岩体剪切强度的半经验准则Hoek-Brown屈服准则。结合岩质边坡工程并基于强度折减法计算边坡安全系数,与等效Mohr-Coulomb屈服准则数值模拟对比,得出Hoek-Brown屈服准则采用基于应力水平的塑性流动法则,考虑了岩体结构、岩块强度、应力状态等多种因素的影响,能很好地反映岩体的非线性破坏特征及机理,符合节理岩体的变形和破坏特点。 相似文献
7.
土体中的水力劈裂破坏机理存在两种观点：张拉破坏机理和剪切破坏机理。前人已经对其进行了大量地研究,并分别基于两种破坏机理提出了起劈压力的理论计算式,但是对其适用性没有进一步探讨。分别分析三向应力状态下土体水力劈裂的张拉破坏和剪切破坏判定准则,并基于Mohr-Coulomb屈服准则分析两种破坏机理的适用条件。分析结果表明,土体的水力劈裂是张拉破坏还是剪切破坏与小主应力σ3和不排水抗剪强度cu的大小密切相关,当小主应力σ3较大且不排水抗剪强度cu较小时,土体的劈裂多受剪切破坏机理控制;反之,土体的劈裂为张拉破坏。最后通过已有文献中的试验结果验证了所提出的两种破坏机理适用条件的合理性。 相似文献
8.
采用不同Drucker-Prager屈服准则得到的边坡安全系数的转换 总被引：3，自引：0，他引：3
研究了边坡稳定分析的有限元强度折减过程中（平面应变条件下）Drucker-Prager（简称D-P）屈服面的变化特点,从而提出了有限元强度折减法采用不同屈服准则计算得到的边坡安全系数之间的转换关系式。在利用有限元强度折减法计算边坡稳定安全系数时,ANSYS有限元软件采用的岩土材料屈服准则为Mohr-Coulomb（简称M-C）六边形外接圆D-P准则,可先求出外接圆D-P准则条件下的安全系数,然后利用所提出的安全系数转换公式就可直接计算出各D-P准则条件下的安全系数。因此,通过转换就可以在ANSYS程序中实现不同M-C准则,而不需要进行二次开发。采用ANSYS软件通过算例分析,比较了由转换关系式得到的安全系数与实际计算的结果,讨论了转算结果的误差。算例结果表明：通过转换关系式得到的安全系数与计算得到的结果非常接近,具有比较高的计算精度,同时也证明所提出的方法是可行的。 相似文献
9.
Drucker-Prager（以下简称D-P）屈服准则因其形式简单、物理意义明确而得到广泛的应用,然而经过多年的应用和研究,其缺点逐渐显现出来,如拉剪区偏大、不具有应力角效应等,针对这些不足,提出了修正的D-P屈服准则。在拉剪区用圆球面和横截面的组合方式代替原来的圆锥面,从而不改变D-P屈服准则压剪区的形式,在考虑岩石实际的拉伸强度的同时,也满足该屈服准则经过单轴抗拉强度点的条件;为考虑岩石在三轴压缩和三轴拉伸状态下不同的强度特性,即应力角效应,引入了一种角隅模型;为保证屈服准则始终通过岩石的单轴抗压强度点,建议了一种材料参数k值的取值方法。为验证屈服准则修正后的适用性和正确性,将修正屈服准则与真三轴试验结果进行了拟合,同时与Mohr-Coulomb（以下简称M-C）准则进行对比,结果表明,修正屈服准则可以很好地描述试验现象,且比M-C准则更接近真实结果。 相似文献
10.
高地应力等复杂环境下大跨度地下洞室群开挖容易诱发高边墙的脆性拉裂破坏,在Mohr-Coulomb准则基础之上,基于拉应变准则,提出大型地下厂房高边墙破坏的复合拉剪屈服准则。根据弹塑性理论,张拉和剪切屈服分别采用相关联和非相关联的塑性流动法则,详细推导了该准则在有限差分法中的计算格式,并对张拉和剪切屈服面的应力空间进行了分界。利用C语言程序编写的动态链接库文件,实现了该准则在通用软件的科学仿真计算。以高应力地区的拉西瓦水电站地下厂房为工程实例,与FLAC3D自带的拉剪屈服准则相比,基于拉应变准则的计算结果表明,大跨度高边墙拉破坏明显,塑性区以及应力松弛区的面积更小,能合理描述硬岩高边墙的脆性拉裂破坏特征。大型硬岩地下洞室群的计算分析验证了该准则的合理性和实用性 相似文献
11.
基于非线性破坏准则的边坡稳定性极限分析 总被引：5，自引：0，他引：5
上限定理是求解岩土工程问题的有效工具。以上限定理为理论基础,分析边坡的稳定性问题,并考虑了岩土材料破坏准则的非线性特性。在非线性Mohr-Coulomb破坏准则下,采用条分法与极限分析上限法相结合的方法,对边坡稳定性进行分析。通过切线法引入非线性强度参数 、 ,推导了岩土材料在非线性破坏准则下的相关联流动法则,建立功能方程,推导了边坡直线滑裂面、折线滑裂面和光滑曲线滑裂面安全系数F的计算方程。采用数学规划法计算后与工程实际中常用的边坡稳定性分析方法进行对比,并获得安全系数F =1.0时的稳定性系数Ns。与已有的研究成果进行比较分析,结果表明了该方法的正确性及优越性 相似文献
12.
In this study, upper bound finite element (FE) limit analysis is applied to stability problems of slopes using a nonlinear criterion. After formulating the upper bound analysis as the dual form of a second-order cone programming (SOCP) problem, the stress field and corresponding shear strength parameters can be determined iteratively. Thus, the nonlinear failure criterion is represented by the shear strength parameters associated with stress so that the analysis of slope stability using a nonlinear failure criterion can be transformed into the traditional upper bound method with a linear Mohr–Coulomb failure criterion. Comparison with published solutions illustrates the accuracy and feasibility of the proposed method for a simple homogeneous slope stability problem. The proposed approach is also applied to a seismic stability problem for a rockfill dam to study the influence of different failure criterions on the upper bound solutions. The results show that the seismic stability coefficients obtained using two different nonlinear failure criteria are similar but that the convergence differs significantly. 相似文献
13.
将极限分析的上限定理、有限元离散思想、随机规划理论和蒙特卡洛方法这四者结合起来,提出了一种土质边坡可靠度分析的上限数值方法。首先采用三节点有限单元离散土质边坡,然后将土体的抗剪参数设为随机变量,根据上限定理构建同时满足三角形单元的塑性流动约束条件、单元公共边的塑性流动约束条件和单元速度边界条件的机动许可速度场,并根据内功功率等于外功功率条件建立目标函数,构建土质边坡可靠度分析的上限法随机规划模型。采用蒙特卡洛方法求解上限法随机规划模型,同时提出了一种基于上限法速度场的边坡失效风险系数估算方法,该方法特别适用于具有多种失效模式的边坡风险分析。对2个经典算例进行了深入分析,验证了方法的正确性。 相似文献
14.
The presence of a weak interlayer has usually an adverse effect on the slope stability. However, the rotational failure mechanism in the conventional upper bound limit analysis cannot rationally describe the sliding of the failure mass along the weak interlayer. Therefore, a new failure mechanism was proposed in this study to evaluate the stability of slopes with weak interlayers using the upper bound limit analysis and the associated factor of safety was determined by the shear strength reduction technique. The new failure mechanism is comprised of rigid blocks undergoing rotational or translational movements, instead of the rotational movement in the conventional failure mechanism. It has also been extended to the stability analysis of slopes in presence of stabilizing piles and pore water pressures. Case studies were carried out on actual slopes with weak interlayers. The proposed rotational–translational failure mechanism was verified by the shear strength reduction finite element method (SSRFEM). Comparisons demonstrate the reliability of limit analysis method with the proposed rotational–translational failure mechanism for slopes with weak interlayers and therefore it can be used as a simple evaluation method for the engineering design. 相似文献
15.
基于非线性Mohr-Coulomb破坏准则,结合极限分析上限法和拟静力分析法,建立功能方程,推导了地震效应下裂缝边坡的安全系数计算方程。采用数学规划方法,计算了不同参数组合条件下的边坡安全系数值,详细分析了非线性条件下一系列参数对边坡稳定性的影响。研究表明,边坡安全系数随非线性参数和地震效应的增大而减小。对比分析可知,在非线性破坏准则下,裂缝深度较大时,裂缝对边坡稳定性影响显著,且边坡越陡影响越大;当裂缝深度超过某个值后,临界破坏面起始端可能不穿过裂缝最底端,而是从裂缝中间某部位穿过。在地震效应作用下,非线性抗剪强度参数对安全系数影响显著。研究成果进一步完善了裂缝边坡稳定性分析内容,所列图表为边坡的设计与施工提供有益参考。 相似文献
16.
Numerical modelling of rock slides is a versatile approach to understand the failure mechanism and the dynamics of rock slopes. Finite element slope stability analysis of three rock slopes in Garhwal Himalaya, India has been carried out using a two dimensional plane strain approach. Two different modelling techniques have been attempted for this study. Firstly, the slope is represented as a continuum in which the effect of discontinuities is considered by reducing the properties and strength of intact rock to those of rock mass. The equivalent Mohr-Coulomb shear strength parameters of generalised Hoek-Brown (GHB) criterion and modified Mohr-Coulomb (MMC) criterion has been used for this continuum approach. Secondly, a combined continuum-interface numerical method has been attempted in which the discontinuities are represented as interface elements in between the rock walls. Two different joint shear strength models such as Barton-Bandis and Patton’s model are used for the interface elements. Shear strength reduction (SSR) analysis has been carried out using a finite element formulation provided in the PHASE2. For blocky or very blocky rock mass structure combined continuum-interface model is found to be the most suitable one, as this model is capable of simulating the actual field scenario. 相似文献
17.
Polygonal finite elements are gaining an increasing attention in the computational mechanics literature, but their application in rock mechanics is very rare. This paper deals with numerical modeling of rock failure under dynamic loading based on polygonal finite elements. For this end, a damage-viscoplastic constitutive model for rock based on the Mohr-Coulomb criterion with the Rankine criterion as a tensile cutoff is employed and implemented with the polygonal finite element method. Moreover, the mineral mesostructure or rock is described by randomly mapping groups of polygonal elements representing the constituent minerals into a global mesh and assigning these groups with the corresponding mineral material properties. The performance of the polygonal elements is compared with that of the linear and quadratic triangular and bilinear quadrilateral elements in numerical simulations of controlled shear band formation under uniaxial compression and lateral splitting failure in the dogbone tension test. Numerical simulations of uniaxial tension and compression tests as well as dynamic Brazilian disc test under increasing loading rates demonstrate that the present approach predicts the correct failure modes as well as the dynamic increase in strength of rock. 相似文献
18.
A micromechanical model is proposed for studying the stability and failure process of slopes based on the gravity increase method (GIM). In this numerical model the heterogeneity of rock at a mesoscopic level is considered by assuming that the material properties conform to the Weibull distribution. Elastic damage mechanics is a method used for describing the constitutive law of the meso-level element, the finite element method (FEM) is employed as the basic stress analysis tool, and the maximum tensile strain criterion and the Mohr–Coulomb criterion are utilised as the damage threshold. The numerical model is implemented into the Realistic Failure Process Analysis (RFPA) code using finite element programming, and an extended version of RFPA, i.e., RFPA-GIM, is developed to analyse the failure process and stability of slopes. In the numerical modelling with RFPA-GIM, the critical failure surface of slopes is obtained by increasing the gravity gradually but keeping material properties constant. The acoustic emission (AE) event rate is employed as the criterion for slope failure. The salient feature of the RFPA-GIM in stability analysis of slopes is that the critical failure surface as well as the safety factor can be obtained without any presumption for the shape and location of the failure surface. Several numerical tests have been conducted to demonstrate the feasibility of RFPA-GIM. Numerical results agree well with experimental results and those predicted using the FEM strength reduction method and conventional limit equilibrium analysis. Furthermore it is shown that selection of the AE rate as the criterion for slope failure is reasonable and effective. Finally, the RFPA-GIM is applied to several more complex cases, including slopes in jointed rock masses and layered rock formations. The results indicate that the RFPA-GIM is capable of capturing the mechanism of slope failure and has the potential for application in a larger range of geo-engineering. 相似文献
19.
In this paper, a nonlinear numerical technique is developed to calculate the limit load and failure mode of structures obeying an ellipsoid yield criterion by means of the kinematic limit theorem, nonlinear programming theory and displacement-based finite element method. Using an associated flow rule, a general yield criterion expressed by an ellipsoid equation can be directly introduced into the kinematic theorem of limit analysis. The yield surface is not linearized and instead a nonlinear purely kinematic formulation is obtained. The nonlinear formulation has a smaller number of constraints and requires less computational effort than a linear formulation. By applying the finite element method, the kinematic limit analysis with an ellipsoid yield criterion is formulated as a nonlinear mathematical programming problem subject to only a small number of equality constraints. The objective function corresponds to the dissipation power which is to be minimized and an upper bound to the plastic limit load of a structure can then be calculated by solving the minimum optimization problem. An effective, direct iterative algorithm has been developed to solve the resulting nonlinear programming formulation. The calculation is based purely on kinematically admissible velocities. The stress field does not need to be calculated and the failure mode of structures can be obtained. The proposed method can be used to calculate the bearing capacity of clay soils in a direct way. Some examples are given to illustrate the validity and effectiveness of the proposed method. 相似文献