COMSOL磁体受力计算问题

这是我研究生论文重要计算的一部分，值此临近毕业之际，讨论一下该问题。先说结论，COMSOL默认的磁体受力计算方法为FEM即有限元法，该方法磁体受力计算结果会出现在附近没有吸引物（铁磁性物体）时受力不为零的现象，且不为零的值较大；而使用BEM即边界元法，在没有吸引物时磁体受力计算结果接近于零；以我毕业论文为例，对这个问题进行了阐述。

 从图中可以明显看到就计算结果而言，边界元法相较于有限元法，计算结果更具有可靠性和规律性。图3-4和图3-5计算结果曲线值相差不大，是因为在作图3-4时我将没有吸引物时磁体受力结果加回去了。因此，就对结果数据处理来说，FEM计算结果还需要进一步处理，而BEM计算结果直接使用也可以保证计算精度。

其实从我对两种方法的浅薄研究，我认为两种方法其实都是BEM，只是有限元法选取的边界就是你选择磁体的表面，而边界元法表面是你重新绘制的表面而已。正如图中我的解释一样，存在一种类似应力集中的现象，越靠近磁体原来表面其存在越严重的应力集中现象，而磁体受力是积分的结果，当应力集中越严重其积分过程中产生的数据误差也越大。COMSOL内置的磁体受力积分方法公式如下

 式中：P表示应力张量、r为你选取的力臂、S为积分表面。

具体实施操作方法可以参考COMSOL官网给出的参考案例，如下：

力计算 3 - 磁矩 BEM FEM

当然在使用边界元法的时候我也发现，虽然边界元法其计算结果优于有限元法，但是其计算精度严重依赖于你建立的封闭曲面形状和大小，需要自己不断尝试才能获得高精度。

国内对于磁体受力研究的论文较少，从理论去深究边界元法计算磁体受力的论文更少，想要深入研究可以看看国外研究文章。