共查询到19条相似文献，搜索用时 437 毫秒
1.
2.
基于非局部化方法的边坡稳定性分析 总被引：1，自引：1，他引：0
应变局部化问题已经成为岩土力学研究的一个重要课题,边坡的稳定性分析必然要考虑这种应变局部化对边坡的稳定性影响。结合Bazant等的非局部化模型方法,将其推广应用到弹黏塑性本构模型中去,考虑两种正则化机制的耦合作用,并采用有限元方法结合滑面应力法对边坡的稳定性进行分析。将数值分析结果与极限平衡法和基于线弹性本构的滑面应力法结果进行了对比,从分析结果来看文中的计算结果是合理的,对实际工程具有一定的指导意义。 相似文献
3.
不同应力路径下剪切带的数值模拟 总被引：4，自引：1，他引：3
采用回映应力更新算法,编写了基于伏斯列夫面的超固结黏土本构关系模型子程序,嵌入非线性有限元软件ABAQUS。通过对单元试验进行三轴压缩、三轴伸长及平面应变等问题的模型预测,再现了超固结黏土在不同初始超固结比和应力路径时的变形和强度特性,从而验证了子程序的正确性。借助该本构模型,对三轴压缩、三轴伸长及平面应变应力路径下超固结黏土体变形局部化问题,进行了三维数值模拟。分析结果表明：超固结黏土在三轴压缩及伸长状态时,土体变形局部化在应力-应变关系软化时出现,而平面应变状态时,在应力-应变关系硬化阶段出现,其超固结黏土的剪胀特性在剪切带的形成过程中起重要作用。 相似文献
4.
土体应变局部化现象的理论解析 总被引：7，自引：6，他引：1
引起土体失稳的应变局部化现象是在特定应力状态下,土体本构产生的分叉特性。基于有限变形理论推导了应变局部化产生的三维解析解。基于应变局部化的理论解析,分析了轴对称和平面应变条件下应变局部化现象在弹塑性硬化阶段的存在性以及剪切带的方向性。 理论分析表明,在轴对称条件下,土体应变局部化产生于土体应力-应变的软化阶段,而平面应变条件下,土体应变局部化一般出现在应力-应变的硬化阶段,其剪切带方向角的理论预测与Arthur等[1]建议值较为一致。 相似文献
5.
针对岩土工程材料应变软化问题及有限元对其数值计算时切线刚度矩阵负定造成求解困难的问题进行研究。建立了基于Drucker-Prager（D-P）强度准则的岩石弹塑性应变软化本构模型,本构积分算法采用一种完全隐式返回映射算法,它具有无条件稳定和精确的特点,详细论述了如何进行本构模型的程序化求解;考虑弧长法在判断切线刚度矩阵正定性导致效率低的缺点,在弹塑性增量有限元方程的迭代计算中尝试采用Newton-Raphson法和arc-length法（NR-AL法）联合迭代求解的思路,即在结构未达到极限荷载前采用NR迭代法,而当结构接近极限荷载时转换为AL法控制迭代,从而使结构越过峰值点进入软化区直至破坏,NR-AL法汲取了2者迭代求解中具有的优势;利用C++语言对所建应变软化模型的本构求解和弹塑性增量有限元方程迭代求解过程给予程序实现,应用所编程序进行数值计算,分析了D-P理想弹塑性模型、应变软化模型、应变硬化模型计算的应力-应变曲线的区别,同时将应变软化模型计算结果与试验数据进行了对比。研究结果表明：所建应变软化本构模型可以较好地模拟岩石材料的峰后软化特性,能够揭示峰后应变软化特性和破坏机制,同时NR-AL法能够求解由于应变软化造成的负刚度问题,也克服了单独使用弧长法时判断切线刚度矩阵正定性效率低的缺点。 相似文献
6.
在季冻区高等级公路路基场地取土,制成重塑土试件进行冻融循环试验。而后对经历不同冻融循环次数的土体进行三轴固结排水（CD）剪切试验,测得应力、应变数据。在大量试验数据的基础上,遵循岩土本构变化规律,以岩土弹塑性力学为理论基础,在p-q空间,以椭圆方程拟合体变屈服面,以抛物线方程拟合剪切屈服面,采用塑性关联流动,得到考虑冻融循环条件下融土的弹塑性本构方程。通过编程计算后,对比计算与试验结果,证明该双屈服面本构模型能够较为准确地预估冻融循环土体的应力-应变关系,该本构模型对季冻区路基土体长期稳定性分析和工程预测具有重要参考意义。 相似文献
7.
利用率相关与梯度塑性耦合本构模型,采用简谐波的分析方法对材料的应变局部化及材料的稳定性进行了研究,得到了二阶耦合模型在二维情况下的内尺度律的变化及其与材料稳定性的关系,得到了波长变化的下界及材料稳定性的条件;在此基础上,结合二阶耦合模型在一维情况下的内尺度律变化的规律,对其进行了对比研究,得到材料变形局部化带在一维与二维条件下的关系公式,这对材料变形局部化带宽的研究有积极意义。 相似文献
8.
岩石屈服之后,弹性参数随塑性变形而变化,即岩石具有弹塑性耦合特征。在弹塑性耦合框架和现有应变分类、定义、本构方程研究的基础上,主要做了以下工作：结合岩石压缩试验,从损伤和塑性变形角度分析了弹性参数随塑性演化的现象;重点研究了加载增量步的应变特征,将与加载应力增量对应的应变增量分解为符合广义Hooke定律的弹性应变增量 、不符合广义Hooke定律的弹性应变增量 和塑性应变增量 ;通过采用加载增量步之后的弹性参数将 与 构造弹性关系,分别在应变空间和应力空间中建立了适用于加载、中性变载和卸载条件的本构方程;对表征物质微观结构变化的内变量进行了讨论。所提出的应变分类和定义方法概念清晰,适合于应变强化阶段和应变软化阶段,且所建立的弹塑性耦合本构方程能够反映不同加载条件对弹性参数变化速率的影响,符合岩石的弹塑性耦合力学特性。 相似文献
9.
10.
基于摩尔-库仑模型的非线性本构模型的开发及其在应变局部化中的应用 总被引：1，自引：0，他引：1
以FLAC中的Mohr-Coulomb（M-C）本构模型为基础,在C++环境下实现了考虑拉伸截断的非线性本构模型的二次开发。针对应力及位移的分布规律,将理想弹塑性的Hoek-Brown（H-B）本构模型、M-C本构模型及自定义本构模型的数值解及H-B本构模型的理论解进行了比较,验证了自定义本构模型的正确性。采用弹-脆-塑性的自定义本构模型,以先加载,后挖洞的方式模拟了圆形巷道围岩的应变局部化过程,在计算模型的边界条件对称及均质本构参数的条件下,模拟出了4个对称的V型坑,模拟结果与静水压力条件下岩爆的现场观测及试验结果吻合,而采用接近脆性的H-B应变软化模型所得结果的对称性较差 相似文献
11.
A finite element model is developed for modelling coupled fluid expulsion/deformation behaviour of dewatering sediments subjected to external loadings under isothermal conditions. The non-linear deformation behaviour of the sediment (soil) skeleton is based on the force equilibrium equation in which the constitutive relationship of stress and strain is implemented by the modified Cam-Clay model in soil plasticity. The fluid flow behaviour in the model is described by the generalized porous media flow equation. The model allows temporal and spatial variations of porosity and permeability. The fluid viscosity and density are assumed to be temperature-dependent. The model also allows the development of single and multiple faults, depending upon the material (sediment and fluid) properties, loading and boundary conditions. Procedures are implemented for (1) updating the material properties such as porosity, permeability, fluid density and viscosity and (2) the development of faults which allow the formation of high-permeability conduits for fluid flow. The solution algorithm for displacements of the sediments and the excess pore (fluid) pressure is based on a residual load technique to handle the non-linear (elastic-plastic) deformation behaviour of the sediment skeleton. The model can be applied to one- and two-dimensional problems. Examples of a plane strain saturated sediment layer subjected to stepwise horizontal tractions versus time are given. 相似文献
12.
《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. 相似文献
13.
In this paper, the novel concept of probabilistic yielding is used for 1‐D cyclic simulation of the constitutive behavior of geomaterials. Fokker–Planck–Kolmogorov equation‐based probabilistic elastic–plastic constitutive framework is applied for obtaining the complete probabilistic (probability density function) material response. Both perfectly plastic and hardening‐type material models are considered. It is shown that when uncertainties in material parameters are taken into consideration, even the simple, elastic‐perfectly plastic model captures some of the important features of geomaterial behavior, for example, modulus reduction with cyclic strain, which, deterministically, is only possible with more advanced constitutive models. Furthermore, it is also shown that the use of isotropic and kinematic hardening rules does not significantly improve the probabilistic material response. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
14.
Kasper D. Fischer 《International Journal of Earth Sciences》2006,95(2):239-249
The Aegean–Anatolian region is characterised by an inhomogeneous deformation pattern with high strain rates and a high seismicity both at the boundaries and in the plate interior. This pattern is controlled by the geometry and rheology of the structural units involved and their tectonic setting. A numerical analysis with a finite-element model of the region is used to quantify the influence of different rheological parameters. Viscoelastic material behaviour is implemented for the mantle lithosphere, whereas the crust is modelled with an elastic–plastic rheology. The variation of the inelastic material properties (viscosity and plastic strength) quantifies the influence of these material parameters on the deformation, stress, and strain patterns. Comparison of the modelled results with geodetic and geophysical observations reveals that the viscosity of the mantle lithosphere is the key to explaining the inhomogeneous deformation pattern. The best-fit model yields a viscosity of 1020 Pa s beneath Anatolia, whereas adjacent regions have viscosities between 1021 and 1023 Pa s. The model also explains the intra-plate seismicity and the stress field as well as its partitioning into regions with strike-slip and normal faulting. The final model is in good agreement with seismological, geodetic, and geological observations. Local deviations can be tracked down to small-scale structures, which are not included in the model. 相似文献
15.
Implementation and applications for a constitutive numerical model on F‐75 silica sand, course silica sand and two sizes of glass beads compressed under plane strain conditions are presented in this work. The numerical model is used to predict the stress versus axial strain and volumetric strain versus axial strain relationships of those materials; moreover, comparisons between measured and predicted shear band thickness and inclination angles are discussed and the numerical results compare well with the experimental measurements. The numerical model is found to respond to the changes in confining pressure and the initial relative density of a given granular material. The mean particle size is used as an internal length scale. Increasing the confining pressure and the initial density is found to decrease the shear band thickness and increase the inclination angle. The micropolar or Cosserat theory is found to be effective in capturing strain localization in granular materials. The finite element formulations and the solution method for the boundary value problem in the updated Lagrangian frame (UP) are discussed in the companion paper. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
16.
The diffused and localized instabilities in sand under drained biaxial loading have been analyzed here following a plane strain bifurcation framework, where the rate independent material is defined using a generalized 3D non-associative constitutive model. This study is focused on how various instability modes emerge with respect to initial density, confining pressure, and the applied boundary conditions. Results from large deformation framework have been compared with those from small deformation approximation and the later is noticed to fail in capturing the emergence of diffused modes and predicts delayed onset of localization. The theoretical predictions compares well with existing experimental observations. 相似文献
17.
A hierarchical thermo‐hydro‐plastic constitutive model for unsaturated soils and its numerical implementation 
下载免费PDF全文
下载免费PDF全文 《国际地质力学数值与分析法杂志》2018,42(15):1785-1805
Unsaturated soils are highly heterogeneous 3‐phase porous media. Variations of temperature, the degree of saturation, and density have dramatic impacts on the hydro‐mechanical behavior of unsaturated soils. To model all these features, we present a thermo‐hydro‐plastic model in which the hydro‐mechanical hardening and thermal softening are incorporated in a hierarchical fashion for unsaturated soils. This novel constitutive model can capture heterogeneities in density, suction, the degree of saturation, and temperature. Specifically, this constitutive model has 2 ingredients: (1) it has a “mesoscale” mechanical state variable—porosity and 3 environmental state variables—suction, the degree of saturation, and temperature; (2) both temperature and mechanical effects on water retention properties are taken into account. The return mapping algorithm is applied to implement this model at Gauss point assuming an infinitesimal strain. At each time step, the return mapping is conducted only in principal elastic strain space, assuming no return mapping in suction and temperature. The numerical results obtained by this constitutive model are compared with the experimental results. It shows that the proposed model can simulate the thermo‐hydro‐mechanical behavior of unsaturated soils with satisfaction. We also conduct shear band analysis of an unsaturated soil specimen under plane strain condition to demonstrate the impact of temperature variation on shear banding triggered by initial material heterogeneities. 相似文献
18.
This paper studies the effects of initial fabric anisotropy of dry sand in simple shear deformation. The effects of anisotropy are taken into consideration through the modification of the mobilized friction in the Mohr–Coulomb‐type yield surface as a function of a fabric parameter. In addition, the constitutive model uses a gradient term that directly incorporates the effects of material length scale. The constitutive formulation is implemented into ABAQUS finite element code and used to simulate shearing of the dry sand under various conditions of simple shear. The numerical simulations show that while the shear stress response is affected by fabric anisotropy, its effects on strain localization in simple shear are minimal. This is in contrast to other devices such as the biaxial shear. The strain localization in simple shear is controlled more by the imposed boundary conditions. The use of material length scale is shown to remove the effects of strain localization in the shearing response. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
19.
Numerical simulations have been performed to investigate the strain-dependent behaviour of rheological and kinematical responses to flow of two-phase rocks using the commercial finite-difference program FLAC2D. It was assumed that the two phases have Maxwell rheology. Plane strain and velocity boundary condition, which produces a simple shear deformation, were also assumed. Two types of geometries were considered: strong phase supported (SPS) and weak phase supported (WPS). We calculated strain-dependent variations of effective viscosity and partitioning of strain rate, vorticity and kinematic vorticity number during deformation in both SPS and WPS structure models.The results show that the strain-dependent behaviour is largely influenced by the geometry of the composite. SPS models show both strain hardening and strain softening during the simulations, with strain hardening preceding strain softening. A critical shear strain is necessary to begin the strain softening behaviour. Strain hardening and strain softening are accompanied by a reduction and an increase of the partition of strain rate into the weak phase, respectively. On the other hand, WPS models show only weak strain hardening and strain softening, being the strain-dependent behaviour close to a steady state flow. In addition, the following results are obtained on vorticity and kinematic vorticity number; (1) in both SPS and WPS models the partition of vorticity into weak phase increases with progressive shear strain, i.e. the strong phase becomes less rotational, (2) in SPS models weak inclusions changes from sub-simple shear to super-simple shear with progressive strain, whereas the strong matrix changes from super-simple shear to sub-simple shear, (3) in WPS models the strong inclusions with high viscosity contrasts are less rotational but can be in super-simple shear condition to high strains.The observed strain-dependent behaviours have been compared with previous proposed analytical models. The degree of agreement is variable. Balshin and Ryshkewitch–Duckworth models are only applicable to SPS models. Ji-generalized mixture rule model is applicable to both models.The results suggest that polyphase rocks with SPS structure during ductile shear deformation respond as strain softening materials, after an initial strain hardening stage that may drive to the strain localization into the material. 相似文献