弹体磁屏蔽效应及目标运动状态对磁探测的影响

2024-01-04程晋伟

兵器装备工程学报 2023年12期

李明,程晋伟

(北京理工大学爆炸科学与技术国家重点实验室, 北京 100081)

0 引言

现代战争中,对敌方关键区域如机场、防空阵地及交通枢纽等进行打击与封控是重要的制胜因素,隐埋式封控弹药作为一种重要的威慑封控方式,可在空间和时间上对敌方重要区域进行压制打击,实现封控效果。钢铁作为现代最重要的材质,被广泛应用于各种军事装备之中,使用铁制目标所固有的磁性特征作为其识别信号,具有准确率高、不受战场烟尘、雨雪环境影响,抗干扰能力强的优势[1-2]。且各种新型材质不断被用于磁探测技术,使其探测灵敏度等有很大提升[3-5]。使用基于磁探测的防排技术可显著提高隐埋式弹药的生存与封控能力,增强其作战效能[6]。

目前对于弹体磁屏蔽效应的主要研究方向为对无磁环境如无磁实验室、线缆屏蔽层等需要隔绝外部磁场,加强磁屏蔽效果的屏蔽层改良,应用于弹药领域的磁屏蔽研究则集中于弹体磁屏蔽对于弹体姿态解算的影响,未对减弱弹体的磁屏蔽的方式进行深入研究。壳体屏蔽问题可简化为旋转椭球模型,用以分析磁引信壳体磁屏蔽效应[7],长径比较大的壳体也可简化为无限长圆柱腔体[8];从磁屏蔽角度可分析其对弹体软磁误差的规律[9];在大型磁屏蔽环境屏蔽层上的孔洞漏磁分析表明面积相同的圆形开孔与方形开孔对磁屏蔽效应的影响接近[10];圆柱形磁屏蔽装置上矩形开孔的长宽比及径向与轴向开孔位置对磁屏蔽效应的影响,仿真分析表明面积一定时,开孔长宽比越大,漏磁越少[11]。

对于车辆等目标的磁特征信号研究主要是进行简化仿真[12],对车辆侧向经过时的信号进行了采集分析[13],车辆检测算法方面集中在车流量识别、车型统计等领域[14-16],本文主要对目标在不同运动状态下的磁特征信号变化规律进行研究。

1 隐埋式弹体全屏蔽数理及仿真分析

由于弹体材质的磁导率远高于大气环境及土壤的磁导率,穿过弹药的磁感线会集中于弹体上,使弹体内部磁感应强度减小,产生磁屏蔽效应,弹体的磁屏蔽效应将对弹体内部的磁探测装置产生巨大的干扰。

磁屏蔽系数是描述铁磁材质磁屏蔽效果的一种量化参数。假设不存在铁磁材质磁屏蔽影响时某处的磁场强度为H0,磁感应强度为B0;存在磁屏蔽时的磁场强度为Hs,磁感应强度为Bs,则该铁磁材质对外加磁场的磁屏蔽系数S可以表示为[17]

(1)

本研究中即采用磁屏蔽系数S来评估全屏蔽状态下隐埋式弹体的磁屏蔽效应,磁屏蔽系数S越大则磁屏蔽效果越强。对于长径比较大的弹体,当外界磁场方向垂直于弹体纵轴时,可以将弹体的磁屏蔽效应简化为研究无限长圆柱腔的磁屏蔽效应,如图1所示。

图1 无限长圆柱腔体剖视图

其磁屏蔽系数为

(2)

式(2)中:r1和r2分别为弹体的内、外半径;μr为弹体材质相对磁导率。

(3)

(4)

式(4)中:弹体壁厚δ=r2-r1;弹体平均半径R=(r1+r2)/2。该式表明在满足上述约束条件的情况下,弹体的静磁屏蔽系数与弹体材质的相对磁导率、弹体壁厚成正相关,与弹体的平均半径成负相关,外界环境磁场的大小不影响弹体的磁屏蔽系数。

隐埋式弹体磁屏蔽效应的理论计算可以为其特性分析提供重要参考,但是其只能应用于计算理想弹体,与实际结果差距较大。通过计算机进行有限元仿真,可以设置更多参数更加贴近真实情况,从而更加接近真实结果,为工程和研究提供更加准确的参考依据。本文中主要使用电磁仿真软件对隐埋式弹体的磁屏蔽效应进行建模仿真。建模如图2所示。

图2 隐埋式弹体仿真模型及空气域

在不同的传播介质中,影响磁场及磁信号传播的主要因素是传播介质的磁导率,空气与土壤的相对磁导率均接近1,在仿真设置中可直接用空气域代替隐埋式弹药所处的土壤环境。空气域设置为弹体直径的5倍,为长宽均750 mm的正方形,通过施加矢量磁位的方法在空气域内添加匀强磁场。其形成的磁感应强度为

(5)

图3 匀强磁场

弹体材质对于弹体磁屏蔽效应影响巨大,常用弹体材质如不锈钢、铸铁、45钢、铝等物质相对磁导率由1到数千不等,根据式(4)可知,圆柱形弹体磁屏蔽系数与其外壳磁导率成正比,因此几何参数相同但材质不同的弹体磁屏蔽系数可以相差数千倍。常用的弹体材质主要包括45钢、铸铁、铝等,本研究中拟使用这3种材质的弹体进行仿真,探寻弹体材质对其磁屏蔽特性的影响。弹体设置为外环外径75 mm,内径71 mm,内环外径22 mm,内径19 mm。

经过查询可知,铝的相对磁导率为1,铸铁为60,45钢常用B-H曲线(表示某种铁磁材质在磁化过程中磁感强度B与磁场强度H之间关系的曲线)来表示其磁导率。其最大相对磁导率一般在600～700之间。其仿真结果如图4所示。

图4 3种材质弹体的磁感应强度分布图

从图4中可知,由于铝的相对磁导率接近1,与空气类似,所以其未对匀强磁场产生扰动,无磁屏蔽效应;45钢在此磁场强度作用下相对磁导率大于铸铁,其对匀强磁场的扰动显著强于铸铁,磁屏蔽效应也更加明显。匀强磁场与隐埋式弹体沿X轴(横轴)完全对称,故取弹体上半区域内均匀分布的6个观测点进行采样,各观测点位置坐标及其磁感应强度如图5及表1所示(3种材质磁场磁感应强度单位:μT)。

图5 各观测点位置示意图

表1 3种材质弹体内部磁感应强度

由表1可得,铸铁材质弹体内部平均磁感应强度为27.84 μT,其磁屏蔽系数S铸铁=H0/H铸铁≈2.395,由公式(4)推算出铸铁弹体理论磁屏蔽系数为2.644,仿真值与理论值接近;45钢材质弹体内部平均磁感应强度为7.003 μT,其磁屏蔽系数为S45钢≈9.52,显著高于铸铁材质。

2 弹体开孔优化与目标磁特征信号仿真

2.1 弹体开孔仿真及优化

隐埋式弹体出于结构强度与成本考虑,其弹体主要材质为钢铁材质,钢铁作为铁磁类材质,相对磁导率远远高于空气、土壤等传播介质。由上一小节可知,钢制弹体对磁场和磁信号具有较强的屏蔽作用,会对磁探测装置造成极大的影响。为了减轻弹体磁屏蔽效应对磁探测的影响,本文中选择在弹体侧壁中线处进行开孔,并使用铝合金或奥氏不锈钢等非铁磁类物质对开孔进行填充,恢复弹体的密封性及部分结构强度;同时使得外界磁场的变化可经由开孔进入弹体内部被磁传感器捕获,从而减轻弹体的磁屏蔽效果。为了使磁传感器捕获更多的外界磁信号,其安放位置为紧贴弹体开孔填充物质内壁的中心位置,本文中所述开孔处磁感应强度即该位置磁传感器所捕获的磁感应强度。

隐埋式弹药由空中平台投放后通过伞降降落至地面附近,而后完成战斗部署。在此过程中,隐埋式弹药的自身姿态处于高度不可控状态,致使战斗部完成部署进入封控状态时自身姿态难以确定。定义弹药进入封控状态后,其弹体纵轴与铅垂线之间的夹角为落角φ,纵轴在水平面上的投影与地磁场水平分量夹角为偏角θ;弹体开孔朝向在水平面上的投影与地磁场水平分量的夹角为旋转角γ,这3个角度限定了弹体在落地后的姿态,如图6所示。

其中落角φ和偏角θ共同作用限定了弹体的纵轴朝向,旋转角γ限定了弹体相对于地磁场水平分量的旋转程度;本文中将弹体纵轴在XOZ平面及垂直面上的投影与z轴之间的夹角简称为铅垂角;将旋转角γ简称为水平角。

图6 弹体姿态示意图

表2 北京地区地磁场标准值

隐埋式弹体总体结构为一圆柱体,中间有一条供射流通过的孔洞。其外径为150 mm,内径38 mm,外壁厚度4 mm,内壁厚度3 mm,高度为220 mm。设置外加匀强磁场方向沿X轴正方向,磁感应强度大小为66.67 μT,弹体建模及磁场分布如图7所示。

对弹体进行三维建模仿真时,要控制的主要变量是开孔的大小、水平面上开孔朝向与外界磁场的夹角以及铅垂面上开孔朝向与外界磁场的夹角。使用45钢作为弹体材质,令偏角θ为90°,旋转角γ为0°,落角φ即为铅垂面上开孔朝向与外界磁场的夹角,该夹角及开孔大小对开孔处磁感应强度的影响仿真结果如图8所示;令偏角θ为90°,落角φ为0°,旋转角γ即为水平面上开孔朝向与外界磁场的夹角,该夹角及开孔大小对开孔处磁感应强度的影响仿真结果如图9所示。

图7 弹体三维建模及磁场分布

图8 不同水平角及开孔直径下开孔处磁感应强度大小

图9 不同铅垂角及开孔直径下开孔处磁感应强度大小

由图8及图9可知,在弹体开孔直径小于20 mm且与外界磁场夹角大于60°时,弹体开孔处磁感应强度受到开孔处弹体边界磁感线扩散及外界磁场的综合作用下,其变化趋势较为混乱;弹体开孔直径大于20 mm后,各夹角下开孔处磁感应强度均随开孔直径增加而增加;各夹角下开孔直径在10～20 mm内均可达到开孔直径130 mm时磁感应强度衰减3 dB强度。

当开孔直径大于20 mm时,开孔朝向对开孔处磁感应强度的影响越发明显,夹角为90°与0°时差值可以达到2倍以上,开孔直径为70 mm时,不同水平面夹角下磁感应强度可拟合为

(6)

不同铅垂面夹角下磁感应强度拟合为

(7)

为保证磁特征信号能够尽可能详细地被磁传感器获取,弹体开孔应达到20 mm以上,为减小开孔朝向对磁探测系统的影响,则应以铝合金或奥氏体不锈钢代替磁传感器处的环状弹体,使得周向磁信号均可进入弹体内部。

2.2 目标运动速度对其磁特征信号的影响仿真

对目标运动速度进行仿真时,由于进行三维瞬态仿真计算量大且难以分析,故使用二维瞬态仿真对三维情况进行近似。仿真工况如图10所示。

从仿真波形可以看出，引入虚拟电容后，牵引过程中直流电压及电池电流的震荡程度明显减小，电压最大震荡幅值为±200 V，系统稳定性得到很好的改善。

图10 仿真工况示意图

目标外形尺寸为2 m×2 m×4 m,弹体开孔直径为70 mm,地磁场方向、目标运动方向与弹体开口朝向相同;由于本弹药所针对的目标速度较低,工作速度大多处于0.1～10 km/h,因此仿真目标速度为0.5、1、1.5、2 m/s,不同速度下目标在弹体开孔处产生的磁感应强度如图11所示。

图11 目标不同速度下弹体开孔处磁感应强度

由图11可知,当目标以不同速度接近弹体时,其磁感应强度仅与其所处位置,即目标与弹体之间的距离有关,因此目标速度的不同导致了弹体开孔处磁感应强度变化的斜率不同,而初始位置与最终位置的同步使得不同速度下弹体开孔处的初始与最终磁感应强度相同。

2.3 目标运动朝向对其磁特征信号的影响实验与仿真

本实验中主要为采集不同运动朝向下目标在其前方产生的磁特征信号变化情况。实验场地为学校西操场,实验装置主要包括推车、铁板、磁力计与计算机。本实验中使用推车及铁板来模拟小尺寸目标,推车尺寸为0.9 m×0.6 m×0.9 m,车身材质主体为塑料,把手、螺栓、加强筋及万向轮转向装置为钢制。单块铁板尺寸为0.3 m×0.3 m×0.05 m,质量约为3.53 kg;实验中共使用5块铁板,为体现出目标前后差异,采用前三后二布置,铁板之间通过气泡膜隔开,以减少铁板之间的碰撞与滑动,如图12所示。

图12 目标磁特征信号采集实验环境

本实验中主要验证的是目标运动朝向对其磁特征信号的影响,实验流程如下所示。

磁力计直接放置在地面,3个敏感轴分别指向西、北、天,另一端接入计算机;目标运动朝向自北向南运动轨迹长度为3 m,运动速度接近0.5 m/s,自远处逐渐接近磁力计探头,至推车前缘稍稍越过磁力计上方为止;之后将方向倒转,推车自南向北重复该过程。

得到实验结果如图13所示。

由于实验进行中无法准确控制推车前进速度,而且每次抵达终点的位置有微小差距,因此实验数据在时间轴上并不完全同步,其磁特征信号波形在横轴上方向有些许位移。

从地磁场总场的变化来看,其初始的标量值为49.91 μT,目标自不同方向接近过程中其变化在47.9～55.75 μT之间。目标与北南方向夹角为0°时,其磁特征信号幅值一直在减少;夹角为180°时,磁特征信号先升后降。

西向分量较为微弱,变化范围在1.567～5.563 μT之间,其变化规律是夹角为180°时磁特征信号先升高后降低,为0°时则持续降低;北向分量变化范围在23.769～29.995 μT之间;天向分量变化范围最大,在-48.016～-41.303 μT之间,夹角为0°时有一个较小的减小过程,之后增大;夹角为180°时幅值磁特征信号持续增大。

图13 各敏感轴分量及总场磁感应强度

由以上实验数据通过以下步骤可计算出当地地磁场的方向。

1) 求取3个敏感轴方向及磁总场量在初始状态下的磁感应强度,求得其均值如表3所示。

表3 三轴及总场初始磁感应强度

2) 计算磁倾角α和磁偏角β

(8)

(9)

由上述计算结果可知,实验场地的背景地磁场与北京地区总体地磁场基本一致,磁倾角为向下53.13°,磁偏角为北偏西7.913°。

由图13可知,当推车目标接近磁力计探头1 m范围时,其三轴磁场分量及总场磁感应强度发生明显变化。由于实验过程中难以控制目标速度始终如一,因此其发生明显变化的距离与时间点难以确定。

为对目标磁特征信号变化原因进行更加精准的分析,现对整个实验过程进行仿真复现。因进行三维瞬态求解所需计算量过大,结合上文所述目标静态磁场与动态磁场高度一致,因此采用静磁场仿真求解器,在原运动轨迹上设置一系列观测点用以采集目标距离磁力计探头不同距离下的磁特征信号,如图14所示。

图14 实验仿真设置

地磁场设置为磁感应强度49.89 μT,方向为沿Z轴负方向60°(磁倾角60°),沿X轴正方向7°(磁偏角7°),目标简化为5块30 mm×30 mm×0.5 mm的45钢铁板,铁板上下间隔5 mm,前后间隔200 mm,观测点位于铁板中线处,以前部铁板前沿下方150 mm处为原点,以100 mm为间隔,向前延伸至1 m处。以Y轴正向为正南方向,变换目标朝向为0°和180°,得到其仿真结果如图15所示。