• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    高土石壩宏細(xì)觀壩水動(dòng)力流固耦合理論研究進(jìn)展

    2013-04-09 08:05:05岑威鈞陳亞南
    水利水電科技進(jìn)展 2013年6期
    關(guān)鍵詞:庫(kù)水細(xì)觀石壩

    岑威鈞,孫 輝,陳亞南

    (河海大學(xué)水利水電學(xué)院,江蘇南京 210098)

    高土石壩宏細(xì)觀壩水動(dòng)力流固耦合理論研究進(jìn)展

    岑威鈞,孫 輝,陳亞南

    (河海大學(xué)水利水電學(xué)院,江蘇南京 210098)

    高土石壩遭遇強(qiáng)烈地震時(shí)會(huì)與壩面庫(kù)水及壩內(nèi)孔隙水發(fā)生動(dòng)力流固耦合相互作用,分別從宏觀壩水耦合系統(tǒng)和細(xì)觀水土耦合系統(tǒng)兩個(gè)角度對(duì)高土石壩壩水動(dòng)力流固耦合理論的研究歷史、主要研究成果和研究趨勢(shì)進(jìn)行了闡述評(píng)價(jià)。對(duì)宏觀尺度的大壩與壩面庫(kù)水流固耦合作用,主要從早期壩水相互作用模型、庫(kù)水運(yùn)動(dòng)精細(xì)分析方法以及流固耦合系統(tǒng)坐標(biāo)描述等方面做出評(píng)述;對(duì)細(xì)觀尺度的壩內(nèi)水土耦合作用,主要從早期解耦或擬耦合的水土動(dòng)力分析方法、基于Boit動(dòng)力固結(jié)理論的細(xì)觀水土動(dòng)力流固耦合,以及基于混合物理論的細(xì)觀水土動(dòng)力流固耦合等方面做出評(píng)述。綜合兩個(gè)尺度的流固耦合作用研究現(xiàn)狀,建議對(duì)高土石壩建立水庫(kù)-土石壩-孔隙水的宏細(xì)觀整體動(dòng)力流固耦合系統(tǒng)進(jìn)行綜合研究,并指出了當(dāng)前高土石壩宏細(xì)觀壩水動(dòng)力流固耦合理論尚待繼續(xù)深入研究的若干相關(guān)問(wèn)題。

    高土石壩;庫(kù)水;孔隙水;宏觀壩水動(dòng)力流固耦合;細(xì)觀水土動(dòng)力流固耦合;動(dòng)力分析

    近年來(lái)我國(guó)高土石壩的數(shù)量和高度均有大幅度提升。據(jù)不完全統(tǒng)計(jì),國(guó)內(nèi)已建、在建及擬建(規(guī)劃)的壩高200 m以上的高土石壩就有10余座,個(gè)別壩高直逼或超過(guò)300 m,如大渡河雙江口(壩高312 m)、雅礱江兩河口(壩高305 m)、怒江松塔(壩高307 m)、瀾滄江如美(壩高315 m)等。我國(guó)是一個(gè)多地震國(guó)家,地震活動(dòng)頻度高、強(qiáng)度大、震源淺、分布廣,震災(zāi)嚴(yán)重,而高土石壩所在的西部地區(qū)又是主要的強(qiáng)地震區(qū)。2008年的“5·12”汶川大地震中有兩座高土石壩受到了地震影響,其中,105.3 m高的碧口心墻土石壩地震損壞輕微,156 m高的紫坪鋪面板堆石壩出現(xiàn)了面板壓碎錯(cuò)位、壩頂較大震陷等較嚴(yán)重的地震損傷。紫坪鋪面板堆石壩也是目前世界上遭遇最強(qiáng)烈地震考驗(yàn)的最高土石壩。相對(duì)高混凝土壩而言,國(guó)內(nèi)外對(duì)高土石壩抗震試驗(yàn)和抗震理論等方面的研究相對(duì)偏少,尤其對(duì)那些擬建的300 m級(jí)超高土石壩,由于壩高、水庫(kù)規(guī)模和設(shè)計(jì)地震加速度又創(chuàng)新高,強(qiáng)震作用下大壩和水體(庫(kù)水和孔隙水)的動(dòng)力相互作用可能非常強(qiáng)烈,對(duì)大壩抗震安全性影響重大,需進(jìn)行深入系統(tǒng)的專項(xiàng)研究。

    高土石壩遭遇強(qiáng)烈地震時(shí)會(huì)和水體(庫(kù)水和孔隙水)發(fā)生動(dòng)力相互作用,即壩水動(dòng)力流固耦合,可分為兩種類(lèi)型:第一類(lèi)是壩水動(dòng)力流固耦合發(fā)生在庫(kù)水與大壩(含地基,下同)的交界面上,兩者在宏觀上是相互分離的,比如巖基上的面板堆石壩;另一類(lèi)是水與壩體土石料無(wú)法宏觀分開(kāi),即壩內(nèi)孔隙水和土骨架之間發(fā)生細(xì)觀水土動(dòng)力耦合作用。地震時(shí),庫(kù)水與壩面土體發(fā)生宏觀水土相互作用后,庫(kù)水進(jìn)一步入滲到土石壩中或孔隙水流出壩體進(jìn)入水庫(kù),形成動(dòng)力非穩(wěn)定滲流場(chǎng),此時(shí)壩體內(nèi)孔隙水與土骨架之間進(jìn)行細(xì)觀水土動(dòng)力耦合作用,兩者間進(jìn)行應(yīng)力(孔壓)的傳遞以及變形的協(xié)調(diào)變化。在排水條件下,由于孔隙水的動(dòng)力滲流作用,土體中所含的孔隙水量要發(fā)生變化,振動(dòng)超靜孔隙水壓力還要進(jìn)一步擴(kuò)散和消散,最終趨于零,土體變形至動(dòng)力固結(jié)穩(wěn)定。(斜)心墻土石壩遭遇地震時(shí)與孔隙水的動(dòng)力相互作用就是屬于這類(lèi)細(xì)觀水土動(dòng)力流固耦合。

    強(qiáng)震作用下,高土石壩自身變形較大,尤其在壩頂附近,地震反應(yīng)明顯。高壩往往形成大水庫(kù),在地震時(shí)庫(kù)水會(huì)發(fā)生激蕩,在水庫(kù)表面尤為激烈。此時(shí),壩庫(kù)兩者間可能會(huì)引發(fā)相互激振,壩體、庫(kù)水和孔隙水會(huì)發(fā)生宏細(xì)觀動(dòng)力相互作用,進(jìn)而可能引發(fā)心墻動(dòng)水劈裂、面板動(dòng)力破損、壩基液化和動(dòng)孔壓引起的壩坡動(dòng)力失穩(wěn)、土骨架的動(dòng)力滲透破壞等一系列大壩安全問(wèn)題,涉及水工結(jié)構(gòu)、巖土力學(xué)、振動(dòng)力學(xué)、流體力學(xué)等多門(mén)學(xué)科。目前,國(guó)內(nèi)外學(xué)者對(duì)土石壩與水體的動(dòng)力流固耦合理論以及分析方法開(kāi)展了較為廣泛的研究,取得了一些成果。本文從宏觀和細(xì)觀兩方面對(duì)土石壩與水體的動(dòng)力流固耦合研究歷史、主要理論研究進(jìn)展和成果進(jìn)行闡述和評(píng)價(jià),主要涉及壩庫(kù)宏觀流固耦合、水庫(kù)運(yùn)動(dòng)分析方法、壩水細(xì)觀流固耦合等內(nèi)容,最后指出尚待深入研究的相關(guān)問(wèn)題。

    1 土石壩-水庫(kù)宏觀動(dòng)力流固耦合

    1.1 早期基本的壩水相互作用模型

    對(duì)于混凝土面板堆石壩,面板與水庫(kù)之間存在宏觀壩水交界面,壩水相互作用機(jī)理比較明確,其動(dòng)力耦合分析理論和方法主要參照混凝土壩。目前最為簡(jiǎn)單實(shí)用的分析方法依舊是Westergaard動(dòng)水壓力附加質(zhì)量法。自1933年Westergaard[1]率先發(fā)表了關(guān)于動(dòng)水壓力的論文以后,許多研究者針對(duì)Westergaard模型的一些限制條件進(jìn)行了討論分析并不斷完善,現(xiàn)常被引入到面板壩壩水動(dòng)力作用的分析中。Westergaard附加質(zhì)量法將宏觀壩水相互作用簡(jiǎn)化為僅水庫(kù)對(duì)大壩的單向作用,忽略大壩變形引起水域形狀的改變,從而省卻了流體計(jì)算域的建模求解問(wèn)題。事實(shí)上,Westergaard法需要滿足壩體和地基為剛性等假定以及用一些近似方法來(lái)模擬水體對(duì)壩體的單向動(dòng)力作用,這與實(shí)際土石壩庫(kù)水相互作用情況有較大差異。另外,對(duì)于高面板壩和強(qiáng)烈地震,只有考慮水庫(kù)運(yùn)動(dòng)特性才能真正體現(xiàn)土石壩-水庫(kù)的雙向動(dòng)力流固耦合作用。早期對(duì)庫(kù)水運(yùn)動(dòng)建模主要采用考慮水體壓縮性的波動(dòng)方程和忽略水體壓縮性的Laplace方程,這兩類(lèi)方程均假定水體是無(wú)黏性的理想流體。張振國(guó)[2]最早開(kāi)展面板堆石壩壩水相互作用的相關(guān)研究,將水體視為不可壓縮流體,采用Laplace方程,對(duì)上游面不同坡度的面板壩受動(dòng)水壓力的作用效應(yīng)進(jìn)行了計(jì)算分析,當(dāng)壩坡較陡時(shí)(坡度1∶1.3),動(dòng)水壓力將使面板動(dòng)應(yīng)力增大40%左右;即使較緩壩坡(坡度1∶2.0),也能使面板動(dòng)應(yīng)力增大10%左右。這種對(duì)地震中水庫(kù)運(yùn)動(dòng)特性的考慮依舊采用將庫(kù)水對(duì)壩體的作用轉(zhuǎn)化成附加質(zhì)量的形式。遲世春等[3]在此基礎(chǔ)上考慮了水體壓縮性的影響,用幾種不同的庫(kù)水模型(不可壓縮的Laplace方程、可壓縮的波動(dòng)方程)對(duì)面板壩自振頻率的影響進(jìn)行對(duì)比分析,并對(duì)水體計(jì)算域長(zhǎng)度的選取等方面進(jìn)行了細(xì)致的分析研究,得到一些有益的結(jié)論。隨后,遲世春等[4]進(jìn)一步分析了不同動(dòng)水壓力模型對(duì)面板壩動(dòng)力反應(yīng)的影響,得到的結(jié)論是不可壓縮水體模型與Westergaard模型動(dòng)力反應(yīng)接近,而可壓縮水體模型與不考慮動(dòng)水壓力情況接近。此外,遲世春等[5]還詳細(xì)研究了不同形式的動(dòng)水壓力附加質(zhì)量矩陣,討論了采用集中陣或分布陣對(duì)百米高面板堆石壩地震動(dòng)力反應(yīng)以及幅頻反應(yīng)的影響。鄧海峰[6]以董箐混凝土面板堆石壩為例,用附加質(zhì)量法分析了150 m級(jí)的高面板堆石壩的動(dòng)力反應(yīng),結(jié)果顯示動(dòng)水壓力使得壩體地震響應(yīng)減小,但是對(duì)面板的動(dòng)應(yīng)力以及殘余應(yīng)力影響較大。

    1.2 更為精細(xì)水體模型的建立

    上述宏觀壩水動(dòng)力流固耦合中,水體建模較為簡(jiǎn)單,沒(méi)有考慮水庫(kù)表面重力波等影響,對(duì)于地震烈度和壩高不大的面板壩,基本能體現(xiàn)“土石壩-庫(kù)水”雙向動(dòng)力相互作用。對(duì)于300 m級(jí)的高土石壩,強(qiáng)震作用下壩頂部位的變形達(dá)數(shù)十厘米甚至1 m以上,且上游壩坡面加速度明顯增大,對(duì)水庫(kù)運(yùn)動(dòng)和局部邊界影響顯著,因此需要更精確地考慮地震時(shí)水庫(kù)的運(yùn)動(dòng)特性,以合理反映水庫(kù)與壩體之間動(dòng)力耦合效應(yīng)。此時(shí)需將水庫(kù)連同大壩一起進(jìn)行細(xì)致建模和耦合求解,普通的Laplace方程與波動(dòng)方程對(duì)水庫(kù)運(yùn)動(dòng)描述就顯得過(guò)于簡(jiǎn)單,需引入更為嚴(yán)格的水體運(yùn)動(dòng)數(shù)學(xué)模型,即Navier-Stokes方程(N-S方程)。由于N-S方程由非線性偏微分方程耦合而成,求解析解是非常困難的,結(jié)合水庫(kù)運(yùn)動(dòng)特點(diǎn)可做適當(dāng)簡(jiǎn)化[7]。目前在分析結(jié)構(gòu)與水體相互作用中應(yīng)用較為廣泛的水體模型是基于勢(shì)的亞音速流體模型[8-10]與N-S流體模型。勢(shì)流體模型需要流體符合無(wú)旋、無(wú)黏、無(wú)熱轉(zhuǎn)化的假定[11]?;谒俣葎?shì)的方程相對(duì)于N-S方程具有更少的未知量,在實(shí)際工程計(jì)算中顯得簡(jiǎn)單快捷,因此更適合于土石壩-庫(kù)水動(dòng)力流固耦合分析。Sussman等[12]首先將基于勢(shì)的亞音速公式應(yīng)用于流固耦合分析;王偉華等[13]以重力壩為例,對(duì)比了傳統(tǒng)附加質(zhì)量模型與勢(shì)流體模型下壩體的動(dòng)力反應(yīng),指出傳統(tǒng)附加質(zhì)量模型結(jié)果偏大,進(jìn)行地震反應(yīng)分析時(shí)偏于保守;N-S方程由于其本身的復(fù)雜性,目前在以求解結(jié)構(gòu)反應(yīng)為主的土石壩壩水動(dòng)力分析中應(yīng)用不多,主要應(yīng)用在潰壩以及溢洪道水流數(shù)值模擬中[14];劉金云等[15]以二維壩水相互作用為例,對(duì)勢(shì)流體模型與N-S流體模型進(jìn)行了比較研究,結(jié)果顯示在滿足基于勢(shì)的亞音速公式假定情況下,兩者的計(jì)算結(jié)果較為相似,且選擇基于勢(shì)的亞音速公式具有更少的自由度,相對(duì)于基于N-S方程的模型更為快捷有效。

    1.3 耦合系統(tǒng)的坐標(biāo)描述及先進(jìn)分析方法的選擇

    大壩與水庫(kù)耦合求解時(shí),大壩固體域采用Lagrange坐標(biāo)系,而流體域采用Lagrange坐標(biāo)系或Euler坐標(biāo)系。但是由于流體自身特性,對(duì)于水庫(kù)的建模無(wú)論是Lagrange坐標(biāo)系還是Euler坐標(biāo)系均有明顯不足。ALE方法綜合了Euler坐標(biāo)系和Lagrange坐標(biāo)系的優(yōu)點(diǎn),可用于水庫(kù)帶自由表面和上游壩面邊界變化較大的水庫(kù)運(yùn)動(dòng),克服了Lagrange方法常見(jiàn)的網(wǎng)格畸變等問(wèn)題。20世紀(jì)80年代發(fā)展起來(lái)的迎風(fēng)格式有限元與分步格式有限元能夠很好地解決水庫(kù)流體數(shù)值解的失真震蕩現(xiàn)象,使基于ALE描述的迎風(fēng)有限元法和分步有限元法應(yīng)運(yùn)而生;它采用ALE描述從而可以精確確定流體邊界的位置,且不會(huì)引起網(wǎng)格糾纏,同時(shí)引入迎風(fēng)格式或分步方法來(lái)消除對(duì)流效應(yīng)引起的非物理振蕩,目前其已在流體結(jié)構(gòu)相互作用方面得到應(yīng)用,但在大壩-庫(kù)水相互作用方面的應(yīng)用目前還不多見(jiàn)。Ramaswamy等[16]運(yùn)用ALE分步有限元來(lái)解決不可壓縮黏性流體自由表面的運(yùn)動(dòng)問(wèn)題。Souli等[17]將ALE技術(shù)和GLS迎風(fēng)有限元法相結(jié)合應(yīng)用到求解帶有自由液面流體的大幅晃動(dòng)問(wèn)題。Takase等[18]將瞬變理論同 SUPG迎風(fēng)有限元結(jié)合起來(lái),利用SUPG瞬變有限元方法來(lái)計(jì)算海岸的淺水波問(wèn)題。岳寶增[19-20]等利用ALE分步有限元算法來(lái)解決三維液體的大幅度晃動(dòng)問(wèn)題。華蕾娜[21]對(duì)水池中的自由表面波利用ALE分步有限元進(jìn)行了模擬。陳文元等[22]運(yùn)用ALE描述將流體域的網(wǎng)格節(jié)點(diǎn)按照自由液面的運(yùn)動(dòng)和耦合面的移動(dòng)不斷更新,模擬了壩體在地震作用下的動(dòng)力特性、庫(kù)水自由表面重力波影響以及庫(kù)水域有效影響范圍的問(wèn)題。上述研究成果可以很好地借鑒應(yīng)用到強(qiáng)震下高土石壩與水庫(kù)的動(dòng)力流固耦合中,以精細(xì)反映水庫(kù)的激振和庫(kù)面運(yùn)動(dòng)等特性。

    壩水動(dòng)力相互作用涉及固體域和流體域的聯(lián)合求解,邊界條件非線性程度高,目前只能采用數(shù)值解法,如有限差分法、有限體積法、邊界元法和有限元法等。王國(guó)輝等[23]詳細(xì)分析了各數(shù)值方法的利弊,其中有限元法相對(duì)于其他方法而言能比較容易處理各種復(fù)雜的幾何邊界條件,在很多情況下能得到較高的精度,因此被ADINA、ABAQUS等大型通用有限元軟件在求解流固耦合問(wèn)題時(shí)廣泛采用。有限元法依舊是目前求解“大壩-水庫(kù)”動(dòng)力耦合問(wèn)題的首選數(shù)值解法。

    2 土石壩土骨架-孔隙水細(xì)觀動(dòng)力流固耦合

    2.1 早期基本的細(xì)觀水土動(dòng)力流固耦合

    筑壩土石料屬于典型的多孔介質(zhì),由土骨架、孔隙水和孔隙氣組成。蓄水后浸潤(rùn)線以下的壩體完全被孔隙水充滿,在地震作用下會(huì)產(chǎn)生振動(dòng)孔隙水壓力。為了模擬地震過(guò)程中孔隙水壓力的變化,Seed等最早建議了一個(gè)用于計(jì)算動(dòng)孔壓增長(zhǎng)的解耦模型,適用于在地震歷時(shí)較短且假設(shè)壩體不對(duì)外排水的情況。這種解耦的動(dòng)孔壓模型初步實(shí)現(xiàn)了孔隙水對(duì)土骨架的細(xì)觀動(dòng)力作用。其后徐志英等[24-25]在此基礎(chǔ)上結(jié)合Biot靜力固結(jié)方程來(lái)考慮動(dòng)孔壓的擴(kuò)散和消散,發(fā)展了土石壩排水有效應(yīng)力動(dòng)力分析方法。這類(lèi)方法將單獨(dú)的動(dòng)孔壓增長(zhǎng)、土石壩運(yùn)動(dòng)與土體靜力固結(jié)方程聯(lián)系起來(lái),在小時(shí)段內(nèi)仿照靜力固結(jié)耦合問(wèn)題來(lái)處理動(dòng)孔壓的擴(kuò)散與消散,通過(guò)與大壩動(dòng)力反應(yīng)分時(shí)段交替計(jì)算來(lái)考慮兩者相互作用,近似實(shí)現(xiàn)孔隙水壓力的增長(zhǎng)、擴(kuò)散和消散及與土骨架間的動(dòng)力相互作用,但未從本質(zhì)上描述土骨架與孔隙水兩者間真正的細(xì)觀動(dòng)力流固耦合過(guò)程。為了更好地解決這一問(wèn)題,Biot動(dòng)力固結(jié)理論和基于連續(xù)介質(zhì)力學(xué)的混合物理論便應(yīng)運(yùn)而生。

    2.2 基于Boit動(dòng)力固結(jié)理論的細(xì)觀水土動(dòng)力流固耦合

    Biot等[26-29]對(duì)水土兩相介質(zhì)的相互作用機(jī)理(完全耦合)進(jìn)行了開(kāi)創(chuàng)性的研究,率先建立了飽和土體線彈性多孔介質(zhì)平衡方程,后又給出了系統(tǒng)動(dòng)能和介質(zhì)衰減函數(shù)表達(dá)式,建立了通過(guò)慣性項(xiàng)和黏性項(xiàng)耦合起來(lái)的系統(tǒng)動(dòng)力方程,發(fā)展了含有可壓縮黏性流體的多孔彈性固體應(yīng)力波傳播理論,后又將其推廣至各向異性、黏彈性以及包含固相熱量耗散的飽和兩相多孔介質(zhì)中。Ghaboussi等[30-31]在Biot動(dòng)力方程的基礎(chǔ)上依據(jù)變分原理建立了有限元方程,分析了動(dòng)荷載作用下飽和多孔半空間土體的瞬態(tài)反應(yīng)和土石壩在平面應(yīng)變情況下的地震瞬態(tài)反應(yīng)。

    Biot波動(dòng)理論可以正確考慮飽和土體中土骨架與孔隙水之間的相互作用,但是Biot動(dòng)力方程中的彈性常數(shù)與慣性耦合系數(shù)難以測(cè)定,限制了該理論的推廣應(yīng)用。為此,門(mén)福錄[32]在假定孔隙水為不可壓縮、固相骨架為彈性的條件下,依據(jù)Biot準(zhǔn)靜力情形下的方程再附加以慣性項(xiàng)建立了動(dòng)力學(xué)方程組。盛虞等[33]根據(jù)有效應(yīng)力原理推導(dǎo)出土體二維動(dòng)力固結(jié)方程,將其與孔隙水壓力計(jì)算結(jié)合,對(duì)土壩進(jìn)行考慮孔隙水壓力產(chǎn)生、擴(kuò)散與消散的有效應(yīng)力動(dòng)力反應(yīng)分析。林本海[34]在文獻(xiàn)[33]的基礎(chǔ)上將其理論推廣到三維問(wèn)題,分析中以動(dòng)力固結(jié)方程為基礎(chǔ),在震動(dòng)過(guò)程中全程跟蹤孔隙水壓力產(chǎn)生、擴(kuò)散和消散的發(fā)展變化,將動(dòng)力滲流與土體動(dòng)力反應(yīng)分析相耦合,較好地反映了土體震動(dòng)過(guò)程中的實(shí)際狀態(tài)。但是由于林本海采用的動(dòng)本構(gòu)模型的限制,在動(dòng)力微分方程中動(dòng)孔壓仍然沿用過(guò)去動(dòng)力反應(yīng)與動(dòng)力固結(jié)分離計(jì)算時(shí)的方法,使得由動(dòng)孔壓模型計(jì)算出的動(dòng)孔壓與動(dòng)力滲流固結(jié)引起的孔壓出現(xiàn)了矛盾。為此,需尋求真正耦合且實(shí)用性強(qiáng)的多孔介質(zhì)動(dòng)力流固耦合理論。Zienkiewicz[35]對(duì) Biot多孔介質(zhì)模型的波動(dòng)問(wèn)題進(jìn)行了進(jìn)一步的研究,考慮了孔隙度以及各相密度的變化,并增加了固體和流體的慣性項(xiàng),能合理反映地震過(guò)程土石壩中土骨架與動(dòng)孔隙水的動(dòng)力相互作用。Zienkiewicz等[36-39]對(duì)飽和多孔介質(zhì)建立了用不同未知量表示的幾種有限元方程形式,即以固相位移u和液相相對(duì)位移w為基本未知量的“u-w”形式,以及固相位移u和孔隙水壓力p為基本未知量的“u-p”形式;由于“u-p”形式計(jì)算結(jié)果的精度略遜于“u-w”形式,又進(jìn)一步給出了以“u-w”形式波動(dòng)方程為基礎(chǔ)的高階有限元法、以“u-w-p”為基本未知量和以“u-w-p-σ”為基本未知量的混合有限元法及基于“u-w”形式波動(dòng)方程的Hermitean方法的有限元方程式。李宏儒[40]在林本海研究成果基礎(chǔ)上采用有效應(yīng)力物態(tài)動(dòng)本構(gòu)關(guān)系,利用瞬態(tài)理論,舍去了孔壓模型的引入,對(duì)動(dòng)力滲流和動(dòng)力固結(jié)相耦合的土體有效應(yīng)力計(jì)算方法進(jìn)行了進(jìn)一步的分析和改進(jìn)。劉凱欣等[41]將飽和多孔介質(zhì)的固相和液相處理成完全獨(dú)立的兩相,通過(guò)二者交界面處的流固耦合作用相互聯(lián)系,給出了固相和液相的基本方程以及二者界面耦合關(guān)系方程,開(kāi)發(fā)了三維流固混合顯式動(dòng)力有限元計(jì)算程序,對(duì)飽和多孔介質(zhì)中應(yīng)力波的傳播進(jìn)行了數(shù)值模擬,并詳細(xì)討論了孔隙率和孔隙形狀等因素對(duì)應(yīng)力波傳播主導(dǎo)波形的影響。上述理論和方法為研究強(qiáng)震作用下高土石壩與孔隙水的動(dòng)力耦合效應(yīng)分析開(kāi)辟了嶄新途徑,勢(shì)必會(huì)不斷得到應(yīng)用和驗(yàn)證提高。

    2.3 基于混合物理論的細(xì)觀水土動(dòng)力流固耦合

    在眾多學(xué)者研究土骨架與孔隙水相互作用的Biot動(dòng)力固結(jié)理論的同時(shí),解決土體與孔隙水耦合問(wèn)題的另一種理論——混合物理論也逐步得到了深入研究。混合物理論以熱力學(xué)理論為基礎(chǔ),對(duì)單一物質(zhì)連續(xù)系統(tǒng)理論進(jìn)行了拓展,具有良好的自適性和系統(tǒng)性。Truesdell[42]提出了任意組分混合物的質(zhì)量、動(dòng)量和能量的局部平衡方程,標(biāo)志著現(xiàn)代混合物多孔介質(zhì)理論(PMT)研究的開(kāi)始。Prevost[43-44]提出了一種飽和多孔介質(zhì)波動(dòng)理論的有限元數(shù)值解法,其中,土骨架可以采用非線性或彈塑性本構(gòu)模型,也可以考慮大變形問(wèn)題,液相可以假定為可壓縮或不可壓縮。為了去掉由于剛性流體的存在而產(chǎn)生的對(duì)時(shí)間步長(zhǎng)的限制,采用隱-顯式積分算法。Yiagos等[45]建立了一種可用于土壩彈塑性地震反應(yīng)分析的簡(jiǎn)單而有效的二維有限元計(jì)算方法,其中,壩剖面近似為對(duì)稱的三角形,土層為飽和多孔介質(zhì)且水平分層,并考慮了水的存在,將土骨架按照非線性滯變體進(jìn)行處理。嚴(yán)波等[46]采用基于混合物理論的兩相多孔介質(zhì)模型,建立了黏性流體飽和兩相多孔介質(zhì)非線性動(dòng)力問(wèn)題的控制方程,利用Galerkin加權(quán)殘值法推導(dǎo)了有限元方程組,并采用隱式Newmark法進(jìn)行求解。秦小軍等[47]根據(jù)流固兩相混合物的連續(xù)介質(zhì)力學(xué)理論,采用Galerkin加權(quán)殘值法,選取固相位移、液相位移和孔壓作為場(chǎng)變量,對(duì)固液兩相耦合方程組進(jìn)行有限元離散,得到解耦的方程組,然后在時(shí)域上采用Wilson-θ法進(jìn)行逐步積分,得到一種分析二維飽和多孔介質(zhì)地震反應(yīng)的三場(chǎng)有限元法。

    基于連續(xù)介質(zhì)力學(xué)的混合物理論將運(yùn)動(dòng)學(xué)、動(dòng)力學(xué)、熱力學(xué)及本構(gòu)理論融為一體,包含了各種復(fù)雜的因素,可以更加全面地反映地震過(guò)程中土骨架、孔隙水甚至是孔隙氣之間的動(dòng)力耦合作用?;诨旌衔锢碚摰亩嗫捉橘|(zhì)模型能夠蛻化為經(jīng)典的Biot模型[48]。兩種理論在物理和數(shù)學(xué)上均有很好的一致性[49],能滿足理論上的精確性。相信兩種理論會(huì)在強(qiáng)震下土石壩與孔隙水的細(xì)觀動(dòng)力耦合問(wèn)題中不斷得到深入研究和應(yīng)用。

    3 水庫(kù)-土石壩-孔隙水宏細(xì)觀動(dòng)力流固耦合

    由于目前尚不能在實(shí)驗(yàn)室或現(xiàn)場(chǎng)對(duì)土石壩和水庫(kù)這一耦合系統(tǒng)進(jìn)行整體動(dòng)力研究,因此理論分析和數(shù)值仿真模擬依然是當(dāng)前主要研究手段。尤其對(duì)于強(qiáng)震區(qū)高土石壩,應(yīng)建立庫(kù)水-高土石壩-地基-孔隙水宏細(xì)觀水土動(dòng)力耦合系統(tǒng),對(duì)其進(jìn)行整體分析、數(shù)學(xué)建模和耦合求解,著重研究水的動(dòng)態(tài)演變過(guò)程及與土骨架之間動(dòng)態(tài)作用關(guān)系,即從庫(kù)水運(yùn)動(dòng)和入滲到超孔隙水的形成,再到超孔隙水的擴(kuò)散和耗散(孔隙水的重分布和流入水庫(kù)),更精確地考察孔隙水可能造成的土體失穩(wěn)、滲透破壞、液化、心墻動(dòng)水劈裂等震害問(wèn)題。

    隨著有限元等數(shù)值求解技術(shù)的飛速發(fā)展,將高土石壩庫(kù)水宏觀流固耦合與土骨架孔隙水細(xì)觀流固耦合完整意義上結(jié)合的壩水動(dòng)力流固耦合成為可能。已有研究人員對(duì)一些實(shí)際土石壩工程開(kāi)展相關(guān)動(dòng)力反應(yīng)分析,如Wang等[50]利用大型商業(yè)有限元軟件ADINA,分別建立了壩體與庫(kù)水網(wǎng)格,對(duì)于流體域分別利用亞音速勢(shì)流體模型與N-S流體模型進(jìn)行模擬,對(duì)于壩體采用基于廣義Biot動(dòng)力固結(jié)的多孔介質(zhì)理論,分析了Sanfernando壩在地震作用下壩體、庫(kù)水以及孔隙水系統(tǒng)動(dòng)力耦合的全過(guò)程,初步證明了該理論應(yīng)用于土石壩的可行性與合理性。Wang等[51]對(duì)新奧爾良17街運(yùn)河大堤的破壞進(jìn)行了分析,同時(shí)采用有效應(yīng)力法和總應(yīng)力法進(jìn)行了對(duì)比,重點(diǎn)考慮了堤壩系統(tǒng)中堤壩填充物和防洪墻所形成的空隙對(duì)堤防性能的影響。牛志偉[52]在研究高壩-庫(kù)水-淤沙系統(tǒng)的動(dòng)力相互作用中,將廣義Biot動(dòng)力固結(jié)理論應(yīng)用于庫(kù)底淤沙,采用P-Z彈塑性本構(gòu)模型,對(duì)庫(kù)水采用簡(jiǎn)化歐拉方程,分析了整個(gè)系統(tǒng)的壩水動(dòng)力響應(yīng)以及淤沙層的液化。李蔚[53]在此基礎(chǔ)上,將該壩水動(dòng)力耦合理論進(jìn)一步應(yīng)用到高土石壩的地震動(dòng)力反應(yīng)分析中,對(duì)雙江口心墻堆石壩進(jìn)行了三維有限元?jiǎng)恿Ψ治?。卞鋒[54]分別采用等價(jià)黏彈性模型與基于廣義Biot動(dòng)力固結(jié)理論的P-Z彈塑性模型對(duì)300 m級(jí)的其宗心墻壩進(jìn)行了動(dòng)力分析,實(shí)現(xiàn)了這一新理論在高土石壩動(dòng)力分析中的應(yīng)用。上述研究為強(qiáng)震作用下庫(kù)水-高土石壩-地基-孔隙水宏細(xì)觀水土動(dòng)力耦合系統(tǒng)研究提供了基礎(chǔ),通過(guò)進(jìn)一步的完善研究,能為科學(xué)合理地建立高土石壩地震反應(yīng)安全評(píng)價(jià)標(biāo)準(zhǔn)及抗震加固措施提供理論依據(jù)。

    4 尚待研究的問(wèn)題

    高土石壩與水體(庫(kù)水和孔隙水)的動(dòng)力流固耦合分析目前尚處于理論發(fā)展階段,還有許多問(wèn)題有待細(xì)致深入研究,如:

    a.對(duì)強(qiáng)震下的高土石壩和大水庫(kù),庫(kù)水運(yùn)動(dòng)需建立更為精細(xì)的控制方程,以反映與壩體運(yùn)動(dòng)之間的雙向動(dòng)力作用。

    b.雖然對(duì)水體采用ALE迎風(fēng)有限元法和ALE分步有限元法已廣泛應(yīng)用于許多流固耦合領(lǐng)域,但將其引入到土石壩庫(kù)水運(yùn)動(dòng)分析中還有待進(jìn)一步的應(yīng)用和驗(yàn)證。

    c.廣義Biot動(dòng)力固結(jié)理論已在海床動(dòng)力特性和混凝土壩壩前庫(kù)底淤沙動(dòng)力特性等方面得到應(yīng)用,但在土石壩地震反應(yīng)分析中應(yīng)用還很少,其中多孔介質(zhì)動(dòng)力本構(gòu)模型以及合理的動(dòng)參數(shù)對(duì)壩水動(dòng)力耦合系統(tǒng)影響很大,還需不斷深入研究。

    d.混合物多孔介質(zhì)模型作為一種較為先進(jìn)的多孔介質(zhì)理論已嶄露頭角,目前多限于靜力滲流固結(jié)等領(lǐng)域,在高土石壩動(dòng)力流固耦合分析中的應(yīng)用還有待深入研究。

    e.土石壩動(dòng)力反應(yīng)分析中壩水宏觀流固耦合與多孔介質(zhì)水土細(xì)觀流固耦合雖已實(shí)現(xiàn)結(jié)合,但庫(kù)水與孔隙水在地震中的交換過(guò)程,即庫(kù)水的動(dòng)力入滲和孔隙水的消散回流至水庫(kù),并未得到全面考慮,也值得進(jìn)一步探索研究。

    f.應(yīng)大力加強(qiáng)高土石壩宏細(xì)觀壩水動(dòng)力耦合效應(yīng)的監(jiān)測(cè)和震后分析等工程應(yīng)用研究,從實(shí)際工程中來(lái)驗(yàn)證和評(píng)價(jià)各種研究理論和分析方法的可行性、重要性和準(zhǔn)確性,以達(dá)到理論與實(shí)踐相互促進(jìn)的目的。

    [1]WESTERGAARD H M.Water pressures on dams during earthquakes[J].Trans,ASCE,1933,98:418-433.

    [2]張振國(guó).考慮壩水相互作用的鋼筋混凝土面板堆石壩三維非線性有限元分析[D].南京:河海大學(xué),1987.

    [3]遲世春,顧淦臣.面板堆石壩壩水系統(tǒng)自振特性研究[J].河海大學(xué)學(xué)報(bào),1995,23(6):104-107.(CHI Shichun,GU Ganchen.Faced rockfilldam vibration characteristics of water systems[J].Journal of Hohai University,1995,23(6):104-107.(in Chinese))

    [4]遲世春,林皋.混凝土面板堆石壩與庫(kù)水動(dòng)力相互作用研究[J].大連理工大學(xué)學(xué)報(bào),1998,38(6):718-723.(CHI Shichun,LIN Gao.Research on hydrodynamic interaction of concrete face rockfill dams and water system[J].Journal of Dalian University of Technology,1998,38(6):718-723.(in Chinese))

    [5]遲世春,林皋.不同動(dòng)水壓力質(zhì)量陣對(duì)面板堆石壩動(dòng)力特性的影響[J].水電站設(shè)計(jì),1999,15(1):46-52.(CHI Shichun,LIN Gao.Influence of different hydrodynamic pressure mass matrix on dynamic characters of concrete facing rockfill dams[J].Design of Hydroelectric Power Station,1999,15(1):46-52.(in Chinese))

    [6]鄧海峰.高面板堆石壩地震三維動(dòng)力反應(yīng)分析[D].宜昌:三峽大學(xué),2010.

    [7]JOHN D A.Computationalfluid dynamics[M].Cambridge:Cambridge University Press,2002.

    [8]HIRSCH C.Numerical computation of internal and external flows[J].Chrchester:John Wiley,1988.

    [9]LASKARIS T E.Finite-elementanalysis ofthreedimensional potential flow in turbomachines[J].AIAA J,1978(16):717.

    [10]OLSIN L G,BATHE K J.Analysis of fluid-structure interactions:a direct symmetric coupled formulation based on the fluid velocity potential[J].Computerand Structure,1985,21(1/2):21-32.

    [11] BATCHELOR G K.An introduction to fluid dynamics[M].Cambridge:Cambridge University Press,1967.

    [12] SUSSMAN T,SUNDQVIST J.Fluid-structure interaction analysis with a subsonic potential-based fluid formulation[J].Computer and Structure,2003,81:949-962.

    [13]王偉華,張燎軍.基于Adina的重力壩地震響應(yīng)分析[J].水電能源科學(xué),2008,26(1):97-99(WANG Weihua,ZHANG Liaojun.ApplicationofADINA to analysis of earthquake response of fravity dam[J].Water Conservancy Science and Technology and Economy,2008,26(1):97-99.(in Chinese))

    [14]劉學(xué)炎.潰壩水流數(shù)值模擬[D].武漢:武漢理工大學(xué),2009.

    [15]劉金云,陳建云.勢(shì)流體與N-S流體的應(yīng)用及比較研究[J].人民長(zhǎng)江,2009,40(20):27-31.(LIU Jinyun,CHEN Jianyun.The comparison between potential fluid and N-S fluid[J].Yangtze River,2009,40(20):27-31.(in Chinese))

    [16]RAMASWAMY B,KWAAHARA M.Arbitrary Lagrangian-Eulerian finite element method for unsteady,convective,incompressible viscousfreesurface fluid flow[J].International Journal for Numerical Methods in Fluids,1987(7):1053-1075.

    [17] SOULI M,ZELESIO J P.Arbitrary Lagrangian-Eluerian and free surface methods in fluid mechanics[J].Computer Methods in Applied Mechanics and Engineering,2001,191:451-466.

    [18] TAKASE S,KASHIYAMA K,TANAKA S.Space-time SUPG finite element computation of shallow-water flows with moving shorelines[J].Computational Mechanics,2011,48(3):293-306.

    [19]岳寶增,劉延柱,王照林.三維液體大幅晃動(dòng)及其抑制的數(shù)值模擬[J].上海交通大學(xué)學(xué)報(bào),2000,34(8):1036-1039.(YUE Baozen,LIU Yanzhu,WANG Zhaolin.Numerical simulation and suppression of three dimension large amplitude liquid sloshing[J].Journal of Shanghai Jiaotong University:2000,34(8):1036-1039.(in Chinese))

    [20]岳寶增,劉延柱,王照林.三維液體非線性晃動(dòng)動(dòng)力學(xué)特性的數(shù)值模擬[J].應(yīng)用力學(xué)學(xué)報(bào),2001,18(1):110-115.(YUE Baozen,LIU Yanzhu,WANG Zhaolin.ALE finite element method for three-dimensional large amplitude liquid sloshing using fractional step method[J].Chinese Journal of Applied Mechanics,2001,18(1):110-115.(in Chinese))

    [21]華蕾娜.ALE分步有限元法研究及其在自由表面水波問(wèn)題中的應(yīng)用[D].天津:天津大學(xué),2005.

    [22]陳文元,趙雷.壩體考慮流固耦合的動(dòng)力特性分析[J].四川建筑科學(xué)研究,2013,39(1):126-129.(CHENG Wenyuan,ZHAO Lei.Analysis dynamic characteristics of the dam considering fluid-structure coupling[J].Sichuan Building Science,2013,39(1):126-129.(in Chinese))

    [23]王國(guó)輝,柴軍瑞,杜成偉.壩庫(kù)系統(tǒng)流固耦合數(shù)值分析方法簡(jiǎn)述[J].水利水電科技進(jìn)展,2009,29(4):89-94.(WANG Guohui,CHAI Junrui,DU Chengwei.Review of numerical analysis method for fluid-solid interaction of dam-reservoir system[J].Advances in Science and Technology of Water Resources,2009,29(4):89-94.(in Chinese))

    [24]徐志英,沈珠江.地震液化的有效應(yīng)力二維動(dòng)力分析方法[J].華東水利學(xué)院學(xué)報(bào),1981,9(3):1-14.(XU Zhiying,SHENG Zhujiang.Twodimensionalanalysis methods of effective stress liquefaction[J].Journal of East China Institute of Water Conservancy,1981,9(3)1-14.(in Chinese))

    [25]周健,徐志英.土(尾礦)壩的三維有效應(yīng)力的動(dòng)力反應(yīng)分析[J].地震工程和工程震動(dòng),1984,4(3):60-70.(ZHOU Jian,XU Zhiyi.Three dimensional dynamic analysis of effective stress in soil(tailings)dams[J].Earthquake Engineering and Engineering Vibration,1984,4(3):60-70.(in Chinese))

    [26]BIOT M A.Theory of elasticity and consolidation for a porous anisotropic solid[J].Appl Phy,1955,26:182-185.

    [27]BIOT M A.Theory of propagation of elastic waves in a fluid-saturated porous solid[J].Acoust Soc Am,1956,28:168-191.

    [28]BIOT M A.Mechanicsofdeformation and acoustic propagation in porous media[J].Appl Phy,1962,33(4):1482-1498.

    [29]BIOT M A,Willis P G.The elastic coefficients of the theory consolidation[J].Appl Phy,1957,24:594-601.

    [30]GHABOUSSI J,WILSON E L.Variational formulation of dynamics of fluid-saturated porous elastic solids[J].Engng Mech Div,ASCE,1972,98(4):947-963.

    [31]GHABOUSSI J,WILSON E L.Seismic analysis of earth dam-reservoir system[J].Soil Mech Found Div,ASCE,1973,99(10):849-862.

    [32]門(mén)福錄.波在飽和孔隙彈性介質(zhì)中的傳播[J].地球物理學(xué)報(bào),1965,14(2):107-114.(MEN Fulu.Waves propagation in saturated porous elastic medium[J].Geophysics,1965,14(2):107-114.(in Chinese))

    [33]盛虞,盧盛松,姜樸.土工建筑物動(dòng)力固結(jié)的耦合振動(dòng)分析[J].水利學(xué)報(bào),1989(12):31-42.(SHEN Yu,LU Shengsong,JIANG Pu.Earth structure coupled vibration analysisofdynamic consolidation[J].Journalof Hydraulic Engineering,1989(12):31-42.(in Chinese))

    [34]林本海.砂土-碎石樁復(fù)合地基的液化檢驗(yàn)理論和數(shù)值分析方法的研究[D].西安:西安理工大學(xué),1997.

    [35]ZIENKIEWICZ O C.Basic formulation of static and dynamic behaviour of soil and other porous media[J].Numerical Methods in Geomechanics,1982,92:39-55.

    [36] ZIENKIEWICZ O C,SHIOMI T.Dynamic behavior of saturated porous media:the generalized Biot formulation and its numerical solution[J].Int Numer Anal Methods and Geometry,1984(8):71-96.

    [37]SIMON B R,ZIENKIEWICZ O C,PAUL D K.An analytical solution for the transient response of saturated porous elastic solids[J].Int Numer Anal Methods and Geometry,1984(8):381-398.

    [38]SIMON B R,WU J S,ZIENKIEWICZ O C,PAUL D K.Evaluation of u-w and u-π finite element methods for the dynamic response of saturated porous media using onedimensional models[J].Int Methods and Geometry,1986(10):461-482.

    [39]SIMON B R,WU J S,ZIENKIEWICZ O C.Evaluation of higherorder,mixed and Hermitian finite element procedures for the dynamic response of saturated porous media using one-dimensional models[J].Int Numer Anal Methods and Geometry,1986(10):483-499.

    [40]李宏儒.土體動(dòng)力反應(yīng)分析方法的分析與改進(jìn)[D].西安:西安理工大學(xué),2005.

    [41]劉凱欣,劉穎.橫觀各向同性含液飽和多孔介質(zhì)中瑞利波的特性分析[J].力學(xué)學(xué)報(bào),2003,35(1):100-104.(LIU Kaixin,LIU Ying.Characteristic of Rayleigh waves in transversely isotropic fluid saturated porous media[J].Acta Mechanica Sinica,2003,35(1):100-104.(in Chinese))

    [42]TRUESDELL C.Exact theory of stress and strain in rods and shells[J].Archive for Rational Mechanics and Analysis,1982,1(1):295-323.

    [43]PREVOST JH.Non-lineartransientphenomenain saturated porous media[J].Comput Methods Appl Mech Engrg,1982,30:3-18.

    [44]PREVOST J H.Wave propagation in fluid-saturated porous media:an efficient finite element procedure[J].Int Soil Dyn Earthquake Engrg,1985,4(4):183-202.

    [45] YIAGOS A N,PREVOST J H.Two-phase elasto-plastic seismic response of earth dams:applications[J].Soil Dynamics and Earthquake Engineering,1991,10(7):371-381.

    [46]嚴(yán)波,張汝清.黏性流體兩相多孔介質(zhì)非線性動(dòng)力問(wèn)題的罰有限元法[J].應(yīng)用力學(xué)和數(shù)學(xué),2000,21(12):1247-1253.(YAN Bo,ZHANG Ruqing.Penalty finite element method for nonlinear dynamic response of viscous fluid-saturated biphasic porous media[J].Applied Mathematics and Mechanics,2000,21(12):1247-1253.(in Chinese))

    [47]秦小軍,陳少林,曾心傳.二維飽和孔隙介質(zhì)的三場(chǎng)有限元方法[J].地殼形變與地震,1999,19(2):60-69.(QING Xiaojun,CHEN Shaolin,ZENG Xinchuan.A three-field finite element procedure for two dimensional fluid-saturated porous media[J].Crustal Deformation and Earthquake,1999,19(2):60-69.(in Chinese))

    [48]BOWEN R M.Plane progressive waves in a heat conducting fluid saturated porous material with relaxing porosity[J].Acta Mechanica,1983,46:189-200.

    [49]de BOER R.多孔介質(zhì)理論發(fā)展史上的重要成果[M].劉占芳,嚴(yán)波,譯.重慶:重慶大學(xué)出版社,1995.

    [50]WANG X,WANG L B.Dynamic analysis of a water-soilpore water coupling system[J].Computers & Structures,2007,85:1020-1031.

    [51]WANG X,CASTAY M.Failure analysis of the breached levee at the 17th Street Canal in New Orleans during Hurricane Katrina[J].Canadian Geotechnical Journal,2012,49(7):812-834.

    [52]牛志偉.庫(kù)底淤沙對(duì)混凝土壩地震響應(yīng)影響分析方法研究[D].南京:河海大學(xué),2008.

    [53]李蔚.考慮壩體-水體-地基相互作用的高土石壩地震動(dòng)力反應(yīng)分析研究[D].南京:河海大學(xué),2010.

    [54]卞鋒.高土石壩地震動(dòng)力響應(yīng)特征彈塑性有限元分析[D].北京:清華大學(xué),2010.

    Recent advances on macro-scale and micro-scale dynamic interaction between high earth-rock dams and water

    CEN Weijun,SUN Hui,CHEN Yanan(College of Water Conservancy and Hydropower Engineering,Hohai University,Nanjing 210098,China)

    The dynamic interaction will happen among high earth-rock dams,reservoir and pore water under a strong earthquake.This paper introduces the main research history,research results and research trends of water-soil coupling system for earth-rock dams from the macro-scale and micro-scale aspects.On the macro-scale dynamic coupling system,the early dam-water interaction model,sophisticated analytical methods of reservoir water and FSI coordinate system description are summarized,and on the micro-scale dynamic coupling system,mainly from the early water-soil decoupled dynamic analysis method,and the water-soil coupling analysis system based on the Boit’s theory and the water-soil mixture theory,ect.Combining the above two research aspects of soil-water coupling system,a macro-scale and micro-scale coupling system of reservoir-dam-pore water is proposed to study the dynamic behavior of high earth-rock dams comprehensively.Finally,the main problems on the macro-scale and micro-scale water-soil coupling theory and further research topics are put forward.

    high earth-rock dams;reservoir water;pore water;macro-scale dam and water interaction;micro-scale soil skeleton and pore water interaction;dynamic analysis

    TV641.1

    A

    1006-7647(2013)06-0010-07

    10.3880/j.issn.1006-7647.2013.06.002

    國(guó)家自然科學(xué)基金(51009055);中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)(2009B07514)

    岑威鈞(1977—),男,浙江慈溪人,副教授,博士,主要從事土石壩工程靜動(dòng)力分析及水工滲流研究。E-mail:hhucwj@163.com

    2013-04-12 編輯:熊水斌)

    猜你喜歡
    庫(kù)水細(xì)觀石壩
    三峽庫(kù)區(qū)舊縣坪滑坡變形機(jī)理及穩(wěn)定性
    庫(kù)水升降對(duì)新疆阿勒泰地區(qū)某土石壩滲流影響的研究
    地下水(2021年3期)2021-07-05 05:25:10
    庫(kù)水可壓縮性對(duì)重力壩動(dòng)力特性和地震響應(yīng)的影響
    基于細(xì)觀結(jié)構(gòu)的原狀黃土動(dòng)彈性模量和阻尼比試驗(yàn)研究
    地震研究(2021年1期)2021-04-13 01:05:24
    庫(kù)水滲流作用下邊坡穩(wěn)定性分析
    歡迎訂閱《碾壓式土石壩設(shè)計(jì)》
    新型加筋土技術(shù)在土石壩除險(xiǎn)加固中的應(yīng)用
    基于四叉樹(shù)網(wǎng)格加密技術(shù)的混凝土細(xì)觀模型
    PBX炸藥的抗壓強(qiáng)度及抗拉強(qiáng)度細(xì)觀尺度的數(shù)值計(jì)算
    開(kāi)裂混凝土中水分傳輸過(guò)程的細(xì)觀模型
    一个人免费在线观看电影| 一级毛片电影观看 | 中文亚洲av片在线观看爽| 精品人妻视频免费看| 床上黄色一级片| videossex国产| 亚洲av成人av| 热99re8久久精品国产| 俄罗斯特黄特色一大片| 免费av毛片视频| 欧美bdsm另类| 欧美bdsm另类| 香蕉av资源在线| 亚洲av.av天堂| 国产精品久久久久久精品电影| 丰满的人妻完整版| 高清日韩中文字幕在线| 久久久久久久久久黄片| 国产私拍福利视频在线观看| 精品午夜福利在线看| 亚洲精品色激情综合| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲精品成人久久久久久| 91麻豆精品激情在线观看国产| 一进一出好大好爽视频| 亚洲精品色激情综合| 少妇裸体淫交视频免费看高清| 在现免费观看毛片| 成年av动漫网址| 日韩在线高清观看一区二区三区| 精品少妇黑人巨大在线播放 | 免费看美女性在线毛片视频| 欧美色视频一区免费| 国产精品一二三区在线看| 亚洲精品乱码久久久v下载方式| 国产在视频线在精品| 日韩欧美国产在线观看| 精品久久久久久成人av| 亚洲在线观看片| 成年av动漫网址| 精品国内亚洲2022精品成人| 成人精品一区二区免费| 亚洲经典国产精华液单| 日韩av不卡免费在线播放| 在线观看66精品国产| 国产精品亚洲美女久久久| 午夜精品一区二区三区免费看| 亚洲av中文字字幕乱码综合| 国产男靠女视频免费网站| 综合色av麻豆| 女人被狂操c到高潮| 日日摸夜夜添夜夜添av毛片| 成人av一区二区三区在线看| 亚洲精品日韩在线中文字幕 | 欧美+日韩+精品| 波多野结衣巨乳人妻| 又爽又黄a免费视频| h日本视频在线播放| 91av网一区二区| 99国产极品粉嫩在线观看| 一级黄色大片毛片| 插逼视频在线观看| 一个人看的www免费观看视频| 久久久久国产网址| 日本-黄色视频高清免费观看| 国产精品综合久久久久久久免费| 亚洲av美国av| 国产高潮美女av| 天天躁日日操中文字幕| 床上黄色一级片| 搡老妇女老女人老熟妇| 日本a在线网址| 午夜免费激情av| 国产成人福利小说| 噜噜噜噜噜久久久久久91| 淫秽高清视频在线观看| 亚洲av第一区精品v没综合| 成人无遮挡网站| 在线看三级毛片| 亚洲人成网站在线播| 99久久中文字幕三级久久日本| 久久精品国产自在天天线| 国产精品国产三级国产av玫瑰| 免费观看人在逋| 禁无遮挡网站| 色尼玛亚洲综合影院| 久久久欧美国产精品| 看黄色毛片网站| 欧美日本视频| 日本成人三级电影网站| 中文字幕av成人在线电影| 一a级毛片在线观看| 天堂动漫精品| 又爽又黄a免费视频| 久久久精品大字幕| 欧美日韩一区二区视频在线观看视频在线 | 99久久中文字幕三级久久日本| 免费搜索国产男女视频| 国产色爽女视频免费观看| 欧美三级亚洲精品| 精品久久国产蜜桃| 日韩 亚洲 欧美在线| 天堂网av新在线| 欧美三级亚洲精品| 国产精品女同一区二区软件| 在线播放国产精品三级| 人妻夜夜爽99麻豆av| 有码 亚洲区| 午夜激情福利司机影院| 亚洲最大成人中文| 免费电影在线观看免费观看| 亚洲国产色片| 内射极品少妇av片p| 一区福利在线观看| 日韩欧美精品免费久久| 亚洲av成人精品一区久久| 国产午夜福利久久久久久| 99在线人妻在线中文字幕| 熟女人妻精品中文字幕| 尾随美女入室| eeuss影院久久| 嫩草影院入口| 国产欧美日韩精品一区二区| 亚洲成人精品中文字幕电影| 久久午夜福利片| 国产 一区 欧美 日韩| 三级国产精品欧美在线观看| 听说在线观看完整版免费高清| 女人被狂操c到高潮| 亚洲经典国产精华液单| 欧美性猛交黑人性爽| 非洲黑人性xxxx精品又粗又长| 特级一级黄色大片| 精品乱码久久久久久99久播| 性欧美人与动物交配| 老师上课跳d突然被开到最大视频| 12—13女人毛片做爰片一| 淫秽高清视频在线观看| 国产乱人偷精品视频| 永久网站在线| 激情 狠狠 欧美| 亚洲欧美精品综合久久99| .国产精品久久| АⅤ资源中文在线天堂| 老司机午夜福利在线观看视频| 三级男女做爰猛烈吃奶摸视频| 久久久久久久久久黄片| 国产精品伦人一区二区| 日本-黄色视频高清免费观看| 欧美一区二区国产精品久久精品| 成年av动漫网址| 不卡一级毛片| 欧美性猛交╳xxx乱大交人| 亚洲欧美精品综合久久99| 国产亚洲精品久久久久久毛片| 欧美精品国产亚洲| 成年女人毛片免费观看观看9| 一级黄色大片毛片| 最近视频中文字幕2019在线8| 91精品国产九色| 国产亚洲精品久久久久久毛片| 国国产精品蜜臀av免费| 内地一区二区视频在线| 寂寞人妻少妇视频99o| 黄色视频,在线免费观看| 欧美成人免费av一区二区三区| 成人亚洲欧美一区二区av| 国产欧美日韩精品亚洲av| 午夜激情福利司机影院| 亚洲国产高清在线一区二区三| 国产黄a三级三级三级人| 国产亚洲欧美98| av天堂中文字幕网| 淫秽高清视频在线观看| 久久久a久久爽久久v久久| 日本免费一区二区三区高清不卡| 男女视频在线观看网站免费| 国产精品野战在线观看| 婷婷六月久久综合丁香| 国产中年淑女户外野战色| 少妇高潮的动态图| a级毛片免费高清观看在线播放| 欧美中文日本在线观看视频| 亚洲成人久久性| 欧美日韩乱码在线| 日日摸夜夜添夜夜添小说| 成人三级黄色视频| 国产成人aa在线观看| 久久鲁丝午夜福利片| 搡老岳熟女国产| 欧美日本亚洲视频在线播放| 直男gayav资源| 一卡2卡三卡四卡精品乱码亚洲| 小蜜桃在线观看免费完整版高清| 女同久久另类99精品国产91| 日本色播在线视频| 91久久精品国产一区二区三区| 乱系列少妇在线播放| 夜夜夜夜夜久久久久| 亚洲性夜色夜夜综合| 国产片特级美女逼逼视频| 国产片特级美女逼逼视频| 老女人水多毛片| 午夜视频国产福利| 人妻夜夜爽99麻豆av| 在线免费观看的www视频| 男插女下体视频免费在线播放| 久久人妻av系列| 亚洲中文字幕一区二区三区有码在线看| 亚洲av免费在线观看| 国产精品嫩草影院av在线观看| 丝袜喷水一区| 午夜福利在线观看免费完整高清在 | 色视频www国产| 亚洲欧美成人精品一区二区| 亚洲精品在线观看二区| 精品久久国产蜜桃| 非洲黑人性xxxx精品又粗又长| 国产精品人妻久久久久久| 国产av在哪里看| 俄罗斯特黄特色一大片| 在线看三级毛片| 搞女人的毛片| 黄片wwwwww| 九九在线视频观看精品| 婷婷亚洲欧美| 日韩精品有码人妻一区| 观看美女的网站| 噜噜噜噜噜久久久久久91| 最近的中文字幕免费完整| 国产综合懂色| 三级男女做爰猛烈吃奶摸视频| 午夜影院日韩av| 欧美激情在线99| 校园春色视频在线观看| 全区人妻精品视频| 亚洲av成人av| 91精品国产九色| 亚洲国产精品合色在线| 婷婷精品国产亚洲av在线| 免费看美女性在线毛片视频| 99久国产av精品| 国产老妇女一区| 色综合站精品国产| 成年女人毛片免费观看观看9| 国产老妇女一区| 国产欧美日韩精品亚洲av| 麻豆精品久久久久久蜜桃| 久久久久国产精品人妻aⅴ院| 97超视频在线观看视频| 三级毛片av免费| 精品少妇黑人巨大在线播放 | 午夜免费男女啪啪视频观看 | 午夜福利成人在线免费观看| .国产精品久久| 九色成人免费人妻av| 亚洲欧美中文字幕日韩二区| 麻豆一二三区av精品| 此物有八面人人有两片| 欧美一区二区亚洲| 日本三级黄在线观看| 成人亚洲欧美一区二区av| 国产欧美日韩精品亚洲av| 亚洲国产色片| 人人妻人人看人人澡| 午夜福利在线在线| 久久久久久国产a免费观看| 三级毛片av免费| 日日摸夜夜添夜夜添小说| 日本五十路高清| a级毛片免费高清观看在线播放| 午夜激情福利司机影院| 日日摸夜夜添夜夜添小说| 亚洲第一区二区三区不卡| av免费在线看不卡| 18禁裸乳无遮挡免费网站照片| 99热这里只有是精品在线观看| 亚洲成av人片在线播放无| 亚洲av中文av极速乱| 一级毛片电影观看 | 中出人妻视频一区二区| 校园春色视频在线观看| 成年版毛片免费区| 国内精品美女久久久久久| 国产一区亚洲一区在线观看| 三级毛片av免费| 啦啦啦观看免费观看视频高清| 成人高潮视频无遮挡免费网站| 亚洲成人中文字幕在线播放| 亚洲国产欧美人成| 国产亚洲精品久久久com| 中文字幕久久专区| 一级黄色大片毛片| 丰满人妻一区二区三区视频av| 亚洲av一区综合| 亚洲欧美日韩高清专用| 日韩在线高清观看一区二区三区| ponron亚洲| 天美传媒精品一区二区| a级毛色黄片| 日韩在线高清观看一区二区三区| 人妻夜夜爽99麻豆av| 最近视频中文字幕2019在线8| 免费看日本二区| 免费一级毛片在线播放高清视频| 91午夜精品亚洲一区二区三区| 免费av毛片视频| 午夜精品在线福利| 欧美日韩国产亚洲二区| 午夜福利18| 国产黄色小视频在线观看| 乱码一卡2卡4卡精品| 午夜福利在线观看免费完整高清在 | 精品久久久久久久久久久久久| 亚洲中文字幕日韩| 久久这里只有精品中国| 国产精品嫩草影院av在线观看| 免费人成在线观看视频色| 亚洲三级黄色毛片| 亚洲国产精品sss在线观看| а√天堂www在线а√下载| 日本三级黄在线观看| 国产真实伦视频高清在线观看| 日本黄色片子视频| 91久久精品电影网| av在线老鸭窝| av在线亚洲专区| 一卡2卡三卡四卡精品乱码亚洲| av专区在线播放| 国产伦一二天堂av在线观看| 久久久国产成人精品二区| 欧美成人精品欧美一级黄| 亚洲久久久久久中文字幕| 大又大粗又爽又黄少妇毛片口| 中文字幕av成人在线电影| 日产精品乱码卡一卡2卡三| 国产精品,欧美在线| 禁无遮挡网站| 99热这里只有是精品50| av.在线天堂| 在线观看一区二区三区| 综合色丁香网| 午夜精品国产一区二区电影 | 又爽又黄无遮挡网站| 日本熟妇午夜| 久久久精品欧美日韩精品| 看片在线看免费视频| 亚洲自偷自拍三级| 精品久久久久久久久久免费视频| 国产精品,欧美在线| 精品久久久久久久久久久久久| 黄片wwwwww| 最近的中文字幕免费完整| 三级国产精品欧美在线观看| 一本一本综合久久| 12—13女人毛片做爰片一| 在现免费观看毛片| 欧美区成人在线视频| 久久韩国三级中文字幕| 成年版毛片免费区| 夜夜爽天天搞| 亚洲精品一卡2卡三卡4卡5卡| 狠狠狠狠99中文字幕| 国产午夜精品论理片| a级毛片免费高清观看在线播放| 成人综合一区亚洲| АⅤ资源中文在线天堂| 天堂网av新在线| 亚洲一级一片aⅴ在线观看| 高清毛片免费看| 久久国产乱子免费精品| 如何舔出高潮| 男女做爰动态图高潮gif福利片| 老司机福利观看| 国产 一区精品| 日本黄色片子视频| 亚洲精品成人久久久久久| 成年女人毛片免费观看观看9| 午夜影院日韩av| 麻豆精品久久久久久蜜桃| 男人舔奶头视频| 精品一区二区三区视频在线| 一进一出抽搐gif免费好疼| 老司机影院成人| 亚洲中文字幕日韩| 午夜精品在线福利| 国产女主播在线喷水免费视频网站 | 一级黄片播放器| 在线观看av片永久免费下载| 午夜激情欧美在线| 亚洲真实伦在线观看| 一级黄片播放器| av中文乱码字幕在线| 久久久国产成人免费| 免费在线观看影片大全网站| 国产真实乱freesex| 12—13女人毛片做爰片一| 最近在线观看免费完整版| 精品熟女少妇av免费看| 日韩亚洲欧美综合| 嫩草影视91久久| 久久精品夜色国产| 国产男靠女视频免费网站| 在线看三级毛片| 久久久午夜欧美精品| 国产中年淑女户外野战色| 嫩草影视91久久| aaaaa片日本免费| 日本五十路高清| 亚洲成a人片在线一区二区| 又粗又爽又猛毛片免费看| 最近的中文字幕免费完整| 人人妻,人人澡人人爽秒播| 国产极品精品免费视频能看的| 69av精品久久久久久| 婷婷精品国产亚洲av| 六月丁香七月| 亚洲人成网站在线播放欧美日韩| 99久国产av精品国产电影| 国产成人精品久久久久久| 一本一本综合久久| 久久欧美精品欧美久久欧美| 欧美中文日本在线观看视频| 搡老岳熟女国产| 精品日产1卡2卡| av国产免费在线观看| 成熟少妇高潮喷水视频| 日韩制服骚丝袜av| 男女那种视频在线观看| 国产成人福利小说| 国产亚洲精品综合一区在线观看| 久久精品综合一区二区三区| 久久久久久久久久成人| 99国产精品一区二区蜜桃av| av天堂中文字幕网| 日本欧美国产在线视频| 天天躁夜夜躁狠狠久久av| 三级毛片av免费| 99国产极品粉嫩在线观看| 国产亚洲欧美98| 少妇被粗大猛烈的视频| 国产av一区在线观看免费| 久久中文看片网| 午夜精品一区二区三区免费看| 欧美日韩国产亚洲二区| 欧美性猛交╳xxx乱大交人| 欧美+日韩+精品| 亚洲欧美日韩卡通动漫| 一个人看的www免费观看视频| 久久午夜福利片| 成人特级黄色片久久久久久久| 男人和女人高潮做爰伦理| 1000部很黄的大片| 国产淫片久久久久久久久| 欧美又色又爽又黄视频| 91精品国产九色| 热99在线观看视频| 久久久久久久亚洲中文字幕| 色哟哟·www| 国产精品不卡视频一区二区| 女人十人毛片免费观看3o分钟| 亚洲色图av天堂| 亚洲aⅴ乱码一区二区在线播放| .国产精品久久| 在线观看午夜福利视频| 色哟哟哟哟哟哟| 狂野欧美白嫩少妇大欣赏| 免费看a级黄色片| 特大巨黑吊av在线直播| 香蕉av资源在线| 中文字幕久久专区| 乱人视频在线观看| 国产精品一区二区三区四区免费观看 | 国产乱人偷精品视频| 色av中文字幕| 男女下面进入的视频免费午夜| 成熟少妇高潮喷水视频| 深夜精品福利| 午夜日韩欧美国产| 成人午夜高清在线视频| 亚洲中文字幕日韩| 99热6这里只有精品| 成人三级黄色视频| 日韩欧美国产在线观看| 精品久久国产蜜桃| 欧美+日韩+精品| 欧美xxxx黑人xx丫x性爽| 国产成人精品久久久久久| 99热全是精品| 男女下面进入的视频免费午夜| 久久久久久久亚洲中文字幕| 最近视频中文字幕2019在线8| 久久鲁丝午夜福利片| 久久久久国内视频| 中文字幕av成人在线电影| 国产v大片淫在线免费观看| 菩萨蛮人人尽说江南好唐韦庄 | 国产淫片久久久久久久久| 亚洲乱码一区二区免费版| 亚洲久久久久久中文字幕| 波多野结衣高清作品| 国产精品人妻久久久久久| 亚洲七黄色美女视频| 99视频精品全部免费 在线| 亚洲国产色片| 日韩欧美精品v在线| 国内久久婷婷六月综合欲色啪| 免费高清视频大片| 国产精品亚洲美女久久久| 成人无遮挡网站| ponron亚洲| 老司机午夜福利在线观看视频| 国产精品久久久久久精品电影| 精品国内亚洲2022精品成人| 校园人妻丝袜中文字幕| 插逼视频在线观看| 中文字幕精品亚洲无线码一区| 五月伊人婷婷丁香| 黄色欧美视频在线观看| 免费大片18禁| 你懂的网址亚洲精品在线观看 | 麻豆久久精品国产亚洲av| 听说在线观看完整版免费高清| 老司机午夜福利在线观看视频| 亚洲专区国产一区二区| 99九九线精品视频在线观看视频| 啦啦啦啦在线视频资源| 欧美不卡视频在线免费观看| 久久久久久国产a免费观看| 成人高潮视频无遮挡免费网站| 一本久久中文字幕| 欧美日本亚洲视频在线播放| 久久久久久久久久久丰满| 精品久久久久久久久av| 久久鲁丝午夜福利片| 日韩亚洲欧美综合| 18禁黄网站禁片免费观看直播| 成人美女网站在线观看视频| 综合色丁香网| 三级国产精品欧美在线观看| 极品教师在线视频| 亚洲熟妇中文字幕五十中出| 国产精品乱码一区二三区的特点| 99久久精品国产国产毛片| 国产午夜精品久久久久久一区二区三区 | 精品国内亚洲2022精品成人| 嫩草影院新地址| 亚洲中文字幕日韩| 又爽又黄无遮挡网站| 国产成年人精品一区二区| 中国美白少妇内射xxxbb| 一本久久中文字幕| 国产一区二区三区av在线 | 成人av一区二区三区在线看| 99久国产av精品| 国产高清视频在线观看网站| 久久鲁丝午夜福利片| 国产真实乱freesex| 最新中文字幕久久久久| 欧美潮喷喷水| 亚洲aⅴ乱码一区二区在线播放| 欧美一区二区精品小视频在线| 特大巨黑吊av在线直播| 国产黄a三级三级三级人| 亚洲国产精品合色在线| 亚洲国产精品成人久久小说 | 99久久无色码亚洲精品果冻| 97碰自拍视频| 直男gayav资源| 嫩草影院精品99| 精品不卡国产一区二区三区| 免费在线观看成人毛片| 久久久久久久久大av| 成人永久免费在线观看视频| 午夜福利视频1000在线观看| 亚洲欧美日韩无卡精品| 欧美最黄视频在线播放免费| а√天堂www在线а√下载| 国产精品一区二区性色av| 99久国产av精品| 九色成人免费人妻av| 欧美激情久久久久久爽电影| 熟女人妻精品中文字幕| 草草在线视频免费看| 午夜免费男女啪啪视频观看 | 97碰自拍视频| 观看免费一级毛片| 观看美女的网站| 国产高清视频在线播放一区| 久久久久久久久久久丰满| 美女xxoo啪啪120秒动态图| 欧美日韩一区二区视频在线观看视频在线 | 久久6这里有精品| ponron亚洲| 国产精品伦人一区二区| 亚洲欧美成人综合另类久久久 | 亚洲一区高清亚洲精品| 色噜噜av男人的天堂激情| 乱码一卡2卡4卡精品| 亚洲av第一区精品v没综合| 成人av一区二区三区在线看| 亚洲第一电影网av| 欧美性猛交╳xxx乱大交人| 不卡一级毛片| 亚洲国产精品久久男人天堂| or卡值多少钱| 成人特级黄色片久久久久久久| 久久精品国产自在天天线| 国产色爽女视频免费观看| 一区二区三区四区激情视频 | 中文字幕熟女人妻在线| 亚洲国产高清在线一区二区三| 欧美另类亚洲清纯唯美| 日韩强制内射视频| 亚洲电影在线观看av|