一种用于锻件裂纹定量识别的相位加权矢量全聚焦成像装置及方法
【专利摘要】本发明涉及一种用于锻件裂纹定量识别的相位加权矢量全聚焦成像装置及方法，首先，对全矩阵数据的相位信息进行分析，从中提取出相位特征参数——相位一致因子；然后，将全阵列划分为若干子阵列，利用提出的相位特征参数对每一子阵列的成像幅值进行加权，得到其相位加权成像幅值矩阵，将两者相乘进一步得到每一子阵列的加权特征向量；最后，对所有子阵列的加权特征向量进行合成，并对合成特征矢量的幅值进行全局化处理，即令合成特征矢量的幅值等于全阵列所得的相位加权成像在任意成像点处的幅值，便可得到全阵列在任意成像点处的幅值矢量；最终根据缺陷的位置提取目标缺陷的局部矢量图，利用局部矢量图中矢量的方向来确定目标缺陷的方向。
【专利说明】
一种用于锻件裂纹定量识别的相位加权矢量全聚焦成像装置 及方法
技术领域
[0001] 本发明涉及一种基于相位加权的矢量全聚焦成像方法，特别是用于锻件裂纹定量 识别的超声阵列检测方法，属于无损检测领域。
【背景技术】
[0002] 锻件是国家重大技术装备和重大工程所必需的重要基础部件，其在核电站的压力 容器、发电机组的低压转子以及各种重型机械的核心部件中均得到了广泛的应用。随着国 民经济的快速发展，大型锻件的需求量也将越来越大。在实际工程应用中，由于工作环境恶 劣，大型锻件常需承受复杂的应力、冲击振动和重负载荷。同时，由于大型锻件的生产工序 多、生产周期长，故影响其质量的因素也较多，这使得大型锻件在生产过程中将不可避免地 出现这样那样的缺陷，锻件中存在缺陷将显著地降低其力学性能，甚至在使用过程中还会 发生断裂，造成严重的安全与质量事故。
[0003] 裂纹类缺陷是锻件中常见的缺陷之一，由于裂纹具有方向性，而不同方向的裂纹 对结构危害性不尽相同，当裂纹方向与结构承载相垂直时，裂纹的危害最大，结构极易发生 突然断裂。因此，及时对结构中裂纹类缺陷进行检测及方向识别，对于预测结构的寿命、及 时维修或更换零部件、保证设备安全具有重要的意义。
[0004] 超声相控阵检测技术采用由多个(一般大于16)阵元组成的阵列换能器，通过电子 技术控制各阵元的超声激励接收延时，实现声束在试件内部的定向偏转及聚焦。利用商用 相控阵检测系统可十分方便地对构件进行实时成像，实现对回波信号的A型、B型、电子B型、 C型、D型及S型实时检测成像，但在保证成像实时性的同时，其检测精度和缺陷表征能力却 极为有限。
[0005] 近年来，随着计算机技术的发展，超声相控阵后处理成像技术得到了长足的发展。 超声相控阵后处理成像技术通过对采集到的阵列数据进行离线处理，来获得高精度的成像 效果，属于一种非实时的检测成像方法。通过对相控阵检测数据进行后处理成像，可以提高 超声相控阵检测的精度和缺陷识别能力。国内外学者对相控阵全矩阵数据后处理方法进行 了大量卓有成效的研究。例如，2009年期刊论文"Phase coherence imaging"通过对阵列数 据的相位分布进行分析，提出了一种抑制旁瓣和栅瓣，提高成像分辨率的超声阵列成像方 法--相位加权成像方法，该方法只能对结构中缺陷进行定位检测，却不能实现缺陷的特 征识别。2014硕士论文"基于矢量全聚焦的超声阵列缺陷识别方法研究及其应用"中对基于 幅值的矢量全聚焦成像及其影响因素进行了深入的研究，优化出了最佳的检测参数，实现 了裂纹的方向识别、长度测定。2008年期刊论文"Defect characterization using an ultrasonic array to measure the scattering coefficientmatrix"通过对缺陷音夂身才系 数矩阵分析，提出了基于散射系数矩阵的缺陷特征识别方法，实现了对孔和裂纹的区分及 裂纹缺陷方向、大小的测定。针对常规矢量全聚焦成像方法仅利用信号幅值信息，其成像质 量受噪声影响大的特点，本发明综合利用全矩阵数据的幅值和相位信息，提出一种基于相 位加权的矢量全聚焦成像方法，实现锻件中裂纹方向识别及长度定量测量。
【发明内容】
[0006] 本发明的目的在于提出一种用于锻件裂纹定量识别的基于相位加权的矢量全聚 焦成像方法。首先，对全矩阵数据的相位信息进行分析，提取出相位特征参数一一相位一致 因子(PCF);然后，将全阵列划分为若干子阵列，利用提出的相位特征参数对各个子阵列的 成像幅值进行加权，求取加权幅值特征向量;最后，对所有子阵列的加权特征向量进行合 成，得到相位加权的矢量全聚焦成像，并从中提取出裂纹方向及尺寸等特征信息。
[0007] 本专利的优势在于综合利用检测信号的幅值和相位信息，克服了仅利用幅值信息 时，成像结果受噪声、旁瓣和栅瓣等的影响大，信噪比和分辨率低，缺陷定量误差大的缺点， 是一种高精度阵列成像方法。
[0008] 利用该方法对锻件进行检测时，其检测系统为超声相控阵检测装置，如图1所示。 该超声相控阵检测装置包括计算机1、相控阵检测仪2、相控阵换能器3、待检测试件4、活动 卡头5、卡槽6、待检测试件平台7、连接块8、丝杠9、活动支座10、连接杆11、驱动电机12、底座 13;其中，相控阵检测仪2-端与计算机1连接，另一端与相控阵换能器3连接，相控阵换能器 3与待检测试件4通过耦合介质进行耦合。
[0009] 待检测试件4放置在待检测试件平台7上，待检测试件平台7上开有垂直相交的两 条卡槽6，卡槽6内安装有活动卡头5，活动卡头5能够沿卡槽6自由调节，待检测试件4通过活 动卡头5进行卡紧固定;待检测试件平台7的底部安装有连接块8,连接块8的底部通过螺纹 配合与丝杠9相连;丝杠9安装在四组对称的活动支座10上，活动支座10与连接杆11相连接， 连接杆11为竖直方向上的升降机构，驱动电机12通过减速器与连接杆11相连接;驱动电机 12、连接杆11均安装在底座13内。
[0010] 所述相控阵换能器3的探头由32个阵元组成。
[0011] 在计算机1的控制下相控阵检测仪2中的激励/接收模块产生激励信号，通过相控 阵换能器3(32个阵元组成的线阵相控阵探头)激励出超声波信号，沿待检测试件4传播，并 通过相控阵换能器3接收反射的超声波信号，然后通过相控阵检测仪2中的信号激励/接收 模块传输到计算机1中，通过计算机中与之配套的采集软件即可采集时域信号。
[0012] 本发明提出的基于相位加权的矢量全聚焦成像方法，其基本原理在于：
[0013] 假设线性阵列换能器阵元的个数为N，那么就会得到NXN组全矩阵数据h^(t)，其 中，i表示激励阵元，j表示接收阵元。
[0014] 通过构造复解析函数HdOzHMtHjHQdt)，便可得到信号的瞬时相位，公式 如下：
[0016] 其中HMt)称为同相分量，是信号h^(t)本身;HQjt)称为正交分量，是原信号的 希尔伯特变换。
[0017] 进一步将吟归一化到[-π,π]，由于归一化后的相位具有周期性，因此相位在穿过-和的边界变化时，会出现不连续的现象。这种不连续性导致了原本在边界附近变化不大 的相位却发生了巨大的变化，为了避免此种情况的出现，在此引入辅助相位科/⑷，其公式如 下：
[0019] 这样处理之后就会使接近±31的相位转变成接近于〇。当激励阵元为i时，对N个接 收阵元计算得到的原始相位羚的和辅助相位分别计算标准差，取两者中的较小者作 为相位变化的评价值。则相位变化的评价函数如下式所示：
[0020] sf, (f) = min[σ(^,(?)),σ(?1 (0}] (3 >
[0021] 根据以上推导，进一步给出相位一致因子的定义式：
[0023] 式中，γ是调节因子;是[-π,π]之间均匀分布的标准差，σ〇~1.8。当成像点为非 聚焦点时，由于得到的Ν个相位值是不一样的，因此奸办）>0,所以PCFdtXl;当成像点恰 好为聚焦点时，得到的N个相位值是一样的，故sfi(t) = 0，进而PCFi(t) = 1，由此可知，相位 一致因子会抑制非聚焦点的幅值，而保留聚焦点的幅值。
[0024] 设F为任意成像点，其坐标为(x，z)，当激励阵元为i时，其对应的PCF因子记作PCFi (X，z)，用得到的PCFi对合成输出",(〃)=1,/λ'Ζ hU.:，进行加权;最终将加权后的幅值 进行叠加。因此，可得任意成像点处的幅值公式IPCF(X，z)，如下式所示：
[0026]式中，tij(x，z)代表声波从第i个阵元激励传播到成像点F，再被第j个阵元接收所 需要的时间，可由特定的延迟时间公式计算获得。t^(X，Z)的计算由两点间的距离除以波速 c即可得到：
[0028] 将上述相位加权成像的基本思想引入到矢量全聚焦成像中，通过构造子阵列来实 现基于相位的矢量全聚焦成像。现假设N个阵元为一个全阵列，将该阵列划分为K个子阵列， 每个子阵列中含有η个阵元(n<N)，相邻两个子阵列间的阵元个数为m(m<N)。则第k个子阵 列对应的阵元在全阵列中的序号最小值为l+m(k-l)，最大值为n+m(k-l)，其中，k=l，2， 3 · · · K〇
[0029] 利用下式可得第k个子阵列在任意成像点的特征矢量为
[0030] (7)
[0031] 其中，/怎(.〇为第k个子阵列的相位加权成像公式，#为第k个子阵列单位 方向矢量，公式如下：
[0034]式中，％是第i个阵元激励超声波入射到任意成像点后经第j个阵元接收形成的法 线方向的单位方向矢量，依据反射定理可知，该单位方向矢量的方向会与反射面垂直，具体 计算公式可表示为：
[0036]在得到每个子阵列的特征矢量后，对所有子阵列的特征矢量进行合成，以 获得准确的缺陷方向信息。若对得到的子阵列特征矢量进行简单地叠加，会使计算得到的 缺陷方向偏离反射能量最强的方向，使得裂纹方向误差较大。为了获得准确的裂纹方向信 息，应使合成后的特征矢量方向接近反射能量最强的方向。为此，在矢量求和过程中引入加 权因子，合成后的特征矢暈表示为：
[0038] 式中，α为子阵列特征矢量合成加权因子。
[0039] 为提高缺陷检测精度及对小缺陷的识别能力，对矢量合成后成像点F处的特征矢 量du二)幅值进行全局化处理，即令合成后特征矢量幅值等于相位加权成像得到的成像点 处的幅值，这样处理后，全阵列在聚焦点(x，z)处的矢量成像特征矢量可表示为：
[0041 ]矢量成像特征矢量p U3)的方向与该点处的反射面能量最强的方向平行，因此，可 提取出缺陷的方向信息。
[0042] 为实现上述目的，本发明采用的技术方案为一种基于相位加权的矢量全聚焦成像 方法，具体可以按照以下步骤实施检测，方法流程如图2所示。
[0043] 步骤一:在如图1所示的检测装置下进行实验，其中，实验所用阵列换能器中心频 率为f，阵元总个数为N，单个阵元的宽度为a，相邻两阵元的中心距离为p，超声波在被测试 件中的传播波速为c，则波长A = c/f。通过全矩阵模式采集得到时域信号h^(t)(i = l，2， 3· · ·Ν;」= 1，2,3· · ·Ν)，其中，下标i表示阵列换能器中第i个阵元激励，j表示阵列换 能器中第j个阵元接收。
[0044] 步骤二:对采集到的时域信号h^(t)进行希尔伯特变换得到信号的包络线，将此时 得到的信号称为包络信号g^(t)。
[0045] 步骤三:建立成像坐标系，如图3所示。其中，图中〇为坐标原点，X轴表示与换能器 位置平行的方向，z轴表示与换能器位置垂直的方向。定义第i个激励阵元到成像点的向量 4，.第j个接收阵元到成像点f的向量〇以向量4.为例，豸表示向量4的模，(?表示向量 爲在1轴方向上的分量，PI表示向量&在2轴方向上的分量，后续步骤公式中的其它向量 也是类似表示。
[0046]步骤四：将全部阵列数据按激励、接收阵元与成像点的传播距离计算声波传播时 间tij(X，z)，在每个成像点位置进行聚焦;并计算每个激励阵元i对应的相位一致因子PCFi， 用得到的PCFi对合成输出?Μχ,ζ)进行加权;最终将加权后的幅值进行叠加，即可利用公式 (5)得到任意成像点的幅值I p c f (X，ζ)。
[0047]步骤五:构造子阵列。将阵元总数为Ν的阵列换能器中多个连续的阵元作为一个子 阵列，对应的时域信号称为子阵列数据。将该阵列换能器划分为Κ个子阵列，每个子阵列中 含有η个阵元(η<Ν)，相邻两个子阵列间的阵元个数为m(m<N)。则第k个子阵列对应的阵元 在全阵列中的序号最小值为l+m(k-l)，最大值为n+m(k-l)，其中，k=l，2,3. . .K。
[0048] 步骤六:计算每个子阵列在任意成像点的幅值矢量
[0049] 根据步骤五划分好的子阵列，计算每个子阵列在任意成像点的幅值矢量，可分为 以下几步：
[0050] 1.根据每个子阵列中对应的激励、接收阵元，重复步骤四，公式中下标i，j的最小 值为l+m(k-l)，最大值为n+m(k-l)，利用公式(8)可得到每个子阵列在任意成像点处矢量的 幅值/#(.〇 ( & = 12,3··I )，其中，上标k表示第k个子阵列。
[0051] 2.根据每个子阵列对应的激励、接收阵元到成像点的位置向量i，利用公式 (9)可得到每个子阵列在每个成像点的单位方向矢量，
[0052] 3.将上两步中得到的《叫χ,ζ)和(u)对应相乘即可得到每个子阵列在任意成像 点处的幅值矢量#^U)，其表达式为(7)。
[0053] 步骤七:对步骤六中k个子阵列的幅值矢量#>以，4进行合成，利用公式（11)得到k 个子阵列在每个成像点的合成幅值矢量,,
[0054] 步骤八:将合成幅值矢量二)进行单位化，然后乘以步骤四中的IPCF(X，Z)，利用 公式(12)可得到在全阵列N下每个成像点的幅值矢量
[0055] 步骤九:将步骤八中的全阵列在任意成像点的幅值矢量进行成像显示，可得 到全阵列在任意成像点的全局矢量成像图。
[0056] 步骤十:根据步骤九中的矢量图，确定目标缺陷的位置，提取目标缺陷的局部矢量 图F ;然后将步骤四中的幅值IPCF(X，z)进行dB处理，找出局部矢量图中幅值的最大值，求解 最大值下降_6dB所对应的成像区域面积。最后根据局部矢量图中矢量的方向与目标缺陷垂 直，依据几何关系，目标缺陷的方向可通过公式(13)计算得出。
[0058]式中，下标X，表示合成后最终的全聚焦幅值矢量f在X方向上的分量;下标z，表示 合成后最终的全聚焦幅值矢量f在z方向上的分量;A-6dB，缺陷成像区域内，由成像幅值最大 值下降-6dB所对应的成像区域面积。arctan函数为数学中的反正切函数，0 m表示矢量与z轴 正向的夹角，即等于缺陷与X轴正向的夹角。
【附图说明】
[0059] 图1为本发明的检测装置框图
[0060] 图2为本发明方法的流程图
[0061 ]图3为本发明方法实施例中建立的成像坐标系 [0062]图4为锻件的局部矢量全聚焦图
[0063] 图5为锻件的局部矢量全聚焦图的放大图
[0064] 图中：1、计算机，2、相控阵检测仪，3、相控阵换能器，4、待检测试件，5、活动卡头， 6、卡槽，7、待检测试件平台，8、连接块，9、丝杠，10、活动支座，11、连接杆，12、驱动电机，13、 底座。
【具体实施方式】
[0065] 本发明提出的基于相位加权的矢量全聚焦成像方法，其基本原理在于：
[0066] 假设线性阵列换能器阵元的个数为N，那么就会得到NXN组全矩阵数据h^(t)，其 中，i表示激励阵元，j表示接收阵元。
[0067]通过构造复解析函数便可得到信号的瞬时相位，公式 如下：
[0069]其中HMt)称为同相分量，是信号h^(t)本身;HQjt)称为正交分量，是原信号的 希尔伯特变换。
[0070] 进一步将％归一化到[-m]，由于归一化后的相位具有周期性，因此相位在穿过-和η的边界变化时，会出现不连续的现象。这种不连续性导致了原本在边界附近变化不大 的相位却发生了巨大的变化，为了避免此种情况的出现，在此引入辅助相位识/(0,其公式如 下：
[0072]这样处理之后就会使接近±31的相位转变成接近于〇。当激励阵元为i时，对N个接 收阵元计算得到的原始相位灼(0和辅助相位仍/(0分别计算标准差，取两者中的较小者作 为相位变化的评价值。则相位变化的评价函数如下式所示：
[0073] = (3)
[0074] 根据以上推导，进一步给出相位一致因子的定义式：
[0076]式中，γ是调节因子;是[-π,π]之间均匀分布的标准差，σ〇~1.8。当成像点为非 聚焦点时，由于得到的N个相位值是不一样的，因此奸办）>0,所以PCFdtXl;当成像点恰 好为聚焦点时，得到的N个相位值是一样的，故sfi(t) = 0，进而PCFi(t) = 1，由此可知，相位 一致因子会抑制非聚焦点的幅值，而保留聚焦点的幅值。
[0077]设F为任意成像点，其坐标为(x，z)，当激励阵元为i时，其对应的PCF因子记作PCFi (X，z)，用得到的PCFi对合成输出乂 = ι/λ Ζ : (L-7))进行加权;最终将加权后的幅值 进行叠加。因此，可得任意成像点处的幅值公式IPCF(X，z)，如下式所示：
[0079]式中，?^(Χ，Ζ)代表声波从第i个阵元激励传播到成像点F，再被第j个阵元接收所 需要的时间，可由特定的延迟时间公式计算获得。t^(X，Z)的计算由两点间的距离除以波速 c即可得到：
[0081] 将上述相位加权成像的基本思想引入到矢量全聚焦成像中，通过构造子阵列来实 现基于相位的矢量全聚焦成像。现假设N个阵元为一个全阵列，将该阵列划分为K个子阵列， 每个子阵列中含有η个阵元(n<N)，相邻两个子阵列间的阵元个数为m(m<N)。则第k个子阵 列对应的阵元在全阵列中的序号最小值为l+m(k-l)，最大值为n+m(k-l)，其中，k=l，2， 3..,K〇
[0082] 利用下式可得第k个子阵列在任意成像点的特征矢量为#
[0083] va\x,z) = n-lt)(x^)i，A'f(^=) (7)
[0084] 其中，为第k个子阵列的相位加权成像公式，^(.w)为第k个子阵列单位 方向矢量，公式如下：
[0087]式中，是第i个阵元激励超声波入射到任意成像点后经第j个阵元接收形成的法 线方向的单位方向矢量，依据反射定理可知，该单位方向矢量的方向会与反射面垂直，具体 计算公式可表示为：
[0089]在得到每个子阵列的特征矢量后，对所有子阵列的特征矢量进行合成，以 获得准确的缺陷方向信息。若对得到的子阵列特征矢量进行简单地叠加，会使计算得到的 缺陷方向偏离反射能量最强的方向，使得裂纹方向误差较大。为了获得准确的裂纹方向信 息，应使合成后的特征矢量方向接近反射能量最强的方向。为此，在矢量求和过程中引入加 权因子，合成后的特征矢量表示为：
[0091] 式中，α为子阵列特征矢量合成加权因子。
[0092] 为提高缺陷检测精度及对小缺陷的识别能力，对矢量合成后成像点F处的特征矢 量幅值进行全局化处理，即令合成后特征矢量幅值等于相位加权成像得到的成像点 处的幅值，这样处理后，全阵列在聚焦点(x，z)处的矢量成像特征矢量可表示为：
[0094] 矢量成像特征矢量fU. 4的方向与该点处的反射面能量最强的方向平行，因此，可 提取出缺陷的方向信息。
[0095] 为实现上述目的，本发明采用的技术方案为一种基于相位加权的矢量全聚焦成像 方法，具体可以按照以下步骤实施检测，方法流程如图2所示。
[0096] 步骤一:在如图1所示的检测装置下进行实验，其中，实验所用阵列换能器中心频 率为f，阵元总个数为N，单个阵元的宽度为a，相邻两阵元的中心距离为p，超声波在被测试 件中的传播波速为c，则波长A = c/f。通过全矩阵模式采集得到时域信号h^(t)(i = l，2， 3...化」=1，2,3...?，其中，下标1表示阵列换能器中第1个阵元激励，」表示阵列换能器中 第j个阵元接收。
[0097] 步骤二:对采集到的时域信号h^(t)进行希尔伯特变换得到信号的包络线，将此时 得到的信号称为包络信号g^(t)。
[0098]步骤三:建立成像坐标系，如图3所示。其中，图中〇为坐标原点，X轴表示与换能器 位置平行的方向，z轴表示与换能器位置垂直的方向。定义第i个激励阵元到成像点的向量 4:，第j个接收阵元到成像点f的向量毛以向量藏为例，4|表示向量.4的模，(?表示向量 在X轴方向上的分量，表示向量^在2轴方向上的分量，后续步骤公式中的其它向量也是 类似表示。
[0099] 步骤四：将全部阵列数据按激励、接收阵元与成像点的传播距离计算声波传播时 间tij(X，z)，在每个成像点位置进行聚焦;并计算每个激励阵元i对应的相位一致因子PCFi， 用得到的PCFi对合成输出?Μχ,ζ)进行加权;最终将加权后的幅值进行叠加，即可利用公式 (5)得到任意成像点的幅值I p c f (X，ζ)。
[0100] 步骤五:构造子阵列。将阵元总数为Ν的阵列换能器中多个连续的阵元作为一个子 阵列，对应的时域信号称为子阵列数据。将该阵列换能器划分为Κ个子阵列，每个子阵列中 含有η个阵元(η<Ν)，相邻两个子阵列间的阵元个数为m(m<N)。则第k个子阵列对应的阵元 在全阵列中的序号最小值为l+m(k-l)，最大值为n+m(k-l)，其中，k=l，2,3. . .K。
[0101] 步骤六:计算每个子阵列在任意成像点的幅值矢量#Μχρ)。
[0102] 根据步骤五划分好的子阵列，计算每个子阵列在任意成像点的幅值矢量，可分为 以下几步：
[0103] 1.根据每个子阵列中对应的激励、接收阵元，重复步骤四，公式中下标i，j的最小 值为l+m(k-l)，最大值为n+m(k-l)，利用公式(8)可得到每个子阵列在任意成像点处矢量的 幅值45dα = 1,2,3…),其中，上标k表示第k个子阵列。
[0104] 2.根据每个子阵列对应的激励、接收阵元到成像点的位置向量^，利用公式 (9)可得到每个子阵列在每个成像点的单位方向矢量
[0105] 3.将上两步中得到的#和CUO对应相乘即可得到每个子阵列在任意成像 点处的幅值矢量# > (.V，其表达式为(7)。
[0106] 步骤七:对步骤六中k个子阵列的幅值矢量进行合成，利用公式（11)得到k 个子阵列在每个成像点的合成幅值矢量。
[0107] 步骤八:将合成幅值矢量ο?αι进行单位化，然后乘以步骤四中的IPCF(X，Z)，利用 公式(12)可得到在全阵列N下每个成像点的幅值矢量Ρ^,:ζΚ
[0108] 步骤九:将步骤八中的全阵列在任意成像点的幅值矢量?、^进行成像显示，可得 到全阵列在任意成像点的全局矢量成像图。
[0109] 步骤十:根据步骤九中的矢量图，确定目标缺陷的位置，提取目标缺陷的局部矢量 图f ;然后将步骤四中的幅值IPCF(X，z)进行dB处理，找出局部矢量图中幅值的最大值，求解 最大值下降_6dB所对应的成像区域面积。最后根据局部矢量图中矢量的方向与目标缺陷垂 直，依据几何关系，目标缺陷的方向可通过公式(13)计算得出。
[0111]式中，下标X，表示合成后最终的全聚焦幅值矢量f在X方向上的分量;下标Z，表示 合成后最终的全聚焦幅值矢量交在Z方向上的分量;A-6dB，缺陷成像区域内，由成像幅值最大 值下降-6dB所对应的成像区域面积。arctan函数为数学中的反正切函数，0 m表示矢量与z轴 正向的夹角，即等于缺陷与X轴正向的夹角。
【主权项】
1. 一种用于锻件裂纹定量识别的相位加权矢量全聚焦成像装置，其特征在于：该装置 包括计算机（1)、相控阵检测仪(2)、相控阵换能器(3)、待检测试件(4)、活动卡头(5)、卡槽 (6)、待检测试件平台（7)、连接块(8)、丝杠(9)、活动支座（10)、连接杆（11)、驱动电机(12)、 底座（13);其中，相控阵检测仪(2)-端与计算机(1)连接，另一端与相控阵换能器(3)连接， 相控阵换能器(3)与待检测试件(4)通过禪合介质进行禪合； 待检测试件(4)放置在待检测试件平台（7)上，待检测试件平台（7)上开有垂直相交的 两条卡槽(6)，卡槽(6)内安装有活动卡头(5)，活动卡头(5)能够沿卡槽(6)自由调节，待检 测试件(4)通过活动卡头(5)进行卡紧固定;待检测试件平台（7)的底部安装有连接块(8)， 连接块(8)的底部通过螺纹配合与丝杠(9)相连;丝杠(9)安装在四组对称的活动支座（10) 上，活动支座（10)与连接杆（11)相连接，连接杆（11)为竖直方向上的升降机构，驱动电机 (12)通过减速器与连接杆(11)相连接，驱动电机(12 )、连接杆(11)均安装在底座(13)内。2. 根据权利要求1所述的一种用于锻件裂纹定量识别的相位加权矢量全聚焦成像装 置，其特征在于:所述相控阵换能器(3)的探头由32个阵元组成。3. 根据权利要求1所述的一种用于锻件裂纹定量识别的相位加权矢量全聚焦成像装 置，其特征在于:在计算机（1)的控制下相控阵检测仪(2)中的激励/接收模块产生激励信 号，通过相控阵换能器(3)激励出超声波信号，沿待检测试件4传播，并通过相控阵换能器 (3)接收反射的超声波信号，然后通过相控阵检测仪(2)中的信号激励/接收模块传输到计 算机(1)中，通过计算机中与之配套的采集软件即可采集时域信号。4. 利用权利要求1所述的成像装置进行的一种用于锻件裂纹定量识别的相位加权矢量 全聚焦成像方法，该方法的基本原理如下： 假设线性阵列换能器阵元的个数为N，那么就会得到NXN组全矩阵数据hu(t)，其中，i 表示激励阵元，j表示接收阵元； 通过构造复解析函数出^*)=^4(*)+^·册u(t)，便可得到信号的瞬时相位，公式如下：(!； 其中Hlu(t)称为同相分量，是信号hu(t)本身;朋1如)称为正交分量，是原信号的希尔 伯特变换； 进一步将辉归一化到[-n，n]，由于归一化后的相位具有周期性，因此相位在穿过-31和31 的边界变化时，会出现不连续的现象;运种不连续性导致了原本在边界附近变化不大的相 位却发生了巨大的变化，为了避免此种情况的出现，在此引入辅助相位皆似，其公式如下：(2) 运样处理之后就会使接近±n的相位转变成接近于0;当激励阵元为i时，对N个接收阵 元计算得到的原始相位％(0和辅助相位巧/的分别计算标准差，取两者中的较小者作为相 位变化的评价值;则相位变化的评价函数如下式所示：(3) 根据w上推导，进一步给出相位一致因子的定义式：(4) 式中，丫是调节因子;〇日是[-K川之间均匀分布的标准差，則>1.8;当成像点为非聚焦 点时，由于得到的N个相位值是不一样的，因此sfi(t)>0,所WP化(t)<l;当成像点恰好为 聚焦点时，得到的N个相位值是一样的，故sf i(t) = 0，进而PCFi(t) = 1，由此可知，相位一致 因子会抑制非聚焦点的幅值，而保留聚焦点的幅值； 设F为任意成像点，其坐标为(x，z)，当激励阵元为i时，其对应的PCF因子记作PCFiU, Z)，用得到的PCFi对合成输出巧行加权;最终将加权后的幅值进 行叠加;因此，可得任意成像点处的幅值公式Ipcf(x，z)，如下式所示：(5J 式中，tu(x，z)代表声波从第i个阵元激励传播到成像点F，再被第j个阵元接收所需要 的时间，可由特定的延迟时间公式计算获得;tu(x，z)的计算由两点间的距离除W波速C即 可得到：梯) 将上述相位加权成像的基本思想引入到矢量全聚焦成像中，通过构造子阵列来实现基 于相位的矢量全聚焦成像;现假设N个阵元为一个全阵列，将该阵列划分为K个子阵列，每个 子阵列中含有η个阵元(n<N)，相邻两个子阵列间的阵元个数为m(m<N);则第k个子阵列对 应的阵元在全阵列中的序号最小值为1+m化-1)，最大值为n+m化-1 )，其中，k=l，2,3...K; 利用下式可得第k个子阵列在任意成像点的特征矢量为V-"'屯，Z):(7) 其中，为第k个子阵列的相位加权成像公式，心?ι(χ，3)为第k个子阵列单位方向矢 量，公式如下：式中，:?是第i个阵元激励超声波入射到任意成像点后经第j个阵元接收形成的法线方 向的单位方向矢量，依据反射定理可知，该单位方向矢量的方向会与反射面垂直，具体计算 公式可表示为：(10、 在得到每个子阵列的特征矢量佔巧后，对所有子阵列的特征矢量进行合成，w获得 准确的缺陷方向信息；若对得到的子阵列特征矢量进行简单地叠加，会使计算得到的缺陷 方向偏离反射能量最强的方向，使得裂纹方向误差较大;为了获得准确的裂纹方向信息，应 使合成后的特征矢量方向接近反射能量最强的方向；为此，在矢量求和过程中引入加权因 子，合成后的特征矢量表示为：(11) 式中，α为子阵列特征矢量合成加权因子； 为提高缺陷检测精度及对小缺陷的识别能力，对矢量合成后成像点F处的特征矢量 6佔对幅值进行全局化处理，即令合成后特征矢量幅值等于相位加权成像得到的成像点处 的幅值，运样处理后，全阵列在聚焦点(x，z)处的矢量成像特征矢量可表示为：(12' 矢量成像特征矢量的方向与该点处的反射面能量最强的方向平行，因此，可提取 出缺陷的方向信息； 其特征在于:该方法流程如下， 步骤一:在搭建该检测装置后进行实验，其中，实验所用阵列换能器中屯、频率为f，阵元 总个数为N，单个阵元的宽度为a，相邻两阵元的中屯、距离为P，超声波在被测试件中的传播 波速为C，则波长λ=c/f;通过全矩阵模式采集得到时域信号hu(t) (i = 1，2，3. . .N; j = 1，2，3...N)，其中，下标i表示阵列换能器中第i个阵元激励，j表示阵列换能器中第j个阵元接 收； 步骤二:对采集到的时域信号hu(t)进行希尔伯特变换得到信号的包络线，将此时得到 的信号称为包络信号gij(t); 步骤Ξ:建立成像坐标系；其中，0为坐标原点，X轴表示与换能器位置平行的方向，Z轴 表示与换能器位置垂直的方向；定义第i个激励阵元到成像点的向量4,第j个接收阵元到 成像点F的向量:拟向量导为例，??,表示向量马的模，(3 )、表示向量4在X轴方向上的分量， (?表示向量马在Ζ轴方向上的分量，后续步骤公式中的其它向量也是类似表示； 步骤四：将全部阵列数据按激励、接收阵元与成像点的传播距离计算声波传播时间tij (X，Z )，在每个成像点位置进行聚焦;并计算每个激励阵元i对应的相位一致因子PCFi，用得 到的PCFi对合成输出出(x，z)进行加权;最终将加权后的幅值进行叠加，即可利用公式(5)得 到任意成像点的幅值ΙροΚχ,ζ); 步骤五：构造子阵列；将阵元总数为Ν的阵列换能器中多个连续的阵元作为一个子阵 列，对应的时域信号称为子阵列数据;将该阵列换能器划分为Κ个子阵列，每个子阵列中含 有η个阵元(η<Ν)，相邻两个子阵列间的阵元个数为m(m<N);则第k个子阵列对应的阵元在 全阵列中的序号最小值为1+m化-1)，最大值为n+m化-1 )，其中，k=l，2,3...K; 步骤六:计算每个子阵列在任意成像点的幅值矢量; 根据步骤五划分好的子阵列，计算每个子阵列在任意成像点的幅值矢量，可分为W下 几步： S6.1根据每个子阵列中对应的激励、接收阵元，重复步骤四，公式中下标i，j的最小值 为1+m化-1)，最大值为n+m化-1)，利用公式(8)可得到每个子阵列在任意成像点处矢量的幅 值趨山，刮化=1,2,3...K)，其中，上标k表示第k个子阵列； S6.2根据每个子阵列对应的激励、接收阵元到成像点的位置向量4、^，利用公式(9) 可得到每个子阵列在每个成像点的单位方向矢量1万"(、·，。，）； S6.3将上两步中得到的。Z)和货f(.T，：)对应相乘即可得到每个子阵列在任意成像点 处的幅值矢量:，其表达式为(7); 步骤屯:对步骤六中k个子阵列的幅值矢量进行合成，利用公式（11)得到k个子 阵列在每个成像点的合成幅值矢量(5仁，Z); 步骤八:将合成幅值矢量巧Χ，ζ)进行单位化，然后乘W步骤四中的IPCF(X，Z)，利用公式 (12)可得到在全阵列N下每个成像点的幅值矢量:？(χ，ζ); 步骤九:将步骤八中的全阵列在任意成像点的幅值矢量fu，=)进行成像显示，可得到全 阵列在任意成像点的全局矢量成像图； 步骤十:根据步骤九中的矢量图，确定目标缺陷的位置，提取目标缺陷的局部矢量图F ; 然后将步骤四中的幅值Ipcf(x，z)进行地处理，找出局部矢量图中幅值的最大值，求解最大 值下降-6地所对应的成像区域面积;最后根据局部矢量图中矢量的方向与目标缺陷垂直， 依据几何关系，目标缺陷的方向可通过公式(13)计算得出；(13) 式中，下标X，表示合成后最终的全聚焦幅值矢量在X方向上的分量;下标Z，表示合成 后最终的全聚焦幅值矢量专在Z方向上的分量;A-6dB，缺陷成像区域内，由成像幅值最大值下 降-6地所对应的成像区域面积;arc化η函数为数学中的反正切函数，目m表示矢量与Z轴正向 的夹角，即等于缺陷与X轴正向的夹角。
【文档编号】G01N29/44GK106093204SQ201610096738
【公开日】2016年11月9日
【申请日】2016年2月22日 公开号201610096738.1, CN 106093204 A, CN 106093204A, CN 201610096738, CN-A-106093204, CN106093204 A, CN106093204A, CN201610096738, CN201610096738.1
【发明人】焦敬品, 杨素方, 何存富, 李先锋
【申请人】南京兰博瑞达检测技术有限公司, 北京工业大学