平面刚架静力分析有限元程序.docVIP

有限单元法 平面刚架静力分析程序 #includeio.h #includestring.h #includestdlib.h #includestdio.h #includemath.h void hbw(); void sncs(int nel); void fix(int np); void trmat(); void fis(int nel); void fpj(); void force(); void stiff(); void addsm(); void restr(); void matmul(); void soleq(); void outdis(); float sm[6][6],tsm[300][30],res[60][2],pj[300],t[6][6],d[6][6],r[300], fo[6],foj[6],pf[200][4],x[100],y[100],ae[10][3],sl,sn,cs,eal,eil; int nj,ne,nt,nr,npj,npf,nn,mbw,jel[100][2],nae[100],is[6]; FILE *infile,*outfile; /************主函数**************/ void main() { char name1[30],name2[30]; int i,j,nel,np; printf(please enter data-filename\n); scanf(%s,name1); printf(please enter result-filename\n); scanf(%s,name2); if((infile=fopen(name1,r))!=NULL) { fscanf(infile,%d%d%d%d%d%d,nj,ne,nt,nr,npj,npf); for(i=0;ine;i++) fscanf(infile,%d%d%d,jel[i][0],jel[i][1],nae[i]); for(i=0;inj;i++) fscanf(infile,%f%f,x[i],y[i]); for(i=0;int;i++) fscanf(infile,%f%f%f,ae[i][0],ae[i][1],ae[i][2]); for(i=0;inr;i++) fscanf(infile,%f%f,res[i][0],res[i][1]); } else { printf(the data-file not exit!); exit(1); } nn=3*nj; outfile=fopen(name2,w); if(outfile==NULL) {printf(the result-file not exist!);exit(1);} fprintf(outfile, statis analysis of plane frame\n); fprintf(outfile,input data\n); fprintf(outfile,************\n); fprintf(outfile,control data\n); fprintf(outfile,the num.of node:%3d\n,nj); fprintf(outfile,the num.of mem:%3d\n,ne); fprintf(outfile,the num.of type of section characteristic:%3d\n,nt); fprintf(outfile,the num.of restricted degrees of freedom:%3d\n,nr); fprintf(outfile,the num.of nodal loads:%3d\n,npj); fprintf(outfile,the num.of non-nodal loads:%3d\n,npf); fprintf(outfile,the num.of nodal degrees of freedom:%3d\n,nn); fprintf(outfile,information of mem.\n); fprintf(outfile,mem. start node end node type\n); for(i=0;ine;i++) fprintf(outfile,%5d%10d%