CN114529528B - 一种胶接结构粘接界面演变太赫兹定量计算方法 - Google Patents
一种胶接结构粘接界面演变太赫兹定量计算方法 Download PDFInfo
- Publication number
- CN114529528B CN114529528B CN202210155752.XA CN202210155752A CN114529528B CN 114529528 B CN114529528 B CN 114529528B CN 202210155752 A CN202210155752 A CN 202210155752A CN 114529528 B CN114529528 B CN 114529528B
- Authority
- CN
- China
- Prior art keywords
- image
- change
- terahertz
- bonding
- bonding interface
- 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
- 239000000853 adhesive Substances 0.000 title claims abstract description 18
- 230000001070 adhesive effect Effects 0.000 title claims abstract description 18
- 238000004364 calculation method Methods 0.000 title claims abstract description 12
- 239000010410 layer Substances 0.000 claims abstract description 15
- 238000003384 imaging method Methods 0.000 claims abstract description 14
- 239000012790 adhesive layer Substances 0.000 claims abstract description 13
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 8
- 230000009466 transformation Effects 0.000 claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000000034 method Methods 0.000 claims description 6
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 claims description 4
- 239000013598 vector Substances 0.000 claims description 4
- 238000011426 transformation method Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 2
- 230000007547 defect Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 6
- 239000003292 glue Substances 0.000 description 6
- 238000012512 characterization method Methods 0.000 description 5
- 238000009659 non-destructive testing Methods 0.000 description 5
- 239000002131 composite material Substances 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000004026 adhesive bonding Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000032683 aging Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000036316 preload Effects 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000035882 stress Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3581—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation
- G01N21/3586—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation by Terahertz time domain spectroscopy [THz-TDS]
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Toxicology (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Quality & Reliability (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种胶接结构粘接界面演变太赫兹定量计算方法,使用太赫兹时域光谱仪分别采集加载前后的胶接结构太赫兹时域信号;从胶接结构太赫兹时域信号中提取粘接界面飞行时间特征,并进行成像,获得胶接结构加载前后的胶层粘接界面飞行时间图像I1及图像I2;利用尺度不变特征变换匹配算法获取图像I1与图像I2匹配的特征点;利用图像I1与图像I2识别出的特征点进行图像匹配;将图像I1与图像I2*的重合区域做差值运算,获得胶接结构加载前后的胶层粘接界面的变化图像I3;对变化图像I3整幅图像的变化值分布占比进行统计,对变化图像I3中大于设定阈值的区域进行统计,计算该区域飞行时间变化的均值,用于表征该区域粘接界面的变化情况。
Description
技术领域
本发明属于太赫兹无损检测领域,具体涉及一种利用太赫兹成像对胶接结构粘接界面演变进行定量计算的方法。
背景技术
复合材料胶接结构在使用过程中,存在温度、应力、冲击等环境载荷,尤其当环境载荷大于胶接结构所能承受的范围时,胶接结构的粘接界面可能会出现变化,甚至出现粘接缺陷,导致局部的传力性能发生改变,使得结构失效的风险增加。因此对于粘接界面的识别以及受载后的演变计算,可以对胶接结构粘接质量的评价提供支撑。
目前在无损检测领域对于胶接结构的质量评价多集中于对粘接缺陷的识别,尚无对胶接结构受载后粘接界面的演变表征方法。在太赫兹粘接界面的识别方法方面,CN108267419B中应用太赫兹三维时域数据矩阵的时间轴数据对胶接结构脱粘缺陷进行成像;其成像方法可以识别粘接界面脱粘缺陷的位置。在损伤演变无损表征方面,多集中于采用CT、红外等方法对试件表面缺陷以及裂纹的动态检测,缺少对于胶接结构粘接界面的损伤演变的无损表征手段。
太赫兹无损检测是复合材料无损检测新的一种检测手段,其中太赫兹时域光谱信号的飞行时间特征参数,可表征材料的折射率、材料密度、材料厚度等信息。太赫兹时域光谱信号的飞行时间特征参数也可对粘接结构,尤其是粘接层厚度、缺陷等信息进行表征。
提出了一种基于尺度不变特征变换匹配算法(SIFT)的太赫兹定量表征方法,可以实现对胶接结构粘接界面的变化监测,并对变化情况进行定量计算。
发明内容
针对胶接结构受载后粘接界面的演变问题,本发明提供了一种胶接结构粘接界面演变太赫兹定量计算方法,基于太赫兹无损检测成像技术对粘接界面进行准确识别,利用尺度不变特征变换匹配算法(SIFT)可以实现对胶接结构粘接界面的变化监测,并对变化情况进行定量计算。
本发明的目的是通过以下技术方案实现的:
一种胶接结构粘接界面演变太赫兹定量计算方法,包括以下步骤:
步骤一、使用太赫兹时域光谱仪以二维逐点扫描方式,分别采集加载前后的胶接结构太赫兹时域信号;
步骤二、从采集到的胶接结构太赫兹时域信号中提取粘接界面飞行时间特征,并进行成像,获得胶接结构加载前的胶层粘接界面飞行时间图像I1以及加载后的胶层粘接界面飞行时间图像I2;
步骤三、利用尺度不变特征变换匹配算法(SIFT)获取所述图像I1与图像I2匹配的特征点;
步骤四、利用图像I1与图像I2识别出的特征点进行图像匹配:通过公式I2*=T·I2,将图像I2旋转变换为图像I2*,使图像I2投影在了图像I1上,式中,T为变换矩阵,通过仿射变换方法计算获得;
步骤五、将图像I1与图像I2*的重合区域做差值运算,获得胶接结构加载前后的胶层粘接界面的变化图像I3;
步骤六、对变化图像I3整幅图像的变化值分布占比进行统计,对变化图像I3中大于设定阈值的区域进行统计,计算该区域飞行时间变化的均值,用于表征该区域粘接界面的变化情况。
本发明其技术效果在于:
通过太赫兹无损检测对胶接结构粘接界面进行识别,选取飞行时间成像方式,可以反映粘接界面的纵向信息,以及脱粘缺陷信息,对胶接结构受载后粘接界面的变化反应较为灵敏。
本发明提供的胶接结构粘接界面演变太赫兹定量表征方法,可以有效的将胶接结构粘接界面以及粘接缺陷变化可视化,并可以定量统计粘接界面任一区域内变化情况。
本发明可以推广到太赫兹无损检测对样件内部缺陷变化的应用中,如复合材料受冲击损伤变化监测、绝缘芯子老化过程监测等。
附图说明
图1为本发明一种胶接结构粘接界面演变太赫兹定量计算方法整体流程图;
图2为本发明实施例所述的胶接结构样件示意图;图中,1-复合材料、2-上层胶层、3-应变隔离垫、4-下层胶层、5-金属基底;
图3为本发明实施例所述胶接结构样件某点的太赫兹时域波形;
图4(a)为本发明实施例所述的加载前胶接结构样件下层胶层飞行时间成像I1;
图4(b)为本发明实施例所述的加载后胶接结构样件下层胶层飞行时间成像I2;
图5(a)为本发明实施例所述的图像I1与I2特征点识别图;
图5(b)为本发明实施例所述的图像I1与I2剔除匹配误差大的特征点后的特征点识别图;
图6为本发明实施例所述图像I1与图像I2匹配结果示意图;
图7为本发明实施例所述胶接结构加载前后的胶层粘接界面的变化图像I3;
图8(a)为本发明实施例所述粘接界面的变化分布直方图;
图8(b)为本发明实施例所述变化图像I3中大于设定阈值的区域示意图;
图8(c)为本发明实施例所述大于设定阈值的区域的粘接界面局部的变化情况。
具体实施方式
为了使发明的目的、技术方案及优点更加地清楚明白,以下结合附图,对本发明进行进一步详细说明。
一种胶接结构粘接界面演变太赫兹定量计算方法,本实施例的胶接结构如图2所示,本发明具体方案如下所述:
步骤一、利用太赫兹时域光谱仪以二维逐点扫描方式,采集加载前后的胶接结构样件太赫兹时域信号,设置行列采集步长分别为a,b,采集步距为1mm,获得样件加载前后的太赫兹数据集合A1、A2。
步骤二、从采集到的胶接结构样件太赫兹时域信号中提取粘接界面特征信息,如图3所示,为获取到的太赫兹数据集中某点(a,b)的太赫兹时域信号Sa,b(t),横轴为飞行时间,纵轴为幅值,其中该点处下层胶层特征波形为b所示区域,计算其飞行时间差作为粘接界面纵向信息的特征值。飞行时间的计算公式为:
其中,tfly为图3波形中下层胶层上表面反射信号Peak1与金属基底反射信号Peak2的飞行时间差tfly=Peak2-Peak1,n为有机胶粘接层的折射率,d为有机胶粘接层的厚度,c为光传播的速度。
通过计算每一点的胶层飞行时间差,获取采集数据集合A1、A2的下层胶层飞行时间值。通过飞行时间成像就可以得到胶接结构受载前后胶层的二维图像信息,如图4(a)为受载前样件下层胶层飞行时间成像图I1,图4(b)为受载后样件下层胶层飞行时间成像图I2,此方法可以将样件的粘接界面可视化。
步骤三、将步骤二获取到的加载前后飞行时间成像图I1、I2,利用尺度不变特征变换匹配算法进行特征点匹配。首先构建图像I1、I2的尺度空间。将图像的多尺度空间定义为一个函数L(x,y,σ),它是由一个可以改变尺度的高斯函数G(x,y,σ)与一个输入图像I(x,y)卷积得到的,如公式(2)所示:
L(x,y,σ)=G(x,y,σ)*I(x,y) (2)
其中(x,y)为坐标空间,σ为可变比例的空间因子。在得到的尺度空间中寻求极值点作为图像的待选特征点。在特征点处提取特征描述符,形成特征向量。当两幅图像的特征向量计算完成后,采用关键点特征向量的欧式距离来作为两幅图像中关键点的相似性判定度量,寻找两幅图像中匹配度较高的特征点,在上诉两幅图像中画上标记,连接对应的匹配的点。结果如图5(a)所示,图像中线条所连接的点为识别出的两幅图像特征相同的对应点。
由于特征空间的高维性,相似的距离可能有大量其他的错误匹配,为了得到准确的图像转换模型,采用随机样本一致性算法(RANSAC),将匹配误差大的特征点进行剔除,获得结果如图5(b)所示。
步骤四、将图像I2中的特征点作为移动点,图像I1中的特征点作为固定点,利用仿射变换,求取变化矩阵,计算公式如(3)所示:
其中图像I1的特征点为(x,y),图像I2的特征点为(u,v),T为变换矩阵。将求取得到的变换矩阵T,利用公式(4)得到I2*,其为图像I2旋转变换后的图像,此时图像I2*图像I1重合,如图6所示。
I2*=T·I2 (4)
步骤五、利用公式(5),获取图像I1与I2*重合区域的差值,得到受载后粘接界面太赫兹飞行时间特征值变化图像I3。结果如图7所示,获得了受载前后粘接界面的变化分布情况。
I3=I2*-I1 (5)
步骤六、在获取了粘接界面的变化图像I3后,对粘接界面变化值进行统计分析,得到该区域内飞行时间变化分布直方图,如图8(a)所示,统计得到不同飞行时间变化区间内所占的比例,从而判断出粘接界面整体的变化情况。同时可以对图像I3中大于0.5ps的区域进行统计,如图8(b)中黑色框内的区域,获取该区域粘接界面局部的变化情况如图8(c)直方图所示,计算其变化均值,表征变化的大小。
Claims (2)
1.一种胶接结构粘接界面演变太赫兹定量计算方法,其特征在于,包括以下步骤:
步骤一、使用太赫兹时域光谱仪以二维逐点扫描方式,分别采集加载前后的胶接结构太赫兹时域信号;
步骤二、从采集到的胶接结构太赫兹时域信号中提取粘接界面飞行时间特征,并进行成像,获得胶接结构加载前的胶层粘接界面飞行时间图像I1以及加载后的胶层粘接界面飞行时间图像I2;
步骤三、利用尺度不变特征变换匹配算法(SIFT)获取所述图像I1与图像I2匹配的特征点;所述步骤三包括:
首先构建所述图像I1及图像I2的尺度空间:将图像的多尺度空间定义为一个函数L(x,y,σ),它是由一个可以改变尺度的高斯函数G(x,y,σ)与一个输入图像I(x,y)卷积得到的,如下式所示:
L(x,y,σ)=G(x,y,σ)*I(x,y)
其中,(x,y)为坐标空间,σ为可变比例的空间因子;
在得到的尺度空间中寻求极值点作为图像的待选特征点;
在特征点处提取特征描述符,形成特征向量;
最后,采用随机样本一致性算法,将匹配误差大的特征点进行剔除;
步骤四、利用图像I1与图像I2识别出的特征点进行图像匹配:通过公式I2*=T·I2,将图像I2旋转变换为图像I2*,使图像I2投影在了图像I1上,式中,T为变换矩阵,通过仿射变换方法计算获得;
步骤五、将图像I1与图像I2*的重合区域做差值运算,获得胶接结构加载前后的胶层粘接界面的变化图像I3;
步骤六、对变化图像I3整幅图像的变化值分布占比进行统计,对变化图像I3中大于设定阈值的区域进行统计,计算该区域飞行时间变化的均值,用于表征该区域粘接界面的变化情况。
2.如权利要求1所述的一种胶接结构粘接界面演变太赫兹定量计算方法,其特征在于,所述步骤四中,仿射变换方法计算变化矩阵T包括:
将所述图像I2中的特征点作为移动点,图像I1中的特征点作为固定点,通过下式求取变化矩阵:
其中,图像I1的特征点为(x,y),图像I2的特征点为(u,v),T为变换矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210155752.XA CN114529528B (zh) | 2022-02-21 | 2022-02-21 | 一种胶接结构粘接界面演变太赫兹定量计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210155752.XA CN114529528B (zh) | 2022-02-21 | 2022-02-21 | 一种胶接结构粘接界面演变太赫兹定量计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114529528A CN114529528A (zh) | 2022-05-24 |
CN114529528B true CN114529528B (zh) | 2023-08-29 |
Family
ID=81625745
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210155752.XA Active CN114529528B (zh) | 2022-02-21 | 2022-02-21 | 一种胶接结构粘接界面演变太赫兹定量计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114529528B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108267419A (zh) * | 2017-12-08 | 2018-07-10 | 山东省科学院自动化研究所 | 太赫兹时域光谱检测复合材料胶接结构脱粘缺陷的方法 |
CN110230987A (zh) * | 2019-05-31 | 2019-09-13 | 四川省派瑞克斯光电科技有限公司 | 一种基于太赫兹时域光谱成像对橡胶复合材料的无损检测方法 |
CN111272690A (zh) * | 2019-12-31 | 2020-06-12 | 长春理工大学 | 有机胶粘剂粘接固化特性的太赫兹表征方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110599558B (zh) * | 2018-06-12 | 2021-06-29 | 清华大学 | 太赫兹光谱成像的数据重构方法、装置、设备和存储介质 |
-
2022
- 2022-02-21 CN CN202210155752.XA patent/CN114529528B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108267419A (zh) * | 2017-12-08 | 2018-07-10 | 山东省科学院自动化研究所 | 太赫兹时域光谱检测复合材料胶接结构脱粘缺陷的方法 |
CN110230987A (zh) * | 2019-05-31 | 2019-09-13 | 四川省派瑞克斯光电科技有限公司 | 一种基于太赫兹时域光谱成像对橡胶复合材料的无损检测方法 |
CN111272690A (zh) * | 2019-12-31 | 2020-06-12 | 长春理工大学 | 有机胶粘剂粘接固化特性的太赫兹表征方法 |
Non-Patent Citations (1)
Title |
---|
Terahertz-Light Induced Structural Transition and Superpermeation of Confined Monolayer Water;Zhi Zhu等;《ACS Photonics 2021》;第781-786页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114529528A (zh) | 2022-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108074231B (zh) | 一种基于卷积神经网络的磁片表面缺陷检测方法 | |
Marani et al. | Modeling and classification of defects in CFRP laminates by thermal non-destructive testing | |
CN110045015B (zh) | 一种基于深度学习的混凝土结构内部缺陷检测方法 | |
Marani et al. | Automatic detection of subsurface defects in composite materials using thermography and unsupervised machine learning | |
CN109919905B (zh) | 一种基于深度学习的红外无损检测方法 | |
CN103940708B (zh) | 一种钢材全形态晶粒的快速测量、精细分类方法 | |
CN111882554B (zh) | 一种基于SK-YOLOv3的电力线故障智能检测方法 | |
Li et al. | Image quality assessment using deep convolutional networks | |
CN115035081B (zh) | 一种基于工业ct的金属内部缺陷危险源定位方法及*** | |
CN111272690B (zh) | 有机胶粘剂粘接固化特性的太赫兹表征方法 | |
CN113838010A (zh) | 基于红外热波及卷积神经网络的复合材料损伤智能检测方法 | |
Wei et al. | Feature Extraction and Retrieval of Ecommerce Product Images Based on Image Processing. | |
CN114529528B (zh) | 一种胶接结构粘接界面演变太赫兹定量计算方法 | |
CN113537236B (zh) | 用于航天器损伤红外检测的热扩散效应缺陷定量识别方法 | |
CN117495852B (zh) | 一种基于图像分析的数字印刷质量检测方法 | |
CN109145770B (zh) | 一种基于多尺度特征融合网络与定位模型相结合的麦蜘蛛自动计数方法 | |
Gao et al. | Estimation of the convolutional neural network with attention mechanism and transfer learning on wood knot defect classification | |
CN117406146A (zh) | 一种基于深度学习的磁场检测方法 | |
CN111948171A (zh) | 基于动态时间规整的复合材料缺陷识别方法 | |
CN109816651B (zh) | 基于变化率与温度差异的热图像缺陷特征提取方法 | |
Liu et al. | Multi-feature vision transformer for automatic defect detection and quantification in composites using thermography | |
CN116047418A (zh) | 基于小样本的多模态雷达有源欺骗干扰识别方法 | |
CN112580568B (zh) | 基于高斯核均值局部线性嵌入的芯片缺陷检测方法 | |
CN109886930B (zh) | 基于变化率与温度差异的热图像缺陷特征提取方法 | |
CN112102394A (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 |