《Computational Particle Mechanics》：Hierarchical Load-Bearing Organization of Heterogeneous Granite: Strong Force Chains, Mineral-Scale Partitioning, and Fracture Localization
编辑推荐：
介观非均质性通过影响力传递、损伤累积和断裂局部化，在结晶岩的脆性失稳中发挥关键作用。然而，接触尺度非均质性、强力链演化与矩张量声发射（AE）特征之间的三维关系仍未被充分认识。本研究采用PFC3D-GBM对多矿物花岗岩的单轴压缩进行模拟，并显式区分晶内接触与晶间
介观非均质性通过影响力传递、损伤累积和断裂局部化，在结晶岩的脆性失稳中发挥关键作用。然而，接触尺度非均质性、强力链演化与矩张量声发射（AE）特征之间的三维关系仍未被充分认识。本研究采用PFC3D-GBM对多矿物花岗岩的单轴压缩进行模拟，并显式区分晶内接触与晶间接触。利用均质度指数m = 2、5、10和20的Weibull分布接触属性以及均质条件进行赋值，并通过多级力链分析、时空微裂纹聚类和矩张量计算进行分析。弹性模量对非均质性几乎不敏感，而峰值强度从59.81 MPa增至153.81 MPa，峰值应变从0.296%增至0.696%。在低均质度条件下，较多弱接触与早期弥散开裂和力传递路径中断相关。随着均质度增加，更连续的晶内强力链骨架得以发育，能够实现更大的峰前荷载累积并在应力峰值附近快速退化。在峰值阶段，绝对强力链荷载顺序为斜长石 > 晶间结构 > 石英 > 碱性长石，而单个强力链接触的平均荷载顺序为石英 > 斜长石 > 碱性长石 > 晶间结构。声发射活动从稀疏、空间分散、低震级事件向更密集的事件簇转变，伴随更强的能量释放和沿优势斜向断裂带的增强局部化。总体而言，结果揭示了接触尺度非均质性、强力链组织与坍塌、裂纹局部化、AE震源演化与宏观脆性失稳之间一致的多尺度对应关系。
研究背景：深部地下工程（超深隧道、地热开发、CO2封存和核废料处置等）不断向高地应力和强扰动环境延伸，脆性岩体的突然失稳与能量释放构成重大安全和经济风险。此类破坏表现为局部裂纹起裂、快速断裂局部化和峰后应力骤降，这些宏观响应不仅取决于岩石的整体平均强度，还取决于其介观结构中局部承载路径的形成、演化和坍塌。花岗岩等结晶岩是由矿物颗粒和晶界组成的多矿物集合体，矿物弹性性质与强度差异、晶界弱化以及矿物相的空间排列共同影响裂纹起裂、扩展和宏观脆性响应。目前存在的问题包括：许多非均质建模方法仍基于等效或二维表征，限制了三维颗粒结构与晶界效应的同时刻画；接触属性非均质性未能在矿物类型和晶内/晶间接触类别间系统设定，限制了对弱单元丰度和承载路径稳定性的定量评价；尽管力链分析已用于研究颗粒体系中的荷载传递，但结晶岩中矿物内部和晶界网络中多级强力链网络的识别与对比演化仍研究不足；许多数值研究将每次键断裂视为独立AE事件或主要以裂纹计数表征损伤，但单次键断裂不一定对应实验检测到的AE事件。本研究旨在探究Weibull控制的接触属性分布变化如何重组跨矿物内部和晶界的强力链骨架，以及该骨架的后续退化如何与断裂局部化相关联。
研究方法：研究人员在固定的PFC3D-GBM框架内引入Weibull均质度指数m作为统计控制变量，在试样、晶内/晶间和矿物尺度上量化强力链骨架的演化，重点关注其形成、强化和坍塌；进一步区分不同结构域的整体荷载贡献与单个强力链接触的平均承载强度；最后将承载结构的重组与从弥散开裂向局部化断裂的转变以及矩张量导出的AE特征演化相关联。
主要技术方法：采用PFC3D-GBM构建圆柱形花岗岩模型（直径50 mm、高100 mm），矿物颗粒尺寸3–5 mm，含26%石英、50%斜长石和24%碱性长石；采用线性平行粘结模型，基于Li et al.的力学参数进行标定，数值与实验弹性模量和强度偏差在5%以内。通过FISH函数遍历接触并按Weibull分布重新赋值接触属性（m = 2、5、10、20及均质条件），采用蒙特卡洛逆变换抽样生成Weibull随机数。以峰值应力状态的平均接触力为固定阈值识别强力链，在试样整体、晶内/晶间和矿物三个层级进行分析。采用时空聚类方法将微裂纹分组为等效AE事件，计算矩张量、标量矩和矩震级，聚类参数为ΔCAE = 10个计算循环和RAE = 2.5d*。每个Weibull条件生成三个独立随机种子实现以评估随机稳健性。
研究结果：
3.1 宏观力学响应：应力–应变曲线在非线性段和峰后段表现出明显的均质度依赖性。弹性模量E在22.8–23.2 GPa间几乎不变，表明初始宏观刚度对非均质性影响较弱。峰值强度从m = 2的59.81 MPa增至m = 5、10、20的103.20、122.70、134.53 MPa，均质条件下达153.81 MPa；峰值应变从0.296%增至0.696%。所有试样峰后应力快速下降，均质度越高应力降越陡峭。裂纹演化在所有条件下均遵循起裂、扩展、局部化和峰后快速增长序列，但m = 2时低荷载水平（0.1–0.3 UCS）即出现大量拉（红）剪（蓝）微裂纹且分布广泛，高均质度则峰前损伤有限、近峰和峰后裂纹快速合并形成集中破坏带。拉伸裂纹在所有条件下占主导。
3.2 力链演化与阶段特征：一般力链网络在整个试样中密集分布，随m增大高力链段更显著。强力链网络在所有条件下均遵循形成—渗流—坍塌的阶段演化：40%峰值荷载时强力链稀疏孤立，80%时更密集互联，峰值时形成近跨试样承载骨架，峰后快速解体为局部残余链段。高m时残余强力链集中于局部化条带，低m时更破碎分散。
3.3 晶内与晶间结构荷载分配：晶内强力链网络形成更早、接触更多、覆盖范围更大，是破坏前的主要承载结构；晶间网络相对稀疏、承载份额有限但辅助局部力传递。峰值阶段晶内强力链累积力从约2.0 MN（m = 2）增至约10.0 MN（均质），晶间从约0.7增至4.2 MN，总强力链荷载从2.77增至14.3 MN。晶内与晶间裂纹数随m呈非单调变化，从m = 2至10均减少，晶内/晶间裂纹比从0.87498降至0.49814；m = 20和均质条件下再次增加（均质时约4055条晶内裂纹和6696条晶间裂纹，比值0.60559），反映从低均质度弥散峰前损伤向高均质度集中近峰断裂的转变，晶间开裂始终占主导。
3.4 矿物尺度选择性承载与失稳机制：峰值阶段晶内强力链网络呈现明显矿物依赖差异。m = 2时三种矿物中强力链均稀疏，斜长石空间覆盖最广，其次石英，碱性长石最少；随m增大三种矿物网络均更密集互联，石英和斜长石尤为明显，碱性长石始终较稀疏。峰值阶段各域强力链荷载随非均质性降低而增大：斜长石从约1.06增至5.09 MN，晶间从0.71增至4.21 MN，石英从0.64增至2.92 MN，碱性长石从0.36增至2.04 MN，荷载层级始终保持斜长石 > 晶间网络 > 石英 > 碱性长石。强力链接触数从m = 2至5显著增加，m = 5后趋于平台；晶间和斜长石域接触数最多（约(3.6–5.0)×10?），石英次之（约(2.0–2.7)×10?），碱性长石最少（约(1.6–2.1)×10?）。单个强力链接触平均荷载随均质度稳步上升，顺序为石英 > 斜长石 > 碱性长石 > 晶间结构，与总荷载排序不同，表明斜长石的总荷载优势部分源于其强力链接触数量更多。各域裂纹数随m非单调变化，排序始终为晶间 > 斜长石 > 碱性长石 > 石英；拉伸裂纹占主导，剪切裂纹对均质度更敏感，m = 5后晶内剪切裂纹急剧减少趋近于零。随机实现评估表明峰值强度、弹性模量、强力链荷载和裂纹数的变异系数均低于1%，强力链接触数变异系数低于3%，主要趋势在三个实现中一致再现。
3.5 均质度对AE特征的影响：m = 2时AE事件稀疏分散，仅少数孤立高能事件，无明显的局部化条带；m = 5和10时事件向试样中部倾斜区域集中，勾勒出优势斜向条带；m = 20时斜向条带连续性和事件密度进一步增强；均质条件下事件密度最高，勾勒出两条相交斜向条带。震级–频率分布显示AE活动以低中震级事件为主，随均质度增加分布向高震级延伸。以Mw ≥ ?6.80为共同比较阈值，其事件比例从m = 2的25.1%增至m = 5、10、20的48.1%、58.1%、61.3%，均质条件下达68.2%。聚类参数敏感性分析表明，虽AE事件计数和震级统计对聚类参数敏感，但两种均质度条件间的差异在所考察范围内保持稳健。
3.6 非均质性、力链、断裂与AE的多尺度对应：结果揭示了接触尺度非均质性、强力链组织、断裂发育、AE震源特征与宏观力学响应之间一致的多尺度对应关系。低均质度下较多弱接触与更早、更分散的键断裂、较不连续的强力链网络和较低峰前荷载累积相关；随均质度增加，强力链骨架更连续、峰前承受更高荷载，开裂日益局部化，AE事件沿优势断裂带集中，震级分布向较高震级偏移。本研究的独特贡献在于将非均质依赖的断裂与AE响应同内部承载骨架演化相联系：多级力链分析表明晶内强力链构成主要峰前承载网络，而晶间开裂仍是最普遍的断裂组分；矿物尺度上承载总强力链荷载最大的结构域不一定具有最高的单接触平均荷载（斜长石总荷载最大，石英单接触平均荷载最高）。
讨论与结论：研究人员指出若干局限性：在所采用的Weibull参数化下，均质度指数m的变化同时改变接触属性的离散度和期望值，m相关效应应解释为对这些组合统计变化的响应而非仅对离散度的响应；非均质数值模型尚未与具有独立表征接触尺度属性分布的物理试样直接验证；定量AE结果依赖聚类准则、模型尺度、颗粒尺寸、接触力学参数和加载条件，所报道差异应主要解释为当前数值框架内的相对趋势。结论如下：(1)接触属性非均质性主要影响峰值和峰后力学响应，对初始弹性刚度影响很小；弹性模量在各工况间变化甚微，峰值强度从m = 2的59.81 MPa增至均质条件的153.81 MPa，峰值应变从0.296%增至0.696%，峰后软化总体更突然；低均质度下较多弱接触与更早、更分散的损伤相关，均质度增加则减少早期损伤，使峰值失稳前能累积更大荷载和变形。(2)强力链网络表现出共同的形成、渗流、强化和峰后坍塌阶段演化；低均质度下承载网络较不连续、接触力较低，均质度增加则发育更连续的强力链骨架并在峰前承受更高荷载，峰值强力链荷载从m = 2的2.77 MN显著增至均质条件的约14.3 MN；峰后连续承载骨架快速解体，仅留局部残余链，伴随宏观承载能力的急剧丧失。(3)荷载传递在结构域和矿物相间强烈分异：峰前晶内强力链构成主要承载骨架，晶间链起次要作用；峰值总强力链荷载排序为斜长石 > 晶间结构 > 石英 > 碱性长石，而单接触平均荷载排序为石英 > 斜长石 > 碱性长石 > 晶间结构，表明各域总荷载不仅取决于强力链接触数还取决于每个接触传递的荷载；尽管晶内强力链主导峰前承载骨架，断裂仍主要受晶间开裂控制，但高均质度下晶内开裂再次增加，表明矿物内部在最终破坏中作用增强。(4)矩张量导出的AE特征一致反映从弥散损伤向局部化断裂的转变；低均质度下AE事件稀疏分散、低震级事件贡献更大，随均质度增加AE活动沿优势斜向断裂带日益局部化并偏向较高矩震级，Mw ≥ ?6.80事件比例从25.1%增至68.2%，较高震级事件贡献增加与AE活动更强的空间局部化同宏观失稳前日益协调的断裂发育相一致。
