亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于級聯(lián)MRAS的PMSM參數(shù)在線辨識方法研究

        2015-07-11 06:10:30尹忠剛張延慶孫向東劉靜鐘彥儒
        電氣傳動 2015年4期
        關(guān)鍵詞:參考模型磁鏈級聯(lián)

        尹忠剛,張延慶,孫向東,劉靜,鐘彥儒

        (1.西安理工大學電氣工程系,陜西西安710048;2.電力設備電氣絕緣國家重點實驗室(西安交通大學),陜西西安710049)

        1 引言

        永磁同步電機由于具有高效率、高功率密度以及高動態(tài)性能等優(yōu)點[1],在高精度、微進給伺服系統(tǒng)中得到了大量應用,但電機的溫升與磁飽和等因素會導致電機參數(shù)發(fā)生變化,進而導致PMSM 控制性能的下降,因此,電機參數(shù)的在線辨識成為當前電機研究領域的熱點[2-4]。在永磁同步電機調(diào)速系統(tǒng)中,轉(zhuǎn)子的速度和位置是不可缺少的參數(shù),利用電流、電壓等物理量并通過一定的方法可以實現(xiàn)無速度傳感器調(diào)速系統(tǒng)中轉(zhuǎn)子速度和位置的估計。由于省去了速度傳感器,無速度傳感器電機調(diào)速系統(tǒng)具有結(jié)構(gòu)簡單、體積小、成本低、可靠性高等優(yōu)點。近年來,無速度傳感器矢量控制系統(tǒng)中PMSM 轉(zhuǎn)速和位置估算的主要方法可分為以下幾類:非線性觀測器[5]、滑模觀測器法[6-7]、高頻信號注入法[8-9]、擴展卡爾曼濾波器[10-11]、模型參考自適應等[12-21]。

        模型參考自適應(MRAS)由Schauder C.提出[12],是首次基于穩(wěn)定性理論設計交流電機轉(zhuǎn)速的辨識方法,狀態(tài)和速度的漸近收斂性由Lyapunov 穩(wěn)定性理論和Popov 超穩(wěn)定性理論保證,但MRAS的速度觀測以參考模型為基礎,而參考模型本身的參數(shù)準確程度直接影響速度辨識的精度[13-14]。文獻[15]對近年來提出的多種估算PMSM轉(zhuǎn)子位置和轉(zhuǎn)速的方法進行了介紹,并比較了各種方法的優(yōu)缺點。文獻[16]在文獻[15]的基礎上提出一種改進MRAS 算法,利用q 軸電流的估算誤差構(gòu)建自適應律對PMSM 進行轉(zhuǎn)速辨識,該算法結(jié)構(gòu)簡單,易于實現(xiàn)。文獻[17]為了實現(xiàn)永磁同步電機參數(shù)的在線實時估算,提出一種在同一模型中對定子電阻、定子電感和永磁體磁鏈進行辨識的自適應算法。文獻[18]提出一種變結(jié)構(gòu)MRAS觀測器,該觀測器結(jié)合變結(jié)構(gòu)控制與自適應控制理論,提高了系統(tǒng)對電機參數(shù)變化與負載擾動的魯棒性。文獻[19]提出一種基于雙參數(shù)MRAS 的感應電機無速度傳感器矢量控制策略,通過對定子電阻進行在線辨識,提高了系統(tǒng)低速帶載能力。該方法雖然目前應用在感應電機上,但對PMSM的多參數(shù)同時在線辨識具有一定指導意義。文獻[20]將擴展卡爾曼濾波器應用于基于MRAS 的PMSM 無傳感器直接轉(zhuǎn)矩控制控制系統(tǒng)中,實現(xiàn)了電機永磁體磁鏈、轉(zhuǎn)子位置和轉(zhuǎn)速的同時在線辨識或估算,整個系統(tǒng)具有較好的動靜態(tài)性能和穩(wěn)定性。

        本文提出了一種基于級聯(lián)模型參考自適應的永磁同步電機參數(shù)在線辨識方法,能夠?qū)崿F(xiàn)電機定子電阻、轉(zhuǎn)子磁鏈、轉(zhuǎn)子速度及其位置的同時在線辨識。首先建立級聯(lián)模型辨識的結(jié)構(gòu),并根據(jù)Popov 超穩(wěn)定性定理推導辨識算法,給出級聯(lián)MRAS 清晰的辨識步驟,最后在Matlab/Simulink的仿真環(huán)境下進行驗證。

        2 基于MRAS的PMSM參數(shù)辨識理論

        MRAS 的主要思想為:設計一個不含未知參數(shù)的參考模型與一個含有待辨識參數(shù)的可調(diào)模型,兩個模型被相同的外部輸入所激勵,并具有相同物理意義的輸出量。當兩個模型同時工作時,其輸出量的差值被送入自適應機構(gòu),通過調(diào)節(jié)可調(diào)模型中的待辨識參數(shù)使差值趨近于零,以達到參數(shù)辨識的目的,基本結(jié)構(gòu)如圖1所示。

        圖1 MRAS基本結(jié)構(gòu)Fig.1 Structure of model reference adaptive scheme

        對于永磁同步電機的參數(shù)辨識,可將圖1 中的可調(diào)模型選取為永磁同步電機dq軸(兩相旋轉(zhuǎn)坐標系)的狀態(tài)方程,待辨識參數(shù)選擇為電機參數(shù),而把電機運行時實際的dq軸電流值作為參考模型的輸出。根據(jù)參考模型與可調(diào)模型輸出量的誤差,通過設計合理的自適應機構(gòu)對可調(diào)模型中的電機參數(shù)值進行調(diào)整,當可調(diào)模型與參考模型的輸出相等時,待辨識的電機參數(shù)就會收斂到實際值,從而實現(xiàn)在線辨識電機參數(shù)的目標。

        基于模型參考自適應的永磁同步電機參數(shù)在線辨識,其主要目的是設計合理的自適應率,該自適應律不僅要保證參考模型與可調(diào)模型的輸出誤差收斂至零,還要保證待辨識參數(shù)的收斂。采用Popov 超穩(wěn)定性理論設計自適應律,其設計方法分為以下幾個步驟:

        1)把MRAS 系統(tǒng)等效轉(zhuǎn)化為非線性反饋系統(tǒng),該系統(tǒng)包含一個線性時不變前饋系統(tǒng)和一個非線性時變反饋系統(tǒng);

        2)設計一部分自適應規(guī)則使非線性時變反饋系統(tǒng)滿足Popov不等式;

        3)設計其余的自適應規(guī)則保證線性時不變前饋系統(tǒng)嚴格正實;

        4)把等效系統(tǒng)再還原成MRAS系統(tǒng)。

        對于基于MRAS 的PMSM 參數(shù)在線辨識系統(tǒng),MRAS 自身就是一個時變的非線性系統(tǒng),其穩(wěn)定性問題是系統(tǒng)固有的也是首要解決的問題,而基于Popov超穩(wěn)定性理論設計的自適應律能夠成功地用來設計穩(wěn)定的MRAS系統(tǒng),故在電機參數(shù)辨識領域得到了廣泛的應用。

        3 基于級聯(lián)MRAS 的PMSM 參數(shù)辨識方法

        3.1 基于級聯(lián)MRAS的PMSM參數(shù)辨識結(jié)構(gòu)

        本文提出一種基于級聯(lián)MRAS 的PMSM 參數(shù)在線辨識方法,對于級聯(lián)MRAS,分析如下:假設M 和N 為兩個不同的函數(shù)集,但其函數(shù)值相同,如果N包含P個參數(shù)而M未包含,則參考模型設為M,可調(diào)模型設為N,通過設計合適的自適應律辨識這P個參數(shù);另一方面,如果M包含Q個參數(shù)而N未包含,則參考模型設為N,可調(diào)模型設為M,通過設計合適的自適應律辨識這Q 個參數(shù)。待辨識參數(shù)在MRAS 中均為估計值,M 與N 級聯(lián)辨識這些參數(shù)。

        由基于傳統(tǒng)模型參考自適應的PMSM 參數(shù)辨識方法可知,辨識電機定子電阻和轉(zhuǎn)子磁鏈時,需要首先獲取電機轉(zhuǎn)速值;而辨識電機轉(zhuǎn)速和位置時,需要首先獲取電機定子電阻與轉(zhuǎn)子磁鏈值。本文提出一種基于級聯(lián)MRAS 的PMSM 參數(shù)在線辨識方法,該方法分為兩個部分,首先采用定子電阻和轉(zhuǎn)子磁鏈的離線辨識值進行轉(zhuǎn)速辨識,當估算轉(zhuǎn)速已經(jīng)穩(wěn)定且給定轉(zhuǎn)速不變時,再進行定子電阻與轉(zhuǎn)子磁鏈的在線辨識。在辨識定子電阻與轉(zhuǎn)子磁鏈時,使用到的辨識轉(zhuǎn)速值要對實際辨識值進行一階濾波,且濾波時間較大,從而保證轉(zhuǎn)速值的平滑穩(wěn)定。

        基于級聯(lián)MRAS 的PMSM 無速度傳感器矢量控制結(jié)構(gòu)圖如圖2 所示,虛線框內(nèi)為基于級聯(lián)MRAS的PMSM參數(shù)辨識模塊,其中,MRAS1辨識電機定子電阻與轉(zhuǎn)子磁鏈,MRAS2 辨識電機轉(zhuǎn)子的速度及其位置。系統(tǒng)采用3個PI調(diào)節(jié)器,其中轉(zhuǎn)速信息作為控制的反饋量與給定轉(zhuǎn)速的偏差經(jīng)過轉(zhuǎn)速PI調(diào)節(jié)器得到給定轉(zhuǎn)矩電流,而電流PI 調(diào)節(jié)器根據(jù)給定電流與反饋電流的偏差調(diào)節(jié)d,q軸定子電壓給定值,通過坐標變換,求得α,β 軸的定子電壓參考值后就可以發(fā)出PWM 波,控制逆變器輸出并驅(qū)動PMSM運轉(zhuǎn)。

        圖2 基于級聯(lián)MRAS的PMSM無速度傳感器矢量控制結(jié)構(gòu)圖Fig.2 The block diagram of speed sensorless vector control system for PMSM based on cascade MRAS

        3.2 轉(zhuǎn)子速度的辨識

        選取PMSM本身作為參考模型,在同步旋轉(zhuǎn)的dq坐標系中,以定子電流為狀態(tài)變量的電流模型可以表示為

        式中:ud,uq為定子電壓在dq 軸的分量;id,iq為定子電流在dq軸的分量;Rs,Ls分別為定子電阻、定子電感;ωr為轉(zhuǎn)子速度;Ψr為轉(zhuǎn)子磁鏈。

        將式(1)改寫為

        則由式(2)可得:

        將式(3)簡寫為

        其中

        根據(jù)式(3),可調(diào)模型描述為

        簡寫為

        將式(7)簡寫為

        式(8)為定子電流矢量的誤差方程,該誤差方程可由圖3所示的標準反饋系統(tǒng)來表示。

        圖3 等效非線性反饋系統(tǒng)Fig.3 Equivalent nonlinear variable feedback system

        圖3中,D 是增益矩陣,它將廣義誤差e 處理為用于自適應控制的另一矢量V,為簡化計算可取D=Ⅰ(單位矢量)。式(1)所示的電流方程原本是非線性方程,但在式(5)中,已將轉(zhuǎn)速ωr處理為一個時變參數(shù),對于數(shù)字化控制系統(tǒng),可以認為在每一采樣周期內(nèi)ωr是不變的。于是,式(5)就成為線性時不變方程,圖3 上半部虛線框內(nèi)即為一線性時不變前饋系統(tǒng)。在尋求自適應矢量V與反饋矢量W的關(guān)系前,先用一個非線性時變反饋環(huán)節(jié)來表示它們之間的關(guān)系,這樣,就得到如圖3所示的等效反饋系統(tǒng)。

        根據(jù)Popov 超穩(wěn)定性理論,該反饋系統(tǒng)穩(wěn)定的條件為:

        1)線性時不變前饋部分傳遞函數(shù)H(s)=(sⅠ-A)-1為嚴格正實;

        在MRAS漸進穩(wěn)定的條件下,通過逆向求解Popov積分不等式可得轉(zhuǎn)速自適應律為

        式中:id,iq為電機實測值由可調(diào)模型計算得到。

        轉(zhuǎn)子速度的自適應律如圖4所示。

        圖4 轉(zhuǎn)速自適應律Fig.4 Speed adaptive law

        圖4中,轉(zhuǎn)速誤差信息εωr經(jīng)過PI調(diào)節(jié)器得到估算速度信號會迫使可調(diào)模型的輸出與參考模型的輸出趨向一致,令估算定子電流逼近實際定子電流,也就使估算轉(zhuǎn)速逼近實際轉(zhuǎn)速ωr,從而達到轉(zhuǎn)速辨識的目的。

        3.3 定子電阻與轉(zhuǎn)子磁鏈的辨識

        基于MRAS 的PMSM 定子電阻與轉(zhuǎn)子磁鏈在線辨識方法,仍然選取PMSM本身作為參考模型,并采用式(1)所示的PMSM 定子電流數(shù)學模型,其辨識算法框圖如圖5所示,參考模型與可調(diào)模型的輸入量都是定子電壓在dq 軸的分量ud和uq,其輸出量均為定子電流在dq軸的分量id和iq。

        圖5 基于MRAS的PMSM定子電阻與轉(zhuǎn)子磁鏈辨識算法框圖Fig.5 Block diagram of PMSM stator resistor and rotor flux parameters identification based on MRAS

        將式(1)所示的PMSM定子電流數(shù)學模型作為參考模型時,可調(diào)模型可以表示如下:

        采用與轉(zhuǎn)速自適應律相似的推導過程,根據(jù)Popov 超穩(wěn)定性理論,可以得到定子電阻與轉(zhuǎn)子磁鏈的辨識算法為

        3.4 基于級聯(lián)MRAS的PMSM參數(shù)辨識步驟

        由式(9)、式(11)、式(12)以及MRAS系統(tǒng)的構(gòu)成可知,辨識PMSM定子電阻與轉(zhuǎn)子磁鏈時需要獲取電機轉(zhuǎn)速值,而另一方面,辨識電機轉(zhuǎn)速時需要獲取定子電阻與轉(zhuǎn)子磁鏈的值?;诩壜?lián)MRAS 理論,電機轉(zhuǎn)速、定子電阻與轉(zhuǎn)子磁鏈的辨識步驟如下:

        1)首先使用定子電阻和轉(zhuǎn)子磁鏈的離線辨識值并通過式(9)計算出轉(zhuǎn)速辨識結(jié)果

        4 仿真驗證及結(jié)論

        本文采用Matlab/Simulink 仿真來驗證所提出的基于級聯(lián)MRAS 的PMSM 參數(shù)在線辨識方法的正確性和有效性,模型中PMSM 的參數(shù)為:額定功率PN=1.5 kW,額定頻率fN=50 Hz,額定轉(zhuǎn)速ωN=1 500 r/min,定子電阻Rs=4.96 Ω,d 軸電感Ld=8.5 mH,q 軸電感Lq=8.5 mH,極對數(shù)pN=2,轉(zhuǎn)子磁鏈Ψr=0.375 Wb。

        圖6是級聯(lián)MRAS輸出的辨識轉(zhuǎn)速與電機實測轉(zhuǎn)速的對比。從圖6中可以看出辨識轉(zhuǎn)速可以很好地跟蹤實測轉(zhuǎn)速,級聯(lián)MRAS在辨識過程中表現(xiàn)出良好的動、靜態(tài)性能。

        圖6 轉(zhuǎn)速的辨識值和實測值Fig.6 Identified and actual value of rotor speed

        圖7和圖8是級聯(lián)MRAS辨識的定子電阻與轉(zhuǎn)子磁鏈波形,在轉(zhuǎn)速辨識穩(wěn)定且轉(zhuǎn)速給定不變的情況下,t=1 s時啟動定子電阻與轉(zhuǎn)子磁鏈的在線辨識算法。從圖7與圖8中可以看出基于級聯(lián)MRAS 的定子電阻辨識值其辨識誤差約為0.2%;轉(zhuǎn)子磁鏈辨識值其辨識誤差約為0.1%,驗證了該在線辨識算法的有效性。

        圖7 定子電阻的辨識值和實測值Fig.7 Identified and actual value of stator resistance

        圖8 轉(zhuǎn)子磁鏈的辨識值和實測值Fig.8 Identified and actual value of rotor flux

        為了驗證級聯(lián)MRAS 在電機參數(shù)變化時的跟蹤辨識能力,對電機參數(shù)給定做如下處理:啟動時轉(zhuǎn)速給定為800 r/min,定子電阻與轉(zhuǎn)子磁鏈采用離線值;t=1 s 時,啟用定子電阻與轉(zhuǎn)子磁鏈在線辨識算法;t=2 s 時,定子電阻以一定規(guī)律遞增至真實值的120%,即5.952 Ω;t=4 s 時,轉(zhuǎn)速給定階躍降至100 r/min;t=6 s 時,轉(zhuǎn)速給定階躍增至1 500 r/min,此后保持轉(zhuǎn)速給定不變;t=7 s時,定子電阻以一定規(guī)律遞減至真實值的80%,即3.968 Ω,同時轉(zhuǎn)子磁鏈以一定規(guī)律遞減至真實值的50%,即0.187 5 Wb,仿真波形如圖9~圖11所示。

        從以上PMSM參數(shù)辨識結(jié)果可以看出,級聯(lián)MRAS 算法辨識的轉(zhuǎn)速可以快速跟蹤電機實測轉(zhuǎn)速,辨識的定子電阻與轉(zhuǎn)子磁鏈也可以較好地跟蹤電機參數(shù)變化,并且系統(tǒng)具有良好的動、靜態(tài)性能。

        圖9 轉(zhuǎn)速變化時的辨識結(jié)果Fig.9 Identified results when given rotor speed varies

        圖10 定子電阻變化時的辨識結(jié)果Fig.10 Identified results when stator resistance varies

        圖11 轉(zhuǎn)子磁鏈變化時的辨識結(jié)果Fig.11 Identified results when rotor flux varies

        提出了一種基于級聯(lián)模型參考自適應的永磁同步電機參數(shù)在線辨識方法,該方法結(jié)構(gòu)簡單,辨識步驟明晰,基于Popov超穩(wěn)定性理論建立了參數(shù)的辨識模型并推導出待辨識參數(shù)的自適應律。仿真結(jié)果表明,級聯(lián)MRAS能同時對電機的轉(zhuǎn)子速度、定子電阻與轉(zhuǎn)子磁鏈進行在線辨識,具有快速收斂性和較小的誤差。

        [1]Pillay P,Krishnan R.Application Characteristics of Permanent Magnet Synchronous and Brushless DC Motor for Servo Drives[J].IEEE Transactions on Industry Application,1991,27(5):986-996.

        [2]胡升,史婷娜.基于模型參考自適應的SPMSM參數(shù)在線辨識策略[J].電氣傳動,2014,44(1):13-16.

        [3]楊永立,張云龍,陳智剛,等.基于參數(shù)辨識的PMSM電流環(huán)在線自適應控制方法[J].電工技術(shù)學報,2012,27(3):86-91.

        [4]梁驕雁,胡育文,魯文其.永磁交流伺服系統(tǒng)的抗擾動自適應控制[J].電工技術(shù)學報,2011,26(10):174-180.

        [5]Lee Junggi,Hong Jinseok,Nam Kwanghee. Sensorless Control of Surface-mount Permanent Magnet Synchronous Motors Based on a Nonlinear Observer[J].IEEE Transactions on Power Electronics,2010,25(2):290-297.

        [6]Chi Song,Zhang Zheng,Xu Longya.Sliding-mode Sensorless Control of Direct-drive PM Synchronous Motors for Washing Machine Applications[J].IEEE Transactions on Industry Applications,2009,45(2):582-590.

        [7]魯文其,黃文新,胡育文. 永磁同步電動機新型滑模觀測器無傳感器控制[J].控制理論與應用,2009,26(4):429-432.

        [8]Bolognani S,Calligaro S,Petrella R.Sensorless Control of IPM Motors in the Low-speed Range and at Standstill by HF Injection and DFT Processing[J].IEEE Transactions on Industry Applications,2011,47(1):96-104.

        [9]Lu K,Vetuschi M,Rasmussen P O.Determination of High-frequency d-and q-axis Inductances for Surface-mounted Permanent Magnet Synchronous Machines[J].IEEE Transactions on Instrumentation and Measurement,2010,59(9):2376-2382.

        [10]陳振,劉向東,靳永強.采用擴展卡爾曼濾波磁鏈觀測器的永磁同步電機直接轉(zhuǎn)矩控制[J].中國電機工程學報,2008,28(33):75-81.

        [11]尹忠剛,張瑞峰,鐘彥儒,等.基于抗差擴展卡爾曼濾波器的永磁同步電機轉(zhuǎn)速估計策略[J].控制理論與應用,2012,29(7):921-927.

        [12]Schauder C. Adaptive Speed Identification for Vector Control of Induction Motors Without Rotational Transducers[J].IEEE TransactionsonIndustryApplications,1992,28(5):1054-1061.

        [13]WEN P,LUTW. Decoupling Control of a Twin Rotor Mimo Systemusing Robust Deadbeat Control Technique[J].Control Theory&Applications,IET,2008,2(1):999-1007.

        [14]Young S K,Sang K K.MRAS Based Sensorless Control of Permanent Magnet Synchronous Motor[C]//IEEE SICE Annual Conference,F(xiàn)ukui,Japan,2003(2):1632-1637.

        [15]Liang Yan,Li Yongdong.Sensorless Control of PM Synchronous Motors Based on MRAS Method and Initial Position Estimation[C]//Electrical Machines and Systems,2003(1):96-99.

        [16]齊放,鄧智泉,仇志堅,等.基于MRAS 的永磁同步電機無速度傳感器[J].電工技術(shù)學報,2007,22(4):53-58.

        [17]安群濤,孫力,趙克.一種永磁同步電機參數(shù)的自適應在線辨識方法[J].電工技術(shù)學報,2008,23(6):31-36.

        [18]王慶龍,張崇巍,張興.基于變結(jié)構(gòu)模型參考自適應系統(tǒng)的永磁同步電機轉(zhuǎn)速辨識[J].中國電機工程學報,2008,28(9):71-75.

        [19]尹忠剛,劉靜,鐘彥儒,等.基于雙參數(shù)模型參考自適應的感應電機無速度傳感器矢量控制低速性能[J].電工技術(shù)學報,2012,27(7):124-130.

        [20]史宇超,孫凱,馬鴻雁,等.內(nèi)埋式永磁同步電機永磁磁鏈的在線辨識[J].電工技術(shù)學報,2011,26(9):48-53.

        [21]Li Zhen,Xu Longya.Sensorless Field Orientation Control of Induction Machines Based on a Mutual MRAS Scheme[J].IEEE Transactions on Industrial Electronics,1998,45(5):824-831.

        猜你喜歡
        參考模型磁鏈級聯(lián)
        級聯(lián)LDPC碼的STBC-OFDM系統(tǒng)
        電子制作(2016年15期)2017-01-15 13:39:09
        基于環(huán)境的軍事信息系統(tǒng)需求參考模型
        語義網(wǎng)絡P2P參考模型的查詢過程構(gòu)建
        基于級聯(lián)MUSIC的面陣中的二維DOA估計算法
        一種弱磁擴速下的異步電機磁鏈觀測和速度辨識
        一種基于簡化MRAS無速度傳感器的永磁電機EKF磁鏈辨識
        基于虛擬磁鏈的STATCOM直接功率控制研究
        LCL濾波器在6kV級聯(lián)STATCOM中的應用
        電測與儀表(2014年1期)2014-04-04 12:00:34
        H橋級聯(lián)型STATCOM的控制策略研究
        電測與儀表(2014年1期)2014-04-04 12:00:28
        基于虛擬磁鏈的APF準諧振控制研究
        亚洲国产精品久久艾草| 久久久精品人妻一区二区三区日本 | h国产视频| av黄片免费在线观看| 麻婆视频在线免费观看| 亚洲av无码专区在线观看成人| 小sao货水好多真紧h视频| 久久久久成人精品免费播放网站| 视频国产一区二区在线| 免费成人电影在线观看| 欧美日韩亚洲tv不卡久久| 国产精品白浆无码流出| 综合成人亚洲网友偷自拍| 欧美国产亚洲日韩在线二区| 欧美freesex黑人又粗又大| 久久久久久久久久免免费精品| 大香蕉视频在线青青草| 亚洲色欲久久久综合网东京热| 丰满爆乳一区二区三区| 99在线无码精品秘 人口| 人妖在线一区二区三区| 中文字幕在线日亚洲9| 久久免费网国产AⅤ| 性色av手机在线观看| 综合国产婷婷精品久久99之一 | 亚洲伊人免费综合网站| 少妇被黑人嗷嗷大叫视频| 欧美俄罗斯40老熟妇| 制服丝袜视频国产一区| 亚洲中文高清乱码av中文| 国产熟妇疯狂4p交在线播放| 精品人妻潮喷久久久又裸又黄| 日本一区免费喷水| 一二三区亚洲av偷拍| 国产午夜精品一区二区三区嫩草| YW亚洲AV无码乱码在线观看| 久久91精品国产91久久跳舞| 一本色道久久88加勒比—综合| 日韩精品一区二区三区视频| 在线播放中文字幕一区二区三区| 日韩精品久久中文字幕|