CN108460220A - 一种基于改进的Jiles-Atherton模型的动态磁滞回线的预测方法 - Google Patents
一种基于改进的Jiles-Atherton模型的动态磁滞回线的预测方法 Download PDFInfo
- Publication number
- CN108460220A CN108460220A CN201810217852.4A CN201810217852A CN108460220A CN 108460220 A CN108460220 A CN 108460220A CN 201810217852 A CN201810217852 A CN 201810217852A CN 108460220 A CN108460220 A CN 108460220A
- Authority
- CN
- China
- Prior art keywords
- atherton
- curve
- jiles
- models
- frequency
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Abstract
本发明公开了一种基于改进的Jiles‑Atherton模型的动态磁滞回线的预测方法。所述的动态磁滞回线预测方法包括如下步骤:通过实验获得的若干组动态磁滞回线的曲线,利用现有的动态Jiles‑Atherton模型获得对应频率下的动态磁滞回线的曲线,获得本发明改进后附加损耗磁场强度Hexcess计算公式的参数,改进现有Jiles‑Atherton模型中有效磁场强度中附加损耗磁场强度Hexcess。本发明的有益之处在于能够预测更准确的任意频率下的动态磁滞回线,改进原有动态Jiles‑Atherton模型预测动态磁滞回线的误差随频率增大而增大的缺点,对铁磁材料铁磁特性的预测提供了更为准确的方法,适用于制造电磁式互感器等涉及非线性铁磁材料瞬态分析的场合。
Description
技术领域
本发明涉及磁滞回线预测领域,特别是涉及一种基于改进的Jiles-Atherton模型的动态磁滞回线的预测方法。
背景技术
铁磁材料的动态磁滞回线是指其在交变磁场磁化下,所得到的B-H关系曲线。描述磁滞现象的理论模型有Preisach(典型如专利CN201410468498.4)、Jiles-Atherton模型(典型如专利CN201710884764.5)等,其中Jiles-Atherton模型因其参数较少、实现方便而被广泛应用在铁磁材料的磁滞建模与仿真中,该模型具有清晰的物理意义,能够真实地描述B-H的非线性关系,通过求解Jiles-Atherton模型方程便能够得到较为准确的B-H曲线。其中Jiles-Atherton模型因其参数较少、实现方便而被广泛应用在铁磁材料的磁滞建模与仿真中,在低频条件下,该模型能够通过求解得到较为准确的B-H曲线,但是在高频条件下,Jiles-Atherton模型的预测B-H曲线误差很大,而且整定参数的算法比较麻烦。所以本发明对Jiles-Atherton模型进行了改良,使其能够更加准确地预测任意频率下的B-H曲线。
发明内容
本发明主要解决的技术问题是Jiles-Atherton模型在动态磁滞回线的预测中误差较随频率增大而增大的问题,在高频情况下预测的B-H曲线误差十分明显。
为解决上述技术问题,本发明采用的一个技术方案是:
一种基于改进的Jiles-Atherton模型的动态磁滞回线的预测方法,包括如下步骤:
1)通过环形铁芯结构的硅钢叠片或者爱泼斯坦方圈测量若干组不同频率下材料的B-H曲线与铁芯损耗;
2)利用步骤1)中得到的铁芯损耗和某一频率下的B-H曲线,用差分进化算法得到现有Jiles-Atherton模型的参数;在步骤1)实验选取的频率下,根据现有Jiles-Atherton模型预测B-H曲线;模型预测的B-H曲线与实验得到的相同频率下的B-H曲线通过做差查找得到对应频率下的最大误差H1与b,其中b为对应频率下磁感应强度随时间变化率的最大值;利用线性回归方法将不同频率下的最大误差H1与b拟合成直线,得斜率k;
3)根据如下公式计算附加损耗磁场强度Hexcess(t):
其中:σ为所测铁磁材料的电导率;
G为一个无量纲系数,大小为0.1356;
S为所测铁磁材料的横截面面积;
H0为额外损耗的系数;
b为磁感应强度随时间变化率的最大值;
λ为方向参数,当λ为+1;当λ为-1;
4)计算总磁场强度Htota;
Htotal=Hhyst(B)+Heddy+Hexcess,其中Heddy为描述瞬态涡流场的磁场强度,Hhyst(B)为磁滞磁场强度。
优选的,步骤2)中所述的某一频率下的B-H曲线是指步骤1)所测得的频率最低的B-H曲线。
本发明还公开了一种基于改进的Jiles-Atherton模型的动态磁滞回线的预测方法,包括如下步骤:
1)通过环形铁芯结构的硅钢叠片或者爱泼斯坦方圈测量材料在两组不同频率下的B-H曲线;
2)利用步骤1)中得到较低频率下的B-H曲线,用差分进化算法得到现有Jiles-Atherton模型的参数;在步骤1)实验选取的频率下,根据现有Jiles-Atherton模型预测另一频率下的B-H曲线,与实验得到的该频率下的B-H曲线通过做差查找得到最大误差H1;
3)根据如下公式计算附加损耗磁场强度Hexcess(t):
其中:σ为所测铁磁材料的电导率;
G为一个无量纲系数,大小为0.1356;
S为所测铁磁材料的横截面面积;
H0为额外损耗的系数;
b为磁感应强度随时间变化率的最大值;
λ为方向参数,当λ为+1;当λ为-1;
4)计算总磁场强度Htota;
Htotal=Hhyst(B)+Heddy+Hexcess,其中Heddy为描述瞬态涡流场的磁场强度,Hhyst(B)为磁滞磁场强度。
本发明的有益效果是:
1、Hexcess(t)更准确
2、可根据若干组低频下的铁芯损耗和动态磁滞回线预测任意频率下的B-H曲线。
3、预测的高频B-H曲线比原模型准确度更高,误差更小(更接近材料实际的动态B-H曲线)。
附图说明
图1计算改进Jiles-Atherton模型的参数及预测B-H曲线的流程图;
图2为实施例的线性拟合结果图;
图3. 500Hz、0.95T下预测的B-H曲线与实验B-H曲线的对比图;
图4. 0.85T下改进Jiles-Atherton模型预测的B-H曲线与实验B-H曲线的对比图。
具体实施方式
下面对本发明的较佳实施例进行详细阐述,以使本发明的优点和特征能更易于被本领域技术人员理解,从而对本发明的保护范围做出更为清楚明确的界定。
下面对本发明的较佳实施例进行详细阐述,以使本发明的优点和特征能更易于被本领域技术人员理解,从而对本发明的保护范围做出更为清楚明确的界定。
本发明实施例如下
1)采用50WW800材料,通过环形铁芯结构的硅钢叠片或者爱泼斯坦方圈测量材料在幅值为0.95T,频率分别为50Hz、100Hz、150Hz、…、450Hz、500Hz下的B-H曲线,并得到铁芯损耗;
2)利用步骤1)中得到的铁芯损耗和50Hz下的B-H曲线,用差分进化算法得到现有Jiles-Atherton模型的参数;在步骤1)实验选取的频率下,根据现有Jiles-Atherton模型预测B-H曲线;模型预测的B-H曲线与实验得到的相同频率下的B-H曲线通过做差查找得到对应频率下的最大误差H1与b,其中b为对应频率下磁感应强度随时间变化率的最大值;利用线性回归方法将不同频率下的最大误差H1与b拟合成直线,得斜率k=0.0226;线性拟合结果如图2所示,
3)根据如下公式计算附加损耗磁场强度Hexcess(t):
其中:σ为所测铁磁材料的电导率;
G为一个无量纲系数,大小为0.1356;
S为所测铁磁材料的横截面面积;
H0为额外损耗的系数;
b为磁感应强度随时间变化率的最大值;
λ为方向参数,当λ为+1;当λ为-1;
4)计算总磁场强度Htota;
Htotal=Hhyst(B)+Heddy+Hexcess,其中Heddy为描述瞬态涡流场的磁场强度,Hhyst(B)为磁滞磁场强度。
在500Hz下,分别用现有JA模型和本方法预测B-H曲线,其结果如图3所示,图3说明了原始Jiles-Atherton模型预测B-H曲线精度较低而改进的Jiles-Atherton模型预测精度较高。
0.85T下,用本方法预测50Hz、100Hz、150Hz下的B-H曲线与实验值进行比较,结果如图4所示,图4说明了在不同频率下本方法预测出的B-H曲线与实验值接近,具有较高精确度。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (3)
1.一种基于改进的Jiles-Atherton模型的动态磁滞回线的预测方法,其特征在于包括如下步骤:
1)通过环形铁芯结构的硅钢叠片或者爱泼斯坦方圈测量若干组不同频率下材料的B-H曲线与铁芯损耗;
2)利用步骤1)中得到的铁芯损耗和某一频率下的B-H曲线,用差分进化算法得到现有Jiles-Atherton模型的参数;在步骤1)实验选取的频率下,根据现有Jiles-Atherton模型预测B-H曲线;模型预测的B-H曲线与实验得到的相同频率下的B-H曲线通过做差查找得到对应频率下的最大误差H1与b,其中b为对应频率下磁感应强度随时间变化率的最大值;利用线性回归方法将不同频率下的最大误差H1与b拟合成直线,得斜率k;
3)根据如下公式计算附加损耗磁场强度Hexcess(t):
其中:σ为所测铁磁材料的电导率;
G为一个无量纲系数,大小为0.1356;
S为所测铁磁材料的横截面面积;
H0为额外损耗的系数;
b为磁感应强度随时间变化率的最大值;
λ为方向参数,当λ为+1;当λ为-1;
4)计算总磁场强度Htota;
Htotal=Hhyst(B)+Heddy+Hexcess,其中Heddy为描述瞬态涡流场的磁场强度,Hhyst(B)为磁滞磁场强度。
2.根据权利要求1所述的基于改进的Jiles-Atherton模型的动态磁滞回线的预测方法,其特征在于步骤2)中所述的某一频率下的B-H曲线是指步骤1)所测得的频率最低的B-H曲线。
3.一种基于改进的Jiles-Atherton模型的动态磁滞回线的预测方法,其特征在于包括如下步骤:
1)通过环形铁芯结构的硅钢叠片或者爱泼斯坦方圈测量材料在两组不同频率下的B-H曲线;
2)利用步骤1)中得到较低频率下的B-H曲线,用差分进化算法得到现有Jiles-Atherton模型的参数;在步骤1)实验选取的频率下,根据现有Jiles-Atherton模型预测另一频率下的B-H曲线,与实验得到的该频率下的B-H曲线通过做差查找得到最大误差H1;
3)根据如下公式计算附加损耗磁场强度Hexcess(t):
其中:σ为所测铁磁材料的电导率;
G为一个无量纲系数,大小为0.1356;
S为所测铁磁材料的横截面面积;
H0为额外损耗的系数;
b为磁感应强度随时间变化率的最大值;
λ为方向参数,当λ为+1;当λ为-1;
4)计算总磁场强度Htota;
Htotal=Hhyst(B)+Heddy+Hexcess,其中Heddy为描述瞬态涡流场的磁场强度,Hhyst(B)为磁滞磁场强度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810217852.4A CN108460220B (zh) | 2018-03-16 | 2018-03-16 | 一种基于改进的Jiles-Atherton模型的动态磁滞回线的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810217852.4A CN108460220B (zh) | 2018-03-16 | 2018-03-16 | 一种基于改进的Jiles-Atherton模型的动态磁滞回线的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108460220A true CN108460220A (zh) | 2018-08-28 |
CN108460220B CN108460220B (zh) | 2020-12-25 |
Family
ID=63236795
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810217852.4A Active CN108460220B (zh) | 2018-03-16 | 2018-03-16 | 一种基于改进的Jiles-Atherton模型的动态磁滞回线的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108460220B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109884564A (zh) * | 2019-03-22 | 2019-06-14 | 华中科技大学 | 一种变压器磁芯磁特性测量方法和装置 |
CN109933914A (zh) * | 2019-03-18 | 2019-06-25 | 北京工业大学 | 双相铁磁性材料磁滞及巴克豪森噪声信号的建模方法 |
CN111044956A (zh) * | 2019-11-28 | 2020-04-21 | 浙江大学 | 一种磁滞损耗的估算方法 |
CN113049998A (zh) * | 2021-02-08 | 2021-06-29 | 华北电力大学(保定) | 多谐波激励作用下铁磁材料损耗预测方法 |
CN114239299A (zh) * | 2021-12-21 | 2022-03-25 | 华北电力大学 | 基于Preisach模型的磁致伸缩确定方法及系统 |
CN114236433A (zh) * | 2021-11-23 | 2022-03-25 | 浙江大学 | 一种铁磁材料磁感应强度在线估测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7502723B1 (en) * | 2005-06-01 | 2009-03-10 | Linear Technology Corporation | Asymmetric minor hysteresis loop model and circuit simulator including the same |
CN107681872A (zh) * | 2017-09-30 | 2018-02-09 | 厦门理工学院 | 一种抑制开关电源温漂的方法 |
-
2018
- 2018-03-16 CN CN201810217852.4A patent/CN108460220B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7502723B1 (en) * | 2005-06-01 | 2009-03-10 | Linear Technology Corporation | Asymmetric minor hysteresis loop model and circuit simulator including the same |
CN107681872A (zh) * | 2017-09-30 | 2018-02-09 | 厦门理工学院 | 一种抑制开关电源温漂的方法 |
Non-Patent Citations (3)
Title |
---|
宋丽娜: "基于智能优化算法磁电复合薄膜特性的研究", 《中国优秀硕士学位论文全文数据库·基础科学辑》 * |
李泽涛: "基于J_A模型的单相变压器建模和偏磁研究", 《中国优秀硕士学位论文全文数据库·工程科技Ⅱ辑》 * |
杨延菊: "基于神经网络的Jiles-Atherton磁滞模型的实现", 《中国优秀硕士学位论文全文数据库·工程科技Ⅱ辑》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109933914A (zh) * | 2019-03-18 | 2019-06-25 | 北京工业大学 | 双相铁磁性材料磁滞及巴克豪森噪声信号的建模方法 |
CN109933914B (zh) * | 2019-03-18 | 2024-01-26 | 北京工业大学 | 双相铁磁性材料磁滞及巴克豪森噪声信号的建模方法 |
CN109884564A (zh) * | 2019-03-22 | 2019-06-14 | 华中科技大学 | 一种变压器磁芯磁特性测量方法和装置 |
CN109884564B (zh) * | 2019-03-22 | 2020-02-14 | 华中科技大学 | 一种变压器磁芯磁特性测量方法和装置 |
CN111044956A (zh) * | 2019-11-28 | 2020-04-21 | 浙江大学 | 一种磁滞损耗的估算方法 |
CN111044956B (zh) * | 2019-11-28 | 2021-05-11 | 浙江大学 | 一种磁滞损耗的估算方法 |
CN113049998A (zh) * | 2021-02-08 | 2021-06-29 | 华北电力大学(保定) | 多谐波激励作用下铁磁材料损耗预测方法 |
CN113049998B (zh) * | 2021-02-08 | 2021-11-02 | 华北电力大学(保定) | 多谐波激励作用下铁磁材料损耗预测方法 |
CN114236433A (zh) * | 2021-11-23 | 2022-03-25 | 浙江大学 | 一种铁磁材料磁感应强度在线估测方法 |
CN114236433B (zh) * | 2021-11-23 | 2022-08-30 | 浙江大学 | 一种铁磁材料磁感应强度在线估测方法 |
CN114239299A (zh) * | 2021-12-21 | 2022-03-25 | 华北电力大学 | 基于Preisach模型的磁致伸缩确定方法及系统 |
CN114239299B (zh) * | 2021-12-21 | 2024-05-03 | 华北电力大学 | 基于Preisach模型的磁致伸缩确定方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN108460220B (zh) | 2020-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108460220A (zh) | 一种基于改进的Jiles-Atherton模型的动态磁滞回线的预测方法 | |
Rezaei-Zare | Enhanced transformer model for low-and mid-frequency transients—Part I: Model development | |
CN107656221A (zh) | 一种基于局部磁滞回线斜率的变压器铁芯剩磁估算方法 | |
Prus et al. | Research of rnergy processes in circuits containing iron in saturation condition | |
CN109444776A (zh) | 三相三柱变压器铁芯剩磁测算方法、系统及存储介质 | |
Hussain et al. | The modified Jiles–Atherton model for the accurate prediction of iron losses | |
CN109884564A (zh) | 一种变压器磁芯磁特性测量方法和装置 | |
Szewczyk | Technical BH saturation magnetization curve models for SPICE, FEM and MoM simulations | |
Sato et al. | Loss computation of soft magnetic composite inductors based on interpolated scalar magnetic property | |
Ducharne et al. | Anomalous fractional magnetic field diffusion through cross-section of a massive toroidal ferromagnetic core | |
Li et al. | Anomalous loss modeling and validation of magnetic materials in electrical engineering | |
Stenglein et al. | Core loss model for arbitrary excitations with DC bias covering a wide frequency range | |
CN102708295B (zh) | 一种电工钢片偏磁特性的分析方法 | |
Sarker et al. | A generalized inverse Preisach dynamic hysteresis model of Fe-based amorphous magnetic materials | |
Zirka et al. | Measurement and Modeling of $ B $–$ H $ Loops and Losses of High Silicon Nonoriented Steels | |
Yamagashira et al. | Vector magnetic properties and 2-D magnetostriction of various electrical steel sheets under rotating flux condition | |
Shimizu et al. | Iron loss eveluation of filter inductor used in PWM inverters | |
Bui et al. | Characterization and modeling of a current transformer working under thermal stress | |
Lu et al. | Dynamic circuit modeling of a high frequency transformer | |
Zhao et al. | A Dynamic Model for Prediction of Hysteresis and Loss Property of Soft Material Under Multiple Harmonic Magnetizations | |
CN107016172A (zh) | Buck变换器磁芯损耗计算方法 | |
Wilk et al. | Simulation of the remanence influence on the transient states of the single-phase transformer including feedback Preisach model | |
Najgebauer | Fractional scaling of magnetic coercivity in electrical steels | |
Tu et al. | Study on excess loss of soft magnetic material considering skin effect | |
CN107688679A (zh) | 一种硅钢片铁芯宽频电路模型的建模方法 |
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 |