三维梁-颗粒模型在岩石材料破坏模拟中的应用

用三维梁-颗粒模型BPM3D(beam-particlemodelinthreedimensions)对岩石类非均质脆性材料的力学性质和破坏过程进行了数值模拟。梁-颗粒模型是在离散单元法基础上,结合有限单元法中的网格模型提出的用于模拟岩石类材料损伤破坏过程的数值模型。在模型中,材料在细观层次上被离散为颗粒单元集合体,相邻颗粒单元由有限单元法中的弹脆性梁单元联结。梁单元的力学性质均按韦伯(Weibull)分布随机赋值,以模拟岩石类材料力学参数的空间变异性。材料内部裂纹通过断开梁单元来模拟。通过自动生成的非均质材料模型对岩石类材料的破坏机理进行研究。岩石类非均质脆性材料在单轴压缩状态下破坏过程细观数值模拟结果显示,岩石材料宏观破坏是由于其内部细观裂纹产生、扩展、贯通的结果。通过数值模拟结果之间的对比分析,揭示出岩石试样宏观破坏模式随细观层次上韦伯分布参数的变化而不同。与实际矿柱破坏形态的对比分析表明了模型的适用性。根据数值模拟结果对岩石类非均质材料的破坏机理进行了探讨。     

弹丸侵彻无钢筋混凝土的数值模拟

采用二维梁-颗粒模型BPM2D（beam-particle model in two dimensions）模拟了刚性弹丸侵彻无钢筋混凝土的过程。离散元法（DEM）和有限元法（FEM）等数值计算方法各有其优势,同时也都存在不足之处：离散元法适于处理由连续介质向非连续介质转化的破坏问题,但对于连续体计算结果精度不高;有限元法适于预测材料破坏的区域,但难以直接用于计算脆性材料破坏过程,因此将两种方法结合可以形成一种较好的混合模型。梁-颗粒模型BPM2D是基于离散元法,结合有限元法开发的二维数值计算模型,采用3种类型梁单元形成混凝土数值试样,每种类型梁单元的力学性质均按韦伯（Weibull）分布随机赋值,以模拟混凝土细观结构的非均匀性,同时梁单元的强度随应变率不同而变化。利用此模型分析了弹丸侵彻下混凝土的破坏过程,并给出侵彻过程弹丸减速度-时间历程曲线。比较计算结果与试验数据表明梁-颗粒模型可有效应用于计算和模拟脆性材料动态破坏问题。     

考虑微裂纹随机分布的堆石颗粒准静态强度统计模型

材料内部裂纹尺寸、形状、转向和空间位置的随机分布对其强度有直接的影响。为了研究地震准静态荷载作用下堆石颗粒强度与尺寸和应变率的统计关系,在最弱链理论的框架下假设颗粒内部裂纹的尺寸和空间位置均服从幂律分布,建立了由颗粒累积破坏概率和体积构成的复合参数与强度的关系。从微观角度应变率效应提高了内部裂纹的强度,降低了单位体积的失效概率,同时使裂纹的空间位置分布更稀疏,减弱了颗粒强度的尺寸效应。不同加载速率的颗粒破碎试验结果表明：随着加载速率的增大,存在逐渐减小的空间分布幂指数使颗粒复合参数汇集在由最弱链统计模型决定的主曲线上。     

一种有限元转化为颗粒离散元的方法及其应用研究

为了充分发挥有限元与颗粒离散元各自的优势,提出了一种由有限元转化为颗粒流的方法。数值模型首先用较粗的有限元网格进行离散,并在单元上引入连续介质本构模型。力学计算开始后,实时跟踪各单元的应力状态。一旦某单元的应力满足Mohr-Coulomb准则或最大拉应力准则,删除该单元,同时创建具有一定数目、随机分布且微嵌入的颗粒簇。其后,该单元所在区域的非连续变形及失稳断裂由颗粒簇演化获得。各颗粒的质量、材料参数、速度、位移、接触力等信息根据插值从有限元单元中继承。为了实现有限元与颗粒流接触面的耦合计算,引入了点-棱（二维）及点-面（三维）接触模型,通过法向及切向弹簧实现接触力的计算。颗粒球与有限元板的碰撞分析、单轴压缩、岩石切割等案例展示了上述方法的精确性及合理性。     

单轴压缩试验下粗糙离散节理网络模型建立及力学特性

针对岩体中随机结构面存在粗糙特性的特点,开展了粗糙离散节理网络(roughness discrete fractures network,RDFN)模型研究。采用数字图像处理技术,开发了生成RDFN模型的算法,并基于Matlab方法建立了与颗粒流离散元PFC方法的接口。采用颗粒流PFC方法研究粗糙度对RDFN模型单轴压缩下的弹性、强度和破坏模式的影响规律,并与直线型DFN模型进行对比。结果表明,RDFN模型的弹性模量低于DFN模型,而峰值强度较高,残余强度较低;DFN模型中微裂纹沿节理方向贯通,局部存在基岩内部裂纹扩展;RDFN模型由破坏模式相对复杂,达到峰值强度后基岩发生显著破坏。     

颗粒形态对离散介质剪切强度的影响

自然条件下,颗粒介质大多以非规则单元形态存在。非规则几何形态对颗粒介质的宏观力学性能有很大影响。针对颗粒单元的不同几何形态,采用团颗粒单元对离散介质的直剪试验过程进行了离散元数值计算,详细地讨论了颗粒形态对离散介质剪切强度的影响。该非规则颗粒由不同形态、不同数目、镶嵌尺寸、组合方位和颗粒大小的球形颗粒进行随机构造,其在局部与整体坐标之间的转动、力矩和方位关系通过4元素方法进行确定,基本球体颗粒之间的作用力采用具有Mohr-Coulomb摩擦定侓的Hertz-Mindlin非线性接触模型,并考虑了非线性法向粘滞力的影响。通过构造7种具有相同的质量概率分布的不同形态的团颗粒,在不同法向应力下,对团颗粒的直剪试验进行了离散元模拟,分析了不同形态颗粒的剪切强度。通过对不同形态颗粒介质剪切强度的数值分析,进一步揭示了非规则颗粒间的咬合互锁效应,为分析非规则颗粒的宏观动力特性提供了依据。     

模拟岩石破裂过程的块体单元离散弹簧模型

在变形体离散元的基础上建立块体单元离散弹簧模型,并应用于岩石破裂过程的数值模拟研究。该模型以连续介质力学理论为基础,将块体单元离散为具有明确物理意义的弹簧系统,通过对弹簧系统的能量泛函求变分获得各弹簧的刚度系数,进而可以直接利用弹簧刚度求解单元的变形和应力,提高计算效率。以重力作用下的岩质边坡计算为例,通过与传统的有限元进行对比,验证该模型弹性计算结果的正确性。在该基础上,引入Mohr-Coulomb与最大拉应力的复合破坏准则,判断单元的破坏状态及破裂方向。当单元的内部破坏面确定后,则通过块体切割的方式实现单元破坏,并建立单元边界和单元内部的双重破裂机制,实现块体由连续到非连续的破裂过程,进而显示的模拟裂纹的形成和扩展。最后,以巴西圆盘劈裂、单轴压缩破裂以及三点弯曲梁等典型算例验证该方法,结果表明该方法可以较好地模拟拉伸、压剪等应力状态下裂纹的形成和扩展,从而可模拟岩石介质由连续到非连续的破裂过程。     

MICP对钙质砂单颗粒的破碎行为影响研究

作为一种特殊的岩土介质材料,钙质砂具有在低压下易破碎的性质。微生物诱导方解石沉淀（MICP）技术得到了广泛的关注和认可,可用来改善钙质砂的破碎特性。文章从室内试验和离散元模拟两个角度分别对钙质砂颗粒MICP固化前后进行单颗粒压碎试验,通过Weibull分布和SEM扫描等探究了MICP对钙质砂颗粒破碎行为的影响。结果表明：离散元模拟得到的生存概率曲线及Weibull模量m值与试验结果均吻合较好,验证了该数值模型的有效性。与室内试验相比,数值模拟可以精确地反映颗粒的裂纹分布及破碎过程,且可以研究同一颗粒MICP固化前后的情形,弥补了室内实验的不足,但其取决于模型参数的选取;经过MICP固化后的钙质砂颗粒表面有明显的方解石结晶生成,颗粒表面及内孔隙分别得到一定程度的包裹和填充,导致颗粒破碎强度有明显的增强且离散性大大降低,破碎模式由“多峰型”向“单峰型”转变,局部裂纹减少,多以表面磨损和直接产生贯穿裂纹为主。     

自由边界形状与近边界裂纹相互作用模型研究

大量的研究表明, 岩爆发生条件与岩石的力学性能有关。岩体中微裂纹的形成、扩展和汇合对岩体的力学性能产生显著影响,可以导致材料的逐渐劣化直至破坏。对于单轴压缩的岩石试件,其破坏型式主要是微裂纹的摩擦滑移、自相似平面扩展和弯折扩展。基于断裂力学机制,对单轴压缩荷载下3种不同边界形状,即直边、凸形、凹形情况下近边界原生斜裂纹的扩展进行了理论研究,并运用FLAC数值模拟技术,建立了基于Interface界面单元的裂纹扩展数值模型。近直边界及凸形边界情况的裂纹扩展使用梁渐近模型来解析薄片岩层的屈曲断裂问题,对于近凹形边界情况给予简化,提出了裂纹－圆－椭圆相互作用的修正模型并进行理论分析。研究结果表明,基于Interface界面单元的裂纹扩展的数值模型正确,解析解与数值解吻合较好,相应的梁渐近模型及裂纹-圆-椭圆模型可以合理地分析应力调整过程中裂纹的扩展和裂纹与不同洞室表面形状的相互作用问题,为岩爆工程及理论研究提出了新的理论依据。     

冻结黏土三轴试验微观变形机理的离散元分析

基于三维颗粒离散单元法,赋予颗粒相应的细观参数,并采用黏结发生在接触颗粒间有限范围内的模型来考虑冻土颗粒中冰的胶结作用,建立了冻结黏土三维离散元数值模型.在相同围压、不同温度和相同温度、不同围压下对冻结黏土的室内三轴试验进行数值模拟,对比了数值试验与室内测试的应力-应变曲线,两者吻合较好.数值模拟结果表明:围压增大会使得接触黏结逐渐失效,在剪切带中胶结冰的破坏区域将增大,而温度的降低则会产生相反结果,这些微观变化都将对冻结黏土的宏观力学变形产生较大影响,同时,细观参数对温度的依赖性也很明显.冻结黏土三轴试验微观变形离散元模拟思路及方法可为今后运用离散单元法研究冻土力学行为提供一定的参考.     

压实黏土的脆性断裂模型及有限元算法

利用弥散裂缝理论,提出了压实黏土拉伸状态下的脆性断裂模型.当压实黏土达到其极限抗拉强度后,通过建立单元的各向异性刚度矩阵,将土体裂缝弥散于实体单元,构造了平面应变条件下考虑压实黏土脆性开裂的有限元计算模式.通过对某压实黏土单轴拉伸试验成果的模拟计算,验证了构建的脆性断裂模型和有限元算法对土体拉伸破坏特性和裂缝发展过程的适用性.本文还进行了模拟软弱面水压“楔劈效应”的简单数值试验,表明压实黏土脆性开裂模型和算法可较好地模拟裂缝扩展行为.     

Meshless numerical modeling of brittle–viscous deformation: first results on boudinage and hydrofracturing using a coupling of discrete element method (DEM) and smoothed particle hydrodynamics (SPH)

We have developed a new approach for the numerical modeling of deformation processes combining brittle fracture and viscous flow. The new approach is based on the combination of two meshless particle-based methods: the discrete element method (DEM) for the brittle part of the model and smooth particle hydrodynamics (SPH) for the viscous part. Both methods are well established in their respective application domains. The two methods are coupled at the particle scale, with two different coupling mechanisms explored: one is where DEM particles act as virtual SPH particles and one where SPH particles are treated like DEM particles when interacting with other DEM particles. The suitability of the combined approach is demonstrated by applying it to two geological processes, boudinage, and hydrofracturing, which involve the coupled deformation of a brittle solid and a viscous fluid. Initial results for those applications show that the new approach has strong potential for the numerical modeling of coupled brittle–viscous deformation processes.     

Numerical simulation of particle plugging in hydraulic fracture by element partition method

Hydraulic fracturing (HF) treatment often involves particle migration and is applied for propping or plugging fractures. Particle migration behaviors, e.g., bridging, packing, and plugging, significantly affect the HF process. Hence, it is crucial to effectively simulate particle migration. In this study, a new numerical approach is developed based on a coupled element partition method (EPM). The EPM is used to model natural and hydraulic fractures, in which a fracture is allowed to propagate across an element, thereby avoiding remeshing in fracture simulations. To characterize the water flow process in a fracture, a fully hydromechanical coupled equation is adopted in the EPM. To model particle transportation in fractures with water flow, each particle is treated as a discrete element. The particles move in the fracture as a result of being dragged by fluid. Their movement, contact, and packing behaviors are simulated using the discrete element method. To reflect the plugging effect, an equivalent aperture approach is proposed. Using this method, the particle migration and its effect on water flow are well simulated. The simulation results show that this method can effectively reproduce particle bridging, plugging, and unblocking in a hydraulic fracture. Furthermore, it is demonstrated that particle plugging significantly affects water flow in a fracture and hence the propagation of hydraulic fracture. This method provides a simple and feasible approach for the simulation of particle migration in a hydraulic fracture.     

A generalized rigid particle contact model for fracture analysis

Rigid particle models taking directly into consideration the physical mechanisms and the influence of the material meso‐structure have recently been developed for fracture studies of quasi‐brittle material such as concrete. The formulation of a generalized contact model for rigid particle simulations is presented in which the contact discretization is a model parameter. The contact model performance for different discretizations is evaluated for uniaxial tensile tests, for uniaxial compression tests and for a notched beam in mode I. Copyright © 2005 John Wiley & Sons, Ltd.     

A 2D mixed fracture–pore seepage model and hydromechanical coupling for fractured porous media

A novel two-dimensional mixed fracture–pore seepage model for fluid flow in fractured porous media is presented based on the computational framework of finite-discrete element method (FDEM). The model consists of a porous seepage model in triangular elements bonded by unbroken joint elements, as well as a fracture seepage model in broken joint elements. The principle for determining the fluid exchange coefficient of the unbroken joint element is provided to ensure numerical accuracy and efficiency. The mixed fracture–pore seepage model provides a simple but effective tool for solving fluid flow in fractured porous media. In this paper, examples of 1D and 2D seepage flow in porous media and porous media with a single fracture or multiple fractures are studied. The simulation results of the model match well with theoretical solutions or results obtained by commercial software, which verifies the correctness of the mixed fracture–pore seepage model. Furthermore, combining FDEM mechanical calculation and the mixed fracture–pore seepage model, a coupled hydromechanical model is built to simulate fluid-driven dynamic propagation of cracks in the porous media, as well as its influence on pore seepage and fracture seepage.

     

基于改进颗粒流声发射片的地下厂房洞室围岩局部损伤细观机制研究

通过改进的PFC-FLAC离散−连续耦合计算方法,提出了一种基于FLAC连续模型耦合区域节点速率双线性插值的PFC颗粒流声发射片模拟方法以及基于声发射强度的围岩破坏区深度细观角度判别方法,并将其应用于地下厂房洞室开挖面、岩锚吊车梁与岩壁接触面等局部围岩损伤细观机制与破化特性研究。研究结果表明：浅层围岩接触力链逐渐稀疏,大量微裂纹发育并汇集成宏观裂隙,最终出现浅层破裂区,呈拉裂破坏;随着距离开挖面的增大,围岩从开挖前期受力较大产生破坏、到开挖后期只有少量微裂纹产生,对应了围岩回弹卸荷区深度。在岩壁围岩劣化和吊车梁超载的情况下,吊车梁与岩壁竖直接触面和倾斜接触面分别发育大量拉裂纹和剪裂纹,宏观表现为拉裂破坏与滑移破坏。上述分析结果与三维有限元方法相一致,并弥补了后者围岩破坏显示方法单一、难以描述围岩损伤程度变化的缺点,为研究地下厂房洞室大变形与应力集中部位的宏细观特性与损伤机制提供参考。     

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.     

大柳树坝址工程断错效应分析

坐落于中卫-同心活动断裂带中间部位的黄河黑山峡大柳树坝址,距F201发震断层垂直距离不足1.5km。统计表明,在青藏高原及其周缘地区,当发生Ⅶ级以上破坏性地震时,距发震断裂两侧3～5km的范围内都有可能出现断层的分支破裂、次生破裂等错断效应。大柳树坝址处于Ⅶ级以上地震时的分支破裂发生带内,因此,存在工程错断效应问题。三维有限元数值模拟结果表明,断层F201发生错动时,坝址地段断层F93、F39、F40相对于F201的错动比率处分别达到14.38%、12.00%和9.84%以上。     

基于离散裂隙网络模型的核素粒子迁移数值模拟研究

高放废物地质处置过程中涉及的核素在围岩裂隙地下水中的迁移问题已引起广泛关注,数值模拟是研究核素粒子运移的重要方法。目前裂隙介质中渗流模型主要是等效连续介质模型、双重介质模型和离散裂隙网络模型。对于岩体尺度裂隙地下水的流动,离散裂隙网络模型能充分表现裂隙介质的各向异性、不连续性等特征。因此,针对裂隙介质准确概化及核素迁移模拟等难点,文章结合Monte Carlo随机生成裂隙方法、裂隙渗流有限元算法和高放射性核素衰变方程等方法,依据花岗岩深钻孔裂隙统计数据,采用离散裂隙网络模型对内蒙古阿拉善高放废物地质处置预选区展开了核素粒子迁移数值模拟研究,并讨论了实例预测分析结果。结果显示:针对设定的地质模型,核素粒子从中心运移到边界的迁移路径长度平均为1293.35 m,粒子运移到边界耗费的时间平均为1.70E+11 d。