美国世贸中心大楼倒塌仿真分析.pdfVIP

土木工程学报,China Civil Engineering Journal,34(6), 2001,8~10 世界贸易中心飞机撞击后倒塌过程的仿真分析 陆新征 江见鲸 清华大学土木工程系 北京 100084 摘要：本文利用动力有限元程序LS-DYNA 对纽约世界贸易中心受飞机撞击后的倒塌进行了 力学分析和仿真，并根据计算结果进行了参数讨论。仿真计算的结果与真实倒塌过程非常接 近，说明通过适当的选取计算参数和计算模型，可以对这种特殊的复杂破坏过程进行模拟分 析和仿真。计算结果说明，世界贸易中心倒塌的直接原因是火灾导致的钢材软化和楼板塌落 冲击荷载引起的连锁反应。如果能够提高结构的抗火能力或者提高结构的延性，将有可能避 免结构倒塌，避免惨剧再次发生。 关键词：动力有限元，世界贸易中心，倒塌 一、引言 2001 年9 月11 日，美国纽约的世界贸易中心双塔陆续受到飞机撞击，导致世界贸易中 心南北两塔完全倒塌，造成5 千余人死亡和重大财产损失。该惨剧受到世界结构工程界的密 切关注。很多人对双塔的倒塌原因、倒塌模式提出了各种各样的见解。一般公认的原因是大 厦由于火灾和堆载而导致倒塌，同时也有一些专家提出了包括二次破坏、燃爆等可能性。由 于大厦倒塌的过程十分复杂，而且很难试验重现，因此目前的分析都是定性的，还没有一个 定量的结论。然而，现代计算机技术的发展为我们提供了一个在计算机上重现整个倒塌过程 的可能，并可以根据计算结果来讨论结构倒塌的原因及如何防止结构的类似倒塌。 二、计算模型 本文采用美国Livemore 软件公司开发的LS-DYNA 动力有限元分析软件进行动力仿真 分析。 由于纽约的世界贸易中心为典型的筒中筒结构。其外筒为密柱深梁结构，钢柱的宽度达 到476.25mm ，而间距仅有558.8mm 。而深梁的高度有 1219.2mm。为了简化计算模型并降 低自由度数量，在计算中将外侧的密柱深梁筒体和内侧的钢桁架筒体都用板单元近似，板单 元厚度设定的原则是按等截面积方法代替。这种近似筒体的整体弯曲和轴向变形与实际情况 吻合，而构件的局部弯曲则与真实情况不尽符合，考虑到筒中筒结构的主要变形是整体弯曲 变形和轴向压缩变形，所以用这种方法近似是合理的。 计算模型中，内、外筒材料设定为钢材，本构关系选用LS-DYNA 中的材料Material 3 ， Plastic Kinematic 模型，材料密度为7800Kg/m3 ，弹性模量为200GPa ，泊松比为0.27，强度 参考美国A440 钢材设定其屈服强度为310MPa，屈服后强化模量为2GPa ，由于涉及到断裂 分析，所以断裂失效塑性变形分别设定为0.5%、1%、5%作为讨论参数。 计算模型的楼板设定为是钢筋混凝土材料，本构关系依然使用 Material 3 ，Plastic Kinematic 模型，材料密度为2500Kg/m3 ，弹性模量为30GPa，泊松比为0.2，30MPa，屈服 后强化模量为0，所有断裂失效塑性变形设定为混凝土最大受压变形，即0.38% 。 由于钢材耐火性差，在撞击后飞机燃油燃烧造成高温对钢材强度有显著影响，对此设定 了火灾后钢材材料，它其他参数均和普通钢材相同，只是其弹性模量和强度都被降低到基本 值的 1/20，即约相当于700℃时钢材的性能。 由于倒塌分析本身非常复杂，因此，计算中选用了LS-DYNA 提供的Single Face Erosion 接触计算模型。该接触计算模型可以自动搜索接触面，判断接触，并可以处理侵蚀、断裂等 复杂边界变化情况。材料的摩擦系数统一设定为0.25 。 1 * 国家自然科学基金资助项目（项目号） 土木工程学报,China Civil Engineering Journal,34(6), 2001,8~10 由于LS-DYNA 是动力计算软件，在施加重力荷载后，结构会上下震动一段时间，这和 实际情况是不相符合的。因此，整个计算需要分为两个阶段：第一个阶段，大厦尚未受到撞 击，施加重力加速度，同时给模型一个比较大的阻尼（这里计