平面四节点等参单元有限元程序

exit(1);
}
//控制信息读入
fscanf(fp1,"%d",&np);
fscanf(fp1,"%d",&ne);
fscanf(fp1,"%d",&nr1);
fscanf(fp1,"%d",&nr2);
fscanf(fp1,"%d",&nd);
exit(1);
}
fprintf(fp2,"位移为：\n"); //输出位移
for(i=0;i<n;i+=2)
fprintf(fp2,"u%d=%f v%d=%f\n",i/2+1,p[i],i/2+1,p[i+1]);
fprintf(fp2,"\n各单元应力：Sx,\tSy,\tSxy,\tS1,\tS2,\tSmises\n"); //输出单元应力
float *x,*y,*p,**pp,u1,u2; //x,y,me,nrr,p为输入信息
float coords[2][4],det,fun[4],pn[2][4],xjac[2][2],rjac[2][2],dnx[2][4]; //求单刚定义的变量
int *LD,**me,*nrr,*nu,*IS,nn; //nn为变带宽一维组中元素个数
void output() //结果输出函数
{ int i,j;
FILE *fp2;
fp2=fopen("output.dat","w"); //运算结果放在output.dat中
if(fp2==NULL)
{
printf("cann't open the file !\n");
}
}
for(i=0;i<nr1;i++) //读入固定位移约束条件
{
fscanf(fp1,"%d",&nrr[i]);//x方向约束的节点号
nrr[i]=(nrr[i]-1)*2; //转化为自由度序号，abaqus模型节点编号从1开始，C语言中从0开始
}
//ns存放各节点的应力
void input() //输入函数定义
{
int i,j;
FILE *fp1;
fp1=fopen("input.dat","r"); //要输入的内容放在input.dat中
if(fp1==NULL)
{
printf("cann't open the file !\n");
}
void tecplot() //结果生成tecplot文件
{ int i,j;
float bl=1.0; //位移显示比例
FILE *fp3;
fp3=fopen("tecplot.plt","w"); //运算结果放在tecplot.plt中
if(fp3==NULL)
{
a[i][j]=0;
}
}
return a;
}
float *float_one_array_malloc(int n) //实型一维动态数组分配
{
float *a;
int j;
a=(float *)malloc(n*sizeof(float));
#include<stdio.h>
#include<math.h>
#include<stdlib.h>
#include<iomanip.h>
#include<iostream.h>
float **float_two_array_malloc(int m,int n) //实型二维动态数组分配
{
int *a;
int j;
a=(int *)malloc(n*sizeof(int));
for (j=0; j<n; j++)
{
a[j]=0;
}
return a;
}
//全局变量定义
int np,ne,nr1,nr2,nd,nf,n,ndf,ld,nm,nu1,nu2;
}
for(i=0;i<nu1;i++) //读入非零位移载荷
{
fscanf(fp1,"%d",&nu[i]);//x方向约束的节点号
nu[i]=(nu[i]-1)*2; //转化为自由度序号，abaqus模型节点编号从1开始，C语言中从0开始
}
for(i=0;i<nu2;i++)
fscanf(fp1,"%d",&nf);
fscanf(fp1,"%d",&ld);
fscanf(fp1,"%d",&nm);
fscanf(fp1,"%d",&nu1);
fscanf(fp1,"%f",&u1);
fscanf(fp1,"%d",&nu2);
fscanf(fp1,"%f",&u2);
for(i=0;i<nm;i++) //读入材料信息
{
for(j=0;j<4;j++)
{
fscanf(fp1,"%f",&mat[j][i]);
}
mat[0][i]=mat[0][i]-1;mat[1][i]=mat[1][i]-1; //abaqus模型节点编号从1开始，C语言中从0开始
fprintf(fp3," ZONE N=%d,E=%d,F=FEPOINT,ET=quadrilateral\n",np,ne);
for(j=0;j<np;j++)
{
fprintf(fp3,"%.4f\t%.4f\t%.4f\t%.4f\t",x[j]+bl*p[j*2],y[j]+bl*p[j*2+1],p[j*2],p[j*2+1]);
{
float **a;
int i,j;
a=(float **)malloc(m*sizeof(float *));
for (i=0;i<m;i++)
{
a[i]=(float *)malloc(n*sizeof(float));
for (j=0; j<n; j++)
fscanf(fp1,"%f",&y[i]);
}
for(i=0;i<ne;i++)
{
for(j=0;j<nd;j++)
{
fscanf(fp1,"%d",&me[j][i]);
me[j][i]=me[j][i]-1; //abaqus模型节点编号从1开始，C语言中从0开始
{
fscanf(fp1,"%d",&nu[nu1+i]);//y方向约束的节点号，abaqus模型节点编号从1开始，C语言中从0开始
nu[nu1+i]=(nu[nu1+i]-1)*2+1;
}
fclose(fp1);
printf("读入信息完成......\n");
}
{
printf("cann't open the file !\n");
exit(1);
}
fprintf(fp3," TITLE = ' test '\n");
fprintf(fp3," VARIABLES=x,y,Ux,Uy,Sx,Sy,Sxy,Smax,Smin,Smises\n");
fscanf(fp1,"%f%f",&pp[i][0],&pp[i][1]);
for(i=0;i<n;i++)
p[i]=0.0;
for(i=0;i<ld;i++)
{
j=int(pp[i][0]);
p[j]=pp[i][1];
}
LD[0]=0;
for(k=0;k<np;k++)
{
ig=np;
me=int_two_array_malloc(nd,ne); //单元节点的总体编号
mat=float_two_array_malloc(4,nm);//单元材料数组，先定义为各向同性材料
//原始数据读入
for(i=0;i<np;i++)
{
fscanf(fp1,"%f",&x[i]);
for (i=0;i<m;i++)
{
a[i]=(int *)malloc(n*sizeof(int));
for (j=0; j<n; j++)
{