有限元作业：三角形单元求解.docVIP

有限元作业：三角形单元求解 有限元作业：三角形单元求解 PAGE 有限元作业：三角形单元求解 《有限元作业》 年 级 2015级 学 院 机电工程学院 专业名称 班级学号 学生姓名 2016年05月 如下图所示为一受集中力P作用的结构，弹性模量E为常量，泊松比V=1/6,厚度为I=1。按平面应力问题计算，运用有限元方法，分别采用三角形及四边形单元求解，求节点位移及单元应力（要求三角形单元数量不少于4个，四边形单元不少于2个） 图(一) 图（二）三角形单元求解 图（三）四边形单元求解 如图划分三角形单元，工分成四个分别为???④ 如图分别进行编号1、2、3、4、5、6，并建立坐标系 编程进行求解，得出结果，其中假设力P=2000N 调用Triangle2D3Node_Stiffness函数，求出单元刚度矩阵 k1 = +06 * 0 0 0 0 0 0 0 0 k2 = +06 * 0 0 0 0 0 0 0 0 k3 = +06 * 0 0 0 0 0 0 0 0 k4 = +06 * 0 0 0 0 0 0 0 0 调用Triangle2D3Node_Assembly函数，求出总体刚度矩阵 求出的节点位移 U = 0 0 0 0 调用Triangle2D3Node_Stress函数，求出应力,S1、S2、S3、中求出的分别为Sx,Sy,Sxy S1 = +03 * S2 = +03 * S3 = +03 * S4 = +03 * 二、 （1）如图划分四边形单元，工分成四个分别为?? （2）如图分别进行编号1、2、3、4、5、6，并建立坐标系 （3）编程进行求解，得出结果，其中假设力P=2000N 调用 Quad2D4Node_Stiffness函数，求出单元刚度矩阵 调用Quad2D4Node_Assembly函数，求出求出总体刚度矩阵 求出节点位移 U = 0 0 0 0 调用Quad2D4Node_Stress函数，求出单元应力中的的S1、S2、S3分别为Sx,Sy,Sxy应力分量 S1 = +03 * S2 = +07 * 程序附录 一、 1、三角形单元总程序： E=1e7; NU=1/6; t=1; ID=1; %调用Triangle2D3Node_Stiffness函数，求出单元刚度矩阵 k1=Triangle2D3Node_Stiffness(E,NU,t,0,1,