共查询到20条相似文献，搜索用时 31 毫秒
1.
In this paper a new finite element formulation for numerical analysis of diffused and localized failure behavior of saturated and partially saturated gradient poroplastic materials is proposed. The new finite element includes interpolation functions of first order (C1) for the internal variables field while classical C0 interpolation functions for the kinematic fields and pore pressure. This finite element formulation is compatible with a thermodynamically consistent gradient poroplastic theory previously proposed by the authors. In this material theory the internal variables are the only ones of non-local character. To verify the numerical efficiency of the proposed finite element formulation, the non-local gradient poroplastic constitutive theory is combined with the modified Cam Clay model for partially saturated continua. Thereby, the volumetric strain of the solid skeleton and the plastic porosity are the internal variables of the constitutive theory. The numerical results in this paper demonstrate the capabilities of the proposed finite element formulation to capture diffuse and localized failure modes of boundary value problems of porous media, depending on the acting confining pressure and on the material saturation degree. 相似文献
2.
In the coupled finite element analysis of the two-phase, saturated soils the pore pressure could decrease below the atmospheric pressure in particular when the plastic dilation of the skeleton is accounted for. The consequent reduction of the bulk modulus of the liquid phase is here modelled in a simplified manner through a tension cut off on the pore pressure. Two alternative approaches are outlined in the framework of the finite element method. The first one introduces a unilateral constraint on the nodal pore pressures and involves manipulations on matrices of the assembled mesh, whilst the second approach requires calculations solely at the element level. The effects of the pore pressure tension cut off are illustrated through some undrained test examples in elastic–plastic regime. 相似文献
3.
An elastoplastic model has been developed for the finite elements modelling of repeated load triaxial tests. This model is based on the shakedown theory established by Zarka for metallic structures. To the previous works, which were based on the Drucker–Prager yield surface and the plastic potential of Von Mises, a compression cap has been added to each one. The model straightforwardly determines the purely elastic state or the elastic shakedown state or the plastic shakedown state and calculates the deviatoric and the volumetric plastic strains. The calibration of the elastoplastic model has been carried out with DEM simulations and an unbound granular material for roads under repeated load triaxial tests using finite element method. The calculations underline the capabilities of the model to take into account, with a unique formalism, the accumulation of the deviatoric and volumetric plastic strains along the loading cycles. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
4.
A coupled elastic–plastic finite element analysis based on simplified consolidation theory for unsaturated soils is used to investigate the coupling processes of water infiltration and deformation. By introducing a reduced suction and an elastic–plastic constitutive equation for the soil skeleton, the simplified consolidation theory for unsaturated soils is incorporated into an in-house finite element code. Using the proposed numerical method, the generation of pore water pressure and development of deformation can be simulated under evaporation or rainfall infiltration conditions. Through a parametric study and comparison with the test results, the proposed method is found to describe well the characteristics during water evaporation/infiltration into unsaturated soils. Finally, an unsaturated soil slope with water infiltration is analyzed in detail to investigate the development of the displacement and generation of pore water pressure. 相似文献
5.
围压循环变化将引起土体孔压累积,加剧土体累积塑性应变,导致地基灾变事故。为研究围压循环变化对饱和粉质黏土变形特性的影响,利用GDS双向振动三轴仪,进行部分排水时饱和粉质黏土的变围压循环三轴试验,研究不同应力路径斜率和循环动应力比时饱和粉质黏土的孔压、轴向累积塑性应变和体应变发展规律,建立部分排水条件时考虑循环围压和循环动应力耦合作用的地基粉质黏土累积应变数学模型。试验结果表明：部分排水时饱和粉质黏土孔压、轴向累积塑性应变和体应变均随循环动应力比和应力路径斜率的增大而增大;动孔压比随着振动次数的增加明显分为急剧增加、快速下降和持续平稳3个阶段;轴向累积塑性应变和体应变均随振动次数的增加而增大。当加载次数超过2 500次时,粉质黏土孔压趋于平稳,变形速率略有降低,但变形持续增加。补充试验结果表明,所建立模型与试验结果具有较高的一致性。研究结果将为交通循环荷载导致的地基灾变控制技术提供理论依据。 相似文献
6.
A three-phase hydro-mechanical model for hydraulic fracturing is proposed. Three phases include: porous solid, fracturing fluid and host fluid. Discontinuity is handled using extended finite element method (XFEM) while cohesive crack model is used as fracturing criterion. Flow through fracture is defined as one-dimensional laminar flow, and flow through porous medium (host reservoir) is defined as two-dimensional Darcy flow. Coupling between two fluids in each space, fracture and pore, is captured through capillary pressure–saturation relationship, while the identical fluids in fracture and pore are coupled through a so-called leak-off mass transfer term. Coupling between fluids and deformation is captured through compatibility of volumetric strain of fluids within fracture and pore, and volumetric strain of the matrix. Spatial and temporal discretisation is achieved using the standard Galerkin method and the finite difference technique, respectively. The model is verified against analytical solutions available from literature. The leaking of fracturing fluid into the medium and suction of porous fluid into the fracture around the tip, are investigated. Sensitivity analyses are carried out for cases with slow and fast injection rates. It is shown that the results by single-phase flow may underestimate the leak-off. 相似文献
7.
对于低路堤公路,合理地预测车辆荷载引起的地基塑性累积变形对于低路堤公路的设计有着重要的意义。从机制上讲,可将地基的塑性应变分为塑性体积应变和塑性剪切应变。其中对于塑性体积应变,提出采用不排水条件下的残余孔压模型与一维固结理论相耦合进行计算;对于塑性剪切应变,考虑了塑性体积应变引起的土体硬化,提出了修正的塑性剪切应变计算模型。然后基于所提模型,对萧山黏土的不排水与排水条件下的动三轴试验进行模拟,初步验证了所提模型的有效性。对于实际工程问题,首先运用数值的方法求得车辆荷载在地基中引起的附加应力,然后运用所提模型计算各土层的塑性应变,再结合分层总和的方法求得地基的累积塑性变形。以日本佐贺机场1号测点为例,说明了模型在工程边值问题中的有效性。 相似文献
8.
Progressive failure in long natural clay slopes with an underlying weak shear zone is studied. The clay slope is assumed to be elastic–perfectly plastic, whereas the shear zone may have an elastic–plastic strain softening behaviour. The failure will occur due to the development of large shear deformations in the weak zone. The general behaviour, including the possibility of progressive failure, can be studied by the use of a one-dimensional finite element model. Other interesting matters which can be studied are critical load, critical disturbance and critical length of the slope. Numerical examples show the overall behaviour due to different residual shear strengths of the weak layer, end disturbance or change of pore pressure in the weak layer. Although simple, the proposed finite element model provides a practically applicable tool for the prediction of whether progressive failure will occur and in which manner the local failure will propagate. 相似文献
9.
10.
This paper presents the theoretical background of an elastic electro-osmosis consolidation model for saturated soils experiencing large strains, which considers volumetric strains induced by changes in both the hydraulic and electric driven pore water flows. Three fully coupled governing equations, considering the soil mechanical behaviour, pore water transport and electrical field, and their numerical implementation within an updated Lagrangian finite element formulation, are presented. The proposed model is first verified against a classical one-dimensional analytical solution for electro-osmosis consolidation to demonstrate its accuracy and efficiency. Then, various numerical examples are investigated to study the deformation characteristics and time dependent evolution of excess pore pressure. Finally, the importance of considering large strains in a consistent and proper way is demonstrated, and differences compared to models based on small strain theory are highlighted. 相似文献
11.
This paper presents a numerical technique for the analysis of the cone penetration test by means of the commercial finite element code ABAQUS. The von Mises yield criterion with its associated flow rule is assumed to model the plastic behaviour of elastoplastic undrained clays. An explicit finite element scheme is used to efficiently perform a large number of loading increments and to simplify the treatment of contact. An Arbitrary Langrangian–Eulerian (ALE) scheme is adopted to preserve the quality of mesh throughout the numerical simulation. A volumetric weighting algorithm adjusts the relative positions of nodes after each loading increment. This prevents mesh over distortion and allows the simulation to run continuously. The variation of the cone resistance is examined in relation to various parameters such as the in situ stress state, shaft and cone face roughness, and the material strength when steady state conditions have been reached. The trends of these variations are highlighted and compared with those found by other researchers. This technique can be extended to analyse the plastic behaviour of elastoplastic sands often modelled using either the Drucker–Prager yield criterion or a critical state model. 相似文献
12.
A multiphase coupled elasto‐viscoplastic finite element analysis formulation, based on the theory of porous media, is used to describe the rainfall infiltration process into a one‐dimensional soil column. Using this framework, we have numerically analyzed the generation of pore water pressure and deformations when rainfall is applied to the soil. A parametric study, including rainfall intensity, soil–water characteristic curves, and permeability, is carried out to observe their influence on the changes in pore water pressure and volumetric strain. From the numerical results, it is shown that the generation of pore water pressure and volumetric strain is mainly controlled by material parameters α and n′ that describe the soil–water characteristic curve. A comparison with the laboratory results shows that the proposed method can describe very well the characteristics observed during the experiments of one‐dimensional water infiltration into a layered unsaturated soil column. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
13.
Akira Murakami Takayuki Shuku Shin‐ichi Nishimura Kazunori Fujisawa Kazuyuki Nakamura 《国际地质力学数值与分析法杂志》2013,37(11):1642-1669
A computational method, incorporating the finite element model (FEM) into data assimilation using the particle filter, is presented for identifying elasto‐plastic material properties based on sequential measurements under the known changing traction boundary conditions to overcome some difficulties in identifying the parameters for elasto‐plastic problems from which the existing inverse analysis strategies have suffered. A soil–water coupled problem, which uses the elasto‐plastic constitutive model, is dealt with as the geotechnical application. Measured data on the settlement and the pore pressure are obtained from a synthetic FEM computation as the forward problem under the known parameters to be identified for both the element tests and the ground behavior during the embankment construction sequence. Parameter identification for elasto‐plastic problems, such as soil behavior, should be made by considering the measurements of deformation and/or pore pressure step by step from the initial stage of construction and throughout the deformation history under the changing traction boundary conditions because of the embankment or the excavation because the ground behavior is highly dependent on the loading history. Thus, it appears that sequential data assimilation techniques, such as the particle filter, are the preferable tools that can provide estimates of the state variables, that is, deformation, pore pressure, and unknown parameters, for the constitutive model in geotechnical practice. The present paper discusses the priority of the particle filter in its application to initial/boundary value problems for elasto‐plastic materials and demonstrates a couple of numerical examples. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
14.
基于卸荷岩体力学基本理论,运用三维有限差分软件,提出以塑性体积应变增量变化区间与不同卸荷区域岩体力学参数相对应的卸荷模拟方法,模拟了不同侧压系数条件下开挖洞周附近围岩最大体积应变增量以及围岩损伤范围的变化情况。数值模拟结果表明：当侧压系数较低时,洞周附近围岩的塑性体积应变增量区域的外边界近似为竖向椭圆形,随着侧压系数的增大,其外边界逐渐向圆形变化,最后变为横向椭圆形。侧压系数越大,围岩的卸载效应越显著,洞周附近围岩的塑性体积应变增量就越大,其所产生的围岩损伤范围也相应越大。利用所获得的结果,可以为深埋隧洞的稳定性分析以及支护设计提供依据。 相似文献
15.
正常固结黏土中平板锚基础的吸力和抗拉力 总被引：3，自引：1，他引：2
平板锚是新近出现的一种系泊深海浮式结构的基础型式。当黏土地基中的平板锚承受上拔力时,平板上、下表面超静孔压差形成的吸力使其抗拉承载力显著增加,对于风浪等快速加载条件尤其如此。利用有限元软件ADINA建立有效应力形式的轴对称动力有限元模型,研究圆形平板锚在缓慢加载与快速加载时的超静孔压分布与地基破坏型式。快速加载时重黏土和高岭土两组典型正常固结土样所得极限承载力系数与塑性极限分析解一致。进而通过变动参数分析,讨论加载速率和埋深对吸力和总抗拉力极限值的影响,并给出排水和不排水加载条件对应的临界加载速率。结果表明,不排水加载条件的总抗拉力可能达到排水总抗拉力的3倍。 相似文献
16.
17.
土工织物与塑料排水板联合处理软基的效果分析 总被引：7，自引：0，他引：7
基于平面应变条件下的Biot固结理论,采用非线性有限元方法对土工织物加筋与塑料板排水联合使用处理软土路基的效果和机制进行了探讨。对比计算表明：在相同固结时间内加设排水板后路堤的竖向沉降大于无排水板情况下的相应值,路堤坡趾处的隆起变小,路堤的侧向位移也变小。塑料排水板极大地缩短了孔压的消散时间,加快了软土的固结进程,使软土地基的承载力得到较快的提高,从而提高了路堤的稳定性。加设排水板后筋材所分担的抗拉力也比无排水板情况下的抗拉力低。为了更好地发挥土工织物的加固效果,土工织物应尽量布置在高应力区域。加设排水板可以缩短路堤各固结阶段的静置时间而达到优化设计。 相似文献
18.
A model of chemo-osmotic consolidation of clays in multi-dimensional domains is presented, with volumetric strains induced by both changes in the chemistry and osmotically driven pore water flow considered. Three fully coupled governing equations considering force equilibrium, pore water transport and chemical transport are presented and solved using the transient finite element method. The proposed approach is verified via consideration of chemo-osmotic consolidation of a compacted clay landfill liner and then applied to investigation of a hypothetical case with a local leachate leak in the compacted clay liner. An assessment of the impact of the two-dimensional nature of the system is made. The consolidation process is found to be dominated by osmotic consolidation in the early stages and subsequently by chemical consolidation. It is found that the surface settlement and the leachate concentration in the compacted clay liner may be highly overestimated by a one-dimensional analysis. Moreover, the peak negative excess pore pressure predicted by the two-dimensional solution remains in a shallow region under the leak while in the one-dimensional solution it moves progressively downwards. 相似文献
19.
为研究近断层脉冲地震动中竖向加速度对砂土场地液化的影响,基于有限元平台OpenSees开发的边界面塑性本构模型,建立了动单剪单元试验模型和饱和砂土三维有限元模型。选取台湾Chi-Chi地震中10条具有速度脉冲特性的地震波,对比分析了水平双向脉冲波与三向脉冲波作用下土柱竖向位移、循环应力比、孔压比及等效循环周数的差异性,继而明确了脉冲地震动中竖向加速度对砂土液化的影响规律。研究表明,三向脉冲地震波中竖向加速度分量对场地永久位移值影响较小,但使永久位移的发展持时明显增大;土柱循环应力比受竖向地震动影响较小,因此分析脉冲地震动对场地剪切特性的影响时,可将三向脉冲地震动简化为水平双向地震动;考虑竖向地震动的三向脉冲地震波引起的孔压比变化幅度较大,孔压消散时间较长;三向脉冲地震波对应的等效循环周数较大,地震动发展持时长,可认为竖向加速度对场地液化有促进作用。 相似文献
20.
不同应力路径下超固结黏土试样变形局部化分析 总被引：2，自引：1，他引：1
基于改进伏斯列夫面超固结黏土三维本构模型,利用有限元软件ABAQUS材料子程序接口,采用回映应力更新算法,实现了该模型在有限元分析中的应用。通过该模型与比奥固结理论的耦合,对超固结比为8的超固结黏土在三轴压缩、三轴伸长及平面应变应力条件下的变形局部化问题,进行了水-土耦合弹塑性有限元分析。分析结果表明：剪切带带内、带外点经历不同应力路径;剪切带带外单元经历了体缩、剪胀及被吸水体缩过程,而剪切带带内单元一直保持剪胀趋势;剪切带的形成伴随着剪胀,剪切带内、外出现了负的孔压,且孔压的分布也具有局部化特性。关于剪切带带内、外的孔隙水压及体变变化趋势与剪切速率有关,而平面应变介于三轴压缩与三轴伸长之间,但平面应变较早出现剪切带。孔隙水的迁移速度影响剪切带带内单元的剪胀,进而影响剪切带的形成及发展;而围压和弱单元位置也对剪切带的形成也有影响。 相似文献