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

    基于混合潤滑模型的軋機工作界面表面形貌及油膜特性

    2021-11-25 11:11:58王橋醫(yī)王瀚崔明超過山
    關(guān)鍵詞:壓下率量綱油膜

    王橋醫(yī),王瀚,崔明超,過山

    (杭州電子科技大學(xué)機械工程學(xué)院,浙江杭州,310018)

    軋制[1-2]是目前最重要的金屬成形工藝,約70%的金屬生產(chǎn)過程涉及軋制。冷軋潤滑[3-4]系統(tǒng)在金屬成形過程中起到了至關(guān)重要的作用,它在工件與軋輥間提供一定的潤滑油膜,以避免嚴(yán)重的金屬間接觸從而導(dǎo)致磨損。當(dāng)軋制速度較低時,油膜厚度也較小,帶鋼通常處于邊界潤滑摩擦狀態(tài)。隨著速度的增加,油膜厚度也越大,軋制過程的狀態(tài)由單一的邊界潤滑摩擦狀態(tài)逐步轉(zhuǎn)變?yōu)榛旌蠞櫥琜5-7]摩擦狀態(tài)。

    在混合潤滑條件下,平均潤滑油膜厚度小于均方根(RMS)表面粗糙度的3倍,且界面存在粗糙接觸。界面載荷由粗糙凹谷中的加壓潤滑劑和粗糙凸峰處的粗糙接觸共同承擔(dān)。在建立潤滑油流動模型時,應(yīng)考慮表面粗糙度排列方向,因為粗糙度縱向排列時允許更多的潤滑油通過粗糙通道,而粗糙度橫向排列時有利于垂直于相對運動的粗糙層攜帶更多的潤滑油?;旌蠣顟B(tài)下的摩擦模型必須考慮在這種復(fù)雜的潤滑條件下不同類型的活動過程,這就是混合潤滑條件下潤滑模型[8-10]更難建立的原因。

    由于粗糙表面大多服從高斯分布,CHRISTENSEN[11]用一個多項式函數(shù)來等效高斯分布密度函數(shù),并基于概率統(tǒng)計理論提出了摩擦潤滑粗糙表面隨機模型。之后,學(xué)者大多采用這種模型來研究表面粗糙度對于摩擦潤滑接觸表面性能的影響。MASJEDI 等[12]首先在軋制入口分析計算混合條件下入口油膜厚度,從而引入了混合狀態(tài)的軋制模型。根據(jù)這個模型以及工件和軋輥表面高斯高度分布的信息,可以估計實際接觸面積,還假設(shè)粗糙接觸引起的摩擦力與接觸面積成正比。通過求解塑性方程,計算了工作區(qū)內(nèi)的壓力。通過考慮入口油膜厚度和軋制金屬加工硬化的熱效應(yīng),進(jìn)一步改進(jìn)了先前的模型。該模型最大的缺陷是沒有考慮帶鋼的整體塑性變形對粗糙度變形的影響。KIMURA 等[13]也對軋制過程的混合潤滑進(jìn)行了建模。然而,只在入口區(qū)分析了粗糙度影響,而在體積變形較大的工作區(qū),沒有進(jìn)一步考慮粗糙度影響。LO 等[14]建立了一個簡化了的混合潤滑模型,其中包括工作區(qū)域內(nèi)的粗糙度影響。HAO 等[15]建立了基于多種簡化和假設(shè)下的雷諾方程,該模型描述了軋制過程的3個不同區(qū)域,即入口區(qū)、過渡區(qū)和工作區(qū)。FUJITA 等[16]針對高速軋制的情況下,入口和過渡區(qū)流體壓力逐漸增大,工作區(qū)內(nèi)的壓力梯度對潤滑油流動的影響可以忽略不計。AKBARZADEH等[17]探討了在低速條件下軋制的混合潤滑,即入口區(qū)產(chǎn)生的壓力可以忽略不計。分析表明,工作區(qū)內(nèi)由于楔塞作用可產(chǎn)生較高的流體壓力,而實際經(jīng)驗表明在低速情況下流體動力效應(yīng)仍然存在,這有助于更好地理解軋制過程中的潤滑作用。CHENG 等[18]基于表面粗糙度對軋制界面變形區(qū)壓強、粗糙峰的擠壓變形等方面的影響規(guī)律進(jìn)行了研究,對塑性變形區(qū)摩擦作用機理以及變形區(qū)各種效應(yīng)對潤滑過程影響進(jìn)行了分析,其研究結(jié)論進(jìn)一步加深了人們對金屬塑性加工混合潤滑機理的理解。為此,本文作者利用表面模式參數(shù)法、粗糙表面隨機模型及流體潤滑平均雷諾方程,建立軋制過程混合潤滑數(shù)學(xué)模型。基于混合潤滑模型數(shù)值模擬在各種工況下軋機工作界面表面形貌及油膜特性,分析表面粗糙度橫向排列方式、表面粗糙度各向同性排列方式和表面粗糙度縱向排列方式下的油膜厚度、摩擦因數(shù)、工作界面接觸面積及界面應(yīng)力等在不同軋制速度和壓下率下的變化情況。

    1 界面表面形貌分析

    圖1所示為軋制界面粗糙表面接觸示意圖,其中,τ1為軋輥表面粗糙度平均線以上特定高度,τ2為軋件表面粗糙度平均線以上特定高度,hn為2個未變形表面的平均距離,U為接觸面積。由圖1可見:在混合潤滑狀態(tài)下,界面載荷由粗糙接觸面和表面凹谷中的加壓潤滑劑承擔(dān)。為了模擬該狀態(tài)下的軋制過程,必須計算出粗糙接觸處的平均粗糙度(峰值)壓力p1、在凹谷的平均谷壓(潤滑劑或流體動力)p2和平均界面壓力p。表面形貌對潤滑劑流量計算的影響也變得至關(guān)重要。準(zhǔn)確的表面特征描述對于潤滑劑流動和粗糙變形的建模都很重要。實際的表面是復(fù)雜的,難以用一個參數(shù)來描述表面特征。

    圖1 軋制界面粗糙表面接觸示意圖Fig.1 Rough surface contact diagram of rolling interface

    高斯曲面是工程曲面最常見的模型。高斯曲面的概率密度PG表達(dá)式為

    式中:τ為τ1和τ2的平均值;Rq為均方根(RMS)表面粗糙度。

    式中:x為沿著測量方向的坐標(biāo);g為從平均線測得的表面局部振幅;L為測量區(qū)間的長度。

    為了近似表示高斯曲面,CHRISTENSEN[11]引入了多項式概率密度函數(shù)PC:

    均方根粗糙度描述了輪廓線與平均線的偏差,但沒有給出關(guān)于輪廓線間距或范圍的信息。這可由定義的自相關(guān)函數(shù)?得到

    式中:α為測量自相關(guān)的距離。對于周期曲面,自相關(guān)函數(shù)也是周期的,而對于隨機曲面,常假定其為指數(shù)形式:

    式中:α*稱為自相關(guān)長度。根據(jù)這個定義,自相關(guān)函數(shù)在自相關(guān)長度處的值為1/e。

    對于各向同性曲面,自相關(guān)函數(shù)與采樣方向無關(guān)。對于由沿特定鋪層方向排列的峰和谷組成的各向異性表面,自相關(guān)函數(shù)依賴于樣品方向與鋪層方向的夾角。PEKLENIK[19]提出了表面模式參數(shù)λs的概念:

    式中:α*x和α*y分別為正交方向上的自相關(guān)長度。在以上計算中,α*x和α*y分別是與表面運動方向平行和與表面運動方向垂直的自相關(guān)長度。

    假設(shè)λs是典型粗糙表面的長寬比。當(dāng)λs=0 時,可表示為粗糙度橫向排列方式;當(dāng)λs=1 時,可表示為粗糙度各向同性排列方式;當(dāng)λs=∞時,可表示為粗糙度縱向排列方式。

    對于高斯表面,平均油膜厚度ha為

    當(dāng)軋錕表面和軋件表面沒有接觸時,ha等于hn;當(dāng)表面有接觸時,ha大于hn。在接觸面積U很大的情況下,hn為負(fù)數(shù),而ha總是正數(shù)。平均薄膜厚度也可以轉(zhuǎn)化為量綱一的形式:

    此外,部分接觸面積U與量綱一的表面間隔hn的關(guān)系為

    2 界面混合潤滑數(shù)學(xué)模型

    通常認(rèn)為軋制界面處于混合潤滑狀態(tài)。圖2所示為考慮油膜厚度的軋件塑性變形示意圖。圖2中:W為軋輥表面速度;V為帶鋼速度;σ1為帶鋼入口厚度;σ2為帶鋼出口厚度;V1為帶鋼入口速度;V2為帶鋼出口速度;b為軋輥半徑。

    圖2 軋制變形區(qū)示意圖Fig.2 Sketch map of rolling deformation

    使用平均雷諾方程可計算表面凹谷中潤滑油的壓力。對于穩(wěn)態(tài)、不可壓縮的流體,平均雷諾方程[20]可表示為

    式中:μ為潤滑油黏度;ηx和ηs分別為壓力流量因子和剪切流量因子,它們被用來補償表面粗糙度的影響。

    均方根粗糙度Rq為

    式中:Rqr和Rqs分別為單個軋輥和帶鋼表面的均方根粗糙度。

    不同粗糙度方向流動因子的半經(jīng)驗方程[21]可表示為

    式中:Cx,rx和As為關(guān)于λs的函數(shù)。

    3 基于粗糙度考慮的界面表面載荷分布

    為了計算平均粗糙度壓力p1,應(yīng)考慮工件塑性變形的影響。通過對具有縱向粗糙度的粗糙表面分析研究,發(fā)現(xiàn)材料的大塑性變形降低了粗糙度硬度。量綱一的粗糙度硬度[22]Yr定義為

    式中:c為平面應(yīng)變下工件的剪切強度;P1和P2分別為粗糙度載荷和潤滑載荷。

    文獻(xiàn)[22]將量綱一的粗糙度硬度與量綱一的應(yīng)變速率M聯(lián)系起來:

    其中:

    量綱一的應(yīng)變速率M定義為

    式中:為體應(yīng)變率;l為粗糙度半節(jié)距;v1和v2分別為接觸面的下降速度和凹谷相對于平均表面的平均上升速度。

    潤滑接觸界面處的總界面載荷由粗糙凸峰的接觸和粗糙凹谷中的加壓潤滑劑共同承擔(dān)。載荷平衡方程將總界面載荷P與粗糙度載荷P1、潤滑載荷P2和接觸面積U聯(lián)系起來,表示為[23-24]

    平均界面壓力用量綱一的形式表示為

    在混合潤滑模型中,壓力變量p1,p2和p可通過式(29)進(jìn)行相互轉(zhuǎn)換。

    4 界面摩擦模型

    圖3所示為假定接觸面為類似于鋸齒面情形的示意圖。在混合潤滑狀態(tài)下,流體動力摩擦應(yīng)力由接觸表面凹谷中的潤滑劑傳遞,而邊界摩擦應(yīng)力由接觸表面凸峰的黏附或剪切傳遞。

    圖3 鋸齒表面示意圖Fig.3 Sketch map of sawtooth surface

    對于給定的滑動速度,局部剪切應(yīng)力大于臨界值cl時,流體表現(xiàn)為可塑性;局部剪切應(yīng)力小于臨界值cl時,潤滑油表現(xiàn)為黏稠性,所以,存在臨界剪切厚度[23-24]ht:

    當(dāng)局部膜厚小于ht時,潤滑油表現(xiàn)為可塑性;而當(dāng)局部膜厚大于ht時,潤滑油表現(xiàn)為黏稠性。

    通常假定粗糙接觸處的界面摩擦應(yīng)力γa為

    式中:k為附著系數(shù)。為了方便起見,假設(shè)潤滑油的剪切強度等于邊界膜的屈服剪切強度,即

    所以,

    平均摩擦應(yīng)力γ可表示為

    式中:Up為所給出的部分塑性剪切面積。

    式中:xt為臨界剪切長度,對應(yīng)于臨界剪切厚度ht。對于鋸齒表面的情況,有

    而微凸體高度d一般與未變形工件的表面均方根粗糙度有關(guān):

    利用式(36)和式(37),潤滑油的量綱一的臨界剪切厚度Ht為

    塑性剪切面積Up為

    計算潤滑油平均黏滯摩擦應(yīng)力γb為

    當(dāng)

    時,定義摩擦因數(shù)m為

    將式(32),(40)和式(41)代入式(35),式(44)可化為

    應(yīng)用鋸齒曲面的幾何關(guān)系:

    式(45)可改寫為

    5 軋制界面工作區(qū)分析

    軋制變形過程的幾何形狀如圖3所示。軋輥和帶鋼之間的界面被劃分為幾個區(qū)域進(jìn)行分析。在入口區(qū),認(rèn)為帶鋼是剛性的,隨著流體動力壓力增大,粗糙接觸所承受的載荷也增大,使得總界面應(yīng)力達(dá)到入口區(qū)與工作區(qū)的交界處工件的屈服應(yīng)力,帶鋼在工作區(qū)內(nèi)發(fā)生大塑性變形。圖1從接觸面微觀角度描述了載荷分布情況。

    假設(shè)軋制過程是等溫的,潤滑劑黏度μ隨潤滑劑壓力p2的變化由Barus公式給出:

    式中:λ為黏度的壓力系數(shù)。

    量綱一的壓力黏度系數(shù)G、量綱一的潤滑劑黏度E、壓下率R、量綱一的速度S、入口區(qū)量綱一的位置I、量綱一的局部帶材厚度Θ和出口速度比Z定義如下:

    式中:b為軋輥半徑;x1為軋輥咬入長度;x′為進(jìn)口帶鋼與軋輥表面的虛擬交點距離;σ1,σ2和σ分別為進(jìn)口、出口和任意局部位置的帶鋼厚度。量綱一的平均雷諾方程[28]可以寫成

    其中:C為與潤滑油物理特性相關(guān)的常數(shù)。

    6 數(shù)值仿真與結(jié)果分析

    為了便于討論基于混合潤滑模型的軋機工作界面表面形貌及油膜特性,仿真分析時將軋輥和帶鋼的粗糙度完全等效到帶鋼上。采用數(shù)學(xué)工具軟件Matlab 軟件對上述方程進(jìn)行迭代求解,數(shù)值模擬在各種工況下混合潤滑模型中軋機工作界面表面形貌及油膜特性。表1所示為數(shù)值仿真中使用的主要參數(shù),仿真結(jié)果如圖4~10 所示。帶鋼表面粗糙度采用3種排列方式進(jìn)行仿真:表面粗糙度橫向排列方式(λs=1/6),表面粗糙度各向同性排列方式(λs=1)和表面粗糙度縱向排列方式(λs=6)。

    表1 數(shù)值仿真中主要參數(shù)Table 1 Main parameters in numerical simulation

    圖4所示為量綱一的速度S=10-4時,摩擦因數(shù)m的分布情況。由于混合潤滑狀態(tài)下摩擦使部分粗糙峰接觸,因此,m的分布形式反映了工作界面接觸面積U的相應(yīng)分布。在粗糙度各種排列情形中,粗糙度縱向排列的摩擦力最大,粗糙度橫向排列的摩擦力最小。摩擦應(yīng)力在中性點附近有明顯的突變。因此,摩擦有助于在入口側(cè)拉進(jìn)帶鋼,并阻礙其在出口側(cè)的運動。當(dāng)粗糙度橫向排列占主導(dǎo)地位時,中性點更靠近出口。

    圖4 不同壓下率和粗糙度方向下的摩擦因數(shù)Fig.4 Friction coefficient under different reductions and roughness directions

    圖5所示為固定壓下率R為0.3 時,量綱一的速度S相對于摩擦因數(shù)m的變化情況。較高的速度往往會降低摩擦水平,并將中性點向出口邊緣移動。從圖5可知:隨著軋制速度增加,中性點向出口一側(cè)移動,因而,粗糙度橫向排列時在高壓下率或高速軋制時會出現(xiàn)打滑的趨勢。

    圖5 不同速度和粗糙度方向下的摩擦因數(shù)Fig.5 Friction coefficient under different roughness directions and speeds

    圖6所示為量綱一的油膜厚度Ha隨量綱一的工作區(qū)位置X在不同速度和粗糙度排列方向組合條件下的變化,比較壓下率R為0.6 時不同速度和粗糙度排列方向?qū)τ湍ず穸鹊挠绊?。從圖6可知:油膜厚度隨著帶鋼從入口到出口的過程而減小。粗糙度橫向排列油膜厚度最大,粗糙度縱向排列油膜厚度最小,粗糙度各向同性排列油膜厚度則處在上面2種情形之間。當(dāng)粗糙度縱向排列、低速情形時(λs=6,S=10-5),油膜膜厚在入口處迅速減小,出口處趨于零。量綱一的速度變化對橫向排列下油膜厚度的影響相對較小。

    圖6 不同速度和粗糙度方向下的量綱一的油膜厚度Fig.6 Film thickness under different roughness directions and speeds

    圖7所示為壓下率R為0.6 時量綱一的平均油膜厚度Ha隨量綱一的速度S的變化。從圖7可知:在各種粗糙度排列方向分布中,縱向排列的油膜厚度最低,橫向排列的油膜厚度最高。當(dāng)S很?。ㄐ∮?0-4)時,Ha基本上是常數(shù)。當(dāng)S較高時,混合潤滑狀態(tài)下的Ha隨S增加而單調(diào)增加,最終所有曲線在對應(yīng)厚膜狀態(tài)S的最高處重合。

    圖7 不同粗糙度方向下油膜厚度隨速度的變化Fig.7 Variation of film thickness with speed under different roughness directions

    圖8所示為2 種典型的量綱一的速度S=10-5和S=10-2(分別對應(yīng)低速和高速)條件下量綱一的油膜厚度Ha隨壓下率R的變化情況。從圖8可知:油膜厚度隨壓下率的增大而減小。粗糙度橫向排列油膜厚度最大,粗糙度縱向排列油膜厚度最小,粗糙度各向同性排列則處在上面2種條件之間。值得注意的是,在相同的壓下率下,由于量綱一的速度從10-5增加到10-2,粗糙度橫向排列的油膜厚度增加不明顯,粗糙度縱向排列的油膜厚度明顯增加。在低速條件下(λs=6,S=10-5),系統(tǒng)處于邊界潤滑狀態(tài)。在高速條件下(λs=6,S=10-2),油膜厚度從0.04 左右跳躍到0.20 左右。在高速(λs=1/6,S=10-2)和低速(λs=1/6,S=10-5)粗糙度橫向排列條件下,薄膜厚度差異不明顯。因此,在量綱一的速度S較低時,可能出現(xiàn)以下2 種狀態(tài):一種是所有界面載荷都由粗糙接觸承擔(dān)的真正邊界潤滑狀態(tài),另一種是界面載荷由粗糙接觸和表面凹谷中的加壓潤滑劑共同承擔(dān)的“流體靜力”混合潤滑狀態(tài)。粗糙度橫向排列的油膜厚度分布對量綱一的速度的變化不敏感,而粗糙度各向同性和粗糙度縱向排列的油膜厚度隨量綱一的速度S從10-5增加到10-2而顯著增加。

    圖8 不同粗糙度方向下油膜厚度隨壓下率的變化Fig.8 Variation of film thickness with reduction under different roughness directions

    圖9所示為壓下率R為0.6 以及粗糙度縱向排列(λs=6)和粗糙度橫向排列(λs=1/6)時的量綱一的界面總載荷P和潤滑載荷P2分布。在量綱一的速度S為10-3時,量綱一的工作區(qū)位置X=0 和X=1分別對應(yīng)工作區(qū)的入口和出口邊緣。在粗糙度橫向排列時,潤滑載荷變化不明顯,它們在中性點處有一個小峰值,在粗糙度縱向排列時,量綱一的界面總載荷P的變化比較大,并且在中性點處有一個明顯的峰值,但這離出口更遠(yuǎn)。這些趨勢都反映出與粗糙度縱向排列相比,粗糙度橫向排列的摩擦要小得多。

    粗糙度縱向分布的量綱一的潤滑載荷P2的分布如圖9所示。它們大致呈拋物線形,在入口和出口處為零,具有粗糙度橫向分布的量綱一的潤滑載荷P2的分布基本上是平的,與P的相應(yīng)分布曲線大致平行。

    圖9 不同粗糙度方向和壓下率下的應(yīng)力變化Fig.9 Stress changes under different roughness directions and reductions

    圖10所示為在量綱一的速度S=10-3時,對于不同的壓下率和粗糙度排列方式,工作界面接觸面積U隨量綱一的工作區(qū)位置X的變化。從圖10可知:工作界面接觸面積U從入口到出口都是單調(diào)增加的。粗糙度縱向排列情況下工作界面接觸面積U遠(yuǎn)大于粗糙度橫向排列情況下工作界面接觸面積U。壓下率對軋制界面接觸面積U的作用不明顯。

    圖10 不同粗糙度方向和壓下率下的工作界面接觸面積Fig.10 Contact area of working interface under different reductions and roughness directions

    7 結(jié)論

    1)在所有粗糙度分布情形下,油膜厚度隨壓下率的增大而減??;粗糙度橫向排列時工作界面油膜厚度最大,粗糙度縱向排列時工作界面油膜厚度最小,粗糙度各向同性排列時工作界面油膜厚度則介于上述2種情形之間。

    2)當(dāng)量綱一的速度很?。ㄐ∮?0-4)時,油膜厚度基本上是常數(shù),在速度較高時,油膜厚度隨速度單調(diào)增加;當(dāng)量綱一的速度從10-5增加到10-2,粗糙度橫向排列時,工作界面的油膜厚度變化隨量綱一的速度的變化不敏感,而粗糙度各向同性和粗糙度縱向排列的油膜厚度隨量綱一的速度從10-5增加到10-2而顯著增加。

    3)對于不同的壓下率和粗糙度排列方式,工作界面接觸面積從入口到出口都是單調(diào)增加的。粗糙度縱向排列情況下工作界面接觸面積遠(yuǎn)比粗糙度橫向排列情況下的大,粗糙度各向同性時的工作界面接觸面積處于兩者之間。

    4)當(dāng)壓下率為定值時,不同速度和粗糙度排列方向下,油膜厚度隨著帶鋼從入口到出口的過程而減小。對于低速條件下的粗糙度縱向排列,油膜膜厚在入口處迅速減小,出口處趨于零。

    5)在粗糙度各種排列情形中,粗糙度縱向排列的摩擦力最大,粗糙度橫向排列的摩擦力最小。摩擦應(yīng)力在軋制過程中性點附件有明顯的突變。當(dāng)粗糙度橫向排列占主導(dǎo)地位時,中性點更靠近出口。

    猜你喜歡
    壓下率量綱油膜
    壓下率對冷軋及退火純鈦板材織構(gòu)的影響
    量綱分析在熱力學(xué)統(tǒng)計物理中的教學(xué)應(yīng)用*
    物理通報(2022年6期)2022-05-27 09:43:14
    長城油膜軸承油在高速棒材生產(chǎn)線的應(yīng)用
    石油商技(2021年1期)2021-03-29 02:36:04
    淺談量綱法推導(dǎo)物理公式的優(yōu)勢
    ——以勻加速直線運動公式為例
    熱軋7075/AZ31B復(fù)合板的顯微組織及結(jié)合性能
    科技論文編輯加工中的量綱問題
    大型數(shù)控立式磨床靜壓轉(zhuǎn)臺油膜熱特性仿真及其實驗分析
    冷軋軋機油膜軸承系統(tǒng)故障分析與對策
    新疆鋼鐵(2016年4期)2016-05-25 00:37:38
    冷軋壓下率對含硼和鈮細(xì)晶高強IF 鋼組織與性能的影響
    上海金屬(2014年5期)2014-12-20 07:58:44
    冷軋壓下率對退火含鈮鈦IF鋼性能及織構(gòu)的影響
    国产亚洲一区二区精品| 欧美日韩国产mv在线观看视频| 亚洲第一青青草原| 中文字幕人妻丝袜制服| 国产精品 国内视频| av视频免费观看在线观看| 日韩中文字幕视频在线看片| 久久久精品区二区三区| 久久午夜综合久久蜜桃| av网站在线播放免费| 日韩视频一区二区在线观看| 久久久久久人人人人人| 久久午夜综合久久蜜桃| 最黄视频免费看| 日韩一卡2卡3卡4卡2021年| 搡老乐熟女国产| 人人妻人人添人人爽欧美一区卜| 又大又爽又粗| 波多野结衣av一区二区av| 午夜福利影视在线免费观看| 免费在线观看黄色视频的| 这个男人来自地球电影免费观看| 欧美黑人欧美精品刺激| 久久 成人 亚洲| 女警被强在线播放| 欧美xxⅹ黑人| 亚洲精品国产精品久久久不卡| 大香蕉久久网| 飞空精品影院首页| 国产av国产精品国产| 国产伦人伦偷精品视频| 国产精品国产av在线观看| 99国产极品粉嫩在线观看| 国产男女内射视频| 国产有黄有色有爽视频| 精品亚洲乱码少妇综合久久| 在线观看一区二区三区激情| 亚洲精品一区蜜桃| 国产免费视频播放在线视频| 亚洲人成电影观看| 国产一卡二卡三卡精品| 免费黄频网站在线观看国产| 久久av网站| 韩国精品一区二区三区| 天堂8中文在线网| 精品熟女少妇八av免费久了| 午夜福利在线免费观看网站| 精品乱码久久久久久99久播| 一级片免费观看大全| av网站免费在线观看视频| 又紧又爽又黄一区二区| 国产男人的电影天堂91| 国产一区二区三区在线臀色熟女 | 高清欧美精品videossex| 成人黄色视频免费在线看| 1024香蕉在线观看| 亚洲精品国产av蜜桃| 成年av动漫网址| 搡老乐熟女国产| 亚洲欧美精品自产自拍| 久久人人爽人人片av| 国产av国产精品国产| 国产高清videossex| 人人妻人人添人人爽欧美一区卜| av国产精品久久久久影院| 淫妇啪啪啪对白视频 | 日韩视频一区二区在线观看| tube8黄色片| 少妇 在线观看| 精品少妇久久久久久888优播| 久久人妻福利社区极品人妻图片| 亚洲全国av大片| 亚洲va日本ⅴa欧美va伊人久久 | 少妇猛男粗大的猛烈进出视频| 狠狠狠狠99中文字幕| 99国产极品粉嫩在线观看| av免费在线观看网站| 天天躁日日躁夜夜躁夜夜| 久久久精品区二区三区| 亚洲精品美女久久久久99蜜臀| 欧美精品人与动牲交sv欧美| 天天操日日干夜夜撸| 蜜桃国产av成人99| av一本久久久久| 黄片大片在线免费观看| 交换朋友夫妻互换小说| 中文字幕色久视频| 日韩中文字幕视频在线看片| 中文字幕另类日韩欧美亚洲嫩草| 国产精品麻豆人妻色哟哟久久| 男人爽女人下面视频在线观看| 日韩一区二区三区影片| 久久狼人影院| 久久人人97超碰香蕉20202| 人人妻人人添人人爽欧美一区卜| 亚洲伊人久久精品综合| 成年人黄色毛片网站| 国产黄频视频在线观看| 国产区一区二久久| 母亲3免费完整高清在线观看| 国产在线免费精品| videos熟女内射| 久久国产精品影院| 亚洲男人天堂网一区| 午夜福利在线免费观看网站| 精品一区二区三卡| 无限看片的www在线观看| 大码成人一级视频| 欧美日韩视频精品一区| 日本猛色少妇xxxxx猛交久久| 伊人亚洲综合成人网| 亚洲av成人不卡在线观看播放网 | 国产亚洲精品第一综合不卡| 极品少妇高潮喷水抽搐| 亚洲第一青青草原| 精品人妻在线不人妻| 精品国产国语对白av| 午夜老司机福利片| 免费观看a级毛片全部| 亚洲精品久久久久久婷婷小说| 亚洲五月婷婷丁香| 亚洲中文av在线| 成年动漫av网址| 日韩一区二区三区影片| 狂野欧美激情性bbbbbb| 亚洲第一av免费看| 超碰97精品在线观看| 欧美久久黑人一区二区| 香蕉国产在线看| av超薄肉色丝袜交足视频| 老鸭窝网址在线观看| 色综合欧美亚洲国产小说| 成人影院久久| 亚洲五月婷婷丁香| 国产伦人伦偷精品视频| 亚洲成人免费av在线播放| 国产欧美日韩一区二区三区在线| 80岁老熟妇乱子伦牲交| 美女高潮喷水抽搐中文字幕| 久久久久国产精品人妻一区二区| 99久久99久久久精品蜜桃| 丰满人妻熟妇乱又伦精品不卡| 久久国产精品人妻蜜桃| www.自偷自拍.com| 亚洲成人国产一区在线观看| 欧美av亚洲av综合av国产av| 在线观看免费高清a一片| 国产男人的电影天堂91| 久久久精品国产亚洲av高清涩受| 免费在线观看日本一区| 两性午夜刺激爽爽歪歪视频在线观看 | 汤姆久久久久久久影院中文字幕| 热99久久久久精品小说推荐| 少妇精品久久久久久久| 久久久久久久大尺度免费视频| 国产精品99久久99久久久不卡| 成年人午夜在线观看视频| 动漫黄色视频在线观看| 国产成人a∨麻豆精品| 老司机午夜福利在线观看视频 | 久久久久久久国产电影| 老熟女久久久| 美女大奶头黄色视频| 亚洲欧洲精品一区二区精品久久久| 窝窝影院91人妻| 中国国产av一级| 精品久久蜜臀av无| 狂野欧美激情性xxxx| 嫁个100分男人电影在线观看| 亚洲熟女精品中文字幕| 狂野欧美激情性xxxx| 别揉我奶头~嗯~啊~动态视频 | 国产男女内射视频| 在线观看免费午夜福利视频| 在线观看免费日韩欧美大片| 亚洲av片天天在线观看| 亚洲精品中文字幕在线视频| 50天的宝宝边吃奶边哭怎么回事| 欧美xxⅹ黑人| 捣出白浆h1v1| 国产在线视频一区二区| 两人在一起打扑克的视频| svipshipincom国产片| 日韩大片免费观看网站| tocl精华| 欧美成狂野欧美在线观看| 国产色视频综合| 少妇的丰满在线观看| 午夜成年电影在线免费观看| 五月开心婷婷网| 伊人久久大香线蕉亚洲五| a级毛片黄视频| 色婷婷久久久亚洲欧美| 国产精品1区2区在线观看. | 少妇猛男粗大的猛烈进出视频| 69精品国产乱码久久久| 国产精品一区二区在线不卡| 国产精品久久久久成人av| 久久青草综合色| 首页视频小说图片口味搜索| 一级,二级,三级黄色视频| 国产一卡二卡三卡精品| 亚洲国产精品一区二区三区在线| 国产精品自产拍在线观看55亚洲 | 午夜免费成人在线视频| 成年av动漫网址| 最近中文字幕2019免费版| 午夜精品久久久久久毛片777| 欧美精品亚洲一区二区| 天堂8中文在线网| 久久久久久亚洲精品国产蜜桃av| 久久香蕉激情| 国产成人免费无遮挡视频| 国产熟女午夜一区二区三区| 亚洲精品国产精品久久久不卡| 国产深夜福利视频在线观看| 亚洲全国av大片| 久久亚洲国产成人精品v| 欧美在线一区亚洲| 男女床上黄色一级片免费看| 精品一品国产午夜福利视频| netflix在线观看网站| 久久久久精品国产欧美久久久 | 免费在线观看影片大全网站| 亚洲三区欧美一区| 亚洲五月色婷婷综合| 国产精品久久久久久精品古装| xxxhd国产人妻xxx| 亚洲精品国产色婷婷电影| 美女扒开内裤让男人捅视频| 国产成人欧美在线观看 | 国产av一区二区精品久久| 成年人免费黄色播放视频| 69av精品久久久久久 | av电影中文网址| 久久国产精品影院| 少妇 在线观看| 又黄又粗又硬又大视频| 动漫黄色视频在线观看| 91成人精品电影| 亚洲成人免费av在线播放| 岛国毛片在线播放| 免费观看a级毛片全部| 女人爽到高潮嗷嗷叫在线视频| 精品人妻在线不人妻| 中文欧美无线码| 欧美大码av| 热99国产精品久久久久久7| 久久ye,这里只有精品| 2018国产大陆天天弄谢| 老司机福利观看| 男女国产视频网站| av不卡在线播放| 考比视频在线观看| 国产男人的电影天堂91| 午夜福利视频在线观看免费| 久9热在线精品视频| 黑人操中国人逼视频| 欧美日韩黄片免| 狠狠精品人妻久久久久久综合| 午夜精品国产一区二区电影| 天天躁狠狠躁夜夜躁狠狠躁| 丝瓜视频免费看黄片| 久久久精品免费免费高清| 久久久久久久久久久久大奶| 三级毛片av免费| 国产男女内射视频| 久久午夜综合久久蜜桃| 国精品久久久久久国模美| netflix在线观看网站| 亚洲欧美精品综合一区二区三区| 日韩欧美一区视频在线观看| 狠狠婷婷综合久久久久久88av| 高潮久久久久久久久久久不卡| 精品视频人人做人人爽| 国产男人的电影天堂91| 人人妻,人人澡人人爽秒播| 人人妻人人澡人人爽人人夜夜| 成在线人永久免费视频| 男女免费视频国产| 大香蕉久久成人网| 欧美日韩中文字幕国产精品一区二区三区 | 久久人妻福利社区极品人妻图片| 亚洲av国产av综合av卡| 欧美精品一区二区大全| 成人影院久久| 天天躁狠狠躁夜夜躁狠狠躁| av片东京热男人的天堂| 嫩草影视91久久| 最新的欧美精品一区二区| 看免费av毛片| www.自偷自拍.com| 99re6热这里在线精品视频| 老司机影院毛片| 99精品久久久久人妻精品| 国产男女内射视频| av在线播放精品| 亚洲精品国产区一区二| 天天添夜夜摸| 国产一区二区 视频在线| 啦啦啦啦在线视频资源| 搡老岳熟女国产| 欧美黑人欧美精品刺激| 色精品久久人妻99蜜桃| 欧美人与性动交α欧美精品济南到| 久久久国产一区二区| 少妇 在线观看| 性色av乱码一区二区三区2| 中文字幕最新亚洲高清| 久久精品亚洲熟妇少妇任你| 国产深夜福利视频在线观看| 国产野战对白在线观看| 性高湖久久久久久久久免费观看| videosex国产| 青草久久国产| av天堂久久9| 午夜91福利影院| 9热在线视频观看99| 亚洲中文字幕日韩| kizo精华| 国产免费福利视频在线观看| e午夜精品久久久久久久| 丝袜脚勾引网站| 亚洲精品中文字幕一二三四区 | 一级黄色大片毛片| 90打野战视频偷拍视频| 操美女的视频在线观看| 久久久久久久精品精品| 国产精品一区二区精品视频观看| 日韩视频在线欧美| 大陆偷拍与自拍| 国产成人av教育| 欧美乱码精品一区二区三区| 大香蕉久久网| 国产成人免费无遮挡视频| 少妇精品久久久久久久| 交换朋友夫妻互换小说| 在线av久久热| 国产成人免费无遮挡视频| 国产亚洲一区二区精品| 欧美老熟妇乱子伦牲交| 欧美精品人与动牲交sv欧美| 成人三级做爰电影| 男人舔女人的私密视频| 18禁黄网站禁片午夜丰满| 精品国产乱码久久久久久小说| 人人妻人人添人人爽欧美一区卜| 亚洲av成人不卡在线观看播放网 | 国产精品秋霞免费鲁丝片| 美国免费a级毛片| 亚洲伊人色综图| 狠狠狠狠99中文字幕| 国产成人一区二区三区免费视频网站| 久久 成人 亚洲| 超碰97精品在线观看| 亚洲精品美女久久av网站| 国产97色在线日韩免费| 国产成人系列免费观看| 99国产综合亚洲精品| 国产极品粉嫩免费观看在线| 黄色a级毛片大全视频| 欧美日韩av久久| a级片在线免费高清观看视频| 18禁裸乳无遮挡动漫免费视频| 国产精品1区2区在线观看. | 亚洲精品粉嫩美女一区| 久久免费观看电影| 久久久精品国产亚洲av高清涩受| 日韩免费高清中文字幕av| www.av在线官网国产| 免费av中文字幕在线| 一边摸一边做爽爽视频免费| 国产淫语在线视频| 亚洲欧美精品自产自拍| av又黄又爽大尺度在线免费看| 91九色精品人成在线观看| 天天躁夜夜躁狠狠躁躁| 成年美女黄网站色视频大全免费| 欧美日韩亚洲国产一区二区在线观看 | 首页视频小说图片口味搜索| 婷婷丁香在线五月| 国产99久久九九免费精品| 宅男免费午夜| 天堂8中文在线网| 日本精品一区二区三区蜜桃| 丁香六月欧美| 美女中出高潮动态图| 51午夜福利影视在线观看| 两人在一起打扑克的视频| 国产在视频线精品| 久久狼人影院| 叶爱在线成人免费视频播放| 欧美人与性动交α欧美精品济南到| 91精品国产国语对白视频| 80岁老熟妇乱子伦牲交| 国产高清国产精品国产三级| 十八禁网站网址无遮挡| 黄片播放在线免费| 在线观看免费视频网站a站| 女警被强在线播放| 国精品久久久久久国模美| 热re99久久精品国产66热6| kizo精华| 超色免费av| 久久人人爽av亚洲精品天堂| 亚洲精品一区蜜桃| 又大又爽又粗| 少妇 在线观看| 国产成人影院久久av| 在线 av 中文字幕| 五月开心婷婷网| 中文精品一卡2卡3卡4更新| 亚洲av成人一区二区三| 老汉色∧v一级毛片| 丝瓜视频免费看黄片| 中文字幕高清在线视频| 男女无遮挡免费网站观看| 免费人妻精品一区二区三区视频| 黄色a级毛片大全视频| 韩国精品一区二区三区| 欧美xxⅹ黑人| 青青草视频在线视频观看| 91成人精品电影| 精品福利永久在线观看| 国产成人免费无遮挡视频| 欧美黑人欧美精品刺激| 精品高清国产在线一区| videosex国产| 高潮久久久久久久久久久不卡| 欧美在线黄色| 免费高清在线观看视频在线观看| 少妇猛男粗大的猛烈进出视频| 午夜福利视频在线观看免费| 亚洲精品日韩在线中文字幕| 少妇猛男粗大的猛烈进出视频| 亚洲国产看品久久| 久9热在线精品视频| 国产精品免费视频内射| 亚洲av欧美aⅴ国产| 久久久久国产精品人妻一区二区| 久9热在线精品视频| 人人妻,人人澡人人爽秒播| 爱豆传媒免费全集在线观看| 亚洲成人免费av在线播放| 久久香蕉激情| 下体分泌物呈黄色| 久久av网站| 国产国语露脸激情在线看| 搡老熟女国产l中国老女人| 午夜免费观看性视频| 99热国产这里只有精品6| 制服诱惑二区| 午夜视频精品福利| 天天操日日干夜夜撸| 久久精品亚洲熟妇少妇任你| 亚洲国产欧美一区二区综合| 狂野欧美激情性bbbbbb| 成年美女黄网站色视频大全免费| 国产成人精品久久二区二区免费| 久久毛片免费看一区二区三区| 老司机午夜十八禁免费视频| 又紧又爽又黄一区二区| 国产真人三级小视频在线观看| 亚洲人成电影观看| 亚洲欧洲精品一区二区精品久久久| 我的亚洲天堂| 国产深夜福利视频在线观看| 可以免费在线观看a视频的电影网站| 久久精品熟女亚洲av麻豆精品| 久久热在线av| 成在线人永久免费视频| 亚洲欧美精品综合一区二区三区| 别揉我奶头~嗯~啊~动态视频 | 免费人妻精品一区二区三区视频| 少妇 在线观看| 99国产精品一区二区蜜桃av | 国产免费av片在线观看野外av| 如日韩欧美国产精品一区二区三区| 欧美人与性动交α欧美软件| 十八禁网站网址无遮挡| 国产男人的电影天堂91| 婷婷成人精品国产| 国产精品久久久久久人妻精品电影 | 欧美xxⅹ黑人| 制服人妻中文乱码| 99国产极品粉嫩在线观看| 国产免费福利视频在线观看| av国产精品久久久久影院| 亚洲国产精品999| 18禁裸乳无遮挡动漫免费视频| www.自偷自拍.com| 免费观看a级毛片全部| 天天躁夜夜躁狠狠躁躁| 国产男女内射视频| 19禁男女啪啪无遮挡网站| 欧美 亚洲 国产 日韩一| 久久精品成人免费网站| 国产深夜福利视频在线观看| 午夜成年电影在线免费观看| 国产三级黄色录像| 国产精品久久久久久精品电影小说| 亚洲色图 男人天堂 中文字幕| 午夜免费成人在线视频| 国产高清videossex| 久久久国产精品麻豆| 老熟妇乱子伦视频在线观看 | 欧美人与性动交α欧美软件| 免费高清在线观看视频在线观看| 精品少妇一区二区三区视频日本电影| 天天操日日干夜夜撸| 青青草视频在线视频观看| 亚洲精品一区蜜桃| 久久性视频一级片| 狂野欧美激情性bbbbbb| 大香蕉久久网| 91大片在线观看| 青春草视频在线免费观看| 国产成人欧美在线观看 | 亚洲国产中文字幕在线视频| 国产高清视频在线播放一区 | 欧美人与性动交α欧美软件| 国产一区二区激情短视频 | 新久久久久国产一级毛片| 国产成人影院久久av| 欧美精品亚洲一区二区| a级毛片在线看网站| 国产亚洲一区二区精品| 久久久久久亚洲精品国产蜜桃av| 久久av网站| 国产在视频线精品| 青草久久国产| 91成人精品电影| 久久久精品94久久精品| 两个人免费观看高清视频| 黑人猛操日本美女一级片| 美女中出高潮动态图| 十分钟在线观看高清视频www| 久久久久久免费高清国产稀缺| 啪啪无遮挡十八禁网站| 亚洲情色 制服丝袜| 欧美黑人精品巨大| 成人国语在线视频| 极品少妇高潮喷水抽搐| 十分钟在线观看高清视频www| 国产精品香港三级国产av潘金莲| www.熟女人妻精品国产| 少妇被粗大的猛进出69影院| 亚洲av男天堂| 一区二区三区激情视频| 伊人亚洲综合成人网| 在线十欧美十亚洲十日本专区| av有码第一页| 一级毛片精品| 日韩熟女老妇一区二区性免费视频| 久久热在线av| 捣出白浆h1v1| 日日摸夜夜添夜夜添小说| 中文字幕另类日韩欧美亚洲嫩草| 亚洲中文字幕日韩| 精品久久蜜臀av无| 国产欧美日韩精品亚洲av| 免费日韩欧美在线观看| 久久香蕉激情| 日本欧美视频一区| 黑人猛操日本美女一级片| 久热这里只有精品99| 夜夜夜夜夜久久久久| 老熟女久久久| 日韩一卡2卡3卡4卡2021年| 久久久精品区二区三区| 国产福利在线免费观看视频| 曰老女人黄片| 欧美精品一区二区大全| 国产视频一区二区在线看| tocl精华| 国产日韩欧美亚洲二区| 亚洲精品乱久久久久久| h视频一区二区三区| 亚洲avbb在线观看| 亚洲国产看品久久| 亚洲第一欧美日韩一区二区三区 | 欧美另类一区| 久久狼人影院| 精品一区二区三区四区五区乱码| 女人高潮潮喷娇喘18禁视频| av又黄又爽大尺度在线免费看| 国产欧美亚洲国产| 九色亚洲精品在线播放| 老司机午夜十八禁免费视频| 男女床上黄色一级片免费看| 老熟妇乱子伦视频在线观看 | 热re99久久国产66热| 大陆偷拍与自拍| 9191精品国产免费久久| 免费高清在线观看日韩| 热re99久久国产66热| 国产一区二区三区在线臀色熟女 | 精品卡一卡二卡四卡免费| 99精国产麻豆久久婷婷| 热99国产精品久久久久久7| 欧美xxⅹ黑人| 老熟妇乱子伦视频在线观看 | 菩萨蛮人人尽说江南好唐韦庄| 久久精品国产a三级三级三级| 日本wwww免费看| 午夜影院在线不卡| av一本久久久久| 亚洲第一青青草原| 窝窝影院91人妻| 精品久久久久久久毛片微露脸 | 一个人免费在线观看的高清视频 | 久久免费观看电影| av在线老鸭窝| 我的亚洲天堂| 亚洲av电影在线进入|