有限元钢架结构分析~手算+matlab+ansys模拟.docxVIP

有限元大作业 钢架结构分析 选题人: 日期：2016年6月2日 目录： TOC \o 1-5 \h \z \o Current Document 第一章：问题重述 1 一、 题目容： 1 \o Current Document 二、 题目要求： 1 \o Current Document 第二章：有限元法手工求解 2 一、 平面两单元离散化 2 二、 单元分析 2 三、 单元组装 5 四、 边界条件引入及组装总体方程 5 五、 求解整体刚度方程，计算节点2的位移和转角 6 六、 求节点1、3支撑反力 6 七、 设定数据，求解结果 7 八、 绘制轴力图、弯矩图、剪力图 8 \o Current Document 第三章、mat lab编程求解： 9 \o Current Document 一、 总体流程图绘制： 9 \o Current Document 二、 输入数据： 9 \o Current Document 三、 计算单元刚度矩阵： 10 \o Current Document 四、 建立总体刚度矩阵： 10 \o Current Document 五、 计算未约束点位移： 10 \o Current Document 六、 计算支反力： 10 \o Current Document 七、 输出数据： 10 \o Current Document 八、 编程： 10 \o Current Document 第四章有限元求解 11 \o Current Document 一、 预处理 11 \o Current Document 二、 模型建立： 12 \o Current Document 二、 分析计算 14 \o Current Document 三、 求解结果 15 \o Current Document 四、 绘制图像 16 \o Current Document 第五章结果比较 19 \o Current Document 第六章心得体会 19 —、王小灿： 19 二、 明哲： 20 三、 国威 20 \o Current Document 第七章附录 22 * mat lab 程序 22 第一章：问题重述 题目容: 图示平面钢架结构 q=P/lPl/1 q=P/l Pl/1 Wu + + + + + + + + 0.6/EI-20EA 0.6/ 1=100cm ■2/ 图1.1题目容 二、题目要求： （1） 采用平面梁单元进行有限元法手工求解，要求写出完整的求解步骤，包括： a） 离散化：单元编号、节点编号； b） 单元分析：单元刚度矩阵，单元节点等效载荷向量； c） 单元组长：总体刚度矩阵，总体位移向量，总体节点等效载荷； d） 边界条件的引入及总体刚度方程的求解； e） B点的位移，A、C处支撐反力，并绘制该结构的弯矩图、剪力图和轴力图。 （2） 编制通用平面钢架分析有限元Mat lab程序，并计算盖提，与手工结果进行比较； （3） 利用Ansys求解，表格列出B点的位移，A、C处支反力，绘制弯矩图、剪力图和轴力图，并与手算和 Mat lab程序计算结果比较。 （4） 攥写报告，利用A4纸打印； （5） 心得体会，并简要说明各成员主要负责完成的工作。 0 0 0 0 第二章：有限元法手工求解 一、平面两单元离散化 将平面梁离散为两个单元，单元编号分别为①和②，节点号分别为1. 2、3； 如图2-1所示： 广P/2 二、单元分析 首先建立整体坐标系与局部坐标系如图所示; 1、求单元刚度矩阵 对于单元①，求局部坐标系的单元刚度矩阵: EA 0 0 / 0 0 12EZ 6EI 0 \2EI 6EI 卩 I1 I3 I2 6EI 4EI 0 6EI 2EI I1 1 I2 1 0 0 EA 1 0 0 12E/ 6EI 0 \2EI 6EI I3 /2 广 I2 6EI 2EI 0 6EI 4EI I2 1 I2 1 EA 0 EA 由于单元①局部坐标系与整体坐标系的夹角为：0 = 90°,则单元①的局部坐标变换矩阵为: COS(P sin (p 0 0 0 o- ■ 0 1 0 0 0 0 一 sin (p cos (p 0 0 0 0 -1 0 0 0 0 0 T；= 0 0 1 0 0 0 — 0 0 1 0 0 0 0 0 0 cos (p sin (p 0 0 0 0 0 1 0 0 0 0 一 sin (p cos (p 0 0 0 0 一 1 0 0 ■ 0 0 0 0 0 1_ .0 0 0 0 0 1 可以得到在总体坐标系下的单元①的刚度矩阵: \2EI F 0 6E1 l2 \2EI l3 0 6EI