引 言
隨著城市化的發(fā)展,地震對經(jīng)濟發(fā)展和人民生活的潛在威脅越來越嚴重,城市建筑震害預測成為防災減災的一個重要方面。城市區(qū)域內(nèi)的建筑群震害預測包含兩個層次:宏觀層次和微觀層次。宏觀層次模型將建筑抽象成一個或少量自由度的結(jié)構(gòu),進行不同烈度下地震損失分析預測。微觀層次模型[2]將建立建筑物更加精細化的分析模型(例如桿系模型等),可以考慮結(jié)構(gòu)的足夠多的自由度和精確的動力特性,進行高精度的震害分析預測。
城市建筑震害預測,需要獲取和存儲城市建筑的基礎數(shù)據(jù)。如果只是建立宏觀層次模型,則只需要建筑面積、年代、層數(shù)、結(jié)構(gòu)類型等少量宏觀數(shù)據(jù),其數(shù)據(jù)來源可以是城市GIS數(shù)據(jù)庫等?;贕IS的城市震害預測已經(jīng)有一定的研究基礎[3]。而建立微觀層次的精細化模型需要更詳細的數(shù)據(jù),如梁柱的尺寸、配筋等,若采用手動方法建立精細化模型,不僅工作量太大,也容易出錯,特別對于某些開源的有限元軟件(如OpenSees),沒有前處理界面,建模極不方便。比較合理的方法是利用建筑已有的數(shù)字化模型,從中提取建筑震害分析所需要的詳細數(shù)據(jù),自動建模供計算模塊進行分析計算。然而不同軟件生成的建筑數(shù)字化模型有著不同的數(shù)據(jù)格式,而計算模塊也可能由不同的結(jié)構(gòu)分析軟件組成(如Marc,OpenSees等)。因此,需要采用一個統(tǒng)一的數(shù)據(jù)標準,以實現(xiàn)不同軟件之間的數(shù)據(jù)交換。
建筑信息模型(Building Information Model,BIM)致力于在建筑的全生命期實現(xiàn)不同建筑工程領域的信息共享和協(xié)同管理,基于BIM的建筑數(shù)字化模型是未來重要的發(fā)展方向。目前,BIM采用IFC標準作為建筑產(chǎn)品的數(shù)據(jù)交換標準,不同軟件之間通過IFC這個統(tǒng)一的標準來交換數(shù)據(jù)。然而研究表明[4],基于IFC標準的建筑模型與結(jié)構(gòu)模型之間的信息傳遞與共享還不是十分成熟。張劍濤等[5],劉照球[6, 7]等開發(fā)了IFC與PKPM結(jié)構(gòu)模型的數(shù)據(jù)轉(zhuǎn)換接口,能初步實現(xiàn)二者的幾何、材料等信息的轉(zhuǎn)換,但缺乏載荷信息等的轉(zhuǎn)換。也有學者嘗試建立通用建筑結(jié)構(gòu)模型轉(zhuǎn)換平臺,以IFC為通用數(shù)據(jù)轉(zhuǎn)換標準,實現(xiàn)任意常用結(jié)構(gòu)分析軟件間的結(jié)構(gòu)分析模型數(shù)據(jù)轉(zhuǎn)換。但這類研究尚處于起步階段,目前僅實現(xiàn)了SATWE、ETABS、SAP2000等軟件與IFC標準的轉(zhuǎn)換接口,所能實現(xiàn)的數(shù)據(jù)轉(zhuǎn)換也僅限于幾何、材料等信息。
因此,能否利用BIM的共享性,解決城市建筑震害預測的數(shù)據(jù)獲取以及自動建模問題,具有重要的研究價值。
1. 建筑震害分析及數(shù)據(jù)需求
對于鋼筋混凝土框架結(jié)構(gòu),通常采用纖維梁模型作為抗震彈塑性分析的精細化計算模型。不少軟件都能具備建筑抗震彈塑性分析能力,其中兩類典型的軟件是Marc和OpenSees。
由于城市建筑信息是以單體建筑為基礎的,因此,本文利用C++編寫了一個模型轉(zhuǎn)換程序,從一個單體鋼筋混凝土建筑的(IFC格式的)數(shù)字化模型中,提取建筑震害分析所需的信息,自動生成Marc和OpenSees軟件模型,驗證轉(zhuǎn)換后的模型能否用于分析建筑結(jié)構(gòu)在地震作用下的響應,如圖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是一款功能齊全的高級非線性通用有限元軟件[8],用Marc進行建筑震害模擬,可以得到較為理想的分析結(jié)果,且適用于大多數(shù)分析需求。
基于纖維模型原理,清華大學土木工程系利用Marc軟件提供的二次開發(fā)接口,編寫了THUFIBER程序[9],使Marc軟件能對鋼筋混凝土結(jié)構(gòu)進行纖維模型計算。THUFIBER將混凝土構(gòu)件的截面分為36個混凝土纖維,以及4個鋼筋纖維(纖維劃分方法也可以由用戶自定義),各纖維的面積由截面尺寸和配筋決定。對每根纖維定義特定的本構(gòu)關系,便可以進行不同精度的有限元分析。
OpenSees(Open System for Earthquake Engineering Simulation),是一個能進行結(jié)構(gòu)抗震彈塑性分析的開源軟件[10]。相比內(nèi)部代碼不透明的商用有限元軟件,OpenSees的開源性使得研究者能全面掌握它的運行機制,并能根據(jù)自身需求進行二次開發(fā),因此,OpenSees是進行建筑震害研究的重要工具。OpenSees中,混凝土梁柱也用纖維模型模擬。
選用上述兩種軟件進行震害分析,可以同時利用商業(yè)軟件的健壯性和開源軟件的透明性。
1.2 建筑震害分析的數(shù)據(jù)需求
為了從IFC文件中提取充分的信息,以自動生成Marc和OpenSees結(jié)構(gòu)模型,需要全面了解這兩種軟件的建筑震害分析建模方法,并總結(jié)震害分析所需的數(shù)據(jù)信息。這些數(shù)據(jù)信息分為兩類,一類與混凝土結(jié)構(gòu)模型有關,如結(jié)構(gòu)構(gòu)件的位置、幾何信息、材料信息、載荷等。另一類則與震害分析的方法和過程有關,如地震波的選取、載荷工況的增量步總數(shù)和每個增量步的時間步長等,這類數(shù)據(jù)隨著分析方法和要求的變化而變化,并且在軟件中可以很方便地實現(xiàn)這類參數(shù)的調(diào)整。因此,應該重點關注第一類數(shù)據(jù)需求,即與混凝土結(jié)構(gòu)模型有關的數(shù)據(jù)。
通過考察這兩種軟件的建筑抗震分析模型,得到抗震分析的數(shù)據(jù)需求如表1所示。
表1 建筑抗震分析的數(shù)據(jù)需求
Tab.1 Data requirements for structural seismic analysis
2. 抗震分析的數(shù)據(jù)需求的IFC表達
IFC的復雜性,也是限制它發(fā)展、導致常用軟件普遍無法導出滿足抗震分析的IFC文件的原因之一。因此本文在仔細考察IFC4體系的基礎之上,選擇其中的常用部分,建立IFC4結(jié)構(gòu)分析子信息模型。例如配筋的描述僅保留了方法(4)。雖然這降低了IFC的全面性,但仍能描述抗震分析所需的關鍵信息,增強了實用性。
描述IFC(子信息)模型的最方便的工具是EXPRESS-G[15],圖3是它的部分圖例。
2.3.2 其他建筑抗震分析數(shù)據(jù)的IFC4描述方法
與2.3.1小節(jié)類似地,建立了結(jié)構(gòu)構(gòu)件、載荷、支座約束、材料、截面幾何特性等信息的IFC4描述方法。以載荷為例,如圖5所示,圖中實體屬性未全部列出。IFC4能描述線性變化的分布載荷,但為了降低復雜性,本文建議將這類載荷轉(zhuǎn)換成等效均布載荷,用IfcStructuralLinearAction實體表達。
圖5 載荷信息在IFC4中的描述方法
Fig.5 Description of load in IFC4 standard
3. 模型轉(zhuǎn)換程序設計
根據(jù)第2章所建立的完整的IFC4結(jié)構(gòu)分析子信息模型框架,可以構(gòu)建滿足抗震分析數(shù)據(jù)需求的IFC4模型文件,進而采用C++語言設計模型轉(zhuǎn)換程序,解析IFC模型文件,自動生成Marc和OpenSees軟件模型。
IFC模型由EXPRESS語言定義,是ASCII格式純文本文件形式,因此IFC文件需要經(jīng)過解析識別,才能供程序代碼處理。本文使用IFCEngine.dll作為IFC文件解析工具,這是一款免費工具包,但要求編程者對IFC體系十分熟悉。
根據(jù)表1所示建筑震害的數(shù)據(jù)需求,設計了節(jié)點、單元、支座約束、截面幾何、材料等5個類(class),并用Map、Vector等C++標準模板庫中的容器將各個類的不同對象有序地組織起來。程序的輸入模塊通過調(diào)用類的成員函數(shù),從IFC模型文件中提取所需數(shù)據(jù),賦值給相應的數(shù)據(jù)成員。例如節(jié)點類Node調(diào)用SetCoordinate函數(shù)提取并設置節(jié)點坐標。模型轉(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)設計,程序的輸出模塊的實現(xiàn)比較方便,分別遍歷相應容器,按照Marc和OpenSees的建模規(guī)則輸出節(jié)點、單元、材料等類的數(shù)據(jù)成員,便可自動生成這兩種軟件模型。
4. 實例驗證模型轉(zhuǎn)換效果
為了驗證模型轉(zhuǎn)換程序的正確性,用IFC4標準建立了一個鋼筋混凝土框架模型。模型為清華大學主樓的一榀10層3跨平面框架[16],邊跨7.5米,中跨3.5米,框架總高度45.05米,模型詳細尺寸及配筋等信息見文獻。IFC模型以及轉(zhuǎn)換后的Marc模型如圖7所示(OpenSees沒有顯示模型的界面),由圖可知單元、節(jié)點、支座約束轉(zhuǎn)換正確。檢查模型參數(shù),截面、材料、載荷等也都轉(zhuǎn)換正確。
圖7 IFC模型(左)及自動生成的Marc模型(右)
Fig.7 IFC model(left) and auto-generated Marc model(right)
分別用Marc和OpenSees對自動生成的模型進行模態(tài)分析,得到的前三階振型所對的頻率相差小于1%,如表3所示:
表3 Marc與OpenSees模態(tài)分析結(jié)果對比
Tab.3 Modal analysis result of Marc and OpenSees
對框架輸入Y方向地震動加速度,分別用Marc和OpenSees進行彈塑性時程分析。選用El Centro地震波,地震峰值加速度設為400Gal,頂點位移的時程曲線如圖8所示,吻合得很好。
圖8 Marc與OpenSees彈塑性時程分析結(jié)果對比
Fig.8 Elasto-plastic time-history analysis result of Marc and OpenSees
兩個軟件的模態(tài)分析和時程分析結(jié)果的一致性,進一步說明了模型的單元、材料等各參數(shù)都得到了正確轉(zhuǎn)換。
5. 結(jié) 論
本文從建筑震害分析的數(shù)據(jù)需求出發(fā),逐一分析了這些數(shù)據(jù)需求在IFC4中的表達方法,建立了一個完整的IFC4結(jié)構(gòu)分析子信息模型。利用該子信息模型,就能通過IFC4標準來描述可進行建筑抗震彈塑性分析的鋼筋混凝土結(jié)構(gòu)模型。
根據(jù)建立的IFC4結(jié)構(gòu)分析子信息模型以及Marc和OpenSees軟件的建模規(guī)則,編寫了模型轉(zhuǎn)換程序,實現(xiàn)了節(jié)點、構(gòu)件、載荷、支座約束、材料、截面幾何等所有抗震分析所必要的信息的正確轉(zhuǎn)換,進行了IFC文件到Marc和OpenSees軟件之間的數(shù)據(jù)轉(zhuǎn)換的新嘗試。
在此基礎上,還有很多研究內(nèi)容可以繼續(xù)開展。例如,擴展本文建立的IFC4結(jié)構(gòu)分析子信息模型,使之能描述樓板、剪力墻等其他結(jié)構(gòu)構(gòu)件;考慮城市建筑群信息的IFC表達方法等。