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

    基于SWCC試驗數(shù)據(jù)的壩體非飽和非穩(wěn)態(tài)滲流與穩(wěn)定性研究*

    2016-07-06 01:12:01陳佳林王英男王育恒
    工程地質學報 2016年2期
    關鍵詞:土水非飽和吸力

    陳佳林 金 鋒 陳 曦 王英男② 王育恒 任 俊

    基于SWCC試驗數(shù)據(jù)的壩體非飽和非穩(wěn)態(tài)滲流與穩(wěn)定性研究*

    陳佳林①金鋒①陳曦①王英男①②王育恒①任?、?/p>

    ( ①北京交通大學土木建筑工程學院北京100044)

    ( ②中國土木工程集團有限公司北京100038)

    摘要飽和-非飽和狀態(tài)是土存在于自然界的真實狀態(tài),描述和解釋這種狀態(tài)的非飽和土力學理論在土石壩滲流、污染物傳輸、凍土滲流相變、邊坡和路基穩(wěn)定性分析等領域有著廣泛的應用。非飽和土的土水特征曲線(SWCC)是非飽和土力學研究的基本內容之一,對非飽和土體滲流和穩(wěn)定性分析至關重要。采用Python語言開發(fā)了非飽和滲流與穩(wěn)定性分析軟件包USSA和非飽和土土水特征曲線試驗數(shù)據(jù)處理和模型擬合界面。對非飽和土的土水特征曲線數(shù)據(jù)進行處理,以模型擬合參數(shù)作為基本輸入進行非飽和滲流場的模擬,再到非飽和土的穩(wěn)定性分析,詳細呈現(xiàn)了非飽和土滲流與穩(wěn)定性分析及軟件開發(fā)的全過程?;赟WCC模型的擬合參數(shù)對壩體滲流穩(wěn)定性進行了分析,分析結果表明水位變化過程中壩體兩側斜坡具有明顯不同的穩(wěn)定性演化規(guī)律。當采用有限元強度折減法進行非飽和土斜坡穩(wěn)定性分析時,最終安全系數(shù)為所有邊坡安全系數(shù)演化曲線的最低包絡線。

    關鍵詞非飽和土滲流和穩(wěn)定性分析土水特征曲線(SWCC)Python腳本語言有限元強度折減法

    陳曦(1977-),男,博士,教授,博士生導師,主要從事計算巖土力學、特殊巖土力學、巖土工程風險方面的教學和科研.Email:xichen.geo@gmail.com

    0引言

    從土力學建立以來,飽和土力學一直主導著土力學的工程應用。其主要原因可以歸結為以下兩方面:(1)從強度角度來講,非飽和土的強度理論通??紤]基質吸力的影響,而飽和土力學強度理論不包括基質吸力的強度貢獻。因此,飽和土力學強度理論通??梢詫е赂颖J?或安全)的工程設計; (2)非飽和土力學理論涉及更加復雜的土-水-氣相互作用,其作用機理還需要深入探究,很多理論還有待完善(Shengetal.,2013)。盡管如此,由于非飽和狀態(tài)具有普遍意義,是土的真實存在狀態(tài),采用非飽和土力學理論和方法進行巖土工程問題研究的意義不僅在于獲得更加精確物理描述,還在于獲得針對巖土體滲流、變形和破壞機理的更加理性的解釋。因此,出于實際需要,土力學研究正在從飽和土力學進入到非飽和土力學。

    土水特征曲線是非飽和土力學研究區(qū)別于飽和土力學研究的主要特征之一。對于非飽和土的土水特征曲線模型,目前主要有Brooks-Corey模型、vanGenuchten模型和Fredlund-Xing模型,Luetal.(2004)的專著對上面3種模型進行了詳細的介紹,并對它們的適用性進行了闡述。毛尚之(2002)對某膨脹土的土水特征曲線進行了測定,并對滲析技術和軸平移技術的結果進行了對比,初步分析了應力及吸力歷史等因素對土水特征曲線的影響。盧應發(fā)等(2008)建立了一種SWCC曲線方程,并以實驗獲得了該方程的各項參數(shù),試驗結果表明SWCC曲線受到物質成分、塑性指數(shù)、孔隙結構、應力狀態(tài)等多種因素影響,并建立了SWCC曲線與非飽和滲流曲線的相互關系。譚曉慧等(2013)認為在對SWCC曲線進行擬合時,試驗數(shù)據(jù)點需要具有較大的吸力范圍,才能很好地給出整條曲線的形狀,否則會造成多解問題。試驗點不必過于密集,但需要具有整條曲線的代表性。劉艷等(2008)研究了砂土SWCC在干濕往復條件下形成的干燥曲線與浸濕曲線滯回圈和相應的滯回行為,提出了基于熱力學理論的具有一般性和普遍適用性的土水特征曲線滯后模型。通過數(shù)值模擬與實驗結果的對比,驗證了所提出模型的有效性。

    通過非飽和土SWCC曲線可以建立基質吸力與土體體積含水率或含水量的關系,通過基質吸力和非飽和土有效應力原理可以建立非飽和土的強度理論,根據(jù)非飽和土的強度可以分析和研究飽和-非飽和土體的穩(wěn)定性。目前大部分非飽和土穩(wěn)定性研究都遵從以上步驟(Ngetal., 2000; 黃潤秋等, 2002;Griffithsetal., 2005),例如,于玉貞等(2008)和林鴻州等(2009)基于所獲得的非飽和滲流場定義了非飽和土的抗剪強度,采用有限元強度折減法對坡體的穩(wěn)定性進行了分析。上述分析方法的主要差別在于所采用的滲流場模型(包括SWCC模型和滲透性模型)不同或是邊坡穩(wěn)定性分析方法不同。例如,張磊等(2014)對于土水特征曲線方程采用了vanGenuchten方程,非飽和滲透系數(shù)模型采用的是Gardner模型,并考慮了滲流潛蝕對孔壓分布和邊坡穩(wěn)定性的影響; 邊坡穩(wěn)定性分析方法則可采用極限平衡法、極限分析法、重力加載法和強度折減法等。例如,Slope/W就采用SWCC試驗數(shù)據(jù)輸入和極限平衡法進行降雨入滲分析。因此土水特征曲線的擬合影響了非飽和土滲流場的質量,進而又影響了非飽和土邊坡的穩(wěn)定性分析。各種因素可能導致SWCC的不確定性,Phoonetal.(2010)對SWCC曲線的不確定性進行了研究,建議采用對數(shù)正態(tài)隨機向量來構建SWCC曲線擬合參數(shù)的概率模型??梢姡琒WCC曲線的不確定性會導致滲流場和邊坡穩(wěn)定性分析的不確定,是未來非飽和土坡穩(wěn)定性分析可能面臨的課題。此外,SWCC應依賴于當前孔隙比,因此隨著孔隙比變化的SWCC會形成一空間曲面。

    盡管SWCC曲線與孔隙比有關,本文基于時間相關的滲流-變形不完全耦合分析方法(陳曦等, 2011),在不考慮孔隙比的情況下,從非飽和土SWCC曲線試驗數(shù)據(jù)采集,到數(shù)據(jù)處理和模型參數(shù)擬合,再到非飽和土體穩(wěn)定性分析,詳細呈現(xiàn)了非飽和土滲流與穩(wěn)定性分析及軟件開發(fā)的全過程。通過所開發(fā)的軟件對水位變化過程中均質壩體和心墻壩體滲流場和穩(wěn)定性進行了分析。

    1非飽和土土水特征曲線試驗、模型與數(shù)據(jù)擬合

    1.1非飽和土土水特征曲線試驗

    測定非飽和土含水量與吸力之間的關系(即土水特征曲線)可采用多種方法,李志清等(2007)介紹了6種常用方法,并評述了這些方法的優(yōu)缺點。圖1 所示為VJTech公司的非飽和土固結儀,該固結儀采用軸平移技術進行土水特征曲線的測量。軸平移技術是測量土體基質吸力us(即us=ua-uw)的主要方法之一,其原理是通過逐漸提升孔隙氣壓力ua,而用飽和高進氣值材料將孔隙水壓力uw維持在某個參考值,實現(xiàn)對基質吸力進行直接控制,并可避免對于孔隙氣壓力為零,基質吸力超過101kPa時陶土板中水的汽化現(xiàn)象。

    圖1 VJ Tech非飽和土固結儀的整體布置Fig. 1 VJ Tech Ltd unsaturated triaxial testing system

    試驗采用VJTech公司所生產(chǎn)的儀器,采集的數(shù)據(jù)類型為第i級基質吸力us條件下,土樣排水量Vw與時間t的關系為:

    (1)

    當前基質吸力條件下,試驗穩(wěn)定后的累計排水量Vw, i+1與試驗開始時刻的累計排水量Vw, i的差值為us基質吸力條件下土樣的排水量ΔVw,根據(jù)試驗初始時的含水量w0與土樣體積V即可計算出us基質吸力條件下土樣的體積含水率w:

    (2)

    本次試驗所施加的基質吸力等級依次為1kPa, 3kPa, 5kPa, 10kPa, 20kPa, 50kPa, 150kPa, 250kPa,所用陶土板的進氣值為5bar(1bar=100kPa), 24h內的排水量不超過0.1cm3時認為達到吸力平衡。圖2 所示為各級吸力下砂土的累計排水量隨時間的變化曲線。對于砂土,在初始較小的基質吸力(如1kPa)作用下,就會產(chǎn)生較大的排水量,而此時的水流驅動力較小,因此排水時間比較長(約17d); 之后各級吸力下的排水量明顯減少,達到吸力平衡的時間也較短(約3~5d)。

    圖2 典型砂土累計排水量隨時間的變化關系Fig. 2 Accumulated quantity of water outflow varying with time for typical sand

    1.2常用的土水特征曲線模型

    常用的土水特征曲線模型有Brooks-Corey(1964)模型,vanGenuchten(1980)模型和Fredlund-Xing(1994)模型。這些模型中部分參數(shù)(如飽和含水量θs和殘余含水量θr等)具有實際物理意義,也有部分參數(shù)沒有物理意義。下面僅詳細給出vanGenuchten模型的具體表達式。

    圖3 SWCC試驗數(shù)據(jù)擬合結果Fig. 3 Fitting results of the SWCC experiment

    vanGenuchten于1980年提出了一個土水特征曲線的模型(簡稱VG模型)。相比BC模型,VG模型具有連續(xù)平滑曲線,在進氣壓力值和趨近殘余含水量狀態(tài)時具有較好的平滑過渡,因此得到了廣泛的應用。VG模型表達式如下:

    (3)

    或表示為標準化的體積含水量Θ

    (4)

    式中,av、nv和mv為經(jīng)驗系數(shù),mv通常簡化為mv=1-1/nv。VG模型是上述3參數(shù)的模型,對于同一干密度條件下的土,θs和θr可由土的基本物理力學試驗直接測出,當沒有θs和θr的實測數(shù)據(jù)時,可將該模型按5參數(shù)模型進行擬合。

    1.3土水特征曲線數(shù)據(jù)擬合與界面開發(fā)

    進行土水特征曲線試驗時可得到一系列的Excel數(shù)據(jù)文件,每個Excel文件對應某一級基質吸力作用下的排水過程,或者說每個Excel文件表示土水特征曲線上的一個數(shù)據(jù)點(us,θ)。點擊“選擇文件夾”選擇存儲SWCC數(shù)據(jù)文件的文件夾; 接下來點擊“獲取所有Excel文件”將該文件夾中的Excel文件導入到列表窗口,并勾選列表窗口的Excel數(shù)據(jù)文件,每個Excel文件對應土水特征曲線上的一個數(shù)據(jù)點。同時需要填寫本試驗的初始水的體積和試樣體積。對數(shù)據(jù)點整理和認定后生成一個新的Excel文件,對應著不同基質吸力條件下的排水量,組合成一個新的數(shù)據(jù)表即可進行擬合。圖3 展示了SWCC可視化界面,具體分為操作區(qū)、圖像顯示區(qū)域和數(shù)據(jù)輸出區(qū)域??梢暬到y(tǒng)所采用的擬合方法為非線性最小二乘法,圖像顯示區(qū)域給出了根據(jù)試驗數(shù)據(jù)擬合的SWCC曲線,同時將擬合的模型參數(shù)傳遞給數(shù)值模擬界面。

    2非飽和滲流的有限元求解與非飽和土穩(wěn)定性分析

    采用Richards方程進行非飽和土滲流模擬已經(jīng)得到了廣泛的認可。Richards方程有3種基本格式,即壓力水頭格式(h-form)、混合格式(mixedform)和體積含水量格式(θ-form)。其中壓力水頭格式在實際工程中的應用更為廣泛。

    非飽和滲流的水力傳導系數(shù)可采用Mualem式(1976),即:

    (5)

    式中: 飽和水力傳導系數(shù)ks為常數(shù);kr(Θ)為相對滲透張量, 為有效飽和度Θ的函數(shù); 而非飽和滲流的水力傳導系數(shù)則隨體積含水量的減小而迅速衰減。

    壓力水頭格式(h-form)的Richards方程如下:

    (6)

    式中,(x,t)∈Ω×(0,T],s為源匯項;η為比存儲量(即土水特征曲線的導數(shù))。

    (7)

    式中: γw為水的單位重度; mw為土水特征曲線的斜率??臻g離散和時間差分后可得h-form的Richards有限元方程的修正Picard迭代公式:

    (8)

    式中,

    (9)

    其中,下標n表示時間步指標; m表示修正Picard的迭代步指標; M,C分別為質量矩陣和水力傳導特征矩陣。

    非飽和土坡穩(wěn)定性分析方法主要分為3大類,即時間無關的定性分析方法,時間相關的滲流與變形不完全耦合分析方法和時間相關的滲流與變形完全耦合分析方法。盡管完全耦合分析方法已經(jīng)成為發(fā)展趨勢,從發(fā)展成熟程度和計算效率角度考慮,時間相關的滲流與變形不完全耦合分析方法仍為當前的主流分析方法,其基本思路是先進行非飽和滲流分析,基于某一時刻的瞬態(tài)滲流場,再采用極限平衡法或強度折減有限元法進行飽和-非飽和土坡的穩(wěn)定性分析。

    非飽和土邊坡穩(wěn)定性評價需要采用非飽和土抗剪強度理論??娏植?1999)曾對非飽和土的抗剪強度進行了研究,提出了基于試驗數(shù)據(jù)擬合的雙曲模型。而常用的非飽和土的抗剪強度理論仍以Bishop和Blight強度理論以及Fedlund強度理論為主。根據(jù)非飽和土的有效應力原理:

    (10)

    式中:χ為與飽和度(Sr)相關的參數(shù)。則Bishop和Blight強度理論可表示為:

    (11)

    Fredlund強度理論則表示為:

    (12)

    式中,c′, φ′分別為土的有效黏聚力和有效內摩擦角; φb是一個描述基質吸力對剪切強度貢獻的角度,可見,F(xiàn)redlund強度理論區(qū)別對待了φb對基質吸力影響,而目前各種非飽和土強度理論的主要差別在于χ和tanφb的選取的不同。為了計算方便,非飽和土的抗剪強度也常記為:

    (13)

    式中,“假黏聚力”c″由c′和基質吸力影響(即ustanφb)兩部分組成。式(13)中,根據(jù)Shengetal.(2011),tanφb具有不同的定義,例如:

    (14)

    式中,當參數(shù)κ取為1.0,F(xiàn)redlund-Fredlund公式退化為Fredlund-Oberg公式。應用強度折減法時,需要進行如下折減:

    (15)

    式中,SRF為強度折減因子。需要強調的是,當在c″考慮基質吸力的影響時,由于積分點上的基質吸力不同,c″對于每個單元并不是常數(shù),需要進行均勻化處理,取積分點的平均值作為該單元的黏聚力指標。

    與傳統(tǒng)的極限平衡法和極限分析法相比,有限元強度折減法具有一些優(yōu)勢,除了能夠給出邊坡的安全系數(shù),還可以得到位移、應力、應變場,塑性屈服點的分布以及邊坡失穩(wěn)破壞的全過程。由于有限元強度折減法需要處理更多的信息,相比極限平衡法和極限分析法,自然需要更多的計算時間,因此如何提高有限元強度折減法的計算效率成為需要關注的一個問題(Chenetal.,2014)。

    3非飽和土滲流與穩(wěn)定性分析程序USSA的計算流程

    圖4 USSA非飽和土滲流與穩(wěn)定性分析軟件包的分析流程Fig. 4 Unsaturated seepage and stability analysis flowchart for USSA software package

    基于上述理論方法,開發(fā)了非飽和土滲流與穩(wěn)定性分析軟件包。圖4 所示為非飽和土滲流與穩(wěn)定性分析軟件包USSA的分析和計算流程。首先給定一系列初始的Excel文件,生成繪制SWCC曲線所需要的兩列數(shù)據(jù)(us,θ),如1.3節(jié)所述。根據(jù)非飽和土的類型選擇合適的SWCC模型進行試驗數(shù)據(jù)擬合,或直接采用原始數(shù)據(jù)(us,θ)進行擬合或插值使用。這里默認選取Mualem公式為非飽和水力傳導系數(shù)公式,當然可進一步擴展并加入其他非飽和水力傳導系數(shù)公式。根據(jù)選用的SWCC模型和擬合得到的模型參數(shù)進行非飽和滲流場模擬,當需要進行土體穩(wěn)定性分析時,可以在進行滲流場計算的同時,進行穩(wěn)定性分析。最后可視化顯示非飽和滲流場,含水量分布,位移場和塑性點分布等。USSA軟件的用戶界面和SWCC界面都采用Python語言開發(fā),計算程序則采用FORTRAN語言開發(fā)。目前,Python作為一種開源的腳本編程語言,已經(jīng)得到了廣泛的應用。

    4數(shù)值應用

    4.1均質土壩算例

    圖5為均質土壩,即壩身土體和壩基土體土性參數(shù)相同。土體有效彈性參數(shù)取E′=10MPa,ν′=0.3; 摩爾-庫侖有效塑性參數(shù)取為c′=10kPa,φ′=20°(此處采用相關塑性流法則); 自然和飽和重度分別為17kN·m-3和20kN·m-3,ks均為1e-6m·s-1; 非飽和滲流土水特征曲線采用VG模型,模型參數(shù)θs,θr,av(1·kPa-1)和nv(mv=1-1/nv)分別為0.367, 0.187, 1.0和1.53。土壩左右邊坡坡度分別為1︰1.8和1︰1.4, 土壩左側初始水位為0m,右側初始水位為7.05m。在此初始條件下達穩(wěn)態(tài)滲流,相應的陰影圖(圖6)。右側水位從初始水位以1m·d-1的速度下降,最后水位為0m。

    圖5 土壩有限元網(wǎng)格劃分與幾何尺寸Fig. 5 Finite element mesh and geometry of soil dam

    圖6 初始高水位時的陰影圖Fig. 6 At the initial water level, contour plotting for(a) Pressure head; (b) Volumetric water content; (c) Shear strain; (d) Plastic zone; (e) Displacement modulus |u|a.壓力水頭; b.體積含水量; c.剪應變; d.塑性區(qū); e.位移場|u|

    從圖6a、圖6b可知,自由液面以下土體處于飽和狀態(tài),自由液面以上壩體的體積含水量隨著高程增加而逐漸減小。在初始水位和滲流條件下,土壩左側坡體的安全系數(shù)最小,為FOS=1.528,失穩(wěn)時的剪應變、塑性點分布和位移場|u|分別如圖6c、圖6d和圖6e所示。

    圖7 最低水位時的陰影圖Fig. 7 At the final water level, contour plotting for(a) Pressure head; (b) Volumetric water content; (c) Shear strain; (d) Plastic zone; (e) Displacement modulus |u|a.壓力水頭; b.體積含水量; c.剪應變; d.塑性區(qū); e.位移場|u|

    圖7所示為右側水位下降為0m時的陰影圖。從圖7a、7b可見,盡管土壩左右兩側水位都為0m,由于時間效應和土體的持水特性,壩體中部自由液面的高度要大于兩側的高度。此外,最危險的滑坡從左側轉移到右側,如圖7c、7d和7e所示,此時坡體安全系數(shù)為FOS=1.398。圖8 完整地呈現(xiàn)了水位下降過程中,土壩坡體安全系數(shù)的變化過程,即安全系數(shù)先穩(wěn)定增加,達到某一值(FOS=1.567)之后開始迅速減小,也正是從這一點開始,最危險的滑坡從左側轉移到右側。因此,可以推斷,當水位從初始高度開始下降時,左右側坡體中的自由液面開始下降,坡體中非飽和土體比例都在增加,而自由液面以下的飽和土體比例逐漸減小。由于飽和土體自重大、強度低,而非飽和土體自重較小,強度較高。自由液面下降導致左側坡體安全系數(shù)逐漸增大,而對于右側坡體則相對復雜,這是由于水位下降過程中,右側水體對坡面的正壓力也在下降,不利于坡體的安全。在水位較高時,右側水體對坡面的正壓力較大,導致右側坡體的安全系數(shù)大于左側坡體的安全系數(shù),但隨著水位的下降,右側坡體的安全系數(shù)迅速減小,并在水位下降超過4m時,開始小于左側坡體的安全系數(shù)。因此,圖8 所示的安全系數(shù)變化曲線可以看成是左側坡體和右側坡體安全系數(shù)的交叉組合變化曲線中安全系數(shù)最小的部分。為了獲得左右邊坡完整的安全系數(shù)變化曲線,可采用極限平衡法進行求解。如何采用有限元強度折減法獲得左右邊坡完整的安全系數(shù)變化曲線而不產(chǎn)生明顯的結果偏差還有待研究。

    圖8 水位降低過程中均質土壩安全系數(shù)的變化曲線Fig. 8 Variation curve of safety factor(FOS)with the fall of water level for homogeneous dam

    4.2非均質土壩算例

    對于非均質土壩,仍采用上述圖5 所示的均質土壩的有限元網(wǎng)格劃分,但壩體和壩基則分別采用非飽和土數(shù)據(jù)庫UNSODA中編號為1183的黏土和編號為3091的砂質壤土(Nemesetal.,1999)。采用圖3 所示的SWCC數(shù)據(jù)處理和模型擬合界面分別對這兩種非飽和土的試驗數(shù)據(jù)進行VG模型擬合,得到的水力參數(shù)如表1所示,擬合曲線如圖9 所示。壩體壩基彈性參數(shù)均為E′=10MPa,ν′=0.3; 采用摩爾-庫侖強度模型,壩體參數(shù)為c′=30kPa,φ′=25°,壩基為c′=20kPa,φ′=35°; 壩體壩基自然和飽和重度分別為17kN·m-3和20kN·m-3。

    表1 非均質土壩水力參數(shù)

    圖9 基于VG模型擬合的SWCC參數(shù)Fig.9 Fitted parameters of SWCC based on VG modela. 編號為1183的粘土; b. 編號為3091的砂質壤土

    非均質土壩左右側初始水位與上述均質土壩相同,但右側水位從初始水位以0.01m·d-1的速度下降,最后水位為0m。圖10 所示為水位下降過程中,土壩坡體安全系數(shù)的演化過程。與均質土壩中安全系數(shù)的變化相似,當水位下降4m左右時,潛在破壞坡體從左側轉移到右側。不過由于壩體和壩基的強度較高,坡體的安全系數(shù)也在較大值附近波動。

    圖10 水位降低過程中非均質土壩安全系數(shù)的變化曲線Fig. 10 Variation curve of safety factor(FOS)with the fall of water level for heterogeneous dam

    5結語

    非飽和狀態(tài)是土的真實存在狀態(tài)。非飽和土體滲流和穩(wěn)定性分析是非飽和土力學研究的重要內容。本文從非飽和土SWCC曲線試驗開始,到SWCC模型參數(shù)擬合,再到滲流和穩(wěn)定性分析,詳細呈現(xiàn)了整個分析過程,開發(fā)了SWCC數(shù)據(jù)處理和模型擬合界面以及滲流穩(wěn)定性分析程序。通過數(shù)值算例,得到了一些結論:

    (1)由于飽和土和非飽和土具有不同的抗剪強度和容重,非飽和滲流場的模擬精度對壩體穩(wěn)定性分析具有顯著的影響??辜魪姸鹊脑u價也是基于不同應力點的“假黏聚力”的定義,根據(jù)不同的原理定義的“假黏聚力”不同,因此也會對壩體穩(wěn)定性的分析結果產(chǎn)生影響。

    (2)水位下降過程中,壩體滲流場發(fā)生相應的變化,兩側壩體斜坡的穩(wěn)定性具有明顯不同的演化規(guī)律。一般來說,壩體下游斜坡的穩(wěn)定性影響因素單一,主要隨著上游水位降低非飽和區(qū)逐漸增加,導致下游斜坡穩(wěn)定性逐漸增加; 而上游斜坡穩(wěn)定性由增加的非飽和區(qū)和撤離的水壓共同支配。當采用有限元強度折減法進行非飽和土斜坡穩(wěn)定性分析時,安全系數(shù)為所有邊坡安全系數(shù)演化曲線的最低包絡線。

    參考文獻

    ChenX,LiuCJ. 2011.Stageddevelopmentoffiniteelementmethodsforstabilityofunsaturatedsoilslopes[J].ChineseJournalofGeotechnicalEngineering, 33(S1): 380~384.

    ChenX,WuYK,YuYZ,etal. 2014.Atwo-gridsearchschemeforlarge-scale3-Dfiniteelementanalysesofslopestability[J].ComputersandGeotechnics,62: 203~215.

    GriffithsDV,LuN. 2005.Unsaturatedslopestabilityanalysiswithsteadyinfiltrationorevaporationusingelasto-plasticfiniteelements[J].InternationalJournalforNumericalandAnalyticalMethodsinGeomechanics,29(3): 249~267.

    HuangRQ,QiGQ. 2002.Theeffectofunsaturatedsoilsuctiononslopestability[J].JournalofEngineeringGeology,10(4): 343~348.

    LiZQ,LiT,HuRL,etal. 2007.MethodsfortestingandpredictingofSWCCinunsaturatedsoilmechanics[J].JournalofEngineeringGeology,15(5): 700~707.

    LinHZ,YuYZ,LiGX,etal. 2009.Onapplicationofsoil-watercharacteristiccurvestolandslideforecast[J].ChineseJournalofRockMechanicsandEngineering, 28(12): 2569~2576.

    LiuY,ZhaoCG. 2008.Hysteresismodelforsoil-watercharacteristiccurves[J].ChineseJournalofGeotechnicalEngineering,30(3): 399~405.

    LuN,LikosWJ. 2004.Unsaturatedsoilmechanics[M].Wiley.

    LuYF,ChenGF,LuoXQ,etal. 2008.Studyofsoil-watercharacteristicalcurveanditsinfluentialfactors[J].RockandSoilMechanics, 29(9): 2481~2486.

    MaoSZ. 2002.Astudyonsoil-watercharacteristiccurveofunsaturatedsoil[J].JournalofEngineeringGeology,10(2): 129~133.

    MiaoLC,YinZZ. 1999.Shearstrengthofunsaturatedsoils[J].RockandSoilMechanics,20(3): 1~6.

    NemesA,SchaapMG,LeijFJ. 1999.TheUNSODAunsaturatedsoilhydraulicdatabase(Version2.0)[R].USSalinityLaboratory,Riverside,CA.

    NgCWW,PangYW. 2000.Influenceofstressstateonsoil-watercharacteristicsandslopestability[J].JournalofGeotechnicalandGeoenvironmentalEngineering,126(2): 157~166.

    PhoonKK,SantosoA,QuekST. 2010.Probabilisticanalysisofsoil-watercharacteristiccurves[J].JournalofGeotechnicalandGeoenvironmentalEngineering,136(3): 445~455.

    ShengDC,ZhouAN,F(xiàn)redlundDG. 2011.Shearstrengthcriteriaforunsaturatedsoils[J].GeotechnicalandGeologicalEngineering,29(2): 145~159.

    ShengDC,ZhangS,YuZW. 2013.Unansweredquestionsinunsaturatedsoilmechanics[J].ScienceChinaTechnologicalSciences,56(5): 1257~1272.

    TanXH,YuW,ShenMF,etal. 2013.Experimentalstudyandcurvefittingofsoil-watercharacteristiccurve[J].RockandSoilMechanics, 34(S2): 51~56.

    YuYZ,LinHC,LiRJ,etal. 2008.Stabilityanalysisofunsaturatedsoilslopeundertransientseepageflowstate[J].RockandSoilMechanics,29(11): 2892~2898.

    ZhangL,ZhangLL,ChengY,etal. 2014.Slopestabilityunderrainfallinfiltrationconsideringinternalerosion[J].ChineseJournalofGeotechnicalEngineering,36(9): 1680~1687.

    陳曦,劉春杰. 2011. 逐步完善的非飽和土邊坡穩(wěn)定性有限元分析方法[J]. 巖土工程學報, 33(S1): 380~384.

    黃潤秋,戚國慶. 2002. 非飽和滲流基質吸力對邊坡穩(wěn)定性的影響[J]. 工程地質學報,10(4): 343~348.

    李志清,李濤,胡瑞林,等. 2007. 非飽和土土水特征曲線(SWCC)測試與預測[J]. 工程地質學報,15(5): 700~707.

    林鴻州,于玉貞,李廣信,等. 2009. 土水特征曲線在滑坡預測中的應用性探討[J]. 巖石力學與工程學報, 28(12): 2569~2576.

    劉艷,趙成剛. 2008. 土水特征曲線滯后模型的研究[J]. 巖土工程學報,30(3): 399~405.

    盧應發(fā),陳高峰,羅先啟,等. 2008. 土-水特征曲線及其相關性研究[J]. 巖土力學, 29(9): 2481~2486.

    毛尚之. 2002. 非飽和膨脹土的土-水特征曲線研究[J]. 工程地質學報,10(2): 129~133.

    繆林昌,殷宗澤. 1999. 非飽和土的剪切強度[J]. 巖土力學,20(3): 1~6.

    譚曉慧,余偉,沈夢芬,等. 2013. 土-水特征曲線的試驗研究及曲線擬合[J]. 巖土力學, 34(S2): 51~56.

    于玉貞,林鴻州,李榮建,等. 2008. 非穩(wěn)定滲流條件下非飽和土邊坡穩(wěn)定分析[J]. 巖土力學,29(11): 2892~2898.

    張磊,張璐璐,程演,等. 2014. 考慮潛蝕影響的降雨入滲邊坡穩(wěn)定性分析[J]. 巖土工程學報,36(9): 1680~1687.

    UNSATURATEDTRANSIENTFLOWINDAMANDITSSTABILITYANALYSISUSINGSWCCEXPERIMENTALDATA

    CHEN Jialin①JIN Feng①CHEN Xi①WANG Yingnan①②WANG Yuheng①REN Jun①

    (①SchoolofCivilEngineering,BeijingJiaotongUniversity,Beijing100044)

    (②ChinaCivilEngineeringConstructionCorporation,Beijing100038)

    AbstractSaturated-unsaturated state is a natural state of soils. Hence, unsaturated soil mechanics has wide applications in seepage analysis of earth and rockfill dam, contaminant transport, seepage and phase transform of frozen soil, and stability analysis of slope and embankment. As a fundamental relationship, the soil-water characteristic curve(SWCC)plays an essential role in the research of unsaturated soil mechanics. Python language is employed to develop the unsaturated seepage and stability analysis(USSA)software as well as the graphic user interface(GUI)for data processing and data fitting of SWCC curve. Based on a practical slope example, the whole analyzing process is described in details. It is from the data processing and data fitting of SWCC curve to unsaturated seepage analysis and stability analysis. The stability of dam under unsaturated seepage is conducted using the fitted parameters of SWCC model. The numerical results show that the stability evolutions of two dam slopes are quite different under the changing process of water level. When applying the shear strength reduction finite element method to the slope stability of unsaturated soil, the final FOS appears to be the lowest envelope of those FOS curves of all slopes.

    Key wordsUnsaturated soil, Seepage and stability analysis, Soil-water characteristic curve(SWCC), Python script language, Shear strength reduction finite element method

    DOI:10.13544/j.cnki.jeg.2016.02.012

    * 收稿日期:2015-02-04; 收到修改稿日期: 2015-04-04.

    基金項目:國家自然科學基金項目(51379103), 國家“973”項目課題(2012CB026104), 教育部博士點新教師基金(20110009120020)資助.

    第一作者簡介:陳佳林(1994-),男,學士,主要從事巖土工程及結構工程方面的科研和軟件開發(fā)工作. Email:12231037@bjtu.edu.cn

    中圖分類號:TV139.14

    文獻標識碼:A

    猜你喜歡
    土水非飽和吸力
    深水大型吸力錨測試技術
    初始孔隙比對高吸力下非飽和土土水特性的影響
    ROV在海上吸力樁安裝場景的應用及安裝精度和風險控制
    化工管理(2022年11期)2022-06-03 07:08:24
    基于粒徑分布曲線的非飽和砂土土水特征曲線概率預測模型*
    深水吸力樁施工技術研究
    非飽和原狀黃土結構強度的試驗研究
    工程與建設(2019年1期)2019-09-03 01:12:24
    干密度對重塑花崗巖殘積土土水特征曲線影響
    山西建筑(2019年3期)2019-01-19 07:15:54
    非飽和多孔介質應力滲流耦合分析研究
    非飽和土基坑剛性擋墻抗傾覆設計與參數(shù)分析
    非飽和地基土蠕變特性試驗研究
    麻豆av在线久日| 在线天堂中文资源库| 午夜福利免费观看在线| 国产深夜福利视频在线观看| 国产亚洲精品第一综合不卡| 在线观看人妻少妇| 巨乳人妻的诱惑在线观看| 亚洲人成网站在线观看播放| 女性生殖器流出的白浆| 国产成人精品久久久久久| 欧美激情高清一区二区三区| 一级片免费观看大全| 日韩熟女老妇一区二区性免费视频| 亚洲人成77777在线视频| 欧美日韩综合久久久久久| 少妇裸体淫交视频免费看高清 | 久久av网站| 两个人免费观看高清视频| 青春草视频在线免费观看| 久久av网站| av网站免费在线观看视频| 亚洲激情五月婷婷啪啪| 亚洲av成人精品一二三区| 色综合欧美亚洲国产小说| √禁漫天堂资源中文www| 黄色怎么调成土黄色| 9191精品国产免费久久| 午夜激情av网站| 超碰97精品在线观看| av天堂在线播放| 成人国产一区最新在线观看 | 日本色播在线视频| 免费不卡黄色视频| 欧美国产精品一级二级三级| 黄色毛片三级朝国网站| 777久久人妻少妇嫩草av网站| 久久久久久久大尺度免费视频| 天天影视国产精品| 成人亚洲精品一区在线观看| 高清黄色对白视频在线免费看| 精品福利永久在线观看| 最新的欧美精品一区二区| 日韩大码丰满熟妇| 高清黄色对白视频在线免费看| 国产高清不卡午夜福利| av有码第一页| 亚洲一区二区三区欧美精品| 日本av免费视频播放| 欧美亚洲 丝袜 人妻 在线| 两人在一起打扑克的视频| 精品第一国产精品| 啦啦啦视频在线资源免费观看| 超碰成人久久| 久久久久国产精品人妻一区二区| 久久久久久久国产电影| 日韩精品免费视频一区二区三区| 在线精品无人区一区二区三| 免费女性裸体啪啪无遮挡网站| 后天国语完整版免费观看| 欧美日韩视频精品一区| 99香蕉大伊视频| 欧美xxⅹ黑人| 国产亚洲精品久久久久5区| 久久99一区二区三区| 国产成人精品久久久久久| 日韩av不卡免费在线播放| 1024视频免费在线观看| 两性夫妻黄色片| 如日韩欧美国产精品一区二区三区| 久久久亚洲精品成人影院| 日韩av不卡免费在线播放| 亚洲国产av新网站| 各种免费的搞黄视频| 亚洲综合色网址| 丰满人妻熟妇乱又伦精品不卡| av不卡在线播放| 丝袜喷水一区| 人成视频在线观看免费观看| 免费av中文字幕在线| 国产片特级美女逼逼视频| 日韩视频在线欧美| 一级,二级,三级黄色视频| 免费久久久久久久精品成人欧美视频| 久久国产精品影院| 中文字幕人妻丝袜一区二区| 精品熟女少妇八av免费久了| 91麻豆av在线| 丝袜美足系列| 日本五十路高清| 国产精品一区二区在线观看99| 女警被强在线播放| 午夜日韩欧美国产| 欧美av亚洲av综合av国产av| 免费观看av网站的网址| 可以免费在线观看a视频的电影网站| 一本久久精品| 国产又爽黄色视频| 大型av网站在线播放| 啦啦啦在线观看免费高清www| a级毛片黄视频| 青春草视频在线免费观看| 国产福利在线免费观看视频| 欧美精品av麻豆av| 九色亚洲精品在线播放| 欧美激情 高清一区二区三区| 亚洲人成电影观看| 电影成人av| e午夜精品久久久久久久| av国产精品久久久久影院| 国产精品一二三区在线看| 深夜精品福利| 日本av免费视频播放| av在线播放精品| 国产欧美日韩一区二区三区在线| 精品国产一区二区三区久久久樱花| 中文字幕人妻熟女乱码| 亚洲中文日韩欧美视频| 日本欧美视频一区| 免费高清在线观看视频在线观看| 男男h啪啪无遮挡| 精品人妻一区二区三区麻豆| 欧美精品人与动牲交sv欧美| 国产精品一区二区免费欧美 | 免费少妇av软件| 亚洲国产精品一区二区三区在线| 搡老乐熟女国产| 午夜两性在线视频| 成人免费观看视频高清| 大码成人一级视频| 亚洲成人免费电影在线观看 | 久久精品国产综合久久久| 人人妻,人人澡人人爽秒播 | 亚洲成人手机| 国产淫语在线视频| 亚洲精品自拍成人| 高清不卡的av网站| 国产淫语在线视频| 成年人黄色毛片网站| 久久久久久久精品精品| 伊人久久大香线蕉亚洲五| 亚洲人成网站在线观看播放| 观看av在线不卡| 9191精品国产免费久久| 久久亚洲国产成人精品v| 国产一区二区三区av在线| 一级毛片女人18水好多 | 国产精品免费大片| 成年av动漫网址| 精品免费久久久久久久清纯 | 岛国毛片在线播放| 99热网站在线观看| 99热全是精品| 老司机影院毛片| 亚洲精品国产一区二区精华液| 丰满饥渴人妻一区二区三| 大香蕉久久网| 多毛熟女@视频| 久久久久久免费高清国产稀缺| 国产成人精品在线电影| 国产精品免费大片| 老熟女久久久| 欧美精品人与动牲交sv欧美| 色播在线永久视频| 久久毛片免费看一区二区三区| 69精品国产乱码久久久| 欧美精品亚洲一区二区| 国产欧美日韩精品亚洲av| 精品久久久久久电影网| 性色av一级| 国产精品一区二区在线观看99| 亚洲精品国产一区二区精华液| 久热爱精品视频在线9| 久久亚洲精品不卡| 色网站视频免费| 别揉我奶头~嗯~啊~动态视频 | 丁香六月欧美| 欧美日韩黄片免| 蜜桃在线观看..| 欧美日韩亚洲综合一区二区三区_| 91国产中文字幕| 黑人猛操日本美女一级片| 午夜91福利影院| 自线自在国产av| 女人久久www免费人成看片| 黄频高清免费视频| 七月丁香在线播放| 亚洲欧美一区二区三区久久| 一二三四在线观看免费中文在| 午夜影院在线不卡| 大陆偷拍与自拍| 国产精品成人在线| 日韩中文字幕视频在线看片| 久久久久久久久免费视频了| 日本欧美国产在线视频| 三上悠亚av全集在线观看| 欧美黑人精品巨大| 性高湖久久久久久久久免费观看| 青草久久国产| av网站在线播放免费| 午夜福利视频精品| 老司机亚洲免费影院| 亚洲专区中文字幕在线| 1024香蕉在线观看| 免费女性裸体啪啪无遮挡网站| e午夜精品久久久久久久| 国产成人欧美在线观看 | 18禁观看日本| 成年人午夜在线观看视频| 热re99久久国产66热| netflix在线观看网站| 老司机深夜福利视频在线观看 | 久久九九热精品免费| 黄色视频不卡| 在线观看免费午夜福利视频| 免费久久久久久久精品成人欧美视频| 午夜激情av网站| 精品视频人人做人人爽| 国产av国产精品国产| 亚洲欧美日韩另类电影网站| 欧美老熟妇乱子伦牲交| 操美女的视频在线观看| av电影中文网址| 赤兔流量卡办理| 人成视频在线观看免费观看| 一本久久精品| 国产淫语在线视频| 欧美黄色淫秽网站| 下体分泌物呈黄色| 久久人妻熟女aⅴ| 大片电影免费在线观看免费| 亚洲图色成人| 成年美女黄网站色视频大全免费| 国产1区2区3区精品| 91字幕亚洲| 这个男人来自地球电影免费观看| 丰满人妻熟妇乱又伦精品不卡| 精品人妻1区二区| 少妇裸体淫交视频免费看高清 | 久久久国产欧美日韩av| 久久精品人人爽人人爽视色| 久久久久视频综合| 侵犯人妻中文字幕一二三四区| 国产精品一区二区在线观看99| 一级毛片 在线播放| 亚洲三区欧美一区| 成年人免费黄色播放视频| 高清av免费在线| 老司机在亚洲福利影院| 欧美日韩综合久久久久久| 欧美中文综合在线视频| 亚洲精品一区蜜桃| a级毛片在线看网站| 男男h啪啪无遮挡| 在线观看人妻少妇| 婷婷色综合大香蕉| 无限看片的www在线观看| www.999成人在线观看| 如日韩欧美国产精品一区二区三区| 校园人妻丝袜中文字幕| 亚洲av在线观看美女高潮| 交换朋友夫妻互换小说| 人人妻人人添人人爽欧美一区卜| 18禁观看日本| 1024视频免费在线观看| 欧美久久黑人一区二区| 国产xxxxx性猛交| 99国产精品一区二区蜜桃av | 一本一本久久a久久精品综合妖精| 欧美成狂野欧美在线观看| 久久中文字幕一级| 国产成人91sexporn| 一本久久精品| 极品人妻少妇av视频| 国产黄频视频在线观看| 大片免费播放器 马上看| 深夜精品福利| 两人在一起打扑克的视频| 老司机靠b影院| 视频区欧美日本亚洲| 亚洲情色 制服丝袜| 国产高清视频在线播放一区 | 嫁个100分男人电影在线观看 | 久久国产精品人妻蜜桃| 美国免费a级毛片| 少妇 在线观看| 免费高清在线观看日韩| www.熟女人妻精品国产| 在线观看一区二区三区激情| 亚洲精品日韩在线中文字幕| 热99国产精品久久久久久7| 久热这里只有精品99| 七月丁香在线播放| 天堂俺去俺来也www色官网| 男女免费视频国产| 一区二区三区四区激情视频| 国产精品国产av在线观看| 制服诱惑二区| 手机成人av网站| 深夜精品福利| 丰满迷人的少妇在线观看| 亚洲国产精品国产精品| 国产国语露脸激情在线看| 热99国产精品久久久久久7| 国产精品一区二区在线观看99| 日本猛色少妇xxxxx猛交久久| 1024视频免费在线观看| 久久精品国产综合久久久| avwww免费| 18禁裸乳无遮挡动漫免费视频| 国产精品一区二区精品视频观看| 国产成人精品久久二区二区91| 欧美另类一区| 久热这里只有精品99| 婷婷丁香在线五月| 一二三四在线观看免费中文在| 人妻一区二区av| 人成视频在线观看免费观看| 9热在线视频观看99| 国产不卡av网站在线观看| 777米奇影视久久| 欧美成人午夜精品| 亚洲精品日韩在线中文字幕| 国产1区2区3区精品| 精品人妻熟女毛片av久久网站| 99热国产这里只有精品6| 老汉色∧v一级毛片| 国产av精品麻豆| 啦啦啦视频在线资源免费观看| 免费少妇av软件| 成人国产av品久久久| cao死你这个sao货| 高潮久久久久久久久久久不卡| a级片在线免费高清观看视频| 99热国产这里只有精品6| 麻豆乱淫一区二区| 亚洲激情五月婷婷啪啪| 欧美人与性动交α欧美软件| 免费观看人在逋| 国产一区二区 视频在线| 日韩欧美一区视频在线观看| 18禁黄网站禁片午夜丰满| 婷婷色av中文字幕| 国产精品一国产av| 国产精品熟女久久久久浪| 中文字幕人妻丝袜制服| 午夜久久久在线观看| 激情视频va一区二区三区| 国产一区二区三区综合在线观看| 久久青草综合色| 涩涩av久久男人的天堂| 18禁黄网站禁片午夜丰满| 国语对白做爰xxxⅹ性视频网站| 国产在线观看jvid| 久久久久久亚洲精品国产蜜桃av| 少妇粗大呻吟视频| 国产精品一国产av| 国产精品99久久99久久久不卡| 久久99一区二区三区| 午夜免费鲁丝| 免费在线观看完整版高清| 最新在线观看一区二区三区 | 亚洲,欧美精品.| 亚洲av男天堂| 国产成人欧美| 1024视频免费在线观看| 曰老女人黄片| 日韩av不卡免费在线播放| 丝袜美足系列| 你懂的网址亚洲精品在线观看| 久久热在线av| 午夜福利,免费看| 国产成人av教育| 另类精品久久| 免费观看av网站的网址| 亚洲国产精品成人久久小说| 国产在线视频一区二区| 精品高清国产在线一区| 久久精品国产a三级三级三级| 婷婷色综合大香蕉| 男人舔女人的私密视频| 18在线观看网站| 日韩中文字幕欧美一区二区 | 欧美日韩精品网址| 久久女婷五月综合色啪小说| 亚洲国产最新在线播放| a级片在线免费高清观看视频| 久久国产精品人妻蜜桃| 精品视频人人做人人爽| 欧美另类一区| 另类亚洲欧美激情| 亚洲国产欧美在线一区| 国产主播在线观看一区二区 | 免费少妇av软件| 一区二区三区精品91| 亚洲熟女毛片儿| 亚洲av日韩精品久久久久久密 | 丝袜在线中文字幕| 纯流量卡能插随身wifi吗| 99国产精品免费福利视频| 国产麻豆69| 亚洲成人手机| 亚洲精品一二三| 99国产精品一区二区蜜桃av | 国产亚洲欧美在线一区二区| 国产精品人妻久久久影院| 91九色精品人成在线观看| 老熟女久久久| 久久久亚洲精品成人影院| 十分钟在线观看高清视频www| www.熟女人妻精品国产| 91九色精品人成在线观看| 国产精品一二三区在线看| 另类亚洲欧美激情| 精品熟女少妇八av免费久了| 男女国产视频网站| 国产精品99久久99久久久不卡| 久久久久精品人妻al黑| 丝袜人妻中文字幕| 真人做人爱边吃奶动态| 国产无遮挡羞羞视频在线观看| 久久久欧美国产精品| 国产麻豆69| 老司机影院成人| 午夜福利影视在线免费观看| 免费一级毛片在线播放高清视频 | 久久久久国产精品人妻一区二区| 男人爽女人下面视频在线观看| 国产成人欧美| 男女床上黄色一级片免费看| 欧美日韩亚洲综合一区二区三区_| 日韩大码丰满熟妇| 久久久国产一区二区| 超碰成人久久| 亚洲国产毛片av蜜桃av| 成年人黄色毛片网站| 日韩av免费高清视频| 亚洲欧美精品综合一区二区三区| 人妻一区二区av| 亚洲色图 男人天堂 中文字幕| 在现免费观看毛片| 超色免费av| 午夜两性在线视频| 亚洲一区二区三区欧美精品| 999精品在线视频| 母亲3免费完整高清在线观看| 一区二区三区激情视频| 久久久久久久久久久久大奶| 国产成人欧美| 日韩av在线免费看完整版不卡| 国产亚洲一区二区精品| 纵有疾风起免费观看全集完整版| 热99国产精品久久久久久7| 久久人人爽av亚洲精品天堂| 午夜福利一区二区在线看| 国产亚洲午夜精品一区二区久久| 国产精品香港三级国产av潘金莲 | 日韩制服丝袜自拍偷拍| 一本—道久久a久久精品蜜桃钙片| 少妇粗大呻吟视频| 可以免费在线观看a视频的电影网站| 国产成人欧美| 女人精品久久久久毛片| 国产深夜福利视频在线观看| 七月丁香在线播放| 精品人妻在线不人妻| 少妇的丰满在线观看| 另类亚洲欧美激情| 2018国产大陆天天弄谢| 精品人妻熟女毛片av久久网站| 久久久国产一区二区| 尾随美女入室| 美女福利国产在线| av有码第一页| a级毛片黄视频| 性色av乱码一区二区三区2| 精品少妇一区二区三区视频日本电影| 久久精品国产亚洲av高清一级| 国产精品国产三级国产专区5o| 亚洲人成电影免费在线| 91字幕亚洲| 又粗又硬又长又爽又黄的视频| 涩涩av久久男人的天堂| 国产人伦9x9x在线观看| 久久亚洲精品不卡| 啦啦啦视频在线资源免费观看| 国产有黄有色有爽视频| 亚洲av在线观看美女高潮| 亚洲熟女精品中文字幕| 国产精品久久久av美女十八| svipshipincom国产片| 80岁老熟妇乱子伦牲交| 亚洲 欧美一区二区三区| 国产伦理片在线播放av一区| 欧美乱码精品一区二区三区| 宅男免费午夜| 国产精品一区二区免费欧美 | 99国产综合亚洲精品| 久久热在线av| 黑丝袜美女国产一区| 亚洲欧美中文字幕日韩二区| 久久免费观看电影| 丝瓜视频免费看黄片| 亚洲美女黄色视频免费看| 99久久人妻综合| 国产男人的电影天堂91| 欧美黑人精品巨大| 真人做人爱边吃奶动态| 色94色欧美一区二区| a 毛片基地| 一级毛片黄色毛片免费观看视频| 一边摸一边抽搐一进一出视频| 日本色播在线视频| 蜜桃国产av成人99| 大码成人一级视频| 婷婷色av中文字幕| 日韩中文字幕视频在线看片| 在线精品无人区一区二区三| 看免费av毛片| 国产女主播在线喷水免费视频网站| 久久人人爽人人片av| 亚洲av国产av综合av卡| 亚洲精品一区蜜桃| 成人国产av品久久久| 日日夜夜操网爽| 国产精品二区激情视频| 一本大道久久a久久精品| 免费女性裸体啪啪无遮挡网站| 丰满饥渴人妻一区二区三| 欧美日韩视频精品一区| 男人舔女人的私密视频| 丰满人妻熟妇乱又伦精品不卡| av视频免费观看在线观看| 中文字幕av电影在线播放| 五月开心婷婷网| 一本久久精品| 午夜av观看不卡| 妹子高潮喷水视频| 桃花免费在线播放| 天天躁夜夜躁狠狠躁躁| 欧美中文综合在线视频| 波野结衣二区三区在线| 女人爽到高潮嗷嗷叫在线视频| 只有这里有精品99| 欧美日韩国产mv在线观看视频| 国产成人影院久久av| 丁香六月天网| 777米奇影视久久| 黄色视频不卡| 国产视频首页在线观看| av有码第一页| 精品少妇一区二区三区视频日本电影| 一区二区三区精品91| 亚洲国产毛片av蜜桃av| 50天的宝宝边吃奶边哭怎么回事| 亚洲中文av在线| 一级毛片我不卡| 黄频高清免费视频| 精品少妇内射三级| 欧美激情 高清一区二区三区| 国产高清国产精品国产三级| 免费观看人在逋| 久久久久久免费高清国产稀缺| 亚洲天堂av无毛| 亚洲精品久久成人aⅴ小说| 大片电影免费在线观看免费| 成人免费观看视频高清| 狂野欧美激情性bbbbbb| 香蕉国产在线看| 天天影视国产精品| 午夜激情久久久久久久| 欧美激情极品国产一区二区三区| 一级a爱视频在线免费观看| 久久久亚洲精品成人影院| 精品人妻一区二区三区麻豆| 国产欧美日韩综合在线一区二区| 男女无遮挡免费网站观看| 国产免费视频播放在线视频| 看免费成人av毛片| 又紧又爽又黄一区二区| 99热网站在线观看| 国产xxxxx性猛交| 中文字幕高清在线视频| cao死你这个sao货| 丰满饥渴人妻一区二区三| 午夜激情久久久久久久| av在线老鸭窝| 啦啦啦视频在线资源免费观看| 亚洲一区二区三区欧美精品| 超碰97精品在线观看| 亚洲国产中文字幕在线视频| 午夜福利免费观看在线| 丰满少妇做爰视频| 交换朋友夫妻互换小说| 亚洲精品在线美女| 又大又爽又粗| 999精品在线视频| 亚洲欧洲国产日韩| 免费在线观看完整版高清| 免费日韩欧美在线观看| 日韩av在线免费看完整版不卡| 丝袜喷水一区| 国产97色在线日韩免费| 91麻豆精品激情在线观看国产 | 大片免费播放器 马上看| 欧美亚洲 丝袜 人妻 在线| 老司机靠b影院| 亚洲精品国产色婷婷电影| 国产亚洲精品久久久久5区| 国产精品一区二区精品视频观看| 久久久久久亚洲精品国产蜜桃av| 精品少妇黑人巨大在线播放| 91九色精品人成在线观看| 99香蕉大伊视频| 热99久久久久精品小说推荐| 丁香六月欧美|