基于matlab桁架分析,基于MATLAB平面桁架有限元分析研究

第 30 卷 第 8 期 嘉应学院学报( 自然科学) Vol． 30 No． 08 2012 年 8 月 JOURNAL OF JIAYING UNIVERSITY( Natural Science) Aug． 2012 基于 MATLAB 平面桁架有限元分析研究 李远瑛，张德生 收稿日期:2012 －05 －26 基金项目:广东省自然科学基金项目( S2011010001679) ; 2011 年梅州市科技计划项目作者简介:李远瑛( 1973 － ) ，女，广东平远人，副教授，硕士，主要研究方向: 结构理论研究与计算． ( 嘉应学院 土木工程学院，广东 梅州 514015) 摘 要:基于平面桁架有限元分析的基本原理，利用 MATLAB 软件对某平面桁架编制了计算程序，得到了结构的内力计算结果、变形图和结构的动力特性，并将计算结果与 ANSYS 有限元分析软件的结果进行比较． 结果表明，二者的计算结果一致． 利用 MATLAB 可以编制结构有限元分析程序，计算速度快，数据准确，并能为复杂结构有限元分析提供参考． 关键词:MATLAB; 平面桁架; 有限元; ANSYS 中图分类号:TU311． 2 文献标识码: A 文章编号: 1006 －642X( 2012) 08 －0029 －05 0 前言 有限元法是根据变分原理求解数学及物理问题的数值计算方法［1］． 随着计算机技术的快速发展，有限元分析得到了广泛应用． MATLAB 是效率高的工程计算软件工具，可以实现编程计算、数据图形化集成． 本文利用 MATLAB 软件，利用有限元分析的基本原理对平面桁架进行受力分析，得到了 结构各杆件的内力、应力、结构变形及结构动力特性，并将计算结果同 ANSYS 有限元分析软件结果进行比较． 结果表明，利用 MATLAB 可以编制结构有限元分析程序，计算速度快，数据准确，并能为复杂 结构有限元分析提供参考． 1 桁架结构有限元分析的基本原理 桁架结构在屋架、桥梁等工程中应用很多． 在桁架结构中，各杆件的连接处近似为铰接，杆件绕结点可以自由转动，结构中的各个杆件只受轴力， 属于一维问题，各杆件内的位移线性变化，杆件内 的位移分布函数是精确的［2］． 在进行结构离散化时，不需要对杆件再细分单元，直接以杆件为分析单元． 1． 1 局部坐标系下的杆件单元刚度矩阵 桁架中任意杆件 ij 为分析单元，建立结构的整体坐标系为 xoy，局部坐标系为 xoy，如图 1 所示． 一维杆件单元在局部坐标中的单元刚度矩阵为 Ke = EA l 1 0 －1 0 0 0 0 0 －1 0 1 0              0 0 0 0 图 1 平面桁架杆件单元整体坐标与局部坐标 1． 2 坐标转换矩阵 对于平面桁架单元，通过坐标系旋转可以完成 从单元局部坐标系到结构整体坐标系的转换． 杆件 局部坐标系与整体坐标系的夹角为 α，转换矩阵为 1． 3 整体坐标系下的杆件单元刚度矩阵 整体坐标系下桁架单元刚度矩阵与局部坐标 系下的单元刚度矩阵关系为 ke = TTkeT，即 1． 4 整体坐标下结构总刚度矩阵 K 在结构的整体分析中涉及许多单元及联接它们的结点，为避免混淆，必须对各单元、结点及结点位移( 水平位移和竖向位移) 进行编号，由单元刚度矩阵通过单元定位向量集成结构总刚度矩阵 K，矩 阵的阶数为结点总数的 2 倍． 单元定位向量由单元 两端的结点位移编号确定． 1． 5 建立结构的结点荷载向量 P 1． 6 建立结构的刚度方程 K△ = P △为结构整体