陳博
(中鐵第五勘察設(shè)計(jì)院集團(tuán)有限公司,北京102600)
盾構(gòu)法隧道施工對(duì)周圍環(huán)境影響小、施工速度快等,已成為城市隧道施工的主要手段。然而掌子面失穩(wěn)造成的地表變形甚至坍塌事故在北京、上海等地都曾發(fā)生過(guò)。因此,掌子面失穩(wěn)研究是盾構(gòu)施工技術(shù)中的關(guān)鍵一環(huán)。
模型試驗(yàn)研究由于可以模擬實(shí)際工況下的掌子面失穩(wěn),為了排除土壓力盒、應(yīng)變片等測(cè)量元件對(duì)模型試驗(yàn)的影響,采用對(duì)稱理論和相似原理并運(yùn)用PIV(粒子圖像測(cè)速技術(shù))技術(shù)獲得干砂地層不同埋深比(C/D=0.5、1.0、2.0、3.0、4.0、5.0。其中,C為覆土厚度;D為隧道直徑)下地表變形規(guī)律及掌子面破壞模式。
本試驗(yàn)圖像處理部分采用基于MATLAB 平臺(tái)開發(fā)的粒子圖像測(cè)速程序MPIV,即(Particle Image Velocimetry in Matlab)。針對(duì)時(shí)間先后順序的數(shù)碼圖像,MPIV 能夠運(yùn)用相關(guān)性算法計(jì)算2 張圖像間同一顆粒位移,進(jìn)一步計(jì)算最大速度、平均速度,繪制顆粒速度場(chǎng)、位移場(chǎng)及位移等值圖。
MPIV 處理圖像時(shí)提供了2 種算法:一種是基于快速傅立葉變換的互相關(guān)法(COR);另一種是最小二次微分法(MQD)。2 種算法均是將圖像看作是隨時(shí)間變化的離散二維信號(hào)變化場(chǎng),進(jìn)行信號(hào)分析,通過(guò)計(jì)算2 幅圖像的相關(guān)性函數(shù)值得到粒子位移。如圖1 所示,在圖像1 和圖像2 上劃分相同區(qū)塊,f1(Xi,Yxj)和f2(Xi+ΔX,Yj+ΔY)分別代表圖像1 和圖像2 上的對(duì)應(yīng)區(qū)塊,大小為L(zhǎng)×L像素,其中L一般為64 或32。利用相關(guān)性函數(shù)計(jì)算f1與f2的匹配峰值,以算法設(shè)定的峰值取值作為該區(qū)塊內(nèi)粒子平均位移。區(qū)別于流體力學(xué)中要求計(jì)算流體速度,本試驗(yàn)只需測(cè)出前后2 個(gè)狀態(tài)的顆粒位移即可,所得位移值即是最終結(jié)果[1]。
圖1 MPI V圖像分析過(guò)程圖
試驗(yàn)?zāi)P拖涑叽鐬?7cm×52cm×91cm,采用10mm 普通玻璃配合鋼之框架黏合。隧道采用4mm 厚、內(nèi)徑10cm、長(zhǎng)20cm的半圓鋼管,盾構(gòu)支撐面板采用10mm 厚鋼制活塞,利用直徑1.5cm、長(zhǎng)25cm 的螺桿控制活塞,活塞與隧道和玻璃間涂抹凡士林消除摩擦影響,填土表面采用百分表代替位移計(jì)測(cè)量最大沉降量,測(cè)試系統(tǒng)如圖2 所示。
圖2 測(cè)試系統(tǒng)示意圖
本試驗(yàn)用土取自北京周邊河砂(粒徑d50=0.38mm),經(jīng)過(guò)水洗、晾曬、篩分,取0.1~2.0mm 粒徑砂,含水率為0.16%,其主要物理性質(zhì)指標(biāo)為:比重GS=2.62,天然密度ρ=1.51g/cm3,最小干密度ρmin=1.46g/cm3,最大干密度ρmax=1.82g/cm3,相對(duì)密實(shí)度Dr=0.55。經(jīng)篩分試驗(yàn)得知該砂土為級(jí)配良好均勻中砂,符合試驗(yàn)要求,砂土進(jìn)行100~400kPa 直剪試驗(yàn)得到摩擦角φ=37.3°,試驗(yàn)時(shí)采用砂雨法制備模擬地層砂土,制得模擬地層密度為1.62~1.65g/cm3,相對(duì)密度約為0.55。
為直觀體現(xiàn)掌子面破壞模式,不埋設(shè)測(cè)量元件,拓展非接觸式測(cè)量室內(nèi)模型試驗(yàn)。本試驗(yàn)運(yùn)用對(duì)稱原理,探討覆土為中砂的干砂地層掌子面失穩(wěn)問(wèn)題。試驗(yàn)采用位移控制法,利用螺桿帶動(dòng)活塞后退模擬掌子面支護(hù)力不足,螺桿加工螺距d=1mm(即操作桿每轉(zhuǎn)一圈,活塞移動(dòng)1mm)。步驟為:(1)固定并調(diào)整數(shù)碼相機(jī),調(diào)整光源、背景,線控拍攝初始圖像,讀取位移計(jì)初始值;(2)采取活塞每后退0.5mm(螺桿轉(zhuǎn)動(dòng)半圈),拍攝照片并讀取沉降數(shù)據(jù);(3)當(dāng)掌子面失穩(wěn)坍塌傳遞至地表,認(rèn)為試驗(yàn)結(jié)束;(4)當(dāng)活塞位移-地表沉降量曲線出現(xiàn)變化較大拐點(diǎn),迅速增大,則視為掌子面破壞。
依上述方法進(jìn)行覆土厚度為0.5D、1.0D、2.0D、3.0D、4.0D、5.0D試驗(yàn),運(yùn)用MPIV 處理圖像,研究不同埋深比下掌子面破壞模式、土拱效應(yīng)及地表沉降發(fā)展規(guī)律。
通過(guò)地表位移計(jì)和數(shù)碼照片得到掌子面失穩(wěn)地表沉降變化規(guī)律、掌子面前方土拱區(qū)域和破壞模式。
如圖3 所示,地表沉降量的變化與P-S曲線(壓力-位移曲線)類似,分為4 個(gè)階段,即無(wú)位移段、線性段、發(fā)展破壞段和破壞段。
第一階段無(wú)位移段,由于土拱效應(yīng)的影響,地表未表現(xiàn)出沉降,土體發(fā)揮抗剪強(qiáng)度。
第二階段線性段,土體處于彈性變形階段,繼續(xù)發(fā)揮抗剪強(qiáng)度,由于掌子面支護(hù)力迅速減小,土體抗剪強(qiáng)度充分發(fā)揮,該階段較為短暫。
第三階段發(fā)展破壞段,沉降量增加速率加快,活塞位移-沉降量曲線出現(xiàn)明顯拐點(diǎn),土體處于塑性變形階段,并逐漸喪失強(qiáng)度。
第四階段破壞段,土體喪失強(qiáng)度,位移量迅速增加,與活塞位移略呈線性增加,掌子面破壞。
對(duì)比0.5D~5.0D對(duì)應(yīng)關(guān)系曲線,隨著覆土厚度的增加,曲線線性段斜率逐漸增大,持續(xù)時(shí)間更長(zhǎng);拐點(diǎn)更加明顯,即第三階段表現(xiàn)更為強(qiáng)烈;第四階段斜率大致相等。
圖3 0.5D~5.0D 沉降量- 活塞位移關(guān)系圖
根據(jù)覆土厚度為2D試驗(yàn)所得掌子面前方土體位移場(chǎng)圖,可看出土拱效應(yīng),活塞位移較小時(shí),掌子面前方土體位移區(qū)域較小,隨著活塞位移,土體位移區(qū)域逐漸擴(kuò)大,最后延伸至地表。
本試驗(yàn)根據(jù)對(duì)稱原理設(shè)計(jì),隧道與盾構(gòu)面板等均采用半圓,利用玻璃箱可直觀了解掌子面破壞模式。通過(guò)粒子圖像測(cè)速技術(shù),對(duì)掌子面失穩(wěn)過(guò)程所得圖像進(jìn)行分析,獲得不同埋深比(C/D=0.5、1.0、2.0、3.0、4.0、5.0)情況下掌子面前方土體破壞模式。
當(dāng)覆土厚度為0.5 倍或1.0 倍隧道直徑時(shí),由位移場(chǎng)可知掌子面前方砂層破壞呈開口狀,破壞模式近似為漏斗狀,地表下1/2 隧道直徑范圍內(nèi),破壞影響范圍較大。當(dāng)覆土厚度2.0 倍隧道直徑以上時(shí),掌子面前方土體破壞模式更接近煙囪狀。地表沉降槽范圍不隨埋深比的變化而變化,受土拱效應(yīng)影響,相同活塞位移下,沉降槽范圍與最大沉降量隨著埋深比增大而減小,最大沉降點(diǎn)位于掌子面前方約0.25D處。
通過(guò)對(duì)不同埋深比(C/D=0.5、1.0、2.0、3.0、4.0、5.0)干砂地層掌子面失穩(wěn)試驗(yàn)地表沉降規(guī)律、土拱效應(yīng)、掌子面破壞模式的分析,得出以下結(jié)論:
1)由于掌子面支護(hù)力不足導(dǎo)致的掌子面失穩(wěn),引起地表沉降。不同埋深比試驗(yàn)規(guī)律相似,活塞位移-沉降量曲線將沉降過(guò)程分為無(wú)位移段、線性段、發(fā)展破壞段和破壞段4 個(gè)階段。
2)土拱效應(yīng)是除了掌子面支護(hù)力、支護(hù)方法、地下水等外另一影響掌子面失穩(wěn)的因素,隨著掌子面位移,前方土體經(jīng)歷土拱形成階段、土拱最大階段、土拱消失階段和破壞階段,當(dāng)?shù)乇沓两颠M(jìn)入發(fā)展破壞段后,判定掌子面失穩(wěn),土體變形加速。
3)當(dāng)埋深比為0.5 和1.0 時(shí),掌子面破壞模式為漏斗狀;當(dāng)埋深比大于1.0 時(shí),掌子面破壞模式更接近煙囪狀。相同活塞位移,地表沉降槽范圍隨埋深比增加而減?。徽谱用嫱耆茐暮蟪两挡鄯秶苈裆畋扔绊?;沉降最大點(diǎn)約位于掌子面前方0.25D處。
4)基于MATLAB 開發(fā)的圖像處理軟件MPIV 能夠成功描繪掌子面前方土體變化。PIV 技術(shù)的運(yùn)用,未在土體中埋設(shè)測(cè)力、測(cè)位移元件,避免因測(cè)量元件的埋入影響土體本身性質(zhì)研究。