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

    永磁同步電機矢量控制系統(tǒng)頻域模型及其應用

    2022-11-11 03:28:46羅學鴻金佳躍孔建會楊向萍陳家新
    東華大學學報(自然科學版) 2022年5期
    關鍵詞:頻率特性傳遞函數(shù)頻域

    羅學鴻, 金佳躍, 孔建會, 楊向萍, 陳家新

    (東華大學 機械工程學院, 上海 201620)

    隨著永磁同步電機(permanent magnet synchronous motor, PMSM)實際應用要求的不斷提高,對其控制系統(tǒng)的性能提出了更高的要求,相關的分析方法也被不斷完善,涌現(xiàn)出大量的PMSM控制系統(tǒng)建模與參數(shù)整定方法。文獻[1-4]將逆變器等效為一個慣性環(huán)節(jié)進行電流環(huán)分析,通過參數(shù)配置降低了速度環(huán)的階次。王莉娜等[5]基于PMSM調(diào)速系統(tǒng)的頻域模型推導出系統(tǒng)等效傳遞函數(shù),并結合工程應用給出了開環(huán)截止頻率與相位裕度的取值范圍。楊淑英等[6]提出以磁鏈為狀態(tài)量的對稱化建模思路,采用復矢量技術以簡化模型描述。文獻[7-9]基于系統(tǒng)傳遞函數(shù)增益尺度和頻域尺度因子設計了一種帶寬參數(shù)化PD反饋控制器。文獻[10-11]基于二自由度控制方式的研究表明,二自由度控制器使得系統(tǒng)抗擾性與快速性實現(xiàn)獨立調(diào)節(jié)。文獻[12-14]通過研究逆變器與四象限變流器拓撲結構的一致性,進而提出逆變器直流側的等效數(shù)學模型,揭示了直流側等效數(shù)學模型與逆變器及交流側電路參數(shù)之間的定量關系。

    上述研究采用時域狀態(tài)方程模型,而在頻域下伺服系統(tǒng)的特性明顯、參數(shù)聯(lián)系明確,因此有研究者通過頻域模型研究伺服系統(tǒng)參數(shù)與性能的關系。目前對伺服系統(tǒng)的頻域分析與建模以文獻 [15-17]中的模型辨識法為主,該方法建立在試驗的基礎上,即根據(jù)試驗測試結果利用模型辨識得到系統(tǒng)相關參數(shù),進而獲得伺服系統(tǒng)頻域模型。該建模方法較為復雜繁瑣,所得的模型物理意義模糊,為此本文提出一種伺服矢量系統(tǒng)頻域模型的理論建模方法。

    本文提出的雙閉環(huán)矢量控制系統(tǒng)的頻域模型理論建模方法將逆變器延遲與采樣延遲直接體現(xiàn)在系統(tǒng)模型中。通過逆變器的輸出特性以及組成逆變橋的場效應管的物理特性給出了逆變器的等效模型,使用傅里葉級數(shù)展開的方式研究延遲環(huán)節(jié)的替代模型,并通過與模型辨識法所建立的模型對比,檢驗理論模型的正確性,同時在所得等效模型的基礎上利用線性優(yōu)化手段對系統(tǒng)進行改進。

    1 頻域模型關鍵環(huán)節(jié)建模

    雙閉環(huán)矢量控制系統(tǒng)框圖如圖1所示。由圖1可知,要得到系統(tǒng)的頻域模型必須解決以下問題:電機的數(shù)學模型在經(jīng)過坐標變換后d、q軸電流雖解耦合,但轉矩電流方程中反電動勢的存在導致電流依舊存在耦合關系,在建模時應如何處理反電動勢對電流環(huán)模型的影響;三相逆變器拓撲結構的精確數(shù)學模型無法獲得,如何建立逆變器的等效數(shù)學模型使其與實際效果吻合;矢量控制系統(tǒng)多為數(shù)字化控制系統(tǒng),數(shù)字系統(tǒng)信號處理過程會帶來延遲作用,而延遲時間的確定將影響系統(tǒng)模型的準確度,應如何確定延遲時間。針對上述問題對頻域模型的關鍵環(huán)節(jié)進行分析,并給出解決方案以期為頻域模型構建提供依據(jù)。

    圖1 雙閉環(huán)系統(tǒng)控制框圖

    1.1 轉矩電流解耦分析

    永磁同步電機在旋轉坐標下的電壓方程為

    (1)

    式中:ud為d軸電壓;uq為q軸電壓;R為定子繞組電阻;Ld為d軸電感;Lq為q軸電感;ωe為電機的電角度;id、iq分別為d、q軸電流;φf為電機永磁體磁鏈。

    當采用id=0的控制策略時,上述電壓方程簡化為式(2)的形式。

    (2)

    式中:ωeφf為反電動勢。

    將式(2)中的q軸電壓方程中的反電動勢項右移得到式(3)。

    (3)

    采用矢量控制策略將交流電機模擬成直流電機,根據(jù)直流電機的電樞電壓方程,可得q軸電流與反電動勢的關系如下:

    (4)

    式中:iqL為對應的負載電流;Tm為伺服系統(tǒng)的機電時間常數(shù);E為電機旋轉坐標下等效反電動勢。

    對式(4)進行拉普拉斯變換得到反電動勢的傳遞函數(shù)如下:

    (5)

    根據(jù)式(5)得到永磁同步電機電樞反應過程如圖2所示。

    圖2 q軸電樞反應流程

    為便于分析永磁同步電機雙閉環(huán)矢量控制系統(tǒng)的等效模型,假設負載為零。通過移動框圖節(jié)點得到圖2的等效控制框圖,如圖3所示。

    圖3 空載下電流等效控制框圖

    根據(jù)圖3得到電流的開環(huán)傳遞函數(shù):

    (6)

    由式(6)可知電流環(huán)是一個二階系統(tǒng),進一步對電流環(huán)進行簡化使系統(tǒng)降階,得到圖3所示系統(tǒng)的開環(huán)傳遞函數(shù)的頻率特性,如式(7)所示。

    (7)

    由于系統(tǒng)的電流環(huán)響應角頻率一般都很高,當TmTlω2>>1時,進一步簡化式(7)得到式(8)所示結果。

    (8)

    由式(8)可知,坐標變換后結合id=0的控制策略,轉矩電流可以獨立控制,在該情況下可忽略耦合的影響,電流環(huán)相當于一階慣性環(huán)節(jié),因此解耦后轉矩電流與電壓關系為

    (9)

    1.2 逆變器等效模型

    由矢量控制的實現(xiàn)過程可知,矢量算法的最終執(zhí)行是通過三相逆變器的開關組合順序給定電機運行所需的繞組電壓。因此,在系統(tǒng)的建模中逆變器模型的構建對系統(tǒng)是至關重要的。為了提高傳遞函數(shù)模型的準確性,需考慮逆變器的物理過程對系統(tǒng)帶來的影響,綜合考慮可得電流環(huán)的控制框圖如圖4所示。

    圖4 實際電流內(nèi)環(huán)結構

    圖4中iqref為輸入的q軸參考電流,G(s)為逆變器的等效傳遞函數(shù)。為得到圖4的電流環(huán)傳遞函數(shù),必須先得到逆變器的等效傳遞函數(shù)模型。逆變器的物理過程是通過開關管的不同狀態(tài)輸出不同的電壓組合,根據(jù)這一物理過程可將逆變器等效為特定的環(huán)節(jié)。三相電機控制中典型的三相逆變器拓撲結構的等效電路圖如圖5所示。

    圖5 三相逆變器等效電路模型

    圖5中L為電機的相電感,a、b、c分別對應三相電機的各個相線,N表示星型連接電機的中性點。根據(jù)圖5可知相電壓是開關狀態(tài)的函數(shù),則電機相電壓的關系式為

    (10)

    式中:ek為三相等效反電勢;ukN為相電壓;ik為三相電流;在分析每一相線時k分別對應a、b、c。

    式(10)表明電壓與開關的狀態(tài)有關,根據(jù)開關平均函數(shù)可知,開關的導通時間由占空比控制,在一個周期內(nèi)假設橋臂上管占空比為dk,此時逆變器輸出的相電壓和母線電壓的關系如式(11)所示。

    (11)

    根據(jù)式(11)可得逆變器的輸入電壓到輸出電壓的等效線性關系為

    (12)

    在脈寬調(diào)制中由于載波頻率(開關頻率)遠高于輸出電壓基波頻率,根據(jù)式(12)在逆變器的輸出電壓幅值上可將逆變器等效為一個比例環(huán)節(jié),比例系數(shù)為Kpwm,根據(jù)調(diào)制原理并在矢量控制中采用電壓前饋解耦控制,可將Kpwm近似處理為1。

    上述分析表明逆變器相當于一個比例環(huán)節(jié),但在實際的PWM(pulse width modulation)調(diào)制過程中由于開關器件的特性以及死區(qū)時間的影響,逆變器的輸出會存在延遲。在矢量控制系統(tǒng)硬件執(zhí)行過程中上述延遲主要的來源是電流的采樣更新過程,而電流采樣到PWM更新過程的延遲時間與電流采樣更新算法有關。本文采用單周期雙次采樣電流算法分析電流采樣到更新延遲的產(chǎn)生過程,如圖6所示。

    圖6 電流采樣時序

    由圖6可知,在第n個采樣時間點完成采樣后,要在下一個采樣命令即第n+0.5時刻更新PWM波,這是由硬件的特性所決定的,因此逆變器輸出PWM波的延遲時間為半個載波周期,記為Tcar。

    1.3 數(shù)字控制系統(tǒng)信號恢復

    矢量控制系統(tǒng)多為數(shù)字化系統(tǒng),在數(shù)字控制系統(tǒng)中經(jīng)過采樣的離散信號需要恢復成連續(xù)的信號,實際工程中一般用零階保持器來恢復信號。在MATLAB的矢量控制仿真系統(tǒng)中于離散控制器的積分輸出后面設置零階保持器,在系統(tǒng)建模時需要考慮零階保持器信號恢復過程的影響,由自動控制理論的知識可得零階保持器的頻率特性,如式(13)所示。

    (13)

    式中:Ts為采樣周期;ω為系統(tǒng)響應頻率。

    根據(jù)采樣定律將式(13)中的采樣時間轉換成與頻率相關的函數(shù),可得到式(13)的另一種表達方式如式(14)所示。

    (14)

    式中:ωs為采樣頻率。

    (15)

    由式(14)與(15)可知,采樣頻率足夠高時,零階保持器恢復的信號與被采樣的信號較為接近,此時零階保持器可等效為一個延遲時間為半個采樣周期的延遲環(huán)節(jié)。

    則電流環(huán)帶來的總延遲時間為

    Tdelay=Tcar+TOH

    (16)

    上述分析根據(jù)逆變器的輸出幅值特性將逆變器等效為一個比例環(huán)節(jié),又從開關器件的硬件特性和數(shù)字采樣系統(tǒng)的實現(xiàn)方式的角度分析了逆變器具有的延遲效應,綜合可得逆變器的等效數(shù)學模型如式(17)所示。

    G(s)=Kpwme-τs

    (17)

    式中:τ為電流采樣到逆變器輸出時總的延遲時間,并且τ=Tdelay。

    2 雙閉環(huán)系統(tǒng)頻率特性模型

    2.1 按典型I型系統(tǒng)設計電流環(huán)的等效傳遞函數(shù)

    傳統(tǒng)的雙閉環(huán)矢量控制系統(tǒng)等效處理方式是將電流環(huán)中逆變器看成具有一定時間常數(shù)的一階慣性環(huán)節(jié)和在幅值效果上看成是放大倍數(shù)為1的放大環(huán)節(jié)。據(jù)傳統(tǒng)的分析方法得到q軸電流環(huán)的結構框圖,如圖7所示。

    圖7 逆變器等效為慣性環(huán)節(jié)電流環(huán)框圖

    圖7中第一個慣性環(huán)節(jié)是濾波帶來的影響,第二個慣性環(huán)節(jié)是逆變器的延遲效果。為分析方便,忽略反電動勢帶來的擾動影響,將PI控制器寫成零極點的形式,可得圖7的閉環(huán)傳遞函數(shù)為

    (18)

    式中:kCP為電流環(huán)比例增益;kCI為電流環(huán)積分增益。

    按照典型I型系統(tǒng)的設計方法對系統(tǒng)模型進行簡化,對電流環(huán)的閉環(huán)系統(tǒng)進行降階處理,通過設置PI參數(shù)使式(18)降階。PI參數(shù)設計如式(19)所示。

    (19)

    (20)

    在分析速度環(huán)傳遞函數(shù)時,用式(20)所示的電流環(huán)函數(shù)代替電流內(nèi)環(huán),可得傳統(tǒng)分析法的雙閉環(huán)矢量控制系統(tǒng)結構框圖,如圖8所示。

    圖8 傳統(tǒng)分析法的雙閉環(huán)系統(tǒng)結構框圖

    圖8中:TL表示外部負載;K為電機角速度與轉速的換算系數(shù),K=30/π。根據(jù)圖8可得按典型I型系統(tǒng)處理電流環(huán)后的雙閉環(huán)矢量控制系統(tǒng)的閉環(huán)傳遞函數(shù),如式(21)所示。

    (21)

    式中:J為電機轉動慣量;B為阻尼系數(shù);KSP為按典型Ⅱ型系統(tǒng)參數(shù)設計的速度環(huán)比例增益,KSP=0.125;KSI為按典型Ⅱ型系統(tǒng)參數(shù)設計的速度環(huán)積分增益,KSI=25;pn為電機的極對數(shù)。

    對式(21)所示雙閉環(huán)矢量控制系統(tǒng)的等效模型頻率特性與原矢量控制系統(tǒng)頻率特性進行對比,結果如圖9所示。圖9中離散點是矢量控制系統(tǒng)的頻率特性,實線是按典型I型系統(tǒng)設計方法處理電流環(huán)后矢量系統(tǒng)的傳遞函數(shù)模型的頻率特性。

    圖9 按照典型I型系統(tǒng)參數(shù)整定下的頻率特性

    圖9表明,由傳統(tǒng)的按I型系統(tǒng)設計電流環(huán)的方法得到的傳遞函數(shù)與實際系統(tǒng)相差較大,對系統(tǒng)的分析與優(yōu)化起到的作用較小,難以發(fā)揮頻域分析的指導作用。

    2.2 改進后的等效傳遞函數(shù)

    由于傳統(tǒng)模型與實際的頻率特性吻合精度的不足,提出一種改進的頻域模型理論建模方法,該方法在基于內(nèi)模原理的參數(shù)整定方式下對電流環(huán)進行等效。根據(jù)第1節(jié)的分析以及內(nèi)模參數(shù)的設計可得不含延遲環(huán)節(jié)的電流內(nèi)環(huán)的閉環(huán)傳遞函數(shù),如式(22)所示。

    (22)

    式中:α為電流環(huán)的期望帶寬。

    將逆變器以及采樣保持的延遲環(huán)節(jié)考慮到系統(tǒng)模型中,再將延遲環(huán)節(jié)進行泰勒級數(shù)展開,得到式(23)所示的等效表達式。

    (23)

    式(23)的階次根據(jù)系統(tǒng)的階次以及實際情況選定,當延遲時間較大時,為提高近似精度,延遲環(huán)節(jié)近似模型的階次應該選得高些,本文電流環(huán)等效傳遞函數(shù)的分母只有一階,這是因為若分子近似階次高于分母,在工程中不易實現(xiàn),所以選取一階近似模型。

    此時電流環(huán)的傳遞函數(shù)變?yōu)?/p>

    (24)

    根據(jù)對電流環(huán)的解耦分析以及數(shù)字控制系統(tǒng)中的延遲分析,可得旋轉坐標系下雙閉環(huán)控制系統(tǒng)的等效控制模型,如圖10所示。

    圖10 改進的矢量控制系統(tǒng)傳遞函數(shù)模型

    根據(jù)圖10的控制結構示意圖,推導出矢量控制系統(tǒng)在旋轉坐標下的閉環(huán)傳遞函數(shù)如式(25)所示。

    Go(s)=[-cKPτs2+c(KP-KIτ)s+cKI]/

    [Jπbs3+(Bbπ+JπαKv-cKPτ)s2+

    (BπαKv+cKP-cKIτ)s+cKI]

    (25)

    式中:KP為速度環(huán)比例增益;KI為速度環(huán)積分增益;Kv為逆變器等效比例放大系數(shù);Ke為轉矩系數(shù),Ke=1.5pnφf。滿足:

    (26)

    3 理論模型頻率特性驗證

    為驗證本文構建的傳遞函數(shù)與實際系統(tǒng)的吻合度,從頻率響應特性與階躍響應兩方面進行驗證試驗。參數(shù)設計依據(jù)內(nèi)模控制原理進行整定,最終選取的參數(shù)如表1所示。

    表1 電機與控制器的主要參數(shù)

    電流環(huán)總延遲時間為零階保持器所帶來的延遲與電流采樣更新延遲時間之和,PWM的開關頻率為20 kHz,據(jù)此得到仿真系統(tǒng)的延遲時間常數(shù)τ=1.25×10-5。

    依據(jù)表1所示參數(shù)以及式(25)得到選定參數(shù)下的等效傳遞函數(shù)模型如式(27)所示。

    Go(s)=(-157.84s2+12 548 292.73s+6 313 606 405)/(s3+2 486.39s2+12 550 055.10s+6 313 606 405)

    (27)

    對矢量控制系統(tǒng)與基于圖10構建的傳遞函數(shù)系統(tǒng)進行階躍仿真對比,結果如圖11所示。圖11中傳遞函數(shù)系統(tǒng)所示的曲線是通過本文理論建模方法即圖10得到的頻域模型的階躍響應結果,矢量控制系統(tǒng)所示的曲線是MATLAB中搭建的矢量控制仿真系統(tǒng)階躍響應結果。

    圖11 理論模型與矢量控制系統(tǒng)的階躍仿真對比

    圖11表明,改進電流環(huán)處理方式下的等效傳遞函數(shù)控制模型與矢量控制系統(tǒng)仿真模型兩者的階躍響應性能幾乎一致。為進一步驗證等效模型的正確性,對比模型辨識建模方式所得模型、理論建模方法所得模型及原矢量控制系統(tǒng)的頻率特性的吻合情況,結果如圖12所示。圖12中離散點是內(nèi)模參數(shù)下矢量控制系統(tǒng)的頻率特性,虛線是基于模型辨識所得系統(tǒng)的頻率特性,實線是本文所提理論建模法所得的系統(tǒng)頻率特性。

    圖12 頻率特性對比圖

    從圖12可以看出,矢量控制系統(tǒng)、模型辨識所得的頻域模型以及理論建模所得的頻域模型的頻率特性基本吻合。圖11與圖12的結果表明,本文構建的模型與實際系統(tǒng)在階躍性能與頻率特性上都相吻合,因此可采用構建的頻域模型對矢量控制系統(tǒng)進行分析與優(yōu)化。

    4 系統(tǒng)優(yōu)化

    根據(jù)三階系統(tǒng)的階躍響應理論可知,系統(tǒng)的快速性與超調(diào)量很難同時滿足最優(yōu),即系統(tǒng)響應較快時超調(diào)量隨之增加,要求系統(tǒng)無超調(diào)時系統(tǒng)的過渡時間又會增加。當然目前有很多抑制超調(diào)的策略,比如采用積分分離式PI控制器、跟蹤微分器等減小超調(diào)量。但這些方法設計復雜,本文根據(jù)第3節(jié)構建的傳遞函數(shù)模型,利用線性分析方法對系統(tǒng)進行改進。由自動控制理論中典型校正方式的特點可知,超前校正可提高系統(tǒng)的帶寬。超前校正環(huán)節(jié)及其參數(shù)設計如下:

    (28)

    10×lga=g

    (29)

    (30)

    式中:a為校正分子系數(shù);g為根據(jù)頻率特性曲線配置截止頻率時校正環(huán)節(jié)所需提供的增益;ωc為設定的開環(huán)截止頻率。

    由圖12的頻率特性曲線可知,系統(tǒng)幅值衰減到-3 dB時的頻率為808.92 Hz,為驗證傳遞函數(shù)描述下的系統(tǒng)正確性,選擇校正環(huán)節(jié)所提供的23.7 dB增益,此時截止頻率為2 117.83 Hz,計算得到α=234.42,T=4.911×10-6。將式(28)的校正環(huán)節(jié)加至矢量控制系統(tǒng)的速度環(huán)調(diào)節(jié)器前對系統(tǒng)性能進行驗證。根據(jù)行業(yè)標準,輸入信號的幅值為額定轉速的0.01倍即30,輸入頻率為2 117.83 Hz。系統(tǒng)的帶寬測試結果如圖13所示。圖13中原矢量系統(tǒng)是在MATLAB軟件中搭建的圖1所示系統(tǒng),根據(jù)頻域模型優(yōu)化系統(tǒng)是指在矢量系統(tǒng)中根據(jù)構建的頻域模型計算加入校正環(huán)節(jié)后的系統(tǒng)。

    圖13 基于頻域模型的帶寬擴展試驗結果

    由圖13表明,在2 117.83 Hz正弦測試信號下,加入校正環(huán)節(jié)的矢量控制系統(tǒng)帶寬得到了有效擴展,與理論分析一致,而未加入校正環(huán)節(jié)的矢量控制系統(tǒng)在該帶寬下幅值明顯衰減。

    根據(jù)構建的頻域模型對已有的矢量控制系統(tǒng)進行優(yōu)化,為驗證系統(tǒng)的階躍響應性能是否得到改善,對比校正前后系統(tǒng)的階躍響應,結果如圖14所示。

    圖14 校正前后階躍仿真對比

    圖14表明,矢量控制系統(tǒng)經(jīng)過校正后階躍響應性能得到大幅改善,在過渡階段幾乎無超調(diào),且調(diào)整時間變短了。圖13與圖14的結果進一步說明,本文的頻域模型在頻率響應與階躍響應上都與原系統(tǒng)吻合,并且基于頻率特性分析的校正方式在原系統(tǒng)中得到有效應用,使系統(tǒng)的校正變得簡單。

    5 結 語

    針對現(xiàn)有雙閉環(huán)矢量控制系統(tǒng)的頻域分析模型存在過程繁瑣以及典型數(shù)學模型頻率特性與實際系統(tǒng)誤差較大的缺點,提出了PMSM矢量控制系統(tǒng)頻域模型的理論建模方法,并驗證了該方法在帶寬擴展和抑制超調(diào)兩個環(huán)節(jié)的作用效果。結果表明,系統(tǒng)帶寬從原來的808.92 Hz提高至2 117.83 Hz。該模型可彌補現(xiàn)有分析方法中頻域模型準確性差的不足,為高性能伺服系統(tǒng)的開發(fā)奠定基礎,同時為系統(tǒng)性能優(yōu)化提供可靠依據(jù)。

    猜你喜歡
    頻率特性傳遞函數(shù)頻域
    基于LabVIEW的火焰?zhèn)鬟f函數(shù)測量系統(tǒng)
    測控技術(2018年7期)2018-12-09 08:58:46
    頻域稀疏毫米波人體安檢成像處理和快速成像稀疏陣列設計
    雷達學報(2018年3期)2018-07-18 02:41:34
    考慮頻率特性的變頻負荷模型研究
    紅外成像系統(tǒng)的調(diào)制傳遞函數(shù)測試
    基于改進Radon-Wigner變換的目標和拖曳式誘餌頻域分離
    DAM中波發(fā)射機電聲指標頻率特性初探
    一種基于頻域的QPSK窄帶干擾抑制算法
    基于頻域伸縮的改進DFT算法
    電測與儀表(2015年3期)2015-04-09 11:37:24
    基于AD9854的簡易頻率特性測試儀的設計
    基于傳遞函數(shù)自我優(yōu)化的BP網(wǎng)絡算法改進
    久久人人爽av亚洲精品天堂| 一边摸一边抽搐一进一出视频| 看片在线看免费视频| 久久久国产成人精品二区| 国产黄a三级三级三级人| 99热只有精品国产| 国产精品久久久av美女十八| 国产不卡一卡二| 1024香蕉在线观看| 亚洲精品久久成人aⅴ小说| 天天躁狠狠躁夜夜躁狠狠躁| 午夜福利高清视频| 国产麻豆69| 一进一出好大好爽视频| 午夜视频精品福利| 岛国在线观看网站| 日韩大尺度精品在线看网址 | 精品一区二区三区四区五区乱码| 日本a在线网址| 亚洲,欧美精品.| 国产精品影院久久| 国产精品99久久99久久久不卡| 久久久久精品国产欧美久久久| 女警被强在线播放| 久久亚洲真实| 成人免费观看视频高清| 欧美日本中文国产一区发布| 亚洲成人国产一区在线观看| 婷婷六月久久综合丁香| 亚洲精品在线观看二区| 国产精品1区2区在线观看.| 午夜福利影视在线免费观看| 精品久久久久久久人妻蜜臀av | 97人妻天天添夜夜摸| 国产免费男女视频| 国产成人精品久久二区二区免费| 99久久综合精品五月天人人| 别揉我奶头~嗯~啊~动态视频| av视频在线观看入口| 日韩av在线大香蕉| 狠狠狠狠99中文字幕| 90打野战视频偷拍视频| 亚洲精品美女久久久久99蜜臀| 久久亚洲精品不卡| 久久国产精品男人的天堂亚洲| 亚洲精品一卡2卡三卡4卡5卡| 午夜日韩欧美国产| 国产精品一区二区免费欧美| 亚洲成国产人片在线观看| 性欧美人与动物交配| 免费不卡黄色视频| 精品欧美国产一区二区三| cao死你这个sao货| 亚洲精品国产精品久久久不卡| 18禁国产床啪视频网站| 97人妻精品一区二区三区麻豆 | 麻豆久久精品国产亚洲av| 精品卡一卡二卡四卡免费| 国产成人av激情在线播放| 久久久久精品国产欧美久久久| av片东京热男人的天堂| 亚洲一区高清亚洲精品| 两性夫妻黄色片| 老鸭窝网址在线观看| 丝袜美腿诱惑在线| 天天添夜夜摸| 一个人免费在线观看的高清视频| 国产成年人精品一区二区| 亚洲视频免费观看视频| 美女大奶头视频| 亚洲av成人不卡在线观看播放网| 亚洲精品在线美女| 日韩国内少妇激情av| 丝袜美足系列| 亚洲第一电影网av| 国产精品亚洲av一区麻豆| 波多野结衣一区麻豆| 99riav亚洲国产免费| 日韩大码丰满熟妇| bbb黄色大片| 91麻豆av在线| 人人澡人人妻人| 欧美乱色亚洲激情| 国产精品1区2区在线观看.| 亚洲中文字幕日韩| 一边摸一边做爽爽视频免费| 一级黄色大片毛片| 91字幕亚洲| 亚洲片人在线观看| 久久精品aⅴ一区二区三区四区| 国产真人三级小视频在线观看| 女人高潮潮喷娇喘18禁视频| 91字幕亚洲| 一进一出抽搐gif免费好疼| 欧美成人午夜精品| 少妇 在线观看| 久久精品亚洲熟妇少妇任你| 国产成人免费无遮挡视频| 成人精品一区二区免费| 99久久99久久久精品蜜桃| 成人精品一区二区免费| 黄色毛片三级朝国网站| 91精品国产国语对白视频| 亚洲 欧美一区二区三区| 亚洲av片天天在线观看| 看黄色毛片网站| 国产免费男女视频| 搞女人的毛片| 国产精品国产高清国产av| 9191精品国产免费久久| 99国产精品99久久久久| 高潮久久久久久久久久久不卡| 免费在线观看日本一区| 搡老熟女国产l中国老女人| 欧美不卡视频在线免费观看 | 美女免费视频网站| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品影院久久| 91字幕亚洲| 亚洲全国av大片| 精品电影一区二区在线| 夜夜看夜夜爽夜夜摸| 国产精品1区2区在线观看.| 亚洲av五月六月丁香网| 欧美精品啪啪一区二区三区| 成人国产综合亚洲| 九色亚洲精品在线播放| 亚洲色图av天堂| 亚洲 欧美 日韩 在线 免费| 亚洲国产欧美日韩在线播放| 少妇粗大呻吟视频| 老司机靠b影院| 国产成人欧美| 国产在线精品亚洲第一网站| 亚洲av五月六月丁香网| 色尼玛亚洲综合影院| 99久久99久久久精品蜜桃| 亚洲中文字幕日韩| 两个人免费观看高清视频| 少妇被粗大的猛进出69影院| 久久精品aⅴ一区二区三区四区| 在线天堂中文资源库| 亚洲性夜色夜夜综合| 久久人妻福利社区极品人妻图片| 琪琪午夜伦伦电影理论片6080| 欧美性长视频在线观看| 香蕉久久夜色| 美国免费a级毛片| 国产亚洲精品一区二区www| 亚洲国产欧美日韩在线播放| 97人妻精品一区二区三区麻豆 | 久久久久国产一级毛片高清牌| 久久久久国产一级毛片高清牌| 亚洲欧美精品综合久久99| 亚洲欧美精品综合久久99| 免费久久久久久久精品成人欧美视频| 国产精品永久免费网站| 99精品在免费线老司机午夜| 又黄又爽又免费观看的视频| 身体一侧抽搐| 免费少妇av软件| 久久国产精品男人的天堂亚洲| 精品人妻在线不人妻| 在线免费观看的www视频| 亚洲男人天堂网一区| videosex国产| 亚洲狠狠婷婷综合久久图片| 男女下面插进去视频免费观看| 91国产中文字幕| 大香蕉久久成人网| 亚洲人成电影免费在线| 人人妻人人澡人人看| 老司机福利观看| 天天一区二区日本电影三级 | 亚洲精华国产精华精| 一卡2卡三卡四卡精品乱码亚洲| 欧美绝顶高潮抽搐喷水| 啦啦啦免费观看视频1| 人成视频在线观看免费观看| 亚洲国产精品久久男人天堂| 精品卡一卡二卡四卡免费| 在线观看www视频免费| 久久人人精品亚洲av| 免费av毛片视频| 久久精品国产亚洲av高清一级| 国产精品自产拍在线观看55亚洲| cao死你这个sao货| 人人妻人人澡人人看| 91大片在线观看| 777久久人妻少妇嫩草av网站| 久久人人97超碰香蕉20202| 国产一卡二卡三卡精品| 国产精品久久久久久人妻精品电影| 国产极品粉嫩免费观看在线| 一级毛片精品| 一区二区三区高清视频在线| 宅男免费午夜| 黄频高清免费视频| 久久久久久亚洲精品国产蜜桃av| 男男h啪啪无遮挡| 999精品在线视频| 一二三四社区在线视频社区8| tocl精华| 欧美一级a爱片免费观看看 | a级毛片在线看网站| 国产97色在线日韩免费| 亚洲 欧美一区二区三区| 久久久久国产精品人妻aⅴ院| 丝袜在线中文字幕| 亚洲精品久久成人aⅴ小说| 久久精品影院6| 免费看美女性在线毛片视频| 欧美成人性av电影在线观看| 后天国语完整版免费观看| 婷婷丁香在线五月| 成人欧美大片| 熟妇人妻久久中文字幕3abv| 动漫黄色视频在线观看| 国产精品久久久久久亚洲av鲁大| 波多野结衣高清无吗| 91成人精品电影| 免费少妇av软件| 亚洲人成77777在线视频| 男女下面插进去视频免费观看| 日本 欧美在线| 桃色一区二区三区在线观看| 午夜免费鲁丝| 日韩精品中文字幕看吧| 亚洲欧美日韩另类电影网站| 天天躁夜夜躁狠狠躁躁| 精品国产超薄肉色丝袜足j| 乱人伦中国视频| 久久精品国产清高在天天线| 午夜福利免费观看在线| 一个人观看的视频www高清免费观看 | 国产亚洲av嫩草精品影院| 亚洲欧美日韩高清在线视频| 自线自在国产av| 国产精品一区二区三区四区久久 | 国产xxxxx性猛交| 国产精品免费一区二区三区在线| 淫秽高清视频在线观看| 又大又爽又粗| 日韩欧美在线二视频| 欧美日本视频| 亚洲成人免费电影在线观看| 波多野结衣av一区二区av| 天堂动漫精品| 亚洲av电影在线进入| 久久久久久久精品吃奶| 人妻久久中文字幕网| 久久久久国产精品人妻aⅴ院| 色在线成人网| av在线天堂中文字幕| 国产精品美女特级片免费视频播放器 | 国产成年人精品一区二区| 狠狠狠狠99中文字幕| 在线观看日韩欧美| 色综合欧美亚洲国产小说| 国产成人精品久久二区二区免费| 亚洲国产高清在线一区二区三 | 精品国产乱子伦一区二区三区| 亚洲av片天天在线观看| 999久久久精品免费观看国产| 中文字幕人成人乱码亚洲影| 啦啦啦 在线观看视频| 免费搜索国产男女视频| 亚洲人成电影免费在线| 伦理电影免费视频| 精品福利观看| 国产av在哪里看| 日韩三级视频一区二区三区| а√天堂www在线а√下载| 亚洲第一欧美日韩一区二区三区| 在线国产一区二区在线| 精品久久久久久,| 欧美日韩乱码在线| bbb黄色大片| 成人精品一区二区免费| 欧美av亚洲av综合av国产av| 一夜夜www| 久久精品国产99精品国产亚洲性色 | 首页视频小说图片口味搜索| 热re99久久国产66热| 午夜影院日韩av| 中亚洲国语对白在线视频| 淫妇啪啪啪对白视频| 午夜视频精品福利| 亚洲av电影在线进入| 国产精品亚洲av一区麻豆| 欧美国产精品va在线观看不卡| 亚洲国产看品久久| 欧美乱妇无乱码| 亚洲精品av麻豆狂野| 亚洲 欧美 日韩 在线 免费| 亚洲第一青青草原| 嫁个100分男人电影在线观看| 亚洲天堂国产精品一区在线| 免费在线观看影片大全网站| 久久中文字幕人妻熟女| 日韩一卡2卡3卡4卡2021年| www日本在线高清视频| 亚洲男人天堂网一区| 久久午夜亚洲精品久久| 国产精品野战在线观看| 夜夜看夜夜爽夜夜摸| 国产蜜桃级精品一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 禁无遮挡网站| 国产激情久久老熟女| 欧美精品啪啪一区二区三区| 999久久久国产精品视频| 精品国产美女av久久久久小说| 一夜夜www| 精品国产一区二区久久| 黑人欧美特级aaaaaa片| 久久久久国产精品人妻aⅴ院| 亚洲专区中文字幕在线| 国产私拍福利视频在线观看| 视频在线观看一区二区三区| 中文字幕最新亚洲高清| 99久久久亚洲精品蜜臀av| 88av欧美| 亚洲人成电影观看| 久久热在线av| 国产亚洲精品久久久久久毛片| 国产免费av片在线观看野外av| 欧美日韩乱码在线| 亚洲成a人片在线一区二区| 黑人欧美特级aaaaaa片| 曰老女人黄片| 色在线成人网| 亚洲人成77777在线视频| 欧美乱妇无乱码| 丰满人妻熟妇乱又伦精品不卡| 亚洲中文日韩欧美视频| www.自偷自拍.com| 中文字幕人妻丝袜一区二区| 亚洲国产精品成人综合色| 久99久视频精品免费| 两个人看的免费小视频| 欧美精品啪啪一区二区三区| ponron亚洲| 丰满的人妻完整版| 51午夜福利影视在线观看| 后天国语完整版免费观看| √禁漫天堂资源中文www| 9191精品国产免费久久| 妹子高潮喷水视频| 男女下面进入的视频免费午夜 | 日韩大尺度精品在线看网址 | 国产精品日韩av在线免费观看 | 99精品在免费线老司机午夜| 老司机在亚洲福利影院| 午夜福利18| 天天添夜夜摸| 99精品久久久久人妻精品| 国产单亲对白刺激| 怎么达到女性高潮| 成年版毛片免费区| 国产精品久久久久久亚洲av鲁大| 91精品三级在线观看| 人人妻人人澡欧美一区二区 | 中文字幕人妻熟女乱码| 欧美中文综合在线视频| 国产精品亚洲美女久久久| 久久人人97超碰香蕉20202| 88av欧美| 好男人在线观看高清免费视频 | 午夜久久久在线观看| 欧美大码av| 老司机在亚洲福利影院| 国产不卡一卡二| 日本a在线网址| 欧美色欧美亚洲另类二区 | 久久狼人影院| 欧美成人一区二区免费高清观看 | 9色porny在线观看| 国产亚洲av高清不卡| 18禁国产床啪视频网站| 麻豆国产av国片精品| 亚洲 欧美一区二区三区| 国产一区二区三区视频了| 国产一区二区激情短视频| 禁无遮挡网站| 高清在线国产一区| 国语自产精品视频在线第100页| 国产精品久久久av美女十八| 国产精品自产拍在线观看55亚洲| 热99re8久久精品国产| 怎么达到女性高潮| 老司机靠b影院| 久久久水蜜桃国产精品网| 香蕉国产在线看| 91老司机精品| 不卡一级毛片| 中文字幕人妻丝袜一区二区| 久久性视频一级片| xxx96com| 亚洲精品国产区一区二| 亚洲第一青青草原| 69精品国产乱码久久久| 精品久久久久久久久久免费视频| 老熟妇仑乱视频hdxx| 露出奶头的视频| 人妻久久中文字幕网| 国产视频一区二区在线看| 国产亚洲欧美在线一区二区| 亚洲欧美精品综合久久99| 女同久久另类99精品国产91| 1024香蕉在线观看| 丝袜人妻中文字幕| 欧美国产日韩亚洲一区| 午夜两性在线视频| 宅男免费午夜| 国产欧美日韩一区二区精品| 一级毛片高清免费大全| 国产成人系列免费观看| av视频在线观看入口| 精品欧美一区二区三区在线| 又黄又粗又硬又大视频| 亚洲国产看品久久| 1024视频免费在线观看| 级片在线观看| 又紧又爽又黄一区二区| 宅男免费午夜| 麻豆国产av国片精品| 搡老熟女国产l中国老女人| 国产亚洲精品第一综合不卡| 国产亚洲精品久久久久5区| 国产精品永久免费网站| 国产亚洲av高清不卡| 女性生殖器流出的白浆| 一进一出抽搐动态| 国产亚洲精品av在线| www.自偷自拍.com| 日韩欧美免费精品| 成年女人毛片免费观看观看9| 女性被躁到高潮视频| 久久精品aⅴ一区二区三区四区| 亚洲五月天丁香| 啦啦啦观看免费观看视频高清 | 香蕉丝袜av| 午夜成年电影在线免费观看| 黑人巨大精品欧美一区二区蜜桃| av超薄肉色丝袜交足视频| 久久热在线av| 午夜精品国产一区二区电影| 久久香蕉国产精品| 日韩一卡2卡3卡4卡2021年| 又黄又粗又硬又大视频| 午夜福利欧美成人| 大型av网站在线播放| 中文字幕最新亚洲高清| 女人高潮潮喷娇喘18禁视频| 高清黄色对白视频在线免费看| 久久久国产成人精品二区| 757午夜福利合集在线观看| 亚洲,欧美精品.| 在线观看免费午夜福利视频| 99国产精品99久久久久| 精品乱码久久久久久99久播| 久久久水蜜桃国产精品网| 国产亚洲欧美98| 精品久久久精品久久久| 男人舔女人下体高潮全视频| 成人手机av| 久久国产亚洲av麻豆专区| 99re在线观看精品视频| 久9热在线精品视频| 老司机深夜福利视频在线观看| 国产免费av片在线观看野外av| 男人操女人黄网站| 咕卡用的链子| 国产精品久久久久久精品电影 | 亚洲一卡2卡3卡4卡5卡精品中文| 一级毛片精品| 91成人精品电影| 一区在线观看完整版| 无限看片的www在线观看| 99国产精品免费福利视频| 99久久精品国产亚洲精品| 手机成人av网站| 国产一区二区激情短视频| 久久午夜综合久久蜜桃| www.自偷自拍.com| 国产av又大| 亚洲午夜精品一区,二区,三区| 国产精品久久久人人做人人爽| 国产精品美女特级片免费视频播放器 | а√天堂www在线а√下载| 成人18禁高潮啪啪吃奶动态图| 波多野结衣高清无吗| 又大又爽又粗| 男男h啪啪无遮挡| 天天躁狠狠躁夜夜躁狠狠躁| 国产熟女午夜一区二区三区| 免费无遮挡裸体视频| 在线天堂中文资源库| 精品欧美一区二区三区在线| 妹子高潮喷水视频| 免费女性裸体啪啪无遮挡网站| www.熟女人妻精品国产| ponron亚洲| 亚洲七黄色美女视频| 国产高清videossex| 男人舔女人的私密视频| av超薄肉色丝袜交足视频| 一级毛片高清免费大全| 久久精品国产亚洲av香蕉五月| 操出白浆在线播放| 超碰成人久久| 国产精品,欧美在线| 桃色一区二区三区在线观看| 日本黄色视频三级网站网址| av在线天堂中文字幕| 91老司机精品| 日本 av在线| 欧美黄色淫秽网站| 91字幕亚洲| 又大又爽又粗| 亚洲国产日韩欧美精品在线观看 | 真人一进一出gif抽搐免费| 成人永久免费在线观看视频| 精品乱码久久久久久99久播| 欧美日韩亚洲国产一区二区在线观看| 久久人人97超碰香蕉20202| 老鸭窝网址在线观看| 一本大道久久a久久精品| 在线免费观看的www视频| 首页视频小说图片口味搜索| 久久国产精品影院| 黄频高清免费视频| 日日爽夜夜爽网站| 国产av一区二区精品久久| 国产97色在线日韩免费| 在线观看66精品国产| 亚洲久久久国产精品| 国产激情久久老熟女| 午夜免费成人在线视频| 精品欧美国产一区二区三| 国产精品1区2区在线观看.| 18禁美女被吸乳视频| 免费一级毛片在线播放高清视频 | 后天国语完整版免费观看| av电影中文网址| 免费在线观看完整版高清| 久久精品国产99精品国产亚洲性色 | 欧美乱码精品一区二区三区| 长腿黑丝高跟| 亚洲天堂国产精品一区在线| 欧美在线黄色| 在线av久久热| 黑人巨大精品欧美一区二区蜜桃| 免费人成视频x8x8入口观看| 久久久久久久久中文| 亚洲,欧美精品.| 亚洲男人的天堂狠狠| 国产xxxxx性猛交| av片东京热男人的天堂| 日韩中文字幕欧美一区二区| av超薄肉色丝袜交足视频| 亚洲va日本ⅴa欧美va伊人久久| 国产av精品麻豆| 中文字幕高清在线视频| 十八禁网站免费在线| 久久午夜综合久久蜜桃| 欧美激情 高清一区二区三区| 级片在线观看| 男女做爰动态图高潮gif福利片 | 看免费av毛片| 亚洲,欧美精品.| 伦理电影免费视频| 一a级毛片在线观看| 久久草成人影院| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲人成伊人成综合网2020| 欧美一区二区精品小视频在线| 国产av又大| 亚洲中文字幕一区二区三区有码在线看 | 悠悠久久av| 男人操女人黄网站| 神马国产精品三级电影在线观看 | 日韩欧美国产一区二区入口| 久久九九热精品免费| 久久久久久人人人人人| 很黄的视频免费| 欧美av亚洲av综合av国产av| 日韩欧美在线二视频| 精品国产美女av久久久久小说| 精品久久久久久久久久免费视频| 欧美日本中文国产一区发布| 可以在线观看毛片的网站| 精品乱码久久久久久99久播| 91老司机精品| 色综合站精品国产| 精品国产国语对白av| 叶爱在线成人免费视频播放| 国产亚洲精品一区二区www| 亚洲一码二码三码区别大吗| 又黄又爽又免费观看的视频| 在线观看免费日韩欧美大片| 久久人人97超碰香蕉20202| 性欧美人与动物交配| 一进一出好大好爽视频| 9色porny在线观看| 亚洲精品美女久久av网站| 99国产精品99久久久久| 国产精品av久久久久免费| 亚洲国产精品合色在线| 免费看十八禁软件| 亚洲精品国产区一区二| 人人妻人人澡人人看| 国产野战对白在线观看| 校园春色视频在线观看|