不同胶结砂土力学特性及胶结破坏的离散元模拟

采用离散元法研究不同胶结砂土的宏观力学特性以及微观胶结破坏,其中一种胶结砂土胶结特性由离散元商业软件PFC2D中的胶结接触模型(Contact bond model)进行控制,另一种则采用蒋明镜等提出的无厚度改进胶结接触模型(改进的蒋氏模型)控制.首先,将改进的蒋氏模型引入PFC2D;其次,使用PFC2D对不同胶结强度和不同围压下的上述两种胶结砂土进行双轴压缩试验模拟,对比分析了两种不同胶结砂土的模拟结果.结果表明蒋氏试样(胶结特性由改进的蒋氏模型控制)应变软化和体积剪胀显著,峰值内摩擦角基本随胶结强度增大而增大;而PFC试样(胶结特性由Contact bond model控制)在高胶结强度时应变软化和体积剪胀性比较显著,低胶结强度时表现为应变硬化和体积剪缩,峰值内摩擦角随胶结强度增加而减小.离散元模拟结果在一定程度上与已有室内实测结果相符,蒋氏试样更能反映胶结砂土的主要力学特性.由胶结破坏微观信息统计可知,蒋氏试样中胶结点受拉破坏率远大于受剪破坏率,而PFC试样中两者相当,且蒋氏试样中的总胶结点破坏速率峰值要大于PFC试样.通过对改进的蒋氏模型参数分析可知,低围压条件下,试样宏观力学特性与胶结破坏形式对切向胶结强度与法向胶结强度的比值较为敏感.     

深海能源土剪切带形成机理离散元分析

天然气水合物的分解开采过程将会劣化深海能源土的力学性能,从而引发一系列岩土工程问题。因此,要实现天然气水合物的安全开采,需要对能源土的强度和变形特性开展研究。结合深海能源土微观胶结模型,通过平面应变双轴试验的离散元模拟,研究了深海能源土剪切带形成机理以及剪切带内外的宏微观变量特征。结果表明：水合物胶结提升了深海能源土的强度,且使其呈现出明显的应变软化特性;剪切带在峰值应力后开始产生,伴随着胶结的大量破坏以及各宏微观变量的局部化;剪切带内外各宏微观变量差异明显,随着轴向应变的增加,土体微观结构也随之发生变化。     

结构性砂土粒间胶结效应的二维数值分析

将理想胶结颗粒接触力学特性的测试结果引入到离散元胶结接触模型中,对结构性砂土粒间胶结效应进行离散元数值模拟。首先,胶结颗粒被理想化为两铝棒在指定部位形成胶结,通过一系列加载试验（拉伸、压缩、压剪）获得胶结铝棒在不同应力路径下的接触力学响应。随后,将测试结果提炼总结后引入到自行开发的二维离散元程序 NS2D 中,用以模拟不同初始密度和胶结强度的结构性砂土等向压缩试验。最后,通过与人工胶结砂土的试验数据进行比较,对文中的数值模拟结果进行验证。研究表明：离散元数值模拟能够有效的捕捉结构性砂土的主要力学特性,即屈服强度和体积模量均随初始密度和胶结强度的变化而变化,且胶结试样的屈服强度与试样内部颗粒间胶结点破坏率密切相关。     

结构性砂土胶结厚度分布特性试验研究

颗粒间胶结物厚度对胶结颗粒接触力学特性有显著影响,进而影响结构性砂土离散元模拟结果。为更好地模拟天然砂土的结构性,需要明晰结构性砂土中颗粒间胶结厚度分布特性,并且能在离散元模拟中反映该分布特性。本文采用人工制备结构性砂土,开展电镜扫描试验,对结构性砂土中颗粒间胶结厚度分布特性进行了研究。研究结果表明:在结构性砂土中,随着水泥含量的增加,平均胶结厚度近似线性增加;但是水泥含量对胶结厚度的概率分布规律影响不大,水泥含量不同的结构性砂土,其颗粒间胶结厚度概率分布均为先线性增大,后指数型衰减。在试验研究结果的基础上,提出了离散元模拟结构性砂土胶结厚度的方法,可以较好地反映结构性砂土中颗粒间胶结厚度指数型衰减的分布特性。     

不同沉积方向各向异性结构性砂土离散元力学特性分析

为探究沉积方向对各向异性结构性砂土力学特性的影响,首先采用椭圆颗粒生成两种不同沉积方向(水平与竖直)的各向异性净砂样,其次,引入一个考虑胶结厚度影响的微观胶接触模型从而生成各向异性结构性砂土。最后,对两种各向异性结构性砂土试样进行双轴压缩试验,并将水平沉积试样试验结果与室内试验结果对比验证该模型的可行性。同时,将两种不同方向试样的试验结果进行对比以探究沉积方向的影响。结果表明:两种试样应力–应变关系均呈软化及剪胀现象,水平沉积试样峰值偏应力较竖直沉积试样大,而二者残余阶段偏应力无明显差别;水平沉积试样临界孔隙比较竖直沉积试样大;胶结破坏速率及胶结破坏率变化与宏观力学特性变化相对应,且变化规律基本相同;水平方向沉积试样的胶结接触主方向始终保持竖直而竖直沉积试样的胶结接触主方向始终保持水平;在水平沉积试样中始终为水平分布颗粒长轴分布主方向,而竖直沉积试样中颗粒长轴逐渐向各向同性分布靠近。     

各向异性砂土宏微观特性三维离散元分析

各向异性对砂土强度和变形特性有显著的影响,为了研究各向异性砂土的宏微观特性,基于三维离散元法,对7个不同沉积角的试样进行了一系列的三轴模拟试验。利用"Clump"命令生成近似椭球形状颗粒,并且采用三维抗转动模型来模拟颗粒间的抗转动能力。离散元模拟结果与已知室内试验结果吻合很好。结果表明:随着沉积角的增大,偏应力和轴向应变的关系逐渐由应变软化向应变硬化发展。沉积角较小的试样剪胀性更强并且容易到达临界状态,颗粒组构–应力联合不变量(表征颗粒长轴组构张量和应力张量的相对角度)的值接近于-1,且颗粒长轴组构各向异性先增大后减小;然而对于沉积角较大的试样,在轴向应变50%处,仍不能达到临界状态,并且联合不变量的值大于-1,颗粒长轴组构各向异性先减小后不断增大。对于法向接触组构,组构主轴方向迅速向应力主轴方向偏转,组构各向异性的演化规律与偏应力随轴向应变的演化规律相似。     

用于离散元分析的胶结材料三维微观接触模型

岩石、结构性土、含水合物深海沉积土(能源土)等岩土材料均可视为由颗粒骨架和粒间胶结物共同构成的胶结型材料。简单实用的微观粒间接触模型是采用离散元数值方法分析胶结体系宏微观力学特性的基础。针对有一定厚度的粒间胶结形式,将胶结物简化为存在于颗粒之间有一定厚度的短圆柱。通过对胶结物内部应力分布进行适当简化推导了三维粒间接触模型,包括法向、切向、弯曲和扭转四个方向的相互作用规律。最终所得模型结合了理论推导和既有室内胶结接触力学试验成果,考虑粒间胶结物的法向压碎过程,能描述切向、弯曲和扭转单独作用的破坏过程以及剪弯扭共同作用下的破坏准则。模型形式简单,易于引入离散元程序。     

结构性砂土力学特性三维离散元分析

利用离散元法对结构性砂土的三轴试验进行了三维数值模拟并对其宏观特性进行了分析。首先将考虑胶结尺寸(宽度和厚度)的三维胶结接触模型导入离散元软件PFC3D中,对结构性砂土数值试样进行三轴试验数值模拟;然后对比分析离散元模拟与室内试验结果;最后从宏观力学角度对试验结果进行了分析。离散元模拟结果表明:结构性砂土与无胶结松散砂土表现不同,其在低围压时表现出应变软化和体积剪胀特征,并随胶结含量的增加或围压的减少而愈发显著,在高围压时则呈应变硬化和体积剪缩现象;低平均应力时,随胶结含量的增加,试样峰值内摩擦角、黏聚力以及内摩擦角均增加,其中黏聚力增加较为明显,随着平均应力的增加,峰值强度包线逐渐趋向于无胶结土。     

深海能源土宏观力学性质离散元数值模拟分析

首先,引入笔者等^[16-17]所提出的微观胶结模型用以反映能源土颗粒之间水合物微观胶结接触力学特性;其次,采用C++语言将模型程序化,建立同商业软件PFC^2D的程序接口,将模型引入离散单元法中;然后,通过简化计算方法确定胶结宽度随水合物浓度的变化规律,进而确定水合物微观胶结参数;最后,根据所确定的胶结参数,针对不同水合物浓度试样进行能源土宏观力学特性离散元双轴试验模拟,并从应力应变、体积应变、水合物对能源土弹性模量的影响等方面与Masui等^[4]所进行的能源土室内三轴试验进行对比分析。结果表明：所选择胶结模型及微观胶结参数能有效反映深海能源土宏观力学规律;能源土峰值强度、弹性模量均随水合物浓度增加而增加,体积膨胀随水合物浓度的增加越来越显著。     

不排水循环荷载条件下胶结砂土宏微观力学性质离散元模拟研究

天然或人工胶结的存在能够提高砂土的抗液化能力,从宏微观尺度对其动力学性质进行研究具有重大意义.将已有的三维完整胶结接触模型引入到三维离散元程序中,对胶结砂土不排水循环三轴剪切试验进行三维离散元模拟,研究颗粒间胶结、循环应力比对离散元试样宏微观力学性质的影响.研究结果表明,胶结的存在能够抑制轴应变和孔压的发展,提高砂土的...     

基于微观破损机理的胶结砂土三维本构模型研究

在岩土破损力学和临界状态土力学框架内,遵循宏微观土力学的研究思路,建立了胶结砂土三维本构模型。定义与重塑砂土屈服面几何相似但尺寸扩大的胶结砂土屈服面;采用经三维离散元验证的Lade-Duncan强度准则作为临界状态强度面;基于胶结材料微观力学理论并结合三维离散元模拟结果,获得具有微观力学机制的胶结破损规律;将胶结破损规律引入到重塑砂土的硬化规律和流动法则,得到胶结砂土的硬化规律和流动法则。将该本构模型应用于人工制备胶结砂土室内常规三轴压缩试验和等平均应力真三轴试验的模拟,初步验证了该模型的适用性。     

用离散单元法分析单剪试验中粒状体的剪切机理(英文)

用铝棒堆积体作为二维粒状体的模型材料进行单剪试验 ,然后对其用离散单元法进行数值模拟 ,对于宏观应力应变关系 ,两者结果基本吻合。基于数据模拟结果 ,对粒状体的剪切机理 ,例如滑动面上粒子接点角的分布及变化规律 ,进行了微观分析。     

Strain localization analyses of idealized sands in biaxial tests by distinct element method

This paper presents a numerical investigation on the strain localization of an idealized sand in biaxial compression tests using the distinct element method (DEM). In addition to the dilatancy and material frictional angle, the principal stress field, and distributions of void ratio, particle velocity, and the averaged pure rotation rate (APR) in the DEM specimen are examined to illustrate the link between microscopic and macroscopic variables in the case of strain localization. The study shows that strain localization of the granular material in the tests proceeds with localizations of void ratio, strain and APR, and distortions of stress field and force chains. In addition, both thickness and inclination of the shear band change with the increasing of axial strain, with the former valued around 10–14 times of mean grain diameter and the later overall described by the Mohr-Coulomb theory.     

不同温压环境下深海能源土力学特性离散元分析

深海能源土宏观力学特性与所处温度与水压环境密切相关,明晰温压对能源土力学特性的影响对水合物的安全开采具有重要意义。首先,介绍深海能源土的温度-水压-力学微观胶结模型,用以描述能源土粒间的水合物胶结接触力学特性;其次,将该模型导入至离散元商业软件PFC2D中,开展不同温度与水压环境下的离散元双轴试验;最后,结合离散元双轴试验结果及同已有室内试验结果的对比分析,探讨温度与水压对深海能源土宏观力学特性的影响规律和微观作用机理。结果表明：引入胶结模型的离散元双轴试验可较好地描述深海能源土强度、变形等力学特性随温度与水压的变化关系;温度与水压影响深海能源土宏观力学特性的微观机理是颗粒间水合物胶结强度与刚度同温度与水压间的相关性;建议采用温压距离参数L（在无量纲化的温度-水压坐标平面内,水合物赋存温度与水压点至其相平衡线的最小距离）评价实际复杂温压场下的深海能源土宏观力学特性。     

平面应变状态下砂土剪切带分析

提出一种砂土剪切带的两相平衡模型,此模型可看作一种新颖的分岔解。把剪切带考虑作高应变相,剪切带外部其它区域考虑作低应变相。在跨越两相界面处,考虑了变形梯度和应力的不连续性,以及拖曳力和位移的连续性,并且强加了Maxwell关系。推导出砂土在平面应变载荷下的控制方程,并且此分析归结为寻找使控制方程得到唯一一组物理上可以接受的实数解的载荷的最小值。Maxwell应力、剪切带内部和外部的应力应变、剪切带倾角等均可确定,并与实验所观察到的结果相一致。     

含不同胶结充填物的碳酸盐岩层面剪切力学特性试验研究

 层面是层状岩体稳定性的主控结构面,层面的剪切力学行为及抗剪强度对建造于层状岩体中的地下工程具有重要影响。由于成岩时期物化作用的差异,同一工程场址常发育不同地质特征(胶结充填物、粗糙度等)的层面。为了深入了解不同类型层面的剪切力学特性,利用乌东德水电站地下厂房区域取得的胶结层面试样开展常法向应力及峰后降法向应力直剪试验。试验结果表明,不同地质特征的层面表现出不同的剪切力学行为,按有无明显峰值可将剪切变形–剪应力曲线分为2类,无明显峰值层面的粗糙度和表面附着物均有别于有明显峰值层面。法向变形与法向应力直接相关,但与层面类别之间的关系不明显。不同类型层面还表现出不同的剪胀行为及相应的剪胀角差异。胶结层面抗剪强度低于完整岩石但高于已分离层面,对层状岩体起弱化作用。按层面地质特征对其进行分类有助于进一步理解层状岩体的复杂力学行为,而选择不同层面力学模型可使计算更贴近实际工程岩体特征。