CN105486244B - 基于超高精度数字图像测量的二维应变测量系统 - Google Patents
基于超高精度数字图像测量的二维应变测量系统 Download PDFInfo
- Publication number
- CN105486244B CN105486244B CN201610061287.8A CN201610061287A CN105486244B CN 105486244 B CN105486244 B CN 105486244B CN 201610061287 A CN201610061287 A CN 201610061287A CN 105486244 B CN105486244 B CN 105486244B
- Authority
- CN
- China
- Prior art keywords
- test specimen
- dimensional strain
- groups
- digital picture
- feature
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明涉及一种基于超高精度数字图像测量的二维应变测量系统,包括中控电脑、标志点模具、摄像头组和显微镜头组;所述标志点模具上具有呈二维布设的三组特征标志孔,其用于贴在被测试件上以确定出被测试件的三组特征标志点;所述摄像头组包括三个摄像头,所述显微镜头组包括三个显微镜头,所述三个显微镜头对应安装在三个摄像头上,以对应测量所述被测试件的三组特征标志点;所述中控电脑与所述三个摄像头通讯连接,以根据所述三个摄像头的测量结果计算所述被测试件的二维应变,并输出二维应变结果。本发明实现了全程自动化,去除了使用移测显微镜测量定标线间距的过程,提高了测量效率,并且从一维应变测量扩展为二维应变测量,提升了准确度。
Description
技术领域
本发明属于应变测量技术领域,具体涉及一种基于超高精度数字图像测量的二维应变测量系统。
背景技术
传统的应变测量方法是在试件表面粘贴应变片,通过应变片的拉伸、压缩所产生的电信号的改变来测量其应变,这种方法是一种直接接触的测量方法,该方法对贴应变片的实验员要求较高,并且在应变片与被测试件之间的粘结剂粘结效果、厚度等对测量结果有多大影响不得而知,特别是在冲击等动态作用下,粘结剂本身的材料响应对被测材料与应变片之间的响应影响是不可以忽视的。再者,很多测试材料是非常不规则的、表面也是凹凸不平的、甚至有的是难以黏贴的,因此传统的应变片在这些情况下难以满足测试要求。而且,试验现场应变片需要粘贴、焊接,还需要配备应变仪、采集仪等,应变片也是一次性的,因此试验成本相对较高。
随着数码硬件的飞速发展,数字图像测量方法逐渐显示出其独特的优势,它是一门新的测量方法,该方法具有测量多样化,测量成本低、操作简单,无需专业的实验员,而且对材料没有任何影响等优点,因此数字图像测量方法在很多领域中都具有广阔的应用前景。
自从1982年,Yamaguchi和Peters各自提出了数字图像相关测量技术以来,数字图像得到了广泛的发展。随后Gong Maohan在对激光散斑相关性和位移的关系的研究基础上指出,白光比激光更适合数字图像相关方面的运算。为了进一步提高计算精度,提出了双线性插值法、样条函数法、多项式法以及形心法等多种方法来得到图像的亚像素。相关算法从理论上来看其精度达到0.01像素甚至更高,但实际的图像由于大量的噪声,其精度难以稳定在0.01像素。数字散斑相关测量系统具有抗噪声干扰能力,该项技术能够消除图像噪声,使得精度稳定在0.01像素。现有的散斑数字图像相关的平面全场应变测量方法已经比较成熟,并且相关的产品在国内外也已经广泛使用。但是,由于其测量的是全场应变,相同的分辨率下所拍摄的区域大,单位区域中的像素太低,难以保证应变测量的精度在1με范围内。杨勇等人在《高精度数字图像相关测量系统及其技术研究》论文中提出了关于一维的应变测量系统,其测量精度在1με范围内。该论文中所测量的两个标志点间距通常为100mm的平均应变,对于尺寸较小的试件难以达到1με精度。再者,该系统中在上下两个标志点处各有两根定标线,该系统在测量前需要使用移测显微镜分别测量出两定标线的间距,由于试验现场难以摆放移测显微镜,因此大大降低了测量效率。
发明内容
基于上述原因,本发明的目的在于,提出一种基于超高精度数字图像测量的二维应变测量系统,去除了使用移测显微镜测量定标线间距的过程,在能够保证测量精度的同时,提高测量效率。
具体地,本发明提出的一种基于超高精度数字图像测量的二维应变测量系统,包括中控电脑、标志点模具、摄像头组和显微镜头组;所述标志点模具上具有呈二维布设的三组特征标志孔,其用于贴在被测试件上以确定出被测试件的三组特征标志点;所述摄像头组包括三个摄像头,所述显微镜头组包括三个显微镜头,所述三个显微镜头对应安装在三个摄像头上,以对应测量所述被测试件的三组特征标志点;所述中控电脑与所述三个摄像头通讯连接,以根据所述三个摄像头的测量结果计算所述被测试件的二维应变,并输出二维应变结果。
在本发明的进一步优选方案中,所述三组特征标志孔中的每一组包括呈“十”字布设的五个特征标志孔。
在本发明的进一步优选方案中,每一组特征标志孔所限定的测量范围为10mm×10mm。
在本发明的进一步优选方案中,各个特征标志孔的直径为0.1mm;布设在同一行或同一列的两组特征标志孔的中心间距为1cm。
在本发明的进一步优选方案中,所述确定出被测试件的三组特征标志点具体是,采用涂料通过标志点模具的三组特征标志孔在被测试件上喷涂出三组特征标志点。
在本发明的进一步优选方案中,计算所述被测试件的二维应变具体是,通过计算所述摄像头组测量到的被测试件样本图像中的灰度值的变化计算所述被测试件的偏移量。
在本发明的进一步优选方案中,计算所述被测试件的二维应变的公式如下:
εα=εxcosα+εysinα
其中,εα为所述被测试件在任意方向的二维应变值,α为应变方向与水平方向的夹角,εx为所述被测试件在水平方向的应变值,εy为所述被测试件在竖直方向的应变值。
在本发明的进一步优选方案中,所述输出二维应变结果包括绘制二维应变曲线。
在本发明的进一步优选方案中,所述摄像头组的图像分辨率为3288×4384。
在本发明的进一步优选方案中,所述中控电脑的图像识别分辨率为0.01像素。
有益效果:本发明提出的超高精度数字图像测量的二维应变测量系统,实现全程自动化,去除了使用移测显微镜测量定标线间距的过程,测量精度能够保证在1με范围内,提高了测量效率,并且从一维应变测量扩展为二维应变测量,相对于一维测量能够获得更全面的数据,提升了精度和准确度。
附图说明
图1是实施例提出的基于超高精度数字图像测量的二维应变测量系统的整体结构原理示意图。
图2是实施例中标志点模具的结构原理示意图。
图3是实施例中利用图2的标志点模具喷涂出的一组特征标志点示意图。
图4是实施例中利用图2的标志点模具喷涂出三组特征标志点过程示意图。
图中:10-中控电脑,21、22、23-摄像头,31、32、33-显微镜头,40-被测试件,50-标志点模具。
具体实施方式
为了便于本领域技术人员理解,下面将结合附图以及实施例对本发明进行进一步描述。
请参阅图1至图4,实施例提出的基于超高精度数字图像测量的二维应变测量系统,包括中控电脑10、标志点模具50、摄像头组和显微镜头组;所述标志点模具50上具有呈二维布设的三组特征标志孔(如图2所示),其用于贴在被测试件40上以确定出被测试件40的三组特征标志点41、42、43;所述摄像头组包括三个摄像头21、22、23,所述显微镜头组包括三个显微镜头31、32、33,所述三个显微镜头31、32、33对应安装在三个摄像头21、22、23上,以对应测量所述被测试件的三组特征标志点41、42、43;所述中控电脑10与所述三个摄像头21、22、23通讯连接,以根据所述三个摄像头21、22、23的测量结果计算所述被测试件40的二维应变,并输出二维应变结果,例如输出绘制二维应变曲线。
本实施例的最佳实施方案中,摄像头组(三个摄像头21、22、23)的图像分辨率为3288×4384,中控电脑的图像识别分辨率为0.01像素;所述三组特征标志孔中的每一组包括呈“十”字布设的五个特征标志孔,如图2所示;各个特征标志孔的直径为0.1mm;布设在同一行或同一列的两组特征标志孔的中心间距为L=1cm;每一组特征标志孔所限定的测量范围(也即摄像头21、22、23的拍摄区域)为Lx×Ly=1.0mm×1.0mm,即水平范围Lx=1.0mm,竖直范围Ly=10mm。
测量原理
本实施例中,基于超高精度数字图像测量来确定测试件40的二维应变主要是,通过计算摄像头组测量到的被测试件样本图像中的灰度值的变化计算所述被测试件40的偏移量。
在测量之前,可将标志点模具50贴在被测试件40上,然后采用涂料通过标志点模具50的三组特征标志孔在被测试件40上喷涂出三组特征标志点,如图4所示。图3为一组特征标志点,为了便于区别,图3中各个特征标志点分别用1至5标出。
五个特征标志点的测量
图3中,以特征标志点5为中心,以特征标志点1、4所在的方向为x轴,以标志点1、2所在的方向为y轴建立直角坐标系。根据数字图像相关方法测量出五个特征标志点在图像中的像素位置(Pj,i,qj,i),其中下标j表示第j组特征标志点,j=1,2,3,下标i为第j组特征标志点中的第i个特征标志点,i=1,2,3,4,5。四周的特征标志点水平方向(x方向)和竖直方向(y方向)间距分别为Lx、Ly,其中Lx=1.0mm,Ly=1.0mm,则:
摄像头21所拍摄图像在水平方向的系统放大率为:
摄像头21所拍摄图像在竖直方向的系统放大率为:
摄像头22所拍摄图像在水平方向的系统放大率为:
摄像头22所拍摄图像在竖直方向的系统放大率为:
摄像头23所拍摄图像在水平方向的系统放大率为:
摄像头23所拍摄图像在竖直方向的系统放大率为:
当图像发生移动后,则变化后的图像像素坐标为(P′j,i,q′j,i),其中i=1,2,3,4,5,j=1,2,3。则可以计算出各组特征标志点的水平方向和竖直方向位移:
应变测量结果
根据三个摄像头21、22、23分别测量的二维位移建立方程,可以计算出三个在各组特征标志点范围内各组特征标志点范围内的二维应变。
所述被测试件在水平方向的应变值εx:
所述被测试件在竖直方向的应变值εy:
所述被测试件在任意方向的二维应变值εα(α为应变方向与水平方向即x方向的夹角):
εα=εxcosα+εysinα
精度分析
本实施例采用数字图像相关算法,其核心计算是样本图像中的灰度值,只要灰度图发生偏移,中控电脑10就可以计算出其偏移量。中控电脑的图像识别分辨率为0.01像素,摄像头组的图像分辨率为3288×4384,拍摄的区域为1mm×1mm。换算成位移后,则各个摄像头测量的水平方向位移分辨率为0.003μm,竖直方向位移分辨率为0.0023μm。水平方向应变分辨率为0.3με,竖向应变分辨率为0.23με。该精度小于1με,能够很好的满足应变测量需要。
图像最大移动的范围是±0.5mm,两个特征标志点的最大相对移动范围为±1.0mm,计算得到其应变测量范围为0~10000με。
具体测量步骤
1、将标志点模具50放置在被测试件40的表面,在其表面三组特征标志孔上喷涂产生三组独立的特征标志点,然后移除标志点模具50。
2、将分别安装有显微镜头31、32、33的三个摄像头21、22、23分别对着三组特征标志点,并保证每组特征标志点的五个特征标志点均在画面中,调整焦距。
3、通过中控电脑10对三个摄像头21、22、23采集的数据进行初始化,得到每个摄像头的系统标定系数。
4、对被测试件40进行加载试验,由中控电脑10实时采集图像并处理,绘制二维应变曲线。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (8)
1.一种基于超高精度数字图像测量的二维应变测量系统,其特征在于,包括中控电脑、标志点模具、摄像头组和显微镜头组;所述标志点模具上具有呈二维布设的三组特征标志孔,其用于贴在被测试件上以确定出被测试件的三组特征标志点;所述三组特征标志孔中的每一组包括呈“十”字布设的五个特征标志孔;所述摄像头组包括三个摄像头,所述显微镜头组包括三个显微镜头,所述三个显微镜头对应安装在三个摄像头上,以对应测量所述被测试件的三组特征标志点;所述中控电脑与所述三个摄像头通讯连接,以根据所述三个摄像头的测量结果计算所述被测试件的二维应变,并输出二维应变结果,计算所述被测试件的二维应变的公式如下:
εα=εxcosα+εysinα
其中,εα为所述被测试件在任意方向的二维应变量,α为应变方向与水平方向的夹角,εx为所述被测试件在水平方向的应变值,εy为所述被测试件在竖直方向的应变值。
2.根据权利要求1所述的基于超高精度数字图像测量的二维应变测量系统,其特征在于,每一组特征标志孔所限定的测量范围为1.0mm×1.0mm。
3.根据权利要求2所述的基于超高精度数字图像测量的二维应变测量系统,其特征在于,各个特征标志孔的直径为0.1mm;布设在同一行或同一列的两组特征标志孔的中心间距为1cm。
4.根据权利要求3所述的基于超高精度数字图像测量的二维应变测量系统,其特征在于,所述确定出被测试件的三组特征标志点具体是,采用涂料通过标志点模具的三组特征标志孔在被测试件上喷涂出三组特征标志点。
5.根据权利要求4所述的基于超高精度数字图像测量的二维应变测量系统,其特征在于,计算所述被测试件的二维应变具体是,通过计算所述摄像头组测量到的被测试件样本图像中的灰度值的变化计算所述被测试件的偏移量。
6.根据权利要求1所述的基于超高精度数字图像测量的二维应变测量系统,其特征在于,所述输出二维应变结果包括绘制二维应变曲线。
7.根据权利要求6所述的基于超高精度数字图像测量的二维应变测量系统,其特征在于,所述摄像头组的图像分辨率为3288×4384。
8.根据权利要求1所述的基于超高精度数字图像测量的二维应变测量系统,其特征在于,所述中控电脑的图像识别分辨率为0.01像素。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610061287.8A CN105486244B (zh) | 2016-01-28 | 2016-01-28 | 基于超高精度数字图像测量的二维应变测量系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610061287.8A CN105486244B (zh) | 2016-01-28 | 2016-01-28 | 基于超高精度数字图像测量的二维应变测量系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105486244A CN105486244A (zh) | 2016-04-13 |
CN105486244B true CN105486244B (zh) | 2018-11-09 |
Family
ID=55673383
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610061287.8A Active CN105486244B (zh) | 2016-01-28 | 2016-01-28 | 基于超高精度数字图像测量的二维应变测量系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105486244B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106558049A (zh) * | 2016-12-08 | 2017-04-05 | 合肥金星机电科技发展有限公司 | 材料变形的高分辨率快速数字图像适配方法 |
CN110763146A (zh) * | 2019-10-31 | 2020-02-07 | 河海大学 | 基于双相机的高精度光学引伸计及测量方法 |
CN112697063B (zh) * | 2021-03-23 | 2021-11-09 | 四川省人民医院 | 一种基于显微视觉的芯片应变测量方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103542815A (zh) * | 2013-09-23 | 2014-01-29 | 西安新拓三维光测科技有限公司 | 大幅面散斑全场应变测量方法 |
CN103994724A (zh) * | 2014-05-13 | 2014-08-20 | 滕军 | 基于数字图像处理技术的结构二维位移及应变监测方法 |
CN104296717A (zh) * | 2014-09-29 | 2015-01-21 | 中国石油天然气股份有限公司 | 一种管道弯曲应变的获取方法及测量设备 |
CN205038089U (zh) * | 2015-08-21 | 2016-02-17 | 四川烟草工业有限责任公司 | 包装膜样品处理装置 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4836067B2 (ja) * | 2005-05-23 | 2011-12-14 | 日立造船株式会社 | 構造物における変形量計測方法 |
JP2007263611A (ja) * | 2006-03-27 | 2007-10-11 | Yokohama Rubber Co Ltd:The | 歪測定装置および歪測定方法 |
CN101033962B (zh) * | 2007-02-12 | 2011-07-20 | 三峡大学 | 一种基于光学的模型试验位移测量方法 |
CN101514890B (zh) * | 2009-03-27 | 2010-11-10 | 清华大学 | 一种基于光学剪切的二维光学应变花测量方法 |
CN101825436B (zh) * | 2010-05-14 | 2011-10-26 | 上海伦恩通用技术有限公司 | 纸张三维变形量检测系统 |
CN103017651A (zh) * | 2011-09-23 | 2013-04-03 | 苏州比特速浪电子科技有限公司 | 可翻转样品六面检测装置 |
CN205785092U (zh) * | 2016-01-28 | 2016-12-07 | 广州大学 | 基于超高精度数字图像测量的二维应变测量系统 |
-
2016
- 2016-01-28 CN CN201610061287.8A patent/CN105486244B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103542815A (zh) * | 2013-09-23 | 2014-01-29 | 西安新拓三维光测科技有限公司 | 大幅面散斑全场应变测量方法 |
CN103994724A (zh) * | 2014-05-13 | 2014-08-20 | 滕军 | 基于数字图像处理技术的结构二维位移及应变监测方法 |
CN104296717A (zh) * | 2014-09-29 | 2015-01-21 | 中国石油天然气股份有限公司 | 一种管道弯曲应变的获取方法及测量设备 |
CN205038089U (zh) * | 2015-08-21 | 2016-02-17 | 四川烟草工业有限责任公司 | 包装膜样品处理装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105486244A (zh) | 2016-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103575227B (zh) | 一种基于数字散斑的视觉引伸计实现方法 | |
CN103714535B (zh) | 双目视觉测量系统中摄像机参数在线调整方法 | |
CN104729665B (zh) | 一种基于双目视觉的柔性悬臂板振动检测装置及方法 | |
WO2015188515A1 (zh) | 基于纵横转换放大的宏微复合光栅尺测量系统 | |
CN105486244B (zh) | 基于超高精度数字图像测量的二维应变测量系统 | |
CN111174961B (zh) | 一种基于模态分析的索力光学测量方法及其测量系统 | |
CN104048744A (zh) | 一种非接触式的基于影像的实时在线振动测量方法 | |
CN108489413A (zh) | 一种混凝土表面非均匀应变测量系统与方法 | |
CN107817044B (zh) | 基于机器视觉的板材振动的测量装置及方法 | |
CN103808567A (zh) | 一种钎焊接头力学性能实验装置及方法 | |
CN105699982B (zh) | 双激光标定的高精度摄像头芯片多点测距装置及方法 | |
CN1776358A (zh) | 一种非接触式测量物体在三维空间位移和形变的新方法 | |
WO2013060218A1 (zh) | 臂架振动位移测量方法、系统及包含其的混凝土泵车 | |
CN106197292A (zh) | 一种建筑物位移监测方法 | |
CN107907064A (zh) | 一种裂隙监测系统及方法 | |
CN204679637U (zh) | 双激光标定的高精度ccd多点测距装置 | |
CN205785092U (zh) | 基于超高精度数字图像测量的二维应变测量系统 | |
CN103278104B (zh) | 一种用于dic测量的双相机系统的标定板及其标定方法 | |
CN109242909B (zh) | 一种面向高精度二维尺寸测量的线阵相机标定算法 | |
CN102012215A (zh) | 基于数字图像的非接触式光学应变测量方法及应变计 | |
CN113884049A (zh) | 一种使用全站仪测量起重机水平弯曲的方法 | |
CN1268892C (zh) | 基于位置传感器件psd的三维测量方法 | |
CN201852566U (zh) | 基于数字图像的非接触式光学应变计 | |
CN207215021U (zh) | 一种基于双目立体视觉的摄像机光轴偏移校正装置 | |
CN108413870A (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 |