身管外壳对电磁轨道发射装置发射性能影响分析

翟小飞1 邹 锟1,2 李配飞1 刘 华1 彭之然1

(1. 海军工程大学舰船综合电力技术国防科技重点实验室 武汉 430033 2. 东南大学电气工程学院 南京 210096)

摘要 电枢在发射过程中,瞬时变化的导轨电流会在金属身管外壳上感应出巨大的涡流,身管涡流会产生额外的损耗,同时外壳涡流的去磁效应会改变电枢、导轨等器部件的应力分布从而影响装置发射性能。该文建立了电磁轨道发射装置的数学模型,并针对整体式、上下分断式以及叠压式三种外壳结构有限元模型进行了电磁场-结构场联合仿真,获得了发射装置电磁参数和各部件应力分布。有限元电磁仿真结果表明,叠压式外壳结构外壳涡流最小、电感梯度最大、器部件应力最大,整体式外壳结构外壳涡流最大、电感梯度最小、器部件应力最小。全系统电气仿真和对比试验证明,在满足外壳支持强度和装置各部件应力条件下,选用高磁导率、低电导率的材料并设计抑制涡流的身管外壳结构,有利于提高发射装置电枢出口速度和系统效率,从而获得优良的发射性能。

关键词:电磁轨道发射装置 外壳结构 发射性能 电流密度分布 应力分布

0 引言

电磁轨道发射装置(Electromagnetic Launcher, EML),简称发射装置,是利用电磁力推动电枢在膛内高速运动的新型发射装置,由于没有传统火药武器的声滞制约,可以突破传统火药武器的出口速度限制,实现超高速的出口速度[1-4]。身管外壳(简称外壳)作为发射装置的一部分,主要对导轨提供预紧力约束导轨外扩变形,同时为导轨提供支持刚度以保障身管的直线度。金属具有良好的强度和刚度,当前一般采用金属作为外壳材质,由于发射装置工作在百kA、ms级短时脉冲工作模式下,电流短时快速变化会在金属外壳产生巨大的涡流[5-6]。外壳涡流磁场对装置内部的磁场具有一定的削弱作用,也会影响导轨电流分布,从而对装置的电气参数和力学性能产生影响,影响电枢出口速度等发射性能指标。由于外壳结构和材料均会影响外壳涡流路径和大小,因此需要分析外壳材料、结构对发射装置电磁参数、电枢和导轨的受力情况的影响,从而获得身管外壳对发射性能的影响规律。

文献[7]对电磁轨道炮上的瞬态涡流场进行了建模和仿真。文献[8]基于二维静态模型分析了不同管身材料、结构和尺寸对电感梯度的影响,表明管身材料电导率越高装置电感梯度越小,但对装置其他部件(如导轨、电枢)受力没有进行分析。文献[9-10]通过仿真分析了身管外壳涡流分布,研究了材料电磁特性和外壳封装厚度对外壳损耗抑制机理。然而,外壳涡流对系统电气参数和应力参数等影响发射性能的综合分析较少。

为此,本文选取整体式、上下分断式及叠压式三种形式的外壳进行了电磁性能和力学性能的分析。建立了考虑外壳涡流的装置磁场控制方程,同时建立电磁场-结构场有限元(Finite Element Method, FEM)联合仿真模型,对三种结构外壳的发射装置电磁参数及应力分布进行了仿真,并建立了全系统电气仿真模型。仿真结果对比表明,减小外壳涡流有利于提高电枢出口速度,但会导致导轨及外壳等部件受力增加。因此,对外壳材料选取及外壳结构进行设计时要进行电磁性能和力学性能综合考虑,才能获得外壳损耗小、器部件应力低等综合性能优良的发射装置。

1 电磁发射装置数学模型

电磁轨道发射装置主要由电枢、上下导轨以及身管外壳组成,导轨与外壳间设置有绝缘层。身管主要用于支撑导轨和电枢,同时提供抵御导轨外扩的预紧力。电磁轨道发射装置的工作原理是通电导轨在导轨间产生磁场,身管内部磁通密度为B1,身管外部磁通密度为B2,通电电枢受到磁场中电磁力作用,从而沿着导轨直线运动。发射装置结构及受力示意图如图1所示。

width=167.8,height=110.2

图1 发射装置结构及受力示意图

Fig.1 Structure and Force diagram of the EML

根据轨道发射装置理论,电枢电磁推力Fa

width=41.95,height=28.8 (1)

其中

width=35.05,height=28.15

式中,width=11.9,height=11.9为装置电感梯度;L为导轨电感;i为导轨电流。可见,在电流i一定的条件下,电磁推力正比于电感梯度,因此电感梯度是装置的核心参数之一[11]

磁场控制方程如下:

忽略位移电流,根据麦克斯韦方程组为

width=63.25,height=73.25 (2)

式中,H为磁场强度(A/m);B为磁感应强度(T);E为电场强度(V/m);J为电流密度(A/m2);D为电位移矢量(C/m2)。

本构关系为

width=80.15,height=33.2 (3)

式中,v为电枢速度;s 为电导率(S/m);m 为磁导率(H/m)。引入矢量磁位A和标量电位函数j,关系为

width=78.25,height=46.95 (4)

考虑到电流密度J

width=40.05,height=15.05 (5)

式中,Js为外加电流密度;Je为涡流密度。根据式(2)~式(5),可以得到带速度项vA-j 磁场控制方程[12-13]

width=202.25,height=65.1(6)

求解式(6)所示的磁场控制方程,可以得到矢量磁位A和磁感应强度B随空间、时间的分布函数,从而获得导轨和电枢上的电流密度J以及外壳上的涡流密度Je等物理量。根据求得各个部件的电流密度J和磁感应强度B,可以计算各部件的体积力密度fa同各个部件的合力F关系为

width=132.1,height=28.15 (7)

式中,W 为各部件的体积积分区域。同理,根据各部件的电流密度J和电导率,可以通过体积分获得相应部件的欧姆损耗P

width=72,height=31.95 (8)

电枢发射过程中,电流密度J和磁感应强度B均随时间和空间不断变化,因此要获得精准的时空分布和各部件的电磁推力和损耗等数值,需要引入有限元模型进行数值仿真[14-15]

2 三种结构形式外壳的发射装置电磁性能

发射装置有限元模型包括了上下导轨、电枢以及外壳组成,其三维结构及其尺寸如图2所示。

为了对比研究不同形式的外壳涡流对发射性能的影响,采用整体式、上下分断式以及叠压式三种外壳结构形式。三种外壳电磁参数见表1。

width=217.25,height=225.4

图2 电磁轨道发射装置三维模型结构示意图

Fig.2 3D structure of the EML

表1 三种外壳材料的参数

Tab.1 Material parameters of the three shell

外壳结构材料电导率/(MS/m)导磁性能 整体式导磁硅钢4.5B-H曲线 上下分断式钛合金0.6不导磁 叠压式导磁硅钢4.5B-H曲线

导轨采用高电导率、高强度的铜合金材料,电枢为铝合金6061。

根据文献[16-17]可知,影响导轨电流密度分布的速度趋肤效应、邻近效应均可以采用涡流场扫频的方式在Ansys软件中对三维有限元仿真模型进行等效分析。根据文献[16]电枢运动过程中速度所引起的电流趋肤所对应的频率fv(称为速度频率),计算公式为

width=41.95,height=28.15 (9)

式中,l 为电枢与导轨接触的长度;v为电枢速度。根据图2所示尺寸l=30 mm,设电枢最高出膛速度v=1 500 m/s时,对应的速度频率fv=5 066 Hz,其作为涡流场频率上限。

利用装置的结构对称特性,建立三维有限元仿真1/4结构模型,采用三维涡流场进行电气性能分析,激励电流峰值为100 kA,频率范围设置为50 Hz~5 kHz,有限元三维仿真模型如图3所示。三种结构的差别主要是通过设置不同的绝缘面和材料电气属性(整体或叠片)分别进行区分实现。

width=226.65,height=104.55

图3 有限元三维仿真模型(1/4部分)

Fig.3 3D FEM model of EML (quarter part)

2.1 整体式外壳的电磁性能

整体式外壳为如图2c所示的截面且沿着轴向延伸的一体式结构,材料为导磁的硅钢为发射装置提供良好导磁路径。外壳内侧(即与导轨接触面)设置为绝缘边界,以实现导轨与外壳间的绝缘,简化有限元计算模型。

采用涡流场扫频作为模型激励源,在最高频率5 kHz时电流密度分布及外壳上感应出来涡流电流密度矢量分布云图如图4所示。

width=179.7,height=147.15

(a)电流密度分布云图

width=179.7,height=144.65

(b)外壳电流密度矢量云图

图4 整体式外壳、电流密度分布及外壳涡流矢量分布云图

Fig.4 Current density and eddy current vector description with the integral structure shell

从图4a可知,由于趋肤效应导轨电流主要集中在导轨上下表层,同时邻近效应使电流更集中在导轨内侧表面。从图4b所示的涡流矢量云图可知,外壳涡流左右对称分布,且与导轨+电枢构成的传导电流回路方向相反,可见外壳涡流对导轨空间磁场具有削弱作用,外壳涡流电流密度最大值约为1.0× 108 A/m2。外壳涡流会反过来影响导轨的电流密度分布和电感梯度,整体式外壳装置电感梯度随频率变化曲线如图5所示。

width=210.35,height=93.9

图5 整体式外壳,电感梯度及外壳损耗随频率变化曲线

Fig.5 Inductance gradient and loss of the shell curves with the integral structure shell

从图5可知,整体式外壳条件下,发射装置电感梯度随频率增加快速下降,在f=5 kHz时电感梯度只有0.242 mH/m左右。外壳损耗随频率增加而增加,外壳损耗在频率f=5 kHz时高达2.3 MW左右,可见,整体式外壳损耗巨大,这将大大降低装置的发射效率。

2.2 上下分断外壳的电磁性能

按照分割涡流回路从而减小涡流的基本思路,将整体式外壳沿水平轴向进行分割,形成上下各半的独立结构形式,其截面如图6a所示。在实际的装置中,采用两个上下独立外壳之间布置绝缘纸的方式进行上下外壳间绝缘。为了保证外壳的支撑刚度,通常利用绝缘后的螺栓将上下两瓣的外壳连接起来,从而实现向导轨提供预紧力的目的[18-19]。外壳上下分断的仿真模型中,外壳上下对接的断面增加了一个绝缘边界面,这样就实现了上下分断外壳涡流的目的。外壳上下分断的仿真结构模型如图6b所示。

width=226.65,height=122.1

图6 上下分断式外壳的有限元模型

Fig.6 The upper-lower separated structure shell FEM model

按照表1所示,外壳采用不导磁、低电导率且低密度的钛合金材料。采用与整体式外壳模型相同的电流峰值和频率范围,对外壳分断模型进行电磁性能分析。频率f=5 kHz时,电流密度分布以及外壳上的电流密度矢量云图如图7所示。

width=179.7,height=321.2

图7 上下分断式外壳、电流密度分布和外壳电流密度矢量云图

Fig.7 Current density and eddy current vector description with the upper-lower separated structure shell

从图7a可以看出,分断式外壳模型中导轨电流分布在导轨上下表面且主要集中在导轨内表面。从图7b可知,由于上下外壳分断并绝缘后,外壳上产生了四个对称涡流,这四个涡流以外壳的四个轴向外沿为中心,且方向与导轨电流方向相反,仍然具有去磁效果。相较于整体外壳,分断的外壳涡流路径几乎增加了一倍,从而使涡流回路总电阻增加了一倍,外壳涡流电流密度最大值约为0.5×108 A/m2,比整体式涡流电流密度最大值降低了近1/2,可见涡流出现了大幅下降。外壳分断时,电感梯度和外壳损耗随频率变化曲线如图8所示。

width=219.75,height=82.65

(a)电感梯度 (b)外壳损耗

图8 分断式外壳电感梯度和外壳损耗随频率f变化曲线

Fig.8 Inductance gradient and loss of the shell curves with the upper-lower separated structure shell

从图8可知,分断式外壳、电感梯度随频率增加不断下降,在f=5 kHz时电感梯度为0.351 mH/m。外壳损耗随频率增加而增加,在f=5 kHz时高达1 100 kW左右,外壳损耗功率虽较整体结构大幅下降但数值仍旧巨大,仍然会降低装置效率。

2.3 叠压式外壳的电磁性能

分断外壳是沿着上下方向将外壳涡流一分为二,由于上半部或下半部外壳仍然是一个整体,故上下外壳涡流依然巨大。借鉴传统电机采用叠压式铁心以减小铁心涡流的思路,采用一定厚度的叠片沿轴向叠压布置,其结构原理如图9所示。同时,外壳采用的导磁金属叠片可以充分利用其导磁特性以降低回路磁阻,从而提高内膛磁场强度和电感梯度。由于叠压式外壳没有轴向支撑刚度,因此实际装置一般情况下还需在叠片外面再布置一个轴向支撑结构以获得优良的轴向支撑刚度。

width=189.1,height=77.65

图9 叠压式外壳结构

Fig.9 Laminate structure shell

叠片截面尺寸如图2c所示,叠片厚度与导轨外扩力以及叠片材料的屈服强度等参数相关,叠片尽量紧密叠压在一起以减小叠片间隙,本文中叠片厚度为1.0 mm。按照表1所示,外壳材料为与整体式外壳相同的导磁硅钢片作为导磁叠片,采用相同的电流峰值和频率范围,对叠压式外壳模型进行电磁性能分析。频率f=5 kHz时,电流密度分布及外壳上的电流密度矢量云图如图10所示。

从图10a可以看出,叠压式外壳结构、导轨电流同样集中分布在导轨上下表面并集中在导轨内侧表面。从图10b可知,外壳中几乎没有沿轴向环流的涡流,且涡流峰值出现数量级式的下降。

width=177.8,height=133.35

(a)电流密度分布云图

width=177.8,height=135.85

(b)外壳电流密度矢量云图

图10 叠压式外壳、电流密度分布和外壳电流密度矢量云图

Fig.10 Current density and eddy current vector description with the laminate structure shell

叠压式外壳结构电感梯度和外壳损耗随频率变化曲线如图11所示。在频率f=5 kHz时,电感梯度为0.424 mH/m。外壳损耗功率随频率增加而增加,在频率f=5 kHz时只有48 kW左右,该数值远远小于前面两种结构形式的身管涡流损耗。

width=216,height=73.25

(a)电感梯度 (b)外壳损耗

图11 叠压式外壳电感梯度及外壳损耗随频率f变化曲线

Fig.11 Inductance gradient and loss of the shell curves with the laminate structure shell

从上述的三种外壳结构形式有限元仿真分析可知,由于速度趋肤效应和邻近效应,导轨电流均集中分布在导轨上下表面且集中在内侧,同时电枢内侧的喉部区域出现了明显的电流密度集中。由于邻近效应,外壳电流密度主要分布在靠近导轨的内表面。三种结构形式的外壳,在f=5 kHz时发射装置电气参数对比见表2。

表2 三种外壳结构形式的发射装置电气参数对比

Tab.2 Electric parameters of the three different shell EML

外壳结构形式电感梯度/(mH/m)外壳损耗/kW 整体式0.2422 300 上下分断式0.3511 100 叠压式0.42448

从表2可以看出,从电气性能来看,采用叠压式外壳发射装置的电感梯度最大且外壳损耗最小,可见叠片外壳的身管结构最有利于装置的发射性能。同时,表2也反映了整体式外壳的发射性能最差,上下分断式外壳电气性能居中。

3 三种结构形式外壳的发射装置应力分析

不同结构形式的外壳对装置磁场分布会产生影响,除了对电磁性能有影响外,还会影响电枢、导轨的应力分布,同时枢轨接触面上的压应力和轨道与外壳接触面的压应力也需要进行应力场仿真计算。

在峰值电流100 kA,频率为5 kHz时,将电磁场仿真获得的各个部件体积力密度导入到结构场模型中,可以获得各部分应力分布和接触面上的应力分布。

3.1 整体式外壳的应力分析

整体式外壳发射装置的主要器部件应力分布云图和接触面接触压力分布云图如图12所示。

3.2 上下分断式外壳的应力分析

上下分断式外壳发射装置的主要器部件应力分布云图和接触面接触压力分布云图如图13所示。

3.3 叠压式外壳的应力分析

叠压式外壳发射装置的主要器部件应力分布云图和接触面接触压力分布云图如图14所示。

从图12~图14可以看出,三种结构形式的应力分布具有相似规律:电枢最大应力出现在喉部内侧,导轨最大应力出现在与电枢接触位置,主要是因为这两个位置均为电流密度高的区域。三种外壳结构形式发射装置受力情况对比见表3。

width=226.65,height=100.15

(a)电枢应力云图 (b)导轨应力云图

width=226.65,height=100.15

(c)枢轨间接触压力 (d)导轨与外壳间接触压力

图12 整体式发射装置器部件和接触面上的应力云图

Fig.12 Stress distribution of the EML with the integral structure shell

width=226.65,height=85.75

(a)电枢应力云图 (b)导轨应力云图

width=226.65,height=89.55

(c)枢轨间接触压力 (d)导轨与外壳间接触压力

图13 上下分断式发射装置器部件和接触面上的应力云图

Fig.13 Stress distribution of the EML with the upper-lower separated structure shell

width=226.65,height=89.55

(a)电枢应力云图 (b)导轨应力云图

width=226.65,height=90.8

(c)枢轨间接触压力 (d)导轨与外壳间接触压力

图14 叠压式发射装置器部件和接触面上的应力云图

Fig.14 Stress distribution of the EML with the laminated structure shell

表3 三种结构形式发射装置受力对比(5 kHz)

Tab.3 Stress of the three different shell EML (5 kHz)

外壳结构形式电枢应力最大值/ MPa导轨应力最大值/ MPa枢轨接触力/N导轨与外壳接触力/N 整体式14.112.091651.3957.5 上下分断式17.193.364884.72 460.6 叠压式22.515.16905.63 298

从表3各个器部件应力及其接触力对比,叠压式最大而整体式最小。这主要是叠压式涡流最小,对装置内部的磁场去磁效应最弱,导致内部磁场最强,从而引起各部件应力以及接触力最大。

4 系统仿真和试验

发射系统主要由脉冲电源和发射装置组成。脉冲电源由电容C、晶闸管SCR、续流二极管VD等主要元器件组成。发射过程中,由于电枢运动将导轨不断接入到电气回路中,发射装置可以等效为电阻R、电感L随电枢位移x不断增加的动态阻感性负载。同时,身管上的涡流区域也会随着导轨通流区域增加而增加,因此也将其等效为随电枢位移x而变化的电阻Rshll。电磁发射系统电气仿真系统框图如图15所示。

width=368.75,height=160.3

图15 电磁发射系统电气仿真系统框图

Fig.15 Diagram of the electrical system simulation

根据文献[20]所提出的仿真参数提取原理,以及速度频率式(9)将电枢速度大小v转化为等效频率f,图15中的参数如电感梯度L'、电感系数Lu、导轨电阻系数Ru及外壳等效电阻系数Rush等参数随频率f变化曲线,均可根据三维模型涡流场扫频获得,本文就不再赘述。利用Simplore软件构建了三种结构形式身管外壳发射装置系统仿真模型,图16为三种仿真结果对比波形。

width=224.75,height=83.9

(a)电流 (b)速度

图16 三种结构形式的身管外壳发射工况仿真结果

Fig.16 Simulation launching results of the three different structure shell

从图16所示的仿真结果可以看出,相同发射工况条件下,整体外壳式电枢出口速度最低为680.8 m/s,分断式外壳出口速度次之为769.5 m/s,叠压式外壳出口速度最高为883.1 m/s,说明了叠压式外壳发射装置发射效果最好。

采用叠压式身管进行了一次发射试验,速度曲线采用沿身管中线轴向布置的B点探头进行数据拟合获得,电流采用罗氏线圈测量,试验与仿真对比结果如图17所示。

width=228.5,height=88.3

(a)电流 (b)速度

图17 叠压式身管发射工况试验与仿真曲线

Fig.17 The comparison of test and simulation

本次试验电枢质量为15.3 g,电容C容值为240 mF,系统仿真和试验结果对比见表4。

表4 仿真与试验对比

Tab.4 Results of the simulation and the experiment

参 数数 值 仿真试验 电容电压/V1 4001 400 峰值电流/kA196.5197.4 出口速度/(m/s)883.1889.9

从表4可以看,仿真和试验数据曲线较为吻合,出口速度误差为0.7%,电流峰值误差为0.4%。

5 结论

根据发射装置的有限元仿真模型结果可以看出,外壳材料和结构对外壳涡流产生了巨大影响,叠压式结构的外壳涡流最小,上下分断式的外壳涡流次之,整体式外壳涡流最大。外壳涡流差别导致装置电磁性能和应力大小存在以下差异:

1)身管涡流越小如叠压式结构,涡流的去磁效应越弱,电感梯度越高,电枢推力越大,同时外壳损耗越小,装置发射效率越高;反之,结论也成立。

2)身管涡流越小如叠压式结构,各元器件应力越大,接触面接触力越大,对器件的强度要求越高;反之,结论也成立。

从装置工作原理可知,导通电流一定条件下,磁场越强将会对位于磁场内的通电导体产生越大的电磁力:电枢就是向前的电磁推力,上下导轨就是外扩力,尾部馈电电缆就是向后的电磁力。导磁性能优良的外壳材料有利于增加装置内部磁场,而外壳涡流有减弱装置内部磁场的作用,这两种作用是相反的。可见,外壳采用高磁导率、低电导率的材料,越能获得较高的电磁推力和出口速度。有效的抑制涡流措施包括:①选取相对磁导率高的铁磁材料并使其工作在非饱和区域;②选择低电导率甚至不导电的材料;③采用降低身管涡流的结构,如本文中采用的上下分断结构或叠片叠压结构。

从上述分析可知,减小外壳涡流有利于增加电枢电磁推力,但也会相应地增加导轨电磁外扩力和尾部馈电装置的受力,需要增加外壳结构强度以提供更高的预紧力。因此,外壳的结构设计和材料选取时要同时兼顾电磁性能和力学性能,才能获得综合性能优良的发射性能。

参考文献

[1] 马伟明, 鲁军勇, 李湘平. 电磁发射超高速一体化弹丸[J]. 国防科技大学学报, 2019, 41(4): 1-10.

Ma Weiming, Lu Junyong, Li Xiangping. Electro- magnetic launch hypervelocity integrated projectile[J]. Journal of National University of Defense Technology, 2019, 41(4): 1-10.

[2] 王莹, 肖峰. 电炮原理[M]. 北京: 国防工业出版社, 1995.

[3] 马伟明. 关于电工学科前沿技术发展的若干思考[J]. 电工技术学报, 2021, 36(22): 4627-4636.

Ma Weiming. Some thoughts on the development of cutting-edge technology in electrotechnics[J]. Transa- ctions of China Electrotechnical Society, 2021, 36(22): 4627-4636.

[4] Ma Weiming, Lu Junyong. Thinking and study of electromagnetic launch technology[J]. IEEE Transa- ctions on Plasma Science, 2017, 45(7): 1071-1077.

[5] 彭之然, 刘华, 汪光森, 等. 基于多模块分时放电半解析模型的脉冲电源触发策略[J]. 电工技术学报, 2021, 36(增刊1): 54-61.

Peng Zhiran, Liu Hua, Wang Guangsen, et al. Trigger strategy of pulse power supply based on semi- analytical model of multi-module time-sharing discharge[J]. Transactions of China Electrotechnical Society, 2021, 36(S1): 54-61.

[6] Landen D, Satapathy S. Eddy current effects in the laminated containment structure of railguns[J]. IEEE Transactions on Magnetics, 2007, 43(1): 150-156.

[7] 林庆华, 栗保明. 有限元/边界元耦合法计算电磁轨道炮三维瞬态涡流场[J]. 南京理工大学学报, 2010, 34(2): 217-221.

Lin Qinghua, Li Baoming. Finite-element/boundary- element coupling method for 3D transient eddy current field calculation in electromagnetic railgun[J]. Journal of Nanjing University of Science and Tech- nology, 2010, 34(2): 217-221.

[8] 刘守豹, 阮江军, 黄道春, 等. 封装对轨道炮电感梯度的影响[J]. 电工电能新技术, 2009, 28(4): 42- 45, 75.

Liu Shoubao, Ruan Jiangjun, Huang Daochun, et al. Influence of shielding on rail Gun induction gradi- ent[J]. Advanced Technology of Electrical Engin- eering and Energy, 2009, 28(4): 42-45, 75.

[9] 楼宇涛, 栗保明. 管身对中口径电磁轨道炮的影响分析[J]. 强激光与粒子束, 2015, 27(9): 198-203.

Lou Yutao, Li Baoming. Influence of containment for medium caliber electromagnetic railgun[J]. High Power Laser and Particle Beams, 2015, 27(9): 198- 203.

[10] 楼宇涛, 栗保明. 电磁轨道炮封装中涡流损耗的抑制方法[J]. 高电压技术, 2016, 42(6): 1935-1941.

Lou Yutao, Li Baoming. Methods for reducing the loss caused by the eddy current in electromagnetic railgun shielding[J]. High Voltage Engineering, 2016, 42(6): 1935-1941.

[11] 彭之然, 汪光森, 翟小飞, 等. 电磁轨道发射装置时变电感梯度建模与分析[J]. 电工技术学报, 2020, 35(23): 4843-4851.

Peng Zhiran, Wang Guangsen, Zhai Xiaofei, et al. Modeling and analysis of time-varying inductance gradient for electromagnetic rail launcher[J]. Transa- ctions of China Electrotechnical Society, 2020, 35(23): 4843-4851.

[12] 李湘平, 鲁军勇, 张晓, 等. 电磁发射弹丸膛口磁场分布特性分析[J]. 电工技术学报, 2021, 36(3): 525-531.

Li Xiangping, Lu Junyong, Zhang Xiao, et al. Analysis of distribution characteristics of electro- magnetic launcher projectile muzzle magnetic field[J]. Transactions of China Electrotechnical Society, 2021, 36(3): 525-531.

[13] 林庆华, 栗保明. 基于瞬态多物理场求解器的电磁轨道炮发射过程建模与仿真[J]. 兵工学报, 2020, 41(9): 1697-1707.

Lin Qinghua, Li Baoming. Modeling and simulation of electromagnetic railgun launching process based on a transient multi-physical field solver[J]. Acta Armamentarii, 2020, 41(9): 1697-1707.

[14] Peng Zhiran, Zhang Xiao, Zhai Xiaofei, et al. Field- circuit coupled analysis of an electromagnetic launcher sliding contact[J]. IEEJ Transactions on Electrical and Electronic Engineering, 2022, 17(8): 1198-1208.

[15] 李白, 鲁军勇, 谭赛, 等. 高速滑动电接触电枢表面动态磨损过程研究[J]. 电工技术学报, 2023, 38(1): 131-139.

Li Bai, Lu Junyong, Tan Sai, et al. Study on dynamic wear process of armature surface in high-speed sliding electrical contact[J]. Transactions of China Electrotechnical Society, 2023, 38(1): 131-139.

[16] 翟小飞, 李鑫航, 刘华, 等. 电磁轨道发射装置动态电感梯度分析[J]. 国防科技大学学报, 2022, 44(3): 156-163.

Zhai Xiaofei, Li Xinhang, Liu Hua, et al. Analysis of dynamic inductance gradient of electromagnetic rail launcher[J]. Journal of National University of Defense Technology, 2022, 44(3): 156-163.

[17] Lü Qingao, Xiang Hongjun, Lei Bin, et al. Physical principle and relevant restraining methods about velocity skin effect[J]. IEEE Transactions on Plasma Science, 2015, 43(5): 1523-1530.

[18] Tosun N, Ceylan D, Polat H, et al. A comparison of velocity skin effect modeling with 2-D transient and 3-D quasi-transient finite element methods[J]. IEEE Transactions on Plasma Science, 2021, 49(4): 1500- 1507.

[19] 张嘉炜, 鲁军勇, 谭赛, 等. 考虑初始接触压力的滑动电接触界面磁扩散模型[J]. 电工技术学报, 2022, 37(2): 488-495.

Zhang Jiawei, Lu Junyong, Tan Sai, et al. A magnetic diffusion model of electromagnetic launcher con- sidering initial contact pressure[J]. Transactions of China Electrotechnical Society, 2022, 37(2): 488- 495.

[20] Zhai Xiaofei, Liu Hua, Peng Zhiran. Research on armature thrust inductance gradient of the elec- tromagnetic rail launcher[J]. IEEE Transactions on Plasma Science, 2022, 50(3): 754-760.

Analyzing Influence of Barrel Shell on Launching Performance of Electromagnetic Rail Launcher

Zhai Xiaofei1 Zou Kun1,2 Li Peifei1 Liu Hua1 Peng Zhiran1

(1. National Key Laboratory of Science and Technology on Vessel Integrated Power System Naval University of Engineering Wuhan 430033 China 2. School of Electrical Engineering Southeast University Nanjing 210096 China)

Abstract As part of the electromagnetic launcher (EML), the barrel shell mainly provides preloading force for the rails to constrain the EML’s expansion and deformation while providing support stiffness for the rail to ensure the straightness of the EML. During the launch process of the armature, the instantaneous change of rail current induces a considerable eddy current on the metal barrel shell of the EML, and the eddy current magnetic field of the shell weakens the magnetic field inside the device. It also affects the current distribution of the rails, thereby affecting the electrical and mechanical properties of the EML and the launching performance, such as the exit speed of the armature. It is necessary to analyze the influence of the shell material and structure on the electromagnetic parameters of the EML and the force of the armature and rails.

Firstly, the magnetic field control equation considering the eddy current of the shell is derived. The finite element model (FEM) of the integral, upper-lower separated, and laminated shell structures is established. After the speed frequency is introduced, the frequency-changed current excitation simulates the velocity skin effect caused by the launching of the armature in the bore. Therefore, the inductance gradient and shell loss with frequency change curve, the electric density distribution of the rails, and the eddy current distribution cloud in the shell can be obtained. Secondly, the simulation of electromagnetic-structural is carried out, and the stress distribution of each component of the EML is obtained. The system electrical simulation is established, and the launching simulation results of the three shell structures are compared. It is shown that the laminated structure shell has the smallest eddy current, the largest inductance gradient and device component stress, and the highest armature exit speed. The integral structure shell has the largest eddy current and shell loss, the smallest inductance gradient and device component stress, and the lowest armature exit speed. Finally, the simulation and test results of the laminated structure shell EML are compared, and the errors of the armature exit speed and the current peak are less than 1%.

The material selection with high permeability and low conductivity and the shell structure design for suppressing eddy currents are conducive to improving the exit speed and system efficiency of the EML. In addition, reducing the eddy current of the shell can increase the electromagnetic thrust of the armature. However, the external expansion force of the rails will be increased, and increasing the structural strength of the shell is required to provide a higher preloading force. Therefore, the shell’s structural design and material selection should consider electromagnetic and mechanical properties for excellent launching performance.

keywords:Electromagnetic rail launcher (EML), barrel structure, launching performance, current density distribution, stress distribution

中图分类号:TM153+.2

DOI: 10.19595/j.cnki.1000-6753.tces.222015

国家自然科学基金项目(92266104)、卓越青年科学基金项目(2021-JCJQ-ZQ-004)和173计划重点项目(2022-JCJQ-ZD-128-11)资助。

收稿日期 2022-10-26

改稿日期 2023-02-06

作者简介

翟小飞 男,1982年生,博士,副研究员,硕士生导师,研究方向为电磁发射技术。E-mail: smartnavy@126.com

李配飞 男,1989年生,博士,助理研究员,研究方向电磁发射技术。E-mail: peifeilee@163.com(通信作者)

(编辑 崔文静)