CN113761691B - 一种电工钢片磁致伸缩模型的建模方法 - Google Patents
一种电工钢片磁致伸缩模型的建模方法 Download PDFInfo
- Publication number
- CN113761691B CN113761691B CN202111133151.0A CN202111133151A CN113761691B CN 113761691 B CN113761691 B CN 113761691B CN 202111133151 A CN202111133151 A CN 202111133151A CN 113761691 B CN113761691 B CN 113761691B
- Authority
- CN
- China
- Prior art keywords
- magnetostriction
- everett
- epsilon
- function
- magnetostrictive
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 229910000976 Electrical steel Inorganic materials 0.000 title claims abstract description 25
- 230000006870 function Effects 0.000 claims description 57
- 238000005259 measurement Methods 0.000 claims description 9
- 230000006698 induction Effects 0.000 claims description 6
- 230000001052 transient effect Effects 0.000 abstract description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000005315 distribution function Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005381 magnetic domain Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- 229910000565 Non-oriented electrical steel Inorganic materials 0.000 description 1
- 102100029469 WD repeat and HMG-box DNA-binding protein 1 Human genes 0.000 description 1
- 101710097421 WD repeat and HMG-box DNA-binding protein 1 Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 239000000696 magnetic material Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Soft Magnetic Materials (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明公开了一种电工钢片磁致伸缩模型的建模方法,包括:测量不同旋转磁场下电工钢片磁致伸缩对称蝴蝶曲线εk;对所测量的对称蝴蝶曲线进行辨识,构建磁致伸缩蝴蝶曲线εk的Everett函数数据库;将输入的矢量磁场B波形用旋转磁场形状参数Bmax,最大值方向的方向角θ和椭圆长短轴的长度比λ拟合,得到拟合后结果;通过线性插值法在磁致伸缩Everett函数数据库中进行插值,获得拟合后结果所对应的磁致伸缩Everett函数值Ek;根据Everett函数值Ek计算得出对应的磁致伸缩蝴蝶曲线εk。应用本发明实施例,不仅可以描述矢量磁致伸缩特性,还可以通过计算估算出旋转磁场下的瞬态磁致伸缩。此外,还能够考虑磁致伸缩的各向异性特性。
Description
技术领域
本发明涉及磁性材料特性建模技术领域,尤其涉及一种电工钢片磁致伸缩模型的建模方法。
背景技术
众所周知,电机(包括旋转电机和变压器)中的大部分振动和噪声是由电工钢片中的磁致伸缩引起的,尤其是变压器中的电工钢片。因此,开发精确的磁致伸缩模型有助于降低电机和变压器设计状态下的振动和噪声。然而,电机和变压器的不同工作状态所对应的各种不同的磁化条件下,如交替和旋转磁化条件下,磁致伸缩特性是相当复杂的。
迄今为止,文献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波形用旋转磁场形状参数Bmax,最大值方向的方向角θ和椭圆长短轴的长度比λ拟合,得到拟合后结果;
通过线性插值法在磁致伸缩Everett函数数据库中进行插值,获得拟合后结果所对应的磁致伸缩Everett函数值Ek;
根据Everett函数值Ek计算得出对应的磁致伸缩蝴蝶曲线εk。
可选的,对所测量的对称蝴蝶曲线进行辨识,构建磁致伸缩蝴蝶曲线εk的Everett函数数据库的步骤为:
在磁感应强度最大值Bmax为Ba的一旋转磁场(θ0,λ0)下,测量对称蝴蝶曲线;
计算矢量磁场点为(Ba,-Ba)的磁致伸缩Everett函数为第一函数值εa;
利用测量曲线计算其他点(Ba,Bb)点的磁致伸缩Everett函数为第二函数值(εa-εb)/2,Bb为不同于-Ba的其他值;其中,εa、εb分别是B=Ba、B=Bb时的磁致伸缩值,εa、εb是测量对称蝴蝶曲线上的值,蝴蝶曲线纵轴是磁致伸缩值ε,横轴是磁感应强度B;
通过磁致伸缩Everett函数关于直线α=-β的对称性,获得对称点的函数值,并且α为磁致伸缩算子的正向翻转阈值,β为磁致伸缩算子的负向翻转阈值;
分别计算每条测量的蝴蝶回线对应的Everett函数值Ek,并基于所计算的函数值建立磁致伸缩Everett函数数据库。
一种实现方式中,根据Everett函数值Ek计算得出对应的磁致伸缩蝴蝶曲线εk的公式为:
其中,Bt为输入波形极值存储器中的极值序列,θ为最大值方向的方向角,λ为椭圆长短轴的长度比。
一种实现方式中,基于Preisach模型的磁致伸缩模型表示为:
εk(t)=∫∫α≥βμk(θ,λ,α,β)γαβ(Bk(t))dαdβ
其中,γαβ为磁致伸缩算子,并且α为其正向翻转阈值,β为其负向翻转阈值,θ为最大值方向的方向角,λ为椭圆长短轴的长度比;μk(θ,λ,α,β)是输入旋转磁场形状参数拟合成规则椭圆形后形状参数θ和λ的函数,Bk(t)为t时刻k方向上的波形。
应用本发明实施例提供的一种电工钢片磁致伸缩模型的建模方法,具备的有益效果如下:
(1)建立的新型磁致伸缩模型,可以计算不同磁场下电工钢片磁致伸缩特性随时间变化的瞬时值。
(2)通过建立磁致伸缩Everett函数数据库,简化了Preisach模型的计算过程。
附图说明
图1是本发明实施例一种电工钢片磁致伸缩模型的建模方法的一种流程示意图。
图2a是35PH105硅钢片在交变磁场下:当θ=20°时,不同Bmax的实测主磁致伸缩应变εx的蝴蝶曲线。
图2b是35PH105硅钢片在交变磁场下:当Bmax=1.0T时,不同θ的实测主磁致伸缩应变εx的蝴蝶曲线。
图3a是在旋转磁场下,当θ=20°时,不同Bmax的实测主磁致伸缩应变εx的蝴蝶曲线。
图3b是在旋转磁场下,当Bmax=1.0T时,不同θ的实测主磁致伸缩应变εx的蝴蝶曲线。
图4是当Bmax=Ba时的对称磁致伸缩蝴蝶曲线示意图。
图5a为εx和Bx的Everett函数:m=εx,n=Bx。
图5b为εy和By的Everett函数:m=εy,n=By。
图5c为εxy和Bx 2+By 2的Everett函数:m=εxy,n=Bx 2+By 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°时,不同Bmax的实测主磁致伸缩应变εx的蝴蝶曲线;如图2b所示,是当Bmax=1.0T时,不同θ的实测主磁致伸缩应变εx的蝴蝶曲线。在旋转磁场下测量得到:如图3a所示,是当θ=20°时,不同Bmax的实测主磁致伸缩应变εx的蝴蝶曲线;如图3b所示,是当Bmax=1.0T时,不同θ的实测主磁致伸缩应变εx的蝴蝶曲线。
S120,对所测量的对称蝴蝶曲线进行辨识,构建磁致伸缩蝴蝶曲线εk的Everett函数数据库。
可以理解的是,磁致伸缩Everett函数数据库的构建需要依靠测量的对称蝴蝶曲线进行辨识,其辨识步骤如下:
(1)在磁感应强度最大值Bmax为矢量磁场点为Ba的某一给定旋转磁场(θ0,λ0)下测量如图4所示的对称蝴蝶曲线;其中,θ为最大值方向的方向角,λ为椭圆长短轴的长度比。θ0λ0表示给定的确定值。
(2)计算点(Ba,-Ba)的磁致伸缩Everett函数为E(θ0,λ0,Ba,-Ba)=εa;
(3)利用测量曲线计算其他点(Ba,Bb)点的磁致伸缩Everett函数为E(θ0,λ0,Ba,Bb)=(εa-εb)/2;
(4)利用磁致伸缩Everett函数关于直线α=-β的对称性,得到对称点的函数值。
(5)重复上述步骤1~4,分别计算每条测量的蝴蝶回线对应的函数值Ek(θ,λ,B,B),并建立磁致伸缩Everett函数数据库。
S130,将输入的矢量磁场B波形用旋转磁场形状参数Bmax,最大值方向的方向角θ和椭圆长短轴的长度比λ拟合,得到拟合后结果。
需要说明的是,本发明使用的拟合方法为最小二乘近似法:
θ=∠(Bx(τmax),By(τmax))
其中,τmax表示矢量B波形最大时的方向角,Bx和By分别是输入矢量B波形的x方向和y方向的分量,τi表示第i时刻的方向角,Bmax表示B波形的最大值,τ表示某个时刻的方向角,B是输入的矢量磁场B波形,为非椭圆波形;b是椭圆波形,NB是输入的矢量磁场B波形的采样点数。τ=ωt,ω表示角速度,t表示时刻,τ表示某个时刻的方向角。
需要说明的是,将输入的B波形拟合成规则椭圆形后形状参数θ和λ的函数,其中θ为最大值方向的方向角,λ为椭圆长短轴的长度比,如图5a、图5b、图5c是施加在35PH105硅钢片上的旋转磁场为(θ=30°,λ=0.3)时,Everett函数数据库的结果:图5a是εx和Bx的Everett函数,m=εx,n=Bx;如图5b是εy和By的Everett函数,m=εy,n=By;如图5c是εxy和Bx 2+By 2的Everett函数,m=εxy,n=Bx 2+By 2。
S140,通过线性插值法在磁致伸缩Everett函数数据库中进行插值,获得拟合后结果所对应的磁致伸缩Everett函数值Ek。
插值方法如图9所示,已知Everett函数值E(Qij)为点Qij(θi,λj)的Everett函数值,其中,θi参数θ第i个数,λj为参数λ的第j个数;
(1)在θ轴方向上进行插值求得点Qa1(θa,λ1)的Everett函数值:
其中,θR,λ1,Q11、Q21、θ1,θ2、E(Q11)、E(Q21)见图示。
(2)在θ轴方向上进行插值求得点Qa2(θa,λ2)的Everett函数值:
其中,λ2,Q12、Q22、E(Q12)、E(Q22)见图示。
(3)在λ轴方向上进行插值求得点Qab(θa,λb)的Everett函数值:
其中,λR,λ1、λ2、R1,R2、E(R2)、E(R1)见图示。
通过上述线性插值法在磁致伸缩Everett函数数据库中进行插值,求得拟合后结果所对应的磁致伸缩Everett函数值Ek。
S150,根据Everett函数值Ek计算得出对应的磁致伸缩蝴蝶曲线εk。
需要说明的是,本发明基于Preisach模型理论,提出一种可以考虑相应矢量磁场下的主磁致伸缩和剪切磁致伸缩的磁致伸缩模型。基于Preisach模型的磁致伸缩模型表示为:
εk(t)=∫∫α≥βμk(θ,λ,α,β)γαβ(Bk(t))dαdβ (1)
式(1)中,下标k分别表示x,y和xy,则εx和εy分别表示x和y方向的主磁致伸缩,εxy表示剪切磁致伸缩。Bx和By分别是输入B波形的x方向和y方向的分量,并且γαβ为磁致伸缩算子,并且α为其正向+1翻转阈值,β为其负向-1翻转阈值,如图6所示,定义如下:
式(2)中,ti是输入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函数Ek(θ,λ,α,β),用来代替分布函数μk。在给定磁场(θ0,λ0),它所对应的定义式的积分区域,如图8所示,在这个区域里给定一点(α0,β0)((α0,β0)是Preisach平面上给定的一个点)的Everett函数定义如下:
通过式(3)和式(4),可以把Preisach磁致伸缩模型的定义式(1)的复杂计算转换成如下公式:
其中,B1为输入波形极值存储器中序号为1的极值序列、B2为输入波形极值存储器中序号为2的极值序列、Bt为输入波形极值存储器中序号为t的极值序列,N表示t的最大值。极值存储器存储输入矢量B波形的极值过程,如图10所示。其大小排序如下:
|B1|>|B3|>···>|B2i-1|
|B2|>|B4|>···>|B2i|
|B1|>|B2|
经过上述计算,可以求得定义式(1)中Preisach平面上的Everett函数值Ek,从而得出对应的磁致伸缩蝴蝶曲线εk。
图11a为交变磁场Bmax=1.0T,θ=80°,λ=0,图11b为旋转磁场Bmax=1.0T,θ=30°,λ=0.3,图11c为旋转磁场Bmax=1.0T,θ=40°,λ=0.1,图11d为旋转磁场Bmax=1.0T,θ=40°,λ=0.2,通过比较了型号为35PH105硅钢片在不同磁场下基于本发明模型得到的磁致伸缩曲线和测量结果,其中虚线代表模型计算曲线,实线代表测量曲线。结果表明本发明提出的模型结果与试验结果基本吻合,本模型可以预测出与实际情况接近的磁致伸缩特性。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。
Claims (2)
1.一种电工钢片磁致伸缩模型的建模方法,其特征在于,包括:
测量不同旋转磁场下电工钢片磁致伸缩对称蝴蝶曲线εk;
对所测量的对称蝴蝶曲线进行辨识,构建磁致伸缩蝴蝶曲线εk的Everett函数数据库;
将输入的矢量磁场B波形用旋转磁场形状参数Bmax,最大值方向的方向角θ和椭圆长短轴的长度比λ拟合,得到拟合后结果;
通过线性插值法在磁致伸缩Everett函数数据库中进行插值,获得拟合后结果所对应的磁致伸缩Everett函数值Ek;
根据Everett函数值Ek计算得出对应的磁致伸缩蝴蝶曲线εk;
对所测量的对称蝴蝶曲线进行辨识,构建磁致伸缩蝴蝶曲线εk的Everett函数数据库的步骤为:
在磁感应强度最大值Bmax为Ba的一旋转磁场下,测量对称蝴蝶曲线;
计算矢量磁场点为(Ba,-Ba)的磁致伸缩Everett函数为第一函数值εa;
利用测量曲线计算其他点(Ba,Bb)点的磁致伸缩Everett函数为第二函数值(εa-εb)/2,Bb为不同于-Ba的其他值,εa、εb是测量对称蝴蝶曲线上的值;
通过磁致伸缩Everett函数关于直线α=-β的对称性,获得对称点的函数值,并且α为磁致伸缩算子的正向翻转阈值,β为磁致伸缩算子的负向翻转阈值;
分别计算每条测量的蝴蝶回线对应的磁致伸缩Everett函数值Ek,并基于所计算的函数值建立磁致伸缩Everett函数数据库;
根据Everett函数值Ek计算得出对应的磁致伸缩蝴蝶曲线εk的公式为:
其中,Bt为输入波形极值存储器中序号为t的极值序列,θ为最大值方向的方向角,λ为椭圆长短轴的长度比,N表示t的最大值。
2.根据权利要求1所述的一种电工钢片磁致伸缩模型的建模方法,其特征在于,基于Preisach模型的磁致伸缩模型表示为:
εk(t)=∫∫α≥βμk(θ,λ,α,β)γαβ(Bk(t))dαdβ
其中,γαβ为磁致伸缩算子,α为其正向翻转阈值,β为其负向翻转阈值,θ为最大值方向的方向角,λ为椭圆长短轴的长度比;μk(θ,λ,α,β)是输入旋转磁场形状参数拟合成规则椭圆形后形状参数θ和λ的函数,Bk(t)为t时刻k方向上的波形。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111133151.0A CN113761691B (zh) | 2021-09-27 | 2021-09-27 | 一种电工钢片磁致伸缩模型的建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111133151.0A CN113761691B (zh) | 2021-09-27 | 2021-09-27 | 一种电工钢片磁致伸缩模型的建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113761691A CN113761691A (zh) | 2021-12-07 |
CN113761691B true CN113761691B (zh) | 2024-04-02 |
Family
ID=78797715
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111133151.0A Active CN113761691B (zh) | 2021-09-27 | 2021-09-27 | 一种电工钢片磁致伸缩模型的建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113761691B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07110369A (ja) * | 1993-10-08 | 1995-04-25 | Iwatsu Electric Co Ltd | 磁気ひずみ測定装置 |
JP2004184234A (ja) * | 2002-12-03 | 2004-07-02 | Nippon Steel Corp | 磁界解析方法、磁界解析装置、コンピュータプログラム及びコンピュータ読み取り可能な記録媒体 |
CN112083362A (zh) * | 2020-09-10 | 2020-12-15 | 华北电力大学 | 基于可逆磁化分量的电工钢片磁滞特性预测方法及系统 |
-
2021
- 2021-09-27 CN CN202111133151.0A patent/CN113761691B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07110369A (ja) * | 1993-10-08 | 1995-04-25 | Iwatsu Electric Co Ltd | 磁気ひずみ測定装置 |
JP2004184234A (ja) * | 2002-12-03 | 2004-07-02 | Nippon Steel Corp | 磁界解析方法、磁界解析装置、コンピュータプログラム及びコンピュータ読み取り可能な記録媒体 |
CN112083362A (zh) * | 2020-09-10 | 2020-12-15 | 华北电力大学 | 基于可逆磁化分量的电工钢片磁滞特性预测方法及系统 |
Non-Patent Citations (1)
Title |
---|
黄雷 ; 包广清 ; 陈俊全 ; .基于Levenberg-Marquardt算法的改进Preisach模型磁特性模拟与验证.中国电机工程学报.(18),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN113761691A (zh) | 2021-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yin et al. | Thickness measurement of metallic plates with an electromagnetic sensor using phase signature analysis | |
Gmyrek | Single sheet tester with variable dimensions | |
Lin et al. | Distributed current source method for modeling magnetic and eddy-current fields induced in nonferrous metallic objects | |
Weise et al. | Optimal magnet design for Lorentz force eddy-current testing | |
Szlagowska-Spychalska et al. | A novel approach for measuring of thickness of induction hardened layers based on the eddy current method and the finite element modeling | |
Baghel et al. | Effect of mechanical stress on different core loss components along orthogonal directions in electrical steels | |
CN106940343A (zh) | 一种基于阵列电磁传感器的材料微损伤检测方法及系统 | |
Wang et al. | Modeling of magnetostrictive property of electrical steel sheet under vectorial excitation | |
Prauzner | Finite Element Method in an analysis of selected parameters of an inductive sensor for protective coatings measurements | |
CN113761691B (zh) | 一种电工钢片磁致伸缩模型的建模方法 | |
Yoon et al. | Finite element implementation of a generalized Chua-type vector hysteresis model and application to iron loss analysis of three-phase transformer | |
Wang et al. | Study for influence of harmonic magnetic fields on vibration properties of core of anode saturable reactor in HVDC converter valve system | |
CN114580235A (zh) | 一种换流变压器振动缩比模型的设计方法 | |
Zhu et al. | A Preisach-based magnetostriction model for highly grain-oriented electrical steel under rotating magnetic field | |
Pluta | Calculating power loss in electrical steel taking into account magnetic anisotropy | |
JP4044891B2 (ja) | 磁界解析方法及びコンピュータプログラム | |
CN114613453A (zh) | 基于ja模型的动态磁致伸缩确定方法及系统 | |
Wang et al. | Magnetostrictive characteristics of the grain-oriented electrical steel in an Epstein frame magnetized with a DC biased magnetic field | |
Miyagi et al. | Usefulness of fixed point method in electromagnetic field analysis in consideration of nonlinear magnetic anisotropy | |
Warzecha et al. | Magnetization measurements in circle-shaped samples of typical dynamo steel sheets | |
Vandevelde et al. | Magnetostriction and magnetic forces in electrical steel: Finite element computations and measurements | |
Zhang et al. | A new hysteresis loss estimation in the induction motor core considering rotating magnetic fields | |
Wakabayashi et al. | Analysis of vector magnetic property and two-dimensional magnetostriction in electrical steel sheets | |
Clénet et al. | Determination of losses’ local distribution for transformer optimal designing | |
Sugahara | Electromagnetic Analysis of Eddy Current Testing With Kelvin Transformation |
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 |