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

    一種改進的內(nèi)嵌式永磁同步電機虛擬信號注入MTPA控制

    2022-10-12 04:16:44嚴晶鋮曾成碧朱奇昊李佩霜邱仲睿
    關(guān)鍵詞:矢量轉(zhuǎn)矩幅值

    嚴晶鋮,曾成碧,苗 虹,朱奇昊,李佩霜,邱仲睿

    (四川大學 電氣工程學院, 成都 610065)

    0 引言

    內(nèi)嵌式永磁同步機(IPMSM)因結(jié)構(gòu)簡單、功率密度高、調(diào)速范圍寬等優(yōu)點被廣泛應(yīng)用于航空航天等領(lǐng)域。最大轉(zhuǎn)矩電流比(maximum torque per ampere,MTPA)控制是在一定的輸出轉(zhuǎn)矩條件下,通過dq軸電流的合理分配降低定子電流及銅耗,提升電機轉(zhuǎn)矩輸出能力和系統(tǒng)效率,在IPMSM的控制中得到廣泛關(guān)注。

    精準的MTPA控制對參數(shù)依賴性極大,然而在不同工況下,受溫度、磁飽和等因素的影響,IPMSM交直軸電感、磁鏈、電阻等均會發(fā)生非線性變化。如何降低MTPA對參數(shù)的依賴,從而減少因電機參數(shù)不確定而產(chǎn)生的的控制誤差,是目前MTPA控制策略的一個研究熱點。

    為了克服電磁參數(shù)非線性的問題,目前的研究大體可以分為離線和在線的方案[1]。通過離線策略最常用的查表法,可得到電機參數(shù)或者最優(yōu)電流參考值[2],但不具有普適性且需占據(jù)大量的存儲空間。在線策略主要包括在線搜索法、代價函數(shù)法和信號注入法。搜索法通過向電流矢量角施加擾動在線搜索最優(yōu)工作點,但其動態(tài)響應(yīng)性能不佳[3]。代價函數(shù)法通過構(gòu)建預(yù)測控制的代價函數(shù),實現(xiàn) MTPA 控制與模型預(yù)測控制的融合,但準確性無法得到保障[4]。信號注入法是向電機的電流或者電壓注入高頻的正弦、方波信號[5-7],對含高頻擾動信號的響應(yīng)信息進行調(diào)制,提取到正比于的信息,通過控制電流矢量角度,實時完成MTPA工作點的跟蹤。但高頻信號的注入會產(chǎn)生高頻的轉(zhuǎn)矩擾動并引入額外損耗,還可能會造成電磁干擾的影響。

    考慮到上述的問題,文獻[8]提出虛擬信號注入(virtual signals injection,VSI)的策略,不同于真實信號注入法,虛擬信號注入法對的提取是在虛擬轉(zhuǎn)矩或虛擬功率模型上,由虛擬信號響應(yīng)的計算得到。虛擬信號注入法具有對電機參數(shù)變化不敏感的特點,并克服了真實信號注入所產(chǎn)生的轉(zhuǎn)矩脈動及損耗問題。為進一步提高虛擬信號注入法的控制精度,降低由于虛擬轉(zhuǎn)矩或虛擬功率模型誤差帶來的控制偏差,文獻[9]提出虛擬多信號注入的控制方法,引入了虛擬功率的概念,避開了濾波器的延時效應(yīng)以提高控制系統(tǒng)的動態(tài)性能,但是虛擬功率算子對電磁參數(shù)的依賴性將降低MTPA控制的準確性。文獻[10]提出向電流矢量角注入虛擬方波信號,依托轉(zhuǎn)矩在當前矢量角度下的泰勒展開式建立了不依賴電磁參數(shù)的虛擬轉(zhuǎn)矩模型,然而由于轉(zhuǎn)矩是當前矢量角度的三角函數(shù),對高階偏導數(shù)的忽略會在一定程度上影響提取信息的準確性。文獻[11]在2種虛擬信號注入法基礎(chǔ)上,提出基于直流虛擬信號注入法的控制策略,但是該方法并沒有考慮到虛擬電流注入對虛擬功率模型中電流幅值的改變所產(chǎn)生的誤差。同時,針對系統(tǒng)在極值處波動的問題,文獻[12]提出自適應(yīng)的可調(diào) 振幅虛擬方波信號,從而得到波動較小的最優(yōu)電流角度,但并提出相應(yīng)閾值的理論依據(jù)。

    本文中提出一種更加精準的虛擬直流信號注入法的MTPA控制,針對一般虛擬直流信號中虛擬功率模型產(chǎn)生的誤差,定義了新的直流信號的注入準則,提高了MTPA判據(jù)的算法效率和精度。為驗證虛擬直流信號的可行性,分析了系統(tǒng)在極值點附近的轉(zhuǎn)矩振蕩和局部穩(wěn)定性問題,探討了控制參數(shù)對MTPA控制性能的影響。通過理論分析和仿真實驗驗證本文提出的控制策略的有效性。

    1 基于虛擬信號注入的MTPA控制

    在正弦定子電流激勵下,IPMSM在dq參考坐標系的數(shù)學模型中電壓方程為

    (1)

    電機電磁轉(zhuǎn)矩方程為

    (2)

    (3)

    式中:vq和vd分別表示交直軸電壓,電機的交直軸電感分別由Lq和Ld表示,永磁磁鏈為ψm。Lq,Ld和ψm隨電機運行的工況變化而改變,可以將其視為id和iq函數(shù)。p為電機的極對數(shù),ωm為電機轉(zhuǎn)子轉(zhuǎn)速。在dq參考軸下的電流幅值為Ia,電流矢量角為β。為滿足電機在穩(wěn)態(tài)運行下MTPA追蹤,通過改變id和iq以達到最大電磁轉(zhuǎn)矩的控制。

    傳統(tǒng)的虛擬直流信號注入法分別向交直軸注入直流信號,如圖1所示。該方法忽略了電流幅值變化產(chǎn)生的誤差,不同直流偏置的注入方法將產(chǎn)生不同的誤差,這將導致?Te/?β的誤差隨電流矢量位置實時變化??紤]電流變化因素,?Te/?β表達式為:

    (4)

    圖1 傳統(tǒng)虛擬直流信號注入法示意圖

    根據(jù)上述分析和電流矢量的關(guān)系,考慮電流矢量變化的影響,對任意的虛擬直流信號注入所產(chǎn)生的誤差可以表示為:

    (5)

    傳統(tǒng)的虛擬直流信號注入法中,?Te/?β是通過分別注入相同大小直流信號后合并計算的結(jié)果,因此無法保證在交直軸電流注入同樣的電流偏置后?id/?Ia,?iq/?Ia偏導項與依據(jù)式(3)得到的偏導相互對應(yīng)。在合理的直流電流A注入的情況下,傳統(tǒng)方法MTPA的誤差可以表示為

    (6)

    其中

    (7)

    式中:A表示注入的直流信號電流;ΔIa1表示單獨向d軸電流注入直流偏置引起電流幅值的改變;ΔIa2表示單獨向q軸電流注入直流偏置引起電流幅值的改變。通過上述分析,分別在交直軸注入直流信號方法將致使ΔIa1,ΔIa2隨不同的id,iq變化而變化,ΔIa1≠ΔIa2導致該方案只能近似估計MTPA控制的誤差。

    2 自適應(yīng)的虛擬直流信號注入法的MTPA控制

    考慮到一般虛擬直流信號注入法產(chǎn)生MTPA控制的穩(wěn)態(tài)誤差,提出一種自適應(yīng)的VSI控制策略,通過注入一個幅值恒定、角度動態(tài)變化的直流信號,求取不同工況下的MTPA下的?Te/?β,利用該導數(shù)項對id,iq進行控制,進而達到MTPA的控制目的。注入電流信號角度的最優(yōu)選擇實現(xiàn)自適應(yīng)控制過程,優(yōu)化?Te/?β計算效率,并保證MTPA的控制精度。

    在利用式(4)求取偏導時,?id/?Ia, ?iq/?Ia對偏導的影響不大。而在實際運行過程中,電流幅值變化更多是受到采樣信號的波動。由于直流信號所引起的電流幅值變化很小,導致對?id/?Ia, ?iq/?Ia進行估測增加了控制算法的復(fù)雜度,同時降低了準確性。因此提出最優(yōu)的注入電流角度自適應(yīng)律,角度應(yīng)滿足如下條件:

    (8)

    動態(tài)的注入信號可得

    (9)

    其中

    (10)

    式中,σ表示注入信號的角度。

    圖2為改進的虛擬直流信號注入示意圖。當負載轉(zhuǎn)矩保持不變,對于給定直軸電流參考值idref,交軸電流參考值iqref由此可以確定,因此對交直流電流的控制可以等效于改變電流矢量角度β。

    圖2 改進的虛擬直流信號注入示意圖

    (11)

    圖3 虛擬直流信號注入所對應(yīng)轉(zhuǎn)矩

    當注入的電流幅值足夠小時,轉(zhuǎn)矩產(chǎn)生的變化分量中一階項起到主導作用,因此在忽略高次項影響的情況下,得到

    (12)

    IPMSM的電磁轉(zhuǎn)矩輸出功率可以表示為

    (13)

    式中,ωe為電機的電角速度。

    在虛擬直流信號注入法的基礎(chǔ)上,根據(jù)式(12),交直軸電流的改變對應(yīng)的虛擬的功率可以表示為

    (14)

    IPMSM中的永磁體磁鏈ψm、Ld、Lq都因為磁飽和效應(yīng)而受到交直軸電流的影響,依據(jù)式(1),由電壓穩(wěn)態(tài)方程可推導出:

    (15)

    交直軸注入電流以后,產(chǎn)生的虛擬功率變化可以表示為

    (16)

    (17)

    將式(16)和式(17)代入式(12)得到

    (18)

    虛擬信號注入法求取的偏導值通過積分環(huán)節(jié)得到idref,如圖4所示。根據(jù)式(19),利用轉(zhuǎn)速環(huán)得到的轉(zhuǎn)矩參考值Tref得到iqref,再通過PI電流控制便可實現(xiàn)虛擬信號注入法的控制。

    (19)

    圖4 虛擬直流信號注入MTPA控制框圖

    3 虛擬直流信號注入法的可行性分析

    虛擬信號注入的MTPA控制自提出以后,研究主要集中于該方法的精度提升和適用范圍,針對調(diào)制電壓、控制方法存在的模型誤差如何補償[9-11,14-17]。目前的研究缺少對虛擬信號穩(wěn)定性的分析,本節(jié)旨在從理論上論證虛擬直流信號注入法的可行性并分析不同設(shè)計參數(shù)選擇對控制的影響,研究轉(zhuǎn)矩在穩(wěn)態(tài)下的振蕩問題,分析系統(tǒng)參數(shù)在極值點附近的局部穩(wěn)定性,并非提供一個全局性的穩(wěn)定性證明。

    在一個穩(wěn)態(tài)時刻,可以將非線性速度控制環(huán)線性化,將電磁參數(shù)和機械參數(shù)視為常值。信號注入法的核心在于提取?Te/?β,并用于獲取電流矢量角度或直軸電流參考值,2種獲取方式因控制回路將造成電流幅值大小和交軸電流的滯后,但分析的方法是一致的。本文研究以獲取直軸電流參考值為例,對于每個控制周期而言,將直軸電流參考值的獲取過程視作為擾動信號Δid,通過系統(tǒng)控制回路改變Δiq,兩者共同決定Δβ,根據(jù)本文的MTPA策略得到?Te/?β,并影響下一個控制周期的Δid,它們之間的關(guān)系如圖5所示。

    圖5 系統(tǒng)狀態(tài)變量與控制變量關(guān)系示意圖

    為簡化虛擬信號注入產(chǎn)生的信號擾動的分析,假設(shè):

    1) 鑒于電流環(huán)回路的響應(yīng)是遠快于轉(zhuǎn)速環(huán),忽略電流控制和電流,電壓采樣的暫態(tài)影響;

    2) 系統(tǒng)運行于穩(wěn)態(tài)過程,此時合成電流矢量位于旋轉(zhuǎn)坐標軸第二象限;

    3) 利用穩(wěn)態(tài)方程求解出的虛擬功率差與轉(zhuǎn)矩關(guān)于交直軸電流的偏導項可以忽略;

    4) 將擾動信號Δid看作以大小任意的控制周期變量,令Δid=Rsin(ωrt+φr)。

    設(shè)在轉(zhuǎn)矩最大的極值點附近存在工作點S0,該工作點的參數(shù)和狀態(tài)集合見式(20)—(22)。同時建立控制回路的小擾動方程。

    S0={id=id0,iq=iq0,Ia=Ia0,

    β=β0,ωm=ωm0,TL=TL0,

    Te0=Te(id0,iq0,ψm0,Ld0,Lq0)}

    (20)

    (21)

    (22)

    其中

    (23)

    式中:ΔTe、Δido、Δiqo、Δeω、Δx表示在系統(tǒng)平衡值的變化量;KPω和KIω表示速度PI控制器的比例和積分增益。由此建立小信號模型,如圖6所示。

    圖6 控制回路的小信號模型結(jié)構(gòu)示意圖

    可以推出關(guān)于擾動的傳遞函數(shù)為

    (24)

    其中

    (25)

    為簡化分析,將系統(tǒng)的id0可以看作控制頻率wr為基礎(chǔ)頻率的周期擾動信號。忽略控制周期對系統(tǒng)的相位的滯后的影響,令

    (26)

    結(jié)合式(24)—(26),得到

    Δiq=kktkGΔid

    (27)

    iq與id的擾動可能會造成系統(tǒng)在穩(wěn)定狀態(tài)下的轉(zhuǎn)矩振蕩。為克服類似問題,控制電流矢量在iq與id的擾動下仍處于恒轉(zhuǎn)矩電流工作區(qū),由此可推得所期望的Δiqe為:

    (28)

    從式(27)可知,k和kt是由電機所處的運行工況決定的,而kG是由系統(tǒng)本身及其控制參數(shù)決定的。因此,通過對系統(tǒng)控制參數(shù)的修正,從控制理論上可避免系統(tǒng)穩(wěn)態(tài)處轉(zhuǎn)矩振蕩問題。聯(lián)立式(27)和式(28),可求得

    (29)

    由式(29)可知,不同轉(zhuǎn)矩和轉(zhuǎn)速所引起的電磁參數(shù)非線性變化,導致最優(yōu)的系統(tǒng)增益kG也隨之改變。因此在控制系統(tǒng)的設(shè)計過程中,可以通過自適應(yīng)調(diào)節(jié)不斷改變系統(tǒng)轉(zhuǎn)速環(huán)的設(shè)定將穩(wěn)態(tài)處的轉(zhuǎn)矩變化控制在較小范圍。 值得一提的是,真實信號注入是在電流上施加實際擾動,而本文分析的電流擾動是建立在控制環(huán)節(jié)本身,系統(tǒng)在穩(wěn)態(tài)下的擾動程度取決于直軸電流參考的獲取方式,因此虛擬信號注入法可以從根本解決在系統(tǒng)在穩(wěn)定狀態(tài)下轉(zhuǎn)矩振蕩問題。

    為簡化系統(tǒng)在極值處的局部穩(wěn)定問題的分析,考慮到電機的轉(zhuǎn)速與定子電流變化所產(chǎn)生交直軸電壓變化是不顯著的,可忽略其對?Te/?β提取的影響,由虛擬功率模型得到的?Te/?id和?Te/?iq是準確的,此時id的控制方法可以表示為:

    (30)

    同時,存在工作點S1滿足

    (31)

    工作點S0與工作點S1之間的系統(tǒng)狀態(tài)為本文局部穩(wěn)定性分析的范圍。單個控制周期內(nèi)的Δid可表示為Δid=-kiλ,為使得Δid收斂于0,因確保對Δid的控制和Δiq的系統(tǒng)響應(yīng)所引起的電流矢量角度的改變小于系統(tǒng)由工作點S0過渡到工作點S1的電流矢量角度變化,故Δid收斂于0的臨界條件可表示為

    Δid=Is0sinβs0-Is1sinβs1

    (32)

    聯(lián)立式(31)和(32),建立電流矢量角的方程為

    (33)

    忽略極值點附近的電流變化情況,得到:

    (34)

    積分系數(shù)的取值是本文控制系統(tǒng)中Δid=0的關(guān)鍵。積分系數(shù)的選取應(yīng)小于式(34)所提出的臨界條件,因此虛擬信號幅值大小只決定了MTPA追蹤精度,并不會影響系統(tǒng)的穩(wěn)定性。同時,系統(tǒng)的負載變化也將改變積分系數(shù)的選取范圍。文獻[12]在積分因子乘上了虛擬信號幅值,因此注入信號的振幅會對電流矢量角度的穩(wěn)態(tài)波動產(chǎn)生影響,但其方法本質(zhì)也是通過調(diào)節(jié)積分系數(shù)去改善載荷變化下的穩(wěn)定性能。

    在滿足上述條件時,系統(tǒng)存在平衡狀態(tài)。當?Te/?β=0時,Δid=0。為進一步分析系統(tǒng)在極值附近的穩(wěn)定性,建立李雅普諾夫方程:

    (35)

    根據(jù)式(3)(30)和(35),李氏函數(shù)關(guān)于時間的導數(shù)可以寫為

    (36)

    dV/dδ與did/dt異號,由式(3)算得dβ/did恒小于0,同時轉(zhuǎn)矩關(guān)于電流矢量的函數(shù)在工作區(qū)間時為凸函數(shù),故dδ/dβ恒小于0,因此可得李氏函數(shù)的一階導數(shù)為負數(shù),故系統(tǒng)在極值處的平衡狀態(tài)是在其鄰域范圍內(nèi)是一致漸進穩(wěn)定的。積分系數(shù)決定了每個控制周期Δid的大小。圖7表示在?Te/?β=0的臨近區(qū)間內(nèi)不同積分系數(shù)對應(yīng)的系統(tǒng)漸進穩(wěn)定和不穩(wěn)定2種情況,積分系數(shù)對MTPA控制的局部穩(wěn)定性起到?jīng)Q定性作用。

    圖7 系統(tǒng)局部穩(wěn)定性判定示意圖

    4 仿真驗證

    為驗證所提出虛擬信號注入法的MTPA控制策略的可行性,在Mtalab中搭建永磁電機控制仿真系統(tǒng)。首先驗證該策略對MTPA工作點的追蹤可靠性,同時通過對比虛擬正弦信號注入法與傳統(tǒng)的虛擬直流信號法,驗證其在變載運行的動態(tài)響應(yīng)和準確性。仿真采用的電機參數(shù)如表1所示,虛擬直流信號幅值為0.01 A,幅值的選取決定了轉(zhuǎn)矩偏導的準確性,并不會改變收斂到MTPA工作點的速度。

    表1 內(nèi)嵌式永磁電機參數(shù)

    MTPA控制的有效性通過兩方面進行驗證。首先是id=0控制到虛擬信號注入法控制的切換運行,電機轉(zhuǎn)速設(shè)定在1 000 r/min,負載轉(zhuǎn)矩為25 N·m,在第4 s時切換為虛擬信號注入法,仿真結(jié)果如圖8所示。此時電流幅值由5.54 A下降到5.40 A,對應(yīng)的電流矢量角度為5.8°。

    考慮到電機可能工作于頻繁的變負載工況,因此在不同負載轉(zhuǎn)矩工況下對虛擬信號注入法進行驗證。電機在1 000 r/min的轉(zhuǎn)速下,負載轉(zhuǎn)矩的變化由5 N·m先逐步增加到30 N·m,然后再減少至10 N·m,此時的負載變化與定子電流響應(yīng)如圖9所示。

    圖8 由id=0切換至虛擬信號MTPA控制仿真結(jié)果

    圖9 電機變載運行的仿真波形

    為測試所提出的MTPA控制策略的動態(tài)性能,將虛擬正弦信號注入法用于電機的控制,同時將速度環(huán)的PI參數(shù)以及決定直軸電流給定的積分器參數(shù)設(shè)置為一致,對比2種控制策略分別在負載轉(zhuǎn)矩由15 N·m變化到25 N·m,再改變到20 N·m的動態(tài)響應(yīng),其結(jié)果如圖10所示。分析結(jié)果可以得到,虛擬正弦信號法的控制策略由于MTPA判據(jù)依賴于帶通濾波器與低通濾波器,直軸電流給定值需要11 s的時間才能達到穩(wěn)定,而本文所提出的控制策略可以在1.5 s時達到穩(wěn)定,具有更好的動態(tài)性能。

    圖10 虛擬正弦信號注入與虛擬直流信號注入控制波形

    虛擬直流信號注入法的動態(tài)性能取決于積分器參數(shù)與電機工況。電機的工作條件不同將導致該時刻所提取的偏導系數(shù)不同,考慮到電機工況的不可控性與電機參數(shù)的交變的影響,積分器不僅在決定了控制的穩(wěn)定性,同時也影響了控制的動態(tài)特性。改變積分器的系數(shù),得到的id電流響應(yīng)如圖11所示,在不同的積分系數(shù)情況下,系統(tǒng)都能在2 s內(nèi)達到穩(wěn)定,同時可以得到最適合當下工況的參數(shù)設(shè)定。

    圖11 不同積分器設(shè)定的d軸電流響應(yīng)

    為驗證本文的MTPA控制策略的準確性,將傳統(tǒng)的虛擬直流信號注入法運用于本文采用的電機中,在變載的運行的穩(wěn)態(tài)條件下,定量分析對應(yīng)的定子電流,利用Hamilton插值的方法得到2種方案的MTPA軌跡,其結(jié)果如圖12所示。為了定量分析傳統(tǒng)虛擬直流信號注入法的理論誤差,忽略電機電磁參數(shù)交變,根據(jù)式(8)推出系統(tǒng)在不同的直流分量與不同負載下的誤差水平,結(jié)果如圖13。MTPA的誤差隨著虛擬直流信號的增大而增大,負載轉(zhuǎn)矩的越小誤差水平越大。傳統(tǒng)的虛擬直流信號與實際的MTPA在理論上存在細微的誤差,而所提出的MTPA控制方法更加準確。

    圖12 虛擬信號注入法的MTPA控制軌跡

    圖13 傳統(tǒng)虛擬直流信號?Te/?β誤差結(jié)果

    5 結(jié)論

    針對IPMSM當前的虛擬信號注入法的MTPA控制存在的問題,提出了一種改進的虛擬直流信號注入MTPA控制策略,彌補了一般虛擬直流信號注入法的誤差,設(shè)定的自適應(yīng)律在保證最優(yōu)工作點跟蹤準確性的同時,簡化了?Te/?β的計算。通過對虛擬直流信號注入法在極值處的穩(wěn)態(tài)分析可以得出,虛擬直流信號注入法可以避免真實信號注入對電磁轉(zhuǎn)矩產(chǎn)生的擾動,通過對獲取直軸電流參考值的改進,有效減少系統(tǒng)在極值處的轉(zhuǎn)矩變化,驗證了控制策略在最大轉(zhuǎn)矩電流比極值處的局部穩(wěn)定性。idref獲取中積分器的設(shè)置決定了控制系統(tǒng)的穩(wěn)定性與動態(tài)性能。虛擬信號的幅值大小決定了?Te/?β的精度并不會影響系統(tǒng)的動態(tài)性能。仿真結(jié)果也驗證了所提控制方法的可行性。

    猜你喜歡
    矢量轉(zhuǎn)矩幅值
    矢量三角形法的應(yīng)用
    卷取機轉(zhuǎn)矩控制技術(shù)優(yōu)化卷形
    四川冶金(2018年1期)2018-09-25 02:39:26
    基于S變換的交流電網(wǎng)幅值檢測系統(tǒng)計算機仿真研究
    電子制作(2017年7期)2017-06-05 09:36:13
    基于矢量最優(yōu)估計的穩(wěn)健測向方法
    容錯逆變器直接轉(zhuǎn)矩控制策略
    正序電壓幅值檢測及諧波抑制的改進
    三角形法則在動態(tài)平衡問題中的應(yīng)用
    低壓電力線信道脈沖噪聲的幅值與寬度特征
    基于零序電壓幅值增量的消弧線圈調(diào)諧新方法
    電測與儀表(2015年7期)2015-04-09 11:40:10
    基于分級變頻的高轉(zhuǎn)矩軟起動器
    www日本黄色视频网| 人妻丰满熟妇av一区二区三区| 国产精品久久电影中文字幕| 国产精品不卡视频一区二区| av中文乱码字幕在线| 波多野结衣高清作品| 性插视频无遮挡在线免费观看| 精品无人区乱码1区二区| 国产大屁股一区二区在线视频| 国产精品乱码一区二三区的特点| 97人妻精品一区二区三区麻豆| 国产高清视频在线观看网站| 国产精品女同一区二区软件| 日本免费a在线| 成熟少妇高潮喷水视频| 人妻少妇偷人精品九色| 色吧在线观看| 亚洲av.av天堂| 亚洲av.av天堂| 国产精品一及| 亚洲18禁久久av| 久久精品久久久久久噜噜老黄 | 免费人成视频x8x8入口观看| 丰满人妻一区二区三区视频av| 91av网一区二区| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av.av天堂| 黑人高潮一二区| 18禁在线播放成人免费| 亚洲婷婷狠狠爱综合网| 精品日产1卡2卡| 日本熟妇午夜| 一a级毛片在线观看| 搡老妇女老女人老熟妇| 成熟少妇高潮喷水视频| av天堂中文字幕网| 99热网站在线观看| 老司机影院成人| av在线老鸭窝| 国产女主播在线喷水免费视频网站 | 日韩精品中文字幕看吧| 国产伦精品一区二区三区视频9| 日韩欧美免费精品| 欧美成人精品欧美一级黄| 国产精品不卡视频一区二区| 人人妻,人人澡人人爽秒播| 九色成人免费人妻av| 91久久精品国产一区二区三区| 亚洲av免费高清在线观看| 淫妇啪啪啪对白视频| 国产极品精品免费视频能看的| 亚洲av美国av| 乱系列少妇在线播放| 成人欧美大片| 久久久午夜欧美精品| 狂野欧美激情性xxxx在线观看| 狂野欧美白嫩少妇大欣赏| 国产激情偷乱视频一区二区| 人妻制服诱惑在线中文字幕| 日韩欧美国产在线观看| 极品教师在线视频| 国产黄a三级三级三级人| 亚洲真实伦在线观看| 少妇人妻一区二区三区视频| 99久久中文字幕三级久久日本| 天堂动漫精品| 俄罗斯特黄特色一大片| 99久久中文字幕三级久久日本| 99国产极品粉嫩在线观看| 国产精品电影一区二区三区| 欧美潮喷喷水| 免费无遮挡裸体视频| 国产高潮美女av| 国产一区二区三区av在线 | 一a级毛片在线观看| 免费观看人在逋| 特大巨黑吊av在线直播| 国产午夜福利久久久久久| 国产亚洲91精品色在线| 桃色一区二区三区在线观看| 亚洲,欧美,日韩| 搡老岳熟女国产| 岛国在线免费视频观看| 国产aⅴ精品一区二区三区波| 日本-黄色视频高清免费观看| 国产成年人精品一区二区| 男人和女人高潮做爰伦理| 国产麻豆成人av免费视频| 日韩高清综合在线| 国产综合懂色| 久久热精品热| 天美传媒精品一区二区| 国产精品久久电影中文字幕| 欧美区成人在线视频| 亚州av有码| 麻豆成人午夜福利视频| 国产午夜精品论理片| 搡老岳熟女国产| 亚洲国产精品成人综合色| 久久精品国产99精品国产亚洲性色| 久久久久国产网址| 亚洲自偷自拍三级| 熟女人妻精品中文字幕| 日韩,欧美,国产一区二区三区 | 久久欧美精品欧美久久欧美| 亚洲欧美日韩卡通动漫| 男女那种视频在线观看| 久久精品人妻少妇| 午夜亚洲福利在线播放| 中文字幕人妻熟人妻熟丝袜美| 三级国产精品欧美在线观看| 我要搜黄色片| 色尼玛亚洲综合影院| 日本黄色视频三级网站网址| 最近视频中文字幕2019在线8| 国产69精品久久久久777片| 人妻制服诱惑在线中文字幕| 亚洲国产欧洲综合997久久,| 最新在线观看一区二区三区| 女同久久另类99精品国产91| 综合色av麻豆| 国产久久久一区二区三区| 成年女人看的毛片在线观看| 欧美性猛交黑人性爽| 亚洲成a人片在线一区二区| 日韩国内少妇激情av| 久久久成人免费电影| 亚洲精品国产av成人精品 | 日韩成人av中文字幕在线观看 | av天堂中文字幕网| av视频在线观看入口| 黄片wwwwww| av天堂中文字幕网| 日韩欧美国产在线观看| 国产精品一区二区免费欧美| 观看免费一级毛片| 日本一本二区三区精品| 免费无遮挡裸体视频| 国产伦一二天堂av在线观看| 熟女人妻精品中文字幕| 国产片特级美女逼逼视频| 久久久色成人| 午夜久久久久精精品| 国产白丝娇喘喷水9色精品| 日韩人妻高清精品专区| 亚洲熟妇中文字幕五十中出| 毛片女人毛片| 欧美国产日韩亚洲一区| 十八禁国产超污无遮挡网站| 亚洲av成人精品一区久久| 亚洲人成网站高清观看| 18禁在线播放成人免费| 听说在线观看完整版免费高清| 天天躁日日操中文字幕| 国产精品三级大全| 亚洲欧美成人综合另类久久久 | 日韩欧美三级三区| 久久精品国产亚洲av香蕉五月| 一个人免费在线观看电影| 91久久精品国产一区二区成人| 色综合亚洲欧美另类图片| 深爱激情五月婷婷| 国产精品久久视频播放| 一级毛片久久久久久久久女| 天堂√8在线中文| 一级毛片电影观看 | 日韩精品青青久久久久久| 成人永久免费在线观看视频| 亚洲欧美日韩东京热| 黄色一级大片看看| 成人永久免费在线观看视频| 国产伦在线观看视频一区| 人妻制服诱惑在线中文字幕| 日本黄色视频三级网站网址| 日韩精品青青久久久久久| 男女啪啪激烈高潮av片| 日本五十路高清| av国产免费在线观看| 黑人高潮一二区| 国产精品三级大全| 亚洲国产欧洲综合997久久,| 丰满乱子伦码专区| 欧美性猛交╳xxx乱大交人| 国产成人a∨麻豆精品| 免费观看的影片在线观看| 国产美女午夜福利| 精品久久久久久成人av| 欧美一级a爱片免费观看看| 久久亚洲国产成人精品v| 国产伦在线观看视频一区| 熟女电影av网| 黄片wwwwww| 亚洲一级一片aⅴ在线观看| 亚洲熟妇中文字幕五十中出| 美女黄网站色视频| 欧美一区二区国产精品久久精品| 免费人成视频x8x8入口观看| 亚洲国产色片| 干丝袜人妻中文字幕| 又爽又黄a免费视频| 国产在线男女| 国产精品乱码一区二三区的特点| 日本黄大片高清| 国产一区二区在线观看日韩| 一区二区三区免费毛片| 桃色一区二区三区在线观看| 免费看光身美女| 国产精品久久久久久精品电影| 春色校园在线视频观看| 99热精品在线国产| 小说图片视频综合网站| 神马国产精品三级电影在线观看| 国产一级毛片七仙女欲春2| 97碰自拍视频| 热99在线观看视频| 久久久久久国产a免费观看| 欧美激情在线99| 少妇熟女欧美另类| 国产久久久一区二区三区| 国产高清不卡午夜福利| 免费av毛片视频| 午夜亚洲福利在线播放| 性色avwww在线观看| 国产日本99.免费观看| 国产av一区在线观看免费| 免费人成在线观看视频色| 人妻制服诱惑在线中文字幕| 精品人妻偷拍中文字幕| 自拍偷自拍亚洲精品老妇| АⅤ资源中文在线天堂| 免费人成视频x8x8入口观看| 51国产日韩欧美| 直男gayav资源| 午夜免费男女啪啪视频观看 | 自拍偷自拍亚洲精品老妇| 99热只有精品国产| 精品免费久久久久久久清纯| 老司机影院成人| 色av中文字幕| 久久精品91蜜桃| 色哟哟·www| 欧美极品一区二区三区四区| 国产伦在线观看视频一区| 一区二区三区高清视频在线| 国产精品无大码| 有码 亚洲区| 我要搜黄色片| 韩国av在线不卡| 国产精品一区二区免费欧美| 国产精品国产高清国产av| 男女啪啪激烈高潮av片| 99热只有精品国产| 欧美激情在线99| 国产视频内射| 91午夜精品亚洲一区二区三区| 国产一区二区激情短视频| 最近2019中文字幕mv第一页| 精品久久久久久久久久免费视频| 伊人久久精品亚洲午夜| 男人舔奶头视频| 久久人人爽人人爽人人片va| 亚洲精品久久国产高清桃花| 如何舔出高潮| 欧美人与善性xxx| 亚洲激情五月婷婷啪啪| 国产精品人妻久久久影院| 国产成人a∨麻豆精品| 欧美中文日本在线观看视频| 美女被艹到高潮喷水动态| 日本黄大片高清| 九九久久精品国产亚洲av麻豆| 女人十人毛片免费观看3o分钟| 国产日本99.免费观看| 成人欧美大片| 国产熟女欧美一区二区| 色噜噜av男人的天堂激情| 免费大片18禁| 亚洲国产色片| av黄色大香蕉| 久久精品国产亚洲av涩爱 | 22中文网久久字幕| 深夜精品福利| 色播亚洲综合网| 亚洲一区高清亚洲精品| 人人妻人人看人人澡| 91久久精品电影网| 国产精品免费一区二区三区在线| 女人十人毛片免费观看3o分钟| 97超级碰碰碰精品色视频在线观看| 亚洲天堂国产精品一区在线| 婷婷亚洲欧美| 亚洲欧美日韩无卡精品| 久久久久久大精品| 久久国内精品自在自线图片| 亚洲久久久久久中文字幕| 天天一区二区日本电影三级| 欧美日本亚洲视频在线播放| 亚洲无线在线观看| 99在线视频只有这里精品首页| 久久久a久久爽久久v久久| 婷婷精品国产亚洲av| 精品人妻视频免费看| 深爱激情五月婷婷| 久久久久性生活片| av在线老鸭窝| 国产亚洲精品久久久久久毛片| 成人三级黄色视频| 舔av片在线| 日韩亚洲欧美综合| 精品人妻一区二区三区麻豆 | 亚洲欧美精品自产自拍| 99热网站在线观看| 一级毛片电影观看 | 日日摸夜夜添夜夜添小说| 亚洲欧美精品自产自拍| 美女黄网站色视频| 高清毛片免费观看视频网站| 亚洲人成网站在线播| 欧美激情久久久久久爽电影| 18禁黄网站禁片免费观看直播| 乱人视频在线观看| 麻豆成人午夜福利视频| 又黄又爽又刺激的免费视频.| 小说图片视频综合网站| 黄片wwwwww| 国产一区二区激情短视频| 成人无遮挡网站| 久久精品国产99精品国产亚洲性色| 亚洲欧美日韩高清在线视频| 亚洲av五月六月丁香网| 看非洲黑人一级黄片| 一夜夜www| 亚洲一区高清亚洲精品| 在线播放无遮挡| 免费无遮挡裸体视频| 国产精品一区二区三区四区免费观看 | 免费看美女性在线毛片视频| 精品久久久久久久久av| 精品乱码久久久久久99久播| 丰满的人妻完整版| 亚洲美女视频黄频| 一夜夜www| 国模一区二区三区四区视频| 国产精品女同一区二区软件| 菩萨蛮人人尽说江南好唐韦庄 | 露出奶头的视频| 亚洲精品久久国产高清桃花| 国产精品久久电影中文字幕| 白带黄色成豆腐渣| 搡老岳熟女国产| 免费观看人在逋| 成年女人毛片免费观看观看9| 色5月婷婷丁香| 深夜a级毛片| 插阴视频在线观看视频| 免费人成在线观看视频色| 村上凉子中文字幕在线| 久久久国产成人免费| 大香蕉久久网| a级毛片免费高清观看在线播放| 国产v大片淫在线免费观看| 午夜福利在线观看免费完整高清在 | 国产综合懂色| 日韩欧美 国产精品| 卡戴珊不雅视频在线播放| 99久国产av精品| 国产精品人妻久久久久久| 日韩精品青青久久久久久| 97热精品久久久久久| 午夜视频国产福利| 国产不卡一卡二| 麻豆成人午夜福利视频| 欧美激情国产日韩精品一区| 蜜臀久久99精品久久宅男| 欧美一区二区国产精品久久精品| 可以在线观看的亚洲视频| 久久久精品大字幕| 一级黄色大片毛片| 久久这里只有精品中国| 欧美成人a在线观看| 国产精品伦人一区二区| 午夜日韩欧美国产| 亚洲电影在线观看av| 中国国产av一级| av在线老鸭窝| 国产精品福利在线免费观看| 老司机影院成人| 伦理电影大哥的女人| 欧美一区二区国产精品久久精品| 级片在线观看| 午夜影院日韩av| 嫩草影院入口| 精品一区二区三区av网在线观看| 国产男靠女视频免费网站| 免费观看的影片在线观看| 欧美绝顶高潮抽搐喷水| 久久久久久久久久黄片| 成人性生交大片免费视频hd| 日韩成人伦理影院| 亚洲av.av天堂| 亚洲人与动物交配视频| 日韩国内少妇激情av| 美女黄网站色视频| 51国产日韩欧美| 亚洲国产精品合色在线| 国产一区亚洲一区在线观看| 成人二区视频| 久久久午夜欧美精品| 欧美色视频一区免费| 又黄又爽又免费观看的视频| 成年女人永久免费观看视频| 久久99热6这里只有精品| 三级男女做爰猛烈吃奶摸视频| 露出奶头的视频| 成年av动漫网址| 国产 一区精品| 乱系列少妇在线播放| 免费av不卡在线播放| 日本免费a在线| a级一级毛片免费在线观看| 精品一区二区免费观看| 国产亚洲精品久久久com| 国产毛片a区久久久久| 亚洲精品色激情综合| 国产精品美女特级片免费视频播放器| 97超碰精品成人国产| 免费av毛片视频| aaaaa片日本免费| 亚洲av二区三区四区| 哪里可以看免费的av片| 欧美成人免费av一区二区三区| 日韩精品青青久久久久久| 久久6这里有精品| 成人综合一区亚洲| 一a级毛片在线观看| 少妇熟女欧美另类| 欧美绝顶高潮抽搐喷水| 精品午夜福利在线看| 精品99又大又爽又粗少妇毛片| 天天躁夜夜躁狠狠久久av| 免费搜索国产男女视频| 免费av观看视频| 国产精品久久久久久精品电影| 秋霞在线观看毛片| 国产熟女欧美一区二区| 中文字幕熟女人妻在线| 老师上课跳d突然被开到最大视频| 深夜a级毛片| 久久精品夜色国产| 久久精品91蜜桃| 一夜夜www| 欧美潮喷喷水| 国产高清激情床上av| 人妻少妇偷人精品九色| 波野结衣二区三区在线| 观看美女的网站| 麻豆一二三区av精品| 国产高清激情床上av| 性欧美人与动物交配| 熟妇人妻久久中文字幕3abv| 18禁黄网站禁片免费观看直播| 男人狂女人下面高潮的视频| 中文字幕av成人在线电影| 噜噜噜噜噜久久久久久91| 又爽又黄无遮挡网站| 国产三级中文精品| 乱人视频在线观看| 欧美成人免费av一区二区三区| 色在线成人网| 综合色丁香网| 国产欧美日韩精品一区二区| 国产成年人精品一区二区| 中文字幕免费在线视频6| h日本视频在线播放| 乱人视频在线观看| 丝袜喷水一区| 最新中文字幕久久久久| 干丝袜人妻中文字幕| 国产亚洲精品久久久久久毛片| 超碰av人人做人人爽久久| 成人综合一区亚洲| 久久这里只有精品中国| 成人鲁丝片一二三区免费| 日韩欧美 国产精品| 97热精品久久久久久| 一个人看视频在线观看www免费| 亚州av有码| 尤物成人国产欧美一区二区三区| 国产 一区精品| 日韩 亚洲 欧美在线| 在线观看66精品国产| 国产男人的电影天堂91| 日本在线视频免费播放| 国产真实伦视频高清在线观看| 人妻少妇偷人精品九色| 精品一区二区三区人妻视频| 亚洲aⅴ乱码一区二区在线播放| 小蜜桃在线观看免费完整版高清| 免费不卡的大黄色大毛片视频在线观看 | 国产av不卡久久| 国产欧美日韩精品亚洲av| 别揉我奶头 嗯啊视频| 国产精品久久久久久精品电影| 秋霞在线观看毛片| 自拍偷自拍亚洲精品老妇| 最近2019中文字幕mv第一页| 此物有八面人人有两片| 久久人人爽人人爽人人片va| 美女被艹到高潮喷水动态| 欧美性猛交黑人性爽| 国产亚洲精品av在线| 我要搜黄色片| 欧美人与善性xxx| 一个人免费在线观看电影| 亚洲在线自拍视频| 99在线人妻在线中文字幕| 直男gayav资源| 久久精品国产清高在天天线| 久久人妻av系列| 高清日韩中文字幕在线| 亚洲在线自拍视频| 99久久精品热视频| 九九爱精品视频在线观看| 色综合亚洲欧美另类图片| 久久人妻av系列| 亚洲一级一片aⅴ在线观看| 亚洲图色成人| 一区福利在线观看| 中文字幕免费在线视频6| 久久亚洲精品不卡| 最近最新中文字幕大全电影3| 色噜噜av男人的天堂激情| 日本在线视频免费播放| 亚洲自拍偷在线| 少妇熟女欧美另类| 人妻夜夜爽99麻豆av| 51国产日韩欧美| 97人妻精品一区二区三区麻豆| 啦啦啦啦在线视频资源| 又粗又爽又猛毛片免费看| 国产欧美日韩一区二区精品| 日韩一本色道免费dvd| 国产黄色视频一区二区在线观看 | 精品乱码久久久久久99久播| 搡老妇女老女人老熟妇| 欧美成人精品欧美一级黄| 国产一区二区在线观看日韩| 欧美色欧美亚洲另类二区| 欧美极品一区二区三区四区| 日韩中字成人| 亚洲不卡免费看| 少妇丰满av| 特大巨黑吊av在线直播| 午夜福利18| 1024手机看黄色片| 成年女人看的毛片在线观看| 中国美女看黄片| 大香蕉久久网| 国产高清有码在线观看视频| 午夜日韩欧美国产| 97在线视频观看| 国产精品电影一区二区三区| 哪里可以看免费的av片| 99热这里只有是精品50| 国产精品久久久久久精品电影| 欧美国产日韩亚洲一区| 国产私拍福利视频在线观看| 免费人成在线观看视频色| 狂野欧美激情性xxxx在线观看| 免费看光身美女| 亚洲高清免费不卡视频| 国产一级毛片七仙女欲春2| 久久精品夜色国产| 无遮挡黄片免费观看| 如何舔出高潮| 国产精品,欧美在线| 精品一区二区三区av网在线观看| 国产一区二区激情短视频| 久久亚洲国产成人精品v| 男女之事视频高清在线观看| 白带黄色成豆腐渣| 18禁在线无遮挡免费观看视频 | 99在线人妻在线中文字幕| 午夜视频国产福利| 狠狠狠狠99中文字幕| 国语自产精品视频在线第100页| 欧美成人一区二区免费高清观看| av天堂中文字幕网| 老司机影院成人| 午夜福利18| 欧美国产日韩亚洲一区| 99热只有精品国产| 啦啦啦韩国在线观看视频| 色哟哟·www| 国产精品不卡视频一区二区| 深爱激情五月婷婷| 亚洲欧美日韩东京热| 久久精品综合一区二区三区| 国产伦一二天堂av在线观看| 日韩成人av中文字幕在线观看 | 欧美zozozo另类| 你懂的网址亚洲精品在线观看 | 日韩欧美三级三区| 久久人人爽人人片av| 身体一侧抽搐| 国产成人影院久久av| 亚洲成人久久性| 中文字幕精品亚洲无线码一区| 国产午夜福利久久久久久| 最近的中文字幕免费完整| 精品人妻偷拍中文字幕| 少妇熟女欧美另类| 美女大奶头视频| 色av中文字幕| 香蕉av资源在线|