CN102945321A - 一种条形非线性永磁体等效磁路模型的建模方法 - Google Patents

一种条形非线性永磁体等效磁路模型的建模方法 Download PDF

Info

Publication number
CN102945321A
CN102945321A CN2012104214138A CN201210421413A CN102945321A CN 102945321 A CN102945321 A CN 102945321A CN 2012104214138 A CN2012104214138 A CN 2012104214138A CN 201210421413 A CN201210421413 A CN 201210421413A CN 102945321 A CN102945321 A CN 102945321A
Authority
CN
China
Prior art keywords
permanent magnet
magnetic
section
leakage
force
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
CN2012104214138A
Other languages
English (en)
Other versions
CN102945321B (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.)
Harbin Institute of Technology
Original Assignee
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201210421413.8A priority Critical patent/CN102945321B/zh
Publication of CN102945321A publication Critical patent/CN102945321A/zh
Application granted granted Critical
Publication of CN102945321B publication Critical patent/CN102945321B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)

Abstract

一种条形非线性永磁体等效磁路模型的建模方法,其包括以下步骤:步骤S1:对永磁体沿长度方向进行饱和充磁,然后放置于空气中自然去磁;步骤S2:对永磁体外中心截面进行磁力线绘制;步骤S3:对磁力线图景中射出或射入永磁体表面的位置进行标记;步骤S4:记录各个标记点对于条形永磁体中心位置的距离,计算每一段的等效长度;步骤S5:计算每一个分段漏磁的漏磁导;步骤S6:计算出总漏磁导并减去各个分段的漏磁导,得到端部漏磁导;步骤S7:结合所述各段永磁的漏磁导、等效长度以及端部漏磁导,建立基于磁力线的非线性条形永磁体的等效磁路模型;步骤S8:建模过程结束。本发明提高了磁路计算条形非线性永磁体截面磁通的计算精度。

Description

一种条形非线性永磁体等效磁路模型的建模方法
技术领域
本发明涉及一种建模方法,尤其涉及一种基于磁力线划分的条形非线性永磁体等效磁路模型的建模方法。
背景技术
目前,条形永磁体截面磁通的求解方法主要有基于磁场和基于磁路两种方法。
基于磁场的求解方法一般是利用商业电磁场仿真软件进行建模求解,该方法的优点是计算精确,但模型建立复杂繁琐,修改不便,计算时间较长,对于非线性永磁体建模较为复杂。
而基于磁路的求解方法模型简单、速度快,但计算精度低。在进行可靠性分析等需要大量计算的场合中,主要以磁路求解的方法为主。分析等效磁路法计算精度低的原因可以发现,以往都是将永磁看成一个整体,等效为一个磁势和一个磁阻,工作在一条回复线上。而实际情况是,永磁内部磁场分布不均匀,永磁各部分并非工作在一条回复线上,这一点在非线性永磁体上体现的尤为突出。
目前已有的分段建模过程中,通常将分段建立在条形永磁体的几何尺寸均分的基础上,即建立等效磁路模型时,模型每一段的长度尺寸一致。
因此,如何将上述问题加以解决,而提供一种基于磁力线划分的条形非线性永磁体等效磁路模型的建模方法,即为本领域技术人员的研究方向所在。
发明内容
本发明的主要目的在于提供一种条形非线性永磁体的等效磁路模型建模方法,以解决上述现有技术中所存在的问题。
为了达到上述目的,本发明提供一种条形非线性永磁体等效磁路模型的建模方法,其包括以下步骤:
步骤S 1:对永磁体沿长度方向进行饱和充磁,然后放置于空气中自然去磁;
步骤S2:对永磁体外中心截面进行磁力线绘制,得到永磁体外中心截面的磁场图景,分段数n由磁力线图景中的磁力线条数确定;
步骤S3:对磁力线图景中射出或射入永磁体表面的位置进行标记;
步骤S4:根据标记的磁力线射出或射入的位置,记录各个标记点对于条形永磁体中心位置的距离Li,计算每一段的等效长度li,li=Li+1-Li,i≥1;
步骤S5:计算每一个分段漏磁的漏磁导;
步骤S6:计算出总漏磁导并减去各个分段的漏磁导,得到端部漏磁导;
步骤S7:结合所述各段永磁的漏磁导、等效长度以及端部漏磁导,建立基于磁力线的非线性条形永磁体的等效磁路模型;
步骤S8:建模过程结束。
其中,在步骤S2中,是使用永磁体磁场图景绘制仪或采用播撒磁粉的方法对永磁体外中心截面进行磁力线绘制。
其中,所述采用播撒磁粉的方法是将永磁体水平放置于纯白色蜡制或纸质平面上,均匀播撒细密铁粉,即显现该平面磁场图景及磁力线。
其中,在步骤S2中,分段数n由磁力线图景中从永磁体上表面的磁力线条数确定,其中,上表面包括永磁端部上端点,不包括侧面及下表面。
其中,在步骤S5中,是根据漏磁导计算公式
Figure BDA00002325180200021
计算每一个分段漏磁的漏磁导,其中μ0为空气磁导率,bi=Li+Li-1,θi1和θi2分别为该段漏磁导磁力线的弧度。
其中,在步骤S6中,计算永磁体总漏磁导的计算公式为:
Figure BDA00002325180200031
其中,μ0为空气磁导率,L、D、H分别为永磁体的长、宽、高。
其中,在步骤S7中,设分段数n为10,十段非线性条形永磁体的等效磁路模型的数学表达为:
B ‾ 1 - B 1 = - 1 4 μ 1 H 1 ‾ B i - 1 - B ‾ i = - G i A ( 1 2 Σ δ = 1 i - 1 H δ ‾ μ δ + 1 4 H i ‾ μ i ) B i = B ‾ i - ( B i - 1 - B ‾ i ) B ‾ i = f ( H ‾ i )
其中,模型的端部校验条件为:
Figure BDA00002325180200033
其中,G为端部漏磁导,Bi,Hi为各段截面的磁场强度和磁感应强度,
Figure BDA00002325180200034
为各段截面的平均磁场强度和平均磁感应强度,A为永磁体截面积,f为非线性永磁体去磁曲线函数,μi为各段相对磁导率,μi=Gili/A。
本发明大大提高了磁路计算条形非线性永磁体截面磁通的计算精度,为非线性条形永磁体的截面磁通计算提供有效的磁路计算方法。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明所述的非线性条形永磁体示意图;
图2为本发明所述的条形永磁体中心面的磁场图景;
图3为本发明所述的方法应用于的条形永磁体的等效磁路模型;
图4为本发明步骤流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明所述的非线性条形永磁体示意图;中心截面为图中虚线所示平面,中心点为O,N、S极在图中标示,长宽高分别为L、D、H,对于条形永磁体有L>D>H。
如图2所示,本发明所述的条形永磁体中心面的磁场图景,分段根据图中标记的磁力线射出或射入的位置(图中永磁体上表面空心圆点位置),各段编号为1,2,3…n,各个标记点相对于条形永磁体中心位置的距离为l1,l2,l3…ln,每一段的等效长度即为l1-l2,l2-l3,…ln-1-ln
如图3所示,本发明所述的方法应用于的条形永磁体的等效磁路模型;其中,L为条形永磁体长度;A为条形永磁体截面积;Φi~Φi+1为通过每分段两个端面的截面磁通;Φsi为每段的漏磁通,Gi~Gn为各段的漏磁导;Gd为永磁体的端部漏磁导;ΔUi为各段的等效磁势差,Ui为永磁各分段的磁势。
如图4所示,为本发明步骤流程图,本发明提供了一种非线性条形永磁体等效磁路模型的建模方法,其包括以下步骤:
步骤S1:对永磁体沿长度方向进行饱和充磁,如图1所示,然后放置于空气中自然去磁,此时条形永磁体处于磁性能相对稳定状态;
步骤S2:使用永磁体磁场图景绘制仪或采用播撒磁粉的方法,对永磁体外中心截面进行磁力线绘制,得到永磁体外中心截面的磁场图景,如图2所示,分段数n由磁力线图景中从永磁体上表面(包括永磁端部上端点,不包括侧面及下表面)的磁力线条数确定,如图2中,永磁体上表面磁力线条数n=10,即分段数为10;
其中,上述采用播撒磁粉的方法是将永磁体水平放置于纯白色蜡制或纸质平面上,均匀播撒细密铁粉,即可显现该平面磁场图景及磁力线。
步骤S3:对磁力线图景中射出或射入永磁体表面的位置进行标记;
步骤S4:根据标记的磁力线射出或射入的位置,记录各个标记点对于条形永磁体中心位置的距离Li(i≥1),如图2中所示,计算每一段的等效长度li(li=Li+1-Li);
步骤S5:根据漏磁导计算公式
Figure BDA00002325180200051
计算每一个分段漏磁的漏磁导,其中μ0为空气磁导率,bi=Li+Li-1,θi1和θi2分别为该段漏磁导磁力线的弧度,如图2中所示;
步骤S6:根据永磁体总漏磁导计算公式计算出总漏磁导并减去各个分段的漏磁导,得到端部漏磁导;
其中,永磁体总漏磁导计算公式为:
Figure BDA00002325180200052
其中,μ0为空气磁导率,L、D、H分别为永磁体的长、宽、高,如图1所示;
步骤S7:结合所述各段永磁的漏磁导、等效长度以及端部漏磁导,建立基于磁力线的非线性条形永磁体的等效磁路模型,如图3所示;其中各参量均由前述步骤计算而得,依据磁路基本原理,图3中的十段非线性条形永磁体的等效磁路模型的数学表达为:
B ‾ 1 - B 1 = - 1 4 μ 1 H 1 ‾ B i - 1 - B ‾ i = - G i A ( 1 2 Σ δ = 1 i - 1 H δ ‾ μ δ + 1 4 H i ‾ μ i ) B i = B ‾ i - ( B i - 1 - B ‾ i ) B ‾ i = f ( H ‾ i )
其中,模型的端部校验条件为:
Figure BDA00002325180200054
其中G为端部漏磁导,Bi,Hi为各段截面的磁场强度和磁感应强度,
Figure BDA00002325180200055
为各段截面的平均磁场强度和平均磁感应强度,A为永磁体截面积,f为非线性永磁体去磁曲线函数,μi为各段相对磁导率,μi=Gili/A;
步骤S8:建模过程结束。
综上所述,本发明所述的基于磁力线的非线性条形永磁体等效磁路模型的建模方法,与实际情况更接近,大大提高了磁路计算条形永磁体磁路计算方法的计算精度,为含有条形永磁体磁系统的计算提供有效的磁路计算基础。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围。

Claims (7)

1.一种条形非线性永磁体等效磁路模型的建模方法,其特征在于,其包括以下步骤:
步骤S1:对永磁体沿长度方向进行饱和充磁,然后放置于空气中自然去磁;
步骤S2:对永磁体外中心截面进行磁力线绘制,得到永磁体外中心截面的磁场图景,分段数n由磁力线图景中的磁力线条数确定;
步骤S3:对磁力线图景中射出或射入永磁体表面的位置进行标记;
步骤S4:根据标记的磁力线射出或射入的位置,记录各个标记点对于条形永磁体中心位置的距离Li,计算每一段的等效长度li,li=Li+1-Li,i≥1;
步骤S5:计算每一个分段漏磁的漏磁导;
步骤S6:计算出总漏磁导并减去各个分段的漏磁导,得到端部漏磁导;
步骤S7:结合所述各段永磁的漏磁导、等效长度以及端部漏磁导,建立基于磁力线的非线性条形永磁体的等效磁路模型;
步骤S8:建模过程结束。
2.根据权利要求1所述的一种条形非线性永磁体等效磁路模型的建模方法,其特征在于,在步骤S2中,是使用永磁体磁场图景绘制仪或采用播撒磁粉的方法对永磁体外中心截面进行磁力线绘制。
3.根据权利要求2所述的一种条形非线性永磁体等效磁路模型的建模方法,其特征在于,所述采用播撒磁粉的方法是将永磁体水平放置于纯白色蜡制或纸质平面上,均匀播撒细密铁粉,即显现该平面磁场图景及磁力线。
4.根据权利要求1所述的一种条形非线性永磁体等效磁路模型的建模方法,其特征在于,在步骤S2中,分段数n由磁力线图景中从永磁体上表面的磁力线条数确定,其中,上表面包括永磁端部上端点,不包括侧面及下表面。
5.根据权利要求1所述的一种条形非线性永磁体等效磁路模型的建模方法,其特征在于,在步骤S5中,是根据漏磁导计算公式
Figure FDA00002325180100021
计算每一个分段漏磁的漏磁导,其中μ0为空气磁导率,bi=Li+Li-1,θi1和θi2分别为该段漏磁导磁力线的弧度。
6.根据权利要求5所述的一种条形非线性永磁体等效磁路模型的建模方法,其特征在于,在步骤S6中,计算永磁体总漏磁导的计算公式为:其中,μ0为空气磁导率,L、D、H分别为永磁体的长、宽、高。
7.根据权利要求6所述的一种条形非线性永磁体等效磁路模型的建模方法,其特征在于,在步骤S7中,设分段数n为10,十段非线性条形永磁体的等效磁路模型的数学表达为:
B ‾ 1 - B 1 = - 1 4 μ 1 H 1 ‾ B i - 1 - B ‾ i = - G i A ( 1 2 Σ δ = 1 i - 1 H δ ‾ μ δ + 1 4 H i ‾ μ i ) B i = B ‾ i - ( B i - 1 - B ‾ i ) B ‾ i = f ( H ‾ i )
其中,模型的端部校验条件为:
Figure FDA00002325180100024
其中,G为端部漏磁导,Bi,Hi为各段截面的磁场强度和磁感应强度,
Figure FDA00002325180100025
为各段截面的平均磁场强度和平均磁感应强度,A为永磁体截面积,f为非线性永磁体去磁曲线函数,μi为各段相对磁导率,μi=Gili/A。
CN201210421413.8A 2012-10-29 2012-10-29 一种条形非线性永磁体等效磁路模型的建模方法 Active CN102945321B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210421413.8A CN102945321B (zh) 2012-10-29 2012-10-29 一种条形非线性永磁体等效磁路模型的建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210421413.8A CN102945321B (zh) 2012-10-29 2012-10-29 一种条形非线性永磁体等效磁路模型的建模方法

Publications (2)

Publication Number Publication Date
CN102945321A true CN102945321A (zh) 2013-02-27
CN102945321B CN102945321B (zh) 2015-12-09

Family

ID=47728264

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210421413.8A Active CN102945321B (zh) 2012-10-29 2012-10-29 一种条形非线性永磁体等效磁路模型的建模方法

Country Status (1)

Country Link
CN (1) CN102945321B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106897487A (zh) * 2017-01-13 2017-06-27 华中科技大学 一种高温超导磁体非线性分析的建模方法
CN107153746A (zh) * 2017-06-02 2017-09-12 山东大学 一种内置式永磁同步电机漏磁系数解析计算方法
CN108052764A (zh) * 2017-12-26 2018-05-18 西安交通大学 一种环形Halbach永磁阵列的磁场建模方法
CN109633281A (zh) * 2018-11-29 2019-04-16 哈尔滨理工大学 一种线性/非线性双层绝缘电介质界面极化特性分析方法
CN111625983A (zh) * 2020-05-30 2020-09-04 哈尔滨工业大学 基于磁滞模型的含非线性永磁电磁机构有限元仿真方法
CN113640717A (zh) * 2021-07-28 2021-11-12 哈尔滨工业大学 一种电器元件内部条形永磁体截面磁通间接测量方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
NIKOLAY KARDJILOV ET.AL: "Three-dimensional imaging of magnetic fields with polarised neutrons", 《NATURE PHYSICS》 *
徐俊: "条形永磁分段回复线模型的研究", 《中国优秀硕士学位论文全文数据库》 *
梁慧敏 等: "磁钢内部磁通测量装置的研制", 《机电元件》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106897487A (zh) * 2017-01-13 2017-06-27 华中科技大学 一种高温超导磁体非线性分析的建模方法
CN106897487B (zh) * 2017-01-13 2018-12-14 华中科技大学 一种高温超导磁体非线性分析的建模方法
CN107153746A (zh) * 2017-06-02 2017-09-12 山东大学 一种内置式永磁同步电机漏磁系数解析计算方法
CN107153746B (zh) * 2017-06-02 2019-12-06 山东大学 一种内置式永磁同步电机漏磁系数解析计算方法
CN108052764A (zh) * 2017-12-26 2018-05-18 西安交通大学 一种环形Halbach永磁阵列的磁场建模方法
CN109633281A (zh) * 2018-11-29 2019-04-16 哈尔滨理工大学 一种线性/非线性双层绝缘电介质界面极化特性分析方法
CN111625983A (zh) * 2020-05-30 2020-09-04 哈尔滨工业大学 基于磁滞模型的含非线性永磁电磁机构有限元仿真方法
CN111625983B (zh) * 2020-05-30 2022-07-15 哈尔滨工业大学 基于磁滞模型的含非线性永磁电磁机构有限元仿真方法
CN113640717A (zh) * 2021-07-28 2021-11-12 哈尔滨工业大学 一种电器元件内部条形永磁体截面磁通间接测量方法
CN113640717B (zh) * 2021-07-28 2023-10-03 哈尔滨工业大学 一种电器元件内部条形永磁体截面磁通间接测量方法

Also Published As

Publication number Publication date
CN102945321B (zh) 2015-12-09

Similar Documents

Publication Publication Date Title
CN102945321A (zh) 一种条形非线性永磁体等效磁路模型的建模方法
US10551350B2 (en) Method for simulating magnetic flux leakage based on loop current
CN104991994A (zh) 一种基于j-a磁滞模型的变压器励磁电流仿真方法
CN111625983B (zh) 基于磁滞模型的含非线性永磁电磁机构有限元仿真方法
DK62481A (da) Fremgangsmaade til flerpolet magnetisering af strimmelmateriale samt apparat til udoevelse af fremgangsmaaden
CN112733429A (zh) 一种基于Jiles-Atherton模型的钕铁硼磁滞回线描述方法
CN103971443B (zh) 对防伪元件进行检测的方法和装置
CN103674349B (zh) 一种基于高斯定理的开环磁通量传感器感应方法
Nielsen et al. Saturation Of Gauge Invariant Schwinger Model Correlation Functions By Instantons
CN204348795U (zh) 钢壳电芯盖帽及钢壳电芯
Samofalov et al. Generation of strong inhomogeneous stray fields by high-anisotropy permanent magnets
DE60114219D1 (de) Vorrichtung und verfahren zur stabilisierung einer bahn oder eines fadens aus ferromagnetischem material, die/der sich in eine richtung bewegt
CN206347947U (zh) 一种工程造价专用的图纸比例尺
CN206489527U (zh) 一种永磁体的磁场调节装置
CN113312746B (zh) 一种预测永磁体服役性能的方法
GB2531948A (en) Semi-analytic inversion method for nuclear magnetic resonance (NMR) signal processing
CN102542163B (zh) 一种晶体织构极图分析方法
CN208141943U (zh) 充磁机极头微调机构
MX2021011193A (es) Codificacion de region de funcion fisica no clonable.
CN102930156A (zh) 一种判断永磁体周围磁场作用区域的方法
CN206281485U (zh) 一种地质剖面测量装置
Kaczér The interaction energy of parallel Bloch walls
CN206003373U (zh) 一种自制教学时轴
RU2018147382A (ru) Магнитная головка
CN106295036B (zh) 一种钕铁硼磁体各向异性场的计算方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant