《有限元分析》课程作业.docVIP

《有限元分析》课程作业 任课教师：徐亚兰 学生姓名：陈新杰 学 号班 级：1304012 时 间：2016-01-05 一、问题描述及分析 问题：如图1所示，有一矩形平板，在右侧受到P=10KN/m的分布力，材料常数为：弹性模量；泊松比；板的厚度为t=0.1m；试按平面应力问题利用三角形与矩形单元分别计算各个节点位移及支座反力。 图 平面矩形结构的有限元分析 二、有限元建模及分析 1、基于3节点三角形单元的有限元建模及分析 （1）结构的离散化与编号 如图，各个节点的位移分别沿方向和方向为 图2 方案一使用两个 （2）各单元的刚度矩阵及刚度方程 a.单元的几何和节点描述 单元；同样，将所有节点上的各个力也组成一个列阵，记作，则有 同理对于单元 b.单元的位移场描述 对于单元①，设位移函数 (1-1) 由节点条件处 (1-2) 将式1)代入节点条件式(1-2)中，可求出式（1-1）中 (1-3) (1-4) (1-5) (1-6) (1-7) (1-8) 在式式-8)中 (1-10) 上式中的符号3)表示下标轮换，如同时更换 将单元代入得 式-8)式代入 (1-11) 令 (1-12) 位移函数式(1-11)写成矩阵形式有 对于单元 形状函数矩阵 (1-14) 位移函数 5) c.单元的应变场描述 对于单元 (1-16) 其中几何矩阵 7) 对于单元 (1-18) 其中几何矩阵 9) d.单元的应力场描述 (1-20) 其中弹性系数矩阵D为 21) (1-22) 其中应力函数矩阵S(1)为 23) 应力方程为 24) 对于单元 弹性系数矩阵D(2)为 (1-25) 应力函数矩阵S(2)为 (1-26) 应力方程为 (1-27) .单元的势能表达 是单元刚度矩阵即 其中薄板厚度式代入式 计算得 同理 将两个单元按节点位移所对应的位置进行组装得到总体刚度矩阵为 节点力 系统的势能 (计算结果在下面呈现 (4)边界条件的处理及方程求解 边界条件为。因此，将针对节点2和节点3的位移求解，节点2和节点3对应总体刚度阵KK中的第3行到第6行、第3列到第6列，则需从KK中提出置给k然后生成对应的载荷列阵p再采用高斯消去法进行求解 k=KK(3:6,3:6); p=[500;0;500;0]; u=k\p u =0.0009 0.0001 0.0010 -0.0002 [将列排成了行] 再计算支反力，再代回原整体刚度方程计算出所有的节点力找出对应的支反力 U=[0;0;u;0;0] U = 0 0 0.0009 0.0001 0.0010 -0.0002 0 0 [将列排成了行] P=KK*U P =-500 -176.4706 500 0 500 0 -500 176.4706 [将列排成了行] 所以，节点1的支反力为，节点。 根据已求得的位移和支反力计算系统的势能 A=0.5*U*KK*U-P*U A = -0.4706 (5）结果分析 上述支反力计算结果 2、基于四节点四边形单元的有限元建模及分析 （1）结构的离散化与编号 如图，各个节点的位移分别沿方向和方向为 图 方案 （2）局部坐标系下单元的描述 a.单元的几何和节点描述 采用无量纲坐标 其中则单元四个节点的几何位置为 q；同样将所有节点上的各个力也组成一个列阵记作F则有 b. 设位移函数为 由节点条件处 将位移试函数代入节点条件中求出待定系数和再代入位移函数中整理后得 其中 如以无量纲坐标来表达可写成 带入上式得到形状函数矩阵 写成矩阵形式，有 c.单元的应变场描述 单元应变为 其中几何矩阵为 d. 应力表达式为 其中应力函数矩阵 e.单元的势能表达 以上已将单元的三大基本变量用基于节点的位移列阵来表示将其代入单元势能表达式中有其中为 其中，t为薄板的厚度，，上式的各