徐威 趙臻真 王岳 張雪才

摘 要：高重力壩地震響應分析是確保建筑物設計合理、評價其抗震特性的重要環(huán)節(jié)。從波動理論出發(fā)，基于Python語言對ABAQUS軟件進行二次開發(fā)，得到黏彈性人工邊界及其相應地震動的輸入方法。通過無質量地基模型和有質量黏彈性人工邊界地基模型的對比分析，說明采用人工邊界模型在考慮地基的輻射阻尼效應后，能顯著降低壩體的動力響應，在一定程度上更符合真實地基情況，而采用無質量地基模型進行壩體的抗震分析偏于保守。通過工程實例驗證了本文黏彈性人工邊界及其相應地震動的輸入方法簡單有效，易于實現，同時能保證一定的精度，適合于高混凝土重力壩抗震分析。

關鍵詞：高重力壩;地震響應分析;黏彈性人工邊界;地震動輸入方法;ABAQUS二次開發(fā)

中圖分類號：TV642.4;TV312

文獻標志碼：A

doi：10.3969/j.issn.1000-1379.2021.12.023

引用格式：徐威，趙臻真，王岳，等.基于黏彈性人工邊界的高重力壩動力特性分析[J].人民黃河，2021，43（12）：118-122.

Abstract： The seismic response analysis of high gravity dam is an important step to ensure the reasonable design of buildings and evaluate their seismic characteristics. Based on the wave theory， the input method of viscoelastic artificial boundary and corresponding ground motion was obtained in the ABAQUS software， which was redeveloped based on Python language. Then the rationality of the method was proved by a simple example and applied to a gravity dam project. The results show that the method is simple， effective， easy to implement and it can ensure certain accuracy， which is suitable for the seismic analysis of a high concrete gravity dam.

Key words： high gravity dam; seismic response analysis; viscoelastic boundary condition; method of earthquake input; secondary development of ABAQUS

1 研究背景

隨著水電事業(yè)的蓬勃發(fā)展及環(huán)保理念的日益深化，我國政府把水能資源作為能源戰(zhàn)略和能源安全的積極發(fā)展領域，強調在貫徹全面協(xié)調、統(tǒng)籌兼顧、保護生態(tài)、發(fā)揮綜合效益原則的基礎上，實現人與自然和諧相處，促進經濟社會可持續(xù)發(fā)展。大壩作為水利水電發(fā)展最重要的標志[1]，其在地震作用下的安全穩(wěn)定至關重要。重力壩因其良好的受力特性及超標洪水漫頂的可靠性而被廣泛采用。近年來，我國水利水電技術發(fā)展較快，高重力壩目前在我國建設多座，其中絕大部分屹立于我國西南部等高地震烈度區(qū)，由此帶來的高重力壩地震安全問題必須加以重視[2]。

在高重力壩抗震分析中，地基的模擬與地震動輸入方法是其關鍵環(huán)節(jié)[3]。有質量地基一致性輸入模型依賴于地基截斷范圍，若不考慮人工邊界，則人為放大了地震反應。廣泛流行的無質量地基一致性輸入模型雖相較于有質量地基削弱了地震動，但因完全不考慮地基質量，故人為放大了系統(tǒng)的自振頻率，且無法考慮地震波的無限逸散問題[4]。

如何真實正確地反映地震動能量向無限地基逸散是眾多學者研究的重要內容。國內外學者考慮用一種輻射阻尼的作用，來實現在有限范圍地基條件下地震動在人工截斷邊界處不發(fā)生反射作用。黏彈性人工邊界正是其中的一種[5]，其不僅能夠吸收散射波能量，而且能夠考慮半無限地基的恢復能力。許多學者對該人工邊界進行了深入研究，劉晶波等[6]給出了二維時域黏彈性邊界，何建濤等[7]對各種實現該邊界的方法進行了系統(tǒng)總結。目前黏彈性人工邊界在地震分析中運用較為廣泛[8-13]。

對于工程設計人員，掌握在通用有限元軟件中實現黏彈性人工邊界的施加方法，是對高重力壩進行地震動力分析的關鍵。本文采用黏彈性人工邊界及其相應地震動的輸入方法對高重力壩動力特性進行了分析。

2 黏彈性人工邊界及其地震動輸入方法

2.1 黏彈性人工邊界

黏彈性人工邊界條件（Visco-Elastic Artifical Boundary Conditions）可等效為模型邊界上連續(xù)分布的并聯(lián)彈簧-阻尼器元件系統(tǒng)。彈簧元件的存在使得該邊界可以考慮介質的彈性恢復作用，這是其優(yōu)于僅含單向黏滯阻尼器、只能考慮介質能量吸收作用的黏性邊界的最重要原因。

將邊界點上的彈簧剛度、阻尼系數和各點作用的等效地震荷載組裝到與結構對應的剛度矩陣、阻尼矩陣及荷載向量組中，即可進行求解。

2.3 通用軟件中的實現方法

在通用有限元軟件中實現黏彈性人工邊界的輸入，常規(guī)的做法是在人工邊界的節(jié)點上設置法向與切向的并聯(lián)彈簧和阻尼器，借助軟件中的彈簧單元和阻尼單元實現，如ANSYS中的combine14單元，ABAQUS中的spring1和dashpot1單元及ADINA中的spring單元。由于要求對邊界上每個節(jié)點賦予彈簧剛度、阻尼系數以及相應的等效荷載，節(jié)點數量眾多，其手動施加過程較為繁瑣，不利于工程設計時采用，因此可結合Python、MATLAB等編程手段進行實現。下面以ABAQUS為例，利用Python編程平臺，介紹實現方法。

由式（5）可知，邊界上節(jié)點B的控制面積是關鍵輸入元素。本文根據潘堅文[16]提出的動水附加質量實現方法，在作用邊界上施加1 Pa的均布荷載，求得節(jié)點反力，其實質即為節(jié)點的控制面積。之后再讀取各節(jié)點的全局坐標，計算出波源至節(jié)點的距離R，再由式（1）、式（2）計算各點的人工邊界系數，由式（3）、式（4）計算各點的等效應力，最后通過.inp文件寫入相應關鍵字，實現批量施加人工邊界彈簧、阻尼器及等效應力的功能。以上施加方法通過編制程序實現。

2.4 典型算例驗證

通過一外源波動算例[17]驗證本文介紹的黏彈性人工邊界施加程序的正確性。選取二維均勻彈性半空間模型，荷載從二維彈性半空間計算模型的底部垂直輸入，有限域模型的尺寸為長400 m、高600 m，模型材料參數見表2。模型采用四節(jié)點平面應變單元（CPE4）進行離散。

該問題的自由地表位移解析解為考慮行波延遲后放大2倍的入射位移時程。由圖1可見，入射波到達地表后與反射波疊加，地表的響應為入射波的2倍;底面觀測點在經歷入射波和經地表反射的波后響應幾乎為0，波傳向遠域，與解析解一致。數值模擬結果與波傳播的理論解一致，說明本文介紹的黏彈性人工邊界及其地震動輸入方式和相關程序是正確且有效的。

3 工程實例

黃登水電站是瀾滄江上游“一庫八級”水電梯級開發(fā)方案的重要組成部分。該工程為大（1）型，工程等別為一等，主要建筑物級別為1級。

選取典型擋水壩段12#壩段進行數值模擬計算，該壩段壩頂高程1 625 m，壩基面高程1 422 m，最大壩高203 m，壩頂寬16 m，正常蓄水位1 619 m。大壩采用四節(jié)點平面應力單元（CPS4）、壩基采用四節(jié)點平面應變單元（CPE4）進行離散。

基礎采用有質量地基并在地基截斷邊界處施加黏彈性人工邊界以模擬無限域輻射阻尼效應，施加過程與方法采用本文介紹的黏彈性人工邊界施加程序。輸入地震動按實測Konya地震波輸入，其時程曲線如圖2所示。計算時水庫水位為正常蓄水位，動水壓力采用westergaard附加質量模型;Rayleigh阻尼系數根據自振分析得到的前兩階頻率計算獲得。

壩體與壩基物理力學參數見表3，不考慮大壩壩體的材料分區(qū)，地基材料假定為均質線彈性介質，忽略材料阻尼。

為分析人工邊界條件下截斷地基的尺寸效應影響，模擬范圍分別選取往上下游及深度方向1.0倍、1.5倍、2.0倍壩高地基;同時采用1.5倍壩高無質量地基設置對照，以考察輻射阻尼效應對大壩地震響應的影響。

圖3給出了無質量地基模型和黏彈性人工邊界地基模型的壩頂水平向相對位移時程曲線?？芍?，考慮輻射阻尼效應后壩頂最大相對位移顯著降低，說明與無質量地基相比，黏彈性人工邊界地基可顯著減小壩頂地震位移響應。

圖4給出了不同地基尺寸黏彈性人工邊界模型的壩頂水平向位移時程?？芍?，隨地基尺寸的增大，壩頂位移逐漸減小，但相對變化較小，數量級僅為0.001 m，壩頂位移對尺寸邊界不敏感。

圖5、圖6為無質量地基和黏彈性人工邊界地基的最大、最小主應力包絡圖。可以看出，兩種地基模型所得的主應力分布規(guī)律基本一致，壩踵處出現最大拉應力，壩趾處出現最大壓應力?？紤]地基輻射阻尼后，壩體的最大主應力及最小主應力均有較大幅度的降低。在現行抗震設計規(guī)范中，按無質量地基考慮，與實際情況相比，設計結果偏于保守，輻射阻尼會大幅降低結構的地震響應，在實際工程抗震設計中應予以適當考慮。

此外，由應力包絡圖可看出，壩體上、下游折坡處及壩踵部位是抗震薄弱環(huán)節(jié)，這些部位主拉應力較大，可能發(fā)生損傷開裂，因此進一步對大壩進行非線性動力分析。

采用ABAQUS中的混凝土塑性損傷模型，通過引入損傷因子反映混凝土材料的拉壓損傷狀態(tài)。混凝土受壓采用硬化模型，初始受壓屈服強度為13.0 MPa，極限抗壓強度為24.1 MPa，極限抗拉強度為2.9 MPa。材料模型通過開裂位移定義混凝土的受拉軟化行為。圖7給出了壩體損傷開裂在地震過程中的演化過程?？芍?，大壩首先在壩踵處發(fā)生損傷，之后在上游直立段、上下游折坡點等處發(fā)生損傷開裂，隨著地震的進行，損傷區(qū)逐漸向壩體內發(fā)展，最終在下游折坡點處發(fā)生交會貫穿整個壩面。在重力壩設計中，特別是高地震烈度區(qū)的高重力壩設計，應重點關注這些部位的抗震性能，并采取必要的工程措施。

4 結 論

（1）基于Python編程平臺，可實現黏彈性人工邊界在大型通用有限元分析軟件中的自動化施加，且典型算例的計算結果證明本文所闡述的黏彈性邊界有較好的準確性與可靠性。

（2）通過無質量地基模型和有質量黏彈性人工邊界地基模型的對比分析，說明采用人工邊界模型在考慮地基的輻射阻尼效應后，壩體的動力響應顯著降低，在一定程度更符合真實地基情況，而采用無質量地基模型進行壩體的抗震分析偏于保守。

（3）本文介紹并編制的黏彈性人工邊界及地震荷載施加程序簡單易行，且容易推廣到三維情況，適合于設計人員在實際工程設計中使用。

參考文獻：

[1] JIA Jinsheng. A Technical Review of Hydro-Project Develop[J]. Engineering，2016（3）：302-312.

[2] 邱流潮，金峰.地震分析中人工邊界處理與地震動輸入方法研究[J].巖土力學，2006，27（9）：1501-1504.

[3] 馬笙杰，遲明杰，陳紅娟，等.黏彈性人工邊界在ABAQUS中的實現及地震動輸入方法的比較研究[J].巖石力學與工程學報，2020，39（7）：1445-1457.

[4] 陳震，徐遠杰.基于波動理論的黏彈性人工邊界內源波動有限元分析[J].武漢大學學報（工學版），2011，44（6）：735-739.

[5] DEEKS A J，RANDOLPH M F. Axisymmetric Time-Domain Transmitting Boundaries[J].Journal of Engineering Mechanics，1994，120（1）：25-42.

[6] 劉晶波，呂彥東.結構—地基動力相互作用問題分析的一種直接方法[J].土木工程學報，1998，31（3）：55-64.

[7] 何建濤，馬懷發(fā)，張伯艷，等.黏彈性人工邊界地震動輸入方法及實現[J].水利學報，2010，41（8）：960-969.

[8] 王進廷，潘堅文，張楚漢.地基輻射阻尼對高拱壩非線性地震反應的影響[J].水利學報，2009，40（4）：413-420.

[9] 陳平，何蘊龍，陳海霞.地基邊界對五嘎沖拱壩地震動力響應的影響[J].武漢大學學報（工學版），2011，44（5）：554-559.

[10] 王銘明，陳健云，徐強.重力壩-庫水-地基相互作用分析方法比較研究[J].大連理工大學學報，2013，53（5）：715-722.

[11] 李渤，辛全才.考慮庫水可壓縮性的重力壩地震動力響應分析[J].人民黃河，2018，40（8）：116-119.

[12] 李笑然，胡曉，邸慶霜，等.ANSYS黏彈性邊界在拱壩抗震分析中的應用[J].人民黃河，2016，38（3）：84-88.

[13] 梁鐘元，王浩，許新勇.基于黏彈性邊界的渡槽結構地震動力特性分析[J].人民黃河，2020，42（2）：81-84，102.

[14] 吳藝，房營光.彈塑性地基中黏性與黏彈性人工邊界條件有效性的驗證[J].巖石力學與工程學報，2006，25（增刊2）：3468-3473.

[15] 劉晶波，谷音，杜義欣.一致粘彈性人工邊界及粘彈性邊界單元[J].巖土工程學報，2006，28（9）：1070-1075.

[16] 潘堅文.ABAQUS水利工程應用實例教程[M].北京：中國建筑工業(yè)出版社，2015：175-180.

[17] 苑舉衛(wèi)，杜成斌，陳燈紅.基于ABAQUS的三維黏彈性邊界單元及地震動輸入方法研究[J].三峽大學學報（自然科學版），2010，32（3）：9-13.

【責任編輯 張華巖】