CN103559376B - 基于等效磁矩的潜航平台磁异特征场表征方法 - Google Patents
基于等效磁矩的潜航平台磁异特征场表征方法 Download PDFInfo
- Publication number
- CN103559376B CN103559376B CN201310602015.0A CN201310602015A CN103559376B CN 103559376 B CN103559376 B CN 103559376B CN 201310602015 A CN201310602015 A CN 201310602015A CN 103559376 B CN103559376 B CN 103559376B
- Authority
- CN
- China
- Prior art keywords
- magnetic
- platform
- field
- submariner
- magnetic moment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
本发明涉及水下平台隐身设计领域,公开了一种基于等效磁矩的潜航平台磁异特征场表征方法,包括以下步骤:(1)联合边界条件下,对潜航平台磁化场数值的建模;(2)潜航平台磁异特征场分布表征:(3)基于场特征的潜航平台磁矩反演:(4)对潜航平台磁矩的磁偶极子进行模拟分析:(5)根据等效磁矩磁偶极子磁场表征潜艇磁异特征场。本发明建立了基于磁异特征场的等效磁矩反演方法,将复杂的磁异特征场分析问题归结为等效磁矩表征问题,解决了复杂因素导致的分析规模受限问题。适用于水下潜航平台的磁隐身设计,也可用于抗磁异探测的磁诱饵模拟。
Description
技术领域
本发明涉及属于水下平台隐身设计领域,尤其涉及一种基于等效磁矩的潜航平台磁异特征场表征方法。
背景技术
随着电磁侦察技术的发展,磁异探测技术将成为对潜侦察的主要手段之一。其基本工作原理是:潜航平台壳体、设备、电缆等铁磁性结构受到地磁场磁化而产生附加静磁场,使地磁场的分布规律产生畸变,形成磁异特征场,通过被动非接触测量,从磁异特征场空间分布曲线中提取相对于地磁场的特征信息,实现对目标的定位和跟踪,磁异探测技术是侦察、监视以及获取潜航目标信息的最可靠手段之一,具有隐蔽性能好、抗干扰能力强、定位精度高等特点。为对抗磁异探测,需要开展水下潜航平台的磁隐身设计,也需要发展抗磁异探测的磁诱饵模拟技术。
磁异特征场表征是潜航平台的磁隐身设计以及磁诱饵模拟的基础,通过获取磁异特征场的发生特性、分布规律,支持磁隐身设计,也可帮助建立抗磁异探测的磁诱饵磁场模拟方法。潜航平台磁异特征场关联因素繁多,主要包括地磁场等磁场源的磁化后产生的感应磁场和磁滞效应引起的固定磁场,由于固定磁场的理论建模与目标的磁化历史直接相关,且平台外形较为复杂、磁化不均匀,直接用解析方法会产生较大的近似误差,无法准确表征磁异特征场机理。
发明内容
本发明要解决的技术问题在于提供一种潜航平台磁异特征场表征方法,基于矢量叠加的磁矩反演和磁偶极子模拟得到潜航平台磁异特征场分布特性。
本发明解决其技术问题所采用的技术方案是:一种基于等效磁矩的潜航平台磁异特征场表征方法,包括以下步骤:
(1)联合边界条件下,对潜航平台磁化场数值的建模;
(2)潜航平台磁异特征场分布表征:
(2.1)对潜航平台被磁化后产生沿磁化方向的磁矩进行矢量分析,得到总磁感应强度B或总磁场强度H的幅度,即|B|或|H|;
(2.2)计算磁场B或H的空间变化率G,
其中,Bxx+Byy+Bzz=0;Bxy-Byx=0,Bxz-Bzx=0,Byz-Bzy=0,Bxy=Byx,Bxz=Bzx,Byz=Bzy,traceG=Bxx+Byy+Bzz=0;
(3)基于场特征的潜航平台磁矩反演:
根据独立点源磁偶极子磁场表达式,获得磁偶极子磁矩:
(4)对潜航平台磁矩的磁偶极子进行模拟分析:
用独立点源磁偶极子模型表示潜航平台的磁场,
其中,为磁矩矢量,为源点指向场点的单位矢量;
(5)根据等效磁矩磁偶极子磁场表征潜艇磁异特征场
其中,Hp为潜艇等效磁矩磁偶极子磁场,Hs为潜艇产生的磁异特征场。
按上述方案,步骤(1)中对潜航平台磁化场数值的建模采用以下步骤:
(1.1)基于有限元方法,对平台结构和不均匀介质进行离散,将体积划分成许多小体积四面体单元,建立插值函数,近似表达每一单元的未知函数;
(1.2)采用极值变分公式,建立磁场等效求解模型;
(1.3)针对等效求解模型建立多未知数的线性方程组,求解各离散点的磁矢位;
(1.4)对感应磁场和固定磁场分别进行求解,矢量叠加后获得磁异特征场。
有限元数值方法静磁求解方法通用性好,计算精度高。
按上述方案,为提高步骤(3)中磁矩的求解精度,采用以下步骤计算最佳磁矩:
(3.1)设有N个点的H值,求出每个点对应的磁矩Mi,最佳磁矩M0满足的准则是:
当磁矩满足下式时:
为最优解;
(3.2)将式(e)代入式(f),求得最优磁矩为
(3.3)将式(g)代入式(b),得到磁偶极子磁矩的表达式:
本发明产生的有益效果是:
(1)建立了基于磁通等效原理的固定磁场理论模型和等效环境场磁化作用下的感应磁场模型,可直观表征潜艇磁异特征场的近、远场分布和梯度张量规律。
(2)建立了基于磁异特征场的等效磁矩反演方法,将复杂的磁异特征场分析问题归结为等效磁矩表征问题,解决了复杂因素导致的分析规模受限问题。适用于水下潜航平台的磁隐身设计,也可用于抗磁异探测的磁诱饵模拟。
(3)有限元数值静磁求解方法具有通用性好,计算精度高、适于计算不规则结构和非均匀介质等优点;
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明实施例的方法流程图;
图2是联合边界条件下磁化场有限元数值计算与文献、实验曲线对比图;
图3是实施例中等效磁矩表征曲线与直接数值计算曲线对比图;
图4是实施例中潜艇等效磁矩相似性量化曲线。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
如图1所示,基于等效磁矩的潜航平台磁异特征场表征方法,其特征在于:(1)联合边界条件下的潜航平台磁化场数值建模
基于有限元方法,对平台结构和不均匀介质进行离散,将体积V划分成许多小体积四面体单元,建立插值函数,近似表达每一单元的未知函数;采用极值变分公式,建立磁场等效求解模型;针对等效求解模型建立多未知数的线性方程组,求解各离散点的磁矢位;对感应磁场和固定磁场分别进行求解,矢量叠加后获得磁异特征场。
(2)潜航平台磁异特征场分布表征
理论上,铁磁体被磁化后会产生沿磁化方向的磁矩,对磁矩进行矢量分析:
1)总磁感应强度B或总磁场强度H的幅度,即|B|或|H|;
2)磁场B或H的空间变化率。由于梯度为矢量,为表征其梯度信息,必须采用磁场梯度张量表示,即
虽然梯度张量有9个要素,由于在无源空间中B的散度和旋度均为零,因此有
divB=0
rotB=0
或写成
所以Bxx+Byy+Bzz=0;Bxy-Byx=0,Bxz-Bzx=0,Byz-Bzy=0,由此在9个要素中,只有5个要素是独立的。梯度张量是对称的,其迹为零,Bxy=Byx,Bxz=Bzx,Byz=Bzy,traceG=Bxx+Byy+Bzz=0。上述参量反映了磁场及其变化量,磁异探测根据这些参量探测潜航平台,超过阈值后,即认为发现目标。这些参量的值越大,被发现的概率越大。
(3)基于场特征的潜航平台磁矩反演
通过独立点源磁偶极子磁场表达式,可得磁偶极子磁矩:
在等效磁矩求解时,已知各点的Hx,Hy,Hz,求Mx,My,Mz。由(2)可知,仅需一个点的H值,即可求出磁矩。但对于不规则结构的潜航平台而言,为提高磁矩的求解精度,须采用多点的磁场数据,本研究采用最小二乘法开展潜艇等效磁矩反演理论研究。
设有N个点的H值,求出每个点对应的磁矩Mi,最佳磁矩M0满足的准则是
磁矩满足下式时
为最优解。将式(3)代入根据式(4),可求得最优磁矩为
代入(2)式,有
(4)潜航平台磁矩的磁偶极子模拟分析
根据级数展开原理,对于复杂结构的磁性目标,其在周围近场空间的磁场不能用独立点源磁偶极子模型精确表示。磁性物体的磁场必须用分布式多极子模型表示:
其中,为磁矩矢量,为源点指向场点的单位矢量。随距离增大,高阶次极子的作用变得很小,可以忽略。模型变成了独立点源偶极子模型:
(5)根据等效磁矩磁偶极子磁场表征潜艇磁异特征场
为定量分析潜艇等效磁矩的适用条件,本发明提出了等效磁矩磁偶极子磁场与潜艇磁异场相似性系数概念,相似性系数定义为
其中,Hp为潜艇等效磁矩磁偶极子磁场,Hs为潜艇产生的磁异特征场。通过相似性系数可以根据等效磁矩磁偶极子磁场表征潜艇磁异特征场
本发明方法有效性的验证:
(1)根据前述磁场模型和数值求解的计算方法,开发了磁场数值计算平台,同IEEE Trans.on Magnetics中的文献数据作了典型圆柱目标的对比验证。对比曲线见图2。
(2)构建潜艇的仿真模型,艇长为60m,直径为10m,壳体厚度为0.1m。采用上文所述的数值有限元方法计算潜艇的磁异特征场,地磁场的磁场强度取为80A/m,固定磁场取为20A/m。为简化问题,设固定磁场与地磁方向同向。潜艇内部设置一电流环,电流为10A,模拟艇内设备的电流。同时考虑到航空磁异探测的实际情况,潜艇磁异场的计算面为潜艇正上方(即XOY面)200m处,纵向(即X方向)的计算区域从-300m到300m,横向(即Y方向)的计算区域从-300m到300m。
(3)将前述潜艇模型作为研究对象,潜艇参数不变,激励参数取为80A/m,沿潜艇长度方向(X向)。采用远场的磁场仿真计算数据作为计算潜艇等效磁矩的依据,利用磁偶极子磁矩计算方法,计算出潜艇的等效磁矩约为1.08×105A·m2,等效磁矩沿磁化方向(X向)。将等效磁矩用磁偶极子模拟,磁场发射模型计算距离其300m处的特征面磁场,再采用前述求解方法计算潜艇在同样位置的磁场,比较二者的磁场分布相似性,如图3所示。
(4)为定量衡量场分布的相似性,计算等效磁偶极子磁场与潜艇磁异特征场的相似性系数。定义为
Hp与Hs分别为磁偶极子和潜艇产生的磁场。S包含三个特征面,S=SX+SY+SZ SX,SY,SZ分别为X,Y,Z特征面的面积。Hp与Hs完全相同时,ρs=1。ρsx,ρsy,ρsz分别为Hx,Hy,Hz的相似性系数,总相似性系数为
沿某一截线l的相似性定义为
图4为潜艇等效磁矩磁偶极子与潜艇磁异特征场的相似性系数随b/L(距离/艇长)的变化曲线。图中可看出,b/L>5时,相似性优于0.9,显示出良好的相似性。
(5)根据前述的计算结果,前述潜艇模型磁矩约为105A·m2量级。可以采用多匝线圈等方式模拟,对于磁诱饵上的小型磁偶极子来说,需要对线圈横截面积S、线圈电流I与匝数N进行综合优化设计。因此,潜艇抗磁异探测的磁诱饵技术从原理上是可行的。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (2)
1.一种基于等效磁矩的潜航平台磁异特征场表征方法,其特征在于,
包括以下步骤:
(1)联合边界条件下,对潜航平台磁化场数值的建模;
(1.1)基于有限元方法,对平台结构和不均匀介质进行离散,将潜航平台按体积划分成许多小体积四面体单元;
(1.2)建立插值函数,近似表达每一单元的未知函数;
(1.3)采用极值变分公式,建立磁场等效求解模型;
(1.4)针对等效求解模型建立多未知数的线性方程组,求解各离散点的磁矢位;
(1.5)对感应磁场和固定磁场分别进行求解,矢量叠加后获得磁异特征场;
(2)潜航平台磁异特征场分布表征:
(2.1)对潜航平台被磁化后产生沿磁化方向的磁矩进行矢量分析,得到总磁感应强度B或总磁场强度H的幅度,即|B|或|H|;
(2.2)计算磁场B或H的空间变化率G,
其中,Bxx+Byy+Bzz=0;Bxy-Byx=0,Bxz-Bzx=0,Byz-Bzy=0,Bxy=Byx,Bxz=Bzx,Byz=Bzy,traceG=Bxx+Byy+Bzz=0;
(3)基于场特征的潜航平台磁矩反演:
根据独立点源磁偶极子磁场表达式,获得磁偶极子磁矩:
(4)对潜航平台磁矩的磁偶极子进行模拟分析:
用独立点源磁偶极子模型表示潜航平台的磁场,
其中,为磁矩矢量,为源点指向场点的单位矢量;
(5)根据等效磁矩磁偶极子磁场表征潜艇磁异特征场
其中,Hp为潜艇等效磁矩磁偶极子磁场,Hs为潜艇产生的磁异特征场。
2.根据权利要求1所述的潜航平台磁异特征场表征方法,其特征在于,步骤(3)中磁矩采用以下步骤计算:
(3.1)设有N个点的H值,求出每个点对应的磁矩Mi,最佳磁矩M0满足的准则是:
当磁矩满足下式时:
为最优解;
(3.2)将式(e)代入式(f),求得最优磁矩为
(3.3)将式(g)代入式(b),得到磁偶极子磁矩的表达式:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310602015.0A CN103559376B (zh) | 2013-11-25 | 2013-11-25 | 基于等效磁矩的潜航平台磁异特征场表征方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310602015.0A CN103559376B (zh) | 2013-11-25 | 2013-11-25 | 基于等效磁矩的潜航平台磁异特征场表征方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103559376A CN103559376A (zh) | 2014-02-05 |
CN103559376B true CN103559376B (zh) | 2016-08-17 |
Family
ID=50013622
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310602015.0A Active CN103559376B (zh) | 2013-11-25 | 2013-11-25 | 基于等效磁矩的潜航平台磁异特征场表征方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103559376B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108072906B (zh) * | 2016-11-18 | 2020-01-10 | 北京自动化控制设备研究所 | 一种分布式磁探测磁目标识别方法 |
CN107122518B (zh) * | 2017-03-23 | 2020-10-30 | 中国人民解放军海军工程大学 | 磁场相似性度量方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5831873A (en) * | 1997-01-27 | 1998-11-03 | Kohnen; Kirk K. | Magnetic dipole target classifier and method |
US6539327B1 (en) * | 1997-09-11 | 2003-03-25 | Commissariat A L'energie Atomique | Process for determining the position of a moving object using magnetic gradientmetric measurements |
CN102928794A (zh) * | 2012-10-22 | 2013-02-13 | 中国船舶重工集团公司第七一〇研究所 | 一种基于三轴矢量磁传感器阵的磁性目标磁矩测量方法 |
-
2013
- 2013-11-25 CN CN201310602015.0A patent/CN103559376B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5831873A (en) * | 1997-01-27 | 1998-11-03 | Kohnen; Kirk K. | Magnetic dipole target classifier and method |
US6539327B1 (en) * | 1997-09-11 | 2003-03-25 | Commissariat A L'energie Atomique | Process for determining the position of a moving object using magnetic gradientmetric measurements |
CN102928794A (zh) * | 2012-10-22 | 2013-02-13 | 中国船舶重工集团公司第七一〇研究所 | 一种基于三轴矢量磁传感器阵的磁性目标磁矩测量方法 |
Non-Patent Citations (5)
Title |
---|
吴楠等.舰船短波电磁环境仿真建模的工程应用研究.《舰船科学技术》.2010,第32卷(第03期),全文. * |
基于磁偶极子等效的潜艇空间磁场分布;张朝阳等;《舰船科学技术》;20130115;第35卷(第1期);全文 * |
张琦等.潜艇磁化场的有限元方法研究.《舰船科学技术》.2009,第31卷(第1期),全文. * |
磁性物体磁偶极子模型适用性的试验研究;张朝阳等;《应用基础与工程科学学报》;20101015;第18卷(第5期);正文第862页第1段第1-3行,第863页第1-3段、倒数第1段,第864页第3节第1-3段 * |
黄玉.地磁场测量及水下定位技术研究.《中国博士学位论文全文数据库 工程科技Ⅱ辑》.2012,(第05期),第21页第2.2节第1-2段,第22页第1段. * |
Also Published As
Publication number | Publication date |
---|---|
CN103559376A (zh) | 2014-02-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Holmes | Modeling a ship’s ferromagnetic signatures | |
Yin et al. | 3D time-domain airborne EM forward modeling with topography | |
US20130018588A1 (en) | Method of real time subsurface imaging using gravity and/or magnetic data measured from a moving platform | |
CN108197389A (zh) | 二维强磁性体磁场的快速、高精度数值模拟方法 | |
Zhou et al. | Spatial signature analysis of submarine magnetic anomaly at low altitude | |
CN109254327B (zh) | 三维强磁性体的勘探方法及勘探系统 | |
Desjardins et al. | Analytical modeling of the transient response of a coil encircling a ferromagnetic conducting rod in pulsed eddy current testing | |
Jian-jun et al. | A method for real-time compensation of moving ferromagnet's magnetic moment | |
CN114236639A (zh) | 基于核磁共振与电阻率联合的矿井灾害水源动态监测方法 | |
Guo et al. | 3D correlation imaging of magnetic total field anomaly and its vertical gradient | |
CN103559376B (zh) | 基于等效磁矩的潜航平台磁异特征场表征方法 | |
Choi et al. | Optimization of degaussing coil currents for magnetic silencing of a ship taking the ferromagnetic hull effect into account | |
US20140056481A1 (en) | Methods of multinary inversion for imaging objects with discrete physical properties | |
Um | Three-dimensional finite-element time-domain modeling of the marine controlled-source electromagnetic method | |
Zhang et al. | Automatic search algorithms for near-field ferromagnetic targets based on magnetic anomaly detection | |
Chun | Method of moving frames to solve time-dependent Maxwell's equations on anisotropic curved surfaces: Applications to invisible cloak and ELF propagation | |
Huang et al. | 3D full-time anisotropic TEM modelling using a mixed BDF2/SAI method | |
Kim et al. | Efficient re-degaussing technique for a naval ship undergoing a breakdown in degaussing coils | |
Zuo et al. | Multi-objective optimization design method of naval vessels degaussing coils | |
Clark | Interpretation of the magnetic gradient tensor and normalized source strength applied to the Tallawang magnetite skarn deposit, New South Wales, Australia | |
BUCHAU | Precise and robust magnetic field computations for high-end smart sensor applications | |
Zhou et al. | Approach for a fast computation of induced magnetic field created by naval vessels | |
Yin et al. | Three-dimensional reconstruction of a small-scale magnetic target from magnetic gradient observations | |
Vijn et al. | Prediction of the Magnetic State of Ferromagnetic Objects by Assimilating Data into a Physical Model | |
YAN et al. | Nonlinear Inversion with Quadratic Function Approaching Method for Magnetotelluric Data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |