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

    考慮尾礦材料參數(shù)空間變異性的尾礦壩穩(wěn)定可靠度分析*

    2020-01-13 03:11:18蔣水華朱明明黃勁松
    關(guān)鍵詞:變異性尾礦安全系數(shù)

    蔣水華,朱明明,黃勁松

    (南昌大學(xué) 建筑工程學(xué)院,江西 南昌 330031)

    0 引言

    隨著國民經(jīng)濟(jì)和工程技術(shù)的快速發(fā)展,我國尾礦庫數(shù)量越來越多,規(guī)模越來越大。目前我國尾礦庫超過了12 000座,其中95%的尾礦庫采用上游法筑壩,尾礦庫同時也是礦業(yè)風(fēng)險的主要來源之一[1]。近年來,地震、洪水等各種原因?qū)е挛覈驳V庫潰壩事故頻繁發(fā)生,危害極其嚴(yán)重。另外,尾礦材料是性質(zhì)極其復(fù)雜的地質(zhì)介質(zhì),在施工過程中人為與非人為因素的作用下,其物質(zhì)組成和內(nèi)部結(jié)構(gòu)會發(fā)生一定的變化,尾礦材料性質(zhì)在空間上相差較大,即呈現(xiàn)一定的空間變異性。尾礦材料參數(shù)空間變異性的客觀存在給尾礦壩穩(wěn)定性分析帶來了一定的困難,常用的確定性分析方法因不能反映這種空間變異性的影響,會導(dǎo)致評價結(jié)果存在一定的偏差。

    目前針對尾礦壩穩(wěn)定可靠度問題開展了大量有益的研究工作,如胡平安等[1]采用JC法探究了不同尾礦材料力學(xué)參數(shù)的變異性對湖北省漁潭尾礦壩可靠度的影響;劉迪等[2]利用貝葉斯網(wǎng)絡(luò)模型建立尾礦壩穩(wěn)定性因果關(guān)聯(lián)分析模型,分析了各因素變化對尾礦壩穩(wěn)定的影響;曹志松等[3]采用二次多項式構(gòu)建響應(yīng)面模型,計算了黏聚力、內(nèi)摩擦角和重度服從正態(tài)分布情況下的尾礦壩穩(wěn)定可靠度;張揚等[4]采用蒙特卡羅模擬(MCS)方法探討了尾礦材料黏聚力、內(nèi)摩擦角、密度等和孔隙水壓力的不確定性對獨木垅尾礦壩抗滑穩(wěn)定的影響。雖然這些研究工作推動了可靠度理論在尾礦庫安全性評價中的應(yīng)用,但是大多采用隨機(jī)變量模型模擬尾礦材料參數(shù)的不確定性,不能有效考慮空間不同點局部和整體尾礦材料物理力學(xué)性質(zhì)的差異,無法客觀地評價尾礦材料參數(shù)空間變異性對尾礦壩穩(wěn)定可靠度的影響。相比之下,隨機(jī)場理論將土層不同位置上的土壤參數(shù)模擬為一個服從某一概率分布的隨機(jī)變量,并與相鄰位置上的隨機(jī)變量相關(guān),非常適合表征尾礦材料的空間變異性,目前在邊坡工程中得到了廣泛應(yīng)用[5-7],但是較少應(yīng)用到尾礦材料參數(shù)空間變異性表征中。此外,對于需考慮尾礦材料參數(shù)空間變異性的復(fù)雜尾礦壩穩(wěn)定問題,MCS等傳統(tǒng)可靠度分析方法的計算精度和效率欠佳,因此需要發(fā)展一種高效的考慮參數(shù)空間變異性的尾礦壩穩(wěn)定可靠度分析方法。

    本文結(jié)合隨機(jī)場理論、代理模型和結(jié)構(gòu)可靠度理論,采用考慮尾礦材料參數(shù)空間變異性尾礦壩可靠度分析的非侵入式隨機(jī)有限元法,通過Karhunen-Loève(K-L)展開方法模擬尾礦材料物理力學(xué)參數(shù)空間變異性,借助Hermite隨機(jī)多項式展開建立尾礦壩安全系數(shù)代理模型,在此基礎(chǔ)上計算尾礦壩邊坡穩(wěn)定可靠度,并以實際尾礦壩工程為例驗證了提出方法的有效性[5-7]。

    1 尾礦材料參數(shù)空間變異性模擬

    采用隨機(jī)場模擬尾礦材料參數(shù)空間變異性時,重要的一步是進(jìn)行隨機(jī)場離散,本文選用計算精度和效率較高的K-L級數(shù)展開方法離散尾礦材料參數(shù)隨機(jī)場。

    以離散尾礦滲透系數(shù)ks與內(nèi)摩擦角φ各向異性對數(shù)正態(tài)參數(shù)隨機(jī)場為例,簡要介紹基于K-L級數(shù)展開方法的尾礦材料參數(shù)空間變異性模擬過程。K-L展開方法離散隨機(jī)場是以譜分解有界,對稱且正定的相關(guān)函數(shù)為基礎(chǔ),將土體參數(shù)隨機(jī)場的離散轉(zhuǎn)化為求解Fredholm積分方程的特征值問題[6],即

    (1)

    式中:(x1,y1)和(x2,y2)表示計算區(qū)域Ω中的任意2個坐標(biāo)點;ρ[(x1,y1),(x2,y2)]為任意2點(x1,y1),(x2,y2)處隨機(jī)場特性值之間的自相關(guān)系數(shù);fi(.)和λi分別為自相關(guān)函數(shù)對應(yīng)的特征函數(shù)和特征值。

    采用描述參數(shù)空間自相關(guān)性的高斯型自相關(guān)函數(shù)所模擬的參數(shù)隨機(jī)場分布平滑度和連續(xù)性較好[8],故本文選取高斯型自相關(guān)函數(shù),表達(dá)式為:

    (2)

    式中:lh和lv分別表示水平和垂直方向上的相關(guān)距離,用于表征參數(shù)空間變異性的自相關(guān)程度,相關(guān)距離越大對應(yīng)的尾礦材料參數(shù)的空間自相關(guān)性越強(qiáng)。

    根據(jù)式(2)可計算得到原始空間不同點處的隨機(jī)場特性值之間的自相關(guān)系數(shù)。一般水平相關(guān)距離是垂直相關(guān)距離的10倍左右,水平相關(guān)距離在10~40 m范圍內(nèi),垂直相關(guān)距離變化范圍為1~3 m[9]。

    當(dāng)采用式(2)高斯型自相關(guān)函數(shù)時,式(1)沒有解析解,為此本文采用wavelet-Galerkin技術(shù)[10]進(jìn)行數(shù)值求解?;谧韵嚓P(guān)函數(shù)的特征解可將隨機(jī)場H(x,y)展開為:

    (3)

    式中:(x,y)為隨機(jī)場區(qū)域中的任意坐標(biāo)點,(x,y)∈Ω;ξXi,j為參數(shù)隨機(jī)場離散的獨立標(biāo)準(zhǔn)正態(tài)隨機(jī)變量;μlnXi和σlnXi分別為相應(yīng)參數(shù)隨機(jī)場lnXi的均值和標(biāo)準(zhǔn)差;n為截斷項數(shù),與隨機(jī)場離散的隨機(jī)變量數(shù)目對應(yīng)。截斷項數(shù)n取決于精度要求和自相關(guān)函數(shù)的形式,文獻(xiàn)[11]~[12]建議采用隨機(jī)場期望能比率因子ε≥95%作為確定截斷項數(shù)n取值的依據(jù)。

    2 尾礦壩可靠度非侵入式隨機(jī)分析方法

    通過K-L展開方法得到參數(shù)隨機(jī)場特性值后,需要進(jìn)行含大量隨機(jī)變量的尾礦壩邊坡可靠度分析。尾礦壩滲流分析和邊坡穩(wěn)定一般借助有限元軟件實現(xiàn),因此尾礦壩安全系數(shù)與尾礦材料參數(shù)之間是非線性隱式函數(shù)關(guān)系。傳統(tǒng)的MCS、一次二階矩等可靠度方法不能有效求解含大量隨機(jī)變量的極限狀態(tài)函數(shù)是隱式的尾礦壩可靠度問題。為此,本文將近幾年發(fā)展起來的非侵入式隨機(jī)有限元法[8]拓展到尾礦庫工程領(lǐng)域,采用Hermite隨機(jī)多項式展開(Hermite Polynomial Chaos Expansion,PCE)代替尾礦壩邊坡安全系數(shù)FS與尾礦材料參數(shù)ξ之間的非線性隱式函數(shù)關(guān)系:

    (4)

    式中:n為隨機(jī)變量的數(shù)目;ξ=(ξ1,ξ2,…,ξn)T為獨立標(biāo)準(zhǔn)隨機(jī)向量,與式(1)參數(shù)隨機(jī)場離散的隨機(jī)變量一一對應(yīng);Γp(.)表示為p階Hermite隨機(jī)多項式展開;a0,ai1,ai1i2,ai1i2i3,…為待定系數(shù),其數(shù)目為(n+p)!/(n!p!)。

    接著,采用拉丁超立方抽樣技術(shù)(Latin Hypercube Sampling,LHS)[13]產(chǎn)生輸入?yún)?shù)隨機(jī)樣本點,代入有限元程序或巖土商業(yè)軟件計算安全系數(shù),再根據(jù)式(4)建立線性代數(shù)方程組求解多項式展開系數(shù)。一旦確定了尾礦壩安全系數(shù)和輸入?yún)?shù)ξ的顯式函數(shù)關(guān)系,即安全系數(shù)代理模型,便可建立新的顯式函數(shù)表達(dá)的尾礦壩穩(wěn)定可靠度分析極限狀態(tài)函數(shù):

    G(ξ)=FS(ξ)-1.0

    (5)

    基于采用Hermite隨機(jī)多項式展開建立的代理模型,利用100 000次MCS方法計算尾礦壩失穩(wěn)破壞概率。顯然,該方法不僅實現(xiàn)了可靠度分析和確定性有限元分析的不耦合,而且直接使用代理模型計算安全系數(shù)FS,無需進(jìn)行大量的有限元計算,提高了計算效率。

    3 工程應(yīng)用

    以西木尾礦庫工程[14]為例,探討尾礦壩材料滲透系數(shù)和內(nèi)摩擦角空間變異性對尾礦壩穩(wěn)定可靠度的影響。西木礦區(qū)位于加拿大魁北克省魯恩-諾蘭達(dá)以東約40 km的鮑斯凱鎮(zhèn),于2004年被發(fā)現(xiàn)。西木尾礦壩壩高15 m,大壩在建設(shè)結(jié)束時的幾何形狀,初始模型中各層材料分布、層厚和層寬如圖1所示。各個土層的物理力學(xué)參數(shù)見表1。土層下限是基巖和過密的冰磧層,在大壩建設(shè)之后保持平衡,被認(rèn)為是穩(wěn)定地層。

    圖1 尾礦壩材料剖面及有限元網(wǎng)格Fig.1 Material profile and finite element mesh of the tailing dam

    表1 尾礦材料物理力學(xué)參數(shù)Table 1 Physical and mechanical parameters of tailings materials

    注:γ,φ,c和ks分別為尾礦材料的重度、內(nèi)摩擦角、黏聚力和滲透系數(shù);a,n和m為土水特征曲線擬合參數(shù)。

    3.1 有限元模型及計算參數(shù)

    西木尾礦壩有限元模型如圖1所示,一共劃分為10 685個節(jié)點和10 407個網(wǎng)格大小為0.5 m的四邊形和三角形混合單元。邊界條件為:上游水位(左側(cè))13 m,下游水位(右側(cè))8 m,并在下游邊坡(右側(cè))設(shè)置了潛在滲流面?;赟EEP/W模塊對西木尾礦壩進(jìn)行飽和滲流分析,再將得到浸潤線和孔隙水壓分布等滲流結(jié)果導(dǎo)入到SLOPE/W模塊,采用摩根斯坦-普萊斯法計算的安全系數(shù)FS=1.931,最危險滑動面如圖2所示。所計算的安全系數(shù)和臨界滑移面與Coulibaly等[14]計算的安全系數(shù)1.967和臨界滑移面基本一致,說明了該確定性有限元分析的有效性。

    因受基質(zhì)吸力、含水量、滲透率等因素的影響,因此對西木尾礦壩進(jìn)行飽和-非飽和滲流分析十分必要。本文將填石和尾細(xì)砂視作非飽和材料,采用目前應(yīng)用廣泛的Fredlund-Xing[15]模型表征的土水特征曲線(SWCC)模擬體積含水量與基質(zhì)吸力之間的函數(shù)關(guān)系。模型擬合參數(shù)a為進(jìn)氣值相關(guān)的土性參數(shù);n為與SWCC斜率相關(guān)的土性參數(shù);m為與殘余含水率相關(guān)的土性參數(shù)[16],參數(shù)取值見表1。圖3(a)和3(b)分別給出了尾細(xì)砂土水特征曲線與滲透系數(shù)函數(shù)。將飽和-非飽和滲流分析結(jié)果導(dǎo)入到SLOPE/W模塊,采用摩根斯坦-普萊斯法計算得到西木尾礦壩安全系數(shù)FS=1.944,大于在所有尾礦材料為飽和狀態(tài)下計算的安全系數(shù)。

    圖2 尾礦壩穩(wěn)定性分析結(jié)果Fig.2 Stability analysis results of the tailings dam

    圖3 尾細(xì)砂土水特征與滲透系數(shù)函數(shù)曲線Fig.3 Soil-water characteristic curve and permeability coefficient function curve of tailings fine sand

    另外,尾礦壩上的局部坍塌和滑移大多是由于地震導(dǎo)致尾部或地基土液化引起的[14]。為此,本文考慮Ⅶ級地震作用,進(jìn)一步對該尾礦壩進(jìn)行穩(wěn)定可靠度分析。根據(jù)建設(shè)部頒布的《關(guān)于統(tǒng)一抗震設(shè)計規(guī)范地面運動加速度設(shè)計取值的通知》規(guī)定,Ⅶ級地震設(shè)計基本地震加速度取0.1 g,豎向地震影響系數(shù)的最大值取水平地震影響系數(shù)最大值的65%[17]。采用擬靜力法模擬Ⅶ級地震作用,同樣通過非飽和滲流穩(wěn)定分析得到安全系數(shù)FS=1.418,圖4給出了最危險滑面位置。根據(jù)尾礦設(shè)施設(shè)計規(guī)范[18]可知,該尾礦壩在Ⅶ級地震作用下仍處于穩(wěn)定狀態(tài)。

    圖4 地震作用下尾礦壩穩(wěn)定性分析結(jié)果Fig.4 Stability analysis results of tailings dam under effect of earthquake

    3.2 尾礦壩穩(wěn)定可靠度分析

    如前所述,受到尾礦壩長期多循環(huán)水力充填以及固結(jié)沉降的作用,尾礦材料強(qiáng)度參數(shù)存在較大的離散性和空間變異性,與尾礦材料抗剪強(qiáng)度參數(shù)空間變異性相比,尾礦材料重度空間變異性較小,一般小于10%,可視作常量。故本例僅考慮尾礦材料滲透系數(shù)ks和內(nèi)摩擦角φ的空間變異性,并假設(shè)尾細(xì)砂、粉土、粉質(zhì)黏土和填石的ks,φ均服從對數(shù)正態(tài)分布,尾礦壩材料參數(shù)統(tǒng)計特征見表2。

    表2 尾礦材料參數(shù)統(tǒng)計特征Table 2 Statistical characteristics of tailings materials parameters

    圖5 尾細(xì)砂滲透系數(shù)隨機(jī)場的典型實現(xiàn)Fig.5 Typical realization of random field for permeability coefficient of tailings fine sand

    在此基礎(chǔ)上,利用LHS法對輸入的變量進(jìn)行抽樣,根據(jù)參數(shù)統(tǒng)計特征生成隨機(jī)變量和隨機(jī)場數(shù)據(jù),分別賦到尾礦壩模型,然后通過滲流有限元和穩(wěn)定分析得到FS,再采用Hermite展開多項式建立FS代理模型(即FS與隨機(jī)變量和隨機(jī)場離散的隨機(jī)變量之間的顯式函數(shù)關(guān)系)。圖6比較了由不同方法計算的極限狀態(tài)函數(shù)累積分布函數(shù)(CDF)曲線。由圖6可見,分別進(jìn)行2 000,4 000和10 000次滲流有限元及穩(wěn)定性分析構(gòu)建FS代理模型,對應(yīng)的二階PCE+MCS方法計算的極限狀態(tài)函數(shù)的CDF曲線非常吻合,并且與10 000次直接LHS方法的計算結(jié)果僅僅在低概率區(qū)域存在微小差別。其原因是10 000次直接LHS方法對于低失穩(wěn)破壞概率(<10-3)水平尾礦壩可靠度問題,計算精度不夠。

    圖6 極限狀態(tài)函數(shù)累積分布曲線的比較Fig.6 Comparison on cumulative distribution curves of limit state functions

    由圖6極限狀態(tài)函數(shù)的CDF曲線可以直接計算尾礦壩失穩(wěn)破壞概率。表3比較了不同方法計算的尾礦壩失穩(wěn)破壞概率。從表3可以看出,西木尾礦壩失穩(wěn)破壞概率水平很低,二階PCE+MCS方法(Np=2 000),二階PCE+MCS方法(Np=4 000)和二階PCE+MCS(Np=10 000)方法計算的破壞概率分別為8.6×10-4,8.1×10-4,7.6×10-4,可見3種方法計算的破壞概率非常吻合,和10 000次直接LHS方法計算的9.9×10-4也基本一致,說明了本文非侵入式隨機(jī)有限元法的有效性。此外,與10 000次直接LHS方法相比,二階PCE+MCS方法只需要進(jìn)行2 000次尾礦壩滲流有限元和穩(wěn)定性分析就可以滿足計算精度要求,計算量是直接LHS方法的1/4,可見本文提出的方法具有較高的計算效率,為解決考慮多個尾礦材料參數(shù)空間變異性的低概率水平尾礦壩穩(wěn)定可靠度問題提供了一條有效路徑。

    表3 尾礦壩可靠度分析結(jié)果Table 3 Reliability analysis results of the tailings dam

    4 結(jié)論

    1)采用隨機(jī)場理論模擬西木尾礦材料參數(shù)空間變異性,采用Hermite隨機(jī)多項式建立了尾礦壩安全系數(shù)與不確定性輸入?yún)?shù)之間的代理模型,將尾礦壩滲流有限元及穩(wěn)定性分析視作黑箱子直接調(diào)用,實現(xiàn)了尾礦壩可靠度分析與確定性穩(wěn)定分析的不耦合。

    2)與直接LHS方法相比,提出方法滿足計算精度要求的同時,計算量減少了80%,為解決考慮多個尾礦材料參數(shù)空間變異性的低概率水平尾礦壩穩(wěn)定可靠度問題提供了一條有效的路徑。

    3)盡管西木尾礦壩在正常和地震工況下的安全系數(shù)很高,失穩(wěn)破壞概率低至10-4量級,與尾礦壩設(shè)計規(guī)范給出的允許安全系數(shù)和破壞概率相比,該尾礦壩在Ⅶ級地震作用下發(fā)生失穩(wěn)破壞的可能性較小,但是仍然需要設(shè)計一些簡單的尾礦壩安全防護(hù)抗震措施。

    猜你喜歡
    變異性尾礦安全系數(shù)
    基于視覺識別的浮選尾礦在線測灰儀的應(yīng)用
    考慮材料性能分散性的航空發(fā)動機(jī)結(jié)構(gòu)安全系數(shù)確定方法
    咳嗽變異性哮喘的預(yù)防和治療
    鐵尾礦資源的研究與應(yīng)用
    江西建材(2018年4期)2018-04-10 12:36:50
    重力式擋土墻抗滑穩(wěn)定性安全系數(shù)的異性分析及經(jīng)驗安全系數(shù)方法
    閘室樁基處理后水平抗滑穩(wěn)定安全系數(shù)提高值的估算范圍研究
    接近物體感測庫顯著提升安全系數(shù)
    汽車零部件(2014年6期)2014-09-20 06:29:36
    咳嗽變異性哮喘的中醫(yī)治療近況
    清肺止咳湯治療咳嗽變異性哮喘40例
    冬病夏治止咳貼貼敷治療小兒咳嗽變異性哮喘40例
    亚洲激情五月婷婷啪啪| 十八禁网站网址无遮挡 | 久久久久国产网址| 18禁动态无遮挡网站| 人人妻人人添人人爽欧美一区卜| 美女主播在线视频| 国产精品久久久久久久电影| 国产在线男女| 下体分泌物呈黄色| 国产国拍精品亚洲av在线观看| 亚洲欧美精品专区久久| 不卡视频在线观看欧美| 久久精品夜色国产| 国产精品免费大片| 国产精品一区二区在线观看99| 精品99又大又爽又粗少妇毛片| 亚洲高清免费不卡视频| 高清av免费在线| 亚洲精品乱码久久久久久按摩| 亚洲精品一二三| 狠狠精品人妻久久久久久综合| 亚洲婷婷狠狠爱综合网| 国产熟女午夜一区二区三区 | 久热久热在线精品观看| 国精品久久久久久国模美| 又大又黄又爽视频免费| 大陆偷拍与自拍| 午夜影院在线不卡| 免费在线观看成人毛片| av福利片在线| 亚洲精华国产精华液的使用体验| 国产淫片久久久久久久久| 在线观看国产h片| 亚洲美女黄色视频免费看| 伦理电影大哥的女人| 欧美精品高潮呻吟av久久| 黄色视频在线播放观看不卡| 青春草亚洲视频在线观看| 91aial.com中文字幕在线观看| 青春草视频在线免费观看| 国产精品99久久99久久久不卡 | 自拍偷自拍亚洲精品老妇| .国产精品久久| 久久av网站| 高清毛片免费看| 在线 av 中文字幕| 天美传媒精品一区二区| 最近最新中文字幕免费大全7| 亚洲一区二区三区欧美精品| 久久久久国产网址| 国产精品不卡视频一区二区| av黄色大香蕉| 国产精品一二三区在线看| 久久午夜福利片| 午夜免费鲁丝| 精品午夜福利在线看| 国产一区二区在线观看日韩| 日韩欧美 国产精品| 夫妻性生交免费视频一级片| 日韩在线高清观看一区二区三区| av在线老鸭窝| 久热久热在线精品观看| 国产成人aa在线观看| 91久久精品国产一区二区成人| 国产精品一区二区在线不卡| 99久久中文字幕三级久久日本| 国产精品福利在线免费观看| 婷婷色麻豆天堂久久| 成人免费观看视频高清| 久久毛片免费看一区二区三区| 久久久国产欧美日韩av| 国产精品久久久久久精品古装| 午夜福利影视在线免费观看| 久久ye,这里只有精品| 又大又黄又爽视频免费| 午夜日本视频在线| 欧美区成人在线视频| 久久久久久久久久成人| 三上悠亚av全集在线观看 | 人人妻人人澡人人爽人人夜夜| 丝袜喷水一区| 国产精品一区www在线观看| 毛片一级片免费看久久久久| 成人漫画全彩无遮挡| 中文精品一卡2卡3卡4更新| 日韩伦理黄色片| 成人毛片a级毛片在线播放| 精品人妻熟女毛片av久久网站| 中文字幕人妻熟人妻熟丝袜美| 亚洲国产欧美日韩在线播放 | 日本黄色日本黄色录像| 秋霞伦理黄片| 亚洲精品乱码久久久v下载方式| 亚洲精品成人av观看孕妇| 最近手机中文字幕大全| 久久国产亚洲av麻豆专区| 2018国产大陆天天弄谢| 欧美日韩视频精品一区| 国产免费福利视频在线观看| 亚洲欧美精品自产自拍| 日本黄色日本黄色录像| 尾随美女入室| 午夜免费观看性视频| 久久久久久久亚洲中文字幕| 国产深夜福利视频在线观看| 免费在线观看成人毛片| 久久久久精品性色| 亚洲精品国产av成人精品| 国精品久久久久久国模美| 亚洲av电影在线观看一区二区三区| 久久人人爽人人爽人人片va| 又大又黄又爽视频免费| 国产熟女午夜一区二区三区 | 美女主播在线视频| 青春草视频在线免费观看| 美女视频免费永久观看网站| 一级黄片播放器| 中文乱码字字幕精品一区二区三区| 亚洲成人手机| 一级二级三级毛片免费看| 国产亚洲午夜精品一区二区久久| 亚洲av成人精品一二三区| 欧美日韩在线观看h| 亚洲激情五月婷婷啪啪| 男女免费视频国产| 人妻制服诱惑在线中文字幕| 日日摸夜夜添夜夜添av毛片| 观看免费一级毛片| 高清在线视频一区二区三区| 人妻系列 视频| a 毛片基地| av女优亚洲男人天堂| 亚洲色图综合在线观看| av一本久久久久| 欧美精品国产亚洲| 国产一区亚洲一区在线观看| 深夜a级毛片| 国产精品福利在线免费观看| 最新的欧美精品一区二区| 亚洲av国产av综合av卡| 中文天堂在线官网| 纵有疾风起免费观看全集完整版| 欧美三级亚洲精品| 男女免费视频国产| 99精国产麻豆久久婷婷| 老司机影院毛片| 最近中文字幕2019免费版| 欧美日韩一区二区视频在线观看视频在线| 中国美白少妇内射xxxbb| 国产午夜精品一二区理论片| 国产成人午夜福利电影在线观看| 日韩精品免费视频一区二区三区 | 中文资源天堂在线| 97超碰精品成人国产| 亚洲三级黄色毛片| 成人美女网站在线观看视频| av线在线观看网站| 又黄又爽又刺激的免费视频.| 在线观看免费视频网站a站| 男人添女人高潮全过程视频| 午夜影院在线不卡| 又爽又黄a免费视频| av黄色大香蕉| 国产一区有黄有色的免费视频| 亚洲三级黄色毛片| 99久久精品一区二区三区| 免费少妇av软件| 国产免费又黄又爽又色| 97在线人人人人妻| 国产亚洲av片在线观看秒播厂| 久久久久久久久久久丰满| 一级爰片在线观看| 丰满饥渴人妻一区二区三| 亚洲国产精品999| 久久av网站| 日本黄色日本黄色录像| 日韩伦理黄色片| av.在线天堂| 亚洲av欧美aⅴ国产| 久久久午夜欧美精品| 欧美日韩在线观看h| 国内揄拍国产精品人妻在线| 国产一区二区在线观看av| 最近最新中文字幕免费大全7| 啦啦啦在线观看免费高清www| 男女边吃奶边做爰视频| 成人18禁高潮啪啪吃奶动态图 | 美女脱内裤让男人舔精品视频| 国产av精品麻豆| 亚洲欧美精品自产自拍| 天天躁夜夜躁狠狠久久av| 亚洲真实伦在线观看| 99热国产这里只有精品6| 成人国产av品久久久| 国产一区二区三区综合在线观看 | 18禁裸乳无遮挡动漫免费视频| 一级a做视频免费观看| 91久久精品国产一区二区三区| 青青草视频在线视频观看| 校园人妻丝袜中文字幕| 国产精品三级大全| 欧美日韩在线观看h| 日本av免费视频播放| 最黄视频免费看| 日韩三级伦理在线观看| 国产国拍精品亚洲av在线观看| 噜噜噜噜噜久久久久久91| 国产免费视频播放在线视频| 人妻少妇偷人精品九色| 免费看日本二区| av在线app专区| 国产淫语在线视频| 人妻一区二区av| 国产成人精品福利久久| av.在线天堂| 男女边摸边吃奶| 哪个播放器可以免费观看大片| 99久久中文字幕三级久久日本| 97在线视频观看| 成人18禁高潮啪啪吃奶动态图 | 亚洲av日韩在线播放| 中文欧美无线码| 看非洲黑人一级黄片| 五月玫瑰六月丁香| 久久午夜综合久久蜜桃| 免费观看无遮挡的男女| 在线 av 中文字幕| 欧美丝袜亚洲另类| 久久久a久久爽久久v久久| 中国美白少妇内射xxxbb| 五月开心婷婷网| 永久网站在线| 街头女战士在线观看网站| 国产成人91sexporn| 少妇的逼好多水| 精品少妇久久久久久888优播| 亚洲高清免费不卡视频| 老司机亚洲免费影院| a级毛色黄片| 日本91视频免费播放| 18+在线观看网站| 亚洲成色77777| a级毛片在线看网站| 久久免费观看电影| av福利片在线| 一级片'在线观看视频| 成人18禁高潮啪啪吃奶动态图 | 国产成人精品久久久久久| 午夜日本视频在线| 亚洲国产精品专区欧美| 欧美少妇被猛烈插入视频| 久久97久久精品| 亚洲人与动物交配视频| 中文字幕制服av| 成人美女网站在线观看视频| 亚洲欧美成人精品一区二区| 精品国产一区二区三区久久久樱花| 午夜日本视频在线| 日本wwww免费看| 国产男女内射视频| 免费观看性生交大片5| 国产精品熟女久久久久浪| 亚洲成人一二三区av| 免费黄频网站在线观看国产| 尾随美女入室| 天堂中文最新版在线下载| 久久99蜜桃精品久久| 免费不卡的大黄色大毛片视频在线观看| 国产亚洲午夜精品一区二区久久| 天堂8中文在线网| 一级毛片我不卡| 亚洲精品456在线播放app| 男女无遮挡免费网站观看| 十八禁网站网址无遮挡 | 一二三四中文在线观看免费高清| 成人无遮挡网站| 秋霞在线观看毛片| 国产色爽女视频免费观看| 亚洲人成网站在线观看播放| 黄色一级大片看看| 交换朋友夫妻互换小说| 少妇人妻一区二区三区视频| 青春草亚洲视频在线观看| 国产精品嫩草影院av在线观看| 亚洲av不卡在线观看| 亚洲人成网站在线观看播放| 免费观看在线日韩| 亚洲精品日韩av片在线观看| 亚洲精品色激情综合| 美女福利国产在线| 人妻系列 视频| 国产成人精品一,二区| 黄色日韩在线| 深夜a级毛片| 久久影院123| 亚洲av免费高清在线观看| 国产精品久久久久久久久免| 一个人免费看片子| 国产精品一区二区三区四区免费观看| 永久免费av网站大全| av免费在线看不卡| 毛片一级片免费看久久久久| 日韩一区二区三区影片| 精品亚洲成国产av| 在线播放无遮挡| 我的老师免费观看完整版| 日日撸夜夜添| 久久国产亚洲av麻豆专区| 日韩一区二区三区影片| 一级爰片在线观看| 熟女人妻精品中文字幕| 成年av动漫网址| 少妇的逼水好多| av.在线天堂| 色婷婷av一区二区三区视频| 少妇裸体淫交视频免费看高清| 日韩 亚洲 欧美在线| 日本91视频免费播放| h日本视频在线播放| 国产熟女午夜一区二区三区 | 亚洲国产精品专区欧美| av国产精品久久久久影院| 乱系列少妇在线播放| 乱码一卡2卡4卡精品| www.av在线官网国产| 国产亚洲5aaaaa淫片| 人人妻人人添人人爽欧美一区卜| 久久精品国产亚洲网站| 一级毛片aaaaaa免费看小| 九九在线视频观看精品| 久久99精品国语久久久| 91久久精品国产一区二区三区| 国产淫片久久久久久久久| 国产精品麻豆人妻色哟哟久久| 欧美最新免费一区二区三区| 深夜a级毛片| 青青草视频在线视频观看| 在线免费观看不下载黄p国产| 日韩电影二区| 亚洲欧美一区二区三区黑人 | 人人妻人人澡人人看| 国产高清有码在线观看视频| 另类精品久久| 高清视频免费观看一区二区| 这个男人来自地球电影免费观看 | 国产免费又黄又爽又色| 欧美日韩一区二区视频在线观看视频在线| 亚洲国产av新网站| 久久精品久久精品一区二区三区| 国产 精品1| 夜夜爽夜夜爽视频| 男女免费视频国产| 亚州av有码| 在现免费观看毛片| 成人综合一区亚洲| 大话2 男鬼变身卡| 国产男人的电影天堂91| 国语对白做爰xxxⅹ性视频网站| 97在线人人人人妻| 国产男女超爽视频在线观看| 国内揄拍国产精品人妻在线| 卡戴珊不雅视频在线播放| 如何舔出高潮| 99re6热这里在线精品视频| 只有这里有精品99| 久久午夜综合久久蜜桃| 一级黄片播放器| a级片在线免费高清观看视频| 亚洲精品久久久久久婷婷小说| 伦理电影免费视频| 国产亚洲一区二区精品| 在线精品无人区一区二区三| 最新的欧美精品一区二区| 一级二级三级毛片免费看| 国产精品一二三区在线看| 亚洲久久久国产精品| 国产午夜精品久久久久久一区二区三区| 久久人人爽av亚洲精品天堂| 韩国高清视频一区二区三区| 国产精品不卡视频一区二区| 免费人妻精品一区二区三区视频| 久久午夜综合久久蜜桃| 国产精品国产三级国产专区5o| 午夜福利,免费看| 成人黄色视频免费在线看| 久久韩国三级中文字幕| 九九爱精品视频在线观看| 三级经典国产精品| 国产一区二区在线观看av| 好男人视频免费观看在线| 中文字幕人妻熟人妻熟丝袜美| 插逼视频在线观看| 亚洲一级一片aⅴ在线观看| 日韩欧美精品免费久久| tube8黄色片| 午夜福利在线观看免费完整高清在| 桃花免费在线播放| 一本大道久久a久久精品| 久久狼人影院| 亚洲精品视频女| 综合色丁香网| 国产爽快片一区二区三区| 久久久久久久久久人人人人人人| 欧美bdsm另类| 男男h啪啪无遮挡| 亚洲美女搞黄在线观看| 精品国产一区二区三区久久久樱花| 久久久国产一区二区| 午夜福利在线观看免费完整高清在| 我的女老师完整版在线观看| 国产亚洲精品久久久com| h视频一区二区三区| 免费人妻精品一区二区三区视频| 校园人妻丝袜中文字幕| 黑人巨大精品欧美一区二区蜜桃 | 在线观看国产h片| 国产无遮挡羞羞视频在线观看| 不卡视频在线观看欧美| 亚洲国产av新网站| 男人爽女人下面视频在线观看| 一级毛片aaaaaa免费看小| 男人舔奶头视频| 偷拍熟女少妇极品色| 亚洲精品乱码久久久v下载方式| 男女免费视频国产| 精品久久国产蜜桃| 嘟嘟电影网在线观看| 国产男女超爽视频在线观看| 一级片'在线观看视频| 国产视频内射| 黄色一级大片看看| 亚洲国产欧美在线一区| 三级国产精品片| 亚洲中文av在线| 啦啦啦在线观看免费高清www| 国产69精品久久久久777片| 国产av码专区亚洲av| 一边亲一边摸免费视频| 99热6这里只有精品| 欧美激情极品国产一区二区三区 | 欧美成人精品欧美一级黄| 九草在线视频观看| 日韩中文字幕视频在线看片| 少妇 在线观看| 欧美激情极品国产一区二区三区 | 三上悠亚av全集在线观看 | 久久久久久伊人网av| 一本大道久久a久久精品| 欧美变态另类bdsm刘玥| 日本午夜av视频| 中文字幕精品免费在线观看视频 | 黑人高潮一二区| 人妻系列 视频| 亚洲成色77777| 免费av中文字幕在线| 亚洲国产欧美日韩在线播放 | 男男h啪啪无遮挡| 久久久久久久亚洲中文字幕| av在线观看视频网站免费| 午夜福利视频精品| 80岁老熟妇乱子伦牲交| 欧美日韩视频精品一区| 最近手机中文字幕大全| 97超碰精品成人国产| 极品人妻少妇av视频| 51国产日韩欧美| 免费大片黄手机在线观看| 亚洲av福利一区| 国产无遮挡羞羞视频在线观看| 99精国产麻豆久久婷婷| 丰满迷人的少妇在线观看| 观看免费一级毛片| 日韩av不卡免费在线播放| 久久影院123| 国产精品一区www在线观看| 中文字幕制服av| 狠狠精品人妻久久久久久综合| 久久久a久久爽久久v久久| 日韩一区二区视频免费看| 国产探花极品一区二区| 高清午夜精品一区二区三区| 国产男人的电影天堂91| 在线观看www视频免费| 欧美区成人在线视频| 美女xxoo啪啪120秒动态图| 乱码一卡2卡4卡精品| 欧美日韩在线观看h| 国产精品久久久久成人av| 国产伦理片在线播放av一区| 亚洲av男天堂| 丝袜在线中文字幕| 亚洲av日韩在线播放| 亚洲,一卡二卡三卡| 国产av一区二区精品久久| 精品久久久久久电影网| 99热这里只有是精品在线观看| 女性生殖器流出的白浆| 国产亚洲最大av| 乱人伦中国视频| 自拍偷自拍亚洲精品老妇| 热99国产精品久久久久久7| www.av在线官网国产| √禁漫天堂资源中文www| av黄色大香蕉| 黄片无遮挡物在线观看| 男人爽女人下面视频在线观看| 成人亚洲精品一区在线观看| 成人影院久久| 国产高清三级在线| 日本午夜av视频| 免费黄色在线免费观看| 人体艺术视频欧美日本| av在线播放精品| 国产午夜精品一二区理论片| 久久久久久久久大av| 国产一级毛片在线| 久久久久久久久久久久大奶| 国语对白做爰xxxⅹ性视频网站| 国产伦理片在线播放av一区| 大码成人一级视频| 精品99又大又爽又粗少妇毛片| 交换朋友夫妻互换小说| 成人无遮挡网站| 国产在线男女| 日韩中文字幕视频在线看片| 国产爽快片一区二区三区| 一区二区三区精品91| 成年人午夜在线观看视频| 成人二区视频| 精品国产一区二区三区久久久樱花| 久久人人爽av亚洲精品天堂| 日韩精品免费视频一区二区三区 | 美女内射精品一级片tv| 亚洲欧美成人精品一区二区| 久久久国产一区二区| 国产乱人偷精品视频| 3wmmmm亚洲av在线观看| 建设人人有责人人尽责人人享有的| 永久免费av网站大全| 日韩欧美精品免费久久| 亚洲中文av在线| 久久精品久久精品一区二区三区| 亚洲久久久国产精品| 亚洲精品国产成人久久av| 在线观看一区二区三区激情| av福利片在线| 免费看日本二区| 国产精品福利在线免费观看| 51国产日韩欧美| 国产成人精品无人区| 夜夜骑夜夜射夜夜干| av福利片在线观看| 91久久精品国产一区二区成人| 亚洲无线观看免费| 一级毛片电影观看| 青春草视频在线免费观看| 爱豆传媒免费全集在线观看| 久久精品国产鲁丝片午夜精品| av播播在线观看一区| 亚洲精品第二区| 久久久久久久久久久久大奶| 成年人午夜在线观看视频| 一级毛片黄色毛片免费观看视频| 亚洲婷婷狠狠爱综合网| 久久久久国产网址| 人人妻人人添人人爽欧美一区卜| 国产精品人妻久久久久久| 丰满乱子伦码专区| 80岁老熟妇乱子伦牲交| 又粗又硬又长又爽又黄的视频| av黄色大香蕉| 国产成人91sexporn| 亚洲不卡免费看| 国产伦精品一区二区三区视频9| 亚洲精品乱码久久久v下载方式| 视频区图区小说| 狠狠精品人妻久久久久久综合| 日韩中字成人| 亚洲精品中文字幕在线视频 | 午夜久久久在线观看| 中文字幕久久专区| 国产国拍精品亚洲av在线观看| 亚洲三级黄色毛片| 少妇的逼好多水| 久久这里有精品视频免费| 2018国产大陆天天弄谢| 久久国产精品男人的天堂亚洲 | 日本91视频免费播放| 国产精品国产三级国产av玫瑰| 精品亚洲成a人片在线观看| 99久久中文字幕三级久久日本| 丝袜喷水一区| 一本色道久久久久久精品综合| 国产亚洲av片在线观看秒播厂| 成人国产av品久久久| 日日爽夜夜爽网站| 岛国毛片在线播放| 在线观看免费高清a一片| 波野结衣二区三区在线| 插阴视频在线观看视频| 国产av码专区亚洲av| 精品久久久精品久久久| 你懂的网址亚洲精品在线观看| 国产精品久久久久成人av| 男女免费视频国产| 下体分泌物呈黄色| 亚洲激情五月婷婷啪啪| 亚洲精品成人av观看孕妇| 日本与韩国留学比较| 免费看日本二区| 免费播放大片免费观看视频在线观看| 国产欧美日韩一区二区三区在线 | 男女啪啪激烈高潮av片| 久久精品夜色国产| av播播在线观看一区|