2022年有限元程序设计终版 .pdf
( 4.5 )《2022年有限元程序设计终版 .pdf》由会员分享，可在线阅读，更多相关《2022年有限元程序设计终版 .pdf（19页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、目录摘要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 - 静定平面桁架影响线绘制摘要
2、：本课程设计介绍了平面静定桁架影响线软件编程作法。它利用有限元的的思想， 以 maple 为编程工具，将有限元思想与世界上最强大的符号计算与无与伦比的数值计算maple软件完美结合，快速、准确的获得结果。关键词： maple 有限元桁架影响线第一章题目的意义有限元方法自 20 世纪 50 年代后期出现以来，就以其理论先进、功能强大以及与工程结构的复杂结构形式、材料性质和作用荷载的高度适应性，在工程结构力学计算中得到越来越广泛的应用。传统的人工计算力学方面，不仅耗费大量的人力、物力，而且计算时间长，精度低。在遇到某些更加困难的力学方程时，只能利用数值分析方法得到近似值。利用有限元思想与计算机结合
3、，不仅能够克服上述困难，更甚将程序举一反三，运用自如。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 19 页 - - - - - - - - - - 2 - 第二章理论分析及程序编制思路2.1 传统方法静定桁架影响线：如图，做 FG内力影响线：、作 1-1 截面，令 P=1在截面左侧移动，研究其右半部：、作 1-1 截面，令 P=1在截面右侧移动，研究其左半部：其他杆件内力影响线做法类似可得。2.2 有限元法编程作影响线：2.2.1 影响线绘制基本思路绘制桁架内力影响线
4、， 常采用静力法。 即先计算出单位荷载作用在位移路径每一节点上时，目标杆件的内力。在绘制影响线时，只要早坐标系中标识出相应内力的大小，并用直线将相邻节点计算得的内力连接起来即可得到该杆件的内力影响线。0DM021rNdRFGBDBFGMrRrdN11120DM021rNdRFGADAFGMrRrdN1112名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 19 页 - - - - - - - - - - 3 - 2.2.2 桁架内力计算桁架结构的内力计算方法主要有：结点法、
5、截面法、联合法。在具体计算时，规定内力符号以杆件受拉为正，受压为负。结点法是截取桁架的一个结点为隔离体计算桁架内力的方法。结点上的荷载、 支座反力和杆件轴力作用线都汇交于一点，组成了平面汇交力系，因此， 结点法是利用平面汇交力系来求解内力的。在每个结点上所受到的力的合力为零。对于任意杆件 K，其两端结点分别为i 、j ，杆件轴力大小为kN，如图所示。将杆件两端的受力分解到X、Y方向上可得对于任意的结点i ，设其联接了杆件1，2，3， n，则nnnninnnnibNbNbNbYaNaNaNaXsincos22112211对于静定平面桁架，设结点数为n，则共有 n-m根杆， m个支座反力。提取系数
6、可得 2n2n 的系数矩阵。 其中第 2i-1 行表示编号为i 结点 X方向分力，2i 行表示 i 号结点 Y方向分力。第 j 列表示编号kkjkkjkkikkiNYNXNYNXsincossincos名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 19 页 - - - - - - - - - - 4 - 为 j 的杆件内力或支座反力在各结点上的投影系数。用K 表示系数矩阵。以各杆内力及支座反力构成一2n1 的矩阵，其中第 j 行表示编号为 j 的杆件内力或支座反力。以R
7、 表示。各结点所受外力也可以构成一2n1 的矩阵， 2i-1行表示编号为 i 结点 X方向的外力，2i 行表示 i 号结点 Y方向外力。以P 表示。如此，则有 P=K R 。若已知外荷载 P, 则可通过桁架的系数矩阵求解桁架各杆件及支座反力 R 2.2.3 桁架内力影响线的绘制内力影响线即表示单位移动荷载作用下结构内力变化规律的图形，是研究移动荷载的基本工具。绘制桁架内力影响线时， 移动荷载有其移动的路径， 一般是桁架的上弦或下弦。 桁架通常承受结点荷载， 任一杆的轴力影响线在相邻结点之间为一直线。采用静力法做桁架影响线时， 需要计算出荷载作用在每一结点时杆件的内力，再用直线将计算得到的内力联
8、接起来即可得到杆件的内力影响线。设移动荷载通过的路径上有m个结点，则当移动荷载作用在每一个结点时，都可以通过P=K R 计算得到该情况下各杆件内力R 。分别以 R1、R2、 Rm表示。在 n个结点的静定桁架中， R是 2n1的矩阵， 将R1、 R2、 、名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 19 页 - - - - - - - - - - 5 - Rm合成起来，可以用一个 2nm的矩阵RR来表示杆件内力的变化情况。则在这个总体内力矩阵中， 第 i 行第 j 列的
9、元素表示当移动荷载作用在路径上第 j 个节点时，编号为 i 的杆件内力。 以此若要绘制任意杆件的影响线， 从该矩阵中提取出对应行的内力之，并以直线连接起来即可。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 19 页 - - - - - - - - - - 6 - 第三章程序设计框图开始静定桁架节点坐标、单元节点编号、约束信息节点坐标矩阵约束矩阵单位荷载路径矩阵单元矩阵各杆件沿荷载路径内力桁架影响线绘图结束计算处理系数矩阵循环名师资料总结 - - -精品资料欢迎下载 -
10、- - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 19 页 - - - - - - - - - - 7 - 第四章程序编制4.1 编制说明（1）程序使用 MAPLE 语言在 maple7.0 中编写（2）本程序适用于平面静定桁架影响线的绘制，不适用于超静定桁架的计算。（3）计算时，由于单位荷载的移动路径不唯一，应通过路径矩阵确定，路径矩阵中元素表示单位荷载依次通过的结点编号。（4）符号说明：参数：XY ：结点坐标矩阵；NE ：杆件单元矩阵；NC ：支座约束矩阵；P：荷载路径矩阵；n1,n2: 输出影响线杆件起止编号，
11、输出编号 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: 名师资料总结
12、- - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 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) d
13、o: 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 n
14、e 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-
- 配套讲稿：
如PPT文件的首页显示word图标，表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制：
部分文档作品中含有的国旗、国徽等图片，仅作为作品整体效果示例展示，禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词：
- 2022年有限元程序设计终版 2022 有限元 程序设计
淘文阁 - 分享文档赚钱的网站所有资源均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。
限制150内