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

    狹縫節(jié)流空氣軸承-轉子系統(tǒng)穩(wěn)定性分析

    2022-09-21 08:28:44于賀春張國慶岳金珂王文博王仁宗張素香時金虎
    機床與液壓 2022年17期
    關鍵詞:氣膜瞬態(tài)承載力

    于賀春,張國慶,岳金珂,王文博,王仁宗,張素香,時金虎

    (中原工學院機電學院,河南鄭州 450007)

    0 前言

    空氣靜壓主軸具有回轉精度高、能耗低、溫升小等優(yōu)點,是現代機床的核心部件之一。隨著加工精度、加工效率等要求的不斷提高,現代加工制造對空氣靜壓主軸的精度和轉速要求也不斷提高。

    軸承-轉子系統(tǒng)的研究常使用線性攝動法和軌跡法。攝動法無法展現轉子在整個系統(tǒng)中的運動軌跡,只適用于轉子在平衡狀態(tài)下的微小擾動;而軌跡法是通過對動力學方程和流體方程直接求解,得到轉子軸心軌跡,進而對整個系統(tǒng)穩(wěn)定性進行分析。CASTELLI和ELROD首先提出了軌跡法,通過研究長徑軸承,得到不同系統(tǒng)工況下轉子軸心軌跡變化,并將分析得到的數值與攝動法的計算結果相比較。LUND等為了研究軸承的動態(tài)特性,提出八系數法,為軸承-轉子系統(tǒng)的動態(tài)特性分析提供了便利,將8個不同方位的剛度阻尼系數作為軸承的動態(tài)參數表征軸承的穩(wěn)定性。韓東江等采用軸心軌跡、分岔圖以及頻譜圖等方式,呈現了系統(tǒng)不同振動下的動力學特征,給出了定性的理論分析以及相應的控制方法。張永芳等通過組裝技術獲得固定瓦-可傾瓦動壓氣體軸承非線性氣膜壓力的分布。HASSINI和 ARGHIR認為旋轉狀態(tài)下的轉子在其平衡位置不斷攝動形成其軸心軌跡的變化,每一個位置的氣膜力都可使其產生線性的動態(tài)變化,采用連續(xù)逼近法對軸承-轉子系統(tǒng)的運動軌跡進行分析。CHEN和WANG分別利用有限差分法和龍格-庫塔方法求得軸承的瞬態(tài)雷諾方程和系統(tǒng)動力學方程,通過計算發(fā)現系統(tǒng)會隨著轉子速度與質量的變化出現倍周期現象。MIYANAGA 和TOMIOKA利用線性攝動理論,研究了人字槽氣體軸承彈性基座的支承剛度和阻尼的變化,得出彈性基座可以提高軸承穩(wěn)定性的結論。陳東菊等通過擾動法對軸承動態(tài)特性系數進行求解,將結果引入建立的軸承-轉子動態(tài)振動模型中,求得了主軸的振動誤差曲線;通過對數據進行頻域處理,得到非線性對主軸徑向誤差的影響,并通過實驗對結果進行分析。任佟等人通過微擾動法和分離變量法對氣體軸承的氣膜阻尼系數建模、求解并分析,將氣膜阻尼系數引入軸承轉子系統(tǒng)進行動力學建模及計算,并將實驗與計算結果進行對比。張強和張霞妹通過研究軸承-轉子系統(tǒng),計算了軸承承載力、臨界轉速以及振型,并通過改變軸承的平均間隙和偏心率對轉子系統(tǒng)的臨界轉速進行研究。

    綜上所述,目前多采用線性方法或攝動方法結合剛度阻尼系數對系統(tǒng)的穩(wěn)定性進行分析,很少考慮氣膜力與轉子非線性運動的相互影響。為探究不同狹縫軸承參數對系統(tǒng)穩(wěn)定性的影響規(guī)律,本文作者利用Gambit軟件建立徑向支撐的雙排連續(xù)性狹縫軸承-轉子系統(tǒng)的動態(tài)耦合模型,考慮軸承與轉子之間的非線性影響。構造系統(tǒng)的穩(wěn)態(tài)環(huán)境,對穩(wěn)態(tài)下的軸心軌跡變化和氣體流場進行分析;對系統(tǒng)載荷方向的軸心軌跡變化進行瞬時響應分析,探究不同軸承參數下系統(tǒng)受到載荷作用后的轉子位移在載荷方向的變化曲線,從而利用時域參數分析不同軸承參數下系統(tǒng)的穩(wěn)定性。

    1 系統(tǒng)耦合方法分析與模型建立

    1.1 流固耦合控制方程

    主軸內的軸承和轉子在工作過程中相互影響,轉子自身的慣性力和氣膜慣性力會使轉子產生非線性運動,轉子軸心軌跡的變化會使氣膜流場分布特性發(fā)生改變,因此可以利用雙向流固耦合方法對狹縫節(jié)流徑向軸承支撐的轉子系統(tǒng)進行求解。圖1所示為系統(tǒng)耦合計算流程。

    圖1 系統(tǒng)耦合計算原理

    通過聯立流體方程與固體方程,可以得到系統(tǒng)的耦合控制方程。流固耦合遵循基本守恒定律,在流固耦合交界面上,流體f與固體s之間變量的傳遞保持守恒,耦合方程滿足公式(1):

    (1)

    式中:為應力,Pa;為位移,m;為熱流量,W;為溫度,K。

    1.2 動網格技術

    動態(tài)網格(動網格)技術是一項網格更新技術,可以改變流場區(qū)域內網格體積大小以及網格的出現和產生,并控制網格質量。流場中的任意網格體積上的一般標量的守恒方程為

    (2)

    式中:為體積,m;?為體積運動邊界;為網格運動速度,m/s;為流體密度,kg/m;為流體速度矢量,m/s;為邊界截面積,m;為耗散系數;為標量的源項。

    1.3 系統(tǒng)耦合參數確定

    以雙排狹縫節(jié)流徑向空氣軸承支撐的轉子系統(tǒng)為研究對象,如圖2所示,軸承參數如表1所示。

    圖2 耦合模型系統(tǒng)簡圖

    表1 狹縫軸承部分參數 單位:mm

    在同一坐標系下對系統(tǒng)進行簡化,軸向為模型的方向,徑向受載荷力方向為方向。轉子模型軸頸長60 mm,軸頸表面與氣膜內表面重合。

    1.4 網格劃分及定義邊界條件

    (1)網格劃分

    在耦合計算過程中,由于流場自身的變化和轉子的運動會使氣膜間隙處網格產生大的變形和重構。在劃分網格過程中,要在保證網格計算精度的基礎上,減少網格數量,且與轉子耦合處的網格一一對應。不同流場區(qū)域網格節(jié)點數如表2所示。

    表2 流場網格節(jié)點數

    利用Gambit軟件對簡化后的流場模型進行網格劃分。由于軸承長度與氣膜間隙尺寸相差較大,為保證網格質量,將流場分為氣膜間隙區(qū)和狹縫區(qū),如圖3所示??芍簹饽らg隙區(qū)網格密度由狹縫至氣體出口方向逐漸減小;狹縫區(qū)網格密度由狹縫入口到狹縫出口逐漸增大。

    圖3 流場網格劃分

    (2)定義邊界條件

    定義邊界條件:①氣體在運動過程中保持溫度恒定;②根據文獻[17]的計算結果,文中采用層流方式對軸承進行分析;③忽略加工誤差對壁面的影響,轉子表面絕對光滑;④系統(tǒng)供氣壓力為0.5 MPa,軸承的出口邊界為壓力出口,出口壓力與環(huán)境壓力相等;⑤設置流場內部表面和轉子外表面為耦合面,氣膜間隙兩側出氣表面定義為變形壁面;⑥對轉子施加軸向位移約束,排除轉子軸向竄動對結果的干擾。

    2 狹縫節(jié)流軸承-轉子系統(tǒng)耦合計算

    2.1 系統(tǒng)流固耦合參數設置

    利用ANSYS Workbench中System Coupling模塊對流體域和固體域的計算結果進行數據交換。根據流體網格大小對計算步長進行調整,設置耦合計算步長為5×10s。

    2.2 系統(tǒng)流固耦合穩(wěn)態(tài)環(huán)境準備

    系統(tǒng)的進氣孔壓力為0.5 MPa,轉子轉速為10 000 r/min,耦合作用下的轉子軌跡如圖4所示??芍恨D子在自身慣性力和氣膜力的作用下逐漸趨于某一穩(wěn)定位置;初始時,軸承與轉子之間的氣流流動較為混亂,轉子受氣流沖擊產生不規(guī)則運動,當氣體流動穩(wěn)定且轉速增加到一定穩(wěn)定值后,轉子的運動趨于穩(wěn)定,最終在這一位置附近進行渦動。

    圖4 耦合作用下的轉子軌跡

    圖5所示為轉子運動穩(wěn)定后,氣膜間隙流場分布情況。此時轉子不受外界因素影響,只在自身慣性力和軸承氣膜力作用下運動,氣膜內流場分布均勻,轉子運動較為平穩(wěn)。

    圖5 穩(wěn)態(tài)環(huán)境下氣體壓力分布

    2.3 系統(tǒng)的位移分析

    圖6所示為狹縫軸承-轉子系統(tǒng)受到階躍載荷后的轉子位移和承載力變化曲線。系統(tǒng)在穩(wěn)態(tài)環(huán)境下運動,當=2.5 ms時,轉子在方向受到階躍載荷作用,轉子運動軌跡發(fā)生改變,氣膜承載力隨之發(fā)生改變,轉子在載荷作用下到達新的平衡位置。

    圖6 耦合作用下的轉子位移及承載力變化關系

    由圖6(a)中可知:轉子位移量增加至最大值7.8 μm后減小,最終趨于平緩,穩(wěn)定在7.1 μm左右;氣膜承載力增大至465 N后減小,趨于平緩后穩(wěn)定在359 N左右,位移最大值與承載力最大值同時出現;當>3.75 ms時,轉子的位移變化和氣膜承載力變化趨于平穩(wěn),系統(tǒng)達到穩(wěn)定狀態(tài)。由圖6(b)可知:系統(tǒng)穩(wěn)定狀態(tài)下轉子的位移變化和氣膜承載力變化相反;當轉子位移增加時,氣膜間隙減小,氣膜承載力增加,反之亦然,此時系統(tǒng)狀態(tài)為動態(tài)平衡。

    2.4 系統(tǒng)的流場分析

    圖7所示為轉子在移動過程中軸承內的氣體壓力分布情況??芍?當=2.5 ms時,轉子未受載荷作用,流場壓力在同一軸向位置的圓周方向均勻分布,轉子在自身慣性和氣膜力作用下處于平衡狀態(tài);當=2.75 ms時,轉子受到方向上350 N階躍載荷,轉子在載荷作用下與狹縫軸承產生偏心,隨著轉子位移的不斷增加,氣膜的高壓面積不斷增加,承載力也不斷增加;當=3.0 ms時,轉子的位移量最大,氣膜高壓面積達到最大,此時氣膜承載力大于載荷,轉子向反方向移動;當>3.75 ms后,載荷與承載力達到平衡,氣膜的高壓區(qū)面積不再發(fā)生明顯變化,系統(tǒng)處于新的動態(tài)平衡。

    圖7 耦合作用下的氣體壓力分布

    2.5 狹縫節(jié)流軸承-轉子系統(tǒng)時域分析

    利用時域分析研究轉子在載荷作用下的運動規(guī)律,在瞬態(tài)響應階段,狹縫軸承-轉子系統(tǒng)受到階躍載荷作用,發(fā)生振蕩,經過一段時間,系統(tǒng)達到新的穩(wěn)態(tài)狀態(tài)。

    轉子在載荷作用下運動,一定時間內,轉子位移首次達到穩(wěn)態(tài)的時間為上升時間;轉子位移首次達到最大值的時間為峰值時間;上升時間與峰值時間的大小反映了系統(tǒng)的快速性。在二階欠阻尼系統(tǒng)中,最大峰值與穩(wěn)態(tài)值差值與穩(wěn)態(tài)值的比值為系統(tǒng)最大超調量,利用式(3)可得出最大超調量和阻尼比之間的關系:

    (3)

    阻尼比決定系統(tǒng)的振蕩特性,越大,系統(tǒng)穩(wěn)定性越大。越小,值越大,振蕩幅值越小,當增大到一定程度時,系統(tǒng)的響應曲線單調上升。

    (1)二階欠阻尼系統(tǒng)

    圖8所示為在=2.5 ms時施加階躍載荷作用下狹縫間距為30 mm的轉子的位移變化曲線。可知:當>2.5 ms時,轉子在方向上不受載荷作用;轉子在方向上的位移分為2個階段:瞬態(tài)響應階段和穩(wěn)態(tài)響應階段,其波動規(guī)律符合二階欠阻尼系統(tǒng)響應。

    圖8 t=2.5 ms施加階躍載荷轉子在x、y方向位移曲線

    (2)二階過阻尼系統(tǒng)

    圖9所示為在=2 ms時施加階躍載荷作用下狹縫間距為10 mm的轉子的位移變化曲線。可知:轉子在方向的位移增長率絕對值逐漸減小,位移變化曲線為二階過阻尼系統(tǒng)響應曲線,此時>1,系統(tǒng)持續(xù)處于瞬態(tài)響應階段,不發(fā)生振蕩,且無最大超調量。

    圖9 t=2 ms施加階躍載荷轉子在x、y方向位移曲線

    3 耦合計算結果與分析

    3.1 狹縫間距對系統(tǒng)瞬態(tài)參數的影響

    在系統(tǒng)供氣壓力和所受載荷(在=2.0 ms時施加)恒定的情況下,依次改變狹縫軸承的狹縫間距為10、20、30、40、50 mm,計算承載力,結果如圖10所示。

    圖10 瞬態(tài)響應下承載力變化(不同狹縫間距)

    由圖10可知:承載力峰值和承載力的變化率隨狹縫間距增加逐漸減小,承載力到達峰值的時間隨著狹縫間距增加而逐漸減小。

    圖11所示為不同狹縫間距下轉子位移變化網格圖。

    圖11 瞬態(tài)響應下y方向位移變化(不同狹縫間距)

    由圖11可知:轉子位移穩(wěn)定值隨著狹縫間距的增加,先減小后增加,即相同供氣壓力下,轉子位移變化反映了系統(tǒng)的承載力變化,即隨著狹縫間距增加,系統(tǒng)承載力先增加后下降;狹縫間距在20、30、40 mm時,轉子位移先增大后減小,最終趨于穩(wěn)定,而狹縫間距在10、50 mm時,轉子位移先增加后減少,未達到穩(wěn)定狀態(tài);當狹縫間距為10、50 mm時,系統(tǒng)的位移變化曲線為二階過阻尼系統(tǒng)響應曲線,此狀態(tài)下系統(tǒng)瞬態(tài)響應時間增加,系統(tǒng)無振蕩,阻尼比>1。狹縫間距在20、30、40 mm時,轉子位移曲線為二階欠阻尼響應曲線。

    3.2 狹縫深度對系統(tǒng)瞬態(tài)參數的影響

    系統(tǒng)其余參數不變,依次改變軸承的狹縫深度為3、5、7、9、11 mm,軸承承載力隨時間變化曲線如圖12所示??芍弘S著狹縫深度增加,承載力最大值逐漸下降,即高壓氣體經過狹縫節(jié)流器,與周圍壁面接觸,流入氣膜間隙后,能量隨著接觸壁面的增加而逐漸減小,導致承載力在瞬態(tài)響應期間的變化不同,當系統(tǒng)達到穩(wěn)態(tài)后,系統(tǒng)所受外界的載荷相同,軸承的承載力相同。

    圖12 瞬態(tài)響應下承載力變化(不同狹縫深度)

    圖13所示為不同狹縫深度下的轉子位移變化網格圖??芍弘S著狹縫深度的不斷增加,轉子位移先增加后減小,由于不同狹縫深度下的氣膜承載能力不同,導致轉子在瞬態(tài)響應階段的位移變化量不同,且最終轉子的穩(wěn)定位置也不相同。

    圖13 瞬態(tài)響應下y方向位移變化(不同狹縫深度)

    不同狹縫深度下系統(tǒng)時域響應參數如表3所示??芍弘S著狹縫深度逐漸增加,轉子位移的上升時間和峰值時間逐漸增加,最大超調量逐漸下降。因此,增加狹縫深度會使系統(tǒng)的響應快速性下降。狹縫深度增加使系統(tǒng)增大,系統(tǒng)振蕩特性減小,系統(tǒng)穩(wěn)定性提高。

    表3 不同狹縫深度下系統(tǒng)瞬態(tài)響應參數

    3.3 狹縫寬度對系統(tǒng)瞬態(tài)參數的影響

    系統(tǒng)其余參數保持不變,改變狹縫寬度,系統(tǒng)承載力變化如圖14所示:狹縫寬度增加會使軸承的承載力發(fā)生改變,隨著狹縫寬度的不斷增加,波動過程中的軸承承載力最大值逐漸增加,承載力最大值與穩(wěn)態(tài)值的比值不斷增加。

    圖14 瞬態(tài)響應下承載力變化(不同狹縫寬度)

    圖15所示為瞬態(tài)響應過程中,轉子在載荷方向下的位移變化網格圖??芍寒敧M縫寬度為0.008 mm時,轉子位移一直增加,前4 ms內未達到穩(wěn)態(tài),系統(tǒng)未產生波動;狹縫寬度在0.01~0.016 mm時,轉子位移有最大值,且隨著寬度的增加,位移絕對值的最大值逐漸增加。穩(wěn)定狀態(tài)下,軸承位移的絕對值隨著狹縫寬度的增加先減小后增加。

    圖15 瞬態(tài)響應下載荷方向位移變化(不同狹縫寬度)

    狹縫寬度為0.008 mm時,轉子位移沒有最大值,為過阻尼系統(tǒng),>1。表4所示為狹縫寬度為0.01~0.016 mm的瞬態(tài)響應參數變化??芍荷仙龝r間和峰值時間在狹縫寬度為0.012 mm時最小,此時位移峰值最大,位移最大值和穩(wěn)態(tài)值的比值最大,系統(tǒng)的反應速度最快;超調量先增加后減少,即先減小后增加,<1。因此,隨著狹縫寬度的增加,系統(tǒng)先減小后增加,系統(tǒng)穩(wěn)定性先減小后增加。

    表4 不同狹縫寬度下系統(tǒng)瞬態(tài)響應參數

    4 結論

    (1)通過流固耦合計算方法可以得到系統(tǒng)動態(tài)耦合下軸承的流場壓力分布與轉子的軸心軌跡變化,可知系統(tǒng)轉子軌跡在穩(wěn)態(tài)時呈不斷波動狀態(tài),流場壓力分布變化與轉子運動軌跡相互影響。

    (2)通過時域分析法可以得出,不同系統(tǒng)對相同載荷下的穩(wěn)態(tài)和瞬態(tài)反應各不相同,不同狹縫參數下的位移變化曲線有2種表現形式:欠阻尼系統(tǒng)和過阻尼系統(tǒng)。在一定范圍內,增加兩狹縫之間的距離可使系統(tǒng)阻尼比和穩(wěn)定性先增加后減?。辉黾营M縫深度可使系統(tǒng)阻尼比和穩(wěn)定性逐漸增加;增加狹縫寬度可使系統(tǒng)阻尼比和穩(wěn)定性先減小后增加。系統(tǒng)穩(wěn)定性增加,系統(tǒng)對外界反應的快速性會下降。

    猜你喜歡
    氣膜瞬態(tài)承載力
    T 型槽柱面氣膜密封穩(wěn)態(tài)性能數值計算研究
    氣膜孔堵塞對葉片吸力面氣膜冷卻的影響
    高壓感應電動機斷電重啟時的瞬態(tài)仿真
    防爆電機(2020年3期)2020-11-06 09:07:36
    靜葉柵上游端壁雙射流氣膜冷卻特性實驗
    火箭推進(2020年2期)2020-05-06 02:53:56
    躲避霧霾天氣的氣膜館
    CFRP-PCP板加固混凝土梁的抗彎承載力研究
    十億像素瞬態(tài)成像系統(tǒng)實時圖像拼接
    中國光學(2015年5期)2015-12-09 09:00:39
    基于瞬態(tài)流場計算的滑動軸承靜平衡位置求解
    DC/DC變換器中的瞬態(tài)特性分析
    耐火鋼圓鋼管混凝土柱耐火極限和承載力
    亚洲欧洲精品一区二区精品久久久 | 黑人猛操日本美女一级片| 青青草视频在线视频观看| 欧美国产精品va在线观看不卡| 久久久久久久大尺度免费视频| 亚洲丝袜综合中文字幕| 亚洲一区二区三区欧美精品| 午夜久久久在线观看| 精品卡一卡二卡四卡免费| 丝袜在线中文字幕| 国产成人精品福利久久| 人人妻人人澡人人看| 看免费av毛片| 亚洲国产欧美日韩在线播放| 久久亚洲国产成人精品v| a级毛片在线看网站| 在线观看三级黄色| 国产 精品1| 日本欧美视频一区| 全区人妻精品视频| 国产成人精品婷婷| 亚洲欧美一区二区三区国产| 久久国产精品男人的天堂亚洲 | 人妻少妇偷人精品九色| 精品国产一区二区三区四区第35| 午夜福利网站1000一区二区三区| 中文精品一卡2卡3卡4更新| 一级毛片我不卡| 自线自在国产av| 国产成人精品福利久久| 日日撸夜夜添| 高清欧美精品videossex| videos熟女内射| 亚洲成色77777| 最黄视频免费看| 黄色怎么调成土黄色| 青春草国产在线视频| 蜜桃在线观看..| a级毛片在线看网站| 有码 亚洲区| www日本在线高清视频| 在线观看人妻少妇| 在线看a的网站| 欧美 亚洲 国产 日韩一| 男人舔女人的私密视频| 插逼视频在线观看| 午夜激情久久久久久久| 夜夜骑夜夜射夜夜干| 伊人久久国产一区二区| 在线天堂最新版资源| 亚洲av男天堂| 人妻人人澡人人爽人人| 日韩视频在线欧美| 汤姆久久久久久久影院中文字幕| 久久久精品94久久精品| 欧美3d第一页| 色94色欧美一区二区| 少妇人妻 视频| 午夜福利视频精品| 又黄又爽又刺激的免费视频.| 99视频精品全部免费 在线| www.色视频.com| 午夜91福利影院| 亚洲少妇的诱惑av| 中文字幕另类日韩欧美亚洲嫩草| 一本—道久久a久久精品蜜桃钙片| 午夜91福利影院| videos熟女内射| 免费观看av网站的网址| 国产精品免费大片| 热re99久久精品国产66热6| 亚洲性久久影院| 国产一区二区三区av在线| 国产成人av激情在线播放| a级毛色黄片| 91精品三级在线观看| 日本午夜av视频| 国精品久久久久久国模美| 晚上一个人看的免费电影| 看免费成人av毛片| 日韩中文字幕视频在线看片| 一级黄片播放器| 日韩在线高清观看一区二区三区| 亚洲色图 男人天堂 中文字幕 | 久久狼人影院| 伊人亚洲综合成人网| 夫妻午夜视频| 美女视频免费永久观看网站| 国产成人精品福利久久| 亚洲一级一片aⅴ在线观看| 高清欧美精品videossex| 欧美97在线视频| 美女大奶头黄色视频| 一级毛片电影观看| 在线免费观看不下载黄p国产| 国产又色又爽无遮挡免| 免费播放大片免费观看视频在线观看| 一本大道久久a久久精品| 免费少妇av软件| 国产永久视频网站| 日本爱情动作片www.在线观看| 在线亚洲精品国产二区图片欧美| videosex国产| 久久人人爽人人爽人人片va| 亚洲av综合色区一区| 捣出白浆h1v1| 一本—道久久a久久精品蜜桃钙片| 色婷婷久久久亚洲欧美| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产极品粉嫩免费观看在线| 久久国内精品自在自线图片| 久久毛片免费看一区二区三区| 免费高清在线观看视频在线观看| 热re99久久国产66热| 亚洲第一区二区三区不卡| 少妇人妻久久综合中文| 熟女人妻精品中文字幕| 亚洲,欧美,日韩| 亚洲熟女精品中文字幕| 国产成人a∨麻豆精品| 十八禁高潮呻吟视频| 亚洲欧美精品自产自拍| 久久综合国产亚洲精品| 在线亚洲精品国产二区图片欧美| 纵有疾风起免费观看全集完整版| 啦啦啦中文免费视频观看日本| 久久精品国产亚洲av天美| 中文字幕另类日韩欧美亚洲嫩草| 男女边吃奶边做爰视频| 18禁观看日本| 天天操日日干夜夜撸| www.熟女人妻精品国产 | 精品一区二区三区四区五区乱码 | 欧美最新免费一区二区三区| 精品久久国产蜜桃| 中文字幕免费在线视频6| 欧美最新免费一区二区三区| 一级毛片电影观看| 亚洲av在线观看美女高潮| 狠狠婷婷综合久久久久久88av| 国产激情久久老熟女| 麻豆精品久久久久久蜜桃| 黄色配什么色好看| av免费在线看不卡| 国产乱来视频区| 亚洲精品美女久久av网站| 亚洲综合色惰| 国产精品一二三区在线看| 女性生殖器流出的白浆| 色哟哟·www| 最近的中文字幕免费完整| 亚洲久久久国产精品| 久久这里只有精品19| 久久国内精品自在自线图片| 国产精品国产av在线观看| 最黄视频免费看| 高清欧美精品videossex| 久久久欧美国产精品| 少妇被粗大猛烈的视频| 水蜜桃什么品种好| 夫妻性生交免费视频一级片| 欧美国产精品一级二级三级| av线在线观看网站| 这个男人来自地球电影免费观看 | 丰满迷人的少妇在线观看| 晚上一个人看的免费电影| 久久人妻熟女aⅴ| 国产高清国产精品国产三级| 久久精品国产a三级三级三级| 纵有疾风起免费观看全集完整版| 男男h啪啪无遮挡| 免费少妇av软件| 在线观看免费日韩欧美大片| 99久国产av精品国产电影| 日韩在线高清观看一区二区三区| 一本—道久久a久久精品蜜桃钙片| 亚洲五月色婷婷综合| 亚洲精品美女久久久久99蜜臀 | 在线免费观看不下载黄p国产| 久久亚洲国产成人精品v| 亚洲国产av影院在线观看| av女优亚洲男人天堂| 欧美+日韩+精品| 菩萨蛮人人尽说江南好唐韦庄| a级毛片黄视频| 99国产精品免费福利视频| 日日撸夜夜添| 久久人人爽人人片av| 久久久久久久久久久免费av| 国产高清国产精品国产三级| 99久久综合免费| 国产精品久久久av美女十八| 欧美精品人与动牲交sv欧美| 久久97久久精品| 国产日韩欧美在线精品| 视频在线观看一区二区三区| 亚洲精品乱码久久久久久按摩| 国产精品偷伦视频观看了| 极品少妇高潮喷水抽搐| 国产探花极品一区二区| 免费高清在线观看视频在线观看| 又粗又硬又长又爽又黄的视频| 美国免费a级毛片| 亚洲av电影在线观看一区二区三区| 欧美3d第一页| 成人免费观看视频高清| 亚洲av男天堂| 99热6这里只有精品| 欧美激情国产日韩精品一区| xxx大片免费视频| 欧美成人午夜精品| 深夜精品福利| 99热网站在线观看| 亚洲精品第二区| 精品人妻熟女毛片av久久网站| 精品一区二区免费观看| 国产乱来视频区| 久久久久久久久久久免费av| 久久人人爽人人爽人人片va| 男人舔女人的私密视频| 亚洲精品乱码久久久久久按摩| 高清不卡的av网站| 午夜影院在线不卡| 亚洲熟女精品中文字幕| a级毛片黄视频| 欧美激情 高清一区二区三区| 亚洲欧洲国产日韩| 亚洲精品第二区| 免费高清在线观看日韩| 中国国产av一级| 久久精品国产亚洲av天美| 午夜福利在线观看免费完整高清在| 亚洲欧洲精品一区二区精品久久久 | 十八禁高潮呻吟视频| 交换朋友夫妻互换小说| 大片免费播放器 马上看| 在线观看免费视频网站a站| 欧美日韩视频精品一区| 欧美日韩av久久| 国产白丝娇喘喷水9色精品| 一边亲一边摸免费视频| 午夜精品国产一区二区电影| 91午夜精品亚洲一区二区三区| 亚洲精品av麻豆狂野| 欧美97在线视频| 亚洲熟女精品中文字幕| 国产精品一区www在线观看| 五月开心婷婷网| 少妇的逼好多水| 18禁动态无遮挡网站| 国产成人精品久久久久久| 免费大片18禁| 精品一区二区三卡| 久久这里只有精品19| 国产黄色视频一区二区在线观看| 亚洲av电影在线进入| 国国产精品蜜臀av免费| 男女午夜视频在线观看 | 三上悠亚av全集在线观看| 亚洲av日韩在线播放| 综合色丁香网| 美女中出高潮动态图| 99香蕉大伊视频| 国产片特级美女逼逼视频| 人人妻人人爽人人添夜夜欢视频| 2018国产大陆天天弄谢| 极品少妇高潮喷水抽搐| 免费播放大片免费观看视频在线观看| 在线观看www视频免费| 中文字幕人妻丝袜制服| 99久久综合免费| 18+在线观看网站| 少妇人妻久久综合中文| 久热久热在线精品观看| 亚洲精品视频女| av国产久精品久网站免费入址| 国产不卡av网站在线观看| 亚洲av.av天堂| 久久久精品区二区三区| 精品人妻偷拍中文字幕| 99久久精品国产国产毛片| 亚洲精品美女久久久久99蜜臀 | 中国国产av一级| 亚洲欧美色中文字幕在线| 国产 一区精品| 国产老妇伦熟女老妇高清| 国产在线免费精品| 久久热在线av| 久久精品aⅴ一区二区三区四区 | 久久热在线av| 青春草视频在线免费观看| 婷婷色麻豆天堂久久| 午夜福利视频精品| 亚洲久久久国产精品| 啦啦啦在线观看免费高清www| 三上悠亚av全集在线观看| 汤姆久久久久久久影院中文字幕| a级毛片黄视频| 亚洲成人av在线免费| 久久综合国产亚洲精品| av视频免费观看在线观看| 久久久久久伊人网av| 超碰97精品在线观看| 纵有疾风起免费观看全集完整版| 国产在视频线精品| 不卡视频在线观看欧美| 欧美精品国产亚洲| 国产成人精品一,二区| 精品国产乱码久久久久久小说| 欧美成人午夜精品| 午夜日本视频在线| 亚洲国产精品成人久久小说| 欧美日韩综合久久久久久| av在线播放精品| a级片在线免费高清观看视频| 日韩三级伦理在线观看| 伊人久久国产一区二区| 国产国语露脸激情在线看| 三级国产精品片| 亚洲av日韩在线播放| av国产久精品久网站免费入址| 欧美日本中文国产一区发布| 美女xxoo啪啪120秒动态图| 久久狼人影院| 大码成人一级视频| 精品亚洲成国产av| 亚洲伊人久久精品综合| 成年av动漫网址| 久久99蜜桃精品久久| 国产精品久久久久久精品电影小说| 国产有黄有色有爽视频| 王馨瑶露胸无遮挡在线观看| 五月开心婷婷网| 日本av手机在线免费观看| 久久久精品免费免费高清| 国产精品国产三级国产av玫瑰| 三上悠亚av全集在线观看| 香蕉国产在线看| 曰老女人黄片| 久久99一区二区三区| av国产久精品久网站免费入址| 岛国毛片在线播放| 啦啦啦中文免费视频观看日本| 免费看不卡的av| 一级片'在线观看视频| tube8黄色片| 日韩不卡一区二区三区视频在线| 亚洲国产成人一精品久久久| 久久精品久久久久久噜噜老黄| 久久精品久久久久久久性| a级毛片在线看网站| 国产老妇伦熟女老妇高清| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩国产mv在线观看视频| 国产欧美亚洲国产| av国产精品久久久久影院| 久久亚洲国产成人精品v| 日本猛色少妇xxxxx猛交久久| 国产亚洲av片在线观看秒播厂| 日韩不卡一区二区三区视频在线| 99久久人妻综合| 看非洲黑人一级黄片| 热re99久久精品国产66热6| 超色免费av| 国产黄色视频一区二区在线观看| 亚洲精品乱码久久久久久按摩| 欧美日韩av久久| 亚洲精品aⅴ在线观看| 久久精品国产综合久久久 | 热99国产精品久久久久久7| 成人亚洲欧美一区二区av| 一二三四中文在线观看免费高清| 国产在线免费精品| 亚洲成人手机| 亚洲人成77777在线视频| 免费观看无遮挡的男女| av有码第一页| 99九九在线精品视频| 成人免费观看视频高清| av国产久精品久网站免费入址| 人人妻人人澡人人爽人人夜夜| 亚洲欧美日韩卡通动漫| 看十八女毛片水多多多| 毛片一级片免费看久久久久| 人体艺术视频欧美日本| 精品国产国语对白av| 看十八女毛片水多多多| 欧美日韩国产mv在线观看视频| 成人亚洲欧美一区二区av| 欧美日韩国产mv在线观看视频| 黑丝袜美女国产一区| 国产精品久久久久久av不卡| 波野结衣二区三区在线| 丰满少妇做爰视频| 国产有黄有色有爽视频| 亚洲精品日本国产第一区| 午夜福利在线观看免费完整高清在| av不卡在线播放| 天堂中文最新版在线下载| 夜夜骑夜夜射夜夜干| 国产av码专区亚洲av| 99九九在线精品视频| 精品国产露脸久久av麻豆| 亚洲国产最新在线播放| 日本免费在线观看一区| 在现免费观看毛片| 母亲3免费完整高清在线观看 | 在线亚洲精品国产二区图片欧美| 国产麻豆69| 乱人伦中国视频| 亚洲国产av影院在线观看| 九九爱精品视频在线观看| 免费高清在线观看视频在线观看| 看免费av毛片| 高清毛片免费看| 夫妻性生交免费视频一级片| 十分钟在线观看高清视频www| 国产精品99久久99久久久不卡 | 免费播放大片免费观看视频在线观看| 中文字幕av电影在线播放| 国产爽快片一区二区三区| 色94色欧美一区二区| 黑人高潮一二区| 一级黄片播放器| 久久久久久久久久久免费av| 一级片'在线观看视频| 性色av一级| 日韩熟女老妇一区二区性免费视频| 一本色道久久久久久精品综合| 日韩一区二区三区影片| 人人妻人人澡人人爽人人夜夜| 午夜激情久久久久久久| 99热这里只有是精品在线观看| 在线精品无人区一区二区三| 久久97久久精品| av有码第一页| 高清欧美精品videossex| 国产成人精品无人区| 国产精品一区www在线观看| 亚洲精品日韩在线中文字幕| 亚洲国产欧美日韩在线播放| 精品一区二区免费观看| 美女中出高潮动态图| 精品一区在线观看国产| 国产一区二区三区综合在线观看 | 一级片免费观看大全| av卡一久久| av有码第一页| 欧美人与性动交α欧美软件 | 国产无遮挡羞羞视频在线观看| 欧美bdsm另类| 亚洲国产毛片av蜜桃av| √禁漫天堂资源中文www| 伦精品一区二区三区| 久久97久久精品| 最近中文字幕2019免费版| 亚洲四区av| 午夜日本视频在线| 国产爽快片一区二区三区| 亚洲精品乱久久久久久| 80岁老熟妇乱子伦牲交| 观看美女的网站| 女性生殖器流出的白浆| 久久国内精品自在自线图片| 男女边摸边吃奶| 另类精品久久| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产毛片av蜜桃av| 中国国产av一级| 女性生殖器流出的白浆| 男人爽女人下面视频在线观看| 1024视频免费在线观看| 狠狠精品人妻久久久久久综合| 一本色道久久久久久精品综合| 考比视频在线观看| 久久 成人 亚洲| 99国产综合亚洲精品| 中文乱码字字幕精品一区二区三区| 亚洲av中文av极速乱| 亚洲精品国产av蜜桃| av免费在线看不卡| 国产av一区二区精品久久| 久久人人爽av亚洲精品天堂| 亚洲成色77777| 曰老女人黄片| 五月伊人婷婷丁香| 在线观看人妻少妇| 激情五月婷婷亚洲| 欧美激情极品国产一区二区三区 | 亚洲av国产av综合av卡| 久久久久久人人人人人| 一本色道久久久久久精品综合| 国产成人午夜福利电影在线观看| 国产精品久久久久久av不卡| 美女大奶头黄色视频| 欧美3d第一页| 一本—道久久a久久精品蜜桃钙片| 少妇的逼水好多| 国产成人免费无遮挡视频| 中文字幕人妻丝袜制服| 久久久精品区二区三区| 女人精品久久久久毛片| 赤兔流量卡办理| 久久国产亚洲av麻豆专区| 久久99一区二区三区| 国产无遮挡羞羞视频在线观看| av免费在线看不卡| 中国美白少妇内射xxxbb| 少妇猛男粗大的猛烈进出视频| 亚洲一码二码三码区别大吗| 99热这里只有是精品在线观看| 精品一区二区免费观看| 午夜免费男女啪啪视频观看| 国产1区2区3区精品| 永久网站在线| 欧美丝袜亚洲另类| 亚洲av.av天堂| 蜜桃在线观看..| videossex国产| 精品人妻偷拍中文字幕| 久久人人97超碰香蕉20202| 中文精品一卡2卡3卡4更新| 国产伦理片在线播放av一区| 卡戴珊不雅视频在线播放| 满18在线观看网站| 婷婷成人精品国产| 亚洲av男天堂| 免费看不卡的av| 99九九在线精品视频| 婷婷色麻豆天堂久久| 亚洲美女搞黄在线观看| 国产精品三级大全| 人人妻人人澡人人爽人人夜夜| 国产黄色视频一区二区在线观看| 精品国产一区二区久久| 国产又色又爽无遮挡免| 好男人视频免费观看在线| 亚洲精品日韩在线中文字幕| 精品亚洲乱码少妇综合久久| 免费久久久久久久精品成人欧美视频 | 成人毛片60女人毛片免费| 免费看av在线观看网站| 777米奇影视久久| 99re6热这里在线精品视频| 最近手机中文字幕大全| 欧美日韩视频精品一区| 久久久久久久亚洲中文字幕| 看免费av毛片| 人成视频在线观看免费观看| 黄色 视频免费看| av免费在线看不卡| 在线精品无人区一区二区三| 中文乱码字字幕精品一区二区三区| 观看av在线不卡| 久久久久久久亚洲中文字幕| av国产精品久久久久影院| 成人毛片a级毛片在线播放| 国产爽快片一区二区三区| 国国产精品蜜臀av免费| 成人亚洲欧美一区二区av| 99热这里只有是精品在线观看| 一本—道久久a久久精品蜜桃钙片| 深夜精品福利| 精品一区二区三区视频在线| 亚洲欧美成人综合另类久久久| 成年动漫av网址| 久久 成人 亚洲| 亚洲激情五月婷婷啪啪| 国产在线一区二区三区精| 精品一区在线观看国产| 熟女电影av网| 午夜激情av网站| 少妇人妻 视频| 欧美变态另类bdsm刘玥| 精品国产国语对白av| 国产精品成人在线| 欧美精品高潮呻吟av久久| av国产精品久久久久影院| 日韩人妻精品一区2区三区| 高清毛片免费看| 成人漫画全彩无遮挡| 国产精品久久久久久久电影| 久热久热在线精品观看| 欧美国产精品一级二级三级| 99精国产麻豆久久婷婷| 黄网站色视频无遮挡免费观看| 大码成人一级视频| 性高湖久久久久久久久免费观看| 国产黄色视频一区二区在线观看| 亚洲av.av天堂| 久久久久精品人妻al黑| 狂野欧美激情性bbbbbb| 亚洲国产精品专区欧美| a级片在线免费高清观看视频| 最后的刺客免费高清国语| 日本色播在线视频| av视频免费观看在线观看| 日韩免费高清中文字幕av| 国产 一区精品| 精品人妻在线不人妻| 国语对白做爰xxxⅹ性视频网站| 成人手机av| 亚洲精品久久成人aⅴ小说| 精品少妇内射三级| 亚洲美女黄色视频免费看| 夜夜爽夜夜爽视频| 中文字幕av电影在线播放| 欧美亚洲 丝袜 人妻 在线| 国产色婷婷99| 久久午夜福利片| 街头女战士在线观看网站| 秋霞伦理黄片| 春色校园在线视频观看| 亚洲天堂av无毛| 国产成人午夜福利电影在线观看|