有限元平面钢架.docVIP

PAGE 1 平面钢架 有限元作业 2017-7-51 2017-7-5 1401s205 — 付小龙 1401s207 — 胡新帅 1401s208 — 贾伟波 目录 TOC \o 1-3 \h \z \u 第一章 问题重述 3 一、题目内容： 3 二、题目要求： 3 第二章 原理说明 4 一、平面刚架的解题思路（步骤）： 4 二、单元刚度矩阵的建立 4 三、总体刚度矩阵 9 四、位移的求解 9 第三章 手工算法求解 11 一、单元离散化 11 二、单元分析 11 三、单元组装 13 四、边界条件引入及组装总体方程 14 五、求解整体刚度方程，计算节点2的位移和转角 14 六、求节点1、3支撑反力 15 七、设定数据，求解结果 15 八、轴力图、弯矩图、剪力图的绘制 16 第四章 matlab编程计算 18 一、流程图 18 二、输入数据 18 三、计算单元刚度矩阵 18 四、建立总体刚度矩阵： 19 五、计算未约束点位移 19 六、计算支反力 19 七、输出数据 19 八、编程 19 第五章 ANSYS分析部分 20 一、参数预处理设置 20 二、建立数学模型 21 三、添加载荷和约束 24 四、分析求解 27 五、绘制图像 29 第六章 结果比较 32 第七章 心得体会 33 胡新帅 33 付小龙 33 贾伟波 33 附录 35 一、matlab原程序 35 二、matlab运行结果 36  第一章 问题重述 一、题目内容： 图示平面钢架结构 图1.1 题目内容 二、题目要求： （1）采用平面梁单元进行有限元法手工求解，要求写出完整的求解步骤，包括： a）离散化：单元编号、节点编号； b）单元分析：单元刚度矩阵，单元节点等效载荷向量； c）单元组装：总体刚度矩阵，总体位移向量，总体节点等效载荷； d）边界条件的引入及总体刚度方程的求解； e）B点的位移，A、C处支撑反力，并绘制该结构的弯矩图、剪力图和轴力图。 （2）编制通用平面钢架分析有限元Matlab程序，并计算盖提，与手工结果进行比较； （3）利用Ansys求解，表格列出B点的位移，A、C处支反力，绘制弯矩图、剪力图和轴力图，并与手算和Matlab程序计算结果比较。 （4）攥写报告，利用A4纸打印； （5）心得体会，并简要说明各成员主要负责完成的工作。  第二章 原理说明 一、平面刚架的解题思路（步骤）： （1）?结构离散化： 将每个杆件作为一个单元，杆件与杆件之间的联结点为节点，并进行编码。 （2）单元刚度矩阵的建立： 具体内容有：建立局部坐标系、单元刚度矩阵的导出、坐标转换、单元刚度矩阵的坐标变换、单元刚度矩阵的特性。 （3）求总体刚度矩阵： 总体刚度矩阵的分块形式、总体刚度矩阵的集成、总体刚度矩阵的特性。 （4）位移的求解： 支撑条件的引入、非节点载荷的处理。 （5）求应变、应力 二、单元刚度矩阵的建立 （1）局部坐标系与整体坐标系： 由于各杆件的取向不同，所以对各杆单元必须先建立各自的局部坐标系。 以节点1(e)为原点，由1(e)到2(e)的方向x轴正向的右手坐标系就叫做该单元的局部坐标系，用英文小写字母x和y表示。如图2-1所示。 图2-1 一维单元（钢架杆单元） 为便于进行整体分析，要求整个结构的所有杆件都采用同一个坐标系，如图所示。该坐标系就称为整体坐标系，用英文大写字母X和Y表示： 在局部坐标系中，杆单元每个节点都有三个位移分量，如图(a)所示。每个节点都有三个力分量，如图（b）所示 图中 ——轴向位移、轴向力 ——横向位移、横向力 ——角位移、弯矩 图中所示各力的方向为正。 每个单元在两端共有六个节点位移和六个节点力分量（2）单元刚度矩阵的导出： 对于刚架杆单元，由于位移和力向量都不再是二维而是六维，所以单元刚度矩阵也不再是2×2阶而是6×6阶，其基本形式为： 求单元刚度矩阵的核心是求各单元刚度系数。下面求单元刚度矩阵。 的第一列对应于只有节点1(e)在轴向方向产生单位位移时的情况，如下图所示。 这时除了外，两节点的其余位移分量都为零。要产生这个变形并能保持该杆单元的平衡，只需要在两个节点处施加两个轴向力 他们分别为刚度系数。很明显，第一列的其它元素均为零。 的第二列对应于只有节点1(e)在横向方向产生单位位移时的情况，如下图所示。 这时除了外，两节点的其余位移分量都为零。它相当于自由端绕度为1、转角为0的一根悬臂梁。要产生这个变形并能保持该梁的平衡，在节点1(e)处要施加横向力 在