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

    空氣靜壓軸承氣旋機(jī)理及其影響因素分析

    2021-10-20 10:58:14唐璐陽張小龍
    機(jī)械設(shè)計(jì)與制造 2021年10期
    關(guān)鍵詞:大渦氣膜氣旋

    唐璐陽,張小龍,寇 盼

    (西安建筑科技大學(xué)機(jī)電工程學(xué)院,陜西 西安710055)

    1 引言

    空氣靜壓電主軸以空氣靜壓軸承支承、以電機(jī)或者渦輪為動(dòng)力驅(qū)動(dòng),帶動(dòng)刀具高速旋轉(zhuǎn),從而實(shí)現(xiàn)高速精密加工,具有回轉(zhuǎn)精度高、低振動(dòng)、轉(zhuǎn)速高等特點(diǎn),在機(jī)床、精密儀器、醫(yī)療器械等領(lǐng)域均得到廣泛應(yīng)用[1]。

    作為超精密加工的重要部件,空氣靜壓軸承的性能影響著被加工精密零件的質(zhì)量[2]。近幾年來,關(guān)于空氣靜壓軸承的微振動(dòng)現(xiàn)象已有很多研究成果。文獻(xiàn)[3]通過實(shí)驗(yàn)發(fā)現(xiàn)了氣體通過供氣孔進(jìn)入軸承氣膜時(shí)具有明顯的壓力降現(xiàn)象。文獻(xiàn)[4]通過對(duì)超精密氣浮工作臺(tái)的振動(dòng)特性分析,發(fā)現(xiàn)了高速氣體通過有腔節(jié)流孔時(shí)產(chǎn)生“氣旋”現(xiàn)象,會(huì)導(dǎo)致氣膜波動(dòng)。文獻(xiàn)[5]利用Star-CD軟件仿真了S型和R型節(jié)流孔的內(nèi)部氣體流動(dòng)情況,并通過實(shí)驗(yàn)驗(yàn)證了不同節(jié)流孔形狀對(duì)軸承微振動(dòng)有影響。文獻(xiàn)[6]進(jìn)一步研究了S型和R型節(jié)流孔對(duì)軸承承載力和穩(wěn)定性的影響。為了抑制微振動(dòng)提高加工精度,文獻(xiàn)[7]實(shí)驗(yàn)研究發(fā)現(xiàn)了增加均壓槽可降低軸承微振動(dòng)80%左右,能顯著提高軸承的穩(wěn)定性。文獻(xiàn)[8]研究認(rèn)為由于空氣軸承節(jié)流孔內(nèi)氣流雷諾數(shù)過大而導(dǎo)致了微振動(dòng),并提出了減小微振動(dòng)的方法。文獻(xiàn)[9]和文獻(xiàn)O[10]考慮是邊界區(qū)和慣性區(qū)、激波與邊界層的相互作用,在供氣壓力較高和氣膜間隙較大的條件下,發(fā)現(xiàn)在節(jié)流孔入口和出口轉(zhuǎn)角處存在激波現(xiàn)象,從而導(dǎo)致氣膜溫度和壓力下降,甚至使節(jié)流孔內(nèi)的氣體流動(dòng)由層流向湍流轉(zhuǎn)化,使軸承出現(xiàn)微振動(dòng)。

    針對(duì)上述研究,目前研究認(rèn)為節(jié)流孔壓力腔內(nèi)氣體流動(dòng)的非均勻性引起的“氣旋”現(xiàn)象是軸承微振動(dòng)的主要原因[11],研究相對(duì)集中在求解氣旋穩(wěn)態(tài)特性下不同工藝參數(shù)和幾何參數(shù)對(duì)氣旋的影響,但以上研究對(duì)氣旋瞬態(tài)發(fā)展過程和微振動(dòng)的產(chǎn)生機(jī)制研究相對(duì)較少。利用STAR-ccm+軟件對(duì)節(jié)流孔內(nèi)氣旋的產(chǎn)生及其瞬態(tài)變化情況進(jìn)行模擬,確定了氣旋導(dǎo)致微振動(dòng)的產(chǎn)生機(jī)制、并從穩(wěn)態(tài)和瞬態(tài)兩方面研究了不同軸承供氣壓力和壓力腔形狀對(duì)氣旋強(qiáng)度的影響,為空氣靜壓軸承的設(shè)計(jì)提供理論依據(jù)。

    2 氣旋及其對(duì)微振動(dòng)的作用

    2.1 大渦模擬

    為了模擬節(jié)流孔內(nèi)流體的動(dòng)態(tài)變化情況,這里利用這里進(jìn)行大渦模擬時(shí),我們利用STAR-ccm+中的大渦模擬進(jìn)行節(jié)流孔內(nèi)流場的瞬態(tài)求解,STAR-ccm+可以更加清楚獲得流場內(nèi)的分布情況,使用方便高效,以其穩(wěn)健精確的數(shù)值算法、覆蓋面廣的理論模型以及易處理的網(wǎng)格體系,是一款優(yōu)秀的CFD模擬軟件,是新一代CFD軟件的最強(qiáng)閃光點(diǎn)[12]。

    大渦模擬[13]將湍流瞬時(shí)脈動(dòng)分解為直接模擬大尺度脈動(dòng)和輸送方程模擬小尺度脈動(dòng),大尺度脈動(dòng)通過N~S方程直接求解,而小尺度脈動(dòng)通過亞格力應(yīng)力模型來求解,計(jì)算量相比直接模擬較小,模型比較容易構(gòu)造,能夠捕捉到雷諾平均法無法得到的許多非穩(wěn)態(tài)和非平穩(wěn)過程中的湍流過程。通過大渦模擬能精確求解節(jié)流孔內(nèi)流體的瞬時(shí)流動(dòng)情況,便于分析氣旋以及微振動(dòng)的機(jī)理。

    利用大渦模擬時(shí)采用的節(jié)流孔形狀及計(jì)算域如圖1所示。節(jié)流孔是空氣靜壓設(shè)計(jì)的關(guān)鍵,節(jié)流孔的直徑一般在(0.1~1.0)mm[14],且計(jì)算域主要幾何尺寸分別為節(jié)流孔直徑d=0.24 mm,氣膜厚度h0=0.013mm,節(jié)流孔深度l=1mm,壓力腔寬度B=2mm,壓力腔深度H=0.1mm。且氣膜間隙區(qū)域的上邊界與軸承內(nèi)表面相接觸,氣膜區(qū)域的下邊界與主軸表面相接觸。另外在分析時(shí)節(jié)省時(shí)間,只取計(jì)算域模型的1/4進(jìn)行分析。

    圖1 計(jì)算域模型圖Fig.1 Computational Domain Model Diagram

    利用STAR-ccm+進(jìn)行大渦模擬時(shí),將雷諾平均法的穩(wěn)態(tài)計(jì)算值作為大渦模擬的計(jì)算初始值,這樣保證為大渦模擬提供準(zhǔn)確的初始流場以減少計(jì)算時(shí)間。這里我們將雷諾平均法的穩(wěn)態(tài)計(jì)算值作為大渦模擬的計(jì)算初始值,可以為大渦模擬提供準(zhǔn)確的初始流場以減少計(jì)算時(shí)間。流場計(jì)算域采用六面體結(jié)構(gòu)進(jìn)行網(wǎng)格劃分,進(jìn)氣口供氣壓力為0.435MPa,出氣口供氣壓力為0.1MPa。Model選擇PISO非穩(wěn)態(tài)計(jì)算,Turbulence選擇大漩渦模擬,選擇氣體動(dòng)力粘度屬性為Sutherland定律,本次大渦模擬選用WALE模型來模擬節(jié)流孔湍流數(shù)值計(jì)算,設(shè)置WALE模型常數(shù)Cw為0.325。由于大渦模擬為非定常計(jì)算,在時(shí)間離散上采用二階隱式模型,設(shè)置迭代時(shí)間步長Δt為10-7s,時(shí)間步長的選擇要滿足CFL條件,即:

    式中:μ-流經(jīng)網(wǎng)格單元的流體速度大小,m/s;Δx-速度方向上的最小網(wǎng)格尺寸,m;Δt-迭代時(shí)間步長。單個(gè)時(shí)間步長內(nèi)迭代20次,最后根據(jù)收斂情況適當(dāng)調(diào)整亞松弛因子保證殘差收斂。

    由于實(shí)際節(jié)流孔內(nèi)氣體流動(dòng)是非常復(fù)雜的,因此我們需要做以下簡單假設(shè):(1)運(yùn)動(dòng)環(huán)境溫度為恒溫過程;(2)設(shè)氣膜厚度h0為常值;(3)在氣膜區(qū)域內(nèi)垂直于軸承間隙的壓力和密度相等。

    2.2 氣旋及微振動(dòng)

    節(jié)流孔段類似于拉伐爾噴管[15],從節(jié)流孔進(jìn)氣口區(qū)域供入的氣體噴射到較大的自由空間,這部分氣體流動(dòng)恰好為流體的沖擊射流,沖擊射流的示意圖,如圖2所示,氣流流動(dòng)按照特征可以分為自由射流區(qū)、阻滯區(qū)以及壁面射流區(qū)。

    圖2 沖擊射流二維示意圖Fig.2 Two-dimensional Schematic Diagram of Impinging Jets

    相比自由射流區(qū)、阻滯區(qū)以及壁面射流區(qū),為了詳細(xì)分析節(jié)流孔內(nèi)流場分布情況,根據(jù)節(jié)流孔形狀和氣體流動(dòng)情況,把流場計(jì)算域分為四個(gè)不同區(qū)域,如圖3所示。

    圖3 內(nèi)部流場過程示意圖Fig.3 Internal Fow Field Process Diagram

    處于進(jìn)氣口下的①自由射流區(qū),射流的邊界接觸著節(jié)流孔壁面,環(huán)境流體與射流邊界之間的相互剪切作用,從而使動(dòng)能、能量在發(fā)生空吸作用,空吸作用使得射流內(nèi)徑向速度的分布發(fā)生變化。氣流將沿①區(qū)域保持高速向下流動(dòng)。

    當(dāng)氣流到達(dá)②阻滯區(qū)域,氣流流線發(fā)生彎曲,氣流產(chǎn)生回流和分離,從而導(dǎo)致氣流紊亂,形成了湍流效應(yīng),湍動(dòng)能強(qiáng)度較大。垂直氣流開始接觸到主軸表面,氣流軸向速度急劇減少且速度方向改變90°,導(dǎo)致氣體速度梯度的產(chǎn)生,對(duì)比圖3中的速度云圖可以清楚看到產(chǎn)生速度梯度的區(qū)域大約是②區(qū)域中的2d范圍內(nèi)。

    在空吸作用下,氣體與周圍的靜止氣體相接觸形成了氣旋,此時(shí)形成氣旋的強(qiáng)度較大,氣旋的中心處是壓力最低的地方,氣旋邊緣與中心產(chǎn)生的壓力差推動(dòng)氣旋運(yùn)動(dòng),氣旋運(yùn)動(dòng)的氣旋受到周圍氣體的粘性作用,使得氣旋強(qiáng)度逐漸減少,隨著從①自由射流區(qū)域氣體的不斷輸入,氣旋不斷的生成和破裂,從而引起劇烈的壓力波動(dòng),導(dǎo)致軸承的微振動(dòng)最大。圖4是利用大渦模擬得到壓力腔內(nèi)部的瞬時(shí)流態(tài)圖,可以看到氣旋生成到破裂消散的過程。

    圖4 不同時(shí)刻壓力腔的流場速度和壓力云圖Fig.4 Flow Field transient Velocity and Pressure Cloud Diagram of the Pressure Chamber

    在進(jìn)入壁面射流區(qū)之前,氣流首先要沿壓力腔徑向位置流動(dòng),氣流會(huì)與壓力腔壁面發(fā)生碰撞,繼續(xù)又在壓力腔內(nèi)產(chǎn)生氣旋,但相比在②阻滯區(qū)域,氣旋強(qiáng)度較小,產(chǎn)生的氣旋影響壓力腔內(nèi)壓力波動(dòng),持續(xù)對(duì)軸承微振動(dòng)產(chǎn)生影響。

    氣流經(jīng)過壁面射流區(qū)時(shí),氣流總體方向沿徑向向外流動(dòng),各個(gè)方向上的氣流在壓力腔出口處相匯合,共同從壓力腔出口流入氣膜區(qū)域,當(dāng)從壓力腔出口流向氣膜區(qū)域時(shí),氣流流動(dòng)區(qū)域面積由厘米變至微米級(jí)別,造成動(dòng)能增強(qiáng),各個(gè)方向上的速度在壓力腔出口處形成湍流,形成的湍動(dòng)能強(qiáng)度和壓力波動(dòng)相比②阻滯區(qū)域較小,也會(huì)對(duì)軸承微振動(dòng)產(chǎn)生影響。在氣膜區(qū)域流動(dòng)的氣體沿軸向流動(dòng),氣流主要受主軸-軸承壁面的摩擦力作用,在摩擦力的作用下氣流速度越來越小,流動(dòng)狀態(tài)為層流,對(duì)軸承的微振動(dòng)幾乎沒有影響。

    另外,分析不同時(shí)刻下計(jì)算域的瞬時(shí)渦量分布圖,氣體在①自由射流區(qū)域高速向下流動(dòng)時(shí),這部分氣體包含有極大的能量,由于不存在空吸作用,此時(shí)的渦量為零;當(dāng)氣流抵達(dá)②阻滯區(qū)域,在“2d”范圍內(nèi)最先形成氣旋,此處的渦量較大,對(duì)比圖5(a)-(c)過程可以看出節(jié)流孔出口形成的氣旋在徑向速度的驅(qū)動(dòng)下,沿徑向開始移動(dòng),但是受氣流粘性作用,空吸作用變?nèi)蹩梢钥吹綒庑龔?qiáng)度明顯減小,從渦量分布圖可以看出渦量有所減少。在壓力腔內(nèi)由于運(yùn)動(dòng)氣流和壓力腔壁面產(chǎn)生碰撞,又有氣旋的產(chǎn)生,壓力腔內(nèi)氣旋持續(xù)影響氣膜壓力波動(dòng)。從圖5(d)看到,渦量在壓力腔出口處集中,壓力和速度波動(dòng)對(duì)軸承微振動(dòng)產(chǎn)生影響。

    圖5 不同時(shí)刻下計(jì)算域的渦量分布云圖Fig.5 Transient Vorticity Distribution Cloud Map of the Computational Domain

    最后進(jìn)入氣膜區(qū)域的氣體沿軸向流動(dòng),氣流受摩擦力作用速度越來越小,不會(huì)產(chǎn)生氣旋,渦量逐漸減為零。通過以上對(duì)氣體在節(jié)流孔內(nèi)的流動(dòng)情況分析,導(dǎo)致軸承產(chǎn)生微振動(dòng)主要有以下幾方面:(1)從沖量角度來講;氣旋與壓力腔不平行且垂直的速度分量會(huì)導(dǎo)致主軸-軸承微振動(dòng);(2)從壓力變化來講;氣旋造成壓力降產(chǎn)生壓力脈動(dòng),這種壓力脈動(dòng)作用在主軸表面產(chǎn)生微振動(dòng);(3)在壁面射流區(qū)各個(gè)方向的氣流在壓力腔出口處匯聚,這種不穩(wěn)定的氣流形態(tài)造成湍動(dòng)能的增加而導(dǎo)致微振動(dòng)。另外空氣靜壓主軸在運(yùn)轉(zhuǎn)時(shí),主要的內(nèi)部熱源為空氣靜壓軸承的摩擦生熱[16]。因?yàn)橄職饽づc運(yùn)轉(zhuǎn)中的高速主軸相接觸,氣膜上下邊界之間形成較大的速度梯度,從而產(chǎn)生大量的剪切摩擦熱,這些剪切摩擦熱會(huì)導(dǎo)致主軸熱變形,進(jìn)而影響氣膜間隙,造成氣膜壓力波動(dòng)。

    3 影響氣旋的因素分析

    氣旋是軸承微振動(dòng)的主要原因,影響氣旋的因素主要由結(jié)構(gòu)幾何參數(shù)和工藝參數(shù)兩部分構(gòu)成,這里對(duì)不同供氣壓力和壓力腔形狀的穩(wěn)態(tài)和瞬態(tài)特性進(jìn)行分析。

    3.1 不同供氣壓力下的穩(wěn)態(tài)分析

    供氣壓力的大小嚴(yán)重影響軸承的工作性能,適當(dāng)增加供氣壓力可以增加軸承的承載力和剛度,但是當(dāng)供氣壓力過大時(shí),可能導(dǎo)致軸承工作不穩(wěn)定,增加軸承氣錘振動(dòng)發(fā)生的可能性。為防止出現(xiàn)“氣錘”現(xiàn)象,仿真了供氣壓力分別等于0.335MPa和0.535MPa兩種情況下節(jié)流孔的穩(wěn)態(tài)內(nèi)部流場。并用壓力降大小和湍動(dòng)能大小表示氣旋強(qiáng)度大小。

    通過繪制不同供氣壓力下節(jié)流孔對(duì)稱中心沿徑向位置的壓力和湍動(dòng)能變化曲線,從圖6壓力分大小分布曲線圖可以看出氣旋的中心為曲線的最低點(diǎn),當(dāng)供氣壓力為0.335MPa時(shí),氣旋導(dǎo)致的壓力降為9134Pa;當(dāng)供氣壓力為0.535MPa時(shí),氣旋導(dǎo)致的壓力降分別為20181Pa。因此隨著供氣壓力的增加,氣旋導(dǎo)致的壓力降增大。從圖7湍動(dòng)能大小分布曲線圖可以看出供氣壓力越大,湍動(dòng)能越大,引起軸承的微振動(dòng)越大。

    圖6 不同供氣壓力下壓力沿徑向位置分布曲線圖Fig.6 Distribution Curve of Pressure Along Radial Position Under Different Air Supply Pressure

    圖7 不同供氣壓力下湍動(dòng)能沿徑向位置分布曲線圖Fig.7 Distribution Curve of Turbulent Kinetic Energy Along Radial Position under Different Air Supply Pressure

    3.2 不同供氣壓力下的瞬態(tài)分析

    如圖8所示,為不同供氣壓力下的計(jì)算域瞬態(tài)速度和壓力云圖,不同時(shí)刻下,另外觀察瞬態(tài)下的計(jì)算域氣流流態(tài)情況,可以看到氣旋形成到消散的過程,在節(jié)流孔出口處的氣旋強(qiáng)度最大,沿著壓力腔徑向位置氣旋強(qiáng)度逐漸減少,引起的壓力波動(dòng)逐漸減少,這里同2.2節(jié)的分析是一致的。在供氣壓力更大的情況下,氣旋強(qiáng)度逐漸增大。隨著供氣壓力增大,氣旋首先形成的位置差別不大,但是氣旋所在區(qū)域的中心與邊緣速度差、壓力差明顯增大,且供氣壓力的增大,輸入系統(tǒng)的能量增加,在氣腔結(jié)構(gòu)不變的情況下,氣體的流動(dòng)狀態(tài)變得紊亂,氣體的最大流速逐漸增大。

    圖8 不同供氣壓力下的瞬態(tài)速度和壓力云圖Fig.8 Transient Speed and Pressure Cloud Diagrams at Different Supply Pressures

    3.3 不同壓力腔形狀下的穩(wěn)態(tài)分析

    為了研究不同壓力腔形狀對(duì)氣旋的影響規(guī)律,選用如圖9所示的壓力腔模型,兩種模型僅壓力腔形狀不同,其具體形狀和參數(shù)如圖所示。

    圖9 壓力腔模型Fig.9 Pressure Chamber Model

    通過繪制不同壓力腔形狀下節(jié)流孔對(duì)稱中心沿徑向位置的壓力和湍動(dòng)能變化曲線,如圖10所示,兩種壓力腔內(nèi)都有氣旋的產(chǎn)生,錐形壓力腔內(nèi)的氣旋所造成的壓力降明顯最大,且氣旋面積相比更廣。而圓角壓力腔起到了過渡緩沖的作用,不僅減緩壓力腔內(nèi)的速度增加,而且也能降低壓力降的大小。

    圖10 不同壓力腔形狀壓力沿徑向位置分布曲線圖Fig.10 Distribution Curve of Pressure Along Radial Pposition under Different Pressure Chamber Shape

    不同壓力腔形狀時(shí)的流體的湍動(dòng)能分布曲線,如圖11所示??梢钥吹綀A角節(jié)流孔的湍動(dòng)能在沿壓力腔徑向分布中明顯最小,錐形節(jié)流孔產(chǎn)生的能量積聚對(duì)壓力腔垂直面產(chǎn)生的沖擊更大,氣旋強(qiáng)度更大。

    圖11 不同壓力腔形狀湍動(dòng)能沿徑向位置分布曲線圖Fig.11 Distribution Curve of Turbulent Kinetic Energy Along Radial Position under Different Pressure Chamber Shape

    3.4 不同壓力腔形狀下的瞬態(tài)分析

    通過對(duì)比圓角形和錐形壓力腔的不同時(shí)刻速度分布圖,如圖12所示可以看到圓角壓力腔由于圓角的存在在節(jié)流孔處的氣流過渡效果總是優(yōu)于錐形壓力腔。且錐形壓力腔內(nèi)氣旋形成面積區(qū)域大于圓角氣旋形成面積區(qū)域。且圓角壓力腔的最大氣流速度明顯小于錐形壓力腔的最大氣流速度。

    圖12 不同壓力腔形狀的瞬態(tài)速度云圖Fig.12 Transient Speed Cloud Diagrams at Different Pressure Chamber Shape

    通過對(duì)比穩(wěn)態(tài)和瞬態(tài)的仿真結(jié)果,因此在壓力腔的設(shè)計(jì)過程中,選擇合適的弧度進(jìn)行過渡,一方面可以減少阻力,另一方面可以使氣體流動(dòng)更加通暢,減小氣旋強(qiáng)度降低微振動(dòng)。

    4 結(jié)論

    利用STAR-ccm+軟件對(duì)空氣靜壓軸承的內(nèi)部流場進(jìn)行了仿真分析,得到以下主要結(jié)論:(1)通過利用沖擊射流的氣流流動(dòng)結(jié)構(gòu)特征劃分了氣流流場區(qū)域,利用大渦模擬研究了不同區(qū)域的氣流瞬時(shí)流動(dòng)情況,研究了氣旋的產(chǎn)生和變化規(guī)律。(2)從沖量原理、壓力變化等方面確定了氣旋導(dǎo)致微振動(dòng)的機(jī)理,即產(chǎn)生于阻滯區(qū)域、壓力腔徑向流動(dòng)區(qū)的氣旋以及在壁面射流區(qū)中的壓力腔出口處形成湍動(dòng)能和壓力波動(dòng)。(3)數(shù)值分析對(duì)比不了同供氣壓力和不同壓力腔形狀下計(jì)算域穩(wěn)態(tài)和瞬態(tài)流場,理論分析結(jié)合仿真結(jié)果表明,隨著供氣壓力的增加,氣旋強(qiáng)度增加,不同形狀的壓力腔對(duì)氣旋強(qiáng)度影響也較大。因此通過數(shù)值模擬氣旋強(qiáng)度,可以選擇較優(yōu)的壓力腔形狀。

    猜你喜歡
    大渦氣膜氣旋
    T 型槽柱面氣膜密封穩(wěn)態(tài)性能數(shù)值計(jì)算研究
    2003年3月北大西洋上兩個(gè)爆發(fā)性氣旋的“吞并”過程及發(fā)展機(jī)制分析?
    氨還原劑噴槍中保護(hù)氣旋流噴射作用的研究
    能源工程(2021年1期)2021-04-13 02:05:56
    氣膜孔堵塞對(duì)葉片吸力面氣膜冷卻的影響
    靜葉柵上游端壁雙射流氣膜冷卻特性實(shí)驗(yàn)
    基于壁面射流的下?lián)舯┝鞣欠€(wěn)態(tài)風(fēng)場大渦模擬
    軸流風(fēng)機(jī)葉尖泄漏流動(dòng)的大渦模擬
    躲避霧霾天氣的氣膜館
    北太平洋上一個(gè)爆發(fā)性氣旋族的結(jié)構(gòu)分析?
    2014年3月大西洋上一個(gè)爆發(fā)性氣旋的研究
    av视频在线观看入口| 国产成年人精品一区二区| 国产精品国产三级专区第一集| 色5月婷婷丁香| 午夜久久久久精精品| 午夜福利视频1000在线观看| 91久久精品国产一区二区三区| 热99在线观看视频| 91午夜精品亚洲一区二区三区| 日日啪夜夜撸| 国语自产精品视频在线第100页| 免费人成在线观看视频色| 久久草成人影院| 大香蕉97超碰在线| 高清av免费在线| 欧美日本亚洲视频在线播放| av福利片在线观看| 日韩一区二区视频免费看| 色视频www国产| videossex国产| 深爱激情五月婷婷| 99久久中文字幕三级久久日本| 国产精品国产三级专区第一集| 在线a可以看的网站| 亚洲av中文字字幕乱码综合| 午夜福利成人在线免费观看| 精品久久久久久久末码| 国产成人精品婷婷| 日产精品乱码卡一卡2卡三| 欧美日本视频| 久久久a久久爽久久v久久| 亚洲精品一区蜜桃| 最近的中文字幕免费完整| 校园人妻丝袜中文字幕| 我的女老师完整版在线观看| 婷婷色麻豆天堂久久 | 国语自产精品视频在线第100页| 日本一本二区三区精品| 少妇熟女欧美另类| 一个人观看的视频www高清免费观看| av又黄又爽大尺度在线免费看 | 欧美精品国产亚洲| 日本黄大片高清| 欧美日韩国产亚洲二区| 天堂影院成人在线观看| 韩国av在线不卡| 能在线免费看毛片的网站| 亚洲精品国产成人久久av| 男女国产视频网站| 在线观看一区二区三区| 国产免费福利视频在线观看| 波多野结衣巨乳人妻| 国产精品久久电影中文字幕| 日韩高清综合在线| 欧美性猛交黑人性爽| 99久久精品热视频| 天天躁夜夜躁狠狠久久av| 国产探花极品一区二区| 久久精品久久久久久久性| 免费观看性生交大片5| 国产欧美另类精品又又久久亚洲欧美| 五月伊人婷婷丁香| 精品久久国产蜜桃| 国产成人a∨麻豆精品| 免费在线观看成人毛片| 久久人妻av系列| 两性午夜刺激爽爽歪歪视频在线观看| 免费观看在线日韩| 五月伊人婷婷丁香| 国产亚洲午夜精品一区二区久久 | 男人狂女人下面高潮的视频| 中国美白少妇内射xxxbb| 激情 狠狠 欧美| 亚洲怡红院男人天堂| 精品人妻偷拍中文字幕| 网址你懂的国产日韩在线| 免费人成在线观看视频色| 亚洲久久久久久中文字幕| 久99久视频精品免费| 国产免费视频播放在线视频 | 国产精品av视频在线免费观看| 男女啪啪激烈高潮av片| 黄色日韩在线| 免费av观看视频| 国产又黄又爽又无遮挡在线| 亚洲国产精品国产精品| 桃色一区二区三区在线观看| 日本一二三区视频观看| 一夜夜www| 久久久精品94久久精品| 国产精品国产三级国产av玫瑰| 国产亚洲5aaaaa淫片| 亚洲欧美精品自产自拍| av线在线观看网站| 乱系列少妇在线播放| 欧美激情在线99| 2022亚洲国产成人精品| 久久久久久久久久成人| 亚洲av电影不卡..在线观看| av女优亚洲男人天堂| 嫩草影院入口| 成人漫画全彩无遮挡| 亚洲精品色激情综合| www日本黄色视频网| 黄色欧美视频在线观看| 亚洲国产成人一精品久久久| 热99re8久久精品国产| 成人无遮挡网站| 校园人妻丝袜中文字幕| 久久精品综合一区二区三区| 亚洲综合色惰| 97人妻精品一区二区三区麻豆| 我的老师免费观看完整版| 国产白丝娇喘喷水9色精品| 91狼人影院| 成人综合一区亚洲| 亚洲精品一区蜜桃| 特大巨黑吊av在线直播| 欧美3d第一页| 久99久视频精品免费| 性色avwww在线观看| 欧美激情在线99| 美女xxoo啪啪120秒动态图| 高清av免费在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | av免费观看日本| 久久久久久九九精品二区国产| 最近中文字幕2019免费版| 99久久九九国产精品国产免费| 日韩欧美在线乱码| 十八禁国产超污无遮挡网站| 久热久热在线精品观看| 国产乱来视频区| 在线免费十八禁| 一本一本综合久久| 一边摸一边抽搐一进一小说| 亚洲国产精品成人综合色| 色播亚洲综合网| 国产黄片视频在线免费观看| 精品免费久久久久久久清纯| 国产精品女同一区二区软件| 激情 狠狠 欧美| 听说在线观看完整版免费高清| 国产亚洲最大av| 床上黄色一级片| 亚洲av.av天堂| 国产成人91sexporn| 亚洲欧洲国产日韩| 国产一区二区在线观看日韩| 成人三级黄色视频| 成人三级黄色视频| 26uuu在线亚洲综合色| 国产一区二区在线观看日韩| 久久亚洲精品不卡| 亚洲伊人久久精品综合 | 免费观看性生交大片5| 99国产精品一区二区蜜桃av| 两个人视频免费观看高清| 日韩av在线大香蕉| 午夜老司机福利剧场| 最近手机中文字幕大全| 亚洲av免费高清在线观看| 国产人妻一区二区三区在| 久久久国产成人精品二区| 国产免费男女视频| 搞女人的毛片| 久久久久久久久久成人| 欧美激情国产日韩精品一区| 黄色配什么色好看| 国产高清不卡午夜福利| 永久网站在线| 三级国产精品欧美在线观看| 国产高潮美女av| 蜜臀久久99精品久久宅男| 欧美日本亚洲视频在线播放| av国产免费在线观看| 日产精品乱码卡一卡2卡三| 国产亚洲5aaaaa淫片| 久久久午夜欧美精品| 亚洲精品日韩在线中文字幕| 午夜日本视频在线| 亚洲电影在线观看av| 人人妻人人看人人澡| 99久久成人亚洲精品观看| 特大巨黑吊av在线直播| 九九爱精品视频在线观看| 国产在线一区二区三区精 | 人人妻人人澡人人爽人人夜夜 | 九九爱精品视频在线观看| 又爽又黄a免费视频| 十八禁国产超污无遮挡网站| 大又大粗又爽又黄少妇毛片口| 精品久久久噜噜| 国产探花在线观看一区二区| 91精品一卡2卡3卡4卡| 久久亚洲国产成人精品v| 欧美97在线视频| 亚洲人成网站在线播| 亚洲人与动物交配视频| 欧美色视频一区免费| 国产真实乱freesex| 麻豆精品久久久久久蜜桃| 亚洲欧美中文字幕日韩二区| 亚洲欧美一区二区三区国产| 狂野欧美白嫩少妇大欣赏| 99热这里只有是精品在线观看| 色视频www国产| 少妇人妻精品综合一区二区| 非洲黑人性xxxx精品又粗又长| 最近的中文字幕免费完整| 国产黄色视频一区二区在线观看 | 男的添女的下面高潮视频| 自拍偷自拍亚洲精品老妇| 少妇被粗大猛烈的视频| 久久久成人免费电影| 国产一区二区在线av高清观看| 99九九线精品视频在线观看视频| 日韩av在线大香蕉| 日本一二三区视频观看| 非洲黑人性xxxx精品又粗又长| 99久国产av精品| 观看免费一级毛片| 在线播放无遮挡| 国产毛片a区久久久久| 永久网站在线| 日本色播在线视频| 国产精品一及| 18禁在线播放成人免费| 久久久久九九精品影院| 日本黄色片子视频| 麻豆一二三区av精品| 亚洲欧美日韩东京热| 精品不卡国产一区二区三区| 天天躁夜夜躁狠狠久久av| 久久亚洲精品不卡| 九九热线精品视视频播放| 久久韩国三级中文字幕| 三级毛片av免费| 色尼玛亚洲综合影院| 国产午夜精品一二区理论片| 久久久久久久久久久免费av| 亚洲精品亚洲一区二区| 久久精品久久久久久噜噜老黄 | 直男gayav资源| 天堂av国产一区二区熟女人妻| 黄色配什么色好看| 国产精品一区二区性色av| 国产精品无大码| 亚洲精品乱久久久久久| 乱人视频在线观看| 亚洲欧美清纯卡通| 欧美3d第一页| 国产免费又黄又爽又色| 日韩,欧美,国产一区二区三区 | 一边摸一边抽搐一进一小说| 热99re8久久精品国产| 99久久成人亚洲精品观看| av在线亚洲专区| 亚洲中文字幕日韩| 日本与韩国留学比较| 亚洲国产欧美在线一区| 最近中文字幕高清免费大全6| 在现免费观看毛片| 精品久久久久久成人av| 免费av观看视频| 国产白丝娇喘喷水9色精品| 国产一区二区在线观看日韩| a级毛片免费高清观看在线播放| 久久久久久久国产电影| 国产精品一区二区性色av| 青春草亚洲视频在线观看| 超碰97精品在线观看| h日本视频在线播放| 爱豆传媒免费全集在线观看| 国产精品久久久久久精品电影小说 | 国产av码专区亚洲av| 国产成年人精品一区二区| 日日干狠狠操夜夜爽| 在线播放国产精品三级| 波野结衣二区三区在线| 七月丁香在线播放| 久久久国产成人免费| 精品久久久久久久人妻蜜臀av| 免费看a级黄色片| 欧美激情在线99| 国产乱人偷精品视频| 亚洲电影在线观看av| 亚洲av成人av| 国产亚洲精品久久久com| 午夜视频国产福利| 国产av码专区亚洲av| 成人特级av手机在线观看| 国产黄色视频一区二区在线观看 | 男的添女的下面高潮视频| 亚洲精品日韩av片在线观看| 久久精品国产自在天天线| 欧美一区二区国产精品久久精品| 欧美性猛交╳xxx乱大交人| 舔av片在线| 午夜免费激情av| 成人二区视频| 国产精品野战在线观看| 三级男女做爰猛烈吃奶摸视频| 人人妻人人澡欧美一区二区| 亚洲av免费在线观看| 岛国在线免费视频观看| 国产中年淑女户外野战色| 精品久久久久久久久久久久久| 亚洲,欧美,日韩| 中文天堂在线官网| 美女脱内裤让男人舔精品视频| 2021天堂中文幕一二区在线观| 日韩三级伦理在线观看| 91久久精品国产一区二区成人| 老司机影院成人| 欧美一级a爱片免费观看看| 日本色播在线视频| 亚洲欧美日韩卡通动漫| 麻豆精品久久久久久蜜桃| 联通29元200g的流量卡| 国产熟女欧美一区二区| 五月伊人婷婷丁香| 一夜夜www| 亚洲国产精品成人综合色| 亚洲人成网站高清观看| 欧美又色又爽又黄视频| 好男人视频免费观看在线| 成人一区二区视频在线观看| av在线观看视频网站免费| 国产精品熟女久久久久浪| 午夜免费男女啪啪视频观看| 深夜a级毛片| 国产亚洲av嫩草精品影院| 尤物成人国产欧美一区二区三区| 免费观看精品视频网站| 亚洲国产高清在线一区二区三| 男女国产视频网站| 国产v大片淫在线免费观看| 精品不卡国产一区二区三区| 国产精品av视频在线免费观看| 麻豆成人午夜福利视频| 国产伦一二天堂av在线观看| 亚洲精品成人久久久久久| 我要搜黄色片| 亚洲av不卡在线观看| 三级经典国产精品| 高清av免费在线| 我要搜黄色片| 亚洲在线自拍视频| 自拍偷自拍亚洲精品老妇| 毛片一级片免费看久久久久| 午夜精品国产一区二区电影 | 91久久精品国产一区二区三区| 国产亚洲一区二区精品| 伦理电影大哥的女人| 成人国产麻豆网| 性色avwww在线观看| 色5月婷婷丁香| 亚洲国产精品专区欧美| 一本一本综合久久| 一个人看视频在线观看www免费| 美女国产视频在线观看| 亚洲最大成人手机在线| 午夜视频国产福利| 亚洲av中文av极速乱| 99久久中文字幕三级久久日本| 久久国内精品自在自线图片| 国产亚洲精品av在线| 国产精品av视频在线免费观看| 亚洲精品国产av成人精品| 亚洲丝袜综合中文字幕| 国产一区二区三区av在线| videossex国产| 国产午夜精品一二区理论片| 看非洲黑人一级黄片| 一边摸一边抽搐一进一小说| 国产色婷婷99| 97人妻精品一区二区三区麻豆| 亚洲综合精品二区| 成人美女网站在线观看视频| 久久99精品国语久久久| 嫩草影院新地址| 亚洲在线自拍视频| 亚洲欧美日韩无卡精品| 国产成人a∨麻豆精品| 亚洲在线自拍视频| 亚洲自拍偷在线| 99久久精品国产国产毛片| 日韩在线高清观看一区二区三区| 久久这里只有精品中国| 亚洲欧美日韩高清专用| 中文资源天堂在线| 七月丁香在线播放| 日韩 亚洲 欧美在线| 亚洲在久久综合| 最近最新中文字幕免费大全7| 简卡轻食公司| 日韩,欧美,国产一区二区三区 | 麻豆乱淫一区二区| 少妇高潮的动态图| 亚洲最大成人av| 国产综合懂色| 日韩欧美精品v在线| 亚洲欧美中文字幕日韩二区| 久久久久国产网址| 国产一区二区在线av高清观看| 啦啦啦啦在线视频资源| 日韩精品青青久久久久久| 嫩草影院新地址| 国产精品人妻久久久久久| 精品久久国产蜜桃| 国产日韩欧美在线精品| 婷婷六月久久综合丁香| 中文字幕免费在线视频6| 国产极品天堂在线| 在线免费十八禁| 99视频精品全部免费 在线| 一个人免费在线观看电影| 亚洲国产精品久久男人天堂| 国产成人精品久久久久久| 黑人高潮一二区| 在现免费观看毛片| 看黄色毛片网站| 免费播放大片免费观看视频在线观看 | 在现免费观看毛片| 欧美区成人在线视频| 亚洲国产精品成人综合色| 免费av观看视频| 国产黄a三级三级三级人| 91精品国产九色| 人体艺术视频欧美日本| 国产高清国产精品国产三级 | 亚洲三级黄色毛片| 99热网站在线观看| 最近2019中文字幕mv第一页| 亚洲国产高清在线一区二区三| 成人漫画全彩无遮挡| 中文字幕人妻熟人妻熟丝袜美| 久久这里只有精品中国| 国产精品久久久久久久久免| 白带黄色成豆腐渣| 亚洲av.av天堂| 国产片特级美女逼逼视频| 免费av观看视频| 精品国产三级普通话版| 亚洲最大成人手机在线| 久久久久免费精品人妻一区二区| 精华霜和精华液先用哪个| 99久久无色码亚洲精品果冻| 男女那种视频在线观看| 人人妻人人澡欧美一区二区| 国产成人精品婷婷| 成人二区视频| 国产女主播在线喷水免费视频网站 | 欧美日韩综合久久久久久| 国产高清国产精品国产三级 | 午夜视频国产福利| 国产不卡一卡二| 狂野欧美白嫩少妇大欣赏| 国产精品一区二区在线观看99 | 国产精品无大码| 搡老妇女老女人老熟妇| 久久久a久久爽久久v久久| 国产真实伦视频高清在线观看| 亚洲最大成人手机在线| 一个人看视频在线观看www免费| 国产高清视频在线观看网站| 18禁动态无遮挡网站| 精品人妻偷拍中文字幕| av国产久精品久网站免费入址| kizo精华| 亚洲18禁久久av| 午夜精品国产一区二区电影 | 日韩欧美三级三区| 嘟嘟电影网在线观看| 一区二区三区高清视频在线| 毛片一级片免费看久久久久| 欧美高清成人免费视频www| 国产精品国产三级国产av玫瑰| 嘟嘟电影网在线观看| 国产爱豆传媒在线观看| 国产精品久久视频播放| 麻豆成人av视频| 国产精品三级大全| 久久精品国产鲁丝片午夜精品| 国产精品福利在线免费观看| 国产老妇伦熟女老妇高清| 精品酒店卫生间| 舔av片在线| 免费人成在线观看视频色| 美女国产视频在线观看| 蜜桃亚洲精品一区二区三区| 中文字幕av成人在线电影| 亚洲精品国产av成人精品| 久久久久国产网址| 乱人视频在线观看| 一级毛片aaaaaa免费看小| 我的女老师完整版在线观看| 极品教师在线视频| 久久精品人妻少妇| 亚洲综合精品二区| 免费无遮挡裸体视频| 超碰av人人做人人爽久久| 秋霞伦理黄片| 免费不卡的大黄色大毛片视频在线观看 | 成人特级av手机在线观看| 嫩草影院精品99| 免费人成在线观看视频色| 七月丁香在线播放| 精品99又大又爽又粗少妇毛片| 亚洲精品,欧美精品| 国产精品一区www在线观看| 亚洲中文字幕日韩| 毛片一级片免费看久久久久| 国产免费男女视频| 国产亚洲最大av| kizo精华| 高清av免费在线| 亚洲三级黄色毛片| 久久久久久久久大av| 亚洲人成网站在线播| 国产 一区精品| 免费在线观看成人毛片| 国产女主播在线喷水免费视频网站 | 精品久久久久久电影网 | 深夜a级毛片| 国产乱来视频区| 男女啪啪激烈高潮av片| 老司机影院成人| 亚洲电影在线观看av| 欧美成人一区二区免费高清观看| 久久久久国产网址| 日本免费a在线| 国产亚洲精品av在线| 九九热线精品视视频播放| 免费人成在线观看视频色| 爱豆传媒免费全集在线观看| 久久久久久久久久久丰满| 婷婷色综合大香蕉| 久久久国产成人精品二区| av在线观看视频网站免费| 免费在线观看成人毛片| 亚洲精品亚洲一区二区| 精品人妻熟女av久视频| 边亲边吃奶的免费视频| 日本欧美国产在线视频| 欧美日韩一区二区视频在线观看视频在线 | 精品久久久久久久末码| 亚洲色图av天堂| 日韩av在线大香蕉| 九草在线视频观看| 女人十人毛片免费观看3o分钟| 两性午夜刺激爽爽歪歪视频在线观看| a级毛片免费高清观看在线播放| 免费一级毛片在线播放高清视频| 久久人人爽人人爽人人片va| 国产成人精品一,二区| 丰满人妻一区二区三区视频av| 中文字幕人妻熟人妻熟丝袜美| 精品人妻视频免费看| 大香蕉97超碰在线| 国产高清不卡午夜福利| 日日摸夜夜添夜夜爱| 看片在线看免费视频| 国产私拍福利视频在线观看| 国产大屁股一区二区在线视频| 国产私拍福利视频在线观看| 国产亚洲av片在线观看秒播厂 | 国产女主播在线喷水免费视频网站 | АⅤ资源中文在线天堂| 国产精品电影一区二区三区| 精品久久国产蜜桃| 午夜免费男女啪啪视频观看| 久久久精品欧美日韩精品| 久久亚洲国产成人精品v| 中文字幕熟女人妻在线| 日韩成人伦理影院| 美女xxoo啪啪120秒动态图| 欧美人与善性xxx| 亚洲欧美日韩卡通动漫| 久久久久久久久久久免费av| 精品久久久久久久末码| 熟女电影av网| 国产精品嫩草影院av在线观看| 麻豆久久精品国产亚洲av| 简卡轻食公司| 欧美又色又爽又黄视频| 欧美3d第一页| 色噜噜av男人的天堂激情| 99九九线精品视频在线观看视频| 色5月婷婷丁香| 久久6这里有精品| 免费av不卡在线播放| 免费观看的影片在线观看| 麻豆久久精品国产亚洲av| 国产精品美女特级片免费视频播放器| 国产爱豆传媒在线观看| 免费在线观看成人毛片| 成人三级黄色视频| 亚洲国产精品久久男人天堂| 日韩大片免费观看网站 | 精品一区二区三区人妻视频| 国产精品一区www在线观看| 日日啪夜夜撸| 99久国产av精品| 久久久成人免费电影| 国产日韩欧美在线精品| 丰满乱子伦码专区| 搞女人的毛片| 亚洲五月天丁香| 亚洲av日韩在线播放| 国产片特级美女逼逼视频| 波多野结衣高清无吗| 99久国产av精品| 久久99热这里只有精品18|