有限元源程序.doc

//**程序开始**#include<>#include<>#defineNE6#defineNJ6#defineNZ6#defineNPJ0#defineNPF1#defineNJ318#defineDD12###-2#*fp;//这是输入参数的初始化和定义全局变量**intjm[NE+1][3]={{0,0,0},{0,1,3},{0,3,2},{0,3,5},{0,6,5},{0,5,4},{0,2,4}};doublegc[NE+1]={,,,,,,};/*杆长数组*/doublegj[NE+1]={,,,,,,};/*杆角数组*/doublemj[NE+1]={,,,,,,};/*面积数组*/doublegx[NE+1]={,I0,I0,I0,I0,I0,I0};/*各杆惯性矩数组*/intzc[NZ+1]={0,1,2,3,16,17,18};doublepf[NPF+1][5]={{0,0,0,0,0},{0,-45,,,}},pj[NPJ+1][3]={,,};doublekz[NJ3+1][NJ3+1],p[NJ3+1];doublepe[7],f[7],f0[7],t[7][7];doubleke[7][7],kd[7][7];//**kz[][]-整体刚度矩阵//**ke[][]-整体坐标下的单元刚度矩阵//**kd[][]-局部坐标下的单元刚度矩阵//**t[][]-坐标变换//**这是函数的声明voidjdugd(int);voidzb(int);voidgdnl(int);voiddugd(int);//**主程序开始voidmain(void){ inti,j,k,e,d,dh,h,ii,jj,hz,a1,b1,m,l,dl,zl,z,j0; doublecl,wy[7]; intIM,IN,jn; //****************************** //文件定义 //****************************** if(fp=fopen("c:\\","w")==NULL) printf("no\n"); else printf("yes\n"); fp=fopen("c:\\","w");//**********************************//功能:形成矩阵p//*********************************** if(NPJ>0) { for(i=1;i<=NPJ;i++) { j=(int)pj[i][2]; p[j]=pj[i][1]; } } if(NPF>0) { for(i=1;i<=NPF;i++) { hz=i; gdnl(hz); e=(int)pf[hz][3]; zb(e); for(j=1;j<=6;j++) { pe[j]=; fo