工程岩体失稳通常是由于非贯通节理的扩展聚合最终导致的剪切破坏，而能量驱动下的状态失稳是岩体破坏的根本原因。目前基于直接剪切试验的含非贯通节理岩石的力学行为研究中，大多集中在含有一条或两条节理的简单样品上，而忽略了天然岩体中的非贯通节理通常会以不同数量和几何形状的多种组合出现，与复杂的实际条件存在较大差异。因此，研究含多种组合形式的非贯通节理岩石的剪切力学性质，并分析其能量的演化过程，有助于更全面，更直观的理解含非贯通节理岩体的剪切破坏过程，在岩土工程和工程地质中具有重要意义。
中国科学院地质与地球物理研究所页岩气与地质工程院重点实验室台大平博士生在导师祁生文研究员的指导下，基于类岩石材料的室内直剪试验结果，使用颗粒流离散元PFC3D构建了一系列含有多种组合形式非贯通节理的岩石模型（BPM）并开展了直接剪切试验的数值模拟，从微裂纹扩展和能量演化的角度研究了节理排数和节理倾角对岩石应力应变特性和破坏模式的影响。
结果表明，导致含非贯通节理岩体破坏的主要是由于剪切-张拉型微裂纹（ST cracks）的扩展贯通，且对于“两列多排型”分布的节理岩石模型（Type B和Type C），当节理排数Nr大于1时，岩石的剪切破坏区域不随着Nr的增加而发生变化（图1，图2）。
图1 含水平多排节理的岩石剪切破坏时的微裂纹扩展形态与位移场分布
图2 含水平多排节理的岩石剪切破坏时的典型宏观裂纹形态
并且，节理会显著增加了模型在峰后的能量耗散，但对峰前弹性应变能与耗散能的在总能量中的占比几乎没有影响（图3）。在剪切过程中，摩擦消耗了大部分能量，而动能消耗只占总能量的不到1%。此外，峰值弹性应变能的变化趋势与峰值剪切位移保持一致。微裂纹的扩展和累积直接影响模型的能量耗散，并且临界微裂纹的累积数量与破坏时的临界耗散能量之间存在显著的线性关系，在相同倾角下，当节理倾向与剪切方向相反时，模型在峰值破坏时更容易积累较多的微裂纹，从而导致更多的能量耗散（图4）。
图3 剪切过程的能量占比变化
图4 剪切破坏时临界微裂纹数量与临界耗散能占比对比
该研究从微裂纹扩展和能量演化的角度定量分析了多种非贯通节理组合对岩石应力应变特性和破坏模式的影响，为理解含非贯通节理岩体的剪切破坏过程提供了新的视角。
研究成果发表于国际学术期刊Journal of Rock Mechanics and Geotechnical Engineering（台大平，祁生文*，郑博文，汪冲浪，郭松峰，罗光明. Shear mechanical properties and energy evolution of rock-like samples containing multiple combinations of non-persistent joints [J]. Journal of Rock Mechanics and Geotechnical Engineering, 2023. DOI: 10.1016/j.jrmge.2022.11.014.）。研究得到国家自然科学杰出青年基金项目（41825018）和第二次青藏高原综合科学考察研究（QZKK20190904）等项目资助。
