1. 引言
近年来随着复合材料在航空航天、军事等领域的广泛应用,对复合材料韧性预测方法的研究越来越重要。在复合材料的韧性研究中,通过断裂力学方法研究复合材料分层断裂的韧性是一种有效方法。在使用过程中复合材料的刚度、强度性能会随着断裂裂纹的扩展而下降。复合材料I型层间断裂韧性的测量方法通常为单向复合材料的末端切口试样的双悬臂梁(Double cantilever beam, DCB)试验 [1] 。姚辽军通过对复合材料分层扩展路径和分层扩展断面进行显微镜观察和SEM电镜扫描分析,明确了复合材料I型静态分层扩展中的典型断面断裂机理 [2] 。张璐以飞机实际构件中的典型分层缺陷为研究对象,预测分层在压缩载荷及疲劳载荷作用下的扩展情况,定量分析分层缺陷对复合材料结构的影响,并通过双柔度法测试复合材料Ⅰ型层间断裂韧性 [3] 。
在进行DCB实验时,为了保证载荷、位移实时对应,实验人员通常在试验过程中每次发生裂纹扩展时中途暂停,再测量此刻的裂纹长度。该方法虽然能够测量出裂纹长度,但由于实验过程中多次暂停,导致实验过程中的裂纹扩展与真实场景存在一些差异,从而造成测量误差 [4] 。因此,为了提高DCB实验过程中的准确性和效率,提出一种基于图像分割的复合材料I型层间断裂裂纹识别方法。首先使用Sobel边缘检测算子提取复合材料截面的刻度值信息,然后使用图像分割模型检测复合材料断裂裂纹区域,最后通过裂纹区域长度与刻度值信息计算出断裂裂纹的真实长度值。
2. 裂纹识别方法原理
2.1. Sobel边缘检测算子
Sobel算子是一种基于图像强度的一阶微分算子,通常用于图像边缘检测 [5] 。与Roberts算子、Laplacian算子和Canny算子 [6] 等经典边缘检测算子相比,Sobel算子具有计算方法简单、处理速度快等优点。传统的Sobel算子边缘检测是利用水平方向(0˚)和垂直方向(90˚)两个方向的边缘检测模板与图像作卷积操作,通过计算图像的一阶导数得到图像的边缘点像素值。在实际使用中Sobel算子通常选取3 * 3尺寸的模板,0˚和90˚的3 * 3尺寸Sobel算子模板S0和S90分别为
假设连续的图像用函数
,其在
处的梯度是一个具有方向和大小的矢量,即
(1)
式中
分别是
方向的单位矢量。
在
处的梯度幅值表示为
(2)
2.2. 图像分割算法
图像语义分割是基于标签数据的像素到像素的分类,将像素(Pixel)按照图像中不同的语义含义进行分割(Segmentation),也称图像语义标注(Image semantic labeling) [7] ,具体实例如图1所示。
(a) 原图
(b) 语义分割图
Figure 1. Image semantic segmentation diagram
图1. 图像语义分割图
Deeplab语义分割系列是Google团队提出的语义分割网络 [8] ,Deeplab v3的网络结构如图2所示,它的网络结构主要分为编码部分和解码部分。在编码部分,原始的图像首先通过以深度卷积网络(Deep Convolutional Network, DCNN)为主干的特征提取网络提取特征信息,从上到下网络越来越深,其中浅层网络卷积的感受野较小,主要用于提取图像中的空间信息,如图像的纹理、轮廓等等;深层网络卷积的感受野大且融合了浅层的特征,主要用于提取图像的高级语义信息。在解码部分,通过上采样的方式还原特征图的尺寸,为了融合网络结构的上下文信息,将编码阶段不同深度的输出拼接到解码部分对应尺寸的特征图上。
空间金字塔池化(Atrous Spatial Pyramid Pooling, ASPP)是DeepLab中用于语义分割的一个模块。由于图像中待检测物体具有不同的尺度,给分割增加了难度。因此,DeepLab设计了一个ASPP的模块,既能得到图像的多尺度信息,同时计算量又比较小。图3即是空间金字塔池化结构,通过使用不同rate大小的空洞卷积(Atrous Convolution)使得每个分支的尺度不同,最后将不同分支的尺度信息进行融合。
![](//html.hanspub.org/file/48-2570871x18_hanspub.png?20230519091716970)
Figure 2. Deeplab V3 image segmentation network structure
图2. Deeplab V3图像分割网络结构
3. 裂纹识别模型构建
3.1. 准备数据集
使用相机拍摄复合材料DCB试验过程,采集复合材料I型层间断裂裂纹图片,图片的分辨率为1920 * 1080,图片的样例图如图4所示。
![](//html.hanspub.org/file/48-2570871x20_hanspub.png?20230519091716970)
Figure 4. Interlayer fracture crack diagram of Type I composite material
图4. 复合材料I型层间断裂裂纹图
3.2. 标注数据
共采集485张断裂裂纹图片,然后使用labelme标注工具对图片进行标注,图像中需要标注的类别有横截面区域和断裂裂纹区域两个类别。标注好的图片如图5所示,其中红色区域为复合材料的横截面,绿色区域为复合材料层间断裂裂纹。
![](//html.hanspub.org/file/48-2570871x21_hanspub.png?20230519091716970)
Figure 5. Effect drawing of fracture crack picture annotation
图5. 断裂裂纹图片标注效果图
3.3. 训练裂纹分割模型
我们将标注好的数据分为训练集和测试集,其中训练集样本数量为385张,测试集数据为100张,训练集和测试集数据随机分配。训练模型时需要设置下面的部分参数,图像输入尺寸inputSize = 512 * 512,初始学习率lr_start = 1e−3,衰减权重weight_decay = 1e−4,最大训练次数max_iter = 40000。分割模型的训练过程如图6所示。
随着模型的训练迭代次数的增加,模型的loss误差越来也越小,loss值从最开始的5.3下降到0.69。训练40,000轮迭代后loss收敛得到最终的模型。训练好的分割模型在测试集中的结果为mIOU = 0.946。
![](//html.hanspub.org/file/48-2570871x22_hanspub.png?20230519091716970)
Figure 6. Training process of crack segmentation model
图6. 裂纹分割模型训练过程
4. 复合材料层间I型层间断裂裂纹识别实验
4.1. 横截面刻度信息提取
图7是复合材料横截面区域,在DCB实验过程中我们需要提取截面区域的刻度信息,方便后续过程中对断裂裂纹的长度进行识别。首先通过图像分割模型提取出复合材料横截面区域,由于横截面区域的颜色信息较多,在进行刻度边缘提取时需要对横截面区域图像进行灰度化处理。原始图像有R、G、B三个通道,通过公式(1)将其转换为单通道灰度图:
(3)
其中,R、G、B分别表示图像的红、绿、蓝通道,(x, y)表示图像像素坐标,f是得到的灰度图像。
![](//html.hanspub.org/file/48-2570871x24_hanspub.png?20230519091716970)
Figure 7. Schematic diagram of cross section of composite material
图7. 复合材料横截面示意图
对复合材料横截面区域图像灰度化处理后,然后使用垂直方向的Sobel边缘检测算法提取该区域中刻度信息,图8是经过Sobel边缘检测算子处理后的图像,横截面的刻度信息被提取出来。整个提取流程如图9所示。
4.2. 断裂裂缝检测
在DCB实验过程中,我们通过训练好的裂纹分割模型对复合材料I型层间断裂裂纹进行实时分割,裂纹的分割结果如图10所示。分割模型能够准确的提取到断裂裂纹。
![](//html.hanspub.org/file/48-2570871x25_hanspub.png?20230519091716970)
Figure 8. Cross section region processed by Sobel edge detection operator
图8. Sobel边缘检测算子处理后的横截面区域
![](//html.hanspub.org/file/48-2570871x26_hanspub.png?20230519091716970)
Figure 9. Flow chart of cross section scale information extraction
图9. 横截面刻度信息提取流程图
![](//html.hanspub.org/file/48-2570871x27_hanspub.png?20230519091716970)
Figure 10. Segmentation effect of crack segmentation model
图10. 裂纹分割模型的分割效果图
4.3. 断裂裂纹长度识别分析
在DCB实验过程中,通过Sobel边缘算子提取到横截面的刻度信息,然后通过裂纹分割模型提取到断裂裂纹的区域,最后将断裂裂纹区域尺寸映射到实际的长度值。这样在DCB实验过程中就能实时的计算出复合材料I型层间断裂裂纹的长度,避免在实验时中途暂停去测量裂纹长度。图11是DCB实验过程中,图像分割方法自动识别出的裂纹长度与实际长度的比较结果。横坐标为DCB实验过程中的时间,纵坐标为断裂裂纹的长度,实验过程中随着时间变化裂纹长度逐渐变大。图像分割方法识别的裂纹长度与实际长度平均误差为2.9%,说明基于图像分割的复合材料I型层间断裂裂纹识别方法在DCB实验过程中能够有效的识别出裂纹长度。
![](//html.hanspub.org/file/48-2570871x28_hanspub.png?20230519091716970)
Figure 11. Comparison of image segmentation crack recognition results
图11. 图像分割裂纹识别结果对比
5. 结论
为了解决DCB实验过程中人工测量复合材料I型层间断裂裂纹长度带来的误差问题,本文提出一种基于图像分割的复合材料I型层间断裂裂纹识别方法。该方法通过Sobel边缘检测算子提取复合材料横截面的刻度信息,然后利用图像分割模型获取断裂裂纹的区域,最终识别出断裂裂纹的真实长度值。通过实验表明,本文提出的方法能够有效地识别出DCB实验过程中复合材料I型层间断裂裂纹的长度,提高了DCB实验过程中裂纹长度测量的准确率和效率,为后面的复合材料分层韧性预测提供了重要参数指标。