有限元程序设计.doc

有限元程序设计.doc有限元程序设计作业程序功能说明本程序包含四个子程序,其主要功能分别是:voidmatmatvoidmatvec调用结点及单元信息判断并循环形成总刚doublematvecdoublegs调用结点位移及结点信息和对应单元信息,回代求出并输出单元杆端力有限元程序流程图输入单元数输入结点数根据单元数形成动态数组存储单元及结点信息以及单元截面信息输入非结点荷载数有形成动态数组存储非结点荷载信息否根据非结点荷载信息判断并形成列向量输入结点荷载数否调用结点及单元信息判断并循环形成总刚形成总荷载列向量输入结点荷载信息并形成结点荷载列向量有形成动态数组存储结构总刚信息引入支承条件,对总刚及荷载列向量进行处理输入铰结点数及其编号 调用结点位移及结点信息和对应单元信息,回代求出并输出单元杆端力解线性方程组,求出结点位移铰结点的处理程序变量说明NE---------------------存储单元号,单元结点号,两端约束情况nelem--------------------单元数nnode---------------------结点数NM--------------------存储单元弹模,面积,惯性矩,长度,与总体X轴夹角nf-----------------------单元自由度nf1------------------------结点自由度K-----------------------结构的总刚度矩阵K0---------------------存储各单元单刚矩阵nelem1-----------------非结点荷载单元数nnode1-------------------有结点荷载结点数k-----------------------单元单刚矩阵nnode2-----------------有约束的结点数nnode3-----------------铰结点数T-----------------------单元坐标转换矩阵Po---------------------非结点荷载列向列量F----------------------单元非结点荷载的固端力nbg-------------------起始结点号nend------------------终止结点号NL--------------------单元的非结点荷载信息LD--------------------单元非结点荷载的大小Pa---------------------所有单元结点荷载的列向量LD1------------------存储各单元荷载信息P-----------------------所有荷载之和列向量NR---------------------存储支座约束的结点号F1----------------------存储各单元的杆端力NLD-----------------存储有结点荷载的结点号JJ-----------------------存储铰结的结点编号V-----------------------返回并存储结点位移列向量V1---------------------调用各单元的结点位移源程序//#include""#include<>#include<>#[6][6],X[6];voidmatmat(introw,intcol,doubleB1[6][6],doubleB2[6][6]);voidmatvec(introw,intcol,doublemat1[6][6],doubleY1[6]);doublematvec(introw,intcol,double**mat1,double*Y1,double*V);doublegs(introw,intcol,double**A,double*B,double*V);intmain(intargc,char*argv[]){doublem=1e8;intnelem,nnode,nelem1,nnode1,nnode2,nnode3,nf,nf1;inti1,i2,i3,nbg,nend,nbg1,nend1,nbg2,nend2,nbg3,nend3,N,N1,N2,n,n1;doubleL,a,E,A,I,q,T[6][6],T1[6][6],k[6][6],mat1[6][6];doubleP1[6],V1[6],F1[6],X1[6];cout<<"输入单元总数:";cin>>nelem;cout<<"输入结点总数:";cin>>nnode;int**NE;NE=newint*[5];for(i1=0;i1<nelem;i1++){NE[i1]=newint[5]; }cout<<"输入:单元编号,