弹塑性结构地震反应分析.ppt

第八章 弹塑性结构地震反应分析 第一节 弹塑性动力分析的一般过程 第二节 串联多自由度体系分析 第三节 平面框架模型 第四节 多维地震波作用下的平扭耦联系统 Date1 第一节 弹塑性动力分析的一般过程 结构弹塑性地震反应问题是地震工程学研究的热点 之一。 一般结构物都在强震中会进入弹塑性变形阶段。 结构的弹塑性反应与线性反应的表现有很大不同： 结构的基本动力特性变化 整体结构的动力反应特征不同 引用弹塑性分析的概念和具体做法，有利于研究结 构地震反应的本质特征，有助于揭示设计结构的最不 利薄弱环节。 Date2 第一节 弹塑性动力分析的一般过程(续） 一、动力方程 二、刚度修正技术 三、一般分析过程 Date3 一、动力方程 结构在多维地震波作用下的一般动力方程为： 结构弹塑性动力分析的基本过程与之相类似： 唯一的变化在于恢复力向量F代替了弹性力向量 KU，这种形式上的替代使我们可以方便地考虑结构 的非线性增量方程 。 Date4 一、动力方程 对单自由度体系，结构在时刻 tj+1的反应可以用tj的反应迭加 一个线形增量: 对多自由度体系，则有 进而得出结构非线性增量方程 : Date5 一、动力方程 用动力分析的逐步积 分法，可以方便地实现结 构弹塑性动力分析计算。 在非线性大变形阶段 ，结构变形可能进入恢复 力下降段，即出现负刚度 。在负刚度条件下各数值 积分方法与正刚度条件有 所不同。 Date6 二、刚度修正技术 结构线性地震反应分析与非线性地震反应分析的主要 差别在于刚度矩阵是否变化。对于弹塑性结构，在每一 步增量反应计算之先，要先行修正刚度矩阵中各元素的 量值，此即刚度修正技术。 修正刚度矩阵的过程实质是重新形成总刚度矩阵的过 程。在这里，区分总刚度矩阵、单元刚度矩阵、刚度系 数、截面抵抗矩等概念十分重要。修正刚度矩阵与应用 恢复力模型的联系途径是通过这些概念转换的。这一途 径可用图6.2加以说明。 Date7 二、刚度修正技术（续） 以常见的三线型刚度退化 型模型介绍刚度修正技术 。骨架曲线包括了开裂点 、屈服点、极限荷载点等 界点。滞回曲线由最大变 形点指向和刚度退化规则 加以规定。在动力计算开 始前要存贮骨架曲线界点 值，在计算中要存贮反向 曾经经历过的变形最大值 和损伤状态值。 Date8 二、刚度修正技术（续） （1）根据变形速度的符号判定变形方向，然后判明本步 变形绝对值是否超过同方向历史最大变形绝对值。当超 过时，则加载点必在骨架曲线上，此时，可将本步累积 变形值与骨架曲线界点变形值相比较。超过界点值时改 变状态标识变量并修正刚度；不超过界点值时，不修正 刚度。而当不超过历史最大变形绝对值时，应进一步判 明相邻时刻内力是否反号，反号时，则修正刚度，否则 不修正刚度； （2）当相邻时刻变形速度值发生变化时，变形反向， 此时，取卸载段退化刚度为本步刚度值。 Date9 二、刚度修正技术（续） Date10 二、刚度修正技术（续） 在刚度修正技术中，还有界点刚度转 换问题，即在前后两时刻刚度发生变化（即 恢复力曲线有转折）时，需将时间步长分割 ，求出刚度发生变化时（即到达恢复力曲线 的转折点）的时刻。在此时刻之前按原刚度 计算，在此时刻之后按改变后的刚度计算。 Date11 三、一般分析过程 弹塑性结构反应分析的思路分为三个基本组 成部分： 数值积分 反应值迭加 刚度修正 一般的分析流程见图6.5。 Date12 三、一般分析过程 Date13 第二节 串联多自由度体系分析 当不考虑结构扭转振动和土 结相互作用时，一般多层结构或高耸结构 可以抽象为一个底部嵌固的串联多自由度 体系。 一、剪切模型 二、弯剪模型 Date14 一、剪切模型 当结构的变形主要表现为集中质量层之间的错动，且这种错动 可视为层间剪切角变位的结果时，则可将结构简化为剪切模型。 一般说来，高宽比不大的多层建筑、强梁弱柱的框架体系等可 以作为剪切结构考虑。 Date15 一、剪切模型（续） 由于不考虑楼板的转角变形，因此，剪切模型的层间 单元刚度矩阵服从以下关系： 其中yi为第I层的位移 ， ， 为剪应力不均匀系数；h为层高；A，J分别为截面 积和惯性矩。 根据Q恢复力关系进行动力分析时，弹性层间刚度 为： 在弹塑性阶段，则有： Date16 二、弯剪模型 高宽比大于4的结构、 强柱弱梁型结构和高耸 结构等，在结构振动时 ，弯曲效应不容忽视。 应采用同时考虑弯曲变 形和剪切变形的弯剪模 型。 Date17 二、弯剪模型 层间单元刚度矩阵服从下述一般关系： 为区别剪切模型，这里以u为水平位移，而为 转角未知量。式中刚度系数的具体形式见公式6.15 。 Date18 二、弯剪模型 弯剪模型区别于剪切模型的根本点在于前者需 考虑楼层处的弯曲转角。引用静态凝聚原理，则 不增加动力方程的自由度数。将总刚度矩阵中与 转角有关的刚度系数并入仅与水平位移有关的刚 度系数项。 自由振动的动平衡方程可以表示为： Date19 二、弯剪模型 从上述方程的后一组n个方程中可以解出： 代入前一组n个方程可得： 其中 称为弯剪模型得等效结构侧移刚度矩阵。 Date20 二、弯剪模型 引用等效侧移刚度矩阵，动力方程如上式 所示，动力自由度数等于集中质量个数。 Date21 第三节 平面框架模型 平面框架模型就是将一般框架结构单独抽出一 榀或经某种假定简化为一榀框架所假定的模型 。 这种模型的弹塑性动力分析反应分析一般以梁 、柱构件作为最小单元进行分析。 典型的构件模型有杆端弹塑性弹簧模型、分割 梁模型、假设变形函数模型等。 Date22 第三节 平面框架模型 一、杆端弹塑性弹簧模型 二、分割梁模型 三、半刚架模型 Date23 一、杆端弹塑性弹簧模型 这种类型的基本思想在于把杆件中的塑性变形全部集 中在杆端，并以杆端等效的弹塑性回转弹簧等价地表示 。弹簧之间的杆件仅发生弹性变形。 模型的基本假定： （1）杆件的弹塑性变形状态可以用图6.8等价地表示 ； （2）杆端塑性转角只与本端弯矩增量有关； （3）采用以单根构件试验为基础的杆端力矩杆端 转角恢复力关系作为基准恢复力曲线，而不考虑各构件 相互联结的影响。 Date24 一、杆端弹塑性弹簧模型 Date25 二、分割梁模型 杆端弹簧模型把沿杆件分布的损伤分别集中于 杆端，并假定杆端塑性转角增量仅与本端弯矩增 量有关。这些假定在反弯点位置偏离构件中点很 远的场合是不适用的。 分割梁模型是把构件分割成若干个沿杆轴线的 假想并列杆件，各杆件仅在杆端相连，而沿杆轴 各点上则有不同的变形。 依据采用的恢复力骨架曲线的不同，分割梁模 型又有双分量模型、三分量模型等类型。 Date26 1、双分量模型 （见图见图 6.11） 克拉夫设想将杆分割为刚度分别为pK和qK的两个平行的 分杆，从而引出双分量模型的概念，图6.11表示这种分 解与合成的全过程。 根据杆端转角的大小，杆件将分别处于两端弹性连接， 一端弹性连接、另一端塑性铰接或两端塑性铰接等不同 状态。 双分量模型的单元刚度矩阵可由各分杆端部连接部分的 变形相容条件和各杆的力平衡条件来确定。（见P178） 双分量模型不考虑刚度退化和反向的最大点指向，卸载 时采用原点加载刚度，反向加载时亦如此。 双分量模型不能表示工程结构构件的刚度退化性质。 Date27 1、双分量模型 （见图见图 6.11） 图c表示第一分杆的恢复力曲线，它表示一根完全弹性杆， 图e表示第二分杆的恢复力曲线，它对应一根弹塑性杆。 Date28 2、三分量模型 （见图见图 6.13） 三分量模型与双分量模型的唯一差别在于它用三 根分杆模拟恢复力模型，图6.13表示其刚度分解过 程和各分杆的恢复力曲线。 构件两端分别有弹性、开裂、屈服等三个不同状 态，因此有九种不同的刚度矩阵。图6.14表示了考 虑刚度退化的恢复力情况。 应当指出，分割梁模型中的单元刚度矩阵是全量 型的，而杆端弹簧模型则为增量型的。 Date29 2、三分量模型 （见图见图 6.13） Date30 三、半刚架模型 采用杆系模型进行弹塑性动力分析中，一 个比较突出的问题就是结构的计算自由度 多，计算工作量大。因此，在上述杆系模 型的基础上，发展了半刚架模型的简化动 力分析方法。 基本的半刚架简化方法有两类： 将各构件特性集成 将各构件特性平均 Date31 三、半刚架模型(续) 如图6.15所示，对于一般的n跨平面刚架结构，依据同一 楼层各节点的转角变形与剪切变形相等的假定，可导出 将平面框架转化为图示半刚架的规则如下： （1）弹性参数 柱弯曲刚度： 梁弯曲刚度： 柱高： 梁长： Date32 三、半刚架模型(续) 如图6.15所示，对于一般的n 跨平面刚架结构，依据同一楼层各节 点的转角变形与剪切变形相等的假定，可导出将平面框架转化为图 示半刚架的规则如下： （2） 恢复力特性参数 柱初始屈服弯矩： 柱极限屈服弯矩： 梁初始屈服弯矩： 梁极限屈服弯矩： 式中引入为反映原框架各层的梁、柱受力不均匀因而屈服 不同时发生所造成的影响，称之为半框架初始屈服折减系数，一般 可取0.70.9。 Date33 三、半刚架模型(续3) 对于图6.16所示的三折线型恢复力骨架曲线，转 化恢复力特性参数尚包括刚度折减系数p1 可取为 ： 柱刚度折减系数： 梁刚度折减系数： Date34 第四节 多维地震波作用下的平扭耦联系统 由于设计要求和随机因素的影响，实际结构都可 视为非对称结构。在多维地震波作用下，非对称结构的 振动一般表现为平移与扭转耦合的振动形式。 根据基本抗侧力构件双向相互作用的强弱，把平 扭耦联振动问题按结构材料类型分为： 弱相互作用模型 强相互作用模型 一、一般平扭耦联系统的动力方程 二、弱相互作用模型 三、强相互作用模型 Date35 一、一般平扭耦联系统的动力方程 一般平扭耦联系统的平面结构形式如 图6.18所示。 对此类结构形式进行平扭耦合振动分析 时，一般沿用楼板平面内无限刚性和平面 外完全柔性的假定。 楼板刚性假定，即不考虑楼板在平面内 的剪切变形与弯曲变形。 楼板平面外完全柔性的假定，则可以不 考虑楼板与框架梁的共同作用。 Date36 一、一般平扭耦联系统的动力方程 一般平扭耦联系统的平面结构形式如图6.18 所示。 对此类结构形式进行平扭耦合振动分析时，一 般沿用楼板平面内无限刚性和平面外完全柔性 的假定。 楼板刚性假定，即不考虑楼板在平面内的剪切 变形与弯曲变形。 楼板平面外完全柔性的假定，则可以不考虑楼 板与框架梁的共同作用。 Date37 一、一般平扭耦联系统的动力方程（续1） Date38 一、一般平扭耦联系统的动力方程（续） 平面子结构在自身主轴方向产生单位层间位移所需力 为k（即层间刚度），则平面子结构i由于位移矢量 而引起的沿子结构轴向的力为: ai为子结构主 轴方向与x轴 的夹角，xi， yi为子结构中 性轴与主轴交 点处的坐标值 ，ri为坐标系 原点到子结构 主轴的距离 Date39 一、一般平扭耦联系统的动力方程（续） 平面子结构的抗弯刚度为： 式中j为单个构件序号；m为子结构中构件个数。 平面子结构的等效抗剪刚度为： 式中第一项为剪力墙平面内的抗剪刚度之和，第二 项为框架平面内等效抗剪刚度之和。推导见PP183 185。 平扭耦联系统的结构动力方程的基本形式为： Date40 二、弱相互作用模型 弱相互作用模型的基本思想是，只考虑一方 向上的强度损失对与其垂直方向的强度及刚度 的影响。例如，砌体结构的抗侧力构件的双向 变形可以处理为平面内剪切变形与平面外弯曲 变形的复合，其非线性关系可分别取图6.21所 示的形式。其中，墙体平面内初始抗剪刚度为 ： 式中，为正应力影响系数；为剪应 力不均匀系数；为弯曲影响系数；为墙体 高宽比；为墙体中点高度平均正应力。 墙体平面外弯曲刚度为： Date41 二、弱相互作用模型（续一） 骨架曲线上的折线刚度由材料试验给出，滞 回曲线可以选用设定模型或复合模型。 （1）墙体初始刚度是平面内抗剪刚度和平面 外抗弯刚度的迭加； （2）当平面内出现屈服进入高度非线性时， 不再考虑平面外刚度的贡献。 Date42 二、内力重分布与变变形集中(续一) 2、变变形集中 塑性铰的形成是非线性变形在构件一小区段内 地集中，与之类似，整体结构层间屈服会造成 非线性变形在该层的集中。 实际结构物，由于工艺、建筑、设计条件等方 面的限制，不可能沿竖向的强度、刚度完全均 匀。在地震中，结构将首