剪力墙结构的静力弹塑性分析
【摘 要】简述了静力弹塑性分析的原理，通过工程实例进行相关的对比和讨论。
【关键词】静力弹塑性分析；剪力墙结构；结构抗震性能评价；EPDA
1、前言
一般建筑结构的抗震设计都需要考虑结构的弹塑性行为。由于时程分析法计算工作量大， 结果处理繁杂， 相比之下， pushover法更方便于进行抗震设计。尤其是上世纪 90年代以后， 随着基于性能的抗震设计思想的提出和发展，pushover 方法得到了深入的研究和日益广泛的应用。我国在新的建筑结构抗震设计规范中也引入了pushover方法。
利用pushover方法进行结构分析的优点在于：既能对结构在多遇地震下的弹性设计进行校核，也能够确定结构在罕遇地震下潜在的破坏机制，找到最先破坏的薄弱环节，从而使设计者仅对局部薄弱环节进行修复和加强，不改变整体结构的性能，就能使整体结构达到预定的使用功能。
对多遇地震的计算，可以与弹性分析的结果进行验证，看总侧移和层间位移角、各杆件是否满足弹性极限要求，各杆件是否处于弹性状态；对罕遇地震的计算，可以检验总侧移和层间位移角、各个杆件是否超过弹塑性极限状态，是否满足大震不倒的要求。
2、原理与实施步骤
2.1 原理
静力非线性分析方法是基于性能评估现有结构和设计新结构的一种方法。它是将静力弹塑性分析和反应谱相结合进行图解的快速计算方法。其原理是使结构分析模型受到一个沿结构高度为单调逐渐增加的侧向力或侧向位移 ，直至控制点达到目标位移或建筑物倾覆为止。
基于结构行为设计使用Pushover分析，包括形成结构近似需求曲线和能力曲线，并确定曲线交点。需求曲线基于反应谱曲线，能力曲线基于静力非线性Pushover分析。在Pushover分析中，结构受到逐渐增加的荷载作用，从而得到需求曲线和能力曲线的交点，即性能点。由于性能点定义了结构的底部剪力和位移，因此通过结构在性能点的行为和现行规范进行比较，从而确定结构是否满足要求。
2.2实施步骤
（1）准备工作：建立结构模型，包括几何尺寸、物理参数以及节点和构件的编号，并输入构件的实配钢筋以便求出各个构件的塑性承载力。
（2）求出结构在竖向荷载作用下的内力。这时还要求出结构的基本自振周期。
（3）施加一定量的水平荷载。水平力大小的确定原则是：水平力产生的内力与第（2）步竖向荷载产生的内力叠加后，恰好能使一个或一批构件进入屈服。
（4）对在上一步进入屈服的构件的端部，设定塑性铰点变更结构的刚度，这样，相当于形成了一个新的结构。求出这个“新”结构的自振周期，在其上再施加一定量的水平荷载，又使一个或一批构件恰好进入屈服。
（5）不断地重复第（4）步，直到结构的侧向位移达到预定的破坏极限。记录每一次有新的塑性铰出现后结构的周期，累计每一次施加的荷载。
（6）成果整理：将每一个不同的结构自振周期及其对应的地震影响系数绘成曲线，也把相应场地的各条反应谱曲线绘在一起，如图1所示。这样如果结构反应曲线能够穿过某条反应谱，就说明结构能够抵抗那条反应谱所对应的地震烈度。还可以在图中绘出相应的变形，更便于评价结构的抗震能力。
3、工程实例
某工程为28层剪力墙结构，用SATWE进行过振型分解法分析，现用EPDA进行push-over法的分析，判断结构的抗震性能，并与SATWE程序的计算结果进行对比。
4、结论及问题讨论
4.1楼板对梁刚度的影响
从以上工程实例与SATWE的对比结果可以看出，在多遇地震下，当SATWE的“中梁刚度增大系数”Bk=1时，结构的自振周期SATWE与push-over的结果基本一致，而当Bk=2时，SATWE计算的自振周期比push-over要短，说明了push-over计算中并未考虑楼板对梁刚度的影响，程序中也没有提供相关的参数设置，建议程序对此作进一步改进。
4.2周期折减系数Tc的影响
从以上的对比中，还可以看出“周期折减系数”Tc的影响，在多遇地震下，当Tc<1时，SATWE计算的基底剪力和位移均大于push-over的计算，只有当Tc=1时，两者的基底剪力和位移才较为接近。Tc在SATWE中的作用是考虑填充墙的抗侧刚度对结构自振周期的影响，这个对比显示push-over忽略了这项因素。简单的修正方法是对求出来的周期乘上一个周期折减系数，然后在用于结果的整理。
周期折减系数的取值要注意的是当结构的侧向变形达到一定范围后，由于填充墙的逐渐破坏，周期折减系数要相应调整，以反映填充墙抗侧刚度的降低甚至消失。
4.3混凝土本构关系模型的假定
本文计算实例的push-over法采用的混凝土应力应变关系模型为Saenz曲线，如图4所示，在应变0～εc范围内（即混凝土的弹性范围），应力应变的关系是非线性的，混凝土的弹性模量E（即图中曲线的切线斜率）随应变、应力的增大而减小。结构的刚度随弹性模量E的减小而减小。自振周期随结构的刚度的减小而增大。因此，自振周期随影响系数的增大而增大。
一般的弹塑性分析程序采用的混凝土应力应变关系模型如图5所示，其中应变0～εc