本发明涉及磁滞回线预测领域，特别是涉及一种基于改进的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为对应频率下磁感应强度随时间变化率
3)根据如下公式计算附加损耗磁场强度hexcess(t)：
其中：σ为所测铁磁材料的电导率；
g为一个无量纲系数，大小为0.1356；
s为所测铁磁材料的横截面面积；
h0为额外损耗的系数；
b为磁感应强度随时间变化率
λ为方向参数，当
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为磁感应强度随时间变化率
λ为方向参数，当
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为对应频率下磁感应强度随时间变化率
3)根据如下公式计算附加损耗磁场强度hexcess(t)：
其中：σ为所测铁磁材料的电导率；
g为一个无量纲系数，大小为0.1356；
s为所测铁磁材料的横截面面积；
h0为额外损耗的系数；
b为磁感应强度随时间变化率
λ为方向参数，当
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曲线与实验值接近，具有较高精确度。
以上所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。