劉天添，　戚承志，　陳昊祥，　楊謹(jǐn)瑞

(北京建筑大學(xué) 土木與交通工程學(xué)院 北京市高校工程結(jié)構(gòu)與新材料工程研究中心， 北京　100044)

?

海底隧道地震動(dòng)力響應(yīng)研究

劉天添，戚承志，陳昊祥，楊謹(jǐn)瑞

(北京建筑大學(xué) 土木與交通工程學(xué)院 北京市高校工程結(jié)構(gòu)與新材料工程研究中心， 北京100044)

摘要：高水壓作用下海底隧道的地震動(dòng)力響應(yīng)是一個(gè)需要深入研究的問(wèn)題. 基于Winkler地基模型，利用動(dòng)力平衡條件和微分關(guān)系推導(dǎo)出了僅由彎曲引起的梁振動(dòng)方程，利用分離變量法求得了海底隧道的振型和自振頻率，并推導(dǎo)出了地震作用下海底隧道與地基之間的動(dòng)力相互作用系數(shù)以及隧道結(jié)構(gòu)位移表達(dá)式. 通過(guò)算例計(jì)算研究了水壓對(duì)于結(jié)構(gòu)的自振頻率、結(jié)構(gòu)位移以及相互作用系數(shù)的影響. 研究表明，高水壓作用下結(jié)構(gòu)的地震動(dòng)力響應(yīng)不可忽略.

關(guān)鍵詞：海底隧道； 水壓力； Winkler地基模型； 自振頻率； 振動(dòng)位移

目前的跨海交通主要有橋梁和隧道兩種形式.也有橋隧并舉的形式. 而海底隧道因?qū)Ｑ蟓h(huán)境及水上航運(yùn)影響小、受惡劣氣候及戰(zhàn)時(shí)破壞影響小等優(yōu)點(diǎn)，在近年來(lái)跨海交通的不同方案比選過(guò)程中越來(lái)越顯現(xiàn)其優(yōu)越性，并日益成為城市、地區(qū)、甚至國(guó)家之間的重要連接手段和交通命脈. 但是相較一般的陸上隧道，海底隧道又有自己的難點(diǎn). 如海底隧道深埋于地下，深水海底地質(zhì)勘察比地面地質(zhì)勘察難度更高、投入更大，不僅影響勘察的準(zhǔn)確性，而且漏勘和失真情況的風(fēng)險(xiǎn)也會(huì)相應(yīng)增加；海底隧道周邊圍巖裂隙中一般會(huì)有較高的孔隙水壓力和滲水壓力，會(huì)造成巖體弱化，使圍巖的有效應(yīng)力和成拱能力降低，進(jìn)而影響地層的穩(wěn)定性；隧道絕大部分位于海床之下，襯砌需要承受長(zhǎng)期的高水壓作用；在穿越斷層構(gòu)造及地震破碎帶時(shí)需要高度關(guān)注地震作用對(duì)隧道安全性的影響等. 這些難點(diǎn)的存在使得海底隧道從設(shè)計(jì)、施工到運(yùn)營(yíng)維護(hù)等一系列過(guò)程都比一般隧道更加復(fù)雜，同時(shí)也推進(jìn)著海底隧道的研究.

在抗震研究早期，人們認(rèn)為地下構(gòu)造物因受到圍巖的約束和限制，在地震時(shí)隨著地基的運(yùn)動(dòng)而運(yùn)動(dòng)，地震對(duì)地下構(gòu)造物的破壞影響相比地表破壞小得多，可以不用考慮. 根據(jù)統(tǒng)計(jì)，過(guò)去地震中很少有地下隧道結(jié)構(gòu)遭受地震破壞似乎也驗(yàn)證了這一點(diǎn). 然而在1995年日本阪神大地震中，神戶市區(qū)間隧道和地鐵車(chē)站等大型地下結(jié)構(gòu)均遭受了重大破壞. 這一事實(shí)迫使人們重新審視地下結(jié)構(gòu)物的破壞，并開(kāi)始懷疑傳統(tǒng)的地下結(jié)構(gòu)抗震理念，自此地下結(jié)構(gòu)抗震研究的熱潮逐漸興起. 目前國(guó)內(nèi)外對(duì)地下結(jié)構(gòu)的抗震研究大多集中在地鐵車(chē)站和山嶺陸地隧道，而對(duì)水壓力影響下海底隧道襯砌結(jié)構(gòu)動(dòng)力響應(yīng)的研究則相對(duì)較少.

總的來(lái)說(shuō)，由于海底隧道建造在海底，屬于一種特殊的地下結(jié)構(gòu)，且它處于非常復(fù)雜的工程地質(zhì)和水文地質(zhì)環(huán)境中，存在很多諸如斷裂構(gòu)造、破碎帶等不確定因素，建成后難以更改，一旦發(fā)生地震破壞，其后果是難以估量的. 因此，對(duì)海底隧道地震時(shí)的動(dòng)力響應(yīng)進(jìn)行研究就顯得尤為重要. 對(duì)于傳統(tǒng)陸上地下線型隧道結(jié)構(gòu)，在研究軸向傳播的地震動(dòng)作用下垂直于軸線的橫向振動(dòng)時(shí)，一般只研究了軸力、內(nèi)外阻尼、轉(zhuǎn)動(dòng)慣性、橫截面剪切變形等一種或幾種因素的影響. 而對(duì)于海底隧道這種線型結(jié)構(gòu)所受到的高水壓力對(duì)結(jié)構(gòu)橫向振動(dòng)的影響研究，主要是通過(guò)一些有限元軟件進(jìn)行數(shù)值模擬，目前還缺乏系統(tǒng)的理論研究，高水壓對(duì)于隧道結(jié)構(gòu)的地震響應(yīng)的影響還不清楚. 鑒于此，本文將基于土與結(jié)構(gòu)的相互作用，結(jié)合結(jié)構(gòu)動(dòng)力學(xué)和彈性力學(xué)相關(guān)知識(shí)，研究高水壓對(duì)于結(jié)構(gòu)自振頻率及橫向振動(dòng)位移的影響. 首先是利用平衡條件和微分關(guān)系推導(dǎo)出僅由彎曲引起的結(jié)構(gòu)自由振動(dòng)方程和強(qiáng)迫振動(dòng)方程. 之后利用分離變量法和輸入地震振動(dòng)形式分別求得了含水壓影響的結(jié)構(gòu)自振頻率和結(jié)構(gòu)振動(dòng)位移的表達(dá)式，并結(jié)合算例研究了地震作用下水壓對(duì)海底隧道結(jié)構(gòu)自振頻率和振動(dòng)位移的影響.

1振動(dòng)方程

1.1模型選取

溫克爾(Winkler，1867)地基梁模型(見(jiàn)圖1)由捷克工程師溫克爾(Winker)提出，是—種常用的線彈性地基模型. 它的基本假設(shè)是土介質(zhì)表面任一點(diǎn)豎向位移w=w(x,t)與該點(diǎn)所受的地基反力p(x,t)成正比，而與其他點(diǎn)地基反力無(wú)關(guān)，其表達(dá)式如下：

p(x,t)=k0w(x,t)

(1)

隧道為線型結(jié)構(gòu)，重點(diǎn)研究沿結(jié)構(gòu)長(zhǎng)度方向的結(jié)構(gòu)振動(dòng)和受力情況，故將式(1)中的單位面積上地基基床系數(shù)k0變換為單位長(zhǎng)度上的地基基床系數(shù)k. 利用關(guān)系式k=k0b，則(1)式變?yōu)椋?/p>

p(x,t)=kw(x,t)

(2)

溫克爾模型簡(jiǎn)單且為線性，在地下結(jié)構(gòu)的設(shè)計(jì)和分析中，常選用此模型來(lái)模擬土與結(jié)構(gòu)的相互作用.

1.2水壓下的軸力

目前，海底隧道水壓力取值通常采用折減系數(shù)法、全水頭法、可能最大水頭值法3種估算方法. 這3種方法應(yīng)用最廣泛的是折減系數(shù)法，而數(shù)值計(jì)算法和解析法有時(shí)也有應(yīng)用. 本文在考慮高水壓影響時(shí)，采用全水頭法計(jì)算水壓力，并將水壓力按照靜水壓力考慮. 對(duì)于海底隧道來(lái)說(shuō)，水壓力是突出問(wèn)題，它不僅影響隧道結(jié)構(gòu)的橫向受力，還會(huì)引起隧道結(jié)構(gòu)縱向軸力變化，而隧道是一個(gè)線型結(jié)構(gòu)，兩端可以看作鉸結(jié)在地基上，故可以將其按照平面應(yīng)變問(wèn)題進(jìn)行考慮. 假設(shè)水壓力用qw表示，則根據(jù)彈性力學(xué)知識(shí)：

(3)

根據(jù)廣義胡克定律由εz=0可得:

(4)

則水壓力引起的軸力可按下式近似求得:

(5)

若設(shè)襯砌結(jié)構(gòu)所受初始軸壓為N0，則結(jié)構(gòu)軸向所受總壓力(以壓力為正)為：

N=N0+Nw=N0+2υπb2qw

(6)

式中，σr、σθ、σz分別為襯砌結(jié)構(gòu)徑向、環(huán)向、軸向的應(yīng)力，a、b分別為襯砌結(jié)構(gòu)的內(nèi)半徑和外半徑.

1.3振動(dòng)方程推導(dǎo)

本文在建立梁橫向振動(dòng)方程時(shí)作了如下幾點(diǎn)假設(shè)：

1)選用線彈性溫克爾地基梁模型；

2)不考慮錯(cuò)縫拼裝和螺栓連接效應(yīng)的影響，結(jié)構(gòu)按等剛度原則考慮；

3)對(duì)于淺埋結(jié)構(gòu)，一般認(rèn)為埋深越大，結(jié)構(gòu)受到圍巖的約束就越強(qiáng)，結(jié)構(gòu)的動(dòng)力響應(yīng)越小，本文在推導(dǎo)過(guò)程中不考慮埋深對(duì)結(jié)構(gòu)動(dòng)力響應(yīng)的這種影響；

4)不考慮梁在運(yùn)動(dòng)過(guò)程中剪切變形和轉(zhuǎn)動(dòng)慣量的影響，即認(rèn)為梁撓曲線的斜率僅由彎矩所引起；

5)忽略阻尼影響.

設(shè)有如圖2所示的地鐵隧道簡(jiǎn)化模型，沿軸向x方向的抗彎剛度為EI，單位長(zhǎng)度質(zhì)量為m(x)=ρA，作用在模型上的橫向荷載包括地基反力和慣性力，模型橫向位移用w(x,t)表示.

如圖3所示，取梁上任一截面x處的微段dx為隔離體，結(jié)構(gòu)軸線在彎矩作用下將會(huì)發(fā)生彎曲，微段右端面與左端面之間存在一個(gè)微小的傾斜角度dφ=dx?2w/?x2. 因?yàn)?2w/?x2=?θ/?x，則有dφ=dx?2w/?x2=dx?θ/?x(θ為彎矩產(chǎn)生的撓曲線斜率). 該微段上除了作用在兩個(gè)截面上的彎矩M、剪力Q、軸力N和地基反力p(x,t)外，根據(jù)達(dá)朗貝爾原理，還有假設(shè)的慣性力fI(x).

根據(jù)豎向力的平衡則有：

(7)

即

(8)

由于不考慮梁在運(yùn)動(dòng)過(guò)程中剪切變形和轉(zhuǎn)動(dòng)慣量的影響，梁撓曲線的斜率僅由彎矩造成. 根據(jù)微分關(guān)系有：

(9)

將式(9)代入式(8)得考慮軸力影響的梁自由振動(dòng)方程：

(10)

將式(2)、式(6)代入上式即可得到考慮高水壓影響的梁自由振動(dòng)方程：

(11)

式中：E——襯砌結(jié)構(gòu)的彈性模量(N/m2)；

Iz——襯砌結(jié)構(gòu)的橫截面慣性矩(m4)；

υ——襯砌結(jié)構(gòu)的泊松比.

事實(shí)上，已有文獻(xiàn)推導(dǎo)了考慮包括軸力、剪切變形、內(nèi)外阻尼、橫截面轉(zhuǎn)動(dòng)慣性以及地基參數(shù)等因素的更為一般的梁振動(dòng)方程，本文只研究包含水壓力影響的僅由彎矩引起的比較簡(jiǎn)單的情況.

若梁作強(qiáng)迫振動(dòng)時(shí)，wg(x,t)為自由場(chǎng)地基側(cè)向位移，則文克爾地基表達(dá)式變?yōu)椋?/p>

p(x,t)=k[w(x,t)-wg(x,t)]

(12)

此時(shí)梁強(qiáng)迫振動(dòng)方程變?yōu)椋?/p>

(13)

2振動(dòng)方程的求解

2.1自振頻率的求解

為求解式(11)表示的梁自由振動(dòng)方程的振型和自振頻率，本文采用分離變量法，假設(shè)解具有如下形式：

w(x,t)=W(x)T(t)

(14)

代入式(11)整理得：

(15)

式(15)等號(hào)兩端分別是x和t的函數(shù)，對(duì)于任意獨(dú)立的x和t，式(15)只能等于一個(gè)常數(shù)，此常數(shù)用λ4表示. 則有：

(16)

(17)

式(16)為常見(jiàn)的無(wú)阻尼單自由度振動(dòng)方程，其解為：

T(t)=Acosωt+Bsinωt=asin(ωt+φ)

W(x)=C1cosαx+C2sinαx+C3coshβx+C4sinhβx

(19)

利用邊界條件：

W(0)=W″(0)=0W(L)=W″(L)=0

可得振型函數(shù)的解為 ：

(20)

(21)

2.2強(qiáng)迫振動(dòng)方程的求解

對(duì)于方程(13)，其解包括兩部分如下：

w(x,t)=wc(x,t)+wp(x,t)

(22)

式中，右邊第一項(xiàng)代表自由振動(dòng)解，反映了初始條件對(duì)結(jié)構(gòu)體系振動(dòng)的作用；第二項(xiàng)代表強(qiáng)迫振動(dòng)特解，表示外來(lái)輸入對(duì)結(jié)構(gòu)體系振動(dòng)的作用. 因?yàn)樽枘嶙饔糜绊懣偸谴嬖?，第一?xiàng)在一段時(shí)間之后將很快衰退而消失，可以不用考慮，只剩下一個(gè)穩(wěn)定的純強(qiáng)迫振動(dòng)(常稱(chēng)為穩(wěn)態(tài)振動(dòng))，其振動(dòng)特性由外界干擾力的性質(zhì)以及變化規(guī)律決定.

強(qiáng)迫振動(dòng)方程在求解時(shí)常作如下假定：地下結(jié)構(gòu)和地基土是各向同性、均勻的線彈性介質(zhì)；不考慮地震波輸入時(shí)的隨機(jī)性和復(fù)雜性，地震波按剪切波考慮，入射方向與隧道軸向之間夾角為零.

如圖4所示，假定自由場(chǎng)土介質(zhì)某一時(shí)刻的側(cè)向位移為wg(x,t)，幅值為D. 隧道結(jié)構(gòu)側(cè)向位移用w(x,t)表示，它與自由場(chǎng)土介質(zhì)對(duì)應(yīng)時(shí)刻側(cè)向位移形狀相似而幅值不同，兩者相差一個(gè)相互作用系數(shù)R.

地震波以簡(jiǎn)諧波形式在地層中傳播，地層在地震波作用下發(fā)生簡(jiǎn)諧位移，其位移數(shù)學(xué)表達(dá)式假定如下：

地震位移場(chǎng)：

wg(x,t)=Dsin(γx-θt)

(23)

結(jié)構(gòu)的振動(dòng)位移：

w(x,t)=RDsin(γx-ωt)

(24)

式中：θ——地震波頻率；

ω——結(jié)構(gòu)振動(dòng)頻率；

λ——地震波波長(zhǎng).

為了求出結(jié)構(gòu)振動(dòng)位移，應(yīng)首先求出相互作用系數(shù)R. 將式(23)和(24)帶入式(13)得：

[EIγ4-(N0+2υπb2qw)γ2+k-ρAω2]RDsin(γx-

ωt)=kDsin(γx-θt)

(25)

整理可以得到相互作用系數(shù)為：

(26)

把式(26)帶入式(24)可以得到結(jié)構(gòu)振動(dòng)位移表達(dá)式為：

(27)

將式(24)和式(27)進(jìn)行比較容易看出，外荷載頻率與結(jié)構(gòu)振動(dòng)頻率相等，即θ=ω，則相互作用系數(shù)表達(dá)式為：

(28)

把式(28)帶入式(24)，得到結(jié)構(gòu)振動(dòng)位移表達(dá)式如下：

(29)

3算例

某盾構(gòu)圓形海底隧道，最低點(diǎn)位于水面下130m，水壓約1.3MPa，隧道長(zhǎng)度L=200m，屬大盾構(gòu)，外徑d1=16.7m，內(nèi)徑d2=15.7m，混凝土襯砌密度ρ=2 600kg/m3，彈性模量E=3.0×104MPa，泊松比υ=0.3，初始軸向壓力N0=50 000kN，取地基基床系數(shù)k0=2.0×107N/m3. 地震強(qiáng)度PGV=0.3m/s(相當(dāng)于8度地震)，地震位移幅值D=2.9cm，地震波長(zhǎng)λ=100m，波速v=180m/s，頻率θ=11.31rad/s，地震動(dòng)位移函數(shù)為wg=0.029×sin(0.062 8x-11.31t).

為了考察水壓對(duì)結(jié)構(gòu)自振頻率的影響，結(jié)構(gòu)單位長(zhǎng)度質(zhì)量、抗彎剛度、地震動(dòng)頻率均保持不變，根據(jù)式(21)可得到不考慮水壓影響和考慮水壓為2MPa影響的結(jié)構(gòu)前10階的自振頻率以及不同水壓下結(jié)構(gòu)前三階的自振頻率，具體計(jì)算結(jié)果如圖5和圖6所示.

從圖5可以看出，水壓對(duì)結(jié)構(gòu)低階自振頻率與高階自振頻率的影響差別很大. 考慮2 MPa水壓時(shí)，前五階自振頻率下降幅度為5.99%～61.56%，六階自振頻率以后下降幅度在3.38%以下. 因此，水壓對(duì)結(jié)構(gòu)低階頻率影響較大，對(duì)結(jié)構(gòu)的高階頻率影響幾乎可以忽略. 根據(jù)圖6可以看出，結(jié)構(gòu)的前三階的自振頻率隨著水壓的逐漸增大而近似線性減小. 并且從圖中可以看出，水壓增大到一定值后，結(jié)構(gòu)的自振頻率為0，此時(shí)結(jié)構(gòu)不振動(dòng)，這也說(shuō)明了水壓對(duì)結(jié)構(gòu)有一定的限制作用，阻礙著結(jié)構(gòu)的振動(dòng). 但從圖5中可以看出，結(jié)構(gòu)的自振頻率為0時(shí)，所受臨界水壓須達(dá)到2.5 MPa以上，目前修建的海底隧道所受水壓最高在2 MPa左右，因此不會(huì)出現(xiàn)自振頻率為0的情況.

在研究水壓對(duì)襯砌結(jié)構(gòu)振動(dòng)位移的影響時(shí)，為使研究具有針對(duì)性，取結(jié)構(gòu)振動(dòng)位移最大的時(shí)刻的位移幅值進(jìn)行研究. 圖7和圖8分別給出了不同水壓力下，結(jié)構(gòu)振動(dòng)位移幅值及地基與結(jié)構(gòu)相互作用系數(shù)的變化情況. 從圖7中可以看出，結(jié)構(gòu)的振動(dòng)位移幅值隨著水壓的增大而逐漸增大. 這可能是因?yàn)殡S著水壓力逐漸增大時(shí)，隧道結(jié)構(gòu)的自振頻率會(huì)減小，表現(xiàn)為結(jié)構(gòu)剛度逐漸降低，從而使結(jié)構(gòu)的振動(dòng)位移幅值增大. 從本例的計(jì)算結(jié)果可以看出，水壓對(duì)結(jié)構(gòu)振動(dòng)位移幅值的變化已超過(guò)8%，水壓力對(duì)結(jié)構(gòu)振動(dòng)位移幅值的影響不可忽略. 從圖8中可以看出，水壓對(duì)相互作用系數(shù)的影響也很顯著. 隨著水壓的逐漸增大，相互作用系數(shù)也在逐漸增大. 這說(shuō)明由于水壓的存在，地震動(dòng)對(duì)結(jié)構(gòu)的動(dòng)力響應(yīng)的影響將加強(qiáng)，這也驗(yàn)證了前述水壓對(duì)結(jié)構(gòu)振動(dòng)位移幅值的增大效應(yīng).

4結(jié)論

本文采用溫克爾地基模型，結(jié)合結(jié)構(gòu)動(dòng)力學(xué)和彈性力學(xué)相關(guān)知識(shí)推導(dǎo)了僅由彎曲引起的梁自由振動(dòng)方程和強(qiáng)迫振動(dòng)方程，并求出了含高水壓影響的結(jié)構(gòu)的自振頻率表達(dá)式以及結(jié)構(gòu)的位移表達(dá)式. 通過(guò)算例得出水壓對(duì)結(jié)構(gòu)的低階自振頻率影響顯著，對(duì)高階自振頻率的影響可以忽略，并且結(jié)構(gòu)的位移

幅值隨水壓增大而增大. 因此，對(duì)處于高水壓情況下的隧道來(lái)說(shuō)，水壓影響應(yīng)該考慮.

參考文獻(xiàn)：

[1]譚忠盛,王夢(mèng)恕,羅時(shí)祥.瓊州海峽鐵路隧道方案初步比選分析[J].中國(guó)工程學(xué),2009,11(7):39-44

[2]石新棟,皇甫明,譚忠盛,等.跨瓊州海峽隧道方案的探討[J].隧道建設(shè)，2010，30(6)：625-628

[3]Dowing C H， Rozen A. Damage to rock tunnels from earthquake shaking[J]. Journal of the Geotechnical Engineering Division， 1978, 104(2):175-191

[4]張維. 日本兵庫(kù)縣南部地震災(zāi)區(qū)現(xiàn)場(chǎng)震害考察[J].工程抗震, 1995(1)：2

[5]楊光,沈繁鑾.日本阪神地震災(zāi)害的一些調(diào)查統(tǒng)計(jì)數(shù)據(jù)[J].華南地震, 2005, 25(1)：83-86

[6]Stephen N G．The second spectrum of Timoshenko beam theory—further assessment[J].Journal of Sound and Vibration, 2006, 292(1-2) :372-389

[7]Majkut L. Free and forced vibrations of Timoshenko beams described by single difference equation[J]. Journal of Theoretical and Applied Mechanics, 2009,47(1): 193-210

[8]Bogacz R, Czyczula W. Response of beam on visco-elastic foundation to moving distributed load[J]. Journal of Theoretical and Applied Mechanics, 2008, 46(4): 763-775

[9]徐英乾,戚承志,沙曼.阻尼對(duì)地下區(qū)間隧道動(dòng)力響應(yīng)的影響[J].巖土力學(xué), 2012(2):196-200

[10]鐵摩辛柯，古地爾.徐芝綸譯.彈性理論[M].北京：人民教育出版社，1964:12

[11]戚承志, 陳燦壽, 陳國(guó)興, 等. 地基彈性參數(shù)對(duì)于長(zhǎng)型大截面地下結(jié)構(gòu)橫向自由振動(dòng)的影響[J]. 防災(zāi)減災(zāi)工程學(xué)報(bào), 2011(6): 622-626

[12]陳燦壽，戚承志，陳劍杰，等.長(zhǎng)型大截面地下結(jié)構(gòu)的橫向自由振動(dòng)方程[J].世界地震工程, 2011, 27(1): 90-95

[13]韋洪剛， 楊秀仁，戚承志，等. 地下結(jié)構(gòu)縱向地震響應(yīng)分析[A]∥全國(guó)第一屆防災(zāi)減災(zāi)工程學(xué)術(shù)研討會(huì)[C], 2011

[14]韋洪剛，楊秀仁，戚承志. 長(zhǎng)型地下結(jié)構(gòu)在地震作用下的橫向響應(yīng)分析 [D].北京：北京建筑大學(xué), 2012

[15]彭俊生，羅永坤，彭地.結(jié)構(gòu)動(dòng)力學(xué)、抗震計(jì)算與SAP2000應(yīng)用[M].成都：西南交通大學(xué)出版社，2006:90-97

[責(zé)任編輯：佟啟巾]

Research on Seismic Dynamic Response of Submarine Tunnel

Liu Tiantian,Qi Chengzhi,Chen Haoxiang,Yang Jinrui

(School of Civil and Traffic Engineering, Beijing Higher Institution Engineering Research Center of Structure and New Materials,Beijing University of Civil Engineering and Architecture, Beijing 100044)

Abstract:Seismic dynamic response of submarine tunnel under high water pressure is a problem in need of further investigation. Based on Winkler’s foundation model, the structure vibration equation with consideration of only bending is derived by using dynamic equilibrium condition and differential relations. Using the method of variable separation, the vibration mode and natural vibration frequencies of submarine tunnel are obtained. And the interaction coefficient between submarine tunnel and foundation and structural displacement is also determined. Numerical example is given and the effect of water pressure on structural natural vibration frequencies and structural displacement and interaction coefficient is studied. Investigation shows that structural seismic dynamic response under effect of high water pressure should not be neglected.

Key words:submarine tunnel; water pressure; Winkler foundation model; natural frequency; vibration displacement

中圖分類(lèi)號(hào)：U459.9; TU435

文獻(xiàn)標(biāo)志碼：A

作者簡(jiǎn)介：劉天添(1988—)，男，碩士研究生，研究方向：結(jié)構(gòu)抗震、巖土力學(xué).

基金項(xiàng)目：國(guó)家自然科學(xué)基金項(xiàng)目(51174012，51478027);北京市屬高等學(xué)校創(chuàng)新團(tuán)隊(duì)建設(shè)與教師職業(yè)發(fā)展計(jì)劃項(xiàng)目(IDHT20130512);國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973)項(xiàng)目(2015CB0578005)

收稿日期：2015-09-06

文章編號(hào)：1004-6011(2016)01-0046-06