有限元程序设计002.doc

有 限 元 程 序 设 计 作 业 程序功能说明 本程序包含四个子程序，其主要功能分别是： void matmat void matvec 调用结点及单元信息判断并循环形成总刚 double matvec double gs 调用结点位移及结点信息和对应单元信息，回代求出并输出单元杆端力 有限元程序流程图 有 否 否 有 程序变量说明 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 stdafx.h #include iostream.h #include math.h #define PI 3double mat3[6][6],X[6]; void matmat(int row,int col,double B1[6][6],double B2[6][6]); void matvec(int row, int col,double mat1[6][6],double Y1[6]); double matvec(int row, int col,double ** mat1,double * Y1,double * V); double gs(int row,int col,double ** A, double * B,double * V); int main(int argc, char* argv[]) { double m=1e8; int nelem,nnode,nelem1,nnode1,nnode2,nnode3,nf,nf1; int i1,i2,i3,nbg,nend,nbg1,nend1,nbg2,nend2,nbg3,nend3,N,N1,N2,n,n1; double L,a,E,A,I,q,T[6][6],T1[6][6],k[6][6],mat1[6][6]; double P1[6],V1[6],F1[6],X1[6]; cout输入单元总数:; cinnelem; cout输入结点总数:; cinnnode; int * *NE; NE=new int * [5]; for(i1=0;i1nelem;i1++)