共查询到20条相似文献，搜索用时 31 毫秒
1.
Piping, as one of the critical patterns of internal erosion, has been reported as a major cause for failures of embankment dams and levees. The fundamental mechanism of piping was traditionally investigated through experimental trials and simplified theoretical methods in macroscale. Nevertheless, the initiation and progressive evolution of piping is a microscale phenomenon in its essence. The current understanding of the micro-mechanism of piping erosion is limited due to a lack of quantitative analysis and visualized evidence. And in fact, seepage flows can affect the soil fabrics and the development of contact forces between particles. But how these fabrics and contact forces evolve under a critical hydraulic gradient is still not fully understood. In this paper, the detailed process of piping erosion is investigated by using a coupled computational fluid dynamics and discrete element method (CFD–DEM) approach. The treatment of soil–flow interactions in CFD–DEM is explained by exchanging the momentum between the two phases. During the simulation, the piping erosion process is initiated by incrementally ascending differential water head across the soil samples. The three main stages of piping erosion (initial movement, continuation of erosion and total heave) can be identified from monitoring the particle velocity and positions. In addition, the evolution of contact force, hydraulic force, coordination number and void fraction is inspected to provide insight into the micro-mechanism of piping erosion. Two cases are simulated, one with a uniform particle size and a relatively uniform porosity distribution and the other with specific particle size and porosity distributions. An interesting finding from this study is that piping does not always initiate from the free surface and the evolution of piping depends heavily on the particle size and porosity distribution. 相似文献
2.
为了充分发挥有限元与颗粒离散元各自的优势,提出了一种由有限元转化为颗粒流的方法。数值模型首先用较粗的有限元网格进行离散,并在单元上引入连续介质本构模型。力学计算开始后,实时跟踪各单元的应力状态。一旦某单元的应力满足Mohr-Coulomb准则或最大拉应力准则,删除该单元,同时创建具有一定数目、随机分布且微嵌入的颗粒簇。其后,该单元所在区域的非连续变形及失稳断裂由颗粒簇演化获得。各颗粒的质量、材料参数、速度、位移、接触力等信息根据插值从有限元单元中继承。为了实现有限元与颗粒流接触面的耦合计算,引入了点-棱（二维）及点-面（三维）接触模型,通过法向及切向弹簧实现接触力的计算。颗粒球与有限元板的碰撞分析、单轴压缩、岩石切割等案例展示了上述方法的精确性及合理性。 相似文献
3.
采用二维梁-颗粒模型BPM2D（beam-particle model in two dimensions）模拟了刚性弹丸侵彻无钢筋混凝土的过程。离散元法（DEM）和有限元法（FEM）等数值计算方法各有其优势,同时也都存在不足之处：离散元法适于处理由连续介质向非连续介质转化的破坏问题,但对于连续体计算结果精度不高;有限元法适于预测材料破坏的区域,但难以直接用于计算脆性材料破坏过程,因此将两种方法结合可以形成一种较好的混合模型。梁-颗粒模型BPM2D是基于离散元法,结合有限元法开发的二维数值计算模型,采用3种类型梁单元形成混凝土数值试样,每种类型梁单元的力学性质均按韦伯（Weibull）分布随机赋值,以模拟混凝土细观结构的非均匀性,同时梁单元的强度随应变率不同而变化。利用此模型分析了弹丸侵彻下混凝土的破坏过程,并给出侵彻过程弹丸减速度-时间历程曲线。比较计算结果与试验数据表明梁-颗粒模型可有效应用于计算和模拟脆性材料动态破坏问题。 相似文献
5.
提出了一种基于颗粒接触的无网格方法（PCMM）,并编制了相应的C++程序,解决了有限元等网格类方法在模拟边坡失稳滑动过程中的网格畸变问题。该方法利用颗粒离散元中的接触拓扑创建连续介质单元,通过颗粒的运动演化实现连续介质单元的自动删除及重建,通过在连续介质单元中引入考虑应变软化效应的Mohr-Coulomb模型及最大拉应力模型,实现边坡的弹塑性分析及失稳滑移过程的模拟。利用PCMM分析了均质边坡的弹塑性场、沙堆的形成过程及土质边坡的失稳过程。计算结果表明,PCMM在小变形下具有足够的精度,且在模拟材料大变形方面具有明显优势,是一种模拟边坡成灾范围的有效方法。 相似文献
6.
Discrete element method has been widely adopted to simulate processes that are challenging to continuum-based approaches. However, its computational efficiency can be greatly compromised when large number of particles are required to model regions of less interest to researchers. Due to this, the application of DEM to boundary value problems has been limited. This paper introduces a three-dimensional discrete element–finite difference coupling method, in which the discrete–continuum interactions are modeled in local coordinate systems where the force and displacement compatibilities between the coupled subdomains are considered. The method is validated using a model dynamic compaction test on sand. The comparison between the numerical and physical test results shows that the coupling method can effectively simulate the dynamic compaction process. The responses of the DEM model show that dynamic stress propagation (compaction mechanism) and tamper penetration (bearing capacity mechanism) play very different roles in soil deformations. Under impact loading, the soil undergoes a transient weakening process induced by dynamic stress propagation, which makes the soil easier to densify under bearing capacity mechanism. The distribution of tamping energy between the two mechanisms can influence the compaction efficiency, and allocating higher compaction energy to bearing capacity mechanism could improve the efficiency of dynamic compaction. 相似文献
7.
The literature review on discrete element (DEM) model analysis of jigging reveals that an idealized fluid behavior is assumed and the damping of the fluid motion across the mineral bed is generally ignored. A microscopic model based on the principles of Computational Fluid Dynamics (CFD) is used to simulate the liquid flow and stratification of coal particles with a wide size range and density distribution in jigging. Fluid motion is calculated by directly solving the Navier–Stokes equations. Coal particles are moved in a Lagrangian frame through the action of forces imposed by the fluid and gravity. Particle effects on fluid motion are fed back at each time step through calculating the velocity disturbance caused by the particle. Particle–particle and particle–wall collisions are also considered. The snapshots of particle configurations for the simulation of stratification in oscillating flow show that the model predicts the macroscopic behavior, such as segregation and stratification, of particles reasonably well. 相似文献
8.
Simulation of fracturing processes in porous rocks can be divided into two main branches: (i) modeling the rock as a continuum enhanced with special features to account for fractures or (ii) modeling the rock by a discrete (or discontinuous) approach that describes the material directly as a collection of separate blocks or particles, e.g., as in the discrete element method (DEM). In the modified discrete element (MDEM) method, the effective forces between virtual particles are modified so that they reproduce the discretization of a first-order finite element method (FEM) for linear elasticity. This provides an expression of the virtual forces in terms of general Hook’s macro-parameters. Previously, MDEM has been formulated through an analogy with linear elements for FEM. We show the connection between MDEM and the virtual element method (VEM), which is a generalization of FEM to polyhedral grids. Unlike standard FEM, which computes strain-states in a reference space, MDEM and VEM compute stress-states directly in real space. This connection leads us to a new derivation of the MDEM method. Moreover, it enables a direct coupling between (M)DEM and domains modeled by a grid made of polyhedral cells. Thus, this approach makes it possible to combine fine-scale (M)DEM behavior near the fracturing region with linear elasticity on complex reservoir grids in the far-field region without regridding. To demonstrate the simulation of hydraulic fracturing, the coupled (M)DEM-VEM method is implemented using the Matlab Reservoir Simulation Toolbox (MRST) and linked to an industry-standard reservoir simulator. Similar approaches have been presented previously using standard FEM, but due to the similarities in the approaches of VEM and MDEM, our work provides a more uniform approach and extends these previous works to general polyhedral grids for the non-fracturing domain. 相似文献
9.
管涌的发生、发展过程是土骨架相在渗流作用下侵蚀为可动细颗粒相,并随水相在孔隙通道中运移流失的过程。在该过程中,渗流与侵蚀相互耦合,相互促进,水相、土相、可动细颗粒相互作用,因此,管涌过程是一个多场、多相耦合的高度非线性的动态过程。现有的管涌试验结果表明,只有当水力梯度大于起始水力梯度时,细颗粒相才会随水相从土体中运移流失,土体才会发生管涌侵蚀,且管涌稳定后土体的孔隙率（稳定孔隙率）和水力梯度之间存在对应关系,根据该结果,提出管涌稳定孔隙率的概念,修正传统的渗流侵蚀本构方程,建立多孔介质中三相耦合的修正的渗流侵蚀管涌控制方程。最后,针对特定应力状态下的土体建立稳定孔隙率和水力梯度之间的对应关系。基于Galerkin有限元法编制有限元程序,在轴对称情况下对该土体的管涌过程进行数值模拟。结果表明,修正后的管涌控制方程能更全面地描述管涌发生、发展直至稳定状态的特性。 相似文献
10.
This paper presents a numerical scheme for fluid‐particle coupled discrete element method (DEM), which is based on poro‐elasticity. The motion of the particles is resolved by means of DEM. While within the proposition of Darcian regime, the fluid is assumed as a continuum phase on a Eulerian mesh, and the continuity equation on the fluid mesh for a compressible fluid is solved using the FEM. Analytical solutions of traditional soil mechanics examples, such as the isotropic compression and one‐dimensional upward seepage flow, were used to validate the proposed algorithm quantitatively. The numerical results showed very good agreement with the analytical solutions, which show the correctness of this algorithm. Sensitivity studies on the effect of some influential factors of the coupling scheme such as pore fluid bulk modulus, volumetric strain calculation, and fluid mesh size were performed to display the accuracy, efficiency, and robustness of the numerical algorithm. It is revealed that the pore fluid bulk modulus is a critical parameter that can affect the accuracy of the results. Because of the iterative coupling scheme of these algorithms, high value of fluid bulk modulus can result in instability and consequently reduction in the maximum possible time‐step. Furthermore, the increase of the fluid mesh size reduces the accuracy of the calculated pore pressure. This study enhances our current understanding of the capacity of fluid‐particle coupled DEM to simulate the mechanical behavior of saturated granular materials. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
11.
The paper deals with two-dimensional (2D) numerical modelling of hydro-fracking (hydraulic fracturing) in rocks at the meso-scale. A numerical model was developed to characterize the properties of fluid-driven fractures in rocks by combining the discrete element method (DEM) with computational fluid dynamics (CFD). The mechanical behaviour of the rock matrix was simulated with DEM and the behaviour of the fracturing fluid flow in newly developed and pre-existing fractures with CFD. The changes in the void geometry in the rock matrix were taken into account. The initial 2D hydro-fracking simulation tests were carried out for a rock segment under biaxial compression with one injection slot in order to validate the numerical model. The qualitative effect of several parameters on the propagation of a hydraulic fracture was studied: initial porosity of the rock matrix, dynamic viscosity of the fracking fluid, rock strength and pre-existing fracture. The characteristic features of a fractured rock mass due to a high-pressure injection of fluid were realistically modelled by the proposed coupled approach. 相似文献
12.
Water retention of clayey soils with wide particle size distributions involves a combination of capillary and adsorbed layers effects that result into suction–saturation relations spanning over multiple decades of matric suction values. The present study provides a physics-based analysis to reproduce the water retention curve of such soils based solely on particle size distribution and porosity. The distribution of inter-particle pore sizes is inferred through a probabilistic treatment of the particle size distribution, which is then used, together with an assigned pore entry pressure, to estimate the inter-particle water volume at a given suction. The contribution to water content from adsorbed layers is also taken into account by considering the balance of electrochemical forces between water and clay material. The total water content is therefore found by summing up the contribution of inter-particle water, as well as adsorbed layers that form around clay particles and around the individual clay platelets. Comparisons with experimental results on nine different soil samples verify the capability of the model in accurately predicting the wide water retention curves without any prior calibration. Additional to capturing the essential features of the water retention curve with remarkable detail, the analytical model also provides insights into the relative contributions of capillary and adsorbed waters to the overall saturation at different suction regimes. Being based upon easily accessible information such as particle size distribution and void ratio, the model can therefore be considered as a substitute for costly and lengthy laboratory and in situ measurements of water retention curve. 相似文献
13.
Channel morphology and bed sediment erodibility are two crucial factors that significantly affect debris flow entrainment processes. Current debris flow entrainment models mostly hypothesize the erodible beds are infinite with uniform slopes. In this study, a series of small-scale flume experiments were conducted to investigate the effects of bed longitudinal inflexion and sediment porosity on basal entrainment characteristics. Experimental observations revealed that sediment entrainment is negligible at early stages and accelerates rapidly as several erosion points appear. Continual evolution of flow-bed interfaces changes interactions between debris flows and bed sediments, rendering the interfacial shear action involved into a mixed shear and frontal collisional action. Lower bed sediment porosity will change the spatial arrangement and orientation of particle mixture, strengthen the interlocking and anti-slide forces of adjacent sediment particles, and promote the formation of particle clusters, all of which will increase bed sediment resistance to erosion. By examining the post-experimental bed morphology, the slope-cutting amounts and topographic reliefs are determined to positively correlate with longitudinal transition angles. These high topographic reliefs may indicate the propensity of triangular slab erosion, rather than strip-shaped slab erosion, in non-uniform channels with relatively steep erodible beds. Empirical formulas are obtained that denote the relationships among bed sediment strength, channel curvature radius, and sediment porosity through a multi-parameter regression analysis. This study may aid in clarifying the complex coupling effects of spatial variations in debris flow dynamics as well as sediment erodibility and bed morphology in non-uniform channels with abundant seismic loose material. 相似文献
15.
把块裂结构岩质边坡看作拟连续介质与块裂介质广义双重介质,运用现代分形理论,建立了1-2级主干裂隙分形网络,并根据两类介质接触处水头相等及节点流量相等建立了合理的块裂结构边坡三维水力学模型,给出了有限元解法,模拟了三峡工程的蓄水位至175 m时边坡渗流场与应力场的变化,在应用上有较好的可靠性与广泛性。 相似文献
16.
页岩作为典型的非常规储层,基质孔隙小,渗透率极低,水平井多级水力压裂为其商业开发的主要手段。准确模拟页岩气产能,应同时考虑水力裂隙和天然裂隙的渗流。基于离散裂隙模型和等效连续模型建立页岩气渗流数学模型,利用有限元分析方法进行数值求解,研究不同走向裂隙组对页岩气井产能的影响。研究认为,页岩基质为气体的生产提供了主要气源,天然裂隙作为渗流的主要通道,将气体输送到水力裂缝,进而到达井筒。模拟结果表征,离散裂隙的渗流特征对于页岩气井的产能有重要影响。根据页岩储层的天然裂隙走向,可以优化相应的水平井方位。对于二维离散裂隙网络模型,水平井沿着2个裂隙组夹角的平分线更有利于生产。 相似文献
17.
采用颗粒离散单元法进行动力计算时,人工截断边界上需设置吸收边界条件,以防止波的反射。鉴于颗粒离散单元数值计算模型的人工边界上颗粒单元半径大小不一、边界面凸凹不平,在连续介质的黏性、黏弹性、自由场边界条件方程基础之上,推导出适用于离散介质的等效方程。在离散介质的黏性边界条件等效方程中引入微调系数,提出比值迭代法以快速确定其最优值,以实现对波的最佳吸收。采用二维颗粒离散单元计算软件PFC2D,分别建立黏性、黏弹性、自由场边界条件相关数值分析模型,探讨颗粒分布模式对黏性边界上颗粒单元半径、速度分布及比值迭代过程的影响;采用外源波动算例及经典Lamb问题算例验证黏弹性边界设置方法的正确性;通过隧洞算例检验提出的自由场边界条件设置方法的正确性。 相似文献
18.
In the present work, a methodology for setting up virgin stress conditions in discrete element models is proposed. The developed algorithm is applicable to discrete or coupled discrete/continuum modeling of underground excavation employing the discrete element method (DEM). Since the DEM works with contact forces rather than stresses there is a need for the conversion of pre-excavation stresses to contact forces for the DEM model. Different possibilities of setting up virgin stress conditions in the DEM model are reviewed and critically assessed. Finally, a new method to obtain a discrete element model with contact forces equivalent to given macroscopic virgin stresses is proposed. The test examples presented show that good results may be obtained regardless of the shape of the DEM domain. 相似文献
19.
区域性渗流场由于范围大,通常忽略掉竖向的水头损失而用水平面二维有限元进行分析,但在城区区域渗流场中悬挂式阻水结构物量多,竖向绕流现象突出,用一般水平面二维有限元进行分析将产生较大的偏差。针对这一问题在悬挂式结构体周边引入一层无厚度界面接触单元。在水平面有限元中利用它考虑因竖向绕流产生的局部水头损失。利用阻力系数的方法和二维状态下截面突变流道流动的解析解,同时考虑防渗墙本身的弱透水性,推导了接触单元的导水矩阵。应用接触单元编制了相关程序,在水平面有限元中利用接触单元对悬挂式阻水结构进行分析,并将结果与三维有限元进行比较,两者的计算结果在流量上保持一致,在水头上前者是后者在同一位置处沿高程方向平均值。提出的该计算方法,对于用水平面二维有限元分析含有悬挂式阻水结构的大区域渗流问题,是一个有效的手段。 相似文献
20.
Internal instability is a phenomenon of fine particle redistribution in granular materials under the seepage action and consequent change in the soil’s internal structure and hydraulic and mechanical properties. It is one of the primary causes of failures of sand-gravel foundations and embankment dams. The criteria establishment is considered the key to solving the erosion problems, so the existing internal stability criteria need a review and classification to study the recent development trends in soil seepage and erosion. Therefore, this paper aims at reviewing the internal stability factors of gap-graded soil with a focus on the internal erosion mechanism and internal stability evaluation based on geometric and hydraulic criteria. Firstly, the paper compared the effect of several commonly used geometric criteria for gap-graded soil evaluation, such as particle size, fine content, void ratio, and fractal dimension. Furthermore, it provided a hydraulic criteria overview and analyzed the effects of the hydraulic gradient, hydraulic shear stress, confining pressure, and pore velocity on internal erosion. The geometric–hydraulic coupling methods were introduced, with a detailed elaboration of the erosion resistance index method based on accumulated dissipated energy. The capabilities and limitations of these criteria were discussed throughout the paper. It was found that combined Kezdi’s criterion and Kenney and Lau’s criterion is more reliable to evaluate internal stability of soil. The gap-graded soil with fine particle content higher than 35% is not necessarily internally stable. Finally, the energy-based method (erosion resistance index method) can effectively reproduce the total amount of erosion mass and the final spatial distribution of fine particles and identifies erosion. The review's outcome can be used as a basis to evaluate the internal erosion risk for gap-graded soils. The evaluation methods discussed here can help identify the zones of relatively high erosion potential. 相似文献
|
