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

基于自適應(yīng)諧波消除的Hexapod平臺微振動(dòng)激勵(lì)控制

2016-03-24 02:00:26周世驥
航天器環(huán)境工程 2016年1期

周世驥,黃 海

(北京航空航天大學(xué) 宇航學(xué)院,北京 100191)

?

基于自適應(yīng)諧波消除的Hexapod平臺微振動(dòng)激勵(lì)控制

周世驥,黃 海

(北京航空航天大學(xué) 宇航學(xué)院,北京 100191)

摘要:針對空間微振動(dòng)環(huán)境模擬的需求,以Hexapod平臺為對象,進(jìn)行正弦振動(dòng)激勵(lì)控制的研究。當(dāng)Hexapod平臺工作在共振頻段時(shí),其輸出的振動(dòng)信號中因含有諧波成分而產(chǎn)生了顯著的控制誤差。為此,提出了一種自適應(yīng)諧波消除算法。該算法以LMS濾波器為基礎(chǔ),將與諧波同頻率的正弦信號和余弦信號作為濾波器的基底信號,將平臺實(shí)際的輸出響應(yīng)作為濾波器的誤差信號,以此實(shí)現(xiàn)諧波分量的自適應(yīng)消除。將基于該算法的控制回路引入傳統(tǒng)的控制器,進(jìn)行了共振頻段的單輸入單輸出和多輸入多輸出的微振動(dòng)激勵(lì)試驗(yàn),結(jié)果表明,該算法可有效地消除諧波失真,大幅提高了Hexapod平臺在共振頻段的控制精度。

關(guān)鍵詞:振動(dòng)激勵(lì);Hexapod平臺;自適應(yīng)控制;共振;諧波消除

0 引言

隨著星載高精度成像設(shè)備的分辨率不斷提高及其尺寸和焦距的不斷增大,由各種內(nèi)外部因素導(dǎo)致的航天器微振動(dòng)對成像質(zhì)量的影響尤為突出[1]。因此,在模擬空間微振動(dòng)環(huán)境中進(jìn)行航天器有效載荷抗干擾檢定變得不可或缺。但是由于微振動(dòng)的頻率較高且振動(dòng)的加速度量級一般不超過10mg,傳統(tǒng)的低頻、大加速度量級的振動(dòng)臺無法滿足微振動(dòng)激勵(lì)的要求。此外,多軸振動(dòng)試驗(yàn)雖然是更有效的評價(jià)手段,但我國由于缺少系統(tǒng)的試驗(yàn)設(shè)備、試驗(yàn)技術(shù)和方法以及試驗(yàn)規(guī)范,尚未將多軸振動(dòng)試驗(yàn)納入航天器設(shè)計(jì)和驗(yàn)證流程[2]。因此,研制滿足航天應(yīng)用的多軸微振動(dòng)激勵(lì)平臺并研究相應(yīng)的振動(dòng)激勵(lì)控制算法具有重要意義。

Hexapod平臺是一種六自由度并聯(lián)機(jī)構(gòu)[3],具有定位精度高、剛度大、結(jié)構(gòu)穩(wěn)定、承載能力強(qiáng)、運(yùn)動(dòng)慣量小、動(dòng)態(tài)特性好等特點(diǎn),被廣泛應(yīng)用于主動(dòng)振動(dòng)控制領(lǐng)域,如美國海軍研究生院的PPH平臺[4],Honeywell公司的VISS[5]、ULB的主動(dòng)隔振平臺[6]等。同樣,許多振動(dòng)激勵(lì)裝置也采用了Hexapod平臺的構(gòu)型,如美國空軍實(shí)驗(yàn)室的高頻振動(dòng)臺[7]及MTS公司的Model 353振動(dòng)臺[8]等。這類平臺不僅具有多軸運(yùn)動(dòng)的能力,而且其工作頻率和精度都能滿足多軸微振動(dòng)激勵(lì)的要求。

目前在單輸入單輸出(SISO)系統(tǒng)的正弦振動(dòng)激勵(lì)的控制策略中主要采用了基于LMS(Linear Minimum Square)濾波器的自適應(yīng)控制技術(shù)[9-10]。對于多輸入多輸出(MIMO)系統(tǒng)的振動(dòng)激勵(lì)控制,需先對控制對象進(jìn)行解耦,再轉(zhuǎn)化為對多個(gè)SISO系統(tǒng)的自適應(yīng)控制[11]。但是,由于LMS濾波器是一種線性濾波器,當(dāng)共振等因素激發(fā)了系統(tǒng)非線性響應(yīng)并產(chǎn)生諧波失真時(shí),采用單一頻率基底信號的LMS濾波器不能消除該諧波失真,使所產(chǎn)生的控制誤差顯著增加。為此,本文基于LMS濾波器的傳統(tǒng)控制方法,提出一種自適應(yīng)諧波消除(Adaptive Harmonic Cancellation,AHC)算法,并進(jìn)行SISO系統(tǒng)和MIMO系統(tǒng)的振動(dòng)激勵(lì)試驗(yàn),以研究所提出算法的有效性。

1 SISO系統(tǒng)的振動(dòng)激勵(lì)控制方法

1.1基于LMS濾波器的振動(dòng)激勵(lì)控制

在SISO系統(tǒng)中進(jìn)行正弦振動(dòng)激勵(lì)所使用的基于LMS濾波器的控制方法如圖1所示。

圖1 基于LMS濾波器的振動(dòng)激勵(lì)控制Fig. 1 Control of vibration excitation based on LMS filter

在圖1中,Aejφ為控制對象在頻率f0處的傳遞函數(shù)值,為了使系統(tǒng)輸出正弦振動(dòng)d(t)=sin(f0t),則輸入信號應(yīng)滿足

其中:系數(shù)w1、w2為

可見,只需用相同頻率的正弦和余弦信號作為基底信號并乘以合適的系數(shù)后輸入系統(tǒng),便可獲得期望的輸出信號。但由于系統(tǒng)傳遞函數(shù)往往是不可確切獲知的,所以w1、w2不能由式(2)直接確定,而需要根據(jù)LMS算法進(jìn)行迭代求解,即

式中:ke是實(shí)際輸出信號ky與期望信號kd之間的誤差;μ是迭代增益,為確保控制收斂,μ須滿足[4]

當(dāng)系統(tǒng)工作在線性狀態(tài)時(shí),上述控制方法可以達(dá)到較高的控制精度。但在試驗(yàn)中發(fā)現(xiàn),當(dāng)系統(tǒng)在共振等因素的作用下工作在非線性狀態(tài)時(shí),輸出信號中除了含有與基底信號同頻的成分外,還會混疊有頻率為基底信號倍頻的諧波成分r(t),即如圖2所示。

圖2 含有諧波擾動(dòng)的情況Fig. 2 The case with harmonic distortion

在這種情況下,因?yàn)椴捎玫膮⒖夹盘柵c諧波的頻率不同,所以并不能消除諧波的影響,輸出信號始終存在諧波失真。為此,需要在控制回路中增加諧波消除回路來減小由諧波失真產(chǎn)生的控制誤差。

1.2SISO系統(tǒng)的自適應(yīng)諧波消除

對于學(xué)生而言,雙語學(xué)習(xí)取決于自身的自覺性和主動(dòng)性[10]。本次調(diào)查結(jié)果顯示,目前全科住院醫(yī)師英語基礎(chǔ)絕大多數(shù)通過了英語四級,近半數(shù)通過了英語六級,具備簡單交流能力,雖然綜合能力仍有不足,但具備自主學(xué)習(xí)基礎(chǔ)。在醫(yī)學(xué)雙語學(xué)習(xí)的安排上,53.19%的住院醫(yī)師沒有固定的學(xué)習(xí)計(jì)劃,只是想起來時(shí)學(xué)習(xí),17.02%的住院醫(yī)師有明確的學(xué)習(xí)計(jì)劃,但僅3.19%的學(xué)員能夠堅(jiān)持按計(jì)劃完成學(xué)習(xí)任務(wù)。說明多數(shù)住院醫(yī)師雖有學(xué)習(xí)意愿,但學(xué)習(xí)的持續(xù)性和計(jì)劃性不足。

當(dāng)參考信號為正弦信號時(shí),LMS算法可以通過在線迭代將誤差信號中與參考信號同頻率的部分削減到0。考慮到諧波的頻率是可知的,本文利用LMS算法的特點(diǎn),提出一種自適應(yīng)諧波消除算法:即將與諧波同頻率的正弦信號和余弦信號作為LMS濾波器的參考信號,以實(shí)際的輸出響應(yīng)作為LMS濾波器的誤差信號,通過LMS算法在線更新權(quán)系數(shù),最終消除輸出響應(yīng)中的諧波。假定輸出信號僅含有n倍頻諧波,即

則加入諧波消除回路后的SISO系統(tǒng)振動(dòng)激勵(lì)控制框圖如圖3所示。該控制回路包括振動(dòng)信號生成回路和諧波消除回路兩部分。

圖3 含自適應(yīng)諧波消除的微振動(dòng)激勵(lì)控制Fig. 3 Control of micro-vibration excitation with adaptive harmonic cancellation

若系統(tǒng)在頻率nf0處的響應(yīng)為,與式(3)類似,自適應(yīng)諧波消除回路中的權(quán)系數(shù)w3、w4的迭代公式及收斂條件為:

2 MIMO系統(tǒng)的振動(dòng)激勵(lì)控制方法

MIMO系統(tǒng)的振動(dòng)激勵(lì)控制包括前饋解耦和自適應(yīng)控制兩部分[11],如圖4所示。其中,對MIMO系統(tǒng)進(jìn)行解耦需先通過辨識獲取系統(tǒng)的頻響矩陣,并進(jìn)一步根據(jù)其逆矩陣構(gòu)造FIR濾波器進(jìn)行前饋補(bǔ)償。

圖4 MIMO系統(tǒng)的微振動(dòng)激勵(lì)控制Fig. 4 Control of micro-vibration excitation in MIMO system

2.1頻響函數(shù)估計(jì)與FIR濾波器構(gòu)造

系統(tǒng)頻響函數(shù)的辨識可采用H1估計(jì)方法,該方法假設(shè)辨識時(shí)的測量噪聲主要來源于輸出噪聲,見圖5。圖中:H(f)為系統(tǒng)的頻響函數(shù);X(f)為輸入信號的頻譜,在進(jìn)行辨識時(shí),輸入信號常選用高斯白噪聲信號以獲得系統(tǒng)較寬頻帶上的頻率響應(yīng);N(f)為噪聲頻譜;Y(f)為測量信號的頻譜;F(f)是線性系統(tǒng)的輸出譜,其無法直接測得。

圖5 含有輸出噪聲時(shí)的系統(tǒng)模型Fig. 5 System model with output noise

設(shè)共進(jìn)行了k次辨識試驗(yàn),則利用H1方法求解系統(tǒng)估計(jì)頻響函數(shù)的公式為

對于MIMO系統(tǒng),需要分別辨識每個(gè)輸入輸出間的頻響函數(shù)再組集成頻響函數(shù)矩陣,從而利用其逆矩陣進(jìn)行解耦。但是,不能在控制器中直接應(yīng)用,而要通過具有相同頻率響應(yīng)的傳遞函數(shù)實(shí)現(xiàn)。對于指定的頻率,通過構(gòu)造有限沖激濾波器(FIR濾波器)來得到這個(gè)傳遞函數(shù)最為簡便,且得到的傳遞函數(shù)不含極點(diǎn),可以保證控制系統(tǒng)的穩(wěn)定。

2.2MIMO系統(tǒng)的自適應(yīng)諧波消除

圖4所示的控制結(jié)構(gòu)中的自適應(yīng)控制器僅采用了基于1.1節(jié)所述的LMS控制,因此無法消除MIMO系統(tǒng)的諧波失真。為此,本文將自適應(yīng)諧波消除控制引入MIMO系統(tǒng)的振動(dòng)激勵(lì)控制,以消除進(jìn)行多自由度激振時(shí)產(chǎn)生的諧波,如圖6所示。首先通過逆頻響函數(shù)矩陣進(jìn)行前饋解耦,將MIMO系統(tǒng)的控制轉(zhuǎn)化為對多個(gè)SISO系統(tǒng)的控制,然后在每個(gè)SISO控制回路中應(yīng)用圖3的控制方法,從而在產(chǎn)生振動(dòng)信號的同時(shí)消除諧波擾動(dòng)。值得注意的是,諧波消除回路中需要使用與諧波頻率相同的逆頻響函數(shù)矩陣進(jìn)行解耦。

圖6 MIMO系統(tǒng)的自適應(yīng)諧波消除Fig. 6 Adaptive harmonic cancellation in MIMO system

3 試驗(yàn)結(jié)果

為驗(yàn)證所引入的AHC算法的效果,在Hexapod平臺上進(jìn)行了單輸入單輸出和多輸入多輸出的振動(dòng)激勵(lì)控制試驗(yàn)。

3.1試驗(yàn)系統(tǒng)

圖7為所搭建的振動(dòng)激勵(lì)系統(tǒng),包括Hexapod平臺和電控系統(tǒng)。

圖7 基于Hexapod平臺的微振動(dòng)激勵(lì)系統(tǒng)Fig. 7 Micro-vibration generation system based on Hexapod platform

Hexapod平臺主要包括上平臺、下平臺和6個(gè)作動(dòng)器。其中,以作動(dòng)器的控制電壓作為輸入信號,以上平臺的加速度作為輸出信號。作動(dòng)器中的主動(dòng)元件為直線音圈電機(jī)。作動(dòng)器與上平臺通過柔性鉸鏈連接,作動(dòng)器的定子與動(dòng)子之間通過膜片彈簧連接。上平臺上方安裝了3根懸吊彈簧以抵消上平臺的重力。

電控系統(tǒng)由加速度計(jì)、信號調(diào)理器、實(shí)時(shí)控制器、A/D和D/A卡、電源和功率放大器組成。其中,安裝于上平臺的6個(gè)壓電加速度計(jì)可用于采集上平臺各自由度的振動(dòng)加速度并轉(zhuǎn)化為電壓信號;該信號經(jīng)過信號調(diào)理器和A/D卡后進(jìn)入實(shí)時(shí)控制器,然后通過所設(shè)計(jì)的控制方法編寫的控制程序得到控制電壓;控制器輸出的電壓信號分別經(jīng)過D/A卡和功率放大器后驅(qū)動(dòng)作動(dòng)器運(yùn)動(dòng),最終使上平臺產(chǎn)生振動(dòng),如圖8所示。

圖8 電控系統(tǒng)構(gòu)架Fig. 8 Structure of electricity-control system

3.2頻率響應(yīng)辨識

通過H1方法對試驗(yàn)平臺各輸入輸出間的頻率響應(yīng)函數(shù)進(jìn)行了辨識,所得的幅頻曲線如圖9所示。由圖可知,平臺存在多個(gè)共振頻段,激振試驗(yàn)可在這些頻段中選取某個(gè)固定頻點(diǎn)進(jìn)行,本文選取了146Hz。

圖9 平臺的幅頻曲線Fig. 9 Amplitude-frequency curves of the platform

3.3單輸入單輸出激振試驗(yàn)

分別以Hexapod平臺中的1號作動(dòng)器和1號加速度計(jì)作為輸入和輸出進(jìn)行振動(dòng)激勵(lì)試驗(yàn),期望響應(yīng)為幅值5mg、頻率146Hz的正弦振動(dòng)。圖10為期望信號波形與未采用諧波消除時(shí)的實(shí)際響應(yīng)信號波形。可見,輸出信號并不是一個(gè)規(guī)則的正弦信號,并且其幅值超過了期望的5mg,接近10mg。對響應(yīng)信號做傅里葉變換可得到其頻譜,如圖11所示。由圖可知:雖然響應(yīng)信號在146Hz的分量已經(jīng)很好地控制在了期望的5mg,但是其中混疊了2倍頻的諧波成分,這個(gè)292Hz的諧波的幅值已超過5mg,因此實(shí)際響應(yīng)信號的誤差很大。

圖10 無自適應(yīng)諧波消除時(shí)的響應(yīng)波形(SISO)Fig. 10 Response waveform without AHC(SISO)

圖11 無自適應(yīng)諧波消除時(shí)的響應(yīng)頻譜(SISO)Fig. 11 Response spectrum without AHC (SISO)

圖12為引入了AHC控制后的實(shí)際響應(yīng)信號波形與期望信號波形,兩者幾乎重合,說明經(jīng)振動(dòng)控制后,理論值與試驗(yàn)值基本吻合。圖13所示的響應(yīng)信號頻譜也進(jìn)一步表明原有的2倍頻諧波已被有效消除。可見,在單輸入單輸出振動(dòng)激勵(lì)試驗(yàn)中引入AHC控制后,控制精度的提升十分顯著。

圖12 引入自適應(yīng)諧波消除后的響應(yīng)波形(SISO)Fig. 12 Response waveform with AHC(SISO)

圖13 引入自適應(yīng)諧波消除后的響應(yīng)頻譜(SISO)Fig. 13 Response spectrum with AHC(SISO)

3.4多輸入多輸出激振試驗(yàn)

以Hexapod平臺中的6個(gè)作動(dòng)器和6個(gè)加速度計(jì)分別作為輸入和輸出進(jìn)行激振試驗(yàn)。在試驗(yàn)中,令1號加速度計(jì)的期望響應(yīng)為幅值5mg、頻率146Hz的正弦振動(dòng),其余加速度計(jì)處的期望響應(yīng)為0。圖14~圖17分別為加入AHC回路前后的試驗(yàn)結(jié)果。

圖14為未采用自適應(yīng)諧波消除時(shí)各加速度計(jì)所測實(shí)際響應(yīng)信號及期望信號的波形。由圖可知,除3、4號加速度計(jì)所測響應(yīng)的幅值接近期望值外,其余加速度計(jì)處的響應(yīng)幅值都超過了期望值,其中1號加速計(jì)約為20mg。圖15所示為各響應(yīng)信號的頻譜,該圖反映了誤差的主要來源仍為292Hz的2倍頻諧波。

圖14 無自適應(yīng)諧波消除時(shí)的響應(yīng)波形(MIMO)Fig. 14 Response waveforms without AHC(MIMO)

圖15 無自適應(yīng)諧波消除時(shí)的響應(yīng)頻譜(MIMO)Fig. 15 Response spectrums without AHC(MIMO)

圖16為在MIMO系統(tǒng)中引入了AHC控制后各加速度計(jì)所測的實(shí)際響應(yīng)信號和期望信號的波形。其中,1號加速度計(jì)所測的響應(yīng)波形幾乎與期望的5mg的波形重合,且其余加速度計(jì)的響應(yīng)信號僅為噪聲水平,滿足所期望的控制目標(biāo)。

圖16 引入自適應(yīng)諧波消除后的響應(yīng)波形(MIMO)Fig. 16 Response waveforms with AHC(MIMO)

圖17所示的各響應(yīng)信號的頻譜表明,引入了AHC控制后,原有的2倍頻諧波已被消除。可見,在MIMO系統(tǒng)中引入AHC控制后,控制精度明顯提高。

圖17 引入自適應(yīng)諧波消除后的響應(yīng)頻譜(MIMO)Fig. 17 Response spectrums with AHC(MIMO)

4 結(jié)束語

為減小Hexapod平臺在共振頻段進(jìn)行微振動(dòng)激勵(lì)時(shí)由諧波失真產(chǎn)生的控制誤差,分別在SISO系統(tǒng)激振和MIMO系統(tǒng)激振的傳統(tǒng)控制方法中引入了一種自適應(yīng)諧波消除算法。通過試驗(yàn)結(jié)果分析,可得到以下結(jié)論: AHC算法在SISO和MIMO系統(tǒng)的振動(dòng)激勵(lì)中均可以有效地消除輸出響應(yīng)中的諧波成分,從而顯著減少了振動(dòng)誤差,使Hexapod平臺產(chǎn)生的振動(dòng)信號滿足期望要求。

參考文獻(xiàn)(References)

[1] 王紅娟, 王煒, 王欣, 等. 航天器微振動(dòng)對空間相機(jī)像質(zhì)的影響[J]. 光子學(xué)報(bào), 2013, 42(10): 1212 Wang Hongjuan, Wang Wei, Wang Xin, et al. Space camera image degradation induced by satellite micro-vibration[J]. Acta Photonica Sinica, 2013, 42(10): 1212

[2] 次永偉, 邱大蘆, 付樂平, 等. 航天器振動(dòng)試驗(yàn)控制技術(shù)進(jìn)展[J]. 動(dòng)力學(xué)與控制學(xué)報(bào), 2014, 12(3): 193 Ci Yongwei, Qiu Dalu, Fu Leping, et al. Progress in spacecraft vibration testing control technology[J]. Journal of Dynamics and Vibration, 2014, 12(3): 193

[3] Stewart D. A platform with six degrees of freedom[J]. Proceedings of the Institution of Mechanical Engineers, 1965, 180(1): 371-386

[4] Taranti C G. A computationally efficient algorithm for disturbance cancellation to meet the requirements for optical payloads in satellites[D]. CA: Research Office, Naval Postgraduate School, 2001

[5] Foshage G K, Davis T, Sullivan J M, et al. Hybrid active/ passive actuator for spacecraft vibration isolation and suppression[C]//SPIE’s 1996 International Symposium on Optical Science, Engineering, and Instrumentation. Denver, 1996

[6] Hanieh A A. A stiff and soft stewart platforms for active damping and active isolation of vibrations[C]//8thInternational Conference on New Actuators. Bremen, 2002: 254- 257

[7] Ewing C M. The advanced guided weapon tested(AGWT) at the Air Force Research Laboratory Munitions Directorate[C]//AIAA Modeling and Simulation Technologies Conference. Chicago, 2009

[8] MTS Systems Corporation. MAST systems[OL]. [2015-09-30]. http://www.mts.com/ucm/groups/public/documents/library/ dev_ 002251.pdf, 2015

[9] 趙勇. 液壓振動(dòng)臺高精度正弦振動(dòng)的控制策略研究[D].哈爾濱: 哈爾濱工業(yè)大學(xué), 2009

[10] 沈剛. 三自由度電液振動(dòng)臺時(shí)域波形復(fù)現(xiàn)控制策略研究[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2011

[11] Wang H. Multiple-degree-of-freedom sinusoidal vibration generation based on a Hexapod platform[J]. Journal of Systems and Control Engineering, 2015, 229(2): 139-148

(編輯:肖福根)

Control of micro-vibration generation using Hexapod platform based on adaptive harmonic cancellation

Zhou Shiji, Huang Hai
(School of Astronautics, Beihang University, Beijing 100191, China)

Abstract:A Hexapod platform driven by voice coil motors is used to simulate the micro-vibration environment in space. The response signal of the Hexapod platform contains harmonics when it works at resonance frequencies, resulting in a significant error between the actual output and the desired output. To reduce the harmonic distortion and improve the control performance of the Hexapod platform, an adaptive harmonic cancellation (AHC) algorithm based on the LMS filter is developed. The reference signals of the LMS filter are a sine signal and a cosine signal with the same frequency of the harmonics. Meanwhile, the output signal is directly used as the error signal of the LMS filter. Experimental results show that the proposed approach can efficiently eliminate the harmonic distortion and greatly enhance the control accuracy in both SISO and MIMO cases.

Key words:vibration excitation; Hexapod platform; adaptive control; resonance; harmonic cancellation

作者簡介:周世驥(1990—),男,碩士研究生,研究方向?yàn)檎駝?dòng)主動(dòng)控制;E-mail: buaazsj@163.com。指導(dǎo)教師:黃 海(1963—),男,博士學(xué)位,教授,研究方向?yàn)轱w行器結(jié)構(gòu)優(yōu)化,空間智能結(jié)構(gòu)及其控制等;E-mail: hhuang@buaa.edu.cn。

收稿日期:2015-10-24;修回日期:2016-01-08

DOI:10.3969/j.issn.1673-1379.2016.01.016

中圖分類號:TB534

文獻(xiàn)標(biāo)志碼:A

文章編號:1673-1379(2016)01-0089-06

主站蜘蛛池模板: 国产成人免费手机在线观看视频 | av在线手机播放| 久久www视频| 亚洲va精品中文字幕| 免费在线a视频| 巨熟乳波霸若妻中文观看免费| 久久国产精品无码hdav| 国模视频一区二区| 秋霞国产在线| 国产精品无码久久久久久| 免费毛片视频| 激情综合五月网| 国产一区二区三区日韩精品| 亚洲视频色图| 国产无码精品在线| 成人精品免费视频| 国产一线在线| 国产黄在线观看| 四虎影视无码永久免费观看| 欧美精品在线免费| 97色伦色在线综合视频| 国产超碰在线观看| 77777亚洲午夜久久多人| 青青青国产免费线在| 国产乱子伦视频三区| 99久久精品国产综合婷婷| 尤物在线观看乱码| 久久国产乱子| 亚洲国产成人麻豆精品| 538精品在线观看| 91成人免费观看在线观看| 亚洲人人视频| 国模粉嫩小泬视频在线观看| 国产久操视频| 国产91精品久久| 久久精品66| 国产h视频在线观看视频| 日韩小视频网站hq| 高潮毛片无遮挡高清视频播放| 亚洲综合极品香蕉久久网| 欧美成人手机在线观看网址| 国产区网址| 91无码网站| 日韩在线网址| 蜜臀AV在线播放| 欧美国产日韩另类| 成人福利在线视频免费观看| 久久综合亚洲鲁鲁九月天| 香蕉在线视频网站| 亚洲第一视频免费在线| 日韩免费毛片| 丝袜久久剧情精品国产| 国产精品亚欧美一区二区三区| 亚瑟天堂久久一区二区影院| 91精品国产一区| 最新国产高清在线| 草草线在成年免费视频2| 2020精品极品国产色在线观看| 国产精品极品美女自在线网站| 亚洲国产日韩一区| 91色爱欧美精品www| 狠狠v日韩v欧美v| 91麻豆精品国产91久久久久| 久久黄色毛片| 亚洲一区免费看| 久久永久免费人妻精品| 91美女在线| 欧美午夜精品| 国产精品久久久久久久久久98| 国模极品一区二区三区| 国产一级视频久久| 喷潮白浆直流在线播放| 操国产美女| jizz亚洲高清在线观看| 色综合婷婷| 亚洲天天更新| 欧美午夜理伦三级在线观看| 国内精品一区二区在线观看| 欧美日韩中文字幕二区三区| 日本高清在线看免费观看| 亚洲欧美精品在线| 麻豆国产在线不卡一区二区|