平面刚架分析.doc

平面框架结构的数值分析 与有限元分析实例 姓 名: 张 建 学 号: 联系方式:指导教师: 王 骑 2011年12月3日 题目:平面框架结构的数值分析与有限元分析 如下图所示平面框架结构,所受荷载如图所示,试分别用MATLAB数值分析软件和ANSYS有限元分析软件分析该结构各杆件的位移和内力。结构中各个截面的参数都为：E=210GPa,I=32×10-5m4，A=1.2×10-2m2。在ANSYS平台上，平面梁单元采用Beam3单元，不考虑剪切影响。 分析： （1）第一种方法：采用MATLAB进行数值分析 一、步骤 1、对结点和单元进行编号，选定整体坐标系和局部坐标系，如下图所示。 2、求出各单元在整体坐标系中的单元刚度矩阵，通过调用MATLAB的PlaneFrameElementStiffness（E，A，I，L，theta）函数实现。 3、将各单刚子块对号入座，形成结构原始刚度矩阵，通过调用MATLAB的PlaneFrameAssemble（K，k，i，j）函数实现。 4、计算非结点荷载作用下的各单元固端力、等效结点荷载及综合结点荷载，根据（《结构力学》上册-李廉锟 主编）表10-3等直杆单元的固端力查表，并在MATLAB中编程计算所需数据。 5、引入支承条件，修改结构的原始刚度方程。结点1和10为固定端，所以结点1和10处的线位移和角位移均为0，在原始刚度矩阵中删去与零位移对应的行和列，同时在结点位移列向量和结点外力列向量中删去相应的行，在MATLAB中通过命令k=K（4:27，4:27）便可得到修改后的结构的原始刚度矩阵。 6、解结构的刚度方程，求出结点位移，在MATLAB中通过命令u=k\f来求得所需数据。 7、计算各单元杆端力，在MATLAB中通过命令PlaneFrameElement Forces(E,A,I,L,theta,u)+各单元在局部坐标系下的固端力来求解。 二、计算结果数据整理 MATLAB计算结构结点位移单元编号ui（m）vi（m）φiuj（m）vj（m）φj10.0000 0.0000 0.0000 -0.0041 -0.0001 0.0032 2-0.0041 -0.0001 0.0032 -0.0096 -0.0002 0.0014 3-0.0096 -0.0002 0.0014 -0.0109 -0.0002 0.0011 4-0.0109 -0.0002 0.0011 -0.0114 0.0001 -0.0007 5-0.0114 0.0001 -0.0007 -0.0115 0.0000 0.0017 6-0.0115 0.0000 0.0017 -0.0096 0.0000 0.0020 7-0.0096 -0.0002 0.0014 -0.0096 -0.0014 -0.0008 8-0.0096 -0.0014 -0.0008 -0.0096 0.0000 0.0020 9-0.0096 0.0000 0.0020 -0.0040 0.0000 0.0031 10-0.0040 0.0000 0.0031 0.0000 0.0000 0.0000  MATLAB计算结构单元杆端力单元编号FNi（kN）FSi（kN）Mi（kN.m）FNj（kN）FSj（kN）Mj（kN.m）1103.8043 -89.2894 -198.4863 -103.8043 89.2894 19.9074 2103.8043 -89.2894 -19.9074 -103.8043 69.2894 -138.6715 341.7878 -44.3236 -3.2707 -41.7878 34.3236 -36.0529 453.3511 -8.8312 36.0529 -53.3511 -6.1688 -40.0464 5-6.1688 53.3511 40.0464 6.1688 24.5912 34.6739 618.2122 -17.6380 -34.6739 -18.2122 47.6380 2.0359 724.9659 62.0165 141.9421 -24.9659 -62.0165 44.1073 824.9659 12.0165 -44.1073 -24.9659 -12.0165