有限元分析期末作业(最终版).doc

有限元分析课程作业任课教师：徐亚兰学生姓名：学号：班级：时间：2017-01大作业要求:双面打印,手写,不准用铅笔。程序附后，双面打印,需包含子程序代码，主程序和子程序每一步都需要有注释。主程序中不要写运行结果，运行结果在“计算结果”中给出。除程序外不得另附页。程序要分别写清楚是两单元还是四单元的。满分 100 分,占总分 70%评分细则：二、有限元建模及分析1、基于两个三节点三角形单元的有限元建模及分析（20 分）（1）结构的离散化与编号（4 分）（2）各单元的刚度矩阵及刚度方程（5 分）（3）整体刚度阵及刚度方程（7 分）（4）边界条件的处理及方程求解（2 分）（5）结果及分析（2 分）2、基于四个三节点三角形单元的有限元建模及分析（20 分）（1）结构的离散化与编号（4 分）（2）各单元的刚度矩阵及刚度方程（5 分）（3）整体刚度阵及刚度方程（7 分）（4）边界条件的处理及方程求解（2 分）（5）结果及分析（2 分）3.两种方案的比较与分析（10 分）三、基于 MATLAB 的编程实现 1. 基于两个三节点三角形单元的有限元建模及分析（20 分）（1）程序编写说明（5 分）（2）程序清单（10 分）（3）计算结果（5 分） 2. 基于四个三节点三角形单元的有限元建模及分析（20 分）（1）程序编写说明（5 分）（2）程序清单（10 分）（3）计算结果（5 分）四、结束语（10 分）打分表：满分得分基于两个三角形单元的有限元建模及分析20基于四个三角形单元的有限元建模及分析20两种方案的比较与分析10基于两个三角形单元的有限元编程实现20基于四个三角形单元的有限元编程实现20结束语10总分100一、问题描述及分析问题：如下图所示矩形平板，不计体力，受如图三角形分布力作用，按平面应力问题，分别利用两个平面三节点三角形单元和四个平面三节点三角形单元对该结构进行有限元分析（求节点位移，应力场，支反力）。其中材料常数为：弹性模量 E = 2 1011 Pa, 泊松比 m = 0.3 ,板厚 t = 0.1m ，板长1.m （ 为每位学生学号的后两位）， 板宽1m 。p0 = 10000 N / m二、有限元建模及分析1、基于两个 3 节点三角形单元的有限元建模及分析（1）结构的离散化与编号(含计算等效节点载荷)（2）各单元的刚度矩阵及刚度方程a.单元的几何和节点描述b.单元的位移场描述c.单元的应变场描述d.单元的应力场描述e.单元的势能表达f.单元刚度阵及刚度方程（3）整体刚度方程(需要具体写出单元刚度矩阵各子块在总体刚度矩阵的对应位置)（4）边界条件的处理及方程求解（5）结果及分析2、基于四个 3 节点三角形单元的有限元建模及分析（1）结构的离散化与编号（2）局部坐标系下单元的描述a. 单元的几何和节点描述(略)b. 单元的位移场描述(略)c. 单元的应变场描述(略)d. 单元的应力场描述(略)e. 单元的势能表达(略)f. 单元刚度阵及刚度方程（3）整体刚度方程(需要具体写出单元刚度矩阵各子块在总体刚度矩阵的对应位置)（4）边界条件的处理及方程求解（5）结果及分析3.两种方案的比较与分析三、基于 MATLAB 的编程实现 1. 基于两个 3 节点三角形单元的有限元编程实现（1）程序编写说明（2）