PFC2D模拟裂隙岩石裂纹扩展特征的研究现状

二维颗粒流数值模拟（PFC2D）是目前研究裂隙岩石裂纹扩展特征的重要手段。在大量已有相关研究文献的基础上做了以下分析和总结:从颗粒接触本构模型、细观参数的标定和裂隙模拟方法3个方面对当前PFC2D的主要模拟方法进行了总结;根据PFC2D模拟裂隙岩石裂纹扩展特征的研究现状,重点对单裂隙、断续双裂隙岩石在不同加载方式下的裂纹扩展特征进行了深入总结。在此基础上,指出当前研究中存在如下不足:裂隙岩石的PFC2D模型未考虑断裂韧度是否符合实际、平行黏结模型模拟结果与室内试验结果存在差异、模拟裂隙与真实裂隙存在差异。结合研究中存在的不足,提出了相应的解决办法并进行展望,以期有助于裂隙岩石PFC2D模拟方法的发展。     

Microparameter Prediction for a Triaxial Compression PFC3D Model of Rock Using Full Factorial Designs and Artificial Neural Networks

This paper presents an integrated approach that predicts the microparameters of the particle flow code in three dimensions (PFC3D) model in triaxial compression simulations. The new approach combines a full factorial design (FFD) with an artificial neural network (ANN). The ANN model maps the input factors (triaxial compressive strength, Poisson’s ratio, and Young’s modulus) onto output variables, which are microparameters that affect the macroscopic responses in a PFC3D model. Emphasis is placed on data collection and optimization of the ANN model using FFD. The data for training and testing the ANN model were obtained from laboratory experiments and numerical simulations of a PFC3D model according to the principles of FFD. Using a backpropagation artificial neural network (BPNN) optimized with FFD principles, the object of the current study (to reliably predict the microparameters for a PFC3D model) has been achieved because the predicting data obtained by the BPNN model were in excellent agreement with the testing data.     

A Discrete Element Model for Predicting Shear Strength and Degradation of Rock Joint by Using Compressive and Tensile Test Data

A discrete element model is proposed to examine rock strength and failure. The model is implemented by UDEC, which is developed for this purpose. The material is represented as a collection of irregular-sized deformable particles interacting at their cohesive boundaries. The interface between two adjacent particles is viewed as a flexible contact whose constitutive law controls the material fracture and fragmentation properties. To reproduce rock anisotropy, an orthotropic cohesive law is developed for the contacts, which allows their shear and tensile behaviors to be different from each other. Using a combination of original closed-form expressions and statistical calibrations, a unique set of the contact microparameters are found based on the uniaxial/triaxial compression and Brazilian tension test data of a plaster. Applying the obtained microparameters, joint specimens, made of the same plaster, are simulated, where the comparison of the obtained results to laboratory data shows a reasonable agreement.     

Calibration of a constitutive model using genetic algorithms

Before any constitutive model can be used in a numerical procedure, the model needs to be calibrated using laboratory test results. The traditional calibration techniques use stress and strain levels at certain states that a material undergoes during certain types of laboratory tests. Sometimes this method of calibrating a constitutive model fails to capture the overall behavior of a material, i.e. behavior at every point in stress/strain path. In this paper, we have shown how a random search technique, genetic algorithm (GA), can be used to calibrate constitutive models. The advantages of using GA are that it considers the overall behavior of a material, not the behavior at some specific states as the traditional method does, and it can work with many types of laboratory tests. The concept is applied to calibrate the hierarchical single surface (HiSS) δ_l model for geologic materials. Three cases have been studied where the difference is in the type of test data used to calibrate the model. These three different test data are: (1) simulated conventional test data, (2) simulated cyclic test data, and (3) real test data. A comparison of two different cross-over schemes, one-point and six-point, has been made.     

不同墙体位移方式下被动土压力的颗粒流模拟

用二维颗粒流程序（PFC2D）对各种移动模式的刚性挡墙被动土压力进行了模拟。再根据已有的模型试验建立颗粒流模型,分别让刚性墙体朝土体一侧平移、绕墙体顶部转动和绕墙体底部转动,给出了这三种情况下墙后被动土压力的分布图形以及相应的竖向应力的分布图形,同时给出了墙后被动土压力随墙体位移的变化规律。最后,将计算结果与实验结果进行了对比表明,颗粒流数值模拟方法不仅使用方便,而且能较好地模拟墙后被动土压力的分布规律。     

Application of unconstrained optimization and sensitivity analysis to calibration of a soil constitutive model

Large sets of soil experimental data (field and laboratory) are becoming increasingly available for calibration of soil constitutive models. A challenging task is to calibrate a potentially large number of model parameters to satisfactorily match many data sets simultaneously. This calibration effort can be facilitated by optimization techniques. The current study aims to explore systematic approaches for exercising optimization and sensitivity analysis in the area of soil constitutive modelling. Analytical, semi‐analytical and numerical optimization techniques are employed to calibrate a multi‐surface‐plasticity sand model. Calibration is based on results from a number of drained triaxial sample tests and a dynamic centrifuge liquefaction test. The analytical and semi‐analytical approaches and associated sensitivity analysis are applied to calibrate the model non‐linear shear stress–strain response. Thereafter, model parameters controlling shear–volume coupling effects (dilatancy) are calibrated using a solid–fluid fully coupled finite element program in conjunction with an advanced numerical optimization code. A related sensitivity study reveals the challenges often encountered in optimizing highly non‐linear functions. Overall, this study demonstrates applicability and limitations of optimization techniques for constitutive model calibration. Copyright © 2003 John Wiley & Sons, Ltd.     

陇西黄土三轴剪切过程微观变化研究

为了解黄土三轴剪切过程中的微观变化,采用PFC3D建立黄土三轴试验模型,模拟围压分别为0、50、150、300 kPa的三轴剪切试验,并与室内试验进行对比分析。分析结果表明：PFC3D能够较好地模拟出不同围压下从开始到破坏到残余变形整个过程中应力-应变的变化规律,且位移场及接触应力场的变化规律与室内三轴试验宏观现象较一致;发现弹性模量、泊松比及峰值强度与数值模型中微观参数有着密切的联系,如法向刚度kn控制试样宏观弹性模量,kn /ks值控制泊松比,摩擦系数控制峰值强度;通过体应变-轴向应变曲线发现,随着围压的增加应变能增大,试样呈现出由体积膨胀到体积减缩的变化规律。其研究结果为进一步探究黄土的应力-应变性状及抗剪强度特性提供参考。     

砂土中等截面桩的上拔机制分析

桩的上拔承载性能的宏观力学现象与桩周土细观结构变化相关,应用细观力学的颗粒流（PFC2D）数值模拟方法对承受上拔荷载作用的桩基进行了分析,数值模拟了上拔荷载作用的桩及桩周土的细观力学特征,研究了桩侧摩阻力的分布、桩周土剪切带的形成过程,较好的再现了桩的荷载－位移关系的实验结果,并与宏观物理实物试验的位移实测结果作了对比分析。分析了土体中剪切带形成过程中的颗粒间的细观变化及其形成过程,当上拔荷载达到极限时,上拔桩的剪切带形成原因是密砂的应变软化效应;颗粒流数值模拟的颗粒接触力与实物物理试验桩侧摩阻力是同一的,数值模拟的荷载-位移曲线与实物物理试验的荷载-位移曲线一致;桩侧摩阻力、桩上拔过程中剪切带的形成过程、桩上拔荷载-位移关系与颗粒流数值模拟的颗粒分布、速度、接触力的细观参数的变化密切相关。桩承受荷载过程中土颗粒细观结构变化的颗粒流仿真,是关于细观力学特征与宏观力学响应的初步研究。     

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.     

Feasible packing of granular materials in discrete-element modelling of cone-penetration testing

This paper explores how the discrete-element method (DEM) was found to play an increasingly important role in cone penetration test (CPT) where continuum-mechanics-based analysis tools are insufficient. We investigated several crucial features of CPT simulations in the two-dimensional DEM. First, the microparameters (stiffness and friction) of discrete material tailored to mimic clean, saturated sand, which is used in cone-penetration tests, were calibrated by curve-fitting drained triaxial tests. Then, three series of cone-penetration simulations were conducted to explore (1) top boundary conditions, (2) reasonable size of discrete particles at different initial porosities, and (3) limit initial porosity of the model for a balance between accurate representation and computational efficiency. Further, we compared the cone-penetration resistance obtained in the laboratory and numerical simulations for the range of relative densities.     

Influence of numerical precision on the calibration of AEM-based groundwater flow models

Groundwater modelers have embraced the use of automated calibration tools based on classical nonlinear regression techniques. While clearly an improvement over trial-and-error calibration, it is not clear to what extent these popular inverse modeling tools yield accurate parameter sets for groundwater flow models. The impact of model configuration and precision upon automated parameter estimation is also unclear. An extensive set of numerical experiments was performed to explore the influence of model configuration on the calibration of a regional groundwater flow model developed using the analytic element method. The results provided insight into the manner in which the specified level of model precision and the location of observation points influence the results of inverse modeling based on nonlinear regression. While the importance of these issues is application-specific, obtaining an accurate model calibration for the case study required both a careful placement of test observations and a greater-than-anticipated level of model precision. The required level of model precision for calibration was more than necessary to produce an acceptable flow solution.     

Strategy for consistent model parameter calibration for soft soils using multi-objective optimisation

Constitutive models for soft soils require a large number of parameters to model the complex material response. One set of parameters should capture the soil response for various laboratory & in situ stress paths. This requires a new method to calibrate a consistent set of model parameters using test data from different load paths of classic geomechanical tests. The feasibility of the proposed method is demonstrated with the recently developed CREEP-SCLAY1S model. After a sensitivity analysis that highlights the model parameters for optimisation, an optimisation process for two different configurations is designed. The latter is successfully verified against artificially generated laboratory data.     

考虑转动阻抗的粗粒土离散元模拟

转动阻抗被定义为作用颗粒接触上的一对对称力偶,用来抵抗颗粒之间的相互转动。将转动阻抗引入到离散元模拟中是对传统离散单元法的重要改进。开发出考虑颗粒转动阻抗的接触模型,并将其嵌入到PFC2D中,利用该模型进行粗粒土的双轴剪切数值模型试验,研究剪切过程中转动阻抗对粗粒土的宏细观力学性质的影响。结果显示,在宏观方面,颗粒转动阻抗对粗粒土的宏观力学行为（应力-应变及体应变-轴应变行为）有重要的影响,随着转动阻抗的增加,粗粒土的剪切强度和最大摩擦角随之增加,这与已有的研究结果一致,证明所建模型是可靠的;在微观方面,考察转动阻抗对粗粒土内部微观结构的影响发现,随着转动阻抗的增加,粗粒土的内部的接触数目减少,而粗粒土的剪切强度增加,表明转动阻抗能够提高粗粒土力链网络的稳定性,同时发现随着转动阻抗的增加,粗粒土的各向异性增加主要是强力链各向异性的增加,说明转动阻抗增强了强力链的传递力的能力以及抵抗力链屈曲破坏能力。数值模拟结果表明,增加颗粒转动阻抗,粗粒土出现组构与轴应变非共轴的现象。     

Discrete Element Modeling of Drop Tests

A discrete element code with impact model has been developed and calibrated to simulate the dynamic behavior of rock materials, with special regard to rock fragmentation upon impact during rock-fall analysis. The paper summarizes the discrete element code, the calibration algorithms developed to identify the model microparameters, and the impact model. Experimental work on drop tests is then used to validate the code on modeling impact fragmentation. It has been found that the developed discrete element code and impact model can reasonably simulate rock fragmentation in drop tests. The use of the discrete element code and impact model can provide good reference results in evaluating impact fragmentation in rock-fall analysis.     

Numerical modelling of rheological properties of landslide debris

Debris flow has caused severe human casualties and economic losses in landslide-prone areas around the globe. A comprehensive understanding of the morphology and deposition mechanisms of debris flows is crucial to delineate the extent of a debris flow hazard. However, due to inherent complex field topography and varying compositions of the flowing debris, coupled with a lack of fundamental understanding about the factors controlling the geomaterial flow, interparticle interactions and its final settlement resulted in a limited understanding of the flow behaviour of the landslide debris. In this study, a physical model was set up in the laboratory to simulate and calibrate the debris flow using PFC, a distinct element modelling-based software. After calibration, a case study of the Varunavat landslide was taken to validate the developed numerical model. Following validation with an acceptable level of confidence, several models were generated to evaluate the effect of slope height, slope angle, slope profile, and grain size distribution of the dislodged geomaterial in the rheological properties of debris flow. Both qualitative and quantitative analysis of the landslide debris flow was performed. Finally, the utility of retaining wall and their effect on debris flow is also studied with different retaining wall positions along the slope surface.

     

砂性土成拱应力释放特性的模型试验及数值模拟研究

通过室内模型试验研究了卸荷尺度及孔隙率对砂性土成拱应力释放特性的影响,并应用离散元分析软件PFC2D进一步对颗粒摩擦性质影响、成拱后应力场以及位移场变化规律和应力释放分区进行了细观研究。研究发现：（1）数值模拟对卸荷尺度影响因素的研究结果与室内试验一致,应用PFC2D研究砂性土应力释放特性是可行的。（2）应力释放区内竖向应力释放大于应力加载区,而水平应力释放小于应力加载区;侧压力系数在应力释放区内距离卸荷位置越近增量越大,而在应力加载区距离卸荷位置越近增量越小。（3）砂性土摩擦性质会减小土体成拱对应力释放以及侧压力系数变化的影响。（4）应力释放随着卸荷尺度及孔隙率的增大而增大;对于侧压力系数,在应力释放区内随着卸荷尺度或孔隙率的增大而增大,而在应力加载区随着卸荷尺度或孔隙率的增大而减小。     

孔隙率、级配参数对粉土双轴压缩性状影响的颗粒流分析

针对颗粒流法中颗粒几何参数等对土体宏观力学性质的影响进行了数值模拟,并将数值模拟结果与室内模型试验验证进行了对比。对二维、三维下土体级配及孔隙率进行了理论分析,并提出二维数值试样颗粒数目及孔隙率转换公式。采用PFC2D软件进行了6组参数验证及24组交叉对比（含4组级配、3组孔隙率及2类接触模型）双轴压缩试验。结果表明：孔隙率对土体力学性质的影响大于级配;随着孔隙率增加,土体强度峰值逐渐降低,对应轴向应变增大,体应变则由剪胀性逐渐向剪缩性过渡。当孔隙率较低时,采用线性及Hertz-Mindlin (H-M)接触模型得到的模拟结果均较接近;而当孔隙率升高时,建议采用H-M接触模型,以更好地反映土体到达强度峰值后产生的后续剪胀效应。     

A Study of the Failure Mechanism of Planar Non-Persistent Open Joints Using PFC2D

Particle flow code 2D (PFC2D) was adopted to simulate the shear behavior of rocklike material samples containing planar non-persistent joints. Direct shear loading was conducted to investigate the effect of joint separation on the failure behavior of rock bridges. Initially calibration of PFC was undertaken with respect to the data obtained from experimental laboratory tests to ensure the conformity of the simulated numerical models response. Furthermore, validation of the simulated models were cross checked with the results of direct shear tests performed on non-persistent jointed physical models. Through numerical direct shear tests, the failure process was visually observed, and the failure patterns were found reasonably similar to the experimentally observed trends. The discrete element simulations demonstrated that the macro-scale shear zone resulted from the progressive failure of the tension-induced micro-cracks. The failure pattern was mostly influenced by joint separation, while the shear strength was linked to the failure pattern and failure mechanism. Furthermore, it was observed that the failure zone is relatively narrow and has a symmetrical pattern when rock bridges occupy a low percentage of the total shear surface. This may be due to the high stress interactions between the subsequent joints separated by a rock bridge. In contrast, when rock bridges are occupying sufficient area prohibiting the stress interactions to occur then the rupture of surface is more complex and turns into a shear zone. This zone was observed to be relatively thick with an unsymmetrical pattern. The shear strength of rock bridges is reduced by increasing the joint length as a result of increasing both the stress concentration at tip of the joints and the stress interaction between the joints.     

石灰岩细观力学特性的颗粒流模拟

岩体地区地质灾害的发生和发展取决于岩石细观组分的运动学行为。研究岩石运动学行为时通常将岩石作为整体研究对象较多,而直接以细观组分为对象的研究较少。以石灰岩为例,根据室内试验获得的岩石力学性质指标,使用基于非连续介质理论的颗粒流方法,将材料离散成刚性颗粒组成的模型,把颗粒细观变化与宏观力学特性联系起来,建立了石灰岩的细观结构模型,获得了颗粒接触力、颗粒接触模量、接触连接强度和连接刚度比等细观力学参数。由于文中直接以细观成分为研究对象、反映了岩石和岩体组成的本质特点,所得结论不仅对含裂隙岩石本构关系研究具有广阔的应用前景,而且对岩体工程性质和地质灾害机制研究也具有重要的理论意义。     

考虑颗粒抗转动的砂土双轴试验离散元模拟

颗粒间抗转动作用是影响砂土力学行为的重要因素。将蒋明镜等（2005年）提出的考虑颗粒抗转动作用的颗粒接触模型植入PFC2D中,对砂土双轴试验进行了模拟,研究颗粒抗转动作用对砂土力学行为的影响。研究显示：颗粒抗转动对砂土力学行为影响显著。抗转作用越强,砂土强度越高;抗转系数为0.4时,松砂亦出现软化及剪胀现象。抗转动作用对孔隙比及配位数变化趋势影响也很显著。