弹性力学平面问题有限元程序.docVIP

有限单元法 弹性力学平面问题有限元程序 #define NODE_NUM 300 #define ELE_NUM 200 #define NODE2 NODE_NUM*2 #includestdlib.h #includestdio.h #includemath.h int nj,ne,nm,nz,npj,jm[ELE_NUM][4],nzc[60],nj2,bw,lxm; float ea,e,um,oul,te,xy[NODE_NUM][4],EA[5],p[NODE2],pj,F[ELE_NUM],etm[5][5]; float ek[7][7]; double **K,*f; int BAND2(int n,int m, double **K,double *f); double **TwoArrayAlloc(int,int); void TwoArrayFree(double **); FILE *outf; /*-------------------------------------------*/ double **TwoArrayAlloc(int r,int c) { double *x,**y; int n; x=(double*)calloc(r*c,sizeof(double)); y=(double**)calloc(r,