outputs.m

来自「基于力密度法的索杆张力结构找形程序」· M 代码 · 共 26 行

function outputs(t,E,A,l,uz,w,cab_num,x,y,z,r,s0,Etype,bianhao,cab_n,filename1)


for i=1:1:cab_num
    if Etype(i)==0
        s(i,1)=t(i)*s0(i)/(E(i)*A(i))+s0(i);
        derts(i,1)=s(i)-s0(i);
    elseif Etype(i)==1
        derts(i,1)=t(i)*l(i)*(l(i)+uz(i)^2/l(i)+w(i)^2*r(i)^2/(12*t(i)^2*l(i)))/(E(i)*A(i)*r(i));
        s(i,1)=fs(l(i),r(i),uz(i),s0(i),E(i),A(i),w(i),t(i))+s0(i)+derts(i,1);
    end
    %     derts(i)=h(i)*(l(i)+uz(i)^2/l(i))/(E(i)*A(i));
end
X=[bianhao,x,y,z]';

smodelfile=strrep(filename1,'.txt','变形后坐标.txt');
fid=fopen(smodelfile,'wt');     %自应力模态文件名smodel
if fid<0
    errordlg({'文件打开出错'});
else
    fprintf(fid,'           单元号    单元初始长度     变形后长度       变化长度       杆件径向力\n');
    fprintf(fid,'%15d%15.5f%15.5f%15.5f%15.5f\n',[cab_n,s0,s,derts,t]');
    fprintf(fid,'\n\n\n          节点号码      X坐标值       Y坐标值         Z坐标值\n');
    fprintf(fid,'%15d%15.4f%15.4f%15.4f\n',X);
end
fclose(fid);