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

面向制動踏板感覺的主缸動力學(xué)模型及其關(guān)鍵影響因素

2015-06-14 07:38:06孟德建張立軍方明霞余卓平
關(guān)鍵詞:踏板

孟德建,張立軍,方明霞,余卓平

(1.同濟大學(xué) 汽車學(xué)院,上海201804;2.同濟大學(xué) 智能型新能源汽車協(xié)同創(chuàng)新中心,上海201804;3.同濟大學(xué) 力學(xué)流動站,上海200092;4.同濟大學(xué),航空航天與力學(xué)學(xué)院,上海200092)

0 引 言

制動踏板感覺作為駕駛員與汽車之間交互的重要組成,直接影響到消費者對車輛品質(zhì)的評價。因此,現(xiàn)在的汽車制動系統(tǒng)不僅要滿足汽車制動安全性,還需要具有良好的制動踏板感覺。主缸是將機械踏板力轉(zhuǎn)化為制動液壓力的關(guān)鍵部件,在制動系統(tǒng)中起到“變壓器”的作用[1-2]。整車制動踏板感覺經(jīng)常使用踏板行程與制動油壓關(guān)系、踏板力與制動油壓關(guān)系進行評價,顯然體現(xiàn)了主缸對制動踏板感覺的重要作用。因此,面向制動踏板感覺的主缸建模與關(guān)鍵影響因素分析,對考慮制動踏板感覺的制動系統(tǒng)正向開發(fā)與設(shè)計具有重要的意義。

中心閥式串列雙腔主缸是現(xiàn)代汽車制動主缸的主要類型,它主要包括第一活塞、第一活塞回位彈簧、第二活塞、第二活塞回位彈簧、中心閥、皮碗和缸體等結(jié)構(gòu),制動液填充在活塞腔內(nèi)。針對這一包含彈簧機械系統(tǒng)和黏性液體的機構(gòu),Khan等[1]采用鍵合圖的方法建立了機械系統(tǒng)的動力學(xué)模型,Gerdes和詹軍等[2-3]則建立了機械系統(tǒng)動力學(xué)數(shù)學(xué)模型。Fisher[4]考慮活塞回位彈簧的預(yù)緊力和系統(tǒng)受到的摩擦力,利用液體體積彈性模量的定義,建立了包含機械系統(tǒng)模型和制動液模型的制動主缸動力學(xué)數(shù)學(xué)模型。在此基礎(chǔ)之上,Yamada等[5]利用一維流體連續(xù)方程建立了制動液的數(shù)學(xué)模型,劉曉[6]重點考慮了補償空位置建立了制動液的數(shù)學(xué)模型,他們都使用Matlab/Simulink軟件對機械系統(tǒng)和制動液數(shù)學(xué)模型進行了仿真。根據(jù)制動主缸的拓撲結(jié)構(gòu),F(xiàn)ortina和賀平 等[7-8]分 別 利 用AMEsim 軟 件 和Modelica/MWorks軟件建立了制動主缸動力學(xué)模型。Day等[9]進一步建立了AMEsim 和Abaqus聯(lián)合仿真的制動主缸動力學(xué)模型,重點分析了橡膠皮腕變形對制動主缸特性的影響。方泳龍等[10]建立了制動主缸行程知識庫,王興東等[11]利用分段函數(shù)建立了汽液混合型制動主缸動力學(xué)模型,為筆者建立制動主缸動力學(xué)模型提供了借鑒。

通過分析發(fā)現(xiàn):①現(xiàn)有的制動主缸模型均是針對制動性能而開發(fā)的,它們往往關(guān)注制動速度、制動距離、制動減速度、制動熱效能等整車制動性能指標[3-4,6],或者是制動油壓、推桿力和推桿行程的時域 特 性[1-2,7-8,10-11]。而 制 動 踏 板 感 覺 往 往 考察主缸推桿力-行程、主缸油壓-行程、主缸油壓-推桿力的關(guān)系[9]。②面向制動踏板感覺的制動主缸特性高度依賴于主缸推桿速度,但是現(xiàn)有研究基本沒有考慮推桿速度對制動踏板特性的影響。③現(xiàn)有的制動主缸的動力學(xué)模型存在不足。例如,在建立機械系統(tǒng)模型時,文獻[6]忽略了系統(tǒng)摩擦力,文獻[1-3,5-6]忽略了回位彈簧的預(yù)緊力,而文獻[4]未給出活塞開始工作的判定條件;在建立制動液模型時,往往忽略體積彈性模量的變化[5-6];而在使用AMEsim 等軟件建模時,沒能給出反應(yīng)系統(tǒng)特性的數(shù)學(xué)表達式[7-9]。

為了完成面向制動踏板感覺的制動系統(tǒng)設(shè)計與開發(fā),本文在考慮回位彈簧的預(yù)緊力、系統(tǒng)摩擦力、閥口間隙以及制動液體積彈性模量的變化的基礎(chǔ)上,建立了進程階段面向制動踏板感覺的主缸動力學(xué)模型。基于試驗辨識了模型的關(guān)鍵參數(shù),并在不同推桿速度下驗證了模型的有效性。以推桿力-行程特性、主缸油壓-行程特性和主缸油壓-推桿力特性形成的3象限圖為評價指標,分析了影響制動主缸特性的關(guān)鍵因素。

1 試驗方案

在制動踏板感覺試驗臺架上,利用伺服電動缸對制動主缸施加推力,使用拉壓力傳感器測量推桿力,拉桿位移傳感器測量推桿位移,油壓傳感器測量活塞腔油壓,傳感器的布置如圖1 所示。為了充分辨識回位彈簧剛度和活塞所受摩擦力,分別在無制動液和有制動液兩種情況下開展試驗。這兩種工況下推桿速度分為5個水平,分別為4、6、8、10、20mm/s。在無制動液工況時不測量油壓,有制動液時將出液口堵住。

圖1 主缸試驗傳感器布置Fig.1 Sensors layout in master cylinder experiment

2 數(shù)學(xué)模型與參數(shù)辨識

2.1 研究對象與基本假設(shè)

中心閥式串列雙腔制動主缸的典型結(jié)構(gòu)及其關(guān)鍵零件如圖2所示。為了建立制動主缸動力學(xué)模型,作基本假設(shè)如下:①只考慮皮腕與缸體之間的摩擦力,忽略制動液及其流動產(chǎn)生的摩擦力;②工作過程中制動液溫度保持不變且為室溫;③忽略制動液流動在活塞腔內(nèi)造成的壓力損失,認為同一腔內(nèi)壓力處處相等。

圖2 主缸結(jié)構(gòu)示意圖Fig.2 Structure schematic diagram of master cylinder

2.2 機械系統(tǒng)動力學(xué)模型

根據(jù)主缸結(jié)構(gòu)和工作原理,可以將其等效為圖3所示的物理模型。主缸的運動過程包含3個階段,分別建立各個階段的動力學(xué)方程,模型中包含的關(guān)鍵零件及其參數(shù)說明如下:mc1、mc2分別為第一、第二活塞質(zhì)量;xc1、xc2分別為第一、第二活塞位移;Fe為第一活塞推力;kc1、kc2分別為第一、第二活塞回位彈簧剛度;cc1、cc2分別為第一、第二活塞阻尼;Fm1、Fm2分別為第一、第二活塞回位彈簧預(yù)緊力;Fu1、Fu2分別為第一、第二活塞摩擦力;P1、P2分別為第一、第二活塞腔油壓。

圖3 主缸物理模型Fig.3 Physical model of master cylinder

第一階段:主缸推桿力增加,克服第二活塞回位彈簧預(yù)緊力后開始運動,第二活塞回位彈簧開始被壓縮,直至第一活塞回位彈簧開始被壓縮為止。該過程中可將第一活塞、第一活塞回位彈簧和第二活塞看作一個整體。制動主缸的動力學(xué)方程為:

式中:xc11為第一階段末第一活塞的位移;FN1為第二階段第一活塞受到的初始力。

第二階段:第一活塞回位彈簧開始被壓縮,直至第二活塞回位彈簧被壓縮至極限。該過程中前、后活塞以不同的速度共同運動,第一活塞動力學(xué)方程為:

第二活塞動力學(xué)方程為:

式中:xc12、xc22分別為第二階段末第一、第二活塞的位移,F(xiàn)N2為第三階段第一活塞受到的初始力。

第三階段:第二活塞回位彈簧被壓縮至極限,直至第一活塞回位彈簧也被壓縮至極限。該過程中只有第一活塞在運動,制動主缸的動力學(xué)方程為:

式中:h1為第一活塞腔有效長度;sc為活塞面積。

2.3 制動液動力學(xué)模型

由一維流體連續(xù)方程可得:

式中:ρ為流體的密度;vx為流體沿方向的流速;x 為流體沿方向的位移;t為時間。

令Q 為流體的體積流量,S 為過流面積,則式(7)可表達為:

對于具有規(guī)則形狀的腔體,假設(shè)流體的體積流量均勻降低,則式(8)可表示為:

式中:Qin為流進的流體流量;Qout為流出的流體流量;V 為腔體體積。

制動液是可壓縮性液體,其壓縮能力用體積彈性模量Ke表示:

式中:P 為流體壓力。

制動液動力學(xué)方程為:

由于實際試驗過程中將主缸出液口堵住,因此Qout=0。所以,第一階段制動液動力學(xué)方程為:

第二階段制動液動力學(xué)方程為:

第三階段制動液動力學(xué)方程為:

式中:hb1、hb2為初始狀態(tài)下第一、第二中心閥口間隙。

制動液體積彈性模量Ke隨著含氣量和壓力變化,工程中制動液的體積彈性模量[12-13]為:

式中:Kl為純制動液的體積彈性模量;δ0為標準大氣壓下空氣的溶解度;γ為氣體絕熱指數(shù);α為制動液中實際的氣體含量;P0為標準大氣壓。

2.4 基于試驗的關(guān)鍵參數(shù)辨識

2.4.1 無制動液工況

在無制動液工況下,主要辨識前回位彈簧的剛度和預(yù)緊力以及活塞受到的摩擦力。圖4為無制動液時主缸的推桿力-行程曲線。由圖4可知:推桿力-行程特性可以分為0~14 mm、14~25 mm、25~35mm 三個階段,其中第三階段中主要壓縮第一活塞回位彈簧,曲線的斜率即為該彈簧剛度;第一階段主要是壓縮第二活塞回位彈簧,曲線的斜率即為該彈簧剛度。

圖4 無制動液工況主缸推桿力-行程曲線Fig.4 Master cylinder pushrod force-travel curves on condition of no brake fluid

推桿速度對其制動主缸機械特性影響可以忽略不計,在相同推桿行程下進程和回程階段的推桿力差值的一半即為摩擦力,但不同推桿位置時摩擦力的大小不同,這主要是由于在第一和第二階段,摩擦力包含推桿和橡膠密封圈摩擦力和3個皮碗與缸壁的摩擦力,而在第三階段摩擦力包含推桿和橡膠密封圈摩擦力和1個皮碗與缸壁的摩擦力。分別在0~25mm、25~35 mm 內(nèi)計算各個行程下進程與回程推桿力的差值的平均值,得到制動主缸摩擦力。

推桿進程中行程為零時的推桿力為第一階段摩擦力和第二活塞回位彈簧的預(yù)緊力之和,由于摩擦力已求,從而可以計算出第二活塞回位彈簧的預(yù)緊力。在推桿行程為14mm 時,第一階段結(jié)束,第二階段開始,此時的推桿力為第一活塞回位彈簧的預(yù)緊力。

2.4.2 有制動液工況

圖5為有制動液工況制動主缸油壓-行程曲線。由圖5可知:當推桿行程為零時的油壓為制動主缸的初始油壓,制動油壓開始迅速增大時的推桿行程為制動主缸空行程,即初始狀態(tài)下第二中心閥口間隙。

經(jīng)試驗辨識的主缸參數(shù)及數(shù)值如下:Fu1=14 N;Fu2=6N;Fm1=94N;Fm2=44N;kc1=4.5N/m;kc2=1.9N/m;hb1=1.6 mm;初始油壓P0=8000Pa。

圖5 有制動液工況主缸油壓-行程曲線Fig.5 Master cylinder pushrod force-travel curves on condition of with brake fluid

3 模型有效性驗證

圖6為推桿速度為4mm/s和20mm/s且無制動液工況時制動主缸特性對比圖。由圖6 可知:無制動液工況下制動主缸模型仿真結(jié)果與試驗結(jié)果具有較好的一致性,僅在起始位置處誤差稍大,這主要是由于皮碗在運動之處的不規(guī)則變形引起的。

圖6 無制動液工況制動主缸推桿力-行程曲線對比Fig.6 Comparison of master cylinder pushrod forcetravel curves on condition of no brake fluid

圖7 為推桿速度為4mm/s和20mm/s且有制動液工況時制動主缸特性對比圖。由圖7 可知:有制動液時制動主缸動力學(xué)模型仿真結(jié)果與試驗結(jié)果的一致性較好。

圖7 有制動液工況制動主缸特性對比Fig.7 Comparison of master cylinder characteristics on condition of with brake fluid

4 關(guān)鍵影響因素

研究已經(jīng)表明[1,10,14-15],主缸結(jié)構(gòu)類型、閥口間隙、橡膠皮碗等因素對制動踏板感覺具有重要影響。本文重點針對活塞內(nèi)徑、彈簧預(yù)緊力和剛度、制動液氣體含量等因素,以推桿速度為4mm/s工況為例,分析它們對制動踏板感覺的影響。

4.1 活塞內(nèi)徑

圖8 活塞內(nèi)徑對制動主缸特性的影響Fig.8 Influence of piston inner diameter on master cylinder characteristics

圖8 為活塞內(nèi)徑對制動主缸特性的影響,其中,20mm 特性曲線為實際制動主缸特性。由圖8可知:隨著活塞內(nèi)徑的增大,主缸油壓-行程曲線沒有變化;制動主缸推桿力-行程曲線變化較小,只有在大行程時主缸推桿力-行程的剛度稍微增加,這將會使制動踏板感覺有變硬的趨勢;主缸油壓-推桿力曲線變化顯著,主缸油壓-推桿力剛度明顯減小,從而導(dǎo)致制動踏板感覺有偏軟的趨勢,這主要是由于活塞內(nèi)徑增大,相同推桿力產(chǎn)生的壓強降低。

4.2 彈簧預(yù)緊力和剛度

圖9為回位彈簧預(yù)緊力對制動主缸特性的影響。其中,前3個工況為第二活塞回位彈簧預(yù)緊力不變,依次增大第一活塞回位彈簧預(yù)緊力。通過分析可知:彈簧預(yù)緊力對主缸油壓-推桿力曲線特性影響較小,對推桿力-行程曲線特性和主缸油壓-行程曲線特性影響較大。隨著第一活塞回位彈簧預(yù)緊力的增大,主缸油壓空行程減小,推桿力-行程曲線拐點提前,從而導(dǎo)致制動踏板感覺油壓空行程減小,制動踏板力-行程曲線拐點提前。

圖9 回位彈簧預(yù)緊力對制動主缸特性的影響Fig.9 Influence of the return spring preload on master cylinder characteristics

在圖9中,第三和第四工況分別使第一、第二活塞回位彈簧預(yù)緊力增加相同倍數(shù),通過分析可知:改變第二活塞回位彈簧預(yù)緊力可以更加有效地改變推桿力-行程特性和主缸油壓-行程特性,進而說明主缸油壓特性對第二活塞回位彈簧預(yù)緊力更加敏感。

圖10為回位彈簧剛度對制動主缸特性的影響,其中,第一工況為第一活塞彈簧剛度增大,第二活塞剛度不變;第二、三、四工況為第一活塞彈簧剛度不變,第二活塞彈簧剛度依次增大。由圖10可知:彈簧剛度對主缸油壓-推桿力曲線特性影響較小,對推桿力-行程曲線特性和主缸油壓-行程曲線特性影響較大。隨著第二活塞回位彈簧剛度的增大,主缸油壓空行程減小,推桿力-行程曲線拐點提前,從而導(dǎo)致制動踏板感覺油壓空行程減小,制動踏板力-行程曲線拐點提前。增大第一活塞回位彈簧剛度時,制動主缸特性變化較小。因此,說明主缸油壓特性對第二活塞回位彈簧剛度更加敏感。

圖10 回位彈簧剛度對制動主缸特性的影響Fig.10 Influence of the return spring stiffness on master cylinder characteristics

4.3 制動液氣體含量

圖11 為制動液氣體含量對制動主缸特性的影響。由圖11 可知:制動液氣體含量對主缸油壓-推桿力特性的影響較小,但是對推桿力-行程曲線特性和主缸油壓-行程曲線特性影響較大。隨著制動液氣體含量的增加,制動主缸推桿力-行程剛度和主缸油壓-行程剛度減小,從而導(dǎo)致制動踏板感覺有變軟的趨勢;制動液氣體含量增大到一定程度后,制動主缸特性變化不明顯。因此,實際制動系統(tǒng)的排氣效果對制動主缸和制動踏板感覺具有重要的影響。

圖11 制動液氣體含量對制動主缸特性的影響Fig.11 Influence of the brake fluid gas content on master cylinder characteristics

5 結(jié)束語

考慮回位彈簧預(yù)緊力、系統(tǒng)摩擦力以及關(guān)鍵結(jié)構(gòu)參數(shù)建立了制動主缸機械動力學(xué)模型,考慮制動液體積彈性模量的可變性建立了制動液動力學(xué)模型,基于試驗辨識了系統(tǒng)的關(guān)鍵參數(shù),進而建立了面向制動踏板感覺的制動主缸動力學(xué)模型。經(jīng)不同推桿速度下的試驗對比與分析,該模型具有較好的有效性和準確性。活塞內(nèi)徑對主缸油壓-推桿力剛度影響較大;回位彈簧預(yù)緊力和剛度對推桿力-行程曲線拐點和主缸油壓空行程影響較大,第二活塞回位彈簧表現(xiàn)的更加敏感;制動液氣體含量對制動主缸特性具有重要的影響。

[1]Khan Y,Kulkarni P,Youceftoumi K.Modeling,experimentation and simulation of a brake apply system[J].Journal of Dynamic Systems Measurement and Control,1994,116(1):111-122.

[2]Gerdes J C,Hedrick J K.Brake system modeling for simulation and control[J].Journal of Dynamic Systems Measurement and Control,1999,121(3):496-503.

[3]詹軍.基于ACC的制動系統(tǒng)模型研究[J].中國機械工程,2005,16(5):76-78.Zhan Jun.Research on model of brake system based on ACC system[J].China Mechanical Engineering,2005,16(5):76-78.

[4]Fisher D K.Brake system component dynamic performance measurement and analysis[C]∥International Automobile Safety Conference,Detroit,Michigan,1970.

[5]Yamada T,Sawada M.Development and implementation of simulation tool for vehicle brake system[C]∥SAE 2001 World Congress,Detroit,Michigan,2001.

[6]劉曉.基于特性的制動系統(tǒng)動態(tài)建模與分析[D].長春:吉林大學(xué)汽車工程學(xué)院,2012.Liu Xiao.Dynamic modeling and analysis of the braking system based on its characteristics[D].Changchun:College of Automotive Engineering,Jilin University,2012.

[7]Fortina A,Velardocchia M,Sorniotti A.Braking system components modelling[C]∥21st Annual Brake Colloquium &Exhibition,Hollywood,F(xiàn)lorida,2003.

[8]賀平.基于Modelica的液壓制動系統(tǒng)多領(lǐng)域建模與仿真[D].武漢:華中科技大學(xué)機械科學(xué)與工程學(xué)院,2011.He Ping.The multi-domain modeling and simulation of hydraulic brake system based on modelica[D].Wuhan:School of Mechanical Science and Engineering,Huazhong University of Science &Technology,2011.

[9]Day A,Ho H,Hussain K,et al.Brake system simulation to predict brake pedal feel in a passenger car[C]∥SAE 2009Brake Colloquium and Exhibition,Tampa,F(xiàn)lorida,2009.

[10]方泳龍,胥向欣,劉璽.制動主缸與踏板行程專家系統(tǒng)知識庫的建立[J].江蘇大學(xué)學(xué)報:自然科學(xué)版,2007,28(6):477-481.Fang Yong-long,Xu Xiang-xin,Liu Xi.Expert system knowledge base construction of brakemaster cylinder stroke and brake pedal stroke[J].Journal of Jiangsu University(Natural Science Edition),2007,28(6):477-481.

[11]王興東,胡于進,李成剛,等.氣液混合型制動系的執(zhí)行機構(gòu)建模試驗和仿真[J].汽車工程,2002,24(6):516-519.Wang Xing-dong,Hu Yu-jin,Li Cheng-gang,et al.Modeling,experiment and simulation of a pneumatic-h(huán)ydraulic brake system[J].Automotive Engineering,2002,24(6):516-519.

[12]朱冬.液壓系統(tǒng)低壓管路瞬態(tài)過程研究[D].哈爾濱:哈爾濱工業(yè)大學(xué)機電工程學(xué)院,2006.Zhu Dong.The research of fluid transients in low pressure hydraulic pipeline[D].Harbin:Mechatronic Engineering,Harbin Institute of Technology,2006.

[13]馮斌,龔國芳,楊華勇.液壓油彈性模量提高方法與試驗[J].農(nóng)業(yè)機械學(xué)報,2010,41(3):219-222.Feng Bin,Gong Guo-fang,Yang Hua-yong.Method and experiment for increasing effective fluid bulk modulus in hydraulic systems[J].Transactions of the Chinese Society for Agricultural Machinery,2010,41(3):219-222.

[14]Goto Y,Yasuda A,Ishida S.Brake Master cylinder for secure brake feel and improved system failure performance[C]∥21st Annual Brake Colloquium &Exhibition,Hollywood,F(xiàn)lorida,2003.

[15]Rajesh A T,Badal B,Pol S.Co-relating subjective and objective brake performance:a case study[C]∥24th Annual Brake Colloquium and Exhibition,Grapevine,Texas,2006.

猜你喜歡
踏板
單踏板不可取
車主之友(2022年6期)2023-01-30 07:58:16
鋼琴演奏中踏板運用的探討
黃河之聲(2020年4期)2020-12-07 10:48:25
汽車制動踏板感匹配方法探討
公路與汽運(2020年6期)2020-12-07 02:22:10
淺談延音踏板在鋼琴演奏中的用法
黃河之聲(2019年23期)2019-12-17 19:08:43
淺談汽車制動踏板的型面設(shè)計
淺談制動踏板機構(gòu)設(shè)計
鋼琴踏板的基本用法
淺談鋼琴踏板的運用
黃河之聲(2017年13期)2017-01-28 13:30:17
論鋼琴踏板的正確使用
淺析鋼琴演奏中踏板的使用
北方音樂(2016年22期)2016-02-04 16:52:16
主站蜘蛛池模板: 亚洲欧美成人网| 波多野结衣一区二区三视频| 亚洲婷婷丁香| 欧美精品导航| 国产日本欧美亚洲精品视| 国产成人无码综合亚洲日韩不卡| 激情网址在线观看| 精品视频一区在线观看| 视频在线观看一区二区| 久久成人18免费| 人妻无码一区二区视频| 亚洲 欧美 中文 AⅤ在线视频| 中文字幕有乳无码| 色网在线视频| 狠狠ⅴ日韩v欧美v天堂| 欧美激情视频二区| av一区二区三区在线观看| 午夜啪啪网| 伊人五月丁香综合AⅤ| 国产成人免费高清AⅤ| 亚洲国产精品日韩av专区| 天堂网亚洲系列亚洲系列| 亚洲精品国产精品乱码不卞 | 狼友av永久网站免费观看| 欧洲精品视频在线观看| 日韩欧美中文字幕在线精品| 2021国产精品自拍| 国产成人禁片在线观看| 亚洲成人福利网站| 午夜免费视频网站| 午夜久久影院| 91无码国产视频| 国产精品黄色片| 91欧洲国产日韩在线人成| 国产精品亚洲欧美日韩久久| 91无码人妻精品一区二区蜜桃| 精品视频一区二区三区在线播| 精品人妻无码区在线视频| 欧美区在线播放| 久久semm亚洲国产| 欧美一区二区人人喊爽| 在线观看无码a∨| 色精品视频| 日本人又色又爽的视频| 97精品伊人久久大香线蕉| 日韩视频免费| 九九热免费在线视频| 精品人妻一区二区三区蜜桃AⅤ| 久久久久久久久久国产精品| 亚洲天堂视频在线免费观看| 黄片一区二区三区| 女人av社区男人的天堂| 成人综合在线观看| 日本午夜网站| 精品国产成人av免费| a网站在线观看| 成人精品在线观看| 毛片网站在线看| 国产91成人| 中文字幕欧美日韩| 宅男噜噜噜66国产在线观看| 国产亚洲精品自在线| 毛片免费视频| 亚洲国模精品一区| 亚洲va在线观看| 中文字幕有乳无码| 婷婷中文在线| 国产a网站| 久久久久九九精品影院| 亚洲午夜天堂| 9久久伊人精品综合| 91小视频在线| av在线无码浏览| 国产原创自拍不卡第一页| 国产经典三级在线| 国产美女叼嘿视频免费看| 国产国语一级毛片| 久久久亚洲色| 亚洲综合经典在线一区二区| 亚洲av综合网| 日韩不卡免费视频| 国产二级毛片|