ANSYS杆单元作业精选.pdf

一、题目描述 ：对于一个复杂的桁架结构，用有限元法，根据结构 的受力示意图，求出各节点的位移和支座反力。杆件横截面积为 2 3250mm ，弹性模量E为2.0*10^4GPa，杆长度均为3.6m。 一、问题分析 1、结构说明 这是求解平面桁架结构的平衡问题，在结构中有包括水: 平杆和倾斜杆，在这个过程中先将本题中的结构离散化，其中有 个7 11 . 节点和 个单元节点编号如下图所示，单元编号由其两端的节点编 号组成。单元结构尺寸列表如下： 单元 始末端节 倾角（°） 单元 始末端节 倾角（°） 点编号 点编号 1 1-2 60 7 4-5 300 2 1-3 0 8 4-6 0 3 2-3 300 9 5-6 60 4 2-4 0 10 5-7 0 5 3-4 60 11 6-7 300 6 3-5 0 2、解决方法: 解题的基本思路是根据单元的参数和平面桁架的刚度矩阵，写出各单 元的刚度矩阵，然后整合进整体刚度矩阵，引入边界条件后得到整体 刚度矩阵方程，解这个方程就可以得到各节点的位移和受力情况,下 面就按照这种解题思路进行求解。 三、求解过程 1 k12 k13 k23 k24 k34 k35 k45 、求整体坐标下单元刚度矩阵 、 、 、 、 、 、 、 k46 k56 k57 k67 、 、 、 已知L1=L2=L3=L4=L5=L6=L7=L8=L9=L10=L11=0.1853； 节点和单元编号如图： 编写并调用调用函数：PlaneTrussElementStiffness(E,A,L,theta) 和PlaneTrussAssemble(K,k,i,j)求单元刚度矩阵和总体刚度矩阵。 (数据较多就不一一写出来了，后面有详细的介绍) 1、单元刚度矩阵叠加求解整体刚度矩阵K (逐步叠加得到整体刚度矩阵) 2、列出平衡方程 K*U=F得到U=K\F 3、带入边界条件 F=[0;-280000;0;0;0;-210000;0;0;0;-280000;0;0;0;-360000]； U=[0;0;u2;v2;u3;v3;u4;v4;u5;v5;u6;v6;u7;0]； 因为位移边界条件中的有零项: U=K(3:13,3:13)\F(3:13) 解出位移矩阵，并将位移回带入平衡方程求得支反力 F=K*U 四、结果说明: 节点受力及位移表及杆内应力表 节点 1 2 3 4 5 6 7 水平 0 1.5871 0.3840 0.8191 1.1903 -0.0256 1.6127 移(mm) 垂直 0 -180310 -338440 -372440 -359870 -192130 0 移(mm) 水平力 0 0 0 0 0 0 0 (N) 垂直力 233330 0 -210000 0 -280000 0 256670 (N) 单元序 1 2 3 4 5 6 7 8 9 10 11 号 应力 82.9016 41.558 82.9016 82.9016 8.2902 87.0467 8.2902 91.