炸药内爆是岩石及混凝土等材料受载破坏的常见形式，其破坏机理为是：在爆炸产生的应力波和高温高压气体共同作用下，介质内部出现裂纹并向自由面扩展，将介质切割成块。但对应力波和高温高压气体的作用分工研究较少，尤其是两者共同作用下的裂纹扩展机理尚不完全清楚。

裂纹扩展过程作为断裂力学的重要内容，一直是学者们研究的重点，测量方法也很多，主要包括电测法、模拟法、红外热成像法和光测法等。电测法主要使用电阻式应变片^[1]和裂纹扩展计（crack propagation gauge, CPG）^[2]等传感器，以金属栅丝断裂时刻作为裂纹发展到此的时间，从而得到裂纹的起裂、传播、止裂时刻。模拟法则主要采用各种动态分析软件，如RFPA-dynamic^[3-4]、LS-DYNA^[5-6]、AUTODYN^[7]、ABAQUS^[8]等进行数值试验，设定爆炸条件，模拟介质的动力响应，由于模拟法具有不受试验条件限制，易设定爆破参数等优势，应用较多。相关研究^[9]表明，岩石材料受到冲击致裂时，伴随体积膨胀和断裂，试件表面会发生能量耗散和温度跃升，由温度变化可得到裂纹发展过程。然而，由于缺乏相关理论研究，基于热成像的方法研究较少。光测法包括激光裂纹成像传感器^[10]、动态焦散线和数字图像相关等方法。杨立云等^[11]、杨仁树等^[12-13]用固体激光器取代多火花式点光源，用数码高速相机代替胶片相机，设计了新型数字激光动态焦散线测试系统，研究了多种爆破模式下裂纹扩展及相互作用过程。数字图像相关（digital image correlation，DIC）具有非接触、全场变形及应变测量等优势，在高速动态过程研究中应用较为广泛。Siviour等^[14]，Zhang等^[15]、申海艇等^[16]将DIC方法应用到SHPB试验中，开展了基于DIC技术的冲击加载下岩石类材料动态响应研究，同时测得了裂纹扩展速度；徐振洋等^[17]应用两台高速相机的DIC-3D方法研究了炸药外爆的聚能射流侵彻过程；杨立云等^[18]应用超高速摄影和DIC技术分析了聚碳酸酯平板内爆炸应力波的传播过程。模拟法虽然不受试验条件限制，可以迅速开展研究，但是在模拟分析过程中，需要对边界条件和材料属性进行简化或近似，难以真实还原试验条件；应变片及裂纹扩展计受载后会破坏，不能测得裂纹扩展全过程；近年来，随着光学成像技术以及图像处理技术的进步，基于光学传感器的光测法被越来越多地应用到爆破裂纹扩展过程研究当中。

本文中应用DIC和数字图像处理（digital image processing, DIP）技术，分析混凝土内爆试验中表面裂纹在长度和宽度方向的扩展过程及表面变形场和应变场的演化过程，讨论应力波和气体对裂纹扩展的分工作用，并分析断裂过程区的发展过程。

1.   裂纹扩展测量方法

1.1   基于数字图像处理技术的裂纹检测算法

高清、高速相机的出现，为研究爆炸、冲击等瞬态过程提供了有力工具，也使数字图像处理技术得到了全面发展。采用千眼狼5KF20高速摄影机、Nikon AF70-200 mm F/2.8镜头采集图像，远程快门触发，拍摄的裂纹图像如图1(a)所示，应用背景减除、图像增强、边缘检测等处理方法，提出了一种新的混凝土表面裂纹检测算法，裂纹检测流程如下。

图  1  裂纹检测过程

Figure  1.  Crack detection process

下载: 全尺寸图片 幻灯片

（1）背景减除

将起爆前拍摄图像作为背景，裂纹发展过程图像作为前景，由于散斑会干扰裂纹边缘检测，因此不能直接进行边缘检测，将前景减除背景可消除散斑影响，如图1(b)所示。

（2）图像增强

采用室外自然光源照明，且帧率较高，图像整体亮度较低，经过背景减除后，虽获得了裂纹对象，但其灰度特征不明显，采用伽马变换增加图像对比度，其基本形式为：

式中： s为输出灰度值，α为灰度缩放系数，r为输入灰度值，γ为灰度输入输出曲线。当γ＜1时，拉伸灰度级较低的区域，压缩灰度级较高的区域；当γ＞1时，则拉伸灰度级较高的区域，压缩灰度级较低的区域。本文中γ取值为0.6，增强后的裂纹特征图像如图1(c)所示。

（3）裂纹边缘检测

边缘检测的本质是对图像进行滤波，通常的做法是求图像一阶或二阶导数以突出图像局部变化最明显的部分，即边缘。求一阶导数常用Sobel算子、Canny算子等，求二阶导数常用Laplacian算子，本文中使用Canny算子检测到的裂纹边缘如图1(d)所示。

（4）临近裂纹连接

将图1(d)中检测到的相近裂纹进行连通，同时删除非裂纹的噪点区域，方法为：对提取的边缘进行腐蚀、开闭等形态学操作，连通临近裂纹，并删除孤立对象，如图1(d)中孤点为一爆破碎石。

（5）提取裂纹骨架

进行腐蚀操作后可得到的裂纹对象如图1(e)所示，以裂纹对象的中心像素为骨架，剔除两侧像素，将裂纹细化为单层像素，采用Zhang等^[19]提出的算法：以某一明点P_i,j为中心，考察其8邻域$\begin{array}{*{20}{c}} {{P_{i - 1,j - 1}}}&{{P_{i - 1,j}}}&{{P_{i - 1,j+1}}} \\ {{P_{i,j - 1}}}&{{P_{i,j}}}&{{P_{i,j + 1}}} \\ {{P_{i + 1,j-1}}}&{{P_{i + 1,j}}}&{{P_{i + 1,j + 1}}} \end{array}$内明点个数及明暗点组合方式，满足删除条件时，则删除该点，每一轮迭代删除一层像素，直至目标只剩一层像素，如图1(f)所示。由裂纹骨架中明点个数得到裂纹的长度，单位为像素，再由图中标尺所代表的实际物理距离，可得到每个像素的物理尺寸，计算裂纹的物理长度。

设在时刻t_i (i=1, 2, ···, n)，第i帧图像中检测的裂纹为$c_i^k\,(k = 1,2, \cdots ,m)$，k为裂纹数，裂纹扩展速度可用下式计算：

式中：$v_i^k$为第k条裂纹扩展速度，$n_i^k$和$n_{i - 1}^k$分别为第k条裂纹在第i帧和第i−1帧中明点累计数，D_p为单个像素点所代表的物理距离，本试验中为0.89 mm/pixel。

1.2   基于数字图像相关方法的全场变形测量

数字图像相关方法的本质是模板匹配。如果将受载变形前后图像分别定义为参考图像I_R和匹配图像I_M，参考图像中某点P₀(x₀,y₀)变形后在x和y方向的位移分别为u和v，那么P₀在匹配图像中的对应点为P(x₀+u,y₀+v)。在参考图像中以P₀为中心建立大小为(2i+1)×(2j+1)的参考子区f(x,y)，在匹配图像中建立同样大小的匹配子区g(x+u,y+v)，使匹配子区在匹配图像中滑动，同时进行相关运算，计算相关系数C_f,g,u,v，当相关系数取得极值时，点P即为P₀的对应点，两点坐标差即为位移矢量(u,v)，此时位移为整像素精度，并采用曲面拟合法求解亚像素位移。采用零均值归一化互相关（zero-based normalized cross correlation, ZNCC）函数计算相关系数，表示如下：

式中：${f_{ij}} = \dfrac{1}{{(2i + 1)(2j + 1)}}\displaystyle\sum\limits_{x = - i}^i {\displaystyle\sum\limits_{y = - j}^j {f(x,y)} }$和${g_{ij}} = \dfrac{1}{{(2i + 1)(2j + 1)}}\displaystyle\sum\limits_{x = - i}^i {\displaystyle\sum\limits_{y = - j}^j {g(x + u,y + v)} }$分别为参考子区和匹配子区灰度均值。对图像中每一像素进行相关运算，可得到位移场分布情况，对位移场采用逐点最小二乘法求解应变场，再扩展到图像序列中每一帧，便可得到位移场和应变场演化过程。

2.   试　验

2.1   测试系统

测试系统由混凝土模型、爆破器材、图像采集处理设备及喷漆、冲击钻、发电机等辅助材料组成，如图2(a)所示。

图  2  混凝土表面裂纹扩展试验系统

Figure  2.  Experiment system for crack propagation of concrete surface

下载: 全尺寸图片 幻灯片

混凝土表面裂纹扩展过程以面内位移为主，故采用单台相机的二维DIC方法，调节三脚架使相机视轴垂直被测表面，高度与被测表面高度一致，拍摄的爆破过程图像如图2(b)所示。

2.2   装药参数

模型采用商用混凝土浇筑而成，尺寸为40 cm×40 cm×30 cm，木制框架，自然养护28 d，经测试其物理力学参数如表1所示。

表  1  混凝土模型物理力学参数

Table  1.  The physical and mechanical parameters of the concrete model

混凝土标号	容重/(kg·m−3)	抗压强度/MPa	弹性模量/GPa	抗拉强度/MPa	泊松比
C30	2.7×103	42.5	31.2	3.2	0.3

下载: 导出CSV 

| 显示表格

爆破器材为8#导爆管雷管加2 g黑索金，先将黑索金倒入内径10 mm、底部封堵的PVC管中，再插入导爆管雷管，压实，顶端用胶带缠绕，总装药长度约4 cm，装药量3 g。

使用冲击钻在模型顶面钻炮孔，孔径12 mm，孔深12 cm，爆破抵抗线约9.5 cm，炮孔距离被测面20 cm，如图2(c)所示。为保证填塞质量，先在孔口填入石英砂，再灌注502胶水固结。被测表面喷涂散斑并描画刻度尺以便于获得像素与物理距离之间的比例。拍摄帧率为20 000 s⁻¹，曝光时间为50 μs，拍摄时长为3 s，分辨率为448×484，远程快门线触发。

3.   结果分析

3.1   起裂点分析

拍摄的裂纹扩展过程如图3所示。起爆后0.50 ms模型侧面开始出现裂纹，起裂点位于模型上表面与侧面棱线处，如图3(b)所示；随后裂纹进一步扩展，于1.75 ms结束，如图3(c)～(e)所示；长度扩展结束后转而向宽度扩张，最终于4.50 ms结束，如图3(f)所示。

图  3  裂纹扩展过程

Figure  3.  Crack propagation progress

下载: 全尺寸图片 幻灯片

炸药在混凝土内部爆炸时，其作用过程分为三种模式^[20]：一是在应力波的作用下，混凝土内部产生微裂隙，并使混凝土本身由于缺陷存在的原生裂隙扩展；二是应力波到达自由面后反射,引起自由面附近发生“剥落”；三是爆炸生成的气体使原生裂隙和次生裂隙进一步发展，并将介质切割成块。基于该理论，并结合装药结构，最小抵抗线约9.5 cm，侧面抵抗线为20 cm，而棱线处抵抗线约22.1 cm，可知模型表面最初的裂纹应形成于炮孔孔口处。侧面起裂点没有出现在侧面而出现在棱线上，说明一旦形成贯通裂纹，该裂纹将成为主裂纹，或称为“优势”裂纹，“优势”裂纹在表面扩展速度极快，也将决定后续裂纹的趋向。

3.2   裂纹扩展过程

3.2.1   长度扩展

应用第1.1节介绍的裂纹检测算法依次对图像序列进行检测，得到裂纹长度扩展过程如图4所示。由图4可知，裂纹长度扩展从0.50 ms开始，到1.75 ms停止，期间最大速度225.95 m/s，平均速度122.27 m/s。依速度值将曲线划为2段，0.50～1.20 ms为第1段，裂纹扩展速度始终大于平均值，起裂后速度跃升快，极值出现在0.60 ms，极值后有一定波动，整体稳定于较高水平，可称为裂纹稳定扩展段，此间裂纹长度为133.15 mm，达到总长的83%；第2段为1.20～1.75 ms，裂纹扩展速度急剧衰减，推测此时应力波作用已十分微弱，甚至已停止，只有爆炸气体继续推动裂纹扩展，最终裂纹长度扩展至159.92 mm停止。黄家蓉等^[21]、冯盼学等^[22]测得内爆条件下应力波在混凝土中传播的速度达10³ m/s级，作用时间为几十至几百微秒，而侧面裂纹起裂于0.5 ms，止裂于1.75 ms，可见，裂纹在长度上的扩展为应力波和爆炸气体共同作用的结果。

图  4  裂纹长度和扩展速度曲线

Figure  4.  Curves of crack length and propagation velocity

下载: 全尺寸图片 幻灯片

3.2.2   表面位移场及宽度扩展过程

应用第1.2节介绍的数字图像相关方法计算裂纹区域水平位移场,如图5所示，横纵坐标分别为裂纹区域像素坐标，位移单位为mm。

图  5  水平位移场的演化过程

Figure  5.  Evolution of horizontal displacement field

下载: 全尺寸图片 幻灯片

由图5可知，裂纹的扩展过程是被测表面位移场由无序到有序不断演化、并在裂纹区域不断集中的过程。起爆后0.25 ms，表面位移普遍较小，左侧为负值，极值为−0.041 mm，右侧为正值，极值为0.052 mm，宏观裂纹尚未形成，裂纹区域位移为0～0.025 mm；0.50 ms时起裂点两侧位移增大，达到−0.052～0.126 mm；0.75～1.75 ms为裂纹两侧位移逐渐增长、裂纹长度逐渐延伸的过程，长度扩展结束时，裂纹右侧位移极值为0.318 mm，左侧为−0.282 mm；1.75～4.50 ms裂纹长度基本不再增加，而转向横向（宽度）发展，裂纹宽度扩展结束时，裂纹两侧位移极值分别为1.067 mm和−0.523 mm。

传统基于应变片、焦散线和裂纹扩展计（CPG）等的方法对裂纹长度方向扩展研究较多，而对裂纹在宽度方向的扩展研究较少，为方便描述，本文中将裂纹长度方向称为扩展，宽度方向称为扩张。

在起裂点左侧和右侧每隔10 pixel (8.9 mm)分别取5对监测点，记为C_Li (i=1, 2, ···, 5)和C_Ri (i=1, 2, ···, 5)，则起裂点裂纹张开位移可用各监测点相对位移表示：

以时间为横轴，以裂纹张开宽度W_i为纵轴，裂纹扩张曲线如图6(a)所示。从起裂点开始，沿着裂纹扩展方向，以5 pixel (4.45 mm)为步长，以裂纹为中心，按照上述方法分别取监测点，得到沿扩展方向的裂纹宽度曲线，如图6(b)所示。

图  6  裂纹张开宽度随时间和裂纹长度的变化曲线

Figure  6.  Curves of crack width with time and crack length

下载: 全尺寸图片 幻灯片

由图6(a)可知，无论监测点距离起裂点多远，其相对位移发展趋势是一致的，裂纹扩展停止后，横向扩张要持续更长时间，可达4.5 ms，在裂纹发展的整个过程中，爆炸气体均起到重要作用，但其对裂纹的扩张作用更明显。如前文所述，裂纹扩展速度于1.20 ms明显下降，这是由于应力波的作用微弱或已结束，但此时爆炸气体的作用仍较强，其作用更主要表现为使裂纹进一步张开，且能维持裂纹稳定扩张。随着距裂纹中心距离的增大，各监测点间的相对位移增大，监测点W₁最大相对位移1.01 mm，W₂最大相对位移1.36 mm，W₃、W₄、W₅最大相对位移较为接近，分别为1.59、1.55、1.54 mm。说明靠近裂纹近区介质内部受应力波拉伸作用较强，表现为因孔隙、微裂隙等增大而变得“疏松”，此“疏松”区宽度可达5.34 cm，由图5中位移最大区域并非裂纹附近而是与裂纹有一段距离，同样可以证实这种“疏松”区的存在。从第3监测点开始远离裂纹区域，介质更多地表现为刚体位移，因此，以W₃、W₄、W₅表示裂纹张开宽度是合适的。由图6(b)可知，随着裂纹长度的增加，其张开宽度逐渐减小，起裂点附近宽度最大，可达1.59 mm，止点附近张开宽度最小，仅为0.4 mm。

3.3   应变场演化过程

最大主应变场的演化过程如图7所示，其中，红色表示应变较大，由图7可知，裂纹区域始终存在应变集中带。在0.25 ms时，宏观裂纹尚未产生，但应变已向裂纹区域集中；0.50 ms时，未来的裂纹附近应变集中带已经形成，起裂点附近最明显，随着宏观裂纹的产生和扩展，应变集中带越发明显，1.75 ms时，应变集中带的总体形态已经确定；在1.75～4.50 ms，在裂纹附近应变集中带上，深色区域不断扩大，宏观表现为此处裂纹正在扩张。

图  7  最大主应变场的演化过程

Figure  7.  Evolutionary process of maximum principle strain field

下载: 全尺寸图片 幻灯片

3.4   断裂过程区分析

断裂过程区（fracture process zone, FPZ）模型是Hillerborg等^[23]在研究混凝土断裂过程中提出的，该模型认为混凝土裂纹的扩展以裂纹前端形成的微裂区为先导，此区域内材料强度降低，但由于骨料和砂浆之间的互锁效应，仍能传递部分应力，称为黏聚力，其作用为抵抗裂纹张开，并使裂纹转向。

FPZ的观测通常在带有预制裂纹混凝土梁的三点弯曲试验中进行，属于准静态加载，对于动态加载过程FPZ的研究则较少。FPZ尺寸的测量也逐步发展为以DIC方法为主，结合声发射（acoustic emission，AE）^[24]测量、引伸计^[25]等。FPZ区域的界定不尽相同，有的以应变作为界定标准，也有的以位移作为标准^[26]。本文中以各帧中宏观裂纹尖端为FPZ起点，应变集中带作为FPZ边界，边界应变值取极限拉伸应变100×10⁻⁶（文献[27]中测得相似配比混凝土极限拉伸应变为(80～120)×10⁻⁶），确定FPZ的尺寸，如图8所示，图中横纵坐标的单位为像素，裂纹尖端和FPZ尖端坐标也为像素，0.80 ms后FPZ已超出观测范围，结合图7，计算各时刻FPZ尺寸如表2所示。

表  2  断裂过程区尺寸

Table  2.  Sizes of fracture process zone

时间/ms	FPZ宽度/mm	FPZ长度/mm
0.50	33.82	132.64
0.75	47.17	129.30
0.80	48.95	132.16
1.00	52.51	−
2.00	50.73	−
3.00	49.84	−
4.00	51.62	−

下载: 导出CSV 

| 显示表格

图  8  断裂过程区发展过程

Figure  8.  Propagation of fracture zone

下载: 全尺寸图片 幻灯片

由表2可知，除起裂时刻FPZ宽度较小，为33.82 mm外，其余时刻宽度变化不大，在47.17～52.51 mm之间，FPZ长度基本不变，为129.30～132.64 mm，约为混凝土骨料最大粒径（15 mm）的8～9倍。赵艳华等^[28]的研究表明，断裂过程区的尺寸与混凝土的配比、骨料粒径等因素有关，为混凝土材料固有性质，但也有研究表明，断裂过程区还表现出与混凝土试件大小相关的尺寸效应。

4.   结　论

结合数字图像相关方法和图像处理等技术，研究了内爆加载下混凝土表面裂纹扩展规律。该技术可观测裂纹的扩展和扩张效应以及断裂过程区的演化规律，同时提供裂纹变化的量化数据，是爆炸裂纹研究的有效方法。通过进行混凝土内爆试验，得到主要结论如下：

（1）侧面裂纹起裂点位于侧面与顶面交线上，扩展平均速度为122.27 m/s，最大速度为225.95 m/s，裂纹总长为159.92 mm，裂纹方向与炮孔轴线方向近似平行；

（2）裂纹长度扩展集中在起爆后0.50～1.75 ms内，为应力波和爆炸气体共同作用结果，扩展停止后，裂纹主要进行横向扩张，持续时间可达4.5 ms，主要动力为爆炸生成的气体楔入裂纹内，使模型被切割，最终裂纹最大张开宽度为1.59 mm；

（3）根据应变场分析了爆破这种超动态加载下混凝土断裂过程区的演化规律，混凝土配比和骨料粒径固定后，其断裂过程区尺寸基本不变，其长度约为骨料最大粒径的8～9倍；

（4）试验采用单药包集中装药，且装药量较小，侧面只产生了一条宏观主裂纹，采用多药包装药，并设计不同装药参数和布置形式，对于全面分析多条裂纹扩展过程大有裨益；另外，电测法擅长对应力波的测量，光测法对裂纹动态扩展过程的测量更方便，将二者结合起来更有助于研究爆炸加载下介质的响应和应力波传播规律，这将是下一步工作的方向。