何小涌，孫 濤，孟 妍，田鎮(zhèn)華

(陸軍勤務(wù)學(xué)院，重慶 401311)

0 前 言

門式剛架輕型房屋鋼結(jié)構(gòu)具有結(jié)構(gòu)體系受力簡(jiǎn)單、傳力路徑明確、構(gòu)件加工快捷、施工周期短、結(jié)構(gòu)布置靈活、經(jīng)濟(jì)效應(yīng)明顯等特點(diǎn)，被廣泛應(yīng)用于廠房、超市、展覽館、庫(kù)房等工業(yè)與民用建筑[1]。內(nèi)爆壓力由于受到結(jié)構(gòu)內(nèi)部空間的影響，沖擊波存在反射與匯聚效應(yīng)，因此與外爆沖擊波的特性存在著較大的差異。目前已有諸多學(xué)者[2-5]對(duì)箱室、艙室、工房、管廊等研究了內(nèi)爆的沖擊波特性、結(jié)構(gòu)失效模式及損傷變形，結(jié)果都表明分析對(duì)象的幾何形狀特征對(duì)沖擊波的傳播特性、結(jié)構(gòu)失效模式及損傷程度存在較大的影響。門式剛架輕型房屋鋼結(jié)構(gòu)相對(duì)于其他類型的結(jié)構(gòu)具有自己的幾何形狀特征，研究?jī)?nèi)爆沖擊波的傳播特性及分布規(guī)律對(duì)結(jié)構(gòu)抗爆設(shè)計(jì)與評(píng)估具有重要的現(xiàn)實(shí)意義。本文以某12 m跨門式剛架結(jié)構(gòu)為例分析特定爆炸位置下梁柱關(guān)鍵構(gòu)件的超壓分布特征。

1 數(shù)值模型

1.1 結(jié)構(gòu)簡(jiǎn)介

某12 m跨門式剛架結(jié)構(gòu)尺寸參數(shù)為跨度12 m，檐口高度5.4 m，柱距6 m。該結(jié)構(gòu)的屋面荷載設(shè)計(jì)值為0.75 kN/m2，基本風(fēng)壓為0.5 kN/m2，梁柱為高頻焊接H型鋼截面，尺寸為L(zhǎng)H200×150×4.5×6，其平面與山墻布置分別如圖1～2所示[6]。圖1中，當(dāng)結(jié)構(gòu)長(zhǎng)度不超過(guò)45 m時(shí)只需要在兩端開間設(shè)置橫向水平支撐，當(dāng)結(jié)構(gòu)長(zhǎng)度大于45 m時(shí)需要在中部開間增設(shè)橫向水平支撐。

圖1 平面布置示意圖(單位：mm)

圖2 山墻面布置示意圖(單位：mm)

1.2 爆源與超壓監(jiān)測(cè)點(diǎn)設(shè)置

圖3給出了爆源與超壓監(jiān)測(cè)點(diǎn)的設(shè)置，為模擬固定位置堆放的易爆物爆炸后門式輕鋼庫(kù)房梁柱關(guān)鍵構(gòu)件的超壓分布規(guī)律，在固定炸藥當(dāng)量為6.33 kg TNT的前提下，距左柱3 m，高度1.5 m處設(shè)置炸藥埋設(shè)點(diǎn)D1。采用LS-DYNA分析埋設(shè)點(diǎn)D1爆炸后在梁柱構(gòu)件邊緣產(chǎn)生的超壓分布規(guī)律，在梁柱關(guān)鍵構(gòu)件的邊緣各設(shè)置4個(gè)超壓監(jiān)測(cè)點(diǎn)以研究超壓時(shí)程曲線的變化規(guī)律。

圖3 爆源與監(jiān)測(cè)點(diǎn)分布

1.3 材料模型

分析模型包括炸藥、空氣兩種材料介質(zhì)，門式剛架的梁柱邊緣則考慮為剛性界面。空氣材料選用理想氣體模型，其模型參數(shù)與狀態(tài)方程參數(shù)分別見(jiàn)表1～2[7]。炸藥選擇TNT高爆炸藥模型，采用JWL狀態(tài)方程來(lái)描述，該方程給出了爆炸產(chǎn)物化學(xué)能所產(chǎn)生壓力的數(shù)學(xué)表達(dá)式：

(1)

式中，E0為單位體積炸藥的內(nèi)能；V為相對(duì)體積；A、B、R1、R2、ω為炸藥的特性參數(shù)，TNT炸藥的特性參數(shù)見(jiàn)表3[7]。

表1 *MAT_NULL模型參數(shù)

表2 *EOS_IDEAL_GAS狀態(tài)方程參數(shù)

表3 TNT炸藥特性參數(shù)

1.4 簡(jiǎn)化處理

采用數(shù)值分析爆炸沖擊波產(chǎn)生的超壓分布規(guī)律需要?jiǎng)澐州^細(xì)的網(wǎng)格，而門式剛架庫(kù)房結(jié)構(gòu)在空間上具有較大的長(zhǎng)度與寬度，如要模擬真實(shí)三維結(jié)構(gòu)的超壓分布則計(jì)算量會(huì)超過(guò)普通計(jì)算機(jī)甚至工作站的極限。本研究的實(shí)際模型是沿著門式剛架長(zhǎng)度堆滿易爆物，但是易爆物在寬度與高度方向的位置不同，考慮到長(zhǎng)度一般遠(yuǎn)大于寬度和高度，且主要研究梁柱關(guān)鍵構(gòu)件的超壓分布，因此可以簡(jiǎn)化為二維模型，如圖4所示。

圖4 模型的簡(jiǎn)化示意

2 結(jié)果與分析

2.1 左立柱超壓特征

圖5給出了左柱測(cè)點(diǎn)G1至G4的超壓時(shí)程曲線。位于柱腳處的監(jiān)測(cè)點(diǎn)G1的超壓時(shí)程曲線具有多次超壓峰值，但第一次超壓峰值遠(yuǎn)大于之后的超壓峰值，該時(shí)刻(0.837 ms)的超壓云圖如圖6所示，第一次超壓峰值是地面與立柱反射的沖擊波沿著柱腳方向傳播并相遇產(chǎn)生匯聚導(dǎo)致的結(jié)果。測(cè)點(diǎn)G2與G3在立柱的三分點(diǎn)處，且超壓時(shí)程曲線較為相似，具有多個(gè)超壓峰值。測(cè)點(diǎn)G4位于梁柱交接處，盡管該處到爆源D1的距離與G2、G3到爆源D1的距離相比更遠(yuǎn)，但是從曲線可以發(fā)現(xiàn)測(cè)點(diǎn)G4的超壓峰值比測(cè)點(diǎn)G2、G3的超壓峰值大很多。圖7給出了測(cè)點(diǎn)G4出現(xiàn)第一次超壓峰值時(shí)刻(1.292 ms)的超壓云圖，可以發(fā)現(xiàn)梁柱交接處的超壓峰值是由梁柱反射的沖擊波向節(jié)點(diǎn)處傳播并相遇產(chǎn)生匯聚導(dǎo)致的結(jié)果，與測(cè)點(diǎn)G1第一超壓峰值的形成機(jī)理相似。

此外，測(cè)點(diǎn)G2至G4的超壓時(shí)程曲線都有多個(gè)超壓峰值且逐漸衰減，但前2次或前3次的超壓峰值具有同樣的數(shù)量級(jí)，能否忽略第二次及以后的超壓峰值對(duì)結(jié)構(gòu)的影響還有待進(jìn)一步研究，尤其是如果第一次超壓峰值已經(jīng)導(dǎo)致構(gòu)件的彎曲或剪切損傷后，后繼超壓峰值能夠繼續(xù)對(duì)構(gòu)件造成損傷具有較大的概率。

圖5 測(cè)點(diǎn)G1至G4的超壓時(shí)程曲線

圖6 測(cè)點(diǎn)G1第一次超壓峰值時(shí)刻的超壓云圖

圖7 測(cè)點(diǎn)G4第一次超壓峰值時(shí)刻的超壓云圖

2.2 左屋梁超壓特征

圖8給出了左屋梁各測(cè)點(diǎn)的超壓時(shí)程曲線。測(cè)點(diǎn)G5至G7的超壓曲線都表現(xiàn)相似的規(guī)律，都具有多次超壓峰值且峰值逐漸衰減，但是前2次甚至是前3次的峰值差距不大，因此對(duì)構(gòu)件可能造成累積的損傷。從圖8中還可發(fā)現(xiàn)，雖然測(cè)點(diǎn)G4至G7到爆源的距離相差不大，但是G5至G7的超壓峰值與G4的超壓峰值相差較大，主要原因在于G5至G7處雖然有沖擊波的反射，但沒(méi)有G4的條件使得反射后的沖擊波再次進(jìn)行匯聚。

圖8 測(cè)點(diǎn)G4至G7的超壓時(shí)程曲線

2.3 右屋梁超壓特征

圖9給出了右側(cè)屋梁各測(cè)點(diǎn)的超壓時(shí)程曲線。測(cè)點(diǎn)G7至G10所表現(xiàn)出的規(guī)律與測(cè)點(diǎn)G4至G7所表現(xiàn)出的規(guī)律相似，各測(cè)點(diǎn)都具有多個(gè)超壓峰值且逐漸減小，但是前3次或前4次的超壓峰值在同一個(gè)數(shù)量級(jí)。如果第一次超壓峰值能夠?qū)?gòu)件造成較大的影響，那么后繼超壓峰值對(duì)構(gòu)件的影響將不能忽略，尤其是測(cè)點(diǎn)G10所表現(xiàn)出的前2次超壓幾乎相近的情況。此外，從圖9可以發(fā)現(xiàn)測(cè)點(diǎn)G10的前2次超壓峰值遠(yuǎn)大于測(cè)點(diǎn)G7至G9的超壓峰值，盡管測(cè)點(diǎn)G10離爆源最遠(yuǎn)，但經(jīng)過(guò)梁柱反射后的沖擊波向角部傳遞時(shí)還會(huì)產(chǎn)生匯聚,該時(shí)刻(2.542 ms)超壓云圖如圖10所示。

圖9 測(cè)點(diǎn)G7至G10的超壓時(shí)程曲線

圖10 測(cè)點(diǎn)G10第一次超壓峰值時(shí)刻的超壓云圖

2.4 右立柱超壓特征

圖11給出了右立柱各測(cè)點(diǎn)的超壓時(shí)程曲線。測(cè)點(diǎn)G10至G13的超壓曲線所表現(xiàn)出的規(guī)律與測(cè)點(diǎn)G7至G10的超壓規(guī)律較為相似，具有多個(gè)超壓峰值且總體趨勢(shì)上峰值逐漸減小。圖11中測(cè)點(diǎn)G12的曲線較為特殊，第2次超壓峰值居然略大于第1次超壓峰值，應(yīng)該值得后續(xù)的深入研究。圖12給出了測(cè)點(diǎn)G13第1次超壓峰值時(shí)刻(2.318 ms)的超壓云圖，可以看出由于地面與立柱形成了直角的幾何特征，沖擊波在此處反射與匯聚形成了第1次超壓峰值且為最大的超壓峰值。

圖11 測(cè)點(diǎn)G10至G13的超壓時(shí)程曲線

圖12 測(cè)點(diǎn)G13第一次超壓峰值時(shí)刻的超壓云圖

2.5 梁柱全長(zhǎng)超壓峰值比較

從測(cè)點(diǎn)G1至G13的超壓時(shí)程曲線可以發(fā)現(xiàn)梁柱不同部位的超壓峰值各不相同，主要原因有以下兩點(diǎn)：①各測(cè)點(diǎn)到爆源的距離不同；②不同位置會(huì)存在著沖擊波的發(fā)射與匯聚，幾何形狀將會(huì)對(duì)超壓峰值產(chǎn)生較大的影響。然而，以上研究規(guī)律若要應(yīng)用到工程設(shè)計(jì)中還有待簡(jiǎn)化，目前的工程防爆設(shè)計(jì)采用統(tǒng)一的超壓施加在構(gòu)件上使得設(shè)計(jì)過(guò)于簡(jiǎn)化且不能反映出結(jié)構(gòu)的形狀特征對(duì)超壓峰值的影響。

沿著構(gòu)件長(zhǎng)度繪制出最大超壓峰值可以反映結(jié)構(gòu)的大小與形狀特征的影響，對(duì)于工程設(shè)計(jì)的簡(jiǎn)化及改進(jìn)將具有參考價(jià)值。通過(guò)加密超壓監(jiān)測(cè)點(diǎn)并提取各點(diǎn)的最大超壓峰值，采用樣條曲線連接各加密監(jiān)測(cè)點(diǎn)的超壓峰值可以得到圖13所示的壁面超壓分布模型且可以發(fā)現(xiàn)以下規(guī)律。

1)在門剛結(jié)構(gòu)柱腳、柱頂?shù)冉屈c(diǎn)位置處，由于沖擊波的反射與匯聚，超壓峰值比構(gòu)件中部大很多，在設(shè)計(jì)或預(yù)防時(shí)值得注意。

2)沿構(gòu)件長(zhǎng)度存在小區(qū)段的超壓峰值急速增長(zhǎng)或衰減，大部分區(qū)段的超壓峰值相差并不大。為簡(jiǎn)化工程設(shè)計(jì)，可將構(gòu)件的超壓峰值沿長(zhǎng)度劃分區(qū)段，將超壓峰值變化取為線性變化或水平不變，這樣可以避免結(jié)構(gòu)內(nèi)力分析時(shí)各點(diǎn)超壓不同帶來(lái)的繁瑣的荷載輸入。

3)圖13僅為固定爆點(diǎn)D1爆炸后所得到的超壓峰值壁面分布模型，改變爆點(diǎn)位置將會(huì)得到不同的分布模型，壁面超壓分布模型更適合評(píng)估特定堆放的易爆物爆炸對(duì)結(jié)構(gòu)的影響。由此得到不同爆點(diǎn)位置對(duì)結(jié)構(gòu)的影響規(guī)律不僅有利于指導(dǎo)如何合理地堆放易爆物，還可在已有堆放方式的情況下采取合理措施防止爆炸產(chǎn)生較大的危害作用。

圖13 壁面超壓分布模型(P0=1.0GPa)

3 結(jié) 論

為模擬固定位置堆放的易爆物在意外爆炸后沖擊波對(duì)門式輕鋼庫(kù)房的影響，以某12 m跨門式剛架結(jié)構(gòu)為例，分析了爆點(diǎn)D1爆炸后沖擊波在梁柱關(guān)鍵構(gòu)件上的超壓分布規(guī)律，主要結(jié)論如下。

1)在沿庫(kù)房長(zhǎng)度規(guī)律堆放的易爆物的情況下，計(jì)算模型可簡(jiǎn)化為平面爆炸模型，能大幅縮小計(jì)算資源的需求。

2)由于門剛結(jié)構(gòu)特殊的幾何特征，沖擊波會(huì)在柱腳區(qū)段及柱頂區(qū)段等角部位置反射與匯聚，超壓峰值相比梁柱構(gòu)件中部區(qū)段有較大的提升，在工程設(shè)計(jì)或采取預(yù)防措施時(shí)值得關(guān)注。

3)如需進(jìn)一步分析結(jié)構(gòu)構(gòu)件在超壓下的響應(yīng)，可先繪制壁面超壓分布模型，依據(jù)分布規(guī)律將超壓分區(qū)段簡(jiǎn)化為線性增強(qiáng)或衰減段以及水平不變段。

[ID:013663]