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

        ?

        基于BP神經(jīng)網(wǎng)絡(luò)和分數(shù)階PIλDμ的 VIENNA整流器控制策略

        2022-09-29 10:28:28楊旭紅方浩旭吳亞雄賈巍
        科學技術(shù)與工程 2022年22期
        關(guān)鍵詞:閉環(huán)控制內(nèi)環(huán)整流器

        楊旭紅, 方浩旭, 吳亞雄, 賈巍

        (1.上海電力大學自動化工程學院, 上海 200090; 2.上海太陽能工程技術(shù)研究中心有限公司, 上海 200241)

        三相Vienna整流器是由學者Kolar等[1]在1994年提出的。與傳統(tǒng)的整流器而言,Vienna整流器具有所使用的功率開關(guān)器件較少,控制電路簡單,電壓應力低和輸出電壓波形質(zhì)量較好等特點。同時,由于Vienna整流器的單向潮流特性,主要用于中、大功率場合的單向整流或兩級變流器的前級升壓整流模塊。由此,Vienna整流器被廣泛的應用于電動汽車、充電樁、航空航天等領(lǐng)域[2-3]。

        目前,中外學者主要的研究方向是Vienna整流器的控制策略問題。根據(jù)內(nèi)外環(huán)的不同,內(nèi)環(huán)主要是對電流的控制或者對功率的控制,外環(huán)則大多為對直流電壓的控制[4]。文獻[5-6]研究了Vienna整流器的大、小信號數(shù)學模型,并且根據(jù)所得模型設(shè)計了傳統(tǒng)比例積分(proportional integral,PI)控制器。盡管PI控制器能夠?qū)ienna整流器進行控制,但是由于Vienna整流器控制系統(tǒng)的非線性和耦合較強,在對誤差進行線性求和以及積分環(huán)節(jié)的飽和造成的電壓超調(diào)以及動態(tài)響應慢等缺點。因此,文獻[7]對內(nèi)環(huán)進行了改進,加入了Bang-bang控制器,電壓外環(huán)則依舊是對直流電壓的偏差信號采用PI控制器,增加了穩(wěn)定性。但是,Bang-bang控制器實際上是對偏差進行了限幅,這樣會降低控制系統(tǒng)的快速性。文獻[8]對電壓外環(huán)進行改進,引入了滑??刂破?,電流內(nèi)環(huán)則是對參考功率進行跟蹤。滑??刂剖且环N變結(jié)構(gòu)的非線性控制,在非線性控制系統(tǒng)中具有快速性好、抗干擾能力強等特點,但文獻[8]中針對滑??刂扑惴ǖ耐茖н^程較為簡單且只分析了直流側(cè)的情況。文獻[9]主要是在文獻[8]上做進一步的改進,將比例諧振(proportional resonant,PR)控制引入內(nèi)環(huán)控制中,去控制電流信號,但實際上與PI控制并無差別。文獻[10]基于滑??刂频碾妷和猸h(huán)和功率內(nèi)環(huán)的新型雙閉環(huán)結(jié)構(gòu)。但由于不是對電流進行直接控制,因此交流電流諧波(THD)較大。文獻[11-12]將無源控制應用于Vienna整流器中。但無源控制需要Vienna控制系統(tǒng)的精確數(shù)學模型,且對電阻電感等器件參數(shù)精度要求較高,而實際工業(yè)現(xiàn)場的電阻電感等參數(shù)由于氣候環(huán)境的變化會存在小幅波動,因此不適用于實際現(xiàn)場。文獻[13]將模型預測控制引入Vienna整流器中,但推導過程較為復雜且對計算機的性能要求較高,不易實現(xiàn)。

        此外,針對Vienna整流器中點電位問題,目前的解決方法主要有:①通過對調(diào)制過程中的小矢量進行控制[14-16];②通過對零序電壓注入的方式進行控制[17-19];③采用模型預測控制[20-21]。

        綜上,由于Vienna整流器控制系統(tǒng)是一個多入多出、耦合較強的非線性系統(tǒng),因此需要采取更復雜的控制技術(shù)對其進行控制。神經(jīng)網(wǎng)絡(luò)由于其原理和實現(xiàn)簡單且能擬合任意非線性函數(shù)被廣泛應用到各個領(lǐng)域當中[22-24]。分數(shù)階PIλ在積分環(huán)節(jié)加入了分數(shù)階理論,控制性能更好。

        因此,現(xiàn)將BP神經(jīng)網(wǎng)絡(luò)和分數(shù)階PIλ分別用于電壓外環(huán)和電流內(nèi)環(huán),該方案對系統(tǒng)的數(shù)學模型要求較低,魯棒性較強,快速性較好并且易于理解方便推導。同時,針對中點電位問題采用基于小矢量SVPWM調(diào)制策略,該調(diào)制策略能有效解決Vienna整流器上下橋臂電壓不平衡的問題。通過與PI控制進行對比,表明本文所提控制策略的優(yōu)越性。

        1 Vienna整流器數(shù)學模型

        圖1為三相Vienna整流器拓撲結(jié)構(gòu)圖,Vienna整流器由6個快速恢復二極管(VD1~VD6)、3個雙向功率開關(guān)管以及上下兩個電容(Cp,Cn)所構(gòu)成。

        ea、eb、ec分別為電網(wǎng)電壓;ia、ib、ic為電網(wǎng)電流;Sip(i=a,b,c)為電網(wǎng)電流;L為三相濾波電感;R為三相濾波電阻;udc為Vienna整流器直流側(cè)輸出電壓;Rf為負載電阻;VD1~VD6為快速恢復二級管;RL為負載電阻;ip、in 分別為輸出直流母線正向電流、負向 電流;iL為負載電流圖1 Vienna整流器拓撲結(jié)構(gòu)圖Fig.1 Vienna rectifier topology diagram

        為了簡化分析,假定所有器件均處于理想狀態(tài)。當電網(wǎng)處于理想狀態(tài)時,假定直流側(cè)上下電容大小相等,即Cp=Cn=C,則有uCp=uCn=1/2udc。根據(jù)Kirchhoff 定律,可得到abc坐標系下的數(shù)學模型:

        (1)

        式(1)中:Sip和Sin分別為正方向電流和負方向電流的開關(guān)函數(shù);uCp和uCn分別為電容Cp和Cn電壓值;uin為整流橋與直流側(cè)母線中點N的電壓;uNO為中點電壓;i=a,b,c。

        將式(1)變換到兩相旋轉(zhuǎn)坐標系,得

        (2)

        (3)

        式中:id和iq為電網(wǎng)電流ia、ib、ic通過Park變換下所得d軸、q軸分量;ed、eq為電網(wǎng)電壓ea、eb、ec在通過Park變換下所得d軸、q軸分量;ω為交流角頻率,ω=2πf;SP,d、SP,q、SN,d、SN,q分別為Sip、Sin在d、q坐標軸下的正方向和負方向電流的開關(guān)函數(shù)。

        根據(jù)式(2)和式(3)可得如圖2所示的等效電路。

        從圖2中可以看出,系統(tǒng)模型中存在電流耦合,因此在dq坐標系下的數(shù)學模型仍然是非線性的。因此,傳統(tǒng)的線性控制策略已經(jīng)不再適用與該系統(tǒng)。在設(shè)計控制策略時應包含解耦控制或者采用相應非線性控制算法。

        圖2 dq坐標系下的等效模型 Fig.2 Equivalent model in dq coordinate system

        2 Vienna整流器電壓電流雙閉環(huán)控制策略

        Vienna整流器控制系統(tǒng)主要是采用雙閉環(huán)控制,其中,外環(huán)對參考電壓進行控制,主要控制目標是使直流輸出電壓能夠快速地跟蹤給定電壓并且能夠使上下橋臂電容電壓值相等,其輸出值作為內(nèi)環(huán)控制系統(tǒng)的參考值。內(nèi)環(huán)控制系統(tǒng)為電流控制環(huán),其主要作用是保證實際電流能夠跟蹤給定電流,同時有功、無功功率達到穩(wěn)定狀態(tài)。

        2.1 電壓外環(huán)控制器設(shè)計

        2.1.1 BP神經(jīng)網(wǎng)絡(luò)

        神經(jīng)網(wǎng)絡(luò)是一種模仿動物神經(jīng)系統(tǒng)的智能算法。BP神經(jīng)網(wǎng)絡(luò)由于其原理簡單,算法較易實現(xiàn)等特點被廣泛運用于非線性系統(tǒng)中。BP神經(jīng)網(wǎng)絡(luò)一般分為三層,包括輸入層、隱含層、輸出層,對于非線性函數(shù)具有很好的逼近效果。

        BP神經(jīng)網(wǎng)絡(luò)輸入層:

        x=[x1x2…xn]T

        (4)

        BP神經(jīng)網(wǎng)絡(luò)權(quán)值:

        w=[w1w2…wm]T

        (5)

        BP神經(jīng)網(wǎng)絡(luò)隱含層神經(jīng)元Netin輸入:

        (6)

        BP神經(jīng)網(wǎng)絡(luò)輸出層:

        (7)

        從式(6)中可看出,BP神經(jīng)網(wǎng)絡(luò)隱含層全部式線性變換,這樣神經(jīng)網(wǎng)絡(luò)的多層就沒有了實際意義,因此,需要在輸出層中增加一個激活函數(shù),使之轉(zhuǎn)變?yōu)榉蔷€性關(guān)系,本文中采用Sigmoid函數(shù),其函數(shù)表達式為

        (8)

        2.1.2 基于BP神經(jīng)網(wǎng)絡(luò)的電壓外環(huán)設(shè)計

        從式(3)中可以看出,Vienna整流器具有較強的非線性,且存在耦合,因此采用BP神經(jīng)網(wǎng)絡(luò)進行控制。

        假定Vienna整流器交流測輸入功率與輸出功率相等,這樣可得

        (9)

        式(9)中:iL為直流側(cè)負載電流值。

        當中點電位平衡時,上下橋臂電容電壓相等,即ucp=ucn=1/2udc。因此,總有一個電容也在向負載供電,則負載電流表達式為

        (10)

        同時,控制系統(tǒng)也滿足:

        (11)

        式(11)中:idref、iqref分別為d、q軸電流的參考量。

        結(jié)合式(3)、式(9)~式(11),可以得非線性函數(shù)關(guān)系為

        (12)

        (13)

        式(13)中:kp和ki為控制器參數(shù)。式(12)可以轉(zhuǎn)變?yōu)?/p>

        (14)

        神經(jīng)網(wǎng)絡(luò)指標函數(shù)為

        (15)

        當系統(tǒng)的實際輸出uCp、uCn和直流電壓設(shè)定值uCp、uCn不相同時,誤差將通過隱含層傳輸?shù)捷斎雽?,并根?jù)每層獲得的誤差信號更新權(quán)重。重復調(diào)整后,誤差沿梯度方向減小,最小誤差對應的電流指令idref自行處理輸出,實現(xiàn)參考電流恒定輸出。

        根據(jù)式(15)的誤差函數(shù),可得隱含層權(quán)值的更新公式為

        網(wǎng)絡(luò)視頻行業(yè)極其燒錢早已不是什么秘密,從過去的群雄混戰(zhàn)到今天的優(yōu)愛騰(優(yōu)酷、愛奇藝和騰訊視頻)三分天下,最重要的競爭籌碼就是資本,這是一個沒錢就沒法玩的游戲,所以,也只能是BAT之間的游戲。

        (16)

        式(16)中:δ1j=ηEf′2w2jf′1;E為神經(jīng)網(wǎng)絡(luò)的指標函數(shù)E(k);f′1和f′2分別為隱含層和輸出層的激活函數(shù);w1j和w2j分別為輸出層和隱含層和權(quán)重。

        同理,輸出層權(quán)重更新公式為

        (17)

        因此,通過BP神經(jīng)網(wǎng)絡(luò)找到非線性函數(shù)g,從而得到電流參考值idref,另外,q軸的參考電流iqref=0,這樣系統(tǒng)中無功功率設(shè)定值則為0。BP神經(jīng)網(wǎng)絡(luò)的輸入如式(18)所示,采用三層前向結(jié)構(gòu),隱含層個數(shù)為12,BP神經(jīng)網(wǎng)絡(luò)模型結(jié)構(gòu)圖如圖3所示。

        (18)

        圖3 BP神經(jīng)網(wǎng)絡(luò)模型結(jié)構(gòu)圖Fig.3 Structure of BP neural network model

        2.2 電流內(nèi)環(huán)控制器設(shè)計

        2.2.1 分數(shù)階理論

        分數(shù)階微積分是指尋找任意階導數(shù)和積分的方法(可以是分數(shù)、無理數(shù),甚至是復數(shù))理論,其一般定義形式為

        (19)

        分數(shù)階微積分的計算方法目前主要有Grünwald-Letnikov定義和Riemann-Liouville定義,本文研究中采用Riemann-Liouville定義,其表達式為

        (20)

        對式(20)進行Laplace變換可得

        (21)

        由此,可以得出分數(shù)階系統(tǒng)的傳遞函數(shù)表達式為

        (22)

        式(2)中:G(s)為系統(tǒng)的傳遞函數(shù);Y(s)和U(s)分別為系統(tǒng)的輸出和輸入;s為微分算子;αi和βj為微分次數(shù);ai和bj為系數(shù)。

        2.2.2 分數(shù)階PIλDμ控制器

        分數(shù)階PIλDμ的傳遞函數(shù)式(23)所示,采用分數(shù)階PIλDμ控制器的結(jié)構(gòu)框圖如圖4所示。

        (23)

        分數(shù)階PIλDμ通過增加兩個可調(diào)參數(shù)λ和μ,實際上是對微分和導數(shù)求任意階次,而非是整數(shù)階導數(shù)和n重積分,從而增加了控制器的自由度和靈活性,同時通過改變參數(shù)λ和μ,也能改變系統(tǒng)的低頻段和高頻段的幅值斜率和相位角,因此也增加了系統(tǒng)的穩(wěn)定性。當λ=μ=0 時,為P控制器;當λ=1、μ=0 時,為PI控制器;當λ=0、μ=1時,為PD控制器;λ=μ=1時,為PID控制器。

        圖4 FOPIλDμ控制器的結(jié)構(gòu)框圖Fig.4 Block diagram of the fractional order PIλDμ controller

        2.2.3 基于分數(shù)階PIλDμ的電流內(nèi)環(huán)設(shè)計

        根據(jù)式(2),令

        (24)

        式(24)中:ud、uq為三相電壓ea、eb、ec在兩相旋轉(zhuǎn)坐標系下d軸、q軸分量。將式(24)代入式(2)可得

        (25)

        根據(jù)Laplace變換,將式(25)轉(zhuǎn)變?yōu)轭l域方程:

        (26)

        從式(26)中可以看出在dq坐標系下電流存在耦合,這樣增加了控制器的復雜性,因此,根據(jù)式(26)同時結(jié)合式(23),可以得到基于電流解耦的分數(shù)階PIλDμ控制器,如式(27)所示:

        (27)

        根據(jù)式(27)可以得到電流內(nèi)環(huán)的控制框圖,如圖5所示。

        圖5 基于分數(shù)階PIλ電流內(nèi)環(huán)的控制框圖Fig.5 Control block diagram based on fractional order PIλ current loop

        3 Vienna整流器控制系統(tǒng)

        3.1 Vienna整流器SVPWM控制

        Vienna整流器輸出側(cè)在上橋臂和下橋臂分別存在一個電容,因此會存在中點電位平衡的問題,若不去考慮中點電位平衡控制,交流電流的總諧波失真(total harmonic distortion, THD)顯著增加且會導致Vienna整流器的損害,故傳統(tǒng)的SVPWM算法并不能適用于該控制系統(tǒng)。

        在SVPWM控制中,小矢量主要是控制直流側(cè)電容的充放電從而直接影響中點電位的平衡。因此采用文獻[25]所提方法對小矢量進行控制,從而使中點電位平衡。

        3.2 Vienna整流器控制系統(tǒng)

        圖6為Vienna整流器控制系統(tǒng)。首先,從交流側(cè)獲取abc坐標系下的三相電壓電流信號uabc和iabc,通過Park變換轉(zhuǎn)變成為dq坐標系下的電流電壓分量udq和idq,其中d軸反應的是有功功率,q軸反應的是無功功率。Park變換所需要的角度信號是將交流測電壓電流信號輸入至鎖相環(huán)中所得。然后將其注入本文所提的BP神經(jīng)網(wǎng)絡(luò)電壓控制環(huán)和分數(shù)階PIλDμ電流控制環(huán)中,使其跟蹤設(shè)定值。接著再將電流環(huán)輸出的參考電壓信號注入SVPWM控制模塊中進行控制,并且在SVPWM中加入中點電位控制,最后使直流輸出電壓達到設(shè)定值且上下橋臂電壓達到平衡。

        4 仿真結(jié)果分析

        為了驗證文中所提算法的正確性和優(yōu)越性,首先根據(jù)圖1建立拓撲結(jié)構(gòu),接著將所提的算法應用于該拓撲結(jié)構(gòu)中并與傳統(tǒng)的PI控制進行對比。所建立仿真系統(tǒng)的主要參數(shù)如表1所示。

        圖6 Vienna整流器控制系統(tǒng)Fig.6 Vienna rectifier control system

        表1 Vienna整流器主要參數(shù)Table 1 Main parameters of Vienna rectifier

        圖7為直流側(cè)電壓波形圖,從圖7中可以看出雙PI控制和BP-FOPI控制均可以達到設(shè)定電壓600 V。BP-FOPI控制的超調(diào)量為3.6%,調(diào)節(jié)時間約為0.02 s;而雙PI閉環(huán)控制的超調(diào)量為17.2%,調(diào)節(jié)時間約為0.05 s,且在0.05 s之前波動較為頻繁。因此,仿真證明了相較于雙PI閉環(huán)控制,BP-FOPI控制具有更快的響應速度和更好的穩(wěn)定性。

        圖7 直流電壓Fig.7 DC voltage

        圖8 直流側(cè)上下橋臂電容電壓差Fig.8 Voltage difference between upper and lower bridge arms on the DC side

        圖8為直流側(cè)上下橋臂電容電壓差udiff,圖8(a)為在BP-FOPI控制下中點電位平衡波形圖,圖8(b)為在雙PI閉環(huán)控制下中點電位平衡波形圖。從圖8(a)中可以看出,在BP-FOPI控制下上下橋臂電容電壓差值大約在[-0.27,0.25]波動,中點電位平衡控制與BP-FOPI控制協(xié)調(diào)較好,波動幅度較??;從圖8(a)中可以看出,在BP-FOPI控制下直流側(cè)上下橋臂電容電壓差值大約在[-0.27,0.25]波動,中點電位平衡控制與BP-FOPI控制協(xié)調(diào)較好,波動幅度較小;然而圖8(b)中可以看出,在雙PI閉環(huán)控制下直流側(cè)上下橋臂電容電壓差值大約在[-0.63,0.33]波動,中點電位平衡控制與雙PI閉環(huán)控制協(xié)調(diào)較差,波動幅度較大。

        圖9為Vienna 整流器輸出功率波形圖,圖9(a)為Vienna整流器輸出的有功功率,圖9(b)為Vienna整流器輸出的無功功率。從圖9(a)中可以看出BP-FOPI控制和雙PI閉環(huán)控制的輸出功率最后都穩(wěn)定在36 kW,BP-FOPI控制的超調(diào)量為3.6%,大約在0.02 s時達到穩(wěn)定狀態(tài);而雙PI控制的超調(diào)量為4.73%,大約在0.05 s時達到穩(wěn)定狀態(tài),且在0.05 s之前波動較大。從圖9(b)中可以看出兩種控制方案下無功功率都能達到設(shè)定值0 Var,而BP-FOPI控制在啟動時的無功功率相較于雙PI控制策略較小,且達到穩(wěn)定時間的速度更快。因此,BP-FOPI控制的快速性和穩(wěn)定性更好,性能更優(yōu)。

        圖9 Vienna 整流器輸出功率Fig.9 Vienna rectifier output power

        圖10 BP-FOPI控制下電網(wǎng)電流波形圖和THDFig.10 THD and grid current under BP-FOPI

        圖10~圖13 展示了Vienna 整流器電網(wǎng)電流以及THD在兩種控制方案下的波形圖。從圖10中可以看出在BP-FOPI控制下的電流,啟動時的峰值電流約為140 A,且在0.02 s時達到了穩(wěn)定狀態(tài);而在PI雙閉環(huán)控制下的電流,啟動時峰值電流約為150 A,比BP-FOPI控制下的啟動電流大,設(shè)備更容易損壞,且電流在0.05 s后才達到穩(wěn)定狀態(tài),快速性不如BP-FOPI控制。通過對電網(wǎng)電流0.1~0.2 s的快速傅立葉變換(fast Fourier transform,F(xiàn)FT)分析可得:在BP-FOPI控制下的THD為1.07%,而在雙PI閉環(huán)控制下的THD值為1.71%。顯然可以看出,文中所提的控制方案可以大幅度降低電網(wǎng)電流的THD,證明了BP-FOPI控制的有效性。

        圖11 雙PI閉環(huán)控制下電網(wǎng)電流波形圖和THDFig.11 THD and grid current under PI control

        圖12 Vienna整流器輸出功率因數(shù)Fig.12 Vienna rectifer output power factor

        圖14為Vienna整流器直流側(cè)輸出功率因數(shù),功率因數(shù)的高低反映了整流器的工作效率,從圖14中可以看出,在兩種控制策略下Vienna整流器的輸出功率因數(shù)都能達到1,啟動時BP-FOPI控制下的功率因數(shù)最低值為98%,但是雙PI閉環(huán)控制下的功率因數(shù)最低值為92%,因此,BP-FOPI控制的下功率因數(shù)波動更小,整流器工作效率更高,動態(tài)性能更好。

        5 結(jié)論

        針對Vienna整流器提出了BP-FOPI控制策略,其中電壓外環(huán)采用BP神經(jīng)網(wǎng)絡(luò)對參考電流進行擬合,電流內(nèi)環(huán)采用分數(shù)階PIλ,提高了系統(tǒng)的穩(wěn)定性。對于中點電位平衡,采用了基于小矢量的SVPWM調(diào)制策略。通過軟件仿真實驗可以得到以下結(jié)論。

        (1)在MATLAB/Simulink仿真實驗中,通過對雙閉環(huán)PI控制進行對比可以得到文中所提的控制策略的優(yōu)越性。

        (2)文中所提的控制策略與中點電位控制上相較于雙PI閉環(huán)控制協(xié)調(diào)性更好,使中點電位波動更小。

        猜你喜歡
        閉環(huán)控制內(nèi)環(huán)整流器
        博物館文創(chuàng)產(chǎn)品設(shè)計的“內(nèi)環(huán)-外循”框架研究
        包裝工程(2023年16期)2023-08-25 11:39:16
        基于LMI的過渡態(tài)主控回路閉環(huán)控制律優(yōu)化設(shè)計
        三電平PWM整流器下的地鐵牽引供電系統(tǒng)探討
        科學家(2021年24期)2021-04-25 17:55:12
        適用于厚度在線測量的水壓閉環(huán)控制系統(tǒng)
        智能車競賽中的閉環(huán)控制算法應用分析
        電子制作(2018年11期)2018-08-04 03:25:58
        經(jīng)臍兩孔法腹腔鏡腹股溝疝內(nèi)環(huán)高位結(jié)扎加臍外側(cè)襞加強術(shù)治療小兒腹股溝斜疝*(附108例報告)
        經(jīng)臍微型腹腔鏡內(nèi)環(huán)高位結(jié)扎術(shù)聯(lián)合包皮環(huán)套術(shù)的臨床應用
        三相電壓型PWM 整流器研究
        PWM整流器啟動瞬時電流過沖抑制策略
        三相PWM整流器解耦與非解耦控制的對比
        欧美一片二片午夜福利在线快| 亚洲中文字幕在线观看| 中文字幕日韩有码国产| 麻豆69视频在线观看| 欧美午夜理伦三级在线观看| 国产又大又硬又粗| 国产免费av片无码永久免费| 亚洲熟女综合一区二区三区| 亚洲天堂在线播放| 免费观看久久精品日本视频| 亚洲一区二区三区在线观看蜜桃 | 国内国外日产一区二区| 中文字幕一区二区三区日日骚| 凌辱人妻中文字幕一区| 又紧又大又爽精品一区二区| 久久久久久久久蜜桃| 亚洲av成人精品日韩一区| 欧美综合区| 中文字幕一二区中文字幕| 加勒比东京热一区二区| 亚洲av色香蕉一区二区三区| 四虎国产精品永久在线国在线| 亚洲欧美精品91| 精品国免费一区二区三区| 精品女同av一区二区三区 | 亚洲免费国产中文字幕久久久| 亚洲av成人片在线观看| 人妻精品动漫h无码网站| japanese无码中文字幕| 亚洲精品一品二品av| av免费在线播放视频| 免费不卡在线观看av| 亚洲男人第一无码av网站| 国产乱妇乱子视频在播放| 久久久久亚洲AV无码专| 日本一区二区三区小视频| 一区二区三区亚洲免费| 国产视频自拍一区在线观看| 国产精品亚洲а∨无码播放不卡| 免费人成无码大片在线观看| 亚洲色欲大片AAA无码|