CN101949683B - 一种涡流位移检测方法 - Google Patents
一种涡流位移检测方法 Download PDFInfo
- Publication number
- CN101949683B CN101949683B CN2010102647179A CN201010264717A CN101949683B CN 101949683 B CN101949683 B CN 101949683B CN 2010102647179 A CN2010102647179 A CN 2010102647179A CN 201010264717 A CN201010264717 A CN 201010264717A CN 101949683 B CN101949683 B CN 101949683B
- Authority
- CN
- China
- Prior art keywords
- signal
- eddy current
- response
- ina
- inb
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Measurement Of Length, Angles, Or The Like Using Electric Or Magnetic Means (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
Abstract
本发明涉及一种涡流位移检测方法,包括步骤1.产生正弦激励信号:通过信号发生模块产生正弦激励信号U0并将其输入涡流探头线圈模块;步骤2.产生响应信号:正弦激励信号输入涡流探头线圈模块后,根据电磁感应原理探测涡流探头线圈模块与被测对象位移变化时所引起的涡流探头模块阻抗的变化,得到响应信号并将其输入到响应信号分离模块;步骤3.响应信号分离:响应信号分离模块将响应信号分为第一路响应子信号和第二路响应子信号等步骤。本发明的有益效果是:可以通过一组涡流位移检测装置完成对多种材料的被测对象的位移检测,减少了标定次数,提高了涡流位移检测装置的互换性,降低了成本。
Description
技术领域
本发明属于无损检测技术领域,具体涉及涡流无损位移检测技术。
背景技术
涡流无损检测技术利用电磁感应原理实现检测、测量和评估。当一随时间谐变的电流输入位于一被测对象上方的探头线圈时,首先在线圈周围产生磁场,该磁场称为一级磁场。根据电磁感应原理,当被测对象与探头线圈之间的距离足够小时,在被测对象中产生涡流,该涡流也会产生一个磁场(二级磁场),阻碍一级磁场的变化,从而引起线圈阻抗的变化。因为一级磁场和二级磁场的耦合作用与诸多因素有关:被测对象和线圈的距离、被测对象参数、检测频率等,因此该技术被开发成检测位移变化的涡流位移检测装置和评估被测对象内不连续性特征的无损检测系统。涡流位移检测装置因具有结构简单、测量范围宽、抗干扰能力强、不受油污污染等优势,在航空航天、输油管道、能源等部门的关键部件(水轮机、风力发电机或其它旋转机械)的振动性能监测或涂层厚度检测(飞机发动机叶片热障涂层、输油管道绝缘层等)中得到广泛应用。
但是根据前述涡流检测理论,现有的涡流位移检测方法的检测结果除了与被测对象表面和探头线圈模块的距离有关外,还与被测对象的电磁特性有关。因此对一种被测对象对应一次标定过程。因此当被测对象的电磁特性发生变化时,同一检测距离下位移检测方法的位移输出会因材料电磁特性的不同而不同。这种现象严重降低了涡流位移检测装置的互换性,限定了其使用范围。同时也大大增加了标定工作量,即每一个涡流位移检测装置都需要根据被测对象的电磁特性逐个标定。例如被测对象的材料为铝(Al)、铜(Cu)、1Cr18Ni9Ti、45#钢、2169#钢和A3#钢六种材料,那么现有的涡流位移检测方法则需要分别 对每种材料对应的被测对象进行标定,因此需要分别设置六组涡流位移检测装置。
发明内容
为了克服现有的涡流位移检测装置的标定工作量大的缺点,提出了一种涡流位移检测方法,采用该方法的涡流位移检测装置可以匹配多种具有不同电磁特性的被测对象。
为了实现上述目的,本发明的技术方案是:一种涡流位移检测方法,包括步骤:步骤1.产生正弦激励信号:通过信号发生模块产生正弦激励信号U0并将其输入涡流探头线圈模块;
步骤2.产生响应信号:正弦激励信号输入涡流探头线圈模块后,根据电磁感应原理探测涡流探头线圈模块与被测对象位移变化时所引起的涡流探头模块阻抗的变化,得到响应信号并将其输入到响应信号分离模块;
步骤3.响应信号分离:响应信号分离模块将响应信号分为第一路响应子信号和第二路响应子信号;
步骤4.响应信号的幅值差信号和相位差信号的提取;
步骤6.模数转换:对响应信号的幅值差信号和相位差信号进行模数转换得到幅值差信号和相位差信号的模拟量对应的数字量。
步骤7.计算得到涡流探头线圈的阻抗向量的实部和虚部;
步骤8.根据涡流探头线圈模块的阻抗向量Z0计算得到被测对象的位移。
上述步骤(3)中第一路响应子信号VinA和第二路响应子信号VinB通过如下公式获得:
上述步骤4中响应信号的幅值差信号VMAG和相位差信号VPHS的提取,通过如下公式得到:
VPHS=φ(VinA)-φ(VinB) ;
上述步骤4中响应信号的幅值差信号VMAG和相位差信号MPHS的提取,通过如下公式得到:
VPHS=φ(VinA)-φ(VinB) ;
上述步骤7中涡流探头线圈的阻抗向量Z0的实部Rb和虚部Xb通过如下公式得到:
上述步骤8具体实现过程包括:步骤81:通过实验方法得到涡流探头线圈模块的等效阻抗Zpj与位移lj之间的函数关系F(lj);步骤82:通过涡流探头线圈模块阻抗模块在任一被测对象下得到的阻抗向量Zij的实部Rb和虚部Xb计算得到其等效阻抗向量Zpj,然后将等效阻抗向量zpj带入位移lj和等效阻抗向量Zpj的函数关系F(lj)进行逆运算得到位移lj。
本发明的有益效果:本发明将通过涡流探头线圈模块阻抗计算模块提取的实部Rb和虚部Xb数据进行矢量投影法处理,消除涡流位移检测装置测量中被测对象电磁特性影响,得到与被测对象电磁特性无关的涡流位移检测装置。因此可以通过一组涡流位移检测装置完成对多种材料的被测对象的位移检测,减少了标定次数,提高了涡流位移检测装置的互换性,降低了成本。
附图说明
图1是本发明的主流程图。
图2是本发明的主结构原理图。
图3是本发明响应信号分离模块的结构原理图。
图4是本发明阻抗向量和位移关系模型坐标图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明。如图1和图2所示,公开了本发明的一种涡流位移检测方法的具体实施例,包括:
步骤1.产生正弦激励信号U0:通过信号发生模块产生正弦激励信号U0并将其输入涡流探头线圈模块,U0由公式(1)表达,
U0=A0cos(ωt+φ0) 公式(1)
公式(1)中,U0是正弦激励信号、A0是正弦激励信号的幅值、ω是角频率、t是时间变量、φ0是正弦激励信号的初始相位,其中A0和ω为预先设定的已知量;φ0为测定得到的已知量。
本实施例中信号发生模块包括了控制芯片、带通滤波器和信号传输电缆,控制芯片产生方波信号,通过带通滤波器,滤掉方波信号中高频奇次谐波信号,得到需要的正弦激励信号,通过信号传输电缆传输到涡流探头线圈模块中。在本实施例中控制芯片、带通滤波器和信号传输电缆具体可以分别采用MSP430F169型号的控制芯片、π型巴特沃斯的带通滤波器和高频传输电缆。
步骤2.产生响应信号:正弦激励信号输入涡流探头线圈模块后,根据电磁感应原理探测涡流探头线圈模块与被测对象位移变化时所引起的涡流探头模块阻抗的变化,得到响应信号并将其输入到响应信号分离模块;
步骤3.响应信号分离:响应信号分离模块通过如下公式(2)和公式(3)将响应信号分为第一路响应子信号和第二路响应子信号,
公式(2)中,VinA是第一路响应子信号、Rb是涡流探头线圈模块的阻抗的实部、jXb是涡流探头线圈模块的阻抗的虚部、R0是信号发生模块的内阻、AinA是第一路响应子信号的幅值、φinA是第一路响应子信号的初始相位;公式(3)中,VinB是第二路响应子信号、I是第二路响应子信号对应的电流值,R1是第二路响应子信号对应的等效电阻、AinB是第二路响应子信号的幅值、φinB是第二路响应子信号的初始相位。式中U0、ω和R0为已知量;AinB、AinA、φinA和φinB为可测定量;R1为预设量;Rb和jXb未知量。
如图3所示,响应信号分离模块将响应信号分为两路,一路经过电压跟随器和反相器后得到第一路响应子信号VinA,另一路经过并联的电阻R1和电流/电压转换器后再经过电压跟随器后得到第二路响应子信号VinB。
步骤4.响应信号的幅值差信号VMAG和相位差信号VPHS的提取:通过如下公式(4)和公式(5)提取幅值差信号VMAG和相位差信号VPHS,
VPHS=φ(VinA)-φ(VinB) 公式(5);
步骤5.滤波:对响应信号的幅值差信号VMAG和相位差信号VPHS进行滤波处理,从响应信号的幅值差信号和相位差信号提取模块获得的幅值差信号VMAG和相位差信号VPHS是直流电压信号,通过直流滤波电路滤除直流电压信号中的畸变信号,减小干扰信号,使电路输出的幅值差信号VMAG和相位差信号VPHS更为 平滑;本步骤为可选步骤。
步骤6.模数转换:对响应信号的幅值差信号VMAG和相位差信号VPHS进行模数转换得到幅值差信号VMAG和相位差信号VPHS的模拟量对应的数字量。
利用两路模/数转换器,分别将滤波模块中的幅值差信号VMAG和相位差信号VPHS的模拟量转化为数字量,供微处理器进行后续的数据处理。
步骤7.计算得到涡流探头线圈的阻抗向量Z0的实部Rb和虚部Xb:得到幅值差信号VMAG和相位差信号VPHS的数字量后,通过如下公式得到实部Rb和虚部Xb,
通过步骤3的公式(2)除以公式(3)推导出如下公式(6)和公式(7),
步骤8.根据涡流探头线圈模块的阻抗向量Z0计算得到被测对象的位移。其具体实现过程包括:
步骤81:通过实验方法得到涡流探头线圈模块的等效阻抗Zpj与位移lj之间的函数关系。
本发明将通过涡流探头线圈模块阻抗计算模块提取的实部Rb和虚部Xb数据进行矢量投影法处理,消除涡流位移检测装置测量中被测对象电磁特性影响,得到与被测对象电磁特性无关的涡流位移检测装置。矢量投影法消除涡流位移检测装置测量中被测对象电磁特性影响的方法原理如图4所示,i和j均为自然数,其中i的取值范围为1,2,3……n,本实施例中n取值6,j为1,2,3……m,本实施例中m取值15,i代表实验过程中具有不同电磁特性的被测对象的组数,j代表实验过程中所取不同位移的组数;lj代表第j组位移, Zij代表在第j组位移下第i种被测对象下涡流探头线圈模块的阻抗向量,此处Zij=Z0。
图中坐标系中,阻抗向量Zij的实部Rb表示坐标系的横坐标,阻抗向量Zij的虚部Xb表示坐标系的纵坐标。经过实验,我们发现在同一组位移lj下具有不同电磁特性的被测对象的阻抗向量Zij在矢量坐标图种满足线性关系(实际测试过程中只是近似满足这种线性关系,但在建立该模型和模型计算过程中则认为完全满足这种线性关系),即可以通过一条斜率直线将阻抗向量Zij连接起来。通过原点并且与位移lj对应的斜率直线垂直的直线在本坐标中被定义为投影平面P,位移lj对应的斜率直线在投影平面P上的投影被定义为投影阻抗向量Zp(lj),位移lj对应的斜率直线与纵坐标Xb的交点被定义为等效阻抗向量Zpj,阻抗向量Zij与横坐标Rb的夹角被定义为阻抗向量Zij的相位角φij,投影平面P与横坐标轴Rb的夹角被定义为夹角θ。
由图可知,阻抗向量Zij与投影阻抗向量Zp(lj)的关系如公式(8),
同理,投影阻抗向量Zp(lj)与等效阻抗向量Zpj的关系如公式(9),
Zp(lj)=Zpjcos(θ-90°) 公式(9)
联合公式(8)和公式(9),可得到如下公式(10),
上述公式中θ值的具体计算过程是采用最小二乘法加黄金分割法对同一位移、不同电磁特性的被测对象的阻抗向量Zij进行拟和得到斜率直线从而计算得到θ值,由于最小二乘法加黄金分割法是现有的数学拟合计算的常用方法,因此不再详细描述。
从图中可知,位移lj和等效阻抗向量Zpj存在一一对应的关系,即在同一位移lj下,不同电磁特性的被测对象的阻抗向量Zij所对应的等效阻抗向量Zpj是相同的,因此可以建立位移lj和等效阻抗向量Zpj的函数关系,本实施例中,位移lj可以取值15组,通过上述公式(10)可以推导得出15组等效阻抗向量Zpj,将这15组位移lj和等效阻抗向量Zpj的数据进行拟合,得到位移lj和等效阻抗向量Zpj的函数关系F(lj),这种函数关系和拟合算法相关,如公式(11)。
Zpj=F(lj) 公式(11)
通过这种方法,同一位移、不同电磁特性的被测对象的阻抗向量Zij可以被等效到同一个等效阻抗向量Zpj,从而消除了因被测对象不同电磁特性对测量结果造成的影响,实现了一次标定条件下检测多种材料。建立如公式(11)的函数关系后,如果已知等效阻抗向量Zpj就可以准确的得到被测对象的位移lj,这就是本发明实现位移检测的基本原理。
步骤82:通过涡流探头线圈模块阻抗模块在任一被测对象下得到的阻抗向量Zij的实部Rb和虚部Xb经过公式(10)计算得到其等效阻抗向量Zpj,然后将等效阻抗向量Zpj带入公式(11)的位移lj和等效阻抗向量Zpj的函数关系F(lj)进行逆运算得到位移lj。
步骤9.位移lj非线性校正:根据等效阻抗向量Zpj与位移lj的非线性关系,对被测对象位移进行非线性校正得到精确的位移值。本步骤为可选步骤,并且可以通过现有的方法进行实施,因此不再详细描述。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的实施方法,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (7)
1.一种涡流位移检测方法,其特征在于,包括:
步骤1.产生正弦激励信号:通过信号发生模块产生正弦激励信号U0并将其输入涡流探头线圈模块;
步骤2.产生响应信号:正弦激励信号输入涡流探头线圈模块后,根据电磁感应原理探测涡流探头线圈模块与被测对象位移变化时所引起的涡流探头模块阻抗的变化,得到响应信号并将其输入到响应信号分离模块;
步骤3.响应信号分离:响应信号分离模块将响应信号分为第一路响应子信号和第二路响应子信号;
步骤4.响应信号的幅值差信号和相位差信号的提取;
步骤6.模数转换:对响应信号的幅值差信号和相位差信号进行模数转换得到幅值差信号和相位差信号的模拟量对应的数字量。
步骤7.计算得到涡流探头线圈的阻抗向量的实部和虚部;
步骤8.根据涡流探头线圈模块的阻抗向量Z0计算得到被测对象的位移。
2.根据权利要求1所述的一种涡流位移检测方法,其特征在于,在步骤(4)和步骤(6)之间还包括:
步骤5.滤波:对响应信号的幅值差信号和相位差信号进行滤波处理,从响应信号的幅值差信号和相位差信号提取模块获得的幅值差信号和相位差信号是直流电压信号,通过直流滤波电路滤除直流电压信号中的畸变信号,减小干扰信号,使电路输出的幅值差信号和相位差信号更为平滑。
3.根据权利要求1或2所述的一种涡流位移检测方法,其特征在于,在步骤(8)之后还包括:
步骤9.位移非线性校正:根据等效阻抗向量与位移的非线性关系,对被测对象位移进行非线性校正得到精确的位移值。
4.根据权利要求1所述的一种涡流位移检测方法,其特征在于,在步骤(3)中第一路响应子信号VinA和第二路响应子信号VinB通过如下公式获得:
上述公式中,VinA是第一路响应子信号,Rb是涡流探头线圈模块的阻抗的实部,jXb是涡流探头线圈模块的阻抗的虚部,R0是信号发生模块的内阻,AinA是第一路响应子信号的幅值,φinA是第一路响应子信号的初始相位;VinB是第二路响应子信号,I是第二路响应子信号对应的电流值,R1是第二路响应子信号对应的等效电阻,AinB是第二路响应子信号的幅值,φinB是第二路响应子信号的初始相位;式中U0,ω和R0为已知量;AinB,AinA,φinA和φinB为可测定量;R1为预设量;Rb和iXb为未知量。
5.根据权利要求1所述的一种涡流位移检测方法,其特征在于,步骤4中响应信号的幅值差信号VMAG和相位差信号VPHS的提取,通过如下公式得到:
VPHS=φ(VinA)-φ(VinB);
上述公式中,VinA是第一路响应子信号,Rb是涡流探头线圈模块的阻抗的实部,jXb是涡流探头线圈模块的阻抗的虚部,AinA是第一路响应子信号的幅值,φinA是第一路响应子信号的初始相位;VinB是第二路响应子信号;φ(VinA)为VinA与相位φ的函数,φ(VinB)为VinB与相位φ的函数;R1是第二路响应子信号对应的等效电阻,AinB是第二路响应子信号的幅值;式中AinB,AinA,φinA和φinB为可测定量;R1为预设量;Rb和iXb为未知量。
7.根据权利要求1所述的一种涡流位移检测方法,其特征在于,步骤8具体实现过程包括:步骤81:通过实验方法得到涡流探头线圈模块的等效阻抗Zpj与位移lj之间的函数关系F(lj);步骤82:通过涡流探头线圈模块阻抗模块在任一被测对象下得到的阻抗向量Zij的实部Rb和虚部Xb计算得到其等效阻抗向量Zpj,然后将等效阻抗向量Zpj带入位移lj和等效阻抗向量Zpj的函数关系F(lj)进行逆运算得到位移lj。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102647179A CN101949683B (zh) | 2010-08-27 | 2010-08-27 | 一种涡流位移检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102647179A CN101949683B (zh) | 2010-08-27 | 2010-08-27 | 一种涡流位移检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101949683A CN101949683A (zh) | 2011-01-19 |
CN101949683B true CN101949683B (zh) | 2012-09-12 |
Family
ID=43453264
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010102647179A Expired - Fee Related CN101949683B (zh) | 2010-08-27 | 2010-08-27 | 一种涡流位移检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101949683B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102809480A (zh) * | 2011-06-03 | 2012-12-05 | 北京理工大学 | 一种用于监测单体泵滑阀运动状况的方法 |
CN102955000B (zh) * | 2012-08-24 | 2015-01-07 | 电子科技大学 | 一种自适应涡流位移检测方法 |
CN102927894B (zh) * | 2012-10-30 | 2015-06-03 | 西安理工大学 | 一种电涡流检测系统及方法 |
CN103257182B (zh) * | 2013-06-07 | 2015-11-18 | 电子科技大学 | 一种脉冲涡流缺陷定量检测方法及检测系统 |
CN103499022B (zh) * | 2013-09-30 | 2015-08-19 | 清华大学 | 一种区分管道内外表面腐蚀缺陷的传感器 |
CN104006731B (zh) * | 2014-06-12 | 2017-07-04 | 株洲时代电子技术有限公司 | 一种对中传感装置及其位移检测方法 |
CN109813211B (zh) * | 2017-07-24 | 2020-09-08 | 华中科技大学 | 一种可同时测量金属工件位移及厚度的方法 |
CN108489370A (zh) * | 2018-02-28 | 2018-09-04 | 天津职业技术师范大学 | 一种适用于铝材的电涡流测距系统及方法 |
CN108534650B (zh) * | 2018-04-04 | 2019-05-14 | 大连理工大学 | 电涡流传感器输出曲线高精标定的线性度优化方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5541510A (en) * | 1995-04-06 | 1996-07-30 | Kaman Instrumentation Corporation | Multi-Parameter eddy current measuring system with parameter compensation technical field |
FR2811076B1 (fr) * | 2000-06-29 | 2002-09-20 | Snecma Moteurs | Procede de mesure d'epaisseur de paroi d'une aube creuse |
US20040239345A1 (en) * | 2001-05-01 | 2004-12-02 | Amini Bijan K. | Measurement of coatings on metal |
JP4646835B2 (ja) * | 2006-03-07 | 2011-03-09 | 株式会社日立製作所 | 渦電流探傷による残存肉厚の評価方法及び評価装置 |
CN101726238B (zh) * | 2009-12-10 | 2011-07-06 | 西安理工大学 | 差动式脉冲电涡流位移检测装置及检测方法 |
-
2010
- 2010-08-27 CN CN2010102647179A patent/CN101949683B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101949683A (zh) | 2011-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101949683B (zh) | 一种涡流位移检测方法 | |
CN101419254B (zh) | 一种特高压输电线路参数测量系统及方法 | |
CN106567708B (zh) | 一种随钻侧向电阻率测井系统及其信号检测方法 | |
CN102590618B (zh) | 一种电网基波正序电压相位的检测方法 | |
CN101846701B (zh) | 不平衡度矢量测试仪及其调整h桥电容器组平衡的方法 | |
CN106546938A (zh) | 一种电流互感器抗直流分量和偶次谐波测试系统 | |
CN107202966B (zh) | 一种变压器绕组相间漏磁场的测量方法和系统 | |
CN102175925A (zh) | 一种用于测量并联电容器组中电容的测试装置 | |
CN102207487B (zh) | 金属材质特性快速鉴别仪系统 | |
KR101346498B1 (ko) | 신호 생성 장치 및 신호 생성 방법 | |
CN103424734A (zh) | 电子式电能表在线校准方法 | |
CN202002977U (zh) | 一种用于测量并联电容器组中电容的测试装置 | |
CN104330623A (zh) | 电力系统中正弦波信号的参数测量方法及系统 | |
CN102445595A (zh) | 一种电力系统时变功率的实时测量方法 | |
CN103605038A (zh) | 一种电子式电流互感器的保护级误差检测校验系统 | |
CN102597791B (zh) | 用于测量电气设备的损耗因子的仪器和方法 | |
CN104749547A (zh) | 新型现场电能表校验仪 | |
Lindahl et al. | Noncontact electrical system monitoring on a US Coast Guard cutter | |
Sadoughi et al. | A new approach for induction motor broken bar diagnosis by using vibration spectrum | |
CN104483577B (zh) | 一种电力线路参数准确度的测定方法 | |
CN202770590U (zh) | 加工机的动平衡刀把及其动平衡控制系统 | |
CN106405465A (zh) | 一种电压互感器测量误差校验方法 | |
CN102354952B (zh) | 半波积分累积式反时限过励磁保护方法 | |
CN103247996B (zh) | 电流互感器饱和引起二次电流畸变的补偿方法 | |
WO2008016274A1 (en) | Measuring instrument for an electric current leakage |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120912 Termination date: 20150827 |
|
EXPY | Termination of patent right or utility model |