徐 偉，李世坤，王 壬，徐建國(guó)

(1.鄭州大學(xué)水利科學(xué)與工程學(xué)院，河南 鄭州 450001；2.中國(guó)人民解放軍陸軍步兵學(xué)院 石家莊校區(qū)，河北 石家莊 050200)

排水管道保證了城市居民生活廢水、工業(yè)廢水以及雨水的正常排放，是衡量現(xiàn)代化城市水平的重要標(biāo)志[1]。近年來(lái)，我國(guó)排水管道建設(shè)發(fā)展迅速，截至2015年全國(guó)大中城市排水管道的總長(zhǎng)已超過(guò)54萬(wàn)km[2]。在交通荷載作用下排水管道災(zāi)害事故日益頻發(fā)，造成其滲漏、脫空和斷裂等情況[3]。因此，研究排水管道在交通荷載下的脫空修復(fù)和動(dòng)力特性十分必要。范新卓等[4]建立了路-管-土系統(tǒng)結(jié)構(gòu)的有限元模型，分析車輛多輪荷載作用下埋地管道的動(dòng)應(yīng)力響應(yīng)特性。吳小剛等[5]建立了交通荷載作用下管道的Eular-Bernoulli彈性地基梁受力模型。徐建國(guó)等[6]采用ABAQUS有限元模擬分析了地下管道經(jīng)高聚物修復(fù)前后在靜力荷載及地震荷載作用下的受力與變形特性。王直民等[7]，將交通荷載簡(jiǎn)化為不連續(xù)的半波正弦荷載，采用有限元法分析了單次和多次交通荷載作用下埋地管道的動(dòng)力響應(yīng)規(guī)律。王復(fù)明等[8]建立了道路結(jié)構(gòu)與帶承插口的埋地排水管相互作用的三維數(shù)值模型，分析了在沖擊荷載作用下管道的動(dòng)力響應(yīng)。蘭國(guó)冠等[9]采用ANSYS軟件分析了穿越公路的輸氣管道在單輪車輛載荷作用下的力學(xué)性狀，探討了橫穿道路管道受車輛荷載作用時(shí)管道中點(diǎn)應(yīng)力及豎向位移隨時(shí)間變化規(guī)律。周正峰等[10]應(yīng)用ABAQUS有限元軟件，考慮管-土相互作用，分析了管周附加應(yīng)力的分布特征，計(jì)算了飛機(jī)、施工用重型車輛雙軸雙輪荷載及壓路機(jī)荷載作用下，管道附加應(yīng)力及其引起的管道結(jié)構(gòu)應(yīng)力、變形隨管道埋深的變化規(guī)律。FANG Hongyuan 等[11]基于有限元探究了排水管道在高聚物注漿修復(fù)前后的力學(xué)特性，分析了排水管道在荷載作用下不同工況的力學(xué)響應(yīng)。但大部分?jǐn)?shù)值分析的結(jié)果未得到現(xiàn)場(chǎng)試驗(yàn)驗(yàn)證?；诖?，筆者以車輛多輪動(dòng)載對(duì)地下排水管道的動(dòng)力特性研究為基礎(chǔ)，采用 ABAQUS有限元分析軟件，建立了帶有承插結(jié)構(gòu)的排水管道三維數(shù)值模型，從數(shù)值分析和現(xiàn)場(chǎng)試驗(yàn)方面探討交通荷載作用下埋地管道的動(dòng)力反應(yīng)，研究成果將為管道注漿技術(shù)、數(shù)值建模分析、工程搶險(xiǎn)及高聚物注漿推廣應(yīng)用等提供相關(guān)指導(dǎo)。

1 地下排水管道足尺試驗(yàn)

1.1 試驗(yàn)設(shè)計(jì)

試驗(yàn)采用小型設(shè)備開挖長(zhǎng)寬深為8 m×2.3 m×2.3 m的試驗(yàn)坑進(jìn)行測(cè)試，試驗(yàn)選取3節(jié)完整管段和兩節(jié)半截管段通過(guò)承插口連接，有效長(zhǎng)度為8 m，管體埋深1.1 m(見圖1)。圖中Pn為第n節(jié)管道，jm為第m個(gè)承插口。

圖1 試驗(yàn)管道Fig.1 Test pipeline

為了評(píng)價(jià)高聚物注漿技術(shù)應(yīng)用于脫空管道的修復(fù)效果，試驗(yàn)進(jìn)行了管道正常埋置，管道底部貫穿脫空和高聚物注漿修復(fù)3種不同工況下的試驗(yàn)。脫空區(qū)設(shè)計(jì)為沿管底(管臀以下約1/4圓環(huán)區(qū)域)縱向底部貫穿式空洞，其截面寬度和高度為0.5 m×0.15 m，車輛各軸左側(cè)車輪沿著管線縱向以10 km/h的速度勻速通過(guò)試驗(yàn)管道。表1為試驗(yàn)用車整車參數(shù)。試驗(yàn)管道應(yīng)變片布置如圖2所示。

表1 試驗(yàn)用車整車參數(shù)Table 1 Parameters of test vehicle

圖2 試驗(yàn)管道應(yīng)變片布置Fig.2 Arrangement of pipeline strain gauge

1.2 試驗(yàn)管道參數(shù)

試驗(yàn)管選取的是柔性接頭C型鋼承口RCPⅢ700×2000GB/T11836，即管段的內(nèi)徑為0.7 m，壁厚為0.07 m，長(zhǎng)度為2.0 m。管道材料為C30混凝土管，抗壓強(qiáng)度為30 MPa，管內(nèi)橡膠圈起到密封防滲的作用，符合標(biāo)準(zhǔn)《混凝土和鋼筋混凝土排水管》(GB/T 11836—2009)中的規(guī)定要求 。

1.3 測(cè)點(diǎn)布置

試驗(yàn)采用80 mm和30 mm長(zhǎng)的單軸測(cè)試應(yīng)變片，通過(guò)粘貼在管道外表面得到鋼筋混凝土管結(jié)構(gòu)的力學(xué)響應(yīng)數(shù)據(jù)。其中較長(zhǎng)的應(yīng)變片沿管線縱向粘貼，較短的沿圓周環(huán)向粘貼[12]，在管身3個(gè)截面以及鋼承口截面，應(yīng)變測(cè)點(diǎn)沿圓周方向布置，管圓周方向0°、90°和180°分別代表管頂、管西測(cè)和管底。數(shù)據(jù)采集采用東華測(cè)試技術(shù)公司生產(chǎn)的DH3821靜態(tài)及動(dòng)態(tài)應(yīng)變測(cè)試分析系統(tǒng)，通過(guò)網(wǎng)線與電腦相連，對(duì)參數(shù)、采樣頻率進(jìn)行設(shè)置，每次測(cè)試開始前，將所有設(shè)備重置歸零。

2 地下排水管道數(shù)值分析模型

2.1 管道有限元模型建立

利用ABAQUS軟件對(duì)管-土-高聚物模型進(jìn)行了三維有限元模擬，三維數(shù)值模型長(zhǎng)寬高為8 m×10 m×10 m，土體采用Mohr-Coulomb彈塑性本構(gòu)，橡膠圈變能函數(shù)采用文獻(xiàn)[13]Mooney-Rivlin相關(guān)參數(shù)，模型材料參數(shù)如表2所示。為了確保數(shù)值模型和現(xiàn)場(chǎng)試驗(yàn)結(jié)構(gòu)間的一致性，垂直邊界設(shè)置為橫向約束，底部邊界設(shè)置為固定約束，考慮管道與土體自重。對(duì)管土模型進(jìn)行網(wǎng)格劃分時(shí)，采用C3D8R六面體單元，整體結(jié)構(gòu)網(wǎng)格和承插口模型如圖3和圖4所示。

表2 材料參數(shù) Table 2 Material parameters

圖3 整體結(jié)構(gòu)網(wǎng)格圖Fig.3 Grid diagram of overall structure

圖4 承插口細(xì)部圖Fig.4 Detailed drawing of socket

2.2 接觸面模型

管周土體不僅是作用在管道上的荷載，還是增強(qiáng)管道強(qiáng)度和剛度的一種介質(zhì)[14]。在管土交界面設(shè)置接觸單元，模擬不同介質(zhì)間產(chǎn)生的相對(duì)錯(cuò)動(dòng)與分離狀態(tài)。切向采用罰函數(shù)，法向采用硬接觸，設(shè)Fs和Fn分別為接觸單元間的摩擦力和法向力。

(1)

(2)

式中:Kt為黏性系數(shù)；u為切向位移；F為接觸單元為發(fā)生滑移時(shí)的摩擦力；Kn為法向剛度，d為法向接觸距離。

2.3 荷載模型

基于ABAQUS用戶子程序平臺(tái)，使用 Fortran 語(yǔ)言編寫子程序Dload，實(shí)現(xiàn)兩個(gè)車軸兩側(cè)車輪平行豎向的移動(dòng)加載過(guò)程，用來(lái)模擬真實(shí)的車輛輪載。需在子程序中定義車輛各軸車輪的起始位置、車輪荷載作用位置、車輛前進(jìn)的速度、車輪荷載大小和作用尺寸等。輪胎接觸面等效為矩形接觸面[15]，前軸單輪輪胎接地面積為0.018 m2，后軸雙輪胎接地面積為0.036 m2，輪載大小為在接觸面內(nèi)的均布荷載。圖5為荷載移動(dòng)過(guò)程，可以看到在四個(gè)輪載作用下土體表面位移的沉降。

圖5 荷載移動(dòng)過(guò)程Fig.5 Load moving process

3 結(jié)果對(duì)比分析

3.1 3種工況下試驗(yàn)值與數(shù)值結(jié)果比較

為了驗(yàn)證模型的有效性，選取J3鋼承口截面底部應(yīng)變數(shù)據(jù)在正常埋置、底部貫穿脫空和高聚物注漿修復(fù)3種不同工況下與數(shù)值模擬結(jié)果進(jìn)行比較(見圖6)。

圖6 試驗(yàn)結(jié)果與數(shù)值結(jié)果對(duì)比Fig.6 Comparison of experimental and numerical results

在車輛多輪動(dòng)載下，應(yīng)變時(shí)程曲線形狀均類似于齒根，呈受壓狀態(tài)，在出現(xiàn)兩次峰值后逐漸趨于平緩，兩個(gè)峰值大小不同，后者大于前者，這是由于軸荷較小的前輪先經(jīng)過(guò)測(cè)點(diǎn)，所以后輪通過(guò)時(shí)產(chǎn)生的應(yīng)變大于前輪。車輛軸距為3.815 m，以10 km/h的速度勻速通過(guò)J3采集點(diǎn)，兩齒根峰值的理論時(shí)間間距約為1.4 s，與試驗(yàn)值相符。對(duì)比結(jié)果如表3所示，可以看出模擬結(jié)果略大于試驗(yàn)值，在同一數(shù)量級(jí)內(nèi)，模擬結(jié)果與試驗(yàn)值吻合良好，說(shuō)明所采用方法的正確性和有效性，模型網(wǎng)格劃分合理，計(jì)算準(zhǔn)確可靠，可以用于地下排水管道結(jié)構(gòu)分析。

表3 試驗(yàn)結(jié)果與數(shù)值結(jié)果對(duì)比Table 3 Comparison of experimental and numerical results

3.2 3種工況下承插口處的應(yīng)力分布

由于承插口往往是引起管道破壞的源頭之一[16]，對(duì)管道的整體性能起著非常重要的作用。當(dāng)車輛后輪經(jīng)過(guò)采集點(diǎn)時(shí)，該點(diǎn)的應(yīng)變值達(dá)到最大，為最不利狀態(tài)，因此將分析對(duì)象選為后輪經(jīng)過(guò)J3采集點(diǎn)時(shí)的J3承插口。取J3承插口環(huán)向路徑，對(duì)3種不同工況下承插口Mises應(yīng)力進(jìn)行分析，其應(yīng)力對(duì)比曲線如圖7所示。0°、90°和180°分別代表管道的頂部、側(cè)面和底部，與試驗(yàn)測(cè)點(diǎn)布置環(huán)向相同。

由于管線結(jié)構(gòu)的對(duì)稱性[17]，管道兩側(cè)對(duì)稱部位的Mises應(yīng)力值較一致。從圖7可以看出，管道在正常埋置和高聚物注漿修復(fù)后的承插口Mises應(yīng)力曲線較吻合，且修復(fù)后的Mises應(yīng)力值較正常埋置值小。承口處底部Mises應(yīng)力值最大，略大于頂部，在管肩和管臀處Mises應(yīng)力值最小，兩種工況之間承插口環(huán)向Mises應(yīng)力差值在管底處最大(差值分別為0.60 MPa和0.05 MPa)，管側(cè)和管頂次之，管肩和管臀最小。可見，管道在高聚物修復(fù)之后可以達(dá)到正常埋置時(shí)的受力狀態(tài)，效果較好。

圖7 不同工況作用下J3承插口環(huán)向Mises應(yīng)力對(duì)比曲線Fig.7 Comparison curve of circumferential Mises stress of J3 socket under different working conditions

當(dāng)管道自管臀以下土體脫空時(shí)，承插口環(huán)向Mises應(yīng)力狀態(tài)發(fā)生很大的變化，相較于正常埋置，承口兩側(cè)管臀的Mises應(yīng)力值分別增大了2.08倍和1.97倍，管底壓應(yīng)力減小很多，插口在管臀和管底處的Mises均減小很多。這是因?yàn)楣艿涝臼苌舷峦馏w約束，在行車荷載作用下產(chǎn)生受壓變形，由于底部土體的脫空致管底失去了土體的約束，使得整個(gè)管道截面受力整體減小，尤其在底部更為明顯，而管臀附近處的土體承擔(dān)大部分的支撐作用，所以在管臀處會(huì)產(chǎn)生較大的應(yīng)力集中，該處的Mises應(yīng)力值自然變大。因此，管道底部脫空，對(duì)其整體受力有很大的影響。

3.3 不同埋深排水管道動(dòng)力響應(yīng)分析

為了研究土體覆蓋深度對(duì)管線力學(xué)性能的影響，進(jìn)行了覆蓋深度分別為0.9 m、1.0 m、1.1 m和1.2 m的數(shù)值模擬。在車載為6.1t作用下，車輛單側(cè)輪壓沿著管線縱向以10 km/h的速度勻速通過(guò)試驗(yàn)管道上方覆土，將J3鋼承口截面底部應(yīng)變時(shí)程最大值、P3插口內(nèi)頂部位移時(shí)程最大值和后輪經(jīng)過(guò)J3采集點(diǎn)時(shí)的J3承插口環(huán)向Mises應(yīng)力最大值進(jìn)行比較(見表4)。

表4 不同埋深作用下J3承插口最大應(yīng)變、位移和環(huán)向Mises應(yīng)力Table 4 Maximum strain,displacement and circumferential Mises stress of J3 bell and spigot under different buried depth

不同埋深作用下J3鋼承口截面底部應(yīng)變和P3插口內(nèi)頂部位移對(duì)比曲線如圖8所示。

圖8 不同埋深作用下底部應(yīng)變和頂部位移對(duì)比曲線Fig.8 Comparison curve between bottom strain and top displacement under different buried depth

可以看出位移與埋深成反比，當(dāng)埋深從0.9 m分別增加到1.0 m、1.1 m和1.2 m時(shí)，J3鋼承口截面底部應(yīng)變最大值分別減小了0.266 MPa、0.212 MPa和0.088 MPa,P3插口內(nèi)頂部位移最大值分別減小了0.005 mm、0.004 mm和0.003 mm。說(shuō)明埋深增量相同時(shí)，隨埋深的增加，J3鋼承口截面底部應(yīng)變和P3插口內(nèi)頂部位移最大值的減少量在減小。

當(dāng)車輛后輪經(jīng)過(guò)J3采集點(diǎn)時(shí)，J3承插口的環(huán)向Mises應(yīng)力如圖9所示。Mises應(yīng)力最大值與埋深成正比，承口部位的環(huán)向Mises應(yīng)力要遠(yuǎn)大于插口部位，這是由于兩個(gè)部位材料屬性的不同，管道部位的Mises應(yīng)力同時(shí)受到土體和管道自重及交通荷載的影響。

圖9 不同埋深作用下J3承插口環(huán)向Mises應(yīng)力對(duì)比曲線Fig.9 Comparison curve of circumferential Mises stress of J3 bell and spigot under different buried depth

圖10為除去了土壓力，即車輛荷載對(duì)J3承插口的環(huán)向Mises應(yīng)力的附加影響?？梢钥闯?，埋深為0.9 m時(shí)交通荷載產(chǎn)生的Mises應(yīng)力增量較大，埋深為1.1 m和1.2 m這兩種情況產(chǎn)生的Mises應(yīng)力增量較小且較為一致，說(shuō)明隨著埋深的增加，交通荷載對(duì)管道施加的應(yīng)力逐漸減小，即土壓力呈線性增加，成為管道的主要荷載，相較于土壓力引起的管道受力變化，交通荷載不再起主導(dǎo)作用[18]。

圖10 車輛荷載對(duì)J3承插口環(huán)向Mises應(yīng)力附加影響的對(duì)比曲線Fig.10 Effect of vehicle load on circumferential Mises stress of J3 bell and spigot

3.4 不同車載排水管道動(dòng)力響應(yīng)分析

在管道埋深為1.1 m，行車速度和車載作用位置不變，計(jì)算對(duì)比4種不同的車輛荷載對(duì)管道力學(xué)響應(yīng)的影響(見表5)。

表5 不同車載作用下管道J3承插口最大應(yīng)變、位移和環(huán)向Mises應(yīng)力Table 5 Maximum strain,displacement and circumferential Mises stress of J3 bell and spigot under different vehicle loads

圖11為J3鋼承口截面底部應(yīng)變時(shí)程和P3插口內(nèi)頂部位移時(shí)程。且隨著車載的增大，J3鋼承口截面底部應(yīng)變和P3插口內(nèi)頂部位移也隨之增加。尤其在車輛后輪經(jīng)過(guò)J3接口時(shí)，車載為6.1 t、8.1 t和10.1 時(shí)的J3鋼承口截面底部應(yīng)變分別是4.1 t時(shí)的1.57、1.98和2.22倍，P3插口內(nèi)頂部位移分別是4.1 t時(shí)的1.54、1.99和2.34倍。

圖11 不同車載作用下底部應(yīng)變和頂部位移對(duì)比曲線Fig.11 Comparison curve between bottom strain and top displacement under different vehicle loads

不同車載作用下J3承插口環(huán)向Mises應(yīng)力對(duì)比曲線如圖12所示，因?yàn)檐囕v載重不一樣，輪胎接觸地面時(shí)壓力的分布也會(huì)有所差異，J3承插口的環(huán)向Mises應(yīng)力差別明顯，所以地面交通荷載越大，對(duì)管道的受力狀態(tài)越不利。

圖12 不同車載作用下J3承插口環(huán)向Mises應(yīng)力對(duì)比曲線Fig.12 Comparison curve of circumferential Mises stress of J3 bell and spigot under different vehicle loads

3.5 不同車速排水管道動(dòng)力響應(yīng)分析

在管道埋深為1.1 m，車載大小和作用位置不變，改變車速為 10 km/h、40 km/h、70 km/h和100 km/h，管道J3鋼承口截面底部應(yīng)變、P3插口內(nèi)頂部位移和環(huán)向Mises應(yīng)力如表6所示。

表6 不同車速作用下管道J3承插口最大應(yīng)變、位移和環(huán)向Mises應(yīng)力Table 6 Maximum strain,displacement and circumferential Mises stress of J3 bell and spigot under different vehicle speeds

不同車速作用下,底部應(yīng)變和頂部位移對(duì)比曲線如圖13所示,J3承插口環(huán)向Mises應(yīng)力對(duì)比曲線圖14所示。當(dāng)車速由10 km/h增加到100 km/h時(shí)，各項(xiàng)指標(biāo)均有所減小，以J3鋼承口截面底部應(yīng)變?yōu)槔S著移動(dòng)車速的增加，其值不斷減小，減幅分別為5.4%、4.7%和 5.2%。這是因?yàn)樾旭偟钠噷?duì)地面施加的荷載有瞬時(shí)性[19]，隨著車輛速度的提高，使輪胎與地面的接觸時(shí)間變短。由于地下結(jié)構(gòu)中應(yīng)力的傳遞是由相鄰顆粒來(lái)完成的，若應(yīng)力出現(xiàn)的時(shí)間很短，則來(lái)不及傳遞分布到地面及管道結(jié)構(gòu)中去，其動(dòng)力特性便不能像速度較慢時(shí)呈現(xiàn)的那樣完全，從而降低了對(duì)埋地管道的影響。

圖13 不同車速作用下底部應(yīng)變和頂部位移對(duì)比曲線Fig.13 Comparison curve between bottom strain and top displacement under different vehicle speeds

圖14 不同車速作用下J3承插口環(huán)向Mises應(yīng)力對(duì)比曲線Fig.14 Comparison curve of circumferential Mises stress of J3 bell and spigot under different vehicle speeds

4 結(jié) 論

(1)數(shù)值模擬結(jié)果略大于試驗(yàn)值，兩者符合較好，說(shuō)明建立的ABAQUS數(shù)值分析模型模擬車輛多輪動(dòng)載的分析方法具有正確性和有效性，可進(jìn)一步開展地下排水管道結(jié)構(gòu)的數(shù)值分析。

(2)管道底部脫空對(duì)管道受力影響較大，高聚物注漿對(duì)底部脫空管道起到了較好的修復(fù)作用，并能夠?qū)⑵涫芰顟B(tài)恢復(fù)到正常埋置狀態(tài)。

(3)隨著排水管道埋深的增加，交通荷載對(duì)管道施加的應(yīng)力逐漸減小，而土壓力呈線性增加，成為管道的主要承重荷載，相較于土壓力引起的管道受力變化，交通荷載不再起主導(dǎo)作用。

(4)荷載大小對(duì)埋地管道的受力特性影響較大，交通荷載越大對(duì)管道越不利。

(5)交通荷載對(duì)管道產(chǎn)生的應(yīng)變和位移影響隨著移動(dòng)車速的增加呈減小趨勢(shì)，說(shuō)明交通荷載的移動(dòng)車速對(duì)地下排水管道的動(dòng)力響應(yīng)具有時(shí)間效應(yīng)。