颗粒形状对粗粒土剪切变形影响的细观研究

颗粒形状是影响粗粒土密实度、力学与渗流等特性的主要因素之一。为了分析颗粒形状对粗粒土剪切特性的影响,采用离散元法生成4种不同形状的颗粒组,进行粗粒土直剪试验模拟与剪切宏细观响应研究,得出了颗粒形状对剪应力-剪位移、体应变-剪位移的影响,分析了粗粒土剪切应力、应变特性与剪胀特性。通过分析剪切带厚度、颗粒旋转量值、平均接触数、孔隙率及接触力系等宏细观参量的演化规律,研究颗粒形状在宏细观尺度上对粗粒土的影响。研究表明:异形颗粒间的咬合自锁作用大于纯圆颗粒,试样的抗剪强度有随形状系数的减小而增大的趋势。试样颗粒在外荷载作用下发生运动,应变主要表现在颗粒运动剧烈、剪胀幅度较大的剪切带内。颗粒形状系数F减小,试样的初始平均接触数增加,内摩擦角φ增大,剪切带内孔隙率增量越大,剪胀幅度越大。剪切过程中强力链聚集于剪切带内并起骨架作用,随着形状系数的减小,力链长度在0～5所占百分比呈增大趋势;剪切带内强力链的数目随着形状系数的减小而增加,峰值含量在30%～35%之间。     

反倾层状岩质边坡倾倒变形机理与影响因素的离散元模拟

为进一步研究层状反倾边坡的弯曲倾倒变形机制,以离心试验为原型,通过离散元数值模拟,研究了层状岩质反倾边坡的变形机理与影响因素。通过预置层内随机裂隙,实现了破裂面的形成和贯通。研究结果表明：模拟结果与试验吻合较好,边坡变形可分为起始蠕变、稳态变形和失稳破坏3个阶段;边坡破裂面在达到破坏荷载（G_f）后瞬间贯通,呈直线型,产状受岩层倾角控制,G_f值与坡角幂函数相关;反倾边坡的破坏需满足倾角和坡角启动条件,且变形破坏与岩层所受弯矩关系密切,当倾角为70°~80°、坡角大于60°时,最易破坏;典型破坏模式有倾倒-折断-块体式、倾倒-弯曲-折断式、倾倒-反折式3种,其受倾角、坡角组合控制;对材料参数的正交试验表明,各参数对G_f的敏感性从大到小依次为密度、层面内摩擦角、层厚、密度比、层面黏聚力,且G_f与层厚、层面内摩擦角及密度比具有良好的线性相关性;层面内摩擦角可影响破裂面产状,从而控制变形体规模,其他参数仅影响G_f的大小。     

An explicitly coupled hydro‐geomechanical model for simulating hydraulic fracturing in arbitrary discrete fracture networks

Modeling hydraulic fracturing in the presence of a natural fracture network is a challenging task, owing to the complex interactions between fluid, rock matrix, and rock interfaces, as well as the interactions between propagating fractures and existing natural interfaces. Understanding these complex interactions through numerical modeling is critical to the design of optimum stimulation strategies. In this paper, we present an explicitly integrated, fully coupled discrete‐finite element approach for the simulation of hydraulic fracturing in arbitrary fracture networks. The individual physical processes involved in hydraulic fracturing are identified and addressed as separate modules: a finite element approach for geomechanics in the rock matrix, a finite volume approach for resolving hydrodynamics, a geomechanical joint model for interfacial resolution, and an adaptive remeshing module. The model is verified against the Khristianovich–Geertsma–DeKlerk closed‐form solution for the propagation of a single hydraulic fracture and validated against laboratory testing results on the interaction between a propagating hydraulic fracture and an existing fracture. Preliminary results of simulating hydraulic fracturing in a natural fracture system consisting of multiple fractures are also presented. Copyright © 2012 John Wiley & Sons, Ltd.     

A bond contact model for methane hydrate‐bearing sediments with interparticle cementation

While methane hydrates (MHs) can be present in various forms in deep seabeds or permafrost regions, this paper deals with MH‐bearing sediments (MHBS) where the MH has formed bonds between sand grains. A bond model based on experimentally validated contact laws for cemented granules is introduced to describe the mechanical behavior of the MH bonds. The model parameters were derived from measured values of temperature, water pressure and MH density. Bond width and thickness adopted for each bond of the MHBS were selected based on the degree of MH saturation. The model was implemented into a 2D distinct element method code. A series of numerical biaxial standard compression tests were carried out for various degrees of MH saturation. A comparison with available experimental data shows that the model can effectively capture the essential features of the mechanical behavior of MHBS for a wide range of levels of hydrate saturation under drained and undrained conditions. In addition, the analyses presented here shed light on the following: (1) the relationship between level of cementation and debonding mechanisms taking place at the microscopic level and the observed macro‐mechanical behavior of MHBS and (2) the relationship between spatial distribution of bond breakages and contact force chains with the observed strength, dilatancy and deformability of the samples. Copyright © 2014 John Wiley & Sons, Ltd.     

Calculating the macrodispersion coefficient of the ensemble averaged solute transport equation in the discrete domain

Recognizing the spatial heterogeneity of hydraulic parameters, many researchers have studied the solute transport by both groundwater and channel flow in a stochastic framework. One of the methodologies used to up‐scale the stochastic solute transport equation, from a point‐location scale to a grid scale, is the cumulant expansion method combined with the calculus for the time‐ordered exponential and the calculus for the Lie operator. When the point‐location scale transport equation is scaled up to the grid scale, using the cumulant expansion method, a new dispersion coefficient emerges in the dispersive term of the solute transport equation in addition to the molecular dispersion coefficient. This velocity driven dispersion is called ‘macrodispersion’. The macrodispersion coefficient is the integral function of the time‐ordered covariance of the random velocity field. The integral is calculated over a Lagrangian trajectory of the flow. The Lagrangian trajectory depends on the following: (i) the spatial origin of the particle; (ii) the time when the macrodispersion is calculated; and (iii) the mean velocity field along the trajectory itself. The Lagrangian trajectory is a recursive function of time because the location of the particle along the trajectory at a particular time depends on the location of the particle at the previous time. This recursive functional form of the Lagrangian trajectory makes the calculation of the macrodispersion coefficient difficult. Especially for the unsteady, spatially non‐stationary, non‐uniform flow field, the macrodispersion coefficient is a highly complex expression and, so far, calculated using numerical methods in the discrete domains. Here, an analytical method was introduced to calculate the macrodispersion coefficient in the discrete domain for the unsteady and steady, spatially non‐stationary flow cases accurately and efficiently. This study can fill the gap between the theory of the ensemble averaged solute transport model and its numerical implementations. Copyright © 2011 John Wiley & Sons, Ltd.     

A discrete thermodynamic approach for anisotropic plastic–damage modeling of cohesive‐frictional geomaterials

A discrete plastic–damage model is developed for cohesive‐frictional geomaterials subjected to compression‐dominated stresses. Macroscopic plastic strains of material are physically generated by frictional sliding along weakness planes. The evolution of damage is related to the evolution of weakness planes physically in connection with the propagation of microcracks. A discrete approach is used to account for anisotropic plastic flow and damage evolution, by introducing two stress invariants and one plastic hardening variable for each family of sliding weakness planes. Plastic flow in each family is coupled with damage evolution. The proposed model is applied to typical geomaterials and comparisons between numerical predictions and experimental data are presented. Copyright © 2009 John Wiley & Sons, Ltd.     

A two‐mechanism elastoplastic model for shakedown of unbound granular materials and DEM simulations

An elastoplastic model has been developed for the finite elements modelling of repeated load triaxial tests. This model is based on the shakedown theory established by Zarka for metallic structures. To the previous works, which were based on the Drucker–Prager yield surface and the plastic potential of Von Mises, a compression cap has been added to each one. The model straightforwardly determines the purely elastic state or the elastic shakedown state or the plastic shakedown state and calculates the deviatoric and the volumetric plastic strains. The calibration of the elastoplastic model has been carried out with DEM simulations and an unbound granular material for roads under repeated load triaxial tests using finite element method. The calculations underline the capabilities of the model to take into account, with a unique formalism, the accumulation of the deviatoric and volumetric plastic strains along the loading cycles. Copyright © 2011 John Wiley & Sons, Ltd.     

基于离散裂缝网络模型裂缝储层的随机模拟——以王窑区块为例

在燕山期、喜山期构造活动作用下,王窑区块发育北东向和南北向共轭剪切裂缝。针对该裂缝进行了统计分析,并利用离散裂缝网络模型对裂缝进行了随机模拟。模拟结果表明:研究区共轭裂缝北东向较为显著,且以高角度裂缝为主。裂缝对渗透率影响较为明显,将研究区的渗透率提高到20×10-3μm2,孔隙度增加约2%。     

Stokes和Hotine积分离散求和的快速算法

提出了一种用于Stokes积分和Hotine积分直接离散求和的快速算法。该算法将积分核表达为计算点纬度、流动点纬度和两点间经度差的函数,充分利用核函数的对称性,相同纬度的所有计算点只需计算一组核函数,计算次数远少于普通离散求和。基于EGM2008地球重力位模型的模拟实验表明,快速算法的计算效率远高于普通算法,有效解决了离散求和计算速度太慢的数值问题,且保留了球面积分的特性,可取代一维FFT用于计算Stokes积分和Hotine积分。     

全球等面积四叉树离散格网建模与编码体系研究

在分析球冠、球带表面积公式的基础上,提出了一种等面积四叉树全球离散格网建模方法,阐述了格网划分机理,制定了格网编号体系,给出了由编号确定格网中心地理坐标及由地理坐标确定格网编号的公式,分析了其实际应用方向。