针对电磁感应测速系统的有限元模拟

吴潇; 张乐; 马小娟; 张明建; 刘福生

doi:10.11858/gywlxb.20230693

针对电磁感应测速系统的有限元模拟

doi: 10.11858/gywlxb.20230693

西南交通大学物理科学与技术学院, 四川成都　610031

基金项目: 国家自然科学基金（41674008，11002120）；四川省自然科学基金（2022NSFSC0340）；西南交通大学专项（XJ2021KJZK055）

详细信息

作者简介:
吴　潇（1996－），男，博士研究生，主要从事动高压下的凝聚态物理研究.E-mail：2215189687@qq.com

通讯作者:
马小娟（1976－），女，博士，教授，博士生导师，主要从事冲击高压下的实验物理研究. E-mail：mxj_swjtu@126.com

中图分类号: O521.3; O44
计量
- 文章访问数: 165
- HTML全文浏览量: 88
- PDF下载量: 23
出版历程
- 收稿日期: 2023-07-12
- 修回日期: 2023-08-29
- 网络出版日期: 2023-12-06
- 刊出日期: 2023-12-15

Finite Element Simulation for Magnet Velocity Induction System

School of Physical Science and Technology, Southwest Jiaotong University, Chengdu 610031, Sichuan, China

摘要

摘要: 在冲击加载实验中，飞片速度是确定样品动态压缩性的重要参量。电磁感应测速系统是针对超高速运动状态的速度测量系统。基于有限元方法和ANSYS Electromagnetics Suite三维建模，模拟了飞片的测速过程，并验证了模拟结果的可靠性。结合实验结果，证实相同条件下模拟信号的感应电动势误差为1.1%，低于系统误差（2.0%）。模拟信号的测速误差小于0.4%，低于系统测速误差（0.9%）。分析了Al、Cu、Ta、W、93W飞片在不同飞片厚度、半径、倾斜角度、运动速度以及线圈孔径条件下检测线圈中感应电动势的变化规律。模拟结果为在冲击实验中更好地测量飞片速度提供了参考。
- 有限元方法 /
- 电磁感应测速系统 /
- 速度测量 /
- 三维模拟 /
- 感应电动势
Abstract: In shockwave experiments, the launch velocity of flyer is a crucial parameter for determining the physical quantities in sample under dynamic compression. Magnet velocity induction system is used to determine the velocity of flyer which can work under hypervelocity conditions. Based on finite element method, ANSYS Electromagnetics Suite module was used to establish a three-dimensional model and the launch process of flyer was simulated. The simulated results well reproduce the experimental signals. Comparing to the experimental results, the relative error of simulated amplitude of induced electromotive force (AIEMF) is 1.1%, lower than the system error of 2.0%. The relative error for flyer velocity determined by simulation is less than 0.4%, also lower than the measurement system error of 0.9%. The dependence of AIEMF with the thickness, radius, tilt angle, launch velocity of flyer and the diameter of pick-up coils is analyzed. The simulation result provides a reference for better signal acquisition in shockwave experiments.
- finite element method /
- magnetic velocity induction system /
- velocity measurement /
- three-dimensional model /
- amplitude of induced electromotive force

HTML全文

随着现代科技的高速发展，特别是机械工业^[1–2]、雷达工程^[3–7]、航空航天^[8–11]以及海洋工程^[12]等领域的快速发展，人们常常需要对物体的运动速度进行高精度测量，据此开展运动参数分析。冲击加载实验是研究物质在高温高压下的物理性质的实验方法^[13]。在轻气炮实验中^[14–18]，精确的速度测量是整个实验的关键。为了分析样品在高温高压下的冲击压强、粒子速度以及雨贡纽参数（C₀和λ）等重要物理量，需要准确测量飞片的速度^[13]。常用的测速方法有光电门测速法^[19]、气垫导轨测速法^[20]、冲击摆^[21]、多普勒效应测速法^[22–26]、电磁感应测速法^[27]。由于高速运动的物体前端存在等离子体发光现象，电磁感应测速法不仅能够克服等离子体发光对测速造成的影响，而且具有装置成本低、速度测量精度高等优点。

1851年，Foucault^[28]提出了涡旋电流的概念。1977年，Kondo等^[29]基于电磁感应原理发明了飞片速度测量技术。1984年，Moody等^[30]设计了基于涡流原理的电磁感应测速系统（magnet velocity induction system，MAVIS），利用检测线圈中的感应电动势（induced electromotive force，IEMF）信号计算出飞片的速度。1991年，施尚春等^[31]改进了MAVIS系统，在永磁环两端设置导磁片，增强了检测线圈中的IEMF信号幅度，提高了飞片速度的测量精度。2017年，张修路等^[32]利用COMSOL对MAVIS进行二维模拟，探究了飞片在不同速度下的IEMF信号波形特征。

MAVIS是一种基于涡流检测^[33–40]原理的速度测量装置。当飞片高速穿过MAVIS时，示波器能否被触发并记录完整的IEMF信号是测速实验成败的关键因素之一。因此，实验前，需对IEMF信号幅值及时间展宽进行预估，据此设置示波器的触发类型、触发电平、触发点位置、电压幅度和时基等。考虑到目前缺乏关于飞片厚度、半径、速度、发射倾斜角度以及线圈孔径对IEMF信号的影响的相关实验数据，因而有必要开展相应的模拟工作。为此，本研究基于有限元方法，利用ANSYS Electromagnetics Suite对MAVIS进行三维建模，分别模拟Al、Cu、Ta、W（金属钨）以及93W（UPGM-93WNIFE金属聚合物复合材料）飞片，结合已有的测速实验数据，验证模拟结果的准确性和可靠性，在此基础上，分析IEMF与飞片的厚度、半径、速度、倾斜角度以及线圈孔径之间的关系，以期为冲击实验中的示波器设置提供参考。

1. MAVIS装置

1.1 MAVIS结构

MAVIS主要由不锈钢壳体、永磁环、导磁片（软铁磁材料）和检测线圈组成。检测线圈一般包括3组甚至更多的测量单元，测量单元之间等距配置，每组测量单元包含1个检测线圈、永磁环和紧贴永磁环的导磁片。检测线圈缠绕于导磁片内侧，与导磁片之间相互绝缘。检测线圈由具有绝缘层的铜导线制成。检测线圈的测量单元之间既可以串联连接，也可以相互独立（即每个检测线圈单独与示波器的一个通道连接）。早期，因为示波器的采集通道有限，所以MAVIS的检测线圈主要采用串联方式连接，其弊端是只要任何一个检测线圈断路，就会导致整个MAVIS失效。以3组检测线圈为例，其MAVIS结构如图1所示。

图 1 MAVIS示意图

Figure 1. Schematic diagram of MAVIS

下载: 全尺寸图片幻灯片

1.2 MAVIS的工作原理

具有3组检测线圈的MAVIS的工作原理如图2所示。永磁环在飞片运动路径上产生非均匀磁场。当飞片高速通过MAVIS系统时，根据法拉第电磁感应定律，飞片表面及内部产生涡旋电流（eddy current，EC）。EC激发的磁场与空间原磁场相互作用，导致整个空间中的磁场发生变化，变化的磁场使检测线圈中产生IEMF，示波器记录检测线圈中的IEMF信号。通过此信号波形，能够获得飞片依次通过相邻两个检测线圈的时间间隔，分别记为∆T₁和∆T₂。相邻两个检测线圈间的距离是已知的，记为S。根据

图 2 MAVIS的工作原理

Figure 2. Working principle of MAVIS

下载: 全尺寸图片幻灯片

${v_{{\rm f}1}} = \frac{S}{{\Delta {T_1}}}$

(1)

${v_{{\rm f}2}} = \frac{S}{{\Delta {T_2}}}$

(2)

$v_{\mathrm{f}}=\frac{v\mathrm{_{f1}}+v_{\mathrm{f2}}}{2}$

(3)

即可得到飞片的速度v_f，其中v_f1和v_f2分别为飞片通过前、后检测线圈的平均速度。

2. 有限元模拟

2.1 模型建立

结合有限元方法，利用ANSYS Electromagnetics Suite对MAVIS的测速过程建模。建立结构、尺寸以及材料与第1节MAVIS相同的三维模型，如图3所示。对于壳体和导磁片，采用自然边界条件；飞片、永磁环以及检测线圈呈轴对称分布。采用对称边界条件，选用最大网格尺寸模式依次对壳体、永磁环、导磁片、检测线圈以及飞片进行划分。将最大迭代步数设为5，模型误差率设为0.1%，非线性残差设为0.001，以期最大限度地提高模型的精度。设置检测线圈的参数，每组测量单元的检测线圈匝数均为50匝，内部初始电流为0 A，考虑到独立的检测线圈能够克服串联存在的问题，因此建模过程中检测线圈采用相互独立的连接方式，假定MAVIS具有3组独立的检测线圈。对于求解步长，通常采用0.01 ms作为时基间隔，如果需要更高精度的实验结果，可以将时基间隔进一步缩短，通过增加求解步数的方法提升计算精度。

图 3 MAVIS的三维模型

Figure 3. Three-dimensional model of MAVIS

下载: 全尺寸图片幻灯片

2.2 静磁场分布

为了理解MAVIS的工作原理，需要掌握MAVIS在静态条件下的磁场分布。在静磁场中，磁场强度满足安培环路定律

${\boldsymbol{B}}\left( {x,y,\textit{z}} \right) = \nabla \times {\boldsymbol{A}}\left( {x,y,\textit{z}} \right)$

(4)

${\boldsymbol{B}}\left( {x,y,\textit{z}} \right) = {\mu _{\rm {r}}}{\mu _0}{\boldsymbol{H}}\left( {x,y,\textit{z}} \right)$

(5)

$\nabla \cdot {\boldsymbol{B}}\left( {x,y,\textit{z}} \right) = 0$

(6)

可以得出

$\nabla \times \left[ {\frac{{{\boldsymbol{B}}\left( {x,y,\textit{z}} \right)}}{{{\mu _{\rm {r}}}{\mu _0}}}} \right] = {\boldsymbol{J}}\left( {x,y,\textit{z}} \right)$

(7)

进一步可以得出

$\nabla \times \left[ {\frac{1}{{{\mu _{\rm {r}}}{\mu _0}}}\nabla \times {\boldsymbol{A}}\left( {x,y,\textit{z}} \right)} \right] = {\boldsymbol{J}}\left( {x,y,\textit{z}} \right)$

(8)

式中：B(x, y, z)为磁场中的磁感应强度，H(x, y, z)为磁场强度，A(x, y, z)为矢量磁位，J(x, y, z)为电流密度，μ₀=4π×10⁻⁷ H/m，μ_r为相对磁导率。对于永磁环，磁感应强度为

$\boldsymbol{B}(x,y,\textit{z})=\mu_0\mu_{\rm r}\boldsymbol{H}(x,y,\textit{z})+\mu_0\boldsymbol M_{\rm p}(x,y,\textit{z})$

(9)

式中：M_p (x, y, z)为永磁环的极化强度。

在MAVIS中，采用各向异性的导磁材料，即

${\mu _{\rm r}} = \left( {\begin{array}{*{20}{c}} {{\mu _{{\rm r}x}}}&{}&{} \\ {}&{{\mu _{{\rm r}y}}}&{} \\ {}&{}&{{\mu _{{\rm r}\textit{z}}}} \end{array}} \right)$

(10)

式中：μ_rx、 μ_ry、μ_rz分别为x、 y、z方向上的相对磁导率。对于磁场强度H(x, y, z)，可以描述如下

$\boldsymbol{H}(x,y,\textit{z})=\boldsymbol H_{\rm p}(x,y,\textit{z})+\nabla\varphi(x,y,\textit{z})+\boldsymbol H_{\rm c}(x,y,\textit{z})$

(11)

式中： $\varphi$ (x, y, z)为磁势，H_p(x, y, z)为四面体网格6条边上的磁场强度，H_c (x, y, z)用于描述永磁环上的磁场强度。理论上，利用上述关系可以计算出空间中任意点的磁场强度^[41]。

2.3 静磁场模拟

MAVIS结构复杂，难以通过理论计算空间中任意一点的磁场强度。利用ANSYS Electromagnetics Suite可以得到磁感应强度分布，模拟结果如图4所示。空间中的磁场集中分布于永磁环两端，说明导磁片能够改变空间中的磁场分布，使磁场在永磁环的两端集中分布。

图 4 静态空间中的磁场分布

Figure 4. Distribution of the static magnetic field

下载: 全尺寸图片幻灯片

2.4 测速模拟

采用动态计算模块模拟MAVIS的测速过程。设飞片以恒定的速度运动，整个飞行时间由飞片速度以及飞片运动路径共同决定。

讨论了沿飞片厚度和半径方向的EC分布形式。对于沿飞片厚度方向的EC，由欧姆定律， ${\boldsymbol{J}} = \sigma {\boldsymbol{E}}$ ，其中E为电场强度， $\sigma$ 为飞片的电导率。两边取散度，可得

$\nabla \cdot {\boldsymbol{J}} = \frac{\sigma }{\varepsilon }q$

(12)

式中： $\varepsilon$ 为飞片的相对介电常数，q为所选区域总的电荷量。结合麦克斯韦方程组，由于飞片上只存在EC，并且EC满足电流扩散方程，因此

${\nabla ^2}{\boldsymbol{J}} - {k^2}{\boldsymbol{J}} = 0$

(13)

式中：k为常数。假定飞片的运动方向沿z轴正方向，利用式(13)，结合EC径向呈轴对称分布的特点，可得飞片厚度d处xy平面内的EC密度的通解

$J_x(d)=J_y(d)=C_1\mathrm{e}^{-kd}+C_2\mathrm{e}^{-kd}$

(14)

式中：J_x(d)和J_y(d)分别为沿x和y方向的EC密度分布，C₁和C₂为常数。实际中，由于J_x(d)和J_y(d)不可能无限大，因此C₂=0。对于飞片表面中心位置d=0处，假定J_x(0)=J_y(0)=J₀，可以得出

$J_x(d)=J_y(d)=J_0\mathrm{e}^{-\alpha d}\mathrm{e}^{-{\rm {j}} \beta d}$

(15)

即飞片的EC密度，其中J₀为飞片表面中心位置的EC密度，α和β为表征飞片表面磁通量变化率的常数。

对于飞片EC的径向分布，根据电磁场理论，EC的电流密度与飞片的电阻率、磁导率、飞片厚度、飞片半径以及飞片表面的磁通量变化率有关。当飞片与永磁环之间的距离z=z₁时，特定厚度处飞片径向分布的EC密度J的大小与径向r以及磁通量变化率Δ $\varphi$ 有关^[42]，即

${\boldsymbol{J}} = {\left( {\frac{{\partial {\boldsymbol{F}}\left( {\textit{z},r,\Delta \varphi } \right)}}{{\partial r}}} \right)_{\textit{z} = {\textit{z}_1}}} = {\boldsymbol{f}}(r)$

(16)

式中：F(z, r, Δ $\varphi$ )为描述飞片不同位置处磁通量的函数，f(r)为电流密度分布函数。空间磁场的非均匀分布导致EC在飞片内非均匀分布。由于难以测定空间磁场分布，因此飞片表面及内部的EC密度的大小与分布难以量化。为此，有必要对MAVIS进行三维建模，通过模拟MAVIS的整个测速过程，获得检测线圈中的AIEMF（amplitude of induced electromotive force）。

2.5 模拟与实验对比

选用Al作为飞片材料，相应的实验数据来源于西南交通大学的二级轻气炮实验，发射管口径为25 mm，电磁感应测速装置安置于炮口末端，允许通过的飞片最大半径为17 mm。所用的Al飞片的厚度d=2.5 mm，半径r=12 mm。实验中，实测飞片速度为2 659.6 m/s。在模拟过程中，将飞片运动速度同样设置为2 659.6 m/s，模拟相同条件下Al飞片通过MAVIS的过程。模拟与实验得到的IEMF波形如图5所示。

图 5 模拟与实验IEMF信号的对比

Figure 5. Comparison of simulated and experimental IEMF signals

下载: 全尺寸图片幻灯片

实验与模拟所得波形的变化趋势基本相同，模拟波形与实验波形存在时间延时，这是由于模拟过程中设置的飞片运动起始时间与实验不同导致的，该延时对测速结果无影响。结合图5可以获得飞片通过MAVIS相邻检测线圈的间隔时间Δt₁和Δt₂。相邻两线圈之间的距离均为50 mm。实验测得Δt_1(exp.)=18.7 μs，Δt_2(exp.)=18.9 μs；模拟计算得到Δt_1(sim.)=18.8 μs，Δt_2(sim.)=18.9 μs。结合式(1)可以计算出实验中的飞片速度v_f(exp.)=2 659.6 m/s，模拟中的飞片速度v_f(sim.)=2 652.5 m/s，二者仅存在0.2%的相对误差。开展了多次不同飞片速度的验证实验，通过模拟波形得到的速度相对误差均小于0.4%，低于MAVIS的系统误差（约0.9%）。将AIEMF定义为最大正向电动势与最大反向电动势的绝对值之和。通过比较实验和模拟波形，得到3组AIEMF（A₁、A₂和A₃），结果见表1。实验得到的AIEMF平均值为3 166.6 mV，模拟的平均值为3 209.0 mV，二者的相对误差约为1.1%，低于MAVIS系统的幅值误差2.0%。

表 1 AIEMF的实验和模拟结果

Table 1. Simulated and experimental results of AIEMF

Method	A₁/mV	A₂/mV	A₃/mV
Simulation	3296.7	3234.1	3096.2
Experiment	3251.4	3208.7	3051.6

下载: 导出CSV

| 显示表格

3. AIEMF的影响因素

3.1 飞片厚度

在冲击实验中，飞片厚度d一般介于2.5～3.5 mm之间。然而，在某些实验中，往往需要减小飞片厚度。例如，超高速发射往往是通过减薄飞片厚度从而降低弹重实现的。为此，模拟时有必要拓宽飞片厚度范围，将飞片的厚度范围扩宽至0.5～3.5 mm，探究不同发射速度（v）下飞片厚度的变化对AIEMF的影响。在模拟与实验中，飞片的半径r=12 mm。模拟与实验得到的不同发射速度下AIEMF随飞片厚度的变化如图6所示。从图6中不难看出，在特定的发射速度下，AIEMF随飞片厚度的增加呈非线性波动，整体呈增大趋势。当飞片速度大于2 km/s时，AIEMF的波动幅度显著；当飞片厚度在3.3 mm附近时，AIEMF达到峰值。整体上看，实验与模拟所得AIEMF的变化趋势一致。

图 6 模拟与实验得到的不同发射速度下AIEMF随飞片厚度的变化

Figure 6. Simulated and experimental AIEMF depending on the thickness of flyer at various launch velocities

下载: 全尺寸图片幻灯片

3.2 飞片半径

对于本研究的二级轻气炮，其飞片半径r介于10～12 mm之间。然而，在某些实验中也可能会使用较小半径的飞片，为此，模拟了不同半径的飞片产生的AIEMF。模拟时，将飞片厚度定为2.5 mm，结果如图7所示。

图 7 模拟与实验得到的不同发射速度下AIEMF随飞片半径的变化

Figure 7. Simulated and experimental AIEMF depending on the flyer radius at various launch velocities

下载: 全尺寸图片幻灯片

研究发现，为了使检测线圈中产生明显的IEMF，飞片半径必须大于某临界值，该临界值随着飞片速度的变化而变化。当飞片厚度为2.5 mm时，飞片速度为1、2、3、4、5和6 km/s对应的飞片半径临界值分别为8.8、9.2、9.4、9.4、9.6和9.6 mm。当飞片半径大于半径临界值时，随着半径的增加，检测线圈中的AIEMF起初呈线性增大趋势，随后呈非线性增大特征。

3.3 飞片倾角

实验中，飞片安装和弹丸制作不可避免地存在偏差，导致飞片在发射过程中出现倾斜，即出现倾斜角θ（见图8），飞片的倾斜会对冲击加载实验产生较大的影响。目前，较为成熟的飞片倾角测量方法有光测法和电探针法，但是，这两种方法并不能适用于所有冲击加载实验。如果能通过比较感应电动势幅值的大小获得飞片倾角，那么该方法将会获得广泛的运用。为了达到这一目标，有必要利用模拟探究飞片倾角θ对测速结果的影响，尤其是倾角θ对AIEMF以及测速时间间隔（飞片经过相邻两组测量单元的时间间隔）的影响。

图 8 飞片倾角示意图

Figure 8. Schematic diagram of tilt angle

下载: 全尺寸图片幻灯片

模拟时，设飞片的厚度d=2.5 mm，飞片半径r=12 mm，倾角θ取0°～20°。模拟得到的不同发射速度下AIEMF和测速时间间隔随θ的变化如图9所示。

图 9 模拟得到的不同发射速度下时间间隔和AIEMF随倾斜角的变化

Figure 9. Simulated time intervals and AIEMF depending on the tilt angles at various launch velocities

下载: 全尺寸图片幻灯片

从图9中可以看出，不同的发射速度下，测速时间间隔随着θ的增加出现微小波动。整体上看，时间间隔并不随θ的增加产生较大的变化。这是因为空间磁场由永磁环和导磁片确定，最大磁通量的变化率由空间磁场确定，这些量在测速过程中不会发生较大的改变，因此时间间隔的波动很小。

从图9中还可以看出，随着θ的增大，AIEMF呈非线性变化，并且整体上随θ增大而减小。当飞片速度为1 km/s时，AIEMF在θ为2°、6°、10°、17°时达到极大值，在θ为5°、7°、16°时达到极小值；当飞片速度为2 km/s时，AIEMF随θ的增加呈波动递减，在θ为15°时出现一个明显的极大值；当飞片速度为3、4、5、6 km/s时，AIEMF的极值点位置较为固定，极大值对应的θ为1°、9°和15°，极小值对应的θ为2°、6°、10°、16°。AIEMF呈现非线性递减的原因在于，随着θ的增加，飞片的有效横截面积减小，同时飞片有效厚度增加，两者相互作用导致AIEMF呈非线性变化。当θ达到一定值时，飞片的有效横截面积成为主要影响因素，飞片有效横截面积的减小使AIEMF整体上趋于减小。

3.4 线圈孔径

由于不同MAVIS系统装配的检测线圈的孔径不同，所以有必要探究不同线圈孔径对AIEMF的影响。选择d=2.5 mm、r=12 mm的Al飞片进行模拟，检测线圈的孔径D从24 mm增加至34 mm。图10显示了v=3 km/s时MAVIS其中一个检测线圈记录的IEMF信号。不同发射速度下检测线圈中的AIEMF随线圈孔径的变化如图11所示。

图 10 模拟得到的不同孔径的检测线圈的IEMF信号（v=3 km/s，d=2.5 mm，r=12 mm）

Figure 10. Simulated IEMF signals produced by the pick-up coils with different diameters (v=3 km/s, d=2.5 mm, r=12 mm)

下载: 全尺寸图片幻灯片

图 11 模拟得到的不同速度下AIEMF随检测线圈孔径的变化（d=2.5 mm，r=12 mm）

Figure 11. Simulated AIEMF at various launch velocities produced by the pick-up coils with different diameters (d=2.5 mm, r=12 mm)

下载: 全尺寸图片幻灯片

如图10所示，检测线圈孔径的改变只会改变AIEMF，并不改变信号所对应的时间宽度。从图11不难看出，随着检测线圈孔径的增加，检测线圈中的AIEMF逐渐衰减，并且衰减速率也随着检测线圈孔径的增加而减小，最终趋于零。飞片速度越大，相同孔径下对应的AIEMF衰减量越大。

3.5 飞片速度

为了探究Al飞片速度对AIEMF的影响，选择多组厚度和半径的飞片参数进行组合模拟，结果如图12所示。对比模拟与实验结果，不难看出，在相同条件下，随着飞片速度的增大，AIEMF逐渐增大。图12(a)中，飞片半径为10 mm，当飞片速度为5 km/s时，AIEMF出现了极小值。图12(b)～图12(f)中，AIEMF随着速度的增大而近似线性增大。

图 12 模拟与实验得到的AIEMF随飞片速度的变化

Figure 12. Dependence of the simulated and experimental AIEMF with the launch velocity of flyer

下载: 全尺寸图片幻灯片

3.6 常用材料的模拟

常用的飞片材料除了Al之外，还有Cu、Ta、W以及93W等一系列材料。为此，模拟了采用Cu、Ta、W和93W飞片时，AIEMF随飞片厚度及半径的变化，结果分别如图13、图14、图15和图16所示。

图 13 模拟与实验得到的AIEMF随铜飞片厚度及半径的变化

Figure 13. Dependence of the simulated and experimental AIEMF with thickness and radius of copper flyer

下载: 全尺寸图片幻灯片

图 14 模拟与实验得到的AIEMF随钽飞片厚度及半径的变化

Figure 14. Dependence of the simulated and experimental AIEMF with thickness and radius of tantalum flyer

下载: 全尺寸图片幻灯片

图 15 模拟与实验得到的AIEMF随钨飞片厚度及半径的变化

Figure 15. Dependence of the simulated and experimental AIEMF with thickness and radius of tungsten flyer

下载: 全尺寸图片幻灯片

图 16 模拟与实验得到的AIEMF随93W飞片厚度及半径的变化

Figure 16. Dependence of the simulated and experimental AIEMF with thickness and radius of 93W flyer

下载: 全尺寸图片幻灯片

3.7 讨　论

对比模拟与实验结果，不难看出，模拟与实验所得AIEMF存在偏差。主要原因是：模拟中未考虑飞片的磁化效应，即飞片接近永磁环时被磁化，磁化效应导致原磁场发生改变，进而使检测线圈产生感应电动势；其次，模拟中只考虑了检测线圈中的EC，未考虑其他部位产生的EC，其他部位产生的EC会导致热损耗。

观察模拟与实验结果，可以发现，飞片依次通过3组测量单元时，AIEMF并不相同，说明飞片通过MAVIS时并非匀速通过。飞片穿过永磁环的过程中，将受到安培力的作用，致使飞片速度减小。通过分析飞片在前后两组测量单元的速度差异，得出飞片速度减小约0.1%，低于MAVIS的系统测量误差（0.9%），因此，飞片速度的减小可以忽略不计。

4. 结　论

结合MAVIS实验装置，利用ANSYS Electromagnetics Suite构建模型，对MAVIS在静态条件下的磁场分布及动态测量过程进行三维模拟。在静态条件下，模拟得到了空间磁场分布，解决了难以通过理论计算空间磁场强度这一难题。在动态过程中，将模拟与实验进行对比，结果表明，通过模拟进行飞片速度测量的相对误差低于0.4%，小于MAVIS的系统测速误差0.9%。实验与模拟得到的AIEMF的相对误差为1.1%，低于系统幅值误差2.0%，说明模拟结果是可靠的。模拟结果表明：AIEMF随飞片厚度的增加呈非线性波动，整体呈增大趋势；在飞片半径大于临界值的条件下，AIEMF随飞片半径的增加先近乎线性增大，后变化趋缓并呈现波动特征；AIEMF随飞片速度的增大而增大；AIEMF随飞片倾角的增大而波动，整体上趋于减小；飞片依次通过3组测量单元时，可以忽略安培力对飞片速度的影响。模拟结果与实验结果存在的误差主要归因于模拟中未考虑测量装置中其他元件产生的EC以及飞片的磁化效应。本研究建立的三维模型及其模拟结果是可靠的，能有效地解决因测速信号不明确而导致测速失效这一难题。

图 1 MAVIS示意图

Figure 1. Schematic diagram of MAVIS

下载: 全尺寸图片幻灯片

图 2 MAVIS的工作原理

Figure 2. Working principle of MAVIS

下载: 全尺寸图片幻灯片

图 3 MAVIS的三维模型

Figure 3. Three-dimensional model of MAVIS

下载: 全尺寸图片幻灯片

图 4 静态空间中的磁场分布

Figure 4. Distribution of the static magnetic field

下载: 全尺寸图片幻灯片

图 5 模拟与实验IEMF信号的对比

Figure 5. Comparison of simulated and experimental IEMF signals

下载: 全尺寸图片幻灯片

图 6 模拟与实验得到的不同发射速度下AIEMF随飞片厚度的变化

Figure 6. Simulated and experimental AIEMF depending on the thickness of flyer at various launch velocities

下载: 全尺寸图片幻灯片

图 7 模拟与实验得到的不同发射速度下AIEMF随飞片半径的变化

Figure 7. Simulated and experimental AIEMF depending on the flyer radius at various launch velocities

下载: 全尺寸图片幻灯片

图 8 飞片倾角示意图

Figure 8. Schematic diagram of tilt angle

下载: 全尺寸图片幻灯片

图 9 模拟得到的不同发射速度下时间间隔和AIEMF随倾斜角的变化

Figure 9. Simulated time intervals and AIEMF depending on the tilt angles at various launch velocities

下载: 全尺寸图片幻灯片

图 10 模拟得到的不同孔径的检测线圈的IEMF信号（v=3 km/s，d=2.5 mm，r=12 mm）

Figure 10. Simulated IEMF signals produced by the pick-up coils with different diameters (v=3 km/s, d=2.5 mm, r=12 mm)

下载: 全尺寸图片幻灯片

图 11 模拟得到的不同速度下AIEMF随检测线圈孔径的变化（d=2.5 mm，r=12 mm）

Figure 11. Simulated AIEMF at various launch velocities produced by the pick-up coils with different diameters (d=2.5 mm, r=12 mm)

下载: 全尺寸图片幻灯片

图 12 模拟与实验得到的AIEMF随飞片速度的变化

Figure 12. Dependence of the simulated and experimental AIEMF with the launch velocity of flyer

下载: 全尺寸图片幻灯片

图 13 模拟与实验得到的AIEMF随铜飞片厚度及半径的变化

Figure 13. Dependence of the simulated and experimental AIEMF with thickness and radius of copper flyer

下载: 全尺寸图片幻灯片

图 14 模拟与实验得到的AIEMF随钽飞片厚度及半径的变化

Figure 14. Dependence of the simulated and experimental AIEMF with thickness and radius of tantalum flyer

下载: 全尺寸图片幻灯片

图 15 模拟与实验得到的AIEMF随钨飞片厚度及半径的变化

Figure 15. Dependence of the simulated and experimental AIEMF with thickness and radius of tungsten flyer

下载: 全尺寸图片幻灯片

图 16 模拟与实验得到的AIEMF随93W飞片厚度及半径的变化

Figure 16. Dependence of the simulated and experimental AIEMF with thickness and radius of 93W flyer

下载: 全尺寸图片幻灯片

表 1 AIEMF的实验和模拟结果

Table 1. Simulated and experimental results of AIEMF

Method A₁/mV A₂/mV A₃/mV

Simulation 3296.7 3234.1 3096.2
Experiment 3251.4 3208.7 3051.6

下载: 导出CSV

参考文献(42)

[1]	孟川舒. 高速磁浮列车测速定位问题综述 [J/OL]. 铁道标准设计(2022-11-28)[2023-07-12]. https://kns.cnki.net/kcms2/article/abstract?v=TzO8JwpG6ugvCh-VGkUDKgKSZPKAKEbopyhecrl89N6uki3CMZ5uouIaNOSIOf8Ya06yC3QPR-3M4jto4JTCewpY3QWpdpyAmYonBoBHY8EntNgbINmCcBTOTp6ELQ3HT9SUsH1gOOQ=&uniplatform=NZKPT&language=CHS. DOI: 10.13238/j.issn.1004-2954.202208010006. MENG C S. Review on problems of speed measurement and positioning of high speed maglev train [J/OL]. Railway Standard Design (2022-11-28)[2023-07-12]. https://kns.cnki.net/kcms2/article/abstract?v=TzO8JwpG6ugvCh-VGkUDKgKSZPKAKEbopyhecrl89N6uki3CMZ5uouIaNOSIOf8Ya06yC3QPR-3M4jto4JTCewpY3QWpdpyAmYonBoBHY8EntNgbINmCcBTOTp6ELQ3HT9SUsH1gOOQ=&uniplatform=NZKPT&language=CHS. DOI: 10.13238/j.issn.1004-2954.202208010006.
[2]	柴博森, 闫东, 王广义, 等. 制动工况桃腔偶合器三维涡特征分析及仿真评价 [J]. 吉林大学学报(工学版), 2023: 1−10. CHAI B S, YAN D, WANG G Y, et al. Three-dimensional vortex characteristic analysis and simulation evaluation of peach cavity hydrodynamic coupling under braking condition [J]. Journal of Jilin University (Engineering and Technology Edition), 2023: 1−10.
[3]	马付远, 张明宇, 徐伟利. 基于多普勒雷达和北斗系统的农业机械测速方法研究 [J]. 农机化研究, 2023, 45(9): 150–154. MA F Y, ZHANG M Y, XU W L. Research on speed measurement methods of agricultural machinery based on Doppler radar and BDS [J]. Journal of Agricultural Mechanization Research, 2023, 45(9): 150–154.
[4]	杨雪亚, 沈显照, 王贤翔. 低空反无雷达的高精度测速方法及波形设计 [J]. 火控雷达技术, 2022, 51(3): 36–39. YANG X Y, SHEN X Z, WANG X X. High-precision speed measurement and waveform design of low-altitude anti-UAV radar [J]. Fire Control Radar Technology, 2022, 51(3): 36–39.
[5]	韩超, 李欣伟, 王凯鑫, 等. 激光测速雷达的性能对比分析 [J]. 测控技术, 2022, 41(9): 90–95. HAN C, LI X W, WANG K X, et al. Performance comparison and analysis of laser velocimetry radar [J]. Measurement and Control Technology, 2022, 41(9): 90–95.
[6]	刘琳, 王彩霞, 田杨萌, 等. 多普勒效应对雷电源辐射定位的影响关系研究 [J]. 电波科学学报, 2022, 37(6): 1049–1056. LIU L, WANG C X, TIAN Y M, et al. Study on the influence of Doppler effect on the radiation localization of lightning sources [J]. Chinese Journal of Radio Science, 2022, 37(6): 1049–1056.
[7]	谢孟龙. 机动车雷达测速仪检定方法探究及不确定度评定 [J]. 计量与测试技术, 2023, 50(1): 108–111. XIE M L. Research on verification method and uncertainty evaluation of vehicle radar velocimeter [J]. Metrology & Measurement Technique, 2023, 50(1): 108–111.
[8]	陈略, 李文潇, 李海涛, 等. 嫦娥五号探测器精细开环测速及空间引力红移验证设计分析 [J]. 中国科学: 物理学力学天文学, 2023, 53(3): 122–132. CHEN L, LI W X, LI H T, et al. High-accuracy open-loop velocity measurement and gravitational redshift verification design analysis based on the Chang’E-5 probe [J]. Scientia Sinica Physica: Mechanica & Astronomica, 2023, 53(3): 122–132.
[9]	王泽杰, 王潜心, 商天文. 基于北斗三号卫星非差多普勒观测信号的优化测速模型 [J]. 电信科学, 2021, 37(12): 42–50. WANG Z J, WANG Q X, SHANG T W. Optimal velocity estimation model based on non-differential Doppler observations from BDS-3 satellites [J]. Telecommunications Science, 2021, 37(12): 42–50.
[10]	刘天琦, 张浩, 焦名. 结合测速测距的卫星定位法研究 [J]. 铁路通信信号工程技术, 2022, 19(3): 17–20, 43. LIU T Q, ZHANG H, JIAO M. Satellite positioning method in combination with speed and distance measurement [J]. Railway Signalling & Communication Engineering, 2022, 19(3): 17–20, 43.
[11]	张伟. 天文光谱测速导航技术与应用思考 [J]. 导航与控制, 2020, 19(Suppl 1): 64–73. ZHANG W. A study of the navigation technology and application based on astronomical spectral velocity measurement [J]. Navigation and Control, 2020, 19(Suppl 1): 64–73.
[12]	刘健. 皮托管流量计在海洋平台天然气计量中的应用 [J]. 石油化工自动化, 2022, 58(2): 84–87. LIU J. Application of pitot tube flowmeter in natural gas metering on offshore platform [J]. Automation in Petro-Chemical Industry, 2022, 58(2): 84–87.
[13]	经福谦, 陈俊祥. 动高压原理与技术 [M]. 北京: 国防教育出版社, 2006: 34. JING F Q, CHEN J X. Dynamic high-pressure generation principle and related technologies [M]. Beijing: National Defense Education Press, 2006: 34.
[14]	MITCHELL A C, NELLIS W J . Diagnostic system of the Lawrence Livermore National Laboratory two-stage light-gas gun [J]. Review of Scientific Instruments, 1981, 52(3): 347–359. doi: 10.1063/1.1136602
[15]	KAWAI N, TSURUI K, HASEGAWA S, et al. Single microparticle launching method using two-stage light-gas gun for simulating hypervelocity impacts of micrometeoroids and space debris [J]. Review of Scientific Instruments, 2010, 81(11): 115105. doi: 10.1063/1.3498896
[16]	廖文强, 孙燕云, 吴浪, 等. 基于电探针及光探针的冲击加载时刻测量技术 [J]. 电子制作, 2022, 30(20): 90–92, 96. LIAO W Q, SUN Y Y, WU L. Impact loading time measurement technology based on electric probe and optical probe [J]. Practical Electronics, 2022, 30(20): 90–92, 96.
[17]	陈源福. 冲击加载下苯的透光特性 [J]. 吉林化工学院学报, 2022, 39(7): 91–93. CHEN Y F. Light transmission characteristics of benzene under shock compression [J]. Journal of Jilin Institute of Chemical Technology, 2022, 39(7): 91–93.
[18]	孙雨荟, 赵传, 王韫泽, 等. 冲击加载下50SiMnVB钢微观断裂机理试验研究 [J]. 兵器装备工程学报, 2022, 43(11): 1–8. SUN Y H, ZHAO C, WANG Y Z, et al. Mechanism of microscopic fracture of 50SiMnVB steel under impact loading [J]. Journal of Ordnance Equipment Engineering, 2022, 43(11): 1–8.
[19]	ROCHON P, GAUTHIER N. Induction transducer for recording the velocity of a glider on an air track [J]. American Journal of Physics, 1982, 50(1): 84–85. doi: 10.1119/1.12992
[20]	CARPENA P. Velocity measurements through magnetic induction [J]. American Journal of Physics, 1997, 65(2): 135–140. doi: 10.1119/1.18526
[21]	FERENCE J R M, LEMON H B, STEPHENSON R J. Analytical experimental physics [M]. Chicago, USA: University of Chicago Press, 1956.
[22]	向永江. 全光纤激光多普勒测速技术的研究 [J]. 光学学报, 1998, 12(11): 1053–1056. XIANG Y J. Investigation of all fiber laser Doppler velocimetry [J]. Acta Optica Sinica, 1998, 12(11): 1053–1056.
[23]	谭朝勇, 叶欣, 胡永明. 全光纤多普勒测速系统的研究 [J]. 光电技术应用, 2011, 32(2): 276–279. TAN C Y, YE X, HU Y M. Study on all fiber-optic laser Doppler velocimeter [J]. Semiconductor Optoelectronics, 2011, 32(2): 276–279.
[24]	胡昌明, 王翔, 刘仓理, 等. 阵列DPS测量技术在材料动态力学性能研究中的应用 [J]. 爆炸与冲击, 2010, 30(1): 105–108. HU C M, WANG X, LIU C L. Applications of DPS arrays testing technique to dynamic properties study of materials [J]. Journal of Experimental Mechanics, 2010, 30(1): 105–108.
[25]	高春峰, 魏国, 王琦, 等. 一维激光多普勒测速仪/单轴旋转惯导组合车载高程测量方法 [J]. 红外与激光工程, 2022, 51(12): 9–17. GAO C F, WEI G, WANG Q, et al. Height measurement method based on one-dimensional laser Doppler velocimeter and integrated navigation system of single-axis rotation inertial navigation system [J]. Infrared and Laser Engineering, 2022, 51(12): 9–17.
[26]	万镇宇. 基于结构光场的多普勒测速研究 [D]. 武汉: 华中科技大学, 2023. WAN Z Y. Research on Doppler velocimetry based on structured light field [D]. Wuhan: Huazhong University of Science and Technology, 2023.
[27]	王金贵. 气体炮原理与技术 [M]. 北京: 国防教育出版社, 1999: 38. WANG J G. Principle and technology of gas gun [M]. Beijing: National Defense Education Press, 1999: 38.
[28]	FOUCAULT M L. Physical demonstration of the rotation of the Earth by means of the pen dulum [J]. Journal of the Franklin Institute, 1851, 51(5): 350–353.
[29]	KONDO K I, SAWAOKA A, SAITO S. Magnetoflyer method for measuring gas-gun projectile velocities [J]. Review of Scientific Instruments, 1977, 48(12): 1581–1582. doi: 10.1063/1.1134948
[30]	MOODY R L, KONRAD C H. Magnetic induction system for two-stage gun projectile velocity measurements: SAND 84-0638 [R]. Albuquerque, NM, USA: Sandia National Laboratories, 1984.
[31]	施尚春, 陈攀森, 黄跃. 高速弹丸的磁感应测速方法 [J]. 高压物理学报, 1991, 5(3): 205–214. SHI S C, CHEN P S, HUANG Y. Velocity measurement of magnet induced system for projectile [J]. Chinese Journal of High Pressure Physics, 1991, 5(3): 205–214.
[32]	张修路, 齐少鹏, 张智, 等. 超高速加载下电磁感应测速的多物理场有限元模拟 [J]. 科技资讯, 2017, 15(32): 1–6. ZHANG X L, QI S P, ZHANG Z, et al. Multi-physical finite element simulation of electromagnetic induction velocity measurement under ultra-high speed loading [J]. Science & Technology Information, 2017, 15(32): 1–6.
[33]	ROCHA T J, RAMOS H G, RIBEIRO A L, et al. Magnetic sensors assessment in velocity induced eddy current testing [J]. Sensors and Actuators A: Physical, 2015, 228: 55–61. doi: 10.1016/j.sna.2015.02.004
[34]	RAO B P C. Eddy current testing: basics [J]. Journal of Non Destructive Testing & Evaluation, 2011, 10: 7–16.
[35]	HE Y, LUO F, PAN M, et al. Pulsed eddy current technique for defect detection in aircraft riveted structures [J]. NDT & E International, 2010, 43(2): 176–181.
[36]	DOGARU T, SMITH S T. Giant magnetoresistance-based eddy-current sensor [J]. IEEE Transactions on Magnetics, 2001, 37(5): 3831–3838. doi: 10.1109/20.952754
[37]	RIPKA P, VOPÁLENSKÝ M, PLATIL A, et al. AMR magnetometer [J]. Journal of Magnetism and Magnetic Materials, 2003, 254/255: 639–641.
[38]	RIBEIRO A L, RAMOS H G. Inductive probe for flaw detection in non-magnetic metallic plates using eddy currents [C]//2008 IEEE Instrumentation and Measurement Technology Conference. Victoria, BC, Canada, 2008: 1447−1451.
[39]	POSTOLACHE, RAMOS H G, RIBEIRO A L. Detection and characterization of defects using GMR probes and artificial neural networks [J]. Computer Standards & Interfaces, 2011, 33(2): 191–200.
[40]	JIAO S, CHENG L, LI X, et al. Monitoring fatigue cracks of a metal structure using an eddy current sensor [J]. EURASIP Journal on Wireless Communications and Networking, 2016: 188.
[41]	郭硕鸿. 电动力学 [M]. 北京: 高等教育出版社, 1997. GUO S H. Electrodynamics [M]. Beijing: China Higher Education Press, 1997.
[42]	吴桔生. 电涡流技术与运用 [M]. 长沙: 中南大学出版社, 2014. WU J S. Eddy current technology and application [M]. Changsha: Central South University Press, 2014.

施引文献

资源附件(0)

访问统计

点击查看大图

图(16) / 表(1)

计量

文章访问数: 165
HTML全文浏览量: 88
PDF下载量: 23

1. MAVIS装置
1.1 MAVIS结构
1.2 MAVIS的工作原理
2. 有限元模拟
2.1 模型建立
2.2 静磁场分布
2.3 静磁场模拟
2.4 测速模拟
2.5 模拟与实验对比
3. AIEMF的影响因素
3.1 飞片厚度
3.2 飞片半径
3.3 飞片倾角
3.4 线圈孔径
3.5 飞片速度
3.6 常用材料的模拟
3.7 讨　论
4. 结　论

1. MAVIS装置
1.1 MAVIS结构
1.2 MAVIS的工作原理
2. 有限元模拟
2.1 模型建立
2.2 静磁场分布
2.3 静磁场模拟
2.4 测速模拟
2.5 模拟与实验对比
3. AIEMF的影响因素
3.1 飞片厚度
3.2 飞片半径
3.3 飞片倾角
3.4 线圈孔径
3.5 飞片速度
3.6 常用材料的模拟
3.7 讨　论
4. 结　论

参考文献(42)

施引文献

资源附件(0)

访问统计

针对电磁感应测速系统的有限元模拟

doi: 10.11858/gywlxb.20230693

作者简介: 吴 潇（1996－），男，博士研究生，主要从事动高压下的凝聚态物理研究.E-mail：2215189687@qq.com

通讯作者: 马小娟（1976－），女，博士，教授，博士生导师，主要从事冲击高压下的实验物理研究. E-mail：mxj_swjtu@126.com

计量

出版历程

Finite Element Simulation for Magnet Velocity Induction System

1. MAVIS装置

1.1 MAVIS结构

1.2 MAVIS的工作原理

2. 有限元模拟

2.1 模型建立

2.2 静磁场分布

2.3 静磁场模拟

2.4 测速模拟

2.5 模拟与实验对比

3. AIEMF的影响因素

3.1 飞片厚度

3.2 飞片半径

3.3 飞片倾角

3.4 线圈孔径

3.5 飞片速度

3.6 常用材料的模拟

3.7 讨 论

4. 结 论

计量

出版历程

目录

1. MAVIS装置

1.1 MAVIS结构

1.2 MAVIS的工作原理

2. 有限元模拟

2.1 模型建立

2.2 静磁场分布

2.3 静磁场模拟

2.4 测速模拟

2.5 模拟与实验对比

3. AIEMF的影响因素

3.1 飞片厚度

3.2 飞片半径

3.3 飞片倾角

3.4 线圈孔径

3.5 飞片速度

3.6 常用材料的模拟

3.7 讨 论

4. 结 论

作者简介:
吴　潇（1996－），男，博士研究生，主要从事动高压下的凝聚态物理研究.E-mail：2215189687@qq.com

通讯作者:
马小娟（1976－），女，博士，教授，博士生导师，主要从事冲击高压下的实验物理研究. E-mail：mxj_swjtu@126.com

3.7 讨　论

4. 结　论

3.7 讨　论

4. 结　论