有限元计算程序源代码

/*  有限元计算程序 */

#include <stdio.h>
#include <math.h>
int nj,ne,nz,ndd,ind,nj2;
int jm[100][3],nzc[200],npj1,npj;
float eo,un,gama,te,ae;
float cjz[100][2],pj[100][2];
float b[3][6],d[3][3],s[3][6],eke[6][6],tkz[200][20],p[200];

void DATA()  /*数据输入函数*/
{ int i,j;
  printf("Qing shu ru 6 ge can shu\n");
  scanf("%d,%d,%d,%d,%d,%d",&nj,&ne,&nz,&ndd,&npj,&ind);
  nj2=nj*2;
  npj1=npj+1;
  getchar();
  printf("Qing shu ru 4 ge ji suan can shu\n");
  scanf("%f,%f,%f,%f",&eo,&un,&gama,&te);
  getchar();
  printf("nj=%d  ne=%d  nz=%d  ndd=%d  npj=%d  ind=%d\n",nj,ne,nz,ndd,npj,ind);
  printf("eo=%f  un=%f  gama=%f  te=%f\n",eo,un,gama,te);
  getchar();
  printf("Qing shu ru JM ju zheng\n");
  for(i=0;i<ne;i++)
   { for(j=0;j<3;j++)
      scanf("%d",&jm[i][j]);
      getchar();
   }
  printf("Qing shu ru CJZ ju zheng\n");

  for(i=0;i<nj;i++)
   { for(j=0;j<2;j++)
      scanf("%f",&cjz[i][j]);
      getchar();
   }
  printf("Qing shu ru NZC ju zheng\n");
   for(i=0;i<nz;i++)
      scanf("%d",&nzc[i]);
      getchar();
  printf("Qing shu ru PJ ju zheng\n");
  for(i=0;i<npj1;i++)
   { for(j=0;j<2;j++)
      scanf("%f",&pj[i][j]);
      getchar();
   }
   getchar();
  printf("CJZ ju zheng ru xia\n");
  for(i=0;i<nj;i++)
   { for(j=0;j<2;j++)
      {
       printf("%f    ",cjz[i][j]);
       if(j%2==1)
       printf("\n");
      }
   }
  getchar();
}

void ELEST(int meo,int iask)/*单元刚度矩阵*/
{ int ie,je,me,i,j,k;