loadff.m

来自「教学程序」· M 代码 · 共 14 行

function f=loadff
%----形成节点荷载
global   nfp_number nfp  node_number npe_number %（定义全局变量）
f = zeros(node_number*6,1) ;%（创建f）
for inf=1:1:nfp_number
  n = nfp( inf, 1 ) ;
  d = nfp( inf, 2 ) ;
  f((n-1)*6 + d ,1) = nfp( inf, 3 ) ;
end
fe=equload;
for i=1:1:6*node_number
    f(i,1)=f(i,1)+fe(i,1);
end
return