Fabric dilatancy and the plasticity modeling of granular media     

B. Muhunthan  N. Sasiharan 《国际地质力学数值与分析法杂志》2012,36(9):1181-1193

A novel conceptual model of the mechanics of sands is developed within an elastic–plastic framework. Central to this model is the realization that volume changes in anisotropic granular materials occur as a result of two fundamentally different mechanisms. The first is purely kinematic, dilative, and is the result of the changes in anisotropic fabric. There is also a second volume change in granular media that occurs as a direct response to changes in stress as in a standard elastic/plastic continuum. The inclusion of the two sources of volume change results in three important datum states. When subjected to isotropic strains, the resulting stress state in granular materials is not isotropic but lies upon the kinematic normal consolidation line. There exists a state at which the fabric‐induced volumetric strain rate becomes equal to the stress‐induced volumetric strain rate making the total plastic volumetric strain rate equal to zero. Granular response changes from contractive to dilative at this phase transformation line. The third datum state is the one in which the stress‐induced volumetric strain rate is zero. The sand, however, continues to dilate at this state with the difference between stress and dilation ratio a constant as predicted by Taylor's stress–dilatancy rule. These predictions are shown in accordance with experimental data from a series of drained tests and undrained on Ottawa sand. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

均质土坡潜在滑面的搜索     

任祥  潘霄  门玉明 《地球科学与环境学报》2002,24(1):46-48

均质土坡潜在滑面的确定 ,是对均质土坡进行稳定性计算的必要前提。传统的以瑞典圆弧法为代表的刚体极限平衡方法 ,虽然能够确定出均质土坡的滑面 ,但工作量较大 ,使用不方便。因此 ,在塑性极限分析方法的基础上 ,提出了一种新的搜索滑面的方法  相似文献   

A Generalized Backward Euler algorithm for the numerical integration of an isotropic hardening elastoplastic model for mechanical and chemical degradation of bonded geomaterials     

Claudio Tamagnini  Riccardo Castellanza  Roberto Nova 《国际地质力学数值与分析法杂志》2002,26(10):963-1004

An extended version of the classical Generalized Backward Euler (GBE) algorithm is proposed for the numerical integration of a three‐invariant isotropic‐hardening elastoplastic model for cemented soils or weak rocks undergoing mechanical and non‐mechanical degradation processes. The restriction to isotropy allows to formulate the return mapping algorithm in the space of principal elastic strains. In this way, an efficient and robust integration scheme is developed which can be applied to relatively complex yield surface and plastic potential functions. Moreover, the proposed algorithm can be linearized in closed form, thus allowing for quadratic convergence in the global Newton iteration. A series of numerical experiments are performed to illustrate the accuracy and convergence properties of the algorithm. Selected results from a finite element analysis of a circular footing on a soft rock layer undergoing chemical weathering are then presented to illustrate the algorithm performance at the boundary value problem level. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

GPS水准的有限元法与多面函数法的加权综合模型   总被引：4，自引：0，他引：4  

刘万林  王利  赵超英 《地球科学与环境学报》2004,26(3):48-51

应用有限元法与多面函数法的加权综合模型改进拟合模型在逼近高程异常曲面时的拟合效果。由大地高求出准确的正常高，通常应用有限元法或多面函数法等单一拟合模型来拟合高程异常曲面。GPS水准的单一模型在逼近高程异常曲面时其拟合效果往往不稳定。通过对实测数据的计算分析可发现：有限元法与多面函数法的加权综合模型的逼近结果比单一模型的逼近结果有了一定的提高。由此可见综合模型对于局部地区解决GPS高程转化为正常高是有使用价值的。  相似文献   

Optimization routine for identification of model parameters in soil plasticity     

Hans Mattsson  Marek Klisinski  Kennet Axelsson 《国际地质力学数值与分析法杂志》2001,25(5):435-472

The paper presents an optimization routine especially developed for the identification of model parameters in soil plasticity on the basis of different soil tests. Main focus is put on the mathematical aspects and the experience from application of this optimization routine. Mathematically, for the optimization, an objective function and a search strategy are needed. Some alternative expressions for the objective function are formulated. They capture the overall soil behaviour and can be used in a simultaneous optimization against several laboratory tests. Two different search strategies, Rosenbrock's method and the Simplex method, both belonging to the category of direct search methods, are utilized in the routine. Direct search methods have generally proved to be reliable and their relative simplicity make them quite easy to program into workable codes. The Rosenbrock and simplex methods are modified to make the search strategies as efficient and user‐friendly as possible for the type of optimization problem addressed here. Since these search strategies are of a heuristic nature, which makes it difficult (or even impossible) to analyse their performance in a theoretical way, representative optimization examples against both simulated experimental results as well as performed triaxial tests are presented to show the efficiency of the optimization routine. From these examples, it has been concluded that the optimization routine is able to locate a minimum with a good accuracy, fast enough to be a very useful tool for identification of model parameters in soil plasticity. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

Modelling mechanical behaviour of limestone under reservoir conditions     

Lúcia Carvalho Coelho  Antonio Claudio Soares  Nelson Francisco F. Ebecken  Jos Luis Drummond Alves  Luiz Landau 《国际地质力学数值与分析法杂志》2006,30(14):1477-1500

High porosity and low permeability limestone has presented pore collapse. As fluid is withdrawn from these reservoirs, the effective stresses acting on the rock increase. If the strength of the rock is overcome, pore collapse may occur, leading to irreversible compaction of porous media with permeability and porosity reduction. It impacts on fluid withdrawal. Most of reservoirs have been discovered in weak formations, which are susceptible to this phenomenon. This work presents a study on the mechanical behaviour of a porous limestone from a reservoir located in Campos Basin, offshore Brazil. An experimental program was undergone in order to define its elastic plastic behaviour. The tests reproduced the loading path conditions expected in a reservoir under production. Parameters of the cap model were fitted to these tests and numerical simulations were run. The numerical simulations presented a good agreement with the experimental tests. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

Delayed plastic model for time‐dependent behaviour of materials     

Tsutomu Namikawa 《国际地质力学数值与分析法杂志》2001,25(6):605-627

A delayed plastic model, based on the theory of plasticity, is proposed to represent the time‐dependent behaviour of materials. It is assumed in this model that the stress can lie outside the yield surface and the conjugate stress called static stress is defined on the yield surface. The stress–strain relation is calculated based on the plastic theory embedding the static stress. Thus, the stress–strain relation of the model practically corresponds to that of the inviscid elastoplastic model under fairly low rate deformation. The delayed plastic model is coupled with the Cam‐clay model for normally consolidated clays. The performance of the model is then examined by comparing the model predictions with reported time‐dependent behaviour of clays under undrained triaxial conditions. It is shown that the model is capable of predicting the effect of strain rate during undrained shear and the undrained creep behaviour including creep rupture. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

Stress–dilatancy based modelling of granular materials and extensions to soils with crushable grains     

Antonio DeSimone  Claudio Tamagnini 《国际地质力学数值与分析法杂志》2005,29(1):73-101

Stress–dilatancy relations have played a crucial role in the understanding of the mechanical behaviour of soils and in the development of realistic constitutive models for their response. Recent investigations on the mechanical behaviour of materials with crushable grains have called into question the validity of classical relations such as those used in critical state soil mechanics. In this paper, a method to construct thermodynamically consistent (isotropic, three‐invariant) elasto‐plastic models based on a given stress–dilatancy relation is discussed. Extensions to cover the case of granular materials with crushable grains are also presented, based on the interpretation of some classical model parameters (e.g. the stress ratio at critical state) as internal variables that evolve according to suitable hardening laws. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

A non‐coaxial critical state soil model and its application to simple shear simulations     

Yunming Yang  H.S. Yu 《国际地质力学数值与分析法杂志》2006,30(13):1369-1390

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.  相似文献   

Finite deformation analysis of geomaterials     

Boris Jeremi&#x;  Kenneth Runesson  Stein Sture 《国际地质力学数值与分析法杂志》2001,25(8):809-840

The mathematical structure and numerical analysis of classical small deformation elasto–plasticity is generally well established. However, development of large deformation elastic–plastic numerical formulation for dilatant, pressure sensitive material models is still a research area. In this paper we present development of the finite element formulation and implementation for large deformation, elastic–plastic analysis of geomaterials. Our developments are based on the multiplicative decomposition of the deformation gradient into elastic and plastic parts. A consistent linearization of the right deformation tensor together with the Newton method at the constitutive and global levels leads toward an efficient and robust numerical algorithm. The presented numerical formulation is capable of accurately modelling dilatant, pressure sensitive isotropic and anisotropic geomaterials subjected to large deformations. In particular, the formulation is capable of simulating the behaviour of geomaterials in which eigentriads of stress and strain do not coincide during the loading process. The algorithm is tested in conjunction with the novel hyperelasto–plastic model termed the B material model, which is a single surface (single yield surface, affine single ultimate surface and affine single potential surface) model for dilatant, pressure sensitive, hardening and softening geomaterials. It is specifically developed to model large deformation hyperelasto–plastic problems in geomechanics. We present an application of this formulation to numerical analysis of low confinement tests on cohesionless granular soil specimens recently performed in a SPACEHAB module aboard the Space Shuttle during the STS‐89 mission. We compare numerical modelling with test results and show the significance of added confinement by the thin hyperelastic latex membrane undergoing large stretching. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献