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

    基于NSGA-Ⅱ算法的ECPT系統(tǒng)PID參數(shù)尋優(yōu)及輸出穩(wěn)壓控制

    2016-11-16 00:44:41蘇玉剛陳苓芷唐春森馬浚豪呼愛國
    電工技術(shù)學報 2016年19期
    關(guān)鍵詞:控制參數(shù)耦合控制器

    蘇玉剛 陳苓芷 唐春森 馬浚豪 呼愛國

    (1.重慶大學自動化學院 重慶 400030 2.奧克蘭大學電子與計算機工程系 奧克蘭 1010)

    ?

    基于NSGA-Ⅱ算法的ECPT系統(tǒng)PID參數(shù)尋優(yōu)及輸出穩(wěn)壓控制

    蘇玉剛1陳苓芷1唐春森1馬浚豪1呼愛國2

    (1.重慶大學自動化學院重慶400030 2.奧克蘭大學電子與計算機工程系奧克蘭1010)

    針對電場耦合型無線電能傳輸(ECPT)系統(tǒng)階數(shù)高、非線性、對系統(tǒng)參數(shù)變化敏感導致其輸出電壓質(zhì)量下降、系統(tǒng)不穩(wěn)定等問題,提出一種最優(yōu)控制器設計方法。該方法首先建立系統(tǒng)廣義狀態(tài)空間平均(GSSA)模型,并基于此模型,利用多目標多約束遺傳算法(NSGA-II)對PID控制參數(shù)進行自動尋優(yōu),有效解決了高階系統(tǒng)PID控制器最優(yōu)參數(shù)難以設計的難題,改善了閉環(huán)系統(tǒng)的上升時間、穩(wěn)態(tài)誤差以及對變化參數(shù)的魯棒性能,并將超調(diào)限定在一定范圍內(nèi),提升系統(tǒng)的穩(wěn)定性。仿真和實驗結(jié)果驗證了GSSA模型的準確性以及在此模型上利用NSGA-II算法進行控制參數(shù)尋優(yōu)的有效性。

    無線電能傳輸電場耦合參數(shù)尋優(yōu)廣義狀態(tài)空間平均模型PID控制

    0 引言

    經(jīng)過多年的發(fā)展,無線電能傳輸(Wireless Power Transfer,WPT)技術(shù)在理論和工程應用方面取得了一系列成果[1-5]。該新型電能傳輸方式可以使電能的傳輸實現(xiàn)從供體到受體的無直接電氣連接,從而消除了傳統(tǒng)接觸式電能傳輸方式存在的布線雜亂、導線磨損、接觸火花等問題。因此,無線電能傳輸技術(shù)尤為適用于安全性要求苛刻且工作環(huán)境復雜的用電場合(如潮濕、易燃、易爆等工作環(huán)境[6,7])。在無線電能傳輸技術(shù)研究領域,基于電磁感應耦合的無線電能傳輸(Inductively Coupled Power Transfer,ICPT)技術(shù)是目前較熱的研究方向,在理論研究與實際應用中都取得了較大的突破且逐步推廣應用。由于電場與磁場具有某些相似特性以及理論分析的對偶性,且基于電場耦合的無線電能傳輸方式具有一些獨有的優(yōu)勢(比如當耦合機構(gòu)之間或周圍存在金屬導體時,不會使導體產(chǎn)生渦流損耗等[8]),因此,基于電場耦合的無線電能傳輸(Electric-field Coupled Power Transfer,ECPT )技術(shù)逐漸被國內(nèi)外相關(guān)團隊與機構(gòu)所重視[9]。近年來,國內(nèi)外相關(guān)學者對ECPT技術(shù)的某些特性進行了研究。文獻[10]對比分析了ECPT系統(tǒng)盤式以及圓柱式耦合機構(gòu)的特性,對基于圓柱式耦合機構(gòu)的ECPT系統(tǒng)的工作頻率、最大傳輸功率以及功率損耗進行了分析。文獻[11]利用頻閃映射建模法建立了ECPT系統(tǒng)模型,并基于此模型推出系統(tǒng)ZVS頻率點以及各交流量的穩(wěn)態(tài)波形,分析系統(tǒng)參數(shù)改變對ZVS頻率點的影響。文獻[12]提出了一種利用電容矩陣方式實現(xiàn)調(diào)諧控制的方法,實現(xiàn)了系統(tǒng)諧振頻率的穩(wěn)定,優(yōu)化了系統(tǒng)能效積。

    在ECPT系統(tǒng)實際應用中,電能發(fā)射極板與電能拾取極板的相對位置經(jīng)常改變。當電場耦合機構(gòu)兩極板的相對位置發(fā)生變化時,電容耦合系數(shù)亦將隨之變化,以至于ECPT系統(tǒng)能量拾取側(cè)輸出電壓的質(zhì)量因動態(tài)波動而下降[8]。此外,在變負載或多負載的ECPT系統(tǒng)中,負載的變化也將很大程度上造成ECPT系統(tǒng)輸出電壓的波動[13]。為了提高ECPT系統(tǒng)輸出電壓的穩(wěn)定性,可靠且實用的輸出穩(wěn)壓控制策略是非常有必要。

    目前,關(guān)于ECPT系統(tǒng)輸出穩(wěn)壓控制方面的研究文獻相對較少,文獻[8]引入極板陣列和可變諧振電感(變壓器二次側(cè)并聯(lián)兩個開關(guān)管),當耦合機構(gòu)的拾取側(cè)極板位置發(fā)生偏移時,通過改變極板陣列的供電極板和可變諧振電感值來保持輸出電壓穩(wěn)定在原有值,但它并不能控制其穩(wěn)定在任意設定值,且沒有考慮負載切換的情況。而在ICPT系統(tǒng)輸出控制研究方面,有諸多文獻可查[14-17],如文獻[14]建立了雙向ICPT系統(tǒng)的動態(tài)模型,并在此模型的基礎上基于多目標遺傳算法設計了一組最優(yōu)PID控制參數(shù),優(yōu)化系統(tǒng)在穩(wěn)壓控制過程的動態(tài)性能。雖然在基本理論分析方面,ICPT技術(shù)與ECPT技術(shù)具有一定的相似性和對偶性,但由于ECPT系統(tǒng)具有本身工作頻率高,對系統(tǒng)參數(shù)變化敏感等特點,使得對其進行輸出控制策略研究變得困難。在對ECPT系統(tǒng)輸出穩(wěn)壓控制研究時,可借鑒ICPT系統(tǒng)已有的研究成果,并結(jié)合ECPT系統(tǒng)自身的特點,研究適合其特性的輸出控制方法。

    本文建立了半橋式ECPT系統(tǒng)的廣義狀態(tài)空間平均(Generalized State Space Averaging,GSSA)模型,并針對ECPT系統(tǒng)結(jié)構(gòu)復雜、階數(shù)較高、對系統(tǒng)參數(shù)變化敏感等問題,基于系統(tǒng)GSSA模型,利用多目標多約束遺傳算法NSGA-Ⅱ自動搜尋最優(yōu)PID控制參數(shù),從而實現(xiàn)系統(tǒng)的輸出穩(wěn)壓控制,保證系統(tǒng)輸出響應快、穩(wěn)態(tài)誤差小,并且當系統(tǒng)參數(shù)在一定范圍內(nèi)變化時其輸出仍能有效跟蹤期望值。由于該方法通過控制系統(tǒng)發(fā)射側(cè)逆變電路開關(guān)管的導通與關(guān)斷來實現(xiàn)輸出電壓穩(wěn)定,故系統(tǒng)無需增加額外的能量調(diào)節(jié)電路,降低了控制電路的復雜性。

    1 半橋式ECPT系統(tǒng)拓撲結(jié)構(gòu)及原理

    半橋式逆變電路因結(jié)構(gòu)簡單、可控性好、工作頻率高等優(yōu)點,在ECPT系統(tǒng)工程中得到較廣泛的應用[9],故本文選用半橋式電路拓撲作為系統(tǒng)的逆變環(huán)節(jié)進行研究。半橋式ECPT系統(tǒng)拓撲結(jié)構(gòu)如圖1所示,該系統(tǒng)由發(fā)射側(cè)和拾取側(cè)組成。

    圖1 基于半橋式ECPT系統(tǒng)電路拓撲Fig.1 The topology of the ECPT system based on half-bridge

    圖1中,Cs1、Cs2為ECPT系統(tǒng)中耦合極板的等效電容,RL為負載等效電阻。該ECPT系統(tǒng)的直流輸入經(jīng)過半橋逆變環(huán)節(jié)和發(fā)射側(cè)諧振網(wǎng)絡后得到高頻交流電,在其作用下,耦合機構(gòu)的兩極板之間形成交互電場,在交互電場的作用下產(chǎn)生位移電流“流過”極板,實現(xiàn)電能的無線傳輸。接收極板接收到的電能經(jīng)過拾取側(cè)諧振網(wǎng)絡與整流環(huán)節(jié)變換后提供給用電設備。若要對系統(tǒng)輸出電壓進行閉環(huán)控制,可以通過能量信號并行傳輸或射頻傳輸?shù)确绞綄⑤敵鲭妷翰蓸有盘枏氖叭?cè)傳輸?shù)桨l(fā)射側(cè)[18,19],采樣獲得的實際值與期望設定值比較后得到偏差量,將此偏差量作為控制器輸入,再根據(jù)具體控制策略生成對應的控制量來控制逆變電路開關(guān)管的導通與關(guān)斷,由此構(gòu)成閉環(huán)控制系統(tǒng),以保證輸出電壓的穩(wěn)定性。

    從拓撲結(jié)構(gòu)上來看,該系統(tǒng)主要由半橋逆變環(huán)節(jié)、阻抗匹配環(huán)節(jié)、單級諧振網(wǎng)絡、LCL復合諧振網(wǎng)絡、整流濾波環(huán)節(jié)以及負載組成。其中,單級諧振網(wǎng)絡中的電感L2a用于對耦合極板的等效電容Cs1、Cs2進行補償;L2b、C2、L3在系統(tǒng)拾取側(cè)構(gòu)成LCL復合諧振網(wǎng)絡,其作用主要有兩點:一是提高A、B端右側(cè)網(wǎng)絡的功率因數(shù);二是在耦合極板電流較小的條件下,提升了流經(jīng)負載RL的輸出電流,且由于耦合極板的電壓差較小,系統(tǒng)的運行安全性得到提高。電容C1b用于補償其右側(cè)網(wǎng)絡的阻抗使其呈現(xiàn)純阻性;由于LCL網(wǎng)絡的引入使系統(tǒng)拾取側(cè)的總阻抗增加,因此在逆變電路輸出端利用電感L1及電容C1a對其右側(cè)所有電路的等效阻抗進行匹配,提升系統(tǒng)輸出功率。

    2 半橋式ECPT系統(tǒng)的廣義狀態(tài)空間平均模型

    ECPT系統(tǒng)因階數(shù)高、開關(guān)非線性等特性使得其對參數(shù)變化較為敏感,而且因為結(jié)構(gòu)復雜,也使其理論分析變得困難。因此,有必要建立一個精確的數(shù)學模型對其內(nèi)在動力學行為進行描述與分析,并在模型的基礎上設計一定的控制律來保證系統(tǒng)輸出的穩(wěn)定性。本文采用廣義狀態(tài)空間平均法對半橋式ECPT系統(tǒng)進行建模,它利用傅里葉變換思想將時域非線性系統(tǒng)轉(zhuǎn)換為頻域線性模型,由此分析其暫態(tài)和穩(wěn)態(tài)過程[17]。

    將圖1所示的系統(tǒng)電路拓撲進行簡化。將電容C1a、C1b簡化為電容C1,即C1=C1a+C1b;將電感L2a、L2b簡化為電感L2,即L2=L2a+L2b;將耦合極板等效電容Cs1、Cs2簡化為Cs,即Cs=Cs1Cs2/(Cs1+Cs2)。并分別定義RL1、RL2、RL3為簡化系統(tǒng)中電感L1、L2、L3的內(nèi)阻值。選取簡化系統(tǒng)中電感電流和電容電壓作為狀態(tài)變量。分別采用零次諧波分量和基波分量近似表示ECPT系統(tǒng)中的直流狀態(tài)變量uCf和交流狀態(tài)變量iL1、uC1、uCs,iL2、uC2、iL3。因此,定義

    (1)

    利用傅里葉變換將ECPT系統(tǒng)時域模型轉(zhuǎn)換為頻域的線性微分模型為

    (2)

    式中,g(t)、p(t)為非線性轉(zhuǎn)換函數(shù),分別表示發(fā)射側(cè)半橋式逆變環(huán)節(jié)和拾取側(cè)整流環(huán)節(jié),函數(shù)的定義為

    (3)

    (4)

    式中,T為逆變器工作周期;θ1、θ2為函數(shù)的初相角;n為周期的整數(shù)倍。展開式(2)中各變量的傅里葉系數(shù)實虛部,依次定義為系統(tǒng)的廣義狀態(tài)變量

    x(t)=[Re〈iL1〉1,Im〈iL1〉1,Re〈uC1〉1,Im〈uC1〉1,

    Re〈uCs〉1,Im〈uCs〉1,Re〈iL2〉1,Im〈iL2〉1,

    Re〈uC2〉1,Im〈uC2〉1,Re〈iL3〉1,Im〈iL3〉1,

    〈uCf〉0]T

    (5)

    設系統(tǒng)的直流輸入Edc為系統(tǒng)控制量,基于式(2)的頻域線性微分方程描述,可建立以系統(tǒng)傅里葉系數(shù)實虛部為變量的廣義狀態(tài)空間平均模型

    (6)

    式中,A為系統(tǒng)矩陣;B為輸入矩陣。通過將時域狀態(tài)變量進行頻域分解,將ECPT系統(tǒng)的時域非線性微分模型轉(zhuǎn)換為頻域線性微分模型,實現(xiàn)非線性系統(tǒng)的線性化。

    3 基于NSGA-Ⅱ算法的PID控制器設計

    在實際工程應用中,ECPT系統(tǒng)輸出電壓的波動將會在一定程度上影響到整個系統(tǒng)的穩(wěn)定性,因此,有必要采取可靠的輸出穩(wěn)壓控制策略,對系統(tǒng)輸出進行實時控制。為了使系統(tǒng)輸出電壓穩(wěn)定,本文采用PID控制器對系統(tǒng)進行輸出控制,閉環(huán)框圖如圖2所示。將GSSA模型的輸出值與設定值相比,經(jīng)過PID控制模塊后得到控制量,將此控制量作為GSSA模型的輸入,由此保證系統(tǒng)輸出電壓穩(wěn)定。

    圖2 系統(tǒng)輸出控制閉環(huán)結(jié)構(gòu)Fig.2 Closed-loop block structure for output voltage regulation

    PID控制器的傳遞函數(shù)一般為

    (7)

    式中,Kp為比例系數(shù);Ti為積分時間常數(shù);Td為微分時間常數(shù)。

    在PID控制器中,Kp、Ti、Td三個參數(shù)的選擇直接影響PID控制器的性能,因此PID控制器設計的關(guān)鍵問題是如何選擇這三個參數(shù)。到目前為止,有許多PID參數(shù)設計方法,其中較常見的是ZN方法[20],然而這些參數(shù)設計方法通常只適用于低階對象,對于ECPT系統(tǒng)這樣高階復雜的系統(tǒng),想要設計出一組能滿足多個控制目標的PID參數(shù),顯然這些方法是不太適用的。針對PID參數(shù)整定問題的多目標性,本文選擇多目標多約束遺傳算法NSGA-Ⅱ?qū)刂破鲄?shù)進行全局優(yōu)化,自動搜尋出一系列非支配最優(yōu)解,并從中找出一組合適的控制參數(shù),使系統(tǒng)在該組控制參數(shù)的作用下,當系統(tǒng)可變參數(shù)發(fā)生變化時其輸出仍能達到預期控制目標(響應快、誤差小、超調(diào)小、對變化參數(shù)魯棒性好)。

    3.1NSGA-Ⅱ算法原理

    圖3 NSGA-Ⅱ算法流程Fig.3 Flowchart of the NSGA-Ⅱ algorithm

    NSGA-Ⅱ算法是由K.Deb提出的一種在多目標遺傳算法界非常有效,并基于(約束性)非支配排序、精英保留策略與多樣性維持機制的優(yōu)化算法[21]。該算法所特有的精英保留策略和多樣性維持機制可以確保其算法的收斂性與多樣性。具體過程如圖3所示:①首先根據(jù)經(jīng)驗估計出PID三個控制參數(shù)的取值范圍,然后在這個范圍內(nèi)隨機生成大小為N的初始種群;②根據(jù)目標函數(shù)和約束函數(shù)分別計算出每個個體所對應的目標函數(shù)值與約束函數(shù)值(即指標計算);③根據(jù)個體的指標值對種群個體進行(約束性)非支配排序分層,并分別計算每個個體的擁擠度距離;④通過選擇(擁擠度錦標賽選擇)、交叉(模擬二元交叉)、變異(多項式變異)操作后得到子代種群;⑤將子代種群與父輩種群合并形成大小為2N的臨時種群;⑥對臨時種群的每個個體進行指標計算;⑦對臨時種群的每個個體進行(約束性)非支配排序分層,并分別計算其擁擠度距離;⑧根據(jù)每個個體與第一非支配前沿的接近程度及其擁擠度距離大小,從臨時種群中選擇前N個個體組成新的種群,作為下一代遺傳操作的父輩,這個過程一直重復到算法滿足結(jié)束要求為止。

    3.2優(yōu)化目標和約束條件

    3.2.1優(yōu)化目標

    所謂PID參數(shù)優(yōu)化,實際就是利用算法來優(yōu)化Kp、Ti、Td三個參數(shù),其本質(zhì)是基于一定的目標函數(shù)和約束條件的參數(shù)尋優(yōu)問題。因此,在遺傳算法中關(guān)鍵的一步就是確定目標函數(shù)與約束條件,它決定了一個染色體(一組PID參數(shù))的好壞。本文選取以下三個目標函數(shù)作為算法的性能評價指標。

    1)上升時間trise定義為系統(tǒng)輸出從穩(wěn)態(tài)值的10%到穩(wěn)態(tài)值的90%所需的時間,即

    J1=t90%-t10%

    (8)

    2)時間乘平方誤差積分ess。

    (9)

    J3=J3,1+J3,2+J3,3+J3,4

    (10)

    式中

    (11)

    式中,J3,1、J3,2、J3,3和J3,4分別表示當可變參數(shù)(RL,Cs)在以下4種切換情況(從標稱值到極端最壞值)下變化時,系統(tǒng)輸出從切換時刻起始的時間乘平方誤差積分,其中,t表示切換時刻。

    3.2.2約束條件

    由于NSGA-II算法是以搜尋整個非支配前沿為目的來進行優(yōu)化的,而非支配前沿的某些邊界極端解(某個目標函數(shù)值非常大,但另一個目標函數(shù)值非常小)實用意義通常不大。因此,相較于第一非支配前沿的某些邊界極端解,設計者更期望第二(第三,等)非支配前沿的中間解能夠保留下來。以兩個目標函數(shù)的最小值優(yōu)化問題為例,對非支配最優(yōu)解的選擇情況分析圖如圖4所示。

    圖4 非支配最優(yōu)解分析Fig.4 Analysis of non-dominated solutions

    從圖4可以看出,優(yōu)化解1、2、3、4屬于第一非支配前沿,優(yōu)化解5、6、7屬于第二非支配前沿,設計者期望的目標函數(shù)值邊界如虛線框所示。顯然,優(yōu)化解1、4、5、7都不是設計者所期望的解,所以,當?shù)谝环侵淝把氐倪吔鐦O端解(如解1)與第二非支配前沿的中間解(如解6)相比較時,設計者更期望將后者(解6)作為較優(yōu)解保留下來。因此,算法需要對非支配最優(yōu)解的目標函數(shù)值作約束處理,若不對其進行約束的話,那些實用意義不大的邊界極端解將會一直被視作最優(yōu)解,并代代遺傳下去,造成最優(yōu)解質(zhì)量下降和計算資源的浪費。為了提高NSGA-II算法的計算效率與非支配最優(yōu)解的實用性,本文引入約束條件對其進行優(yōu)化處理,旨在將每一代的非支配最優(yōu)解都限制到其所有目標函數(shù)值均接近各自平均值的那些非支配優(yōu)化解。

    因此,本文選取以下4個經(jīng)驗條件作為算法的約束評價指標

    Cons1: trise≤0.006

    (12)

    Cons2: ess≤0.000 3

    (13)

    Cons3: err≤0.1

    (14)

    (15)

    式中,Cons1、Cons2、Cons3分別對三個目標函數(shù)值進行邊界條件約束;Cons4對系統(tǒng)輸出在響應過程中的超調(diào)量進行約束。

    在基于NSGA-Ⅱ算法的ECPT系統(tǒng)的PID控制參數(shù)尋優(yōu)問題中,對上述目標函數(shù)值與約束函數(shù)值的計算即圖3中的指標計算部分,它是根據(jù)ECPT系統(tǒng)的GSSA模型來進行計算的。具體的指標計算流程如圖5所示。

    圖5 指標計算流程Fig.5 Flowchart of indexes evaluation of objectives and constraints

    3.3優(yōu)化結(jié)果

    應用帶約束的NSGA-Ⅱ算法對基于GSSA模型的PID控制參數(shù)優(yōu)化模型進行參數(shù)尋優(yōu),具體遺傳算法的參數(shù)設置為:種群大小100,遺傳代數(shù)30,交叉概率0.9,變異概率0.1。圖6為當遺傳代數(shù)為30時100個種群個體的目標值分布情況。

    由圖6可看出,當種群進化到30代時,其非支配最優(yōu)解均勻分布在最優(yōu)解集中,種群中每個個體的三個目標值均達到較小值,控制性能效果較好。在100個種群個體中選擇等級為1,擁擠度距離較大的一個非支配最優(yōu)解作為最終選定的最優(yōu)解Sopt

    Kp=0.887,Ti=0.002、Td=0

    即選用(0.887,0.002,0)作為此次PID控制參數(shù)優(yōu)化的最終解。

    圖6 NSGA-Ⅱ算法優(yōu)化結(jié)果Fig.6 Optimization result of NSGA-Ⅱalgorithm

    4 仿真結(jié)果與實驗分析

    為了驗證基于NSGA-Ⅱ算法優(yōu)化的PID控制器對ECPT系統(tǒng)輸出控制效果,本文基于上述所提出的PID控制器參數(shù)尋優(yōu)方法及目標函數(shù)和約束條件設計出一組最優(yōu)PID控制參數(shù)(0.887,0.002,0),將其帶入閉環(huán)控制系統(tǒng)進行仿真分析和實驗驗證。在圖1所示的系統(tǒng)電路拓撲基礎上,加入PID控制回路,根據(jù)系統(tǒng)期望輸出電壓值與采樣得到的實際值比較后獲得偏差及偏差變化率,經(jīng)過PID控制器模塊運算后得到相應的控制量,占空比計算模塊將其轉(zhuǎn)換為逆變電路開關(guān)管的控制信號,從而構(gòu)成閉環(huán)控制系統(tǒng),保證輸出電壓穩(wěn)壓控制。

    令輸出電壓設定值為10 V,在Matlab/Simulink平臺下對閉環(huán)系統(tǒng)進行仿真分析,系統(tǒng)電路參數(shù)取值見表1。

    表1 系統(tǒng)參數(shù)取值Tab.1 System parameters

    在仿真運行中利用開關(guān)器件同時切換可變參數(shù)RL、Cs,使負載RL變化±20%,Cs變化±10%,分別模擬4種參數(shù)變化的最壞情況。負載兩端電壓u0的仿真波形如圖7所示。

    圖7 ECPT系統(tǒng)輸出電壓仿真結(jié)果Fig.7 Simulation results of the ECPT system

    由圖7可以看出,當系統(tǒng)耦合機構(gòu)等效電容以及負載電阻在允許變化范圍內(nèi)發(fā)生4種最壞情況變化時,閉環(huán)系統(tǒng)輸出經(jīng)過最快8 ms、最長20 ms左右的調(diào)節(jié)過程后均能穩(wěn)定下來。而在系統(tǒng)起動過程,輸出電壓經(jīng)過8 ms左右從穩(wěn)態(tài)值的10%上升到穩(wěn)態(tài)值的90%,與圖6中的優(yōu)化目標值一致。

    基于同樣的控制條件,現(xiàn)利用傳統(tǒng)ZN參數(shù)整定方法對控制參數(shù)進行整定。由于ZN算法是一種處于經(jīng)驗和計算之間的中間方法,它只是以系統(tǒng)輸出瞬態(tài)響應特性為優(yōu)化目標,建立一套經(jīng)驗公式以獲得一組瞬態(tài)響應較優(yōu)的控制參數(shù),而難以在多目標優(yōu)化條件下整定出一組控制參數(shù)使閉環(huán)系統(tǒng)實現(xiàn)多個目標值最優(yōu)的效果。為了比較兩者的瞬態(tài)響應特性,將ZN算法整定出的控制參數(shù)(3.3,1.52×10-3,3.96×10-3)帶入閉環(huán)系統(tǒng)中,得到響應結(jié)果見表2。從表中數(shù)據(jù)可以看出,閉環(huán)系統(tǒng)在ZN參數(shù)整定方法得到的控制器的作用下上升時間優(yōu)于NSGA-Ⅱ算法,但其穩(wěn)定時間(輸出從0上升到穩(wěn)態(tài)值)明顯劣于NSGA-Ⅱ算法。當系統(tǒng)可變參數(shù)RL、Cs發(fā)生變化時,利用ZN算法的閉環(huán)系統(tǒng)能夠穩(wěn)定下來,但調(diào)節(jié)時間最短為12 ms、最長為22 ms,且每種跳變情況的調(diào)節(jié)時間均劣于或等于NSGA-Ⅱ算法。由此可見,經(jīng)過NSGA-Ⅱ算法優(yōu)化后的PID控制器的控制效果在動態(tài)性能以及魯棒性能方面優(yōu)于傳統(tǒng)的ZN算法,且NSGA-Ⅱ算法是基于多個優(yōu)化目標來對參數(shù)進行整定的,它具有很強的延伸性,可以通過調(diào)整不同的優(yōu)化目標實現(xiàn)不同的控制效果,這是傳統(tǒng)參數(shù)優(yōu)化算法不能比擬的。

    表2 PID控制器輸出控制效果比較Tab.2 Comparison of different PID controllers

    基于圖1的拓撲結(jié)構(gòu)和表1的系統(tǒng)參數(shù),搭建了半橋式ECPT系統(tǒng)實驗平臺。本文選用STM32F407單片機作為系統(tǒng)的控制器實現(xiàn)PID控制器運算與占空比計算。單片機輸出逆變電路開關(guān)管的控制信號,驅(qū)動電路將其轉(zhuǎn)換為實際的驅(qū)動脈沖。系統(tǒng)起動過程的負載輸出電壓u0、Q1管柵極觸發(fā)脈沖uGS、逆變器輸出電壓u1實驗波形如圖8所示。從圖中可看出,Q1管柵極觸發(fā)脈沖uGS作為閉環(huán)系統(tǒng)的直接控制信號,當系統(tǒng)開始運行時便根據(jù)PID控制器輸出的控制量轉(zhuǎn)換為占空比可變的觸發(fā)脈沖信號。系統(tǒng)在半橋逆變環(huán)節(jié)的作用下,逆變器輸出電壓u1為占空比可變的方波輸出,其占空比與觸發(fā)脈沖uGS的占空比一致。與仿真結(jié)果類似,PID控制器經(jīng)過約12 ms的動態(tài)調(diào)節(jié),使系統(tǒng)輸出電壓u0從1 V上升到9 V,且進入穩(wěn)態(tài)后,其穩(wěn)態(tài)誤差基本為零。在整個動態(tài)過程中系統(tǒng)超調(diào)很小,調(diào)節(jié)時間與仿真基本一致。由此說明,在實際應用中,該PID控制器使系統(tǒng)有較好的動態(tài)特性和靜態(tài)特性。

    圖8 ECPT系統(tǒng)在起動過程的實驗波形Fig.8 Experimental waveform of the ECPT system with step response

    當系統(tǒng)可變參數(shù)(Cs,RL)發(fā)生變化時,系統(tǒng)輸出電壓u0與流經(jīng)L3的電流iL3的實驗波形如圖9所示。實驗結(jié)果與圖7中的仿真結(jié)果基本一致,當系統(tǒng)參數(shù)發(fā)生4種最壞情況跳變時,輸出電壓u0均能經(jīng)過一個較短的調(diào)節(jié)過程后進入10 V穩(wěn)態(tài)階段,調(diào)節(jié)時間最快為12 ms、最長為30 ms左右。當負載阻值發(fā)生變化時,流經(jīng)L3的電流iL3會發(fā)生相應變化來保證系統(tǒng)輸出電壓值恒定。當負載阻值增大時,iL3的幅值會相應減小,反之,當負載阻值減小時,iL3的幅值會相應增大以保證u0值不變。

    圖9 ECPT系統(tǒng)在參數(shù)跳變下的實驗波形Fig.9 Transient response of the ECPT system when parameter variations

    5 結(jié)論

    本文采用PID控制器對半橋式ECPT系統(tǒng)進行輸出穩(wěn)壓控制,建立了系統(tǒng)的廣義狀態(tài)空間平均模型,并基于此模型提出了一種利用多目標多約束遺傳算法對控制參數(shù)進行自動尋優(yōu)的方法。該優(yōu)化方法在多目標多約束條件下自動搜尋出一組最優(yōu)PID控制參數(shù),保證閉環(huán)系統(tǒng)在該控制器作用下,當系統(tǒng)可變參數(shù)(Cs,RL)在一定范圍內(nèi)發(fā)生變化時,輸出仍能保持上升時間短、穩(wěn)態(tài)誤差小、魯棒性能好的效果,從而提高系統(tǒng)的魯棒性與傳輸性能。本文利用Matlab/Simulink平臺建立了閉環(huán)系統(tǒng)的仿真模型,給出了仿真結(jié)果,并進行了分析;基于理論分析和仿真結(jié)果搭建了系統(tǒng)實驗平臺,給出實驗結(jié)果。通過仿真結(jié)果和實驗驗證了上述研究內(nèi)容的有效性。

    [1]蘇玉剛,陳龍,王智慧,等.基于最小二乘法的感應電能無線傳輸系統(tǒng)負載辨識方法[J].電工技術(shù)學報,2015,30(5):9-14.

    Su Yugang,Chen Long,Wang Zhihui,et al.A load identification method for inductive power transfer system based on least squares algorithm[J].Transactions of China Electrotechnical Society,2015,30(5):9-14.

    [2]張獻,楊慶新,崔玉龍,等.大功率無線電能傳輸系統(tǒng)能量發(fā)射線圈設計、優(yōu)化與驗證[J].電工技術(shù)學報,2013,28(10):12-18.

    Zhang Xian,Yang Qingxin,Cui Yulong,et al.Design optimization and verification on the power transmitting coilin the high-power wireless power transmission system[J].Transactions of China Electrotechnical Society,2013,28(10):12-18.

    [3]Wang Z H,Li Y P,Sun Y,et al.Load detection model of voltage-fed inductive power transfer system[J].IEEE Transactions on Power Electronics,2013,28(11):5233- 5243.

    [4]王智慧,呂瀟,孫躍,等.諧振式無線電能傳輸系統(tǒng)損耗模型[J].電工技術(shù)學報,2014,29(9):17-21.

    Wang Zhihui,Lü Xiao,Sun Yue,et al.Modeling of power loss in resonant wireless power transfer system [J].Transactions of China Electrotechnical Society,2014,29(9):17-21.

    [5]程志遠,朱春波,魏國,等.串-并補償結(jié)構(gòu)大功率感應充電系統(tǒng)諧振變換器[J].電工技術(shù)學報,2014,29(9):44-48.

    Cheng Zhiyuan,Zhu Chunbo,Wei Guo,et al.Resonant converter for high power ICPT system withseries-parallel compensation[J].Transactions of China Electrotechnical Society,2014,29(9):44-48.

    [6]楊慶新,章鵬程,祝麗花,等.無線電能傳輸技術(shù)的關(guān)鍵基礎與技術(shù)瓶頸問題[J].電工技術(shù)學報,2015,30(5):1-8.

    Yang Qingxin,Zhang Pengcheng,Zhu Lihua,et al.Key fundamental problems and technical bottlenecks of the wireless power transmission technology[J].Transactions of China Electrotechnical Society,2015,30(5):1-8.

    [7]侯佳,陳乾宏,嚴開沁,等.新型S/SP補償?shù)姆墙佑|諧振變換器分析與控制[J].中國電機工程學報,2013,33(33):1-9.

    Hou Jia,Chen Qianhong,Yan Kaiqin,et al.Analysis and control of S/SP compensation contactless resonant converters[J].Proceedings of the CESS,2013,33(33):1-9.

    [8]Liu C,Hu A P,Wang B,et al.Capacitively coupled contactless matrix charging platform with soft switched transformer control[J].Industrial Electronics,2013,60(1):249-260.

    [9]Theodoridis M P.Effective capacitive power transfer[J].IEEE Transactions on Power Electronics,2012,27(12):4906-4913.

    [10]Liu C,Hu A P,Nair N K C.Coupling study of a rotary capacitive power transfer system[C]//IEEE International Conference on Industrial Technology,Melbourne,Australia,2009:1-6.

    [11]Liu C,Hu A P.Steady state analysis of a capacitively coupled contactless power transfer system[C]//IEEE Energy Conversion Congress and Exposition,San Jose,CA,2009:3233-3238.

    [12]蘇玉剛,徐健,謝詩云,等.電場耦合型無線電能傳輸系統(tǒng)調(diào)諧技術(shù)[J].電工技術(shù)學報,2013,28(11):189-194.

    Su Yugang,Xu Jian,Xie Shiyun,et al.A tuning technology of electrical-field coupledwireless power transfer system[J].Transactions of China Electrotechnical Society,2013,28(11):189-194.

    [13]Liu C,Hu A P.Power flow control of a capacitively coupled contactless power transfer system[C]//2009 35th Annual Conference of IEEE Industrial Electronics,Porto,Portugal,2009:743-747.

    [14]Neath M,Swain A K,Madawala U K,et al.An optimal PID controller for a bidirectional inductive power transfer system using multiobjective genetic algorithm[J].IEEE Transactions on Power Electronics,2014,29(3):1523-1531.

    [15]張炯,樓佩煌,錢曉明,等.基于雙LCL補償?shù)臒o接觸供電系統(tǒng)研究[J].電工技術(shù)學報,2013,28(10):19-24.

    Zhang Jiong,Lou Peihuang,Qian Xiaoming,et al.Study of contactless power transfer system based on double LCL Compensation[J].Transactions of China Electrotechnical Society,2013,28(10):19-24.

    [16]蔡華,史黎明,李耀華,等.感應耦合電能傳輸系統(tǒng)輸出功率調(diào)節(jié)方法[J].電工技術(shù)學報,2014,29(1):215-220.

    Cai Hua,Shi Liming,Li Yaohua,et al.Output power adjustment in inductively coupled power transfer system[J].Transactions of China Electrotechnical Society,2014,29(1):215-220.

    [17]Li Y L,Sun Y,Dai X.μ-Synthesis for frequency uncertainty of the ICPT system[J].IEEE Transactions on Industrial Electronics,2013,60(1):291-300.

    [18]劉洋.ICPT系統(tǒng)中信號雙向傳輸機理研究[D].重慶:重慶大學,2013.

    [19]孫小凌.可分區(qū)供電電源負載識別與系統(tǒng)建模仿真[D].重慶:重慶大學,2009.

    [20]Gude J J,Kahoraho E.Modified ziegler-nichols method for fractional PI controllers[C]//IEEE Conference on Emerging Technologies & Factory Automation,Bilbao,2010:1-5.

    [21]Dai Xin,Zou Yang,Sun Yue.Uncertainty modeling and robust control for LCL resonant inductive power transfer system[J].Journal of Power Electronics,2013,5(13):814-828.

    Evolutionary Multi-Objective Optimization of PID Parameters for Output Voltage Regulation in ECPT System Based on NSGA-Ⅱ

    Su Yugang1Chen Lingzhi1Tang Chunsen1Ma Junhao1Hu Aiguo2

    (1.Automation College of Chongqing UniversityChongqing400030China 2.Department of Engineering the University of AucklandAuckland1010New Zealand)

    Several issues including output voltage distortion and system instability might be caused by the high-order resonant networks,nonlinearity,and parameter-sensitivity in electrical-field coupled power transfer(ECPT) systems.An optimal controller design method is proposed in order to overcome these problems.A generalized state space averaging(GSSA) model is built to approximate the ECPT system firstly.Basing on this model,an evolutionary optimization framework is utilized to search for the optimal controllers automatically with the evolutionary multi-objective algorithm NSGA-II.As the results,the difficulty of designing a proportion-integration-differentiation(PID) controller of a high-order systems is solved effectively.And the system performance,indicated by rise time,steady state error,and robustness of the system,are optimized.Meanwhile,the system overshoot is limited within an acceptable range and the stability of the system is improved.At last,the simulation and experiments verify the accuracy of the GSSA model and the validity of the utilized NSGA-II optimization method.

    Wireless power transfer,electrical-field coupled,parameter optimization,generalized state space averaging (GSSA) model,PID control

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

    2015-06-10改稿日期2015-08-06

    TM724

    蘇玉剛男,1962年生,博士,教授,研究方向為無線電能傳輸技術(shù)、電力電子技術(shù)、控制理論應用與自動化系統(tǒng)集成。

    E-mail:420888970@qq.com(通信作者)

    陳苓芷女,1991年生,碩士研究生,研究方向為電力電子技術(shù)和無線電能傳輸技術(shù)。

    E-mail:368729089@qq.com

    猜你喜歡
    控制參數(shù)耦合控制器
    非Lipschitz條件下超前帶跳倒向耦合隨機微分方程的Wong-Zakai逼近
    高超聲速飛行器滑??刂茀?shù)整定方法設計*
    飛控與探測(2022年6期)2022-03-20 02:16:14
    Birkhoff系統(tǒng)穩(wěn)定性的動力學控制1)
    力學學報(2020年4期)2020-08-11 02:32:12
    基于PI與準PR調(diào)節(jié)的并網(wǎng)逆變器控制參數(shù)設計
    黑龍江電力(2017年1期)2017-05-17 04:25:08
    基于“殼-固”耦合方法模擬焊接裝配
    大型鑄鍛件(2015年5期)2015-12-16 11:43:20
    模糊PID控制器設計及MATLAB仿真
    MOXA RTU控制器ioPAC 5542系列
    自動化博覽(2014年9期)2014-02-28 22:33:17
    倍福 CX8091嵌入式控制器
    自動化博覽(2014年4期)2014-02-28 22:31:15
    求解奇異攝動Volterra積分微分方程的LDG-CFEM耦合方法
    非線性耦合KdV方程組的精確解
    国产一区二区激情短视频 | 亚洲精品美女久久久久99蜜臀| 国产精品国产三级国产专区5o| 人人妻人人添人人爽欧美一区卜| 后天国语完整版免费观看| 老司机影院毛片| 久久久久视频综合| 人妻 亚洲 视频| videos熟女内射| 丝袜在线中文字幕| 搡老熟女国产l中国老女人| 亚洲精品国产一区二区精华液| 免费观看a级毛片全部| 国产免费福利视频在线观看| 亚洲国产精品999| 国产精品久久久久久精品古装| av不卡在线播放| 久久青草综合色| 亚洲精品乱久久久久久| 亚洲国产精品一区二区三区在线| 国产欧美亚洲国产| 久久天堂一区二区三区四区| 日本五十路高清| 久久久精品区二区三区| 在线十欧美十亚洲十日本专区| 国产人伦9x9x在线观看| 在线av久久热| 伊人亚洲综合成人网| 成人亚洲精品一区在线观看| 十八禁高潮呻吟视频| 国产一级毛片在线| 老熟妇乱子伦视频在线观看 | 亚洲国产精品成人久久小说| 国产1区2区3区精品| 欧美大码av| 夜夜夜夜夜久久久久| 亚洲色图 男人天堂 中文字幕| 久久人妻福利社区极品人妻图片| 久久久水蜜桃国产精品网| 蜜桃在线观看..| 精品视频人人做人人爽| 不卡一级毛片| 国产精品久久久久久精品电影小说| 9色porny在线观看| 日本猛色少妇xxxxx猛交久久| 欧美黄色淫秽网站| 国产色视频综合| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲国产精品一区三区| 久久久国产成人免费| 中亚洲国语对白在线视频| 香蕉丝袜av| 成年人免费黄色播放视频| 日韩精品免费视频一区二区三区| 搡老岳熟女国产| 精品亚洲成国产av| 日韩欧美国产一区二区入口| 国产亚洲av高清不卡| 女性被躁到高潮视频| 亚洲精品国产av成人精品| 亚洲精品久久午夜乱码| 国产精品一区二区免费欧美 | 在线亚洲精品国产二区图片欧美| 麻豆乱淫一区二区| 一二三四在线观看免费中文在| 免费日韩欧美在线观看| 在线观看舔阴道视频| 搡老熟女国产l中国老女人| 久久久久久久精品精品| 我的亚洲天堂| 9色porny在线观看| 人人澡人人妻人| 91老司机精品| 国产成人av激情在线播放| 国产精品一区二区免费欧美 | 久久久久久免费高清国产稀缺| 老司机在亚洲福利影院| 午夜两性在线视频| 欧美日韩精品网址| 在线观看舔阴道视频| 看免费av毛片| xxxhd国产人妻xxx| 男女高潮啪啪啪动态图| 精品久久蜜臀av无| 久久久国产成人免费| 久久性视频一级片| 超色免费av| 十八禁人妻一区二区| 人成视频在线观看免费观看| 香蕉丝袜av| 亚洲久久久国产精品| 亚洲精品美女久久久久99蜜臀| 久久亚洲国产成人精品v| 久久精品aⅴ一区二区三区四区| 波多野结衣av一区二区av| 国产真人三级小视频在线观看| 午夜影院在线不卡| 女人高潮潮喷娇喘18禁视频| 日本五十路高清| 国产真人三级小视频在线观看| avwww免费| 午夜福利一区二区在线看| 久久中文看片网| 啦啦啦 在线观看视频| 亚洲精品粉嫩美女一区| 久久精品人人爽人人爽视色| 美女高潮喷水抽搐中文字幕| 99精品欧美一区二区三区四区| 精品少妇黑人巨大在线播放| 国产精品九九99| 欧美日韩中文字幕国产精品一区二区三区 | 国产在线观看jvid| videos熟女内射| 国产深夜福利视频在线观看| 精品人妻一区二区三区麻豆| 自线自在国产av| 国产日韩欧美亚洲二区| 男女午夜视频在线观看| 在线观看人妻少妇| 丝瓜视频免费看黄片| 少妇精品久久久久久久| 久热这里只有精品99| 在线精品无人区一区二区三| 中文字幕人妻丝袜一区二区| 性高湖久久久久久久久免费观看| 岛国毛片在线播放| 肉色欧美久久久久久久蜜桃| 久久中文字幕一级| 日韩视频一区二区在线观看| 日韩视频一区二区在线观看| 亚洲欧美一区二区三区黑人| 日本a在线网址| 国产亚洲精品久久久久5区| 99国产综合亚洲精品| 肉色欧美久久久久久久蜜桃| 亚洲av日韩在线播放| 大型av网站在线播放| 国产有黄有色有爽视频| 母亲3免费完整高清在线观看| 亚洲精品中文字幕一二三四区 | 午夜福利,免费看| 欧美日本中文国产一区发布| av有码第一页| 在线看a的网站| 久久精品成人免费网站| 后天国语完整版免费观看| 制服诱惑二区| 成人手机av| 国产黄色免费在线视频| 这个男人来自地球电影免费观看| 欧美日韩成人在线一区二区| 最近最新免费中文字幕在线| 波多野结衣一区麻豆| 美女主播在线视频| 18禁黄网站禁片午夜丰满| 日韩制服骚丝袜av| 日韩制服骚丝袜av| 国产激情久久老熟女| 91老司机精品| 啦啦啦视频在线资源免费观看| 国产成人免费无遮挡视频| 天天躁日日躁夜夜躁夜夜| 久久久久网色| 国产日韩欧美视频二区| 久久久精品区二区三区| 女人被躁到高潮嗷嗷叫费观| 99久久人妻综合| 热99久久久久精品小说推荐| www.999成人在线观看| 国内毛片毛片毛片毛片毛片| 最近最新中文字幕大全免费视频| 电影成人av| 免费观看a级毛片全部| 各种免费的搞黄视频| 久久影院123| 99国产精品一区二区三区| 亚洲精品第二区| 亚洲精品国产av蜜桃| 国产精品香港三级国产av潘金莲| 91精品三级在线观看| 欧美激情高清一区二区三区| 国产精品成人在线| 老司机福利观看| 淫妇啪啪啪对白视频 | 激情视频va一区二区三区| 人人妻,人人澡人人爽秒播| avwww免费| 日本精品一区二区三区蜜桃| 久久精品人人爽人人爽视色| 黄色视频,在线免费观看| 国产黄色免费在线视频| 成人黄色视频免费在线看| 久久性视频一级片| 国产黄色免费在线视频| 99热网站在线观看| 黑人猛操日本美女一级片| 99久久综合免费| 久热爱精品视频在线9| 欧美日韩亚洲高清精品| av在线老鸭窝| tube8黄色片| 亚洲精品一二三| 岛国毛片在线播放| 在线看a的网站| 欧美日韩福利视频一区二区| 99热网站在线观看| 精品第一国产精品| 免费一级毛片在线播放高清视频 | 欧美黑人精品巨大| 国产精品香港三级国产av潘金莲| 美女高潮喷水抽搐中文字幕| 国产成人精品在线电影| 日韩人妻精品一区2区三区| 最黄视频免费看| 美女扒开内裤让男人捅视频| 啦啦啦中文免费视频观看日本| 老司机亚洲免费影院| 亚洲欧美日韩另类电影网站| 亚洲精品国产区一区二| 天天躁狠狠躁夜夜躁狠狠躁| 1024香蕉在线观看| 亚洲性夜色夜夜综合| 国产色视频综合| 亚洲av男天堂| 制服人妻中文乱码| 又黄又粗又硬又大视频| 久久久久久久大尺度免费视频| 国产成人系列免费观看| 精品人妻熟女毛片av久久网站| 久久99一区二区三区| 精品国产乱码久久久久久小说| 精品国产一区二区久久| 欧美人与性动交α欧美精品济南到| 69精品国产乱码久久久| 十八禁人妻一区二区| 国产精品亚洲av一区麻豆| 50天的宝宝边吃奶边哭怎么回事| 免费日韩欧美在线观看| 婷婷色av中文字幕| 久久亚洲精品不卡| 国产精品一区二区在线观看99| 狠狠狠狠99中文字幕| 亚洲人成电影免费在线| 日本撒尿小便嘘嘘汇集6| 国产高清国产精品国产三级| 日韩免费高清中文字幕av| 午夜精品国产一区二区电影| 久久亚洲精品不卡| 日韩视频一区二区在线观看| 国产成人影院久久av| 久久久久国产一级毛片高清牌| 色视频在线一区二区三区| 亚洲自偷自拍图片 自拍| 每晚都被弄得嗷嗷叫到高潮| 久久狼人影院| 午夜福利免费观看在线| a在线观看视频网站| 高清av免费在线| 国产免费一区二区三区四区乱码| 每晚都被弄得嗷嗷叫到高潮| 欧美黑人欧美精品刺激| 嫁个100分男人电影在线观看| 国产高清国产精品国产三级| av不卡在线播放| 午夜福利视频在线观看免费| 国产精品一区二区在线观看99| av福利片在线| 亚洲av电影在线观看一区二区三区| 国产精品99久久99久久久不卡| 丝袜喷水一区| 久久国产精品大桥未久av| 成年人黄色毛片网站| 国产高清videossex| 丝袜在线中文字幕| 免费在线观看完整版高清| 美女大奶头黄色视频| 国产区一区二久久| 夜夜骑夜夜射夜夜干| 女人高潮潮喷娇喘18禁视频| 亚洲国产日韩一区二区| 免费观看av网站的网址| 成人国产av品久久久| 日韩有码中文字幕| 久久国产亚洲av麻豆专区| 亚洲欧美色中文字幕在线| 一区二区三区精品91| 精品卡一卡二卡四卡免费| 成年动漫av网址| 制服诱惑二区| 美女福利国产在线| 纯流量卡能插随身wifi吗| 亚洲精品久久成人aⅴ小说| 老司机影院毛片| 最黄视频免费看| 熟女少妇亚洲综合色aaa.| 搡老熟女国产l中国老女人| 涩涩av久久男人的天堂| 国产男女超爽视频在线观看| 国产欧美日韩一区二区精品| 啦啦啦中文免费视频观看日本| 女人高潮潮喷娇喘18禁视频| 国产精品国产av在线观看| 一级黄色大片毛片| 午夜精品久久久久久毛片777| 国产日韩欧美在线精品| 亚洲国产av新网站| 精品久久久久久电影网| 一本—道久久a久久精品蜜桃钙片| 啦啦啦视频在线资源免费观看| 国产无遮挡羞羞视频在线观看| 一本一本久久a久久精品综合妖精| 亚洲激情五月婷婷啪啪| 亚洲精品国产色婷婷电影| 少妇 在线观看| 宅男免费午夜| 欧美精品人与动牲交sv欧美| 男女高潮啪啪啪动态图| 亚洲五月色婷婷综合| 美女高潮喷水抽搐中文字幕| 麻豆国产av国片精品| 天堂8中文在线网| 女人被躁到高潮嗷嗷叫费观| 午夜福利免费观看在线| 国产成人av教育| 精品少妇内射三级| 日韩欧美国产一区二区入口| 国产精品自产拍在线观看55亚洲 | 国产精品一区二区精品视频观看| 精品一品国产午夜福利视频| 久久精品国产亚洲av高清一级| 国产成人啪精品午夜网站| 黄色视频在线播放观看不卡| 精品视频人人做人人爽| cao死你这个sao货| 精品一区在线观看国产| 国产精品国产三级国产专区5o| 99精品欧美一区二区三区四区| 国产精品偷伦视频观看了| 黑人欧美特级aaaaaa片| 三上悠亚av全集在线观看| 国产精品一区二区在线观看99| 日韩有码中文字幕| 无限看片的www在线观看| 日韩欧美国产一区二区入口| www.999成人在线观看| 麻豆国产av国片精品| 看免费av毛片| 老司机靠b影院| 免费高清在线观看日韩| 成人三级做爰电影| 韩国高清视频一区二区三区| 日韩,欧美,国产一区二区三区| 操出白浆在线播放| 亚洲欧美一区二区三区久久| 国产精品麻豆人妻色哟哟久久| 大型av网站在线播放| 91大片在线观看| 亚洲精品中文字幕在线视频| 日韩欧美一区二区三区在线观看 | 亚洲欧美一区二区三区黑人| 免费av中文字幕在线| 亚洲国产欧美日韩在线播放| 老司机影院毛片| 欧美黄色淫秽网站| 久热爱精品视频在线9| 国产一区二区三区综合在线观看| 三级毛片av免费| 国产老妇伦熟女老妇高清| 久久毛片免费看一区二区三区| 成人黄色视频免费在线看| 黑人操中国人逼视频| 久久国产精品男人的天堂亚洲| av超薄肉色丝袜交足视频| 国产精品一区二区免费欧美 | 啪啪无遮挡十八禁网站| 黄片播放在线免费| 成人国产一区最新在线观看| 国产福利在线免费观看视频| 国产成人啪精品午夜网站| 天堂俺去俺来也www色官网| 女性生殖器流出的白浆| 99热网站在线观看| 狂野欧美激情性bbbbbb| 亚洲第一青青草原| 日韩制服骚丝袜av| 亚洲国产精品成人久久小说| 黄片小视频在线播放| 中文字幕人妻熟女乱码| 国产成人精品久久二区二区免费| 欧美在线黄色| 亚洲美女黄色视频免费看| 日日爽夜夜爽网站| 涩涩av久久男人的天堂| 国产精品久久久久久人妻精品电影 | 国产亚洲欧美在线一区二区| 色婷婷av一区二区三区视频| 在线 av 中文字幕| 蜜桃在线观看..| 三级毛片av免费| 黄色视频在线播放观看不卡| 日本黄色日本黄色录像| 老司机午夜十八禁免费视频| 97人妻天天添夜夜摸| 久久久国产精品麻豆| 2018国产大陆天天弄谢| 飞空精品影院首页| 最近中文字幕2019免费版| 精品欧美一区二区三区在线| 久热爱精品视频在线9| 天天添夜夜摸| 国产人伦9x9x在线观看| 亚洲专区中文字幕在线| 免费高清在线观看视频在线观看| 亚洲精品国产色婷婷电影| 国产精品国产三级国产专区5o| 久久精品人人爽人人爽视色| 一级毛片电影观看| 精品少妇一区二区三区视频日本电影| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲第一青青草原| 菩萨蛮人人尽说江南好唐韦庄| 黄色视频在线播放观看不卡| 久久久国产精品麻豆| 国产精品免费大片| 午夜视频精品福利| 91av网站免费观看| 一级片'在线观看视频| 18禁裸乳无遮挡动漫免费视频| 秋霞在线观看毛片| 日韩欧美免费精品| 淫妇啪啪啪对白视频 | 女性生殖器流出的白浆| 亚洲综合色网址| 男人舔女人的私密视频| 大型av网站在线播放| 亚洲 欧美一区二区三区| 国产日韩欧美在线精品| 亚洲熟女精品中文字幕| 色婷婷久久久亚洲欧美| 国产精品二区激情视频| 人妻久久中文字幕网| 91字幕亚洲| 国产黄色免费在线视频| 国产免费一区二区三区四区乱码| 精品人妻在线不人妻| 国产极品粉嫩免费观看在线| 高清欧美精品videossex| 亚洲成av片中文字幕在线观看| 婷婷丁香在线五月| 好男人电影高清在线观看| 中亚洲国语对白在线视频| 黑人操中国人逼视频| 欧美黄色片欧美黄色片| 老司机福利观看| 狠狠精品人妻久久久久久综合| 国产精品av久久久久免费| 中文字幕最新亚洲高清| 丝瓜视频免费看黄片| 久久九九热精品免费| 美国免费a级毛片| 黄片小视频在线播放| 午夜久久久在线观看| av线在线观看网站| 欧美成狂野欧美在线观看| 欧美精品啪啪一区二区三区 | 一区在线观看完整版| 成年人午夜在线观看视频| 国产一区二区三区综合在线观看| 亚洲国产av影院在线观看| 中文字幕精品免费在线观看视频| 99精品欧美一区二区三区四区| e午夜精品久久久久久久| 日韩 亚洲 欧美在线| 岛国毛片在线播放| 可以免费在线观看a视频的电影网站| 欧美精品亚洲一区二区| 99国产精品免费福利视频| 欧美日韩亚洲国产一区二区在线观看 | 欧美亚洲 丝袜 人妻 在线| 一本久久精品| 99精品欧美一区二区三区四区| 人妻久久中文字幕网| 在线亚洲精品国产二区图片欧美| 亚洲av欧美aⅴ国产| 日日摸夜夜添夜夜添小说| 天天影视国产精品| 国产黄频视频在线观看| 色婷婷久久久亚洲欧美| 欧美激情极品国产一区二区三区| 国产高清videossex| 精品人妻一区二区三区麻豆| 青草久久国产| 男人爽女人下面视频在线观看| 午夜福利,免费看| 日韩中文字幕视频在线看片| 老熟妇仑乱视频hdxx| 国产成人精品无人区| 99国产精品99久久久久| 亚洲,欧美精品.| 在线av久久热| 每晚都被弄得嗷嗷叫到高潮| 久久免费观看电影| 亚洲精品国产av成人精品| 精品人妻熟女毛片av久久网站| 老司机午夜福利在线观看视频 | 成年人黄色毛片网站| 啦啦啦视频在线资源免费观看| 久久中文看片网| 桃花免费在线播放| 啦啦啦啦在线视频资源| 天堂中文最新版在线下载| 亚洲精品av麻豆狂野| 午夜福利乱码中文字幕| 国产精品秋霞免费鲁丝片| 婷婷色av中文字幕| 母亲3免费完整高清在线观看| 久久久久久久久免费视频了| 欧美在线黄色| 老司机午夜十八禁免费视频| 美女视频免费永久观看网站| 麻豆av在线久日| 夜夜骑夜夜射夜夜干| 午夜福利免费观看在线| 国产一区二区 视频在线| 国产xxxxx性猛交| 黄色片一级片一级黄色片| 97精品久久久久久久久久精品| 久久人妻熟女aⅴ| 真人做人爱边吃奶动态| 亚洲精品一二三| 国产精品免费大片| 97在线人人人人妻| 久久久国产成人免费| 麻豆乱淫一区二区| 香蕉丝袜av| 日日爽夜夜爽网站| 国产精品国产av在线观看| 夜夜夜夜夜久久久久| 亚洲精品国产av蜜桃| 亚洲精品中文字幕一二三四区 | 亚洲伊人久久精品综合| 欧美 日韩 精品 国产| 亚洲五月婷婷丁香| 亚洲欧美日韩另类电影网站| 伊人亚洲综合成人网| 一区二区三区乱码不卡18| 一二三四社区在线视频社区8| 免费av中文字幕在线| 美女中出高潮动态图| av又黄又爽大尺度在线免费看| 99国产极品粉嫩在线观看| 99热国产这里只有精品6| 久久久久久免费高清国产稀缺| 男女边摸边吃奶| 国产精品久久久久成人av| 国产精品九九99| 亚洲av电影在线观看一区二区三区| 欧美97在线视频| 日韩制服骚丝袜av| 亚洲国产成人一精品久久久| 两性夫妻黄色片| 国产日韩欧美视频二区| 久久久久久久久久久久大奶| 免费在线观看视频国产中文字幕亚洲 | 亚洲一区二区三区欧美精品| 一二三四在线观看免费中文在| a级毛片在线看网站| 一二三四在线观看免费中文在| 亚洲一区二区三区欧美精品| 秋霞在线观看毛片| 成年人黄色毛片网站| 在线亚洲精品国产二区图片欧美| 三级毛片av免费| 亚洲专区国产一区二区| 狠狠婷婷综合久久久久久88av| 欧美乱码精品一区二区三区| 91av网站免费观看| 五月天丁香电影| 亚洲精品自拍成人| 在线观看www视频免费| av福利片在线| 精品人妻一区二区三区麻豆| 国产精品亚洲av一区麻豆| 精品国产一区二区三区四区第35| 国产男女内射视频| 久久亚洲国产成人精品v| 色视频在线一区二区三区| 岛国毛片在线播放| 王馨瑶露胸无遮挡在线观看| 国产有黄有色有爽视频| 亚洲va日本ⅴa欧美va伊人久久 | av网站在线播放免费| 亚洲精品久久成人aⅴ小说| 欧美大码av| 久久久久精品人妻al黑| 国产成人a∨麻豆精品| 免费观看人在逋| 下体分泌物呈黄色| 久久久久国内视频| 国产一级毛片在线| 亚洲国产av影院在线观看| 少妇精品久久久久久久| 91麻豆av在线| 国产精品国产av在线观看| 色视频在线一区二区三区| 国产三级黄色录像| 色播在线永久视频| 肉色欧美久久久久久久蜜桃| 亚洲午夜精品一区,二区,三区| 一级,二级,三级黄色视频| 久久久精品区二区三区| 中文字幕人妻熟女乱码| videos熟女内射| 欧美另类亚洲清纯唯美| 国产一卡二卡三卡精品|