999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

微重力模擬裝置中懸吊彈簧局部模態(tài)動力學(xué)建模與試驗研究

2023-02-27 13:14:44楊鴻杰
振動與沖擊 2023年4期
關(guān)鍵詞:模態(tài)振動模型

馬 超,劉 磊,楊鴻杰

(1. 西北工業(yè)大學(xué) 航天學(xué)院,西安 710072; 2. 陜西省空天飛行器設(shè)計重點實驗室,西安 710072)

空間微振動隔離技術(shù)已經(jīng)成功用于空間站以及空間高分辨率遙感衛(wèi)星。微振動隔振器可以解決空間精密載荷對微振動的抑制與隔離問題,從而滿足精密載荷對微重力加速度水平的要求[1]。模擬空間微重力環(huán)境進行地面試驗是驗證隔振器能否在軌可靠工作的重要措施之一。目前進行微重力模擬[2]的主要方法有:水浮法[3]、氣浮法[4-5]、落塔法[6]、拋物線飛行法[7]以及懸吊法[8]等,其中落塔法進行微重力模擬的精度能夠達到10-4g~10-5g,然而受到落塔高度的限制,進行地面試驗的時間很難超過10 s。拋物線飛行法為部分自由落體方式,微重力模擬的精度能夠達到10-2g~10-3g,該方法能夠?qū)⒃囼灂r間提高到20~30 s。然而落塔法和拋物線飛行法進行地面試驗的成本較高,時間較短,使得對于微振動隔振器的各種性能指標(biāo)無法較好的考核,在地面微重力模擬試驗應(yīng)用方面受到很大限制。水浮法利用水的浮力模擬微重力環(huán)境,模擬的時長以及尺寸不受限,然而對于帶電設(shè)備,需要考慮水阻和密封,很大程度限制了地面試驗的應(yīng)用。氣浮法利用氣體的推力實現(xiàn)微重力模擬,實現(xiàn)簡單,力可調(diào),然而時長和載荷質(zhì)量受限,考慮到氣瓶的容量要進行長時間地面試驗只能允許較小的載荷,并且氣浮法會產(chǎn)生較大的噪聲,對微振動隔振器的隔振性能產(chǎn)生一定的影響。懸吊法因其原理簡單,制造方便,微重力模擬時間長等優(yōu)點在微重力模擬技術(shù)中被廣泛應(yīng)用[9-11]。

美國的專家學(xué)者從20世紀(jì)70年代就開始將固定長度懸吊法用于航天員的微重力模擬,NASA[12-13]采用多根固定長度的吊索對航天員身體多個部位進行懸吊,模擬月球微重力下走、跑和跳。劉榮強等[14]使用固定長度的吊索為空間索桿鉸接式伸展臂設(shè)計了一套重力平衡裝置,對伸展臂展開和收攏進行了零重力地面模擬試驗。這種方法構(gòu)型簡單,適用于運動自由度少,運動形式簡單,幾乎無豎直位移的情況。然而吊索長度不可調(diào)節(jié)或不可在線調(diào)節(jié),造成豎直力不可調(diào)節(jié),因此適用范圍小,微振動模擬精度低。卡耐基梅隆大學(xué)的Brown等[15]在1994年提出一種利用吊索配重克服重力來模擬太空低重力環(huán)境的方法,并且可在地面實現(xiàn)三維空間的運動。彭浩等[16]采用吊絲配重懸吊法通過一圈帶有配重的吊絲成功對環(huán)形天線進行重力補償。配重式懸吊法實現(xiàn)簡單,適用于運動自由度少,豎直位移小,豎直加速度小的微低重力模擬,當(dāng)懸吊繩慢速移動時補償效果較好。然而快速移動時配重的慣量嚴(yán)重影響拉力精度。為了解決配重懸吊的問題,彈簧懸吊開始被大量研究,Pyatibratov等[17]利用彈簧懸吊法對零重力和低重力條件下的模擬器進行研制,用于航天員在月球和火星空間的低重力訓(xùn)練。Liu等[18]采用扭桿和彈簧的組合形式對剛度進行調(diào)整,使整個機構(gòu)的剛度較低,并使用魯棒控制進行了高精度低重力懸浮模擬試驗。文獻[19]通過四個懸吊點實現(xiàn)懸掛,每個懸吊點包括兩個呈V形安裝的螺旋彈簧和兩個用于低頻定位的位置傳感器組成,通過彈簧可以實現(xiàn)低頻懸吊。然而以上研究都只關(guān)注了如何獲得剛度低的特性,而沒有考慮彈簧自身模態(tài)振動對載荷的影響。

懸吊彈簧為質(zhì)量不可忽略的長彈簧,描述有質(zhì)量彈簧動力學(xué)特性的模型包括離散模型和連續(xù)模型,離散模型是將懸吊彈簧離散成一系列小的彈簧振子進行分析。通過有限個小的彈簧振子來對整個彈簧的性能進行近似時會引入系統(tǒng)誤差,不能準(zhǔn)確描述懸吊彈簧的動力學(xué)特性。為了提高懸吊彈簧動力學(xué)模型的準(zhǔn)確性,本文采用連續(xù)模型對懸吊彈簧進行動力學(xué)建模,搭建了一種微重力模擬裝置,對所建立的動力學(xué)模型進行試驗驗證。

為了降低隔振器在微重力模擬過程中懸吊彈簧局部模態(tài)振動對隔振器的干擾,通過PVC膠與懸吊彈簧并聯(lián)的方式增加彈簧的模態(tài)阻尼,利用連續(xù)模型對加膠帶懸吊彈簧進行動力學(xué)建模,并開展了被動和主動隔振試驗。試驗結(jié)果顯示PVC膠帶可以減小懸吊彈簧局部模態(tài)處的振動幅值,緩解對隔振器控制性能的限制,更準(zhǔn)確的反應(yīng)在軌工作狀態(tài)。

1 微重力模擬裝置理論分析

使用懸吊方法對隔振器進行微重力模擬的總體方案如圖1所示,微重力模擬裝置由懸吊架、懸吊彈簧、PVC阻尼膠帶、懸吊繩以及輔助機構(gòu)組成。懸吊架布置在試驗臺兩側(cè),由長2 m,寬2 m,高3 m的可拆卸鋁型材搭建,起到支撐的作用;懸吊繩通過滑輪纏繞在絞盤上,用于升降懸吊彈簧;PVC膠帶與懸吊彈簧并聯(lián)后直接與隔振上平臺相連,用于微重力模擬;輔助機構(gòu)包括絞盤和滑輪裝置。隔振下平臺固連在試驗臺上,隔振上平臺與下平臺通過支撐彈簧連接在一起。由于在地面受重力影響隔振上平臺會一直壓縮支撐彈簧,懸吊彈簧進行拉伸時會產(chǎn)生向上的拉力,當(dāng)懸吊彈簧提供的拉力與隔振上平臺的重力平衡,使得支撐彈簧保持在原長的狀態(tài),此時隔振器完成微重力模擬。

圖1 微重力模擬裝置總體方案Fig.1 The overall scheme of microgravity simulation device

為了滿足空間精密載荷對空間微重力環(huán)境的要求,隔振器在1 Hz~100 Hz需要將空間站mg量級的環(huán)境振動降低到μg量級。因此懸吊架在豎直方向設(shè)計基頻高于100 Hz,從而降低懸吊架對隔振器隔振性能的影響。為了達到以上要求,懸吊架采用橫截面積100×100高剛度可拆卸鋁型材作為支撐結(jié)構(gòu)進行搭建,增加了微重力模擬裝置的剛度。在每個豎直梁的底面連接一塊底板,每個豎直梁與橫梁以及與底板之間均用斜梁進行加固,底板采用不銹鋼,保證懸吊架具有較高的穩(wěn)定性和較高豎向的基頻。懸吊架設(shè)計模型如圖2所示。

圖2 懸吊架模型Fig.2 Hanging shelf model

懸吊彈簧為質(zhì)量不可忽略的長彈簧,阻尼較小,在對隔振器進行微重力模擬時,懸吊彈簧在局部模態(tài)處的振動會引入到隔振上平臺,這種模態(tài)振動被加速度傳感器測量后進入主動控制回路,容易導(dǎo)致隔振控制失穩(wěn),對隔振器的隔振性能產(chǎn)生影響,并且這種狀態(tài)與在軌狀態(tài)不符。

懸吊彈簧與隔振器內(nèi)部彈簧為并聯(lián)關(guān)系,懸吊彈簧和隔振器內(nèi)部彈簧的關(guān)系見式(1),懸吊彈簧和隔振器內(nèi)部彈簧的連接示意如圖3所示。

kz=k+kt

(1)

式中:kz為懸吊彈簧和隔振器內(nèi)部彈簧的總剛度;k為懸吊彈簧的剛度;kt為隔振器內(nèi)部彈簧的剛度。

圖3 懸吊彈簧和隔振器內(nèi)部彈簧的連接示意圖Fig.3 Schematic diagram of connection between suspension spring and internal spring of vibration isolator

由式(1)可以看出,懸吊彈簧的剛度越大,總剛度中懸吊彈簧占有的剛度越大,因此懸吊彈簧引入到隔振器的基頻就越大。為了避免懸吊彈簧在1 Hz~100 Hz與隔振器產(chǎn)生共振,懸吊彈簧的基頻應(yīng)低于1 Hz。受到實驗室高度限制懸吊架總高度為3 m,試驗臺高度為0.8 m,因此懸吊彈簧工作狀態(tài)總高度應(yīng)小于2.2 m,設(shè)計懸吊彈簧的剛度k=0.08 N/mm,懸吊彈簧的基頻為0.83 Hz,小于設(shè)計值1 Hz的要求。隔振器內(nèi)部彈簧的剛度kt=7.20 N/mm,因此隔振器的基頻為7.86 Hz,懸吊彈簧與隔振器整體的基頻為7.91 Hz。懸吊彈簧的設(shè)計參數(shù)如表1所示。

表1 懸吊彈簧的設(shè)計參數(shù)Tab.1 Design parameters of suspension spring

2 懸吊彈簧的動力學(xué)建模

微重力模擬裝置主要包括懸吊架、四根懸吊繩以及四根懸吊彈簧。四根完全相同的懸吊彈簧以并聯(lián)的方式連接,在形變量相同的情況下四根懸吊彈簧的剛度可以等效成一根主彈簧的剛度,主彈簧的剛度為四根懸吊彈簧剛度之和。懸吊彈簧連接在懸吊架上可以理解為懸吊彈簧一端固支。將懸吊彈簧看做均勻的連續(xù)介質(zhì),一端固定;另一端系一質(zhì)量為M的彈簧振子系統(tǒng),等效模型如圖4所示。

圖4 懸吊彈簧的等效模型Fig.4 Equivalent model of suspension spring

彈簧的剛度為k,長為L,質(zhì)量為m,阻尼為C1。在t時刻離固定端距離為x點處的彈簧位移量為u(x,t)。觀察其中原長為dx的一小段彈簧,假設(shè)在某一時刻t整個彈簧處于拉伸狀態(tài),則這一小段彈簧的剛度為kL/dx,形變是du,受到的拉力為Fk=kLdu/dx=kLux。這一小段彈簧的阻尼為c1dx/L,彈簧受到的阻尼力為Fc=(c1/L)dx·du/dt,這一小段彈簧的質(zhì)量為mdx/L,加速度為utt,則可以導(dǎo)出彈簧的波動方程為

(2)

為了減小懸吊彈簧對隔振器的影響,通過在懸吊彈簧外部粘貼一層PVC膠帶的方式增加懸吊彈簧阻尼,從而減小彈簧在局部模態(tài)的振動,減小對隔振平臺的影響。

PVC膠帶屬于黏彈性材料,兼具黏性和彈性的性質(zhì)。在外力作用下,其應(yīng)力大小既依賴于應(yīng)變也依賴于應(yīng)變速率,其力學(xué)行為介于理想彈性體(彈簧)和理想黏性體(黏壺)之間。當(dāng)受到外力作用時,類似于彈簧和黏壺共同來應(yīng)對外力的作用,當(dāng)外力消失后,彈簧會恢復(fù)到原來的狀態(tài),而黏壺?zé)o法恢復(fù)到原來的狀態(tài),從而耗散了振動能量,起到阻尼作用。采用彈簧與黏壺并聯(lián)的Kelvin模型對PVC膠帶的黏彈性阻尼進行建模,PVC膠帶與懸吊彈簧通過并聯(lián)方式連接。PVC膠帶與懸吊彈簧的模型如圖5所示。

圖5 PVC膠帶與懸吊彈簧的等效模型Fig.5 Equivalent model of PVC tape and suspension spring

膠帶的剛度為kl,阻尼為C2。觀察其中原長為dx的一小段彈簧和一小段膠帶,設(shè)在某一時刻t整個彈簧和膠帶處于拉伸狀態(tài),則這一小段彈簧受到的拉力為Fk=(k+kl)Ldu/dx=(k+kl)Lux。這一小段膠帶的阻尼為c2dx/L,彈簧受到的阻尼力為Fc=[(c1+c2)/L]dx·du/dt,在式(2)的基礎(chǔ)上彈簧加膠帶的波動方程為

(3)

隔振上平臺和載荷的質(zhì)量為M,懸吊彈簧上端固定,下端與隔振上平臺相連,則邊界條件可表示為

(4)

t=0時懸吊彈簧被隔振上平臺拉伸位移u0,此時懸吊彈簧處于靜止?fàn)顟B(tài)。

則初始條件可表示為

(5)

令u(x,t)=X(x)T(t)代入式(4)和式(5)并分離變量,得

(6)

(7)

(8)

式中:ωn為無阻尼彈簧的頻率;ωd為有阻尼彈簧的頻率;ε為有阻尼彈簧的阻尼比;φ為彈簧有阻尼振動的相位。

(9)

且可以得到

(10)

ytany=β

(11)

解此方程可得到彈簧振子的本征頻率。以yn表示第n個周期所對應(yīng)的本征值y,在|y|<π/2內(nèi),將tany展開成麥克勞林級數(shù),可得

將tany在區(qū)間π/2

(13)

依此方法可求得更高的本征頻率ωn。

解出頻率ωn代入式(8)中可以得到懸吊彈簧的振動解為

e-εωntsin(ωdt+φ)

(15)

隔振上平臺(振子)的振動解為

e-εωntsin(ωdt+φ)

(16)

由式(15)和式(16)可以看出,懸吊彈簧的振動頻率與隔振上平臺的振動頻率相同,說明微重力模擬裝置中懸吊彈簧的振動會引入到隔振上平臺,從而對隔振器產(chǎn)生影響。

隔振上平臺的振動振幅為

(17)

式中,ε為有阻尼彈簧的阻尼比,由式(17)可以看出,阻尼越大隔振上平臺的振動振幅越小,增加阻尼會減小隔振上平臺的振動。

懸吊彈簧為剛性彈簧,阻尼比在0.001~0.005內(nèi),阻尼與阻尼比的關(guān)系為

(18)

式中,ε1為懸吊彈簧的阻尼比。根據(jù)式(18),懸吊彈簧計算阻尼在0.014~0.068內(nèi)。

(19)

式中:E為膠帶的等效彈性模量;B為膠帶的寬度;L為膠帶的長度;τ為流變系數(shù)。

膠帶的性能參數(shù)如表2所示。

表2 PVC膠帶的性能參數(shù)Tab.2 Performance parameters of PVC tape

經(jīng)過計算,膠帶的阻尼為0.510,懸吊彈簧加膠帶后阻尼增加到0.524以上,加膠帶后整體阻尼可以提高8.5倍以上。由此說明加膠帶可以減小隔振上平臺的振動,從而減小微重力模擬裝置中懸吊彈簧對隔振器的影響。

3 微重力模擬裝置的試驗驗證

為了驗證懸吊彈簧局部模態(tài)處的振動對隔振上平臺的影響,本文搭建了地面試驗驗證系統(tǒng)。地面試驗驗證系統(tǒng)包括試驗臺、微重力模擬裝置以及單自由度隔振器。地面試驗驗證系統(tǒng)如圖6所示。

圖6 地面試驗驗證系統(tǒng)Fig.6 The ground experimental system

單自由度隔振器由隔振上平臺和下平臺通過四個支撐彈簧進行連接,隔振下平臺固連在試驗臺上,隔振上平臺與微重力模擬裝置相連。加速度傳感器采用PCB公司生產(chǎn)的393B05加速度傳感器,該傳感器能夠測量0.2 Hz~1 700 Hz頻率范圍內(nèi)的加速度信號,靈敏度為10 V/g,可以采集到隔振上平臺μg量級的低幅微振動。加速度信號調(diào)理盒能夠?qū)⒓铀俣葴y量信號經(jīng)過信號放大、信號調(diào)理以及A/D轉(zhuǎn)換傳遞給dspace實時控制系統(tǒng),并傳輸給主控計算機從而完成低幅微振動的測量。主控計算機可以根據(jù)隔振平臺采集的加速度信號實現(xiàn)主動隔振。

在懸吊狀態(tài)下對隔振器進行被動隔振試驗。隔振上平臺在頻域下的加速度響應(yīng)如圖7所示。

圖7 隔振上平臺振動頻率Fig.7 Vibration frequency of upper platforms

由圖7可以看出,懸吊彈簧實測第一階頻率為0.75 Hz,第二階頻率實測值為5.85 Hz。加速度響應(yīng)的最高峰值為7.6 Hz,該頻率為懸吊彈簧和隔振器整體的實測基頻,懸吊彈簧和隔振器整體基頻理論值為7.91 Hz,理論值與實測值的誤差為3.9%。

通過式(9)、式(12)~式(14)對懸吊彈簧模態(tài)振動頻率進行求解,得到彈簧模態(tài)振動頻率的理論值。將懸吊彈簧模態(tài)振動頻率理論值與實測值(見圖7)進行對比,見表3所示。

優(yōu)點:延續(xù)并兼容PTN技術(shù),依托以太網(wǎng)產(chǎn)業(yè)鏈低成本優(yōu) 勢。基于以太網(wǎng)的FlexE輕量級增強,簡單易實現(xiàn)。

表3 懸吊彈簧振動理論值與實測值的對比Tab.3 Comparison between theoretical and measured vibration values of suspension spring

從表3可以看出,懸吊彈簧的理論基頻為0.83 Hz,與實測第一階頻率誤差為9.6%。懸吊彈簧的第二階頻率理論值5.95 Hz與實測值的誤差為1.21%。懸吊彈簧振動頻率的理論值與實測值除了第一階頻率誤差為9.64%,其他階頻率誤差在4%以內(nèi),說明懸吊彈簧局部模態(tài)處的振動引入到隔振上平臺,驗證了微重力模擬裝置中懸吊彈簧動力學(xué)模型的正確性。

為了減小懸吊彈簧局部模態(tài)處的振動對隔振器性能的影響,使微重力模擬裝置更接近在軌狀態(tài),在懸吊彈簧側(cè)面以粘貼膠帶的方式增加彈簧的阻尼,PVC膠帶與懸吊彈簧的安裝方式如圖8所示。

圖8 PVC膠帶與懸吊彈簧的安裝方式Fig.8 Installation of PVC tape and suspension spring

對四根懸吊彈簧粘貼膠帶前后進行隔振試驗,通過被動隔振后上平臺加速度頻域?qū)Ρ热鐖D9所示。

圖9 懸吊彈簧有無膠帶被動隔振后上平臺加速度對比圖Fig.9 Comparison diagram of acceleration of upper platform after passive vibration isolation of suspension spring with or without tape

由圖9可以看出,懸吊彈簧沒有加膠帶進行被動隔振試驗時,懸吊彈簧在1 Hz~100 Hz的頻段內(nèi)會將自身的局部模態(tài)引入到隔振上平臺,在5.85 Hz,12.16 Hz,18.15 Hz,24.23 Hz,42.35 Hz處的振動均超過隔振要求的1 μg,加上膠帶后隔振上平臺1 Hz~100 Hz頻段內(nèi)由懸吊彈簧引入的振動幅值均小于1 μg,說明在懸吊彈簧外部加膠帶可以增加彈簧的阻尼,并且可以減小懸吊彈簧引入到隔振上平臺的振動,從而減小微重力模擬裝置中懸吊彈簧對隔振器的影響。

為了驗證懸吊彈簧局部模態(tài)處的振動對主動隔振的影響,在微重力模擬狀態(tài)下對隔振器進行主動隔振試驗。主動隔振采用PID控制,當(dāng)增益升高到100時,控制器失穩(wěn),音圈作動器出現(xiàn)嚴(yán)重漂移。通過主動隔振后上下平臺加速度在頻域上的對比如圖10所示。

圖10 上下平臺加速度在頻域下的對比Fig.10 Comparison of acceleration of upper and lower platforms in frequency domain

由圖10可看出,對隔振器進行主動隔振后,由懸吊彈簧局部模態(tài)(5.81 Hz,12.14 Hz,54.52 Hz,60.09 Hz)引入到隔振上平臺的振動無法被隔離到1 μg以下,使得主動隔振無法滿足性能要求。說明懸吊彈簧局部模態(tài)的振動會影響到隔振器的性能。

通過以上試驗可以說明,懸吊彈簧局部模態(tài)的振動已經(jīng)引入到隔振上平臺,這種模態(tài)振動被加速度傳感器測量后進入主動控制回路,導(dǎo)致隔振控制失穩(wěn)。

對四根懸吊彈簧粘貼膠帶后進行主動隔振試驗,主動隔振采用PID控制,增益升高到205時,控制器失穩(wěn),相對于懸吊彈簧無膠帶的微重力模擬裝置,增益提高了一倍。通過主動隔振后上下平臺加速度在頻域上的對比如圖11所示。

圖11 上下平臺加速度在頻域下的對比Fig.11 Comparison of acceleration of upper and lower platforms in frequency domain

由于48.8 Hz為試驗臺的基頻,通過隔振器引入到隔振上平臺,因此與懸吊彈簧沒有關(guān)系。從圖11可以看出,在對隔振器進行主動隔振后,微重力模擬裝置中的懸吊彈簧加膠帶后引入到隔振上平臺的振動均降低到1 μg以下,說明加膠帶可以降低懸吊彈簧局部模態(tài)處的振動,減小對隔振器主動隔振性能的影響。

通過以上試驗可以說明,在懸吊彈簧側(cè)面粘貼PVC膠帶可以起到增加阻尼的效果,降低了懸吊彈簧局部模態(tài)處的振動,從而有效緩解了模態(tài)振動對主動隔振控制器控制性能的限制,這種狀態(tài)更加接近在軌狀態(tài)。

4 結(jié) 論

在利用懸吊法對隔振器進行微重力模擬時,懸吊彈簧局部模態(tài)處的振動很容易引入隔振上平臺,限制隔振控制性能,并且與在軌狀態(tài)不符。本文利用連續(xù)模型對懸吊彈簧和PVC膠帶進行動力學(xué)建模并對純彈簧和加膠帶彈簧的阻尼進行對比,從而對加膠帶懸吊彈簧的模態(tài)振動和阻尼進行深入研究,通過被動和主動隔振試驗對所建立的動力學(xué)模型進行了試驗驗證,研究了加膠帶懸吊彈簧對隔振器隔振性能的影響,從而得到以下結(jié)論:

(1)通過連續(xù)模型進行動力學(xué)建模能夠準(zhǔn)確描述微重力模擬裝置中懸吊彈簧加膠帶的動力學(xué)特性。

(2)懸吊彈簧局部模態(tài)處的振動會引入到隔振上平臺,通過被動和主動隔振都無法隔離,這種狀態(tài)與在軌狀態(tài)不符。采用PVC膠帶可以起到增加阻尼的效果,這種方法能夠?qū)椈傻哪B(tài)振動降低到1 μg以下,從而減小模態(tài)振動對被動隔振的影響,并有效緩解了模態(tài)振動對主動隔振控制器控制性能的限制,這種狀態(tài)能夠更準(zhǔn)確的反應(yīng)在軌工作狀態(tài)。

猜你喜歡
模態(tài)振動模型
一半模型
振動的思考
重要模型『一線三等角』
振動與頻率
重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
中立型Emden-Fowler微分方程的振動性
3D打印中的模型分割與打包
國內(nèi)多模態(tài)教學(xué)研究回顧與展望
基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識別
UF6振動激發(fā)態(tài)分子的振動-振動馳豫
計算物理(2014年2期)2014-03-11 17:01:44
主站蜘蛛池模板: 香蕉久人久人青草青草| 日韩a级片视频| 国产毛片网站| 日韩黄色精品| 亚洲男人的天堂在线| 国产va免费精品| 老司机精品99在线播放| 欧美综合成人| 欧美色亚洲| 中文字幕不卡免费高清视频| 本亚洲精品网站| 精品五夜婷香蕉国产线看观看| 国产精品伦视频观看免费| 一级福利视频| 在线日韩日本国产亚洲| 自拍欧美亚洲| 无码福利日韩神码福利片| 99久久99视频| 亚洲中文无码h在线观看| 四虎永久免费地址| 亚洲人成网站色7777| 国产精品成人第一区| 国产最新无码专区在线| 99精品高清在线播放| 成人午夜视频在线| 国产福利免费视频| 最新国语自产精品视频在| 人妻21p大胆| 亚洲日韩精品无码专区97| 日韩久草视频| 亚洲香蕉久久| 中文字幕 91| 91久久天天躁狠狠躁夜夜| 国产成人午夜福利免费无码r| 毛片网站免费在线观看| 一级毛片在线免费看| 免费一级无码在线网站| 欧美亚洲另类在线观看| 欧美在线中文字幕| 毛片在线看网站| 久久黄色视频影| 国产福利大秀91| 日本高清免费不卡视频| av在线手机播放| 久久久91人妻无码精品蜜桃HD | 秘书高跟黑色丝袜国产91在线| 广东一级毛片| 国产精品吹潮在线观看中文| 国内a级毛片| 高清久久精品亚洲日韩Av| 亚洲永久视频| 成人在线综合| 亚洲欧美一区二区三区图片 | 青青青国产视频| 久久天天躁夜夜躁狠狠| 99热这里只有精品国产99| 综合人妻久久一区二区精品 | 欧美不卡二区| 凹凸国产分类在线观看| lhav亚洲精品| 中文字幕首页系列人妻| 亚洲婷婷六月| 婷婷综合缴情亚洲五月伊| 国产精品亚欧美一区二区三区| 精品五夜婷香蕉国产线看观看| 国产黄在线免费观看| 2020国产免费久久精品99| 91丝袜乱伦| 国产成人调教在线视频| 国产成人高清精品免费软件| 国产美女主播一级成人毛片| 91精品人妻互换| 日韩在线2020专区| 欧美在线视频不卡第一页| 日韩精品成人在线| 久久永久精品免费视频| 狠狠干欧美| 亚洲性一区| 久久永久视频| 熟女日韩精品2区| 在线国产91| 伊人久久大香线蕉aⅴ色|