CN107748813A - 非晶丝在非轴向磁场作用下的巨磁阻抗建模方法 - Google Patents

非晶丝在非轴向磁场作用下的巨磁阻抗建模方法 Download PDF

Info

Publication number
CN107748813A
CN107748813A CN201710890013.4A CN201710890013A CN107748813A CN 107748813 A CN107748813 A CN 107748813A CN 201710890013 A CN201710890013 A CN 201710890013A CN 107748813 A CN107748813 A CN 107748813A
Authority
CN
China
Prior art keywords
amorphous wire
magnetic
equation
impedance
magnetic field
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.)
Granted
Application number
CN201710890013.4A
Other languages
English (en)
Other versions
CN107748813B (zh
Inventor
许洪光
路小菲
胡琴芳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Graduate School Harbin Institute of Technology
Original Assignee
Shenzhen Graduate School Harbin Institute of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenzhen Graduate School Harbin Institute of Technology filed Critical Shenzhen Graduate School Harbin Institute of Technology
Priority to CN201710890013.4A priority Critical patent/CN107748813B/zh
Publication of CN107748813A publication Critical patent/CN107748813A/zh
Application granted granted Critical
Publication of CN107748813B publication Critical patent/CN107748813B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/12Measuring magnetic properties of articles or specimens of solids or fluids
    • G01R33/1223Measuring permeability, i.e. permeameters

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Measuring Magnetic Variables (AREA)
  • Hall/Mr Elements (AREA)

Abstract

本发明提出了一种非晶丝在非轴向磁场作用下的巨磁阻抗建模方法,首先利用LLG方程,对磁化矢量旋转方程进行约束,建立非晶丝内部总能量的表达式;然后,根据体系自由能最小原理,建立磁化矢量运动方程,通过对运动方程求解得到非轴向磁场作用下非晶丝磁导率张量,利用磁场强度和磁感应强度的本征关系以及麦克斯韦方程,建立关于非晶丝环向和轴向磁场参量的电磁场方程;最后针对高频和低频的应用条件,对非晶丝电磁场方程的解进行简化,使用高频近似法和低频近似法得到了非晶丝巨磁阻抗的一般表达式。本发明能精确建立非晶丝在非轴向磁场下的巨磁阻抗传感模型,给出巨磁阻抗的计算方法。

Description

非晶丝在非轴向磁场作用下的巨磁阻抗建模方法
技术领域
本发明属于磁性材料传感技术领域,具体涉及一种非轴向磁场作用下的巨磁阻抗建模方法。
背景技术
非晶丝磁导率高,矫顽力较小,具有显著的磁阻抗效应以及低场下的高灵敏度,和现有的一些半导体器件、磁电阻器件、磁通门比较起来,非晶丝磁场传感器兼具稳定性、高灵敏度、微型化和低功耗的优势。非晶丝传感器在工业、导航、军事、安全、无损检测、地磁及生物医学等领域具有广泛应用,例如,高密度磁记录设备、微磁设备、军用传感器、交通监控、卫星飞行器检测、工业设备控制、地磁传感器等都应用了非晶丝传感器。
非晶丝的阻抗张量包含直流阻抗部分以及交流阻抗部分,直流阻抗可以从经典电动力模型中得到,但为了得到交流阻抗,除了分析静态磁化外还需要考虑高频交流电激发的交变磁场引发的交流磁化。磁矩的电磁场需要满足Maxwell方程,磁化转动需要满足Landau-Lifshitz方程。由于非轴向磁场下非晶丝磁化过程相当复杂,磁矩方程求解通常难以进行,并且几乎不可能得到解析解。只有在某些限定条件下,比如忽略退磁场,交换作用以及内部应力,才可以对方程进行近似线性化,此时可以将磁导率可看作仅由磁矩的旋转确定。
发明内容
为解决现有技术中存在的问题,本发明提出了一种非晶丝在非轴向磁场作用下的巨磁阻抗建模方法,给出了非轴向磁场作用下非晶丝巨磁阻抗的计算方法。根据经典电动力模型,非晶丝的轴向巨磁阻抗正比于环向磁导率的开方,环向磁化过程主要是通过直流外磁场对环向磁导率的变化来体现,在相对较高频率下,使用磁化矢量旋转理论计算环向磁导率。
本发明具体通过如下技术方案实现:
一种非晶丝在非轴向磁场作用下的巨磁阻抗建模方法,所述方法用于非轴向磁场作用下非晶丝巨磁阻抗的计算,所述方法包括:首先利用LLG方程,对磁化矢量旋转方程进行约束,建立非晶丝内部总能量的表达式;根据体系自由能最小原理,建立磁化矢量运动方程,通过对运动方程求解得到非轴向磁场作用下非晶丝磁导率张量;然后利用磁场强度和磁感应强度的本征关系以及麦克斯韦方程,建立关于非晶丝环向和轴向磁场参量的电磁场方程;最后使用高频近似法和低频近似法得到非晶丝巨磁阻抗表达式。
本发明的有益效果是:本发明提供了一种对非晶丝在非轴向外磁场下的巨磁阻抗作用机理的建模方法,使用该方法计算的结果可以分析非晶丝在非轴向外磁场下的巨磁阻抗特性。本发明能够得到非晶丝最佳的工作环境,在特定的外界参数范围内可以制备灵敏度高、低功耗、微型化的磁敏传感器。本发明可以极大提高非晶丝磁场传感器的检测精度,强化其对空间方向磁场的分辨能力,扩大巨磁阻抗的应用范围。可以用于高密度磁记录设备、工业设备控制、交通测量和监控、电子磁罗盘、微磁医疗设备等方面。
附图说明
图1是非晶丝工作原理图;
图2是磁化强度与外磁场以及等效各向异性场的角度示意图。
具体实施方式
下面结合附图说明及具体实施方式对本发明进一步说明。
本发明结合人体工程学与体育健康测量等先验知识提出全新的特征点标记方法,标记结果具有更高的准确性。特征点可根据实际需求增加或减少特征点数目与位置。具体的各个部位特征点定义规则及标定方法如下。
本发明的方法采用的材料是接近于零的负向磁致伸缩系数的Co基非晶丝。Co基非晶丝在无外磁场施加时具有芯-壳磁畴结构,内芯具有轴向各向异性,外壳是交替反向的圆周各向异性。
非晶丝工作原理如图1所示,其工作条件是,直流电流源和高频交流电流源驱动非晶丝,在非晶丝内部和表面产生环向直流和交流磁场,直流电流源驱动两轴亥姆霍兹线圈产生任意方向均匀磁场。其中,Hex是交流磁场强度,Hdc是环向直流磁场强度,Idc是直流电流源,r是电阻,Iac是环向电流,ew是两轴亥姆霍兹线圈产生的电场。
非晶丝的阻抗是一个对称张量形式,分为对角分量和非对角分量,可以用非晶丝表面的电磁场信息表示:
其中,Zzz是对角阻抗分量,是非对角阻抗分量,c是光速,N表示耦合线圈的匝数,l是导线的长度ez是轴向电场强度,是轴向电场强度,是轴向磁场强度
非晶丝复阻抗的经典电动力模型,当驱动交流电流的频率为几MHz以上,阻抗近似为:
其中,a是非晶丝的半径,ω是激励交流电流频率,ρ是非晶丝的电阻率,是环向磁导率。
在相对较高频率下,使用磁化矢量旋转理论计算环向磁导率。在经典电动力模型中,是在假设磁感应强度和磁场强度是线性变化的情况下得到的,实际情况下非晶丝在受到非轴向外磁场的磁化时,磁感应强度和磁场强度的关系不是线性的,这时我们需要增加适合的磁化矢量旋转方程加以约束。使用朗道-栗弗席兹-吉尔伯特(Landau-Lifshitz-Gilber,LLG)方程对磁化矢量进行约束,其表达式为:
其中,M是磁化强度,Heff是磁矩感受到的总的有效磁场,Ms是饱和磁化强度,γ是旋磁比、γ=|e|/me,α是定性的无量纲的阻尼系数。
由于非轴向磁场下非晶丝磁化过程相当复杂,磁矩方程求解通常难以进行,并且几乎不可能得到解析解。只有在某些限定条件下,比如忽略退磁场,交换作用以及内部应力,才可以对方程进行近似线性化,此时可以将磁导率可看作仅由磁矩的旋转确定。
体系自由能的能量项主要包括磁各向异性能、外磁场能、环向磁场能、磁应力能、退磁场能和畴壁交换能等,忽略退磁场和交换作用场等的影响,只考虑有各向异性场的情况,如图2,非晶丝内的总能量可以表示为:
其中,μ0是真空磁导率,Hk是各向异性场的等效磁场强度,Ms是饱和磁化强度,ψ是各向异性轴与横截面方向的夹角,θex是外磁场与非晶丝轴线的夹角,θ是饱和磁化磁矩方向与圆周方向的夹角。
根据体系自由能最小原理,可建立LLG运动方程。对LLG方程求解,得到磁导率张量。
其中
ωm=γ4πM
ω1=γ[Ha cos2(θ-ψ)+He sinθ]
ω2=γ[Ha cos2{2(θ-ψ)}+Hesinθ]
下标r表示非晶丝的径向,下标表示非晶丝的环向,下标z表示非晶丝的轴向,θ是饱和磁化磁矩方向与圆周方向的夹角,表示环向-轴向磁导率,其他的以此类推。M是饱和磁化强度,Ha是各项异性磁场的强度,γ是回旋常数,κ是吉尔伯特阻尼参数。
利用磁场强度和磁感应强度的本征关系,以及麦克斯韦方程,建立关于非晶丝环向和轴向磁场参量的电磁场方程。
其中,,r表示圆柱坐标系的径向坐标,σ表示电导率,ω是激励交流电流频率。
磁场分量的解为
其中,A,B,C,D是系数,可以通过边界条件求得,这里特征根:
a是非晶丝的半径。
非晶丝的表面阻抗张量可以由电磁场方程与磁导率张量联立得到,高频近似法下的复阻抗,对角阻抗分量Zzz和非对角阻抗分量分别为:
在低频率情况下,解可以表示为以下序列的线性组合:
hz=AN1+BN2+CN3+DN4
这里函数MiNj表示磁场分量解的贝塞尔函数的泰勒展开,由下面的式子决定:
其中的系数可以通过递归方式计算,
整理可得得到低频情况下的对角阻抗分量Zzz和非对角阻抗分量
本发明提供的非晶丝在非轴向磁场下的巨磁阻抗的建模方法,使用该方法计算的结果可以分析非晶丝在非轴向外磁场下的巨磁阻抗特性。本发明能够得到非晶丝最佳的工作环境,在特定的外界参数范围内可以制备灵敏度高、低功耗、微型化的磁敏传感器。本发明可以极大提高非晶丝磁场传感器的检测精度,强化其对空间方向磁场的分辨能力,扩大巨磁阻抗的应用范围。本发明还可以用于高密度磁记录设备、工业设备控制、交通测量和监控、电子磁罗盘、微磁医疗设备等方面。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。

Claims (9)

1.一种非晶丝在非轴向磁场作用下的巨磁阻抗建模方法,所述方法用于非轴向磁场作用下非晶丝巨磁阻抗的计算,所述非晶丝的巨磁阻抗是指使用高频交流电流源驱动非晶丝时,非晶丝两端的交流阻抗;其特征在于:所述方法包括:S1、利用LLG方程,对磁化矢量旋转方程进行约束,建立非晶丝内部总能量的表达式;S2、根据体系自由能最小原理,建立磁化矢量运动方程,通过对运动方程求解得到非轴向磁场作用下非晶丝磁导率张量;S3、利用磁场强度和磁感应强度的本征关系以及麦克斯韦方程,建立关于非晶丝环向和轴向磁场参量的电磁场方程;S4、使用高频近似法和低频近似法得到非晶丝巨磁阻抗表达式。
2.根据权利要求1所述的方法,其特征在于:所述非晶丝是有接近于零的负向磁致伸缩系数的Co基非晶丝。
3.根据权利要求1所述的方法,其特征在于:使用两轴亥姆霍兹线圈,直流电流源驱动两轴亥姆霍兹线圈在所述非晶丝周围产生二维任意方向均匀磁场。
4.根据权利要求1所述的方法,其特征在于:所述非晶丝的阻抗是一个对称张量形式,分为对角分量和非对角分量,用非晶丝表面的电磁场信息表示:
其中,Zzz是对角阻抗分量,是非对角阻抗分量,c是光速,N表示线圈的匝数,l是导线的长度ez是轴向电场强度,是轴向电场强度,是轴向磁场强度。
5.根据权利要求1所述的方法,其特征在于:使用LLG方程对磁化矢量进行约束,其表达式为:
其中,M是磁化强度,Heff是磁矩感受到的总的有效磁场,Ms是饱和磁化强度,γ是旋磁比,α是定性的无量纲的阻尼系数。
6.根据权利要求1所述的方法,其特征在于:只考虑有各向异性场的情况,非晶丝内的总能量可以表示为:
其中,μ0是真空磁导率,Hk是各向异性场的等效磁场强度,Ms是饱和磁化强度,ψ是各向异性轴与横截面方向的夹角,θex是外磁场与非晶丝轴线的夹角,θ是饱和磁化磁矩方向与圆周方向的夹角。
7.根据权利要求1所述的方法,其特征在于,所述步骤S2包括对运动方程求解,得到磁导率张量:
其中,
ωm=γ4πM
ω1=γ[Hacos2(θ-ψ)+Hesinθ]
ω2=γ[Hacos2{2(θ-ψ)}+Hesinθ]
下标r表示非晶丝的径向,下标表示非晶丝的环向,下标z表示非晶丝的轴向,θ是饱和磁化磁矩方向与圆周方向的夹角,表示环向-轴向磁导率,其他的以此类推。M是饱和磁化强度,Ha是各项异性磁场的强度,γ是回旋常数,κ是吉尔伯特阻尼参数。
8.根据权利要求7所述的方法,其特征在于,所述步骤S3包括:利用磁场强度和磁感应强度的本征关系,以及麦克斯韦方程,建立关于非晶丝环向和轴向磁场参量的电磁场方程:
其中,r表示圆柱坐标系的径向坐标,ω是激励交流电流频率,σ是非晶丝的电导率,磁场分量的解为
其中,A,B,C,D是系数,可以通过边界条件求得。
这里特征根:ω是激励交流电流频率。
9.根据权利要求8所述的方法,其特征在于:非晶丝的表面阻抗张量由电磁场方程与磁导率张量联立得到,高频近似法下的复阻抗,对角阻抗分量Zzz和非对角阻抗分量分别为:
在低频率情况下,可以表示为以下序列的线性组合:
hz=AN1+BN2+CN3+DN4
CN201710890013.4A 2017-09-27 2017-09-27 非晶丝在非轴向磁场作用下的巨磁阻抗建模方法 Active CN107748813B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710890013.4A CN107748813B (zh) 2017-09-27 2017-09-27 非晶丝在非轴向磁场作用下的巨磁阻抗建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710890013.4A CN107748813B (zh) 2017-09-27 2017-09-27 非晶丝在非轴向磁场作用下的巨磁阻抗建模方法

Publications (2)

Publication Number Publication Date
CN107748813A true CN107748813A (zh) 2018-03-02
CN107748813B CN107748813B (zh) 2021-01-01

Family

ID=61255916

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710890013.4A Active CN107748813B (zh) 2017-09-27 2017-09-27 非晶丝在非轴向磁场作用下的巨磁阻抗建模方法

Country Status (1)

Country Link
CN (1) CN107748813B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112380737A (zh) * 2020-09-02 2021-02-19 南京理工大学 一种基于表面阻抗边界的薄层电磁结构的时域分析方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104532174A (zh) * 2014-12-31 2015-04-22 哈尔滨工业大学 一种非晶微丝可调控阻抗线性响应量程的方法
CN105891743A (zh) * 2014-12-19 2016-08-24 中国人民解放军军械工程学院 一种巨磁阻抗传感器的激励方法及装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105891743A (zh) * 2014-12-19 2016-08-24 中国人民解放军军械工程学院 一种巨磁阻抗传感器的激励方法及装置
CN104532174A (zh) * 2014-12-31 2015-04-22 哈尔滨工业大学 一种非晶微丝可调控阻抗线性响应量程的方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
L.V.PANINA等: "Giant magneto-impedance (GMI) in amorphous wire, single layer film and sandwich film", 《PHYSICA A:STATISTICAL MECHANICS AND ITS APPLICATIONS》 *
V.RAPOSO等: "Internal Magnetoimpedance of Amorphous Wires", 《IEEE TRANSACTIONS ON MAGNETICS》 *
姜婷婷: "非晶丝各向异性对巨磁阻抗的影响研究", 《中国优秀硕士论文全文数据库 工程科技II辑(月刊)电力工程》 *
柴秀丽等: "带材非对角非对称巨磁阻抗效应传感机理建模", 《仪器仪表学报》 *
胡含凯等: "基于非晶丝巨磁阻抗效应的微磁场探测系统", 《电子设计工程》 *
高科军: "基于巨磁阻抗效应的地磁导航研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑(月刊)航空航天科学与工程》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112380737A (zh) * 2020-09-02 2021-02-19 南京理工大学 一种基于表面阻抗边界的薄层电磁结构的时域分析方法

Also Published As

Publication number Publication date
CN107748813B (zh) 2021-01-01

Similar Documents

Publication Publication Date Title
Ripka et al. Advances in magnetic field sensors
CN107315150A (zh) 一种正交基模磁通门传感器
Lu et al. Responsivity and noise of a wire-bonded CMOS micro-fluxgate sensor
Jordán et al. Detection of magnetic moment in thin films with a home-made vibrating sample magnetometer
Panina Asymmetrical giant magneto-impedance (AGMI) in amorphous wires
CN207181651U (zh) 一种正交基模磁通门传感器
Das et al. Detection of low-frequency magnetic fields down to sub-pT resolution with planar-Hall effect sensors
Hrakova et al. Sensitivity and noise of parallel fluxgate sensor with amorphous wire cores
CN107748813A (zh) 非晶丝在非轴向磁场作用下的巨磁阻抗建模方法
Ye et al. A quantitative model for the sensitivity of untuned voltage output fluxgate sensors
Kennedy et al. Applications of nanoparticle-based fluxgate magnetometers for positioning and location
Beato-López et al. Contactless magnetic nanoparticle detection platform based on non-linear GMI effect
Yao et al. High precision MI sensor with low energy consumption driven by low-frequency Wiegand pulse
Chen et al. Sensitivity model for residence times difference fluxgate magnetometers near zero magnetic field
Shah Vibrating sample magnetometery: analysis and construction
Jeng et al. Enhancement in sensitivity using multiple harmonics for miniature fluxgates
Coillot et al. On-board hybrid magnetometer of NASA CHARM-II rocket: principle, design and performances
CN210864010U (zh) 宽频带感应式磁场传感器
US20120101772A1 (en) Device for the magnetic measurement of the rotation of a magnetised ball and method for measuring the rotation of the ball
CN113932939A (zh) 基于扫场法的铁磁共振测温方法
Sinibaldi et al. NMR Detection at 8.9 mT with a GMR based sensor coupled to a superconducting Nb Flux Transformer
Deak et al. Dynamic calculation of the responsivity of monodomain fluxgate magnetometers
Fodil et al. Model calculation of the magnetic induction generated by magnetic nanoparticles flowing into a microfluidic system: Performance analysis of the detection
Dong et al. Simulation and measurements of transient fields from conductive plates of shielded room for SQUID-based ultralow field magnetic resonance imaging
Weiss et al. Excess magnetic noise in orthogonal fluxgates employing discontinuous excitation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant