粘彈性邊界輸入及大直徑隧道地震反應(yīng)數(shù)值模擬

[摘要] """為了確定粘彈性邊界輸入地震動時底部波場的選擇，并在此基礎(chǔ)上研究兩條并行大直徑盾構(gòu)隧道的地震響應(yīng)，建立了相應(yīng)的數(shù)值模型進行分析。首先通過對粘彈性邊界地震動輸入波動法的研究，使用自由場和散射場算例對比了求解模型底部等效節(jié)點力時幾種不同的波場選擇方式。然后，使用粘彈性邊界和前面的地震動輸入方法，建立了大直徑不同間距的土-隧道有限元模型，輸入不同頻譜特性的地震動，對土-隧道體系地震響應(yīng)問題進行分析。研究結(jié)果表明：粘彈性邊界以等效節(jié)點力的方式輸入地震動時，采用自由波場求解等效節(jié)點力精度最高，采用基于人工邊界子結(jié)構(gòu)的地震動輸入方法也有較高精度；地震動頻譜特性對隧道的動力響應(yīng)影響較大，隧道間距對隧道的動力響應(yīng)影響相對較??；隧道間距對隧道襯砌的彎矩和軸力的影響存在差異。研究結(jié)果可為使用粘彈性邊界時地震動的準(zhǔn)確輸入和大直徑盾構(gòu)隧道抗震設(shè)計提供參考。

[關(guān)鍵詞] 粘彈性邊界; 大直徑隧道; 地震響應(yīng); 動力相互作用

[DOI] 10.19987/j.dzkxjz.2024-013

基金項目："中國地震局工程力學(xué)研究所基本科研業(yè)務(wù)費專項資助項目（2019B10）資助。

0 "引言

隨著工程技術(shù)的快速發(fā)展，大量的隧道建設(shè)開始向大直徑、大深度、長距離和高烈度場地發(fā)展^[1-2]。同時，隨著城市地下空間利用率日益提高，雙線隧道或隧道相互鄰近的情況不可避免。大量的研究和震害現(xiàn)場調(diào)查表明，穿越復(fù)雜場地的隧道襯砌往往存在開裂、坍塌、滲水等破壞情況^[3-4]，這些已引起了 工程設(shè)計人員和研究者的高度重視。

目前，有很多學(xué)者對地震作用下土-隧道之間的相互作用進行研究，主要研究方法有理論推導(dǎo)、原型觀測、模型試驗和數(shù)值模擬等。理論推導(dǎo)形式復(fù)雜，且往往局限于各向同性彈性場地和線彈性襯砌等理想情況^[5]，限制了其實際應(yīng)用。原型觀測^[3-4]結(jié)果可靠，其缺點是往往只能觀測結(jié)果而不能獲得破壞時的響應(yīng)。模型試驗主要有擬靜力試驗、振動臺試驗和離心機振動臺試驗。其中，鄒炎等^[6]通過振動臺試驗研究了穿過土層分界面的隧道地震反應(yīng)，認(rèn)為此類隧道產(chǎn)生破壞的主要原因是隧道相對位移的出現(xiàn)；Wang等^[7]通過振動臺試驗研究了土-隧道-結(jié)構(gòu)相互作用系統(tǒng)的地震反應(yīng)規(guī)律，發(fā)現(xiàn)隧道結(jié)構(gòu)的存在對周圍土體和地表建筑的影響不同，且與地震動頻譜特性有關(guān)；王彤輝等^[8]通過振動臺試驗分析了地震動特性對場地-地下結(jié)構(gòu)-地上結(jié)構(gòu)相互作用體系的影響；Baziar等^[9]通過離心機振動臺試驗及數(shù)值模擬分析，發(fā)現(xiàn)隧道的存在對地表加速度影響有正有負(fù)；李延濤等^[10]通過振動臺試驗研究了在上部結(jié)構(gòu)存在時土-上下平行隧道體系的地震響應(yīng)；Liang等^[11]通過振動臺試驗和數(shù)值模擬研究了超大直徑盾構(gòu)隧道上軟下硬場地的響應(yīng)。在數(shù)值模擬方面的研究也很多，有直接對土-隧道體系進行分析的，如王國波等^[12]，黃忠凱和張冬梅^[13]，張季等^[14]；還有通過反應(yīng)位移法等簡化方法的，如Kontoe等^[15]對圓形隧道抗震設(shè)計的4個解析解進行了驗證，并驗證了使用等效線性化來模擬非線性的近似方法；Zhao等^[16]通過數(shù)值模擬采用改進反應(yīng)位移法研究了大直徑海底盾構(gòu)隧道混凝土襯砌在地震荷載作用下的損傷模式。當(dāng)?shù)卣饎幼饔糜谕?隧道體系時，其動力相互作用是復(fù)雜的，受到場地和地震動頻譜特性、有無并行隧道及其間距等的影響，有必要對該體系的動力響應(yīng)進行進一步研究。

地下隧道的地震反應(yīng)計算，是半無限域中土-結(jié)動力相互作用問題，在進行涉及到地基半無限性影響的動力計算時，需要使用人工邊界模擬半無限域?qū)τ邢抻嬎銋^(qū)域的作用，使得散射波穿過人工邊界或被人工邊界吸收而不反射回到計算區(qū)域。正確的地震動輸入方法是數(shù)值模擬得到準(zhǔn)確結(jié)果的前提。本文首先通過SV波垂直入射下的自由場和散射場算例給出了使用粘彈性邊界時合理的底邊界等效節(jié)點力計算方法；然后，采用粘彈性邊界和前面提到的地震動輸入方法，建立了實際工程中典型的上軟下硬場地雙線大直徑盾構(gòu)隧道的數(shù)值計算模型，研究地震作用下隧道襯砌的受力情況并探討了地震動頻譜特性和隧道間距對隧道襯砌動力響應(yīng)的影響。本研究對復(fù)雜場地大直徑隧道的抗震設(shè)計和性能分析有一定的參考作用。

1 "粘彈性邊界與地震動輸入

1.1 "粘彈性邊界

ABAQUS自帶的無限元在應(yīng)用于動力計算時，本質(zhì)上是粘性邊界^[17]，因此也具有粘性邊界在低頻動力荷載作用下模型容易整體漂移的缺點，粘彈性邊界克服了粘性邊界這一缺點，能吸收散射波并模擬半無限地基的恢復(fù)能力^[18]。粘彈性邊界首先由Deeks和Randolph^[19]提出；隨后劉晶波等^[18]在此基礎(chǔ)上給出了二維時域粘彈性邊界和外源荷載輸入方法，并發(fā)展了三維時域粘彈性邊界^[20]；后來，劉晶波等^[21]又提出了一種基于子結(jié)構(gòu)的外源地震動輸入方法，避免了輸入地震動時對自由場應(yīng)力進行求解。粘彈性邊界概念清晰，實現(xiàn)簡單，穩(wěn)定性良好，因此被廣泛應(yīng)用^{[14，16，22]}。

粘彈性邊界是在有限的土體模型邊界上施加彈簧和阻尼器元件來模擬半無限土體的，其彈簧系數(shù)K和阻尼系數(shù)C分別如下：

將地震動合理輸入到計算區(qū)域是取得正確數(shù)值模擬結(jié)果的前提。使用粘彈性邊界，需要將地震動轉(zhuǎn)化為人工邊界節(jié)點上的等效節(jié)點力實現(xiàn)地震動輸入，如果使用強制的位移或加速度邊界條件輸入會導(dǎo)致粘彈性邊界失效^{[22， 24]}。根據(jù)劉晶波等的研究^[18]，施加在粘彈性邊界上一點P的力為：

1.2 "自由場和散射場工況及計算模型

根據(jù)馬笙杰等^[22]，譚輝等^[24]的研究，在數(shù)值模擬時應(yīng)當(dāng)在土體側(cè)面和底面截斷處都施加粘彈性邊界，側(cè)邊界輸入波場應(yīng)為自由波場?；诖?，在本節(jié)中，保持側(cè)邊界均為自由波場，僅考慮底部波場的不同，再與劉晶波等^[21]提出的基于人工邊界子結(jié)構(gòu)的地震動輸入方法對比。具體的計算方法見表1。

1.3 "不同方法下的結(jié)果分析

圖3為自由場各方法的計算結(jié)果，從自由場的計算結(jié)果中可以看到，求解底部邊界等效節(jié)點荷載時速度采用入射波場或自由波場，位移采用自由波場能獲得最佳的模擬精度，如果位移選擇了入射波場，則無論速度輸入自由波場還是入射波場都存在較大的誤差，采用子結(jié)構(gòu)法輸入地震動的模擬精度也較高。在粘彈性邊界中，反射波在邊界上引起的荷載一部分由阻尼器吸收，該部分荷載不需要額外輸入來平衡，剩余部分由彈簧來平衡。當(dāng)?shù)走吔缟系牡刃Ч?jié)點力采用方法1，3的公式進行計算時，位移項僅考慮了入射波場，沒有考慮反射波對底邊界位移的影響，因此計算結(jié)果與理論值存在較大差異。同時，發(fā)現(xiàn)方法1和方法3，方法2和方法4結(jié)果完全相同，這是因為在底部邊界，當(dāng)二維SV波垂直入射時，自由場只有水平位移。因此，根據(jù)胡克定律，自由場應(yīng)力表示為：

圖4為散射場A，B，C，D四個點在不同方法下水平方向位移與遠(yuǎn)置邊界計算結(jié)果的對比，從圖4可以看出，方法2，4，5與遠(yuǎn)置邊界擬合度較高，但都不是完全重合，這是因為粘彈性邊界是一種近似的人工邊界，在推導(dǎo)過程中采用了一系列假設(shè)，但是從計算結(jié)果來看，完全滿足工程精度要求；方法1，3如前文所述，其計算結(jié)果存在較大偏差。

綜上所述，當(dāng)采用波動法作為粘彈性邊界輸入方法時，底邊界建議采用方法2，4所列公式，可以高精度地模擬計算域有散射波的情況；基于人工邊界子結(jié)構(gòu)的地震動輸入方法也具有較高的精度?；谌斯み吔缱咏Y(jié)構(gòu)的地震動輸入方法使用方便，可用于不同的應(yīng)力型人工邊界，但需要對人工邊界子結(jié)構(gòu)模型進行一次動力計算，速度較慢。

2 "隧道-土體模型數(shù)值模擬

2.1 "模型參數(shù)

為了分析不同地震動作用下不同間距的大直徑隧道的響應(yīng)，在ABAQUS中建立了不同隧道間距的土-隧道模型，計算模型見圖5，其中，B表示并行隧道的間距。選擇90 m厚的上軟下硬土層進行分析，為了減小邊界對計算結(jié)果的影響，土體水平向?qū)挾热?00 m。隧道外直徑為15 m，襯砌厚度為0.65 m，隧道最大埋深為27.5 m，采用4節(jié)點平面應(yīng)變單元離散土體和襯砌。為了滿足波傳播的條件^[25]，土體單元的最大尺寸為1 m×1 m，并在隧道附近適當(dāng)加密，土體劃分了約45000個單元（不同工況方法單元數(shù)略有不同），每個隧道襯砌劃分了1920個單元。隧道襯砌和土體之間為有限滑移接觸，其中，法向為“硬”接觸，切向為罰剛度算法，摩擦系數(shù)為0.6。為了更貼近隧道的真實受力情況，參考王蘇等^[26]的研究，對模型進行了地應(yīng)力平衡。地震動以SV波形式從土體底部入射，為了考慮地震動頻譜特性對隧道動力響應(yīng)的影響，分別輸入El Centro波，Northridge波，Kobe波，并將其調(diào)幅到0.1g，輸入地震動的加速度時程曲線見圖6。

采用等效線性化模型考慮土體的非線性，利用一維土層反應(yīng)分析程序EERA計算自由場響應(yīng)，獲得土體的等效阻尼比和等效剪切模量。將土體的等效阻尼比以瑞利阻尼的形式施加，以自由場的一階和二階頻率計算瑞利阻尼參數(shù)。土層信息見表2，各類土的剪切模量比、阻尼比與剪應(yīng)變的關(guān)系曲線見圖7。隧道襯砌選擇為C60混凝土，采用基于規(guī)范的混凝土損傷塑性（CDP）模型來考慮混凝土的彈塑性，具體實現(xiàn)方式參考方自虎等^[27]的研究。混凝土密度為2400 kg/m³，彈性模量為36000 MPa，泊松比為0.18。根據(jù)土層參數(shù)和自由場響應(yīng)，采用自編Python腳本批量施加粘彈性邊界和等效節(jié)點荷載。

2.2 "數(shù)值模擬結(jié)果

分別考察在輸入不同地震動，隧道間距為20 m、 25 m、 30 m、40 m和50 m的雙線隧道以及單孔隧道情況下襯砌彎矩（圖8）和軸力（圖9）的不同，來分析地震動頻譜特性和隧道間距對隧道受力的影響。彎矩以使襯砌外側(cè)受拉為正，軸力以受壓為正，為方便起見隧道襯砌的環(huán)寬取1 m。

2.2.1 "代表性計算結(jié)果

由于土-隧道模型為對稱結(jié)構(gòu)，當(dāng)?shù)卣饎訌牡撞看怪比肷鋾r，左右隧道受力也是對稱的，因此，只分析一側(cè)的隧道。圖8給出了在0.1g的El Centro波入射時單一隧道以及不同間距并行隧道中左側(cè)隧道襯砌的總彎矩沿著圓周方向的分布圖。從圖8可以看出，襯砌所受的彎矩在水平和豎直方向上最大，水平位置附近為外側(cè)受拉，豎直方向附近為內(nèi)側(cè)受拉。并行隧道的彎矩與單一隧道有較大區(qū)別，且與并行隧道的間距關(guān)系較大，如當(dāng)B=20 m、50 m時，并行左側(cè)隧道襯砌的最大彎矩分別為548.5 kN·m、603.2 kN·m，比單一隧道襯砌的最大彎矩625.1 kN·m分別下降了約12.3%、3.5%。和單一隧道相比，在并行隧道拱頂拱底附近和2個隧道內(nèi)側(cè)相對的水平位置附近，單一隧道襯砌的彎矩大于并行隧道，而在并行隧道內(nèi)側(cè)上下45°方向則相反，并行隧道的間距越小，這一趨勢越明顯。在并行隧道靠外的一側(cè)，單一隧道襯砌的彎矩和并行隧道襯砌的彎矩相差不大。

圖9給出了在0.1g的El Centro波入射時單一隧道以及不同間距并行隧道中左側(cè)隧道襯砌的總軸力沿著圓周方向的分布圖。由圖9可知，隧道襯砌軸力以壓力為主，最小值出現(xiàn)在隧道頂部位置，在底部和中部也較小，最大值出現(xiàn)在2個隧道下半部相對方向的45°位置。并行隧道襯砌所受的軸力比單一隧道要大一些，兩個隧道的間距越小這種趨勢越明顯。

2.2.2 "場地的地震響應(yīng)

為了說明地震動頻譜特性對場地響應(yīng)的影響，提取了不同入射波作用下自由場中隧道所在土層相對于隧道底部土體的最大位移的絕對值（圖10）。由圖10可知，地震動頻譜特性對場地的最大相對位移影響較大。在0.1g的El Centro波、Northridge波和Kobe波作用下，隧道頂部土層相對于隧道底部土層的最大相對位移分別為6.71 mm、5.51 mm和3.84 mm，由于隧道的受力和變形主要由土體的相對位移引起，因此可以預(yù)見，在El Centro波作用下，隧道襯砌的動彎矩和動軸力最大，在Kobe波作用下，隧道襯砌的動彎矩和動軸力最小。

2.2.3 "隧道所受動力荷載與隧道間距和地震動頻譜特性的關(guān)系

用隧道襯砌受力的最大最小值減去地應(yīng)力平衡分析步的隧道襯砌受力，即可得到隧道襯砌的最大動彎矩和動軸力，以此來分析并行隧道的動力相互作用。

首先考慮地震動頻譜特性的影響，在El Centro波、Northridge波和Kobe波作用下，單孔隧道襯砌所受的最大正動彎矩分別為206.58 kN·m，162.40 kN·m，107.63 kN·m?？梢?，不同地震動入射時，隧道襯砌所受的動軸力和動彎矩有很大的區(qū)別。

圖11為不同間距，不同地震動入射時，隧道襯砌所受動彎矩的包絡(luò)圖。可以看出，隧道襯砌動彎矩的最大值出現(xiàn)在襯砌上下共軛45°方向上，且襯砌所受正負(fù)動彎矩的最大值位置不同。隨著隧道間距的減小，雙線隧道襯砌所受正負(fù)動彎矩在2個隧道相對側(cè)的上下45°方向增加明顯，有逐漸向水平方向移動的趨勢，而在另一側(cè)變化不大。以左側(cè)隧道的負(fù)動彎矩為例，當(dāng)B=20 m、25 m、30 m時，El Centro波入射下的并行左側(cè)隧道的負(fù)動彎矩相較單孔隧道提高了15.6%、9.5%、5.8%，在Northridge波入射時，提高值為8.6%、3.3%、3.4%，而在Kobe波入射時，分別提高了12.9%、5.4%、2.0%。這說明并行隧道間距對地震動作用下襯砌動彎矩的影響很大，這種影響主要體現(xiàn)在2個隧道相靠近的一側(cè)，且與地震動頻譜特性的關(guān)系也較大。

圖12為不同間距，不同地震動入射時，隧道所受動軸力的包絡(luò)圖。可以看出，隧道所受動軸力的最大值位置和動彎矩類似，出現(xiàn)在襯砌上下共軛45°方向上。隨著隧道間距的減小，襯砌所受動軸力逐漸增大，增大位置出現(xiàn)在并行隧道相對側(cè)的上部45°和另一側(cè)下部45°位置，而在其他位置變化不大，甚至有可能略有減小。以右側(cè)隧道的正動軸力為例，當(dāng)B=20 m、25 m、30 m時，El Centro波入射下并行右側(cè)隧道的正動軸力相較單孔隧道提高了9.5%、5.7%、3.2%，在Northridge波入射時，提高值為5.5%、3.6%、3.7%，而在Kobe波入射時，提高值為10.4%、4.8%、1.8%。這說明并行隧道間距對地震動作用下襯砌動軸力的影響很大，這種影響主要體現(xiàn)在2個隧道相靠近的一側(cè)的上半部和另一側(cè)的下半部，且受地震動頻譜特性的影響，地震動頻譜特性對襯砌動軸力和動彎矩的影響存在不同。

3 "結(jié)論

本文從地震動輸入波動法節(jié)點等效荷載的計算公式入手，在使用中發(fā)現(xiàn)理論公式的計算結(jié)果有一定偏差，通過自由場和散射場算例驗證了底邊界幾種輸入方法的計算精度，并通過分析計算結(jié)果與理論解和遠(yuǎn)置邊界計算結(jié)果的差異，得出精確程度最高的地震動輸入方式。然后，利用有限元軟件建立了土-大直徑隧道模型，采用粘彈性邊界和前面提出的地震動輸入方式，考慮地震動頻譜特性和隧道間距的影響，得出了以下結(jié)論：

（1）粘彈性邊界采用等效節(jié)點荷載方式輸入地震動時，側(cè)面邊界輸入自由波場，底部位移輸入自由波場，速度輸入入射波場或自由波場，能得到最高的模擬精度，基于人工邊界子結(jié)構(gòu)的地震動輸入方法精度也較高。

（2）隧道在地震動作用下的響應(yīng)受地震動頻譜成分的影響較大，在設(shè)計時應(yīng)當(dāng)綜合考慮多條地震動。

（3）地震作用下，隧道襯砌所受的動彎矩、動軸力峰值位置出現(xiàn)在襯砌的共軛45°方向，與靜力作用下襯砌彎矩和軸力最大值的位置不同。

（4）地震動作用下雙線隧道和單一隧道的受力存在差異，隨著雙線隧道的間距逐漸減小，2個隧道之間的動力相互作用會越來越強烈，且這種作用對動彎矩和動軸力的影響不同。

參考文獻

洪開榮，馮歡歡. 近2年我國隧道及地下工程發(fā)展與思考（2019—2020年）[J]. 隧道建設(shè)（中英文），2021，41（8）：1259-1280 """Hong K R，F(xiàn)eng H H. Development and thinking of tunnels and underground engineering in China in recent 2 years （from 2019 to 2020）[J]. Tunnel Construction，2021，41（8）：1259-1280

《中國公路學(xué)報》編輯部. 中國交通隧道工程學(xué)術(shù)研究綜述：2022[J]. 中國公路學(xué)報，2022，35（4）：1-40 """Editorial Department of China Journal of Highway and Transport. Review on China’s traffic tunnel engineering research：2022[J]. China Journal of Highway and Transport，2022，35（4）：1-40

李天斌. 汶川特大地震中山嶺隧道變形破壞特征及影響因素分析[J]. 工程地質(zhì)學(xué)報，2008，16（6）：742-750 """Li T B. Failure characteristics and influence factor analysis of mountain tunnels at epicenter zones of great Wenchuan earthquake[J]. Journal of Engineering Geology，2008，16（6）：742-750

崔光耀，王明年，于麗，等. 汶川地震公路隧道洞口結(jié)構(gòu)震害分析及震害機理研究[J]. 巖土工程學(xué)報，2013，35（6）：1084-1091 """Cui G Y，Wang M N，Yu L，et al. Seismic damage and mechanism of portal structure of highway tunnels in Wenchuan earthquake[J]. Chinese Journal of Geotechnical Engineering，2013，35（6）：1084-1091

Lee V W，Trifunac M D. Response of tunnels to incident SH-waves[J]. Journal of the Engineering Mechanics Division，1979，105（4）：643-659

鄒炎，景立平，李永強. 隧道穿過土層分界面振動臺模型試驗研究[J]. 巖石力學(xué)與工程學(xué)報，2014，33（增刊1）：3340-3348 """Zou Y，Jing L P，Li Y Q. Study of shaking table model test of tunnel through soil interface[J]. Chinese Journal of Rock Mechanics and Engineering，2014，33（S1）：3340-3348

Wang G B，Yuan M Z，Yu M，et al. Experimental study on seismic response of underground tunnel-soil-surface structure interaction system[J]. Tunnelling and Underground Space Technology，2018，76：145-159

王彤輝，李小軍，陳紅娟，等. 場地-隧道-地上建筑結(jié)構(gòu)體系地震響應(yīng)的振動臺試驗[J]. 防災(zāi)減災(zāi)工程學(xué)報，2023，43（4）：742-751 """Wang T H，Li X J，Chen H J，et al. Shaking table test of seismic response of site-tunnel-aboveground building structure system[J]. Journal of Disaster Prevention and Mitigation Engineering，2023，43（4）：742-751

Baziar M H，Moghadam M R，Kim D S，et al. Effect of underground tunnel on the ground surface acceleration[J]. Tunnelling and Underground Space Technology，2014，44：10-22

李延濤，田野，宗金輝，等. 受相鄰上部結(jié)構(gòu)影響的隧道-土體系振動臺試驗研究[J]. 振動與沖擊，2020，39（3）：233-241，259 """Li Y T，Tian Y，Zong J H，et al. Shaking-table tests for seismic response of tunnel-soil system affected by adjacent upper structure[J]. Journal of Vibration and Shock，2020，39（3）：233-241，259

Liang J W，Xu A Q，Ba Z N，et al. Seismic response of ultra-large diameter shield tunnel in upper-soft and lower-hard site：Shaking table tests and numerical simulations[J]. Soil Dynamics and Earthquake Engineering，2023，164：107634

王國波，王亞西，陳斌，等. 隧道-土體-地表結(jié)構(gòu)相互作用體系地震響應(yīng)影響因素分析[J]. 巖石力學(xué)與工程學(xué)報，2015，34（6）：1276-1287 """Wang G B，Wang Y X，Chen B，et al. Analysis of factors influencing seismic responses of tunnel-soil-ground structural system[J]. Chinese Journal of Rock Mechanics and Engineering，2015，34（6）：1276-1287

黃忠凱，張冬梅. 軟土地層淺埋盾構(gòu)隧道地震易損性分析[J]. 現(xiàn)代隧道技術(shù)，2020，57（5）：51-60 """Huang Z K，Zhang D M. Seismic vulnerability analysis of the shallow-buried shield tunnel in soft soil stratum[J]. Modern Tunnelling Technology，2020，57（5）：51-60

張季，譚燦星，許開成. 地震作用下軟土-隧道-地上框架體系非線性動力反應(yīng)分析[J]. 振動與沖擊，2021，40（12）：159-167 """Zhang J，Tan C X，Xu K C. Nonlinear dynamic response analyses of a soft soil-tunnel-aboveground frame structure system under earthquake[J]. Journal of Vibration and Shock，2021，40（12）：159-167

Kontoe S，Avgerinos V，Potts D M. Numerical validation of analytical solutions and their use for equivalent-linear seismic analysis of circular tunnels[J]. Soil Dynamics and Earthquake Engineering，2014，66：206-219

Zhao Z Y，Cui J，Liu C，et al. Seismic damage characteristics of large-diameter shield tunnel lining under extreme-intensity earthquake[J]. Soil Dynamics and Earthquake Engineering，2023，171：107958

Lysmer J，Kuhlemeyer R L. Finite dynamic model for infinite media[J]. Journal of the Engineering Mechanics Division，1969，95（4）：859-877

劉晶波，呂彥東. 結(jié)構(gòu)-地基動力相互作用問題分析的一種直接方法[J]. 土木工程學(xué)報，1998，31（3）：55-64 """Liu J B，Lü Y D. A direct method for analysis of dynamic soil-structure interaction[J]. China Civil Engineering Journal，1998，31（3）：55-64

Deeks A J，Randolph M F. Axisymmetric time-domain transmitting boundaries[J]. Journal of Engineering Mechanics，1994，120（1）：25-42

劉晶波，王振宇，杜修力，等. 波動問題中的三維時域粘彈性人工邊界[J]. 工程力學(xué)，2005，22（6）：46-51 """Liu J B，Wang Z Y，Du X L，et al. Three-dimensional visco-elastic artificial boundaries in time domain for wave motion problems[J]. Engineering Mechanics，2005，22（6）：46-51

劉晶波，譚輝，寶鑫，等. 土-結(jié)構(gòu)動力相互作用分析中基于人工邊界子結(jié)構(gòu)的地震波動輸入方法[J]. 力學(xué)學(xué)報，2018，50（1）：32-43 """Liu J B，Tan H，Bao X，et al. The seismic wave input method for soil-structure dynamic interaction analysis based on the substructure of artificial boundaries[J]. Chinese Journal of Theoretical and Applied Mechanics，2018，50（1）：32-43

馬笙杰，遲明杰，陳紅娟，等. 黏彈性人工邊界在ABAQUS中的實現(xiàn)及地震動輸入方法的比較研究[J]. 巖石力學(xué)與工程學(xué)報，2020，39（7）：1445-1457 """Ma S J，Chi M J，Chen H J，et al. Implementation of viscous-spring boundary in ABAQUS and comparative study on seismic motion input methods[J]. Chinese Journal of Rock Mechanics and Engineering，2020，39（7）：1445-1457

劉晶波，杜義欣，閆秋實. 粘彈性人工邊界及地震動輸入在通用有限元軟件中的實現(xiàn)[J]. 防災(zāi)減災(zāi)工程學(xué)報，2007，27（增刊1）：37-42 """Liu J B，Du Y X，Yan Q S. Implementation of viscous-elastic artificial boundary and earthquake dynamic input in common finite element software[J]. Journal of Disaster Prevention and Mitigation Engineering，2007，27（S1）：37-42

譚輝，劉晶波，王東洋，等. 地下結(jié)構(gòu)地震反應(yīng)分析中人工邊界條件和地震波動輸入方法對比研究[J]. 工程力學(xué)，2018，35（增刊1）：212-216，222 """Tan H，Liu J B，Wang D Y，et al. Comparison on artificial boundaries and seismic wave input methods in seismic analysis of underground structures[J]. Engineering Mechanics，2018，35（S1）：212-216，222

廖振鵬. 工程波動理論導(dǎo)論[M]. 北京：科學(xué)出版社，2002：7-9 """Liao Z P. Introduction to wave motion theories in engineering[M]. Beijing：Science Press，2002：7-9

王蘇，路德春，杜修力. 地下結(jié)構(gòu)地震破壞靜-動力耦合模擬研究[J]. 巖土力學(xué)，2012，33（11）：3483-3488 """Wang S，Lu D C，Du X L. Research on underground structure seismic damage using static-dynamic coupling simulation method[J]. Rock and Soil Mechanics，2012，33（11）：3483-3488

方自虎，周?？?，賴少穎，等. ABAQUS混凝土應(yīng)力-應(yīng)變關(guān)系選擇[J]. 建筑結(jié)構(gòu)，2013，43（增刊2）：559-561 """Fang Z H，Zhou H J，Lai S Y，et al. Choose of ABAQUS concrete stress-strain curve[J]. Building Structure，2013，43（S2）：559-561

Seismic motion input methods of viscous-elastic boundary and numerical simulation of seismic response on large diameter tunnels

Jing Weibiao， Wang Zhan， Jing Liping Qi Wenhao

Institute of Engineering Mechanics， China Earthquake Administration， Heilongjiang Harbin 150080， China

[Abstract] """"To determine the selection of input wave field at the bottom of the viscous-elastic boundary， and study the seismic response of two parallel large-diameter shield tunnels， corresponding numerical models were established for analysis. First， through the study of wave method for seismic motion input of viscous-elastic boundary， several different wave fields selection methods for calculating the equivalent nodal forces at the bottom of the model were compared using free field and scattering field models. Then， using viscous-elastic boundary and the previous seismic motion input method， established soil-tunnels system finite element models with different interval between two large-diameter tunnels and input seismic waves with different spectral characteristics analyzed the dynamic response of the soil-tunnels system. The results show that the simulation accuracy is highest when calculating the equivalent nodal loads by the free wave field， and the seismic input method based on the substructure of artificial boundaries also has high accuracy; spectral characteristics of seismic motion has a significant impact on the dynamic response of tunnels， while the influence of two tunnels interval is relatively small; the influence of tunnels interval on the bending moment and axial force of tunnels lining is different. The research results can provide reference for accurately seismic motion input of viscous-elastic boundary and seismic design of large-diameter shield tunnels.

[Keywords] viscous-elastic boundary; large diameter tunnel; seismic response; dynamic interaction