米中陽，張智宇，黃永輝，雷 振

(1. 昆明理工大學(xué) 國土資源工程學(xué)院，昆明 650093；2. 昆明理工大學(xué)爆破新技術(shù)應(yīng)用研究所，昆明 650093；3. 貴州理工學(xué)院 礦業(yè)工程學(xué)院，貴陽 550003)

隨著社會的發(fā)展，爆破技術(shù)得到了廣泛的應(yīng)用，爆破誘發(fā)的一系列的負(fù)面效應(yīng)尤其是爆破地震效應(yīng)問題受到了人們普遍的關(guān)注和重視[1]。曹孝君等通過對各種爆破條件下的爆破地震波的特性分析和建筑物破壞現(xiàn)象、結(jié)構(gòu)動力響應(yīng)參數(shù)和結(jié)構(gòu)動力特性參數(shù)的測試工作，對爆破振動在爆破地震波的特性、傳播規(guī)律及爆破地震波與建筑物動力響應(yīng)的關(guān)系等方面做了大量的研究，并且取得了一定的成果[2]。

爆破地震波對建筑物的影響本質(zhì)上是結(jié)構(gòu)的動力響應(yīng)問題，建筑物在爆破地震波荷載作用下，由于基礎(chǔ)運(yùn)動而引起結(jié)構(gòu)振動，爆破振動破壞效應(yīng)不僅與地表振動強(qiáng)度參數(shù)有關(guān)，同時與爆破地震波頻率及振動持續(xù)時間也有較大的相關(guān)性，爆破振動破壞效應(yīng)是結(jié)構(gòu)體動態(tài)破壞問題，其對結(jié)構(gòu)體的破壞效應(yīng)實(shí)際上是一個動力響應(yīng)過程[3]。高富強(qiáng)[4]等基于結(jié)構(gòu)動力學(xué)理論，利用ANASYS/LS-DYNA軟件構(gòu)建了數(shù)值模型，綜合了地震波的振幅、頻率及結(jié)構(gòu)特征等對結(jié)構(gòu)振動的影響，對結(jié)構(gòu)在爆破振動作用下的響應(yīng)過程進(jìn)行了系統(tǒng)分析。黃志強(qiáng)[5]等通過采用結(jié)構(gòu)動力學(xué)原理，對結(jié)構(gòu)在爆破振動激勵下的動力響應(yīng)進(jìn)行了計算分析，并運(yùn)用ANSYS軟件對煙囪結(jié)構(gòu)在爆破振動下的響應(yīng)情況進(jìn)行了數(shù)值模擬，結(jié)果表明，隨著煙囪高度的增加，其位移響應(yīng)也越大。魏曉林[6]的研究表明，炸藥量和爆心距是影響振動的主要因素，建筑物自身頻率對爆破振動的放大程度起著決定性作用。

由于礦山爆破裝藥量大、爆破次數(shù)頻繁，與其他爆破相比爆破的地震效應(yīng)較強(qiáng)，對周圍建筑的影響較大。隨著某礦礦體開采逐步向礦辦公區(qū)及住宅區(qū)和周邊村莊推進(jìn)，以及多次爆破振動會對建筑結(jié)構(gòu)產(chǎn)生累積振損，該礦山爆破對周邊建筑結(jié)構(gòu)產(chǎn)生的危害越來越大，有必要對礦山爆破對鄰近框架建筑結(jié)構(gòu)的影響進(jìn)行研究，從而為建筑物的現(xiàn)場監(jiān)測提供一定依據(jù)。

本文采用動力有限元分析方法，對礦山爆破振動載荷下的建筑結(jié)構(gòu)動力響應(yīng)進(jìn)行了研究，通過有限元軟件ANSYS/LS-DYNA建立該礦辦公樓五層混凝土框架結(jié)構(gòu)有限元模型，將現(xiàn)場實(shí)測的爆破振動波形輸入建筑基礎(chǔ)，分析了建筑結(jié)構(gòu)的動力響應(yīng)情況，并參照相關(guān)標(biāo)準(zhǔn)[7]判斷建筑結(jié)構(gòu)的破壞程度。

1 辦公樓數(shù)值模型構(gòu)建

1.1 材料模型的選取及約束條件

對筋混凝土的密度取為2 500 kg/m3，梁、柱采用C30混凝土，彈性模量取值E=30 GPa，樓板采用C25混凝土，彈性模量取值E=28 GPa，混凝土材料泊松比取為0.2。對于鋼筋混凝土結(jié)構(gòu)，結(jié)構(gòu)阻尼比ζ取0.05。計算中未考慮土與上部結(jié)構(gòu)共同作用，基礎(chǔ)與地基剛性連接，即約束與地面接觸的框架柱以及剪力墻上節(jié)點(diǎn)的所有自由度。

1.2 三維有限元模型

建立辦公樓有限元模型如圖1所示。在該模型中，節(jié)點(diǎn)總數(shù)為22 330個，單元總數(shù)為26 838個。模型水平面中以模型的長邊方向?yàn)閤軸方向，短邊方向?yàn)閥軸方向，以豎直向上方向?yàn)閦軸方向。此外，該計算分析模型未考慮填充墻的影響。

圖1 某辦公樓有限元模型Fig.1 Finite element model of phosphate office building

2 辦公樓三維有限元動力分析

2.1 提取計算分析用爆破振動波

二采區(qū)2050-2040水平和三采區(qū)2110-2100水平爆破時，對礦辦樓1樓、3樓和5樓進(jìn)行了兩次實(shí)測。對比二采區(qū)開采爆破時所監(jiān)測的辦公樓1樓振動波速曲線和三采區(qū)開采爆破時所監(jiān)測的辦公樓1樓振動波速曲線，發(fā)現(xiàn)二采區(qū)開采爆破時所監(jiān)測的辦公樓1樓振動波速相對較大。因此，選取如圖2所示二采區(qū)開采爆破時所監(jiān)測的辦公樓1樓振動波速曲線作為該礦辦公樓動力響應(yīng)計算分析用爆破振動波。

圖2 該礦辦公樓動力響應(yīng)分析用爆破振動速度波Fig.2 Blasting vibration velocity wave for dynamic response analysis of office building in this mine

該礦辦公樓動力響應(yīng)分析采用10 s時間歷程爆破振動加速度時程曲線分別如圖3至圖5所示。

圖3 水平徑向加速度時程曲線Fig.3 Horizontal radial acceleration time-history curve

圖4 水平切向加速度時程曲線Fig.4 Horizontal tangential acceleration time-history curve

圖5 垂直方向加速度時程曲線Fig.5 Vertical acceleration time-history curve

將實(shí)際的爆破振動波輸入爆破振動波加速度時程曲線，并利用公式(1)和公式(2)將加速度時程轉(zhuǎn)換為應(yīng)力時程。

σs=-2(ρCs)νs

(1)

σn=-2(ρCp)νn

(2)

式中，σs、σn分別為施加在靜態(tài)邊界上的法向應(yīng)力和切向應(yīng)力，ρ為介質(zhì)密度，Cp和Cs分別為介質(zhì)P波和S波的波速，νn和νs分別為結(jié)構(gòu)模型邊界上法向和切向的速度分量。

2.2 辦公樓模型的自振特性

為確定該礦辦公樓模型的自振特性，研究應(yīng)用ANSYS軟件，采用子空間迭代方法計算獲得模型的頻率對辦公樓模型進(jìn)行了模態(tài)分析，并提取了結(jié)構(gòu)前40階頻率，由頻率和周期的對應(yīng)關(guān)系，得到辦公樓結(jié)構(gòu)的自振周期如圖6所示。

圖6 結(jié)構(gòu)自振周期隨階數(shù)變化曲線Fig.6 The curve of natural vibration period with order of structure

從圖6可知，結(jié)構(gòu)前四階自振周期逐步下降，第五階至第七階自振周期較接近，7階以后自振周期逐漸變短，第16階自振周期以后趨于平穩(wěn)變化，階數(shù)越大，自振周期緩慢變短。其中，模型的高階振型出現(xiàn)扭轉(zhuǎn)效應(yīng)。

2.3 爆破振動作用下辦公樓振動響應(yīng)

2.3.1 動力響應(yīng)時程分析

根據(jù)在進(jìn)行動力響應(yīng)時程分析時，采用NEWMARK-β法進(jìn)行數(shù)值積分，迭代求解。為提高迭代精度，保證迭代收斂性，將地震波或震動波記錄的時間間隔0.125 s，改為取迭代步長為0.0125 s；在加載過程中采用逐級加載的方法進(jìn)行加載。

在進(jìn)行動力響應(yīng)時程分析時，鑒于結(jié)構(gòu)阻尼對結(jié)構(gòu)的動力響應(yīng)影響很大，阻尼比的影響基本上反比于位移反應(yīng)。因此，在選取阻尼時，采用輸入Rayleigh阻尼(瑞利阻尼)的方法。即選取瑞利阻尼為：

C=αM+βK

(3)

式中：C為動力方程中的阻尼矩陣，K為剛度矩陣，M為質(zhì)量矩陣，α為與質(zhì)量成比例的阻尼常數(shù)；β為與剛度成比例的阻尼常數(shù)，αM為質(zhì)量分量，βK為剛度分量。α、β可由公式(4)確定，即：

(4)

式中：?i、?j分別為結(jié)構(gòu)第i、j振型的圓頻率；ξi、ξj為結(jié)構(gòu)第i、j振型的阻尼比。對于混凝土結(jié)構(gòu)，兩個振型假設(shè)具有相同的阻尼比ξi=ξj=0.05。通過模態(tài)分析提取了結(jié)構(gòu)的頻率ω，根據(jù)式(4)，然后可以確定阻尼參數(shù)α和β。

2.3.2 爆破振動作用下的結(jié)構(gòu)加速度、速度、位移時程曲線分析

由于篇幅的原因加速度、速度時程曲線省略，分析過程三者基本相同，此處僅以位移時程曲線為例進(jìn)行了分析。結(jié)構(gòu)受爆破振動水平徑向加速度作用，在0.875 s時刻(爆破振動水平徑向最大加速度時刻)結(jié)構(gòu)頂層的位移時程曲線如圖7、8所示，結(jié)構(gòu)徑向位移和切向位移最大值分別為1.59×10-6m和2.17×10-6m。

圖7 水平徑向結(jié)構(gòu)頂層位移時程曲線Fig.7 The top displacement time-history curve of the horizontal radial direction structure

圖8 水平切向結(jié)構(gòu)頂層位移時程曲線Fig.8 The top displacement time-history curve of the horizontal tangential direction structure

結(jié)構(gòu)受爆破振動水平切向加速度作用，在9.25 s時刻(爆破振動水平切向最大加速度時刻)結(jié)構(gòu)頂層的位移時程曲線如圖9、10所示，結(jié)構(gòu)徑向位移和切向位移最大值分別為2.21×10-6m和1.55×10-6m，均小于最小安全值。

圖9 水平徑向結(jié)構(gòu)頂層位移時程曲線Fig.9 The top displacement time-history curve of the horizontal radial direction structure

圖10 水平切向結(jié)構(gòu)頂層位移時程曲線Fig.10 The top displacement time-history curve of the horizontal tangential direction structure

結(jié)構(gòu)受爆破振動垂直加速度作用，在5 s時(爆破振動水平切向最大加速度時刻)結(jié)構(gòu)頂層的位移時程曲線如圖11所示，結(jié)構(gòu)垂直位移和切向位移最大值為4.92×10-6m。

圖11 豎直方向結(jié)構(gòu)頂層位移時程曲線Fig.11 The top displacement time-history curve of the upright direction structure

通過對爆破振動作用下辦公樓受爆破振動各個方向上最大加速度作用下結(jié)構(gòu)頂層的加速度、速度和位移的時程曲線進(jìn)行分析，得出最大值，均小于安全值。

通過對加速度、速度和位移時程曲線分析，可以對加速度、速度和位移隨時間的變化過程有一個充分的認(rèn)識和理解，通過對其分析可以對框架結(jié)構(gòu)建筑物對爆破震動波的響應(yīng)過程有一個更清晰的認(rèn)識。

3 辦公樓結(jié)構(gòu)動力響應(yīng)程度分析

在二采區(qū)開采爆破時所監(jiān)測的辦公樓1樓振動波速(水平切向)作用下，計算的監(jiān)測點(diǎn)處響應(yīng)加速度時程曲線如圖12所示。

圖12 二采區(qū)開采爆破時計算的監(jiān)測點(diǎn)處響應(yīng)加速度時程曲線Fig.12 Response time curves of monitoring points at the time of blasting in the two mining areas

根據(jù)結(jié)構(gòu)動力響應(yīng)規(guī)律，常定義加速度放大系數(shù)η為結(jié)構(gòu)加速度動峰值與結(jié)構(gòu)基礎(chǔ)加速度動峰值的比值。因此，如果結(jié)構(gòu)體內(nèi)任意一點(diǎn)E的動力響應(yīng)的加速度峰值為aE，基礎(chǔ)處的動力響應(yīng)的加速度峰值為a，則E點(diǎn)的加速度放大系數(shù)η可用式(5)表示，即有：

(5)

由圖12和式(5)可得辦公樓2層、3層、4層、5層、6層對應(yīng)監(jiān)測點(diǎn)處的加速度放大系數(shù)，如表1所示。

由表1可知，除2層對應(yīng)監(jiān)測點(diǎn)處的爆破振動響應(yīng)加速度縮小外，其余樓層對應(yīng)監(jiān)測點(diǎn)處的爆破振動響應(yīng)加速度均放大。3層、4層、5層、6層對應(yīng)監(jiān)測點(diǎn)處的爆破振動響應(yīng)加速度放大系數(shù)分別為26%、17%、34%和69%。但各樓層對應(yīng)監(jiān)測點(diǎn)處的爆破振動響應(yīng)加速度均較小，最大值均未超過6.00×10-4m/s2。

表1 結(jié)構(gòu)相關(guān)部位加速度放大系數(shù)一覽表(有阻尼)Table 1 List of acceleration magnification coefficients for structure related parts (damping)

二采區(qū)開采爆破在辦公樓引起的振動響應(yīng)總體較小，對辦公樓整體穩(wěn)定性基本無影響。除2層對應(yīng)監(jiān)測點(diǎn)處的爆破振動響應(yīng)加速度縮小外，其余樓層對應(yīng)監(jiān)測點(diǎn)處的爆破振動響應(yīng)加速度放大效應(yīng)并隨著樓層的增高逐漸增大的趨勢，6層放大效應(yīng)達(dá)到最大。

4 結(jié)論

1)在二采區(qū)開采爆破振動作用下，辦公樓結(jié)構(gòu)前4階自振周期逐步下降，第5階至第7階自振周期較接近， 7階以后隨著階數(shù)越大自振周期逐漸變短，第16階自振周期以后趨于平穩(wěn)變化，模型的高階振型出現(xiàn)扭轉(zhuǎn)效應(yīng)，使得整個建筑物的安全性和穩(wěn)定性都大幅下降。

2)通過對爆破振動作用下辦公樓受爆破振動各個方向上最大加速度作用下結(jié)構(gòu)頂層的加速度、速度和位移的時程曲線進(jìn)行分析，得出最大值，均小于安全值。通過時程曲線可以對加速度、速度和位移進(jìn)行分析，從而為判斷建筑物的安全性提供一定的參考。

3)除2層外框架結(jié)構(gòu)對礦山爆破的加速度放大效應(yīng)隨樓層有逐漸變大的趨勢，頂層的加速度放大效應(yīng)最大。因此，礦山爆破對大于2層框架建筑物有影響，應(yīng)加強(qiáng)對頂層的監(jiān)測。對于2層加速度的放大效應(yīng)出現(xiàn)異常變小現(xiàn)象，需要做進(jìn)一步研究和分析。