基于BIM数据的建筑抗震弹塑性分析建模技术

曾翔,陆新征,许镇

(清华大学土木工程系,北京市  100084

计算机辅助工程/Computer Aided Engineering, 2014, 23(4): 5-11.

PDF 全文下载

摘要:利用建筑信息模型(BIM)数据自动生成结构分析模型可提高建模效率,BIM的共享数据标准为IFC标准,但现有软件导出的IFC格式模型缺少关键数据,不能满足结构非线性抗震分析的数据需求。因此,本研究在深入了解IFC标准的基础上,建立了建筑抗震弹塑性分析数据基于IFC标准的表达方法,从而有望将BIM的共享性应用于结构分析领域。并编写了C++转换程序,将BIM包含的数据信息正确地转换成不同软件(如MarcOpenSees)的模型,从而可以使用不同的结构分析软件实现建筑抗震弹塑性分析。该研究成果为基于BIM数据实现建筑抗震弹塑性分析提供了参考。

关键词: IFC标准;MarcOpenSees;前处理;建筑抗震分析;纤维模型;结构信息转换;IFC子信息模型

 

Modeling Technology of Structural Seismic Elasto-plastic Analysis Based on BIM Data

ZENG Xiang, LU Xinzheng, XU Zhen

(Department of Civil Eng., Tsinghua Univ., 100084, Beijing, China)

Abstract: It is more efficient if the structural analytical model is generated automatically from BIM (Building Information Model), which uses IFC (Industry Foundation Classes) as the data sharing standard. However, IFC structural analysis models exported by conventional used software lack some key information, thus cannot satisfy the data requirements for nonlinear structural seismic analysis. Therefore, on the basis of a deep understanding of IFC standard, the corresponding expression of those data in the IFC standard is established for the application of BIM in structural analysis domain. A C++ code is developed to correctly convert the BIM data into different structural analytical models (e.g. Marc, OpenSees, etc.), so that the seismic analysis can be performed based on different kinds of structural analysis software. The research serves as a reference for structural seismic elasto-plastic analysis based on BIM data.

Key words: IFC standard; Marc; OpenSees; Pre-processing; structural seismic analysis; fiber model; conversion of structural information; IFC sub-model


 

随着城市化的发展,地震对经济发展和人民生活的潜在威胁越来越严重,城市建筑震害预测成为防灾减灾的一个重要方面。城市区域内的建筑群震害预测包含两个层次:宏观层次和微观层次。宏观层次模型[1]将建筑抽象成一个或少量自由度的结构,进行不同烈度下地震损失分析预测。微观层次模型[2]将建立建筑物更加精细化的分析模型(例如杆系模型等),可以考虑结构的足够多的自由度和精确的动力特性,进行高精度的震害分析预测。

城市建筑震害预测,需要获取和存储城市建筑的基础数据。如果只是建立宏观层次模型,则只需要建筑面积、年代、层数、结构类型等少量宏观数据,其数据来源可以是城市GIS数据库等。基于GIS的城市震害预测已经有一定的研究基础[3]。而建立微观层次的精细化模型需要更详细的数据,如梁柱的尺寸、配筋等,若采用手动方法建立精细化模型,不仅工作量太大,也容易出错,特别对于某些开源的有限元软件(如OpenSees),没有前处理界面,建模极不方便。比较合理的方法是利用建筑已有的数字化模型,从中提取建筑震害分析所需要的详细数据,自动建模供计算模块进行分析计算。然而不同软件生成的建筑数字化模型有着不同的数据格式,而计算模块也可能由不同的结构分析软件组成(如MarcOpenSees等)。因此,需要采用一个统一的数据标准,以实现不同软件之间的数据交换。

建筑信息模型(Building Information ModelBIM)致力于在建筑的全生命期实现不同建筑工程领域的信息共享和协同管理,基于BIM的建筑数字化模型是未来重要的发展方向。目前,BIM采用IFC标准作为建筑产品的数据交换标准,不同软件之间通过IFC这个统一的标准来交换数据。然而研究表明[4],基于IFC标准的建筑模型与结构模型之间的信息传递与共享还不是十分成熟。张剑涛等[5],刘照球[6, 7]等开发了IFCPKPM结构模型的数据转换接口,能初步实现二者的几何、材料等信息的转换,但缺乏载荷信息等的转换。也有学者尝试建立通用建筑结构模型转换平台,以IFC为通用数据转换标准,实现任意常用结构分析软件间的结构分析模型数据转换。但这类研究尚处于起步阶段,目前仅实现了SATWEETABSSAP2000等软件与IFC标准的转换接口,所能实现的数据转换也仅限于几何、材料等信息。

因此,能否利用BIM的共享性,解决城市建筑震害预测的数据获取以及自动建模问题,具有重要的研究价值。

建筑震害分析及数据需求

对于钢筋混凝土框架结构,通常采用纤维梁模型作为抗震弹塑性分析的精细化计算模型。不少软件都能具备建筑抗震弹塑性分析能力,其中两类典型的软件是MarcOpenSees

由于城市建筑信息是以单体建筑为基础的,因此,本文利用C++编写了一个模型转换程序,从一个单体钢筋混凝土建筑的(IFC格式的)数字化模型中,提取建筑震害分析所需的信息,自动生成MarcOpenSees软件模型,验证转换后的模型能否用于分析建筑结构在地震作用下的响应,如图1所示。


图1 编写模型转换程序,将IFC格式模型转换成Marc模型和OpenSees模型

图1  编写模型转换程序,将IFC格式模型转换成Marc模型和OpenSees模型

Fig.1  Convert IFC-format files into Marc and OpenSees models using a model converter


1.1  MarcOpenSees软件

Marc是一款功能齐全的高级非线性通用有限元软件[8],用Marc进行建筑震害模拟,可以得到较为理想的分析结果,且适用于大多数分析需求。

基于纤维模型原理,清华大学土木工程系利用Marc软件提供的二次开发接口,编写了THUFIBER程序[9],使Marc软件能对钢筋混凝土结构进行纤维模型计算。THUFIBER将混凝土构件的截面分为36个混凝土纤维,以及4个钢筋纤维(纤维划分方法也可以由用户自定义),各纤维的面积由截面尺寸和配筋决定。对每根纤维定义特定的本构关系,便可以进行不同精度的有限元分析。

OpenSeesOpen System for Earthquake Engineering Simulation),是一个能进行结构抗震弹塑性分析的开源软件[10]。相比内部代码不透明的商用有限元软件,OpenSees的开源性使得研究者能全面掌握它的运行机制,并能根据自身需求进行二次开发,因此,OpenSees是进行建筑震害研究的重要工具。OpenSees中,混凝土梁柱也用纤维模型模拟。

选用上述两种软件进行震害分析,可以同时利用商业软件的健壮性和开源软件的透明性。

1.2  建筑震害分析的数据需求

为了从IFC文件中提取充分的信息,以自动生成MarcOpenSees结构模型,需要全面了解这两种软件的建筑震害分析建模方法,并总结震害分析所需的数据信息。这些数据信息分为两类,一类与混凝土结构模型有关,如结构构件的位置、几何信息、材料信息、载荷等。另一类则与震害分析的方法和过程有关,如地震波的选取、载荷工况的增量步总数和每个增量步的时间步长等,这类数据随着分析方法和要求的变化而变化,并且在软件中可以很方便地实现这类参数的调整。因此,应该重点关注第一类数据需求,即与混凝土结构模型有关的数据。

通过考察这两种软件的建筑抗震分析模型,得到抗震分析的数据需求如表1所示。

表1  建筑抗震分析的数据需求

Tab.1  Data requirements for structural seismic analysis

项目

震害分析的数据需求

单元

元编号、总数、类型,单元所包含的节点。

节点

节点编号,节点总数,节点坐标。

钢筋混凝土材料

(核心区和保护层的)混凝土弹性模量、混凝土的应力应变关系,混凝土保护层厚度,截面配筋情况,钢筋弹性模量、屈服强度,配筋截面对应的单元编号。

截面几何性质

截面尺寸,截面局部坐标x轴在整体坐标系下的方向,截面几何性质对应的单元编号。

支座约束

支座自由度约束状况,支座约束名称,支座约束对应的单元编号。

载荷

恒、活载荷、其他载荷,载荷对应的单元编号。

抗震分析的数据需求的IFC表达

2.1  IFC标准概述

IFCIndustry Foundation Classes)标准是针对建筑信息模型的开源的产品模型标准,最新版本为IFC4[11],于20133月正式发布。

IFC4包含类型、实体、函数、规则、(预定义)属性集、量集等元素。其中,实体对应于面向对象思想中的“类”, 每个实体都拥有一条或多条属性,对应于“类”中的数据成员。然而IFC体系中,为了更灵活、全面地表达实体的属性,引入了“属性集”概念[12],即属性的集合。通过关联实体的关联机制,将某个属性集和某个实体相关联,就实现了对该实体的属性的扩展。

IFC4标准既能描述建筑模型,又能描述结构分析模型,但这两类模型的内容有很大差异。例如同一根梁,在IFC4建筑模型中由IfcBeam实体表示,是由一个封闭截面拉伸而形成的实体(或者其他实体生成方式),并且具有颜色、渲染材质等属性;而在IFC4结构分析模型中由IfcStructuralCurveMember表示,是由两个顶点和之间的连线组成的实体。图2是一个建筑模型和对应结构分析模型的对比示意图[13]

图2 建筑模型和结构分析模型的对比

图2  建筑模型和结构分析模型的对比

Fig.2  Comparison between architectural model and corresponding structural analysis model

2.2  IFC标准在结构分析中的应用

目前,IFC在建筑设计和施工等领域的研究和应用十分丰富,而在结构分析领域相对缺乏。

Wan[14]通过对比结构分析软件SAP2000所需信息和IFC2x2标准所能描述的信息之间的差异,分析了IFC标准对结构分析模型的支持程度,发现IFC可以提供简单的静力载荷、静力工况、部分材料力学属性等数据,但无法描述弯矩释放、材料剪切性能等部分信息。

虽然IFC标准本身具备描述基本结构分析模型的能力,但其应用并不成熟,表2是常用建筑设计软件和结构分析软件的IFC模型导出效果,导出的IFC文件最高版本为IFC2x3。从表中可以看出,现有软件的IFC结构分析模型导出能力较差,不能满足结构抗震弹塑性分析的要求。

表2  部分软件的IFC文件导出能力测试

Tab.2 IFC file exporting test for conventional used software

测试软件

几何属性

材料属性

结构分析模型实体

配筋信息

载荷

支座约束

Revit Structure 2013

×

×

×

×

Etabs 9

×

×

×

×

ArchiCAD 16

×

×

×

×

SAP2000 v15

×

×

×

×

×

Tekla Structures 18.1

×

×

×

本文建立的IFC4结构分析子信息模型

2.3  IFC4结构分析子信息模型

由于现有常用软件导出的IFC文件中缺少建筑结构抗震弹塑性分析的关键数据,因此有必要建立一个完整的IFC结构分析模型,如表2最后一行所示。子信息模型是指整个IFC标准体系的一个子集。作为一个建筑产品通用数据标准,IFC提供了丰富的建筑描述方法,以适应各种描述需求。但同时也导致IFC标准非常复杂,例如描述一根梁的配筋至少有四种方法:(1)用IfcMaterialLayer实体定义层状材料,将一个构件的材料分成若干钢筋层、混凝土层,这种方法更适合表达剪力墙、楼板等的配筋。(2)用IfcReinforcementDefinitionProperties实体定义配筋,这种方法通常用于早期设计阶段。(3)用IfcReinforcingBar实体描述钢筋,并通过关系实体将钢筋与梁相关联,然而钢筋的位置既可能用全局坐标表达,又可能用局部坐标表达;既可以为每根钢筋定义一个IfcReinforcingBar实体,又可以仅定义一根钢筋,再用映射方法IfcMappedItem生成任意数量的相同类型钢筋;既可以为每根钢筋指定屈服强度等属性,又可以为IfcReinforcingBarType指定属性从而一次性给所有同类型钢筋指定了该属性。因此使用这种方法定义配筋,可以描述大多数类型的配筋方式,但若要从中提取配筋信息,实现起来将非常复杂。(4)用自定义属性集的方式为梁添加配筋属性,它表达灵活,也便于编程提取信息。

IFC的复杂性,也是限制它发展、导致常用软件普遍无法导出满足抗震分析的IFC文件的原因之一。因此本文在仔细考察IFC4体系的基础之上,选择其中的常用部分,建立IFC4结构分析子信息模型。例如配筋的描述仅保留了方法(4)。虽然这降低了IFC的全面性,但仍能描述抗震分析所需的关键信息,增强了实用性。

描述IFC(子信息)模型的最方便的工具是EXPRESS-G[15],图3是它的部分图例。

图3 EXPRESS-G部分图例

图3  EXPRESS-G部分图例

Fig.3  Legend for EXPRESS-G (partial)

2.3.1 IFC4结构分析模型的描述

任何一个IFC4模型文件中,有且仅有一个IfcProject实体。因此,任何IFC子信息模型中都包含IfcProject实体,用于表征一个建筑项目。IFC4中,用IfcStructuralAnalysisModel实体表示一个结构分析模型,该实体有结构分析类型、全局坐标系、所包含的载荷工况等属性。载荷工况由IfcStructuralLoadGroup实体的子类,即IfcStructuralLoadCase表达,该实体的主要属性有载荷类型(永久载荷、可变载荷、偶然载荷等)与载荷模式(恒载荷、活载荷、风载荷等)。一个IfcProject可对应几个IfcStructuralAnalysisModel,二者之间通过IfcRelDeclares关系实体建立关联关系,其EXPRESS-G表达如图4所示,图中实体属性未全部列出,黄色标出的是主要实体。

图4 IFC4结构分析模型的描述方法

图4  IFC4结构分析模型的描述方法

Fig.4  Description of IFC4 structural analysis model

2.3.2其他建筑抗震分析数据的IFC4描述方法

2.3.1小节类似地,建立了结构构件、载荷、支座约束、材料、截面几何特性等信息的IFC4描述方法。以载荷为例,如图5所示,图中实体属性未全部列出。IFC4能描述线性变化的分布载荷,但为了降低复杂性,本文建议将这类载荷转换成等效均布载荷,用IfcStructuralLinearAction实体表达。


图5 载荷信息在IFC4中的描述方法

图5  载荷信息在IFC4中的描述方法

Fig.5  Description of load in IFC4 standard


模型转换程序设计

根据第2章所建立的完整的IFC4结构分析子信息模型框架,可以构建满足抗震分析数据需求的IFC4模型文件,进而采用C++语言设计模型转换程序,解析IFC模型文件,自动生成MarcOpenSees软件模型。

IFC模型由EXPRESS语言定义,是ASCII格式纯文本文件形式,因此IFC文件需要经过解析识别,才能供程序代码处理。本文使用IFCEngine.dll作为IFC文件解析工具,这是一款免费工具包,但要求编程者对IFC体系十分熟悉。

根据表1所示建筑震害的数据需求,设计了节点、单元、支座约束、截面几何、材料等5个类(class),并用MapVectorC++标准模板库中的容器将各个类的不同对象有序地组织起来。程序的输入模块通过调用类的成员函数,从IFC模型文件中提取所需数据,赋值给相应的数据成员。例如节点类Node调用SetCoordinate函数提取并设置节点坐标。模型转换程序的IFC解析及提取数据部分的流程如图6所示。


图6 模型转换程序的IFC解析及提取数据部分的流程图

图6  模型转换程序的IFC解析及提取数据部分的流程图

Fig.6  Flow chart of the IFC parse and data obtaining part of the model converter


由于采用了合理的数据结构设计,程序的输出模块的实现比较方便,分别遍历相应容器,按照MarcOpenSees的建模规则输出节点、单元、材料等类的数据成员,便可自动生成这两种软件模型。

实例验证模型转换效果

为了验证模型转换程序的正确性,用IFC4标准建立了一个钢筋混凝土框架模型。模型为清华大学主楼的一榀103跨平面框架[16],边跨7.5米,中跨3.5米,框架总高度45.05米,模型详细尺寸及配筋等信息见文献[16]IFC模型以及转换后的Marc模型如图7所示(OpenSees没有显示模型的界面),由图可知单元、节点、支座约束转换正确。检查模型参数,截面、材料、载荷等也都转换正确。

图7 IFC模型(左)及自动生成的Marc模型(右)

图7  IFC模型(左)及自动生成的Marc模型(右)

Fig.7  IFC model(left) and auto-generated Marc model(right)

分别用MarcOpenSees对自动生成的模型进行模态分析,得到的前三阶振型所对的频率相差小于1%,如表3所示:

表3  Marc与OpenSees模态分析结果对比

Tab.3 Modal analysis result of Marc and OpenSees

分析软件

1阶圆频率

2阶圆频率

3阶圆频率

基本周期/s

Marc

4.07

11.48

21.69

1.53

OpenSees

4.11

11.57

21.82

1.54

对框架输入Y方向地震动加速度,分别用MarcOpenSees进行弹塑性时程分析。选用El Centro地震波,地震峰值加速度设为400Gal,顶点位移的时程曲线如图8所示,吻合得很好。

图8 Marc与OpenSees弹塑性时程分析结果对比

图8  Marc与OpenSees弹塑性时程分析结果对比

Fig.8  Elasto-plastic time-history analysis result of Marc and OpenSees

两个软件的模态分析和时程分析结果的一致性,进一步说明了模型的单元、材料等各参数都得到了正确转换。

 

本文从建筑震害分析的数据需求出发,逐一分析了这些数据需求在IFC4中的表达方法,建立了一个完整的IFC4结构分析子信息模型。利用该子信息模型,就能通过IFC4标准来描述可进行建筑抗震弹塑性分析的钢筋混凝土结构模型。

根据建立的IFC4结构分析子信息模型以及MarcOpenSees软件的建模规则,编写了模型转换程序,实现了节点、构件、载荷、支座约束、材料、截面几何等所有抗震分析所必要的信息的正确转换,进行了IFC文件到MarcOpenSees软件之间的数据转换的新尝试。

在此基础上,还有很多研究内容可以继续开展。例如,扩展本文建立的IFC4结构分析子信息模型,使之能描述楼板、剪力墙等其他结构构件;考虑城市建筑群信息的IFC表达方法等。


参考文献:

[1] 韩博, 陆新征,许镇, . 基于高性能GPU计算的城市建筑群震害模拟[J]. 自然灾害学报, 2012,21(5):16-22.

    HAN Bo, LU Xinzheng, XU Zhen, et al. Seismic damage simulation of urban buildings based on high performance GPU computing[J]. Journal of Natural Disasters, 2012,21(5):16-22.

[2] 卢啸, 陆新征, 张万开, . 特大地震下超高层建筑的倒塌模拟[J]. 中国科学:技术科学, 2011(11):1405-1416.

    LU Xiao, LU Xinzheng, ZHANG Wankai, et al. Collapse simulation of a super high-rise building subjected to extremely strong earthquakes. Sci China Tech Sci, 2011, 54: 2549-2560.

[3] XU Feng, CHEN Xuping, REN Aizhu, et al. Earthquake disaster simulation for an urban area, with GIS, CAD,FEA, and VR integration[J]. Tsinghua Science and Technology, 2008,13(S1):311-316.

[4] Applied Technology Council. Project ATC-75[EB/OL]. [2012.11.12]. https://www.atcouncil.org/Projects/atc-75-project.html.

[5] 张剑涛, 刁波, 唐春风, . IFC标准在PKPM结构软件中的实现[J]. 建筑科学, 2006,22(4):103-106.

    ZHANG Jiantao, DIAO Bo, TANG Chunfeng, et al. The exchange of information between IFC and PKPM[J]. Building Science, 2006,22(4):103-106.

[6]   刘照球, 李云贵, 吕西林, . 建筑结构信息集成的程序实现[J]. 沈阳建筑大学学报(自然科学版), 2009,25(3):467-473.

    LIU Zhaoqiu, LI Yungui, LV Xilin, et al. Program realization for integrated information in architectural and structural design field[J]. Journal of Shenyang Jianzhu University ( Natural Science), 2009,25(3):467-473.

[7]   刘照球, 李云贵, 吕西林, . 基于BIM 建筑结构设计模型集成框架应用开发[J]. 同济大学学报(自然科学版), 2010,38(7):948-953.

    LIU Zhaoqiu, LI Yungui, LV Xilin, et al. BIM-based integrated information framework for architectural and structural design model[J]. Journal of Tongji University ( Natural Science), 2010,38(7):948-953.

[8]   蒋学武, 吴新跃, 朱石坚. 综合应用UG, HyperMeshMSC Marc软件进行有限元分析[J]. 计算机辅助工程, 2007(02):11-14.

JIANG Xuewu, WU Xinyue, ZHU Shijian. Comprehensive application of UG, HyperMesh and MSC.Marc in finite element analysis[J]. Computer Aided Engineering, 2007(02):11-14.

[9]  陆新征, 叶列平, 缪志伟, . 建筑抗震弹塑性分析——原理、模型与在ABAQUSMSC.MARCSAP2000上的实践[M]. 北京: 中国建筑工业出版社, 2009: 143-162.

    LU Xinzheng, YE Lieping, MIU Zhiwei, et al. Elasto-plastic analysis of building against earthquaketheory, model and implementation on ABAQUS, MSC.MARC and SAP2000[M]. Beijing: China Architecture & Building Press, 2009: 143-162.

[10] McKenna F. OpenSees: A Framework for Earthquake Engineering Simulation[J]. COMPUTING IN SCIENCE & ENGINEERING, 2011,13(4):58-66.

[11]  BuildingSMART. Industry Foundation Classes IFC4 Official Release [EB/OL]. [2013.03.20]. http://www.buildingsmart-tech.org/ifc/IFC4/final/html/.

[12] 张洋. 基于BIM的建筑工程信息集成与管理研究[D]. 北京: 清华大学, 2009.

    ZHANG Yang. Research on BIM-based building information integration and management[D]. Beijing: Tsinghua University, 2009.

[13] 秦领, 刘西拉. 建筑物理模型与结构分析模型的数据映射研究[J]. 土木建筑工程信息技术, 2010,2(2):28-36.

    QIN Lin, LIU Xila. Research on data mapping between architectural physical models and structural analysis models[J]. Journal of Information Technology in Civil Engineering and Archietecture, 2010,2(2):28-36.

[14] Wan Caiyun, Chen Po-Han, Tiong Robert L. K. Assessment of IFCS for Structural Analysis Domain[J]. Information Technology in Construction, 2004,9:75-95.

[15] Data Modelling Using EXPRESS-G for IFC Development[EB/OL]. http://www.stanford.edu/group/narratives/classes/08-09/CEE215/ReferenceLibrary/Industry%20Foundation%20Classes%20(IFC)/IFC%20General/Data_Modelling_Using_EXPRESS-G_for_IFC_Development.pdf.

[16] 邹积麟. 空间RC框架结构全过程静力弹塑性分析[D]. 北京: 清华大学, 2001.

    ZOU Jilin. 3D RC frame structure full-range pushover analysis[D]. Beijing: Tsinghua University, 2001.

个人信息
研究工作
实际工程
论文工作
教学工作
资料下载
专题
其他

 

我们的实验室

抗倒塌专业委员会