• <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)矩軟起動器
    久久精品国产清高在天天线| 一a级毛片在线观看| 亚洲国产欧洲综合997久久, | 天堂影院成人在线观看| 搡老熟女国产l中国老女人| 亚洲中文字幕日韩| 国产又色又爽无遮挡免费看| 免费看美女性在线毛片视频| 非洲黑人性xxxx精品又粗又长| 精品午夜福利视频在线观看一区| 好看av亚洲va欧美ⅴa在| 麻豆久久精品国产亚洲av| 精品国产一区二区三区四区第35| a级毛片a级免费在线| 国产精品久久久久久人妻精品电影| 日韩有码中文字幕| 夜夜躁狠狠躁天天躁| 欧美绝顶高潮抽搐喷水| 日韩成人在线观看一区二区三区| 国产精品爽爽va在线观看网站 | 久久久久国产一级毛片高清牌| 婷婷丁香在线五月| av片东京热男人的天堂| 亚洲色图 男人天堂 中文字幕| 淫秽高清视频在线观看| 一进一出好大好爽视频| 神马国产精品三级电影在线观看 | 欧美日本视频| 中文字幕高清在线视频| 日日夜夜操网爽| 日本撒尿小便嘘嘘汇集6| 久久人人精品亚洲av| 午夜免费成人在线视频| 久久国产乱子伦精品免费另类| av在线播放免费不卡| 国产精品免费一区二区三区在线| 国产精品一区二区三区四区久久 | 久久天堂一区二区三区四区| 久久香蕉国产精品| 黑丝袜美女国产一区| 免费看美女性在线毛片视频| 日本a在线网址| 欧美激情高清一区二区三区| 国内少妇人妻偷人精品xxx网站 | 久久人妻福利社区极品人妻图片| 老熟妇仑乱视频hdxx| 91在线观看av| 国产又黄又爽又无遮挡在线| 免费看日本二区| 国产亚洲欧美精品永久| 十分钟在线观看高清视频www| 青草久久国产| 亚洲av熟女| 成人免费观看视频高清| or卡值多少钱| 久久久国产精品麻豆| 欧美性猛交╳xxx乱大交人| 免费电影在线观看免费观看| 久久久久久久久免费视频了| 熟女少妇亚洲综合色aaa.| 国产在线精品亚洲第一网站| 久久久久久久久中文| 欧美性长视频在线观看| 欧美 亚洲 国产 日韩一| 亚洲自偷自拍图片 自拍| 欧美激情久久久久久爽电影| 一级黄色大片毛片| 日韩欧美 国产精品| 波多野结衣高清无吗| 村上凉子中文字幕在线| 午夜福利欧美成人| avwww免费| 宅男免费午夜| 亚洲精品美女久久av网站| 亚洲无线在线观看| 免费女性裸体啪啪无遮挡网站| 亚洲av第一区精品v没综合| 欧美成人性av电影在线观看| 成人午夜高清在线视频 | 法律面前人人平等表现在哪些方面| 两人在一起打扑克的视频| 最近在线观看免费完整版| 后天国语完整版免费观看| 麻豆国产av国片精品| 久久九九热精品免费| 久久性视频一级片| 日本a在线网址| 免费看十八禁软件| 黄片大片在线免费观看| 亚洲精华国产精华精| 欧美性猛交╳xxx乱大交人| 日韩欧美三级三区| 一进一出好大好爽视频| 久久精品国产亚洲av香蕉五月| 亚洲自拍偷在线| 久久欧美精品欧美久久欧美| 一进一出抽搐gif免费好疼| 亚洲精品国产区一区二| 精品少妇一区二区三区视频日本电影| 欧美 亚洲 国产 日韩一| 好看av亚洲va欧美ⅴa在| 国产伦一二天堂av在线观看| 中文字幕高清在线视频| 午夜成年电影在线免费观看| 夜夜夜夜夜久久久久| 亚洲一区高清亚洲精品| 免费在线观看亚洲国产| 亚洲色图 男人天堂 中文字幕| 国产精品 国内视频| 满18在线观看网站| 少妇裸体淫交视频免费看高清 | 黄色 视频免费看| 精品不卡国产一区二区三区| 老熟妇乱子伦视频在线观看| 国产亚洲精品一区二区www| 免费在线观看日本一区| 成人18禁高潮啪啪吃奶动态图| 日韩欧美国产在线观看| 精品久久蜜臀av无| 亚洲免费av在线视频| 国产精品美女特级片免费视频播放器 | 亚洲熟妇熟女久久| 一二三四社区在线视频社区8| 精品国产美女av久久久久小说| 亚洲精品美女久久av网站| 性色av乱码一区二区三区2| 一本久久中文字幕| 最近在线观看免费完整版| 看黄色毛片网站| 午夜免费激情av| cao死你这个sao货| 大型av网站在线播放| 午夜亚洲福利在线播放| 国产一区二区三区在线臀色熟女| 成人亚洲精品一区在线观看| 久久久国产精品麻豆| 色综合站精品国产| 亚洲av日韩精品久久久久久密| 久久香蕉国产精品| 99在线人妻在线中文字幕| 99精品欧美一区二区三区四区| 亚洲久久久国产精品| 免费电影在线观看免费观看| 波多野结衣高清作品| 国产一卡二卡三卡精品| 久久精品夜夜夜夜夜久久蜜豆 | 国产国语露脸激情在线看| 久久亚洲真实| 两性夫妻黄色片| x7x7x7水蜜桃| 精品久久蜜臀av无| 免费在线观看成人毛片| 一本久久中文字幕| 日日干狠狠操夜夜爽| 一个人免费在线观看的高清视频| 一二三四社区在线视频社区8| 男女下面进入的视频免费午夜 | 美女大奶头视频| 欧美在线黄色| 日韩精品免费视频一区二区三区| 久久久久亚洲av毛片大全| 欧美亚洲日本最大视频资源| 悠悠久久av| 亚洲精品中文字幕在线视频| 国产亚洲精品久久久久久毛片| 国产精品一区二区精品视频观看| 岛国在线观看网站| 两性午夜刺激爽爽歪歪视频在线观看 | 在线国产一区二区在线| 一卡2卡三卡四卡精品乱码亚洲| 国产v大片淫在线免费观看| 91麻豆av在线| 久久伊人香网站| 午夜福利成人在线免费观看| 精品日产1卡2卡| 国产成人欧美在线观看| 精品久久蜜臀av无| 夜夜夜夜夜久久久久| www国产在线视频色| 亚洲成国产人片在线观看| 少妇的丰满在线观看| 99久久国产精品久久久| 三级毛片av免费| 好男人在线观看高清免费视频 | 9191精品国产免费久久| 成人永久免费在线观看视频| 国产伦在线观看视频一区| 亚洲国产欧美网| 久久精品国产综合久久久| 中文字幕高清在线视频| 亚洲av中文字字幕乱码综合 | 国产成人系列免费观看| av中文乱码字幕在线| 少妇被粗大的猛进出69影院| 亚洲精华国产精华精| 国产精品永久免费网站| 精品福利观看| 中文字幕人成人乱码亚洲影| 最近最新中文字幕大全免费视频| 91麻豆av在线| 日韩精品中文字幕看吧| 中文字幕人成人乱码亚洲影| 亚洲精品中文字幕一二三四区| 成年女人毛片免费观看观看9| 亚洲国产日韩欧美精品在线观看 | 国产成年人精品一区二区| 免费看十八禁软件| 制服人妻中文乱码| 黑人巨大精品欧美一区二区mp4| 久久婷婷人人爽人人干人人爱| 99久久无色码亚洲精品果冻| 熟女少妇亚洲综合色aaa.| 欧美不卡视频在线免费观看 | av在线播放免费不卡| 国产精品久久久久久人妻精品电影| 99热这里只有精品一区 | 国产黄色小视频在线观看| 成年女人毛片免费观看观看9| 国产伦一二天堂av在线观看| 国产一区二区三区视频了| 美女午夜性视频免费| 久久精品国产亚洲av高清一级| 亚洲av电影不卡..在线观看| 婷婷精品国产亚洲av| 18禁国产床啪视频网站| 欧美乱妇无乱码| 啦啦啦观看免费观看视频高清| 国产精品综合久久久久久久免费| 视频区欧美日本亚洲| 日本撒尿小便嘘嘘汇集6| 中文字幕精品免费在线观看视频| 啦啦啦免费观看视频1| 亚洲欧美一区二区三区黑人| 成人三级做爰电影| 欧美激情久久久久久爽电影| 国产蜜桃级精品一区二区三区| 久久精品91无色码中文字幕| 欧美久久黑人一区二区| 亚洲国产欧美日韩在线播放| 国产黄a三级三级三级人| 午夜老司机福利片| 热re99久久国产66热| 亚洲精品在线观看二区| 国产精品av久久久久免费| 国产v大片淫在线免费观看| 婷婷亚洲欧美| 午夜免费观看网址| 国产亚洲av嫩草精品影院| 免费在线观看日本一区| 大型黄色视频在线免费观看| 亚洲性夜色夜夜综合| 欧美黄色淫秽网站| www日本在线高清视频| 国产又色又爽无遮挡免费看| 久久精品成人免费网站| 人人澡人人妻人| 老熟妇仑乱视频hdxx| 亚洲狠狠婷婷综合久久图片| 十分钟在线观看高清视频www| 国产精品乱码一区二三区的特点| 亚洲第一av免费看| 淫秽高清视频在线观看| av超薄肉色丝袜交足视频| av福利片在线| 久久草成人影院| 免费在线观看日本一区| 欧美日韩瑟瑟在线播放| 亚洲一区二区三区色噜噜| 无遮挡黄片免费观看| 久久中文看片网| 国产亚洲精品久久久久5区| 两个人视频免费观看高清| 真人一进一出gif抽搐免费| 女人被狂操c到高潮| 国内精品久久久久精免费| 亚洲一区中文字幕在线| 亚洲专区中文字幕在线| 久久亚洲精品不卡| 欧美激情极品国产一区二区三区| cao死你这个sao货| 日韩有码中文字幕| 高清在线国产一区| 最近在线观看免费完整版| 欧美黄色淫秽网站| 两个人看的免费小视频| 亚洲自偷自拍图片 自拍| 国产精品综合久久久久久久免费| 亚洲av日韩精品久久久久久密| 久久这里只有精品19| 国产激情欧美一区二区| 国产男靠女视频免费网站| 黄频高清免费视频| 国产成人啪精品午夜网站| 激情在线观看视频在线高清| 婷婷六月久久综合丁香| 给我免费播放毛片高清在线观看| 91老司机精品| 999久久久精品免费观看国产| 哪里可以看免费的av片| 大型av网站在线播放| 久久久久国产一级毛片高清牌| 欧美 亚洲 国产 日韩一| 国产男靠女视频免费网站| 国产成年人精品一区二区| 国产成人精品久久二区二区91| 1024手机看黄色片| 午夜成年电影在线免费观看| 成人精品一区二区免费| 午夜两性在线视频| 国产精品久久久av美女十八| 久久香蕉国产精品| 大香蕉久久成人网| 99久久综合精品五月天人人| 国产亚洲精品一区二区www| 侵犯人妻中文字幕一二三四区| 久久人人精品亚洲av| 国产av不卡久久| 视频在线观看一区二区三区| e午夜精品久久久久久久| 国产99久久九九免费精品| 在线观看免费视频日本深夜| 午夜福利18| 午夜福利免费观看在线| 俺也久久电影网| 免费看日本二区| 亚洲午夜理论影院| 老熟妇仑乱视频hdxx| 一进一出好大好爽视频| 免费搜索国产男女视频| 十分钟在线观看高清视频www| 在线观看日韩欧美| 俺也久久电影网| 日本三级黄在线观看| 欧美乱色亚洲激情| 免费在线观看影片大全网站| 久久精品国产亚洲av香蕉五月| 免费观看人在逋| 美女 人体艺术 gogo| 脱女人内裤的视频| 亚洲熟女毛片儿| 国产三级在线视频| 亚洲av成人av| 人成视频在线观看免费观看| 欧美成狂野欧美在线观看| 伦理电影免费视频| 国产亚洲精品一区二区www| 亚洲av片天天在线观看| 久久婷婷成人综合色麻豆| 香蕉av资源在线| 日日干狠狠操夜夜爽| 国产精品久久电影中文字幕| 精品国内亚洲2022精品成人| 少妇 在线观看| 久久狼人影院| 日韩av在线大香蕉| 免费在线观看成人毛片| 真人做人爱边吃奶动态| 亚洲精品在线美女| 日韩欧美一区视频在线观看| 国产精品自产拍在线观看55亚洲| 麻豆国产av国片精品| 欧洲精品卡2卡3卡4卡5卡区| 99久久99久久久精品蜜桃| 精品电影一区二区在线| 欧美黄色片欧美黄色片| 国产亚洲av嫩草精品影院| 成人18禁在线播放| 亚洲国产精品合色在线| 神马国产精品三级电影在线观看 | 国产伦一二天堂av在线观看| 老司机靠b影院| 亚洲成av片中文字幕在线观看| 国产成人影院久久av| 伦理电影免费视频| 精品福利观看| 欧美黑人精品巨大| cao死你这个sao货| 熟女少妇亚洲综合色aaa.| 超碰成人久久| 99re在线观看精品视频| 两个人视频免费观看高清| www国产在线视频色| 免费在线观看视频国产中文字幕亚洲| 村上凉子中文字幕在线| 亚洲人成77777在线视频| av片东京热男人的天堂| 91九色精品人成在线观看| 久久人人精品亚洲av| 亚洲av第一区精品v没综合| 老司机福利观看| 亚洲熟女毛片儿| 亚洲成人国产一区在线观看| 一边摸一边抽搐一进一小说| 日韩欧美在线二视频| 国产精品乱码一区二三区的特点| 午夜影院日韩av| 成人免费观看视频高清| 精品一区二区三区四区五区乱码| 成人一区二区视频在线观看| 国产爱豆传媒在线观看 | 成人三级黄色视频| 久久中文字幕人妻熟女| 成年免费大片在线观看| 久久久久久久午夜电影| 在线观看免费日韩欧美大片| 制服人妻中文乱码| 熟女电影av网| 国产私拍福利视频在线观看| www.www免费av| 亚洲天堂国产精品一区在线| 国语自产精品视频在线第100页| 国产一区在线观看成人免费| 69av精品久久久久久| 午夜福利成人在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 久久国产乱子伦精品免费另类| 搡老妇女老女人老熟妇| 首页视频小说图片口味搜索| 国产亚洲精品av在线| 男女午夜视频在线观看| 亚洲人成77777在线视频| 日韩一卡2卡3卡4卡2021年| 亚洲自偷自拍图片 自拍| 午夜a级毛片| 久久天躁狠狠躁夜夜2o2o| 久久伊人香网站| 日韩中文字幕欧美一区二区| 两个人免费观看高清视频| 99riav亚洲国产免费| 人人澡人人妻人| 免费看a级黄色片| 中国美女看黄片| 中文资源天堂在线| 变态另类丝袜制服| 老司机午夜十八禁免费视频| 夜夜爽天天搞| 色播亚洲综合网| 成人18禁在线播放| 日韩中文字幕欧美一区二区| 哪里可以看免费的av片| 99久久国产精品久久久| 香蕉久久夜色| 国产蜜桃级精品一区二区三区| 欧美成狂野欧美在线观看| 丁香六月欧美| 亚洲精品在线观看二区| 国产成人精品无人区| 国产精品久久电影中文字幕| 精品国产超薄肉色丝袜足j| 哪里可以看免费的av片| e午夜精品久久久久久久| 亚洲午夜精品一区,二区,三区| 国产蜜桃级精品一区二区三区| 亚洲五月婷婷丁香| 亚洲精品av麻豆狂野| 日韩精品免费视频一区二区三区| 亚洲国产精品成人综合色| 怎么达到女性高潮| 婷婷精品国产亚洲av在线| cao死你这个sao货| 淫妇啪啪啪对白视频| av天堂在线播放| 久久亚洲真实| 精品国内亚洲2022精品成人| 精品人妻1区二区| 欧美黑人巨大hd| 99精品在免费线老司机午夜| 亚洲自拍偷在线| 欧美色视频一区免费| 香蕉久久夜色| 亚洲激情在线av| 久久这里只有精品19| 国产精品99久久99久久久不卡| 国产黄a三级三级三级人| 色综合欧美亚洲国产小说| 亚洲中文字幕日韩| 国产精品精品国产色婷婷| 在线播放国产精品三级| 亚洲av成人不卡在线观看播放网| 日本免费a在线| 中文在线观看免费www的网站 | 黄片播放在线免费| 黄色毛片三级朝国网站| 亚洲片人在线观看| 欧美成人免费av一区二区三区| 99国产极品粉嫩在线观看| 成年人黄色毛片网站| 老熟妇乱子伦视频在线观看| 国产亚洲精品久久久久5区| 国产激情偷乱视频一区二区| 波多野结衣高清作品| 久久午夜综合久久蜜桃| 国产主播在线观看一区二区| 夜夜夜夜夜久久久久| 男人舔女人的私密视频| 好男人在线观看高清免费视频 | 精品免费久久久久久久清纯| 麻豆一二三区av精品| 国产精品二区激情视频| 丰满人妻熟妇乱又伦精品不卡| 久久欧美精品欧美久久欧美| 一级作爱视频免费观看| 亚洲精品在线美女| 中亚洲国语对白在线视频| 757午夜福利合集在线观看| 久久久久久大精品| 久久青草综合色| 亚洲自拍偷在线| a在线观看视频网站| 在线观看舔阴道视频| 国产精品免费一区二区三区在线| 熟女少妇亚洲综合色aaa.| 美女大奶头视频| 十分钟在线观看高清视频www| 国内少妇人妻偷人精品xxx网站 | 亚洲在线自拍视频| 国内揄拍国产精品人妻在线 | 热re99久久国产66热| 中文字幕人妻丝袜一区二区| 亚洲成国产人片在线观看| 精品不卡国产一区二区三区| a级毛片a级免费在线| 欧美久久黑人一区二区| 精品国产一区二区三区四区第35| 国产av又大| 婷婷精品国产亚洲av| 天天添夜夜摸| 国产精品98久久久久久宅男小说| 国产一区二区激情短视频| 欧美另类亚洲清纯唯美| 国产精品久久久久久亚洲av鲁大| a在线观看视频网站| 99国产精品一区二区三区| 亚洲专区国产一区二区| 久久久水蜜桃国产精品网| 熟女少妇亚洲综合色aaa.| 国产精品久久久av美女十八| 亚洲成人久久爱视频| 国产亚洲精品久久久久久毛片| 国产成人精品无人区| 99久久99久久久精品蜜桃| 久久精品影院6| 日韩欧美国产在线观看| 国产又爽黄色视频| 99久久精品国产亚洲精品| 黄色 视频免费看| 久久精品成人免费网站| 韩国av一区二区三区四区| 国产欧美日韩精品亚洲av| 午夜久久久久精精品| 啦啦啦免费观看视频1| 亚洲成a人片在线一区二区| 少妇粗大呻吟视频| 给我免费播放毛片高清在线观看| 777久久人妻少妇嫩草av网站| 欧美性猛交╳xxx乱大交人| 亚洲欧美日韩高清在线视频| 他把我摸到了高潮在线观看| 成人国语在线视频| 午夜免费成人在线视频| 在线观看日韩欧美| 麻豆久久精品国产亚洲av| 一区二区三区精品91| 777久久人妻少妇嫩草av网站| 一夜夜www| 韩国av一区二区三区四区| 久久精品影院6| 亚洲在线自拍视频| 国产免费av片在线观看野外av| 亚洲av电影不卡..在线观看| 2021天堂中文幕一二区在线观 | 欧美日韩精品网址| 午夜亚洲福利在线播放| 亚洲av成人一区二区三| 哪里可以看免费的av片| 黄色a级毛片大全视频| 国产精品,欧美在线| 成人18禁在线播放| 国产成人系列免费观看| 国产成人一区二区三区免费视频网站| 精品一区二区三区视频在线观看免费| www.www免费av| 国产精品免费视频内射| 国产精品 国内视频| 亚洲国产精品合色在线| 亚洲欧美一区二区三区黑人| 国产伦一二天堂av在线观看| 亚洲狠狠婷婷综合久久图片| 久久国产亚洲av麻豆专区| 啦啦啦 在线观看视频| 国产99久久九九免费精品| 波多野结衣高清作品| 男女那种视频在线观看| 草草在线视频免费看| 国产区一区二久久| 欧美日本亚洲视频在线播放| 国产麻豆成人av免费视频| 好男人电影高清在线观看| 国产极品粉嫩免费观看在线| 天天躁夜夜躁狠狠躁躁| 午夜精品久久久久久毛片777| 亚洲精品国产区一区二| 亚洲成av人片免费观看| 国产午夜福利久久久久久| 亚洲色图 男人天堂 中文字幕| 国产99白浆流出| 久久久国产成人免费| 久久久久久人人人人人| 亚洲一卡2卡3卡4卡5卡精品中文| 色哟哟哟哟哟哟| 给我免费播放毛片高清在线观看| 熟女电影av网| 国产色视频综合| 日本免费一区二区三区高清不卡| 亚洲一卡2卡3卡4卡5卡精品中文|