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

    多場耦合作用下高混凝土壩水力劈裂研究綜述

    2013-08-15 00:51:44沈振中徐力群
    水利水電科技進(jìn)展 2013年2期
    關(guān)鍵詞:水力滲流裂紋

    甘 磊,沈振中,徐力群

    (1.河海大學(xué)水利水電學(xué)院,江蘇南京 210098; 2.河海大學(xué)水文水資源與水利工程科學(xué)國家重點實驗室,江蘇南京 210098)

    多場耦合作用下高混凝土壩水力劈裂研究綜述

    甘 磊1,2,沈振中1,2,徐力群1,2

    (1.河海大學(xué)水利水電學(xué)院,江蘇南京 210098; 2.河海大學(xué)水文水資源與水利工程科學(xué)國家重點實驗室,江蘇南京 210098)

    闡述了混凝土結(jié)構(gòu)水力劈裂的研究現(xiàn)狀,結(jié)合多場耦合作用下高混凝土壩水力劈裂研究中存在的問題,從水工混凝土結(jié)構(gòu)水力劈裂理論、試驗技術(shù)、數(shù)值模擬和不同耦合條件下裂縫擴(kuò)展研究等方面歸納了國內(nèi)外的研究成果,提出了多場耦合作用下高混凝土壩水力劈裂試驗儀器設(shè)備研制、試驗測試技術(shù)、多場耦合作用下裂縫擴(kuò)展機(jī)理及耦合機(jī)理、考慮結(jié)構(gòu)裂縫任意性的數(shù)值模擬技術(shù)等尚待進(jìn)一步研究的問題。

    高混凝土壩;水力劈裂;多場耦合;裂縫擴(kuò)展;綜述

    水工混凝土內(nèi)部存在著隨機(jī)分布的初始微觀裂紋,當(dāng)混凝土結(jié)構(gòu)受力時,這些缺陷極易產(chǎn)生應(yīng)力集中并迅速起裂、擴(kuò)展、匯聚形成宏觀裂縫,誘發(fā)混凝土結(jié)構(gòu)的宏觀開裂或破壞[5]。同時,混凝土中的裂紋通常不是單一的,而是由多種裂紋共同擴(kuò)展形成。對于水工混凝土水力劈裂問題,由于水壓力和溫度荷載作用下的混凝土裂縫擴(kuò)展試驗難度大,裝置復(fù)雜,導(dǎo)致目前的相關(guān)試驗研究還不多。

    水工混凝土水力劈裂問題的實質(zhì)就是具有高勢能的水對混凝土的破壞,由于混凝土內(nèi)孔隙水的滲流過程及孔隙水壓力的存在,使得對滲流與應(yīng)力耦合作用下的混凝土水力劈裂機(jī)理過程的描述更為復(fù)雜和困難。同時,混凝土的溫度改變將引起混凝土應(yīng)力發(fā)生變化,當(dāng)約束拉應(yīng)變超過混凝土的極限拉伸值時,混凝土結(jié)構(gòu)發(fā)生開裂破壞。溫度變化引起熱應(yīng)力和熱應(yīng)變使得應(yīng)力場和位移場發(fā)生改變,從而引起水工混凝土裂縫特性發(fā)生改變,甚至引起裂縫的進(jìn)一步擴(kuò)展。滲流-應(yīng)力-溫度等多場耦合作用下高混凝土壩裂縫產(chǎn)生、發(fā)展的機(jī)理以及模擬分析方法的研究,具有重要的理論意義和應(yīng)用價值,研究成果可廣泛應(yīng)用于水工混凝土結(jié)構(gòu)的設(shè)計、施工和安全評價。

    1 混凝土裂縫擴(kuò)展理論研究

    材料的破壞過程是兩大力學(xué)難題之一。傳統(tǒng)的力學(xué)強(qiáng)度計算理論以材料力學(xué)和結(jié)構(gòu)力學(xué)為基礎(chǔ),通常假定材料為均勻連續(xù)體,而且認(rèn)為材料各向同性,避開了材料客觀存在的缺陷、裂縫等,無法解釋理論強(qiáng)度遠(yuǎn)高于實際強(qiáng)度的原因,也不能描述裂縫尖端存在應(yīng)力集中嚴(yán)重的現(xiàn)象。

    斷裂力學(xué)主要研究對象為裂紋,通過考慮裂紋尖端的應(yīng)力場和位移場,計算裂紋尖端的應(yīng)力強(qiáng)度因子,從而分析荷載作用下裂紋的擴(kuò)展規(guī)律[6-7]。1913年Inglis采用應(yīng)力集中觀點解釋了由于材料存在缺陷使得材料的實際強(qiáng)度遠(yuǎn)低于理論強(qiáng)度。1921 年,Griffith提出了脆性材料裂紋擴(kuò)展的能量準(zhǔn)則,用彈性體能量平衡的觀點研究了玻璃、陶瓷等脆性材料中的裂紋擴(kuò)展問題。1957年Irwin提出了應(yīng)力場強(qiáng)度觀點及應(yīng)力強(qiáng)度因子斷裂準(zhǔn)則。應(yīng)力強(qiáng)度因子斷裂準(zhǔn)則與裂縫擴(kuò)展的能量準(zhǔn)則構(gòu)成了線彈性斷裂力學(xué)的核心內(nèi)容。彈塑性斷裂研究的發(fā)展起步于20世紀(jì)60年代,英國科學(xué)家認(rèn)為彈塑性斷裂過程聚集于裂縫前方的條狀屈服區(qū),并建立了BCS連續(xù)位錯模型、Cottrell斷裂過程區(qū)模型、Dugdale簡化條狀塑性區(qū)模型、Barenblatt內(nèi)聚斷裂模型等彈塑性斷裂模型。1965年Wells根據(jù)上述模型提出以裂縫尖端張開位移COD(crack opening displacement)作為表征斷裂的參量,建立了裂縫的斷裂準(zhǔn)則。

    從應(yīng)力強(qiáng)度因子出發(fā),1968年美國科學(xué)家提出以J積分和HRR理論作為斷裂力學(xué)參量。1991年美國電力研究院以J積分為斷裂判據(jù),推出了較完整的延性斷裂手冊,標(biāo)志著彈塑性斷裂力學(xué)的完善。1961年Kaplan[8]將斷裂力學(xué)的概念引入到混凝土中。Lajtail[9]建立了壓剪斷裂的強(qiáng)變理論,指出在受壓剪作用下,裂縫端部形成不均勻應(yīng)力場,存在拉應(yīng)力和壓應(yīng)力集中現(xiàn)象,產(chǎn)生垂直受力方向的正剪切裂縫。在國內(nèi),徐道遠(yuǎn)等[10]在分析大體積混凝土破壞機(jī)理并進(jìn)行大量試驗的基礎(chǔ)上,提出了一種損傷-斷裂耦合模型,該模型能反映如混凝土壩這樣的大體積混凝土結(jié)構(gòu)開裂至破壞過程的實際情況。張鏡劍等[11]提出了一種適用于分析含壓剪裂縫結(jié)構(gòu)的非線性斷裂力學(xué)數(shù)值分析模型——壓剪鈍裂縫帶模型;楊延毅[12]運用損傷力學(xué)理論,從微裂縫損傷演化發(fā)展的角度研究混凝土裂縫的形成與擴(kuò)展機(jī)理,提出了相應(yīng)的理論判據(jù)?;炷两Y(jié)構(gòu)中不可避免地存在缺陷和裂縫,用連續(xù)介質(zhì)來反映應(yīng)混凝土材料的力學(xué)性能不符合實際情況。現(xiàn)今,許多學(xué)者將混凝土斷裂力學(xué)和損傷力學(xué)引入到混凝土裂縫的擴(kuò)展研究工作中。

    2 混凝土水力劈裂試驗研究

    從20世紀(jì)末開始,我國進(jìn)入了水利建設(shè)的黃金時期,一大批處于高水頭、大埋深等惡劣水文地質(zhì)條件下的高壩、深埋隧洞相繼開工建設(shè)。相對于復(fù)雜多樣的工程建設(shè)背景,巖體和混凝土水力劈裂的研究卻還處于初始階段,工程實踐缺少理論支持,急需深入開展水力劈裂的試驗以及計算方法研究。水力劈裂實驗是揭示水力劈裂機(jī)理的最佳方法,有關(guān)研究巖體單裂紋滲流特性與周圍應(yīng)力之間關(guān)系和裂紋水力劈裂機(jī)理的試驗較多[13-21],但研究混凝土材料斷裂細(xì)觀物理機(jī)制方面的試驗以及裂紋水力劈裂機(jī)理的試驗卻較少。由于混凝土的力學(xué)性能與試件尺寸有關(guān),要確切描述混凝土的力學(xué)性能就必須研究混凝土的尺寸效應(yīng)。很多實際的混凝土結(jié)構(gòu)尺寸均比較大,難以進(jìn)行實際結(jié)構(gòu)的系統(tǒng)試驗,僅能在實驗室進(jìn)行小尺寸試件的模擬,據(jù)此得到的結(jié)果對實際結(jié)構(gòu)的指導(dǎo)意義和實用性有限,成為目前研究的難題之一。吳智敏等[22-23]的試驗表明,混凝土的斷裂韌度隨試件尺寸的增加而增加,但當(dāng)試件尺寸增加到一定值時,斷裂韌度趨于恒定,即混凝土斷裂韌度只在一定范圍內(nèi)存在尺寸效應(yīng)。同時,徐世烺等[24]采用體積、跨度、厚度和高度四系列的三點彎曲梁試件研究了混凝土斷裂韌度的尺寸效應(yīng)規(guī)律,發(fā)現(xiàn)試件高度是影響混凝土斷裂韌度的主要原因。

    在混凝土水力劈裂機(jī)理研究方面,Bruhwiler 等[25-27]采用混凝土楔形劈裂試件對混凝土裂縫進(jìn)行了試驗研究和數(shù)值模擬,研究了裂縫中靜水壓力對混凝土表觀斷裂能和表觀斷裂韌度的影響,測試了裂紋內(nèi)水壓力分布,結(jié)果表明,裂紋內(nèi)水壓力會導(dǎo)致裂紋尖端損傷,降低裂尖部位的斷裂韌度。賈金生等[3]采用預(yù)設(shè)裂縫的全級配混凝土大試件進(jìn)行試驗,在裂縫面上直接施加水壓力,直至混凝土被高壓水擊穿。由于混凝土的力學(xué)行為與巖石材料有許多相似之處,故可以將有關(guān)巖石水力劈裂方面的試驗方法引入混凝土水力劈裂試驗研究工作中,李宗利等[28]針對巖體混凝土類材料從單裂紋的水力劈裂模型、試驗技術(shù)、數(shù)值方法、裂紋內(nèi)水壓力分布等方面介紹了國內(nèi)外的理論研究成果。

    1.4 腦電圖檢測 FS組患兒在入院第1天采用意大利ENG_GNT2SW腦電圖機(jī)行常規(guī)腦電圖描記。結(jié)果按照黃遠(yuǎn)貴編制的《臨床腦電圖學(xué)》分為正常、異常,異常包括非特異性異常(廣泛性慢波)及發(fā)作性異常[棘波、尖波、棘(尖)-慢復(fù)合波]。

    3 混凝土水力劈裂數(shù)值方法研究

    3.1 連續(xù)介質(zhì)力學(xué)分析法

    有限元法是當(dāng)今應(yīng)用最為廣泛的數(shù)值方法。為了解決裂縫的擴(kuò)展問題,Benzley[29]提出采用奇異單元法,Atluri等[30-31]提出采用增強(qiáng)單元法和動態(tài)奇異單元法。Swenson等[32]采用節(jié)點釋放技術(shù)模擬裂縫沿單元邊界擴(kuò)展,Bittencourt等[33]通過調(diào)整局部網(wǎng)格以實現(xiàn)對線彈性材料裂縫擴(kuò)展的模擬;Lim 等[34]應(yīng)用上述相近方法實現(xiàn)對彈塑性材料的復(fù)合型斷裂的模擬;楊慶生等[35]提出采用自適應(yīng)有限元法對裂縫擴(kuò)展進(jìn)行模擬分析;Poston等[36]通過引入接觸單元實現(xiàn)裂縫擴(kuò)展模擬,該方法在可能產(chǎn)生裂縫的區(qū)域布置厚度為零的接觸面單元,通過接觸面單元的應(yīng)力判斷裂縫的開裂情況及路徑。

    Belyschko等[37-38]等采用擴(kuò)展有限元法(extended finite element method,XFEM)解決了網(wǎng)格不斷調(diào)整的問題,將裂縫和有限元網(wǎng)格分開考慮,通過在位移函數(shù)中加入增強(qiáng)項以解決裂縫面的不連續(xù)性和裂尖的奇異性,從而解決裂縫擴(kuò)展問題。目前國內(nèi)有學(xué)者將XFEM應(yīng)用于混凝土重力壩的開裂分析[39-40],方修君等[41]還利用XFEM對裂隙水壓作用下混凝土試件的開裂過程進(jìn)行了模擬,首次將XFEM應(yīng)用于混凝土水力劈裂研究??梢钥闯?XFEM的應(yīng)用已經(jīng)由最初的梁、板等構(gòu)件拓寬到了大型水工結(jié)構(gòu)的開裂破壞研究,并且可以用于研究混凝土的水力劈裂破壞。董玉文等[42]研究了采用XFEM進(jìn)行重力壩水力劈裂數(shù)值模擬的方法。

    邊界元法[43]可以分別通過直接法和間接法建立邊界方程。為了處理裂縫擴(kuò)展問題,Sladek[44]采用拉普拉斯變化法分析了在沖擊荷載作用下,無限彈性內(nèi)含圓片裂縫的動態(tài)應(yīng)力強(qiáng)度因子;Portela[45]采用雙重邊界元法研究二維和三維的裂縫擴(kuò)展問題;Shou等[46]采用位移不連續(xù)法求解裂縫的不連續(xù)位移場和應(yīng)力場;Dwyer等[47]采用邊函數(shù)法,通過不同的勢函數(shù)考慮不同的邊界,如內(nèi)部空洞、裂縫等邊界;程玉民等[48-49]分別采用相似邊界元法和無限相似邊界元法,對瞬態(tài)載荷作用下的動態(tài)應(yīng)力強(qiáng)度因子進(jìn)行了分析計算。盡管邊界元法的研究已經(jīng)開展了很多年,但是邊界元法還是不適合用來處理非線性材料、多介質(zhì)等復(fù)雜問題。

    3.2 非連續(xù)介質(zhì)力學(xué)分析法

    數(shù)值流形方法(numerical manifold method, NMM)以拓?fù)淞餍魏臀⒎至餍螢榛A(chǔ),采用數(shù)學(xué)覆蓋和物理網(wǎng)格兩套獨立網(wǎng)格,將有限元網(wǎng)格作為數(shù)學(xué)覆蓋,邊界條件、裂縫、交界面等作為物理網(wǎng)格,通過物理網(wǎng)格對數(shù)學(xué)覆蓋的再剖分,形成覆蓋材料全域的物理覆蓋系統(tǒng)。此方法可以較好地模擬裂縫擴(kuò)展問題,Tsay等[50]采用數(shù)值流形方法和斷裂力學(xué)相結(jié)合的方法研究裂縫擴(kuò)展及尖端應(yīng)力場;王水林等[51]利用這種將連續(xù)變形與非連續(xù)變形分析相結(jié)合的數(shù)值流形法對裂縫擴(kuò)展進(jìn)行了數(shù)值模擬;Chiou 等[52]用虛位移擴(kuò)展法和數(shù)值流形方法相結(jié)合來研究混合裂縫擴(kuò)展問題。

    無單元法(EFM)是一種較為新穎的數(shù)值方法,它只需計算域的幾何邊界及計算點,不需要單元信息,因此具有邊界元法的優(yōu)點,且可以在開裂點附近布置可移動的加密節(jié)點以方便地跟蹤裂縫的擴(kuò)展過程,同時無單元法基本方程的數(shù)學(xué)基礎(chǔ)與有限元法類似,所以無單元法也具有有限元法的優(yōu)點,又比邊界元法具有更廣泛的應(yīng)用范圍。應(yīng)用無單元法追蹤裂紋擴(kuò)展,目前還存在著不少難題有待解決,在工程中的應(yīng)用也很少,國內(nèi)外許多學(xué)者都在進(jìn)行著這方面的研究:Belytschko等[53]提出了無單元伽遼金法(EFGM),這種方法在計算形函數(shù)導(dǎo)數(shù)時保留了被Nayroles等[54]忽略的所有項,并利用拉格朗日乘子法引入本質(zhì)邊界條件;在國內(nèi),清華大學(xué)張雄等[55]提出了基于子域法的緊支函數(shù)加權(quán)殘量法的無網(wǎng)格格式和最小二乘配點的無網(wǎng)格法,避免采用背景網(wǎng)格形式以消除其對數(shù)值積分的限制,并將此法應(yīng)用于求解彈塑性及波傳動等問題;欒茂田等[56-57]提出了采用有限覆蓋無單元法解決連續(xù)與非連續(xù)問題,并將此法應(yīng)用于裂縫擴(kuò)展追蹤問題、巖土類弱拉型材料摩擦接觸問題和脆性材料損傷斷裂演化行為問題等方面的研究;劉欣等[58]提出了流形覆蓋思想的無單元法,該法在邊界奇異性半解析無單元法研究的基礎(chǔ)上,采用四象限法以確定覆蓋范圍大小;周維垣等[59]對無單元伽遼金法進(jìn)行了詳細(xì)介紹,提出了一種拱壩開裂分析的近似方法(代表體積法),隨后又將無單元伽遼金法用于三維斷裂力學(xué)[60],并應(yīng)用于拱壩三維開裂分析中[61]。但是無單元法計算中系數(shù)矩陣帶寬比較大,致使計算量也相當(dāng)大,從計算精度及方便性來說,無單元法優(yōu)于有限元法,但就計算效率而言,無單元法不及有限元法,因此,將無單元法與有限元法耦合起來處理問題效果更好、更為方便,為解決高混凝土壩水力劈裂問題提供了全新的途徑。

    4 混凝土多場耦合模型研究

    4.1 滲流-應(yīng)力耦合模型

    滲流-應(yīng)力耦合模型主要包括滲流場和應(yīng)力場的數(shù)值模擬及裂縫特性的模擬兩個方面的內(nèi)容,其中耦合模型中的關(guān)鍵問題就是如何反映滲流過程中裂縫的擴(kuò)展和萌生。一直以來,關(guān)于混凝土水力劈裂方面的研究都只考慮裂縫的內(nèi)水壓力或者整體滲流場對結(jié)構(gòu)應(yīng)力、變形狀態(tài)的作用以及對裂縫擴(kuò)展的影響。混凝土結(jié)構(gòu)的應(yīng)力狀態(tài)和裂縫位置、長度、寬度以及形狀的變化也都會改變結(jié)構(gòu)的滲流場和裂縫的內(nèi)水壓力,因此滲流場與應(yīng)力場之間的相互作用在研究高混凝土壩水力劈裂時是不容忽視的。高壓水致使混凝土裂紋張開、擴(kuò)展、貫通的同時,伴隨著水流在裂紋中運動。裂紋內(nèi)水力特性及壓力直接影響混凝土的應(yīng)力應(yīng)變狀態(tài),尤其是對于在深水壓力下運行的高混凝土壩而言,建立其流固耦合模型存在一定的困難:①裂紋的失穩(wěn)擴(kuò)展也并非穩(wěn)定的發(fā)展,而是跳躍式前進(jìn),水流也并非恒定流。裂紋流體運動方程的求解邊界條件注入流量和外部水頭是恒定的,這樣的邊界條件不符合深埋地下結(jié)構(gòu)、高壩結(jié)構(gòu)等結(jié)構(gòu)中水流所產(chǎn)生的水力劈裂問題。②流固耦合分析時,涉及裂紋的張開寬度、張開速率及加速度、裂紋的滲透系數(shù)等參數(shù),實際應(yīng)用較為困難。

    目前滲流場與應(yīng)力場的耦合研究多半都以巖體,尤其是含裂隙的巖體作為研究對象。鑒于混凝土與巖體同為多孔介質(zhì),它們的內(nèi)部結(jié)構(gòu)在很多方面具有類似的特征,將巖體的滲流場與應(yīng)力場耦合模型應(yīng)用到混凝土結(jié)構(gòu)上也不失為一種有效的解決辦法。沈振中等[62-63]將無單元伽遼金法引入滲流場計算,運用無單元法求解重力壩以及巖體的滲流場與應(yīng)力場。將無單元法引入滲流-應(yīng)力耦合分析模型,研究耦合條件下高壩的水力劈裂問題是一種比較可行的方法。

    4.2 滲流-溫度耦合模型

    溫度荷載是混凝土壩的主要荷載之一,當(dāng)混凝土的溫度發(fā)生改變時,例如冷卻水的侵入、水泥水化熱以及環(huán)境溫度的改變,都將引起混凝土的應(yīng)力改變,當(dāng)約束拉應(yīng)變超過混凝土的極限拉伸值時,混凝土結(jié)構(gòu)發(fā)生開裂破壞。溫度變化引起熱應(yīng)力和熱應(yīng)變使得應(yīng)力場和位移場發(fā)生改變,同時隨著溫度的變化,材料的自身物理力學(xué)性質(zhì)也會發(fā)生改變,從而引起水工混凝土裂縫的特性發(fā)生改變,甚至引起裂縫的進(jìn)一步擴(kuò)展。在溫度場分析中,由于水等滲透介質(zhì)在滲流場中的運移過程中,必然存在與壩體的熱量傳遞與交換,導(dǎo)致水工結(jié)構(gòu)內(nèi)的溫度場分布發(fā)生一定的改變,因此滲流場對溫度場的作用是不可忽視的。當(dāng)前滲流-溫度耦合作用大多都是針對地下巖體進(jìn)行研究的,Harlan[64]考慮滲流場和溫度場的耦合作用,提出了凍結(jié)過程中的熱質(zhì)遷移數(shù)學(xué)模型;柴軍瑞等[65]采用等效連續(xù)介質(zhì)法,提出滲流和溫度耦合模型,并給出了有限元法計算方法;黃濤等[66]通過將裂隙圍巖體滲透性能與熱物理性能的等效連續(xù)化處理,建立了圍巖體溫度場與滲流場耦合作用數(shù)學(xué)模型,計算了某隧洞的溫度和滲流場;賴遠(yuǎn)明等[67]應(yīng)用伽遼金法導(dǎo)出相變溫度場和滲流場的有限元方程,并計算了寒區(qū)大壩和擋土墻的滲流場和溫度場;周志芳等[68]建立了考慮熱對流、熱傳導(dǎo)、熱機(jī)械彌散及含水介質(zhì)骨架與水之間的熱量交換的地下水熱量運移數(shù)學(xué)模型,運用BEM-FAM耦合法求解熱量運移的水流-傳熱耦合模型,具有較好的計算精度;張樹光[69]提出滲流改變了溫度場的對稱分布,熱交換隨著滲流速度的增加而向順滲流方向移動;張玉軍[70]采用裂隙巖體的滲透性能和熱學(xué)性能等效連續(xù)處理方法,指出由于裂隙中的滲流通過與巖體的熱量傳遞使得溫度場的分布發(fā)生改變。到目前為止,關(guān)于混凝土裂縫的滲流-溫度耦合分析研究的報道還很少。

    4.3 溫度滲流-應(yīng)力耦合模型

    目前溫度-滲流-應(yīng)力三場耦合大多應(yīng)用于核能、地?zé)嵋约昂藦U料填埋等方面,三場耦合理論的發(fā)展過程中,國內(nèi)外研究人員提出了多種形式的三場耦合數(shù)學(xué)模型,大致可分為兩類:飽和介質(zhì)三場耦合數(shù)學(xué)模型與非飽和-飽和介質(zhì)三場耦合數(shù)學(xué)模型。Noorishad等[71]在擴(kuò)展的Biot固結(jié)理論基礎(chǔ)上,提出了飽和孔隙介質(zhì)的三場耦合控制方程組;Hart 等[72]基于混合物理論,推導(dǎo)了飽和孔隙介質(zhì)的熱力水流完全耦合方程組;Gatmiri等[73]提出了飽和土的熱液力耦合微分方程組;Bower等[74]提出了飽和雙重介質(zhì)的三場耦合數(shù)學(xué)模型;Neaupane 等[75]推導(dǎo)了各向異性非飽和飽和-孔隙介質(zhì)的三場耦合控制方程;賴遠(yuǎn)明等[76]根據(jù)凍土力學(xué)、滲流理論和傳熱學(xué)提出了帶相變的三場耦合數(shù)學(xué)力學(xué)模型;黃濤等[77]通過等效性能場之間耦合作用的研究建立了巖體地下水滲流場、應(yīng)力場與溫度場耦合作用的簡單數(shù)學(xué)模型;劉亞晨等[78]根據(jù)不可逆熱力學(xué)基本理論,推導(dǎo)了核廢料貯庫裂隙巖體介質(zhì)熱液力耗散過程的定解方程;梁冰等[79]在塊裂介質(zhì)理論的基礎(chǔ)上,導(dǎo)出了三維裂隙網(wǎng)絡(luò)的塊裂介質(zhì)高溫巖體熱-流-固耦合方程;Rutqvist等[80]導(dǎo)出了對于飽和-非飽和介質(zhì)三場耦合控制方程的一般形式。

    三場耦合模型的數(shù)值模擬由于能更好地體現(xiàn)結(jié)構(gòu)的運行工作狀態(tài),已經(jīng)成為很多研究者的重要研究方向,但是有關(guān)水工混凝土結(jié)構(gòu)裂縫擴(kuò)展方面的研究卻很少。水工混凝土結(jié)構(gòu)大多都是在水壓力作用和溫度變化條件下運行的,而溫度場、滲流場和應(yīng)力場的分布是評價水工混凝土結(jié)構(gòu)工程運行狀況的重要方面,同時三者之間并非是孤立存在的,它們之間是相互作用、相互影響和相互制約的,當(dāng)某一場發(fā)生擾動變化時,另外兩場也將受到一定程度的影響。要想準(zhǔn)確、客觀和全面評價水工混凝土結(jié)構(gòu)的裂縫萌生、擴(kuò)展以及可能的失穩(wěn)狀況,應(yīng)當(dāng)考慮溫度、滲流和應(yīng)力之間的相互作用關(guān)系,即考慮溫度場、滲流場、應(yīng)力場三場間的耦合作用。因此,需要研究多場耦合條件下的高混凝土壩水力劈裂問題就必須建立合理的混凝土裂縫溫度-滲流-應(yīng)力多場耦合模型。

    5 尚待研究的問題

    目前,多場耦合作用下水工混凝土結(jié)構(gòu)水力劈裂試驗、發(fā)生機(jī)理以及高壓水力劈裂數(shù)值模擬等方面的研究不多,尚有許多有待進(jìn)一步研究的問題:

    a.研究水工混凝土裂縫在外荷載和水壓力的雙重影響下,水工混凝土在非恒定水頭和流量時裂縫擴(kuò)展與流體的耦合機(jī)理,研究滲流-應(yīng)力耦合條件下裂縫的萌生、發(fā)展以及失穩(wěn)機(jī)理和特性。

    b.在溫度滲流應(yīng)力耦合作用條件下,其耦合機(jī)理復(fù)雜,需要深入研究水工混凝土結(jié)構(gòu)裂縫開度、擴(kuò)展長度、裂縫內(nèi)的內(nèi)水壓力、起裂韌度、裂縫在不同水壓情況下的起裂寬度、外部水溫等對混凝土結(jié)構(gòu)裂縫擴(kuò)展的影響。

    c.采用數(shù)值模擬方法研究多場耦合作用下水工混凝土裂縫產(chǎn)生和發(fā)展過程,需要重點解決多裂縫、交叉裂縫以及任意裂縫開展的計算技術(shù),同時其斷裂準(zhǔn)則也有待進(jìn)一步深入研究。

    d.針對高混凝土壩水力劈裂問題,需要研制先進(jìn)的測試技術(shù)和試驗設(shè)備,以及研發(fā)混凝土裂縫起裂判斷和擴(kuò)展過程的相關(guān)計算程序,注重多場耦合作用和結(jié)構(gòu)裂縫的任意性,研究高混凝土壩水力劈裂的萌生、發(fā)展和失穩(wěn)破壞的演化規(guī)律。

    [1]朱伯芳.混凝土壩溫度控制與防止裂縫的現(xiàn)狀與展望[J].水利學(xué)報,2006,37(12):1424-1432.(ZHU Bofang.Current situation and prospect of temperature control and cracking prevention technology for concrete dam[J].JournalofHydraulicEngineering,2006,37 (12):1424-1432.(in Chinese))

    [2]徐道遠(yuǎn),王向東,朱為玄,等.混凝土壩的損傷及損傷仿真計算[J].河海大學(xué)學(xué)報:自然科學(xué)版,2002,30(4): 14-17.(XU Daoyuan,WANG Xiangdong,ZHU Weixuan. Damage to concrete dams and its simulation calculation [J].Journal of Hohai University:Natural Sciences,2002, 30(4):14-17.(in Chinese))

    [3]賈金生,李新宇,鄭璀瑩.特高重力壩考慮高壓水劈裂影響的初步研究[J].水利學(xué)報,2006,37(12):1509-1515.(JIA Jinsheng,LI Xinyu,ZHENG Cuiying.Studies on safety problem of high gravity dams higher than 200 m with consideration of hydraulic fracturing under water pressure[J].Journal of Hydraulic Engineering,2006,37 (12):1509-1515.(in Chinese))

    [4]李雪紅.重大水工混凝土結(jié)構(gòu)裂縫演變規(guī)律及轉(zhuǎn)異診斷方法研究[D].南京:河海大學(xué),2003.

    [5]朱萬成,趙啟林,唐春安,等.混凝土斷裂過程的力學(xué)模型與數(shù)值模擬[J].力學(xué)進(jìn)展,2002,32(4):579-598. (ZHU Wancheng,ZHAO Qilin,TANG Chun’an,et al. Mechanical model and numerical simulation of fractur eprocess of concrete[J].Advances in Mechanics,2002,32 (4):579-598.(in Chinese))

    [6]尹雙增.斷裂判據(jù)與在混凝土工程中應(yīng)用[M].北京:科學(xué)出版社,1996:15-18.

    [7]蔡四維,蔡敏.混凝土的損傷斷裂[M].北京:人民交通出版社,2000:55-65.

    [8]KAPLAN M F.Crack propagation and the fracture of concrete[J].ACI Journal,1961,58(11):591-610.

    [9]LAJTAL E Z.Brittle fracture in compression[J]. International Journal of Fracture,1977,10(4):195-208.

    [10]徐道遠(yuǎn),符曉陵,朱為玄,等.壩體混凝土損傷斷裂模型[J].大連理工大學(xué)學(xué)報,1997,37(增刊1):1-6.(XU Daoyuan,WANG FU Xiaoling,ZHU Weixuan,et al. Damage-fracture model for concrete of dam[J].Journal of Dalian University of Technology,1997,37(Sup1):1-6. (in Chinese))

    [11]張鏡劍,涂金良,翟奇愚,等.龍灘碾壓混凝土重力壩水平層面穩(wěn)定安全的斷裂分析[J].華北水利水電學(xué)院學(xué)報,1996,17(4):2-7.(ZHANG Jingjian,TU Jinliang, ZHAI Qiyu,et al.Fracture analysis of the stability safety of the horizontal construction joint in longtan rollcrete gravity dam[J].JournalofNorthChinaInstituteofWater Conservancy and Hydroelectric Power,1996,17(4):2-7. (in Chinese))

    [12]楊延毅.混凝土損傷斷裂過程研究[J].浙江大學(xué)學(xué)報:自然科學(xué)版,1993,27(5):654-661.(YANG Yanyi.The research of progressive cracking in concrete structure with damage theory[J].Journal of Zhejiang University:Natural Science,1993,27(5):654-661.(in Chinese))

    [13]BARTONN,BANDIS,BAKHTURK.Strength deformation and conductivity coupling of rock joints[J]. Int J Rock Mechanics Mining Science and Geomechanics, 1985,22:121-140.

    [14]耿克勤,陳鳳翔,劉光廷,等.巖體裂隙滲流水力特性的實驗研究[J].清華大學(xué)學(xué)報:自然科學(xué)版,1996,36(1):102-106.(GENG Keqin,CHEN Fengxiang,LIU Guangting,etal.Experimentalresearchofhydraulic properties of seepage flow in fracture[J].Journal of Tsinghua University:Science and Technology,1996,36 (1):102-106.(in Chinese))

    [15]ESAKI T,DU S,MITANI Y,et al.Development of a shear flow test apparatus and determination of coupled properties for a single rock joint[J].Int J Rock Mechanics and Mining Science,1999,36:641-650.

    [16]LEE H S,CHO T F.Characteristics of rough fractures in linear flow under normal and shear load[J].Int J Rock Mechanics and Rock Engineering,2002,35(4):299-318.

    [17]速寶玉,詹美禮,王媛.裂隙滲流與應(yīng)力耦合特性的試驗研究[J].巖土工程學(xué)報,1997,19(4):1-3.(SU Baoyu,ZHAN Meili,WANG Yuan.Experimental research of coupling properties of seepage and stress in fracture [J].Chinese Jounal of Geotechnical Engineering,1997,19 (4):1-3.(in Chinese))

    [18]朱珍德,徐衛(wèi)亞,張愛軍.脆性巖石損傷斷裂機(jī)理分析與試驗研究[J].巖石力學(xué)與工程學(xué)報,2003,22(9): 1411-1416.(ZHU Zhende,XU Weiya,ZHANG Aijun. Mechanism analysis and testing study on damage and fracture of brittle rock[J].Chinese Journal of Rock Mechanics and Engineering,2003,22(9):1411-1416.(in Chinese))

    [19]謝興華.巖體水力劈裂機(jī)理試驗及數(shù)值模擬研究[D].南京:河海大學(xué),2004.

    [20]王國慶.巖體水力劈裂試驗及裂紋擴(kuò)展的無單元法計算[D].南京:河海大學(xué),2004.

    [21]陳衛(wèi)忠,朱維申,羅超文.萬家寨引黃工程總干一、二級泵站水力劈裂試驗研究[J].巖土力學(xué),2001,22(1): 26-28.(CHEN Weizhong,ZHU Weishen,LUO Chaowen. Hydraulic jacking tests at GMPS1 and GMPS2 of Shanxi YellowRiverdiversionproject[J].RockandSoil Mechanics,2001,22(1):26-28.(in Chinese))

    [22]吳智敏,趙國藩,徐世烺.大尺寸混凝土試件的斷裂韌度[J].水利學(xué)報,1997(6):67-76.(WU Zhimin,ZHAO Guofan,XU Shilang.Fracture toughness of large concrete specimen[J].Journal of Hydraulic Engineering,1997 (6):67-76.(in Chinese))

    [23]吳智敏,趙國藩.混凝土斷裂韌度的尺寸效應(yīng)研究[J].工業(yè)建筑,1995,25(4):20-22.(WU Zhimin,ZHAO Guofan.Studies on size effect of fracture toughness of concrete[J].Industrial Construction,1995,25(4):20-22. (in Chinese))

    [24]徐世烺,趙國藩.混凝土斷裂力學(xué)研究[M].大連:大連理工大學(xué)出版社,1991.

    [25]BRUHWILER E,SAOUMA V.Water fracture interaction in concrete:partⅠ:fracture properties[J].ACI Materials Journal,1995,92(3):296-303.

    [26]BRUHWILER E,SAOUMA V.Water fracture interaction in concrete:partⅡ:hydrostatic pressure cracks[J].ACI Materials Journal,1995,92(4):383-390.

    [27]VOLKER V,SAOUMA V.Water pressure in propagating concrete cracks[J].Journal of Structure Engineering, 2000,126(2):235-242.

    [28]李宗利,任青文.巖石混凝土類材料單裂紋水力劈裂研究述評[J].水利水運工程學(xué)報,2005(1):67-74.(LI Zongli,REN Qingwen.Review of hydraulic fracturing on a single fracture of rock and concrete materials[J].Hydro-Science and Engineering,2005(1):67-74.(in Chinese))

    [29]BENZLEYSE.Representationofsingularitieswith isoperimetric finite elements[J].International Journal for Numerical Methods in Engineering,1974(8):537-545.

    [30]ATLURI S N.An assumed displacement hybrid finite element modelforlinearfracturemechanics[J]. International Journal of Fracture,1975,11(2):257-271.

    [31]NISHIOKA T,ATLURI S N.Path-independent integral and movingisoperimetricelementsfordynamicscrack propagation[J].AIAA Journal,1984,22(3):409-414.

    [32]SWENSON D V,INGRAFFEA A R.Modeling mixed-mode dynamic crack propagation using finite element:theory and applications[J].Computational Mechanics,1988(3):381-397.

    [33]BITTENCOURT T U,WAWRZYNEK P A,INGRAFFEA A R.Quasi-automatic simulation of crack propagation for 2D LEFM problems[J].Engineering Fracture Mechanics, 1996,55(2):321-334.

    [34]LIM I L,JOHNSTON I W,CHOI S K.A finite element code for fracture propagation analysis within elasto-plastic continuum[J].Engineering Fracture Mechanics,1996,53 (2):193-211.

    [35]楊慶生,楊衛(wèi).斷裂過程的有限元模擬[J].計算力學(xué)學(xué)報,1997,14(4):79-88.(YANG Qingsheng,YANG Wei. Finite element simulation of fracture process[J].Chinese Journal of Computational Mechanics,1997,14(4):79-88. (in Chinese))

    [36]POSTON T,STEWART I.Catastrophe theory and its applications[M].New York:Pitman,1978.

    [37]BELYSCHKO T,BLACK T.Elastic crack growth in finite elements with minimal remeshing[J].International Journal for Numerical Methods in Engineering,1999,45:601-620.

    [38]BELYSCHKOT,MOESN,USUIS.Arbitrary discontinuities in finite elements[J].Numer Methods Engrg,2001,50:993-1013.

    [39]杜效鵠,段云嶺,王光綸.重力壩斷裂數(shù)值分析研究[J].水利學(xué)報,2005,36(9):1035-1042.(DU Xiaohu, DUAN Yunling,WANG Guanglun.Numerical analysis of fracture ingravitydam[J].JournalofHydraulic Engineering,2005,36(9):1035-1042.(in Chinese))

    [40]方修君,金峰,王進(jìn)廷.基于擴(kuò)展有限元法的Koyna重力壩地震開裂過程模擬[J].清華大學(xué)學(xué)報:自然科學(xué)版,2008,48(12):2065-2069.(FANG Xiujun,JIN Feng, WANG Jinting.Seismic fracture simulation of the Koynagravity dam using an extended finite element method[J]. Journal of Tsinghua University:Science and Technology, 2008,48(12):2065-2069.(in Chinese))

    [41]方修君,金峰.裂隙水流與混凝土開裂相互作用的耦合模型[J].水利學(xué)報,2007,38(12):1466-1474.(FANG Xiujun,JIN Feng.Coupling model for interaction between fissure water and cracking in concrete[J].Journal of Hydraulic Engineering,2007,38(12):1466-1474.(in Chinese))

    [42]董玉文,任青文.重力壩水力劈裂分析的擴(kuò)展有限元法[J].水利學(xué)報,2011,42(11):1361-1367.(DONG Yuwen,REN Qingwen.An extended finite element method for modeling hydraulic fracturing in gravity dam[J]. Journal of Hydraulic Engineering,2011,42(11):1361-1367.(in Chinese))

    [43]姚振漢,杜慶華.邊界元法應(yīng)用的若干近期研究及國際新進(jìn)展[J].清華大學(xué)學(xué)報:自然科學(xué)版,2001,41(4/ 5):89-93.(YAO Zhenhan,DU Qinghua.Some recent investigations and new progresses in the application of boundary element methods[J].Journal of Tsinghua University:Science and Technology,2001,41(4/5):89-93.(in Chinese))

    [44]SLADEK J.Dynamics stress intensity fractures studied by boundary integro-differential equations[J].International Journal Numerical Methods in Engineering,1986,23(5): 919-928

    [45]PORTELA A.The dual BEM:effective implementation for crack problem[J].International Journal for Numerical Methods in Engineering,1992,33:269-287.

    [46]SHOU K J,CROUCH S L.A higher order displacement discontinuity method for analysis of crack problems[J]. International Journal of Rock Mechanics and Mining Sciences&Geomechanics Abstracts,1995,32(1):49-55.

    [47]DWYER J F,AMADAI B.Application of the edge function method to rock mechanics problems[J].Rock Mechanics and Rock Engineering,1995,28(4):185-209.

    [48]程玉民,彭妙娟.彈性動力學(xué)的相似邊界元法[J].重慶建筑大學(xué)學(xué)報,2000,22(6):1-4.(CHENG Yumin, PENG Miaojuan.Similar boundary element method in elastodynamics[J].Journal of Chongqing Jianzhu University, 2000,22(6):1-4.(in Chinese))

    [49]程玉民,嵇醒,賀鵬飛.動態(tài)斷裂力學(xué)的無限相似邊界元法[J].力學(xué)學(xué)報,2004,36(1):43-48.(CHENG Yumin,JI Xing,HE Pengfei.Infinite similar boundary element method for dynamic fracture mechanics[J]. Chinese Journal of Theoretical and Applied Mechanics, 2004,36(1):43-48.(in Chinese))

    [50]TSAY R J,CHIOU Y J,CHUANG W L.Crack growth prediction by manifold method[J].Journal of Engineering Mechanics,1999,125(8):884-890.

    [51]王水林,葛修潤.流形元方法在模擬裂紋擴(kuò)展中的應(yīng)用[J].巖石力學(xué)與工程學(xué)報,1997,16(5):405-410. (WANG Shuilin,GE Xiurun.Application of manifold method in simulating crack propagation[J].Chinese Journal of Rock mechanics and Engineering,1997,16 (5):405-410.(in Chinese))

    [52]CHIOU Y J,LEE Y M,TSAY R J.Mixed mode fracture propagation by manifold method[J].International Journal of Fracture,2002,114:327-347.

    [53]BELYTSCHKO T,LU Y Y,GU L.Element free Galerkin methods[J].Int J Num Meth Engng,1994,37:229-256.

    [54]NAYROLES B,TOUZOT G,VILLON P.Generalizing the finite element method:diffuse approximation and diffuse elements[J].Computational Mechanics,1992,10(5): 307-318.

    [55]張雄,宋康祖,陸明萬.緊支試函數(shù)加權(quán)殘值法[J].力學(xué)學(xué)報,2003,35(1):43-49.(ZHANG Xiong,SONG Kangzu,LU Mingwan.Weighted residual method with compactly supported trial functions[J].Chinese Journal of Theoretical and Applied Mechanics,2003,35(1):43-49. (in Chinese))

    [56]欒茂田,田榮,楊慶,等.有限覆蓋無單元法在巖土類弱拉型材料摩擦接觸問題中的應(yīng)用[J].巖土工程學(xué)報, 2002,24(2):137-141.(LUAN Maotian,TIAN Rong, YANGQing,etal.Applicationsofthefinite-cover element-free methodinfrictioncontactproblemsof tension-weakgeomaterials[J].ChineseJounalof Geotechnical Engineering,2002,24(2):137-141.(in Chinese))

    [57]田榮.連續(xù)與非連續(xù)變形分析的有限覆蓋無單元方法及其應(yīng)用研究[D].大連:大連理工大學(xué),2000.

    [58]劉欣,朱德懋,陸明萬,等.平面裂紋問題的h,p,hp型自適應(yīng)無網(wǎng)格方法的研究[J].力學(xué)學(xué)報,2000,32 (3):308-318.(LIU Xin,ZHU Demao,LU Mingwan,et al.h,p,hp adaptive meshless method for plane crack problem[J].Chinese Journal of Theoretical and Applied Mechanics,2000,32(3):308-318.(in Chinese))

    [59]周維垣,寇曉東.無單元法及其工程應(yīng)用[J].力學(xué)學(xué)報,1998,30(2):193-202.(ZHOUWeiyuan,KOU Xiaodong.Element free method and its application in engineering[J].ChineseJournalofTheoreticaland Applied Mechanics,1998,30(2):193-202.(in Chinese))

    [60]胡云進(jìn),周維垣,林鵬.無單元法在三維斷裂力學(xué)中的應(yīng)用[J].巖土力學(xué),2003,24(增刊1):21-24.(HU Yunjin,ZHOU Weiyuan,LIN Peng.Application of EFG method to three-dimensional fracture mechanics[J].Rock andSoilMechanics,2003,24(Sup1):21-24.(in Chinese))

    [61]周維垣,黃巖松,林鵬.三維無單元伽遼金法及其在拱壩分析中的應(yīng)用[J].水利學(xué)報,2005,36(6):645-649. (ZHOU Weiyuan,HUANG Yansong,LIN Peng.Threedimensionalelement-freeGalerkinmethodandits application in analysis of arch dams[J].Journal of Hydraulic Engineering,2005,36(6):645-649.(inChinese))

    [62]沈振中,陳小虎,徐力群.重力壩應(yīng)力-滲流相互作用的無單元耦合分析[J].巖土力學(xué),2008,29(增刊1):74-78.(SHENZhenzhong,CHENXiaohu,XULiqun. Coupled analysis of stress-seepage interaction for gravity dambyelement-freemethod[J].RockandSoil Mechanics,2008,29(Sup1):74-78.(in Chinese))

    [63]沈振中,張鑫,孫粵琳.巖體水力劈裂的應(yīng)力-滲流-損傷耦合模型研究[J].計算力學(xué)學(xué)報,2009,26(4):523-528.(SHEN Zhenzhong,ZHANG Xin,SUN Yuelin. Research on stress-seepage-damage coupling model of hydraulic fracturing for rock mass[J].Chinese Journal of Computational Mechanics,2009,26(4):523-528.(in Chinese))

    [64]HARLAN R L.Analysis of coupled heat-fluid transport in partial frozen soil[J].Water Resources Research,1973,9 (5):1314-1323.

    [65]柴軍瑞,韓群柱.巖體滲流場與溫度場耦合的連續(xù)介質(zhì)模型[J].地下水,1997,19(2):59-62.(CHAI Junrui, HAN Qunzhu.Continuous medium model for the coupled problem of temperature and seepage fields for rock mass [J].Groundwater,1997,19(2):59-62.(in Chinese))

    [66]黃濤,楊立中.隧道裂隙巖體溫度-滲流耦合數(shù)學(xué)模型研究[J].巖土工程學(xué)報,1999,21(5):554-558. (HUANG Tao,YANG Lizhong.A study of mathematical modeloncouplingbetweentemperature-seepagein fractured rock mass surrounding tunnel[J].Chinese Jounal of Geotechnical Engineering,1999,21(5):554-558.(in Chinese))

    [67]賴遠(yuǎn)明,劉松玉,鄧華鈞.寒區(qū)大壩溫度場和滲流場耦合問題的非線性數(shù)值模擬[J].水利學(xué)報,2001(8):26-31. (LAI Yuanming,LIU Songyu,DENG Xuejun.Numerical simulation for the coupled problem of temperature and seepage fields in dams located in cold area[J].Journal of Hydraulic Engineering,2001(8):26-31.(in Chinese))

    [68]周志芳,王錦國.裂隙介質(zhì)水動力學(xué)[M].北京:中國水利水電出版社,2004.

    [69]張樹光.深埋巷道圍巖溫度場的數(shù)值模擬分析[J].科學(xué)技術(shù)與工程,2006,6(14):2194-2196.(ZHANG Shuguang.Numerical simulation analysis of temperature field in surrounding rock of deeply buried roadways[J]. Science Technology and Engineering,2006,6(14):2194-2196.(in Chinese))

    [70]張玉軍.核廢料地質(zhì)處置近場熱水應(yīng)力遷移耦合二維有限元分析[J].巖土工程學(xué)報,2007,29(10):1553-1557.(ZHANG Yujun.2D FEM analysis for coupled thermo-hydro-mechanical-migratory process in near field of geological disposal of nuclear waste[J].Chinese Journal of Geotechnical Engineering,2007,29(10):1553-1557. (in Chinese))

    [71]NOORISHAD J,TSANG C F,WITHERSPOON P A. Coupledthermal-hydraulic-mechanicalphenomenain saturated fractured porous rocks:numerical approach[J].J Geophysical Res,1984,89(B12):10365-10373.

    [72]HART R D,JOHN C M.Formulation of a fully-coupled thermal-mechanical-fluid model for non-linear geologic systems[J].Int J Rock Mech Min Sic&Geornech Abstr, 1986,23(3):213-224.

    [73]GATMIRI B,DELAGE P.A formulation of fully coupled thermal-hydro-mechanical behavior of saturated porous media-numerical approach[J].Int J Num&Anal Methods in Geomech,1997,21(3):199-225.

    [74]BOWER K M,ZYVOLOSKI G.A numerical model for thermo-hydro-mechanical coupling in fractured rock[J]. Int J Rock Mech Min Sic,1997,34(8):1201-1211.

    [75]NEAUPANEKM,YAMABET,YOSHINAKAR. Simulation of a fully coupled thermo-hydro-mechanical system in freezing and thawing rock[J].Int J Rock Mech Min Sic,1999,36(5):563-580.

    [76]賴遠(yuǎn)明,吳紫汪,朱元林,等.寒區(qū)隧道溫度場、滲流場和應(yīng)力場耦合問題的非線性分析[J].巖土工程學(xué)報, 1999,21(5):529-533.(LAI Yuanming,WU Ziwang, ZHU Yuanlin,et al.Nonlinear analyses for the couple problem of temperature,seepage and stress fields in cold region tunnels[J].ChineseJounalofGeotechnical Engineering,1999,21(5):529-533.(in Chinese))

    [77]黃濤,楊立中,陳一立.工程巖體地下水滲流-應(yīng)力-溫度耦合作用數(shù)學(xué)模型的研究[J].西南交通大學(xué)學(xué)報, 1999,34(1):11-15.(HUANG Tao,YANG Lizhong, CHEN Yili.A mathematical model for the coupling among groundwater seepage stress temperature in engineering rock mass[J].Journal of Southwest Jiaotong University, 1999,34(1):11-15.(in Chinese))

    [78]劉亞晨,劉泉聲,吳玉山,等.核廢料貯庫圍巖介質(zhì)不可逆過程熱力學(xué)和熱彈性[J].巖石力學(xué)與工程學(xué)報, 2000,19(3):361-365.(LIU Yachen,LIU Quansheng, WU Yushan,etal.Irreversiblethermodynamicsand thermoelasticity of fractured rock mass surrounding nuclear waste repositories[J].Chinese Journal of Rock Mechanics and Engineering,2000,19(3):361-365.(in Chinese))

    [79]梁冰,孫可明,薛強(qiáng).地下工程中的流-固耦合問題的探討[J].遼寧工程技術(shù)大學(xué)學(xué)報.2001,20(2):129-34. (LIANG Bing,SUN Keming,XUE Qiang.The research of fluid-solid coupling in the ground engineering[J].Journal of Liaoning Technical University:Natural Science Edition, 2001,20(2):129-34.(in Chinese))

    [80]RUTQVIST J,BBRGESSON L,CHIJIMATSU M,et al. Thermo-hydro-mechanics of partially saturated geological media:governing equations and formulation of four finite element models[J].Int J Rock Mech Min Sic,2001,38 (1):105-127.

    Review on hydraulic fracture in high concrete dam under multi-field coupling conditions

    //GAN Lei1,2,SHEN Zhenzhong1,2,XU Liqun1,2(1.College of Water Conservancy and Hydropower Engineering,Hohai University,Nanjing 210098,China;2.State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering,Hohai University, Nanjing 210098,China)

    The research status of hydraulic fractures in concrete structures is reviewed.With regard to the existing problems in the research on hydraulic fractures in high concrete dams under multi-field coupling conditions,the worldwide results of hydraulic fractures in concrete structures are summarized from the perspectives of hydraulic fracturing theories,test techniques,numerical simulation methods,and crack propagation studies under different coupling conditions.Several prospects are proposed:test equipment and techniques for hydraulic fractures in high concrete dams under multi-field coupling conditions,crack propagation and coupling mechanisms for different coupling conditions,and numerical techniques considering arbitrary structural cracks.

    high concrete dam;hydraulic fracture;multi-filed coupling;crack propagation;review

    10.3880/j.issn.10067647.2013.02.019

    TV642

    A

    10067647(2013)02008708

    2012-06-14 編輯:熊水斌)

    國家自然科學(xué)基金(51179062);江蘇省普通高校研究生科研創(chuàng)新計劃(CXLX11_0447)

    甘磊(1987—),男,博士研究生,江西撫州人,主要從事壩工設(shè)計計算理論與試驗技術(shù)研究。E-mail:ganlei2015@hhu.edu.cn

    猜你喜歡
    水力滲流裂紋
    水力全開
    Epidermal growth factor receptor rs17337023 polymorphism in hypertensive gestational diabetic women: A pilot study
    微裂紋區(qū)對主裂紋擴(kuò)展的影響
    球墨鑄鐵管的水力計算
    水力噴射壓裂中環(huán)空水力封隔全尺寸實驗
    預(yù)裂紋混凝土拉壓疲勞荷載下裂紋擴(kuò)展速率
    簡述滲流作用引起的土體破壞及防治措施
    河南科技(2014年12期)2014-02-27 14:10:26
    關(guān)于渠道滲流計算方法的選用
    河南科技(2014年11期)2014-02-27 14:09:48
    低水力停留時間氧化溝的改造與調(diào)控
    低合金鋼焊接裂紋簡述
    国产亚洲精品第一综合不卡| 成年女人毛片免费观看观看9 | 日日爽夜夜爽网站| 亚洲中文日韩欧美视频| 99九九在线精品视频| 高清欧美精品videossex| 精品国产国语对白av| 手机成人av网站| 成人av一区二区三区在线看 | 亚洲国产av影院在线观看| 色婷婷av一区二区三区视频| 可以免费在线观看a视频的电影网站| 自拍欧美九色日韩亚洲蝌蚪91| 国产成人啪精品午夜网站| 国产日韩一区二区三区精品不卡| 丝袜美足系列| www日本在线高清视频| 99国产精品一区二区三区| 曰老女人黄片| 在线观看www视频免费| 男人舔女人的私密视频| 三上悠亚av全集在线观看| 国产亚洲一区二区精品| 久久久久久免费高清国产稀缺| 亚洲成人免费电影在线观看| 国产成人免费无遮挡视频| 亚洲va日本ⅴa欧美va伊人久久 | 欧美变态另类bdsm刘玥| 宅男免费午夜| 90打野战视频偷拍视频| a级毛片在线看网站| 视频区图区小说| 爱豆传媒免费全集在线观看| 正在播放国产对白刺激| 日日摸夜夜添夜夜添小说| 我要看黄色一级片免费的| 亚洲中文日韩欧美视频| 欧美性长视频在线观看| 纯流量卡能插随身wifi吗| 久久毛片免费看一区二区三区| 久久精品人人爽人人爽视色| av网站免费在线观看视频| 精品熟女少妇八av免费久了| 久久久精品94久久精品| 国产成+人综合+亚洲专区| 老汉色∧v一级毛片| av欧美777| 国产精品自产拍在线观看55亚洲 | 亚洲第一青青草原| 51午夜福利影视在线观看| 国产成人影院久久av| 久热这里只有精品99| 蜜桃在线观看..| 亚洲av电影在线进入| 国产男女超爽视频在线观看| 日韩人妻精品一区2区三区| 国产成人免费观看mmmm| 精品国产乱子伦一区二区三区 | 亚洲精品久久久久久婷婷小说| 菩萨蛮人人尽说江南好唐韦庄| 亚洲情色 制服丝袜| 老司机亚洲免费影院| 成人国产一区最新在线观看| 亚洲av美国av| av网站免费在线观看视频| 亚洲成国产人片在线观看| 国产1区2区3区精品| 精品国产超薄肉色丝袜足j| 人妻 亚洲 视频| 中文字幕另类日韩欧美亚洲嫩草| 在线观看免费午夜福利视频| 黄片播放在线免费| 超碰成人久久| 乱人伦中国视频| 男女床上黄色一级片免费看| 国产精品久久久久久人妻精品电影 | 久久久久久久久久久久大奶| 精品一品国产午夜福利视频| 久久精品aⅴ一区二区三区四区| 丝袜喷水一区| 久久久国产精品麻豆| 99精国产麻豆久久婷婷| 三上悠亚av全集在线观看| 伊人亚洲综合成人网| 国产色视频综合| 国产亚洲午夜精品一区二区久久| 亚洲va日本ⅴa欧美va伊人久久 | 51午夜福利影视在线观看| 国产一区二区三区av在线| 国产高清国产精品国产三级| 色婷婷久久久亚洲欧美| 亚洲 欧美一区二区三区| 国产精品一区二区精品视频观看| 曰老女人黄片| 亚洲国产成人一精品久久久| 欧美97在线视频| 蜜桃在线观看..| 男女无遮挡免费网站观看| 狠狠婷婷综合久久久久久88av| 老鸭窝网址在线观看| 国产精品亚洲av一区麻豆| 亚洲精品国产一区二区精华液| 亚洲欧洲日产国产| 高潮久久久久久久久久久不卡| 热re99久久国产66热| 国产老妇伦熟女老妇高清| 久久久精品区二区三区| 日韩,欧美,国产一区二区三区| 中文字幕人妻丝袜制服| 2018国产大陆天天弄谢| 久久国产精品人妻蜜桃| 热re99久久精品国产66热6| 精品亚洲乱码少妇综合久久| 欧美激情 高清一区二区三区| 一级,二级,三级黄色视频| 成人黄色视频免费在线看| 动漫黄色视频在线观看| 一个人免费看片子| 一二三四社区在线视频社区8| av有码第一页| 又紧又爽又黄一区二区| 国产成人欧美在线观看 | 亚洲国产精品一区二区三区在线| 久久久精品国产亚洲av高清涩受| 国产成人免费无遮挡视频| 欧美国产精品va在线观看不卡| 婷婷丁香在线五月| 欧美另类亚洲清纯唯美| 国产在视频线精品| 美女视频免费永久观看网站| 欧美激情 高清一区二区三区| 男女免费视频国产| 在线av久久热| 黑人操中国人逼视频| 久久久久久人人人人人| 亚洲国产日韩一区二区| 免费看十八禁软件| 国产精品欧美亚洲77777| av不卡在线播放| 伊人久久大香线蕉亚洲五| 91麻豆精品激情在线观看国产 | √禁漫天堂资源中文www| 亚洲欧美一区二区三区久久| www.自偷自拍.com| 人人妻人人澡人人爽人人夜夜| 色精品久久人妻99蜜桃| 高清av免费在线| 国产精品一区二区精品视频观看| 国产精品香港三级国产av潘金莲| 久久久国产精品麻豆| 欧美成狂野欧美在线观看| av天堂久久9| 不卡一级毛片| 久久久久精品人妻al黑| 一本—道久久a久久精品蜜桃钙片| 看免费av毛片| 蜜桃国产av成人99| 日韩三级视频一区二区三区| 女警被强在线播放| 国产亚洲精品一区二区www | 我要看黄色一级片免费的| 妹子高潮喷水视频| 99精国产麻豆久久婷婷| 亚洲精品乱久久久久久| 亚洲中文字幕日韩| 国产91精品成人一区二区三区 | 国产一区二区三区av在线| 一级毛片女人18水好多| 不卡一级毛片| 国产在线视频一区二区| av又黄又爽大尺度在线免费看| 美女国产高潮福利片在线看| 日韩电影二区| 69av精品久久久久久 | 精品少妇内射三级| 男人舔女人的私密视频| 久久久久精品人妻al黑| 成人黄色视频免费在线看| 久久久久视频综合| 国产精品久久久久成人av| 亚洲 国产 在线| 一二三四社区在线视频社区8| 视频在线观看一区二区三区| 男女高潮啪啪啪动态图| 99久久人妻综合| av有码第一页| 国产麻豆69| 啦啦啦视频在线资源免费观看| 美女视频免费永久观看网站| 国产欧美亚洲国产| 国产极品粉嫩免费观看在线| 欧美日韩精品网址| 精品一品国产午夜福利视频| 中国美女看黄片| 精品熟女少妇八av免费久了| 亚洲精品av麻豆狂野| 老熟妇仑乱视频hdxx| 亚洲人成77777在线视频| 国产成人一区二区三区免费视频网站| 日日夜夜操网爽| 亚洲av成人不卡在线观看播放网 | 成年女人毛片免费观看观看9 | 亚洲激情五月婷婷啪啪| 亚洲欧洲精品一区二区精品久久久| 国产一级毛片在线| 成人亚洲精品一区在线观看| 性少妇av在线| 99国产精品一区二区三区| 男人添女人高潮全过程视频| 国产成人影院久久av| 黄色a级毛片大全视频| 精品亚洲成a人片在线观看| 大陆偷拍与自拍| 老司机福利观看| 亚洲国产精品成人久久小说| 一本一本久久a久久精品综合妖精| 麻豆av在线久日| 亚洲av成人一区二区三| 欧美精品高潮呻吟av久久| cao死你这个sao货| 少妇猛男粗大的猛烈进出视频| 黑人巨大精品欧美一区二区蜜桃| 狠狠狠狠99中文字幕| 精品一区二区三区av网在线观看 | 女性被躁到高潮视频| 在线亚洲精品国产二区图片欧美| 丝袜美腿诱惑在线| 亚洲欧洲日产国产| 亚洲成人免费电影在线观看| 大型av网站在线播放| 国产精品欧美亚洲77777| 韩国高清视频一区二区三区| 久久久国产精品麻豆| 老司机午夜十八禁免费视频| 亚洲国产欧美网| 欧美老熟妇乱子伦牲交| 黄色片一级片一级黄色片| netflix在线观看网站| 精品亚洲成国产av| 国产福利在线免费观看视频| 国产成人欧美在线观看 | 亚洲伊人久久精品综合| av不卡在线播放| 19禁男女啪啪无遮挡网站| 视频区图区小说| 一本久久精品| 国产欧美日韩一区二区三 | 自拍欧美九色日韩亚洲蝌蚪91| 黑人操中国人逼视频| 午夜福利,免费看| 日韩有码中文字幕| 精品卡一卡二卡四卡免费| 丁香六月天网| 午夜福利在线观看吧| 免费久久久久久久精品成人欧美视频| 美女视频免费永久观看网站| 欧美激情高清一区二区三区| 成年女人毛片免费观看观看9 | 桃花免费在线播放| 亚洲一区中文字幕在线| 99re6热这里在线精品视频| 亚洲激情五月婷婷啪啪| 91精品国产国语对白视频| 日韩欧美一区二区三区在线观看 | 色老头精品视频在线观看| h视频一区二区三区| 国产成人影院久久av| www.熟女人妻精品国产| 无遮挡黄片免费观看| 国产精品久久久久久人妻精品电影 | 69av精品久久久久久 | 少妇的丰满在线观看| 午夜激情av网站| 男女床上黄色一级片免费看| 十分钟在线观看高清视频www| 精品福利观看| av视频免费观看在线观看| 午夜老司机福利片| 国产精品1区2区在线观看. | 免费一级毛片在线播放高清视频 | 欧美 日韩 精品 国产| 国产一区二区激情短视频 | 12—13女人毛片做爰片一| 免费看十八禁软件| 狂野欧美激情性bbbbbb| 交换朋友夫妻互换小说| 免费在线观看视频国产中文字幕亚洲 | 青草久久国产| 王馨瑶露胸无遮挡在线观看| 国产野战对白在线观看| 亚洲精品日韩在线中文字幕| 亚洲自偷自拍图片 自拍| 十八禁高潮呻吟视频| 每晚都被弄得嗷嗷叫到高潮| 精品一区二区三区av网在线观看 | 99热网站在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲少妇的诱惑av| 91av网站免费观看| 成年人黄色毛片网站| 看免费av毛片| 精品少妇一区二区三区视频日本电影| 日本91视频免费播放| 亚洲精品美女久久久久99蜜臀| 肉色欧美久久久久久久蜜桃| 在线观看免费午夜福利视频| xxxhd国产人妻xxx| 色播在线永久视频| 欧美老熟妇乱子伦牲交| 老司机午夜十八禁免费视频| 99国产极品粉嫩在线观看| 亚洲av欧美aⅴ国产| 超碰成人久久| 性高湖久久久久久久久免费观看| 别揉我奶头~嗯~啊~动态视频 | 狠狠狠狠99中文字幕| 欧美大码av| 日韩精品免费视频一区二区三区| 亚洲一区二区三区欧美精品| 欧美乱码精品一区二区三区| 久久狼人影院| 高清av免费在线| 国产高清视频在线播放一区 | 水蜜桃什么品种好| 狠狠婷婷综合久久久久久88av| 正在播放国产对白刺激| 国产av一区二区精品久久| 久9热在线精品视频| 国产精品国产三级国产专区5o| 久久精品人人爽人人爽视色| 日日夜夜操网爽| 久久精品熟女亚洲av麻豆精品| 视频在线观看一区二区三区| 少妇 在线观看| 精品人妻熟女毛片av久久网站| 亚洲av国产av综合av卡| 欧美人与性动交α欧美软件| 少妇被粗大的猛进出69影院| 母亲3免费完整高清在线观看| av在线播放精品| 国产淫语在线视频| 各种免费的搞黄视频| 久久久精品区二区三区| 男女下面插进去视频免费观看| 无限看片的www在线观看| 国产91精品成人一区二区三区 | 日韩制服骚丝袜av| 久久国产亚洲av麻豆专区| 涩涩av久久男人的天堂| 欧美日韩成人在线一区二区| 欧美精品亚洲一区二区| 日韩熟女老妇一区二区性免费视频| 高潮久久久久久久久久久不卡| 精品国产一区二区三区四区第35| 国产国语露脸激情在线看| 少妇粗大呻吟视频| 电影成人av| 少妇猛男粗大的猛烈进出视频| 亚洲第一av免费看| 精品国产一区二区久久| 91精品伊人久久大香线蕉| 18禁观看日本| 丁香六月欧美| 婷婷色av中文字幕| 成年人免费黄色播放视频| 老司机午夜十八禁免费视频| 久久精品亚洲熟妇少妇任你| 国产av一区二区精品久久| 欧美亚洲 丝袜 人妻 在线| www.熟女人妻精品国产| 午夜激情av网站| 久久亚洲国产成人精品v| 午夜福利在线观看吧| 人妻人人澡人人爽人人| 日韩电影二区| 久久精品成人免费网站| 黑人欧美特级aaaaaa片| 99热全是精品| 91成年电影在线观看| 亚洲av男天堂| 精品一区二区三区av网在线观看 | 热99久久久久精品小说推荐| 精品久久久久久久毛片微露脸 | 午夜福利在线观看吧| 一进一出抽搐动态| 丁香六月欧美| 人妻 亚洲 视频| 啪啪无遮挡十八禁网站| 欧美少妇被猛烈插入视频| 亚洲精品成人av观看孕妇| 男女床上黄色一级片免费看| 日韩有码中文字幕| 国产精品一区二区精品视频观看| 国产不卡av网站在线观看| a级毛片在线看网站| 欧美精品av麻豆av| 狂野欧美激情性bbbbbb| 一二三四社区在线视频社区8| 亚洲人成电影免费在线| 一本大道久久a久久精品| 老熟妇仑乱视频hdxx| 老司机靠b影院| 国产精品.久久久| 最近最新免费中文字幕在线| 国产激情久久老熟女| 下体分泌物呈黄色| 我要看黄色一级片免费的| 亚洲国产欧美日韩在线播放| 他把我摸到了高潮在线观看 | 精品久久久久久电影网| 国产成人免费无遮挡视频| 自线自在国产av| 国产一区二区三区在线臀色熟女 | 人妻久久中文字幕网| 精品国内亚洲2022精品成人 | av线在线观看网站| 国产精品九九99| 亚洲色图 男人天堂 中文字幕| 久久久久精品国产欧美久久久 | 亚洲中文字幕日韩| 免费观看av网站的网址| 午夜久久久在线观看| 欧美日韩精品网址| 欧美黑人欧美精品刺激| 十分钟在线观看高清视频www| 国产精品久久久久久精品古装| 午夜福利免费观看在线| bbb黄色大片| 欧美 亚洲 国产 日韩一| 91麻豆精品激情在线观看国产 | 天天添夜夜摸| 中文字幕精品免费在线观看视频| 永久免费av网站大全| 黄片大片在线免费观看| 黑人猛操日本美女一级片| 成年av动漫网址| 制服诱惑二区| 97人妻天天添夜夜摸| 91av网站免费观看| 午夜日韩欧美国产| 国产一区二区激情短视频 | 男人舔女人的私密视频| 中文欧美无线码| 精品少妇一区二区三区视频日本电影| 天堂中文最新版在线下载| 日日摸夜夜添夜夜添小说| 国产精品久久久av美女十八| 一级,二级,三级黄色视频| 免费在线观看视频国产中文字幕亚洲 | 天天躁狠狠躁夜夜躁狠狠躁| 极品人妻少妇av视频| 啦啦啦免费观看视频1| 国产97色在线日韩免费| 久久人人97超碰香蕉20202| 男人操女人黄网站| 麻豆乱淫一区二区| 亚洲黑人精品在线| 免费在线观看视频国产中文字幕亚洲 | 日韩视频一区二区在线观看| 老司机深夜福利视频在线观看 | 免费在线观看日本一区| 亚洲第一青青草原| 欧美另类一区| 97在线人人人人妻| 1024视频免费在线观看| 97人妻天天添夜夜摸| 最近最新中文字幕大全免费视频| 亚洲欧美一区二区三区久久| 日韩三级视频一区二区三区| 亚洲av电影在线观看一区二区三区| 久久天堂一区二区三区四区| 国产精品欧美亚洲77777| 91精品三级在线观看| 正在播放国产对白刺激| 精品一品国产午夜福利视频| 日韩中文字幕视频在线看片| 一个人免费看片子| 中文字幕人妻熟女乱码| 精品人妻一区二区三区麻豆| 99久久精品国产亚洲精品| 国产欧美日韩一区二区三区在线| 免费观看a级毛片全部| 亚洲精品成人av观看孕妇| 美女主播在线视频| 女人被躁到高潮嗷嗷叫费观| 亚洲一卡2卡3卡4卡5卡精品中文| 国产日韩欧美视频二区| 一级片免费观看大全| 欧美人与性动交α欧美软件| 精品国产乱码久久久久久男人| 国产成人精品久久二区二区免费| 国产精品.久久久| 91精品国产国语对白视频| 欧美黑人欧美精品刺激| 脱女人内裤的视频| 亚洲国产精品999| 国产1区2区3区精品| 精品国产乱码久久久久久男人| 一级a爱视频在线免费观看| 亚洲av国产av综合av卡| 好男人电影高清在线观看| 久9热在线精品视频| 不卡一级毛片| 一本色道久久久久久精品综合| 精品人妻一区二区三区麻豆| 桃花免费在线播放| 97在线人人人人妻| 久久久久精品国产欧美久久久 | 一区在线观看完整版| av免费在线观看网站| 天天操日日干夜夜撸| 午夜福利影视在线免费观看| 日韩欧美免费精品| 国产精品免费视频内射| 这个男人来自地球电影免费观看| 女性被躁到高潮视频| 69av精品久久久久久 | 97精品久久久久久久久久精品| 亚洲精品中文字幕一二三四区 | 国产精品秋霞免费鲁丝片| 国产亚洲av片在线观看秒播厂| 国产成人av教育| 一区二区日韩欧美中文字幕| 精品国内亚洲2022精品成人 | 不卡av一区二区三区| 91av网站免费观看| 免费日韩欧美在线观看| 伊人亚洲综合成人网| 超碰成人久久| 亚洲av日韩在线播放| 国产精品久久久久久精品电影小说| 成人亚洲精品一区在线观看| 久热爱精品视频在线9| 18在线观看网站| 亚洲少妇的诱惑av| 久久久久国产一级毛片高清牌| 搡老岳熟女国产| 美女扒开内裤让男人捅视频| 蜜桃在线观看..| 免费在线观看日本一区| 80岁老熟妇乱子伦牲交| 午夜免费鲁丝| 国产在线视频一区二区| 一边摸一边抽搐一进一出视频| 老鸭窝网址在线观看| 免费在线观看视频国产中文字幕亚洲 | 成人18禁高潮啪啪吃奶动态图| 午夜福利在线免费观看网站| 亚洲中文日韩欧美视频| 亚洲精华国产精华精| 国产精品影院久久| 菩萨蛮人人尽说江南好唐韦庄| 日本91视频免费播放| 2018国产大陆天天弄谢| 国产亚洲欧美在线一区二区| 午夜91福利影院| 在线观看免费日韩欧美大片| 老司机靠b影院| 后天国语完整版免费观看| 国产成人系列免费观看| 日本撒尿小便嘘嘘汇集6| 国产深夜福利视频在线观看| 久久久久久亚洲精品国产蜜桃av| 免费少妇av软件| 中文字幕av电影在线播放| 国产亚洲欧美在线一区二区| 色播在线永久视频| 嫁个100分男人电影在线观看| 老汉色∧v一级毛片| 国产97色在线日韩免费| 久久久精品免费免费高清| 18禁裸乳无遮挡动漫免费视频| 一二三四社区在线视频社区8| 日日摸夜夜添夜夜添小说| 满18在线观看网站| 在线观看www视频免费| 麻豆国产av国片精品| 一区二区日韩欧美中文字幕| 亚洲情色 制服丝袜| 中文字幕另类日韩欧美亚洲嫩草| 动漫黄色视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产在线免费精品| 91九色精品人成在线观看| 久久久久国内视频| 亚洲精品成人av观看孕妇| 午夜精品久久久久久毛片777| 91字幕亚洲| 日韩中文字幕视频在线看片| 日韩精品免费视频一区二区三区| 99久久99久久久精品蜜桃| 成年美女黄网站色视频大全免费| 午夜精品国产一区二区电影| 免费不卡黄色视频| 人妻一区二区av| 啦啦啦中文免费视频观看日本| 中亚洲国语对白在线视频| 日韩一区二区三区影片| 一本—道久久a久久精品蜜桃钙片| 日韩欧美一区视频在线观看| 久久久久久久久免费视频了| 精品福利永久在线观看| 黑人巨大精品欧美一区二区mp4| 国产成人影院久久av| 日本五十路高清| 91老司机精品| 纵有疾风起免费观看全集完整版| av天堂在线播放| 91老司机精品| 欧美激情高清一区二区三区| 三上悠亚av全集在线观看| 中文字幕色久视频| www.av在线官网国产| 中文字幕av电影在线播放| 人人妻人人爽人人添夜夜欢视频|