• <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)特性分析
    耐火鋼圓鋼管混凝土柱耐火極限和承載力
    黄色女人牲交| 麻豆久久精品国产亚洲av| 午夜福利高清视频| 91麻豆精品激情在线观看国产| or卡值多少钱| 97超级碰碰碰精品色视频在线观看| 国产精品98久久久久久宅男小说| 一进一出抽搐gif免费好疼| 欧美黑人精品巨大| 免费女性裸体啪啪无遮挡网站| 成年人黄色毛片网站| 午夜成年电影在线免费观看| 一二三四在线观看免费中文在| 波多野结衣av一区二区av| 国产亚洲欧美在线一区二区| 国产精品久久久av美女十八| 可以在线观看毛片的网站| 我的亚洲天堂| 黄片播放在线免费| √禁漫天堂资源中文www| www.精华液| 夜夜躁狠狠躁天天躁| 一区二区三区国产精品乱码| 欧美乱色亚洲激情| 亚洲精品一卡2卡三卡4卡5卡| 亚洲中文字幕一区二区三区有码在线看 | 日本熟妇午夜| 日本成人三级电影网站| 老熟妇仑乱视频hdxx| 亚洲国产看品久久| 天堂动漫精品| 99精品欧美一区二区三区四区| a级毛片在线看网站| 亚洲欧美精品综合久久99| 久久精品国产综合久久久| 18禁观看日本| 黄色视频,在线免费观看| 午夜老司机福利片| 丁香六月欧美| 午夜免费成人在线视频| 国产精品九九99| 久9热在线精品视频| 久久久久久久久免费视频了| 动漫黄色视频在线观看| 日韩三级视频一区二区三区| 精品一区二区三区四区五区乱码| 亚洲一区二区三区不卡视频| 大型黄色视频在线免费观看| 黑丝袜美女国产一区| 三级毛片av免费| 99国产精品99久久久久| 好男人电影高清在线观看| 日韩三级视频一区二区三区| 一级作爱视频免费观看| 丁香六月欧美| 国产精品久久久av美女十八| 日本 欧美在线| 美女扒开内裤让男人捅视频| 日本成人三级电影网站| 亚洲av成人av| 99在线视频只有这里精品首页| 精品一区二区三区av网在线观看| 男人舔女人的私密视频| 国内毛片毛片毛片毛片毛片| av片东京热男人的天堂| 国产精品乱码一区二三区的特点| 无遮挡黄片免费观看| 亚洲九九香蕉| 亚洲精品国产一区二区精华液| 成人一区二区视频在线观看| xxx96com| 在线观看一区二区三区| 亚洲成人国产一区在线观看| 欧美黄色片欧美黄色片| 亚洲人成77777在线视频| 午夜免费激情av| 别揉我奶头~嗯~啊~动态视频| 最近在线观看免费完整版| 欧美乱妇无乱码| 中文字幕精品亚洲无线码一区 | 精品国内亚洲2022精品成人| 欧美绝顶高潮抽搐喷水| 国产一级毛片七仙女欲春2 | 黄色片一级片一级黄色片| 国产精品久久电影中文字幕| 成人欧美大片| 欧美不卡视频在线免费观看 | 亚洲 欧美 日韩 在线 免费| 亚洲av成人av| 亚洲一区高清亚洲精品| 搡老熟女国产l中国老女人| 免费看a级黄色片| 欧美激情高清一区二区三区| 国产精品 欧美亚洲| 一边摸一边做爽爽视频免费| 免费高清视频大片| av片东京热男人的天堂| 亚洲精品一卡2卡三卡4卡5卡| 香蕉av资源在线| 在线观看免费日韩欧美大片| 久久国产亚洲av麻豆专区| 亚洲精品在线观看二区| 亚洲aⅴ乱码一区二区在线播放 | 国产精品综合久久久久久久免费| 一个人免费在线观看的高清视频| 精品日产1卡2卡| 欧美国产精品va在线观看不卡| 精品国产一区二区三区四区第35| 男人的好看免费观看在线视频 | 久久精品人妻少妇| 久久中文字幕一级| 日本三级黄在线观看| 夜夜看夜夜爽夜夜摸| 99在线人妻在线中文字幕| 一进一出抽搐动态| 日韩国内少妇激情av| 日韩 欧美 亚洲 中文字幕| 夜夜爽天天搞| 欧美最黄视频在线播放免费| 亚洲精品一区av在线观看| 制服丝袜大香蕉在线| 人妻久久中文字幕网| 欧美乱码精品一区二区三区| 欧美日本亚洲视频在线播放| 久久精品国产亚洲av高清一级| 一卡2卡三卡四卡精品乱码亚洲| 日韩精品中文字幕看吧| 1024香蕉在线观看| 热99re8久久精品国产| 中文资源天堂在线| 午夜福利一区二区在线看| 1024视频免费在线观看| 日日摸夜夜添夜夜添小说| 午夜免费成人在线视频| a在线观看视频网站| xxx96com| 午夜影院日韩av| 成人午夜高清在线视频 | 国产一区二区在线av高清观看| aaaaa片日本免费| 国产精品乱码一区二三区的特点| 免费高清视频大片| 老司机福利观看| 亚洲人成网站在线播放欧美日韩| 黑丝袜美女国产一区| 亚洲男人的天堂狠狠| 亚洲人成网站在线播放欧美日韩| 国产一区二区激情短视频| 欧美激情 高清一区二区三区| 亚洲欧美激情综合另类| 久久亚洲真实| 日本免费a在线| 亚洲av第一区精品v没综合| 免费看十八禁软件| 亚洲成国产人片在线观看| 欧美最黄视频在线播放免费| 日本成人三级电影网站| 国产不卡一卡二| 嫩草影院精品99| 99久久精品国产亚洲精品| 亚洲国产欧美网| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品国产一区二区精华液| 亚洲一卡2卡3卡4卡5卡精品中文| 久久草成人影院| 国产成人一区二区三区免费视频网站| 国产精品亚洲av一区麻豆| 成人国产综合亚洲| 国产一区在线观看成人免费| videosex国产| 三级毛片av免费| 国产激情欧美一区二区| 18禁观看日本| 波多野结衣高清作品| 亚洲一卡2卡3卡4卡5卡精品中文| 久久 成人 亚洲| 国产精品98久久久久久宅男小说| 一级a爱片免费观看的视频| 国产精品av久久久久免费| 人人妻人人澡人人看| 天堂影院成人在线观看| 又黄又爽又免费观看的视频| 久久婷婷人人爽人人干人人爱| 色精品久久人妻99蜜桃| 欧美黄色淫秽网站| 亚洲国产欧美网| 老熟妇仑乱视频hdxx| 国产黄a三级三级三级人| 国产精品久久视频播放| 日韩大尺度精品在线看网址| 18禁观看日本| 无遮挡黄片免费观看| 宅男免费午夜| 亚洲欧美日韩无卡精品| 一进一出好大好爽视频| 一进一出好大好爽视频| cao死你这个sao货| 欧美日韩亚洲国产一区二区在线观看| 亚洲中文日韩欧美视频| 免费在线观看视频国产中文字幕亚洲| 手机成人av网站| 男女下面进入的视频免费午夜 | 最好的美女福利视频网| 久久99热这里只有精品18| 一级片免费观看大全| 午夜免费激情av| 国产精华一区二区三区| 日本在线视频免费播放| 国产亚洲精品综合一区在线观看 | 亚洲人成网站在线播放欧美日韩| 国产亚洲av高清不卡| 精品不卡国产一区二区三区| 在线视频色国产色| 午夜a级毛片| 日韩欧美免费精品| 无人区码免费观看不卡| 亚洲欧美日韩高清在线视频| 一区二区三区激情视频| 岛国视频午夜一区免费看| 青草久久国产| 成熟少妇高潮喷水视频| 亚洲国产日韩欧美精品在线观看 | 国产三级在线视频| 香蕉av资源在线| 日本免费a在线| 久久久精品国产亚洲av高清涩受| 人妻丰满熟妇av一区二区三区| 成年版毛片免费区| 丁香欧美五月| 日本 欧美在线| 成人18禁高潮啪啪吃奶动态图| 婷婷精品国产亚洲av| 在线视频色国产色| 亚洲人成77777在线视频| 中文字幕人妻熟女乱码| 免费观看精品视频网站| 男人舔奶头视频| 日本a在线网址| 少妇的丰满在线观看| 亚洲一区二区三区不卡视频| av欧美777| 麻豆一二三区av精品| 精品福利观看| 国产又爽黄色视频| 精品久久久久久久毛片微露脸| 久久久国产成人免费| 久久热在线av| 女同久久另类99精品国产91| www.999成人在线观看| 首页视频小说图片口味搜索| 国产激情欧美一区二区| 久久久久久人人人人人| 欧美激情高清一区二区三区| 国产激情久久老熟女| 国产成+人综合+亚洲专区| 久久久久亚洲av毛片大全| 国产黄色小视频在线观看| 精品久久久久久成人av| 国产黄片美女视频| 男女做爰动态图高潮gif福利片| 亚洲国产欧美一区二区综合| 久久精品国产99精品国产亚洲性色| 国产成+人综合+亚洲专区| 日韩欧美国产一区二区入口| 88av欧美| 久久久精品国产亚洲av高清涩受| 亚洲av五月六月丁香网| 欧美激情 高清一区二区三区| 一级作爱视频免费观看| 在线观看舔阴道视频| 亚洲国产毛片av蜜桃av| 极品教师在线免费播放| 亚洲男人的天堂狠狠| 国产一区二区三区视频了| 免费高清在线观看日韩| 成年女人毛片免费观看观看9| 久热爱精品视频在线9| 亚洲精品美女久久av网站| 国产亚洲欧美98| 亚洲精品一卡2卡三卡4卡5卡| 久久精品影院6| 桃色一区二区三区在线观看| 精品少妇一区二区三区视频日本电影| 久久久水蜜桃国产精品网| 国产av又大| 不卡av一区二区三区| 亚洲久久久国产精品| 99久久精品国产亚洲精品| 黄色成人免费大全| 欧美日韩中文字幕国产精品一区二区三区| 首页视频小说图片口味搜索| 精品不卡国产一区二区三区| 久久久久久国产a免费观看| 99久久精品国产亚洲精品| 啦啦啦韩国在线观看视频| 欧美激情久久久久久爽电影| 国产av在哪里看| 精品国产乱码久久久久久男人| 中出人妻视频一区二区| 看免费av毛片| 成人国产一区最新在线观看| 久久中文字幕人妻熟女| 欧洲精品卡2卡3卡4卡5卡区| 757午夜福利合集在线观看| 午夜福利一区二区在线看| 亚洲成人久久爱视频| 国产免费男女视频| 久久亚洲精品不卡| 在线观看免费日韩欧美大片| 91九色精品人成在线观看| 丰满的人妻完整版| 久久精品91无色码中文字幕| 中文字幕av电影在线播放| 欧美黑人精品巨大| 久久九九热精品免费| 91老司机精品| 国产精品久久久久久亚洲av鲁大| 久久亚洲真实| 午夜久久久在线观看| 欧美成人午夜精品| 久久中文看片网| 国内精品久久久久久久电影| 国产精品98久久久久久宅男小说| 亚洲久久久国产精品| 国产精品亚洲av一区麻豆| 色老头精品视频在线观看| 哪里可以看免费的av片| 国产v大片淫在线免费观看| 国产精品国产高清国产av| 欧美黑人欧美精品刺激| 欧美日韩一级在线毛片| 国产精品一区二区精品视频观看| 丰满的人妻完整版| 18禁黄网站禁片午夜丰满| 国产精品亚洲一级av第二区| 女生性感内裤真人,穿戴方法视频| 99久久精品国产亚洲精品| 激情在线观看视频在线高清| 国产精品香港三级国产av潘金莲| 亚洲五月婷婷丁香| 亚洲专区字幕在线| 免费电影在线观看免费观看| 可以在线观看的亚洲视频| 香蕉久久夜色| 久久精品夜夜夜夜夜久久蜜豆 | 国产欧美日韩一区二区三| 亚洲成人久久爱视频| 成年版毛片免费区| 久久中文字幕人妻熟女| 99久久综合精品五月天人人| 国产精品亚洲一级av第二区| 免费看美女性在线毛片视频| 丰满人妻熟妇乱又伦精品不卡| 久久久久久九九精品二区国产 | 久久人妻福利社区极品人妻图片| 亚洲欧洲精品一区二区精品久久久| 男人舔奶头视频| 嫁个100分男人电影在线观看| 精品国产国语对白av| 欧美黄色片欧美黄色片| 在线av久久热| 亚洲中文字幕日韩| 日本免费a在线| 久久青草综合色| 欧美午夜高清在线| 国产精品亚洲av一区麻豆| 欧美亚洲日本最大视频资源| 亚洲精品国产精品久久久不卡| www国产在线视频色| 亚洲欧美精品综合一区二区三区| 亚洲av日韩精品久久久久久密| 黄色丝袜av网址大全| 亚洲精品一区av在线观看| cao死你这个sao货| 免费在线观看日本一区| 神马国产精品三级电影在线观看 | 69av精品久久久久久| 九色国产91popny在线| 亚洲成人精品中文字幕电影| 18禁裸乳无遮挡免费网站照片 | 亚洲自拍偷在线| 亚洲aⅴ乱码一区二区在线播放 | 99久久综合精品五月天人人| 国产午夜福利久久久久久| 99久久无色码亚洲精品果冻| 脱女人内裤的视频| 亚洲av成人不卡在线观看播放网| 亚洲精华国产精华精| 热re99久久国产66热| 国产私拍福利视频在线观看| 欧美一级毛片孕妇| 亚洲成人久久性| 91国产中文字幕| 美国免费a级毛片| 亚洲成人国产一区在线观看| 很黄的视频免费| 国产爱豆传媒在线观看 | 真人做人爱边吃奶动态| 一个人免费在线观看的高清视频| 欧美乱色亚洲激情| 一个人观看的视频www高清免费观看 | a在线观看视频网站| 久久精品91蜜桃| 亚洲成av人片免费观看| 成人国产一区最新在线观看| 女人被狂操c到高潮| 欧美在线黄色| 中文在线观看免费www的网站 | 最近在线观看免费完整版| 日韩免费av在线播放| 精品一区二区三区四区五区乱码| 村上凉子中文字幕在线| 色综合欧美亚洲国产小说| 精品国产一区二区三区四区第35| 精品久久久久久成人av| 亚洲,欧美精品.| 久久青草综合色| 亚洲av日韩精品久久久久久密| 久久精品国产综合久久久| 搡老妇女老女人老熟妇| 国产精品一区二区精品视频观看| 久久九九热精品免费| 久久久久久人人人人人| 亚洲一区高清亚洲精品| 亚洲色图av天堂| 色播亚洲综合网| 国产视频一区二区在线看| tocl精华| 亚洲av成人不卡在线观看播放网| 久久草成人影院| 亚洲精品久久成人aⅴ小说| 午夜久久久久精精品| 欧美黄色片欧美黄色片| 一区二区三区国产精品乱码| 精品国产乱码久久久久久男人| 色老头精品视频在线观看| 欧美成人一区二区免费高清观看 | 老熟妇乱子伦视频在线观看| www.www免费av| 黄色 视频免费看| 很黄的视频免费| 18禁黄网站禁片免费观看直播| 日本在线视频免费播放| 变态另类丝袜制服| 又黄又粗又硬又大视频| 欧美成人一区二区免费高清观看 | 女同久久另类99精品国产91| 欧美日韩一级在线毛片| 亚洲国产精品999在线| 亚洲精品国产一区二区精华液| 亚洲成人国产一区在线观看| 男男h啪啪无遮挡| 免费看美女性在线毛片视频| 国产精品1区2区在线观看.| 嫁个100分男人电影在线观看| 欧美一区二区精品小视频在线| 夜夜夜夜夜久久久久| av在线播放免费不卡| 国产av一区在线观看免费| 深夜精品福利| 韩国av一区二区三区四区| 久久精品国产亚洲av高清一级| 女同久久另类99精品国产91| 18禁黄网站禁片午夜丰满| 亚洲国产毛片av蜜桃av| 久久午夜综合久久蜜桃| 女生性感内裤真人,穿戴方法视频| 搡老熟女国产l中国老女人| 亚洲国产精品成人综合色| 精品国产超薄肉色丝袜足j| 欧美日韩福利视频一区二区| 亚洲第一av免费看| 欧美大码av| 国产午夜精品久久久久久| 国产免费男女视频| 午夜免费激情av| 韩国精品一区二区三区| 欧美成人免费av一区二区三区| 91麻豆精品激情在线观看国产| 丁香六月欧美| 欧美zozozo另类| 妹子高潮喷水视频| 成年免费大片在线观看| 欧美午夜高清在线| 国产成人精品久久二区二区免费| 最近最新中文字幕大全免费视频| 亚洲熟女毛片儿| 国产极品粉嫩免费观看在线| 一a级毛片在线观看| 人人澡人人妻人| 听说在线观看完整版免费高清| 欧美精品亚洲一区二区| 免费在线观看黄色视频的| 亚洲自拍偷在线| 中文字幕最新亚洲高清| 精品久久久久久久毛片微露脸| 老熟妇乱子伦视频在线观看| 国产熟女午夜一区二区三区| 久久国产乱子伦精品免费另类| 在线播放国产精品三级| 亚洲精品国产精品久久久不卡| 欧美国产精品va在线观看不卡| 国产精华一区二区三区| 91成人精品电影| 天堂影院成人在线观看| 久久久精品欧美日韩精品| av免费在线观看网站| 免费在线观看黄色视频的| 国产真人三级小视频在线观看| 国产又黄又爽又无遮挡在线| 俄罗斯特黄特色一大片| 黄网站色视频无遮挡免费观看| 国内精品久久久久久久电影| 国产免费av片在线观看野外av| 亚洲av中文字字幕乱码综合 | 久久精品91蜜桃| 桃色一区二区三区在线观看| 久久中文看片网| 搡老熟女国产l中国老女人| 在线观看午夜福利视频| 精品人妻1区二区| 久久久久精品国产欧美久久久| 亚洲精品国产精品久久久不卡| 岛国在线观看网站| 中文字幕人妻熟女乱码| 男人舔女人下体高潮全视频| 中文字幕av电影在线播放| 国产成人精品久久二区二区91| 老汉色∧v一级毛片| 亚洲一区高清亚洲精品| 亚洲av美国av| 不卡av一区二区三区| 欧美成狂野欧美在线观看| 亚洲欧美日韩无卡精品| 欧美一级a爱片免费观看看 | 精品午夜福利视频在线观看一区| 亚洲精品色激情综合| 搡老熟女国产l中国老女人| 一本久久中文字幕| 男人舔奶头视频| 亚洲欧洲精品一区二区精品久久久| or卡值多少钱| 亚洲成人免费电影在线观看| 欧美丝袜亚洲另类 | 国产在线观看jvid| 啦啦啦观看免费观看视频高清| 狠狠狠狠99中文字幕| 999久久久精品免费观看国产| 狂野欧美激情性xxxx| 亚洲精品色激情综合| 亚洲va日本ⅴa欧美va伊人久久| 好男人在线观看高清免费视频 | 亚洲欧美日韩高清在线视频| 最近最新中文字幕大全免费视频| 悠悠久久av| 99热这里只有精品一区 | 国产亚洲精品第一综合不卡| 日本一本二区三区精品| 久久精品国产亚洲av高清一级| 日韩欧美 国产精品| 久久精品国产99精品国产亚洲性色| www.999成人在线观看| 又黄又爽又免费观看的视频| 日本一区二区免费在线视频| 在线永久观看黄色视频| av电影中文网址| 真人一进一出gif抽搐免费| 午夜福利欧美成人| 真人做人爱边吃奶动态| 丝袜人妻中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 50天的宝宝边吃奶边哭怎么回事| 69av精品久久久久久| 韩国精品一区二区三区| 国产又黄又爽又无遮挡在线| 51午夜福利影视在线观看| 欧美亚洲日本最大视频资源| 亚洲精品国产一区二区精华液| 91九色精品人成在线观看| 国产色视频综合| 国产在线观看jvid| 黄色毛片三级朝国网站| 久久久久国内视频| 美女大奶头视频| 免费一级毛片在线播放高清视频| 国产精品 国内视频| 此物有八面人人有两片| 黑人欧美特级aaaaaa片| 国产成人精品无人区| 999精品在线视频| 法律面前人人平等表现在哪些方面| 久久久精品欧美日韩精品| 国产99白浆流出| 国产成人啪精品午夜网站| 一个人免费在线观看的高清视频| 久久精品91蜜桃| 亚洲avbb在线观看| 成在线人永久免费视频| 老汉色∧v一级毛片| 国产精品 欧美亚洲| 日韩欧美免费精品| 一进一出抽搐gif免费好疼| 国产高清videossex| 亚洲九九香蕉| 在线观看66精品国产| 好男人在线观看高清免费视频 | 国产精品98久久久久久宅男小说| 欧美成人一区二区免费高清观看 | 最新美女视频免费是黄的| 在线天堂中文资源库| 久久精品亚洲精品国产色婷小说| 国产精品久久久久久人妻精品电影| 久久香蕉精品热| 美女 人体艺术 gogo| 国产精品香港三级国产av潘金莲|