含2条断续裂隙试件直接拉伸下的颗粒流模拟

采用颗粒流程序PFC构建了试件直接拉伸模型,模拟了含2条断续裂隙试件的直接拉伸试验,通过完整试件的模拟结果与室内试验结果的对比,验证了细观参数的准确性和可靠性,同时分析了裂隙倾角和岩桥倾角对抗拉强度和裂纹扩展的影响。研究表明:在直接拉伸情况下,含2条断续裂隙试件以萌生次生倾斜反翼裂纹为主,少量次生共面裂纹和翼裂纹;岩桥倾角和裂隙倾角均较小时,次生共面裂纹与原始裂隙之间发生搭接、贯通,随着岩桥倾角和裂隙倾角的加大,次生倾斜反翼裂纹扩展使试件断裂;裂纹扩展的方向与最大拉应力的方向基本保持垂直;岩桥倾角对于试件的抗拉强度和起裂应力影响不大;裂隙倾角对于试件抗拉强度和起裂应力影响明显,裂隙倾角越大,试件抗拉强度和起裂应力越大。     

考虑脆性损伤和渗流的圆形水工隧洞解析解

针对脆性岩石的损伤特性引入脆性损伤模型,导出三维各向同性损伤情况下的损伤变量表达式。在考虑损伤、不同工况主应力对应关系和渗流作用影响下,推导出圆形水工隧洞洞周围岩的弹塑性应力和塑性半径的解析解。对岩石的脆性特性的研究表明,考虑脆性损伤时的洞周围岩塑性半径比不考虑时有一定程度提升;围岩脆性特征越明显,塑性半径越大。无论洞周围岩应力状态如何,考虑渗流影响后,洞周围岩塑性区范围均会有所扩展。在内水压从0逐渐增加过程中,塑性半径随内水压增加而逐渐减小,直至塑性区消失,在超出临界内压后,塑性区重新出现,塑性半径逐渐增大。     

含２条断续裂隙试件直接拉伸下的颗粒流模拟

采用颗粒流程序 ＰＦＣ构建了试件直接拉伸模型，模拟了含２条断续裂隙试件的直接拉伸试验，通过完整试件的模拟结果与室内试验结果的对比，验证了细观参数的准确性和可靠性，同时分析了裂隙倾角和岩桥倾角对抗拉强度和裂纹扩展的影响。研究表明:在直接拉伸情况下，含２条断续裂隙试件以萌生次生倾斜反翼裂纹为主，少量次生共面裂纹和翼裂纹；岩桥倾角和裂隙倾角均较小时，次生共面裂纹与原始裂隙之间发生搭接、贯通，随着岩桥倾角和裂隙倾角的加大，次生倾斜反翼裂纹扩展使试件 断裂；裂纹扩展的方向与最大拉应力的方向基本保持垂直；岩桥倾角对于试件的抗拉强度和起裂应力影响不大；裂隙倾角对于试件抗拉强度和起裂应力影响明显，裂隙倾角越大，试件抗拉强度和起裂应力越大。      

离散裂隙网络对岩石力学性质和声发射特性影响的颗粒流分析

基于颗粒流软件平台,采用图像处理技术和Monte Carlo法,生成离散裂隙网络。结合等效晶质模型进行花岗岩细观结构的数值建模,从细观角度阐述了离散裂隙网络对岩石强度与变形性质、声发射特性的影响。主要结论如下：（1）裂隙平均长度和密度影响着岩石的力学性质,岩石的单轴压缩强度、弹性模量随裂隙密度的增加而减小。（2）单轴压缩下岩石破坏时均是晶间-拉裂纹和晶内-拉裂纹为主。（3）离散裂隙网络与岩石的声发射事件密切相关。随着裂隙密度的增加,声发射b值和分形维数D值均呈下降趋势。当b值和D值急剧下降时,预示着岩石失稳破坏。（4）离散裂隙网络对声发射事件的破裂强度有一定的影响。在离散裂隙网络附近易产生破裂强度大的声发射事件。综上所述,基于离散裂隙网络和等效晶质模型联合建模获得的细观角度上岩石破裂机制和声发射现象等结论,可有效弥补现有室内试验中研究手段的不足,并对野外具有更为复杂的裂隙网络的岩体稳定性评价与工程施工提供有力支撑。     

水工隧洞围岩裂纹扩展的临界水压解析解

通过把断裂力学扩展判据引入水工隧洞解析解中,使宏观的水工隧洞受力分析与细观的围岩裂纹扩展研究相结合;将隧洞围岩裂纹分为拉剪型裂纹和压剪型裂纹,并针对不同类型裂纹的开展机制,分别推导了水工隧洞围岩拉剪型和压剪型裂纹扩展时的临界内水压力表达式。在此基础上,分析了裂纹的走向、长度、埋深等参数对临界水压的影响规律。算例显示,在高内水压力作用下,或者由于地应力释放的影响,裂纹既可能发生拉剪型扩展,也可能发生压剪型扩展,但后者发生的几率更大。在内水压力控制工况下,当裂纹走向垂直于洞壁法线方向时最不容易发生扩展;裂纹长度越长,埋深越浅,围岩孔隙水的作用面积越大,裂纹发生扩展的可能性越大。     

含顺层断续节理岩质边坡地震作用下的破坏模式与动力响应研究

基于二维颗粒流软件PFC2D的人工合成岩体技术（SRM）,研究了岩桥倾角和节理间距不同组合形式的含顺层断续节理岩质边坡在地震作用下的破坏模式与动力响应规律。研究结果显示：在地震动力作用下,含单潜在滑动面的顺层断续节理岩质边坡呈现出滑移-倾倒的混合破坏特征,含多潜在滑动面的顺层断续节理岩质边坡则主要发生倾倒破坏;由顺层断续节理以及岩桥交替连接所组成的潜在滑动面是控制边坡动力稳定性的关键因素。在地震动力作用下,最靠近坡脚的岩桥段首先萌生翼裂纹,使得拉应力得到释放,随后各节理相继萌生裂纹并扩展、贯通,最终导致坡体发生阶梯状整体失稳。裂纹扩展受顺层断续节理控制,萌生裂纹中以张拉裂纹为主,且裂纹数量与输入地震波的加速度曲线具有同步性。另一方面,节理面的存在对边坡动力响应产生明显影响,沿坡表以及沿水平方向上的峰值速度、峰值位移随着岩桥倾角的增大、节理间距的减小而增大,同时节理间距和岩桥倾角对于峰值加速度（PGA）放大系数的影响范围主要集中在坡表、坡肩;沿竖直方向上,峰值位移随着岩桥倾角、节理间距的增大而减小,PGA放大系数曲线随高程变化总体呈现U型分布特征。     

基于MeshPy的3DGIS与三维有限元数值计算无缝耦合方法

3DGIS与有限元数值计算耦合方法及其系统开发,是目前GIS与岩土工程交叉学科的重要研究内容。通过对MeshPy、GRASS GIS、vtkPython以及Scipy等多个开源库的探索与编程实践,以Python为“黏合剂”,构建了GIS与三维有限元数值计算无缝耦合的平台框架。在满足几何形状、Delaunay特性、角度、面积以及体积等约束条件下,通过MeshPy包含的Triangle与TetGen库获得能用于有限元数值计算的三角形与四面体网格,并阐述了限定Delaunay三角化（CDT）节点插入算法与四面体网格质量控制方法。以排土场为例,针对CDT网格建立的地表模型,结合vtkPython封装的模块与算法,提出了根据节点分布位置综合采用线性与三角形重心坐标的插值方法,有效消除了其他插值手段导致的部分地表高程畸变的现象。借助Scipy提供的矩阵运算模块,开发了邓肯-张E-? 模型,给出了采用中点增量法进行非线性求解的算法,并通过对室内大型三轴剪切试验的数值模拟验证了程序的正确性。最后,对准东露天煤矿北部排土场堆排过程开展三维有限元数值计算,结合GRASS GIS对排弃物料高度与厚度空间分布规律的分析结果,对排土场边坡位移变化规律进行初步研究,结果表明,排土场垂直位移与排弃料高度及厚度相关,而最大水平位移则出现在各台阶坡体边缘位置,与实际情况一致。