CN113406039A - 一种太赫兹主动式风力叶片泡沫芯材无损检测方法及装置 - Google Patents
一种太赫兹主动式风力叶片泡沫芯材无损检测方法及装置 Download PDFInfo
- Publication number
- CN113406039A CN113406039A CN202110644852.4A CN202110644852A CN113406039A CN 113406039 A CN113406039 A CN 113406039A CN 202110644852 A CN202110644852 A CN 202110644852A CN 113406039 A CN113406039 A CN 113406039A
- Authority
- CN
- China
- Prior art keywords
- terahertz
- core material
- foam core
- image
- defect
- 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.)
- Pending
Links
- 239000011162 core material Substances 0.000 title claims abstract description 81
- 239000006260 foam Substances 0.000 title claims abstract description 79
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000009659 non-destructive testing Methods 0.000 title claims description 25
- 230000007547 defect Effects 0.000 claims abstract description 99
- 238000001514 detection method Methods 0.000 claims abstract description 53
- 238000003384 imaging method Methods 0.000 claims abstract description 16
- 238000002592 echocardiography Methods 0.000 claims abstract description 6
- 238000003325 tomography Methods 0.000 claims abstract description 4
- 238000003709 image segmentation Methods 0.000 claims description 18
- 239000000463 material Substances 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 11
- 239000002184 metal Substances 0.000 claims description 9
- 229910052751 metal Inorganic materials 0.000 claims description 9
- 238000001914 filtration Methods 0.000 claims description 6
- 238000005192 partition Methods 0.000 claims description 6
- 239000000758 substrate Substances 0.000 claims description 6
- 239000013077 target material Substances 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 3
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 3
- 239000000919 ceramic Substances 0.000 claims description 3
- 239000002131 composite material Substances 0.000 claims description 3
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims description 3
- 239000010931 gold Substances 0.000 claims description 3
- 229910052737 gold Inorganic materials 0.000 claims description 3
- 229910002804 graphite Inorganic materials 0.000 claims description 3
- 239000010439 graphite Substances 0.000 claims description 3
- 239000004033 plastic Substances 0.000 claims description 3
- 229920003023 plastic Polymers 0.000 claims description 3
- 229920000642 polymer Polymers 0.000 claims description 3
- 238000003672 processing method Methods 0.000 claims description 3
- 238000011002 quantification Methods 0.000 claims description 3
- 229910052710 silicon Inorganic materials 0.000 claims description 3
- 239000010703 silicon Substances 0.000 claims description 3
- 239000011152 fibreglass Substances 0.000 claims description 2
- 238000012804 iterative process Methods 0.000 claims description 2
- 229910052755 nonmetal Inorganic materials 0.000 claims description 2
- 238000011144 upstream manufacturing Methods 0.000 claims description 2
- 230000002950 deficient Effects 0.000 abstract description 4
- 239000002699 waste material Substances 0.000 abstract description 4
- 238000010248 power generation Methods 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 3
- 230000001066 destructive effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 241001391944 Commicarpus scandens Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
Images
Classifications
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Toxicology (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提供一种太赫兹主动式风力叶片泡沫芯材无损检测方法及装置，属于太赫兹检测技术领域，该方法是利用连续太赫兹波检测系统的太赫兹发射源对样品风力叶片泡沫芯材发射太赫兹电磁波，根据被泡沫芯材透射和/或反射后的太赫兹波强度得到样品的太赫兹图像；连续太赫兹波检测系统对风力叶片芯材进行扫描式成像；连续太赫兹波检测系统的一次二维扫描可以同时采集风力叶片泡沫芯材不同深度的太赫兹反射回波，形成层析成像，实现太赫兹无损检测。本发明利用连续太赫兹波通过对泡沫芯材的高分辨率成像分析各种缺陷的成像特征进而有效的分辨缺陷泡沫芯材，更加直观、更易于辨别缺陷种类缺陷位置，提高检测效率的同时也减少对泡沫芯材资源的浪费。
Description
技术领域
本发明涉及太赫兹检测技术领域，具体地说是一种太赫兹主动式风力叶片泡沫芯材无损检测方法及装置。
背景技术
一般的，随着风力发电投资的逐年增加，风力发电装备的制备研发逐年增强，但是实际的风力发电叶片在制备研发、使用过程中发现的问题也尤为突出。风力发电叶片泡沫芯材由于加工过程中不可控因素的影响导致裂缝、空隙、空心等缺陷的产生；带有缺陷泡沫芯材的风力叶片在实际运行中由于持续受到风力的多方向剪切力作用极易发生断裂、变形甚至发生整个发电机组的倾倒毁损。
现有的X射线检测技术不适用于泡沫芯材缺陷的难点，因为X射线属于高能射线，对于密度相对较低的泡沫芯材会完全穿透，不产生任何成像效果。
现有的超声无损检测技术不能准确提供缺陷的种类、缺陷的定位时间慢等问题不能有效指导缺陷的查找、修补工作。
现有的其他检测或探测技术对风力叶片泡沫芯材造成破坏性影响的主要是较大的缺陷。如何在保证风力叶片泡沫芯材不被破坏的前提下对风力叶片泡沫芯材进行无损检测是现阶段需要克服的重要课题。
发明内容
本发明的技术任务是解决现有技术的不足，提供一种太赫兹主动式风力叶片泡沫芯材无损检测方法及装置。
本发明的技术方案是按以下方式实现的，本发明的一种太赫兹主动式风力叶片泡沫芯材无损检测方法，该方法是利用连续太赫兹波检测系统的太赫兹发射源对样品风力叶片泡沫芯材发射太赫兹电磁波，根据被泡沫芯材透射和/或反射后的太赫兹波强度得到样品的太赫兹图像；
连续太赫兹波检测系统对风力叶片芯材进行扫描式成像；连续太赫兹波检测系统使太赫兹波聚焦在风力叶片泡沫芯材上，尤其粘接层处，风力叶片泡沫芯材会对太赫兹波产生损耗，而粘接缺陷的边缘对太赫兹波的散射效应，会影响太赫兹波的强度分布，获得层析图像，尤其粘接层处的层析图像，根据焦平面附近的图像，可直观判断出粘接面有无缺陷；
连续太赫兹波检测系统的一次二维扫描可以同时采集风力叶片泡沫芯材不同深度的太赫兹反射回波，形成层析成像，实现太赫兹无损检测。
得到泡沫芯材太赫兹图像后，将所得原始RGB图像转换为相应的灰度图像，并进行中值滤波，消除孤立的噪声点；根据连续太赫兹图像中缺陷的灰度值分布特点，选择图像分割算法对其进行分割，从而提取缺陷前景，实现对缺陷面积定量。
通过将强度信号线性量化成RGB图像可以得出样品物体内部的形状、缺陷及边缘位置；不同频率的太赫兹波被透镜汇聚到样件的不同深度，太赫兹波穿过泡沫芯材与粘接层后被金属基板反射，在混频器中与入射信号进行混频并被太赫兹探测器接收，其中含有被测量的相对距离信息，经过锁相放大器的信号放大传输至上位机进行图像处理。
连续太赫兹波检测系统的工作运算是：
信号发射接收所用时间为τ，
太赫兹源发射的初始频率f0，
带宽B，
周期T，
调频斜率μ＝B/T，
设太赫兹发射头到目标的距离为R，
ct为太赫兹波在空气中的传播速度，则τ＝2R/ct，
分别用实线和虚线表示太赫兹源发射频率ft和经金属基板反射后太赫兹波频率fr，
f0为ft的初始值，
f1则为fr的终值，
t为检测时间，
发射信号频率为:ft＝f0+μt 式(1)
反射信号频率为:fr＝ft(t－τ)＝f0+μ(t－τ) 式(2)
混频器的作用是使发射信号频率与反射信号频率相减得到中频信号频率fif，表示为:fif＝ft－fr＝μτ＝2BR/(TcT) 式(3)
于是有:R＝TcTfif/(2B) 式(4)
频带宽B和周期T在特定情况下，混频器中频信号输出fif与到材料的距离R成正比，根据中频信号fif可得出太赫兹透镜到材料的距离，中频信号的振幅反映了被测材料的折射率，通过混频器发出的差拍信号，即可实现目标材料的振幅与相位成像；锁相放大器将信号放大后交送上位机，再经由图像处理后，可显示Z厚度方向0～1000mm纵深范围内任一深度样件的方向X-Y二维截面太赫兹图像。
对于所获得的太赫兹图像，结合中值滤波，采用模糊c均值FCM聚类算法对缺陷图像进行处理，分割得到含有目标缺陷信息的二值图像，定位目标缺陷信息所示的缺陷位置，并计算像素点得到缺陷面积。
对于缺陷图像处理方法，分为[1]缺陷图像分割和[2]缺陷面积估算检测两步；
[1]缺陷图像分割
缺陷图像分割算法通过迭代运算实现目标函数：类内误差平方和函数的不断优化，在图像分割的应用就是根据图像中像素和c个聚类中心的每一个中心间的加权隶属度，对目标函数进行多次迭代优化，得到聚类目标函数数J最小时的聚类中心矩阵P与模糊划分矩阵U；
uij表示第j个样本属于第i个聚类的隶属度，U由uij构成；在模糊聚类算法中，
其聚类目标函数J为：
式(5)中:
n为图像中总像素个数；
m为加权指数，决定分类模糊程度；
dij＝||pi-xj||是样本xj与聚类中心pi两者间的欧式距离，所有pi构成P；
为约束
引入拉格朗日乘子λ，则聚类目标函数可以写成：
为得到聚类目标函数J最小时的模糊划分矩阵U，对目标函数J求关于λ和uij的偏导数，令其为0，则有:
由式(7)、式(8)得：
令目标函数关于聚类中心pi的偏导数为0，则有:
最终得到pi的更新公式：
该算法的主要实现步骤为：
(4)确定聚类中心个数c和加权指数m，设定迭代阈值ε＞0，初始迭代次数为0，并初始化聚类中心矩阵P；
(5)计算更新U(N),P(N)；
(6)若||P(N+1)-P(N)||＜ε停止迭代，否则将置迭代系数N＝N+1，并返回步骤(2)继续迭代；
针对风力叶片泡沫芯材缺陷的连续太赫兹图像，目标缺陷即为一类；该算法的主要思想：评价聚类性能的准则函数在多次迭代过程后达到最优，并将数据集划分为不同的聚类，从而使生成的每个聚类达到内紧凑与类间独立；粘接层中存在的缺陷使反射回波的强度存在差异，在连续太赫兹波图像上，粘接缺陷体现为一种单独的聚类，像素聚类与图像分割一致，从而将缺陷图像与其它类分隔开；
[2]缺陷面积估算
将粘接层太赫兹图像进行图像分割得到含有缺陷信息的二值图像，遍历图像中所有缺陷像素点，即可求出缺陷类像素点的总和；
由于连续太赫兹系统的的X轴步进精度和Y轴的步进精度均为1mm，则每个像素点的所代表的真实面积为1mm2；所以缺陷类像素点总和的数值即为缺陷面积的数值。
连续太赫兹波检测系统的太赫兹发射源发射的太赫兹波采用频率范围为0.3～1.5THz的反射式连续太赫兹波；能够通过频率源、信号源的整合发射稳定频率范围的太赫兹信号，太赫兹波功率300mW～5W；
连续太赫兹波检测系统的太赫兹探测器，其探测频率范围0.1～2.0THZ。
风力叶片芯材采用陶瓷、石墨、高分子复合材料、塑料或泡沫的非金属芯材。风力发电叶片使用的泡沫主要是玻璃钢夹芯结构的芯材。
一种实现太赫兹主动式风力叶片泡沫芯材无损检测方法的装置，其结构是：
在太赫兹发射源的发射轴向上配置有喇叭天线、混频器、太赫兹透镜、样品风力叶片泡沫芯材、反射板，
于太赫兹发射轴向上在太赫兹透镜的上游配置有混频器，
在混频器与发射轴交点处的垂直方向配置有太赫兹探测器，
太赫兹探测器连接锁相放大器，锁相放大器连接上位机；
以上构成连续太赫兹波检测系统。
连续太赫兹波检测系统的太赫兹发射源发射的太赫兹波采用频率范围为0.3～1.5THz的反射式连续太赫兹波；能够通过频率源、信号源的整合发射稳定频率范围的太赫兹信号，太赫兹波功率300mW～5W；
连续太赫兹波检测系统的太赫兹探测器，其探测频率范围0.1～2.0THZ；
喇叭天线材质采用金属金，发射角度范围-60°～60°；
混频器：使发射信号频率与反射信号频率相减得到中频信号，中频信号的振幅反映了被测材料的折射率，通过混频器发出的差拍信号，即可实现目标材料的振幅与相位成像；
锁相放大器是从干扰环境中分离出特定载波频率信号的放大器；
太赫兹透镜采用高纯硅太赫兹透镜，直径40cm，焦距30cm。
本发明与现有技术相比所产生的有益效果是：
本发明的一种太赫兹主动式风力叶片泡沫芯材无损检测方法及装置利用太赫兹主动式风力叶片泡沫芯材无损检测装置对泡沫芯材进行扫描式成像；太赫兹检测泡沫芯材的分辨率为10-20μm，通过对泡沫芯材的高分辨率成像分析各种缺陷的成像特征进而有效的分辨缺陷泡沫芯材，提高风力发电叶片运行的整体安全性。
太赫兹主动式风力叶片泡沫芯材无损检测装置因为采用二维图像成像的方式相对于超声波无损检测技术更加直观、更易于辨别缺陷种类缺陷位置，提高检测效率的同时也减少对泡沫芯材资源的浪费。
本发明能够通过确定风力发电泡沫芯材中的缺陷位置进行泡沫芯材缺陷的检测，提高了风力叶片泡沫芯材的检测效率；同时也提升了风力发电叶片整体的安全性能。通过将缺陷种类、位置信息反馈给缺陷处理单元，实现了缺陷泡沫板材的重复利用，减少了资源的浪费。
本发明的一种太赫兹主动式风力叶片泡沫芯材无损检测方法及装置设计合理、结构简单、安全可靠、使用方便、易于维护，具有很好的推广使用价值。
附图说明
附图1是本发明的检测装置的结构示意图；
附图2是本发明的连续太赫兹波检测系统的锯齿波调制连续太赫兹发射与回波时频曲线示意图。
附图中的标记分别表示：
1、太赫兹发射源，2、喇叭天线，3、混频器，4、太赫兹探测器，5、锁相放大器，6、上位机，7、太赫兹透镜，8、风力叶片泡沫芯材，9、反射板。
具体实施方式
下面结合附图对本发明的一种太赫兹主动式风力叶片泡沫芯材无损检测方法及装置作以下详细说明。
如附图所示，本发明的一种太赫兹主动式风力叶片泡沫芯材无损检测方法及装置，采用频率范围为0.3～1.5THz的反射式连续太赫兹波检测系统，对风力叶片泡沫芯材进行了无损检测。使太赫兹波聚焦在粘接层处，获得层析图像，根据焦平面附近的图像，可直观判断出粘接面有无缺陷。结合中值滤波，采用模糊c均值(FCM)聚类算法对缺陷图像进行分割得到含有目标缺陷信息的二值图像，计算像素点得到缺陷面积。
发射频率范围为0.1～10THz的太赫兹发射源发射的太赫兹电磁波，不能穿透金属材料，但可穿透陶瓷、石墨、高分子复合材料、塑料、泡沫等材料。
本发明是根据被泡沫芯材透射或反射后的太赫兹波强度得到样品的太赫兹图像。风力叶片泡沫芯材会对太赫兹波产生损耗，而粘接缺陷的边缘对太赫兹波的散射效应，会影响太赫兹波的强度分布。
通过将强度信号线性量化成RGB图像可得出物体内部的形状、缺陷及边缘位置。不同频率的太赫兹波被透镜汇聚到样件的不同深度，太赫兹波穿过泡沫芯材与粘接层后被金属基板反射，在混频器中与入射信号进行混频并被太赫兹探测器接收，其中含有被测量的相对距离信息，经过锁相放大器的信号放大传输至上位机进行图像处理。装置一次二维扫描，可同时采集风力叶片泡沫芯材不同深度的太赫兹反射回波，形成层析成像，实现太赫兹无损检测。
连续太赫兹波检测系统：
信号发射接收所用时间为τ，太赫兹源发射的初始频率f0，带宽B，周期T，调频斜率μ＝B/T，设太赫兹发射头到目标的距离为R，ct为太赫兹波在空气中的传播速度，则τ＝2R/ct。
图2中，分别用实线和虚线表示太赫兹源发射频率ft和经金属基板反射后太赫兹波频率fr，f0为ft的初始值，f1则为fr的终值，t为检测时间。
发射信号频率为:ft＝f0+μt 式(1)
反射信号频率为:fr＝ft(t－τ)＝f0+μ(t－τ) 式(2)
混频器的作用是使发射信号频率与反射信号频率相减得到中频信号频率fif，表示为:fif＝ft－fr＝μτ＝2BR/(TcT) 式(3)
于是有:R＝TcTfif/(2B) 式(4)
频带宽B和周期T在特定情况下，混频器中频信号输出fif与到材料的距离R成正比，根据中频信号fif可得出太赫兹透镜到材料的距离，中频信号的振幅反映了被测材料的折射率，通过混频器发出的差拍信号，即可实现目标材料的振幅与相位成像。锁相放大器将信号放大后交送上位机，再经由图像处理后，可显示Z(厚度方向)方向0～1000mm纵深范围内任一深度样件的X-Y二维截面太赫兹图像。
太赫兹发射源，频率范围为0.3THZ～1.5THZ，能够通过频率源、信号源的整合发射稳定频率范围的太赫兹信号，太赫兹波功率300mW～5W。
喇叭天线，材质采用金属金，发射角度范围-60°～60°
混频器，使发射信号频率与反射信号频率相减得到中频信号，中频信号的振幅反映了被测材料的折射率，通过混频器发出的差拍信号，即可实现目标材料的振幅与相位成像。
太赫兹探测器，探测频率范围0.1-2.0THZ。
锁相放大器，也称为相位检测器，是一种可以从干扰极大的环境(信噪比可低至-60dB，甚至更低)中分离出特定载波频率信号的放大器。
太赫兹透镜，材质选用高纯硅，直径40cm，焦距30cm。
实际应用中，对风力叶片泡沫芯材造成破坏性影响的主要是较大的缺陷。得到泡沫芯材太赫兹图像后，将所得原始RGB图像转换为相应的灰度图像，并进行中值滤波，从而消除孤立的噪声点。根据连续太赫兹图像中缺陷的灰度值分布特点，选择图像分割算法对其进行分割，从而提取缺陷前景，实现对缺陷面积定量。所述图像处理方法，分为缺陷图像分割和缺陷面积检测两步。
1)缺陷图像分割
所述算法通过迭代运算实现目标函数(类内误差平方和函数)的不断优化，在图像分割的应用就是根据图像中像素和c个聚类中心的每一个中心间的加权隶属度，对目标函数进行多次迭代优化，得到聚类目标函数数J最小时的聚类中心矩阵P与模糊划分矩阵U。
uij表示第j个样本属于第i个聚类的隶属度，U由uij构成。在模糊聚类算法中，
其聚类目标函数J为：
式中:
n为图像中总像素个数；
m为加权指数，决定分类模糊程度；
dij＝||pi-xj||是样本xj与聚类中心pi两者间的欧式距离，所有pi构成P。
为约束
引入拉格朗日乘子λ，则聚类目标函数可以写成：
为得到聚类目标函数J最小时的模糊划分矩阵U，对目标函数J求关于λ和uij的偏导数，令其为0，则有:
由(7)(8)得：
令目标函数关于聚类中心pi的偏导数为0，则有:
最终得到pi的更新公式：
算法的主要实现步骤为：
(1)确定聚类中心个数c和加权指数m,设定迭代阈值ε＞0，初始迭代次数为0，并初始化聚类中心矩阵P；
(2)计算更新U(N),P(N)；
(3)若||P(N+1)-P(N)||＜ε停止迭代，否则将置迭代系数N＝N+1，并返回步骤(2)继续迭代。
针对风力叶片泡沫芯材缺陷的连续太赫兹图像，目标缺陷即为一类。
算法的主要思想:评价聚类性能的准则函数在多次迭代过程后达到最优，并将数据集划分为不同的聚类，从而使生成的每个聚类达到内紧凑与类间独立。粘接层中存在的缺陷使反射回波的强度存在差异，在连续太赫兹波图像上，粘接缺陷体现为一种单独的聚类，像素聚类与图像分割一致，从而将缺陷图像与其它类分隔开。
2)缺陷面积估算
将粘接层太赫兹图像进行图像分割得到含有缺陷信息的二值图像，遍历图像中所有缺陷像素点，即可求出缺陷类像素点的总和。由于连续太赫兹系统的的X轴步进精度和Y轴的步进精度均为1mm，则每个像素点的所代表的真实面积为1mm2。所以缺陷类像素点总和的数值即为缺陷面积的数值。
本发明能够通过确定风力发电泡沫芯材中的缺陷位置进行泡沫芯材缺陷的检测，提高了风力叶片泡沫芯材的检测效率；同时也提升了风力发电叶片整体的安全性能。通过将缺陷种类、位置信息反馈给缺陷处理单元实现了缺陷泡沫板材的重复利用，减少了资源的浪费。
Claims (10)
1.一种太赫兹主动式风力叶片泡沫芯材无损检测方法，其特征在于包括
该方法是利用连续太赫兹波检测系统的太赫兹发射源对样品风力叶片泡沫芯材发射太赫兹电磁波，根据被泡沫芯材透射和/或反射后的太赫兹波强度得到样品的太赫兹图像；
连续太赫兹波检测系统对风力叶片芯材进行扫描式成像；连续太赫兹波检测系统使太赫兹波聚焦在风力叶片泡沫芯材上，尤其粘接层处，风力叶片泡沫芯材会对太赫兹波产生损耗，而粘接缺陷的边缘对太赫兹波的散射效应，会影响太赫兹波的强度分布，获得层析图像，尤其粘接层处的层析图像，根据焦平面附近的图像，可直观判断出粘接面有无缺陷；
连续太赫兹波检测系统的一次二维扫描可以同时采集风力叶片泡沫芯材不同深度的太赫兹反射回波，形成层析成像，实现太赫兹无损检测。
2.根据权利要求1所述的一种太赫兹主动式风力叶片泡沫芯材无损检测方法，其特征在于：
得到泡沫芯材太赫兹图像后，将所得原始RGB图像转换为相应的灰度图像，并进行中值滤波，消除孤立的噪声点；根据连续太赫兹图像中缺陷的灰度值分布特点，选择图像分割算法对其进行分割，从而提取缺陷前景，实现对缺陷面积定量。
3.根据权利要求1所述的一种太赫兹主动式风力叶片泡沫芯材无损检测方法，其特征在于：
通过将强度信号线性量化成RGB图像可以得出样品物体内部的形状、缺陷及边缘位置；不同频率的太赫兹波被透镜汇聚到样件的不同深度，太赫兹波穿过泡沫芯材与粘接层后被金属基板反射，在混频器中与入射信号进行混频并被太赫兹探测器接收，其中含有被测量的相对距离信息，经过锁相放大器的信号放大传输至上位机进行图像处理。
4.根据权利要求3所述的一种太赫兹主动式风力叶片泡沫芯材无损检测方法，其特征在于：
连续太赫兹波检测系统的工作方式和运算是：
信号发射接收所用时间为τ，
太赫兹源发射的初始频率f0，
带宽B，
周期T，
调频斜率μ＝B/T，
设太赫兹发射头到目标的距离为R，
ct为太赫兹波在空气中的传播速度，则τ＝2R/ct，
分别用实线和虚线表示太赫兹源发射频率ft和经金属基板反射后太赫兹波频率fr，
f0为ft的初始值，
f1则为fr的终值，
t为检测时间，
发射信号频率为:ft＝f0+μt 式(1)
反射信号频率为:fr＝ft(t－τ)＝f0+μ(t－τ) 式(2)
混频器的作用是使发射信号频率与反射信号频率相减得到中频信号频率fif，表示为:fif＝ft－fr＝μτ＝2BR/(TcT) 式(3)
于是有:R＝TcTfif/(2B) 式(4)
频带宽B和周期T在特定情况下，混频器中频信号输出fif与到材料的距离R成正比，根据中频信号fif可得出太赫兹透镜到材料的距离，中频信号的振幅反映了被测材料的折射率，通过混频器发出的差拍信号，即可实现目标材料的振幅与相位成像；锁相放大器将信号放大后交送上位机，再经由图像处理后，可显示Z厚度方向0～1000mm纵深范围内任一深度样件的方向X-Y二维截面太赫兹图像。
5.根据权利要求2、3或4所述的一种太赫兹主动式风力叶片泡沫芯材无损检测方法，其特征在于：
对于所获得的太赫兹图像，结合中值滤波，采用模糊c均值FCM聚类算法对缺陷图像进行处理，分割得到含有目标缺陷信息的二值图像，定位目标缺陷信息所示的缺陷位置，并计算像素点得到缺陷面积。
6.根据权利要求5所述的一种太赫兹主动式风力叶片泡沫芯材无损检测方法，其特征在于：
对于缺陷图像处理方法，分为[1]缺陷图像分割和[2]缺陷面积估算检测两步；
[1]缺陷图像分割
缺陷图像分割算法通过迭代运算实现目标函数：类内误差平方和函数的不断优化，在图像分割的应用就是根据图像中像素和c个聚类中心的每一个中心间的加权隶属度，对目标函数进行多次迭代优化，得到聚类目标函数数J最小时的聚类中心矩阵P与模糊划分矩阵U；
uij表示第j个样本属于第i个聚类的隶属度，U由uij构成；在模糊聚类算法中，
其聚类目标函数J为：
式(5)中:
n为图像中总像素个数；
m为加权指数，决定分类模糊程度；
dij＝||pi-xj||是样本xj与聚类中心pi两者间的欧式距离，所有pi构成P；
为约束
引入拉格朗日乘子λ，则聚类目标函数可以写成：
为得到聚类目标函数J最小时的模糊划分矩阵U，对目标函数J求关于λ和uij的偏导数，令其为0，则有:
由式(7)、式(8)得：
令目标函数关于聚类中心pi的偏导数为0，则有:
最终得到pi的更新公式：
该算法的主要实现步骤为：
(1)确定聚类中心个数c和加权指数m，设定迭代阈值ε＞0，初始迭代次数为0，并初始化聚类中心矩阵P；
(2)计算更新U(N),P(N)；
(3)若||P(N+1)-P(N)||＜ε停止迭代，否则将置迭代系数N＝N+1，并返回步骤(2)继续迭代；
针对风力叶片泡沫芯材缺陷的连续太赫兹图像，目标缺陷即为一类；该算法的主要思想：评价聚类性能的准则函数在多次迭代过程后达到最优，并将数据集划分为不同的聚类，从而使生成的每个聚类达到内紧凑与类间独立；粘接层中存在的缺陷使反射回波的强度存在差异，在连续太赫兹波图像上，粘接缺陷体现为一种单独的聚类，像素聚类与图像分割一致，从而将缺陷图像与其它类分隔开；
[2]缺陷面积估算
将粘接层太赫兹图像进行图像分割得到含有缺陷信息的二值图像，遍历图像中所有缺陷像素点，即可求出缺陷类像素点的总和；
由于连续太赫兹系统的的X轴步进精度和Y轴的步进精度均为1mm，则每个像素点的所代表的真实面积为1mm2；所以缺陷类像素点总和的数值即为缺陷面积的数值。
7.根据权利要求1～4任一所述的一种太赫兹主动式风力叶片泡沫芯材无损检测方法，其特征在于：
连续太赫兹波检测系统的太赫兹发射源发射的太赫兹波采用频率范围为0.3～1.5THz的反射式连续太赫兹波；能够通过频率源、信号源的整合发射稳定频率范围的太赫兹信号，太赫兹波功率300mW～5W；
连续太赫兹波检测系统的太赫兹探测器，其探测频率范围0.1～2.0THZ。
8.根据权利要求1～4任一所述的一种太赫兹主动式风力叶片泡沫芯材无损检测方法，其特征在于：
风力叶片芯材采用玻璃钢、陶瓷、石墨、高分子复合材料、塑料或泡沫的非金属芯材。
9.一种实现如权利要求1～4任一所述的太赫兹主动式风力叶片泡沫芯材无损检测方法的装置，其特征在于该装置的结构是：
在太赫兹发射源的发射轴向上配置有喇叭天线、混频器、太赫兹透镜、样品风力叶片泡沫芯材、反射板，
于太赫兹发射轴向上在太赫兹透镜的上游配置有混频器，
在混频器与发射轴交点处的垂直方向配置有太赫兹探测器，
太赫兹探测器连接锁相放大器，锁相放大器连接上位机；
以上构成连续太赫兹波检测系统。
10.根据权利要求9所述的一种实现太赫兹主动式风力叶片泡沫芯材无损检测方法的装置，其特征在于：
连续太赫兹波检测系统的太赫兹发射源发射的太赫兹波采用频率范围为0.3～1.5THz的反射式连续太赫兹波；能够通过频率源、信号源的整合发射稳定频率范围的太赫兹信号，太赫兹波功率300mW～5W；
连续太赫兹波检测系统的太赫兹探测器，其探测频率范围0.1～2.0THZ；
喇叭天线材质采用金属金，发射角度范围-60°～60°；
混频器：使发射信号频率与反射信号频率相减得到中频信号，中频信号的振幅反映了被测材料的折射率，通过混频器发出的差拍信号，即可实现目标材料的振幅与相位成像；
锁相放大器是从干扰环境中分离出特定载波频率信号的放大器；
太赫兹透镜采用高纯硅太赫兹透镜，直径40cm，焦距30cm。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110644852.4A CN113406039A (zh) | 2021-06-09 | 2021-06-09 | 一种太赫兹主动式风力叶片泡沫芯材无损检测方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110644852.4A CN113406039A (zh) | 2021-06-09 | 2021-06-09 | 一种太赫兹主动式风力叶片泡沫芯材无损检测方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113406039A true CN113406039A (zh) | 2021-09-17 |
Family
ID=77683336
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110644852.4A Pending CN113406039A (zh) | 2021-06-09 | 2021-06-09 | 一种太赫兹主动式风力叶片泡沫芯材无损检测方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113406039A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117969449A (zh) * | 2024-03-29 | 2024-05-03 | 三峡金沙江云川水电开发有限公司 | 一种太赫兹检测发电机定子线棒绝缘缺陷的方法及系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106525862A (zh) * | 2016-10-27 | 2017-03-22 | 首都师范大学 | 利用太赫兹成像检测层状绝缘材料内部缺陷的方法和装置 |
| CN107677608A (zh) * | 2016-08-02 | 2018-02-09 | 北京遥感设备研究所 | 一种太赫兹准光主动扫描反射成像系统 |
| WO2019127949A1 (zh) * | 2017-12-28 | 2019-07-04 | 深圳市太赫兹科技创新研究院 | 多频太赫兹检测装置、系统及其方法 |
-
2021
- 2021-06-09 CN CN202110644852.4A patent/CN113406039A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107677608A (zh) * | 2016-08-02 | 2018-02-09 | 北京遥感设备研究所 | 一种太赫兹准光主动扫描反射成像系统 |
| CN106525862A (zh) * | 2016-10-27 | 2017-03-22 | 首都师范大学 | 利用太赫兹成像检测层状绝缘材料内部缺陷的方法和装置 |
| WO2019127949A1 (zh) * | 2017-12-28 | 2019-07-04 | 深圳市太赫兹科技创新研究院 | 多频太赫兹检测装置、系统及其方法 |
Non-Patent Citations (1)
| Title |
|---|
| 郑晓康: "基于连续太赫兹波的航天复合材料检测方法研究", 《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117969449A (zh) * | 2024-03-29 | 2024-05-03 | 三峡金沙江云川水电开发有限公司 | 一种太赫兹检测发电机定子线棒绝缘缺陷的方法及系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Doctor et al. | SAFT—the evolution of a signal processing technology for ultrasonic testing | |
| US7236564B2 (en) | Linear array detector system and inspection method | |
| JP6018047B2 (ja) | 複合材料を検査するためのｒｆ反射 | |
| Zimdars et al. | Technology and applications of terahertz imaging non‐destructive examination: inspection of space shuttle sprayed on foam insulation | |
| Abou‐Khousa et al. | Comparison of X‐ray, millimeter wave, shearography and through‐transmission ultrasonic methods for inspection of honeycomb composites | |
| WO1988008530A1 (en) | Gamma ray flaw detection system | |
| JP4369915B2 (ja) | 電磁波イメージングシステム | |
| CN113406039A (zh) | 一种太赫兹主动式风力叶片泡沫芯材无损检测方法及装置 | |
| CN108872277A (zh) | 基于龙虾眼透镜的x射线无损探伤装置 | |
| Friederich et al. | Continuous wave terahertz inspection of glass fiber reinforced plastics with semi-automatic 3-D image processing for enhanced defect detection | |
| CN113686960B (zh) | 声场阈值分割的相控阵曲面全聚焦成像优化方法及系统 | |
| CN104422715A (zh) | 基于动态子窗口技术的激光扫描热波成像方法及装置 | |
| US5247302A (en) | Entropy-based signal receiver | |
| Ravuri et al. | Microwave and millimeter wave near‐field methods for evaluation of radome composites | |
| CN111721672B (zh) | 基于多普勒和声呐图像技术的多相测试方法 | |
| Fang et al. | Visualization and quantitative evaluation of delamination defects in GFRPs via sparse millimeter-wave imaging and image processing | |
| CN108459088A (zh) | 一种相控阵超声波信号数据的三维转换和融合方法 | |
| CN114449431B (zh) | 一种使用太赫兹波对扬声器振膜进行无损检测的方法及系统 | |
| Liu et al. | Manifold learning and segmentation for ultrasonic inspection of defects in polymer composites | |
| Matheis et al. | Terahertz FMCW inspection of GFRP composites: Comparison with conventional NDT techniques and enhanced defect detection capability through | |
| CN114487115B (zh) | 一种基于Canny算子与超声平面波成像相结合的高分辨缺陷无损检测方法 | |
| Kim et al. | Quantitative characterization of multiple delaminations in laminated composites using the compton backscatter technique | |
| CN117554320A (zh) | 一种基于太赫兹调频连续波合成孔径成像的无损检测系统及其使用方法 | |
| Martin et al. | Signal processing of leaky Lamb wave data for defect imaging in composite laminates | |
| US5280291A (en) | Thermodynamics-based signal receiver |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210917 |