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

    基于離散元方法的豎直群樁冰荷載分析

    2017-02-27 00:04:37李輝輝張寶森翟必垚季順迎
    南水北調(diào)與水利科技 2017年1期

    李輝輝+張寶森+翟必垚+季順迎

    摘要:在高緯度海域或高寒地區(qū)的江河中,碎冰對海洋平臺樁腿和橋墩等直立群樁結(jié)構(gòu)產(chǎn)生很大的沖擊作用,并在一定條件下造成不同程度的破壞。依據(jù)碎冰在自然條件下的離散分布特性,采用Voronoi切割算法構(gòu)造隨機(jī)分布和具有非規(guī)則幾何形態(tài)的碎冰初始狀態(tài),在此基礎(chǔ)上采用若干個具有黏結(jié)-破碎功能的球形顆粒構(gòu)造具有一定厚度、尺寸和形狀等物理參數(shù)的河冰離散單元模型,將豎直群樁簡化為剛體結(jié)構(gòu),從而對碎冰的運(yùn)移及其對群樁的沖擊作用進(jìn)行數(shù)值模擬。在此基礎(chǔ)上可確定碎冰撞擊下的群樁冰荷載,并分析冰塊尺寸、冰速和冰厚對樁結(jié)構(gòu)所受冰荷載的影響。

    關(guān)鍵詞:冰荷載;離散元方法;豎直群樁

    中圖分類號:TU473文獻(xiàn)標(biāo)志碼:A文章編號:1672-1683(2017)01-0126-06

    Abstract:In the seas and rivers in high-latitude regions,the broken ice floes can generate obvious ice load on offshore platform piles and bridge piers,and even cause damage of the vertical piles.Based on the discrete distribution characteristics of ice floes under natural conditions,the Voronoi cutting algorithm was used to construct the initial state of randomly-distributed broken ice floes with irregular geometric shapes.The discrete element method was adopted to simulate the drifting process of ice floes and the interaction between ice floes and vertical piles.The ice floe model with certain thickness,size,shape and other physical parameters was constructed with bonded spheres with consideration to their failure effect.The vertical piles were simplified as rigid structures.A numerical simulation was conducted to determine the ice load on the pile group under the impact of ice floes,and analyze the effects of ice size,ice velocity and ice thickness on the ice load of pile structures.

    Key words: ice load;discrete element method;vertical pile group

    對冰的認(rèn)識和研究最初始于北方內(nèi)河航道交通的需要。黃河以北的地區(qū),包括新疆、內(nèi)蒙和東北三省,秋末冬初,氣溫逐漸降低,太陽輻射減弱,江河出現(xiàn)水凍結(jié)成冰的現(xiàn)象,浮冰隨流水順流而下形成秋季流冰[1]。春季由于氣候變暖、太陽輻射加強(qiáng),封凍期形成的大面積的冰蓋融化破碎形成尺寸大小不等的冰塊,順流而下,形成春季流冰[2]。流冰在經(jīng)過水工建筑物時,會造成對水工建筑物有嚴(yán)重作用的流冰荷載。

    冰荷載是寒區(qū)內(nèi)河水工建筑物的主要荷載之一。水工建筑物因冰凌的靜力動力作用遭到破壞的實例時有發(fā)生[3]。20世紀(jì)50年代官廳水庫中的橋梁墩臺被冰力剪斷;20世紀(jì)60年代渤海灣中海洋平臺的烽火臺被流冰撞倒。國外,20世紀(jì)60年代美國阿拉斯加庫克灣兩座采油平臺在劇烈的冰振下倒塌;日本雅內(nèi)港的聲向崎燈標(biāo)被冰推倒;20世紀(jì)70年代芬蘭波茲尼亞灣的鋼制燈塔在流冰撞擊下倒塌;20世紀(jì)90年代美國緬因州的一座鋼筋混凝土大橋被流冰撞毀[4]。隨著社會的發(fā)展,流冰引發(fā)的災(zāi)害現(xiàn)象也越來越得到人們的重視[5-7]。

    目前,國內(nèi)外對冰荷載的研究包括試驗和理論分析兩種方法。試驗方法主要是在原型結(jié)構(gòu)上進(jìn)行現(xiàn)場測量和室內(nèi)實驗室試驗,如加拿大Beaufort海域Molikpaq平臺[8]、加拿大聯(lián)邦大橋橋墩以及Baltic海Kemi-I燈塔上的冰荷載測量[9-10],國內(nèi)主要集中在渤海的海洋平臺的測量工作[11-13]。理論分析主要有理論模型分析和數(shù)值模型分析兩種方法,理論模型有Matlock首次提出的Matlock模型,Karna引入失效區(qū)域?qū)ζ溥M(jìn)行了改進(jìn)[14],Withalm等其他學(xué)者也對其進(jìn)行不斷的發(fā)展與應(yīng)用[15-16]。數(shù)值模型主要包括有限元方法和離散元方法,其更加關(guān)注冰作為復(fù)雜的工程材料這一特點(diǎn),更多的考慮到冰的力學(xué)性質(zhì),因而對理論分析法而言有了較大改進(jìn)。

    近年來,離散元方法得到廣泛的關(guān)注,尤其在對冰的研究領(lǐng)域內(nèi)發(fā)揮著重要作用。該方法在細(xì)觀上能夠合理地描述冰材料的結(jié)構(gòu)性質(zhì),在宏觀上可以準(zhǔn)確地模擬出冰與結(jié)構(gòu)作用中的破碎過程[17]。如Hopkins采用擴(kuò)展圓盤[18-19]來模擬近似圓形的浮冰,并將其應(yīng)用于對Alaska北岸浮冰的運(yùn)動、斷裂情況的分析,Hansen采用二維圓盤離散單元模型對錨系船舶在矩形浮冰中的航行進(jìn)行了分析[20],另外還有很多學(xué)者利用離散元方法對河冰或者海冰進(jìn)行了不斷深入的研究[21-22]。

    為此,本文采用離散元方法,將冰離散為若干個具有黏結(jié)-破碎功能的球體顆粒單元,將豎直群樁簡化為剛體結(jié)構(gòu),對冰運(yùn)移以及冰對群樁的作用過程進(jìn)行數(shù)值模擬。分析不同工況(冰厚、尺寸、冰速)對樁結(jié)構(gòu)所受冰荷載的影響。

    1 冰的離散單元模型

    在離散元模擬中,顆粒單元之間視為在法向上由彈簧與阻尼器連接,在切向上由彈簧、阻尼器和滑動摩擦器連接,見圖1(a)。法向彈性力為Fn=Knxn-Cnn,切向彈性力為Fs=min(Ksxs-Csx,μFn)。式中Kn和Ks分別是法向和切向剛度系數(shù),Cn和Cs是法向和切向阻尼系數(shù),x和分別為顆粒的變形和變形率。法向阻尼Cn=ζ2MKn,其中ζn為無量綱黏滯系數(shù),ζn=-lneπn+ln2e,e為回彈系數(shù)。切向與法向的剛度、阻尼系數(shù)有如下關(guān)系:Ks=αKn,Cs=βCn。這里取α=0.5,β=0.5。

    考慮單元間的凍結(jié)作用,建立了顆粒單元的黏結(jié)-破碎功能。兩個黏結(jié)單元之間采用平行黏結(jié)模型,如圖1(b)所示,即在兩個顆粒單元之間設(shè)定彈性黏結(jié)圓盤,傳遞軸向力、剪力、彎矩和扭矩,力和力矩都用法向和切向分量來表示,i=ni+si,i=ni+si。式中,ni、ni、si、si分別為力與力矩的法向分量和切向分量。

    在顆粒相互作用時,相對位置和相對轉(zhuǎn)角的變化引起黏彈性力的變化,并通過判斷最大軸向力和剪切力分別與給定的法向黏結(jié)強(qiáng)度和切向黏結(jié)強(qiáng)度做比較,判斷顆粒間的黏結(jié)是否發(fā)生破壞,黏結(jié)作用是否消失,宏觀的表現(xiàn)是冰塊是否出現(xiàn)裂紋,是否發(fā)生破碎。

    2 豎直樁冰荷載的離散元分析

    本文計算模擬的矩形水道見圖2,長500 m寬100 m,左側(cè)冰區(qū)長300 m寬100 m。群樁位于冰區(qū)前方100 m處,共有5個直立樁,樁的直徑為2 m,樁間距20 m。浮冰是采用Voronoi切割算法構(gòu)造的,隨機(jī)分布且形狀非規(guī)則,由單層球體單元黏結(jié),單元直徑為0.3 m,豎直樁采用圓柱體單元。整個冰區(qū)的浮冰密集度為80%,厚度為0.3 m,平均尺寸是25 m2,水面高度是20 m,水流速度為1 m/s。整個冰區(qū)有1 600個浮冰塊,共由172 695個球體顆粒組成。具體離散元參數(shù)見表1。

    碎冰在水流速度的拖曳力作用下基本保持勻速運(yùn)動,到達(dá)群樁處時撞擊樁腿,有冰堆積,部分冰塊出現(xiàn)破碎,在直立樁前堆積的冰塊能迅速繞過樁腿繼續(xù)向前運(yùn)動,整個過程如圖3所示。初始時刻碎冰在直水道的分布情況見圖3(a),不同顏色代表不同冰塊。浮冰會發(fā)生三個相互垂直方向上的平動和轉(zhuǎn)動,其中規(guī)定浮冰前進(jìn)的方向為x方向,垂直于水面向上為z方向,x、y、z三個方向符合右手定則。經(jīng)過計算可以提取得到5個豎直樁的三個方向的冰荷載加以研究。

    統(tǒng)計剩下4個豎直樁的x方向的冰力時程曲線如圖5所示,5個樁腿冰力峰值分別是860.77 kN、919.98 kN、1 208.59 kN、1 201.78 kN、1 203.47 kN,相應(yīng)的均值為85.40 kN、69.22 kN、91.18 kN、94.94 kN、86.14 kN。因為水道邊界形狀規(guī)則,碎冰尺寸相近,碎冰運(yùn)移時冰速保持恒定,所以5個豎直樁所受冰力具有很高的相似性,冰力峰值和均值也都很相近。因此本文在分析冰荷載時,只取水道中間的1號樁腿所受的冰荷載作為研究對象,這會和研究其他樁腿的冰力得到相同的結(jié)論。

    3 豎直樁上冰荷載的影響因素分析

    3.1 冰速對冰荷載的影響

    冰速對冰荷載的影響顯而易見,不論是動冰能量理論和動冰強(qiáng)度極限理論,都將冰速作為影響冰力的關(guān)鍵因素加以考慮。如工程中流動冰塊對豎直樁沖擊荷裁P=VhARa,V為冰塊流速,h為冰厚,A為冰蓋面積,Ra為壓縮強(qiáng)度,冰速V作為單獨(dú)的一項加以研究[4]。而且冰的強(qiáng)度值與加載速率緊密相關(guān),高應(yīng)變率下表現(xiàn)為脆性材料,低應(yīng)變率下表現(xiàn)為彈性材料,因而冰速又間接影響冰的壓縮強(qiáng)度與彎曲強(qiáng)度。為研究冰速與冰荷載的關(guān)系,本文將保持冰厚為0.4 m,尺寸為25 m2,冰速分別設(shè)置成1 m/s、2 m/s、3 m/s、4 m/s、5 m/s、6 m/s對比,得到冰荷載時程。

    計算提取不同速度下冰荷載最大值和相應(yīng)的均值,由此做出曲線見圖6,對比可以看出冰速對于冰荷載的影響很大,冰力隨著冰速增加大體呈線性增長的趨勢。根據(jù)動冰能量理論推導(dǎo)的流冰對河道上的橋墩作用冰載荷計算公式:P=0.4VhmsΩR,式中,V是流冰速度;h是冰厚,ms是冰塊質(zhì)量;Ω冰塊面積;R是局部擠壓強(qiáng)度。在工程應(yīng)用中,將冰速與冰荷載大小視為成線性相關(guān),數(shù)值模擬結(jié)果與工程規(guī)范很好的對應(yīng)。這主要是因為冰速的增加會導(dǎo)致浮冰與樁腿碰撞時的動量增大,產(chǎn)生的冰力也隨之增大。模擬結(jié)果很好的驗證了根據(jù)動冰能量理論所建立規(guī)范的合理性。

    3.2 冰塊尺寸對冰荷載的影響

    浮冰尺寸作為流冰動能的重要影響要素,在動冰與樁撞擊過程中也會影響冰荷載的大小。本文為研究尺寸與冰荷載的關(guān)系,將冰速設(shè)置為1 m/s,冰厚為0.4 m,生成冰塊平均面積分別為15 m2、20 m2、25 m2、30 m2、35 m2、40 m2,尺寸大小見圖7,其他計算參數(shù)同前。

    計算模擬了浮冰250 s運(yùn)移過程,記錄了其x方向的冰力,計算提取得到不同尺寸下冰荷載最大值和相應(yīng)的均值,并由此做出曲線見圖8。結(jié)果表明,冰塊尺寸對冰荷載有一定影響,隨著冰塊尺寸的增加,樁所受冰荷載隨之增加。很多工程應(yīng)用中都只將冰與建筑物接觸的水平投影寬度作為影響因素考慮,冰塊尺寸間接影響碎冰碰撞樁腿的接觸面積,從而影響冰力大??;另一方面浮冰尺寸的增加會導(dǎo)致整塊浮冰質(zhì)量增加,動量也隨之增大,產(chǎn)生的冰荷載也隨之增大。由冰荷載均值可以看出,在一定尺寸范圍內(nèi),冰荷載遵循動冰能量理論,即冰荷載受碎冰動能決定,隨著尺寸的增大,冰塊的動能增大,冰荷載均值也隨之增大;尺寸超過某個臨界值之后,冰荷載遵循動冰強(qiáng)度極限理論,即統(tǒng)計的冰荷載是在碎冰碰撞樁腿發(fā)生破壞的時刻產(chǎn)生,由碎冰的破壞強(qiáng)度決定,尺寸的變化并不影響其破壞強(qiáng)度,因此在圖8(b)后半段的冰荷載均值受尺寸的影響很小。

    3.3 冰厚對冰荷載的影響

    冰厚與冰速一樣,是影響冰力的重要因素。為了研究冰厚與冰荷載的關(guān)系,將保持冰塊尺寸為25 m2,冰速為4 m/s,冰厚分別設(shè)置成0.3 m、0.4 m、0.5 m、0.6 m做了四個算例對比,得到的冰荷載時程,計算提取得到每個尺寸下冰荷載最大峰值和相應(yīng)的均值,由此做出曲線見圖9。可以看出,冰厚對冰荷載影響較大,隨著冰塊尺寸的增大,樁腿所受垂向冰荷載也隨之增大,之間近似成線性關(guān)系。主要原因是冰厚增加,冰塊質(zhì)量增大,從而引起碎冰對樁腿的沖擊動量的增加;另一方面,隨著冰厚增加,模擬冰塊的顆粒粒徑相應(yīng)增加,因而顆粒之間的黏結(jié)力也增大,宏觀上即表現(xiàn)為冰塊破壞強(qiáng)度的增大,冰塊破壞前給樁腿造成的冰荷載的增大。

    4 結(jié)論

    本文采用了離散元方法對冰的運(yùn)移以及冰對群樁的作用過程進(jìn)行數(shù)值模擬,討論了不同影響因素下,冰荷載的變化規(guī)律。以上研究表明,本文建立的離散元法能夠合理計算流凌在水域中的運(yùn)移過程,可對冰與樁的動力作用過程進(jìn)行模擬,能得到很好的冰力時程曲線。經(jīng)過分析可以發(fā)現(xiàn)冰厚、冰塊尺寸和冰速都是影響冰力大小的重要因素。

    由于影響冰荷載的因素諸多,加之工程規(guī)范中的計算公式各系數(shù)項的意義及取值原則,在實際工程中不可照搬,應(yīng)該根據(jù)具體情況設(shè)置。用離散元法對冰的模擬更加接近實際情況,根據(jù)實際工況設(shè)置參數(shù)進(jìn)行數(shù)值模擬的計算結(jié)果可以為規(guī)范計算提供很好的參考,也可為實際工程問題提供很好的參照。

    參考文獻(xiàn)(References):

    [1] 陸欽年,段忠東,歐進(jìn)萍,等.黑龍江省的凌汛災(zāi)害及其區(qū)域劃分[J].自然災(zāi)害學(xué)報,2000(3):54-58.(LU Qin-nian,DUAN Zhong-dong,OU Jin-ping,et al.Ice flood and its zonation in Heilongjiang Province [J].Journal of Natural Disasters,2000(3):54-58.(in Chinese))DOI:10.3969/j.issn.1004-4574.2000.03.008

    [2] 于天來,張礫元,雷俊卿,等.春季流冰對橋墩產(chǎn)生動冰壓力時冰抗壓強(qiáng)度取值探討[J].中外公路,2010,30(3):168-172.(YU Tian-lai,ZHANG Li-yuan,LEI Jun-qing,et al.Value of ice compressive strength when spring ice produces dynamic ice pressure on the pier [J].Journal of China & Foreign Highway,2010,30(3):168-172.(in Chinese)) DOI:10.14048/j.issn.1671-2579.2010.03.028

    [3] 蘇盛奎.冰荷載對水工建筑物的作用[J].水運(yùn)工程,1980(5).(SU Sheng-kui.The effect of ice load on hydraulic structures [J].Port & Waterway Engineering,1980(5).(in Chinese)) DOI:10.16233/j.cnki.issn1002-4972.1980.05.008

    [4] 蔡之瑞,孫柏濤,郭世榮,等.冰荷載的實驗研究與計算方法[J].地震工程與工程振動,1997(4):49-56.(CAI Zhi-rui,SUN Bai-tao,GUO Shi-rong,et al.Experimental research on ice load and its calculation method [J].Earthquake Engineering and Engineering Dynamics,1997(4):49-56.(in Chinese)) DOI:10.13197/j.eeev.1997.04.007

    [5] 宋波,牛立超,黃帥,等.強(qiáng)震作用下冰水域橋墩動力響應(yīng)研究分析[J].巖土力學(xué),2014(S1):285-291.(SONG Bo,NIU Li-chao,HUANG Shuai,et al.Research on dynamic response of pier column structures in ice-plagued area under highly seismic action [J].Rock and Soil Mechanics,2014(S1):285-291.(in Chinese)) DOI:10.16285/j.rsm.2014.s1.041

    [6] 賈玲玲.波浪力對深水橋墩結(jié)構(gòu)在地震和冰荷載作用下的響應(yīng)影響分析[J].震災(zāi)防御技術(shù),2010,05(2):263-269.(JIA Ling-ling.Seismic and ice response analysis of bridge pier in deep water with the wave effect [J].Technology for Earthquake Disaster Prevention,2010,05(2):263-269.(in Chinese)) DOI:10.3969/j.issn.1673-5722.2010.02.015

    [7] 吳云鳳.黃河海勃灣水利樞紐施工期防凌設(shè)計[J].人民黃河,2013(11):90-92.(WU Yun-feng.Ice prevention design during construction of Haibowan water control project on the Yellow River [J].Yellow River,2013(11):90-92.(in Chinese))DOI:10.3969/j.issn.1000-1379.2013.11.032

    [8] Timco G W,Johnston M.Ice Loads on the Molikpaq in the Canadian Beaufort Sea [J].Cold Regions Science and Technology,2003,37:51-68.DOI:10.1016/S0165-232X(03)00035-1

    [9] Brown T G,Mttnen M.Comparison of Kemi-I and Confederation Bridge Cone Ice Load Measurement Results [J].Cold Regions Science and Technology,2009,55:3-13.DOI:10.1016/j.coldregions.2008.04.005

    [10] Brown T G,Tibbo J S,Tripathi D,et al.Extreme ice load events on the Confederation Bridge [J].Cold Regions Science and Technology,2010,60:1-14.DOI:10.1016/j.coldregions.2009.08.004

    [11] Qu Y,Yue Q,Bi X,et al.A random ice force model for narrow conical structures[J].Cold Regions Science and Technology,2006,45:148-157.DOI:10.1016/j.coldregions.2006.05.008

    [12] Yue Q,Bi X.Ice-induced jacket structure vibrations in Bohai Sea [J].Journal of Cold Regions Engineering,2000,14:81-92.DOI:10.1061/(ASCE)0887-381X(2000)14:2(81)

    [13] 岳前進(jìn),畢祥軍,于曉,等.錐體結(jié)構(gòu)的冰激振動與冰力函數(shù)[J].土木工程學(xué)報,2003,36:16-19.(YUE Qian-jin,BI Xiang-jun,YU Xiao,et al.Ice-induced vibration and ice force function of conical structure [J].China Civil Engineering Journal,2003,36:16-19.(in Chinese)) DOI:10.3321/j.issn:1000-131X.2003.02.004

    [14] Krn T,Kamesaki K,Tsukuda H.A numerical model for dynamic ice-structure interaction [J].Computers & Structures,1999,72:645-658.DOI:10.1016/S0045-7949(98)00337-X

    [15] Withalm M,Hoffmann N P.Simulation of full-scale ice-structure-interaction by an extended matlock-model [J].Cold Regions Science and Technology,2010,60:130-136.DOI:10.1016/j.coldregions.2009.09.006

    [16] Yue Q,Guo F,Krn T.Dynamic ice forces of slender vertical structures due to ice crushing [J].Cold Regions Science & Technology,2009,56(2-3):77-83.DOI:10.1016/j.coldregions.2008.11.008

    [17] 季順迎,李紫麟,李春花,等.碎冰區(qū)海冰與船舶結(jié)構(gòu)相互作用的離散元分析[J].應(yīng)用力學(xué)學(xué)報,2013(4):520-526.(JI Shun-ying,LI Zi-lin,LI Chun-hua,et al.Discrete element modeling of ice loads on ship hulls in broken ice fields [J].Chinese Journal of Applied Mechanics,2013(4):520-526.(in Chinese)) DOI:10.11776/cjam.30.04.D032

    [18] Hopkins M A.Discrete element modeling with dilated particles [J].Engineering Computations,2004,21:422-430.DOI:10.1108/02644400410519866

    [19] Hopkins M A,Shen H H.Simulation of pancake-ice dynamics in a wave field [J].Annals of glaciology,2001,33:355-360.DOI:10.1016/j.coldregions.2012.02.003

    [20] Hansen E H,Lset S.Modelling floating offshore units moored in broken ice:model description [J].Cold Regions Science and Technology,1999,29:97-106.DOI:10.1016/S0165-232X(99)00023,3

    [21] Lau M,Lawrence KP,Rothenburg L.Discrete element analysis of ice loads on ships and structures [J].Ships and Offshore Structures,2011,6(3):211-221.DOI:10.1080/17445302.2010.544086

    [22] Lu W,Lubbad R,Lset S.Simulating ice-sloping structure interactions with the cohesive element method [J].Journal of Offshore Mechanics and Arctic Engineering,2014,136:031501.

    日韩欧美 国产精品| 男女床上黄色一级片免费看| 久久婷婷成人综合色麻豆| 欧美乱码精品一区二区三区| 老司机福利观看| 中文字幕人妻熟女乱码| 久久中文字幕一级| 亚洲自拍偷在线| 免费一级毛片在线播放高清视频| 色综合欧美亚洲国产小说| 国产欧美日韩一区二区精品| 19禁男女啪啪无遮挡网站| 国产又爽黄色视频| av电影中文网址| 国产精品美女特级片免费视频播放器 | 两个人免费观看高清视频| 免费无遮挡裸体视频| 丰满人妻熟妇乱又伦精品不卡| 精品熟女少妇八av免费久了| 少妇裸体淫交视频免费看高清 | av在线天堂中文字幕| 亚洲中文字幕日韩| 亚洲av成人一区二区三| 国产精品永久免费网站| 热99re8久久精品国产| 午夜免费激情av| 18禁美女被吸乳视频| 久久99热这里只有精品18| 国产不卡一卡二| 12—13女人毛片做爰片一| 久久久久久人人人人人| 大型黄色视频在线免费观看| 日韩欧美 国产精品| 日韩成人在线观看一区二区三区| 自线自在国产av| 老汉色∧v一级毛片| 成年版毛片免费区| 亚洲欧美精品综合一区二区三区| 妹子高潮喷水视频| 身体一侧抽搐| 国产黄片美女视频| 欧美三级亚洲精品| 精品不卡国产一区二区三区| 曰老女人黄片| 免费无遮挡裸体视频| 国产精品国产高清国产av| 亚洲av五月六月丁香网| 成人国产综合亚洲| 成人永久免费在线观看视频| 91大片在线观看| 国产高清videossex| 特大巨黑吊av在线直播 | 国产精品1区2区在线观看.| 手机成人av网站| 99国产精品一区二区三区| 9191精品国产免费久久| 日本撒尿小便嘘嘘汇集6| 女人被狂操c到高潮| x7x7x7水蜜桃| 日韩高清综合在线| 91av网站免费观看| 亚洲黑人精品在线| 丝袜人妻中文字幕| av视频在线观看入口| 禁无遮挡网站| 日韩大码丰满熟妇| 十八禁人妻一区二区| 大香蕉久久成人网| 一本久久中文字幕| 制服诱惑二区| 亚洲成av人片免费观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国内毛片毛片毛片毛片毛片| 久久精品影院6| 国产成年人精品一区二区| 久久精品国产亚洲av香蕉五月| 亚洲国产看品久久| 日韩成人在线观看一区二区三区| 精华霜和精华液先用哪个| 精品久久久久久,| 亚洲中文字幕日韩| 精品福利观看| 美女国产高潮福利片在线看| 热re99久久国产66热| 国产aⅴ精品一区二区三区波| 国产亚洲欧美精品永久| 午夜精品久久久久久毛片777| 亚洲aⅴ乱码一区二区在线播放 | 波多野结衣高清无吗| 在线观看66精品国产| 777久久人妻少妇嫩草av网站| 韩国av一区二区三区四区| 看片在线看免费视频| 亚洲在线自拍视频| 精品久久久久久久久久免费视频| 国产在线精品亚洲第一网站| 欧美一级a爱片免费观看看 | 成年女人毛片免费观看观看9| 国产真实乱freesex| 真人一进一出gif抽搐免费| 一边摸一边抽搐一进一小说| 亚洲成av人片免费观看| 高清在线国产一区| 精品高清国产在线一区| 中文资源天堂在线| 精品午夜福利视频在线观看一区| 久久国产精品人妻蜜桃| www.精华液| 脱女人内裤的视频| 黄色毛片三级朝国网站| 看免费av毛片| 91九色精品人成在线观看| 中文亚洲av片在线观看爽| av天堂在线播放| 国产男靠女视频免费网站| 国产精品久久电影中文字幕| 国产精品九九99| 久久久久久久午夜电影| 久久国产精品男人的天堂亚洲| 好男人电影高清在线观看| 美女国产高潮福利片在线看| 成人特级黄色片久久久久久久| 国产精品综合久久久久久久免费| 黄色 视频免费看| 亚洲专区国产一区二区| 男人舔女人的私密视频| 午夜免费观看网址| 成人免费观看视频高清| 国产亚洲精品久久久久5区| 一夜夜www| 欧美成人一区二区免费高清观看 | 女人爽到高潮嗷嗷叫在线视频| 国产黄片美女视频| 又紧又爽又黄一区二区| 美女国产高潮福利片在线看| 韩国精品一区二区三区| 亚洲美女黄片视频| 一进一出好大好爽视频| 长腿黑丝高跟| 亚洲狠狠婷婷综合久久图片| 欧美成人一区二区免费高清观看 | 久久国产精品男人的天堂亚洲| 99热只有精品国产| 亚洲一码二码三码区别大吗| 麻豆成人午夜福利视频| 婷婷丁香在线五月| 看免费av毛片| 亚洲国产精品合色在线| 日韩免费av在线播放| 精品国产一区二区三区四区第35| 琪琪午夜伦伦电影理论片6080| 色播亚洲综合网| tocl精华| 亚洲成国产人片在线观看| 天堂影院成人在线观看| 亚洲国产日韩欧美精品在线观看 | 亚洲专区中文字幕在线| 久久亚洲精品不卡| 国产一区二区激情短视频| av超薄肉色丝袜交足视频| 国产麻豆成人av免费视频| 精品高清国产在线一区| 成人特级黄色片久久久久久久| 亚洲国产看品久久| 欧美黄色片欧美黄色片| 亚洲av第一区精品v没综合| 两性夫妻黄色片| 成年女人毛片免费观看观看9| 我的亚洲天堂| 色婷婷久久久亚洲欧美| 欧美三级亚洲精品| 婷婷六月久久综合丁香| 国产精品,欧美在线| 久久香蕉国产精品| 亚洲aⅴ乱码一区二区在线播放 | 亚洲黑人精品在线| tocl精华| 亚洲成国产人片在线观看| 日韩欧美国产在线观看| 国产成人av激情在线播放| 国产熟女xx| 一级a爱视频在线免费观看| 国产一级毛片七仙女欲春2 | 亚洲色图 男人天堂 中文字幕| 日韩一卡2卡3卡4卡2021年| 精品欧美一区二区三区在线| 男男h啪啪无遮挡| 韩国精品一区二区三区| 黄色视频不卡| 色播亚洲综合网| 亚洲avbb在线观看| 99国产精品一区二区三区| 国产精品久久久av美女十八| 欧美激情高清一区二区三区| 麻豆国产av国片精品| 日韩有码中文字幕| 女人被狂操c到高潮| 欧美绝顶高潮抽搐喷水| 中文亚洲av片在线观看爽| 国内揄拍国产精品人妻在线 | 又大又爽又粗| 夜夜夜夜夜久久久久| 午夜日韩欧美国产| 久久久久久久久中文| 精品福利观看| ponron亚洲| 一a级毛片在线观看| 久久精品影院6| a级毛片在线看网站| 大型av网站在线播放| 99精品久久久久人妻精品| 国产区一区二久久| 中文在线观看免费www的网站 | 丰满人妻熟妇乱又伦精品不卡| 真人做人爱边吃奶动态| 少妇熟女aⅴ在线视频| 国产午夜福利久久久久久| 欧美日韩瑟瑟在线播放| 一区二区日韩欧美中文字幕| 午夜免费激情av| 亚洲 欧美一区二区三区| 中文字幕精品免费在线观看视频| 日韩精品免费视频一区二区三区| 18禁美女被吸乳视频| 婷婷丁香在线五月| 黑人巨大精品欧美一区二区mp4| 久久国产乱子伦精品免费另类| 满18在线观看网站| 两人在一起打扑克的视频| 老司机深夜福利视频在线观看| 日韩精品免费视频一区二区三区| 天堂√8在线中文| 国产黄a三级三级三级人| 欧美黄色淫秽网站| 香蕉久久夜色| 日韩成人在线观看一区二区三区| 丝袜美腿诱惑在线| av天堂在线播放| 亚洲七黄色美女视频| 亚洲一区高清亚洲精品| 亚洲国产中文字幕在线视频| 日本免费a在线| 给我免费播放毛片高清在线观看| 欧美久久黑人一区二区| 日本黄色视频三级网站网址| 久久精品国产亚洲av香蕉五月| 亚洲精华国产精华精| 日韩大尺度精品在线看网址| 91字幕亚洲| 伊人久久大香线蕉亚洲五| 亚洲精品美女久久av网站| 国产一区二区激情短视频| tocl精华| 最新美女视频免费是黄的| 色哟哟哟哟哟哟| 免费高清在线观看日韩| 亚洲成国产人片在线观看| 亚洲人成网站高清观看| 亚洲色图av天堂| 久久国产乱子伦精品免费另类| 亚洲精品美女久久久久99蜜臀| 色精品久久人妻99蜜桃| 麻豆成人午夜福利视频| 啦啦啦韩国在线观看视频| 日本免费一区二区三区高清不卡| 90打野战视频偷拍视频| 国产欧美日韩一区二区精品| 亚洲 欧美一区二区三区| 精品一区二区三区视频在线观看免费| 天堂动漫精品| 99riav亚洲国产免费| 欧美中文综合在线视频| 无限看片的www在线观看| 亚洲精品在线观看二区| 一个人免费在线观看的高清视频| 精品熟女少妇八av免费久了| 最新在线观看一区二区三区| 国产精品一区二区免费欧美| 亚洲电影在线观看av| 在线永久观看黄色视频| 亚洲黑人精品在线| 人妻久久中文字幕网| 黑人操中国人逼视频| 伊人久久大香线蕉亚洲五| 日韩三级视频一区二区三区| 桃红色精品国产亚洲av| 最近最新中文字幕大全电影3 | 欧美性长视频在线观看| 中文亚洲av片在线观看爽| 日韩欧美国产一区二区入口| 中文字幕人妻丝袜一区二区| cao死你这个sao货| 中文字幕人妻熟女乱码| 女人被狂操c到高潮| 黄片播放在线免费| 欧美日韩精品网址| 久久精品人妻少妇| 天堂动漫精品| 午夜福利在线观看吧| 免费看日本二区| 麻豆一二三区av精品| 999精品在线视频| 日韩欧美一区二区三区在线观看| 免费观看人在逋| 久9热在线精品视频| 亚洲av电影不卡..在线观看| 搞女人的毛片| 国产成人欧美| 亚洲成人国产一区在线观看| 日日干狠狠操夜夜爽| 很黄的视频免费| 国产精品电影一区二区三区| 国产精品1区2区在线观看.| 国产av在哪里看| 久久精品影院6| 中亚洲国语对白在线视频| 日韩一卡2卡3卡4卡2021年| 精品一区二区三区av网在线观看| 又黄又爽又免费观看的视频| 欧美+亚洲+日韩+国产| 深夜精品福利| 日本 欧美在线| 亚洲欧美精品综合久久99| 久久久久久久久中文| 不卡一级毛片| 日韩欧美一区视频在线观看| 亚洲五月天丁香| 亚洲国产中文字幕在线视频| 亚洲av片天天在线观看| 在线观看免费视频日本深夜| 国产不卡一卡二| 欧美性猛交╳xxx乱大交人| 一区二区三区精品91| 少妇的丰满在线观看| 免费观看精品视频网站| 不卡一级毛片| 久久天堂一区二区三区四区| 精品久久久久久久末码| 久久人人精品亚洲av| 亚洲久久久国产精品| 香蕉丝袜av| 久9热在线精品视频| 久久久久精品国产欧美久久久| 日本熟妇午夜| 熟女电影av网| 欧美黑人精品巨大| 精品免费久久久久久久清纯| 两人在一起打扑克的视频| 欧美最黄视频在线播放免费| 老司机靠b影院| 免费高清视频大片| 欧美色视频一区免费| 亚洲欧美激情综合另类| 人人妻人人看人人澡| 少妇粗大呻吟视频| 亚洲久久久国产精品| 久久午夜亚洲精品久久| 99久久综合精品五月天人人| 无限看片的www在线观看| 欧美日本亚洲视频在线播放| 久久久久久大精品| 国产真实乱freesex| 高潮久久久久久久久久久不卡| 伊人久久大香线蕉亚洲五| 黄片小视频在线播放| 一级片免费观看大全| 中文在线观看免费www的网站 | 最近最新中文字幕大全电影3 | 97碰自拍视频| 亚洲精品国产精品久久久不卡| 欧美成狂野欧美在线观看| 级片在线观看| 午夜激情av网站| 午夜视频精品福利| 美女大奶头视频| 久久久国产精品麻豆| 禁无遮挡网站| 精品第一国产精品| 欧美激情久久久久久爽电影| 99国产极品粉嫩在线观看| 19禁男女啪啪无遮挡网站| 久久午夜亚洲精品久久| 精品久久久久久久久久免费视频| 夜夜看夜夜爽夜夜摸| 美女高潮喷水抽搐中文字幕| 久99久视频精品免费| 中文字幕高清在线视频| 国产精品综合久久久久久久免费| 天天添夜夜摸| 亚洲男人的天堂狠狠| 精品不卡国产一区二区三区| 久久久久久久久久黄片| bbb黄色大片| 又黄又粗又硬又大视频| 午夜福利一区二区在线看| 又大又爽又粗| 啦啦啦韩国在线观看视频| 精品人妻1区二区| 国产黄色小视频在线观看| 亚洲欧美精品综合久久99| 亚洲精品一卡2卡三卡4卡5卡| 国产成年人精品一区二区| 人人妻人人澡人人看| 国产精品亚洲美女久久久| 久久狼人影院| 黄片大片在线免费观看| 50天的宝宝边吃奶边哭怎么回事| 少妇粗大呻吟视频| 满18在线观看网站| 天堂影院成人在线观看| 免费在线观看影片大全网站| 天天添夜夜摸| 日韩欧美在线二视频| 欧美三级亚洲精品| videosex国产| 久久久国产成人精品二区| 岛国在线观看网站| 宅男免费午夜| bbb黄色大片| 色婷婷久久久亚洲欧美| 亚洲精品久久成人aⅴ小说| 国产精华一区二区三区| 久久人妻av系列| 一进一出抽搐gif免费好疼| 亚洲狠狠婷婷综合久久图片| 精品日产1卡2卡| 欧美精品亚洲一区二区| 一区二区日韩欧美中文字幕| 国产熟女午夜一区二区三区| 无限看片的www在线观看| 久久精品人妻少妇| 国产精品久久视频播放| 看黄色毛片网站| 黄片小视频在线播放| 久久草成人影院| 999精品在线视频| 黄频高清免费视频| 成人18禁高潮啪啪吃奶动态图| 91字幕亚洲| 男女之事视频高清在线观看| 亚洲成国产人片在线观看| 午夜激情av网站| 成人亚洲精品一区在线观看| 久久草成人影院| 欧美三级亚洲精品| 免费观看人在逋| 国产精品久久电影中文字幕| 在线观看www视频免费| 亚洲av成人av| 三级毛片av免费| 亚洲国产欧美网| 99re在线观看精品视频| 成在线人永久免费视频| 亚洲男人的天堂狠狠| 欧美日韩黄片免| 亚洲av电影在线进入| 三级毛片av免费| 亚洲五月色婷婷综合| 日韩欧美三级三区| 免费在线观看亚洲国产| 麻豆成人av在线观看| 女人高潮潮喷娇喘18禁视频| 人人澡人人妻人| 一a级毛片在线观看| 久久亚洲真实| 热re99久久国产66热| 日韩欧美一区二区三区在线观看| 精品无人区乱码1区二区| 亚洲熟女毛片儿| 亚洲性夜色夜夜综合| 久久久久国产一级毛片高清牌| 美女国产高潮福利片在线看| 欧美激情 高清一区二区三区| x7x7x7水蜜桃| 亚洲成人免费电影在线观看| 黑丝袜美女国产一区| 精品一区二区三区av网在线观看| 91成人精品电影| 91国产中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲国产精品成人综合色| 日韩大尺度精品在线看网址| 亚洲成人精品中文字幕电影| 日韩高清综合在线| netflix在线观看网站| 免费av毛片视频| 女人被狂操c到高潮| 成人国产一区最新在线观看| 18禁裸乳无遮挡免费网站照片 | www.自偷自拍.com| 欧美成人一区二区免费高清观看 | 视频区欧美日本亚洲| 久久午夜综合久久蜜桃| 亚洲自拍偷在线| 操出白浆在线播放| 国产99白浆流出| 又大又爽又粗| av免费在线观看网站| 搡老妇女老女人老熟妇| 亚洲 欧美一区二区三区| 日韩欧美三级三区| 高清在线国产一区| 国产又色又爽无遮挡免费看| 最新在线观看一区二区三区| 伊人久久大香线蕉亚洲五| 欧美色视频一区免费| 亚洲人成网站在线播放欧美日韩| 啪啪无遮挡十八禁网站| 欧美乱码精品一区二区三区| 久久久国产精品麻豆| 精品国内亚洲2022精品成人| 国产一卡二卡三卡精品| 人妻久久中文字幕网| 最好的美女福利视频网| 制服诱惑二区| 国产视频一区二区在线看| 亚洲av熟女| 一级毛片精品| 88av欧美| 天堂影院成人在线观看| 国产亚洲精品一区二区www| 美女午夜性视频免费| 国内精品久久久久久久电影| 侵犯人妻中文字幕一二三四区| 国产午夜精品久久久久久| 国产人伦9x9x在线观看| 一区二区日韩欧美中文字幕| 午夜成年电影在线免费观看| av在线播放免费不卡| 欧美中文综合在线视频| 免费电影在线观看免费观看| 十分钟在线观看高清视频www| 波多野结衣av一区二区av| 91麻豆av在线| 亚洲国产看品久久| 午夜精品久久久久久毛片777| 欧美一级毛片孕妇| 久久婷婷人人爽人人干人人爱| 欧美国产日韩亚洲一区| 亚洲成人精品中文字幕电影| 久久久精品国产亚洲av高清涩受| 韩国精品一区二区三区| 久久这里只有精品19| 国产亚洲精品第一综合不卡| 精品久久蜜臀av无| 草草在线视频免费看| svipshipincom国产片| 男人舔女人下体高潮全视频| 特大巨黑吊av在线直播 | 夜夜爽天天搞| 美女国产高潮福利片在线看| 亚洲国产欧美一区二区综合| 一二三四在线观看免费中文在| 免费在线观看影片大全网站| 久久精品91无色码中文字幕| 亚洲 欧美 日韩 在线 免费| 最近在线观看免费完整版| 999久久久精品免费观看国产| 久久九九热精品免费| 国产成+人综合+亚洲专区| 精品国产超薄肉色丝袜足j| 久久天躁狠狠躁夜夜2o2o| 99国产精品一区二区三区| 少妇裸体淫交视频免费看高清 | 欧美黑人巨大hd| 中文字幕av电影在线播放| 免费高清视频大片| 精品高清国产在线一区| 久久久水蜜桃国产精品网| 国产激情欧美一区二区| 欧美日韩中文字幕国产精品一区二区三区| 国产精品电影一区二区三区| 欧美黑人欧美精品刺激| 正在播放国产对白刺激| 久久99热这里只有精品18| 国产激情久久老熟女| 久久亚洲精品不卡| 欧美另类亚洲清纯唯美| 国产成人系列免费观看| aaaaa片日本免费| 国产免费男女视频| 日韩欧美免费精品| 欧美不卡视频在线免费观看 | 淫秽高清视频在线观看| 波多野结衣巨乳人妻| 欧美激情极品国产一区二区三区| 黄片播放在线免费| 亚洲国产欧洲综合997久久, | 午夜福利一区二区在线看| 中文资源天堂在线| 日韩欧美三级三区| 中文字幕最新亚洲高清| 久久精品国产亚洲av香蕉五月| 桃色一区二区三区在线观看| 国产精品久久久久久亚洲av鲁大| 97超级碰碰碰精品色视频在线观看| 日本免费一区二区三区高清不卡| 一级作爱视频免费观看| 12—13女人毛片做爰片一| 别揉我奶头~嗯~啊~动态视频| 999久久久精品免费观看国产| 黄片大片在线免费观看| 精品少妇一区二区三区视频日本电影| 丝袜人妻中文字幕| 白带黄色成豆腐渣| 亚洲一码二码三码区别大吗| 中文字幕精品免费在线观看视频| 男女之事视频高清在线观看| 91在线观看av| 国产精品av久久久久免费| 久久精品亚洲精品国产色婷小说| 精品久久久久久久毛片微露脸| 热re99久久国产66热| 午夜福利免费观看在线| 国产成人精品久久二区二区免费|