共查询到20条相似文献，搜索用时 78 毫秒
1.
J. A. M. Teunissen 《国际地质力学数值与分析法杂志》2007,31(1):23-51
This paper evaluates the mechanical behaviour of yielding frictional geomaterials. The general Double Shearing model describes this behaviour. Non‐coaxiality of stress and plastic strain increments for plane strain conditions forms an important part of this model. The model is based on a micro‐mechanical and macro‐mechanical formulation. The stress–dilatancy theory in the model combines the mechanical behaviour on both scales. It is shown that the general Double Shearing formulation comprises other Double Shearing models. These models differ in the relation between the mobilized friction and dilatancy and in non‐coaxiality. In order to describe reversible and irreversible deformations the general Double Shearing model is extended with elasticity. The failure of soil masses is controlled by shear mechanisms. These shear mechanisms are determined by the conditions along the shear band. The shear stress ratio of a shear band depends on the orientation of the stress in the shear band. There is a difference between the peak strength and the residual strength in the shear band. While peak stress depends on strength properties only, the residual strength depends upon the yield conditions and the plastic deformation mechanisms and is generally considerably lower than the maximum strength. It is shown that non‐coaxial models give non‐unique solutions for the shear stress ratio on the shear band. The Double Shearing model is applied to various failure problems of soils such as the direct simple shear test, the biaxial test, infinite slopes, interfaces and for the calculation of the undrained shear strength. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
2.
Numerous constitutive models built on coaxial theory and validated under axi‐symmetric condition often describe the stress–stain relationships and predict the inceptions of shear banding in sands inaccurately under true triaxial condition. By adopting an elaborated Mohr–Coulomb yield function and using non‐coaxial non‐associated flow rule, a 3D non‐coaxial elasto‐plasticity model is proposed and validated by a series of true triaxial tests on loose sands. The bifurcation analysis of true triaxial tests on dense sands predicts the influence of the intermediate principal stress ratio on the onset of shear band accurately. The failure of soils is shown to be related to the formation of shear band under most intermediate principal stress ratio conditions except for those which are close to the axi‐symmetric compression condition. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
3.
In this paper, the non‐coaxial relation between the principal plastic strain increments and the principal stresses, which results from the internal friction in geomaterials, is analyzed, and the phenomenon of the unbalanced development of plastic flow in two conjugate directions is discussed. A non‐coaxial, unbalanced plastic flow model for Coulomb frictional materials is developed and used to determine the orientation of shear band in geomaterials. It is shown that the unbalanced index r of plastic flow has important effect on the orientation of the shear band, and the orientation determined by the conventional plastic flow theory is only a special case of the proposed model when r=0. This result soundly explains the reason that the geomaterials with the same internal friction angle and dilatancy angle can have very different shear band orientations. In addition, the difference between the intrinsic and apparent dilatancy angles is analyzed, and it is emphasized that the dilatancy angle commonly used in practice is indeed the apparent dilatancy angle. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
4.
A constitutive model for the simulation of non‐coaxiality, an aspect of anisotropic behavior of sand subjected to the rotation of the principal stress axes, is presented in this paper. Experimental studies have shown that non‐coaxiality or non‐coincidence of principal plastic strain increments with principal stress axes under loadings involving the rotation of principal stress axes may be considerable. Besides, the rotation of the principal stress axes results in dramatic effects on stiffness and dilatant behavior of sand. Therefore, the consequences of principal stress axes rotation on deformational behavior, dilatancy and soil stiffness must be taken into account in theoretical and practical problems. To this aim, the following steps are taken: (1) A general relationship for flow direction with respect to possibility of non‐coaxial flow is developed. Moreover, special circumstances linking non‐coaxiality to instantaneous interaction between loading and soil fabric are proposed. (2) Proposing novel expressions for plastic modulus and dilatancy function, the model is enforced to provide realistic simulations when sand is subjected to the rotation of the principal stress axes. Finally, with numerous examples and comparisons, the model capabilities are shown under various stress paths and drainage conditions. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
5.
6.
We investigate the quasi‐static simple shear flow of a two‐dimensional assembly of cohesionless particles using discrete element method (DEM) simulations. We focus on the unsteady flow regime where the solid would experience significant evolution of stresses, mobilised shear strength and dilation. We construct the DEM model using a discretised‐wall confined granular cell where the apparent boundary is allowed to dilate or contract synchronously with the confined solid. A rather uniform simple shear field is achieved across the whole assembly, which benefits rheological studies in generalising constitutive laws for continuum methods. We examine two aspects of the simple shear behaviour: macroscopic stress and strain rate evolution, particularly the non‐coaxiality between the principal directions of the two; and micromechanics such as evolution of fabric. For an initially anisotropic specimen sheared under constant normal pressure, the direction of principal stress rotates towards that of the principal strain rate, gradually reducing the degree of non‐coaxiality from about 45° to fluctuating around 0°. The rate in approaching coaxiality is slower in samples with larger initial porosity, stress ratio and mean stress. Generally, a faster rate in approaching coaxiality in simple shear is observed in a more dilatant sample, which often shows a larger degree of mobilised fabric anisotropy, suggesting the possible important role of instantaneous internal friction angle. The evolution of principal fabric direction resembles that of the principal stress direction. © 2013 The Authors. International Journal for Numerical and Analytical Methods in Geomechanics published by John Wiley & Sons, Ltd. 相似文献
7.
Euripides Papamichos 《国际地质力学数值与分析法杂志》2010,34(6):581-603
Borehole failure under anisotropic stresses in a sandstone is analyze numerically for various borehole sizes using a nonlinear elastic–plastic constitutive model for a Cosserat continuum. Borehole failure is identified as macroscopic failure of the borehole through the development of shear bands and breakouts. The results compare well both qualitatively and quantitatively with experimental results from polyaxial tests on Red Wildmoor sandstone. They show that the hole size effect of the borehole failure strength is independent of the far‐field stress anisotropy and follows a ? power law of the hole size. A similar scale effect equation with a ? power law is proposed for the scale effect of the maximum plastic shear strain at failure. This equation can be useful for better predicting hole‐size‐dependent failure with standard codes based on classical continua. The effect of stress anisotropy on the borehole failure stress is found to be independent of the hole size. The failure stress decreases linearly to 40% as the stress anisotropy increases. However, the maximum plastic shear strain at failure is stress anisotropy independent and therefore the critical plastic shear strain for failure is only hole‐size dependent. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
8.
The yield vertex non‐coaxial theory is implemented into a critical state soil model, CASM (Int. J. Numer. Anal. Meth. Geomech. 1998; 22 :621–653) to investigate the non‐coaxial influences on the stress–strain simulations of real soil behaviour in the presence of principal stress rotations. The CASM is a unified clay and sand model, developed based on the soil critical state concept and the state parameter concept. Without loss of simplicity, it is capable of simulating the behaviour of sands and clays within a wide range of densities. The non‐coaxial CASM is employed to simulate the simple shear responses of Erksak sand and Weald clay under different densities and initial stress states. Dependence of the soil behaviour on the Lode angle and different plastic flow rules in the deviatoric plane are also considered in the study of non‐coaxial influences. All the predictions indicate that the use of the non‐coaxial model makes the orientations of the principal stress and the principal strain rate different during the early stage of shearing, and they approach the same ultimate values with an increase in loading. These ultimate orientations are dependent on the density of soils, and independent of their initial stress states. The use of the non‐coaxial model also softens the shear stress evolutions, compared with the coaxial model. It is also found that the ultimate shear strengths by using the coaxial and non‐coaxial models are dependent on the plastic flow rules in the deviatoric plane. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
9.
In this work, the interface behavior between an infinite extended narrow granular layer and a rough surface of rigid body is investigated numerically, using finite element method in the updated Lagrangian (UL) frame. In this regard, the elasto‐plastic micro‐polar (Cosserat) continuum approach is employed to remove the limitations caused by strain‐softening of materials in the classical continuum. The mechanical properties of cohesionless granular soil are described with Lade's model enhanced by polar terms, including Cosserat rotations, curvatures, and couple stresses. Furthermore, the mean grain diameter as the internal length is incorporated into the constitutive relations accordingly. Here, the evolution and location of shear band, within the granular layer in contact with the rigid body, are mainly focused. In this regard, particular attention is paid to the effects of homogeneous distribution and periodic fluctuation of micro‐polar boundary conditions, prescribed along the interface. Correspondingly, the effects of pressure level, mean grain diameter, and stratified soil are also considered. The finite element results demonstrate that the location and evolution of shear band in the granular soil layer are strongly affected by the non‐uniform micro‐polar boundary conditions, prescribed along the interface. It is found that the shear band is located closer to the boundary with less restriction of grain rotations. Furthermore, the predicted thickness of shear band is larger for higher rotation resistance of soil grains along the interface, larger mean grain diameter, and higher vertical pressure. Regarding the stratified soil, comprising a thin layer with slightly different initial void ratio, the shear band moves towards the layer with initially higher void ratio. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
10.
Wei Wu 《国际地质力学数值与分析法杂志》2000,24(3):245-263
The paper investigates the incipience of shear band with an incrementally non‐linear constitutive equation. Necessary conditions for the emergence of shear band are derived. The lower bound solution is obtained by taking the strain rate inside and outside the shear band into consideration. Numerical results of localized bifurcation for general stress and strain are presented and compared with experiments. In the principal stress space, the stresses at the onset of shear band form a surface, which is partially enclosed by the failure surface for homogeneous straining. The significance of the analysis for identification of the material parameters and verification of the constitutive model against experiments is discussed. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
11.
It has been known that classical continuum mechanics laws fail to describe strain localization in granular materials due to the mathematical ill‐posedness and mesh dependency. Therefore, a non‐local theory with internal length scales is needed to overcome such problems. The micropolar and high‐order gradient theories can be considered as good examples to characterize the strain localization in granular materials. The fact that internal length scales are needed requires micromechanical models or laws; however, the classical constitutive models can be enhanced through the stress invariants to incorporate the Micropolar effects. In this paper, Lade's single hardening model is enhanced to account for the couple stress and Cosserat rotation and the internal length scales are incorporated accordingly. The enhanced Lade's model and its material properties are discussed in detail; then the finite element formulations in the Updated Lagrangian Frame (UL) are used. The finite element formulations were implemented into a user element subroutine for ABAQUS (UEL) and the solution method is discussed in the companion paper. The model was found to predict the strain localization in granular materials with low dependency on the finite element mesh size. The shear band was found to reflect on a certain angle when it hit a rigid boundary. Applications for the model on plane strain specimens tested in the laboratory are discussed in the companion paper. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
12.
Geotechnical experiments show that Lode angle‐dependent constitutive formulations are appropriate to describe the failure of geomaterials. In the present study, we have adopted one such class of failure criteria along with a versatile constitutive relationship to theoretically analyze the effects of Lode angle on localized shear deformation or shear band formation in loose sand for both drained and undrained conditions. We determine the variation in the possible stress states for shear localization due to the introduction of Lode angle by considering the localized deformation as a bifurcation problem. Further, similar bifurcation analysis is performed for the stress states along a specific loading path, namely, plane strain compression at the constitutive level. In addition, the plane strain compression tests have been simulated as a boundary value finite element problem to see how Lode angle affects the post‐localization response. Results show that the inclusion of a Lode angle parameter within the failure criterion has considerable effects on the onset, plastic strain, and propagation of shear localization in loose sand specimens. For drained condition, we notice early inception of shear localization and multiple band formation when the Lode angle‐dependent failure criterion is used. Undrained localization characteristics, however, found to be independent of Lode angle consideration. 相似文献
13.
This paper deals with FE investigations of shear localization in dilatant granular bodies. The calculations were carried out with a hypoplastic constitutive law enhanced by micro‐polar terms to properly model the shear zone evolution. The behaviour of an initially medium dense sand specimen with very smooth and very rough horizontal boundaries was analyzed during a plane strain compression test. A stochastic distribution of the initial void ratio was assumed to be spatially correlated. Attention was focused on the non‐coaxiality of the directions of the principal strain increments and principal stresses in the shear zone and on the stress–dilatancy rule. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
14.
重力坝深层滑动主要表现为沿缓倾角的软弱结构面形成滑移通道,滑移通道内应变积聚且应变梯度急剧不连续,是典型的应变局部化现象,采用经典连续介质理论进行数值模拟时存在病态的有限元网格依赖性。引入Cosserat连续体理论作为正则化机制,提出了基于Cosserat理论的Mohr-Coulomb弹塑性模型,考虑非关联的流动法则,在经典塑性理论框架下采用向后Euler隐式积分算法进行应力更新。采用ABAQUS的自定义单元接口（UEL）进行二次开发,进行了平面应变条件下单轴受压的数值验证。数值模拟结果表明,该模型能保证应变局部化问题的正定性。基于Cosserat理论的重力坝深层抗滑稳定分析结果表明,采用超载法进行重力坝渐进破坏过程模拟时,基于经典连续体理论的模拟结果有较大的网格依赖性,而且结果偏于安全,而采用Cosserat连续体理论的结果对网格密度不敏感。 相似文献
15.
The paper reports on the results of theoretical and experimental investigations on the spontaneous formation of shear bands in sand bodies. The phenomenon is considered as a bifurcation problem. Consequently, material response and configuration-dependent loading determine the bifurcation mode. Both Coulomb's and Roscoe's solutions of inclination of the shear band can be correct theoretically and experimentally. The first one holds for non-rotating stress axes, the second one for co-rotating stress and strain increment axes during failure. Values in between can occur if the rotation of principal stress axes is not equal to one of these limits. If Coulomb's inclination of shear band occurs, there is a thin deforming material layer separating rigid bodies. Inside the shear band non-coaxiality of strain increment and stress holds from the beginning. If Roscoe's inclination of shear band occurs, it is separating two deforming bodies. Inside the shear band strain increment and stress are coaxial at peak. 相似文献
16.
与 Cosserat 理论相比,偶应力理论在一定程度上可以降低数值框架的复杂度,已逐渐应用于岩土体应变局部化分析中。然而,一般的偶应力有限元法需要满足 C1连续性,即单元内部和单元交界面上的应变都需要具有连续性。为了避免开发较为复杂的C1型偶应力单元,在 Cosserat 连续体理论框架下,通过借助罚函数方法对 C1连续性进行松弛来获得偶应力理论的逼近解,建立了基于罚函数的偶应力有限元方法 PCS-FEM。通过平面应变条件下的弹性圆孔应力集中问题对 PCS-FEM 方法的有效性进行了验证,并应用于土体应变局部化分析中。通过对Ottawa砂的平面应变试验进行数值模拟,发现 PCS-FEM 方法获得的应力−应变曲线及剪切带破坏形态与试验结果基本一致,且能够克服经典连续体理论病态的网格敏感性问题,保证应变局部化问题的适定性;通过对承受偏心荷载作用下的土坡应变局部化经典算例进行分析,发现 PCS-FEM 方法同样可以克服土坡应变软化阶段的网格敏感性问题,展现土体的渐进破坏过程。 相似文献
17.
Numerical studies of shear banding in interface shear tests using a new strain calculation method 总被引：1，自引：0，他引：1
Strain localization is closely associated with the stress–strain behaviour of an interphase system subject to quasi‐static direct interface shear, especially after peak stress state is reached. This behaviour is important because it is closely related to deformations experienced by geotechnical composite structures. This paper presents a study using two‐dimensional discrete element method (DEM) simulations on the strain localization of an idealized interphase system composed of densely packed spherical particles in contact with rough manufactured surfaces. The manufactured surface is made up of regular or irregular triangular asperities with varying slopes. A new simple method of strain calculation is used in this study to generate strain field inside a simulated direct interface shear box. This method accounts for particle rotation and captures strain localization features at high resolution. Results show that strain localization begins with the onset of non‐linear stress–strain behaviour. A distinct but discontinuous shear band emerges above the rough surface just before the peak stress state, which becomes more expansive and coherent with post‐peak strain softening. It is found that the shear bands developed by surfaces with smaller roughness are much thinner than those developed by surfaces with greater roughness. The maximum thickness of the intense shear zone is observed to be about 8–10 median particle diameters. The shear band orientations, which are mainly dominated by the rough boundary surface, are parallel with the zero extension direction, which are horizontally oriented. Published in 2007 by John Wiley & Sons, Ltd. 相似文献
18.
This paper investigates the effects of a non‐coaxial model on simulated stress–strain behaviour of granular materials subject to simple shearing under various initial conditions. In most cases, a significant difference of predictions between coaxial and non‐coaxial modelling is found during the early stage in shearing. With the increase in shearing, non‐coaxial simulations approach and tend to coincide with coaxial simulations. It is also found that the roles of non‐coaxial modelling in simulating simple shear behaviour are considerably influenced by hardening rules, flow rules, initial static lateral pressure coefficients. In some cases, the non‐coaxial modelling gives a similar simulation as the coaxial modelling. In other cases, the non‐coaxial modelling decreases the hardening response or softening response of materials, compared with the coaxial modelling. Under certain conditions, the predicted peak strength of materials with non‐coaxial modelling is larger than that for coaxial modelling. Some of these observations can be attributed to the amount of principal stress rotation in various cases analysed. Others can be attributed to the difference between the directions of the non‐coaxial plastic flow and those for coaxial plastic flow. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
19.
We consider discontinuous bifurcations as the indicator of a localized failure for a class of composites that are characterized by elastic fibres reinforcing an elastic–plastic matrix. A macroscopic tangent stiffness tensor for the fibre‐reinforced composite is developed by consistently homogenizing the contribution of fibres in a spherical representative volume element. Analytical solutions are derived for the critical hardening modulus and corresponding bifurcation directions for the case of plane strain loading. Properties of the solutions are further illustrated on the example of the non‐associated Drucker–Prager model at onset of yielding. Results show that presence of fibres decreases the critical hardening modulus, thus inhibiting the onset of strain localization. The rate of decrease in the critical hardening modulus is the highest for pure shear, followed by uniaxial tension, uniaxial compression, biaxial tension and biaxial compression. The main fibre parameters that control the onset of strain localization are their volumetric content and their stiffness modulus whereby very stiff fibres can produce the most significant decrease in the critical hardening modulus, especially for the state of biaxial tension. The critical hardening modulus for the non‐associated Drucker–Prager model exhibits a full range of localization modes including compaction bands, dilation bands, and transition in the form of shear bands regardless of the presence of fibres. Presence of fibres affects bifurcation directions, except in the case when Poisson's ratio of the matrix is equal to 0.25. The results demonstrate stabilizing effects of fibres by which they provide the control against the onset of strain localization. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
20.
考虑塑性应变率梯度的单轴压缩岩样轴向响应 总被引：9，自引：7，他引：2
基于剪切应变率梯度格式,采用解析方式研究了岩石材料在单轴压缩条件下的应变软化的结构响应。根据非局部连续介质模型,提出了一维二阶剪切应变率梯度格式。非局部剪切应变率与局部剪切应变率及其二阶梯度有关。将经典塑性理论中的局部剪切应变率替换为非局部剪切应变率,可以直接得到局部剪切应变率的封闭解析解,而不必通过将局部剪切应变对时间求导获得。通过对局部剪切应变率积分,得到了沿剪切带方向的相对剪切速度。试件峰值强度后的端部速度由弹性及塑性两部分构成。前一部分由虎克定律描述;后一部分与相对剪切速度有关。对弹性及塑性两部分速度求和,得到了单轴压缩岩样剪切破坏问题轴向响应的解析式。研究表明：试样高度越大、内部长度越小、剪切软化模量越大及泊松比越小,则岩样的轴向响应倾向于脆性。根据岩样与矿柱的相似性,岩样响应倾向于脆性,意味着矿柱将失去稳定性,发生矿柱岩爆。目前的基于剪切应变率梯度格式的主要优点是简洁。 相似文献