基于力磁耦合效应的铁磁材料修正磁化模型*

罗旭 朱海燕 丁雅萍

1) (西南石油大学机电工程学院,成都 610500)

2) (成都理工大学能源学院,成都 610059)

3) (四川建筑职业技术学院交通与市政工程系,成都 610399)

1 引 言

应力-应变对铁磁材料磁化过程的影响规律,即力磁耦合效应,一直是国内外学者的研究热点[1−9].随着应力-应变的改变,特别是应力集中或塑性变形的出现,铁磁材料内部的磁畴结构和钉扎点密度等会发生显著变化,从而影响材料的磁化曲线、矫顽力、剩余磁化强度等磁化特性[10−14].基于上述现象,发展出多种磁无损检测技术,如: 金属磁记忆检测技术、巴克豪森磁噪声技术、磁声发射技术.建立准确的磁化模型来描述应力、塑性变形对铁磁材料磁化的影响规律,是利用上述磁无损检测技术实现应力集中、塑性变形检测与反演的基础.

1984年,Jiles和Atherton[15,16]根据接近定理,将磁化强度分为可逆磁化和不可逆磁化两部分,建立了著名的J-A模型; 在该模型的基础上,Sablik等[6]将应力等效为附加磁场,建立了反映均匀弹性应力对磁化过程影响的Jiles-Atherton-Sablik (J-A-S)模型; 初期,Sablik等[6,17]认为磁致伸缩系数只与磁化强度的二次项有关,并给出了三种不同形式的磁致伸缩应变关系式; 此后,Sablik等[7,18]基于磁畴理论和能内定理,建立了新的磁致伸缩应变关系式; 由于实际材料中磁畴结构很难提前预知,Jiles[8]提出利用经验公式来拟合磁致伸缩应变与磁化强度之间的变化关系,并将其Taylor展开式简化为只包含磁化强度二次项和四次项的式子,由此形成了J-A模型的中经典磁致伸缩应变关系式; 此后,文献[19−21]针对原磁致伸缩系数关系式不能很好反映低应力条件下磁致伸缩系数的变化规律,提出利用神经网络来拟合磁致伸缩系数与磁化强度关系.为描述拉应力和压应力对磁化特性影响的非对称性,Li和Xu[22]以及 Sablik[23]在原有模型基础上引入退磁项,并基于Gaussian方程建立了饱和磁致伸缩系数与应力的变化关系; Li和Jiles[24]则利用Rayleigh定理对J-A模型中的接近定理进行了修正; 针对循环应力作用下零应力附近的磁化强度不符合逐渐趋于稳定状态的问题,Sablik等[25]提出了两种方法对接近定理进行改进.为实现利用J-A模型描述塑性变形对材料磁化的影响,文献[3,26−33]分析了塑性变形对材料位错密度的影响,建立了钉扎系数、Weiss分子场耦合系数、形状系数等与晶粒尺寸、位错密度之间的关系,推导了基于J-A模型的不同形式的塑性磁化模型.

J-A模型中,应力对磁化的影响是通过应力改变磁致伸缩系数实现的.因此,磁致伸缩应变关系式是决定J-A模型准确性的关键.目前,不同形式J-A模型中的磁致伸缩应变关系式都是采用Jiles提出的经典拟合公式,其中磁化强度二次项和四次项的系数与应力均为线性关系,并不能很好地描述铁磁材料磁致伸缩应变随应力、磁化强度的非线性变化规律.针对这一问题,Zhou等[34,35]通过在Zheng Xiao-Jing-Liu Xing-En (Z-L)模型中的磁致伸缩应变关系中引入磁化强度四次项,更好地反映了铁磁材料磁致伸缩应变随应力、磁化强度的非线性变化规律; 但是,周浩淼等改进的磁化模型中没有考虑Weiss分子场、退磁场和磁滞效应的影响,不能反映循环磁化过程中的磁滞现象.Shi等[36,37]在考虑退磁场和磁滞效应的基础上,对Z-L模型进行了改进,但仍未考虑Weiss分子场作用,同时由于采用了基于弹性能的接近定理,无法描述塑性变形对磁化曲线的影响.

针对上述问题,本文在结合Z-L模型中的非线性磁致伸缩应变关系式以及J-A模型中的磁滞理论基础上,考虑应力、塑性变形对模型参数的影响,建立了能够反映弹-塑性应力与应变对铁磁材料磁化曲线影响的修正磁化模型.通过与试验结果及原有模型计算结果进行对比,验证改进模型的有效性、准确性; 同时,讨论了弹性拉、压应力和塑性拉、压应变对磁化曲线、矫顽力和剩余磁化强度的影响,可为分析应力、应变对铁磁材料磁化的影响规律提供新选择.

2 非线性力磁耦合模型

2.1 非磁滞条件下的力磁耦合模型

当铁磁材料所受应力方向与外加磁场方向同轴时,根据热力学第一定律,材料单位体积内能U(σ,Man,T)满足[35]

式中,Q为材料从外部吸收的热能;W为外力对材料做所的功;T为外界温度;S为体积熵的变化量;σ为应力;ε为应变;µ0为真空磁导率,µ0=4π×10−7;H为外磁场强度;Man为非磁滞磁化强度.

根据可逆体系的热力学平衡原理,理想磁化过程中铁磁材料的Gibbs自由能Gm(σ,Man,T) 可表示为[36]

式中,α为Weiss分子场耦合系数,Nd为退磁系数.则Gm(σ,Man,T) 的全微分可写为

当不考虑温度变化时,即 dT=0 ,则可得到只包含应力σ与非磁滞磁化强度Man的Gibbs自由能Gm(σ,Man) 全微分表达式,即

由上述热力学关系可得到非磁滞条件下铁磁材料应变ε(σ,Man) 和有效场He(σ,Man) 关系为[34−36]

为获得多项式形式的本构关系,将Gm(σ,Man)在 (σ,Man)=(0,0) 进行泰勒展开,并根据试验现象对泰勒展开式进行简化.根据Kurzar和Cullity[38],Yamasaki等[39]的试验结果及Jiles、周浩淼等的理论分析,铁磁材料的磁致伸缩曲线是关于y轴对称的偶函数,因此磁致伸缩应变ε(σ,Man)中的耦合项只保留含有项,则

式中,Es为材料固有弹性模量;λ0(σ) 为应力单独作用下磁畴移动导致的磁致伸缩应变量;λmax(σ)为应力单独作用时由于磁畴移动导致的最大磁致伸缩应变量;Mws为无应力时的饱和壁移磁化强度;λws为无应力条件下材料的饱和磁致伸缩应变;M0(σ)为应力作用时的饱和壁移磁化强度;ϑ为阶跃函数,当Man

其中,Ms为饱和磁化强度;a为形状系数; tanh(x)是双曲正切函数;σ0表示将λ0(σ) 在原点附近线性化为斜率为 1/E′的直线时,使λ0(σ)=λws时的应力值,1/E′=1/E0−1/Es,E0为初始弹性模量;β为比例系数,影响曲线斜率.周浩淼改进的Z-L模型采用了分段函数来表示压应力和拉应力对λ0(σ)的不同影响,为简化公式,本文通过平移变换,将Z-L模型中λ0(σ) 的表达式变为了连续函数,其中−arctanh(1/3)的引入是为保证当应力σ=0 MPa时,λ0(σ)=0.

将(8)式代入(7)式中,即可得到非磁滞条件下材料磁致伸缩应变和磁化强度随应力、磁场的变化规律.

2.2 磁滞效应对磁化强度的影响

实际磁化过程中,材料内部不可避免地存在缺陷、位错、晶格滑移等钉扎点,从而阻碍磁畴运动.钉扎效应的存在,导致材料在磁化过程中不可避免地存在能量损耗.因此,根据能量守恒原理,实际磁化过程中的磁化能量等于无磁滞磁化能减去钉扎点造成的磁滞损耗[33],即

式中,M为实际磁化强度,M=Mrev+Mirr,Mirr为不可逆磁化强度,Mrev为可逆磁化强度;keff为有效钉扎系数;δ为方向系数,当dHe(M)/dt>0时,δ=1 ,当 dHe(M)/dt<0 时,δ=−1 ,以保证钉扎点的作用总是阻碍磁场的变化;Be(M) 表示实际磁感应强度,Be(M)=µ0He(M) ;He(M) 为实际磁化强度条件下的等效磁场强度,He(M) 可将(7)式中的Man换成M得到,即

整理(9)式可得

由M=Mrev+Mirr ,Mrev=c(Man−Mirr) 可得

(12)式对外磁场强度H取微分可得

(10)式对外磁场强度H取微分可得

联立公(11),(13)和(14)式,并整理可得

其中,dMan/dH可联立(7)和(8)式得到,即

联立(15)和(16)式即可计算不同应力条件下磁化强度M随外磁场H的变化规律.

2.3 应力-应变对模型参数的影响

在应力、应变,特别是塑性变形作用下,材料内部磁畴结构的形状和尺寸以及钉扎点密度都会发生改变,从而改变Weiss分子场耦合系数α、有效钉扎系数keff以及形状系数a的取值.

在一定范围内,应变硬化应力σF与塑性应变εr之间的关系可根据Ludwik定律得到[40]:

式中A为常数;m为Ludwik指数,与材料性质有关.

根据Astie-Degauque关系,位错密度ζd与材料的应变硬化应力σF呈比例关系[28]:

式中,ζd0为无应力时材料初始位错密度;G为剪切模量,其与固有弹性模量Es和泊松比υ之间的关系为G=Es/2(1+υ) ;b为柏氏矢量,是描述位错实质的物理量;αk为常数,取0.76.

Lo等[28,29]的试验结果表明: 拉应力使钉扎密度减小,压应力使钉扎密度增大; 而位错会以钉扎的形式阻碍畴壁运动.因此,有效钉扎系数keff和形状系数a与应力和位错密度之间的关系可表示为[30]

由于应力-应变改变了磁畴的形状和尺寸,Weiss分子场耦合系数α会随着应力-应变的增加而发生变化,其与应力-应变的关系为[23,31,41]

式中,α0为初始Weiss分子场耦合系数;q1ln(q2εr+q3)表示塑性变形对α改变量;

表示弹性应力对α改变量;q1,q2,q3均为拟合常数.

3 模型理论计算结果与试验结果对比

3.1 单次磁化过程中应力对磁化曲线的影响

为说明修正模型相对于Z-L模型、J-A模型的改进情况,利用Jiles等[8]进行的–200,–100,0,100和200 MPa五组不同预应力条件下的磁化测试结果进行验证,试验结果如文献[8]图3虚线所示.Jiles模型中参数与文献[8]中的图2一致; Z-L模型中相关参数取值与文献[35]中图9一致; 修正模型中β=1.4 ,ϑ=0.75 ,α0=1.44×10−3,其他参数与Z-L模型参数取值一致,不同应力条件下不同模型计算结果与试验结果对比如图1所示.

由试验结果可知: 与无应力条件下的磁化曲线相比,应力将阻碍材料的磁化,且在低磁场阶段压应力对磁化的阻碍作用明显高于拉应力; 当外磁场增强至中高磁场阶段后,压应力条件下的磁化曲线和拉应力条件下的磁化曲线将发生反转,最终压应力条件下的磁化曲线将趋于自由状态下的磁化曲线,而拉应力条件下的磁化曲线会比无应力条件下磁化小很多.对比三种模型计算结果可知: J-A模型能够定性反映应力对磁化的影响,但是定量上却存在一定差别,例如试验结果的无应力饱和磁感应强度在1.6 T左右,而J-A模型的无应力饱和磁感应强度达到了1.8 T左右; 相对于J-A模型,Z-L模型在定性和定量上都有很大的改进,能够很好地反映应力对磁化影响,但是定量上仍然存在一定差异; 而本文修正后的模型无论在定性上还是定量上都能十分准确地反映试验结果,不同应力条件下的理论预测结果与试验结果基本重合.

图1 单次磁化条件下不同模型计算得到的 Be - H 曲线对比 (a)修正模型计算结果; (b) J-A模型计算结果; (c) Z-L模型计算结果; (d)不同模型计算结果与试验结果相关系数对比Fig.1.Comparison of initial Be - H curves calculated by different models: (a) Our theoretical model; (b) J-A model; (c) Z-L model;(d) correlation coefficients of different models.

为定量地说明上述修正模型、Z-L模型和JA模型对单次磁化曲线的计算精度,验证修正模型的改进效果,选用相关系数R2作为模型拟合度优劣的评价标准,其计算公式为

式中,Pi为试验结果,为模型计算结果,i表示测试数据编号;n为数据点数量.R2值越接近1,则表明模型拟合度越高.不同模型计算结果的相关系数R2如表1和图1(d)所示.

表1 不同模型的相关系数 R2 比较Table 1.Correlation coefficients R2 of initial magnetization curve predicated by different models.

结果表明: 相对于J-A模型,Z-L模型和本文修正模型的理论预测结果的准确性都有极大的改进,模型计算得到的磁化曲线与试验测试曲线的相关系数R2远远高于J-A模型的计算结果,这种准确性的提高主要源于Z-L模型和本文修正模型都采用了更好的非线性磁致伸缩应变关系式; 同时,相对于Z-L模型,本文修正模型的准确性更好,不同应力条件下预测结果与试验结果之间的相关系数比Z-L模型的相关系数要高0.05以上.

3.2 加载条件下应力对磁滞回线的影响

利用Makar等[13]在不同含碳量的碳钢棒材试件上进行的试验进一步验证修正模型的改进情况.由于,Makar等进行的试验包含了弹-塑性两个阶段,而原Z-L模型只能描述弹性阶段应力的影响,本节中只将修正模型、J-A模型计算结果与Makar等的试验结果进行对比.选择含碳量为0.003%和0.15%两种碳钢棒材试件的测试结果进行验证,修正模型及J-A模型的部分参数可根据参考文献[13]中的试验数据直接得到,如: 含碳量为0.003%试件的屈服强度σs=90 MPa,平均晶粒直径d=1.25×10−3m,λws=3.85×10−6m,Mws=1.25×106A·m–1,−Ms=1.45×106A·m–1,所受应力分别为0,33和160 MPa; 含碳量为0.15%试件的屈服强度σs=122 MPa,平均晶粒直径d=1.9×10−4m,λws=4.17×10−6m,Mws=1.35×106A·m–1,Ms=1.6×106A·m–1,所受应力分别为0,36和182 MPa; 修正模型其他参数则需要根据磁化曲线通过拟合得到,如: 当含碳量为0.003%时,

当含碳量为0.15%时,c=0.25 ,k0=4.522×10−3A·m−1,其他参数与含碳量为0.003%时相同.J-A模型中磁致伸缩系数关系式中的参数可根据参考文献[13]中给出的磁致伸缩应变曲线拟合得到: 当含碳量为0.003%时,

当含碳量为0.15%时,

其他参数与修正模型中的对应参数取值一致.将修正模型和J-A模型计算得到的磁滞回线与文献[13]中的图3的试验结果进行对比,结果如图2所示.

图2结果表明: 随着应力的增大,磁滞回线均会产生顺时针转动,且塑性应力对磁滞回线的影响更为显著; 随着含碳量的增大,磁滞损耗增大,磁滞回线面积增大; 此外,在弹性应力(0.003 wt% C,33 MPa; 0.15 wt% C,36 MPa)作用下,铁磁材料的饱和磁化强度变化很小,而塑性应力(0.003 wt%C,160 MPa; 0.15 wt% C,182 MPa)作用下,铁磁材料的饱和磁化强度却显著降低,其原因在于当材料发生塑性变形时,材料内部的位错、晶格滑移增大,钉扎点数量急剧增多,使得材料磁化过程中的磁滞损耗显著增大.对比不同模型计算结果和试验结果可知: 在弹性应力范围内,J-A模型能够较好地反映应力对磁滞回线的影响趋势,但当进入塑性阶段后,J-A模型计算结果与试验结果存在明显差距; 而修正模型能够同时反映弹性阶段和塑性阶段应力对不同含碳量试件磁滞回线的影响规律,不同应力条件下理论计算得到的磁滞回线与试验测试得到的磁滞回线基本重合.

进一步,计算不同模型理论预测的磁滞回线与试验曲线的相关系数R2,结果如表2所列.

通过表2的对比可知: 在加载条件下,修正模型计算得到的磁滞回线与试验结果相关系数均在0.98以上,明显高于J-A模型的相关系数,能够更准确地预测弹性阶段和塑性阶段不同应力条件下的磁滞回线.

3.3 卸载条件下残余塑性变形对磁滞回线的影响

在上述试验基础上,Makar和Tanner[14]进一步分析了在塑性载荷卸载后碳钢棒材试件存在的残余塑性变形对磁滞回线的影响,同样选择含碳量为0.003%和0.15%两种试件的分别加载至160 MPa和182 MPa再卸载后测试得到的磁滞回线进行验证,结果如图3所示.

由图3试验结果可知: 当不存在残余塑性变形时,即σ=0 MPa,不同含碳量试件的磁滞回线均为标准的“S”型; 当存在残余塑性变形时,为保证材料内部晶粒间的总体应力平衡,试件内部将产生残余压应力,从而使得磁滞回线的形状显著改变:磁滞回线发生顺时针旋转,磁滞回线面积增大,饱和磁化强度明显降低,低磁场条件下磁滞回线的斜率显著减小,说明残余塑性变形增加了磁化难度,使得试件更难磁化.对比理论计算结果和试验数据可知: 与J-A模型相比,修正模型的计算结果与试验结果更为吻合,能够更好地反映残余塑性变形对材料磁滞回线的影响.计算不同模型预测的磁滞回线与试验曲线的相关系数,结果如表3所列.

图2 不同应力条件下修正模型和J-A模型磁滞回线计算结果的对比 (a)含碳量为0.003%时修正模型计算结果; (b)含碳量为0.003%时J-A模型计算结果; (c)含碳量为0.15%时修正模型计算结果; (d)含碳量为0.15%时J-A模型计算结果Fig.2.Hysteresis loops predicated by modified model and J-A model under different loading stresses: (a) Our modified model for 0.003 wt% C sample; (b) J-A model for 0.003 wt% C sample; (c) our modified model for 0.15 wt% C sample; (d) J-A model for 0.15 wt% C sample.

表2 加载条件下不同模型计算得到的磁滞回线与试验曲线相关系数 R2 比较Table 2.Correlation coefficients R2 of hysteresis loops predicated by different models under loading condition.

图3 不同残余塑性变形条件下修正模型和J-A模型磁滞回线计算结果的对比 (a)含碳量0.003%时修正模型计算结果; (b)含碳量0.003%时J-A模型计算结果; (c)含碳量0.15%时修正模型计算结果; (d)含碳量0.15%时J-A模型计算结果Fig.3.Hysteresis loops predicated by modified model and J-A model under different residual plastic deformation: (a) Our modified model for 0.003 wt% C sample; (b) J-A model for 0.003 wt% C sample; (c) our modified model for 0.15 wt% C sample; (d) J-A model for 0.15 wt% C sample.

表3 卸载条件下不同模型计算得到的磁滞回线与试验曲线相关系数 R2 比较Table 3.Correlation coefficients R2 of hysteresis loops predicated by different models under different residual plastic deformation.

表3数据表明: 卸载后,修正模型计算得到的磁滞回线与试验测试得到的磁滞回线都吻合得非常好,相关系数均在0.98以上,明显高于J-A模型计算得到的相关系数.

综上可以看出: 相比于J-A模型和Z-L模型,本文的修正磁化模型能够更准确地预测不同材料在不同应力、应变条件下的磁滞回线变化规律,具有更好的适用性.

4 应力、塑性变形对磁滞回线的影响

由第3节的分析可知,在弹性变形和塑性变形阶段,磁滞曲线形状将发生显著改变.为更好地掌握弹性拉、压应力和塑性拉、压应变对材料磁化过程的不同影响,本节进一步进行具体讨论.

4.1 弹性应力对磁滞回线的影响

在弹性应阶段,拉伸应力取σ=0,50,100,150和200 MPa,压应力取σ=–50,–100,–150和–200 MPa,不同应力条件下磁滞回线、矫顽力Hc和剩余磁感应强度Br随外磁场的变化规律如图4所示.

由图4(a)和图4(b)可知: 在弹性阶段,随着拉、压应力的增大,初始磁化曲线斜率降低,磁滞回线发生顺时针旋转,饱和磁化强度将发生极小幅度的衰减,磁滞回线所包围的面积增大,说明应力增大了材料磁化的难度; 相比于拉应力,压应力条件下的上述变化更为明显,说明压应力对材料磁化的阻碍作用更强.由图4(c)和图4(d)可知:Hc随着拉应力的增大呈非线性增大,而随着压应力的增大呈线性增大,相同幅值条件下拉应力对应的Hc远远小于压应力;Br随拉应力的增大先小幅增大后迅速减小,而随压应力的增大则是呈非线性减小趋势.需要指出的是,Hc和Br随弹性拉应力的上述变化规律与文献[13]中图5、图6中试验得到的变化规律是一致的.

4.2 塑性应变对磁滞回线的影响

当存在塑性变形时,拉伸塑性变形分别取εr=0,0.001 ,0.002 ,0.003和0.004 ,压缩塑性变形分别取εr=−0.001 ,−0.002 ,−0.003和−0.004 ,不同塑性变形条件下,磁滞回线、矫顽力Hc和剩余磁感应强度Br随外磁场的变化如图5所示.

图5 塑性应变对磁滞回线、矫顽力及剩余磁化强度的影响 (a)拉伸塑性应变对磁滞回线的影响; (b)压缩塑性变形对磁滞回线的影响; (c)塑性变形对矫顽力的影响; (d)塑性变形对剩余磁感应强度的影响Fig.5.Effects of plastic deformation on hysteresis loops,coercivity and remanence: (a) Effect of plastic tensile deformation on hysteresis loop; (b) effect of plastic compressive deformation on hysteresis loop; (c) effect of plastic tensile and compressive deformation on coercivity; (d) effect of plastic tensile and compressive deformation on remanence.

由图5(a)和图5(b)可知: 当发生塑性变形时,材料内部的位错、晶格滑移增大,钉扎点数量急剧增多,增大了材料磁化的阻力,导致初始磁化曲线的斜率大幅降低,磁滞回线形状将发生显著改变且沿顺时针旋转; 随着塑性变形的增大,磁滞回线面积增大,饱和磁化强度显著降低; 相比于拉伸塑性变形,压缩塑性变形条件下的上述变化更为明显,说明压缩塑性变形对材料磁化的阻碍作用更强; 此外,较小的塑性变形都将极大地改变磁化曲线形状,这也与文献[31]得到的结论相符.由图5(c)和图5(d)可知:Hc随着塑性变形的增大呈非线性而增大,而Br随着塑性变形的增大呈非线性减小,且压缩塑性变形条件下Hc和Br的变化幅度明显大于拉伸塑性变形.需要特别指出的是: 参考文献[28]中改进的J-A模型计算分析结果显示Hc和Br均是随塑性变形的增大而增大; 而对比Makar等[13,14]的试验结果可以看出: 在塑性变形条件下,矫顽力随塑性变形的增大呈线性增大趋势,而剩余磁化强度随塑性的增大呈非线性减小的趋势.因此,文献[28]改进的J-A模型计算得到的剩余磁化强度随塑性变形的变化规律与Makar等的试验结果相比存在偏差,而本文修正模型所反映的Hc和Br随塑性变形的变化趋势更符合试验结果.

5 结 论

本文针对J-A模型和Z-L的不足,在结合两种模型各自优势的基础上,建立了能够更好反映弹-塑性应力/应变对铁磁材料磁化曲线影响的修正磁化模型,通过与试验结果及原有模型计算结果进行对比,验证了本文修正磁化模型的有效性、准确性; 在此基础上,进一步分析了弹性拉、压应力及塑性拉、压应变对磁化曲线、矫顽力和剩余磁化强度的影响规律.

1)在单次磁化条件下: 在低磁场阶段,拉、压应力将阻碍材料磁化,使应力条件下磁化强度小于无应力条件下的磁化强度,且压应力对磁化的阻碍作用强于拉应力; 在高磁场阶段,拉、压应力对磁化强度强度的影响将发生逆转,使得压应力条件下的磁化强度高于拉应力条件下的磁化强度.

2)在循环磁化条件下: 随着应力和塑性应变的增大,初始磁化曲线斜率大幅降低,磁滞回线沿顺时针旋转,磁滞回线面积增大,饱和磁化强度降低; 相比于拉应力和拉伸塑性变形,压应力和压缩塑性变形对磁化曲线、Hc和Br的影响更为明显,说明压应力和压缩塑性变形对材料磁化的阻碍作用更强; 此外,相对于弹性变形,塑性变形对磁化的影响更大,较小的塑性变形都将显著改变磁化曲线形状、Hc和Br,这为利用Hc和Br检测塑性变形提供了理论基础.

3)通过与试验结果及原有模型计算结果对比,结果表明: 相对于J-A模型和Z-L模型,本文的修正模型无论在定性还是定量上都能更准确地反映单次磁化、循环磁化过程中应力、塑性变形对磁化曲线的影响规律,不同应力条件下的理论预测结果与试验结果的相关系数均在0.98以上,高于J-A模型和Z-L模型的计算结果,可为分析力磁耦合效应对铁磁材料磁化过程的影响规律提供更准确的理论分析模型.