基于三维并行SPH模型的土体流滑特性研究

土体的大变形流滑导致了许多地质灾害的发生,对人们的生命财产安全构成了极大的威胁,因此越来越多的研究开始关注土体的大变形流滑特性。其中,光滑粒子流体动力学（SPH）方法是常用的模拟方法之一,但SPH方法的粒子特性导致其计算时间过长,影响了在工程地质领域的进一步应用。对此,本研究基于SPH方法的基本原理、非牛顿流体理论和等效黏度概念,提出了适用于土体大变形流滑分析的三维SPH仿真模型。结合OpenMP并行计算原理,实现了SPH算法的并行优化。在此基础上,对土体流滑模型试验进行了二维和三维分析,得到了滑动距离、滑动冲击力和冲击力峰值等动力学参数,分析了计算维数和边界条件对流滑特性的影响机制。通过不同线程数下计算时间的对比,获得了计算效率随线程数的变化规律。结果证明了本文的OpenMP并行优化具有较高的计算效率,显著降低了三维SPH模拟的计算耗时,对工程地质数值方法的效率提升具有重要的借鉴意义。     

Flow analysis of liquefied soils based on smoothed particle hydrodynamics

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.     

Lagrangian meshfree particles method (SPH) for large deformation and failure flows of geomaterial using elastic–plastic soil constitutive model

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.     

Numerical simulation of flow processes in liquefied soils using a soil–water-coupled smoothed particle hydrodynamics method

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.     

A comprehensive study on the parameters setting in smoothed particle hydrodynamics (SPH) method applied to hydrodynamics problems

Smoothed particle hydrodynamics (SPH) is a meshfree, Lagrangian particle method which has advantages in handling solids with extremely large deformation. Like any other numerical methods, cares must be taken to ensure its desirable accuracy and stability through considering several correction techniques in calculation. The selection of values for parameters in those correction approaches is a key step in SPH simulation, which is always difficult for new beginners to deal well with effectively. This paper examines the common inconsistency and instability problems in SPH method and studies its computational efficiency when applied to hydrodynamics problems with material strength like soil column collapse. We analyzed in detail how the correction techniques mitigate these inconsistency and instability problems. Also, the numerical testing results associate with different values for the parameters used in the correction techniques are provided for better understanding the influence of these parameters and for finding out the desirable values. It is found that (1) the SPH method is easily subjected to an inconsistency problem in the boundary area due to the boundary deficiency, and it can be treated well by adopting “virtual particles” contributing to the particle summations. (2) The numerical oscillation in SPH simulation can be mitigated effectively by artificial viscosity with the suggested parameter values. (3) The tension cracking treatment, artificial viscosity and artificial stress work well in removing the tensile instability problem in SPH method. In addition, the nearest neighboring particle searching (NNPS) algorithm, spacing ratio, smoothing length and time step influence the efficiency and accuracy of SPH method significantly. It is shown that SPH method with suggested parameters values can produce a very good result compared with the experimental result.     

Simulation of earthquake‐induced slope deformation using SPH method

This paper presents the development, calibration, and validation of a smoothed particle hydrodynamics (SPH) model for the simulation of seismically induced slope deformation under undrained condition. A constitutive model that combines the isotropic strain softening viscoplasticity and the modified Kondner and Zelasko rule is developed and implemented into SPH formulations. The developed SPH model accounts for the effects of wave propagation in the sliding mass, cyclic nonlinear behavior of soil, and progressive reduction in shear strength during sliding, which are not explicitly considered in various Newmark‐type analyses widely used in the current research and practice in geotechnical earthquake engineering. Soil parameters needed for the developed model can be calibrated using typical laboratory shear strength tests, and experimental or empirical shear modulus reduction curve and damping curve. The strain‐rate effects on soil strength are considered. The developed SPH model is validated against a readily available and well‐documented model slope test on a shaking table. The model simulated slope failure mode, acceleration response spectra, and slope deformations are in excellent agreement with the experimental data. It is thus suggested that the developed SPH model may be utilized to reliably simulate earthquake‐induced slope deformations. This paper also indicates that if implemented with appropriate constitutive models, SPH method can be used to model large‐deformation problems with high fidelity. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical analysis on seepage failures of dike due to water level-up and rainfall using a water–soil-coupled smoothed particle hydrodynamics model

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.     

Parallel simulations of pore-scale flow through porous media

Smoothed particle hydrodynamics (SPH) is a versatile technique which can be applied to single and multiphase flow through porous media. The versatility of SPH is offset by its computational expense which limits the practicability of SPH for large problems involving low Reynolds number flow. A parallel pore-scale numerical model based on SPH is described for modeling flow phenomena in porous media. Aspects of SPH which complicate parallelization are emphasized. The speed of the method is demonstrated to be proportional to the number of processors for test cases where load balance was achieved. The parallel algorithm permits the application of SPH to more complicated porous media problems than previously considered. For such problems, best performance is achieved when several soil grains are simulated by each processor. Finally, future applications of the method and possible extensions are discussed.     

Stable axisymmetric SPH formulation with no axis singularity

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.     

A novel computational approach for large deformation and post‐failure analyses of segmental retaining wall systems

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.     

水势理论在四水转化研究中应用时代特征与趋势:纪念ZFP方法引进我国30周年

20世纪80年代零通量面方法在我国应用中解决了"四水"转化研究中参数不确定性带来的问题,发现在降水入渗补给地下水过程中土壤总水势梯度大于1.0cm H2O/cm,且逐渐降低,流入、流出被监测土层的水量相等时土壤总水势梯度趋近于1.0cm H2O/cm。20世纪90年代,利用土壤水势与含水量之间量化关系,指导了农业节水灌溉,提出在灌溉过程中土壤水势梯度等于1.0cm H2O/cm的时间持续愈长,表明过剩灌溉而浪费的水量愈多的认识。进入21世纪以来,水势理论较广泛地用来解决土壤水盐分运移数值模拟与入渗模型中水文地质参数问题和降水入渗土壤水势运移微观机理研究,并发现表聚型、中聚型和底聚型土壤盐分剖面的水动力学特征。通过土壤水动力场调控改变土壤水盐(养分或污染物)运移是未来重要研究方向。     

地震动加速度幅值对土坡失稳后滑坡体大小及致灾范围的影响

地震滑坡的致灾范围是判断滑坡能否会对已有建构筑物造成损失、确定预警疏散范围的重要依据,因此对地震土坡破坏后的滑坡体大小和致灾范围进行研究具有重要的意义.本研究基于SPH动力分析方法,结合弹塑性本构模型和固体力学控制方程建立了地震土坡破坏的动力分析模型;通过设置振动边界粒子和自由场边界粒子,实现了地震动加速度的施加以及自...     

An improved SPH method for saturated soils and its application to investigate the mechanisms of embankment failure: Case of hydrostatic pore‐water pressure

The method of smoothed particle hydrodynamics (SPH) has recently been applied to computational geomechanics and has been shown to be a powerful alternative to the standard numerical method, that is, the finite element method, for handling large deformation and post‐failure of geomaterials. However, very few studies apply the SPH method to model saturated or submerged soil problems. Our recent studies of this matter revealed that significant errors may be made if the gradient of the pore‐water pressure is handled using the standard SPH formulation. To overcome this problem and to enhance the SPH applications to computational geomechanics, this article proposes a general SPH formulation, which can be applied straightforwardly to dry and saturated soils. For simplicity, the current work assumes hydrostatic pore‐water pressure. It is shown that the proposed formulation can remove the numerical error mentioned earlier. Moreover, this formulation automatically satisfies the dynamic boundary conditions at a submerged ground surface, thereby saving computational cost. Discussions on the applications of the standard and new SPH formulations are also given through some numerical tests. Furthermore, techniques to obtain the correct SPH solution are also proposed and discussed throughout. As an application of the proposed method, the effect of the dilatancy angle on the failure mechanism of a two‐sided embankment subjected to a high groundwater table is presented and compared with that of other solutions. Finally, the proposed formulation can be considered a basic formulation for further developments of SPH for saturated soils. Copyright © 2011 John Wiley & Sons, Ltd.     

Run-out analysis of flow-like landslides triggered by the Ms 8.0 2008 Wenchuan earthquake using smoothed particle hydrodynamics

Flow-like landslides have caused significant damage and casualties worldwide. However, studying such phenomena with traditional simulation methods is made difficult by their complex fluidization characteristics. In this paper, we use smoothed-particle hydrodynamics (SPH) for the run-out analysis of flow-like landslides. Compared with conventional methods, the proposed SPH modeling technique is the combination of a Bingham flow model and Navier?CStokes equations in the framework of computational fluid dynamics. At first, two benchmark problems of dam break and granular flow are simulated and verified to evaluate the accuracy of the SPH model. Then, run-out analyses are performed for flow-like landslides triggered by the Ms 8.0 Wenchuan earthquake that occurred on 12 May 2008 in Sichuan Province, China. Run-out analyses of the Tangjiashan, Wangjiayan, and Donghekou landslides are conducted by the application of SPH models to real flow-like landslides. All simulations show good agreement with characteristics of flow-like landslides observed in the field. We have found that numerical modeling can capture the fundamental dynamic behavior of these flow-like landslides and produce preliminary results for hazard assessment and site selection for reconstruction in earthquake-prone areas.     

生态脆弱矿区土层中导水裂缝带发育高度研究

为指导生态脆弱区"保水采煤",掌握土层中导水裂缝带发育高度,采用微电阻率扫描成像技术、物理模拟、数值模拟、钻孔冲洗液消耗量观测等手段综合研究,结果表明:物理模拟、数值模拟辨别出的裂采比为28.1~29.1,微电阻率扫描成像探测出裂采比28.3~28.5,微电阻率扫描成像识别与其他手段所得结果吻合。微电阻率扫描成像技术可推广应用于其他矿区土层中导水裂缝带发育高度的探测。      

SPH-DEM coupling method based on GPU and its application to the landslide tsunami. Part I: method and validation

Landslide-induced tsunami is a complex fluid–solid coupling process that plays a crucial role in the study of a disaster chain. To simulate the coupling behaviors between the fluid and solid, a graphics processing unit-based coupled smoothed particle hydrodynamics (SPH)-discrete element method (DEM) code is developed. A series of numerical tests, which are based on the laboratory test by Koshizuka et al. (Particle method for calculating splashing of incompressible viscous fluid, 1995) and Kleefsman et al. (J Comput Phys 206:363–393, 2005), are carried out to study the influence of the parameters, and to verify the accuracy of the developed SPH code. To ensure accurate results of the SPH simulation, the values for the diffusion term, particle resolution (1/25 characteristic length), and smoothing length (1.2 times of particle interval) are suggested. The ratio of the SPH particle size and the DEM particle’s diameter influences the accuracy of the coupling simulation between solid particles and water. For the coupling simulation of a single particle or a loose particle assembly (not contact each other) with fluid, this ratio should be smaller than 1/20; for a dense particle assembly, a ratio of smaller than 1/6 will be good.

     

A smoothed particle hydrodynamics modelling of soil–water mixing and resulting changes in average strength

Soil–water interaction is a pivotal process in many underwater geohazards such as underwater landslides where soil sediments gradually evolve into turbidity currents after interactions with ambient water. Due to the large deformations, multiphase interactions and phase changes this involves, investigations from numerical modelling of the transition process have been limited so far. This study explores a simple numerical replication of such soil–water mixing with respect to changes in average strength using smoothed particle hydrodynamics (SPH). A uniform viscoplastic model is used for both the solid-like and fluid-like SPH particles. The proposed numerical solution scheme is verified by single-phase dam break tests and multiphase simple shear tests. SPH combinations of solid-like and fluid-like particles can replicate the clay–water mixture as long as the liquidity index of the solid-like particles is larger than unity. The proposed numerical scheme is shown to capture key features of an underwater landslide such as hydroplaning, water entrainment and wave generation and thus shows promise as a tool to simulate the whole process of subaquatic geohazards involving solid–fluid transition during mass transport.     

Frictional contact algorithms in SPH for the simulation of soil–structure interaction

Simulation of frictional contact between soils and rigid or deformable structure in the framework of smoothed particle hydrodynamics (SPH) is presented in this study. Two algorithms are implemented into the SPH code to describe contact behavior, where the contact forces are calculated using the law of conservation of momentum based on ideal plastic collision or using the criteria of partial penetrating. In both algorithms, the problem of boundary deficiency inherited from SPH is properly handled so that the particles located at contact boundary can have precise acceleration, which is critical for contact detection. And the movement and rotation of the rigid structure are taken into account so that it is easy to simulate the process of pile driving or movement of a retaining wall in geotechnical engineering analysis. Furthermore, the capability of modeling deformability of a structure during frictional contact simulations broadens the fields of SPH application. In contrast to previous work dealing with contact in SPH, which usually use particle‐to‐particle contact or ignoring sliding between particles and solid structure, the method proposed here is more efficient and accurate, and it is suitable to simulate interaction between soft materials and rigid or deformable structures, which are very common in geotechnical engineering. A number of numerical tests are carried out to verify the accuracy and stability of the proposed algorithms, and their results are compared with analytical solutions or results from finite element method analysis. Good agreement obtained from these comparisons suggests that the proposed algorithms are robust and can be applied to extend the capability of SPH in solving geotechnical problems. Copyright © 2013 John Wiley & Sons, Ltd.     

灵新煤矿导水裂隙带发育高度数值模拟研究

煤矿的开采利用给国民经济发展带来巨大的收益,但也引发了许多环境地质问题,特别在煤层开采过程中,煤层上覆基岩变形破坏形成的裂隙通道极易发生矿井涌（突）水事故,时刻威胁着井下工人的生命安全。本文以灵新煤矿051505工作面为研究对象,利用Flac3D数值模拟软件,对14号主采煤层上覆基岩导水裂隙带高度进行了模拟研究。模拟结果表明：当煤层开采厚度为2.5 m时,导水裂隙带发育最大高度为59.5 m。同时选取经验公式法对导水裂隙带高度进行了计算。最终通过钻孔实测法得到的结果与前两种方法对比分析,数值模拟结果与钻孔实测结果基本吻合,认为数值模拟方法能够高效、简单、合理达到预测导水裂隙带高度的目的,也为同类矿井安全、绿色生产提供一定的借鉴。     

An elasto‐plastic three phase model for partially saturated soil for the finite element simulation of compressed air support in tunnelling

This paper presents a fully coupled finite element formulation for partially saturated soil as a triphasic porous material, which has been developed for the simulation of shield tunnelling with heading face support using compressed air. While for many numerical simulations in geotechnics use of a two‐phase soil model is sufficient, the simulation of compressed air support demands the use of a three‐phase model with the consideration of air as a separate phase. A multiphase model for soft soils is developed, in which the individual constituents of the soil—the soil skeleton, the fluid and the gaseous phase—and their interactions are considered. The triphasic model is formulated within the framework of the theory of porous media, based upon balance equations and constitutive relations for the soil constituents and their mixture. An elasto‐plastic, cam–clay type model is extended to partially saturated soil conditions by incorporating capillary pressure according to the Barcelona basic model. The hydraulic properties of the soil are described via DARCY 's law and the soil–water characteristic curve after VAN GENUCHTEN . Water is modelled as an incompressible and air as a compressible phase. The model is validated by means of selected benchmark problems. The applicability of the model to geotechnical problems is demonstrated by results from the simulation of a compressed air intervention in shield tunnelling. Copyright © 2009 John Wiley & Sons, Ltd.