结构力学大作业 matlab 矩阵位移.doc

结构力学大作业 二、输出结果的过程（完程序后） 矩阵不限） 2.通过函数A)得到结点位移矩阵函数W,A)得到各单元的的杆端力。的过程详见附录详细程序见附录二） 程序设计过程 上述的输出结果过程，可得，工作量较小，工作中，的单元刚度矩阵以及单元等效结点荷载都仅仅是码字工作，唯一难度的是定位向量的，将的单元刚度矩阵以及单元等效结点荷载成刚度矩阵整体等效荷载。 思路：作用是定位，试想定位向量得到行列变换，得到刚度矩阵假设共有位移编号，则建立=(N+6)×(N+6)的单位矩阵。kk=)×(N+6)的空矩阵为rr(1,2,5,13,14,0)系下单元刚度矩阵为ke其中kk左上角该刚度矩阵ke 目标：矩阵列变换，将的右下角变成刚度矩阵即得到相应的行变换矩阵RR)以及列变换矩阵RR’)。过程见附件 rr.m） 计算结果求解器对比）刚度矩阵 结点等效荷载 得出： 与求解器结果的比较）过程： 结点,2,0.6,3 结点,3,1,5 结点,4,4,5 结点,5,4,3 结点,6,4,0 单元,1,2,1,1,1,1,1,0 单元,2,3,1,1,0,1,1,1 单元,3,4,1,1,1,1,1,1 单元,4,5,1,1,1,1,1,1 单元,5,6,1,1,1,1,1,1 单元,2,5,1,1,0,1,1,1 单元,2,6,1,1,0,1,1,1 单元荷载,3,3,1,0,1,90 单元荷载2,4,5,3,0.5,0,1,90 单元荷载,6,1,4,1/2,90 结点支承,1,6,0,0,0,0 结点支承,6,6,0,0,0,0 单元材料性质,1,7,0.5,1,0,0,-1 杆端 1 杆端 2 ---------------------------------------- ------------------------------------------ 单元码 u -水平位移 v -竖直位移 ?-转角 u -水平位移 v -竖直位移 ?-转角 ----------------------------------------------------------------------------------------------- 1 0 0 0 -5 -17.3064174 1 2 -5 -17.3064174 -0 -6 -24.7602310 -0 3 -6 -24.7602310 -0 -9 -24.3680038 1 4 -9 -24.3680038 1 -6 -19.3457527 2 5 -6 -19.3457527 2 0 0 0 6 -5 -17.3064174 -3 -6 -19.3457527 2 7 -5 -17.3064174 5 0 0 0即，整体得到 结构力学求解器 MATLAB计算 易看出，精度低了一些，计算结果完全一致。 、 MATLAB计算的 下面是中结果 弯矩图 剪力图 轴力图 六、总结结果与求解器完全一致，但的说的程序远远不及求解器编写的程序单元的定位需要手动输入，而求解器是自动完成的，方便了很多但，在求解杆端力，结点位移上，应该是比较接近求解器的，的数没有限制程序，因素以及时间安排不合理，导致太多的时间在上面，可以臻于完美的功能没能实现。编写程序，也让我极大的矩阵位移法，也让我土木人也是应该学程序的。 附录一 1.输入原始数据矩阵ac 2.运行函数k(ac)，得到结点位移de 3.运行函数li(de,ac)，得到