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.     

基于颗粒识别分析系统的碎屑流堆积物颗粒识别和统计方法研究

碎屑流堆积物颗粒识别和统计是碎屑流灾害的研究重点。文章基于图像处理孔隙（颗粒）及裂隙图像识别与分析系统（PCAS）,以贵州纳雍普洒村崩塌-碎屑流为例,结合纳雍崩塌堆积物粒径实测结果,通过阐释识别过程中阈值、孔喉封闭半径、最小孔隙面积的参数意义,研究PCAS软件在碎屑流颗粒识别与统计中的应用,并提出了颗粒识别时这些参数的选取方法。分析结果表明:（1）PCAS能自动准确地识别碎屑流堆积物颗粒与孔隙,相比人工计数更精细,所识别堆积物各区小颗粒比重较大,0～2 m颗粒粒径各区占比都在50%以上;（2）当阈值为170（像素）时能获得精细的二值图像,颗粒与孔隙得到了准确地区分;（3）不同参数取值下获得堆积物颗粒粒径分布结果不同,碎屑流堆积物颗粒识别宜采用较大的孔喉封闭半径和较小孔隙面积,当二者比值为3/30（像素）时能更好地反应颗粒粒径分布情况;（4）PCAS具有较高的可行性,统计结果显示,各粒径含量变化趋势与人工统计相近,两种统计方法各粒径占比、分布规律基本吻合,说明利用PCAS可以实现对崩塌碎屑流颗粒粒径分布的高效便捷分析。     

粗颗粒悬移质测验误差分析

本文根据流速与含沙量分布,研究了悬移质积点法对粗颗粒含沙量测验误差.首先指出无论从输沙率还是从悬浮高度看,粗颗粒悬移质集中于河底,致使水文测验中的一些测验方法(如一点法、二点法、三点法等)带有较大的误差.进而分析了各种测法对粗颗粒的相对误差,并且随着颗粒变粗,误差愈来愈大,以致不能接受.这种情况在水库下游河道冲刷过程中,会在一定程度上出现,必须重视研究和调整测验方法.最后通过与长江委水文局在沙市及监利两水文站试验结果对比,证明本文理论上的分析是可信的.     

苏州澄湖湖底硬粘土粒度特征及成因初探

澄湖是太湖平原上的一个浅水湖泊,湖底十分平坦,主要由硬粘土组成,同时发育有NW—SE横穿澄湖的古河道和泥沼洼地。在接近湖泊中心部位的古河道中打钻获得了4 m长的岩芯,并进行了详细的粒度分析。结果表明,该岩芯下部的硬粘土主要由细粉砂、中粉砂和粘土组成,分别占26.7%~38.1%,26.9%~34.3%和19.7%~26.2%,可命名为粘土质粉砂;众数粒径为10~50 μm,平均占57.9%;平均粒径为6.3~6.9。将其与上部的湖泊沉积、镇江下蜀黄土和中国北方各地黄土的粒度比较表明,苏州澄湖湖底的硬粘土原本也是风尘沉积。     

全尾砂-废石混合回填膏体渗透特性试验研究

塌陷区混合回填膏体是由全尾砂、废石人工制备得到的复合宽粒级散体,其渗透性能直接影响到塌陷区回填体稳定性及次生灾害发生的可能性。采用自制的渗透性试验装置,研究了废石含量、废石粒径对混合回填体渗透系数的影响,并对回填体固有特征参数（粒径小于0.075、0.02 mm颗粒含量、不均匀系数 、平均粒径）与渗透系数之间关联性进行了研究。试验结果表明,回填膏体渗透系数随废石掺量及粒径增大而增加,细粒效应显著。回填体中粒径小于0.075、0.02 mm颗粒含量对其渗透性影响至关重要,与渗透系数呈负指数关系。渗透系数随 值增大而增大。当 20时,渗透系数增加趋于稳定。给出了混合回填体渗透系数定量方程,其结果与试验结果较为一致。     

崩积混合体直剪试验与PFC2D数值模拟分析

通过野外崩积混合体结构特征分析,考虑原样级配,开展对崩积混合体重塑样的室内大尺度直剪试验研究,结合PFC2D颗粒离散元仿真试验,分析了不同含石量情况下崩积混合体变形与强度变化规律及内在机理。将数值试验的成果与室内试验结果进行对比分析。结果表明,随着含石量的增加,崩积混合体的应变硬化效应显著。低法向应力下崩积混合体表现为剪胀,高法向应力下则表现为剪缩;当含石量低于40%时,崩积混合体的力学性质由土体控制,超过80%以后,其力学性质基本上完全由块石控制。     

考虑土石料颗粒破碎和密度变化的次塑性本构模型建模方法

试验和现场观测都发现高土石坝的坝料在高压及湿化作用下会发生显著的颗粒破碎,颗粒破碎会改变土石料的级配曲线和密度,因而影响其后继力学行为,因此,土石料的颗粒破碎是当前岩土工程领域的研究热点。为了模拟土石料在高压及湿化作用下发生显著的颗粒破碎现象,以及循环加载中的颗粒破碎与应力诱导各向异性随动硬化共同影响下土石料的变形,本文提出了一个建模方法,考虑土石料颗粒破碎和密度变化的影响。所建议的次塑性本构模型在经过试验资料的验证后可用于动力有限元数值计算。     

锦屏大理岩蠕变损伤演化细观力学特征的数值模拟研究

岩石的蠕变损伤和断裂是岩石流变效应的重要表现形式,但其损伤演化过程往往难以直观观测,为此,采用二维颗粒流数值模拟方法(PFC2D)对岩石的蠕变损伤和断裂的细观力学机制进行了分析。在锦屏大理岩室内试验基础上,利用颗粒流应力腐蚀模型(PSC),建立了能反映其短期和长期强度特征的柱状岩样数值模型,并开展了大量数值试验。结果表明,蠕变损伤的演化过程与暂态的损伤演化过程具有明显的差异。在岩样蠕变损伤过程中,其内部微裂纹多沿加载方向开裂且分布均匀,先快速增加再稳定扩展,最后则发生快速断裂。当荷载较小时,岩样宏观上呈现劈裂破坏特征,当荷载较大时,岩样呈现剪切破坏特征。在岩样蠕变损伤初始和稳定演化阶段的前期,荷载大小对岩样的损伤演化过程影响不大;在稳定演化阶段的后期至断裂过程中,低荷载下岩样的损伤增速比高荷载下快。     

盾构切刀顺次破岩机制的数值模拟研究

在对盾构切刀顺次破岩实际工况提出合理简化的基础上,从岩土细观角度出发,采用颗粒离散元法建立了切刀破岩的二维数值模型,研究了两把盾构切刀顺次作用下的破岩机制和影响因素,并通过试验对切刀破岩过程中切屑堆积现象及破岩力学特性进行了验证。研究表明,刀具尖端的破坏作用最为明显;切削岩石时切削力随着切削行程不断波动,水平切削力大于垂直切削力;不同于单刀切削,切刀顺次作用时前刀刮过的岩面留下了大量残余裂纹,使得后刀所受的切削力减小;从切削性能来看,随着切削深度的增加,岩石破碎块度不断增大,切削力和裂纹数迅速增加,说明切深与切削力密切相关;切削试验观测到了切屑在前刀面堆积和切削力的波动现象,与数值模拟具有较好的一致性。     

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.