有限元分析作业1精选.docxVIP

《有限元分析》课程作业任课教师：徐亚兰学生姓名：林声佳学号级：1304012时间：2016-01-10问题描述及分析如图1所示为一矩形薄平板，在右端部受分布力P=10KN作用，材料常数为：弹性模量E =1×107Pa 、泊松比μ=1/3，板的厚度为t = 0.1m，试按平面应力问题计算各个节点位移及支座反力 y 1m 1②① 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=[u1v1u2v2u3v3]TP=[ 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中有，边界条件的处理及方程求解该问题的位移边界条件是u3=0,v3=0,u4=0,v4=0,将其带入上式中，划去已知节点位移对应的第5行至第8行（列），则所以：节点位移u1=0.0020,v1=-0.0004,u2=0.0018,v2=-0.0000结果分析支反力Rx3=-500.