基于MATLAB的平面刚架静力分析.doc

  1. 1、本文档共17页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
基于MATLAB的平面刚架静力分析

基于MATLAB的平面刚架静力分析 为了进一步理解有限元方法计算的过程,本文根据矩阵位移法的基本原理应用MATLAB编制计算程序对以平面刚架结构进行了静力分析。本文还利用ANSYS大型商用有限元分析软件对矩阵位移法的计算结果进行校核,发现两者计算结果相当吻合,验证了计算结果的可靠性。 问题描述 如图1所示的平面刚架,各杆件的材料及截面均相同,E=210GPa,截面为0.12×0.2m的实心矩形,现要求解荷载作用下刚架的位移和内力。 图1 二、矩阵位移法计算程序编制 为编制程序方便考虑,本文计算中采用“先处理法”。具体的计算步骤如下。 对结构进行离散化,对结点和单元进行编号,建立结构(整体)坐标系和单元(局部) 对结点位移分量进行编码,形成单元定位向量; 建立按结构整体编码顺序排列的结点位移列向量,计算固端力、等效结点荷载 及综合结点荷载列向量; 计算个单元局部坐标系的刚度矩阵,通过坐标变换矩阵形成整体坐标系下的单元刚度矩阵 ; 利用单元定位向量形成结构刚度矩阵; 按式 求解未知结点位移; 计算各单元的杆端力。 根据上述步骤编制了平面刚架的分析程序。程序中单元刚度矩阵按下式计算。 转换矩阵则按下式计算。 计算程序框图如图2所示,具体的程序代码见附录1。 图2 MATLAB矩阵分析法程序框图 三、解题步骤 取整体坐标系如图3所示,对结构进行离散化,对结点和单元进行编号如图4所示,局部坐标系用单元中箭头的方向表示,原始数据如下: 图3 图4 刚架结点输入矩阵为, x=[0 0;0 -5;1.63 -6.37;4 -5;4 -1;4 2]; 各单元定位向量为, locvec1=[1 2 3 0 0 0]; locvec2=[1 2 3 4 5 6]; locvec3=[4 5 6 7 8 9]; locvec4=[7 8 9 10 11 12]; locvec5=[10 11 12 0 0 0]; 输入截面参数, E=2.1e11;%E=210GPa a=0.12; %矩形截面长0.12m b=0.2; %矩形截面宽0.2m 输入整体坐标系下各单元结点荷载列阵, f(1,:)=zeros(1,6); f(2,:)=[0 0 0 0 40e3 0]; f(3,:)=zeros(1,6); f(4,:)=[0 0 0 -50e3 0 0]; f(5,:)=zeros(1,6); 输入整体坐标系下单元1等效节点荷载 q=10e3; %10kN/m fe=[0.5*q*l(1),0,-q*l(1)^2/12,0.5*q*l(1),0,q*l(1)^2/12]; 由此计算得到平面刚架整体坐标系下的结点位移(m), d= 0.0035 0.0000 -0.0004 0.0030 -0.0005 -0.0004 0.0027 0.0000 0.0016 -0.0051 0.0000 -0.0006 各个单元的杆端力如表1所示, 表1 各单元杆端力   单元 1 2 3 4 5 i端 Fx(kN) -17917.047 17917.05 17917.05 17917.05 -32083 Fy(kN) 17507.3731 -17507.4 22492.63 22492.63 22492.63 Mz(kN·m) 1897.83076 -1897.83 2092.833 -26668.3 44999.85 j端 Fx(kN) -32082.953 -17917 -17917 32082.95 32082.95 Fy(kN) -17507.373 -22492.6 -22492.6 -22492.6 -22492.6 Mz(kN·m) -37312.596 -2092.83 26668.34 -44999.8 51249.01 四、计算结果校核 在ANSYS中使用beam3单元,按照如图4所示的离散结构建立平面刚架模型施加约束和荷载,得到的有限元模型如图5所示。计算分析后得到结构的轴力图、剪力图和弯矩图如图6、7、8所示,命令流见附录2。 图5 有限元模型 图6 轴力图(kN) 图7 剪力图(kN) 图8 弯矩图(kN·m) 从ANSYS计算结果中提取各结点位移、内力,并与矩阵位移法分析的结果比较,得到表2、3。 表2 各节点位移比较 节点号 项目 矩阵位移法 ANSYS 误差 1 Ux(m) 0 0 0 Uy(m) 0 0 0 R

文档评论(0)

haihang2017 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档