張益瑄,林文凱,晏啟祥,朱洪雪,劉 陽

(西南交通大學(xué)交通隧道工程教育部重點(diǎn)實(shí)驗(yàn)室,成都 610031)

近些年來，隨著我國在交通隧道領(lǐng)域和地下工程建設(shè)方面的不斷發(fā)展，TBM法具有環(huán)保、掘進(jìn)速度快等優(yōu)點(diǎn)，在許多地下工程中被采用，而在TBM法施工過程中，圍巖蠕變對(duì)隧道開挖的影響越來越大，已經(jīng)成為不可忽視的因素，近些年來，許多的學(xué)者研究了圍巖蠕變對(duì)隧道開挖的影響[1-9]，蔡燕燕等[10]考慮圍巖蠕變的全過程，推導(dǎo)出了深埋隧洞非線性位移解，侯公羽等[11]在支護(hù)反力變化條件下提出被動(dòng)支護(hù)反力連續(xù)變化下的迭代算法，并且結(jié)合工程實(shí)例對(duì)圍巖的蠕變變形進(jìn)行了分析。林文凱等[12]基于Burgers的蠕變模型對(duì)圓形隧道的內(nèi)力進(jìn)行了研究，對(duì)比分析了地層結(jié)構(gòu)模型與荷載結(jié)構(gòu)模型兩種分析方法，并得出了地層結(jié)構(gòu)模型分析更為有效的結(jié)論。王明年等[13]基于有限元的分析方法，結(jié)合二郎山隧道的工程實(shí)例，分析了高地應(yīng)力軟巖隧道的蠕變特性。同時(shí)，也有不少學(xué)者在蠕變位移方面進(jìn)行了研究，余東明等[14]為了描述深埋圓形隧道的開挖支護(hù)，基于Burgers與Drucker-parger組合的黏彈塑性模型，推導(dǎo)出了深埋圓形隧道考慮剪脹性能的圍巖蠕變位移解析解，結(jié)合實(shí)際算例，分析了圍巖剪脹力與支護(hù)反力對(duì)深埋圓形圍巖蠕變位移的影響規(guī)律。

目前在國內(nèi)，基于Burgers本構(gòu)模型，利用有限元分析軟件ANSYS等建立荷載結(jié)構(gòu)模型和地層結(jié)構(gòu)模型，分析圍巖蠕變作用下隧道襯砌結(jié)構(gòu)的內(nèi)力的研究已經(jīng)很多，本文做進(jìn)一步研究分析，建立地形復(fù)雜，地表崎嶇不平的三維地質(zhì)模型，在初始應(yīng)力作用下，實(shí)施隧道開挖數(shù)值模擬，以隧道開挖模擬的結(jié)果為基礎(chǔ)，分析在圍巖蠕變作用下，支護(hù)結(jié)構(gòu)拱頂、拱底及左右拱腰處的位移變化，從而得到圍巖的蠕變特性。

1 圍巖蠕變理論

綜合國內(nèi)外的巖體蠕變理論研究，目前常用組合模型的方式來模擬巖體的流變關(guān)系。組合模型的基本原理是根據(jù)巖體的彈性、塑性及黏滯性質(zhì)來設(shè)定一些基本元件，通過這些基本元件的串并聯(lián)組合，以求能夠反映巖體的本構(gòu)關(guān)系，并以此模擬巖體的應(yīng)力-應(yīng)變關(guān)系。常用基本元件基礎(chǔ)上組合的模型有:麥克斯韋體(Maxweii)模型和開爾文體(Keivin)組合模型，以及由M體和K體串聯(lián)所組合的伯格斯(Burgers)(M-K)體組合模型[6]。結(jié)合高黎貢山隧道所處的砂巖層地段開展研究，根據(jù)地層巖性特點(diǎn)，選用黏彈性Burgers蠕變模型開展隧道圍巖蠕變分析。Burgers模型的本構(gòu)方程為

σ+p1σ1+p2σ2=q1ε1+q2ε2

(1)

式中，E1、η1分別為Maxweii模型的彈性模量、黏滯系數(shù)；E2、η2分別為Keivin模型的彈性模量、黏滯系數(shù)；σ1、σ2分別為正應(yīng)力對(duì)時(shí)間的一階、二階導(dǎo)數(shù)；ε1、ε2分別為正應(yīng)變對(duì)時(shí)間的一階、二階導(dǎo)數(shù)。由于Burgers模型的4個(gè)參數(shù)不能作為Prony級(jí)數(shù)輸入值，因此需將E1、E2、η1、η2轉(zhuǎn)化為剪切模量的Prony級(jí)數(shù)參數(shù)[15]。對(duì)于Burgers模型，Prony剪切模量表達(dá)式為[12]

(2)

式中，g∞為松弛時(shí)間無限大時(shí)的剪切模量；g1，τ1分別為Maxweii模型的相對(duì)剪切模量和松弛時(shí)間；g2，τ2分別為Keivin模型的相對(duì)剪切模量和剪切時(shí)間；G0為松弛時(shí)間為0的剪切模量。

將Burgers模型中的彈性模量系數(shù)進(jìn)行轉(zhuǎn)換，并將獲得的表達(dá)式進(jìn)行拉普拉斯變換和拉普拉斯逆變換，便可以得到剪切模量的Prony級(jí)數(shù)表達(dá)式，通過對(duì)比式(2)，便可得到Burgers模型在ANSYS中需要的4個(gè)Prony級(jí)數(shù)g1，g2，τ1，τ2[16-17]，可用于下文的蠕變分析。

2 三維地質(zhì)模型與應(yīng)力場(chǎng)回歸方法

建立高黎貢山隧道砂巖層地段三維地質(zhì)模型。為了模擬地形復(fù)雜，地表崎嶇不平的模型，結(jié)合地質(zhì)平面圖中的等高線資料，導(dǎo)入CAD中得到dxf文件，然后用GoCAD軟件調(diào)用dxf文件，提取等高線上的相關(guān)擬合點(diǎn)坐標(biāo)數(shù)據(jù)。由于數(shù)據(jù)的不均勻性，為了便于這些提取數(shù)據(jù)在ANSYS中建模，借用Surfer軟件使數(shù)據(jù)均勻化，使數(shù)據(jù)達(dá)到均勻值，然后將處理后的數(shù)據(jù)導(dǎo)入ANSYS中，就可以得到地形復(fù)雜，崎嶇不平的三維地形地質(zhì)模型。結(jié)合高黎貢山實(shí)際條件，建立三維地質(zhì)山體模型，模型下部是一個(gè)長(zhǎng)方體模型，長(zhǎng)、寬、高尺寸為250 m×150 m×120 m，在山體下部模型中包含有隧道襯砌結(jié)構(gòu)模型，其中模型下部位置進(jìn)行隧道開挖模擬，襯砌的內(nèi)、外半徑分別為4、4.5 m。從而獲得ANSYS分析的整體模型如圖1所示。

圖1 數(shù)值模型

將圍巖視為各向同性連續(xù)介質(zhì)，隧道圍巖、襯砌力學(xué)參數(shù)見表1。

表1 周圍土體與襯砌參數(shù)

圍巖蠕變參數(shù)分別為：Maxweii模型的彈性模量、黏滯系數(shù)分別為E1=5.546 GPa，η1=32.697 GPa·h，E2=3.211 GPa，η2=2.378 GPa·h，由前面公式可得到4個(gè)Prony級(jí)數(shù)，分別為g1=0.118 2，g2=0.883 6，τ1=17.265 4，τ2=0.267 6。

地應(yīng)力場(chǎng)是隧道從開挖到支護(hù)等仿真分析的計(jì)算基礎(chǔ)，對(duì)合理進(jìn)行地下工程設(shè)計(jì)與施工具有很大的現(xiàn)實(shí)意義[16-17]，初始地應(yīng)力場(chǎng)也是影響隧道圍巖穩(wěn)定性的重要影響因素[18]，以現(xiàn)場(chǎng)部分實(shí)測(cè)地應(yīng)力數(shù)據(jù)為基礎(chǔ)，同時(shí)考慮地形、地質(zhì)方面的因素，結(jié)合理論研究與數(shù)值分析進(jìn)行反演[19-20]。本文考慮用多元線性回歸分析法對(duì)高黎貢山隧道砂巖層地段進(jìn)行初始地應(yīng)力回歸。

(3)

為了方便統(tǒng)計(jì)與計(jì)算，取m個(gè)相關(guān)點(diǎn)，由最小二乘法殘差平方和的方法得到

(4)

由最小二乘法的基本原理，可以得到回歸系數(shù)，為使得S殘最小，令S殘對(duì)Li的偏導(dǎo)數(shù)為0，即

(5)

整理后可以得到如下矩陣

(6)

(7)

j=1，2，3，…，6表示6個(gè)初始應(yīng)力分量。

山體原地應(yīng)力測(cè)試孔CS1、CS2分別位于模型x=40 m，y=75 m和x=140 m，y=75 m兩處附近，兩應(yīng)力測(cè)試孔的實(shí)測(cè)數(shù)據(jù)如表2所示，其中數(shù)據(jù)已轉(zhuǎn)化為模型坐標(biāo)方向的應(yīng)力分量。

進(jìn)行初始地應(yīng)力場(chǎng)的反演，分別對(duì)該復(fù)雜地形三維模型施加構(gòu)造應(yīng)力場(chǎng)，應(yīng)力場(chǎng)1、3分別為YZ面擠壓應(yīng)力場(chǎng)與豎向的剪切應(yīng)力場(chǎng)，應(yīng)力場(chǎng)2、3、5分別為XZ面的擠壓應(yīng)力場(chǎng)，水平與豎向剪切應(yīng)力場(chǎng)。

用相關(guān)軟件分別計(jì)算上述各構(gòu)造應(yīng)力場(chǎng)和自重應(yīng)力場(chǎng)，可以得到相關(guān)點(diǎn)深度的應(yīng)力分量值，以表2實(shí)際測(cè)量得到的應(yīng)力值為基礎(chǔ)，代入式(3)～式(7)進(jìn)行分析，計(jì)算可以得到該處圍巖在上述假設(shè)各構(gòu)造應(yīng)力場(chǎng)下的初始狀態(tài)，具體如式(8)回歸方程

表2 測(cè)試孔應(yīng)力分量

σ=1.073σ自重+6.374σ構(gòu)1+1.628σ構(gòu)2+
42.336σ構(gòu)3+0.433σ構(gòu)4+0.245σ構(gòu)5

(8)

σ構(gòu)1、σ構(gòu)2、σ構(gòu)3、σ構(gòu)4、σ構(gòu)5、σ自重分別對(duì)應(yīng)5個(gè)構(gòu)造應(yīng)力場(chǎng)和自重應(yīng)力場(chǎng)。對(duì)該處山體的回歸效果進(jìn)行分析，可以計(jì)算得到該處山體的初始應(yīng)力場(chǎng)回歸方程的相關(guān)系數(shù)r=0.964，接近于1.0。殘差平方和S殘=12.687，說明初始地應(yīng)力場(chǎng)回歸具有較好的相似度。

3 結(jié)果分析

3.1 復(fù)雜地形下構(gòu)造應(yīng)力場(chǎng)

運(yùn)用多元線性回歸法進(jìn)行分析，得到山體在構(gòu)造應(yīng)力場(chǎng)以及自重應(yīng)力場(chǎng)下的回歸方程，以式(8)和相關(guān)數(shù)據(jù)為基礎(chǔ)，可以得到相應(yīng)的邊界條件，運(yùn)用有限元軟件分析，復(fù)雜地形下山體地層的豎向位移和豎向應(yīng)力如圖2、圖3所示。

圖2 初始地層豎向位移云圖(單位：mm)

圖3 初始地層豎向應(yīng)力云圖(單位：Pa)

由圖2和圖3可知，受到山體地形復(fù)雜、崎嶇不平的影響，山頭地表處有較大的初始豎向位移，且在底層位置山體有較大的初始豎向應(yīng)力。

3.2 隧道開挖支護(hù)狀態(tài)

為了得到復(fù)雜地形下隧道開挖后的支護(hù)狀態(tài)，用有限元軟件進(jìn)行隧道開挖數(shù)值模擬，在上文假定的構(gòu)造應(yīng)力場(chǎng)與自重應(yīng)力場(chǎng)作用下，隧道開挖完支護(hù)結(jié)構(gòu)的豎向位移云圖見圖4。

圖4 隧道開挖支護(hù)襯砌豎向位移云圖(單位：mm)

從圖4可以看出，受到山頂復(fù)雜地形的影響，隧道縱向的中部范圍內(nèi)豎向位移較大，在隧道頂部按照梭形分布，在靠近隧道兩端的位置，豎向位移逐漸趨于平緩，這種分布與圖2初始地層豎向位移較為類似。

隧道開挖支護(hù)結(jié)構(gòu)豎向應(yīng)力云圖如圖5所示。

圖5 隧道開挖支護(hù)襯砌豎向應(yīng)力云圖(單位：Pa)

由圖5可以看出，隧道頂部沿縱向有較大的應(yīng)力分布，且應(yīng)力分布比較均勻，隧道頂部的應(yīng)力分布與山頭處的應(yīng)力分布較為相似。

3.3 圍巖蠕變下支護(hù)結(jié)構(gòu)行為分析

下文研究地形復(fù)雜三維地質(zhì)模型下，隧道開挖后，圍巖的蠕變對(duì)支護(hù)結(jié)構(gòu)的影響。為了更好地得到結(jié)果進(jìn)行分析，在建模過程中，取隧道模型的縱向長(zhǎng)度為300 m，在模型的縱向方向，分別選取位置不同，埋深不同的5個(gè)數(shù)值采集面，并按照順序進(jìn)行編號(hào)，數(shù)值檢測(cè)斷面的位置與埋深見表3。

表3 各數(shù)值采集斷面分布

其中，斷面1和斷面5位于模型縱向的兩側(cè)，位置上位于隧道埋設(shè)較淺處，斷面2、斷面3、斷面4位于模型縱向中部，相對(duì)于斷面1、斷面5，埋深較大。各數(shù)值采集斷面分布如圖6所示。

支護(hù)結(jié)構(gòu)的拱頂部位也設(shè)置相應(yīng)隧道縱向的數(shù)值采集斷面，各支護(hù)結(jié)構(gòu)數(shù)值采集斷面示意見圖6。支護(hù)結(jié)構(gòu)橫斷面監(jiān)測(cè)點(diǎn)布置見圖7，如監(jiān)測(cè)斷面2的拱頂為A2，拱底為B2，左側(cè)拱腰為C2，右側(cè)拱腰為D2。

圖6 數(shù)值監(jiān)測(cè)點(diǎn)斷面示意

圖7 監(jiān)測(cè)點(diǎn)橫向示意

以5個(gè)隧道橫向數(shù)值采集斷面和隧道縱向數(shù)值采集斷面的數(shù)據(jù)為基礎(chǔ)，分析研究隧道在圍巖蠕變荷載作用下，支護(hù)結(jié)構(gòu)的響應(yīng)特點(diǎn)。

3.3.1隧道各監(jiān)測(cè)斷面Ai處的位移

以隧道5個(gè)橫向數(shù)值采集斷面支護(hù)結(jié)構(gòu)Ai處位移隨圍巖蠕變變化的數(shù)據(jù)為基礎(chǔ)，得到Ai處位移隨時(shí)間變化的趨勢(shì)如圖8所示。

圖8 隧道各數(shù)值采集斷面Ai位移曲線

由圖8可知，在圍巖蠕變荷載作用下，位于隧道中部的斷面2、斷面3、斷面4橫向數(shù)值采集斷面的Ai處位移不斷增加，而位于隧道兩端的1、5數(shù)值采集斷面隨著時(shí)間的增加位移不斷減小。

圖8中，圍巖蠕變產(chǎn)生的初始階段，A3處沉降為38.71 mm，是產(chǎn)生最大沉降位置，A1處沉降為28.87 mm，是所有拱頂位置處的最小沉降值，這與在初始應(yīng)力場(chǎng)作用下支護(hù)結(jié)構(gòu)的位移分布是相似的，因?yàn)锳3相比于A1埋深較大。由于A1相對(duì)于A5有更大的埋深，A1處的位移減小率更大，在斷面2、斷面3、斷面4中，在圍巖蠕變作用下，可以發(fā)現(xiàn)A3處位移增加的速率是最快的，原因是斷面3相對(duì)于其他斷面埋深較大。

3.3.2隧道各數(shù)值采集斷面Bi處的位移

通過對(duì)5個(gè)橫向數(shù)值采集斷面得到的數(shù)據(jù)進(jìn)行匯總分析，得到B1處位移隨時(shí)間變化的趨勢(shì)如圖9所示。

圖9 隧道各數(shù)值采集斷面Bi處位移曲線

由圖9可知，受到蠕變荷載的影響，1與5數(shù)值采集斷面支護(hù)結(jié)構(gòu)Bi的位移隨著天數(shù)增加不斷減小， 隧道斷面2、斷面3、斷面4橫向數(shù)值采集斷面的B1處位移是不斷增加的。Bi處的位移變化趨勢(shì)與Ai處的位移變化趨勢(shì)是類似的。

在圖9中，圍巖蠕變產(chǎn)生的初始階段，B3處位移的大小為22.01 mm，是產(chǎn)生位移的最大位置，B1沉降為15.98 mm，是所有拱頂位置處沉降的最小值。同時(shí)可以發(fā)現(xiàn)，在斷面2、斷面3、斷面4中，斷面3擁有較大的埋深，因此A3處位移最大，在2、3、4斷面中，在圍巖蠕變隨時(shí)間作用下，埋深越大的位置處的位移增長(zhǎng)速率也越大。

3.3.3隧道各監(jiān)測(cè)斷面Ci與Di的收斂位移

通過對(duì)5個(gè)橫向數(shù)值采集斷面得到的數(shù)據(jù)進(jìn)行匯總分析，得到Ci與Di的收斂位移隨時(shí)間變化的趨勢(shì)如圖10所示。

圖10 隧道各監(jiān)測(cè)斷面Ci與Di的收斂位移曲線

由圖10可以看出，隧道橫向數(shù)值采集斷面2、斷面3、斷面4隨著時(shí)間增加，左右拱腰收斂位移呈減小趨勢(shì)，其中埋深較大的斷面3處收斂位移減小最大，斷面1和斷面5處左右拱腰收斂位移呈現(xiàn)增大趨勢(shì)，由于斷面1和斷面5埋深接近，增長(zhǎng)趨勢(shì)也較為接近。在圍巖蠕變作用下，通過對(duì)比分析Ai與Bi處的位移，可以發(fā)現(xiàn)Ci與Di之間的收斂位移的變化趨勢(shì)與拱頂和拱底處的位移變化趨勢(shì)恰好是相反的。

在圖10中，圍巖蠕變產(chǎn)生的初始階段，Ci與Di位置的收斂位移4.62 mm是收斂位移最大值處，相對(duì)地，C3與D3位置的收斂位移大小為3.89 mm，是收斂位移的最小值處。由圖中各曲線的變化趨勢(shì)，在圍巖蠕變的作用下，隨著斷面埋深的不斷增大，相應(yīng)左右拱腰的收斂位移不斷減小，與上述各斷面Ai、Bi處位移隨蠕變變化相比，Ci與Di位置的收斂位移數(shù)值變化幅度較小，并且變化的數(shù)值也比較小。

3.3.4 隧道縱向各數(shù)值采集斷面拱頂位移分析

為了具體分析研究縱向拱頂位移隨蠕變的變化趨勢(shì)，這里沿著隧道的縱向方向，在拱頂位置取一系列采集點(diǎn)，采集點(diǎn)的范圍為10～240 m，每隔10 m共取24個(gè)點(diǎn)，研究蠕變作用天數(shù)下位移的變化趨勢(shì)，結(jié)果如圖11所示。

圖11 隧道橫向數(shù)值采集斷面拱底位移曲線

結(jié)合圖11進(jìn)行分析，在地形崎嶇復(fù)雜的山體中間部位，沿著隧道縱向進(jìn)行分析，相對(duì)埋深較大處，支護(hù)結(jié)構(gòu)的拱頂位移大于山體兩邊處的拱頂位移。隨著圍巖蠕變作用不斷增加，將圖11根據(jù)變化趨勢(shì)分為3個(gè)區(qū)間，其中區(qū)間I，在隧道縱向x=10 m到x=60 m范圍內(nèi)，支護(hù)結(jié)構(gòu)拱頂位移隨著蠕變時(shí)間的增加不斷減??；區(qū)間Ⅱ，沿隧道縱向區(qū)間x=60 m到x=200 m，支護(hù)結(jié)構(gòu)拱頂位移在蠕變荷載作用下，隨著時(shí)間是不斷增加的；區(qū)間Ⅲ，沿隧道縱向區(qū)間x=200 m到x=240 m范圍內(nèi)，隨著荷載作用時(shí)間的增加，拱頂位移不斷減小，區(qū)間Ⅱ范圍內(nèi)形成了一個(gè)拱形曲線，觀察圖2發(fā)現(xiàn)，支護(hù)結(jié)構(gòu)拱頂位移隨著時(shí)間逐漸形成一個(gè)山頭，與地形崎嶇復(fù)雜的三維地質(zhì)模型的地表形狀是類似的。

其中，區(qū)間Ⅰ處的拱頂位移隨著蠕變時(shí)間增加，變化幅度大于區(qū)間Ⅲ的位移變化，在圖11中表現(xiàn)為區(qū)間Ⅰ處的曲線比區(qū)間Ⅲ的曲線更為稀疏，在三維地質(zhì)模型中表現(xiàn)為區(qū)間Ⅰ處的地表地形更加復(fù)雜陡峭。

3.3.5同一數(shù)值采集斷面下支護(hù)結(jié)構(gòu)位移隨時(shí)間變化分析

為了研究圍巖蠕變對(duì)支護(hù)結(jié)構(gòu)的響應(yīng)，選取上文的5個(gè)監(jiān)測(cè)斷面，分別提取它們?cè)贏i，Bi處的位移以及Ci，Di兩處的收斂位移進(jìn)行分析。

數(shù)值采集斷面1支護(hù)結(jié)構(gòu)的位移隨時(shí)間變化如圖12所示。

圖12 數(shù)值采集斷面1支護(hù)結(jié)構(gòu)位移曲線

由圖12可以看出，數(shù)值采集斷面1在整個(gè)蠕變分析過程中，拱頂處的位移要大于拱底處的位移，拱頂處和拱底處的位移都隨蠕變時(shí)間的增加而減小，但減小的速率不斷增大，左右拱腰處的位移隨時(shí)間緩慢增加，增加速率較為穩(wěn)定。

數(shù)值采集斷面3支護(hù)結(jié)構(gòu)的位移隨時(shí)間變化曲線如圖13所示。

圖13 數(shù)值采集斷面3支護(hù)結(jié)構(gòu)位移曲線

由圖13可以看出，數(shù)值采集斷面3在蠕變時(shí)間增加的過程中，拱頂?shù)奈灰埔笥诠暗椎奈灰疲绊敽凸暗孜灰剖遣粩嘣黾拥?，增加的趨?shì)較為穩(wěn)定，而左右拱腰的位移隨蠕變時(shí)間增加不斷減小，整體上監(jiān)測(cè)斷面3左右拱腰的位移較小。

數(shù)值采集斷面5支護(hù)結(jié)構(gòu)的位移隨時(shí)間變化如圖14所示。

圖14 數(shù)值采集斷面5支護(hù)結(jié)構(gòu)位移曲線

由圖14可以看出，在蠕變時(shí)間增加的過程中，拱頂位移數(shù)值要大于拱底位移，但拱頂位移和拱底位移都隨時(shí)間不斷減小，左右拱腰位移隨蠕變時(shí)間呈平緩增加的趨勢(shì)。

4 結(jié)論

通過建立地形復(fù)雜，崎嶇不平的三維地質(zhì)模型，在初始應(yīng)力作用下，實(shí)施隧道開挖數(shù)值模擬，分析隧道開挖模擬后在圍巖蠕變作用下，支護(hù)結(jié)構(gòu)不同部位的位移變化，對(duì)支護(hù)結(jié)構(gòu)5個(gè)橫向數(shù)值采集斷面的Ai，Bi，Ci，Di處，以及1個(gè)縱向數(shù)值采集斷面進(jìn)行相關(guān)數(shù)據(jù)提取，從而得出以下結(jié)論。

(1)在山體中間部位，支護(hù)結(jié)構(gòu)的拱頂Ai處的位移在蠕變荷載作用下，位移隨時(shí)間變化的趨勢(shì)是不相同的，其中，埋深相對(duì)較淺的斷面A1和A5處拱頂位移在圍巖蠕變荷載作用下，位移呈逐漸減小的趨勢(shì)，而位于隧道中部的A2，A3，A4隨著蠕變時(shí)間的增加，位移不斷增大，原因是上述3個(gè)斷面埋深相對(duì)較大。

(2)支護(hù)結(jié)構(gòu)Bi處的拱底位移受到圍巖蠕變荷載的影響，隨著時(shí)間的增加，Bi處的位移變化趨勢(shì)類似于Ai處的拱頂位移變化，同時(shí)Bi處的位移數(shù)值要小于Ai處。

(3)Ci與Di之間的收斂位移的變化趨勢(shì)也隨著斷面位置的變化而不同，其中C2與D2，C3與D3，C4與D4之間的收斂位移在蠕變荷載作用下，隨時(shí)間呈現(xiàn)減小的趨勢(shì)，C1與D1，C5與D5之間的收斂位移則隨著時(shí)間的增大不斷增大。

(4)通過分析同一數(shù)值采集斷面下支護(hù)結(jié)構(gòu)位移隨時(shí)間變化曲線，發(fā)現(xiàn)拱頂處與拱底處的位移變化趨勢(shì)類似，左右拱腰處的收斂位移變化趨勢(shì)相反，位于隧道兩端的斷面1和5各位置的位移變化趨勢(shì)是相似的，而位于中部的數(shù)值采集斷面各位置處則與兩端位置有著相反的變化趨勢(shì)。