电磁超声无损检测技术的ANSYS仿真研究

电子测量技术

ＥＬＥＣＴＲｏＮＩＣＭＥＡＳＵＲＥＭ［ＥＮＴＴＥＣＨＮｏＬＯＧＹ第３１卷第７期２００８年７月

电磁超声无损检测技术的ＡＮＳＹＳ仿真研究

任晚可李健

（天津大学精密测试技术及仪器国家重点实验室天津３０００７２）

摘要：利用有限元分析软件ＡＮＳＹＳ对电磁超声无损检测技术进行了研究，通过仿真，比较了电磁超声检测过程中缺陷对被测导体表面磁流密度、磁感应强度和涡流等的影响，并介绍了ＡＮＳＹＳ软件进行瞬态电磁场仿真的过程，为电磁超声的研究提供了帮助。

关键词：电磁超声；有限元；无损检测；管道

中图分类号：ＴＨ８７８文献标识码：Ａ

ＳｉｍｕｌａｔｉｏｎｒｅｓｅａｒｃｈｏｎｅｌｅｃｔｒｏｍａｇｎｅｔｉｃａｃｏｕｓｔｉｃＮＤＴｂｙＡＮＳＹＳ

ＲｅｎＸｉａｏｋｅＬｉＪｉａｎ

（ＳｔａｔｅＫｅｙＬａｂｏｒａｔｏｒｙｏｆＰｒｅｃｉｓｉｏｎＭｅａｓｕｒｉｎｇＴｅｃｈｎｏｌｏｇｙａｎｄＩｎｓｔｒｕｍｅｎｔｓ，ＴｉａｎｊｉｎＵｎｉｖｅｒｓｉｔｙ，Ｔｉａｎｊｉｎ３０００７２）

Ａｂｓｔｒａｃｔ：Ｔｈｉｓｐａｐｅｒｒｅｓｅａｒｃｈｅｓｔｈｅｔｅｃｈｎｏｌｏｆｙｏｆｎｏｎｄｅｓｔｒｕｃｔｉｖｅｔｅｓｔｉｎｇ（ＮＤＴ）ｉｎｅｌｅｃｔｒｏｍａｇｎｅｔｉｃａｃｏｕｓｔｉｃｔｒａｎｓｄｕｃｅｒ（ＥＭＡＴ）ｕｓｉｎｇｆｉｎｉｔｅｅｌｅｍｅｎｔｓｏｆｔｗａｒｅＡＮＳＹＳ．Ｉｔｓｉｍｕｌａｔｅｓｔｈｅｉｎｆｌｕｅｎｃｅｗｈｉｃｈｃｏｍｅｓｆｒｏｍｄｅｆｅｃｔｉｎｍａｇｎｅｔｉｃｆｌｕｘｄｅｎｓｉｔｙ，ｍａｇｎｅｔｉｃｉｎｄｕｃｔｉｏｎｉｎｔｅｎｓｉｔｙａｎｄｅｄｄｙｃｕｒｒｅｎｔｄｕｒｉｎｇｔｈｅｅｌｅｃｔｒｏｍａｇｎｅｔｉｓｍｕｌｔｒａｓｏｕｎｄｔｅｓｔｉｎｇ．ＴｈｅｓｔｅｐｓｏｆｕｓｉｎｇＡＮＳＹＳｔＯｓｉｍｕｌａｔｅｔｒａｎｓｉｅｎｔｓｔａｔｅｏｆｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｆｉｅｌｄａｒｅｉｎｔｒｏｄｕｃｅｄ，ｗｈｉｃｈｃｏｕｌｄｈｅｌｐｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｕｌｔｒａｓｏｎｉｃｔｅｓｔｉｎｇｔｅｃｈｎｏｌｏｇｙｒｅｓｅａｒｃｈｄｉｒｅｃｔｌｙ．

Ｋｅｙｗｏｒｄｓ：ＥＭＡＴ；ｆｉｎｉｔｅｅｌｅｍｅｎｔｍｅｔｈｏｄ；ｎｏｎ－ｄｅｓｔｒｕｃｔｉｖｅｔｅｓｔｉｎｇ；ｐｉｐｅｌｉｎｅ

０引言

管道作为能源运输的主要方式，其安全性是相当重要的。油气管道无损检测技术就是预测隐患、为安全评价提供基础资料的一种重要手段。管道无损检测技术的原理有很多种，目前应用较多的是超声法、漏磁法和涡流法。新兴的电磁超声技术作为一种激发原理特殊的超声，因其不需要耦合介质，易于选定所激发的超声模式等优点，正引起无损检测领域越来越高的重视［１］。

电磁超声虽然具有众多优点，但由于其转换效率低，衰减严重，受噪声污染大，检测结果常常不理想。因此，利用有限元分析软件，建立ＥＭＡＴ的二维模型，通过仿真比较了缺陷对被测物体表面磁流密度、磁感应强度和涡流等情况的影响，提出了设计ＥＭＡＴ检测装置和进行检测数据分析时值得注意的方面，对指导ＥＭＡＴ无损检测的研究具有重要的意义。

１电磁超声检测技术的理论依据

在被测物体表面设置一曲折线圈，某一时刻，通过往线圈中加载一高压窄脉冲或时谐电流源激励产生一个交变的电磁场Ｈ，会在被测物体内部产生感应电流，即涡流Ｉ。如果同时施加一由永磁体产生的稳定磁场曰，就会与金属内部涡流相互作用，产生洛伦兹力。被测物体的质点在此力的作用下产生机械振动，从而产生沿被测件辐射或沿表面传播的超声波。这就是电磁超声激发的基本原理［２］。电磁超声的接收是激发的逆过程。当被测物体表面有超声自内部投射时，质点发生位移，带正电荷的晶格在偏置磁场的作用下受力，产生交变电流。这个交变电流将导致被测导体的表层出现交变的磁场，这个交变磁场漏出导电体，在被测导体上方的线圈中感生出电动势。这样，就可以进行检测了。在电磁超声无损检测中，被测物体作为龟磁超声传感器的一部分，必须是电导体或磁导体。若被测物体是铁磁性材料，除洛仑兹力外，还受到磁致伸缩力的作用‘３ｌ。

２仿真方法

有限元分析是将物体划分成有限个单元，单元看作是不可变形的刚体，通过有限个节点相互连接，单元之间的力通过节点传递，然后利用能量原理建立各单元矩阵，在输入材料特性、载荷和约束等边界条件后，利用计算机进行物体特性计算和分析。ＡＮＳＹＳ软件是目前最常用的有限元分析软件之一，融结构、流体、电场、磁场、声场分析于一体‘“。软件主要包括３个模块：前处理模块，求解模块和后处理模块。

任晓可等：电磁超声无损检测技术的ＡＮＳＹＳ仿真研究

第７期

２．１前处理

根据电磁超声无损检测装置的结构和工作特点，取单元类型为空气、线圈和被测导体，定义材料属性。设被测导体的相对磁导率胁一５００，电阻率Ｐ＝９．１７×１０－８０·ｒｎ，线圈的相对磁导率雎一１，电阻率Ｐ＝１．６８×１０１

Ｑ·ｍ。

考虑到模型具有对称性，且假设线圈导线的长度比线间距要大得多，可将模型简化成二维情况来进行有限元分析口８］。设被测物体为半无限大导体，在其上方放置曲折线圈，线圈由一根导线来回沿直角弯曲成一组６条平行线，相邻导线间距ｂ一１ｍｍ，提离距离ｈ一０．５ｍｍ，导线横截面为正方形，边长ａ一０．２１ＴＬｒｎ，相邻两导线电流方向相反。导体处在磁感应强度Ｂ。＝１Ｔ的偏置磁场中，方向竖直向下。为比较缺陷对电磁超声检测过程中各因素的影响，在导线３下方设置了一个０．１

ｍｍ×０．２ｌＴＬｒｎ的裂缝，ＥＭＡＴ

模型如图１所示。建立模型后，给面赋予特性，确定单元

关系，并进行智能网格划分。

图１

ＥＭＡＴ激发系统二维平面模型

２．２求解

在电磁场微分方程的求解过程中，只有限制了边界条件和初始条件，才有确定解，因此，为被测导体设置边界条件，并规定各项力标志。根据电磁超声的激发机理，采用时变电流对线圈施加激励。设激励电流的形式为ｉ（ｆ）一

１００ｓｉｎｚ（７ｃ×１０５ｔ）ｃｏｓ（ｎ×１０６ｔ），取分析时间为１０脚，此时激励电流可视为衰减到０，时间步长设为０．２卢ｓ，信号波形如图２所示。根据所加激励的类型，耦合线圈自由度，进行瞬态磁场分析。

蠢

时间，。

图２激励电流信号波形

２．３后处理

求解结束后，可选择相应时间步，获得分析过程中相关的等值线图，云图和矢量图等。可通过选择节点，读取其在分析过程中某一参数随时间或频率连续变换的情况啪。还可以以麦克斯韦方程组为出发点，求解磁通、力和能量等参数。

３仿真实验结果

当线圈受时谐电流源激励时，会产生一个交变的电磁

场Ｈ，交变磁场的磁流密度如图３和图４所示。磁流密度的具体数值可通过等值线读取。可见当不存在缺陷时，各条导线下方的磁流密度大致相等，当存在缺陷时，缺陷两侧钢板表面的磁流密度就比其它地方大，而缺陷深处磁流密度却变得非常小。

彳／’＼一』

一生／＼＼

／⑨⑩◎⑩⑩⑥＼

Ｄ

Ｅ

卫

。』Ｄ一

￡一Ｃ

ＭＸ

Ｅ

妥滔滞司◇

，—气

，、

一

．一

，、＞＜

Ａ＝．６０５Ｅ埘

Ｂ＝．１８２Ｅ．０３Ｃ＝．３０３Ｅ－０３Ｄ＝．４２４Ｅ．０３Ｅ＝．５４５Ｅ－０３Ｆ＝．６６６Ｅ－０３Ｇ＝．７８７Ｅ．０３

日＝．９０８Ｅ．０３

Ｉ＝．００１０２９

图３被测导体表面不存在缺陷时交变磁场的磁流密度

衙画』画否一、、。＾Ｂ九Ｃ

ｄ．Ｋ

Ｄ，．、爿』、。．、／．

可一溺—～

图４被测导体表面存在缺陷时交变磁场的磁流密度

＝．８００Ｅ．０４＝．２４０Ｅ．０３＝．４００Ｅ．０３＝．５６０Ｅ．０３＝．７２０Ｅ．０３＝．８８０Ｅ．０３。．００１０４＝．００１２

一．００１３６

爿口ＣＤＥＦＧ

Ｈ，