静力弹塑性分析原理

第4章静力弹塑性分析原理4.1概要4.1.1非线性分析的目的非线性抗震分析方法可分为非线性静力分析方法和非线性动力分析方法。非线性动力分析方法可以认为是比较准确的方法，但是因为分析时间较长并对技术人员理论水准有较高的要求，所以在实际工程上的普及应用受到了限制。相反静力分析方法虽然在反映结构动力特性方面有所不足，但是因为计算效率较高和操作简单、理论概念清晰等原因被广大设计人员所普遍使用。静力弹塑性分析又被称为Pushover分析，是基于性能的抗震设计(Performance-BasedSeismicDesign,PBSD)中最具代表性的分析方法。所谓基于性能的抗震设计是以某种目标性能(targetperformance)为设计控制目标，而不是单纯的满足规范要求的极限承载能力的设计方法。其步骤是先按照规范要求进行抗震分析和构件设计，然后通过Pushover分析获得结构的极限承载能力，最后通过非线性位移结果评价结构是否满足目标性能要求。目前规范中推荐的基底剪力法和反应谱分析方法均为弹性分析方法，其评价标准是地震作用下的抗力不小于地震作用下产生的内力，这些方法也被称为基于荷载的设计方法。而基于性能的设计方法则是使用与结构损伤直接相关的位移来评价结构的变形能力(耗能能力)，所以又被称为基于位移的设计(displacement-baseddesign)方法。通过Pushover分析可得如图4.1.1所示的荷载-位移关系曲线(能力谱)，根据结构耗能情况可得到非线性需求谱。能力谱与需求谱的交点就是结构对于地震作用的性能点(performancepoint)。性能点意味着结构对于地震作用所拥有的最大的非线性承载力和最大位移，该点在控制目标性能范围内则表示该结构满足了性能要求。通过非线性分析可以了解结构具有的的极限承载能力和安全度。5%SmoothElasticDesignSpectraPerformancePointDemandSpectrumCapacitySpectrumSpectralDisplacementSpectralAccelerationSaSdDdesign图4.1.1基于位移设计法的结构抗震性能评价4.1.2静力弹塑性分析的抗震设计原理基于性能的抗震分析方法有下列四种。线性静力分析法(LinearStaticProcedure,LSP)线性动力分析法(LinearDynamicProcedure,LDP)非线性静力分析法(NonlinearStaticProcedure,NSP)非线性动力分析法(NonlinearDynamicProcedure,NDP)其中Pushover分析方法属于非线性静力分析法，又被称为塑性铰分析法。该分析方法主要PushoverAnalysisCapacitySpectrumdSaSSDOFSystemroofroofCapacityCurveMDOFSystemtransformFbaseVbaseV被应用于受高阶振型和动力特性影响较小的结构。Pushover分析就是按照指定的加载模式逐渐加载至控制目标并获得结构的荷载-位移曲线(capacitycurve)，然后将其转换为单自由度体系的加速度响应和位移响应的能力谱，同时将加速度-周期格式的加速度反应谱转换为加速度-位移格式(ADRS，Acceleration-DisplacementResponseSpectrum)的需求谱(demandspectrum)，将需求谱和能力谱反映在同一个坐标系中，两条谱曲线的交点(性能点)就是满足该水准地震作用的极限承载能力和变形能力点。因此可通过定义不同的需求谱(小震、中震、大震)，通过验算不同性能水准下的承载力和变形，实现“小震不坏、中震可修、大震不倒”的三水准(performancelevel)抗震设计原则。结构大师中使用了ATC-40(1996)和FEMA-273(1997)等报告中的能力谱法(CapacitySpectrumMethod,CSM)和推荐的参数对构件的抗震性能进行评价。能力谱法的原理如图4.1.2所示。(a)结构的能力曲线(capacitycurve)和能力谱(capacityspectrum)(b)需求谱(demandspectrum)(c)性能点(performancepoint)图4.1.2能力谱法(CapacitySpectrumMethod,CSM)Pushover分析的目的是要了解结构具有的承载能力和变形能力，钢筋砼结构在进行Pushover分析前必须先进行线弹性分析和构件设计以获得结构的配筋结果。然后才能进一步进行非线性分析。Pushover分析的优点如下：可获得结构屈服后的响应和极限承载能力可获得结构耗能能力和位移需求可获得结构构件的出铰顺序在维修加固工程中事先了解需要加固的构件4.1.3静力弹塑性分析方法如下图4.1.3所示，结构在横向荷载作用的初期处于弹性状态，当内力超过构件的开裂或屈服内力时部分构件将发生开裂或屈服，构件和结构的刚度和阻尼都将发生变化，荷载和位移的相关关系显示非线性特性。由弹性进入屈服阶段的点A被称为弹性极限，部分构件屈服后随着荷载的增加结构的位移会显著增加，到达B点后较小的外力增量也会发生较大的位移，最后在C点后即使不再增加外力位移也会增加，C点被称为极限承载能力点。DemandSpectrumn,2Ta22ndS4TSn,1TnTdSaSaSResponseSpectrumtransform5%ElasticSpectrumPerformancePointDemandSpectrumCapacitySpectrumaSdSmaxDmaxAInternalForceDisplacementAElasticRangeCInelasticRangeBPlasticHinge图4.1.3内力和位移的关系点C是通过荷载增量进行分析的荷载控制法所能得到稳定解的极限点，要想获得C点之后的曲线只能通过位移增量进行分析，即采用位移控制法。结构大师中既提供荷载控制法又提供位移控制法。4.1.4迭代分析方法Pushover分析中由于发生裂缝和屈服造成结构的刚度变化，在分析过程中会产生不平衡力也叫残余力(ResidualForce)，为了消除不平衡力需要进行迭代计算使不平衡力达到可以忽略的程度(满足收敛条件)。结构大师的迭代计算方法使用了完全牛顿-拉普森法(FullNewton-RaphsonMethod)，该方法具有收敛速度快的特点。使用完全牛顿-拉普森法的非线性分析过程如下图所示，分析过程如下。(1)nF(0)nK1nU(3)nU(1)nU(2)nU(4)nUU;ReferenceLoad;Displacement;ResidualForce0Pn10Pn(2)nK(1)nK0P(1)R(2)R()iRn(1)Un(2)Un(3)U0PnAB;InternalForce()inFn(2)UC图4.1.4完全牛顿-拉普森法(1).在当前步骤(n)增加荷载向量0nP可得图4.1.4所示的A点，此时的平衡方程式如下。10nnnnKUFP(4.1.4a)其中，nK:当前步骤(n)的结构切线刚度矩阵nU:当前步骤(n)的位移增量1nF:前次步骤(n-1)的内力向量n:当前步骤(n)的加载系数0P:荷载向量0nP:当前步骤(n)的荷载向量可将式(4.1.1)用增量形式表达如下，010nnnnnnnKUPFKUP(4.1.4b)其中，0nP:当前步骤(n)的荷载增量向量解式(4.1.2)得位移增量nU。(2).利用位移增量nU计算各单元的切线刚度和内力，将各单元的内力组合构成切线刚度矩阵()inK。将各单元的内力与节点力组合构成内力向量()inF。此时结构的内力和位移的关系满足图4.1.4点B上的平衡条件。(3).荷载增加0nP时如果单元发生屈服则单元产生残余力()inR，可通过下面的迭代计算消除残余力。()()()0()()()iiinnnniiinnnKUPFKUR(4.1.4c)其中，()inK:当前步骤(n)内的第i次迭代计算时的切线刚度矩阵()inU:当前步骤(n)内的第i次迭代计算时的位移向量()inF:当前步骤(n)内的第i次迭代计算时的内力向量()inR:当前步骤(n)内的第i次迭代计算时的残余力解方程(4.1.3)得位移向量()inU。计算各单元的内力和切线刚度后可得残余力()inR，重复(1)~(3)步骤直到满足收敛条件。(4).满足收敛条件时(在点C)将进行下一个增量步骤的分析。残余力和收敛计算塑形铰的出现造成了单元刚度的变化，单元刚度的变化又引起了单元内力的变化，从而使外力和单元内力之间产生了不平衡力(残余力)。程序中消除残余力的方法如下：进行收敛迭代计算时(在时程荷载工况对话框中勾选了迭代计算选项)使用完全牛顿-拉普森法进行迭代收敛计算直至满足收敛条件。但是仍有下面的残余力累计到下一个增分步骤的外力中。当最大迭代次数输入1时等同于不进行迭代计算。a.到最大迭代次数时仍未满足收敛条件时的残余力b.满足了收敛条件但仍残留的不平衡力不进行收敛迭代计算时(在时程荷载工况对话框中未勾选迭代计算选项)各增量步骤的残余力将累计到下一个增量步骤中的外力中。因此即便是某个增量步骤中没有收敛只要下一个步骤中收敛时，可以认为最终分析结果收敛。收敛判断条件因为不可能完全消除残余力，所以为了既满足计算结果的精确度又保证计算效率，需要设置适当的收敛判断条件。迭代计算的收敛判断采用范数标准，有位移范数、荷载范数、能量范数，可选择其一也可多选作为收敛判断标准。位移范数()()()()iTinnDiTinnUUUU(4.1.4d)荷载范数()()()()iTinnFiTinnFFFF(4.1.4e)能量范数()()()()iTinnEiTinnFUFU(4.1.4f)其中，D:位移范数F:荷载范数E:能量范数()inU:当前步骤(n)内的第i次迭代计算累计的位移增量向量()inU:当前步骤(n)内的第i次迭代计算的位移向量()inF:当前步骤(n)内的第i次迭代计算的累计内力增量向量()inF:当前步骤(n)内的第i次迭代计算的内力向量4.1.5初始荷载因为地震作用前结构的竖向荷载是始终存在的，所以有Pushover分析有必要考虑竖向荷载作用下的初始内力状态，这样计算的杆件的内力才是接近真实的。特别是考虑轴力和弯矩相关的柱构件在计算屈服面时需要考虑竖向荷载引起的轴力。结构大师的Pushover分析对初始荷载也进行非线性分析，以获得更接近于实际情况的初始内力。4.1.6加载模式Pushover的横向荷载应该能相对准确地反映实际地震作用，即实际地震力在各楼层的惯性力分布状态，这样才能保证分析结果更接近于实际状态。一般来说，在Pushover分析中推荐使用两种以上的横荷载分布模式进行分析，通过比较取不利的结果进行判断。在结构大师中提供了四种横向荷载加载模式，特别是其中的层剪力模式使用了反应谱分析得到的层剪力分布模式，更接近于实际的地震力分布。静力荷载模式：按用户定义的静力荷载分布加载。振型模式：可按振型形状分布模式加载，也可以将几个振型线性组合。按各楼层的质量分布比例加载按反应谱分析的层剪力分布模式加载4.1.7P-Delta效应结构大师的Pushover分析支持梁、柱、支撑的P-Delta效应，非线性墙单元目前不支持P-Delta效应。P-Delta分析属于几何非线性分析，单元刚度矩阵中要考虑几何刚度，考虑P-delta效应的平衡方程如下。10GnnnnKKUFP(4.1.7a)其中，K:结构的弹性刚度矩阵GK:结构的几何刚度矩阵在Pushover分析过程中刚度矩阵行列式应大于零，当刚度矩阵的行列式为零或负值时将忽略几何刚度的影响。0GKK(