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

基于變論域模糊控制的車輛半主動懸架控制方法

2017-03-01 09:27:19王大勇
中國機械工程 2017年3期
關(guān)鍵詞:振動模型

王大勇 王 慧

遼寧工程技術(shù)大學(xué)機械工程學(xué)院,阜新,123000

基于變論域模糊控制的車輛半主動懸架控制方法

王大勇 王 慧

遼寧工程技術(shù)大學(xué)機械工程學(xué)院,阜新,123000

針對車輛半主動懸架系統(tǒng),提出了一種基于變論域模糊PID控制方法,目標(biāo)是提高車輛在隨機路面激勵作用下的平順性。通過將變論域方法與模糊PID控制器相結(jié)合來解決模糊PID存在的因模糊規(guī)則制定盲目性而產(chǎn)生的在線調(diào)節(jié)時間過長等問題。由仿真和實驗研究對比可知,變論域模糊PID控制下的半主動懸架系統(tǒng)中的車身垂直振動速度和加速度比常規(guī)PID控制下的車身垂直振動速度和加速度分別減小了46.56%和29.21%,相比被動懸架系統(tǒng)的車身垂直振動速度和加速度分別減小了58.05%和49.74%。使用該車輛半主動懸架模糊控制方法可提高車輛的平順性。

半主動懸架;變論域;模糊PID算法;車輛振動

0 引言

隨著科學(xué)技術(shù)和社會的發(fā)展,人們對車輛行駛過程中的平順性和舒適性提出了更高的要求。車輛懸架系統(tǒng)是影響車輛平順性和舒適性最重要的部件[1-2]。經(jīng)過長時間的探索和研究,人們研究出了通過調(diào)節(jié)阻尼大小從而減小車輛振動的主動懸架系統(tǒng)和半主動懸架系統(tǒng)。

文獻[3]建立了主動懸架二自由度動力模型,并將白噪聲和正弦信號作為路面激勵,使用PID控制器對主動懸架系統(tǒng)進行控制,仿真得到主動懸架作用下車輛的振動撓度和加速度響應(yīng)。文獻[4-5]使用模糊PID控制器對主動懸架系統(tǒng)進行控制,解決了常規(guī)PID控制器的缺陷,通過MATLAB仿真和實驗方法對所研究的主動懸架控制方法的優(yōu)勢進行了驗證。文獻[6-7]使用遺傳算法對模糊PID控制器的量化因子和參數(shù)修正系數(shù)進行優(yōu)化,期望得到最優(yōu)的控制器參數(shù),通過建立的主動懸架模型進行仿真和實驗分析。文獻[8]建立了半主動懸架的1/4車輛模型,使用PID閉環(huán)控制方法進行阻尼力的實時調(diào)節(jié),實現(xiàn)了車輛懸架的半主動控制,通過仿真方法對PID閉環(huán)控制方法進行研究。文獻[9]建立了半主動懸架的1/2車輛模型,并使用模糊PID控制算法進行懸架阻尼和剛度調(diào)節(jié),實現(xiàn)了半主動懸架控制。文獻[10]在ADAMS軟件中建立了車輛的半主動懸架模型,使用模糊控制器作為半主動懸架的控制算法,并使用遺傳算法對模糊控制器的參數(shù)進行優(yōu)化,以C級路面作為路面激勵,仿真研究結(jié)果表明,使用該控制方法可明顯地提高汽車的行駛平順性。

常規(guī)PID控制器已不能滿足車輛懸架系統(tǒng)這類復(fù)雜系統(tǒng)的控制要求,而模糊PID控制器的固定基本論域和參數(shù)會直接影響模糊PID控制器的控制性能。相對于主動懸架系統(tǒng),半主動懸架系統(tǒng)具有結(jié)構(gòu)簡單、易于維護、能耗低等優(yōu)勢。在建立車輛懸架模型時,1/4模型能夠簡化求解過程,降低計算量。本文針對車輛半主動懸架系統(tǒng)建立1/4車輛簡化模型,提出了一種基于變論域模糊PID控制方法,目標(biāo)是提高車輛在隨機路面激勵作用下的平順性。

1 車輛半主動懸架系統(tǒng)數(shù)學(xué)模型

1.1 路面激勵數(shù)學(xué)模型

路面的隨機激勵是研究車輛懸架系統(tǒng)動力學(xué)行為的主要外部激勵之一。路面不平度是十分復(fù)雜的隨機過程,許多文獻對此進行了報道,隨機路面輪廓的生成一般使用積分高斯白噪聲或者通過對高斯白噪聲得到。使用比較廣泛的隨機路面輪廓輸入表達式為

(1)

式中,f0為濾波器截止頻率下限,當(dāng)f0=0時,路面的隨機輪廓用積分白噪聲表示;n0為參考空間頻率;Gq(n0)為不平度系數(shù);v為車輛速度;w(t)為高斯白噪聲。

當(dāng)車輛速度v為定值時,速度功率譜為定值。因此,可用積分后速度功率譜密度的高斯白噪聲來表示隨機的路面激勵[11]:

(2)

1.2 車輛懸架數(shù)學(xué)模型

1.2.1 車輛懸架模型的簡化

為了簡化求解過程,降低計算量,通常將七自由度的整車懸架模型進行簡化,在滿足特定條件的情況下,先簡化為1/2懸架模型,進而簡化為1/4模型。首先對懸架系統(tǒng)作如下假設(shè):①將輪胎簡化為等效剛度的彈簧,忽略其阻尼的影響;②將車軸和車輛視為剛性體,車輪和路面視為點接觸;③忽略車輛內(nèi)部振動的影響,即將車身與車座等部件的連接視為剛性連接;④假設(shè)路面對車輛左右兩側(cè)的垂直激勵是相同的[12]。

本文使用1/4車輛懸架模型進行研究,車輛懸架模型如圖1所示,其中虛線框內(nèi)即1/4車輛懸架模型。

圖1 車輛懸架模型Fig.1 Vehicle suspension model

為了使簡化后的模型計算不失真,需要使簡化模型滿足特定的條件,將車輛模型簡化為1/4模型需滿足下式:

(3)

式中,m為半車身(左半車身或右半車身)簧載質(zhì)量;a為車輛質(zhì)心到前車輪距離;b為車輛質(zhì)心到后車輪距離;I為半車身(左半車身或右半車身)俯仰運動的轉(zhuǎn)動慣量;mA為半車身的前1/4懸架的簧載質(zhì)量;mC為半車身的后1/4懸架的簧載質(zhì)量[13]。

1.2.2 1/4車輛半主動懸架動力學(xué)模型

1/4車輛的半主動懸架簡化模型如圖2所示。圖2中,mw為非彈簧負載質(zhì)量;Ks為半主動懸架的剛度;Kt為輪胎的等效剛度;xw為非彈簧負載位移;FMR為磁流變減振器的可控阻尼力,即本文進行半主動懸架控制的對象;ce為半主動懸架的阻尼系數(shù);xb為彈簧負載的位移;xg為路面激勵的位移;mb為彈簧負載等效質(zhì)量。

圖2 1/4車輛的半主動懸架簡化模型Fig.2 Simplified model of semi-active suspension for 1/4 vehicle

對1/4車輛的半主動懸架簡化模型(圖2)進行動力學(xué)分析,可得其動力學(xué)微分方程:

(4)

(5)

本文建立的狀態(tài)方程為

(6)

式中,X(t)為狀態(tài)變量;A為狀態(tài)矩陣;B為輸入矩陣;C為輸出矩陣;D為狀態(tài)反饋矩陣;F為擾動輸出矩陣;W(t)為擾動輸入矩陣;U(t)為控制矩陣,即本文控制方法的控制對象。

A=

B=[1/mb-1/mw0 0 0]T

D=[1/mb0 0]TW(t)=[w(t)]

2 變論域模糊PID控制

模糊PID控制器中設(shè)置有固定的量化和比例因子,以及固定的基本論域,然而固定的基本論域會直接影響模糊PID控制器的控制性能。假設(shè)當(dāng)輸入基本論域過大而輸入量較小時,大部分模糊規(guī)則閑置而沒有被利用,當(dāng)輸入基本論域過小而輸入量較大時,控制器的輸出調(diào)整量無法滿足系統(tǒng)需求,從而影響模糊PID控制器的控制性能。對于輸出論域同理。因此,本文使用變論域模糊PID控制算法對半主動懸架進行控制,變論域模糊PID控制器原理如圖3所示。

圖3 變論域模糊PID控制器原理框圖Fig.3 Diagram of variable universe fuzzy PID controller

變論域模糊PID控制器即在模糊PID控制器基礎(chǔ)上增加一個論域調(diào)整模塊,論域調(diào)整模塊的輸入為系統(tǒng)被控對象反饋的誤差和誤差變化率,輸出為對輸入和輸出基本論域的調(diào)整比例系數(shù),從而能夠在不改變輸入和輸出基本論域的情況下調(diào)節(jié)論域范圍,提高控制器的自適應(yīng)能力[16]。

隨著誤差和誤差變化率的不斷變化,論域也隨之變化,變化后的論域表示為

式中,αi(t)為輸入論域的伸縮因子;β(t)為輸出論域的伸縮因子;Ei為第i個輸入的論域;U為輸出的論域。

輸入量誤差e和誤差變化率ec的論域伸縮因子表示為

(7)

式中,τ1、τ2為0~1之間的常數(shù);E為輸入誤差的論域;EC為輸入誤差變化率的論域。

輸出論域伸縮因子表示為

(8)

式中,KI為比例因子;β(0) 為輸出論域伸縮因子的初值;pi為輸入變量權(quán)重常數(shù)向量,pi=[p1p2…pi];ei(τ)為誤差矢量。

變論域的模糊PID控制器的總輸出為

(9)

由此,將變論域與常規(guī)的模糊PID控制器進行結(jié)合,通過變論域方法自適應(yīng)調(diào)節(jié)論域,從而提高模糊PID控制器輸出參數(shù)的選取精確性,提高模糊PID控制器的自適應(yīng)能力[17]。

3 仿真與實驗分析

本文通過仿真和實驗的方法對基于變論域模糊PID車輛半主動懸架模糊控制方法進行分析。使用MATLAB中的Simulink仿真模塊建立仿真模型,其中隨機路面輸入仿真模型如圖4所示。

圖4 隨機路面輸入模塊Fig.4 Random road input module

為研究不同路面激勵作用下使用不同控制方法的作用效果,以B級和C級路面作為路面激勵類型。B級、C級路面的不平度系數(shù)Gq(n0)分別取64×10-6m2/m-1和2.56×10-6m2/m-1。研究兩種路面激勵以及40m/s、60m/s車速情況下的車身和速度響應(yīng)。

設(shè)置本文研究的車輛半主動懸架系統(tǒng)參數(shù)如下:車輛的彈簧負載等效質(zhì)量mb=350kg;半主動懸架的剛度Ks=18kN/m;半主動懸架的阻尼系數(shù)ce=550N·s/m;非彈簧負載質(zhì)量mw=50kg;輪胎的等效剛度Kt=200kN/m。

PID和模糊控制器中參數(shù)設(shè)置如下:PID的初始參數(shù)KP0=2.1,KI0=0.05,KD0=0.5;τ1和τ2均取0.5;輸入和輸出論域均歸一化到[-1,1]范圍內(nèi)。

仿真得到車輛以 40m/s速度在B級和C級路面上行駛時車身的速度和加速度響應(yīng),如圖5和圖6所示。

(a)車身的速度響應(yīng)

(b)車身的加速度響應(yīng)圖5 B級路面車身的速度和加速度響應(yīng)仿真曲線Fig.5 Speed and acceleration response simulation curve of B level road vehicle body

(a)車身的速度響應(yīng)

(b)車身的加速度響應(yīng)圖6 C級路面車身的速度和加速度響應(yīng)仿真曲線Fig.6 Speed and acceleration response simulation curve of C level road vehicle body

以表格數(shù)據(jù)方式列出B級和C級路面激勵以及40m/s、60m/s車速情況,不同控制方法作用下的車身垂直振動速度和加速度的均方根值見表1。

表1 不同懸架系統(tǒng)和控制方法作用下車身垂直振動仿真響應(yīng)對比Tab.1 Comparison of vertical vibration response of vehicle under different suspension systems and control methods

使用實驗的方法對基于變論域模糊PID車輛半主動懸架模糊控制方法進行研究,實驗裝置簡圖見圖7。懸架系統(tǒng)的參數(shù)與仿真模型參數(shù)一致,采用康泰電子的ULT2051/V型加速度傳感器,工控機采用美國 PXI-6509型NI 工控機,運動控制卡采用NI PCI-7390型運動卡,數(shù)據(jù)采集卡采用NI PCI-6221 型采集卡[18]。

圖7 實驗裝置簡圖Fig.7 Schematic diagram of experimental apparatus

使用實驗方法得到車輛以40 m/s速度在B級和C級路面上行駛時車身的速度和加速度響應(yīng),如圖8和圖9所示。

以表格數(shù)據(jù)方式列出B級和C級路面激勵以及40 m/s、60 m/s車速情況,不同控制方法作用下的車身垂直振動速度和加速度實測數(shù)據(jù)的均方根值見表2。

(a)車身的速度響應(yīng)

(b)車身的加速度響應(yīng)圖8 B級路面車身的速度和加速度響應(yīng)實驗曲線Fig.8 Speed and acceleration response test curve of B level road vehicle body

(a)車身的速度響應(yīng)

(b)車身的加速度響應(yīng)圖9 C級路面車身的速度和加速度響應(yīng)實驗曲線Fig.9 Speed and acceleration response test curve of C level road vehicle body

車速(m/s)監(jiān)測量(均方根值)被動懸架PID控制半主動懸架模糊控制半主動懸架B級路面40振動速度0.04150.03010.0169振動加速度0.41030.30820.198960振動速度0.05160.03920.0243振動加速度0.51220.38610.3281C級路面40振動速度0.04920.03520.0198振動加速度0.42950.31090.216260振動速度0.05620.04310.0265振動加速度0.56180.41560.3654

由仿真和實驗對比可知,兩者的數(shù)據(jù)趨勢基本一致;比較B級和C級路面激勵以及40 m/s、60 m/s車速情況下不同控制方法作用下的車身垂直振動速度和加速度的均方根值可知,變論域模糊PID控制下的半主動懸架系統(tǒng)中的車身垂直振動速度和加速度比常規(guī)PID控制下的振動速度和加速度分別減小了46.56%和29.21%,相比被動懸架系統(tǒng)的振動速度和加速度分別減小了58.05%和49.74%。這說明使用本文提出的車輛半主動懸架模糊控制方法可提高車輛的平順性。

4 結(jié)語

使用1/4車輛模型建立半主動懸架的動力學(xué)模型,求解簡便,計算量降低,且不會使計算結(jié)果失真。提出了一種基于變論域模糊PID的模糊控制方法,目標(biāo)是提高車輛在隨機路面激勵作用下的平順性。由仿真和實驗對比可知,兩者的數(shù)據(jù)趨勢基本一致。

[1] 喻凡, 林逸. 汽車系統(tǒng)動力學(xué)[M].北京: 機械工業(yè)出版社, 2005. YU Fan, LIN Yi.Automotive System Dynamics [M].Beijing: Machinery Industry Press, 2005.

[2] LIU Yanqing, MATSUHISA H, UTSUNO H. Semi-active Vibration Isolation System with Variable Stiffness and Damping Control[J]. Journal of Sound and Vibration,2008,313(1/2): 16-28.

[3] 杜常清,常曉瑞. PID控制策略下主動懸架系統(tǒng)的動態(tài)仿真[J].武漢理工大學(xué)學(xué)報(信息與管理工程版),2015,37(6):680-683. DU Changqing, CHANG Xiaorui. Dynamic Simulation of Active Suspension System Based on PID Control Strategy[J].Journal of Wuhan University of Technology(Information &Management Engineering),2015,37(6):680-683.

[4] 劉海潮,劉夫云,張驥,等.剛?cè)狁詈戏蔷€性主動懸架系統(tǒng)的FUZZY-PID控制[J].機械研究與應(yīng)用,2016,29(1):231-234. LIU Haichao, LIU Fuyun, ZHANG Ji, et al. FUZZY-PID Control of Nonlinear Active Suspension for Rigid-flexible Coupling Model[J].Mechanical Research & Application,2016, 29(1):231-234.

[5] 陳學(xué)文,曹曉政,任治.1/2汽車主動懸架模糊控制器設(shè)計及Simulink仿真研究[J].遼寧工業(yè)大學(xué)學(xué)報(自然科學(xué)版),2016,36(2):117-121. CHEN Xuewen, CAO Xiaozheng, REN Zhi.Research on Fuzzy Controller Design and Simulink Simulation of a Half Automobile Active Suspension[J].Journal of Liaoning University of Technology(Natural Science Edition) ,2016,36 (2):117-121.

[6] 鄭泉, 陳黎卿, 王繼先, 等. 基于ADAMS/Car 和Simulink 的主動懸架遺產(chǎn)模糊控制[J]. 農(nóng)業(yè)機械學(xué)報,2009, 40(5): 7-11. ZHENG Quan, CHEN Liqing, WANG Jixian, et al. Control of ASS Genetic Algorithms and Fuzzy Based on ADAMS/Car and Simulink [J]. Transactions of the Chinese Society for Agricultural Machinery, 2009, 40(5): 7-11.

[7] 趙強,何法,王鑫,等.基于遺傳算法優(yōu)化的車輛主動懸架模糊PID控制[J].重慶理工大學(xué)學(xué)報(自然科學(xué)),2016,30(2):6-11. ZHAO Qiang, HE Fa, WANG Xin, et al. Fuzzy-PID Control of Vehicle Active Suspension Based on Genetic Algorithm Optimization[J].Journal of Chongqing University of Technology( Natural Science) ,2016,30(2):6-11.

[8] 操林林,曹旭陽,王殿龍,等.半主動油氣懸架PID閉環(huán)控制研究[J].機械設(shè)計與制造,2014,12(12):9-12. CAO Linlin, CAO Xuyang, WANG Dianlong, et al. Research on PID Looped Control System of Semi-active Hydro-pneumatic Suspension[J].Machinery Design & Manufacture,2014,12(12):9-12.

[9] 張慶慶,潘公宇,陳立付. 可變剛度和阻尼的半主動懸架控制方法研究[J]. 拖拉機與農(nóng)用運輸車,2010,37( 3): 62-67 ZHANG Qingqing, PAN Gongyu, CHEN Lifu. Research on Control of Semi-active Suspension System with Variable Damping and Stiffness[J].Tractor & Farm Transporter,2010,37 (3):62-67.

[10] 高瑞貞,張京軍,趙子月,等.基于改進遺傳算法的半主動懸架系統(tǒng)模糊控制優(yōu)化研究[J].工程力學(xué),2012,29(1):240-248. GAO Ruizhen, ZHANG Jingjun, ZHAO Ziyue, et al.Research on the Fuzzy Control of Semi-active Suspension Systems Based on Improved Genetic Algorithms[J].Engineering Mechanics,2012,29(1):240-248.

[11] 王國敬.磁流變半主動懸架系統(tǒng)控制策略的研究[D].哈爾濱:哈爾濱工業(yè)大學(xué),2013. WANG Guojing. Research on Control Strategy of Magnetorheological Semi-active Suspension System [D]. Harbin :Harbin Institute of Technology,2013.

[12] 楊小衛(wèi). 磁流變減振器磁路分析及磁流變半主動懸架控制策略研究[D].鎮(zhèn)江:江蘇大學(xué),2007. YANG Xiaowei. Analyse of MR Damper’s Magnetic Circuit and Study on the Control Policy of MR Semi-active Suspension [D].Zhenjiang:Jiangsu University,2007.

[13] 朱一波.半主動懸架數(shù)學(xué)模型與自適應(yīng)神經(jīng)模糊控制研究[D].杭州:浙江工業(yè)大學(xué),2009. ZHU Yibo. Study on Mathematical Model and Adaptive Neuro-fuzzy Control of Semi-active Suspension [D].Hangzhou:Zhejiang University of Technology,2009.

[14] 蔣東升. 基于磁流變減振器的汽車半主動懸架的控制研究[D].錦州:遼寧工業(yè)大學(xué),2012. JIANG Dongsheng. Study on Control of Vehicle Semi-active Suspension Based on MR Damper[D].Jinzhou:Liaoning University of Technology,2012.

[15] 王祺明. 磁流變減振器半主動懸架控制仿真及試驗研究[D]. 合肥:合肥工業(yè)大學(xué),2009. WANG Qiming. The Study of SASS with MR Damper Control Simulation & Experiment[D]. Hefei :Hefei University of Technology,2009.

[16] 徐三玲.基于變論域模糊PID的數(shù)控機床伺服系統(tǒng)研究[D].哈爾濱:哈爾濱工業(yè)大學(xué),2013. XU Sanling. Research on Servo System of CNC Machine Tool Based on Variable Universe Fuzzy PID Control[D]. Harbin :Harbin Institute of Technology,2013.

[17] 崔家瑞,李擎,祝喬.飛航導(dǎo)彈變論域自適應(yīng)模糊PID控制方法[J].控制工程,2014,21(6):848-852. CUI Jiarui, LI Qing, ZHU Qiao. Variable Universe Adaptive Fuzzy PID Control Method in Maneuverable Missile[J].Control Engineering of China,2014,21(6):848-852.

[18] 魏慶福.飛機牽引車磁流變半主動懸架自適應(yīng)控制[D].哈爾濱:哈爾濱工業(yè)大學(xué),2012.WEI Qingfu. Adaptive Control for Aircraft Tractor’s Magnetorheological Semi-active Suspension[D]. Harbin :Harbin Institute of Technology,2012.

(編輯 陳 勇)

Control Method of Vehicle Semi Active Suspensions Based on Variable Universe Fuzzy Control

WANG Dayong WANG Hui

School of Mechnical Engineering,Liaoning Technical University,F(xiàn)uxin,Liaoning,123000

A fuzzy control method was proposed based on variable universe fuzzy PID for vehicle semi-active suspension systems herein. The goals were to improve the ride comfort of vehicles under random road excitations. The fuzzy PID controller was combined with variable universe method to solve the problems that the fuzzy rules were blind and the online regulation time was too long. The simulations and experimental results show that, variable universe fuzzy PID control semi-active suspension systems in the body vertical vibration velocity and acceleration are reduced by 46.56% and 29.21% respectively in comparison with conventional PID controls, compared to the passive suspension systems they are reduced by 58.05% and 49.74% respectively. The vehicle semi-active suspension fuzzy control method has better effects on vehicle ride comfort.

semi-active suspension; variable universe; fuzzy PID algorithm; vehicle vibration

2016-03-31

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

TP273

10.3969/j.issn.1004-132X.2017.03.019

王大勇,男,1979年生。遼寧工程技術(shù)大學(xué)機械工程學(xué)院副教授、博士研究生。主要研究方向為機械系統(tǒng)設(shè)計、建模及動態(tài)分析、仿真。E-mail:651725085@qq.com。王 慧,男,1960年生。遼寧工程技術(shù)大學(xué)機械工程學(xué)院副院長、教授、博士研究生導(dǎo)師。

猜你喜歡
振動模型
一半模型
振動的思考
噴水推進高速艇尾部振動響應(yīng)分析
重要模型『一線三等角』
This “Singing Highway”plays music
重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
振動攪拌 震動創(chuàng)新
中國公路(2017年18期)2018-01-23 03:00:38
中立型Emden-Fowler微分方程的振動性
3D打印中的模型分割與打包
FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
主站蜘蛛池模板: 911亚洲精品| 国产97区一区二区三区无码| 欧美综合成人| 国产区免费| 欧美激情视频一区| 欧美中文字幕一区| 国产日韩精品欧美一区灰| 国产精品亚洲精品爽爽| 国产精品v欧美| 国产第一色| 亚洲国产成人综合精品2020| 成人毛片在线播放| 日韩毛片免费| 99中文字幕亚洲一区二区| 一级毛片基地| 国产免费福利网站| 激情综合五月网| 免费毛片视频| 伊人久久青草青青综合| 亚洲欧美另类色图| 国产99在线| 9久久伊人精品综合| 亚洲国产精品日韩av专区| 亚洲午夜国产精品无卡| 久久不卡精品| 天堂成人在线视频| 亚洲福利片无码最新在线播放| 日韩精品久久久久久久电影蜜臀| 免费国产一级 片内射老| 亚瑟天堂久久一区二区影院| 朝桐光一区二区| 伊人激情综合网| 国产91在线|日本| 热re99久久精品国99热| 国产欧美日韩另类精彩视频| 免费看av在线网站网址| 精品国产一区二区三区在线观看| 国产在线自乱拍播放| 黄色网址免费在线| 日韩在线中文| 制服丝袜 91视频| 免费在线不卡视频| 不卡视频国产| 国产精品jizz在线观看软件| 精品国产免费观看| 欧美丝袜高跟鞋一区二区| 欧美成人一区午夜福利在线| 曰AV在线无码| 日韩精品成人在线| 一区二区三区高清视频国产女人| 成人综合网址| 久久99国产精品成人欧美| 青青青国产视频手机| 亚洲人成网站观看在线观看| 欧美激情网址| 亚洲欧美国产五月天综合| 国产一级毛片yw| 国产专区综合另类日韩一区 | 91麻豆久久久| 国产精女同一区二区三区久| 91精品视频在线播放| 91久久精品国产| 国产亚洲欧美日韩在线观看一区二区| 欧美亚洲另类在线观看| 亚洲综合久久成人AV| …亚洲 欧洲 另类 春色| 色丁丁毛片在线观看| 亚洲综合色吧| 国产精品嫩草影院视频| 国产精品一老牛影视频| 久久综合成人| 日韩在线观看网站| 免费毛片全部不收费的| 欧洲亚洲一区| 99精品伊人久久久大香线蕉| 国产亚洲男人的天堂在线观看| AV无码无在线观看免费| 国产在线高清一级毛片| 99视频在线精品免费观看6| 国产幂在线无码精品| 亚洲第一天堂无码专区| 秋霞国产在线|