变分原理及有限元大作业.docx

八、附录?附录一：程序代码??#includestdio.h?#includemath.h?#includestdlib.h?#includeiomanip.h?#includeiostream.h?float?**float_two_array_malloc(int?m,int?n)???//实型二维动态数组分配?{?float?**a;?????int?i,j;?a=(float?**)malloc(m*sizeof(float?*));?????for?(i=0;im;i++)?????{???????? ?a[i]=(float?*)malloc(n*sizeof(float));for?(j=0;?jn;?j++)???{?? ??a[i][j]=0;?????}?????}??return?a;}float?*float_one_array_malloc(int?n)???//实型一维动态数组分配?{?????float?*a;?????int?j;????? a=(float?*)malloc(n*sizeof(float));??for?(j=0;?jn;?j++)??{?? ?a[j]=0;} return a ;}int?**int_two_array_malloc(int?m,int?n)???//整型二维动态数组分配{???? ?int?**a;?????int?i,j;?????a=(int?**)malloc(m*sizeof(int?*));?????for?(i=0;im;i++)?{?a[i]=(int?*)malloc(n*sizeof(int));???for?(j=0;?jn;?j++)???{????a[i][j]=0;???}?????}??return?a;?}?int?*int_one_array_malloc(int?n)???//整型一维动态数组分配{?int?*a;?????int?j;?a=(int?*)malloc(n*sizeof(int));??for?(j=0;?jn;?j++)??{?a[j]=0;??}?????return?a;?}?//全局变量定义int?np,ne,nr1,nr2,nd,nf,n,ndf,ld,nm,nu1,nu2;?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为变带宽一维组中元素个数???????????????????????//IS表示单元节点位移在结构位移阵列中对应的序号float?**EK,**s,*A,**ns,**mat;???//EK为总体坐标系下单元的刚度矩阵//A存储结构刚度矩阵的一维变带宽数组???????????????????????//s存放各单元应力中心点的应力???????//ns存放各节点的应力?void?input()????//输入函数定义?{??int?i,j;??FILE?*fp1;?? fp1=fopen(input.dat,r);?//要输入的内容放在input.dat中??if(fp1==NULL)??{??? printf(cannt?open?the?file?!\n);???exit(1);??}?//控制信息读入?fscanf(fp1,%d,np);??fscanf(fp1,%d,ne);?fscanf(fp1,%d,nr1);????fscanf(fp1,%d,nr2);???fscanf(fp1,%d,nd);?????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);??????n=nf*np;?ndf=nf*nd;?//?printf(%d\t%f\t%d\t%f\n,nu1,u1,nu2,u2);?//动态数组生成??x=float_one_array_malloc(np);?????y=float_one_array_malloc(np);??nrr=int_one_array_malloc(nr1+nr2);???//位移约束为零的节点，nr1--x方向约束，nr2---y方向约束??nu=int_one_array_mall