999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

船舶電力推進系統(tǒng) VC 與 DTC 調(diào)速控制策略建模與仿真研究

2017-06-19 19:18:59劉恩東羅永吉朱漢華司宇航
艦船科學(xué)技術(shù) 2017年5期
關(guān)鍵詞:控制策略船舶系統(tǒng)

劉恩東,高 嵐,羅永吉,朱漢華,司宇航,胡 佳

(武漢理工大學(xué) 船舶動力工程技術(shù)交通行業(yè)重點實驗室,湖北 武漢 430063)

船舶電力推進系統(tǒng) VC 與 DTC 調(diào)速控制策略建模與仿真研究

劉恩東,高 嵐,羅永吉,朱漢華,司宇航,胡 佳

(武漢理工大學(xué) 船舶動力工程技術(shù)交通行業(yè)重點實驗室,湖北 武漢 430063)

為研究不同永磁同步電機調(diào)速控制策略對船舶大功率推進電機調(diào)速性能的影響和在電力推進系統(tǒng)中的適用性,在推進電機及其螺旋槳負載數(shù)學(xué)模型的基礎(chǔ)上,建立永磁同步電機矢量控制(VC)系統(tǒng)和直接轉(zhuǎn)矩控制(DTC)系統(tǒng)。根據(jù)實船系統(tǒng)建立船舶電力推進系統(tǒng)并進行仿真,分析在船舶加速工況下 2 種控制策略的調(diào)速性能以及船舶電站的穩(wěn)定性。仿真結(jié)果表明;2 種控制策略都有很好的控制效果,在保證良好的調(diào)速性能同時,保持整個系統(tǒng)的穩(wěn)定性,二者在船舶電力推進系統(tǒng)中都具有一定的適用性。

電力推進系統(tǒng);永磁同步電機;矢量控制;直接轉(zhuǎn)矩控制

0 引 言

永磁同步電動機(PMSM)以其效率高、功率因數(shù)高、力矩慣量比大、定子電流和定子電阻損耗小且轉(zhuǎn)子參數(shù)可測、控制性能好等優(yōu)點,成為船舶推進電機的最優(yōu)選擇。作為船舶電力推進的核心,交流電機調(diào)速目前有2種比較成熟且應(yīng)用廣泛的控制策略,分別是矢量控制(Vector Control,VC)和直接轉(zhuǎn)矩控制(Direct Torque Control,DTC)[1]。雖然這 2 種控制方法在其它領(lǐng)域已得到廣泛應(yīng)用,但由于船舶推進電機具有大功率、大轉(zhuǎn)動慣量等特點[2],且運行工況較為復(fù)雜,所以 VC 和 DTC 兩種控制策略能否適用于船舶電力推進系統(tǒng)的永磁同步電動機調(diào)速控制以及兩種策略的調(diào)控效果優(yōu)劣仍有待深入研究。

本文在船舶永磁同步推進電機模型的基礎(chǔ)上分別建立矢量控制和直接轉(zhuǎn)矩控制系統(tǒng),在 Matlab/Sim-ulink 仿真平臺上進行全船電力推進系統(tǒng)建模與仿真分析,對永磁同步電機的矢量控制、直接轉(zhuǎn)矩控制效果進行對比分析。

1 推進電機及負載建模

1.1 永磁同步電機模型

永磁同步電機(PMSM)的磁場是由永磁體產(chǎn)生,既節(jié)約了能耗,又簡化了電機結(jié)構(gòu),提高了電機的質(zhì)能比。與交流電機相比,永磁同步電機沒有激磁損耗,不需要從電網(wǎng)中吸收無功功率,提高了功率因素。PMSM 電機的物理模型如圖 1所示。

為了方便分析數(shù)學(xué)方程,運用坐標變換原理做線性化處理,將圖 1 分解成不同坐標系,實現(xiàn)電機解耦控制,通常所用到的坐標有 ABC 三相靜止坐標系、αβ 兩相靜止坐標系和 dq 兩相旋轉(zhuǎn)坐標系 3 種,如圖 2所示。

通常情況下,在永磁同步電機建立數(shù)學(xué)模型時一般忽略鐵芯飽和,不考慮轉(zhuǎn)子上的阻尼繞組,不計渦流損耗,不考慮諧波且感應(yīng)電動勢為理想正弦波,兩相旋轉(zhuǎn)坐標系中的 PMSM 的電壓方程,即數(shù)學(xué)模型為[3]:

式中:ud,uq分別為定子在 dq 坐標下直軸電壓和交軸電壓;id,iq為定子在 dq 坐標系下的電流分量;ω 為轉(zhuǎn)子磁鏈旋轉(zhuǎn)角速度;ψd,ψq為定子在 dq 坐標系下的磁鏈分量。磁鏈方程為:

式中:Ld,Lq分別為 PMSM 的直軸電感和交軸電感;ψf為永磁體的磁鏈。

電機轉(zhuǎn)矩方程:

電機機械運動方程:

式中:TL為負載轉(zhuǎn)矩;J 為轉(zhuǎn)動慣量;F 為阻尼系數(shù);pn為電機磁極對數(shù);ωe為轉(zhuǎn)子電角度。

1.2 螺旋槳負載模型

螺旋槳與船體之間具有強耦合性,一般將二者作為一個整體建模,數(shù)學(xué)模型如下:

螺旋槳轉(zhuǎn)矩

螺旋槳推力

船槳系統(tǒng)運動方程

船體與螺旋槳間的關(guān)系

有界形式的進速比

式中:vs為船速;vp為螺旋槳前進速度;ω 為伴流系數(shù);ρ 為海水密度;D 為螺旋槳直徑;n 為螺旋槳轉(zhuǎn)速;t 為推力減額系數(shù);R 為船舶阻力;m 為船體重量;Δm 為隨船一起運動的附著水質(zhì)量;K’M和 K’p分別為修正扭矩系數(shù)和推力系數(shù),是 J′ 的函數(shù)。

2 推進電機調(diào)速控制策略

2.1 矢量控制(VC)

推進電機的速度調(diào)節(jié)通常是通過控制電磁轉(zhuǎn)矩來實現(xiàn),對于表貼式三相永磁同步電機,電磁轉(zhuǎn)矩表達式(3)可變?yōu)椋?/p>

式(10)表明,決定電磁轉(zhuǎn)矩的是定子電流在 q軸上的分量,即轉(zhuǎn)矩電流 iq。在 dq 坐標系沿轉(zhuǎn)子磁場定向的前提下,若控制 id= 0,則在 dq 坐標系中定子電流全是轉(zhuǎn)矩電流且 iq與 ψf始終相對靜止。此時從轉(zhuǎn)矩生成的角度,可將 PMSM 等效為他勵直流電動機,控制 iq可獲得與他勵直流電動機同樣的轉(zhuǎn)矩控制效果[4]。本文擬用 PMSM 矢量控制系統(tǒng)(VC),如圖 3 所示。

圖3 中,PMSM 矢量控制系統(tǒng)為雙閉環(huán)控制系統(tǒng),外環(huán)為速度環(huán),內(nèi)環(huán)為電流環(huán)。系統(tǒng)由 1 個轉(zhuǎn)速調(diào)節(jié)器、2 個電流調(diào)節(jié)器、3 個坐標變換模塊、1 個空間電壓矢量脈寬調(diào)制(SVPWM)模塊、角度傳感器、電壓源逆變器和 PMSM[5]組成。系統(tǒng)將檢測到的定子三相電流通過坐標變換轉(zhuǎn)化為 dq 坐標系下的 2 個分量,分別與電流勵磁分量給定值和轉(zhuǎn)矩分量給定值做差,通過電流調(diào)節(jié)器生成參考電壓,經(jīng)坐標變換轉(zhuǎn)化為兩相靜止坐標系下的分量,輸入到 SVPWM 模塊,然后產(chǎn)生 PWM 信號控制逆變器開關(guān),驅(qū)動電機工作。其中定子電流勵磁分量給定值設(shè)為 0,即 id*= 0[6];轉(zhuǎn)矩分量給定值由實際轉(zhuǎn)速與給定轉(zhuǎn)速的差值經(jīng)速度調(diào)節(jié)器得到。該過程通過轉(zhuǎn)子磁場定向和矢量變換,實現(xiàn)了定子電流勵磁分量和轉(zhuǎn)矩分量的解耦,分別控制 2 個分量可獲得與直流電動機同樣的控制品質(zhì)。

2.2 直接轉(zhuǎn)矩控制(DTC)

電磁轉(zhuǎn)矩的生成可看成是定子磁場和轉(zhuǎn)子磁場相互作用的結(jié)果,電磁轉(zhuǎn)矩表達式(3)可變形為

式中: ψs為定子磁鏈;δ 為定子磁鏈和轉(zhuǎn)子磁鏈之間的夾角。轉(zhuǎn)子磁鏈 ψf為常數(shù),若能控制定子磁鏈的幅值不變,則電磁轉(zhuǎn)矩就只和負載角 δ 有關(guān),通過控制δ 可控制電磁轉(zhuǎn)矩的大小,這是直接轉(zhuǎn)矩控制的基本原理[7]。

PMSM 定子電壓方程為:

船用永磁同步電機的定子電阻很小,若忽略定子電阻的影響,則有

所以可通過外加電壓 us來控制定子磁鏈的幅值和轉(zhuǎn)速。在很短時間內(nèi)依靠 us的作用使定子磁鏈加速旋轉(zhuǎn),而這期間轉(zhuǎn)子速度來不及變化,由此可加大負載角 δ,若保持定子磁鏈的幅值不變,就可以使電磁轉(zhuǎn)矩增大;若使定子磁鏈反向旋轉(zhuǎn),可使電磁轉(zhuǎn)矩減小[8]。

本文擬用 PMSM 直接轉(zhuǎn)矩控制系統(tǒng)(DTC)如圖 4 所示,系統(tǒng)由速度調(diào)節(jié)器、滯環(huán)比較器、電壓矢量開關(guān)表、Clark 變換、磁鏈和轉(zhuǎn)矩觀測器,角度傳感器、電壓源逆變器和 PMSM組成。系統(tǒng)將檢測的定子三相電流 iABC經(jīng)過 Clark 變換后輸入磁鏈和轉(zhuǎn)矩觀測器,估算出定子磁鏈和轉(zhuǎn)矩的大小分別與磁鏈和轉(zhuǎn)矩給定值做比較,其差值分別經(jīng)過滯環(huán)比較器得到磁鏈控制信號 Dψ和轉(zhuǎn)矩控制信號 DT,然后根據(jù)定子磁鏈位置 θ(N)、轉(zhuǎn)矩控制信號 DT和磁鏈控制信號 Dψ,利用電壓矢量開關(guān)表選擇響應(yīng)的電壓矢量控制逆變器開關(guān)的通斷,驅(qū)動電機運行。其中磁鏈給定信號為常數(shù),轉(zhuǎn)矩給定信號由實際轉(zhuǎn)速與給定轉(zhuǎn)速的差值經(jīng)速度調(diào)節(jié)器得到的。該過程將轉(zhuǎn)矩和磁鏈作為反饋信號,省去復(fù)雜的矢量變換,直接控制電磁轉(zhuǎn)矩,從而達到調(diào)速的目的[9]。

3 船舶電力推進系統(tǒng)仿真分析

除了永磁同步推進電機及螺旋槳負載建模和 2 種

控制系統(tǒng)模型外,船舶電力推進系統(tǒng)中還有柴油發(fā)電機組、高低壓負載等模塊建模,參見文獻[10]。

3.1 仿真對象及系統(tǒng)參數(shù)

以某船為母型船,基于 VC 和 DTC 兩種控制策略,在 Matlab/Simulink 仿真平臺中建立船舶電力推進系統(tǒng)模型,如圖 5 和圖 6 所示。該船采用的是 6 600 V、50 Hz 船舶中壓電力系統(tǒng),由 2 個中心電站組成,每個電站有 2 臺中壓柴油發(fā)電機組;2 套推進系統(tǒng),本文采用 1 個電站和 1 套推進器進行仿真分析。發(fā)電機組額定容量 3 388 kVA,額定電壓 6 600 V,額定電流 296 A,功率因素 0.85,頻率 50 Hz;推進電機的額定功率 4.088 MW,電機額定電壓 660 V,額定轉(zhuǎn)速 200 r/min,額定轉(zhuǎn)矩 Tn = 195 200 N·m;螺旋槳直徑 3.6 m,槳葉數(shù)為 4,船體質(zhì)量 16 229 t。

3.2 仿真設(shè)置

仿真時間長度為 8 s,0 s 時發(fā)電機組啟動,達到額定電壓后接入電網(wǎng);1 s 時船上其它負載投入運行;3 s 時船舶啟動航行,初始轉(zhuǎn)速設(shè)定為 100 r/min;4 s 時加速,設(shè)定速度變?yōu)?150 r/min;5 s 時加入負載轉(zhuǎn)矩 100 kN·m擾動,6 s 時擾動消失。電力系統(tǒng)中觀測的主要參數(shù)有發(fā)電機有功功率、發(fā)電機端電壓、勵磁電壓、頻率。推進系統(tǒng)中觀測的參數(shù)有推進電機轉(zhuǎn)速、電磁轉(zhuǎn)矩、定子三相電流、船速等。

3.3 仿真結(jié)果分析

3.3.1 船舶電站特性參數(shù)分析

船舶電站主要特性參數(shù)的系統(tǒng)仿真結(jié)果如下:

2 種控制策略下船舶電站及負載的變化趨勢基本相同,船舶電站的參數(shù)變化曲線如圖 7 所示。取其中一種情況加以說明,1 s 時輔機及日常用電設(shè)備投入運行,發(fā)電機有功功率增大,經(jīng)過短暫的過渡時間后進入穩(wěn)定狀態(tài),勵磁電壓立刻上升然后回落至額定值,使得發(fā)電機端電壓基本無變化,向電網(wǎng)穩(wěn)定供電,同時發(fā)電機頻率經(jīng)過小范圍下降后回到額定值;3 s 時大功率推進系統(tǒng)投入運行,發(fā)電機功率增大,勵磁電流瞬時大范圍增大以調(diào)節(jié)并補償大功率負載造成發(fā)電機輸出電壓的壓降,發(fā)電機輸出電壓下降 0.5%,對于電力系統(tǒng)影響可以忽略不計;4 s 加速時,發(fā)電機組的勵磁電壓、輸出電壓和頻率與上一秒變化相同,但變化范圍減小;5 s 加入負載擾動時,發(fā)電機組有功功率繼續(xù)增大,推進系統(tǒng)的轉(zhuǎn)矩呈現(xiàn)波動狀態(tài),此時勵磁電壓增大,發(fā)電機輸出電壓有少量的波動,電壓波動不到 0.7%,頻率降低;6 s 負載擾動消失,發(fā)電機組有功功率減小,勵磁電壓降低,發(fā)電機輸出電壓波動逐漸消失并穩(wěn)定在額定值,頻率增大然后穩(wěn)定在額定值。整個過程中,頻率相對保持穩(wěn)定,最大變化率不超過0.2%,在允許范圍內(nèi)。由此可見,基于 2 種控制策略的推進系統(tǒng)在運行時對于船舶電力系統(tǒng)的影響不大,都在可調(diào)節(jié)的范圍內(nèi),故 2 種控制策略在船舶電力推進系統(tǒng)都有一定的適用性。

3.3.2 控制策略比較分析

推進系統(tǒng)在最后 1 s 時間內(nèi)穩(wěn)定運行,性能指標沒有變化,所以這里取 2~7 s 的時間段對推進系統(tǒng)性能指標進行觀測,VC 和 DTC 控制策略仿真結(jié)果分別如圖 8 和圖 9 所示。

圖8 和圖 9 分別為矢量控制(VC)和直接轉(zhuǎn)矩控制(DTC)下的推進系統(tǒng)參數(shù)變化曲線,仿真結(jié)果顯示,2種控制算法都有良好的動態(tài)響應(yīng),二者具體性能比較見表 1。在啟動和加速時,2 種控制策略都有快速響應(yīng),其中 DTC 響應(yīng)更快。因為 VC 的輸出電壓由電流調(diào)節(jié)器產(chǎn)生,存在時間上的滯后,而 DTC 通過滯環(huán)比較器和電壓開關(guān)選擇表輸出電壓,沒有電流控制回路,可獲得較大的電壓變化率和較快的電流響應(yīng),所以 DTC 的轉(zhuǎn)矩響應(yīng)比 VC 快些。VC 系統(tǒng)的起動轉(zhuǎn)矩比 DTC 系統(tǒng)大,且啟動電流相對較小。VC 系統(tǒng)的速度環(huán)輸出采取了電流限制,而 DTC 系統(tǒng)為防止過電流,速度環(huán)輸出采取轉(zhuǎn)矩限制措施,因此 DTC 的轉(zhuǎn)矩被限制在一定范圍,而起動電流比 VC 大。

系統(tǒng)加入負載擾動時,VC 系統(tǒng)的過渡時間稍長,轉(zhuǎn)速的動態(tài)偏差也較大,但系統(tǒng)穩(wěn)定后無靜態(tài)偏差;而 DTC 系統(tǒng)轉(zhuǎn)速動態(tài)偏差雖小,但系統(tǒng)穩(wěn)定后存在靜態(tài)偏差。而且在穩(wěn)態(tài)時 VC 系統(tǒng)的轉(zhuǎn)矩和定子電流變化相對平滑;DTC 系統(tǒng)的轉(zhuǎn)矩波動范圍相對較大,較為粗糙,定子電流脈動大,不是標準的正弦電流。在VC 系統(tǒng)中,對定子電流的勵磁和轉(zhuǎn)矩分量進行解耦,并通過 PI 調(diào)節(jié)器分別對電流勵磁分量和轉(zhuǎn)矩分量進行調(diào)節(jié),進而對轉(zhuǎn)矩進行控制。該過程通過 PI 調(diào)節(jié)器實現(xiàn)了電流和轉(zhuǎn)矩的連續(xù)控制,使得電流和轉(zhuǎn)矩波形較為平滑,在穩(wěn)態(tài)時由于轉(zhuǎn)矩誤差小并與負載轉(zhuǎn)矩達到高度平衡,使得轉(zhuǎn)速幾乎沒有靜態(tài)偏差。而 DTC 系統(tǒng)采用滯環(huán)比較器,使得轉(zhuǎn)矩和定子磁鏈總存在誤差,穩(wěn)態(tài)時轉(zhuǎn)矩和電流脈動大;由于轉(zhuǎn)矩波動大,與負載轉(zhuǎn)矩之間總存在差值,使得轉(zhuǎn)速難以精確達到期望值,存在靜態(tài)偏差。

表1 矢量控制與直接轉(zhuǎn)矩控制系統(tǒng)性能比較Tab. 1 Vector control and direct torque control system performance comparison

在整個過程中,2 種控制策略下的船舶速度都較為平穩(wěn),加速時船速響應(yīng)快,負載擾動下,船速幾乎沒有變化,可見兩種控制策略下的船舶電力推進系統(tǒng)靈活性好,操縱性強,而且具有很好的魯棒性。

4 結(jié) 語

本文將永磁同步電機的矢量控制(VC)和直接轉(zhuǎn)矩控制(DTC)策略分別應(yīng)用于船舶電力推進系統(tǒng),仿真結(jié)果表明:

1)2 種控制策略都有良好的動態(tài)響應(yīng),在工況變化時可以保證推進電機調(diào)速的準確性、快速性和穩(wěn)定性。

2)矢量控制系統(tǒng)的算法相對復(fù)雜,但可以實現(xiàn)轉(zhuǎn)矩的平滑控制,轉(zhuǎn)矩和電流波動小;直接轉(zhuǎn)矩控制算法相對簡單,動態(tài)響應(yīng)快,但轉(zhuǎn)矩和電流波動較大。

3)2 種控制策略在船舶電力推進系統(tǒng)中都有一定的適用性,目前矢量控制策略已相對完善,直接轉(zhuǎn)矩控制策略在減小轉(zhuǎn)矩波動等方面仍有改進空間。

[1]池飛飛. 船用永磁同步電機直接轉(zhuǎn)矩控制仿真研究[D]. 武漢: 武漢理工大學(xué), 2013.

[2]李崇堅. 交流同步電機調(diào)速系統(tǒng)[M]. 北京: 科學(xué)出版社, 2013: 292–294.

[3]KRISHNAN R. 永磁無刷電機及其驅(qū)動技術(shù)[M]. 柴鳳, 譯.北京: 機械工業(yè)出版社, 2013: 182–187.

[4]王成元, 夏加寬, 孫宜標. 現(xiàn)代電機控制技術(shù)[M]. 北京: 機械工業(yè)出版社, 2014.

[5]柯常國, 還芳, 王勁. 電力推進系統(tǒng)控制策略仿真[J]. 船電技術(shù), 2016, 36(4): 15–18.

[6]JIA G X, ZHAO H C, SHAO H J. Simulation research on PMSM vector control system based on SVPWM[C]//Electrical and Control Engineering (ICECE), IEEE Conference Publications, 2010: 1936–1940.

[7]楊建飛, 胡育文. 永磁同步電機直接轉(zhuǎn)矩控制轉(zhuǎn)矩調(diào)節(jié)器設(shè)計[J]. 中國電機工程學(xué)報, 2011, 31(9): 76–81.

[8]徐艷平, 鐘彥儒. 扇區(qū)細分和占空比控制相結(jié)合的永磁同步電機直接轉(zhuǎn)矩控制[J]. 中國電機工程學(xué)報, 2009, 29(3): 102–108.

[9]郭昊昊, 劉彥呈, 任俊杰. 電力推進船舶中永磁電機直接轉(zhuǎn)矩控制[J]. 電機與控制應(yīng)用, 2011, 38(11): 31–33.

[10]胡淼. 船舶電力推進永磁同步電機控制系統(tǒng)的研究[D]. 大連: 大連海事大學(xué), 2013.

Research on modeling & simulation analysis of marine electric propulsion system of VC/DTC speed control strategy

LIU En-dong, GAO Lan, LUO Yong-ji, ZHU Han-hua, SI Yu-hang, HU Jia
(Key Laboratory of Marine Power Engineering and Technology, Ministry of Transportation, Wuhan University of Technology, Wuhan 430063, China)

To study of different permanent magnet synchronous motor speed control strategy of power Marine propulsion motor speed control performance and applicability in the electric propulsion system, on the basis of mathematical model of the propulsion motor and propeller load, established the permanent magnet synchronous motor vector control (VC) and direct torque control (DTC) system. Establishment of marine electric propulsion system and simulation based on the real ship, while analyzing two kinds of control strategy of speed control performance and stability of ship power station, under the working condition of ship speed.The simulation results show that the two kinds of control strategy has good control effect, while ensuring good speed regulating performance at the same time, maintain the stability of the whole system, both in the ship electric propulsion system has certain applicability.

electric propulsion system;permanent magnet synchronous motor;vector control;direct torque control

U664

A

1672–7619(2017)05–0088–06

10.3404/j.issn.1672–7619.2017.05.017

2016–08–22;

2016–10–20

國家自然科學(xué)基金資助項目(51179144)

劉恩東(1990–),男,碩士研究生,主要研究方向為輪機仿真及自動化。

猜你喜歡
控制策略船舶系統(tǒng)
計算流體力學(xué)在船舶操縱運動仿真中的應(yīng)用
Smartflower POP 一體式光伏系統(tǒng)
《船舶》2022 年度征訂啟事
船舶(2021年4期)2021-09-07 17:32:22
WJ-700無人機系統(tǒng)
考慮虛擬慣性的VSC-MTDC改進下垂控制策略
能源工程(2020年6期)2021-01-26 00:55:22
ZC系列無人機遙感系統(tǒng)
北京測繪(2020年12期)2020-12-29 01:33:58
船舶!請加速
工程造價控制策略
山東冶金(2019年3期)2019-07-10 00:54:04
現(xiàn)代企業(yè)會計的內(nèi)部控制策略探討
連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
主站蜘蛛池模板: 日韩欧美亚洲国产成人综合| 久视频免费精品6| 欧美一级特黄aaaaaa在线看片| 久久国语对白| 欧美一级高清片久久99| 2022精品国偷自产免费观看| 丰满人妻被猛烈进入无码| 国产日韩精品一区在线不卡| 免费无码AV片在线观看国产| 中文字幕亚洲乱码熟女1区2区| 久久久国产精品无码专区| 日韩欧美国产另类| 日韩人妻少妇一区二区| 成人综合网址| 亚洲欧美在线综合一区二区三区| 欧美在线导航| 在线亚洲天堂| 69视频国产| 国产精品私拍99pans大尺度| 国产成人精品一区二区不卡| 欧美日韩一区二区三区四区在线观看| 国产激爽爽爽大片在线观看| 亚洲免费三区| 99久视频| 狠狠色婷婷丁香综合久久韩国| 亚洲午夜福利在线| 欧美精品v日韩精品v国产精品| 色老头综合网| 国产精品无码AV中文| 成人av专区精品无码国产| 亚洲欧美一级一级a| 亚洲视屏在线观看| 老司机aⅴ在线精品导航| 国产一级妓女av网站| 第一页亚洲| 欧美另类一区| 情侣午夜国产在线一区无码| 91色爱欧美精品www| 精品一區二區久久久久久久網站 | 亚洲国产清纯| 无码中文字幕加勒比高清| 亚洲综合日韩精品| 亚洲AV成人一区国产精品| 91麻豆精品视频| 久久久久久尹人网香蕉| 中文字幕在线观| 一级香蕉人体视频| 秋霞一区二区三区| 日韩色图在线观看| 精品一区二区无码av| 午夜视频免费试看| 国产夜色视频| 一级毛片在线播放免费观看| 美女亚洲一区| 在线精品亚洲一区二区古装| 黄色在线网| 精品国产免费第一区二区三区日韩| 国产成人成人一区二区| 波多野结衣一二三| 亚洲精品视频网| 中国国产A一级毛片| 国产乱人视频免费观看| 日本一本正道综合久久dvd | 国产精品免费电影| 色天堂无毒不卡| 666精品国产精品亚洲| 久久久成年黄色视频| 毛片视频网| 日韩精品亚洲人旧成在线| 亚洲男女在线| 最新加勒比隔壁人妻| 毛片免费视频| 麻豆精品久久久久久久99蜜桃| 啦啦啦网站在线观看a毛片 | 中文字幕无码中文字幕有码在线 | 欧美日韩成人在线观看| 91午夜福利在线观看| 国产精品精品视频| 99精品在线看| 一本大道香蕉高清久久| 乱人伦99久久| av午夜福利一片免费看|