目 录
摘要……………………………………………………………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
静定平面桁架影响线绘制
摘要：本课程设计介绍了平面静定桁架影响线软件编程作法。
它利用有限元的的思想，以maple为编程工具，将有限元思想与世界上最强大的符号计算与无与伦比的数值计算maple软件完美结合，快速、准确的获得结果。
关键词：maple 有限元 桁架影响线
第一章 题目的意义
有限元方法自20世纪50年代后期出现以来，就以其理论先进、功能强大以及与工程结构的复杂结构形式、材料性质和作用荷载的高度适应性，在工程结构力学计算中得到越来越广泛的应用。传统的人工计算力学方面，不仅耗费大量的人力、物力，而且计算时间长，精度低。在遇到某些更加困难的力学方程时，只能利用数值分析方法得到近似值。利用有限元思想与计算机结合，不仅能够克服上述困难，更甚将程序举一反三，运用自如。
- 1 -
第二章 理论分析及程序编制思路
2.1传统方法静定桁架影响线：
如图，做FG内力影响线：
①、作1-1截面，令P=1在截面左 侧移动，研究其右半部：
?M ?0DRB?2d?NFG?r1?0?NFG??2d1?RB???MDr1r1②、作1-1截面，令P=1在截面右 侧移动，研究其左半部：
?M ?0DRA?2d?NFG?r1?0?NFG??其他杆件内力影响线做法类似可得。
2d1?RA???MDr1r12.2有限元法编程作影响线：
2.2.1 影响线绘制基本思路
绘制桁架内力影响线，常采用静力法。即先计算出单位荷载作用在位移路径每一节点上时，目标杆件的内力。在绘制影响线时，只要早坐标系中标识出相应内力的大小，并用直线将相邻节点计算得的内力连接起来即可得到该杆件的内力影响线。
- 2 -
2.2.2 桁架内力计算
桁架结构的内力计算方法主要有：结点法、截面法、联合法。在具体计算时，规定内力符号以杆件受拉为正，受压为负。
结点法是截取桁架的一个结点为隔离体计算桁架内力的方法。结点上的荷载、支座反力和杆件轴力作用线都汇交于一点，组成了平面汇交力系，因此，结点法是利用平面汇交力系来求解内力的。在每个结点上所受到的力的合力为零。
对于任意杆件K，其两端结点分别为i、j，杆件轴力大小为Nk，如图所示。将杆件两端的受力分解到X、Y方向上可得
?Xi??Nk?cos?k??Yi??Nk?sin?k??Xj??Nk?cos?k???Yj??Nk?sin?k
对于任意的结点i，设其联接了杆件1，2，3，…，n，则
?Xi?a1N1?a2N2???anNnan??cos?n??Yi?b1N1?b2N2???bnNnbn??sin?n
对于静定平面桁架，设结点数为n，则共有n-m根杆，m个支座反力。提取系数可得2n×2n的系数矩阵。其中第2i-1行表示编号为i结点X方向分力，2i行表示i号结点Y方向分力。第j列表示编号
- 3 -
为j的杆件内力或支座反力在各结点上的投影系数。用[K]表示系数矩阵。
以各杆内力及支座反力构成一2n×1的矩阵，其中第j行表示编号为j的杆件内力或支座反力。以{R}表示。
各结点所受外力也可以构成一2n×1的矩阵， 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}是2n×1的矩阵，将{R1}、{R2}、…、
- 4 -