您好，欢迎访问三一刀客
当前位置：首页 > 建筑/环境 > 建筑制度 > MATLAB计算平面桁架体系的程序
globall;E=1e6;A=1;p=fopen('INPUTDATA.txt','r');NP=fscanf(p,'%d',1);NF=fscanf(p,'%d',1);MP=zeros(NP,2);fori=1;NPMP(i,:)=fscanf(p,'%d',2);endNE=fscanf(p,'%d',1);ME=zeros(NE,2);fori=1:NEME(i,:)=fscanf(p,'%d',2)endNR=fscanf(p,'%d',1);NRR=zeros(NR,2);fori=1:NRNRR(i,:)=fscanf(p,'%d',2);endNL=fscanf(p,'%d',1);LL=zeros(NL,2);fori=1:NLLL(i,:)=fscanf(p,'%d',2);endNDF=NF*NP;K=zeros(NDF,NDF);fori=1:NEl(i)=sqrt((MP(ME(i,2),1)-MP(ME(i,1),1))^2+(MP(ME(i,2),2)-MP(ME(i,1),2))^2);v=[MP(ME(i,2),1)-MP(ME(i,1),1),MP(ME(i,2),2)-MP(ME(i,1),2)]/l(i);ax=v*[1;0];ay=v*[0;1];T(:,:,i)=[ax,ay,0,0;0,0,ax,ay];ke=(E*A/l(i))*[1,-1;-1,1];Ke=T(:,:,i)'*ke*T(:,:,i);K(2*ME(i,1)-1:2*ME(i,1),2*ME(i,1)-1:2*ME(i,1))=...K(2*ME(i,1)-1:2*ME(i,1),2*ME(i,1)-1:2*ME(i,1))+Ke(1:2,1:2);K(2*ME(i,1)-1:2*ME(i,1),2*ME(i,2)-1:2*ME(i,2))=...K(2*ME(i,1)-1:2*ME(i,1),2*ME(i,2)-1:2*ME(i,2))+Ke(1:2,3:4);K(2*ME(i,2)-1:2*ME(i,2),2*ME(i,1)-1:2*ME(i,1))=...K(2*ME(i,2)-1:2*ME(i,2),2*ME(i,1)-1:2*ME(i,1))+Ke(3:4,1:2);K(2*ME(i,2)-1:2*ME(i,2),2*ME(i,2)-1:2*ME(i,2))=...K(2*ME(i,2)-1:2*ME(i,2),2*ME(i,2)-1:2*ME(i,2))+Ke(3:4,3:4);endP=zeros(2*NP,1);fori=1:NLP(LL(i,1))=LL(i,2);endfori=1:NRj=NRR(i,1);K(j,j)=10^30;endfori=1:NRj=NRR(i,1);P(j)=10^30*NRR(i,2);endu=K^(-1)*P;fprintf('ёчЧФУЙ¶ИО»ТЖОЄ:/n');fprintf('%1.4f/n',u);fprintf('ёчёЛµҐФЄДЪБ¦ОЄ:/n');fori=1:NED=(E/(1-0.3^2))*[1,0.3,0;0.3,1,0;0,0,(1-0.3)/2];B=[-1/l(i),0,1/l(i),0;0,0,0,0;0,-1/l(i),0,1/l(i)];uii=[u(2*ME(i,1));u(2*ME(i,1));u(2*ME(i,2)-1);u(2*ME(i,2))];ui=T(:,:,i)*uii;F=E*(ui(2)-ui(1))/l(i);fprintf('%3.4f/n',F);end
三一刀客所有资源均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。
扫描二维码
笨笨的乖乖兔
本文标题：MATLAB计算平面桁架体系的程序
链接地址：https://www.111doc.com/doc-6339403 .html