本发明公开了一种电工钢片磁致伸缩模型的建模方法，包括：测量不同旋转磁场下电工钢片磁致伸缩对称蝴蝶曲线εk；对所测量的对称蝴蝶曲线进行辨识，构建磁致伸缩蝴蝶曲线εk的Everett函数数据库；将输入的矢量磁场B波形用旋转磁场形状参数Bmax,最大值方向的方向角θ和椭圆长短轴的长度比λ拟合，得到拟合后结果；通过线性插值法在磁致伸缩Everett函数数据库中进行插值，获得拟合后结果所对应的磁致伸缩Everett函数值Ek；根据Everett函数值Ek计算得出对应的磁致伸缩蝴蝶曲线εk。应用本发明实施例，不仅可以描述矢量磁致伸缩特性，还可以通过计算估算出旋转磁场下的瞬态磁致伸缩。此外，还能够考虑磁致伸缩的各向异性特性。
基本信息
申请号：CN202111133151.0
申请日期：20210927
公开号：CN202111133151.0
公开日期：20211207
申请人：上海海事大学
申请人地址：201306 上海市浦东新区临港新城海港大道1550号
发明人：朱立勋;刘仕一;高宁;吴卫民
当前权利人：上海海事大学
代理机构：上海互顺专利代理事务所(普通合伙) 31332
代理人：成秋丽
主权利要求
1.一种电工钢片磁致伸缩模型的建模方法，其特征在于，包括：测量不同旋转磁场下电工钢片磁致伸缩对称蝴蝶曲线εk；对所测量的对称蝴蝶曲线进行辨识，构建磁致伸缩蝴蝶曲线εk的Everett函数数据库；将输入的矢量磁场B波形用旋转磁场形状参数Bmax,最大值方向的方向角θ和椭圆长短轴的长度比λ拟合，得到拟合后结果；通过线性插值法在磁致伸缩Everett函数数据库中进行插值，获得拟合后结果所对应的磁致伸缩Everett函数值Ek；根据Everett函数值Ek计算得出对应的磁致伸缩蝴蝶曲线εk。
权利要求
1.一种电工钢片磁致伸缩模型的建模方法，其特征在于，包括：
测量不同旋转磁场下电工钢片磁致伸缩对称蝴蝶曲线ε k；
对所测量的对称蝴蝶曲线进行辨识，构建磁致伸缩蝴蝶曲线ε k的Everett函数数据库；
将输入的矢量磁场B波形用旋转磁场形状参数B max,最大值方向的方向角θ和椭圆长短轴的长度比λ拟合，得到拟合后结果；
通过线性插值法在磁致伸缩Everett函数数据库中进行插值，获得拟合后结果所对应的磁致伸缩Everett函数值E k；
根据Everett函数值E k计算得出对应的磁致伸缩蝴蝶曲线εk。
2.根据权利要求1所述的一种电工钢片磁致伸缩模型的建模方法，其特征在于，对所测量的对称蝴蝶曲线进行辨识，构建磁致伸缩蝴蝶曲线ε k的Everett函数数据库的步骤为：
在磁感应强度最大值B max为B a的一旋转磁场下，测量对称蝴蝶曲线；
计算矢量磁场点为(B a,-B a)的磁致伸缩Everett函数为第一函数值εa；
利用测量曲线计算其他点(B a,B b)点的磁致伸缩Everett函数为第二函数值(εa-εb)/2，B b为不同于-B a的其他值，εa、εb是测量对称蝴蝶曲线上的值；
通过磁致伸缩Everett函数关于直线α＝-β的对称性，获得对称点的函数值，并且α为磁致伸缩算子的正向翻转阈值，β为磁致伸缩算子的负向翻转阈值；
分别计算每条测量的蝴蝶回线对应的磁致伸缩Everett函数值E k,并基于所计算的函数值建立磁致伸缩Everett函数数据库。
3.根据权利要求1所述的一种电工钢片磁致伸缩模型的建模方法，其特征在于，根据Everett函数值E k计算得出对应的磁致伸缩蝴蝶曲线εk的公式为：
其中，B t为输入波形极值存储器中序号为t的极值序列，θ为最大值方向的方向角，λ为椭圆长短轴的长度比，N表示t的最大值。
4.根据权利要求1所述的一种电工钢片磁致伸缩模型的建模方法，其特征在于，基于Preisach模型的磁致伸缩模型表示为：
ε k(t)＝∫∫ α≥βμ k(θ,λ,α,β)γ αβ(B k(t))dαdβ
其中，γ αβ为磁致伸缩算子，α为其正向翻转阈值，β为其负向翻转阈值，θ为最大值方向的方向角，λ为椭圆长短轴的长度比；μ k(θ,λ,α,β)是输入旋转磁场形状参数拟合成规则椭圆形后形状参数θ和λ的函数，B k(t)为t时刻k方向上的波形。
说明书
一种电工钢片磁致伸缩模型的建模方法
技术领域
本发明涉及磁性材料特性建模技术领域，尤其涉及一种电工钢片磁致伸缩模型的建模方法。
背景技术
众所周知，电机(包括旋转电机和变压器)中的大部分振动和噪声是由电工钢片中的磁致伸缩引起的，尤其是变压器中的电工钢片。因此，开发精确的磁致伸缩模型有助于降低电机和变压器设计状态下的振动和噪声。然而，电机和变压器的不同工作状态所对应的各种不同的磁化条件下，如交替和旋转磁化条件下，磁致伸缩特性是相当复杂的。
迄今为止，文献IEEE Transactions on Magnetic,56(3):7514604,2020(Animproved magnetostriction model for electrical steel sheet based on Jiles–Atherton model)介绍了基于Jiles-Atherton模型改进的电工钢片磁致伸缩模型，文献IEEE Transactions on Magnetic,33(5):pp.3931-3933,1997(Utilizing anisotropicPreisach-type models in the accurate simulation of magnetostriction)介绍了在磁致伸缩的精确模拟中利用的各向异性Preisach型模型，文献IEEE Transactions onMagnetic,51(5):6000407,2015(Modeling 2-Dmagnetostriction in nonorientedelectrical steels using a simple magnetic domain model)介绍了模拟无取向电工钢的二维磁致伸缩的磁畴模型，文献International Conference on Electrical Machines(ICEM),2014(Two-dimensional magnetostriction analysis using E&S-W model ininduction motor model core)介绍了在感应电机模型铁芯中进行二维磁致伸缩分析的E&S-W模型。在上述模型中，由于Preisach模型具有精度高、易于实现、具有物理意义等优点，在磁滞领域得到了广泛的研究。由于磁致伸缩特性被认为与磁滞特性相似，因此已经开发了一些基于Preisach模型的磁致伸缩模型。而且有的只能模拟磁致伸缩的最大值和最小值，有的只能模拟交变磁场下的磁致伸缩。然而，旋转磁场分布在电机中，在分析瞬态振动和噪声特性时需要瞬态磁致伸缩特性。上述几种模型不能同时满足模拟矢量磁致伸缩和估算瞬态磁致伸缩的需求。
发明内容
本发明的目的在于提供一种电工钢片磁致伸缩模型的建模方法，不仅可以描述矢量磁致伸缩特性，还可以通过计算估算出旋转磁场下的瞬态磁致伸缩。
为了实现上述目的，本发明提供一种电工钢片磁致伸缩模型的建模方法，包括：
测量不同旋转磁场下电工钢片磁致伸缩对称蝴蝶曲线ε k；
对所测量的对称蝴蝶曲线进行辨识，构建磁致伸缩蝴蝶曲线ε k的Everett函数数据库；
将输入的矢量磁场B波形用旋转磁场形状参数B max,最大值方向的方向角θ和椭圆长短轴的长度比λ拟合，得到拟合后结果；
通过线性插值法在磁致伸缩Everett函数数据库中进行插值，获得拟合后结果所对应的磁致伸缩Everett函数值E k；
根据Everett函数值E k计算得出对应的磁致伸缩蝴蝶曲线εk。
可选的，对所测量的对称蝴蝶曲线进行辨识，构建磁致伸缩蝴蝶曲线ε k的Everett函数数据库的步骤为：
在磁感应强度最大值B max为B a的一旋转磁场(θ 0,λ 0)下，测量对称蝴蝶曲线；
计算矢量磁场点为(B a,-B a)的磁致伸缩Everett函数为第一函数值εa；
利用测量曲线计算其他点(B a,B b)点的磁致伸缩Everett函数为第二函数值(εa-εb)/2，B b为不同于-B a的其他值；其中，εa、εb分别是B＝B a、B＝B b时的磁致伸缩值，εa、εb是测量对称蝴蝶曲线上的值，蝴蝶曲线纵轴是磁致伸缩值ε，横轴是磁感应强度B；
通过磁致伸缩Everett函数关于直线α＝-β的对称性，获得对称点的函数值，并且α为磁致伸缩算子的正向翻转阈值，β为磁致伸缩算子的负向翻转阈值；
分别计算每条测量的蝴蝶回线对应的Everett函数值E k,并基于所计算的函数值建立磁致伸缩Everett函数数据库。
一种实现方式中，根据Everett函数值E k计算得出对应的磁致伸缩蝴蝶曲线εk的公式为：
其中，B t为输入波形极值存储器中的极值序列，θ为最大值方向的方向角，λ为椭圆长短轴的长度比。
一种实现方式中，基于Preisach模型的磁致伸缩模型表示为：
ε k(t)＝∫∫ α≥βμ k(θ,λ,α,β)γ αβ(B k(t))dαdβ
其中，γ αβ为磁致伸缩算子，并且α为其正向翻转阈值，β为其负向翻转阈值，θ为最大值方向的方向角，λ为椭圆长短轴的长度比；μ k(θ,λ,α,β)是输入旋转磁场形状参数拟合成规则椭圆形后形状参数θ和λ的函数，B k(t)为t时刻k方向上的波形。
应用本发明实施例提供的一种电工钢片磁致伸缩模型的建模方法，具备的有益效果如下：
(1)建立的新型磁致伸缩模型，可以计算不同磁场下电工钢片磁致伸缩特性随时间变化的瞬时值。
(2)通过建立磁致伸缩Everett函数数据库，简化了Preisach模型的计算过程。
附图说明
图1是本发明实施例一种电工钢片磁致伸缩模型的建模方法的一种流程示意图。
图2a是35PH105硅钢片在交变磁场下：当θ＝20°时，不同B max的实测主磁致伸缩应变ε x的蝴蝶曲线。
图2b是35PH105硅钢片在交变磁场下：当B max＝1.0T时，不同θ的实测主磁致伸缩应变ε x的蝴蝶曲线。
图3a是在旋转磁场下，当θ＝20°时，不同B max的实测主磁致伸缩应变ε x的蝴蝶曲线。
图3b是在旋转磁场下，当B max＝1.0T时，不同θ的实测主磁致伸缩应变ε x的蝴蝶曲线。
图4是当B max＝B a时的对称磁致伸缩蝴蝶曲线示意图。
图5a为ε x和B x的Everett函数：m＝ε x,n＝B x。
图5b为ε y和B y的Everett函数：m＝ε y,n＝B y。
图5c为ε xy和B x 2+B y 2的Everett函数：m＝ε xy,n＝B x 2+B y 2。
图6是为磁滞算子示意图。
图7为输入B波形拟合成旋转磁场椭圆轨迹示意图。
图8为Preisach平面示意图。
图9为本发明线性插值法的插值过程。
图10为极值存储器存储输入矢量B波形的极值过程。
图11a、图11b、图11c、图11d为35PH105硅钢片在不同输入B波形下模拟和测量的磁致伸缩结果的比较。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式，本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本发明的精神下进行各种修饰或改变。
如图1本发明提供一种电工钢片磁致伸缩模型的建模方法，包括：
S110，测量不同旋转磁场下电工钢片磁致伸缩对称蝴蝶曲线ε k。
需要说明的是，在圆形二维单片应变仪里测量电工钢片在多种矢量磁场条件下的磁致伸缩ε k，输入矢量磁场的x轴和y轴分量,可以通过圆形二维单片应变仪里固定在定子上的两个相线圈产生，通过控制线圈的电流可以产生任意旋转磁场。通过测量得到的一部分磁致伸缩蝴蝶曲线如图2a、图2b、图3a、图3b，图中的磁滞特性曲线是基于韩国浦项钢铁公司型号为35PH105的硅钢片测量得到的。
具体的，在交变磁场下测量得到：如图2a所示，是当θ＝20°时，不同B max的实测主磁致伸缩应变ε x的蝴蝶曲线；如图2b所示，是当B max＝1.0T时，不同θ的实测主磁致伸缩应变ε x的蝴蝶曲线。在旋转磁场下测量得到：如图3a所示，是当θ＝20°时，不同B max的实测主磁致伸缩应变ε x的蝴蝶曲线；如图3b所示，是当B max＝1.0T时，不同θ的实测主磁致伸缩应变ε x的蝴蝶曲线。
S120，对所测量的对称蝴蝶曲线进行辨识，构建磁致伸缩蝴蝶曲线ε k的Everett函数数据库。
可以理解的是，磁致伸缩Everett函数数据库的构建需要依靠测量的对称蝴蝶曲线进行辨识，其辨识步骤如下：
(1)在磁感应强度最大值B max为矢量磁场点为B a的某一给定旋转磁场(θ 0,λ 0)下测量如图4所示的对称蝴蝶曲线；其中，θ为最大值方向的方向角，λ为椭圆长短轴的长度比。θ 0λ 0表示给定的确定值。
(2)计算点(B a,-B a)的磁致伸缩Everett函数为E(θ 0,λ 0,B a,-B a)＝ε a；
(3)利用测量曲线计算其他点(B a,B b)点的磁致伸缩Everett函数为E(θ 0,λ 0,B a,B b)＝(ε a-ε b)/2；
(4)利用磁致伸缩Everett函数关于直线α＝-β的对称性，得到对称点的函数值。
(5)重复上述步骤1～4，分别计算每条测量的蝴蝶回线对应的函数值E k(θ,λ,B,B),并建立磁致伸缩Everett函数数据库。
S130，将输入的矢量磁场B波形用旋转磁场形状参数B max,最大值方向的方向角θ和椭圆长短轴的长度比λ拟合，得到拟合后结果。
需要说明的是，本发明使用的拟合方法为最小二乘近似法：
θ＝∠(B x(τ max),B y(τ max))
其中，τ max表示矢量B波形最大时的方向角，B x和B y分别是输入矢量B波形的x方向和y方向的分量，τ i表示第i时刻的方向角，B max表示B波形的最大值，τ表示某个时刻的方向角，B是输入的矢量磁场B波形，为非椭圆波形；b是椭圆波形，N B是输入的矢量磁场B波形的采样点数。τ＝ωt，ω表示角速度，t表示时刻，τ表示某个时刻的方向角。
需要说明的是，将输入的B波形拟合成规则椭圆形后形状参数θ和λ的函数，其中θ为最大值方向的方向角，λ为椭圆长短轴的长度比，如图5a、图5b、图5c是施加在35PH105硅钢片上的旋转磁场为(θ＝30°,λ＝0.3)时，Everett函数数据库的结果：图5a是ε x和B x的Everett函数,m＝ε x,n＝B x；如图5b是ε y和B y的Everett函数，m＝ε y,n＝B y；如图5c是ε xy和B x 2+B y 2的Everett函数，m＝ε xy,n＝B x 2+B y 2。
S140，通过线性插值法在磁致伸缩Everett函数数据库中进行插值，获得拟合后结果所对应的磁致伸缩Everett函数值E k。
插值方法如图9所示，已知Everett函数值E(Q ij)为点Q ij(θ i,λ j)的Everett函数值，其中，θ i参数θ第i个数，λ j为参数λ的第j个数；
(1)在θ轴方向上进行插值求得点Q a1(θ a,λ 1)的Everett函数值：
其中，θ R,λ 1，Q 11、Q 21、θ 1,θ 2、E(Q 11)、E(Q 21)见图示。
(2)在θ轴方向上进行插值求得点Q a2(θ a,λ 2)的Everett函数值：
其中，λ 2，Q 12、Q 22、E(Q 12)、E(Q 22)见图示。
(3)在λ轴方向上进行插值求得点Q ab(θ a,λ b)的Everett函数值：
其中，λ R，λ 1、λ 2、R 1,R 2、E(R 2)、E(R 1)见图示。
通过上述线性插值法在磁致伸缩Everett函数数据库中进行插值，求得拟合后结果所对应的磁致伸缩Everett函数值E k。
S150，根据Everett函数值E k计算得出对应的磁致伸缩蝴蝶曲线ε k。
需要说明的是，本发明基于Preisach模型理论，提出一种可以考虑相应矢量磁场下的主磁致伸缩和剪切磁致伸缩的磁致伸缩模型。基于Preisach模型的磁致伸缩模型表示为：
ε k(t)＝∫∫ α≥βμ k(θ,λ,α,β)γ αβ(B k(t))dαdβ (1)
式(1)中，下标k分别表示x,y和xy，则ε
x和ε
y分别表示x和y方向的主磁致伸缩，ε
xy表示剪切磁致伸缩。B
x和B
y分别是输入B波形的x方向和y方向的分量，并且
γ
αβ为磁致伸缩算子，并且α为其正向+1翻转阈值，β为其负向-1翻转阈值，如图6所示，定义如下：
式(2)中，t i是输入B波形的时间步长，μ x、μ y、μ xy分别为主磁致伸缩和剪切磁致伸缩算子对应的分布函数，该分布函数根据不同磁场激励下的磁致伸缩测量曲线计算得到。μ x、μ y和μ xy不仅是阈值α和β的函数，而且是输入B波形拟合成规则椭圆形后形状参数θ和λ的函数，其中θ为最大值方向的方向角，λ为椭圆长短轴的长度比，如图7所示。
从Preisach磁致伸缩模型的定义式(1)可知，式中的积分区域是一个三角形区域，称为Preisach平面，如图8所示。
在Preisach平面中，每个点都代表一种磁滞算子。根据输入的矢量磁场，历史依赖阶梯线L(t)将Preisach平面分为T +(t)和T -(t)部分，其磁滞算子的值分别为+1和-1。所以Preisach磁致伸缩模型的定义式(1)可以表示为：
本发明所提出的Preisach磁致伸缩模型，引入磁致伸缩Everett函数E k(θ,λ,α,β)，用来代替分布函数μ k。在给定磁场(θ 0,λ 0)，它所对应的定义式的积分区域，如图8所示，在这个区域里给定一点(α 0,β 0)((α 0,β 0)是Preisach平面上给定的一个点)的Everett函数定义如下：
通过式(3)和式(4)，可以把Preisach磁致伸缩模型的定义式(1)的复杂计算转换成如下公式:
其中，B 1为输入波形极值存储器中序号为1的极值序列、B 2为输入波形极值存储器中序号为2的极值序列、B t为输入波形极值存储器中序号为t的极值序列，N表示t的最大值。极值存储器存储输入矢量B波形的极值过程，如图10所示。其大小排序如下：
|B 1|＞|B 3|＞···＞|B 2i-1|
|B 2|＞|B 4|＞···＞|B 2i|
|B 1|＞|B 2|
经过上述计算，可以求得定义式(1)中Preisach平面上的Everett函数值E k，从而得出对应的磁致伸缩蝴蝶曲线εk。
图11a为交变磁场B max＝1.0T,θ＝80°,λ＝0，图11b为旋转磁场B max＝1.0T,θ＝30°,λ＝0.3，图11c为旋转磁场B max＝1.0T,θ＝40°,λ＝0.1，图11d为旋转磁场B max＝1.0T,θ＝40°,λ＝0.2，通过比较了型号为35PH105硅钢片在不同磁场下基于本发明模型得到的磁致伸缩曲线和测量结果，其中虚线代表模型计算曲线，实线代表测量曲线。结果表明本发明提出的模型结果与试验结果基本吻合，本模型可以预测出与实际情况接近的磁致伸缩特性。
上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下，对上述实施例进行修饰或改变。因此，举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变，仍应由本发明的权利要求所涵盖。
