共查询到20条相似文献，搜索用时 15 毫秒
1.
A Physical and Numerical Investigation on the Stability of Shallow Tunnels in Strain Softening Media 总被引：11，自引：0，他引：11
Summary The behaviour up to failure of shallow underground openings is discussed on the basis of some laboratory, small-scale model tests and of finite element simulation. The experimental results are first illustrated. They were obtained from two-dimensional (plane strain) and three-dimensional tunnel models tested under standard gravity conditions. Then, the phenomenon of strain localisation that characterizes the medium surrounding the model tunnels is discussed, recalling two alternative approaches for its numerical interpretation. On this basis, a finite element procedure for strain softening analyses is outlined and applied to the simulation of the tests in both two- and three-dimensional conditions. The comparison between experimental and numerical results leads to some conclusions on the influence of strain localisation on the overall behaviour of shallow tunnels and on the stability of their headings. 相似文献
2.
E. Alonso L. R. Alejano F. Varas G. Fdez‐Manin C. Carranza‐Torres 《国际地质力学数值与分析法杂志》2003,27(13):1153-1185
A literature review has shown that there exist adequate techniques to obtain ground reaction curves for tunnels excavated in elastic‐brittle and perfectly plastic materials. However, for strain‐softening materials it seems that the problem has not been sufficiently analysed. In this paper, a one‐dimensional numerical solution to obtain the ground reaction curve (GRC) for circular tunnels excavated in strain‐softening materials is presented. The problem is formulated in a very general form and leads to a system of ordinary differential equations. By adequately defining a fictitious ‘time’ variable and re‐scaling some variables the problem is converted into an initial value one, which can be solved numerically by a Runge–Kutta–Fehlberg method, which is implemented in MATLAB environment. The method has been developed for various common particular behaviour models including Tresca, Mohr–Coulomb and Hoek–Brown failure criteria, in all cases with non‐associative flow rules and two‐segment piecewise linear functions related to a principal strain‐dependent plastic parameter to model the transition between peak and residual failure criteria. Some particular examples for the different failure criteria have been run, which agree well with closed‐form solutions—if existing—or with FDM‐based code results. Parametric studies and specific charts are created to highlight the influence of different parameters. The proposed methodology intends to be a wider and general numerical basis where standard and newly featured behaviour modes focusing on obtaining GRC for tunnels excavated in strain‐softening materials can be implemented. This way of solving such problems has proved to be more efficient and less time consuming than using FEM‐ or FDM‐based numerical 2D codes. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
3.
针对岩土工程材料应变软化问题及有限元对其数值计算时切线刚度矩阵负定造成求解困难的问题进行研究。建立了基于Drucker-Prager（D-P）强度准则的岩石弹塑性应变软化本构模型,本构积分算法采用一种完全隐式返回映射算法,它具有无条件稳定和精确的特点,详细论述了如何进行本构模型的程序化求解;考虑弧长法在判断切线刚度矩阵正定性导致效率低的缺点,在弹塑性增量有限元方程的迭代计算中尝试采用Newton-Raphson法和arc-length法（NR-AL法）联合迭代求解的思路,即在结构未达到极限荷载前采用NR迭代法,而当结构接近极限荷载时转换为AL法控制迭代,从而使结构越过峰值点进入软化区直至破坏,NR-AL法汲取了2者迭代求解中具有的优势;利用C++语言对所建应变软化模型的本构求解和弹塑性增量有限元方程迭代求解过程给予程序实现,应用所编程序进行数值计算,分析了D-P理想弹塑性模型、应变软化模型、应变硬化模型计算的应力-应变曲线的区别,同时将应变软化模型计算结果与试验数据进行了对比。研究结果表明：所建应变软化本构模型可以较好地模拟岩石材料的峰后软化特性,能够揭示峰后应变软化特性和破坏机制,同时NR-AL法能够求解由于应变软化造成的负刚度问题,也克服了单独使用弧长法时判断切线刚度矩阵正定性效率低的缺点。 相似文献
4.
The aim of the present paper is to numerically analyse the behaviour of frozen sand by using a viscoplastic constitutive model with strain softening. A constitutive model has been developed introducing the stress history tensor which is a functional of the stress history, with respect to a generalized time measure. It is shown that Adachi and Oka's model is applicable to the results of triaxial tests on a frozen Toyoura sand at different strain rates. First, the instability of the model is discussed within the framework of bifurcation theory. The model is then implemented into a FEM code to numerically simulate the behaviour under plane strain conditions. From the numerical results, it is revealed that the formation of shear bands is possible and the characteristics of strain localization, such as shear banding, depend on the strain rates. 相似文献
5.
地基渐进破坏及极限承载力的Cosserat连续体有限元分析 总被引：2，自引：2，他引：2
利用Cosserat连续体理论和所发展的有限元数值方法,模拟了地基由应变软化引起的以应变局部化为特征的渐进破坏过程,并从等价塑性应变的发展变化,阐述了渐进破坏过程对所能发挥的极限承载能力的影响。结果表明,Cosserat连续体模型能有效地模拟由应变软化引起以应变局部化为特征的渐进破坏现象,对地基等土工结构物有必要进行渐进破坏分析。同时指出,在求解软化型土体地基的极限承载力时,如果仍按传统的极限平衡或极限分析理论进行分析,可能得出偏于危险的结果。 相似文献
6.
不同应力路径下剪切带的数值模拟 总被引：4，自引：1，他引：3
采用回映应力更新算法,编写了基于伏斯列夫面的超固结黏土本构关系模型子程序,嵌入非线性有限元软件ABAQUS。通过对单元试验进行三轴压缩、三轴伸长及平面应变等问题的模型预测,再现了超固结黏土在不同初始超固结比和应力路径时的变形和强度特性,从而验证了子程序的正确性。借助该本构模型,对三轴压缩、三轴伸长及平面应变应力路径下超固结黏土体变形局部化问题,进行了三维数值模拟。分析结果表明：超固结黏土在三轴压缩及伸长状态时,土体变形局部化在应力-应变关系软化时出现,而平面应变状态时,在应力-应变关系硬化阶段出现,其超固结黏土的剪胀特性在剪切带的形成过程中起重要作用。 相似文献
7.
砂土的渐进破坏及其数值模拟 总被引：1，自引：1，他引：0
紧密砂土在排水剪切过程中由于材料的应变软化特性,会出现变形局部化并伴随剪切带的产生。对此类问题进行数值模拟时会遇到很多困难,比如如何模拟紧砂的应变软化特性,如何处理有限元数值分析过程中刚度矩阵的非正定导致计算结果发散问题,以及如何处理应力不连续的问题,拟针对以上问题进行了探讨。 相似文献
8.
The paper mainly concerns the mechanical response of 2D dry dense sand specimens under shock loading. The problem is numerically analysed by means of a SEM dynamic code, within which an already conceived non‐local viscoplastic constitutive model characterized by a non‐associated flow rule and by an anisotropic strain hardening has been implemented. In particular the strain localization and time dependency of the material mechanical response are taken into consideration. Both rapid/static loading and dynamic histories are numerically simulated. In the first case, the time dependency of the material mechanical response can be captured by neglecting inertial effects, while in the second one the two factors are superimposed and the propagation of the stress waves within the specimen is discussed. The interest of these analyses derives from the fact that the diffusion phenomenon takes place within a specimen already localized. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
9.
Material behaviour that exhibits characteristics of creep induced by a spontaneous mineral dissolution enhanced by material damage is studied. It is believed that the characteristic rates of the chemical processes involved determine the time‐rate dependence of the resulting strain. A basic model of a combined chemo‐plastic softening and chemically enhanced deviatoric strain hardening for saturated geomaterials is presented. Chemical softening is postulated to occur as a consequence of the net mass removal resulting from dissolution and precipitation of specific minerals occurring at the damage‐generated inter‐phase interfaces. Closed and open systems are discussed. In the former case, deformation at constant stress results entirely from a local compensation mechanism between the chemical softening and strain hardening. The classical three stages of creep are interpreted in terms of mechanisms of dissolution and precipitation, as well as the variation in the reaction surface areas involved in the mass exchange. In an open system, the above local mechanism is enhanced by the removal of mass via diffusion of species affecting the mass balance. Such a system is addressed via a boundary value problem as shown in an example. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
10.
This paper studies the excavation of a spherical cavity subjected to hydrostatic initial stresses in the infinite homogeneous and isotropic rock mass with strain‐softening Mohr–Coulomb (M‐C) and Hoek–Brown (H‐B) behaviors. Numerical solutions of the spherical cavity are obtained and the application to determining stress–strain curve of strain‐softening M‐C and H‐B rock mass is studied. A closed‐form solution for the elastic–brittle–plastic medium is introduced first, and then a numerical procedure that simplifies the strain‐softening process into a series of brittle–plastic ones is presented. The approach is validated against the facts that the strain‐softening process evolves into a brittle–plastic one when the softening slope is very steep, whereas it evolves into an elasto‐plastic one when the softening slope approaches zero. Numerical solutions for the prediction of displacements and stresses around the spherical cavity in the strain‐softening M‐C and H‐B rock mass are presented. On the basis of the analysis of the spherical cavity in strain‐softening rock mass, the stress–strain relationship at an infinitesimal cube around the cavity is obtained and discussed with different evolution laws for the strength parameters considered. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
11.
《Geomechanics and Geoengineering》2013,8(4):235-251
The bearing capacity of footing has been studied by both conventional and numerical methods by many researchers. However, degradation of the microstructure of material, that is, a change in the microstructure of the soil, has not been adequately taken into account. Degradation of microstructure causes strain softening of materials and it leads to strain localization such as shear bands and slip bands. From an engineering point of view the strain localization is crucial because it is a precursor of failure. In the present study, finite element analyses of the bearing capacity of a shallow foundation on homogeneous and inhomogeneous saturated clay strata have been conducted using an elasto-viscoplastic soil constitutive model of microstructure change. A series of analyses of footing on clay deposit with different microstructure parameters have been carried out. Numerical results show that strain localization can be predicted during the loading of rigid footing on highly structured soil and strain localization affects the footing–soil interaction. The effects of footing roughness on the failure mechanism are also discussed in the study. 相似文献
12.
Constitutive models for rocks and soils typically incorporate some form of strain softening. Moreover, many plasticity models for frictional materials use a nonassociated flow rule. Strain softening and nonassociated flow rules can cause loss of well-posedness of the initial-value problem, which can lead to a severe mesh dependence in simulations and poor convergence of the iterative solution procedure. The inclusion of viscosity, which is a common property of materials, seems a natural way to restore well-posedness, but the mathematical properties of a rate-dependent model, and therefore the effectiveness with respect to the removal of mesh dependence, can depend strongly on how the viscous element is incorporated. Herein, we show that rate-dependent models, which are commonly applied to problems in the Earth's lithosphere, such as plate tectonics, are very different from the approach typically adopted for more shallow geotechnical engineering problems. We analyse the properties of these models under dynamic loadings, using dispersion analyses and one-dimensional finite difference analyses, and complement them with two-dimensional simulations of a typical strain localisation problem under quasi-static loading conditions. Finally, we point out that a combined model, which features two viscous elements, may be the best way forward for modelling time-dependent failure processes in the deeper layers of the Earth, since it not only enables modelling of the creep characteristics typical of long-term behaviour but also regularises the initial/boundary-value problem. 相似文献
13.
非共轴本构模型在地基承载力数值计算中若干影响因素的探讨 总被引：1，自引：0，他引：1
土体在剪切变形过程中产生主应力方向的旋转时,主应变增量方向与主应力方向之间存在着非共轴现象,然而传统的弹塑性本构模型未能考虑该现象的影响。通过在屈服面的切线方向增加一项非共轴塑性应变增量,即可实现对非共轴现象的反映。采用显式积分算法和自动分步方法,将非共轴本构模型运用到桶形基础地基承载力问题的有限元计算中,并讨论了流动法则、内摩擦角、膨胀角等因素与非共轴模型的联系。计算结果表明：采用有限元程序默认容许误差时,该本构模型可达到理想的收敛精度,并且,该模型对关联、非关联流动法则均适用。采用共轴模型进行数值计算时,不同流动法则对计算结果的影响可以忽略;采用非共轴模型时,不同流动法则的计算结果之间存在差异。非共轴现象对地基承载力-位移曲线具有软化作用,并且,该软化作用在采用非关联流动法则时变得更加明显 相似文献
14.
《Journal of Structural Geology》2002,24(6-7):1045-1053
The development of shear zones initiating on random weaker initial perturbations is modelled numerically for low Deborah number viscoelastic materials, considering the influence of effective viscosity contrast, power law rheology, strain softening, and different imposed bulk deformation geometries, ranging from pure to simple shear. Conjugate shear zones initiate at ∼90° to one another, and rotate with increasing bulk deformation, the basic pattern not being markedly influenced by the vorticity of imposed deformation. The rate of propagation of individual conjugate shear zones is little affected by increased effective viscosity contrast between matrix and inclusion but is promoted by power-law rheology. However, the most marked effect is observed for strain softening behaviour, where rapid propagation produces straighter and narrower shear zones. The localisation of strain is reflected in a correspondingly heterogeneous stress distribution. In particular, mean stress or pressure is higher in the extending, near planar, weaker zones of localised shear. Melting of gneissic or pelitic compositions is pressure dependent. With free water present, increased pressure promotes melting, whereas the opposite is true for water-absent melting. For water-present conditions, a positive feedback could develop between localised shearing, increased pressure and partial melting. This is potentially more effective in concentrating melt in shear zones than shear heating, where melt-related softening has a negative feedback effect. 相似文献
15.
We investigate the internal deformation of orogenic wedges growing by frontal accretion with a two-dimensional numerical model. Our models are limited to crustal deformation and assume a horizontal detachment as observed for various natural orogens (e.g. Alaska and Costa Rica). The model wedges develop as a result of convergence of a brittle sediment layer in front of a strong backstop. We find that our reference model develops in-sequence forward-thrusts which propagate upward from the basal detachment. For this reference model we investigate the sensitivity of shear zone activity to surface processes and strain softening. Model results show that diffusive or slope dependent erosion enhances material transport across the wedge and slows down forward propagation of the deformation front. Frictional strain softening focuses deformation into narrow shear zones and enhances displacement along them. This has also been postulated for natural thrusts such as the Glarus thrust in the Swiss Alps and the Moine thrust in the Scottish Caledonides. A second series of models investigates the effects of regularly spaced weak inclusions within the sediment layer which simulate remnants of previous deformation phases. These inclusions facilitate and focus internal deformation, influence the thrust dip and thrust vergence and enable thrust reactivation in the internal part of the wedge. Our results show that inactive thrusts in the internal part of the wedges may be reactivated in models with diffusive surface processes, strain softening or weak inclusions. Thrust reactivation occurs as models seek to maintain their critical taper angle. First order characteristics of our numerical models agree well with natural orogenic wedges and results from other numerical and analogue models. 相似文献
16.
A distinct element method numerical investigation of compaction processes in highly porous cemented granular materials 总被引：1，自引：0，他引：1
In this paper, the results of an oedometric numerical test campaign, performed by means of a 3D Discrete Element Code on idealised cemented granular cylindrical specimens, are illustrated. The idealised microstructure taken into account is characterised by the following: (i) rigid grains bonded to one another; (ii) a high void ratio; and (iii) two different families of voids: the micro and the macro‐voids. The compaction process developing within the specimens, as well as the localization along tabular zones of pure compressive deformation (compaction banding) that in some cases takes place, are discussed. The influence on the evolution of this peculiar strain localization process of many microstructural/numerical parameters like material porosity, macro‐void size, the constitutive relationship adopted for the bonds and the bond damage rate is analysed. Tests for different values of porosity were run. Below a certain porosity threshold value, the onset of mixed modes of localisation was detected whereas the increase in the macro‐void size is observed to favour the onset of instability. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
17.
膨胀土边坡受降雨影响产生膨胀变形,是典型的非饱和土多场耦合问题。为探究降雨入渗对其渐进性破坏的失稳过程,基于饱和-非饱和渗流理论、膨胀土弹塑性本构关系和应变软化理论,利用应变软化模型、FLAC3D二次开发平台和内置FISH语言,提出了一种综合考虑非饱和渗流、膨胀变形和应变软化的多场耦合数值分析法。结合工程实例,通过该方法探讨了降雨入渗条件下膨胀土边坡非饱和渗流、位移响应及渐进性破坏的变化规律。结果表明:膨胀变形和应变软化受控于非饱和渗流的时空分布,对边坡位移响应过程影响显著,也易导致饱和-非饱和分界带形成剪应力集中区。膨胀土边坡渐进性破坏由局部破坏转变为整体性失稳,其塑性破坏区首先随悬挂型暂态饱和区的变化向坡内扩展,雨后逐渐形成第二条由坡脚向坡顶扩展的滑动带,呈现出多重滑动性和后退牵引式的破坏特征。 相似文献
18.
19.
Guang Y. Sheu 《国际地质力学数值与分析法杂志》2006,30(8):823-842
This study presents a preliminary development of a direct back analysis procedure by the meshless local Petrov–Galerkin (MLPG) method and Bayesian statistics and the application of resulting procedure to characterize soil properties using laboratory results. As compared to direct back analyses that are based on the finite element method (FEM), it is intended to show that MLPG‐based direct back analyses may be more suitable for some kinds of characterization problems; for example, involving a complex subsurface stratification or the characterization of soil properties of just an inclusion of a soil profile. The existing MLPG method is first slightly modified to analyse time‐dependent problems. Using the resulting method, quantities to be characterized are evaluated so that they give numerical results as close to measured data as possible. The Akaike information criterion is introduced for simplifying the evaluation. A one‐dimensional finite strain consolidation problem is introduced to do an error analysis for prediction by the proposed MLPG method. Another example illustrates experiences of performing an MLPG‐based direct back analysis. Comparison of MLPG‐based and FEM‐based direct back analyses is taken. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
20.
A plain strain problem of an isotropic elastic liquid-saturated porous medium in poroelasticity has been studied. The eigenvalue
approach using the Laplace and Fourier transforms has been employed and these transforms have been inverted by using a numerical
technique. An application of infinite space with concentrated force at the origin has been presented to illustrate the utility
of the approach. The displacement and stress components in the physical domain are obtained numerically. The results are shown
graphically and can be used for a broad class of problems related to liquid-saturated porous media. 相似文献