耿少波，羅 干，陳佳龍，趙 洲

（中北大學土木工程學科部，山西 太原 030051）

易爆物品儲存運輸不當、燃氣爆炸及暴恐襲擊等時有發(fā)生，建筑結構在其服役期內遭受空爆荷載的概率逐漸增大。在進行結構抗爆分析時，抗爆設計規(guī)范普遍推薦采用無阻尼結構動力學體系，實現(xiàn)基于動力系數的等效靜載抗爆分析。空爆荷載作用時長很短，可近似簡化為三角形衰減荷載，這種簡化可使不熟悉空爆荷載的結構設計人員完成抗爆計算。目前，民用建筑抗爆設計采用延性比等參數完成彈塑性抗爆分析，延性比為結構振動彈塑性變形最大值與彈性變形最大值的比。無阻尼自由振動體系沒有能量耗散，是一種無休止的簡諧振動。忽略阻尼作用，將放大結構振動各個階段的位移幅值，而對延性比、動力系數等防爆設計參數的影響，還沒有明確的理論結論。

認識到阻尼對結構振動響應存在影響，已開展了一些研究。在理論方面：Biggs采用等效單自由度（single degree of freedom, SDOF）進行抗爆分析時，提及阻尼對塑性階段結構振動存在一定影響；Gantes 等分析空爆作用結構彈塑性振動時，指出阻尼對結構響應前幾個振動周期存在一些影響；Riedel 等認為對于空爆作用下結構失效的情況，阻尼可忽略不計；方秦等建立并求解了端部有阻尼支承的梁體空爆作用下的振動方程，表明空爆荷載作用時長越短，端部的阻尼支撐對梁體的抗力提高作用越顯著；郭東等采用杜哈梅積分方法，求解了空爆作用下彈性階段含阻尼體系的等效單自由度動力方程，表明阻尼對反彈階段的彈性位移振動影響顯著，建議按30%進行位移折減，但未解決阻尼參數對塑性階段響應的影響；陳萬祥等求解了含阻尼的柔性邊界支承下淺梁的振動方程，認為邊界阻尼對結構的破壞模式不會發(fā)生變化；董彬等通過數值方法分析了含阻尼體系的梁體振動，表明加設阻尼器能有效控制空爆作用下的動力響應；杜志鵬等將船體結構簡化為梁模型，分析了水下爆炸船體鞭狀運動的阻尼效應，結果表明不考慮阻尼效應將高估運動響應幅值。由空爆作用結構試驗可知，阻尼對沖擊波荷載結束后的自由振動階段確實存在影響：Liu 等分析了尺寸0.15 m×0.15 m×1.7 m 的鋼筋混凝土梁在近場空爆作用下的破壞特征，實測位移顯示，阻尼影響下結構在4～5 個振動周期后靜止；Nassr 等完成了5 組工字型鋼梁在遠場空爆下的動力響應，4～6 個振動周期后結構靜止；Zhang 等完成了50 kg炸藥近場空爆作用下長2.5 m 鋼管混凝土梁構件的振動效應，結果顯示3～7 個振動周期后構件靜止；Liu 等進行了尺寸0.22 m×0.3 m×2 m 的鋼筋-玻璃纖維-混凝土梁構件在0.3～4 kg 當量炸藥下近場空爆試驗，發(fā)現(xiàn)梁體彈性振動、塑性振動、截面斷裂3 種響應類型均在2～4 個振動周期后靜止，塑性振動及截面斷裂對應的阻尼明顯偏大；Nagata 等、Syed 等、Ritchie 等、Shi 等和Foglar 等完成的結構空爆試驗也充分說明阻尼對結構振動的顯著影響。阻尼對空爆作用結構彈性階段和塑性階段強迫振動的影響、影響動力系數的程度、進而對抗爆設計的影響，試驗上無法直接判別，也缺乏理論上的研究探索。

本文中，建立含阻尼的等效自由度體系振動模型，考慮空爆荷載作用時長與結構進入塑性振動時長的關系，將結構分為柔性結構、剛性結構、臨界結構等3 種類型，進行各種類型下結構彈塑性振動方程的求解，結合延性比、動力系數等抗爆設計參數定義，以典型的阻尼比設計驗算工況，對比現(xiàn)行抗爆設計規(guī)范的推薦公式，考查阻尼對動力系數的影響。

1 柔性結構空爆作用動力響應求解

1.1 彈塑性階段動力方程

采用延性比的概念進行梁式及單向板結構空爆作用動力分析，且結構響應以彎曲振動分析為主時，等效單自由度方法具有良好的計算精度、簡易的計算流程，被廣泛采用，如圖1 所示。

圖1 理想彈塑性的含阻尼等效單自由度體系Fig. 1 Elastic-perfectly plastic SDOF vibration system with damping

考慮阻尼參數的彈性階段等效單自由度振動方程為：

式中：為彈性階段等效結構質量，為彈性階段等效結構阻尼，為等效結構剛度，為與等效結構相等的、真實結構在典型位置處的振動位移，Δ()為等效結構承受的隨時間變化的等效空爆荷載。各等效系數分別為：

式中：為結構每延米的質量，為結構跨長，ξ 為結構阻尼比，為真實結構剛度，為彈性階段質量變換系數，為彈性階段荷載變換系數?？毡奢d持續(xù)時間非常短，可簡化為等沖量的線性荷載，GB 50038–2005《人民防空地下室設計規(guī)范》中，推薦采用等效空爆荷載：

式中：為空爆荷載作用時長，為沖擊波超壓峰值，為彈性階段荷載變換系數或塑性階段荷載變換系數。

設結構彈性位移最大值對應的時刻為，對應的振動速度為，此為結構進入塑性振動的區(qū)分點，將彈性參數替換為塑性參數后，結構塑性階段方程為：

式中：為塑性階段等效結構質量，為塑性階段等效結構阻尼，為結構最大抗力。各等效系數分別為：

式中：為塑性階段質量變換系數。

1.2 結構分類及柔性結構彈性階段動力響應求解

柔性結構指的是空爆荷載作用時長小于該結構振動從0 至最大彈性位移的時長，即＜；類似地，剛性結構指的是＞對應的結構，臨界結構指的是=對應的結構。

求解空爆作用結構動力方程時，常采用杜哈梅積分方法，在求解過程中，會出現(xiàn)多次分部積分，計算過程復雜。根據微分方程理論：任意微分方程解答均可表示為通解與特解之和的形式，這個方法可一定程度簡化本文動力方程的求解。

對于柔性結構，在彈性階段且在荷載作用時長范圍0≤≤，由初始位移和初始速度均為0，結合式(1)、(3)可求該強迫振動階段解：

式中：ω 為無阻尼振動等效頻率，ω為含阻尼振動等效頻率，為超壓峰值為靜載時的結構位移。各參數計算公式為：

將代替代入式(6)～(7)，即可得到強迫振動結束時的位移、速度。當空爆荷載消去，這個階段結構外荷載為0、以位移及速度為初始條件的含阻尼彈性自由振動（即當＜≤），可由式(1)、(3)求解：

式中：γ 為阻尼綜合降低系數，θ、θ為結構時長參數。計算公式為：

1.3 柔性結構塑性階段動力響應求解

根據理想彈塑性理論假設，結構完成彈性振動后結構達到最大的抗力，在塑性階段該抗力保持不變，此時為外荷載為0、以位移及速度為初始條件的含阻尼塑性階段自由振動。由式(4)可得：

在達到彈塑性位移最大值時，此時振動速度=0。令式（15）右為0，得出：

將式(16)對應的時刻代入式(14)，得出結構彈塑性振動最大位移為：

代入等效單自由度體系對應參數，即將式(5)中、及代入(17)，可得：

式中：、分別為彈性、塑性階段質量變換系數與荷載變換系數的比。即：

由結構彈塑性理論及抗爆設計規(guī)范，彈塑性階段抗力動力系數和延性比β 分別為：

將式(11)、(18)代入式(20)，且令：

則得到柔性結構延性比β 關于抗力動力系數的表達式：

2 剛性及臨界結構空爆作用動力響應求解

2.1 剛性結構彈性階段動力響應求解

根據定義及理論分析可知，剛性結構彈性階段（0＜≤）與柔性結構0＜≤時段的振動方程相同，將代入式(6)～(7)，可得時刻的位移和速度：

2.2 剛性結構塑性階段動力響應求解

由定義可知，剛性結構從彈性振動進入塑性振動后，第1 個塑性響應階段為外荷載不為0、以位移和速度為初始條件的含阻尼塑性階段強迫振動，即在＜＜內，由式(4)求出其動力響應解答后，可得到時刻結構振動位移和速度分別為：

剛性結構振動＞時，為無外荷載、位移和速度為初始條件的含阻尼塑性階段自由振動，即當＜≤時，該振動方程與柔性結構塑性階段響應方程求解一致，將式(14)～(16)中用替換為后，便得出其解答，其中為：

同理，將式(5)中、和代入式(27)后，可將等效單自由體系轉變?yōu)樵Y構參數。為精簡篇幅、清晰顯示延性比對應的參數關系，本文中省略該化簡過程。令：

由式(20)、(23)和(27)可知，基于動力系數的延性比為：

2.3 臨界狀態(tài)結構動力響應求解

根據臨界狀態(tài)定義可知，空爆荷載作用時長與結構完成彈性振動時長恰好相等，即=。在彈性階段且在荷載作用時長范圍0＜＜，由彈性振動式(23)～(24)可得：

振動時刻大于時，外荷載為0、以和為初始條件的含阻尼塑性階段自由振動，即當＜＜時，與柔性結構塑性階段振動一樣，即利用式(18)、(20)、(23)～(24)，且令：

可得臨界結構基于動力系數的延性比：

由式(22)、(31)、(35)可以看出，將動力系數化解為關于阻尼比、延性比等參數的顯性表達式難度很大，幾乎不可能。因此，采用本文公式進行考慮阻尼的動力系數分析時，可根據結構類型的界定條件迭代計算。

3 算例驗證

3.1 算例工況遴選

GB 50038–2005《人民防空地下室設計規(guī)范》中，推薦采用等效單自由振動動力系數計算方法。無阻尼、不含塑性參數的簡化公式為：

為了校核本文公式的精準性，以式（36）為對比前提，以簡支梁為結構選型，以文獻[1]中推薦的鋼筋混凝土構件允許延性比1～4 為計算范圍，以文獻[23]中推薦的空爆結構-荷載參數θ≤2.2 為參數范圍；為了獨立觀察阻尼對動力系數的影響程度，選用典型的塑性參數，即彈性、塑性階段等效質量系數與等效荷載系數之比、，由文獻[5]取常數值0.78、0.66；在阻尼比為0.0001～0.1 時，選擇5 種典型阻尼比，即0.000 1（接近無阻尼）、0.001（極小阻尼）、0.01（小阻尼，建筑鋼構件可采用數值）、0.05（常用阻尼比，鋼筋混凝土及砌體構件常采用數值）、0.1（較大阻尼，塑性階段可能性數值）作為建筑物構件典型的阻尼比。共計20 種典型計算工況，見表1。

表1 典型工況Table 1 Typical calculation cases

為了查看動力系數與延性比、阻尼比的關系，考慮阻尼比后的20 種工況計算結果，與文獻[1]中公式在延性比為1～4 時計算結果的比較，如圖2 所示；各工況計算結果與文獻[1]中公式計算結果的相對誤差如圖3(a)所示，與阻尼比為0.000 1 工況計算結果的相對誤差如圖3(b)所示。

圖2 本文工況計算結果與文獻公式的比較Fig. 2 Comparison of the results from the calculation cases and from the code formula

圖3 本文工況計算結果的相對誤差Fig. 3 Relative errors of the calculation cases

延性比β=1 表示彈塑性抗爆設計退化為彈性設計，此時各阻尼比的動力系數曲線均為光滑曲線。延性比增加表示結構的塑性變形比例增加，隨著β 從2 遞增到4，動力系數對參數θ的斜率從0.8 至0.6 附近發(fā)生轉折。相同延性比數值下，阻尼比越大，動力系數越小，阻尼比為0.000 1（接近無阻尼）、0.001（極小阻尼）、0.01（小阻尼）的差異非常小。相同阻尼比情況下，隨著延性比β 的增加，相同θ對應的動力系數降低。規(guī)范計算公式與柔性結構、臨界結構計算結果差異很小，與剛性結構動力系數計算結果差異較大。

3.2 結果分析與討論

由圖2 可知，延性比β 為1 時：文獻[1]中公式計算結果與阻尼比為0.0001～0.1 的計算結果差異均較小，低于阻尼比0.0001 的計算結果，最小為0.06%，最大為4.23%，平均為2.43%；文獻[1]中公式計算結果位于阻尼比0.01 與0.05（常用阻尼比）之間，相差最小值為0.07%，最大值為2.63%，平均值為0.86%；阻尼比為0.1 時，差異變得顯著，平均值為11.73%。可見，用文獻[1]中公式進行彈性階段抗爆設計具有很高計算精度和優(yōu)勢，可觀察到增加5%以上的阻尼比，對抗爆設計具有明顯的經濟效益。

由圖3(a)可知：延性比β 為2～4 時，文獻[1]中公式計算結果與柔性結構、臨界結構計算公式的計算結果差異很小，略低于阻尼比0.0001、0.001 的，略高于阻尼比0.01～0.1 的，其差異為0～4%，用文獻[1]中公式進行柔性結構抗爆設計仍具有較高精度；而較大幅度低于本文剛性結構的計算結果，在阻尼比0.0001～0.1 下，延性比β=2 時最大差異為31.88%～14.52%，延性比β=3 時最大差異為38.83%～18.25 %，延性比β=4 時最大差異為41.49%～18.09%。這種差異與文獻[1]中忽略阻尼比及塑性參數有關，也與文獻[1]中擬合公式時選擇的類型有關。

由圖3(b)（以阻尼比0.0001 對應的動力系數為基準）可知：阻尼比為0.001 時，動力系數結果降低有限，其差異為?0.11%～?0.20%，基本可以忽略，說明結構阻尼比在0.001 以下，可按無阻尼進行結構抗爆設計；阻尼比為0.01 時，動力系數值略有降低，其差異為?1.49%～?2.08%，數值較小，在結構設計允許范圍內，可忽略該阻尼比對結構的影響；阻尼比為0.05 時，動力系數值降低明顯，平均值為?7.33%～?9.92%；阻尼比為0.1 時，對應的動力系數降低非常明顯，動力系數值降低平均值為?13.42%～?19.43%；阻尼比大于0.05 時，臨界結構的動力系數差異性明顯低于柔性結構、剛性結構數值，即對于大多數抗爆結構，不應忽略此時的阻尼影響，可增設阻尼構造來降低抗爆設計的工程造價。

4 結 論

推導了空爆作用下柔性結構、剛性結構和臨界結構等效靜載動力系數的隱式函數表達式，通過典型計算工況，分析了阻尼比對動力系數的影響，得到以下結論。

（1） 阻尼比0.001、0.01 較阻尼比0.0001 的動力系數值降低幅度約0.20%、2.08%，即阻尼比小于0.01 時，數值差異性很小，即可忽略阻尼比小于0.01 對抗爆設計動力系數影響作用。

（2） 阻尼比0.05、0.1 較阻尼比0.0001 的動力系數值降低幅度約9.92%、19.43%，說明將抗空爆結構的阻尼比提高至0.05 以上，將產生有較大的經濟效益，是一種良好的抗空爆設計手段。

（3） GB 50038–2005《人民防空地下室設計規(guī)范》中，空爆動力系數計算公式忽略阻尼比、塑性參數，對完成彈性階段抗爆設計具有很高計算精度和優(yōu)勢。當進行彈塑性抗爆設計時，規(guī)范中公式較適用于柔性結構、臨界結構抗爆設計，且誤差在4%以內。運用于剛性結構時，規(guī)范中公式均小于本文中計算公式計算結果，且阻尼比越小誤差越大，延性比β=4 時相對差異值可達18.09%～41.49%。