todate.m

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

function todate(node_all,nodenum,nodefree,nodefix,cab_num,Xall,Yall,Zall,bianhao,nfix_n,cab_all,filename)

for i=1:1:nodenum
    for j=1:1:nodenum
        if node_all(i,1)==bianhao(j)
            node_all(i,2)=Xall(j);
            node_all(i,3)=Yall(j);
            node_all(i,4)=Zall(j);
        end
    end
end
    
smodelfile=strrep(filename,'date.txt','update.txt');
fid=fopen(smodelfile,'wt');     %自应力模态文件名smodel
if fid<0
    errordlg({'文件打开出错'});
else
    fprintf(fid,'%d\n%d\n%d\n%d\n',nodenum,nodefree,nodefix,cab_num);
     fprintf(fid,'%3d%15.5f%15.5f%15.5f%15.5f%15.5f%15.5f\n',node_all');
    fprintf(fid,'%d\n',nfix_n);
    fprintf(fid,'%5d%5d%5d%10.4g%10f%15.8f%10d\n',cab_all');
end
fclose(fid);