基于单元一向量技术的结构矩阵分析研究.pdfVIP

基于单元定位向2技术的结构矩阵分析 赵超燮 赵东松 (北京工业大学建工学院)(北京市市政工程设计研究总院) [提要]本文总结单元定位向量的使用经验 首次全面系统地沦述采用单元定位向量进行结构矩阵 分析的全过程。并为忽略轴向变形、主从位移关系、各种中间连接结点以及边界条件的处理指出统一、 简捷的途径 一、前 官 结构矩阵分析是采用矩阵和单元定位向量技术分析杆件结构受力状态的方法，因此 也有人称为杆系结构有限元法。结构矩阵分析法可分为矩阵位移法和矩阵力法，这里只讨 论前者。 杆系结构有限元法的要点是先把整个结构分拆为许多单元，根据这些单元在有限个结 点处的平衡条件或势能原理，装配成结构刚度方程，求解这组联立方程得到结点位移。这 样一分一合，就把复杂的结构计算转化为简单的单元分析和整体求解 从而可以利用计算 机按照单元定位向量 自动地进行计算。 采用矩阵表达整个结构受力分析过程是十分简洁的，按照矩阵运算规则进行计算可得 到准确的结果。然而为了提高计算效率，在编写程序时需要选择适当的算法，这样就在用 矩阵表达的数学模型与计算机程序之间出现 一定的距离和障碍。为了跨越这个距离，在编 程中就不得不加以研究。如果在结构矩阵分析中采用单元定位向量技术，则装配结构刚度 矩阵、结构结点力列阵以及求取单元结点位移列阵，均可直接采用单元定位向v所描述的 装配逻辑关系式来表达。同时可以直接按照装配逻辑关系式编写程序段的语句，这样就填 补了从数学模型到程序编写之间的缝隙。此外，对于一些特殊的问题，如忽略轴向变形、 主从位移关系、各种中间连接结点以及边界条件的处理等，利用单元定位向量来处理具有 方法统一、简单和便于程序编写等优点。 二、求解步砚 卜面以平面刚架为例来说明矩阵位移法的求解步骤: (1)有关结点、单元和荷载信息 根据结构计算简图确定结点编号、结点坐标和各结 点的约束代码;单元的连接结点、单元的材料性质分组号和材性值 (EA,El)以及荷载工 况、给定位移和弹性支座等信息; (2)形成在局部坐标系 卜的单元刚度矩阵el，并变换到裕体坐标系 卜的单元刚度矩阵 r: (3)利用单元定位向量把单元刚度矩阵K`装配成结构刚度矩阵K ，348 ， c4)利用单元定位向量把结点荷载或单元等效结点力列阵 厂装配成结构结点力列M (5)求解结构刚度方程Kd=F,得到结构结点位移列阵d; (6}利用单元刚度矩阵和单元结点位移列阵计算杆端力和支座反力 (7)校核 三、单元定位向t 1.坐标系和结点位移未知量 约定采用右手直角坐标系: (L)局部坐标系— 用于单元分析和计算单元内力; cz整体坐标系— 用于装配结构刚度矩阵和计算结点位移和支座反力 在空间任一个点的可能位移为(o,V,We.,a.,B.)ou为沿 :坐标轴方向的线位移 护为 沿r坐标轴方向的线位移;，为沿Z坐标轴方向的线位移 a为绕x坐标轴方向的角位移 a,为绕Y坐标轴方向的角位移;9为绕2坐标轴方向的角位移 常见结构的结点位移未知量列在表1中 班1常见结构的结点位移未知， 结构名称 u v w Bx 氏 a- N 架 (书s架 召 J 连 续 梁 了 J 甲而 刚 架 J J J 交 义 梁 系 J J J 空 问 析 架 了 J J 空 A(I刚 契 J J 丫 J