弹塑性结构地震反应分析

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