结构程序分析设计桁架c语言

《结构程序分析设计桁架c语言》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、#include#includemain(){FILE*fp;inti,j,NN,NE,NF,n,m,z[20],a[20]⑵,g[20][4],c[20];floatE,A,Cx,Cy,l,T,Xi,Xj,Yi,Yj,K[20][20]={0},k[4][4]={0},F[4][4]={0},y,t,d,ui,uj,vi,vj,Ui,Uj,Vi,Vj,fxi,fxj,fyi,fyj;clrscr();fp=fopen(Hdata3.txt,,,HrH);fscanf(fp,“％d”,&NN);fscanf(fp,”％d”,&NE);fscanf(fp,“％d

2、”,&NF);for(i=0;i<=2*NF-1;i++)fscanf(fp,”％d”,&z[i]);for(i=();i<=NN-l;i++)for(j=0;j<=l;j++)fscanfCfp/^d'^&afiirj]);for(i=0;i<=NN-l;i++)for(j=0;j<=3;j++)fscanf(fp；'%dM,&g[i][j]);for(i=0;i<=2*NN-1;i++)fscanf(fp；'%dn,&c[i]);(close(fp);fp=fopen(nanswer3.txt",nw");fprintf(fp,HNNtNEtNFn);fprintf(fp,”％d

3、t%dt%d”,NN,NE,NF);fprintf(fp,“SUPPRESSEDDISPLACIMENTPOINTu);for(i=0;i<=2*NF-l;i++){fprintf(fp,"%dt",z[i]);if(i%2!=0)fprintf(fp,MH);}fprintf(fp,HNO.tX-COORD.tY-COORDM);for(i=0;i<=NN-l;i++){fprintf(fp,"%dt",i+1);for(j=0;j<=l;j++)fprintf(fp,M%dtM,a[iJ[jJ)；fprintf(fp,MH);}fprintf(fp,HIt

4、JtAtEM);for(i=0;i<=NN-l;i++){for(j=O;j<=3;j++)fprintf(fp,M%dfg[i][j]);fprintf(fp,MH);}fprintf(fp,HVECTOROFNODALLOADS');fprintf(fp,MNO.tX-LOAD.tY-LOAD.M);for(i=0;i<=NN-1;i++)fprintf(fp,"%dt%dH,c[2*i],c[2*i+1]);fprintf(fp,MH);}fprintf(fp,HTHENODALDISPLACEMENTARE:u);for(i=l;i<=NN;i++)f

5、or(j=i+l;j<=NN;j++)for(t=0;t<=3;t++){if(i==g[t][O]&&j==g[t][l]){Xi=a[i-l][01;Yi=a[i-l][l];Xj=a

6、j-l][O];Yj=aU-l][l];l=sqrt((Xi-Xj)*(Xi-Xj)+(Yi-Yj)*(Yi-Yj));Cx=(Xj-Xi)/l;Cy=(Yj-Yi)/l;E=g[纽3]；A=g[t][21;T二E*A/1;k[0][0]=T*Cx*Cx;k[O]LlJ=T*Cx*Cy;k[0]⑵=T*(-l)*Cx*Cx;k[0]⑶=T*(-l)*Cx*Cy;k[l][0]=T*Cx*Cy;k[l][l]

7、=T*Cy*Cy;k[l][2]=T*(-l)*Cx*Cy;k[l]⑶=T*(-l)*Cy*Cy;k[2][0]=T*(-l)*Cx*Cx;k[2][l]=T*(-l)*Cx*Cy;k⑵[2]=T*Cx*Cx;k[2]⑶=T*Cx*Cy;k[3][0]=T*(-l)*Cx*Cy;k[3][l]=T*(-l)*Cy*Cy;k⑶L2J=T*Cx*Cy;k⑶⑶二T*Cy*Cy;K[2Sl][2*i-l]=K[2*i-l][2*i-l]+k[0][0];K[2*i-l][2*i]=K[2*i-l][2*i]+k[0][l];K[2*i][2*i-1]=K[2*i][2*i-1]+k[1][0];K[2

8、*i][2*i]=K[2*i][2*i]+k[1][1];K[2*j-1][2*j-1]=K[2*j-1][2*j-1]+k[2][2];K[2*j-l]r2*jl=K[2*j-H[2*j]+kR]⑶；K[2*j][2*j-l]=K[2*j][2*j-l]+k[3][2];K0j][2*j]=K[2*j][2*j]+k⑶⑶;K[2*i-l][2*j-l]=K[2*i-l][2*j-l]+k[0J[