Erosion of Schist due to Solid Particle Impingement

The erosion of schist due to impinging quartz and glass particles is investigated. The experimental set-up used to perform the erosion tests is briefly described. The erosion process is characterized by two failure modes, namely intercrystalline fracture and lateral cracking. The power exponent for particle velocity effects is 2.44. The schist shows the features of plastic deformation, and a model proposed by Ruff and Wiederhorn (Erosion. Academic Press, New York, pp 69–126, 1979) can be applied to the erosion process. The formation of glassy microfibers is noted, and it is attributed to localized heating due to the particle impact process.     

流态地聚物固化土强度特性及其强度预测

地聚物胶凝材料能够替代水泥基胶凝材料作为固化剂应用于狭窄肥槽回填等工程问题中,有效降低水泥生产过程中的污染及能耗,但目前对于流态地聚物固化土胶凝材料的研究较少。采用3种新型绿色胶凝材料联合碱激发剂固化工程渣土形成流态地聚物固化土,通过对比其无侧限抗压强度,探究每种胶凝材料对于固化土强度特性的影响,同时建立强度预测模型,分析不同因素对于强度的影响程度。研究结果表明：固化土的强度随着碱激发剂模数的增加先提高后降低;固化土强度随着高炉矿渣（GGBS）、粉煤灰、稻壳灰掺量的增加均呈上升趋势,随着稻壳灰粒径的增长呈下降趋势;碱激发剂模数增至1.2、GGBS掺量增至10%、粉煤灰掺量增至8%和稻壳灰掺量增至11%时,固化土强度提升最为显著;强度预测模型预测结果的平均相对误差仅为5.57%,预测结果较为精准;预测模型中各层权值的计算结果表明养护龄期对于固化土强度影响最大,稻壳灰粒径影响程度最小。研究结果可以为固化土在实际工程的应用提供理论支持。     

输沙能力模式在坡面非恒定输沙过程中的应用

在坡面土壤侵蚀输沙计算中,可选用的输沙能力模式较多,但这些模式多运用于恒定输沙,在非恒定输沙中运用不多。建立了坡面土壤侵蚀非恒定输沙数学模型,模型中的输沙能力运用水流切应力、水流功率、单位水流功率3种模式,对模型过程采用有限差分格式离散求解。根据实测水沙资料进行模型参数率定,运用3种输沙能力模式于不同降雨强度、不同坡度的非恒定坡面输沙过程中。结果表明:在坡面非恒定输沙计算中,不同输沙能力模式对计算结果有明显影响,在降雨强度较小时,单位水流功率模式结果较其他两个模式为好,而在雨强较大时,切应力模式计算结果较好。     

砂岩类脆性无序介质连续破坏过程的细观模拟

为了分析砂岩类无序介质在各种荷载作用下的连续破坏过程,本文提出了梁-颗粒细观模型。在此模型中,介质被离散为一系列的颗粒单元,这些颗粒单元均由弹脆性的梁单元联结。颗粒的力学行为由离散单元法和有限单元法确定。当某个梁单元所受应力超过其强度时,就随即将它从计算网络中剔除,以模拟介质的破坏行为。为了验证数值模型的可靠性,分别从破坏模式和荷载-位移曲线两方面与物理模型实验进行了对比,结果表明两者基本一致。此外,探讨了无序性对介质裂纹分布形式的影响,并用分维表述了裂纹分布与荷载强度之间的内在联系。      

3D finite element modeling of sand particle fracture based on in situ X‐Ray synchrotron imaging

Compressive loading of granular materials causes inter‐particle forces to develop and evolve into force chains that propagate through the granular body. At high‐applied compressive stresses, inter‐particle forces will be large enough to cause particle fracture, affecting the constitutive behavior of granular materials. The first step to modeling particle fracture within force chains in granular mass is to understand and model the fracture of a single particle using actual three‐dimensional (3D) particle shape. In this paper, the fracture mode of individual silica sand particles was captured using 3D x‐ray radiography and Synchrotron Micro‐computed Tomography (SMT) during in situ compression experiments. The SMT images were used to reconstruct particle surfaces through image processing techniques. Particle surface was then imported into Abaqus finite element (FE) software where the experimental loading setup was modeled using the extended finite element method (XFEM) where particle fracture was compared to experimental fracture mode viewed in radiograph images that were acquired during experimental loading. Load‐displacement relationships of the FE analysis were also compared with experimental measurements. 3D FE modeling of particle fracture offers an excellent tool to map stress distribution and monitors crack initiation and propagation within individual sand particles. Copyright © 2015 John Wiley & Sons, Ltd.     

Pore‐scale modeling of surface erosion in a particle bed

A fully coupled transient two‐dimensional model was employed to study fundamentals of flood‐induced surface erosion in a particle bed. The interaction of the liquid and solid phases is the key mechanism related to surface erosion. The solid phase was idealized at a particle scale by using the discrete element method. The fluid phase was modeled at a mesoscale level and solved using the lattice Boltzmann method. The fluid forces applied on the particles were calculated on the basis of the momentum the fluid exchanges with the particle. The proposed approach was used to model both single particles and particle beds subjected to Couette flow conditions. The behavior of both the single particle and the particle bed depended on particle diameter and surface shear fluid velocity. The conducted simulations show that the fluid flow profile penetrates the bed for a small distance. This penetration initiates sheet‐flow and surface erosion as the fluid interacts with particles. The effect of suppressing particle rotation on the fluid‐induced forces on the particle was also examined. Suppressing particle spinning may lead to underestimated erosion rate. Results of fluid and particle velocities were compared against experimental results and appeared to agree with the observed trends.Copyright © 2013 John Wiley & Sons, Ltd.     

岩石化学风化时效效应的离散元模拟

微观上,岩石由晶粒和晶粒之间的联接组成。在岩石的风化过程中,晶粒之间的联接会被削弱,导致岩石微观缺陷的发展,从而表现为显著的时效效应。在已有胶结模型的基础上,将岩石的化学风化过程等效为颗粒之间胶结物质不断溶蚀进而影响颗粒之间胶结强度和刚度的过程,基于此建立了考虑化学风化的微观接触模型。随后将模型植入离散元软件,运用离散元方法（DEM）,模拟和对比分析了室内腐蚀试验与DEM模拟化学风化试验,结果表明：该模型可以反映岩石强度的化学风化时效效应,即随着反应时间增加,强度减小;可以反映单轴压缩试验中裂隙压密阶段以后变形的时效效应,即随着反应时间的增加,弹性模量、破坏时应变减小;粒间胶结的破坏,即微观缺陷的发育,是导致岩石化学风化时效效应的微观原因。同时DEM模拟风化与加载耦合的试验现象表明,施工过程中的施工速度会影响岩石的强度,施工应该选择合适的速率。     

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

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

Stress‐dependent elastic wave velocity of microfractured sandstone

A model for the stress‐dependent elastic wave velocity response of fractured rock mass is proposed based on experimental evidence of stress‐dependent fracture normal and shear stiffness. Previously proposed models and previous experimental studies on stress‐dependent fracture stiffness have been reviewed to provide a basis for the new model. Most of the existing stress‐dependent elastic wave velocity models are empirical, with model parameters that do not have clear physical meanings. To propose the new model, the rock mass is assumed to have randomly oriented microscopic fractures. In addition, the characteristic length of microfractures is assumed to be sufficiently short compared to the rock mass dimensions. The macroscopic stress‐dependent elastic wave velocity response is assumed to be attributed to the stress dependency of fracture stiffness. The stress‐dependent fracture normal stiffness is defined as a generalized power law function of effective normal stress, which is a modification of the Goodman's model. On the other hand, the stress dependency of fracture shear stiffness is modeled as a linear function of normal stress based on experimental data. Ultrasonic wave velocity responses of a dry core sample of Berea sandstone were tested at effective stresses ranging from 2 to 55 MPa. Visual observation of thin sections obtained from the Berea sandstone confirms that the assumptions made for microstructure of rock mass model are appropriate. It is shown that the model can describe the stress‐dependent ultrasonic wave velocity responses of dry Berea sandstone with a set of reasonable material parameter values. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

The scaling law in the hole erosion test with a constant pressure drop

A process called ‘piping’, which often occurs in the soil at dams, levees, and dykes, involves the formation and development of a continuous tunnel between upstream and downstream ends. The hole erosion test is commonly used to quantify the critical stress and the rate of piping erosion progression. The aim of this study is to draw up a model for interpreting the results of this test. A characteristic internal erosion time is defined and expressed as a function of the initial hydraulic gradient and the coefficient of surface erosion. It is established here that the product of the coefficient of erosion and the flow velocity is a significant dimensionless number: when this number is small, the kinetics of erosion are low, and the particle concentration does not have any effect on the flow. This finding applies to most of the available test results. Theoretical and experimental evidence is presented showing that the evolution of the pipe radius during erosion with a constant pressure drop obeys a scaling exponential law. Copyright © 2008 John Wiley & Sons, Ltd.     

On the heat transfer coefficient between rock fracture walls and flowing fluid

Two analytical solutions are derived to model the heated flow-through experiments for granite fractures in the literature. The first model, which assumes an identical/continuous temperature between the bulk fluid and fracture surfaces, represents an upper bound solution of water temperature in rock fractures. The second model including the empirical parameter of heat transfer coefficient is used to calculate the average heat transfer coefficient based on the available experimental data. The obtained heat transfer coefficients are smaller than that from the thermal boundary layer theory for flat plates, but larger than the previous estimates. A power function is fitted to describe the relation between heat transfer coefficient and flow velocity. Both models show that water temperature increases non-linearly along fracture plane.     

满流条件下管道破损诱发渗流侵蚀的试验研究

针对富水砂层中地下管道破损诱发地面沉降的问题,设计了一套可视化试验装置,对骨架粒径d90 =1.45～8.45 mm的11种砂土样、5种满流流速下土体渗流侵蚀诱发地表沉降的规律进行了研究。研究表明：（1）管道破损诱发的土体渗流侵蚀有3种模式,分别为：只突水无沉降、形成土拱并发生沉降和溃砂沉降;（2）土体骨架粒径、破损口尺寸、厚跨比是决定土体产生何种渗流侵蚀模式的主要因素,满流流速主要影响土体沉降范围;（3）管道破损口上方土体形成土拱或溃砂时,土体骨架粒径d90与厚跨比r的关系为：当8.0≥r≥4.2时,d90随r的增大呈抛物线下降;当12.5≥r≥8.0时,d90保持不变;（4）管道破损口上方土体形成土拱或溃砂时,流速为0的初始沉降半径、沉降深度由破损口直径D与土体平均粒径d50的比值（D/d50）决定;沉降半径、沉降深度均随满流流速的增大而线性增大;形成土拱时,沉降半径、沉降深度随流速增大的扩展速度（VL、VH）与D/d50呈对数关系;溃砂时,沉降半径随流速增大的扩展速度VL规律如下：当23.0≥D/d50≥6.0时,VL随D/d50的增加而线性增大,当42.0≥D/d50≥23.0时,VL随D/d50的增大呈对数减小。     

电渗排水固结中电极材料的对比试验

采用辽宁沿海地区淤泥质软土对铁、铜、铝和新型复合电极等不同电极材料进行室内电渗排水固结试验,从有效电势、电流、排水量与排水速率、能耗、电渗后土体含水率和承载力等对电渗效果进行分析。新型复合电极材料是由碳纤维、塑料排水板与铁片组成,其充分利用各种材料的优点,既解决了铁与碳纤维直接相连通电后易断的问题,又解决了电极腐蚀严重的问题。多组试验结果表明：新型复合电极材料的电渗排水固结使土体的加固强度提升较高,且相对均匀,加固强度达到160 kPa以上的土体面积占加固总面积的75%,说明新型复合电极材料的优越性。金属电极的腐蚀较为严重,且主要发生在阳极。根据电渗后阳极的腐蚀比,新型复合电极的腐蚀量明显低于金属电极,且新型复合电极可提高加固效果,节约电极材料成本。     

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.     

多尺度内聚颗粒模型破碎过程研究

岩石破碎过程涉及到多个变量（应力、应变和孔隙率等）变化和裂纹的产生、拓展和聚集,是研究岩石破碎机制的重要途径。为了克服现有研究中岩石颗粒模型未考虑岩石内部特征问题,针对岩石内部颗粒非均匀分布、聚集的特点,开展了岩石轴压破碎试验和岩石岩相分析试验,并在此基础上,构建符合真实岩石内部特征的多尺度内聚颗粒模型。根据离散元的颗粒黏结模型（BPM）理论,求解了多尺度内聚颗粒模型不同粒级颗粒间黏结键的力学关系,发现与二级颗粒形成的黏结键断裂判据为 ≥ 2 GPa,三级颗粒之间形成的黏结键断裂判据为 ≥ 6 GPa,并基于该判据建立了用于模拟颗粒模型破碎的演化模型。通过模拟轴压破碎试验,破碎演化模型可以从细观角度得到颗粒模型各颗粒间黏结键承受力的实时变化和岩石破碎过程中黏结键从外而内的断裂顺序;从宏观角度得到岩石内部裂纹呈V形从上表面两端延伸并相交于岩石中部。通过与岩石轴压破碎试验结果对比发现,模拟试验得到的岩石裂纹特征与岩石轴压破碎试验结果相似,验证模型的可靠性,实现了从细观和宏观角度分析岩石破碎过程。     

The settling of consecutive spheres in viscoplastic fluids

One of the factors contributing to the uncertainties involved in the estimation of particle settling velocity in viscoplastic fluids is the time-dependent effect where the viscous parameters of the fluid change as a particle flows through and shears the medium. These changes, particularly at low shear Reynolds numbers, are reflected in the settling velocity of a following sphere that is released some time after an initial one, with the following sphere having a significantly greater velocity. This study found that changes in both fall velocity and equivalent viscosity can be correlated satisfactorily by a power law equation to the dimensionless form of the time interval between releases, and the rheogram shape factor for the fluid. A collision of particles occurs in cases where the time interval between releases is small, after which the particles combine and travel at a terminal velocity. A new variable, β, which takes into account the different surficial stress of the combined spheres, was introduced to the correlation of Wilson et al. [Wilson, K.C., Horsley, R.R., Kealy, T., Reizes, J.A., Horsley, M.R., 2003. Direct prediction of fall velocities in non-Newtonian materials. Int. J. Miner. Process. 71, 17–30] β was found to depend on the rheogram shape factor for the fluid and the shear Reynolds number for the particle. The validity of this approach was supported by experimental data.     

Rill erosion: Exploring the relationship between experiments,modelling and field observations

Predicting rill erosion rates for a given discharge and slope minimally requires a model for rill hydraulics that allows the prediction of hydraulic parameters and a model for sediment detachment. Several relationships that describe rill hydraulics and/or sediment detachment within an eroding rill have been proposed and are incorporated into state of the art soil erosion models.In this paper a critical review of the theoretical concepts that are underpinning current rill flow and sediment detachment models is made in the light of recent experimental results. Approaches to define detachment–hydraulics relationship are generally based on developments in alluvial river hydraulics. However, experimental evidence to support the use of these concepts in models of rill erosion is scarce and recent experimental findings suggest that the basic assumptions used to model rill erosion are to some extent flawed. An analysis of empirically collected data on rill hydraulics conclusively shows that the empirical Manning equation does not hold for rill flow and should therefore not be used in rill erosion models. An empirical power law relationship relating velocity to discharges is much better in agreement with available experimental data, both for soils with and without rock fragments. In the absence of vegetation residue and/or other macroscopic, immobile elements such as rock fragments, total shear stress and unit length shear force can be used to predict soil detachment. The use of unit length shear force has the advantage that no information about rill geometry is necessary.The evidence for sediment load and rill flow detachment interaction is somewhat conflicting: the presence of a heavy sediment load appears to restrict rill flow detachment, but the exact form of the relationship between detachment rate and sediment load remains unclear. The effect of the presence of a sediment load on flow detachments under natural conditions is also limited by the nature of the detachment and transporting capacity relationships: on a rectilinear hillslope, transporting capacity increases much more rapidly with discharge than detachment capacity.We propose modifications to the theoretical formulations used in rill erosion models so that they are in better agreement with experimental evidence. Finally, we illustrate the potential of simplified models and conclude that the combination of empirical equations for flow detachment and rill hydraulics leads to results that are consistent with empirical data relating rill erosion rates to topography.     

Modeling the particle breakage of rockfill materials with the cohesive crack model

A practical combined finite–discrete element method was developed to simulate the breakage of irregularly shaped particles in granular geomaterials, e.g., rockfill. Using this method, each particle is discretized into a finite element mesh. The potential fracture paths are represented by pre-inserted cohesive interface elements (CIEs) with a progressive damage model. The Mohr–Coulomb model with a tension cut-off is employed as the damage initiation criterion to rupture the predominant failure mode occurs at the particle scale. Two series of biaxial tests were simulated for both the breakable and unbreakable particle assemblies. The two assemblies have identical configurations, with the exception that the former is inserted with CIEs and is breakable. The simulated stress–strain–dilation responses obtained for both assemblies are in agreement with experimental observations. We present a comprehensive study of the role of particle breakage on the mechanical behavior of rockfill materials at both the macroscopic and microscopic scales. The underlying mechanism of particle breakage can be explained by the force chain in the assemblies.     

粒子重复冲击破岩细观损伤及破碎特征试验研究

粒子冲击破岩作为一种新型辅助破岩技术已在钻探和油气开采等领域得到有效运用,其在极坚硬岩层中隧（巷）道辅助掘进方面具有广阔的应用前景。从试验角度研究了粒子的冲击次数、颗粒强度及冲击速度等因素对极坚硬花岗岩表面冲击坑损伤破碎特征的影响,对冲击坑三维形貌、岩石碎屑及坑内矿物破碎特征进行了定量分析。结果表明：成坑最大深度随冲击次数增加呈抛物线趋势增大,而成坑体积和坑顶面积随冲击次数增加呈线性增长;成坑体积随冲击速度呈先增大后减小的规律,该临界冲击速度约为82.5 m/s;冲击坑中心与外部的细观破碎机制差异导致岩石碎屑平均尺寸具有明显的双峰特征。从能量角度分析发现,成坑体积、坑顶面积和最大深度随粒子动能在双对数坐标系下线性增长。基于图像处理方法获得了冲击坑周边主要矿物内部裂纹分布随冲击速度和次数增加的分形维数变化规律。破岩效果表明,增大粒子冲击速度和冲击次数能有效扩大岩石冲击坑的损伤范围,但冲击速度的影响程度显著大于冲击次数。     

矿渣胶凝材料固化软土的力学性状及机制

利用矿渣胶凝材料固化软土,既可利用工业废渣,又能减少水泥的用量。以矿渣胶凝材料固化黏土、砂土二种软土。发现矿渣胶凝材料加固软土的效果远好于水泥、石灰,其9 %掺量的固化土28 d的无侧限强度达到2.0 MPa以上,普遍高于15 %掺量的水泥固化土,且其28 d固化土的软化系数普遍高于90 %以上,固化黏土后CBR值远高于同掺量的石灰固化土。X衍射结构分析表明,矿渣胶凝材料水化时产生的高强难溶的矿物晶体是其固化软土效果好的主要原因。因此,矿渣胶凝材料是一种性能优异的软土加固材料。