夏 晨，戚承志，利 学，周卓群

(1. 北京建筑大学土木与交通工程学院，北京100044；2. 北京未来城市设计高精尖创新中心，北京100044；3. 城市交通基础设施建设北京国际科技合作基地，北京100044)

在外部加载条件下，脆性材料将产生新的裂纹，而初始裂纹将会进一步发展，最后贯穿并破坏材料。为了研究这类材料中裂纹的演化机理，已有多位学者提出并发展了I 型裂纹模型[1−8]。随着轴向荷载的增大，主裂纹顶部会产生次生裂纹，而次生裂纹将呈现曲线扩展，最终平行于加载方向。

HUANG 等[9]分析了侧限以及应变率效应对于岩石材料动力强度的影响。由于缺乏对材料力学性质退化以及惯性效应的影响，其研究结果不符合实验结果。BHAT 等[10]在前人模型[4− 11]的基础上，引入损伤演化率，扩大了该模型的适用范围，但是该模型没有考虑动力加载时裂纹间的相互作用。

在单轴压缩加载条件下，I 型裂纹模型可以表达具有稳定裂纹扩展的非线性应力-应变关系[5−7]。该模型不能确定材料的压缩强度极限，为了改进这个缺陷，部分学者提出了一些可考虑裂纹相互作用的模型[12−19]，试件在动力荷载作用下，裂纹间的相互作用强度取决于裂纹的初始尺寸、走向、密度及其空间分布。ASHBY 和SAMMIS[4]通过引入裂纹间的桥连附加内力以考虑裂纹相互作用效应，然而，该方法是一种近似的处理方式。KACHANOV 等 [12 − 14]以及FABRIKANT[15− 16]提出了更一般的裂纹间相互作用的处理方法，KACHANOV 的方法基于应力叠加原理来求解相互作用函数。根据KACHANOV 考虑裂纹间相互作用的方法，ZHENG 等[17]揭示了材料的整体破坏和率效应的机理，但是其研究未考虑裂纹扩展速度对断裂韧性的影响以及弹性模型退化的影响。QI 等[20]考虑裂纹扩展速度、弹性模量退化以及裂纹相互作用等因素，提出了岩石类材料动力加载下的应力-应变本构关系。此模型未考虑裂纹面滑动过程中，滑动速度对裂纹面摩擦系数的影响。研究表明：裂纹面摩擦滑动对裂纹尖端应力强度因子和裂纹发展有重要影响[21−24]。

本文将考虑裂纹面滑动过程对裂纹面摩擦行为的影响，提出一种考虑裂纹相互作用、裂纹密度以及裂纹面滑动摩擦的裂纹动力扩展模型，用以研究轴向动力加载过程中脆性材料的裂纹扩展和动力行为。

1 轴向加载下岩石的三维动力模型

1.1 惯性效应引起的附加应力

1.2 轴向加载作用下的裂纹扩展三维模型

假定体积为V的试样中，含有N个币状的初始裂纹，其平均直径为2a。根据ASHBY 和SAMMIS[4]，有效弹性模量可以表达为：

式中：Eeff和E分别为材料的有效弹性模量和弹性模量；D为损伤因子，可由式(4)求解：

式中：l为次生裂纹的长度； α=cosψ ， ψ为主裂纹平面与加载方向的夹角；Cr为单位体积内的裂纹数，Cr=N/V。

根据胡克定理，试件z方向的应变可以表达为：

式中，νeff为有效泊松比。假定有效弹性模量与体积模量相等，则有：

将式(1)～式(4)代入式(8)，试件的轴向应力-应变关系如下：

式(9)右侧第一项为动力加载项，第二项为应变加速度的贡献。根据式(9)，可研究动力加载过程中，裂纹相互作用、裂纹密度以及裂纹面滑动摩擦行为等不同因素对岩石材料的裂纹扩展、试件受力变形和强度的影响。

1.3 裂纹尖端运动方程

通过求解式(14)，可以求得裂纹扩展速度以及裂纹长度。

假定主裂纹和次生裂纹的半径分别为a和l，裂纹中心的间距为2b，如图2 所示。依据ASHBY和SAMMIS[4]的研究，次生裂纹的应力强度因子表达为：

图2 I 型裂纹受力状态及单个I 型裂纹示意图[27]Fig. 2 Stress condition for mode I crack and diagram of one isolated crack[27]

1.4 裂纹的动力相互作用模型

不同于ASHBY 和SAMMIS[4]引入裂纹间的桥连附加内应力来近似考虑裂纹间相互作用，LEKESIZ等[28]基于Kachanov 的研究，提出针对I 型裂纹的裂纹相互作用计算方法：

1.5 裂纹壁摩擦系数与滑移速度的相关性

在研究断层活动以及与摩擦相关的岩石破裂问题时，常采用速度及状态依赖性摩擦本构模型(Rate-and state-frictional law)[29− 31]描 述 其 摩 擦 行为。次生裂纹受到初始裂纹滑动的驱动，亦可以采用速度及状态依赖性摩擦本构关系进行研究。裂纹摩擦系数可以表达为：

实验研究结果[32−33]表明：当滑移速度大于0.1 m/s 后，摩擦系数会快速降低，这与速度弱化和非稳态滑移的初动是密切相关的。岩石类别对摩擦系数与滑动速度的依赖性影响不大，可以将滑动分为三个速度区间，如图3 所示。

图3 不同岩石材料摩擦系数与滑移速度相关性[32−33]Fig. 3 Dependence of friction coefficient on slip velocity[32−33]

为了统一这三个速度区间，SPAGNUOLO 等[34]提出一种修正的摩擦定律(Modified friction law, MFL)：

为了确定摩擦对裂纹扩展的影响，需要先求解裂纹面的相对滑移量，进一步可求得相应的摩擦面滑动速度。图4 展示了单一I 型裂纹的张开。N、n、T以及t分别是单位矢量，其方向如图4所示。AA'为初始裂纹，其长度为2a；次生裂纹长度为l。驱动初始裂纹滑动的力如图4(b)所示，其中：F为驱动力；Q为抵抗力；而FC为弹性闭合力；R为沿裂纹均匀分布的开裂力。滑动会沿初始裂纹AA′产生相对位移uW。而裂纹间相互作用会引起I 型裂纹上均匀分布的附加荷载，进而引起初始裂纹的相对位移uF，如图5 所示[35]。

图4 I 型裂纹初始裂纹面滑移时的受力状态[35]Fig. 4 Slip over pre-existing crack as a driving mechanism for wing crack growth[35]

图5 I 型裂纹张开时的相互作用效应[35]Fig. 5 Effect of crack interaction on wing crack opening[35]

因此，初始裂纹面的总相对滑动量可表达为：

对式(22)所求解的总位移量对时间求导，即可得到初始裂纹面的滑移速度，进而通过式(21)求解裂纹面的动摩擦系数。

2 不同裂纹面摩擦行为对动力加载过程的影响

根据ZHANG 等[36]的霍普金森杆试验(SHPB)，试样的半径为37 mm，弹性模量E和泊松比ν分别取17.2 GPa 和0.19，初始裂纹长度为3 mm。入射应力脉冲持续时间与撞击杆长相关，可以表达为：

式中：t/µs 为时间； σi为入射应力脉冲的最大幅值，与撞击速度相关。对于半径37 mm 试件的SHPB 试验，入射应力脉冲的最大幅值与撞击速度的关系如下：

当撞击速度为15.7 m/s 时，入射应力脉冲的最大幅值为128.04 MPa。假设试件的轴向应变值与入射应力脉冲值成正比，即：

根据上述方程及参数，可求得试件的轴向应变、应变率和应变加速度，其结果如图6～图8 所示。

图6 轴向应变随时间的变化Fig. 6 Axial strain vs time

图7 轴向应变率随着时间的变化Fig. 7 Axial strain rate of the specimen vs time

图8 轴向应变加速度随时间的变化Fig. 8 Variation of axial strain acceleration of the specimen with time

QI 等[20]研究了在摩擦系数为常数情况下，不同裂纹密度对裂纹的启裂应力、动力加载过程、动强度以及惯性效应的影响。但是实际观测和试验结果表明：裂纹面上的摩擦系数与滑移速度相关，在裂纹面滑动过程中并非常数。下面将采用第1 节提出的模型，以ZHANG 等[36]的试验实验数据为参考，采用Matlab 进行数值计算，分析动力摩擦过程对动力加载过程中岩石材料行为的影响。

2.1 动摩擦条件下裂纹密度的影响

裂纹面上发生相对滑动，其摩擦行为可以采用速度及状态依赖型摩擦本构模型进行研究。本节将采用式(21)描述的摩擦本构模型，以正方形分布的裂纹阵列(如图9 所示)为例，考虑裂纹间相互作用，分析不同裂纹密度条件对裂纹启裂应力、动力加载过程、最高强度以及惯性效应的影响。

图9 正方形分布裂纹阵列示意图[31]Fig. 9 Periodic square configuration cracks[31]

a2−b2、p和Vc分别取0.01、0.8 以及0.1。V0=1 µm/s，μ0(V0)≈0.7。为考虑不同裂纹密度的影响，假定 λ=a/b分别为0.15、0.30、0.50、0.70、0.8、0.9 和0.95。其余主要计算参数如表1 所示。联合式(9)和式(18)，可计算不同裂纹密度条件时试件的动力加载过程、破坏强度以及惯性效应的影响等，计算结果如图10～图12 和表2 所示。在图10 中，“S-0.9”代表裂纹密度参数为 λ=0.9的正方形阵列裂纹结果，文中其他相关表示均遵守这种规则。

图1 圆柱体试件坐标示意图Fig. 1 Diagram of coordinate axes for cylindrical sample.

图10 裂纹密度参数为0.9 时的裂纹面滑移速度Fig. 10 Slip velocity vs time for 'S-0.9' cracks

表1 数值计算采用的参数Table 1 List of the parameters for the numerical calculations

从图11 和表2 可以看出，随着裂纹密度从0.15 增大至0.95，裂纹互相作用逐渐增强，相互作用系数从1.0008 逐步增长至1.4218。而裂纹的启裂时刻从64×10−6s 逐步提前至59×10−6s，启裂时的加载应力由3.46 MPa 降低至2.29 MPa。这个过程中，试件破坏的时刻亦不断提前，破坏时的加载应力从199.18 MPa 大幅下降至22.21 MPa。当 λ=0.15时，裂 纹 扩 展 持 续 了60×10−6s；当λ=0.95时，裂纹扩展仅持续了29×10−6s，为裂纹密度0.15 时的一半。由图12 可见，惯性效应引起的试件轴向附加应力随裂纹密度的增大而逐渐增大。惯性作用引起的附加轴向应力峰值出现在为裂纹密度0.5 的试件中，最大值为1.83 MPa。相对于QI 等[20]摩擦系数为常数的结果，动摩擦过程会促进裂纹发展，并降低试样的强度。

图11 不同裂纹密度条件下轴向应力随加载时间变化Fig. 11 Axial stress vs time for different crack concentration

图12 不同裂纹密度条件下惯性引起附加轴向应力随加载时间变化Fig. 12 Inertia induced additional axial stress vs time for different crack density

表2 不同裂纹密度条件下启裂和破坏时刻以及相应的加载应力Table 2 The moments of crack growth initiation, sample failure and corresponding loading stresses for different crack density

2.2 动强度因子结果对比

动强度因子(Dynamic increase factor, DIF，为试件破坏时的动强度与精力强度之比)常用于研究应变率效应对岩石材料动力强度的影响。目前，已有众多学者通过实验研究了应变率处于1 s−1～103s−1之间时岩石材料的DIF。为验证上述提出的模型，将第2.1 节的理论结果与SHPB 实验结果[36− 39]对比在图13 中。根据ZHENG 等[17]的研究，当材料的受压强度十倍于受拉强度时，裂纹密度参数为0.3。因此，选取了裂纹密度参数为0.3 时，采用MFL 的理论值，在图13 中用黑色星代表。在图13 中，动强度因子均随应变增大而增大，表现出明显的应变率效应。同时，在400 s−1～700 s−1应变区间，随着试件直径从12.7 mm 增大到74 mm，动强度因子亦呈现增大趋势，表现出明显的动力尺寸效应。

图13 不同尺寸试样的动强度因子对比Fig. 13 DIFs of concrete-like materials obtained from this study and SHPB specimens with different diameters

2.3 不同摩擦系数的结果对比

为进一步研究摩擦效应对动强度的影响，将摩擦系数假定为常数，分别取值0.05 和0.7。裂纹密度参数λ 分别取0.3 和0.95，以考虑不同裂纹密度的影响。采用MFL 时，裂纹面开始滑动前，摩擦系数假定为0.7。其他参数与第2.1 小节保持一致，对比结果如图14 所示。

图14 不同摩擦系数条件下轴向应力随加载时间变化Fig. 14 Axial stress vs time for different constant friction coefficients for the model of wing crack interaction

随着摩擦系数的增加，轴向应力亦有明显增加，而试件的破坏时间会延迟。相对于摩擦系数为常数0.7 的结果，MFL 破坏出现得更早(分别为113×10−6s 和123×10−6s)，破坏时试样轴向应力更低(分别为63.8 MPa 和91.6 MPa)。这是由于摩擦系数增大后，裂纹面需更高应力克服摩擦进行滑动，进而影响有效模量和试件的强度。此外，当初始裂纹开始滑动后，MFL 里的摩擦系数随着速度的增加而降低，促使裂纹增长并降低试件的强度。在地震事件中，尤其是大地震，断层面摩擦滑移产生的热能会引起近1000 ℃的温度增长，导致断层接触面数毫米到数厘米的融化[40]，断层面的摩擦系数急剧下降。因此，在研究动力加载过程中岩石材料的摩擦行为时，MFL 相对于静摩擦系数更为合理，可以考虑接触面滑动引起的摩擦弱化现象。

3 结论

基于岩石类材料的I 型裂纹模型，考虑正方形裂纹阵列情况下的裂纹相互作用，采用MFL 模型考虑裂纹面的动力摩擦过程，提出了一种考虑裂纹相互作用、裂纹密度以及裂纹面滑动摩擦的脆性材料动力模型。主要结论如下：

(1) 根据此模型和数值计算，定量研究了MFL模型以及常摩擦系数条件下，裂纹密度对裂纹扩展过程、试件受力和破坏的影响。结果表明，随着裂纹密度的增大，裂纹的相互作用增强，试件破坏时的加载应力降低，而惯性效应引起的试件轴向附加应力增大。

(2) 不同裂纹摩擦系数的结果表明，随着摩擦系数的增大，试样轴向应力增加，而破坏时刻会延迟。采用MFL 模型计算裂纹的动力摩擦过程时，假定在裂纹滑动前的摩擦系数为0.7。MFL 模型结果表明，随着裂纹开始滑动，裂纹面摩擦系数逐渐降低。相对于摩擦系数为常数0.7 的结果，其破坏出现得更早(分别为113×10−6s 和123×10−6s)，破坏时试件轴向应力更低(分别为6.38×107Pa 和9.16×107Pa)。动力摩擦过程会促进岩石材料中的裂纹发展，并降低材料强度。相对于常数摩擦系数，考虑速度及状态依赖型摩擦模型更为合理。

(3) 本模型计算的动强度因子与其他实验结果对比显示，脆性材料在动力加载过程中表现出明显的应变率效应和动力尺寸效应。本文提出的模型有助于进一步研究动力加载条件下，脆性材料的动力破坏过程。