共查询到20条相似文献，搜索用时 484 毫秒
1.
Block theory has been widely applied to stability analysis of rock engineering due to its clear concept and elegant geometrical theory. For a general block with multiple discontinuity planes, it is assumed that contact is only maintained on a single plane (single-plane sliding) or two intersecting planes (double-plane sliding) in block theory analysis. Since the normal forces and shear resistances acting on the other discontinuity planes are omitted, it can cause unreasonable estimations of block failure modes and incorrect calculation of factors of safety. In this study, a new method is presented that permits to consider the contribution of the normal forces and shear resistances acting on each discontinuity plane to the block stability. The proposed method meets all of the force-equilibrium and moment-equilibrium conditions and provides a rigorous solution for stability of general blocks with any number of faces and any shape. Some typical polyhedral blocks in rock slopes are analyzed using block theory and the proposed method. The results indicate that the traditional block theory may give a misleading conclusion for the predictions of stability and sliding direction of rock blocks when contact occurs on more than two discontinuity planes. 相似文献
2.
A Boundary Element based Discontinuous Deformation Analysis (BE‐DDA) method is developed by implementing the improved dual reciprocity boundary element method into the open close iterations based DDA. This newly developed BE‐DDA is capable of simulating both the deformation and movement of blocks in a blocky system. Based on geometry updating, it adopts an incremental dynamic formulation taking into consideration initial stresses and dealing with external concentrated and contact forces conveniently. The boundaries of each block in the discrete blocky system are discretized with boundary elements while the domain of each block is divided into internal cells only for the integration of the domain integral of the initial stress term. The contact forces among blocks are treated as concentrated forces and the open–close iterations are applied to ensure the computational accuracy of block interactions. In the current method, an implicit time integration scheme is adopted for numerical stability. Three examples are used to show the effectiveness of the algorithm in simulating block movement, sliding, deformation and interaction of blocks. At last, block toppling and tunnel stability examples are conducted to demonstrate that the BE‐DDA is applicable for simulation of blocky systems. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
3.
深部岩体具有块状层次结构,深部动载造成岩块发生相互间的振动脱离产生低摩擦效应,从而极易诱发原先处于平衡状态的岩体的动力变形破坏。在前人研究基础上,将块系岩体振动简化为等效质量-黏弹性模型,引入岩石摩擦滑移速率弱化模式,最终得到块系岩体滑移失稳计算模型。通过计算分析块系岩体自身特性及外荷载特性对岩块间低摩擦效应的影响。理论计算表明：水平静力及外扰动保持不变,增大岩块间弹性系数或者减小黏性系数,更容易引发岩体低摩擦滑移。随着冲击扰动、水平拉力幅值的增加,岩块的水平残余位移量值增加,当它们幅值超过一临界值时,岩块发生自持续滑移失稳运动。冲击扰动诱发岩块间不可逆位移、动力滑移失稳的临界能量与剪切力水平密切相关,在较大的剪切内力条件下,极其微弱的动力扰动即可诱发较大的岩块间不可逆位移甚至岩块的动力滑移失稳,随着剪切内力的减小,诱发岩块滑移失稳的能量阈值不断增大,当剪切内力低于岩块动摩擦强度时,单次冲击扰动只能诱发岩块间的不可逆位移。初步开展扰动诱发含初应力紫砂岩块体滑移试验,试验结果与理论计算基本符合,证明该模型的可行性。 相似文献
4.
F. Tonon 《国际地质力学数值与分析法杂志》2007,31(14):1567-1608
After describing the kinematics of a generic rigid block subjected to large rotations and displacements, the Udwadia's General Principle of Mechanics is applied to the dynamics of a rigid block with frictional constraints to show that the reaction forces and moments are indeterminate. Thus, the paper presents an incremental‐iterative algorithm for analysing general failure modes of rock blocks subject to generic forces, including non‐conservative forces such as water forces. Consistent stiffness matrices have been developed that fully exploit the quadratic convergence of the adopted Newton–Raphson iterative scheme. The algorithm takes into account large block displacements and rotations, which together with non‐conservative forces make the stiffness matrix non‐symmetric. Also included in the algorithm are in situ stress and fracture dilatancy, which introduces non‐symmetric rank‐one modifications to the stiffness matrix. Progressive failure is captured by the algorithm, which has proven capable of detecting numerically challenging failure modes, such as rotations about only one point. Failure modes may originate from a limit point or from dynamic instability (divergence or flutter); equilibrium paths emanating from bifurcation points are followed by the algorithm. The algorithm identifies both static and dynamic failure modes. The calculation of the factor of safety comes with no overhead. Examples show the equilibrium path of a rock block that undergoes slumping failure must first pass through a bifurcation point, unless the block is laterally constrained. Rock blocks subjected to water forces (or other non‐conservative forces) may undergo flutter failure before reaching a limit point. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
5.
《Geomechanics and Geoengineering》2013,8(4):293-302
This paper presents a comparative study of two methods, Sarma's method and the discontinuous deformation analysis (DDA), for rock slope stability analysis. The comparison concerns the stability analysis of two classic rock slopes. The study shows that the DDA, which accounts for the block kinematics, provides a very different factor of safety as compared with Sarma's method. More realistic reaction forces around each rock block can be obtained by the DDA, including the thrust forces between rock blocks and the forces between the base and the blocks. The DDA's result shows two possible directions for the relative movement between two contiguous blocks at the initiation of slope failure. It also indicates that the limit equilibrium condition may not occur along the interfaces of rock blocks at the initiation of slope failure. The determination of realistic interaction forces around each block will be very important in rock slope stability analysis if nonlinear failure criteria are considered. 相似文献
6.
An original theoretical model has been devised to simulate mass flow over hill slopes due to gravitational sliding. The sliding mass is discretized into a sequence of contiguous blocks which are subjected to gravitational forces, to bottom friction and to surface resistance stresses that are generally negligible for subaerial flows, but are relevant for submarine slides. The blocks interact with each other while sliding down the hill flanks because of internal forces that dissipate mechanical energy and produce a momentum exchange between the individual blocks, yet conserving the total momentum of the mass. Internal forces are expressed in terms of interaction coefficients depending on the instantaneous distance between the block centers of mass, which is a measure of the deformation experienced by the blocks: the functional dependence includes three parameters, namely the interaction intensity ¯, the deformability parameter and the shape parameter , by means of which a wide range of interaction types can be fully accounted for. The time integration is performed numerically by solving the equations for the block velocities and positions at any time ti by means of the block accelerations at the previous time ti-1, and by subsequently updating the block accelerations, which allows to proceed iteratively to the following times. The model has been tested against laboratory results available from literature and by means of several numerical experiments involving a simplified geometry both for the sliding body and the basal surface, with the purpose of clarifying the influence of the model parameters on the slide dynamics. The model improves the performance of the existing kinematic models for slides, moreover preserving an equivalent numerical simplicity. Future applications and possible improvements of this model are suggested. 相似文献
7.
Summary Discontinuities in a rock mass can intersect an excavation surface to form discrete blocks (keyblocks) which can be unstable. Once a potentially unstable block is identified, initial normal and shear stresses on each block face are calculated using elastic theory, and are then modified by discontinuity deformations as the keyblock displaces. The modified stresses are summed into resultant forces to evaluate block stability. Since the resultant forces change with displacement, successive increments of block movement are examined to see whether the block ultimately becomes stable or fails. Calculated keyblock stability increases with larger in situ stress magnitudes, larger lateral stress ratios, and larger shear strengths. Discontinuity stiffness controls block displacement more strongly than it does stability itself. Large keyblocks are less stable than small ones, and stability increases as blocks become more slender. 相似文献
8.
经典块体理论在进行块体可移动性分析和稳定性计算时是分步独立进行的,默认不可移动块体是稳定的,对不可移动块体不再进行稳定性分析。实际工程中,某些阻碍岩石块体运动的“岩桥”可能发生破坏,使得不可移动块体发生滑动。当块体的规模比较大时,在重力作用下破坏局部岩桥的可能性会增高。文中提出一种考虑岩桥破坏的块体稳定性分析方法,不再默认不可移动块体绝对稳定,计算块体稳定系数时,将阻碍块体移动的岩桥的抗剪切力纳入到块体的抗滑力中。该方法首先找出块体所有可能的滑动面和裂隙交线,识别块体沿每一个方向运动时需要破坏的岩桥,基于摩尔-库仑强度准则计算出块体沿每个方向运动时的稳定系数,然后,选取最小的稳定系数作为块体的实际稳定系数,最小稳定系数对应的方向即为块体最有可能的运动方向。 相似文献
9.
土石混合体是由软弱的土颗粒与刚度较大、强度较高的块石混合而成,因此土石混合体边坡失稳时,其滑动面受块石-土形成的细观结构特征影响,导致边坡稳定性分析存在诸多困难。本文基于一典型的土石混合体边坡,建立颗粒离散元数值模型,通过室内试验标定土体细观参数,结合强度折减法分析土石混合体边坡与纯土体边坡破坏模式及滑面形态的差别。结果表明:块石的存在使得土石混合体边坡的接触力链分布极不均匀,块石的存在使接触力链形成闭环,增大了滑动面积,对边坡有加固作用;土石混合体边坡的滑面发展趋势总倾向于穿过块石间的缝隙,滑面的发展与块石分布状况相关联;当块石含量达到20%以上时,边坡稳定性显著提高;块石含量达到30%~50%时,滑坡以局部块体运动为主,无法形成连续的滑动面。研究成果可为土石混合体边坡的变形破坏机理分析提供依据与参考。 相似文献
10.
The paper considers a plane joint or interface element suitable for implementation into a standard non-linear finite element code. The element is intended to model discontinuities with rough contact surfaces, such as rock joints, where dilatant behaviour is present. Of particular concern is the formulation of a constitutive model which fully caters for all possible histories of opening, closing and sliding (accompained by dilation or contraction) in any direction. The non-linear incremental constitutive equations are formulated in a manner appropriate for a back-ward difference discretization in time along the path of loading. The advantage of such an approach is that no essential distinction need be drawn between opening, closing and sliding. Further, a convenient formulation of the constitutive equations is facilitated by representing the different contact conditions in relative displacement space. The state diagram in relative displacement space, however, changes from one time step to the next, and evolution equations for the updating must be formulated. These concepts are illustrated for two rock-joint models: a sawtooth asperity model and a limited dilation model. The models are based on a penalty formulation to enforce the contact constraints, and explicit equations for the tangent stiffness matrix and for the corrector step of the standard Newton–Raphson iterative algorithm are derived. These equations have been implemented as an user element into the finite element code ABAQUS7. Three examples are presented to illustrate the predictions of the formulation. 相似文献
11.
Summary Although wedge and plane sliding stability analyses are well established in the geotechnical literature, certain geologic environments produce blocks which cannot be adequately modelled as either wedges or plane slides. An example is blocks forming in cylindrically folded sedimentary rocks, where the surface of sliding is neither a single plane nor a double plane but is curved. This type of block may be idealized as a prismatic block with multiple sliding planes, all with parallel lines of intersection. If the sliding planes number three or more, the distribution of normal forces, and hence the factor of safety, is indeterminate. A new analytical model for sliding stability analysis is described in which the distribution of normal forces on the contact planes is chosen to minimize the potential energy of the system. The classic wedge and plane solutions are shown to be special cases of this more general model, which allows determination of the safety factor for any shape of prismatic contact surface. An example from Tennessee concerning a block with a curved sliding surface is described and the factor of safety compared with the standard wedge analysis. It is shown that with three or more contact planes, the safety factor may be significantly lower than that calculated from the wedge model, which provides an upper limit on stability. 相似文献
12.
To improve the computational efficiency of the numerical manifold method for discontinuous deformation simulations, a spatial-domain coupled explicit-implicit time integration algorithm is proposed. A subdomain partition algorithm based on a super manifold element is developed for the numerical manifold method to simulate dynamic motions of blocky rock mass. In different subdomains, explicit or implicit time integration method is employed respectively based on its contact and motion status. These subdomains interact through assembling the corresponding explicit or implicit time integration-based matrices of different rock blocks. The computational efficiency of the discontinuity system under dynamic loading is improved by partially diagonalizing the global matrices. Two verification examples of a sliding block along an inclined plane under a horizontal acceleration excitation and a multiblock system acted on by dynamic forces are studied to examine the accuracy of the proposed numerical method, respectively. A highly fractured rock mass situated on an inclined slope subjected to seismic excitations is then studied to show the computational efficiency of the developed algorithm. The simulated results are in good agreement with those from the versions using purely implicit or explicit time integration algorithm for the numerical manifold method. The computational efficiency is shown to be higher using the proposed algorithm, which demonstrates its potential for application in dynamic analysis of highly fractured rock masses. 相似文献
13.
14.
《Journal of Structural Geology》2004,26(6-7):1317-1339
Classifying and assessing geotechnical aspects of rock masses involves combining parameters in various ways, guided by empirical considerations, to derive quantitative geotechnical parameters. Geological structures and the deformation history of rocks underpin the nature of rock masses. The kinematics of a deforming rock mass may occur as sliding along throughgoing discontinuities or as distributed sliding on block faces. Distributed sliding will tend to disrupt the continuity of planar structures such that data on the size and shape of blocks is needed, rather than relying on discontinuity orientation data alone. Orientation and spacing data can be combined to provide a geometric analysis of block systems from linear samples, such as drill core. Dihedral angles and spacing of sequential pairs of discontinuities provides a sample of the size and shape of blocks that can be interpreted stereologically. Further detail can be derived by combining neighbouring intersections that enclose or partially enclose individual blocks. The shape and size of a block can be represented on a stereograph with the enclosing faces shown as poles and their perpendicular distance from an arbitrary point inside the block shown as a number. Identifying the size and shape of specific blocks rather than relying on statistical methods is beneficial to critical aspects of design such as analysing keyblocks that would be exposed during excavations. The detailed characterization of block size and shape is also a step toward interpreting the kinematics of rock mass deformation and the analysis of rock masses as ultra-close packed dilatant granular systems. 相似文献
15.
Discontinuous deformation analysis (DDA) is a numerical approach used to simulate the post-failure behavior of a blocky assembly. Three available algorithms incorporate seismic impacts into DDA simulations for earthquake-induced slope failure. The following methods are used: directly applying time-dependent accelerations to falling/sliding blocks (Method 1); adding time-dependent accelerations to base block (Method 2); and time-dependently constraining seismic displacements of the base block (Method 3). However, incorrect absolute movements of falling/sliding and base blocks were obtained using Method 1. Additionally, relative movements between falling/sliding blocks and the base block are opposite to those simulated by the other two algorithms—Methods 2 and 3. Since locating an earthquake-induced landslide before an earthquake is extremely difficult, the seismic movements of base rock are recorded. Method 1 applies recorded seismic data to sliding blocks in conflict with d’Alembert’s principle of mechanics. Additionally, in Method 2, when the computation time step must be longer than the time in seismic data, computational results reveal abnormal base block displacements due to the non-zero velocity recorded at the end time of seismic data in seismic DDA. In this study, a novel algorithm to diminish the velocity of the base rock in the seismic analysis is utilized to modify Method 2. Furthermore, this work confirms that DDA with the modified Methods 2 and 3 is a practical approach for earthquake-induced landslide simulations. 相似文献
16.
G.H. Shi 《Geomechanics and Geoengineering》2014,9(2):80-96
Continuous computation and limit equilibrium computation are the two independent computations for practical rock engineering. For global stability analysis, limit equilibrium is still the fundamental method. For any numerical method, reaching limit equilibrium requires large displacements, discontinuous contacts, precise friction law, multistep computation and stabilised time-step dynamic computation. Therefore three convergences are unavoidable: convergence of equilibrium equations, convergence of open-close iterations for contacts and convergence of the contact forces of dynamic computations. This paper utilises mainly two dimensional discontinuous deformation analysis (DDA) and an available simple version of three dimensional DDA. The applications show DDA has the ability to reach the limit equilibrium of block systems. 相似文献
17.
An interface constitutive model is presented accounting for slip and sliding effects and also for dilatancy phenomena. The microslip effects are described by considering spherical asperity interaction with variation of contact area and generation of progressive or reverse slip zones. The incremental constitutive equations are derived with proper memory rules accounting for generation and annihilation of particular slip zones during the process of variable loading. It is further assumed that sliding of spherical contacts occurs along large asperities whose slope varies due to the wear process. The predicted shear and dilatancy curves are shown to provide close quantitative simulation of available experimental data. The strain ratchetting effect for non-symmetric cyclic loading was exhibited using the asperity wear model. The model presented could be applied to simulate rock joints, masonry, or concrete cracked interfaces, under monotonic and cyclic loading. 相似文献
18.
P. P. Oreste 《Geotechnical and Geological Engineering》2009,27(1):53-69
The analysis of dowels (non pre-stressed passive reinforcements) for the stabilization of potentially unstable rock blocks
due to sliding is one of the most interesting problems concerning the static of underground rock chambers. The stabilizing
force produced by dowels is not in fact known a priori and it depends on the dowel–rock interaction. The parameters that influence
the problem are not only the geometrical ones of the dowel, but also the mechanical characteristics of the rock and the orientation
of the displacement vector of the block. A detailed study of the dowel–rock interaction has been carried out in this work.
This analysis is able to lead to the evaluation of the axial and shear forces in the dowel in correspondence to the crossed
discontinuity of the block. The vectorial composition of these two forces constitutes the stabilization force of the dowel.
An extensive parametric analysis then made it possible to determine the great variability of the stabilization force of the
dowel in function of the influential parameters of the problem. The graphics that were obtained can be considered a useful
design instrument as they quickly allow the dimensioning of the dowels to be done to reach the required safety factor for
the rock block. 相似文献
19.
在现行条分法的基础上,建立了一种全新的边坡条块力-位移分析方法。现行条分法中条块底边的力学参数采用极限平衡状态力学参数,即每个条块底边均处于极限平衡状态,这种极限平衡状态对于残余应力是较为适宜的。采用理想弹塑性模型和一种全新的本构模型描述土条块底边力学特性,并以不平衡推力法为例,提出了临界状态条块的确定方法,针对可能的破坏模式进行了分析,并提出了边坡整体破坏各条块的应力和位移的确定方法,计算出相应的稳定系数,如传统强度折减法、综合下滑力-抗滑力法、主推力法、综合位移法和富余位移法稳定系数。通过计算分析论证了理想弹塑性模型在不改变力学参数的情况下是难以描述边坡渐近破坏过程的,而一种全新的本构模型可以描述渐近破坏各条块的力学行为。提出的条块力-位移法可以确定边坡在不同荷载和位移条件下的稳定性,也可以获得边坡的推力变化和滑面移动特征,进而实施边坡的应力、位移和稳定性的初步预测预报。 相似文献
20.
裂隙岩体一般块体理论初步 总被引：9，自引：0，他引：9
本文在前人工作的基础上提出了一种通用方法解决任意大小裂隙、任意形状工程岩体的岩石块体识别及稳定性计算问题。在这一方法中裂隙既可以是实测裂隙也可以是通过随机模拟方法生成的随机裂隙,工程岩体可以是任意由多面体组合成的形状,如复杂边坡或地下洞室,而且岩体和裂隙面可以是非均质的。这一通用过程被称为一般块体理论。它克服了关键块体理论(或楔形体法)中存在的弱点:即把裂隙假设为无限大的不连续面,因此无法预测岩石块体的数量和位置,而且一般只能识别出简单形状的凸形体。它在块体识别的同时将复杂的块体分解为几个简单的凸形块体,从而使复杂块体的几何描述、体积及重力计算、力学分析、可视化表示等一系列问题大大简化。在对一公路边坡的裂隙进行了详细观测之后,利用本文的方法进行了块体识别和稳定性分析,对本文的研究进行初步的验证。 相似文献