首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
For seepage failures of dike due to water level-up and rainfall, surface infiltration and strength change induced by suction reduction are important factors; thus, numerical analysis should consider the coupling of water and soil, as well as the effect of saturation to obtain more precise failure mechanism. Based on the advanced smoothed particle hydrodynamics (SPH) method, this work proposed a two-phase-coupled SPH model in coordination with a novel constitutive model for unsaturated soils. Then, a triaxial compression test is simulated to check the applicability of the SPH method on the soil phase. After that, the failure test of a dike due to water level-up is discretized and simulated, from which the seepage process, the distribution of maximum shear strain, the slip surface, and pore water pressure are obtained. The two-phase-coupled SPH model is also applied to a slope failure test of heavy rainfall, and the results are compared to the model test. Finally, a dike failure test due to rainfall is analyzed using the proposed SPH model to reproduce the surface infiltration and suction reduction. The proposed SPH model provides several insights of seepage failures and can be a helpful tool for the analysis of dike failures induced by water level-up and rainfall.  相似文献   

2.
Segmental retaining wall (SRW) systems are commonly used in geotechnical practice to stabilize cut and fill slopes. Because of their flexibility, these systems can tolerate minor movements and settlements without incurring damage or crack. Despite these advantages, very few numerical studies of large deformations and post‐failure behavior of SRW systems are found in the current literature. Traditional numerical methods, such as the finite element method, suffer from mesh entanglement, thus are unable to simulate large deformations and flexible behavior of retaining wall blocks in SRW systems. To overcome the above limitations, a novel computational framework based on the smoothed particle hydrodynamics (SPH) method was developed to simulate large deformations and post‐failure behavior of soils and retaining wall blocks in SRW systems. The proposed numerical framework is a hybrid continuum/discontinuum approach that can model soil as an elasto‐plastic material and retaining wall blocks as independent rigid bodies associated with both translational and rotational degrees of freedom. A new contact model is proposed within the SPH framework to simulate the interaction between the soil and the blocks and between the blocks. As an application of the proposed numerical method, a two‐dimensional simulation of an SRW collapse was simulated and compared to experimental results conducted under the same conditions. The results showed that the proposed computational approach provided satisfactory agreement with the experiment. This suggests that the new framework is a promising numerical approach to model SRW systems. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

3.
In this paper, a new approach to applying confining stress to flexible boundaries in the smoothed particle hydrodynamics (SPH) method is developed to facilitate its applications in geomechanics. Unlike the conventional SPH methods that impose confining boundary conditions by creating extra boundary particles, the proposed approach makes use of kernel truncation properties of SPH approximations that occur naturally at free-surface boundaries. Therefore, it does not require extra boundary particles and, as a consequence, can be utilised to apply confining stresses onto any boundary with arbitrary geometry without the need for tracking the curvature change during the computation. This enables more complicated problems that involve moving confining boundaries, such as confining triaxial tests, to be simulated in SPH without difficulties. To further enhance SPH applications in elasto-plastic computations of geomaterials, a robust numerical procedure to implement Mohr-Coulomb plasticity model in SPH is presented for the first time to avoid difficulties associated with corner singularities in Mohr-Coulomb model. The proposed approach was first validated against two-dimensional finite element (FE) solutions for confining biaxial compression tests to demonstrate its predictive capability at small deformation range when FE solutions are still valid. It is then further extended to three-dimensional conditions and utilised to simulate triaxial compression experiments. Simulation results predicted by SPH show good agreement with experiments, FE solutions, and other numerical results available in the literature. This suggests that the proposed approach of imposing confining stress boundaries is promising and can handle complex problems that involve moving confining boundary conditions.  相似文献   

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

5.
It is a well-known fact that addressing hydromechanical problems in saturated soils with the finite element method and equal-order interpolation formulations in displacements and pore pressures produces unstable results. Classically, stabilization has been achieved by increasing the interpolation degree of displacement with respect to pore pressure, hence fulfilling the Babuska–Brezzi condition. However, the use of quadratic elements involves high computational costs. From that point of view, the use of stabilized low-order elements is a more desirable option. Much research has been carried out in different directions in the stabilization of low-order formulations for saturated soils in quasistatic conditions, among others with the technique based on strain field enhancement through internal degrees of freedom. This article presents an alternative displacement–pore pressure formulation for saturated soil dynamics based on the enhancement of the displacement field through incompatible modes.  相似文献   

6.
Liquefaction can result in the damage or collapse of structures during an earthquake and can therefore be a great threat to life and property. Many site investigations of liquefaction disasters are needed to study the large-scale deformation and flow mechanisms of liquefied soils that can be used for performance assessments and infrastructure improvement. To overcome the disadvantages of traditional flow analysis methods for liquefied soils, a soil–water-coupled smoothed particle hydrodynamics (SPH) modeling method was developed to analyze flow in liquefied soils. In the proposed SPH method, water and soil were simulated as different layers, while permeability, porosity, and interaction forces could be combined to model water-saturated porous media. A simple shear test was simulated using the SPH method with an elastic model to verify its application to solid phase materials. Subsequently, the applicability of the proposed SPH modeling method to the simulation of interaction forces between water and soil was verified by a falling-head permeability test. The coupled SPH method produced good simulations for both the simple shear and falling-head permeability tests. Using a fit-for-purpose experimental apparatus, a physical flow model test of liquefied sand has been designed and conducted. To complement the physical test, a numerical simulation has been undertaken based on the soil–water-coupled SPH method. The numerical results correspond well with the physical model test results in observed configurations and velocity vectors. An embankment failure in northern Sweden was selected so that the application of the soil–water-coupled SPH method could be extended to an actual example of liquefaction. The coupled SPH method simulated the embankment failure with the site investigation well. They have also estimated horizontal displacements and velocities, which can be used to greatly improve the seismic safety of structures.  相似文献   

7.
Simulation of large deformation and post‐failure of geomaterial in the framework of smoothed particle hydrodynamics (SPH) are presented in this study. The Drucker–Prager model with associated and non‐associated plastic flow rules is implemented into the SPH code to describe elastic–plastic soil behavior. In contrast to previous work on SPH for solids, where the hydrostatic pressure is often estimated from density by an equation of state, this study proposes to calculate the hydrostatic pressure of soil directly from constitutive models. Results obtained in this paper show that the original SPH method, which has been successfully applied to a vast range of problems, is unable to directly solve elastic–plastic flows of soil because of the so‐called SPH tensile instability. This numerical instability may result in unrealistic fracture and particles clustering in SPH simulation. For non‐cohesive soil, the instability is not serious and can be completely removed by using a tension cracking treatment from soil constitutive model and thereby give realistic soil behavior. However, the serious tensile instability that is found in SPH application for cohesive soil requires a special treatment to overcome this problem. In this paper, an artificial stress method is applied to remove the SPH numerical instability in cohesive soil. A number of numerical tests are carried out to check the capability of SPH in the current application. Numerical results are then compared with experimental and finite element method solutions. The good agreement obtained from these comparisons suggests that SPH can be extended to general geotechnical problems. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

8.
The development of a numerical procedure for the finite element analysis of anchors dynamically penetrating into saturated soils is outlined, highlighting its unique features and capabilities. The mechanical behaviour of saturated porous media is predicted using mixture theory. An algorithm is developed for frictional contact in terms of effective normal stress. The contact formulation is based on a mortar segment-to-segment scheme, which considers the interpolation functions of the contact elements to be of order N, thus overcoming a numerical deficiency of the so-called node-to-segment (NTS) contact algorithm. The nonlinear behaviour of the solid constituent is captured by the Modified Cam Clay soil model. The soil constitutive model is also adapted so as to incorporate the dependence of clay strength on strain rate. An appropriate energy-absorbing boundary is used to eliminate possible wave reflections from the artificial mesh boundaries. To illustrate the use of the proposed computational scheme, simulations of dynamically penetrating anchors are conducted. Results are presented and discussed for the installation phase followed by ‘setup’, i.e., pore pressure dissipation and soil consolidation. The results, in particular, reveal the effects of strain rate on the generation of excess pore pressure, bearing resistance and frictional forces. The setup analyses also illustrate the pattern in which pore pressures are dissipated within the soil domain after installation. Hole closure behind a dynamic projectile is also illustrated by an example.  相似文献   

9.
方志  陈育民  何森凯 《岩土力学》2018,39(5):1851-1857
减饱和法是一种通过减小饱和砂土地基中水的饱和度来处置可液化砂土地基的方法。基于单相流-固耦合理论,将减饱和砂中水-气两相流体等效为单相流体,提出一种可以考虑加载过程中孔隙流体体积模量变化的减饱和砂土静态液化的单相流改进算法,用其进行单调荷载作用下三轴不排水压缩试验数值模拟研究,分析了不同饱和度条件下的减饱和砂土的不排水反应并与饱和砂土三轴不排水试验结果进行对比。研究结果表明,单相流改进算法能够很好地反映减饱和法的抗液化特性。此外,对比不同数值分析方法模拟结果,得出以下结论：采用单相流简化算法分析减饱和砂土的不排水反应时,因为不考虑加载过程中的孔隙流体体积模量变化,所以初始的流体体积模量取值会影响减饱和砂土的强度,初始围压为100 kPa、饱和度为96%的减饱和砂土在单调荷载作用下,气体体积模量取值从100 kPa增加至200 kPa时,减饱和砂试样的峰值偏应力会减小约30%,孔隙压力会增加约40%;通过对比同等条件下的单相流改进算法和两相流算法的应力-应变关系曲线以及饱和度和体积应变变化曲线,两者结果相近,误差在5%以内。综上所述,单相流改进算法是一种较为合理而简洁地模拟减饱和砂土静态液化的计算方法。  相似文献   

10.
The axisymmetric formulation of the governing equations for geomechanics in the framework of smoothed particle hydrodynamics (SPH) is presented in this study. Two forms of SPH discretization for the motion equations, which are labeled as form I and form II, are proposed, and the methods to compute the hoop stress and strain terms including hoop strain rate and the acceleration introduced by the hoop stress are compared. To avoid possible singularity problem near the axis of symmetry, a perfectly smooth contact along with ghost particles are applied to prevent the real particles from overly approaching the axis of symmetry to remove this potential singularity. In addition, the Mohr–Coulomb constitutive model is implemented into the SPH formulation in describing soil behavior. Four numerical tests are carried out to validate and compare the accuracy and stability of the proposed algorithms, and their results are compared with analytical solutions and results from FEM analysis. The performance in these comparisons suggests that SPH II with hoop terms computed through direct hoop method is more stable than the others, and the adoption of contact for the symmetric axis is efficient in eliminating the singularity problem. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

11.
It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

12.
Flow analysis of liquefied soils based on smoothed particle hydrodynamics   总被引：7，自引：6，他引：1  
To overcome the disadvantages of traditional flow analysis methods for liquefied soils that exhibit fluidization and large deformation characteristics, Smoothed particle hydrodynamics (SPH) is adopted in this study to analyze the flow processes of liquefied soils. Bingham model with the use of the Mohr–Coulomb yield criterion, the concepts of equivalent Newtonian viscosity, and the Verlet neighbor list method are introduced into the framework of SPH to build an algorithm for the analysis of flowing liquefied soils. This modeling involves a simulation of physical model test of flowing liquefied soils that can be compared with numerical results. In addition, a shaking table test is selected from the literature for SPH analysis to verify the validation of the SPH method and extend its applications. The SPH simulation can reproduce the flow processes of liquefied soils and constrain estimates of the horizontal displacement, vertical displacement, and velocity of soils after liquefaction. According to the dynamic behaviors of the materials involved, designs can be implemented to improve the seismic safety of structures.  相似文献   

13.
ABSTRACT

Location of failure slip surfaces plays a critical role in landslide risk assessment and mitigation, particularly for unstable slopes, because it is a key input to design of stabilisation measures for unstable slopes and it determines the volume of the sliding soil mass (i.e. landslide consequence). The failure slip surfaces in the numerical analysis (e.g. finite element/different method, FEM/FDM) are often identified using shear strength reduction (SSR) method. A careful examination of FEM results showed that, although the SSR method performs well for stable slopes, it might provide misleading results for unstable slopes. To properly locate failure slip surfaces for unstable slopes, this paper presents a particle-based numerical method called smoothed particle hydrodynamics (SPH), which is mesh-free, immune to the mesh distortion problem in FEM/FDM, and able to directly simulate large deformation of soils that occurs during landslides. A series of slope stability analyses is performed using an in-house SPH programme. Failure slip surfaces are properly identified by SPH for both stable and unstable slopes. Furthermore, because SPH provides a spatial distribution of the post-landslide large displacement of soils, the failure slip surfaces can be identified conveniently using soil displacement. A displacement-based criterion is proposed to locate the failure slip surfaces.  相似文献   

14.
Reservoir depletion results in rock failure, wellbore instability, hydrocarbon production loss, oil sand production, and ground surface subsidence. Specifically, the compaction of carbonate reservoirs with soft rocks often induces large plastic deformation due to rock pore collapse. On the other hand, following the compaction of reservoirs and failure of rock formations, the porosity and permeability of formations will, in general, decrease. These bring a challenge for reservoir simulations because of high nonlinearity of coupled geomechanics and fluid flow fields. In this work, we present a fully implicit, fully coupled, and fully consistent finite element formulation for coupled geomechanics and fluid flow problems with finite deformation and nonlinear flow models. The Pelessone smooth cap plasticity model, an important material model to capture rock compaction behavior and a challenging material model for implicit numerical formulations, is incorporated in the proposed formulation. Furthermore, a stress-dependent permeability model is taken into account in the formulation. A co-rotational framework is adopted for finite deformation, and an implicit material integrator for cap plasticity models is consistently derived. Furthermore, the coupled field equations are consistently linearized including nonlinear flow models. The physical theories, nonlinear material and flow models, and numerical formulations are the focus of part I of this work. In part II, we verify the proposed numerical framework and demonstrate the performance of our numerical formulation using several numerical examples including a field reservoir with soft rocks undergoing serious compaction.  相似文献   

15.
A new formulation of the element‐free Galerkin (EFG) method is developed for solving coupled hydro‐mechanical problems. The numerical approach is based on solving the two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Spatial variables in the weak form, i.e. displacement increment and pore water pressure increment, are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on a penalty method. Numerical stability of the developed formulation is examined in order to achieve appropriate accuracy of the EFG solution for coupled hydro‐mechanical problems. Examples are studied and compared with closed‐form or finite element method solutions to demonstrate the validity of the developed model and its capabilities. The results indicate that the EFG method is capable of handling coupled problems in saturated porous media and can predict well both the soil deformation and variation of pore water pressure over time. Some guidelines are proposed to guarantee the accuracy of the EFG solution for coupled hydro‐mechanical problems. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

16.
Risk analysis of existing slopes in catchment areas requires quantification of their stability. This quantification becomes particularly difficult when dealing with larger areas under 3D conditions and including saturated and unsaturated water flow. This paper proposes the use of an effective numerical procedure to solve three-dimensional slope stability problems in large areas subjected to pore pressure effects. This numerical approach, numerical limit analysis, utilizes the finite element method and mathematical programming techniques. Mathematical programming is needed because the basic plasticity theorems for limit analysis can be cast as optimization problems. The generated optimization problem is formulated under a second-order cone programming framework, which is known to solve large-scale problems with great computational efficiency. The main objective of this work was to determine the slope safety factor and the collapse mechanism of soils governed by the Drucker–Prager yield criterion for large-scale 3D problems including pore pressure effects. This approach is applied to an experimental catchment in the Oregon Coast Range that failed after an intense rainfall. The results were compared with a previous stability analysis of the area available in the literature that used a novel 3D limit equilibrium method.  相似文献   

17.
The smoothed particle hydrodynamics (SPH) method was recently extended to simulate granular materials by the authors and demonstrated to be a powerful continuum numerical method to deal with the post-flow behaviour of granular materials. However, most existing SPH simulations of granular flows suffer from significant stress oscillation during the post-failure process, despite the use of an artificial viscosity to damp out stress fluctuation. In this paper, a new SPH approach combining viscous damping with stress/strain regularisation is proposed for simulations of granular flows. It is shown that the proposed SPH algorithm can improve the overall accuracy of the SPH performance by accurately predicting the smooth stress distribution during the post-failure process. It can also effectively remove the stress oscillation issue in the standard SPH model without having to use the standard SPH artificial viscosity that requires unphysical parameters. The predictions by the proposed SPH approach show very good agreement with experimental and numerical results reported in the literature. This suggests that the proposed method could be considered as a promising continuum alternative for simulations of granular flows.  相似文献   

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

19.
Development of robust numerical solutions for poro‐elasticity is an important and timely issue in modern computational geomechanics. Recently, research in this area has seen a surge in activity, not only because of increased interest in coupled problems relevant to the petroleum industry, but also due to emerging applications of poro‐elasticity for modelling problems in biomedical engineering and materials science. In this paper, an original mixed least‐squares method for solving Biot consolidation problems is developed. The solution is obtained via minimization of a least‐squares functional, based upon the equations of equilibrium, the equations of continuity and weak forms of the constitutive relationships for elasticity and Darcy flow. The formulation involves four separate categories of unknowns: displacements, stresses, fluid pressures and velocities. Each of these unknowns is approximated by linear continuous functions. The mathematical formulation is implemented in an original computer program, written from scratch and using object‐oriented logic. The performance of the method is tested on one‐ and two‐dimensional classical problems in poro‐elasticity. The numerical experiments suggest the same rates of convergence for all four types of variables, when the same interpolation spaces are used. The continuous linear triangles show the same rates of convergence for both compressible and entirely incompressible elastic solids. This mixed formulation results in non‐oscillating fluid pressures over entire domain for different moments of time. The method appears to be naturally stable, without any need of additional stabilization terms with mesh‐dependent parameters. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

20.
The chemical‐dissolution front propagation problem exists ubiquitously in many scientific and engineering fields. To solve this problem, it is necessary to deal with a coupled system between porosity, pore‐fluid pressure and reactive chemical‐species transport in fluid‐saturated porous media. Because there was confusion between the average linear velocity and the Darcy velocity in the previous study, the governing equations and related solutions of the problem are re‐derived to correct this confusion in this paper. Owing to the morphological instability of a chemical‐dissolution front, a numerical procedure, which is a combination of the finite element and finite difference methods, is also proposed to solve this problem. In order to verify the proposed numerical procedure, a set of analytical solutions has been derived for a benchmark problem under a special condition where the ratio of the equilibrium concentration to the solid molar density of the concerned chemical species is very small. Not only can the derived analytical solutions be used to verify any numerical method before it is used to solve this kind of chemical‐dissolution front propagation problem but they can also be used to understand the fundamental mechanisms behind the morphological instability of a chemical‐dissolution front during its propagation within fluid‐saturated porous media. The related numerical examples have demonstrated the usefulness and applicability of the proposed numerical procedure for dealing with the chemical‐dissolution front instability problem within a fluid‐saturated porous medium. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号