CN110221230B - 一种交变电流磁光调制下调制器磁场分析方法 - Google Patents

一种交变电流磁光调制下调制器磁场分析方法 Download PDF

Info

Publication number
CN110221230B
CN110221230B CN201910562101.0A CN201910562101A CN110221230B CN 110221230 B CN110221230 B CN 110221230B CN 201910562101 A CN201910562101 A CN 201910562101A CN 110221230 B CN110221230 B CN 110221230B
Authority
CN
China
Prior art keywords
magnetic field
magneto
modulator
alternating current
internal
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
Application number
CN201910562101.0A
Other languages
English (en)
Other versions
CN110221230A (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.)
Rocket Force University of Engineering of PLA
Original Assignee
Rocket Force University of Engineering of PLA
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 Rocket Force University of Engineering of PLA filed Critical Rocket Force University of Engineering of PLA
Priority to CN201910562101.0A priority Critical patent/CN110221230B/zh
Publication of CN110221230A publication Critical patent/CN110221230A/zh
Application granted granted Critical
Publication of CN110221230B publication Critical patent/CN110221230B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/02Measuring direction or magnitude of magnetic fields or magnetic flux

Landscapes

  • Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Magnetic Variables (AREA)
  • Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)

Abstract

本发明涉及一种交变电流磁光调制下调制器磁场分析方法,通过分析调制器结构特点,结合麦克斯韦方程组分别建立调制器内部电磁场的基本模型、空心调制器内部磁场模型、无松弛极化磁光材料内部磁场模型,最终构建交变电流磁光调制下含无松弛极化磁光材料的调制器内部磁场分布模型。通过所设置的基本仿真参数,由交变电流磁光调制下调制器内部电磁场的基本模型,求解空心螺线管的磁场构建无松弛极化磁光材料内部磁场模型,最终构建的磁场分布模型反映了介质内外磁场幅值的频率特性,本发明同现有技术比较的最大优越性在于:通过对磁光调制过程中含无松弛介质的调制器内部磁场分布计算,为详细研究磁光调制过程提供了磁场分析方法。

Description

一种交变电流磁光调制下调制器磁场分析方法
技术领域
本发明属于磁光调制技术领域,涉及一种交变电流磁光调制下含无松弛极化磁光材料的调制器磁场分析方法。
背景技术
磁光调制是以法拉第磁致旋光效应为基础的物理现象,当调制信号为交流信号时,能够实现对传递信号的高精度调制,使其在高精度方位测量、工业、医疗等领域应用日益广泛。
传统磁光调制器内部中轴线位置放置一定长度的磁光材料,当给调制器通以交流调制信号时,调制器内部产生磁场,磁光材料磁化,此时一束光通过磁化后的磁光材料时,光的振动方向旋转一个角度θ,这就是法拉第磁致旋光效应。法拉第旋转角计算方法为θ=VBL,其中B为磁光材料内部的磁场,V为磁光材料的维尔德常数。当调制信号为直流信号时,磁场的表达式清晰,但是当调制信号为交流信号时,调制器内部的磁场分布情况复杂。
在本发明以前的现有技术中,关于交流信号调制下调制器内部磁场的精确描述的文献很少,仅有文献《通以交变电流的长直螺线管内部磁场和电场的分布》葛松华.通以交变电流的长直螺线管内部磁场和电场的分布[J].物理与工程,2003,13(6):6-8,以及申请人团队发表的期刊文章《基于麦克斯韦方程的交变电流长螺线管磁场》蔡伟,伍樊成,杨志勇,等.基于麦克斯韦方程的交变电流长螺线管磁场[J].强激光与粒子束,2015,27(12):123201.
,但是此两篇文献仅仅是交流信号调制下调制器内部磁场的概略分析,分析过程中简化了部分因素,此外关于交流信号调制下调制器内部中轴线位置放置的磁光材料内部磁场(近轴区磁场)的分析,目前没有相关文献报道。
发明内容
针对上述现有技术现状,本发明的目的在于:利用麦克斯韦方程,结合调制系统参数,通过精确计算交流信号调制下调制器内部磁场,以及调制器内部中轴线位置放置的无松弛极化磁光材料的内部磁场(近轴区磁场),提供一种可实现交变电流磁光调制下含无松弛极化磁光材料的调制器磁场内部磁场精确表示的方法。
为达到上述目的,本发明提出一种交变电流磁光调制下的调制器磁场分析方法,其特征在于:通过分析调制器结构特点,结合麦克斯韦方程组分别建立交变电流磁光调制下调制器内部电磁场的基本模型、空心调制器内部磁场模型、无松弛极化磁光材料内部磁场模型,最终构建交变电流磁光调制下含无松弛极化磁光材料的调制器内部磁场分布模型,包括以下步骤:
步骤1:根据调制器结构特点,结合麦克斯韦方程组构建交变电流磁光调制下调制器内部电磁场的基本模型
Figure BDA0002108557420000021
Figure BDA0002108557420000022
式中,Bz为磁场的轴向分量,
Figure BDA0002108557420000023
为电场的圆周方向分量,r代表距中轴线的距离,
Figure BDA0002108557420000024
为真空中的波矢,ε0为真空介电常数,μ0为真空磁导率,ω为加载的正弦交变电流信号的频率;
步骤1.1:给调制器的螺线管线圈通以正弦交变电流Is=Ime-iωt,即一个幅值为Im,频率为ω的正弦交变电流信号。螺线管线的线电流可等效为柱面电流,基于调制器结构的对称性,螺线管内部的电磁场表达式为
Figure BDA0002108557420000025
B、E分别代表调制器的磁场和电场,Br为磁场的径向分量,z代表中轴线坐标,t为时间变量,r、
Figure BDA0002108557420000026
z分别为柱坐标系中沿径向、周向、轴向方向的单位矢量;
步骤1.2:将电磁场表达式代入麦克斯韦方程组,得到真空中磁场、电场各分量的关系式为
Figure BDA0002108557420000027
Figure BDA0002108557420000028
Figure BDA0002108557420000029
式中,Br为磁场的径向分量;
步骤1.3:忽略Br分量,即令Br=0,结合式(5)和式(6)得到交变电流磁光调制下调制器内部电磁场的基本模型为
Figure BDA0002108557420000031
Figure BDA0002108557420000032
步骤2:由交变电流磁光调制下调制器内部电磁场的基本模型求解空心螺线管的磁场为
Bz(r,t)=B0J0(k0r)·e-iωt (7)
式中,J0(k0r)是第一类零阶贝塞尔函数,B0为直流电信号激励时螺线管内部磁场幅值大小,即:
Figure BDA0002108557420000033
式中,
Figure BDA0002108557420000034
Figure BDA0002108557420000035
Figure BDA0002108557420000036
Figure BDA0002108557420000037
分别为第一类和第三类全椭圆积分;
步骤2.1:由式(1)得到螺线管内部轴向磁场Bz的通解为
Bz(r)=AJ0(k0r)+CN0(k0r) (9)
式中,A、C为待求积分常数,N0(k0r)是第二类零阶贝塞尔函数。
步骤2.2:确定常数A、C的值
当ω→0时因为Bz(ω→0)=AJ0(k0r→0)+CN0(k0r→0)=B0,而J0(k0r→0)=1,N0(k0r→0)→-∞,所以A=B0,C=0;于是得到式7中空心螺线管的磁场模型
Bz(r,t)=B0J0(k0r)·e-iωt (7)
步骤3:构建无松弛极化磁光材料内部磁场模型
对于无松弛极化的介质,在外电场作用下,其内部极化强度P和电位移矢量D以及电场强度E的关系为
D=ε0E+P=(1+αe0E=εrε0E=εE (10)
式中,ε为介质的介电常数,εr为介质的相对介电常数,αe为介质的极化率。将式(10)代入到麦克斯韦方程组中,重复步骤1、2得到介质中的磁场表达式为
Bz(r,t)=B0J0(k0r)·e-iωt (7)
式中,
Figure BDA0002108557420000041
为光在介质中的波矢。
步骤4:构建交变电流磁光调制下调制器内部电磁场的最终模型
Figure BDA0002108557420000042
式中,R为介质的半径。
步骤4.1:根据步骤2、3汇总得到交变电流磁光调制下调制器内部磁场可分为两部分
Bzr(r,t)=B0J0(kr)·e-iωt,r<R (12)
Bzo(r,t)=AJ0(k0r)·e-iωt,R<r<a (13)
式中,Bzr表示介质内部的磁场,Bzo表示介质外部的磁场,A为待求系数。
步骤4.2:根据步骤4.1的结果,结合安培环路定律得到磁光材料边界处微小闭合回路内的介质内磁场和介质外的磁场满足
Bzr(R)-Bzo(R)=0 (14)
由式(12)、(13)、(14)解得A=B0·J0(kR)/J0(k0R),于是得到交变电流磁光调制下调制器内部电磁场的最终模型为
Figure BDA0002108557420000043
同现有技术比较,本发明的最大优越性在于:通过对磁光调制过程中含无松弛介质的调制器内部磁场分布计算,为详细研究磁光调制过程提供了磁场分析方法。
附图说明
图1:磁光调制器示意图
图2磁场幅值Bz随频率变化关系
具体实施方式
现结合附图对本发明具体实施方式阐述如下:
设定基础仿真参数:电流信号幅值Im=1A,线圈密度n=1000匝/m,圆螺线管半径a=3.5cm,半长b=5.5cm,介质半径R=0.5cm,采用磁光介质TGG,其相对介电常数εr=4.1,相对磁导率μr=1,空气的介电常数和磁导率均设为1。
步骤1:根据圆调制器轴对称的结构特点,结合麦克斯韦方程组构建交变电流磁光调制下调制器内部电磁场的基本模型
Figure BDA0002108557420000051
Figure BDA0002108557420000052
步骤2:由交变电流磁光调制下调制器内部电磁场的基本模型求解空心螺线管的磁场为
Bz(r,t)=B0J0(k0r)·e-iωt (7)
其中
Figure BDA0002108557420000053
式中
Figure BDA0002108557420000054
Figure BDA0002108557420000055
Figure BDA0002108557420000056
分别为第一类和第三类全椭圆积分。
步骤3:构建无松弛极化磁光材料内部磁场模型
Figure BDA0002108557420000057
步骤4:构建交变电流磁光调制下调制器内部电磁场的最终模型
Figure BDA0002108557420000058
分别选取r=0.4R、4.0R两个点位,研究不同位置处调制信号平率变化对磁场幅值Bz的影响,如图2所示。图2反映了介质内外磁场幅值的频率特性,其中虚线表示介质内部r=0.4R处磁场,实线表示无介质时r=0.4R处磁场,o线表示介质外部r=4R处磁场,+线表示无介质时r=4R处磁场。

Claims (1)

1.一种交变电流磁光调制下的调制器磁场分析方法,其特征在于:通过分析调制器结构特点,结合麦克斯韦方程组分别建立交变电流磁光调制下调制器内部电磁场的基本模型、空心调制器内部磁场模型、无松弛极化磁光材料内部磁场模型,最终构建交变电流磁光调制下含无松弛极化磁光材料的调制器内部磁场分布模型,具体包括以下步骤:
步骤1:根据调制器结构特点,结合麦克斯韦方程组构建交变电流磁光调制下调制器内部电磁场的基本模型
步骤1.1:给调制器的螺线管线圈通以正弦交变电流Is=Ime-iωt,即一个幅值为Im,频率为ω的正弦交变电流信号,螺线管线的线电流可等效为柱面电流,基于调制器结构的对称性,螺线管内部的电磁场表达式为
Figure FDA0002662659440000011
B、E分别代表调制器的磁场和电场,Br为磁场的径向分量,Bz为磁场的轴向分量,
Figure FDA0002662659440000019
为电场的圆周方向分量,z代表中轴线坐标,t为时间变量,r、
Figure FDA0002662659440000012
z分别为柱坐标系中沿径向、周向、轴向方向的单位矢量;
步骤1.2:将电磁场表达式(1)代入麦克斯韦方程组,得到真空中磁场、电场各分量的关系式为
Figure FDA0002662659440000013
Figure FDA0002662659440000014
Figure FDA0002662659440000015
步骤1.3:忽略Br分量,即令Br=0,结合式(3)和式(4)得到交变电流磁光调制下调制器内部电磁场的基本模型为
Figure FDA0002662659440000016
Figure FDA0002662659440000017
式中,r代表距中轴线的距离,
Figure FDA0002662659440000018
为真空中的波矢,ε0为真空介电常数,μ0为真空磁导率,ω为加载的正弦交变电流信号的频率;
步骤2:由交变电流磁光调制下调制器内部电磁场的基本模型求解空心螺线管的磁场
步骤2.1:由式(5)得到螺线管内部轴向磁场Bz的通解为
Bz(r)=AJ0(k0r)+CN0(k0r) (7)
式中,A、C为待求积分常数,N0(k0r)是第二类零阶贝塞尔函数;
步骤2.2:确定常数A、C的值
当ω→0时因为Bz(ω→0)=AJ0(k0r→0)+CN0(k0r→0)=B0,而J0(k0r→0)=1,N0(k0r→0)→-∞,所以A=B0,C=0;于是得到空心螺线管的磁场
Bz(r,t)=B0J0(k0r)·e-iωt (8)
式中,J0(k0r)是第一类零阶贝塞尔函数,B0为直流电信号激励时螺线管内部磁场幅值大小,即:
Figure FDA0002662659440000021
式中,
Figure FDA0002662659440000022
Figure FDA0002662659440000023
Figure FDA0002662659440000024
Figure FDA0002662659440000025
分别为第一类和第三类全椭圆积分;
步骤3:构建无松弛极化磁光材料内部磁场模型
对于无松弛极化的介质,在外电场作用下,其内部极化强度P和电位移矢量D以及电场强度E的关系为
D=ε0E+P=(1+αe0E=εrε0E=εE (10)
式中,ε为介质的介电常数,εr为介质的相对介电常数,αe为介质的极化率,将式(10)代入到麦克斯韦方程组中,重复步骤1、2得到介质中的磁场表达式为
Bz(r,t)=B0J0(kr)·e-iωt (11)
式中,
Figure FDA0002662659440000026
为光在介质中的波矢;
步骤4:最终构建交变电流磁光调制下含无松弛极化磁光材料的调制器内部磁场分布模型
步骤4.1:根据步骤2、3汇总得到交变电流磁光调制下调制器内部磁场可分为两部分
Bzr(r,t)=B0J0(kr)·e-iωt,r<R (12)
Bzo(r,t)=AJ0(k0r)·e-iωt,R<r<a (13)
式中,Bzr表示介质内部的磁场,Bzo表示介质外部的磁场,A为待求系数;
步骤4.2:根据步骤4.1的结果,结合安培环路定律得到磁光材料边界处微小闭合回路内的介质内磁场和介质外的磁场满足
Bzr(R)-Bzo(R)=0 (14)
由式(12)、(13)、(14)解得A=B0·J0(kR)/J0(k0R),于是最终得到交变电流磁光调制下含无松弛极化磁光材料的调制器内部磁场分布模型为
Figure FDA0002662659440000031
式中,R为介质的半径。
CN201910562101.0A 2019-06-26 2019-06-26 一种交变电流磁光调制下调制器磁场分析方法 Active CN110221230B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910562101.0A CN110221230B (zh) 2019-06-26 2019-06-26 一种交变电流磁光调制下调制器磁场分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910562101.0A CN110221230B (zh) 2019-06-26 2019-06-26 一种交变电流磁光调制下调制器磁场分析方法

Publications (2)

Publication Number Publication Date
CN110221230A CN110221230A (zh) 2019-09-10
CN110221230B true CN110221230B (zh) 2020-12-11

Family

ID=67815062

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910562101.0A Active CN110221230B (zh) 2019-06-26 2019-06-26 一种交变电流磁光调制下调制器磁场分析方法

Country Status (1)

Country Link
CN (1) CN110221230B (zh)

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2514114A (en) * 2013-05-13 2014-11-19 Univ Bath Apparatus and method for measuring electromagnetic properties
CN106405457B (zh) * 2016-08-29 2018-11-13 中国科学院武汉物理与数学研究所 一种用于材料铁磁性及磁化性能检测的装置和方法
US10948556B2 (en) * 2017-11-30 2021-03-16 Electronics And Telecommunications Research Institute Method for modifying and controlling magnetic field and apparatus for the same
CN109490798B (zh) * 2018-09-26 2020-02-21 上海交通大学 基于光纤磁场传感技术的高温超导线圈磁场测量方法
CN109933911B (zh) * 2019-03-15 2023-04-07 中国人民解放军陆军装甲兵学院 密绕螺线管内金属圆筒电磁场有限元分析方法

Also Published As

Publication number Publication date
CN110221230A (zh) 2019-09-10

Similar Documents

Publication Publication Date Title
CN102565728B (zh) 一种圆柱形软磁材料磁导率测量方法
CN109933911A (zh) 密绕螺线管内金属圆筒电磁场有限元分析方法
Wolski Maxwell's equations for magnets
CN110221230B (zh) 一种交变电流磁光调制下调制器磁场分析方法
US6377148B1 (en) Asymmetric zonal shim coils for magnetic resonance
Rosenthal et al. Torque measurements in spin-up flow of ferrofluids
Skumiel et al. The use of Gramme coils in a 2-phase system for generation of a high frequency rotating magnetic field
Li et al. Rapid prototyping eddy current sensors using 3D printing
Ghadirian et al. Restricted diffusion in annular geometrical pores
Weir et al. Approximate shape factors for soft magnetic composites
Pankrac et al. Simple algorithms for the calculation of the intensity of the magnetic field of current loops and thin-wall air coils of a general shape using magnetic dipoles
JP2003512629A (ja) 永久磁石とシムの設計
RU2617723C2 (ru) Способ определения намагниченности вещества
Gerus et al. Effect of nonuniform magnetic field That magnetizes a ferrite film on the measurement accuracy for characteristics of spin waves
Mohdeb et al. A new approximation for calculating the attraction force in cylindrical permanent magnets arrays and cylindrical linear single-axis-actuator
Yating et al. Two approaches to coil impedance calculation of eddy current sensor
CN107748813B (zh) 非晶丝在非轴向磁场作用下的巨磁阻抗建模方法
Wei et al. Magnetic field analysis of solenoid driven by alternating current
Yang et al. Excitation structure design and magnetic field analysis of a new electromagnetic flowmeter based on magnetically permeable material
Wu et al. Analytical expression to impedance for solenoid coil with a coaxial cylindrical ferrite core of finite length
Emelyanov et al. Estimation of physical parameters of magnetic nanoparticles
Yao et al. Numerical study on magnetic field characteristics of electromagnetic flowmeter with small excitation module
Kosek et al. Skin effect in conductor of rectangular cross-section–approximate solution
Junxia et al. Design and Implementation of Portable NMR Probe Magnet
Ma et al. Analytical calculation method of free static magnetic field point in magnetic signal measurement

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