紫坪铺面板坝堆石料颗粒破碎试验研究

采用大型三轴仪对紫坪铺面板坝堆石料进行了单调和循环荷载下的固结排水剪切试验,研究了不同孔隙比情况下颗粒破碎及剪胀的变化规律。试验表明：（1）单调和循环荷载条件下,堆石料颗粒破碎率与塑性功之间存在一致的双曲线关系;（2）峰值应力处剪胀率与颗粒破碎率在半对数坐标中呈近似线性关系;（3）峰值应力处主应力比与相应的剪胀率呈近似线性关系,且上述结果受初始孔隙比的影响不大。研究成果有助于进一步了解堆石料的颗粒破碎特点,对建立复杂应力条件下考虑颗粒破碎和状态相关性的弹塑性本构模型,分析紫坪铺面板堆石坝汶川地震破损机制是十分有益的。     

DEM study on the effect of particle breakage on the macro- and micro-behavior of rockfill sheared along different stress paths

The behavior of crushable rockfill sheared along different stress paths is studied using the discrete element method. Rockfill particles are modeled as breakable agglomerates, and reasonable consistency is found between the predicted and experimental results. The simulation highlights the influence of agglomerate breakage on both the macro- and micro-behavior of the assembly. The evolution of both the number and mode of agglomerate breakage during shearing is significantly influenced by the confining pressure, deviator stress ratio, and loading direction. The relationship between the deviator stress ratio on the macro-scale and the deviator fabric on the micro-scale along different stress paths is explored.     

Discrete modelling of time‐dependent rockfill behaviour

This modelling study deals with the time‐dependent behaviour of rockfill media, which is of particular interest during the life of rockfill dams. Breakage of rock blocks and crack propagation are the main processes responsible for rockfill creep and collapse. The modelling procedure presented here is performed on two scales: on the rock block scale, where the grain is taken to be an assembly of rigid particles initially endowed with cohesive bonds, and on the rockfill scale, which is taken to involve a set of breakable grains interacting via contact and friction processes. The grain breakage process is described in term of a thermodynamically consistent damage interface model, where the damage is a gradual delayed process. This model was implemented in a non‐smooth contact dynamics code. The effects of the main parameters involved were analysed by performing numerical studies. The ability of the model to predict the creep behaviour of rockfill media is confirmed by presenting several simulations. Copyright © 2008 John Wiley & Sons, Ltd.     

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

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

A discrete element method (DEM)-based approach to simulating particle breakage

An approach for particle breakage simulation based on the framework of discrete element method was proposed in the current study. Convex polyhedron blocks were adopted as elementary particles for the complex particle shapes, and the variability of particle breakage strength is modeled using the invertible function method. Additionally, the traditional modified “Brazilian” criterion was adopted as the breakage criterion. Under the assumption that the eventual fractures within a particle can be determined according to the contact points and the centroid of the particle, once a target particle fulfilled the breakage criterion, it was cut into several fragments by a series of virtual cutting faces, which are consistent with the eventual fractures. With this, the production of local stress and the non-conservation of mass and volume can be avoided. A pre-defined fragmentation mode was also unnecessary for this approach. A series of numerical triaxial tests adopting this new presented approach was then conducted according to the configurations reported in the literature and comparisons made with experimental results. It revealed that while the presented approach is capable of reproducing the macroscopic shear responses and particle breakage characteristics of breakable particle assemblies, some fragmentation modes of particles such as surface grinding and corner abrasion cannot be captured using this approach, presenting an area for future investigation.

     

A random virtual crack DEM model for creep behavior of rockfill based on the subcritical crack propagation theory

The post-construction settlement of rockfill dams and high filled ground of airport, which is a phenomenon of much significance, is mainly caused by the time-dependent breakage of the rockfill material. In this paper, a random virtual crack DEM model is proposed for creep behavior of rockfill in PFC2D according to the theory of subcritical crack propagation induced by stress corrosion mechanisms. The bonded clusters are adopted to represent the rockfill particles so as to simulate their irregular shapes. Virtual cracks are set at the bonds of the clusters, and the length of the crack is considered as a random value, which leads the crushing strength of a single particle to follow the Weibull’s statistical model and the corresponding size rules. Oedometric creep tests for rockfill are simulated by using this proposed model. The results show that the model, validated preliminarily by some test data, can reflect qualitatively the creep mechanism as well as the size effects reasonably. Particles can develop various breakage patterns during creep, including global breakage, local breakage and even complex mixed breakage. The increase in stress levels and particle size will lead to an obvious growth of the creep strain and creep rate of the rockfill. The scale effects on the creep behavior of rockfill are analyzed through 35 specimens, and formulas including the effects of scales and stress levels are tentatively proposed.     

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

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

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

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

Discrete numerical modelling of rockfill dams

The aim of this study is to obtain quantitative information on the behaviour of rockfill used in embankment dams, and particularly on the influence of block breakage on the displacement field, from a numerical analysis using the Distinct element method. A methodology is set up to define the resistance of the 2D particles so that the same probability of breaking blocks may be reproduced as in a 3D material. The model uses the discrete element code PFC^2D (Itasca Consulting Group Inc., PFC2D (Particle Flow Code in Two Dimensions), Version 3.0, 2002) and considers breakable clusters of 2D balls. The different parameters are determined from experimental data obtained from laboratory tests performed on rock blocks. The model is validated by comparing the results of the simulation of shearing tests with actual triaxial tests on rockfill material published in the literature. The numerical analysis of block crushing in an actual dam is proposed in the last part of this paper. Copyright © 2006 John Wiley & Sons, Ltd.     

高聚物胶凝堆石料静力特性试验研究

采取加固措施可提高土体抗剪强度、减小变形等。采用聚氨酯高聚物胶凝堆石料加固土石坝以改良堆石料强度、变形特性。开展了聚氨酯胶凝堆石料三轴剪切试验研究,结果表明：聚氨酯胶凝堆石料偏应力-轴向应变曲线具有明显的非线性,抗剪强度提高;采用抛物线拟合偏应力-轴向应变曲线,可反映其应变软化和应变硬化,并得到初始弹性模量随聚氨酯含量的提高而提高;聚氨酯泡沫填充堆石颗粒间孔隙并将堆石颗粒黏结在一起,减少颗粒破碎、颗粒重排列,使堆石料剪缩特征减弱,当聚氨酯含量高时表现出剪胀;聚氨酯增加了堆石料准黏聚力,对内摩擦角无明显影响。通过有限元方法对堆石坝变形和稳定性进行分析,结果表明,聚氨酯加固堆石坝可减少边坡浅层滑动,增大滑弧深度,提高坝坡稳定安全系数。     

面板堆石坝垂直缝破坏下三维渗流场有限元模拟

采用改进节点虚流量法求解无压稳定渗流场,并引入无厚度的裂缝模型对破坏的垂直缝渗流行为进行模拟,得到面板堆石坝裂缝渗流问题的有限元分析方法,并编制Fortran程序。以某混凝土面板堆石坝为例,计算了面板单一垂直缝破坏和多条垂直缝破坏条件下的三维渗流场,得到不同条件下渗流场的水头分布、浸润线以及渗漏量,系统分析了面板堆石坝在垂直缝破坏条件下的稳定渗流场规律和特点。结果表明,该方法能对渗流逸出点和浸润线进行准确定位,还能很好地模拟面板垂直缝破坏对坝区渗流场的影响,可以为面板堆石坝的接缝设计提供参考。     

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.     

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

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

Numerical simulation of the seismic response of the Zipingpu concrete face rockfill dam during the Wenchuan earthquake based on a generalized plasticity model

The strong ground motion of the 2008 Wenchuan earthquake in China caused considerable damage to the Zipingpu concrete face rockfill dam (CFRD). The maximum crest settlement was approximately 1.0 m, and compressive failure and joint dislocations were observed in the face slabs. This damage has made it necessary to understand the damage pattern and safety of high CFRDs subjected to strong earthquake shaking, and the response of the Zipingpu CFRD during the Wenchuan earthquake can be used as a benchmark for this purpose. In this study, a 3D dynamic procedure was employed to simulate the dynamic responses of the Zipingpu CFRD. The rockfill materials were described using a generalized plasticity model, while the interfaces between the face slabs and cushions were modeled using zero-thickness interface elements that follow a perfect elasto-plastic stress–strain model in the tangential direction using Coulomb’s friction law. Dam deformation, face-slab stress, and face joint dislocations were simulated, and the results were compared with the field measurements. Using the generalized plastic model, the residual deformation of the dam during the earthquake could be directly obtained without being complemented by separate, semi-empirical procedures. The rockfill materials shrank to the center of the valley due to the strong shaking, causing crushing damage in the zone of the slabs. The dislocation of construction joints was also duplicated by the numerical procedure. The results of this study indicate that a 3D finite element procedure based on a generalized plasticity model can be used to evaluate the dynamic responses of CFRDs during strong earthquakes.     

高填石路堤蠕变本构模型及其参数反演分析与应用

结合高填石路堤工后沉降变形机理与工程特点，基于工程实测数据，提出了高填石路堤工后沉降蠕变变形的双曲线型三参数本构模型，并引进遗传算法与有限元分析理论，建立了该模型参数的反演分析方法。在此基础上，利用蠕变有限元分析手段，深入探讨了高填石路堤工后沉降的分析计算方法，并开发了相应的分析计算软件。结合某高填石路堤工程实践，探讨了高填石路堤双曲线三参数蠕变本构模型和高填石路堤工后沉降计算的应用方法。工程实例分析表明，提出的本构模型及其工后沉降计算方法简单，可以满足工程要求，初步建立了一种新的高填石路堤工后沉降计算方法。     

Modelling size effect on rock aggregates strength using a DEM bonded-cell model

Empirical evidence has shown that particle breakage affects the mechanical behaviour of granular materials. The source of this mechanism takes place at the particle scale, and the main consequence on the macromechanical behaviour is increasing compressibility. Due to the inverse correlation between particle size and particle crushing strength, coarse rockfill materials are particularly vulnerable to mechanical degradation due to particle breakage. However, such coarse materials do not fit in standard laboratory devices, and the alternative of large sample testing is usually unavailable or too expensive. Alternatively, recent works have proposed multi-scale approaches using the discrete element method (DEM) to carry out numerical testing of coarse crushable materials, although few studies have focused on size effects. This article presents the application of a DEM bonded-cell model to study particle size-strength correlation on angular rock aggregates. Each particle is modelled by a cluster of perfectly rigid polyhedral cells with Mohr–Coulomb contact law. Constant cell density within particles implies that the presence of potential fragmentation planes increases with size. Therefore, particle strength decreases with size. A comprehensive sensitivity analysis was carried out through 1477 particle crushing simulations in a given particle size. Based on published experimental data on calcareous rock aggregates, part of the simulations were used for calibration, and 97 additional simulations of a coarser size fraction were performed for validation. The results show a good agreement with the empirical data in terms of size effect and data scatter through Weibull statistics.

     

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

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

Energy distribution during the quasi-static confined comminution of granular materials

A method is proposed to calculate the distribution of energy during the quasi-static confined comminution of particulate assemblies. The work input, calculated by integrating the load-displacement curve, is written as the sum of the elastic deformation energy, the breakage energy and the redistribution energy. Experimental results obtained on samples subjected to compression stresses ranging between 0.4 and 92 MPa are used to calibrate the model. The elastic energy stored in the samples is obtained by simulating the compression test on the final particle size distributions (PSDs) with the discrete element method and by extracting the contact forces. A PSD evolution law is proposed to account for particle breakage. The PSD is related to the total particle surface in the sample, which allows calculating the breakage energy. The redistribution energy, which comprises the kinetic energy of particles being rearranged and the friction energy dissipated at contacts, is obtained by subtracting the elastic energy and breakage energy from the work input. Results show that: (1) at least 60% of the work input is dissipated by particle redistribution; (2) the fraction of elastic deformation energy increases, and the fraction of redistribution energy decreases as the compression stress increases; (3) the breakage energy accounts for less than 5% of the total input energy, and this value is independent of the compressive stress; (4) the energy dissipated by redistribution is between 14 and 30 times larger than the breakage energy.     

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

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

基于堆石坝竣工期沉降观测数据的材料非线性本构模型参数反演

为了高效和准确地确定堆石料的非线性本构模型参数,提出了基于响应面方法的参数反演方法。采用有限元方法数值模拟了堆石坝的分层填筑过程。建立了堆石坝变形观测点沉降的响应面函数,确定了多项式响应面函数的系数。根据堆石坝竣工期变形观测数据和确定的响应面函数,采用优化方法反演确定了堆石料本构模型参数。工程实际应用结果表明,该方法具有较高的计算效率,预测的堆石坝沉降变形与现场观测值吻合较好。