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

    基于LSTM神經(jīng)網(wǎng)絡(luò)的壓電執(zhí)行器位移遲滯建模

    2023-05-31 04:51:16時(shí)夢(mèng)想徐希瀟
    壓電與聲光 2023年2期
    關(guān)鍵詞:正弦波執(zhí)行器壓電

    時(shí)夢(mèng)想,胡 泓,吳 浩,徐希瀟

    (1. 哈爾濱工業(yè)大學(xué)(深圳)機(jī)電工程與自動(dòng)化學(xué)院,廣東 深圳 518000;2. 深圳市西渥智控科技有限公司,廣東 深圳 518000)

    0 引言

    LED芯片生產(chǎn)過(guò)程中,在經(jīng)過(guò)前道加工步驟后,晶圓片(Wafer)上芯片的數(shù)目從幾千到數(shù)萬(wàn)不等。為了控制質(zhì)量和降低后續(xù)工序生產(chǎn)成本,需要對(duì)每顆芯片通電后進(jìn)行光電參數(shù)測(cè)試。在半導(dǎo)體測(cè)試領(lǐng)域,生產(chǎn)線上大多采用探針結(jié)構(gòu)通電。目前工業(yè)上探針與芯片電極的接觸和脫離一般通過(guò)傳統(tǒng)步進(jìn)或伺服電機(jī)帶動(dòng)。壓電執(zhí)行器作為一種新型的執(zhí)行機(jī)構(gòu),具有工作精度高,響應(yīng)速度快及驅(qū)動(dòng)力大的優(yōu)點(diǎn),使其適用于各種應(yīng)用場(chǎng)景[1]。本文研究的壓電執(zhí)行器主要應(yīng)用于LED芯片測(cè)試過(guò)程中。使用末端帶有探針的壓電陶瓷執(zhí)行器實(shí)現(xiàn)探針和芯片電極接觸并進(jìn)行測(cè)量。提高芯片測(cè)試效率和探針位移精度。但壓電執(zhí)行器位移與驅(qū)動(dòng)電壓存在非線性關(guān)系,主要源于壓電陶瓷的遲滯效應(yīng)。遲滯效應(yīng)導(dǎo)致壓電執(zhí)行器的輸出位移不僅與當(dāng)前驅(qū)動(dòng)電壓有關(guān),還與之前的驅(qū)動(dòng)電壓有關(guān)。一般遲滯效應(yīng)會(huì)隨著驅(qū)動(dòng)電壓的頻率改變而改變,這給精確控制壓電執(zhí)行器的位移帶來(lái)了挑戰(zhàn)。為了解決這個(gè)問(wèn)題。一方面使用閉環(huán)控制方法,可以達(dá)到較理想的位移控制效果,但閉環(huán)控制需要為壓電執(zhí)行器額外配置位移傳感器,增加了使用成本和系統(tǒng)復(fù)雜度。另一方面使用基于位移遲滯模型的開環(huán)控制方法,但此方法的遲滯模型對(duì)于位移預(yù)測(cè)需要有足夠精度。因此,研究者提出各種模型用于遲滯效應(yīng)的建模。其中數(shù)學(xué)解析模型一般具有大量的參數(shù),識(shí)別過(guò)程較難。靜態(tài)的遲滯模型包括Preisach模型、Prandtl-lshliskii(PI)模型及Maxwell模型等[2]。

    深度學(xué)習(xí)可模擬非線性映射關(guān)系[3],已有很多研究者使用神經(jīng)網(wǎng)絡(luò)對(duì)壓電遲滯效應(yīng)進(jìn)行建模[4]。郭岱宗等[5]使用循環(huán)神經(jīng)網(wǎng)絡(luò)(GRU)對(duì)壓電執(zhí)行器位移遲滯進(jìn)行建模,與傳統(tǒng)的數(shù)學(xué)模型相比,該模型具有較好的模擬遲滯效果。鄒守睿等[4]設(shè)計(jì)了一種循環(huán)神經(jīng)網(wǎng)絡(luò)對(duì)壓電驅(qū)動(dòng)器的遲滯特性進(jìn)行建模,進(jìn)而得到能夠準(zhǔn)確模擬輸出位移和輸入電壓之間關(guān)系的逆模型,并據(jù)此對(duì)壓電驅(qū)動(dòng)器進(jìn)行前饋補(bǔ)償。本文研究了使用長(zhǎng)短期記憶(LSTM)神經(jīng)網(wǎng)絡(luò)模擬遲滯的可行性。

    本研究基于一種杠桿式壓電執(zhí)行器[6],使用3層LSTM神經(jīng)網(wǎng)絡(luò)對(duì)執(zhí)行器位移遲滯進(jìn)行建模,并將所提出的基于LSTM神經(jīng)網(wǎng)絡(luò)遲滯補(bǔ)償模型與傳統(tǒng)Prandtl-Ishlinskii補(bǔ)償模型進(jìn)行對(duì)比,證明了神經(jīng)網(wǎng)絡(luò)模型具有更好的模擬壓電執(zhí)行器壓電位移遲滯的效果。

    1 實(shí)驗(yàn)方法

    1.1 LSTM神經(jīng)網(wǎng)絡(luò)

    壓電遲滯現(xiàn)象主要由于壓電執(zhí)行器當(dāng)前時(shí)刻的位移輸出不僅受當(dāng)前時(shí)刻電壓輸入的影響,還與過(guò)去時(shí)刻的電壓輸入有關(guān)。LSTM是一種特殊的循環(huán)神經(jīng)網(wǎng)絡(luò)(RNN),LSTM的提出主要是為了解決序列預(yù)測(cè)的問(wèn)題,在預(yù)測(cè)當(dāng)前時(shí)間點(diǎn)輸出時(shí)考慮以往時(shí)間點(diǎn)的輸入,這種預(yù)測(cè)方法與遲滯現(xiàn)象發(fā)生的本質(zhì)吻合。

    LSTM由一系列的LSTM細(xì)胞(Cell)組成,每個(gè)細(xì)胞含有若干個(gè)LSTM單元(Unit)。每個(gè)時(shí)間步的細(xì)胞狀態(tài)ct存儲(chǔ)了經(jīng)過(guò)本時(shí)間步及之前的時(shí)間步篩選處理后所留存的信息,從而長(zhǎng)記憶得以建立。

    圖1 LSTM神經(jīng)元

    ft=σ(Wfxt+Ufht-1+bf)

    (1)

    it=σ(Wixt+Uiht-1+bi)

    (2)

    ot=σ(Woxt+Uoht-1+bo)

    (3)

    (4)

    (5)

    ht=ot⊙tanh(ct)

    (6)

    式中:σ(·)為sigmoid激活函數(shù),值域[0,1];tanh(·)為雙曲正切激活函數(shù),值域[-1,1];⊙為矩陣的哈達(dá)瑪積(Hadamard product,按元素相乘);W,U為權(quán)重矩陣;b為偏置向量sigmoid函數(shù)和雙曲正切函數(shù)的定義如下:

    (7)

    (8)

    假設(shè)每次訓(xùn)練的樣本數(shù)(batch size)為b,每個(gè)樣本的維度為d,每個(gè)LSTM Cell隱藏神經(jīng)元數(shù)(hidden size)為h。xt∈d,W,U∈h×d,ht∈h,b∈h。

    1.2 神經(jīng)網(wǎng)絡(luò)超參數(shù)的選擇

    本文所設(shè)計(jì)神經(jīng)網(wǎng)絡(luò)由3層LSTM層和1層全連接層組成,如圖2所示。

    圖2 神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)

    由于壓電執(zhí)行器非線性遲滯由多方面因素決定,故選取3層LSTM網(wǎng)絡(luò)。對(duì)于多層神經(jīng)網(wǎng)絡(luò),上一層每個(gè)Cell的輸入為下面一層對(duì)應(yīng)Cell的隱藏狀態(tài)ht,根據(jù)實(shí)驗(yàn)和訓(xùn)練經(jīng)驗(yàn)選取層數(shù),更多的層數(shù)不僅可能導(dǎo)致較差的擬合效果,還會(huì)增加計(jì)算時(shí)間。

    根據(jù)其他研究者對(duì)壓電遲滯的研究,LSTM窗口大小選取為50個(gè),隱藏單元數(shù)為64個(gè)[8],可以獲得較好的擬合效果及較小的計(jì)算要求。窗口大小表示由ti時(shí)刻及其之前49個(gè)采樣時(shí)刻所組成的激勵(lì)電壓序列ut-49,ut-48,…,ut輸入到神經(jīng)網(wǎng)絡(luò)中預(yù)測(cè)ti時(shí)刻的位移yi。

    由于LSTM層最終輸出為ht維度(h維)的意義與實(shí)際所要預(yù)測(cè)的位移輸出有差別,因此,在最后增加一層全連接層進(jìn)行必要的變換。神經(jīng)網(wǎng)絡(luò)超參數(shù)如表1所示。

    表1 神經(jīng)網(wǎng)絡(luò)超參數(shù)

    1.3 經(jīng)典PI模型

    經(jīng)典PI遲滯模型是一種數(shù)學(xué)解析模型,該模型采用雙邊Play算子Frj(t)和權(quán)重系數(shù)pj加權(quán)疊加來(lái)描述壓電陶瓷遲滯非線性特性[2]。經(jīng)典PI模型在離散處理后的表達(dá)式為

    (9)

    式中:rj為離散閾值(r0=0);pj為每個(gè)遲滯算子對(duì)應(yīng)的權(quán)重系數(shù),兩者通過(guò)遲滯實(shí)驗(yàn)數(shù)據(jù)識(shí)別得出。

    對(duì)每個(gè)離散的rj值,相應(yīng)的遲滯算子具體計(jì)算方法為

    (10)

    式中z0j為各個(gè)遲滯算子初值。

    為獲得模型,繼續(xù)對(duì)時(shí)間進(jìn)行離散,每個(gè)rj值對(duì)應(yīng)的遲滯算子計(jì)算方法為

    (11)

    式中:k=0,1,…,Ns;t0=0;tE=NsTs;Ts為采樣間隔,共使用Ns+1個(gè)采用數(shù)據(jù)。

    將離散后的遲滯算子和權(quán)重系數(shù)疊加,計(jì)算出壓電執(zhí)行器的位移輸出為

    (12)

    2 實(shí)驗(yàn)系統(tǒng)

    為了收集數(shù)據(jù)進(jìn)行LSTM神經(jīng)網(wǎng)絡(luò)的訓(xùn)練和測(cè)試,整體的實(shí)驗(yàn)框圖如圖3所示。

    圖3 實(shí)驗(yàn)系統(tǒng)

    圖3包含一個(gè)用于研究位移遲滯的壓電執(zhí)行器,其功能原理如圖4所示。該執(zhí)行器由杠桿機(jī)構(gòu)將壓電疊堆的唯一輸出放大輸出到末端,該壓電疊堆型號(hào)為P-885.91,最大工作電壓為110 V,壓電執(zhí)行器組裝時(shí),通過(guò)彈簧對(duì)壓電疊堆施加一定的預(yù)緊力,以保證其穩(wěn)定運(yùn)行。

    圖4 壓電執(zhí)行器功能原理圖

    上位機(jī)使用Labview軟件結(jié)合數(shù)據(jù)采集卡廠商的驅(qū)動(dòng)程序,編寫數(shù)據(jù)采集和電壓控制信號(hào)發(fā)出的程序。電壓控制信號(hào)通過(guò)與上位機(jī)相連的USB-5621A數(shù)據(jù)采集卡發(fā)出,輸出采樣頻率為10 kHz。

    ATG-2042型號(hào)的功率放大器將數(shù)據(jù)采集卡輸出的控制信號(hào)電壓進(jìn)行放大,輸出給壓電執(zhí)行器,驅(qū)動(dòng)其產(chǎn)生位移。

    位移傳感器探測(cè)壓電執(zhí)行器的位移量,傳感器選取型號(hào)為米朗LVDT8筆式傳感器,該位移傳感器精度為2.5 μm。最后USB-5621A將傳感器位移信號(hào)的采集輸入到上位機(jī)保存,位移信號(hào)和電壓控制信號(hào)同步采集,采樣頻率均為10 kHz,最終獲取壓電執(zhí)行器的位移遲滯信息。

    3 模型辨識(shí)

    3.1 模型識(shí)別訓(xùn)練數(shù)據(jù)獲取

    為了提供訓(xùn)練LSTM神經(jīng)網(wǎng)絡(luò)所需的多種遲滯數(shù)據(jù),一方面使用多頻率的正弦信號(hào)作為激勵(lì)電壓,另一方面使用多頻率的衰減正弦信號(hào)作為激勵(lì)電壓。得到壓電執(zhí)行器的多個(gè)遲滯環(huán)數(shù)據(jù),即:

    (13)

    (14)

    式中:τ=0.15為衰減正弦波的時(shí)間常數(shù);f為正弦波的頻率;Vmax,Vmin分別為最大、最小激勵(lì)電壓。

    由于功率放大器只能輸出正電壓,為防止之后的補(bǔ)償電壓為負(fù)值,Vmin=30 V,由于壓電疊堆限制,Vmax=110 V。兩種波形在1~60 Hz內(nèi)各采樣60組數(shù)據(jù),傳感器采樣頻率為10 kHz,每種信號(hào)持續(xù)20 s,共120組電壓-位移數(shù)據(jù),每種信號(hào)進(jìn)行降采樣后共40 000個(gè)采樣點(diǎn)用于訓(xùn)練和驗(yàn)證神經(jīng)網(wǎng)絡(luò)模型。通過(guò)兩種激勵(lì)信號(hào)產(chǎn)生了多種遲滯環(huán)。圖5為正弦波激勵(lì)信號(hào)所產(chǎn)生的遲滯環(huán)。圖6為衰減正弦波產(chǎn)生的遲滯環(huán)。

    圖5 正弦波激勵(lì)信號(hào)產(chǎn)生的遲滯環(huán)

    圖6 衰減正弦波產(chǎn)生的遲滯環(huán)

    由圖5、6可看出,壓電執(zhí)行器的遲滯環(huán)與頻率相關(guān),驅(qū)動(dòng)頻率越高,遲滯非線性越明顯,且是非對(duì)稱的。這是由于壓電執(zhí)行器的非線性特性不僅由壓電疊堆本身決定,還與機(jī)械零件之間的摩擦(如在杠桿與支點(diǎn)間的摩擦,頂針和頂針套間的摩擦)及零件的彈性形變有關(guān)。

    3.2 LSTM訓(xùn)練過(guò)程

    訓(xùn)練過(guò)程中,為提高學(xué)習(xí)效率,所有訓(xùn)練數(shù)據(jù)被歸一化到[-1,1]區(qū)間內(nèi),與神經(jīng)網(wǎng)絡(luò)中的激活函數(shù)相匹配。使用Pytorch開源機(jī)器學(xué)習(xí)框架實(shí)現(xiàn)神經(jīng)網(wǎng)絡(luò)。訓(xùn)練迭代次數(shù)選取為80。本研究選用均方損失函數(shù)(MSE),并選擇Adam算法作為反向傳播優(yōu)化算法,學(xué)習(xí)率設(shè)定為0.001。使用前述的神經(jīng)網(wǎng)絡(luò)及超參數(shù)訓(xùn)練神經(jīng)網(wǎng)絡(luò)模型。將訓(xùn)練集和驗(yàn)證集按7∶3劃分。本文選取正弦波和衰減正弦波的頻率為2 Hz、10 Hz、20 Hz、30 Hz、50 Hz,以此作為模型訓(xùn)練集和驗(yàn)證集數(shù)據(jù),其他數(shù)據(jù)作為測(cè)試集數(shù)據(jù)。

    訓(xùn)練使用計(jì)算機(jī)配置:中央處理器為Intel Xeon E5-2680 v4 CPU,主頻2.4 GHz;圖形處理器(GPU)為:NVIDIA Tesla P40(12TFLOPS 單精度浮點(diǎn)計(jì)算,47INT8 TOPS);內(nèi)存為:DDR4,內(nèi)存帶寬達(dá)2 666 MT/s。網(wǎng)絡(luò)需要訓(xùn)練的參數(shù)共有83 777個(gè)。共完成訓(xùn)練所需時(shí)間為30 min。訓(xùn)練集損失函數(shù)和驗(yàn)證集準(zhǔn)確率值隨迭代次數(shù)的變化曲線如圖7所示。

    圖7 損失曲線

    3.3 PI模型辨識(shí)

    為驗(yàn)證LSTM神經(jīng)網(wǎng)絡(luò)的有效性,以傳統(tǒng)PI模型作為對(duì)照組。由于經(jīng)典PI模型為靜態(tài)模型,當(dāng)研究分析遲滯靜態(tài)特性時(shí),輸入信號(hào)頻率不高,頻率過(guò)高將導(dǎo)致模型誤差過(guò)大。使用5 Hz正弦波遲滯數(shù)據(jù)進(jìn)行模型參數(shù)識(shí)別。

    識(shí)別PI模型首先確定遲滯算子的數(shù)量。選用遲滯算子數(shù)量越多時(shí),模型擬合精度越高,計(jì)算量越大。本文選擇遲滯算子個(gè)數(shù)為9個(gè),閾值和遲滯算子初始值[2]為

    (15)

    式中電壓最大值max{|v(k·Ts)|}為實(shí)際輸入電壓最大值110 V減去基值電壓30 V,即80 V。

    遲滯算子的初值均設(shè)置為0,即:

    z0j=0

    (16)

    使用Matlab編寫PI遲滯模型的識(shí)別和驗(yàn)證程序。權(quán)重系數(shù)計(jì)算使用最小二乘法,即通過(guò)最小化PI模型計(jì)算位移和實(shí)際位移間的誤差來(lái)辨識(shí)9個(gè)權(quán)重系數(shù)q值。模型辨識(shí)結(jié)果如表2所示。

    表2 PI模型辨識(shí)結(jié)果

    4 實(shí)驗(yàn)結(jié)果分析

    本文使用3種準(zhǔn)確性的指標(biāo)來(lái)對(duì)模型預(yù)測(cè)結(jié)果進(jìn)行評(píng)估,分別是均方根誤差(ERMS)、ERMS在觀測(cè)值范圍內(nèi)的位置(ENRMS)和最大絕對(duì)誤差(EMA)。

    EMA表示預(yù)測(cè)值和觀測(cè)值之間絕對(duì)誤差的最大值。計(jì)算方法如下:

    (i=1,2,…,N)

    (17)

    ERMS表示預(yù)測(cè)值和觀測(cè)值之間差異(稱為殘差)的樣本標(biāo)準(zhǔn)差。計(jì)算方法如下:

    (i=1,2,…,N)

    (18)

    ENRMS表示ERMS在觀測(cè)值范圍內(nèi)的位置:

    (19)

    由實(shí)驗(yàn)結(jié)果可知:

    1) 由PI模型的預(yù)測(cè)誤差結(jié)果(見表3)可知,經(jīng)典PI模型對(duì)于該執(zhí)行器的遲滯仍具有較好的靜態(tài)建模能力。雖然最大絕對(duì)誤差(EMA)值仍有6.50 μm,但是誤差保證在[-5,+5]內(nèi),誤差對(duì)稱,ENRMS僅為2.15%。圖8為對(duì)用于模型識(shí)別的5 Hz正弦波所對(duì)應(yīng)的位移遲滯預(yù)測(cè)結(jié)果。

    表3 PI模型的位移預(yù)測(cè)誤差

    圖8 PI模型5 Hz正弦波位移預(yù)測(cè)

    2) PI模型動(dòng)態(tài)預(yù)測(cè)能力很差,對(duì)于9 Hz信號(hào)ENRMS已達(dá)2.34%。對(duì)于1 Hz信號(hào)預(yù)測(cè)(見圖9)ENRMS誤差已超過(guò)3%,很難達(dá)到工業(yè)上的精度要求。對(duì)于同頻率衰減正弦波預(yù)測(cè)誤差更大,EMA已達(dá)10.95 μm,模型預(yù)測(cè)基本失效。

    圖9 PI模型1 Hz正弦波位移預(yù)測(cè)

    將上述的經(jīng)典PI模型作為對(duì)照組,由LSTM正弦波和衰減正弦波實(shí)驗(yàn)預(yù)測(cè)誤差(見表4、5)可知,LSTM神經(jīng)網(wǎng)絡(luò)模型具有廣泛的動(dòng)態(tài)預(yù)測(cè)能力。

    表4 LSTM神經(jīng)網(wǎng)絡(luò)正弦波預(yù)測(cè)誤差

    表5 LSTM神經(jīng)網(wǎng)絡(luò)衰減正弦波預(yù)測(cè)誤差

    3) 對(duì)于多頻率正弦波信號(hào),神經(jīng)網(wǎng)絡(luò)的最大絕對(duì)誤差可以保證在±5 μm內(nèi)基本上可滿足LED測(cè)試要求。ENRMS可控制在2%內(nèi),精度較高。LSTM模型對(duì)5 Hz正弦波激勵(lì)位移預(yù)測(cè)如圖10所示。

    圖10 LSTM模型5 Hz正弦波位移預(yù)測(cè)

    4) 對(duì)于衰減系數(shù)τ=0.15的各頻率衰減正弦波遲滯,ENRMS基本保證在3%內(nèi),說(shuō)明神經(jīng)網(wǎng)絡(luò)不僅對(duì)正弦波具有動(dòng)態(tài)建模能力,對(duì)于其他波形同樣具有良好的預(yù)測(cè)能力。LSTM模型對(duì)5 Hz衰減正弦波激勵(lì)位移預(yù)測(cè)如圖11所示。

    圖11 LSTM模型5 Hz衰減正弦波位移預(yù)測(cè)

    5 結(jié)束語(yǔ)

    本文為了提高壓電執(zhí)行器的定位精度,分別使用經(jīng)典PI模型和LSTM神經(jīng)網(wǎng)絡(luò)對(duì)其靜態(tài)和動(dòng)態(tài)遲滯進(jìn)行建模。使用數(shù)據(jù)采集卡同時(shí)采集壓電執(zhí)行器的驅(qū)動(dòng)電壓控制信號(hào)和位移信號(hào),從而獲取壓電執(zhí)行器的多頻率正弦波和衰減正弦波的遲滯數(shù)據(jù)。設(shè)計(jì)了一種由3層LSTM層和1層全連接層構(gòu)成LSTM神經(jīng)網(wǎng)絡(luò)遲滯模型,使用GPU加速,通過(guò)80輪迭代,通過(guò)已獲取的遲滯數(shù)據(jù)計(jì)算神經(jīng)網(wǎng)絡(luò)參數(shù),并對(duì)經(jīng)典PI模型參數(shù)進(jìn)行識(shí)別。最后分別考察兩種遲滯模型的位移預(yù)測(cè)效果,發(fā)現(xiàn)傳統(tǒng)PI模型對(duì)于靜態(tài)遲滯有較好的建模效果,動(dòng)態(tài)效果較差,對(duì)于衰減正弦波的建?;臼?。以PI模型作對(duì)照組,LSTM模型不僅在動(dòng)態(tài)遲滯建模方面具有優(yōu)良的性能,整體誤差保證在2%內(nèi),而且對(duì)于衰減正弦波同樣可獲得較良好的建模效果,基本上可達(dá)到工業(yè)應(yīng)用要求。通過(guò)對(duì)壓電執(zhí)行器的遲滯效應(yīng)進(jìn)行有效建模,可以為L(zhǎng)ED晶圓檢測(cè)提供一種更加高效的方法。

    猜你喜歡
    正弦波執(zhí)行器壓電
    單相正弦波變頻電源設(shè)計(jì)與實(shí)現(xiàn)
    采用BC5016S的純正弦波逆變器設(shè)計(jì)及制作
    電子制作(2019年7期)2019-04-25 13:18:18
    《壓電與聲光》征稿啟事
    壓電與聲光(2019年1期)2019-02-22 09:46:06
    測(cè)控技術(shù)(2018年12期)2018-11-25 09:37:50
    新型壓電疊堆泵設(shè)計(jì)及仿真
    飛機(jī)裝配預(yù)連接緊固件自動(dòng)化安裝末端執(zhí)行器設(shè)計(jì)
    考慮執(zhí)行器飽和的改進(jìn)無(wú)模型自適應(yīng)控制
    一類具有執(zhí)行器飽和的非線性系統(tǒng)抗飽和方法研究
    基于嵌入式技術(shù)的電網(wǎng)同步正弦波發(fā)生裝置研究
    基于壓電激振的彈性模量測(cè)量方法
    国产免费视频播放在线视频| 在线观看免费日韩欧美大片 | 久久久a久久爽久久v久久| 特大巨黑吊av在线直播| 女的被弄到高潮叫床怎么办| 久热久热在线精品观看| 麻豆乱淫一区二区| 国内精品宾馆在线| 大香蕉97超碰在线| 国产老妇伦熟女老妇高清| 亚洲精品日本国产第一区| 亚洲欧美清纯卡通| 婷婷色综合www| 国产淫片久久久久久久久| 成年人午夜在线观看视频| 简卡轻食公司| 热99国产精品久久久久久7| 五月天丁香电影| 亚洲精品第二区| 日韩一区二区三区影片| 欧美丝袜亚洲另类| 国产中年淑女户外野战色| 亚洲精品成人av观看孕妇| 三级国产精品片| 中文字幕亚洲精品专区| 高清毛片免费看| 欧美 日韩 精品 国产| 涩涩av久久男人的天堂| 一个人免费看片子| 国产伦精品一区二区三区四那| 国产伦精品一区二区三区四那| 欧美老熟妇乱子伦牲交| 日韩一区二区视频免费看| 久久久欧美国产精品| 亚洲最大成人中文| av天堂中文字幕网| 寂寞人妻少妇视频99o| 欧美日韩一区二区视频在线观看视频在线| 欧美区成人在线视频| 日本黄大片高清| av在线观看视频网站免费| 哪个播放器可以免费观看大片| av免费观看日本| 日日撸夜夜添| 国产精品免费大片| 在线观看美女被高潮喷水网站| 人人妻人人添人人爽欧美一区卜 | 狂野欧美激情性bbbbbb| 国产无遮挡羞羞视频在线观看| 久久久久久人妻| 我要看黄色一级片免费的| 在线 av 中文字幕| 啦啦啦视频在线资源免费观看| 男男h啪啪无遮挡| 久久久色成人| 国产视频内射| 男人狂女人下面高潮的视频| 美女高潮的动态| 九九爱精品视频在线观看| 老熟女久久久| 国产精品麻豆人妻色哟哟久久| 联通29元200g的流量卡| 91精品一卡2卡3卡4卡| 免费大片黄手机在线观看| 国产高潮美女av| 亚洲精品aⅴ在线观看| 日本av免费视频播放| 九九久久精品国产亚洲av麻豆| 大码成人一级视频| 深爱激情五月婷婷| av国产久精品久网站免费入址| 亚洲国产精品成人久久小说| 如何舔出高潮| 久久久久久久国产电影| 亚洲精品国产av蜜桃| 久久久a久久爽久久v久久| 国产色爽女视频免费观看| 国产精品人妻久久久久久| 三级经典国产精品| 纵有疾风起免费观看全集完整版| 爱豆传媒免费全集在线观看| 亚洲一级一片aⅴ在线观看| 少妇的逼好多水| 免费大片18禁| 精品久久久精品久久久| 日韩亚洲欧美综合| 亚洲av免费高清在线观看| 久久人人爽人人片av| 国产精品蜜桃在线观看| 国产淫片久久久久久久久| 插阴视频在线观看视频| 亚洲欧美中文字幕日韩二区| 啦啦啦啦在线视频资源| 亚洲av不卡在线观看| 亚洲图色成人| 日本欧美视频一区| 久久久久久久久久成人| 国产精品久久久久久久电影| 成人毛片60女人毛片免费| 久久99精品国语久久久| 色婷婷av一区二区三区视频| 18禁裸乳无遮挡动漫免费视频| 成人国产av品久久久| 丝袜脚勾引网站| 一区二区三区免费毛片| 成年人午夜在线观看视频| 中文字幕亚洲精品专区| av在线老鸭窝| 免费黄频网站在线观看国产| 国产又色又爽无遮挡免| 亚洲中文av在线| 久久久久久九九精品二区国产| 国产精品偷伦视频观看了| 视频区图区小说| 日韩av免费高清视频| 视频中文字幕在线观看| 日本一二三区视频观看| 亚洲人成网站在线播| 久久精品国产鲁丝片午夜精品| 2018国产大陆天天弄谢| 国产亚洲午夜精品一区二区久久| 777米奇影视久久| 亚洲精品,欧美精品| 久久久久久久亚洲中文字幕| 啦啦啦中文免费视频观看日本| 国产成人免费无遮挡视频| 夜夜骑夜夜射夜夜干| 久久精品久久久久久久性| 中文在线观看免费www的网站| 中文字幕久久专区| 欧美人与善性xxx| 国产伦在线观看视频一区| 建设人人有责人人尽责人人享有的 | 国产成人精品久久久久久| 日韩国内少妇激情av| 国产av码专区亚洲av| 一区二区三区精品91| 日韩不卡一区二区三区视频在线| 麻豆国产97在线/欧美| 欧美日韩视频高清一区二区三区二| av国产久精品久网站免费入址| 亚洲人与动物交配视频| av在线app专区| 成年女人在线观看亚洲视频| 日韩国内少妇激情av| 日本黄大片高清| 亚洲av成人精品一二三区| 蜜臀久久99精品久久宅男| 亚洲成人中文字幕在线播放| 国产精品一区二区性色av| 日产精品乱码卡一卡2卡三| h日本视频在线播放| 日本免费在线观看一区| 99热全是精品| 国产精品国产三级国产av玫瑰| 亚洲精品,欧美精品| 亚洲精品,欧美精品| 黑人高潮一二区| 天天躁日日操中文字幕| 国产无遮挡羞羞视频在线观看| 成人午夜精彩视频在线观看| 国产高清不卡午夜福利| 中文欧美无线码| 国产精品欧美亚洲77777| 高清欧美精品videossex| 热re99久久精品国产66热6| 男人添女人高潮全过程视频| 成人免费观看视频高清| av播播在线观看一区| 国产亚洲午夜精品一区二区久久| 美女cb高潮喷水在线观看| 免费黄网站久久成人精品| 亚洲在久久综合| freevideosex欧美| 黑人高潮一二区| 成人二区视频| 午夜免费男女啪啪视频观看| 欧美xxⅹ黑人| 少妇人妻久久综合中文| 黄色怎么调成土黄色| 下体分泌物呈黄色| 日日啪夜夜爽| 亚洲国产日韩一区二区| 一区二区三区精品91| 亚洲av男天堂| 丰满乱子伦码专区| 日韩亚洲欧美综合| 国产精品麻豆人妻色哟哟久久| 久久久久久久久久久免费av| 观看av在线不卡| 亚洲国产精品一区三区| 秋霞在线观看毛片| 精品久久久噜噜| 你懂的网址亚洲精品在线观看| 美女国产视频在线观看| 国产黄片视频在线免费观看| 色视频www国产| 国产av一区二区精品久久 | 国产一区二区三区av在线| 国国产精品蜜臀av免费| 高清欧美精品videossex| 国产在线一区二区三区精| 国产精品一区二区三区四区免费观看| 日本色播在线视频| 美女xxoo啪啪120秒动态图| 在线 av 中文字幕| 一个人免费看片子| 夜夜看夜夜爽夜夜摸| 日本wwww免费看| 久久韩国三级中文字幕| 国产精品久久久久久久久免| 人人妻人人爽人人添夜夜欢视频 | 久久久久久久久久久免费av| 观看av在线不卡| 久久精品久久精品一区二区三区| 国产无遮挡羞羞视频在线观看| 尾随美女入室| 汤姆久久久久久久影院中文字幕| 欧美少妇被猛烈插入视频| 丝袜脚勾引网站| av黄色大香蕉| 欧美激情极品国产一区二区三区 | 99九九线精品视频在线观看视频| 国产黄色视频一区二区在线观看| 久热这里只有精品99| 男女免费视频国产| 中文精品一卡2卡3卡4更新| 午夜福利影视在线免费观看| 久久国产精品男人的天堂亚洲 | 日本vs欧美在线观看视频 | av不卡在线播放| 久久久久久久亚洲中文字幕| 国产熟女欧美一区二区| 少妇人妻 视频| 国产黄片视频在线免费观看| 啦啦啦啦在线视频资源| 亚洲精华国产精华液的使用体验| 在线观看av片永久免费下载| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲国产色片| 在线观看一区二区三区激情| 三级国产精品欧美在线观看| 久久久久视频综合| 国产精品嫩草影院av在线观看| 极品少妇高潮喷水抽搐| 国产成人精品久久久久久| 男女国产视频网站| 国产人妻一区二区三区在| 成年女人在线观看亚洲视频| 纵有疾风起免费观看全集完整版| 国产成人精品婷婷| 亚洲国产色片| 秋霞伦理黄片| 久久久亚洲精品成人影院| 免费av中文字幕在线| 少妇人妻精品综合一区二区| 精品视频人人做人人爽| 男女免费视频国产| 自拍欧美九色日韩亚洲蝌蚪91 | 99热6这里只有精品| 赤兔流量卡办理| 黄片无遮挡物在线观看| 视频中文字幕在线观看| 欧美成人精品欧美一级黄| 亚洲久久久国产精品| 欧美日韩精品成人综合77777| 一二三四中文在线观看免费高清| 国产爽快片一区二区三区| 亚洲成人av在线免费| 亚洲av.av天堂| 色婷婷久久久亚洲欧美| 亚洲丝袜综合中文字幕| 亚洲精品自拍成人| 国产精品女同一区二区软件| 国产精品久久久久久精品古装| 久久鲁丝午夜福利片| 亚洲精品乱久久久久久| 免费黄色在线免费观看| 日韩欧美一区视频在线观看 | 欧美bdsm另类| 精品久久久精品久久久| 夜夜爽夜夜爽视频| 男女边摸边吃奶| 一二三四中文在线观看免费高清| av在线app专区| 我的老师免费观看完整版| 有码 亚洲区| 久久精品夜色国产| 日本黄色日本黄色录像| 亚洲伊人久久精品综合| a 毛片基地| a级一级毛片免费在线观看| 久久热精品热| 美女xxoo啪啪120秒动态图| 日韩亚洲欧美综合| 成年女人在线观看亚洲视频| 只有这里有精品99| 天天躁日日操中文字幕| av在线老鸭窝| 日韩人妻高清精品专区| 欧美少妇被猛烈插入视频| 亚洲精品中文字幕在线视频 | 久久亚洲国产成人精品v| 欧美日韩在线观看h| 亚洲精品国产av成人精品| 国产又色又爽无遮挡免| 少妇丰满av| 亚洲精品久久午夜乱码| 亚洲欧洲日产国产| 国产亚洲午夜精品一区二区久久| 建设人人有责人人尽责人人享有的 | 观看免费一级毛片| 纵有疾风起免费观看全集完整版| 我要看日韩黄色一级片| 熟女人妻精品中文字幕| 性高湖久久久久久久久免费观看| 丰满人妻一区二区三区视频av| 欧美日韩视频高清一区二区三区二| 久久国产精品男人的天堂亚洲 | 这个男人来自地球电影免费观看 | 观看免费一级毛片| 成人毛片a级毛片在线播放| 丰满迷人的少妇在线观看| 久久久久久久大尺度免费视频| 久久国产精品大桥未久av | 国产国拍精品亚洲av在线观看| 久久久久性生活片| 久久久久国产网址| 亚州av有码| 狂野欧美白嫩少妇大欣赏| 免费大片黄手机在线观看| 日韩中文字幕视频在线看片 | 视频区图区小说| 国产精品一区www在线观看| 三级经典国产精品| 亚洲成人手机| 亚洲成人av在线免费| av卡一久久| 性高湖久久久久久久久免费观看| 美女福利国产在线 | 青春草亚洲视频在线观看| 国产在线男女| 一级毛片电影观看| 在线播放无遮挡| 国产亚洲最大av| a级一级毛片免费在线观看| 国产精品免费大片| 久久毛片免费看一区二区三区| 色哟哟·www| 青春草视频在线免费观看| xxx大片免费视频| 成人美女网站在线观看视频| 亚洲av福利一区| 又黄又爽又刺激的免费视频.| 菩萨蛮人人尽说江南好唐韦庄| 内地一区二区视频在线| 少妇人妻 视频| 成年av动漫网址| 久久久久精品久久久久真实原创| 亚洲av在线观看美女高潮| 九九久久精品国产亚洲av麻豆| 亚洲美女搞黄在线观看| 日韩制服骚丝袜av| 少妇精品久久久久久久| 亚洲自偷自拍三级| 久久婷婷青草| 九色成人免费人妻av| 男女无遮挡免费网站观看| 自拍欧美九色日韩亚洲蝌蚪91 | 纵有疾风起免费观看全集完整版| 国产成人aa在线观看| 午夜福利在线观看免费完整高清在| 99久久精品国产国产毛片| 亚洲国产高清在线一区二区三| 亚洲欧美清纯卡通| 欧美日韩一区二区视频在线观看视频在线| 大香蕉97超碰在线| 国产精品一区二区性色av| 中文乱码字字幕精品一区二区三区| 男女免费视频国产| 最近的中文字幕免费完整| 亚洲激情五月婷婷啪啪| 国产69精品久久久久777片| 亚洲精品久久久久久婷婷小说| 久久久久性生活片| 免费观看a级毛片全部| 亚洲av不卡在线观看| 男女下面进入的视频免费午夜| 纯流量卡能插随身wifi吗| 久久久久精品久久久久真实原创| 中文字幕人妻熟人妻熟丝袜美| 七月丁香在线播放| 联通29元200g的流量卡| 噜噜噜噜噜久久久久久91| 91精品一卡2卡3卡4卡| 欧美一区二区亚洲| 色综合色国产| 最近手机中文字幕大全| 亚洲精品aⅴ在线观看| 欧美精品亚洲一区二区| 自拍偷自拍亚洲精品老妇| 婷婷色av中文字幕| 亚洲怡红院男人天堂| 精品亚洲成a人片在线观看 | 免费看光身美女| 观看美女的网站| 国产精品一区二区三区四区免费观看| 免费观看性生交大片5| 亚洲国产日韩一区二区| 能在线免费看毛片的网站| 免费观看av网站的网址| 国产黄色视频一区二区在线观看| 成人毛片a级毛片在线播放| 国产精品久久久久成人av| 国产伦在线观看视频一区| 激情 狠狠 欧美| 久久国产精品男人的天堂亚洲 | 国产精品国产av在线观看| 青春草亚洲视频在线观看| 精品少妇黑人巨大在线播放| 卡戴珊不雅视频在线播放| 国产av精品麻豆| 男人和女人高潮做爰伦理| av网站免费在线观看视频| 欧美xxxx性猛交bbbb| 久久久午夜欧美精品| 一级av片app| 男女下面进入的视频免费午夜| 日日啪夜夜爽| 国产精品免费大片| 国产精品99久久久久久久久| 六月丁香七月| 成人毛片60女人毛片免费| 少妇人妻精品综合一区二区| 婷婷色综合大香蕉| 国产 一区 欧美 日韩| 久久久国产一区二区| 亚洲人成网站在线观看播放| 国产高清有码在线观看视频| 久久国产精品大桥未久av | 青春草视频在线免费观看| av一本久久久久| 亚洲久久久国产精品| 欧美成人a在线观看| 久久人妻熟女aⅴ| 国产黄色视频一区二区在线观看| 日本色播在线视频| 精品亚洲成a人片在线观看 | 黄色一级大片看看| 亚洲va在线va天堂va国产| 涩涩av久久男人的天堂| 性色av一级| 亚洲精品自拍成人| 日韩 亚洲 欧美在线| 日本黄大片高清| 我要看黄色一级片免费的| 国产男女超爽视频在线观看| 一级毛片久久久久久久久女| 青春草视频在线免费观看| 亚洲美女视频黄频| 蜜桃亚洲精品一区二区三区| 久久久亚洲精品成人影院| 亚洲精品国产av成人精品| 女性生殖器流出的白浆| 我要看日韩黄色一级片| 80岁老熟妇乱子伦牲交| 一本色道久久久久久精品综合| 国产v大片淫在线免费观看| 只有这里有精品99| 久久久久网色| 日本-黄色视频高清免费观看| 午夜老司机福利剧场| 久久鲁丝午夜福利片| 日本爱情动作片www.在线观看| 亚洲欧美日韩东京热| 欧美zozozo另类| 国产精品一区二区三区四区免费观看| av国产精品久久久久影院| 亚洲丝袜综合中文字幕| 十八禁网站网址无遮挡 | 十分钟在线观看高清视频www | 国产一级毛片在线| 亚洲精品,欧美精品| 国产一级毛片在线| 欧美bdsm另类| 亚洲精品中文字幕在线视频 | 亚洲久久久国产精品| 久久久精品94久久精品| 99久久精品一区二区三区| 国产精品欧美亚洲77777| 中文字幕久久专区| av免费在线看不卡| 涩涩av久久男人的天堂| 国产成人91sexporn| 国产深夜福利视频在线观看| 亚洲欧美成人综合另类久久久| 欧美成人午夜免费资源| 少妇人妻 视频| 丝瓜视频免费看黄片| 日韩不卡一区二区三区视频在线| 纯流量卡能插随身wifi吗| 精品久久久久久久久亚洲| 国产精品久久久久久久久免| 精华霜和精华液先用哪个| 亚洲av不卡在线观看| 视频区图区小说| 久久精品国产a三级三级三级| 草草在线视频免费看| 国产无遮挡羞羞视频在线观看| 日本av免费视频播放| 91精品国产九色| 最新中文字幕久久久久| 成年av动漫网址| 亚洲国产av新网站| 伊人久久国产一区二区| 中国三级夫妇交换| 国产色婷婷99| 国产男女内射视频| 少妇的逼水好多| 国产亚洲午夜精品一区二区久久| 国产精品国产av在线观看| 成人漫画全彩无遮挡| 日韩在线高清观看一区二区三区| 我要看日韩黄色一级片| 亚洲国产精品999| 日本av手机在线免费观看| 欧美人与善性xxx| 日韩av免费高清视频| 国产 一区 欧美 日韩| 久久久久久久久大av| 熟妇人妻不卡中文字幕| 亚洲美女视频黄频| av黄色大香蕉| 各种免费的搞黄视频| 黄色欧美视频在线观看| 久久久国产一区二区| 国产伦精品一区二区三区四那| 男女国产视频网站| 久久av网站| 久久6这里有精品| 亚洲av二区三区四区| 一级毛片我不卡| 男女免费视频国产| 狂野欧美激情性bbbbbb| 免费高清在线观看视频在线观看| 日本黄色片子视频| 婷婷色综合大香蕉| 老熟女久久久| 国产欧美另类精品又又久久亚洲欧美| 联通29元200g的流量卡| 一级毛片电影观看| 人人妻人人爽人人添夜夜欢视频 | 自拍欧美九色日韩亚洲蝌蚪91 | 国产探花极品一区二区| 久久99热这里只频精品6学生| 免费看光身美女| 我要看黄色一级片免费的| 成人毛片a级毛片在线播放| 亚洲精品成人av观看孕妇| 欧美亚洲 丝袜 人妻 在线| 成年免费大片在线观看| av国产免费在线观看| 在线亚洲精品国产二区图片欧美 | 精品久久久久久久末码| 国产av码专区亚洲av| 久久韩国三级中文字幕| 日韩成人伦理影院| 18禁在线播放成人免费| 免费看光身美女| 国产免费一区二区三区四区乱码| 久久人妻熟女aⅴ| 国产成人aa在线观看| av播播在线观看一区| 国产爽快片一区二区三区| 七月丁香在线播放| 联通29元200g的流量卡| 成人影院久久| 黄色配什么色好看| 午夜福利影视在线免费观看| 久久久久久九九精品二区国产| 国产精品一二三区在线看| 99热国产这里只有精品6| av在线app专区| 一边亲一边摸免费视频| 91精品一卡2卡3卡4卡| av国产精品久久久久影院| 99视频精品全部免费 在线| 热99国产精品久久久久久7| 国产极品天堂在线| 亚洲av不卡在线观看| 国产精品国产三级国产av玫瑰| 人妻 亚洲 视频| 99九九线精品视频在线观看视频| 精品酒店卫生间| 久久人人爽人人片av| 爱豆传媒免费全集在线观看| 黄色一级大片看看| 秋霞在线观看毛片| xxx大片免费视频| 国产黄片美女视频| 色视频www国产| 九色成人免费人妻av| 亚洲成人av在线免费| 嘟嘟电影网在线观看| a级毛片免费高清观看在线播放| 你懂的网址亚洲精品在线观看| 久久久久久久大尺度免费视频| 国产精品爽爽va在线观看网站| 日本-黄色视频高清免费观看| 成年av动漫网址| 一级二级三级毛片免费看| 国产免费一级a男人的天堂| 久久精品国产鲁丝片午夜精品| 校园人妻丝袜中文字幕| 日韩av在线免费看完整版不卡|