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

    基于QPSO 和雙線性插值的混合儲能系統(tǒng)改進(jìn)解耦控制

    2022-02-25 14:05:58曹以龍涂少博帥祿瑋江友華
    電源學(xué)報(bào) 2022年1期
    關(guān)鍵詞:插值穩(wěn)態(tài)表格

    曹以龍,涂少博,帥祿瑋,江友華

    (上海電力大學(xué)電子與信息工程學(xué)院,上海 200090)

    隨著能源與環(huán)境問題的日益加劇,可再生能源存在巨大的發(fā)展空間。新能源發(fā)電系統(tǒng)中,為減小太陽能、風(fēng)能等可再生能源由于發(fā)電不連續(xù)、不穩(wěn)定等問題對系統(tǒng)造成的干擾,需要配備一定容量的儲能單元。蓄電池與超級電容組成的混合儲能系統(tǒng)相較于單一的儲能單元,具有特性互補(bǔ)、可延長電池使用壽命等優(yōu)勢,在實(shí)際中應(yīng)用較廣。傳統(tǒng)的混合儲能單元通過多個(gè)分立的DC/DC 變換器接入發(fā)電系統(tǒng),這樣導(dǎo)致系統(tǒng)結(jié)構(gòu)復(fù)雜,不易集中管理。多端口DC/DC變換器可將多個(gè)儲能單元有效結(jié)合在一起,提高系統(tǒng)的功率密度,具有電路拓?fù)浜唵巍⒊杀据^低、易于實(shí)現(xiàn)系統(tǒng)集中控制和能量管理等優(yōu)點(diǎn)[1]。

    完全隔離型多端口DC/DC 變換器通過一個(gè)多繞組高頻變壓器將各端口進(jìn)行連接,所有端口均實(shí)現(xiàn)電氣隔離,因此實(shí)際應(yīng)用更加廣泛。但是由于漏感的存在,使得各端口功率通過繞組磁通疊加后產(chǎn)生耦合關(guān)系,導(dǎo)致各個(gè)變量的控制環(huán)路設(shè)計(jì)較為復(fù)雜[2]。文獻(xiàn)[3]通過設(shè)置控制環(huán)路的帶寬,引入了新的解耦控制變量,大大降低了兩個(gè)控制變量之間的耦合,但是該方法并沒有實(shí)現(xiàn)完全解耦。文獻(xiàn)[4]推導(dǎo)了解耦條件下功率傳輸?shù)牡葍r(jià)關(guān)系,進(jìn)而得到解耦變量的關(guān)系表達(dá)式,實(shí)現(xiàn)解耦,該方法模型過于簡單,控制精度差。文獻(xiàn)[5]通過離線訓(xùn)練得到解耦模型,采用神經(jīng)網(wǎng)絡(luò)在線求得控制參數(shù)。但是神經(jīng)網(wǎng)絡(luò)運(yùn)算量過大,難以實(shí)現(xiàn)實(shí)時(shí)控制。

    本文針蓄電池與超級電容混合儲能三端口變換器的耦合效應(yīng),采用改進(jìn)解耦控制策略,為提高系統(tǒng)解耦精度,在傳統(tǒng)解耦控制的基礎(chǔ)上考慮變換器軟開關(guān)范圍以及器件損耗對解耦參數(shù)的影響。由于在解耦過程中需要實(shí)時(shí)求解一個(gè)非線性方程組,為避免實(shí)時(shí)計(jì)算的壓力,采用量子粒子群優(yōu)化QPSO(quantum-behaved particle swarm optimization)算法對非線性解耦方程組進(jìn)行離線求解,得到解耦參數(shù)表,存入控制系統(tǒng)當(dāng)中。為降低控制系統(tǒng)的內(nèi)存消耗,利用雙線性插值算法減少解耦參數(shù)的存儲量,以適用于實(shí)際工程應(yīng)用。

    1 混合儲能三端口DC/DC 變換器

    圖1 為三端口有源全橋TAB(triple active bridge)變換器。圖中:T 為三繞組高頻變壓器,其匝數(shù)比為N1∶N2∶N3,根據(jù)實(shí)際變壓器模型可將其等效為漏感L與勵磁電感為變換器串聯(lián)電感,由于實(shí)際情況下變壓器漏感遠(yuǎn)小于串聯(lián)電感,故可將L1、L2、L3視為該變換器的功率傳輸電感;S1~S12為MOSFET 開關(guān)管為電感電流;C1、C2、C3為濾波電容為橋臂中點(diǎn)電壓;V1為蓄電池電壓,V2為母線電壓,V3為超級電容電壓;i1、i2、i3為端口電流;R 為負(fù)載。變換器所有開關(guān)管頻率均為f,每個(gè)端口的驅(qū)動信號占空比均為50%,但存在一定相位差(即移相角)。設(shè)定端口1、2 之間的相位差為φ12,端口1、3 之間的相位差為φ13,端口2、3 之間的相位差為φ23,其范圍均為(-π/2~π/2)。通過改變φ12、φ13的大小和方向,即可調(diào)節(jié)功率的大小和流向。

    圖1 混合儲能三端口DC/DC 變換器Fig.1 Hybrid energy storage three-port DC/DC converter

    1.1 等效電路

    圖2 TAB 變換器等效電路Fig.2 Equivalent circuits of TAB converter

    根據(jù)Y-△轉(zhuǎn)換原理,可得△型等效電路,如圖2(b)所示,各等效電感可表示為

    當(dāng)φ12>φ13時(shí),令t0=0、t1=φ13/(2πf)、t2=φ12/(2πf)、t3=1/(2f),由電磁感應(yīng)定律,變換器半個(gè)周期內(nèi)電感電流可表示為

    1.2 軟開關(guān)特性

    文獻(xiàn)[6]根據(jù)TAB 的Y 型等效電路,將一個(gè)開關(guān)周期分為12 個(gè)子區(qū)間,并對每個(gè)區(qū)間內(nèi)電路開關(guān)管的狀態(tài)和電流的流向進(jìn)行了分析。由其工作模態(tài)分析可知,在移相控制下該電路易于實(shí)現(xiàn)零電壓開通,但無法實(shí)現(xiàn)零電流關(guān)斷。零電壓開通的具體條件可表述為,在開關(guān)管導(dǎo)通之前其體二極管已經(jīng)導(dǎo)通,將開關(guān)管兩端的電壓箝位在0。令iL1(t0)、iL2(t2)、(t1)分別表示橋臂中點(diǎn)電壓上升沿時(shí)刻的各端口電感電流的表達(dá)式。根據(jù)△型等效電路,軟開關(guān)條件可表示為

    圖3 軟開關(guān)范圍三維圖Fig.3 3D maps of soft switching range

    1.3 開關(guān)損耗

    以蓄電池所在端口1 為例進(jìn)行分析。開關(guān)管在工作過程中的損耗包括導(dǎo)通損耗、關(guān)斷損耗和開關(guān)損耗3 種。針對本文研究的MOSFET 而言,其導(dǎo)通狀態(tài)下可以等效為一個(gè)電阻Ron,因此其導(dǎo)通損耗可表示為

    式中,ID為流過開關(guān)管的電流。

    根據(jù)式(2),令L12=L13=LΔ,則開關(guān)管導(dǎo)通時(shí)間內(nèi)電感電流可表示為

    進(jìn)一步推導(dǎo)可得端口1 的開關(guān)管導(dǎo)通損耗為

    對于關(guān)斷損耗,在MOSFET 處于關(guān)斷狀態(tài)時(shí),該器件可以等效為一個(gè)無窮大的電阻,此時(shí)漏極電流很小,關(guān)斷損耗就很小,因此可忽略不計(jì)。根據(jù)文獻(xiàn)[7],開關(guān)管的開關(guān)損耗等于其開關(guān)過程中兩端電壓與電流的乘積,計(jì)算公式可表示為

    式中:VDS(on)、VDS(off)、IDS(on)、IDS(off)為開關(guān)管導(dǎo)通與關(guān)斷時(shí)兩端電壓與電流;tup、tfall為開關(guān)管上升、下降時(shí)間。根據(jù)對稱性,端口1 的開關(guān)損耗可表示為

    由軟開關(guān)特性分析,變換器能夠在一定范圍內(nèi)實(shí)現(xiàn)軟開關(guān),此時(shí)開關(guān)管的開關(guān)損耗可忽略不計(jì)。則開關(guān)損耗可進(jìn)一步表示為

    通過以上分析,可將端口1 的開關(guān)管損耗最終表示為

    同理可得,端口2、3 的開關(guān)管損耗為

    2 改進(jìn)解耦模型

    根據(jù)文獻(xiàn)[8],采用基波分析法可得到各端口之間功率傳輸公式為

    以i1、i3為系統(tǒng)控制變量,可得系統(tǒng)在穩(wěn)態(tài)工作點(diǎn)R 處的小信號模型為

    式中:Δi1、Δi3為端口小信號電流;G 為傳輸矩陣,其G11、G12、G21、G22的表達(dá)式為

    文獻(xiàn)[9]給出了傳統(tǒng)TAB 變換器解耦控制方法下系統(tǒng)穩(wěn)態(tài)工作點(diǎn)的電流i1R、i3R表達(dá)式,即

    由式(15)知,該方法在推導(dǎo)系統(tǒng)穩(wěn)態(tài)點(diǎn)的電流表達(dá)式時(shí),僅考慮理想情況下的能量傳輸,導(dǎo)致系統(tǒng)穩(wěn)態(tài)工作點(diǎn)求解精度不高,系統(tǒng)整體控制效果較差。為提高控制效果,本文將軟開關(guān)范圍和開關(guān)管損耗考慮到系統(tǒng)解耦控制中,根據(jù)式(10)~式(12)可得改進(jìn)后穩(wěn)態(tài)工作點(diǎn)處的電流分別為

    由式(13)可知,端口1 和端口3 之間存在耦合,控制變量i1、i3不能夠獨(dú)立控制,需要進(jìn)行解耦。本文采用對角解耦法來解除耦合,在控制回路中增加G的逆矩陣H 作為解耦環(huán)節(jié),解耦矩陣可表示為

    由式(16)、式(17)可知,系統(tǒng)穩(wěn)態(tài)工作點(diǎn)方程是一個(gè)非線性方程組,求解過于復(fù)雜,且系統(tǒng)的解耦矩陣與當(dāng)前的穩(wěn)態(tài)工作點(diǎn)相關(guān),當(dāng)變換器負(fù)載發(fā)生變化或者控制變量參考值發(fā)生變化時(shí),系統(tǒng)的穩(wěn)態(tài)工作點(diǎn)也會發(fā)生變化,相應(yīng)的解耦矩陣H 也發(fā)生變化。因此,需要實(shí)時(shí)求解非線性方程組,進(jìn)一步增加了控制算法的復(fù)雜度。

    3 解耦參數(shù)求解與控制策略

    3.1 基于QPSO 的解耦參數(shù)求解

    傳統(tǒng)的求解非線性方程組的方法有蒙特卡洛法、牛頓迭代法、遺傳算法等。在實(shí)際求解過程中,蒙特卡洛法迭代次數(shù)過多,牛頓迭代法對初值過于敏感,遺傳算法求解精度較低。粒子群優(yōu)化PSO(particle swarm optimization)算法能夠更快地尋找到最優(yōu)解,且實(shí)現(xiàn)較為簡單,在實(shí)際中應(yīng)用較多,但是PSO 算法容易陷入局部最優(yōu)解。QPSO 算法是對PSO 算法的改進(jìn),該算法中的粒子具有量子行為,可以在整個(gè)解空間中進(jìn)行搜索,能夠克服PSO 算法的局限性[10]。

    QPSO 算法的粒子位置更新計(jì)算公式為

    式中:vmax為常數(shù),vi∈[-vmax,vmax];u 為服從[0,1]均勻分布的隨機(jī)數(shù);xi(t)為粒子當(dāng)前最優(yōu)位置;T 為最大迭代次數(shù);t 為當(dāng)前迭代次數(shù)為t 次迭代后最優(yōu)位置平均值,可表示為

    BR作為一種重要的植物激素,在植物生長發(fā)育的各個(gè)方面都發(fā)揮重要的作用。研究表明,BR可以調(diào)控?cái)M南芥種子的發(fā)育[12, 31],并提高種子的產(chǎn)量[5, 6]。油菜是與擬南芥親緣關(guān)系相近的一種重要油料作物,已有研究發(fā)現(xiàn),使用體外油菜素內(nèi)酯處理可以增加油菜單株籽粒的產(chǎn)量[23];外源施加eBL可以增加油菜幼苗根的細(xì)胞呼吸速率,增強(qiáng)其抗逆性[24, 25];甘藍(lán)型油菜過表達(dá)擬南芥的AtDWF4基因能夠提高油菜的種子產(chǎn)量,并增強(qiáng)油菜的抗逆性[26]。以上研究結(jié)果表明BR在油菜發(fā)育過程中也發(fā)揮重要作用。

    由于實(shí)際應(yīng)用中,控制器的運(yùn)算能力有限。本文采用QPSO 算法離線求解式(16)的非線性方程組,求得改進(jìn)穩(wěn)態(tài)工作點(diǎn)的電流值,得到i1、i3與φ12、φ13的對應(yīng)表格,存入控制系統(tǒng)中,通過查表法獲得當(dāng)前的穩(wěn)態(tài)工作點(diǎn),避免了實(shí)時(shí)運(yùn)算的壓力。算法求解流程如圖4 所示。

    圖4 QPSO 算法流程Fig.4 Flow chart of QPSO algorithm

    本文設(shè)定QPSO 算法優(yōu)化目標(biāo)函數(shù)為

    其中

    式中,l、m 為穩(wěn)態(tài)電流,根據(jù)實(shí)際控制情況,本文中其取值范圍為[0,3],求解時(shí)應(yīng)取其離散值,設(shè)定間隔為0.15,即l、m∈{0,0.15,0.30,0.45,…,2.55,2.70,2.85,3.00}。

    圖5 為l=2.25、m=1.50 時(shí),PSO 與QPSO 兩種算法的迭代結(jié)果對比,可知:PSO 算法在多次迭代過程中雖然收斂速度較快,但是計(jì)算結(jié)果容易出現(xiàn)差異,存在局部收斂;QPSO 算法雖然收斂速度較慢,但是能最大程度上克服局部收斂的情況。圖6為QPSO 算法迭代一次的誤差曲線,可知隨著迭代次數(shù)的增加誤差趨向于0,算法收斂。多次求解,最終可得到i1、i3與φ12、φ13的關(guān)系表格。為方便分析,假設(shè)移相控制器寄存器的最大值為1 125,將結(jié)果轉(zhuǎn)化為移相控制寄存器的值。表1 展示了部分表格值,表中(φ12,φ13)為對應(yīng)結(jié)果。

    表1 QPSO 算法求解結(jié)果(分辨率0.15)Tab.1 Solution results of QPSO algorithm(at resolution of 0.15)

    圖5 PSO 與QPSO 算法迭代結(jié)果Fig.5 Iterative results of PSO and QPSO algorithms

    圖6 QPSO 算法迭代誤差曲線Fig.6 Iterative error curve of QPSO algorithm

    3.2 解耦控制策略

    由式(20)可知,i1、i3的分辨率決定了表格數(shù)據(jù)的存儲量,數(shù)據(jù)量的多少與系統(tǒng)的控制效果直接相關(guān),數(shù)據(jù)太少不僅達(dá)不到解耦的目的,還可能出現(xiàn)系統(tǒng)不穩(wěn)定的情況;數(shù)據(jù)過多會消耗大量的系統(tǒng)內(nèi)存。

    雙線性插值法是一種二次插值法,通過在x 和y 兩個(gè)方向分別計(jì)算2 次得到目標(biāo)點(diǎn)的值,如圖7所示,設(shè)定P(x,y)為所求的目標(biāo)點(diǎn),4 個(gè)插值點(diǎn)為P(x1,y1)=Z1、P(x2,y1)=Z2、P(x1,y2)=Z3、P(x2,y2)=Z4,令a=x2-x1、b=y2-y1,則目標(biāo)點(diǎn)的值[11]為

    圖7 雙線性插值算法示意Fig.7 Schematic of bilinear interpolation algorithm

    該算法可以在i1、i3取低分辨率的情況下提高控制的精度,保證系統(tǒng)控制穩(wěn)定的同時(shí),降低表格的數(shù)據(jù)量,減少內(nèi)存的消耗。但是在過低的分辨率情況下,即使采用該算法,系統(tǒng)的控制精度依然會很低,系統(tǒng)穩(wěn)定性較差。因此,定義φ12的插值誤差為Δe12,其表達(dá)式為

    式中:φH(k)為i1、i3取較高分辨率時(shí)生成的表格值;φL(k)為i1、i3取較低分辨率時(shí)生成的表格值;N 為插值的次數(shù)。同理可得φ13的插值誤差Δe13。

    根據(jù)表1 提供的求解結(jié)果,為方便分析,令分辨率取0.15 時(shí)生成的為高分辨率表格,則分辨率取0.30、0.45、0.60 時(shí)得到的為低分辨率表格。對比高低分辨率表格,可通過式(23)插值得到低分辨率表格中所有未知的高分辨率表格中的值,即為φL(k);插值得到的φL(k)對應(yīng)在高分辨率表格中的值即為φH(k)。最后根據(jù)式(24)可得插值誤差。表2 展示了4 種分辨率的插值誤差與內(nèi)存消耗對比。

    由表2 可得,由于本文以分辨率取0.15 時(shí)生成的表格為基準(zhǔn),因此沒有插值誤差;當(dāng)分辨率為0.30 時(shí),插值誤差小于1%,而且內(nèi)存消耗相較于分辨率為0.15 時(shí)降低了73%,內(nèi)存使用大大減小。分辨率為0.45、0.60 時(shí),有相同結(jié)果。綜合考慮本文實(shí)際情況,以插值誤差小于1%,內(nèi)存消耗小于200 word 為界限,由以上分析可知分辨率為0.30 時(shí)符合條件。故本文最終采用的i1、i3分辨率為0.30。

    表2 不同分辨率插值誤差與內(nèi)存消耗對比Tab.2 Comparison of interpolation error and memory consumption at different resolutions

    系統(tǒng)的整體控制框圖如圖8 所示。實(shí)驗(yàn)中端口2 連接負(fù)載;端口1 接蓄電池,采用電壓電流雙閉環(huán)控制,用來穩(wěn)定端口2 的母線電壓;端口3 接超級電容,采用電流單環(huán)控制,并向母線傳輸能量??刂葡到y(tǒng)將端口1、3 的電流參考值輸入離線計(jì)算得到的表格中,并代入式(23)進(jìn)行雙線性插值運(yùn)算,得到系統(tǒng)解耦控制參數(shù),進(jìn)而實(shí)現(xiàn)解耦控制。

    圖8 系統(tǒng)控制框圖Fig.8 Control block diagram of system

    4 仿真與實(shí)驗(yàn)驗(yàn)證

    4.1 仿真分析

    為驗(yàn)證所提方案可行性,在Matlab/Simulink 中搭建了仿真電路。仿真參數(shù)為V1=150 V、V3=150 V、R=30 Ω;開關(guān)頻率40 kHz;L1、L2、L3分別為1.5 μH、5 μH、1.5 μH;變壓器匝數(shù)比為1∶2∶1。仿真過程中,端口1 的電壓外環(huán)給定電壓為300 V,端口3的電流環(huán)給定電流為15 A。在0.1 s 時(shí),端口2 的負(fù)載由30 Ω 突變?yōu)?5 Ω;在0.2 s 時(shí),負(fù)載由15 Ω突變?yōu)?0 Ω。控制算法分別采用未解耦控制、文獻(xiàn)[9]中的傳統(tǒng)解耦控制以及本文的改進(jìn)解耦控制。

    圖9 為TAB 變換器仿真結(jié)果。由圖9(a)可知,負(fù)載在突增突減過程中,端口1 的電流i1跟隨負(fù)載變化而變化,三種控制方式均能保持母線電壓的穩(wěn)定,達(dá)到閉環(huán)控制效果。從圖9(b)的局部放大圖可知,由于耦合作用,未解耦時(shí)i1、i3之間存在較大影響關(guān)系,端口輸出電流穩(wěn)態(tài)誤差的峰峰值較大,不能無誤差跟隨給定值,且負(fù)載突變時(shí)響應(yīng)速度較慢;采用解耦控制能夠改善系統(tǒng)控制效果,并且改進(jìn)解耦控制方法相較于傳統(tǒng)的解耦控制方法,系統(tǒng)負(fù)載突變時(shí)響應(yīng)速度更快,端口電流基本能夠無誤差的跟隨指定值,變換器的耦合關(guān)系進(jìn)一步減弱。

    圖9 解耦后負(fù)載突變仿真結(jié)果Fig.9 Simulation result of load mutation after decoupling

    4.2 實(shí)驗(yàn)驗(yàn)證

    為進(jìn)一步驗(yàn)證改進(jìn)解耦控制的實(shí)際控制效果,搭建了原理樣機(jī),如圖10 所示。實(shí)驗(yàn)所用控制器為德州儀器公司的TMS320F28069,電路參數(shù)如表3所示。實(shí)驗(yàn)過程中,系統(tǒng)從50%額定負(fù)載跳變到80%額定負(fù)載,再跳回50%額定負(fù)載。3 種控制算法與仿真相同。

    圖10 混合儲能三端口變換器硬件平臺Fig.10 Hardware platform of hybrid energy-storage three-port converter

    表3 電路參數(shù)Tab.3 Circuit parameters

    圖11 為TAB 電感電流波形。圖12~圖14 分別展示了3 種TAB 變換器控制方式的實(shí)驗(yàn)結(jié)果。從圖12~圖14 可知,實(shí)驗(yàn)結(jié)果與仿真基本一致。從圖12(b)~圖14(b)可以看出,3 種控制方法下,1、3 端口的電流穩(wěn)態(tài)誤差峰峰值分別約為400、250、150 mA;負(fù)載突增時(shí),電流的調(diào)節(jié)時(shí)間分別約為30、20、12 ms。因此,改進(jìn)解耦控制相對于傳統(tǒng)解耦控制,系統(tǒng)響應(yīng)速度更快,端口電流擾動更小,控制效果更好。

    圖11 TAB 電感電流波形Fig.11 Waveforms of TAB inductor current

    圖12 未解耦控制實(shí)驗(yàn)結(jié)果Fig.12 Experimental results under undecoupled control

    圖13 傳統(tǒng)解耦控制實(shí)驗(yàn)結(jié)果Fig.13 Experimental results under traditional decoupling control

    圖14 改進(jìn)解耦控制實(shí)驗(yàn)結(jié)果Fig.14 Experimental results under improved decoupling control

    5 結(jié)論

    為減小混合儲能三端口變換器之間的功率耦合,本文提出改進(jìn)解耦控制方法,將TAB 的軟開關(guān)特性和開關(guān)管損耗考慮到解耦網(wǎng)絡(luò)中,利用QPSO算法得到了解耦參數(shù)表格,采用雙線性插值算法,降低了系統(tǒng)的內(nèi)存消耗。仿真與實(shí)驗(yàn)結(jié)果表明:

    (1)與傳統(tǒng)解耦控制相比,改進(jìn)解耦控制的控制效果更好,端口間電流穩(wěn)態(tài)誤差更小。

    (2)在負(fù)載突增突減過程中,改進(jìn)解耦控制整體響應(yīng)更快。

    猜你喜歡
    插值穩(wěn)態(tài)表格
    可變速抽水蓄能機(jī)組穩(wěn)態(tài)運(yùn)行特性研究
    《現(xiàn)代臨床醫(yī)學(xué)》來稿表格要求
    碳化硅復(fù)合包殼穩(wěn)態(tài)應(yīng)力與失效概率分析
    電廠熱力系統(tǒng)穩(wěn)態(tài)仿真軟件開發(fā)
    煤氣與熱力(2021年4期)2021-06-09 06:16:54
    元中期歷史劇對社會穩(wěn)態(tài)的皈依與維護(hù)
    中華戲曲(2020年1期)2020-02-12 02:28:18
    統(tǒng)計(jì)表格的要求
    統(tǒng)計(jì)表格的要求
    統(tǒng)計(jì)表格的要求
    基于Sinc插值與相關(guān)譜的縱橫波速度比掃描方法
    一種改進(jìn)FFT多譜線插值諧波分析方法
    国产精品一区www在线观看| 下体分泌物呈黄色| 18+在线观看网站| 国产欧美另类精品又又久久亚洲欧美| 狠狠精品人妻久久久久久综合| 亚洲精品久久久久久婷婷小说| av国产久精品久网站免费入址| 在线播放无遮挡| 一级av片app| 91久久精品国产一区二区三区| 亚洲美女黄色视频免费看| 伦精品一区二区三区| 亚洲在久久综合| 网址你懂的国产日韩在线| kizo精华| 人妻少妇偷人精品九色| 免费高清在线观看视频在线观看| 干丝袜人妻中文字幕| 国产精品偷伦视频观看了| 国产精品一二三区在线看| 日本av免费视频播放| 特大巨黑吊av在线直播| 久久精品国产自在天天线| 亚洲av在线观看美女高潮| 亚洲伊人久久精品综合| 亚洲av在线观看美女高潮| 18禁裸乳无遮挡动漫免费视频| 亚洲人成网站高清观看| 久久人人爽人人爽人人片va| 日韩电影二区| 国产亚洲最大av| 国产伦理片在线播放av一区| 日韩成人av中文字幕在线观看| 亚洲性久久影院| 啦啦啦视频在线资源免费观看| h日本视频在线播放| 国产欧美日韩精品一区二区| 久久精品国产亚洲av天美| 一本色道久久久久久精品综合| 99热这里只有是精品50| 国产精品免费大片| av视频免费观看在线观看| 亚洲av电影在线观看一区二区三区| 伊人久久国产一区二区| 日日啪夜夜撸| 如何舔出高潮| 最近最新中文字幕免费大全7| 久久婷婷青草| 国产亚洲最大av| 不卡视频在线观看欧美| 欧美日韩国产mv在线观看视频 | 国产爱豆传媒在线观看| 韩国av在线不卡| 精品国产乱码久久久久久小说| 国产在线一区二区三区精| 久久人人爽人人爽人人片va| 亚洲精品一二三| 一区二区三区四区激情视频| 国产有黄有色有爽视频| 日韩 亚洲 欧美在线| 一区在线观看完整版| 精品亚洲成a人片在线观看 | 18禁动态无遮挡网站| 色吧在线观看| 欧美变态另类bdsm刘玥| 啦啦啦中文免费视频观看日本| 全区人妻精品视频| 一级黄片播放器| 菩萨蛮人人尽说江南好唐韦庄| 在线观看免费视频网站a站| 午夜免费鲁丝| 你懂的网址亚洲精品在线观看| 99热6这里只有精品| 国内揄拍国产精品人妻在线| 黄色欧美视频在线观看| 美女xxoo啪啪120秒动态图| 人人妻人人澡人人爽人人夜夜| 在线免费十八禁| 亚洲成人中文字幕在线播放| 又黄又爽又刺激的免费视频.| 美女xxoo啪啪120秒动态图| 熟女电影av网| 日韩一区二区三区影片| 菩萨蛮人人尽说江南好唐韦庄| 另类亚洲欧美激情| 六月丁香七月| 国产黄色免费在线视频| 又爽又黄a免费视频| 成人漫画全彩无遮挡| 人体艺术视频欧美日本| 久久人人爽人人片av| 日韩制服骚丝袜av| 男女无遮挡免费网站观看| 国产精品久久久久久久久免| 性色av一级| 大香蕉97超碰在线| 国精品久久久久久国模美| 麻豆国产97在线/欧美| 亚洲自偷自拍三级| 欧美最新免费一区二区三区| 国产69精品久久久久777片| 日韩三级伦理在线观看| 女性生殖器流出的白浆| 三级国产精品片| 中文精品一卡2卡3卡4更新| 国产精品一区二区三区四区免费观看| 黄色视频在线播放观看不卡| 久久99热这里只有精品18| 最近最新中文字幕大全电影3| 熟女av电影| 99精国产麻豆久久婷婷| 国内揄拍国产精品人妻在线| 美女xxoo啪啪120秒动态图| 久久久色成人| 亚洲不卡免费看| 久久毛片免费看一区二区三区| 亚洲成色77777| 2022亚洲国产成人精品| 少妇丰满av| 国产精品国产三级国产av玫瑰| 97在线视频观看| 亚洲成人手机| 亚洲国产成人一精品久久久| 十八禁网站网址无遮挡 | 国产成人freesex在线| 日韩人妻高清精品专区| 国产亚洲欧美精品永久| 日韩欧美精品免费久久| 亚洲精品中文字幕在线视频 | 午夜福利视频精品| 国产 精品1| 亚洲国产毛片av蜜桃av| 欧美激情极品国产一区二区三区 | 亚洲av男天堂| 啦啦啦视频在线资源免费观看| 国产乱来视频区| 久热这里只有精品99| 精品国产一区二区三区久久久樱花 | 青春草视频在线免费观看| 中文在线观看免费www的网站| 国产成人一区二区在线| 亚洲成色77777| 精品熟女少妇av免费看| 亚洲av国产av综合av卡| 搡女人真爽免费视频火全软件| av免费观看日本| 日本vs欧美在线观看视频 | 亚州av有码| 成人漫画全彩无遮挡| 亚洲欧美日韩另类电影网站 | 亚洲性久久影院| 岛国毛片在线播放| 一级毛片黄色毛片免费观看视频| 自拍偷自拍亚洲精品老妇| 色婷婷av一区二区三区视频| 多毛熟女@视频| 嫩草影院新地址| 久久99热这里只频精品6学生| 丝袜脚勾引网站| 国产 一区 欧美 日韩| 99国产精品免费福利视频| 我的老师免费观看完整版| 亚洲av.av天堂| 综合色丁香网| 久久精品国产亚洲网站| 国产欧美日韩精品一区二区| 伊人久久精品亚洲午夜| 交换朋友夫妻互换小说| 欧美日韩精品成人综合77777| 大码成人一级视频| 啦啦啦中文免费视频观看日本| 丝袜脚勾引网站| 国产欧美日韩一区二区三区在线 | 精品久久久久久电影网| 激情 狠狠 欧美| 国产 精品1| 国产成人精品一,二区| 国产av精品麻豆| 久久99热这里只有精品18| 最后的刺客免费高清国语| 99热这里只有精品一区| 亚洲欧美日韩东京热| 国产精品无大码| 美女xxoo啪啪120秒动态图| 偷拍熟女少妇极品色| 国产精品一区二区在线不卡| 一级毛片 在线播放| 国产亚洲91精品色在线| 欧美成人a在线观看| 尤物成人国产欧美一区二区三区| 内射极品少妇av片p| 国产精品国产三级专区第一集| 综合色丁香网| 久久久久久人妻| av天堂中文字幕网| 精品一区二区免费观看| 精品久久久久久久末码| 男男h啪啪无遮挡| 亚洲欧美一区二区三区国产| 美女脱内裤让男人舔精品视频| 国产亚洲最大av| 久久精品国产a三级三级三级| 日韩一本色道免费dvd| 少妇猛男粗大的猛烈进出视频| 免费人成在线观看视频色| 精华霜和精华液先用哪个| 国产精品熟女久久久久浪| 久久精品国产亚洲av涩爱| 国产一区亚洲一区在线观看| 亚洲av不卡在线观看| 制服丝袜香蕉在线| 一区二区三区免费毛片| 精品少妇黑人巨大在线播放| 九色成人免费人妻av| 特大巨黑吊av在线直播| 欧美亚洲 丝袜 人妻 在线| 午夜老司机福利剧场| 尾随美女入室| 国产国拍精品亚洲av在线观看| 国产一级毛片在线| 亚洲激情五月婷婷啪啪| av不卡在线播放| 成人18禁高潮啪啪吃奶动态图 | 久久人人爽av亚洲精品天堂 | 伦理电影免费视频| 日本欧美视频一区| 日韩制服骚丝袜av| 欧美日韩亚洲高清精品| 只有这里有精品99| 亚洲欧美中文字幕日韩二区| 国产精品久久久久久久电影| 91精品国产国语对白视频| av专区在线播放| 国产伦精品一区二区三区视频9| 成年人午夜在线观看视频| 伊人久久精品亚洲午夜| 久久久成人免费电影| 六月丁香七月| 国产精品不卡视频一区二区| 久久99精品国语久久久| 我要看黄色一级片免费的| 在线精品无人区一区二区三 | 国产伦在线观看视频一区| 亚洲精品色激情综合| 久久久久久久大尺度免费视频| 人体艺术视频欧美日本| 毛片一级片免费看久久久久| 国产伦在线观看视频一区| 国产男人的电影天堂91| 又爽又黄a免费视频| 国产精品福利在线免费观看| 亚洲欧美日韩另类电影网站 | av卡一久久| 久久久欧美国产精品| 一本—道久久a久久精品蜜桃钙片| 亚洲成人一二三区av| 中文天堂在线官网| 五月天丁香电影| 成人亚洲精品一区在线观看 | 精品人妻一区二区三区麻豆| 久久精品国产自在天天线| 少妇人妻精品综合一区二区| 国产女主播在线喷水免费视频网站| 91精品伊人久久大香线蕉| 亚洲欧美成人精品一区二区| 精品久久国产蜜桃| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品日韩av片在线观看| 国产视频内射| 日韩av在线免费看完整版不卡| 插阴视频在线观看视频| 大香蕉久久网| 久久99热这里只频精品6学生| 80岁老熟妇乱子伦牲交| 日本猛色少妇xxxxx猛交久久| 久久精品国产a三级三级三级| 亚洲美女黄色视频免费看| 国产人妻一区二区三区在| 美女脱内裤让男人舔精品视频| 久久精品久久精品一区二区三区| 久久av网站| 午夜激情久久久久久久| 免费黄频网站在线观看国产| 91aial.com中文字幕在线观看| 久久精品久久精品一区二区三区| 免费观看性生交大片5| 午夜免费鲁丝| 97在线人人人人妻| 亚洲精品乱码久久久久久按摩| 免费观看性生交大片5| 国产在线一区二区三区精| 国产高潮美女av| 偷拍熟女少妇极品色| 最近最新中文字幕大全电影3| 免费人妻精品一区二区三区视频| 天天躁日日操中文字幕| av在线老鸭窝| 久久99精品国语久久久| 亚洲高清免费不卡视频| 蜜桃在线观看..| 亚洲欧洲国产日韩| 在线免费观看不下载黄p国产| 久久国产精品大桥未久av | 亚洲av成人精品一二三区| 国产亚洲最大av| 美女视频免费永久观看网站| 女人久久www免费人成看片| 亚洲精品亚洲一区二区| 偷拍熟女少妇极品色| 九九爱精品视频在线观看| 亚洲经典国产精华液单| 国产亚洲91精品色在线| 一级黄片播放器| 联通29元200g的流量卡| 亚洲精品视频女| 啦啦啦视频在线资源免费观看| 国产成人a区在线观看| 亚洲精品自拍成人| 天天躁日日操中文字幕| 亚洲人成网站高清观看| 亚州av有码| 新久久久久国产一级毛片| 久久久久精品久久久久真实原创| 国产爽快片一区二区三区| 国产乱人视频| 国产高清有码在线观看视频| 秋霞在线观看毛片| 久久99热这里只频精品6学生| 在线免费观看不下载黄p国产| 久久久欧美国产精品| 九草在线视频观看| 亚洲怡红院男人天堂| 国产美女午夜福利| 国产精品一区二区三区四区免费观看| 爱豆传媒免费全集在线观看| 亚洲欧美一区二区三区国产| 看免费成人av毛片| 这个男人来自地球电影免费观看 | 九色成人免费人妻av| 高清av免费在线| 99热全是精品| 91精品一卡2卡3卡4卡| 视频区图区小说| a 毛片基地| 久久久久国产网址| 国产国拍精品亚洲av在线观看| 蜜桃亚洲精品一区二区三区| 日本黄大片高清| 日日摸夜夜添夜夜爱| 免费黄网站久久成人精品| 免费观看a级毛片全部| 国产无遮挡羞羞视频在线观看| 一区二区三区精品91| 国产 精品1| 久久久久久九九精品二区国产| av不卡在线播放| 欧美三级亚洲精品| 免费观看的影片在线观看| 日本免费在线观看一区| 2021少妇久久久久久久久久久| 你懂的网址亚洲精品在线观看| 国语对白做爰xxxⅹ性视频网站| av福利片在线观看| 日本一二三区视频观看| 1000部很黄的大片| 最近2019中文字幕mv第一页| 久久久午夜欧美精品| 韩国av在线不卡| 国产极品天堂在线| 日本黄色片子视频| 日韩中字成人| 国产在线一区二区三区精| 亚洲国产av新网站| 国内揄拍国产精品人妻在线| 蜜桃久久精品国产亚洲av| 精品国产三级普通话版| 国产深夜福利视频在线观看| 精品久久久久久电影网| 最近最新中文字幕大全电影3| 一本色道久久久久久精品综合| tube8黄色片| 我的女老师完整版在线观看| 又黄又爽又刺激的免费视频.| 国产一区亚洲一区在线观看| av专区在线播放| 久久久精品94久久精品| 男女无遮挡免费网站观看| 少妇被粗大猛烈的视频| 嘟嘟电影网在线观看| 欧美日韩视频精品一区| 精品人妻熟女av久视频| .国产精品久久| 欧美xxxx黑人xx丫x性爽| 99视频精品全部免费 在线| 一级毛片我不卡| 久久人妻熟女aⅴ| 毛片女人毛片| 两个人的视频大全免费| 黄色视频在线播放观看不卡| 一区二区三区精品91| 成人二区视频| 亚洲美女视频黄频| 在线观看一区二区三区激情| 国产精品免费大片| 国产伦理片在线播放av一区| 高清不卡的av网站| 色网站视频免费| 久久国内精品自在自线图片| 交换朋友夫妻互换小说| 青春草国产在线视频| 精品人妻视频免费看| 亚洲国产成人一精品久久久| 国产精品一及| 免费人成在线观看视频色| 久久久精品94久久精品| 日韩国内少妇激情av| 久久韩国三级中文字幕| 熟女电影av网| 亚洲欧美一区二区三区国产| 深夜a级毛片| 欧美日韩亚洲高清精品| av免费在线看不卡| 日韩一区二区视频免费看| 国产高清国产精品国产三级 | 天天躁日日操中文字幕| 亚洲欧美一区二区三区黑人 | 亚洲人成网站在线观看播放| 日本av手机在线免费观看| 看十八女毛片水多多多| 最近手机中文字幕大全| 日韩一本色道免费dvd| 日本黄色片子视频| 草草在线视频免费看| 一个人免费看片子| 欧美精品国产亚洲| 国产成人午夜福利电影在线观看| 国产av国产精品国产| 人人妻人人澡人人爽人人夜夜| 最近手机中文字幕大全| 久久精品人妻少妇| 国产精品蜜桃在线观看| 丰满人妻一区二区三区视频av| 免费人成在线观看视频色| 中国国产av一级| 欧美激情极品国产一区二区三区 | 亚洲精品乱码久久久久久按摩| 一本色道久久久久久精品综合| 乱码一卡2卡4卡精品| 大片电影免费在线观看免费| 91狼人影院| 美女福利国产在线 | 91久久精品国产一区二区成人| 精品国产露脸久久av麻豆| 午夜福利影视在线免费观看| 日韩一本色道免费dvd| 男人狂女人下面高潮的视频| 国产成人免费观看mmmm| 日韩欧美一区视频在线观看 | 欧美丝袜亚洲另类| 亚洲天堂av无毛| 99久久人妻综合| 国产亚洲91精品色在线| 国产成人一区二区在线| 一区二区三区免费毛片| 国产永久视频网站| 特大巨黑吊av在线直播| 国产欧美另类精品又又久久亚洲欧美| 97超碰精品成人国产| 亚洲,一卡二卡三卡| 美女脱内裤让男人舔精品视频| 下体分泌物呈黄色| av国产免费在线观看| 精品亚洲成国产av| 噜噜噜噜噜久久久久久91| 黄色一级大片看看| 日本欧美国产在线视频| 久久99蜜桃精品久久| 久久精品国产a三级三级三级| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美日韩在线观看h| 高清毛片免费看| 大片电影免费在线观看免费| 内射极品少妇av片p| 一级毛片我不卡| 国产乱人视频| 亚洲人成网站在线播| xxx大片免费视频| 日韩在线高清观看一区二区三区| 嫩草影院入口| 午夜福利高清视频| 国产欧美日韩一区二区三区在线 | 国产伦精品一区二区三区视频9| 日韩视频在线欧美| 精品少妇久久久久久888优播| 在线天堂最新版资源| 在线免费观看不下载黄p国产| 99热国产这里只有精品6| 国产日韩欧美在线精品| 少妇裸体淫交视频免费看高清| 成人亚洲欧美一区二区av| 亚洲欧美一区二区三区国产| 免费播放大片免费观看视频在线观看| 免费人成在线观看视频色| 日本爱情动作片www.在线观看| 最近中文字幕2019免费版| 亚洲欧美一区二区三区国产| 最新中文字幕久久久久| 日韩电影二区| av国产精品久久久久影院| 久久人妻熟女aⅴ| 韩国高清视频一区二区三区| 国产成人91sexporn| 成人午夜精彩视频在线观看| 欧美zozozo另类| 人妻夜夜爽99麻豆av| 亚洲成人一二三区av| 97超视频在线观看视频| 欧美+日韩+精品| 最近手机中文字幕大全| 国产精品成人在线| 亚洲av在线观看美女高潮| 欧美另类一区| av天堂中文字幕网| 少妇熟女欧美另类| 国产成人免费观看mmmm| 国产伦理片在线播放av一区| 春色校园在线视频观看| 久久99蜜桃精品久久| 在线观看人妻少妇| 青青草视频在线视频观看| 免费av不卡在线播放| 午夜激情福利司机影院| 国产精品av视频在线免费观看| a级一级毛片免费在线观看| 久久国产亚洲av麻豆专区| 我的老师免费观看完整版| 欧美 日韩 精品 国产| 在线观看三级黄色| 看免费成人av毛片| 乱系列少妇在线播放| 日韩不卡一区二区三区视频在线| a级一级毛片免费在线观看| 国产成人a∨麻豆精品| 日韩av不卡免费在线播放| 高清毛片免费看| 精品一区在线观看国产| 久久久久久久久久人人人人人人| 97精品久久久久久久久久精品| 日本wwww免费看| 国产成人午夜福利电影在线观看| 精品久久久精品久久久| 欧美精品一区二区免费开放| 免费少妇av软件| 久久久精品94久久精品| 熟女av电影| 一区二区三区四区激情视频| h日本视频在线播放| 大香蕉97超碰在线| 中文字幕av成人在线电影| 久久久久性生活片| 在线观看一区二区三区| 纯流量卡能插随身wifi吗| 久久久久视频综合| 看免费成人av毛片| 亚洲欧美一区二区三区国产| 少妇人妻久久综合中文| 亚洲精品一二三| av国产久精品久网站免费入址| 一区二区三区精品91| 亚洲av成人精品一区久久| 久久青草综合色| 人妻夜夜爽99麻豆av| 久久韩国三级中文字幕| 亚洲av中文av极速乱| 97在线人人人人妻| 少妇裸体淫交视频免费看高清| 人妻夜夜爽99麻豆av| 内地一区二区视频在线| 国产精品爽爽va在线观看网站| 在线天堂最新版资源| 国产v大片淫在线免费观看| 国产精品三级大全| 最近中文字幕2019免费版| 亚洲av综合色区一区| 亚洲内射少妇av| av福利片在线观看| 国产无遮挡羞羞视频在线观看| 午夜激情久久久久久久| 一个人看的www免费观看视频| 麻豆精品久久久久久蜜桃| 亚洲国产欧美在线一区| 久久久久久久久久成人| 下体分泌物呈黄色| 99热网站在线观看| 国产探花极品一区二区| 狂野欧美激情性bbbbbb| 一本色道久久久久久精品综合| 久久久久国产网址| 少妇熟女欧美另类| 在线观看三级黄色| 久久女婷五月综合色啪小说| 国产视频内射| 国产亚洲欧美精品永久| 成人亚洲欧美一区二区av| 免费看日本二区| 一级毛片电影观看| 免费观看在线日韩| 亚洲人成网站在线播| 18禁在线播放成人免费| 国产高潮美女av| 国产日韩欧美亚洲二区| 久久99热这里只频精品6学生| 亚洲精品乱码久久久久久按摩| 国产精品麻豆人妻色哟哟久久| 日韩成人av中文字幕在线观看|