2022年有限元程序设计终版 .pdf

目录摘要2 第一章题目的意义2 第二章理论分析及程序编制思路3 2.1传统方法静定桁架影响线3 2.2有限元法编程作影响线3 2.2.1 影响线绘制基本思路3 2.2.2桁架内力计算4 2.2.3桁架内力影响线的绘制5 第三章程序设计框图7 第四章程序编制8 4.1编制说明8 4.2程序代码8第五章计算实例10 第六章结论12 参考文献12 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 19 页 - - - - - - - - - - 1 - 静定平面桁架影响线绘制摘要：本课程设计介绍了平面静定桁架影响线软件编程作法。它利用有限元的的思想， 以 maple 为编程工具，将有限元思想与世界上最强大的符号计算与无与伦比的数值计算maple软件完美结合，快速、准确的获得结果。关键词： maple 有限元桁架影响线第一章题目的意义有限元方法自 20 世纪 50 年代后期出现以来，就以其理论先进、功能强大以及与工程结构的复杂结构形式、材料性质和作用荷载的高度适应性，在工程结构力学计算中得到越来越广泛的应用。传统的人工计算力学方面，不仅耗费大量的人力、物力，而且计算时间长，精度低。在遇到某些更加困难的力学方程时，只能利用数值分析方法得到近似值。利用有限元思想与计算机结合，不仅能够克服上述困难，更甚将程序举一反三，运用自如。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 19 页 - - - - - - - - - - 2 - 第二章理论分析及程序编制思路2.1 传统方法静定桁架影响线：如图，做 FG内力影响线：、作 1-1 截面，令 P=1在截面左侧移动，研究其右半部：、作 1-1 截面，令 P=1在截面右侧移动，研究其左半部：其他杆件内力影响线做法类似可得。2.2 有限元法编程作影响线：2.2.1 影响线绘制基本思路绘制桁架内力影响线， 常采用静力法。 即先计算出单位荷载作用在位移路径每一节点上时，目标杆件的内力。在绘制影响线时，只要早坐标系中标识出相应内力的大小，并用直线将相邻节点计算得的内力连接起来即可得到该杆件的内力影响线。0DM021rNdRFGBDBFGMrRrdN11120DM021rNdRFGADAFGMrRrdN1112名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 19 页 - - - - - - - - - - 3 - 2.2.2 桁架内力计算桁架结构的内力计算方法主要有：结点法、截面法、联合法。在具体计算时，规定内力符号以杆件受拉为正，受压为负。结点法是截取桁架的一个结点为隔离体计算桁架内力的方法。结点上的荷载、 支座反力和杆件轴力作用线都汇交于一点，组成了平面汇交力系，因此， 结点法是利用平面汇交力系来求解内力的。在每个结点上所受到的力的合力为零。对于任意杆件 K，其两端结点分别为i 、j ，杆件轴力大小为kN，如图所示。将杆件两端的受力分解到X、Y方向上可得对于任意的结点i ，设其联接了杆件1，2，3， n，则nnnninnnnibNbNbNbYaNaNaNaXsincos22112211对于静定平面桁架，设结点数为n，则共有 n-m根杆， m个支座反力。提取系数可得 2n2n 的系数矩阵。 其中第 2i-1 行表示编号为i 结点 X方向分力，2i 行表示 i 号结点 Y方向分力。第 j 列表示编号kkjkkjkkikkiNYNXNYNXsincossincos名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 19 页 - - - - - - - - - - 4 - 为 j 的杆件内力或支座反力在各结点上的投影系数。用K 表示系数矩阵。以各杆内力及支座反力构成一2n1 的矩阵，其中第 j 行表示编号为 j 的杆件内力或支座反力。以R 表示。各结点所受外力也可以构成一2n1 的矩阵， 2i-1行表示编号为 i 结点 X方向的外力，2i 行表示 i 号结点 Y方向外力。以P 表示。如此，则有 P=K R 。若已知外荷载 P, 则可通过桁架的系数矩阵求解桁架各杆件及支座反力 R 2.2.3 桁架内力影响线的绘制内力影响线即表示单位移动荷载作用下结构内力变化规律的图形，是研究移动荷载的基本工具。绘制桁架内力影响线时， 移动荷载有其移动的路径， 一般是桁架的上弦或下弦。 桁架通常承受结点荷载， 任一杆的轴力影响线在相邻结点之间为一直线。采用静力法做桁架影响线时， 需要计算出荷载作用在每一结点时杆件的内力，再用直线将计算得到的内力联接起来即可得到杆件的内力影响线。设移动荷载通过的路径上有m个结点，则当移动荷载作用在每一个结点时，都可以通过P=K R 计算得到该情况下各杆件内力R 。分别以 R1、R2、 Rm表示。在 n个结点的静定桁架中， R是 2n1的矩阵， 将R1、 R2、 、名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 19 页 - - - - - - - - - - 5 - Rm合成起来，可以用一个 2nm的矩阵RR来表示杆件内力的变化情况。则在这个总体内力矩阵中， 第 i 行第 j 列的元素表示当移动荷载作用在路径上第 j 个节点时，编号为 i 的杆件内力。 以此若要绘制任意杆件的影响线， 从该矩阵中提取出对应行的内力之，并以直线连接起来即可。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 19 页 - - - - - - - - - - 6 - 第三章程序设计框图开始静定桁架节点坐标、单元节点编号、约束信息节点坐标矩阵约束矩阵单位荷载路径矩阵单元矩阵各杆件沿荷载路径内力桁架影响线绘图结束计算处理系数矩阵循环名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 19 页 - - - - - - - - - - 7 - 第四章程序编制4.1 编制说明（1）程序使用 MAPLE 语言在 maple7.0 中编写（2）本程序适用于平面静定桁架影响线的绘制，不适用于超静定桁架的计算。（3）计算时，由于单位荷载的移动路径不唯一，应通过路径矩阵确定，路径矩阵中元素表示单位荷载依次通过的结点编号。（4）符号说明：参数：XY ：结点坐标矩阵；NE ：杆件单元矩阵；NC ：支座约束矩阵；P：荷载路径矩阵；n1,n2: 输出影响线杆件起止编号， 输出编号 n1 到 n2 的杆件的影响线，当 n1=n2=a时，输出第 a 根杆件的内力影响线。局部变量：ne：循环变量；kk：杆件 cos、sin 的值；te ：外力矩阵 P ；st ：系数矩阵 K; np：内力矩阵 R ；NP ：各荷载作用点下 R的合成；4.2 程序代码work:=proc(XY,NE,NC,P,n1,n2) local ne,kk,te,st,np,NP,x,pp: st:=Matrix(2*rowdim(XY),2*rowdim(XY),0): np:=Matrix(1,2*rowdim(XY),0): for ne to rowdim(NE) do: 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 19 页 - - - - - - - - - - 8 - kk:=normalize(XYNEne,2,1-XYNEne,1,1,XYNEne,2,2-XYNEne,1,2): #kk=cosa,sina st2*NEne,1-1,ne:=-kk1: st2*NEne,2-1,ne:=+kk1: st2*NEne,1,ne:=-kk2: st2*NEne,2,ne:=+kk2; end do; for ne to rowdim(NC) do: st2*NCne,1+NCne,2-2,rowdim(NE)+ne:=-1: end do; NP:=matrix(rowdim(P),2*rowdim(XY),0): for ne to rowdim(P) do: te:=Vector(2*rowdim(XY),0): xne:=XYPne,1,1: te2*Pne,1:=-1: np:=matrix(linsolve(evalm(st*1.0),te): copyinto(np,NP,ne,1): end do: NP:=transpose(NP): print(NP): pp:=matrix(n2-n1+1,1,0): for ne from 1 to n2-n1+1 do: ppne,1:=display( plot(XYP1,1,1,0,XYProwdim(P),1,1,0,style=line,thickness=1,color=black,axes=none), plot(XYP1,1,1,0,x1,NPn1-1+ne,1,style=line,thickness=2,color=blue), plot(xcoldim(NP),NPn1-1+ne,coldim(NP),XYProwdim(P),1,1,0,style=line,thickness=2,color=blue), plot(seq(xi,NPn1-1+ne,i,xi+1,NPn1-1+ne,i+1,i=1.coldim(NP)-1),style=line,thickness=2,color=blue), textplot(seq(xi,NPn1-1+ne,i,NPn1-1+ne,i,i=1.coldim(NP),color=red), textplot(seq(XYP1,1,1-1,0,i,i=cat(No ,n1-1+ne),color=black,thickness=2): od: display(pp); end proc: 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 19 页 - - - - - - - - - - 9 - 第五章计算实例实例一图示静定桁架，共18 个结点， 33 根杆， 3 个支座反力。单位荷载分别沿上弦（ 2、5、8、10、13、16、18 各点）和下弦（ 1、3、6、9、11、14、17 各点）移动，各杆件的内力影响线。程序数据输入见 附录 1，计算结果见 附录 2。从计算得到的影响线中提出任一杆件的影响线如移动荷载沿桁架下弦移动时编号为14 的杆件的内力影响线。再采用传统方法做该杆的影响线，结果如图。可以看到程序计算结果与传统方式计算得到的结果一致，可知程序计算结果是正确的。35（）名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 19 页 - - - - - - - - - - 10 - 实例二图示静定桁架，共14 个结点， 25 根杆， 3 个支座反力。单位荷载分别沿上弦（ 2、4、6、8、10、12、14 各点）和下弦（ 1、3、5、7、9、11、13 各点）移动，计算并绘制各杆件的内力影响线。桁架数据输入如下，结果略。XY1:=Matrix(0,0,0,3,3,0,3,3,6,0,6,3,9,0,9,3,12,0,12,3,15,0,15,3,18,0,18,3): NE1:=Matrix(1,2,1,3,2,3,2,4,3,4,3,5,4,5,4,6,5,6,5,7,6,7,6,8,7,8,7,9,7,10,8,10,9,10,9,11,9,12,10,12,11,12,11,13,11,14,12,14,13,14): NC1:=Matrix(1,1,1,2,13,2): 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 19 页 - - - - - - - - - - 11 - 第六章结论有限元法作为结构分析的一种有效方法已经被广泛采用，而有限元法与计算机程序的结合使得大型结构分析得以实现。通过以上实例，可以看出，通过有限元的方法编程计算桁架影响线的结果准确，与传统方式比较具有相当程度的优越性。效率高，结果精度高。在处理传统方式很难或无法求解的大型、复杂结构时， 通过有限元程序也能够进行计算。当然程序并不是完美的，仍有许多地方可以完善，这些将在以后的工作中完成。参考文献1 科学出版社maple9.0 符号处理及应用黎捷编著2 高等教育出版社结构力学 I 龙驭求 包世华 主编名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 19 页 - - - - - - - - - - 12 - 附录 1 XY:=Matrix(0,0,0,6,3,0,3,3,3,7,6,0,6,6,6,8,9,0,9,9,12,0,12,6,12,8,15,0,15,3,15,7,18,0,18,6): NE:=Matrix(1,3,1,4,1,2,2,5,3,4,4,5,3,6,4,6,4,7,5,8,6,7,7,8,6,9,7,9,7,10,8,10,9,10,9,11,9,12,10,12,10,13,11,12,12,13,11,14,11,15,12,15,13,16,14,15,15,16,14,17,15,17,16,18,17,18): NC:=Matrix(1,1,1,2,17,2): P1:=Matrix(2,5,8,10,13,16,18): P2:=Matrix(1,3,6,9,11,14,17): := XY00063033376066689099120126128150153157180186 := NE131412253445364647586778697971081091091191210121013111212131114111512151316141515161417151716181718 := NC1112172 := P125810131618 := P21369111417名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 13 页，共 19 页 - - - - - - - - - - 13 - 附录 2 P1路径下各杆件影响线总体内力矩阵 RR：0.8333.6667.5000.3333.16670.0.-1.179-.9430-.7072-.4715-.23570.-1.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.-1.0.0.0.0.0.0.8333.6667.5000.3333.16670.0.-.70720.0.0.0.0.0.-.4715-.9430-.7072-.4715-.23570.0.0.0.0.0.0.0.0.50000.0.0.0.0.0.0.-1.0.0.0.0.0.3333.6667.5000.3333.16670.0.-.3727-.74540.0.0.0.0.-.2357-.4715-.7072-.4715-.23570.0.0.0.0.0.0.0.0.3333.66670.6667.33330.0.1667.3333.5000.6667.33330.0.0.0.0.-.7454-.37270.0.-.2357-.4715-.7072-.4715-.23570.0.0.0.0.0.0.0.0.0.0.0.0.50000.0.0.0.0.-1.0.0.0.1667.3333.5000.6667.83330.0.0.0.0.0.-.70720.0.-.2357-.4715-.7072-.9430-.47150.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.-1.0.0.1667.3333.5000.6667.83330.0.-.2357-.4715-.7072-.9430-1.1790.0.0.0.0.0.0.0.0.0.0.0.0.0.-1.0.0.0.0.0.0.0.1.8333.6667.5000.3333.16670.0.1667.3333.5000.6667.83331.名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 14 页，共 19 页 - - - - - - - - - - 14 - 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 15 页，共 19 页 - - - - - - - - - - 15 - 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 16 页，共 19 页 - - - - - - - - - - 16 - P2路径下各杆件影响线总体内力矩阵 RR：0.8333.6667.5000.3333.16670.0.-1.179-.9430-.7072-.4715-.23570.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.1.0.0.0.0.0.0.0.0.0.0.0.0.0.8333.6667.5000.3333.16670.0.-.70720.0.0.0.0.0.-.4715-.9430-.7072-.4715-.23570.0.0.0.0.0.0.0.0.50001.0.0.0.0.0.0.0.0.0.0.0.0.3333.6667.5000.3333.16670.0.-.3727-.74540.0.0.0.0.-.2357-.4715-.7072-.4715-.23570.0.0.0.0.0.0.0.0.3333.66671.6667.33330.0.1667.3333.5000.6667.33330.0.0.0.0.-.7454-.37270.0.-.2357-.4715-.7072-.4715-.23570.0.0.0.0.0.0.0.0.0.0.0.1.50000.0.0.0.0.0.0.0.0.1667.3333.5000.6667.83330.0.0.0.0.0.-.70720.0.-.2357-.4715-.7072-.9430-.47150.0.0.0.0.0.0.0.0.0.0.0.0.1.0.0.0.0.0.0.0.0.0.1667.3333.5000.6667.83330.0.-.2357-.4715-.7072-.9430-1.1790.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.0.1.8333.6667.5000.3333.16670.0.1667.3333.5000.6667.83331.名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 17 页，共 19 页 - - - - - - - - - - 17 - 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 18 页，共 19 页 - - - - - - - - - - 18 - 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 19 页，共 19 页 - - - - - - - - -