《有限元基础教程》_【MATLAB算例】3.2.5(2)四杆桁架结构的有限元分析(Bar2D2Node).docx

【MATLAB  算例】  3.2.5(2)  四杆桁架结构的有限元分析  (Bar2D2Node) 如图  3-8  所示的结构，各个杆的弹性模量和横截面积都为  E 29.5  104 N / mm2  , 100mm2 。试基于 MATLAB 平台求解该结构的节点位移、单元应力以及支反力。 图 3-8 四杆桁架结构 解答 ：对该问题进行有限元分析的过程如下。 （1） 结构的离散化与编号 对该结构进行自然离散， 节点编号和单元编号如图 3-8 所示，有关节点和单元的信息见 表 3-1~ 表 3-3。 （2）计算各单元的刚度矩阵 ( 基于国际标准单位 ) 建立一个工作目录，将所编制的用于平面桁架单元分析的 4 个 MATLAB 函数放置于该 工 作 目 录 中 ， 分 别 以 各 自 函 数 的 名 称 给 出 文 件 名 ， 即 ： Bar2D2Node_Stiffness ， Bar2D2Node_Assembly ，Bar2D2Node_Stress，Bar2D2Node_Forces 。然后启动 MATLAB ，将 工作目录设置到已建立的目录中，在 MATLAB 环境中，输入弹性模量 E、横截面积 A，各 点坐标 x1,y1,x2,y2,x3,y3,x4,y4, 角度 alpha 1, alpha 2 和 alpha 3,然后分别针对单元 1，2， 3 和 4，调用 4 次 Bar2D2Node_Stiffness ，就可以得到单元的刚度矩阵。相关的计算流程如下。 E=2.95e11; A=0.0001; x1=0; y1=0; x2=0.4; y2=0; x3=0.4; y3=0.3; x4=0; y4=0.3; alpha1=0; alpha2=90; alpha3=atan(0.75)*180/pi; k1=Bar2D2Node_Stiffness (E,A,x1,y1,x2,y2,alpha1) k1 =0 0 0 0 0 0 00 0 0 0 0 k2=Bar2D2Node_Stiffness (E,A,x2,y2,x3,y3,alpha2) k2 = 1.0e+007 * 0.0000 0.0000 -0.0000 -0.0000 0.0000 9.8333 -0.0000 -9.8333 -0.0000 -0.0000 0.0000 0.0000 -0.0000 -9.8333 0.0000 9.8333 k3=Bar2D2Node_Stiffness (E,A,x1,y1,x3,y3,alpha3) k3 = 1.0e+007 * 3.77602.8320 -3.7760 -2.8320 2.83202.1240 -2.8320 -2.1240 -3.7760 -2.83203.77602.8320 -2.8320 -2.1240 2.8320 2.1240 k4=Bar2D2Node_Stiffness (E,A,x4,y4,x3,y3,alpha1) k4 = 0 0 0 0 0 0 00 0 0 0 0 （3） 建立整体刚度方程 由于该结构共有 4 个节点，因此，设置结构总的刚度矩阵为 KK(8× 8)，先对 KK 清零， 然后四次调用函数 Bar2D2Node _Assembly 进行刚度矩阵的组装。相关的计算流程如下。 KK=zeros(8,8); KK=Bar2D2Node_Assembly (KK,k1,1,2); KK=Bar2D2Node_Assembly (KK,k2,2,3); KK=Bar2D2Node_Assembly (KK,k3,1,3); KK=Bar2D2Node_Assembly (KK,k4,4,3) KK= 1.0e+008 * 1.1151 0.2832 -0.7375 0 -0.3776 -0.2832 0 0 0.2832 0.2124 0 0 -0.2832 -0.2124 0 0 -0.7375 0 0.7375 0.0000 -0.0000 -0.0000 0 0 0 0 0.0000 0.9833 -0.0000 -0.9833 0 0 -0.3776 -0.2832 -0.0000 -0.0000 1.1151 0.2832 -0.7375 0 -0.2832 -0.2124 -0.0000 -0.9833 0.2832 1.1957 0 0 0 0 0 0 -0.7375 0 0.7375 0 0 0 0 0