纳米晶永磁Pr-%2c2-Fe-%2c14-B的微磁学有限元法的模拟计算的研究.pdfVIP

苎±=星全里壁兰塑壁丝苎壁垒望垒奎叁———— 纳米晶永磁Pr：Fe，．8的微磁学有限元法的模拟计算研究 张宏伟 荣传兵 张健 张绍英 沈保根 中国科学院物理研究所凝聚态物理中心，磁学国家重点实验室，北京100080 摘 要根据实验数据，我们构造了接近实际纳米永磁Pr。Fe。．B的样品，用傲磁学有限元法进行了 模拟计算．计算结果表明，晶界处各向异性的下降会导致矫顽力降低、剩磁增加，而晶界处交换作用常数 的减小则会使朝硅降低、矫顽力增加；通过对实验样品的模拟研究发现，晶界处备向异性和交换作用常数 的共同下降能够同时拟合出真实的矫顽力和剩磁．模拟计算与实验在退磁曲线形状上的差距则说明，我们 曲横拟还存在不足． 1引言 自从最大磁能积高达1MJ／m3的纳米复合永磁材料被预言以来[1】，纳米晶永磁研究成为近十年磁 MJ／m3【21。理论预言与实 学研究的重点之一㈤。但实验上各向同性磁体的最大磁能积尚未能超过0．2 验差距巨大的主要原因体现在：(1)实验上不能够同时保证较高的矫顽力和剩磁；(2)理论预言最大 磁能积1MJim3的模型过于简单。因此，人们对微磁学有限元法的模拟计算研究寄以希望。 永磁材料的硬磁性能既取决于材料的内禀特性又决定于材料的微结构。关于三维磁体磁滞回线的 学有限元法。微磁学有限元法己在磁体设计方面取得了很大的成功。近几年，这一方法又应用于纳米 晶永磁的磁滞回线计算研究，并且取得了一定的进展i1-tl】。如果微磁学有限元法能够模拟实际样品的 磁滞回线，那么这一算法对提高磁体性能的指导价值是显而易见的。然而，迄今为止没有发现用这种 算法模拟实验上磁滞回线的研究报道。因此，本文将报道我们采用微磁学有限元法来模拟真实退磁曲 线的尝试。 2．实验和计算方法 合金快淬条带样品是各向同性的。HRTEM测试结果显示，快淬样品的晶粒约为20nm，晶粒边界处 的缺陷层厚度约3脚f12】。磁性测量是用SQUID在室温下进行的。 模拟计算的样品是由64个多面体晶粒构成的立方体，在其相邻晶粒间构造了3nm晶间层$】。计 算退磁曲线的有限元算法(以下简称FEM)参见文献17,11。我们使用了40000个有限元对样品进行了 剖分：模拟所采用的磁饱和极化强度Js为1．6T、各向异性常数K为5．567MJ／m3、交换积分常数A为 x 7．7 10-12J／m。 3结果与讨论 计算的曲线为理想晶界条件。尽管计算的退磁曲线与实验的相差较大，但FEM得到的结果要优于 粒尺寸的较宽分布使矫顽力附近的反磁化分步进行，这可能是源于大尺寸晶粒与小尺寸晶粒之间的磁 化反转场的不周。 ． 对于理想晶界条件。FEM计算的矫顽力和剩磁均大于实验结果。鉴于实际样品中晶粒边界处存 在缺陷的事实，我们对晶界处的参数K和A做了调整．图2给出了晶界处参数K岫”对剩磁和矫顽力 的影响。由图可以看出。晶界处K蛔一的下降使矫顽力降低。并在K岍r砷．15K时计算值与实验值相等。 但是．计算的剩磁随K-一的下降而增加，这使得剩磁的计算值与实验值差距加大。 第十一届全国磋学和磁性材料会议论文篥 图3给出了晶界处参数Al-y。r对剩磁和矫顽力的影响。可以看出，随着A1”的下降，剩磁是减少 的，而矫顽力是增大的。这一趋势与文献【8·1q的结果相近。但在数值变化速率上均小于文献tto]的值。 我们认为，这主要是在样品构造和有限元形状方面的不同造成的。与文献‘8】的计算结果相比，在剩磁 与A]．YcVA近似呈线性关系上是一致的，但我们计算的矫顽力却与A1。Y“／A不呈线性。我们认为矫顽力 近，但矫顽力却随着A1ayer的下降与实验值相差变大。 由图2和图3的结果得出，单独调整K1一或A1一的数值是很难同时拟合实验上的矫顽力和剩磁。 我们进一步同时调接了K1ayes和A1ayes的数值，图4分别给出了晶粒大小分布较窄的样品的典型计算结 果。对于适当的K“”r和A1一数值，我们可以拟合出接近实验的矫顽力和剩磁值，但是在退磁曲线的 形状上有较大的差距。这说明我们在模拟实验曲线上还存在关键问题尚未找到。进一步的研究仍在进 行中。