CN110221230B - 一种交变电流磁光调制下调制器磁场分析方法 - Google Patents
一种交变电流磁光调制下调制器磁场分析方法 Download PDFInfo
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring 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:根据调制器结构特点,结合麦克斯韦方程组构建交变电流磁光调制下调制器内部电磁场的基本模型
步骤1.1:给调制器的螺线管线圈通以正弦交变电流Is=Ime-iωt,即一个幅值为Im,频率为ω的正弦交变电流信号。螺线管线的线电流可等效为柱面电流,基于调制器结构的对称性,螺线管内部的电磁场表达式为
步骤1.2:将电磁场表达式代入麦克斯韦方程组,得到真空中磁场、电场各分量的关系式为
式中,Br为磁场的径向分量;
步骤1.3:忽略Br分量,即令Br=0,结合式(5)和式(6)得到交变电流磁光调制下调制器内部电磁场的基本模型为
步骤2:由交变电流磁光调制下调制器内部电磁场的基本模型求解空心螺线管的磁场为
Bz(r,t)=B0J0(k0r)·e-iωt (7)
式中,J0(k0r)是第一类零阶贝塞尔函数,B0为直流电信号激励时螺线管内部磁场幅值大小,即:
步骤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+αe)ε0E=εrε0E=εE (10)
式中,ε为介质的介电常数,εr为介质的相对介电常数,αe为介质的极化率。将式(10)代入到麦克斯韦方程组中,重复步骤1、2得到介质中的磁场表达式为
Bz(r,t)=B0J0(k0r)·e-iωt (7)
步骤4:构建交变电流磁光调制下调制器内部电磁场的最终模型
式中,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),于是得到交变电流磁光调制下调制器内部电磁场的最终模型为
同现有技术比较,本发明的最大优越性在于:通过对磁光调制过程中含无松弛介质的调制器内部磁场分布计算,为详细研究磁光调制过程提供了磁场分析方法。
附图说明
图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:根据圆调制器轴对称的结构特点,结合麦克斯韦方程组构建交变电流磁光调制下调制器内部电磁场的基本模型
步骤2:由交变电流磁光调制下调制器内部电磁场的基本模型求解空心螺线管的磁场为
Bz(r,t)=B0J0(k0r)·e-iωt (7)
其中
步骤3:构建无松弛极化磁光材料内部磁场模型
步骤4:构建交变电流磁光调制下调制器内部电磁场的最终模型
分别选取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,频率为ω的正弦交变电流信号,螺线管线的线电流可等效为柱面电流,基于调制器结构的对称性,螺线管内部的电磁场表达式为
步骤1.2:将电磁场表达式(1)代入麦克斯韦方程组,得到真空中磁场、电场各分量的关系式为
步骤1.3:忽略Br分量,即令Br=0,结合式(3)和式(4)得到交变电流磁光调制下调制器内部电磁场的基本模型为
步骤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为直流电信号激励时螺线管内部磁场幅值大小,即:
步骤3:构建无松弛极化磁光材料内部磁场模型
对于无松弛极化的介质,在外电场作用下,其内部极化强度P和电位移矢量D以及电场强度E的关系为
D=ε0E+P=(1+αe)ε0E=εrε0E=εE (10)
式中,ε为介质的介电常数,εr为介质的相对介电常数,αe为介质的极化率,将式(10)代入到麦克斯韦方程组中,重复步骤1、2得到介质中的磁场表达式为
Bz(r,t)=B0J0(kr)·e-iωt (11)
步骤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),于是最终得到交变电流磁光调制下含无松弛极化磁光材料的调制器内部磁场分布模型为
式中,R为介质的半径。
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)
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 | 中国人民解放军陆军装甲兵学院 | 密绕螺线管内金属圆筒电磁场有限元分析方法 |
-
2019
- 2019-06-26 CN CN201910562101.0A patent/CN110221230B/zh active Active
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 |