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

    高鐵轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼比的協(xié)同優(yōu)化

    2016-11-17 03:17:44周長城于曰偉趙雷雷
    大連交通大學學報 2016年4期
    關鍵詞:阻尼比轉(zhuǎn)向架高鐵

    周長城,于曰偉,趙雷雷,2

    (1.山東理工大學 交通與車輛工程學院,山東 淄博 255049;2.北京郵電大學 自動化學院,北京 100876)*

    ?

    高鐵轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼比的協(xié)同優(yōu)化

    周長城1,于曰偉1,趙雷雷1,2

    (1.山東理工大學 交通與車輛工程學院,山東 淄博 255049;2.北京郵電大學 自動化學院,北京 100876)*

    根據(jù) 1/4車體 4自由度垂向振動模型,利用 MATLAB/Simulink,建立了高鐵轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼比協(xié)同優(yōu)化設計仿真模型;以人體振動舒適性最佳為目標,建立了阻尼比協(xié)同優(yōu)化設計數(shù)學模型.在此基礎上,以軌道高低不平順作為輸入激勵,以一系及二系懸掛垂向行程和一系懸掛垂向動作用力為約束條件,創(chuàng)建了轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼比協(xié)同優(yōu)化設計方法.通過實例對阻尼比進行了優(yōu)化設計,可知阻尼比的優(yōu)化設計值能夠改善高鐵的乘坐舒適性,表明所建立的高鐵轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼比的協(xié)同優(yōu)化設計方法是正確的,該研究為高鐵一系和二系垂向懸掛系統(tǒng)阻尼比的初始設計提供了重要的理論指導.

    高鐵 ;轉(zhuǎn)向架垂向懸掛 ;阻尼比 ;協(xié)同優(yōu)化

    0 引言

    高鐵必須具有良好的運行平穩(wěn)性和安全性,否則將會影響到乘客的乘坐舒適性,甚至會引起脫軌安全事故[1].一系和二系垂向懸掛系統(tǒng)作為轉(zhuǎn)向架懸掛系統(tǒng)的重要組成部分,其阻尼匹配對高鐵的運行平穩(wěn)性和安全性具有重要影響[2-3].然而,由于受高速列車懸掛最優(yōu)阻尼匹配理論的制約,目前國內(nèi)外對于高鐵轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼比的設計,大都是將一系和二系垂向懸掛系統(tǒng)分別單獨進行研究,并根據(jù)經(jīng)驗在可行性設計區(qū)間選擇某一一系和二系垂向懸掛系統(tǒng)阻尼比值,然后利用多體動力學軟件SIMPACK或ADAMS/Rail通過實體建模仿真確定其設計值[4-8].利用該方法所得到的一系和二系垂向懸掛系統(tǒng)阻尼比值,盡管可使車輛滿足當前行駛工況的要求,然而并非是轉(zhuǎn)向架垂向懸掛系統(tǒng)的最佳阻尼匹配值.隨著高鐵行駛速度的不斷提高,對車輛的乘坐舒適性和運行安全性提出了更高的設計要求,目前轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼匹配的設計方法不能給出具有指導意義的創(chuàng)新理論.

    近年來,已有許多國內(nèi)外學者對軌道車輛轉(zhuǎn)向架垂向懸掛系統(tǒng)進行了大量的研究,但這些研究主要是針對其減振器阻尼系數(shù)進行的,主要采用的方法有控制設計理論法、智能優(yōu)化設計法和建模仿真優(yōu)化法.例如,Sugahara Y,Takigami T等人[9-10],分別利用H∞和LOG控制算法對軌道車輛的一系垂向懸掛系統(tǒng)進行研究,給出了基于控制設計理論的一系垂向懸掛系統(tǒng)減振器的最佳阻尼系數(shù)優(yōu)化設計值;Nguyen H C,Zuo L等人[11-12]分別利用H∞和分散控制技術給出了基于控制設計理論的二系垂向懸掛系統(tǒng)減振器的最佳阻尼系數(shù)優(yōu)化設計值;Mei T X,Kim Y G,曾京等人[13-20],分別利用遺傳算法,神經(jīng)網(wǎng)絡算法,全局優(yōu)化算法,穩(wěn)健性設計,多目標優(yōu)化方法等給出了基于智能優(yōu)化方法的一系和二系垂向懸掛系統(tǒng)減振器的最佳阻尼系數(shù)優(yōu)化設計值;Zhang Y W,Nishimura K等人[21-22],分別利用ANSYS軟件和MATLAB軟件給出了基于建模仿真的一系和二系垂向懸掛系統(tǒng)懸掛參數(shù)的優(yōu)化設計值.雖然這些研究能夠給出一系或二系垂向懸掛系統(tǒng)減振器的最佳阻尼系數(shù)設計值,但這些研究所建立的振動模型或仿真模型中未考慮減振器端部連接結(jié)構(gòu)的彈性作用,且未曾給出具有指導意義的一系及二系垂向懸掛系統(tǒng)阻尼比的設計方法.

    本文根據(jù)1/4車體4自由度垂向振動模型,通過高鐵行駛平穩(wěn)性和安全性分析,對轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼比進行研究,并結(jié)合實例,對一系和二系垂向懸掛系統(tǒng)阻尼比進行優(yōu)化設計.

    1 1/4車體4自由度垂向振動模型

    1.1 模型的建立

    為了研究轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼對鐵道車輛頻率響應的影響,同時,又能比較直接地反映出減振器的特性,不受到整車復雜大系統(tǒng)的影響,本文采用1/4車體4自由度垂向振動模型,對高鐵一系及二系垂向懸掛系統(tǒng)阻尼比進行研究,如圖1所示.其中,模型考慮了一系及二系垂向減振器的端部連接剛度,坐標原點位于各自靜平衡位置處.

    圖1 1/4車體4自由度垂向振動模型

    圖中,m1為單個轉(zhuǎn)向架構(gòu)架質(zhì)量的一半,m2為單節(jié)車體滿載質(zhì)量的1/4;K1,K2分別為每臺轉(zhuǎn)向架單側(cè)一系和二系垂向懸掛彈簧的等效剛度;C1,C2分別為每臺轉(zhuǎn)向架單側(cè)一系和二系垂向減振器的等效阻尼系數(shù);Kd1,Kd2分別為每臺轉(zhuǎn)向架單側(cè)一系和二系垂向減振器的端部連接等效剛度;zd1,zd2分別為一系和二系垂向減振器的活塞桿垂向位移;z1,z2分別為轉(zhuǎn)向架構(gòu)架和車體的垂向位移;zv為軌道高低不平順隨機輸入.

    1.2 車體及轉(zhuǎn)向架構(gòu)架垂向振動微分方程

    根據(jù)1/4車體4自由度垂向振動模型,在不計軌道耦合振動作用的影響及減振器質(zhì)量情況下,利用牛頓第二定律,可建立車體及轉(zhuǎn)向架構(gòu)架垂向振動微分方程,即

    (1)

    利用上述振動微分方程,可對在軌道激勵下的高鐵行駛振動響應及一系和二系垂向懸掛系統(tǒng)最佳阻尼比的協(xié)同優(yōu)化設計進行研究.

    2 軌道高低不平順激擾模型

    2.1 軌道高低不平順

    軌道高低不平順是引起機車車輛產(chǎn)生垂向振動的主要原因.近年來,國內(nèi)、外對于軌道隨機不平順已進行了大量的研究,其中應用較為成熟、廣泛的主要有美國六級軌道譜和德國高速軌道譜[23].本文采用對高速軌道線路擬合程度較高的德國高低不平順作為高鐵垂向振動的軌道輸入,其中,高低不平順以空間頻率形式表示的解析表達式為

    (2)

    式中,Sv(Ω)為軌道高低不平順功率譜密度,Ω為軌道不平順空間頻率,Av為軌道粗糙度系數(shù);Ωc、Ωr為截斷空間頻率.其中,各已知參數(shù)值如表1所示.

    表1 德國高低不平順參數(shù)值

    其中,低干擾譜適合250 km/h及以上車速,高干擾譜適合250 km/h以下車速.

    2.2 軌道高低不平順時域樣本的模擬合成

    目前國內(nèi)外對于軌道不平順時域樣本的模擬合成,主要采用的方法有二次濾波法、三角級數(shù)法、白噪聲濾波法、以及基于功率譜密度采樣的軌道不平順數(shù)值模擬新方法等[23].本文采用對軌道線路擬合程度高且能夠真實反映線路實際情況的基于功率譜密度采樣的軌道不平順數(shù)值模擬新方法對軌道高低不平順時域樣本進行合成,其中,所合成的車輛運行速度為300 km/h時的軌道高低不平順時域模擬序列如圖2所示.

    圖2 軌道高低不平順時域模擬序列

    3 轉(zhuǎn)向架懸掛系統(tǒng)阻尼比協(xié)同優(yōu)化

    為使高鐵轉(zhuǎn)向架垂向懸掛系統(tǒng)的阻尼匹配達到最佳,本文利用MATLAB/Simulink和多島遺傳算法,通過構(gòu)建阻尼比協(xié)同優(yōu)化設計Simulink仿真模型,以人體振動舒適性最佳為目標,以軌道高低不平順作為輸入激勵,以一系及二系懸掛垂向行程和一系懸掛垂向動作用力為約束條件,對高鐵轉(zhuǎn)向架一系及二系垂向懸掛系統(tǒng)阻尼比進行協(xié)同優(yōu)化.

    3.1 舒適性評價指標

    為了對高鐵一系和二系垂向懸掛系統(tǒng)的阻尼比進行優(yōu)化設計,首先必須明確人體振動舒適性評價指標.目前,國內(nèi)外最常用的人體舒適性和健康評價指標為ISO2631標準[24]中提出的加權加速度均方根值,該標準認為人體對不同頻率及不同方向振動的敏感程度不同,由于本研究主要針對垂向振動對舒適性的影響進行評價,因此,僅計算垂向振動加權加速度均方根值.其中,在不同頻率下的振動頻率加權加速度均方根值的頻率加權值w(f)的計算公式為

    (3)

    3.2 阻尼比協(xié)同優(yōu)化設計目標函數(shù)

    根據(jù)每臺轉(zhuǎn)向架單側(cè)一系垂向減振器的等效阻尼系數(shù)C1、二系垂向減振器的阻尼系數(shù)C2,分別與各系懸掛系統(tǒng)參數(shù)及待優(yōu)化設計阻尼比之間的關系,可得

    (4)

    (5)

    式中,ξ1為一系垂向懸掛系統(tǒng)阻尼比,ξ2為二系垂向懸掛系統(tǒng)阻尼比.

    (6)

    3.3 阻尼比協(xié)同優(yōu)化設計約束條件

    為了保證高鐵高速行駛時具有良好的運行平穩(wěn)性和安全性,在對其轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼比進行優(yōu)化設計時,應滿足以下約束條件:

    (1)一系及二系懸掛垂向行程

    為了減小車輛高速行駛過程中撞擊限位的概率,使車輛具有良好的運行平穩(wěn)性和安全性,因此,一系及二系懸掛的垂向行程不應超出其垂向限位行程,即

    (7)

    (8)

    式中,z1-zv為一系懸掛垂向行程,z2-z1為二系懸掛垂向行程,[fd1]為一系懸掛垂向限位行程,[fd2]為二系懸掛垂向限位行程.

    (2)一系懸掛垂向動作用力

    為了使車輪不抬離軌道表面,以保障車輛的運行安全性,輪對所受的一系懸掛垂向動作用力不應超出其靜作用力,即

    (9)

    3.4 阻尼比協(xié)同優(yōu)化

    (1)阻尼比優(yōu)化設計仿真模型

    根據(jù)所建立的振動微分方程式(1)及式(4)、(5),利用MATLAB軟件的Simulink工具箱,構(gòu)建高鐵轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼比協(xié)同優(yōu)化設計Simulink仿真模型,如圖3所示.

    圖3 高鐵轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼比協(xié)同優(yōu)化設計Simulink仿真模型

    (2)一系及二系垂向懸掛系統(tǒng)阻尼比優(yōu)化設計

    多島遺傳算法作為一種偽并行遺傳算法可有效避免早熟和加快收斂速度,可以很好地在優(yōu)化域中尋找全局最優(yōu)解,為此,本文采用多島遺傳算法,利用MATLAB對高鐵轉(zhuǎn)向架一系和二系垂向懸掛系統(tǒng)阻尼比進行協(xié)同優(yōu)化,優(yōu)化設計流程圖如圖4所示,其中,參數(shù)設定如下:高鐵一系及二系垂向懸掛系統(tǒng)阻尼比的優(yōu)化設計初始值ξ1=0、ξ2=0,優(yōu)化范圍ξ1∈(0,0.5)、ξ2∈(0,0.5);多島遺傳算法的子群規(guī)模為10,島個數(shù)為10,進化代數(shù)為10,交叉概率為1,變異概率為0.01,遷移概率為0.01,遷移的間隔代數(shù)為5.

    圖4 優(yōu)化設計流程圖

    根據(jù)車輛參數(shù),所建立的Simulink仿真模型及式(6)~(9),以一系垂向懸掛系統(tǒng)阻尼比ξ1和二系垂向懸掛系統(tǒng)阻尼比ξ2為設計變量,以德國軌道高低不平順作為軌道輸入激勵,依據(jù)上述

    優(yōu)化設計流程,利用所編寫的優(yōu)化設計程序求目標函數(shù)Jo(ξ1,ξ2)的最小值,便可得到高鐵一系和二系垂向懸掛系統(tǒng)的最優(yōu)阻尼比ξ1、ξ2.

    4 設計實例

    5 仿真驗證

    在相同車輛結(jié)構(gòu)參數(shù)和軌道激勵下,分別對轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼比優(yōu)化設計前、后的該車輛1/4車體進行了模型仿真,其中,仿真所得到的該車輛優(yōu)化設計前、后車體垂向振動加速度的時域信號及功率譜密度對比曲線,分別如圖5、圖6所示.

    圖5 車體垂向振動加速度時域信號

    圖6 車體垂向振動加速度功率譜密度曲線

    分析圖5和圖6可知,高鐵一系和二系垂向懸掛系統(tǒng)阻尼比優(yōu)化設計后,其車體垂向振動加速度和功率譜密度值比優(yōu)化設計前有所降低.其中,優(yōu)化設計前的車體垂向振動加權加速度均方根值為0.24 m/s2,優(yōu)化設計后的為0.19 m/s2,舒適性提高了20.8%.對比可知,優(yōu)化設計后車輛的乘坐舒適性得到了明顯的提高,表明所設計的一系和二系垂向懸掛系統(tǒng)的最優(yōu)阻尼比值是可靠的.

    6 結(jié)論

    通過高鐵轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼比協(xié)同優(yōu)化方法的理論研究、實例設計及仿真驗證,可知:

    (1)根據(jù)1/4車體4自由度垂向振動模型,利用MATLAB/Simulink,可建立高鐵轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼比協(xié)同優(yōu)化設計仿真模型;

    (2)利用高鐵轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼比協(xié)同優(yōu)化設計仿真模型,以人體振動舒適性最佳為目標,以軌道高低不平順作為輸入激勵,以一系及二系懸掛垂向行程和一系懸掛垂向動作用力為約束條件,可建立轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼比協(xié)同優(yōu)化設計方法;

    (3)實例設計及對比驗證結(jié)果表明:所建立的高鐵轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼比協(xié)同優(yōu)化設計方法是可行的,能夠改善高鐵的乘坐舒適性,該研究為高鐵一系和二系垂向懸掛系統(tǒng)阻尼比的初始設計提供了重要的理論指導,為基于1/4車輛模型的高鐵轉(zhuǎn)向架垂向懸掛系統(tǒng)阻尼比的優(yōu)化設計提供了一種行之有效的指導方法.

    [1]王福天.車輛系統(tǒng)動力學[M].北京:中國鐵道出版社,1994.

    [2]AUCIELLO J, MELI E, FALOMI S, et al. Dynamic behaviour of tramways with different kinds of bogies[J]. Veh Syst Dyn, 2009, 47(7): 867-899.

    [3]孟宏,翟婉明,王開云.二系懸掛對機車動力學性能的影響[J].鐵道機車車輛,2005,25(5):1-4.

    [4]楊國楨,王福天.機車車輛液壓減振器[M].北京:中國鐵道出版社,2002.

    [5]VERROS G, GOUDAS H, NATSIAVAS S. Dynamics of large scale vehicle models using ADAMS/FLEX[C]. International ADAMS User Conference, 2000.

    [6]EOM B G, LEE H S. Assessment of running safety of railway vehicles using multibody dynamics[J]. Int. J. Precis. Eng. Man, 2010, 11(2): 315-320.

    [7]EICHBERGER A, HOFMANN G. TMPT: multi-body package SIMPACK[J]. Veh Syst Dyn, 2007, 45(Suppl): 207-216.

    [8]WALLRAPP O. Review of past developments in multi-body system dynamics at DLR-from FADYNA to SIMPACK[J]. Veh Syst Dyn, 2004, 41(5): 339-348.

    [9]SUGAHARA Y, TAKIGAMI T, KAZATO A. Suppression vertical vibration in railway vehicles through air spring damping control[J]. Journal of System Design and Dynamics, 2007, 1(2): 213-223.

    [10]SUGAHARA Y, TAKIGAMI T, KAZATO A, et al. Suppression vertical vibration in railway vehicles by damping force control of primary suspension using an LQG controller[J]. Journal of System Design and Dynamics, 2008, 2(1): 251-262.

    [11]NGUYEN H C, SONE A, IBA D, et al. Design of Passive Suspension System of Railway Vehicles via Control Theory[J]. Journal of System Design and Dynamics, 2008, 2(2): 518-527.

    [12]ZUO L, NAYFEH S A. Design of Passive Mechanical Systems via Decentralized Control Techniques[C]. 43rd AIAA/ASME/ ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, 2002, AIAA2002-1282: 1-9.

    [13]MEI T X, GOODALL R M. Use of multiobjective genetic algorithms to optimize inter-vehicle active suspensions[J]. Journal of Rail and Rapid Transit , 2002, 216(1): 53-63.

    [14]MASTINU R M, GOBBI M. On the optimal design of railway passenger vehicles[J]. Journal of Rail and Rapid Transit, 2001, 215: 111-124.

    [15]KIM Y G, PARK C K, PARK T W. Design optimization for suspension system of high speed train using neural network[J]. JSME International Journal Series C, 2003, 46(2): 727-735.

    [16]KUZNETSOV A, MAMMADOV M, SULTAN I, et al. Optimization of a quarter-car suspension model coupled with the driver biomechanical effects[J]. Journal of Sound and Vibration, 2011, 330(12): 2937-2946.

    [17]BATOU A, SOIZE C, CHOI C K, et al. Robust design in multibody dynamics-application to vehicle ride-comfort optimization[J]. Procedia IUTAM, 2015, 13: 90-97.

    [18]GEORGIOU G, VERROS G, NATSIAVAS S. Multi-objective optimization of quarter-car models with a passive or semi-active suspension System[J]. Veh Syst Dyn, 2007, 45 (1): 77-92.

    [19]郝建華,曾京,鄔平波.鐵道客車垂向隨機減振及懸掛參數(shù)優(yōu)化[J].鐵道學報,2006,28(6):35-40.

    [20]SUGAHARA Y, KAZATO A, TAKIGAMI T, et al. Suppression of vertical vibration in railway vehicles by controlling the damping force of primary and secondary suspensions[J]. QR of RTRI, 2008, 49(1): 7-15.

    [21]ZHANG Y W, ZHAO Y, ZHANG Y H, et al. Riding comfort optimization of railway trains based on pseudo-excitation method and symplectic method[J]. Journal of Sound and Vibration, 2013, 332(21): 5255-5270.

    [22]NISHIMURA K, PERKINS N C, ZHANG W M. Suspension dynamics and design optimization of a high speed railway vehicle[C]. In Proceedings of the 2004 ASME/IEEE Joint Rail Conference, Maryland, USA, 2004.

    [23]翟婉明.車輛-軌道耦合動力學[M].4版,北京:科學出版社,2015.

    [24]ISO 2631-1 . Mechanical Vibration and Shock-Evaluation of Human Exposure to Whole-driver Vibration-Part1: General Requirements, International Organization for Standardization[P].Geneva [s.n.], 1997.

    Collaborative Optimization of Bogie Vertical Suspension Damping Ratio for High-Speed Rail

    ZHOU Changcheng1, YU Yuewei1, ZHAO Leilei1,2

    (1. School of Transportation and Vehicle Engineering, Shandong University of Technology, Zibo 255049, China; 2. School of Automation, Beijing University of Posts and Telecommunications, Beijing 100876, China)

    According to the 1/4 vehicle body four-degree-of-freedom vertical vibration model of high-speed rail, a collaborative optimal design simulation model for damping ratio of bogie vertical suspension was established using MATLAB/Simulink. Taking optimal ride comfort as target, a collaborative optimal design mathematical model for damping ratio of bogie vertical suspension was built. Using track vertical profile irregularity as input and using the primary and secondary suspension vertical stroke and the vertical dynamic force of primary suspension as constraint conditions, a collaborative optimization design method for damping ratio of the bogie vertical suspension was presented. With a practical example of high-speed rail, the damping ratio was designed. The results show that the damping ratio value designed can significantly improve ride comfort. Thus, the collaborative optimization method of bogie vertical suspension damping ratio for high-speed rail is correct. This research provides important theoretical basis for the initial design of the damping ratio of the primary and secondary vertical suspension system of high-speed rail.

    high-speed rail; bogie vertical suspension; damping ratio; collaborative optimization

    1673-9590(2016)04-0055-06

    2015-12-14

    國家自然科學基金資助項目(51575325);山東省自然科學基金資助項目(ZR2013EEM007);山東省重點研發(fā)計劃資助項目(2015GGX105006)

    周長城(1962-),男,教授,博士,主要從事汽車和軌道車輛懸架設計及理論的研究E-mail:greatwall@sdut.edu.cn.

    A

    猜你喜歡
    阻尼比轉(zhuǎn)向架高鐵
    一起坐高鐵
    第一次乘高鐵
    基于細觀結(jié)構(gòu)的原狀黃土動彈性模量和阻尼比試驗研究
    地震研究(2021年1期)2021-04-13 01:05:24
    20t軸重米軌轉(zhuǎn)向架裝用搖枕的研制
    基于SPS模式的轉(zhuǎn)向架軸箱裝配線仿真研究
    高鐵會飛嗎
    黏滯阻尼器在時程分析下的附加有效阻尼比研究
    振動與沖擊(2019年4期)2019-02-22 02:33:34
    基于結(jié)構(gòu)應力的轉(zhuǎn)向架構(gòu)架焊縫疲勞強度研究
    波形分析法求解公路橋梁阻尼比的探討
    上海公路(2018年3期)2018-03-21 05:55:48
    人地百米建高鐵
    學與玩(2017年12期)2017-02-16 06:51:16
    久久久午夜欧美精品| 男插女下体视频免费在线播放| 欧美性感艳星| 性色avwww在线观看| 在线观看午夜福利视频| 非洲黑人性xxxx精品又粗又长| 禁无遮挡网站| aaaaa片日本免费| 国产精品久久久久久久久免| 亚洲三级黄色毛片| 女人十人毛片免费观看3o分钟| 一级毛片我不卡| 国产精品一区www在线观看| 成年免费大片在线观看| 乱码一卡2卡4卡精品| 欧美潮喷喷水| 亚洲人成网站在线播放欧美日韩| 午夜福利在线在线| 我的女老师完整版在线观看| 亚洲内射少妇av| 99九九线精品视频在线观看视频| 亚洲最大成人手机在线| 国产日本99.免费观看| 中文字幕免费在线视频6| 又黄又爽又刺激的免费视频.| 亚洲丝袜综合中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 国产免费一级a男人的天堂| 国产精品久久久久久亚洲av鲁大| 亚洲av.av天堂| 亚洲国产欧洲综合997久久,| 麻豆国产97在线/欧美| av在线天堂中文字幕| 天天躁日日操中文字幕| 午夜激情欧美在线| 精品人妻熟女av久视频| 亚洲久久久久久中文字幕| 淫妇啪啪啪对白视频| 精品久久久噜噜| 欧美日韩国产亚洲二区| 欧美最黄视频在线播放免费| 精品少妇黑人巨大在线播放 | 村上凉子中文字幕在线| 日韩欧美三级三区| а√天堂www在线а√下载| 国产精华一区二区三区| 久久国内精品自在自线图片| 午夜精品国产一区二区电影 | 久久精品91蜜桃| 欧美+亚洲+日韩+国产| 少妇丰满av| 亚洲av成人av| 亚洲婷婷狠狠爱综合网| 成人鲁丝片一二三区免费| 老师上课跳d突然被开到最大视频| 成人精品一区二区免费| 亚洲人成网站在线播放欧美日韩| 少妇丰满av| 女生性感内裤真人,穿戴方法视频| 亚洲无线在线观看| 午夜免费激情av| 少妇熟女欧美另类| 真实男女啪啪啪动态图| 国内揄拍国产精品人妻在线| 国产高清有码在线观看视频| 亚洲不卡免费看| 中国美白少妇内射xxxbb| 精品不卡国产一区二区三区| 一级a爱片免费观看的视频| 亚洲精品色激情综合| 男女做爰动态图高潮gif福利片| 长腿黑丝高跟| 中国美女看黄片| 中文字幕精品亚洲无线码一区| 在线播放国产精品三级| 一本精品99久久精品77| 黄色一级大片看看| 国产三级在线视频| 日韩三级伦理在线观看| 色综合站精品国产| 麻豆成人午夜福利视频| 久久久久久久久久久丰满| 国产精品野战在线观看| 亚洲av五月六月丁香网| 村上凉子中文字幕在线| 淫秽高清视频在线观看| 黄色视频,在线免费观看| 黄色视频,在线免费观看| 日韩av在线大香蕉| 久久久国产成人精品二区| 亚洲av成人精品一区久久| 国产 一区 欧美 日韩| 国产精品不卡视频一区二区| 国产高清视频在线播放一区| 长腿黑丝高跟| 免费大片18禁| 午夜福利在线观看免费完整高清在 | 亚洲精品国产av成人精品 | 少妇的逼好多水| 禁无遮挡网站| 国产伦精品一区二区三区四那| 99久久久亚洲精品蜜臀av| 久久精品影院6| 午夜爱爱视频在线播放| av专区在线播放| 成人av一区二区三区在线看| 丰满乱子伦码专区| 97超级碰碰碰精品色视频在线观看| 国产精品人妻久久久久久| 亚洲精品亚洲一区二区| 免费观看人在逋| 精品午夜福利视频在线观看一区| 高清午夜精品一区二区三区 | 99久国产av精品国产电影| 日本一本二区三区精品| 免费在线观看影片大全网站| a级毛片a级免费在线| 非洲黑人性xxxx精品又粗又长| 在线播放无遮挡| 搡老岳熟女国产| 91久久精品国产一区二区成人| 99热这里只有精品一区| 一级黄色大片毛片| 男女边吃奶边做爰视频| 亚洲一区高清亚洲精品| 黄色配什么色好看| 久久精品国产清高在天天线| 成人av在线播放网站| 波多野结衣巨乳人妻| 久久精品人妻少妇| 亚洲精品国产成人久久av| 色播亚洲综合网| 人妻久久中文字幕网| 婷婷亚洲欧美| 香蕉av资源在线| 久久久久久国产a免费观看| 久久久国产成人免费| 性欧美人与动物交配| 成人精品一区二区免费| 成人三级黄色视频| 在线观看午夜福利视频| 精品久久久久久久人妻蜜臀av| 久久久午夜欧美精品| 日韩欧美免费精品| 在线播放国产精品三级| 婷婷精品国产亚洲av在线| 国产精品精品国产色婷婷| 啦啦啦啦在线视频资源| av天堂中文字幕网| 亚洲av成人av| 国产v大片淫在线免费观看| 大香蕉久久网| 国产精品一区二区三区四区久久| 女人被狂操c到高潮| 国产毛片a区久久久久| 久久亚洲精品不卡| 乱系列少妇在线播放| 看黄色毛片网站| 欧美又色又爽又黄视频| 男女视频在线观看网站免费| 精品久久久噜噜| 国产精品综合久久久久久久免费| 一个人免费在线观看电影| 内地一区二区视频在线| 国产精品亚洲一级av第二区| 国产精品爽爽va在线观看网站| 成年女人毛片免费观看观看9| 国产黄片美女视频| 国产精品国产三级国产av玫瑰| 久久午夜亚洲精品久久| 欧美日韩国产亚洲二区| 国产一区亚洲一区在线观看| 久99久视频精品免费| 午夜日韩欧美国产| 欧美日韩精品成人综合77777| 性插视频无遮挡在线免费观看| 国产大屁股一区二区在线视频| 亚洲一区二区三区色噜噜| 久久久久久大精品| 久久精品国产亚洲av天美| 在线天堂最新版资源| 女生性感内裤真人,穿戴方法视频| 久久久国产成人精品二区| 亚洲欧美日韩东京热| 国内揄拍国产精品人妻在线| 久久久久九九精品影院| 国产精品不卡视频一区二区| 国产精品免费一区二区三区在线| 亚洲国产高清在线一区二区三| 国产黄a三级三级三级人| 欧美色欧美亚洲另类二区| 少妇熟女aⅴ在线视频| 久久婷婷人人爽人人干人人爱| 好男人在线观看高清免费视频| 人人妻人人看人人澡| 91av网一区二区| 亚洲高清免费不卡视频| 中文资源天堂在线| 久久亚洲精品不卡| 欧美色欧美亚洲另类二区| 国产精品1区2区在线观看.| 国产成年人精品一区二区| 熟女人妻精品中文字幕| 国产高清不卡午夜福利| 最近视频中文字幕2019在线8| 人妻久久中文字幕网| 亚洲欧美精品自产自拍| 精品久久国产蜜桃| 亚洲av不卡在线观看| 成年女人看的毛片在线观看| 成人漫画全彩无遮挡| 国产69精品久久久久777片| 国产白丝娇喘喷水9色精品| 深夜精品福利| 欧美3d第一页| 国产视频一区二区在线看| 色噜噜av男人的天堂激情| 热99在线观看视频| 少妇熟女aⅴ在线视频| 成人毛片a级毛片在线播放| 色吧在线观看| 精品欧美国产一区二区三| 精品久久久噜噜| 嫩草影视91久久| 日韩欧美三级三区| 欧美zozozo另类| 麻豆av噜噜一区二区三区| av专区在线播放| ponron亚洲| 国产国拍精品亚洲av在线观看| 久久午夜亚洲精品久久| 99热这里只有是精品50| 九九在线视频观看精品| 日韩制服骚丝袜av| 69人妻影院| 免费在线观看成人毛片| 免费人成在线观看视频色| 我的老师免费观看完整版| 深夜精品福利| 好男人在线观看高清免费视频| 国内精品久久久久精免费| 欧美日韩一区二区视频在线观看视频在线 | 久久午夜福利片| 国产av一区在线观看免费| 大香蕉久久网| 日产精品乱码卡一卡2卡三| 在线免费观看的www视频| 在线观看午夜福利视频| 老司机福利观看| 欧美+日韩+精品| 精品人妻偷拍中文字幕| 九色成人免费人妻av| 欧美色视频一区免费| 国产色婷婷99| 久久精品91蜜桃| 亚洲国产欧洲综合997久久,| 简卡轻食公司| 看片在线看免费视频| 国产精品亚洲一级av第二区| 久久亚洲精品不卡| 国产成人a∨麻豆精品| 久久精品国产亚洲av天美| 黄色日韩在线| 久久99热这里只有精品18| 亚洲va在线va天堂va国产| 久久精品夜色国产| 人妻丰满熟妇av一区二区三区| 日本黄大片高清| 有码 亚洲区| 三级毛片av免费| 日韩欧美在线乱码| 亚洲真实伦在线观看| 超碰av人人做人人爽久久| 一进一出抽搐动态| av在线亚洲专区| 又黄又爽又免费观看的视频| 日韩欧美在线乱码| 18禁黄网站禁片免费观看直播| 我要搜黄色片| 国产69精品久久久久777片| 亚洲精品一卡2卡三卡4卡5卡| 搡女人真爽免费视频火全软件 | 成人无遮挡网站| 性色avwww在线观看| 高清毛片免费看| 亚洲av五月六月丁香网| 国产真实伦视频高清在线观看| 男女边吃奶边做爰视频| 久久人人爽人人片av| 亚洲中文日韩欧美视频| 午夜影院日韩av| 日日摸夜夜添夜夜爱| 国产在视频线在精品| 熟女电影av网| 免费高清视频大片| 在线免费观看不下载黄p国产| 中国美白少妇内射xxxbb| 日韩中字成人| 亚洲高清免费不卡视频| 国产高清三级在线| 在线观看av片永久免费下载| 偷拍熟女少妇极品色| 中文字幕免费在线视频6| 国产高清不卡午夜福利| 国产国拍精品亚洲av在线观看| 少妇熟女欧美另类| 全区人妻精品视频| 亚洲精品影视一区二区三区av| 国产亚洲精品av在线| 五月玫瑰六月丁香| 丝袜喷水一区| 男人舔女人下体高潮全视频| 欧美3d第一页| 又粗又爽又猛毛片免费看| 俺也久久电影网| 日韩,欧美,国产一区二区三区 | 久久精品91蜜桃| 日日摸夜夜添夜夜爱| 日韩成人伦理影院| 国产精品久久久久久精品电影| 亚洲精品在线观看二区| 老熟妇乱子伦视频在线观看| eeuss影院久久| 久久韩国三级中文字幕| 国内久久婷婷六月综合欲色啪| 你懂的网址亚洲精品在线观看 | a级一级毛片免费在线观看| 久久人人精品亚洲av| 欧美zozozo另类| 深爱激情五月婷婷| 国内揄拍国产精品人妻在线| 一级毛片我不卡| 看黄色毛片网站| 国内揄拍国产精品人妻在线| 在线免费观看不下载黄p国产| 白带黄色成豆腐渣| 少妇被粗大猛烈的视频| 夜夜爽天天搞| 国产又黄又爽又无遮挡在线| 精品熟女少妇av免费看| 色综合亚洲欧美另类图片| 国产精品一区www在线观看| 少妇的逼水好多| 日韩中字成人| 最近手机中文字幕大全| 淫妇啪啪啪对白视频| 免费观看在线日韩| 狂野欧美激情性xxxx在线观看| 男人狂女人下面高潮的视频| 老熟妇乱子伦视频在线观看| 精品午夜福利在线看| 最近的中文字幕免费完整| 国产精品乱码一区二三区的特点| 久久久久久伊人网av| 午夜老司机福利剧场| 久久久久久久久久久丰满| 麻豆久久精品国产亚洲av| 久久久久久大精品| 欧美成人精品欧美一级黄| 精品不卡国产一区二区三区| 日韩三级伦理在线观看| 日本撒尿小便嘘嘘汇集6| 一个人看的www免费观看视频| 精品乱码久久久久久99久播| av卡一久久| 国产欧美日韩精品一区二区| 黄片wwwwww| 国产精品一区www在线观看| 国内精品一区二区在线观看| 国产精品一区二区三区四区久久| 成人鲁丝片一二三区免费| 国产爱豆传媒在线观看| 六月丁香七月| 亚洲av中文字字幕乱码综合| 香蕉av资源在线| 99热全是精品| 国产白丝娇喘喷水9色精品| 欧美性猛交黑人性爽| 久久精品国产亚洲av天美| av在线观看视频网站免费| 我要搜黄色片| 国产91av在线免费观看| 97碰自拍视频| 久久天躁狠狠躁夜夜2o2o| 91在线观看av| 日韩三级伦理在线观看| 成人亚洲精品av一区二区| 国产高清视频在线播放一区| 欧美zozozo另类| 可以在线观看毛片的网站| 亚洲精品日韩av片在线观看| 日韩高清综合在线| 国产亚洲精品综合一区在线观看| www.色视频.com| 一级黄片播放器| h日本视频在线播放| 国产精品日韩av在线免费观看| 人妻久久中文字幕网| 国产高清视频在线播放一区| 卡戴珊不雅视频在线播放| 国产探花在线观看一区二区| 一a级毛片在线观看| 麻豆国产av国片精品| 日韩精品中文字幕看吧| 国产亚洲91精品色在线| 成人鲁丝片一二三区免费| 成人特级黄色片久久久久久久| 色噜噜av男人的天堂激情| 日本爱情动作片www.在线观看 | 国产精品亚洲一级av第二区| 久久中文看片网| 色播亚洲综合网| 日韩欧美在线乱码| 亚洲欧美精品自产自拍| 免费高清视频大片| 91精品国产九色| 精品一区二区三区人妻视频| 噜噜噜噜噜久久久久久91| 成年av动漫网址| 国产女主播在线喷水免费视频网站 | 欧美日韩精品成人综合77777| 亚洲最大成人手机在线| 日本熟妇午夜| 我的女老师完整版在线观看| 夜夜夜夜夜久久久久| 亚洲中文字幕日韩| 麻豆成人午夜福利视频| 午夜福利高清视频| 乱系列少妇在线播放| 久久热精品热| 久久99热6这里只有精品| 亚洲成人精品中文字幕电影| 免费看av在线观看网站| 久久精品久久久久久噜噜老黄 | 最近在线观看免费完整版| 人妻夜夜爽99麻豆av| 午夜爱爱视频在线播放| 日产精品乱码卡一卡2卡三| 在线看三级毛片| 日本五十路高清| 男人狂女人下面高潮的视频| 成年女人看的毛片在线观看| 一级毛片久久久久久久久女| 又爽又黄a免费视频| 99国产极品粉嫩在线观看| 亚洲最大成人av| 桃色一区二区三区在线观看| 欧美三级亚洲精品| 噜噜噜噜噜久久久久久91| 村上凉子中文字幕在线| 91久久精品国产一区二区成人| 一级毛片aaaaaa免费看小| 亚洲高清免费不卡视频| 日韩一区二区视频免费看| 亚洲综合色惰| 欧美3d第一页| 亚洲自拍偷在线| 欧美色视频一区免费| 中文亚洲av片在线观看爽| 高清毛片免费观看视频网站| 插逼视频在线观看| av卡一久久| 久久久久国产精品人妻aⅴ院| 噜噜噜噜噜久久久久久91| 午夜精品一区二区三区免费看| 亚洲一级一片aⅴ在线观看| 日本在线视频免费播放| 男人舔女人下体高潮全视频| 国产视频内射| 国产精品亚洲一级av第二区| 少妇猛男粗大的猛烈进出视频 | 男人和女人高潮做爰伦理| 欧美区成人在线视频| 女人十人毛片免费观看3o分钟| 日韩在线高清观看一区二区三区| 波多野结衣高清作品| 国内久久婷婷六月综合欲色啪| 少妇被粗大猛烈的视频| 99热这里只有是精品在线观看| 嫩草影视91久久| 国产真实乱freesex| 高清午夜精品一区二区三区 | 男人的好看免费观看在线视频| 老熟妇乱子伦视频在线观看| 午夜亚洲福利在线播放| 国产精品亚洲美女久久久| 日本a在线网址| 一个人看的www免费观看视频| 日韩亚洲欧美综合| 亚洲精品一区av在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲熟妇中文字幕五十中出| 欧美人与善性xxx| 性插视频无遮挡在线免费观看| 亚洲一区二区三区色噜噜| 少妇熟女欧美另类| 男女做爰动态图高潮gif福利片| 天堂网av新在线| 噜噜噜噜噜久久久久久91| 又黄又爽又免费观看的视频| 亚洲欧美日韩高清专用| 99国产极品粉嫩在线观看| 高清毛片免费看| 亚洲熟妇中文字幕五十中出| 搞女人的毛片| 97热精品久久久久久| 一区二区三区高清视频在线| 成年免费大片在线观看| 99久久中文字幕三级久久日本| 国产欧美日韩精品一区二区| 岛国在线免费视频观看| 国产一级毛片七仙女欲春2| 最新在线观看一区二区三区| 美女黄网站色视频| 欧美高清成人免费视频www| 毛片一级片免费看久久久久| 不卡一级毛片| 亚洲成a人片在线一区二区| 在线观看一区二区三区| 国产欧美日韩精品一区二区| 亚洲aⅴ乱码一区二区在线播放| 插逼视频在线观看| 久久99热6这里只有精品| 一a级毛片在线观看| 天堂网av新在线| av在线天堂中文字幕| 国产精品1区2区在线观看.| 亚洲国产色片| 婷婷色综合大香蕉| 久久人妻av系列| 中文字幕熟女人妻在线| 在线观看免费视频日本深夜| 国产成人aa在线观看| 国产一区二区在线观看日韩| 国产高清三级在线| 国产黄色小视频在线观看| 全区人妻精品视频| 日韩亚洲欧美综合| av天堂在线播放| 久久精品国产亚洲av天美| 国产一区二区在线av高清观看| 伦精品一区二区三区| 亚洲av一区综合| 日韩精品中文字幕看吧| 国内精品宾馆在线| 欧美日韩精品成人综合77777| 97超级碰碰碰精品色视频在线观看| 成人国产麻豆网| 欧美国产日韩亚洲一区| 久久精品国产亚洲av香蕉五月| 毛片女人毛片| 亚洲欧美日韩无卡精品| 老师上课跳d突然被开到最大视频| or卡值多少钱| 欧美日本视频| 插阴视频在线观看视频| 波多野结衣高清无吗| 一级av片app| 亚洲18禁久久av| 欧美高清性xxxxhd video| 日韩精品有码人妻一区| 国产伦精品一区二区三区四那| 九九热线精品视视频播放| 亚洲色图av天堂| 亚洲av熟女| 91久久精品电影网| 亚洲精品乱码久久久v下载方式| 国产成人91sexporn| 色5月婷婷丁香| 免费av毛片视频| 国产伦一二天堂av在线观看| 日本免费a在线| 亚洲第一区二区三区不卡| 丰满人妻一区二区三区视频av| 国产激情偷乱视频一区二区| 激情 狠狠 欧美| 国产av麻豆久久久久久久| 日本爱情动作片www.在线观看 | 我要搜黄色片| 寂寞人妻少妇视频99o| 精品午夜福利视频在线观看一区| 黄片wwwwww| 中文字幕av在线有码专区| 成人无遮挡网站| 亚洲色图av天堂| 在线免费观看不下载黄p国产| 色吧在线观看| 亚洲在线观看片| av视频在线观看入口| 午夜福利高清视频| 精品无人区乱码1区二区| 91麻豆精品激情在线观看国产| 久久精品国产亚洲av涩爱 | 亚洲欧美清纯卡通| 国产毛片a区久久久久| 国产麻豆成人av免费视频| 国产精品国产三级国产av玫瑰| 天堂网av新在线| av卡一久久| 亚洲真实伦在线观看| 国产单亲对白刺激| 国产精品一区二区三区四区久久| 精品福利观看| videossex国产| 搡老岳熟女国产| 看黄色毛片网站| 一夜夜www| 久久欧美精品欧美久久欧美| 国产 一区 欧美 日韩| 亚洲av第一区精品v没综合| 亚洲精品一卡2卡三卡4卡5卡| 在线看三级毛片| 久久人人爽人人爽人人片va| 国内精品一区二区在线观看| 亚洲精品影视一区二区三区av|