夜雨聆风一、引言：循环加卸载下的岩石破裂，你真的搞懂了吗？
在矿山开采、隧道掘进、深部资源开发这些工程现场，岩石往往不是一次性被压垮的，而是经历反复加卸载的"折磨"后才最终破坏。这种循环荷载作用下的岩石力学行为，跟单次加载完全不是一回事儿。
很多做数值模拟的同行在这块儿踩过不少坑：模拟出来的应力-应变曲线跟试验对不上、裂纹演化规律看不清楚、破坏形态跟实际试件差十万八千里……这些问题归根结底，往往是建模思路和参数标定没到位。
▎ 循环路径影响大 — 恒下限和变下限两种循环方式，砂岩的破坏模式截然不同
▎ 裂纹演化难捕捉 — 张拉裂纹和剪切裂纹的萌生扩展规律需要精细化追踪
▎ 模拟验证不可少 — PFC模拟结果必须与室内试验对比验证才能放心用
▎ 位移场分析关键 — 颗粒位移场是揭示细观破裂机制的核心抓手
今天这篇干货，我就结合一篇PFC3D数值模拟成果，把砂岩在单轴压缩、恒下限循环、变下限循环三种工况下的破裂演化过程，从头到尾给你盘一遍。
二、研究概况与模型建立
▶ 2.1 数值模型搭建
本研究采用PFC3D（Particle Flow Code 3D）离散元软件，建立尺寸为50mm×100mm的圆柱形砂岩数值模型。模型共包含约52544个颗粒，颗粒间微键数约301825个。颗粒粒径控制在0.8mm~0.96mm均匀分布，最短边与平均粒径之比大于5，有效消除了粒径效应。
PFC3D数值模型与颗粒接触示意
关键建模参数一览 模型尺寸：φ50mm × H100mm（与室内试验一致） 颗粒数量：约52544个 颗粒粒径：0.8mm ~ 0.96mm 接触模型：平行黏结模型（Parallel Bond Model） 颗粒密度：2396 kg/m³ 颗粒有效模量：7.61 GPa |
接触模型选用平行黏结模型，这种模型在颗粒黏结时能抵抗扭矩并表现为线弹性，当受力超过极限后黏结破坏产生微裂纹，非常契合砂岩这种质地均匀、硬度较大的刚性颗粒集合体特征。
▶ 2.2 细观参数标定
颗粒微观参数与材料宏观力学性能之间没有定量关系，因此需要通过"试错法"进行标定。标定思路是：调整pb_emod（平行黏结模量）、emod（有效接触模量）、kratio（法向/切向刚度比）等参数，使模拟得到的抗压强度、弹性模量、泊松比与室内试验匹配。
细观参数 | 数值 | 细观参数 | 数值 |
最小颗粒半径Rmin | 0.80 mm | 平行黏结半径系数λ | 1.00 |
Rmax/Rmin | 1.20 | 平行黏结有效模量pb_emod | 0.78 GPa |
颗粒体密度ρ | 2396 kg/m³ | 法向黏结强度pb_ten | 15.69 MPa |
粒间摩擦系数μ | 0.70 | 切向黏结强度pb_coh | 7.80 MPa |
颗粒有效模量emod | 7.61 GPa | 法向-切向刚度比kn/ks | 0.70 |
标定完成后，单轴压缩模拟强度为63.482MPa，与室内试验的63.140MPa非常接近；巴西劈裂模拟也以张拉破坏为主，裂纹贯穿圆盘，与试验结果基本一致。
单轴压缩模拟与室内试验破坏形态对比
巴西劈裂模拟与室内试验结果对比
三、单轴压缩模拟：裂纹演化的基准线
在进行循环加卸载模拟之前，先做单轴压缩模拟作为基准对照。设置24120个加载步，记录应力、应变、裂纹数等变量的演化过程。
▶ 3.1 应力-应变响应
从应力-应变曲线可以看出，模拟得到的抗压强度63.482MPa与室内试验吻合。但峰前区斜率变化不大，没有出现明显的压密阶段——这是因为PFC模型中颗粒接触紧密、分布均匀，没有真实砂岩内部的孔隙和缺陷。
单轴压缩应力-应变对比曲线
▶ 3.2 裂纹演化规律
裂纹演化是PFC模拟最精彩的部分。单轴压缩下，总裂纹数在加载初期几乎为0，当应力达到约42.28MPa（峰前66%）时裂纹开始迅速增加，峰值应力时裂纹数达到8457500个。其中剪切裂纹占绝对主导（8410000个），张拉裂纹仅47500个。
单轴压缩下应力与裂纹数随时步关系曲线
单轴压缩下应力与剪切裂纹随时步关系曲线
单轴压缩核心结论 模拟抗压强度63.482MPa，与室内试验63.140MPa高度一致 裂纹以剪切破坏为主，剪切裂纹占比超过99% 裂纹在峰前约66%应力水平开始加速扩展 峰值应力时总裂纹数约846万个，试件瞬间贯通破坏 |
▶ 3.3 位移场与破坏模式分类
根据颗粒相对位移的大小和方向，将裂纹形成规律分为六种模式：直接剪切（DS）、间接剪切（RS）、直接拉伸（DT）、间接拉伸（RT）、混合破坏MF-1~MF-4。
不同裂纹模式下颗粒位移场分类示意
单轴压缩峰值应力状态的位移场显示，主裂纹产生位置以混合破坏（MF-1）和间接剪切（RS）为主，边缘区域以间接拉伸（RT）和间接剪切为主。这种位移场分布直接决定了最终的X型剪切破坏形态。
单轴压缩峰值应力状态颗粒位移场分布
四、恒下限循环加卸载：强度提升与X型贯通破坏
恒下限循环加卸载是指每次循环都从0加载到逐步增大的上限应力，下限始终保持为0。这种路径下，砂岩经历了7个循环阶段才最终破坏。
▶ 4.1 应力-应变与滞回特性
模拟得到的抗压强度为68.89MPa，比单轴压缩提高了约8.5%。应力-应变曲线表现出明显的滞回现象：随着循环次数增加，滞回环面积逐渐增大，反映出岩样内部损伤持续加深。这一演化特征与室内试验完全吻合。
▶ 4.2 裂纹演化特征
恒下限循环的裂纹演化非常有规律。第1~4循环阶段裂纹增长较缓，从第5循环开始裂纹数目增长速度明显加快。当达到峰值应力68.89MPa时，总裂纹数达到1735万个，其中剪切裂纹1730万个，张拉裂纹仅5万个。
恒下限循环加载下应力与裂纹数随时步关系曲线
值得注意的是，在卸载阶段，新裂纹的产生会受到抑制，但在已有裂纹尖端由于应力集中，仍会出现新的裂纹。这说明循环加载的损伤累积效应非常显著。
▶ 4.3 破坏形态与位移场
恒下限循环加载下，试件呈现典型的X型贯通剪切破坏，与室内试验结果高度一致。从位移场可以看出，空间分布呈典型"X"形状，主裂纹尖端法向以混合破坏（MF-1）为主，切向以间接剪切（RS）为主。
恒下限循环加载模拟与室内试验破坏形态对比
恒下限循环加载峰值应力状态位移场分布
恒下限循环核心结论 抗压强度68.89MPa，比单轴压缩提高约8.5% 滞回环面积随循环次数逐渐增大，损伤持续累积 第5循环后裂纹增速明显加快，呈加速破坏特征 破坏模式为典型X型贯通剪切破坏，与试验吻合 |
五、变下限循环加卸载：未完全贯通的X型破坏
变下限循环加卸载是指上下限应力都随循环逐步增大，应力幅值不断变化。这种路径下，砂岩的破坏特征与恒下限有明显差异。
▶ 5.1 应力-应变响应差异
变下限循环模拟得到的抗压强度为65.7MPa，比单轴压缩提高约4.1%，但低于恒下限循环。应力-应变曲线整体形态与室内试验相近，但数值模型在达到峰值后呈现快速应力跌落，而室内试验由于碎屑填充增强了界面摩擦，强度有所提升。
从滞回特性来看，室内试验在初始加载阶段就表现出显著的滞回特性，而数值模型由于理想均质假设，初始阶段损伤积累有限，滞回环面积始终维持在较低水平。
变下限循环加卸载应力-应变曲线对比
▶ 5.2 裂纹演化与破坏特征
变下限循环下，试件在第7个加卸载循环发生破坏，总裂纹数714万个，其中剪切裂纹710万个，张拉裂纹4万个。与恒下限相比，变下限的总裂纹数明显更少，说明其损伤累积程度相对较轻。
变下限循环加载下应力与裂纹数随时步关系曲线
破坏形态上，变下限呈现未完全贯通的"X"型破坏，主裂纹没有彻底贯穿整个试件。这与室内试验观察到的现象一致。
变下限循环加载模拟与室内试验破坏形态对比
▶ 5.3 位移场分析
变下限位移场呈现未完全贯通的"X"型分布模式。在主裂纹方向上，位移场以混合破坏（MF-1）和间接剪切（RS）为主导；在边缘处主要以间接拉伸（RT）为主。试样中部颗粒的横向位移与上下加载纵向位移形成破裂面，这种位移场的空间分布直接决定了主裂纹的萌生与扩展路径。
变下限循环加载峰值应力状态位移场分布
变下限循环核心结论 抗压强度65.7MPa，比单轴压缩提高约4.1% 滞回环面积始终较低，损伤累积程度轻于恒下限 破坏模式为未完全贯通的X型破坏 边缘区域以间接拉伸为主，中部以混合破坏和间接剪切为主 |
六、避坑指南与实战建议
结合这次PFC3D模拟的实战经验，我给大家总结几条做岩石离散元模拟时最容易踩的坑，以及对应的解决思路。
实战建议一：参数标定要耐心 平行黏结模型的细观参数与宏观力学性能没有解析关系，必须通过试错法标定。 建议先固定kratio（刚度比）和λ（黏结半径系数），再微调pb_emod和pb_ten。 每次只改一个参数，观察对应宏观指标的变化，建立参数敏感性认知。 |
实战建议二：循环路径设置要合理 恒下限和变下限的加载策略完全不同，伺服控制代码要分别编写。 建议在每个循环的峰值点和谷值点都保存模型状态，方便后续分析。 加载步长不宜过大，否则容易错过裂纹萌生和扩展的关键时刻。 |
常见误区提醒 不要指望PFC模拟的应力-应变曲线跟室内试验完全重合，峰前压密阶段很难复现。 颗粒数量不是越多越好，要在计算精度和时间成本之间找平衡。 位移场分析时，注意区分MF-1~MF-4不同混合破坏模式的判别标准。 |
建模效率提升技巧 利用PFC的history函数自动记录应力、位移、裂纹数等变量，减少后处理工作量。 裂纹分类统计时，建议同时输出张拉裂纹和剪切裂纹的时程曲线，便于对比分析。 位移场可视化时，可将颗粒按位移大小着色，配合箭头方向展示运动趋势。 |
如有其他需要，欢迎关注我的咸鱼号:Orangetu
本推送内容基于公开文献进行梳理与解读，著作权归原作者及出版机构所有。内容仅代表个人对原文的理解与归纳，旨在学术交流与知识传播，不作商业用途。如需引用或进一步探讨，请务必参阅原始文献，若涉及版权事宜请联系我们处理。
感谢您的关注，麻烦多多对推文点赞、收藏和转发，并衷心希望您多多指教！帮助CAE数值仿真模拟进步提升！
