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

雙三相異步電機電流預測控制算法

2015-04-06 03:15:32柴建云孫旭東陸海峰
電工技術學報 2015年9期
關鍵詞:模塊化實驗

張 杰 柴建云 孫旭東 陸海峰

(電力系統及發電設備安全控制和仿真國家重點實驗室(清華大學) 北京 100084)

?

雙三相異步電機電流預測控制算法

張 杰 柴建云 孫旭東 陸海峰

(電力系統及發電設備安全控制和仿真國家重點實驗室(清華大學) 北京 100084)

針對一種半對稱式雙三相異步電機,提出兩種電流預測控制算法。一種為基本算法,另一種為在基本算法的基礎上通過引入合理近似而實現的模塊化算法。相對而言,后者更為簡潔,且易于在其他一些結構更為復雜的多相電機中進行推廣。此外,針對離散控制的延時對上述兩種算法帶來的穩定性問題進行分析,并給出改進方案。實驗結果表明,兩種電流預測控制算法均能實現良好的電流控制效果,同時該文針對穩定性問題的改進方法也是有效的。

雙三相異步電機 電流預測控制 模塊化 穩定性

0 引言

多相電機因具有轉矩脈動小、單個功率器件電流容量要求小、可容錯運行等優點,得到越來越廣泛的應用[1,2]。多相電機種類眾多,其中有一類為多Y型結構,即電機定子繞組由多組相互間錯開一定角度的三相繞組單元構成,如半對稱式雙三相電機[3]、四Y移15°十二相電機[4]等。

為了充分發揮多相電機的優勢,學者們從多相電機的高性能控制算法的角度開展了大量工作[5,6]。然而,對于控制算法中的電流內環,大多數學者只是簡單地給出PI調節器的方案。對于三相電機,盡管已有PI調節器參數的設計方法[7],但設計出的參數應用于實際系統時,往往需要根據不同工況做較大幅度的調整。而對于多相電機,尚未有學者給出系統化的PI調節器參數的設計方法。此外,隨著多相電機系統復雜程度的提高和PI調節器數目的增加,通過實際調試獲得一組性能良好的PI參數也變得較為困難。為此,一些學者考慮了其他形式的電流調節器。以半對稱式雙三相異步電機為例,文獻[8]采用了滯環比較器的方法,但存在開關頻率不固定的缺點。文獻[9]提出了使用PR調節器的思路,但仍面臨PR參數的整定問題。F.Barrero等[10-15]將電流預測控制(predictive current control,PCC)算法應用于電流內環,在獲得良好的電流控制效果的同時,降低了參數整定困難。PCC的數學本質是從可行解中尋找最優解的優化問題,按可行解的構成大致可分為兩類:一類是由逆變器開關狀態決定的有限個電壓矢量組成可行解;另一類是考慮SVPWM方案后,將逆變器電壓輸出限幅內所有的電壓值作為可行解。F.Barrero等選擇的是第一類思路。在雙三相異步電機領域,尚未有學者采用第二類思路。對于雙三相異步電機,每個控制周期需要考慮的電壓矢量從三相系統的8個增加至64個,運算壓力大大增加。文獻[10-15]的主要工作之一即為降低算法運算量。文獻[10,11]采用動態減少每個周期的可選電壓矢量數目的方案,文獻[12-15]直接將可選電壓矢量子集限定為12個最大電壓矢量和零矢量。文獻[10-15]的策略降低了PCC算法的運算量,獲得了良好的電流動態響應,但其電流高頻紋波還較大,尤其是諧波子空間電流的紋波水平還有待改善。

實際上,如果采用PCC算法的第二類思路,算法的運算量將近似隨電機相數增加呈線性增長,而非第一類思路的指數增長。并且對于諧波子空間的電流,也易于通過合理計算電壓指令而實現紋波更小的控制效果。本文首先基于PCC算法的第二類思路,提出雙三相異步電機的一種基本PCC算法。在此基礎上,考慮到雙三相異步電機繞組構成方面具有模塊化特點,為使算法更為簡潔,對電機模型引入了模塊化的合理近似,將雙三相異步電機近似按兩個獨立運行的三相異步電機模塊進行建模,通過分別實現每個模塊的PCC來實現雙三相異步電機整體的PCC。此外,還對這兩種PCC算法的穩定性問題進行了一處分析和改進。最后通過實驗研究,對比本文提出的兩種PCC算法和文獻[15]的PCC算法的電流控制效果,并驗證本文對PCC算法穩定性的一處改進的有效性。

1 雙三相異步電機的基本PCC算法

1.1 雙三相異步電機的數學模型

雙三相異步電機的繞組示意圖如圖1所示。在定子側,abc三相繞組Y型聯結,uvw三相繞組亦Y型聯結。abc三相繞組在空間位置上超前uvw三相繞組30°電角度,且兩套繞組的中性點相互隔離。兩套繞組沒有電的耦合,再忽略兩套定子繞組的互漏感,則兩套繞組僅通過氣隙磁鏈發生磁的耦合。

圖1 雙三相異步電機繞組示意圖Fig.1 Schematic of the windings of the dual three phase induction machine

采用文獻[3]提出的空間矢量分解方法對雙三相異步電機進行建模。其中,所用坐標變換式如式(1)所示

(1)

式中,x表示電壓u、電流i和磁鏈ψ等電氣量;xsa、xsb、xsc、xsu、xsv、xsw分別表示定子a、b、c、u、v、w各相的電氣量;xsα、xsβ、xsz1、xsz2分別表示在定子基波子空間的α軸和β軸及諧波子空間z1軸和z2軸分量。

將文獻[3]中的電機模型改寫為電阻,電感和反電勢串聯的形式[16]為

(2)

其中

(3)

式中,Rs、Rr、Lm、Lr、σLs、Lsl、Tr分別為電機的定子電阻、轉子電阻、互感、轉子電感、定子暫態電感、定子漏感和轉子時間常數;ψrα、ψrβ、ωr分別為電機的轉子αβ軸磁鏈和轉速。

1.2 基本PCC算法的實現

對式(2)所示的電機方程采用梯形公式,以Ts為控制周期,假定當前為n時刻,進行離散化,其結果如式(4)所示。其中,為了便于對反電勢的估計,將梯形公式離散化時的[eα(n)+eα(n-1)]/2和[eβ(n)+eβ(n-1)]/2分別用eα(n)和eβ(n)代替。

(4)

(5)

(6)

采用上述方法處理后,當前n時刻的電流采樣值將完全不出現在表達式中。這樣,綜合式(4)~式(6),并用反電勢的預測值代替實際值,最終可得到本文提出的基本PCC算法中電壓指令的表達式為

(7)

文獻[15]中基于有限開關狀態的PCC算法需要對每個可選電壓矢量進行一次計算,而本文方法只需要依據式(7)進行一次運算,運算量得到降低。由于式(7)是以下一周期實際電流跟蹤指令電流為依據進行電壓指令的計算,因而從電流的控制效果看,該方法優于以電流偏差最小為目標,在有限個電壓矢量中進行選擇的方法。

2 雙三相異步電機的模塊化PCC算法

2.1 雙三相異步電機的模塊化模型

分別取a相繞組軸線和u相繞組軸線為α軸和α′軸,β軸和β′軸分別由α軸和α′軸逆時針旋轉90°得到。兩套繞組分別以αβ軸和α′β′軸做clarke變換,如式(8)、式(9)

(8)

(9)

式中,x表示電壓u、電流i和磁鏈ψ等電氣量;xsα1、xsβ1、xsα′2、xsβ′2分別表示第一套繞組的電氣量x在定子αβ軸的分量和第二套繞組的電氣量x在α′β′軸的分量。對式(9)取式(10)所示的park變換,將第二套繞組在α′β′軸的電氣量變換到αβ軸。

(10)

雙三相異步電機穩態對稱運行時,時間上abc三相電流分別超前uvw三相30°,幅值上各相電流均相等,電壓和磁鏈的情況相同。而式(10)的另一個物理含義是分別將uvw相的電氣量的相位增加了30°,從而使得變換后的結果xsα2、xsβ2分別與xsα1、xsβ1的相位相同,幅值相同。綜合式(1)、式(8)、式(10)可得到

(11)

將式(11)的結果代入式(2)中可得

(12)

定義下述兩個變量

(13)

則式(12)可轉換為

(14)

(15)

考慮到電機兩套繞組結構參數相同且運行對稱的條件,近似有isα1=isα2,isβ1=isβ2,且pisα1=pisα2,pisβ1=pisβ2,則式(14)和式(15)分別轉換為

(16)

(17)

式(16)、式(17)即為雙三相異步電機的模塊化模型。下面給出上述近似處理的合理性:

2)近似處理帶來的偏差可歸結為電機模型的外界干擾項。PCC算法是一種電流閉環控制算法,對外界擾動具有良好的抑制能力,在實現電流閉環控制時,能夠抑制模型近似化處理所帶來的偏差。

2.2 模塊化PCC算法的實現

(18)

模塊化PCC算法使得PCC算法在由三相繞組為單元組成的多相電機中實現更為簡潔,且軟件上的模塊化將為硬件上的模塊化奠定基礎。本文的后續實驗將表明,盡管采用了一定的近似處理,模塊化的PCC算法能獲得與基本PCC算法基本一致的電流控制效果。

3 考慮穩定性問題后對PCC算法的改進

PCC算法的穩定性一般會受到離散控制的延時和參數偏差等兩方面因素的影響。雙三相異步電機的PCC算法主要使用的是定轉子的漏感參數,其在電機運行過程中,受鐵心飽和程度變化的影響較小,因而本文主要考慮離散控制的延時帶來的穩定性問題。討論穩定性問題一般需得到系統的傳遞函數,且傳遞函數中涉及的變量應能連續取值。對于文獻[15]中基于有限開關狀態的PCC算法,推導其傳遞函數的過程較為復雜,且涉及的電壓變量取值不連續,難以進行穩定性分析和改進,而本文提出的兩種PCC算法則不存在上述困難。

在分析PCC算法的穩定性問題時,一般認為電機的反電勢得到了準確預測,電壓指令計算式,即控制器模型和電機模型中的反電勢部分可相互抵消,不再考慮。在此基礎上,靜止坐標系中電機模型和控制器模型在不同坐標軸中分別具有相同的形式,因而僅以其中一個坐標軸的方程為例進行分析即可。

不考慮反電勢,取其中一個坐標軸的電機模型進行Z變換后為

(19)

對一個坐標軸的控制器模型進行Z變換后為

(20)

綜合式(19)和式(20)可得

(21)

一般有ReTs?Le[17],因而式(21)的極點趨近于0,為穩定系統,綜合式(19)和式(21),考慮U(z)與I*(z)的關系為

(22)

式(22)的極點除式(21)外的兩個接近0的極點外,還存在一個單位圓上的極點-1,因而系統臨界穩定,實際運行時則較易進入不穩定狀態。這就意味著盡管I(z)可以跟隨I*(z),但U(z)會出現幅值不斷增大的振蕩情況。

為解決該穩定性問題,本文首先從穩態入手,考慮到u(n+1),u(n)和u(n-1)之間在穩態時相差較小的特點,提出一種在PCC算法中用u(n+1)部分代替u(n)和u(n-1)的策略,來保證系統的穩定性。在確保算法的穩定性后,再考慮電流的控制效果是否仍能滿足要求。

采用該方法后,PCC算法的控制器模型由式(20)轉換為

(23)

式中,k1、k2分別表示u(n+1)代替u(n)和u(n-1)的比例。綜合式(19)和式(23),并忽略包含ReTs/Le的項,可得此時U(z)與I*(z)的關系式為

(24)

圖2 模最大極點的模與k1、k2的關系圖Fig.2 Relationship between the magnitude of the largest magnitude pole and k1、k2

觀察圖2可知,理論上只要取k1略大于0,模最大極點的模就小于1,系統進入穩定區域。實際中,需要為系統留有一定穩定裕度,因而k1、k2取值不能過小。圖3a給出了k2取0.1,k1取0.1~0.6時,H(z)=[Ts/(2Le)]×(U(z)/I*(z))的幅頻特性,其中離散周期Ts取 0.2 ms。可看出,在頻率較高的區域,H(z)幅值較大,且在2.5kHz處存在一個尖峰。H(z)在高頻區域較高的幅值會將反饋回路中的高頻擾動進行放大,對系統的穩定產生不利影響。隨著k1的增大,該尖峰逐漸降低,當k1大于0.4后,尖峰的減小已不再明顯。圖3b給出了k1取0.4,k2取0~0.4時,H(z)的幅頻特性。可看出,隨著k2的增加,尖峰有所變窄,但不改變2.5 kHz時的取值,且當k2大于0.1后,對尖峰的改善已不是很明顯。綜合上述考慮,本文在下文的實驗中取k1=0.4,k2=0.1。在實際中,該取值還可進行調整,但總體取值范圍不會變化很大。

圖3 (Ts/2Le)(U(z)/I*(z))的幅頻特性與k1、k2的關系圖Fig.3 Relationship between the magnitude-frequency characteristic of (Ts/2Le)(U(z)/I*(z)) and k1、k2

采用上述改進策略后,實際電流與指令電流間的關系由式(21)轉換為式(25),該式已忽略包含ReTs/Le的項。其伯德圖如圖4所示。圖4表明,采用上述方法改善穩定性后,仍能保證電流響應的帶寬在1 kHz以上,且對于最為關注的0~150 Hz范圍內的電流,其幅值衰減和相位滯后幾乎為0,因而總體上可認為電流響應仍滿足要求。

(25)

圖4 采用本文策略后,I(z)/I*(z)關系式的伯德圖Fig.4 The Bode diagram of I(z)/I*(z) with the proposed method

4 實驗結果

實驗在一臺2.5 kW雙三相異步電機機組上完成。電機參數見表1。實驗平臺如圖5所示。在圖5a中,該電機同軸連接了一臺由安川H1000型四象限變頻器控制的三相異步電機作為負載。在圖5b中,由兩套三相逆變器主電路構成雙三相逆變器的主電路,控制板選用TMS320C2812DSP,開關頻率5 kHz。實驗過程中利用CAN總線以5 kHz速率讀取并記錄DSP芯片內部的電流電壓等相關數據。

表1 實驗用雙三相異步電機參數Tab.1 Parameters of dual three phase induction machine

圖5 實驗平臺Fig.5 Experimental platform

本文實驗采用間接矢量控制方案,圖6給出了模塊化PCC算法的控制框圖。兩套繞組均分整個系統的轉矩電流指令和勵磁電流指令。兩套系統各自獨立實現自己的電流預測控制算法。對于基本PCC算法的控制框圖,其與傳統的雙三相異步電機矢量控制框圖類似,只需將電流PI調節器用預測電流控制器代替即可。

圖6 實驗系統框圖Fig.6 Diagram of the experimental system

圖7 3種電流預測控制算法的對比實驗結果Fig.7 Experimental results with the three PCC strategies

圖8 3種電流預測控制算法對比實驗結果(局部放大圖)Fig.8 Experimental results with the three PCC strategies (partial enlarged view)

圖9 參數有偏差時,模塊化PCC實驗結果Fig.9 Experimental results of the proposed modularized PCC strategy with inaccurate machine parameter

圖10 參數有偏差時,模塊化PCC實驗結果(局部放大圖)Fig.10 Experimental results of the proposed modularized PCC strategy with inaccurate parameter (partial enlarged view)

為了驗證第3節中關于穩定性問題的分析的正確性,本文給出了從穩定狀態向不穩定狀態過渡的實驗結果,如圖11所示。在該實驗中,電流指令不變,修正系數k1和k2分別由0.4和0.1逐漸減少至0。為使在系統失控狀態電流不會嚴重過電流,對輸出電壓取0.3 pu的限幅值。在0~0.2 s范圍內,k1和k2取前文的設定值,電流控制效果良好。在0.2~0.8 s范圍內,隨著k1和k2逐漸減小,圖11a的abc繞組的α軸電流波形isα1和圖11b的abc繞組的α軸電壓波形usα1均出現了明顯的高頻脈動,但電流并未完全失控。在0.8 s之后,k1和k2均減小為0,電流完全失控。該實驗結果表明,如果取k1=k2=0,即不采用改善穩定性的方法,PCC算法將不穩定。而k1取略大于0的值即可使算法穩定,但為確保足夠的穩定裕度,k1、k2取值不可過小。

圖11 隨k1,k2變化,從穩定過渡到不穩定的對比實驗結果Fig.11 Comparison experimental results of transition from stable to unstable conditions with the change of k1 and k2

5 結論

本文針對雙三相異步電機,提出了基本PCC算法和模塊化PCC算法兩種電流控制算法。在運算量有所減少的前提下,兩種算法實現了與文獻[15]中性能最優的一種PCC算法基本一致的動態性能,并且在減小電流紋波,尤其是諧波子空間的電流紋波方面有明顯改進。本文的模塊化PCC算法相對于基本PCC算法,采用了含有一定近似處理的模塊化電機模型。相關分析和實驗表明,該近似不會明顯影響PCC算法的性能。而模塊化PCC算法的思路更為簡潔,并可在其他一些繞組結構具有模塊化特點的多相電機中加以推廣。此外,本文針對基本PCC算法和模塊化PCC算法的穩定性問題的一處分析和改進的有效性也通過一組實驗得到了驗證。

[1] Levi E.Multiphase electric machines for variable-speed applications[J].IEEE Transactions on Industrial Electronics,2008,55(5):1893-1909.

[2] 莊朝暉,熊有倫,馬挺.多相感應電機變頻調速系統——回顧、現狀及展望[J].電氣傳動,2001,31(2):3-7. Zhuang Zhaohui,Xiong Youlun,Ma Ting.Multiphase variable-speed induction machine drives system—state of arts and trends[J].Electric Drive,2001,31(2):3-7.

[3] Zhao Yifan,Lipo Thomas A.Space vector PWM control of dual three-phase induction machine using vector space decomposition[J].IEEE Transactions on Industry Applications,1995,31(5):1100-1109.

[4] 謝衛,朱軍.十二相四Y移15°繞組同步電動機的數學模型及動態仿真[J].電工技術學報,2004,19(1):2-6. Xie Wei,Zhu Jun.Modeling and simulation of 12-phase synchronous motor with four Y-connected 3-phase symmetrical windings displaced in turn by 15°[J].Transactions of China Electrotechnical Society,2004,19(1):2-6.

[5] 楊金波,李鐵才,楊貴杰.一相開路雙三相永磁同步電機建模與控制[J].電工技術學報,2011,26(10):167-173. Yang Jinbo,Li Tiecai,Yang Guijie.Modeling and control of dual three-phase PMSM with one open phase[J].Transactions of China Electrotechnical Society,2011,26(10):167-173.

[6] Hatua K,Ranganathan V T.Direct torque control schemes for split-phase induction machine[J].IEEE Transactions on Industry Applications,2005,41(5):1243-1254.

[7] 李永東.交流電機數字控制系統[M].北京:機械工業出版社,2012.

[8] Kianinezhad R,Nahid B,Baghli L,et al.Aspects of current regulation in indirect field oriented control of dual three phase induction machines[C].IEEE International Conference on Industrial Technology,Mumbai,India,2006:933-938.

[9] Bojoi R,Levi E,Farina F,et al.Dual three-phase induction motor drive with digital current control in the stationary reference frame[J].IEE Proceedings-Electric Power Application,2006,153(1):129-139.

[10]Duran M J,Barrero F,Prieto J,et al.Predictive current control of dual three-phase drives using restrained search techniques and multi level voltage source inverters[C].2010 IEEE International Symposium on Industrial Electronics,Bari,Italy,2010:3171-3176.

[11]Duran M J,Prieto J,Barrero F,et al.Predictive current control of dual three-phase drives using restrained search techniques[J].IEEE Transactions on Industrial Electronics,2011,58(8):3253-3263.

[12]Barrero F,Arahal M R,Gregor R,et al.A proof of concept study of predictive current control for VSI-driven asymmetrical dual three-phase AC machines[J].IEEE Transactions on Industrial Electronics,2009,56(6):1937-1954.

[13]Barrero F,Arahal M R,Gregor R,et al.One-step modulation predictive current control method for the asymmetrical dual three-phase induction machine[J].IEEE Transactions on Industrial Electronics,2009,56(6):1974-1983.

[14]Gregor R,Barrero F,Toral S L.Predictive-space vector PWM current control method for asymmetrical dual three-phase induction motor drives[J].IET Electric Power Applications,2010,4(1):26-34.

[15]Barrero F,Prieto J,Levi E.An enhanced predictive current control method for asymmetrical six-phase motor drives[J].IEEE Transactions on Industrial Electronics,2011,58(8):3242-3252.

[16]Quang N P,Dittrich J A.Vector control of three-phase AC machines[M].Berlin,Heidelberg:Springer,2008.

[17]Guzinski J,Krzeminski Z.Advanced current regulated PWM inverter with simplified load model[J].Electric Power Components and Systems,2004,32(10):1031-1042.

Predictive Current Control Methods for Dual Three Phase Induction Machine

ZhangJieChaiJianyunSunXudongLuHaifeng

(State Key Laboratory of Security Control and Simulation of Power Systems and Large Generation Equipments Tsinghua University Beijing 100084 China)

For an asymmetrical dual three phase induction machine,two predictive current control methods are proposed.One of them is the basic method.Based on this,the other modularized method is suggested with reasonable assumption.Comparing to the basic method,the latter is more concise and easy to be applied to some other more sophisticated multi-phase machines.Besides,the stability issue of these two predictive current control strategies caused by the delay of the discrete control is analyzed,and one improvement solution is provided.Experimental results show that the two proposed predictive current control strategies can achieve a good current performance.Meanwhile,the effectiveness of the improvement on the stability issue is also validated.

Dual three phase induction machine,predictive current control,modularized,stability

國家自然科學基金(51277102)資助項目。

2014-11-15 改稿日期2015-03-12

TM321

張 杰 男,1987年生,博士研究生,研究方向為電動汽車驅動控制。(通信作者)

柴建云 男,1961年生,教授,博士生導師,研究方向為電機及其控制系統。

猜你喜歡
模塊化實驗
記一次有趣的實驗
模塊化自主水下機器人開發與應用
微型實驗里看“燃燒”
基于模塊化控制系統在一體化教學中的應用
模塊化住宅
現代裝飾(2020年5期)2020-05-30 13:01:56
做個怪怪長實驗
ACP100模塊化小型堆研發進展
中國核電(2017年2期)2017-08-11 08:00:56
從模塊化中得到的二氧化碳
流程工業(2017年4期)2017-06-21 06:29:52
模塊化VS大型工廠
流程工業(2017年4期)2017-06-21 06:29:50
NO與NO2相互轉化實驗的改進
主站蜘蛛池模板: 真实国产乱子伦高清| 亚洲综合亚洲国产尤物| 亚洲乱码在线播放| 草草线在成年免费视频2| 97成人在线视频| 午夜小视频在线| 亚洲天堂久久| 久久人人爽人人爽人人片aV东京热 | 亚洲成a人片| 成人国产小视频| 亚洲最大看欧美片网站地址| 一区二区日韩国产精久久| 欧美黄色网站在线看| 亚洲精品无码抽插日韩| 中文字幕丝袜一区二区| 92午夜福利影院一区二区三区| 99人妻碰碰碰久久久久禁片| 亚洲第一黄色网| 亚洲AV无码乱码在线观看裸奔 | 中文字幕 日韩 欧美| 欧美亚洲综合免费精品高清在线观看| 福利在线一区| 国产日本一区二区三区| 亚洲欧美日韩中文字幕在线| 欧美精品亚洲精品日韩专区| 亚洲高清日韩heyzo| 丁香婷婷综合激情| 午夜成人在线视频| 亚洲国产中文综合专区在| 欧美一级黄片一区2区| 日韩精品高清自在线| аⅴ资源中文在线天堂| 人妻少妇久久久久久97人妻| 久久久久国产精品嫩草影院| 精品无码国产自产野外拍在线| 亚洲人成色在线观看| 国产乱子伦精品视频| 国产无码在线调教| 精品一區二區久久久久久久網站 | 欧美在线视频不卡第一页| 国产国语一级毛片在线视频| 国产小视频免费观看| 热九九精品| 国产女人爽到高潮的免费视频 | 人人艹人人爽| 国产激情第一页| 影音先锋亚洲无码| 国产亚洲欧美在线中文bt天堂| 久久国产免费观看| 欧洲熟妇精品视频| a在线亚洲男人的天堂试看| 日韩中文字幕免费在线观看| 欧美色亚洲| 一本大道东京热无码av| 欧美不卡二区| 亚洲婷婷六月| 91免费观看视频| 成人欧美日韩| 午夜性刺激在线观看免费| 亚洲国产成人在线| 国产亚洲精| 国产毛片高清一级国语| 91精品国产91久久久久久三级| 亚洲精品国产综合99| 在线免费观看AV| 98精品全国免费观看视频| 国产自无码视频在线观看| 亚洲人成网站观看在线观看| 国产毛片一区| 三级毛片在线播放| 亚洲第一区欧美国产综合| 亚洲欧美日韩动漫| 国产a网站| 久久久久九九精品影院| 中国特黄美女一级视频| 免费一级成人毛片| 欧美日韩动态图| 精品国产成人国产在线| 一本一道波多野结衣一区二区| 国模极品一区二区三区| 亚洲二区视频| 九九久久99精品|