CN115950371A - 基于光学显微的一点应变测量方法 - Google Patents
基于光学显微的一点应变测量方法 Download PDFInfo
- Publication number
- CN115950371A CN115950371A CN202310167033.4A CN202310167033A CN115950371A CN 115950371 A CN115950371 A CN 115950371A CN 202310167033 A CN202310167033 A CN 202310167033A CN 115950371 A CN115950371 A CN 115950371A
- Authority
- CN
- China
- Prior art keywords
- point
- line segment
- strain
- 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.)
- Granted
Links
- 238000000399 optical microscopy Methods 0.000 title claims abstract description 25
- 238000000691 measurement method Methods 0.000 title claims abstract description 18
- 238000005259 measurement Methods 0.000 claims abstract description 12
- 230000014509 gene expression Effects 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 10
- 230000003287 optical effect Effects 0.000 claims description 5
- 239000003086 colorant Substances 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 2
- 238000012360 testing method Methods 0.000 abstract description 8
- 239000007788 liquid Substances 0.000 abstract description 4
- 230000007613 environmental effect Effects 0.000 abstract description 2
- 230000008859 change Effects 0.000 description 3
- 239000000853 adhesive Substances 0.000 description 2
- 230000001070 adhesive effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种基于光学显微的一点应变测量方法,包括:基于光学显微设备对介质被测点进行高倍率放大摄像或拍照,获得被测点局部范围高像素细观特征图像,在图像中选择若干个特征点,将特征点相连获得特征线段;测量特征线段在初始时刻的线段长度及其与坐标轴正方向夹角,计算夹角的余弦值;测量特征线段在后续某一时刻的线段长度,计算特征线段相比于初始时刻长度的相对变化量;基于一点的应变状态与该点任意方向正应变之间的关系式计算被测点的应变,获得应变测量结果。本发明解决了电阻式应变测量方法应变片覆盖面积大、环境适应性差、测试精度低、不适用用导电液体环境等问题。
Description
技术领域
本发明属于应变测量领域,特别是涉及一种基于光学显微的一点应变测量方法。
背景技术
一点的应变测量通常由电阻式应变花进行测量,每个应变花包含不少于3个不同方向的应变片,为了保证测量精度电阻应变片的长度普遍大于1cm。因此,每个应变片测得的应变应为应变片覆盖范围沿应变片长度方向的平均正应变,进而由应变计算得到的一点应变近似为应变花覆盖区域的平均应变值,因此,电阻应变测量方法无法精确获得微小区域的真实应变。
由于电阻式应变测量需将应变片与被测介质采用粘结剂进行紧密粘贴,介质在变形过程中一定程度上会受到粘结剂和电阻式应变片的束缚,其并非是完全自由变形状态,尤其是测量对象为质地较软时,这种影响将非常显著,导致错误的测量结果。若介质处于潮湿环境或介质不致密时,电阻式应变片粘贴效果难以保障,进而产生错误的测试数据,甚至致使测试失败。当被测介质浸于导电的液体中时(如充水钻孔、海上轮船底部等),电阻式应变无法进行直接测量,必须临时封隔出干燥腔体进行测量,给测试工作带来的很大的困难。即便如此,由于腔体内介质不受液体压力作用,其应变状态与正常工作条件时并不相同。电阻式应变花为一次性测量元件,大量测量时其成本相对较大。
发明内容
为解决上述问题,本发明提供了如下方案:基于光学显微的一点应变测量方法,包括:
基于光学显微设备对介质被测点进行高倍率放大摄像或拍照,获得被测点局部范围高像素细观特征图像,在图像中选择若干个特征点,将所述特征点相连获得特征线段,基于所述特征线段获得特征线段集;
测量所述特征线段在初始时刻的线段长度及所述特征线段与坐标轴正方向夹角,计算夹角的余弦值;
基于所述特征线段在后续某一时刻的线段长度与所述特征线段在初始时刻的线段长度,计算特征线段相比于初始时刻长度的相对变化量;
基于所述特征线段长度的相对变化量、一点的应变状态与该点任意方向正应变之间的关系,计算被测点的应变,获得应变测量结果。
优选地,所述光学显微设备的放大倍率在1倍~100000倍之间连续可调。
优选地,所述特征图像像素不低于100万。
优选地,所述特征线段的长度及所述特征线段与坐标轴正方向之间的夹角基于所述光学显微设备或图像分析软件的虚拟标尺、标记和像素点位置识别功能实现。
优选地,所述特征点为与被测介质固接且同步变形的可分辨物理特征点,至少包括尖点、凹点、散斑、颜色。
优选地,所述特征线段为从一个特征点中心位置至另一个特征点中心位置的直线段。
优选地,所述若干个特征点为特征点数量满足两两相连所得的特征线段中不平行线段数量不少于3条;
所述特征线段集为根据需要连接任意两个特征点得到的特征线段的集合,且不平行线段数量不少于3条。
优选地,所述相对变化量在数值上等于被测点应变在线段方向上的正应变,其计算表达式为:
优选地,所述一点的应变状态为:ε=[εx,εy,γxy];
所述一点应变与任意方向正应变关系为:
优选地,所述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]。
进一步地优化方案,所述特征线段长度相比于初始时刻的相对变化量指线段当前长度与其初始长度的差值与线段初始长度的比值。
设一点的应变状态为ε=[εx,εy,γxy],则该点在任意方向的正应变为:
采用光学显微方法对被测点进行放大摄像或拍照,获得被测点高像素特征图像并选择待测目标微区域内的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 (10)
1.一种基于光学显微的一点应变测量方法,其特征在于,包括:
基于光学显微设备对介质被测点进行高倍率放大摄像或拍照,获得被测点局部范围高像素细观特征图像,在图像中选择若干个特征点,将所述特征点相连获得特征线段,基于所述特征线段获得特征线段集;
测量所述特征线段在初始时刻的线段长度及所述特征线段与坐标轴正方向夹角,计算夹角的余弦值;
基于所述特征线段在后续某一时刻的线段长度与所述特征线段在初始时刻的线段长度,计算特征线段相比于初始时刻长度的相对变化量;
基于所述特征线段长度的相对变化量、一点的应变状态与该点任意方向正应变之间的关系,计算被测点的应变,获得应变测量结果。
2.根据权利要求1所述的基于光学显微的一点应变测量方法,其特征在于,所述光学显微设备的放大倍率在1倍~100000倍之间连续可调,特征图像像素不低于100万。
3.根据权利要求1所述的基于光学显微的一点应变测量方法,其特征在于,所述特征线段的长度及所述特征线段与坐标轴正方向之间的夹角基于所述光学显微设备或图像分析软件的虚拟标尺、标记和像素点位置识别功能实现。
4.根据权利要求1所述的基于光学显微的一点应变测量方法,其特征在于,所述特征点为与被测介质固接且同步变形的可分辨物理特征点,至少包括尖点、凹点、散斑、颜色。
5.根据权利要求1所述的基于光学显微的一点应变测量方法,其特征在于,所述特征线段为从一个特征点中心位置至另一个特征点中心位置的直线段。
6.根据权利要求1所述的基于光学显微的一点应变测量方法,其特征在于,
所述若干个特征点为特征点数量满足两两相连所得的特征线段中不平行线段数量不少于3条;
所述特征线段集为根据需要连接任意两个特征点得到的特征线段的集合,且不平行线段数量不少于3条。
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 true CN115950371A (zh) | 2023-04-11 |
CN115950371B 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)
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 | 大连理工大学 | 一种基于机器视觉的最大拉压应变检测传感器 |
-
2023
- 2023-02-27 CN CN202310167033.4A patent/CN115950371B/zh active Active
Patent Citations (14)
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)
Title |
---|
MIKAMI, YOSHIKI: "Effect of reference point selection on microscopic stress measurement using EBSD", MATERIALS SCIENCE AND ENGINEERING A-STRUCTURAL MATERIALS PROPERTIES MICROSTRUCTURE AND PROCESSING, pages 256 - 264 * |
程国玺: "叙永地区五峰组-龙马溪组页岩变形及物性演变机理", 中国博士学位论文全文数据库基础科学辑, pages 011 - 131 * |
罗钧文: "微尺度MEMS上层芯片表面应变测量的方法研究", 《机械设计与制造》, pages 196 - 198 * |
Also Published As
Publication number | Publication date |
---|---|
CN115950371B (zh) | 2023-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107042528B (zh) | 一种工业机器人的运动学标定系统及方法 | |
US6694634B2 (en) | Position error evaluation method of moving device and method of improving movement accuracy based on evaluation results thereof | |
CN110948522B (zh) | 一种基于拉线旋转传感器的工业机器人空间位姿测量机构及测量方法 | |
US8502871B2 (en) | Gauge line position measuring device, program for measuring a gauge line position, and gauge line mark | |
CN101900531B (zh) | 双目视觉位移测量误差的测算方法及测量系统 | |
CN109099883A (zh) | 高精度大视场机器视觉测量与标定装置及方法 | |
CN108801166A (zh) | 基于悬臂梁理论的光纤光栅机翼形变测量建模及标定方法 | |
CN104729481A (zh) | 一种基于pnp透视模型的合作目标位姿精度测量方法 | |
CN108627104A (zh) | 一种零部件高度尺寸的点激光测量方法 | |
CN105865349A (zh) | 一种大型建筑物位移监测方法 | |
CN114184931A (zh) | 探针调整方法、装置、电子设备和存储介质 | |
CN112792814A (zh) | 基于视觉标志的机械臂零点标定方法 | |
CN106671081B (zh) | 一种基于单目视觉的少自由度机器人运动学标定方法 | |
CN109342008B (zh) | 基于单应性矩阵的风洞试验模型迎角单相机视频测量方法 | |
CN105865401B (zh) | 用于可伸缩臂架的形变量检测方法和装置 | |
CN109990801B (zh) | 基于铅垂线的水平仪装配误差标定方法 | |
CN115950371A (zh) | 基于光学显微的一点应变测量方法 | |
CN111145247B (zh) | 基于视觉的位置度检测方法及机器人、计算机存储介质 | |
CN110672094B (zh) | 一种分布式pos多节点多参量瞬间同步标校方法 | |
CN116216631A (zh) | 一种高温硅压力传感器的点胶装配方法 | |
CN110954081A (zh) | 一种磁罗盘快速校准装置及方法 | |
CN114061472B (zh) | 基于标靶的测量坐标误差修正的方法 | |
CN115674271A (zh) | 一种基于多工位测量的机器人标定方法 | |
DE112018005152T5 (de) | Druckmesseinrichtung und druckmessverfahren | |
KR100735494B1 (ko) | 방위각 측정장치 및 그 측정방법 |
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 |