nodefenli.m

基于力密度法的索杆张力结构找形程序

function [Xfl,Xl,Yfl,Yl,Zfl,Zl,node_fixedno,node_freeno]=nodefenli(nodefix,nodenum,node_no,nfix_n,node_x,node_y,node_z)
nfix_n(nodefix+1)=-1;
j=1;fixed=1;free=1;
for i=1:1:nodenum
    if node_no(i,1)==nfix_n(j,1)
        node_fixedno(fixed)=node_no(i);
        Xfl(fixed,1)=node_x(i,1);
        Yfl(fixed,1)=node_y(i,1);
        Zfl(fixed,1)=node_z(i,1);
        fixed=fixed+1;
        j=j+1;
    else
        node_freeno(free)=node_no(i);
        Xl(free,1)=node_x(i,1);
        Yl(free,1)=node_y(i,1);
        Zl(free,1)=node_z(i,1);
        free=free+1;
    end
end
nfix_n(nodefix+1,:)=[];
fixed=fixed-1;free=free-1;   %fixed表示约束节点的数目，free表示自由节点的数目