引  言

    隨著城市化的發(fā)展，地震對(duì)經(jīng)濟(jì)發(fā)展和人民生活的潛在威脅越來越嚴(yán)重，城市建筑震害預(yù)測(cè)成為防災(zāi)減災(zāi)的一個(gè)重要方面。城市區(qū)域內(nèi)的建筑群震害預(yù)測(cè)包含兩個(gè)層次：宏觀層次和微觀層次。宏觀層次模型將建筑抽象成一個(gè)或少量自由度的結(jié)構(gòu)，進(jìn)行不同烈度下地震損失分析預(yù)測(cè)。微觀層次模型[2]將建立建筑物更加精細(xì)化的分析模型（例如桿系模型等），可以考慮結(jié)構(gòu)的足夠多的自由度和精確的動(dòng)力特性，進(jìn)行高精度的震害分析預(yù)測(cè)。

    城市建筑震害預(yù)測(cè)，需要獲取和存儲(chǔ)城市建筑的基礎(chǔ)數(shù)據(jù)。如果只是建立宏觀層次模型，則只需要建筑面積、年代、層數(shù)、結(jié)構(gòu)類型等少量宏觀數(shù)據(jù)，其數(shù)據(jù)來源可以是城市GIS數(shù)據(jù)庫(kù)等?；贕IS的城市震害預(yù)測(cè)已經(jīng)有一定的研究基礎(chǔ)[3]。而建立微觀層次的精細(xì)化模型需要更詳細(xì)的數(shù)據(jù)，如梁柱的尺寸、配筋等，若采用手動(dòng)方法建立精細(xì)化模型，不僅工作量太大，也容易出錯(cuò)，特別對(duì)于某些開源的有限元軟件（如OpenSees），沒有前處理界面，建模極不方便。比較合理的方法是利用建筑已有的數(shù)字化模型，從中提取建筑震害分析所需要的詳細(xì)數(shù)據(jù)，自動(dòng)建模供計(jì)算模塊進(jìn)行分析計(jì)算。然而不同軟件生成的建筑數(shù)字化模型有著不同的數(shù)據(jù)格式，而計(jì)算模塊也可能由不同的結(jié)構(gòu)分析軟件組成（如Marc，OpenSees等）。因此，需要采用一個(gè)統(tǒng)一的數(shù)據(jù)標(biāo)準(zhǔn)，以實(shí)現(xiàn)不同軟件之間的數(shù)據(jù)交換。

    建筑信息模型（Building Information Model，BIM）致力于在建筑的全生命期實(shí)現(xiàn)不同建筑工程領(lǐng)域的信息共享和協(xié)同管理，基于BIM的建筑數(shù)字化模型是未來重要的發(fā)展方向。目前，BIM采用IFC標(biāo)準(zhǔn)作為建筑產(chǎn)品的數(shù)據(jù)交換標(biāo)準(zhǔn)，不同軟件之間通過IFC這個(gè)統(tǒng)一的標(biāo)準(zhǔn)來交換數(shù)據(jù)。然而研究表明[4]，基于IFC標(biāo)準(zhǔn)的建筑模型與結(jié)構(gòu)模型之間的信息傳遞與共享還不是十分成熟。張劍濤等[5]，劉照球[6, 7]等開發(fā)了IFC與PKPM結(jié)構(gòu)模型的數(shù)據(jù)轉(zhuǎn)換接口，能初步實(shí)現(xiàn)二者的幾何、材料等信息的轉(zhuǎn)換，但缺乏載荷信息等的轉(zhuǎn)換。也有學(xué)者嘗試建立通用建筑結(jié)構(gòu)模型轉(zhuǎn)換平臺(tái)，以IFC為通用數(shù)據(jù)轉(zhuǎn)換標(biāo)準(zhǔn)，實(shí)現(xiàn)任意常用結(jié)構(gòu)分析軟件間的結(jié)構(gòu)分析模型數(shù)據(jù)轉(zhuǎn)換。但這類研究尚處于起步階段，目前僅實(shí)現(xiàn)了SATWE、ETABS、SAP2000等軟件與IFC標(biāo)準(zhǔn)的轉(zhuǎn)換接口，所能實(shí)現(xiàn)的數(shù)據(jù)轉(zhuǎn)換也僅限于幾何、材料等信息。

    因此，能否利用BIM的共享性，解決城市建筑震害預(yù)測(cè)的數(shù)據(jù)獲取以及自動(dòng)建模問題，具有重要的研究?jī)r(jià)值。

1. 建筑震害分析及數(shù)據(jù)需求

    對(duì)于鋼筋混凝土框架結(jié)構(gòu)，通常采用纖維梁模型作為抗震彈塑性分析的精細(xì)化計(jì)算模型。不少軟件都能具備建筑抗震彈塑性分析能力，其中兩類典型的軟件是Marc和OpenSees。

    由于城市建筑信息是以單體建筑為基礎(chǔ)的，因此，本文利用C++編寫了一個(gè)模型轉(zhuǎn)換程序，從一個(gè)單體鋼筋混凝土建筑的（IFC格式的）數(shù)字化模型中，提取建筑震害分析所需的信息，自動(dòng)生成Marc和OpenSees軟件模型，驗(yàn)證轉(zhuǎn)換后的模型能否用于分析建筑結(jié)構(gòu)在地震作用下的響應(yīng)，如圖1所示。

 

圖1 編寫模型轉(zhuǎn)換程序，將IFC格式模型轉(zhuǎn)換成Marc模型和OpenSees模型
Fig.1  Convert IFC-format files into Marc and OpenSees models using a model converter

1.1 Marc與OpenSees軟件

    Marc是一款功能齊全的高級(jí)非線性通用有限元軟件[8]，用Marc進(jìn)行建筑震害模擬，可以得到較為理想的分析結(jié)果，且適用于大多數(shù)分析需求。

    基于纖維模型原理，清華大學(xué)土木工程系利用Marc軟件提供的二次開發(fā)接口，編寫了THUFIBER程序[9]，使Marc軟件能對(duì)鋼筋混凝土結(jié)構(gòu)進(jìn)行纖維模型計(jì)算。THUFIBER將混凝土構(gòu)件的截面分為36個(gè)混凝土纖維，以及4個(gè)鋼筋纖維（纖維劃分方法也可以由用戶自定義），各纖維的面積由截面尺寸和配筋決定。對(duì)每根纖維定義特定的本構(gòu)關(guān)系，便可以進(jìn)行不同精度的有限元分析。

    OpenSees（Open System for Earthquake Engineering Simulation），是一個(gè)能進(jìn)行結(jié)構(gòu)抗震彈塑性分析的開源軟件[10]。相比內(nèi)部代碼不透明的商用有限元軟件，OpenSees的開源性使得研究者能全面掌握它的運(yùn)行機(jī)制，并能根據(jù)自身需求進(jìn)行二次開發(fā)，因此，OpenSees是進(jìn)行建筑震害研究的重要工具。OpenSees中，混凝土梁柱也用纖維模型模擬。

    選用上述兩種軟件進(jìn)行震害分析，可以同時(shí)利用商業(yè)軟件的健壯性和開源軟件的透明性。

1.2  建筑震害分析的數(shù)據(jù)需求

    為了從IFC文件中提取充分的信息，以自動(dòng)生成Marc和OpenSees結(jié)構(gòu)模型，需要全面了解這兩種軟件的建筑震害分析建模方法，并總結(jié)震害分析所需的數(shù)據(jù)信息。這些數(shù)據(jù)信息分為兩類，一類與混凝土結(jié)構(gòu)模型有關(guān)，如結(jié)構(gòu)構(gòu)件的位置、幾何信息、材料信息、載荷等。另一類則與震害分析的方法和過程有關(guān)，如地震波的選取、載荷工況的增量步總數(shù)和每個(gè)增量步的時(shí)間步長(zhǎng)等，這類數(shù)據(jù)隨著分析方法和要求的變化而變化，并且在軟件中可以很方便地實(shí)現(xiàn)這類參數(shù)的調(diào)整。因此，應(yīng)該重點(diǎn)關(guān)注第一類數(shù)據(jù)需求，即與混凝土結(jié)構(gòu)模型有關(guān)的數(shù)據(jù)。

    通過考察這兩種軟件的建筑抗震分析模型，得到抗震分析的數(shù)據(jù)需求如表1所示。
 

表1  建筑抗震分析的數(shù)據(jù)需求
Tab.1  Data requirements for structural seismic analysis

2. 抗震分析的數(shù)據(jù)需求的IFC表達(dá)

圖5 載荷信息在IFC4中的描述方法
Fig.5  Description of load in IFC4 standard

3. 模型轉(zhuǎn)換程序設(shè)計(jì)

    根據(jù)第2章所建立的完整的IFC4結(jié)構(gòu)分析子信息模型框架，可以構(gòu)建滿足抗震分析數(shù)據(jù)需求的IFC4模型文件，進(jìn)而采用C++語言設(shè)計(jì)模型轉(zhuǎn)換程序，解析IFC模型文件，自動(dòng)生成Marc和OpenSees軟件模型。

    IFC模型由EXPRESS語言定義，是ASCII格式純文本文件形式，因此IFC文件需要經(jīng)過解析識(shí)別，才能供程序代碼處理。本文使用IFCEngine.dll作為IFC文件解析工具，這是一款免費(fèi)工具包，但要求編程者對(duì)IFC體系十分熟悉。

    根據(jù)表1所示建筑震害的數(shù)據(jù)需求，設(shè)計(jì)了節(jié)點(diǎn)、單元、支座約束、截面幾何、材料等5個(gè)類（class），并用Map、Vector等C++標(biāo)準(zhǔn)模板庫(kù)中的容器將各個(gè)類的不同對(duì)象有序地組織起來。程序的輸入模塊通過調(diào)用類的成員函數(shù)，從IFC模型文件中提取所需數(shù)據(jù)，賦值給相應(yīng)的數(shù)據(jù)成員。例如節(jié)點(diǎn)類Node調(diào)用SetCoordinate函數(shù)提取并設(shè)置節(jié)點(diǎn)坐標(biāo)。模型轉(zhuǎn)換程序的IFC解析及提取數(shù)據(jù)部分的流程如圖6所示。

圖6  模型轉(zhuǎn)換程序的IFC解析及提取數(shù)據(jù)部分的流程圖
Fig.6  Flow chart of the IFC parse and data obtaining part of the model converter

    由于采用了合理的數(shù)據(jù)結(jié)構(gòu)設(shè)計(jì)，程序的輸出模塊的實(shí)現(xiàn)比較方便，分別遍歷相應(yīng)容器，按照Marc和OpenSees的建模規(guī)則輸出節(jié)點(diǎn)、單元、材料等類的數(shù)據(jù)成員，便可自動(dòng)生成這兩種軟件模型。

4.  實(shí)例驗(yàn)證模型轉(zhuǎn)換效果

    為了驗(yàn)證模型轉(zhuǎn)換程序的正確性，用IFC4標(biāo)準(zhǔn)建立了一個(gè)鋼筋混凝土框架模型。模型為清華大學(xué)主樓的一榀10層3跨平面框架[16]，邊跨7.5米，中跨3.5米，框架總高度45.05米，模型詳細(xì)尺寸及配筋等信息見文獻(xiàn)。IFC模型以及轉(zhuǎn)換后的Marc模型如圖7所示（OpenSees沒有顯示模型的界面），由圖可知單元、節(jié)點(diǎn)、支座約束轉(zhuǎn)換正確。檢查模型參數(shù)，截面、材料、載荷等也都轉(zhuǎn)換正確。

圖7  IFC模型（左）及自動(dòng)生成的Marc模型（右）
Fig.7  IFC model(left) and auto-generated Marc model(right)

    分別用Marc和OpenSees對(duì)自動(dòng)生成的模型進(jìn)行模態(tài)分析，得到的前三階振型所對(duì)的頻率相差小于1%，如表3所示：

表3  Marc與OpenSees模態(tài)分析結(jié)果對(duì)比
Tab.3 Modal analysis result of Marc and OpenSees

    對(duì)框架輸入Y方向地震動(dòng)加速度，分別用Marc和OpenSees進(jìn)行彈塑性時(shí)程分析。選用El Centro地震波，地震峰值加速度設(shè)為400Gal，頂點(diǎn)位移的時(shí)程曲線如圖8所示，吻合得很好。

圖8  Marc與OpenSees彈塑性時(shí)程分析結(jié)果對(duì)比
Fig.8  Elasto-plastic time-history analysis result of Marc and OpenSees

 

    兩個(gè)軟件的模態(tài)分析和時(shí)程分析結(jié)果的一致性，進(jìn)一步說明了模型的單元、材料等各參數(shù)都得到了正確轉(zhuǎn)換。

5. 結(jié)  論

    本文從建筑震害分析的數(shù)據(jù)需求出發(fā)，逐一分析了這些數(shù)據(jù)需求在IFC4中的表達(dá)方法，建立了一個(gè)完整的IFC4結(jié)構(gòu)分析子信息模型。利用該子信息模型，就能通過IFC4標(biāo)準(zhǔn)來描述可進(jìn)行建筑抗震彈塑性分析的鋼筋混凝土結(jié)構(gòu)模型。

    根據(jù)建立的IFC4結(jié)構(gòu)分析子信息模型以及Marc和OpenSees軟件的建模規(guī)則，編寫了模型轉(zhuǎn)換程序，實(shí)現(xiàn)了節(jié)點(diǎn)、構(gòu)件、載荷、支座約束、材料、截面幾何等所有抗震分析所必要的信息的正確轉(zhuǎn)換，進(jìn)行了IFC文件到Marc和OpenSees軟件之間的數(shù)據(jù)轉(zhuǎn)換的新嘗試。

    在此基礎(chǔ)上，還有很多研究?jī)?nèi)容可以繼續(xù)開展。例如，擴(kuò)展本文建立的IFC4結(jié)構(gòu)分析子信息模型，使之能描述樓板、剪力墻等其他結(jié)構(gòu)構(gòu)件；考慮城市建筑群信息的IFC表達(dá)方法等。