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

1、#include <stdio.h>#include <math.h>main()file *fp;int i,j,nn,ne,nf,n,m,z20,a20,g204,c20;floate,a,cx,cy,l,t,xi,xj,yi,yj,k2020=0,k44=0,f44=0,y,t,d,ui,uj,vi,vj,ui,uj,vi,vj,fx i,fxj,fyi,fyj;clrscr();fp=fopen(hdata3.txt,hrh);fscanf(fp,“d”,&nn);fscanf(fp,”d”,&ne);fscanf(fp,“d”,&nf)

2、;for(i=0;i<=2*nf-1 ;i+)fscanf(fp,”d”,&zi);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,&gij);for(i=0;i<=2*nn-1 ;i+)fscanf(fp；'%dn,&ci);(close(fp);fp=fopen(nanswer3.txt",nw");fprintf

3、(fp,hnntnetnfnn);fprintf(fp,”dt%dt%dn”,nn,ne,nf);fprintf(fp,“suppressed displaciment pointnu);for(i=0;i<=2*nf-l;i+)fprintf(fp,"%dt",zi);if(i%2!=0)fprintf(fp,mnh);fprintf(fp,hno.tx-coord.ty-coordnm);for(i=0;i<=nn-l;i+)fprintf(fp,"%dt",i+1);for(j=0;j<=l;j+)fprintf(fp,m%dtm

4、,aijjj)；fprintf(fp,mnh);fprintf(fp,hitjtatenm);for(i=0;i<=nn-l;i+)for(j=o;j<=3;j+)fprintf(fp,m%dfgij);fprintf(fp,mnh);fprintf(fp,hvector of nodal loads');fprintf(fp,mno.tx-load.ty-load.nm);for(i=0;i<=nn-1 ;i+)fprintf(fp,"%dt%dh,c2*i,c2*i+1); fprintf(fp,mnh);fprintf(fp,hthe nodal di

5、splacement are:nu);for(i=l;i<=nn;i+)for(j=i+l ;j<=nn;j+)for(t=0;t<=3;t+)if(i=gto&&j=gtl)xi=ai-l01;yi=ai-ll;xj=a|j-lo;yj=au-ll; l=sqrt(xi-xj)*(xi-xj)+(yi-yj)*(yi-yj);cx=(xj-xi)/l;cy=(yj-yi)/l;e=g纽3；a=gt21;t 二 e*a/1;k00=t*cx*cx;kollj=t*cx*cy;k0=t*(-l)*cx*cx;k0=t*(-l)*cx*cy;kl0=t*cx*cy

6、;kll=t*cy*cy;kl2=t*(-l)*cx*cy;kl=t*(-l)*cy*cy;k20=t*(-l)*cx*cx;k2l=t*(-l)*cx*cy;k 2 =t*cx*cx;k2=t*cx*cy;k30 =t*(-l)*cx*cy;k3 l=t*(-l)*cy*cy;k l2j=t*cx*cy;k 二 t*cy*cy;k2sl2*i-l=k2*i-l2*i-l+k00;k2*i- l2*i=k2*i- l2*i+k 0l;k2*i 2*i-1 =k2*i 2*i-1 +k 1 0;k2*i 2*i=k2*i 2*i+k 1 1 ;k2*j-1 2*j-1 =k2*j-1 2*j-1

7、 +k2 2;k2*j-lr2*jl=k2*j-h2*j+kr；k2*j2*j-l=k2*j2*j-l+k32;k0j2*j=k2*j2*j+k ;k2*i-l2*j-l=k2*i-l2*j-l+k0j2;k2*il2*j=k2怜 l2*j+k0;k2*i2*j-l=k2*i2*j-l+kl2;k2*i2*j=k2*i2*j+kl3;k2*j-l2*i-l=k2*j-l 2*i-l+k20;k2*j-1 2衍=k2*j 12呵+k1 ;k2*j 2*i-1 =k2*j 2*i-1 +k 30j;k2*j2*i=k2*j2*i+k 1;for(i=0;i<=7;i+)ki+l9=cij;f

8、or(i=l;i<=8;i+)for(j=l；j<=9；j+) fprintf(fp,%3f役 kij);fprintf(fp,hnh);for(i=2;i<=8;i=i+2) if(ki9=0)kii=l;for(j=i-l;j>=l;j-)kiu=o;for(j=i+1 ;j<=9;j+)kijuj=o； for(j=i-l;j>=l;j-) kji=o;for(j=i+l;j<=9;j+) k|ji二 0;t=i-1;for(j=t-l;j>=l;j-) krtjl=0;for(j=t+l;j<=9;j+) ktj二 0;for(j=

9、t-l;j>=l;j-) kut=0;for(j=t+l ；j<=9;j+) kjt=0; for(i=l;i<=8;i+)for(j=l;j<=9;j+)fprintf(fp,h%2.3f skij); fprintf(fp,"nh);*/n=2*nn;m=2*nn+l; for(i=l;i<=n;i+)if(i=n)for(j=i;j<=m;j+)kiju=(kiu/y); for(t=i-l;t>=l;t-) d=kgi;for(j=i;j<=m;j+)ktu=(ktuj-kiu*d);elsey=kiij；for(j=i;j&l

10、t;=m;j+)kiu=(kiu/y);for(t=i+1 ;t<=n;t+)d=kti;for(j=i；j<=m;j+)ktj=(krtj.kij*d);for(t=i-l;t>=l;t)d=kti;for(j=i;j<=m;j+)ktiu=(krtrj-kiui*d)； for(i=l;i<=8; i+)for(j=l;j<=9;j+)fprintf(fp,n%.3f kij)； fprintf(fp,mnh);*/fprintf(fp,mnodet u:tv;nh);for(i= 1 ;i<=4;i+)fprintf(fp；,(%d)t%ft%f

11、n,i,k2*i-l9,k2*i9);fprintf(fp,hforces in each members');fprintf(fp,mmembertforces');for(i=l;i<=nn;i+)for(j=i+1 ;j<=nn;j+)for(t=0;t<=3;t+)if(i=gto&&j=gtl)xi=ai-l0;yi=ai-ll;xj=au-lo;yj=au-ulu；lrsqrt(xt-xj)mxrxj)+(yi,yj)£y7-yj)；cxaxjxx一；cyayy<ehg 三2;aug 三 e; the 沃 a、一；k

12、o=oitcx- so=lcy;20=2110;s0=3llp ku=on(')© 竺二三hex; k=2pss3h0;s2=0p keullp竺 2=2ncx;km=3ncy竺3=0p竺3=l0;53=2 (二 fey kl3=3itcx-un-kt3*=9k vh-k-2*=9k un;k2*h=9k vkej=9kuli-ko=o*ui+so=l*vi+ko=2*uj+k©=3*vj; vn-k=oui+k二=1 = vi+辽 1 =2匚uj+k= =3匚 vj - ujus2=0*ui+s2=l 匚 vi+k-2=2uj+s2=3*vj 一 vjh 竺 3=0=ui+竺 3=二*<+竺3=2=£+龙 3=3 匸 vh f2oiihfoeuo;fo=2lr(qt;fc=3hphfhpf三三ho;fu=2n0;fu=3uo;feolret;fellrpfe2ht;22=3o;33=0