• <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例
    黄色欧美视频在线观看| 99久久人妻综合| 亚洲第一av免费看| 色5月婷婷丁香| 亚洲av二区三区四区| 免费看av在线观看网站| 大码成人一级视频| 国产成人精品久久久久久| 丁香六月天网| 精品亚洲乱码少妇综合久久| 精品人妻熟女av久视频| 久久韩国三级中文字幕| 免费观看无遮挡的男女| 国产国语露脸激情在线看| 亚洲熟女精品中文字幕| 精品国产一区二区三区久久久樱花| 国产精品国产av在线观看| 亚洲久久久国产精品| 欧美亚洲 丝袜 人妻 在线| 亚洲内射少妇av| 五月天丁香电影| videos熟女内射| 国产精品一区二区在线观看99| 精品国产乱码久久久久久小说| 欧美日韩视频高清一区二区三区二| 最近的中文字幕免费完整| 久久国产精品男人的天堂亚洲 | 成人手机av| 国产亚洲午夜精品一区二区久久| 在线免费观看不下载黄p国产| 97在线人人人人妻| 国产免费福利视频在线观看| 日韩大片免费观看网站| 美女国产高潮福利片在线看| 国产精品人妻久久久久久| 日韩av在线免费看完整版不卡| 久久人人爽人人爽人人片va| 狂野欧美白嫩少妇大欣赏| 久久这里有精品视频免费| 高清黄色对白视频在线免费看| av又黄又爽大尺度在线免费看| 高清不卡的av网站| 久热久热在线精品观看| 亚洲成色77777| 日本wwww免费看| 热re99久久国产66热| 女的被弄到高潮叫床怎么办| 欧美激情国产日韩精品一区| 亚洲欧美成人综合另类久久久| 成人18禁高潮啪啪吃奶动态图 | 成人午夜精彩视频在线观看| 你懂的网址亚洲精品在线观看| 视频区图区小说| 免费大片18禁| 九色亚洲精品在线播放| 天天影视国产精品| 久久精品人人爽人人爽视色| 哪个播放器可以免费观看大片| 国产精品人妻久久久影院| 成人亚洲欧美一区二区av| 99热这里只有是精品在线观看| 国产视频首页在线观看| 一本色道久久久久久精品综合| 黄色怎么调成土黄色| 我要看黄色一级片免费的| 熟妇人妻不卡中文字幕| 日韩在线高清观看一区二区三区| 夫妻性生交免费视频一级片| 久久精品人人爽人人爽视色| 街头女战士在线观看网站| xxx大片免费视频| 一本大道久久a久久精品| 在线观看人妻少妇| 亚洲精品视频女| 欧美另类一区| 插阴视频在线观看视频| 成人无遮挡网站| 99热全是精品| 久久久久久久久久久丰满| 亚洲久久久国产精品| 51国产日韩欧美| 美女视频免费永久观看网站| 啦啦啦啦在线视频资源| 亚洲五月色婷婷综合| 亚洲欧美清纯卡通| av在线播放精品| 久久青草综合色| 亚洲精品色激情综合| 只有这里有精品99| 色婷婷av一区二区三区视频| 国产免费福利视频在线观看| 精品国产露脸久久av麻豆| 久久久久久久久久成人| 在线观看一区二区三区激情| 99re6热这里在线精品视频| 18在线观看网站| 一级毛片电影观看| 2021少妇久久久久久久久久久| 丰满饥渴人妻一区二区三| 国产成人一区二区在线| 精品国产乱码久久久久久小说| 大片免费播放器 马上看| 国产乱人偷精品视频| 一级毛片黄色毛片免费观看视频| 国产亚洲午夜精品一区二区久久| 飞空精品影院首页| 国产免费福利视频在线观看| 男女高潮啪啪啪动态图| 在线亚洲精品国产二区图片欧美 | 热re99久久国产66热| 在线观看人妻少妇| 两个人的视频大全免费| 中文精品一卡2卡3卡4更新| 一本久久精品| av女优亚洲男人天堂| 国产精品欧美亚洲77777| 日日撸夜夜添| 国产日韩一区二区三区精品不卡 | 亚洲欧美中文字幕日韩二区| 精品久久久噜噜| 成年人免费黄色播放视频| 日日摸夜夜添夜夜添av毛片| 最近中文字幕2019免费版| 久久人妻熟女aⅴ| 男女边摸边吃奶| 最近最新中文字幕免费大全7| 在线观看三级黄色| 免费播放大片免费观看视频在线观看| 日韩av在线免费看完整版不卡| 免费高清在线观看日韩| 亚洲人成网站在线观看播放| 一区二区日韩欧美中文字幕 | 亚洲精品成人av观看孕妇| 亚洲国产精品成人久久小说| 日韩免费高清中文字幕av| 久久人人爽人人片av| 天堂俺去俺来也www色官网| 亚洲精品第二区| 国产日韩欧美亚洲二区| 在线精品无人区一区二区三| 极品人妻少妇av视频| 在线 av 中文字幕| 国产伦理片在线播放av一区| 爱豆传媒免费全集在线观看| 久久精品久久久久久久性| videosex国产| 精品久久国产蜜桃| 亚洲人成网站在线播| 日本-黄色视频高清免费观看| 丰满乱子伦码专区| 久久久久久久久久久丰满| 我的老师免费观看完整版| 欧美性感艳星| 久久久久久久国产电影| 制服丝袜香蕉在线| 久久99热这里只频精品6学生| 亚洲不卡免费看| 亚洲欧美精品自产自拍| 久久精品久久精品一区二区三区| 男女无遮挡免费网站观看| 在线观看免费日韩欧美大片 | 在线看a的网站| 91久久精品国产一区二区三区| 嘟嘟电影网在线观看| 熟女电影av网| 亚洲精品久久午夜乱码| 在线观看免费视频网站a站| 亚洲美女黄色视频免费看| av卡一久久| 一级毛片aaaaaa免费看小| 日韩av不卡免费在线播放| 国产视频首页在线观看| 在线观看人妻少妇| 哪个播放器可以免费观看大片| 亚洲综合精品二区| 亚洲在久久综合| 一本—道久久a久久精品蜜桃钙片| 嫩草影院入口| 欧美 日韩 精品 国产| 欧美三级亚洲精品| 中文精品一卡2卡3卡4更新| 色婷婷久久久亚洲欧美| 久久精品人人爽人人爽视色| 成人漫画全彩无遮挡| 精品99又大又爽又粗少妇毛片| 亚洲国产精品一区二区三区在线| 91久久精品国产一区二区三区| 久久久久国产精品人妻一区二区| 中文字幕免费在线视频6| 中文字幕亚洲精品专区| 亚洲国产毛片av蜜桃av| 黑丝袜美女国产一区| 国产精品成人在线| av一本久久久久| 亚洲欧洲精品一区二区精品久久久 | 大话2 男鬼变身卡| 一本久久精品| av女优亚洲男人天堂| 午夜免费观看性视频| 日韩av在线免费看完整版不卡| 成人国产麻豆网| 成人漫画全彩无遮挡| 97精品久久久久久久久久精品| 中文字幕制服av| 国产精品不卡视频一区二区| 久热这里只有精品99| 黄色欧美视频在线观看| 亚洲精品久久成人aⅴ小说 | 永久免费av网站大全| 久久久久久久久久久久大奶| 18禁裸乳无遮挡动漫免费视频| 国产精品99久久99久久久不卡 | 97在线视频观看| 人人澡人人妻人| 国产精品国产三级专区第一集| 国产精品国产三级国产av玫瑰| 最近手机中文字幕大全| 日本91视频免费播放| 久久久久久久久久久久大奶| 国产亚洲最大av| 亚洲欧美一区二区三区国产| 日韩一区二区视频免费看| 国产精品国产三级专区第一集| 欧美日韩视频高清一区二区三区二| 麻豆精品久久久久久蜜桃| kizo精华| 91午夜精品亚洲一区二区三区| 丰满迷人的少妇在线观看| 亚洲国产欧美日韩在线播放| 成年人免费黄色播放视频| 最黄视频免费看| 在线观看美女被高潮喷水网站| 在线免费观看不下载黄p国产| 精品亚洲乱码少妇综合久久| 免费av不卡在线播放| 日本色播在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 一本一本综合久久| 国产免费又黄又爽又色| 中文字幕人妻丝袜制服| 色婷婷av一区二区三区视频| 精品一区在线观看国产| 久久97久久精品| 亚洲精品中文字幕在线视频| 日日摸夜夜添夜夜添av毛片| 伦理电影大哥的女人| 免费高清在线观看视频在线观看| 一级二级三级毛片免费看| 黄色视频在线播放观看不卡| 亚洲第一av免费看| 99视频精品全部免费 在线| 精品国产一区二区久久| 久久精品国产亚洲网站| 欧美xxⅹ黑人| 精品国产国语对白av| 久久这里有精品视频免费| 视频区图区小说| 狂野欧美激情性bbbbbb| 久久久久国产网址| 国产日韩一区二区三区精品不卡 | 免费观看在线日韩| 天堂8中文在线网| 亚洲三级黄色毛片| 97超视频在线观看视频| 午夜免费鲁丝| 热99久久久久精品小说推荐| 中文精品一卡2卡3卡4更新| 欧美成人午夜免费资源| 亚洲熟女精品中文字幕| 国产视频内射| 久久久午夜欧美精品| 精品久久久精品久久久| 色94色欧美一区二区| 国产综合精华液| 黄色欧美视频在线观看| 一本—道久久a久久精品蜜桃钙片| 日韩一区二区三区影片| 寂寞人妻少妇视频99o| 卡戴珊不雅视频在线播放| 亚洲欧美成人精品一区二区| 免费观看av网站的网址| 国产亚洲最大av| 欧美精品高潮呻吟av久久| 国产极品天堂在线| 久久 成人 亚洲| 欧美成人午夜免费资源| 少妇高潮的动态图| 久久久亚洲精品成人影院| 国产成人精品婷婷| 国产无遮挡羞羞视频在线观看| 国产在线视频一区二区| 人人妻人人澡人人看| 亚洲av免费高清在线观看| 亚洲成色77777| 成人漫画全彩无遮挡| 成人免费观看视频高清| 2018国产大陆天天弄谢| 又黄又爽又刺激的免费视频.| 视频区图区小说| 免费看av在线观看网站| 午夜福利网站1000一区二区三区| av一本久久久久| 妹子高潮喷水视频| 亚洲精品一二三| 午夜日本视频在线| 亚洲精品国产色婷婷电影| 国产精品国产三级国产专区5o| 黄色配什么色好看| 人妻夜夜爽99麻豆av| 大片电影免费在线观看免费| 最后的刺客免费高清国语| 老熟女久久久| 满18在线观看网站| 国产精品三级大全| 人妻系列 视频| 亚洲,欧美,日韩| 99热网站在线观看| 国语对白做爰xxxⅹ性视频网站| 自线自在国产av| 欧美人与善性xxx| 人人妻人人澡人人看| 欧美日韩视频高清一区二区三区二| 蜜桃久久精品国产亚洲av| 亚洲精品美女久久av网站| 日产精品乱码卡一卡2卡三| 在线观看美女被高潮喷水网站| 纯流量卡能插随身wifi吗| 人妻夜夜爽99麻豆av| 搡老乐熟女国产| 免费看av在线观看网站| 国产欧美亚洲国产| 久久国产精品大桥未久av| 免费人妻精品一区二区三区视频| 亚洲综合色惰| 丝袜美足系列| 国产精品久久久久久久久免| 岛国毛片在线播放| 爱豆传媒免费全集在线观看| 中文字幕亚洲精品专区| 99热国产这里只有精品6| 美女视频免费永久观看网站| 一区在线观看完整版| 妹子高潮喷水视频| 男人操女人黄网站| 嘟嘟电影网在线观看| 国产黄频视频在线观看| 建设人人有责人人尽责人人享有的| 日韩av在线免费看完整版不卡| av女优亚洲男人天堂| 亚洲欧美一区二区三区黑人 | 中文精品一卡2卡3卡4更新| 老女人水多毛片| 美女国产视频在线观看| 乱人伦中国视频| 亚洲欧美成人综合另类久久久| 女的被弄到高潮叫床怎么办| 一区在线观看完整版| 婷婷色av中文字幕| 国产精品一国产av| 国产精品久久久久久久久免| 99热国产这里只有精品6| 午夜精品国产一区二区电影| 国产极品粉嫩免费观看在线 | 日本黄大片高清| 欧美 亚洲 国产 日韩一| 一级毛片黄色毛片免费观看视频| 亚洲四区av| 精品国产一区二区三区久久久樱花| 欧美日韩国产mv在线观看视频| 成人亚洲精品一区在线观看| 亚洲国产最新在线播放| 99国产综合亚洲精品| 人妻人人澡人人爽人人| 国产熟女欧美一区二区| 在现免费观看毛片| 天天躁夜夜躁狠狠久久av| av女优亚洲男人天堂| 亚洲精品aⅴ在线观看| 国产老妇伦熟女老妇高清| 满18在线观看网站| 丝袜脚勾引网站| 母亲3免费完整高清在线观看 | 国产欧美另类精品又又久久亚洲欧美| 欧美成人午夜免费资源| 色婷婷av一区二区三区视频| 国产一区有黄有色的免费视频| 爱豆传媒免费全集在线观看| 夜夜爽夜夜爽视频| 国产亚洲午夜精品一区二区久久| 另类亚洲欧美激情| 国产日韩一区二区三区精品不卡 | 十分钟在线观看高清视频www| videos熟女内射| 成人国产av品久久久| 精品国产露脸久久av麻豆| 最后的刺客免费高清国语| 精品国产露脸久久av麻豆| 最近手机中文字幕大全| 午夜福利网站1000一区二区三区| 一级黄片播放器| 欧美亚洲 丝袜 人妻 在线| 久久久精品区二区三区| 久久精品久久久久久噜噜老黄| 精品99又大又爽又粗少妇毛片| 搡女人真爽免费视频火全软件| av线在线观看网站| 国产精品国产三级国产av玫瑰| 黄色视频在线播放观看不卡| 一级毛片黄色毛片免费观看视频| 国产老妇伦熟女老妇高清| 天天躁夜夜躁狠狠久久av| 亚洲五月色婷婷综合| 只有这里有精品99| 在线亚洲精品国产二区图片欧美 | 亚洲第一av免费看| 欧美国产精品一级二级三级| 欧美+日韩+精品| 久久97久久精品| 亚洲经典国产精华液单| 免费日韩欧美在线观看| 久久热精品热| 在线观看www视频免费| 观看av在线不卡| 亚洲精品乱码久久久v下载方式| 国产成人免费观看mmmm| 日本欧美视频一区| 国产视频首页在线观看| 欧美激情国产日韩精品一区| 高清不卡的av网站| 亚洲,欧美,日韩| 国国产精品蜜臀av免费| 精品亚洲成a人片在线观看| 免费播放大片免费观看视频在线观看| 三级国产精品片| 久久久久久久久久久久大奶| 美女cb高潮喷水在线观看| .国产精品久久| 中文字幕最新亚洲高清| 精品人妻一区二区三区麻豆| 搡女人真爽免费视频火全软件| 免费看av在线观看网站| 久久免费观看电影| 亚洲欧美清纯卡通| 国产精品.久久久| 亚洲精品国产色婷婷电影| 亚洲av成人精品一区久久| 国产精品一国产av| 少妇人妻精品综合一区二区| 免费高清在线观看日韩| 久久综合国产亚洲精品| 精品人妻在线不人妻| 免费播放大片免费观看视频在线观看| 91成人精品电影| 高清av免费在线| 日韩欧美一区视频在线观看| 一级二级三级毛片免费看| 午夜福利,免费看| 午夜福利视频精品| 日本黄大片高清| 国产黄色视频一区二区在线观看| 午夜av观看不卡| 久久久久久久国产电影| 亚洲精品色激情综合| 国产精品国产av在线观看| 久久久久久久精品精品| 国产有黄有色有爽视频| 亚洲精品第二区| 免费人妻精品一区二区三区视频| 亚洲人与动物交配视频| 欧美激情 高清一区二区三区| 在线精品无人区一区二区三| 97超视频在线观看视频| 亚洲美女搞黄在线观看| 欧美精品一区二区免费开放| 久久久久国产精品人妻一区二区| 永久免费av网站大全| 婷婷色综合www| 国产有黄有色有爽视频| 精品少妇久久久久久888优播| 精品亚洲乱码少妇综合久久| 国产av精品麻豆| 日日摸夜夜添夜夜爱| 久热这里只有精品99| 一级毛片 在线播放| 三上悠亚av全集在线观看| 18在线观看网站| a级毛色黄片| 男人操女人黄网站| 女的被弄到高潮叫床怎么办| 久久99精品国语久久久| 九色亚洲精品在线播放| 亚洲一区二区三区欧美精品| 简卡轻食公司| 午夜福利网站1000一区二区三区| 日韩免费高清中文字幕av| 三级国产精品欧美在线观看| 亚洲人与动物交配视频| 亚洲激情五月婷婷啪啪| 欧美日韩av久久| 日本黄大片高清| 欧美日韩av久久| 男的添女的下面高潮视频| 国产精品99久久99久久久不卡 | 婷婷色麻豆天堂久久| 在线免费观看不下载黄p国产| 天天躁夜夜躁狠狠久久av| 色婷婷久久久亚洲欧美| 免费看光身美女| 嫩草影院入口| 满18在线观看网站| 国产成人aa在线观看| 欧美精品人与动牲交sv欧美| 曰老女人黄片| 国产永久视频网站| 在线免费观看不下载黄p国产| 高清黄色对白视频在线免费看| 久久国内精品自在自线图片| 人妻制服诱惑在线中文字幕| 精品国产国语对白av| 男女免费视频国产| 在现免费观看毛片| 满18在线观看网站| 国产亚洲精品久久久com| 夜夜看夜夜爽夜夜摸| 精品国产国语对白av| 高清视频免费观看一区二区| 免费观看无遮挡的男女| 欧美日韩一区二区视频在线观看视频在线| 老司机亚洲免费影院| 99精国产麻豆久久婷婷| 最近手机中文字幕大全| 高清午夜精品一区二区三区| 日韩亚洲欧美综合| 人妻制服诱惑在线中文字幕| 日本爱情动作片www.在线观看| 极品少妇高潮喷水抽搐| 国产国拍精品亚洲av在线观看| 少妇被粗大的猛进出69影院 | 亚洲精品,欧美精品| 亚洲成人一二三区av| 国产熟女午夜一区二区三区 | 在线观看免费日韩欧美大片 | 亚洲精品久久久久久婷婷小说| 亚洲精品成人av观看孕妇| 一区在线观看完整版| 老女人水多毛片| 2018国产大陆天天弄谢| 一级毛片aaaaaa免费看小| a级毛片免费高清观看在线播放| 精品视频人人做人人爽| 国产亚洲午夜精品一区二区久久| 国产精品人妻久久久久久| 欧美另类一区| 尾随美女入室| 国产黄片视频在线免费观看| 精品少妇黑人巨大在线播放| 一区二区三区四区激情视频| 日本免费在线观看一区| 精品卡一卡二卡四卡免费| 一级毛片黄色毛片免费观看视频| a级毛片黄视频| 亚洲av中文av极速乱| 在线免费观看不下载黄p国产| 免费人妻精品一区二区三区视频| 亚洲精品中文字幕在线视频| 亚洲怡红院男人天堂| 日日啪夜夜爽| 亚洲伊人久久精品综合| 国产一区有黄有色的免费视频| 高清av免费在线| 一级毛片aaaaaa免费看小| 国国产精品蜜臀av免费| 日本黄色日本黄色录像| av国产久精品久网站免费入址| 亚洲国产精品一区二区三区在线| 一本一本综合久久| 亚洲精品亚洲一区二区| 三级国产精品片| 婷婷色综合大香蕉| 九九爱精品视频在线观看| 国产精品国产三级国产av玫瑰| 秋霞伦理黄片| 黑人欧美特级aaaaaa片| 日韩熟女老妇一区二区性免费视频| 高清不卡的av网站| 人人妻人人澡人人爽人人夜夜| 在现免费观看毛片| freevideosex欧美| 免费黄网站久久成人精品| 精品久久蜜臀av无| 亚洲国产精品999| 丝袜美足系列| 国产永久视频网站| 如何舔出高潮| 久久久久久伊人网av| 日韩欧美一区视频在线观看| 日本av手机在线免费观看| 免费观看av网站的网址| 亚洲五月色婷婷综合| av免费观看日本| 丁香六月天网| 国产亚洲最大av| 国产免费现黄频在线看| 久久久久网色| 国产黄频视频在线观看| 亚洲国产精品国产精品| 九色成人免费人妻av| 欧美日本中文国产一区发布| 久久国内精品自在自线图片| 亚洲精品视频女| 国产爽快片一区二区三区| 亚洲av福利一区| 在线精品无人区一区二区三| 成人无遮挡网站| 亚洲精品自拍成人|