基于抗震规范的Abaqus隐式算法的高层建筑结构时程分析

详细信息

作者优秀: 优秀教师/意见领袖/博士学历/特邀专家/独家讲师
平台推荐: 主编推荐

3年前浏览3481

一、写在前面

2021年5月18日下午，位于深圳市华强北商圈的赛格大厦出现强烈晃动现象，一位当时在赛格大厦电子企业工作者坦言大厦出现明显晃动后，他们没法在高楼里安心工作。正当人们还不明具体原因时，5月19日中午和20日中午大楼再次出现了晃动。虽然晃动的感觉没有18日强烈，但依然引发了一定的恐慌情绪，有些公司将自己的重点文件和设备打包带离了赛格大厦，另觅地点存放。

为了安全起见，5月21日，深圳赛格集团发布通知称，自21日零时起，暂停所有业主、商户、租户进出赛格大厦写字楼和电子市场，待相关检测工作完成后再有序开放。

大楼晃动的原因是什么呢？有人说是地震，也有人说是风振动引起的，还有人说是大楼内的大量机械振动引起的。不过根据深圳气象局和应急管理局的监测数据可以发现，事发当天并未出现明显极端天气，大楼内的机械振动现象也不太可行，因为晚上大楼会断电。

在前期对大楼的房屋沉降、倾斜率、加速度等情况持续进行实时监测后，检测机构判断大楼的主体结构处于安全范围内，而且自6月1日起将会对大楼进行激振试验，以便获得大楼的振动频率，振型和阻尼等重要信息，力图通过试验方法找出大楼晃动的真实原因。

设计院的工程师在进行建筑结构分析设计时，不可能对每栋大楼进行试验，一般是先进行模态分析获取频率和振型，按照归一法处理各阶模态并按照由小到大的方法列出各个频率或周期和振型（当外界荷载与某一阶频率相近时会发生共振现象，振动形态即为与此时频率对应的振型），在此基础上设计人员开展动力学分析获取在动力荷载作用下的内力，位移时程曲线，为结构设计提供数据支持。

动力学分析一般有显式分析和隐式分析两类方法，隐式分析相比显式分析的最大好处是可以做到无条件稳定，隐式算法对线性系统可以不用考虑时间步长的选择，从而避免了显式算法为了满足稳定性要求，需要选择过小的时间步长所带来的计算成本。本文将给大家介绍隐式分析在高层建筑中的应用。

二、原理介绍

1、ABAQUS/Standard采用的隐式算法H.H.T算法(Hilber, Hughes, and Taylor (1977)),是Newmark-β的延伸，Newmark-β法计算原理如下：

对于一个承受动力的系统而言，假设在时间步内，加速度a是介于 i和 i 1之间的一个常量，则有

或是

通过在时间步

间对加速度a的积分可以得到

时刻的速度和位移：

联立以上各式可以解得两个基本递推公式

进而可以解得

时刻的速度与加速度公式：

在计算过程中，由于外荷载

为

时刻的位移，速度加速度及t 1时刻的外荷载函数，故是已知的，可以由KU=P先求出

，进而求出

时刻的速度与加速度，循环以上步骤（采用牛顿法）得到所有后续离散点上的位移速度及加速度。在最后一步（伪动力分析）利用向后欧拉法对计算结果进行修正，采用下列公式对结果进行调整

由此对应不同的α=0.0（无阻尼系统），有多种逐步积分法。取值γ=0.5(a=0)可以达到二阶精度，为了达到无条件稳定性分析效果，一般采用γ=2β的参数组合来实现。

对应a的其它情况则有：

a=-0.05可以消除高阶频率,并对低阶频率无影响；

a=-0.333，阻尼最大。

2、同时为了评价对所选择的时间步的准确性，引入了半增量误差这一参数(Hibbitt and Karlsson (1979))。其计算公式如下:

数值分析结果表明当这个参数是外力幅值的0.1倍时，结果精度高。

三、工程案例实操

《抗规》5.1.2规定：3 特别不规则的建筑、甲类建筑和表5.1.2-1所列高度范围的高层建筑，应采用时程分析法进行多于地震下的补充计算；

所谓“补充”，主要指对计算结果的底部剪力、楼层剪力、和层间位移进行比较，当时程分析法大于振型分解反应谱法时，相关部位的构件和配筋做相应的调整。

同时在地震波选择时需要注意地震动三要素（频谱特性、有效峰值和持续时间）需满足相关规定：工程所在场地决定频谱特性；加速度有效峰值按照表5.1.2-2采用；有效持续时间可取为结构基本周期的5~10倍。

现有长沙市天心区一小区主体结构为30层剪力墙住宅，标准层层高2.9m，总高为87.0m，其平面图如下所示：

在YJK或第三方软件建模完成后导入至ABAQUS如下图所示（为保证两个软件模型的一致性，可以对比两个模型的质量，和基本周期及振动形态，ABAQUS的结果数据在dat文件查找）：

导入模型后，在ABAQUS中的主要参数设置如下：

材料采用混凝土，阻尼参数：α=0.38，β=0.006；

分析步：先进行静力分析，再进行隐式分析；静力分析时间可取为1e-05，隐式分析时长为输入地震波的时间，增量步长同地震波数据，一般取为0.02，其它数据默认（γ，β数值可根据α进行调整；半增量误差参数设置参考上节内容）。

荷载输入：分别输入不同方向的地震波，放大系数为规范表5.1.2-2所采用的加速度峰值除以地震波数据的最大值。

本文选择帝谷波ELCENTRO作为地震动输入，其图形如下图所示：

分别输入不同方向下的地震波数据后得到各个方向的计算指标：

X向基底剪力如下图所示：

X向顶点时程曲线如下图所示：

X向中间层某点的位移曲线如下图所示：

Y向的计算结果与X向一样的处理方法：输入地震波数据，提取基底剪力与位移时程曲线，后续工作主要为与振型分解法的结果进行比较，当大于反应谱法结果时需要进行内力及配筋的调整。

四、小结

本文从ABAQUS的隐式分析算法原理，软件设置中的关键参数理解，国家规范及实际案例这四大方面介绍了H.H.T算法在高层建筑结构分析中的应用，对实际复杂建筑结构分析中具有很好借鉴作用，同时可以作为超限结构分析的补充方法。

对本文案例模型感兴趣的朋友，请在文章末尾点赞并好看，且分享本文到朋友圈后截图发到本公众号（仿真秀App），我会在24小时内把本文相关模型和学习资料发给你。

在重要的超高层建筑结构或大跨空间结构设计时，往往需采用节点试验与有限元分析相结合来确定节点的安全。由于节点实验或非常规构件的实验比较费时，业主不愿因为局部的不确定来影响工期，故采取可靠成熟的商业软件来进行数值模拟是一个不错的选择。

《基于ABAQUS钢结构复杂节点分析9讲—获得钢结构节点FEA分析流程和核心关键技能》就是笔者利用ABAQUS软件对钢结构节点FEA流程和技巧进行的详细讲解，希望对学习型仿真工程师有些许帮助。

以下是大纲安排（点击下图可以试看）：

1、课程适应人群

学习型仿真工程师
理工科院校学生
从事钢结构设计及分析的结构工程师
不满足于专用结构分析软件，想提升有限元分析水平的资深结构师

2、课程学习后的效果

理解钢结构节点的设计分析流程；
掌握常见钢结构节点的建模、荷载施加、结果校核；
掌握实体单元的网格划分；
掌握平面桁架屋盖结构的节点设计；
掌握钢管桁架结构的焊接节点设计；
掌握钢框架的连接节点设计
掌握空间网架结构的连接节点设计
掌握门式刚架结构连接节点设计
掌握通用软件ABAQUS在节点分析中的应用
订阅用户可以获得讲师的答疑专栏（可公开提问或私问）
购买全套课程的用户可以获得全课程模型，ppt 等资料下载服务。

所有订阅用户可以获得讲师的90天答疑专栏服务（可公开提问或私问），并且提供VIP群交流，后续可以根据用户需求开展免费的加餐直播。购买全套课程的用户可以获得全课程模型，PDF等资料下载服务。进群和开具发票请联系平台任意仿真小助手。

（完）

作者：赵老师，11年的结构分析经验，6年大型设计院实战经验，工学硕士学历，仿真秀科普作者。曾就职于北京市建筑设计院、中国市政西南设计院研究总院。项目经历：湖南岳阳大酒店超限设计报告。

声明：本文首发仿真秀App，部分图片和内容源自网络，如有不当请联系我们，欢迎分享，禁止私自转载，转载请联系我们。

喜欢作者，请点赞和在看

道桥岩土建筑结构基础 Abaqus

著作权归作者所有，欢迎分享，未经许可，不得转载

首次发布时间：2021-06-02