function fe=equload
%--- 計算等效結點力
global  element  nfe  nfe_number nf  node_number %（定義全局變量）
fe = zeros(node_number*3,1) ;%（創建fe）
for i=1:1:nfe_number
    ee=nfe(i,1);
    ni=element(ee,1);
    nj=element(ee,2);
    t=trans(ee);
    FF=efix(i);
    wf=t*transpose(FF);
    fe((ni-1)*3 + 1 ) = fe((ni-1)*3 + 1 )-wf(1,1);
    fe((ni-1)*3 + 2 ) = fe((ni-1)*3 + 2 )-wf(2,1);
    fe((ni-1)*3 + 3 ) = fe((ni-1)*3 + 3 )-wf(3,1);
    fe((nj-1)*3 + 1 ) = fe((nj-1)*3 + 1 )-wf(4,1);
    fe((nj-1)*3 + 2 ) = fe((nj-1)*3 + 2 )-wf(5,1);
    fe((nj-1)*3 + 3 ) = fe((nj-1)*3 + 3 )-wf(6,1);
end