CN115950371B - 基于光学显微的一点应变测量方法 - Google Patents

基于光学显微的一点应变测量方法 Download PDF

Info

Publication number
CN115950371B
CN115950371B CN202310167033.4A CN202310167033A CN115950371B CN 115950371 B CN115950371 B CN 115950371B CN 202310167033 A CN202310167033 A CN 202310167033A CN 115950371 B CN115950371 B CN 115950371B
Authority
CN
China
Prior art keywords
line segment
strain
point
characteristic
characteristic line
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
Application number
CN202310167033.4A
Other languages
English (en)
Other versions
CN115950371A (zh
Inventor
张强
代万里
吴培楠
王红英
韩贵雷
蒋斌松
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN202310167033.4A priority Critical patent/CN115950371B/zh
Publication of CN115950371A publication Critical patent/CN115950371A/zh
Application granted granted Critical
Publication of CN115950371B publication Critical patent/CN115950371B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明公开了一种基于光学显微的一点应变测量方法,包括:基于光学显微设备对介质被测点进行高倍率放大摄像或拍照,获得被测点局部范围高像素细观特征图像,在图像中选择若干个特征点,将特征点相连获得特征线段;测量特征线段在初始时刻的线段长度及其与坐标轴正方向夹角,计算夹角的余弦值;测量特征线段在后续某一时刻的线段长度,计算特征线段相比于初始时刻长度的相对变化量;基于一点的应变状态与该点任意方向正应变之间的关系式计算被测点的应变,获得应变测量结果。本发明解决了电阻式应变测量方法应变片覆盖面积大、环境适应性差、测试精度低、不适用用导电液体环境等问题。

Description

基于光学显微的一点应变测量方法
技术领域
本发明属于应变测量领域,特别是涉及一种基于光学显微的一点应变测量方法。
背景技术
一点的应变测量通常由电阻式应变花进行测量,每个应变花包含不少于3个不同方向的应变片,为了保证测量精度电阻应变片的长度普遍大于1cm。因此,每个应变片测得的应变应为应变片覆盖范围沿应变片长度方向的平均正应变,进而由应变计算得到的一点应变近似为应变花覆盖区域的平均应变值,因此,电阻应变测量方法无法精确获得微小区域的真实应变。
由于电阻式应变测量需将应变片与被测介质采用粘结剂进行紧密粘贴,介质在变形过程中一定程度上会受到粘结剂和电阻式应变片的束缚,其并非是完全自由变形状态,尤其是测量对象为质地较软时,这种影响将非常显著,导致错误的测量结果。若介质处于潮湿环境或介质不致密时,电阻式应变片粘贴效果难以保障,进而产生错误的测试数据,甚至致使测试失败。当被测介质浸于导电的液体中时(如充水钻孔、海上轮船底部等),电阻式应变无法进行直接测量,必须临时封隔出干燥腔体进行测量,给测试工作带来的很大的困难。即便如此,由于腔体内介质不受液体压力作用,其应变状态与正常工作条件时并不相同。电阻式应变花为一次性测量元件,大量测量时其成本相对较大。
发明内容
为解决上述问题,本发明提供了如下方案:基于光学显微的一点应变测量方法,包括:
基于光学显微设备对介质被测点进行高倍率放大摄像或拍照,获得被测点局部范围高像素细观特征图像,在图像中选择若干个特征点,将所述特征点相连获得特征线段,基于所述特征线段获得特征线段集;
测量所述特征线段在初始时刻的线段长度及所述特征线段与坐标轴正方向夹角,计算夹角的余弦值;
基于所述特征线段在后续某一时刻的线段长度与所述特征线段在初始时刻的线段长度,计算特征线段相比于初始时刻长度的相对变化量;
基于所述特征线段长度的相对变化量、一点的应变状态与该点任意方向正应变之间的关系,计算被测点的应变,获得应变测量结果。
优选地,所述光学显微设备的放大倍率在1倍~100000倍之间连续可调。
优选地,所述特征图像像素不低于100万。
优选地,所述特征线段的长度及所述特征线段与坐标轴正方向之间的夹角基于所述光学显微设备或图像分析软件的虚拟标尺、标记和像素点位置识别功能实现。
优选地,所述特征点为与被测介质固接且同步变形的可分辨物理特征点,至少包括尖点、凹点、散斑、颜色。
优选地,所述特征线段为从一个特征点中心位置至另一个特征点中心位置的直线段。
优选地,所述若干个特征点为特征点数量满足两两相连所得的特征线段中不平行线段数量不少于3条;
所述特征线段集为根据需要连接任意两个特征点得到的特征线段的集合,且不平行线段数量不少于3条。
优选地,所述相对变化量在数值上等于被测点应变在线段方向上的正应变,其计算表达式为:
式中,为特征线段的初始长度,/>为被测点产生变形后特征线段在T时刻的线段长度;/>为T时刻特征线段的相对变化量(正应变)。
优选地,所述一点的应变状态为:ε=[εxyxy];
所述一点应变与任意方向正应变关系为:
式中,l和m分别为线段与x轴和y轴正方向夹角的余弦值;为一点应变ε在该方向的正应变。
优选地,所述T时刻各特征线段方向上的正应变可用被测点的应变状态分量εx、εy、γxy表示为:
式中,li和mi分别为在初始时刻各线段与x轴和y轴正向夹角的余弦值。
优选地,当特征线段集中不平行特征线段数量等于3时,可构建一点应变计算方程组,求解该方程组可得一点应变分量εx、εy、γxy;当特征线段集中不平行特征线段数量大于3时,所构建方程组为超静定方程组,利用最小二乘法可求解一点应变的分量εx、εy、γxy的最优解;
所述一点应变分量εx、εy、γxy计算表达式可统一表示为:
本发明公开了以下技术效果:
本发明提供的基于光学显微的一点应变测量方法,基于光学显微设备对介质被测点进行高倍率放大摄像或拍照,获得被测点局部范围高像素细观特征图像,在图像中选择若干个特征点,将特征点相连获得特征线段,测量特征线段在初始时刻的线段长度及其与坐标轴正方向夹角,并计算夹角的余弦值,由特征线段在后续某一时刻的线段长度可计算其相比于初始时刻长度的相对变化量,基于一点的应变状态与该点任意方向正应变之间的关系式计算被测点的应变,获得应变测量结果。解决了电阻式应变测量方法应变片覆盖面积大、环境适应性差、测试精度低、不适用用导电液体环境等问题。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的方法流程图;
图2为本发明实施例的初始时刻特征点空间关系示意图;
图3为本发明实施例的T时刻特征点空间关系示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1所示,本发明提供了一种基于光学显微的一点应变测量方法,包括:基于光学显微设备对介质被测点进行高倍率放大摄像或拍照,获得被测点局部范围高像素细观特征图像,在图像中选择若干个特征点,将特征点相连获得特征线段,测量特征线段在初始时刻的线段长度及其与坐标轴正方向夹角,并计算夹角的余弦值,由特征线段在后续某一时刻的线段长度可计算其相比于初始时刻长度的相对变化量,基于一点的应变状态与该点任意方向正应变之间的关系式计算被测点的应变,获得应变测量结果。
进一步地优化方案,所述光学显微设备放大倍率可在10倍~100000倍之间连续可调,特征图像像素不低于100万。
进一步地优化方案,所述光学显微设备或图像分析软件具有虚拟标尺、标记和像素点位置识别功能,可测量两点之间的距离、两线段夹角。
进一步地优化方案,所述特征点指与被测介质固接且同步变形的可分辨物理特征点,至少包括尖点、凹点、散斑、颜色。
进一步地优化方案,所述特征线段指从一个特征点中心位置至另一个特征点中心位置的直线段。
进一步地优化方案,所述若干个特征点指特征点数量满足两两相连所得的特征线段中不平行线段数量不少于3条;
进一步地优化方案,所述一系列特征线段集指根据需要连接任意两个特征点得到的特征线段的集合,且不平行线段数量不少于3条。
进一步地优化方案,所述特征线段与坐标轴正方向夹角指以线段任一端点为起点,指向线段另一端的射线分别与x轴和y轴正方向的夹角,其取值范围为[-π,π]。
更进一步,特征线段与坐标轴正方向夹角的余弦值值域范围为[-1,1]。
进一步地优化方案,所述特征线段长度相比于初始时刻的相对变化量指线段当前长度与其初始长度的差值与线段初始长度的比值。
设一点的应变状态为ε=[εxyxy],则该点在任意方向的正应变为:
式中,l和m分别为该线段与x轴和y轴正方向夹角的余弦;为一点应变ε在该方向的正应变。
采用光学显微方法对被测点进行放大摄像或拍照,获得被测点高像素特征图像并选择待测目标微区域内的k个可识别特征点,记做p1、p2、……、pk,将特征点相连可得n条特征线段,特征线段数量n不超过k(k-1)/2,在初始时刻各特征线段与x轴和y轴正向夹角的余弦值为分别为li和mi、初始长度为(i=1,2,…,n)。当测点区域受到外部扰动产生变形后,特征点位置产生变化,若T时刻测量得到的上述特征线段长度为/>则特征线段相比于初始时刻长度的相对变化量可表示为:
鉴于特征图像为光学显微得到的极小区域的微视场,被测点应变也较小,由于特征点位置改变引起的线段方向的相对变化量相比于应变为高阶小量,可忽略不计,在变形过程中可认为各线段方向不发生变化,各线段相对变化量即为被测点应变在线段方向上的正应变。因此,在T时刻各线段的正应变可用被测点应变状态分量εx、εy、γxy表示为:
当特征线段集中不平行特征线段数量等于3时,可构建3个线性无关方程,进而获得测点的三个应变分量;当特征线段集中不平行特征线段数量大于3时,所构建方程组为超静定方程组,利用最小二乘法可求解该点应变分量的最优解,具体的法方程组为:
则该点的应变分量可直接统一表示为
所述不在同一直线的特征点数量不少于3个,宜取不少于5个特征点以提高计算精度,避免客观测量误差导致错误计算结果。
实施例一
本发明通过以下实施例对本技术方案进行详细描述。
(1)在对被测点进行100倍放大、采用5000万像素CCD相机拍照,在待测点区域附近选取k=4个特征点分别为p1、p2、p3、p4,且任意三个特征点不在同一直线,由四个测点两两相连可得6条特征线段,如图2所示。在具体测试分析时可取所有或部分特征线段,特征线段数量满足不平行线段数量不少于3条即可。为提高测试精度本案例中取所有线段进行测试分析。在初始时刻,各特征线段长度分别为L1=0.5545mm、L2=0.6194mm、L3=0.6953mm、L4=0.6300mm、L5=0.8807mm、L6=0.8812mm,其方向余弦(与x轴和y轴夹角的余弦值)依次为(0.2853,0.9584)、(0.9853,-0.1708)、(-0.2685,-0.9633)、(-0.9221,0.3869)、(0.8742,0.4856)、(0.4782,-0.8782)。
(2)结构受到外力作用后产生变形,导致待测区域选取的特征点产生位移,在T时刻各特征点当前位置为p1′、p2′、p3′、p4′,相应的特征点两两相连特征线段长度变为L1′=0.5541mm、L2′=0.6194mm、L3′=0.6946mm、L4′=0.6329mm、L5=0.8807mm、L6=0.8812mm,如图3所示。
(3)六条线段当前长度相比于初始时刻长度相对变化量分别为:
上述相对变化量即为待测点区域应变沿着线段方向的正应变。
(4)根据一点应变分量εx、εy、γxy与该点在任意方向正应变的理论关系式可得如下方程:
(5)上述方程组为超静定方程组,其对应的求解法方程组为:
求解该方程组可得待测点的应变,各应变分量分别为
以上所述的实施例仅是对本发明的优选方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。

Claims (7)

1.一种基于光学显微的一点应变测量方法,其特征在于,包括:
基于光学显微设备对介质被测点进行高倍率放大摄像或拍照,获得被测点局部范围高像素细观特征图像,在图像中选择若干个特征点,将所述特征点相连获得特征线段,基于所述特征线段获得特征线段集;
测量所述特征线段在初始时刻的线段长度及所述特征线段与坐标轴正方向夹角,计算夹角的余弦值;
基于所述特征线段在后续某一时刻的线段长度与所述特征线段在初始时刻的线段长度,计算特征线段相比于初始时刻长度的相对变化量;
基于所述特征线段长度的相对变化量、一点的应变状态与该点任意方向正应变之间的关系,计算被测点的应变,获得应变测量结果;
所述光学显微设备的放大倍率在1倍~100000倍之间连续可调,特征图像像素不低于100万;
所述特征线段的长度及所述特征线段与坐标轴正方向之间的夹角基于所述光学显微设备或图像分析软件的虚拟标尺、标记和像素点位置识别功能实现;
当特征线段集中不平行特征线段数量等于3时,构建一点应变计算方程组,求解该方程组可得应变状态分量εx、εy、γxy;当特征线段集中不平行特征线段数量大于3时,所构建方程组为超静定方程组,利用最小二乘法求解应变状态分量εx、εy、γxy的最优解;
所述应变状态分量εx、εy、γxy计算表达式可统一表示为:
式中,l和m分别为线段与x轴和y轴正方向夹角的余弦值;ε为一点应变ε在该方向的正应变。
2.根据权利要求1所述的基于光学显微的一点应变测量方法,其特征在于,所述特征点为与被测介质固接且同步变形的可分辨物理特征点,至少包括尖点、凹点、散斑、颜色。
3.根据权利要求1所述的基于光学显微的一点应变测量方法,其特征在于,所述特征线段为从一个特征点中心位置至另一个特征点中心位置的直线段。
4.根据权利要求1所述的基于光学显微的一点应变测量方法,其特征在于,
所述若干个特征点为特征点数量满足两两相连所得的特征线段中不平行线段数量不少于3条;
所述特征线段集为根据需要连接任意两个特征点得到的特征线段的集合,且不平行线段数量不少于3条。
5.根据权利要求1所述的基于光学显微的一点应变测量方法,其特征在于,
所述相对变化量在数值上等于被测点应变在线段方向上的正应变,其计算表达式为:
式中,为特征线段的初始长度,/>为被测点产生变形后特征线段在T时刻的线段长度;/>为T时刻特征线段的相对变化量。
6.根据权利要求1所述的基于光学显微的一点应变测量方法,其特征在于,
所述一点的应变状态为:ε=[εxyxy];
所述一点应变与任意方向正应变关系为:
7.根据权利要求6所述的基于光学显微的一点应变测量方法,其特征在于,
T时刻各特征线段方向上的正应变用被测点的应变状态分量εx、εy、γxy表示为:
式中,li和mi分别为在初始时刻各线段与x轴和y轴正向夹角的余弦值。
CN202310167033.4A 2023-02-27 2023-02-27 基于光学显微的一点应变测量方法 Active CN115950371B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310167033.4A CN115950371B (zh) 2023-02-27 2023-02-27 基于光学显微的一点应变测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310167033.4A CN115950371B (zh) 2023-02-27 2023-02-27 基于光学显微的一点应变测量方法

Publications (2)

Publication Number Publication Date
CN115950371A CN115950371A (zh) 2023-04-11
CN115950371B true CN115950371B (zh) 2023-10-03

Family

ID=87297665

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310167033.4A Active CN115950371B (zh) 2023-02-27 2023-02-27 基于光学显微的一点应变测量方法

Country Status (1)

Country Link
CN (1) CN115950371B (zh)

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4969106A (en) * 1989-02-27 1990-11-06 Camsys, Inc. Computerized method of determining surface strain distributions in a deformed body
EP0641992A2 (en) * 1993-09-08 1995-03-08 Texas Instruments Incorporated Apparatus systems and methods for measuring defects in semiconductor processing
EP0857944A1 (en) * 1997-02-10 1998-08-12 Xerox Corporation Vibration detection and control system for printers
JP2008020330A (ja) * 2006-07-13 2008-01-31 National Institute Of Advanced Industrial & Technology 光学分光測定方法及び装置
CN104279975A (zh) * 2014-10-17 2015-01-14 中国科学院武汉岩土力学研究所 一种光学显微测量孔径变形的地应力测试方法
CN104807563A (zh) * 2015-05-04 2015-07-29 中国矿业大学 基于钻孔显微数字摄像的非接触式地应力测试装置及方法
CN205785092U (zh) * 2016-01-28 2016-12-07 广州大学 基于超高精度数字图像测量的二维应变测量系统
CN107014304A (zh) * 2017-04-17 2017-08-04 西安交通大学 一种高精度拱坝变形监测设备和测量方法
CN109175179A (zh) * 2018-11-23 2019-01-11 中国航发南方工业有限公司 一种高温合金锻件锻造拔长方法
CN109211153A (zh) * 2018-11-16 2019-01-15 中冶建筑研究总院有限公司 一种结构表面应变的测量方法
CN109855970A (zh) * 2019-04-08 2019-06-07 大连理工大学 一种基于三轴土样全表面变形的应变状态计算方法
CN112212797A (zh) * 2020-10-26 2021-01-12 大连理工大学 一种显微视觉应变传感器
CN112697063A (zh) * 2021-03-23 2021-04-23 四川省人民医院 一种基于显微视觉的芯片应变测量方法
CN114739310A (zh) * 2022-04-19 2022-07-12 大连理工大学 一种基于机器视觉的最大拉压应变检测传感器

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4969106A (en) * 1989-02-27 1990-11-06 Camsys, Inc. Computerized method of determining surface strain distributions in a deformed body
EP0641992A2 (en) * 1993-09-08 1995-03-08 Texas Instruments Incorporated Apparatus systems and methods for measuring defects in semiconductor processing
EP0857944A1 (en) * 1997-02-10 1998-08-12 Xerox Corporation Vibration detection and control system for printers
JP2008020330A (ja) * 2006-07-13 2008-01-31 National Institute Of Advanced Industrial & Technology 光学分光測定方法及び装置
CN104279975A (zh) * 2014-10-17 2015-01-14 中国科学院武汉岩土力学研究所 一种光学显微测量孔径变形的地应力测试方法
CN104807563A (zh) * 2015-05-04 2015-07-29 中国矿业大学 基于钻孔显微数字摄像的非接触式地应力测试装置及方法
CN205785092U (zh) * 2016-01-28 2016-12-07 广州大学 基于超高精度数字图像测量的二维应变测量系统
CN107014304A (zh) * 2017-04-17 2017-08-04 西安交通大学 一种高精度拱坝变形监测设备和测量方法
CN109211153A (zh) * 2018-11-16 2019-01-15 中冶建筑研究总院有限公司 一种结构表面应变的测量方法
CN109175179A (zh) * 2018-11-23 2019-01-11 中国航发南方工业有限公司 一种高温合金锻件锻造拔长方法
CN109855970A (zh) * 2019-04-08 2019-06-07 大连理工大学 一种基于三轴土样全表面变形的应变状态计算方法
CN112212797A (zh) * 2020-10-26 2021-01-12 大连理工大学 一种显微视觉应变传感器
CN112697063A (zh) * 2021-03-23 2021-04-23 四川省人民医院 一种基于显微视觉的芯片应变测量方法
CN114739310A (zh) * 2022-04-19 2022-07-12 大连理工大学 一种基于机器视觉的最大拉压应变检测传感器

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Effect of reference point selection on microscopic stress measurement using EBSD;Mikami, Yoshiki;MATERIALS SCIENCE AND ENGINEERING A-STRUCTURAL MATERIALS PROPERTIES MICROSTRUCTURE AND PROCESSING;256-264 *
叙永地区五峰组-龙马溪组页岩变形及物性演变机理;程国玺;中国博士学位论文全文数据库基础科学辑;A011-131 *
微尺度MEMS上层芯片表面应变测量的方法研究;罗钧文;《机械设计与制造》;196-198 *

Also Published As

Publication number Publication date
CN115950371A (zh) 2023-04-11

Similar Documents

Publication Publication Date Title
CN110689579B (zh) 基于合作目标的快速单目视觉位姿测量方法及测量系统
CN109696356B (zh) 一种土工合成材料拉伸试样全域应变场测量装置及方法
CN107256568B (zh) 一种高精度机械臂手眼相机标定方法及标定系统
JP5082941B2 (ja) 標線位置測定装置、標線位置測定用プログラム、および標線マーク
CN107042528A (zh) 一种工业机器人的运动学标定系统及方法
CN110779797B (zh) 一种测量金属拉伸试验过程中塑性应变比的方法
CN101900531B (zh) 双目视觉位移测量误差的测算方法及测量系统
CN112991460B (zh) 一种用于获取汽车零件尺寸的双目测量系统、方法及装置
CN107907064A (zh) 一种裂隙监测系统及方法
CN109952176A (zh) 一种机器人标定方法、系统、机器人及存储介质
CN115950371B (zh) 基于光学显微的一点应变测量方法
CN102012215B (zh) 基于数字图像的非接触式光学应变测量方法及应变计
CN115540775A (zh) 一种ccd单相机的3d视频引伸计
CN109990801B (zh) 基于铅垂线的水平仪装配误差标定方法
CN111145247B (zh) 基于视觉的位置度检测方法及机器人、计算机存储介质
CN109342008B (zh) 基于单应性矩阵的风洞试验模型迎角单相机视频测量方法
Ehrenfried Processing calibration-grid images using the Hough transformation
Li et al. Normal strain measurement by machine vision
CN114061472B (zh) 基于标靶的测量坐标误差修正的方法
CN111829439B (zh) 一种高精度平移测量方法及装置
Li et al. F-touch sensor for three-axis forces measurement and geometry observation
CN114571506B (zh) 工业机器人性能测量的姿态对齐方法
DE112018005152T5 (de) Druckmesseinrichtung und druckmessverfahren
CN116061196B (zh) 一种多轴运动平台运动学参数标定方法和系统
CN114087988B (zh) 一种单目单线结构光传感器的精度评价方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant