外爆荷載下K iew itt 8型單層球面網(wǎng)殼的動(dòng)力響應(yīng)

馬加路，支旭東，范 峰，Stewart G.

（1.哈爾濱工業(yè)大學(xué)結(jié)構(gòu)工程災(zāi)變與控制教育部重點(diǎn)實(shí)驗(yàn)室，哈爾濱 150090；2.哈爾濱工業(yè)大學(xué)土木工程學(xué)院，哈爾濱 150090；3.Centre for Infrastructure Performance and Reliability，The University of Newcastle，New South Wales，2308，Australia）

外爆荷載下K iew itt 8型單層球面網(wǎng)殼的動(dòng)力響應(yīng)

馬加路1，2，支旭東1，2，范 峰1，2，Stewart G.3

（1.哈爾濱工業(yè)大學(xué)結(jié)構(gòu)工程災(zāi)變與控制教育部重點(diǎn)實(shí)驗(yàn)室，哈爾濱 150090；2.哈爾濱工業(yè)大學(xué)土木工程學(xué)院，哈爾濱 150090；3.Centre for Infrastructure Performance and Reliability，The University of Newcastle，New South Wales，2308，Australia）

為研究Kiewitt8型單層球面網(wǎng)殼在爆炸荷載下的動(dòng)力響應(yīng)，基于LS-DYNA建立了爆炸中考慮附屬結(jié)構(gòu)（檁托，檁條，鉚釘和屋面板）影響的40 m單層球面網(wǎng)殼結(jié)構(gòu)有限元模型。通過(guò)分析相同爆炸荷載下不同結(jié)構(gòu)矢跨比，構(gòu)件尺寸和材料屬性等參數(shù)對(duì)結(jié)構(gòu)響應(yīng)的影響，得到了球面網(wǎng)殼的矢跨比、桿件截面、屋面板厚度及屋面活荷載等是影響結(jié)構(gòu)動(dòng)力響應(yīng)的主要因素。在此基礎(chǔ)上，通過(guò)討論相同結(jié)構(gòu)在不同爆炸作用下的結(jié)構(gòu)破壞及倒塌情況，提出了單層球面網(wǎng)殼在外爆荷載下安全距離的計(jì)算方法，為更好地建立網(wǎng)殼結(jié)構(gòu)的抗爆設(shè)計(jì)方法和安全評(píng)估理論提供了依據(jù)。

外部爆炸；單層球面網(wǎng)殼；動(dòng)力響應(yīng)；安全距離

美國(guó)俄克拉荷馬州聯(lián)邦大樓爆炸案以及國(guó)內(nèi)外大量的汽車(chē)炸彈襲擊案例，已經(jīng)引起對(duì)于民用建筑如何有效抵御外部爆炸襲擊等問(wèn)題的持續(xù)關(guān)注。文獻(xiàn)［1－3］對(duì)爆炸試驗(yàn)及爆炸作用簡(jiǎn)化計(jì)算理論的大量研究；近年來(lái)，Wu等［4］利用數(shù)值模擬和試驗(yàn)結(jié)合的方法對(duì)爆炸荷載的一些經(jīng)驗(yàn)公式進(jìn)行了改進(jìn)。在此基礎(chǔ)上，國(guó)內(nèi)外學(xué)者對(duì)于受到?jīng)_擊載荷的方形勁肋鋼板［5－6］、鋼結(jié)構(gòu)梁柱、桁架等構(gòu)件及節(jié)點(diǎn)［7－9］的抗爆性能也進(jìn)行了探討。但是，上述多是集中于單獨(dú)構(gòu)件或者尺寸較小的結(jié)構(gòu)，對(duì)于更易受到恐怖襲擊關(guān)注的大型公共建筑或者體育場(chǎng)館中的大跨空間結(jié)構(gòu)抗爆的研究卻十分少見(jiàn)。文獻(xiàn)［10］通過(guò)POD（本征正交分解）方法對(duì)某柱殼結(jié)構(gòu)的爆炸荷載及爆炸響應(yīng)進(jìn)行研究，但受計(jì)算條件所限，研究對(duì)象仍是主體結(jié)構(gòu)，對(duì)于受爆炸影響較大的屋面板，填充墻等附屬結(jié)構(gòu)所承受的爆炸沖擊波的影響卻沒(méi)有合理考慮。

本文對(duì)于工程中典型的跨度40 m Kiewitt8型單層球面網(wǎng)殼建立了考慮屋面板等附屬結(jié)構(gòu)的有限元模型，采用ConWep軟件包，對(duì)不同大當(dāng)量外部爆炸下網(wǎng)殼結(jié)構(gòu)的動(dòng)力響應(yīng)進(jìn)行了分析，討論了結(jié)構(gòu)矢跨比、構(gòu)件截面、材料屬性等參數(shù)對(duì)結(jié)構(gòu)響應(yīng)的影響規(guī)律；并提出單層球面網(wǎng)殼遭受外爆時(shí)安全距離的計(jì)算方法，為該類(lèi)結(jié)構(gòu)的抗爆設(shè)計(jì)及安全評(píng)估提供了基礎(chǔ)數(shù)據(jù)。

1 有限元模型及結(jié)構(gòu)的動(dòng)力響應(yīng)指標(biāo)

1.1 爆炸荷載

對(duì)于爆炸問(wèn)題，傳統(tǒng)的分析方法是建立完整的空氣流場(chǎng)作為沖擊波的傳播媒介，并使用任意拉格朗日歐拉方法（ALE）進(jìn)行流固耦合（FSI）分析［11］。但由于網(wǎng)殼結(jié)構(gòu)跨度大，模擬完整的空氣流場(chǎng)所需要的空氣單元數(shù)量會(huì)導(dǎo)致計(jì)算難以進(jìn)行。針對(duì)這種大尺寸復(fù)雜結(jié)構(gòu)的獨(dú)特性質(zhì)，采用基于Kingery-Bulmash方程建立的ConWep方法對(duì)結(jié)構(gòu)直接施加爆炸荷載更為高效實(shí)用［12］。ConWep能夠自動(dòng)考慮入射角度對(duì)于反射超壓的影響，從而實(shí)現(xiàn)對(duì)結(jié)構(gòu)表面施加爆炸載荷的過(guò)程。其在無(wú)遮攔反射面上的加載壓力表達(dá)式［13］：

式中，Pl為加載壓力，Pi為入射超壓，Pr為反射超壓。LSTC公司在其顯式動(dòng)力學(xué)軟件LS-DYNA中集成了ConWep編碼，并通過(guò)*LOAD＿BLAST等關(guān)鍵字把爆炸效應(yīng)直接加載到目標(biāo)結(jié)構(gòu)上［14］。此方法不需要空氣單元作為傳播介質(zhì)，大大節(jié)約計(jì)算資源。計(jì)算中只需提供等效TNT當(dāng)量，起爆點(diǎn)坐標(biāo)以及爆炸作用方式（地表爆炸或空中爆炸）。

1.2 有限元模型

爆炸荷載作用下，材料的應(yīng)變率可達(dá)到100～10 000／s［5］，所以此時(shí)材料的應(yīng)變率效應(yīng)不可忽略。本文網(wǎng)殼有限元模型材料采用24號(hào)材料（關(guān)鍵字*MAT＿PIECEWISE＿LINEAR＿PLASTICITY［14］），可以根據(jù)用戶需要自定義材料的應(yīng)力應(yīng)變關(guān)系，并能夠考慮應(yīng)變率效應(yīng)的影響。網(wǎng)殼模型中所有構(gòu)件均采用Q235鋼材，其屈服強(qiáng)度為235 MPa。為了對(duì)比彈性模量及泊松比變化對(duì)結(jié)構(gòu)響應(yīng)的影響，其取值在下文的計(jì)算中略有變化。

爆炸荷載首先作用在結(jié)構(gòu)屋面板等附屬結(jié)構(gòu)上，然后通過(guò)連接構(gòu)件把荷載傳遞到結(jié)構(gòu)的主要受力構(gòu)件上。為了更加準(zhǔn)確地模擬這些動(dòng)力行為，需建立較為精細(xì)的網(wǎng)殼結(jié)構(gòu)有限元模型。本文建立的網(wǎng)殼結(jié)構(gòu)計(jì)算模型包含了附屬結(jié)構(gòu)及連接件等，如圖1（a）所示。有限元模型由網(wǎng)殼肋桿、斜桿、檁托、檁條、屋面板及鉚釘?shù)冉M成，各構(gòu)件之間相對(duì)位置如圖1（b）所示，其詳細(xì)尺寸如表1。

圖1 Kiewitt8型單層球面網(wǎng)殼有限元模型Fig.1 FEmodel of the single layer reticulated dome

單層網(wǎng)殼結(jié)構(gòu)中各構(gòu)件通過(guò)焊接與球節(jié)點(diǎn)相連，因此有限元模型中用剛接節(jié)點(diǎn)模擬；計(jì)算過(guò)程中各個(gè)構(gòu)件之間可能產(chǎn)生的相互接觸用自動(dòng)接觸關(guān)鍵字*CONTACT＿AUTOMATIC＿GENERAL［14］處理，其接觸算法為罰函數(shù)算法，其原理是：在每一個(gè)時(shí)間步首先檢查各個(gè)節(jié)點(diǎn)是否穿透主面，如果沒(méi)有穿透，認(rèn)為接觸仍沒(méi)有發(fā)生，不作任何處理。如果穿透，則在該從節(jié)點(diǎn)與被穿透的主面之間引入一個(gè)較大的接觸面力，其大小與穿透的深度、主面的剛度等成正比［14］。從而解決了爆炸過(guò)程中某些構(gòu)件發(fā)生大變形之后產(chǎn)生相互穿透的問(wèn)題。

表1 Kiewitt－8型單層球面網(wǎng)殼構(gòu)件尺寸Tab.1 Details of com ponent dimensions

1.3 結(jié)構(gòu)的動(dòng)力響應(yīng)指標(biāo)

不同于往復(fù)作用的地震荷載中網(wǎng)殼變形特點(diǎn)，爆炸荷載作用下，網(wǎng)殼結(jié)構(gòu)的整體變形及塑性分布情況比較相似，都主要集中在網(wǎng)殼迎爆面一側(cè)。因此，本文僅選擇網(wǎng)殼中節(jié)點(diǎn)最大相對(duì)位移和不同屈服程度桿件比例兩個(gè)物理量作為單層球面網(wǎng)殼動(dòng)力響應(yīng)的衡量標(biāo)準(zhǔn)［15］。

（1）網(wǎng)殼中節(jié)點(diǎn)最大相對(duì)位移。為便于分析，本文定義矢跨比1／5的標(biāo)準(zhǔn)桿件截面網(wǎng)殼結(jié)構(gòu)遭受50 m外等效當(dāng)量1 000 kgTNT炸藥作用的算例中結(jié)構(gòu)的節(jié)點(diǎn)最大位移0.279 m作為節(jié)點(diǎn)位移基準(zhǔn)值Ds，max。節(jié)點(diǎn)最大相對(duì)位移DR就可以表示為任意工況中的節(jié)點(diǎn)最大位移Di，max與基準(zhǔn)值的比值（如式（2）所示），可以從宏觀上體現(xiàn)結(jié)構(gòu)的剛度水平。

（2）不同屈服程度桿件比例。網(wǎng)殼結(jié)構(gòu)有限元模型，每個(gè)Beam單元橫截面上分別有4個(gè)積分點(diǎn)，如圖2 （a），在爆炸過(guò)程中，按照截面的屈服程度不同定義為1P～4P。不同屈服程度桿件比例Ri可以用1P～4P桿件數(shù)量nip占總桿件數(shù)n的百分率進(jìn)行衡量，

式中，i＝1，2，3，4，1P表示4個(gè)積分點(diǎn)中有且僅有1個(gè)進(jìn)入塑性，相應(yīng)的，4P表示全截面進(jìn)入塑性（如圖2（b）所示）。

圖2 網(wǎng)殼桿件截面示意圖Fig.2 The sketches of dome components

2 相同爆炸荷載下的結(jié)構(gòu)參數(shù)分析

假定發(fā)生相同的爆炸時(shí)，分別考慮結(jié)構(gòu)矢跨比、構(gòu)件尺寸以及材料屬性變化對(duì)網(wǎng)殼結(jié)構(gòu)動(dòng)力響應(yīng)的影響。對(duì)40 m跨度的Kiewitt 8型單層球面網(wǎng)殼結(jié)構(gòu)在遭受到比例距離為2 m／kg1／3（距離50 m，TNT等效當(dāng)量1 000 kg）的外部爆炸時(shí)的結(jié)構(gòu)響應(yīng)進(jìn)行了如表2所示的參數(shù)分析。

表2 單層球形網(wǎng)殼結(jié)構(gòu)分析參數(shù)Tab.2 Summary of parametric studies for K iew itt 8 domes

2.1 結(jié)構(gòu)矢跨比的影響

圖3 矢跨比的影響Fig.3 Effect of ratios of rise to span

由圖3可知，網(wǎng)殼結(jié)構(gòu)的節(jié)點(diǎn)最大相對(duì)位移受矢跨比影響很大，差別最大可達(dá)到142%。在結(jié)構(gòu)矢跨比較大時(shí)，節(jié)點(diǎn)最大位移隨網(wǎng)殼矢跨比的減小而下降，在矢跨比為1／5時(shí)節(jié)點(diǎn)相對(duì)位移最??；然后又隨著結(jié)構(gòu)矢跨比減小而增大。為了僅僅對(duì)比矢跨比變化對(duì)結(jié)構(gòu)響應(yīng)的影響，圖3中所涉及的算例是均采用了相同的桿件截面Φ114×4.0 mm。所以，小矢跨比時(shí)桿件截面因低于正常的設(shè)計(jì)標(biāo)準(zhǔn)而出現(xiàn)了位移過(guò)大；而大矢跨比時(shí)雖然桿件相對(duì)較粗，但是同時(shí)也承受更多的爆炸波。從而矢跨比為1／5時(shí)節(jié)點(diǎn)的相對(duì)位移最小。

矢跨比的變化對(duì)于網(wǎng)殼不同屈服程度桿件比例的影響不十分明顯，尤其是1P桿件百分率，結(jié)構(gòu)矢跨比改變對(duì)其幾乎沒(méi)有太大的影響，最大差別也僅有5%；對(duì)于2P、3P和4P百分率影響略微明顯，最大可達(dá)到19%，這說(shuō)明矢跨比的改變并不能明顯改變發(fā)生塑性屈服的桿件數(shù)量。

2.2 結(jié)構(gòu)構(gòu)件尺寸的影響

桿件截面的變化對(duì)于網(wǎng)殼結(jié)構(gòu)的節(jié)點(diǎn)最大相對(duì)位移和不同屈服程度桿件比例都有很顯著的影響，最大差別分別為74%和32%（圖4）。桿件截面越大，結(jié)構(gòu)的整體剛度也隨之增強(qiáng)，網(wǎng)殼結(jié)構(gòu)在爆炸中的響應(yīng)越小。

如圖5所示，屋面板厚度的變化對(duì)節(jié)點(diǎn)最大相對(duì)位移以及不同屈服程度桿件比例也有較明顯的影響。隨著屋面板厚度的增加結(jié)構(gòu)在爆炸荷載下的響應(yīng)逐漸減弱。節(jié)點(diǎn)最大相對(duì)位移差別可達(dá)79%。

本文有限元模型中，屋面活荷載改變的影響是通過(guò)調(diào)整屋面板的密度的方法實(shí)現(xiàn)。如圖6所示，屋面活荷載的大小對(duì)節(jié)點(diǎn)最大相對(duì)位移以及不同屈服程度桿件比例也有很顯著的影響。隨著屋面活荷載的增加結(jié)構(gòu)響應(yīng)迅速減小，節(jié)點(diǎn)最大位移的差別可達(dá)到69%。

圖4 桿件截面的影響Fig.4 Effect of component sections

圖5和圖6中，網(wǎng)殼屋面板厚度及屋面活荷載不斷增加，爆炸荷載作用下網(wǎng)殼的動(dòng)力響應(yīng)呈現(xiàn)出與常規(guī)動(dòng)力分析（例如地震響應(yīng)分析）中截然相反的規(guī)律。其原因是，有限元模型中屋面板厚度和屋面活荷載兩個(gè)參數(shù)的變化的實(shí)質(zhì)都是由屋面板質(zhì)量改變控制的。所以，在爆炸過(guò)程中，屋面板質(zhì)量增加，致使發(fā)生相同的變形所需要的能量增多，從而導(dǎo)致節(jié)點(diǎn)位移及不同屈服程度桿件比例隨之減小。

圖6 屋面活荷載的影響Fig.6 Effect of live loads

2.3 材料屬性的影響

如圖7所示，隨著鋼材彈性模量增加，節(jié)點(diǎn)最大相對(duì)位移變化較小，幾乎可以忽略。彈性模量變化對(duì)網(wǎng)殼結(jié)構(gòu)不同屈服程度桿件比例的影響也很小，1P～4P桿件的數(shù)目幾乎不變。

圖7 彈性模量的影響Fig.7 Effect of Elastic modulus

同樣，在圖8中泊松比的改變對(duì)節(jié)點(diǎn)最大相對(duì)位移以及不同屈服程度桿件比例的影響也都很小，均不超過(guò)5%。因此，爆炸荷載下泊松比變化對(duì)于結(jié)構(gòu)響應(yīng)的影響也可以忽略。

圖8 泊松比的影響Fig.8 Effect of Poisson's ratios

3 網(wǎng)殼結(jié)構(gòu)抗爆安全距離的確定方法

網(wǎng)殼結(jié)構(gòu)具有良好的整體工作性能，部分桿件的破壞或結(jié)構(gòu)局部的大變形未必能引發(fā)結(jié)構(gòu)的整體倒塌。但是，從結(jié)構(gòu)保持穩(wěn)定振動(dòng)（具有一定的損傷）到發(fā)生整體倒塌卻是一個(gè)突然的瞬變過(guò)程。對(duì)于網(wǎng)殼一類(lèi)相對(duì)重要的結(jié)構(gòu)來(lái)說(shuō)，研究如何避免結(jié)構(gòu)發(fā)生整體倒塌顯得十分必要。圖9是網(wǎng)殼結(jié)構(gòu)在相同炸藥當(dāng)量不同爆炸距離時(shí)的結(jié)構(gòu)破壞情況。本文將裝藥量一定時(shí)能夠使網(wǎng)殼結(jié)構(gòu)發(fā)生倒塌的臨界狀態(tài)所對(duì)應(yīng)的爆距定義為網(wǎng)殼結(jié)構(gòu)的抗爆安全距離。

圖9 節(jié)點(diǎn)位移云圖Fig.9 Contours of nodal displacements

圖10 不同炸藥當(dāng)量時(shí)1 s與1.5 s位移對(duì)比Fig.10 Results comparison between 1 s and 1.5 s

在網(wǎng)殼結(jié)構(gòu)爆炸響應(yīng)分析中，由于部分桿件發(fā)生屈服或斷裂，導(dǎo)致材料非線性和整體結(jié)構(gòu)變形持續(xù)加劇，數(shù)值計(jì)算的效率會(huì)受到相當(dāng)嚴(yán)重的影響；同時(shí)，為了準(zhǔn)確計(jì)算毫秒級(jí)的爆炸作用過(guò)程中的結(jié)構(gòu)響應(yīng)，又需設(shè)置較為精細(xì)的計(jì)算參數(shù)，造成后續(xù)結(jié)構(gòu)自由振動(dòng)過(guò)程的分析也相當(dāng)漫長(zhǎng)。因此，本文提出一個(gè)簡(jiǎn)化的判斷結(jié)構(gòu)是否發(fā)生失穩(wěn)倒塌的衡量標(biāo)準(zhǔn)：將外部爆炸發(fā)生之后1 s與1.5 s時(shí)刻節(jié)點(diǎn)最大位移進(jìn)行對(duì)比，若二者的節(jié)點(diǎn)最大位移沒(méi)有明顯區(qū)別，則可認(rèn)為結(jié)構(gòu)不會(huì)發(fā)生倒塌并能夠一直保持現(xiàn)有的穩(wěn)定狀態(tài)；反之，若1.5 s時(shí)刻的節(jié)點(diǎn)最大位移明顯增大，則可認(rèn)為結(jié)構(gòu)發(fā)生整體失穩(wěn)倒塌。對(duì)于網(wǎng)殼結(jié)構(gòu)遭受爆炸荷載來(lái)說(shuō)，1 s時(shí)間已經(jīng)足夠使得結(jié)構(gòu)的變形及材料塑性充分發(fā)展，因此，如果結(jié)構(gòu)在1 s后的變形持續(xù)增大，則是處于振動(dòng)不穩(wěn)定的狀態(tài)。

對(duì)上節(jié)中矢跨比1／5跨度40 m的單層球面網(wǎng)殼（桿件截面Φ114×4.0 mm）進(jìn)行不同爆距（1～20 m）不同當(dāng)量（900～1 200 kg）的爆炸荷載全過(guò)程分析。并將1 s與1.5 s時(shí)刻的計(jì)算結(jié)果同時(shí)繪制于圖10中，則各圖中分叉點(diǎn)所對(duì)應(yīng)的距離就可以認(rèn)為是某一當(dāng)量下該網(wǎng)殼結(jié)構(gòu)抗爆的安全距離。

由圖10可知，在四種不同的炸藥當(dāng)量下，網(wǎng)殼結(jié)構(gòu)在1 s與1.5 s時(shí)刻的節(jié)點(diǎn)最大位移對(duì)比均出現(xiàn)了分叉點(diǎn)，則可以確定此時(shí)所對(duì)應(yīng)的爆距就是該爆炸荷載下網(wǎng)殼結(jié)構(gòu)的安全距離。但是由于恐怖襲擊中炸藥當(dāng)量的不確定性，而通常認(rèn)為，恐怖襲擊中汽車(chē)炸彈的TNT炸藥當(dāng)量最大約為1 000 kg，因此可確定該網(wǎng)殼結(jié)構(gòu)的抗爆安全距離是12 m。

抗爆安全距離可以為擬建建筑抗爆性能設(shè)計(jì)提供重要依據(jù)，而且還可作為已建建筑加固或防御措施的評(píng)估依據(jù)。例如，對(duì)于本文網(wǎng)殼結(jié)構(gòu)應(yīng)在12 m外設(shè)置具有阻隔功能的建筑構(gòu)件（防撞柱、防爆墻或者防爆柵欄等［16］），從而有效地控制外部爆炸對(duì)結(jié)構(gòu)的危害，避免結(jié)構(gòu)發(fā)生整體倒塌。

4 結(jié) 論

本文建立了爆炸荷載下考慮附屬結(jié)構(gòu)影響的Kiewitt8型單層球面網(wǎng)殼精細(xì)化有限元模型，通過(guò)LS-DYNA對(duì)其在外部爆炸荷載下的動(dòng)力響應(yīng)進(jìn)行研究。并得到以下結(jié)論：

（1）針對(duì)大跨空間結(jié)構(gòu)的特點(diǎn)，從結(jié)構(gòu)的全局響應(yīng)出發(fā)，提出爆炸荷載下衡量網(wǎng)殼結(jié)構(gòu)整體動(dòng)力響應(yīng)的指標(biāo)，即節(jié)點(diǎn)最大相對(duì)位移和不同屈服程度桿件比例程度。

（2）對(duì)比相同爆炸荷載下單層球面網(wǎng)殼結(jié)構(gòu)不同的結(jié)構(gòu)矢跨比、構(gòu)件尺寸、材料屬性等參數(shù)對(duì)于網(wǎng)殼節(jié)點(diǎn)最大相對(duì)位移及不同屈服程度桿件比例程度的影響。結(jié)果表明結(jié)構(gòu)矢跨比變化對(duì)于網(wǎng)殼動(dòng)力響應(yīng)的影響最為明顯；其次是桿件截面，屋面板厚度和屋面活荷載的影響；彈性模量，泊松比等材料屬性的影響最小，幾乎可以忽略不計(jì)。

（3）通過(guò)對(duì)特定的網(wǎng)殼結(jié)構(gòu)在不同爆炸荷載時(shí)的動(dòng)力響應(yīng)對(duì)比，提出網(wǎng)殼結(jié)構(gòu)抗爆安全距離計(jì)算方法。并明確了跨度40 m的單層球面網(wǎng)殼應(yīng)至少在距離結(jié)構(gòu)12 m處設(shè)置防爆墻或者防爆柵欄以有效的抵御外部爆炸對(duì)結(jié)構(gòu)的影響。

［1］Baker W E，Explosions in Air［M］.Austin：University of Texas Press，1973.

［2］Henrych J.The Dynamics of Explosion and Its Use［M］.Amsterdam：Scientific Pub.Co.，1979.

［3］UFC 3－340－02.Structures to Resist the Effects of Accidental Explosions［S］.Departments of the Army，the Navy，and the Air Force，2008.

［4］Wu C，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.

［5］Yuen S C K，Nurick G N.Nurick，Experimental and numerical studies on the response of quadrangular stiffened plates.Part I：subjected to uniform blast load［J］.International Journalof Impact Engineering，2005，31（1）：55 －83.

［6］Langdon G S，Yuen SC K，Nurick G N.Experimental and numerical studies on the response of quadrangular stiffened plates.Part II：localised blast loading［J］.International Journal of Impact Engineering，2005，31（1）：85－111.

［7］Guzas E L，Earls C J.Simulating blasteffects on steel beamcolumn members：Methods［J］.Computers＆Structures，2011，89（23／24）：2133－2148.

［8］Guzas E L，Earls C J.Simulating blasteffects on steel beamcolumn members：Applications［J］.Computers＆Structures，2011，89（23／24）：2149－2161.

［9］Sabuwala T，Linzell D，Krauthammer T.Finite element analysis of steel beam to column connections subjected to blast loads［J］.International Journal of Impact Engineering，2005，31（7）：861－876.

［10］高軒能，王書(shū)鵬.大空間柱殼結(jié)構(gòu)爆炸動(dòng)力響應(yīng)Ritz-POD數(shù)值模擬［J］.土木建筑與環(huán)境工程，2010.32（2）：64－70.

GAO Xuan-neng，WANG Shu-peng.Numerical simulation for dynamic response of large-space cylindrical reticulated shell under internal explosion by ritz-POD method［J］.Journal of Civil，Architectural＆Environment Engineering，2010，32 （2）：64－70.

［11］Remennikov A M.A review ofmethods for predicting bomb blast effects on buildings［J］.Journal of Battlefield Technology，2003，6（3）：5－10.

［12］Hyde D.User's guide for microcomputer program conwep，application of TM5－855－1，fundamental of protective design for conventional weapons［R］.Instruction Report SL-88-1，Structures Laboratory，US Army Waterways Experiment Station，Vicksburg，Mississippi，1992.

［13］Randers Pehrson G，Bannister K A.Airblast Loading Model for DYNA2D and DYNA3D［R］.1997.

［14］LSTC.LS-DYNA Keyword User'sManual，Version971［M］.2007.

［15］Zhi xu-dong，F(xiàn)an feng，Shen shi-zhao.Elasto-plastic instability of single-layer reticulated shells under dynamic actions［J］.Thin-Walled Structures，2010，48（10／11）：837－845.

［16］Rose T A，Smith PD，Mays G C.The effectiveness of walls designed for the protection of structures against airblast from high explosives［J］.Structures＆Buildings，1995，110（1）：78－85.

Dynam ic responses of K iew itt 8 single-Layer reticulated domes subjected to outside exp losion loads

MA Jia-lu1，2，ZHIXu-dong1，2，F(xiàn)AN Feng1，2，Stewart G3

（1.MOE Key Laboratory of Structural Engineering Disaster and Control，Harbin Institute of Technology，Harbin 150090，China；2.School of Civil Engineering，Harbin Institute of Technology，Harbin 150090，China；3.Centre for Infrastructure Performance and Reliability，The University of Newcastle，New South Wales 2308，Australia）

In order to study the design methods of single-layer reticulated domes，the FE model for a Kiewitt 8 dome of 40 m span under explosion loads was built with the explicit FE packages LS-DYNA to investigate its dynamic responses considering the effects of non-structural components including purlins，purlin props，rivets and roof panels.The blast loading was generated with ConWep，and applied on the dome directly.The effects of different structural ratios of rise to span，component sizes and material properties，etc parameters on the dome's dynamic responses were also analyzed.The results showed that the ratios of rise to span，sizes of cross sections，thicknesses of roof panels and live loads are themain factors to affect the dome's dynamic responses.At last，the concept and calculation method of safety distance for single-layer reticulated domes subjected to outside explosion were proposed.

outside explosion；single-layer reticulated dome；dynamic response；safety distance

TU393；TU393.3

A

10.13465／j.cnki.jvs.2015.21.012

國(guó)家自然科學(xué)基金項(xiàng)目（51478144，51378151）

2014－06－24 修改稿收到日期：2014－10－17

馬加路男，博士生，1984年生

支旭東男，博士，教授，1977年生