某電站取水隧洞與熱水回流溝銜接段抗震分析

趙 杰，王桂萱，程 翔

(大連大學土木工程技術研究與開發中心，遼寧大連 116622)

某電站取水隧洞與熱水回流溝銜接段抗震分析

趙 杰，王桂萱，程 翔

(大連大學土木工程技術研究與開發中心，遼寧大連 116622)

運用動力時程反應分析方法，利用非線性有限差分軟件FLAC3D對某電站取水隧洞與熱水回流溝銜接段結構進行抗震分析。針對場地地質特征和回填條件分析結構地震響應，通過編制FLAC3D程序并結合五點公式計算襯砌結構的內力，給出了襯砌內力包絡圖。數值分析結果表明:在相同的荷載效應組合情況下，溫度的改變對內力的影響較為顯著，溫度荷載不僅顯著改變了隧洞和熱水回流溝銜接段的內力大小，而且改變了某些部位的內力方向;銜接段采用圓形接口優于采用方形接口。

取水隧洞;熱水回流溝;抗震分析;內力包絡圖;FLAC3D軟件

近年來，地下結構在能源、交通、通訊、城市建設和國防工程中獲得廣泛應用，地下結構的重要性也日益顯現。隨著地下結構震害的頻繁發生，地下結構抗震問題日益受到世界各國地震工作者的高度重視。研究地下結構抗震性能的主要途徑有原型觀測、模型試驗和數值模擬［1］。由于問題極其復雜，目前還沒有哪一種手段能夠完全實現對地下結構動力反應進行全面而真實的解釋和模擬。傳統的地下結構抗震分析方法主要包括地震系數法、反應位移法、圍巖應變傳遞法等［2-4］。這些方法實際上都是擬靜力法［5-7］，具有計算簡便的優點，但假定條件較多，無法精確考慮材料非線性非均質性和復雜邊界變化等因素的影響，難以考慮地震作用下土體與結構之間的相互作用。隨著計算機和計算理論的發展，以動力有限元法為代表的數值方法應運而生，它彌補了以上方法的缺陷，為各種復雜情況下地下結構抗震特性的全面深入研究提供了重要手段。本文在給定的地質和回填條件下開展某電站取水隧洞和熱水回流溝銜接段三維地震響應分析，編制三維有限差分FLAC3D程序并結合五點公式計算襯砌結構的內力，繪出襯砌實體單元內力包絡圖，在此基礎上研究銜接段混凝土結構的內力變化規律，為地下結構抗震設計提供依據。

1 工程概況

某電站機組通過兩條直徑為5.5 m的取水隧洞將水引至泵房，兩隧洞的中心距為25.00～27.40m。隧洞圍巖為Ⅴ類片麻巖。根據工程需要，在隧洞出口與泵房前池閘門井間增設一道熱水回流溝(圖1)。熱水通過熱水回流溝，經隧洞頂部預留孔進入取水隧洞，隨后進入泵房前池閘門井。銜接段采用碎石回填。鑒于取水隧洞工程的重要性，需專門對該段取水隧洞進行抗震分析。

圖1 熱水回流溝(單位:mm)

2 銜接段抗震分析模型的建立

2.1 邊界條件的設置

為準確模擬實際場地中地震波的傳播過程，需在模型中設置人工邊界。黏性邊界是發展最早且簡單實用的人工邊界［8］。黏性邊界最早由Kuhlmeyer和Lysmer［9-10］提出，是指在邊界上施加法向和切向與邊界無關的獨立阻尼器，由阻尼器提供法向和切向黏性阻力，以此消除地震波在人為設置的模型邊界上的反射效應，吸收或消耗傳往邊界外的波動能量，能比較真實地反映地震波的傳播過程。黏性邊界的黏性阻力計算公式為

式中:tn、ts分別為法向和切向黏性阻力;ρ為材料密度;Cp、Cs分別為縱波和橫波速度;vn、vs分別為邊界上速度法向和切向分量。

在模型各側面的邊界條件必須考慮沒有地面結構的自由場。通過在模型四周生成二維和一維網格的方法來實現這種自由場的邊界條件，柱體網格的側邊界通過阻尼器與自由場網格進行耦合，將自由場網格的不平衡力施加到主體網格的邊界上。

在FLAC3D動力分析中，為準確模擬地震波的傳播過程，空間單元尺寸必須小于輸入波最高頻率成分所對應波長的1/10～1/8。本次計算模型最大的網格寬度為9.25 m，可以考慮地震波主要頻譜成分的影響。同時，數值模擬中采用局部阻尼，局部阻尼系數取0.157。

2.2 模型的建立

三維動力分析模型與靜力模型計算范圍一致，水平方向取5倍隧洞洞徑，基巖深度自隧洞底部向下取50 m。模型寬135 m，長76 m，高82 m。計算時模型底部設為黏性邊界，兩側采用能量透射邊界。模型采用六面體單元，巖體本構模型采用 mohrcoulomb彈塑性模型，采用null模型模擬隧洞的開挖，采用CABLE單元模擬系統錨桿和土釘。取水隧洞與熱水回流溝銜接段結構模型見圖2，三維靜力和動力分析模型見圖3。

圖2 取水隧洞與熱水回流溝銜接段結構模型

圖3 三維靜力和動力分析模型

2.3 地震波的選取

場地基巖輸入采用該電站地震安全評價報告中的場址地震波，地震動荷載持時25 s，地震波時程曲線見圖4。根據該電站地震安全性評價報告批準文件，場地基巖50a超越概率10%的地震動荷載加速度峰值為0.1g，豎向加速度取水平向加速度的2/3。

2.4 計算參數

抗震分析計算參數根據該電站施工圖設計階段巖土工程詳勘報告取值，列于表1中。根據該電站取水泵房溫排水回熱系統試驗研究報告，內外溫差取6℃。

表1 銜接段抗震分析計算參數

圖4 場址地震波時程曲線

2.5 計算方案

具體計算時，綜合考慮隧洞洞口的開挖回填以及熱水回流溝的施工，建立整體三維分析模型。根據給定的設計方案開展抗震分析，分析地震動荷載作用下溫度變化時結構的內力分布，為設計提供依據。利用FLAC3D計算分析地下結構截面上的內力分布規律時，如采用三維實體單元，截面形狀雖能滿足工程需求，但實體單元只給出了單元對應的位移及應力結果，無法直接得出相應截面上的內力(如彎矩、剪力和軸力)，因此，需要對FLAC3D進行二次開發并結合五點公式計算襯砌結構的內力。式(3)～(5)分別為彎矩、剪力和軸力計算公式。式中應力一律用同一標號來代替，即軸力計算公式中σ1～σ5為所求軸力方向的應力，彎矩計算公式中σ1～σ5為所求彎矩方向的應力，剪力計算公式中σ1～σ5為所求剪力方向的應力。

式中:M、F、N分別為彎矩、剪力和軸力;σ1、σ2、σ3、σ4、σ5分別為截面高度方向第1個節點、第2個節點、第3個節點、第4個節點、第5個節點的應力;b為截面寬度;h為截面高度。

選取兩個控制截面(圖5)，獲取不同位置處的結構內力及變形變化。控制截面1和2的控制點布置見圖6。

圖5 控制截面位置示意圖

圖6 控制截面控制點布置(單位:m)

3 結構抗震分析

在給定的設計地震動荷載作用下，進行銜接段接口處的地震響應分析，在此基礎上繪制襯砌內力包絡圖(圖7)，研究結構的內力分布。

3.1 抗震計算結果

a.工況1:設計地震作用，不考慮溫度荷載。從圖7可以看出，控制截面1最大彎矩絕對值為267.8kN·m，最大軸力絕對值為4833.0kN，最大剪力絕對值為693.6 kN;控制截面2最大彎矩絕對值為100.4 kN·m，最大軸力絕對值為1 247.4 kN，最大剪力絕對值為507.6 kN。

b.工況2:設計地震作用+溫度荷載(內外溫差6℃)。從圖7可以看出，控制截面1最大彎矩絕對值為485.4 kN·m，最大軸力絕對值為5 348.2 kN，最大剪力絕對值為823.8 kN;控制截面2最大彎矩絕對值為314.9kN·m，最大軸力絕對值為1557.5 kN，最大剪力絕對值為778.2 kN。

對比可以看出:在相同的荷載效應組合情況下，溫度的改變對內力的影響較為顯著，隨著隧洞內水溫的升高，銜接段內力顯著增大;溫度荷載不僅顯著改變了隧洞和熱水回流溝銜接段的內力大小，而且改變了某些部位的內力方向。

圖7 襯砌內力包絡圖

3.2 接口形狀優化

在上述抗震分析中，取水隧洞與熱水回流溝銜接口形狀為方形，在優化設計分析中將接口形狀改為圓形，將圓形接口的內力計算結果與方形接口的內力計算結果進行比較，為設計提供參考。計算結果見表2，從表2可以看出，圓形接口的受力優于方形接口，方形接口的彎矩值比圓形接口的彎矩值高約15%，說明圓形接口更合理。

4 結語

利用三維有限差分軟件FLAC3D開展某電站取水隧洞與熱水回流溝銜接段地震響應分析，闡述了利用FLAC3D進行地下結構動力分析的關鍵性問題，針對場地地質特征和回填條件進行結構地震響應分析，在此基礎上給出了襯砌內力包絡圖。結果表明:①在相同的荷載效應組合情況下，溫度的改變對內力的影響較為顯著;②溫度荷載不僅顯著改變了隧洞和熱水回流溝銜接段的內力大小，而且改變了某些部位的內力方向;③圓形接口形狀優于方形接口形狀，方形接口的彎矩值比圓形接口的彎矩值高約15%。該成果對類似地下工程抗震設計具有一定的參考價值。

表2 兩種不同接口形狀銜接段內力最大值比較

［1］劉晶波，李彬.地鐵地下結構抗震分析及設計中的幾個關鍵問題［J］.土木工程學報，2006，39(6):106-110.(LIU Jingbo，LI Bin.Issues on the seismic analysis and design of subway structures［J］.China Civil EngineeringJournal，2006，39(6):106-110.(in Chinese))

［2］趙杰，王桂萱，裴強，等.核電廠取水隧洞抗震分析［J］.世界地震工程，2009，25(3):135-139.(ZHAO Jie，WANG Guixuan，PEI Qiang，et al.Seismic analysis of an intake tunnel for a nuclear power plant［J］.World Earthquake Engineering，2009，25(3):135-139.(in Chinese))

［3］趙杰，王桂萱，裴強，等.某核電廠取水隧洞不同設計方案抗震性能分析［J］.防災減災學報，2010，26(1):31-35.(ZHAO Jie，WANG Guixuan，PEI Qiang，et al.Seismic analysis of nuclear power project for intake tunnel under different design scheme［J］.Journal of Disaster Prevention and Reduction，2010，26(1):31-35.(in Chinese))

［4］趙杰，王桂萱，裴強，等.影響核電廠取水隧洞抗震性能的幾個因素分析［J］.遼寧工程技術大學學報，2011，30(1):46-49.(ZHAO Jie，WANG Guixuan，PEI Qiang，et al.Seismic analysis of nuclear power project for intake tunnel［J］.Journalof Liaoning Technical University，2011，30(1):46-49.(in Chinese))

［5］SHUKLA D K，RIZZO P C，STEPHENSON D E.Earthquake load analysis of tunnels and shafts［C］//Proceeding of the Seventh World Conference on Earthquake Engineering.Michigan:［s.n.］，1980:20-28.

［6］孫鈞，侯學淵.地下結構:下冊［M］.北京:科學出版社，1988:783-794.

［7］鄭永來，楊林德，李文藝，等.地下結構抗震［M］.上海:同濟大學出版社，2005.

［8］FLAC-3D(fast lagrangian analysis of continua in 3 dimensions)users manual［R］.Minneapolis，Minnesota:Itasca Consulting Group Inc.，2006.

［9］KUHLMEYER R L，LYSMER J.Finite element method accuracy for wave propagation problems［J］.Journal of the Soil Mechanics and Foundations Division，ASCE，1973，99:421-427.

［10］LYSMER J，KUHLMEYER R L.Finite dynamic model for infinite media［J］.Journal of the Soil Mechanics and Foundations Division，ASCE，1969，95:859-877.

Seismic analysis of the joint section between intake tunnel and hot-water circumfluence groove for a hydropower

station

ZHAO Jie，WANG Guixuan，CHENG Xiang(Research and Development Center of Civil Engineering Technology，Dalian University，Dalian116622，China)

By using the dynamic time-history analysis method，seismic analysis was performed on the joint section between the intake tunnel and the hot-water circumfluence groove for a hydropower station by FLAC3D software.Aiming at the geological characteristics and the backfill condition of the site，analysis of seismic response for the joint section was carried on，and the internal forces of the lining structure were calculated by combining secondary development of FLAC3D with the five-point formula.Based on it，the internal force envelop diagram for the lining structure was presented.Numerical simulation shows that，under identical load effect combination，the effect of temperature on internal forces is remarkable:the variety on temperature load not only significantly changes the value but also changes the orientations of the internal forces for certain parts.Meanwhile，it can be concluded that the circular interface is superior to the square interface for the joint section.

intake tunnel;hot-water circumfluence groove;seismic analysis;internal force envelop diagram;FLAC3D software

TV671;P315.9

A

1006-7647(2013)02-0073-05

10.3880/j.issn.1006-7647.2013.02.016

遼寧省教育廳科學技術研究項目(L2011216);大連市科學技術基金(2010J21DW013)

趙杰(1980—)，男，河北邢臺人，講師，博士，主要從事地下結構穩定與工程抗震研究。E-mail:zhaojie_gd@163.com

2012-05-28 編輯:駱超)