永磁轴承计算

1、汇报人：XXX2022-3-20目录1、磁学基本概念2、永磁轴承分类3、磁场的计算4、轴承磁力计算磁学基础知识 静磁现象 材料的磁化 磁性与磁性材料的分类磁化曲线磁滞回线物质的磁性分类磁性材料分类磁矩m磁化强度M磁场强度H和磁感应强度B磁化率和磁导率退磁场静磁能电场电场磁场磁场从从N N极出发，进入与其最邻近的极出发，进入与其最邻近的S S极，并形成闭合回路；极，并形成闭合回路；通常呈直线或曲线，不存在呈直通常呈直线或曲线，不存在呈直角拐弯的磁力线；角拐弯的磁力线；任意二条同向磁力线之间相互排任意二条同向磁力线之间相互排斥，因此不存在相交的磁力线斥，因此不存在相交的磁力线；电荷之间存在相互作用

2、221rqqkF 电荷周围存在电场，可以用电力线来表示磁力线特点：如果一个小磁体能够用无限小的电流回路如果一个小磁体能够用无限小的电流回路来表示，我们就称为来表示，我们就称为磁偶极子磁偶极子。用磁偶极。用磁偶极矩矩j jm m表示：表示：磁偶极子和磁矩与磁偶极子等效的平面回路的电流和回路与磁偶极子等效的平面回路的电流和回路面积的乘积定义为面积的乘积定义为磁矩磁矩，用，用m表示：表示：jm=mlm=iS i+m-mliSm 磁偶极矩和磁矩具有相同的物理意义，存在关系：磁偶极矩和磁矩具有相同的物理意义，存在关系：jm=0m ，o=4=41010-7-7H Hm m-1-1 ，真空磁导率，真空磁导率

3、m为磁极强度磁极化强度磁极化强度J Jm m、磁化强度M定义单位体积磁体内磁偶极子具有的磁偶极矩矢量和称为定义单位体积磁体内磁偶极子具有的磁偶极矩矢量和称为磁极化强度磁极化强度，用，用J Jm m 表示；表示； 单位体积磁体内磁偶极子具有的磁矩矢量和称为单位体积磁体内磁偶极子具有的磁矩矢量和称为磁化强度磁化强度，用用M M表示。表示。J J m和和M M亦有如下关系：亦有如下关系： Jm=0M VjJmm2/mWbVMmmA /磁化强度可以看成是磁偶极子的集合磁化强度又可以看成是闭合电流环的集合等效磁荷模型等效电流模型磁化强度M和磁场强度H存在如下关系： Mc H 或 cM/Hc称为磁体的磁化

4、率，表征磁体磁性强弱的参量。由此 B=0(H+c H)0(1+c)H 定义（1c）为相对磁导率，即B/ 0 H。 磁导率是表征磁体的磁性、导磁性及磁化难易程度的磁学量。磁化率和磁导率磁化曲线在非铁磁材料中，磁通密度 B 与磁场强度 H 成比即B = 0H 式中，0 真空磁导率。0 = 410-7H /m，B 与H 呈线性关系。铁磁材料的磁通密度（即磁感应强度）B 与磁场强度 H 呈非线性关系，即 即 B =f (H)是一条曲线，称磁化曲线。 将铁磁材料磁化一个循环，得到一个闭合回线 abcdefa，称为铁磁材料的磁滞回线。磁滞回线Br，Mr表示剩磁M=0时的矫顽力，称为内禀矫顽力B=0时的矫顽

5、力，称为磁通矫顽力矫顽力是表征材料在磁化以后保持磁化状态的能力通常将BHC81038105的材料称为硬磁材料硬磁材料；介于120kA m-1之间的为半硬磁材料退磁曲线退磁曲线上每一点的B和H的乘积(BH)为磁能积，表征永磁材料中能量大小的物理量。(BH)的最大值为最大磁能积(BH) max按磁化率的大小，可将物质磁性分为五个种类：抗磁性 ：原子无磁矩，磁化率 ，且顺磁性：磁矩混乱，当外加磁场时，呈现弱磁性，磁化强度与外磁场相同， 。反铁磁性：当负交换时，若磁矩相互抵消，为反铁磁；若外加磁场时，显弱磁性， 。亚铁磁性：当负交换时，若磁矩不相互抵消，显一定的磁性，为亚铁磁材料；易被磁化， 。铁磁性

6、：当正交换时，磁矩平行排列，显强磁性；极易被磁化 ， 。物质的磁性分类0dc0c581010在c0c510c非铁磁性材料非铁磁性材料0c0c511010在c铁磁材料铁磁材料HexeWABWAAWBBABWiiWiiWijBA（a）抗磁性（b）顺磁性（c）反铁磁性（d）铁磁性（e）亚铁磁性五种磁性材料的基本磁结构磁学中单位统一磁学中单位统一径向永磁轴承：1、径向充磁2、轴承充磁径向永磁轴承的分类3、天津大学张大卫新型磁气混合永磁轴承永磁轴承中重要理论1、根据Earnshaw理论，仅由永磁材料构成的磁轴承，不可能实现稳定的完全悬浮，至少在1个自由度上是不稳定的。径向永磁轴承，径向上稳定，而轴向上不

7、稳定；轴向上稳定，径向上不稳定。2、轴向上固定，有利于提高径向承载力和刚度。3、永磁轴承的研究与设计过程中，承载能力和刚度的计算是非常重要的。永磁体的磁场计算方法：磁化磁体的物理计算模型磁化磁体磁场的计算空间磁场计算的数值方法等效磁荷模型等效电流模型有限差分法有限元法边界元法Mm0为真空磁导率0m等效磁荷模型认为，磁体内部存在密度为 的磁荷分布，并有如下关系:磁化强度M 永磁体沿磁化方向均匀充磁，M为常矢量， 恒等于零。但是在边界上由于M并不连续，因而在磁体边界上仍然存在面磁荷密度msmMns0m位矢量磁体边界的外法向单n永磁体空间磁场的分析计算及其在永磁磁力轴承中的应用-汤双清用等效磁荷法计

8、算永磁体磁场-李景天根据 Maxwell 方程组中的高斯定律, 可得由永磁体产生的空间任意点上标量磁位为:dSRdVRRSsVm+m00m4141可以求出永磁体周围的空间磁场强度:+dSRRdVRRHSsVmm3m041dSRRHSsm3m041根据标量磁位 与磁场强度的关系:mmHm因为 恒等于零，所以磁场强度可表示为：m同轴环形磁铁磁作用力计算的等效磁荷法修世超度永磁环剩余磁感应强真空磁导率；平面上投影；在两点间的距离；内外环端面；磁荷端面磁荷面密度外环内外半径；内环内外半径；内外磁环的轴向高度位移；内环相对外环的轴向方向一致；径向外载荷，假定与方向的位移；内环相对外环沿径向上任一点的极坐

9、标；端面上任一点的极坐标；端面rmsByzrrPQrRRRRHxyFerr014321032,3,2,永磁向心轴承承载能力与刚度的计算谭庆昌Permanent Magnet Bearings and CouplingsYonnet内环端面2上任一点P处的磁荷量为:222drdrdSqmsms根据电磁学理论电磁学理论,P 点处的点磁荷在外环端面 3上任一点Q处产生的磁场强度为：0222002204141rrdrdrrrqHmsQ点处的磁荷量为333drdrdSqmsms023232023334rrdrdrddrrdrdrHqHFdmsmsQ点处的点磁荷所受到P点产生磁场的的磁场力为：由于rrr

10、0则rrdrdrddrrdrdrHqHFdmsms43323202333根据静磁学理论公式：nMnJs0m()MHB+0磁极化强度矢量J由于磁化方向与磁环端面法线重合，因此综合上面公式可得出：当时，0HrBB所以MB0rsBM 0m位矢量磁体边界的外法向单n024rBC 其中C 是与永磁材料的磁性参数有关的系数, 称为磁性系数因此， 可表示成rrdrdrddrrBFdr43323202rrdrdrddrrC33232Fd向y方向投影，可得到PQ两点处点磁荷沿径向外载荷方向的作用力jrrdrdrddrrCdF33232其中根据几何关系可得：()()22322320sinsincoscosrrer

11、rxr+对整个端面积分,得端面2与3的磁作用力沿径向外载方向上的大小：()()()2322322320323223202032sinsincoscoscoscos4321rrerrxdrdrddrrerrCFRRRR+此磁作用力的数学表达式包含了十分复杂的积分。用MATLAB, Maple和Mathematic。这三大数学软件进行数值积分，结果不仅耗费大量时间，而且很难得出计算结果。物理学上常用蒙特卡洛方法来计算高维积分，但这种方法积出的结果是随机的，只有通过统计分析后才能找到大致的结果。模型简化()byxP,11P点坐标()ebyxQ+,22Q点坐标112dydxqmsP点处的点磁荷量223

12、dydxqmsQ点处的点磁荷量rrdydxdydxqHFdms432211023rrdydxrrqHms41413110320是P指向Q的矢量r()()212212212eyyxxPQr+r的大小是P点处的点磁荷在Q点处产生的磁场强度矢量：P点产生的磁场对Q点的磁力：()322110232211024e-eb4rdydxdydxerdydxdydxdFmsms+()()+a02322122122211l0a0l002234eyyxxdydxdydxeFms()()+a02322122122211l0a0l023eyyxxdydxdydxeCF所以024rBC 其中rsBM 0m令由于对称性，沿

13、X，沿Y方向的力相互抵消。只求沿Z轴方向的力对整个端面积分,得端面2与3的磁作用力沿Z轴方向上的力：1xN1yNlajiP ,111 yxS111yxNaNlS 1111dydxdSS()()+a02322122122211l0a0l0eyyxxdydxdydx直接积四重积分，很难解出真实解，从原式物理意义出发，利用有限元思想，给出数值计算方法。把2面沿X轴方向分成 份，沿Y轴方向分成 份1xN1yNjiP ,表示为 +2111iNlixx +2111jNajyy()11, yxP同理，面3可以看成：222yxNaNlS 222 yxS2222dydxdSS121NNNxx221NNNyy令(

14、)kiNliNlkNlxxxx+112212112()jhNljNahNayyyy+212212112()()+a02322122122211l0a0l0eyyxxdydxdydxhkQ , +2122kNlkxx +2122hNajyy可表示为()11, yxQ()()()()23222122111111111a02322122122211l0a0l02121+jhNlkiNlNaNleyyxxdydxdydxNhNkNjNi通过分割把原积分化成如下算式：通过分割把原积分化成如下算式：()()232221221111111112121+jhNlkiNlNaNlfNhNkNjNi令()()feCeyyxxdydxdydxeCFa02322122122211l0a0l023+则永磁体规格：厚度b=5mm;宽度a=10mm;长度l=20mm;