删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

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

本站小编 Free考研考试/2021-12-29

摘要:Jiles-Atherton (J-A)模型和Zheng Xiao-Jing-Liu Xing-En (Z-L)模型在分析应力对铁磁材料磁化的影响方面应用十分广泛. 目前, J-A模型中的磁致伸缩应变与应力和磁化强度的关系式采用Jiles给出的经典拟合公式, 该拟合公式中磁化强度的二次项和四次项系数与应力均为线性关系, 不能准确描述铁磁材料磁致伸缩系数随应力、磁化强度的非线性变化规律; Z-L模型中磁致伸缩应变与应力和磁化强度的关系式采用了双曲正切函数tanh(x), 更好地描述了铁磁材料磁致伸缩应变和磁化强度随应力的非线性变化规律, 但Z-L模型却没有考虑Weiss分子场、钉扎效应的作用, 且由于采用了基于弹性能的接近定理, 只能描述弹性应力对磁化过程的影响. 针对上述问题, 本文结合Z-L模型中的非线性磁致伸缩应变关系式以及J-A模型中的磁滞理论, 考虑弹性应力、塑性变形对模型参数的影响, 建立了能够反映弹-塑性阶段应力与塑性变形对铁磁材料磁化曲线影响的修正磁化模型, 分析了弹性拉、压应力及塑性拉、压变形对磁化曲线、矫顽力和剩余磁化强度的影响规律. 通过与试验结果及原有模型的计算结果进行对比, 发现修正模型能够更好地反映单次磁化、循环磁化过程中应力、塑性变形对磁化曲线的影响规律, 理论预测结果与试验结果之间的相关系数均在0.98以上, 可为分析力磁耦合效应对铁磁材料磁化影响规律提供更准确的理论模型.
关键词: 力磁耦合/
磁致伸缩/
磁滞回线/
Jiles-Atherton模型/
Zheng Xiao-Jing-Liu Xing-En模型

English Abstract


--> --> -->
应力-应变对铁磁材料磁化过程的影响规律, 即力磁耦合效应, 一直是国内外****的研究热点[19]. 随着应力-应变的改变, 特别是应力集中或塑性变形的出现, 铁磁材料内部的磁畴结构和钉扎点密度等会发生显著变化, 从而影响材料的磁化曲线、矫顽力、剩余磁化强度等磁化特性[1014]. 基于上述现象, 发展出多种磁无损检测技术, 如: 金属磁记忆检测技术、巴克豪森磁噪声技术、磁声发射技术. 建立准确的磁化模型来描述应力、塑性变形对铁磁材料磁化的影响规律, 是利用上述磁无损检测技术实现应力集中、塑性变形检测与反演的基础.
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模型的中经典磁致伸缩应变关系式; 此后, 文献[1921]针对原磁致伸缩系数关系式不能很好反映低应力条件下磁致伸缩系数的变化规律, 提出利用神经网络来拟合磁致伸缩系数与磁化强度关系. 为描述拉应力和压应力对磁化特性影响的非对称性, Li和Xu [22]以及 Sablik[23]在原有模型基础上引入退磁项, 并基于Gaussian方程建立了饱和磁致伸缩系数与应力的变化关系; Li和Jiles [24]则利用Rayleigh定理对J-A模型中的接近定理进行了修正; 针对循环应力作用下零应力附近的磁化强度不符合逐渐趋于稳定状态的问题, Sablik等[25]提出了两种方法对接近定理进行改进. 为实现利用J-A模型描述塑性变形对材料磁化的影响, 文献[3, 2633]分析了塑性变形对材料位错密度的影响, 建立了钉扎系数、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\left( {\sigma,{M_{{\rm{an}}}},T} \right)$满足[35]
${\rm{d}}U = {\rm{d}}Q - {\rm{d}}W = T{\rm{d}}S + \sigma {\rm{d}}\varepsilon + {\mu _0}H{\rm{d}}{M_{{\rm{an}}}}, $
式中, $Q$为材料从外部吸收的热能; $W$为外力对材料做所的功; T为外界温度; $S$为体积熵的变化量; $\sigma $为应力; $\varepsilon $为应变; ${\mu _0}$为真空磁导率, ${\mu _0} = 4{\text{π}} \times {10^{ - 7}}$; $H$为外磁场强度; ${M_{{\rm{an}}}}$为非磁滞磁化强度.
根据可逆体系的热力学平衡原理, 理想磁化过程中铁磁材料的Gibbs自由能${G_{\rm{m}}}\left( {\sigma,{M_{{\rm{an}}}},T} \right)$可表示为[36]
$\begin{split}{G_{\rm{m}}}\left( {\sigma,{M_{{\rm{an}}}},T} \right) = &\; U + {{{\mu _0}\alpha {M_{{\rm{an}}}}^2} / 2}\\ & - {{{\mu _0}{N_{\rm{d}}}{M_{{\rm{an}}}}^2} / 2} - TS - \sigma \varepsilon,\end{split}$
式中, $\alpha $为Weiss分子场耦合系数, ${N_{\rm{d}}}$为退磁系数. 则${G_{\rm{m}}}\left( {\sigma,{M_{{\rm{an}}}},T} \right)$的全微分可写为
$\begin{split}{\rm{d}}{G_m}\,& = {\rm{d}}U - T{\rm{d}}S - S{\rm{d}}T - \sigma {\rm{d\varepsilon }} - \varepsilon {\rm{d\sigma }}\\ & = {\mu _0}\left( {H + \alpha {M_{{\rm{an}}}} - {N_{\rm{d}}}{M_{{\rm{an}}}}} \right){\rm{d}}{M_{{\rm{an}}}} - \varepsilon {\rm{d\sigma }} - S{\rm{d}}T.\end{split}$
当不考虑温度变化时, 即${\rm{d}}T = 0$, 则可得到只包含应力$\sigma $与非磁滞磁化强度${M_{{\rm{an}}}}$的Gibbs自由能${G_{\rm{m}}}\left( {\sigma,{M_{{\rm{an}}}}} \right)$全微分表达式, 即
${\rm{d}}{G_{\rm{m}}} = {\mu _0}\left( {H + \alpha {M_{{\rm{an}}}} - {N_{\rm{d}}}{M_{{\rm{an}}}}} \right){\rm{d}}{M_{{\rm{an}}}} - \varepsilon {\rm{d\sigma }}.$
由上述热力学关系可得到非磁滞条件下铁磁材料应变$\varepsilon \left( {\sigma,{M_{{\rm{an}}}}} \right)$和有效场${H_{\rm{e}}}\left( {\sigma,{M_{{\rm{an}}}}} \right)$关系为[3436]
$\varepsilon \left( {\sigma,{M_{{\rm{an}}}}} \right) = - \frac{{\partial {G_{\rm{m}}}}}{{\partial \sigma }},$
$\begin{split}{H_{\rm{e}}}\left( {\sigma,{M_{{\rm{an}}}}} \right) =\,& \frac{1}{{{\mu _0}}}\frac{{\partial {G_{\rm{m}}}}}{{\partial {M_{{\rm{an}}}}}} = H + \alpha {M_{{\rm{an}}}} - {N_{\rm{d}}}{M_{{\rm{an}}}} \\ &- \frac{1}{{{\mu _0}}}\frac{{\partial \left( {\int {\varepsilon \left( {\sigma,{M_{{\rm{an}}}}} \right){\rm{d\sigma }}} } \right)}}{{\partial {M_{{\rm{an}}}}}}.\end{split}$
为获得多项式形式的本构关系, 将${G_{\rm{m}}}\left( {\sigma,{M_{{\rm{an}}}}} \right)$$\left( {\sigma,\;{M_{{\rm{an}}}}} \right) = \left( {0,\;0} \right)$进行泰勒展开, 并根据试验现象对泰勒展开式进行简化. 根据Kurzar和Cullity[38], Yamasaki等[39]的试验结果及Jiles、周浩淼等的理论分析, 铁磁材料的磁致伸缩曲线是关于$y$轴对称的偶函数, 因此磁致伸缩应变$\varepsilon \left( {\sigma,{M_{{\rm{an}}}}} \right)$中的耦合项只保留含有$M_{\rm{an}}^2$$M_{\rm{an}}^4$项, 则

$\left\{ {\begin{aligned}& {\varepsilon \left( {\sigma,{M_{{\rm{an}}}}} \right) = \frac{\sigma }{{{E_{\rm{s}}}}} + {\lambda _0}\left( \sigma \right) + {\lambda _{\max }}\left( \sigma \right){{\left( {\frac{{{M_{{\rm{an}}}}}}{{{M_{{\rm{ws}}}}}}} \right)}^2} - \vartheta {\lambda _{{\rm{ws}}}}\frac{{\left( {{M_{{\rm{an}}}}^4 - {M_0}{{\left( \sigma \right)}^4}} \right)}}{{{M_{{\rm{ws}}}}^4}},} \\ & {{H_{{\rm{e}}\left( {{M_{{\rm{an}}}}} \right)}} = H + \alpha {M_{{\rm{an}}}} - {N_{\rm{d}}}{M_{{\rm{an}}}} + \frac{1}{{{\mu _0}}}\Bigg[ {\frac{{{M_{{\rm{an}}}}}}{{{M_{{\rm{ws}}}}^2}}\int_0^\sigma {{\lambda _{\max }}\left( \sigma \right){\rm{d}}\sigma } - 4\vartheta \sigma {\lambda _{{\rm{ws}}}}\frac{{\left( {{M_{{\rm{an}}}}^3 - {M_0}{{\left( \sigma \right)}^3}} \right)}}{{{M_{{\rm{ws}}}}^4}}} \Bigg],} \end{aligned}} \right.$

式中, ${E_{\rm{s}}}$为材料固有弹性模量; ${\lambda _0}\left( \sigma \right)$为应力单独作用下磁畴移动导致的磁致伸缩应变量; ${\lambda _{\max }}\left( \sigma \right)$为应力单独作用时由于磁畴移动导致的最大磁致伸缩应变量; ${M_{{\rm{ws}}}}$为无应力时的饱和壁移磁化强度; ${\lambda _{{\rm{ws}}}}$为无应力条件下材料的饱和磁致伸缩应变; ${M_0}\left( \sigma \right)$为应力作用时的饱和壁移磁化强度; $\vartheta $为阶跃函数, 当${M_{{\rm{an}}}} < {M_0}\left( \sigma \right)$时, $\vartheta = 0$, 当${M_{{\rm{an}}}} \geqslant {M_0}\left( \sigma \right)$时, $\vartheta = n$, 以保证只有当磁化强度${M_{{\rm{an}}}}$超过${M_0}\left( \sigma \right)$后, 磁化强度的四次项${M_{{\rm{an}}}}^4$才对磁致伸缩应变产生影响. 在文献[35, 36]理论分析的基础上, 结合文献[32, 33]对J-A模型的讨论, ${M_{{\rm{an}}}}$, ${M_0}\left( \sigma \right)$, ${\lambda _0}\left( \sigma \right)$${\lambda _{\max }}\left( \sigma \right)$可表示为

$\left\{ {\begin{aligned}& {{M_{{\rm{an}}}} = {M_{\rm{s}}}\left( {\coth \left( {\frac{{{H_{{\rm{e}}\left( {{M_{{\rm{an}}}}} \right)}}}}{a}} \right) - \frac{a}{{{H_{{\rm{e}}\left( {{M_{{\rm{an}}}}} \right)}}}}} \right),} \\ &{{M_0}\left( \sigma \right) = \frac{3}{4}{M_{{\rm{ws}}}}\left[ {{\rm{1}} - {\rm{tanh}}\left( {\beta \frac{\sigma }{{{\sigma _0}}} - \rm{arctanh} \left( {{1 / 3}} \right)} \right)} \right],} \\ & {{\lambda _0}\left( \sigma \right) = {\lambda _{{\rm{ws}}}}\left[ {\frac{3}{4}{\rm{tanh}}\left( {\beta \frac{\sigma }{{{\sigma _0}}} - \rm{arctanh} \left( {{1 / 3}} \right)} \right) + \frac{1}{4}} \right],} \\ & {{\lambda _{\max }} = {\lambda _{{\rm{ws}}}} - {\lambda _0}\left( \sigma \right) = \frac{3}{4}{\lambda _{{\rm{ws}}}}\left[ {{\rm{1}} - {\rm{tanh}}\left( {\beta \frac{\sigma }{{{\sigma _0}}} - \rm{arctanh} \left( {{1 / 3}} \right)} \right)} \right],} \\ & {\int_0^\sigma {{\lambda _{\max }}\left( \sigma \right){\rm{d}}\sigma } = \frac{3}{4}{\lambda _{{\rm{ws}}}}\left[ {\sigma - \frac{{{\sigma _0}}}{\beta }{\rm{ln cosh}}\left( {\beta \frac{\sigma }{{{\sigma _0}}} - \rm{arctanh} \left( {{1 / 3}} \right)} \right)} \right],} \end{aligned}} \right.$

其中, ${M_{\rm{s}}}$为饱和磁化强度; $a$为形状系数; ${\rm{tanh}}\left( x \right)$是双曲正切函数; ${\sigma _0}$表示将${\lambda _0}\left( \sigma \right)$在原点附近线性化为斜率为${1 / {{E'}}}$的直线时, 使${\lambda _0}\left( \sigma \right) = {\lambda _{{\rm{ws}}}}$时的应力值, ${1 / {{E'}}} = {1 / {{E_0}}} - {1 / {{E_{\rm{s}}}}}$, ${E_0}$为初始弹性模量; $\beta $为比例系数, 影响曲线斜率. 周浩淼改进的Z-L模型采用了分段函数来表示压应力和拉应力对${\lambda _0}\left( \sigma \right)$的不同影响, 为简化公式, 本文通过平移变换, 将Z-L模型中${\lambda _0}\left( \sigma \right)$的表达式变为了连续函数, 其中$ - \rm{arctanh} \left( {{1 / 3}} \right)$的引入是为保证当应力$\sigma = $0 MPa时, ${\lambda _0}\left( \sigma \right) = 0$.
将(8)式代入(7)式中, 即可得到非磁滞条件下材料磁致伸缩应变和磁化强度随应力、磁场的变化规律.
2
2.2.磁滞效应对磁化强度的影响
-->实际磁化过程中, 材料内部不可避免地存在缺陷、位错、晶格滑移等钉扎点, 从而阻碍磁畴运动. 钉扎效应的存在, 导致材料在磁化过程中不可避免地存在能量损耗. 因此, 根据能量守恒原理, 实际磁化过程中的磁化能量等于无磁滞磁化能减去钉扎点造成的磁滞损耗[33], 即
$\begin{split} &\int {M{\rm{d}}{B_{\rm{e}}}\left( M \right)} \\ =\, & \int {{M_{{\rm{an}}}}{\rm{d}}{B_{\rm{e}}}\left( M \right)} - {k_{{\rm{eff}}}}\delta \int {\left( {\frac{{{\rm{d}}{M_{{\rm{irr}}}}}}{{{\rm{d}}{B_{{\rm{e}}\left( M \right)}}}}} \right)} {\rm{d}}{B_{\rm{e}}}\left( M \right),\end{split}$
式中, $M$为实际磁化强度, $M = {M_{{\rm{rev}}}} + {M_{{\rm{irr}}}}$, ${M_{{\rm{irr}}}}$为不可逆磁化强度, ${M_{{\rm{rev}}}}$为可逆磁化强度; ${k_{{\rm{eff}}}}$为有效钉扎系数; $\delta $为方向系数, 当${{{\rm{d}}{H_{\rm{e}}}\left( M \right)} / {{\rm{d}}t}} > 0$时, $\delta = 1$, 当${{{\rm{d}}{H_{\rm{e}}}\left( M \right)} / {{\rm{d}}t}} < 0$时, $\delta = - 1$, 以保证钉扎点的作用总是阻碍磁场的变化; ${B_{\rm{e}}}\left( M \right)$表示实际磁感应强度, ${B_{\rm{e}}}\left( M \right) = {\mu _0}{H_{\rm{e}}}\left( M \right)$; ${H_{\rm{e}}}\left( M \right)$为实际磁化强度条件下的等效磁场强度, ${H_{\rm{e}}}\left( M \right)$可将(7)式中的${M_{{\rm{an}}}}$换成$M$得到, 即

${H_e}\left( M \right) \!=\! H + \alpha M \!-\! {N_d}M \!+ \frac{1}{{{\mu _0}}}\!\left( {\frac{3}{2}{\lambda _{\rm ws}}\left\{ {\sigma \!-\! \frac{{{\sigma _0}}}{\beta }{\rm{ln}}\left[ {\cosh \left( {\beta \frac{\sigma }{{{\sigma _0}}} \!-\! \rm{arctanh} \left( {{\dfrac13}} \right)} \right)} \right]} \right\}\frac{M}{{{M_{{\rm{ws}}}}^2}} \!-\! 4\vartheta \sigma {\lambda _{{\rm{ws}}}}\frac{{\left( {{M^3} \!-\! {M_0}{{\left( \sigma \right)}^3}} \right)}}{{{M_{{\rm{ws}}}}^4}}} \right). $
整理(9)式可得
$\frac{{{\rm{d}}{M_{{\rm{irr}}}}}}{{{\rm{d}}{H_{\rm{e}}}\left( M \right)}} = \frac{{{\mu _0}}}{{{k_{{\rm{eff}}}}\delta }}\left( {{M_{{\rm{an}}}} - M} \right).$
$M = {M_{\rm rev}} + {M_{\rm irr}}$, ${M_{{\rm{rev}}}} = c\left( {{M_{{\rm{an}}}} - {M_{{\rm{irr}}}}} \right)$可得
$M = c{M_{{\rm{an}}}} + \left( {1 - c} \right){M_{{\rm{irr}}}}.$
(12)式对外磁场强度$H$取微分可得
$\frac{{{\rm{d}}{M_{{\rm{irr}}}}}}{{{\rm{d}}H}} = \frac{1}{{1 - c}}\left( {\frac{{{\rm{d}}M}}{{{\rm{d}}H}} - c\frac{{{\rm{d}}{M_{{\rm{an}}}}}}{{{\rm{d}}H}}} \right).$
(10)式对外磁场强度$H$取微分可得

$\frac{{{\rm{d}}{H_{\rm{e}}}}}{{{\rm{d}}H}} \!=\! 1 \!+\! \left[ {\left( {\alpha \!-\! {N_{\rm{d}}}} \right) \!+\! \frac{1}{{{\mu _0}}}\left( {\frac{3}{2}{\lambda _{{\rm{ws}}}}\left\{ {\sigma \!-\! \frac{{{\sigma _0}}}{\beta }{\rm{ln}}\left[ {\cosh \left( {\beta \frac{\sigma }{{{\sigma _0}}} \!-\! \rm{arctanh} ( {{\dfrac13}} )} \right)} \right]} \right\}\frac{1}{{{M_{{\rm{ws}}}}^2}} \!-\! 12\vartheta \sigma {\lambda _{{\rm{ws}}}}\frac{{\left( {{M^2} \!-\! {M_0}{{( \sigma )}^2}} \right)}}{{{M^4_{{\rm{ws}}}}}}} \right)} \right]\frac{{{\rm{d}}M}}{{{\rm{d}}H}}. $
联立公(11), (13)和(14)式, 并整理可得
$\begin{split}\dfrac{{{\rm{d}}M}}{{{\rm{d}}H}} ={}& \bigg[\dfrac{{{\mu _0}}}{{{k_{{\rm{eff}}}}\delta }}\left( {{M_{{\rm{an}}}} - M} \right) + \dfrac{c}{{1 - c}}\dfrac{{{\rm{d}}{M_{{\rm{an}}}}}}{{{\rm{d}}H}}\bigg]\times\Bigg\{\dfrac{1}{{1 - c}} - \dfrac{{{\mu _0}}}{{{k_{{\rm{eff}}}}\delta }}\left( {{M_{{\rm{an}}}} - M} \right) \bigg[\left( {\alpha - {N_{\rm{d}}}} \right) + \dfrac{1}{\mu _0} \\ & \times \left( \dfrac{3}{2}{\lambda _{{\rm{ws}}}}\left\{ {\sigma - \dfrac{{{\sigma _0}}}{\beta }{\rm{ln}}\left[ {\cosh \left( {\beta \dfrac{\sigma }{{{\sigma _0}}} - \rm{arctanh} \left( {{\dfrac13}} \right)} \right)} \right]} \right\}\dfrac{1}{{{M_{{\rm{ws}}}}^2}} - 12\vartheta \sigma {\lambda _{{\rm{ws}}}}\dfrac{{\left( {{M^2} - {M_0}{{\left( \sigma \right)}^2}} \right)}}{{{M_{{\rm{ws}}}}^4}} \right) \bigg]\Bigg\}^{-1}.\end{split}$
其中, ${{{\rm{d}}{M_{{\rm{an}}}}} / {{\rm{d}}H}}$可联立(7)和(8)式得到, 即
$\begin{split}\dfrac{{{\rm{d}}{M_{{\rm{an}}}}}}{{{\rm{d}}H}} =\;& - {M_{\rm{s}}}\left[ {\dfrac{1}{a}{{\rm{csch} }^2}\left( {\dfrac{{{H_{\rm{e}}}}}{a}} \right) - \dfrac{a}{{{H_{\rm{e}}}^2}}} \right]\times \Bigg\{ 1 + {M_{\rm{s}}}\left( {\dfrac{1}{a}{{\rm{csch} }^2}\left( {\dfrac{{{H_{\rm{e}}}}}{a}} \right) - \dfrac{a}{{{H_{\rm{e}}}^2}}} \right)\Bigg[ \left( {\alpha - {N_{\rm{d}}}} \right) + \dfrac{1}{{{\mu _0}}} \\& \times \left( {\dfrac{3}{2}{\lambda _{{\rm{ws}}}}\left\{ {\sigma \!-\! \dfrac{{{\sigma _0}}}{\beta }{\rm{ln}}\left[ {\cosh \left( {\beta \dfrac{\sigma }{{{\sigma _0}}} \!-\! \rm{arctanh} \left( {{\dfrac13}} \right)} \right)} \right]} \right\}\dfrac{1}{{{M_{{\rm{ws}}}}^2}} \!-\! 12\vartheta \sigma {\lambda _{{\rm{ws}}}}\dfrac{{\left( {{M^2_{{\rm{an}}}} \!-\! {M_0}{{\left( \sigma \right)}^2}} \right)}}{{{M_{{\rm{ws}}}}^4}}} \right) \Bigg] \Bigg\}^{-1}.\end{split}$

联立(15)和(16)式即可计算不同应力条件下磁化强度$M$随外磁场$H$的变化规律.
2
2.3.应力-应变对模型参数的影响
-->在应力、应变, 特别是塑性变形作用下, 材料内部磁畴结构的形状和尺寸以及钉扎点密度都会发生改变, 从而改变Weiss分子场耦合系数$\alpha $、有效钉扎系数${k_{{\rm{eff}}}}$以及形状系数$a$的取值.
在一定范围内, 应变硬化应力${\sigma _{\rm{F}}}$与塑性应变${\varepsilon _{\rm{r}}}$之间的关系可根据Ludwik定律得到[40]:
${\sigma _{\rm{F}}} = A{\varepsilon _{\rm{r}}}^m,$
式中$A$为常数; $m$为Ludwik指数, 与材料性质有关.
根据Astie-Degauque关系, 位错密度${\zeta _{\rm{d}}}$与材料的应变硬化应力${\sigma _{\rm{F}}}$呈比例关系[28]:
${\zeta _{\rm{d}}} = {\left[ {\left( {\frac{{\left| {{\sigma _{\rm{F}}}} \right|}}{{{\alpha _k}Gb}}} \right) + \zeta _{{\rm{d}}0}^{{1 / 2}}} \right]^2},$
式中, ${\zeta _{{\rm{d}}0}}$为无应力时材料初始位错密度; $G$为剪切模量, 其与固有弹性模量${E_{\rm{s}}}$和泊松比$\upsilon $之间的关系为$G = {{{E_{\rm{s}}}} / {2\left( {1 + \upsilon } \right)}}$; $ b$为柏氏矢量, 是描述位错实质的物理量; ${\alpha _k}$为常数, 取0.76.
Lo等[28,29]的试验结果表明: 拉应力使钉扎密度减小, 压应力使钉扎密度增大; 而位错会以钉扎的形式阻碍畴壁运动. 因此, 有效钉扎系数${k_{{\rm{eff}}}}$和形状系数$a$与应力和位错密度之间的关系可表示为[30]
${k_{{\rm{eff}}}} = \left( {{G_1} + \frac{{{G_2}}}{d}} \right)\zeta _{\rm{d}}^{{1 / 2}}{k_0} - \frac{3}{2}\eta {\lambda _{{\rm{ws}}}}\sigma,$
$a = \left( {{G_3} + \frac{{{G_4}}}{d}} \right)\zeta _{\rm{d}}^{{1 / 2}}{a_0},$
式中, $\left[ {{G_1} + {{{G_2}} / d}} \right]\zeta _{\rm{d}}^{{1 / 2}}{k_0}$表示塑性变形对钉扎系数的改变量, $ - {{3\eta {\lambda _{{\rm{ws}}}}\sigma } / 2}$表示$\sigma $${k_{{\rm{eff}}}}$的改变量; ${\zeta _{\rm{d}}}$为错密度, ${k_0}$为初始钉扎系数; ${G_1}$, ${G_2}$, ${G_3}$, ${G_4}$, $\eta $为拟合常数, 一般取${G_1} = {G_3}$, ${G_2} = {G_4}$; ${a_0}$为初始形状系数; $d$为材料晶粒尺寸. 在弹性阶段, 材料内部位错密度不会发生改变, 取$\left[ {{G_1} \!+\! {{{G_2}} / d}} \right]\zeta _{{{\rm{d}_{\rm{0}}}}}^{{1 / 2}}\!=\! 1$, $\left( {{G_3} + {{{G_4}} / d}} \right)\zeta _{{{\rm{d}}_0}}^{{1 / 2}} = 1$.
由于应力-应变改变了磁畴的形状和尺寸, Weiss分子场耦合系数$\alpha $会随着应力-应变的增加而发生变化, 其与应力-应变的关系为[23,31,41]
$\begin{split}\alpha =\, & {\alpha _0} + {q_1}\ln \left( {{q_2}{\varepsilon _{\rm{r}}} + {q_3}} \right) \\ &+ \frac{{21}}{4}\frac{{{E_{\rm{s}}}}}{{{\mu _0}}}{\left\{ {{{{\lambda _{{\rm{ws}}}}\exp \left[ { - \frac{1}{2}{{\left( {\frac{\sigma }{{{\sigma _0}}}} \right)}^2}} \right]} / {{M_{\rm{s}}}}}} \right\}^2},\end{split}$
式中, ${\alpha _0}$为初始Weiss分子场耦合系数; ${q_1}\ln \left( {{q_2}{\varepsilon _{\rm{r}}} + {q_3}} \right)$表示塑性变形对$\alpha $改变量;
$\frac{{21}}{4}\frac{{{E_{\rm{s}}}}}{{{\mu _0}}}{\left\{ {{{{\lambda _{{\rm{ws}}}}\exp \left[ { - \frac{1}{2}{{\left( {\frac{\sigma }{{{\sigma _0}}}} \right)}^2}} \right]} / {{M_{\rm{s}}}}}} \right\}^2}$
表示弹性应力对$\alpha $改变量; ${q_1}$, ${q_2}$, ${q_3}$均为拟合常数.
2
3.1.单次磁化过程中应力对磁化曲线的影响
-->为说明修正模型相对于Z-L模型、J-A模型的改进情况, 利用Jiles等[8]进行的–200, –100, 0, 100和200 MPa五组不同预应力条件下的磁化测试结果进行验证, 试验结果如文献[8]图3虚线所示. Jiles模型中参数与文献[8]中的图2一致; Z-L模型中相关参数取值与文献[35]中图9一致; 修正模型中$\beta = 1.4$, $\vartheta = 0.75$, ${\alpha _0} = 1.44 \times {10^{ - 3}}$, 其他参数与Z-L模型参数取值一致, 不同应力条件下不同模型计算结果与试验结果对比如图1所示.
图 1 单次磁化条件下不同模型计算得到的${B_{\rm{e}}}$-$H$曲线对比 (a)修正模型计算结果; (b) J-A模型计算结果; (c) Z-L模型计算结果; (d)不同模型计算结果与试验结果相关系数对比
Figure1. Comparison of initial${B_{\rm{e}}}$-$H$curves calculated by different models: (a) Our theoretical model; (b) J-A model; (c) Z-L model; (d) correlation coefficients of different models

由试验结果可知: 与无应力条件下的磁化曲线相比, 应力将阻碍材料的磁化, 且在低磁场阶段压应力对磁化的阻碍作用明显高于拉应力; 当外磁场增强至中高磁场阶段后, 压应力条件下的磁化曲线和拉应力条件下的磁化曲线将发生反转, 最终压应力条件下的磁化曲线将趋于自由状态下的磁化曲线, 而拉应力条件下的磁化曲线会比无应力条件下磁化小很多. 对比三种模型计算结果可知: J-A模型能够定性反映应力对磁化的影响, 但是定量上却存在一定差别, 例如试验结果的无应力饱和磁感应强度在1.6 T左右, 而J-A模型的无应力饱和磁感应强度达到了1.8 T左右; 相对于J-A模型, Z-L模型在定性和定量上都有很大的改进, 能够很好地反映应力对磁化影响, 但是定量上仍然存在一定差异; 而本文修正后的模型无论在定性上还是定量上都能十分准确地反映试验结果, 不同应力条件下的理论预测结果与试验结果基本重合.
为定量地说明上述修正模型、Z-L模型和J-A模型对单次磁化曲线的计算精度, 验证修正模型的改进效果, 选用相关系数${R^2}$作为模型拟合度优劣的评价标准, 其计算公式为
$S{S_{{\rm{tot}}}} = \sum\nolimits_1^n {{{\left( {{P_i} - \bar P} \right)}^2}},$
$S{S_{{\rm{err}}}} = \sum\nolimits_1^n {{{\left( {{{\hat P}_i} - {P_i}} \right)}^2}},$
${R^2} = 1 - S{S_{{\rm{err}}}}/S{S_{{\rm{tot}}}},$
式中, ${P_i}$为试验结果, ${\hat P_i}$为模型计算结果, $i$表示测试数据编号; $n$为数据点数量. ${R^2}$值越接近1, 则表明模型拟合度越高. 不同模型计算结果的相关系数${R^2}$表1图1(d)所示.
模型类型应力值/MPa
–200–1000100200
修正模型${R^2}$0.99040.99470.99730.98590.9724
Z-L模型${R^2}$0.94210.92900.93330.93000.8956
J-A模型${R^2}$0.62890.0394–0.02890.26880.6440


表1不同模型的相关系数${R^2}$比较
Table1.Correlation coefficients ${R^2}$ of initial magnetization curve predicated by different models.

结果表明: 相对于J-A模型, Z-L模型和本文修正模型的理论预测结果的准确性都有极大的改进, 模型计算得到的磁化曲线与试验测试曲线的相关系数${R^2}$远远高于J-A模型的计算结果, 这种准确性的提高主要源于Z-L模型和本文修正模型都采用了更好的非线性磁致伸缩应变关系式; 同时, 相对于Z-L模型, 本文修正模型的准确性更好, 不同应力条件下预测结果与试验结果之间的相关系数比Z-L模型的相关系数要高0.05以上.
2
3.2.加载条件下应力对磁滞回线的影响
-->利用Makar等[13]在不同含碳量的碳钢棒材试件上进行的试验进一步验证修正模型的改进情况. 由于, Makar等进行的试验包含了弹-塑性两个阶段, 而原Z-L模型只能描述弹性阶段应力的影响, 本节中只将修正模型、J-A模型计算结果与Makar等的试验结果进行对比. 选择含碳量为0.003%和0.15%两种碳钢棒材试件的测试结果进行验证, 修正模型及J-A模型的部分参数可根据参考文献[13]中的试验数据直接得到, 如: 含碳量为0.003%试件的屈服强度${\sigma _{\rm{s}}} = $ 90 MPa, 平均晶粒直径$d = 1.25 \times\! {10^{ - 3}}$ m, ${\lambda _{{\rm{ws}}}} = 3.85 \times \!{\rm{1}}{{\rm{0}}^{ - 6}}$ m, ${M_{{\rm{ws}}}} = $ $ 1.25 \times {10^6}$ A·m–1, $-{M_{\rm{s}}} = 1.45 \times {10^6} $ A·m–1, 所受应力分别为0, 33和160 MPa; 含碳量为0.15%试件的屈服强度${\sigma _{\rm{s}}} = $ 122 MPa, 平均晶粒直径$ d = $$ 1.9 \times {10^{ - 4}}$ m, ${\lambda _{{\rm{ws}}}} = 4.17 \times {\rm{1}}{{\rm{0}}^{ - 6}}$ m, ${M_{{\rm{ws}}}} = 1.35 \;\times$ 106 A·m–1, ${M_{\rm{s}}} = 1.6 \times {10^6}$ A·m–1, 所受应力分别为0, 36和182 MPa; 修正模型其他参数则需要根据磁化曲线通过拟合得到, 如: 当含碳量为0.003%时,
$\begin{split}& \beta = 5,\;\vartheta = 0.75,\;{\sigma _0} = 100\;{\rm{ MPa}},\; {a_0} = 900,\\ &{N_{\rm{d}}} = 5 \times {10^{ - 5}},\;{\alpha _0} = 1.44 \times {10^{ - 3}},\; c = 0.65,\\ &{k_0} = 1.52 \times {10^{ - 3}}\;{\rm{ A}}\cdot{{\rm{m}}^{{\rm{ - 1}}}},\;b = 2.5 \times {10^{ - 10}}\;{\rm{ m}},\\ & {G_1} = {G_3} = 0.81 \times {10^{ - 6}}\;{\rm{m}},\;{G_2} = {G_4} = 3.8 \times {10^{ - 12}}\;{{\rm{m}}^2},\\& A = 2.5 \times {10^{10}}\;{\rm{J}}\cdot{{\rm{m}}^{{\rm{ - 3}}}},\;{q_1} = - 2.15566 \times {10^{ - 4}},\\& {q_2} = 62.5447,\;{q_3} = 0.0618315, \;\eta = 1.5 \times {10^{ - 7}};\end{split}$
当含碳量为0.15%时, $c = 0.25$, ${k_0} = 4.522\; \times $$ {10^{ - 3}}\;{\rm{ A}}\cdot{{\rm{m}}^{{\rm{ - 1}}}}$, 其他参数与含碳量为0.003%时相同. J-A模型中磁致伸缩系数关系式中的参数可根据参考文献[13]中给出的磁致伸缩应变曲线拟合得到: 当含碳量为0.003%时,
$\begin{split}& {r_{1(0)}} = 6.05 \times {10^{ - 18}}\;{{\rm{A}}^{{\rm{ - 2}}}}\cdot{{\rm{m}}^{\rm{2}}},\\ &r_{1(0)}' = - 1 \times {10^{ - 25}}\;{{\rm{A}}^{{\rm{ - 2}}}}\cdot{{\rm{m}}^{\rm{2}}}\cdot{\rm{P}}{{\rm{a}}^{{\rm{ - 1}}}},\\ &{r_{2(0)}} = - 2.45 \times {10^{ - 30}}\;{{\rm{A}}^{{\rm{ - 4}}}}\cdot{{\rm{m}}^4},\\ &r_{_{2(0)}}' = 4.53 \times {10^{ - 38}}\;{{\rm{A}}^{{\rm{ - 4}}}}\cdot{{\rm{m}}^4}\cdot{\rm{P}}{{\rm{a}}^{{\rm{ - 1}}}};\end{split}$
当含碳量为0.15%时,
$\begin{split}& {r_{1(0)}} = 5.54 \times {10^{ - 18}}\;{{\rm{A}}^{{\rm{ - 2}}}}\cdot{{\rm{m}}^{\rm{2}}},\\ &r_{1(0)}' = - 3.04 \times {10^{ - 23}}\;{{\rm{A}}^{{\rm{ - 2}}}}\cdot{{\rm{m}}^{\rm{2}}}\cdot{\rm{P}}{{\rm{a}}^{{\rm{ - 1}}}},\\ &{r_{2(0)}} = - 1.99 \times {10^{ - 30}}\;{{\rm{A}}^{{\rm{ - 4}}}}\cdot{{\rm{m}}^4},\\ &r_{_{2(0)}}' = - 5.20 \times {10^{ - 38}}\;{{\rm{A}}^{{\rm{ - 4}}}}\cdot{{\rm{m}}^4}\cdot{\rm{P}}{{\rm{a}}^{{\rm{ - 1}}}};\end{split}$
其他参数与修正模型中的对应参数取值一致. 将修正模型和J-A模型计算得到的磁滞回线与文献[13]中的图3的试验结果进行对比, 结果如图2所示.
图 2 不同应力条件下修正模型和J-A模型磁滞回线计算结果的对比 (a)含碳量为0.003%时修正模型计算结果; (b)含碳量为0.003%时J-A模型计算结果; (c)含碳量为0.15%时修正模型计算结果; (d)含碳量为0.15%时J-A模型计算结果
Figure2. 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结果表明: 随着应力的增大, 磁滞回线均会产生顺时针转动, 且塑性应力对磁滞回线的影响更为显著; 随着含碳量的增大, 磁滞损耗增大, 磁滞回线面积增大; 此外, 在弹性应力(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模型计算结果与试验结果存在明显差距; 而修正模型能够同时反映弹性阶段和塑性阶段应力对不同含碳量试件磁滞回线的影响规律, 不同应力条件下理论计算得到的磁滞回线与试验测试得到的磁滞回线基本重合.
进一步, 计算不同模型理论预测的磁滞回线与试验曲线的相关系数${R^2}$, 结果如表2所列.
模型类型试件含碳量0.003 wt%试件含碳量0.15 wt%
0 MPa33 MPa160 MPa0 MPa36 MPa182 MPa
修正模型${R^2}$0.98450.98940.98980.98400.98080.9937
J-A模型${R^2}$0.91330.91880.89840.94850.94960.9153


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

通过表2的对比可知: 在加载条件下, 修正模型计算得到的磁滞回线与试验结果相关系数均在0.98以上, 明显高于J-A模型的相关系数, 能够更准确地预测弹性阶段和塑性阶段不同应力条件下的磁滞回线.
2
3.3.卸载条件下残余塑性变形对磁滞回线的影响
-->在上述试验基础上, Makar和Tanner [14]进一步分析了在塑性载荷卸载后碳钢棒材试件存在的残余塑性变形对磁滞回线的影响, 同样选择含碳量为0.003%和0.15%两种试件的分别加载至160 MPa和182 MPa再卸载后测试得到的磁滞回线进行验证, 结果如图3所示.
图 3 不同残余塑性变形条件下修正模型和J-A模型磁滞回线计算结果的对比 (a)含碳量0.003%时修正模型计算结果; (b)含碳量0.003%时J-A模型计算结果; (c)含碳量0.15%时修正模型计算结果; (d)含碳量0.15%时J-A模型计算结果
Figure3. 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试验结果可知: 当不存在残余塑性变形时, 即$\sigma = 0$ MPa, 不同含碳量试件的磁滞回线均为标准的“S”型; 当存在残余塑性变形时, 为保证材料内部晶粒间的总体应力平衡, 试件内部将产生残余压应力, 从而使得磁滞回线的形状显著改变: 磁滞回线发生顺时针旋转, 磁滞回线面积增大, 饱和磁化强度明显降低, 低磁场条件下磁滞回线的斜率显著减小, 说明残余塑性变形增加了磁化难度, 使得试件更难磁化. 对比理论计算结果和试验数据可知: 与J-A模型相比, 修正模型的计算结果与试验结果更为吻合, 能够更好地反映残余塑性变形对材料磁滞回线的影响. 计算不同模型预测的磁滞回线与试验曲线的相关系数, 结果如表3所列.
模型类型试件含碳量0.003 wt%试件含碳量0.153 wt%
0 MPa160 MPa0 MPa182 MPa
修正模型${R^2}$0.98450.99430.98400.9858
J-A模型${R^2}$0.91330.98260.94850.9765


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

表3数据表明: 卸载后, 修正模型计算得到的磁滞回线与试验测试得到的磁滞回线都吻合得非常好, 相关系数均在0.98以上, 明显高于J-A模型计算得到的相关系数.
综上可以看出: 相比于J-A模型和Z-L模型, 本文的修正磁化模型能够更准确地预测不同材料在不同应力、应变条件下的磁滞回线变化规律, 具有更好的适用性.
由第3节的分析可知, 在弹性变形和塑性变形阶段, 磁滞曲线形状将发生显著改变.为更好地掌握弹性拉、压应力和塑性拉、压应变对材料磁化过程的不同影响, 本节进一步进行具体讨论.
2
4.1.弹性应力对磁滞回线的影响
-->在弹性应阶段, 拉伸应力取$\sigma = $0, 50, 100, 150和200 MPa, 压应力取$\sigma = $–50, –100, –150和–200 MPa, 不同应力条件下磁滞回线、矫顽力${H_{\rm{c}}}$和剩余磁感应强度${B_{\rm{r}}}$随外磁场的变化规律如图4所示.
图 4 弹性应力对磁滞回线、矫顽力及剩余磁感应强度的影响 (a)弹性拉应力对磁滞回线的影响; (b)弹性压应力对磁滞回线的影响; (c)弹性拉、压应力对矫顽力的影响; (d)弹性拉、压应力对剩余磁感应强度的影响
Figure4. Effects of elastic stress on hysteresis loops, coercivity and remanence: (a) Effect of elastic tensile stress on hysteresis loop; (b) effect of elastic compressive stress on hysteresis loop; (c) effect of elastic tensile and compressive stress on coercivity; (d) effect of elastic tensile and compressive stress on remanence

图4(a)图4(b)可知: 在弹性阶段, 随着拉、压应力的增大, 初始磁化曲线斜率降低, 磁滞回线发生顺时针旋转, 饱和磁化强度将发生极小幅度的衰减, 磁滞回线所包围的面积增大, 说明应力增大了材料磁化的难度; 相比于拉应力, 压应力条件下的上述变化更为明显, 说明压应力对材料磁化的阻碍作用更强. 由图4(c)图4(d)可知: ${H_{\rm{c}}}$随着拉应力的增大呈非线性增大, 而随着压应力的增大呈线性增大, 相同幅值条件下拉应力对应的${H_{\rm{c}}}$远远小于压应力; ${B_{\rm{r}}}$随拉应力的增大先小幅增大后迅速减小, 而随压应力的增大则是呈非线性减小趋势. 需要指出的是, ${H_{\rm{c}}}$${B_{\rm{r}}}$随弹性拉应力的上述变化规律与文献[13]中图5、图6中试验得到的变化规律是一致的.
2
4.2.塑性应变对磁滞回线的影响
-->当存在塑性变形时, 拉伸塑性变形分别取${\varepsilon _{\rm{r}}} = 0$, $0.001$, $0.002$, $0.003$$0.004$, 压缩塑性变形分别取${\varepsilon _{\rm{r}}} = - 0.001$, $ - 0.002$, $ - 0.003$$ - 0.004$, 不同塑性变形条件下, 磁滞回线、矫顽力${H_{\rm{c}}}$和剩余磁感应强度${B_{\rm{r}}}$随外磁场的变化如图5所示.
图 5 塑性应变对磁滞回线、矫顽力及剩余磁化强度的影响 (a)拉伸塑性应变对磁滞回线的影响; (b)压缩塑性变形对磁滞回线的影响; (c)塑性变形对矫顽力的影响; (d)塑性变形对剩余磁感应强度的影响
Figure5. 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)可知: ${H_{\rm{c}}}$随着塑性变形的增大呈非线性而增大, 而${B_{\rm{r}}}$随着塑性变形的增大呈非线性减小, 且压缩塑性变形条件下${H_{\rm{c}}}$${B_{\rm{r}}}$的变化幅度明显大于拉伸塑性变形. 需要特别指出的是: 参考文献[28]中改进的J-A模型计算分析结果显示${H_{\rm{c}}}$${B_{\rm{r}}}$均是随塑性变形的增大而增大; 而对比Makar等[13,14]的试验结果可以看出: 在塑性变形条件下, 矫顽力随塑性变形的增大呈线性增大趋势, 而剩余磁化强度随塑性的增大呈非线性减小的趋势. 因此, 文献[28]改进的J-A模型计算得到的剩余磁化强度随塑性变形的变化规律与Makar等的试验结果相比存在偏差, 而本文修正模型所反映的${H_{\rm{c}}}$${B_{\rm{r}}}$随塑性变形的变化趋势更符合试验结果.
本文针对J-A模型和Z-L的不足, 在结合两种模型各自优势的基础上, 建立了能够更好反映弹-塑性应力/应变对铁磁材料磁化曲线影响的修正磁化模型, 通过与试验结果及原有模型计算结果进行对比, 验证了本文修正磁化模型的有效性、准确性; 在此基础上, 进一步分析了弹性拉、压应力及塑性拉、压应变对磁化曲线、矫顽力和剩余磁化强度的影响规律.
1)在单次磁化条件下: 在低磁场阶段, 拉、压应力将阻碍材料磁化, 使应力条件下磁化强度小于无应力条件下的磁化强度, 且压应力对磁化的阻碍作用强于拉应力; 在高磁场阶段, 拉、压应力对磁化强度强度的影响将发生逆转, 使得压应力条件下的磁化强度高于拉应力条件下的磁化强度.
2)在循环磁化条件下: 随着应力和塑性应变的增大, 初始磁化曲线斜率大幅降低, 磁滞回线沿顺时针旋转, 磁滞回线面积增大, 饱和磁化强度降低; 相比于拉应力和拉伸塑性变形, 压应力和压缩塑性变形对磁化曲线、${H_{\rm{c}}}$${B_{\rm{r}}}$的影响更为明显, 说明压应力和压缩塑性变形对材料磁化的阻碍作用更强; 此外, 相对于弹性变形, 塑性变形对磁化的影响更大, 较小的塑性变形都将显著改变磁化曲线形状、${H_{\rm{c}}}$${B_{\rm{r}}}$, 这为利用${H_{\rm{c}}}$${B_{\rm{r}}}$检测塑性变形提供了理论基础.
3)通过与试验结果及原有模型计算结果对比, 结果表明: 相对于J-A模型和Z-L模型, 本文的修正模型无论在定性还是定量上都能更准确地反映单次磁化、循环磁化过程中应力、塑性变形对磁化曲线的影响规律, 不同应力条件下的理论预测结果与试验结果的相关系数均在0.98以上, 高于J-A模型和Z-L模型的计算结果, 可为分析力磁耦合效应对铁磁材料磁化过程的影响规律提供更准确的理论分析模型.
相关话题/塑性 材料 计算 文献 应力

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 重费米子材料与物理
    摘要:作为典型的强关联电子体系,重费米子材料表现出丰富的量子基态,如反铁磁序、铁磁序、非常规超导、非费米液体、自旋液体、轨道序和拓扑态等.相比其他强关联电子体系,重费米子体系的特征能量尺度低,可以通过压力、磁场或掺杂等参量对不同量子态进行连续调控,因而是研究量子相变、超导及其相互作用的理想体系.本文 ...
    本站小编 Free考研考试 2021-12-29
  • 半导体黄光发光二极管新材料新器件新设备
    摘要:在可见光范围内,半导体发光二极管(LED)发展很不平衡,黄光LED光效(光功率效率)长期远低于其他颜色光效.本文基于GaN/Si体系,从材料生长、芯片制造、器件物理和专用装备等方面进行了系统研究,解决了外延膜龟裂、位错过多、量子阱应力过大、InGaN黄光阱材料相分离、空穴浓度不足、阱材料生长温 ...
    本站小编 Free考研考试 2021-12-29
  • 钡作为掺杂元素调控铅基钙钛矿材料的毒性和光电特性
    摘要:近年来,有机-无机杂化卤化物钙钛矿材料(ABX3)由于具有优异的光电性质,受到了材料、能源等领域的广泛关注.但是,卤化物钙钛矿存在两个明显阻碍其商业化应用的问题:热稳定性差和含有有毒的铅(Pb)元素.相比于有机-无机杂化卤化物钙钛矿,全无机卤化物钙钛矿通常拥有更好的热稳定性.同时,采用一些无毒 ...
    本站小编 Free考研考试 2021-12-29
  • 有机-无机杂化钙钛矿材料的本征稳定性
    摘要:有机-无机杂化钙钛矿太阳能电池的光电转换效率已逾24%,效率的飞速提升加之可低成本溶液法制备的市场优势,使人们越来越期待钙钛矿太阳能电池的商业化.目前钙钛矿太阳能电池商业化所面临最大的障碍是材料乃至器件的长期不稳定性,这使其无法在使用寿命上与已商品化的硅基等太阳能电池匹敌.本文从化学不稳定性和 ...
    本站小编 Free考研考试 2021-12-29
  • 非互易旋电材料硅基矩形波导的色散特性研究
    摘要:研究设计了基于光通信C波段旋电材料的矩形波导,利用有效折射率法对波导有效折射率及横向电场分布进行求解,得到矩形波导中${m{E}}_{mn}^x$导模的色散方程.研究了在外磁场作用下表面磁等离子体激元的非互易传播特性.还研究了结构参数和材料折射率对非互易色散关系、时延特性的影响.结果表明: ...
    本站小编 Free考研考试 2021-12-29
  • 锂离子电池正极材料Li<sub>2</sub>FeO<sub>2</sub>的电子结构性质和Li扩散
    摘要:采用基于密度泛函理论的第一性原理方法计算了锂离子电池正极材料Immm-Li2FeO2的声子谱、电子结构性质和Li扩散系数并与Li2MO2(M=Co,Ni,Cu)材料进行对比.计算结果显示,Immm-Li2FeO2材料具有结构稳定性,计算结果呈铁磁性,能带结构具有半金属的特征.Fe离子外层d电子 ...
    本站小编 Free考研考试 2021-12-29
  • 金刚石-碳化硅超硬复合材料的冲击强度
    摘要:不同于延性介质,脆性介质的失效破坏严重制约着材料的强度.本文采用一种定量描述脆性介质力学性质的格点-弹簧模型,研究了金刚石-碳化硅超硬复合材料的冲击强度及其细观损伤机理,有助于避免灾变破坏、提高冲击强度.在模型中,通过构建不同体积分数比的金刚石和碳化硅两相复合材料,模拟获得了经受冲击波压缩形变 ...
    本站小编 Free考研考试 2021-12-29
  • 基于有效介质理论的物理性能计算模型的软件实现
    摘要:在改进的有效介质理论的基础上采用C++/Qt混合编程,设计并开发出一套复合材料物理性能模拟计算软件—CompositeStudio.该软件通过格林函数对本构方程进行求解,计算体积分数、颗粒长径比、取向分布、宏观位向对复合材料有效性能的影响.目前软件开发了弹性模量和介电常数两个模块,提供了友好的 ...
    本站小编 Free考研考试 2021-12-29
  • 纳米氧化锡负极材料锂化反应机理的原位透射电镜研究
    摘要:二氧化锡(SnO2)材料因具有储量丰富、理论容量高、嵌脱锂电位安全等一系列优点,在锂离子电池负极材料研究中受到广泛关注.然而,SnO2纳米材料在锂化反应过程中的机理,尤其是第一步转化反应是否可逆尚存在争议.本文利用常规水热法成功制备了平均粒径为4.4nm的SnO2纳米颗粒,并在透射电子显微镜中 ...
    本站小编 Free考研考试 2021-12-29
  • 等离激元材料和器件的动态调控研究进展
    摘要:等离激元学连接着光子学和电子学,在光产生、显微显示、数据存储、光集成和光子芯片、传感技术和纳米制造技术等方面展示出重要应用,正极大地促进既拥有纳米电子学的尺寸又兼有介电光子学速度的新一代信息材料和器件的发展.但是,传统上绝大部分等离激元材料和器件都是基于静态的设计,即一旦被制备,其性能也就确定 ...
    本站小编 Free考研考试 2021-12-29