有限元分析作业1.docxVIP

《有限元分析》课程作业 任课教师：徐亚兰 学生姓名：林声佳 学号班级：1304012 时间：2016-01-10 问题描述及分析 如图1所示为一矩形薄平板，在右端部受分布力P=10KN作用，材料常数为：弹性模量E =1×107Pa 、泊松比μ=1/3，板的厚度为t = 0.1m，试按平面应力问题计算各个节点位移及支座反力 y 1m 1 1m ② ① 3 ② ① P P 1m 4 2 x (a)问题描述 (b)有限元分析模型 图1 右端部分受均布力作用的平面问题 二、有限元建模及分析 1、基于3节点三角形单元的有限元建模及分析 （1）结构的离散化与编号 对该结构进行离散，单位编号及节点编号如图（1-b）所示，即有两个3节点三角形单元。载荷F按静力等效原则向节点1，节点2移置等效。 节点位移列阵 q=[u1v1u2v2u3v3u4v4]T 节点外载荷列阵F=[ 0 0 0 0 0 0]T 约束的支反力列阵 R=F=[0 0 0 0 Rx3 Ry3 Rx4 Ry4]T 总的节点载荷列阵 P=F+R=F=[ 0 0 Rx3 Ry3 Rx4 Ry4]T 其中的Rx3、Ry3、Rx4、 Ry4分别为节点3和节点4的两个方向的支反力 （2）各单元的刚度矩阵及刚度方程 a.单元的几何和节点描述 图2 如图2所示，一个单元体有6个节点位移自由度（DOF），将所有节点上的位移组成列阵q；同样，将所有节点上个力也组成列阵P，那么 q=[u1v1u2v2u3 v3]T P=[ Px1 Py1 Px2 Py2 Px3 Py3]T 当两个单元取图a示中的局部编码（i,j,m)时，其各单元的位移场，应力场，应变场，势能，刚度矩阵完全相同。 b.单元的位移场描述 就如图 2所示的平面 3 节点三角形单元，由于有 3 个节点，每一个节点有两个位移，因此共有 6 个节点位移，考虑到简单性、完备性、连续性及待定数的唯一确定性原则，分别选取单元中各个方向的位移模式为 N(x,y)= Ni=(a+bix+ciy),i=1,2,3 其中a=x2y3-x3y2,b=y2-y3,c=-x2+x3 上式中的符号（1，2，3）表示下标轮换，如 1→2，2→3，3→1 同时更换。 u(x,y)=N（x,y)·q c.单元的应力场描述 （x,y)= = =[?]u 其中[?]为几何方程的算子矩阵(operator matrix)，即 [?]= （x,y)=[?]·N(x,y)·q=B(x,y)·q B=[B1 B2 B3] 其中 Bi=·,i=1,2,3 d.单元的应变场描述 由弹性力学中平面问题的物理方程，将其写成矩阵形式 其中平面应力问题的弹性系数矩阵 D 为 D= 若为平面应变问题，则将上式中的系数(E,μ)换成平面应变问题的系数。 e.单元的势能表达 单元的势能的表达以上已将单元的三大基本变量(u, ,ε σ) 用基于节点位移列阵 来进行表达q来进行表达，将其带入单位势能的表达式，有 其中K是单元刚度矩阵，即 t为平面问题的厚度，这时B为常系数矩阵，因此上式可以写成 其中各个矩阵为 其中 f.单元刚度阵及刚度方程 将单元的势能式对节点位移取一阶极值,可得到单元的刚度方程q （3）整体刚度阵及刚度方程 a.单元贡献矩阵 各个单元的描述当两个单元取图示中的局部编码(i,j,m)时，其单元刚度矩阵完全相同，即 b.整体刚度阵及刚度方程 按单元的位移自由度所对应的位置进行组装可以得到整体刚度矩阵，该组装过程可以写成 具体写出单元刚度矩阵的各个子块在总刚度矩阵的对应位置如下 带入整体刚度Kq=p中有， 边界条件的