堆石料大三轴试验的细观模拟

采用三维颗粒流计算程序,对堆石料的大三轴排水剪切试验过程进行了数值模拟。在数值模型的生成过程中,使堆石料集合体达到规定密度和级配;采用动态松弛算法迭代求解,并引入粒间阻尼,以吸收颗粒填装的多余动能使其稳定;颗粒间引入黏结力来提高试样的峰值强度;确定堆石料的细观力学参数,拟合室内试验应力应变曲线。模拟不同围压下堆石料的试验,得到的应力-应变曲线与室内试验基本一致,说明颗粒流方法可以较好的模拟堆石料的大三轴试验过程。由于没有考虑颗粒形状及颗粒破碎的影响,造成大应变剪胀偏大。     

堆石料碾压试验的颗粒流模拟新方法

针对堆石料现场碾压试验,研究基于颗粒流的仿真模拟新方法。首先,以随机多边形颗粒簇形式,建立堆石料的颗粒流模型,充分体现了堆石料的级配、细观组构和破碎效应,从而可以较为真实地模拟堆石料在碾压过程中的宏、细观力学响应;然后,提出振动碾压荷载及碾压过程的模拟新方法,实现振动碾压荷载及碾压过程的仿真模拟;最后,结合水布垭主堆石料的现场碾压试验,利用该方法进行模拟验证,从而论证该方法的可靠性。该方法可为仿真模拟堆石料碾压试验以及研究碾压效应对堆石料宏观力学特性的影响奠定基础,具有重要工程指导意义。     

颗粒劣化效应的堆石料流变细观数值模拟

堆石料流变的主要机理是由于水位变化、降雨入渗、日晒雨淋等环境因素导致堆石料发生性质劣化,与此同时,颗粒发生高接触应力-破碎和重新排列-应力释放、调整和转移,这一过程由于颗粒的持续劣化而不断重复。在考虑颗粒破碎的堆石体不连续变形分析方法(SGDD)中,引入颗粒强度劣化模型反映颗粒强度随外界环境的持续劣化。应用该方法进行堆石料三轴流变数值试验,模拟结果与室内试验所观察到的规律一致,表明考虑流变效应的SGDD方法抓住堆石料流变的主要机理,适合模拟堆石料的流变变形这一复杂的、非线性演化问题。数值试验结果表明,堆石料随外界环境的劣化程度、劣化速率、母岩强度对宏观流变变形有较大影响。     

堆石料的临界状态与考虑颗粒破碎的本构模型

高应力水平时堆石料的颗粒破碎对其强度和变形机制有重要影响。临界状态土力学理论对重塑土的应力-应变关系的描述较为成功,但目前颗粒破碎对堆石料的临界状态的影响及其数学描述鲜有研究。对堆石料进行了固结应力从0.4 MPa到4 MPa的18组固结排水和固结不排水常规三轴压缩试验,以及6组等向压缩试验。试验结果表明：在排水条件和不排水条件下,不同的固结应力试样都趋于临界状态;堆石料的临界状态在q-p′平面和e-lgp′平面均为非线性变化。基于此试验结果,通过引入状态参数,在广义塑性力学的理论框架下,建立了考虑颗粒破碎的堆石料本构模型,并给出了模型参数的确定方法。与长河坝料的试验进行了对比,结果表明所建议的本构模型可以较好地模拟堆石料从低围压0.4 MPa到高围压 3.5 MPa下的应力-应变特性     

堆石料单轴流变试验的颗粒流模拟

岩石断裂力学的亚临界裂缝扩展理论认为微裂缝扩展可导致岩石破碎,即岩石颗粒破碎具有时间效应。根据亚临界裂缝扩展理论,提出了考虑微裂缝扩展导致堆石颗粒破碎时间效应的数值流变模拟新方法,并进行了考虑颗粒不同典型破碎模式的单轴流变颗粒流数值试验。在对比数值与室内流变试验曲线的基础上,分析了数值流变过程中颗粒破碎情况与颗粒体内部结构发展过程等。研究成果表明,两种试验手段得到的堆石流变的发展趋势基本一致,由微裂缝扩展引起的颗粒延时破碎是堆石流变的主要原因之一,深化了对堆石料变形机制的认识。     

堆石料单粒强度尺寸效应的颗粒流模拟方法研究

围绕堆石料单粒强度尺寸效应的颗粒流模拟方法展开研究。首先,基于FISH二次开发建立了堆石料的随机不规则单粒模型,充分考虑堆石料的形状特征和破碎现象;然后,建立了堆石料单粒强度尺寸效应的等效模拟方法,以单粒强度随其粒径的变化规律为基础,推导了堆石料模型中细观黏结强度与堆石料等效粒径的负指数经验公式;其次,基于建立的数值模型对堆石料的室内单粒压缩试验进行仿真模拟,验证数值模型的正确性和合理性,并对较大粒径堆石料的单粒强度进行模拟预测,突出数值试验的优势;最后,基于建立的数值模型对相同粒径不同形状特征堆石料的单粒强度分布特征进行模拟研究。研究结果表明：（1）堆石料内部缺陷含量和尺寸随粒径增加对其单粒强度所产生的尺寸效应,可通过堆石料模型中细观强度参数随粒径折减进行等效模拟;（2）形状特征对堆石料的破裂机制具有重要影响,方形颗粒为压剪破裂,单粒强度较高,而随机不规则颗粒和圆形颗粒为拉剪或劈裂,单粒强度相对较低;（3）拉剪或劈裂条件下,堆石料形状越不规则,其单粒强度的离散程度越高,反之则离散程度越低。相关研究成果可为进一步研究荷载作用下堆石体内各粒径段堆石料的破碎量奠定基础,从而更加真实地反映堆石体的级配演化规律。     

砂土力学性质的细观模拟

基于颗粒流理论,对砂土的室内双轴试验和砂土剪切带形成与发展进行了数值模拟,基本再现了砂土试样应力-应变关系。主要研究了颗粒粒径、颗粒摩擦系数等细观参数变化时试样宏观性质的变化情况。对比分析了室内试验和颗粒流数值模拟试验的砂土剪切带特征。同时对比研究了模型试样颗粒粒径、颗粒刚度和摩擦系数等细观参数的变化,对剪切带的形成与发展过程的影响。结果表明：通过颗粒流数值模型试验可以有效模拟砂土剪切带的形成与发展机理。     

考虑颗粒破碎效应的堆石料静动力本构模型

为合理反映颗粒破碎对堆石料力学特性的影响,基于试验结果分析,得出了堆石料在压缩和剪切作用下的颗粒破碎特性规律。通过引入压缩破碎和剪切破碎的相关参数,借鉴已有本构模型的合理定义,吸收临界状态理论和边界面理论的优点,发展了考虑颗粒破碎和状态相关的堆石料静动力统一弹塑性本构模型,并阐述了模型参数的确定方法。该模型不仅能够反映堆石料在静力荷载作用下的低压剪胀、高压剪缩、应变软化和硬化等特性,还能够反映在循环荷载作用下应力-应变的滞回特性和残余变形的累积效应。最后为验证模型的合理性,分别对堆石料的静力三轴和循环三轴试验进行了数值模拟预测,结果表明：模型预测与试验数据吻合良好,所提出的本构模型能够合理地描述颗粒破碎对堆石料静动力变形特性的影响。     

基于颗粒流的堆石料三轴剪切试验研究

近年来频发的地质灾害使得深入研究大坝堆石体受力变形机制显得愈发迫切和重要。针对堆石颗粒具有形状极其不规则、咬合作用突出的特点,编写了一个简单易用的程序来随机生成形状不规则的数值颗粒。以正四面体为核,按照晶胞繁衍的方式,生成堆石的仿真颗粒集合体。模拟排水剪切条件下的大三轴试验,通过分析数值试验过程中细观参量的变化规律,探讨三轴剪切条件下堆石体变形的细观机制。分析表明,按照晶胞繁衍法生成的数值颗粒与圆颗粒相比能更好地模拟堆石,由颗粒簇形成的咬合力能形成更真实的力-变形关系;数值试验得到的应力-应变曲线和体变曲线与室内三轴试验结果基本一致;颗粒的破碎速率是联系堆石体宏、细观力学性质的重要纽带,通过分析堆石体颗粒破碎4个阶段细观参量的变化规律,深化了对堆石体变形机制的认识。     

粗粒料颗粒破碎数值模拟研究

粗粒料在外力作用下存在明显的颗粒破碎特性,研究颗粒破碎过程是当前研究的热点问题之一。基于粗粒料单颗粒破碎机制,考虑单颗粒破碎强度与直径的变化规律,采用非线性接触H-Z模型和密度控制法建立了粗粒料颗粒破碎数值模型。开展粗粒料双轴剪切试验数值模拟研究,并与室内试验结果进行对比分析。研究表明：建立的粗粒料颗粒破碎数值模型能够较好地模拟粗粒料偏应力与轴向应变和体积应变与轴向应变的关系;数值模拟获得的粗粒料颗粒破碎率与室内试验结果基本一致;去除试样制备过程和固结过程引起的颗粒破碎,不同围压条件下的颗粒破碎率归一化后基本重合,且可以近似地采用双曲线函数进行拟合。颗粒破碎率随着围压的增大,逐渐增大,试验级配趋于Einav提出的颗粒破碎的最终级配（分维数等于2.6）     

堆石料三维形状量化及其对颗粒破碎的影响

筑坝堆石料由于尺寸较大,必须对其按一定比例缩尺后才能用来开展室内三轴试验。但缩尺前后颗粒形状难免会有差异,如何评价颗粒形状变化对堆石料变形特性的影响是十分重要的。引入了高精度的三维激光扫描技术对紫坪铺面板坝筑坝堆石料2.5～5、5～10、10～20 mm以及20～40 mm 4个粒径组的颗粒进行了空间形状分析,在此基础上进一步开展了单一粒组的三轴试验,研究了4个粒组的颗粒形状指标与颗粒破碎率的相关性。试验表明,紫坪铺堆石料颗粒破碎率随着其平均球度的增加而减小,并且呈近似半对数线性关系;随着围压的增加,颗粒形状对颗粒破碎的影响逐渐减弱,颗粒强度的影响逐渐增大。紫坪铺堆石料的颗粒强度随着颗粒尺寸的增加逐渐减小,但其破碎率反而随着颗粒尺寸的增加而逐渐减小,主要是因为试验所采用的紫坪铺堆石料颗粒尺寸越小时,其形状越不规则。因此,研究缩尺效应对颗粒破碎率的影响时,要同时考虑颗粒尺寸和颗粒形状。     

Study on the particle breakage of ballast based on a GPU accelerated discrete element method

Breakage of particles will have greatly influence on mechanical behavior of granular material(GM)under external loads,such as ballast,rockfill and sand.The discrete element method(DEM)is one of the most popular methods for simulating GM as each particle is represented on its own.To study breakage mechanism of particle breakage,a cohesive contact mode is developed based on the GPU accelerated DEM code-Blaze-DEM.A database of the 3D geometry model of rock blocks is established based on the 3D scanning method.And an agglomerate describing the rock block with a series of non-overlapping spherical particles is used to build the DEM numerical model of a railway ballast sample,which is used to the DEM oedometric test to study the particles’breakage characteristics of the sample under external load.Furthermore,to obtain the meso-mechanical parameters used in DEM,a black-analysis method is used based on the laboratory tests of the rock sample.Based on the DEM numerical tests,the particle breakage process and mechanisms of the railway ballast are studied.All results show that the developed code can better used for large scale simulation of the particle breakage analysis of granular material.     

颗粒破碎效应对堆石坝应力变形的影响

颗粒破碎对堆石体的变形有重大影响,颗粒破碎除与母岩强度有关外,与应力状态密切相关,高堆石坝堆石体的大部分变形正是高应力场下颗粒破碎的反映。沈珠江院士提出的南水双屈服面模型虽可以较全面地反映堆石体的变形特性,但未充分考虑颗粒破碎引起的剪缩特性,因此坝体变形计算结果往往偏小,适用于中低堆石坝的分析中。本文的计算方法在南水双屈服面模型的基础上,考虑了堆石体在一定应力状态下颗粒破碎引起的应变,较真实地反映了堆石体的变形特性,可为堆石坝尤其是高堆石坝的应力变形研究提供更为合理的结果。     

应力路径和干湿状态对堆石料颗粒破碎的影响研究

堆石料颗粒破碎是引起高土石坝变形的重要因素。在大坝填筑、蓄水期,堆石料的应力路径和干湿状态均是变化的。通过大型三轴试验,系统地研究了应力路径和干湿状态对堆石料颗粒破碎规律的影响。试验结果表明：（1）相同初始条件下,按不同应力路径达到同一轴向应变停机时测定的颗粒破碎率是不同的,等围压?3试验产生的颗粒破碎最大,等平均主应力p试验的次之,等最大主应力?1试验的最小,但不同应力路径下的颗粒级配演化规律是一致的。（2）相同初始条件下,湿样的颗粒破碎率明显高于干样,且二者的差距随着围压的增大而增大,不同干湿条件下的颗粒级配演化规律同样是一致的。（3）建立的考虑母岩强度的颗粒破碎率与塑性功的关系可以较好地统一描述不同应力路径及干湿状态下的颗粒破碎。该研究成果可为建立复杂应力路径及干湿变化条件下考虑颗粒破碎的堆石料弹塑性本构模型提供依据。