电工钢片磁特性模型研究的几个问题.docxVIP

电工钢片磁特性模型研究的几个问题 计算电磁学是一门综合着磁力学理论、数值计算方法和计算机软件技术等学科的新兴学科。在过去30年中，计算机技术的快速进步和普及，计算数学和恶意软件的快速发展，以及电气设备在工业中的应用，为该学科的快速发展带来了很大的推动。对磁性材料在磁值分析中的精确模拟，与电气产品的能耗和优化设计密切相关。这是当前计算磁体学的热点问题之一。在磁性材料中，电磁工具片（也称为硅钢片）被广泛用于各种电机和设备的磁极材料，形成主磁路。因此，其磁特性的模拟受到了很多关注。研究问题包括两个方面。一是创造合适的材料特征数学模型。其次，采用研究与材料特征模型相结合的干磁值计算方法。研究的困难可以概括如下。 1)铁磁材料的磁特性具有特殊的非线性性质,包括与磁化过程相关的多值性(即包含极限磁滞回线和局部磁滞回线在内的磁滞特性)和各向异性. 2)若要在工程电磁场数值分析中将上述材料特性尽可能精确地结合进去,需要冗长的计算时间. 3)建立任何一种磁特性模型,都需要对表征材料特性的若干组参数进行实际测量,通常用于测量的实验装置和测量方法是非标准的,并且不容易操作. 由于问题本身的复杂性,目前这一领域的研究仍然处于发展阶段.本文将以电工钢片材料为例,介绍磁特性模型研究的新进展,并针对上述关键和难点问题展开讨论. 1 真空磁导率检测 磁性材料的宏观磁性质可用下式表征 式中:μ0——真空磁导率; B、H和M——磁通密度、磁场强度和磁化强度. 当磁场变化时,材料呈现复杂的磁滞特性. 1.1 磁滞特性的模拟 在传统的电磁场数值计算中,通常并未考虑材料的磁滞特性,磁通密度B与磁场强度H之间用一个单值函数来表示,在不同的激励条件下该单值函数可采用平均磁化曲线、交流磁化曲线或有效磁导率.即使考虑磁滞特性,通常也只利用极限磁滞回线而不考虑局部回线.但是在需要高精度计算的场合,这样的简化处理就不能满足要求.磁滞特性的模拟本是一个经典问题,描述磁滞特性的Mandelung定则于1905年就已提出.与图1相对照,该定则可表述为: 1)磁滞回线上的任一回转(reversal)曲线由反向点坐标唯一确定. 2)如果曲线3-4-1上的任意一点4成为新的反向点,那么从点4出发的曲线4-5-3将回到起始点3,这一特性称为“返回点记忆”(return-pointmemory). 3)如果曲线4-5-3上的点5成为新的返回点,并且曲线5-4超过了点4,那么工作点将沿曲线3-4-1的4-1段移动,就像是封闭的回环4-5-4根本不存在一样,这称为“擦除特性”(wiping-outproperty). 上述定则描述了宏观磁滞现象的几个基本特性,但并未包括复杂工况下电机电器中磁性材料磁滞特性的一些其它重要特征. 1.2 回转时磁场强度的阈值 Preisach类磁滞模型是目前采用最多的一种磁滞模型.德国物理学家F.Preisach于1935年基于磁化机理提出了标量Preisach模型,以后M.Krasnoselskii将这一模型从它的物理思想中分离出来,用一种纯数学形式来描述,在此基础上形成了经典Preisach磁滞模型.其基本原理是将磁性物质表示为一组磁偶极子的集合,这些磁偶极子都具有如图2a所示的矩形磁滞特性,材料的局部宏观磁特性被看作是这些磁偶极子磁滞特性的总和,磁通密度B与磁场强度H的关系可以表示为 式中:S——图2b中的三角形区域; α和β——磁偶极子极性反转时磁场强度的阈值; μ(α,β)——磁偶极子的分布函数,且μ(α,β)=μ(-α,-β),当(α,β)∈/S时,μ(α,β)=0,并且在S+区域内有γαβ=1,在S-区域内有γαβ=-1;Hsat为饱和磁场强度,且Hsat≥α≥β≥-Hsat,在S+和S-之间的阶梯形分界线的每一个顶点对应于过去外磁场变化的极值. 应用Preisach模型的关键是确定其分布函数μ(α,β),这需要大量的实验测试和统计工作.许多研究者提出了不同的方法,其中Everett函数是广为人知的一种,可以应用Everett函数的二阶导数与测试值相结合来计算μ(α,β), 式中,Ev代表Everett函数,求导运算需要采用数值方法.在许多确定分布函数的其它方法中也包含数值微分和积分运算,这将引入额外的数值误差并增加计算的复杂性.文献提出了一种方法,在确定Everett函数时需要测定一系列局部磁滞回线,每个局部回线要求精确地对中(输入磁场的正、负绝对值相等),这一方法的优点是在将Everett函数与Preisach模型相结合计算磁化强度时不需要进行数值微分和积分.为了减少实验工作量,有些研究者提出,可以按照一定的函数变换,由B-H或H-M极限磁滞回线上的数据计算出μ(α,β),这种变换只须利用极限磁滞回线数据和工作点变化的历史,就可表述包含局部磁滞回线在