王樂華,鄧華鋒,湯開宇,朱 敏,彭云楓

(1.三峽大學(xué)三峽庫區(qū)地質(zhì)災(zāi)害教育部重點實驗室,湖北宜昌 443002;2.三峽地區(qū)地質(zhì)災(zāi)害與生態(tài)環(huán)境湖北省協(xié)同創(chuàng)新中心,湖北宜昌 443002)

強(qiáng)震作用下蘇家河口面板堆石壩動力響應(yīng)

王樂華1,2,鄧華鋒1,2,湯開宇1,朱 敏1,彭云楓2

(1.三峽大學(xué)三峽庫區(qū)地質(zhì)災(zāi)害教育部重點實驗室,湖北宜昌 443002;2.三峽地區(qū)地質(zhì)災(zāi)害與生態(tài)環(huán)境湖北省協(xié)同創(chuàng)新中心,湖北宜昌 443002)

為了分析蘇家河口面板堆石壩壩體在強(qiáng)震作用下的安全性能,建立了面板堆石壩的整體三維有限元模型和面板、趾板子模型,采用Hardin-Drnevich本構(gòu)模型,對壩體在強(qiáng)震作用下的動力響應(yīng)規(guī)律進(jìn)行了詳細(xì)分析。結(jié)果表明,在0.38g的強(qiáng)震作用下,面板堆石壩的加速度、動位移、動應(yīng)力以及垂直縫和周邊縫的動變形分布符合一般規(guī)律,未發(fā)現(xiàn)有特殊不利的現(xiàn)象,整體上滿足抗震要求,但是壩體頂部振動的“鞭梢”效應(yīng)比較明顯,在3/4壩高以上壩頂區(qū)域以及靠近下游坡面處存在較小的拉應(yīng)力區(qū)域。建議加強(qiáng)3/4壩高以上區(qū)域面板、壩頂及下游護(hù)坡的抗震工程措施,以保證大壩在強(qiáng)震下的安全。

面板堆石壩;三維有限元模型;強(qiáng)震;動力響應(yīng);蘇家河口

混凝土面板堆石壩以其安全性、經(jīng)濟(jì)性、適應(yīng)性等優(yōu)點受到國內(nèi)外壩工界的普遍重視[1-2]。目前很多面板堆石壩位于強(qiáng)震區(qū),使得面板堆石壩抗震問題研究的重要性和迫切性越來越突出[3-5]。蘇家河口面板堆石壩位于云南省境內(nèi),為一梯級開發(fā)電站,最大壩高為131.49m,水庫正常蓄水位為1590.00m,總庫容為2.25億m3。根據(jù)DL5180—2003《水電樞紐工程等級劃分及設(shè)計安全標(biāo)準(zhǔn)》,該水電站工程等別為Ⅱ等工程,工程規(guī)模為大(2)型,根據(jù)場地地震安全性評價,工程區(qū)50 a超越概率10％的地震動基巖峰值加速度為0.19g,50 a超越概率5％的地震動基巖峰值加速度為0.25g,100 a超越概率2％的地震動基巖峰值加速度為0.38g,基本烈度為Ⅷ度,工程按8度設(shè)防。對于該類型的工程,一般按照基本烈度進(jìn)行設(shè)計計算(有的工程在此基礎(chǔ)上考慮50 a超越概率5％的校核地震),但很少考慮強(qiáng)震的作用。汶川地震時,紫坪鋪面板堆石壩的壩基加速度在0.5g以上,達(dá)到了Ⅸ度[6],因此,有必要對強(qiáng)震區(qū)的壩體進(jìn)行強(qiáng)震作用下動力分析。本文研究的工程場地位于地震活動強(qiáng)烈的北西向騰沖—龍陵以西約20km,屬于區(qū)域構(gòu)造穩(wěn)定性較差地區(qū),因此,本文對壩體在強(qiáng)震作用下的動力響應(yīng)進(jìn)行了詳細(xì)計算分析,對其安全性做出評價,并提出加強(qiáng)壩體抗震安全性能的工程措施,為面板堆石壩的抗震設(shè)計提供科學(xué)依據(jù)。

1 動力分析方法

1.1 堆石料動力本構(gòu)模型及參數(shù)

地震作用是一種非等幅、非等周期的不規(guī)則荷載,在一次地震中,土石料將經(jīng)歷數(shù)十次甚至上百次卸載和再加載過程,并且它們之間無規(guī)律可循,給動力分析造成了一定困難。為了解決分析困難,比較常用的方法是應(yīng)用Mashing規(guī)則制定一個應(yīng)力應(yīng)變關(guān)系的骨架曲。在此基礎(chǔ)上的本構(gòu)模型主要有雙線性、黏彈性和彈塑性等模型,其中等效線性模型因具有概念明確、應(yīng)用方便的優(yōu)點,在有限元動力計算中得到廣泛運用[7-8]。本文采用Hardin-Drnevich雙曲線本構(gòu)模型進(jìn)行計算,相關(guān)的動力計算參數(shù)見表1。

表1 Hardin-Drnevich模型動力計算參數(shù)

地基計算條件:在動力計算中,基巖均采用線彈性模型,動彈性模量在靜彈性模量的基礎(chǔ)上提高30％,即靜彈性模量乘以系數(shù)1.3。只考慮壩基的彈性作用,消除壩基對振動的放大作用,采用無質(zhì)量壩基方案進(jìn)行分析。

1.2 加速度時程曲線

根據(jù)地震安全性評價報告,壩址區(qū)100a超越概率2％的地震動基巖峰值加速度為0.38g,地震動反應(yīng)譜特征周期0.35 s。計算中,根據(jù)抗震規(guī)范反應(yīng)譜擬合的人工地震波作為輸入地震波,并按3個方向輸入(順河向、橫河向和豎直向),計算地震歷時30 s,順河向、橫河向最大加速度為0.38g,豎直向最大加速度為水平向的2/3。地震波各方向加速度時程曲線如圖1所示。

1.3 計算流程

建立三維有限元模型,首先進(jìn)行靜力計算,壩體與覆蓋層材料采用鄧肯E-B雙曲線模型,并模擬分層填筑與蓄水過程,將靜力結(jié)果作為時程動力計算的初始狀態(tài)。動力計算中,壩體與覆蓋層材料采用Hardin-Drnevich模型,地震波三向基底輸入。大地震屬小概率事件,因此一般情況下,在地震計算中,將地震作用與水庫的正常蓄水位組合,在本文的動水壓力分析中,采用附加質(zhì)量法進(jìn)行考慮[9]。

2 有限元計算模型

圖1 輸入地震加速度時程曲線

面板堆石壩的垂直縫、周邊縫以及面板和墊層之間的接觸是大壩三維動力分析的關(guān)鍵位置,但是這些部位構(gòu)造尺寸與大壩整體尺寸相差懸殊,而且結(jié)構(gòu)尺寸、材料屬性和結(jié)構(gòu)邊界的變異性均較大。目前,對面板堆石壩的應(yīng)力變形分析通常采用整體模型進(jìn)行計算。分析時,對接觸部位及壩體同時進(jìn)行保證應(yīng)力計算精度的網(wǎng)格剖分難度較大,此時若均采用較密的網(wǎng)格劃分,會導(dǎo)致有限元模型規(guī)模非常龐大,求解困難;若采用較稀疏的網(wǎng)格則無法保證細(xì)節(jié)部位的結(jié)構(gòu)特征,上述關(guān)鍵部位的真實應(yīng)力變形狀態(tài)不能得到保證。如何既快又好地計算局部結(jié)構(gòu)的真實應(yīng)力和變形是一個值得研究的問題,以往的研究表明采用子模型進(jìn)行局部結(jié)構(gòu)的細(xì)化分析可以較好地解決此類問題[10-12]。

子模型方法又稱切割邊界位移法或特定邊界位移法,其主要思路是從整體模型中把重點研究的區(qū)域切割出來,作為子模型單獨進(jìn)行網(wǎng)格細(xì)分,并以切割邊界的位移計算結(jié)果作為子模型的邊界條件。由于子模型和整體模型的相對獨立性,因此可以增加子模型中的網(wǎng)格密度,以便對指定區(qū)域進(jìn)行更高精度計算,具體分析流程如圖2所示。在本文的分析中,建立了大壩的有限元整體模型、面板和趾板的子模型。

2.1 有限元整體模型

為考慮壩基及兩岸山體變形對混凝土面板堆石壩應(yīng)力變形的影響,本文計算范圍除了混凝土面板堆石壩以及壩基覆蓋層外,還包括了覆蓋層下的巖石基礎(chǔ)和兩岸山體,地基范圍在沖積層底部以下延伸150m,上下游方向和左右岸方向也各延伸150m。順?biāo)鞣较驗閤軸的正向;沿壩軸方向從右到左為y軸正向;豎直向上為z軸的正向(圖3)。三維網(wǎng)格剖分時主要采用8節(jié)點六面體單元和少數(shù)的6節(jié)點三棱柱過渡單元。單元總數(shù)為10762,節(jié)點總數(shù)為9253,模型底部、四周均為法向約束。

圖2 有限元子模型分析流程

圖3 面板堆石壩整體三維有限元網(wǎng)格

2.2 面板、趾板子模型

面板、趾板子模型如圖4所示,該子模型包括面板、垂直縫、周邊縫、趾板、面板和墊層的接觸面、墊層、過渡層及下部部分堆石料和基巖。在子模型中采用基于薄層單元的有厚度動力接觸單元進(jìn)行接觸面的模擬[13-14],Desai等[15]建議薄層單元的厚度t=(0.01～0.10)B(B表示薄層單元的寬度),根據(jù)殷宗澤等[16]的研究,在不產(chǎn)生顯著誤差的情況下,應(yīng)使厚度盡可能的小。本文計算中薄層單元的厚度取5cm,子模型單元總數(shù)為19 349,節(jié)點總數(shù)為17423。

圖4 子模型范圍和有限元網(wǎng)格

3 計算結(jié)果及分析

3.1 壩體動力響應(yīng)分析

3.1.1 壩頂動加速度反應(yīng)

壩頂沿壩軸線方向上各點的峰值絕對動加速度分布如圖5所示。壩頂3個方向的峰值絕對動加速度變化趨勢基本一致,以壩中線為對稱軸基本呈對稱分布,分布圖近似呈馬鞍形、有雙峰,壩軸線中點峰值絕對動加速度不是最大值。壩頂順河向峰值絕對動加速度在河床最大斷面附近達(dá)到最大值,最大值為8.24m/s2,動力放大倍數(shù)為2.21;壩頂豎直向峰值絕對動加速度最大值發(fā)生在左岸靠近河床部位,最大值為6.12m/s2,動力放大倍數(shù)為1.64;壩頂橫河向峰值絕對動加速度最大值發(fā)生在最大斷面右岸部位,最大值為6.26m/s2,動力放大倍數(shù)為1.68。

圖5 壩頂峰值絕對動加速度

壩體最大斷面的峰值絕對動加速度分布如圖6所示,順河向峰值絕對動加速度隨壩高增加而逐漸增大,最大值為8.13m/s2,動力放大倍數(shù)為2.18,發(fā)生在壩頂部位;豎直向峰值絕對動加速度亦隨壩高的增加而逐漸增大,最大值為6.78m/s2,動力放大倍數(shù)為1.82,發(fā)生在3/4壩高下游面;橫河向峰值絕對動加速度也是隨壩高增加而逐漸增大,在壩體中部上下游面產(chǎn)生了2個加速度放大區(qū)域,最大值為6.44m/s2,動力放大倍數(shù)為1.73,發(fā)生在3/4壩高下游面。

圖6 壩體最大斷面峰值絕對動加速度(單位:m/s2)

3.1.2 壩體地震動位移

壩體最大剖面的順河向、豎直向和橫河向峰值絕對動位移等值線如圖7所示。由圖7可見,最大剖面順河向峰值絕對動位移隨著壩體的升高而增大,最大位移發(fā)生在壩頂位置,最大值為13.70cm;最大剖面豎直向峰值絕對動位移也隨著壩體的升高而增大,最大位移發(fā)生在3/4壩高下游面,最大值為6.30cm;最大剖面橫河向峰值絕對動位移分布規(guī)律和豎直向峰值絕對動位移分布規(guī)律相似,最大動位移發(fā)生在壩頂靠下游面,最大值為10.10cm。

3.1.3 壩體地震動應(yīng)力

壩體最大剖面順河向、豎直向和橫河向峰值動拉應(yīng)力等值線如圖8所示。順河向峰值動拉應(yīng)力在2/3壩高以下的分布規(guī)律是由壩坡向壩軸線逐漸增大,在2/3壩高以上的拉應(yīng)力區(qū)的分布規(guī)律是由壩坡向壩軸線逐漸減小,最大值也發(fā)生在此區(qū)域,最大值為191.55kPa;豎直向峰值動拉應(yīng)力越靠近壩基越大,最大值為329.52kPa,發(fā)生在壩基靠下游的地方;橫河向峰值動拉應(yīng)力也基本是越靠近壩基越大,最大值為143.07kPa,發(fā)生在壩基靠下游的部位。

圖7 壩體最大剖面峰值絕對動位移(單位:cm)

圖8 壩體最大剖面峰值動拉應(yīng)力等值線(單位:kPa)

壩體動拉應(yīng)力最大值與靜應(yīng)力疊加后的等值線如圖9所示。與靜應(yīng)力疊加后,在順河向和豎直向壩體應(yīng)力基本呈受壓趨勢,但在順河向,在壩頂附近的上游堆石區(qū)以及靠近下游坡面處有少量拉應(yīng)力出現(xiàn)。

3.2 面板動力響應(yīng)分析

3.2.1 面板地震動加速度反應(yīng)

面板順坡向、垂直向和橫河向峰值絕對動加速度等值線如10所示。面板順坡向峰值絕對動加速度由面板頂部向底部逐漸減小,面板最大順坡向峰值絕對動加速度發(fā)生在面板頂部偏右岸處,最大值為6.50m/s2,放大倍數(shù)為1.74;豎直向峰值絕對動加速度也是由面板頂部向底部逐漸減小,最大值為6.37m/s2,放大倍數(shù)為1.71,發(fā)生在面板頂部偏左岸處;面板橫河向峰值絕對動加速度分別在面板兩岸形成2個加速度放大區(qū),最大值發(fā)生在面板頂部偏右岸處,最大值為6.11m/s2,放大倍數(shù)為1.64。

3.2.2 面板地震動位移

面板的動位移很大程度上取決于堆石體的動位移,3個方向的峰值絕對動位移等值線如圖11所示。順坡向峰值絕對動位移分布規(guī)律比較對稱,形成規(guī)則的水波紋,在面板頂中部位置達(dá)到最大,最大值為10.90cm;垂直向峰值絕對動位移最大值同樣發(fā)生在面板頂中部位置,為9.16cm;橫河向峰值絕對動位移在面板頂部偏左岸達(dá)到最大,為9.84cm。

圖9 壩體最大剖面動拉應(yīng)力與靜應(yīng)力疊加后等值線(單位:kPa)

圖10 面板峰值絕對動加速度(單位:m/s2)

圖11 面板峰值絕對動位移(單位:cm)

3.2.3 面板地震動應(yīng)力

面板順坡向、橫河向峰值動拉應(yīng)力等值線如圖12所示,順坡向峰值動拉應(yīng)力在面板中部3/4壩高的地方達(dá)到最大,向四周逐漸減小,最大動拉應(yīng)力為2295.71kPa;橫河向峰值動拉應(yīng)力在面板頂部右岸位置達(dá)到最大,最大動拉應(yīng)力為1630.78kPa。

圖12 面板峰值動拉應(yīng)力(單位:kPa)

3.3 面板垂直縫、周邊縫動變形響應(yīng)分析

面板垂直縫、周邊縫的動變形見表2,垂直縫、周邊縫的動張拉變形如圖13所示(圖中為了表示垂直縫不同高程張拉變形的變化,在面板上繪制了多條橫向輔助線,圖中所標(biāo)數(shù)值為垂直縫與輔助線交點位置處的張拉變形)。

由圖13可以看出,垂直縫、周邊縫的動張拉變形均以面板中心軸線呈對稱形式分布,其中垂直縫張拉變形在面板中心軸線一側(cè)同一高程上變形有先增大后減小的趨勢,且隨著高程的增加張拉位移不斷變大,最大張拉變形為7.287mm,出現(xiàn)在面板中心軸線右側(cè)頂部;周邊縫動張拉變形隨高程的增加呈現(xiàn)先變大后減小的趨勢,在面板頂部和底部張拉變形均較小,最大張拉變形為3.576mm,出現(xiàn)在面板中心軸線右側(cè)約1/2高程處。

表2 面板垂直縫、周邊縫的動變形

圖13 面板垂直縫、周邊縫動張拉變形(單位:mm)

4 結(jié)論及建議

a.在0.38g強(qiáng)震作用下,面板堆石壩的加速度、動位移、動應(yīng)力以及垂直縫和周邊縫的動變形分布符合一般規(guī)律,而其量值也比較合理,未發(fā)現(xiàn)有特殊不利的現(xiàn)象,整體滿足抗震要求,說明該面板堆石壩的設(shè)計合理。

b.地震作用下頂部振動的“鞭梢”效應(yīng)比較明顯,壩體頂部地震慣性力可能導(dǎo)致壩頂區(qū)堆石體的松動,堆石顆粒間咬合力喪失,從而有可能使壩頂和下游側(cè)發(fā)生壩面堆石表層滑動。

c.地震過程中,在壩高3/4以上區(qū)域面板的加速度、動位移、動拉應(yīng)力均較大,因此,建議加強(qiáng)壩高3/4以上區(qū)域面板、壩頂及下游護(hù)坡的抗震工程措施,以保證大壩抗震安全。

[1]陳厚群,徐澤平,李敏.汶川大地震和大壩抗震安全[J].水利學(xué)報,2008,39(10):1158-1167.(CHEN Houqun,XU Zeping,LI Min.Wenchuan Earthquake and seismic safety of large dams[J].Journal of Hydraulic Engineering,2008,39(10):1158-1167.(in Chinese))

[2]朱亞林,孔憲京,丁克偉,等.高土石壩的抗震研究進(jìn)展[J].安徽建筑工業(yè)學(xué)院院報:自然科學(xué)版,2009,17(3):1-6.(ZHU Yalin,KONG Xianjing,DING Kewei,et al.Research progress on earthquake resistant studies of high earth-rockfill dams[J].Journal of Anhui Institute of Architecture＆amp;Industry:Natural Science,2009,17(3):1-6.(in Chinese))

[3]趙劍明,常亞屏,陳寧.高心墻堆石壩地震變形與穩(wěn)定分析[J].巖土力學(xué),2004,25(增刊2):423-428.(ZHAO Jianming,CHANG Yaping,CHEN Ning.Study onearthquake-inducedpermanentdeformationand dynamic stability of high core rockfill dam[J].Rock and Soil Mechanics,2004,25(Sup2):423-428.(in Chinese))

[4]顧淦臣.土石壩地震工程[M].南京:河海大學(xué)出版社,1989.

[5]孔憲京,鄒德高,鄧學(xué)晶,等.高土石壩綜合抗震措施及其效果的驗算[J].水利學(xué)報,2006,37(12):1489-1495.(KONG Xianjing,ZOU Degao,DENG Xuejing,et al.Comprehensive earthquake resistant measure of high earth-rockfill dams and effectiveness verification[J].Journal of Hydraulic Engineering,2006,37(12):1489-1495.(in Chinese))

[6]張楚漢.汶川地震工程震害的啟示[J].水利水電技術(shù),2009(1):1-3.(ZHANG Chuhan.The Wenchuan Earthquake engineering disaster revelation[J].Water Resources and Hydropower Engineering,2009(1):1-3.(in Chinese))

[7]李宏男,肖詩云,霍林生.汶川地震震害調(diào)查與啟示[J].建筑結(jié)構(gòu)學(xué)報,2008,29(4):10-19.(LI Hongnan,XIAOShiyun,HUOLinsheng.Damage investigation and analysis of engineering structures in the WenchuanEarthquake[J].JournalofBuilding Structures,2008,29(4):10-19.(in Chinese))

[8]何蘊龍,劉俊林,熊堃.汶川地震冶勒大壩動力響應(yīng)規(guī)律分析[J].四川大學(xué)學(xué)報:工程科學(xué)版,2009,41(3):157-164.(HE Yunlong,LIU Junlin,XIONG Kun.SeismicresponseofYeleDamduringWenchuan Earthquake[J].JournalofSichuanUniversity:Engineering Science Edition,2009,41(3):157-164.(in Chinese))[9]趙劍明,常亞屏,陳寧.龍首二級面板堆石壩三維真非線性地震反應(yīng)分析和評價[J].巖土力學(xué),2004(增刊2):388-392,396.(ZHAO Jianming,CHANG Yaping,CHEN Ning.3-D authentic nonlinear dynamic analysis and evaluation of Longshou concrete faced rockfill dam[J].Rock and Soil Mechanics,2004(Sup2):388-392,396.(in Chinese))

[10]汪明元,程展林,林紹忠,等.高面板堆石壩應(yīng)力變形的三維子模型法研究[J].長江科學(xué)院院報,2005,22(5):49-51,69.(WANG Mingyuan,CHENG Zhanlin,LIN Shaozhong,et al.Study on 3D submodel method for stress and deformation analysis for high concrete raced rockfill dam[J].Journal of Yangtze River Scientific Research Institute,2005,22(5):49-51,69.(in Chinese))

[11]武甲慶,楊令強(qiáng),楊克坤,等.用子模型法研究拱壩裂縫的產(chǎn)狀[J].中國農(nóng)村水利水電,2009(9):108-111.(WU Jiaqing,YANG Lingqiang,YANG Kekun,et al.Submodels used to study strength cracking generation of arch dams[J].China Rural Water and Hydropower,2009(9):108-111.(in Chinese))

[12]BAZANT Z P,CEDOLIN L.Blunt crack band propagation infinite element analysis[J].Journal of the Engineering Mechanics Division,ASCE,1979,105(EM2):297-315.

[13]沈振中,徐志英.V形河谷中非均勻土石壩振動的簡化分析[J].水利學(xué)報,2002(3):74-79.(SHEN Zhenzhong,XU Zhiying.Simplified analysis of transverse vibration for inhomogeneous earth-rock dams in triangular canyons[J].Journal of Hydraulic Engineering,2002(3):74-79.(in Chinese))

[14]李遠(yuǎn)略,李中原.觀音巖水壩地震動力分析[J].河南科學(xué),2009,27(10):1301-1304.(LIYuanlue,LI Zhongyuan.Guanyinyan dynamic analysis on the dam[J].Henan Science,2009,27(10):1301-1304.(in Chinese))

[15]DESAI C S,ZAMAN M M,LIGHTNER J G,et al.Thinlayer element for interfaces and joints[J].International JournalforNumericalandAnalyticalMethodsin Geomechanics,1984,8:19-43.

[16]殷宗澤,朱泓,許國華.土與結(jié)構(gòu)材料接觸面的變形及其數(shù)學(xué)模擬[J].巖土工程學(xué)報.1994,16(3):14-22(YINZongze,ZHUHong,XUGuohua.Numerical simulation of the deformation in the interface between soil andstructuralmaterial[J].ChineseJournalof Geotechnical Engineering,1994,16(3):14-22.(in Chinese))

Dynamic response analysis of Sujiahekou concrete-face rock-fill dam under strong earthquake

//WANG Lehua1,2,DENG Huafeng1,2,TANG kaiyu1,ZHU Min1,PENG Yunfeng2(1.China Three Gorges University Kay Laboratory of Geological Hazards on Three Gorges Reservoir Area,Ministry of Education,Yichang 443002,China;2.Collaborative Innovation Center for Geo-Hazards and Eco-Environment in Three Gorges Area,Yichang 443002,China)

As the Sujiahekou concrete-face rock-fill dam(CFRD)is located next to the strong earthquake activity belt,the analyze and evaluation of the safety performance of the dam should be assessed.To archive this goal,this study established the three-dimensional finite element model of the CFRD and corresponding sub model of concrete face and toe board.We used the Hardin-Drnevich constitutive model to calculate and analyze the dynamic response rules of dam's body under strong earthquakes in detail.The results show that under the action of 0.38g earthquake,the acceleration,dynamic displacement,dynamic stress,and dynamic deformation of vertical seam and peripheral joint distribution conform to the general rules.There is no special adverse phenomenon;The CFRD can meet seismic requirements.However,the whipping effect of the dam body vibration is obvious and there is a small tensile stress area in 3/4 dam high above crest and areas close to the downstream slope.Therefore,seismic engineering measures,strengthening the 3/4 dam high above area panel,dam crest and downstream slope,have been suggested to ensure the safety of the dam in earthquake.

concrete-faced rock-fill dam;three-dimensional finite element model;strong earthquake;dynamic response;Sujiahekou

TV641.4

:A

:1006-7647(2014)06-0007-06

10.3880/j.issn.1006-7647.2014.06.002

2013-09-05 編輯:劉曉艷)

湖北省教育廳科學(xué)研究計劃基金(Q20091301);國家重點基礎(chǔ)研究發(fā)展計劃(973計劃)(2012CB426502);水利部公益性行業(yè)科研專項(201401029)

王樂華(1977—),男,安徽懷寧人,副教授,博士,主要從事巖土工程研究。E-mail:lehuatg@126.com