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

    糧食平房倉底部豎向壓力現(xiàn)場試驗及離散元模擬

    2022-01-26 06:36:04曾長女
    關鍵詞:倉壁縮尺糧堆

    曾長女,李 鑫,范 量,張 斌

    河南工業(yè)大學 土木工程學院,河南 鄭州 450001

    糧食是關系國家安全的戰(zhàn)略物資,是人民群眾的基本生活保障。我國建立了特殊的糧食長期儲存制度,建設了大量的糧食儲備倉以確保國家糧食安全,這些糧倉中平房倉占很大部分[1]。糧食壓力的確定對于糧食平房倉的設計和安全運營具有重要意義?!都Z食平房倉設計規(guī)范》(GB 50320—2014,以下簡稱《規(guī)范》)中提出,糧食引起的底部豎向壓力分布均勻,僅與糧堆高度和糧食容重相關。該理論將糧食作為一種理想的連續(xù)性介質(zhì),采用流體壓力理論計算糧堆豎向壓力,但忽略了糧食顆粒的離散特性[2],因而難以反映真實的糧堆豎向壓力分布規(guī)律。

    目前對筒倉側壓力的大小和分布研究較多[3-9]。糧食作為一種典型的散粒體介質(zhì),其底部豎向壓力分布與預測儲糧數(shù)量、平房倉設計等都相關。為了能準確掌握國家糧食儲備情況,張德賢等[10]通過現(xiàn)場試驗建立糧倉底部壓力與儲糧數(shù)量之間的關系,從而對儲糧數(shù)量進行監(jiān)測,確保國家儲糧安全。在對實倉[11-12]和模擬倉[13-14]的研究中,提出了糧堆底部豎向壓力呈非均勻分布,大致呈中間大、邊界小的趨勢,其不均勻程度隨糧堆高度增加而增大,糧堆單點底部壓力與糧堆高度之間近似呈線性關系。但由于現(xiàn)場試驗具有規(guī)模大、費用高、周期長、不易重復進行等缺點,且現(xiàn)場存在較多不確定因素會影響試驗結果,導致試驗數(shù)據(jù)有限,難以得到更多實用成果。為了彌補現(xiàn)場試驗的不足,部分學者采用離散元數(shù)值模擬方法研究。PFC(particle flow code)是一種基于離散元思想模擬不連續(xù)物體運動及相互作用的顆粒離散元法,近年來被應用于研究糧倉底部壓力。許啟鏗等[15]、陳家豪等[16]、楊真真等[17]利用 PFC3D軟件對散糧堆底部壓力進行了模擬分析,研究了糧堆豎向壓力、水平壓力等底部壓力的分布規(guī)律。這些研究大多針對糧堆堆積高度、糧堆內(nèi)部摩擦應力等因素對倉底壓力進行分析,而且都采用的是圓形顆粒模擬小麥,不符合實際顆粒形狀;此外,對倉底豎向壓力影響因素的系統(tǒng)研究也較少。因此,作者以現(xiàn)場試驗模型為原型,利用PFC3D離散元方法重現(xiàn)小麥顆粒,并將現(xiàn)場試驗結果與模擬現(xiàn)場裝料的數(shù)值模擬結果作對比,系統(tǒng)分析影響糧堆底部豎向壓力的內(nèi)因和外因,為糧倉底部壓力分布的正確描述提供理論支持。

    1 糧食顆粒單元改進

    2 離散元數(shù)值模擬

    2.1 數(shù)值縮尺模型相似比分析

    由于實際模型尺寸大,在數(shù)值建模時,若按1∶ 1進行建模會導致數(shù)值模型生成的顆粒數(shù)量龐大,導致數(shù)值模擬無法進行。因此,利用Zhao等[18]提出相似理論來反映不同模型尺寸之間的內(nèi)在本質(zhì)聯(lián)系,建立縮尺模型可較準確地表示原始物理模型的力學行為。

    保證建立的原物理模型和縮尺模型的幾何尺寸相似,即兩個模型的顆粒半徑、模型尺寸、顆粒重疊量等按照等比例進行縮放。

    式中:Rw和Rs分別為原物理模型和縮尺模型的粒子半徑;Lw和Ls分別為原物理模型和縮尺模型的幾何長度;uw和us分別為原物理模型和縮尺模型的顆粒重疊量。

    建立縮尺模型時,縮尺模型與原物理模型的應力、應變必須保持一致[19]。應力(σ)、應變(ε)可表示:

    式中: F(u,R)為相鄰顆粒之間所有接觸力的合力;D為顆粒特征長度,D=2R;A為顆粒特征面積,A=D2。由式(1)可知,可保證縮尺模型與原物理模型中顆粒應變相等。在三維中,縮尺模型與原物理模型特征面積比值為λA=λ2,則可知縮尺模型與原物理模型的顆粒間接觸力的合力比值λF=λ2。

    式中:F(u,R)w為原物理模型顆粒之間所有接觸力的合力;F(u,R)s為縮尺模型顆粒之間所有接觸力的合力。

    在離散元中,顆粒的運動方程滿足牛頓第二定律[19]。

    (4)

    在準靜態(tài)狀態(tài)下,式(4)中的慣性項可以忽略不計[19],即:

    F(u,R)=Q(t)。

    (5)

    為了保證建立的縮尺模型與原物理模型具有相同的力,則兩個模型中對應的力必須相等或相差一個固定倍數(shù)[20]。

    式中:t為常數(shù);Q(t)w為原物理模型顆粒外部施加的力;Q(t)s為縮尺模型顆粒外部施加的力。

    故在縮尺模型中,若將試驗模型尺寸等比例縮小λ,則模擬模型顆粒密度應該為原試驗模型顆粒密度的λ-1。

    2.2 PFC3D模型建立

    2.2.1 工程概況

    本試驗倉是方形平房倉實倉,倉體長45 m,寬28 m,高12 m。根據(jù)對稱性,選取整個糧倉的1/4作為試驗區(qū)域,糧食最終裝滿全倉。在試驗區(qū)域倉底安裝44個振弦式壓力傳感器和6個重復檢測試點監(jiān)測糧堆底部豎向壓力[21]。壓力傳感器分為4行,在倉壁2距離0.1、9、16 m位置處分別布置間隔為1 m或2 m的3列壓力傳感器,分別標注為R3、R2和R1。在倉壁1距離7 m位置處布置間隔為1 m或2 m的一行壓力傳感器,標注為R4。傳感器布置如圖2所示。在倉底安裝振弦式壓力傳感器,讀取穩(wěn)定初始值,然后開始裝糧。裝糧時間為每天8:00至23:00,歷時5 d完成。裝糧完成后靜置儲存。在試驗過程中,每隔0.5 h讀取并記錄每個壓力傳感器上的豎向壓力值。

    圖2 現(xiàn)場試驗壓力傳感器的布置

    2.2.2 模型細觀參數(shù)確定

    基于上述現(xiàn)場試驗以及相似理論,將試驗倉尺寸縮小10%進行數(shù)值建模,即所建模型的長、寬、高分別為4.5、2.8、1.2 m。為使模型倉達到原型倉所需的重力場,按上述縮尺模型相似比結論,可通過增大散體顆粒的密度來實現(xiàn),即模型倉中散體顆粒密度為原型倉中顆粒密度的10倍,可使模型倉和原型倉內(nèi)應力狀態(tài)一致。

    在PFC3D中,顆粒與顆粒、顆粒與墻體相互接觸都會產(chǎn)生相互作用[22]。顆粒之間的相互作用通過模擬顆粒間的接觸來實現(xiàn),在每一個特定接觸上的本構模型可以分為線性接觸模型、接觸黏結模型等。線性接觸模型適用于無黏性模型,根據(jù)糧食顆粒離散性特點,本文中采用線性接觸模型,主要細觀參數(shù):顆粒法向接觸剛度kn、法向與切向剛度比kratio和摩擦系數(shù)μ。采用參數(shù)標定方法確定顆粒試樣的細觀參數(shù)[23],當模擬值與試驗值的應力-應變吻合時,認為當前選擇的參數(shù)比較合適。表1為PFC3D模型所選用的細觀參數(shù)。

    表1 PFC3D模型細觀參數(shù)Table 1 Mesoscopic parameters of PFC3D model

    2.2.3 儲糧過程模擬

    建模時,倉壁和倉底采用墻單元表示,糧食顆粒采用顆粒簇clump單元表示。PFC3D提供的測量球命令作為虛擬傳感器對底部豎向壓力進行測量,測量球半徑為0.15m,測量球放置位置如圖3a所示,可方便地監(jiān)測全倉的糧食底部豎向壓力。為模擬人工平糧狀態(tài),該模型倉內(nèi)散體顆粒的生成采用在模型底部生成顆粒,在重力作用下讓顆粒達到平衡狀態(tài),再生成下一層顆粒,如此循環(huán)直到達到預定填料高度,散體顆粒上部為平頂。本次數(shù)值模擬分4次裝滿,每層高度分別為0.20、0.40、0.60、0.75m。先在模型頂部生成半徑為0.020~0.024m的小球顆粒,再利用顆粒簇clump單元進行替換,最后設定重力場,運行一定時步后,觀察最大不平衡力,當其趨近于零時,顆粒堆積穩(wěn)定密實,開始下次裝料,循環(huán)重復此過程4次,達到預定高度時停止裝料,結果如圖3b所示。

    圖3 PFC3D建模模型Fig.3 PFC3D modeling model diagram

    3 糧食底部豎向壓力影響因素分析

    在驗證模型的基礎上,改變影響因素,進一步分析糧堆底部豎向壓力分布影響因素。所建模型具有對稱性,糧堆豎向壓力分布也具有對稱性,因此取坐標軸一半進行分析。

    3.1 顆粒形狀影響

    圖4所示為球形顆粒、橢球顆粒的底部豎向壓力分布。由圖4可知,球形顆粒整體趨勢與橢球顆粒相近,但橢球顆粒在距倉壁最遠位置處的底部豎向壓力值較大,在靠近倉壁位置處底部壓力值較小,而且隨著糧堆高度的增加,距倉壁最遠位置與倉壁位置處的底部壓力差值越大。對比2種顆粒形狀的力鏈數(shù),球形顆粒的力鏈數(shù)為355 547,橢球顆粒的力鏈數(shù)為562 822,這是由于橢球顆粒比球形顆粒更容易發(fā)生自鎖現(xiàn)象[24],使橢球顆粒在內(nèi)部形成能夠有效抵御外界的顆粒鏈,不易受外界影響。

    圖4 不同顆粒形狀對糧堆底部壓力

    配位數(shù)為某一顆粒與其他顆粒接觸的顆粒數(shù)量,是表示堆積密實度的細觀結構參數(shù)。不同顆粒形狀配位數(shù)變化如圖5所示。由圖5a可知,配位數(shù)在X軸上距倉壁距離的增大而增大,與糧堆底部豎向壓力分布相對應。將倉底所有測點處配位數(shù)取平均值,由圖5b可知,隨著堆糧高度的增加,顆粒平均配位數(shù)隨之增加,使糧堆更加密實。模型達到平衡狀態(tài)時,球形顆粒配位數(shù)為6.12,Aparicio等[25]從理論上證明了單一粒徑球體準靜態(tài)下堆積的配位數(shù)為6,與本文模擬值相近;橢球顆粒配位數(shù)為10.36。由此可知,不同顆粒形狀在不同堆糧高度下,其配位數(shù)不同,橢球顆??梢杂行г龃箢w粒的配位數(shù),可以使糧堆在堆積過程中更加密實,不易發(fā)生波動。

    圖5 配位數(shù)的變化

    3.2 裝料方式影響

    為模擬裝料方式的影響,采用2種顆粒生成方式:一種在模型底部一層一層生成顆粒,如此循環(huán)直到達到預定填料高度;另一種在模型頂部生成顆粒,在重力作用下讓顆粒下落至倉底,待顆粒平衡后,再生成下一層顆粒,如此循環(huán)直到達到預定填料高度。2種顆粒生成方式的橢球顆粒在不同糧堆高度下豎向壓力分布曲線如圖6所示。由圖6可知,頂部生成顆粒方式使底部豎向壓力曲線浮動較小,底部豎向壓力分布較為均勻;底部生成顆粒方式對底部豎向壓力分布影響較大,且受倉壁摩擦系數(shù)影響較大。

    注:底部生成顆粒簡稱底顆粒,頂部生成顆粒簡稱頂顆粒。圖7同。

    孔隙率、配位數(shù)隨堆糧高度變化如圖7所示。由圖7可知,孔隙率隨堆糧高度增加而減小,配位數(shù)隨堆糧高度增加而增加。堆糧高度達到0.75m時底部生成顆粒的平均孔隙率與平均配位數(shù)分別為0.392、10.360,頂部生成顆粒的平均孔隙率與平均配位數(shù)分別為0.349、11.689。頂部生成顆粒在重力作用下下落至模型倉底部,導致糧堆的孔隙率減小、配位數(shù)增大,孔隙率與配位數(shù)變化成反比,與齊陽等[26]得到的孔隙率與配位數(shù)變化趨勢相似??勺C明頂部生成顆粒的生成方式使糧堆更加密集、接觸力分布更均勻,對糧堆底部豎向壓力影響更小。對比這2種生成方式的接觸力鏈數(shù),底部生成顆粒的生成填充方式接觸力鏈數(shù)為573 923,而頂部生成顆粒的生成方式接觸力鏈數(shù)為688 852,也可證明頂部生成顆粒的生成方式使糧堆更加密集。因此,顆粒生成方式不同,導致糧堆結構的密實狀態(tài)不同,會影響糧堆底部豎向壓力的分布。

    圖7 孔隙率、配位數(shù)隨堆糧高度變化

    3.3 倉壁摩擦系數(shù)影響

    考慮倉壁摩擦系數(shù)為0.0、0.2和0.4對糧堆底部豎向壓力影響,不同倉壁摩擦系數(shù)糧堆底部豎向壓力分布如圖8所示。由圖8可知,倉壁摩擦系數(shù)為0.0的底部豎向壓力曲線波動較大,呈多個波峰和多個波谷狀態(tài);倉壁摩擦系數(shù)為0.4的底部豎向壓力隨著距倉壁距離的增加呈增加趨勢。由于不同摩擦系數(shù)倉壁對顆粒的摩擦作用,導致不同倉壁對糧堆底部豎向壓力的影響不同。

    圖8 不同倉壁摩擦系數(shù)糧堆底部豎向壓力分布

    X向對稱軸中心位置、倉壁位置處底部豎向壓力在不同倉壁摩擦系數(shù)作用下隨堆糧高度變化曲線如圖9所示。由圖9可知,在模型中心位置處,隨著倉壁摩擦系數(shù)增大,中心底部豎向壓力值增大;在倉壁位置處,隨著倉壁摩擦系數(shù)越大,邊界底部豎向壓力值越小。由此表明,倉壁摩擦系數(shù)對邊界底部壓力的分布影響顯著,對中心位置處的底部壓力影響較小,底部壓力分布與倉壁摩擦系數(shù)、距離倉壁距離有關,與文獻[16]的結論相似。

    圖9 X向對稱軸中心位置、倉壁位置底部壓力隨堆糧高度變化

    4 現(xiàn)場試驗與數(shù)值模擬結果

    結合上述糧堆底部豎向壓力的影響因素分析,在倉壁摩擦系數(shù)為0.4的模型倉頂部生成橢球顆粒,在重力作用下下落至模型倉底部,并將數(shù)值模擬與現(xiàn)場試驗結果進行對比分析。數(shù)值模擬與現(xiàn)場測點數(shù)據(jù)結果對比如圖10所示。模擬采用的測量球半徑較大,并非與現(xiàn)場試驗布置的傳感器位置一一對應,因此會產(chǎn)生偏差。利用《規(guī)范》中糧堆豎向壓力公式pvk=γh( γ 為糧食容重,kN/m3;h為堆糧高度,m)計算最大糧堆高度下底部豎向壓力值作為理論值,由圖10可知,《規(guī)范》獲得的理論值均比現(xiàn)場試驗測試值和數(shù)值模擬值大。隨著距垂直倉壁距離的增加,數(shù)值模擬與試驗監(jiān)測的糧堆底部豎向壓力均先逐漸增加,之后達到穩(wěn)定狀態(tài)。就整體趨勢而言,數(shù)值模擬的底部豎向壓力值與試驗值具有相同的變化趨勢,可利用本文建立的模型進行糧堆底部豎向壓力變化趨勢研究。

    圖10 數(shù)值模擬與現(xiàn)場試驗對比分析

    為了定量地表述倉底糧食豎向壓力分布規(guī)律,采用壓力分布系數(shù)表述不同位置處糧堆底部豎向壓力。定義壓力分布系數(shù)λV為各測點位置糧堆底部豎向壓力(V)與對應列中最大的糧堆底部豎向壓力值(Vmax)之比。

    由圖11可知,壓力分布系數(shù)λV隨距離垂直對應倉壁距離的增加而增加,在垂直對應倉壁0~4m范圍內(nèi)壓力分布系數(shù)增長較快。在對應倉壁大于8m后,兩個測點之間壓力分布系數(shù)大部分在0.9~1.0之間,增長幅度較小,因此該位置處倉壁對糧堆底部豎向壓力影響較小。需要說明的是,本現(xiàn)場試驗中靠近中部3個點有異常,主要是受中間柱等因素的影響。

    圖11 現(xiàn)場試驗底部糧食豎向壓力分布系數(shù)

    將現(xiàn)場試驗R3、R2、R1和R4列的數(shù)值模擬結果與文獻[16]數(shù)值模擬結果對比,如圖12所示。壓力分布系數(shù)λV隨著距離相應壁面距離的增加而逐漸增大,隨后又緩慢下降逐漸趨于穩(wěn)定。距倉壁0~5m范圍內(nèi),壓力分布系數(shù)增長較快,與文獻[16]增長趨勢相同;距倉壁8m左右,增長幅度較小,逐步趨于穩(wěn)定,與現(xiàn)場試驗結論相似。

    圖12 數(shù)值模擬底部糧食豎向壓力分布系數(shù)

    由現(xiàn)場試驗和數(shù)值模擬試驗可知,距倉壁距離對倉底糧食壓力分布影響較大,根據(jù)上述分析,糧倉與糧食之間摩擦系數(shù)、裝料方式、顆粒形狀等都會對糧堆底部豎向壓力產(chǎn)生影響,顯然裝料方式、倉壁摩擦系數(shù)影響更大,尤其是倉壁摩擦系數(shù)的影響獲得了眾多學者的關注。由本文的試驗可得,在試驗范圍內(nèi),距倉壁距離大于8m后,墻體對糧食底部豎向壓力的影響就很小了。整體而言,現(xiàn)場試驗結果比數(shù)值模擬結果大,且這兩者結果均小于《規(guī)范》推薦的理論值,有必要積累更多的現(xiàn)場數(shù)據(jù),對現(xiàn)有規(guī)范相關內(nèi)容進行適當?shù)男拚?/p>

    5 結論

    以現(xiàn)場平房倉裝料試驗為基礎,利用PFC3D軟件對平房倉進行數(shù)值模擬,研究了顆粒形狀、倉壁摩擦系數(shù)和裝料方式對糧堆底部豎向壓力分布的影響。

    離散元軟件PFC3D建立的橢球顆粒能更好地模擬小麥在填料過程中糧堆底部豎向壓力變化趨勢。球形顆粒配位數(shù)為6.12,橢球顆粒配位數(shù)為10.36,橢球顆粒更容易發(fā)生自鎖現(xiàn)象。顆粒不同的生成方式會影響糧堆結構的密實狀態(tài),從而影響糧堆底部豎向壓力的分布。頂部生成顆粒方式會導致糧堆孔隙率小、配位數(shù)大,使得糧堆更加密集。

    糧堆底部豎向壓力分布與倉壁摩擦系數(shù)以及距離倉壁距離有關。不同倉壁摩擦系數(shù)對倉底豎向壓力分布影響不同,倉壁摩擦系數(shù)越大,對邊界底部壓力的分布影響越顯著,對中心位置處的底部壓力影響越小。

    壓力分布系數(shù)能很好地表示不同位置處糧堆底部豎向壓力的分布情況。壓力分布系數(shù)在倉壁處較小,隨距倉壁距離的增加而逐漸增大。現(xiàn)場試驗和模擬試驗均表明在距倉壁距離小于5m的范圍內(nèi),壓力分布系數(shù)增長較快,到達8m后壓力分布系數(shù)趨于穩(wěn)定。

    本文獲得的現(xiàn)場試驗和數(shù)值分析結果,將為今后規(guī)范修訂倉底壓力提供可靠的理論依據(jù)。

    致謝:感謝云南省糧油工業(yè)有限公司為試驗提供幫助。

    猜你喜歡
    倉壁縮尺糧堆
    地下裝配式大直徑組合筒倉倉壁穩(wěn)定性分析
    河南科技(2023年16期)2023-10-10 09:11:08
    爆炸荷載作用下鋼筋混凝土構件縮尺效應的數(shù)值模擬研究
    結構工程師(2022年4期)2022-10-14 06:47:12
    箱梁渦振的縮尺效應及振幅修正研究
    裝配式地下糧倉鋼板-混凝土組合倉壁整體結構力學性能分析
    尺度效應對噴水推進系統(tǒng)進出口流場及推力影響分析
    鋼筋混凝土土體地下糧倉倉壁結構選型初探
    儲料豎向壓力對糧倉中小麥糧堆濕熱傳遞的影響
    溫濕度對稻谷糧堆結露的影響及實倉結露預警
    地下庫糧堆充氮技術研究
    土石堆積體介質(zhì)力學特性縮尺試驗研究
    日本精品一区二区三区蜜桃| 亚洲av成人不卡在线观看播放网| 成人高潮视频无遮挡免费网站| 人妻丰满熟妇av一区二区三区| 啦啦啦观看免费观看视频高清| av天堂中文字幕网| 国产真人三级小视频在线观看| 听说在线观看完整版免费高清| 人妻丰满熟妇av一区二区三区| 黄色丝袜av网址大全| 两人在一起打扑克的视频| 在线免费观看不下载黄p国产 | 哪里可以看免费的av片| 国产精品久久视频播放| 日本在线视频免费播放| 国产精品98久久久久久宅男小说| 精品国产超薄肉色丝袜足j| 国产精品1区2区在线观看.| 国产黄片美女视频| 亚洲专区中文字幕在线| 色综合站精品国产| 最新在线观看一区二区三区| 亚洲av片天天在线观看| 老熟妇仑乱视频hdxx| 日韩欧美三级三区| 在线观看午夜福利视频| www.www免费av| 国产精品久久久久久人妻精品电影| 岛国在线免费视频观看| 最近最新中文字幕大全免费视频| 午夜福利在线观看吧| 亚洲精华国产精华精| 成人av在线播放网站| 午夜激情福利司机影院| 久久草成人影院| 欧美中文综合在线视频| 国产v大片淫在线免费观看| 日本撒尿小便嘘嘘汇集6| 中出人妻视频一区二区| 国产av一区在线观看免费| 国内久久婷婷六月综合欲色啪| 韩国av一区二区三区四区| 国产一区二区三区视频了| 成年女人看的毛片在线观看| 在线观看免费视频日本深夜| 老汉色∧v一级毛片| 深夜精品福利| 午夜激情欧美在线| 国产亚洲精品综合一区在线观看| 日本与韩国留学比较| 综合色av麻豆| 少妇的丰满在线观看| 淫妇啪啪啪对白视频| 国内精品久久久久精免费| 久久久精品大字幕| 91老司机精品| 人人妻人人看人人澡| 免费在线观看影片大全网站| 国产一区二区激情短视频| 真人一进一出gif抽搐免费| 午夜久久久久精精品| 亚洲精品一卡2卡三卡4卡5卡| 男人舔女人下体高潮全视频| 宅男免费午夜| 天堂√8在线中文| 国产一区二区激情短视频| 精品久久久久久久末码| 亚洲欧美日韩高清在线视频| 伊人久久大香线蕉亚洲五| 国产午夜精品久久久久久| 禁无遮挡网站| 日韩欧美在线乱码| 人人妻人人澡欧美一区二区| 精品国产三级普通话版| 91麻豆精品激情在线观看国产| 日韩欧美一区二区三区在线观看| 夜夜爽天天搞| 精品久久久久久久久久免费视频| 男人舔女人的私密视频| 嫩草影院精品99| 精品无人区乱码1区二区| 日韩欧美在线乱码| 在线视频色国产色| 草草在线视频免费看| 午夜福利欧美成人| 精品国产三级普通话版| 色精品久久人妻99蜜桃| 99久久成人亚洲精品观看| 在线观看免费视频日本深夜| 久久精品国产清高在天天线| 欧美成狂野欧美在线观看| 国产在线精品亚洲第一网站| 国产乱人视频| av天堂在线播放| 精品久久久久久,| 国产一区二区在线观看日韩 | 免费观看人在逋| 免费看日本二区| 日本熟妇午夜| 免费电影在线观看免费观看| 国产真实乱freesex| 精品久久久久久久久久久久久| 亚洲欧美日韩高清在线视频| 久久久久久久久中文| 国产精品一区二区三区四区免费观看 | 国产一区二区三区视频了| 亚洲18禁久久av| 国产淫片久久久久久久久 | 女人高潮潮喷娇喘18禁视频| 欧美在线黄色| 99久久精品热视频| 日韩欧美一区二区三区在线观看| 最近最新免费中文字幕在线| 99精品在免费线老司机午夜| 色在线成人网| 中亚洲国语对白在线视频| 欧美大码av| 亚洲成人久久性| 九九热线精品视视频播放| 别揉我奶头~嗯~啊~动态视频| 色在线成人网| 亚洲成人久久性| 亚洲乱码一区二区免费版| 亚洲无线在线观看| 人妻久久中文字幕网| 国产精品99久久99久久久不卡| 国产私拍福利视频在线观看| 亚洲国产色片| 国产免费男女视频| 男女做爰动态图高潮gif福利片| 一二三四社区在线视频社区8| 亚洲av电影在线进入| 国产成人精品无人区| 非洲黑人性xxxx精品又粗又长| 国产亚洲精品av在线| 成人三级黄色视频| 国产亚洲精品久久久com| 看免费av毛片| 欧美又色又爽又黄视频| 叶爱在线成人免费视频播放| 在线免费观看不下载黄p国产 | 精品久久久久久成人av| 12—13女人毛片做爰片一| 十八禁网站免费在线| 麻豆成人午夜福利视频| 午夜日韩欧美国产| 九九久久精品国产亚洲av麻豆 | 手机成人av网站| 午夜福利高清视频| 久久久国产欧美日韩av| 国产精品一及| 日韩高清综合在线| 99热这里只有是精品50| 2021天堂中文幕一二区在线观| 日本免费a在线| 亚洲第一电影网av| 国产极品精品免费视频能看的| 亚洲精品美女久久久久99蜜臀| 男插女下体视频免费在线播放| 欧美黄色淫秽网站| 日本三级黄在线观看| 观看免费一级毛片| 国产av在哪里看| 母亲3免费完整高清在线观看| 伦理电影免费视频| 国产av在哪里看| 成人午夜高清在线视频| 怎么达到女性高潮| 天天躁狠狠躁夜夜躁狠狠躁| 69av精品久久久久久| 亚洲成人久久爱视频| 蜜桃久久精品国产亚洲av| 中文在线观看免费www的网站| 日本熟妇午夜| 嫩草影视91久久| 久久亚洲精品不卡| 午夜精品在线福利| 熟妇人妻久久中文字幕3abv| 亚洲av美国av| 女人被狂操c到高潮| 婷婷精品国产亚洲av| 一二三四社区在线视频社区8| 亚洲成a人片在线一区二区| 亚洲成av人片在线播放无| 成年女人看的毛片在线观看| 少妇人妻一区二区三区视频| 美女黄网站色视频| 黄色片一级片一级黄色片| 两性午夜刺激爽爽歪歪视频在线观看| 草草在线视频免费看| 亚洲色图 男人天堂 中文字幕| 精品久久久久久久人妻蜜臀av| 亚洲成人中文字幕在线播放| 两性夫妻黄色片| 高潮久久久久久久久久久不卡| 一个人看的www免费观看视频| 999久久久国产精品视频| 成人三级黄色视频| 亚洲成av人片免费观看| 在线十欧美十亚洲十日本专区| 国产精品久久久久久精品电影| 国产精华一区二区三区| 欧美乱色亚洲激情| 国产精品,欧美在线| 亚洲人与动物交配视频| 在线观看66精品国产| 在线观看一区二区三区| 国产精品av久久久久免费| 久久精品91无色码中文字幕| 99久久精品国产亚洲精品| 欧美不卡视频在线免费观看| 舔av片在线| 午夜成年电影在线免费观看| 男人舔女人的私密视频| 男人的好看免费观看在线视频| 国产亚洲精品综合一区在线观看| 亚洲激情在线av| 1024手机看黄色片| 国产综合懂色| 亚洲av五月六月丁香网| 神马国产精品三级电影在线观看| 国产精品免费一区二区三区在线| 精品一区二区三区视频在线 | 国内少妇人妻偷人精品xxx网站 | 我的老师免费观看完整版| 黄频高清免费视频| 成人欧美大片| 大型黄色视频在线免费观看| 美女黄网站色视频| 亚洲欧美日韩卡通动漫| 国产日本99.免费观看| 久久久久免费精品人妻一区二区| 亚洲av美国av| 国产精品 欧美亚洲| 一进一出抽搐gif免费好疼| 日本a在线网址| 啦啦啦观看免费观看视频高清| www日本黄色视频网| 久久精品91无色码中文字幕| 国产人伦9x9x在线观看| 国产毛片a区久久久久| 一区二区三区高清视频在线| 欧美午夜高清在线| 99精品久久久久人妻精品| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲国产欧美人成| 一本精品99久久精品77| 日韩 欧美 亚洲 中文字幕| av女优亚洲男人天堂 | 日韩高清综合在线| 国产精品电影一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 97超视频在线观看视频| 俄罗斯特黄特色一大片| 成人鲁丝片一二三区免费| 亚洲 欧美 日韩 在线 免费| а√天堂www在线а√下载| 精品一区二区三区av网在线观看| 久久久国产欧美日韩av| 岛国在线免费视频观看| 日本撒尿小便嘘嘘汇集6| 99热这里只有是精品50| 午夜两性在线视频| 免费av毛片视频| 免费在线观看成人毛片| 久久香蕉精品热| 舔av片在线| 天堂网av新在线| 亚洲成人精品中文字幕电影| 麻豆成人午夜福利视频| 动漫黄色视频在线观看| 亚洲狠狠婷婷综合久久图片| 国产亚洲欧美98| 一本综合久久免费| 欧美+亚洲+日韩+国产| 欧美一区二区国产精品久久精品| 亚洲人成伊人成综合网2020| 一区二区三区激情视频| 亚洲精品久久国产高清桃花| 无限看片的www在线观看| АⅤ资源中文在线天堂| 国产高清三级在线| 18禁国产床啪视频网站| 欧美日韩亚洲国产一区二区在线观看| 小说图片视频综合网站| 深夜精品福利| 后天国语完整版免费观看| 亚洲欧洲精品一区二区精品久久久| 一区二区三区激情视频| 我要搜黄色片| 亚洲av片天天在线观看| 哪里可以看免费的av片| 色尼玛亚洲综合影院| 最好的美女福利视频网| 操出白浆在线播放| 不卡一级毛片| 99久久无色码亚洲精品果冻| 亚洲 欧美一区二区三区| 日韩免费av在线播放| www.999成人在线观看| 性色av乱码一区二区三区2| 伦理电影免费视频| 一级a爱片免费观看的视频| 可以在线观看的亚洲视频| 九九在线视频观看精品| 老鸭窝网址在线观看| 搡老熟女国产l中国老女人| 亚洲精品美女久久久久99蜜臀| 蜜桃久久精品国产亚洲av| 亚洲精品在线观看二区| 69av精品久久久久久| 天堂影院成人在线观看| 欧美日韩黄片免| 长腿黑丝高跟| 欧美日韩一级在线毛片| 亚洲国产精品合色在线| 男女之事视频高清在线观看| 久久精品亚洲精品国产色婷小说| 国内毛片毛片毛片毛片毛片| 国产成+人综合+亚洲专区| 国产一区二区激情短视频| 精品久久久久久久久久免费视频| 97超视频在线观看视频| 黑人欧美特级aaaaaa片| 国产91精品成人一区二区三区| 久9热在线精品视频| 天堂av国产一区二区熟女人妻| 一进一出抽搐gif免费好疼| 国产一区二区激情短视频| 国产精品久久电影中文字幕| 色噜噜av男人的天堂激情| 亚洲人成伊人成综合网2020| 亚洲真实伦在线观看| 日韩有码中文字幕| 日韩欧美 国产精品| 欧美在线一区亚洲| 法律面前人人平等表现在哪些方面| 国产精品1区2区在线观看.| 久久中文看片网| 欧美成狂野欧美在线观看| 少妇的逼水好多| 免费高清视频大片| 国产精品一区二区免费欧美| 国产私拍福利视频在线观看| 国产 一区 欧美 日韩| 国产精品久久久久久精品电影| 亚洲avbb在线观看| 国产一区在线观看成人免费| 午夜激情福利司机影院| 黄色成人免费大全| 变态另类丝袜制服| 精品99又大又爽又粗少妇毛片 | 在线免费观看不下载黄p国产 | 亚洲男人的天堂狠狠| 久久精品人妻少妇| 久久久久久久精品吃奶| 国产高清三级在线| 国产精品综合久久久久久久免费| 国产精品女同一区二区软件 | 婷婷六月久久综合丁香| 欧美中文综合在线视频| 日本 av在线| 色av中文字幕| 婷婷六月久久综合丁香| 亚洲欧洲精品一区二区精品久久久| 国产一区在线观看成人免费| 女警被强在线播放| 欧美激情久久久久久爽电影| 国产欧美日韩一区二区三| 亚洲欧美日韩高清专用| 亚洲欧美日韩东京热| 国产精品久久久人人做人人爽| 99久国产av精品| 在线视频色国产色| 久久草成人影院| 91av网一区二区| 免费看日本二区| 国产三级在线视频| 成年免费大片在线观看| 日韩欧美精品v在线| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美日韩东京热| 香蕉av资源在线| 美女cb高潮喷水在线观看 | 久久久久久久久久黄片| 一进一出抽搐动态| 在线观看免费视频日本深夜| 男插女下体视频免费在线播放| 女同久久另类99精品国产91| 嫁个100分男人电影在线观看| 国产黄片美女视频| 桃红色精品国产亚洲av| 人妻夜夜爽99麻豆av| 男人的好看免费观看在线视频| 午夜两性在线视频| 观看美女的网站| 日韩欧美国产在线观看| 黄色 视频免费看| 日本免费a在线| 日本三级黄在线观看| 国产精品综合久久久久久久免费| 母亲3免费完整高清在线观看| 免费在线观看视频国产中文字幕亚洲| 欧美成人性av电影在线观看| 在线观看午夜福利视频| 色哟哟哟哟哟哟| 99精品久久久久人妻精品| 国产99白浆流出| cao死你这个sao货| 精品一区二区三区av网在线观看| 亚洲精品一区av在线观看| 母亲3免费完整高清在线观看| 国产精品野战在线观看| 男女床上黄色一级片免费看| 精品久久久久久久人妻蜜臀av| 国产一区二区三区在线臀色熟女| 免费看光身美女| 91久久精品国产一区二区成人 | 国产一区二区三区视频了| 两性夫妻黄色片| 欧美日韩黄片免| 国产成人精品无人区| 天堂av国产一区二区熟女人妻| 99re在线观看精品视频| 18美女黄网站色大片免费观看| 亚洲av片天天在线观看| 国产免费av片在线观看野外av| 最新在线观看一区二区三区| 国产真人三级小视频在线观看| 99热这里只有精品一区 | 很黄的视频免费| 女人高潮潮喷娇喘18禁视频| 欧美极品一区二区三区四区| 免费电影在线观看免费观看| 精品国产乱子伦一区二区三区| 美女黄网站色视频| 国产激情久久老熟女| 一进一出抽搐gif免费好疼| 少妇裸体淫交视频免费看高清| 在线免费观看不下载黄p国产 | 亚洲人成伊人成综合网2020| 亚洲国产精品999在线| 精品久久久久久成人av| 久久天堂一区二区三区四区| 国产成人欧美在线观看| 成年女人看的毛片在线观看| 嫩草影院精品99| 久久人妻av系列| 日本与韩国留学比较| 真人一进一出gif抽搐免费| 我的老师免费观看完整版| 欧美日韩福利视频一区二区| 婷婷精品国产亚洲av在线| 亚洲成a人片在线一区二区| 国产男靠女视频免费网站| 国产单亲对白刺激| 婷婷六月久久综合丁香| 中文字幕久久专区| 精品久久久久久久毛片微露脸| 国产熟女xx| 噜噜噜噜噜久久久久久91| 制服人妻中文乱码| 亚洲乱码一区二区免费版| 嫁个100分男人电影在线观看| 欧美日本视频| 99在线人妻在线中文字幕| 波多野结衣巨乳人妻| 极品教师在线免费播放| 国产精品98久久久久久宅男小说| 香蕉国产在线看| 久久欧美精品欧美久久欧美| 狂野欧美激情性xxxx| 国产伦精品一区二区三区视频9 | 国产精品野战在线观看| 国产一区二区三区在线臀色熟女| 精品一区二区三区av网在线观看| 午夜福利欧美成人| bbb黄色大片| 一级黄色大片毛片| 黄频高清免费视频| 在线观看66精品国产| 两性夫妻黄色片| 91字幕亚洲| 国产v大片淫在线免费观看| 日本 av在线| 精品国产乱子伦一区二区三区| 视频区欧美日本亚洲| 熟妇人妻久久中文字幕3abv| 又大又爽又粗| 久久久久久人人人人人| 99久久综合精品五月天人人| 男人舔女人下体高潮全视频| aaaaa片日本免费| 黄色 视频免费看| 色精品久久人妻99蜜桃| 国产亚洲精品av在线| 精品久久久久久久毛片微露脸| 757午夜福利合集在线观看| 亚洲五月婷婷丁香| 欧美又色又爽又黄视频| 一卡2卡三卡四卡精品乱码亚洲| 国产私拍福利视频在线观看| 看黄色毛片网站| 91在线精品国自产拍蜜月 | 国产蜜桃级精品一区二区三区| 99热精品在线国产| 搡老妇女老女人老熟妇| 久久久久久久久中文| 天堂√8在线中文| 亚洲第一电影网av| 久久久国产欧美日韩av| 亚洲人成电影免费在线| 亚洲人成网站在线播放欧美日韩| 成人无遮挡网站| 一区二区三区高清视频在线| 少妇的丰满在线观看| 他把我摸到了高潮在线观看| 国产伦精品一区二区三区四那| 精品国产乱码久久久久久男人| 曰老女人黄片| 亚洲欧美日韩高清专用| 精品人妻1区二区| 99久久精品国产亚洲精品| 日韩欧美国产一区二区入口| 午夜激情欧美在线| 国产男靠女视频免费网站| 国产97色在线日韩免费| 成人三级黄色视频| 午夜福利成人在线免费观看| 久久久国产欧美日韩av| 国产视频内射| 国产三级中文精品| 午夜精品久久久久久毛片777| 国内少妇人妻偷人精品xxx网站 | 99热只有精品国产| 香蕉久久夜色| a在线观看视频网站| 日本五十路高清| 国产av一区在线观看免费| 综合色av麻豆| 久久国产精品影院| 99久久精品国产亚洲精品| 欧美午夜高清在线| 91九色精品人成在线观看| 成人无遮挡网站| 99久久久亚洲精品蜜臀av| 成人av在线播放网站| 一进一出抽搐gif免费好疼| 一本一本综合久久| 国产99白浆流出| 黄片大片在线免费观看| 身体一侧抽搐| 欧美中文日本在线观看视频| 天堂动漫精品| 亚洲国产精品合色在线| 在线观看日韩欧美| 国产高清激情床上av| 99国产精品99久久久久| 国产精品久久视频播放| 免费看a级黄色片| 九色国产91popny在线| 免费大片18禁| 欧美日韩综合久久久久久 | 久久久精品大字幕| 男女床上黄色一级片免费看| 精品一区二区三区av网在线观看| 午夜免费观看网址| 久久久久久久久免费视频了| 日韩欧美免费精品| 我的老师免费观看完整版| 熟女少妇亚洲综合色aaa.| 亚洲av日韩精品久久久久久密| 又紧又爽又黄一区二区| 亚洲一区高清亚洲精品| 亚洲av成人av| 国产乱人伦免费视频| 精品无人区乱码1区二区| 亚洲,欧美精品.| 亚洲欧美日韩东京热| 欧美日韩福利视频一区二区| 久久久成人免费电影| 亚洲国产欧洲综合997久久,| 精品国产乱子伦一区二区三区| 波多野结衣高清作品| 亚洲在线观看片| 19禁男女啪啪无遮挡网站| 免费看光身美女| 国产美女午夜福利| 宅男免费午夜| 在线十欧美十亚洲十日本专区| 18禁黄网站禁片午夜丰满| 91麻豆精品激情在线观看国产| 老司机深夜福利视频在线观看| 99精品久久久久人妻精品| 九九久久精品国产亚洲av麻豆 | 久久久水蜜桃国产精品网| 国产精品1区2区在线观看.| 夜夜夜夜夜久久久久| 老汉色av国产亚洲站长工具| 国产成+人综合+亚洲专区| 一个人观看的视频www高清免费观看 | 少妇的逼水好多| 久久中文字幕人妻熟女| 激情在线观看视频在线高清| 三级毛片av免费| 88av欧美| 白带黄色成豆腐渣| 桃色一区二区三区在线观看| 99热只有精品国产| 老熟妇乱子伦视频在线观看| 国产视频内射| 亚洲熟妇熟女久久| 99视频精品全部免费 在线 | 天堂动漫精品| 亚洲精品一卡2卡三卡4卡5卡|