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

    基于顆粒流程序的廣義Kelvin 模型及其應(yīng)用

    2015-02-04 12:17:48金愛兵張秀鳳孟新秋楊振偉
    巖土力學(xué) 2015年9期
    關(guān)鍵詞:模型

    金愛兵,王 凱,張秀鳳,孟新秋,楊振偉

    (1.北京科技大學(xué) 土木與環(huán)境工程學(xué)院,北京 100083;2.溫州市城建設(shè)計(jì)院,浙江 溫州 325000)

    1 引 言

    土體作為工程材料的一種,在外荷載作用下的力學(xué)響應(yīng)不能完全由經(jīng)典力學(xué)描述。特別是軟土,由于其內(nèi)部結(jié)構(gòu)特殊,其應(yīng)力-應(yīng)變關(guān)系呈現(xiàn)出一定的非線性流變特征。陳宗基[1]通過將Maxwell 體與Hook 彈性體并聯(lián),得到一種新模型來描述土體變形規(guī)律;孫鈞[2]對(duì)上海地區(qū)3 種典型軟土進(jìn)行了大量的單軸剪切與三軸壓縮流變?cè)囼?yàn),并對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行詳細(xì)分析后,提出了上海軟土蠕變特性的方程式。粉質(zhì)黏土作為一種常見軟土,也常常表現(xiàn)出一定的流變特性。王坤等[3]對(duì)深圳地區(qū)深層粉質(zhì)黏土進(jìn)行了低圍壓作用下的三軸流變?cè)囼?yàn),并得到了粉質(zhì)黏土的三軸流變規(guī)律。而對(duì)于服役時(shí)間較長的露天礦山邊坡,粉質(zhì)黏土的流變效應(yīng)對(duì)邊坡工程的安全影響更為顯著。

    目前,對(duì)巖土流變特性的研究主要集中在室內(nèi)試驗(yàn)[4-6]和流變模型開發(fā)[7-9]等方面。室內(nèi)試驗(yàn)研究主要有單軸試驗(yàn)[6]和不同圍壓下的三軸試驗(yàn)[10-14]。熊良宵等[15]將修正后的Burgers 模型應(yīng)用于蠕變?cè)囼?yàn)的擬合分析;李娜等[16]用改進(jìn)的西原模型對(duì)循環(huán)加、卸載試驗(yàn)結(jié)果進(jìn)行擬合,所得效果良好。但室內(nèi)試驗(yàn)研究往往存在操作復(fù)雜、花費(fèi)時(shí)間長、試驗(yàn)費(fèi)用高等缺點(diǎn)。

    作為巖土工程穩(wěn)定性研究的常用工具,數(shù)值模擬方法已經(jīng)被許多科研工作者應(yīng)用于巖土流變特性的研究。張志沛等[17]利用FLAC3D中的Burgers 模型對(duì)不同形狀的試樣進(jìn)行了蠕變特性模擬分析;熊良宵等[18]在FLAC3D中開發(fā)出新的流變模型,并將其應(yīng)用到巖土應(yīng)變軟化的相關(guān)工程計(jì)算中;劉珊珊等[19]在FLAC3D中二次開發(fā)出黏彈性廣義Kelvin 模型,并用算例驗(yàn)證了所開發(fā)模型的正確性;王濤等[20]開發(fā)出適用于PFC2D的廣義Kelvin 接觸本構(gòu)模型,并將其應(yīng)用到了隧洞工程問題中。

    本文在對(duì)PFC2D中自定義本構(gòu)模型進(jìn)行研究的基礎(chǔ)上,利用Microsoft Visual Studio 平臺(tái),開發(fā)出廣義Kelvin 模型的動(dòng)態(tài)鏈接庫(DLL)文件。通過算例驗(yàn)證程序編制的準(zhǔn)確性,并在此基礎(chǔ)上,將其應(yīng)用到露天煤礦含粉質(zhì)黏土邊坡的失穩(wěn)研究中。

    2 廣義Kelvin 模型本構(gòu)關(guān)系

    廣義Kelvin 模型由一個(gè)Maxwell 體和任意個(gè)Kelvin 體串聯(lián)組成,其模型如圖1 所示。其中,Maxwell 體由一個(gè)彈簧和一個(gè)阻尼器串聯(lián)而成,Kelvin 體由一個(gè)彈簧和一個(gè)阻尼器并聯(lián)而成,Kn和Ks代表彈簧法向和切向剛度,Cn和Cs代表阻尼器法向和切向阻尼,下標(biāo)數(shù)字1,2…,n為原件編號(hào)。

    圖1 廣義Kelvin 模型Fig.1 Generalized Kelvin model

    由圖1 可知,廣義Kelvin 模型的總位移u 等于一個(gè)Maxwell 體的位移 um1與若干個(gè)Kelvin 體的位移ukr(下標(biāo)r為開爾文體的編號(hào),r=2,3,…,n)之和,如式(1)所示。其中,um1由Maxwell 體中彈簧的位移 uk1和阻尼器的位移 uc1共同組成。

    對(duì)廣義Kelvin 模型的本構(gòu)關(guān)系推導(dǎo)如下,下列各式中的“±”(或“?”)分別對(duì)應(yīng)著法向和切向的情況;Kr和Cr(r=1,2,…,n) 分別為廣義Kelvin 模型中彈簧的剛度和阻尼器的阻尼;f為模型中顆粒間的接觸力;Ar、Br、C和D 是分別與Kr、Cr和Δt 有關(guān)的中間變量。

    對(duì)式(1)進(jìn)行一階求導(dǎo),得

    接觸力f 可表示為

    將式(3)、(5)、(6)代入式(2)得

    利用有限差分格式的中心差分法近似對(duì)時(shí)間t進(jìn)行差分,其中對(duì)t 取平均值得

    對(duì)于Kelvin 體,則有

    整理后,接觸力ft+1可表示為

    ft+1可由已知量ut+1、ut、和ft求出。

    3 廣義Kelvin 模型開發(fā)

    PFC 程序中的用戶自定義接觸模型的編寫包括基類描述、模型注冊(cè)、成員函數(shù)描述、模型和PFC程序間的信息交換。軟件還提供了自帶本構(gòu)模型的源代碼,用戶在進(jìn)行本構(gòu)模型開發(fā)時(shí),可以利用和修改這些代碼。下面介紹如何修改本構(gòu)模型的頭文件和源文件,從而創(chuàng)建用戶自己的本構(gòu)模型。

    在頭文件中進(jìn)行新的本構(gòu)模型派生類的聲明,修改模型的編號(hào)、名稱和版本,修改派生類的私有成員,包括模型的基本變量及主要中間變量。為避免沖突,新自定義模型的編號(hào)需不小于100。

    在源文件中修改以下各個(gè)成員函數(shù):

    (1)修改構(gòu)造函數(shù)Kelvin::Kelvin(),對(duì)其私有成員進(jìn)行初始化。

    (2)修改Kelvin::Name()函數(shù),返回用戶自定義模型的名稱。

    (3)根據(jù)派生類私有成員的名稱字符串,依次修改Kelvin::PropNames()、Kelvin::ReturnProp()和Kelvin::AcceptProp()3 個(gè)函數(shù)。PFC 的PROPERTY命令可用來調(diào)用Kelvin::AcceptProp()函數(shù),對(duì)模型的參數(shù)進(jìn)行賦值。

    (4)修改Kelvin::FDlaw()函數(shù),此函數(shù)在計(jì)算離散元的力與位移關(guān)系時(shí)被調(diào)用。用戶根據(jù)推導(dǎo)的模型中力與位移數(shù)學(xué)關(guān)系,編寫C++語句,由應(yīng)變?cè)隽坑?jì)算得到應(yīng)力增量,從而獲得新的應(yīng)力。本構(gòu)模型將在每個(gè)循環(huán)中更新顆粒間接觸的力與力矩,以及法向和切向剛度。

    (5)修改USERCM1::SaveRestore()函數(shù),方法同步驟(3),此函數(shù)用來保存與恢復(fù)實(shí)型、整型和邏輯型的內(nèi)部數(shù)據(jù),PFC 程序給出save和restore命令時(shí)調(diào)用此函數(shù)。

    另外,針對(duì)本文自定義的廣義Kelvin本構(gòu)模型,程序中設(shè)置控制變量i來靈活控制Kelvin體的個(gè)數(shù)。編寫Kelvin.h和Kelvin.cpp 文件后,編譯動(dòng)態(tài)鏈接庫文件。至此,所開發(fā)的廣義Kelvin 模型以動(dòng)態(tài)鏈接庫文件的形式存儲(chǔ)起來,可供PFC2D程序調(diào)用。

    4 廣義Kelvin 模型驗(yàn)證

    當(dāng)控制變量i=0 時(shí),即廣義Kelvin 模型中只含有Maxwell 體,無Kelvin 體時(shí),模型即為Maxwell模型;當(dāng)控制變量i=1 時(shí),即廣義Kelvin 模型中只含有Maxwell 體和1 個(gè)Kelvin 體時(shí),模型即為Burgers 模型。本文分別采用i=0和1 兩種情況下的Maxwell 模型和Burgers 模型進(jìn)行松弛試驗(yàn)的驗(yàn)證。

    (1)Maxwell 模型理論值求解

    當(dāng)i=0,無Kelvin 體時(shí),式(7)可表示為

    (2)Burgers 模型理論值求解

    當(dāng)i=1 時(shí),對(duì)于單位輸入函數(shù) u( t)=1,Burgers模型的理論值為

    兩種模型中Maxwell 體彈簧法向、切向剛度:K1n=K1s=100MPa;阻尼器法向、切向阻尼:C1n=C1s=100MPa·s;Kelvin體彈簧法向、切向剛度:K2n=K2s=100MPa,阻尼器法向、切向阻尼:C2n=C2s=100MPa·s。

    通過上述計(jì)算,得到兩種接觸模型松弛過程的理論值后,采用PFC2D命令流和Fish 語言編寫應(yīng)力松弛程序,然后在PFC2D軟件中調(diào)用該程序,分別加載i為0和1 時(shí)的廣義Kelvin 本構(gòu)模型動(dòng)態(tài)鏈接庫文件,模擬兩個(gè)相互接觸球體(兩球重疊量為0.01 m)的應(yīng)力松弛過程,得出球體間法向接觸力隨時(shí)間的變化關(guān)系,并將數(shù)值模擬解與理論值進(jìn)行對(duì)比分析。其對(duì)比結(jié)果如圖2 所示。

    由圖2 可知,兩球接觸力隨時(shí)間的增加呈指數(shù)趨勢(shì)減小,本構(gòu)模型的數(shù)值模擬解與理論解所對(duì)應(yīng)曲線相互吻合,從而驗(yàn)證了所開發(fā)模型的準(zhǔn)確性。

    圖2 法向接觸力模擬計(jì)算結(jié)果與理論結(jié)果Fig.2 Analytical and numerical results of normal contact forces

    5 工程應(yīng)用

    5.1 邊坡失穩(wěn)規(guī)律

    含軟弱粉質(zhì)黏土邊坡廣泛分布于山西、陜西等黃土高原地區(qū),尤其是在山西地區(qū),煤炭資源豐富,露天開采形成的含粉質(zhì)黏土邊坡普遍存在。在山西某露天煤礦,經(jīng)長期實(shí)地調(diào)研,發(fā)現(xiàn)此類邊坡的滑坡形式較為特殊,其失穩(wěn)規(guī)律有以下幾個(gè)特征:①邊坡的滑坡區(qū)域含有粉質(zhì)黏土層;②巖土接觸帶傾向坡內(nèi),傾角較小(6°~9°),且有地下水滲出;③滑坡體前緣平緩,局部有底鼓隆起現(xiàn)象;④滑坡體中部較為平緩;⑤滑坡體后緣較陡,張裂縫嚴(yán)重,下落量??;⑥ 滑動(dòng)面形式表現(xiàn)異常,為“圓弧+折線”的組合型,最終輪廓表現(xiàn)為帽檐形。

    此類含粉質(zhì)黏土邊坡失穩(wěn)破壞的一般過程與現(xiàn)場情況如圖3 所示。

    5.2 模型計(jì)算參數(shù)

    巖土體的物理力學(xué)參數(shù)由巖石壓縮試驗(yàn)、土工直剪試驗(yàn)及粉質(zhì)黏土三軸流變?cè)囼?yàn)獲得,見表1。參照室內(nèi)試驗(yàn),進(jìn)行數(shù)值模擬試驗(yàn),粉質(zhì)黏土部分調(diào)用本文自主開發(fā)的廣義Kelvin 模型,其他巖土體均采用PFC2D軟件默認(rèn)線彈性接觸模型;數(shù)值模型尺寸均為100 mm×50 mm 的標(biāo)準(zhǔn)試樣;巖土體顆粒間均采用能夠抗拉、抗剪、承受彎矩的平行黏結(jié)。

    模擬試驗(yàn)過程中,對(duì)廣義Kelvin 模型中的控制變量i 及模型的細(xì)觀力學(xué)參數(shù)進(jìn)行反復(fù)調(diào)試。最終結(jié)果表明,當(dāng)控制變量i=2,即廣義Kelvin 模型中含1 個(gè)Maxwell 體和2 個(gè)Kelvin 體時(shí),粉質(zhì)黏土的模擬試驗(yàn)與室內(nèi)試驗(yàn)結(jié)果較吻合(見圖4、5)。此時(shí),在廣義Kelvin 模型中,Maxwell 體彈簧法向、切向剛度:K1n=K1s=100 MPa;阻尼器法向、切向阻尼:C1n=C1s=300 MPa·s;兩個(gè)Kelvin 體相同,彈簧法向、切向剛度:K2n=K2s=K3n=K3s=100 MPa;阻尼器法向、切向阻尼:C2n=C2s=C3n=C3s=300 MPa·s 。邊坡模型中巖土體的其他細(xì)觀力學(xué)參數(shù)見表2。

    圖3 邊坡失穩(wěn)破壞過程與現(xiàn)場Fig.3 Process and scene of slope failure

    表1 巖土體力學(xué)參數(shù)Table 1 Mechanical parameters of rock and soils

    圖4 數(shù)值試驗(yàn)與直剪試驗(yàn)結(jié)果Fig.4 Results of numerical test and direct shear test

    圖5 數(shù)值試驗(yàn)與流變?cè)囼?yàn)結(jié)果Fig.5 Results of numerical test and rheological test

    5.3 邊坡模型構(gòu)建

    根據(jù)參數(shù)調(diào)試過程中獲取的廣義Kelvin模型和細(xì)觀力學(xué)參數(shù),建立邊坡開挖前的PFC2D數(shù)值計(jì)算模型,見圖6。根據(jù)工程實(shí)際情況,邊坡模型長×高尺寸為:217.51 m×82.66 m。模型內(nèi)部共生成36 261 個(gè)巖土體顆粒,模型的外部由1 道水平墻體和2 道豎直墻體構(gòu)成邊界限制條件,墻體法向剛度取 kn=1 000 MN/m,切向剛度取 ks=1 000 MN/m,摩擦系數(shù)取1.0。模型中不同顏色顆粒分別代表不同地層,邊坡自上而下依次為粉土、粉質(zhì)黏土(天然)、粉質(zhì)黏土(飽和)、接觸帶及泥巖。

    表2 顆粒細(xì)觀力學(xué)參數(shù)Table 2 Meso mechanical parameters of particles

    圖6 邊坡數(shù)值模型Fig.6 Numerical model of slope

    對(duì)于粉質(zhì)黏土,考慮到吸水-失水循環(huán)對(duì)其強(qiáng)度的影響,根據(jù)現(xiàn)場實(shí)際情況,在邊坡模型的不同地層位置和計(jì)算的不同時(shí)段,對(duì)粉質(zhì)黏土分別采用天然強(qiáng)度與飽和強(qiáng)度。無降水時(shí)下層粉質(zhì)黏土采用天然強(qiáng)度,降水時(shí)下層粉質(zhì)黏土采用飽和強(qiáng)度。后續(xù)數(shù)值模擬中,邊坡開挖后通過在計(jì)算的不用時(shí)步對(duì)該層粉質(zhì)黏土賦予相應(yīng)參數(shù)來實(shí)現(xiàn)這一過程,根據(jù)當(dāng)?shù)赜昙镜慕邓?guī)律:模擬中每40 萬時(shí)步采用天然強(qiáng)度,每10 萬時(shí)步采用飽和強(qiáng)度,循環(huán)多次至邊坡最終破壞形成,共計(jì)約240 萬時(shí)步。

    5.4 計(jì)算結(jié)果及分析

    在PFC2D中編寫命令流,實(shí)現(xiàn)對(duì)邊坡數(shù)值模型的計(jì)算與開挖。對(duì)于邊坡中的粉質(zhì)黏土地層,分別采用線彈性接觸模型和本文所開發(fā)的具有流變特性的廣義Kelvin 模型兩種情況進(jìn)行模擬計(jì)算,兩算例中吸水-失水循環(huán)模擬過程相同。最終將其計(jì)算結(jié)果進(jìn)行對(duì)比分析。

    計(jì)算中采用顆粒間平行黏結(jié)的拉裂或剪斷來表示邊坡巖土體結(jié)構(gòu)的破壞[21],為便于觀察,將模型中黏結(jié)破壞的顆粒設(shè)置為深藍(lán)色。模型中粉質(zhì)黏土地層采用線彈性接觸模型的計(jì)算結(jié)果如圖7 所示,邊坡最終破壞模式呈常見的圓弧形土體滑坡模式,但與現(xiàn)場邊坡實(shí)際破壞情況不符。

    圖7 線彈性接觸模型邊坡破壞結(jié)果Fig.7 Slope failure predicted by a linear elastic contact model

    模型中粉質(zhì)黏土地層采用用戶自定義廣義Kelvin 模型的計(jì)算結(jié)果,如圖8 所示,邊坡最終破壞模式比較符合現(xiàn)場的實(shí)際情況。計(jì)算過程中,邊坡的失穩(wěn)破壞可分為5 個(gè)階段。

    第1階段:邊坡施工開挖,如圖8(a)所示,各臺(tái)階開挖臨空面小范圍內(nèi)巖土體受到擾動(dòng),深藍(lán)色顆粒開始出現(xiàn),并零星散布于邊坡面的淺層。

    第2階段:邊坡局部破壞,如圖8(b)所示,邊坡內(nèi)部整體變化微小,因賦存地下水,粉質(zhì)黏土層強(qiáng)度降低,其坡腳處開始有深藍(lán)色顆粒聚集,局部位移增大,破壞初顯。

    第3階段:破壞范圍擴(kuò)大,如圖8(c)所示,飽和粉質(zhì)黏土層邊坡坡腳處深藍(lán)色顆粒增多,孔隙率增加,破壞范圍擴(kuò)大,少量黏結(jié)已被破壞的顆粒在重力作用下開始滑移出坡面,邊坡后緣開始呈現(xiàn)出較明顯的圓弧滑面。

    第4階段:邊坡出現(xiàn)失穩(wěn),如圖8(d)所示,坡面周邊大量深藍(lán)色顆粒翻滾下滑,孔隙率大幅增加,邊坡破壞區(qū)域輪廓成“折線+圓弧”的帽檐形。

    第5階段:帽檐形滑坡形成,如圖8(e)所示,模型中的顆粒停止下滑,邊坡模型已經(jīng)破壞。

    圖8 廣義Kelvin 模型邊坡破壞過程Fig.8 Slope failure process derived from the generalized Kelvin model

    通過上述計(jì)算,并結(jié)合現(xiàn)場滑坡情況,邊坡失穩(wěn)過程主要表現(xiàn)為:

    (1)天然狀態(tài)下,粉質(zhì)黏土層強(qiáng)度較高,所在邊坡具有足夠的自穩(wěn)能力,隨著開采的進(jìn)行,粉質(zhì)黏土層中顆粒黏結(jié)遭到破壞,導(dǎo)致裂隙產(chǎn)生,計(jì)算過程中坡體共產(chǎn)生498 條裂隙。

    (2)微觀裂隙貫通產(chǎn)生宏觀裂縫,受水影響,飽和粉質(zhì)黏土在接觸帶處強(qiáng)度較低,顆粒黏結(jié)容易遭到破壞,所以該處產(chǎn)生裂隙最多,成為邊坡破壞的主要區(qū)域。

    (3)受開采影響,失水-吸水循環(huán)過程中,裂隙進(jìn)一步擴(kuò)張,松散體影響范圍逐漸擴(kuò)大,并在邊坡后緣開始形成圓弧滑面。

    (4)因松散體自身承載力不夠,內(nèi)部土體顆粒開始不斷下滑,在粉質(zhì)黏土流變作用下,一段時(shí)間后,粉質(zhì)黏土不能承受上部土體的荷載,松散體呈折線形式向上蔓延,聯(lián)通至邊坡后緣圓弧滑面,以“折線+圓弧”的組合滑動(dòng)面形式破壞,破壞區(qū)最終輪廓呈帽檐形。

    含粉質(zhì)黏土邊坡失穩(wěn)模擬結(jié)果表明,該類邊坡失穩(wěn)與普通土質(zhì)邊坡的失穩(wěn)形態(tài)以及失穩(wěn)過程并不一致,粉質(zhì)黏土的流變性對(duì)該類邊坡失穩(wěn)具有重要影響;而水是在粉質(zhì)黏土流變過程中導(dǎo)致邊坡失穩(wěn)的重要原因,因此,實(shí)際工程中可以采取壓坡阻水的方式處治該類邊坡,以保證邊坡穩(wěn)定。

    6 結(jié) 論

    (1)自定義廣義Kelvin 本構(gòu)模型包含1 個(gè)Maxwell 體和若干個(gè)Kelvin 體,實(shí)際應(yīng)用中可以根據(jù)具體情況選擇合適數(shù)量的Kelvin 體,以便更好地模擬巖土體流變特性。

    (2)在PFC2D軟件中調(diào)用基于C++語言開發(fā)的自定義廣義Kelvin 本構(gòu)模型動(dòng)態(tài)鏈接庫,進(jìn)行算例驗(yàn)證,得到的數(shù)值模擬解與理論解吻合,模型具有較好的準(zhǔn)確性。

    (3)利用PFC 調(diào)用自定義廣義Kelvin 本構(gòu)模型對(duì)國內(nèi)某含粉質(zhì)黏土邊坡模擬表明,含1個(gè)Maxwell體和2 個(gè)Kelvin 體的用戶自定義廣義Kelvin 模型可很好地模擬粉質(zhì)黏土的流變特性,得到的邊坡失穩(wěn)破壞過程與現(xiàn)場基本一致。

    (4)由于流變特性的影響,該類粉質(zhì)黏土邊坡破壞形式較為特殊,與常規(guī)土質(zhì)邊坡破壞呈圓弧滑動(dòng)不一致,其破壞區(qū)最終輪廓成“折線+圓弧”的帽檐形。

    [1]陳宗基.固結(jié)及時(shí)間效應(yīng)的單維問題[J].土木工程學(xué)報(bào),1958,5(1):1-10.CHEN Zong-ji.Unidimensional problems of consolidation and time effect[J].China Civil Engineering Journal,1958,5(1):1-10.

    [2]孫鈞.巖土材料流變及其工程應(yīng)用[M].北京:中國建筑工業(yè)出版社,1999.

    [3]王坤,蔡永昌,魏赟,等.深圳地區(qū)深層粉質(zhì)黏土的流變特性試驗(yàn)研究[J].巖土工程界,2009,12(12):86-89.WANG Kun,CAI Yong-chang,WEI Yun,et al.Experimental study on the rheological properties of deep silty clay in Shenzhen[J].Geotechnical Engineering World,2009,12(12):86-89.

    [4]ZHANG ZHILIANG,XU WEIYA,WANG WEI,et al.Triaxial creep tests of rock from the compressive zone of dam foundation in Xiangjiaba Hydropower Station[J].International Journal of Rock Mechanics and Mining Sciences,2012,50(1):133-139.

    [5]胡華,鄭曉栩.動(dòng)載作用頻率對(duì)海相沉積軟土動(dòng)態(tài)流變特性影響試驗(yàn)研究[J].巖土力學(xué),2013,34(增刊1):9-13.HU Hua,ZHENG Xiao-xu.Experimental research on dynamic rheological characteristics of marine deposit soft soil under different frequencies of dynamic loading[J].Rock and Soil Mechanics,2013,34(Supp.1):9-13.

    [6]LI YONGSHENG,XIA CAICHU.Time-dependent tests on intact rocks in uniaxial compression[J].International Journal of Rock Mechanics and Mining Sciences,2000,37(3):467-475.

    [7]ZHANG HUABIN,WANG ZHIYIN,ZHENG YALI,et al.Study on triaxial creep experiment and constitutive relation of different rock salts[J].Safety Science,2012,50(4):801-805.

    [8]LIU LANG,WANG GEMIN,CHEN JIANHONG,et al.Creep experiment and rheological model of deep saturated rock[J].Transactions of Nonferrous Metals Society of China,2013,23(2):478-483.

    [9]SHAO J F,ZHU Q Z,SU K.Modeling of creep in rock materials in terms of material degradation[J].Computers and Geotechnics,2003,30(7):549-555.

    [10]ZHANG YU,XU WEIYA,GU JINJIAN,et al.Triaxial creep tests of weak sandstone from fracture zone of high dam foundation[J].Journal of Central South University,2013,20(9):2528-2536.

    [11]YANG SHENGQI,JIANG YUZHOU.Triaxial mechanical creep behavior of sandstone[J].Mining Science and Technology,2010,20(3):339-349.

    [12]ZHAO YANLIN,CAO PING,WANG WEIJUN,et al.Viscoelasto-plastic rheological experiment under circular increment step load and unload and nonlinear creep model of soft rocks[J].Journal of Central South University of Technology,2009,16(3):484-491.

    [13]李亞麗,于懷昌,劉漢東.三軸壓縮下粉砂質(zhì)泥巖蠕變本構(gòu)模型研究[J].巖土力學(xué),2012,33(7):2035-2040.LI Ya-li,YU Huai-chang,LIU Han-dong.Study of creep constitutive model of silty mudstone under triaxial compression[J].Rock and Soil Mechanics,2012,33(7):2035-2040.

    [14]徐進(jìn),張家生,趙同順,等.軟弱路基土體三軸蠕變?cè)囼?yàn)及蠕變模型研究[J].中南大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,42(10):3136-3142.XU Jin,ZHANG Jia-sheng,ZHAO Tong-shun,et al.Creep triaxial tests and constitutive model of embankment foundation soil[J].Journal of Central South University(Science and Technology),2011,42(10):3136-3142.

    [15]熊良宵,楊林德,張堯.巖石的非定常Burgers 模型[J].中南大學(xué)學(xué)報(bào)(自然科學(xué)版),2010,41(2):679-684.XIONG Liang-xiao,YANG Lin-de,ZHANG Yao.Non-stationary Burgers model for rock[J].Journal of Central South University(Science and Technology),2010,41(2):679-684.

    [16]李娜,曹平,衣永亮,等.分級(jí)加卸載下深部巖石流變實(shí)驗(yàn)及模型[J].中南大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,42(11):3465-3471.LI Na,CAO Ping,YI Yong-liang,et al.Creep properties experiment and model of deep rock with step loading and unloading[J].Journal of Central South University(Science and Technology),2011,42(11):3465-3471.

    [17]張志沛,王芝銀,彭惠.陜南泥巖三軸壓縮蠕變?cè)囼?yàn)及其數(shù)值模擬研究[J].水文地質(zhì)工程地質(zhì),2011,38(1):53-58.ZHANG Zhi-pei,WANG Zhi-yin,PENG Hui.Study on triaxial compression creep test and numerical simulation about mud rock in southern Shaanxi[J].Hydrogeology &Engineering Geology,2011,38(1):53-58.

    [18]熊良宵,楊林德,張堯.硬巖的復(fù)合黏彈塑性流變模型[J].中南大學(xué)學(xué)報(bào)(自然科學(xué)版),2010 41(4):1540-1548.XIONG Liang-xiao,YANG Lin-de,ZHANG Yao.Composite viscoelasto-plastic rheological model for hard rock[J].Journal of Central South University(Science and Technology),2010,41(4):1540-1548.

    [19]劉姍姍,趙同彬.黏彈性廣義Kelvin 模型的FLAC3D二次開發(fā)[J].山東科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2010,29(4):20-23.LIU Shan-shan,ZHAO Tong-bin.Secondary development on generalized viscoelastic Kelvin model with FLAC3D[J].Journal of Shandong University of Science and Technology(Natural Science),2010,29(4):20-23.

    [20]王濤,呂慶,李楊,等.顆粒離散元方法中接觸模型的開發(fā)[J].巖石力學(xué)與工程學(xué)報(bào),2009,28(增刊2):4040-4045.WANG Tao,Lü Qing,LI Yang,et al.Development of contact model in particle discrete element method[J].Chinese Journal of Rock Mechanics and Engineering,2009,28(Supp.2):4040-4045.

    [21]張曉平,吳順川,王思敬.類土質(zhì)路塹邊坡動(dòng)態(tài)監(jiān)測及數(shù)值模擬分析[J].巖石力學(xué)與工程學(xué)報(bào),2008,27(增刊2):3431-3440.ZHANG Xiao-ping,WU Shun-chuan,WANG Si-jing.Dynamic monitoring and numerical analysis of soil-like cut slope[J].Chinese Journal of Rock Mechanics andEngineering,2008,27(Supp.2):3431-3440.

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    国产高清视频在线播放一区| 老熟妇仑乱视频hdxx| 亚洲精品久久午夜乱码| 中出人妻视频一区二区| 99精品久久久久人妻精品| 日本免费a在线| 欧美日韩中文字幕国产精品一区二区三区 | 久久香蕉国产精品| 亚洲美女黄片视频| 大型黄色视频在线免费观看| 国产人伦9x9x在线观看| 亚洲男人的天堂狠狠| 校园春色视频在线观看| 国产麻豆69| 人人妻人人澡人人看| www.www免费av| 国产精品久久视频播放| 亚洲中文av在线| av视频免费观看在线观看| 亚洲国产毛片av蜜桃av| 9热在线视频观看99| 美女福利国产在线| 久久国产亚洲av麻豆专区| 夫妻午夜视频| 亚洲第一欧美日韩一区二区三区| tocl精华| 这个男人来自地球电影免费观看| 丝袜美腿诱惑在线| 国产深夜福利视频在线观看| 三上悠亚av全集在线观看| 精品一区二区三区视频在线观看免费 | 高清欧美精品videossex| 欧美色视频一区免费| av有码第一页| 在线观看一区二区三区激情| 妹子高潮喷水视频| 又黄又粗又硬又大视频| 成人特级黄色片久久久久久久| 中文字幕精品免费在线观看视频| 在线播放国产精品三级| 亚洲一区高清亚洲精品| 亚洲人成电影免费在线| 日韩高清综合在线| 国产精品偷伦视频观看了| 性少妇av在线| 国产高清国产精品国产三级| 两性午夜刺激爽爽歪歪视频在线观看 | 黄色a级毛片大全视频| 国产熟女xx| 一级黄色大片毛片| 国产精品久久久av美女十八| videosex国产| 日韩av在线大香蕉| 久久久精品国产亚洲av高清涩受| 黄色丝袜av网址大全| 国产在线精品亚洲第一网站| 亚洲一区中文字幕在线| 成人国产一区最新在线观看| 男女下面进入的视频免费午夜 | 日韩av在线大香蕉| 国产成人欧美在线观看| 丁香六月欧美| 国产精品99久久99久久久不卡| 国产精品综合久久久久久久免费 | 成在线人永久免费视频| 免费在线观看黄色视频的| 国产亚洲欧美精品永久| 亚洲专区字幕在线| 欧美+亚洲+日韩+国产| 国产高清激情床上av| av有码第一页| 又大又爽又粗| 亚洲一区高清亚洲精品| 色播在线永久视频| 欧美日韩亚洲高清精品| a级片在线免费高清观看视频| 女生性感内裤真人,穿戴方法视频| 一级片免费观看大全| 女警被强在线播放| 夜夜躁狠狠躁天天躁| 免费在线观看视频国产中文字幕亚洲| 黑人巨大精品欧美一区二区蜜桃| 中文字幕人妻丝袜一区二区| 视频在线观看一区二区三区| 操出白浆在线播放| av福利片在线| 国产区一区二久久| 精品国内亚洲2022精品成人| 高潮久久久久久久久久久不卡| 精品高清国产在线一区| 两性午夜刺激爽爽歪歪视频在线观看 | 搡老乐熟女国产| 咕卡用的链子| 热re99久久精品国产66热6| 国产在线观看jvid| 亚洲av日韩精品久久久久久密| av超薄肉色丝袜交足视频| 久久这里只有精品19| 正在播放国产对白刺激| 一边摸一边抽搐一进一出视频| 97碰自拍视频| 国产一区二区三区在线臀色熟女 | 午夜影院日韩av| 丝袜美腿诱惑在线| 中文字幕最新亚洲高清| 午夜精品国产一区二区电影| 91av网站免费观看| svipshipincom国产片| 亚洲精品一卡2卡三卡4卡5卡| 80岁老熟妇乱子伦牲交| 这个男人来自地球电影免费观看| 成人精品一区二区免费| 成人三级黄色视频| 人成视频在线观看免费观看| 黄色a级毛片大全视频| 一进一出抽搐gif免费好疼 | 久久香蕉精品热| 亚洲熟妇中文字幕五十中出 | 最新美女视频免费是黄的| 看片在线看免费视频| 又紧又爽又黄一区二区| 国产高清激情床上av| 中文字幕人妻丝袜制服| 国产极品粉嫩免费观看在线| 精品国内亚洲2022精品成人| 女人被躁到高潮嗷嗷叫费观| 亚洲激情在线av| 天堂中文最新版在线下载| 一级作爱视频免费观看| 欧美最黄视频在线播放免费 | 侵犯人妻中文字幕一二三四区| 精品人妻1区二区| 88av欧美| 日韩人妻精品一区2区三区| av片东京热男人的天堂| 91精品国产国语对白视频| 国产精品亚洲一级av第二区| 午夜a级毛片| 黄色视频,在线免费观看| 一区二区三区精品91| 国产成人av教育| 老司机午夜十八禁免费视频| 美女福利国产在线| 黄色成人免费大全| 美女高潮喷水抽搐中文字幕| 日本撒尿小便嘘嘘汇集6| 免费观看精品视频网站| 精品国产亚洲在线| 91麻豆精品激情在线观看国产 | 最新美女视频免费是黄的| 色综合婷婷激情| 一边摸一边抽搐一进一出视频| 免费看a级黄色片| 一个人免费在线观看的高清视频| 18禁裸乳无遮挡免费网站照片 | 亚洲男人天堂网一区| 一进一出抽搐动态| 丁香欧美五月| av免费在线观看网站| 国产伦一二天堂av在线观看| 色播在线永久视频| 国产高清国产精品国产三级| 50天的宝宝边吃奶边哭怎么回事| 老司机午夜福利在线观看视频| 久久久久国产一级毛片高清牌| 日本wwww免费看| 人人妻人人澡人人看| 精品久久久久久电影网| 欧美日韩一级在线毛片| 国产成人av教育| 精品国产一区二区久久| 久99久视频精品免费| 久久九九热精品免费| 成年人黄色毛片网站| 满18在线观看网站| 精品熟女少妇八av免费久了| 人人妻人人添人人爽欧美一区卜| 丰满饥渴人妻一区二区三| 另类亚洲欧美激情| 亚洲熟妇熟女久久| a在线观看视频网站| 精品久久久久久成人av| 色老头精品视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 1024香蕉在线观看| 一级片免费观看大全| 国产精品久久电影中文字幕| 色老头精品视频在线观看| 色在线成人网| 99热只有精品国产| 久久久精品欧美日韩精品| 国产成人av教育| 亚洲国产中文字幕在线视频| 十八禁网站免费在线| 老司机亚洲免费影院| 精品久久久精品久久久| 亚洲一区中文字幕在线| 欧美一区二区精品小视频在线| 丝袜在线中文字幕| 成人av一区二区三区在线看| 亚洲国产看品久久| 亚洲国产毛片av蜜桃av| 亚洲精品中文字幕在线视频| 99久久久亚洲精品蜜臀av| 日韩国内少妇激情av| 在线观看舔阴道视频| 色播在线永久视频| 国产伦一二天堂av在线观看| 又黄又爽又免费观看的视频| 水蜜桃什么品种好| 无限看片的www在线观看| 老熟妇乱子伦视频在线观看| 97人妻天天添夜夜摸| av免费在线观看网站| 国产精品98久久久久久宅男小说| 国产精品1区2区在线观看.| 精品一区二区三区视频在线观看免费 | 亚洲av第一区精品v没综合| 69精品国产乱码久久久| 操美女的视频在线观看| 色播在线永久视频| 女同久久另类99精品国产91| 一进一出好大好爽视频| 桃红色精品国产亚洲av| xxx96com| 亚洲国产欧美日韩在线播放| 国产精品自产拍在线观看55亚洲| 桃色一区二区三区在线观看| 国产精品99久久99久久久不卡| av在线播放免费不卡| 亚洲第一青青草原| 99香蕉大伊视频| 国产99久久九九免费精品| 三上悠亚av全集在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 国产成人精品久久二区二区91| 人成视频在线观看免费观看| 久久影院123| 亚洲人成网站在线播放欧美日韩| 亚洲第一青青草原| 99riav亚洲国产免费| 欧美日韩亚洲国产一区二区在线观看| 久热爱精品视频在线9| 久久天堂一区二区三区四区| 午夜福利欧美成人| 精品一区二区三卡| 1024香蕉在线观看| 欧美最黄视频在线播放免费 | 国产单亲对白刺激| 丰满迷人的少妇在线观看| 欧美人与性动交α欧美软件| 在线观看www视频免费| 久久国产乱子伦精品免费另类| 日本精品一区二区三区蜜桃| 免费女性裸体啪啪无遮挡网站| 成年人黄色毛片网站| 久99久视频精品免费| 成人黄色视频免费在线看| 中出人妻视频一区二区| 人妻丰满熟妇av一区二区三区| 亚洲专区字幕在线| 变态另类成人亚洲欧美熟女 | 99久久综合精品五月天人人| 亚洲精品国产色婷婷电影| 老熟妇乱子伦视频在线观看| 在线观看免费视频网站a站| 国产av又大| 中文字幕精品免费在线观看视频| 亚洲国产精品999在线| 亚洲av熟女| 亚洲精品av麻豆狂野| 少妇被粗大的猛进出69影院| 91麻豆av在线| 97碰自拍视频| 波多野结衣一区麻豆| 国产精品99久久99久久久不卡| 亚洲狠狠婷婷综合久久图片| 日日摸夜夜添夜夜添小说| 一进一出好大好爽视频| √禁漫天堂资源中文www| 黑人巨大精品欧美一区二区蜜桃| 中文字幕人妻熟女乱码| 天堂俺去俺来也www色官网| 国产成人精品久久二区二区91| 成人精品一区二区免费| 国产一区二区激情短视频| av天堂久久9| 国产xxxxx性猛交| x7x7x7水蜜桃| 久久久久久人人人人人| 久久精品亚洲av国产电影网| 欧美激情久久久久久爽电影 | 成人手机av| 50天的宝宝边吃奶边哭怎么回事| 80岁老熟妇乱子伦牲交| 校园春色视频在线观看| 999久久久国产精品视频| 日韩欧美免费精品| 男男h啪啪无遮挡| 亚洲激情在线av| 日本撒尿小便嘘嘘汇集6| 麻豆成人av在线观看| 露出奶头的视频| 成人18禁高潮啪啪吃奶动态图| 国产男靠女视频免费网站| 国产精品成人在线| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲人成电影观看| 亚洲午夜理论影院| 深夜精品福利| 欧美成人免费av一区二区三区| 日韩精品青青久久久久久| 日韩高清综合在线| 欧美精品啪啪一区二区三区| 久久精品人人爽人人爽视色| 我的亚洲天堂| 国产真人三级小视频在线观看| 日本撒尿小便嘘嘘汇集6| 午夜精品国产一区二区电影| 69精品国产乱码久久久| 国产精品免费视频内射| 亚洲一码二码三码区别大吗| 午夜福利影视在线免费观看| 琪琪午夜伦伦电影理论片6080| 国产亚洲欧美精品永久| 身体一侧抽搐| 一个人观看的视频www高清免费观看 | 亚洲黑人精品在线| 超碰成人久久| 久久久精品欧美日韩精品| 国产成人av激情在线播放| 午夜a级毛片| 91国产中文字幕| 精品少妇一区二区三区视频日本电影| www.www免费av| 老熟妇乱子伦视频在线观看| 精品免费久久久久久久清纯| 757午夜福利合集在线观看| 久久久精品欧美日韩精品| 成人精品一区二区免费| 精品久久久久久久毛片微露脸| 日韩大尺度精品在线看网址 | 欧美色视频一区免费| 欧美精品亚洲一区二区| 老鸭窝网址在线观看| 欧美性长视频在线观看| 男女之事视频高清在线观看| 岛国在线观看网站| 男女下面插进去视频免费观看| 日韩国内少妇激情av| 淫秽高清视频在线观看| 亚洲在线自拍视频| 日韩三级视频一区二区三区| 亚洲欧洲精品一区二区精品久久久| 两个人免费观看高清视频| 午夜久久久在线观看| 天天影视国产精品| 国产精品久久电影中文字幕| 亚洲 国产 在线| 如日韩欧美国产精品一区二区三区| 国产精品1区2区在线观看.| 一夜夜www| 99国产精品一区二区三区| 天天添夜夜摸| 亚洲专区字幕在线| 国产亚洲精品久久久久5区| 一本大道久久a久久精品| 99精国产麻豆久久婷婷| 午夜福利影视在线免费观看| 午夜亚洲福利在线播放| 欧美性长视频在线观看| 在线看a的网站| 成人永久免费在线观看视频| 国产伦人伦偷精品视频| 免费在线观看视频国产中文字幕亚洲| 国产成人精品无人区| a级毛片黄视频| 国产精品久久久久久人妻精品电影| www日本在线高清视频| 高潮久久久久久久久久久不卡| 一a级毛片在线观看| 精品欧美一区二区三区在线| 欧美黑人精品巨大| 久久久久精品国产欧美久久久| 久久精品aⅴ一区二区三区四区| 成人亚洲精品av一区二区 | 成年女人毛片免费观看观看9| 久久人人精品亚洲av| 两人在一起打扑克的视频| 女人高潮潮喷娇喘18禁视频| 亚洲成国产人片在线观看| av中文乱码字幕在线| 国产97色在线日韩免费| 欧美av亚洲av综合av国产av| www.精华液| 欧美黑人精品巨大| 久久国产精品人妻蜜桃| 日韩三级视频一区二区三区| 人人妻人人添人人爽欧美一区卜| 琪琪午夜伦伦电影理论片6080| 丝袜美足系列| 成年人免费黄色播放视频| 看免费av毛片| 国产一区二区三区在线臀色熟女 | 国产黄色免费在线视频| 欧美日韩黄片免| 国产av又大| 纯流量卡能插随身wifi吗| 国产成人av教育| 国产男靠女视频免费网站| 久久精品国产99精品国产亚洲性色 | 国产成人欧美| x7x7x7水蜜桃| 亚洲国产毛片av蜜桃av| 国产激情久久老熟女| 日韩欧美一区视频在线观看| 黑人猛操日本美女一级片| 黄色女人牲交| 制服人妻中文乱码| 国产一区二区三区综合在线观看| 精品高清国产在线一区| 一夜夜www| 欧美不卡视频在线免费观看 | 妹子高潮喷水视频| 日韩有码中文字幕| 女人高潮潮喷娇喘18禁视频| 大码成人一级视频| 精品人妻在线不人妻| 国产高清videossex| 亚洲第一青青草原| 啦啦啦在线免费观看视频4| 亚洲人成77777在线视频| 国产亚洲精品综合一区在线观看 | 亚洲中文字幕日韩| 如日韩欧美国产精品一区二区三区| 三上悠亚av全集在线观看| 精品国产一区二区久久| 免费少妇av软件| av网站在线播放免费| 精品国产美女av久久久久小说| 亚洲精品在线观看二区| 中文字幕人妻熟女乱码| av天堂在线播放| 欧美一级毛片孕妇| 国产蜜桃级精品一区二区三区| 黄片播放在线免费| 久久久国产精品麻豆| 老司机深夜福利视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 别揉我奶头~嗯~啊~动态视频| 国产精品影院久久| 亚洲avbb在线观看| 免费av毛片视频| 淫妇啪啪啪对白视频| 女生性感内裤真人,穿戴方法视频| 欧美乱码精品一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 国产97色在线日韩免费| 亚洲少妇的诱惑av| 国产精品秋霞免费鲁丝片| 在线观看www视频免费| 极品教师在线免费播放| 久久久国产欧美日韩av| 99香蕉大伊视频| 麻豆成人av在线观看| 日本a在线网址| 在线十欧美十亚洲十日本专区| 国内毛片毛片毛片毛片毛片| 久久 成人 亚洲| www.999成人在线观看| xxxhd国产人妻xxx| www.999成人在线观看| 精品久久蜜臀av无| av欧美777| 国产高清国产精品国产三级| 黄片小视频在线播放| 老司机福利观看| 午夜影院日韩av| 久久久国产成人精品二区 | 久久亚洲真实| 757午夜福利合集在线观看| 国产av精品麻豆| 淫妇啪啪啪对白视频| 午夜福利一区二区在线看| 校园春色视频在线观看| 精品久久久久久,| 亚洲免费av在线视频| 日本免费一区二区三区高清不卡 | 亚洲免费av在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 老司机午夜福利在线观看视频| 波多野结衣av一区二区av| 日韩大码丰满熟妇| 国产精品一区二区精品视频观看| 啦啦啦免费观看视频1| 操出白浆在线播放| 丝袜在线中文字幕| 国产熟女午夜一区二区三区| 色哟哟哟哟哟哟| 久久性视频一级片| 亚洲三区欧美一区| 超色免费av| 亚洲 欧美 日韩 在线 免费| 国产黄色免费在线视频| 老司机在亚洲福利影院| 69av精品久久久久久| 久久精品亚洲熟妇少妇任你| 日韩欧美免费精品| 老司机深夜福利视频在线观看| 久久精品国产亚洲av高清一级| 成人永久免费在线观看视频| 黄色视频不卡| 丰满迷人的少妇在线观看| 亚洲男人天堂网一区| 高清av免费在线| 丰满的人妻完整版| 最近最新中文字幕大全电影3 | 日韩欧美一区二区三区在线观看| 老司机午夜福利在线观看视频| 亚洲美女黄片视频| 精品久久蜜臀av无| 日韩精品青青久久久久久| 国产成人一区二区三区免费视频网站| 每晚都被弄得嗷嗷叫到高潮| 麻豆国产av国片精品| 欧美日韩国产mv在线观看视频| 精品乱码久久久久久99久播| 黑人欧美特级aaaaaa片| 国产成年人精品一区二区 | 深夜精品福利| 高潮久久久久久久久久久不卡| 亚洲精品美女久久久久99蜜臀| 精品久久久久久久久久免费视频 | 99热国产这里只有精品6| 在线观看免费视频日本深夜| 日韩大码丰满熟妇| 亚洲精品成人av观看孕妇| 男女之事视频高清在线观看| 女人爽到高潮嗷嗷叫在线视频| 最近最新中文字幕大全免费视频| 亚洲va日本ⅴa欧美va伊人久久| 久久精品国产亚洲av高清一级| 国产亚洲精品一区二区www| 18禁美女被吸乳视频| 在线观看66精品国产| 久久精品亚洲av国产电影网| 黄色 视频免费看| 99香蕉大伊视频| 自线自在国产av| 国产欧美日韩综合在线一区二区| 亚洲一区二区三区色噜噜 | 国产精品一区二区免费欧美| 亚洲成人国产一区在线观看| 人人妻人人添人人爽欧美一区卜| 久久香蕉国产精品| 91九色精品人成在线观看| 久久中文字幕一级| 欧美激情久久久久久爽电影 | 99在线视频只有这里精品首页| 悠悠久久av| 国产国语露脸激情在线看| 国产精品影院久久| 中文字幕色久视频| 国产不卡一卡二| av天堂久久9| 久久久水蜜桃国产精品网| 午夜免费鲁丝| 国产亚洲欧美在线一区二区| 亚洲av片天天在线观看| 久久精品国产亚洲av高清一级| 精品福利观看| 亚洲色图 男人天堂 中文字幕| 国产xxxxx性猛交| 久久 成人 亚洲| 国产伦人伦偷精品视频| 少妇 在线观看| 亚洲色图 男人天堂 中文字幕| 欧美av亚洲av综合av国产av| 日日夜夜操网爽| 一级,二级,三级黄色视频| 99久久人妻综合| 亚洲情色 制服丝袜| 最近最新中文字幕大全免费视频| 国产精品综合久久久久久久免费 | 亚洲欧美日韩高清在线视频| 99国产精品一区二区蜜桃av| 99久久人妻综合| 久久人人爽av亚洲精品天堂| 18禁美女被吸乳视频| 精品电影一区二区在线| 国产精品一区二区在线不卡| 久久久精品欧美日韩精品| 国产精品1区2区在线观看.| 桃色一区二区三区在线观看| 老汉色av国产亚洲站长工具| 亚洲自偷自拍图片 自拍| 欧美日韩一级在线毛片| 欧美午夜高清在线| 亚洲va日本ⅴa欧美va伊人久久| 女警被强在线播放| 国产一区在线观看成人免费| 黄片播放在线免费| 一个人免费在线观看的高清视频| 超碰成人久久| 亚洲熟妇中文字幕五十中出 | 亚洲精品一卡2卡三卡4卡5卡| 色精品久久人妻99蜜桃| 999久久久精品免费观看国产| x7x7x7水蜜桃| 99re在线观看精品视频| 最近最新中文字幕大全电影3 | 亚洲男人的天堂狠狠| 国产成年人精品一区二区 | 宅男免费午夜| 国产1区2区3区精品|