下雨天的天.docVIP

黄河科技学院 结构力学论文 结构力学中矩阵位移法的 院系：工学院建筑系 班级：级土木二班 学号： 姓名： 结构力学中矩阵位移法的 \ 摘　要 :根据矩阵位移法的原理 ,讨论了其计机程序实现中需要注意的几个问题 ,旨在为广大师生提供一些关于如何利用程序解决求解平面刚架、桁架、连续梁的位移、内力等问题的帮助。 关键词 :矩阵位移法 ;单元刚度 ;结构刚度。 Programming Implementation of Matrix Displacement Method Abstract : This paper discusses items on comp uter p rogramming implementation of matrix displacement method based on its. 结构矩阵分析方法是以传统结构力学理论为基础、以矩阵作为数学表述形式、以电子计算机作为计算手段大规模的计算方法。矩阵位移法是结构力学中用来求解杆系结构位移、内力的重要方法 ,它的理论基础是位移法。 矩阵位移法思路 取结构结点位移作为基本未知量 ,将待分析的结构离散成为单个杆件 ,然后通过对各单元的分析和组装 ,建立结构刚度方程 ,求得结点位移和单元杆端力。它实质上是位移法的一种拓展 ,就是将矩阵运算引入到位移法中 ,从而利用计算机程序进行数值计算。由于计算机的强大运算能力 ,使得原本用手工难以完成的较大数目未知量求解 ,变的易于实现。正是矩阵位移法的存在 ,才使得结构力学更具活力 ,能够解决更为复杂的问题。 矩阵位移法之所以具备如此强大的功能 ,需归功于其程序化的结果。这里我们介绍一下矩阵位移法计算机程序实现的思路 ,学习它的使用方法。 在使用矩阵位移法解题时 ,我们首先要进行结构的离散化 ,这一步是需要手工完成的。对用结点将结构进行划分所得到的单元集合体 ,按一定顺序对结点、单元分别加以编号 (如图)。接着用数字描述结点坐标、单元材料与截面特性以及支承信息和荷载信息等等 ,为矩阵位移法的分析、计算 包括结点个数、坐标位置、单元数、单元码单元的截面性质和单元的弹性模量。对于刚架的铰结点 ,由于相对转动位置不确定 , 需要在同一位置编多个号 ,如图 中的 3、4结点。一般按结点 和单元编号的顺序依次输入。其中需要注意的是 ,在进行连续梁和刚架分析中不考虑轴向变形时 ,往往不对程序做改动 ,而是在结构数据输入中将相应单元的设为较大值 ,可扩大 1000倍 ;对于桁架分析 ,没有弯曲变形存在可将设为 0。约束数、约束所在结点和方向。对于桁架的各结点 ,即使是自由端 ,也可视其转角位移为 刚度矩阵 单元刚度矩阵可分为在局部 坐标下 ,和在整体坐标下两种。 其中各个局部坐标下的单元 刚度矩阵具有统一的基 本形式 (这里假设杆端弯矩 顺时针为正) : 称为局部坐标系中的单元刚度矩阵（简称单刚）。 的行数等于杆端力向量的分量数, 列数等于杆端位移向量的分量数，的每一个元素称为单元刚度系数，其表示了一个力。 总刚度矩阵反映的是整个结构的刚度 ,根据引入边界条件的先后 ,形成总刚度矩阵的方法分为后处理法和先处理法。在形成结构刚度矩阵后引入位移条件的方法称为后处理法 ;在形成结构刚度矩阵之前已经引入位移条件的方法称为先处理法。在程序实现时一般选用后处理法 ,采取对号入座 ,由单元刚度矩阵形成总刚度矩阵。 单元刚度矩阵 单一位移时的单元杆端力 平面桁架单元单元刚度矩阵为: 连续梁的整体刚度矩阵 若把连续梁两支座间的一跨取单元，杆端位移条件为： 杆端位移向量与单元杆端力向量为: 单元刚度方程为 单元刚度矩阵的坐标变换: 整体分析时必须建立一个统一的坐标系，称为整体坐标系，其作用是把各单元上不同方向的量值统一到整体坐标系方向上来。整体坐标系中，单元杆端位移向量记为{δe },单元杆端力向量记为{Fe} 单元坐标转换矩阵: 局部坐标系 与整体坐标 系为xoy的夹角α以x轴 逆时针转到与局部坐标系 为正。 端点i 处的杆端力分量，有下列转换关系： j 端点杆端力转换关系 联合起来写成矩阵形式 简记为 [T]称为单元坐标转换矩阵, [T]是一正交矩阵。 [I]为与[T ]同阶的单位矩阵。 例：试求图示刚架的原始刚度矩阵。 已知各杆