有限元上机报告.doc

《有限元上机报告.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、在以常应变三角形单元计算弹性力学平面问题原理的基础上，编写可计算平面应力问题，亦可计算平面应变问题的通用程序。给出程序设计的流程图，程序源代码和程序应用的工程算例及其结果分析。一、程序流程图此程序采用本人比较熟悉C语言编写。C语言相比于FORTRAN语言拥有更强大的实用编程环境，同时由于提供了指针变量、动态内存分配函数和结构体变量等，使其编写的软件在维护性、可读性和内存利用效率方面具有明显的优势。1.程序总框图图1程序总框图整个程序由一个主调主程序(主函数main（）)和七个子程序组成，其中数据子程序DATA()用来接受输入的参数和变量。

2、单刚子程序ELEST()为含有参数的函数，功能控制参数iask可取1、2、3，分别计算相应单元的面积、应力矩阵和单元刚度矩阵，主程序不直接调用它，而通过总刚子程序、载荷子程序和求应力子程序间接调用。总刚子程序TOTSTI()用来合成总刚矩阵，载荷子程序LOAD()用来计算合成载荷，支承子程序SUPPOR()用来引入约束，解方程子程序SOLVEQ()用来求解并输出各个结点的位移，求应力子程序STRESS()用来计算和输出应力、主应力及主平面角。2.主函数及各子程序流程图图2voidmain()和voidDATA()图3voidELEST(i

3、ntmeo,intiask)图4voidTOTSTI()图5voidLOAD()图6voidSUPPOR()图7voidSOLVEQ()图8voidSTRESS()二、程序源代码/*有限元计算程序*/#include#includeintnj,ne,nz,ndd,ind,nj2;intjm[100][3],nzc[200],npj1,npj;floateo,un,gama,te,ae;floatcjz[100][2],pj[100][2];floatb[3][6],d[3][3],s[3][6],eke[

4、6][6],tkz[200][20],p[200];voidDATA()/*数据输入函数*/{inti,j;printf("Qingshuru6gecanshu");scanf("%d,%d,%d,%d,%d,%d",&nj,&ne,&nz,&ndd,&npj,&ind);nj2=nj*2;npj1=npj+1;getchar();printf("Qingshuru4gejisuancanshu");scanf("%f,%f,%f,%f",&eo,&un,&gama,&te);getchar();printf("nj=%dne=%

5、dnz=%dndd=%dnpj=%dind=%d",nj,ne,nz,ndd,npj,ind);printf("eo=%fun=%fgama=%fte=%f",eo,un,gama,te);getchar();printf("QingshuruJMjuzheng");for(i=0;i

6、)scanf("%f",&cjz[i][j]);getchar();}printf("QingshuruNZCjuzheng");for(i=0;i

7、r(j=0;j<2;j++){printf("%f",cjz[i][j]);if(j%2==1)printf("");}}getchar();}voidELEST(intmeo,intiask)/*单元刚度矩阵*/{intie,je,me,i,j,k;floatcm,bm,cj,bj;ie=jm[meo-1][0];je=jm[meo-1][1];me=jm[meo-1][2];cm=cjz[je-1][0]-cjz[ie-1][0];bm=cjz[ie-1][1]-cjz[je-1][1];cj=cjz[ie-1][0]-cjz[m

8、e-1][0];bj=cjz[me-1][1]-cjz[ie-1][1];ae=(bj*cm-bm*cj)/2;if(iask>1){for(i=0;i<3;i++){for(j=0;j<6;