基于薄板理論的混凝土橋梁近場抗爆方法研究*

胡志堅 劉 超 方建橋

(武漢理工大學交通學院1) 武漢 430063) (襄陽市交通工程質量監督處2) 襄陽 441000)

基于薄板理論的混凝土橋梁近場抗爆方法研究*

胡志堅1)劉 超1, 2)方建橋1)

(武漢理工大學交通學院1)武漢 430063) (襄陽市交通工程質量監督處2)襄陽 441000)

以近場爆炸條件下的混凝土橋梁為研究對象，結合現有公路橋梁多采用梁板式混凝土，其橋面板具有高寬比和高跨比都比較小的特點，將橋面板做整體薄板考慮，基于動力薄板理論提出具體的動態控制方程，開展爆炸荷載作用下混凝土梁板式橋梁動力特性研究，并通過理論分析結果與實測數據的對比分析，驗證了基于薄板理論開展常規梁板式混凝土橋梁抗爆性能研究的有效性與可行性.同時提出合理控制比例距離是實現橋梁結構抗爆設計的有效途徑.

薄板理論；近場爆炸荷載；混凝土橋梁；比例距離

國內外許多學者對結構或構件在爆炸荷載作用下的動力響應分析方法進行了研究，Rong等[1]提出了基于性能的RC梁(柱)抗爆設計方法NES，并通過實驗梁(柱)驗證了該方法的有效性；Fang等[2]開展了運用薄板彎曲理論計算了防護土墻體的抗爆研究，獲得了爆炸荷載作用下防護土墻體的p-I曲線.數值分析方法可以綜合模擬爆炸沖擊波荷載的復雜時程曲線、材料在爆炸沖擊波荷載作用下的復雜本構關系、結構構件的復雜邊界條件、材料的局部損傷以及結構整體動力分析等，是進行結構抗爆分析的有力工具[3].但對于近場爆炸，數值模擬方法計算精度偏低，與試驗實測結果存在明顯差異[4].試驗研究是對理論分析和數值模擬的有效驗證手段，Williamson等[5-6]以1∶2的相似比開展了RC柱的抗爆試驗研究；Fujikura等[7]開展了鋼管混凝土柱式排架的比例模型試驗研究；孫文彬[8]選擇鋼筋混凝土簡支單向板進行爆炸試驗,觀測了鋼筋混凝土板在彈性區域的動態響應和塑性區域的破壞特征.本文擬將橋面板做整體薄板分析，提出相應的求解方法，并通過具體的實橋爆炸測試數據，驗證該方法的有效性，進而研究比例距離變化對橋梁結構爆炸荷載響應的影響.

1 動態響應分析

(1)

(2)

(3)

式中：Mx，My為板側面正應力σx和σy產生的分布彎矩；Mxy為板側面切應力τxy產生的分布轉矩；w(x,y,t)為撓度；ρ為密度；C為板的阻尼；h為板的厚度.

如果將板側面上的應力分量向中面簡化，得到作用在板側面中線上的分布彎矩和轉矩如下關系式.

(4)

(5)

(6)

將等式(4)和式(5)代入式(2)，得

(7)

將式(5)和(6)代入式(3)，得

(8)

最后將式(7)和式(8)代入式(1)，得到一般動態響應方程

(9)

在式(9)中，令q(x,y,t)=0并忽略阻尼項，可以得到如下無阻尼自由振動方程

(10)

假設它有一種如下形式的解

(11)

此處Z(x,y)是振型函數.

再將式(11)代入式(10)，得

(12)

(13)

式(13)左邊是關于(x,y)的函數，而右邊只與t有關.欲使等式兩邊相等，則其必定等于某一常數，用ω2表示.故而式(13)可以化為2個等式，一個是時間函數φ(t);另一個是空間位置函數Z(x,y).

(14)

(15)

此處α2=D/ρh.

考慮橋面板在縱向兩端簡支和橫向兩端自由，令Z(x,y)=X(x)Y(y)，邊界條件如下.

縱向兩端簡支邊界條件：

(16)

(17)

橫向兩端自由邊界條件：

(18)

假設板豎向振型為

(19)

式(19)中δ為大于零的不定常量，可根據實際情況進行調整.式(19)滿足邊界條件等式(18).于是得到如下等式

(20)

此處

(21)

將式(20)代入式(15),得到

(22)

分情況解此微分方程.

(23)

代入邊界條件式(16)和式(17)，則A,B,C和D解出現奇異，不符合要求.

(24)

代入邊界條件等式(16)和等式(17)，欲使A,B,C和D有非零解.則有

(25)

所以

sin (γL)=0

故而

(26)

所以，根據式(24)可以得到固有頻率

(27)

此處γm和θn相應的由等式(21)和(26)給出.

與頻率ωmn相一致的振型函數如下

(28)

假設前述的特征值問題可以由頻率和振型解決，則橋面板豎向撓度可以由下式給出

(29)

此處Zmn(x,y)由式(28)給出.

將式(29)代入式(9),由于模態正交而忽略阻尼項，即令C=0.則模態控制方程為

(30)

式(30)給出的廣義質量和廣義剛度可由下述式子得到

(31)

(32)

此處γm和θn相應值由式(26)和(21)給出.廣義力Pmn(t)可以通過以下積分式得到

(33)

則式(30)可寫成以下形式

(34)

若考慮阻尼影響，即C≠0時，式(30)可寫成以下形式

(35)

3 工程案例

某高速公路橋梁東西雙線結構分幅獨立設置，橋梁上部構造為4孔20m跨徑的先張法預應力混凝土簡支空心板，20°斜交，單幅每跨上部構造由12片梁組成.空心板預制寬99cm、高85cm；橋梁下部構造為雙柱式墩和三肋式臺配以鉆孔灌注樁基礎.

2009年8月22日凌晨，該橋西線幅橋第四跨橋面發生滿載煙花爆竹的貨車劇烈爆炸，造成西線第四跨結構嚴重損壞，經炸后檢測，除外邊梁斷裂外，其余梁體出現嚴重的下撓現象，各梁下撓與開裂程度基本吻合.特別是2,3號梁下撓值達到756mm和658mm.

該橋面板的外觀尺寸為單跨長L=20.0 m、寬H=12.0 m、高h=1.0 m，密度ρ=2 500 kg/m3，楊氏模量E=34.5 GPa，泊松比υ=0.20.由于實際爆炸荷載的復雜性和實載重量的不確定性，按車輛型號和文獻[3]分析結果，擬定炸藥TNT當量為450kg，爆心距橋面1.2m.

3.1 固有頻率與振型

將上述參數代入式(26)和(21)可得到γm和θn，再將γm和θn計算值代入式(27)得到橋面板的固有頻率.其中ω11與δ取值的具體關系見圖1.

圖1 ω11與δ取值關系

再綜合考慮橋面板振型圖y方向兩端滿足自由的邊界條件和x方向兩端滿足簡支的邊界條件，經試算，當δ值大于或等于35時，振型圖y方向y=H端幅值與y=0端幅值相差很小，即兩端變形基本對稱，故最終確定δ取值為35.此時x方向(縱向)上的振型和y方向(橫向)上的振型的頻率ωmn見表1，其中下標m對應x方向第m階振型;n對應y方向第n階振型.振型Zjj(j=1,2,3,4)見圖2.

表1 頻率ωmn Hz

3.2 橋面板豎向撓度

在計算豎向撓度w(x,y,t)時，取q(x,y,t)為線性下降的三角形荷載，其中Pmax和t0分別為文獻[11]中超壓公式和正壓作用時間公式計算所得到的超壓峰值和正壓時間.由式(34)可得

圖2 橋面板振型圖

(36)

(37)

將式(20)和(37)代入式(29)并進行整理得豎向撓度

(38)

為求得跨中位置的豎向位移，應確定各變量的取值.則δ=35，x=10m，y=6m，又當m，n取值超過10以后，w(x,y,t)基本趨于穩定，故計算時m，n均取1～10即可.當炸藥距離橋面為1.2m至9m時，豎向位移w(x,y,t)隨t變化的關系曲線見圖3.

圖3 不同爆心距離時豎向位移時程曲線

由圖3可知，當爆心距橋面1.2m時，基于本文理論計算方法所得的跨中撓度為678mm，與實際梁體最大跨中撓度為756，658mm分別相差11.5%和2.9%，造成上述偏差的原因主要是由于對炸藥當量和爆心位置的判斷存在偏差所致.表明該分析方法對于常規的梁板式混凝土橋梁抗爆分析具有較好的精度.

另外，在炸藥當量不變時，隨著爆心距離增加，爆炸荷載所產生的跨中最大位移將迅速降低，且這種降勢呈先急后緩態勢.當爆心距離由1.2m提升至3m，此時比例距離由0.156m/kg1/3變為0.315m/kg1/3，最大跨中撓度將由678mm降至48mm.由此可見，合理設計比例距離是進行結構抗爆設計的有效途徑.

4 結束語

由于爆炸荷載的不確定性，有關橋梁結構抗爆研究的理論方法相對較少.本文基于動力薄板理論，將混凝土橋梁的橋面板做整體薄板考慮，提出具體的動態控制方程，開展實橋結構動力特性分析，并通過理論分析結果與實測數據的對比分析，驗證了基于薄板理論開展常規梁板式混凝土橋梁抗爆性能研究的有效性與可行性.

同時，由本文分析可知，隨著比例距離增大，爆炸荷載效應會迅速減小，合理控制比例距離是實現橋梁結構抗爆設計的有效途徑.

[1]RONG H ,LI B.Deformation-controlled design of reinforced concrete flexural members subjected to blast loadings [J]. Journal of Structural Engineering,2008,134(10):1598-1610.

[2]FANG J Q,CHUNG P W.WOLFE R.Analysis of a blast-loaded protective wall for bridge columns [J].Bridge Structures, 2008(3):135-141.

[3]WINGET D G,MARCHAND K A,WILLIAMSON E B.Analysis and design of critical bridges subject to blast loads [J].Journal of Structural Engineering,2005,131(8):1243-1255.

[4]LI B,HUANG Z,LIM C.Verification of nondimensional energy spectrum-based blast design for reinforced concrete members through actual blast tests[J].Journal of Structural Engineering,2010,136(6):627-636

[5]WILLIAMSON E B,BAYRAK O,DAVIS C,et al.Performance of bridge columns subjected to blast loads: experimental program [J].Journal of Bridge Engineering,ASCE,2011,16(6):693-702

[6]WILLIAMSON E B,BAYRAK O,DAVIS C，et al.Performance of bridge columns subjected to blast loads: results and recommendations [J].Journal of Bridge Engineering,ASCE,2011，16(6):703-710

[7]FUJIKURA S,MICHEL B,DIEGO L.Experimental investigation of multihazard resistant bridge piers having concrete-filled steel tube under blast loading[J].Journal of Bridge Engineering,ASCE,2008,13(6):586-594

[8]孫文彬．鋼筋混凝土對邊簡支板的爆炸試驗[J]．力學與實踐．2008,30(4):58-60.

[9]畢繼紅，王 繪．工程彈塑性力學[M]．天津：天津大學出版社，2008.

[10]LI Q S,FANG J Q，JEARY A P.Free vibration analysis of multi-storey buildings with narrow rectangular configuration[J].Structure Engineering,1999:135-141.

[11]LIU C, HU C,WANG Y,et al.Numerical simulation of peak overpressure and positive pressure time of explosion shock wave in the air[J].Advanced Materials Research,2012,472(2):1500-1504.

Method of Close Blast Response Analysis for Concrete Bridges Based on Shear-flexural Plate Theory

HU Zhijian1)LIU Chao1,2)FANG Jianqiao1)

(SchoolofTransportation,WuhanUniversityofTechnology,Wuhan430063,China)1)(XiangyangQualitySupervisionBureauofTrafficEngineering,Xiangyang441000 ,China)2)

In this paper the research focuses on close blast loadings. And the analytical method of analyzing the dynamic response of a bridge deck subjected to the blast loadings is presented with the deck be treated as a shear-flexural plate based on the deformation observations and analyses, since beam/slab structures are widely applied for highway bridges with low values for both depth to width ratio and depth to span ratio. Then the governing equation is derived, and the method of the blast response analysis is proposed. For a real bridge project, the numerical results, including dynamic characteristics and structural response are given to verify and validate the proposed shear-flexural plate model. The obtained data have also indicated that scaled distance is one of the key parameters for blast resistant design for structures and should be chosen appropriately.

shear-flexural plate theory; close blast loadings; concrete bridges; scaled distance

2014-12-10

*江西省科技支撐項目(批準號：20132BBG7010)、中央高校基本科研業務費專項資金項目(批準號：2013V0016)資助

U445.75

10.3963/j.issn.2095-3844.2015.02.011

胡志堅(1974- )：男，博士，副教授，主要研究領域為橋梁結構理論與安全技術