多自由度体系的地震反应分析.pptVIP

2003-8-9 多自由度体系的地震反应分析 计算简图 对于多层房屋，可把质量集中在每个楼层处。如多层框架、多层砖房。 计算简图——续 多跨不等高单层厂房，可以把质量集中到每个柱顶或屋盖标高处。 计算简图——续 烟囱，根据计算要求可分成若干段，每段集中成一个质点。 运动方程 取每个集中质量作为隔离体 运动方程——续 根据质点的平衡条件： 对于所有质点（i=1，2，…，n），可以得到与式（2）相似的一组方程。该方程组可简化为矩阵的形式： 多自由度体系的自由振动 频率方程 为使计算简化，按无阻尼自由振动方程来计算结构体系的自振频率及主振型。 无阻尼自由振动方程： 设方程的解为： 将式（b）对时间求二次导数得： 将式（b）和式（c）代入式（a）得： 显然，{X}不能为零，否则体系无振动。 为使得方程有非零解，则方程的系数行列式必须等于零： 即为频率方程，解之可得体系的自振频率： 2. 主振型 对于不同的频率?j ，可以写出n个主振型向量方程： 由上式可求得体系的n个主振型。 振型的正交性 求频率举例 计算图示二层框架结构的自振频率和振型。 解答 计算刚度系数 由上图可以写出刚度系数： 质量矩阵为 由上式得 解上式一元二次方程得自振频率 由振型向量方程求得振型 主振型图 第1振型、第2振型分别如下图示 自振频率和振型的实用计算 已知： [M]，[K] 求： 频率、振型 方法： 广义雅可比法、矩阵迭代法、能量法、等效质量法等 1.广义雅可比法 本方法是计算机程序，不适合手算 已知[M],[K]，利用程序可以求出全部的 1.广义雅可比法 广义雅可比法很成熟，网上可以找到多种语言版本 Matlab软件中的eig函数也是广义雅可比法 2.矩阵迭代法 较繁琐，适合求第一振型和频率 矩阵迭代法是采用逐步逼近的计算方法来确定结构的频率与振型。 由结构力学知，刚度矩阵与柔度矩阵互逆，即 则式 改写成： 展开后得： 迭代法步骤 先假定一个振型并代入上式等号的右边，计算后得频率?2和主振型第一次近似值。 再将第一次近似值代入上式的右边，计算后得频率?2和主振型的第二次近似值。 如此计算，直至前后两次结果接近为止。 当一个振型求得后，则利用振型正交性质，求出更高阶的频率和振型。 迭代法举例 某三层框架结构图，假定横梁的刚度为无限大。各参数如图示。用矩阵迭代法求结构的频率和振型。 已知质量为： m1=2561t，m2=2545t，m3=559t 层间刚度为： k1=5.43×105kN/m， k2=9.03×105kN/m， k3=8.23×105kN/m 柔度系数：单位力作用下的位移。 解答 求柔度系数 第1振型 设第1振型的第一次近似值为： 上式运算后得 再将此第三次近似值代入得 由（a）与（b）比较可以看出，第三次近似值很接近第二次近似值，因此取第一振型为： 第二振型 对于第二振型 将上式展开得： 现在对（e）式进行迭代，假设 再由（d）式得 第三振型 根据主振型的正交性，则有 将上两式展开，得： 求第三频率， 最后得相应的第三振型： 矩阵迭代法讨论 采用矩阵迭代法求频率和振型时，因求高频及振型时需取已求得的较低振型，故计算误差随振型提高而增加。但分析中，一般只需采用前几个振型，故误差积累对结构反应分析影响不大。 采用矩阵迭代法，当质点较多时，计算太繁。 3.能 量 法 能量法（Reyleigh法），求解第一频率 要点： 设法找到一个近似的振型{X}，越接近越好 结构的基本振型可以近似取为当重力荷载水平作用于质点上时的结构弹性曲线。 能量法讨论 因采用了近似的振型曲线，故基本频率也是近似的。 若要提高计算的精度，应提高振型的精度。可采用迭代法进行修正。 能量法算例 已知质量为： m1=2561t m2=2545t m3=559t 层间刚度为： k1=5.43×105kN/m， k2=9.03×105kN/m， k3=8.23×105kN/m 解答（采用能量法） 结构在重力荷载Gi水平作用下的弹性曲线如图。 层间相对位移 各质点的水平位移 相应的基本振型 为了提高精度，进行修正。 由惯性力引起的层间相对位移 各层位移为 基本频率为 基本振型为 4.等效质量法 也称折算质量法、转移质量法 在求多自由度