超高层静力弹塑性分析

超高层建筑结构静力弹塑性分析

于定富

西安建筑科技大学建筑深圳设计分院，深圳(518000)

E-mail：yudfu@

摘要：上海某北塔楼是带加强层巨型柱核心筒超高层结构, 塔楼从室外地面至主屋面装饰体之高度约为260m，地上共58层。单体建筑面积 ( 地上 )141,890 平方米。主屋面上设约 15 至 20m 高钢桁架构筑物。塔楼范围设 3 层地下室 ( 其中局部有二夹层 ) 。地下1，2，3层之层高分别为6.2m，8.3m，3.55m（不含疏水层厚度）,存在加强层及伸臂桁架的钢-混凝土组合结构.静力弹塑性分析采用中国建筑科学研究院的程序PKPM系列PUSH进行。该程序是一个三维有限元空间弹塑性静力分析程序，利用PKPM的PUSH建立了上海某北塔楼的三维非线性有限元模型,并进行了推覆分析,得到了结构的层间位移和层间位移角,找出了结构

的薄弱部位,分析了构件的屈服和破坏规律,对同类工程的分析有较高的参考价值。

关键词:结构工程；弹塑性分析；push- over；反应谱；结构抗震性能评价

1. 引言

现行的《建筑抗震设计规范》（G B50011-2001）中，3.6.2 条为：“......罕遇地震作用下的弹塑性变形分析。此时，可根据结构特点采用静力弹塑性分析或弹塑性时程分析方法[2]。”这里的静力弹塑性分析，即静力非线性分析，除了指一般的与反应谱结合不密切的非线性静力分析外，也包括了push-over方法。《抗规》条文说明5.5.3 明确提出“......较为精确的结构弹塑性分析方法，可以是三维的静力弹塑性（如push-over方法）......”。因为弹塑性时程分析对计算机软硬件和分析人员要求较高，工作量也较大，在一段时期内不容易成为一种被广泛采用的方法。因此逐步推广push-over这种较一般静力分析有许多改进而且相对简便易行的方法，在目前是一种可行的方向[3]。

2. 原理与实施步骤

2.1 原理

Push-over方法是近年来在国外得到广泛应用的一种结构抗震能力评价的新方法，其应用范围主要集中于对现有结构或设计方案进行抗震能力的估计。这种方法从本质上说是一种静力非线性计算方法，与以往的抗震静力计算方法不同之处主要在于它将设计反应谱引入了计算过程和计算成果的工程解释。这种方法的优点在于：水平力的大小是根据结构在不同工作阶段的周期由设计反应谱求得，而分布则根据结构的振型变化求得[1]。

2.2 实施步骤

（1）准备工作：建立结构模型，包括几何尺寸、物理参数以及节点和构件的编号，并输入构件的实配钢筋以便求出各个构件的塑性承载力。（2）求出结构在竖向荷载作用下的内力。这时还要求出结构的基本自振周期。（3）施加一定量的水平荷载。水平力大小的确定原则是：水平力产生的内力与第(2)步竖向荷载产生的内力叠加后，恰好能使一个或一批构件进入屈服。（4）对在上一步进入屈服的构件的端部，设定塑性铰点变更结构的刚度，这样，相当于形成了一个新的结构。求出这个“新”结构的自振周期，在其上再施加一定量的水平荷载，又使一个或一批构件恰好进入屈服。（5）不断地重复第(4)步，直到结构的侧向位移达到预定的破坏极限。记录每一次有新的塑性铰出现后结构的周期，累计每一次施加的荷载。（6）成果整理：将每一个不同的结构自振周期及其对应的地震影响系数绘成曲线，也

把相应场地的各条反应谱曲线绘在一起.这样如果结构反应曲线能够穿过某条反应谱，就说明结构能够抵抗那条反应谱所对应的地震烈度。还可以在图中绘出相应的变形，更便于评价结构的抗震能力。

3. 工程实例

北塔楼位于场地之西北角。塔楼从室外地面至主屋面装饰体之高度约为260m，地上共58层。单体建筑面积 ( 地上 )141,890 平方米。主屋面上设约 15 至 20m 高钢桁架构筑物。塔楼范围设 3 层地下室 ( 其中局部有二夹层 ) 。地下1，2，3层之层高分别为6.2m，8.3m，3.55m（不含疏水层厚度）,存在加强层及伸臂桁架的钢-混凝土组合结构.静力弹塑性分析采用中国建筑科学研究院的程序PKPM系列PUSH进行。该程序是一个三维有限元空间弹塑性静力分析程序，程序单元库包括梁（柱）元和剪力墙元两种非线性单元。弹塑性梁（柱）单元采用标准的有限元方法（所谓的微观方法）构造，单元切线刚度直接基于混凝土材料微元和钢筋材料微元的本构关系，这种模型通常也被称为纤维束模型。弹塑性墙元面内刚度的力学模型采用平面应力膜，并可以考虑开洞，与梁柱元一样，它的单元切线刚度也直接基于混凝土材料微元和钢筋材料微元的本构关系。墙元的面外刚度用简化的弹塑性板元进行考虑。采用PUSH程序进行静力弹塑性分析的计算模型与SATWE分析模型一致。为了提高计算效率，计算中去掉了作为上部结构嵌固端的地下室部分。考虑到北塔楼结构X向、Y 向在地震作用下抗侧性能的相似性，且X向相对较弱，故现设计仅对北塔楼的X向进行推覆分析[4]。

3.2 结构材料模型及假设

程序计算中采用的混凝土应力应变关系模型为Saenz曲线，如图1所示。在应变0~εc 范围内（即混凝土的弹性范围），应力应变的关系是非线性的，混凝土的弹性模量随应变、应力的增大而减小，程序忽略了混凝土的抗拉能力；钢筋的应力应变曲线采用双折线模拟，如图1所示，其中塑流段有一微小斜率，其切线模量Ep=0.01E 。

在本次计算中，设定当杆件的刚度退化为初始刚度的20%时,显示为出现塑性铰。

3.3 结构分析方法

a.施加竖向静力荷载，结构的计算内力作为侧推分析的初始应力。

b.给定侧推荷载形式，采用弧长控制增量法进行非线性静力推覆分析。

形分布与第一振型地震力分布形式接近,如图2所示。使用程序提供的需求谱与能力谱关系确定对应的多遇地震、设防烈度地震和罕遇地震下的结构反应特性和弹塑性发展情况。

推覆力分布方式程序中提供了两种侧向力分布形式，一种是倒三角形，一种是矩形.静力推覆分析作为结构一种非线性地震反应分析的简化方法,可以对已建、待建或初步设计已完成的结构进行抗震性能的评估,已在实际工程中得到不同程度的应用.水平力加载模式分别为均匀分布、倒三角形分布和抛物线分布.目前对不同加载模式的优劣,研究者尚未达成共识,大多数工程应采用倒三角形分布的加载模式,一般认为分布形状在加载过程中恒定不变[6].一般认为倒三角形加载模式只适合于刚度变化规则的结构.这也是目前国内外抗震设计规范中大多采用的侧向力分布形式,即假定水平侧向力沿高度为倒三角形.由于本塔楼刚度沿竖向分布不均匀,本设计时用了倒三角形和均匀分布两种加载模式进行计算.由于倒三角形分布与矩形分布的推覆分析结果非常接近，所以本论文选用倒三角形分布形式进行静力推覆分析。推覆力分布方式与反应谱地震分析楼层地震力分布比较图见图2.

图2a推覆力倒三角形分布方式图2b 第一震型楼层地震力分布

3.4 结果分析

3.4.1 结构推覆曲线

按倒三角形加载方式分析得到的顶点位移-荷载曲线见图3。

推覆曲线（顶点位移-荷载曲线）显示了推覆全过程中塔楼顶点位移的变化过程。从图中可以看出，当水平推覆力由零增加到相当于设防烈度下基底剪力时，位移-荷载曲线基本保持直线，仅靠近设防烈度点时开始显示斜率微小的下降，此时第一批塑性铰开始出现，少量连梁进入塑性状态。随着推覆力的进一步加大，曲线的斜率明显下降，表现出明显的弹塑性性能。但即使到了罕遇地震对应点，该曲线仍然保持向上的趋势，说明结构仍具有稳定的抗侧推能力。此时塔楼顶点的位移角为1/215。

本次推覆计算中塔楼没有达到结构的极限荷载，从推覆曲线可以看出塔楼具有良好的抗推覆和抗倒塌能力。