有限元程序设计.doc

目录摘要……………………………………………………………2第一章题目的意义…………………………………………2第二章理论分析及程序编制思路……………………………3 ……………………………3 ………………………………3 ……………………………3 ……………………………………4 ……………………………5第三章程序设计框图…………………………………………7第四章程序编制………………………………………………8 ………………………………………………8 ………………………………………………8第五章计算实例……………………………………………10第六章结论…………………………………………………12参考文献………………………………………………………12静定平面桁架影响线绘制摘要:本课程设计介绍了平面静定桁架影响线软件编程作法。它利用有限元的的思想,以maple为编程工具,将有限元思想与世界上最强大的符号计算与无与伦比的数值计算maple软件完美结合,快速、准确的获得结果。关键词:maple有限元桁架影响线第一章题目的意义有限元方法自20世纪50年代后期出现以来,就以其理论先进、功能强大以及与工程结构的复杂结构形式、材料性质和作用荷载的高度适应性,在工程结构力学计算中得到越来越广泛的应用。传统的人工计算力学方面,不仅耗费大量的人力、物力,而且计算时间长,精度低。在遇到某些更加困难的力学方程时,只能利用数值分析方法得到近似值。利用有限元思想与计算机结合,不仅能够克服上述困难,更甚将程序举一反三,运用自如。:如图,做FG内力影响线:①、作1-1截面,令P=1在截面左侧移动,研究其右半部:②、作1-1截面,令P=1在截面右侧移动,研究其左半部:其他杆件内力影响线做法类似可得。:,常采用静力法。即先计算出单位荷载作用在位移路径每一节点上时,目标杆件的内力。在绘制影响线时,只要早坐标系中标识出相应内力的大小,并用直线将相邻节点计算得的内力连接起来即可得到该杆件的内力影响线。:结点法、截面法、联合法。在具体计算时,规定内力符号以杆件受拉为正,受压为负。结点法是截取桁架的一个结点为隔离体计算桁架内力的方法。结点上的荷载、支座反力和杆件轴力作用线都汇交于一点,组成了平面汇交力系,因此,结点法是利用平面汇交力系来求解内力的。在每个结点上所受到的力的合力为零。对于任意杆件K,其两端结点分别为i、j,杆件轴力大小为,如图所示。将杆件两端的受力分解到X、Y方向上可得对于任意的结点i,设其联接了杆件1,2,3,…,n,则对于静定平面桁架,设结点数为n,则共有n-m根杆,m个支座反力。提取系数可得2n×2n的系数矩阵。其中第2i-1行表示编号为i结点X方向分力,2i行表示i号结点Y方向分力。第j列表示编号为j的杆件内力或支座反力在各结点上的投影系数。用[K]表示系数矩阵。以各杆内力及支座反力构成一2n×1的矩阵,其中第j行表示编号为j的杆件内力或支座反力。以{R}表示。各结点所受外力也可以构成一2n×1的矩阵,2i-1行表示编号为i结点X方向的外力,2i行表示i号结点Y方向外力。以{P}表示。如此,则有{P}=[K]×{R}。若已知外荷载{P},则可通过桁架的系数矩阵求解桁架各杆件及支座反力{R},是研究移动荷载的基本工具。绘制桁架内力影响线时,移动荷载有其移动的路径,一般是桁架的上弦或下弦。桁架通常承受结点荷载,任一杆的轴力影响线在相邻结点之间为一直线。采用静力法做桁架影响线时,需要计算出荷载作用在每一结点时杆件的内力,再用直线将计算得到的内力联接起来即可得到杆件的内力影响线。设移动荷载通过的路径上有m个结点,则当移动荷载作用在每一个结点时,都可以通过{P}=[K]×{R}计算得到该情况下各杆件内力{R}。分别以{R1}、{R2}、…、{Rm}表示。在n个结点的静定桁架中,{R}是2n×1的矩阵,将{R1}、{R2}、…、{Rm}合成起来,可以用一个2n×m的矩阵{RR}来表示杆件内力的变化情况。则在这个总体内力矩阵中,第i行第j列的元素表示当移动荷载作用在路径上第j个节点时,编号为i的杆件内力。以此若要绘制任意杆件的影响线,从该矩阵中提取出对应行的内力之,并以直线连接起来即可。第三章程序设计框图开始静定桁架节点坐标、单元节点编号、约束信息节点坐标矩阵约束矩