沿磁化方向的单元组合缺陷漏磁信号计算方法技术

本发明专利技术提出一种沿磁化方向的单元组合缺陷漏磁信号计算方法，包括以下步骤：S1：根据预设的目标缺陷尺寸，选取单元缺陷尺寸；S2：获取单元缺陷漏磁信号；S3：对单元缺陷漏磁信号进行三维傅里叶变换，得到变换后的频域信号；S4：对变换后的频域信号沿磁化方向进行平移变换操作，得到平移变换后的两个频域信号；S5：将沿磁化方向平移变换后的两个频域信号进行组合操作，得到组合后的频域信号；S6：将组合后的频域信号进行三维傅里叶反变换，得到目标缺陷的漏磁信号。本发明专利技术不需要对目标缺陷再次进行有限元等复杂计算，具有计算模型简单、速度快的优点。

Calculation method of magnetic flux leakage signal of unit combination along magnetization direction

The invention provides a unit of defects along the magnetization direction of the magnetic flux leakage signal calculation method comprises the following steps: S1, according to the target defect size preset, select unit flaw size; S2: acquisition unit of MFL signal; S3: the unit of MFL signals for three-dimensional Fu Liye transform, obtain the frequency signal after conversion S4: to transform domain signal; the translation operation along the direction of magnetization, get two frequency signal after translation; S5: the two frequency signal the direction of magnetization translation after combination operation, the frequency domain signal combination; S6: the frequency domain signal after the combination of three-dimensional Fu Liye inverse transform, magnetic flux leakage signal to obtain the target defects. The invention does not need to carry out complex calculation of target defects again by finite elements, and has the advantages of simple calculation model and fast speed.

【技术实现步骤摘要】
沿磁化方向的单元组合缺陷漏磁信号计算方法
本专利技术涉及无损检测
，特别涉及一种沿磁化方向的单元组合缺陷漏磁信号计算方法。
技术介绍
缺陷漏磁信号计算是漏磁检测与缺陷评估的基础。由于缺陷漏磁信号对缺陷尺寸的不完全映射，通常需要构建根据已知缺陷尺寸求解漏磁信号的正向模型，对模型中的缺陷尺寸参数进行迭代更新来逼近目标漏磁信号，从而实现缺陷尺寸的反演。因此如何快速准确地计算缺陷漏磁信号，对于漏磁检测与缺陷评估具有重要的意义。在现有的相关技术中，主要采用磁偶极子法和有限元法直接对目标缺陷进行漏磁信号的计算，但是在漏磁信号迭代计算过程中存在大量的重复性冗余工作，大大降低了漏磁信号的计算效率。为了提高目标缺陷漏磁信号的计算效率，可以考虑利用预先计算的已知缺陷漏磁信号去表示目标缺陷漏磁信号，从而在原理上能够大大提高目标缺陷漏磁信号的计算效率。然而，目标缺陷尺寸与其漏磁信号间存在严重的非线性关系，因此无法利用预先计算的已知缺陷漏磁信号或其组合直接去表示目标缺陷漏磁信号。目前尚未有一种方法可以根据预先计算的已知缺陷漏磁信号计算目标缺陷的漏磁信号。
技术实现思路
本专利技术旨在至少解决上述技术问题之一。为此，本专利技术的目的在于提出一种沿磁化方向的单元组合缺陷漏磁信号计算方法，该方法不需要对目标缺陷再次进行有限元等复杂计算，具有计算模型简单、速度快的优点。为了实现上述目的，本专利技术的实施例提出了一种沿磁化方向的单元组合缺陷漏磁信号计算方法，包括以下步骤：S1：根据预设的目标缺陷尺寸l，w，d，选取单元缺陷尺寸l0，w0，d0；S2：获取单元缺陷漏磁信号HE(x,y,z)；S3：对单元缺陷漏磁信号HE(x,y,z)进行三维傅里叶变换，得到变换后的频域信号FE(α,β,γ)；S4：对变换后的频域信号沿磁化方向进行平移变换操作，得到平移变换后的两个频域信号FEα-(α,β,γ)和FEα+(α,β,γ)；S5：将沿磁化方向平移变换后的两个频域信号FEα-(α,β,γ)和FEα+(α,β,γ)进行组合操作，得到组合后的频域信号FEcombine(α,β,γ)；S6：将组合后的频域信号FEcombine(α,β,γ)进行三维傅里叶反变换，得到目标缺陷的漏磁信号HT(x,y,z)。另外，根据本专利技术上述实施例的沿磁化方向的单元组合缺陷漏磁信号计算方法还可以具有如下附加的技术特征：在一些示例中，所述目标缺陷是单元缺陷沿磁化方向的组合，单元缺陷尺寸l0，w0，d0与目标缺陷尺寸l，w，d之间满足下述条件：l0＝l/2，w0＝w，d0＝d其中，l0和l是缺陷的长度，沿磁化方向，w0和w是缺陷的宽度，垂直磁化方向，d0和d是缺陷的深度，沿被测材料的厚度方向。在一些示例中，所述单元缺陷漏磁信号HE(x,y,z)是根据单元缺陷尺寸l0，w0，d0通过预设算法计算得到的，其中，所述预设算法至少包括：磁偶极子法、有限元法或神经网络法。在一些示例中，在所述S3中，在进行三维傅里叶变换时，满足如下公式：其中，α，β和γ分别是x，y和z方向上的空间频率变量。在一些示例中，在所述S4中，在沿磁化方向进行平移变换操作时，满足如下公式：在一些示例中，在所述S5中，在进行组合操作时，满足如下公式：FEcombine(α,β,γ)＝FEα-(α,β,γ)+FEα+(α,β,γ)。在一些示例中，在所述S6中，在进行三维傅里叶反变换时，满足如下公式：根据本专利技术实施例的沿磁化方向的单元组合缺陷漏磁信号计算方法，可以根据预先计算的已知缺陷的漏磁信号通过变换、组合、反变换等操作，得到沿磁化方向组合后的目标缺陷的漏磁信号，不需要重复进行有限元等正向计算，提高了缺陷漏磁信号的计算效率，具有计算模型简单、速度快的优点。本专利技术的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本专利技术的实践了解到。附图说明本专利技术的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：图1是根据本专利技术实施例的沿磁化方向的单元组合缺陷漏磁信号计算方法的流程图；图2是根据本专利技术一个具体实施例中长宽深为24mm×12mm×3.6mm的目标缺陷漏磁信号求解结果对比图；图3是根据本专利技术另一个具体实施例中长宽深为48mm×24mm×4.8mm的目标缺陷漏磁信号求解结果对比图。具体实施方式下面详细描述本专利技术的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本专利技术，而不能理解为对本专利技术的限制。在本专利技术的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本专利技术和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本专利技术的限制。此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性。在本专利技术的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本专利技术中的具体含义。以下结合附图描述根据本专利技术实施例的沿磁化方向的单元组合缺陷漏磁信号计算方法。图1是根据本专利技术一个实施例的沿磁化方向的单元组合缺陷漏磁信号计算方法的流程图。如图1所示，该方法包括以下步骤：步骤S1：根据预设的目标缺陷尺寸l，w，d，选取单元缺陷尺寸l0，w0，d0。在本专利技术的一个实施例中，目标缺陷是单元缺陷沿磁化方向的组合，单元缺陷尺寸l0，w0，d0与目标缺陷尺寸l，w，d之间满足下述条件：l0＝l/2，w0＝w，d0＝d其中，l0和l是缺陷的长度，沿磁化方向，w0和w是缺陷的宽度，垂直磁化方向，d0和d是缺陷的深度，沿被测材料的厚度方向。步骤S2：获取单元缺陷漏磁信号HE(x,y,z)。在本专利技术的一个实施例中，单元缺陷漏磁信号HE(x,y,z)是根据单元缺陷尺寸l0，w0，d0通过预设算法计算得到的，其中，预设算法至少包括：磁偶极子法、有限元法或神经网络法，即预设算法包括但不限于以上列举的几种计算方法。步骤S3：对单元缺陷漏磁信号HE(x,y,z)进行三维傅里叶变换，得到变换后的频域信号FE(α,β,γ)。具体地，在步骤S3中，在进行三维傅里叶变换时，满足如下公式：其中，α，β和γ分别是x，y和z方向上的空间频率变量。步骤S4：对变换后的频域信号沿磁化方向进行平移变换操作，得到平移变换后的两个频域信号FEα-(α,β,γ)和FEα+(α,β,γ)。具体地，在步骤S4中，在沿磁化方向进行平移变换操作时，满足如下公式：步骤S5：将沿磁化方向平移变换后的两个频域信号FEα-(α,β,γ)和FEα+(α,β,γ)进行组合操作，得到组合后的频域信号FEcombine(α,β,γ)。具体地，在步骤S5中，在进行组本文档来自技高网...

【技术保护点】
一种沿磁化方向的单元组合缺陷漏磁信号计算方法，其特征在于，包括以下步骤：S1：根据预设的目标缺陷尺寸l，w，d，选取单元缺陷尺寸l

【技术特征摘要】
1.一种沿磁化方向的单元组合缺陷漏磁信号计算方法，其特征在于，包括以下步骤：S1：根据预设的目标缺陷尺寸l，w，d，选取单元缺陷尺寸l0，w0，d0；S2：获取单元缺陷漏磁信号HE(x,y,z)；S3：对单元缺陷漏磁信号HE(x,y,z)进行三维傅里叶变换，得到变换后的频域信号FE(α,β,γ)；S4：对变换后的频域信号沿磁化方向进行平移变换操作，得到平移变换后的两个频域信号FEα-(α,β,γ)和FEα+(α,β,γ)；S5：将沿磁化方向平移变换后的两个频域信号FEα-(α,β,γ)和FEα+(α,β,γ)进行组合操作，得到组合后的频域信号FEcombine(α,β,γ)；S6：将组合后的频域信号FEcombine(α,β,γ)进行三维傅里叶反变换，得到目标缺陷的漏磁信号HT(x,y,z)。2.根据权利要求1所述的沿磁化方向的单元组合缺陷漏磁信号计算方法，其特征在于，所述目标缺陷是单元缺陷沿磁化方向的组合，单元缺陷尺寸l0，w0，d0与目标缺陷尺寸l，w，d之间满足下述条件：l0＝l/2，w0＝w，d0＝d其中，l0和l是缺陷的长度，沿磁化方向，w0和w是缺陷的宽度，垂直磁化方向，d0和d是缺陷的深度，沿被测材料的厚度方向。3.根据权利要求1所述的沿磁化方向的单元组合缺陷漏磁信号计算方法，其特征在于，所述单元缺陷漏磁信号HE(x,y,z)是根据单元缺陷尺寸l0，w0，d0通过预设算法计算得到的，其中，所述预设算法至少包括：磁偶极子法、有限元法或神经网络法。4.根据权利要求1所述的沿磁化方向的单元组合缺陷漏磁信号计算方法，其特征在于，在所述S3中，在进行三维傅里叶变换时，满足如下公式：

【专利技术属性】
技术研发人员：黄松岭，赵伟，彭丽莎，于佳，
申请(专利权)人：清华大学，
类型：发明
国别省市：北京,11