资源描述:
《有限元程序设计作业》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
1、长沙理工大学有限元程序设计作业姓名:陆文林王子昱学号:071053205071054236班级:07级研究生专业:结构工程院系:土木与建筑学院第21页长沙理工大学有限元程序流程图输入单元数输入结点数根据单元数形成动态数组存储单元及结点信息以及单元截面信息输入非结点荷载数有形成动态数组存储非结点荷载信息否根据非结点荷载信息判断并形成列向量输入结点荷载数否调用结点及单元信息判断并循环形成总刚形成总荷载列向量输入结点荷载信息并形成结点荷载列向量有形成动态数组存储结构总刚信息第21页长沙理工大学引入支承条件,对总刚及荷载列向量进行处理输入铰结点数及其编号调用结点位移及结点信息和对
2、应单元信息,回代求出并输出单元杆端力解线性方程组,求出结点位移铰结点的处理第21页长沙理工大学程序变量说明NE---------------------存储单元号,单元结点号,两端约束情况nelem--------------------单元数nnode---------------------结点数NM--------------------存储单元弹模,面积,惯性矩,长度,与总体X轴夹角nf-----------------------单元自由度nf1------------------------结点自由度K-----------------------结构的总刚度
3、矩阵K0---------------------存储各单元单刚矩阵nelem1-----------------非结点荷载单元数nnode1-------------------有结点荷载结点数k-----------------------单元单刚矩阵nnode2-----------------有约束的结点数nnode3-----------------铰结点数T-----------------------单元坐标转换矩阵Po---------------------非结点荷载列向列量F----------------------单元非结点荷载的固端力nbg---
4、----------------起始结点号nend------------------终止结点号NL--------------------单元的非结点荷载信息LD--------------------单元非结点荷载的大小Pa---------------------所有单元结点荷载的列向量LD1------------------存储各单元荷载信息P-----------------------所有荷载之和列向量NR---------------------存储支座约束的结点号F1----------------------存储各单元的杆端力NLD----------
5、-------存储有结点荷载的结点号JJ-----------------------存储铰结的结点编号V-----------------------返回并存储结点位移列向量第21页长沙理工大学V1---------------------调用各单元的结点位移源程序//#include"stdafx.h"#include#include#definePI3.14159265doublemat3[6][6],X[6];voidmatmat(introw,intcol,doubleB1[6][6],doubleB2[6][6]);vo
6、idmatvec(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,n
7、bg2,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