許有俊 車彥文 田治旺 崔廣琴

?

基于改進(jìn)反應(yīng)位移法的盾構(gòu)隧道橫向抗震規(guī)律研究1

許有俊1）車彥文1）田治旺2）崔廣琴1）

1）內(nèi)蒙古科技大學(xué)土木工程學(xué)院，內(nèi)蒙古自治區(qū)，包頭 014010 2）中國中鐵二院工程集團(tuán)有限責(zé)任公司，成都 610000

在地震動作用下，小間距盾構(gòu)隧道夾土層區(qū)域的動力反應(yīng)影響結(jié)構(gòu)的內(nèi)力與變形，而反應(yīng)位移法未考慮夾土層彈簧參數(shù)的變化，本文針對這一誤差來源，提出了改進(jìn)反應(yīng)位移法。一方面在反應(yīng)位移法合理內(nèi)容的基礎(chǔ)上完善計算模型，另一方面通過數(shù)值計算將改進(jìn)反應(yīng)位移法與動力時程法進(jìn)行對比分析，并擬合理論模型對改進(jìn)的方法進(jìn)行修正。結(jié)果表明，兩種算法下結(jié)構(gòu)的地震響應(yīng)規(guī)律一致；間距≥8m時，結(jié)構(gòu)地震動反應(yīng)變化不明顯，可按反應(yīng)位移法計算；理論模型提高了改進(jìn)反應(yīng)位移法計算結(jié)果的精度。

盾構(gòu)隧道 夾土層彈簧 改進(jìn)反應(yīng)位移法 地震響應(yīng) 理論模型

引言

隨著社會經(jīng)濟(jì)的快速發(fā)展，城市化程度的不斷提高，城市人口的急劇增加，城市地上空間的開發(fā)與利用已接近極限，由此促進(jìn)了城市地下空間的開發(fā)與利用，促進(jìn)了以地鐵盾構(gòu)隧道為代表的地下結(jié)構(gòu)工程的大規(guī)模興建，下層盾構(gòu)隧道（簡稱下層隧道）近距離穿越上層盾構(gòu)隧道（簡稱上層隧道）的工程也越來越多。如北京地鐵10號線三元橋站至亮馬橋站盾構(gòu)隧道左右線并行，最小凈距僅為1.7m，6號線南鑼鼓巷站至東四站區(qū)間上下行隧道疊落，間距僅為2.1m。其抗震問題成為城市工程抗震和防災(zāi)減災(zāi)研究的重要組成部分（權(quán)登州等，2015）。

針對地下結(jié)構(gòu)抗震問題的設(shè)計和計算方法，從力學(xué)特性上可分為動力時程法和擬靜力計算法兩大類。動力時程法能夠模擬地震波輸入過程中各個時刻結(jié)構(gòu)的內(nèi)力和變形狀態(tài)，準(zhǔn)確地反映地下結(jié)構(gòu)的地震響應(yīng)特性。當(dāng)然，在進(jìn)行動力分析時需要注意動力加載方式和人工邊界、力學(xué)阻尼及波的傳播，以便完整地模擬土體與結(jié)構(gòu)在地震過程的慣性作用以及相互作用。擬靜力法主要有地震系數(shù)法（施仲衡，1997）、自由場變形法（Hashash等，2001）、反應(yīng)位移法（川島一彥，1994）、Pushover分析方法（劉晶波等，2008）以及整體式反應(yīng)位移法（劉晶波等，2010）等。其中反應(yīng)位移法的基本原理是通過地層位移來計算地下結(jié)構(gòu)的地震響應(yīng)，在一定程度上反映地下結(jié)構(gòu)地震反應(yīng)的特點。該方法理論清晰、計算簡便，在眾多的設(shè)計規(guī)范和工程中得到了廣泛應(yīng)用。

目前反應(yīng)位移法只適用于單個結(jié)構(gòu)，但在強(qiáng)震作用下，小間距盾構(gòu)隧道夾土層的地震動反應(yīng)對結(jié)構(gòu)的內(nèi)力和變形影響很大，現(xiàn)采用的計算模型存在一定缺陷，會增大計算結(jié)果的誤差率。本文通過改進(jìn)反應(yīng)位移法的計算模型，對與下層隧道呈不同角度、間距、埋深等多種工況下的上層隧道橫斷面進(jìn)行抗震研究，并與準(zhǔn)確的動力有限元的模擬數(shù)值進(jìn)行對比分析，同時對兩種算法生成的數(shù)據(jù)進(jìn)行曲線擬合，得到的擬合理論模型對提高反應(yīng)位移法計算結(jié)果的精度具有重要的實際意義，也為后續(xù)深入研究奠定了基礎(chǔ)。

1 反應(yīng)位移法概述

在地震動作用下，盾構(gòu)隧道變形主要受周圍地層的約束。反應(yīng)位移法的基本原理是將地震荷載作用下地層周圍的剪力、結(jié)構(gòu)自身慣性力等施加于結(jié)構(gòu)，同時把地層在地震時產(chǎn)生的位移差通過地基彈簧以靜載的形式作用于結(jié)構(gòu)上，從而求得結(jié)構(gòu)的內(nèi)力和變形等。

反應(yīng)位移法最重要的一步是在結(jié)構(gòu)周圍設(shè)置地基彈簧，以模擬隧道與地層之間的相互作用，其系數(shù)可按計算公式或通過有限元方法求解。但隧道間距較小時，夾土層在強(qiáng)震作用下對結(jié)構(gòu)產(chǎn)生不可忽略的影響，若采用上述方法只分析單個隧道結(jié)構(gòu)，未考慮夾土層區(qū)域彈簧參數(shù)變化的影響，將使結(jié)構(gòu)內(nèi)力的計算結(jié)果產(chǎn)生誤差（唐金良，2014）。

2 改進(jìn)反應(yīng)位移法

2.1 改進(jìn)模型

為克服反應(yīng)位移法的計算模型未考慮夾土層彈簧變化這一不足，在反應(yīng)位移法內(nèi)容合理的基礎(chǔ)上，改進(jìn)反應(yīng)位移法的計算模型。

基于有限元結(jié)構(gòu)分析軟件SAP84，改進(jìn)計算模型：兩條隧道夾土層等效為連接結(jié)構(gòu)上兩點的直彈簧，如圖1所示。結(jié)構(gòu)的受力分析與其周圍其他地基彈簧按反應(yīng)位移法設(shè)置（蔣英禮等，2013）。

圖1 改進(jìn)反應(yīng)法計算模型

2.2 夾土層彈簧剛度

結(jié)合地基彈簧常數(shù)的近似取值法（林皋，1990），擬定夾土層彈簧的計算公式如下：

式中，為夾土層彈簧剛度；為與地震震動最大應(yīng)變幅度相應(yīng)的地基土的剪切模量；為夾土層彈簧反映的土體面積；為兩條隧道結(jié)構(gòu)的間距。

2.3 擬合理論模型

采用上述改進(jìn)模型對兩個結(jié)構(gòu)進(jìn)行計算，不可避免會產(chǎn)生一些偏差，為簡單有效地解決這個問題，通過與準(zhǔn)確的動力時程法進(jìn)行對比分析，提出了引入擬合理論模型來提高計算精度，使改進(jìn)方法與動力時程法的結(jié)果相吻合。

在Matlab軟件中對兩種算法生成的數(shù)據(jù)進(jìn)行擬合時，對函數(shù)對話框“Results”中的參數(shù)SSE（擬合誤差方差的平方和，其值越小，擬合結(jié)果越好）與-square（置信區(qū)間，越趨近1，擬合結(jié)果越好）進(jìn)行不斷對比和分析，同時參考Adjusted-square（調(diào)整置信區(qū)間，越趨近1，擬合度越高），選出較合理的理論模型（唐家德等，2008；史立新等，2007）。

3 算例

地鐵隧道的相對位置主要存在并行、斜穿與下穿這3類情況，本章利用改進(jìn)模型，在多種埋深、間距工況下，研究隧道的橫向地震響應(yīng)特性，并與動力時程法的數(shù)值結(jié)果進(jìn)行對比分析。

基于Midas-NX有限元分析軟件對算例進(jìn)行分析。為提高計算精度（蘭景巖等，2012），動力時程法所使用的地震波依據(jù)《北京地鐵16號線工程場地地震安全性評價報告》，采用提供的50年超越概率為10%的加速度時程曲線，見圖2。土體遵從經(jīng)典的Mohr-Coulomb破壞準(zhǔn)則，盾構(gòu)隧道結(jié)構(gòu)采用彈性模型。進(jìn)行靜力計算時，土層兩側(cè)采用水平約束的人工天然邊界；進(jìn)行動力計算時，在模型底部豎直方向施加地震動的加速度邊界，四周為粘彈性動力人工邊界；動力分析是在靜力計算獲得應(yīng)力場后進(jìn)行的。

圖2 輸入地震動加速度時程和加速度傅氏譜

3.1 計算工況

以北京地鐵17號線03標(biāo)段區(qū)間盾構(gòu)隧道為背景，襯砌外徑為6m，內(nèi)徑為5.4m，厚度為0.3m。通過查閱并統(tǒng)計北京市地鐵盾構(gòu)隧道資料，設(shè)定上下層隧道所呈角度＝0°、30°、60°、90°，間距＝2m、4m、6m、8m、10m、12m、14m，上層隧道埋深＝8m、10m、12m、14m、16m。隧道尺寸與計算簡圖詳見圖3。盾構(gòu)隧道周圍土層簡化為代表北京地區(qū)常見的均一粉質(zhì)黏土（施春花等，2009），結(jié)構(gòu)與土體的參數(shù)如表1所示。

圖3 隧道尺寸及計算簡圖

表1 土體、隧道計算參數(shù)

3.2 數(shù)值結(jié)果及分析

盾構(gòu)隧道設(shè)計主要依據(jù)結(jié)構(gòu)的彎矩值、軸力值及直徑變化率，其值越大，結(jié)構(gòu)變形越明顯。在進(jìn)行變化規(guī)律分析時，選取最大值進(jìn)行數(shù)據(jù)處理。

3.2.1 彎矩分析

（1）對比分析

呈不同角度的兩個隧道的變形模態(tài)基本一致，以上層隧道為例進(jìn)行論述。上層隧道彎矩的最大值隨埋深和間距呈規(guī)律性變化，見圖4，其中以動力時程法的模擬結(jié)果作參照，驗證改進(jìn)反應(yīng)位移法的準(zhǔn)確性。由圖4易知兩種算法下彎矩最大值的變化規(guī)律基本一致：

①通過圖4（a）與（b）可知，埋深對彎矩值的影響最大，間距次之；與下層隧道位于同一角度、間距工況下的上層隧道，其彎矩最大值與單孔隧道的最大值均隨埋深的增加而增加（楊寧波，2013）。

②通過圖4（c）與（d）可知，在0°并行的情況下，兩隧道的上層隧道在同一埋深。隨著間距增大，彎矩最大值由先小于單孔隧道，然后增大到與單孔隧道近乎相同。在地震動作用下，當(dāng)兩條隧道間距較小時，由于下層隧道-夾土層-上層隧道慣性作用增強(qiáng)，提高整體結(jié)構(gòu)穩(wěn)定性，從而減小了上層隧道的差異位移，對上層隧道起到抗震的效果。隨間距逐漸增大，抗震效果也隨之減弱。隨埋深的增加，上層隧道彎矩最大值趨于單孔隧道的速率越來越快。

③通過圖4（f）—（j）可知，在30°、60°并行與90°疊落的情況下，當(dāng)上層隧道在同一埋深下，隨著間距增大，彎矩最大值由先大于單孔隧道，然后減小到與單孔隧道近乎相同。下層隧道的開挖造成臨近上層隧道的一定范圍內(nèi)土層應(yīng)力的釋放，在強(qiáng)震中，松弛的地層明顯放大了上層隧道的地震反應(yīng)特性，使其彎矩值出現(xiàn)大于單孔隧道的可能。彎矩最大值在隧道間距處于4—6m時出現(xiàn)峰值，該范圍是設(shè)計雙線隧道的薄弱區(qū)域。

在橫斷面抗震分析中應(yīng)考慮內(nèi)力值的分布區(qū)域（以30°并行的上層隧道為例，如圖5所示），彎矩最大值主要分布在上部拱腰兩側(cè)。

采用改進(jìn)反應(yīng)位移法計算時，在同一埋深與角度下，當(dāng)隧道間距≥8m時，彎矩值的變化不大，并且趨近于單孔隧道，這說明當(dāng)兩條隧道間距≥8m時，可忽略夾土層地基彈簧的變化對上層隧道彎矩值的影響，按現(xiàn)有的反應(yīng)位移法進(jìn)行計算。

圖4 上層隧道彎矩最大值變化圖（二）

圖5 彎矩分布區(qū)域圖（30°并行，埋深12m，間距6m，單位：kN·m）

（2）擬合理論模型選取

基于Matlab軟件進(jìn)行計算，經(jīng)殘差曲線的對比分析，彎矩值擬合理論模型采用函數(shù)Cubic Polynomial（公式2）進(jìn)行擬合，效果良好，如圖6所示。

式中，f（x）為動力時程法彎矩值；x為改進(jìn)反應(yīng)位移法彎矩值；p1、p2、p3、p4為修正系數(shù)，其均值與取值范圍分別為－1.00×10－5（－4.88×10－5，2.88×10－5）、1.16×10－3（－1.42×10－2，1.65×10－2）、1.20（－0.76，3.15）、－8.58（－88.69，71.53）。

擬合優(yōu)度方面，擬合誤差平方和（SSE）為0.04；置信區(qū)間（-square）為0.95；調(diào)整置信區(qū)間（Adjusted-square）為0.95。

3.2.2 軸力分析

（1）對比分析

軸力與彎矩的分析方法一致，即對每種工況下的最大值進(jìn)行研究。兩種算法下軸力最大值的變化規(guī)律基本一致，同時與改進(jìn)位移法彎矩最大值的變化規(guī)律也一致，以30°并行的上層隧道為例，見圖7。

在地震動作用下，軸力最小值分布在隧道頂部區(qū)域，最大值基本上分布在結(jié)構(gòu)下部拱腰附近，集中現(xiàn)象比較明顯，以30°并行的上層隧道為例，見圖8。

（2）擬合理論模型選取

基于Matlab軟件進(jìn)行計算分析，軸力值擬合理論模型采用Sum of Sin Functions函數(shù)（公式3）進(jìn)行擬合，效果良好，曲線見圖9。

式中，（）為動力時程法軸力值；為改進(jìn)反應(yīng)位移法軸力值；1、1、1、2、2、2為修正系數(shù)，其均值與取值范圍分別為1.51×103（5.69×102，2.45×103）、6.75×10－4（－5.366×10－5，1.40×10－3）、6.25（5.98，6.51）、28.46（14.42，42.5）、1.41×10－2（9.02×10－3，1.92×10－2）、－4.90（－10.42，0.63）。

在擬合優(yōu)度方面，擬合誤差平方和（SSE）為0.07；置信區(qū)間（-square）為0.96；調(diào)整置信區(qū)間（Adjusted-square）為0.96。

圖8 軸力分布區(qū)域圖（30°并行，埋深12m，間距6m，單位：kN）

圖9 軸力最大值擬合曲線

3.2.3 直徑變化率分析

（1）對比分析

易知直徑變化率受間距影響不大，隨埋深增加而減小，且變化幅度很小；單孔隧道直徑變化率最大，而在與下層隧道呈90°疊落、60°、30°、0°并行的情況下，上層隧道的直徑變化率依次減小。該變化規(guī)律與動力時程法模擬的上層隧道的地震動反應(yīng)特性類似，見圖10。

（2）擬合理論模型選取

基于Matlab軟件與上述理論分析，用擬合理論模型的Fourier函數(shù)來擬合直徑變化率，效果良好。曲線見圖11。

式中，（）為動力時程法直徑變化率；為改進(jìn)反應(yīng)位移法直徑變化率；0、1、1、2、2、為修正系數(shù)，其均值與取值范圍分別為6.91×10－2（－5.83×10－2，1.97×10－1）、5.75×10－2（－2.46×10－2，1.40×10－1）、－2.30×10－2（－0.61，0.56）、1.76×10－2（－0.26，0.29）、1.51× 10－2（－0.21，0.24）、32.39（－14.13，78.91）。

在擬合優(yōu)度方面，擬合誤差平方和（SSE）為2.78×10－3；置信區(qū)間（-square）為0.97；調(diào)整置信區(qū)間（Adjusted-square）為0.96。

圖10 直徑變化率（動力時程法）

圖11 直徑變化率擬合曲線

表2 直徑變化率（‰）

續(xù)表

4 結(jié)語

（1）本文改進(jìn)了反應(yīng)位移法計算模型，通過在雙線隧道之間設(shè)置夾土彈簧來反映土體-結(jié)構(gòu)間的相互作用，提出了該彈簧剛度的計算公式。在此基礎(chǔ)上，將上層隧道橫斷面的計算結(jié)果與準(zhǔn)確的動力時程法模擬數(shù)據(jù)進(jìn)行對比分析。

（2）在與下層隧道呈不同角度、間距、埋深的多種工況下，上層隧道橫斷面的抗震研究結(jié)果表明，兩種方法所得的彎矩與軸力的最大值的變化規(guī)律基本一致：埋深對內(nèi)力值的影響最大，間距次之；隧道水平并行時，上層隧道在同一埋深下，隨著間距增大，內(nèi)力值由先小于單孔隧道，然后增大到與單孔隧道近乎相同；隧道由斜穿變化到豎直下穿時，上層隧道在同一埋深下，隨著間距增大，內(nèi)力值由先大于單孔隧道，然后減小到與單孔隧道近乎相同，并且在4—6m時出現(xiàn)峰值，該范圍是設(shè)計雙線隧道的薄弱區(qū)域。同時，上層隧道直徑最大變化率受間距影響不大，隨埋深增加而減小，且隧道所呈角度越大變化率越大。

（3）彎矩最大值主要分布在上部拱腰兩側(cè)；軸力最小值分布在隧道頂部區(qū)域，最大值基本上分布在結(jié)構(gòu)下部拱腰附近，集中現(xiàn)象比較明顯。

（4）當(dāng)隧道間距≥8m時，內(nèi)力最大值變化不大，并趨近于單孔隧道，這說明隧道間距≥8m時，可忽略夾土層地基彈簧的變化對上層隧道內(nèi)力值的影響，按現(xiàn)有的反應(yīng)位移法進(jìn)行計算。

（5）基于Matlab軟件與準(zhǔn)確的動力時程法，對改進(jìn)模型的彎矩、軸力最大值和直徑最大變化率分別采用Cubic Polynomial、Sum of Sin Functions、Fourier擬合理論模型進(jìn)行修正來減少誤差，擬合效果良好。

綜上所述，本文提出的改進(jìn)反應(yīng)位移法較為合理，適用于小間距盾構(gòu)隧道的抗震分析與設(shè)計。

蔣英禮，王勁松，2013．基于反應(yīng)位移法的盾構(gòu)隧道橫斷面抗震計算．現(xiàn)代隧道技術(shù)，50（6）：52—57．

蘭景巖，呂悅軍，劉紅帥，2012．地震動強(qiáng)度及頻譜特征對場地地震反應(yīng)分析結(jié)果的影響．震災(zāi)防御技術(shù)，7（1）：37—45．

林皋，1990．地下結(jié)構(gòu)抗震分析綜述（下）．世界地震工程，（3）：1—10，42．

劉晶波，劉祥慶，李彬，2008．地下結(jié)構(gòu)抗震分析與設(shè)計的Pushover分析方法．土木工程學(xué)報，41（4）：73—80．

劉晶波，王文暉，趙冬冬，2010．地下結(jié)構(gòu)地震反應(yīng)計算反應(yīng)位移法的改進(jìn)．土木建筑與環(huán)境工程，32（增刊2）：211—213．

權(quán)登州，王毅紅，井彥林等，2015．黃土地區(qū)地鐵地下結(jié)構(gòu)抗震研究綜述．震災(zāi)防御技術(shù)，10（3）：565—574．

施春花，呂悅軍，彭艷菊等，2009．北京地區(qū)粉質(zhì)粘土土動力學(xué)參數(shù)的統(tǒng)計分析．震災(zāi)防御技術(shù)，4（1）：69—79．

史立新，聶信天，季明，2007．基于Matlab曲線擬合工具箱的列表曲線擬合．新技術(shù)新工藝，（7）：39—41．

施仲衡，1997．地下鐵道設(shè)計與施工．西安：陜西科學(xué)技術(shù)出版社．

唐家德，2008．基于Matlab的非線性曲線擬合．計算機(jī)與現(xiàn)代化，（6）：15—19．

唐金良，2014．地鐵交叉隧道地震響應(yīng)特性研究．成都：西南交通大學(xué)．

楊寧波，2013．盾構(gòu)隧道地震反應(yīng)和影響因素分析．哈爾濱：中國地震局工程力學(xué)研究所．

姚毅超，董正方，王君杰，2012．盾構(gòu)隧道的反應(yīng)位移法與抗震性能指標(biāo)．中國科技論文，（5）：334—338．

川島一彥，1994．地下構(gòu)筑物的耐震設(shè)計．日本：鹿島出版社．

Hashash Y. M. A., Hook J. J., Schmidt B., et al., 2001. Seismic design and analysis of underground structures. Tunnelling and Underground Space Technology, 16(4): 247—293.

許有俊，車彥文，田治旺，崔廣琴，2017．基于改進(jìn)反應(yīng)位移法的盾構(gòu)隧道橫向抗震規(guī)律研究．震災(zāi)防御技術(shù)，12（3）：613—624．

Study of Transverse Anti-seismic Law in Shield Tunnel Based on the Improved Response Deformation Method

Xu Youjun1), Che Yanwen1), Tian Zhiwang2)and Cui Guangqin1)

1) The School of Civil Engineering, Inner Mongolia University of Science and Technology, Baotou 014010, Inner Mongolia, China 2) China Railway Eryuan Engineering Group Co. ltd, Chengdu 610000, China

The seismic response of soil layer, in shield tunneling with small interval, influences structural internal force and deformation. However, the change of soil layer spring was not considered in the response deformation method. The improved method is proposed to reduce the error source of the original method. In our study, the model is completed with improved model, which is based on rational content of the original method. Moreover, the improved method is compared with the dynamic analysis method by numerical calculation, and theory model is fit out to modify the method further. The analysis results show that seismic law under two methods is basically consistent; when the interval between two tunnels is equal or greater than 8m, the variation of seismic response of shield tunnel is little, and this can be done by the original method. This theory model increases the accuracy of the results from improved method.

Shield tunnel; Soil layer spring; The improved response deformation method; Seismic response; Theory model

10.11899/zzfy20170317

內(nèi)蒙古自治區(qū)高等學(xué)?？茖W(xué)技術(shù)研究項目（NJZY14167）

2016-11-17

許有俊，男，生于1979年。博士，教授，碩士生導(dǎo)師。主要從事隧道、地下工程等方面的教學(xué)和科研工作。E-mail：xyoujun@163.com