三相耦合渗流侵蚀管涌机制研究及有限元模拟

管涌的发生、发展过程是土骨架相在渗流作用下侵蚀为可动细颗粒相,并随水相在孔隙通道中运移流失的过程。在该过程中,渗流与侵蚀相互耦合,相互促进,水相、土相、可动细颗粒相互作用,因此,管涌过程是一个多场、多相耦合的高度非线性的动态过程。现有的管涌试验结果表明,只有当水力梯度大于起始水力梯度时,细颗粒相才会随水相从土体中运移流失,土体才会发生管涌侵蚀,且管涌稳定后土体的孔隙率（稳定孔隙率）和水力梯度之间存在对应关系,根据该结果,提出管涌稳定孔隙率的概念,修正传统的渗流侵蚀本构方程,建立多孔介质中三相耦合的修正的渗流侵蚀管涌控制方程。最后,针对特定应力状态下的土体建立稳定孔隙率和水力梯度之间的对应关系。基于Galerkin有限元法编制有限元程序,在轴对称情况下对该土体的管涌过程进行数值模拟。结果表明,修正后的管涌控制方程能更全面地描述管涌发生、发展直至稳定状态的特性。     

均质土坝管涌发展过程的渗流场空间性状研究

堤坝发生管涌破坏时,渗流场具有强烈的空间特性,管涌发展过程的数值模拟对于堤坝工程管涌预防与治理都具有重要的参考价值。提出以一维通道嵌入三维块体的方法,克服了在三维问题有限元模拟方法上存在的困难。在探讨均质土坝管涌发展过程模拟方法的基础上,使用数值方法分析了土坝管涌发展过程的渗流场时空分布特性;揭示了管涌发展过程中渗流力集中状况及其变化规律。研究表明,改善渗流场的形态是防治堤坝工程管涌的关键措施。     

细砂层埋深对堤基管涌影响的试验研究

不同土层结构的堤基,其管涌发生和发展的情形不同。以前的研究多集中于双层堤基的情况。事实上,含夹砂层堤基中,位于砂砾石层内的细砂层会对堤基管涌的发生产生重大影响,并且这些影响随细砂层在砂砾石层中的位置不同而不尽相同。利用室内砂槽试验,模拟当细砂层位于砂砾石层内部不同深度时含夹砂层堤基管涌的发展过程,研究细砂层埋深对 堤基管涌的影响。试验结果表明,细砂层埋深存在一个临界深度,当细砂层埋深较浅时,堤基的临界水力梯度较小,易于发生管涌,且管涌发生以后,堤基被侵蚀的速率较大,出砂较多;堤基破坏的临界水力梯度随着细砂层埋深的增加而增大,当细砂层埋深大于临界埋深时,堤基的临界水力梯度基本不变,同时堤基抵抗管涌破坏的能力较强,其破坏形式与双层堤基 相类似。     

黄土缓坡片蚀过程及其水力参数适宜性试验研究

为确定黄土缓坡片蚀过程中侵蚀限制条件的变化及水力参数的适宜性,采用模拟降雨试验,在60 mm/h雨强、3种坡度（7.5°、10°、15°）和2种坡长（5 m和10 m）条件下,验证含沙量作为表征参量判定片蚀过程侵蚀限制条件的可行性,提出检验水力参数与产沙量作用关系适宜性的评价方法。结果表明:①依据含沙量随径流的变化规律,片蚀过程依次经历输沙能力限制阶段、剥蚀能力限制阶段Ⅰ和剥蚀能力限制阶段Ⅱ共3个侵蚀阶段,其中剥蚀能力限制阶段Ⅰ的侵蚀量占绝对优势。②从水力参数的适宜性来看,径流功率和雷诺数均能很好地表征与产沙量的关系,径流剪切力和弗劳德数则不能有效反映与产沙量的关系。③在不同侵蚀限制阶段,径流功率和雷诺数与产沙量的作用关系是不同的。输沙能力限制阶段和剥蚀能力限制阶段Ⅱ时,径流功率和雷诺数与产沙量分别呈指数关系和线性正相关关系;剥蚀能力限制阶段Ⅰ时,该2种水力参数均与产沙量呈线性负相关关系,但作用关系不唯一,且存在雷诺数临界值,即在输沙能力限制阶段和剥蚀能力限制阶段Ⅱ,径流功率和雷诺数与产沙量的作用关系不受坡长效应影响;在剥蚀能力限制阶段Ⅰ,需要同时考虑坡长和雷诺数临界值的共同影响。该研究结果可为深入理解片蚀过程的复杂性提供参考依据。     

Quantitative analysis of piping erosion micro-mechanisms with coupled CFD and DEM method

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.     

An auto-adaptive moving mesh method for the numerical simulation of piping erosion

A numerical approach for simulating piping erosion is proposed based on an auto-adaptive moving mesh. This approach adapts a flownet mesh (FM) to discretize the domain. It reduces the simulation of two-dimensional transport of moving particles to one-dimensional problem by taking advantages of flow tube concept. Furthermore, it employs a stepwise procedure to decouple the solid-liquid interaction equations for modeling piping erosion. Specifically, the FM is updated and evolved automatically in accordance to simulated erosion since the flownet is refined in high-velocity regions, where the piping erosion progresses. Therefore, the FM is auto-adaptive to the piping erosion paths.     

基坑防渗体局部失效渗流场变化特征

采用三维有限元方法分析了基坑开挖工程在防渗体出现局部失效情况下，引起渗流在失效部位的集中及渗流场空间状况。提出了以一维通道嵌入三维块体的方法模拟管涌发展过程及渗流场变化特征。讨论了管涌通道渗透性对基坑管涌发展规律的影响。     

A finite element method for localized erosion in porous media with applications to backward piping in levees

This paper presents a computational method able to effectively model both the simultaneous processes typically observed in backward erosion piping, ie, the pipe tip propagation and the conduit cross-section enlargement. The numerical method is based on the novel formulation of a problem of localized erosion along a line propagating in a multidimensional porous medium. In this line, a conduit with evolving transverse size is embedded, which conveys a multiphase flow. The two systems, porous medium and pipe, are bridged by exchange terms of multiphase fluid mass and by a shared fluid pressure field. On the contrary, different fields are considered to describe flows, which are assumed as Darcian in the porous medium and turbulent in the conduit. These two flows drive pipe propagation and enlargement, respectively, as modeled by means of proper erosion kinetic laws. The corresponding numerical formulation is based on the combination between one- and multidimensional finite elements, to model the erosion conduit and the porous medium, respectively. Several simulations are proposed to demonstrate the ability of the proposed approach in reproducing available experimental data of real-scale tests on levees. Our results point out the crucial role played by the combined influence of pipe propagation and enlargement, as well as of three-dimensional (3D) effects. We also assess the mesh independence of the proposed numerical solution, particularly as concerns the calculated pipe propagation history.     

三层堤基中细砂层厚度对管涌影响的试验研究

对由弱透水黏土层、细砂层和强透水砂砾层组成的三层堤基进行了管涌发展的砂槽模型试验,为了便于观察分析,细砂层由各种颜色的细彩砂依次排列在砂砾石层上表面,通过改变彩砂层的厚度分析研究了不同细砂层厚度对管涌发生、发展机制及过程的影响。试验结果表明,三层堤基细砂层厚度的不同使管涌发生的临界水力梯度、涌砂量和通道发展的速度不同,与双层堤基有很大区别。临界水力梯度是由多种元素决定的,包括破坏土体的性质及其整体性等;细砂层的存在使流量在渗透变形初期对涌砂不敏感;在试验中发生的相同水位下多次间歇性涌砂,其原因一方面是颗粒在运动过程中发生堵塞,另一方面是通道边界的土体失去支撑发生应力释放,抵抗力随着时间逐渐减小。     

缺级粗粒土管涌类型的判别方法

利用大型渗透变形仪,对不同最大粒径和细粒含量的缺级粗粒土的渗透破坏特性进行了系统的试验研究。研究结果表明,缺级粗粒土在细料含量小于30 %的情况下都发生管涌型渗透破坏,基本符合前人提出的判别流土和管涌的准则。管涌破坏的类型可划分为发展性管涌和非发展性管涌两种,采用刘杰和康德拉且夫对发展性和非发展性管涌的判别方法都不能区分本次试验的管涌破坏类型。在康德拉且夫提出的判别粗粒土管涌类型的方法的基础上,考虑细料对粗料孔隙的填充程度及试样中缺级颗粒百分含量的影响,建立了缺级粗粒土管涌类型的判别式,为缺级粗粒土管涌类型的判别提供了新的思路和方法。     

Numerical simulation of the developing course of piping

Piping is the main cause of dangerous situations and one of the most serious threats to levees in China. It is unlikely that Darcy’s law is suitable for evaluation of piping zones; therefore, the ordinary seepage computation based method of simply increasing the permeability coefficient of the piping zone is not appropriate. The discrete element method (DEM) may be a promising alternative to the use of Darcy’s law, but the large number of particles hampers its effective application in practice. In this study, a new approach based on combination of pipe flow in the erosion channel and seepage analysis in the remainder of the levee is proposed. Based on the changeable boundary between the two types of computational domains, the element free Galerkin method (EFG) was employed to facilitate the efficiency of coupling iteration. Examples show that this approach can simulate the main characteristics of developing processes in the erosion channel, which implies that this simplified method can be employed practically.     

Modeling of fluid–solid interaction in granular media with coupled lattice Boltzmann/discrete element methods: application to piping erosion

In this article, we present a numerical method to deal with fluid–solid interactions and simulate particle–fluid systems as encountered in soils. This method is based on a coupling between two methods, now widely used in mechanics of granular media and fluid dynamics respectively: the discrete element (DE) method and the lattice Boltzmann (LB) method. The DE method is employed to model interactions between particles, whereas the LB method is used to describe an interstitial Newtonian fluid flow. The coupling presented here is a full one in the sense that particle motions act on fluid flow and reciprocally. This article presents in details each of the two methods and the principle of the coupling scheme. Determination of hydrodynamic forces and torques is also detailed, and the treatment of boundaries is explained. The coupled method is finally illustrated on a simple example of piping erosion, which puts in evidence that the combined LB–DE scheme constitutes a promising tool to study coupled problems in geomechanics. Copyright © 2011 John Wiley & Sons, Ltd.     

硅酸钾材料加固土遗址耐风蚀颗粒元模拟分析

西北干旱地区土遗址受风化、风蚀等破坏严重,大量土质文物亟待加固抢修。加固后土遗址的各耐环境因素及加固机制研究是土遗址加固的理论基础。首次引入颗粒元程序PFC,通过改变模型中颗粒间平行连接强度,对硅酸钾（简称PS）加固前后的土样进行数值模拟。在考虑实际土样颗粒粒径和密度的前提下,拟合了生土PS加固前后的抗压和抗拉强度,并将拟合后的颗粒元模型应用于风蚀模拟。通过随机生成挟沙风颗粒,以一定的速度撞向土体,模拟挟沙风的吹蚀作用。挟沙风颗粒数与循环步数成正比例,因此,可以用挟沙风颗粒数来代表吹蚀时间的长短。挟沙风颗粒的速度则代表挟沙风风速。模拟结果表明,在20 m/s的挟沙风吹蚀作用下,风蚀程度随吹蚀时间的增加而增大,未加固土样的风蚀程度增幅度远大于加固土样;同样吹蚀时间条件下,加固土样的抗风蚀强度明显高于未加固土样。这些模拟结论与风洞试验结果的统计规律一致。本研究拟合的颗粒流模型可进一步应用于PS加固机制研究及耐风蚀、雨蚀、冻融等诸环境影响分析研究。     

Analysis of development and depth of backward erosion pipes in the presence of a coarse sand barrier

Backward erosion piping (BEP) is a failure mechanism that can affect the safety of water-retaining structures. It can occur when a local anomaly on the downstream side of an embankment causes a concentration of seepage flow at that location. Shallow pipes may then form, progressing in the upstream direction and leading to a collapse of the water-retaining structure. A novel and economically appealing measure against BEP is the coarse sand barrier (CSB), which is now being developed in a multiscale experimental programme in the Netherlands. The method involves placing a trench filled with coarse sand below the blanket layer on the downstream side of the embankment. The CSB prevents the upstream progression of the pipe and significantly enhances resistance to BEP. This paper presents medium-scale laboratory tests involving a range of sands, barrier depths and relative densities. The piping process and the observations of pipe progression in the presence of a CSB are presented, followed by a conceptual model. The presence of a CSB changed the erosion pattern. It resulted in pipe formation perpendicular to the flow direction over the entire width of the barrier before the barrier was damaged. The findings also demonstrate the effect of material properties on pipe initiation, progression and pipe depth. Measurements of the pipe depth are presented and analysed, revealing the significance of pipe depth for understanding the piping process. This analysis shows considerable erosion in the downstream background sand and demonstrates that erosion profiles and measured pipe depths are significantly larger than in BEP tests without a CSB.

     

沙漠公路风沙土路基风蚀破坏试验研究

以沙漠公路风沙土路基为研究对象,通过室内风蚀风洞试验研究路基的风蚀破坏规律,以及路基不同断面对风沙流运动的影响。以路基高度、路基边坡坡率和路基宽度作为路基断面主要设计参数,研究不同路基断面下风沙流扰动、增速、减速、恢复的过程,以及路基周围风速流场的变化特征,分析路基病害较未病害时路基周围流场的变化。试验结果表明：路基高度和边坡坡率对风沙流运动的影响较大。随路基高度增加,路基对风沙流流场扰动增强,迎风坡坡顶处吹蚀破坏和背风坡坡底处堆蚀破坏越显著,在确定的路基边坡坡率下,路基模型高度为250 mm较模型高度为60 mm时,迎风坡坡顶风速增加1.13倍,背风坡坡底风速减小2.53倍,建议沙漠公路路基高度宜小于2.5 m。进一步,在确定的路基高度下,比较不同的边坡坡率对路基沿程风速的影响,发现当路基边坡坡率为1:1.75时,路基沿程风速变化不明显,沙漠公路风沙土路基不宜被风蚀破坏。     

Volumetric sand production model and experiment

A sand production model was developed for volumetric sand production predictions that take into account the effects of the external stresses and fluid flow rate. The model couples the poro‐mechanical behaviour of the solid–fluid system with the erosion behaviour of the solids due to fluid flow. It predicts reasonably experimental volumetric sand production data from a hollow cylinder test on a weak sandstone. The test results show that in weak and compactive sandstones, sand production is associated with decohesioning and plasticification of a zone around the inner hole which can then be mobilized by the hydrodynamic forces of the fluid flow. The sand production rate increases both with external applied stress and fluid flow rate but it is constant with time under constant external stress and fluid flow rate. In both cases a critical lower limit has to be exceeded for sand production initiation. Copyright © 2001 John Wiley & Sons, Ltd.     

土质陡坡降雨侵蚀的数学模型及求解方法

建立了土质陡坡降雨侵蚀的数学模型,模型包括陡坡细沟间坡面流控制方程、变沟宽陡坡细沟流控制方程、陡坡细沟间坡面侵蚀泥沙连续方程及陡坡变沟宽细沟侵蚀泥沙连续方程。通过紊流冲击分布概率确定土壤剥蚀方程,确定了泥沙源/汇项的表达形式及细沟高程与宽度的表达式。通过有限元和有限差分相结合的方法求解该数学模型,利用有限单元法对数学模型的水动力方程、运动方程及泥沙连续方程进行空间上的离散,利用有限差分法解决时间域的问题,得出了顺序求解的数值计算公式及模型数值求解的具体方法和步骤。     

堤坝管涌的室内试验与颗粒流细观模拟研究

本文对水流垂直方向的堤坝管涌现象分别进行了室内模型模拟和颗粒流模拟实验。在室内模型试验模拟中采用先进的数码成像跟踪设备在细观的尺度上对土样在管涌发生发展的全过程进行了观察和记录,为颗粒流细观模拟提供了必要的细观参数及宏观依据。提出的二维孔隙率与三维孔隙率转换的关系公式很好的将室内模型试验的孔隙率同颗粒流数值模拟结合,运用颗粒流程序采用固液耦合的分析模型对管涌现象作进一步的数值模拟分析。根据模拟结果,分析了管涌产生、发展在颗粒细观尺度的原因,对管涌产生发展中的试样细观变化规律及流速的变化规律进行了分析。这些结果为PFC2D在渗流和渗透破坏等方面的深入研究应用提供了一定的实用依据,同时也为在细观尺度上解释和分析管涌现象提供了一条新的途径。     

青海省格拉沟形成演化机制及其主要工程地质问题

格拉沟位于青海省内陆滩地，其形成演化反映了西北地区典型的水土流失过程。文章在收集大量野外资料的基础上，研究了格拉沟的形成演化机制，并对该沟道主要工程地质问题进行了探讨。分析认为，格拉沟形成演化是该地区特殊的内因外因共同作用的结果，基于此对其发育演化阶段进行了划分与预测。该沟主要工程地质问题为边坡稳定性、渗透稳定性、河流侵蚀等。其中边坡稳定性主要包括天然岸坡整体稳定性及崩塌、坠溜、滑塌等；渗透稳定性主要包括潜蚀与管涌；河流下切与溯源侵蚀则成为岸坡失稳破坏的重要诱发因素。探讨格拉沟形成演化机制及其主要工程地质问题，对于研究西北地区水土流失机理及其产生的环境问题具有一定的借鉴作用。     

银西高铁董志塬地区边坡侵蚀特性分析

黄土坡面侵蚀特性研究对于铁路边坡及路基防护具有重要的意义。通过银西高铁董志塬段某路基护坡坡面冲刷试验,获得了不同冲刷历时、冲刷流量、坡度等条件下的坡面冲刷结果,并对坡面流水动力学特性、坡面产沙规律、坡面产沙机理进行了分析,得出以下结论:（1）坡顶和坡底比坡面中部更易受侵蚀;30°～60°斜坡在较小的冲刷强度下也能产生较明显的侵蚀沟,宜采取45°左右的多级矮坡来减弱侵蚀强度。（2）坡面流水动力学特性分析表明,试验工况水流主要处于过渡流区;平均流速与冲刷流量、坡度呈幂函数关系;达西阻力系数与冲刷流量、坡度呈负相关,且与雷诺数相关性较低。（3）平均含沙量随冲刷流量与坡度的增大而增大,随历时近似线性增加,约20 min以后,含沙量基本稳定,此过程为坡面沟道发展阶段。（4）坡面侵蚀产沙量与侵蚀切应力、有效水流功率都呈正相关,与前者近似呈线性增大关系,而与后者近似呈幂函数关系。