張 峰,裴華富
(大連理工大學(xué)巖土工程研究所,遼寧 大連 116024)
滑坡是一種嚴(yán)重的全球性自然災(zāi)害,給人們的生命財(cái)產(chǎn)造成了巨大損失。邊坡監(jiān)測是滑坡災(zāi)害防治研究的重要組成部分[1-2]。根據(jù)地表位移、地下水壓力、樁身位移、錨固應(yīng)力、深部位移等監(jiān)測對象的不同,大致可分為多種監(jiān)測技術(shù)[3-5]。深部位移監(jiān)測是其中必不可少的內(nèi)容,可以獲取邊坡任意深度的水平位移,確定滑裂面和薄弱帶的位置,有效判斷邊坡的穩(wěn)定性。此外,該方法可用于監(jiān)測所有滑坡和基坑坍塌變形。
測斜儀是監(jiān)測邊坡深部位移的主要設(shè)備。傳統(tǒng)的測斜儀包括振弦式、伺服加速度計(jì)、應(yīng)變式和電子式[6-7]。然而,這些傳統(tǒng)的測斜儀存在穩(wěn)定性差、電磁干擾大等缺點(diǎn),并且大多屬于點(diǎn)式監(jiān)測,采集的數(shù)據(jù)量有限,難以滿足邊坡監(jiān)測的要求。在過去的幾十年里,分布式光纖傳感(DFOS)技術(shù)得到了發(fā)展,包括光纖布拉格光柵(FBG)、光時(shí)域反射(OTDR)、布里淵光時(shí)域反射/分析(BOTDR/A)和布里淵光頻域反射計(jì)/分析[8-11]。Ho[12]等應(yīng)用時(shí)域反射(TDR)技術(shù)檢測九峰山滑坡位移并測量其變形大小。該技術(shù)可有效監(jiān)測滑坡滑動面的位置,邊坡位移可以通過TDR反射波形的變化來確定。裴華富[3]等人研制了一種基于FBG的新型原位測斜儀。測斜儀的位移由光纖光柵傳感器測得的應(yīng)變經(jīng)過有限差分法計(jì)算得到。Sun[8]等人開發(fā)了基于歐拉梁理論的BOTDR測斜儀,并將其安裝在實(shí)際邊坡上進(jìn)行滑坡監(jiān)測。然而這些技術(shù)在提高空間分辨率、監(jiān)測精度、應(yīng)變位移轉(zhuǎn)換算法等方面仍存在一定的局限性。例如,TDR技術(shù)空間分辨率低,誤差大。FBG屬于準(zhǔn)分布式技術(shù),測點(diǎn)間距離較大,難以準(zhǔn)確捕捉位移突變點(diǎn)處滑移面的位置。BOTDR技術(shù)雖然可以實(shí)現(xiàn)分布式測量,但其監(jiān)測精度低于FBG技術(shù)。OFDR技術(shù)結(jié)合了分布式監(jiān)測和高精度的優(yōu)點(diǎn),與傳統(tǒng)的傳感設(shè)備相比,具有明顯的優(yōu)勢[13-14],主要體現(xiàn)在以下幾個(gè)方面:1)傳感器的選擇更簡單,更廣泛。溫度和應(yīng)變測量可以用普通的通信光纖進(jìn)行,既便宜又簡單;2)工作原理先進(jìn),空間分辨率可達(dá)1 mm,測量精度可達(dá)1 με;3)使用的光纜為普通光纜,價(jià)格便宜,單根光纖可同時(shí)測量數(shù)千個(gè)傳感點(diǎn),測量精度高;4)分布式測量,傳感器結(jié)構(gòu)簡單,測試性能穩(wěn)定。
因此,本研究將OFDR技術(shù)應(yīng)用于邊坡監(jiān)測。首先,研制了一種基于共軛梁算法的OFDR原位測斜儀。然后,通過數(shù)值模擬驗(yàn)證了共軛梁算法的有效性,通過對比實(shí)驗(yàn)驗(yàn)證了OFDR原位測斜儀的精度,通過模型實(shí)驗(yàn)驗(yàn)證了測斜儀的實(shí)用性。最后,對這些結(jié)果進(jìn)行了討論和分析,以驗(yàn)證基于共軛梁算法的OFDR原位測斜儀可以有效地用于邊坡位移監(jiān)測中。
當(dāng)被測光纖置于溫度場或應(yīng)變場中時(shí),光纖內(nèi)部的折射率分布會在溫度或應(yīng)變的影響下發(fā)生變化,瑞利散射信號光的相應(yīng)頻率也會發(fā)生變化。通過測量瑞利散射信號光的頻率可以得到對應(yīng)外部溫度場或應(yīng)變場的變化,從而實(shí)現(xiàn)分布式光纖傳感測量。
OFDR的基本原理如圖1所示。來自光源的線性掃描連續(xù)光被耦合器分成兩路,其中一路作為參考光。另一個(gè)用作檢測光并發(fā)射到要測量的光纖中。當(dāng)檢測光在光纖中向前傳播時(shí),會連續(xù)產(chǎn)生瑞利散射信號。這些信號光和反射參考光通過耦合器發(fā)生拍頻干涉,被光電檢測器檢測到。將參考光和檢測光得到的瑞利散射信號分成若干個(gè)信號窗口,通過互相關(guān)運(yùn)算計(jì)算每個(gè)信號窗口的光譜偏移量。譜位移與應(yīng)變和溫度變化的關(guān)系可用下式表示[13-14]:
圖1 OFDR技術(shù)示意圖(改自文獻(xiàn)[13])
其中Cε和CT是校準(zhǔn)常數(shù)。
基于上述OFDR技術(shù),設(shè)計(jì)了測斜儀。將一根橫截面直徑為20 mm,長度為1 000 mm的PVC棒上刻出一個(gè)深度為2 mm的凹槽,將光纖預(yù)先拉直放入凹槽中,然后用環(huán)氧樹脂膠進(jìn)行粘貼制作測斜儀,如圖2所示。
圖2 OFDR原位測斜儀設(shè)計(jì)(單位:mm)
環(huán)氧樹脂可以很好地將光纖和PVC棒粘合在一起,使光纖精確測量PVC棒表面的應(yīng)變分布。OFDR技術(shù)可每1 mm取一個(gè)測點(diǎn),獲得一個(gè)應(yīng)變值。對滑坡進(jìn)行監(jiān)測時(shí),預(yù)先在邊坡相應(yīng)位置鉆孔,將OFDR原位測斜儀垂直插入鉆孔內(nèi),如圖3(a)所示。當(dāng)邊坡在地震、降雨或外部荷載作用下產(chǎn)生滑動時(shí),由于PVC棒具有良好的彈性,能隨土體協(xié)調(diào)變形。根據(jù)PVC棒表面光纖測得的應(yīng)變,結(jié)合相應(yīng)算法可計(jì)算出OFDR原位測斜儀的位移,從而反映邊坡的變形,判斷邊坡的穩(wěn)定狀態(tài)。
圖3 位移計(jì)算示意圖
測斜儀垂直插入邊坡中,可視為懸臂梁,下端固定,上端自由。假設(shè)測斜儀的水平位移相對于測斜儀的長度較小,即滿足材料力學(xué)中的小變形假設(shè),忽略剪切變形的影響,則測斜儀可視為小變形歐拉梁,滿足曲率k與彎矩M的物理關(guān)系:
式中:E——材料的彈性模量;
I——截面的慣性矩。
在共軛梁理論中[5],實(shí)梁的撓度和角度對應(yīng)于共軛梁的彎矩和剪力。實(shí)梁的曲率在數(shù)值上等于共軛梁的荷載集度?;谠摾碚摚梢詫?shí)現(xiàn)懸臂梁的撓度計(jì)算。結(jié)合材料力學(xué),可以建立如下關(guān)系:
式中: ε(x)——光纖測得的測斜儀各點(diǎn)的應(yīng)變;
y——彈性梁的半徑;
假設(shè)測斜儀的長度為L,按照OFDR技術(shù)的采樣間隔,L被分為n個(gè)單元,則每個(gè)分段的長度為l=L/n。各單元上下表面分別取兩個(gè)應(yīng)變監(jiān)測點(diǎn),簡化模型如圖3(b)所示。各截面載荷分布可表示為:
根據(jù)共軛梁的荷載分布,可以計(jì)算出共軛梁上任意一點(diǎn)的彎矩。
根據(jù)實(shí)梁與共軛梁的關(guān)系,可以得到實(shí)梁上任意一點(diǎn)的撓度。為方便計(jì)算,將式(7)轉(zhuǎn)換為如下矩陣:
采用有限元數(shù)值模擬方法驗(yàn)證共軛梁法用于測斜儀中計(jì)算位移的準(zhǔn)確性。在Abaqus中建立長度為2 m,半徑分別為2 cm、3 cm和4 cm的懸臂梁模型。梁的彈性模量和泊松比分別為3 000 MPa和0.35。將梁的一端進(jìn)行完全固定約束,在懸臂梁上施加不同的載荷和位移,分別模擬:1)不同半徑懸臂梁在相同荷載作用下的工況,見圖4(a);2)同一半徑懸臂梁在不同荷載作用下的工況,如圖4(b)所示;3)對懸臂梁施加雙向位移的工況,見圖4(c)。
圖4 不同工況下的數(shù)值模擬
通過有限元分析,獲得懸臂梁表面沿長度方向的應(yīng)變值,采用共軛梁法將應(yīng)變值計(jì)算為位移值。計(jì)算出的位移值與Abaqus中直接輸出得到的位移值進(jìn)行比較。對比結(jié)果如圖5、圖6、圖7所示。
圖5 相同荷載作用下不同半徑懸臂梁的計(jì)算位移與模擬位移比較
圖6 不同荷載作用下同一半徑懸臂梁的計(jì)算位移與模擬位移比較
圖7 雙向彎曲下計(jì)算位移與模擬位移的比較
結(jié)果表明,共軛梁法計(jì)算的位移與有限元法計(jì)算的位移一致。在荷載不變、半徑變化的工況下,三組最大誤差分別為 0.01 mm/m、0.07 mm/m、0.64 mm/m。隨著半徑的減小,誤差會逐漸增大,因此在該測斜儀設(shè)計(jì)制作時(shí)應(yīng)對半徑進(jìn)行控制,但該工況的對比結(jié)果仍可說明該算法的可靠性。在半徑不變、荷載變化的工況下,三組最大誤差分別為0.15 mm/m、0.22 mm/m、0.29 mm/m,具有較高的精度。在施加雙向彎曲位移時(shí)計(jì)算結(jié)果與模擬結(jié)果吻合較好。不同工況下模擬和計(jì)算的對比結(jié)果說明共軛梁法將應(yīng)變轉(zhuǎn)化為位移的準(zhǔn)確性,并驗(yàn)證了可將該方法用于OFDR原位測斜儀,根據(jù)測得的應(yīng)變計(jì)算測斜儀的變形。
基于共軛梁法的OFDR原位測斜儀在進(jìn)行邊坡模型實(shí)驗(yàn)前,需要提前標(biāo)定。標(biāo)定裝置如圖8所示。該標(biāo)定裝置主要由鋼架、OFDR原位測斜儀、位移螺母、夾具和位移計(jì)組成。鋼架起整體支撐作用。鋼架底部橫梁上有一個(gè)固定夾,可以固定OFDR原位測斜儀(PVC棒)的一端,模擬固定端的邊界條件。中間橫梁有導(dǎo)向槽,用于控制OFDR原位測斜儀(PVC棒)在二維平面內(nèi)的變形。右端的位移螺母可以上下移動。實(shí)驗(yàn)過程中,轉(zhuǎn)動螺母使位移桿左右移動,從而對測斜儀的不同部位施加力和位移。用磁性底座將誤差為0.01 mm的位移計(jì)固定在規(guī)定的高度,測出測斜儀的實(shí)際水平位移。不考慮溫度與應(yīng)變的交叉敏感性,在測試過程中,OFDR原位測斜儀長度1 m,每隔1 cm取一個(gè)點(diǎn)進(jìn)行測量,傳感桿上1 m長的光纖共有100個(gè)應(yīng)變測點(diǎn)。通過共軛梁算法將得到的應(yīng)變值轉(zhuǎn)換為位移,并與百分表各點(diǎn)測量的實(shí)際位移值進(jìn)行比較。
圖8 標(biāo)定實(shí)驗(yàn)
除此之外還采用了FBG原位測斜儀與OFDR原位測斜儀進(jìn)行了對比,F(xiàn)BG測斜儀長度為1 m,每隔25 cm一個(gè)測點(diǎn),共有4個(gè)測點(diǎn)。將FBG測斜儀測得的位移值與OFDR測斜儀測得的位移進(jìn)行對比,驗(yàn)證OFDR測斜儀的測量精度。
為方便測試結(jié)果的直觀對比,將測量結(jié)果繪制為測點(diǎn)位置與水平位移的關(guān)系,如圖9所示。從圖9可以看出,基于共軛梁算法的OFDR原位測斜儀的測量值與百分表提供的參考值非常接近,并且相比于FBG原位測斜儀OFDR原位測斜儀測量結(jié)果更加精確。對測量的結(jié)果進(jìn)行統(tǒng)計(jì),采用相對誤差(RE)、平均相對誤差(MRE)和均方根誤差(RMSE)對標(biāo)定結(jié)果進(jìn)行了分析,結(jié)果見表1。OFDR原位測斜儀測量結(jié)果的標(biāo)準(zhǔn)偏差為0.441 mm,相對誤差為-9.66% ~ -1.18%,在允許的測量要求范圍內(nèi)。綜上,在實(shí)驗(yàn)室條件下驗(yàn)證了所提測量技術(shù)的可行性與可靠性。這種技術(shù)在實(shí)際條件下的可用性與在實(shí)驗(yàn)室條件下的差別不大,但是光纖的連接、斷裂、光彎曲損耗以及光纖的粘接或埋入工藝和傳感棒材料等因素都會影響光纖的信噪比和測量不確定度,這將是下一步的研究重點(diǎn)。
表1 實(shí)驗(yàn)測量結(jié)果統(tǒng)計(jì)
圖9 位移計(jì)測得的變形與使用共軛梁算法得到的計(jì)算位移的比較
為驗(yàn)證OFDR原位測斜儀的可靠性和實(shí)用性,進(jìn)行了1g條件下的邊坡模型試驗(yàn)。此外,還使用北京中弘泰科生產(chǎn)的誤差0.1 mm的RT-1010測斜儀測量邊坡的內(nèi)部水平位移。最后,對兩種測斜儀的監(jiān)測結(jié)果進(jìn)行了比較。
試驗(yàn)所用土樣為標(biāo)準(zhǔn)砂與高嶺土的混合物,按9∶1的混合比配制而成。本實(shí)驗(yàn)中使用的試驗(yàn)箱由鋼制成,尺寸為長1 440 mm,寬640 mm,高800 mm,如圖10所示。在填充模型邊坡的過程中,采用分層壓實(shí)方法。每次鋪設(shè)厚度控制在5 cm,并進(jìn)行壓實(shí),確保試驗(yàn)土具有良好的均質(zhì)性。OFDR原位測斜儀和RT-1010測斜儀垂直安裝在邊坡的同一截面指定位置,下端固定,如圖10所示。在試驗(yàn)過程中,為了研究OFDR原位測斜儀在邊坡變形監(jiān)測中的有效性,將液壓千斤頂放置在寬160 mm、長600 mm的坡頂鋼板上,利用液壓千斤頂向坡頂逐級施加載荷。
圖10 邊坡模型試驗(yàn)
在邊坡模型試驗(yàn)中,采用兩種測斜儀對邊坡內(nèi)水平位移進(jìn)行監(jiān)測,監(jiān)測結(jié)果如圖11所示,兩種測斜儀的測量結(jié)果具有良好的一致性。表明了OFDR原位測斜儀在邊坡監(jiān)測中的可靠性和實(shí)用性。隨著高度的增加,位移逐漸增加。這主要是因?yàn)榛瑒用嬉韵碌耐馏w相對穩(wěn)定,變形較小,而滑動面以上的土體變形更明顯。當(dāng)產(chǎn)生較大變形時(shí),OFDR原位測斜儀的測量結(jié)果與RT-1010測斜儀的測量結(jié)果會有一些偏差,主要是因?yàn)镺FDR原位測斜儀的位移算法基于小變形假定,隨著變形增大測量會有一定誤差,但是測量精度仍能夠滿足工程測量要求,具有一定的工程實(shí)用性。如何解決大變形測量將是下一步的研究重點(diǎn)。
圖11 兩類測斜儀的測量結(jié)果
本文借助分布式光纖測量的技術(shù)優(yōu)勢,開發(fā)了一種基于OFDR技術(shù)可用于邊坡深部水平位移監(jiān)測的測斜儀?;诠曹椓核惴ㄍ茖?dǎo)了光纖應(yīng)變與測斜儀位移的關(guān)系。通過有限元模擬,驗(yàn)證了共軛梁算法用于測斜儀位移計(jì)算中的有效性,并通過標(biāo)定校準(zhǔn)實(shí)驗(yàn)驗(yàn)證了OFDR原位測斜儀的有效性和精度。具體結(jié)論如下:
1)基于梁應(yīng)變和撓度之間的關(guān)系,可使用分布式光纖技術(shù)測量土體內(nèi)部的變形。室內(nèi)實(shí)驗(yàn)表明,利用分布式光纖計(jì)算梁撓曲變形的方法簡單,并可實(shí)現(xiàn)1 mm一個(gè)監(jiān)測點(diǎn)的分布式測量。測試得到的位移與百分表測量的實(shí)際位移非常接近,該方法測量精度高,輸出穩(wěn)定,實(shí)驗(yàn)相對誤差在±10%以內(nèi)。
2)相比于FBG原位測斜儀,基于分布式光纖的OFDR原位測斜儀每毫米取一個(gè)測點(diǎn),測點(diǎn)間隔更小,應(yīng)變測量精度為1 με,位移測量精度更高。同時(shí)較小的測點(diǎn)間隔更能準(zhǔn)確捕捉邊坡滑裂面的位置。
3)將OFDR原位測斜儀應(yīng)用于邊坡模型實(shí)驗(yàn)中,監(jiān)測結(jié)果與傳統(tǒng)測斜儀監(jiān)測結(jié)果一致,說明了OFDR原位測斜儀于實(shí)際監(jiān)測中具有一定的可靠性和工程實(shí)用性。