基于FDM-DEM耦合的冲击损伤大理岩静态断裂力学特征研究

为探究循环冲击损伤后大理岩的静态断裂力学特征,基于有限差分（finite difference method,FDM）-离散元（discrete element method,DEM）耦合的建模技术构建了三维分离式霍普金森压杆（split Hopkinson pressure bar,SHPB）数值模型,其中杆件系统和岩石试件分别采用FLAC3D和PFC3D程序建模。利用该模型对中心直切槽半圆盘（NSCB）试样进行了恒定子弹速度下的循环冲击,随后对受损试样进行静态三点弯曲断裂实验。通过编写Fish程序,提取试样断裂面数据,对断裂面进行重构并定量计算表面粗糙度。通过与相关室内实验结果的对比分析,验证了本文数值分析的合理性与可靠性。模拟结果表明,随着循环冲击次数的增加,试样内部微裂纹、破碎颗粒均增加。连接力场分布混乱,部分力链发生断裂。力链的变化是试样力学性能劣化的根本原因。在静态三点弯曲断裂实验中,冲击5次后试样的静态断裂韧度较天然试样产生一定程度的降低。试样在静载过程中产生的微裂纹和碎块的数量随循环冲击次数的增加而增加,断裂面粗糙度随循环冲击次数的增加而增加。     

循环加卸载损伤大理岩的动力学特性

利用MTS 815电液伺服岩石实验系统进行上限应力为80%、85%、90%、95%单轴抗压强度的大理岩单轴压缩循环加卸载实验，每种上限应力条件分别设置20、40、60、80次循环。再利用分离式Hopkinson压杆对损伤岩样进行动力学实验。分析了循环加卸载上限应力及循环次数对大理岩塑性应变的影响，揭示了大理岩动态力学参数和破碎吸收能随损伤变量的演化规律。实验结果表明:塑性应变与循环次数呈正相关，且上限应力越大，塑性应变趋于稳定所需的循环次数也会增大；动态单轴抗压强度、动态弹性模量随损伤变量增加呈指数衰减；破碎吸能占比以损伤变量D=0.343为临界点分为两个阶段，D＜0.343时，破碎吸能占比稳定在10%左右，数值约为13 J，当D＞0.343时破碎吸能占比随损伤变量增加不断增大。研究结果可为岩体工程的设计、施工及支护参数的选取提供参考。