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

次級斷續時直線感應牽引電機的等效電路

2021-03-25 00:12:46羅志昆
電工技術學報 2021年6期
關鍵詞:區域

呂 剛 羅志昆

次級斷續時直線感應牽引電機的等效電路

呂 剛 羅志昆

(北京交通大學電氣工程學院 北京 100044)

直線感應電機作為軌道交通車輛的牽引電機,在道岔、轉彎等路段會出現次級感應板缺失的情況,次級連續時的等效電路不再適用。針對此特殊工況,提出一種基于氣隙磁通密度的分布模型,考慮次級斷續的等效電路。首先,建立直線感應電機次級斷續時的分析模型,根據初、次級耦合情況把斷續工況分成三個階段。然后,利用線性近似化法修正隨耦合區域長度變化的初級漏感、勵磁電感等參數。結合氣隙磁通密度分布模型與次級斷續工況,考慮縱向端部效應的影響,修正勵磁支路上的勵磁電感和表征端部效應渦流損耗的等效電阻,推導出三個階段的等效電路來表示初級通過次級感應板斷續區域這一動態過程。最后,根據等效電路求出電機推力和效率的動態變化曲線,并與實驗臺的測量結果對比,驗證了所提等效電路的正確性。

直線感應電機 次級斷續 等效電路 縱向端部效應 氣隙磁通密度

0 引言

單邊短初級直線感應電機在交通系統中的運用主要有直線輪軌交通和磁懸浮交通兩種方式。直線輪軌交通的優勢是爬坡、轉彎路段通過性強,線路規劃靈活,拆遷成本低,隧道截面小,建設成本低;磁懸浮交通的優勢則為靜音、舒適、速度快[1-3]。目前由中車集團研制的時速600km磁懸浮試驗車于同濟大學試驗線上成功試跑。直線電機牽引系統克服了傳統軌道交通方式依靠輪軌之間黏著牽引的限制,有望極大地提高陸地交通工具的速度上限。隨著人類對于交通工具速度要求的不斷提高,以及電力電子技術和計算機控制技術的進步,直線電機將會在軌道交通中得到越來越廣泛的應用[4-6]。

在軌道交通系統中,直線感應牽引電機的初級安裝在列車上,次級板沿軌道鋪設。當初級繞組中通入對稱的三相交流電時,列車在電磁推力的牽引下做直線運動。為了保證運行過程中初、次級不發生碰撞,直線電機的氣隙一般較大,從而導致功率因數和能量傳遞效率較低;同時,由于初級開斷造成的端部效應和法向力的存在,使直線感應電機的運行性能進一步惡化。

為提高電機運行效率、減少能量損失,國內外學者提出了許多高性能的優化控制策略[7-8]。一種思路是根據電機的等效電路得到損耗的數學模型,推導出損耗最小時勵磁磁通的表達式。文獻[9]綜合考慮了直線地鐵車輛運行時電機和逆變器的損耗,推導出總損耗最小時的勵磁磁通。文獻[10]根據電機數學模型推導出推力、銅耗和法向力的表達式,結合最優化理論考慮在推力滿足要求的條件下,電機銅耗和法向力造成的損耗之和最小。以上兩種控制策略均需要準確的數學模型,一個精確的等效電路至關重要,因此直線感應電機的等效電路一直是研究熱點。

上述對于直線感應電機等效電路的研究已經相當成熟,但研究對象都是次級連續工況下的電機。由于工程實際的限制,直線感應電機的次級感應板在道岔、轉彎和列車進出庫等路段不能連續鋪設,次級感應板不連續路段如圖1所示,此時會造成電機初級繞組中頻繁過電流,導致保護系統誤動作和牽引傳動系統的IGBT器件損壞。同時推力和法向力的突變導致列車輪對的不規則磨損,嚴重影響列車的安全穩定運行。目前對次級斷續時的相關研究較少,文獻[17]提出一種列車運行過程中直線感應電機互感在線辨識的算法,判斷列車是否進入次級斷續區域,并調整電機dq軸電流的指令值,來抑制初級電流和推力的突變,但使用的仍是基于次級連續工況的等效電路。初級通過次級斷續區域時,初、次級的電阻和漏感等參數會發生變化,導致次級時間常數發生變化,因此次級斷續時傳統針對次級連續的穩態電路不再適用。

圖1 次級感應板不連續路段

針對上述問題,本文提出了一種基于氣隙磁通密度模型的次級斷續時直線感應電機的等效電路。首先建立了次級斷續時直線感應牽引電機的分析模型,將初級通過次級斷續區域的過程分成三個階段。其次,使用線性近似化法修正次級斷續時初級漏感、次級電阻、次級漏感和勵磁電感等參數。然后根據次級斷續時氣隙磁通密度和去磁電流的分布,修正了勵磁支路的電感和電阻。最后在此基礎上,用三個階段的等效電路表示初級通過斷續區域的動態過程,并得到推力和效率的變化曲線,通過與實驗臺測量結果對比,驗證了等效電路的正確性。

1 次級斷續時直線感應電機的分析模型

次級斷續時直線感應牽引電機的分析模型如圖2所示。圖2a是直線感應電機運行的起始位置,初級長度為p,次級斷續區域的長度為s,起始時刻初級入口端與斷續區域距離為c,笛卡爾坐標系的軸正方向與初級運動方向一致,軸正方向垂直于次級板向上,原點位于起始時刻初級的入口端,且坐標系固定不動。初級先后經過次級板1、次級缺失區域和次級板2,如圖2b所示。根據軌道交通中次級板實際的鋪設情況,設置初級長度大于缺失區域的長度。

圖2 次級斷續時直線感應電機的分析模型

為了簡化等效電路參數的求解過程,合理地做出以下假設:

(1)斷續區域的氣隙磁通密度不影響初、次級耦合區域的磁通密度分布。

(2)初級通過次級斷續區域時,速度保持恒定,因為直線感應電機牽引的列車慣性很大。

(3)忽略鐵磁材料的飽和。

(4)初級繞組的三相電流保持對稱。

次級的缺失會對氣隙磁場和次級電流的分布產生很大的影響,從而導致推力、次級渦流損耗、氣隙傳遞的無功功率等物理量發生變化,對電機的運行情況產生較大的影響。綜合考慮耦合區域的長度變化和縱向端部效應的影響,把初級通過斷續區域的過程分成三個階段:階段1,初級僅與次級板1耦合;階段2,初級與次級板1、2同時耦合;階段3,初級僅與次級板2耦合。

2 次級斷續時的等效電路

2.1 考慮轉差率的新型氣隙磁通密度分布模型

在直線感應電機中,初級長度有限并不斷與新的次級區域耦合,氣隙磁場會產生畸變。根據楞次定律,初級入口端的感應渦流會削弱氣隙磁場,出口端的氣隙磁場則會被加強,初級電流和次級感應渦流共同決定了氣隙磁場的分布,這種由于初級不斷進入新的次級區域導致的氣隙磁場畸變稱為動態縱向端部效應,其示意圖如圖3所示。氣隙磁場和次級去磁渦流作用產生的電磁力始終阻止初級進入新的次級區域,縱向端部效應會對電機的運行狀況產生較大影響,而且速度越大,影響越明顯。本文根據縱向氣隙磁通密度和次級去磁渦流的分布,在傳統的T型等效電路上根據運行速度,動態修正并聯支路上的電感和電阻,來考慮縱向端部效應的影響。

圖3 縱向端部效應示意圖

J. Duncan提出用e指數函數擬合縱向氣隙磁通密度的分布(簡稱Duncan模型)[15],如圖4a所示,21、22、23為不同的次級運動速度。為了表示不同速度下縱向氣隙磁通密度和次級去磁渦流的分布情況,/2,2為次級時間常數,2=(m2)/2,m、2和2分別為勵磁電感、次級漏感和次級電阻。當初級運動速度為時,初級在次級上方前進一個自身長度所需時間T=/,定義一個無量綱的參數為

式中,為時間域上電機長度的標幺值,速度越大,電機的標幺長度則越小;反之,則越大。當速度為零時,電機長度趨近于無窮大,縱向端部效應不發生。

圖4b為縱向氣隙磁通密度分布的有限元計算結果,Duncan模型僅考慮了速度對磁通密度分布的影響,然而不同轉差率下磁通密度分布情況相差較大,Duncan模型僅在轉差率較小時比較準確,無法適用于直線感應電機運行的全部工況。為了考慮轉差率的影響并提高等效電路模型的計算精度,使用一種與轉差率有關的函數來模擬不同工況下的縱向氣隙磁通密度分布(簡稱新型磁通密度分布模 型)[16],定義為

圖4 縱向氣隙磁通密度的分布情況

初級供電頻率一定時,不同轉差率下縱向氣隙磁通密度和次級去磁渦流的分布模型如圖5所示。氣隙磁場在初級入口端迅速增強,到達高位后振蕩趨向于穩定,相比于Duncan模型,新型磁通密度分布模型更接近于有限元計算結果,說明新型磁通密度模型更接近于實際情況,能反映不同頻率、不同轉差率下各種運行工況的磁場分布。

圖5 氣隙磁通密度、去磁渦流分布模型

2.2 次級斷續時的等效電路

2.2.1 初級漏感等參數的修正

為了降低分析的復雜程度,使用線性近似化法來修正與耦合區域面積有關的等效電路參數,表達式為

圖6 直線感應電機的等效電路

式中,p、r分別為初級長度和耦合區域長度。上標1為次級完整時的參數值,可通過靜止狀態下的參數辨識求得[18-20];上標0為次級全部缺失時的值,為簡化分析,認為無次級板時初、次級無能量傳遞,故勵磁電感、次級漏感和次級電阻均為零,即

2.2.2 勵磁支路的修正系數

斷續工況階段1初級僅與次級板1耦合,初級入口端進入缺失區域向次級板2靠近。次級去磁渦流的分布情況如圖7所示,1為次級板1的邊緣與初級的交點。初級漏感修正后,認為初、次級非耦合區域氣隙磁通密度為零,次級去磁渦流的分布與氣隙磁場相對應,且缺失區域不產生次級渦流和 推力。

圖7 階段1去磁渦流的分布

縱向端部效應產生的去磁渦流會削弱氣隙磁場,此影響通過修正勵磁電感來考慮,按照圖6a所示的演化過程,有

式中,2ea為去磁渦流的平均值;mea為勵磁電流的平均值;并聯電感me起分流作用,故勵磁電感的修正系數為Km=1-1。

勵磁支路的電阻表征感應板渦流損耗e,由兩部分組成:①初、次級耦合區域的鋁板上的渦流損耗e1;②初級出口端外側感應渦流產生的損耗e2。去磁渦流的方均根值為

故勵磁支路電阻的修正系數為

隨著初級往前運動,次級缺失區域在初級下方逐漸向出口端移動,從階段1過渡到階段2、3。去磁渦流分布曲線的中斷部分隨之往后移動,階段2、3去磁渦流的分布如圖8所示,21、22、3分別為運行階段2、3次級板邊緣與初級的交點。

圖8 階段2、3去磁渦流的分布

根據階段1修正系數的推導過程,結合階段2、3次級去磁渦流的分布,可得到勵磁支路電感和電阻的修正系數Km和K2。至此,次級斷續工況各運行階段勵磁支路的修正系數已全部推導出,見表1,結合圖6的T型電路,次級斷續時的等效電路構建完成。

根據構建的等效電路,可求出次級斷續時推力、效率和功率因數等物理量。推力F由兩部分組成:①次級轉差電流2產生的正推力F1,對應于等效電路中次級支路電阻的損耗;②次級去磁渦流產生的、阻礙初級前進的負推力F2,對應于勵磁支路電阻的損耗。

表1 各運行階段勵磁支路的修正系數

Tab.1 Correction coefficient of excitation branch in each operation stage

效率為

3 等效電路的計算結果與實驗驗證

為了驗證所提出的等效電路的正確性,按照表2所示的電機參數,使用Matlab軟件計算初級通過斷續區域時,等效電路的勵磁支路參數和推力、效率的動態變化曲線。

表2 直線感應電機參數

Tab.2 Parameters of linear induction motor

圖10為直線感應牽引電機實驗臺,電機參數同表2。接觸軌提供1 500V直流電壓,經由VVVF逆變器向電機初級供電,能夠改變電壓(0~DC 1 100V)和頻率(0~35Hz)來驅動。轉向架上安裝了用于測量力、電壓、電流和速度的傳感器,數據采集器用Labview編程,獲取并處理被測數據。分別使用基于新型磁通密度模型和Duncan模型推導的等效電路(下文分別稱為方法1和方法2),計算得到推力和效率的動態變化曲線,與實驗測量結果的對比如圖11所示。

從圖11中可以看出,初級經過斷續區域時,推力和效率均先減小后增大,大體上與耦合區域面積呈正相關,為V型曲線。方法1考慮了整個轉差范圍,得到了更精確的結果,推力、效率動態曲線相比于方法2更接近于實驗的測量結果,特別是在推力方面,計算量也沒有明顯增加。基于等效電路計算得到的推力、效率變化曲線,與實驗測量結果大體上是吻合的,無論是在次級連續還是斷續區域,二者均有相同的變化趨勢,但以下兩個因素導致了一些偏差:①基于縱向氣隙磁通密度分布模型推導等效電路,未充分考慮橫向端部效應的影響;②對初級漏感等參數使用近似線性化處理的方法存在一定誤差。

圖10 直線感應牽引電機實驗臺

圖11 推力和效率變化曲線

4 結論

針對直線輪軌交通線路中初級經過次級感應板斷續路段這一特殊工況,本文綜合考慮了初、次級耦合區域長度變化和縱向端部效應的影響,基于氣隙磁通密度和去磁渦流的分布模型,構建了次級斷續時的等效電路。初級漏感、次級電阻等參數變化根據耦合區域線性長度進行近似化處理,基于縱向氣隙磁通密度的分布推導出勵磁支路電感和電阻的修正系數。根據等效電路計算得到推力和效率的動態變化曲線,與實驗臺測量結果進行對比,驗證了等效電路的準確性。

[1] 盧琴芬, 孔浩, 石佳蒙, 等. 高速列車永磁同步牽引電機基于聯合仿真模型的單脈沖控制研究[J]. 電工技術學報, 2015, 30(14): 61-66.

Lu Qinfen, Kong Hao, Shi Jiameng, et al. Research on single-pulse control of traction PMSM in high speed train based on co-simulation model[J]. Transa- ctions of China Electrotechnical Society, 2015, 30(14): 61-66.

[2] 朱進權, 葛瓊璇, 孫鵬琨, 等. 基于自抗擾的高速磁浮列車牽引控制策略[J]. 電工技術學報, 2020, 35(5): 1065-1074.

Zhu Jinquan, Ge Qiongxuan, Sun Pengkun, et al. Traction-system research of high-speed maglev based on active disturbance rejection control[J]. Transa- ctions of China Electrotechnical Society, 2020, 35(5): 1065-1074.

[3] 孫曉, 史黎明, 張志華, 等. 并聯雙邊直線感應電機次級定向控制[J]. 中國電機工程學報, 2019, 39(4): 1181-1189.

Sun Xiao, Shi Liming, Zhang Zhihua, et al. The secondary field oriented control of parallel connected double sided linear induction motors[J]. Proceedings of the CSEE, 2019, 39(4): 1181-1189.

[4] 周天豪, 陳磊, 祝長生, 等. 基于自適應變步長最小均方算法的磁懸浮高速電機不平衡補償[J]. 電工技術學報, 2020, 35(9): 1900-1911.

Zhou Tianhao, Chen Lei, Zhu Changsheng, et al. Unbalance compensation for magnetically levitated high-speed motors based on adaptive variable step size least mean square algorithm[J]. Transactions of China Electrotechnical Society, 2020, 35(9): 1900- 1911.

[5] 鄭曉欽, 徐杰, 陳春濤, 等. 超高速磁浮渦流裝置制動力的解析分析[J]. 電工技術學報, 2020, 35(9): 1891-1899.

Zheng Xiaoqin, Xu Jie, Chen Chuntao, et al. Analytical calculation of braking force of super high speed maglev eddy current device[J]. Transactions of China Electrotechnical Society, 2020, 35(9): 1891- 1899.

[6] 張邦富, 程明, 王颯颯, 等. 基于改進型代理模型優化算法的磁通切換永磁直線電機優化設計[J]. 電工技術學報, 2020, 35(5): 1013-1021.

Zhang Bangfu, Cheng Ming, Wang Sasa, et al. Optimal design of flux-switching permanent magnet linear machine based on improved surrogate-based optimization algorithm[J]. Transactions of China Electrotechnical Society, 2020, 35(5): 1013-1021.

[7] 李超, 魯軍勇, 馬偉明, 等. 電磁發射用多級混合儲能充電策略優化[J]. 電工技術學報, 2017, 32(13): 118-124.

Li Chao, Lu Junyong, Ma Weiming, et al. Charging strategy amelioration of multilevel hybrid energy storagefor electromagnetic launch[J]. Transactions of China Electrotechnical Society, 2017, 32(13): 118- 124.

[8] 王金兵, 沈艷霞. 基于增量模型的永磁同步直線電機魯棒預測電流控制[J]. 電力系統保護與控制, 2020, 48(8): 69-77.

Wang Jinbing, Shen Yanxia. Robust predictive current control for a permanent magnet synchronous linear motor based on an incremental model[J]. Power System Protection and Control, 2020, 48(8): 69-77.

[9] Hu Dong, Xu Wei, Dian Renjun, et al. Loss mini- mization control of linear induction motor drive for linear metros[J]. IEEE Transactions on Industrial Electronics, 2018, 65(9): 6870-6880.

[10] 呂剛, 孫守光, 馬云雙. 城軌交通中直線感應牽引電機的效率最優控制[J]. 電機與控制學報, 2009, 13(4): 490-495.

Lü Gang, Sun Shouguang, Ma Yunshuang. Efficiency optimal control of linear induction motor for urban rail transit[J]. Electric Machines and Control, 2009, 13(4): 490-495.

[11] Pai R M, Boldea I, Nasar S A. A complete equivalent circuit of a linear induction motor with sheet secondary[J]. IEEE Transactions on Magnetics, 1988, 24(1): 639-654.

[12] Gieras J F. Linear induction drives[M]. New York: Clarendon Press, 1994.

[13] 龍遐令. 直線感應電動機等效電路的通用推導方法[J]. 電工技術學報, 1993, 8(4): 55-60.

Long Xialing. A general method driving the equiva- lent circuit of linear induction motor[J]. Transactions of China Electrotechnical Society, 1993, 8(4): 55-60.

[14] Xu Wei, Zhu Jianguo, Zhang Yongchang, et al. Equivalent circuits for single-sided linear induction motors[J]. IEEE Transactions on Industry Appli- cations, 2010, 46(6): 2410-2423.

[15] Duncan J. Linear induction motor-equivalent-circuit model[J]. IEE Proceedings. Part B: Electric Power Applications, 1983, 130(1): 51-57.

[16] Lü Gang, Zeng Dihui, Zhou Tong. An advanced equivalent circuit model for linear induction motors[J]. IEEE Transactions on Industrial Electronics, 2018, 65(9): 7495-7503.

[17] 劉可安, 田紅旗, 劉勇. 直線感應電機過無次級感應板區檢測[J]. 浙江大學學報: 工學版, 2016, 50(6): 9-18.

Liu Kean, Tian Hongqi, Liu Yong. Detection of passing area of no secondary reaction plate for linear induction motor[J]. Journal of Zhejiang Univerisity: Engineering Science, 2016, 50(6): 9-18.

[18] 邸珺, 范瑜, 劉亞靜. 基于等效次級的直線感應電機的電磁分析與參數辨識[J]. 電工技術學報, 2017, 32(11): 145-154.

Di Jun, Fan Yu, Liu Yajing. Electromagnetic analysis and parameter estimation for the linear induction motor based on equivalent secondary[J]. Transactions of China Electrotechnical Society, 2017, 32(11): 145-154.

[19] 姜杰, 王學斌, 殷家敏, 等. 基于感應電動機復數簡化模型的參數辨識研究[J]. 電力系統保護與控制, 2014, 42(19): 87-92.

Jiang Jie, Wang Xuebin, Yin Jiamin, et al. Research on identification of induction motor based on its simplified complex quantity models[J]. Power System Protection and Control, 2014, 42(19): 87-92.

[20] 呂剛, 楊琛, 羅志昆. 考慮死區時間的直線感應電機參數靜止自整定方法[J]. 北京交通大學學報, 2020, 44(1): 91-97, 105.

Lü Gang, Yang Chen, Luo Zhikun. Parameter static self-tuning method of linear induction motor considering dead-time[J]. Journal of Beijing Jiaotong University, 2020, 44(1): 91-97, 105.

[21] Freeman E M. Travelling wave in induction machines: input impedance and equivalents[J]. Proceedings of the Institution of Electrical Engineers, 1968, 115(12): 1772-1776.

An Equivalent Circuit of Linear Induction Traction Motor with Discontinuous Secondary

(School of Electrical Engineering Beijing Jiaotong University Beijing 100044 China)

When linear induction motor is used as traction motor of rail transit vehicles, the secondary induction plate will be missing in turnouts, turns and other sections, and the equivalent circuit with continuous secondary is no longer applicable. For this special working condition, a equivalent circuit based on the new air gap flux density distribution model and considering secondary discontinuity is proposed. Firstly, the analysis model of linear induction motor with discontinuous secondary is established. According to the coupling conditions of the primary and secondary, the discontinuous working condition is divided into three stages. Then, the linear approximation method is used to correct the primary leakage inductance, excitation inductance and other parameters which vary with the length of the coupling region. Combined with the new air gap flux density distribution model and the discontinuous condition of the secondary, considering the influence of longitudinal end effect, the excitation inductance on the excitation branch and the equivalent resistance representing the eddy current loss of end effect are modified. The equivalent circuit of three stages is derived to represent the dynamic process of the primary passing through the discontinuous region of the secondary induction plate. Finally, the dynamic curve of the thrust and efficiency of the motor are obtained according to the equivalent circuit, and the correctness of the proposed equivalent circuit is verified by comparing with the experimental results.

Linear induction motor, discontinuous secondary, equivalent circuit, longitudinal end effect, air gap flux density

TM359.4

10.19595/j.cnki.1000-6753.tces.200804

國家自然科學基金資助項目(51777009,52077003)。

2020-07-07

2020-12-24

呂 剛 男,1976年生,博士,教授,研究方向為軌道交通直線電機與磁懸浮技術。E-mail: ganglv@bjtu.edu.cn(通信作者)

羅志昆 男,1996年生,碩士研究生,研究方向為直線電機及其控制。E-mail: 18121474@bjtu.edu.cn

(編輯 陳 誠)

猜你喜歡
區域
分割區域
探尋區域創新的密碼
科學(2020年5期)2020-11-26 08:19:22
基于BM3D的復雜紋理區域圖像去噪
軟件(2020年3期)2020-04-20 01:45:18
小區域、大發展
商周刊(2018年15期)2018-07-27 01:41:20
論“戎”的活動區域
敦煌學輯刊(2018年1期)2018-07-09 05:46:42
區域發展篇
區域經濟
關于四色猜想
分區域
公司治理與技術創新:分區域比較
主站蜘蛛池模板: 亚洲天堂久久| 被公侵犯人妻少妇一区二区三区| 毛片视频网| 五月天婷婷网亚洲综合在线| 福利姬国产精品一区在线| 国产一区在线视频观看| 欧美国产综合色视频| 精品人妻一区无码视频| 激情网址在线观看| 91色在线观看| 国产毛片高清一级国语 | 国产超薄肉色丝袜网站| 小蝌蚪亚洲精品国产| 青草免费在线观看| 18禁黄无遮挡网站| 五月天久久综合| 欧美区日韩区| 中文字幕精品一区二区三区视频| 99这里只有精品免费视频| 国产乱子伦视频三区| 一本色道久久88亚洲综合| 国产内射一区亚洲| 亚洲AV电影不卡在线观看| 波多野结衣中文字幕一区二区 | 国产另类乱子伦精品免费女| 国产网站免费| 亚洲精品第一页不卡| 欧美福利在线观看| 97在线国产视频| 日本午夜精品一本在线观看 | 亚洲香蕉在线| 中国毛片网| 成人毛片免费在线观看| 九九精品在线观看| 专干老肥熟女视频网站| 亚洲综合婷婷激情| 午夜丁香婷婷| 国产微拍一区二区三区四区| 国产小视频在线高清播放| 狠狠色丁香婷婷| 在线看AV天堂| 国产91视频观看| 这里只有精品免费视频| 九色国产在线| 国产精品精品视频| 国产无码精品在线| 欧美日韩亚洲国产| 免费欧美一级| 亚洲欧美另类中文字幕| 麻豆国产在线观看一区二区| 一本大道无码日韩精品影视| 国产二级毛片| 亚洲欧洲自拍拍偷午夜色| 国产成人亚洲精品蜜芽影院| 国产亚洲精久久久久久久91| 欧美色亚洲| 丝袜无码一区二区三区| 精品国产污污免费网站| 欧美一区二区三区欧美日韩亚洲 | 成人欧美日韩| 亚洲精品无码日韩国产不卡| 免费看一级毛片波多结衣| 亚洲一级无毛片无码在线免费视频 | 在线欧美日韩| 69精品在线观看| 亚洲三级色| 亚洲视频免费播放| 午夜老司机永久免费看片| 伦精品一区二区三区视频| 精品久久高清| 国产在线一区视频| 亚洲精品自拍区在线观看| 亚洲Av综合日韩精品久久久| 国产美女久久久久不卡| 亚洲国产系列| 欧美啪啪网| 久久99国产精品成人欧美| 蜜臀AVWWW国产天堂| 欧美午夜视频在线| 久久99国产精品成人欧美| 国产正在播放| 91久久偷偷做嫩草影院|