震作用下钢筋混凝土结构倒塌反应的数值模拟计算分析.pdf

地震作用下钢筋混凝土结构倒塌反应的数值模拟计 算分析 张颖，倪镇国 辽宁工程技术大学结构工程系，辽宁阜新（123000 ） E-mail ：lncathy@ 摘 要：本文综述了钢筋混凝土结构倒塌反应的国内外研究现状，分析了钢筋混凝土结构倒 塌反应的计算机仿真中的几种数值方法的优缺点，指出了该研究中目前存在的问题和解决问 题的方向。 关键词： 数值模拟；倒塌；有限单元法；离散单元法 中图分类号：TU375.4 1. 引 言 近年来钢筋混凝土结构在地震作用下的内力及变形可能达到极限承载能力而遭受破损， 如结构裂缝的出现与扩展、混凝土保护层的剥落、钢筋的屈服与断裂、连接锚固的破坏等。 构件的局部破坏可引发建筑物的层间破坏或整体倒塌。钢筋混凝土结构[1]是地震作用下破坏 最严重的结构形式。因此，改善钢筋混凝土结构的抗震性能是非常重要的。我国和世界各国 的学者在这方面做了大量的试验与数值分析的研究工作，特别是数值分析法，近年来得到了 很大的重视。 [1]性质本身的很大的离散性及各向异性性、手工操作过程中质量的难以 由于混凝土材料 控制等诸多因素的影响，使混凝土结构试验的数量增大，成本增高，且又有模型与原型的相 似理论的许多问题。应用可视化技术在计算机上模拟[2]钢筋混凝土在地震作用下的破坏全过 程依据可作为一个重要的辅助研究手段，为新结构的设计和旧结构的改造与加固提供理论。 另外，结构在地震中的倒塌方式将直接决定它的影响范围，应用计算机仿真技术分析结构在 地震中的倒塌方式亦可为小区的规划、震后的救援和恢复等提供决策依据， 本文对这一领 域的成果和存在的问题进行了综合分析，提出了解决问题的方向，为进一步完善钢筋混凝土 结构倒塌反应的仿真理论[3]作必要的准备。 2. 钢筋混凝土框架结构倒塌过程数值模拟的研究现状 钢筋混凝土结构的计算机数值模拟技术是研究结构地震反应的重要方法之一。早期的结 构非线性动力响应[4]的数值模拟，大都集中在结构弹塑性阶段的分析。其中比较典型的结构 分析程序有 Park 等人提出的 IDARC 程序，该程序属于平面分析程序[6]，采用梁杆单元、剪 力墙单元、连梁单元、弹簧约束单元等单元[5]，利用三线形恢复力模型，来模拟框架结构强 度、刚度退化等特性。Li 等人建立了 CANNY 程序，该程序有多种恢复力模型、多种数值 计算方法、实用的前后处理功能，是一个比较全面的三维非线性分析程序[7] 。我国的许多学 者在这一方面取得了研究成果，顾祥林、孙飞飞等人运用 Visual Basic 和 Fortran 程序设计语 言编制[8]一套钢筋混凝土平面框架结构地震反应仿真软件，该软件能纪实仿真结构模型上各 单元(梁、柱、节点)裂缝的开展过程。倪强、唐家祥等人在这一方面也做了一定的研究工作。 同济大学的秦东、范立础等人在模拟结构倒塌的过程[9]中提出了一种扩展散体单元模型，这 种模型可以解决结构从局部破坏到整体破坏的全过程、混凝土的开裂、混凝土块体的碰撞与 刚体运动、钢筋的屈服等现象，但这一单元的刚度退化模型(恢复力模型)与实际恢复力模型 - 1 - 有一定的差距，还需改进处理。焦双健、冯启民等人开发了 DS5·0(Damage Simulation 5·0) 系统程序，建立刚度沿杆长变化的破坏单元模型[10]，采用考虑变形和滞回耗能的双重破坏 指数作为单元的破坏准则，对于弯矩释放的破坏单元。采用以变形作为判断单元破坏的准则， 建立了判断整体结构是否为机构的机动分析方法，提出贯穿牵连方法，来区分框架结构的局 部倒塌和整体倒塌。目前随着计算机技术的发展，特别是大型有限元软件如 ABAQUS 、 ANSYS 、SAP、PKPM 等软件的发展，使得可以利用这些软件丰富的单元类型、众多的材 料本构模型以及强大的后处理功能等特点，在三维实体中分别对结构的梁柱、节点乃至整体