CN110188480A - 一种直流偏磁条件下铁磁材料的磁滞特性模拟分析***和方法 - Google Patents

一种直流偏磁条件下铁磁材料的磁滞特性模拟分析***和方法 Download PDF

Info

Publication number
CN110188480A
CN110188480A CN201910471820.1A CN201910471820A CN110188480A CN 110188480 A CN110188480 A CN 110188480A CN 201910471820 A CN201910471820 A CN 201910471820A CN 110188480 A CN110188480 A CN 110188480A
Authority
CN
China
Prior art keywords
curve
point
single order
magnetic field
field strength
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
CN201910471820.1A
Other languages
English (en)
Other versions
CN110188480B (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.)
North China Electric Power University
Original Assignee
North China Electric Power University
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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN201910471820.1A priority Critical patent/CN110188480B/zh
Publication of CN110188480A publication Critical patent/CN110188480A/zh
Application granted granted Critical
Publication of CN110188480B publication Critical patent/CN110188480B/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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Magnetic Variables (AREA)
  • Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)

Abstract

本发明涉及一种直流偏磁条件下铁磁材料的磁滞特性模拟分析***和方法,所述方法包括:生成一阶回转曲线,并显示所述生成的一阶回转曲线;利用生成的一阶回转曲线进行模型参数辨识和磁滞特性模拟分析。本发明仅需极限磁滞回线的实验数据,所需实验数据少,避免了复杂的实验设计和测量工作以及由此造成的低效和误差,有效提高生成的直流偏磁条件下铁磁材料一阶回转曲线及Preisach磁滞模型模拟结果的准确性。同时还可以预测高阶回转曲线,使得磁滞特性的模拟分析更加准确和多样化。

Description

一种直流偏磁条件下铁磁材料的磁滞特性模拟分析***和 方法
技术领域
本发明属于铁磁材料磁滞特性研究领域,尤其涉及一种直流偏磁条件下铁磁材料的磁滞特性模拟分析***和方法。
背景技术
磁滞现象是磁性材料固有的重要特征之一,不同磁性材料的磁滞特性主要通过该材料的磁滞回线形状和影响参数来表现。而磁性材料几乎存在于所有的电气设备中,如电力***中变压器的铁心,发电机、电子线路中的电感线圈等,由于在铁磁材料内部存在磁滞现象和涡流现象,当内部磁场随时间变化时,便会产生铁耗,从而影响电气设备的运行性能。
实际运行中的电力变压器,由于磁暴以及直流输电的影响,会发生直流偏磁现象,从而导致变压器运行性能降低、运行噪声加大、结构件和箱体的局部过热,对变压器的稳定安全运行更为不利。因此,深入分析铁磁材料在直流偏磁条件下的磁滞特性及损耗特性,有助于直流偏磁问题的研究和解决。
Preisach模型是研究磁性材料磁滞特性、损耗特性的重要数学模型,一阶回转曲线是进行其参数辨识的重要数据,一般通过测量得到。相关的实验研究表明,与无偏磁条件下铁磁材料的磁滞特性相比,直流偏磁磁滞回线不再对称,直流偏磁磁化曲线在一、三象限不对称且不再通过原点,因此直流偏磁一阶回转曲线也必然与无偏磁一阶回转曲线不同。而在实际应用中,一阶回转曲线的测量工作十分复杂繁琐,并且测量实验本身也会对结果造成误差,因此找到一个利用尽可能少量的实验数据来获得准确有效的直流偏磁一阶回转曲线,建立Preisach模型,将使直流偏磁条件下磁滞模拟结果更加准确,对实现Preisach模型的参数辨识十分具有实际意义。
发明内容
本发明的目的在于克服上述不足之处,提供一种直流偏磁条件下铁磁材料的磁滞特性模拟分析***和方法,通过本发明提供的模拟分析***和方法,能够利用少量的实验数据来获得准确有效的直流偏磁一阶回转曲线,磁滞模拟分析结果更加准确,同时还可以预测高阶回转曲线,使得磁滞特性的模拟分析更加准确和多样化。
为实现上述目的本发明提供一种直流偏磁条件下铁磁材料的磁滞特性模拟分析***,所述***包括:
生成显示模块,用于生成一阶回转曲线,并显示所述生成的一阶回转曲线;
模拟分析模块,利用生成的一阶回转曲线进行模型参数辨识和磁滞特性模拟分析。
根据本发明的另一个方面,所述生成显示模块包括回转点确定子模块,用于确定一阶回转曲线的回转点位置,得到曲线上任一点处磁场强度,所述一阶回转曲线回转点位置如下:
所述一阶回转曲线回转点R对应的极限磁滞回线水平宽度,即磁场强度差值ΔHrev如下:
ΔHrev=Ha(BR)-Hd(BR) (1)
其中,Ha(B)为极限磁滞回线上升支磁场强度;Hd(B)为极限磁滞回线下降支磁场强度;Ha(BR)、Hd(BR)分别为回转点R对应的极限磁滞回线上升支和下降支磁场强度;
所述一阶回转曲线回转点R与正深度饱和点T之间的垂直宽度,即磁感应强度差值ΔBrev如下:
ΔBrev=BT-BR (2)
其中,BT为正深度饱和点T处磁感应强度;BR为回转点R处磁感应强度;
所述一阶回转曲线上任一点P处磁场强度HP如下:
HP=Ha(BP)-ΔH (3)
其中,BP为P点磁感应强度;Ha(BP)为BP对应的极限磁滞回线上升支磁场强度;ΔH为BP对应的极限磁滞回线上升支磁场强度与P点磁场强度差值。
根据本发明的另一个方面,所述生成显示模块包括相对位置确定子模块,用于引入无量纲参数,用比值表征一阶回转曲线回转点及曲线上任一点在曲线中的相对位置,所述无量纲参数表示如下:
用比值β表征回转点R在一阶回转曲线中的相对位置:
其中,ΔBout为极限磁滞回线正、负深度饱和点之间磁感应强度差值;
用比值x表征点P在一阶回转曲线中的相对位置:
其中,ΔB为点P与点T之间磁感应强度差值。
根据本发明的另一个方面,所述生成显示模块包括变换求解子模块,用于将同一磁感应强度对应的极限磁滞回线上升支磁场强度与曲线上任一点磁场强度差值按变化速率用两项分量叠加表示,并利用均方根逼近得出上述表达式中系数最优值,其中其各项参数取值决定变化速率;所述同一磁感应强度BP对应的极限磁滞回线上升支磁场强度Ha(BP)与P点磁场强度HP差值ΔH(x)表示如下:
ΔH(x)=ΔHrev·(1-b)xe-a(1-x)+ΔHout(BP)·bxc (6)
ΔHout(BP)=Ha(BP)-Hd(BP) (7)
其中,a、b、c为参数,其取值决定ΔH(x)变化速率;ΔHout(BP)为BP对应的极限磁滞回线水平宽度,即磁场强度差值。
根据本发明的另一个方面,所述生成显示模块包括曲线生成子模块,用于选取所得系数最优值,生成直流偏磁条件下铁磁材料一阶回转曲线;
所模拟分析模块进行磁滞特性模拟分析包括:利用生成的一阶回转曲线模拟预测高阶回转曲线,进行磁滞特性分析。
本发明还提供了一种直流偏磁条件下铁磁材料的磁滞特性模拟分析方法,所述方法包括以下步骤:
A1:生成一阶回转曲线,并显示所述生成的一阶回转曲线;
A2:利用生成的一阶回转曲线进行模型参数辨识和磁滞特性模拟分析。
根据本发明的另一个方面,所述步骤A1还包括步骤:
S1:确定一阶回转曲线回转点位置,得到曲线上任一点处磁场强度;所述一阶回转曲线回转点位置如下:
所述一阶回转曲线回转点R对应的极限磁滞回线水平宽度,即磁场强度差值ΔHrev如下:
ΔHrev=Ha(BR)-Hd(BR) (1)
其中,Ha(B)为极限磁滞回线上升支磁场强度;Hd(B)为极限磁滞回线下降支磁场强度;Ha(BR)、Hd(BR)分别为回转点R对应的极限磁滞回线上升支和下降支磁场强度;
所述一阶回转曲线回转点R与正深度饱和点T之间的垂直宽度,即磁感应强度差值ΔBrev如下:
ΔBrev=BT-BR (2)
其中,BT为正深度饱和点T处磁感应强度;BR为回转点R处磁感应强度;
所述一阶回转曲线上任一点P处磁场强度HP如下:
HP=Ha(BP)-ΔH (3)
其中,BP为P点磁感应强度;Ha(BP)为BP对应的极限磁滞回线上升支磁场强度;ΔH为BP对应的极限磁滞回线上升支磁场强度与P点磁场强度差值。
根据本发明的另一个方面,所述步骤A1还包括步骤:
S2:引入无量纲参数,用比值表征一阶回转曲线回转点及曲线上任一点在曲线中的相对位置,所述无量纲参数表示如下:
用比值β表征回转点R在一阶回转曲线中的相对位置:
其中,ΔBout为极限磁滞回线正、负深度饱和点之间磁感应强度差值;
用比值x表征点P在一阶回转曲线中的相对位置:
其中,ΔB为点P与点T之间磁感应强度差值。
根据本发明的另一个方面,所述步骤A1还包括步骤:
S3:将同一磁感应强度对应的极限磁滞回线上升支磁场强度与曲线上任一点磁场强度差值按变化速率用两项分量叠加表示,并利用均方根逼近得出上述表达式中系数最优值,其中其各项参数取值决定变化速率;所述同一磁感应强度BP对应的极限磁滞回线上升支磁场强度Ha(BP)与P点磁场强度HP差值ΔH(x)表示如下:
ΔH(x)=ΔHrev·(1-b)xe-a(1-x)+ΔHout(BP)·bxc (6)
ΔHout(BP)=Ha(BP)-Hd(BP) (7)
其中,a、b、c为参数,其取值决定ΔH(x)变化速率;ΔHout(BP)为BP对应的极限磁滞回线水平宽度,即磁场强度差值。
根据本发明的另一个方面,所述方法还包括以下步骤:
选取所得系数最优值,生成直流偏磁条件下铁磁材料一阶回转曲线;
利用生成的一阶回转曲线模拟预测高阶回转曲线,进行磁滞特性分析。
本发明的有益效果是:
本发明将直流偏磁条件下铁磁材料一阶回转曲线中同一磁感应强度对应的极限磁滞回线上升支磁场强度与曲线上任一点磁场强度差值采用两项分量叠加表示,以反映所述曲线在不同分段趋近极限磁滞回线的变化速率,利用均方根拟合逼近的方法选取系数最优解,算法简单;其次,本发明仅需极限磁滞回线的实验数据,所需实验数据少,避免了复杂的实验设计和测量工作以及由此造成的低效和误差,有效提高生成的直流偏磁条件下铁磁材料一阶回转曲线及Preisach磁滞模型模拟结果的准确性。同时还可以预测高阶回转曲线,使得磁滞特性的模拟分析更加准确和多样化。
附图说明
图1为本发明的磁滞特性模拟分析***结构示意图;
图2为本发明的原理示意图;
图3为本发明中的磁滞特性模拟分析方法流程示意图;
图4为本发明的一个优选实施例中磁滞特性模拟分析方法流程示意图;
图5为本发明的由计算机执行的生成一阶回转曲线的方法流程示意图;
图6A和图6B为本发明实例中的一阶回转曲线生成图。
具体实施方式
下面结合附图并通过具体实施例对本发明作进一步详述,以下实施例只是描述性的,不是限定性的,不能以此限定本发明的保护范围。
如图1和图2所示,本发明提供了一种直流偏磁条件下铁磁材料的磁滞特性模拟分析***,所述***包括:
生成显示模块,用于生成一阶回转曲线,并显示所述生成的一阶回转曲线;
模拟分析模块,利用生成的一阶回转曲线进行模型参数辨识和磁滞特性模拟分析。
优选地,所述生成显示模块包括回转点确定子模块,用于确定一阶回转曲线的回转点位置,得到曲线上任一点处磁场强度,所述一阶回转曲线回转点位置如下:
所述一阶回转曲线回转点R对应的极限磁滞回线水平宽度,即磁场强度差值ΔHrev如下:
ΔHrev=Ha(BR)-Hd(BR) (1)
其中,Ha(B)为极限磁滞回线上升支磁场强度;Hd(B)为极限磁滞回线下降支磁场强度;Ha(BR)、Hd(BR)分别为回转点R对应的极限磁滞回线上升支和下降支磁场强度;
所述一阶回转曲线回转点R与正深度饱和点T之间的垂直宽度,即磁感应强度差值ΔBrev如下:
ΔBrev=BT-BR (2)
其中,BT为正深度饱和点T处磁感应强度;BR为回转点R处磁感应强度;
所述一阶回转曲线上任一点P处磁场强度HP如下:
HP=Ha(BP)-ΔH (3)
其中,BP为P点磁感应强度;Ha(BP)为BP对应的极限磁滞回线上升支磁场强度;ΔH为BP对应的极限磁滞回线上升支磁场强度与P点磁场强度差值。
优选地,所述生成显示模块包括相对位置确定子模块,用于引入无量纲参数,用比值表征一阶回转曲线回转点及曲线上任一点在曲线中的相对位置,所述无量纲参数表示如下:
用比值β表征回转点R在一阶回转曲线中的相对位置:
其中,ΔBout为极限磁滞回线正、负深度饱和点之间磁感应强度差值;
用比值x表征点P在一阶回转曲线中的相对位置:
其中,ΔB为点P与点T之间磁感应强度差值。
优选地,所述生成显示模块包括变换求解子模块,用于将同一磁感应强度对应的极限磁滞回线上升支磁场强度与曲线上任一点磁场强度差值按变化速率用两项分量叠加表示,并利用均方根逼近得出上述表达式中系数最优值,其中其各项参数取值决定变化速率;所述同一磁感应强度BP对应的极限磁滞回线上升支磁场强度Ha(BP)与P点磁场强度HP差值ΔH(x)表示如下:
ΔH(x)=ΔHrev·(1-b)xe-a(1-x)+ΔHout(BP)·bxc (6)
ΔHout(BP)=Ha(BP)-Hd(BP) (7)
其中,a、b、c为参数,其取值决定ΔH(x)变化速率;ΔHout(BP)为BP对应的极限磁滞回线水平宽度,即磁场强度差值。
优选地,所述生成显示模块包括曲线生成子模块,用于选取所得系数最优值,生成直流偏磁条件下铁磁材料一阶回转曲线;
所模拟分析模块进行磁滞特性模拟分析包括:利用生成的一阶回转曲线模拟预测高阶回转曲线,进行磁滞特性分析。
如图3所示,本发明还提供了一种直流偏磁条件下铁磁材料的磁滞特性模拟分析方法,所述方法包括以下步骤:
A1:生成一阶回转曲线,并显示所述生成的一阶回转曲线;
A2:利用生成的一阶回转曲线进行模型参数辨识和磁滞特性模拟分析。
优选地,所述步骤A1还包括步骤:
S1:确定一阶回转曲线回转点位置,得到曲线上任一点处磁场强度;所述一阶回转曲线回转点位置如下:
所述一阶回转曲线回转点R对应的极限磁滞回线水平宽度,即磁场强度差值ΔHrev如下:
ΔHrev=Ha(BR)-Hd(BR) (1)
其中,Ha(B)为极限磁滞回线上升支磁场强度;Hd(B)为极限磁滞回线下降支磁场强度;Ha(BR)、Hd(BR)分别为回转点R对应的极限磁滞回线上升支和下降支磁场强度;
所述一阶回转曲线回转点R与正深度饱和点T之间的垂直宽度,即磁感应强度差值ΔBrev如下:
ΔBrev=BT-BR (2)
其中,BT为正深度饱和点T处磁感应强度;BR为回转点R处磁感应强度;
所述一阶回转曲线上任一点P处磁场强度HP如下:
HP=Ha(BP)-ΔH (3)
其中,BP为P点磁感应强度;Ha(BP)为BP对应的极限磁滞回线上升支磁场强度;ΔH为BP对应的极限磁滞回线上升支磁场强度与P点磁场强度差值。
优选地,所述步骤A1还包括步骤:
S2:引入无量纲参数,用比值表征一阶回转曲线回转点及曲线上任一点在曲线中的相对位置,所述无量纲参数表示如下:
用比值β表征回转点R在一阶回转曲线中的相对位置:
其中,ΔBout为极限磁滞回线正、负深度饱和点之间磁感应强度差值;
用比值x表征点P在一阶回转曲线中的相对位置:
其中,ΔB为点P与点T之间磁感应强度差值。
优选地,所述步骤A1还包括步骤:
S3:将同一磁感应强度对应的极限磁滞回线上升支磁场强度与曲线上任一点磁场强度差值按变化速率用两项分量叠加表示,并利用均方根逼近得出上述表达式中系数最优值,其中其各项参数取值决定变化速率;所述同一磁感应强度BP对应的极限磁滞回线上升支磁场强度Ha(BP)与P点磁场强度HP差值ΔH(x)表示如下:
ΔH(x)=ΔHrev·(1-b)xe-a(1-x)+ΔHout(BP)·bxc (6)
ΔHout(BP)=Ha(BP)-Hd(BP) (7)
其中,a、b、c为参数,其取值决定ΔH(x)变化速率;ΔHout(BP)为BP对应的极限磁滞回线水平宽度,即磁场强度差值。
优选地,所述方法还包括以下步骤:
选取所得系数最优值,生成直流偏磁条件下铁磁材料一阶回转曲线;
利用生成的一阶回转曲线模拟预测高阶回转曲线,进行磁滞特性分析。
如图4所示,本实施例给出了一个具体的磁滞特性模拟分析方法流程示意图。该方法包括以下步骤:
S1:确定一阶回转曲线回转点R的位置:
一阶回转曲线回转点R对应的极限磁滞回线水平宽度,即磁场强度差值ΔHrev
ΔHrev=Ha(BR)-Hd(BR) (1)
其中,Ha(B)为极限磁滞回线上升支磁场强度;Hd(B)为极限磁滞回线下降支磁场强度;Ha(BR)、Hd(BR)分别为回转点R对应的极限磁滞回线上升支和下降支磁场强度;
一阶回转曲线回转点R与正深度饱和点T之间的垂直宽度,即磁感应强度差值ΔBrev
ΔBrev=BT-BR (2)
其中,BT为深度饱和点T处磁感应强度;BR为回转点R处磁感应强度;
一阶回转曲线上任一点P处磁场强度HP
HP=Ha(BP)-ΔH (3)
其中,BP为P点磁感应强度;Ha(BP)为BP对应的极限磁滞回线上升支磁场强度;ΔH为BP对应的极限磁滞回线上升支磁场强度与P点磁场强度差值。
S2:引入无量纲参数β及x:
比值β表征回转点R在一阶回转曲线中的相对位置:
其中,ΔBout为极限磁滞回线正、负深度饱和点之间磁感应强度差值;
比值x表征比值x表征点P在一阶回转曲线中的相对位置:
其中,ΔB为点P与点T之间磁感应强度差值。
S3:将同一磁感应强度BP对应的极限磁滞回线上升支磁场强度Ha(BP)与P点磁场强度HP差值ΔH(x)按变化速率用两项分量叠加表示:
ΔH(x)=ΔHrev·(1-b)xe-a(1-x)+ΔHout(BP)·bxc (6)
ΔHout(BP)=Ha(BP)-Hd(BP) (7)
其中,a、b、c为参数,其取值决定ΔH(x)变化速率;ΔHout(BP)为BP对应的极限磁滞回线水平宽度,即磁场强度差值。
S4:利用均方根逼近得出上述表达式中系数最优值。
系数a、b均与β有关,因此用关于β的多项式表示系数a、b,拟合结果如下:
a=ΔBrev(7.73+2.76β-28.63β2+28.36β3) (8)
b=0.22(1-β) (9)
c=0.125 (10)
S5:选取所得系数最优值,生成直流偏磁条件下铁磁材料一阶回转曲线,实现Preisach模型的参数辨识及磁滞特性模拟分析。
优选地,本步骤中,可以根据生成的一阶回转曲线模拟预测高阶回转曲线,进行磁滞特性分析。
一阶回转曲线定义为回转点在主磁滞回线或初始磁化曲线上的各回转曲线,二阶回转曲线定义为回转点在一阶回转曲线上的各回转曲线,依此类推。将二阶或高于二阶的回转曲线统称为高阶回转曲线。如何得到一阶回转曲线是数值Preisach模型的首要工作,一阶回转曲线数据测量的准确与否直接关系到对高阶回转曲线预测结果的准确性。一阶回转曲线应具备以下两个特征:第一,一阶回转曲线应满足描述磁滞特性的Madelung定则,即每条一阶回转曲线由回转点坐标唯一确定;第二,一阶回转曲线应该具备铁磁材料动态磁滞特性的基本信息,能够反映高阶回转曲线磁滞特性的基本特征。
进行Preisach模型的参数辨识时,Preisach分布函数可由极限磁滞回线上的回转点数据和一阶回转曲线上的数据确定。由一阶回转曲线点的数据可以得到每个网格关于一阶回转曲线点数据的输出函数,由该函数可计算出P平面内任意点所对应的一阶回转曲线上的输出值,从而求得整个P平面各个点对应的一阶回转曲线对应的输出值。根据Preisach模型的同余性和输入值的磁化历史,可以得到任意输入值H对应的输出值B,从而可以得到B-H磁滞回线。Preisach模型的同余性可以表述:当输入变量在相同的两个回转点之间来回变化时,所得的局部回线对应的垂直弦长相等。
模拟预测高阶回转曲线时,如果一阶回转曲线的同心磁滞回线符合高阶回转曲线的特点,则对同心磁滞回线的部分数据进行提取整理作为测量值,得到B-H高阶回转曲线。利用一阶回转曲线Preisach模型以高阶回转曲线磁场强度H为输入值进行计算得到磁感应强度B,从而得到B-H高阶回转曲线。
本发明实施例只需一阶回转线和极限磁滞回线的实验测量数据就能得到满足准确性要求的高阶回转曲线,大大提高了工作效率,节约了实验测量时间。
图5为本发明的由计算机执行的生成一阶回转曲线的方法流程示意图。
首先,读取数据及初始化变量N,M;然后,确定回转点R位置,对BR由Bmin到Bmax均匀采样N个值,即BR(i);对i赋值,令i=1,确定第i个回转点R对应的
β(i),a(i),b(i),c及ΔHrev(i)。判断条件1≤j≤M-1是否成立,如果成立,则根据ΔH表达式,确定第i个回转点对应的一阶回转曲线上第j个采样点P的位置,即P点坐标(HP,BP);j的值加1后返回继续判断条件1≤j≤M-1是否成立。判断i≤N-1是否成立,如果成立,则i的值加1返回继续迭代,否则,绘制所述N-1条一阶回转曲线。本流程可以通过编写代码导入到计算机中,由计算机运行该代码来自动生成一阶回转曲线。
图6A和图6B为本发明实例中的一阶回转曲线生成图。其中,图6A为上升支一阶回转曲线,图6B为下降支一阶回转曲线。
综上所述,本发明将直流偏磁条件下铁磁材料一阶回转曲线中同一磁感应强度对应的极限磁滞回线上升支磁场强度与曲线上任一点磁场强度差值采用两项分量叠加表示,以反映所述曲线在不同分段趋近极限磁滞回线的变化速率,利用均方根拟合逼近的方法选取系数最优值,算法简单;其次,本发明仅需极限磁滞回线的实验数据,所需实验数据少,避免了复杂的实验设计和测量工作以及由此造成的低效和误差,有效提高生成的直流偏磁条件下铁磁材料一阶回转曲线及Preisach磁滞模型模拟结果的准确性。同时还可以预测高阶回转曲线,使得磁滞特性的模拟分析更加准确和多样化。
尽管为说明目的公开了本发明的实施例和附图,但是本领域的技术人员可以理解:在不脱离本发明及所附权利要求的精神和范围内,各种替换、变化和修改都是可能的,因此,本发明的范围不局限于实施例和附图所公开的内容。

Claims (10)

1.一种直流偏磁条件下铁磁材料的磁滞特性模拟分析***,其特征在于:所述***包括:
生成显示模块,用于生成一阶回转曲线,并显示所述生成的一阶回转曲线;
模拟分析模块,利用生成的一阶回转曲线进行模型参数辨识和磁滞特性模拟分析。
2.根据权利要求1所述的模拟分析***,其特征在于:所述生成显示模块包括回转点确定子模块,用于确定一阶回转曲线的回转点位置,得到曲线上任一点处磁场强度,所述一阶回转曲线回转点位置如下:
所述一阶回转曲线回转点R对应的极限磁滞回线水平宽度,即磁场强度差值ΔHrev如下:
ΔHrev=Ha(BR)-Hd(BR) (1)
其中,Ha(B)为极限磁滞回线上升支磁场强度;Hd(B)为极限磁滞回线下降支磁场强度;Ha(BR)、Hd(BR)分别为回转点R对应的极限磁滞回线上升支和下降支磁场强度;
所述一阶回转曲线回转点R与正深度饱和点T之间的垂直宽度,即磁感应强度差值ΔBrev如下:
ΔBrev=BT-BR (2)
其中,BT为正深度饱和点T处磁感应强度;BR为回转点R处磁感应强度;
所述一阶回转曲线上任一点P处磁场强度HP如下:
HP=Ha(BP)-ΔH (3)
其中,BP为P点磁感应强度;Ha(BP)为BP对应的极限磁滞回线上升支磁场强度;ΔH为BP对应的极限磁滞回线上升支磁场强度与P点磁场强度差值。
3.根据权利要求2所述的模拟分析***,其特征在于:所述生成显示模块包括相对位置确定子模块,用于引入无量纲参数,用比值表征一阶回转曲线回转点及曲线上任一点在曲线中的相对位置,所述无量纲参数表示如下:
用比值β表征回转点R在一阶回转曲线中的相对位置:
其中,ΔBout为极限磁滞回线正、负深度饱和点之间磁感应强度差值;
用比值x表征点P在一阶回转曲线中的相对位置:
其中,ΔB为点P与点T之间磁感应强度差值。
4.根据权利要求3所述的模拟分析***,其特征在于:所述生成显示模块包括变换求解子模块,用于将同一磁感应强度对应的极限磁滞回线上升支磁场强度与曲线上任一点磁场强度差值按变化速率用两项分量叠加表示,并利用均方根逼近得出上述表达式中系数最优值,其中其各项参数取值决定变化速率;所述同一磁感应强度BP对应的极限磁滞回线上升支磁场强度Ha(BP)与P点磁场强度HP差值ΔH(x)表示如下:
ΔH(x)=ΔHrev·(1-b)xe-a(1-x)+ΔHout(BP)·bxc (6)
ΔHout(BP)=Ha(BP)-Hd(BP) (7)
其中,a、b、c为参数,其取值决定ΔH(x)变化速率;ΔHout(BP)为BP对应的极限磁滞回线水平宽度,即磁场强度差值。
5.根据权利要求4所述的模拟分析***,其特征在于:所述生成显示模块包括曲线生成子模块,用于选取所得系数最优值,生成直流偏磁条件下铁磁材料一阶回转曲线;
所模拟分析模块进行磁滞特性模拟分析包括:利用生成的一阶回转曲线模拟预测高阶回转曲线,进行磁滞特性分析。
6.一种直流偏磁条件下铁磁材料的磁滞特性模拟分析方法,应用于如权利要求1-5任一项所述的直流偏磁条件下铁磁材料的磁滞特性模拟分析***,其特征在于:所述方法包括以下步骤:
A1:生成一阶回转曲线,并显示所述生成的一阶回转曲线;
A2:利用生成的一阶回转曲线进行模型参数辨识和磁滞特性模拟分析。
7.根据权利要求6所述的模拟分析方法,其特征在于:所述步骤A1还包括步骤:
S1:确定一阶回转曲线回转点位置,得到曲线上任一点处磁场强度;所述一阶回转曲线回转点位置如下:
所述一阶回转曲线回转点R对应的极限磁滞回线水平宽度,即磁场强度差值ΔHrev如下:
ΔHrev=Ha(BR)-Hd(BR) (1)
其中,Ha(B)为极限磁滞回线上升支磁场强度;Hd(B)为极限磁滞回线下降支磁场强度;Ha(BR)、Hd(BR)分别为回转点R对应的极限磁滞回线上升支和下降支磁场强度;
所述一阶回转曲线回转点R与正深度饱和点T之间的垂直宽度,即磁感应强度差值ΔBrev如下:
ΔBrev=BT-BR (2)
其中,BT为正深度饱和点T处磁感应强度;BR为回转点R处磁感应强度;
所述一阶回转曲线上任一点P处磁场强度HP如下:
HP=Ha(BP)-ΔH (3)
其中,BP为P点磁感应强度;Ha(BP)为BP对应的极限磁滞回线上升支磁场强度;ΔH为BP对应的极限磁滞回线上升支磁场强度与P点磁场强度差值。
8.根据权利要求7所述的模拟分析方法,其特征在于:所述步骤A1还包括步骤:
S2:引入无量纲参数,用比值表征一阶回转曲线回转点及曲线上任一点在曲线中的相对位置,所述无量纲参数表示如下:
用比值β表征回转点R在一阶回转曲线中的相对位置:
其中,ΔBout为极限磁滞回线正、负深度饱和点之间磁感应强度差值;
用比值x表征点P在一阶回转曲线中的相对位置:
其中,ΔB为点P与点T之间磁感应强度差值。
9.根据权利要求8所述的模拟分析方法,其特征在于:所述步骤A1还包括步骤:
S3:将同一磁感应强度对应的极限磁滞回线上升支磁场强度与曲线上任一点磁场强度差值按变化速率用两项分量叠加表示,并利用均方根逼近得出上述表达式中系数最优值,其中其各项参数取值决定变化速率;所述同一磁感应强度BP对应的极限磁滞回线上升支磁场强度Ha(BP)与P点磁场强度HP差值ΔH(x)表示如下:
ΔH(x)=ΔHrev·(1-b)xe-a(1-x)+ΔHout(BP)·bxc (6)
ΔHout(BP)=Ha(BP)-Hd(BP) (7)
其中,a、b、c为参数,其取值决定ΔH(x)变化速率;ΔHout(BP)为BP对应的极限磁滞回线水平宽度,即磁场强度差值。
10.根据权利要求6或9所述的模拟分析方法,其特征在于:所述方法还包括以下步骤:
选取所得系数最优值,生成直流偏磁条件下铁磁材料一阶回转曲线;
利用生成的一阶回转曲线模拟预测高阶回转曲线,进行磁滞特性分析。
CN201910471820.1A 2019-05-31 2019-05-31 一种直流偏磁条件下铁磁材料的磁滞特性模拟分析***和方法 Active CN110188480B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910471820.1A CN110188480B (zh) 2019-05-31 2019-05-31 一种直流偏磁条件下铁磁材料的磁滞特性模拟分析***和方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910471820.1A CN110188480B (zh) 2019-05-31 2019-05-31 一种直流偏磁条件下铁磁材料的磁滞特性模拟分析***和方法

Publications (2)

Publication Number Publication Date
CN110188480A true CN110188480A (zh) 2019-08-30
CN110188480B CN110188480B (zh) 2020-12-08

Family

ID=67719582

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910471820.1A Active CN110188480B (zh) 2019-05-31 2019-05-31 一种直流偏磁条件下铁磁材料的磁滞特性模拟分析***和方法

Country Status (1)

Country Link
CN (1) CN110188480B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114239299A (zh) * 2021-12-21 2022-03-25 华北电力大学 基于Preisach模型的磁致伸缩确定方法及***
CN114282422A (zh) * 2021-12-30 2022-04-05 哈尔滨工业大学 一种基于局部磁滞回线模型的铁损计算方法
CN117272683A (zh) * 2023-11-13 2023-12-22 江苏南方永磁科技有限公司 一种铁磁材料磁滞效应仿真的模型参数优化方法及***

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120179403A1 (en) * 2009-09-16 2012-07-12 Bae Systems Plc Compensating for hysteresis
US20120207002A1 (en) * 2010-03-01 2012-08-16 Katsnelson Esfir Z Magneto-optical device with an optically induced magnetization
CN102708295A (zh) * 2012-05-14 2012-10-03 广东工业大学 一种电工钢片偏磁特性的分析方法
CN104200038A (zh) * 2014-09-15 2014-12-10 河北工业大学 电工钢片磁滞损耗的Preisach模型优化设计方法
CN104657599A (zh) * 2015-01-30 2015-05-27 国家电网公司 等效微分电(磁)路原理计算直流偏磁单相变压器模型
CN104777384A (zh) * 2015-04-17 2015-07-15 华北电力大学(保定) 直流偏磁状态下变压器磁滞特性及损耗特性确定方法
CN104849587A (zh) * 2015-04-30 2015-08-19 国网四川省电力公司电力科学研究院 一种分析变压器在直流偏磁影响下的励磁特性变化方法
CN106709102A (zh) * 2015-11-17 2017-05-24 黑龙江恒能自控科技有限公司 一种磁悬浮控制力矩陀螺***谐波减速器迟滞建模方法
CN109388883A (zh) * 2018-10-08 2019-02-26 重庆邮电大学 面向emtp的磁通-电流型ja磁滞电感获取方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120179403A1 (en) * 2009-09-16 2012-07-12 Bae Systems Plc Compensating for hysteresis
US20120207002A1 (en) * 2010-03-01 2012-08-16 Katsnelson Esfir Z Magneto-optical device with an optically induced magnetization
CN102708295A (zh) * 2012-05-14 2012-10-03 广东工业大学 一种电工钢片偏磁特性的分析方法
CN104200038A (zh) * 2014-09-15 2014-12-10 河北工业大学 电工钢片磁滞损耗的Preisach模型优化设计方法
CN104657599A (zh) * 2015-01-30 2015-05-27 国家电网公司 等效微分电(磁)路原理计算直流偏磁单相变压器模型
CN104777384A (zh) * 2015-04-17 2015-07-15 华北电力大学(保定) 直流偏磁状态下变压器磁滞特性及损耗特性确定方法
CN104849587A (zh) * 2015-04-30 2015-08-19 国网四川省电力公司电力科学研究院 一种分析变压器在直流偏磁影响下的励磁特性变化方法
CN106709102A (zh) * 2015-11-17 2017-05-24 黑龙江恒能自控科技有限公司 一种磁悬浮控制力矩陀螺***谐波减速器迟滞建模方法
CN109388883A (zh) * 2018-10-08 2019-02-26 重庆邮电大学 面向emtp的磁通-电流型ja磁滞电感获取方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
SAJID HUSSAIN等: "An Efficient Implementation of the Classical Preisach Model", 《IEEE TRANSACTIONS ON MAGNETICS》 *
李慧奇 等: "基于一阶回转曲线的Preisach磁滞模型研究", 《磁性材料及器件》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114239299A (zh) * 2021-12-21 2022-03-25 华北电力大学 基于Preisach模型的磁致伸缩确定方法及***
CN114239299B (zh) * 2021-12-21 2024-05-03 华北电力大学 基于Preisach模型的磁致伸缩确定方法及***
CN114282422A (zh) * 2021-12-30 2022-04-05 哈尔滨工业大学 一种基于局部磁滞回线模型的铁损计算方法
CN117272683A (zh) * 2023-11-13 2023-12-22 江苏南方永磁科技有限公司 一种铁磁材料磁滞效应仿真的模型参数优化方法及***
CN117272683B (zh) * 2023-11-13 2024-01-26 江苏南方永磁科技有限公司 一种铁磁材料磁滞效应仿真的模型参数优化方法及***

Also Published As

Publication number Publication date
CN110188480B (zh) 2020-12-08

Similar Documents

Publication Publication Date Title
CN110188480A (zh) 一种直流偏磁条件下铁磁材料的磁滞特性模拟分析***和方法
AU2013201595B2 (en) Method and apparatus for harmonic state estimation in power system
CN104316894B (zh) 一种电能表实际运行环境的仿真校验方法
CN106443285B (zh) 基于总体最小二乘法的多谐波源谐波责任量化分析方法
CN107317331A (zh) 一种特高压变压器轴对称直流偏磁仿真模型
Yuan et al. Research on dynamic load modelling based on power quality monitoring system
Huang et al. A fast and stable method for modeling generalized nonlinearities in power electronic circuit simulation and its real-time implementation
CN109884402B (zh) 一种三维非对称结构高温超导磁体交流损耗的获取方法
CN103279625A (zh) 不同结构变压器铁心等值电路模型建立方法
Zhang et al. Calculation of DC bias reactive power loss of converter transformer via finite element analysis
Zhu et al. Multi‐objective optimisation design of air‐cored axial flux PM generator
CN108170885A (zh) 一种配电网中多谐波源识别方法
Kontis et al. Dynamic equivalencing of active distribution grids
Chandrasena et al. Simulation of hysteresis and eddy current effects in a power transformer
Jin et al. Modified complementary sliding mode control with disturbance compensation for permanent magnet linear synchronous motor servo system
CN103424621A (zh) 一种谐波电流的人工神经网络检测方法
Ovalle et al. GICharm: A system level analysis tool for geomagnetic disturbance related harmonics
CN109086533B (zh) 一种三相变压器直流偏磁频域计算方法
Shen et al. Deep neural network for magnetic core loss estimation using the magnet experimental database
Basmanov et al. Developing a methodology for calculating asymmetric modes of power supply systems at enterprises if source information is incomplete
Duan et al. Modeling and experimental validation of a dynamic regional saturation JA model for protective current transformer
Jain et al. Adaptive wavelet neural network based harmonic estimation of single-phase systems
CN116911090A (zh) 基于磁导-电容类比法的高频变压器电磁模型建模方法
Shu et al. Novel estimation method of utility harmonic impedance based on short‐term impedance minimum variance criterion
CN115618797A (zh) 基于有限元法的st电磁暂态模型和场路耦合计算方法

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