CN104950039B - 基于非线性磁饱和脉冲涡流的铁磁管道定量无损评价方法 - Google Patents
基于非线性磁饱和脉冲涡流的铁磁管道定量无损评价方法 Download PDFInfo
- Publication number
- CN104950039B CN104950039B CN201510350720.5A CN201510350720A CN104950039B CN 104950039 B CN104950039 B CN 104950039B CN 201510350720 A CN201510350720 A CN 201510350720A CN 104950039 B CN104950039 B CN 104950039B
- Authority
- CN
- China
- Prior art keywords
- mrow
- magnetic
- pulse
- msub
- eddy current
- 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
Landscapes
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
Abstract
基于非线性磁饱和脉冲涡流的铁磁管道定量无损评价方法，首先在电磁铁线圈中通入较强直流电、形成强静磁场使被检测对象磁饱和，然后在被检测件处于磁饱和状态下对TR型脉冲涡流探头完成脉冲激励和脉冲信号检出，其次基于磁饱和脉冲涡流非线性高效正问题信号模拟方法和缺陷重构反问题算法，实现铁磁管道中局部减薄缺陷的定量无损评价；本方法具有非接触、易实现、易操作、检测效率高等优点，可广泛运用于电站、化工等结构中大量使用的铁磁材料管道容器的局部减薄缺陷的定量无损评价。
Description
技术领域
本发明涉及基于电磁方法的缺陷定量无损检测技术领域，具体涉及一种基于非线性磁饱和脉冲涡流的铁磁管道定量无损评价方法。
背景技术
为保障能源供应、调整能源结构和降低温室气体排放，我国正在大力发展核电，近几年即将达到40余座，迈入国际前列。日本福岛核电站堆芯熔毁事故凸显了核电结构安全的重要性。除全电源丧失可能导致堆芯无法冷却而引发堆芯事故外，冷却所必需的大量核级管道的局部减薄等缺陷的发生和发展也可能导致核级管道的泄漏甚至大破裂，对核能结构安全造成严重威胁。例如，日本关西电力属下的美滨核电站2004年发生爆管事故，导致4人死亡，7人严重受伤，其原因是主蒸汽管道局部管壁壁厚由10mm锐减为1.4mm，导致强度不足发生事故。依据核电标准，主蒸汽管道容许安全壁厚为4.7mm，由于管道减薄缺陷未能及时检出是导致此次重大事故发生的主要原因。在核级管道中大量存在流速测定部位、L弯头、T形接头等特殊部位，由于流动的复杂性其下游管道在服役过程中发生减薄缺陷不可避免。核级管道内面产生的减薄缺陷从机理上主要可分为两类，即流动加速腐蚀(FAC:FlowAccelerated Corrosion)和液滴冲蚀(LDI:Liquid Droplet Impingement)。虽然管道初始厚度有足够的裕度，FAC和LDI可能导致管道严重减薄，可能造成如美滨核电站甚至更加严重的安全事故。LDI缺陷容易在管道L弯部位产生，严重时甚至可形成通孔缺陷。
因此，核级管道减薄缺陷对核电站的安全造成严重威胁，需要定期进行定量无损检测以保障其安全性和经济性。
核级管道材料主要为不锈钢和碳素钢，核电结构定量无损检测主要采用超声检测(UT)方法。由于管道外部的隔热/绝缘层的存在，常规超声方法无法实现直接对内壁缺陷的定量检测，而常规涡流检测方法不需要接触及耦合剂，同时脉冲涡流检测方法(PECT)对于厚壁结构无损检测具有优势。对于不锈钢管道，理论和实验研究已证明了PECT方法对于管道局部减薄缺陷检测的有效性。但对于强磁性碳素钢管道，由于磁导率影响大大降低涡流的趋肤深度，导致PECT方法无法直接应用于厚壁管道检测，需要开发新的方法。
磁饱和现象可降低材料磁导率，磁饱和涡流检测已在磁性材料检测中广泛使用，但无法实现厚板内面缺陷检测。结合脉冲涡流和磁饱和涡流的特点，磁饱和脉冲涡流检测有望克服上述困难，实现碳素钢核级管道内面减薄缺陷的定量检测，国内外尚无相关研究报道。
鉴于此，本发明提出了基于磁饱和脉冲涡流的铁磁管道定量无损评价方法，利用产生强静磁场的电磁铁、TR型脉冲涡流探头和非线性磁饱和脉冲涡流正反问题算法，实现铁磁材料管道的缺陷定量无损检测。
发明内容
为了解决上述现有技术存在的问题，本发明的目的在于提供一种基于非线性磁饱和脉冲涡流的铁磁管道定量无损评价方法，通过在电磁铁的线圈中通入较强直流电、形成强静磁场，然后在被检测件处于磁饱和状态下由脉冲涡流装置对TR型脉冲涡流探头完成脉冲激励和脉冲信号检出；再附加以非线性磁饱和脉冲涡流法的高效正问题信号模拟方法和基于逆问题的缺陷重构算法，实现铁磁管道中局部减薄缺陷的定量无损评价，可广泛运用于电站、化工等结构中大量使用的铁磁材料管道容器的局部减薄缺陷的定量检测。具有非接触、易实现、易操作、检测效率高等优点。
为达到以上目的，本发明采用如下技术方案：
基于非线性磁饱和脉冲涡流的铁磁管道定量无损评价方法，包括如下步骤：
步骤1：选取饱和磁感应强度大于1.5T的材料制作成U型电磁铁磁轭1，然后在U型电磁铁磁轭1的两腿部绕制线径为1mm～2mm的磁轭线圈2作为电磁铁的磁轭线圈；
步骤2：绕制激励线圈3和检出线圈4分别形成激励探头和检出探头，置于U型电磁铁磁轭1内部，根据信噪比取优标准优化探头间距，制作TR型脉冲涡流探头；
步骤3：搭建磁饱和脉冲涡流实验***，主要包括四部分，首先是脉冲涡流装置，由脉冲信号发生器和功率放大器组成，脉冲信号发生器用于产生方波脉冲激励信号，功率放大器用来放大脉冲激励信号；其次是磁饱和脉冲涡流探头，包括电磁铁和TR型脉冲涡流探头，其中TR型脉冲涡流探头又包括激励探头和检出探头，在电磁铁的磁轭线圈中通入强直流电，形成强静磁场，使被检测铁磁性材料处于磁饱和状态，降低其磁导率，增强涡流的趋肤深度，在此状态下，通过脉冲涡流装置对TR型脉冲涡流探头的激励探头进行脉冲信号激励，并通过检出探头接收脉冲检出信号；另外还包括提供强直流电的直流电源，和采集与分析检出探头信号的数据采集分析***；
步骤4：在被检测铁磁性材料处于磁饱和状态下，对TR型脉冲涡流探头进行方波形式的脉冲信号激励和脉冲信号检出，测量不同缺陷大小的磁饱和脉冲涡流信号，方波脉冲激励信号基频f0为100Hz，占空比为50％；对TR型脉冲涡流探头进行二维平面或曲面扫描，得到C扫图进行成像显示；
步骤5：首先基于模拟极化法计算电磁铁产生的强静磁场空间分布，然后根据被检测材料的B-H曲线计算被检测对象磁导率的三维空间分布；基于退化磁矢量位法Ar及Crank-Nicholson时域积分法计算具有三维磁导率分布的非线性脉冲涡流信号，并结合模拟极化法计算得到的强静磁场空间分布，采用非线性磁饱和脉冲涡流法对铁磁性材料检测的正问题信号进行数值模拟；
首先，基于模拟极化法计算电磁铁产生的静磁场分布的方法如下：
应用库仑规范条件模拟极化法的控制方程写作：
其中，ΩFEM是指包含铁磁性材料和磁轭的铁磁性区域，ΩBEM是指包含激励线圈的空气区域，J0是激励线圈中通入的电流源密度，M是铁磁性区域的磁化场，A是矢量磁位，μ0是空气磁导率；
方程(3)给出了铁磁性体的非线性磁化本构关系，
其中，B0是磁感应强度，H是磁场强度；
基于上述方程(1)、(2)、(3)和有限元-边界元混合法能够计算得到电磁铁产生的静磁场分布即磁场强度H的空间分布，进而得到磁感应强度B0的空间分布；
其次，利用Crank-Nicholson时域积分法求解脉冲瞬态激励的响应问题，基于时间域的逐步积分方法如下：
基于退化磁矢量位法的低频涡流场控制方程简写成，
式中：
[K]，[C]，[M]——全局系数矩阵；
I(t)——依赖于时间函数的脉冲瞬态激励电流；
对于瞬态问题，矢量磁位A关于时间的导数由关于时间步的差分(Ak-Ak-1)/Δt近似计算，其中Ak＝A(t0+kΔt)，k表示第k次时间步，Δt表示时间步长，t0表示初始时间；为了提高计算稳定性，矢量磁位A由Crank-Nicholson直接积分法
A＝θAk-1+(1-θ)Ak (5)
替换，其中(0≤θ≤1)是控制积分稳定性的参数；方程(5)代入(4)，当前步的矢量磁位A通过式(6)计算出来；
[[K](1-θ)Δt+[C]]{Ak}＝ (6)
Δt{M}I(t0+kΔt)+[[C]-[K]θΔt]{Ak-1}
通过公式(6)计算得到矢量磁位A，然后计算得到瞬态磁感应强度B′；
最后，通过结合模拟极化法计算得到的具有空间分布的强静磁感应强度B0和脉冲瞬态激励产生的瞬态磁感应强度B′，得到每一个时间积分步的磁感应强度的总和B，即B＝B0+B′；
其中，控制方程中的单元刚度矩阵与单元磁导率μ相关，因此将所有导体单元的单元磁导率μ定义为一个具有三维空间分布的数组，以每个单元中心的磁感应强度总和B作为判定参数，满足一定条件则对该单元的单元磁导率μ进行重新赋值，然后进入下一个积分时间步；在每一个积分时间步中，导入碳素钢磁导率和磁感应强度总和B的非线性磁化曲线，通过非线性磁化曲线对单元磁导率μ进行动态赋值，重新赋值后的单元磁导率μ进行单元刚度阵组装并进行下一次的积分时间步的信号计算，直至所有积分时间步计算完成，最后得到并输出所有时间步的检出信号；
基于上述步骤能够计算得到磁饱和条件下的脉冲瞬态激励的响应信号；
步骤6：结合步骤5，根据反问题中的共轭梯度算法，把求解缺陷尺寸的优化问题转化为求解目标函数的最小值问题，即通过公式(7)，
得到管道局部缺陷的尺寸c；其中，l和m代表二维扫描方向上的扫描点，L和M分别为二维扫描方向上的总扫描点个数，c是缺陷的尺寸参数，ε(c)是目标函数或者称为残差函数，Pl,m(c)是当缺陷尺寸为c时从第(l,m)个扫描点的计算得到的脉冲瞬态响应信号中提取的特征参数，是对应的从第(l,m)个扫描点的实验得到的脉冲瞬态响应信号中提取的特征参数；
利用步骤3搭建的磁饱和脉冲涡流实验***，测量不同管道缺陷的实验信号，对其进行反演，得到管道局部缺陷的实际尺寸c。
所述U型电磁铁磁轭1的饱和磁感应强度大于被检测铁磁性材料的饱和磁感应强度，其材料采用坡莫和金。
所述TR型脉冲涡流探头的激励探头和检出探头的间距为40mm。
和现有技术相比，本发明的优点如下：
1)本发明方法将磁饱和的原理引入到了脉冲涡流无损检测中。针对电站中铁磁性材料结构的局部减薄缺陷问题，结合脉冲涡流和磁饱和涡流的优点，解决了铁磁性材料的涡流的较小趋肤深度的难点。而且本方法操作简单、容易实现，因此可广泛的运用于电站结构或压力容器中大量使用的铁磁性材料的减薄缺陷的检测。
2)本发明基于开发的非线性磁饱和脉冲涡流正问题信号模拟方法和改良的常规反问题算法，实现了铁磁性材料缺陷的定量重构。
附图说明
图1为本发明方法所用***各组件连接示意图。
图2为本发明方法所用磁饱和脉冲涡流探头示意图。
图3为本发明中用到的模拟大口径管道的平板局部减薄矩形沟缺陷示意图(即被检测铁磁性材料)。
具体实施方式
如图1、图2和图3所示，本发明方法的检测步骤为：在磁轭线圈2中通入强直流电、形成强静磁场，使被检测铁磁性材料5处于磁饱和状态，降低其磁导率，增强涡流的趋肤深度。在被检测铁磁性材料5处于磁饱和状态下，通过脉冲涡流装置对TR型脉冲涡流探头的激励探头即激励线圈3进行脉冲信号激励，根据电磁感应定律，会于被检测铁磁性材料5内产生瞬态涡流，缺陷的存在会影响涡流的分布，从而影响涡流产生的瞬态次生磁场的分布，检出探头即检出线圈4两端电压会发生相应的变化，由数据采集分析***采集得到该线圈的检出电压信号，分析检出信号，结合开发的非线性磁饱和脉冲涡流正问题信号模拟方法和缺陷重构的逆问题算法，从而得到被检测铁磁性材料5内部的情况。
下面结合图1、图2、图3和具体实施例对本发明作进一步的详细描述。
本发明基于非线性磁饱和脉冲涡流的铁磁管道定量无损评价方法，包括如下步骤：
步骤1：如图2所示，选取具有较大饱和磁感应强度的材料制作U型电磁铁磁轭1，然后在U型电磁铁磁轭1的两腿部绕制线径较粗(线径选为1mm)的磁轭线圈2；
步骤2：如图2所示，绕制激励线圈3和检出线圈4，线圈的内直径均为5mm、外直径均为10mm、高均为5mm，根据信噪比取优标准优化探头间距，制作TR型脉冲涡流探头；
步骤3：如图1所示，搭建磁饱和脉冲涡流实验***，主要包括四部分，首先是脉冲涡流装置，由脉冲信号发生器和功率放大器组成，脉冲信号发生器用于产生方波脉冲激励信号，功率放大器用来放大脉冲激励信号；其次是磁饱和脉冲涡流探头，包括电磁铁和TR型脉冲涡流探头，其中TR型脉冲涡流探头又包括激励探头和检出探头，在电磁铁的磁轭线圈中通入直流电(电流选为20安培)，形成强静磁场，使被检测铁磁性材料5处于磁饱和状态，降低其磁导率，增强涡流的趋肤深度，在此状态下，通过脉冲涡流装置对TR型脉冲涡流探头的激励探头进行脉冲信号激励，并通过检出探头接收脉冲检出信号；另外还包括提供强直流电的直流电源，和采集/分析检出探头信号的数据采集与分析***；
步骤4：在被检测铁磁性材料5处于磁饱和状态下，对TR型脉冲涡流探头进行方波形式的脉冲信号激励和脉冲信号检出，测量不同缺陷大小的磁饱和脉冲涡流信号，方波脉冲激励信号的幅值为1安培，基频f0为100Hz，占空比为50％；对脉冲涡流探头进行二维平面或曲面扫描，得到C扫图进行成像显示；
步骤5：首先基于模拟极化法计算电磁铁产生的强静磁场空间分布，然后根据被检测材料的B-H曲线计算被检测对象磁导率的三维空间分布；基于退化磁矢量位法Ar及Crank-Nicholson时域积分法计算具有三维磁导率分布的非线性脉冲涡流信号，并结合模拟极化法计算得到的强静磁场空间分布，采用非线性磁饱和脉冲涡流法对铁磁性材料检测的正问题信号进行数值模拟；
首先，基于模拟极化法计算电磁铁产生的静磁场分布的方法如下：
应用库仑规范条件模拟极化法的控制方程写作：
其中，ΩFEM是指包含铁磁性材料和磁轭的铁磁性区域，ΩBEM是指包含激励线圈的空气区域，J0是激励线圈中通入的电流源密度，M是铁磁性区域的磁化场，A是矢量磁位，μ0是空气磁导率；
方程(3)给出了铁磁性体的非线性磁化本构关系，
其中，B0是磁感应强度，H是磁场强度；
基于上述公式(1)、(2)(3)和有限元-边界元混合法能够计算得到电磁铁产生的静磁场分布即磁场强度H的空间分布，进而得到磁感应强度B0的空间分布；
其次，利用Crank-Nicholson时域积分法求解脉冲瞬态激励的响应问题，基于时间域的逐步积分方法如下：
基于退化磁矢量位法的低频涡流场控制方程简写成，
式中：
[K]，[C]，[M]——全局系数矩阵；
I(t)——依赖于时间函数的脉冲瞬态激励电流；
对于瞬态问题，矢量磁位A关于时间的导数由关于时间步的差分(Ak-Ak-1)/Δt近似计算，其中Ak＝A(t0+kΔt)，k表示第k次时间步，Δt表示时间步长，t0表示初始时间；为了提高计算稳定性，矢量磁位A由Crank-Nicholson直接积分法
A＝θAk-1+(1-θ)Ak (5)替换，其中(0≤θ≤1)是控制积分稳定性的参数；方程(5)代入(4)，当前步的矢量磁位A通过式(6)计算出来；
[[K](1-θ)Δt+[C]]{Ak}＝ (6)
Δt{M}I(t0+kΔt)+[[C]-[K]θΔt]{Ak-1}
通过公式(6)计算得到矢量磁位A，然后计算得到瞬态磁感应强度B′；
最后，通过结合模拟极化法计算得到的具有空间分布的强静磁感应强度B0和脉冲瞬态激励产生的瞬态磁感应强度B′，得到每一个时间积分步的磁感应强度的总和B，即B＝B0+B′；
其中，控制方程中的单元刚度矩阵与单元磁导率μ相关，因此将所有导体单元的单元磁导率μ定义为一个具有三维空间分布的数组，以每个单元中心的磁感应强度总和B作为判定参数，满足一定条件则对该单元的单元磁导率μ进行重新赋值，然后进入下一个积分时间步；在每一个积分时间步中，导入碳素钢磁导率和磁感应强度总和B的非线性磁化曲线，通过非线性磁化曲线对单元磁导率μ进行动态赋值，重新赋值后的单元磁导率μ进行单元刚度阵组装并进行下一次的积分时间步的信号计算，直至所有积分时间步计算完成，最后得到并输出所有时间步的检出信号；
基于上述步骤能够计算得到磁饱和条件下的脉冲瞬态激励的响应信号。
步骤6：结合步骤5，根据反问题中的共轭梯度算法，把求解缺陷尺寸的优化问题转化为求解目标函数的最小值问题，即通过公式(7)，
得到管道局部缺陷的尺寸c。其中，l和m代表二维扫描方向上的扫描点，L和M分别为二维扫描方向上的总扫描点个数，c是缺陷的尺寸参数，ε(c)是目标函数(或者称为残差函数)。Pl,m(c)是当缺陷尺寸为c时从第(l,m)个扫描点的计算得到的脉冲瞬态响应信号中提取的特征参数，是对应的从第(l,m)个扫描点的实验得到的脉冲瞬态响应信号中提取的特征参数。
利用步骤3搭建的磁饱和脉冲涡流实验***，测量不同管道缺陷的实验信号，对其进行反演，得到管道局部缺陷的实际尺寸c。
作为本发明的优选实施方式，所述U型电磁铁磁轭1的饱和磁感应强度大于被检测铁磁性材料的饱和磁感应强度，其材料采用坡莫和金。
作为本发明的优选实施方式，所述TR型脉冲涡流探头的激励探头和检出探头的间距需要根据信噪比取优标准优化选取，这里为40mm。
Claims (3)
1.基于非线性磁饱和脉冲涡流的铁磁管道定量无损测量方法，其特征在于：包括如下步骤：
步骤1：选取饱和磁感应强度大于1.5T的材料制作成U型电磁铁磁轭(1)，然后在U型电磁铁磁轭(1)的两腿部绕制线径为1mm～2mm的磁轭线圈(2)作为电磁铁的磁轭线圈；
步骤2：绕制激励线圈(3)和检出线圈(4)分别形成激励探头和检出探头，置于U型电磁铁磁轭(1)内部，根据信噪比取优标准优化探头间距，制作TR型脉冲涡流探头；
步骤3：搭建磁饱和脉冲涡流实验***，该实验***包括四部分：脉冲涡流装置，由脉冲信号发生器和功率放大器组成，脉冲信号发生器用于产生方波脉冲激励信号，功率放大器用来放大脉冲激励信号；磁饱和脉冲涡流探头，包括电磁铁和TR型脉冲涡流探头，其中TR型脉冲涡流探头包括激励探头和检出探头，在电磁铁的磁轭线圈中通入强直流电，形成强静磁场，使被检测铁磁性材料处于磁饱和状态，降低其磁导率，增强涡流的趋肤深度，在此状态下，通过脉冲涡流装置对TR型脉冲涡流探头的激励探头进行脉冲信号激励，并通过检出探头接收脉冲检出信号；提供强直流电的直流电源，和采集与分析检出探头信号的数据采集分析***；
步骤4：在被检测铁磁性材料处于磁饱和状态下，对TR型脉冲涡流探头进行方波形式的脉冲信号激励和脉冲信号检出，测量不同缺陷大小的磁饱和脉冲涡流信号，方波脉冲激励信号基频f0为100Hz，占空比为50％；对TR型脉冲涡流探头进行二维平面或曲面扫描，得到C扫图进行成像显示；
步骤5：首先基于模拟极化法计算电磁铁产生的强静磁场空间分布，然后根据被检测材料的B-H曲线计算被检测对象磁导率的三维空间分布；基于退化磁矢量位法Ar及Crank-Nicholson时域积分法计算具有三维磁导率分布的非线性脉冲涡流信号，并结合模拟极化法计算得到的强静磁场空间分布，采用非线性磁饱和脉冲涡流法对铁磁性材料检测的正问题信号进行数值模拟；
首先，基于模拟极化法计算电磁铁产生的静磁场分布的方法如下：
应用库仑规范条件模拟极化法的控制方程写作：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mi>&Delta;</mi>
<mi>A</mi>
<mo>=</mo>
<mo>&dtri;</mo>
<mo>&times;</mo>
<mi>M</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</mtd>
<mtd>
<msub>
<mi>&Omega;</mi>
<mrow>
<mi>F</mi>
<mi>E</mi>
<mi>M</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mi>&Delta;</mi>
<mi>A</mi>
<mo>=</mo>
<msub>
<mi>J</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</mtd>
<mtd>
<msub>
<mi>&Omega;</mi>
<mrow>
<mi>B</mi>
<mi>E</mi>
<mi>M</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，ΩFEM是指包含铁磁性材料和磁轭的铁磁性区域，ΩBEM是指包含激励线圈的空气区域，J0是激励线圈中通入的电流源密度，M是铁磁性区域的磁化场，A是矢量磁位，μ0是空气磁导率；
方程(3)给出了铁磁性体的非线性磁化本构关系，
<mrow>
<mi>M</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>B</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msub>
<mi>B</mi>
<mn>0</mn>
</msub>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>-</mo>
<mi>H</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>B</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，B0是磁感应强度，H是磁场强度；
基于上述方程(1)、(2)、(3)和有限元-边界元混合法能够计算得到电磁铁产生的静磁场分布即磁场强度H的空间分布，进而得到磁感应强度B0的空间分布；
其次，利用Crank-Nicholson时域积分法求解脉冲瞬态激励的响应问题，基于时间域的逐步积分方法如下：
基于退化磁矢量位法的低频涡流场控制方程简写成，
<mrow>
<mo>&lsqb;</mo>
<mi>K</mi>
<mo>&rsqb;</mo>
<mo>{</mo>
<mi>A</mi>
<mo>}</mo>
<mo>+</mo>
<mo>&lsqb;</mo>
<mi>C</mi>
<mo>&rsqb;</mo>
<mo>{</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>A</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>}</mo>
<mo>=</mo>
<mo>&lsqb;</mo>
<mi>M</mi>
<mo>&rsqb;</mo>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：
[K]，[C]，[M]——全局系数矩阵；
I(t)——依赖于时间函数的脉冲瞬态激励电流；
对于瞬态问题，矢量磁位A关于时间的导数由关于时间步的差分(Ak-Ak-1)/Δt近似计算，其中Ak＝A(t0+kΔt)，k表示第k次时间步，Δt表示时间步长，t0表示初始时间；为了提高计算稳定性，矢量磁位A由Crank-Nicholson直接积分法
A＝θAk-1+(1-θ)Ak (5)
替换，其中θ是控制积分稳定性的参数，0≤θ≤1；方程(5)代入(4)，当前步的矢量磁位A通过公式(6)计算出来；
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mo>&lsqb;</mo>
<mo>&lsqb;</mo>
<mi>K</mi>
<mo>&rsqb;</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>+</mo>
<mo>&lsqb;</mo>
<mi>C</mi>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
<mo>{</mo>
<msup>
<mi>A</mi>
<mi>k</mi>
</msup>
<mo>}</mo>
<mo>=</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>{</mo>
<mi>M</mi>
<mo>}</mo>
<mi>I</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mi>k</mi>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mo>&lsqb;</mo>
<mo>&lsqb;</mo>
<mi>C</mi>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>&lsqb;</mo>
<mi>K</mi>
<mo>&rsqb;</mo>
<mi>&theta;</mi>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>&rsqb;</mo>
<mo>{</mo>
<msup>
<mi>A</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
通过公式(6)计算得到矢量磁位A，然后计算得到瞬态磁感应强度B′；
最后，通过结合模拟极化法计算得到的具有空间分布的强静磁感应强度B0和脉冲瞬态激励产生的瞬态磁感应强度B′，得到每一个时间积分步的磁感应强度的总和B，即B＝B0+B′；
其中，控制方程中的单元刚度矩阵与单元磁导率μ相关，因此将所有导体单元的单元磁导率μ定义为一个具有三维空间分布的数组，以每个单元中心的磁感应强度总和B作为判定参数，满足一定条件则对该单元的单元磁导率μ进行重新赋值，然后进入下一个积分时间步；在每一个积分时间步中，导入碳素钢磁导率和磁感应强度总和B的非线性磁化曲线，通过非线性磁化曲线对单元磁导率μ进行动态赋值，重新赋值后的单元磁导率μ进行单元刚度阵组装并进行下一次的积分时间步的信号计算，直至所有积分时间步计算完成，最后得到并输出所有时间步的检出信号；
基于上述步骤能够计算得到磁饱和条件下的脉冲瞬态激励的响应信号；
步骤6：结合步骤5，根据反问题中的共轭梯度算法，把求解缺陷尺寸的优化问题转化为求解目标函数的最小值问题，即通过公式(7)，
<mrow>
<mi>&epsiv;</mi>
<mrow>
<mo>(</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>L</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mi>m</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>c</mi>
<mo>)</mo>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mi>m</mi>
</mrow>
<mrow>
<mi>o</mi>
<mi>b</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>/</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>L</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mi>m</mi>
</mrow>
<mrow>
<mi>o</mi>
<mi>b</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
得到管道局部缺陷的尺寸c；其中，l和m代表二维扫描方向上的扫描点，L和M分别为二维扫描方向上的总扫描点个数，c是缺陷的尺寸参数，ε(c)是目标函数或者称为残差函数，Pl,m(c)是当缺陷尺寸为c时从第(l,m)个扫描点的计算得到的脉冲瞬态响应信号中提取的特征参数，是对应的从第(l,m)个扫描点的实验得到的脉冲瞬态响应信号中提取的特征参数；
利用步骤3搭建的磁饱和脉冲涡流实验***，测量不同管道缺陷的实验信号，对实验信号进行反演，得到管道局部缺陷的实际尺寸c。
2.根据权利要求1所述的基于非线性磁饱和脉冲涡流的铁磁管道定量无损测量方法，其特征在于：所述U型电磁铁磁轭(1)的饱和磁感应强度大于被检测铁磁性材料的饱和磁感应强度，其材料采用坡莫合金。
3.根据权利要求1所述的基于非线性磁饱和脉冲涡流的铁磁管道定量无损测量方法，其特征在于：所述TR型脉冲涡流探头的激励探头和检出探头的间距为40mm。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510350720.5A CN104950039B (zh) | 2015-06-23 | 2015-06-23 | 基于非线性磁饱和脉冲涡流的铁磁管道定量无损评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510350720.5A CN104950039B (zh) | 2015-06-23 | 2015-06-23 | 基于非线性磁饱和脉冲涡流的铁磁管道定量无损评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104950039A CN104950039A (zh) | 2015-09-30 |
| CN104950039B true CN104950039B (zh) | 2017-09-12 |
Family
ID=54164852
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510350720.5A Active CN104950039B (zh) | 2015-06-23 | 2015-06-23 | 基于非线性磁饱和脉冲涡流的铁磁管道定量无损评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104950039B (zh) |
Families Citing this family (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104359389B (zh) * | 2014-11-05 | 2017-04-26 | 北京航空航天大学 | 一种测量铁磁构件壁厚相对变化量的脉冲涡流检测方法 |
| CN106770632B (zh) * | 2015-11-24 | 2023-06-16 | 核动力运行研究所 | 一种适用于ω焊缝的基于收发式线圈的直流磁化探头 |
| CN105427324B (zh) * | 2015-12-07 | 2018-03-09 | 电子科技大学 | 基于二值化阈值自动搜索的磁光成像缺陷检测方法 |
| CN105546359B (zh) * | 2016-01-20 | 2017-10-03 | 合肥中大检测技术有限公司 | 便携式承压工业管道在线恒磁检测装置 |
| CN105823797B (zh) * | 2016-03-21 | 2019-07-12 | 电子科技大学 | 一种基于共生式磁轭线圈的感应热像无损检测装置 |
| CN106290561A (zh) * | 2016-10-11 | 2017-01-04 | 武汉华宇目检测装备有限公司 | 一种基于涡流磁导率测量的钢管外壁裂纹检测方法 |
| CN106289042A (zh) * | 2016-10-11 | 2017-01-04 | 武汉华宇目检测装备有限公司 | 一种钢管剩余壁厚磁化脉冲涡流测量方法与装置 |
| CN106442713A (zh) * | 2016-10-11 | 2017-02-22 | 武汉华宇目检测装备有限公司 | 一种基于表层磁导率差动测量的厚壁钢管内裂纹检测方法 |
| CN106404900A (zh) * | 2016-11-02 | 2017-02-15 | 中国计量大学 | 一种钢板表面缺陷检测装置 |
| CN108267504B (zh) * | 2016-12-30 | 2021-08-24 | 爱德森（厦门）电子有限公司 | 一种铁磁性发动机壳体内叶片动态原位监测方法 |
| CN106940343A (zh) * | 2017-03-28 | 2017-07-11 | 中国人民解放军国防科学技术大学 | 一种基于阵列电磁传感器的材料微损伤检测方法及*** |
| CN107505388B (zh) * | 2017-07-25 | 2018-07-03 | 西安交通大学 | 一种柔性磁饱和脉冲涡流检测探头及检测方法 |
| CN109490406B (zh) * | 2017-09-11 | 2020-11-20 | 清华大学 | 动态磁检测***、检测方法及电磁控阵方法 |
| CN109491306B (zh) * | 2017-09-11 | 2024-01-23 | 清华大学 | 动态磁检测探头及电磁控阵方法 |
| JP2021502558A (ja) * | 2017-11-13 | 2021-01-28 | エクソンモービル リサーチ アンド エンジニアリング カンパニーＥｘｘｏｎ Ｒｅｓｅａｒｃｈ Ａｎｄ Ｅｎｇｉｎｅｅｒｉｎｇ Ｃｏｍｐａｎｙ | 非破壊材料検査の方法及びシステム |
| CN108152365B (zh) * | 2017-11-20 | 2019-07-23 | 西安交通大学 | 基于小波分析的脉冲涡流电磁超声复合无损检测方法 |
| CN109975406A (zh) * | 2017-12-28 | 2019-07-05 | 西南石油大学 | 一种油气管道风险评估方法 |
| CN108303461A (zh) * | 2017-12-29 | 2018-07-20 | 西安交通大学 | 引入多介质单元的指套管不规则磨损缺陷的涡流检测方法 |
| CN109324085B (zh) * | 2018-11-20 | 2021-08-17 | 四川沐迪圣科技有限公司 | 基于新型高效开口磁轭线圈结构的感应热像无损检测*** |
| CN110389058B (zh) * | 2019-06-21 | 2022-04-08 | 中国特种设备检测研究院 | 确定奥氏体不锈钢材料中马氏体含量的方法 |
| CN111982967A (zh) * | 2020-08-22 | 2020-11-24 | 核动力运行研究所 | 一种基于永磁铁的磁饱和脉冲涡流红外无损评价方法 |
| CN112782273B (zh) * | 2020-12-31 | 2024-02-20 | 天津因科新创科技有限公司 | 基于脉冲涡流的管道裂纹检测方法及装置 |
| CN113640373A (zh) * | 2021-06-30 | 2021-11-12 | 清华大学 | 基于预磁化和扫频脉冲涡流的工件缺陷检测方法及装置 |
| CN113671022B (zh) * | 2021-07-14 | 2023-06-23 | 南京航空航天大学 | 基于脉冲涡流检测探头线圈间距交叉点的提离测量方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4188577A (en) * | 1975-12-11 | 1980-02-12 | Magnetic Analysis Corporation | Pulse eddy current testing apparatus for magnetic materials, particularly tubes |
| CN101520435B (zh) * | 2009-04-01 | 2011-08-24 | 华中科技大学 | 对具有导磁材料保护层的构件腐蚀检测方法及装置 |
| CN103175891A (zh) * | 2013-02-28 | 2013-06-26 | 厦门大学 | 一种永磁与脉冲涡流复合的漏磁检测方法 |
| CN103868986A (zh) * | 2012-12-13 | 2014-06-18 | 上海海事大学 | 一种金属管道内表面缺陷的涡流检测探头及其检测方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101551359B (zh) * | 2009-05-05 | 2011-10-19 | 重庆大学 | 封闭铁磁壳体内部三维作用力在线检测方法 |
| CN104316594B (zh) * | 2014-11-16 | 2017-03-29 | 吉林大学 | 一种钢件缺陷的电磁无损检测装置 |
-
2015
- 2015-06-23 CN CN201510350720.5A patent/CN104950039B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4188577A (en) * | 1975-12-11 | 1980-02-12 | Magnetic Analysis Corporation | Pulse eddy current testing apparatus for magnetic materials, particularly tubes |
| CN101520435B (zh) * | 2009-04-01 | 2011-08-24 | 华中科技大学 | 对具有导磁材料保护层的构件腐蚀检测方法及装置 |
| CN103868986A (zh) * | 2012-12-13 | 2014-06-18 | 上海海事大学 | 一种金属管道内表面缺陷的涡流检测探头及其检测方法 |
| CN103175891A (zh) * | 2013-02-28 | 2013-06-26 | 厦门大学 | 一种永磁与脉冲涡流复合的漏磁检测方法 |
Non-Patent Citations (1)
| Title |
|---|
| 带包覆层管道壁厚减薄脉冲涡流检测理论与方法;徐志远;《万方数据企业知识服务平台》;20130516;1-101 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104950039A (zh) | 2015-09-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104950039B (zh) | 基于非线性磁饱和脉冲涡流的铁磁管道定量无损评价方法 | |
| CN107505388B (zh) | 一种柔性磁饱和脉冲涡流检测探头及检测方法 | |
| Xie et al. | A novel magnetic force transmission eddy current array probe and its application for nondestructive testing of defects in pipeline structures | |
| Xie et al. | Features extraction and discussion in a novel frequency-band-selecting pulsed eddy current testing method for the detection of a certain depth range of defects | |
| Xie et al. | Quantitative non-destructive evaluation of wall thinning defect in double-layer pipe of nuclear power plants using pulsed ECT method | |
| Ali et al. | Review on system development in eddy current testing and technique for defect classification and characterization | |
| CN106596712B (zh) | 一种基于缺陷深度的选频带脉冲涡流无损检测方法 | |
| Chang et al. | Effects of excitation system on the performance of magnetic-flux-leakage-type non-destructive testing | |
| Christen et al. | Three-dimensional localization of defects in stay cables using magnetic flux leakage methods | |
| CN108195931B (zh) | 金属构件损伤检测的放置式探头及可视化定量评估方法 | |
| Jin et al. | Electromagnetic stimulation of the acoustic emission for fatigue crack detection of the sheet metal | |
| CN106769771A (zh) | 一种基于低场核磁共振技术的非饱和土渗透系数的测量方法 | |
| Yusa et al. | Numerical modeling of general cracks from the viewpoint of eddy current simulations | |
| Su et al. | Experiment and simulation on testing steel plate with corrosion defects via magnetic flux leakage method | |
| Karimian et al. | Differential permeability behaviour of P9 and T22 power station Steels | |
| Kim et al. | An application of a magnetic camera for an NDT system for aging aircraft | |
| Zakaria et al. | Simulation of magnetic flux leakage (MFL) analysis using FEMM software | |
| Ou et al. | Magnetic flux leakage testing for steel plate using pot-shaped excitation structure | |
| Suresh et al. | Characterization of defects on ferromagnetic tubes using magnetic flux leakage | |
| Wu et al. | A novel magnetic testing method for the loss of metallic cross-sectional area of bridge cables | |
| Li et al. | A novel 3D simulation prediction model of mechanical properties of ferromagnetic materials via incremental permeability method | |
| Xing et al. | Magnetic evaluation of fatigue damage in train axles without artificial excitation | |
| Wang et al. | A multi-frequency strategy for reconstruction of deep stress corrosion cracks from ECT signals at multiple liftoffs | |
| Kikuchi et al. | Applicability of magnetic flux leakage method for wall thinning monitoring in nuclear power plants | |
| Cheng | Pulsed eddy current characterization of local wall-thinning |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |