静力弹塑性分析

1、静力弹塑性分析(Pushover分析)简介Pushover分析是考虑构件的材料非线性特点，分析构件进入弹塑性状态直至到达极限状态时结构响应的方法。Pushover分析是最近在地震研究及耐震设计中经常采用的基于性能的耐震设计(Performance-BasedSeismicDesign,PBSD)方法中最具代表性的分析方法。所谓基于性能的耐震设计就是由用户及设计人员设定结构的目标性能(targetperformance),并使结构设计能满足该目标性能的方法。Pushover分析前要经过一般设计方法先进行耐震设计使结构满足小震不坏、中震可修的规范要求，然后再通过pushover分析评价结构在大震作

2、用下是否能满足预先设定的目标性能。计算等效地震静力荷载一般采用如图2.24所示的方法。该方法是通过反应修正系数(R)将设计荷载降低并使结构能承受该荷载的方法。在这里使用反应修正系数的原因是为了考虑结构进入弹塑性阶段时吸收地震能量的能力，即考虑结构具有的延性使结构超过弹性极限后还可以承受较大的塑性变形，所以设计时的地震作用就可以比对应的弹性结构折减很多，设计将会更经济。目前我国的抗震规范中的反应谱分析方法中的小震影响系数曲线就是反应了这种设计思想。这样的设计方法可以说是基于荷载的设计(force-baseddesign)方法。一般来说结构刚度越大采用的修正系数R越大，一般在110之间。但是这种基

3、于荷载与抗力的比较进行的设计无法预测结构实际的地震响应，也无法从各构件的抗力推测出整体结构的耐震能力，设计人员在设计完成后对结构的耐震性能的把握也是模糊的。基于性能的耐震设计中可由开发商或设计人员预先设定目标性能，即在预想的地震作用下事先设定结构的破坏程度或者耗能能力，并使结构设计满足该性能目标。结构的耗能能力与结构的变形能力相关，所以要预测到结构的变形发展情况。所以基于性能的耐震设计经常通过评价结构的变形来实现，所以也可称为基于位移的设计(displacement-baseddesign)。rasnsesaDDisplacement图2.24基于荷载的设计方法中地震作用的计算Pushover

4、分析是评价结构的变形性能的方法之一，分析后会得到如图2.25所示的荷载一位移能力谱曲线。另外，根据结构耗能情况会得到弹塑性需求谱曲线。两个曲线的交点就是针对该地震作用结构所能发挥的最大内力以及最大位移点。当该交点在目标性能范围内，则表示该结构设计满足了目标性能要求。5%SmoothElasticDesignSpectanuThrFe8AlhtceDS分析方法结构设计必须满足规范的一系列规定和要求，在完成满足规范要求的设计之后，结构的目标性能具体控制在哪个水准上，则由建筑物的使用者和设计者决定。为了评价结构性能需要进行结构分析，基于性能的耐震设计方法一般有下列四种。线性静力分析方法(Linear

5、StaticProcedure,LSP)线性动力分析方法(LinearDynamicProcedure,LDP)非线性静力分析方法(NonlinearStaticProcedure,NSP)非线性动力分析方法(NonlinearDynamicProcedure,NDP)MIDAS/Gen中提供了上述四种分析方法，其中Pushover分析属于非线性静力分析方法。Pushover分析又称为静力弹塑性分析,是评价结构进入弹塑性状态后的结构极限状态和稳定状态的有效而简捷的方法。该方法主要适用于低频结构影响较大的结构中Pushover分析中可以考虑材料和几何非线性，材料非线性特性是通过定义构件截面的荷载

6、一位移的非线性特性实现的。Pushover分析是通过逐渐加大预先设定的荷载直到最大性能控制点位置，获得荷载一位移能力曲线(capacitycurve)。多自由度的荷载一位移关系转换为使用单自由度体系的加速度-位移方式表现的能力谱(capacityspectrum),地震作用的响应谱转换为用ADRS(Acceleration-DisplacementResponseSpectrum)方式表现的需求谱(demandspectrum)o通过比较两个谱曲线，评价结构在弹塑性状态下的最大需求内力和变形能力，通过与目标性能的比较，决定结构的性能水平(performancelevel)。在MIDAS/Gen

7、中使用ATC-40(1996)和FEMA-273(1997)中提供的能力谱法(CapacitySpectrumMethod,CSM)评价结构的耐震性能。能力谱法(CSM羽原理如图2.26所示。(a)计算结构物的能力曲线(capacitycurve)和能力谱(capacityspectrum)CapacitySpectrumSDOFSystem(b)计算需求谱(demandspectrum)(c)评价性能点(performancepoint)图2.26能力谱法(CapacitySpectrumMethod,CSM)的原理Pushover分析是为了评价结构所拥有的耐震性能，具前提条件是已经完成了初

8、步的分析和设计，即对于混凝土结构必须已经完成了配筋设计。Pushover分析的优点如下：可以评价结构进入弹塑性阶段的响应以及所拥有的抵抗能力可以掌握结构的耗能能力和位移需求可以掌握各构件屈服的顺序对确定需要维修和加固的构件提供计算依据分析中适用的单元类型MIDAS/Gen中Pushover分析中适用的单元类型有二维梁单元(2-dimensionalbeamelement)、三维梁一柱单元(3-dimensionalbeam-columnelement)、三维墙单元(3-dimensionalwallelement)、桁架单元(trusselement)。各单元的特性如下。二维梁单元和三维梁-柱

9、单元梁单元和梁-柱单元采用的模型如图2.27所示，其位移和荷载如下，其中适用于梁单元时无轴力项(1.a)PT=FxrMx1,Fy1,IM*Fz1,Mz,Fx2,Mx2,Fy2,My2,Fz2,Mz2(1.b)u=ux1，71x1，vy1，?y1，''z1，71z1，Ux2，入2，vy2，71y2，',z2，入2Shearspring图2.27二维梁单元和三维梁-柱单元模型三维墙单元模型如图2.28所示墙单元模型由中间的线单元，上下两端的刚性杆构成。中间的线单元与三维梁一柱单元相同，刚性杆在xz平面内做Axialspring图2.28墙单元的节点力和节点位移刚体运动。桁架

10、单元模型如图2.29所示，桁架单元采用轴向(x方向)的弹簧模型。非线性弹簧的特性在各单元模型中表现的弹簧并非表示弹簧的存在，而是表现分析的方法，即在弹簧位置将发生塑性变形。弹簧具有的特性如下。梁单元模型的弹簧特性用荷载一位移、轴力-单向弯矩-位移角、剪力-剪切变形、扭矩-扭转角等关系来表现；柱以及墙体单元模型的弹簧特性用荷载一位移轴力一双向弯矩一位移角、剪力一剪切变形、扭矩-扭转角等关系来表现；桁架单元模型的弹簧特性用荷载一位移关系来表现。单元的变形可用下面的各式来表现。弯曲变形节点上构件的变形角可用下列三项之和来表现。r-erP.rs(2)在此，小、e、&分别为弹性弯曲变形角、塑性弯

11、曲变形角、因剪切产生的弯曲变形角。另外，如图2.30所示弯矩引起的塑性变形将假设集中在aL区段内。图形中阴影部分表示发生塑性变形的区段。因此包含塑性变形和剪切变形的柔度矩阵(flexibilitymatrix)如下LaL"211"21111f11=+l(3-3ot+a)-+豆-+3EI。31ElJ<EI2ElJJGAL(3.a)L12f216EI1111EI。“GALaL-2f11)ia+(3-3a+a31<EI1EI。1GAL(3.b)(3.c)图2.30弯曲刚度的分布假定构件的荷载一位移关系可用柔度矩阵表现如下=fM在此，epsf=fff(5)如图2.31所

12、示，式(5)中各项分别表示弹性弯曲变形角、塑性变形角、因剪切引起的弯曲变形角。剪切变形弹簧轴向变形、扭转变形、图2.31弯矩一变形角关系8斗0口+守在MIDAS/Genfi勺Pushover分析中假定轴力、扭矩、剪力在构件内不变，塑性较发生在构件中央。其荷载一位移关系可参照弯曲变形中的各式。双向弯曲弹簧双向受弯且受轴力作用时，先计算各向的屈服弯矩后使用下面关系式建立双向受弯相关公式。Mnx<MnoxMnyMnoyr=1.0(6)上式适用与钢筋混凝土和钢材等所有构件。塑性较(plastichinge)特性随着荷载的增加，结构构件将产生塑性较，结构的刚度会发生变化，横向位移也将逐渐加大。MI

13、DAS/Ge冲采用的塑性较特性如下。较特性：多折线类型(Multi-LinearType)-采用切向刚度矩阵(tangentstiffnessmatrix)-荷载控制(loadcontrol)和位移控制(displacementcontrol)均可-可考虑P-Delta效果较特性：FEMA类型(FEMAType)时-割线刚度矩阵(secantstiffnessmatrix)- 采用位移控制(displacementcontrol)- 可考虑P-Delta效果和大位移(largedeformation)效果因为结构承受的荷载大小为已知条件，所以一般采用荷载控制方法。荷载控制方法就是将荷载从零开始

14、逐渐加载到极限荷载的方法。位移控制是在基于性能的耐震设计中采用比较多的方法。虽然不知道加载的荷载大小，但是可以通过预先设定满足目标性能的位移进行分析。分析过程中可以获得荷载传递能力(load-carryingcapacity)和失稳(unstable)状态。采用位移控制和割线刚度矩阵(secantstiffnessmatrix)时，在最大荷载附件可以获得稳定的解。多折线较类型多折线较类型可以用于荷载控制和位移控制方法中。- 荷载一位移关系采用双折线(Bilinear)和三折线(Trilinear)两中形式- 屈服后刚度和抗裂刚度用与初始刚度的比值(stiffnessratio)来表现- 能表现

15、构件的刚度降低，但不能表现材料的强度降低daoLLaroraFEM酸类型FEMA较类型是将钢筋混凝土构件和钢构件的循环加载试验(reversedcyclicload)获得的资料理想化的结果，其特性如下。MIDAS/GenH勺FEM段特性只能使用位移控制方法。LateralDeformation图2.33FEMA较类型的塑性较特性点A位置：未加载状态,由材料、构件AB区段：具有构件的初始刚度(initialstiffness)尺寸、配筋率、边界条件、应力和变形水准决定.点B位置：公称屈服强度(nominalyieldstrength)状态BC区段：强度硬化(strainhardening)区段，

16、刚度一般为初始刚度的5-10%,对相邻构件间的内力重分配有较大影响。点C位置：由公称强度(nominalstrength)开始构件抵抗能力开始下降CD区段：构件的初始破坏(initialfailure)状态，钢筋混凝土构件的主筋断裂(fracture)或混凝土压碎(spalling)状态，钢构件的抗剪能力急剧下降区段。DE区段：残余抵抗(residualresitance)状态，公称强度的20%E右点E位置：最大变形能力位置，无法继续承受重力荷载的状态。Pushover分析方法MIDAS/Gen中提供两种Pushover分析方法，即基于荷载增分的荷载控制法和基于目标位移的位移控制法。基于荷载增

17、分的荷载控制法MIDAS/Gen的荷载控制法采用全牛顿拉普森(FullNewton-Raphson)方法。牛顿拉普森方法是采用微分原理求解的方法，其优点是速度快。采用荷载增的Pushover分析方法的图形接介绍如下。荷载图2.34基于荷载增分法的Pushover分析荷载采用具有一定分布模式的横向荷载。荷载分布模式既可以采用地震荷载（Qud）也可以采用任意的荷载分布模式。另外，也可以采用包含节点荷载在内的用户定义的任何荷载工况。第1阶段：计算弹性极限首先使用用户定义的水平荷载计算构件的应力，然后计算各构件的应力与屈服应力的比值九。将各构件的比值中的最小值乘以加载的荷载工况重新定义荷载。'

18、=U-L/M（7）在此：匕各校计算的屈服荷载系数（最大值0.5）U:钱的屈服应力L：初始应力M：荷载工况计算的钱的应力第2阶段：基于等差级数的增分分析由弹性极限到预估的坍塌荷载（Qud*X）之间的荷载增量由下面的等差级数计算。P=（n+1）-il（PxX心i（8）在此：i:等差增分步骤数pi：第i步的荷载增量P:总荷载X：预估的坍塌荷载与总荷载的比值（基本值为0.4）第3阶段：预估坍塌荷载之后的荷载增分使用最终计算的（n+1）步骤的荷载增分。终止分析的条件- 到达最大增分步时- 层间位移角到达极限层间位移角时- 分析中计算的水平内力到达指定的大小时- 刚度矩阵为负（negative）时基于目标

19、位移的位移控制法MIDAS/Gen的位移控制法是由用户定义目标位移，然后逐渐增加荷载直到达到目标位移的方法。目标位移分为整体控制和主节点控制两种，整体控制是所有节点的位移都要满足用户输入最大位移，位移也是整体位移，不设置某一方向的位移控制。主节点控制是用户指定特定节点的特定方向上的最大位移的方法。基于性能的耐震设计大部分是先确定可能发生最大位移的节点和位移方向后给该节点设定目标位移的方法。初始的目标位移一般可假定为结构总高度的1%2%4%这些数值一般相当于最大层间位移值，于结构的破坏情况相关。ATC-40或FEMA-273中将层间位移为1%时定义为直接居住水准(ImmediateOccupan

20、tLevel),2%寸定义为生命安全水准(LifeSafetyLevel),4%时定义为坍塌预防水准(CollapsePreventionLevel)。这些值在构件级别上的意义可能会稍有不同。作用荷载作用荷载应该采用能反映各层惯性力的横向荷载。所以推荐至少使用两种以上的横向荷载分布模式。在MIDAS/Ger提供了三种类型的荷载分布模式，即静力荷载分布模式、振型形状分布模式、与各层质量成比例关系的等加速度分布模式。采用静力荷载分布模式时，用户可以定义任意形状的静力荷载分布。采用振型形状荷载分布模式时必须先做特征值分析基于性能点的耐震性能评价在MIDAS/Gen中使用能力谱（CSM）M理评价结构的

21、保有内力和耐震性能。结构的保有内力可通过Pushover分析获得的能力曲线和能力谱进行评价。对地震作用的需求谱可以适用有效阻尼的弹性设计谱来评价。将这两个谱表现在相同的坐标系上将获得意味着结构非线性最大需求内力的交点，即性能点（performancepoint）。利用性能点位置的变形程度和保有内力来评价结构的耐震性能和性能水准。能力谱和需求谱评价结构的耐震性能和性能水准时会使用能力谱和需求谱的概念。通过Pushover分析将获得荷载一位移关系（V-U）,响应谱也可获得加速度周期（A-T）的相关关系。为了比较两个谱，需要将其转换为力口速度一位移谱（acceleration-displacemen

22、tresponsespectrum,ADRS）图2.35将荷载-位移关系转换为加速度-位移谱DemarrispectruinDTnS图2.36将加速度-周期谱转换为加速度-位移谱如图2.35所示，荷载一位移关系转换为加速度-位移关系的方法如下:(9)在此，rk和Mk为各自方向的(10)k阶振型的振型参与系数和有效质量系数，计算方法如下:N''mj'jkj1振型参与系数N“mj'jkj£(11)N2二mjjk.jaN振型参与质量(12).二mjji)2kj上式(9)和(10)为动力学理论的多自由度(MDOF体系和单自由度(SDOF林系之间的关系。即A和D

23、为单自由度体系响应谱上的响应加速度和响应位移，V和U为多自由度体系的基底剪力和位移。如图2.36,弹性响应谱可以利用单自由度体系的位移和加速度关系式(13)进行转换。(13)性能点(performancepoint)的评价能力谱和需求谱的交点称为性能点。在MIDAS/Gen提供的计算性能点的方法为ATC-40的能力谱(CSM)中提供的Procedure-A和Procedure-B两种方法。两种方法的基本原理相同，通过计算有效阻尼反复计算获得性能点的方法为Procedure-A方法，利用延性比和有效周期原理计算性能点的方法为Procedure-B。(1)计算等效阻尼(equivalentdamp

24、ing)在能力谱法(CSM)中，通过pushover分析获得能力谱后如下图所示使用具有相同面积的双折线(biloinear)曲线来表现。在CSM中使用具有5%fi尼的弹性响应谱和能力谱计算结构的等效阻尼。因为结构的阻尼而耗散的能量等于双折线滞回曲线的面积，可按式(14)计算。UOW网"4u*Eiuuo.M图2.37利用滞回曲线计算等效阻尼-eq+0.050(14)apidpi1Ed63.7(aydpi-dyapi)4二Eso在此，ED=结构阻尼引起的耗散能eso=结构的最大变形能将式(14)使用百分率的形式表现如下。：eq(15)05=63.7(9dpi-dyapi)5apidpi在

25、此，表示阻尼比(),在ATC-40中阻尼比超过25%寸，需要谨慎的判断，且不许超过50%(2)计算有效阻尼(effectivedamping)地震作用作用下的钢筋混凝土结构的滞回曲线中没有表现刚度退化(stiffnessdegradation)和强度退化(strengthdeterioration)、滑移或握裹(sliporpinching)的特性的理想化的滞回模型。所以在ATC-40中为了反映钢筋混凝土的这些滞回特性，使用阻尼调整系数(dampingmodificationfactor)来调整等效阻尼。调整后的等效阻尼称为有效阻尼系数，按下式计算。(16)，:63.7-(aydpi-dyap

26、i)匚eq=05=5apidpi式左侧的阻尼比5恸弹性体系的阻尼，不同结构特性的阻尼调整系数如下表。结构特性等效阻尼国(%)阻尼调整系数()口A(完全滞回特性)<16.25>16.251.01130.51(aydpidyapi)apidpi口B(一般滞回特性)<25>250.670845_0.446(aydpi-dyapi)apidpi口C(弱滞回特性)所有值0.33(3)非线性需求谱使用前面计算的有效阻尼系数决定非线性响应谱。即利用有效阻尼系数计算响应谱的谱折减系数(spectrumreductionfactor,SR)。如图2.27所示加速度和速度的谱折减系数不同。

27、谱折减系数采用了Newmark和Hall(1982)的地基运动扩大系数，加速度的谱折减系数(SRA)和速度速度的谱折减系数(SRV)的计算式如下。根据结构的滞回特性，ATC-40中给出了谱折减系数的下限值。图2.38根据谱折减系数计算的非线性响应谱表2.2结构的滞回特性对应的谱折减系数下限值项目SRAsrv类型A（完全滞回特性）1.000.330.50类型B（一般滞回特性）0.670.440.56类型C（弱滞回特性）0.330.560.67根据上述的计算过程可以获得设计地震作用或线弹性反应谱对3.21-0.68lnSRa63.7k(aydpi-dy3pi)+51apidpi2.120.33fo

28、rTypeA二；0.44forTypeB0.56(forTypeC)SR63.7*(aydpi-dyapi2.31-0.41lnapidpi1.655)0.50forTypeA,0.56forTypeB0.67forTypeC(17)应的非线性需求谱。将获得的非线性地震需求谱和通过Pushover分析获得的结构的能力谱进行比较，可以获得结构的性能点。（4）计算性能点利用Pushover分析得到的结构的能力谱和非线性设计响应谱的比较，可以获得表现结构的非线性最大位移和保有内力的性能点，并且利用其来评价结构的性能水准。确定性能点的方法MIDAS/GerfN艮据能力谱（CSM）!定性能电白方法采用A

29、TC-40中提供的两种方法。其基本原理为使用有效阻尼系数评价需求谱并求其与能力谱的交点作为性能点Procedure-A是ATC-40中提供的基本方法，首先将能力谱中斜率为初始刚度的切线和阻尼比为5%勺弹性设计响应谱的交点作为初始的性能点。然后确定初始性能点位置的等效阻尼，然后求使用有效阻尼系数的非线性设计响应谱，然后重新计算交叉点作为性能点。重复上述过程，直到在使用有效阻尼系数的非线性设计响应谱和能力谱的的交点位置上位移响应和加速度响应的变化量在误差范围内，将此时的交点视为性能点。采用Procedure-A方法确定性能点的方法参见图2.39。SpectralDisplacement,inche

30、s图2.39使用Procedure-A方法计算性能点(ATC-40)Procedure-BATC-40中计算性能点的第二种方法是首先假设位移延性比，然后计算对应延性比的结构的结构的有效周期，将有效周期直线和5%弹性设计响应谱的交点作为初始的性能点。对弈于假定的位移延性比的放射线状的有效周期和非线性设计响应谱的交点将形成一个轨迹线，该轨迹线与结构的能力谱的交点为最终的性能点。利用Procedure-B方法计算性能点的原理如图2.40所示。DemandCurvesfor%=5%,10%,15%,20%,25%and30%CepacHy.pectrorn5%rwpan»LlgginaCfl

31、gpoint*plotttd1tpacHladinBPirformanctpointInttrMCtktnolctpackty«p«ctrumundMn«pfdttedjpAClfladinstapTSpectralDisplacementinches图2.40利用Procedure-B方法计算性能点(ATC-40)RCOIW里aucsegn该方法是首先假定位移延性比，然后逐步计算有效阻尼系数,所以在交叉点计算的响应误差发散的概率较低。前面介绍的Procedure-A方法在寻找性能点的过程当中收敛性不是很好，而Procedure-B方法不仅收敛性能好，而且不必建立

32、针对多个阻尼比的弹性反应谱，而是根据变化的阻尼比和振动周期获得响应谱的轨迹即可获得性能点，所以Procedure-B方法是相对比较简单的方法。在MIDAS/Ge冲提供的两种方法的操作界面如下图所示。PushoverCurvePushoverLo#Case|PUSH1(PUSHPAT)PlolTypeCBdseShearvs.DlsplacernemShearCoetficienlvs.DisplacementCShearCoetficienlvs.DriftAdditionalCurvesatOtherNodesCapacitySpectrumDemandSpectrumDefineDesig

33、nSpeelrum7DemandSpactraatDampingRatios(%)_TO15po|7Constan!PeriodLiinosatPeriods(sec)pfPF-EvaluallanofPerfarmancePoIaif*1Procedure-ArProcedure-BCapacitySpedrumvs,DemamdSnectrumJ6J5。D.ID.行CDHJnKk心石一岳Jl=<DtlmSpeclralDisplacemenKSd)DescrlpilonforPrintedOutput|f=vailuaiion-cfP&rlomnaA«PointCS

34、MProcedure-A(ATC-40)144.4,16i0Sa.SdCL21比3.MTeff.DeffI.3U1E.2Displ,ControlNode：19Dir：DXMaterLoadPattern；StaticLoad-BackgroundColorrBlack/WhileChangeGraphTitleChangeGraphRangeSS谑'箱而16而福，田i丽"JAdditionalPushoverStepforStoryDrift|TextOutput|Redraw|Close2.41利用Procedure-A的方法决定性能点(MIDAS/Gen)Pushove

35、rCurvePushoverLoadCase|PUSH1(PUSHPAT)PloiTypeCBaseShear特.Dis悟cementCShearCoetficienlvs.DisplacementCShearCmtficMlvs.DriftAdditionalCurvesatOtherNodesFPFF<*'CapacliySpectrumDemandSpectrumDsflneDesignSpecirumi|7DemandSpectraatDampingRatios曲)FFFIFCanstantPeriodUnesatPeriDds(sec)Evaluallon口tPerfoirmancePolinirProcedur&-A国Procedure-BDampingParamelerInherent+AddNonalDamping4)StructuralBehaviorType佰三SpeclralDisplacemenKSd)ChangsGraphTitlePerforrriiaricePointDispl,ControlNode：ISDir：DXMaterLoad