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

航空航天结构中的有限元方法 编程大作业 姓名熊蕾班级110514学号航空航天结构中的有限元方法 编程大作业 PAGE \* MERGEFORMATi 目录  TOC \h \z \t 标题 2,1,标题 3,2  HYPERLINK \l _Toc388918797 一、 使用程序介绍  PAGEREF _Toc388918797 \h 1  HYPERLINK \l _Toc388918798 二、 数据文件生成程序  PAGEREF _Toc388918798 \h 1  HYPERLINK \l _Toc388918799 三、 计算程序  PAGEREF _Toc388918799 \h 2  HYPERLINK \l _Toc388918800 四、 例题求解  PAGEREF _Toc388918800 \h 4  HYPERLINK \l _Toc388918801 例1:  PAGEREF _Toc388918801 \h 4  HYPERLINK \l _Toc388918802 例2：  PAGEREF _Toc388918802 \h 5  HYPERLINK \l _Toc388918803 例3：  PAGEREF _Toc388918803 \h 6  HYPERLINK \l _Toc388918804 五、 程序不足之处与改进方法  PAGEREF _Toc388918804 \h 7  HYPERLINK \l _Toc388918805 1、 数据文件的可读性  PAGEREF _Toc388918805 \h 7  HYPERLINK \l _Toc388918806 2、 数据长度混乱  PAGEREF _Toc388918806 \h 7  HYPERLINK \l _Toc388918807 3、 输入量很难从图形得出  PAGEREF _Toc388918807 \h 8  HYPERLINK \l _Toc388918808 4、 算法的性质  PAGEREF _Toc388918808 \h 8  HYPERLINK \l _Toc388918809 六、 参考文献  PAGEREF _Toc388918809 \h 8 PAGE \* MERGEFORMAT8 使用程序介绍 由于有限元刚度方程求解过程是一个矩阵运算的过程，使用matlab编写程序会比C语言更加方便。 MATLAB是近年来得到迅速发展的数学软件，它将高性能的数值计算和可视化集成在一起，并提供了大量内置函数，被广泛地应用于科学计算、控制系统、信息处理等领域的分析、仿真和实际工作。 数据文件生成程序 为了方便计算程序读取数据，我编写了以下程序来生成txt数据文件，这样所得的数据文件中数据储存的格式正是计算时需要的格式。 format long nn = 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)%边界条件，有位移约束的置1 data=zeros(11,2*nn) data(1,1)=nn data(2,1)=nm for i=1:nm data(3,i)=L(i,1) end data(4,1)=E data(5,1)=A for i=1:nm data(6,i)=Q(i,1) end for i=1:nm data(7,i)=JM(i,1) end for i=1:nm data(8,i)=JM(i,2) end data(9,1)=nf for i=1:2*nn data(10,i)=F(1,i) end for i=1:2*nn data(11,i)=bc(1,i) end filename=input(请输