非经典振型分解法求混凝土房屋及其上钢塔的地震响应.pptVIP

非经典振型分解法求混凝土房屋及其上钢塔的地震响应 王依群，王华明 （天津大学建筑工程学院,天津市300072） 引言 钢结构与混凝土结构阻尼比不同，混凝土房屋与其顶上钢塔组成了非比例阻尼结构系统。为避免采用直接积分法耗时多及结构高阶自振周期对积分步长要求较小的限制，本文探讨用非经典振型分解法求解该类结构系统的线弹性地震响应。 1 非经典振型分解法 一般线弹性结构的动力方程为： 2 顶上钢塔的混凝土结构 钢结构塔的阻尼比一般为0.01、混凝土结构的阻尼比一般为0.05，两者的组合构成了非比例阻尼结构系统。以下举一算例表明上述方法的有效性。 按有限元的方法，先建立各子结构自由度编号与结构系统总自由度编号之间的关系，据此关系将各子结构质量、刚度矩阵扩充为阶数为总自由度数的各子结构质量、刚度扩充矩阵[Mi ]、[Ki ]。则具有n子结构的整体结构系统的质量、刚度及阻尼矩阵分别为： 2.1 算例简况 某10层钢筋混凝土框架及其顶上钢塔如图1～3示。钢筋混凝土构件尺寸及材料性质见文[3]、钢塔构件尺寸及材料性质见文[4]。 在结构基础X向施加峰值调整为70gal的EL-Centro波，对该结构系统按二部分不同阻尼比，混凝土框架阻尼比取0.05，钢塔阻尼比取0.01进行计算。为比较起见，还计算了混凝土框架和钢塔阻尼比均取0.05，即比例阻尼的情况。框架楼面采用刚性平面假定后，系统共有792个自由度。采用Newmark逐步积分方法进行了三维弹性动力时程计算。 2.2 算例结果分析 直接积分法和取不同振型反应数目迭加的钢塔顶、框架顶X向位移响应计算结果如图３、４、５示，图中direct代表直接积分、m代表是m阶振型响应迭加的结果。由三图可见，5个振型位移响应迭加结果可与直接积分结果接近。第一阶振型为Y向平动振型，所以在三图示的X向位移响应中，第一阶振型响应所占比重很小。 图3 钢塔顶X向位移响应时程 图4 钢筋混凝土框架顶X向位移响应时程 图5 钢塔顶按0.05比例阻尼X向位移响应时程 由三图比较可见，对于刚度几乎不变且阻尼比也不变的框架部分，仅两个振型反应的迭加就与直接积分的结果相当接近（图5），这是由于二阶以上振型主要是钢塔的变形。而对刚度突变后又逐渐减小且阻尼比不变（图6）、和阻尼比减小（图4）的钢塔，有两点值得注意：1、要有相对多的振型反应参与迭加才能达到与直接积分结果接近的程度，2、较少振型反应迭加的时程曲线波峰、波谷的位置也与直接积分结果差距较大。从图形上看，非比例阻尼情况（图4）相对于比例阻尼情况（图6），这两点更显著，这种差别显然是非比例阻尼特性(振动相位差)造成的。为弥补此影响，对非比例阻尼系统应增加参与响应迭加的振型阶数，由于算例较小，这点显示不明显。 如将此问题当成比例阻尼系统计算（阻尼比0.05的结果见图6），则钢塔顶X向位移最大值为0.366m，而按非比例阻尼系统计算相应值为0.388m；即按比例阻尼计算的结果偏小（此例小了5.7%）。考察钢塔最下面单元弯矩计算精度（前几阶与前30阶振型迭加结果比）随参与迭加振型阶数的变化，发现此例前五阶时达89%，前十阶时达99.8%。此例在pentium4微机Windows XP操作系统下，直接积分位移计算耗时82秒，取5阶振型非经典振型迭加计算耗时11秒，可见效率很高。 3 结语 混凝土楼顶钢塔“鞭梢效应”是由于钢塔刚度和阻尼均小于混凝土房屋的双重因素造成的，忽视二者中任一因素都将造成安全隐患。 高阶振型对顶塔地震响应有较大影响，要较准确地得到顶塔的地震响应，须增加参与迭加的振型阶数；考虑非比例阻尼特性时，参与迭加的振型阶数还应有所增加。 参与迭加的振型阶数应比满足位移精度要求的振型阶数至少多一倍以上时才能达到满意的内力结果，这点与丁幼亮文[5] 按比例阻尼计算结果相同。 震害The toppled crown atop the 95-foot high tower of a 1920 Art Deco building on Hollywood Boulevard. (photo: Ken Fowler) 欢迎请意见 谢谢！ * 对于一实际工程有上千个自由度，直接求解方程要花费大量的计算机时间。对于线弹性结构可采用振型分解法计算 这里 为振型矩阵。 当阻尼满足比例阻尼条件时， 称为经典的振型迭加法。 是满阵，缩减后的方程仍耦连。对此降阶的动力方程组积分，将有限阶振型反应迭加就得到整个结构系统的动力响应[1]，此即非经典的振型分解法。 是对角阵， 非比例阻尼时， 其广义阻尼矩阵为： 将这种不同阻尼材料构成的