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

考慮氣動系統的高速受電弓分層控制

2022-05-20 07:44:50劉志剛
電工技術學報 2022年10期
關鍵詞:控制策略模型系統

楊 鵬 張 靜 金 偉 劉志剛

考慮氣動系統的高速受電弓分層控制

楊 鵬1張 靜2金 偉1劉志剛1

(1. 西南交通大學電氣工程學院 成都 611756 2. 西南交通大學機械工程學院 成都 610031)

0 引言

隨著高速列車運行速度的不斷提升,弓網的垂向振動、縱向沖擊、橫向擺動及耦合振動等問題越發顯著[1]。列車的動力來源于弓網滑動接觸過程傳遞的電能,穩定的弓網受流是保證高速列車正常運行的重要條件,弓網系統如圖1所示。

圖1 弓網系統

弓網系統的接觸波動不僅增加弓網間接觸磨損、離線燃弧等問題[2],也使得列車運行時的受流質量急劇惡化[3-5]。因此,有效降低弓網間的耦合振動,提高弓網受流質量,是當下亟須解決的問題之一。

目前,弓網受流質量的優劣主要采用弓網接觸力波動程度表征[6]。為降低接觸力波動,文獻[7-9]提出對受電弓采取參數優化的方法,但能降低的接觸力標準差極為有限。另外,文獻[10-11]開展接觸網吊弦、懸掛方式等結構優化方面的研究,建議對既有線路改造,但是改造工程耗資巨大。受電弓主動控制能夠解決優化設計存在的局限性問題,文獻[12-13]以接觸力的偏差和控制器的能量約束等信息,建立多個控制目標,分別采用模糊控制和PID控制策略,并通過弓網模型驗證了控制器的有效性。但是,上述研究皆僅選取弓網接觸力作為反饋信號,且接觸力信號獲取過程中常存在過多噪聲干擾。文獻[14]提出了一種阻抗控制器,選取受電弓框架的位移量為反饋信息,并與PID控制策略比較。在抑制接觸力波動的問題上,所提方法的控制效果更加明顯。文獻[15]以受電弓各結構的位移和速度為反饋控制量,設計了一種基于Backstepping方法的控制器,減小接觸力實際值與參考值之間的偏差。文獻[16]針對模型的不確定性和難以準確建模等問題,提出了一種模糊自適應終端滑模控制方法。在模型建模存在誤差的情況下,仍能夠保證良好的控制效果。文獻[17]基于擴展狀態變量下的受電弓非線性模型,采用微分幾何理論進行線性化處理,并增設干擾觀測器來提高控制系統的精確度。但在高速運行下,部分反饋量因傳感器靈敏度難以達到控制器的要求而造成數據缺失,影響到控制系統的有效性。因此,文獻[18]在采用二階滑模控制策略的基礎上,采用代數可觀測性理論,根據受電弓上、下框架的位移測量值估計接觸力,避免了直接對接觸力信號的采集。文獻[19]考慮到電磁干擾下的數據丟失問題,利用魯棒卡爾曼濾波器對受電弓狀態進行估計。但是,控制過程中出現時滯情況也會對控制效果造成一定影響。在此基礎上,文獻[20]采用最優控制策略,解決作動器輸入延遲和接觸力時變問題,該方法不僅能夠調整狀態和控制指標的相對權重,而且構建的穩定性判據兼顧系統穩定性。

綜上所述,現有高速受電弓控制策略方面的研究主要是根據受電弓控制器計算主動控制力大小,并直接施加于受電弓模型。雖然部分文獻中考慮了控制器造成的通信時滯和作動器輸入時滯問題,但并未涉及高速受電弓氣動系統動態特性的影響。在受電弓實際控制過程中,氣動系統達到期望壓強的響應時間遠大于通信時滯和輸入時滯此類純時滯情況,較長的響應時間將會嚴重影響控制器精度。因此,在研究受電弓主動控制策略時,考慮氣動系統的動態特性,以減少響應時間的影響具有重要理論意義和工程實用價值。

1 弓網及氣壓驅動系統模型

1.1 弓網系統模型

為更好表述受電弓的動態特性,在高速受電弓主動控制研究中,受電弓多是采用歸算質量模 型[21]。本文選取受電弓三質量塊模型,質量塊之間通過等效彈簧剛度和阻尼建立連接關系。受電弓三質量塊模型如圖2所示。

圖2 受電弓三質量塊模型

受電弓三質量塊的運動學方程為

接觸網結構較為復雜,包括接觸線、承力索、定位器和吊線等裝置。為建立更精確的模型,本文采用文獻[22]中建立的接觸網系統模型。使用有限元方法,完成接觸網的建模。接觸網有限元模型動力學方程為

其中

1.2 受電弓氣動系統模型

在受電弓主動控制中,受電弓控制器的輸出值以電信號傳遞給氣動系統,受電弓氣動系統原理如圖3所示。

在氣動系統工作時,兩位三通電磁換向閥開啟,使絕對壓力為0.8MPa的氣源裝置向氣動系統提供氣壓。當氣體到達氣控先導式精密調壓閥后,根據控制器傳遞的控制信號,將氣源壓強轉變為控制壓強,并經過降弓調速閥、絕緣軟管等輸送至氣囊內,完成氣囊內部壓力的調節,進而實現弓網接觸力的控制。而對于較高的氣壓值,安全閥會根據設定值進行溢流動作,保護整個系統安全運作。同時,為減少氣源的輸出壓力對氣控先導式精密調壓閥的沖擊,所調壓閥輸入壓力范圍為0.005~0.8MPa。根據文獻[23],假設氣動系統內部恒溫,根據氣體流量和器件之間的壓強差關系,建立氣動系統的數學模型。針對各器件的功能,對整個系統進行虛擬容器劃分,氣動系統模型如圖4所示。圖中,ADD支路為受電弓氣壓驅動系統,為保證受電弓在特殊情況下能快速自動降落而設計的。電氣比例閥根據控制回路產生的脈沖寬度調制(Pulse Width Modu- lation, PWM)電信號,由電磁閥控制產生控制壓強。該過程對整個氣動系統運行時的影響較小,可將電氣比例閥視為理想元件,忽略控制器到氣動系統之間信號轉變過程造成的延遲影響。

圖3 受電弓氣動系統原理

1—氣源 2—截止閥 3—過濾器 4—兩位三通電磁換向閥 5—升弓調速閥 6—電氣比例閥 7—壓力傳感器 8—氣控先導式精密調壓閥 9—降弓調速閥 10—安全閥 11—絕緣軟管 12—快速排氣閥 13—試驗閥

Fig.3 Principle of pneumatic system of pantograph

圖4 氣動系統模型

2 受電弓分層

圖5 分層控制結構

2.1 上層控制器設計

在設計上層控制器時,除了考慮弓網接觸力波動外,還分析了主動控制力和弓頭位移等因素的影響。因此,控制目標包含的具體性能指標選定為:

由于接觸網有限元模型的復雜性,需要對模型進行簡化,將接觸網有限元模型等效為面向控制的靜態剛度系數模型。目前,弓網間耦合關系主要通過罰函數關系建立[26],弓網之間的接觸力為

(1)閉環系統內部穩定。

針對閉環系統的約束條件,控制器的矩陣不等式滿足條件如下,其中“*”表示矩陣對稱位置變量的轉置。

2.2 下層控制器設計

針對氣動系統響應時間過長的問題,采用結構簡單且適應性強的PID控制器。但考慮到PID控制器的整定參數過多,為保證良好的控制性能,需經過大量的調試,這一過程極大地增加了控制器設計的復雜程度。因此,本文將內模理論應用到PID控制器設計中[29],對氣動系統采用IMC-PID控制策略。該方法將PID多參數調試難題簡化為單一參數整定問題。通過對濾波常數的整定,完成PID參數的調節,從而實現響應時間的控制。內模控制結構及其等效形式如圖6所示。

圖6 內模及其等效反饋結構

由于氣壓驅動系統的結構復雜,難以獲取整個系統的內模模型傳遞函數。因此,在考慮到控制器通信過程的時滯時間基礎上,將氣動系統模型化簡為一階時滯模型[30]為

將引入低通濾波器后的內模控制器和簡化的內模模型代入式(18),可得到系統的反饋控制器為

為方便控制器參數的整定,可將式(23)轉換為

相比于響應時間,氣壓驅動系統的純時滯時間較短。為保證氣動系統的穩定性,在IMC-PID控制器中不考慮微分環節,故忽略二階以上高階項。將式(24)用麥克勞林公式展開,內模控制器轉變為所需的PI控制器,具體表達式為

3 仿真分析

3.1 氣動系統響應控制仿真

為驗證下層控制器對縮短氣動系統響應時間的效果,建立氣動系統模型,參數見表1。

表1 氣動系統元件參數

Tab.1 Parameters of pneumatic system components

在安全閥設定的安全壓力值下,通過對模型進行仿真,受電弓氣囊壓強與靜態接觸壓力之間近似呈一定線性關系,其關系如圖7所示。在文獻[31]中通過對氣壓驅動系統的試驗數據采集擬合,氣囊壓強對應靜態接觸壓力間的大致線性關系進一步驗證了仿真結論的可靠性。

圖7 靜態接觸壓力與氣囊壓強關系

比較簡化模型與氣動系統的階躍響應,氣壓驅動系統和簡化模型響應輸出曲線具體如圖8所示。對于不同的階躍信號,兩者之間的響應曲線都始終保持大致重合。因此,氣動系統的簡化模型滿足內模控制傳遞函數的要求,可替代氣動系統模型進行下層控制器設計。

圖9 氣動系統輸出曲線

仿真結果表明,當控制信號改變時,相比于無控制狀態下的輸出曲線,采用IMC-PID控制的氣動系統到達控制目標的響應時間明顯縮短,顯著加快了氣動系統內部的響應速度,且調節過程中幾乎無超調量產生,也不會對系統的穩定性造成過多的影響,有效改善了輸出壓強滯后問題。

3.2 弓網控制仿真分析

表2 DSA380受電弓質量塊模型參數

Tab.2 DSA380 pantograph mass model parameters

表3 京津線接觸網結構參數

Tab.3 Catenary structure parameters of Beijing-Tianjin line

為驗證受電弓主動控制器的普適性,分別在250km/h、300km/h和350km/h 3個運行工況下,對接觸網5跨內的受電弓控制效果進行了仿真證明。比較無控制、傳統單層控制策略和分層控制策略的控制效果,弓網接觸力控制結果如圖10所示。對于不同控制方法的弓網接觸力標準差的計算值,具體結果見表4。

圖10 各時速的接觸力仿真結果

表4 不同控制策略的弓網接觸壓力標準差

Tab.4 Standard deviation of pantograph catenary contact pressure with different control strategies

由圖10可知,根據本文提出的控制策略所設計的傳統控制器和分層控制器,都降低了弓網接觸力的標準差。由表4可知,傳統單層控制策略雖然能降低接觸力波動,但控制效果并不明顯,這是因為氣動系統自身的響應時間過長,導致氣囊壓強變化難以跟蹤單層控制器計算出的目標控制值。而在采用分層控制策略后,下層控制器的作用使響應時間明顯減少,加快了氣動系統的響應速度。圖11為各速度的接觸力標準差下降百分比。

圖11 各速度的接觸力標準差下降百分比

由圖11可知,針對不同運行時速工況,在考慮了氣動系統的動態特性影響后,相比傳統單層控制的效果,受電弓分層控制的接觸力標準差下降百分比提高10%左右,表明了分層控制方法在解決弓網耦合振動問題上具有明顯優越性。針對弓網接觸的弓頭位移指標,具體仿真結果如圖12所示。

圖12 弓頭位移波動

兩種控制策略對弓頭位移均有一定程度減小。但分層控制下的控制量變化實時性更快,對弓頭位移的約束效果優于傳統控制策略。因此,弓頭位移的仿真結果進一步證明了受電弓分層控制策略的有效性。

4 結論

本文針對受電弓氣動系統中存在的響應時間較長的難題,在傳統單層受電弓控制的基礎上引入氣動系統控制器,提出了一種基于分層結構的受電弓主動控制策略。在考慮氣動系統動態響應影響的情況下,通過非線性弓網模型仿真計算,驗證了分層控制策略降低接觸力波動的有效性,可得到以下結論:

2)氣動系統控制采用內模理論,根據簡化后的氣動系統模型,設計了內模控制器,并利用麥克勞林展開公式得到了僅與單一參數相關的PID控制器。在IMC-PID控制下,氣動系統的響應時間明顯縮短,大幅提升了對控制目標值的跟蹤速率。

[1] 張靜, 劉志剛, 魯小兵, 等. 高速弓網空氣動力學研究進展[J]. 鐵道學報, 2015, 37(1): 7-15.

Zhang Jing, Liu Zhigang, Lu Xiaobing, et al. Study on aerodynamics development of high-speed pantograph and catenary[J]. Journal of the China Railway Society, 2015, 37(1): 7-15.

[2] 陳忠華, 唐俊, 時光, 等. 弓網強電流滑動電接觸摩擦振動分析與建模[J]. 電工技術學報, 2020, 35(18): 3869-3877.

Chen Zhonghua, Tang Jun, Shi Guang, et al. Analysis and modeling of high current sliding electrical contact friction dynamics in pantograph-catenary system[J]. Transactions of China Electrotechnical Society, 2020, 35(18): 3869-3877.

[3] Lee J H, Kim Y G, Paik J S, et al. Performance evaluation and design optimization using differential evolutionary algorithm of the pantograph for the high-speed train[J]. Journal of Mechanical Science and Technology, 2012, 26(10): 3253-3260.

[4] 劉方林. 電氣化鐵路動態弓網接觸電阻研究[J]. 電氣技術, 2018, 19(9): 69-72.

Liu Fanglin. Study on dynamic contact resistance between pantograph and catenary in electrified rail- way[J]. Electrical Engineering, 2018, 19(9): 69-72.

[5] 陳忠華, 吳迪, 回立川, 等. 波動載荷下弓網滑動電接觸失效研究[J]. 電工技術學報, 2019, 34(21): 4492-4500.

Chen Zhonghua, Wu Di, Hui Lichuan, et al. Research on failure of pantograph-catenary sliding electrical contact under fluctuation load[J]. Transactions of China Electrotechnical Society, 2019, 34(21): 4492- 4500.

[6] Wang Hongrui, Liu Zhigang, Song Yang, et al. Dete- ction of contact wire irregularities using a quadratic time-frequency representation of the pantograph- catenary contact force[J]. IEEE Transactions on Instrumentation and Measurement, 2016, 65(6): 1385- 1397.

[7] Kim J W, Yu S N. Design variable optimization for pantograph system of high-speed train using robust design technique[J]. International Journal of Precision Engineering and Manufacturing, 2013, 14(2): 267- 273.

[8] 周寧. 350km/h及以上弓網動態行為研究[D]. 成都: 西南交通大學, 2013.

[9] 吳孟臻, 劉洋, 許向紅. 高速弓網系統動力學參數敏度分析及優化[J]. 力學學報, 2021, 53(1): 75-83.

Wu Mengzhen, Liu Yang, Xu Xianghong. Sensitivity analysis and optimization on parameters of high speed pantograph-catenary system[J]. Chinese Journal of Theoretical and Applied Mechanics, 2021, 53(1): 75-83.

[10] 蔣先國, 古曉東, 鄧洪, 等. 基于微動理論的整體吊弦損傷機理及優化研究[J]. 鐵道學報, 2019, 41(6): 40-45.

Jiang Xianguo, Gu Xiaodong, Deng Hong, et al. Research on damage mechanism and optimization of integral dropper string based on fretting theory[J]. Journal of the China Railway Society, 2019, 41(6): 40-45.

[11] Zhang Weihua, Zhou Ning, Li Ruiping, et al. Panto- graph and catenary system with double pantographs for high-speed trains at 350 km/h or higher[J]. Journal of Modern Transportation, 2011, 19(1): 7-11.

[12] 時光, 劉健辰, 陳忠華, 等. 基于DE-EDA多目標優化的受電弓模糊控制[J]. 計算機工程與應用, 2016, 52(1): 229-232, 253.

Shi Guang, Liu Jianchen, Chen Zhonghua, et al. Control of pantograph based on DE-EDA multi- objective optimization[J]. Computer Engineering and Applications, 2016, 52(1): 229-232, 253.

[13] Cristina S R, Jimenez J O, Carnicero A. Active control strategy on a catenary-pantograph validated model[J]. Vehicle System Dynamics, 2013, 51(4): 554-569.

[14] Yoshitaka Y, Mitsuru I. Advanced active control of contact force between pantograph and catenary for high-speed trains[J]. Quarterly Report of RTRI, 2012, 53(1): 28-33.

[15] Zhu Bing, Ren Zhiling, Xie Wenjing, et al. Active nonlinear partial-state feedback control of contacting force for a pantograph-catenary system[J]. ISA Transa- ctions, 2019, 91(1): 78-89.

[16] 時光, 陳忠華, 郭鳳儀, 等. 基于最優載荷的受電弓自適應終端滑模控制[J]. 電工技術學報, 2017, 32(4): 140-146, 153.

Shi Guang, Chen Zhonghua, Guo Fengyi, et al. Adaptive terminal sliding mode control of pantograph based on optimal load[J]. Transactions of China Elec- trotechnical Society, 2017, 32(4): 140-146, 153.

[17] 時光, 陳忠華, 郭鳳儀, 等. 弓網接觸力反饋線性化控制[J]. 控制理論與應用, 2016, 33(1): 85-91.

Shi Guang, Chen Zhonghua, Guo Fengyi, et al. Feedback linearization control of load between panto- graph and catenary[J]. Control Theory and Appli- cations, 2016, 33(1): 85-91.

[18] Pisano A, Usai E. Contact force estimation and regulation in active pantographs: an algebraic obser- vability approach[J]. Asian Journal of Control, 2011, 13(6): 761-772.

[19] 張靜, 宋寶林, 謝松霖, 等. 基于狀態估計的高速受電弓魯棒預測控制[J]. 電工技術學報, 2021, 36(5): 1075-1083.

Zhang Jing, Song Baolin, Xie Songlin. Robust predictive control of high-speed pantograph based on state estimation[J]. Transactions of China Electro- technical Society, 2021, 36(5): 1075-1083.

[20] Lin Yuchen, Shieh N C, Liu V T. Optimal control for rail vehicle pantograph systems with actuator delays[J]. IET Control Theory and Applications, 2015, 9(13): 1917-1926.

[21] Lu Xiaobing, Liu Zhigang, Song Yang, et al. Estimator- based multi-objective robust control strategy for an active pantograph in high-speed railways[J]. Pro- ceedings of the Institution of Mechanical Engineers Part F-Journal of Rail and Rapid Transit, 2018, 232(4): 1064-1077.

[22] Song Yang, Liu Zhigang, Wang Hongrui, et al. Nonlinear modelling of high-speed catenary based on analytical expressions of cable and truss elements[J]. Vehicle System Dynamics, 2015, 53(10): 1455-1479.

[23] 任輝文. 氣壓驅動系統特性對受電弓動力學的影響[D]. 長沙: 湖南大學, 2019.

[24] 黃文龍, 胡海濤, 陳俊宇, 等. 樞紐型牽引變電所再生制動能量利用系統能量管理及控制策略[J]. 電工技術學報, 2021, 36(3): 588-598.

Huang Wenlong, Hu Haitao, Chen Junyu, et al. Energy management and control strategy of regener- ative braking energy utilization system in hub traction substation[J]. Transactions of China Electrotechnical Society, 2021, 36(3): 588-598.

[25] 郭偉, 趙洪山. 基于事件觸發機制的直流微電網多混合儲能系統分層協調控制方法[J]. 電工技術學報, 2020, 35(5): 1140-1151.

Guo Wei, Zhao Hongshan. Coordinated control method of multiple hybrid energy storage system in DC microgrid based on event-triggered mechanism[J]. Transactions of China Electrotechnical Society, 2020, 35(5): 1140-1151.

[26] 王勖成. 有限單元法[M]. 北京: 清華大學出版社, 2003.

[27] 楊京, 王彤, 畢經天, 等. 含直驅風電機組的電力系統次同步振蕩魯棒阻尼控制[J]. 電力系統自動化, 2020, 44(3): 56-65.

Yang Jing, Wang Tong, Bi Jingtian, et al. Robust damping control of subsynchronous oscillation in power system with direct-drive wind turbines[J]. Automation of Electric Power Systems, 2020, 44(3): 56-65.

[29] 周天豪, 楊智, 祝長生, 等. 電磁軸承高速電機轉子系統的內模-PID控制[J]. 電工技術學報, 2020, 35(16): 3414-3425.

Zhou Tianhao, Yang Zhi, Zhu Changsheng, et al. Internal model control-PID control of an active magnetic bearing high-speed motor rotor system[J]. Transactions of China Electrotechnical Society, 2020, 35(16): 3414-3425.

[30] Jin Qibing, Liu Qi. IMC-PID design based on model matching approach and closed-loop shaping[J]. ISA Transaction, 2014, 53(2): 462-473.

[31] 苗海濤. 基于空氣彈簧高速受電弓的氣壓伺服主動控制研究[D]. 成都: 西南交通大學, 2009.

[32] 寇發榮, 李冬, 許家楠, 等. 車輛電動靜液壓主動懸架內模PID控制研究[J]. 液壓與氣動, 2018(6): 1-7.

Kou Farong, Li Dong, Xu Jianan, et al. Study on internal model PID control of vehicle active suspen- sion with electro-hydrostatic actuator[J]. Chinese Hydraulics and Pneumatics, 2018(6): 1-7.

Hierarchical Control of High-Speed Pantograph Considering Pneumatic System

1211

(1. School of Electrical Engineering Southwest Jiaotong University Chengdu 611756 China 2. School of Mechanical Engineering Southwest Jiaotong University Chengdu 610031 China)

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

TM571

楊 鵬 男,1996年生,碩士,研究方向為受電弓主動控制。E-mail: 18234136381@163.com

張靜 女,1979年生,副教授,碩士生導師,研究方向為高速弓網動力學特性、參數優化及匹配及控制技術、外骨骼機器人優化設計及控制技術。E-mail: sdzj2006@126.com(通信作者)

2021-04-19

2021-08-25

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

(編輯 陳 誠)

猜你喜歡
控制策略模型系統
一半模型
Smartflower POP 一體式光伏系統
工業設計(2022年8期)2022-09-09 07:43:20
WJ-700無人機系統
考慮虛擬慣性的VSC-MTDC改進下垂控制策略
能源工程(2020年6期)2021-01-26 00:55:22
ZC系列無人機遙感系統
北京測繪(2020年12期)2020-12-29 01:33:58
重要模型『一線三等角』
重尾非線性自回歸模型自加權M-估計的漸近分布
工程造價控制策略
山東冶金(2019年3期)2019-07-10 00:54:04
現代企業會計的內部控制策略探討
消費導刊(2018年10期)2018-08-20 02:57:02
連通與提升系統的最后一塊拼圖 Audiolab 傲立 M-DAC mini
主站蜘蛛池模板: 久久久久无码精品| 国产浮力第一页永久地址 | 国产丝袜无码精品| 国内老司机精品视频在线播出| 欧美激情综合一区二区| 97国产精品视频人人做人人爱| 欧美视频在线不卡| 久久视精品| 亚洲精品视频网| 性喷潮久久久久久久久| 免费一级全黄少妇性色生活片| 国产xx在线观看| 玖玖精品在线| 国产成人喷潮在线观看| 午夜国产不卡在线观看视频| 国产电话自拍伊人| 成人午夜福利视频| 国模粉嫩小泬视频在线观看| 久久91精品牛牛| 青青草91视频| 亚洲天堂在线免费| a级毛片网| 伦精品一区二区三区视频| 欧美伦理一区| 全部免费毛片免费播放| 性视频久久| 日本精品视频| 91成人在线观看| 国产自无码视频在线观看| 三上悠亚在线精品二区| 精品视频福利| 伊人福利视频| 久久9966精品国产免费| 91久久偷偷做嫩草影院电| 国产成人av一区二区三区| 玩两个丰满老熟女久久网| 欧美翘臀一区二区三区| 国产精品视频猛进猛出| 欧美日本激情| 国产三级毛片| 毛片网站观看| 国产成人麻豆精品| 国产精品污污在线观看网站| 日韩视频免费| 久久免费观看视频| 成人韩免费网站| 久久中文电影| 人人妻人人澡人人爽欧美一区| AV无码无在线观看免费| 国产高颜值露脸在线观看| 99国产精品国产| 欧美性猛交xxxx乱大交极品| 欧美日韩中文国产va另类| 人妻丰满熟妇av五码区| 中文字幕在线观看日本| 亚洲国产日韩在线成人蜜芽| 久久久久国产精品免费免费不卡| 久草热视频在线| 亚洲欧洲日产无码AV| 国产成人做受免费视频| 全午夜免费一级毛片| 国产在线观看人成激情视频| 奇米精品一区二区三区在线观看| 天天躁夜夜躁狠狠躁躁88| 午夜精品影院| 中国一级特黄大片在线观看| 最新日韩AV网址在线观看| 韩国v欧美v亚洲v日本v| 一级一级一片免费| 天堂在线www网亚洲| 国产自在线播放| 国产夜色视频| a毛片在线| 亚洲午夜福利精品无码不卡| 免费可以看的无遮挡av无码| 综合色在线| 美女被操91视频| 欧美翘臀一区二区三区| 国产成人高精品免费视频| 国产在线98福利播放视频免费 | 国产又色又刺激高潮免费看| 大学生久久香蕉国产线观看|