范鵬賢，王明洋，2，馮淑芳，王德榮，李 杰，2

（1．解放軍理工大學爆炸沖擊防災減災國家重點實驗室，南京 210007；2．南京理工大學機械工程學院，南京 210094；3．沈陽軍區(qū)司令部工程科研設計所，沈陽 110162）

直墻拱頂?shù)叵陆Y(jié)構受爆炸地震波作用的動力響應

范鵬賢1，王明洋1，2，馮淑芳3，王德榮1，李 杰1，2

（1．解放軍理工大學爆炸沖擊防災減災國家重點實驗室，南京 210007；2．南京理工大學機械工程學院，南京 210094；3．沈陽軍區(qū)司令部工程科研設計所，沈陽 110162）

研究有彈性墊層的地下直墻拱頂結(jié)構在爆炸地震波作用下動力響應，采用共同變形理論及矩陣力法考慮襯砌與圍巖的相互作用建立爆炸地震波作用的地下防護結(jié)構動力響應快速計算方法。分析圍巖等級、墊層參數(shù)及荷載作用形式等對典型直墻圓拱結(jié)構動力響應影響。結(jié)果表明，圍巖等級對結(jié)構動力響應具有顯著影響，圍巖質(zhì)量越好結(jié)構位移、內(nèi)力越小；墊層有延遲響應、降低響應峰值作用，荷載作用形式也有一定影響。

爆炸地震波；防護工程；直墻拱頂結(jié)構；動力響應

隨新型鉆地武器打擊能力的提高，地下爆炸效應已對防護工程結(jié)構構成嚴重威脅。直墻拱頂?shù)叵陆Y(jié)構為目前防護工程最常用工程結(jié)構形式。開展直墻拱頂?shù)叵陆Y(jié)構在爆炸地震波作用下動力響應問題研究具有重要意義。利用試驗、理論分析及數(shù)值模擬等對爆炸荷載作用的地下結(jié)構動力響應進行廣泛研究［1－6］。我國的此項研究開展較晚，且主要集中于圓形結(jié)構，對直墻拱頂結(jié)構試驗研究［7］及數(shù)值模擬成果［8－10］較少，難以滿足工程設計與評估需求。

范鵬賢等［11］在矩陣力法理論基礎上用集中質(zhì)量方法將圓形隧道離散為有限自由度體系，將爆炸地震波作用視為作用于集中質(zhì)量的動力建立地下圓形結(jié)構整體動力響應快速計算方法。本文采用矩陣力法，考慮襯砌與圍巖間存在彈性墊層影響，建立地下直墻拱頂結(jié)構計算模型。利用算例分析圍巖、墊層、荷載等參數(shù)對結(jié)構動力響應影響規(guī)律，為實際工程設計、評估提供參考。

1 計算模型及結(jié)構運動方程

1.1 計算模型

據(jù)防護工程常用坑道截面形式，選直墻拱頂結(jié)構為研究對象。設結(jié)構為巖土介質(zhì)中鋼筋混凝土封閉的拱形框架，存在彈性墊層時結(jié)構可視為雙層彈性介質(zhì)條件的工作框架。介質(zhì)工作據(jù)彈性半空間模型考慮，結(jié)構承受幾種爆炸波組合的聯(lián)合作用。計算模型為含巖土介質(zhì)、彈性墊層及鋼筋混凝土框架相互作用的三部分，見圖1，圖中：1為鋼筋混凝土框架，2為彈性墊層，3為動載荷。不考慮爆炸地震波在巖土介質(zhì)中傳播特性，將其對彈性墊層及結(jié)構作用視作動荷載加載到墊層上研究結(jié)構動力響應。

為研究結(jié)構強迫振動，用集中質(zhì)量法將結(jié)構簡化成有限自由度的離散體系見圖2。每種地震爆炸波作用均以一系列作用于集中質(zhì)量的動力代替，每個力均用波到達時間toi、波中最大壓力σmi及隨時間變化規(guī)律τi（τ）表示。

圖1 結(jié)構模型簡圖Fig．1 The schematic diagram of structuralmodel

圖2 單元及受力簡圖Fig．2 The elements and their force diagram

圖3 鏈桿法計算簡圖Fig．3 The schematic diagram of Chain-bar Method

1.2 結(jié)構運動方程

對雙層介層采用集中質(zhì)量體系，考慮障礙物柔性強迫振動矩陣形式可寫為

式中：P＝｛Pi｝26為動載荷列向量；M＝［mij］26×26為集中質(zhì)量對角線矩陣；F＝［fij］26×26為集中質(zhì)量單位位移矩陣；C為阻尼矩陣；y＝｛yi｝26，y··＝｛y··｝26分別為集中質(zhì)量的位移及加速度向量。

1．2．1 荷載確定方式

動載荷列向量P以每個波k的列矩陣總和形式確定，即

式中：σkm，αk分別為第k個波到達時最大壓力及角度；k0為波反射系數(shù)：v為圍巖泊松比。

式（3）中向量Fk＝［Fki］26各元素代表第k個波產(chǎn)生的動載荷隨時間變化規(guī)律，可寫為

式中：U為結(jié)構自身變形；y1，ij為彈性半空間變形；y3，ij為墊層變形；J為截面慣性距；A為截面面積；G2為材料剪切模量；γ為截面修正系數(shù)；Fij為未知力Xi作用下j點位移；h3為墊層材料泊松比；E3為墊層材料彈性模量；li為計算段長度。

由此可得確定系數(shù)矩陣B與自由項矩陣D的計算式為

1.3 求解方法

據(jù)以上各式用數(shù)值方法求解有限自由度體系強迫振動。在每個積分步長中對某個參數(shù)進行限制，使原始二階方程式變?yōu)橐来吻蠼獾囊浑A方程組，利用Fortran語言用四階Rung-Kutta方法求解。方程組可寫為

2 圍巖等級對結(jié)構響應影響

2.1 算例參數(shù)設定

矩陣力法對求解地下圓形結(jié)構在爆炸地震波作用下動力響應中已獲得成功應用［11］，此處不重復驗證。為考察圍巖、墊層、荷載等對直墻拱頂結(jié)構動力響應影響，用本文方法計算典型直墻拱頂結(jié)構動力響應。直墻拱頂結(jié)構斷面尺寸寬l＝3．0 m，直墻高h＝1．8 m，拱頂高f＝1．2 m，襯砌結(jié)構厚h0＝0．4 m。據(jù)資料取自由場荷載升壓時間為28．1 ms，正壓作用時間為200 ms，荷載峰值為3．02×105N／m2；鋼筋混凝土結(jié)構彈性模量E＝2．5×1010N／m2，結(jié)構材料密度ρ＝2．55×103kg／m3。如無特殊說明，本文均采用此參數(shù)。

據(jù)《公路隧道設計規(guī)范》，選取典型圍巖進行分析。其物性參數(shù)見表1。

表1 典型圍巖物性參數(shù)Tab.1 Parameters of rocks

2.2 圍巖參數(shù)不同時結(jié)構動力響應

圍巖等級不同時拱頂位移時程曲線見圖4。由圖4看出，荷載與直墻垂直時拱頂產(chǎn)生向下的位移，圍巖等級對隧道位移影響較顯著，圍巖特性越好阻抗比越大，拱頂位移越小，限制隧道位移能力越好；隨圍巖剛度的減小結(jié)構產(chǎn)生位移峰值時間明顯延后。

圖4 拱頂豎向位移時程曲線Fig．4 The time-history curves of displacement of the arch top

圖5 拱頂截面彎矩時程曲線Fig．5 The time-history curves of themoment of the arch top cross-section

圖6 拱頂截面軸力時程曲線Fig．6 The time-history curves of the axial force of the arch top cross-section

圍巖等級不同時拱頂截面彎矩、軸力的時程曲線見圖5、圖6。由兩圖可知，圍巖特性越好襯砌結(jié)構截面承受的內(nèi)力越小。與襯砌拱頂位移時程（圖4）比較知，圍巖等級較差（圍巖5）時襯砌拱頂位移峰值約為圍巖等級較好（圍巖1）的2倍，但相應截面軸力值相差約6倍，而彎矩相差超10倍。由此可見，爆炸地震波對地下結(jié)構內(nèi)力影響遠大于對結(jié)構位移影響。結(jié)構的彎矩、軸力為結(jié)構設計選擇截面及配筋的重要依據(jù)，因此對結(jié)構截面動內(nèi)力進一步研究十分必要。

3 墊層參數(shù)變化對結(jié)構響應影響

在結(jié)構與圍巖間增加軟墊層為防護工程常用隔振減震方法。為說明墊層對結(jié)構動力響應作用，本文模型考慮墊層變形影響，計算不同厚度墊層的結(jié)構動力響應。設襯砌厚度不變，分別計算襯砌厚度h＝0．6 m，0．8 m，1．0 m，1．2 m，1．4 m時襯砌拱頂位移隨時間變化曲線，見圖7。由圖7看出，設置墊層后拱頂位移明顯減小，隨墊層厚度增大襯砌達到位移峰值的時間延遲，峰值減小。此因設置墊層與圍巖間波阻抗不匹配，沖擊波在墊層與圍巖之間產(chǎn)生反射，部分能量被反射，部分能量被墊層吸收，從而減小對結(jié)構的動力作用。而當墊層厚度與襯砌厚度比大于3時，增加墊層厚度對結(jié)構影響效果明顯減小。

圖7 墊層厚度變化時襯砌拱頂位移時程曲線Fig．7 The time-history curves of the displacement of the arch top cross-section under different subcrust conditions

有墊層襯砌拱頂截面環(huán)向應力與無墊層時最大應力值之比的時程曲線見圖8。由圖8看出，襯砌結(jié)構加入彈性墊層后因彈性墊層的消波作用，截面應力明顯減?。?0%～17%）。隨墊層厚度增加影響越明顯，且應力最大值出現(xiàn)時間明顯滯后；但墊層厚度與襯砌厚度比大于3時，應力變化不再明顯?？梢妷|層厚度并非越厚越好。設計中應綜合考慮經(jīng)濟因素及抗爆效果。通常墊層彈性模量遠小于襯砌結(jié)構。墊層較軟時分析表明，墊層彈性模量對墊層抗爆隔振效果影響較小。

圖8 拱頂截面環(huán)向應力與襯砌無墊層時最大應力比值時程曲線Fig．8 The time-history curves of the non-dimensional axial stress of the arch top cross-section under different subcrust conditions

4 荷載作用對結(jié)構響應影響

4.1 爆炸波入射角變化對結(jié)構響應影響

對襯砌、圍巖荷載作用位置不變、爆炸波入射角由0°～90°變化時結(jié)構動力響應進行分析。

爆炸地震波入射角度變化時拱頂位移變化時程曲線見圖9。由圖9看出，入射角度增加拱頂位移越大，入射角大于45°后拱頂位移逐漸變??；入射角θ＝90°時結(jié)構位移最?。蝗肷浣铅龋?5°時對結(jié)構位移影響最大。

圖9 拱頂位移時程曲線Fig．9 The time-history curves of the displacementof the arch top cross-section

圖10 拱頂彎矩時程曲線Fig．10 The time-history curves of the momentof the arch top cross-section

圖11 拱頂軸力時程曲線Fig．11 The time-history curves of the axial force of the arch top cross-section

荷載作用位置變化時拱頂截面彎矩、軸力時程曲線見圖10、圖11。由二圖看出，θ＝0°時荷載對拱頂截面彎矩影響最大；θ＝45°時拱頂產(chǎn)生反向彎矩，但拱頂截面軸力此時達最大。說明此時截面產(chǎn)生的應力最大，即同等荷載對結(jié)構產(chǎn)生的動力響應最大。

4.2 荷載正壓作用時間對結(jié)構響應影響

正壓作用時間對結(jié)構位移、應力亦有影響。增加正壓作用時間結(jié)構位移、應力峰值會增加約10%，但位移、應力出現(xiàn)峰值的時間基本一致。增加荷載作用時間升壓作用時間不變。正壓作用時間增加時結(jié)構位移、應力變化基本一致。

5 結(jié) 論

據(jù)矩陣力法，考慮結(jié)構、墊層及圍巖的相互作用，建立直墻圓拱結(jié)構動力響應計算方法；研究圍巖特性、墊層參數(shù)及荷載作用形式對結(jié)構動力響應影響，結(jié)論如下：

（1）圍巖等級對隧道動力響應影響顯著，圍巖特性越好，拱頂位移、內(nèi)力越小；爆炸地震波對地下結(jié)構內(nèi)力影響遠大于對其位移影響。

（2）鋪設墊層利于減小襯砌的位移及內(nèi)力，增加墊層厚度襯砌達位移峰值時間有一定延遲；墊層厚度并非越厚越好。

（3）荷載作用位置不變、爆炸波入射角的變化對結(jié)構動力響應影響較大。本文所提方法及分析結(jié)果可為防護工程抗爆設計、評估提供參考。

［1］Akiyoloi T，F(xiàn)uschida K．Soil-pipelines interaction through a frictional interface during earthquake［J］．Soil Dynamic and Earthquake Engineering，1984，3（1）：77－81．

［2］NateghiR，Kiany M，Gholipouri O．Control negative effects of blasting waves on concrete of the structures by analyzing of parameters of ground vibration［J］．Tunneling and Underground Space Technology，2009，24（6）：508－515．

［3］沈俊，顧金才，陳安敏，等．巖土工程抗爆結(jié)構模型試驗裝置研制及應用［J］．地下空間與工程學報，2007，3（6）：1077－1132．

SHEN Jun，GU Jin-cai，CHEN An-m in，et al．Evelopment and applications of themodel test apparatus on anti-explosion structures in geotechnical engineering［J］．Chinese Journal of Underground Space and Engineering，2007，3（6）：1077－1132．

［4］辛凱，姜忻良，吳祥云．爆炸荷載作用下兩相飽和土中結(jié)構響應試驗研究［J］．巖石力學與工程學報，2009，28（S2）：4065－4070．

XIN Kai，JIANG Xin-liang，WU Xiang-yun．Experimental study on structure response in two-phase saturated soil under blasting load［J］．Chinese Journal of Rock Mechanics and Engineering，2009，28（S2）：4065－4070．

［5］高盟，高廣運，王瀅，等．均布突加荷載作用下圓柱形襯砌振動響應的解析解［J］．巖石力學工程學報，2010，29（2）：237－242．

GAOMeng，GAO Guang-yun，WANG Ying，etal．Analytical solution on dynamic response of lining subjected to sudden internal uniform loading［J］．Chinese Journal of Rock Mechanics and Engineering，2010，29（2）：237－242．

［6］Feldgun V R，Kochetkov A V，Karinski Y S，et al．Blast response of a lined cavity in a porous saturated soil［J］．International Journal of Impact Engineering，2008，35：953－966．

［7］楊自友，顧金才，陳安敏，等．錨桿長度對洞室抗爆性能影響的模型試驗研究［J］．振動與沖擊，2009，28（1）：24－27．

YANG Zi-you，GU Jin-cai，CHEN An-min，et al．Model experiment study on the influence of anti-blast characteristics of length of rock bolts in wall rock ofunderground cavern［J］．Journal of Vibration and Shock，2009，28（1）：24－27．

［8］國勝兵，王明洋，趙躍堂，等．爆炸地震波作用下地下結(jié)構動力響應數(shù)值分析［J］．世界地震工程，2004，20（4）：137－142．

GUO Sheng-bing，WANG Ming-yang，ZHAO Yue-tang，et al．Dynamic numerical analysis of underground structures under action of exp losion seism ic wave［J］．World Information on Earthquake Engineering，2004，20（4）：137－142．

［9］孫惠香，許金東．爆炸荷載作用下地下結(jié)構與圍巖的動力響應［J］．土木工程學報，2011，44（S）：148－152．

SUN Hui-xiang，XU Jin-dong．Dynamic action study between structure and rock of underground structure subjected to blast load［J］．China Civil Engineering Journal，2011，44（S）：148－152．

［10］鄧春梅，許金余，沈劉軍．裝藥爆炸下地下拱形結(jié)構變形及破壞特征分析［J］．解放軍理工大學學報（自然科學版），2007，8（5）：534－537．

DENG Chun-mei，XU Jin-yu，SHEN Liu-jun．Deformation and damage characteristics analysis of underground arch structure subjected to subsurface blast［J］．Journal of PLA University of Science and Technology（Natural Science Edition），2007，8（5）：534－537．

［11］范鵬賢，王明洋，馮淑芳，等．爆炸地震波作用下深埋圓形隧道的動力響應分析［J］．巖石力學與工程學報，2013，32（4）：671－680．

FAN Peng-xian，WANG Ming-yang，F(xiàn)ENG Shu-fang，et al．Dynamic response of deep-buried circular tunnel to explosion seismic wave［J］．Chinese Journal of Rock Mechanics and Engineering，2013，32（4）：671－680．

［12］王瑞民，華毅杰，武建勛．半空間地下結(jié)構與土相互作用的鏈桿法研究［J］．同濟大學學報，1999，27（4）：495－498．

WANG Rui-m in，HUA Yi-jie，WU Jian-xun．Chain-bar method analysis of soil and underground structure interaction in semi-space［J］．Journal of Tongji University，1999，27（4）：495－498．

［13］Kerr A D．Elastic and viscoelastic foundation models［J］．Journal of Applied Mechanics，1964，31：491－498．

Dynam ic response of underground straight-wall-round-arch structure subjected to explosion seism ic wave

FAN Peng-xian1，WANGMing-yang1，2，F(xiàn)ENG Shu-fang3，WANGDe-rong1，LIJie1，2
（1．State Key Laboratory of Disaster Prevention＆Mitigation of Explosion＆Impact，PLA University of Science and Technology，Nanjing，Jiangsu 210007，China；2．School of Mechanical Engineering，Nanjing University of Science and Technology，Nanjing，Jiangsu 210007，China；3．Institute of Engineering Research and Design，Headquarters of Shenyang Military Area，Shenyang 110162，China）

The dynamic response of underground structure impacted by explosion is of great importance in the designing and evaluating of the protective engineering．The computational method for the dynamic response of straightwall-round-arch structure subjected to explosion seismic wave was investigated on the foundation of common deformation theory and matrix forcemethod，and a rapid computation method was established．By the presented method，the dynamic responses of a typical straight-wall-round-arch structure were studied under different surrounding rock，subcrust and impact loading conditions．The results show that the rock grade，the thickness of subcrust and the loadingmode all have notable influences on the structure's dynamic response．The better quality the rock is of the lower the displacementand the internal force will be．And as the thickness of the subcrust increases，the peak value of the dynamic response decreases and delays．

explosion seismic wave；protective engineering；straight-wall-round-arch structure；dynamic response

TU352；TU 928

：A

10．13465／j．cnki．jvs．2014．22．033

國家自然科學創(chuàng)新研究群體科學基金（51021001）；博士后科學基金面上項目

2013－01－17 修改稿收到日期：2013－11－28

范鵬賢男，博士，講師，1983年生