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 Name
ModelName = 'Woofer.fem';

-- Maximum excursion + /- from the centered position:
Xlim = 10;

-- Movement increments used during the analysis
dX = 0.1;

--------------------------------------
-- Analyze BL and incremental inductance at 1 mm steps between 
- Xlim and + Xlim
open(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);
end

mi_close();
remove('temp.fem');
remove('temp.ans');

print('');
print('DC coil resistance = ',R);

更多技術內容，關注：揚聲器系統設計與仿真