基于超声测量模型的镁合金弹壳自动检测方法

刘希玲 ++吕铎 ++李雄兵 ++倪培君

摘要：镁合金弹壳是具有变厚度变曲率特征的组合回转体，根据其结构特点提出了一种镁合金弹壳超声自动检测方案，设计了相应的超声自动检测系统。对六自由度超声自动检测台进行了运动学建模，实现了超声扫描运动路径与探头位姿的精确控制。以壳体内缺陷散射模型和多元高斯声束模型为基础，建立了壳体缺陷超声测量模型，并绘制了壳体缺陷定量表征曲线。利用六自由度超声自动检测台对镁合金弹壳试块进行超声扫描实验，实现了对壳体人工缺陷的检出与定量，从而验证了本文方法的有效性。

关键词：超声；镁合金；弹壳；自动检测；缺陷定量

中图分类号：TB553 文献标识码：A

新型镁合金导弹弹壳相比传统钢质或铝合金弹壳具有明显的减重优势，其应用是实现武器装备轻量化的重要技术途径＼[1＼]。受原材料及加工工艺的影响，镁合金弹壳易出现夹杂、裂纹等缺陷，如不能及时有效检出上述缺陷，将严重影响武器装备的性能和安全。超声无损检测具有检测灵敏度高、检测深度大、缺陷定位准确、对人体无害等优点＼[2＼]，是进行镁合金弹壳缺陷检测的较理想方法。弹壳是具有变曲率变厚度特征的组合回转体，受曲面的影响，超声波在壳体中传播时的聚焦、散射特性增强＼[3＼]，不同部位的壳体厚度也将产生不同的声压衰减＼[4＼]，因此缺陷的检出与定量难度较高。本文利用六自由度超声自动检测台实现超声扫描运动及探头位姿的精确控制，将超声测量模型应用于镁合金弹壳缺陷检测，建立了缺陷定量表征曲线，考虑了壳体厚度、表面曲率对超声响应的影响，通过检测镁合金弹壳试块中的平底孔人工缺陷，验证了本文方法的有效性。

镁合金弹壳超声自动检测技术方案见图2，可分为2个步骤。步骤1：对六自由度检测台进行运动学建模，得出运动学方程正解与逆解，结合镁合金弹壳的CAD模型控制壳体与超声探头的预定路径相对运动，实现全壳体完整超声扫描，从而获取所需的超声A波信息。步骤2：针对镁合金材质特性建立壳体内部多元高斯声束模型，结合由基尔霍夫近似得到的缺陷散射模型，建立壳体缺陷超声测量模型，基于该模型建立壳体缺陷定量表征曲线。最终通过扫描获得的超声A波波形特征来获取缺陷信息，实现缺陷的检出；通过获取相对缺陷波高，利用超声测量模型反求出理论缺陷半径值，实现缺陷的定量。

2检测过程的运动控制

镁合金弹壳一般为组合回转体结构，回转体表面为平面、锥面或者复杂曲面，弹壳在不同部分有不同的曲率和厚度。在对镁合金弹壳进行自动检测时，随着弹壳表面曲率和厚度的变化，信号中表面波、缺陷波、底波的位置和幅度变化很大，要想获得准确有效的超声信号，以便于后续信号处理，需要利用多自由度机械手精确控制超声探头的位姿。

采用CYS100型六自由度超声自动检测台来实现超声探头位姿与扫描运动的精确控制。检测台的结构和运动原理见图3，检测台由一个五自由度空间开式运动链（五轴机械手）和一个独立旋转轴（底部转盘）组合，开式运动链末端装有超声探头。检测台的五轴机械手包含3个移动关节（A1，A2，A3）和2个转动关节（A4，A5），5个关节变量分别为x，y，z，θ1，θ2。转动关节A5对应的连杆长度为l1，探头的夹持长度为l2。检测台的底部转盘包含一个转动关节A6，关节变量为θ3。

根据DenavitHartenberg方法＼[5-6＼]建立五轴机械手的运动学方程，运动方程正解见式（1），运动方程逆解见式（2）。

αx=-sinθ2，

αy=-sinθ1cosθ2，

αz=cosθ1cosθ2，

px=l2sinθ2+x，

py=l2sinθcosθ2-l1cosθ1+l1+y，

pz=-l2cosθ1cosθ2-l1sinθ1+l2+z；（1）

θ1=tan-1（-αyα-1z），

θ2=tan-1（αxα-1ysinθ1），

x=px-l2sinθ2，

y=py-l2sinθ1cosθ2+l1cosθ1-l1，

z=pz+l2cosθ1cosθ1+l1sinθ1-l2。（2）

式中：（αx，αy，αz）为检测点表面法向矢量；（px，py，pz）为检测点坐标值。

图3检测台结构和运动原理图

Fig。3Schematic diagram of testing platform

and its movement principle

3镁合金弹壳缺陷定量检测

3。1超声测量模型的建立

以镁合金弹壳圆锥面部分为例，介绍超声测量模型的建立过程。如图4所示，ρ1，ρ2分别为水和镁合金弹壳材质的密度；c1，c2分别为水和镁合金弹壳中的声速；a，b分别为超声探头端面和平底孔缺陷的半径；s1为水声距；s2为缺陷深度。

镁合金弹壳材质内部各向均一同性，超声传播过程中遇到平底孔缺陷时形成散射回波，其幅值为＼[7＼]：

A（ω）=ik2b22。（3）

式中：k2为超声波在镁合金材质中的波数。超声检测系统获取的平底孔缺陷回波电压频域响应为＼[8＼]：

VR（ω）=s（ω）A（ω）v（ω）24πρ2c2-ik2Z。（4）

式中：s（ω）为脉冲激励系统函数，可通过超声脉冲激励实验获取＼[9＼]；v（ω）为探头辐射声场中质点振动相对速度；Z=ρ1c1S为探头辐射阻抗，S为探头端面面积。

图4声束传播示意图

Fig。4Schematic diagram of ultrasonic transmission

当超声垂直入射，且仅考虑纵波声场时，利用多元高斯模型，探头辐射声场中质点振动速度为＼[10＼]：

V（ω）d=∑15r=1[v0（0）]rdet[M2（s2）]rdet[M2（0）]rdet[Μ1（s1）]rdet[Μ1（0）]r×

ArT12dexpik1s1+ik2s2+ik12XTc1Μ2（s2）rX。（5）

式中：v0（0）为探头表面处质点振动速度；Ar为高斯系数＼[11＼]；T12为声压透射系数；d为单位向量；k1为超声波在水中的波数；X=[x，z]为探头辐射声场中振动质点位置坐标；M1（0）和M1（s1）分别为超声波在水中起始和传播相位矩阵，M2（0）和M2（s2）分别为超声波在试块中起始和传播相位矩阵。利用Huang提出的ABCD转换矩阵法＼[12＼]，可得

M1（0）=iBr/c1D00iBr/c1D； （6）

M1（s1）=1/c1（s1-iD/Br）001/c1（s1-iD/Br）；（7）

M2（0）=M/c1（s1-iD/Br）00M/c1（s1-iD/Br）；（8）

M2（s2）=1c1（s1-D/Br）/M+s2c2001c1（s1-D/Br）+s2c2； （9）

M=1+（s1-iD/Br）（1-c1/c2）/R。（10）

式中：Br为高斯系数；D=k1a2/2为瑞利距离；R为超声入射界面曲率半径。探头辐射声场中质点振动相对速度为：

v（ω）=V（ω）/v0（0）。（11）

超声检测系统获取的底面回波电压频域响应为：

VR（ω）=s（ω）t（ω）。（12）

式中：t（ω）为声/弹性转化函数。最终超声测量模型表达式为：

VR（ω）=2s（ω）1S∫Sv（ω）dS（X）。（13）

将式（3）和式（11）及脉冲激励系统函数s（ω）代入式（4）中，可得到平底孔缺陷回波电压频域响应，将式（11）及脉冲激励系统函数s（ω）代入式（13）中，可得到底面回波电压频域响应，对频域响应进行傅里叶逆变换可得到相应回波信号时域响应。

3。2缺陷定量表征曲线的绘制

在缺陷定量检测中，通常用缺陷回波电压与底面回波电压的相对幅值来判定缺陷当量大小。缺陷定量表征曲线描述了缺陷/底面回波电压相对幅值随超声入射界面曲率、缺陷深度、缺陷尺寸等参数的变化规律。为对镁合金弹壳进行缺陷定量，可分两步绘制平底孔缺陷定量表征曲线：1）利用所建立的超声测量模型计算底面回波与缺陷回波的电压幅值，并以底面回波电压幅值为归一化基准，计算电压相对幅值；2）改变超声测量模型中超声入射界面曲率、缺陷深度、缺陷尺寸等参数，通过拟合运算可得到缺陷定量表征曲线。

图5（a）为利用超声测量模型绘制的相同曲率，不同缺陷半径及缺陷深度的缺陷定量表征曲线。图中横坐标为缺陷深度，纵坐标为缺陷/底面回波电压相对幅值。模型中试块材质为镁合金，设定试块厚度为50 mm，缺陷深度为1～45 mm，平底孔缺陷半径b依次为1。0，1。5，2。0，2。5和3。0 mm，超声入射界面曲率半径R为70 mm。超声探头频率为2。25 MHz，直径为12。7 mm，水声距为75 mm。

图5（b）为在相同条件下利用超声测量模型绘制的相同缺陷半径，不同曲率、不同缺陷深度的缺陷定量表征曲线。平底孔缺陷半径均为1 mm，R为超声入射界面曲率半径。

综上所述，利用缺陷定量表征曲线进行缺陷定量的步骤如下：1）利用所建立的超声测量模型，绘制缺陷定量表征曲线；2）获取实验测试底面回波与缺陷回波幅值，计算缺陷/底面回波电压相对幅值；3）根据缺陷回波和表面回波的时间差，确定缺陷深度；4）对比缺陷/底面回波电压相对幅值与模型预测结果。如果实验与预测结果吻合良好，其误差在允许范围之内，则定量表征曲线对应缺陷半径即为实际检测缺陷当量值。

4检测实例

针对图6（a）所示的镁合金弹壳进行缺陷检测，在该试块的壳体内表面利用电火花法加工了5个平底孔缺陷，缺陷所在位置的壳体外圆半径为70 mm，1#～5#缺陷半径依次为1，2，3，3，3 mm；缺陷深度依次为15，15，15，10，5 mm。采用OLYMPUS5072P/R脉冲发射接收仪，直径12。7 mm，频率2。25 MHz的超声探头，以脉冲激励方式进行实验，增益-3 dB，检测水声距为75 mm，利用CYS1100型六自由度超声自动检测台精确控制扫描运动。

图6（b）为实验测试得到的平底孔缺陷/底面回波电压相对幅值与定量表征曲线预测结果的对比，详细实验测试数据以及模型预测数据如表1所示。

结合图6与表1的数据分析可知，模型预测的相对缺陷波高值与实验测试结果的绝对误差小于1。1 dB，且利用超声测量模型判定缺陷半径的当量尺寸，与实际缺陷半径值之间的相对误差不超过6。5%。可见，基于超声测量模型的缺陷定量表征曲线，能有效辅助缺陷定量。

（a）镁合金弹壳人工缺陷示意图

s2/mm（b）实验测试与定量表征曲线对比

图6检测实例

Fig。6Test example

表1预测及实验数据

Tab。1Prediction and experimental data

缺陷

序号

缺陷

深度

/mm

预测

声压

/dB

实验

声压

/dB

声压

误差

/dB

预测缺

陷半径

/mm

实际缺

陷半径

/mm

半径相

对误差

/%

1#

15

-27。9

-29。0

1。1

1。06

1

6。0

2#

15

-21。4

-22。5

1。1

2。13

2

6。5

3#

15

-16。8

-17。5

0。7

3。16

3

5。3

4#

10

-16。0

-15。3

-0。7

2。85

3

5。0

5#

5

-16。4

-16。7

0。3

3。12

3

4。0

实验测试数据与模型预测数据的差异主要是建立超声测量模型以及绘制缺陷定量曲线时对超声传播过程的理想假设，此外还有超声检测系统带来的误差。

5结论

1）对六自由度超声自动检测台进行了运动学建模，实现了超声扫描运动与探头位姿的精确控制，满足镁合金弹壳等回转体的超声检测要求。

2）结合多元高斯声束模型和缺陷散射模型，建立了缺陷超声测量模型，并绘制了缺陷定量表征曲线。根据超声测量模型进行反求计算所得的缺陷半径理论值与实际值误差不超过6。5%。本文方法为实现镁合金弹壳缺陷的检出与定量提供了一种新的有效途径。

参考文献

[1]肖冰，康凤，胡传凯，等。 国外轻质结构材料在国防工业中的应用＼[J＼]。 兵器材料科学与工程， 2011， 34（1）： 94-97。

XIAO Bing， KANG Feng， HU Chuankai， et al。 Application of the light structure material in the defense industry abroad＼[J＼]。Ordnance Material Science and Engineering， 2011， 34（1）： 94-97。（In Chinese）

[2]李喜孟。 无损检测＼[M＼]。 北京：机械工业出版社，2011：18-30。

LI Ximeng。 Nondestructive testing＼[M＼]。 Beijing： China Machine Press， 2011：18-30。 （In Chinese）

[3]LEE J R， TAKATSUBO J J， TOYAMA N， et al。 Health monitoring of complex curved structures using an ultrasonic wave field propagation imaging system＼[J＼]。 Measurement Science and Technology， 2007， 18（12）： 3816-3825。

[4]张杨，周晓军，杨辰龙，等。 变厚度曲面构件超声检测灵敏度补偿＼[J＼]。 浙江大学学报：工学版，2013， 47（1）： 116-121。

ZHANG Yang， ZHOU Xiaojun， YANG Chenlong， et al。 Sensitivity compensation method in ultrasonic inspection curved components with variable thickness＼[J＼]。 Journal of Zhejiang University：Engineering Science， 2013， 47（1）：116-121。 （In Chinese）

[5]孙树栋。 工业机器人技术基础＼[M＼]。 西安：西北工业大学出版社，2006：31-47。

SUN Shudong。 Fundamentals of robotics＼[M＼]。 Xi′an： Northwestern Polytechnical University Press， 2006：31-47。 （In Chinese）

[6]谢宏，王朝辉，邹帆，等。 基于遗传Hopfield混合神经网络的多轴控制系统设计与实现＼[J＼]。 湖南大学学报：自然科学版，2012，39（7）：44-48。

XIE Hong， WANG Zhaohui， ZOU Fan， et al。 Design and realization of multiaxis control system based on modified GAHopfiled neutral network＼[J＼]。 Journal of Hunan University：Natural Sciences， 2012，39（7）：44-48。 （In Chinese）

[7]KIM H J， SCHMERR L W， SEDOV A。 Transferring distanceamplitude correction curves using ultrasonic modeling＼[C＼]//Review of Quantitative Nondestructive Evaluation。Green Bay， Wisconsin，2004：753-756。

[8]THOMPSON R B， GRAY T A。A model relating ultrasonic scattering measurements through liquidsolid interfaces to unbounded medium scattering amplitudes＼[J＼]。 Journal of the Acoustical Society of America， 1983， 74（4）： 1279-1290。

[9]DANG C J。Electromechanical characterization of ultrasonic NDE systems＼[D＼]。 Ames：Center for Nondestructive Evaluation，Iowa State University，2001。

[10]SCHMERR L W，SONG S J。Ultrasonic nondestructive evaluation systems models and measurements＼[M＼]。 New York： Springer Press，2007：152-167。

[11]WEN J J，BREAZEALE M A。A diffraction beam field expressed as the superposition of gaussian beams＼[J＼]。 Journal of the Acoustical Society of America，1988， 83（5）：1752-1756。

[12]HUANG R。 Ultrasonic modeling for complex geometries and materials＼[D＼]。 Ames： Center for Nondestructive Evaluation，Iowa State University，2006。