• <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
    黑人巨大精品欧美一区二区蜜桃| 一级毛片女人18水好多| 99国产精品一区二区蜜桃av | 在线观看舔阴道视频| 久久九九热精品免费| 成人亚洲精品一区在线观看| 丝袜在线中文字幕| 午夜激情av网站| 欧美精品啪啪一区二区三区 | 动漫黄色视频在线观看| 国产日韩欧美视频二区| 岛国在线观看网站| 人妻人人澡人人爽人人| 天天躁日日躁夜夜躁夜夜| 精品少妇一区二区三区视频日本电影| 老汉色∧v一级毛片| 黑人巨大精品欧美一区二区蜜桃| 亚洲欧美色中文字幕在线| av视频免费观看在线观看| 国产精品久久久av美女十八| 日韩欧美一区视频在线观看| 波多野结衣av一区二区av| 丁香六月天网| 青春草视频在线免费观看| 国产精品国产av在线观看| 国产精品国产av在线观看| 十八禁网站网址无遮挡| 国产精品影院久久| 亚洲国产欧美网| av在线播放精品| 黄色怎么调成土黄色| 免费久久久久久久精品成人欧美视频| 极品人妻少妇av视频| 日韩制服骚丝袜av| 国产区一区二久久| 久久综合国产亚洲精品| 久9热在线精品视频| 日本一区二区免费在线视频| 精品熟女少妇八av免费久了| 人人妻人人澡人人看| 亚洲人成电影免费在线| 一本—道久久a久久精品蜜桃钙片| 国产麻豆69| av又黄又爽大尺度在线免费看| 欧美另类亚洲清纯唯美| 大香蕉久久网| 欧美精品一区二区免费开放| 真人做人爱边吃奶动态| 在线观看人妻少妇| 大型av网站在线播放| 成年女人毛片免费观看观看9 | 亚洲七黄色美女视频| 久久久久久亚洲精品国产蜜桃av| 欧美激情久久久久久爽电影 | 国产福利在线免费观看视频| 精品一区在线观看国产| 国产精品久久久久久精品电影小说| 亚洲第一av免费看| 亚洲七黄色美女视频| 中文精品一卡2卡3卡4更新| 夜夜骑夜夜射夜夜干| videos熟女内射| 99精国产麻豆久久婷婷| 久久久久国产一级毛片高清牌| 国产成人a∨麻豆精品| 好男人电影高清在线观看| 美女高潮到喷水免费观看| 最黄视频免费看| av国产精品久久久久影院| 老汉色∧v一级毛片| 国产一区二区三区综合在线观看| 午夜精品久久久久久毛片777| 久久99热这里只频精品6学生| 麻豆乱淫一区二区| 久久久水蜜桃国产精品网| 成人黄色视频免费在线看| 色视频在线一区二区三区| 999久久久精品免费观看国产| 国产精品二区激情视频| 日韩制服骚丝袜av| av电影中文网址| 无遮挡黄片免费观看| 国产男女超爽视频在线观看| av网站在线播放免费| 老熟妇乱子伦视频在线观看 | 成人国语在线视频| 国产精品一区二区精品视频观看| 99精品欧美一区二区三区四区| 国产精品成人在线| 午夜福利视频精品| 国产深夜福利视频在线观看| 日韩视频一区二区在线观看| 亚洲国产中文字幕在线视频| 国产视频一区二区在线看| 黄片小视频在线播放| 制服诱惑二区| 少妇猛男粗大的猛烈进出视频| 国产欧美日韩一区二区三区在线| 精品国产乱码久久久久久男人| av又黄又爽大尺度在线免费看| 老熟妇乱子伦视频在线观看 | 可以免费在线观看a视频的电影网站| av在线播放精品| 久久久久网色| 每晚都被弄得嗷嗷叫到高潮| 亚洲欧美色中文字幕在线| 考比视频在线观看| 岛国毛片在线播放| 国产福利在线免费观看视频| 欧美人与性动交α欧美精品济南到| 最新在线观看一区二区三区| 亚洲欧美日韩高清在线视频 | 亚洲中文字幕日韩| 丁香六月天网| 国产伦理片在线播放av一区| 男女免费视频国产| 黄色片一级片一级黄色片| 午夜久久久在线观看| 欧美日韩国产mv在线观看视频| 热re99久久精品国产66热6| 一区二区三区激情视频| 一级片免费观看大全| 人人妻人人澡人人爽人人夜夜| 亚洲一码二码三码区别大吗| 99热国产这里只有精品6| 法律面前人人平等表现在哪些方面 | 90打野战视频偷拍视频| 汤姆久久久久久久影院中文字幕| 在线av久久热| av国产精品久久久久影院| 欧美精品亚洲一区二区| 18禁观看日本| 久久亚洲国产成人精品v| 亚洲av欧美aⅴ国产| 最新在线观看一区二区三区| 欧美另类一区| 久久女婷五月综合色啪小说| 成年美女黄网站色视频大全免费| 女人久久www免费人成看片| 欧美激情极品国产一区二区三区| 欧美精品人与动牲交sv欧美| 日本vs欧美在线观看视频| 淫妇啪啪啪对白视频 | av线在线观看网站| 欧美午夜高清在线| 精品久久久久久久毛片微露脸 | 日日摸夜夜添夜夜添小说| 亚洲专区国产一区二区| 亚洲欧美色中文字幕在线| 亚洲av电影在线进入| 欧美日韩一级在线毛片| 搡老乐熟女国产| 免费在线观看视频国产中文字幕亚洲 | 亚洲成人免费av在线播放| 日本精品一区二区三区蜜桃| 人成视频在线观看免费观看| 久久这里只有精品19| 淫妇啪啪啪对白视频 | av视频免费观看在线观看| 俄罗斯特黄特色一大片| 两人在一起打扑克的视频| 777久久人妻少妇嫩草av网站| 高清视频免费观看一区二区| 老司机午夜十八禁免费视频| 一级黄色大片毛片| 最近中文字幕2019免费版| 国产激情久久老熟女| 曰老女人黄片| 老汉色av国产亚洲站长工具| 最新在线观看一区二区三区| 亚洲欧美激情在线| 欧美日韩亚洲综合一区二区三区_| 日本猛色少妇xxxxx猛交久久| 久久久精品区二区三区| 国产福利在线免费观看视频| 免费在线观看日本一区| 久热这里只有精品99| 欧美亚洲 丝袜 人妻 在线| 久久久久精品国产欧美久久久 | 12—13女人毛片做爰片一| 久久精品aⅴ一区二区三区四区| 国产福利在线免费观看视频| 99热网站在线观看| 国产成人精品久久二区二区免费| 精品视频人人做人人爽| 老司机午夜十八禁免费视频| 少妇 在线观看| 岛国毛片在线播放| 国产人伦9x9x在线观看| 成人av一区二区三区在线看 | 国产精品一区二区在线不卡| 男人爽女人下面视频在线观看| 亚洲伊人色综图| 国产在线视频一区二区| 欧美亚洲 丝袜 人妻 在线| 法律面前人人平等表现在哪些方面 | 久久精品国产a三级三级三级| 人妻一区二区av| 美女大奶头黄色视频| 成人国语在线视频| 免费在线观看视频国产中文字幕亚洲 | 亚洲中文av在线| 超色免费av| 国产激情久久老熟女| 国产精品久久久久成人av| 好男人电影高清在线观看| 91精品三级在线观看| 久久久久国内视频| 午夜激情久久久久久久| 国产成+人综合+亚洲专区| 欧美+亚洲+日韩+国产| 精品人妻熟女毛片av久久网站| 一级片免费观看大全| 777米奇影视久久| 午夜91福利影院| 999久久久国产精品视频| 亚洲国产精品999| 亚洲国产欧美一区二区综合| 国产男人的电影天堂91| 欧美国产精品va在线观看不卡| 天堂俺去俺来也www色官网| 国产日韩一区二区三区精品不卡| 黄色怎么调成土黄色| 美女福利国产在线| 成年av动漫网址| 亚洲成av片中文字幕在线观看| 亚洲精品日韩在线中文字幕| 国产男人的电影天堂91| 免费在线观看视频国产中文字幕亚洲 | 久久天堂一区二区三区四区| 最近最新免费中文字幕在线| 国产成人啪精品午夜网站| 一本久久精品| 精品亚洲成a人片在线观看| 国产免费现黄频在线看| 日本vs欧美在线观看视频| 欧美日韩国产mv在线观看视频| 好男人电影高清在线观看| 免费高清在线观看日韩| 一级毛片电影观看| 久久国产精品男人的天堂亚洲| 午夜精品久久久久久毛片777| 老熟妇乱子伦视频在线观看 | 国产一区二区 视频在线| 在线观看人妻少妇| 日韩有码中文字幕| 韩国高清视频一区二区三区| 亚洲国产欧美一区二区综合| 老司机亚洲免费影院| 日韩,欧美,国产一区二区三区| 亚洲国产欧美在线一区| 女人高潮潮喷娇喘18禁视频| 免费在线观看黄色视频的| 久热爱精品视频在线9| 国内毛片毛片毛片毛片毛片| 国产片内射在线| 建设人人有责人人尽责人人享有的| 美女扒开内裤让男人捅视频| 黄色片一级片一级黄色片| 亚洲avbb在线观看| av不卡在线播放| 美女午夜性视频免费| 美女大奶头黄色视频| 婷婷丁香在线五月| 亚洲国产精品成人久久小说| 美女中出高潮动态图| 交换朋友夫妻互换小说| 自拍欧美九色日韩亚洲蝌蚪91| 一边摸一边抽搐一进一出视频| 国产主播在线观看一区二区| 午夜福利在线观看吧| 久久久久久久久免费视频了| 国产高清视频在线播放一区 | 精品久久蜜臀av无| 九色亚洲精品在线播放| 国产欧美日韩精品亚洲av| 国产一区二区在线观看av| 久久精品亚洲熟妇少妇任你| 热re99久久精品国产66热6| 亚洲一区中文字幕在线| 国产亚洲av片在线观看秒播厂| a级毛片在线看网站| 欧美少妇被猛烈插入视频| 中文字幕精品免费在线观看视频| 色婷婷久久久亚洲欧美| 国产三级黄色录像| 久久久久久久久免费视频了| 亚洲第一青青草原| 欧美精品av麻豆av| 国产免费av片在线观看野外av| 亚洲精品乱久久久久久| 成年av动漫网址| 免费人妻精品一区二区三区视频| 亚洲一区中文字幕在线| 亚洲七黄色美女视频| 王馨瑶露胸无遮挡在线观看| 性高湖久久久久久久久免费观看| 丰满迷人的少妇在线观看| 国产精品香港三级国产av潘金莲| 亚洲av男天堂| 女人高潮潮喷娇喘18禁视频| 亚洲人成电影免费在线| 久久精品亚洲熟妇少妇任你| 久久毛片免费看一区二区三区| 欧美成狂野欧美在线观看| 一级毛片精品| 亚洲欧美色中文字幕在线| 母亲3免费完整高清在线观看| √禁漫天堂资源中文www| 丝袜脚勾引网站| 亚洲第一青青草原| 精品国产一区二区久久| 777久久人妻少妇嫩草av网站| 欧美黄色片欧美黄色片| 最新在线观看一区二区三区| 久久国产精品影院| 国产在线免费精品| 国产精品免费大片| 伊人久久大香线蕉亚洲五| 国产免费福利视频在线观看| 黄色怎么调成土黄色| 精品人妻1区二区| 亚洲av欧美aⅴ国产| 亚洲欧美日韩高清在线视频 | 久久国产精品大桥未久av| 啦啦啦中文免费视频观看日本| a级片在线免费高清观看视频| 国产精品免费大片| 成年动漫av网址| 成人国产av品久久久| 美女午夜性视频免费| 极品人妻少妇av视频| 欧美精品亚洲一区二区| 男女边摸边吃奶| 一边摸一边抽搐一进一出视频| 51午夜福利影视在线观看| 丰满少妇做爰视频| 亚洲av电影在线进入| 久久 成人 亚洲| 十八禁网站网址无遮挡| 亚洲精品久久午夜乱码| 国产精品九九99| av国产精品久久久久影院| 亚洲久久久国产精品| 国产一区二区 视频在线| 考比视频在线观看| 亚洲欧洲精品一区二区精品久久久| 午夜福利乱码中文字幕| 午夜成年电影在线免费观看| av视频免费观看在线观看| √禁漫天堂资源中文www| 午夜福利乱码中文字幕| www.熟女人妻精品国产| 午夜福利乱码中文字幕| 国产精品 欧美亚洲| 欧美精品高潮呻吟av久久| 黄色片一级片一级黄色片| 久久久久久久国产电影| 视频区图区小说| 纯流量卡能插随身wifi吗| 别揉我奶头~嗯~啊~动态视频 | 亚洲免费av在线视频| 国产精品影院久久| 菩萨蛮人人尽说江南好唐韦庄| 乱人伦中国视频| 亚洲欧美成人综合另类久久久| 超碰97精品在线观看| 欧美xxⅹ黑人| 成人三级做爰电影| 制服诱惑二区| 脱女人内裤的视频| 国产成人系列免费观看| 免费女性裸体啪啪无遮挡网站| svipshipincom国产片| 蜜桃国产av成人99| 90打野战视频偷拍视频| 日韩免费高清中文字幕av| 午夜福利免费观看在线| 国产在视频线精品| 最近最新中文字幕大全免费视频| 下体分泌物呈黄色| 十八禁人妻一区二区| 激情视频va一区二区三区| 久久九九热精品免费| 蜜桃国产av成人99| 国产成人一区二区三区免费视频网站| 涩涩av久久男人的天堂| 自拍欧美九色日韩亚洲蝌蚪91| 日韩视频在线欧美| 日本一区二区免费在线视频| 新久久久久国产一级毛片| 欧美黄色淫秽网站| 日韩精品免费视频一区二区三区| 成人影院久久| 十八禁网站免费在线| 日日夜夜操网爽| 成人免费观看视频高清| 美国免费a级毛片| 免费女性裸体啪啪无遮挡网站| 久久精品亚洲熟妇少妇任你| 色精品久久人妻99蜜桃| 国产无遮挡羞羞视频在线观看| 香蕉国产在线看| 飞空精品影院首页| 日韩制服骚丝袜av| 久久国产精品男人的天堂亚洲| 美女扒开内裤让男人捅视频| 欧美激情 高清一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 国产免费现黄频在线看| 超碰成人久久| 亚洲国产欧美一区二区综合| 妹子高潮喷水视频| 久久这里只有精品19| 亚洲少妇的诱惑av| 欧美日韩亚洲国产一区二区在线观看 | 电影成人av| av电影中文网址| 欧美成人午夜精品| 亚洲精品一区蜜桃| 国产成人精品无人区| 久久久久久久精品精品| 女人高潮潮喷娇喘18禁视频| 美女大奶头黄色视频| 成人影院久久| 天天躁狠狠躁夜夜躁狠狠躁| bbb黄色大片| 999精品在线视频| 欧美日韩视频精品一区| 国产在线免费精品| 黑人猛操日本美女一级片| 青春草视频在线免费观看| 国产一区二区在线观看av| 好男人电影高清在线观看| 视频区图区小说| 成人国产一区最新在线观看| 成人手机av| 亚洲三区欧美一区| 永久免费av网站大全| 久久久久网色| 日本a在线网址| 9191精品国产免费久久| 十八禁网站网址无遮挡| 亚洲精品国产色婷婷电影| 久久国产精品男人的天堂亚洲| 三级毛片av免费| 亚洲欧美激情在线| 亚洲精品中文字幕一二三四区 | 国产片内射在线| av欧美777| 可以免费在线观看a视频的电影网站| 十八禁高潮呻吟视频| 成人黄色视频免费在线看| 国产免费一区二区三区四区乱码| 18禁国产床啪视频网站| 国产亚洲欧美精品永久| 丰满饥渴人妻一区二区三| 免费观看人在逋| 中文字幕高清在线视频| 国产伦理片在线播放av一区| 亚洲国产中文字幕在线视频| av在线老鸭窝| 性少妇av在线| 三上悠亚av全集在线观看| 国产精品偷伦视频观看了| 一二三四在线观看免费中文在| 男女下面插进去视频免费观看| 欧美 亚洲 国产 日韩一| 男女边摸边吃奶| 丝袜在线中文字幕| 在线av久久热| 亚洲视频免费观看视频| 91精品伊人久久大香线蕉| 51午夜福利影视在线观看| 亚洲精品成人av观看孕妇| 日本av免费视频播放| av天堂在线播放| 2018国产大陆天天弄谢| 国产欧美日韩精品亚洲av| 精品久久久精品久久久| 久久精品成人免费网站| 亚洲专区中文字幕在线| 亚洲第一欧美日韩一区二区三区 | 十八禁网站免费在线| 国产成人精品久久二区二区免费| 黄频高清免费视频| 久久久精品区二区三区| 99久久精品国产亚洲精品| 男女无遮挡免费网站观看| 久久午夜综合久久蜜桃| 免费高清在线观看日韩| 一级a爱视频在线免费观看| 免费观看a级毛片全部| 无遮挡黄片免费观看| 欧美精品av麻豆av| 麻豆国产av国片精品| 亚洲欧美清纯卡通| 啦啦啦啦在线视频资源| 久久精品亚洲av国产电影网| 亚洲欧美精品综合一区二区三区| 在线看a的网站| 岛国在线观看网站| 中文字幕另类日韩欧美亚洲嫩草| 王馨瑶露胸无遮挡在线观看| 久久中文字幕一级| 女人爽到高潮嗷嗷叫在线视频| 91精品三级在线观看| 亚洲中文日韩欧美视频| 成人亚洲精品一区在线观看| 制服人妻中文乱码| 日韩欧美一区二区三区在线观看 | 免费在线观看日本一区| 欧美变态另类bdsm刘玥| 中文精品一卡2卡3卡4更新| 欧美在线黄色| 国产成人欧美在线观看 | 免费不卡黄色视频| 999久久久精品免费观看国产| 中文字幕另类日韩欧美亚洲嫩草| 丝袜人妻中文字幕| 99国产精品一区二区三区| 少妇裸体淫交视频免费看高清 | 国产精品 国内视频| 9热在线视频观看99| 亚洲精华国产精华精| 纯流量卡能插随身wifi吗| av网站在线播放免费| 亚洲精品美女久久av网站| 欧美精品啪啪一区二区三区 | 亚洲五月色婷婷综合| 一进一出抽搐动态| 久久亚洲精品不卡| 在线观看免费高清a一片| 亚洲av欧美aⅴ国产| 天天躁夜夜躁狠狠躁躁| 欧美另类一区| 久久国产精品大桥未久av| 一边摸一边做爽爽视频免费| 日本vs欧美在线观看视频| 91麻豆精品激情在线观看国产 | 亚洲精华国产精华精| 成人国产一区最新在线观看| 男女国产视频网站| 久久国产精品大桥未久av| 99热全是精品| 一本—道久久a久久精品蜜桃钙片| 亚洲男人天堂网一区| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品久久久久成人av| 一区二区三区精品91| 波多野结衣av一区二区av| 久久人人爽av亚洲精品天堂| 中文字幕另类日韩欧美亚洲嫩草| 国产精品久久久久久精品古装| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产日韩一区二区| 欧美激情久久久久久爽电影 | 搡老熟女国产l中国老女人| 纵有疾风起免费观看全集完整版| 狂野欧美激情性bbbbbb| 亚洲五月婷婷丁香| 中国国产av一级| av天堂在线播放| 99久久精品国产亚洲精品| tube8黄色片| 中文精品一卡2卡3卡4更新| 青春草视频在线免费观看| svipshipincom国产片| 亚洲精品一二三| 青青草视频在线视频观看| 国产三级黄色录像| 中文字幕精品免费在线观看视频| 国产成人精品久久二区二区91| 亚洲欧美精品综合一区二区三区| 欧美日韩成人在线一区二区| 12—13女人毛片做爰片一| 国产成人av教育| av视频免费观看在线观看| 国产亚洲精品一区二区www | 国产男女超爽视频在线观看| 韩国高清视频一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| av有码第一页| 国产视频一区二区在线看| 成年女人毛片免费观看观看9 | 乱人伦中国视频| 两个人看的免费小视频| 亚洲一码二码三码区别大吗| 黑人巨大精品欧美一区二区蜜桃| 亚洲国产精品一区三区| 国产免费福利视频在线观看| 亚洲全国av大片| 搡老乐熟女国产| 亚洲七黄色美女视频| 自线自在国产av| 国产一区二区 视频在线| 久久99一区二区三区| 亚洲av成人不卡在线观看播放网 | 欧美日韩亚洲综合一区二区三区_| 老汉色∧v一级毛片| 国产亚洲午夜精品一区二区久久| 最近中文字幕2019免费版| 999精品在线视频| 夫妻午夜视频| 国内毛片毛片毛片毛片毛片| 午夜视频精品福利| 国产精品 欧美亚洲| 狠狠婷婷综合久久久久久88av| 三上悠亚av全集在线观看| 欧美大码av| 亚洲欧美一区二区三区黑人| 在线永久观看黄色视频| 中文字幕精品免费在线观看视频| 操出白浆在线播放| 国产伦理片在线播放av一区|