王 楓， 胡 豐， 黃 鵬， 顧 明

(同濟(jì)大學(xué) 土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海 200092)

近年來(lái)建筑主體及其圍護(hù)結(jié)構(gòu)的抗風(fēng)研究受到普遍關(guān)注.當(dāng)建筑結(jié)構(gòu)設(shè)計(jì)需考慮圍護(hù)結(jié)構(gòu)的易損性問(wèn)題時(shí)，風(fēng)致飛射物造成的破壞就成為一個(gè)至關(guān)重要的因素.風(fēng)致飛射物是指在極端風(fēng)氣候條件下(如臺(tái)風(fēng))被吹起的高速飛行的物體[1]，而我國(guó)是世界上熱帶氣旋登陸最多、受其影響最嚴(yán)重的國(guó)家之一[2].風(fēng)致飛射物對(duì)高層建筑立面及低矮房屋屋面造成的沖擊破壞作用[3]，是臺(tái)風(fēng)給登陸地區(qū)帶來(lái)巨大損失的主要原因之一.

目前，國(guó)內(nèi)外學(xué)者對(duì)沖擊效應(yīng)的研究多集中于建筑玻璃、壓型鋼板及復(fù)合板的數(shù)值模擬和沖擊試驗(yàn)[4-6]，風(fēng)致飛擲物對(duì)低矮建筑屋面瓦片的沖擊破壞未得到關(guān)注.此外，在已有風(fēng)致飛射物損傷預(yù)測(cè)的統(tǒng)計(jì)模型中，大多數(shù)是針對(duì)“區(qū)域”尺度(如城市)進(jìn)行評(píng)估[7]，很少有研究關(guān)注“單體”尺度(如房屋)的損傷評(píng)估.而多數(shù)風(fēng)災(zāi)破壞往往發(fā)生于低矮房屋，對(duì)低矮房屋單體風(fēng)災(zāi)破壞的評(píng)估與預(yù)防是亟待解決的問(wèn)題.為有效預(yù)防風(fēng)致飛射物引發(fā)的次生災(zāi)害，有必要對(duì)其造成的沖擊破壞機(jī)理進(jìn)行深入研究.

多項(xiàng)調(diào)查顯示[8-9]，風(fēng)致飛射物對(duì)低矮房屋屋面瓦片破壞作用顯著.選取一種典型的陶土屋面瓦片，參考落球沖擊試驗(yàn)數(shù)據(jù)[10]，基于ANSYS/LS-DYNA平臺(tái)，深入研究球狀飛射物對(duì)瓦片的沖擊碰撞機(jī)理，建立相應(yīng)損傷評(píng)估方程，對(duì)瓦片的抗沖擊性能做相應(yīng)分析，并計(jì)算陶土瓦片在球狀混凝土飛射物沖擊下的風(fēng)速閾值.

1 落球沖擊試驗(yàn)

落球沖擊試驗(yàn)是測(cè)試物體抗沖擊性能最直接有效的方法[11].筆者參與文獻(xiàn)[10]中對(duì)陶土瓦片的落球沖擊試驗(yàn)，該試驗(yàn)選取一個(gè)質(zhì)量為0.549 kg的混凝土球塊對(duì)3組陶土瓦片進(jìn)行落球沖擊試驗(yàn)，瓦片幾何尺寸如圖1所示.試驗(yàn)中，瓦片放置于平臺(tái)上，落球的高度有規(guī)律地增加，直至瓦片被沖擊破壞.試驗(yàn)中獲取的平均落球高度反映了該陶土瓦片的破壞強(qiáng)度.由于試驗(yàn)中落球的質(zhì)量一定，故瓦片的破壞能量與落球的高度成正比.圖2為整個(gè)試驗(yàn)過(guò)程，試驗(yàn)結(jié)果如表1所示.

a 瓦片截面

b 瓦片俯視圖

a 落球b 試驗(yàn)裝置c 黏土瓦片破壞情況

圖2落球沖擊試驗(yàn)

Fig.2Dropballimpacttest

瓦片破壞時(shí)裂紋縱向開(kāi)展，呈現(xiàn)明顯的脆性性質(zhì).

表1 陶土瓦片落球沖擊試驗(yàn)結(jié)果

注：表中“√”表示瓦片完好，“×”表示瓦片破壞.

2 數(shù)值模擬

2.1 數(shù)值模型的建立

考慮到碰撞機(jī)理的復(fù)雜性，為更深入地研究球狀飛射物對(duì)陶土瓦片的沖擊作用，基于ANSYS/LS-DYNA平臺(tái)，對(duì)上述落球沖擊試驗(yàn)建立相應(yīng)的有限元模型.模型中采用3D SOLID164實(shí)體單元建立球體及瓦片模型，將球體與瓦片分為2個(gè)部分，落球撞擊位置為瓦片拱形區(qū)域中心，為使計(jì)算更加精確，對(duì)撞擊區(qū)域附近瓦片網(wǎng)格進(jìn)行加密處理，經(jīng)過(guò)網(wǎng)格無(wú)關(guān)性檢驗(yàn)，最終確定計(jì)算模型如圖3所示.模型共18 032個(gè)單元，其中球體單元864個(gè)，瓦片單元17 168個(gè).約束瓦片底部的平動(dòng)自由度即邊界條件為簡(jiǎn)支.

圖3 落球沖擊試驗(yàn)有限元模型

由于僅關(guān)注沖擊過(guò)程，故設(shè)定初始狀態(tài)為球即將碰撞瓦片的時(shí)刻，落球與瓦片的初始距離為0.01 m，定義接觸類型為面面接觸(surface to surface)，接觸過(guò)程為侵蝕接觸(eroding contact).采用Lagrange與單點(diǎn)積分算法并控制沙漏現(xiàn)象[12]，計(jì)算單位為：cm·g·μs.

2.2 材料模型的選擇

數(shù)值模型中材料本構(gòu)模型的選取對(duì)于計(jì)算結(jié)果的準(zhǔn)確性有決定性的影響[13].Johnson-Holmquist-Concrete模型能較好模擬混凝土球的特性；由于瓦片在試驗(yàn)中體現(xiàn)出脆性性質(zhì)，而Johnson-Holmquist-Ceramic模型在模擬脆性材料方面具有獨(dú)特的優(yōu)勢(shì)，故選用該模型模擬瓦片.

混凝土落球模型參數(shù)的選取參考文獻(xiàn)[14]，由于較少研究資料介紹陶土瓦片的材料參數(shù)，因此設(shè)計(jì)試驗(yàn)獲取了瓦片的相關(guān)材料參數(shù).彈性模量是有限元分析中一個(gè)重要的材料參數(shù)，其定義如式(1)所示：

(1)

式中：E為彈性模量；σ為試件應(yīng)力；P為軸向壓力；A為試件橫截面積；ε為試件軸向應(yīng)變.

根據(jù)美國(guó)材料與實(shí)驗(yàn)協(xié)會(huì)(ASTM)E111-04標(biāo)準(zhǔn)進(jìn)行試驗(yàn)[15]，裝置如圖4a所示，2組試件的橫截面尺寸為0.05 m×0.01 m、高為0.08 m.試驗(yàn)數(shù)據(jù)如圖4b所示，取平均值得到陶土瓦片的彈性模量為1.40×1010Pa.總結(jié)現(xiàn)有文獻(xiàn)[8-11]數(shù)據(jù)及上述試驗(yàn)數(shù)據(jù)，確定混凝土落球及陶土瓦片的材料參數(shù)如表2所示,表中D1、D2為損傷指數(shù).

2.3 數(shù)值計(jì)算結(jié)果

根據(jù)以上材料參數(shù)構(gòu)建的本構(gòu)模型即可進(jìn)行相應(yīng)數(shù)值計(jì)算，本文計(jì)算時(shí)長(zhǎng)為104 μs，經(jīng)過(guò)反復(fù)試算發(fā)現(xiàn)，當(dāng)沖擊速率設(shè)置為u=2.17 m·s-1時(shí)，瓦片僅拱形區(qū)域下側(cè)及水平連接處上端部分開(kāi)裂，并未發(fā)生整體完全破壞，開(kāi)裂情況如圖5a所示.

a 試驗(yàn)裝置

b 試驗(yàn)數(shù)據(jù)

物體密度/(kg·m-3)彈性模量/GPa泊松比剪切模量/GPaD1D2抗拉強(qiáng)度/GPa混凝土落球2 440300.214.860.0401.00.004 0陶土瓦片1 380140.21.960.0010.50.002 6

當(dāng)沖擊速率增大至u=2.30 m·s-1時(shí)裂紋逐漸開(kāi)展至邊緣直至瓦片完全破壞，破壞情況如圖5b所示，瓦片裂紋呈十字對(duì)稱狀，縱向裂紋形狀與落球沖擊試驗(yàn)結(jié)果吻合良好.計(jì)算結(jié)果表明，數(shù)值試驗(yàn)與落球試驗(yàn)中臨界沖擊破壞速率(2.17 m·s-1)誤差為6.0%，這有效證明了數(shù)值模型的正確性.

a 瓦片部分開(kāi)裂b 瓦片完全開(kāi)裂

圖5瓦片損傷破壞情況模擬

Fig.5Tiledamagesimulation

當(dāng)沖擊速率u=2.30 m·s-1時(shí)，進(jìn)一步分析瓦片在碰撞歷程中的von mises應(yīng)力狀態(tài)如圖6所示，落球與瓦片在t=850 μs時(shí)開(kāi)始接觸，瓦片應(yīng)力集中于碰撞點(diǎn)附近區(qū)域，隨著沖擊過(guò)程的進(jìn)行，瓦片應(yīng)力逐漸開(kāi)展，t=1 100 μs時(shí)瓦片內(nèi)部出現(xiàn)最大等效應(yīng)力31.31 MPa，最終在t=1 650 μs時(shí)擴(kuò)散到整體，如圖6a所示.t=1 900 μs時(shí)瓦片邊緣開(kāi)始開(kāi)裂，t=6 650 μs時(shí)裂紋完全開(kāi)展至邊緣，其內(nèi)部應(yīng)力分別如圖6b和6c所示.在整個(gè)時(shí)間歷程中，瓦片最大應(yīng)力出現(xiàn)在13 129號(hào)單元，選取其時(shí)間歷程示于圖6d中，由圖可知，自落球與瓦片接觸時(shí)刻開(kāi)始，該單元內(nèi)部應(yīng)力激增，到達(dá)峰值后迅速衰減，此時(shí)單元已達(dá)到破壞條件而失效.

2.4 參數(shù)分析

以上計(jì)算證明了數(shù)值方法的正確性，在此基礎(chǔ)上對(duì)沖擊速率u、球塊質(zhì)量m、瓦片傾角α這幾個(gè)主要影響參數(shù)進(jìn)行分析，以避免試驗(yàn)操作的復(fù)雜性，更方便快捷地得到可靠的結(jié)果.

以u(píng)=2.30 m·s-1、m=0.549 kg、α=0°為基準(zhǔn)工況，假定各變量之間相互獨(dú)立，以控制變量法分別改變每個(gè)變量的取值，計(jì)算瓦片在各沖擊過(guò)程中的最大等效應(yīng)力.例如在研究沖擊速率u時(shí)，固定球塊質(zhì)量m=0.549 kg、瓦片傾角α=0°，僅改變沖擊速率的取值.具體變量參數(shù)選擇見(jiàn)表3.

a 應(yīng)力完全開(kāi)展時(shí)(t=1 650 μs)

b 開(kāi)始開(kāi)裂時(shí)(t=1 900 μs)

c 完全破壞時(shí)(t=6 650 μs)

d 第13 129號(hào)單元應(yīng)力時(shí)程

圖6u=2.30m·s-1時(shí)沖擊歷程中瓦片應(yīng)力狀態(tài)

Fig.6Tilestressstateduringimpactprocessatu=2.30m·s-1

表3 參數(shù)分析

注：*為基準(zhǔn)工況參數(shù)取值.

對(duì)以上21個(gè)參數(shù)對(duì)應(yīng)的19個(gè)工況(1個(gè)基準(zhǔn)工況加上其余18組參數(shù)對(duì)應(yīng)的18個(gè)工況)分別進(jìn)行計(jì)算，提取每種工況下瓦片的最大等效應(yīng)力，計(jì)算結(jié)果見(jiàn)圖7.由圖可知，瓦片最大等效應(yīng)力與瓦片傾角α相關(guān)性不大，與落球質(zhì)量m和落球速率u近似成正比.當(dāng)u=3.00 m·s-1時(shí)，瓦片最大等效應(yīng)力為35.77 MPa；當(dāng)m=0.560 kg時(shí)，瓦片最大等效應(yīng)力為31.38 MPa；當(dāng)α=0°時(shí)，瓦片最大等效應(yīng)力為31.31 MPa.

3 損傷方程

3.1 飛射物形成條件

在進(jìn)行損傷估計(jì)之前首先應(yīng)確定建筑周邊環(huán)境中的碎片是否具備成為飛射物的條件，根據(jù)Wills等[16]的研究，物體被風(fēng)吹起成為飛射物的臨界條件為

a 沖擊速率

b 球塊質(zhì)量

c 瓦片傾角

(2)

根據(jù)上式可得球狀混凝土飛射物直徑D與來(lái)流風(fēng)速U的關(guān)系如圖8所示，該圖反映了一定風(fēng)速U下能被吹起的最大飛射物直徑D.

圖8 球狀混凝土飛射物最大起飛直徑與來(lái)流風(fēng)速的關(guān)系

以上海地區(qū)為例，參考我國(guó)建筑結(jié)構(gòu)荷載規(guī)范(GB50009—2012)[18]，上海地區(qū)10年、50年和100年重現(xiàn)期空曠地貌10 m高度處10 min平均基本風(fēng)壓分別為0.40 kN·m-2、0.55 kN·m-2和0.60 kN·m-2，換算成基本風(fēng)速分別為24.90 m·s-1、29.20 m·s-1和30.50 m·s-1.假設(shè)該風(fēng)速即為來(lái)流風(fēng)速，對(duì)應(yīng)的球狀混凝土飛射物的最大起飛直徑分別為10 mm、14 mm和16 mm(對(duì)應(yīng)最大起飛質(zhì)量分別為0.001 3 kg、0.003 5 kg和0.005 2 kg).

3.2 能量損傷方程

第2節(jié)的研究表明，當(dāng)沖擊速率較小、球塊質(zhì)量較大時(shí)，瓦片內(nèi)部最大等效應(yīng)力與二者均近似成正比，因而多數(shù)研究采用落球試驗(yàn)得到的臨界動(dòng)量作為材料抗沖擊性能的判別標(biāo)準(zhǔn).試驗(yàn)所得的臨界動(dòng)量值Rp=m·v=0.549 kg×2.17 m·s-1=1.19 kg·m·s-1.進(jìn)一步研究表明，當(dāng)沖擊速率較大、球塊質(zhì)量較小時(shí)，材料能承受的最大沖擊動(dòng)量較落球試驗(yàn)(相對(duì)低速情況)不同.例如選取m=0.070 kg、v=15.00 m·s-1的工況進(jìn)行計(jì)算，其沖擊動(dòng)量為Pdeb=m·v=0.070 kg×15.00 m·s-1=1.05 kg·m·s-1

(3)

式中：Ddeb為損傷指數(shù)；K為飛射物速率與來(lái)流風(fēng)速比，根據(jù)宋芳芳[9]和Holmes[17]的研究成果，球狀飛射物在風(fēng)速小于40 m·s-1時(shí)其最終落地速率約為風(fēng)速的一半，即K=0.5；RD為臨界能量損傷指數(shù).

a 不同速率下球狀飛射物動(dòng)能

b 不同質(zhì)量下球狀飛射物動(dòng)能

Fig.9Kineticenergytimehistoriesofsphericaldebrisduringimpactprocess

由圖9a中v=2.30 m·s-1的標(biāo)準(zhǔn)工況計(jì)算結(jié)果可知，球塊的沖擊動(dòng)能為1.45 J，將此動(dòng)能作為使得瓦片破壞的臨界能量值，即取RD= 1.45 J.

將式(2)代入式(3)，消去飛射物特征長(zhǎng)度，得

(4)

圖10給出了能量損傷指數(shù)與風(fēng)速之間的關(guān)系，由該圖可知能量判別準(zhǔn)則下陶土瓦片能承受不超過(guò)34.35 m·s-1的來(lái)流風(fēng)速下球狀混凝土飛射物的沖擊破壞作用，其含義為34.35 m·s-1的風(fēng)速能吹起最大直徑為20 mm的球狀混凝土飛射物(圖8)，這些飛射物造成的最大沖擊效應(yīng)剛好使得陶土瓦片破壞.

由式(4)可知，能量損傷指數(shù)與風(fēng)速的八次方相關(guān)，這表明能量損傷指數(shù)對(duì)風(fēng)速極為敏感.如圖10所示，當(dāng)風(fēng)速U增大至35 m·s-1以上時(shí)，Ddeb迅速增大，說(shuō)明風(fēng)速越大，其破壞程度顯著增大.

圖10 能量損傷指數(shù)與風(fēng)速的關(guān)系

仍以上海地區(qū)為例，規(guī)范[18]給出上海地區(qū)空曠地貌10 m高度處的陣風(fēng)系數(shù)βgz=1.70，則10年、50年和100年的陣風(fēng)風(fēng)壓(βgzw0，w0為基本風(fēng)壓)分別為0.68 kN·m-2、0.94 kN·m-2和1.02 kN·m-2，換算成陣風(fēng)風(fēng)速分別為32.47 m·s-1、38.18 m·s-1和39.77 m·s-1.因此若該陶土瓦片(能承受風(fēng)速不超過(guò)34.35 m·s-1時(shí)球狀混凝土飛射物的沖擊作用)應(yīng)用于上海地區(qū)，在10年重現(xiàn)期的陣風(fēng)風(fēng)速下是安全的，而在50年和100年重現(xiàn)期陣風(fēng)風(fēng)速下將受到飛射物的沖擊破壞.

3.3 易損性分析

無(wú)論是從動(dòng)量角度還是從能量角度判別材料的抗沖擊性能，沖擊物的質(zhì)量和速度都是至關(guān)重要的2個(gè)物理量.由這2個(gè)物理量可繪制相應(yīng)的易損性包絡(luò)曲線.

已有文獻(xiàn)表明，臺(tái)風(fēng)最大陣風(fēng)風(fēng)速常分布在20到40 m·s-1，對(duì)應(yīng)的球狀飛射物沖擊速率約為10到20 m·s-1.取沖擊速率依次為10、12、14、16、18和20 m·s-1.基于第2節(jié)的數(shù)值模型，取不同質(zhì)量的球狀飛射物反復(fù)進(jìn)行數(shù)值模擬試驗(yàn)，記錄瓦片的沖擊破壞情況.試驗(yàn)結(jié)果見(jiàn)表4.

根據(jù)上表的臨界破壞點(diǎn)即可給出陶土瓦片在風(fēng)致球狀飛射物沖擊下的易損性曲線，如圖11所示.

由能量損傷方程Ddeb≥RD及瓦片破壞的臨界能量值RD=1.45 J推導(dǎo)出瓦片在混凝土飛射物沖擊作用下的易損性判別式如下：

表4 陶土瓦片碰撞結(jié)果

圖11 陶土瓦片的易損性曲線

(5)

(6)

將RD代入式(5)，給出瓦片在一定風(fēng)速的飛射物沖擊作用下所能承受飛射物的臨界質(zhì)量如式(6)，從而可以根據(jù)式(6)繪制易損性曲線.

在易損性曲線的左下方區(qū)域，瓦片是安全的，即瓦片在球狀飛射物的沖擊下不會(huì)發(fā)生破壞.而在易損性曲線的右上方，瓦片將受到飛射物的沖擊而發(fā)生破壞.將數(shù)值模擬結(jié)果與能量損失方程判別結(jié)果對(duì)比可知，數(shù)值模擬得到的易損性曲線與能量判別曲線基本一致，數(shù)值模擬結(jié)果略小于能量判別公式結(jié)果，能量損傷方程可作為材料沖擊破壞的判別標(biāo)準(zhǔn).另一方面，2個(gè)結(jié)果比較一致,也說(shuō)明了本文數(shù)值模擬模型的合理性.

本文數(shù)值模擬模型能較好模擬混凝土球塊沖擊陶土瓦片的沖擊破壞過(guò)程，數(shù)值模擬得到的陶土瓦片的易損性曲線可用于低矮房屋屋面陶土瓦片風(fēng)災(zāi)評(píng)估與設(shè)計(jì)，且其偏于安全.

4 結(jié)論

針對(duì)風(fēng)致飛射物對(duì)建筑圍護(hù)結(jié)構(gòu)的破壞問(wèn)題，基于數(shù)值模擬方法深入研究了球狀混凝土飛射物對(duì)陶土瓦片的破壞損傷機(jī)理，主要結(jié)論如下：

(1) 球狀飛射物-瓦片碰撞數(shù)值模型計(jì)算結(jié)果與落球試驗(yàn)結(jié)果吻合良好，數(shù)值模型能準(zhǔn)確反映沖擊過(guò)程中應(yīng)力、應(yīng)變、位移、能量等關(guān)鍵物理量.

(2) 參數(shù)分析表明瓦片最大等效應(yīng)力與沖擊速率及球塊質(zhì)量成正比，而與瓦片傾角相關(guān)性不大.

(3) 沖擊過(guò)程中球狀飛射物的動(dòng)能呈階梯式迅速衰減，說(shuō)明該過(guò)程是一個(gè)能量顯著轉(zhuǎn)移的過(guò)程，為能量判別標(biāo)準(zhǔn)提供了依據(jù).在此基礎(chǔ)上建立了球狀混凝土飛射物對(duì)陶土瓦片沖擊的能量損傷方程.

(4) 該典型陶土瓦片能承受不超過(guò)34.35 m·s-1的風(fēng)速下球狀混凝土飛射物的沖擊作用.若將該瓦片應(yīng)用于上海地區(qū)，能滿足10年重現(xiàn)期下陣風(fēng)風(fēng)速所致球狀混凝土飛射物的沖擊作用，不能滿足50年和100年重現(xiàn)期下的設(shè)計(jì)要求.

(5) 在大量數(shù)值模擬試驗(yàn)的基礎(chǔ)上，并與能量判別方程對(duì)比，得到了陶土瓦片的易損性曲線，供設(shè)計(jì)與評(píng)估以參考.