楊濤春 羅堯治

摘 要： 爆炸超壓是描述爆炸荷載的重要指標(biāo)，不同方法的超壓計(jì)算結(jié)果具有較高的離散性。通過分析試驗(yàn)數(shù)據(jù)圖表擬合確定入射超壓與反射超壓關(guān)系的反射系數(shù)計(jì)算式，搜集轉(zhuǎn)化得到大量爆炸超壓理論計(jì)算公式與爆炸試驗(yàn)數(shù)據(jù)，從而對(duì)爆炸超壓在不同比例距離下的分布特征進(jìn)行分析。研究結(jié)果表明：當(dāng)比例距離小于0.5? m/kg1/3時(shí)，爆炸超壓概率密度服從指數(shù)分布;當(dāng)比例距離大于0.5? m/kg1/3時(shí)，爆炸超壓概率密度服從正態(tài)分布。當(dāng)比例距離小于0.5? m/kg1/3時(shí)，爆炸超壓變異系數(shù)達(dá)最大值1，比例距離在1.5～6.0? m/kg1/3之間時(shí)，變異系數(shù)較小，在0.13～0.2? m/kg1/3之間;反射超壓變異系數(shù)較入射超壓略大。依據(jù)不同比例距離下爆炸超壓分布期望數(shù)據(jù)，擬合得到爆炸超壓的計(jì)算公式與具有95%保證率條件下的爆炸超壓分布范圍計(jì)算公式。

關(guān)鍵詞： 爆炸超壓;比例距離;變異系數(shù);分布特征

中圖分類號(hào)：TU312? ? 文獻(xiàn)標(biāo)志碼：A? ?文章編號(hào)：2096-6717（2020）02-0115-10

Distribution characteristic and determination of overpressure for blast resistant study of buildings

Yang Taochun1，2， Luo Yaozhi1

（1.College of Civil Engineering and Architecture， Zhejiang University， Hangzhou 310027，P.R.China; 2.School of Civil Engineering and Architecture， University of Jinan， Jinan 250022， P.R.China）

Abstract： Explosion overpressure is an important indicator to describe blast load.And the results obtained by different methods have high discreteness.Based on this， the reflection coefficient formula of the relationship between incident overpressure and reflection overpressure is determined firstly by analyzing the experimental data chart fitting.Then a large number of theoretical calculation formulas and explosion test data were collected and transformed to analyze the distribution characteristics of explosion overpressure at different proportional distances.The results show that the probability density of explosion overpressure obeys exponential distribution when the scaled distance less than 0.5? m/kg1/3.And obeys normal distribution when the scaled distance greater than 0.5 m/kg1/3.The variation coefficient reaches the maximum value 1 when the scaled distance less than 0.5 m/kg1/3.While variation coefficient is in the range of 0.13 to 0.2 when the scaled distance within the range of 1.5 m/kg1/3 to 6 m/kg1/3.The variation coefficient of reflect overpressure is a little bigger than that of incident overpressure.According to the expected data of explosion overpressure distribution at different proportional distances， the calculation formula of explosion overpressure and the calculation formula of explosion overpressure distribution range with 95% guarantee rate are obtained by fitting.

Keywords： explosion overpressure; scale distance; variation coefficient; distribution characteristic

隨著恐怖襲擊事件和偶然爆炸事故的不斷發(fā)生，建筑結(jié)構(gòu)防爆、抗爆研究已成為土木工程領(lǐng)域的熱點(diǎn)問題，特別是從“9·11”事件以來，世界多國(guó)學(xué)者已開展很多相關(guān)建筑結(jié)構(gòu)的抗爆研究工作，如爆炸荷載、結(jié)構(gòu)動(dòng)力響應(yīng)、破壞模式及簡(jiǎn)化計(jì)算、抗爆分析方法及抗爆加固措施等，得出很多非常有意義的結(jié)論，對(duì)指導(dǎo)工程結(jié)構(gòu)防爆、抗爆安全有重要參考價(jià)值。

在搜集的試驗(yàn)數(shù)據(jù)中，既有垂直入射數(shù)據(jù)，也有非垂直入射數(shù)據(jù)，為分析數(shù)據(jù)的統(tǒng)一和嚴(yán)謹(jǐn)性，僅保留垂直入射的數(shù)據(jù)，即上述數(shù)據(jù)共計(jì)125組。而在這125組試驗(yàn)數(shù)據(jù)中，有的僅有入射超壓，有的僅有反射超壓，入射超壓與反射超壓同時(shí)存在的只有11組，對(duì)于同批試驗(yàn)，有22組數(shù)據(jù)存在較大差異，在分析中未考慮，因此，實(shí)際使用的試驗(yàn)統(tǒng)計(jì)數(shù)據(jù)有103組。

分析爆炸超壓計(jì)算公式發(fā)現(xiàn)，入射超壓的計(jì)算公式主要有兩種形式，如式（4）、式（5）所示，在每個(gè)公式中，均存在3個(gè)待定系數(shù)。其中，式（4）通過爆炸力學(xué)理論求得，并通過試驗(yàn)確定相關(guān)系數(shù)，此方法在已有入射超壓的計(jì)算公式中得到更多應(yīng)用，式（5）主要通過試驗(yàn)數(shù)據(jù)擬合而得，在入射超壓的計(jì)算公式中也有一定應(yīng)用。對(duì)于統(tǒng)計(jì)的僅有入射超壓或反射超壓的試驗(yàn)數(shù)據(jù)，通過反射系數(shù)式（2）和入射超壓計(jì)算公式可得到該組試驗(yàn)條件下的相關(guān)入射超壓公式和反射超壓相關(guān)數(shù)據(jù)，為分析得到爆炸超壓的概率分布特征提供統(tǒng)計(jì)樣本。

P0= a Z + b Z2 + c Z3? （4）

P0=d·Ze+f （5）

式中：Z為比例距離。

2 超壓分布特征

基于統(tǒng)計(jì)的爆炸超壓數(shù)據(jù)，對(duì)不同比例距離下的離散超壓值進(jìn)行分析，畫出對(duì)應(yīng)不同比例距離的超壓分布直方圖，結(jié)合超壓直方圖分布特征，并通過K-S檢驗(yàn)和Lilliefors檢驗(yàn)分別開展不同比例距離下的超壓概率分布擬合優(yōu)度檢驗(yàn)。結(jié)果表明，當(dāng)比例距離小于0.5 m/kg1/3時(shí)，入射超壓和反射超壓分布均服從指數(shù)分布，當(dāng)比例距離大于0.5 m/kg1/3時(shí)，超壓分布均服從正態(tài)分布。在確定概率分布模型基礎(chǔ)上，計(jì)算得到不同比例距離條件下的超壓均值和標(biāo)準(zhǔn)差的極大似然估計(jì)值，同時(shí)，得到超壓均值在95%保證率條件下的置信區(qū)間，如表2、表3所示，從而得到不同比例距離下的超壓分布概率密度曲線，如圖4、圖5所示（指數(shù)分布和正態(tài)分布各兩組）。從表2、表3中可以看出，比例距離較小時(shí)，標(biāo)準(zhǔn)差最大，超壓分布越分散。

為對(duì)比不同比例距離條件下入射超壓與反射超壓的分散程度，通過表2、表3中超壓期望和標(biāo)準(zhǔn)差的極大似然估計(jì)值得到超壓分布的變異系數(shù)，如圖6所示。從圖6可以看出，比例距離小于0.5 m/kg1/3時(shí)，超壓變異系數(shù)達(dá)到最大，為指數(shù)分布的常值1;當(dāng)比例距離約在1.5～6 m/kg1/3之間時(shí)，得益于較多的試驗(yàn)數(shù)據(jù)，變異系數(shù)較小，在0.13～0.2之間，且反射超壓的變異系數(shù)較入射超壓略大。

3 爆炸超壓公式確定

根據(jù)表2、表3中超壓期望的極大似然估計(jì)值，可得到入射超壓和反射超壓95%置信區(qū)間上、下限隨比例距離的變化關(guān)系，分別如圖7、圖8所示，因不同比例距離的超壓值相差較大，故將比例距離分3段分別繪制，從圖中也可看出超壓分散程度隨比例距離的變化趨勢(shì)，比例距離越小，超壓分布越分散。取比例距離和超壓值的自然對(duì)數(shù)，再通過最小二乘法對(duì)自然對(duì)數(shù)超壓進(jìn)行多項(xiàng)式擬合，如圖9、圖10所示，入射超壓和反射超壓擬合曲線的回歸系數(shù)均大于0.99，最終得到入射超壓的計(jì)算公式為

ln P0=0.158 3ln2Z-2.342ln Z-0.097 7 （6）

入射超壓95%置信區(qū)間上、下限的計(jì)算公式為

ln P0=0.214 9ln2Z-2.486 7ln Z+0.071 1（上限） （7）

ln P0=0.111 1ln2Z-2.242 8ln Z-0.250 6（下限） （8）

同理，可得到反射超壓的計(jì)算公式為

ln Pr=0.208 7ln2Z-2.926 3ln Z+1.564 7 （9）

反射超壓95%置信區(qū)間上、下限的計(jì)算公式為

ln Pr=0.264 3ln2Z-3.065 2ln Z+1.751 2（上限） （10）

ln Pr=0.162ln2Z-2.833 9ln Z+1.389 8（下限） （11）

其中，比例距離Z在0.1～15 m/kg1/3范圍內(nèi)。

4 結(jié)論

通過搜集文獻(xiàn)獲取大量爆炸超壓的試驗(yàn)與理論數(shù)據(jù)，并從不確定性角度出發(fā)，研究了爆炸超壓的概率分布特征，主要得到以下結(jié)論：

1）爆炸超壓試驗(yàn)數(shù)據(jù)受炸藥類型、當(dāng)量、形狀及試驗(yàn)環(huán)境等因素影響明顯，超壓分布具有較高的離散性，且試驗(yàn)數(shù)據(jù)多以小當(dāng)量炸藥為主，比例距離多集中于0.4～2.0 m/kg1/3之間。

2）基于試驗(yàn)數(shù)據(jù)，針對(duì)垂直入射情況擬合給出反射系數(shù)公式，并得到根據(jù)入射超壓獲取反射超壓的計(jì)算公式。

3）比例距離小于0.5 m/kg1/3時(shí)，爆炸超壓概率密度服從指數(shù)分布;比例距離大于0.5 m/kg1/3時(shí)，爆炸超壓概率密度服從正態(tài)分布。

4）比例距離小于0.5 m/kg1/3時(shí)，爆炸超壓變異系數(shù)達(dá)最大值1;比例距離約在1.5～6 m/kg1/3間時(shí)，變異系數(shù)較小，在0.13～0.2之間;反射超壓變異系數(shù)較入射超壓略大。

5）根據(jù)不同比例距離下爆炸超壓分布期望數(shù)據(jù)，擬合得到爆炸超壓的計(jì)算公式與具有95%保證率條件下的超壓分布范圍計(jì)算公式。

參考文獻(xiàn)：

[1]? ? Department of the Army. Fundamentals of protective design for conventional weapons：TM 5-855-1 [M]. Department of the Army， Washington， DC， USA， 1986.

[2]? ? Department of the Army. Structures to resist the effects of accidental explosion：TM 5-1300 [M]. Department of the Army， 1969.

[3]? ? Department of Defense Explosives Safety Board. Structures to Resist the Effects of Accidental Explosions：UFC 3-340-02 [M]. Department of Defense Explosives Safety Board （DDESB）， 2008.

[4]? ? Design and assessment of buildings subjects to blast loads：CSA/S 850-12 [S]. Ontario， Canada：Canadian Standards Association， 2012.

[5]? ?爆破安全規(guī)程：GB 6722—2014 [S]. 北京：中國(guó)標(biāo)準(zhǔn)出版社， 2014.

Safety regulations for blasting：GB 6722-2014 [S]. Beijing：China Standard Press， 2014. （in Chinese）

[6]? ?人民防空地下室設(shè)計(jì)規(guī)范：GB 50038—2005 [S]. 北京：中國(guó)計(jì)劃出版社， 2005.

Code for design of civil air defence basement：GB 50038-2005 [S]. Beijing：China Planning Press， 2014. （in Chinese）

[7]? ? LOW H Y， HAO H. Reliability analysis of reinforced concrete slabs under explosive loading [J]. Structural Safety， 2001， 23（2）：157-178.

[8]? ?BOGOSIAN D， FERRITTO J， SHI Y J. Measuring uncertainty and conservatism in simplified blast models [C] // 30th Explosives Safety Seminar， August 2002， Atlanta， Georgia， USA.

[9]? ?李忠獻(xiàn)，路建輝，師燕超，等. 不確定爆炸荷載作用下鋼梁的可靠度分析[J]. 工程力學(xué)， 2014， 31（4）：112-118， 133.

LI Z X， LU J H， SHI Y C， et al. Reliability analysis of steel beam under uncertain blast loads [J]. Engineering Mechanics， 2014， 31（4）：112-118， 133. （in Chinese）

[10]? ?路建輝. 不確定爆炸荷載作用下鋼梁的可靠度分析[D]. 天津：天津大學(xué)， 2012.

LU J H. Reliability analysis of steel beam under uncertain blast loads [D]. Tianjin：Tianjin University， 2012. （in Chinese）

[11]? ?李忠獻(xiàn)， 任其武， 師燕超， 等. 重要建筑結(jié)構(gòu)抗恐怖爆炸設(shè)計(jì)爆炸荷載取值探討[J]. 建筑結(jié)構(gòu)學(xué)報(bào)， 2016， 37（3）：51-58.

LI Z X， REN Q W， SHI Y C， et al. Research on blast load value in design of important building structures against terrorist explosions [J]. Journal of Building Structures， 2016， 37（3）：51-58. （in Chinese）

[12]? ?MILLS C A. The design of concrete structure to resist explosions and weapon effects [C] // Proceedings of the 1st International Conference on Concrete for Hazard Protections， Edinburgh， UK：1987：61-73.

[13]? ? BRODE H L. Blast wave from a spherical charge [J]. Physics of Fluids， 1959， 2（2）：217.

[14]? ?BRODE H L. Numerical solutions of spherical blast waves [J]. Journal of Applied physics. 1955， 26（6）：766-775.

[15]? ?BRODE H L. A calculation of the blast wave from a spherical charge of TNT [M]. Austin：University of Texas at Austin， 1961.

[16]? ? HENRYCH J. 爆炸動(dòng)力學(xué)及其應(yīng)用[M]. 熊建國(guó). 譯. 北京：科學(xué)出版社， 1987.

HENRYCH J. The Dynamics of explosion and its use [M]. Translated by XIONG J G. Beijing：Science Press， 1987. （in Chinese）

[17]? ?HENRYCH J， ABRAHAMSON G R. The dynamics of explosion and its use [M]. Amsterdam and New York：Elsevier Scientific Pub. Co.， 1979.

[18]? ?KINNEY G F， GRAHAM K J. Explosive shocks in air [M]. Berlin：Springer-Verlag， 1985.

[19]? ?NEWMARK N M， HANSEN R J. Design of blast resistant structures [M]. New York：Shock and Vibration Handbook. Harris and Crede. McGraw-Hill， 1961：3.

[20]? ?BAKER W E， COX P， WESTINE P， et al. Explosion hazards and evaluation [M]. Amsterdam：Elsevier Scientific Pub. Co.， 1983.

[21]? ? WU C Q， HAO H. Modeling of simultaneous ground shock and airblast pressure on nearby structures from surface explosions [J]. International Journal of Impact Engineering， 2005， 31（6）：699-717.

[22]? ?NAUMYENKO I A， PETROVSKY I G. The shock wave of a nuclear explosion [M]. Moscow：вoен. издат. мии. сборокы СССР， 1956：183-184.

[23]? ?SADOVSKII M A. Mechanical effects of air shock waves from explosions according to experiments [M]. Moscow：физика взрыва， No.1， AH СССР， 1952.

[24]? ?葉曉華. 軍事爆破工程[M]. 北京：解放軍出版社，1999.

YE X H， Military blasting engineering [M]. Beijing：the Peoples Liberation Army Press，1999.（in Chinese）

[25]? ? SIDDIQUI J I， AHMAD S. Impulsive loading on a concrete structure [J]. Proceedings of the Institution of Civil Engineers - Structures and Buildings， 2007， 160（4）：231-241.

[26]? ? AHMAD S， ELAHI A， IQBAL J， et al. Impulsive loading on reinforced concrete wall [J]. Proceedings of the Institution of Civil Engineers - Structures and Buildings， 2013， 166（3）：153-162.

[27]? ?IQBAL J， AHMAD S. Improving safety provisions of structural design of containment against external explosion [C] // International Conference on Opportunities and Challenges for Water Cooled Reactors in the 21st Century， International Atomic Energy Agency（IAEA）， 2011.

[28]? ? HELD M. Blast waves in free air [J]. Propellants， Explosives， Pyrotechnics， 1983， 8（1）：1-7.

[29]? ? GHANI RAZAQPUR A， TOLBA A， CONTESTABILE E. Blast? loading response of reinforced concrete panels reinforced with externally bonded GFRP laminates [J]. Composites Part B：Engineering， 2007， 38（5/6）：535-546.

[30]? ? WANG H W， WU C Q， ZHANG F R， et al. Experimental study of large-sized concrete filled steel tube columns under blast load [J]. Construction and Building Materials， 2017， 134：131-141.

[31]? ? LI J， WU C Q， HAO H， et al. Post-blast capacity of ultra-high performance concrete columns [J]. Engineering Structures， 2017， 134：289-302.

[32]? ?金何偉，劉中憲，劉申永，等. 鋼管超高強(qiáng)鋼纖維混凝土柱抗爆性能試驗(yàn)研究[J]. 建筑結(jié)構(gòu)， 2016， 46（4）：45-49.

JIN H W， LIU Z X， LIU S Y， et al. Experimental study of ultra-high performance fiber reinforced concrete filled steel tube columns under blast loading [J]. Building Structure， 2016， 46（4）：45-49. （in Chinese）

[33]? ? LIANG X X， WANG Z Q， WANG R N. Deformation model and performance optimization research of composite blast resistant wall subjected to blast loading [J]. Journal of Loss Prevention in the Process Industries， 2017， 49：326-341.

[34]? ?李國(guó)強(qiáng)，瞿海雁，楊濤春，等. 鋼管混凝土柱抗爆性能試驗(yàn)研究[J]. 建筑結(jié)構(gòu)學(xué)報(bào)， 2013， 34（12）：69-76.

LI G Q， QU H Y， YANG T C， et al. Experimental study of concrete-filled steel tubular columns under blast loading [J].Journal of Building Structures， 2013， 34（12）：69-76. （in Chinese）

[35]? ? KIM J H J， KIM S B， CHOI J K. Experiment study on blast loading response of FRP-retrofitted RC slab structures [C] // The Second Official International Conference of International Institute for FRP in Construction for Asia-Pacific Region， 2009.

[36]? ? BARBHUIYA S， AHMAD S， ELAHI A， et al. Experimental study of masonry wall exposed to blast loading [J]. Materiales de Construcción， 2014， 64（313）：e007.

[37]? ? RITCHIE C B， PACKER J A， SEICA M V， et al. Behavior of steel rectangular hollow sections subject to blast loading [J]. Journal of Structural Engineering， 2017， 143（12）：04017167.

[38]? ? AUNE V， FAGERHOLT E， HAUGE K O， et al. Experimental study on the response of thin aluminium and steel plates subjected to airblast loading [J]. International Journal of Impact Engineering， 2016， 90：106-121.

[39]? ? TABATABAEI Z S， VOLZ J S， BAIRD J， et al. Experimental and numerical analyses of long carbon fiber reinforced concrete panels exposed to blast loading [J]. International Journal of Impact Engineering， 2013， 57：70-80.

[40]? ? CHOI J H， CHOI S J， CHO C M， et al. Experimental evaluation of Bi-directionally unbonded prestressed concrete panel blast resistance behavior under blast loading scenario [J]. Journal of the Korea Concrete Institute， 2016， 28（6）：673-683.

[41]? ? CHOI J H， CHOI S J， KIM J H J， et al. Evaluation of blast resistance and failure behavior of prestressed concrete under blast loading [J]. Construction and Building Materials， 2018， 173：550-572.

[42]? ? SAFARI H， ZAMANI J， KHALILI R， et al. Experimental， theoretical， and numerical studies on the response of square plates subjected to blast loading [J]. The Journal of Strain Analysis for Engineering Design， 2011， 46（8）：805-816.

[43]? ? ZHANG F R， WU C Q， ZHAO X L， et al. Experimentalstudy of CFDST columns infilled with UHPC under close-range blast loading [J]. International Journal of Impact Engineering， 2016， 93：184-195.

[44]? ? HAO Y F， HAO H， SHI Y C， et al. Field testing of fence type blast wall for blast load mitigation [J]. International Journal of Structural Stability and Dynamics， 2017， 17（9）：1750099.

[45]? ? JACINTO A C， AMBROSINI R D， DANESI R F. Experimental and computational analysis of plates under air blast loading [J]. International Journal of Impact Engineering， 2001， 25（10）：927-947.

[46]? ?崔瑩，趙均海，張常光，等. 爆炸沖擊波在鋼管混凝土柱表面壓力分布試驗(yàn)研究及數(shù)值模擬[J]. 北京工業(yè)大學(xué)學(xué)報(bào)， 2014（12）：1828-1836.

CUI Y， ZHAO J H， ZHANG C G， et al. Research on pressure distribution of blast wave on the surface of CFST column based on explosion test and numerical simulation [J]. Journal of Beijing University of Technology， 2014， 40（12）：1828-1836. （in Chinese）

（編輯 胡玥）

收稿日期：2019-05-15

基金項(xiàng)目：? 國(guó)家自然科學(xué)基金（51608229）;山東省高?？蒲邪l(fā)展計(jì)劃（J18KA206）

作者簡(jiǎn)介：? 楊濤春（1983- ），男，副教授，博士，主要從事爆炸災(zāi)害分析評(píng)估研究，E-mail：yangtaochun@126.com。

Received： 2019-05-15

Foundation items：? National Natural Science Foundation of China （No.51608229）; Project of Shandong Province Higher Educational Science and Technology Program （No. J18KA206）.

Author brief：? Yang Taochun （1983- ）， associate professor， PhD， main research interest： anlysis and assessment of explosion disaster， E-mail： yangtaochun@126.com.