一种应用复合裂纹位移场函数的数字图像相关方法与流程

1.本发明涉及实验力学、非接触式全场位移测量、数字图像相关方法，具体涉及一种应用复合裂纹位移场函数的数字图像相关方法。


背景技术：

2.在航空航天、建筑桥梁等诸多领域中，保证结构可靠性永远是第一要务。随着科学技术的发展，能够在结构发生不可逆损伤之前，对其施以控制修缮，成为诸多科学家的重要目标。而由实验力学提供的材料响应载荷表征的测量结果，对于完善设计，改进结构以及监测损伤都具有丰富的指导意义。位移与应变，作为实验力学重要参量，对于分析模型力学特性、验证基础假设及特征识别等都具有举足轻重的作用。面对日益复杂的测量需求，传统的位移应变测量方法就显得捉襟见肘。与此同时，有限处位移应变测量结果远不能胜任复杂模型的分析工作。并且针对各种测量条件的约束，能够实现全场应变及非接触式测量的新兴测量技术应运而生。
3.能够进行全场测量的不同方法，各有其准确性及适用性优势。而又能够同时满足非接触测量条件的数字图像相关法（digital image correlation，dic），由于其较低的实验成本及数据采集过程相对简单等优势已被广泛用于实验力学中，其对测量环境和隔振要求较低，普适性强，广泛应用于各种工程和学科领域，如航空航天、生物医疗等。
4.传统的数字图像相关法（dic）中，子区的位移表征函数为一阶或二阶连续函数。而在断裂实验进行过程中，由于裂纹的产生，导致原本连续的参考子区会被裂纹切断，从而导致含裂纹区域的参考子区位移计算不准确甚至完全无法计算的情况。


技术实现要素：

5.本发明针对现有技术中存在的问题，公开了一种应用复合裂纹位移场函数的数字图像相关方法，该方法应用i、ii型复合裂纹位移场函数作为数字图像相关方法中计算子区的位移表征函数。本发明的方法能保证数字图像相关方法在处理含裂纹图像时的计算精度。
6.本发明是这样实现的：一种应用复合裂纹位移场函数的数字图像相关方法，其特征在于，所述的方法为使用i、ii型复合裂纹的位移场函数作为数字图像相关方法的子区位移表征函数，具体为：步骤一，搭建数字图像相关测量系统，获取试件变形前的参考图像及变形后的带有裂纹的变形图像；步骤二、使用裂纹尖端检测方法计算得到变形图像中裂纹尖端的初始位置（，）；步骤三、对该裂纹尖端初始位置点进行整像素匹配，计算其在参考图像上的位置（，）；
步骤四、在参考图像上布置计算点，为每个计算点分配计算子区；步骤五、以复合裂纹位移场函数为各子区内的位移表征函数，以（，）为裂纹尖端初值，使用数字图像相关法计算各子区内位移场函数的待定系数及裂纹尖端终值（，），并以此计算得到各计算点的位移值；步骤六、由各计算点的位移值插值得到全场位移值。
7.进一步，所述的使用i、ii型复合裂纹的位移场函数的表达式如下：（1）其中，u为以裂纹尖端为中心沿裂纹方向的位移，v为以裂纹尖端为中心垂直于裂纹方向的位移，，，，为i型裂纹的待定系数，， ，，为ii型裂纹的待定系数，，表示待定的u方向和v方向的刚体平移，r和θ为以裂纹尖端为原点的极坐标，m为待定系数的个数。
8.进一步，所述的位移函数能够有效描述单一i型裂纹、单一ii型裂纹或i、ii型复合裂纹的位移场。由使用者决定位移函数为单一i型裂纹、单一ii型裂纹或i、ii型复合裂纹，并决定各待定系数的个数，即m的数值，位移场函数中，通过决定m值的大小来权衡计算精度及计算速度，m值越大，计算精度越高，同时计算速度越慢。
9.进一步，所述的步骤二中裂纹尖端检测方法包括：基于图像处理的检测方法、基于机器学习的检测方法及基于深度卷积神经网络的检测方法。
10.进一步，所述的步骤三中裂纹尖端整像素匹配具体为：在变形图像上以裂纹尖端初始位置（，）为中点，分配变形子区，寻找参考图像上与该变形子区相关性最高的参考子区，该参考子区中点坐标即为（，）。
11.进一步，所述的步骤五中计算裂纹尖端终值及位移场函数待定系数的方法具体为：以公式（1）作为参考图像各计算点子区内的位移表征函数，以（，）为裂纹尖端初值，使用数字图像相关法中的牛顿迭代法进行迭代求解与参考子区内灰度信息相关性匹配度最高的变形图像子区信息，从而计算得到公式（1）的待定系数，，，，，，，，，及裂纹尖端终值（，）；以裂纹尖端终值为极坐标原点计算参考子区中点的极坐标，然后根据公式（1）求得该子区中点的位移值u，v。
12.本发明与现有技术的相比的有益效果在于：本发明的方法首先通过智能识别算法识别出变形图像中裂纹尖端的位置，并通过
整像素匹配计算裂纹尖端在参考图像中的位置。然后应用i、ii型复合裂纹位移场函数作为数字图像相关方法中计算子区的位移表征函数进行子区匹配从而迭代求得裂纹尖端终值及位移场函数的待定系数。根据裂纹尖端终值及位移场函数的待定系数可求得子区中点的位移值。
13.本发明应用i、ii型复合裂纹位移场函数作为数字图像相关方法中计算子区的位移表征函数，使用者可以根据实际情况决定位移场函数为单一i型裂纹、单一ii型裂纹或i、ii型复合裂纹。并决定各待定系数的个数，即m的数值，通过决定m值的大小来权衡计算精度及计算速度，m值越大，计算精度越高，同时计算速度越慢。该复合裂纹位移场函数可以有效表征计算子区由裂纹导致的子区断裂情况下的位移场，从而保证数字图像相关方法在处理含裂纹图像时的计算精度。
附图说明
14.图1为本发明一种应用复合裂纹位移场函数的数字图像相关方法的实施例中参考图像与变形图像及裂纹位置；图2为本发明一种应用复合裂纹位移场函数的数字图像相关方法的实施例中真实位移场分布图；图3为本发明一种应用复合裂纹位移场函数的数字图像相关方法的实施例中计算位移场分布图；图4为本发明一种应用复合裂纹位移场函数的数字图像相关方法的实施例中计算位移场与真实位移场的误差分布图。
具体实施方式
15.为使本发明的目的、技术方案及效果更加清楚，明确，以下列举实施例对本发明进一步详细说明。应当指出此处所描述的具体实施仅用以解释本发明，并不用于限定本发明。
16.本发明的方法所述的方法为使用i、ii型复合裂纹的位移场函数作为数字图像相关方法的子区位移表征函数，表达式如下：（1）其中，u为以裂纹尖端为中心沿裂纹方向的位移，v为以裂纹尖端为中心垂直于裂纹方向的位移，，，，为i型裂纹的待定系数，， ，，为ii型裂纹的待定系数，，表示待定的u方向和v方向的刚体平移，r和θ为以裂纹尖端为原点的极坐标。
17.所述的方法具体步骤为：步骤一、搭建数字图像相关测量系统，获取试件变形前的参考图像及变形后的带有裂纹的变形图像；
步骤二、使用裂纹尖端检测方法计算得到变形图像中裂纹尖端的初始位置（，）；步骤三、对该裂纹尖端初始位置点进行整像素匹配，计算其在参考图像上的位置（，）；步骤四、在参考图像上布置计算点，为每个计算点分配计算子区；步骤五、以公式（1）作为参考图像各计算点子区内的位移表征函数，以（，）为裂纹尖端初值，使用数字图像相关法中的牛顿迭代法进行迭代求解与参考子区内灰度信息相关性匹配度最高的变形图像子区信息，从而计算得到公式（1）的待定系数，，，，，，，，，及裂纹尖端终值（，）；以裂纹尖端终值为极坐标原点计算参考子区中点的极坐标，然后根据公式（1）求得该子区中点的位移值u，v；步骤六、由各计算点的位移值插值得到全场位移值。
18.以下列举具体的实施例进行叙述本发明的方法：模拟散斑图像由计算机高斯散斑算法生成，像素为501
×
501，散斑数目为1500，生成参考图像；为了简便，设定i型裂纹位移场函数如下： （2）其中裂纹尖端位置为（251，251），生成变形图像，如图1所示。通过上述公式（2）计算位移场，如图2所示，然后将计算结果与真实解进行对比。
19.具体实施步骤如下：步骤一：通过图像处理的方法获得变形图像裂纹尖端初始位置（249，250）；步骤二：通过整像素匹配计算，裂纹尖端在参考图像中的位置为（250，250）；步骤三：在参考图像上以2个像素的间隔布置计算点，每个计算点的计算子区大小为181
×
181；步骤四：假设子区内i型裂纹位移场函数为公式（2）；步骤五：对每个计算点进行位移迭代计算，并与真实位移值进行比较；步骤六：插值得到全场位移场。
20.图3所示为计算得到的全场位移场。
21.图4所示为计算位移场与真实位移场的误差分布。
22.将计算位移场与真实位移场进行比较发现，u向位移场误差均值为-0.0013，方差为0.0112，v向位移场误差均值为-0.0086，方差为0.0367。可见本方法计算所得位移场误差较小，可有效解决传统数字图像相关法对于裂纹区域存在的位移计算不准确甚至计算缺失的问题。
23.最后应当说明的是，以上实施方案仅用以说明本发明的实现方式而非对其限制；人们应该理解，对该发明的实施过程进行修改或者部分算法过程进行同等替换，而不会脱离本发明技术方案的精神，其均应涵盖在本发明请求保护的技术方案范围内。