【JY】力荐 | 区域建筑地震安全性有限元分析示例
一、研究背景
城市重要建筑的地震安全水平及韧性水平由于现代抗震及减震控制技术的大量应用正在获得实质性提高，但其水平到底如何，最终提高到什么水平才能满足经济发展对于地震安全性的诉求，需要一个客观的评价标准和技术体系。
因此，基于韧性的抗震设计可以进一步指导工程结构的利益相关者采取更加合适的措施，去改进或提高其韧性能力，或者采取其它的防灾备灾预案，从而达到更加有效的减灾目标。而本实例主要从多尺度下进行抛砖引玉，给出区域建筑地震仿真分析的实例。本贴仅探究城市区域性建筑的建筑结构系统抗震分析，暂不讨论地震后城市引发的次生灾害及非结构构件等其他评估。旨在对于区域地震进行可用于韧性评估较准确的精细化分析。
由于土木工程的研究对象是材料、构件、结构、建筑、社区、城市等多尺度的范围，地震灾害会对各类受灾体造成损伤，因此无论是小到一种材料、一个个体，还是大到一座城市区域，都存在韧性的问题。从多尺度角度进行问题的拆分层次看，将区域建筑拆分成单体建筑，进而拆分成构件，进而研究材料本构。本文的验证研究思路将从本构模型验证到构件验证，进而到单体建筑最后到区域建筑模拟分析。
图 1‑1 城市建筑地震安全性分析层次2.1材料模型
当遭遇强烈地震作用时，结构将进入非线性阶段，其材料特性发生较大改变，只进行线弹性阶段的分析会造成较大误差，因此在进行建筑结构时程分析中应考虑材料的非线性。
ABAQUS中，软件自带的有两种常用的适用于混凝土的本构材料模型，分别为混凝土弥散裂缝模型（Concrete Smeared Cracking）和混凝土塑性损伤模型（Concrete Damaged Plastici-ty），前者以裂缝模型为基础，一般适用于带有钢筋特性的混凝土，仅适用于单调加载分析；后者以损伤模型为基础，考虑了损伤效应，可用于往复荷载作用分析，适用于模拟地震工况下的混凝土力学行为。
在使用ABAQUS结合《混凝土结构设计规范》GB50010-2010对结构进行弹塑性分析时，其中规范中的“损伤演化参数”Dc与ABAQUS中的“损伤因子”dc并不是一个概念，Dc为应力应变曲线上的割线损伤，而dc为卸载刚度损伤。因而在根据规范中的本构模型计算输入ABAQUS塑性损伤参数的时候，要进行相应的转换。
一维混凝土材料模型采用规范指定的单轴本构模型,能反应混凝土滞回、刚度退化和强度退化等特性，其轴心抗压和轴心抗拉强度标准值按《混凝土结构设计规范》表4.1.3采用。混凝土单轴受拉的应力-应变曲线方程按附录C公式C.2.3-1~C.2.3-4计算。
混凝土材料进入塑性状态伴随着刚度的降低。如应力-应变及损伤示意图所示，其刚度损伤分别由受拉损伤参数dt和受压损伤参数dc来表达,dt和dc由混凝土材料进入塑性状态的程度决定。
二维混凝土本构模型采用弹塑性损伤模型，该模型能够考虑混凝土材料拉压强度差异、刚度及强度退化以及拉压循环裂缝闭合呈现的刚度恢复等性质。
当荷载从受拉变为受压时,混凝土材料的裂缝闭合,抗压刚度恢复至原有抗压刚度;当荷载从受压变为受拉时,混凝土的抗拉刚度不恢复，如下图所示。
关于混凝土损伤模型及Abaqus-CDP参数详解，可参看：
【JY】浅谈混凝土损伤模型及Abaqus中CDP的应用
【JY】ABAQUS混凝土CDP插件分享
【JY】混凝土分析工具箱：CDP模型插件与滞回曲线数据
由于混凝土的CDP模型无法在杆系单元（B31/B32）中使用，所以这里采用了Kent-Park的本构关系，利用Fortran将该本构写入UMAT进行验证分析，并将UMAT的本构修改成VUMAT格式，对模型进行显式动力弹性分析。Kent-Park的本构关系其考虑了箍筋作用对混凝土强度和应变的提高及不同加载速率对应力-应变关系的影响，数学表达式如下：
式中fc、εc分别为任一点的单轴压应力和相应的应力；f'c为圆柱体抗压强度；ρs为体积配箍率，fyh为箍筋的屈服强度；h''为箍筋约束核心区混凝土的高度；sh为箍筋间距。
ABAQUS软件自带的适用于动力分析下的钢材模型有等向强化模型与随动强化模型。等向强化即假设加载曲面是屈服曲面在形状、中心位置均不变的情况下做等向膨胀得到的；随动强化则假设加载曲面、屈服曲面只是发生中心位置的移动，两个曲面大小、形状保持一致。通过大量实验表明，材料的加载曲面与屈服曲面实际上既有中心位置的移动，也有曲面大小的改变，即介于等向强化和随动强化之间。
本次模型中，钢材的本构选用随动强化模型，需要弹性模量、屈服强度以确定本构关系。该模型可以考虑包辛格效应，能够较好地描述钢材在往复作用下的弹塑性发展，并有较高的计算效率。钢材的强屈比设定为1.2，极限应力所对应的极限塑性应变为0.025。
二维墙板采用分层壳单元，钢筋层采用钢材本构，混凝土采用CDP，内置本构已经过大量验证，这里不再验证。一维杆系单元采用自主开发UMAT/VUMAT子程序进行计算分析，为验证本构正确性，取清华大学钢筋混凝土框架及关键构件试验数据库中的混凝土框架柱来进行验证模型分析。
由图 2-10可知，模拟计算分析与试验的峰值与趋势吻合较好。模拟时，曲线负向向正向加载时刚度略大，原因是采用的钢筋本构为双折线，而钢筋存在软化现象。所以，将以上材料、构件模型综合起来即可建立超高层建筑与城市防灾仿真。
对于城市防灾（地震）的分析，对结构进行的弹塑性分析将采用基于显式积分的动力弹塑性分析方法，这种分析方法未作任何理论的简化，直接模拟结构在地震力作用下的非线性反应，具有如下优越性：
目前，我国超高层建筑工程的建设规模已位居世界前列。超高层建筑体型优美，同时能为人们提供舒适办公和生活环境。由于其层数多，使用面积大，容纳人数较多，分析超高层结构在罕遇地震下的弹塑性响应，对于抗震设计有重要意义,并能实现城市建筑地震安全性水平及韧性的科学评价。
由于该分析的数据量较多，本文采用Matlab Python的形式辅助Abaqus进行一系列分析及数据提取。
本工程为某超高层写字楼，地面以上33层，标准层层高，采用现浇钢筋混凝土主次梁楼盖形式，主次梁和板的混凝土强度为C30，柱和剪力墙的混凝土强度为C45，有限元模型见图 3-1。
图3-2、3-3分别给出了所输入的El波地震动加速度时程曲线、加速度反应谱。根据GB 50011-2010《建筑抗震设计规范》(2016年版)的规定，截取该地震波的有效波长前30s进行计算分析。修正X 0.85Y输入烈度为8度罕遇400Gal。
关于SignalData软件及使用教程，可参看：
在该地震波情况下最薄弱的地方是5号楼层，且XY方向的最大层间位移角为1/91和1/34，最大顶点位移XY分别为0.144m和0.323m。
由下图可得到，在该地震波情况下，X向的基地剪力为128.3MN，Y向的基地剪力为97.5MN。
由上可验证得到材料→构件→单体建筑的分析验证，进一步的将构建城市区域建筑进行抛砖引玉的浅析。建立的区域建筑群模型如下，仍采用El波8度罕遇地震进行输入。
本文基于有限元软件Abaqus，建立了一套从材料、构件、建筑再到城市的地震安全性的有限元技术体系。通过与既有试验数据进行对比，表明本文的本构模型具有较好的精度。并建立某超高层结构以及城市的动力弹塑性分析模型，通过施加了El波，研究了某超高层结构以及区域建筑在遭受罕遇地震作用下的弹塑性时程响应。
结果表明：该方法可用于建筑地震安全性水平及结构层面水平的分析。基于Abaqus为总结建筑地震韧性的评价流程做一个铺垫分析方法，对后续重要建筑地震韧性评价研究给了方法，为后续建筑地震韧性的定义及等级划分标准，以及对地震韧性的相关评价指标，建立损伤指数与功能损失之间的关系，做好铺垫分析。
