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

    基于混合粒子群算法和能量法主動(dòng)土壓力計(jì)算

    2012-01-08 07:12:16陳昌富唐仁華梁冠亭
    巖土力學(xué) 2012年6期
    關(guān)鍵詞:螺線擋土墻對(duì)數(shù)

    陳昌富,唐仁華,梁冠亭

    (湖南大學(xué) 巖土工程研究所,長(zhǎng)沙 410082)

    1 引 言

    自Terzaghi[1]使用對(duì)數(shù)螺線滑動(dòng)面模型求解擋土墻土壓力后,此模型經(jīng)歷了幾十年試驗(yàn)和實(shí)踐的考證,被普遍認(rèn)為是最接近實(shí)際的破壞機(jī)構(gòu)[2]。Chen[3]較早地系統(tǒng)研究了極限分析法在土壓力中的應(yīng)用,隨后許多學(xué)者運(yùn)用塑性極限上限方法來(lái)計(jì)算土壓力,取得了較多成果。范文等[4]和楊建民[5]利用多三角形破壞機(jī)構(gòu),分別基于統(tǒng)一強(qiáng)度理論和非線性強(qiáng)度準(zhǔn)則求得了土壓力的上限解。楊小禮[6]運(yùn)用非線性強(qiáng)度準(zhǔn)則分析了墻后填土面水平時(shí)動(dòng)、靜主動(dòng)土壓力上限解。Chen[3]分別采用6 種破壞機(jī)構(gòu)結(jié)合最速下降法求得了土壓力系數(shù),其中2 條直線中間夾對(duì)數(shù)螺線破壞機(jī)構(gòu)的計(jì)算結(jié)果與廣泛接受的Sokoloskii[7]解誤差最小,Chen 認(rèn)為,是對(duì)數(shù)螺線夾層起到了關(guān)鍵作用。但上述機(jī)構(gòu)假定對(duì)數(shù)螺線的中心固定于坡頂不變,而實(shí)際情況并非如此。本文采用接近實(shí)際情況的單純對(duì)數(shù)螺線滑動(dòng)面破壞機(jī)構(gòu),不約束對(duì)數(shù)螺線中心的位置,推導(dǎo)出土壓力及其系數(shù)計(jì)算公式,然后運(yùn)用基于自然選擇的混合粒子群算法對(duì)最危險(xiǎn)滑動(dòng)面進(jìn)行全局搜索并計(jì)算主動(dòng)土壓力,最后將本文方法計(jì)算結(jié)果與已有解答和實(shí)測(cè)值進(jìn)行對(duì)比分析。

    2 基于能量法和對(duì)數(shù)螺線滑裂面的主動(dòng)土壓力計(jì)算

    能量法計(jì)算土壓力需用到以下兩個(gè)結(jié)論[8]:

    (1)流動(dòng)法則:土體處于塑性流動(dòng)或剪切滑動(dòng)狀態(tài)時(shí),滑動(dòng)面上任一點(diǎn)的應(yīng)變速度矢量v 與該點(diǎn)處的滑動(dòng)線成φ 角。

    (2)土體在塑性流動(dòng)中的能量消散為cvcosφ。c、φ 分別為土體的黏聚力和內(nèi)摩擦角。

    如圖1 所示的擋土墻,墻高為H,墻背與水平面成β 角,墻后填土與水平面成α 角。填土上作用均布荷載q。以坡頂處A 點(diǎn)為坐標(biāo)原點(diǎn),向左和向上分別定為x,y 坐標(biāo)軸的正方向。

    圖1 對(duì)數(shù)螺線滑動(dòng)面計(jì)算圖 Fig.1 Calculation diagram of log-spiral sliding surface

    假定滑動(dòng)面為通過(guò)墻底的對(duì)數(shù)螺線滑動(dòng)面,其動(dòng)徑方程為

    式中:rθ為0r 繞中心旋轉(zhuǎn)θ 角后的動(dòng)徑長(zhǎng)度(m);0θ為0r 與水平面的夾角(°);0r 為對(duì)數(shù)螺線通過(guò)邊坡的旋入點(diǎn)到螺線中心的長(zhǎng)度(m)。

    當(dāng)θ =hθ 時(shí),

    式中:hr 為對(duì)數(shù)螺線通過(guò)邊坡的旋出點(diǎn)至螺線中心的長(zhǎng)度(m);hθ 為hr 與水平面的夾角(°)。

    過(guò)O 點(diǎn)作AB 的平行線與CA 的延長(zhǎng)線交與點(diǎn)E;作AC 的平行線與BA 的延長(zhǎng)線交于點(diǎn)F,則

    由式(3)可以推出

    此外,對(duì)數(shù)螺線參數(shù)還滿足式(5)中的幾何關(guān)系

    式中:xo為對(duì)數(shù)螺線中心橫坐標(biāo)(m);yo為對(duì)數(shù)螺線中心縱坐標(biāo)(m);L 為滑裂面起點(diǎn)到坡頂?shù)木嚯x;xAB為L(zhǎng) 在x 軸負(fù)方向的投影。

    作用在滑動(dòng)土體上的外力有均布荷載載q,滑動(dòng)土體的重力W 已及擋土墻對(duì)滑動(dòng)土體的反作用力p。假設(shè)滑動(dòng)土體繞對(duì)數(shù)螺線中心產(chǎn)生微小的角速度ω,由假定的對(duì)數(shù)螺線方程可知,螺線上每點(diǎn)的切線方向與動(dòng)徑成(π/2 φ+ )角,結(jié)合流動(dòng)法則,則滑動(dòng)面上每點(diǎn)的應(yīng)變速度方向均與動(dòng)徑垂直,即

    如圖2、3 所示,由于超載和重力的方向都是豎直向下,故其與速度方向夾角等于動(dòng)徑與水平線的夾角。即超載或重力做的功等于力與角速度ω 的乘積,再乘以超載或者單元土體形心到對(duì)數(shù)螺線中心O 點(diǎn)的水平距離。

    圖2 滑體發(fā)生微小轉(zhuǎn)動(dòng)計(jì)算圖 Fig.2 Calculation diagram of small rotation on slide soil

    圖3 滑體發(fā)生微小轉(zhuǎn)動(dòng)計(jì)算圖 Fig.3 Calculation diagram of small rotation on slide soil

    超載q 所作的功為

    如圖4 所示,滑動(dòng)土體重力作的功可由疊加法求出。以W1表示OBC 對(duì)假定轉(zhuǎn)動(dòng)中心O 作的功,W2表示OAB 對(duì)假定轉(zhuǎn)動(dòng)中心O 作的功,W3表示OAC 對(duì)假定轉(zhuǎn)動(dòng)中心O 作的功。OAB 形心到O 點(diǎn)的水平距離可由OAM 與ABM 組合求出,OAC 形心到O 點(diǎn)的水平距離可由OCN 與ACN 組合求出。則W1,W2,W3分別為

    沿對(duì)數(shù)螺線微分長(zhǎng)度 rθdθ / cosφ 上的消耗功為( rθdθ / cosφ ) cv cosφ = crθvd θ,內(nèi)部總消耗功為

    圖4 AC 面上單元土體速度矢量圖 Fig.4 Soil element’s velocity vector on surface AC

    如圖4(a),由余弦定理知

    由正弦定理知

    式 中:A=sin-1[ H sin(β + θh)/(2 ozsin β)];G = sin(β +θh)/(2sin β);J = ( r0/ H )2e2(θh-θ0)tanφ;U=cos(β +θ )/ H sin β。 h

    圖4(b)右為AC 面上單元土體速度矢量合成圖,可見(jiàn)其平均速度即為中點(diǎn)Z 處的速度。土壓力合力P 與墻面法線夾角為δ (外摩擦角),位于法線下方,與Z 點(diǎn)速度方向夾角為π-(hθ -zθ )-β -δ ,則擋土墻對(duì)AC 面土體的反力所做功為

    由于擋土墻對(duì)土體的反力與土體的應(yīng)變速率方向相反,故其所做之功為負(fù)值。由外力作的功等于土體內(nèi)部消耗的能量,消去角速度可得:

    式(17)即為主動(dòng)土壓力計(jì)算公式。顯然,這是一個(gè)二元函數(shù)的極值問(wèn)題,對(duì)于主動(dòng)土壓力,應(yīng)取其最大值。

    特別地,當(dāng)α =0°、β =90°,即墻背豎直、墻后填土面水平時(shí),整理公式后可得:

    式中:

    不同深度處的主動(dòng)土壓力強(qiáng)度pH為

    主動(dòng)土壓力的作用點(diǎn)距墻底的高度為

    式中:Pa為主動(dòng)土壓力,kaγ,kac,kaq分別為單獨(dú)考慮重度,黏聚力,超載時(shí)的土壓力系數(shù)。

    3 基于混合粒子群算法主動(dòng)土壓力優(yōu)化計(jì)算

    基于能量法主動(dòng)土壓力計(jì)算可歸結(jié)為如下優(yōu)化模型:

    式中:X 為優(yōu)化設(shè)計(jì)變量;f (X)為目標(biāo)函數(shù),計(jì)算土壓力時(shí)f (X)=-Pa,計(jì)算土壓力系數(shù)時(shí)f (X)=-kaγ,其中Pa與kaγ分別按式(17)與式(20)計(jì)算;gj(X)( j=1, 2, 3,4)為約束條件,其中xAB按式(5)計(jì)算。

    式(25)是一個(gè)含有復(fù)雜的目標(biāo)函數(shù)和約束條件的非線性規(guī)劃問(wèn)題,本文采用適用性強(qiáng)且具有良好的全局搜索能力的粒子群算法[9](PSO)求解。

    在PSO 算法中,每個(gè)粒子都具有3 個(gè)參數(shù):位置矢量X,代表目標(biāo)函數(shù)中自變量的值;速度矢量V,代表粒子飛翔的力向和距離,用來(lái)引導(dǎo)粒子不斷朝著更好的解空間改變當(dāng)前位置;適應(yīng)度f(wàn)itness,由被優(yōu)化的目標(biāo)函數(shù)決定,目標(biāo)函數(shù)值越小,則適應(yīng)度越好。在每一次迭代中,粒子通過(guò)跟蹤2 個(gè)最優(yōu)解來(lái)更新自己,一是粒子本身所找到的最優(yōu)解,另一個(gè)是整個(gè)種群目前所找到的最優(yōu)解。

    混合粒子群算法是將自然選擇機(jī)制與粒子群算法結(jié)合,在每次迭代過(guò)程中將整個(gè)粒子群按適應(yīng)值排序,用群體中最好的一半的粒子的速度和位置替換最差的一半的位置和速度,同時(shí)保留原來(lái)每個(gè)個(gè)體所記憶的歷史最優(yōu)值。

    基于以上原理和模型,主動(dòng)土壓力混合粒子群優(yōu)化算法計(jì)算步驟如下:

    (1)設(shè)定粒子群數(shù)目為40,即i∈[1,40]。對(duì)每個(gè)粒子的當(dāng)前位置矢量 Xi=[xi,1,xi,2]和速度矢量Vi= [vi,1,vi,2]在變量范圍內(nèi)[0, π-β ]隨機(jī)賦予初始值;

    (2)根據(jù)每個(gè)粒子當(dāng)前的iX 計(jì)算其適應(yīng)度。若gi(X)>0 成立,則fitness(i)=f (iX );否則,將一個(gè)遠(yuǎn)大于正常土壓值的大數(shù)賦給適應(yīng)度函數(shù)作為此粒子的適應(yīng)值以使其受到懲罰,例如fitness(i)=10 000;

    (3)將各微粒的位置和適應(yīng)值存儲(chǔ)在各微粒的個(gè)體最好位置pbest 中,將所有pbest 中適應(yīng)值最優(yōu)個(gè)體的位置和適應(yīng)值存儲(chǔ)于群體最好位置gbest 中;

    (4)按式(26)和(27)更新每個(gè)微粒的速度和位置:

    式中:t 為迭代的步數(shù);pi為個(gè)體最優(yōu)適應(yīng)值對(duì)應(yīng)的[xi,1, xi,2]值,即第i 個(gè)粒子最優(yōu)適應(yīng)值對(duì)應(yīng)的的[0θ ,hθ ]值;pg為群體最優(yōu)適應(yīng)值對(duì)應(yīng)的[0θ ,hθ ]值;w 為慣性權(quán)重因子,文中取w=0.5;c1、c2為正的學(xué)習(xí)因子,文中取c1= c2=2;r1、r2為[0,1]范圍內(nèi)的均勻隨機(jī)數(shù)。

    (5)對(duì)每個(gè)微粒,將其適應(yīng)值與其經(jīng)歷過(guò)的最好適應(yīng)值比較,如果較好,則將其對(duì)應(yīng)的位置作為當(dāng)前的最好位置,然后比較當(dāng)前所有pbest 和gbest的值,更新gbest;

    (6)將整個(gè)粒子群按適應(yīng)值排序,用群體中最好的一半的粒子的速度和位置替換掉最差的一半的位置和速度,保持pbest 和gbest 不變。

    (7)設(shè)定最大迭代次數(shù)為100 步,達(dá)到迭代次數(shù)后則輸出結(jié)果,即土壓力最大時(shí)的0θ 、hθ 值以及土壓力或者土壓力系數(shù)的值。

    4 結(jié)果對(duì)比分析

    當(dāng)墻背豎直時(shí),對(duì)于無(wú)超載的砂性土,將c =0,q=0,β =90°代入式(20),運(yùn)用上述優(yōu)化方法計(jì)算土壓力系數(shù)kaγ隨填土面傾角α 的變化結(jié)果,并與文獻(xiàn)[3]中結(jié)果進(jìn)行對(duì)比。

    由表1 知,在不同的填土面傾角下,本文結(jié)果與文獻(xiàn)[3]結(jié)果幾乎一致(差值在2%以內(nèi))。說(shuō)明兩者的差別對(duì)填土面傾斜角度不敏感。

    當(dāng)填土面水平時(shí)(α =0°),土壓力系數(shù)kaγ隨墻背傾角β 的變化與文獻(xiàn)[3]結(jié)果對(duì)比如表2 所示。

    表1 墻背豎直時(shí)與文獻(xiàn)[3]結(jié)果對(duì)比 Table 1 Comparison between results obtained by reference [3] and with a vertical wall

    表2 填土面水平時(shí)與文獻(xiàn)[3]結(jié)果對(duì)比 Table 2 Comparison between results obtained reference [3] and with horizontal backfills

    從表2 中選擇部分?jǐn)?shù)據(jù)繪成更直觀的圖5。

    圖5 填土面水平時(shí)差值隨墻面傾角的變化(δ =φ /2) Fig.5 Difference varying with wall angle for horizontal backfills(δ =φ /2)

    從圖可以看出,當(dāng)墻面與水平面的夾角為90°、80°、70°時(shí),兩者差值接近于0。而當(dāng)墻背向內(nèi)傾斜角度大于30°時(shí),兩者存在明顯的差別,尤其是填土摩擦角較大時(shí),本文比文獻(xiàn)[3]的計(jì)算結(jié)果大50%左右。將兩者結(jié)果與朱大勇等[10]的極限平衡解對(duì)比如表3 所示。

    表3 墻背向內(nèi)傾斜30°時(shí)土壓力系數(shù)對(duì)比 Table 3 Comparison of active earth pressure coefficients while the wall angle is 30°

    從表3 中可以看出,當(dāng)墻背傾角較大時(shí),本文比文獻(xiàn)[3]更接近文獻(xiàn)[10]的計(jì)算結(jié)果,筆者認(rèn)為這是由不同的破壞機(jī)構(gòu)導(dǎo)致的。由于文獻(xiàn)[3]沒(méi)有對(duì)對(duì)數(shù)螺線中心進(jìn)行搜索,所以在一些特定情況下不能搜索到土壓力系數(shù)的最優(yōu)值(最大值),從而小于實(shí)際值。顯然,此時(shí)的庫(kù)倫理論精度更小,因?yàn)閴Ρ硟A角較大時(shí),真實(shí)滑動(dòng)面與平面滑動(dòng)面的假定存在較大的差異。

    5 實(shí)例驗(yàn)證

    四川建筑科學(xué)研究所在四川簡(jiǎn)陽(yáng)做的重力式擋土墻試驗(yàn)[11]:墻背豎直,墻高為4 m,土天然重度γ =19.34 kN/m3,不排水三軸快剪指標(biāo)φ = 16.63°、δ =8.32°、c=4.7 kPa。

    計(jì)算結(jié)果見(jiàn)表4。同時(shí)給出了建筑地基基礎(chǔ)規(guī)范[12]推薦方法的計(jì)算結(jié)果與實(shí)測(cè)結(jié)果。

    表4 本文計(jì)算結(jié)果 Table 4 Result of the article

    表5 本文計(jì)算結(jié)果與實(shí)測(cè)值對(duì)比 Table 5 Comparison of results obtained by the article and measured date

    顯然,本文對(duì)黏性土的計(jì)算結(jié)果比規(guī)范推薦方法計(jì)算結(jié)果更接近實(shí)測(cè)值。

    6 結(jié) 論

    (1)本文基于能量法和混合離子群優(yōu)化算法,假定滑動(dòng)面為對(duì)數(shù)螺線滑裂面,得到了一種新的主動(dòng)土壓力及其系數(shù)計(jì)算方法。

    (2)對(duì)于黏性土,通過(guò)與一重力式擋墻試驗(yàn)結(jié)果對(duì)比分析發(fā)現(xiàn),本文結(jié)果與實(shí)測(cè)值誤差僅為5.4%,表明本文方法計(jì)算結(jié)果可靠。

    (3)對(duì)于無(wú)超載砂性土的土壓力系數(shù),本文方法計(jì)算結(jié)果與經(jīng)典的極限分析上限解在墻面傾角較小時(shí)基本一致,但當(dāng)墻面傾角較大時(shí),經(jīng)典解明顯偏小,與本文解最大差別達(dá)54.71%,而本文解與基于最優(yōu)性原理的極限平衡解較接近。

    [1] TERZAGHI K. Theoretical soil mechanics[M]. New York: [s. n.], 1943.

    [2] 趙樹(shù)德, 趙樹(shù)惠. 用對(duì)數(shù)螺線滑動(dòng)面計(jì)算擋土墻主動(dòng)土壓力[J]. 西安建筑科技大學(xué)學(xué)報(bào), 2002, 34(4): 342-345. ZHAO Shu-de, ZHAO Shu-hui. Calculation of earth pressure of active retaining wall with logarithmic spiral sliding surface[J]. Xi'an Univ. of Arch.& Tech., 2002, 34(4): 342-345.

    [3] CHEN W F. Limit analysis and soil plasticity[M]. New York: Elsevier Scientific Publishing Company, 1975.

    [4] 范文, 沈珠江, 愈茂宏. 基于統(tǒng)一強(qiáng)度理論的土壓力極限上限分析[J]. 巖土工程學(xué)報(bào), 2005, 27(10): 1147-1153. FAN Wen, SHENG Zhu-jiang, YU Mao-hong. Upper-bound limit analysis of earth pressure based on unified strength theory[J]. Chinese Journal of Geotechnical Engineering, 2005, 27(10): 1147-1153.

    [5] 楊建民. 基于非線性破壞準(zhǔn)則的主動(dòng)土壓力上限計(jì) 算[J]. 巖土力學(xué), 2009, 30(2): 503-508. YANG Jian-min. Upper bound solution for active earth pressure based on a nonlinear failure criterion[J]. Rock and Soil Mechanics, 2009, 30(2): 503-508.

    [6] YANG Xiao-li. Upper bound limit analysis of active earth pressure with different fracture surfaces and nonlinear yield criterion[J]. Theoretical and Applied Fracture Mechanics, 2007, 47(1): 46-56.

    [7] SOKOLOSKII V V. Statics of granular media[M]. New York: Pergamon Press, 1965.

    [8] 顧慰慈. 擋土墻土壓力計(jì)算手冊(cè)[M]. 北京: 中國(guó)建材出版社, 2005.

    [9] 龔純, 王正林. 精通Matlab 最優(yōu)化計(jì)算[M]. 北京: 電子工業(yè)出版社, 2009.

    [10] 朱大勇, 錢七虎. 極限平衡法計(jì)算土壓力系數(shù)的新途徑[J]. 土木工程學(xué)報(bào), 2002, 33(1): 63-68. ZHU Da-yong, QIAN Qi-hu. New approach for computation of earth pressure coefficients based on limit equilibrium method[J]. China Civil Engineering Journal, 2002, 33(1): 63-68.

    [11] 四川省建筑科學(xué)研究院. 重力式擋土墻的試驗(yàn)研究[R]. 成都: 四川省建筑科學(xué)研究院, 1982.

    [12] 中華人民共和國(guó)國(guó)家標(biāo)準(zhǔn). 建筑地基基礎(chǔ)設(shè)計(jì)規(guī)范[S].北京: 中國(guó)建筑工業(yè)出版社, 2002.

    猜你喜歡
    螺線擋土墻對(duì)數(shù)
    含有對(duì)數(shù)非線性項(xiàng)Kirchhoff方程多解的存在性
    指數(shù)與對(duì)數(shù)
    指數(shù)與對(duì)數(shù)
    三維Minkowski空間中的k-型偽零螺線
    走近等角螺線
    探秘等角螺線
    走近等角螺線
    對(duì)數(shù)簡(jiǎn)史
    淺談選廠常見(jiàn)擋土墻設(shè)計(jì)
    仰斜式重力擋土墻穩(wěn)定計(jì)算復(fù)核
    不卡一级毛片| 女警被强在线播放| 成人国语在线视频| 国产片内射在线| 久久99热这里只频精品6学生| 中文字幕高清在线视频| 久久人妻熟女aⅴ| 极品教师在线免费播放| 丰满饥渴人妻一区二区三| av不卡在线播放| 9热在线视频观看99| 国产一区二区三区视频了| 12—13女人毛片做爰片一| 大型黄色视频在线免费观看| avwww免费| 欧美日韩中文字幕国产精品一区二区三区 | 国产成人一区二区三区免费视频网站| 亚洲国产中文字幕在线视频| 美女福利国产在线| 黄网站色视频无遮挡免费观看| 久久精品国产a三级三级三级| 久久狼人影院| 女人久久www免费人成看片| 日本a在线网址| 欧美成狂野欧美在线观看| 亚洲精品乱久久久久久| 脱女人内裤的视频| av免费在线观看网站| 国产一区有黄有色的免费视频| 黄片小视频在线播放| 国产有黄有色有爽视频| 久久香蕉激情| 伦理电影免费视频| 欧美国产精品va在线观看不卡| 精品福利观看| av线在线观看网站| 成在线人永久免费视频| 亚洲国产欧美一区二区综合| 成人国产av品久久久| 久久婷婷成人综合色麻豆| 欧美+亚洲+日韩+国产| 亚洲精品久久成人aⅴ小说| 久热这里只有精品99| 国产高清激情床上av| 国产三级黄色录像| a级毛片在线看网站| 国产极品粉嫩免费观看在线| 黄片播放在线免费| 亚洲欧美精品综合一区二区三区| 午夜福利视频在线观看免费| 亚洲精品美女久久久久99蜜臀| 香蕉丝袜av| 91国产中文字幕| 母亲3免费完整高清在线观看| 久久精品国产99精品国产亚洲性色 | 亚洲五月婷婷丁香| 久久久国产一区二区| 亚洲精品av麻豆狂野| 亚洲一区中文字幕在线| 无遮挡黄片免费观看| 亚洲第一青青草原| 欧美成狂野欧美在线观看| 熟女少妇亚洲综合色aaa.| 亚洲精品av麻豆狂野| 一级黄色大片毛片| 精品视频人人做人人爽| 精品久久久久久久毛片微露脸| 久久精品91无色码中文字幕| 一级片'在线观看视频| 久久精品亚洲熟妇少妇任你| 极品少妇高潮喷水抽搐| 欧美激情极品国产一区二区三区| 18禁裸乳无遮挡动漫免费视频| 色精品久久人妻99蜜桃| 久久精品成人免费网站| 久久精品成人免费网站| 亚洲色图 男人天堂 中文字幕| 成人亚洲精品一区在线观看| 国产福利在线免费观看视频| 纵有疾风起免费观看全集完整版| 在线看a的网站| 一个人免费在线观看的高清视频| 久久人妻福利社区极品人妻图片| 90打野战视频偷拍视频| 日韩熟女老妇一区二区性免费视频| 天堂动漫精品| 国产午夜精品久久久久久| 国产91精品成人一区二区三区 | 国产av精品麻豆| 97在线人人人人妻| 国产日韩欧美亚洲二区| 99在线人妻在线中文字幕 | 日韩欧美免费精品| 国产av又大| 亚洲第一欧美日韩一区二区三区 | 一级黄色大片毛片| 国产不卡一卡二| 岛国毛片在线播放| 国产色视频综合| 成人黄色视频免费在线看| 国产aⅴ精品一区二区三区波| 精品久久久精品久久久| 美女视频免费永久观看网站| 精品人妻在线不人妻| 人人妻人人澡人人爽人人夜夜| 日韩大码丰满熟妇| 高清av免费在线| 18禁裸乳无遮挡动漫免费视频| 99久久99久久久精品蜜桃| 天堂俺去俺来也www色官网| 国产淫语在线视频| 国产有黄有色有爽视频| 99精品久久久久人妻精品| 免费av中文字幕在线| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美激情在线| 国产日韩欧美在线精品| 久久精品成人免费网站| svipshipincom国产片| 99九九在线精品视频| 亚洲成人免费电影在线观看| 国产真人三级小视频在线观看| 国产一卡二卡三卡精品| 1024视频免费在线观看| 精品第一国产精品| 亚洲三区欧美一区| 国产精品免费一区二区三区在线 | 亚洲精品一二三| 露出奶头的视频| a在线观看视频网站| 亚洲国产欧美在线一区| 香蕉久久夜色| 日韩制服丝袜自拍偷拍| 人妻一区二区av| 桃红色精品国产亚洲av| av天堂久久9| 无限看片的www在线观看| 日韩免费av在线播放| kizo精华| 午夜精品久久久久久毛片777| 大片电影免费在线观看免费| 9191精品国产免费久久| 国产精品影院久久| 日韩人妻精品一区2区三区| 老司机午夜十八禁免费视频| 热re99久久国产66热| 男人舔女人的私密视频| 午夜福利欧美成人| 亚洲av电影在线进入| 亚洲熟女精品中文字幕| 一本综合久久免费| 国产在线观看jvid| 大陆偷拍与自拍| 一区二区三区精品91| 欧美另类亚洲清纯唯美| 老司机在亚洲福利影院| 欧美午夜高清在线| 国产精品偷伦视频观看了| 一级毛片电影观看| 亚洲专区国产一区二区| 午夜福利乱码中文字幕| 亚洲国产精品一区二区三区在线| 亚洲国产欧美在线一区| 亚洲欧美日韩另类电影网站| 欧美老熟妇乱子伦牲交| 电影成人av| 老司机福利观看| 一本综合久久免费| 中文字幕人妻丝袜一区二区| 日本av免费视频播放| av一本久久久久| 在线观看免费视频日本深夜| 久久毛片免费看一区二区三区| 国产在线观看jvid| 久久久久久久久久久久大奶| 亚洲avbb在线观看| 一级黄色大片毛片| 十八禁网站网址无遮挡| 国产欧美日韩综合在线一区二区| 国产99久久九九免费精品| 久久久久久免费高清国产稀缺| 另类亚洲欧美激情| 亚洲精品国产精品久久久不卡| 欧美日韩亚洲综合一区二区三区_| 国产精品成人在线| 日本wwww免费看| 久久人妻av系列| av免费在线观看网站| 精品国产一区二区三区四区第35| 久久国产精品大桥未久av| 香蕉国产在线看| av天堂久久9| 国产精品98久久久久久宅男小说| 丁香欧美五月| 国内毛片毛片毛片毛片毛片| 操出白浆在线播放| 黑人欧美特级aaaaaa片| 热99国产精品久久久久久7| 高潮久久久久久久久久久不卡| av国产精品久久久久影院| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品 国内视频| 午夜激情久久久久久久| 9热在线视频观看99| 国产在线观看jvid| 亚洲一码二码三码区别大吗| 国产人伦9x9x在线观看| 久久婷婷成人综合色麻豆| 国产有黄有色有爽视频| 夫妻午夜视频| 亚洲精品成人av观看孕妇| 亚洲久久久国产精品| 国产色视频综合| 50天的宝宝边吃奶边哭怎么回事| 激情视频va一区二区三区| 宅男免费午夜| 69精品国产乱码久久久| 一本久久精品| 高清视频免费观看一区二区| 最黄视频免费看| 一级毛片精品| 国产不卡一卡二| 一级片免费观看大全| 国产精品亚洲av一区麻豆| 一级毛片女人18水好多| 老司机靠b影院| 亚洲av第一区精品v没综合| 熟女少妇亚洲综合色aaa.| 亚洲精品一卡2卡三卡4卡5卡| 99国产精品免费福利视频| 亚洲一卡2卡3卡4卡5卡精品中文| 女人爽到高潮嗷嗷叫在线视频| 亚洲五月婷婷丁香| 亚洲色图av天堂| 啦啦啦 在线观看视频| 午夜激情久久久久久久| 国产成人欧美| 丝瓜视频免费看黄片| 精品国产一区二区三区久久久樱花| av又黄又爽大尺度在线免费看| 欧美日韩精品网址| av有码第一页| 午夜久久久在线观看| 国产精品一区二区精品视频观看| 狠狠婷婷综合久久久久久88av| 午夜成年电影在线免费观看| 日韩免费av在线播放| 日韩熟女老妇一区二区性免费视频| 女人精品久久久久毛片| 久热爱精品视频在线9| 久久中文看片网| 免费观看a级毛片全部| 99国产精品99久久久久| 国产深夜福利视频在线观看| 如日韩欧美国产精品一区二区三区| 日日摸夜夜添夜夜添小说| 精品人妻1区二区| 一级,二级,三级黄色视频| 久久久精品94久久精品| 久久久久久久国产电影| 日韩视频一区二区在线观看| 啦啦啦中文免费视频观看日本| 人人妻,人人澡人人爽秒播| 超色免费av| 日日夜夜操网爽| 在线永久观看黄色视频| 国产精品美女特级片免费视频播放器 | 精品高清国产在线一区| 国产精品久久久人人做人人爽| 制服诱惑二区| 亚洲欧美一区二区三区久久| 女人久久www免费人成看片| 在线av久久热| 亚洲中文字幕日韩| 久久人妻福利社区极品人妻图片| 国产在线精品亚洲第一网站| 新久久久久国产一级毛片| 国产av又大| 99九九在线精品视频| 一区福利在线观看| tocl精华| 热99国产精品久久久久久7| 国产91精品成人一区二区三区 | av网站在线播放免费| 在线十欧美十亚洲十日本专区| 黄片小视频在线播放| 国产精品秋霞免费鲁丝片| 不卡av一区二区三区| 黄片小视频在线播放| 国产麻豆69| 国产在视频线精品| 亚洲少妇的诱惑av| 国产极品粉嫩免费观看在线| 曰老女人黄片| 老汉色av国产亚洲站长工具| 中文字幕色久视频| 天堂俺去俺来也www色官网| 国产一区二区三区在线臀色熟女 | 亚洲 国产 在线| 午夜福利在线免费观看网站| www.精华液| 国产欧美日韩一区二区三区在线| 男女边摸边吃奶| 日日爽夜夜爽网站| av视频免费观看在线观看| 欧美精品一区二区免费开放| 成人永久免费在线观看视频 | 亚洲精品国产一区二区精华液| 日韩中文字幕欧美一区二区| 肉色欧美久久久久久久蜜桃| 人人澡人人妻人| 十八禁高潮呻吟视频| 午夜福利在线观看吧| 啦啦啦视频在线资源免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品亚洲av国产电影网| 怎么达到女性高潮| 欧美在线黄色| tocl精华| 五月天丁香电影| 亚洲精品美女久久久久99蜜臀| 久久精品国产99精品国产亚洲性色 | 操出白浆在线播放| 日韩一区二区三区影片| 欧美激情极品国产一区二区三区| 久久久久视频综合| 中文亚洲av片在线观看爽 | 国产欧美亚洲国产| 国产亚洲午夜精品一区二区久久| 亚洲 国产 在线| 日本黄色日本黄色录像| 麻豆av在线久日| 正在播放国产对白刺激| 777久久人妻少妇嫩草av网站| 老汉色∧v一级毛片| 欧美激情高清一区二区三区| 国产日韩一区二区三区精品不卡| 巨乳人妻的诱惑在线观看| 国产黄频视频在线观看| av天堂久久9| 麻豆乱淫一区二区| 亚洲精华国产精华精| 丝袜美足系列| 欧美激情 高清一区二区三区| 久久久久久亚洲精品国产蜜桃av| 日本黄色视频三级网站网址 | 免费在线观看完整版高清| 可以免费在线观看a视频的电影网站| 亚洲精品国产区一区二| 亚洲精华国产精华精| 大片免费播放器 马上看| av视频免费观看在线观看| 国产亚洲午夜精品一区二区久久| 国产精品国产av在线观看| 操出白浆在线播放| 一级黄色大片毛片| 肉色欧美久久久久久久蜜桃| 国产成人免费无遮挡视频| 久久亚洲精品不卡| 一区二区av电影网| 亚洲熟女毛片儿| 中文字幕另类日韩欧美亚洲嫩草| 久久青草综合色| 人妻 亚洲 视频| 两人在一起打扑克的视频| 一区二区三区激情视频| 欧美精品人与动牲交sv欧美| 免费在线观看日本一区| av不卡在线播放| 狠狠婷婷综合久久久久久88av| 香蕉久久夜色| 亚洲av日韩在线播放| svipshipincom国产片| 久久久久久久国产电影| 免费看十八禁软件| 亚洲国产欧美网| 欧美变态另类bdsm刘玥| 亚洲第一欧美日韩一区二区三区 | 在线看a的网站| 欧美乱码精品一区二区三区| 视频区欧美日本亚洲| 搡老乐熟女国产| 丁香六月天网| 国产黄频视频在线观看| 欧美 亚洲 国产 日韩一| 亚洲精品国产精品久久久不卡| 亚洲精品国产色婷婷电影| 亚洲第一av免费看| 一进一出抽搐动态| 成人特级黄色片久久久久久久 | 国产在线一区二区三区精| 亚洲色图综合在线观看| 国产在视频线精品| 91成人精品电影| 国产成人精品久久二区二区免费| 超色免费av| 精品久久久精品久久久| 亚洲成av片中文字幕在线观看| 亚洲精品国产一区二区精华液| 一个人免费看片子| 日韩中文字幕视频在线看片| 国产欧美日韩一区二区三区在线| 日本一区二区免费在线视频| 自线自在国产av| 午夜精品久久久久久毛片777| 一级a爱视频在线免费观看| 亚洲av欧美aⅴ国产| 手机成人av网站| 69av精品久久久久久 | 18禁黄网站禁片午夜丰满| 每晚都被弄得嗷嗷叫到高潮| 国产欧美日韩一区二区精品| 成人av一区二区三区在线看| 视频区图区小说| 大型黄色视频在线免费观看| 黄色片一级片一级黄色片| 天天躁夜夜躁狠狠躁躁| 精品久久久久久电影网| 精品亚洲成国产av| 久久久久视频综合| 国产又色又爽无遮挡免费看| 成人18禁在线播放| 久久精品国产亚洲av香蕉五月 | 老熟妇乱子伦视频在线观看| 啪啪无遮挡十八禁网站| 九色亚洲精品在线播放| 高清毛片免费观看视频网站 | 人人澡人人妻人| 亚洲成av片中文字幕在线观看| 黄色a级毛片大全视频| 天堂动漫精品| 一级毛片精品| 久久中文字幕人妻熟女| 99国产精品一区二区三区| 亚洲精品中文字幕一二三四区 | 后天国语完整版免费观看| 蜜桃国产av成人99| 别揉我奶头~嗯~啊~动态视频| 国产高清视频在线播放一区| 国产精品美女特级片免费视频播放器 | 一进一出好大好爽视频| 9色porny在线观看| 99在线人妻在线中文字幕 | 韩国精品一区二区三区| netflix在线观看网站| videosex国产| 午夜福利影视在线免费观看| 在线看a的网站| 黄色视频在线播放观看不卡| av天堂在线播放| 菩萨蛮人人尽说江南好唐韦庄| 亚洲伊人色综图| 最新在线观看一区二区三区| 国产日韩欧美在线精品| 午夜成年电影在线免费观看| 1024香蕉在线观看| 脱女人内裤的视频| 久久国产精品大桥未久av| 91九色精品人成在线观看| 国产亚洲精品久久久久5区| 欧美精品啪啪一区二区三区| 母亲3免费完整高清在线观看| 国产精品免费视频内射| 亚洲熟女精品中文字幕| 韩国精品一区二区三区| 精品欧美一区二区三区在线| 国产主播在线观看一区二区| 最新美女视频免费是黄的| 高清视频免费观看一区二区| 中文字幕高清在线视频| 日本a在线网址| 国产精品99久久99久久久不卡| 亚洲人成电影观看| 纯流量卡能插随身wifi吗| 亚洲伊人久久精品综合| 在线观看一区二区三区激情| 免费av中文字幕在线| 国产精品免费一区二区三区在线 | 亚洲中文日韩欧美视频| 国内毛片毛片毛片毛片毛片| 久久精品亚洲熟妇少妇任你| av超薄肉色丝袜交足视频| 动漫黄色视频在线观看| 精品国产一区二区久久| 国产一区有黄有色的免费视频| 多毛熟女@视频| 亚洲,欧美精品.| 脱女人内裤的视频| 叶爱在线成人免费视频播放| 成年人午夜在线观看视频| 中国美女看黄片| 国产视频一区二区在线看| 色精品久久人妻99蜜桃| 激情视频va一区二区三区| 久9热在线精品视频| 一级黄色大片毛片| 人成视频在线观看免费观看| 9热在线视频观看99| 色精品久久人妻99蜜桃| 手机成人av网站| 亚洲欧美一区二区三区黑人| 黄色视频在线播放观看不卡| 欧美乱妇无乱码| 久久人人97超碰香蕉20202| 成人国产av品久久久| 热99re8久久精品国产| 91成人精品电影| 国产亚洲午夜精品一区二区久久| 国产伦理片在线播放av一区| 真人做人爱边吃奶动态| 久久精品国产亚洲av高清一级| 国产精品1区2区在线观看. | 黄片大片在线免费观看| 久久久久久久大尺度免费视频| 性少妇av在线| avwww免费| 91老司机精品| 99国产综合亚洲精品| 下体分泌物呈黄色| 欧美成狂野欧美在线观看| 久久精品国产a三级三级三级| 欧美成人免费av一区二区三区 | 性高湖久久久久久久久免费观看| 精品久久久精品久久久| 亚洲精品自拍成人| 成人av一区二区三区在线看| 欧美黄色淫秽网站| av视频免费观看在线观看| 国产99久久九九免费精品| 叶爱在线成人免费视频播放| 亚洲专区国产一区二区| 99精品久久久久人妻精品| 欧美精品高潮呻吟av久久| 国产片内射在线| 亚洲精品国产精品久久久不卡| 超碰成人久久| 久久精品成人免费网站| 日本黄色视频三级网站网址 | 午夜福利视频精品| 大陆偷拍与自拍| 欧美日韩av久久| 欧美日韩视频精品一区| 女人高潮潮喷娇喘18禁视频| 法律面前人人平等表现在哪些方面| 国产伦理片在线播放av一区| 欧美日韩国产mv在线观看视频| 天天操日日干夜夜撸| 热re99久久精品国产66热6| 亚洲人成伊人成综合网2020| 午夜福利乱码中文字幕| 中文字幕av电影在线播放| 性色av乱码一区二区三区2| 精品少妇一区二区三区视频日本电影| av福利片在线| 国产又爽黄色视频| 人人澡人人妻人| 少妇裸体淫交视频免费看高清 | 国产精品一区二区在线不卡| 欧美日韩中文字幕国产精品一区二区三区 | 99国产极品粉嫩在线观看| 久久久国产一区二区| 日本黄色日本黄色录像| 在线十欧美十亚洲十日本专区| 香蕉国产在线看| 国产在线观看jvid| 精品一区二区三区av网在线观看 | 精品国产乱子伦一区二区三区| 俄罗斯特黄特色一大片| 日韩一区二区三区影片| 精品久久久久久电影网| 一进一出好大好爽视频| 精品国产国语对白av| 高清黄色对白视频在线免费看| 精品国产超薄肉色丝袜足j| av片东京热男人的天堂| 国产91精品成人一区二区三区 | 久久精品成人免费网站| av天堂在线播放| 制服人妻中文乱码| 黑丝袜美女国产一区| 中文亚洲av片在线观看爽 | 美国免费a级毛片| 欧美亚洲 丝袜 人妻 在线| videos熟女内射| 国产精品一区二区在线不卡| 欧美日韩亚洲综合一区二区三区_| 精品少妇一区二区三区视频日本电影| 亚洲伊人久久精品综合| 亚洲,欧美精品.| av不卡在线播放| 亚洲熟女精品中文字幕| 午夜福利欧美成人| 日日摸夜夜添夜夜添小说| 人成视频在线观看免费观看| 好男人电影高清在线观看| av一本久久久久| 精品国产一区二区三区久久久樱花| 欧美成人午夜精品| 亚洲精品乱久久久久久| 天堂动漫精品| 国产精品国产高清国产av | 精品高清国产在线一区| 一边摸一边抽搐一进一小说 | 国产亚洲精品久久久久5区| 婷婷成人精品国产| 欧美日韩亚洲国产一区二区在线观看 | 精品国产一区二区三区四区第35| 看免费av毛片| 亚洲精华国产精华精| 一二三四在线观看免费中文在| 日韩精品免费视频一区二区三区| 狠狠婷婷综合久久久久久88av| 精品欧美一区二区三区在线| 热re99久久国产66热| 老熟妇仑乱视频hdxx|