本发明涉及实验力学数字图像处理技术领域，特别涉及一种非连续变形测量方法。
背景技术：
当在外力作用下，混凝土、岩石等材料未产生很大变形时，其表面通常会伴随有裂缝的产生，在断裂力学领域中，将具有上述破坏特征的材料称为准脆性材料。由于准脆性材料在产生裂缝后，其表面将不再保持为连续状态，因此将准脆性材料的这种破坏方式称为非连续变形破坏。与准脆性材料相对应，能够承受很大变形，且一般不会产生裂缝的材料，称为延性材料，如金属、机械等，延性材料产生的破坏称为连续变形破坏。
裂缝作为准脆性材料受力变形后所产生的一种外观表象，已成为广大科研工作者与工程技术人员关注的热点问题。由于准脆性材料的断裂破坏过程始终伴随有裂缝的产生，因此借助于一种可靠的试验测量手段，观测准脆性材料裂缝的扩展过程，准确获取裂缝扩展尺寸、裂缝扩展尖端位置以及该位置处的真实位移和应变信息，对于理解准脆性材料的断裂破坏机理，具有十分重要的科研意义与工程实用价值。
在实验力学领域中，数字图像相关方法以其非接触、全场测量且操作简便等优良特性，逐渐得到了推广与应用。黄才政等(黄才政.基于数字图像相关方法的混凝土断裂过程区力学特性研究[d].重庆交通大学,2016.)通过开展混凝土三点弯曲预制缺口梁断裂试验，结合数字图像相关方法，对裂缝扩展过程进行了研究。图2为其试件破坏前的图像及在不同荷载阶段下由数字图像相关方法计算得出的全场水平应变云图，其中(a)给出了该试验中试件最终破坏时的图形，由图可看出在预制裂尖前端有一条竖向扩展裂缝产生，(b)～(d)为不同荷载阶段下由数字图像相关方法计算得出的全场水平应变云图(图中pmax代表峰值荷载，(↑)代表荷载上升段，(↓)代表荷载下降段)，可看出在裂缝产生后，裂缝扩展位置处仍存在应变信息，未能将裂缝区域进行剥离。吴智敏等(wuzm,rongh,zhengjj,etal.anexperimentalinvestigationonthefpzpropertiesinconcreteusingdigitalimagecorrelationtechnique[j].engineeringfracturemechanics,2011,78(17):2978-2990.)通过开展混凝土三点弯曲梁断裂试验，采用数字图像相关方法，对混凝土断裂过程区进行了研究。图3为试件破坏前的图像及在不同荷载阶段下由数字图像相关方法计算得出的全场水平应变云图，其中(a)给出了试件最终破坏时的图形，由图可知该试件在最终破坏时有一条竖向裂缝产生，(b)～(c)为不同荷载阶段下由数字图像相关方法计算得出的全场水平应变云图，可看出在裂缝产生后，裂缝扩展位置处仍存在应变信息，未能将裂缝区域进行剥离。
由以上两个实例可看出，现有的数字图像相关方法，在处理准脆性材料的非连续变形破坏问题时，存在以下不足之处：①无法自动检测裂缝扩展边界，获取裂缝扩展尺寸信息；②无法对裂缝扩展区域进行剥离，获取裂缝扩展尖端所在位置以及该位置处的真实位移和应变信息。因此有必要结合现有的数字图像处理技术，提出一种非连续变形测量方法，解决准脆性材料的非连续变形测量问题，服务于科学研究与实际工程领域。
技术实现要素：
针对现有数字图像相关方法的不足，本发明拟解决的技术问题是，针对准脆性材料的非连续变形问题，提供一种非连续变形测量方法，该方法基于目前已有的数字图像处理技术，能够对裂缝扩展区域进行边界检测与剥离，获取裂缝扩展尺寸、裂缝扩展尖端所在位置以及该位置处的真实位移和应变信息。
本发明解决所述技术问题采用的技术方案是：提供一种非连续变形测量方法，该测量方法的步骤如下：
第一步：在被测物体表面绘制人工散斑图案；
第二步：安装装置；将两部平行线光源固定于固定支座上，根据周围环境光线明暗程度手动调节光照强度；将采集相机固定于三角支架上，同时采集相机通过导线与计算机相连，安装好所用装置；
第三步：开动试验机，对被测物体进行加载实验，按照预先设定的采集帧率进行图像采集；
第四步：根据实验荷载数据信息，均匀挑选图像，进行扩展裂缝边界检测，并标识出扩展裂缝的边界，得到裂缝扩展区域，根据扩展裂缝边界的灰度坐标信息，获取裂缝扩展尺寸；若没有检测到扩展裂缝，则通过数字图像处理得到图像的全场位移和应变信息；
第五步：对第四步中已经检测出的裂缝扩展区域进行剥离，获取裂缝剥离之后的图像：即将裂缝扩展区域内的灰度数据全部转化为0或255，以与图像中其他非裂缝扩展区域进行区分，并突显裂缝剥离效果，得到裂缝剥离后的图像；
第六步：对裂缝剥离后的图像采用数字图像处理方法进行分析，获取裂缝扩展尖端所在位置以及该位置处的位移和应变信息，进而得到剥离裂缝扩展区域后的全场位移和应变信息。
与现有技术相比，本发明的有益效果为：
与现有数字图像相关方法相比，该方法的特色在于能够解决混凝土、岩石等准脆性材料的非连续变形测量问题，获取裂缝扩展尺寸、裂缝扩展尖端所在位置以及该位置处的真实位移和应变信息，能够在科学研究领域与实际工程应用中得到广泛应用。
附图说明
图1为本发明非连续变形测量方法所用装置立面示意图；
图2为背景文件1中试件破坏前的图像及在不同荷载阶段下由数字图像相关方法计算得出的全场水平应变云图；
图3为背景文件2中试件破坏前的图像及在不同荷载阶段下由数字图像相关方法计算得出的全场水平应变云图；
图4(a)为混凝土沿竖直方向产生了一条贯穿型裂缝的图像；
图4(b)为图4(a)的图像经本申请方法进行边界检测得到的图像；
图4(c)为图4(b)的图像经本申请方法进行裂缝剥离后得到的图像；
图5(a)为混凝土产生了多条裂缝的图像；
图5(b)为图5(a)的图像经本申请方法进行边界检测得到的图像；
图5(c)为图5(b)的图像经本申请方法进行裂缝剥离后得到的图像；
图6为实施例所用装置的实物图；
图7(a)为实施例中试件未受外力时拍摄的图像；
图7(b)为实施例中试件在峰值荷载时刻对应的带扩展裂缝的图像；
图7(c)为采用本发明非连续变形测量方法计算得到的全场水平位移云图；
图7(d)为采用现有数字图像方法计算得到的全场水平位移云图；
图中，人工散斑1，被测物体2，线光源3，采集相机4，数据传输导线5，计算机6，裂缝识别模块7、图像处理模块8、三角支架9，隔震垫10，固定支座11，固定支架12，卡扣13。
具体实施方式
为了使本发明的目的及技术方案更加清晰，以下结合附图及实施例，对本发明做进一步详细介绍。本发明的具体实施方式并不限于此，任何本领域的普通技术人员在未付出实质性创新前提下所做出的同等修改与完善，均应在本发明的保护范围之内。
本发明提供一种非连续变形测量方法，该测量方法的具体步骤如下：
第一步：在被测物体表面绘制人工散斑图案；
第二步：安装装置；将两部平行线光源固定于固定支座11上，根据周围环境光线明暗程度手动调节光照强度；将采集相机4固定于三角支架9上，同时采集相机4通过导线5与计算机6相连，安装好所用装置；
第三步：开动试验机，对被测物体进行加载实验，按照预先设定的采集帧率进行图像采集；
第四步：根据实验荷载数据信息，均匀挑选图像，假设图像的灰度值在0～255之间，进行扩展裂缝边界检测，并标识出扩展裂缝的边界，得到裂缝扩展区域，根据扩展裂缝边界的灰度坐标信息，获取裂缝扩展尺寸；若没有检测到扩展裂缝则通过数字图像处理得到图像的全场位移和应变信息；
第五步：对第四步中已经检测出的裂缝扩展区域进行剥离，获取裂缝剥离之后的图像：即将裂缝扩展区域内的灰度数据全部转化为0或255，以与图像中其他非裂缝扩展区域进行区分，并显示裂缝剥离效果，得到裂缝剥离后的图像；
第六步：对裂缝剥离后的图像采用数字图像处理方法进行分析，获取裂缝扩展尖端所在位置以及该位置处的位移和应变信息，进而得到剥离裂缝扩展区域后的全场位移和应变信息。
上述非连续变形测量方法所用装置(参见图1)包括用于加载被测物体的试验机、线光源3、采集相机4、用于对图像中裂缝扩展区域进行边缘检测与剥离并获取裂缝扩展尺寸信息的裂缝识别模块7以及用于接收裂缝识别模块7传输来的去裂缝扩展区域的图像数据并进行图像处理分析获取全场位移和应变信息的图像处理模块8；被测物体上涂有人工散斑1，所述线光源3采用两部平行线光源，光源的两端垂直设置固定支架13，固定支架13固定在固定支座11上，在固定支架朝向光源的内侧上从上至下设有若干数量的卡扣12，线光源通过卡扣12能固定在固定支架的任意位置上，便于手动调节两部线光源之间的距离，以达到最佳的补光效果，光照强度可根据周围环境光线明暗程度手动进行调节；所述采集相机4固定于三角支架9上，三角支架底部设有隔震垫，用于消除地面震动对相机造成的干扰；所述采集相机4通过数据传输导线5与计算机相连，实时传输图像数据；计算机6内加载有裂缝识别模块7和图像处理模块8，可以接收通过数据传输导线5传来的图像数据，同时对接收端的图像数据进行分析处理。
所述裂缝识别模块7用于对图像中裂缝扩展区域进行边缘检测与剥离，进而获取裂缝扩展尺寸信息。若材料未有裂缝产生，则直接送入到图像处理模块8中。
所述图像处理模块8用于接收裂缝识别模块7传输来的去裂缝扩展区域的图像数据，使用现有数字图像技术对图像进行处理分析，获取全场位移和应变信息。
所述人工散斑1采用哑光黑与哑光白两种喷漆绘制而成。首先均匀喷涂一层哑光白漆作为底色，再反复喷洒哑光黑漆形成高反差散斑图案。
所述线光源3采用两部平行线光源，光源放于固定支架13上，光照强度和光照间距可根据周围光线条件手动进行调节。
所述三角支架9底部设有隔震垫10，用于消除地面震动对采集相机4的影响。
所述采集相机4用于采集试验加载过程前后的散斑图像，应根据试验需要选用静态采集相机或者高速采集相机，精准捕捉被测物体裂缝开裂瞬间。若试验被测物体2从开始加载到破坏比较突然，则宜选用高速相机；若破坏过程比较缓慢，则使用静态采集相机即可满足。
本发明对图像中的裂缝区域进行检测剥离时，采用了目前常用的裂缝识别处理算法，其中包括直方图均衡化、中值滤波、图像二值化操作，以下分别对上述操作步骤进行描述：
1)直方图均衡化
受周围环境因素(如大气、光照等)的影响，在试验现场采集到的图像有可能存在整体或局部偏暗、过亮等现象，影响裂缝识别检测效果。通过对图像进行直方图均衡化处理可以增强图像的整体对比度，突显裂缝扩展区域，提高裂缝检测效果。
2)中值滤波
采用中值滤波方法，可以消除图像采集和传输过程中因成像设备与传输介质干扰而产生的噪声，有效保护图像的边缘细节信息，该方法可通过对像素邻域向量化取中值的方法来实现。
3)图像二值化
图像二值化是指通过设定一个灰度阈值t来分割目标与背景，大于阈值t的像素属于目标即记为1，其他则属于背景即记为0。假设一幅灰度裂缝图像用f(x,y)表示，其中(x,y)为图像中像素的坐标位置，若t为阈值，则阈值分割后的二值图像b(x,y)满足以下关系
在本非连续变形测量方法中，采用迭代优化的方法来计算阈值t，以提高裂缝边界的检测效果，其具体操作步骤如下：
1)设定阈值
统计带裂缝图像所有像素点位置的最小灰度值tmin、最大灰度值tmax，计算两者平均值作为初始阈值
2)阈值分割
根据步骤1)得到的初始阈值t，采用公式(1)对图像进行分割，按照式(2)得到目标像素集合g1和背景像素集合g2
g1＝{f(x,y)≥t},g2＝{f(x,y)＜t}(2)
3)均值
按照式(3)计算目标像素集合g1和背景像素集合g2的灰度平均值t1和t2，
其中，num(g1)表示目标像素集合g1中的像素点的总数量；num(g2)表示目标像素集合g2中的像素点的总数量。
4)迭代
根据t1和t2计算新的阈值
图4(a)所示混凝土沿竖直方向产生了一条贯穿型裂缝，图5(a)中混凝土产生了多条裂缝，将原始带裂缝图像传输给裂缝识别模块，其裂缝边界检测效果如图4(b)、图5(b)所示，该模块会根据边界检测结果自动在裂缝周围生成两条边界线，此时可以获得两条边界线上所有像素点的坐标信息，根据对应点的坐标信息可以获得任意位置处的裂缝扩展尺寸信息；图4(c)、图5(c)给出了裂缝的剥离效果图(将两条边界线区域内的所有像素点的灰度值设定为255)。
实施例
本实施例的实验在河北省土木工程技术研究中心开展，试验对象为混凝土三点弯曲梁试件，试验目的是为了获取任意荷载阶段下混凝土裂缝的扩展尺寸、裂缝扩展尖端所在位置以及该位置处的位移和应变信息。
图6为本实施例测量方法所用装置的实物图，其中i为电液伺服压力试验机，用于提供试验所需压力，使三点弯曲梁试件ii产生劈裂破坏；iv为采集相机，用于对试验全过程进行拍摄；iii为人工随机散斑图案，散斑作为变形信息载体，与三点弯曲梁试件保持为同步变形状态。图7(a)为混凝土试件未受外力时拍摄的图像，图7(b)为峰值荷载时刻对应的一张带扩展裂缝的图像，可看到在该时刻下产生了一条比较明显的竖向扩展裂缝。图7(c)和图7(d)给出了采用非连续变形测量方法与现有数字图像方法计算得到的全场水平位移云图(单位：mm)，可以看出，当有扩展裂缝产生时，该非连续变形测量方法能够对裂缝扩展边界进行检测剥离，根据边界像素点坐标，可以获取任意位置处的裂缝扩展尺寸信息，同时获取裂缝扩展尖端位置以及该位置处的位移和应变信息，而现有数字图像处理方法无法实现以上所述功能。
本发明未述及之处适用于现有技术。