蔣莉， 王磊， 周曉嵐， 許新勇

(1.華北水利水電大學 水利學院，河南 鄭州 450046； 2.水資源高效利用與保障工程河南省協(xié)同創(chuàng)新中心，河南 鄭州 450046；3.河南省水工結構安全工程技術研究中心，河南 鄭州 450046； 4.黃河水利職業(yè)技術學院，河南 開封 475001)

我國西南地區(qū)水能資源豐富，水電作為清潔能源得到了大力開發(fā)。在水電站建設中，隨著施工開挖機械不斷改進和施工技術的不斷提高，陸續(xù)出現(xiàn)了一大批地下廠房。地下廠房抗震性能較好，受到水利行業(yè)眾多專家的青睞。西南地區(qū)屬于地震高烈度區(qū)，地下廠房在地震動荷載作用下的響應特征及受力狀況成為不容忽視的問題。

目前，水電站廠房抗震領域研究已較為成熟，其中：喻虎圻等[1]運用三維有限元分析技術，對廠房地震持時過程及結構抗震安全性進行了分析，發(fā)現(xiàn)地基阻尼可以有效降低廠房地震響應；朱少坤等[2]運用耐震時程法合成耐震加速度曲線，研究了廠房在不同地震強度作用下的結構響應情況，得出耐震時程法可避免調幅計算，提高了計算效率；谷音等[3]借用經典Lamb算例與成層半無限空間算例，進行了不同邊界條件下的結構響應分析，得出黏彈性邊界施加簡便且計算精度高的結論；韓濤[4]運用黏彈性邊界以及脈沖波方案完成了廠房振動時域分析，根據共振校核和動力系數(shù)復核驗證了地下廠房總體設計的合理性；范書立等[5]通過多個模型試驗對不同變量進行對比分析，探究了龍開口水電站大壩動力破壞特征，找到了不同變量對壩體開裂形式和裂縫發(fā)展的影響；付杰等[6]通過研究黏彈性邊界條件下的拱壩結構流固耦合作用機理，找到了降低拱壩響應的兩個因素，即水可壓縮性和無限地基輻射阻尼；歐陽金惠等[7]對抽水蓄能電站進行現(xiàn)場測試與數(shù)值計算的對比分析，探明了主振源及振動原因并提供了減振措施。

雖然許多學者對于黏彈性邊界和地震動有了很深入的研究，但是將地下廠房水流系統(tǒng)、發(fā)電系統(tǒng)及支撐結構作為一個整體，在黏彈性邊界條件下對其結構抗震性能的研究相對較少。本文以瀾滄江干流河段某大型水電站為例，建立包括水流系統(tǒng)、發(fā)電系統(tǒng)及支撐結構的數(shù)值仿真模型，采用黏彈性邊界進行人工地震波條件下的結構動力特性研究，為水電站地下廠房建設及結構性態(tài)規(guī)律研究提供工程借鑒與指導建議。

1 研究理論及方法

在應力、應變問題上，結構的變形與彈性本構之間的相互影響使結構荷載的累加發(fā)生變化。而這種累加變化正適合于復雜的廠房結構應力狀態(tài)[8]，則可以證明：

εp=[1+ks]-1[kD](ε-εz)，

(1)

εs=(I+k*C)-1k*W(ε-εz)。

(2)

水電站廠房樓板采用殼體結構，則：

式中：s為單位矩陣;k、k*分別為彈性系數(shù)、塑性系數(shù)；I為與C同階的矩陣；C為結構的彈性矩陣；D、W分別為彈性單位矩陣、塑性單位矩陣；η為地基黏滯系數(shù)；參數(shù)β的取值范圍為[0.0,0.2]，根據具體情況而定；Δt為迭代的時間增量；ε、εp、εs、εz分別代表總應變向量、彈性應變向量、塑性應變向量、屈服極限應變向量；Ex、Ey彈性勢能在x、y方向的分量；Gxy為重力在x、y方向上的合力；Sx、Sy、Sxy為彈性勢能在x、y、xy方向上的增量；δx、δy、δxy分別為x、y、xy方向上強化系數(shù)的倒數(shù)。

三維人工邊界的彈簧-阻尼元件參數(shù)為：

法向：

(3)

切向：

(4)

式中：KBN為邊界節(jié)點B的法向彈簧剛度值；CBN為該點的法向剛度阻尼值；KBT為該點的切向彈簧剛度值；CBT為該點的切向剛度阻尼值；Ab為節(jié)點B的控制面積；R為結構中心到邊界的距離；Ep為結構彈性模量；G為剪切模量；v為泊松比。

進行廠房地震動計算時，將地基視為彈簧剛度矩陣，以更好地模擬遠域地震對結構的影響。利用MATLAB編程語言完成三維彈簧單元的COMBIN14實現(xiàn)，其實現(xiàn)方法如圖1所示。

圖1 黏彈性邊界實現(xiàn)方法

2 黏彈性邊界效應驗證

為了驗證黏彈性邊界在地下廠房地震動計算中的適用性，建立三維彈性半無限空間模型，尺寸為1.0 m×1.0 m×0.5 m。在模型邊界位置設置黏彈性邊界，并與常規(guī)固定邊界相對比，假定模型材料為各向同性，模型剪切模量G取值為16，泊松比v為0.25，材料密度ρ取為1。其固定邊界和黏彈性邊界數(shù)值模型如圖2所示。

圖2 固定邊界和黏彈性邊界數(shù)值模型示意圖

為研究波在彈性介質中的傳播過程及在邊界位置的能量彌散現(xiàn)象，選取點M(圖2)作為監(jiān)測點進行分析，計算總時長為1 s。圖3(a)為根據脈沖波實測數(shù)據施加在模型中心點O的脈沖荷載時程曲線，圖3(b)為固定邊界下的三維半無限空間模型與黏彈性邊界下的三維半無限空間模型在脈沖波荷載作用下M點的位移時程曲線。由圖3(b)可知：在相同荷載作用下，點M在固定邊界條件下產生周期性震蕩現(xiàn)象，說明固定邊界沒有能量的耗散；在黏彈性邊界條件下，M點位移波動幅值較小，加載結束以后，位移波動極小，說明荷載到達邊界后，未返回彈性介質，能量被邊界幾乎完全吸收，說明黏彈性邊界較好地模擬了波的傳播過程。

圖3 脈沖荷載及M點位移時程曲線勢

3 數(shù)值模型

瀾滄江干流某大型水電站地下廠房樞紐由攔河壩、擋水建筑物、泄水建筑物、水力發(fā)電系統(tǒng)、升船機等組成。地下廠房位于引水系統(tǒng)首部，廠房整體尺寸(長×寬×高)為160 m×118 m×187.5 m，流道全長約110 m，其中主廠房(長×寬×高)為36 m×34 m×42 m，副廠房(長×寬×高)為18 m×34 m×50 m，風罩半徑為13 m，機墩半徑為7 m，單臺機組發(fā)電機層樓板尺寸為34 m×33 m。選取地下廠房為研究對象，建立數(shù)值仿真模型，如圖4所示。

圖4 廠房結構數(shù)值模型圖

模型采用實體單元(Solid 60)、桿單元(Link 180)、殼單元(Shell 63)等。工程忽略溫度、短暫荷載、機械振動等對水電站廠房結構的影響，結構材料參數(shù)見表1。

表1 水電站廠房結構材料參數(shù)表

4 黏彈性邊界條件下的廠房結構動力響應分析

4.1 地震波的選取與輸入

根據地下廠房所在場地類型與地震區(qū)域劃分，該地區(qū)地震設防烈度為Ⅶ度，場地類型為Ⅰ類，其峰值加速度為0.15g，地震持續(xù)時間為20 s，特征周期為0.35 s。根據《水電工程水工建筑物抗震設計規(guī)范》(NB 35047—2015)有關規(guī)定，廠房結構設計反應譜最大代表值為2.25。利用廠房設計反應譜擬合人工地震波，場地及地震波反應譜如圖5(a)所示；其中豎向地震波加速度取水平向的2/3，擬合的水平向地震波歸一化時程曲線如圖5(b)所示。通過積分軟件Seismosignal進行積分，得到整個自由場的速度和位移譜。根據波的傳播原理與黏彈性邊界施加方法，將人工地震波施加在地下廠房模型底部以及四周共5個面的黏彈性人工邊界上，實現(xiàn)地震波輸入，等效替代地震作用。

圖5 反應譜及地震波加速度時程曲線

4.2 水電站廠房時程演變規(guī)律分析

已知地下廠房在有圍巖壓力條件下，進行黏彈性邊界下的人工地震波輸入，利用時程分析法分析廠房整體結構加速度響應，并完成地下廠房地震動計算。圖6分別給出了發(fā)電機層吊物孔左角點、相應的母線層吊物孔左角點、風罩上端內環(huán)某點的3向地震動時程曲線圖。模型順水流方向為縱向，垂直水流方向為橫向，重力方向為豎向。分析圖6得到如下結論：

1)從整體來看，縱向加速度值最大，豎向加速度次之，橫向加速度最小，并且隨地震持時變化，各部位響應規(guī)律大致相同。在14 s左右各結構加速度逐漸變小，最終歸于零點，廠房各部位加速度曲線隨時程變化趨勢與地震波加速度曲線時程變化趨勢相符，且符合地震波相應情況。

2)由圖6(a)可知，在3 s左右，發(fā)電機層出現(xiàn)峰值加速度0.98 m/s2，結構響應出現(xiàn)滯后現(xiàn)象。說明結構受到慣性力及黏彈性邊界消能的共同作用，響應以縱向為主。

3)風罩峰值加速度為0.98 m/s2(圖6(c))，母線層峰值加速度為0.90 m/s2(圖6(b))，風罩結構的峰值加速度高于母線層的。說明風罩受到地震荷載作用下，結構響應明顯，上部結構響應大于下部結構的。

圖6 廠房結構地震響應加速度時程曲線

圖7給出了發(fā)電機層吊物孔左角點、相應的母線層吊物孔左角點和風罩上端內環(huán)某點的3向位移時程曲線，分析圖7得到如下結論：

1)由圖7(a)可知，發(fā)電機層縱向最大位移接近0.06 m，豎向最大位移為0.03 m，橫向最大位移為0.02 m，最大總位移可達0.07 m；較加速度響應峰值，特征點位移峰值分別出現(xiàn)5、10、15 s的滯后。

2)圖7(b)、圖7(c)分別給出了母線層和風罩的位移時程曲線，其位移曲線與發(fā)電機層的相似。母線層和風罩的最大位移分別為0.052、0.056 m，即Δu=0.004 m，分別根據風罩中部及母線層距基礎底部的距離98、92 m，可以得到層間位移角σ=0.000 67(層間位移角σ<1/1 000)，結構抗震安全，滿足設計要求。

圖7 廠房結構位移時程曲線圖

3)從整體來看，結構縱向位移大于豎向和橫向的，并且發(fā)生周期性擺動，縱向位移與豎向和橫向的相差約0.5個周期，說明結構響應以縱向為主，應加強對結構縱向位移的安全檢測。

4.3 不同邊界下廠房地震響應時程分析

選取模型邊界點A、B，其中A點為廠房下游左側邊界頂部角點，B點在蝸殼層與地平線同高的廠房下游邊界棱上，如圖4(a)所示。圖8為廠房邊界特征點A、B在不同邊界條件下的位移曲線，分析圖8得到如下結論：

圖8 不同邊界條件下的特征點A、B的位移時程曲線圖

1)從整體來看，邊界點位移隨著地震持時的加速度變化而發(fā)生相應變化，其趨勢與廠房結構位移發(fā)展規(guī)律大致相同，先增加后回落，且最大位移值都發(fā)生在10 s前后。

2)荷載作用相同時，固定邊界條件下A、B兩點的位移值大于黏彈性邊界條件下的；隨著地震持時增加，固定邊界出現(xiàn)6個波峰，而黏彈性邊界出現(xiàn)4個波峰，說明固定邊界在地震波荷載下，位移反應敏感性大于黏彈性邊界的。

3)在波峰位置，固定邊界位移變化速率高于黏彈性邊界的，說明黏彈性邊界具有一定的吸能作用，應用于地下廠房地震動分析具有一定的適用性及可行性。

4.4 廠房結構承載能力分析

經計算，發(fā)電機層樓板、母線層樓板和風罩的地震動最終時刻應力云圖如圖9所示，其中樓板邊緣缺少的矩形結構為吊物孔。分析圖9得到如下結論：

1)由圖9(a)可知：發(fā)電機層樓板應力最大值為4.17 MPa，處于樓板內緣邊界135°位置；吊物孔處樓板應力值較小，說明吊物孔受到地震作用力較?。欢黧w結構與樓板之間相互作用，在邊沿位置產生較大作用力，應該加強該處結構監(jiān)測。

2)由圖9(b)可知，母線層結構所承受應力大于發(fā)電機層結構的，應力值大于3.06 MPa的區(qū)域占比為整個樓板區(qū)域的5/8左右，說明在地震力的作用下，樓板上層結構所承受應力小于下層結構的。

3)由圖9(c)可知：風罩的應力最大值位于內環(huán)邊緣位置的上部，而內環(huán)的下部結構應力值相對較??；外環(huán)應力值較大的區(qū)域較廣，且最大應力值為7.97 MPa；風罩最大應力與發(fā)電機層樓板最大應力出現(xiàn)的位置相對應。

圖9 廠房結構各部位應力云圖(單位：Pa)

圖10為地下廠房座環(huán)、蝸殼以及尾水管的應力云圖，由圖可知：

圖10 廠房水流系統(tǒng)結構應力云圖(單位：Pa)

1)從整體來看，水流系統(tǒng)結構的應力值大于發(fā)電機系統(tǒng)結構的，說明水流系統(tǒng)結構在承受上部荷載及地震力的作用下，對結構的抗震性能要求更高。

2)固定導葉的最大應力值為167 MPa，且最大應力值位于導葉的外邊緣位置；蝸殼應力最大值為75.4 MPa，位于蝸殼最小斷面隔舌處，應力值較大區(qū)域在蝸殼包角為255°～300°的范圍內；尾水管應力最大值位于肘管段內側及尾水管擴散段與側板相交處，最大值為23.6 MPa。根據《鋼結構設計標準》(GB 50017—2017)，水工金屬結構滿足抗震設計要求。

5 結論

本文研究了黏彈性邊界條件下廠房在地震動荷載作用下的結構動力特性，結論如下：

1)通過邊界算例驗證，在進行地震動計算時，黏彈性邊界對于地震波的傳遞具有一定的耗散作用，符合大地遠域地基對波的傳遞規(guī)律；黏彈性邊界方法可以模擬廠房地震動反應，驗證了黏彈性邊界應用于水電站廠房地震動數(shù)值模擬計算的可行性。

2)通過施加黏彈性邊界和應用地震波輸入方法，對地下廠房整體結構進行地震時程分析，對發(fā)電機層樓板、母線層樓板及風罩結構進行研究，結果發(fā)現(xiàn)：廠房結構以縱向震動為主；廠房結構在地震荷載作用下，上部結構的加速度、位移均大于下部結構的，廠房下部結構應力大于上部結構的。

3)經計算，發(fā)電機層與母線層之間的層間位移角σ<1/1 000，證明結構抗震安全，滿足抗震設計要求；同時可適當增加樓板厚度，并加強樓板與風罩銜接處結構、吊物孔附近結構的配筋，提高結構抗震性能。