永磁装置中磁场力的计算

1、工艺·技术·应用永磁装置中磁场力的计算王瑜(山西金山磁材有限公司山西太原 030024摘要:提供力学服务是永磁装置重要的用途之一,因而磁场力的计算是磁力机械设计、应用的重要内容。磁场力的计算有公式法和数值算法,本文结合两个具体实例对这两种算法进行了介绍,给出了计算结果,并与实测值进行比较。结果表明公式算法简单、方便,但计算误差较大;数值算法虽计算复杂,但精确、可靠。同时,文中对实例中的有限元数值算法提供了源程序,可供参考。关键词:永磁装置;磁力计算;公式算法;数值算法中图分类号:TM153+.3 文献标识码:B 文章编号:1001-3830(200705-0049-04Ca

2、lculation of Magnetic Force of Permanent Magnet DevicesWANG YuShanxi Jinshan Magnetic Material Co, Ltd, Taiyuan 030024, ChinaAbstract:Providing mechanic service is one of the important uses for the permanent magnet devices, thus the computation magnetic field strength is an important task for the de

3、sign and application of magnetic force machine.There are two methods for the magnetic field strength computation, ie. formula method and mumerical algorithm. This article unifies two concrete examples to carry on the introduction to these two algorithms, gives the computing results, and comparies th

4、e computing results with the measured values. The result indicates that the formula algorithm is simple, convenient, but the error is big; although the numercal algorithm is complex for calculating, its result is precise, reliable.At the same time, the article has provided all the FE source program

5、to the example that may supply reference to the readers.Key words: permanent magnet device; calculation of magnetic force; formula method; numerical algorithm1 引言以永磁材料为核心的永磁装置或用来产生磁场,或用来提供力学服务,所以磁场力的计算是永磁装置设计和使用中的一项重要内容。尤其是在磁选、吸合工具等行业常常要进行以永磁材料为磁场源的磁场力的计算,而这方面的相关资料较少。这主要是由于实用化的永磁材料尤其是第三代永磁材料钕铁硼出现相对较

6、晚,基于传统电磁理论的磁场力计算主要针对的是电磁力,即电磁铁磁场力的计算。同时用传统电磁理论进行磁场力计算时,如要保证精度则计算相当复杂或根本无法计算,如使用简化公式计算精度又不高,实用价值有限;新发展起来的数值算法精度高,但掌握数值分析的人员收稿日期:2006-07-10 修回日期:2007-03-31作者通信:E-mail: wy 对永磁行业又比较陌生,而永磁磁路设计人员大多又不熟悉数值计算,故涉及永磁装置磁场力计算的资料较少。经常有客户或相关行业的人员迫切需要这方面的知识,许多磁性材料网站论坛上也经常涉及这方面的问题,故笔者根据自己的一些经验及了解掌握的一些情况,对这方面内容作一较详细的

7、介绍。2 磁场力的公式算法2.1 磁场力计算公式的推导我们知道,载流导体和导磁材料在磁场中会受到力的作用,我们把这种力的作用称为磁场力。经典电磁理论认为,导磁材料在磁场中所受的力可归结为分子电流所受的力1。导磁材料在磁场中被磁化后,内部存在磁化电流,材料表面存在表面磁化电流,其磁化电流体密度和表面磁化电流面密度分别为v 、s : v =×M (1 s =-n ×M (2 式中M 为介质磁化强度,n 为表面法向矢量。则磁场对导磁材料的作用力为:v s vsvsd d (d (d (3V sV s =×+×=××+××

8、w w F B B M B n M B 式中B 为磁感应强度。对于各向同性介质有:V(d B v =××F M (4又01r rµµµ=M B ,0µ为真空磁导率,r µ为磁介质相对磁导率,经矢量运算可得2r 0rv1d 2v µµµ=F B (5由矢量梯度积分公式vd d v =s 可得2r 0rs1d 2µµµ=F B s (6该式即为磁场力计算公式。由于磁场在导磁材料所在区域分布的复杂性,故直接使用该式积分来求磁场力往往比较困难,实际使用中,我们一般假设磁场在

9、导磁材料所在区域分布均匀一致,又由于导磁材料磁导率较大(即r µ>>1,这样根据磁路设计和计算实际我们把求解磁场力大小的公式简化表示为22r 0r 0111222F B S B S BHS µµµµ= (7式中B 为磁场与导磁材料作用面处的磁感应强度,H 为磁场与导磁材料作用面处的磁场强度,S 为磁场与导磁材料作用面的面积。公式采用SI 单位制,即式中F 、B 、H 、S 单位分别为N 、T 、A/m 、m 2。实际工程应用中,经常使用其它单位制,为方便使用,上式在其它常用单位中的表达形式表述如下100.397898F BHS BH

10、S = (8 使用该式F 、B 、H 、S 单位分别为N 、kG 、kOe 、cm 2;5108F BHS = (9 使用该式F 、B 、H 、S 单位分别为N 、G 、Oe 、cm 2;62101(85000F BHS BHS = (10使用该式F 、B 、H 、S 单位分别为kgf 、G 、Oe 、cm 2。在磁场力的求解时,也经常根据具体情况使用虚位移法来求导磁材料在磁场中的受力:即有n 个回路构成的系统,如仅有一个广义坐标x 发生变化,则首先写出磁回路系统能量表达式W m ,则在各回路磁通保持不变的情况下,导磁材料在磁场中的受力为mW F x= (112.2 计算实例如图1所示,规格为

11、50×50×25mm ,牌号为N35 的NdFeB 磁体(B r =1.25T ,磁体和衔铁之间的距离为6mm 。我们来计算一下磁体对衔铁的吸引力。从公式可知,首先需计算衔铁处的B ,由于磁体与衔铁所构成的磁路漏磁太大,计算误差过大,不适合用磁路建立方程求解,故我们简化处理,直接用距孤立磁体表面中心点x 处的磁场强度公式来求B ,结果相对更准确些(工程实践中直接用高斯计测量也可。 r tan tan (12B B =式中L 、W 、H 分别为磁体长、宽和厚度。据公式可得离磁体6mm 处磁感应强度B 的值约为0.42T, 则201F B S µ=176N如图2所示,

12、磁体、软磁回路和衔铁构成的简 图1 磁体对衔铁的吸引示意图单磁路,磁体规格为30×30×30mm ,牌号为N35 的NdFeB 磁体(B r 1.25T ,磁体和衔铁之间的距离为2mm 。我们来计算一下磁体对衔铁的吸引力。首先说明一点的是,磁体矫顽力H cb 对磁路及气隙磁场也会有影响,这里我们假定磁体矫顽力H cb 足够高,足以保证磁体工作点(B m , H m 始终能处于磁体B r -H cb 退磁曲线弯曲点以上,则计算时我们就可以不再考虑磁体矫顽力的影响。我们来计算衔铁处的B ,对于磁路,我们可建立下面的方程组m m f g g m m f g g m r 0r m

13、(13B S K B S H L K H L B B H µµ×=×××=××=式中K f 、K r 分别为漏磁系数和磁阻系数,K f 取值为2.5,K r 取值为1.2;S m =9cm 2、L m =3cm 分别为磁体截面、长度;S g =9cm 2、L g =0.2cm 分别为气隙截面和长度;B m 、H m 为磁体的工作点;r µ=1.05为磁体相对磁导率。上述三个公式联立可求出磁体B g =0.48T ,每个磁极磁场力为2g g 012F B S µ=,则总磁场力为2g g 01216

14、6N 2F B S µ=×=使用公式法时, K f 、K r 的确定更多靠的是经验,如没有实际磁路设计经验, 要准确确定K f 、K 比较困难。3 磁场力的数值计算3.1 磁场力的数值计算方法介绍公式法计算磁场力简单方便,但相关参数难以准确估计,误差较大,因而大大限制了它的实用价值。磁场力的精确计算需应用数值分析方法,如有限元法、有限差分法等,其中发展较成熟的是有限元法。有限元法是基于建立起来的数字模型,用现代数学方法求解有关微分方程定解问题,并对求解结果进行处理和解释的一种数学方法。可供磁场有限元分析的计算机软件很多,有专用软件,有商业化通用软件。可进行电磁场分析的比较有

15、影响的商业化有限元软件有:美国SwansonAnalysis System 的Ansys ;美国MacNeal Schwendlert Corp. 的MSC/EMS 软件;美国Ansoft Corp.的Maxwell 软件等,其中使用Ansys 的客户最多,分析领域也很广泛,下面将以Ansys 为例进行介绍。Ansys 有限元软件处理问题的一般步骤为:建立几何模型;将问题中涉及到的材料参数进行定义;根据不同的模型要求选定所需的单元类型;给各个几何模型模块赋予材料属性和单元属性;对几何模型进行网格划分;根据物理问题的不同给问题施加载荷(如电流密度并指定边界条件(定解条件;选用合适的求解器和求解参

16、数对问题进行求解;运用后处理程序按需要查看不同的求解结果。具体学习时可参考这方面的相关书籍24。3.2 计算实例下面我们仍以上面提到的例子,对如何应用用Ansys 软件处理问题,有图形模式和命令流模式两种方法,为方便介绍,我们采用命令流模式。对第一个例子,用Ansys 软件求解该问题的命令流程序如下(Ansys 有限元程序语句中,以!、/com 、C*开头的内容为注释性文字,具体求解时其内容可以不必输入。/COM, Magnetic-Nodal !选用基于节点的磁场分析环境/PREP7 !进入前处理器BLOCK,-0.025,0.025,-0.025,0.025,-0.028,-0.003 !

17、建立磁体几何模型BLOCK,-0.025,0.025,-0.025,0.025,0.003,0.028 !建立衔铁几何模型BLOCK,-0.05,0.05,-0.05,0.05,-0.05,0.05 !建立外部区域几何模型VOVLAP,all !对几何模型进行布尔操作 ET,1,SOLID98,10 !单元选择MP,MURX,1,1.05 !定义磁体(1#材料磁导率 MP,MGXX,1,0 !定义磁体X 轴向H cb MP,MGYY ,1,0 !定义磁体Y 轴向H cb MP,MGZZ,1,950000 !定义磁体Z 轴向H cb(钕图2 磁体回路对衔铁的吸引示意图铁硼材料退磁曲线按线性处理:

18、即B r =0r µµH cb 。 MP,MURX,2,2000 !定义2#材料衔铁(电工纯铁磁导率。定义磁导率严格来说要定义其B -H 曲线,这里直接取一常量,从简。MP,MURX,3,1 !定义空气(3#材料磁导率 MSHAPE,1,3D !划分单元参数选取 MSHKEY ,0 ! 设置划分单元网格大小 type,1 !选取单元类型 mat,1 !选取1#材料 vmesh,1 !对1#体积单元划分 mat,2 ! 选取2#材料vmesh,2 ! 对2#体积单元划分 mat,3 ! 选取3#材料 vmesh,4 ! 对4#体积单元划分ESEL,S,MAT,2 !选取材料号

19、为2的单元 CM,xiantie,ELEM !定义为衔铁组元 FINISH !退出前处理器 /SOL !进入求解器FMAGBC,'XIANTIE' !加载衔铁磁场力标志 ALLSEL,ALL ! 选中所有单元MAGSOLV ,2,0.001,25,0 !选用静态求解器求解FINISH !退出求解器 /POST1 ! 进入后处理器 /VSCALE,1,1,0 ! 单位制PLVECT,B,VECT,ELEM,ON,0 !显示磁感应强度B 分布FMAGSUM,'XIANTIE' !磁场力计算结果显示上述有限元程序计算表明,用虚功原理计算磁吸合力为119.8N ,用Ma

20、xwell 应力法计算磁吸合力为120.4N ,图3为磁场力计算结果显示情况。对于第二个例子,用有限元数值计算结果为:用虚功原理计算磁吸合力为156.5N ,用Maxwell 应力法计算磁吸合力为155.7N 。图4为磁场力计算结果显示情况。同时仍将源程序附后,便于参考。 /COM, Magnetic-edge /PREP7 SMRT,ONBLOCK,0,0.03,0,0.06,-0.015,0.015 BLOCK,0,0.03,0.06,0.09,-0.015,0.015 BLOCK,0,0.03,0.09,0.15,-0.015,0.015 BLOCK,0.03,0.06,0,0.03,-

21、0.015,0.015BLOCK,0.03,0.06,0.12,0.15,-0.015,0.015 BLOCK,0.062,0.092,0,0.15,-0.015,0.015 VGLUE,ALLBLOCK,-0.03,0.122,-0.03,0.18,-0.045,0.045 VOVLAP,ALL ET,1,SOLID117 MP,MURX,1,1.05 MP,MGXX,1,0MP,MGYY ,1,950000 MP,MGZZ,1,0MP,MURX,2,200 !A3钢的磁导率约为200 MP,MURX,3,200 MP,MURX,4,1 MSHAPE,1,3D MSHKEY ,0 TYPE,

22、1 MAT,1 VMESH,7 MAT,2VMESH,8,11,1 MAT,3 VMESH,6 MAT,4 VMESH,2(下转60页图3 例1磁场力有限元计算结果图4 例2磁场力有限元计算结果但大量的商用电子设备正朝着小型化和高频化方向发展,目前使用频率6GHz以下(甚至低至30kHz的电磁波为多,因此有必要对适用于这一频段的抗电磁干扰铁氧体吸波材料进行特别的研究,以达到在轻型薄层化、宽频带的同时做到高效吸收。粘结类吸波制品特别是柔性橡胶类材料具有密度小、厚度薄、工艺简单无污染、能制作复杂形状、易于回收等优点,正被越来越多的研究人员所关注,预计将有更为广阔的实用前景和发展空间。参考文献:1

23、TDK产品目录. 2006. 2 邓龙江, 等. J. 功能材料, 2001, 32(2: 144.3 FDK产品目录. 2006.4 过璧君, 邓龙江. J. 电子科技大学学报, 1992, 21(2:158.5 过璧君, 等. J. 电子科技大学学报, 1992, 21(3: 312.6 刘素琴, 等. J. 磁性材料及器件, 2000, 31(2: 12.7 Ruan SP, et al. J. J Magn Magn Mater, 2000, 212: 175.8 张雄, 等. J. 材料导报, 2003, 17: 69.9 孟凡君, 等. J. 无机化学学报, 2002, 18(10

24、: 1067.10 张海军, 等. J. 功能材料, 2003, 34(1: 39.作者简介:陶振声(1978-,男,吉林省吉林市人,从事磁性材料及器件的研究与制造。(上接52页ESEL,S,MAT,3CM,XIANTIE,ELEMFINISH/SOLFMAGBC,'XIANTIE'ALLSEL,ALLANTYPE,STATICMAGSOLVFINISH/POST1/VSCALE,1,1,0PLVECT,B,VECT,ELEM,ON,0FMAGSUM,'XIANTIE'3.3 实验验证及结果说明为了验证上述求解结果,笔者对计算结果进行了实际验证。具体方法是:按上述实例中磁体、软磁轭及衔铁部分尺寸加工模型,充磁组装,磁块部分与衔铁之间按气