崔 瑩，趙均海，張常光，孫珊珊，陳 兵

（1.西安石油大學機械工程學院，西安 710065；2.長安大學建筑工程學院，西安 710061；3.西北工業(yè)大學航空學院，西安 710072）

復式空心鋼管混凝土柱抗爆性能及損傷研究

崔 瑩1，3，趙均海2，張常光2，孫珊珊2，陳 兵1

（1.西安石油大學機械工程學院，西安 710065；2.長安大學建筑工程學院，西安 710061；3.西北工業(yè)大學航空學院，西安 710072）

通過爆炸試驗，對復式空心鋼管混凝土柱抗爆性能及損傷程度進行研究，并結合所得試驗數據進行數值模擬，分析并建立柱的損傷評估準則。結果表明：在折合距離為0.14 m／kg1／3試驗條件的爆炸荷載作用下，復式空心鋼管混凝土柱迎爆面柱中發(fā)生了明顯的塑性彎曲變形；柱迎爆面壓力峰值柱中最高、柱底次之、柱頂最小，且爆炸沖擊波對迎爆面柱中造成的破壞最為強烈、柱底次之、柱頂最?。恢瞬抗?jié)點強度與構件整體強度協調對提升復式空心鋼管混凝土柱的抗爆性能尤為重要。最終建立了基于固端約束下復式空心鋼管混凝土柱中撓度的超壓－沖量（P－I）損傷準則及判定公式。

爆炸力學；復式空心鋼管混凝土柱；爆炸試驗；數值模擬；超壓－沖量（P－I）曲線

由于全球范圍內的爆炸事故和恐怖襲擊頻繁發(fā)生，常規(guī)建筑結構受到爆炸沖擊作用的可能性不斷增加。要使結構在爆炸沖擊荷載作用下避免破壞，柱本身的抗爆性能起著非常重要的作用。復式空心鋼管混凝土柱具有承載力高、塑性好、抗彎剛度大的特點，被廣泛用做工程結構中的重要受力構件［1－4］。研究復式空心鋼管混凝土柱的爆炸沖擊破壞效應，進而建立其損傷評估方法具有重要的理論意義和工程應用價值。近年來，國內外學者對于爆炸沖擊荷載作用下柱的動力響應開展了部分研究工作，主要采用數值模擬和試驗分析的方法研究了流固相互作用下的有缺陷簡支柱的動態(tài)彎曲特性［5－6］；柱端約束條件對框架柱在爆炸沖擊波作用下的影響及框架柱抗爆分析簡化計算［7］；鋼筋混凝土柱和實心鋼管混凝土柱在爆炸沖擊荷載下的動力響應、抗爆性能、破壞形式及影響要素等［8－13］。但是對于本文所提及的復式空心鋼管混凝土柱在爆炸荷載作用下動態(tài)響應的研究工作相對較少。本文以所開展的復式空心鋼管混凝土柱爆炸試驗結果為依據，對爆炸荷載下復式空心鋼管混凝土柱的動態(tài)響應進行數值模擬，建立了固端約束下基于構件最大撓度的復式空心鋼管混凝土柱損傷評估準則。

1 復式空心鋼管混凝土柱爆炸試驗

1.1 試件設計

復式空心鋼管混凝土柱試件以文獻［14］中所應用的鋼管混凝土柱為原型進行設計。按照縮放比例1∶2及《無縫鋼管尺寸、外形、重量及允許偏差》（GB／T17395－2008）中鋼管規(guī)格的要求，外鋼管采用國標外徑為273 mm，壁厚為7 mm的無縫鋼管，芯鋼管選擇為外徑50 mm，壁厚為3 mm的普通直縫焊管。內填混凝土選擇為C40的細石混凝土。柱腳參照《鋼管混凝土結構設計與施工規(guī)程》（CECS28－2012）以及《多、高層民用建筑鋼結構節(jié)點構造詳圖》（GJBT－543－01SG519）中端承式柱腳的構造進行設計，肋板與外鋼管采用焊接連接，焊腳尺寸不小于10 mm。具體試件參數見表1所示。

1.2 試驗準備

試件柱腳底板通過10.9級的M30高強螺栓與反力架完成固結。

試驗所用的TNT藥量為50 kg，爆心距為0.5m，爆心保持與柱迎爆面柱中等高，折合距離為0.14 m／kg1／3。現場示意圖如圖1所示。

表1 復式空心鋼管混凝土柱試件設計參數表Tab.1 Design parameters of dup lex hollow CFST column

圖1 試驗現場示意圖Fig.1 Schematic view of test site

1.3 試驗結果分析

爆炸荷載作用下反力架整體完好，未發(fā)生斷裂破壞，柱宏觀變形如圖2所示。

由圖2（a）可以看出，在折合距離為0.14 m／kg1／3的爆炸荷載作用下，復式空心鋼管混凝土柱迎爆面柱中發(fā)生了明顯的塑性彎曲變形，經量測柱中撓度為195 mm。同時柱腳底板迎爆面右下角有輕微翹曲變形，這是因為右下角與地面固結的螺栓較其余螺栓連接緊密，螺栓提供的約束較好，導致柱向后移動趨勢受阻，因而發(fā)生底板翹曲現象。

需要引起重視的是柱迎爆面柱底肋板上緣柱表面有一條細微未貫通裂縫存在，如圖2（b）所示。分析出現此條細微裂縫的原因是由于柱中發(fā)生彎曲凹陷，而柱底節(jié)點因為有肋板的設計而強度高于構件，柱底迎爆面受到爆炸沖擊波的剪切及拉伸作用所致。因此要使復式空心鋼管混凝土柱的抗爆性能得到提升，柱端部節(jié)點強度如何與構件整體強度協調尤為重要。

圖2 復式空心鋼管混凝土柱宏觀變形Fig.2 Macroscopic deformation of duplex hollow CFST column

2 復式空心鋼管混凝土柱數值模擬

2.1 材料本構模型

數值模擬對象包括鋼管混凝土柱、空氣及炸藥，均采用SOLID164單元模擬。單位制取為mm－ms－MPa。

（1）混凝土本構模型

選用LS-DYNA材料庫中的111號*MAT＿JOHNSON＿HOLMGUIST＿CONCRETE材料模型［15］（簡稱HJ－C模型）來模擬混凝土材料在高應變率下的動力特性。

（2）鋼材本構模型

選用LS-DYNA材料庫中的3號*MAT＿PLASTIC＿KINEMATIC材料模型［16］來模擬鋼管在高應變率下的動力特性。

（3）炸藥和空氣本構模型

采用LS-DYNA提供的*MAT＿HIGH＿EXPLOSIVE＿BURN高能炸藥燃燒模型來模擬TNT炸藥的爆轟。選擇LS-DYNA提供的*MAT＿NULL材料模型對空氣進行模擬。

2.2 材料狀態(tài)方程

在LS-DYNA中，對炸藥和空氣等流體材料的處理，需要同時使用材料的本構模型和狀態(tài)方程兩種方式來描述材料的特性：采用本構模型來描述應力和應變的關系，即Δσij和Δεij的關系，采用狀態(tài)方程來描述炸藥和空氣在爆炸發(fā)生后體積變形率與壓力的關系，即ΔV／V和ΔP的關系。

（1）炸藥的狀態(tài)方程

采用*EOS＿JWL狀態(tài)方程可以對爆炸過程中產生的壓力做出與試驗結果較為相近的預測。*EOS＿JWL狀態(tài)方程的表達式為：

式中，P為爆轟產物的壓力；V為爆轟產物的相對體積；E為爆轟產物單位體積初始內能；w、A、B、R1、R2為材料常數。

（2）空氣的狀態(tài)方程

采用線性多項式*EOS＿LINEAR＿POLYNOMIAL狀態(tài)方程來模擬空氣［16－18］。線性多項式*EOS＿LINEAR＿POLYNOMIAL狀態(tài)方程為：

2.3 有限元模型

（1）幾何模型及邊界條件

為了真實模擬柱受爆炸沖擊荷載后的動態(tài)響應，邊界條件確定如下：鋼管混凝土柱上下端為固端約束。空氣區(qū)域按照文獻［17］的建議及試驗的具體要求確定如下：柱背爆面及兩個側面空氣邊界按照各超過柱邊緣300 mm確定，上部邊界按照柱頂預留200 mm確定，下部邊界按照與地面相接處予以確定。同時模型在底部加設剛性反射面以模擬真實試驗時沖擊波在地面的反射，空氣邊界設定為透射邊界［18］。建立的整體模型幾何關系示意圖如圖3所示。

（2）有限元模型

空氣和炸藥考慮結點共用以及網格劃分的尺寸效應［19－20］，單元大小確定為20 mm。共劃分相應的SOLID164單元為534 050，計算時間設定為20 ms。所建立的有限元數值模型如圖4所示。

圖3 幾何模型示意圖Fig.3 Sketch of geometric model

圖4 復式空心鋼管混凝土柱有限元模型Fig.4 Numericalmodel of duplex hollow CFST column

2.4 數值模擬結果

復式空心鋼管混凝土柱的數值模擬結果如圖5所示。

圖5 復式空心鋼管混凝土柱數值模擬結果（折合距離＝0.14 m／kg1／3）Fig.5 Result of numerical simulation model for duplex hollow CFST column （Scale distance＝0.14 m／kg1／3）

以試驗所測得柱中的彎曲撓度來衡量數值模擬結果的準確性。經比較可以發(fā)現，數值模擬迎爆面柱中彎曲撓度為109.44mm，現場試驗結果為195 mm，兩者在同一數量級內，現場試驗結果大于數值模擬結果，兩者誤差為43.9%。分析兩者誤差產生是由于數值模擬中柱的約束條件以及地面反射條件與現場試驗均有較大差異，且空氣域模擬有限。同時數值模擬是理想狀態(tài)，而現場試驗的條件復雜，加之爆炸試驗采樣率有限。文獻［21］研究表明：當折合距離在0.1 m／kg1／3時，兩者的誤差會達到90%左右，且數值模擬結果均會小于實際爆炸試驗結果。因此綜合考慮上述各個方面因素，可以認為針對爆炸試驗的數值模擬結果是合理的。

為了掌握柱在爆炸荷載下的動態(tài)響應情況，從數值模擬結果中提取柱迎爆面柱底、四分之一柱高、柱中、四分之三柱高和柱頂五個位置的柱位移時程曲線如圖6所示。圖中位移負值表明柱發(fā)生變形的方向與爆炸沖擊波傳播方向一致。

從圖6中可以看出，柱迎爆面五個位置的變形均在6 ms左右達到穩(wěn)定。其中，柱迎爆面變形由于約束和爆心位置兩者的共同影響呈現柱中變形最大，兩端變形最小的特征。從曲線的斜率來看，迎爆面五個位置變形均在2～4 ms的時間段范圍內急劇增加，而后增加減緩趨于穩(wěn)定。圖6（a）中迎爆面柱底位移時程曲線在0.5ms左右的突變分析是由于剛性地面反射所致。將五個位置的最終塑性變形提取出來如表2所示。從表2中的數據可以看出，柱迎爆面的塑性變形量關于爆心近似成上下對稱，同時柱底和柱頂均有少量的位移量，這與實際試驗中實際構件的變形情況是一致的。

圖6 柱迎爆面不同位置單元位移時程曲線Fig.6 Displacement-time curve of different position element on the column surface facing the explosive

表2 柱迎爆面不同位置最終變形量表Tab.2 Displacement value of different position element on the column surface facing the exp losive

在數值模擬結果中提取復式空心鋼管混凝土柱迎爆面柱底、柱中和柱頂三個位置單元的壓力時程曲線如圖7所示。

從圖7中可以看出，在50 kg藥量，爆心距為0.5 m，折合距離為0.14 m／kg1／3的條件下，柱迎爆面柱底、柱中和柱頂的壓力峰值分別為18.72 MPa，176.05 MPa 和4.17 MPa。柱中的壓力峰值分別為柱底和柱頂的9.4倍和42.2倍。通過比對時間數據可以發(fā)現迎爆面上、中、下三個位置壓力峰值產生的時間分別為0.3 ms，0.2 ms和0.3 ms，這說明爆炸沖擊波到達迎爆面柱底、柱中和柱頂的時間幾乎是同時的，而柱底的壓力峰值大于柱頂的原因是因為在模型中添加了剛性面以模擬地面反射的緣故。

在迎爆面正壓持續(xù)時間方面，柱頂、柱中和柱底三個位置的正壓持時分別為2.5ms、2.2ms及1.3ms。因此綜合考慮上述三個位置的壓力峰值可以得出結論：爆炸沖擊波對迎爆面柱中造成的破壞最為強烈、柱底次之、柱頂最小。

圖7 柱迎爆面不同位置單元壓力時程曲線Fig.7 Pressure-time curve of different position element on the column surface facing the explosive

3 復式空心鋼管混凝土柱的損傷評估

3.1 損傷評估準則

目前針對柱構件損傷評估準則可以分為柱中最大位移、最大應力、最大應變以及柱豎向剩余承載力等［10］。確定損傷評估準則的主要原則是：與準則相關的復式空心鋼管混凝土柱的整體特性應較容易通過試驗或數值模擬方法得到，同時又便于針對實際工程開展。因此依據《石油化工控制室抗爆設計規(guī)范》（GB50779－2012）中5.6.3規(guī)定：在爆炸荷載作用下，柱構件支座處彈塑性轉角允許值為2°。結合上述柱構件損傷評估準則要求，可以將控制柱構件支座轉角等效為控制柱中塑性變形量，即柱中撓度不能超過柱計算高度的1／60。本文建立基于固端約束下復式空心鋼管混凝土柱中撓度的超壓－沖量損傷評估準則。復式空心鋼管混凝土柱損傷評估準則可以簡單表述為：柱中撓度小于柱計算高度的1／60，判定柱為安全；柱中撓度大于等于柱計算高度的1／60，判定柱為危險。柱中撓度示意圖如圖8所示。

3.2 復式空心鋼管混凝土柱P－I曲線建立

以建立的爆炸荷載下復式空心鋼管混凝土柱有限元模型為研究對象進行多次有限元數值模型的調整。每一次數值模擬后，均選擇柱迎爆面幾何中心的柱中單元，得到該單元的位移時程曲線。定義柱中撓度為f，柱計算高度為L，將進行了大量數值模擬試算后所得的壓力和沖量數據點繪制在P－I平面內，從數據點中找出符合模型柱中撓度f＝L／60損傷等級分界點的超壓沖量組合臨界值，得到本文中復式空心鋼管混凝土柱損傷等級分界線，即P－I曲線，如圖9所示。

圖8 柱中撓度示意圖Fig.8 Schematic view of column deflection

圖9中的P－I曲線的含義是：曲線將P－I空間劃分為兩個部分，若P、I組合的數據點落在該曲線的左下方，則表明復式空心鋼管混凝土柱迎爆面柱中撓度沒有超過L／60，變形仍在安全范圍內，可以判定該柱為安全；若P、I組合的數據點落在該曲線的右上方，則表明復式空心鋼管混凝土柱迎爆面柱中撓度超過L／60，變形在危險范圍內，可以判定該柱為危險；若P、I組合的數據點落在該曲線上，則表明處于臨界狀態(tài)，基于安全考慮應判定該柱為危險。例如依據本文爆炸試驗，進行數值模擬得到的柱中P＝176.05 MPa、I＝22.70 MPa·ms，該數據點落在P－I曲線的右上方，則可以判定該復式空心鋼管混凝土柱危險。

3.3 P－I曲線擬合公式

為了使經過大量數值模擬后獲得的復式空心鋼管混凝土柱P－I曲線具有明確的數學含義，對圖9中的P－I曲線進行數學表達式擬合。通過對圖9中的數據進行分析可以發(fā)現，I和P近似成自然對數關系，因此經過反復嘗試確定出P－I曲線可以用式（3）進行擬合表達。

式中，P為爆炸沖擊波峰值超壓，I為爆炸沖擊波正向階段沖量，A、B、C為三個實常數，與復式空心鋼管混凝土柱的損傷情況有關。

公式（3）擬合效果如圖10所示。

由圖10可以看出，采用本文推導擬合出的P－I曲線公式所繪曲線與數值模擬曲線基本吻合良好，因此，可以將公式（3）推廣到一般情況，建立以迎爆面柱中撓度為指標的復式空心鋼管混凝土柱損傷評估準則。例如以迎爆面柱中撓度為L／60＝30 mm作為判定指標，則公式（3）可以變換為式（4）：

此時A＝180.5；B＝92.08；C＝－0.5。

圖9 復式空心鋼管混凝土柱P－I曲線Fig.9 P－I curve of duplex hollow CFST column

圖10 公式擬合P－I曲線與數值模擬P－I曲線對比Fig.10 Comparison between fitting P－I curve by formula and P－I curve by numerical simulation

公式（4）適用于本文設計參數的、固端約束條件下復式空心鋼管混凝土柱在爆炸荷載下的損傷評估。其他不同約束條件和設計參數的復式空心鋼管混凝土柱可以依據文中所歸納的方法予以相應推導。

4 結 論

本文采用爆炸試驗和數值模擬的方法研究爆炸荷載下復式空心鋼管混凝土柱的抗爆性能及其損傷評估，主要結論如下：

（1）加強柱端部節(jié)點強度與構件整體強度協調可以提升復式空心鋼管混凝土柱的抗爆性能。

（2）當爆心與柱中等高時，爆炸沖擊波對迎爆面柱中造成的破壞最為強烈、柱底次之、柱頂最小。

（3）依據建立的復式空心鋼管混凝土柱超壓—沖量（P－I）損傷評估準則可以進行爆炸荷載下不同計算高度的、基于迎爆面柱中撓度固端約束復式空心鋼管混凝土柱損傷評估。評價爆炸荷載下復式空心鋼管混凝土柱的安全性能。

［1］蔡紹懷.現代鋼管混凝土結構［M］.北京：人民交通出版社，2003.

［2］鐘善桐.鋼管混凝土結構［M］.3版.北京：清華大學出版社，2003.

［3］王志浩，成戎.復合方鋼管混凝土短柱的軸壓承載力［J］.清華大學學報：自然科學版，2005，45（12）：1596－1599.WANG Zhi-hao，CHENG Rong.Axial bearing capacity of composite-sectioned square concrete-filled steel tubes［J］.Journal ofTsinghua University：Science and Technology，2005，45（12）：1596－1599.

［4］黃宏，陶忠，韓林海.圓中空夾層鋼管混凝土柱軸壓工作機理研究［J］.工業(yè)建筑，2006，36（11）：11－14.

HUANG Hong，TAO Zhong，HAN Lin-hai.Mechanism of concrete-filled double-skin steel tubular columns subjected to axial compression［J］.Industrial Construction，2006，36 （11）：11－14.

［5］Hao H，Cheong H K，Cui S J.Numerical study of dynamic buckling of steel columns subjected to underground explosion ［J］.Engineering Materials，2002，233（1）：211－216.

［6］Cui SJ，Cheong H K，Hao H.Elastic-plastic dynamic response and buckling of steel columns under strong vertical ground motion［J］.Engineering Materials，236（1）：217－222.

［7］李國強，孫建運，王開強.爆炸沖擊荷載作用下框架柱簡化分析模型研究［J］.振動與沖擊，2007，26（1）：8－11.LIGuo-qiang，SUN Jian-yun，WANG Kai-qiang.Research on a simplified frame column model to resist blast load［J］.Journal of Vibration and Shock，2007，26（1）：8－11.

［8］馬云玲，謝冰，何大治，等.圓鋼管混凝土爆炸響應的數值分析［J］.武漢理工大學學報，2010，32（14）：120－124.

MA Yun-ling，XIE Bing，HE Da-zhi，et al.Numerical analysis on blasting response of circular concrete-filled steel tubes［J］.Journal of Wuhan University of Technology，2010，32（14）：120－124.

［9］吳賽.爆炸荷載下復式鋼管混凝土柱動力響應研究［D］.西安：長安大學，2012.

［10］師燕超.爆炸荷載作用下鋼筋混凝土結構的動態(tài)響應行為與損傷破壞機理［D］.天津：天津大學，2009.

［11］孫建運.爆炸沖擊荷載作用下鋼骨混凝土柱性能研究［D］.上海：同濟大學，2006.

［12］Fujikura S，Bruneau M.Experimental investigation of seismically resistant bridge piers under blast loading［J］.Journal of Bridge Engineering，2011，16（1）：63－71.

［13］Fujikura S，Bruneau M，Diego lopez-garcia.Experimental investigation of multihazard resistant bridge piers having concrete-filled steel tube under blast loading［J］.Journal of Bridge Engineering，2008，13（6）：586－594.

［14］程寶坪.深圳賽格廣場鋼管混凝土柱－鋼結構施工特點［J］.施工技術，2000，29（6）：2－5.

CHENG Bao-ping.Constructional features of steel-pipe concrete column and steel structure of shenzhen saiger plaza ［J］.Construction Technology，2000，29（6）：2－5.

［15］Randers Pehrson G，Bannister K A.Airblast loading model for DYNA2D and DYNA3D［R］.Army Research Laboratory，Rept.ARL-TR－1310，1997.

［16］石少卿，康建功，汪敏，等.ANSYS／LS-DYNA在爆炸與沖擊領域內的工程應用［M］.北京：中國建筑工業(yè)出版社，2011.

［17］白金澤.LS-DYNA3D理論基礎與實例分析［M］.北京：科學出版社，2005.

［18］盧紅琴，劉偉慶.空中爆炸沖擊波的數值模擬研究［J］.武漢理工大學學報，2009，31（19）：105－108.

LU Hong-qin，LIU Wei-qing.Research on numerical simulation of blast wave in air［J］.Journal of Wuhan University of Technology，2009，31（19）：105－108.

［19］李順波，東兆星，齊燕軍，等.爆炸沖擊波在不同介質中傳播衰減規(guī)律及數值模擬［J］.振動與沖擊，2009，28（7）：115－117.

LIShun-bo，DONG Zhao-xing，QIYan-jun，et al.Numerical simulation on decay spread of blasting shock wave in different media［J］.Journal of Vibration and Shock，2009，28（7）：115－117.

［20］Shi Yan-chao，Li Zhong-xian，Hao Hong.Mesh size effect in numerical simulation of blastwave propagation and interaction with structures［J］.Transactions of Tianjin University，2008，14（6）：396－402.

［21］楊鑫，石少卿，程鵬飛.空氣中TNT爆炸沖擊波超壓峰值的預測及數值模擬［J］.爆破，2008，25（1）：15－18.

YANG Xin，SHI Shao-qing，CHEN Peng-fei.Forecast and simulation of peak overpressure of TNT explosion shock wave in the air［J］.Blasting，2008，25（1）：15－18.

Explosion-resistance behavior and damage assessment of a dup lex hollow CFST column subjected to blast loading

CUIYing1，3，ZHAO Jun-hai2，ZHANG Chang-guang2，SUN Shan-shan2，CHEN Bing1

（1.School of Mechanical Engineering，Xi'an Shiyou University，Xi'an 710065，China；2.School of Civil Engineering，Chang'an University，Xi'an 710061，China；3.School of Aeronautics，Northwestern Polytechnical University，Xi'an 710072，China）

Test analysis and numerical simulationmethod were integrated to evaluate dynamic response and damage of a duplex hollow CFST column subjected to blast loading.A pressure-impulse damage evaluation criterion for the duplex hollow CFST column was defined by analyzing the test data.The results showed that under the blast loading for the test condition of the converted distance of 0.14m／kg1／3，the plastic bending deformation occurs on the surface facing the explosion of the duplex hollow CFST column；the peak pressure in the middle of the column is the highest，the peak pressure at the bottom of the column is higher than that at the top of the column，the latter is the lowest；with the analysis of the duration of the positive pressure，the damage of themiddle of the column is themost serious；so，it is important to coordinate the ends'strength and the whole column's strength in order to enhance the explosion-resistance performance of the duplex hollow CFST column.Finally，a pressure-impulse damage criterion based on the deflection of the fixed-end duplex hollow CFST column was defined，and amathematical formula to generate pressure-impulse diagram（P-I curve）was also established.

explosion mechanics；duplex hollow CFST column；explosion test；numerical simulation；pressureimpulse diagram

TU398+.9

A

10.13465／j.cnki.jvs.2015.21.033

國家自然科學基金項目（41202191）；教育部博士點基金項目（20110205130001）；中國博士后科學基金面上項目（2012M520079、2013T60868）；陜西省自然科學基金項目（2011JM7002）；陜西省教育廳專項科研計劃項目（15JK1501）

2014－12－31 修改稿收到日期：2015－03－19

崔瑩男，博士后，講師，1979年生