自编程序计算空间框架周期，与sap2000对比

本帖最后由 silyvin 于 2013-6-6 12:20 编辑


跨度4m（局部8m），首层高度6m，其他4m
弹性模量取2×10^11 N/m^2
密度7850kg/m^3，所有杆件的截面积为2m^2（抱歉，为了方便）

分析为侧刚模型，引入刚性楼板假定，位移自由度为3×节点数+3×楼层数
动力分析时，将刚度矩阵静力凝聚，自由度减为3×楼层数，
质量矩阵为对角阵，阶数为3×楼层数，对于层来说，前两个为楼层质量，包括梁、柱质量（目前没有板、隔墙、恒活等），第三个为楼层转动惯量， J=mr^2，这里比较纠结，我是将梁柱取一半的自重加到节点，然后算出质心，然后计算各节点到质心的距离 然后平方再乘以节点的质量（包括与该节点相连杆件的一半质量）

这是我的 这是sap2000的

sap2000 的模态好像超过9，难道用的是总刚模型，sap2000设置了三层节点的constraints约束其z轴角位移为同层牵连
质量源定义为第一项，来自对象和附加质量，反正现在没有恒活之类的荷载

对了，xoy平面唯一的那根斜杆，没有计入质量矩阵，不知道是不是这个问题