CN105157588B - 一种应变局部化带间距演变规律的多维同步优化测量方法 - Google Patents
一种应变局部化带间距演变规律的多维同步优化测量方法 Download PDFInfo
- Publication number
- CN105157588B CN105157588B CN201510557390.7A CN201510557390A CN105157588B CN 105157588 B CN105157588 B CN 105157588B CN 201510557390 A CN201510557390 A CN 201510557390A CN 105157588 B CN105157588 B CN 105157588B
- Authority
- CN
- China
- Prior art keywords
- strain
- localization band
- strain localization
- band
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提供了一种应变局部化带间距随应力或应变演变规律的观测方法,包括:利用拍摄设备获取试样在加载过程中的数字图像;利用亚像素数字图像相关方法和双三次样条插值方法获取光滑性较好的应变场;获取每条被测应变局部化带的应变数据;设置迭代的初值,任一初值的维数比应变局部化带的数目多1,利用群智能算法,通过不断同步迭代,逐渐缩小各条应变局部化带上的数据与试探解答之间的偏差,最终获得唯一的应变局部化带角度和各自的不同截距;利用两条平行线之间的距离公式,计算任意两条相邻应变局部化带的间距。本发明的优点为:能够保证测量结果的准确性和唯一性,能够获得应变局部化带间距的演变规律,能够实现自动化、批量、高效率测量。
Description
技术领域
本发明涉及一种应变局部化带间距测量技术,属于岩土力学、岩土工程、实验力学、地质力学、工程材料力学性能测试等领域。
背景技术
应变局部化是在材料破坏之前观察到的应变集中于狭窄的带状区域内的现象,应变局部化按载荷类型可以划分为剪切应变局部化、拉伸应变局部化和压缩应变局部化。应变局部化现象可以在不同的层次上被观察到,例如,晶粒尺度上的网状滑移线、工程尺度上的交叉的小断层和地壳尺度上的网状地质构造和地震带都可归结于应变局部化现象(王学滨,潘一山.地质灾害中的应变局部化现象.地质灾害与环境保护,2001,12(4):1-5;王学滨,赵扬锋,代树红等.地震块体模型的共轭剪切破裂带数值模拟.防灾减灾工程学报,2004,24(2):119-125)。应变局部化现象往往以网络形式存在,这样,应变局部化带间距的测量就成为应变局部化现象研究中不可回避的重要问题。
目前,应变局部化带间距的测量主要采用两种方法:采用尺子和数字图像技术测量肉眼可见的裂纹面之间的间距,以此作为应变局部化带的间距。实际上,这种测量方法的误差很大,测量结果因人而异,工作量大,难于实现自动化、批量、高效率测量,所获取的数据量有限,而且,测量的是裂纹面的间距,并非应变局部化带的间距。应变局部化带在启动及发展过程中,由于变形阶段的不同和应力场的调整,应变局部化的方向及位置将发生适当的变化。没有任何理论证明,应变局部化带间距等同于裂纹面间距。欲准确测量应变局部化带间距,必须着眼于应变场,在裂纹出现之前,应变场已变得不均匀,这种变化仅凭肉眼一般不容易观测到,不断发展的数字图像相关方法(王学滨,杜亚志,潘一山.基于DIC粗-细搜索方法的单轴压缩砂样的应变分布及应变梯度的实验研究.岩土工程学报,2012,34(11):2050-2057;王学滨,杜亚志,潘一山.考虑一阶和二阶位移梯度的数字图像相关方法在剪切带测量中的比较.工程力学,2013,30(7):282-287;王学滨,杜亚志,潘一山.单轴压缩湿砂样局部及整体体积应变的数字图像相关方法观测.岩土工程学报,2014,36(9):1648-1656)为应变场的精确测量提供了便利条件。数字图像相关方法是一种光测力学方法,具有测量设备简单、对测量环境要求低、测量精度高的优点。
本发明提出采用数字图像相关方法测量应变局部化带的间距。以亚像素数字图像相关方法和双三次样条插值方法获得的光滑性较好的应变场作为依据,将唯一应变局部化带间距的测量问题归结于一个多维同步优化问题,利用群智能算法,通过不断同步迭代,逐渐缩小各条应变局部化带数据与试探解答(各线段)之间的偏差,最终实现兼顾所有数据的唯一的应变局部化带角度和应变局部化带间距的测量。
发明内容
为了解决现有的应变局部化带间距测量方法的精度低、效率低、获取数据量有限及所得结果不唯一的问题,本发明提供了一种基于数字图像相关方法的应变局部化带间距演变规律的多维同步优化测量方法,提高了测量的效率、精度,可以获取丰富的数据量,所得结果唯一。
本发明的特征在于,包括如下步骤:
步骤1:利用拍摄设备(数码相机或CCD摄像头)获取在加载条件下试样一个包含散斑表面的数字图像;
步骤2:利用亚像素数字图像相关方法,获取试样表面的变形场:水平线应变εx、垂直线应变εy、剪切应变γxy及最大剪切应变γmax;
步骤3:设置计算区域,获取每条应变局部化带中央及其附近的应变数据,该计算区域可以是包含多条应变局部化带的1个四边形区域,也可是包含1条应变局部化带的四边形区域的集合;
对于计算区域是包含多条应变局部化带的1个四边形区域的情形:
首先,设置包含多条应变局部化带的1个四边形计算区域;
之后,根据应变场,确定走向一致的应变局部化带的数目m;
之后,利用插值方法,获取计算区域内光滑性较好的各种应变场;
最后,利用粗估的被测应变局部化带角度和宽度,设置多个狭长四边形区域,每个狭长四边形区域仅包含1条应变局部化带,获取每条应变局部化带中央及其附近的应变数据;
对于计算区域是包含1条应变局部化带的四边形区域的集合的情形:
首先,根据应变场,确定走向一致的应变局部化带的数目m;
之后,利用粗估的被测应变局部化带角度和宽度,设置多个狭长四边形区域作为计算区域,每个狭长四边形区域仅包含1条应变局部化带;
之后,利用插值方法,获取计算区域内光滑性较好的各种应变场;
最后,获取每条应变局部化带中央及其附近的应变数据;
步骤4:对于每条应变局部化带,设置n个迭代初值和有关的迭代参数,对于任一初值(一系列平行线段),初值的维数为m+1,其中包括1个斜率和m个截距,利用群智能算法,通过不断同步迭代,获取各条线段共同的斜率和各自不同的截距;
步骤5:利用平行线之间距离公式,计算任意两条相邻应变局部化带的间距。
进一步地,其中,所述获取在加载条件下试样一个包含散斑表面的数字图像进一步为:若试样表面的天然纹理可作为散斑场,则不必人工制作散斑场,否则,利用油漆、颜料、墨水等喷涂人工散斑场,利用试验机或加载装置对试样进行加载,同时,利用拍摄设备记录散斑场的图像。
进一步地,其中,所述利用亚像素数字图像相关方法,获取试样表面的变形场进一步为:选择若干记录的图像,设置有关的计算参数和计算模式,使用亚像素数字图像相关方法计算图像的位移场和应变场,应变场通过对位移场进行中心差分获得,最大剪切应变γmax由εx、εy及γxy获得
进一步地,其中,所述走向一致的应变局部化带是指应变局部化带网络中其角度基本相同的若干应变局部化带的集合。
进一步地,其中,所述获取每条应变局部化带中央及其附近的应变数据进一步为:设置临界应变参数,将超过该参数的数据作为被测应变局部化带的数据。
进一步地,其中,所述狭长四边形区域的宽度方向尺寸可取为材料平均颗粒直径的20~60倍,即应变局部化带宽度的1~3倍;所述临界应变参数由经验确定,一般可取为计算区域内γmax最大值的60%以上。
进一步地,其中,所述迭代初值包括各条应变局部化带所对应线段的共同斜率的初值和不同截距的初值,故初值的维数比应变局部化带的数目m多1,设任一个迭代初值为Ωi={ki,b1i,b2i,...,bmi},其中,ki为共同斜率的初值,{b1i,b2i,...,bmi}为不同截距的初值。
进一步地,其中,所述群智能算法可在粒子群优化算法、遗传算法、差分进化算法等算法中进行选择。
进一步地,其中,所述各条线段共同的斜率和各自不同的截距是指由这些参数确定的m条线段与各条应变局部化带中央及其附近的应变数据最为接近,即偏差达到最小
其中,下标l代表n个迭代初值中的任一个,l=1~n;ci代表第i条应变局部化带的数据数目,通常,不同应变局部化带的数据数目不同,i=1~m;sij代表第i条应变局部化带的第j个数据的纵坐标,i=1~m,j=1~ci,sij=sij(xij),xij为第i条应变局部化带的第j个数据的横坐标;fij代表与第i条应变局部化带逼近的线段在第j个数据处的纵坐标,由于fij是线性函数,所以,fij(kl,bil)=klxij+bil,kl是任意1个初值的共同斜率,bil是1个初值的不同截距,(sij-fij)2代表1个应变数据点与1条线段之间的偏差的平方,对一条应变局部化带上的所有数据ci求和,代表了一条应变局部化带的应变数据与1条线段之间的偏差的平方,再对每条应变局部化带求和,代表各条应变局部化带的应变数据与各条线段之间的偏差的平方,通过不断同步迭代,使其达到最小的向量Ω*={k*,b1*,b2*,…,bm*}即为所述的唯一且最佳结果。
进一步地,其中,所述同步迭代是迭代初值或试探解答不断更新的过程,对于不同的群智能算法,试探解答的更新原理并不相同,但是基本上都考虑了不同试探解答之间的信息共享、相互学习和影响以及优胜劣汰和适者生存等策略,迭代结束的条件包括两种:1)迭代次数达到最大代数;2)迭代结果已趋于稳定(尽管迭代次数未达到最大代数)。
进一步地,其中,所述迭代参数包括:最大代数、迭代初值的数目、迭代结束的阈值条件以及不同群智能算法所要求的迭代参数,例如,对于粒子群优化算法,需要粒子飞行的最大速度、学习因子、惯性常数等参数。
进一步地,其中,所述任意两条相邻应变局部化带的间距di公式为
di=|b(i+1)*-bi*|sinθ*=|b(i+1)*-bi*|sin[arctan(k*)],i=1~(m-1) (3)
本发明所述的一种基于数字图像相关方法的应变局部化带间距演变规律的多维同步优化测量方法,能够实现应变局部化带间距的自动化、快速、准确测量,获得的信息丰富,能保证测量结果的唯一性。
附图说明
图1为一种应变局部化带间距演变规律的多维同步优化测量方法的流程图;
图2为在加载条件下数字图像采集装置的示意图,图中,1—试样,2—包含散斑的表面,3—光源,4—拍摄设备,5—计算机;
图3为自主开发的基于粒子群优化和Newton-Raphson迭代方法的数字图像相关方法程序界面图;
图4为多维同步优化过程的示意图,其中图4-a为设置计算区域;图4-b为确定走向一致的应变局部化带的范围的示意图;图4-c为获取每条应变局部化带中央及其附近的应变数据的示意图;图4-d为设置迭代初值和迭代参数的示意图;图4-e为迭代的最终结果(仅给出了使目标函数最小的真实解答)的示意图;图中6—计算区域,7—应变局部化带,7-1—走向一致的应变局部化带(被测应变局部化带),7-2—走向一致的应变局部化带,8—狭长四边形区域,9—任一数据点,10—某一迭代初值(这样的迭代初值有n个),11—唯一的真实解答;
图5为基于粒子群优化算法的多维同步优化过程的流程图;
图6为计算任意两条相邻应变局部化带的间距的原理图;
图7为不同纵向应变εa时的单向压缩砂样的最大剪切应变γmax的分布图,其中7-a为εa=0.08时的γmax的分布图;图7-b为εa=0.09时的γmax的分布图;图7-c为εa=0.11时的γmax的分布图;图7-d为εa=0.13时的γmax的分布图;图7-e为εa=0.14时的γmax的分布图;图7-f为εa=0.15时的γmax的分布图;图7-g为εa=0.16时的γmax的分布图;12—限制3条应变局部化带的3个四边形区域;
图8为纵向应变εa=0.16时走向一致的3条应变局部化带的数据和迭代初值分布图,其中,图8-a为迭代代数N=1时的分布图;图8-b为迭代代数N=10时的分布图;图8-c为迭代代数N=30时的分布图;图8-d为迭代代数N=100时的分布图;图中,13-15—迭代初值,16—应变局部化带的数据,17—真实解答,17-1—针对第1条带的真实解答,17-2—针对第2条带的真实解答,17-3—针对第3条带的真实解答;
图9为εa=0.16且不同迭代代数时目标函数J的演化图;
图10为不同纵向应变εa时被测应变局部化带的间距的演化图。
具体实施方式
以下结合附图对该方法的具体实施方式进行说明。
本发明所述的一种应变局部化带间距演变规律的多维同步优化测量方法,其流程图如图1所示,包括如下步骤:
步骤1:如图2所示,利用拍摄设备4(数码相机或CCD摄像头)获取在加载条件下试样1一个包含散斑的表面2的数字图像,光源3设置于拍摄设备4的前方,拍摄的图像保存在计算机5的硬盘中;
这一步骤实施具体如下:首先,在试样1一个表面制作人工散斑场,若试样表面的天然纹理可达到作为散斑场的条件,则不必人工制作散斑场;然后,利用试验机或加载装置对试样1进行加载,与此同时,利用拍摄设备4记录散斑场的图像。
所述人工散斑场是指在试样表面采用油漆、颜料、墨水等喷涂的随机斑点状图案。
步骤2:利用亚像素数字图像相关方法(图3),获取试样表面的变形场:水平线应变εx、垂直线应变εy、剪切应变γxy及最大剪切应变γmax;
这一步骤具体实施如下:首先,从记录的大量图像中,选择若干典型的用于应变场计算;之后,设置有关的计算参数和计算模式,使用亚像素数字图像相关方法计算图像的位移场;之后,通过对位移场进行中心差分获得应变场;最后,根据水平线应变、垂直线应变和面内剪切应变计算最大剪切应变。
所述计算参数包括:子区尺寸、测点间距、测点数目、粒子数目及粒子最大飞行速度等。
所述计算模式包括:增量模式和全量模式,所谓增量模式是指以前一张图像作为变形前图像,所得的位移场和应变场是增量;所谓的全量模式是指以第一张图像作为变形前图像,所得的位移场和应变场是全量。
所述亚像素数字图像相关方法可选择基于粒子群优化及Newton-Raphson迭代方法的数字图像相关方法(杜亚志,王学滨.基于Newton-Raphson迭代与PSO算法的数字图像相关方法.计算机工程与应用,2012,48(34):184-189),该方法可以避免传统Newton-Raphson迭代方法易于陷入局部最优的弊端,也可以避免迭代的初值不容易确定的弊端。尽管该方法可以同时对位移场和应变场进行相关搜索,但考虑到应变场的精度一般不高,因而弃用应变场的数据,通常采用对位移场进行中心差分获得新的应变场。
所述最大剪切应变场由3种应变场获得,该最大剪切应变场总大于零,随着变形的继续,应变局部化带内部的变形愈加明显,而带外的变形增加较慢,因而可以很好地表征应变局部化过程。
步骤3:设置计算区域,获取每条应变局部化带中央及其附近的应变数据;
这一步骤实施具体如下:首先,通过观察应变场,设置一个包含多条应变局部化带7的计算区域(这里,以包含多条应变局部化带的1个四边形区域作为计算区域6,如图4-a所示);之后,根据应变场,确定走向一致的应变局部化带7的数目m(图4-b);之后,对上述计算区域的数据量有限的应变场进行双三次样条插值,以提高应变场的数据量和光滑性,不需对整个试样进行,这可降低插值工作量,提高插值的效率;最后,如图4-b和图4-c所示,利用粗估的被测应变局部化带7-1的角度和宽度,剔除被测应变局部化带7-1之外的数据,获取每条被测应变局部化带7-1中央及其附近的应变数据。
所述走向一致的应变局部化带是指应变局部化带网络中其角度基本相同的若干局部化带的集合,通常,应变局部化带网络中的应变局部化带7是以共轭的方式存在的,即两簇,其中的任一簇即为走向一致的应变局部化带7-1或7-2,此处,以走向一致的应变局部化带7-1作为被测应变局部化带。
所述获取每条被测应变局部化带中央及其附近的应变数据进一步为:首先,设置适当的临界应变参数,将超过该参数的数据作为被测应变局部化带7-1的数据;之后,根据经验,确定各条被测应变局部化带7-1的范围,舍弃在这些范围之外的数据。
所述各条被测应变局部化带的范围由1个狭长四边形区域8限定,狭长四边形区域8均位于计算区域6之内,狭长四边形区域8的两个较长的边平行,其角度凭借经验设置,对于岩土材料,剪切应变局部化带的方向与最大主应力(σ3)方向之间的角度一般在55°~80°;狭长四边形区域8的两个较长的边的距离,即其宽度方向的尺寸可取为剪切应变局部化带宽度的1~3倍,即材料平均颗粒直径的20~60倍。
所述临界应变参数可取为计算区域6内γmax的最大值的60%以上。
步骤4:如图4-d所示,对于每条被测应变局部化带7-1,设置n个迭代初值和有关迭代参数,对于某一迭代初值10,即一系列平行线段,初值的维数为m+1,其中包括1个斜率和m个截距,利用群智能算法,通过不断同步迭代,获取各条线段共同的斜率和各自不同的截距。
所述迭代初值是指计算区域6之内的一系列线段,每条线段针对1条应变局部化带7,迭代初值由这些线段的斜率和截距完全确定,坐标原点可取在计算区域6的中心处;对于某一迭代初值10,其包含的线段彼此平行,这样,初值的维数将比应变局部化带的数目m多1;对于不同的初值,通常,其包含的全部线段并不满足相互平行的要求。
所述任一迭代初值设为Ωi={ki,b 1i,b2i,...,bmi},其中,i=1~n,n为迭代初值的数目;应当指出,该迭代初值的数目不等于初值的维数,维数和应变局部化带7的数目有关,而初值的数目与问题的难易程度有关,对于容易求解的优化问题,初值的数目不需太多;ki为任一迭代初值的共同斜率,{b1i,b2i,...,bmi}为任一迭代初值的不同截距;在迭代过程中,迭代初值将被试探解答更新,但试探解答仍满足迭代初值的特点:对于任一个试探解答,其包含的线段彼此平行。
所述群智能算法包括粒子群优化算法、遗传算法、差分进化算法等。不同的群智能算法有不同的表现和起源,这决定了试探解答不同的更新原理,但这些群智能算法具有下列共性:具有较强的全局并行搜索能力,这根源于考虑了不同试探解答之间的信息共享、相互学习和影响以及优胜劣汰和适者生存等策略。
所述同步迭代过程是指试探解答不断向真实解答逼近的过程,在此过程中,m条线段与m条应变局部化带中央及附近的应变数据之间的偏差同步逐渐缩小,上述迭代或寻优过程可归结于使目标函数J达到最小,即
其中,下标l代表n个迭代初值中的任一个,l=1~n;ci代表第i条应变局部化带的数据数目,通常,不同应变局部化带的数据数目不同,i=1~m;sij代表第i条应变局部化带的第j个数据的纵坐标,j=1~cj,sij=sij(xij),xij为第i条应变局部化带的第j个数据的横坐标;fij代表与第i条应变局部化带逼近的线段在第j个数据处的纵坐标,由于fij是线性函数,所以,fij(kl,bil)=klxij+bil,kl是任意1个初值的共同斜率,bil是1个初值的不同截距,sij-fij代表1个数据点9与1个试探解答中1条线段之间的偏差,欲获取所有数据点9与一个试探解答中所有线段的偏差的平方,需要进行两次求和计算:第1次,对1条应变局部化带7中的所有数据的偏差的平方求和,第2次,对各条应变局部化带7的偏差的平方求和;通过不断同步迭代,使上述偏差的平方达到最小的试探解答记为Ω*,Ω*={k*,b1*,b2*,…,bm*},其具有唯一性,即为所求的真实解答13。如果对各条应变局部化带7的数据分别拟合,则无法保证每条应变局部化带7的角度相同,这样,对于若干角度稍有差别的若干线段,难于采用1个确定的公式获取应变局部化带间距的唯一结果。通过设置迭代初值的维数比应变局部化带的数目m多1,而且,其中仅包含一个斜率,这是为了保证应变局部化带7间距的迭代结果的唯一性和准确性,其唯一性和准确性得到保障是由于应变局部化带7角度的迭代结果兼顾了所有走向一致的应变局部化带7中央及其附近的所有数据。
所述迭代参数包括:最大代数、迭代初值的数目、迭代结束的阈值条件及不同群智能算法的具体参数。下面,给出采用粒子群优化算法进行迭代的流程图(图5)和计算步骤:
(1)初始化粒子在计算区域内6,根据被测应变局部化带7的形态,给出其中心线角度和截距的估计值,以此产生n个粒子的初值,每个初值的维数为m+1,当前粒子为第1代粒子,N=1;
(2)评价各粒子利用式(1),计算各粒子的优劣,即评价各粒子与应变局部化带数据的偏离程度;
(3)寻找任一粒子最优试探解答pbest对于任一粒子,在所有迭代代数中,寻找最优试探解答pid,即pbest,下标i为粒子编号,i=1~n,d为某一维数,d=1~m+1;
(4)寻找所有粒子中的最优试探解答gbest在所有迭代代数中,寻找所有粒子中的最优试探解答gbest;
(5)更新试探解答利用下式更新粒子的飞行速度和新位置
其中,代表第i个粒子在第N+1代时第d维上的速度;c1、c2为加速常数,一般取c1=c2=2;r1和r2分别为取值在[0,1]之间的随机数;w为惯性常数,取值随着迭代代数N的增加而线性递减,其最大值wmax和最小值wmin分别为1.4和0;Nmax为预设的最大迭代次数;代表第i个粒子在第N+1代时第d维上的坐标。如果粒子更新后的速度在某维上超出了范围[-vmax,vmax],则被限定在该维的边界上。如果粒子更新后的坐标超出了搜索域,也将其限定在边界上。
(6)如果满足迭代结束的条件,即迭代结果趋于稳定,则停止迭代。否则,迭代代数增加1,回到(2)继续计算,直到满足迭代结束的条件。
步骤5:如图6所示,利用两条平行线之间距离公式,计算任意两条相邻应变局部化带7的间距di,其公式容易推得
di=|b(i+1)*-bi*|sinθ*,i=1~(m-1) (5)
其中,θ*为线段的角度,θ*与k*有关,所以
di=|b(i+1)*-bi*|sin[arctan(k*)],i=1~(m-1) (6)
当i=1时,d1代表第1条与第2条应变局部化带的间距,以此类推。亦可采用点到直线距离公式计算di,但要求点在一条直线上。
图7-a至图7-g为不同纵向应变εa时的单向压缩砂样的最大剪切应变γmax的分布图,这一结果是采用自主开发的基于粒子群优化和Newton-Raphson迭代方法的数字图像相关方法获得的。可以发现,随着εa的增加,应变局部化带越来越清晰,网络格局越来越明显,从中选取3条走向一致的应变局部化带作为被测对象,其中,将限制3条应变局部化带的3个四边形区域12作为计算区域。图8-a至图8-d为εa=0.16时走向一致的3条应变局部化带的数据、迭代初值及多维同步优化过程图,其中粒子数取为3,即迭代初值13-15共计3个,限于篇幅,未给出其它时的过程图。可以发现,随着迭代代数N的增加,各试探解答不断向各条应变局部化带的数据16靠拢,最终,各试探解答趋向同一目标,即为所求的真实解答17。图9为εa=0.16且不同迭代代数时目标函数J的演化图,可以发现,随着迭代代数N的增加,目标函数J不断降低,这说明试探解答不断向真实解答逼近,当N=60~100时迭代结果已趋于稳定。图10为不同纵向应变εa时被测应变局部化带的间距的演化图,可以发现,不同εa时,应变局部化带之间的距离有所不同,随着εa的增加,针对第1条带的真实解答17-1与针对第2条带的真实解答17-2的间距d1有减小的趋势,针对第2条带的真实解答17-2与针对第3条带的真实解答17-3的间距d2有增大的趋势。
Claims (8)
1.一种应变局部化带间距演变规律的多维同步优化测量方法,该方法具体步骤如下:
步骤1:利用拍摄设备获取在加载条件下试样一个包含散斑表面的数字图像;
步骤2:利用亚像素数字图像相关方法,获取试样表面的变形场:水平线应变ε x、垂直线应变ε y、剪切应变γ xy最大剪切应变γ max;
步骤3:设置计算区域,获取每条应变局部化带中央及其附近的应变数据,该计算区域可以是包含多条应变局部化带的1个四边形区域,也可是包含1条应变局部化带的四边形区域的集合;
对于计算区域是包含多条应变局部化带的1个四边形区域的情形:
首先,设置包含多条应变局部化带的1个四边形计算区域;
之后,根据应变场,确定走向一致的应变局部化带的数目m;
之后,利用插值方法,获取计算区域内光滑性较好的各种应变场;
最后,利用粗估的被测应变局部化带角度和宽度,设置多个狭长四边形区域,每个狭长四边形区域仅包含1条应变局部化带,获取每条应变局部化带中央及其附近的应变数据;
对于计算区域是包含1条应变局部化带的四边形区域的集合的情形:
首先,根据应变场,确定走向一致的应变局部化带的数目m;
之后,利用粗估的被测应变局部化带角度和宽度,设置多个狭长四边形区域作为计算区域,每个狭长四边形区域仅包含1条应变局部化带;
之后,利用插值方法,获取计算区域内光滑性较好的各种应变场;
最后,获取每条应变局部化带中央及其附近的应变数据;
步骤4:对于每条应变局部化带,设置n个迭代初值和有关迭代参数,对于任一初值,即一系列平行线段,初值的维数为m+1,其中包括1个斜率和m个截距,利用群智能算法,通过不断同步迭代,获取各条线段共同的斜率和各自不同的截距;
步骤5:利用平行线之间距离公式,计算任意两条相邻应变局部化带的间距。
2.根据权利要求1所述的一种应变局部化带间距演变规律的多维同步优化测量方法,其特征在于,将应变局部化带间距测量问题归结于多维同步优化问题,通过不断同步迭代,逐渐缩小应变局部化带数据与试探解答之间的偏差,直到满足迭代结束条件;
所述同步迭代是指试探解答向真实解答逼近的过程,在此过程中,m条线段与m条应变局部化带中央及其附近的应变数据之间的偏差同步逐渐缩小,上述迭代或寻优过程可归结于使目标函数J达到最小,即
(1)
其中,下标l代表n个迭代初值中的任一个,l=1~n;c i 代表第i条应变局部化带的数据数目,不同应变局部化带的数据数目不同,i=1~m;s ij 代表第i条应变局部化带的第j个数据的纵坐标,j=1~c j ,s ij =s ij (x ij ),x ij 为第i条应变局部化带的第j个数据的横坐标;f ij 代表与第i条应变局部化带逼近的线段在第j个数据处的纵坐标,由于f ij 是线性函数,所以, f ij (k l ,b il )=k l x ij +b il ,k l 是任意1个初值的共同斜率,b il 是1个初值的不同截距,s ij -f ij 代表1个数据点与1个试探解答中1条线段之间的偏差;欲获取所有数据点与一个试探解答中所有线段的偏差的平方,需要进行两次求和计算:第1次,对1条应变局部化带中的所有数据的偏差的平方求和,第2次,对各条应变局部化带的偏差的平方求和;使上述偏差的平方达到最小的试探解答记为Ω * ,Ω * ={k * ,b 1*,b 2*,…,b m*},其具有唯一性,即为所求的真实解答。
3.根据权利要求1所述的一种应变局部化带间距演变规律的多维同步优化测量方法,其特征在于,对于任一个试探解答,其维数比应变局部化带数目m多1,其中包含各条线段共同的斜率和各自不同的截距。
4.根据权利要求1所述的一种应变局部化带间距演变规律的多维同步优化测量方法,其特征在于,利用群智能算法的全局并行搜索能力,对试探解答的各条平行线段的不同截距和共同的斜率同步搜索,兼顾了所有应变局部化带的数据,确保了所得结果的正确性和唯一性。
5.根据权利要求1所述的一种应变局部化带间距演变规律的多维同步优化测量方法,其特征在于,利用插值方法对所述计算区域内的由亚像素数字图像相关方法获得的数据量有限的应变场进行插值,以获取光滑的应变场,根据应变场的特点,获取被测应变局部化带的数据;
所述应变场的特点包括:(1)最大剪切应变场、线应变场或剪切应变场的统计信息;(2)各条应变局部化带的范围及角度的估计值,上述估计值由直接观察或凭借经验确定。
6.根据权利要求1所述的一种应变局部化带间距演变规律的多维同步优化测量方法,其特征在于,任意两条相邻应变局部化带的间距d i 的计算公式为
(2)
其中,参数的下标*代表使目标函数达到最小的参数,亦可采用点到直线距离公式计算两条相邻应变局部化带的间距,但要求点在一条直线上。
7.根据权利要求1所述的一种应变局部化带间距演变规律的多维同步优化测量方法,其特征在于,被测应变局部化带是走向一致的,即其角度基本相同。
8.根据权利要求1所述的一种应变局部化带间距演变规律的多维同步优化测量方法,其特征在于,通过计算在不同应力或应变条件下试样中的相邻应变局部化带之间的间距,即可获得应变局部化带间距随应力或应变的演变规律。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510557390.7A CN105157588B (zh) | 2015-09-05 | 2015-09-05 | 一种应变局部化带间距演变规律的多维同步优化测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510557390.7A CN105157588B (zh) | 2015-09-05 | 2015-09-05 | 一种应变局部化带间距演变规律的多维同步优化测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105157588A CN105157588A (zh) | 2015-12-16 |
CN105157588B true CN105157588B (zh) | 2018-01-26 |
Family
ID=54798532
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510557390.7A Expired - Fee Related CN105157588B (zh) | 2015-09-05 | 2015-09-05 | 一种应变局部化带间距演变规律的多维同步优化测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105157588B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105654444B (zh) * | 2016-01-27 | 2018-07-13 | 辽宁工程技术大学 | 一种变形数字图案的复原方法 |
CN106979752B (zh) * | 2017-03-24 | 2019-01-25 | 辽宁工程技术大学 | 一种应变局部化带宽度及间距的统一光学测量方法 |
CN107144230B (zh) * | 2017-05-22 | 2019-05-03 | 辽宁工程技术大学 | 一种基于长方形子区的应变局部化带应变的光学测量方法 |
CN107747916B (zh) * | 2017-10-27 | 2019-07-30 | 辽宁工程技术大学 | 一种应变局部化带剪胀角的光学测量方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5065331A (en) * | 1981-05-18 | 1991-11-12 | Vachon Reginald I | Apparatus and method for determining the stress and strain in pipes, pressure vessels, structural members and other deformable bodies |
CN102494614A (zh) * | 2011-11-14 | 2012-06-13 | 河南科技大学 | 一种高精度数字散斑相关测量方法 |
CN104848838A (zh) * | 2015-04-29 | 2015-08-19 | 辽宁工程技术大学 | 两种构形条件下岩土试样剪切带倾角演变规律的观测方法 |
-
2015
- 2015-09-05 CN CN201510557390.7A patent/CN105157588B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5065331A (en) * | 1981-05-18 | 1991-11-12 | Vachon Reginald I | Apparatus and method for determining the stress and strain in pipes, pressure vessels, structural members and other deformable bodies |
CN102494614A (zh) * | 2011-11-14 | 2012-06-13 | 河南科技大学 | 一种高精度数字散斑相关测量方法 |
CN104848838A (zh) * | 2015-04-29 | 2015-08-19 | 辽宁工程技术大学 | 两种构形条件下岩土试样剪切带倾角演变规律的观测方法 |
Non-Patent Citations (4)
Title |
---|
《单轴压缩湿砂样局部及整体体积应变的数字图像相关方法观测》;王学滨 等;《岩土工程学报》;20140605;第36卷(第9期);第1648-1656页 * |
《基于微粒子群优化算法的数字散斑图像相关方法》;梁智锦 等;《激光技术》;20140925;第38卷(第5期);第603-607页 * |
《基于粒子群算法的整像素数字图像相关方法》;杜亚志 等;《计算机工程与应用》;20110517;第48卷(第6期);第200-204页 * |
《考虑一阶和二阶位移梯度的数字图像相关方法在剪切带测量中的比较》;王学滨 等;《工程力学》;20130725;第30卷(第7期);第282-287页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105157588A (zh) | 2015-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105157588B (zh) | 一种应变局部化带间距演变规律的多维同步优化测量方法 | |
CN105761312B (zh) | 一种微地形表面重建方法 | |
CN105157594B (zh) | 一种基于半子区分割法的数字图像相关方法 | |
CN109598714A (zh) | 一种基于图像三维重建及网格曲面的隧道超欠挖检测方法 | |
CN102903145B (zh) | 植物群体形态结构三维重建方法 | |
CN110409369A (zh) | 边坡开挖数字化施工与质量控制方法 | |
CN109033538A (zh) | 一种基于实测结构面参数的裂隙岩体渗透张量的计算方法 | |
CN103983230B (zh) | 一种地面叶面积指数间接测量验证方法 | |
CN105678757B (zh) | 一种物***移测量方法 | |
CN103838852B (zh) | 一种快速查找多块结构化网格对接关系的方法 | |
CN104848838B (zh) | 两种构形条件下岩土试样剪切带倾角演变规律观测方法 | |
CN109086649A (zh) | 卫星遥感图像水体识别方法 | |
CN106530397B (zh) | 一种基于稀疏剖面地质轮廓线的地质面三维重建方法 | |
CN104766335A (zh) | 一种岩土材料变形数字图像相关分析优化方法 | |
CN106157368B (zh) | 一种大范围区域重力场精确建模和重构方法 | |
CN106767481B (zh) | 一种应变局部化带内应变场的半子区相关光学测量方法 | |
Liu et al. | Recent movement changes of main fault zones in the Sichuan-Yunnan region and their relevance to seismic activity | |
CN106443782B (zh) | 一种断层和裂缝发育密度、均匀性以及组合样式评价方法 | |
CN108256218A (zh) | 一种基于实测腰线数据的地下连通巷道精细建模方法 | |
CN105043390A (zh) | 基于泛克里金法的重力场插值方法 | |
CN113552637B (zh) | 一种航空-地面-井中磁异常数据协同三维反演方法 | |
CN110569606B (zh) | 峡谷高应力区硬岩大型地下洞室轴线布置量化设计方法 | |
CN104121864B (zh) | 一种岩石错动裂缝面的间隙评价方法 | |
CN106526667B (zh) | 叠前时间偏移方法和装置 | |
CN107239559A (zh) | 一种基于矢量网格进行空间动目标位置计算的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180126 Termination date: 20200905 |
|
CF01 | Termination of patent right due to non-payment of annual fee |