考慮防火門易損性的震后火災模擬

王海東 袁佳欣

摘? ?要：探討一種基于BIM技術的震后次生火災模擬技術，定量評估防火門損傷對室內火災溫度場分布的影響. 通過選取3條場地類型及地震設計分組一致的地震波對9層Benchmark鋼結構模型進行非線性結構分析，并以層間位移角作為非結構構件門損傷的評估指標，結合FEMA P-58中門易損性曲線建立防火門損傷模型，進而基于BIM且利用得到的防火門損傷模型快速準確地建立建筑火源信息模型，導入Pyrosim模擬分析其對火災溫度分布的影響，并得到基于燃燒物與著火點位置關聯下的真實溫度場分布. 研究結果表明：防火門損傷數量的增加直接影響框架柱的溫度，且數量越多框架柱溫度越高;除位于著火點及未燃燒的防火分區外，位于其他防火分區的框架柱考慮門損傷的溫度明顯高于未考慮門損傷的情況，最大達483.98 ℃;且著火點位于不同層數時，各層達到溫度邊界值的框架柱數量存在差異.

關鍵詞：震后火災;建筑信息模型;易損性;建筑火源信息模型;溫度場分布

中圖分類號：TU391;TU243.2 ? ? ? ? ? ? ? ? ?文獻標志碼：A

Fire Following Earthquake Simulation Analysis

Considering of Fire Door Fragility

WANG Haidong1，2，YUAN Jiaxin2

（1. Key Laboratory of Building Safety and Energy Efficiency of Ministry of Education（Hunan University），Changsha 410082，China;

2. College of Civil Engineering，Hunan University，Changsha 410082，China）

Abstract：Post-earthquake fire simulation technology based on BIM was discussed to quantitatively evaluate the impact of fire door damage on the fire temperature distribution. Three seismic waves consistent with the site type and its seismic design group were selected to analyze the nonlinear structure of the 9-story Benchmark steel structure model，and the inter-story displacement angle? was used as the evaluation index of the non-structural component door damage，combined with the door vulnerability curve of FEMA P-58 to establish the door damage model. Then，using the fire door damage model obtained by BIM，the information model of the fire was quickly and accurately established. The Pyrosim was imported to simulate and analyze its impact on the fire temperature distribution，and the real temperature field distribution based on the correlation between the burning substance and the fire location was obtained. The results show that the increase in the number of door damages directly affects the temperature of the frame columns，and the larger the number，the higher the temperature of the frame columns; except for the fire and unburned fire compartments，the temperature of frame columns in other fire compartments with regard to door damage is significantly higher than that without considering door damage，with a maximum of 483.98 ℃. When the fire is located at different layers，the number of frame columns that reach the temperature boundary value of each layer is different.

Key words：post-earthquake fire;Building Information Modelling（BIM）;fragility;the information model of the fire;fire temperature distribution

地震過后往往易發生次生災害，難以控制且損失巨大. 例如：1923年日本關東大地震、1999年臺灣集集大地震、2008年汶川地震、2011年日本地震和2014年智利地震等，震后引發大火，火災造成建筑物損壞甚至倒塌，因而震后火災結構的性能已成為國內外學者研究的一個熱點. Jelinek等[1]對鋼框架在震后火災下的反應與僅受火災影響時結構的反應進行對比研究;Memari等[2]提出一種基于概率性能的地震火災分析框架，并利用該框架分析鋼結構在地震和火災連鎖危害下的反應，同時研究火災荷載密度等因素對結構失效概率的影響;Suwondo等[3]研究了不同的火災位置和地震破壞對荷載重分布路徑和結構倒塌的影響;張文[4]通過模糊數學方法建立震后建筑火災危險性評估模型，研究震后火災對建筑結構的影響;毛小勇等[5]對輕鋼-混凝土組合梁在標準升溫下進行了抗火性能的研究;李強等[6]對先后經歷地震和火災作用的破壞形態、損傷機理進行分析，研究損傷指數對圓鋼管混凝土柱震后耐火極限的影響. 然而，從已有的研究中發現，研究者對于火災發生的位置、溫度分布，尤其是與火災溫度密切相關的燃燒物數量、分布、防火分區等等都是根據自己的假定而并不與建筑本身的信息相關，缺乏一定的合理性. 與此同時，基于性能的抗震設計開始得到廣泛的重視和發展[7]，筆者發現部分學者將建筑信息模型（Building Information Modelling，BIM）技術引入到火災分析中，具有參考意義.

王婷等[8]將BIM模型導入火災模擬軟件PyroSim中，實現簡化火災建模;陸揚[9]將BIM 技術運用到性能化防火設計中，實現了消防性能相關模擬;杜二峰等[10]設計建造縮尺比為1 ∶ 4的門式剛架廠房模型，測量其模型內部各關鍵位置的熱空氣和鋼構件溫度. 但僅考慮消防性能并未將非結構構件考慮于結構有限元模型中，故由此評估建筑整體防火性能可能缺乏安全性. 且隨著人們對建筑使用性的要求不斷提高，利用基于FEMA P-58的抗震性能評估方法，開始專注于非結構構件的評估是抗震性能評估的趨勢. 而防火門這一非結構構件是建筑在發生火災時防止火災蔓延的重要措施，地震發生將導致門有不同程度的破壞，從而影響結構溫度場分布，最終影響結構的整體性能，因此通過有限元模型與該非結構構件進行聯結，并對建筑整體的防火性能評估應值得關注. Xu等[11]利用BIM和美國抗震設計標準FEMA P-58開展了建筑消防噴淋系統的震害評價，表明其建筑抗火性能與未考慮地震作用時有明顯區別. 但是，假定火災位置進行火災模擬，未將建筑內燃燒物與火災位置進行關聯，且未明確防火分區，難以得出震后火災溫度場的真實分布，可能缺乏合理性.

本文探討一種利用BIM技術和FEMA P-58的考慮非結構構件破壞的次生火災模擬方法，定量評估由于門損傷對火溫分布的影響. 該方法基于BIM快速準確地建立建筑火源信息模型，并通過FEMA P-58的抗震性能評估方法評估震后門損傷的情況，進而通過Pyrosim模擬分析防火門破損程度對火災溫度分布的影響，得到基于燃燒物與火災位置關聯下的真實溫度場分布，為后續震后火災結構整體性能研究提供基礎.

1 建模分析

1.1? ?建模思路

考慮震后防火門破損的火災溫度模擬主要包括2種不同的模型：①研究地震作用下非結構構件性能的有限元模型;②反映建筑物實際防火分區和燃燒物信息的火災模擬模型.

如何將2種獨立分割的模型信息統一，本文提出了基于BIM技術建立綜合的建筑信息模型的方法，將以上2種模型有效聯結，具體思路見圖1. 該思路核心為建筑BIM信息模型，即為實現火源信息化的關鍵，并利用FEMA P-58得出非結構構件的易損性分析，同時結合有限元模型中結構時程分析得出門損傷概率模型，并利用火災模擬軟件FDS模擬震后火災，分析得出該模型溫度的影響，更有利于將地震及震后火災結合考慮. 其中，有限元模型主要解決地震作用下非結構構件的破損分析問題;FDS火災模型主要解決火災模擬分析以及溫度實時分析問題;BIM建筑信息模型主要提供結構信息、建筑防火分區、防火門破損位置、狀態以及建筑火源分布情況，并貫穿于整個分析流程，為該建模思路的核心模型.

對于體量較大、非結構構件較多且需要建立建筑火源信息的模型，直接建立含有建筑火源信息的火災數值模型的工作量龐大，故將建筑內火源信息建立于BIM模型中，從而為整體火災模型的火源信息提供精細化的數據基礎，進而通過格式的轉換將其導入火災模擬軟件中得到建筑火源信息FDS模型. 將BIM模型轉換為FDS模型，就是利用BIM軟件將模型導出為fbx文件，再導入至圖形化操作軟件PyroSim中，軟件將正確識別文件格式，并按構件完成分組，此時對模型設置網格數量并對各構件材料屬性進行定義，最終運行.

該模型充分實現BIM技術“一模多用”的特點，減少了直接利用PyroSim建立FDS模型的復雜流程，而許鎮等[12]，采取直接建模及轉換生成FDS模型的方式分別進行火災模擬，模擬得出同一位置的煙氣濃度曲線圖，兩種方式得出的結果基本一致，證明該轉換方式具備較高的準確性.

1.2? ?模型建立

1.2.1? ?BIM模型建立

本文將以9層Benchmark模型為例，利用Autodesk Revit2017進行建筑建模，圖2為該建筑的三維BIM模型，模型體現建筑功能分區、建筑火源信息、結構全部信息，且各層采用相同的平面布局. 圖3為該建筑的平面布置圖. 根據美國防火規范NFPA 5000-2012將在該建筑每層設置6個防火分區，圖4為各層平面防火分區分布圖.

1.2.2? ?有限元模型創建

本文選用9層Benchmark模型，平面尺寸為45.73 m × 45.73 m，X方向、Y方向對稱，各5跨，每跨9.15 m，首層高度為5.49 m，標準層高度為3.96 m，梁、柱單元材料分別采用的是屈服強度為248 MPa和345 MPa的工字鋼，表1給出了9層Benchmark鋼框架結構模型的構件尺寸，詳細參數見文獻[13]. 利用SAP2000有限元軟件對9層Benchmark三維框架模型進行結構分析，阻尼比ζ = 0.05. 抗震設防烈度為8度（0.20g），場地類別為Ⅱ類（D類），場地分組為第1組.

梁和柱單元塑性鉸定義為SAP2000中基于FEMA365中的默認鉸屬性的集中塑性鉸模型. 其中柱單元采用的是考慮彎矩和軸力相互作用的P-M2-M3耦合塑性鉸，梁單元采用的是考慮彎矩的M3鉸，分別布置在梁柱0.1倍和0.9倍的長度位置處. 圖5為該9層Benchmark模型的三維框架立面及平面布置圖，圖6為利用SAP2000建立的有限元模型.

1.2.3? ?火災模擬模型創建

將該9層Benchmark建筑BIM模型導出為fbx文件，并將導出的文件導入至Pyrosim中，補充應有的信息并設置相應的參數可完成模型轉換，圖7為該9層Benchmark模型轉換完成的火災數值FDS模型.

2? ?地震后防火門破損評估分析

2.1? ?易損性數據

本文將利用基于FEMA P-58來考慮震后門損傷對火災的影響. 在美國規范FEMA P-58中提供了門這一非結構構件的易損性情況的數據，其中包含門的易損性曲線以及兩種不同破壞程度的說明. 表2為FEMA P-58中提供的門易損性數據表.

圖8為FEMA P-58中提供的門的易損性曲線，將其與地震損傷狀態相關聯，該曲線中包括兩種不同的破壞狀態分別為DS1和DS2，其中橫坐標為層間位移角數據，縱坐標為對應層間位移角下的損傷概率. DS1指門受到一定程度的變形，導致門的開合受到影響，而當結構變形消失時，門的變形也隨之消失，為彈性變形，據研究顯示，一般發生于層間位移角為0.001左右時;而DS2指門受到不可自行恢復的損傷且需要修理，據研究顯示，一般發生于層間位移角為0.004～0.01時.

2.2? ?門損傷模型

2.2.1? ?門損傷概率模型

本文選取9層Benchmark三維框架模型，以最大層間位移角θmax作為結構的損傷判定指標，利用王海東等[14]結合我國抗震規范及美國HAZUS[15]耐震規范中對不同損傷的定義劃分. 對應于FEMA P-58提供的門的易損性曲線，由于4種損傷狀態下結構均已處于結構非線性階段，故選用DS2為門損傷狀態計算其門損傷概率. 表3為結合門的易損性曲線及各破壞階段下對應的層間位移角限值，得出在不同損傷狀態下建筑內門的整體損傷概率.

2.2.2? ?門損傷評估

采用非線性評估理論，依照Memari等[16]在研究過程中選取的地震波，并參考我國《建筑抗震設計規范》的標準，選取表4中3條場地類型為Ⅱ類，參照之前學者研究成果[17]，依照震中距及震級得出各地震波的相應設計分組，以便工程應用.

將以上選取的3條地震波按規范要求，做最大值歸一化處理，按照抗震設防烈度為8度（0.20g）的罕遇地震情況，取PGA = 0.4g進行調幅. 作用于9層Benchmark三維框架模型中進行非線性結構分析，以各層的層間位移角這一結構評估值作為非結構構件門損傷的評估指數. 表5為最終得出的在各地震記錄影響下模型各層層間位移角的平均值及層占比，從而得出建筑模型下，在不同損傷狀態時，各層門的損傷情況.

2.2.3? ?9層Benchmark門損傷模型

該建筑火源信息模型中，門總數為289樘，其中每層防火門為11樘，其他為普通門，各層防火門位置一致，一層防火門位置如圖9所示，其中防火門1-1～1-7防火門尺寸均為M1521，1-8～1-11均為M1821. 依照前文所述，將選取的地震記錄作用于該建筑的有限元模型中，得到各層平均層間位移角并進行層占比計算. 通過門的易損性曲線對應各狀態下門的損傷比例，得出不同狀態下各層門損傷數量，從而得出該模型的門損傷模型，表6為地震記錄下，結構不同損傷狀態時各層門損傷數量.

其中，當結構處于嚴重破壞及倒塌狀態時，門的損傷概率已到達100%，則門損傷數量不依照各層層間位移角平均值的層占比分布，故不考慮. 而結構在中等破壞狀態下，考慮最不利布置原則，由于各層門損傷數量均大于11樘（各層防火門數量），則最不利情況為各層防火門均損傷，此時各層火災均蔓延至各防火分區，均已達最不利狀態，故本文不考慮.

本文后續將考慮在結構輕微破壞下建筑的火災模擬分析，并考慮最不利布置原則，將建筑模型中各層損傷的門均設置為防火門，由此更大程度地擴大室內火災區域，從而考慮各種火災發生的最不利情況.

3? ?火災模擬分析

本文以結構在輕微破壞狀態下為例，利用PyroSim火災模擬軟件，分別將有無門損傷兩種情況對火災的影響進行模擬實驗，模擬時間取為500 s. 由于同層門損傷狀態下發生火災，對其他層溫度的影響較小，故對各層分別進行模擬實驗. 本文選取模型的1層、5層及9層為模擬層數，得出最低、正中間層、最高層各柱處最高溫度曲線.

3.1? ?火災場景設置

圖10中所示的起火點區域位于①防火分區（見圖11），該區域在建筑設計中設置了茶水間，內有大功率用電設備，且與各防火分區均有接觸面，當防火門損傷時，火災具有最大的擴散面積，為最不利，且圖10表示各層考慮門未損傷時的火災場景. 圖11為各層柱編號所對應框架柱所在位置及防火分區. 同時，由于該建筑功能為辦公樓，考慮反應過程、燃燒產物等因素，設定與現實場景相符的起火場景及相應的化學反應. 本文則根據辦公室的起火特性，選取辦公室標準的Polyurethane GM37（聚氨酯）反應起火. 由于辦公樓可能存在大量的沙發坐墊以及可能存放的紙箱，依據最不利原則，設定火災為快速增長型，火災增長系數α取0.046 89，火源最大熱釋放速率為1 000.0 kW，按照t2火模型計算得到火災到達最大速率的時間為146 s，在本文中，由于設置的燃燒物有限，模擬火災過程影響，發現各柱在500 s內均已達到火災溫度下降段，故選取該模型模擬時間為500 s.

本文主要研究該模型結構在輕微破壞下，1、5、9層在有無考慮門損傷狀態兩種情況下對結構溫度的影響. 由前文可知，5層門損傷的數量為3，1、9層門損傷數量為4，由此進行模擬計算.

由于1、5、9層的門損傷數量以平面視角分類，則共有2種損傷類型，分別為門損傷數量為3和4. 考慮火災現場設計的最不利原則，應將損傷的防火門設置在火源物質影響更大的位置處，故選擇圖12、圖13所示的設置方式，門損傷數量分別為3和4.

3.2? ?溫度數據

在以上火災模擬中，共有6種工況，分別為工況1：9層4防火門損傷;工況2：9層門未損傷;工況3：5層3防火門損傷;工況4：5層門未損傷;工況5：1層4防火門損傷 ;工況6：1層門未損傷. 由于模型中各層框架柱共計36個，模擬500 s內各層柱中部溫度變化情況，并記錄其最高溫度，按各防火分區進行統計，分別得出圖14～圖19各防火分區6種不同工況的最高溫度情況，表7為各防火分區內包含的柱編號.

3.3? ?模擬結果分析

由各層柱在6種工況、不同門損傷情況下統計得到的最高溫度可知，門損傷數量的增多將對框架柱溫度有直接影響，且門損傷數量越多，框架柱溫度越高. 統計可得，在火災模擬500 s內，當框架柱最高溫度高于300 ℃時，則100 s內基本達到構件的最高溫度;當框架柱最高溫度低于300 ℃時，則在模擬時間為300 s后基本達到構件的最高溫度.

通過統計各工況下各框架柱的最高溫度可知：框架柱位于著火點及在所有工況下均未燃燒的防火分區時，考慮地震作用和未考慮地震作用的框架柱燃燒溫度基本吻合;而框架柱位于其他燃燒點的防火分區時，未考慮地震作用下門損傷情況的燃燒溫度將明顯低于考慮地震作用下門損傷情況的燃燒溫度. 且以本文1、5、9層對應考慮地震作用與未考慮地震作用的所有工況下，統計可得考慮門損與未考慮門損的框架柱溫度差值最大達483.98 ℃，該點為9層23號柱.

根據《建筑鋼結構防火技術規范》[18]可知，當溫度達到300 ℃以上時鋼材的屈服強度將進行折減，故本文以計算屈服強度折減系數為1.0、0.95、0.9時鋼材的溫度作為統計邊界值，分別為300 ℃、375 ℃、409 ℃. 圖20為統計各工況下框架柱最高溫度達到該三者統計邊界值以上的數量，為后續分析溫度對柱承載能力的影響提供基礎. 且由各框架柱最高溫度的統計結果可知，考慮地震作用門損傷時達到邊界溫度值的框架柱數量比不考慮地震作用門損傷時多達4.25～7.33倍;當考慮門損傷情況下，框架柱溫度及達到統計邊界值的數量排序均為：9層>1層>5層;當不考慮門損傷的情況下各層達統計邊界值數量趨于穩定.

4? ?結? ?論

本文利用BIM技術和FEMA P-58提出了一種基于性能的考慮非結構構件破損的震后次生火災模擬方法，定量評估由于門損傷對火災溫度場的影響，并反映于框架柱的溫度，相關結論如下.

1）提出了基于BIM的建筑火源信息模型，并將其作為信息傳遞中心，將結構模型與火災模型相結合.

2）將地震作用下的層間位移角與非結構構件門的易損性進行關聯：

①當結構處于嚴重破壞或倒塌狀態時，各層門損傷概率達100%，此時震后火災溫度場與防火門無關;

②當結構處于中等破壞時，門損傷概率為56%，此時各門為有概率的損傷，若考慮最不利情況為防火門全數損傷，火災蔓延至全部防火分區，將達到火災最不利影響;

③在結構處于輕微破壞時，門損傷的概率為10%，防火門是否損傷具有隨機性，考慮地震作用下防火門損傷情況對震后火災溫度場分布具有顯著的影響，故在中低地震設防區具有重要意義.

3）當結構處于輕微破壞時，通過對比6種工況下各框架柱的溫度變化情況及最高溫度統計可得：

①當框架柱位于非著火點及引燃的防火分區時，未考慮地震作用門損傷的燃燒溫度將明顯低于考慮地震作用門損傷的燃燒溫度，且框架柱溫度差值最大為483.98 ℃;

②通過統計鋼材屈服強度折減系數為1.0、0.95、0.9所對應的邊界溫度值300 ℃、375 ℃、409 ℃的框架柱數量可知，考慮門損時達到邊界溫度值的框架柱數量將遠遠高于未考慮地震作用門損時的框架柱數量，其比值最高可達7.33倍;

③達統計邊界值的數量具有差異但有一定規律性，在考慮地震作用門損傷時，各層達統計邊界值數量為9層>1層>5層.未考慮地震作用門損傷時，各層達統計邊界值數量趨于穩定.

參考文獻

[1]? ? JELINEK T，ZANIA V，GIULIANI L. Post-earthquake fire resistance of steel buildings[J]. Journal of Constructional Steel Research，2017，138（8）：774—782.

[2]? ? MEMARI M，MAHMOUD H. Framework for a performance-based analysis of fires following earthquakes[J]. Engineering Structures，2018，171（5）：794—805.

[3]? ? SUWONDO R，CUNNINGHAM L，GILLIE M. Progressive collapse analysis of composite steel frames subject to fire following earthquake[J]. Fire Safety Journal，2019，103（1）：49—58.

[4]? ? 張文. 基于模糊數學的震后建筑火災危險性評估模型研究[J].地震工程學報，2018，40（6）：1372—1377.

ZHANG W. A model of fire risk assessment of post-earthquake buildings based on fuzzy mathematics[J]. China Earthquake Engineering Journal，2018，40（6）：1372—1377. （In Chinese）

[5]? ? 毛小勇，肖巖. 標準升溫下輕鋼-混凝土組合梁的抗火性能研究[J]. 湖南大學學報（自然科學版），2005，32（2）：64—70.

MAO X Y，XIAO Y. Behavior of lightweight steel-concrete composite beams subjected to standard fire[J]. Journal of Hunan University（Natural Sciences），2005，32（2）：64—70. （In Chinese）

[6]? ? 李強，王文達，毛文婧. 震損后圓鋼管混凝土柱耐火性能有限元分析[J]. 建筑結構學報，2017，38（S1）：118—125.

LI Q，WANG W D，MAO W J. Finite element analysis of concrete-filled circular steel tubular columns subjected to post-earthquake fire[J]. Journal of Building Structures，2017，38（S1）：118—125. （In Chinese）

[7]? ? 陸新征，曾翔，許鎮. 建設地震韌性城市所面臨的挑戰[J]. 城市與減災，2017（4）：29—34.

LU X Z，ZENG X，XU Z. The challenge of building a seismically resilient city[J]. City and Disaster Reduction，2017（4）：29—34. （In Chinese）

[8]? ? 王婷，杜慕皓，唐永福，等. 基于BIM的火災模擬與安全疏散分析[J]. 土木建筑工程信息技術，2014，6（6）：102—108.

WANG T，DU M H，TANG Y F，et al. Analysis on fire model and safety evacuation based on BIM[J]. Journal of Information Technology in Civil Engineering and Architecture，2014，6（6）：102—108. （In Chinese）

[9]? ? 陸揚. 基于BIM的性能化分析手段在建筑防火設計中的研究與實踐[J]. 土木建筑工程信息技術，2011，3（4）：63—71.

LU Y. BIM-based performance analysis tools in building fire safety design of the research and practice[J]. Journal of Information Technology in Civil Engineering and Architecture，2011，3（4）：63—71. （In Chinese）

[10]? 杜二峰，舒贛平，呂曉. 1 ∶ 4縮尺廠房中功率火災試驗研究[J]. 湖南大學學報（自然科學版），2020，47（3）：82—88.

DU E F，SHU G P，L？譈 X.Medium power fire test study on a 1 ∶ 4 scaled portal frame building[J]. Journal of Hunan University （Natural Sciences），2020，47（3）：82—88. （In Chinese）

[11]? XU Z，ZHANG Z C，LU X Z，et al. Post-earthquake fire simulation considering overall seismic damage of sprinkler systems based on BIM and FEMA P-58[J]. Automation in Construction，2018，90（2）：9—22.

[12]? 許鎮，張宗才，陸新征，等. 基于BIM和FEMAP-58的建筑地震次生火災模擬[J]. 消防科學與技術，2018，37（6）：745—749.

XU Z，ZHANG Z C，LU X Z，et al. Building post-earthquake fire simulation based on BIM and FEMA P-58[J]. Fire Science and Technology，2018，37（6）：745—749. （In Chinese）

[13]? OHTORI Y，CHRISTENSON R E，SPENCER B F J，et al. Benchmark control problems for seismically excited nonlinear buildings[J]. Journal of Engineering Mechanics，2004，130（4）：366—385.

[14]? 王海東，武戰戰，劉巴黎. 地震動輸入對中長周期結構抗倒塌性能評估的影響[J]. 地震工程與工程振動，2019，39（3）：176—184.

WANG H D，WU Z Z，LIU B L. Influence of ground motion input on evaluation of anti-collapse performance of medium-long period structure[J].Earthquake Engineering and Engineering Dynamics，2019，39（3）：176—184. （In Chinese）

[15]? Users manual：HAZUS99[S]. Washington D C：Federal Emergency Management Agency，l999：236—268.

[16]? MEMARI M，MAHMOUD H，ELLINGWOOD B. Post-earthquake fire performance of moment resisting frames with reduced beam section connections[J]. Journal of Constructional Steel Research，2014，103（9）：215—229.

[17]? 冀昆，溫瑞智，任葉飛. 中國抗震規范強震記錄選取的初選條件研究[J]. 防災減災工程學報，2016，36（1）：44—49.

JI K，WEN R Z，REN Y F. Study on the initial of selection criteria strong motion records for Chinese seismic code[J]. Journal of Disaster Prevention and Mitigation Engineering，2016，36（1）：44—49. （In Chinese）

[18]? 建筑鋼結構防火技術規范：GB 51249—2017[S]. 北京：中國計劃出版社，2017：19—20.

Code for fire safety of steel structures in buildings：GB 51249—2017[S]. Beijing：China Planning Press，2017：19—20. （In Chinese）