CN103969337B - 一种基于矢量全聚焦成像的超声阵列裂纹类缺陷方向识别方法 - Google Patents
一种基于矢量全聚焦成像的超声阵列裂纹类缺陷方向识别方法 Download PDFInfo
- Publication number
- CN103969337B CN103969337B CN201410191263.5A CN201410191263A CN103969337B CN 103969337 B CN103969337 B CN 103969337B CN 201410191263 A CN201410191263 A CN 201410191263A CN 103969337 B CN103969337 B CN 103969337B
- Authority
- CN
- China
- Prior art keywords
- rightarrow
- vector
- array
- imaging point
- imaging
- 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.)
- Withdrawn - After Issue
Links
- 239000013598 vector Substances 0.000 title claims abstract description 186
- 238000003384 imaging method Methods 0.000 title claims abstract description 152
- 230000007547 defect Effects 0.000 title claims abstract description 77
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000003491 array Methods 0.000 claims abstract description 29
- 238000010586 diagram Methods 0.000 claims abstract description 27
- 230000005284 excitation Effects 0.000 claims description 33
- 238000001514 detection method Methods 0.000 claims description 31
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000002474 experimental method Methods 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000005242 forging Methods 0.000 claims description 7
- 230000003247 decreasing effect Effects 0.000 claims description 6
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 229910000831 Steel Inorganic materials 0.000 claims description 2
- 239000002131 composite material Substances 0.000 claims description 2
- 230000008878 coupling Effects 0.000 claims description 2
- 238000010168 coupling process Methods 0.000 claims description 2
- 238000005859 coupling reaction Methods 0.000 claims description 2
- 239000010959 steel Substances 0.000 claims description 2
- 230000002194 synthesizing effect Effects 0.000 claims description 2
- 230000001902 propagating effect Effects 0.000 claims 1
- 238000009659 non-destructive testing Methods 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 235000019687 Lamb Nutrition 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明涉及一种基于矢量全聚焦成像的超声阵列裂纹类缺陷方向识别方法,属于无损检测领域。该方法通过对换能器阵列接收到的信号进行矢量全聚焦处理得到全局矢量图,由缺陷处反射信号能量最强来确定缺陷的位置,然后根据缺陷的位置提取缺陷的局部矢量图,最后利用局部矢量图中矢量的方向来确定缺陷的方向。该矢量全聚焦成像的原理是在换能器阵列中构造多个子阵列,通过子阵列计算任意成像点处的单位方向矢量,利用该单位方向矢量对换能器阵列在任意成像点处补偿后的回波幅值进行矢量化,得到任意成像点处的幅值矢量。
Description
技术领域
本发明涉及一种基于矢量全聚焦成像的超声阵列裂纹类缺陷方向识别方法,该方法涉及裂纹类缺陷的方向识别技术,属于无损检测技术领域。
背景技术
锻件作为组成各种机械设备的基础部件,在航天、核电、船舶等关键领域中应用广泛。由于锻件形状和生产工艺复杂,在生产过程中容易产生裂纹类缺陷。锻件中存在缺陷将显著地降低其力学性能,甚至在使用过程中还会发生断裂,造成严重的安全与质量事故。对于裂纹类缺陷都存在一定的方向性,获取缺陷的方向性能够有效地改进锻件生产工艺过程,提高生产效率。
超声波探伤作为五大常规检测手段之一,对锻件内部裂纹等危害性缺陷检测灵敏度较高,该方法是锻件缺陷检测的主要发展方向。但传统超声检测存在检测效率低、分辨率低、检测结果直观性差等缺点,并且无法对锻件内部裂纹等危害性缺陷的方向、形状等特征进行识别。为克服上述缺点,可以利用多个阵元组成的阵列换能器超声相控阵***进行无损检测,实现对缺陷的方向等特征识别。
超声相控阵检测技术是通过电子***控制阵列换能器中各个阵元,按照一定的延迟时间,规则地激励和接收超声波,控制超声波束的偏转和聚焦,来实现对锻件内部缺陷的无损检测。因此,超声相控阵技术具有快速、准确、适应性强等优点,而且对缺陷的检出率也高。同时,通过对采集到的数据进行后处理,可得到高精度和大范围的成像图,根据成像图中缺陷处局部图像信息的特点,可进一步提取缺陷的方向、形状等特征,从而实现对缺陷的特征识别。
专利《基于超声Lamb波的储罐底板腐蚀检测***及方法》(专利号:201110281076)中所述的全聚焦成像方法对时域信号数据进行处理,实现了对储罐底板结构中缺陷的定位。2013硕士论文“脉冲压缩技术及其在锻件相控阵检测中的应用研究”中对全聚焦成像及其优化方法进行了研究,实现对缺陷的定位检测,并通过优化方法提高了成像的精度。但全聚焦成像方法无法对裂纹类缺陷的方向进行识别。本专利的创新点在于提出一种矢量全聚焦成像方法,通过在阵列中构造多个子阵列,求解出任意成像点处的矢量,然后根据成像点矢量的方向来判断裂纹类缺陷的方向。此方法不仅实现了对裂纹类缺陷的定位检测,而且能有效地识别出裂纹类缺陷的方向。此外,对缺陷方向、形状等特征识别也存在其它的方法,2012年发表的博士论文“基于散射系数矩阵法的超声兰姆波与典型缺陷交互作用研究”中利用散射系数矩阵成功地对圆孔和裂纹类缺陷进行识别,并能较好地识别出缺陷的方向和尺寸。
发明内容
本发明的目的在于提出一种基于矢量全聚焦成像的超声阵列裂纹类缺陷方向识别方法。该方法用于对裂纹类缺陷的方向识别,首先利用超声相控阵***采集得到全矩阵数据,然后对采集到的全矩阵数据进行矢量全聚焦处理得到被测试件的全局矢量图,由目标缺陷处反射信号能量最强来确定缺陷的位置,根据缺陷的位置提取目标缺陷的局部矢量图,最后利用局部矢量图中矢量的方向来确定目标缺陷的方向。该矢量全聚焦成像的原理是在阵列换能器中构造多个子阵列,通过子阵列计算任意成像点处的单位方向矢量,利用该单位方向矢量对全阵列在任意成像点处补偿后的回波幅值进行矢量化,得到全阵列在任意成像点处的幅值矢量。
该方法需要的检测装置包括超声波激励/接收模块、采集模块、计算机、线性阵列换能器和被测试件,其中,超声波激励/接收模块一端与计算机、采集模块连接,另一端与线性阵列换能器连接,采集模块另一端与计算机相连,线性阵列换能器与被测试件通过耦合介质进行耦合,如图1所示。所述的超声波激励/接收模块采用的是Multi2000系列相控阵仪器,所述的采集模块采用的是与Multi2000系列相控阵仪器相配套的采集软件,所述的线性阵列换能器是由多个阵元组成的线性换能器。在进行检测实验时,首先在计算机的控制下超声波激励/接收模块产生激励信号,通过阵列换能器激励出超声波信号沿被测试件发射出去,并通过阵列换能器接收反射回来的超声波信号;然后又经过超声波激励/接收模块传输给采集模块;最后通过计算机控制采集模块里的采集软件进行数据采集,即可获得检测的时域信号。
为实现上述目的,本发明采用的技术方案为一种基于矢量全聚焦成像的超声阵列裂纹类缺陷方向识别方法,具体可以按照以下步骤实施检测,方法流程如图2所示。
步骤一:在如图1所示的检测装置下进行实验,其中,实验所用阵列换能器中心频率为f,阵元总个数为N,单个阵元的宽度为a,相邻两阵元的中心距离为p,超声波在被测试件中的传播波速为c,则波长通过全矩阵模式采集得到时域信号f(i)j(t)(i=1,2,3,…,N,j=1,2,3,…,N),其中,下标(i)表示阵列换能器中第i个阵元激励,j表示阵列换能器中第j个阵元接收。
步骤二:对采集到的时域信号f(i)j进行希尔伯特变换得到信号的包络线,将此时得到的信号称为包络信号g(i)j(t)。
步骤三:建立成像坐标系,如图3所示。定义表示原点到成像点P的向量,表示第i个激励阵元的位置向量,表示第j个接收阵元的位置向量,则第i个激励阵元到成像点P的向量第j个接收阵元到成像点P的向量其中,图中x轴表示与换能器位置平行的方向,z轴表示与换能器位置垂直的方向,以向量为例,表示向量的模,表示向量在x轴方向上的分量,表示向量在z轴方向上的分量,后续步骤公式中的其它向量也是类似表示。
步骤四:将全部阵列数据按激励、接收阵元与成像点的距离计算声波传播时间,在每个成像点位置进行聚焦,并将信号幅值叠加。因此,可通过公式(1)计算全阵列在每个成像点的幅值
式中,c为声波在被测试件中的传播波速;
步骤五:将步骤四中的幅值利用波束补偿因子和衰减补偿因子进行补偿,得到补偿后的幅值可通过公式(2)计算得出。
式中,P为波束补偿因子,B衰减补偿因子
其中,激励、接收阵元在任意成像点的波束补偿因子可分别表示为:
式中,a为单个阵元的宽度,λ为波长,sinc函数为数学中的辛格函数,当用sinc(x)表示时,其非归一化的sinc函数表示为
激励、接收阵元在任意成像点的衰减补偿因子可表示为:
步骤六:构造子阵列。将阵元总数为N的阵列换能器中多个连续的阵元作为一个子阵列,对应的时域信号称为子阵列数据。将该阵列换能器划分为K个子阵列,每个子阵列中含有N1个阵元(N1<N),相邻两个子阵列间的阵元个数为M(M<N)。则第k个子阵列对应的阵元在全阵列中的序号最小值为1+M(k-1),最大值为1+M(k-1)+N1,其中,k=1,2,3,…,K。
步骤七:计算每个子阵列在每个成像点的幅值矢量。
根据步骤六划分好的子阵列,计算每个子阵列在每个成像点的幅值矢量,可分为以下几步:
(1)根据每个子阵列中对应的激励、接收阵元,重复步骤四、五,公式中下标i、j的最小值为1+M(k-1),最大值为1+M(k-1)+N1,可得到每个子阵列在每个成像点处矢量的幅值其中,上标k表示第k个子阵列。
(2)根据每个子阵列对应的激励、接收阵元到成像点的位置向量可得到每个子阵列在每个成像点的单位方向矢量具体可通过公式(6)计算得出。
其中,是第i个阵元激励超声波入射到任意成像点后经第j个阵元接收形成的法线方向的单位方向矢量,依据反射定理可知,该单位方向矢量的方向会与反射面垂直,具体计算公式可表示为:
(3)将上两步中得到的和对应相乘即可得到每个子阵列在每个成像点处的幅值矢量可通过公式(8)计算得出。
步骤八:对步骤七中K个子阵列的幅值矢量进行合成,得到K个子阵列在每个成像点的合成幅值矢量可通过公式(9)计算得出。
式中,α取1到正无穷的任意数,当α取值越大,任意成像点的合成幅值矢量的方向越接近成像点处反射信号能量最强的方向。
步骤九:将合成幅值矢量进行单位化,然后乘以步骤四中的可得到在全阵列N下每个成像点的幅值矢量可通过公式(10)计算得出。
根据上述计算可知,全阵列在每个成像点的幅值矢量的大小与步骤四中的全阵列在每个成像点的幅值相等,并且矢量的方向会与成像点处的反射面垂直。
步骤十:将步骤九中的全阵列在每个成像点的幅值矢量进行成像显示,可得到全阵列在每个成像点的全局矢量成像图。
步骤十一:根据步骤十中的矢量图,确定目标缺陷的位置,提取目标缺陷的局部矢量图然后将步骤四中的幅值进行分贝化,找出局部矢量图中幅值的最大值,求解最大值下降-6dB所对应的成像区域面积。最后根据局部矢量图中矢量的方向与目标缺陷垂直,依据几何关系,目标缺陷的方向可通过公式(11)计算得出。
式中,和表示与阵列换能器平行和垂直的单位向量,A表示在目标缺陷局部矢量图中,由幅值的最大值下降-6dB所对应的成像区域面积,arctan函数为数学中的反正切函数,θm表示矢量与z轴正向的夹角,即等于缺陷与x轴正向的夹角。
与现有检测方法相比,本发明具有以下优点:(1)对阵列数据进行幅值叠加后处理,实现波束在任意成像点处的聚焦,并利用波束补偿因子和衰减补偿因子,提高缺陷检测精度,对小缺陷的检出能力有明显的提高。(2)在阵列中构造多个子阵列,计算子阵列中各成像点的矢量,根据矢量的方向确定缺陷的方向,获取更多的缺陷特征信息。
附图说明
图1为本发明的检测装置框图
图2为本发明方法的流程图
图3为本发明方法实施例中建立的成像坐标系
图4为锻件的局部矢量全聚焦图
图5为锻件的局部矢量全聚焦图的放大图
具体实施方式
以下结合附图和实施例对本发明作进一步详细说明。
本实例检测仪器为Multi2000系列相控阵设备,其中包括超声波激励/接收模块、采集模块及相控阵线性阵列换能器。实验中使用的线性阵列换能器中心频率f=5MHz,阵元总个数N=32,单个阵元宽度a=0.5mm,相邻两阵元的中心距离p=0.6mm,声波在被测试件中的传播波速c=5900m/s,则波长检测对象为圆形钢制锻件,直径为300mm,厚度为40mm。
本发明提出的一种基于矢量全聚焦成像的超声阵列裂纹类缺陷方向识别方法的具体步骤为:
步骤一:在如图1所示的检测装置下进行实验,通过全矩阵模式采集得到时域信号f(i)j(t)(i=1,2,3,…,32,j=1,2,3,…,32),其中,下标(i)表示阵列换能器中第i个阵元激励,j表示阵列换能器中第j个阵元接收。
步骤二:对采集到的时域信号f(i)j进行希尔伯特变换得到信号的包络线,将此时得到的信号称为包络信号g(i)j(t)。
步骤三:建立成像坐标系。
以换能器的正中心为原点建立直角坐标系,对向量的计算将其分解为x、z方向来求解。定义两个一维数组X(s)、Z(s),数组大小均为32,s为阵列换能器中阵元的编号,则s=1,2,3,…,32,其中第s个元素(X(s),Z(s))表示第s个阵元的x和z轴坐标,其表达式为:X(s)=((s-1)-(32-1)/2)×0.6,Z(s)=0。则第i个激励阵元的位置向量在x、z方向上的分量可表示为: 第j个接收阵元的位置向量在x、z方向上的分量可表示为:
设置成像区域x方向长度为150mm,z方向长度为150mm,成像精度为0.2mm。定义两个一维数组Px(m)、Pz(n),数组大小分别为MX、MZ,其中MX=150/0.2=750,MZ=150/0.2=750,m、n分别为x、z方向离散点编号,则m=1,2,3,…,750,n=1,2,3,…,750。当m取特定值m0时,Px(m)表示x方向上编号为m0对应的点的x轴坐标;当n取特定值n0时,Pz(n)表示z方向上编号为n0对应的点的z轴坐标,则任意编号(m0,n0)对应的任意成像点的x、z轴坐标可表示为:Px(m0)=-0.75-0.0002×(m0-1),Pz(n0)=0.0002×(n0-1)。
则任意编号(m0,n0)对应的任意成像点P的位置向量在x、z方向上的分量为: 相应地,第i个激励阵元到成像点P的位置向量在x、z方向上的分量 可表示为:
第j个激励阵元到成像点P的位置向量在x、z方向上的分量 可表示为:
因此,第i个激励阵元、第j个接收阵元到任意成像点的距离可表示为:
后续步骤中的均可通过此步中的表达式计算得出。
步骤四:将全部阵列数据按激励、接收阵元与成像点的距离计算声波传播时间,在每个成像点位置进行聚焦,并将信号幅值叠加。因此,可以通过公式(1)计算全阵列在每个成像点的幅值
式中,c为声波在被测试件中传播的波速。
步骤五:将步骤四中的幅值利用波束补偿因子和衰减补偿因子进行补偿,得到补偿后的幅值可通过公式(2)计算得出。
式中P为波束补偿因子,B衰减补偿因子;
其中,激励、接收阵元在任意成像点的波束补偿因子可分别表示为:
式中a为单个阵元的宽度,λ为波长,sinc函数为数学中的辛格函数,当用sinc(x)表示时,其非归一化的sinc函数表示为
激励、接收阵元在任意成像点的衰减补偿因子可表示为:
步骤六:构造子阵列。将阵元总数为32的阵列换能器中多个连续的阵元作为一个子阵列,对应的时域信号作为子阵列数据。将该阵列换能器划分为7个子阵列,每个子阵列中含有8个阵元,相邻两个子阵列间的阵元个数为4个。则第k个子阵列对应的阵元在全阵列中的序号最小值为1+4(k-1)=4k-3,最大值为1+4(k-1)+7=4k+4,其中,k=1,2,3,…,K。
步骤七:计算每个子阵列在每个成像点的幅值矢量。
根据步骤六划分好的子阵列,计算每个子阵列在每个成像点的矢量,可分为以下几步:
(1)根据每个子阵列中对应的激励、接收阵元,重复步骤四、五,公式中下标i、j的最小值为4k-3,最大值为4k+4,可得到每个子阵列在每个成像点处矢量的幅值其中,上标k表示第k个子阵列。
(2)根据每个子阵列对应的激励、接收阵元到成像点的位置向量可得到每个子阵列在每个成像点的单位方向矢量量即为每个子阵列在每个成像点处矢量的方向。该单位方向矢量方向上的分量可通过公式(6)计算得出。
其中,是第i个阵元激励超声波入射到任意成像点后经第j个阵元接收形成的法线方向的单位方向矢量,依据反射定理可知,该单位方向矢量的方向会与反射面垂直,可通过公式(7)计算得出。该单位方向矢量在x、z方向上的分量 可表示:
式中,为矢量的模,具体表达式为:
(3)将上两步中得到的分向量和对应相乘即可得到每个子阵列在每个成像点的幅值矢量可通过公式(8)计算得出。该幅值矢量方向上的分向量可表示为:
步骤八:取α=4,将步骤七中7个分向量分别在x、z方向上进行合成,得到7个子阵列在每个成像点的合成幅值矢量可通过公式(9)计算得出。该合成幅值矢量方向上的分向量可表示为:
步骤九:将步骤八中得到的进行单位化,然后乘以步骤四中的得到在32阵元下每个成像点的幅值矢量可通过公式(10)计算得出,该32阵元下的幅值矢量方向上的分向量可表示为:
式中,为矢量的模,具体表达式为:
根据上述计算可知,32阵元在每个成像点的幅值矢量的大小与步骤四中的32阵元在每个成像点幅值相等,并且幅值矢量的方向会与成像点处的反射面垂直。
步骤十:将步骤九中的全阵列在每个成像点的幅值矢量进行成像显示,可得到全阵列在每个成像点的全局矢量成像图。
步骤十一:根据步骤十中的全局矢量图,确定目标缺陷的位置,提取目标缺陷的局部矢量图和局部矢量图的放大图,如图4、图5所示。然后将步骤四中的幅值进行分贝化,找出局部矢量图中幅值的最大值,求解最大值下降-6dB所对应的成像区域面积。最后根据局部矢量图中矢量的方向与目标缺陷垂直,依据几何关系,目标缺陷的方向可通过公式(11)计算得出。
式中,A表示在目标缺陷局部矢量图中,由幅值的最大值下降-6dB所对应的成像区域面积,arctan函数为数学中的反正切函数,θm表示矢量方向与z轴正向的夹角,即等于缺陷与x轴正向的夹角;
最后,通过局部矢量图可以看出,缺陷处矢量幅值较明显,用箭头表现的形式,而其它不含缺陷的地方由于幅值很小,看不出来箭头,只呈现出圆点的形式,由此可以判断出缺陷的位置;并且通过计算求得缺陷处的角度θm=3.21°,即表示缺陷与x轴正向的夹角为3.21°。
上述步骤只是本发明的一个典型实施例,本发明的实施不限于此。
Claims (3)
1.一种基于矢量全聚焦成像的超声阵列裂纹类缺陷方向识别方法,其特征在于:该方法按照以下步骤实施检测,
步骤一:在检测装置下进行实验,其中,实验所用阵列换能器中心频率为f,阵元总个数为N,单个阵元的宽度为a,相邻两阵元的中心距离为p,超声波在被测试件中的传播波速为c,则波长通过全矩阵模式采集得到时域信号f(i)j(t),i=1,2,3,…,N,j=1,2,3,…,N,其中,下标(i)表示阵列换能器中第i个阵元激励,j表示阵列换能器中第j个阵元接收;
步骤二:对采集到的时域信号f(i)j进行希尔伯特变换得到信号的包络线,将此时得到的信号称为包络信号g(i)j(t);
步骤三:建立成像坐标系;定义表示原点到成像点P的向量,表示第i个激励阵元的位置向量,表示第j个接收阵元的位置向量,则第i个激励阵元到成像点P的向量第j个接收阵元到成像点P的向量其中,x轴表示与换能器位置平行的方向,z轴表示与换能器位置垂直的方向,表示向量的模,表示向量在x轴方向上的分量,表示向量在z轴方向上的分量,后续步骤公式中的其它向量也是同样表示;
步骤四:将全部阵列数据按激励、接收阵元与成像点的距离计算声波传播时间,在每个成像点位置进行聚焦,并将信号幅值叠加;因此,可通过公式(1)计算全阵列在每个成像点的幅值
式中,c为声波在被测试件中的传播波速;
步骤五:将步骤四中的幅值利用波束补偿因子和衰减补偿因子进行补偿,得到补偿后的幅值可通过公式(2)计算得出;
式中,P为波束补偿因子,B衰减补偿因子;
其中,激励、接收阵元在任意成像点的波束补偿因子可分别表示为:
式中,a为单个阵元的宽度,λ为波长,sin c函数为数学中的辛格函数,当用sin c(x)表示时,其非归一化的sin c函数表示为
激励、接收阵元在任意成像点的衰减补偿因子可表示为:
步骤六:构造子阵列;将阵元总数为N的阵列换能器中多个连续的阵元作为一个子阵列,对应的时域信号称为子阵列数据;将该阵列换能器划分为K个子阵列,每个子阵列中含有N1个阵元(N1<N),相邻两个子阵列间的阵元个数为M,M<N;则第k个子阵列对应的阵元在全阵列中的序号最小值为1+M(k-1),最大值为1+M(k-1)+N1,其中,k=1,2,3,…,K;
步骤七:计算每个子阵列在每个成像点的幅值矢量;
根据步骤六划分好的子阵列,计算每个子阵列在每个成像点的幅值矢量,可分为以下几步:
(1)根据每个子阵列中对应的激励、接收阵元,重复步骤四、五,公式中下标i、j的最小值为1+M(k-1),最大值为1+M(k-1)+N1,可得到每个子阵列在每个成像点处矢量的幅值k=1,2,3,…,K,其中,上标k表示第k个子阵列;
(2)根据每个子阵列对应的激励、接收阵元到成像点的位置向量可得到每个子阵列在每个成像点的单位方向矢量具体可通过公式(6)计算得出;
其中,是第i个阵元激励超声波入射到任意成像点后经第j个阵元接收形成的法线方向的单位方向矢量,依据反射定理可知,该单位方向矢量的方向会与反射面垂直,具体计算公式可表示为:
(3)将上两步中得到的和对应相乘即可得到每个子阵列在每个成像点处的幅值矢量可通过公式(8)计算得出;
步骤八:对步骤七中K个子阵列的幅值矢量进行合成,得到K个子阵列在每个成像点的合成幅值矢量可通过公式(9)计算得出;
式中,α取1到正无穷的任意数,当α取值越大,任意成像点的合成幅值矢量的方向越接近成像点处反射信号能量最强的方向;
步骤九:将合成幅值矢量进行单位化,然后乘以步骤四中的可得到在全阵列N下每个成像点的幅值矢量可通过公式(10)计算得出;
根据上述计算可知,全阵列在每个成像点的幅值矢量的大小与步骤四中的全阵列在每个成像点的幅值相等,并且矢量的方向会与成像点处的反射面垂直;
步骤十:将步骤九中的全阵列在每个成像点的幅值矢量进行成像显示,可得到全阵列在每个成像点的全局矢量成像图;
步骤十一:根据步骤十中的全局矢量成像图,确定目标缺陷的位置,提取目标缺陷的局部矢量图然后将步骤四中的幅值进行分贝化,找出局部矢量图中幅值的最大值,求解最大值下降-6dB所对应的成像区域面积;最后根据局部矢量图中矢量的方向与目标缺陷垂直,依据几何关系,目标缺陷的方向可通过公式(11)计算得出;
式中,和表示与阵列换能器平行和垂直的单位向量,A表示在目标缺陷局部矢量图中,由幅值的最大值下降-6dB所对应的成像区域面积,arctan函数为数学中的反正切函数,θm表示矢量与z轴正向的夹角,即等于缺陷与x轴正向的夹角。
2.根据权利要求1所述的一种基于矢量全聚焦成像的超声阵列裂纹类缺陷方向识别方法,其特征在于:该方法需要的检测装置包括超声波激励/接收模块、采集模块、计算机、线性阵列换能器和被测试件,其中,超声波激励/接收模块一端与计算机、采集模块连接,另一端与线性阵列换能器连接,采集模块另一端与计算机相连,线性阵列换能器与被测试件通过耦合介质进行耦合;所述的超声波激励/接收模块采用的是Multi2000系列相控阵仪器,所述的采集模块采用的是与Multi2000系列相控阵仪器相配套的采集软件,所述的线性阵列换能器是由多个阵元组成的线性换能器;在进行检测实验时,首先在计算机的控制下超声波激励/接收模块产生激励信号,通过阵列换能器激励出超声波信号沿被测试件发射出去,并通过阵列换能器接收反射回来的超声波信号;然后又经过超声波激励/接收模块传输给采集模块;最后通过计算机控制采集模块里的采集软件进行数据采集,即可获得检测的时域信号。
3.根据权利要求1所述的一种基于矢量全聚焦成像的超声阵列裂纹类缺陷方向识别方法,其特征在于:检测仪器为Multi2000系列相控阵设备,其中包括超声波激励/接收模块、采集模块及相控阵线性阵列换能器;实验中使用的相控阵线性阵列换能器中心频率f=5MHz,阵元总个数N=32,单个阵元宽度a=0.5mm,相邻两阵元的中心距离p=0.6mm,超声波在被测试件中的传播波速c=5900m/s,则波长检测对象为圆形钢制锻件,直径为300mm,厚度为40mm;
一种基于矢量全聚焦成像的超声阵列裂纹类缺陷方向识别方法的具体步骤为:
步骤一:在检测装置下进行实验,通过全矩阵模式采集得到时域信号f(i)j(t),i=1,2,3,…,32,j=1,2,3,…,32,其中,下标(i)表示阵列换能器中第i个阵元激励,j表示阵列换能器中第j个阵元接收;
步骤二:对采集到的时域信号f(i)j进行希尔伯特变换得到信号的包络线,将此时得到的信号称为包络信号g(i)j(t);
步骤三:建立成像坐标系;
以换能器的正中心为原点建立直角坐标系,对向量的计算将其分解为x、z方向来求解;定义两个一维数组X(s)、Z(s),数组大小均为32,s为阵列换能器中阵元的编号,则s=1,2,3,…,32,其中第s个元素(X(s),Z(s))表示第s个阵元的x和z轴坐标,其表达式为:X(s)=((s-1)-(32-1)/2)×0.6,Z(s)=0;则第i个激励阵元的位置向量在x、z方向上的分量可表示为:第j个接收阵元的位置向量在x、z方向上的分量可表示为:
设置成像区域x方向长度为150mm,z方向长度为150mm,成像精度为0.2mm;定义两个一维数组Px(m)、Pz(n),数组大小分别为MX、MZ,其中MX=150/0.2=750,MZ=150/0.2=750,m、n分别为x、z方向离散点编号,则m=1,2,3,…,750,n=1,2,3,…,750;当m取特定值m0时,Px(m)表示x方向上编号为m0对应的点的x轴坐标;当n取特定值n0时,Pz(n)表示z方向上编号为n0对应的点的z轴坐标,则任意编号(m0,n0)对应的任意成像点的x、z轴坐标可表示为:Px(m0)=-0.75-0.0002×(m0-1),Pz(n0)=0.0002×(n0-1);
则任意编号(m0,n0)对应的任意成像点P的位置向量在x、z方向上的分量为: 相应地,第i个激励阵元到成像点P的位置向量在x、z方向上的分量可表示为:
第j个激励阵元到成像点P的位置向量在x、z方向上的分量可表示为:
因此,第i个激励阵元、第j个接收阵元到任意成像点的距离
可表示为:
后续步骤中的均可通过此步中的表达式计算得出;
步骤四:将全部阵列数据按激励、接收阵元与成像点的距离计算声波传播时间,在每个成像点位置进行聚焦,并将信号幅值叠加;因此,通过公式(1)计算全阵列在每个成像点的幅值
式中,c为声波在被测试件中传播的波速;
步骤五:将步骤四中的幅值利用波束补偿因子和衰减补偿因子进行补偿,得到补偿后的幅值可通过如下公式计算得出;
式中P为波束补偿因子,B衰减补偿因子;
其中,激励、接收阵元在任意成像点的波束补偿因子可分别表示为:
式中a为单个阵元的宽度,λ为波长,sin c函数为数学中的辛格函数,当用sin c(x)表示时,其非归一化的sin c函数表示为
激励、接收阵元在任意成像点的衰减补偿因子可表示为:
步骤六:构造子阵列;将阵元总数为32的阵列换能器中多个连续的阵元作为一个子阵列,对应的时域信号作为子阵列数据;将该阵列换能器划分为7个子阵列,每个子阵列中含有8个阵元,相邻两个子阵列间的阵元个数为4个;则第k个子阵列对应的阵元在全阵列中的序号最小值为1+4(k-1)=4k-3,最大值为1+4(k-1)+7=4k+4,其中,k=1,2,3,…,K;
步骤七:计算每个子阵列在每个成像点的幅值矢量;
根据步骤六划分好的子阵列,计算每个子阵列在每个成像点的矢量,可分为以下几步:
(1)根据每个子阵列中对应的激励、接收阵元,重复步骤四、五,公式中下标i、j的最小值为4k-3,最大值为4k+4,可得到每个子阵列在每个成像点处矢量的幅值k=1,2,3,…,7,其中,上标k表示第k个子阵列;
(2)根据每个子阵列对应的激励、接收阵元到成像点的位置向量可得到每个子阵列在每个成像点的单位方向矢量量即为每个子阵列在每个成像点处矢量的方向;该单位方向矢量在x、z方向上的分量可通过公式(6)计算得出;
其中,是第i个阵元激励超声波入射到任意成像点后经第j个阵元接收形成的法线方向的单位方向矢量,依据反射定理可知,该单位方向矢量的方向会与反射面垂直,可通过公式(7)计算得出;该单位方向矢量在x、z方向上的分量可表示:
式中,为矢量的模,具体表达式为:
(3)将上两步中得到的分向量和对应相乘即可得到每个子阵列在每个成像点的幅值矢量可通过公式(8)计算得出;该幅值矢量在x、z方向上的分向量 可表示为:
步骤八:取α=4,将步骤七中7个分向量分别在x、z方向上进行合成,得到7个子阵列在每个成像点的合成幅值矢量可通过公式(9)计算得出;该合成幅值矢量在x、z方向上的分向量可表示为:
步骤九:将步骤八中得到的进行单位化,然后乘以步骤四中的得到在32阵元下每个成像点的幅值矢量可通过公式(10)计算得出;该32阵元下的幅值矢量在x、z方向上的分向量可表示为:
式中,为矢量的模,具体表达式为:
根据上述计算可知,32阵元在每个成像点的幅值矢量的大小与步骤四中的32阵元在每个成像点幅值相等,并且幅值矢量的方向会与成像点处的反射面垂直;
步骤十:将步骤九中的全阵列在每个成像点的幅值矢量进行成像显示,可得到全阵列在每个成像点的全局矢量成像图;
步骤十一:根据步骤十中的全局矢量成像图,确定目标缺陷的位置,提取目标缺陷的局部矢量图和局部矢量图的放大图;然后将步骤四中的幅值进行分贝化,找出局部矢量图中幅值的最大值,求解最大值下降-6dB所对应的成像区域面积;最后根据局部矢量图中矢量的方向与目标缺陷垂直,依据几何关系,目标缺陷的方向可通过公式(11)计算得出;
式中,A表示在目标缺陷局部矢量图中,由幅值的最大值下降-6dB所对应的成像区域面积,arctan函数为数学中的反正切函数,θm表示矢量方向与z轴正向的夹角,即等于缺陷与x轴正向的夹角;
最后,通过局部矢量图看出,缺陷处矢量幅值较明显,用箭头表现的形式,而其它不含缺陷的地方由于幅值很小,看不出来箭头,只呈现出圆点的形式,由此判断出缺陷的位置;并且通过计算求得缺陷处的角度θm=3.21°,即表示缺陷与x轴正向的夹角为3.21°。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410191263.5A CN103969337B (zh) | 2014-05-07 | 2014-05-07 | 一种基于矢量全聚焦成像的超声阵列裂纹类缺陷方向识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410191263.5A CN103969337B (zh) | 2014-05-07 | 2014-05-07 | 一种基于矢量全聚焦成像的超声阵列裂纹类缺陷方向识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103969337A CN103969337A (zh) | 2014-08-06 |
CN103969337B true CN103969337B (zh) | 2017-02-22 |
Family
ID=51239088
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410191263.5A Withdrawn - After Issue CN103969337B (zh) | 2014-05-07 | 2014-05-07 | 一种基于矢量全聚焦成像的超声阵列裂纹类缺陷方向识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103969337B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NO346618B1 (en) * | 2015-02-24 | 2022-10-31 | Halfwave As | An apparatus and method for inspecting a pipeline |
CN104730146A (zh) * | 2015-03-11 | 2015-06-24 | 北京工业大学 | 一种用于板结构裂纹定量评价的超声阵列复合成像方法 |
CN104698089A (zh) * | 2015-03-20 | 2015-06-10 | 大连理工大学 | 一种适用于倾斜裂纹定量和成像的超声相对时间传播技术 |
CN104792874A (zh) * | 2015-04-30 | 2015-07-22 | 南京迪威尔高端制造股份有限公司 | 一种实心轴锻件内部缺陷的超声波截面图识别方法 |
EP3396359B1 (en) * | 2015-12-25 | 2021-04-28 | FUJIFILM Corporation | Damage information processing device and damage information processing method |
CN106093204A (zh) * | 2016-02-22 | 2016-11-09 | 南京兰博瑞达检测技术有限公司 | 一种用于锻件裂纹定量识别的相位加权矢量全聚焦成像装置及方法 |
CN105572230B (zh) * | 2016-02-22 | 2018-10-12 | 北京工业大学 | 一种用于裂纹类缺陷定量识别的极性加权矢量全聚焦成像方法 |
CN106093206A (zh) * | 2016-07-15 | 2016-11-09 | 国网浙江省电力公司电力科学研究院 | 一种基于斜入射纵波的焊缝超声阵列全聚焦成像方法 |
CN106802323A (zh) * | 2017-01-22 | 2017-06-06 | 北京工业大学 | 一种基于全矩阵数据的超声全聚焦成像*** |
CN108053479B (zh) * | 2017-11-23 | 2020-05-12 | 河海大学 | 一种基于扩展有限元法的多缺陷群无损识别方法 |
CN108693252B (zh) * | 2018-03-31 | 2021-04-16 | 大连交通大学 | 三角矩阵聚焦成像的工件探伤检测方法 |
CN109239191B (zh) * | 2018-09-29 | 2020-11-17 | 中国特种设备检测研究院 | 一种超声导波缺陷定位成像方法及*** |
CN111487321A (zh) * | 2020-04-14 | 2020-08-04 | 广州多浦乐电子科技股份有限公司 | 基于超声反射提升聚焦能量的全聚焦成像方法 |
CN112798686B (zh) * | 2020-12-24 | 2022-04-05 | 西安交通大学 | 一种基于频率响应函数的Lamb波相控阵聚焦成像方法 |
CN113899816B (zh) * | 2021-09-10 | 2022-06-17 | 国营芜湖机械厂 | 一种t型复合结构的超声无损检测装置及方法和r区检测方法及装置 |
CN114624241B (zh) * | 2022-03-12 | 2024-07-09 | 北京工业大学 | 一种基于改进的全聚焦算法的风机主轴裂纹成像方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1813934A1 (en) * | 2004-11-02 | 2007-08-01 | Hara, Tooru | Concrete structure crack inspection device and crack inspection method |
CN101806778A (zh) * | 2010-03-05 | 2010-08-18 | 北京工业大学 | 金属材料疲劳早期损伤非线性超声在线检测方法 |
CN101819182A (zh) * | 2010-03-18 | 2010-09-01 | 安徽理工大学 | 重构非均匀介质中缺陷形状的方法 |
CN102226783A (zh) * | 2011-03-25 | 2011-10-26 | 北京工业大学 | 一种基于振动声调制技术的管道闭合裂纹检测装置及方法 |
CN102980945A (zh) * | 2012-10-11 | 2013-03-20 | 北京工业大学 | 基于双谱分析的结构微裂纹混频非线性超声检测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8942465B2 (en) * | 2011-12-13 | 2015-01-27 | General Electric Company | Methods and systems for processing images for inspection of an object |
-
2014
- 2014-05-07 CN CN201410191263.5A patent/CN103969337B/zh not_active Withdrawn - After Issue
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1813934A1 (en) * | 2004-11-02 | 2007-08-01 | Hara, Tooru | Concrete structure crack inspection device and crack inspection method |
CN101806778A (zh) * | 2010-03-05 | 2010-08-18 | 北京工业大学 | 金属材料疲劳早期损伤非线性超声在线检测方法 |
CN101819182A (zh) * | 2010-03-18 | 2010-09-01 | 安徽理工大学 | 重构非均匀介质中缺陷形状的方法 |
CN102226783A (zh) * | 2011-03-25 | 2011-10-26 | 北京工业大学 | 一种基于振动声调制技术的管道闭合裂纹检测装置及方法 |
CN102980945A (zh) * | 2012-10-11 | 2013-03-20 | 北京工业大学 | 基于双谱分析的结构微裂纹混频非线性超声检测方法 |
Non-Patent Citations (1)
Title |
---|
周进节等.基于非轴对称激励的管道裂纹时反导波检测研究.《机械工程学报》.2012,第48卷(第16期), * |
Also Published As
Publication number | Publication date |
---|---|
CN103969337A (zh) | 2014-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103969337B (zh) | 一种基于矢量全聚焦成像的超声阵列裂纹类缺陷方向识别方法 | |
Liu et al. | Damage localization in aluminum plate with compact rectangular phased piezoelectric transducer array | |
Yu et al. | Crack imaging and quantification in aluminum plates with guided wave wavenumber analysis methods | |
Yu et al. | In situ 2-D piezoelectric wafer active sensors arrays for guided wave damage detection | |
US20090048789A1 (en) | Optimized Embedded Ultrasonics Structural Radar System With Piezoelectric Wafer Active Sensor Phased Arrays For In-Situ Wide-Area Damage Detection | |
Leleux et al. | Long range detection of defects in composite plates using Lamb waves generated and detected by ultrasonic phased array probes | |
Ostachowicz et al. | Damage localisation in plate-like structures based on PZT sensors | |
Yu et al. | In-situ optimized PWAS phased arrays for Lamb wave structural health monitoring | |
Bulavinov et al. | Sampling phased array a new technique for signal processing and ultrasonic imaging | |
CN102692453B (zh) | 一种基于非线性声学的材料无损检测方法和装置 | |
Ciampa et al. | Acoustic emission localization in complex dissipative anisotropic structures using a one-channel reciprocal time reversal method | |
CN103901109A (zh) | 一种复合绝缘子内部缺陷的相控阵超声检测装置及方法 | |
CN111323784B (zh) | 一种近场分布式声源定位方法 | |
CN103575808B (zh) | 基于多角度立体匹配的高实时定量超声检测方法 | |
CN104730146A (zh) | 一种用于板结构裂纹定量评价的超声阵列复合成像方法 | |
CN110542723B (zh) | 一种基于导波信号稀疏分解及损伤定位的两阶段损伤位置识别方法 | |
Hua et al. | Modified minimum variance imaging of Lamb waves for damage localization in aluminum plates and composite laminates | |
CN113281416B (zh) | 一种用于各向异性焊缝缺陷检测的改进相位相干成像方法 | |
CN105572230B (zh) | 一种用于裂纹类缺陷定量识别的极性加权矢量全聚焦成像方法 | |
IL272819A (en) | A method of image reconstruction based on a model | |
CN106645420B (zh) | 基于费马原理的棒材超声线型阵列成像检测方法 | |
CN103424475B (zh) | 基于相控阵超声检测的被测面轮廓提取方法 | |
CN113777166A (zh) | 一种基于超声平面波成像和时间反转算子相结合的高分辨缺陷无损检测方法 | |
Liu et al. | Ultrasonic phased array total focusing method of imaging with Rayleigh waves based on principal component analysis | |
Song et al. | Accelerated noncontact guided wave array imaging via sparse array data reconstruction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
AV01 | Patent right actively abandoned | ||
AV01 | Patent right actively abandoned |
Granted publication date: 20170222 Effective date of abandoning: 20171011 |