CN109086561A - 一种计及各向异性的ja磁滞模型参数提取方法 - Google Patents
一种计及各向异性的ja磁滞模型参数提取方法 Download PDFInfo
- Publication number
- CN109086561A CN109086561A CN201811168788.1A CN201811168788A CN109086561A CN 109086561 A CN109086561 A CN 109086561A CN 201811168788 A CN201811168788 A CN 201811168788A CN 109086561 A CN109086561 A CN 109086561A
- Authority
- CN
- China
- Prior art keywords
- frog
- hysteresis
- model parameter
- magnetic
- anisotropic
- 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/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明公开了一种计及各向异性的JA磁滞模型参数提取方法,该方法包括以下步骤:S1.建立考虑各向异性的无磁滞磁化曲线模型;S2.求解JA磁滞电感;S3.采用差分‑混合蛙跳复合算法获取JA磁滞模型参数。本发明提出的计及各向异性的JA磁滞模型参数提取方法,从原理上解决了JA磁滞模型参数传统提取方法所存在的技术经济难题,且计算速度快、精度高,具有十分重要的应用价值。
Description
技术领域
本发明涉及电力技术领域,特别是一种计及各向异性的JA(Jiles-Atherton)磁滞模型参数提取方法。
背景技术
变压器饱和特性和铁芯损耗在低频电磁暂态过程中具有十分重要的作用。在磁化过程中存在磁感应强度变化滞后于磁场强度变化的现象,即磁滞效应。现有研究通常采用单值非线性电感和电阻的并联组合表征磁滞特性,不能准确模拟光滑的磁滞回线;现有JA磁滞模型采用朗之万函数描述无磁滞磁化曲线,其结构简单、计算方便,但在模拟各向异性特性方面还存在明显不足。此外,对经典JA磁滞模型进行修正后,模型复杂程度显著增加,而在JA模型参数的优化算法方面,现有研究表明GA和SA在对经典JA磁滞模型参数辨识时能够取得一定的求解精度但在求解速度和收敛性能方面还存在不足。因此,开展计及各向异性的JA磁滞模型参数提取方法研究对于保障变压器电磁暂态模型的应用具有极其重要的理论意义。
发明内容
有鉴于此,为了解决上述问题,本发明提供一种计及各向异性的JA磁滞模型参数提取方法,为变压器电磁暂态模型研究提供了基础。
为实现上述目的及其他目的,本发明提供一种计及各向异性的JA磁滞模型参数提取方法,该方法包括以下步骤:
S1.建立考虑各向异性的无磁滞磁化曲线模型;
其中,Ms为饱和磁化强度,Em(θ)表示磁矩能量,kB表示Boltzman常数,θ表示原子磁矩mat和磁场方向之间的夹角;
S2.求解JA磁滞电感;
S3.采用差分-混合蛙跳复合算法获取JA磁滞模型参数。
可选地,所述求解JA磁滞电感,具体包括以下步骤:
①当t=(t0+Δt)时,该时刻电感两端电压u(t0+Δt)已知,由式dB/dt=u(t0+Δt)/(N·S)对电压进行积分变换,依据B(t0+Δt)=B(t0)+ΔB/Δt·Δt≈B(t0)+dB/dt·Δt得到磁感应强度B(t0+Δt);
N为绕组匝数,S为铁芯横截面积,Δt为时间增量,ΔB为磁感应强度增量;
②将磁感应强度B(t0),磁场强度H(t0),磁化强度M(t0),dB/dt,以及dH/dt作为Matlab Function模块输入,计算得到dM/dt;
③利用表达式M(t0+Δt)=M(t0)+dM/dt·Δt,计算(t0+Δt)时刻磁化强度M(t0+Δt),由式H(t0+Δt)=B(t0+Δt)/μ0–M(t0+Δt)可以计算得到磁场强度H(t0+Δt),根据安培环路定律N·i=H·l,可以输出此时电流i(t0+Δt)。
可选地,所述采用差分-混合蛙跳复合算法获取JA磁滞模型参数,具体包括以下步骤:
S31.获取试验磁场强度数据HT;
S32.随机生成参数饱和磁化强度Ms,无磁滞磁化曲线形状参数a,平均场参数α,磁场方向和晶体特征轴方向夹角ψ,牵制系数k,磁化因数c,单轴各向异性的平均能量密度Kan;
S33.求解JA微分方程;
S34.根据试验磁感应强度BT,计算每一青蛙适应度值E,并据此排序;
S35.划分mp个青蛙子组族群,并在每一个子族群内进行局部深度搜索;
S36.混合整个族群并重新排序;
S37.生成新一代青蛙种群;
S38.判断是否达到最大迭代次数,如果是则输出最优青蛙个数;如果否,则重复步骤S35~步骤S37。
可选地,所述在每一个子族群内进行局部深度搜索,具体包括以下步骤;
S351.记录本子族群的局部最优青蛙Ub和局部最差青蛙Uw以及整个种群的全局最优青蛙Ug;
S352.进行最差青蛙第一次更新;
S353.判断更新前后青蛙性能是否提高;若是,则取代本子族群局部最差青蛙Uw;若否,则用全局最优青蛙Ug替代局部最优青蛙Ub进行第二次更新;
S354.再次比较更新前后青蛙性能是否提高,如果是,则取代本子族群局部最差青蛙Uw,如果否,则随机产生一个青蛙取代本子族群局部最差青蛙Uw;
S355.判断是否达到局部搜索次数,若是则局部搜索结速,若否则重复步骤S351~S354。
由于采用了上述技术方案,本发明具有如下的优点:
本发明提出的计及各向异性的JA磁滞模型参数提取方法,从原理上解决了JA磁滞模型参数传统提取方法所存在的技术经济难题,且计算速度快、精度高,具有十分重要的应用价值。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步的详细描述:
图1为本发明所述方法的流程图;
图2为JA磁滞电感在MATLAB中求解;
图3为JA参数DE-SFLA优化流程图;
图4为取向和无取向硅钢磁滞回线测量的原理图;
图5为GO Bmax=1.0T的试验与优化后仿真磁滞回线;
图6为GO Bmax=1.5T的试验与优化后仿真磁滞回线;
图7为NGO Bmax=1.0T的试验与优化后仿真磁滞回线;
图8为NGO Bmax=1.5T的试验与优化后仿真磁滞回线。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。
需要说明的是,本实施例中所提供的图示仅以示意方式说明本发明的基本构想,遂图式中仅显示与本发明中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制,其实际实施时各组件的型态、数量及比例可为一种随意的改变,且其组件布局型态也可能更为复杂。
如图1所示,本发明提供一种计及各向异性的JA磁滞模型参数提取方法,包括:
S1:建立考虑各向异性的无磁滞磁化曲线模型。
S2:求解JA磁滞电感。
S3:采用差分-混合蛙跳复合算法获取本文JA磁滞模型参数,并开展取向硅钢和无取向硅钢磁滞回线测量,分别比较不同算法在优化速度、精度以及收敛速度等方面的性能差异。
下面以一具体实施方式对本发明的计及各向异性的JA磁滞模型参数提取方法进行说明。具体步骤为:
S1:建立考虑各向异性的无磁滞磁化曲线模型。
当施加外磁场后,磁性材料磁化强度可能从易磁化轴方向转出,表现出磁晶各向异性,晶体磁晶各向异性能Ek可以用自发磁化与相互正交的晶体学主轴的方向余弦函数表示:
EK=K1(α1 2α2 2+α2 2α3 2+α3 2α1 2)+K2α1 2α2 2α3 2 (1)
式中:K1,K2为立方晶体磁晶各向异性常数,α1,α2和α3为方向余弦。
根据式(1)建立一种考虑各向异性的无磁滞磁化曲线模型,对于各向异性材料,由于各向异性能在不同方向不一致,此时无磁滞磁化曲线只能用磁矩在外磁场方向的总和进行计算,可采用Boltzman分布表示为:
式中:Ms为饱和磁化强度,θ表示原子磁矩mat和磁场方向之间的夹角,kB为Boltzman常数。磁矩能量Em(θ)可以表示为:
Em(θ)=-μ0·mat·H·cosθ (3)
式中:H表示磁场强度;μ0为真空磁导率。
对于各向异性材料,磁矩与磁场方向夹角不尽相同。在二维平面上,设磁矩m1,m2与磁场夹角为θ,简化后磁晶各向异性能为:
式中:K0和K1为各向异性常数,可通过试验测量得到,ψ为磁场方向和晶体特征轴方向夹角。
S2:求解JA磁滞电感。
比较各向异性函数与传统朗之万函数的性能差异,分析ode45,ode23以及四阶定步长Runge-Kutta算法对复合JA磁滞电感求解时间的区别;
在Matlab中求解上述JA磁滞电感的过程。如图2所示,具体求解步骤如下:
①当t=(t0+Δt)时,该时刻电感两端电压u(t0+Δt)已知,由式dB/dt=u(t0+Δt)/(N·S)对电压进行积分变换,依据B(t0+Δt)=B(t0)+ΔB/Δt·Δt≈B(t0)+dB/dt·Δt得到磁感应强度B(t0+Δt);
N为绕组匝数,S为铁芯横截面积,Δt为时间增量,ΔB为磁感应强度增量;
②将磁感应强度B(t0),磁场强度H(t0),磁化强度M(t0),dB/dt,以及dH/dt作为Matlab Function模块输入,计算得到dM/dt;
③利用表达式M(t0+Δt)=M(t0)+dM/dt·Δt,计算(t0+Δt)时刻磁化强度M(t0+Δt),由式H(t0+Δt)=B(t0+Δt)/μ0–M(t0+Δt)可以计算得到磁场强度H(t0+Δt),根据安培环路定律N·i=H·l,可以输出此时电流i(t0+Δt)。
以典型的单轴晶体各向异性材料模型参数为例,即k1=0.011,k2=0.049,a=1000A/m,k=5000A/m,c=0.1,Ms=1.3×106A/m,α=0.001,Kan=4×104J/m3,ψ=45o为例,分别采用ode45,ode23,以及四阶定步长Runge-Kutta法通过Matlab编程求解JA磁滞电感,,计算机操作***为:Windows 7Enterprise,***类型:64位,CPU:Intel Core i5–46703.40GHz。不同求解算法下单次JA磁滞模型求解所需时间如表1所示,四阶定步长Runge-Kutta算法求解效率最高,其求解时间仅为ode23算法的26.8%,ode45求解时间介于四阶定步长Runge-Kutta算法和ode23算法之间。因此,在对JA磁滞模型进行求解时采用四阶定步长Runge-Kutta算法最合适。
表1求解算法对本文JA磁滞模型计算时间影响
S3:采用差分-混合蛙跳复合算法获取本文JA磁滞模型参数,并开展取向硅钢和无取向硅钢磁滞回线测量,分别比较不同算法在优化速度、精度以及收敛速度等方面的性能差异:
引入DE/target-to-best/1型差分变异算子后的青蛙更新步长为:
S1=rand·(Ub-Uw)+rand·(Ur1-Ur2) (5)
式中:Ur1和Ur2为本族群中与最差和最优青蛙均不同的两随机个体。
提出如下自适应步长系数:
式中:M1和M2分别为当前迭代次数和最大迭代次数。
引入自适应步长系数后的青蛙更新步长变为:
S2=k1·(Ub-Uw)+k2·(Ur1-Ur2) (7)
引入惯性权值α,可提高原始算法的局部搜索和全局搜索能力,从而得到新的更新公式:
Uw new=α·Uw+S2 (8)
式中:α为惯性权值,其取值范围为[0,1]。
依据上述原理,对经典SFLA中青蛙更新策略分别引入差分变异思想、自适应步长系数以及惯性权值后即可得到DE-SFLA的优化原理。采用DE-SFLA对JA模型参数优化的流程如图3所示。
S31.获取试验磁场强度数据HT;
S32.随机生成参数饱和磁化强度Ms,无磁滞磁化曲线形状参数a,平均场参数α,磁场方向和晶体特征轴方向夹角ψ,牵制系数k,磁化因数c,单轴各向异性的平均能量密度Kan;
S33.求解JA微分方程;
式中:δ表示方向系数,Man表示无磁滞磁化强度。
S34.根据试验磁感应强度BT,计算每一青蛙适应度值E,并据此排序;
式中:N为试验数据长度,BTj为试验磁感应强度,BSj为计算所得磁感应强度,BTmax为试验磁感应强度最大值。
S35.划分mp个青蛙子组族群,并在每一个子族群内进行局部深度搜索;
S36.混合整个族群并重新排序;
S37.生成新一代青蛙种群;
S38.判断是否达到最大迭代次数,如果是则输出最优青蛙个数;如果否,则重复步骤S35~步骤S37。
可选地,所述在每一个子族群内进行局部深度搜索,具体包括以下步骤;
S351.记录本子族群的局部最优青蛙Ub和局部最差青蛙Uw以及整个种群的全局最优青蛙Ug;
S352.按式(11)进行最差青蛙第一次更新;
α=(α2-α1)·(N2-N1)/N2+α1 (11)
式中:α1和α2分别为惯性权值最小值和最大值;N1和N2分别为当前局部搜索迭代次数和局部搜索最大迭代次数。
S353.判断更新前后青蛙性能是否提高;若是,则取代本子族群局部最差青蛙Uw;若否,则用全局最优青蛙Ug替代式(7)局部最优青蛙Ub进行第二次更新;
S354.再次比较更新前后青蛙性能是否提高,如果是,则取代本子族群局部最差青蛙Uw,如果否,则随机产生一个青蛙取代本子族群局部最差青蛙Uw;
S355.判断是否达到局部搜索次数,若是则局部搜索结速,若否则重复步骤S351~S354。
本实施例采用湖南联众科技有限公司的MATS–2010M硅钢测量装置对取向(GO)和无取向(NGO)硅钢环形铁芯(20VA,10V)进行测量,试验原理如图4所示。1为电流探头,2为绕组,3为样品,4为计算机。相同尺寸的GO和NGO铁芯几何参数如表2所示。采用MATS–2010M硅钢测量装置固定最大磁感应强度Bm模式,即设定Bm值,励磁电源自动加压直到B=Bm,在频率为50Hz下,分别对Bm=1.0T和Bm=1.5T的磁滞回线进行测量。
表2铁芯几何尺寸
分别采用上节所述优化算法,在Matlab 2013a中依据实测B-H磁滞回线对JA模型参数进行辨识(计算机操作***:64位Windows 7 Enterprise,CPU:Intel Core i5–46703.40GHz)。
优化后的磁滞回线结果如图5~8所示,从图中可以看出,DE、GA和DE-SFLA三种算法在GO和NGO两种材料不同磁感应强度下均能达到全局最优解附近,而SA在磁感应强度未饱和时(Bmax=1.0T),其结果能够收敛于全局最优解,而当磁感应强度饱和时(Bmax=1.5T),其优化结果为局部最优。在对SA的温度函数以及参数邻域矩阵进行反复调整后,其结果依旧陷入局部最优。因此,SA在磁滞回线饱和情况下对JA模型参数进行优化时容易得到局部最优解。
本发明提出的计及各向异性的JA磁滞模型参数提取方法,从原理上解决了JA磁滞模型参数传统提取方法所存在的技术经济难题,且计算速度快、精度高,具有十分重要的应用价值。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本技术方案的宗旨和范围,其均应涵盖在本发明的保护范围当中。
Claims (4)
1.一种计及各向异性的JA磁滞模型参数提取方法,其特征在于,该方法包括以下步骤:
S1.建立考虑各向异性的无磁滞磁化曲线模型;
其中,Ms为饱和磁化强度,Em(θ)表示磁矩能量,kB表示Boltzman常数,θ表示原子磁矩mat和磁场方向之间的夹角;
S2.求解JA磁滞电感;
S3.采用差分-混合蛙跳复合算法获取JA磁滞模型参数。
2.根据权利要求1所述的一种计及各向异性的Jiles-Atherton磁滞模型参数提取方法,其特征在于,所述求解JA磁滞电感,具体包括以下步骤:
①当t=(t0+Δt)时,该时刻电感两端电压u(t0+Δt)已知,由式dB/dt=u(t0+Δt)/(N·S)对电压进行积分变换,依据B(t0+Δt)=B(t0)+ΔB/Δt·Δt≈B(t0)+dB/dt·Δt得到磁感应强度B(t0+Δt);
N为绕组匝数,S为铁芯横截面积,Δt为时间增量,ΔB为磁感应强度增量;
②将磁感应强度B(t0),磁场强度H(t0),磁化强度M(t0),dB/dt,以及dH/dt作为MatlabFunction模块输入,计算得到dM/dt;
③利用表达式M(t0+Δt)=M(t0)+dM/dt·Δt,计算(t0+Δt)时刻磁化强度M(t0+Δt),由式H(t0+Δt)=B(t0+Δt)/μ0–M(t0+Δt)可以计算得到磁场强度H(t0+Δt),根据安培环路定律N·i=H·l,可以输出此时电流i(t0+Δt)。
3.根据权利要求1所述的一种计及各向异性的Jiles-Atherton磁滞模型参数提取方法,其特征在于,所述采用差分-混合蛙跳复合算法获取JA磁滞模型参数,具体包括以下步骤:
S31.获取试验磁场强度数据HT;
S32.随机生成参数饱和磁化强度Ms,无磁滞磁化曲线形状参数a,平均场参数α,磁场方向和晶体特征轴方向夹角ψ,牵制系数k,磁化因数c,单轴各向异性的平均能量密度Kan;
S33.求解JA微分方程;
S34.根据试验磁感应强度BT,计算每一青蛙适应度值E,并据此排序;
S35.划分mp个青蛙子组族群,并在每一个子族群内进行局部深度搜索;
S36.混合整个族群并重新排序;
S37.生成新一代青蛙种群;
S38.判断是否达到最大迭代次数,如果是则输出最优青蛙个数;如果否,则重复步骤S35~步骤S37。
4.根据权利要求3所述的一种计及各向异性的Jiles-Atherton磁滞模型参数提取方法,其特征在于,所述在每一个子族群内进行局部深度搜索,具体包括以下步骤;
S351.记录本子族群的局部最优青蛙Ub和局部最差青蛙Uw以及整个种群的全局最优青蛙Ug;
S352.进行最差青蛙第一次更新;
S353.判断更新前后青蛙性能是否提高;若是,则取代本子族群局部最差青蛙Uw;若否,则用全局最优青蛙Ug替代局部最优青蛙Ub进行第二次更新;
S354.再次比较更新前后青蛙性能是否提高,如果是,则取代本子族群局部最差青蛙Uw,如果否,则随机产生一个青蛙取代本子族群局部最差青蛙Uw;
S355.判断是否达到局部搜索次数,若是则局部搜索结速,若否则重复步骤S351~S354。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811168788.1A CN109086561B (zh) | 2018-10-08 | 2018-10-08 | 一种计及各向异性的ja磁滞模型参数提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811168788.1A CN109086561B (zh) | 2018-10-08 | 2018-10-08 | 一种计及各向异性的ja磁滞模型参数提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109086561A true CN109086561A (zh) | 2018-12-25 |
CN109086561B CN109086561B (zh) | 2019-09-24 |
Family
ID=64843243
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811168788.1A Active CN109086561B (zh) | 2018-10-08 | 2018-10-08 | 一种计及各向异性的ja磁滞模型参数提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109086561B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110008645A (zh) * | 2019-04-30 | 2019-07-12 | 广东工业大学 | 一种变压器损耗计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103279625A (zh) * | 2013-06-15 | 2013-09-04 | 国家电网公司 | 不同结构变压器铁心等值电路模型建立方法 |
CN107515992A (zh) * | 2017-09-06 | 2017-12-26 | 武汉大学 | 一种铁磁材料磁滞特性仿真中的ja模型参数优化自适应遗传算法 |
CN107609298A (zh) * | 2017-09-26 | 2018-01-19 | 广东电网有限责任公司电力科学研究院 | 一种Jiles‑Atherton模型参数辨识方法和装置 |
CN108519569A (zh) * | 2018-05-07 | 2018-09-11 | 上海交通大学 | 含应力的ja磁滞模型参数识别方法 |
-
2018
- 2018-10-08 CN CN201811168788.1A patent/CN109086561B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103279625A (zh) * | 2013-06-15 | 2013-09-04 | 国家电网公司 | 不同结构变压器铁心等值电路模型建立方法 |
CN107515992A (zh) * | 2017-09-06 | 2017-12-26 | 武汉大学 | 一种铁磁材料磁滞特性仿真中的ja模型参数优化自适应遗传算法 |
CN107609298A (zh) * | 2017-09-26 | 2018-01-19 | 广东电网有限责任公司电力科学研究院 | 一种Jiles‑Atherton模型参数辨识方法和装置 |
CN108519569A (zh) * | 2018-05-07 | 2018-09-11 | 上海交通大学 | 含应力的ja磁滞模型参数识别方法 |
Non-Patent Citations (3)
Title |
---|
DIANHAI ZHANG: ""Analysis of Magnetizing Process of a New Anisotropic Bonded NdFeB Permanent Magnet Using FEM Combined With Jiles-Atherton Hysteresis Model"", 《 IEEE TRANSACTIONS ON MAGNETICS》 * |
王洋 等: ""基于蛙跳模糊算法的Jiles Atherton 铁心磁滞模型参数确定"", 《电工技术学报》 * |
郭子政: ""Jiles-Atherton模型及其应用"", 《信息记录材料》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110008645A (zh) * | 2019-04-30 | 2019-07-12 | 广东工业大学 | 一种变压器损耗计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109086561B (zh) | 2019-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Caux et al. | Computation of dynamical correlation functions of Heisenberg chains in a magnetic field | |
Grinstein | Critical reanalysis of CP asymmetries in B0 decays to CP eigenstates | |
Shannon et al. | Nematic, vector-multipole, and plateau-liquid states in the classical O (3) pyrochlore antiferromagnet with biquadratic interactions in applied magnetic field | |
Tocchio et al. | Spin-liquid and magnetic phases in the anisotropic triangular lattice: The case of κ-(ET) 2 X | |
CN105302975B (zh) | 一种电磁式电流互感器谐波传变建模方法 | |
CN110399677B (zh) | 基于偏磁状态下改进j-a公式的变压器直流偏磁仿真模拟方法 | |
CN109086561B (zh) | 一种计及各向异性的ja磁滞模型参数提取方法 | |
CN104991994A (zh) | 一种基于j-a磁滞模型的变压器励磁电流仿真方法 | |
CN209232531U (zh) | 一种便携式核磁共振永磁体 | |
Silveyra et al. | On the modelling of the anhysteretic magnetization of homogeneous soft magnetic materials | |
Duan et al. | Comparison of limiting loop model and elemental operator model for magnetic hysteresis of ferromagnetic materials | |
Dodrill et al. | High-temperature first-order-reversal-curve (FORC) study of magnetic nanoparticle based nanocomposite materials | |
Meesuk et al. | Magnetic field analysis for a distribution transformer with unbalanced load conditions by using 3-d finite element method | |
Batista et al. | Electron spin resonance of defects in the Haldane system Y 2 BaNiO 5 | |
CN104899441A (zh) | 一种适用于计算磁场vrms均匀度的数值算法 | |
Pang et al. | Analysis of current and magnetic field distributions in rail launcher with Peaceman–Rachford finite-difference method | |
Najgebauer | Fractional scaling of magnetic coercivity in electrical steels | |
CN105528506B (zh) | 一种饱和电抗器等值电路模型的参数获取方法 | |
Richter et al. | Magnetic susceptibility of frustrated spin-s J1-J2 quantum Heisenberg magnets: High-temperature expansion and exact diagonalization data | |
Li et al. | An Improved Jiles-Atherton Model for DC Bias of Transformer | |
Wang et al. | 3-D FEM analysis in electromagnetic system considering vector hysteresis and anisotropy | |
Zheng et al. | Mixed-spin S=(1 2, 1) XXZ ferrimagnet at zero temperature: Cluster series expansions | |
Zhang et al. | Monte Carlo simulation of the Ising ferromagnetic iron nitride system | |
Wang et al. | Effect of ball scribing on relative permeability of grain-oriented electrical steel | |
Coffey | Analyzing mutual inductance measurements to determine the London penetration depth |
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 |