樊玉光，方軍，崔瑩，黨文波，王嘉樂(lè)

西安石油大學(xué)機(jī)械工程學(xué)院，陜西 西安 710065

近年來(lái)各種意外爆炸事故頻發(fā)，加之戰(zhàn)爭(zhēng)與恐怖襲擊的威脅，使得爆炸荷載已成為影響埋地油氣管道安全的重要第三方荷載[1-4]。為了提高管道的抗爆性能，合理包裹吸能材料是一個(gè)有效的思路。碳纖維復(fù)合材料(carbon fiber reinforced polymer, CFRP)由于其重量輕、耐腐蝕性、耐沖擊性等優(yōu)點(diǎn)很適合用于油氣管道加固[5-8]。因此，研究埋地管道在地面爆炸作用下的力學(xué)響應(yīng)以及防護(hù)很有意義。

許多學(xué)者對(duì)油氣管道在爆炸荷載下的力學(xué)響應(yīng)及防護(hù)問(wèn)題進(jìn)行了大量研究，夏宇磬等[9]結(jié)合爆炸應(yīng)力波傳播理論，推導(dǎo)了可表征爆炸應(yīng)力波穿過(guò)土-巖界面后施加至埋地管道的爆破振動(dòng)荷載衰減公式。王亞光等[10]針對(duì)飽和土中鋼頂管在地面爆炸荷載作用下的動(dòng)力響應(yīng)問(wèn)題展開(kāi)研究，探究了炸藥量、起爆位置、鋼頂管埋置深度以及頂管半徑等參數(shù)對(duì)于鋼頂管位移響應(yīng)之間的關(guān)系。龔相超等[11]設(shè)計(jì)和實(shí)施了系列爆炸波作用下管道的動(dòng)態(tài)響應(yīng)實(shí)驗(yàn)，分析管道在不同的爆心距和爆源埋深條件下，不同位置測(cè)點(diǎn)峰值應(yīng)變的衰減規(guī)律。RIGAS等[12]針對(duì)爆炸荷載下加壓天然氣管道地面-管道相互作用以及管道的響應(yīng)展開(kāi)研究，提出了評(píng)估爆破安全距離的實(shí)用計(jì)算方法。MOKHTARI等[13]通過(guò)數(shù)值模擬對(duì)CFRP加固埋地管道展開(kāi)研究，研究表明CFRP包覆可以顯著提高鋼管在爆炸載荷下的性能。HAJIAZIZI等[14]通過(guò)數(shù)值模擬對(duì)GFRP加固埋地管道展開(kāi)研究，結(jié)果表明，GFRP可以顯著減少爆炸引起的管道變形，使管道的損傷區(qū)長(zhǎng)度、截面變形和最大縱向應(yīng)變分別降低59%、48%和73%。

通過(guò)上述國(guó)內(nèi)外相關(guān)研究可知，雖然驗(yàn)證了CFRP加固的有效性，但目前針對(duì)油氣管道與土耦合響應(yīng)數(shù)值模擬方法，以及CFRP防護(hù)優(yōu)化的研究相對(duì)較少。本文基于流固耦合的方法表述管道與土的相互作用，通過(guò)試驗(yàn)研究和數(shù)值模擬，探討了埋地管道在爆炸荷載作用下的力學(xué)響應(yīng)以及CFRP的防護(hù)作用，并進(jìn)一步討論了CFRP層數(shù)、炸藥量、爆心距等因素對(duì)管道響應(yīng)的影響。

1 淺埋爆炸試驗(yàn)

1.1 試驗(yàn)材料

試驗(yàn)選擇與實(shí)際工程中所采用的X80管線鋼屈服強(qiáng)度一致的N80油管作為試驗(yàn)對(duì)象，TNT炸藥量為9.35kg，如圖1所示。N80管道試樣的具體參數(shù)如表1所示。將預(yù)拌和好的碳纖維浸漬膠涂抹在N80管道上，選擇300g一級(jí)CFRP布進(jìn)行粘貼加固，CFRP具體參數(shù)如表2所示。

圖1 普通管道與CFRP管道Fig.1 Plain pipe and CFRP pipe

表1 N80管道試件基本參數(shù)表

表2 CFRP基本參數(shù)表

1.2 試驗(yàn)布置

如圖2所示，在試驗(yàn)場(chǎng)預(yù)先開(kāi)挖1200mm×900mm×800mm(長(zhǎng)×寬×深)的埋設(shè)用坑，將兩組試樣及反力架一同放入坑內(nèi)，土體回填包括兩個(gè)階段，第一次回填界限依據(jù)爆心距為440mm的標(biāo)準(zhǔn)確定以埋設(shè)炸藥；第二次回填界限為自然地面，通過(guò)插入鋼棒方便量測(cè)。試件實(shí)際就位如圖3所示。

圖2 試驗(yàn)剖面圖 圖3 試件就位圖Fig.2 Test section Fig.3 Test piece in position

1.3 試驗(yàn)結(jié)果分析

對(duì)淺埋爆炸荷載作用后的普通管道和CFRP管道變形進(jìn)行測(cè)量，如圖4所示。由圖4可知，由于爆炸超壓的作用，普通管道與CFRP管道的迎爆面均發(fā)生了較大的凹陷變形，CFRP布未發(fā)生破損。對(duì)兩組管道迎爆面的凹陷變形分別進(jìn)行了量測(cè)，其中，普通管道的變形量為120mm，CFRP管道變形量為70mm。

圖4 普通管道和CFRP管道變形量Fig.4 Deformation of plain pipe and CFRP pipe

2 數(shù)值模擬

2.1 本構(gòu)模型

采用LS-DYNA有限元軟件進(jìn)行數(shù)值模擬，有限元模型涉及五種材料：土體、TNT炸藥、管道、CFRP、空氣。

2.1.1 土體材料模型

試驗(yàn)場(chǎng)地土質(zhì)為黏性土，密度1.73×10-3g/cm3，剪切模量6.93×103MPa，體積模量2.1×104MPa，選用LS-DYNA的*MAT_SOIL_AND_FOAM材料模型對(duì)淺埋爆炸荷載下土體的動(dòng)態(tài)力學(xué)響應(yīng)進(jìn)行數(shù)值模擬，其理想塑性屈服函數(shù)為：

φ=J2-(a0+a1σ+a2σ2)

(1)

2.1.2 炸藥材料模型與狀態(tài)方程

采用*MAT_HIGH_EXPLOSIVE_BURN高能炸藥燃燒模型和JWL狀態(tài)方程來(lái)模擬TNT炸藥爆轟。炸藥密度為1.64×10-3g/cm3，爆速6.93×103m/s。JWL狀態(tài)方程對(duì)于爆炸過(guò)程中爆轟產(chǎn)物的壓力、體積和能量特性可以達(dá)到精準(zhǔn)控制，可表示為：

(2)

式中:p為爆轟產(chǎn)物的壓力，MPa；V為爆轟產(chǎn)物的相對(duì)體積；E0為爆轟產(chǎn)物單位體積初始內(nèi)能，J；w、A、B、R1、R2為材料常數(shù)。

2.1.3 管道材料模型

選用LS-DYNA材料庫(kù)中的*MAT_PLASTIC_ KINEMATIC材料模型來(lái)模擬鋼管道在爆炸荷載作用下的動(dòng)力特性。

2.1.4 CFRP材料模型

CFRP屬于單向受拉材料且在破壞前可以視為線彈性材料，因此選用LS-DYNA材料庫(kù)中的*MAT_ ELASTIC材料模型來(lái)模擬可以忽略應(yīng)變率效應(yīng)的CFRP材料。

2.1.5 空氣材料模型和狀態(tài)方程

選取LS_DYNA提供的*MAT_NULL材料模型和來(lái)模擬空氣，空氣密度設(shè)置為1.29 × 10-6g/cm3?？諝鉅顟B(tài)方程選用*EOS_LINEAR_POLYNOMIAL，可表示為:

P=C0+C1μ+C2μ2+C3μ3+(C4+C5μ+C6μ2)E1

(3)

(4)

圖5 有限元模型(尺寸單位：mm)Fig.5 Finite element model(size unit：mm)

式中:E1是單位初始體積的內(nèi)能，J；ρ是空氣質(zhì)量密度，g/cm3；ρ0是參考質(zhì)量密度，g/cm3；C0、C1、C2、C3、C4、C5和C6為狀態(tài)方程參數(shù)。

2.2 管土耦合有限元模型

炸藥在土中爆炸對(duì)埋地管道的響應(yīng)問(wèn)題是一個(gè)非線性問(wèn)題，利用解析計(jì)算方法研究埋地管道在淺埋爆炸作用下的破壞是非常困難的。因此，數(shù)值模擬更適合解決這個(gè)問(wèn)題。為了防止大變形和非線性結(jié)構(gòu)分析中的單元變形，本文采用了任意拉格朗日-歐拉(arbitrary Lagrange-Euler, ALE)算法建立了管-土三維有限元模型，有限元模型的示意圖如圖5所示。為了減少計(jì)算量，考慮到問(wèn)題對(duì)稱性，建立二分之一幾何結(jié)構(gòu)和荷載的分析模型。其中，管道外徑為73mm，管道壁厚為5.5mm。整個(gè)模型在X、Y、Z方向的尺寸分別為2000mm×1800mm×1500mm(長(zhǎng)×寬×高)。有限元網(wǎng)格大小為15mm，空氣、土壤、炸藥、管道使用Soild164單元，共劃分Solid164單元696087個(gè)， CFRP布使用Shell單元，共劃分Shell單元1480個(gè)。約束垂直于對(duì)稱平面(XY和YZ平面)的節(jié)點(diǎn)的過(guò)渡位移，土體單元四周及底部對(duì)比實(shí)際試驗(yàn)條件為透射邊界，頂部為自由邊界，以描述半無(wú)限土體。

2.3 數(shù)值模擬結(jié)果驗(yàn)證

數(shù)值模擬所獲取的管道超壓分布及變形如圖6所示，由圖6可知，兩類管道迎爆面均發(fā)生了凹陷變形。高應(yīng)力區(qū)出現(xiàn)在迎爆面正對(duì)爆心單元，然后向管道端部延伸衰減。由圖6(a)至圖6(d)可知，不同時(shí)刻兩組管道受到爆炸超壓后的變形情況有所不同，炸藥爆炸1ms兩組管道開(kāi)始受到爆炸超壓影響，至2ms高應(yīng)力區(qū)出現(xiàn)在迎爆面正對(duì)爆心單元，然后向管道端部延伸衰減，迎爆面發(fā)生較為明顯塑性變形，至3ms迎爆面塑性變形進(jìn)一步加大，背爆面開(kāi)始變形，管道中部開(kāi)始呈現(xiàn)微小彎曲變形，最后在5ms時(shí)管道背爆面發(fā)生明顯塑性變形，管道中部彎曲進(jìn)一步加大。圖7是普通管道與CFRP管道迎爆面凹陷時(shí)程曲線圖，由圖7可知，首先兩組管道迎爆面變形起始均未在零時(shí)刻，表明管道周邊土體使得爆炸超壓作用在管道上的時(shí)間得到了延緩，同時(shí)，CFRP管道的變形開(kāi)始時(shí)間略晚于普通管道，表明CFRP的存在一定程度上延緩了管道變形的發(fā)生。普通管道迎爆面凹陷變形數(shù)值模擬結(jié)果為125.98mm，試驗(yàn)結(jié)果為120mm，兩者誤差為4.9%，CFRP管道迎爆面凹陷變形數(shù)值模擬結(jié)果為76.43mm，試驗(yàn)結(jié)果為70mm，兩者誤差為9.1%。兩組管道凹陷變形的數(shù)值模擬結(jié)果實(shí)際誤差均未超出10%。綜上所述，數(shù)值模擬結(jié)果是合理的。

圖6 不同時(shí)刻淺埋爆炸荷載下兩類管道超壓分布及變形Fig.6 Overpressure distribution and deformation of two kinds of pipelines under shallow explosion load at different time

圖8為土體變形圖，由圖8可知，看出在爆炸荷載作用下土體發(fā)生變形，并逐漸向兩邊拋擲，形成爆坑。為了研究爆炸沖擊波超壓在土中規(guī)律，選取沿深度及距離爆心位置折合距離為0.38、0.41、0.44、0.49m/kg1/3的四個(gè)測(cè)點(diǎn)，對(duì)土中超壓進(jìn)行觀測(cè)。繪制土中爆炸產(chǎn)生的沖擊波壓力-時(shí)間曲線如圖9所示，由圖9可知，各個(gè)測(cè)點(diǎn)的爆炸波壓力隨著時(shí)間的增進(jìn)而快速衰減，最終趨近于零，且不同測(cè)點(diǎn)的爆炸波峰值壓力也存在著一定差異性，埋置深度和折合距離的增加，沖擊波的峰值壓強(qiáng)相應(yīng)迅速降低。

3 結(jié)果與討論

3.1 CFRP層數(shù)對(duì)管道的影響

圖7 普通管道與CFRP管道迎爆面凹陷時(shí)程曲線圖Fig.7 Time-course diagram of burst-front depression for plain and CFRP pipes

為了在證實(shí)CFRP布提升管道抗爆性能的基礎(chǔ)上，進(jìn)一步探索不同CFRP層數(shù)對(duì)管道抗爆性能的影響，建立三種不同的CFRP層數(shù)的管道在同等淺埋爆炸荷載下的數(shù)值模型，并獲取管道超壓分布及變形如圖10所示，由圖10可知，CFRP層數(shù)的增多有效地提高了防爆性能?；贏LE數(shù)值模擬方法，獲取三種不同層數(shù)的CFRP包裹條件下管道的最大等效應(yīng)變變化情況，結(jié)果如圖11所示。由圖11可知，隨著CFRP層數(shù)的增多，管道受爆炸荷載后，變形量減小。同時(shí)，用6層CFRP包裹加固的管道的等效應(yīng)變要比1層CFRP包裹加固的管道的等效應(yīng)變要小28.3%。用10層CFRP包裹加固的管道的等效應(yīng)變要比6層CFRP包裹加固的管道的等效應(yīng)變要小56.4%，可見(jiàn)在一定程度上增加CFRP層數(shù)可以有效提升管道抗爆性能。

圖8 土體變形圖 圖9 土中不同測(cè)點(diǎn)壓力-時(shí)間曲線Fig.8 Soil deformation diagram Fig.9 Pressure-time curves at different measurement points in the soil

3.2 炸藥量的影響

進(jìn)一步討論不同炸藥量條件(5、10、20、35kg)下，不同CFRP包裹層數(shù)的管道變形情況如圖12所示。由圖12可知，隨著TNT量級(jí)的增加，不同CFRP包裹的管道變形量均相應(yīng)增加。當(dāng)炸藥量小于5kg時(shí)，纏繞1層CFRP布的管道斜率較大，而纏繞6層和10層CFRP布的管道斜率變化較小。分析其原因，主要是因?yàn)橥馏w吸收了部分爆炸能量，傳遞到管道上的超壓較小，同時(shí)包裹的CFRP也起到了一定的緩沖作用。

3.3 爆心距的影響

爆心距是影響管道變形的一個(gè)重要因素?；跀?shù)值模擬結(jié)果，獲取不同CFRP包裹管道在不同爆心距條件下的迎爆面最大等效應(yīng)力變化情況如圖13所示。由圖13可知，隨著爆心距的增加，最大等效應(yīng)力逐漸下降。同時(shí)，而在同一爆心距條件下，CFRP層數(shù)越多，相應(yīng)的最大等效應(yīng)力越小。在最大等效應(yīng)力衰減方面，爆心距在0.3～0.6m時(shí)，最大等效應(yīng)力衰減較緩慢，大于0.6m時(shí)衰減較快。因此，小于0.6m的爆心距條件可以視為淺埋爆炸荷載下影響管道變形的重要因素之一。

圖10 淺埋爆炸荷載下不同層數(shù)CFRP管道超壓分布情況Fig.10 Overpressure distribution of CFRP pipes of different layers under shallow buried blast load

4 結(jié)論

圖11 不同CFRP層數(shù)管道最大等效應(yīng)變變化Fig.11 Variation of maximum equivalent variation for pipes with different number of CFRP layers

通過(guò)設(shè)計(jì)開(kāi)展淺埋爆炸試驗(yàn)和數(shù)值模擬，研究了爆炸荷載下埋地管道的動(dòng)力響應(yīng)特征，并探討了不同CFRP層數(shù)對(duì)管道抗爆性能的影響，所得結(jié)論如下：

1)基于淺埋爆炸試驗(yàn)結(jié)果，驗(yàn)證了ALE方法在模擬土中淺埋爆炸荷載與管道相互作用的有效性和可行性，該方法可以很好地描述爆炸波在土中的傳播規(guī)律，以及管-土間的耦合作用。

2)CFRP包裹可以有效提升油氣管道的抗爆性能。從同等淺埋爆炸荷載下，不同包裹層數(shù)的CFRP管道最大等效應(yīng)變的變化可知，在一定范圍內(nèi)增加CFRP層數(shù)可以有效提升埋地管道的抗爆性能。

3)隨著炸藥量的提升，管道迎爆面的凹陷變形量相應(yīng)提升。同時(shí)，5kg炸藥量對(duì)于超過(guò)1層CFRP包裹的管道迎爆面變形量降低有限，說(shuō)明在小于5kg藥量條件下，如果對(duì)管道抗爆無(wú)特殊要求，可以適當(dāng)減少CFRP包裹的層數(shù)。

4)隨著爆心距的增加，管道迎爆面最大等效應(yīng)力逐漸下降。爆心距大于0.6m時(shí)，管道最大有效應(yīng)力衰減較快。因此，小于0.6m的爆心距條件會(huì)對(duì)管道安全形成較大的威脅，需要盡量規(guī)避。

圖12 不同CFRP層數(shù)炸藥量與管道凹形 圖13 不同CFRP層數(shù)爆心距與迎爆面最大 變形量的關(guān)系 等效應(yīng)力的關(guān)系Fig.12 Relationship between the amount of explosive charge Fig.13 Relationship between the blast centre distance and the amount of concave deformation of and the maximum equivalent force on the face the pipe for different CFRP layers of the blast for different CFRP layers