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

    變剛度支承下電磁軸承-轉子系統(tǒng)振動特性研究

    2022-06-08 09:33:18魏子航宋春生李俊李民輝
    軸承 2022年3期
    關鍵詞:剛性振型固有頻率

    魏子航,宋春生, 2,李俊,李民輝

    (1.武漢理工大學 機電工程學院,武漢 430070; 2.湖北省磁懸浮工程技術研究中心,武漢 430070)

    電磁軸承具有無接觸,無摩擦,發(fā)熱少,無需潤滑系統(tǒng),長壽命以及可在極端環(huán)境下工作等優(yōu)點。此外,電磁軸承具有極佳的主動振動控制功能,可以根據(jù)實際工況通過控制策略在一定范圍內(nèi)調(diào)節(jié)自身的支承參數(shù),從而抑制轉子的不平衡振動。文獻[1]對電磁輔助支承的變剛度減振進行研究,在一階臨界轉速附近取得了較好的減振效果;文獻[2]在不同轉速區(qū)間采用不同控制參數(shù),能夠有效減小轉子的不平衡振動;文獻[3]研究了基于控制參數(shù)切換的磁懸浮旋轉機械隔振技術,并給出了控制參數(shù)的設計思路及確定參數(shù)切換轉速的方法:以上文獻大多聚焦于變參數(shù)控制的設計及實現(xiàn),未分析電磁軸承等效剛度變化對轉子系統(tǒng)振動特性產(chǎn)生的影響。

    本文較為系統(tǒng)地對變剛度電磁軸承-轉子系統(tǒng)的動力學及振動傳遞特性進行研究:建立電磁軸承支承的剛性轉子動力學模型并驗證電磁軸承等效剛度與控制參數(shù)之間的關系,分析等效剛度的影響因素;將電磁軸承-剛性轉子、柔性轉子數(shù)學模型與有限元仿真軟件結合,研究等效剛度對轉子系統(tǒng)模態(tài)頻率、振型的影響;分析電磁軸承支承下的轉子振動傳遞特性。

    1 電磁軸承電磁力及廣義剛度

    本試驗裝置的電磁軸承共有8個磁極,相鄰2個磁極通入相同的控制電流并采用NNSSNN順序進行繞線,以45°偏置安裝,其結構參數(shù)見表1。

    表1 電磁軸承結構參數(shù)Tab.1 Structural parameters of AMB

    電磁軸承的控制系統(tǒng)包括控制器、傳感器、功率放大器。考慮一般情況,設電磁軸承的控制規(guī)律為

    G(s)=Gs(s)Ga(s)Gc(s),

    (1)

    式中:Gs(s),Ga(s),Gc(s)分別為傳感器、功率放大器、控制器的傳遞函數(shù),傳遞函數(shù)G(s)中的復變量s=σ+jω,令σ=0可得到系統(tǒng)的頻率特性G(jω)。

    電磁軸承等效剛度、阻尼的一般表達式為[4]

    (2)

    式中:ω為轉子角速度;P(ω),Q(ω)分別為G(jω)的實部和虛部。

    對于PID控制,控制器的傳遞函數(shù)為

    (3)

    式中:cP,cI,cD分別為比例系數(shù)、積分系數(shù)、微分系數(shù);TD為微分時間常數(shù)。

    功率放大器與傳感器的傳遞函數(shù)以一階慣性環(huán)節(jié)表示,可得PID控制下等效剛度、阻尼表達式為

    (4)

    式中:Aa,As分別為功率放大器、傳感器的放大倍數(shù);Ta,Ts分別為功率放大器、傳感器的時間常數(shù)。

    由于系統(tǒng)固有頻率受電磁軸承等效剛度的影響,因此分別分析PID控制參數(shù)及轉子轉頻對等效剛度的影響,如圖1—圖3所示。

    圖1 等效剛度隨比例系數(shù)cP的變化Fig.1 Variation of equivalent stiffness with cP

    由圖1—圖3可知:在給定的轉速區(qū)間,比例系數(shù)cP對等效剛度的影響最為明顯,兩者之間存在近似線性的關系;積分系數(shù)cI對等效剛度的影響非常??;等效剛度受微分系數(shù)cD與轉子轉速共同作用的影響,當微分系數(shù)較小時,等效剛度能夠在較大轉速范圍內(nèi)保持不變。

    圖2 等效剛度隨積分系數(shù)cI的變化Fig.2 Variation of equivalent stiffness with cI

    圖3 等效剛度隨微分系數(shù)cD的變化Fig.3 Variation of equivalent stiffness with cD

    綜上可知,在較大的轉速范圍內(nèi),比例系數(shù)cP為等效剛度的主要影響因素且兩者體現(xiàn)出近似線性關系,而轉速和微分系數(shù)對等效剛度的影響相對較小,因此可將等效剛度簡化為與比例系數(shù)cP有關的線性化表達式以便于分析。

    2 電磁軸承-轉子系統(tǒng)動力學分析

    2.1 電磁軸承支承剛性轉子動力學模型

    磁懸浮轉子受力如圖4所示,θy,θx分別為轉子繞y,x軸的轉角。

    圖4 磁懸浮轉子受力Fig.4 Force of magnetic suspension rotor

    由于積分參數(shù)對等效剛度和固有頻率的影響不大,若采用PD控制,并取各通道控制參數(shù)相同,只考慮穩(wěn)態(tài)響應下的情況,依據(jù)轉子動力學相關理論[5]可得

    (5)

    q=(θy,x,θx,y)T,

    KsS=BKsBT,

    式中:q為轉子的位移矢量;P為比例系數(shù)矩陣;J,Jz分別為轉子的直徑轉動慣量、極轉動慣量;KsS為軸承負剛度矩陣;M,G,D分別為轉子的質(zhì)量、陀螺、阻尼矩陣;B,C分別為電磁軸承輸入矩陣和傳感器輸出矩陣。

    令Kc=BKiPC為PD控制下的剛度矩陣,Dc=BKiDC為PD控制下的阻尼矩陣,其中Ki為2套電磁軸承在2個方向上的電流剛度矩陣。并將該微分方程轉換為狀態(tài)空間描述,可得狀態(tài)空間矩陣為

    (6)

    令行列式det(λI-A)=0可求得狀態(tài)空間矩陣A的特征值λ,進而得到電磁軸承-剛性轉子系統(tǒng)的固有頻率,得到其坎貝爾圖如圖5所示:隨著轉速的增加,轉子特征頻率線開始出現(xiàn)分叉現(xiàn)象,正向渦動頻率隨轉速的增加而增加,反向渦動頻率隨轉速的增加而減?。晦D子的平動、錐動頻率隨轉速的增加沒有發(fā)生明顯的分叉,說明轉子本身是細長的,陀螺效應對臨界轉速的影響不明顯。

    圖5 轉子系統(tǒng)的坎貝爾圖Fig.5 Campbell diagram of rotor system

    2.2 等效剛度與等效阻尼的有限元仿真驗證

    參考電磁軸承的常見等效剛度范圍,在確保系統(tǒng)穩(wěn)定的前提下,選取一系列方案的控制參數(shù),見表2。

    表2 電磁軸承控制參數(shù)Tab.2 Control parameters of AMB

    在不考慮時滯的條件下,根據(jù)(7)式計算得到各方案下等效剛度和等效阻尼,見表3,阻尼比ζ可以通過(8)式計算,即

    表3 電磁軸承等效剛度和等效阻尼Tab.3 Equivalent stiffness and damping of AMB

    (7)

    (8)

    式中:m為轉子質(zhì)量。

    為與電磁軸承-剛性轉子數(shù)學模型保持一致,在Workbench中將轉子設置為剛體,2組徑向電磁軸承簡化為對地彈簧阻尼單元,連接到轉子疊壓硅鋼片位置,約束其徑向4個自由度并做模態(tài)分析,平動模態(tài)頻率計算結果見表4。

    表4 剛性轉子平動模態(tài)頻率Tab.4 Parallel modal frequency of rigid rotor

    由表4可知,剛性轉子的平動模態(tài)頻率有限元仿真與電磁軸承-剛性轉子數(shù)學模型計算的結果基本一致,此時電磁軸承支承下轉子振型為明顯的彈性支承-剛性轉子形式,如圖6所示。說明電磁軸承的等效剛度和等效阻尼能夠很好地反映其支承特性,觀察方案4,7,8的結果可知,隨著等效阻尼的增大,固有頻率降低但影響不大。

    圖6 等效剛度50.00×106 N/m時平動模態(tài)振型圖Fig.6 Parallel mode shape with equivalent stiffness of50.00×106 N/m

    2.3 電磁軸承支承柔性轉子動力學模型

    由于需要采集諧響應分析中電磁軸承處的加速度信號以便觀察電磁軸承的支承位置,將電磁軸承定子簡化為薄壁圓筒,與轉子之間建立彈簧阻尼約束,同時認為電磁軸承基座部分不發(fā)生變形,約束電磁軸承外表面的位移,修改后的模型如圖7所示。

    圖7 簡化電磁軸承支承轉子模型Fig.7 Simplified model of rotor supported by AMB

    根據(jù)控制方案1,2,3,4,5,6進行約束模態(tài)分析,并與2.2節(jié)電磁軸承-剛性轉子模型得到的固有頻率進行比較,結果見表5。

    表5 轉子系統(tǒng)平動模態(tài)頻率Tab.5 Parallel modal frequency of rotor system Hz

    由表5可知,當?shù)刃偠葹?.15×106,0.30×106,0.50×106,1.00×106N/m時,兩者結果較為一致,而當?shù)刃偠葹?0.00×106,50.00×106N/m時,固有頻率出現(xiàn)了明顯差異,此時等效剛度較大,轉子的一階平動模態(tài)出現(xiàn)了剛性支承-柔性轉子的特點,柔性轉子模型更能反映此時的轉子振動特性。

    利用有限元理論建立柔性轉子的數(shù)學模型,將轉子分為56個單元,依據(jù)每一單元的尺寸與材料屬性,利用MATLAB計算得到每一單元的質(zhì)量、剛度、陀螺矩陣,將所有單元矩陣組合得到系統(tǒng)質(zhì)量、剛度、陀螺矩陣[6],將電磁軸承剛度疊加到剛度矩陣的對應節(jié)點上,得到運動方程,即

    (9)

    y=Cu,

    式中:Ω為轉子自轉角速度;B1為外部擾動力分布矩陣;f為外部擾動力;u為各節(jié)點在x與y方向組成的位移與轉動矢量。

    對轉子模型進行模態(tài)截斷降階,在模態(tài)矩陣Φ中保留前4階模態(tài)并進行模態(tài)變換,令u=Φη,代入(9)式可得

    (10)

    y=CΦη,

    式中:Φ為模態(tài)矩陣;η為模態(tài)坐標。

    令Mre=ΦTMΦ,Kre=ΦTKΦ,Cre=CΦ,Gre=ΦT(D+ΩG)Φ可得

    (11)

    y=Creη。

    金屬轉子取阻尼比ζ=0.3%,并忽略其陀螺效應,經(jīng)質(zhì)量量綱一化處理后可得轉子狀態(tài)空間運動方程為

    (12)

    (13)

    轉子降階后為

    (14)

    Cr=[Cre0]。

    2.4 電磁軸承-柔性轉子系統(tǒng)模態(tài)計算

    通過計算電磁軸承-柔性轉子系統(tǒng)的固有頻率驗證理論模型的準確性,利用MATLAB中的eig函數(shù)求解系統(tǒng)矩陣得到系統(tǒng)固有頻率及振型,用Workbench模態(tài)分析得到轉子在電磁軸承支承下的前3階固有頻率及振型,繪制固有頻率隨等效剛度的變化如圖8所示。

    圖8 固有頻率隨等效剛度的變化Fig.8 Variation of natural frequency with equivalent stiffness

    由圖8可知,利用MATLAB與Workbench計算的固有頻率基本一致,驗證了理論與仿真模型的準確性。在不同的等效剛度區(qū)間,剛度變化對各階固有頻率的影響程度不同。

    文獻[7]指出,簡支梁在彈性支承下的一階彈性模態(tài)可由其在自由狀態(tài)下的一階剛性模態(tài)與一階彎曲模態(tài)線性組合得到。電磁軸承等效剛度對轉子系統(tǒng)動態(tài)特性的影響分為3個階段:當?shù)刃偠容^低,為0.50×106N/m時,轉子相對更接近自由邊界,因而一階剛性模態(tài)貢獻更大,對應振型為如圖9所示的彈性支承-剛性轉子形式,等效剛度僅對前2階剛性固有頻率有明顯的影響,而對轉子一階彎曲固有頻率影響微?。划?shù)刃偠葹?.00×106N/m時,等效剛度對前3階固有頻率均有比較明顯的影響,彎曲模態(tài)對一階振型的影響逐漸增加,一階平動振型出現(xiàn)明顯彈性模態(tài),且節(jié)點最先出現(xiàn)在轉子兩端部,如圖10所示;當?shù)刃偠容^高,為10.00×106N/m時,一階平動固有頻率變化不再明顯,對應振型的節(jié)點位置隨著等效剛度的增加向電磁軸承支承位置不斷移動,如圖11所示。由于控制穩(wěn)定性以及結構尺寸限制的因素,實際電磁軸承的等效剛度在(0.10~10.00)×106N/m之間[8],理論上可以通過調(diào)節(jié)控制參數(shù),實現(xiàn)一階平動模態(tài)或轉子一階彎曲模態(tài)的固有頻率偏移,減小跨臨界引起的振動。對剛性、柔性轉子數(shù)學模型與有限元仿真計算得到的一階平動模態(tài)頻率進行對比,見表6。

    圖9 等效剛度0.50×106 N/m時模態(tài)振型圖Fig.9 Mode shape with equivalent stiffness of0.50×106 N/m

    圖10 等效剛度5.00×106 N/m時模態(tài)振型圖Fig.10 Mode shape with equivalent stiffness of5.00×106 N/m

    圖11 等效剛度10.00×106 N/m時模態(tài)振型圖Fig.11 Mode shape with equivalent stiffness of10.00×106 N/m

    表6 轉子系統(tǒng)平動模態(tài)頻率理論與仿真結果對比Tab.6 Comparison of theory and simulation results of parallel modal frequency of rotor system Hz

    由表6可知:剛性轉子與柔性轉子數(shù)學模型在電磁軸承等效剛度較低時固有頻率相差不大,柔性轉子模型誤差主要來源于低等效剛度下阻尼對固有頻率的影響,此時建立電磁軸承-剛性轉子模型具有建模簡單,方便求解含阻尼固有頻率的優(yōu)勢;而當?shù)刃偠容^高時,阻尼對固有頻率影響不再明顯,由于轉子在一階平動模態(tài)即出現(xiàn)了明顯的彈性模態(tài),會導致剛性轉子模型求解不準確,應當利用柔性轉子數(shù)學模型進行分析。因此,對于數(shù)學模型的選取除了依據(jù)轉速外,也需要參考電磁軸承的等效剛度及轉子本身的彈性剛度。

    3 電磁軸承支承振動傳遞特性分析

    在旋轉機械運行過程中,轉子系統(tǒng)會在不平衡力的作用下產(chǎn)生振動,因此要通過諧響應分析得到系統(tǒng)在簡諧載荷下的穩(wěn)態(tài)響應。為模擬飛輪1的不平衡量對轉子系統(tǒng)引入的不平衡激勵,在飛輪1(圖4)位置施加554.28 N的簡諧力,求解力傳遞率、基座加速度響應以研究電磁軸承在不同等效剛度和等效阻尼下的振動傳遞特性。

    3.1 電磁軸承柔性轉子系統(tǒng)位移響應與力傳遞率

    飛輪1存在的不平衡量將為轉子系統(tǒng)引入不平衡激振力,在MATLAB中定義外部擾動力分布矩陣B1與外部擾動力f,在飛輪1中心位置對應節(jié)點施加554.28 N的擾動力,等效剛度不同取值時電磁軸承處的轉子位移響應如圖12所示。

    圖12 不同等效剛度下驅(qū)動端電磁軸承處的位移響應Fig.12 Displacement response of AMB at driving endunder different equivalent stiffnesses

    由圖12可知,隨著等效剛度的增加,位移響應有下降的趨勢,而一階固有頻率(不同剛度下一階、二階固有頻率可以通過圖中不同剛度下各自的第1,2個峰值對應橫坐標頻率獲得)有逐漸增加的趨勢。當?shù)刃偠葟?.30×106N/m增大到1.00×106N/m時,主要改變的是前2階剛性模態(tài)固有頻率,而幾乎不影響一階彎曲固有頻率,當?shù)刃偠葟?0.00×106N/m增大到50.00×106N/m時,一階固有頻率不再明顯增加,與利用有限元分析軟件得到的結果一致,即此時由于轉子在高剛度支承下一階平動模態(tài)即出現(xiàn)了彈性模態(tài),固有頻率增加不明顯。

    定義傳遞力為電磁軸承處位移響應乘以此時的等效剛度[9],則力傳遞率為

    (15)

    式中:Ff(ω)為傳遞到磁軸承的力;f(ω)為外部擾動力。

    等效剛度不同取值時電磁軸承處的力傳遞率如圖13所示。

    由圖13可知:在轉速較低接近靜態(tài)懸浮時,每套電磁軸承處力傳遞率略小于零,這是由于外界擾動力分散作用在2套電磁軸承上;增加電磁軸承的等效剛度,一階共振頻率向高頻區(qū)移動,力傳遞率一階峰值先增大后趨于穩(wěn)定,電磁軸承等效剛度降低會使振動峰值向低頻區(qū)移動,力傳遞率整體降低。

    3.2 不同支承參數(shù)下基座的加速度響應

    在轉子的飛輪1位置施加2個互相垂直的徑向不平衡激振力,幅值為554.28 N,相位相差90°,邊界條件與模態(tài)分析一致,取等效剛度為0.30×106,1.00×106,10.00×106,50.00×106N/m,設置求解頻率范圍為0~500 Hz,頻率間隔為2 Hz,采用完全法做諧響應分析,輸出驅(qū)動端電磁軸承基座處的加速度響應。

    驅(qū)動端電磁軸承基座處的加速度響應如圖14所示,隨著等效剛度的降低,基座的加速度響應也有減小的趨勢,說明隨著等效剛度的降低,轉子不平衡振動傳遞到基座的振動也隨之減小,但由于等效剛度的降低,在低頻段內(nèi)出現(xiàn)了一階共振峰值,因此等效剛度要根據(jù)旋轉設備的運行狀況進行選擇。

    圖14 不同等效剛度下驅(qū)動端電磁軸承基座的加速度響應Fig.14 Acceleration response of AMB housing withdifferent equivalent stiffnesses (DE)

    相比于傳統(tǒng)的機械軸承,電磁軸承的等效剛度和等效阻尼由控制參數(shù)決定,可以通過調(diào)整控制參數(shù),即低轉速、高等效剛度或高轉速、低等效剛度的方式減小傳遞到基座的振動,針對頻率區(qū)間內(nèi)出現(xiàn)振動峰值的問題,也可以通過調(diào)整局部等效剛度使共振峰值移動,并將其移出關心的頻率區(qū)間。

    等效阻尼分別取100,500,1 000 N·s/m,采集驅(qū)動端電磁軸承基座處的加速度響應如圖15所示:隨著等效阻尼的增加出現(xiàn)了明顯的阻尼削峰效果,當c=1 000 N·s/m時,頻率段內(nèi)響應峰值由3個減少到2個,響應峰值得到了有效控制;但在整個頻段內(nèi)轉子振動傳遞到基座引起的加速度響應增大,這是由于峰值處的能量被分散到整個頻率段,使整個頻率段加速度幅值響應增大。

    圖15 不同等效阻尼下驅(qū)動端電磁軸承基座的加速度響應Fig.15 Acceleration response of AMB housing withdifferent equivalent dampings (DE)

    4 結論

    針對變參數(shù)控制對電磁軸承等效剛度的影響,對轉子系統(tǒng)振動特性的影響進行分析,得出以下結論:

    1)將傳感器、功率放大器等效為一階慣性環(huán)節(jié),在PID控制下等效剛度主要受比例系數(shù)的影響,同時也受轉速與微分系數(shù)共同作用的影響,積分系數(shù)的影響可忽略不計。

    2)等效剛度對前3階固有頻率的影響分3個階段。當?shù)刃偠葹?.10×106N/m時,調(diào)節(jié)的是剛性固有頻率,而轉子一階彎曲固有頻率幾乎不受影響,此時振型為明顯的柔性支承-剛性轉子形式;當?shù)刃偠葹?.00×106N/m時,轉子的一階平動模態(tài)振型表現(xiàn)為明顯的彈性模態(tài),最先在轉子兩端出現(xiàn)節(jié)點,前3階固有頻率同時升高;當?shù)刃偠雀邥r,節(jié)點隨著等效剛度的升高向電磁軸承支承位置移動,一階平動固有頻率變化不再明顯,僅錐動與一階彎曲固有頻率有明顯上升。

    3)從力傳遞率及加速度響應的角度來看,降低等效剛度有利于減小傳遞到基座的振動,但共振峰值也會向低頻區(qū)移動,因此可采用低轉速、高剛度或高轉速、低剛度的方式,依據(jù)轉子轉速調(diào)節(jié)固有頻率以減小過臨界共振的影響。

    猜你喜歡
    剛性振型固有頻率
    關于模態(tài)綜合法的注記
    力學與實踐(2022年5期)2022-10-21 08:10:34
    縱向激勵下大跨鋼桁拱橋高階振型效應分析
    自我革命需要“剛性推進”
    當代陜西(2022年5期)2022-04-19 12:10:46
    現(xiàn)場測定大型水輪發(fā)電機組軸系的固有頻率
    大電機技術(2021年2期)2021-07-21 07:28:38
    加權p-Laplace型方程的剛性
    河南科學(2020年3期)2020-06-02 08:30:10
    塔腿加過渡段輸電塔動力特性分析
    特種結構(2019年2期)2019-08-19 10:05:52
    剛性兌付的法律治理
    金融法苑(2018年2期)2018-12-07 00:59:52
    結構振型幾何辨識及應用研究
    山西建筑(2015年14期)2015-06-05 09:37:07
    總溫總壓測頭模態(tài)振型變化規(guī)律研究
    A novel functional electrical stimulation-control system for restoring motor function of post-stroke hemiplegic patients
    videosex国产| 咕卡用的链子| 99精品久久久久人妻精品| 亚洲国产日韩一区二区| 美女视频免费永久观看网站| 亚洲欧美清纯卡通| 天堂8中文在线网| 亚洲免费av在线视频| 精品国产乱码久久久久久男人| 欧美在线黄色| 亚洲国产成人一精品久久久| 中文字幕精品免费在线观看视频| xxxhd国产人妻xxx| 欧美少妇被猛烈插入视频| av片东京热男人的天堂| 在线观看免费日韩欧美大片| 性高湖久久久久久久久免费观看| 久久婷婷青草| 亚洲专区中文字幕在线 | 水蜜桃什么品种好| 男女国产视频网站| 女人精品久久久久毛片| 国产精品 国内视频| 亚洲人成77777在线视频| 欧美黄色片欧美黄色片| 丁香六月天网| 美女国产高潮福利片在线看| 久久国产亚洲av麻豆专区| 国产精品 欧美亚洲| 99热网站在线观看| 大香蕉久久成人网| 国产高清不卡午夜福利| 欧美黑人精品巨大| 久久久久久人妻| 青春草国产在线视频| 咕卡用的链子| 久久97久久精品| 亚洲av电影在线观看一区二区三区| av又黄又爽大尺度在线免费看| 国产精品人妻久久久影院| 精品久久久久久电影网| 黄片无遮挡物在线观看| 制服丝袜香蕉在线| 老汉色av国产亚洲站长工具| 亚洲av成人不卡在线观看播放网 | 高清视频免费观看一区二区| 丝袜喷水一区| 在线亚洲精品国产二区图片欧美| 日本欧美国产在线视频| 国产1区2区3区精品| 伦理电影免费视频| 爱豆传媒免费全集在线观看| 亚洲欧洲国产日韩| 成年美女黄网站色视频大全免费| 人人妻人人澡人人爽人人夜夜| 不卡视频在线观看欧美| 人人妻人人爽人人添夜夜欢视频| 欧美激情 高清一区二区三区| 国产午夜精品一二区理论片| 人人妻,人人澡人人爽秒播 | 狂野欧美激情性xxxx| 桃花免费在线播放| 一区二区av电影网| 国产亚洲欧美精品永久| 操出白浆在线播放| 波多野结衣一区麻豆| 美女主播在线视频| 在线观看免费午夜福利视频| 黄片无遮挡物在线观看| 久久性视频一级片| 亚洲av国产av综合av卡| 国产精品 国内视频| 亚洲七黄色美女视频| 日韩一本色道免费dvd| 亚洲七黄色美女视频| 大陆偷拍与自拍| 色网站视频免费| 久久国产精品大桥未久av| 久久鲁丝午夜福利片| 男女床上黄色一级片免费看| 亚洲第一青青草原| 亚洲欧美成人精品一区二区| 成人午夜精彩视频在线观看| 欧美黑人精品巨大| 国产毛片在线视频| 在线精品无人区一区二区三| 丁香六月天网| 亚洲美女搞黄在线观看| 成年av动漫网址| 久久久久久人妻| 韩国av在线不卡| 国产日韩欧美视频二区| 日韩 亚洲 欧美在线| 美女视频免费永久观看网站| 成年av动漫网址| 九九爱精品视频在线观看| 久久久久久久大尺度免费视频| 亚洲国产精品成人久久小说| 欧美乱码精品一区二区三区| 丰满乱子伦码专区| av一本久久久久| 黑人欧美特级aaaaaa片| 99热网站在线观看| 日本wwww免费看| 国产精品成人在线| 视频在线观看一区二区三区| 只有这里有精品99| 国产精品一国产av| 国产黄色免费在线视频| 日韩av在线免费看完整版不卡| 中文字幕另类日韩欧美亚洲嫩草| 成年人免费黄色播放视频| 亚洲国产日韩一区二区| 成人亚洲精品一区在线观看| 国产精品.久久久| 免费观看av网站的网址| 大话2 男鬼变身卡| 欧美激情高清一区二区三区 | av片东京热男人的天堂| 最新在线观看一区二区三区 | 午夜福利乱码中文字幕| 亚洲第一av免费看| 久久久精品免费免费高清| 午夜福利,免费看| 久久久精品免费免费高清| 午夜免费观看性视频| av视频免费观看在线观看| 男女床上黄色一级片免费看| kizo精华| 亚洲人成77777在线视频| 电影成人av| 丝袜人妻中文字幕| 亚洲精品自拍成人| 日韩人妻精品一区2区三区| 天天躁夜夜躁狠狠久久av| 街头女战士在线观看网站| 亚洲美女黄色视频免费看| 国产精品成人在线| 69精品国产乱码久久久| 桃花免费在线播放| 久久久久国产一级毛片高清牌| 男女边吃奶边做爰视频| 91精品伊人久久大香线蕉| 9191精品国产免费久久| 超碰成人久久| 建设人人有责人人尽责人人享有的| 国产精品久久久久久精品古装| 三上悠亚av全集在线观看| 人人妻,人人澡人人爽秒播 | 女人爽到高潮嗷嗷叫在线视频| 最近手机中文字幕大全| 一边摸一边抽搐一进一出视频| 18禁动态无遮挡网站| 高清在线视频一区二区三区| 午夜福利一区二区在线看| 国产毛片在线视频| 一级毛片黄色毛片免费观看视频| 精品福利永久在线观看| 观看av在线不卡| 操美女的视频在线观看| 亚洲成人国产一区在线观看 | 午夜福利视频在线观看免费| 日本欧美视频一区| 亚洲视频免费观看视频| 国产成人一区二区在线| 亚洲av电影在线进入| 男的添女的下面高潮视频| 丝袜人妻中文字幕| 亚洲婷婷狠狠爱综合网| 十八禁人妻一区二区| 久久精品亚洲av国产电影网| 亚洲天堂av无毛| 哪个播放器可以免费观看大片| 看非洲黑人一级黄片| 最新在线观看一区二区三区 | videosex国产| 国产伦人伦偷精品视频| 久久精品国产亚洲av高清一级| 欧美亚洲 丝袜 人妻 在线| 免费观看a级毛片全部| 哪个播放器可以免费观看大片| 精品一区二区三区四区五区乱码 | 男女床上黄色一级片免费看| 亚洲成色77777| 男的添女的下面高潮视频| 亚洲熟女精品中文字幕| 极品人妻少妇av视频| 男女国产视频网站| 9191精品国产免费久久| 国产老妇伦熟女老妇高清| 老司机靠b影院| 一区二区日韩欧美中文字幕| 国产在线免费精品| 中文天堂在线官网| 男人操女人黄网站| 亚洲成国产人片在线观看| 黄色怎么调成土黄色| 国产片特级美女逼逼视频| 亚洲av综合色区一区| 精品亚洲乱码少妇综合久久| 国产日韩欧美在线精品| 久久精品国产亚洲av涩爱| a级毛片在线看网站| 日本午夜av视频| 久久 成人 亚洲| 久久久精品区二区三区| 久久精品国产亚洲av涩爱| 亚洲精品久久久久久婷婷小说| 精品一区二区免费观看| 一级毛片电影观看| 久久精品久久精品一区二区三区| 亚洲欧美清纯卡通| 久久精品久久久久久久性| 国产色婷婷99| 久久女婷五月综合色啪小说| 久久99精品国语久久久| 中文字幕色久视频| 精品一区二区三区av网在线观看 | 天天操日日干夜夜撸| 国产精品无大码| 国产激情久久老熟女| 亚洲国产最新在线播放| 水蜜桃什么品种好| 亚洲少妇的诱惑av| 亚洲精品一二三| 中文字幕另类日韩欧美亚洲嫩草| 欧美精品一区二区免费开放| 一本色道久久久久久精品综合| 欧美久久黑人一区二区| 国产精品久久久久成人av| 欧美日韩亚洲高清精品| 男女无遮挡免费网站观看| 日韩,欧美,国产一区二区三区| 亚洲精品aⅴ在线观看| 在线观看免费午夜福利视频| 国产精品欧美亚洲77777| 天天躁日日躁夜夜躁夜夜| 丰满饥渴人妻一区二区三| 日韩欧美一区视频在线观看| 免费观看人在逋| 精品国产一区二区久久| 成人18禁高潮啪啪吃奶动态图| a级毛片黄视频| 亚洲精品国产av蜜桃| 亚洲精华国产精华液的使用体验| 香蕉丝袜av| 欧美人与性动交α欧美精品济南到| 午夜激情久久久久久久| 嫩草影院入口| 亚洲成人一二三区av| 宅男免费午夜| 国产女主播在线喷水免费视频网站| 极品人妻少妇av视频| av一本久久久久| 国产在视频线精品| 久久久久久人人人人人| 国产不卡av网站在线观看| 永久免费av网站大全| 亚洲精品,欧美精品| 国产亚洲欧美精品永久| 欧美日韩国产mv在线观看视频| 亚洲一区中文字幕在线| 久久久久久人妻| 99国产精品免费福利视频| 久久久久久久久久久久大奶| 毛片一级片免费看久久久久| 一本久久精品| 悠悠久久av| 久久精品aⅴ一区二区三区四区| 国产有黄有色有爽视频| 精品国产一区二区久久| 久久久久人妻精品一区果冻| 一二三四在线观看免费中文在| 午夜影院在线不卡| 深夜精品福利| 中文天堂在线官网| 可以免费在线观看a视频的电影网站 | 十八禁人妻一区二区| 夫妻午夜视频| 一级黄片播放器| 看十八女毛片水多多多| xxx大片免费视频| 午夜福利视频精品| 欧美在线黄色| 欧美人与性动交α欧美精品济南到| 亚洲成人手机| 午夜福利一区二区在线看| 国产精品秋霞免费鲁丝片| 午夜激情av网站| 熟女av电影| 亚洲欧洲国产日韩| 高清欧美精品videossex| 欧美黑人精品巨大| 国产精品一国产av| 精品少妇久久久久久888优播| 婷婷成人精品国产| 熟女av电影| 黄片无遮挡物在线观看| 欧美日韩视频高清一区二区三区二| 国产日韩欧美视频二区| 欧美变态另类bdsm刘玥| 国产精品 国内视频| 国产激情久久老熟女| 成年人午夜在线观看视频| www.自偷自拍.com| 亚洲精品美女久久av网站| 日韩一本色道免费dvd| 色综合欧美亚洲国产小说| 亚洲自偷自拍图片 自拍| 伊人久久国产一区二区| 五月开心婷婷网| 91aial.com中文字幕在线观看| 日韩中文字幕视频在线看片| 午夜福利影视在线免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲欧美一区二区三区黑人| 国产精品人妻久久久影院| 亚洲自偷自拍图片 自拍| 亚洲精品一区蜜桃| 久久久国产精品麻豆| 亚洲国产av新网站| 美女福利国产在线| 亚洲精品国产av蜜桃| 不卡视频在线观看欧美| 最近最新中文字幕大全免费视频 | 香蕉国产在线看| 天堂8中文在线网| 国产一级毛片在线| 男女床上黄色一级片免费看| 午夜av观看不卡| 国产免费视频播放在线视频| 999久久久国产精品视频| 夫妻性生交免费视频一级片| 日韩大片免费观看网站| 亚洲欧美清纯卡通| 亚洲伊人色综图| 久久鲁丝午夜福利片| 两性夫妻黄色片| 久久精品熟女亚洲av麻豆精品| 久久人妻熟女aⅴ| 精品午夜福利在线看| 一区二区三区乱码不卡18| 丝瓜视频免费看黄片| 国产精品av久久久久免费| 一本—道久久a久久精品蜜桃钙片| 精品久久久精品久久久| 日日啪夜夜爽| 国产精品av久久久久免费| 精品一区二区免费观看| 人妻人人澡人人爽人人| 精品国产一区二区三区四区第35| 国产成人精品无人区| 国产精品亚洲av一区麻豆 | 熟女av电影| 中文欧美无线码| 中文字幕av电影在线播放| 18禁动态无遮挡网站| 成年人免费黄色播放视频| 人人妻人人澡人人看| 热re99久久国产66热| 在线观看人妻少妇| 美女大奶头黄色视频| 国产人伦9x9x在线观看| 亚洲欧美激情在线| 国产伦理片在线播放av一区| 国产精品一国产av| 天天躁日日躁夜夜躁夜夜| 在线观看免费视频网站a站| 午夜日本视频在线| 免费观看性生交大片5| 伊人久久国产一区二区| 国产成人精品在线电影| 色婷婷av一区二区三区视频| 午夜福利影视在线免费观看| 午夜福利乱码中文字幕| 麻豆av在线久日| 青春草亚洲视频在线观看| 激情视频va一区二区三区| 看免费成人av毛片| 久久久久视频综合| 韩国av在线不卡| 大香蕉久久网| 日韩不卡一区二区三区视频在线| 我要看黄色一级片免费的| 少妇被粗大的猛进出69影院| 婷婷色综合大香蕉| 精品久久久精品久久久| 日韩 欧美 亚洲 中文字幕| 国产精品国产三级专区第一集| 黄色毛片三级朝国网站| tube8黄色片| 国产精品国产三级专区第一集| 久久青草综合色| 五月开心婷婷网| 99精品久久久久人妻精品| 好男人视频免费观看在线| 极品少妇高潮喷水抽搐| 精品国产一区二区三区久久久樱花| 2018国产大陆天天弄谢| 婷婷成人精品国产| 欧美日韩一区二区视频在线观看视频在线| 国产激情久久老熟女| 99热全是精品| 九草在线视频观看| 中文字幕色久视频| 免费看av在线观看网站| av有码第一页| av网站在线播放免费| 亚洲五月色婷婷综合| 黄色一级大片看看| 爱豆传媒免费全集在线观看| 中文乱码字字幕精品一区二区三区| 精品少妇黑人巨大在线播放| 日本vs欧美在线观看视频| 亚洲欧美中文字幕日韩二区| 国产成人精品福利久久| 熟女少妇亚洲综合色aaa.| 婷婷成人精品国产| 国产熟女欧美一区二区| av一本久久久久| 亚洲精品一区蜜桃| 亚洲国产欧美一区二区综合| 国产成人av激情在线播放| 极品少妇高潮喷水抽搐| 久久免费观看电影| 丁香六月欧美| 两个人免费观看高清视频| 免费在线观看黄色视频的| 高清欧美精品videossex| 尾随美女入室| 日韩一卡2卡3卡4卡2021年| 黄色视频在线播放观看不卡| 高清视频免费观看一区二区| 少妇人妻久久综合中文| 国产精品99久久99久久久不卡 | 大码成人一级视频| 悠悠久久av| 男女下面插进去视频免费观看| 日韩成人av中文字幕在线观看| 国产野战对白在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲国产最新在线播放| 18禁国产床啪视频网站| 99香蕉大伊视频| 国产精品熟女久久久久浪| 老汉色av国产亚洲站长工具| 亚洲av福利一区| 妹子高潮喷水视频| 叶爱在线成人免费视频播放| √禁漫天堂资源中文www| 18禁观看日本| 亚洲精品自拍成人| 色播在线永久视频| 久久久欧美国产精品| 亚洲一卡2卡3卡4卡5卡精品中文| 热99国产精品久久久久久7| 人妻一区二区av| 免费不卡黄色视频| 精品国产国语对白av| a 毛片基地| 国产爽快片一区二区三区| 狠狠婷婷综合久久久久久88av| 久久影院123| av在线app专区| 黄片无遮挡物在线观看| 超碰成人久久| 久久久久久久久久久免费av| 欧美精品av麻豆av| 无遮挡黄片免费观看| 观看av在线不卡| 最近2019中文字幕mv第一页| 99国产综合亚洲精品| 人人澡人人妻人| 亚洲国产看品久久| 免费在线观看黄色视频的| 国产精品国产三级国产专区5o| 精品国产露脸久久av麻豆| 亚洲精品久久午夜乱码| 精品国产一区二区三区四区第35| 亚洲欧洲国产日韩| 亚洲三区欧美一区| 天天操日日干夜夜撸| 咕卡用的链子| 日韩,欧美,国产一区二区三区| 欧美中文综合在线视频| 男女午夜视频在线观看| 街头女战士在线观看网站| 美女国产高潮福利片在线看| 久久久国产精品麻豆| 色视频在线一区二区三区| 亚洲精品一二三| √禁漫天堂资源中文www| 精品亚洲乱码少妇综合久久| 国产精品成人在线| 老司机在亚洲福利影院| 国产精品av久久久久免费| xxxhd国产人妻xxx| 考比视频在线观看| av又黄又爽大尺度在线免费看| 亚洲精品国产av蜜桃| 大香蕉久久成人网| 狂野欧美激情性bbbbbb| 国产成人精品久久久久久| 91精品国产国语对白视频| 久久久久精品人妻al黑| 久久久久人妻精品一区果冻| 成年动漫av网址| 成人国语在线视频| 亚洲精品一区蜜桃| 日本午夜av视频| 亚洲欧美中文字幕日韩二区| 香蕉丝袜av| 9191精品国产免费久久| 老汉色av国产亚洲站长工具| 国产片内射在线| h视频一区二区三区| 国产成人欧美在线观看 | 国产成人欧美在线观看 | 多毛熟女@视频| 精品一区二区三区av网在线观看 | 色播在线永久视频| 亚洲视频免费观看视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲少妇的诱惑av| 免费在线观看黄色视频的| 激情视频va一区二区三区| 亚洲一区二区三区欧美精品| 少妇人妻精品综合一区二区| 天天操日日干夜夜撸| 国产在视频线精品| 国产av精品麻豆| 99久久综合免费| 欧美激情高清一区二区三区 | 一级毛片黄色毛片免费观看视频| 9191精品国产免费久久| 成人亚洲欧美一区二区av| 国语对白做爰xxxⅹ性视频网站| 男女下面插进去视频免费观看| 国产精品国产三级国产专区5o| 伊人亚洲综合成人网| 我要看黄色一级片免费的| 丝袜美腿诱惑在线| 99久久综合免费| 伦理电影大哥的女人| 国产一区二区 视频在线| 国产成人午夜福利电影在线观看| √禁漫天堂资源中文www| 一区二区日韩欧美中文字幕| av网站在线播放免费| 美女中出高潮动态图| 精品午夜福利在线看| 亚洲第一av免费看| 国产男人的电影天堂91| 日本91视频免费播放| 久久久久视频综合| 波多野结衣av一区二区av| 亚洲av男天堂| 亚洲少妇的诱惑av| 99久久精品国产亚洲精品| 国产乱人偷精品视频| 99热网站在线观看| 丁香六月天网| 成年女人毛片免费观看观看9 | 伊人亚洲综合成人网| 久久精品亚洲熟妇少妇任你| 亚洲欧美一区二区三区黑人| 2021少妇久久久久久久久久久| 欧美日韩亚洲综合一区二区三区_| 午夜福利影视在线免费观看| 成年美女黄网站色视频大全免费| 老司机影院成人| 欧美在线一区亚洲| 久久久久国产精品人妻一区二区| 久久久亚洲精品成人影院| 水蜜桃什么品种好| 久久久亚洲精品成人影院| 亚洲精品第二区| 一区二区三区精品91| 精品少妇一区二区三区视频日本电影 | 亚洲久久久国产精品| 我的亚洲天堂| 不卡视频在线观看欧美| 黄色毛片三级朝国网站| 久久久国产精品麻豆| 热re99久久国产66热| 亚洲av电影在线进入| netflix在线观看网站| 中文字幕另类日韩欧美亚洲嫩草| 我要看黄色一级片免费的| 亚洲欧美精品自产自拍| 女性被躁到高潮视频| 精品国产露脸久久av麻豆| 美女福利国产在线| 久久久久久久久久久免费av| 精品国产国语对白av| 美女脱内裤让男人舔精品视频| 午夜福利影视在线免费观看| 亚洲七黄色美女视频| 黄色视频在线播放观看不卡| 日韩一本色道免费dvd| 国产日韩欧美视频二区| 青草久久国产| 久久久久久人人人人人| 亚洲美女视频黄频| 九九爱精品视频在线观看| 99久久精品国产亚洲精品| 天天躁狠狠躁夜夜躁狠狠躁| 中国三级夫妇交换| 亚洲第一av免费看| 久久ye,这里只有精品| 国产伦理片在线播放av一区| 久久久久久久精品精品| 成人免费观看视频高清| 亚洲欧美一区二区三区久久| av.在线天堂|