CN109668506B - 一种基于涡流稳态特性的磁性金属材料厚度检测方法 - Google Patents

一种基于涡流稳态特性的磁性金属材料厚度检测方法 Download PDF

Info

Publication number
CN109668506B
CN109668506B CN201910129231.5A CN201910129231A CN109668506B CN 109668506 B CN109668506 B CN 109668506B CN 201910129231 A CN201910129231 A CN 201910129231A CN 109668506 B CN109668506 B CN 109668506B
Authority
CN
China
Prior art keywords
eddy current
tested piece
piece
magnetic
thickness
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
Application number
CN201910129231.5A
Other languages
English (en)
Other versions
CN109668506A (zh
Inventor
田露露
白利兵
程玉华
林科培
王晓杰
周权
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201910129231.5A priority Critical patent/CN109668506B/zh
Publication of CN109668506A publication Critical patent/CN109668506A/zh
Application granted granted Critical
Publication of CN109668506B publication Critical patent/CN109668506B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B7/00Measuring arrangements characterised by the use of electric or magnetic techniques
    • G01B7/02Measuring arrangements characterised by the use of electric or magnetic techniques for measuring length, width or thickness
    • G01B7/06Measuring arrangements characterised by the use of electric or magnetic techniques for measuring length, width or thickness for measuring thickness
    • G01B7/10Measuring arrangements characterised by the use of electric or magnetic techniques for measuring length, width or thickness for measuring thickness using magnetic means, e.g. by measuring change of reluctance
    • G01B7/102Height gauges

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
  • Measurement Of Length, Angles, Or The Like Using Electric Or Magnetic Means (AREA)

Abstract

本发明公开了一种基于涡流稳态特性的磁性金属材料厚度检测方法,先通过激励线圈的磁场分布来建立测试模型,然后求取激励线圈产生的激励磁场在被测试件中产生的涡流密度分布与趋肤深度的关系,再通过涡流密度分布求解出被测试件中涡流产生的磁感应强度,最后在被测试件被磁化的状态下,建立涡流磁感应强度与被测试件的厚度关系,并通过传感器实时测量涡流磁感应强度用于求取出被测试件的厚度。

Description

一种基于涡流稳态特性的磁性金属材料厚度检测方法
技术领域
本发明属于涡流电磁场技术领域,更为具体地讲,涉及一种基于涡流稳态特性的磁性金属材料厚度检测方法。
背景技术
目前,在工业生产和设备检测中,对于金属板材或者金属设备的厚度测量,通常有超声、射线、电磁涡流等不同方法。其中,涡流无损检测技术具有灵敏度高、检测速度快、非接触等特点而广泛应用于金属厚度检测,它的原理是当载有交变电流的检测线圈靠近金属导体时在导体中产生涡流,该涡流又影响原磁场,使得线圈的阻抗和感应电压发生变化,通过分析阻抗或感应电压的变化来获得被测导体的信息.。
其中,针对基于涡流特性的测厚方法主要依据涡流的瞬态特性,即利用趋肤效应,求取趋肤深度与频率之间的关系来检测被测材料的厚度,这种方法技术比较成熟,测量精度达到μm级别,但是对基于涡流稳态特性的测厚方法研究鲜有报道或文献记载。因此,对基于涡流稳态特性的磁性金属材料的厚度检测研究,成为工业测量领域内亟待解决的问题。
发明内容
本发明的目的在于克服现有技术的不足,提供一种基于涡流稳态特性的磁性金属材料厚度检测方法,利用涡流动态特性来实时测量被测试件的厚度。
为实现上述发明目的,本发明一种基于涡流稳态特性的磁性金属材料厚度检测方法,其特征在于,包括以下步骤:
(1)、建立测试模型
在激励线圈中心插入圆柱型探针,在探针垂线方向,距离探针下端探头距离为d处设为点a,在点a处放置被测试件,在探针下端探头处设置HALL传感器;
在点a处由激励线圈产生的磁感应强度B1为:
Figure GDA0002227749950000021
其中,μ为探针的磁导率,I1为激励线圈的激励电流,r为激励线圈的半径,rp为探头的半径,n为激励线圈的匝数,
Figure GDA0002227749950000022
为绕制激励线圈材料的直径;
(2)、计算被测试件的涡流密度分布
(2.1)、计算被测试件在点a处感生出的电场E1
Figure GDA0002227749950000023
其中,S为探头下一个圆形闭合平面的面积;
(2.2)、计算点a处被测试件表面的涡流密度J;
Figure GDA0002227749950000024
其中,σ为被测试件的电导率,M为常数;
当激励线圈在正弦信号I1=Acos(ωt)的激励下,那么在点a处被测试件表面的涡流密度J表示为:
J=σMSAω·sin(ωt)
其中,A为正弦信号的幅度,ω为正弦信号的角频率;
(2.3)、计算被测试件内部的涡流密度J(h');
J(h')=J·eh'=σMSAω·sin(ωt)·eh'
其中,J(h')表示被测试件内部h'深度处的涡流密度,δ表示在角频率ω下被测试件中涡流的趋肤深度;
(3)、根据被测试件内部的涡流密度J(h')计算被测试件内部的涡流磁感应强度B2
Figure GDA0002227749950000025
其中,B2方向与B1相反,μ2为被测试件的磁导率,rs为被测试件中涡流的等效最大半径,J(z)表示试件z深度处的涡流密度;
(4)、建立稳态特性与被测试件的厚度关系
(4.1)、估算被测试件在完全被磁化时的磁导率;
被测试件在完全被磁化时满足:μH>εH0s0m0h,其中,h为被测试件厚度,s0为单位体积的底面积,ε为磁极子间摩擦衰减因子,H为稳态磁场强度,H0为单个磁极子提供的磁场强度,m0为单位体积可磁化磁极子数;
稳态特性下,稳态磁场强度与单磁极子的关系为:H=kH0,其中,k为磁场合成因子;
那么被测试件在完全被磁化时的磁导率估计值
Figure GDA0002227749950000031
为:
Figure GDA0002227749950000032
(4.2)、用
Figure GDA0002227749950000033
替换μ2,估算出被测试件内部的涡流磁感应强度B2
Figure GDA0002227749950000034
(4.3)、估算被测试件的厚度h;
Figure GDA0002227749950000035
(5)、计算待测试件的厚度
(5.1)、当只有激励线圈工作时,点a处未放置被测试件,利用HALL传感器测出激励线圈在点a处的磁感应强度B1
(5.2)、当磁场稳定时,在点a处放置被测试件,然后利用HALL传感器测出点a处的磁感应强度,记为B;B是激励线圈磁感应强度B1与涡流磁感应强度B2叠加而成,且B2与B1方向相反,故B2=B-B1
(5.3)、将测出的涡流磁感应强度B2代入至步骤(4.3)的公式中,计算出被测试件的厚度。本发明的发明目的是这样实现的:
本发明一种基于涡流稳态特性的磁性金属材料厚度检测方法,先通过激励线圈的磁场分布来建立测试模型,然后求取激励线圈产生的激励磁场在被测试件中产生的涡流密度分布与趋肤深度的关系,再通过涡流密度分布求解出被测试件中涡流产生的磁感应强度,最后在被测试件被磁化的状态下,建立涡流磁感应强度与被测试件的厚度关系,并通过传感器实时测量涡流磁感应强度用于求取出被测试件的厚度。
同时,本发明一种基于涡流稳态特性的磁性金属材料厚度检测方法还具有以下有益效果:
(1)、本发明通过建立涡流磁场与被测试件的厚度关系,推导出涡流磁感应强度与被测试件的厚度关系式,由此非常有效、快捷、简单的测量处被测试件的厚度;
(2)、本发明采用了磁场的磁化动态特性,能够实现对所有利用磁化特性的材料厚度测量;
(3)、相比其他的涡流测厚方法,本发明提供的测厚方法完成材料厚度测量,操作比较方便,误差比较小,测量结果更加准确。
附图说明
图1是本发明一种基于涡流稳态特性的磁性金属材料厚度检测方法流程图;
图2是本发明搭建的测试模型示意图;
图3是测试值与实际值的曲线图;
图4是均值的拟合曲线。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
实施例
图1是本发明一种基于涡流稳态特性的磁性金属材料厚度检测方法流程图。
在本实施例中,如图1所示,本发明一种基于涡流稳态特性的磁性金属材料厚度检测方法,包括以下步骤:
S1、建立测试模型
如图2所示,在激励线圈中心插入圆柱型探针,在探针垂线方向,距离探针下端探头距离为d处设为点a,在点a处放置被测试件,在探针下端探头处设置HALL传感器,用于测量点a处磁感应强度;在本实施例中,探针采用圆柱型铁芯。
在点a处由激励线圈产生的磁感应强度B1为:
Figure GDA0002227749950000051
其中,μ为探针的磁导率,I1为激励线圈的激励电流,r为激励线圈的半径,rp为探头的半径,n为激励线圈的匝数,
Figure GDA0002227749950000052
为绕制激励线圈材料的直径;
S2、计算被测试件的涡流密度分布
S2.1、计算被测试件在点a处感生出的电场E1
Figure GDA0002227749950000053
其中,S为探头下一个圆形闭合平面的面积,上标箭头表示向量,在实际情况下,由于该面积很小,可以将S平面内通过的磁感应强度看作均匀电场。
S2.2、计算点a处被测试件表面的涡流密度J;
Figure GDA0002227749950000054
其中,σ为被测试件的电导率,M为常数;
当激励线圈在正弦信号I1=Acos(ωt)的激励下,那么在点a处被测试件表面的涡流密度J表示为:
J=σMSAω·sin(ωt)
其中,A为正弦信号的幅度,ω为正弦信号的角频率;
S2.3、在无限导体平面内,导体内的电流密度会从表面到体内逐渐衰减,那么被测试件内部的涡流密度J(h');
J(h')=J·eh'=σMSAω·sin(ωt)·eh'
其中,J(h')表示被测试件内部h'深度处的涡流密度,δ表示在角频率ω下被测试件中涡流的趋肤深度;
S3、由于探头相对于被测试件的尺寸很小,因此被测试件可以被看成是无限大导体表面,且涡流在导体中呈环状,同环形电流理论一样,涡流产生的磁感应强度B2可表示:
Figure GDA0002227749950000061
其中,B2方向与B1相反,μ2为被测试件的磁导率,rs为被测试件中涡流的等效最大半径,J(z)表示试件z深度处的涡流密度;
在本实施例中,由于探头尺寸很小,因此等效涡流半径rs在这里认为是一个很小的常数,ω=2πf,f为该激励下的频率,则B2可表示为:
Figure GDA0002227749950000062
考虑到涡流的衰减和传感器的尺寸,涡流的有效半径rs远大于d+z,为简化计算,令
Figure GDA0002227749950000067
则经积分计算可得:
Figure GDA0002227749950000063
通过上式,可以计算出在同一激励频率下,从导体表面到h深度的涡流在点a处产生的磁感应强度。
S4、建立稳态特性与被测试件的厚度关系
S4.1、估算被测试件在完全被磁化时的磁导率;
被测试件在完全被磁化时满足:μH>εH0s0m0h,其中,h为被测试件厚度,s0为单位体积的底面积,ε为磁极子间摩擦衰减因子,H为稳态磁场强度,H0为单个磁极子提供的磁场强度,m0为单位体积可磁化磁极子数;
稳态特性下,稳态磁场强度与单磁极子的关系为:H=kH0,其中,k为磁场合成因子;
那么被测试件在完全被磁化时的磁导率估计值
Figure GDA0002227749950000064
为:
Figure GDA0002227749950000065
S4.2、用
Figure GDA0002227749950000066
替换μ2,估算出被测试件内部的涡流磁感应强度B2
Figure GDA0002227749950000071
S4.3、估算被测试件的厚度h;
Figure GDA0002227749950000072
在本实施例中,当被测试件的材料为颗粒状高导磁材料时,h远大于δ,那么指数的负无穷就趋向零,则e-h/δ=0,此时
Figure GDA0002227749950000073
上式等号右边中除了厚度h,其他都是可测或已知常数,因此厚度h和涡流中产生的磁感应强度B2存在一一对应的关系,可用此式测得B2进而求取材料厚度h值。经变换可得:
Figure GDA0002227749950000074
S5、计算待测试件的厚度
S5.1、当只有激励线圈工作时,点a处未放置被测试件,利用HALL传感器测出激励线圈在点a处的磁感应强度B1
S5.2、当磁场稳定时,在点a处放置被测试件,然后利用HALL传感器测出点a处的磁感应强度,记为B;B是激励线圈磁感应强度B1与涡流磁感应强度B2叠加而成,且B2与B1方向相反,故B2=B-B1
S5.3、将测出的涡流磁感应强度B2代入至步骤S4.3的公式中,计算出被测试件的厚度。
在本实施例中,对厚度已知的材料试件根据上述方法进行测试,记录下实际值与测试值,并计算其误差。测试数据如表1所示,测试值与实际值的曲线图如图3所示,从图中可以看出,此方法的误差很小,在理想的测量范围之内。
Figure GDA0002227749950000081
表1
对不同厚度的材料试件进行多次测试,记录下测试值,求厚度的均值和标准差,并对均值进行拟合。测试数据如表2所示,,可以看出,标准差是比较小的,说明数据的波动不大,也就说这种测厚方法的稳定性很好。均值的拟合曲线如图4所示,可以看出测试的厚度均值与涡流磁感应强度有着良好的线性关系,这也方便了厚度测量。
实际厚度 112 220 336 438 497 598
测量点1 123 210 354 443 513 586
测量点2 113 226 338 418 517 590
测量点3 108 213 362 436 486 613
测量点4 100 206 336 426 509 602
测量点5 111 211 343 425 516 594
平均值 111 213.2 346.6 429.6 500.2 597
标准差 8.33667 7.5961 11.082 9.8641 12.795 10.72
表2
在本实施例中,通过上述方法和操作,从得到的结果来看,用本发明进行材料测厚时,误差比较小,准确度很高,材料厚度均值与涡流磁感应强度有着良好的线性关系,操作很方便;因此,本发明完全可以应用于实际测厚中。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。

Claims (1)

1.一种基于涡流稳态特性的磁性金属材料厚度检测方法,其特征在于,包括以下步骤:
(1)、建立测试模型
在激励线圈中心插入圆柱型探针,在探针垂线方向,距离探针下端探头距离为d处设为点a,在点a处放置被测试件,在探针下端探头处设置HALL传感器;
在点a处由激励线圈产生的磁感应强度B1为:
Figure FDA0002227749940000011
其中,μ为探针的磁导率,I1为激励线圈的激励电流,r为激励线圈的半径,rp为探头的半径,n为激励线圈的匝数,
Figure FDA0002227749940000012
为激励线圈的直径;
(2)、计算被测试件的涡流密度分布
(2.1)、计算被测试件在点a处感生出的电场E1
Figure FDA0002227749940000013
其中,S为探头下一个圆形闭合平面的面积;
(2.2)、计算点a处被测试件表面的涡流密度J;
Figure FDA0002227749940000014
其中,σ为被测试件的电导率,M为常数;
当激励线圈在正弦信号I1=Acos(ωt)的激励下,那么在点a处被测试件表面的涡流密度J表示为:
J=σMSAω·sin(ωt)
其中,A为正弦信号的幅度,ω为正弦信号的角频率;
(2.3)、计算被测试件内部的涡流密度J(h');
J(h')=J·eh'/δ=σMSAω·sin(ωt)·eh'/δ
其中,J(h')表示被测试件内部h'深度处的涡流密度,δ表示在角频率ω下被测试件中涡流的趋肤深度;
(3)、根据被测试件内部的涡流密度J(h')计算被测试件内部的涡流磁感应强度B2
Figure FDA0002227749940000021
其中,B2方向与B1相反,μ2为被测试件的磁导率,rs为被测试件中涡流的等效最大半径,J(z)表示试件z深度处的涡流密度;
(4)、建立稳态特性与被测试件的厚度关系
(4.1)、估算被测试件在完全被磁化时的磁导率;
被测试件在完全被磁化时满足:μH>εH0s0m0h,其中,h为被测试件厚度,s0为单位体积的底面积,ε为磁极子间摩擦衰减因子,H为稳态磁场强度,H0为单个磁极子提供的磁场强度,m0为单位体积可磁化磁极子数;
稳态特性下,稳态磁场强度与单磁极子的关系为:H=kH0,其中,k为磁场合成因子;
那么被测试件在完全被磁化时的磁导率估计值
Figure FDA0002227749940000022
为:
Figure FDA0002227749940000023
(4.2)、用
Figure FDA0002227749940000024
替换μ2,估算出被测试件内部的涡流磁感应强度B2
Figure FDA0002227749940000025
(4.3)、估算被测试件的厚度h;
Figure FDA0002227749940000026
(5)、计算待测试件的厚度
(5.1)、当只有激励线圈工作时,点a处未放置被测试件,利用HALL传感器测出激励线圈在点a处的磁感应强度B1
(5.2)、当磁场稳定时,在点a处放置被测试件,然后利用HALL传感器测出点a处的磁感应强度,记为B;B是激励线圈磁感应强度B1与涡流磁感应强度B2叠加而成,且B2与B1方向相反,故B2=B-B1
(5.3)、将测出的涡流磁感应强度B2代入至步骤(4.3)的公式中,计算出被测试件的厚度。
CN201910129231.5A 2019-02-21 2019-02-21 一种基于涡流稳态特性的磁性金属材料厚度检测方法 Active CN109668506B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910129231.5A CN109668506B (zh) 2019-02-21 2019-02-21 一种基于涡流稳态特性的磁性金属材料厚度检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910129231.5A CN109668506B (zh) 2019-02-21 2019-02-21 一种基于涡流稳态特性的磁性金属材料厚度检测方法

Publications (2)

Publication Number Publication Date
CN109668506A CN109668506A (zh) 2019-04-23
CN109668506B true CN109668506B (zh) 2020-03-24

Family

ID=66152092

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910129231.5A Active CN109668506B (zh) 2019-02-21 2019-02-21 一种基于涡流稳态特性的磁性金属材料厚度检测方法

Country Status (1)

Country Link
CN (1) CN109668506B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111999378B (zh) * 2020-08-18 2023-08-29 中国人民解放军空军工程大学 一种基于tmr传感器的金属材料电导率和厚度的测量方法
CN112179261B (zh) * 2020-09-24 2022-09-23 桂林理工大学 一种基于电磁响应的钢轨脱碳层厚度的检测方法
CN112432588B (zh) * 2020-12-08 2021-09-14 电子科技大学 一种通过磁饱和特性测量吸波涂层厚度的方法
CN112666395B (zh) * 2020-12-17 2024-05-31 无锡学院 非接触式的金属材质电导率测量方法及系统
CN116576807B (zh) * 2023-07-11 2023-09-05 中国石油大学(华东) 一种无线能量和信号传输的电磁超声体波测厚装置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105300267B (zh) * 2015-11-17 2017-10-31 电子科技大学 基于电磁涡流和数据融合算法的金属板厚度测量方法
DE102015120067B4 (de) * 2015-11-19 2018-08-16 Sick Ag Verfahren zur Erkennung eines Profils einer einen Detektionsbereich eines Abstandssensors durchlaufenden Anordnung
JP6780538B2 (ja) * 2017-02-23 2020-11-04 日本製鉄株式会社 形状測定装置及び形状測定方法
CN107990820B (zh) * 2017-11-28 2019-09-24 四川元匠科技有限公司 一种基于脉冲涡流的金属板厚度信息检测方法

Also Published As

Publication number Publication date
CN109668506A (zh) 2019-04-23

Similar Documents

Publication Publication Date Title
CN109668506B (zh) 一种基于涡流稳态特性的磁性金属材料厚度检测方法
Li et al. A thickness measurement system for metal films based on eddy-current method with phase detection
CN101131314B (zh) 一种铁基体上镍镀层的无损测厚方法
Zhao et al. Uniaxial ACFM detection system for metal crack size estimation using magnetic signature waveform analysis
CN103499404A (zh) 铁磁构件交变应力测量装置及其测量方法
Wang et al. Design of ultrastable and high resolution eddy-current displacement sensor system
CN104316594B (zh) 一种钢件缺陷的电磁无损检测装置
CN107990820A (zh) 一种基于脉冲涡流的金属板厚度信息检测方法
CN105300266B (zh) 一种汽车用镀锌板镀层厚度的电涡流检测装置与方法
CN106524892A (zh) 一种基于涡流磁导率测量的钢管壁厚测量方法
CN107677888A (zh) 一种金属电导率涡流检测方法
CN102175131B (zh) 一种利用漏磁场测量钢板厚度的方法
CN108267700A (zh) 一种永磁体抗退磁能力的检测系统及其检测方法
CN104483530B (zh) 高温超导块材表层临界电流密度不均匀性测量方法及装置
Krause et al. A pulsed eddy current probe for inspection of support plates from within alloy-800 steam generator tubes
CN113607214B (zh) 一种金属管道参数确定方法及系统
CN114764086B (zh) 基于偏置磁化下涡流检测差动磁导率的管道内检测方法
Ge et al. Development of a velocity-adaptable alternating current field measurement device for crack inspection in rails
CN103940902B (zh) 利用涡流阻抗平面检测仪检测非金属材料不连续性方法
Cheng et al. Absorbing coating thickness measurement based on lift-off effect of eddy current testing
Aguila-Munoz et al. Crack detection in steel using a GMR-based MFL probe with radial magnetization
Long et al. A characteristic approximation approach to defect edge detection in magnetic flux leakage testing
Yang et al. Excitation structure design and magnetic field analysis of a new electromagnetic flowmeter based on magnetically permeable material
Gotoh et al. Measurement of thickness of nickel-layer on steel using electromagnetic method
Wang et al. Inductance calculation of planar eddy-current sensor coils in grating-type displacement measurement system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant