• <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
    嫩草影院入口| 啦啦啦啦在线视频资源| 精品久久久噜噜| 熟妇人妻久久中文字幕3abv| 亚洲av中文av极速乱| 国产成人精品婷婷| 久久精品久久精品一区二区三区| 午夜激情欧美在线| 联通29元200g的流量卡| ponron亚洲| 中文欧美无线码| 最近最新中文字幕免费大全7| 国产免费视频播放在线视频 | 中国美白少妇内射xxxbb| 国产高潮美女av| 国产成人91sexporn| 一个人免费在线观看电影| 一级黄片播放器| 精品国内亚洲2022精品成人| ponron亚洲| 一边亲一边摸免费视频| 亚洲精品日本国产第一区| 乱码一卡2卡4卡精品| 欧美激情国产日韩精品一区| 国产精品一区二区三区四区免费观看| 精品不卡国产一区二区三区| 亚洲在久久综合| av在线亚洲专区| 欧美激情久久久久久爽电影| 国产亚洲精品av在线| 欧美人与善性xxx| 国产免费又黄又爽又色| 天堂av国产一区二区熟女人妻| 午夜免费男女啪啪视频观看| 欧美日本视频| 人人妻人人看人人澡| 十八禁网站网址无遮挡 | 街头女战士在线观看网站| 夫妻午夜视频| 日本黄色片子视频| www.av在线官网国产| 国产 一区精品| 一级二级三级毛片免费看| 非洲黑人性xxxx精品又粗又长| 美女被艹到高潮喷水动态| 亚洲国产高清在线一区二区三| 精品久久久久久电影网| 精品欧美国产一区二区三| 老女人水多毛片| 全区人妻精品视频| 亚洲18禁久久av| 精品久久久精品久久久| 伊人久久国产一区二区| 成人亚洲精品一区在线观看 | av福利片在线观看| av又黄又爽大尺度在线免费看| 午夜福利视频1000在线观看| 国产欧美日韩精品一区二区| 婷婷色麻豆天堂久久| 91久久精品国产一区二区三区| 国产精品美女特级片免费视频播放器| 国产伦理片在线播放av一区| 国产黄频视频在线观看| 最近中文字幕高清免费大全6| 一级av片app| 久久久久久久久久成人| 免费播放大片免费观看视频在线观看| 最后的刺客免费高清国语| 久久久久久久久久久丰满| 成人亚洲精品一区在线观看 | 中文字幕av成人在线电影| 亚洲国产欧美在线一区| 免费黄色在线免费观看| 特级一级黄色大片| 亚洲国产av新网站| 高清欧美精品videossex| 天堂网av新在线| 国产女主播在线喷水免费视频网站 | 国产一区二区亚洲精品在线观看| 美女高潮的动态| 一个人看视频在线观看www免费| 国产真实伦视频高清在线观看| 国产精品女同一区二区软件| 伊人久久国产一区二区| 精品人妻偷拍中文字幕| 搡老妇女老女人老熟妇| 精品亚洲乱码少妇综合久久| 边亲边吃奶的免费视频| 超碰av人人做人人爽久久| 亚洲一区高清亚洲精品| 三级经典国产精品| 亚洲av福利一区| 久久久久精品性色| 91狼人影院| 精华霜和精华液先用哪个| 丝袜美腿在线中文| 乱码一卡2卡4卡精品| 寂寞人妻少妇视频99o| 精品久久国产蜜桃| 久久草成人影院| 国产精品.久久久| 日韩制服骚丝袜av| 成人av在线播放网站| 中文字幕av在线有码专区| 久久久午夜欧美精品| 超碰av人人做人人爽久久| 午夜福利成人在线免费观看| 国产午夜福利久久久久久| 亚洲av福利一区| 午夜激情欧美在线| 女人久久www免费人成看片| 日本午夜av视频| 综合色av麻豆| 卡戴珊不雅视频在线播放| 久久久久性生活片| 日本熟妇午夜| 在线免费观看的www视频| 久久久久网色| 人人妻人人看人人澡| 精品人妻视频免费看| 国精品久久久久久国模美| 亚洲av在线观看美女高潮| 插逼视频在线观看| 男女国产视频网站| 特大巨黑吊av在线直播| av.在线天堂| 精品久久国产蜜桃| 中文字幕免费在线视频6| 老司机影院毛片| 少妇丰满av| 亚洲综合色惰| 高清午夜精品一区二区三区| 夜夜看夜夜爽夜夜摸| 22中文网久久字幕| 人人妻人人澡欧美一区二区| 日韩欧美三级三区| 久久韩国三级中文字幕| 中文字幕久久专区| 人人妻人人澡欧美一区二区| 18禁裸乳无遮挡免费网站照片| 精品一区二区三卡| 亚洲精品日韩av片在线观看| 一区二区三区免费毛片| 80岁老熟妇乱子伦牲交| 热99在线观看视频| 一级黄片播放器| 久久久久国产网址| 久久久久久久久久黄片| 成人鲁丝片一二三区免费| 黄色日韩在线| 免费播放大片免费观看视频在线观看| 亚洲伊人久久精品综合| 成人欧美大片| av黄色大香蕉| 少妇的逼水好多| 欧美精品一区二区大全| 成人漫画全彩无遮挡| 亚洲av不卡在线观看| 人体艺术视频欧美日本| 少妇裸体淫交视频免费看高清| 欧美高清性xxxxhd video| 寂寞人妻少妇视频99o| 中国国产av一级| 午夜福利成人在线免费观看| 国产爱豆传媒在线观看| 亚洲国产精品成人综合色| 天堂网av新在线| 国产三级在线视频| 国产老妇伦熟女老妇高清| 97超碰精品成人国产| 日韩伦理黄色片| 一本一本综合久久| 日韩,欧美,国产一区二区三区| 亚洲综合色惰| 婷婷色av中文字幕| 22中文网久久字幕| 午夜老司机福利剧场| 一级毛片aaaaaa免费看小| 黄片无遮挡物在线观看| 成人av在线播放网站| 日韩,欧美,国产一区二区三区| 国产精品国产三级国产av玫瑰| 国产黄色视频一区二区在线观看| 三级男女做爰猛烈吃奶摸视频| 国国产精品蜜臀av免费| 国产在线一区二区三区精| 纵有疾风起免费观看全集完整版 | 99热6这里只有精品| www.av在线官网国产| 欧美精品国产亚洲| 亚洲av电影不卡..在线观看| 精品人妻偷拍中文字幕| 一级毛片我不卡| 亚洲人成网站在线观看播放| 亚洲av免费在线观看| 美女高潮的动态| 夜夜看夜夜爽夜夜摸| 一级二级三级毛片免费看| 中文字幕人妻熟人妻熟丝袜美| 国产爱豆传媒在线观看| 一本久久精品| 午夜免费观看性视频| 精品一区二区三区人妻视频| 毛片一级片免费看久久久久| 国产色爽女视频免费观看| 亚洲电影在线观看av| 青青草视频在线视频观看| 久久精品国产鲁丝片午夜精品| 国产精品一区二区性色av| 91在线精品国自产拍蜜月| 卡戴珊不雅视频在线播放| 久久精品国产鲁丝片午夜精品| 久久久亚洲精品成人影院| 欧美性猛交╳xxx乱大交人| 又爽又黄a免费视频| 久久久午夜欧美精品| videossex国产| 18+在线观看网站| 3wmmmm亚洲av在线观看| 精品欧美国产一区二区三| 久久鲁丝午夜福利片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲人与动物交配视频| 人体艺术视频欧美日本| 丝袜美腿在线中文| 久久久久久久久久成人| 久久久久久国产a免费观看| 九九爱精品视频在线观看| 22中文网久久字幕| 2021少妇久久久久久久久久久| 欧美人与善性xxx| 日韩三级伦理在线观看| 亚洲不卡免费看| 日日啪夜夜撸| 精品人妻熟女av久视频| 大话2 男鬼变身卡| 丰满少妇做爰视频| 亚洲欧美日韩东京热| 亚洲av成人av| av在线老鸭窝| 日本色播在线视频| 久久久色成人| 夜夜看夜夜爽夜夜摸| 久久99热这里只频精品6学生| 亚洲不卡免费看| 日产精品乱码卡一卡2卡三| 午夜福利在线观看吧| 91狼人影院| 国产国拍精品亚洲av在线观看| 欧美极品一区二区三区四区| 欧美精品一区二区大全| 日本欧美国产在线视频| 国产午夜精品论理片| 97超视频在线观看视频| 熟妇人妻不卡中文字幕| 免费高清在线观看视频在线观看| 日本一二三区视频观看| 国产精品一区二区在线观看99 | 久久精品国产亚洲av天美| 国产精品国产三级专区第一集| 能在线免费看毛片的网站| 亚洲在久久综合| 丰满人妻一区二区三区视频av| 极品教师在线视频| 国产精品福利在线免费观看| 久久人人爽人人片av| 高清午夜精品一区二区三区| 好男人视频免费观看在线| 日韩成人伦理影院| 国产黄a三级三级三级人| 别揉我奶头 嗯啊视频| 只有这里有精品99| 中文乱码字字幕精品一区二区三区 | 亚洲内射少妇av| 亚洲不卡免费看| 91aial.com中文字幕在线观看| 日韩亚洲欧美综合| 成年版毛片免费区| 日日摸夜夜添夜夜爱| 亚洲欧美成人精品一区二区| 国产亚洲午夜精品一区二区久久 | 免费观看av网站的网址| 亚洲高清免费不卡视频| 国产成人freesex在线| 亚洲精品国产成人久久av| 天堂网av新在线| 日韩制服骚丝袜av| 亚洲不卡免费看| 亚洲熟女精品中文字幕| 51国产日韩欧美| 日本-黄色视频高清免费观看| 亚洲欧美日韩卡通动漫| 久久久久国产网址| 国产高清有码在线观看视频| 日韩,欧美,国产一区二区三区| 亚洲最大成人av| 日韩视频在线欧美| 久久久久久久久大av| 99久国产av精品| 亚洲av国产av综合av卡| 国产淫片久久久久久久久| 国产有黄有色有爽视频| 国模一区二区三区四区视频| 亚洲成人av在线免费| 亚洲欧美日韩卡通动漫| 欧美激情久久久久久爽电影| a级毛片免费高清观看在线播放| 超碰97精品在线观看| www.av在线官网国产| 搡女人真爽免费视频火全软件| 最后的刺客免费高清国语| 精品久久久噜噜| 天天躁夜夜躁狠狠久久av| 精品久久国产蜜桃| 最近2019中文字幕mv第一页| 麻豆成人av视频| 久久国产乱子免费精品| 亚洲高清免费不卡视频| 日韩视频在线欧美| 蜜臀久久99精品久久宅男| 观看美女的网站| 免费黄频网站在线观看国产| 男女啪啪激烈高潮av片| 一个人免费在线观看电影| 亚洲国产精品国产精品| 国产女主播在线喷水免费视频网站 | 亚洲精品第二区| 边亲边吃奶的免费视频| 简卡轻食公司| 欧美三级亚洲精品| 国产精品熟女久久久久浪| 久久久久久久久久人人人人人人| 欧美激情在线99| 久久精品国产亚洲av涩爱| 日韩强制内射视频| 欧美变态另类bdsm刘玥| 搡老乐熟女国产| 国产激情偷乱视频一区二区| 精品久久国产蜜桃| 久久久久久久久中文| 亚洲内射少妇av| 精品一区二区三卡| 亚洲国产欧美人成| 国产午夜精品一二区理论片| 嫩草影院精品99| 欧美高清成人免费视频www| 国产亚洲午夜精品一区二区久久 | 22中文网久久字幕| 最近的中文字幕免费完整| 国产高潮美女av| 国产淫片久久久久久久久| 人妻夜夜爽99麻豆av| 国产精品嫩草影院av在线观看| 肉色欧美久久久久久久蜜桃 | 久久久久久久亚洲中文字幕| 午夜激情久久久久久久| 久久精品熟女亚洲av麻豆精品 | 久久久久久久久久久丰满| 成人国产麻豆网| 亚洲美女搞黄在线观看| 亚洲国产精品专区欧美| videos熟女内射| 一级毛片 在线播放| 尤物成人国产欧美一区二区三区| 高清日韩中文字幕在线| 亚洲在线观看片| 国产精品人妻久久久影院| 最近最新中文字幕大全电影3| 99久久中文字幕三级久久日本| 国产在线一区二区三区精| 久久久成人免费电影| 亚洲不卡免费看| 九九久久精品国产亚洲av麻豆| 精品久久久久久成人av| 99久久精品热视频| 日韩 亚洲 欧美在线| 纵有疾风起免费观看全集完整版 | 全区人妻精品视频| 国产精品嫩草影院av在线观看| 午夜福利在线观看免费完整高清在| 99久国产av精品国产电影| 精品熟女少妇av免费看| 在线免费十八禁| 中文天堂在线官网| 国产乱来视频区| 免费高清在线观看视频在线观看| 亚洲精品456在线播放app| 精华霜和精华液先用哪个| 女的被弄到高潮叫床怎么办| 久久久成人免费电影| 国产综合懂色| 国产精品一区二区三区四区久久| 秋霞伦理黄片| 国产精品精品国产色婷婷| 久久久久免费精品人妻一区二区| 国产综合精华液| 99re6热这里在线精品视频| 91在线精品国自产拍蜜月| 少妇高潮的动态图| 能在线免费看毛片的网站| 夜夜看夜夜爽夜夜摸| 一本一本综合久久| 亚洲精品视频女| 三级国产精品片| 国产黄频视频在线观看| 午夜精品国产一区二区电影 | 久久精品综合一区二区三区| 国产精品不卡视频一区二区| 国产av国产精品国产| 亚洲成色77777| 国语对白做爰xxxⅹ性视频网站| 欧美三级亚洲精品| 国产麻豆成人av免费视频| 神马国产精品三级电影在线观看| 免费看av在线观看网站| 国产伦精品一区二区三区视频9| 蜜桃久久精品国产亚洲av| 久久99热这里只有精品18| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品色激情综合| 成人欧美大片| 视频中文字幕在线观看| .国产精品久久| 禁无遮挡网站| 男女啪啪激烈高潮av片| 国产精品久久久久久久久免| 久久久久久久久大av| 天堂网av新在线| 久久久精品94久久精品| 不卡视频在线观看欧美| 亚洲高清免费不卡视频| 我的老师免费观看完整版| 免费黄色在线免费观看| 黄片wwwwww| 午夜久久久久精精品| 中文精品一卡2卡3卡4更新| 小蜜桃在线观看免费完整版高清| 久久99热这里只频精品6学生| 国产黄片美女视频| 五月天丁香电影| 日日啪夜夜撸| 免费黄网站久久成人精品| 日韩大片免费观看网站| 精品久久久久久久人妻蜜臀av| 免费看美女性在线毛片视频| av天堂中文字幕网| 男人爽女人下面视频在线观看| 亚洲精品视频女| 久热久热在线精品观看| av免费在线看不卡| 精品久久久久久久人妻蜜臀av| 欧美日韩视频高清一区二区三区二| 色综合色国产| 老司机影院成人| 久久热精品热| 亚洲精品中文字幕在线视频 | 99久国产av精品| 国产成人91sexporn| 免费少妇av软件| 久久久精品欧美日韩精品| 九草在线视频观看| 99久国产av精品国产电影| 亚洲精品亚洲一区二区| 亚洲av不卡在线观看| 最近最新中文字幕免费大全7| 午夜免费观看性视频| 乱人视频在线观看| 国产午夜精品论理片| 在线播放无遮挡| 亚洲,欧美,日韩| 亚洲精品国产成人久久av| 欧美潮喷喷水| 国产成人精品婷婷| 亚洲av成人精品一区久久| 国产有黄有色有爽视频| 丰满乱子伦码专区| 国产淫语在线视频| 亚洲欧洲国产日韩| 简卡轻食公司| 国产精品久久久久久久电影| 亚洲综合精品二区| 成人亚洲欧美一区二区av| av国产久精品久网站免费入址| 国产精品一区二区在线观看99 | 少妇人妻精品综合一区二区| ponron亚洲| 乱系列少妇在线播放| 两个人的视频大全免费| 好男人在线观看高清免费视频| 大话2 男鬼变身卡| 国产精品嫩草影院av在线观看| 国产成人精品久久久久久| 久久久久九九精品影院| 最近中文字幕2019免费版| 成人特级av手机在线观看| 久久午夜福利片| 特级一级黄色大片| 深爱激情五月婷婷| 80岁老熟妇乱子伦牲交| 国产一区亚洲一区在线观看| 伊人久久国产一区二区| 99久久精品热视频| 伊人久久国产一区二区| 99热网站在线观看| 久久久久久久久久久丰满| 草草在线视频免费看| 亚洲精品一区蜜桃| 久久这里只有精品中国| 九九爱精品视频在线观看| 热99在线观看视频| 亚洲精品乱码久久久v下载方式| 亚洲欧洲国产日韩| 色尼玛亚洲综合影院| 亚洲综合色惰| 国产 一区 欧美 日韩| 国产精品人妻久久久影院| 一级爰片在线观看| 黑人高潮一二区| 一级爰片在线观看| 2018国产大陆天天弄谢| 日韩一区二区视频免费看| 精品一区在线观看国产| 国产成人午夜福利电影在线观看| 少妇的逼好多水| 尾随美女入室| 男插女下体视频免费在线播放| 免费高清在线观看视频在线观看| videos熟女内射| 国产精品精品国产色婷婷| 美女脱内裤让男人舔精品视频| 最近的中文字幕免费完整| 草草在线视频免费看| 波多野结衣巨乳人妻| 亚洲高清免费不卡视频| 久久久午夜欧美精品| av黄色大香蕉| 国产精品女同一区二区软件| 国产午夜福利久久久久久| 又爽又黄无遮挡网站| 久久久久久久久久黄片| 久久精品国产亚洲av涩爱| 日韩亚洲欧美综合| 看非洲黑人一级黄片| 精品国内亚洲2022精品成人| 91久久精品国产一区二区成人| 精品人妻偷拍中文字幕| 色尼玛亚洲综合影院| 狂野欧美激情性xxxx在线观看| av天堂中文字幕网| 国产视频内射| 国产 亚洲一区二区三区 | 国产大屁股一区二区在线视频| 亚洲最大成人av| 国产有黄有色有爽视频| 成人二区视频| 国产av不卡久久| 少妇熟女aⅴ在线视频| 在线免费十八禁| 亚洲欧美成人综合另类久久久| 国产午夜精品一二区理论片| 一区二区三区乱码不卡18| 亚洲自拍偷在线| 国产成人a区在线观看| 两个人视频免费观看高清| 22中文网久久字幕| 伊人久久国产一区二区| 超碰av人人做人人爽久久| 精品一区二区三区人妻视频| 亚洲人成网站在线播| 黑人高潮一二区| 成人高潮视频无遮挡免费网站| 人妻制服诱惑在线中文字幕| 简卡轻食公司| 九九久久精品国产亚洲av麻豆| 91精品一卡2卡3卡4卡| 国产在视频线精品| 精华霜和精华液先用哪个| av专区在线播放| 97热精品久久久久久| 国产精品福利在线免费观看| 欧美成人午夜免费资源| 免费不卡的大黄色大毛片视频在线观看 | 人人妻人人澡欧美一区二区| 午夜精品国产一区二区电影 | 淫秽高清视频在线观看| 免费电影在线观看免费观看| 青春草视频在线免费观看| 淫秽高清视频在线观看| 少妇高潮的动态图| 丰满少妇做爰视频| 国产永久视频网站| 欧美 日韩 精品 国产| 日韩三级伦理在线观看| 久久久久免费精品人妻一区二区| 777米奇影视久久| 亚洲,欧美,日韩| 乱系列少妇在线播放| 亚洲av免费在线观看| 亚洲国产精品成人综合色| 最后的刺客免费高清国语| 欧美潮喷喷水| 中文天堂在线官网| 婷婷色麻豆天堂久久| 七月丁香在线播放| 成年人午夜在线观看视频 | 成年免费大片在线观看| 少妇熟女aⅴ在线视频| 欧美bdsm另类| 午夜爱爱视频在线播放| 伊人久久精品亚洲午夜| 久久99热这里只有精品18| 久久久久精品性色| 18+在线观看网站| 看非洲黑人一级黄片| 日本一二三区视频观看|