共查询到20条相似文献，搜索用时 531 毫秒
1.
We present a Lagrangian formulation for simulating the continuum hydrodynamics of dry granular flows based on multiplicative elastoplasticity theory for finite deformation calculations. The formulation is implemented within the smoothed particle hydrodynamics (SPH) method along with a variant of the usual dynamic boundary condition. Three benchmark simulations on dry sands are presented to validate the model: (a) a set of plane strain collapse tests, (b) a set of 3D collapse tests, and (c) a plane strain simulation of the impact force generated by granular flow on a rigid wall. Comparison with experimental results suggests that the formulation is sufficiently robust and accurate to model the continuum hydrodynamics of dry granular flows in a laboratory setting. Results of the simulations suggest the potential of the formulation for modeling more complex, field-scale scenarios characterized by more elaborate geometry and multi-physical processes. To the authors’ knowledge, this is the first time the multiplicative plasticity approach has been applied to granular flows in the context of the SPH method. 相似文献
2.
Fávero Neto Alomir H. Askarinejad Amin Springman Sarah M. Borja Ronaldo I. 《Acta Geotechnica》2020,15(10):2757-2777
We present an updated Lagrangian continuum particle method based on smoothed particle hydrodynamics (SPH) for simulating debris flow on an instrumented test slope. The site is a deforested area near the village of Ruedlingen, a community in the canton of Schaffhausen in Switzerland. Artificial rainfall experiments were conducted on the slope that led to failure of the sediment in the form of a debris flow. We develop a 3D mechanistic model for this test slope and conduct numerical simulations of the flow kinematics using an SPH formulation that captures large deformation, material nonlinearity, and the complex post-failure movement of the sediment. Two main simulations explore the impact of changes in the mechanical properties of the sediment on the ensuing kinematics of the flow. The first simulation models the sediment as a granular homogeneous material, while the second simulation models the sediment as a heterogeneous material with spatially varying cohesion. The variable cohesion is meant to represent the effects of root reinforcement from vegetation. By comparing the numerical solutions with the observed failure surfaces and final free-surface geometries of the debris deposit, as well as with the observed flow velocity, flow duration, and hot spots of strain concentration, we provide insights into the accuracy and robustness of the SPH framework for modeling debris flows.
相似文献3.
For research on granular materials, establishing a method to calculate continuum strain from particle displacements is necessary for understanding the material behaviour at macro-level and developing continuum constitutive models. Existing methods are generally based on constructing a mesh or background grid to calculate strain from particle motions. These methods offer rigorous ways to measure strain for granular materials; however, they suffer from several problems such as mesh distortion and lacking grid-to-particle strain mapping procedure, which hinders their capability of calculating strain accumulation during large deformation processes of granular media. To address this issue, this study proposes a new strain calculation method for discrete element simulations of granular materials. This method describes a particle assembly as an equivalent continuum system of material points, each of which corresponds to a particle centre and represents a continuous region with its initial volume/area presumably equal to the volume/area of Voronoi cells generated in accordance with the particle assembly configuration. Smooth Particle Hydrodynamics (SPH) interpolation functions are then employed to calculate strain for these material points. This SPH-based method does not require any mesh or background grid for computation, leading to advantages in calculating strain accumulation under large deformation. Simulations of granular materials in both uniform and heterogeneous gradations were carried out, and strain results obtained by the proposed method indicate good agreements with analytical and numerical solutions. This demonstrates its potential for strain calculations in discrete element simulations of granular materials involving large deformations and/or large displacements. 相似文献
4.
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. 相似文献
5.
土体的大变形流滑导致了许多地质灾害的发生,对人们的生命财产安全构成了极大的威胁,因此越来越多的研究开始关注土体的大变形流滑特性。其中,光滑粒子流体动力学（SPH）方法是常用的模拟方法之一,但SPH方法的粒子特性导致其计算时间过长,影响了在工程地质领域的进一步应用。对此,本研究基于SPH方法的基本原理、非牛顿流体理论和等效黏度概念,提出了适用于土体大变形流滑分析的三维SPH仿真模型。结合OpenMP并行计算原理,实现了SPH算法的并行优化。在此基础上,对土体流滑模型试验进行了二维和三维分析,得到了滑动距离、滑动冲击力和冲击力峰值等动力学参数,分析了计算维数和边界条件对流滑特性的影响机制。通过不同线程数下计算时间的对比,获得了计算效率随线程数的变化规律。结果证明了本文的OpenMP并行优化具有较高的计算效率,显著降低了三维SPH模拟的计算耗时,对工程地质数值方法的效率提升具有重要的借鉴意义。 相似文献
6.
The paper presents Cauchy stress tensor computation over parallel grids of message passing interface (MPI) parallel three-dimensional (3D) discrete element method (DEM) simulations of granular materials, considering spherical and nonspherical particles. The stress tensor computation is studied for quasi-static and dynamic conditions, and its resulting symmetry or asymmetry is discussed within the context of classical continuum mechanics (CCM), granular materials mechanics (GMM), and micropolar continuum mechanics (MCM). The average Cauchy stress tensor computation follows Bagi's and Nicot's formulations and is verified within MPI parallel 3D DEM simulations involving dynamically adaptive compute grids. These grids allow calculation of temporal and spatial distributions of stress across granular materials under static and dynamic conditions. The vertical stress component in gravitationally deposited particle assemblies exhibits nonuniform spatial distributions under static equilibrium, and its zone of maximum value changes during the process of gravitational pluviation and collapse. These phenomena reveal a microstructural effect on stress distribution within granular materials that is attributed to their discrete particulate nature (particle size, shape, gradation, boundary conditions, etc). 相似文献
7.
Shaohan Zhao Ha H. Bui Vincent Lemiale Giang D. Nguyen Felix Darve 《国际地质力学数值与分析法杂志》2019,43(5):1005-1031
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. 相似文献
8.
One of the purposes in this study is to develop a modified micromorphic continuum model for granular materials on the basis of a micromechanics approach. A symmetric curvature tensor is proposed in this model, and a symmetric couple stress tensor is derived conjugating the symmetric curvature tensor. In addition, a correct derivation is presented to obtain the symmetric stress tensor conjugated with the symmetric strain tensor. The modified model provides a complete deformation mode for granular materials by considering the decomposition for motions (displacement and rotation) of particles. Consequently, the macroscopic constitutive relationships and constitutive moduli are derived in expressions of the microstructural information. Furthermore, the balance equations and boundary conditions are obtained for the modified micromorphic model. By considering the extended Drucker-Prager yield criterion, the micromorphic elastoplastic model is developed. Another purpose of this study is to derive the finite element formulation for the developed micromorphic elastoplastic model. Based on the ABAQUS user element (UEL) interface, numerical simulations investigated the load-displacement relationship and the strain localization behavior of granular materials and investigated the influence of microscopic parameters in the micromorphic model on these macroscopic mechanical responses. Numerical results illustrate the presented model's capability of simulating the strain-softening and strain localization behaviors, and the capability of considering the influence of microstructural information on the macroscopic mechanical behaviors of granular materials. 相似文献
9.
An assessment of the material point method for modelling large scale run-out processes in landslides 总被引：1，自引：1，他引：0
This paper demonstrates the predictive capabilities of a numerical model based on continuum mechanics for the simulation of run-out processes during landslides. It assesses a particle-based method that takes advantage of a double Lagrangian-Eulerian discretization and known as the material point method (MPM). Attention is given to the post-failure behaviour and, in particular, to the computation of important quantities such as run-out distance, maximum velocity and energy release. The MPM is a step forward in computational solid mechanics and has the potential to simulate large deformations such as those occurring during landslides. A validation is conducted based on simulations of two case studies of different scales, namely the Tokai-Hokuriku expressway failure in Japan and the Vajont landslide in Italy. The results show a very good agreement with field and other numerical observations. 相似文献
10.
A new approach for calculating strain for particulate media 总被引：1，自引：0，他引：1
Discrete element modelling is a viable alternative to conventional continuum‐based analysis for analysing problems involving localized deformations of particulate media. However, to aid in the interpretation of the results, it is useful to express the results of discrete element analyses in terms of the continuum parameters of stress and strain. A number of homogenization methods have been proposed to calculate strain in discrete systems; however, two significant limitations of these methods remain. First, none of these methods incorporate particle rotation effects satisfactorily, although significant particle rotation occurs in shear bands in both physical tests and numerical simulations of granular materials. Additionally, observations of the particle displacement fields in shear bands in granular materials indicate that the displacements within the localizations are erratic. Consequently, existing linear, local interpolation approaches produce substantial variations in the strain values calculated in adjacent elements in the region of localization, hindering clear visualization of the strain localization as it evolves. A new method of domain discretization for calculating strain is proposed. This method is capable of capturing particle rotation and employs a non‐local meshfree interpolation procedure capable of smoothing the erratic displacements in strain localizations, which better defines their evolution. The proposed method is validated for problems involving both two and three dimensions. A number of methods are compared with the proposed method and pertinent insights are made. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
11.
对基于平均场理论的颗粒材料多尺度方法实施过程中的基本理论问题进行了探讨。公式推导表明,考虑微观颗粒间滚动矩相互作用时宏观需要采用Cosserat类连续体模型。基于平均场理论,给出了离散颗粒表征元等效连续介质应力、偶应力表达式。基于Hill-Mandel宏观均质化条件,给出了离散颗粒微结构波动场的边界条件,并且讨论了内部长度尺度参数的选择问题。该工作是颗粒材料多尺度计算的理论部分研究。 相似文献
12.
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. 相似文献
13.
This paper is devoted to numerical analysis of strength and deformation of cohesive granular materials. The emphasis is put on the study of effects of confining pressure and loading path. To this end, the three-dimensional discrete element method is used. A nonlinear failure criterion for inter-granular interface bonding is proposed, and it is able to account for both tensile and shear failure for a large range of normal stress. This criterion is implemented in the particles flow code. The proposed failure model is calibrated from triaxial compression tests performed on representative sandstone. Numerical results are in good agreement with experimental data. In particular, the effect of confining pressure on compressive strength and failure pattern is well described by the proposed model. Furthermore, numerical predictions are studied, respectively, for compression and extension tests with a constant mean stress. It is shown that the failure strength and deformation process are clearly affected by loading path. Finally, a series of numerical simulations are performed on cubic samples with three independent principal stresses. It is found that the strength and failure mode are strongly influenced by the intermediate principal stress.
相似文献14.
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. 相似文献
15.
F. Molenkamp 《国际地质力学数值与分析法杂志》1986,10(2):151-176
In this paper it is shown that the Jaumann stress rate is not accurate for moderate deformations with deviatoric strains of more than 10 per cent; in such cases the objective stress rate is better calculated using the material rotation rate as calculated by polar decomposition. For very large deformations of granular materials this approach becomes invalid too because then the fabric tends to rotate together with principal stresses. Apparently in general the fabric rotation of a granular material cannot be described by continuum mechanics alone. For the numerical calculation of the material rotation for moderate deformations a numerical subroutine, based on polar decomposition, has been developed. 相似文献
16.
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. 相似文献
17.
Equivalent continuum strain calculations based on 3D particle kinematic measurements of sand 
下载免费PDF全文
下载免费PDF全文 Siavash Amirrahmat Khalid A. Alshibli Maha F. Jarrar Boning Zhang Richard A. Regueiro 《国际地质力学数值与分析法杂志》2018,42(8):999-1015
Constitutive modeling of granular materials has been a subject of extensive research for many years. While the calculation of the Cauchy stress tensor using the discrete element method has been well established in the literature, the formulation and interpretation of the strain tensor are not as well documented. According to Bagi, 1 researchers mostly adopt well‐known continuum or discrete microstructural approaches to calculate strains within granular materials. However, neither of the 2 approaches can fully capture the behavior of granular materials. They are considered complementary to each other where each has its own strengths and limitations in solving granular‐mechanics problems. Zhang and Regueiro 2 proposed an equivalent continuum approach to calculating finite strain measures at the local level in granular materials subjected to large deformations. They used three‐dimensional discrete element method results to compare the proposed strains measures. This paper presents an experimental application of the Zhang and Regueiro 2 approach using three‐dimensional synchrotron microcomputed tomography images of a sheared Ottawa sand specimen. Invariant Eulerian finite strain measures were calculated for representative element volumes within the specimen. The spatial maps of Eulerian octahedral shear and volumetric strain were used to identify zones of intense shearing within the specimen and compared well with maps of incremental particle translation and rotation for the same specimen. The local Eulerian volumetric strain was compared to the global volumetric strains, which also can be considered as an averaging of all local Eulerian volumetric strains. 相似文献
18.
存在临界状态是颗粒材料的一个重要特性。基于孔隙胞元的颗粒离散元方法对二维颗粒体进行双轴加载数值试验,在详细分析数值模拟结果的基础上,从微观几何组构的角度揭示了临界状态的存在机制。基于剪胀性原理,提出了以接触价键表征的微观临界状态理论模型,得到了接触价键与塑性剪切应变的关系表达式,理论模型的结果和二维离散元数值模拟得到的结果吻合较好。通过比较不同情况下数值结果和理论模型中的参数,得到以下结论：表征微观临界状态的参数（临界接触价键和达到临界状态所需要的塑性剪切应变）依赖于颗粒体的微观特性,如颗粒形状、表面摩擦性质、颗粒体的围压和初始孔隙比。 相似文献
19.
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. 相似文献
20.
Run-out analysis of flow-like landslides triggered by the Ms 8.0 2008 Wenchuan earthquake using smoothed particle hydrodynamics 总被引：8，自引：5，他引：3
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. 相似文献