单轴循环冲击下弱风化花岗岩的损伤演化*

闫 雷，刘连生,2，李仕杰，杨道学，刘 伟

（1. 江西理工大学资源与环境工程学院，江西 赣州 341000；2. 江西理工大学江西省矿业工程重点实验室，江西 赣州 341000）

赣南地区的花岗岩风化壳是我国南部红土风化壳的一部分[1]，爆破开挖活动不可避免的经过风化壳。爆破作业过程中岩体在多次动力扰动作用[2]下极易发生失稳破坏，这是动态损伤累积的结果[3]。因此，研究风化岩石在爆炸应力波反复作用下的损伤演化规律，对揭示风化岩石动态损伤力学特性、提高风化岩体稳定性及对风化壳矿体爆破增渗技术研究[4]具有重要的理论和实际意义。

分离式霍普金森压杆（split Hopkinson pressure bar, SHPB）实验系统广泛应用于岩石动态力学特性研究[5]，许多学者开展了循环冲击下岩石力学特性的研究：金解放等[6-8]研究了砂岩在静载与循环冲击组合加载下的破坏模式，用波阻抗定义了损伤变量，并结合Logistic 方程的逆函数建立了损伤演化模型，分析了各参数的物理意义及围压和轴压对其的影响。王春等[9]对深部岩石在三维高静载下进行频繁动力扰动实验，推演了损伤演化方程和岩石本构模型并进行了实验验证。不考虑轴压和围压组合加载的情况下，朱晶晶等[10]对花岗岩进行单轴循环冲击实验，基于Weibull 分布的动态损伤本构模型分析了岩石累积损伤的演化规律，变形模量、屈服应变等随冲击次数增加均会发生变化，可以用来描述损伤演化的过程。Li 等[11]利用摆锤驱动的SHPB 装置对绿砂岩反复冲击，从声学特性、能量耗散、变形特征和微裂纹演变的角度揭示了绿砂岩损伤演化的内在机制。王志亮等[12]对黑云母花岗岩以4 种不同应力幅值进行等幅值循环冲击研究其损伤演化机理，结合裂纹起裂应力讨论了循环冲击岩石的损伤应力阈值问题。上述研究主要以深部致密岩石或完整性、均质性好的岩块为主，对具有一定风化程度的岩石，在循环冲击载荷下的力学性质和损伤演化规律的研究鲜有报道。

本文中，以弱风化花岗岩为研究对象，利用改进的SHPB 装置，以不同的冲击速度单次冲击同组岩样，并结合准静态应力应变曲线，从能量角度分析损伤应力阈值的确定方法；在损伤阈值内，以3 种速度对岩样进行等速循环冲击，利用每次冲击后的轴向最大应变研究损伤演化规律，用变形模量定义的损伤变量建立损伤演化模型，并讨论各个参数的物理意义，以期为花岗岩风化壳岩土爆破破岩或防护设计提供参考。

1 实 验

1.1 试样制备

岩样取自赣南安远县某稀土矿半风化层，为粗粒黑云母弱风化花岗岩。按照岩石力学实验规定，将岩样粗坯经过钻、切、磨加工成直径为50 mm 的圆柱体试样，静载试件长度取100 mm，动载试件长度取50 mm[13]，端面不平行度和不垂直度均小于0.02 mm。根据《GB50021—2001 岩土工程勘察标准》对花岗岩风化程度进行划分，并利用超声波纵波波速和强制饱和吸水法测得的有效孔隙度再对弱风化岩样进行筛选，保证实验样本的物理力学性质一致。为避免水对实验结果的影响，将所选试样在108 ℃下烘干48 h，置于干燥器中冷却至室温待用。采用江西省矿业工程重点实验室的RMT-150C 试验机对干燥花岗岩进行单轴压缩实验。图1 为加工好的部分试样。试样的基本物理力学参数：密度ρ = 2 389 kg/m3，有效孔隙度neff= 2.55%，纵波波速vp= 3 658 m/s，应变率 ε ˙=5.0×10−5s−1，准静态单轴抗压强度σqc= 25.15 MPa，弹性模量E = 6.942 GPa，泊松比µ = 0.257。

1.2 实验装置及原理

冲击实验在江西理工大学SHPB 实验系统上[14]完成，实验系统如图2 所示。该装置的冲头和杆件的材质均为40Cr 合金钢，密度为7 810 kg/m3，纵波波速为5 400 m/s，弹性模量为240 GPa，波阻抗为42 TPa/s；杆件直径为50 mm，入射、透射和吸收杆的长度分别为2.0、1.5 和0.5 m；采用纺锤形冲头消除P-C 振荡[15]实现恒应变率加载，冲头直径为50 mm，长度为265 mm。

图 1 加工好的部分岩石试样Fig. 1 Processed weakly-weathered granite specimens

图 2 SHPB 实验系统Fig. 2 SHPB experimental system

基于一维应力波和动态平衡假定，结合采集到的入射应变信号εi(t)、反射应变信号εr(t)和透射应变信号εt(t)，对脆性岩石用三波法[16]计算其动态应力σ(t)、应变ε(t)和应变率 ε ˙(t) 能够得出准确结果，分别为：

式中：A、L 分别为试样的横截面积和长度；Eb、Ab分别为压杆的弹性模量和横截面积；cb为压杆纵波波速。

1.3 实验方案

为了研究弱风化花岗岩在爆炸应力波作用下的损伤演化机理，采用SHPB 实验系统对其施加一维循环冲击载荷。由于风化壳内岩石埋深较浅，可不考虑围压和轴压。通过调节气室气压和冲头在发射腔的位置来控制冲头速度。将所有试样以5 个/组，分为单次冲击实验组A 到E，冲击速度由低到高，组内以字母加数字的形式编号，例如A1、D2；等速循环冲击组F、G 和H，组内以字母加数字的形式编号，例如F1-4 为F 组1 号岩样第4 次冲击。（1）以不同的速度对一组试样进行单次冲击直到发生宏观破裂，确定产生有效损伤的冲击速度范围；（2）以不同的有效损伤冲击速度对试样进行等速循环冲击至裂实验。在试样/压杆接触面涂抹润滑油来降低端面摩擦效应[17]。每次冲击后的动态应力平衡检验是必要的[18]，并以此来筛选有效实验数据，图3 为E4 岩样在7.89 m/s 冲击下的应力平衡图。

图 3 岩样E4 动态力平衡检验Fig. 3 Dynamic stress balance check for specimen E4

2 岩石损伤阈值确定

岩石在循环冲击荷载下，存在一个开始出现损伤的阈值。葛修润等[19]通过岩石周期疲劳实验和CT 细观分析得出疲劳损伤门槛值对应裂纹起裂应力。王宇等[20]结合声发射事件对单轴压缩起裂机制进行了分析，并用裂纹体积应变拐点确定起裂应力，具有明确的物理意义。梁玉昌等[21]研究得出在中低应变率范围内花岗岩的静态起裂应力和峰值强度的比值与应变率无明显关系，王志亮等[12]以此为基础利用Nicksiar 等[22]提出的“侧向应变响应法”确定起裂应力，提出动态起裂应力近似等于静态起裂应力与动态增长因子[23]的乘积，作为循环冲击的损伤阈值。Li 等[24]使用大直径SHPB 对花岗岩反复冲击研究其损伤特性时，以该花岗岩静态应力-应变曲线线性偏离点，即损伤强度作为损伤冲击门槛值。由上可知，起裂应力在中低应变率范围内作为岩石损伤阈值是可行的，在高应变率加载下，动态起裂应力与动态抗压强度的比值的率无关性有待进一步实验验证；直接用准静态损伤强度作为动态冲击损伤阈值，忽略了动态峰值应力、应变的率效应。以静态损伤加载能量为阈值，可避免高应变率下岩样具有高动动态峰值强度σdc对应低动态峰值应变εdc时，分别以应力或应变为阈值判定的损伤结果相矛盾。

单轴压缩条件下，岩样达到静态损伤强度σcg吸收的体积能量Wcg为其应力应变曲线下所包含的面积[25-26]，如图4（a）所示；动态体积加载能Wdc为其动态应力-应变曲线加载段面积，如图4（b）所示。可通过定积分的方法计算，即：

式中：σi和εi分别为应力应变曲线上的任一点的应力及其对应的轴向应变（其中i=1, 2, 3, ···），当i=0时σi=εi=0；将x=cg, dc 代入式(4)中，分别计算Wcg和Wdc，单位为J/cm3。

图 4 加载段单位体积能量计算Fig. 4 Volume energy calculation for load segment

由上述分析得到，Wcg=0.041 0 J/cm3为循环冲击损伤下限，由图4（b）可得，当冲头速度v ≥ 4.06 m/s时，Wdc＞Wcg；对一组试样以不同速度进行单次冲击，以冲击后产生宏观裂纹所对应的速度为损伤上限，实验结果见图5 和表1。表1 中neff为有效孔隙度，L/D 为试样长径比，ρ 为试样密度，v 为冲击速度，vl为纵波波速，σdc第一次冲击动态峰值应力，为第一次冲击平均应变率，N 为累积冲击次数，Wi为入射能量。

图5(a)给出了不同应变率下岩样的应力-应变曲线。当冲击速度为6.94 m/s 时，岩石破裂，并以此速度为损伤上限。在冲击载荷下，弱风化岩石损伤大致可分为弹性阶段、裂纹扩展阶段和卸载阶段3 个阶段；相对准静态加载而言，其压密阶段不明显但存在。其中，卸载阶段以最大应变为界又可分为第1 和第2 卸载段[11]，第2 卸载段均出现应变回弹，是由于岩样在卸载过程中未完全破碎(图5(b)、(c))仍具有一定的承载能力，并伴随弹性应变能释放而造成的。裂纹扩展阶段有两种表现形式，即近似线性的稳定增长和下凹型的非稳定增长，且随应变率升高下凹段越明显。这是由于风化岩石内部空隙发育，加载初期微裂隙压密，骨架和孔洞共同承力进入弹性变形阶段；随应力增高，部分孔隙坍塌并再次平衡，进入稳定扩展阶段；随着总的入射能增大，超过大部分孔洞的屈服极限，在弹性阶段后大量孔洞坍塌出现应力松弛平台段，耗散大部分能量[27]。同时，在发生脆性断裂前出现短暂的应变硬化，可理解为孔洞坍塌后应力迅速在其附近集中，达到岩样整体的屈服强度后突然断裂破坏。

图 5 不同速度单次冲击的应力-应变曲线及破坏形式Fig. 5 Stress-strain curves and failure modes of the specimens subjected to single impact at different velocities

表 1 单次冲击和首次循环冲击实验结果Table 1 Single impact and first cycle impact test results

3 损伤演化

3.1 循环冲击实验结果

由上节可以确定有效冲击速度在4.06 ～ 6.94 m/s 之间，分别以近似4、5 和6 m/s 的速度对岩样进行等速循环冲击，实验结果见表1，图6 为岩样在不同冲击速度下的动态应力-应变曲线及其破裂形态。

从图6 可以看出，随着冲击速度的升高，初始孔隙度相近的岩样达到破裂的冲击总次数明显减少，最后一次冲击破裂的最大应变也有显著的增长。经过多次冲击后，岩样最终出现一条贯穿裂纹，表现出随着冲击次数增加的渐进性扩展（图6(a)）和突然贯通（图6(b)和(c)）两种。裂纹扩展段的应力松弛平台，也会随着冲击速度和次数的增加而凸显。动态峰值应力和最大应变与冲击速度正相关；随冲击次数的增多，峰值应力递减，应变率呈先升后降的趋势；峰值应力的变化是由于每一次冲击都会对岩样造成有效损伤，使其力学性质不断劣化，第二次冲击后峰值应力偶有增大，是冲击速度控制误差造成的；应变率表示应变的速率，与变形对时间的一阶导数。冲击前期力作用时间一定，变形不断增大，应变率升高；冲击后期岩石内部积累了大量的微裂隙，卸载后变形不会立即停止，即变形时间延长，应变率降低。从最后一次冲击的应力-应变曲线也可以看出，在卸载阶段出现应变软化，且末尾有回弹现象，说明卸载后变形依然持续，使得应变率降低。由此可见，应变率不仅对加载应力敏感，也与损伤程度相关。

图 6 不同速度循环冲击岩样的应力-应变曲线及其破坏形式Fig. 6 Stress-strain curves and failure modes of different specimens subjected to cyclic impact at different velocities

3.2 损伤演化规律

为了描述循环荷载下岩石的疲劳损伤演化规律，Xiao 等[28]推导了完整岩石的倒S 形疲劳损伤模型，可分为3 个阶段：初始加速、稳定发展和加速扩展阶段；Liu 等[29]对间隔节理岩石进行循环单轴压缩，不可逆应变也以3 个阶段的倒S 形发展；金解放等[8]研究了岩石在静载荷和循环冲击组合加载下的损伤演化规律，同样可分为3 个阶段。由上可知，岩石的疲劳损伤到破坏基本遵循倒S 演化模型。本文中，每次冲击后轴向最大应变εmax和循环冲击次数n 也遵循这一规律，如图7 所示。

图7 显示了不同冲击速度冲击下，最大应变表示下岩石的疲劳损伤演化过程。可以看出，在低速多次冲击下岩样破坏时的最大应变最小，且随冲击速度增高，稳定发展阶段的速度越快、历时越短，最大应变值增大。用三次函数拟合可获得良好的效果，但没有明确的物理意义，需要建立简单合理的累积损伤演化模型。

图 7 不同冲击速度下不同岩样最大轴向应变的演化Fig. 7 Evolution of the maximum axial strain with cyclic-impact number at different impact velocities for different specimens

3.3 损伤演化模型

建立可靠的疲劳累积损伤模型来评估、预测岩石在频繁动力扰动下的寿命，首先要定义合理的损伤变量。Xiao 等[30]列举了多种疲劳损伤变量的定义方法，包括弹性模量法、最大应变法、能量耗散法、残余应变法、超声波速法和声发射累积计数法等，都能用来描述损伤的演化。本文中选用弹性模量定义损伤变量，通常取峰值应力的40%～60%之间的割线斜率，或峰值应力的50 %处切线斜率作为弹性模量计算损伤。从图6 中可以观察到，风化花岗岩的应力应变曲线具有明显的分段性，且随冲击次数的增长变化明显，因此不能用上述方法直接计算弹性模量。这里采用弹性阶段的切线模量Een和裂纹扩展阶段割线斜率Epn(见图5(a)，其中n 表示冲击次数)来描述损伤的演化，即：

式中：Ee0和Ep0分别取同批岩样中孔隙度最小、完整性最好的动态应力应变曲线[31]上的相应值，也是同批岩样中Een和Epn的最大值，分别为53 GPa 和21.5 GPa，以此作为弱风化花岗岩的初始变形模量。

岩石在循环加载下的损伤演化基本满足倒S 形，Xiao 等[28]推导的疲劳损伤累积模型计及了岩石初始和第一次加载到最大应力所造成的损伤总和，记为D0，不能单一量化首次加载所产生的损伤。金解放等[8]用Logistic 方程的逆函数建立了损伤累积演化模型，并分析了各个参数的物理意义，但没有讨论过零点的意义，且模型在初始和加速阶段损伤累积速率太快，难以反映低冲击载荷扰动下损伤累积初始和破坏前渐进的演化过程。假设岩样冲击前损伤量D0= 0，这里讨论的初始损伤仅由冲击扰动产生。由其他非冲击扰动因素造成岩样的初始损伤，可通过孔隙度、密度、纵波波速等无损检测得到的宏观指标来保持实验组岩样未冲击前损伤程度近似相同，进而可不考虑；冲击损伤量用式（5）计算，破坏时DN= 1。

由3.2 节知，岩石循环冲击下累积损伤演化可用一元三次函数f (x) = αx3+βx2+ηx+C 表示，当α＞0 时，f (x) 在(−∞, +∞) 内呈倒S 形。自变量x 用相对循环次数n/N 表示，因变量y 用D 表示；由D(0)=0 得C=0，由D(1)=α+β+η=1 得η=1−α−β，整理得：

其中，令D"=0，可得n/N=β/(3α)，即该点为损伤减速累积与加速累积的拐点，将其代入式（6）得到拐点对应的损伤值近似等于|β|/α−1，相对差值在4 %之内。因此，可用β/(3α)表示达到中值点的相对循环次数，|β|/α−1 示中值点累积损伤程度，图8 显示了两种损伤变量表示的损伤演化模型。

由图8 可以看出，Dp的值明显高于De的值，当Dp＞0.8 或De＞0.4 时均进入加速破坏阶段，不同的损伤变量表示的演化模式差异较大，对于风化岩石选取统一的损伤变量是预测疲劳寿命的必要前提。不同冲击速度下岩样损伤演化拟合相关系数均大于0.97，拟合效果理想。

图9 分析了参数α 和β 对损伤演化模型的影响，随α 值增加，初始阶段斜率变缓，加速阶段斜率变陡，控制着初始和加速阶段损伤累积的速率，故参数α 可定义为始、末段损伤累积速率因子；随β 值减小，(β/α)−1 增大，因此，参数β 可表征稳定发展阶段损伤累积程度。综上所述，本文提出的损伤演化模型可以较好地模拟不同冲击速度下岩石损伤演化的3 个阶段；双参数控制的模型，关系式简单、物理意义明确，可以估算中值点的损伤程度及其对应的冲击次数，指导工程实践。

图 8 不同损伤定义下的损伤演化模型Fig. 8 The damage evolution models represented by different damage variables

图 9 参数α 和β 对损伤累积模型的影响Fig. 9 Effects of parameters α and β on the damage accumulation model

4 结 论

（1）根据单次冲击作用下弱风化花岗岩力学特性和破坏形式，结合准静态加载时达到损伤强度吸收的单位体积能和动态单位体积加载能的大小，讨论了用能量确定损伤阈值的方法，并通过确定产生有效损伤的冲击速度范围在4.06～6.96 m/s 之间。

（2）不同应变率下弱风化花岗岩应力应变曲线可分为弹性、裂纹扩段和卸载3 个阶段；随应变率增加，弹性模量增大、裂纹扩展阶段割线斜率呈先增后减趋势，与新鲜致密岩石不同该阶段随应变率增大下凹越明显，是由于孔洞坍塌出现的应力松弛；峰值应力呈递增的趋势，峰值应变后期增加不明显。

（3）不同速度循环冲击下，弱风化花岗岩应力应变曲线也具有3 个阶段，随冲击次数的增加弹模、割线模量均减小；在等速循环冲击过程中，由于损伤累积使得峰值应力递减，应变率呈先增后减的趋势；裂纹扩展段的应力松弛平台，随冲击次数和冲击速度的增加而愈发明显；最终形成一条贯穿的张开裂纹，岩样整体仍保持完整。

（4）循环冲击作用下弱风化岩石损伤演化具有初始加速、稳定发展、加速扩展3 个阶段，且随冲击速度增大，这种趋势越不明显；用弹性阶段的切线模量、裂纹扩展阶段割线斜率定义的损伤变量与相对循环次数建立了的双参数损伤演化模型；讨论了待定参数的物理意义，其中α 为始、末段损伤累积速率因子，β 为稳定发展阶段损伤累积程度因子，且中值点的相对循环次数β/(3α)和中值点累积损伤度(β/α)−1与冲击速度呈正相关。