matlab 矩阵位移法编程 结构力学,matlab 矩阵位移法编程 结构力学

矩阵位移法编程大作业

(091210211)

一、编制原理

本程序的原理是基于结构力学矩阵位移法原理，以结构结点位移作基本未知量，将要分析的结构拆成已知节点力—结点力位移关系的单跨梁集合，通过强令结构发生待定的基本未知位移，在各个单跨梁受力分析结果的基础上通过保证结构平衡建立位移法的线性方程组，从而求得基本未知量。

二、程序说明

本程序是计算3层11跨框架右侧结点的位移和弯矩的程序，编译过程是按照矩阵位移法的先处理法进行的。首先将结构杆件的交汇点作为结点，共有36个结点和108个位移编号，然后根据梁、柱、斜杆的不同分别建立单元刚度矩阵，然后转换为整体坐标系下的刚度矩阵，然后将所有杆件的单元刚度矩阵整合成为总体刚度矩阵，在进行整合时连续运用for函数，最终形成108阶的总体刚度矩阵。然后通过对荷载的分析自己确定出荷载矩阵，直接写进程序。这样就可以把36个结点的108个位移求得，然后再利用各个单元的单元刚度矩阵和所得的位移求得单元杆件的内力。 离散化编号如下图：

三、算法流程

离散结构化编号 确定梁、柱、斜杆在整体 坐标系下的刚度矩阵 先处理法把各个杆件的单元刚度矩阵整合成总体刚度矩阵 利用杆件单元刚度矩阵和所求位移求内力 单元分析 建立梁、柱、斜杆的局部坐标系下的单元刚度矩阵 确定综合结点荷载矩阵 建立方程，求解位移

输出结果 四、源代码

结束 %结构力学大作业 3层11跨框架矩阵位移法编程 王贝 091210211 h=input('输入单层高h：'); L=input('输入单跨度L：');

EIc=input('输入柱子的抗弯刚度EIc：'); EAc=input('输