航空航天结构中的有限元方法编程作业.docx

航空航天结构中的有限元方法编程大作业姓名熊蕾班级110514学录一、使用程序介绍1二、数据文件生成程序1三、计算程序2四、例题求解4例1:4例2：5例3：6五、程序不足之处与改进方法71、数据文件的可读性72、数据长度混乱73、输入量很难从图形得出84、算法的性质8六、参考文献8使用程序介绍由于有限元刚度方程求解过程是一个矩阵运算的过程，使用matlab编写程序会比C语言更加方便。MATLAB是近年来得到迅速发展的数学软件，它将高性能的数值计算和可视化集成在一起，并提供了大量内置函数，被广泛地应用于科学计算、控制系统、信息处理等领域的分析、仿真和实际工作。数据文件生成程序为了方便计算程序读取数据，我编写了以下程序来生成txt数据文件，这样所得的数据文件中数据储存的格式正是计算时需要的格式。format longnn = input(请输入节点总数\n)%节点数nm = input(请输入单元数\n)%单元数L = input(请输入各杆件的长度，用分号隔开，单位为:m\n)%各单元长度（单位：米）E = input(请输入弹性模量，单位为:GPa\n)%弹性模量A = input(请输入杆件截面积，单位为:m^2\n)%截面积Q = input(请输入各杆倾角，用分号隔开，单位为:度\n)%各杆倾角JM = input(请输入各杆对应的节点(同一杆的节点序号用空格隔开，不同杆之间用分号隔开）\n)%各杆对应的节点(同一杆的节点序号用空格隔开，不同杆之间用分号隔开）nf=4%外载荷数量F=input(请输入外载荷，单位为:kN\n)%外载荷矩阵bc=input(请输入边界条件，被约束自由度序号的元素置1，未被约束置0)%边界条件，有位移约束的置1data=zeros(11,2*nn)data(1,1)=nndata(2,1)=nmfor i=1:nmdata(3,i)=L(i,1)enddata(4,1)=Edata(5,1)=Afor i=1:nmdata(6,i)=Q(i,1)endfor i=1:nmdata(7,i)=JM(i,1)endfor i=1:nmdata(8,i)=JM(i,2)enddata(9,1)=nffor i=1:2*nndata(10,i)=F(1,i)endfor i=1:2*nn data(11,i)=bc(1,i)endfilename=input(请输入保存的文件名\n,s)fid=fopen(filename,w);fprintf(fid,%f %f %f %f %f %f %f %f %f %f %f %f\n,data)fclose(fid)运行这个.m文件，在指引下输入各项数据，最后能在工作目录下生成文件名可以自定义的txt数据文件。由于时间仓促，编写时设置的输入量写成了杆件长度和杆件的倾斜角，今后可对这些进行改进，输入节点坐标，然后由可得出相应杆件的长度和倾角。另外，也可以不使用这个小程序生成数据文件，只要存储数据时按照格式存成一个矩阵，计算程序就能读取相关数据进行计算。计算程序编写的计算程序如下，matlab对矩阵的操作十分方便，首先读入数据并将它们存入一个名为data的矩阵，再用这个矩阵对各参数进行赋值，组装总体刚度矩阵也可以很轻松地进行，因为编程时设置的字长比较短，为确保计算精度，在引入边界条件时用了置1法。filename=input(请输入生成的数据文件的文件名\n,s)fid=fopen(filename,r);[data,count]=fscanf(fid,%f %f %f %f %f %f %f %f %f %f %f,[11 inf]);fprintf(1,%f %f %f %f %f %f %f %f %f %f %f,data);fclose(fid)%读入数据nn=data(1,1)nm=data(2,1)for i=1:nmL(i,1)=data(3,i)endE=data(4,1)A=data(5,1)for i=1:nmQ(i,1)=data(6,i)endfor i=1:nmJM(i,1)=data(7,i)endfor i=1:nmJM(i,2)=data(8,i)endnf=data(9,1)for i=1:2*nnF(i,1)=data(10,i)endfor i=1:2*nn bc(i,1)=data(11,i)endndof = 2;%平面桁架节点自由度数为2q = Q*pi/180;a = cos(q);b = sin(q);K2=zeros(nn*ndof);for e=1:nm%对每一个单元（杆）k=A*E/L(e)*[a(e)^2,a(e)*b(e);a(e)*b(e),b(e)^2];K=[k,-k;-k,