梁 瑞，包 娟，周文海，王樹江

(蘭州理工大學(xué) 石油化工學(xué)院，蘭州 730050)

為滿足城市化規(guī)模擴大與基礎(chǔ)設(shè)施建設(shè)，難免在埋地管線附近進(jìn)行爆破施工作業(yè)[1,2]。爆破產(chǎn)生的瞬時高強度沖擊波會對臨近的地下管線造成嚴(yán)重的沖擊，使管道承受較強激振作用而產(chǎn)生應(yīng)力損傷，且在振動頻率等于管道固有頻率時，會引起共振位移破壞[3-6]。張黎明等利用現(xiàn)場爆破試驗，結(jié)合應(yīng)變強度準(zhǔn)則擬合分析得到管道合成峰值振速與埋地管道軸向應(yīng)變呈指數(shù)函數(shù)關(guān)系[7]；王棟等對埋地管道動力響應(yīng)進(jìn)行數(shù)值模擬研究，并通過現(xiàn)場試驗驗證了管道應(yīng)變與地表變化具有比例關(guān)系[8]。Yilmaz T等通過改變爆破參數(shù)對托普恰姆水電站進(jìn)行大量的振動分量監(jiān)測研究[9]。鄭爽英等利用有限元模擬和數(shù)值試驗對影響埋地管道應(yīng)力峰值的管道直徑、壁厚、管道內(nèi)壓等5個爆破影響因素進(jìn)行探究[10]；李鏡培等利用LS-DYNA模擬和現(xiàn)場試驗研究了炸藥藥量和管道壁厚對鋼管破壞的影響情況[11]。郝郁清等結(jié)合屈服強度理論，建立了管道應(yīng)變與對應(yīng)單元振速的關(guān)系，并通過試驗得到爆破振速及爆心距的函數(shù)關(guān)系及最小施工安全距離[12]。馬華原等結(jié)合數(shù)值模擬和實地試驗研究管土耦合模型的作用機理和應(yīng)力應(yīng)變規(guī)律，發(fā)現(xiàn)裂紋尖端的應(yīng)力集中會導(dǎo)致管道開裂[13]。梁博等基于SPH-FEM法建立爆炸模型，分析了管體的擾動和土壤變形情況，得到管體迎爆面和背爆面的應(yīng)變特征[14]。張震等通過現(xiàn)場監(jiān)測和數(shù)值模擬，并結(jié)合最大拉應(yīng)力強度理論確定了管道安全的地表爆破振速[15]。

目前，國內(nèi)外研究學(xué)者對在不同的爆破載荷下埋地管道動力響應(yīng)規(guī)律的研究，主要集中在管道整體的應(yīng)力應(yīng)變和振動速度，但就管道具體迎爆面、背爆面的應(yīng)力應(yīng)變、振速和管道能量的衰減研究較少。利用數(shù)值分析建立管土耦合的爆破荷載有限元模型，并通過管土體拉壓本構(gòu)關(guān)系分析其動力響應(yīng)規(guī)律；在此基礎(chǔ)上引入爆破振動波能量衰減公式，討論埋地管道在爆破載荷下能量的波動衰減規(guī)律，提高了在爆破荷載作用下對埋地管道監(jiān)測的可行性。

1 爆破載荷作用下埋地管道模型的建立

以南京4號地鐵中段兩站區(qū)間的南北兩線隧道開挖掘進(jìn)作為工程背景參考，隧道位于北京西路下方深10～18 m，站間總長約650 m，隧道呈馬蹄形斷面。巖體主要以沉積巖為主，在隧道南北兩線上方1.1～9.7 m近區(qū)范圍內(nèi)，分布有較多的磚混建筑結(jié)構(gòu)。地下1～2 m埋置有地下管線，考慮減小爆破振動對其影響，采用小距離掘進(jìn)的松動爆破。

本文中模型包含的炸藥、巖土體與管道均采用SOLID164類型單元，擬建立與實際工程相似的三維數(shù)值仿真模型。土壤層深4 m，下部巖層深16 m，巖土體整體尺寸為20 m×20 m×2 m；管道在土壤層中埋深2 m，巖土深度10 m處左右兩側(cè)各布置孔徑φ70 mm，炮孔深1 m的圓柱孔炸藥，單次起爆總藥量Q為8.0 kg。兩端同時起爆掘進(jìn)作業(yè)，并對模型均采取實體網(wǎng)格進(jìn)行劃分。數(shù)值模型中管道為空管狀態(tài)，炸藥為連續(xù)耦合裝藥結(jié)構(gòu)并設(shè)置藥包中心位置起爆。在巖體外側(cè)延無限伸區(qū)域及炸藥邊界面進(jìn)行了位移約束與非反射邊界條件，管道兩端設(shè)置固定約束，模型結(jié)構(gòu)如圖1。管道與炸藥、土體接觸界面采取流固耦合算法，炸藥采用ALE算法對單元進(jìn)行流變控制，巖土體及管道均采用LAGRANGE算法控制進(jìn)行求解，求解總時長設(shè)置為300 ms。炸藥選取2#巖石乳化炸藥，其材料類型為*MAT-HIGH-EXPLOSIVE-BURN，炸藥狀態(tài)方程選取不考慮炸藥產(chǎn)物成分的JWL方程[16]。

圖 1 爆破掘進(jìn)與淺埋管道結(jié)構(gòu)模型(單位：m)Fig. 1 Blasting driving and shallow pipe structure model(unit：m)

(1)

式中：ρe為炸藥密度；D為爆速；A、B、R1、R2均為材料參數(shù)；其中R1、R2為無量綱參數(shù)；ω為格林愛森常量；PCJ為爆轟壓力；V0為比體積；e為初始比內(nèi)能。見表1。

假定巖體為滿足各向同性的彈塑性材料，選用*MAT-PLASTIC-KINEMATIC；土體選用材料類型為*MAT-SOIL-AND-FOAM。見表2。

表 1 炸藥類型及狀態(tài)方程參數(shù)Table 1 Types of explosives and parameters of equation of state

表 2 巖土體材料參數(shù)Table 2 Material parameters of rock and soil

表中：ρ為密度；E為彈性模量；G為剪切模量；μ為泊松比；σy為屈服極限；β為巖體硬化系數(shù)。

管道選用*MAT-PLASTIC-KINEMATIC作為低碳合金鋼制管道的材料類型。見表3。

表中，ρ為密度；E為彈性模量；μ為泊松比；σy為屈服應(yīng)力；Etan為剪切模量。

2 管道動力響應(yīng)研究

2.1 模擬試驗設(shè)計

見表4。

表 4 模型對照組Table 4 Model comparison group

2.2 管道受載應(yīng)力分布分析

為了更好地分析爆破載荷下埋地管道的動力響應(yīng)，管道的應(yīng)力應(yīng)變簡化為理想線彈性與理想剛塑性兩種不同狀態(tài)下的本構(gòu)關(guān)系。取管道上任意微小單元，其應(yīng)力分量的大小和方向均已知如圖2所示。

圖 2 四面體微元受力圖Fig. 2 Stress diagram of tetrahedron element

物體的剪切破壞主要由應(yīng)力偏張量引起，定義Si為偏主應(yīng)力，故偏應(yīng)力主平面方程為

(2)

得到應(yīng)力偏張量不變量

(3)

由上式可知，J2由三個方向的應(yīng)力變化決定，考慮了中間主應(yīng)力和平均應(yīng)力，并且靜水壓力不僅對平均應(yīng)力有影響，還會引起屈服破壞。由于管道的屈服破壞主要由偏應(yīng)力第二不變量J2引起，而偏應(yīng)力第二不變量主要反映剪應(yīng)力的大小，故根據(jù)J2函數(shù)形式定義管道的應(yīng)力破壞強度指標(biāo)，即有效應(yīng)力或廣義剪應(yīng)力σi為

(4)

(1)管道底部-迎爆面軸向應(yīng)力分布規(guī)律

應(yīng)力波于8～10 ms時作用于管道迎爆面兩端并迅速沿管道壁面的軸向及環(huán)向傳播，如圖3有效應(yīng)力云圖可以看到軸向有效應(yīng)力大于環(huán)向有效應(yīng)力。考慮到管道的破壞主要由軸向和環(huán)向的應(yīng)力波的動力響應(yīng)造成，故下文對管道的軸向和環(huán)向動力響應(yīng)進(jìn)行研究。

圖 3 管道有效應(yīng)力云圖Fig. 3 Cloud chart of effective stress of pipeline

管道底部作為迎爆面，為爆破沖擊載荷最強的部位。沿軸線方向從管道底部左側(cè)至右側(cè)依次選取5個等距5 m的單元作為測點如圖4(a)。提取各單元有效應(yīng)力數(shù)據(jù)，繪制有效應(yīng)力時程曲線。

對比3組不同管壁厚度的管組，分析受載管道軸向的應(yīng)力波分布。如圖5所示，(a)、(b)、(c)分別表示壁厚為4 mm、6 mm、9 mm管道的有效應(yīng)力Effective Stress。

圖 4 壁厚為4 mm的監(jiān)測單元示意圖Fig. 4 Schematic diagram of monitoring unit with 4 mm wall thickness

圖 5 不同壁厚管道迎爆面軸向單元有效應(yīng)力時程曲線Fig. 5 Time history curve of effective stress of axial element on the blasting face of pipeline with different wall thickness

分析迎爆面-管道底部有效應(yīng)力時程曲線，發(fā)現(xiàn)隨著管道壁厚的增大，有效應(yīng)力峰值逐漸減小，說明管道壁厚在一定范圍內(nèi)，壁厚的增大可以有效增大管道的波阻抗，使管壁單元的有效應(yīng)力減小，管道越不易發(fā)生破壞。由于B、D兩單元距離爆源最近，在15～20 ms時的初始應(yīng)力峰值均為最大值。

(2)管道頂部-背爆面軸向應(yīng)力分布規(guī)律

同理提取背爆面監(jiān)測單元，繪制不同壁厚管道組的有效應(yīng)力時程曲線。如圖6所示，(a)、(b)、(c)分別表示壁厚為4 mm、6 mm、9 mm管道背爆面軸向單元有效應(yīng)力時程曲線。

分析背爆面-管道頂部有效應(yīng)力時程曲線可以得到，隨著管道壁厚的增大，有效應(yīng)力峰值仍然逐漸減小，管道頂部作為背爆面，為管道最后受到爆破沖擊載荷的部位。同時，對比迎爆面的管底單元有效應(yīng)力峰值發(fā)現(xiàn)，背爆面受到應(yīng)力波正向傳播和地面反射波共同作用，產(chǎn)生應(yīng)力局部放大效應(yīng)。

2.3 管道受載速度分布分析

由于管道單元和巖體受沖擊載荷發(fā)生非線性振動，且瞬時沖擊載荷下應(yīng)力變化過快，監(jiān)測難度大，故選擇分析更易監(jiān)測的振動強度因子振動速度。

管道底部-頂部(迎爆面-背爆面)軸向速度分布規(guī)律。

同方法選取了沿管道軸線方向頂部單元，提取各單元振動速度數(shù)據(jù)，繪制振動速度時程曲線，分別表示不同管道壁厚迎爆面軸向的振動速度時程曲線，如圖7、圖8所示。

圖 6 不同壁厚管道背爆面軸向單元有效應(yīng)力時程曲線Fig. 6 Time history curve of effective stress of axial element on back blasting surface of pipeline with different wall thickness

圖 7 不同壁厚管道軸向迎爆面單元速度時程曲線Fig. 7 Time history curve of unit velocity on the axial blast surface of pipeline with different wall thickness

圖 8 不同壁厚管道軸向背爆面單元速度時程曲線Fig. 8 Unit velocity time history curve of axial back blasting surface for different wall thickness pipeline

分析管道迎爆面-背爆面上監(jiān)測單元的振動合成速度時程曲線，得到在10～21 ms時，振動速度迅速達(dá)到峰值，并逐漸發(fā)生衰減。管道兩側(cè)1/4處(B、D監(jiān)測單元)由于與爆源的直線距離最短，其振動速度峰值為時程內(nèi)最高值。隨著管壁厚度的增大，峰值振速呈減小趨勢；管道壁厚為4 mm的峰值速度的29%，而管道壁厚為6 mm的衰減速率為7%，管道壁厚越大，振動速率衰減速率越小。對比迎背爆面振速分布得到，不同壁厚管道振速由迎爆面至背爆面具有不同的衰減，但振速衰減相比應(yīng)力衰減具有一定線性關(guān)系。

2.4 管道受載速度分布分析

為了直觀對比探究管道在爆破振動波載荷作用下的應(yīng)力分布、響應(yīng)狀態(tài)及傳播規(guī)律，現(xiàn)分別繪制沿著管道軸向的迎爆面-背爆面，不同爆心距R下的有效應(yīng)力峰值、振動合速度峰值的變化曲線圖分析不同壁厚的管道軸向整體的受載振動響應(yīng)如圖9、圖10所示。

(1)管道迎-背爆面有效應(yīng)力峰值分布規(guī)律

分別以管道軸向長度作為x坐標(biāo)，管道的迎爆面-背爆面單元的有效應(yīng)力峰值為y坐標(biāo)，分析三種不同壁厚的軸向迎-背爆面應(yīng)力峰值變化規(guī)律如圖9所示。

由圖可知管壁厚度越大，有效應(yīng)力峰值越小，衰減程度越小。管道兩端受爆破載荷同時作用，應(yīng)力波由兩端向中間傳播的過程中逐漸衰減，在管道中部應(yīng)力峰值最小，薄壁管道應(yīng)力峰值的衰減速率大于厚壁管道應(yīng)力峰值的衰減速率。

圖 9 軸向迎-背爆面應(yīng)力峰值Fig. 9 Peak value of stress on the axial front back blasting surface

(2)管道迎-背爆面速度峰值分布規(guī)律

以管道軸向長度作為x坐標(biāo)，管道的迎爆面-背爆面單元的振動合速度峰值為y坐標(biāo)，分析三種不同壁厚的軸向迎-背爆面速度峰值變化規(guī)律如圖10所示。

圖 10 軸向迎-背爆面速度峰值Fig. 10 Peak value of axial velocity of front and back blasting surface

由圖10可知迎爆面單元振動速度峰值較背爆面振動速度峰值更大，管道兩側(cè)1/4處振動速度最大，中間位置最小。結(jié)合圖9的有效峰值應(yīng)力分布，管道兩端距離爆源最近，應(yīng)力峰值最大，但速度峰值最大位置為管道兩端的1/4處。說明振動速度的傳播規(guī)律不是嚴(yán)格按照應(yīng)力的線性傳播進(jìn)行衰減，更易受其他因素影響，如波的入射角，振動波以球面波傳播至管道壁面上時，根據(jù)距離計算，速度峰值最大處與爆源的夾角為32°。

3 爆破振動波能通量計算

對管道的動力響應(yīng)研究，除了對應(yīng)力和振速的分析，更應(yīng)考慮能量的改變，振動波能量傳播規(guī)律對受載結(jié)構(gòu)的響應(yīng)分析能夠綜合考慮爆破振動強度、振動頻率、持續(xù)時間的多因素共同作用。

爆破振動波傳播是一個行進(jìn)的點對點擾動過程，同時點的積分伴隨著能量的傳遞。爆破振動強度的可量化因素通常以爆破振動強度因子進(jìn)行運動描述，如質(zhì)點振動位移、速度、加速度等。對于爆破強度A，其廣義表達(dá)式為

A=KQαRβ

(5)

式中：A為爆破地震強度物理量；K為場地系數(shù)；Q為炸藥量；α、β為衰減指數(shù)；R為測點至爆心距離。

對于某監(jiān)測單元的爆破振動波能量的計算可以通過爆破振動波傳播過程中經(jīng)過該監(jiān)測單元的能通量的積分求得。假設(shè)爆破振動波為無限連續(xù)介質(zhì)中傳播的球面波，在球面坐標(biāo)系中應(yīng)力張量的分量表示為

(6)

式中：u1為振動位移；r為爆心距；λ和μ為拉姆常數(shù)。由求和定則，能通量求解為

Φ=τijnjvj

(7)

而單元振動速度可由后處理中提取。設(shè)通過爆心距r處的能通量為定值，通過該單元截面的總能量P為

P=4πr2Φ

(8)

(9)

(10)

爆破振動波的能量隨著波陣面繼續(xù)擴展，圍巖介質(zhì)的開裂消耗了部分能量，所經(jīng)過的質(zhì)點振速逐漸降低，其微分表達(dá)式為

(11)

式中:Es為監(jiān)測單元的爆破振動波能量；β表示振動波能量衰減量系數(shù)，負(fù)號表示能量的衰減。再對表達(dá)式積分得

Es=E0e-βr

(12)

式中:E0為爆源炸藥能量轉(zhuǎn)換為爆破振動波的初始能量；r為爆心距。又聯(lián)合式E0=ηEe，η代表炸藥能與轉(zhuǎn)化成振動波初始能的百分比。得到爆破振動能量衰減公式

Es=ηEee-βr

(13)

對比式(5)，可以發(fā)現(xiàn)式(13)中爆破振動波能量Es對應(yīng)強度物理量A；η對應(yīng)場地系數(shù)K，而η與Ee的積表達(dá)為裝藥量Q的關(guān)系；式中的指數(shù)部分e-βr表達(dá)為能量的衰減形式，兩公式間具有一致性。

通過爆破振動波衰減式(13)對表4中3組公稱直徑相同壁厚不同的管道環(huán)向監(jiān)測單元分析，如圖4(b)所示，沿管壁面底部至頂部依次選取8個單元并以字母A～H排序，監(jiān)測單元具體數(shù)據(jù)的見表5。計算爆心距r并與式(10)計算得到的爆破振動波能量Es進(jìn)行擬合回歸。

首先將式(13)兩邊取對數(shù)

lnEs=lnη+lnEe+(-βr)

(14)

令y=lnEs，x=r，α=-β，b=lnη+lnEe，線性化表示為

表 5 管壁單元峰值振速與振動波能量Table 5 Peak vibration velocity and vibration wave energy of tube wall element

y=ax+b

(15)

表中：r為爆心距；∑v2(ti)為單元時程速度平方和。

將表5數(shù)據(jù)代入式(15)進(jìn)行線性回歸擬合，得到管道中部環(huán)向監(jiān)測單元的線性回歸分析結(jié)果如圖11。分析發(fā)現(xiàn)3組管道的振動波能量Es隨爆心距r增大而減??；不同壁厚的管道能量衰減程度不同：壁厚為4 mm、6 mm、9 mm的管道，爆破振動波能量衰減量分別為142.71 MJ、125.7 MJ、66.88 MJ，兩者近似呈負(fù)相關(guān)關(guān)系。分析3組管道振動波能量的衰減變化，發(fā)現(xiàn)在同種爆破條件下，管道壁厚增大2 mm和3 mm的爆破振動波能量衰減速率分別減少約12%和47%，隨著管壁厚度增大，振動波能量在管道上的衰減幅度也下降。

嘗試探尋管道監(jiān)測單元振動波能量和對應(yīng)振速的關(guān)系，將振動波能量與單元峰值速度進(jìn)行擬合，得到爆破振動波能量Es與振動速度的平方v2近似正相關(guān)，且相關(guān)性均大于0.85，如圖12。

圖 11 能量與爆心距線性回歸結(jié)果Fig. 11 Results of linear regression of energy and blast center distance

圖 12 能量與速度平方線性回歸結(jié)果Fig. 12 Energy and velocity squared linear regression results

為比較由衰減公式(13)擬合得到的衰減系數(shù)β與數(shù)值模擬結(jié)果計算得到的衰減率貼近程度，定義迎-背爆面振動波能量衰減率κ

κ=(EsA-EsH)/EsA

(16)

式中:EsA、EsH分別為管道迎爆面單元和背爆面單元的瞬時振動波能量。對管道迎爆面與背爆面振動波能量數(shù)據(jù)計算得到κ作為能量衰減率，通過能量衰減公式對三組不同壁厚的管道單元擬合回歸得到η、β值，令i=T/D作為管道厚徑比參數(shù)，如表6，κ值由表5計算得到。

表 6 管道回歸結(jié)果與能量衰減公式Table 6 Pipeline regression results and energy attenuation formula

根據(jù)表6，繪制η、κ、β和厚徑比的變化關(guān)系圖如圖13。

圖 13 i與η、κ、β變化關(guān)系Fig. 13 Relationship between i and changes of η,κ,β

分析得管道厚徑比i<0.0525的管道，炸藥能的衰減更為顯著；厚徑比0.0525≤i≤0.05375，管道能量衰減系數(shù)β、炸藥能轉(zhuǎn)換率η與迎-背爆面能量衰減率κ的能量衰減程度均處于26.7%；厚徑比i≥0.05375，管道振動波能量的衰減更加明顯；管道厚徑比i≥0.06時，管道環(huán)向的振動波能量衰減規(guī)律基本相似。隨著管道厚徑比i增大，管道環(huán)向監(jiān)測單元的η、κ與β均逐漸減小，管道壁面的振動波阻抗增大且管道對振動的響應(yīng)敏感度降低，振動波能量衰減程度變?nèi)酢１普駝硬芰克p模型關(guān)于爆心距回歸得到的衰減率β，其數(shù)值及變化趨勢與測點數(shù)據(jù)計算得到的κ貼近程度較高，表明在爆破振動波對管道的動力響應(yīng)分析中建立的能量衰減模型較為準(zhǔn)確。

4 結(jié)論

通過模擬在爆破載荷作用下，不同壁厚埋地管道的動力響應(yīng)，對其有效應(yīng)力、振動速度和能量變化規(guī)律進(jìn)行分析：

(1)通過分析迎爆面-背爆面軸向監(jiān)測單元應(yīng)力時程曲線和速度時程曲線，發(fā)現(xiàn)管道迎爆面1/4(B、D單元)處存在應(yīng)力集中現(xiàn)象，在20～50 ms時管道有效應(yīng)力達(dá)到峰值；管道背爆面有局部應(yīng)力放大效應(yīng)；在10～21 ms時管道迎爆面和背爆面的振動速度達(dá)到峰值。

(2)隨著管道壁厚的增加，管道迎爆面和背爆面同一位置處的軸向有效應(yīng)力和速度均逐漸減小。在管道兩端，迎爆面和背爆面的有效應(yīng)力峰值最大；距管道兩端1/4處，迎爆面和背爆面的振動速度峰值最大。

(3)建立管道振動波能量衰減模型，擬合得到管道環(huán)向振動波能量衰減量與管道壁厚呈近似負(fù)相關(guān)關(guān)系，管壁振速的平方與振動波能量呈一次函數(shù)關(guān)系；當(dāng)管道厚徑比i<0.0525時，炸藥能的衰減最快，0.0525≤i≤0.05375時，衰減指標(biāo)η、β、κ數(shù)值均為26.7%，且i≥0.06時管道環(huán)向監(jiān)測單元各項衰減指標(biāo)規(guī)律趨于一致。

(4)在本文工況下，模擬所得的管道振速峰值分布于1.0～2.25 cm/s范圍內(nèi)，將《爆破安全規(guī)程》中爆破振動波的振動速度與振動頻率作為安全判據(jù)，管道處于安全狀態(tài)。