使用Femm进行轴对称磁路非线性仿真
Femm是一个免费的有限元软件，可以进行2维平面或2维轴对称的电、磁、热的有限元仿真分析。通常会用femm做磁路的仿真，优化磁路设计，或者计算Bl值。实际上femm的功能还有不少拓展的空间，可以求解Bl(X)、Le(x)等等磁路的非线性，短路环对Le（x）的影响等等。其自带了Lua脚本语言的输入窗口和编译器，可以直接执行Lua命令。或者也可以耦合Matlab、Mathematica、Octave（类似Matlab的免费软件）Femm软件的帮助文件中有详细讲解如何进行接口参数调用。Femm官网也有对应的案例参考。Femm官网网址，可能需要翻墙http://www.femm.info/wiki/HomePage附录一款产品求解磁路Bl（x）和Le（x）的Lua命令，注释已经写得很详细清晰：-- Model NameModelName = 'Woofer.fem';-- Maximum excursion + /- from the centered position:Xlim = 10;-- Movement increments used during the analysisdX = 0.1;---------------------------------------- Analyze BL and incremental inductance at 1 mm steps between - Xlim and + Xlimopen(ModelName);mi_saveas('temp.fem');mi_selectgroup(1);mi_movetranslate(0, -Xlim);mi_clearselected();showconsole();clearconsole();print('Disp(mm)','BL(N/A)','Inductance(uH)');for k=-Xlim,Xlim,dX do mi_modifycircprop('icoil',1,1); mi_analyze(1); mi_loadsolution(); mo_groupselectblock(1); fz = mo_blockintegral(12); parm1,R,fl1 = mo_getcircuitproperties('icoil'); mo_close(); mi_modifycircprop('icoil', 1, 0); mi_analyze(1); mi_loadsolution(); parm1,parm2,fl0 = mo_getcircuitproperties('icoil'); L = (fl1 - fl0)*10^6; print(k, fz, L); mi_selectgroup(1); mi_movetranslate(0, 1);endmi_close();remove('temp.fem');remove('temp.ans');print('');print('DC coil resistance = ',R);来源：声学号角
