王小威，陳俊濤，肖明，鄧建

?

地震作用下地下廠房混凝土襯砌結(jié)構(gòu)穩(wěn)定分析

王小威1, 2，陳俊濤1, 2，肖明1, 2，鄧建1, 2

(1. 武漢大學(xué) 水資源與水電工程科學(xué)國家重點(diǎn)實驗室，湖北 武漢，430072；2. 武漢大學(xué) 水工巖石力學(xué)教育部重點(diǎn)實驗室，湖北 武漢，430072)

針對混凝土循環(huán)荷載作用下的率相關(guān)性和疲勞損傷特性，提出考慮拉、剪損傷各向異性的混凝土動力本構(gòu)模型，并推導(dǎo)出損傷變量動態(tài)演化方程；針對圍巖與襯砌動力相互作用特點(diǎn)，提出考慮圍巖與襯砌非連續(xù)變形的顯式動力有限元積分方法；由此，構(gòu)建地下廠房混凝土襯砌結(jié)構(gòu)動力響應(yīng)分析方法，并將該方法應(yīng)用于映秀灣水電站地下廠房抗震穩(wěn)定計算。襯砌結(jié)構(gòu)各部位在地震下處于同步震動狀態(tài)，其中襯砌頂拱應(yīng)力和位移波動幅度較大；主廠房襯砌頂拱與邊墻最大相對位移達(dá)到1.6 cm，說明主廠房發(fā)生了較為明顯的結(jié)構(gòu)變形；襯砌結(jié)構(gòu)損傷區(qū)主要分布在頂拱處，其中拉裂損傷造成的損傷破壞較為嚴(yán)重；圍巖對襯砌結(jié)構(gòu)的約束作用有助于襯砌結(jié)構(gòu)的抗震穩(wěn)定性。應(yīng)用上述方法所得的計算結(jié)果較為合理地反映了襯砌結(jié)構(gòu)動力響應(yīng)特征，可為地下廠房混凝土襯砌結(jié)構(gòu)抗震設(shè)計提供參考。

地下廠房；混凝土襯砌；動力響應(yīng)；拉、剪損傷各向異性；動力本構(gòu)；非連續(xù)變形；顯式動力有限元

中國西南高地震烈區(qū)分布著多座大型水電站，這些水電站地下廠房在地震作用下的穩(wěn)定性對水電站安全運(yùn)行至關(guān)重要。人們對地下廠房動力響應(yīng)進(jìn)行了大量研究[1?4]，取得了豐富的研究成果，但是對廠房襯砌結(jié)構(gòu)動力響應(yīng)特征研究較少。張雨霆等[5?6]的研究表明，地下廠房襯砌結(jié)構(gòu)是地下廠房抗震設(shè)計的重要組成部分，因此研究其地震動力響應(yīng)具有重要意義。地下廠房襯砌結(jié)構(gòu)地震動力響應(yīng)分析主要包括2個方面的內(nèi)容：一是襯砌混凝土材料的動態(tài)響應(yīng)分析，二是地震作用下洞室圍巖與襯砌相互作用分析?；炷敛牧系膭恿Ρ緲?gòu)模型是地下廠房襯砌結(jié)構(gòu)動力分析的基礎(chǔ)，許多研究者基于室內(nèi)實驗或理論推導(dǎo)建立了各種混凝土動態(tài)本構(gòu)模型以研究混凝土材料的動力響應(yīng)規(guī)律，主要包括彈塑性動力本構(gòu)模型[7?9]、非線性彈性動力本構(gòu)模型[10?11]、黏塑性動力本構(gòu)模型[12?13]和動力損傷本構(gòu)模型[14?16]，其中動力損傷本構(gòu)模型應(yīng)用最為廣泛。由于自身結(jié)構(gòu)特征，混凝土在加載前其內(nèi)部就分布著微小缺陷，在地震循環(huán)荷載作用下，這些缺陷不斷發(fā)展、貫通。由此可見，混凝土的破壞過程是損傷累積發(fā)展的過程，要全面了解混凝土的動力響應(yīng)特征就需要研究其內(nèi)部缺陷的演化過程[17]，這也使得損傷理論被廣泛應(yīng)用于混凝土的動態(tài)響應(yīng)分析[18?19]。關(guān)于圍巖與襯砌相互作用的研究，人們展開了大量工作。朱維申[20]建立了黏彈塑性巖體中圍巖與襯砌共同作用時一個軸對稱問題的解析表達(dá)式；肖明等[21]采用彈塑性有限元研究了洞室開挖變形對襯砌作用的影響；CHEN等[22]采用三維非線性薄層單元研究了圍巖與襯砌相互作用?？梢钥闯觯@些研究主要針對靜力情況，對地震荷載作用下圍巖與襯砌的動響應(yīng)問題關(guān)注較少。事實上，圍巖與襯砌作為一個共同作用體系，其在地震荷載作用下的響應(yīng)與靜力條件下的有很大不同；靜力條件下可認(rèn)為圍巖與襯砌變形協(xié)調(diào)，兩者在接觸面處位移是連續(xù)的，但是在動力條件下，由于兩者材料的差異性，加上接觸面對地震波傳播的影響，兩者在接觸面處的位移是不連續(xù)的。因此，在構(gòu)建圍巖與襯砌動力條件下的相互作用模型時，需要考慮到圍巖與襯砌在地震荷載作用下的非連續(xù)變形?；谏鲜鏊悸?，本文作者建立混凝土拉、剪損傷各向異性動力本構(gòu)模型，推導(dǎo)出損傷變量動態(tài)演化方程，提出考慮圍巖與襯砌非連續(xù)變形的顯式動力有限元積分方法，并在此基礎(chǔ)上構(gòu)建地下廠房混凝土襯砌結(jié)構(gòu)動力響應(yīng)分析方法。將該方法應(yīng)用到映秀灣水電站地下廠房動力計算中，研究地下廠房混凝土襯砌結(jié)構(gòu)震損演化過程，為地下廠房襯砌結(jié)構(gòu)抗震設(shè)計提供參考。

1 混凝土襯砌動力損傷本構(gòu)模型

前期研究表明[23?26]，在動力循環(huán)荷載作用下，混凝土材料表現(xiàn)出明顯的率相關(guān)性和疲勞損傷特性，其率相關(guān)性主要體現(xiàn)在物理力學(xué)參數(shù)隨應(yīng)變率的提高而增加?；炷潦且环N非均質(zhì)的多相材料，在地震循環(huán)荷載作用下，其內(nèi)部固有的微裂隙和微空洞不斷地擴(kuò)展和累積，最終形成宏觀上的損傷破壞。在這個過程中，混凝土材料的強(qiáng)度和剛度不斷退化，表現(xiàn)出疲勞損傷特性?；谝陨鲜聦嵖芍獦?gòu)建混凝土襯砌動力本構(gòu)模型，則需要包含這2個方面特征。因此，本文作者在損傷理論的基礎(chǔ)上，建立考慮應(yīng)變率效應(yīng)的混凝土動力損傷本構(gòu)模型。

1.1 基于拉、剪分離的混凝土動力損傷模型

混凝土動力損傷本構(gòu)需要反映混凝土的率敏感性和損傷的產(chǎn)生與演化，其一般形式可以寫為

對損傷變量描述的合理性很大程度上決定了損傷模型的準(zhǔn)確性，目前工程計算采用的模型中，損傷變量形式主要為標(biāo)量形式，其形式簡單，便于理論推導(dǎo)和大規(guī)模數(shù)值計算。然而混凝土作為一種非均質(zhì)材料，在荷載作用下，特別是在地震循環(huán)荷載作用下，其損傷演化是具有方向性的，因此損傷的度量應(yīng)該采用張量的形式[17]?；炷恋膿p傷形態(tài)主要有受拉損傷和受剪損傷，2種損傷形態(tài)的產(chǎn)生機(jī)制各不相同，損傷對強(qiáng)度和剛度的影響也表現(xiàn)迥異。因此將動態(tài)損傷變量分為受拉和受剪2個部分[27]。

1.2 混凝土拉、剪損傷變量演化方程

由前面推導(dǎo)過程可知，要構(gòu)建完整的混凝土動態(tài)損傷模型，還需要建立動力作用下混凝土拉、剪損傷變量演化方程。動力作用下的拉、剪損傷演化過程與靜力作用下不同之處在于：動力損傷的演化過程不僅與應(yīng)變有關(guān)，而且與應(yīng)變率相關(guān)。因此，建立率相關(guān)的損傷演化方程是關(guān)鍵。

假設(shè)在三維復(fù)雜應(yīng)力狀態(tài)下，單元的受拉破壞滿足最大拉應(yīng)變法則。按照文獻(xiàn)[28]中的方法，采用等效應(yīng)變來表示三維應(yīng)力狀態(tài)下的最大拉應(yīng)變，則受拉損傷變量演化方程可以表示為

式(8)~(12)構(gòu)成了完整的動力受拉損傷變量演化方程。類似地，三維應(yīng)力狀態(tài)下的剪切損傷變量演化方程可以表示為

同理可知，式(13)~(16)構(gòu)成了完整的動力剪切損傷變量演化方程。由以上分析可知，本文的動態(tài)受拉損傷變量與剪切損傷變量是直接在最大拉應(yīng)變準(zhǔn)則和Mohr-Coulomb準(zhǔn)則的基礎(chǔ)上，考慮應(yīng)變率效應(yīng)而得出的，形式簡單，易于有限元計算。

2 圍巖與襯砌非連續(xù)變形的動力分析方法

地下廠房可以看作是混凝土結(jié)構(gòu)與圍巖的聯(lián)合體?；炷烈r砌與圍巖在地震波作用下發(fā)生受迫振動，由于混凝土介質(zhì)與巖石介質(zhì)的差異性，加上圍巖與襯砌間存在非連續(xù)接觸面，圍巖與襯砌對于地震動的響應(yīng)存在明顯差異。這種差異造成了襯砌與圍巖非連續(xù)變形，即襯砌與圍巖間存在相對位移差，由相對位移差引起的應(yīng)力差則是圍巖與襯砌間發(fā)生損傷破壞的主要原因。本文作者基于顯式積分法，結(jié)合圍巖與襯砌接觸條件，建立圍巖與襯砌非連續(xù)變形動力分析方法。

2.1 圍巖與襯砌介質(zhì)運(yùn)動方程的時域顯式積分方法

根據(jù)襯砌與圍巖的有條件聯(lián)合承載特征，可以將襯砌與圍巖地震作用下的聯(lián)合承載過程看作是混凝土襯砌與巖體這2種材料有條件單獨(dú)承載過程的疊加。依據(jù)圍巖與襯砌的動力有限元離散方程，并結(jié)合襯砌與圍巖接觸面的初始條件，可以推導(dǎo)出兩介質(zhì)的時域有限元逐步積分格式。

經(jīng)過有限元離散后，地下廠房結(jié)構(gòu)模型可以看作是多質(zhì)點(diǎn)組成的多自由度體系，地震作用下圍巖的動力平衡方程可以表示為

根據(jù)時域加權(quán)殘數(shù)的顯式逐步積分格式求解巖體運(yùn)動方程，可以得到巖體介質(zhì)運(yùn)動方程的時域顯式積分表達(dá)式為

依此類推，襯砌混凝土介質(zhì)動力平衡方程為

同樣，襯砌混凝土介質(zhì)運(yùn)動方程的時域顯式積分表達(dá)式為

假設(shè)圍巖介質(zhì)和襯砌混凝土介質(zhì)靜力開挖后的狀態(tài)為初始條件，則有

2.2 圍巖與襯砌動力相互作用分析

圍巖與襯砌間相互作用關(guān)系可以看作是兩介質(zhì)在接觸面處變形不連續(xù)，而應(yīng)力連續(xù)。兩介質(zhì)在接觸面對應(yīng)節(jié)點(diǎn)處的相對位移可以表示為

由接觸面對應(yīng)節(jié)點(diǎn)處的相對位移引起的節(jié)點(diǎn)力增量可以表示為

由+1時刻接觸面節(jié)點(diǎn)力增量可以得到此時的節(jié)點(diǎn)力為

由于式(36)和式(37)中法向和切向接觸力是在圍巖與襯砌彈性接觸狀態(tài)下求得，而圍巖與襯砌接觸面是地下廠房結(jié)構(gòu)中薄弱部位，其在地震循環(huán)荷載作用下的損傷是不可忽視的，因此需要對接觸面損傷后的接觸力進(jìn)行修正。

接觸面的破壞形式主要是沿法向的張拉開裂和沿切向的剪切滑移。本文作者從這2個方面對接觸狀態(tài)進(jìn)行校核并修正接觸力。

圖1 圍巖與襯砌非連續(xù)變形的動力分析流程圖

3 工程實例

3.1 工程介紹

映秀灣水電站位于四川省汶川縣映秀鎮(zhèn)的岷江干流上，其地下廠房埋深約200 m，主廠房長×寬×高為52.8 m×17.0 m×37.2 m，主廠房邊墻襯砌厚0.5 m，頂拱厚0.9 m左右，主變洞長×寬×高為59.4 m× 7.2 m×27.9 m。水電站工程場地距汶川地震震中約 8 km，地震影響烈度高達(dá)Ⅸ度。由于其距離震中較近，相當(dāng)于經(jīng)受了一次地震原型實驗，因此其震后地下廠房襯砌震損情況具有很強(qiáng)的代表性。本文作者以映秀灣水電站地下廠房為工程實例，采用課題組自主研發(fā)的動力有限元程序[35]，結(jié)合本文的混凝土襯砌動力響應(yīng)分析模塊，研究地下廠房混凝土襯砌地震作用下的響應(yīng)特征。

3.2 有限元模型與計算條件

選取映秀灣水電站主廠房和主變洞來建立數(shù)值模型。模型采用8節(jié)點(diǎn)6面體單元進(jìn)行離散，共劃分為97 652個單元，其中混凝土襯砌單元7 298個，整體三維有限元模型剖面如圖2所示。取=0平面為典型斷面，取該斷面主廠房和主變洞襯砌上共6個監(jiān)測點(diǎn)為研究對象。襯砌有限元模型和監(jiān)測點(diǎn)分布如圖3 所示。

在襯砌結(jié)構(gòu)進(jìn)行動力分析之前，需要對地下廠房進(jìn)行開挖和支護(hù)模擬，將靜力計算的結(jié)果作為動力響應(yīng)分析計算的初始條件。計算采用的材料物理力學(xué)參數(shù)如表1所示。

圖2 有限元模型剖面

(a) 主廠房襯砌；(b) 主變洞襯砌

表1 模型材料力學(xué)參數(shù)表

模型四周采用自由場人工邊界，底部采用黏彈性人工邊界，頂部模擬至自由面。地震波采用距汶川地震較近的臥龍臺測得的加速度時程，截取其中變化劇烈、幅值較大的20~40 s時段，經(jīng)過濾波、基線校正和幅值折減、坐標(biāo)轉(zhuǎn)換[5]后，作為地震輸入波從模型底部輸入。輸入地震波在3個方向的加速度時程曲線如圖4所示。

(a) X方向；(b) Y方向；(c) Z方向

3.3 地下廠房襯砌結(jié)構(gòu)動力響應(yīng)分析

3.3.1 襯砌監(jiān)測點(diǎn)應(yīng)力與位移時程分析

圖5所示為襯砌結(jié)構(gòu)6個監(jiān)測點(diǎn)的位移時程曲線。從圖5可以看出：6條時程曲線的波形和波動規(guī)律較為相似，均出現(xiàn)了3個明顯的波峰，波峰和波谷幾乎同時出現(xiàn)，說明襯砌結(jié)構(gòu)各監(jiān)測點(diǎn)處于同步震動狀態(tài)。在0~6 s時間段內(nèi)，各監(jiān)測點(diǎn)位移變化量很小，在6 s之后，監(jiān)測點(diǎn)位移曲線開始大幅度波動，波動范圍為1.5~5.5 cm，3個波峰分別出現(xiàn)在10，13.2和16.5 s附近，其對應(yīng)的峰值分別為4.8，5.4和5.0 cm。地震波輸入完成后，監(jiān)測點(diǎn)位移與震前相比增加了0.5~ 1.0 cm。各個監(jiān)測點(diǎn)的位移時程曲線間的差異主要體現(xiàn)在波動幅值上，主廠房邊墻和處位移波動幅值較小，與主廠房頂拱處的最大位移差達(dá)到1.6 cm，表明主廠房襯砌結(jié)構(gòu)有較大相對變形；主變洞3個監(jiān)測點(diǎn)位移幅值均較大，但相互之間的差值較小，說明主變洞襯砌結(jié)構(gòu)主要是整體變形。

混凝土材料抗拉強(qiáng)度遠(yuǎn)小于抗壓強(qiáng)度，因此拉裂破壞是混凝土材料常見的破壞形式。取，和這3個監(jiān)測點(diǎn)的最大主應(yīng)力(拉應(yīng)力為正)時程進(jìn)行分析，時程曲線如圖6所示。開挖完成后，3個監(jiān)測點(diǎn)的拉應(yīng)力范圍為0.8~0.9 MPa。隨著地震荷載的施加，在0~4 s，拉應(yīng)力波動幅度較小，主要在0.8~1.0 MPa范圍內(nèi)上下波動；在4~15 s，拉應(yīng)力劇烈波動，波動范圍為0.45~1.64 MPa，其中點(diǎn)和點(diǎn)的拉應(yīng)力最大值分別為1.64 MPa和1.52 MPa，均超過混凝土襯砌抗拉強(qiáng)度；在15~20 s，拉應(yīng)力波動幅度有所減小，最終拉應(yīng)力為0.8~1.0 MPa，與靜力開挖后的應(yīng)力相比有所增加。從圖6可以看出：監(jiān)測點(diǎn)最大主應(yīng)力時程曲線與輸入地震波波形較為相似，拉應(yīng)力的波動規(guī)律主要受到輸入地震波的影響。點(diǎn)和點(diǎn)波動幅度比點(diǎn)的大，這表明地震波荷載對襯砌頂拱影響較大，加劇了頂拱拉應(yīng)力，使得頂拱損傷破壞嚴(yán)重，這也與文獻(xiàn)[32]中的調(diào)查結(jié)果相吻合。

1—A監(jiān)測點(diǎn)；2—B監(jiān)測點(diǎn)；3—C監(jiān)測點(diǎn)；4—D監(jiān)測點(diǎn)；5—E監(jiān)測點(diǎn)；6—F監(jiān)測點(diǎn)。

1—A監(jiān)測點(diǎn)；2—C監(jiān)測點(diǎn)；3—D監(jiān)測點(diǎn)。

4.3.2 襯砌結(jié)構(gòu)損傷情況分析

為了直觀地顯示襯砌結(jié)構(gòu)損傷程度，將損傷張量d轉(zhuǎn)化為標(biāo)量損傷系數(shù)b，b為損傷張量主值的平均值。圖7所示為襯砌結(jié)構(gòu)震前和震后損傷系數(shù)分布圖。由圖7可知：震前(=0 s)損傷區(qū)分布在主廠房襯砌頂拱與邊墻交界處、主變洞襯砌頂拱和頂拱與邊墻交界處，損傷區(qū)范圍較小，損傷系數(shù)主要為0.1~0.2，最大損傷系數(shù)小于0.3。在地震波輸入完成后(=20 s)，襯砌結(jié)構(gòu)損傷區(qū)范圍和損傷系數(shù)均有較大程度增加，損傷區(qū)主要向襯砌頂拱和邊墻擴(kuò)展，主廠房襯砌結(jié)構(gòu)頂拱大部分區(qū)域及其與邊墻交界處均有不同程度損傷，主變洞襯砌結(jié)構(gòu)頂拱及其與邊墻交界處損傷區(qū)進(jìn)一步擴(kuò)大，邊墻部分區(qū)域也出現(xiàn)了損傷情況；主廠房頂拱部分區(qū)域和頂拱與邊墻交界處損傷嚴(yán)重，損傷系數(shù)最大值接近于1，主變洞襯砌結(jié)構(gòu)損傷系數(shù)為0~ 0.8，最大損傷系數(shù)位于襯砌端部頂拱與邊墻交界處。

(a) t=0 s; (b) t=20 s

震后襯砌結(jié)構(gòu)拉、剪損傷分布情況如圖8所示。損傷區(qū)主要分布在襯砌頂拱及其與邊墻交界處，并且有向邊墻擴(kuò)展的趨勢，其中拉裂損傷分布在主廠房頂拱頂部和頂拱與邊墻交界處、主變室頂拱與邊墻交界處，剪切損傷主要分布在主廠房和主變室頂拱與邊墻處。結(jié)合圖7~8可知：拉裂損傷區(qū)域損傷系數(shù)較大，襯砌損傷破壞嚴(yán)重，說明在進(jìn)行襯砌結(jié)構(gòu)抗震設(shè)計時，需要采取措施限制襯砌結(jié)構(gòu)拉裂破壞。

圖8 襯砌結(jié)構(gòu)拉、剪損傷分布

4.3.3 襯砌與圍巖接觸面滑移開裂情況分析

襯砌結(jié)構(gòu)是附著在圍巖上的，如果將圍巖看作是一個整體，那么襯砌與圍巖的整體性明顯要弱于圍巖自身的整體性。因此，在地震循環(huán)荷載作用下，圍巖與襯砌變形的一致性較差，圍巖與襯砌接觸面出現(xiàn)滑移、開裂情況。襯砌與圍巖接觸面滑移與開裂區(qū)分布如圖9所示。由接觸面滑移、開裂分布情況可知，滑移開裂區(qū)主要分布在襯砌頂拱、邊墻、頂拱與邊墻交界處。結(jié)合圖8和圖9可知：接觸面滑移、開裂區(qū)與襯砌結(jié)構(gòu)損傷區(qū)分布情況較為一致，說明接觸面滑移開裂區(qū)襯砌結(jié)構(gòu)更容易發(fā)生損傷情況，圍巖對襯砌結(jié)構(gòu)的約束作用可以有效減小地震荷載對襯砌結(jié)構(gòu)的損傷破壞。

圖9 襯砌與圍巖接觸面滑移與開裂區(qū)分布

4 結(jié)論

1) 基于混凝土拉、剪損傷各向異性動力本構(gòu)模型和圍巖與襯砌動力接觸分析方法，建立地下廠房混凝土襯砌結(jié)構(gòu)動力響應(yīng)分析模型，并將該模型應(yīng)用于映秀灣地下廠房襯砌結(jié)構(gòu)震損演化過程分析。

2) 各監(jiān)測點(diǎn)位移時程曲線和應(yīng)力時程曲線較為相似，襯砌結(jié)構(gòu)各監(jiān)測點(diǎn)在地震荷載作用下同步震動。主廠房襯砌頂拱與邊墻相比，位移和拉應(yīng)力波動幅度均較大，頂拱拉應(yīng)力超過了襯砌抗拉強(qiáng)度，其位移與邊墻位移最大相差1.6 cm，表明主廠房襯砌結(jié)構(gòu)有較大相對變形；主變室襯砌各點(diǎn)位移較為接近，說明主變洞襯砌結(jié)構(gòu)主要是整體變形。

3) 襯砌結(jié)構(gòu)損傷區(qū)主要分布在襯砌頂拱及其與邊墻交界處，并且有向邊墻擴(kuò)展的趨勢，這與映秀灣地下廠房震害調(diào)查結(jié)果相吻合。主廠房頂拱部分區(qū)域和頂拱與邊墻交界處損傷最為嚴(yán)重，拉裂損傷區(qū)域損傷系數(shù)較大，襯砌損傷破壞嚴(yán)重，需要采取措施限制襯砌結(jié)構(gòu)拉裂破壞。

4) 襯砌與圍巖接觸面滑移開裂區(qū)與襯砌結(jié)構(gòu)損傷區(qū)分布情況較為一致，說明接觸面滑移開裂區(qū)的襯砌結(jié)構(gòu)更容易損傷破壞，圍巖對襯砌結(jié)構(gòu)的約束作用有助于襯砌結(jié)構(gòu)的抗震穩(wěn)定性，這也印證了地下工程與地面工程相比有較強(qiáng)抗震性。

[1] 李海波, 朱蒞, 呂濤, 等. 考慮地震動空間非一致性的巖體地下洞室群地震反應(yīng)分析[J]. 巖石力學(xué)與工程學(xué)報, 2008, 27(9): 1758?1766. LI Haibo, ZHU Li, Lü Tao. Seismic response analysis of an underground cavern groups in rock subjected to spatially non-uniform seismic ground motion[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(9): 1758?1766.

[2] 李小軍, 盧淘. 水電站地下廠房洞室群地震反應(yīng)顯示有限元分析[J]. 水利發(fā)電學(xué)報, 2009, 28(5): 1877?1882. LI Xiaojun,LU Tao. Explicit finite element analysis of earthquake response for underground caverns of hydropower stations[J]. Journal of Hydroelectric Engineering, 2009, 28(5): 1877?1882.

[3] 李海波, 馬行東, 李俊如, 等. 地震荷載作用下地下巖體洞室位移特征的影響因素分析[J]. 巖土工程學(xué)報, 2006, 28(3): 358?362. LI Haibo, MA Xingdong, LI Junru, et al. Study on influence factors of rock cavern displacement under earthquake[J]. Chinese Journal of Geotechnical Engineering, 2006, 28(3): 358?362.

[4] 陳健云, 胡志強(qiáng), 林皋. 超大型地下洞室群的隨機(jī)地震響應(yīng)分析[J]. 水利學(xué)報, 2002(1): 71?75.CHEN Jianyun, HU Zhiqiang, LIN Gao. Random seismic response analysis of large-scale underground caverns[J]. Journal of Hydraulic Engineering, 2002(1): 71?75.

[5] 張雨霆, 肖明, 李玉婕. 汶川地震對映秀灣水電站地下廠房的震害影響及動力響應(yīng)分析[J]. 巖石力學(xué)與工程學(xué)報, 2010, 29(增2): 3663?3671. ZHANG Yuting, XIAO Ming, LI Yujie. Effect of Wenchuan earthquake on earthquake damage and dynamic response analysis of underground powerhouse of Yingxiuwan hydropower station[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(S2): 3663?3671.

[6] 崔臻, 盛謙, 冷先倫, 等. 大型地下洞室群地震動力災(zāi)變研究綜述[J]. 防災(zāi)減災(zāi)工程學(xué)報, 2013, 33(5): 606?616. CUI Zhen, SHENG Qian, LENG Xianlun, et al. A review of study on seismic catastrophe of large-scale underground cavern group[J]. Journal of Disaster Prevention & Mitigation Engineering, 2013, 33(5): 606?616.

[7] 吳建營, 李杰. 考慮應(yīng)變率效應(yīng)的混凝土動力彈塑性損傷本構(gòu)模型[J]. 同濟(jì)大學(xué)學(xué)報(自然科學(xué)版), 2006, 34(11): 1427?1430. WU Jianying, LI Jie. Elastoplastic damage constitutive model for concrete considering strain rate effect under dynamic loading[J]. Journal of Tongji University(Natural Science), 2006, 34(11): 1427?1430.

[8] 吳建營, 李杰. 反映阻尼影響的混凝土彈塑性損傷本構(gòu)模型[J]. 工程力學(xué), 2006, 23(11): 116?121. WU Jianying, LI Jie. Elastoplastic damage model for concrete considering damping effects[J]. Engineering Mechanics, 2006, 23(11): 116?121.

[9] WINNICKI A, CICHON C. Plastic model for concrete in plane stress state[J]. Journal of Engineering Mechanics, 1998, 124(6): 591?602.

[10] TEDESCO J W, POWELL J C, ROSS C A, et al. A strain-rate-dependent concrete material model for ADINA[J]. Computers and Structures, 1997, 64(5/6): 1053?1067.

[11] SHKOLNIK I E. Influence of high strain rates on stress strain relationship, strength and elastic modulus of concrete[J]. Cement & Concrete Composites, 2008, 30(10): 1000?1012.

[12] 李杰, 任曉丹, 黃橋平. 混凝土黏塑性動力損傷本構(gòu)關(guān)系[J]. 力學(xué)學(xué)報, 2011, 43(1): 193?201. LI Jie, REN Xiaodan, HUANG Qiaoping. A viscoplastic damage constitutive model for concrete[J]. Chinese Journal of Theoretical and Applied Mechanics, 2011, 43(1): 193?201.

[13] 肖詩云, 李宏男, 杜榮強(qiáng), 等. 混凝土四參數(shù)動態(tài)本構(gòu)模型[J]. 哈爾濱工業(yè)大學(xué)學(xué)報, 2006, 38(10): 1754?1757. XIAO Shiyun, LI Hongnan, DU Rongqiang, et al. Four-parameter dynamic constitutive model of concrete[J]. Journal of Harbin Institute of Technology,2006, 38(10): 1754?1757.

[14] 李慶斌, 張楚漢, 王光綸. 單軸狀態(tài)下混凝土的動力損傷本構(gòu)模型[J]. 水利學(xué)報, 1994(12): 55?60. LI Qingbin, ZHANG Chuhan, WANG Guanglun. Dynamic damage constitutive model of concrete in uniaxial tension and compression[J]. Journal of Hydraulic Engineering, 1994(12): 55?60.

[15] SHKOLNIK I E. Influence of high strain rates on stress-strain relationship, strength and elastic modulus of concrete[J]. Cement and Concrete Composites, 2008, 30(10): 1000?1012.

[16] 李慶斌, 鄧宗才, 張立翔. 考慮初始彈模變化的混凝土動力損傷本構(gòu)模型[J]. 清華大學(xué)學(xué)報(自然科學(xué)版), 2003, 43(8): 1088?1091. LI Qingbin, DENG Zongcai, ZHANG Lixiang. Dynamic damage constitutive model for concrete accounting for the difference of the initial elastic modulus[J]. Journal of Tsinghua University (Science and Technology), 2003, 43(8): 1088?1091.

[17] 林皋, 劉軍, 胡志強(qiáng). 混凝土損傷類本構(gòu)關(guān)系研究現(xiàn)狀與進(jìn)展[J]. 大連理工大學(xué)學(xué)報, 2010, 50(6): 1055?1064. LIN Gao, LIU Jun, HU Zhiqiang. Current situation and progress of research on damage constitutive relation of concrete[J]. Journal of Dalian University of Technology, 2010, 50(6): 1055?1064.

[18] 杜成斌, 蘇擎柱. 混凝土壩地震動力損傷分析[J]. 工程力學(xué), 2003, 20(10): 170?175. DU Chengbin, SU Qingzhu.Dynamic damage of concrete gravity dams under earthquake excitation[J]. Engineering Mechanics,2003, 20(10): 170?175.

[19] 陳建云, 李靜, 林皋. 基于率性損傷模型的混凝土重力壩地震響應(yīng)分析[J]. 哈爾濱工業(yè)大學(xué)學(xué)報, 2005, 37(6): 786?789. CHEN Janyun, LI Jing, LIN Gao. Seismic response analysis of concrete gravity dam based on strain rate-dependant damage model[J]. Journal of Harbin Institute of Technology,2005, 37(6): 786?789.

[20] 朱維申. 粘彈塑性介質(zhì)中圍巖與襯砌的應(yīng)力狀態(tài)[J]. 力學(xué)學(xué)報, 1981, 17(1): 56?67. ZHU Weishen. Stress state of rock and lining in viscoelastic plastic rock media[J]. Chinese Journal of Theoretical and Applied Mechanics, 1981, 17(1): 56?67.

[21] 肖明, 楊建東, 賴旭. 大型長廊阻抗式尾水調(diào)壓井結(jié)構(gòu)受力分析[J]. 巖石力學(xué)與工程學(xué)報, 2000, 19(4): 476?480. XIAO Ming, YANG Jiandong, LAI Xu. Mechanics on large long corridor of surge tank[J]. Chinese Journal of Rock Mechanics and Engineering, 2000, 19(4): 476?480.

[22] CHEN Juntao, YANG Yang, YE Chao Three-dimensional numerical analysis of compound lining in complex underground surge-shaft structure[J]. Mathematical Problems in Engineering, 2015, 2015(1):1-10.

[23] RAPHAEL J M. Tensile strength of concrete[J]. ACI Journal, 1984, 81(2): 158?165.

[24] 林皋, 陳健云, 肖詩云. 混凝土的動力特性與拱壩的非線性地震響應(yīng)[J]. 水利學(xué)報, 2003, 6: 30?36. LIN Gao, CHEN Janyun, XIAO Shiyun. Dynamic behavior of concrete and nonlinear seismic response of arch dam[J]. Shui Li Xue Bao, 2003, 6: 30?36.

[25] 肖詩云, 林皋. Drucker-Prager材料一致率型本構(gòu)模型[J]. 工程力學(xué), 2003, 20(4): 147?151. XIAO Shiyun, LIN Gao. A Drucker-Prager consistent rate-dependent model[J]. Engineering Mechanics, 2003, 20(4): 147?151.

[26] 葛修潤, 蔣宇, 盧允德, 等. 周期荷載作用下巖石疲勞變形特性試驗研究[J]. 巖石力學(xué)與工程學(xué)報, 2003, 22(10): 1581?1585. GE Xiurun, JIANG Yu, LU Yunde, et al. Testing study on fatigue deformation law of rock under cyclic loading[J]. Chinese Journal of Rock Mechanics and Engineering, 2003, 22(10): 1581?1585.

[27] 李杰, 任曉丹. 混凝土靜力與動力損傷本構(gòu)模型研究進(jìn)展述評[J]. 力學(xué)進(jìn)展, 2010, 40(3): 284?297. LI Jie, REN Xiaodan. A review on the constitutive model for static and dynamic damage of concrete[J]. Advances Mechanics, 2010, 40(3): 284?297.

[28] WANG S Y, SUN L, YANG C, et al. Numerical study on static and dynamic fracture evolution around rock cavities[J]. Journal of Rock Mechanics and Geotechnical Engineering, 2013, 5(4): 262?276.

[29] ZHOU X Q, HAO H. Modelling of compressive behaviour of concrete-like materials at high strain rate[J]. International Journal of Solids and Structures, 2008, 45(17): 4648?4661.

[30] 張志國, 肖明, 陳俊濤. 大型地下洞室地震災(zāi)變過程三維動力有限元模擬[J]. 巖石力學(xué)與工程學(xué)報, 2011, 30(3): 509?523. ZHANG Zhiguo, XIAO Ming,CHEN Juntao. Simulation of earthquake disaster process of large-scale underground caverns using three-dimensional dynamic finite element method[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(3): 510?523.

[31] 楊陽. 水電站地下廠房圍巖與結(jié)構(gòu)地震響應(yīng)分析[D]. 武漢: 武漢大學(xué)水利水電學(xué)院, 2015: 18. YANG Yang. Seismic response analysis of surrounding rock and structure for underground powerhouse of hydropower station[D]. Wuhan:Wuhan University. School of Water Resources and Hydropower Engineering, 2015: 18.

[32] 尹顯俊, 王光綸. 巖體結(jié)構(gòu)面切向循環(huán)加載本構(gòu)關(guān)系研究[J]. 巖石力學(xué)與工程學(xué)報, 2005, 24(7): 1158?1163. YIN Xianjun, Wang Guanglun. Study of consititutive model for rock interfaces and joins under normal cyclic loading. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(7): 1158?1163.

[33] GRASSELLI G, EGGER P. Constitutive law for the shear strength of rock joints based on three-dimensional surface parameters[J]. International Journal of Rock Mechanics and Mining Sciences, 2003, 40: 25?40.

[34] BANDIS S C, LUMSDEN A C, BARTON N R. Fundamentals of rock joint deformation[J]. International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts, 1983, 20(6): 249?268.

[35] 張志國, 肖明, 張雨霆, 等. 大型地下洞室三維彈塑性損傷動力有限元分析[J]. 巖石力學(xué)與工程學(xué)報, 2010, 29(5): 982?989. ZHANG Zhiguo, XIAO Ming,ZHANG Yuting, et al. Dynamic finite element analysis of large-scale complex underground caverns with three-dimensional elastoplastic damage model[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(5): 982?989.

(編輯 伍錦花)

Stability analysis of concrete lining structure in underground powerhouse with seismic action

WANG Xiaowei1, 2, CHEN Juntao1, 2, XIAO Ming1, 2, DENG Jian1, 2

(1. State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, China;2. Key Laboratory of Rock Mechanics in Hydraulic Structural Engineering, Ministry of Education, Wuhan University, Wuhan430072, China)

Based on rate-dependant and fatigue damage characteristics of concrete under cyclic loading, a dynamic constitutive model that considers tension and shear damage anisotropic was proposed, and the evolution equations of damage variables were derived. According to dynamic interaction characteristics of surrounding rock and lining, a explicit dynamic finite element analysis method that considers the discontinuous deformation between surrounding rock and lining was presented. Thus, a dynamic response analysis method of concrete lining structure in underground powerhouse was constructed and then applied to calculate seismic stability of Yingxiuwan Hydropower Station. The results show that various positions of lining structure are in a synchronous vibration state, but the stress and displacement response of the top arch are larger than that ofother positions. The maximum relative displacement of the top arch and side wall of the main powerhouse reaches 1.6 cm, indicating that the main building has a obvious structural deformation. Damage zone of lining structure is mainly distributed in its top arch, and the damage caused bytension is more serious. The confinement effects of the surrounding rock to the lining structure is helpful for seismic stability of the lining structure. Calculation results from the proposed method reflect seismic response characteristics of concrete lining reasonably, and thus provide references for seismic design of concrete lining in underground powerhouse.

underground powerhouse; concrete lining; dynamic response; tension and shear damage anisotropic; dynamic constitutive; discontinuous deformation; explicit dynamic finite element

10.11817/j.issn.1672-7207.2018.05.023

TV45

A

1672?7207(2018)05?1203?10

2017?05?08；

2017?06?28

國家重點(diǎn)基礎(chǔ)研究發(fā)展計劃(973計劃)項目(2015CB057904)；國家自然科學(xué)基金資助項目(51579191) (Project(2015CB057904) supported by the National Program on Key Basic Research Project (973 Program); Project(51579191) supported by the National Natural Science Foundation of China)

陳俊濤，博士，副教授，從事地下工程和工程軟件開發(fā)等方面的研究；E-mail: chenjt000@163.com