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

    多角度彎曲管磨粒流加工數(shù)值模擬研究

    2021-06-23 09:36:12尹洪超翟鎮(zhèn)德
    中國機(jī)械工程 2021年11期
    關(guān)鍵詞:磨料磨粒湍流

    尹洪超 劉 宵 翟鎮(zhèn)德 穆 林

    大連理工大學(xué)能源與動力學(xué)院,大連,116024

    0 引言

    隨著工業(yè)的進(jìn)步與發(fā)展,人們對機(jī)械零件的加工制造有了更高要求,其中表面粗糙度是零件加工過程的一項重要指標(biāo)。降低零件的表面粗糙度有利于減小摩擦、降低磨損、提高機(jī)械工作效率等。傳統(tǒng)的減小表面粗糙度值的方法主要有機(jī)械法和化學(xué)法。機(jī)械法是通過砂紙、砂輪等硬度較高的固體打磨被加工材料的表面[1],從而達(dá)到表面拋光的目的。化學(xué)法是通過化學(xué)試劑與被加工零件表面發(fā)生化學(xué)反應(yīng),使表面不平整的凸出的部分優(yōu)先溶解,從而達(dá)到降低表面粗糙度的目的。無論是機(jī)械法還是化學(xué)法,都存在加工工藝的局限性。對于復(fù)雜曲面、狹窄空間等零件表面,采用傳統(tǒng)的加工工藝很難達(dá)到給定的表面粗糙度要求,而產(chǎn)生于20世紀(jì)80年代的磨粒流加工工藝可有效解決此類表面的拋光加工問題。

    磨粒流加工本質(zhì)上是一個固液兩相流問題,在許多工業(yè)過程中都涉及到兩相流引起的磨損問題。如何雅玲等[2]研究了鍋爐內(nèi)高溫?zé)煔鈱Q熱器壁面的磨損現(xiàn)象;王坤等[3]、李昳[4]、殷華鋒等[5]分別研究了固液兩相流對異徑管、離心泵、攪拌罐的侵蝕磨損問題。研究管道內(nèi)固液兩相流的流動特性,對于預(yù)防和減少壁面磨損、粗糙壁面拋光等有重要意義。

    磨粒流的加工特性與流速、壓力、磨粒濃度等諸多因素有關(guān)[6-8]。李俊燁等[9]研究了共軌管、噴嘴小孔、整體葉輪等不同零件的磨粒流加工過程特性;ZHANG等[10]研究了夾具表面形狀對磨粒流加工的影響;CHENG等[11]研究了磨粒流加工的磨削機(jī)理,提出了改進(jìn)的Preston方程??傮w來說,磨粒流加工過程的主要影響因素包括流體特性、磨粒性質(zhì)、力學(xué)特性、流道及夾具特性。

    細(xì)長圓管作為輸送高速流體、高壓氣體等介質(zhì)的通道,在工業(yè)生產(chǎn)中應(yīng)用廣泛。本文以某船用高壓氧氣管道為研究對象,對不同彎曲角度的細(xì)長圓管內(nèi)表面拋光過程進(jìn)行數(shù)值模擬,重點分析了流體特性和磨粒性質(zhì)兩類因素對磨粒流加工過程的影響,從而為降低高壓氧氣輸送時的壓力和速度損耗提供理論支持。

    1 固液兩相流的數(shù)學(xué)模型

    計時鳴等[12]以基于拉格朗日方法的VOF模型進(jìn)行了低濃度磨粒流的數(shù)值模擬。為研究較高濃度的磨粒流加工特性,本文選擇歐拉方法,將流體相與固體相均視作連續(xù)相進(jìn)行模擬計算,假設(shè)流體為不可壓縮流體,連續(xù)方程和動量守恒方程分別為

    (1)

    式中,ρ為流體密度;t為時間;u為速度矢量;i為規(guī)定方向,j為垂直于i的方向,ij表示i方向與j方向的相互作用;xi為i方向的坐標(biāo);xj為j方向的坐標(biāo);p為壓力;τ為剪切應(yīng)力張量;gi為i方向的重力加速度;Fi為其余的體積力。

    湍流的求解方程分為湍流動能方程和湍流耗散率方程:

    式中,k為湍流動能;μ為流體黏度;μt為湍動黏度;σk為湍流動能的普朗特數(shù);σε為湍流耗散率的普朗特數(shù);Gk為由平均速度梯度產(chǎn)生的湍流動能;Gb為浮力產(chǎn)生的湍流動能;ε為耗散率;YM為可壓縮湍流脈動膨脹對總耗散率的影響;C1ε、C2ε、C3ε為經(jīng)驗系數(shù),本文取C1ε=1.44,C2ε=1.92,C3ε=0.09。

    湍動黏度μt和平均速度梯度產(chǎn)生的湍流動能Gk由下式計算:

    式中,Cμ為經(jīng)驗常數(shù),本文取Cμ=0.09。

    能量輸運(yùn)方程用來求解流體及固體域的溫度場,形式如下:

    式中,E為單位質(zhì)量的能量;keff為有效傳熱系數(shù),keff=k0+kt,k0為介質(zhì)傳熱系數(shù),kt為根據(jù)湍流模型確定的湍流傳熱系數(shù);T為溫度;h為焓;Jm為物質(zhì)m的擴(kuò)散通量;τeff為有效剪切應(yīng)力張量;Sh為熱源。

    等號右邊四項分別表示能量輸運(yùn)的傳導(dǎo)項、擴(kuò)散項、黏性耗散項以及熱源項。

    單位質(zhì)量的能量E和管壁的固體域?qū)岱匠谭謩e為

    式(9)中,等號左邊第一項為非穩(wěn)態(tài)項,第二項表示固體旋轉(zhuǎn)或平移運(yùn)動的熱傳輸,右邊兩項分別為固體導(dǎo)熱項和熱源項。

    2 物理模型及模擬設(shè)置

    2.1 物理模型

    本文的模擬對象為6種不同彎曲角度的彎管,彎曲角度α分別為15°、30°、45°、60°、75°、90°,管長1 m,外徑D=18 mm,內(nèi)徑d=9 mm。圖1僅展示了α為30°、60°、90°的管型。流體的第一相為水,第二相為SiC固體顆粒。選擇的模擬軟件為FLUENT18.2,對幾何模型劃分網(wǎng)格如圖2、圖3所示,分別為流體域和固體域,并對流體域及彎管肘部進(jìn)行網(wǎng)格加密,網(wǎng)格數(shù)約為204萬。

    圖1 不同彎曲角度的彎管物理模型

    圖2 肘部網(wǎng)格劃分

    圖3 入口處網(wǎng)格劃分

    2.2 模擬參數(shù)設(shè)置

    首先估算雷諾數(shù),Re=ρvdh/μ,ρ取998.2 kg/m3,v為流體的速度,取v=10 m/s,dh為圓管的當(dāng)量直徑,取dh=9 mm,μ取1.003×10-3kg/(m·s),計算得Re=8.9×104。當(dāng)速度增大時,雷諾數(shù)會隨之增大,該雷諾數(shù)為本文模擬工況中的最小雷諾數(shù),故判斷管內(nèi)流動狀態(tài)為湍流狀態(tài)。選擇標(biāo)準(zhǔn)k-ε湍流模型,并開啟能量方程。多相流方程選擇了Mixture兩相流模型,這是由于模擬中所用的兩相流流體的第二相磨料濃度為10%~30%,且模擬過程未發(fā)生相變,Mixture兩相流模型是基于歐拉方法的流動模型,更適合高濃度的兩相流模擬,而VOF模型、DPM模型均不符合計算條件,模擬結(jié)果的準(zhǔn)確性較差。

    穩(wěn)態(tài)模擬的入口邊界條件為速度入口(velocity inlet),主相為水,第二相為SiC顆粒,磨料濃度用顆粒體積分?jǐn)?shù)η表示,磨料參數(shù)見表1。出口邊界條件為自由流出口(outflow)。管壁的物理參數(shù)取值見表2,內(nèi)壁開啟黏性熱計算,用以計算流體與壁面摩擦的產(chǎn)熱。外壁面設(shè)置為對流換熱表面,表面?zhèn)鳠嵯禂?shù)取20 W/(m2·K),用來模擬暴露在空氣中的自然對流換熱。

    表1 磨粒流參數(shù)

    表2 管壁參數(shù)

    為研究磨粒流加工過程中的熱量積聚情況,本文選取α= 45°的管型,在入口流速vin=30 m/s的工況下進(jìn)行非穩(wěn)態(tài)的模擬研究,時間步長為10 s,模擬時間為1 h。

    2.3 無關(guān)性及準(zhǔn)確性驗證

    使用有限元方法進(jìn)行數(shù)值模擬計算時,網(wǎng)格數(shù)量會影響計算的精度,這里對網(wǎng)格無關(guān)性進(jìn)行了驗證,結(jié)果如表3所示,表明計算偏差在可接受范圍內(nèi),可以認(rèn)為計算結(jié)果與網(wǎng)格數(shù)量無關(guān)。此外,F(xiàn)LUENT軟件中的不同湍流模型是對實際湍流作用的簡化,因此驗證本文使用模型的準(zhǔn)確性是必要的。這里選用李俊燁等[9]對噴嘴小孔的模擬及實驗結(jié)果作為對比,驗證結(jié)果見圖4,可以認(rèn)為本文采用的湍流模型是準(zhǔn)確的。

    表3 網(wǎng)格無關(guān)性驗證

    圖4 模型準(zhǔn)確性驗證

    3 結(jié)果分析

    3.1 入口流速對局部壓差的影響

    Preston方程[13]用于描述工件表面材料去除量:

    (10)

    其中,kP為Preston系數(shù),與加工過程的諸多因素有關(guān),切削量(即材料去除量)Δz與加工時間t、速度v、壓力p成正比。設(shè)置入口流速為定值,以p的分布反映切削量的大小。

    分別對局部壓力和進(jìn)出口壓力進(jìn)行分析,這里局部壓力是指管子的彎曲部位內(nèi)壁的壓力,見圖5。將內(nèi)壁曲率小處稱為內(nèi)側(cè)(此處壓力記為pin),曲率大處稱為外側(cè)(此處壓力記為pout),這里截取了彎曲部位前后各25 mm的長度和彎曲部位,在彎曲部位使用流道長度描述坐標(biāo)位置是不準(zhǔn)確的,為了保證管壁的內(nèi)側(cè)與外側(cè)一一對應(yīng),圖5的橫坐標(biāo)在彎曲部位使用了彎角β來描述流道長度,圖1中標(biāo)注了在α=60°時,彎角β的角度劃分。由伯努利原理可知,流體流速快的部位壓力小,管子內(nèi)表面靠近內(nèi)側(cè)的流體流速大,而靠近外側(cè)的流體流速小,因此在內(nèi)壁內(nèi)側(cè)產(chǎn)生了收縮效應(yīng),導(dǎo)致內(nèi)、外側(cè)出現(xiàn)了較大的壓力差。在實際的磨粒流加工過程中,為保證壁面粗糙度盡可能均勻,需要避免較大的壓力差出現(xiàn)。

    (a) α=15°

    根據(jù)Preston方程,材料的切削量與pvin的整體作用效果有關(guān)。本文對pvin的整體效果也進(jìn)行了研究,定義λ=pvin來反映磨削程度,管壁內(nèi)側(cè)記為λin,外側(cè)記為λout。以45°彎管為研究對象,結(jié)果如圖6所示,與p-L圖不同的是,λ-L圖在彎管肘部的前、后部也出現(xiàn)了差值,見圖6中A、C區(qū)域。這是由于管壁外側(cè)的流速先減小后增大,管壁內(nèi)側(cè)的流速先增大后減小造成的。但管子內(nèi)外側(cè)的磨削差異主要仍發(fā)生在肘部,即B區(qū)域,與圖5對比發(fā)現(xiàn),這主要是由于p的局部差異導(dǎo)致的,因此主要對p展開討論。此外,根據(jù)λ的變化趨勢可知,對彎管的進(jìn)出口進(jìn)行加工可以使A、C區(qū)域的差異抵消,從而僅在彎管肘部的B區(qū)域有較大的磨削差異。

    圖6 彎管局部λ的變化趨勢

    圖7所示為磨料濃度η為10%時的局部壓力差隨入口流速的變化曲線。由圖7可看出,入口流速是影響局部壓差的主要因素,這是由于水的初始流速越大,動能越大,從而壓力越大。從管子的彎曲角度看,在管子彎曲角度小于30°時,增大管子的彎曲角度會導(dǎo)致局部壓差有較為明顯的增大。管子的彎曲角度在30°~90°之間時,局部壓差受局部彎曲角度的影響較小。在管子彎曲處,由于流道形狀改變造成的局部阻力使主相水和第二相SiC顆粒在內(nèi)壁的內(nèi)、外側(cè)發(fā)生不同形式的碰撞,這種碰撞的不均勻性是導(dǎo)致內(nèi)外側(cè)壓差的根本原因。

    圖7 磨料濃度η=10%的局部壓差

    磨料濃度η為10%時的管子進(jìn)出口壓差見圖8。局部壓差最高可達(dá)進(jìn)出口壓差的23.94%,且局部壓差的最大絕對值為0.56 MPa,在實際的磨粒流加工過程中不會產(chǎn)生較大的不均勻性。

    圖8 磨料濃度η=10%的進(jìn)出口壓差

    3.2 磨料濃度對局部壓差的影響

    湍流動能反映了流動發(fā)展成為湍流或維持湍流狀態(tài)的能力,其表達(dá)式為k=1.5(uI)2,湍流強(qiáng)度I=0.16Re-1/8,u是平均流速。層流流動中,只有靠近壁面的流層起到了磨削作用,而湍流流動的相鄰流層之間會出現(xiàn)滑動、混合以及渦團(tuán),流體微團(tuán)的流動表現(xiàn)出不確定性。湍流的產(chǎn)生使得遠(yuǎn)離壁面的流層也有機(jī)會與壁面接觸,因此會達(dá)到比層流更好的磨削效果。

    表4~表6分別為磨料濃度η為10%、20%、30%時的局部壓差及進(jìn)出口壓差。壓差隨著磨料濃度的增大而增大,這是因為固體顆粒的增多會導(dǎo)致兩相流體更具黏滯性[14],對磨粒運(yùn)動的約束能力提高,使顆粒與壁面碰撞的次數(shù)增多,導(dǎo)致壁面受到的壓力增大,局部壓差和進(jìn)出口壓差也同樣增大。總體來看,在磨料濃度明顯增大時,局部壓差的增長并不顯著。從湍流動能的角度分析,湍流動能描述了顆粒與壁面碰撞時具有的能量,其值越大,說明磨粒對壁面的磨削作用越強(qiáng)。

    表4 磨料濃度η=10%的壓差

    表5 磨料濃度η=20%的壓差

    表6 磨料濃度η=30%的壓差

    圖9為湍流動能的分布云圖。這里以入口流速30 m/s的45°彎管為例進(jìn)行分析。磨料濃度由10%提高至20%時,壁面處的湍流動能有一定提高,而濃度由20%提高至30%時,湍流動能的提高并不明顯。這是由于隨著磨料濃度的增大,

    (a) η=10%

    流體的黏滯性增加,流動性則會降低,雖然顆粒與壁面的碰撞次數(shù)增多了,但是每次碰撞所傳遞的動能在降低,從而可以預(yù)測,高濃度下的磨粒流加工效果會變?nèi)酢?/p>

    3.3 加工過程中的產(chǎn)熱分析

    本文選取的模擬壁面均為非光滑壁面,固體顆粒與壁面的碰撞會產(chǎn)生一定熱量。圖10 為不同磨料濃度下管子肘部的溫度分布云圖。仍以45°彎管、入口流速30 m/s為例進(jìn)行分析。由圖10可看出,磨削過程中有熱量產(chǎn)生,流體流經(jīng)彎曲處時局部的磨料濃度分布不均勻,速度也發(fā)生改變。摩擦產(chǎn)熱量由磨料濃度、壓力、流速共同決定。磨料濃度越高,與壁面接觸摩擦的幾率越大、時間越長,從而產(chǎn)生更多熱量;壓力會影響固體顆粒對壁面的擠壓力大小,擠壓力越大,摩擦力越大,從而壓力升高,產(chǎn)熱量增多。流速大小會影響顆粒與壁面摩擦的劇烈程度,流速越大,產(chǎn)生的熱量越多,但是積聚在管壁的熱量并不多,這是由于導(dǎo)熱介質(zhì)水的比熱容較大,流體流速較快時,顆粒與壁面、流體與壁面摩擦產(chǎn)生的熱量迅速隨水的流動被帶走,因而熱量積聚使管壁溫度上升不足1℃。在實際的磨粒流加工中,可以忽略磨料濃度變化導(dǎo)致的管壁溫度的變化,以及此溫度變化對加工過程的影響。

    (a) η=10%

    圖11為45°彎管不同流速及磨料濃度下的管壁溫度升高量云圖。由圖可見,流速是影響熱量積聚的主要因素,因為流速不僅影響了產(chǎn)熱量的多少,還影響了熱量的散失趨勢;而磨料濃度僅僅影響了熱量的產(chǎn)生過程。此外,產(chǎn)熱的原因不僅是顆粒與壁面的摩擦,還存在一部分由流體的黏性產(chǎn)生的熱量,因此提高磨料濃度,管壁溫度的升高并不顯著。

    圖11 45°彎管溫度升高量云圖

    4 結(jié)論

    (1)建立了磨粒流打磨不同彎曲角度細(xì)長圓管的固液兩相流模型,開展不同條件下的磨粒流數(shù)值模擬研究,得到入口流速、磨料濃度、彎曲角度等參數(shù)對管內(nèi)壓力分布、顆粒跡線以及溫度分布的影響規(guī)律。

    (2)入口流速是影響流體在壁面處壓力變化的主要原因。入口流速小于30 m/s時,形成的管內(nèi)局部壓差較小,在實際的磨粒流加工過程中應(yīng)選擇較大的流速以提高加工效率。

    (3)磨料濃度也會影響流體在壁面處的壓力。磨料濃度在10%~30%范圍內(nèi)增大時,局部壓差和湍流動能的增長均不顯著。但是增大磨料濃度會提高加工成本,同時也會增加兩相流體的循環(huán)裝置堵塞、磨損加重的風(fēng)險,可以使用低濃度的磨粒流進(jìn)行實際加工。

    (4)在300 K的加工環(huán)境下,以水為介質(zhì)、SiC為磨粒進(jìn)行磨粒流加工,加工過程中積聚的熱量和溫度變化均較小,由溫度變化導(dǎo)致的磨料性質(zhì)的變化較小。在實際加工過程中,可以忽略溫度變化對加工過程的不利影響。

    猜你喜歡
    磨料磨粒湍流
    基于凸多面體碰撞檢測的虛擬砂輪建模研究
    超精密表面研拋磨粒的研究進(jìn)展
    硅酸鹽通報(2021年3期)2021-04-18 11:01:42
    單個鐵氧體磨粒尺寸檢測電磁仿真
    重氣瞬時泄漏擴(kuò)散的湍流模型驗證
    微晶剛玉磨粒磨削20CrMnTi鋼的數(shù)值模擬研究
    40Cr熱噴涂件抗植物磨料磨損性能試驗
    磨料水射流技術(shù)及其在水下結(jié)構(gòu)物切割中的應(yīng)用
    焊接(2015年5期)2015-07-18 11:03:40
    “青春期”湍流中的智慧引渡(三)
    “青春期”湍流中的智慧引渡(二)
    磨料水射流作用下混凝土損傷場的數(shù)值模擬
    九草在线视频观看| 国产免费视频播放在线视频| 亚洲av免费高清在线观看| 亚洲精品中文字幕在线视频| 美女中出高潮动态图| 尾随美女入室| 最新的欧美精品一区二区| xxx大片免费视频| 看免费av毛片| 久久久精品区二区三区| 国产精品久久久av美女十八| 免费女性裸体啪啪无遮挡网站| 最新中文字幕久久久久| 91精品国产国语对白视频| 亚洲欧美一区二区三区国产| 成人午夜精彩视频在线观看| 午夜福利在线观看免费完整高清在| 久久99蜜桃精品久久| 国产精品免费大片| 欧美日韩av久久| 午夜福利视频在线观看免费| 全区人妻精品视频| 各种免费的搞黄视频| 精品少妇内射三级| 午夜精品国产一区二区电影| 久久久久精品人妻al黑| 欧美日韩成人在线一区二区| 午夜久久久在线观看| 嫩草影院入口| 极品少妇高潮喷水抽搐| 国产熟女欧美一区二区| 最近最新中文字幕免费大全7| 国产免费视频播放在线视频| 成年av动漫网址| 亚洲成av片中文字幕在线观看 | 欧美人与善性xxx| 国产精品 国内视频| 18禁裸乳无遮挡动漫免费视频| 国产精品久久久久久av不卡| 久久人人97超碰香蕉20202| 高清视频免费观看一区二区| 男人操女人黄网站| 日韩成人av中文字幕在线观看| av线在线观看网站| 黑丝袜美女国产一区| 日韩中字成人| 少妇猛男粗大的猛烈进出视频| 妹子高潮喷水视频| 亚洲精品美女久久久久99蜜臀 | 国产精品成人在线| 国产精品欧美亚洲77777| 亚洲四区av| 国产日韩欧美亚洲二区| 中文字幕人妻丝袜制服| 国产男人的电影天堂91| 国产精品久久久久久精品电影小说| tube8黄色片| 黑丝袜美女国产一区| 久久这里只有精品19| 色视频在线一区二区三区| 国产日韩欧美视频二区| 午夜免费男女啪啪视频观看| 亚洲一区二区三区欧美精品| 亚洲情色 制服丝袜| 97精品久久久久久久久久精品| 男人添女人高潮全过程视频| 在线天堂最新版资源| 亚洲国产精品专区欧美| 精品久久久精品久久久| 精品国产一区二区三区四区第35| 国产女主播在线喷水免费视频网站| 日韩三级伦理在线观看| 中文字幕精品免费在线观看视频 | 亚洲欧洲日产国产| 人妻少妇偷人精品九色| 久久久欧美国产精品| 日韩一本色道免费dvd| 国产亚洲午夜精品一区二区久久| 一边摸一边做爽爽视频免费| 80岁老熟妇乱子伦牲交| 亚洲激情五月婷婷啪啪| 国产成人av激情在线播放| 天堂俺去俺来也www色官网| 18+在线观看网站| 亚洲av中文av极速乱| 一级毛片电影观看| 免费看av在线观看网站| 90打野战视频偷拍视频| 黑人高潮一二区| 我要看黄色一级片免费的| 丝袜喷水一区| 全区人妻精品视频| 少妇的丰满在线观看| 高清欧美精品videossex| 精品亚洲成国产av| 边亲边吃奶的免费视频| 青春草亚洲视频在线观看| 免费黄色在线免费观看| 精品少妇黑人巨大在线播放| 草草在线视频免费看| 91精品伊人久久大香线蕉| 久热久热在线精品观看| 伦理电影大哥的女人| 久久久久精品久久久久真实原创| 99久久精品国产国产毛片| av在线老鸭窝| 男女边摸边吃奶| 欧美精品av麻豆av| 国产69精品久久久久777片| 人人妻人人爽人人添夜夜欢视频| 精品国产乱码久久久久久小说| 黄色毛片三级朝国网站| 精品亚洲乱码少妇综合久久| 人成视频在线观看免费观看| 精品第一国产精品| 亚洲国产成人一精品久久久| 国产有黄有色有爽视频| 激情五月婷婷亚洲| 女人被躁到高潮嗷嗷叫费观| 肉色欧美久久久久久久蜜桃| 久久久精品94久久精品| 婷婷色av中文字幕| 成人亚洲精品一区在线观看| av片东京热男人的天堂| 中文字幕精品免费在线观看视频 | 青青草视频在线视频观看| 亚洲精品456在线播放app| 看免费成人av毛片| 亚洲精品乱久久久久久| 青春草国产在线视频| 欧美人与善性xxx| 国产一区二区在线观看日韩| www日本在线高清视频| 亚洲国产看品久久| freevideosex欧美| 91精品伊人久久大香线蕉| 男女午夜视频在线观看 | 少妇熟女欧美另类| 人妻一区二区av| 欧美人与性动交α欧美精品济南到 | 国产成人精品无人区| 国产精品欧美亚洲77777| 建设人人有责人人尽责人人享有的| 欧美人与善性xxx| 久久精品久久久久久噜噜老黄| 9191精品国产免费久久| 久久狼人影院| 观看av在线不卡| 91精品伊人久久大香线蕉| 久久久国产欧美日韩av| 亚洲一区二区三区欧美精品| 一区二区三区四区激情视频| 国产精品久久久av美女十八| 两性夫妻黄色片 | 成年女人在线观看亚洲视频| 国产熟女午夜一区二区三区| 热re99久久国产66热| 国产日韩欧美亚洲二区| 成人毛片a级毛片在线播放| av一本久久久久| 亚洲av综合色区一区| 五月开心婷婷网| 各种免费的搞黄视频| 欧美精品亚洲一区二区| 国产探花极品一区二区| 亚洲精品乱久久久久久| 另类精品久久| 黑丝袜美女国产一区| 一级片'在线观看视频| 少妇人妻精品综合一区二区| 国产一区二区激情短视频 | 国产免费视频播放在线视频| 成年女人在线观看亚洲视频| 最黄视频免费看| 嫩草影院入口| 国产淫语在线视频| 五月伊人婷婷丁香| 亚洲四区av| kizo精华| 菩萨蛮人人尽说江南好唐韦庄| 99久久精品国产国产毛片| 国产av一区二区精品久久| av播播在线观看一区| 国产激情久久老熟女| 日韩一区二区三区影片| 国产男女超爽视频在线观看| 91久久精品国产一区二区三区| 亚洲国产欧美日韩在线播放| 插逼视频在线观看| 在线观看国产h片| 熟妇人妻不卡中文字幕| 天天躁夜夜躁狠狠久久av| 人人澡人人妻人| 极品少妇高潮喷水抽搐| 一区二区三区精品91| 美女福利国产在线| 久久这里有精品视频免费| 国产 一区精品| 性色avwww在线观看| 日韩欧美一区视频在线观看| 亚洲熟女精品中文字幕| 春色校园在线视频观看| 另类亚洲欧美激情| 最新中文字幕久久久久| 下体分泌物呈黄色| 深夜精品福利| 午夜激情久久久久久久| 纯流量卡能插随身wifi吗| 中文字幕人妻丝袜制服| 久久青草综合色| 国产精品久久久av美女十八| 少妇被粗大猛烈的视频| a级毛片黄视频| 精品人妻熟女毛片av久久网站| av片东京热男人的天堂| 蜜桃国产av成人99| 国产成人精品久久久久久| 亚洲情色 制服丝袜| av在线app专区| 亚洲激情五月婷婷啪啪| 国语对白做爰xxxⅹ性视频网站| 人妻 亚洲 视频| 五月伊人婷婷丁香| 亚洲精品第二区| www.av在线官网国产| 日韩中文字幕视频在线看片| 午夜精品国产一区二区电影| 亚洲精品久久久久久婷婷小说| 亚洲高清免费不卡视频| 熟女av电影| 亚洲国产av影院在线观看| 丝袜脚勾引网站| 黄色配什么色好看| 狂野欧美激情性xxxx在线观看| 高清毛片免费看| 成人毛片a级毛片在线播放| 国产av国产精品国产| 夜夜爽夜夜爽视频| 欧美人与性动交α欧美精品济南到 | 狂野欧美激情性xxxx在线观看| 老熟女久久久| 中文欧美无线码| av电影中文网址| 午夜福利,免费看| 18禁裸乳无遮挡动漫免费视频| 丰满饥渴人妻一区二区三| 九九在线视频观看精品| 精品少妇久久久久久888优播| 一个人免费看片子| av在线观看视频网站免费| 最后的刺客免费高清国语| 亚洲精华国产精华液的使用体验| 久久久久视频综合| 亚洲美女搞黄在线观看| 成人手机av| 又黄又粗又硬又大视频| 捣出白浆h1v1| 中国国产av一级| 日日啪夜夜爽| 69精品国产乱码久久久| 久久人人爽人人片av| 大陆偷拍与自拍| 中国国产av一级| tube8黄色片| 亚洲天堂av无毛| a级毛片黄视频| videosex国产| 人人妻人人添人人爽欧美一区卜| 午夜福利乱码中文字幕| 一本—道久久a久久精品蜜桃钙片| 激情视频va一区二区三区| 免费女性裸体啪啪无遮挡网站| 人人妻人人添人人爽欧美一区卜| 亚洲精品久久成人aⅴ小说| 亚洲美女视频黄频| 亚洲成人一二三区av| 亚洲一码二码三码区别大吗| 一级毛片我不卡| 久久久久久久久久成人| 国产亚洲最大av| 在线观看人妻少妇| 又大又黄又爽视频免费| 国产无遮挡羞羞视频在线观看| 午夜福利网站1000一区二区三区| 男男h啪啪无遮挡| av视频免费观看在线观看| 一级毛片 在线播放| 亚洲av综合色区一区| 高清黄色对白视频在线免费看| 91精品国产国语对白视频| 午夜免费男女啪啪视频观看| 国产男女内射视频| 97在线人人人人妻| 久久人妻熟女aⅴ| 啦啦啦在线观看免费高清www| 一边摸一边做爽爽视频免费| 看十八女毛片水多多多| 精品久久蜜臀av无| 日韩一区二区三区影片| 五月开心婷婷网| 亚洲精品国产色婷婷电影| tube8黄色片| 美女主播在线视频| 啦啦啦啦在线视频资源| 男女啪啪激烈高潮av片| 在线 av 中文字幕| 免费女性裸体啪啪无遮挡网站| 久久 成人 亚洲| 日韩在线高清观看一区二区三区| 亚洲av福利一区| 国产在线一区二区三区精| 如日韩欧美国产精品一区二区三区| 国产片特级美女逼逼视频| 亚洲av电影在线进入| 国产精品女同一区二区软件| 久久久精品免费免费高清| av免费在线看不卡| 亚洲精品一区蜜桃| 精品国产一区二区久久| 内地一区二区视频在线| 街头女战士在线观看网站| 男人爽女人下面视频在线观看| videosex国产| 久久国产亚洲av麻豆专区| 亚洲精华国产精华液的使用体验| www.熟女人妻精品国产 | 亚洲欧洲精品一区二区精品久久久 | 国产男人的电影天堂91| 久久精品国产a三级三级三级| 91aial.com中文字幕在线观看| 欧美日韩av久久| 午夜福利在线观看免费完整高清在| av卡一久久| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 五月玫瑰六月丁香| av在线老鸭窝| 亚洲精品乱久久久久久| av在线老鸭窝| 黑人高潮一二区| 极品少妇高潮喷水抽搐| av片东京热男人的天堂| av黄色大香蕉| www.色视频.com| 狠狠精品人妻久久久久久综合| 国产极品天堂在线| 纯流量卡能插随身wifi吗| 亚洲激情五月婷婷啪啪| 韩国av在线不卡| 国产女主播在线喷水免费视频网站| 午夜福利影视在线免费观看| 国产熟女欧美一区二区| 一本久久精品| 插逼视频在线观看| 夜夜骑夜夜射夜夜干| 国产男女超爽视频在线观看| av国产久精品久网站免费入址| 国产欧美亚洲国产| 色视频在线一区二区三区| 国产不卡av网站在线观看| 九色成人免费人妻av| 欧美97在线视频| 777米奇影视久久| av线在线观看网站| 亚洲综合色惰| 亚洲欧美色中文字幕在线| 日韩大片免费观看网站| av在线老鸭窝| 又粗又硬又长又爽又黄的视频| 天天影视国产精品| 国产男女内射视频| 午夜免费观看性视频| 国产精品免费大片| 国产又色又爽无遮挡免| 9色porny在线观看| 精品少妇黑人巨大在线播放| 黄色毛片三级朝国网站| 国产男女超爽视频在线观看| 久久国产亚洲av麻豆专区| 精品酒店卫生间| 亚洲欧美日韩卡通动漫| 色94色欧美一区二区| 男男h啪啪无遮挡| 亚洲成色77777| 美女内射精品一级片tv| 波野结衣二区三区在线| 9191精品国产免费久久| 日本-黄色视频高清免费观看| 一本久久精品| 永久免费av网站大全| av女优亚洲男人天堂| av网站免费在线观看视频| 久久久久精品性色| a级毛片在线看网站| 九色亚洲精品在线播放| 亚洲精品,欧美精品| 另类精品久久| 人人妻人人澡人人看| 黑人巨大精品欧美一区二区蜜桃 | 成年女人在线观看亚洲视频| 久久ye,这里只有精品| 尾随美女入室| 日日摸夜夜添夜夜爱| 色94色欧美一区二区| 亚洲国产最新在线播放| 亚洲第一区二区三区不卡| 中文字幕精品免费在线观看视频 | 母亲3免费完整高清在线观看 | av网站免费在线观看视频| 免费黄网站久久成人精品| 日韩精品免费视频一区二区三区 | freevideosex欧美| 国产极品粉嫩免费观看在线| 五月玫瑰六月丁香| 18禁在线无遮挡免费观看视频| 国产精品久久久久久精品古装| 亚洲精品美女久久久久99蜜臀 | 国产精品.久久久| 哪个播放器可以免费观看大片| 一二三四中文在线观看免费高清| 最近2019中文字幕mv第一页| h视频一区二区三区| 亚洲精品,欧美精品| 激情五月婷婷亚洲| 欧美3d第一页| 精品人妻熟女毛片av久久网站| 欧美亚洲 丝袜 人妻 在线| 日韩成人av中文字幕在线观看| 18禁动态无遮挡网站| 亚洲久久久国产精品| 七月丁香在线播放| 亚洲高清免费不卡视频| 亚洲图色成人| 亚洲一码二码三码区别大吗| 18禁裸乳无遮挡动漫免费视频| 国产精品嫩草影院av在线观看| 午夜视频国产福利| 日韩一本色道免费dvd| 国精品久久久久久国模美| 18禁观看日本| 成人亚洲精品一区在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品久久久久久精品电影小说| 欧美+日韩+精品| 一区在线观看完整版| 久久精品久久精品一区二区三区| 久久精品国产a三级三级三级| 性高湖久久久久久久久免费观看| 最近最新中文字幕免费大全7| 国产亚洲精品久久久com| 少妇精品久久久久久久| 最新的欧美精品一区二区| 制服人妻中文乱码| 久久久久久人妻| 亚洲一码二码三码区别大吗| 国产亚洲午夜精品一区二区久久| 人妻一区二区av| 少妇高潮的动态图| 日日撸夜夜添| 韩国av在线不卡| 99re6热这里在线精品视频| 精品久久久久久电影网| 日本欧美国产在线视频| 亚洲欧美色中文字幕在线| 亚洲精品乱码久久久久久按摩| 欧美性感艳星| 精品亚洲成国产av| 久久99蜜桃精品久久| 国产免费一区二区三区四区乱码| 日韩视频在线欧美| freevideosex欧美| 亚洲欧美清纯卡通| 精品人妻在线不人妻| 久久久久久伊人网av| 老司机亚洲免费影院| 亚洲精华国产精华液的使用体验| 国产一区二区在线观看日韩| 大片免费播放器 马上看| 国产精品不卡视频一区二区| av免费观看日本| 一本久久精品| 久久久欧美国产精品| 最黄视频免费看| 欧美日韩av久久| 国产成人一区二区在线| 中国三级夫妇交换| 青春草国产在线视频| 国产免费现黄频在线看| 免费av中文字幕在线| 欧美 日韩 精品 国产| 日韩人妻精品一区2区三区| 综合色丁香网| av黄色大香蕉| 婷婷色麻豆天堂久久| 国产av国产精品国产| 少妇的逼水好多| 成人黄色视频免费在线看| 国产精品熟女久久久久浪| 99国产精品免费福利视频| 亚洲综合色惰| 欧美日韩亚洲高清精品| 男女啪啪激烈高潮av片| 女性被躁到高潮视频| videosex国产| 九色成人免费人妻av| 97在线人人人人妻| 国产永久视频网站| 国产精品久久久久久av不卡| 美女中出高潮动态图| 国产精品秋霞免费鲁丝片| 亚洲精品成人av观看孕妇| 亚洲中文av在线| 国产亚洲最大av| 九九爱精品视频在线观看| av在线观看视频网站免费| 丝袜美足系列| 午夜福利,免费看| 少妇的丰满在线观看| 婷婷成人精品国产| 母亲3免费完整高清在线观看 | 亚洲欧美一区二区三区黑人 | 中文欧美无线码| 精品一区二区三区视频在线| 久久久久精品人妻al黑| 久久人人爽人人爽人人片va| 在线亚洲精品国产二区图片欧美| 寂寞人妻少妇视频99o| av卡一久久| 久久这里有精品视频免费| www.色视频.com| 97精品久久久久久久久久精品| 久久av网站| 男女边吃奶边做爰视频| 亚洲国产欧美日韩在线播放| 国产精品国产三级专区第一集| 丰满少妇做爰视频| 亚洲国产日韩一区二区| 有码 亚洲区| av黄色大香蕉| 国产欧美日韩综合在线一区二区| 国产精品99久久99久久久不卡 | 多毛熟女@视频| 国产日韩欧美在线精品| 亚洲熟女精品中文字幕| 狂野欧美激情性bbbbbb| 午夜福利在线观看免费完整高清在| 国产成人av激情在线播放| 热99久久久久精品小说推荐| 十分钟在线观看高清视频www| av线在线观看网站| 欧美日韩成人在线一区二区| 少妇 在线观看| 国产熟女欧美一区二区| 午夜福利视频精品| 亚洲国产精品999| 久久久国产一区二区| av在线老鸭窝| 久久久久精品久久久久真实原创| 成人二区视频| 日韩av不卡免费在线播放| 成年人午夜在线观看视频| 精品酒店卫生间| 久久97久久精品| 97超碰精品成人国产| 97在线视频观看| 久久久a久久爽久久v久久| 极品人妻少妇av视频| 婷婷色av中文字幕| 亚洲美女黄色视频免费看| 黄色视频在线播放观看不卡| 亚洲av欧美aⅴ国产| 亚洲国产最新在线播放| 99精国产麻豆久久婷婷| 91在线精品国自产拍蜜月| 不卡视频在线观看欧美| 欧美另类一区| 国产午夜精品一二区理论片| 少妇的丰满在线观看| 国产一区有黄有色的免费视频| 26uuu在线亚洲综合色| 狠狠精品人妻久久久久久综合| 汤姆久久久久久久影院中文字幕| 久久精品夜色国产| 欧美日韩一区二区视频在线观看视频在线| 大片免费播放器 马上看| 久久久久国产网址| 国产综合精华液| 国产伦理片在线播放av一区| 激情五月婷婷亚洲| 久久毛片免费看一区二区三区| 国产亚洲欧美精品永久| 精品一区二区三卡| 久久久久久伊人网av| 97在线人人人人妻| 日本wwww免费看| 高清在线视频一区二区三区| av福利片在线| 熟妇人妻不卡中文字幕| 制服诱惑二区| 欧美日韩亚洲高清精品| 国产精品偷伦视频观看了| 久久99热6这里只有精品| 男女边摸边吃奶| 秋霞在线观看毛片| 人成视频在线观看免费观看| 九色亚洲精品在线播放| 亚洲国产精品一区二区三区在线| www日本在线高清视频| 你懂的网址亚洲精品在线观看| 日本欧美国产在线视频| 九九在线视频观看精品| 国语对白做爰xxxⅹ性视频网站| 亚洲成色77777| 中文字幕精品免费在线观看视频 | 久久人人爽av亚洲精品天堂| 国产精品人妻久久久久久|