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

    基于Excel的直剪試驗數(shù)據(jù)處理方法

    2015-03-13 08:44:54莫林利趙秀紹
    鐵道建筑 2015年9期
    關(guān)鍵詞:黏聚力土樣抗剪

    莫林利,趙秀紹,王 旭,鄭 偉

    (1.華東交通大學(xué) 軟件學(xué)院,江西 南昌 330013;2.華東交通大學(xué) 土木建筑學(xué)院,江西 南昌 330013)

    基于Excel的直剪試驗數(shù)據(jù)處理方法

    莫林利1,趙秀紹2,王 旭2,鄭 偉2

    (1.華東交通大學(xué) 軟件學(xué)院,江西 南昌 330013;2.華東交通大學(xué) 土木建筑學(xué)院,江西 南昌 330013)

    回歸分析中的最小二乘法是處理直剪試驗數(shù)據(jù)誤差最小的方法,但計算工作量大且容易出錯。應(yīng)用Excel軟件給出可用于直剪試驗計算的5種方法,可顯著減少計算工作量并降低出錯率。實際工程數(shù)據(jù)計算對比表明,除基于點線距的規(guī)劃求解法外,其余4種方法均為基于擬合誤差的計算方法,兩種原理的計算結(jié)果差距可忽略不計。對于無黏性土情況,除INTERCEPT函數(shù)與SLOPE函數(shù)聯(lián)合求解法無法實現(xiàn)強制黏聚力為0外,其余方法均可實現(xiàn)。實際應(yīng)用表明Excel圖解法、TREND函數(shù)法及截距和斜率函數(shù)聯(lián)合求解法操作簡單,易于應(yīng)用。

    直剪試驗 Excel內(nèi)部函數(shù) 最小二乘法 規(guī)劃求解

    邊坡土的抗剪強度是鐵路路基邊坡穩(wěn)定性定量評價和計算滑坡推力的重要參數(shù),有時內(nèi)摩擦角相差1°~2°,滑坡推力就可能成倍增加[1]。馬顯春等[2]采用正交分析法對滑坡穩(wěn)定性的影響因子進(jìn)行了分析,得出滑坡影響因素敏感性順序為內(nèi)摩擦角>黏聚力>地下水位>重度>坡高??梢姾侠泶_定巖土材料的抗剪強度指標(biāo)對于工程設(shè)計、施工等十分重要[3-5]。

    直剪試驗的數(shù)據(jù)處理實際上是求出土的抗剪強度參數(shù),即黏聚力 c和內(nèi)摩擦角 φ。斜截式計算法是在試驗中測得N,M兩點的垂直壓應(yīng)力σN,σM,其對應(yīng)的抗剪強度為τN,τM。通過 N,M的直線方程可以獲得土的抗剪強度指標(biāo)[6]。此法的不足之處是非常明顯的,計算時僅能使用2組數(shù)據(jù)?!锻凉ぴ囼灱夹g(shù)手冊》[7]建議的方法為繪制抗剪強度與垂直壓力關(guān)系曲線,直線的傾角為內(nèi)摩擦角φ,縱坐標(biāo)的截距為黏聚力c。但是在繪庫侖破壞直線時,一般是繪圖者根據(jù)經(jīng)驗和直觀判斷,盡量讓落在直線兩邊的點數(shù)大致相等,然后求出截距和斜率。由于它是人為地根據(jù)經(jīng)驗和直觀判斷得出的,c,φ值往往帶有一定的人為性和不確定性。采用最小二乘法對試驗數(shù)據(jù)進(jìn)行擬合求解[8-11],計算量大且容易出錯。本文探討如何巧用Excel的各項功能簡化計算過程,并考慮直剪試驗中c值為負(fù)或強制c=0(無黏性土)的情況,使數(shù)據(jù)處理達(dá)到全自動化。

    1 最小二乘法求解原理

    為了便于說明各種求解方法的原理與實現(xiàn)方法,以江西撫州八景鎮(zhèn)火車站路基工程一組土樣實測數(shù)據(jù)(見表1)為例進(jìn)行分析。

    表1 撫州八景鎮(zhèn)路基土樣測試數(shù)據(jù)

    直剪試驗求解中,當(dāng)繪制測點的破壞包線時,由于土樣不均勻,試驗讀數(shù)誤差等原因,往往繪制的直線不能通過所有的測點,總是存在一些誤差,根據(jù)這些誤差的平方和為最小,可以確定一條近似的破壞包線。設(shè)直剪試驗中每組n個土樣的垂直應(yīng)力、抗剪強度分別為σi,τi,利用庫侖公式 τ=σtanφ+c進(jìn)行擬合,其中c,φ被視為待定系數(shù)。以偏差+c)]2最小為條件來確定c,φ值。運用高等數(shù)學(xué)的多元函數(shù)求極值的方法即可得到當(dāng) M為最小時的 c,φ值。分別對c,φ求一階偏導(dǎo)并令其等于0得[8]

    根據(jù)式(1),式(2)解得

    如果用上述數(shù)學(xué)法求解,計算工作量大且易出錯。Excel軟件中有豐富的回歸分析函數(shù)和規(guī)劃求解方法,只要挖掘其功能,就可以簡單有效地完成上述計算。

    2 直剪試驗的Excel求解方法

    2.1 Excel圖解法

    Excel軟件有強大的繪圖工具,先把試驗數(shù)據(jù)填入數(shù)據(jù)表格中,然后插入圖表,并為數(shù)據(jù)添加趨勢線和顯示公式,就可以得到直剪試驗圖。具體實現(xiàn)方法為:選擇正應(yīng)力和剪應(yīng)力數(shù)據(jù),選擇“插入”→“圖表”→“散點圖”,然后選擇圖形中的點后“添加趨勢線”,在“選項”欄中選擇“顯示公式”,可直接得到擬合公式 y= 0.430 6x+3.260 9,進(jìn)而得 c=3.260 9 kPa,tanφ= 0.430 6,即φ=23.297 1°。

    2.2 TREND函數(shù)求解法

    2.2.1 TREND函數(shù)求解的基本原理

    由于直剪試驗各測點存在線性關(guān)系,這樣就可以方便地利用 Excel中的線性插值函數(shù) TREND。TREND函數(shù)用途:返回一條線性回歸擬合線的一組縱坐標(biāo)值(y值)。即找到適合給定的數(shù)組known_y's和known_x's的直線,并返回指定數(shù)組new_x's值在直線上對應(yīng)的y值。TREND函數(shù)的語法為TREND(known _y's,known_x's,new_x's,const)。其中known_y's為已知的τ=σtanφ+c=σk+c中的τi值集合,known_x's為已知關(guān)系τ=σk+c中可選的σi值的集合,new_x's為需要函數(shù)TREND返回對應(yīng)新σ值的新τ值,const參數(shù)為邏輯值指明是否強制常數(shù)項c為0。當(dāng)擬合發(fā)現(xiàn)c值為負(fù)時,由于黏聚力不可能為負(fù)值,因此可以強制c值為0。

    TREND函數(shù)求解原理如圖1所示,很顯然,對于抗剪強度包線τ=σtanφ+c,當(dāng)new_x's=σ=0時,其返回的函數(shù)值為τ0=c,即為黏聚力值。當(dāng)new_x's= σ=1時,函數(shù)返回值 τ1=1×tanφ+c,根據(jù)圖1中的幾何關(guān)系可得tanφ=(τ1-τ0)/1=τ1-τ0。

    2.2.2 TREND函數(shù)求解法的實現(xiàn)

    將表1中的原始數(shù)據(jù)填入Excel表中,如圖2(屏幕截圖)中所示。

    圖1 TREND函數(shù)求解原理

    圖2 TREND函數(shù)求解抗剪強度指標(biāo)示例

    在C2和D2單元格中,分別輸入“=TREND(B2: B5,A2:A5,0)”,“=ATAN(TREND(B2:B5,A2:A5,1)-TREND(B2:B5,A2:A5,0))*180/PI()”,即可求得該土樣的黏聚力c和內(nèi)摩擦角φ。

    TREND(B2:B5,A2:A5,0)在σ=0時可得τ0=c; TREND(B2:B5,A2:A5,1)在 σ=1 kPa時可得 τ1=c +k;TREND(B2:B5,A2:A5,1)-TREND(B2:B5,A2: A5,0)為系數(shù)k。

    當(dāng)土樣為無黏性土或計算出的黏聚力c值為負(fù)值時,可以強制TREND()函數(shù)中const項參數(shù)為0。對表1的測試數(shù)據(jù)求解結(jié)果為:c=3.260 9 kPa,φ= 23.297 1°。

    2.3 INTERCEPT函數(shù)與SLOPE函數(shù)聯(lián)合求解法

    2.3.1 INTERCEPT函數(shù)與SLOPE函數(shù)聯(lián)合求解原理

    在Excel軟件中,函數(shù)INTERCEPT和SLOPE利用現(xiàn)有的x值和y值計算直線與y軸的截距和斜率。其中INTERCEPT(known_y's,known_x's)函數(shù)可以利用現(xiàn)有的σ值和τ值計算直線與τ軸的截距,即抗剪強度指標(biāo)黏聚力 c;SLOPE(known_y's,known_x's)函數(shù)可以利用現(xiàn)有的σ值和τ值計算直線的斜率 k,即抗剪包線斜率tanφ。known_y's表示因變的觀察值或數(shù)據(jù)集合;known_x's表示自變的觀察值或者數(shù)據(jù)集合。

    2.3.2 INTERCEPT函數(shù)與 SLOPE函數(shù)聯(lián)合求解法的實現(xiàn)

    設(shè)數(shù)據(jù)在 Excel表中的排列仍如圖2所示,則在C2和D2單元格中分別輸入“=INTERCEPT(B2:B5,A2:A5)”,“=ATAN(SLOPE(B2:B5,A2:A5))* 180/PI()”。求解出的結(jié)果為 c=3.260 9 kPa,φ= 23.297 1°。

    2.4 基于點線距的規(guī)劃求解法

    2.4.1 基本原理

    如圖3所示,對抗剪強度庫侖公式τ=σtanφ+c,令tanφ=k,則任意測點 A(σi,τi)到直線的距離為di[12],見式(5)。

    圖3 點線距求解基本原理示意

    根據(jù)最小二乘法的原理,回歸分析時就是令 d2i和達(dá)到最小,即

    式(6)中有k和c兩個變量,實際為線性規(guī)劃求解問題,求解目標(biāo)函數(shù)的最小值。Excel規(guī)劃求解一般步驟為建立目標(biāo)函數(shù)f(x)min,確定約束條件將決策變量用數(shù)學(xué)公式表達(dá)出來,使需要優(yōu)化的目標(biāo)及評判標(biāo)準(zhǔn)量化為這組變量的函數(shù)。

    2.4.2 直剪試驗的規(guī)劃求解實現(xiàn)

    Excel規(guī)劃求解首先要設(shè)計好一個工作表(如圖4所示),根據(jù)建立的數(shù)學(xué)模型,將決策變量、約束條件、目標(biāo)函數(shù)等依次排列,并在工作表中將相關(guān)固定信息初始化,具體實現(xiàn)步驟如下。

    圖4 基于點線距的規(guī)劃求解示例

    1)打開Microsoft Excel軟件后,按圖4基于點線距的規(guī)劃求解所示的格式建立工作表,將有關(guān)數(shù)據(jù)輸入工作表中(其中C2和D2為空白,有待計算),E2單元格輸入“=((C2*A2-B2+D2)/SQRT(C2^2+1)) ^2+((C2*A3-B3+D2)/SQRT(C2^2+1))^2+ ((C2*A4-B4+D2)/SQRT(C2^2+1))^2+((C2* A5-B5+D2)/SQRT(C2^2+1))^2”。

    2)在“工具”欄下拉菜單中找到“規(guī)劃求解”,如果沒有此選項,可以選擇“加載宏”將“規(guī)劃求解”加載進(jìn)來。單擊“規(guī)劃求解”,出現(xiàn)“規(guī)劃求解參數(shù)”對話框,設(shè)置目標(biāo)單元格:$E$2,在“等于”選項中選“最小值”,在“可變單元格”中輸入$F$2:$G$2,在“約束”中單擊“添加”后,逐個輸入下式:“$C$2>=0”,“$D$2>=0”,最后按“確定”。當(dāng)需要強制 c值為0時,則條件變?yōu)椤?C$2>=0”,“$D$2=0”。求解結(jié)果為c=3.254 6 kPa,k=0.430 6,在C3單元格中求得φ=23.298 8°,此時目標(biāo)函數(shù)最小,f(x)min=5.618 9。

    2.5 基于測試值與擬合值誤差的規(guī)劃求解

    直剪試驗數(shù)據(jù)σi,τi利用庫侖公式擬合直剪試驗中每組n個土樣的試驗數(shù)據(jù)時,c,φ被視為待定系數(shù)。以偏差最小為條件從而確定c,k值。這仍為線性規(guī)劃求解問題,求解目標(biāo)函數(shù)的最小值。其實現(xiàn)方法同2.4中方法類似,僅需要根據(jù)公式更改其目標(biāo)函數(shù)為“=(B2-A2*C2-D2)^2+(B3-A3*C2-D2)^2+ (B4-A4*C2-D2)^2+(B5-A5*C2-D2)^2”。

    運行規(guī)劃求解后,選擇“保存極限值”即可同時求得c值和k值,然后根據(jù)k值可計算出內(nèi)摩擦角φ值。根據(jù)上述方法,表1數(shù)據(jù)的求解結(jié)果為c=3.261 0 kPa,φ=23.297 1°,此時目標(biāo)函數(shù)最小,M=6.660 9。

    3 各種直剪試驗求解方法的對比

    為了對比5種求解方法的有效性,以表1中測試數(shù)據(jù)作為第1組數(shù)據(jù),以福建省福清市江陰鎮(zhèn)海邊固結(jié)淤泥測試數(shù)據(jù)作為第2組數(shù)據(jù),第2組數(shù)據(jù)如表2所示。

    表2 福清江陰鎮(zhèn)淤泥土樣測試數(shù)據(jù)

    通過在Excel中求解,5種方法求解結(jié)果如表3所示。根據(jù)表3對2組試驗的計算結(jié)果,分析如下:

    1)常規(guī)數(shù)據(jù)計算結(jié)果對比。雖然各種方法均采用了最小二乘法,在不強制黏聚力為0的情況下,除基于點線距的Excel規(guī)劃求解法外,其余方法求解結(jié)果均相同,說明其余4種方法均是基于擬合誤差的求解方法。從數(shù)據(jù)對比可知5種方法計算結(jié)果非常接近,

    均可作為直剪試驗的求解結(jié)果。

    表3 Excel中5種直剪試驗求解方法計算結(jié)果對比

    2)黏聚力為負(fù)值情況處理對比。表3中福州省福清市江陰鎮(zhèn)由于取樣撓動等原因,多組測試數(shù)據(jù)測試結(jié)果擬合后,黏聚力均為負(fù)值,這與理論并不相符。因此對于c值為負(fù)的情況或無黏性土的黏聚力為0的情況,此時需要強制黏聚力為0。對于強制c=0的情況,INTERCEPT與SLOPE函數(shù)聯(lián)合求解法無法完成強制c=0,其余方法均可強制c=0。通過對比可得,對于無黏性土或黏聚力出現(xiàn)負(fù)值情況,不宜采用INTERCEPT與SLOPE函數(shù)聯(lián)合求解法。

    3)操作簡易程度對比。從5種方法在數(shù)據(jù)更換后Excel操作繁簡比較,規(guī)劃求解需要重新求解,特別是試驗組數(shù)更改后,目標(biāo)函數(shù)應(yīng)隨著更改,而且會生成中間報告,故對計算機(jī)操作不太熟悉人員有些難度。其它方法均操作簡便,且計算結(jié)果會自動完成更新。

    4 結(jié)論

    根據(jù)最小二乘法的基本原理及前人的研究成果,提出了基于Excel圖解法、INTERCEPT與SLOPE函數(shù)聯(lián)合求解法、TREND函數(shù)法、基于點線距的規(guī)劃求解法和基于擬合誤差的規(guī)劃求解法,重點分析了TREND函數(shù)與基于點線距規(guī)劃求解法2種方法的求解原理,并針對每種求解方法在Excel軟件中設(shè)計了相應(yīng)公式或?qū)崿F(xiàn)方法。通過幾種直剪試驗的求解方法對比分析,可以得出以下結(jié)論:

    1)Excel中的 TREND函數(shù)求解法,截距函數(shù)INTERCEPT與斜率函數(shù) SLOPE聯(lián)合求解法均為Excel內(nèi)部函數(shù)法,大大地減小了編輯公式的工作量,是一種操作簡便計算方法,在無黏性土或計算黏聚力<0的情況下,TREND函數(shù)求解法最為方便。

    2)5種求解方法基本原理均采用基于最小二乘法的回歸分析,其擬合誤差保證其與實測值差值的平方和最小,其精度較采用兩點計算的斜截求解法和手工繪圖更能反映各測點對試驗結(jié)果的影響,更能代表土的實際抗剪強度指標(biāo)。

    3)在5種求解方法中,當(dāng)計算的黏聚力 c值為正值時,各種方法求解結(jié)果相近,均可作為試驗計算結(jié)果。當(dāng)由于某種原因而引起擬合c值為負(fù)值或土樣為無黏性土?xí)r,c=0才正確。Excel圖解法、TREND函數(shù)法、兩種規(guī)劃求解法均可實現(xiàn)強制黏聚力為0,而INTERCEPT與SLOPE函數(shù)聯(lián)合求解法無法實現(xiàn)強制截距為0,故不能采用。

    [1]李向國.滑帶土抗剪強度參數(shù)的確定[J].西部探礦工程,2010(2):37-38.

    [2]馬顯春,王雷,趙法鎖.滑坡穩(wěn)定影響因子敏感性分析及治理方案探討[J].地質(zhì)力學(xué)學(xué)報,2008,14(4):381-388.

    [3]劉遠(yuǎn)明,劉杰,夏才初.不同節(jié)理表面形貌下非貫通節(jié)理巖體強度特性直剪試驗研究[J].巖土力學(xué),2014,35(5): 1269-1273.

    [4]衡帥,楊春和,曾義金,等.基于直剪試驗的頁巖強度各向異性研究[J].巖石力學(xué)與工程學(xué)報,2014,33(5):874-883.

    [5]陳建勝,陳從新,魯祖德,等.強風(fēng)化角巖力學(xué)—變形特性的直剪試驗研究[J].巖土力學(xué),2010,31(9):2869-2874.

    [6]王曉斌,張鐵男,董志國.直剪試驗中的斜截式計算方法[J].水利科技與經(jīng)濟(jì),1999,5(2):82-83.

    [7]南京水利科學(xué)研究院土工研究所.土工試驗技術(shù)手冊[M].北京:人民交通出版社,2003.

    [8]談云志,王世梅.直剪試驗數(shù)據(jù)處理方法的對比分析[J].三峽大學(xué)學(xué)報(自然科學(xué)版),2005,27(2):132-134.

    [9]林立軍.二次曲線擬合法在無機(jī)結(jié)合料穩(wěn)定材料擊實試驗數(shù)據(jù)處理中的應(yīng)用[J].福建建筑,2012,30(5):83-85.

    [10]余東明,姚海林,吳少鋒.三軸試驗抗剪強度參數(shù)值回歸分析法的區(qū)別與修正[J].巖土力學(xué),2012,33(10):3037-3042.

    [11]陳輝,顏勝才,孟軍濤.K30試驗數(shù)據(jù)的曲線擬合分析方法[J].鐵道建筑,2012(11):67-70.

    [12]阮波,張向京,彭意.Excel規(guī)劃求解三軸試驗抗剪強度指標(biāo)[J].鐵道科學(xué)與工程學(xué)報,2009,6(5):57-60.

    (責(zé)任審編 李付軍)

    TU411

    :ADOI:10.3969/j.issn.1003-1995.2015.09.29

    2014-11-03;

    :2015-02-10

    國家自然科學(xué)基金項目(51208198);江西省教育廳項目(JXJG-13-5-15);華東交通大學(xué)科研基金(12RJ02)

    莫林利(1977— ),女,河南安陽人,講師,碩士。

    1003-1995(2015)09-0102-04

    猜你喜歡
    黏聚力土樣抗剪
    季凍黃土地區(qū)框錨邊坡預(yù)應(yīng)力錨桿內(nèi)力計算方法
    灌區(qū)渠道基土工程水敏性試驗研究
    檸檬酸對改良紫色土中老化銅的淋洗研究
    配合比對三合土抗剪強度影響的試驗研究
    膨脹土干濕交替作用下殘余強度試驗方案分析
    治淮(2018年6期)2018-01-30 11:42:44
    槽和黏層油對瀝青混合料層間抗剪性能的影響
    黏聚力強度對滑面作用的差異分析
    科技視界(2016年7期)2016-04-01 20:15:33
    PVA-ECC抗剪加固帶懸臂RC梁承載力計算研究
    鋼-混凝土組合梁開孔板連接件抗剪承載力計算研究
    基于準(zhǔn)黏聚力原理的加筋黏土強度理論淺析
    四川建筑(2014年6期)2014-09-03 11:00:08
    欧美三级亚洲精品| 又爽又黄a免费视频| 看黄色毛片网站| 熟女电影av网| 国产精品免费一区二区三区在线| 免费人成视频x8x8入口观看| 欧美xxxx性猛交bbbb| 又黄又爽又刺激的免费视频.| 又黄又爽又免费观看的视频| 亚洲精品一区av在线观看| 亚洲熟妇熟女久久| 亚洲国产色片| 欧美潮喷喷水| 在线观看av片永久免费下载| 长腿黑丝高跟| 99riav亚洲国产免费| 国产69精品久久久久777片| 啦啦啦韩国在线观看视频| 免费无遮挡裸体视频| 久久国产乱子伦精品免费另类| 精品久久国产蜜桃| 熟女人妻精品中文字幕| 国产色婷婷99| 免费在线观看影片大全网站| 五月玫瑰六月丁香| 欧美成人a在线观看| 婷婷丁香在线五月| 久久久久久久久久成人| 亚洲国产精品sss在线观看| 国产精品爽爽va在线观看网站| 亚洲欧美清纯卡通| 午夜激情福利司机影院| АⅤ资源中文在线天堂| 亚洲av免费在线观看| 久久这里只有精品中国| 久久久久久久久久黄片| 看片在线看免费视频| 精品熟女少妇八av免费久了| 三级毛片av免费| 成年免费大片在线观看| 免费看日本二区| 九色国产91popny在线| 在线播放无遮挡| 中文资源天堂在线| 国产一级毛片七仙女欲春2| 亚洲国产精品999在线| 少妇的逼好多水| 久久久久久久精品吃奶| 亚洲av免费高清在线观看| 国产美女午夜福利| 国产成人啪精品午夜网站| 成人一区二区视频在线观看| 白带黄色成豆腐渣| 好男人电影高清在线观看| 搡女人真爽免费视频火全软件 | 99热6这里只有精品| 99久久成人亚洲精品观看| 免费人成在线观看视频色| 久久久久免费精品人妻一区二区| 国产黄片美女视频| av天堂在线播放| 亚洲自拍偷在线| 少妇人妻一区二区三区视频| 亚洲乱码一区二区免费版| eeuss影院久久| www.色视频.com| 长腿黑丝高跟| 日本a在线网址| 狠狠狠狠99中文字幕| 性色av乱码一区二区三区2| 国产黄色小视频在线观看| 俄罗斯特黄特色一大片| 一区二区三区免费毛片| 一级av片app| 国产伦精品一区二区三区四那| 国产免费男女视频| 亚洲电影在线观看av| 亚洲精品一区av在线观看| 18禁黄网站禁片免费观看直播| 搡女人真爽免费视频火全软件 | 高清毛片免费观看视频网站| 亚洲欧美激情综合另类| 国产成人啪精品午夜网站| 在线播放无遮挡| 我的老师免费观看完整版| 丰满的人妻完整版| 欧美zozozo另类| 少妇裸体淫交视频免费看高清| 国产成人影院久久av| 精品99又大又爽又粗少妇毛片 | 久久久久久久久久成人| 亚洲电影在线观看av| 亚洲国产精品sss在线观看| 搡女人真爽免费视频火全软件 | 最新在线观看一区二区三区| 免费在线观看成人毛片| 日本黄色片子视频| 日韩精品中文字幕看吧| 天天一区二区日本电影三级| 变态另类成人亚洲欧美熟女| 好男人电影高清在线观看| 亚洲真实伦在线观看| 亚洲激情在线av| 香蕉av资源在线| 美女高潮喷水抽搐中文字幕| 十八禁网站免费在线| 琪琪午夜伦伦电影理论片6080| 观看美女的网站| 国产精品久久久久久人妻精品电影| 欧美成人a在线观看| 久久久久久久精品吃奶| 嫁个100分男人电影在线观看| 久久性视频一级片| 夜夜夜夜夜久久久久| 成人性生交大片免费视频hd| 午夜福利成人在线免费观看| 91麻豆av在线| 中文字幕久久专区| 国产毛片a区久久久久| 欧美一区二区精品小视频在线| 欧美bdsm另类| 亚洲av日韩精品久久久久久密| 一本久久中文字幕| 精品一区二区三区视频在线| bbb黄色大片| 免费在线观看成人毛片| 蜜桃亚洲精品一区二区三区| 日韩欧美在线二视频| 亚洲av不卡在线观看| 窝窝影院91人妻| 亚洲av五月六月丁香网| 亚洲精品亚洲一区二区| 国产主播在线观看一区二区| 色吧在线观看| 九九热线精品视视频播放| 一边摸一边抽搐一进一小说| 一区二区三区四区激情视频 | 婷婷色综合大香蕉| 精品一区二区免费观看| 丰满乱子伦码专区| 亚洲国产精品成人综合色| 无人区码免费观看不卡| 午夜久久久久精精品| 欧美绝顶高潮抽搐喷水| 神马国产精品三级电影在线观看| 亚洲精品色激情综合| 免费观看人在逋| 男人和女人高潮做爰伦理| 精品99又大又爽又粗少妇毛片 | 欧美乱色亚洲激情| 成人毛片a级毛片在线播放| 国产欧美日韩一区二区三| 欧美成人性av电影在线观看| 99热这里只有是精品在线观看 | 国产av不卡久久| 日韩高清综合在线| 人妻久久中文字幕网| 在线国产一区二区在线| 日韩免费av在线播放| 91麻豆精品激情在线观看国产| 欧美成人性av电影在线观看| 五月伊人婷婷丁香| 五月玫瑰六月丁香| 色精品久久人妻99蜜桃| 免费观看的影片在线观看| 久久精品人妻少妇| 国产毛片a区久久久久| 午夜两性在线视频| 精品人妻1区二区| 国产主播在线观看一区二区| 欧美zozozo另类| 校园春色视频在线观看| 伊人久久精品亚洲午夜| 精品一区二区三区视频在线观看免费| 成年女人看的毛片在线观看| av欧美777| 露出奶头的视频| 国产黄色小视频在线观看| 欧美乱妇无乱码| 久久草成人影院| 亚洲男人的天堂狠狠| 丁香六月欧美| 99久久久亚洲精品蜜臀av| 久久久久久久久久黄片| av视频在线观看入口| 亚洲欧美精品综合久久99| 国产精品人妻久久久久久| 国产精品综合久久久久久久免费| 91字幕亚洲| 国产91精品成人一区二区三区| 男女下面进入的视频免费午夜| 淫秽高清视频在线观看| 最后的刺客免费高清国语| 蜜桃久久精品国产亚洲av| .国产精品久久| 久久精品国产亚洲av香蕉五月| 成人无遮挡网站| 国产精品久久电影中文字幕| 免费无遮挡裸体视频| 3wmmmm亚洲av在线观看| 女人十人毛片免费观看3o分钟| 非洲黑人性xxxx精品又粗又长| 看免费av毛片| 日本黄色视频三级网站网址| 麻豆成人av在线观看| 搡女人真爽免费视频火全软件 | 极品教师在线视频| 国产成+人综合+亚洲专区| 国产午夜精品久久久久久一区二区三区 | 如何舔出高潮| 久久伊人香网站| 白带黄色成豆腐渣| 久久国产乱子伦精品免费另类| 床上黄色一级片| 99在线视频只有这里精品首页| 久久久久精品国产欧美久久久| 国产午夜精品论理片| 日本黄大片高清| 亚洲自偷自拍三级| 美女cb高潮喷水在线观看| 亚洲天堂国产精品一区在线| 在线观看免费视频日本深夜| 久久精品国产亚洲av香蕉五月| 亚洲精品久久国产高清桃花| 中文在线观看免费www的网站| av福利片在线观看| 婷婷六月久久综合丁香| 一区二区三区免费毛片| 91麻豆av在线| 成人性生交大片免费视频hd| 中文在线观看免费www的网站| 欧美日韩综合久久久久久 | 天堂网av新在线| 好男人电影高清在线观看| 欧美日韩乱码在线| 男人狂女人下面高潮的视频| 黄色一级大片看看| 免费av毛片视频| 一边摸一边抽搐一进一小说| 亚洲人成伊人成综合网2020| 哪里可以看免费的av片| 亚洲自偷自拍三级| 久99久视频精品免费| 国产精品人妻久久久久久| 久久国产乱子伦精品免费另类| 啪啪无遮挡十八禁网站| 最近最新中文字幕大全电影3| 国产爱豆传媒在线观看| 欧美精品国产亚洲| 国产成人aa在线观看| 成人精品一区二区免费| 精品人妻1区二区| 男女之事视频高清在线观看| 精品99又大又爽又粗少妇毛片 | 又黄又爽又免费观看的视频| 国产不卡一卡二| 色播亚洲综合网| 中出人妻视频一区二区| 精品久久久久久久人妻蜜臀av| 久久久久亚洲av毛片大全| a级毛片a级免费在线| 热99re8久久精品国产| 午夜免费男女啪啪视频观看 | 午夜久久久久精精品| 夜夜看夜夜爽夜夜摸| 日本精品一区二区三区蜜桃| 国产欧美日韩一区二区精品| 99精品在免费线老司机午夜| 波多野结衣高清无吗| 久久久久久国产a免费观看| 国产亚洲欧美98| 757午夜福利合集在线观看| 久久精品国产亚洲av香蕉五月| 欧美乱妇无乱码| 国产精品野战在线观看| 亚洲成人精品中文字幕电影| 久久久久久久精品吃奶| 国产aⅴ精品一区二区三区波| 午夜福利在线观看免费完整高清在 | 国产成人av教育| 亚洲一区高清亚洲精品| 日日夜夜操网爽| 亚洲成人免费电影在线观看| 精品久久久久久久久av| 欧美在线一区亚洲| 国产精品久久电影中文字幕| 99视频精品全部免费 在线| av在线蜜桃| 乱人视频在线观看| 亚洲人成网站在线播| 99国产精品一区二区三区| 在线观看一区二区三区| 亚洲五月天丁香| 国产色婷婷99| 日本黄色视频三级网站网址| 亚洲精品成人久久久久久| 最后的刺客免费高清国语| 午夜免费男女啪啪视频观看 | 国产精品一区二区三区四区久久| 午夜激情欧美在线| 久久精品国产亚洲av涩爱 | 国产精品精品国产色婷婷| 国产欧美日韩一区二区精品| 精品国内亚洲2022精品成人| 两个人的视频大全免费| 免费无遮挡裸体视频| 中文字幕av成人在线电影| 欧美在线一区亚洲| 日本三级黄在线观看| 国产精品98久久久久久宅男小说| 少妇的逼好多水| 成人鲁丝片一二三区免费| 老熟妇乱子伦视频在线观看| 此物有八面人人有两片| av在线天堂中文字幕| 99热6这里只有精品| 欧美精品啪啪一区二区三区| 午夜a级毛片| 国产亚洲精品久久久com| 日本成人三级电影网站| 日韩人妻高清精品专区| 免费搜索国产男女视频| 久久久精品大字幕| 最好的美女福利视频网| 亚洲黑人精品在线| 亚洲三级黄色毛片| 亚洲成人精品中文字幕电影| 国产精品精品国产色婷婷| 亚洲av电影不卡..在线观看| 一级黄色大片毛片| 91狼人影院| 国产av在哪里看| 天美传媒精品一区二区| 美女 人体艺术 gogo| 日本 欧美在线| 欧美国产日韩亚洲一区| 国产一区二区三区视频了| 午夜亚洲福利在线播放| 99久久成人亚洲精品观看| 人妻制服诱惑在线中文字幕| 村上凉子中文字幕在线| 欧美xxxx黑人xx丫x性爽| 国产色爽女视频免费观看| 日韩免费av在线播放| a级毛片a级免费在线| 欧美黄色淫秽网站| 全区人妻精品视频| 亚洲最大成人av| 欧美色欧美亚洲另类二区| 国产欧美日韩精品一区二区| 久久国产精品影院| 在线观看66精品国产| 床上黄色一级片| 特大巨黑吊av在线直播| 熟女电影av网| 国产一区二区在线av高清观看| 欧美色欧美亚洲另类二区| 亚洲av二区三区四区| 变态另类丝袜制服| 久久精品国产99精品国产亚洲性色| .国产精品久久| 国产一区二区亚洲精品在线观看| 国产单亲对白刺激| 全区人妻精品视频| 国产精品爽爽va在线观看网站| 亚洲人成网站在线播| 12—13女人毛片做爰片一| 午夜亚洲福利在线播放| 久久人妻av系列| 最好的美女福利视频网| 成人无遮挡网站| 日韩欧美在线乱码| 草草在线视频免费看| 一本一本综合久久| 老司机深夜福利视频在线观看| 亚洲,欧美精品.| 久久99热6这里只有精品| 99视频精品全部免费 在线| 男女视频在线观看网站免费| 亚洲av五月六月丁香网| 国模一区二区三区四区视频| 国产视频一区二区在线看| 亚洲精华国产精华精| 亚洲色图av天堂| 成人特级av手机在线观看| 亚洲avbb在线观看| 免费在线观看亚洲国产| 在线国产一区二区在线| 波多野结衣高清无吗| 免费看美女性在线毛片视频| 国产成人啪精品午夜网站| 夜夜看夜夜爽夜夜摸| 欧美日韩乱码在线| 中文字幕人妻熟人妻熟丝袜美| 日韩免费av在线播放| 亚洲精品影视一区二区三区av| 丁香六月欧美| 精品久久久久久久久亚洲 | 国产成年人精品一区二区| 色视频www国产| 国产精品乱码一区二三区的特点| 色吧在线观看| 国内精品久久久久久久电影| 又紧又爽又黄一区二区| 亚洲熟妇中文字幕五十中出| 成年女人毛片免费观看观看9| 一区二区三区激情视频| 我要搜黄色片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 91麻豆精品激情在线观看国产| 18禁黄网站禁片午夜丰满| 欧美精品啪啪一区二区三区| 国产国拍精品亚洲av在线观看| 麻豆av噜噜一区二区三区| 欧美日韩黄片免| 欧美成人a在线观看| 亚洲专区国产一区二区| 亚洲av免费高清在线观看| 精品一区二区三区人妻视频| 热99re8久久精品国产| 大型黄色视频在线免费观看| 欧美丝袜亚洲另类 | 看免费av毛片| 免费看日本二区| 免费一级毛片在线播放高清视频| 亚洲片人在线观看| 欧美xxxx性猛交bbbb| 女同久久另类99精品国产91| 99在线人妻在线中文字幕| 国产精品永久免费网站| 99国产极品粉嫩在线观看| 亚洲精品久久国产高清桃花| 成人永久免费在线观看视频| 久久精品影院6| 黄色丝袜av网址大全| 人人妻,人人澡人人爽秒播| 欧美成人a在线观看| 91在线观看av| 亚洲av美国av| 欧美另类亚洲清纯唯美| 嫩草影视91久久| 欧美成人一区二区免费高清观看| 午夜免费男女啪啪视频观看 | 变态另类丝袜制服| 国产精品精品国产色婷婷| 一个人看视频在线观看www免费| 欧美zozozo另类| 麻豆成人av在线观看| 亚洲最大成人手机在线| 久久欧美精品欧美久久欧美| 蜜桃亚洲精品一区二区三区| 婷婷色综合大香蕉| 久久草成人影院| 国产精品影院久久| 午夜精品久久久久久毛片777| 成人性生交大片免费视频hd| 日日摸夜夜添夜夜添小说| av欧美777| 国产中年淑女户外野战色| 亚洲中文字幕日韩| 亚洲黑人精品在线| 国产精品99久久久久久久久| 97超视频在线观看视频| 国产精品爽爽va在线观看网站| 一本精品99久久精品77| 嫩草影院精品99| 欧美丝袜亚洲另类 | 日本一二三区视频观看| 亚洲国产日韩欧美精品在线观看| 极品教师在线免费播放| 91狼人影院| 国产一区二区激情短视频| 国产大屁股一区二区在线视频| 男女下面进入的视频免费午夜| 国产精品久久久久久久久免 | 免费大片18禁| 免费在线观看日本一区| 床上黄色一级片| 在线播放国产精品三级| 免费看a级黄色片| 日本撒尿小便嘘嘘汇集6| 亚洲专区中文字幕在线| 青草久久国产| .国产精品久久| 最后的刺客免费高清国语| 在线天堂最新版资源| 夜夜爽天天搞| 亚洲精品一卡2卡三卡4卡5卡| 久久天躁狠狠躁夜夜2o2o| 伦理电影大哥的女人| 国产精品自产拍在线观看55亚洲| 可以在线观看毛片的网站| 少妇的逼好多水| 久久久久久久亚洲中文字幕 | 又黄又爽又刺激的免费视频.| 精品国产三级普通话版| 日本黄色片子视频| ponron亚洲| 禁无遮挡网站| ponron亚洲| 国产美女午夜福利| 夜夜看夜夜爽夜夜摸| 久久午夜福利片| 亚洲色图av天堂| 中国美女看黄片| 在线观看66精品国产| 俺也久久电影网| 国产精品1区2区在线观看.| 99riav亚洲国产免费| 波多野结衣巨乳人妻| 亚州av有码| 十八禁国产超污无遮挡网站| 久久精品国产自在天天线| 91麻豆av在线| 亚洲美女视频黄频| 九色国产91popny在线| 草草在线视频免费看| 精品午夜福利视频在线观看一区| 99国产极品粉嫩在线观看| 亚洲人与动物交配视频| 亚洲国产高清在线一区二区三| 免费观看的影片在线观看| 日韩免费av在线播放| 18禁黄网站禁片午夜丰满| 日韩欧美 国产精品| 成人av在线播放网站| 不卡一级毛片| 男插女下体视频免费在线播放| 中文亚洲av片在线观看爽| 欧美乱妇无乱码| 精品国产亚洲在线| 欧美xxxx性猛交bbbb| 在线观看66精品国产| 成人无遮挡网站| 成年版毛片免费区| 国产av不卡久久| 国产单亲对白刺激| x7x7x7水蜜桃| 欧美黄色淫秽网站| 97人妻精品一区二区三区麻豆| 成年免费大片在线观看| 久久久久亚洲av毛片大全| 国产成人aa在线观看| 国产综合懂色| 少妇熟女aⅴ在线视频| 我要看日韩黄色一级片| 午夜日韩欧美国产| 国产精品精品国产色婷婷| av黄色大香蕉| 国产精品99久久久久久久久| 色哟哟·www| 成人一区二区视频在线观看| 哪里可以看免费的av片| 可以在线观看毛片的网站| 波多野结衣高清作品| 简卡轻食公司| 国产毛片a区久久久久| 亚洲精华国产精华精| 欧美三级亚洲精品| av女优亚洲男人天堂| av天堂在线播放| 麻豆av噜噜一区二区三区| 国产不卡一卡二| 一本精品99久久精品77| 久久精品影院6| 91九色精品人成在线观看| 欧美zozozo另类| 久久精品91蜜桃| 精品人妻熟女av久视频| 精品人妻1区二区| 亚洲成a人片在线一区二区| 国产野战对白在线观看| 99精品久久久久人妻精品| 免费观看精品视频网站| 我的老师免费观看完整版| 老鸭窝网址在线观看| 亚洲va日本ⅴa欧美va伊人久久| 欧美zozozo另类| 波野结衣二区三区在线| 女人十人毛片免费观看3o分钟| 老熟妇乱子伦视频在线观看| 中亚洲国语对白在线视频| 女同久久另类99精品国产91| 一本一本综合久久| 大型黄色视频在线免费观看| 在线看三级毛片| 3wmmmm亚洲av在线观看| 99久久精品国产亚洲精品| 久久精品国产亚洲av天美| av国产免费在线观看| 在线观看午夜福利视频| 午夜福利在线观看免费完整高清在 | 最后的刺客免费高清国语| 中文亚洲av片在线观看爽| 精品久久久久久久末码| 麻豆成人av在线观看| 久9热在线精品视频| 大型黄色视频在线免费观看| 国产精品嫩草影院av在线观看 | 99热这里只有是精品50| АⅤ资源中文在线天堂| 天堂网av新在线| 超碰av人人做人人爽久久| 一个人看的www免费观看视频| 97热精品久久久久久| 国产免费一级a男人的天堂| 美女 人体艺术 gogo| 天天一区二区日本电影三级| 午夜a级毛片| 精品人妻一区二区三区麻豆 | 757午夜福利合集在线观看| 男女之事视频高清在线观看| 少妇高潮的动态图| 国产精品嫩草影院av在线观看 | 两性午夜刺激爽爽歪歪视频在线观看|