It is of great significance to accurately evaluate the magnetization force and Lorentz force etc.electromagnetic loads of the vacuum vessel, divertor and blanket of tokamak fusion reactors to optimize the design and to ensure the structural safety. A numerical method for magnetization force calculation is introduced utilizing the Maxwell stress tensor approach based on the edge element provided by ANSYS. This method can calculate not only the total magnetization force, but also can obtain the equivalent distribution of the loads on the studying object. The results show that the force calculated by this method can be affected by the size and deformation of the elements. Fortunately, if the force is output only through the corner nodes, the influence is negligible and the total force agrees well with the analytical one. Moreover, the qualitative analysis shows that the result of magnetization force distribution is consistent with the relevant theory of electromagnetism. This numerical method can be applied to the magnetization force analysis and subsequent electromagnetic-structural coupling analysis of the blankets in fusion reactors.

Tokamak fusion reactor, Magnetization force, ANSYS, Edge element, Maxwell stress tensor approach

为了保证托卡马克聚变堆的真空室、偏滤器、包层等结构或部件的安全性，需要在设计过程中准确评估其磁化力和洛伦兹力等电磁力载荷。基于 ANSYS 提供的棱边元，介绍了一种采用麦克斯韦应力张量法的磁化力数值计算方法。该方法不仅可以计算总磁化力，还可以获得载荷在研究对象上的等效分布。研究表明，该方法的计算结果会受到单元尺寸和变形程度的影响，但当仅通过角节点输出力的计算结果时，该影响较小，且计算结果非常接近解析法的近似计算结果。此外，定性分析表明计算的磁化力分布与电磁学相关理论一致。该数值计算方法可应用于针对聚变堆包层的磁化力分析及后续的电磁-结构耦合分析。

托卡马克聚变堆, 磁化力, ANSYS, 棱边元, 麦克斯韦应力张量法