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

    某越野車車身穩(wěn)定系統(tǒng)匹配設(shè)計(jì)與試驗(yàn)研究

    2015-08-19 09:07:06方德廣張雨張邦基張農(nóng)
    汽車工程學(xué)報(bào) 2015年4期
    關(guān)鍵詞:試驗(yàn)設(shè)計(jì)

    方德廣++張雨++張邦基++張農(nóng)

    摘 要: 過(guò)硬的橫向穩(wěn)定桿會(huì)限制車輛的越野能力,因此某些車型為了保障越野能力而選用較軟的橫向穩(wěn)定桿,便導(dǎo)致了越野車側(cè)傾剛度不足的問(wèn)題。針對(duì)這一問(wèn)題設(shè)計(jì)了一套可取代橫向穩(wěn)定桿的車身穩(wěn)定系統(tǒng),用于在不影響車輛越野性能的同時(shí)增加車輛的側(cè)傾剛度。通過(guò)建立整車及液壓系統(tǒng)動(dòng)力學(xué)模型進(jìn)行數(shù)值運(yùn)算仿真,以求得能和整車匹配的車身穩(wěn)定系統(tǒng)的關(guān)鍵參數(shù)。然后根據(jù)該車的底盤結(jié)構(gòu)對(duì)系統(tǒng)進(jìn)行結(jié)構(gòu)設(shè)計(jì),并開(kāi)發(fā)出了原型樣機(jī)進(jìn)行裝車試驗(yàn)。通過(guò)對(duì)原裝越野車和改裝越野車進(jìn)行懸架性能試驗(yàn)、蛇行試驗(yàn)和平順性試驗(yàn),并對(duì)結(jié)果作對(duì)比分析,驗(yàn)證了車身穩(wěn)定系統(tǒng)可以大幅度提高車輛動(dòng)態(tài)側(cè)傾剛度,改善其操縱穩(wěn)定性而并不影響平順性。

    關(guān)鍵詞:車身穩(wěn)定系統(tǒng);設(shè)計(jì);試驗(yàn);操縱穩(wěn)定性;平順性

    中圖分類號(hào):U461.6 文獻(xiàn)標(biāo)識(shí)碼:A

    Body Stabilization System Design and Experimental Investigation Based on an Off-Road Vehicle

    Fang Deguang1,Zhang Yu2,Zhang Bangji2,Zhang Nong2

    (1.Nanjing Iveco Motor Company Ltd, Nanjing 210028, Jiangsu,China;

    2. State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body, Hunan University, Changsha 410082, Hunan,China)

    Abstract:

    Stiffer anti-roll bars limit the off-road performance of vehicles, so to ensure the off-road capability, some vehicles are installed with soft anti-roll bars which are lack of the roll-plane stiffness. A body stabilizer system is designed based on an off-road vehicle to replace anti-roll bars to improve roll stiffness and have no effect on the off-road performance. A hydraulic system model and a full-car model are employed to conduct numerical study to obtain the key parameters for the body stabilizer system to match the vehicle. Furthermore, based on the obtained parameters, a protoyphy of the system is then designed and installed to the vehicle chassis. A set of tests are performed to evaluate the effects of the body stabilizer system on suspension performance, vehicle handling and ride performance. Overall, the experimental investigation indicated that the body stabilizer system can significantly improve the vehicle dynamic roll stability and handling performance without compensating ride performance.

    Key words: Body stabilization system; design; experiment; handling stability; ride performance

    懸架是車輛最重要的子系統(tǒng)之一,直接決定車輛的操縱穩(wěn)定性和行駛平順性。而車輛的操縱穩(wěn)定性和行駛平順性是相互制約的,難以同時(shí)滿足,因此如何在兩者之間取得合理的平衡以達(dá)到最好的效果一直是懸架研究的重要方向[1-3]。對(duì)于汽車操縱穩(wěn)定性和行駛平順性的優(yōu)化問(wèn)題也成為目前底盤開(kāi)發(fā)設(shè)計(jì)中的重要課題[4-6]。橫向穩(wěn)定桿以機(jī)械式互聯(lián)方式將左右車輪互聯(lián)提高車輛側(cè)傾穩(wěn)定性,但是受其原理限制,橫向穩(wěn)定桿的剛度不可以設(shè)計(jì)的過(guò)高。最重要的是,在越野路面穩(wěn)定桿會(huì)限制車輛對(duì)角車輪的運(yùn)動(dòng),進(jìn)而影響車輛的越野性能[5]。因此,發(fā)展可提高車輛操縱性又不影響車輛平順性和越野能力的懸架系統(tǒng)是越野車輛亟待解決的難題[7-8]。

    本文所涉及的特種越野車在實(shí)際投入使用后發(fā)現(xiàn)車輛經(jīng)過(guò)彎道時(shí)車身的傾斜角較大,車輛的側(cè)傾穩(wěn)定性有待提高。在原越野車現(xiàn)有底盤基礎(chǔ)上改裝基于液壓互聯(lián)懸架原理的車身穩(wěn)定系統(tǒng)[9-12],通過(guò)對(duì)比裝有車身穩(wěn)定系統(tǒng)的樣車與原車的懸架性能試驗(yàn)、蛇形試驗(yàn)及平順性試驗(yàn),驗(yàn)證了車身穩(wěn)定系統(tǒng)能夠有效地提高車輛操縱穩(wěn)定性。

    1 越野車相關(guān)性能試驗(yàn)方法

    研究對(duì)象為某一改進(jìn)型的越野車樣車,如圖1所示。該越野車的前懸架采用扭桿彈簧,后懸架采用鋼板彈簧,前后懸架分別匹配有橫向穩(wěn)定桿。

    圖1 某越野車樣車

    1.1越野車懸架性能試驗(yàn)與分析

    車輛懸架性能的一個(gè)評(píng)價(jià)指標(biāo)為車輛的側(cè)傾剛度 [1]。根據(jù)汽車懸掛系統(tǒng)固有頻率測(cè)試方法,參照國(guó)標(biāo)GB/T 4783-1984《汽車懸掛系統(tǒng)的固有頻率和阻尼比測(cè)定方法》,試驗(yàn)過(guò)程中記錄了車輛的自由衰減響應(yīng)曲線,通過(guò)模態(tài)參數(shù)識(shí)別得到車輛的側(cè)傾固有頻率。汽車懸架系統(tǒng)的側(cè)傾固有頻率與懸架系統(tǒng)的剛度為正比關(guān)系,反映汽車的抗側(cè)傾性能,關(guān)系到汽車的操縱穩(wěn)定性[2]。

    1.2越野車操縱穩(wěn)定性試驗(yàn)與分析

    車身側(cè)傾角的大小與懸架系統(tǒng)的剛度密切相關(guān)。車輛在彎道上行駛受到側(cè)向力作用,車身側(cè)傾角過(guò)大會(huì)使乘員產(chǎn)生恐懼感,一旦側(cè)傾角達(dá)到臨界值,車輛側(cè)翻事故便不可避免。車輛操縱穩(wěn)定性蛇行試驗(yàn)屬于駕駛員—汽車—外界環(huán)境組合而成的閉路系統(tǒng)性能試驗(yàn),可以反映車輛急劇的轉(zhuǎn)向能力和乘員的舒適性和安全性[2-4]。參照國(guó)標(biāo)GB/T 6323.1-1994《汽車操穩(wěn)性試驗(yàn)方法-蛇行試驗(yàn)》,試驗(yàn)過(guò)程中記錄了行駛車速、側(cè)向加速度和車身側(cè)傾角數(shù)值。

    1.3越野車平順性試驗(yàn)與分析

    車輛的平順性由安裝于車輛特定部位的加速度傳感器的測(cè)量信號(hào)進(jìn)行評(píng)估,根據(jù)隨機(jī)過(guò)程理論,某一時(shí)域信號(hào)的均方根值等于其功率譜密度在整個(gè)頻率范圍的積分的開(kāi)方值,因此,加速度均方根值可由頻域積分的方法求出[4-6]。參照國(guó)標(biāo)GB/T 4970-2009 《汽車平順性隨機(jī)輸入行駛試驗(yàn)方法》,選取車輛駕駛員座椅處、質(zhì)心、后右座椅處和后左座椅處為振動(dòng)加速度測(cè)點(diǎn),用加速度均方根值作為整車振動(dòng)的評(píng)價(jià)指標(biāo)。

    2車身穩(wěn)定系統(tǒng)與整車的匹配設(shè)計(jì)

    2.1匹配設(shè)計(jì)流程

    液壓系統(tǒng)關(guān)鍵物理元件性能參數(shù)的取值來(lái)源于樣C件的性能參數(shù)與設(shè)計(jì)參數(shù)的組合,是離散數(shù)據(jù)組合。該離散數(shù)據(jù)組合構(gòu)成了車身穩(wěn)定系統(tǒng)設(shè)計(jì)的設(shè)計(jì)變量集。因此車身穩(wěn)定系統(tǒng)參數(shù)設(shè)計(jì)的目標(biāo)是為尋求該變量集內(nèi)的最佳參數(shù)組合,實(shí)現(xiàn)最佳車輛性能。車身穩(wěn)定系統(tǒng)設(shè)計(jì)實(shí)質(zhì)為多目標(biāo)設(shè)計(jì)優(yōu)化問(wèn)題,提出一種針對(duì)車身穩(wěn)定系統(tǒng)的設(shè)計(jì)方法來(lái)獲取懸架參數(shù)的最優(yōu)組合解。

    針對(duì)整車性能需要達(dá)到的目標(biāo),根據(jù)車身穩(wěn)定系統(tǒng)提供的抗側(cè)傾力矩初步計(jì)算液壓缸無(wú)桿腔和有桿腔的截面積,蓄能器預(yù)充壓力和預(yù)充氣體體積,并預(yù)設(shè)蓄能器的工作壓力,計(jì)算在給定的等效垂向運(yùn)動(dòng)、側(cè)傾運(yùn)動(dòng)轉(zhuǎn)角下蓄能器的最大壓力。從蓄能器的結(jié)構(gòu)強(qiáng)度和抗疲勞角度出發(fā),將蓄能器工作時(shí)的工作壓力限定在0.9倍預(yù)充壓力和4倍的預(yù)充壓力范圍之間[8],先更改工作壓力以滿足整車性能要求,若改變?nèi)匀粺o(wú)法達(dá)到此要求,則改變蓄能器的預(yù)充壓力和預(yù)充體積參數(shù);若還是沒(méi)法滿足,則改變液壓缸有桿腔和無(wú)桿腔的截面積。最后取滿足整車性能目標(biāo)的液壓系統(tǒng)參數(shù)作為最后匹配確定的設(shè)計(jì)參數(shù)。

    車身穩(wěn)定系統(tǒng)的參數(shù)設(shè)計(jì)流程圖如圖2所示,

    圖2 車身動(dòng)態(tài)穩(wěn)定系統(tǒng)相關(guān)參數(shù)設(shè)計(jì)流程

    通過(guò)參數(shù)識(shí)別方法[13]得到原越野車的整車動(dòng)力學(xué)參數(shù)以便進(jìn)行參數(shù)匹配模型的仿真計(jì)算,獲得的整車動(dòng)力學(xué)參數(shù)見(jiàn)表1。

    表1車身動(dòng)態(tài)穩(wěn)定系統(tǒng)部分參數(shù)

    參數(shù)

    簧上質(zhì)量 ms/kg

    2 768

    俯仰轉(zhuǎn)動(dòng)慣量 Iyy/(kg·m2)

    2 923.1

    側(cè)傾轉(zhuǎn)動(dòng)慣量 Ixx/(kg·m2)

    1 783.1 kg.

    前簧下質(zhì)量 muf/kg

    105

    后簧下質(zhì)量 mur/kg

    130

    前懸剛度 ksf y/(kNm/rad)

    95.75

    后懸剛度 ksr(kNm/rad)

    112.665

    前輪剛度 ktf/(N·m)

    410 000

    后輪剛度 ktr f/(N·m)

    410 000

    前懸左右距離 lff /m 0.9

    后懸左右距離 lr/m 0.9

    前軸到質(zhì)心距離 a/m 1.3

    后軸到質(zhì)心距離 b/m 1.5

    在原越野車現(xiàn)有底盤基礎(chǔ)上改裝車身穩(wěn)定系統(tǒng)[9-12],需要提高車輛動(dòng)態(tài)側(cè)傾剛度,改善其操縱穩(wěn)定性而并不影響平順性,通過(guò)對(duì)比裝有車身穩(wěn)定系統(tǒng)的樣車與原車的懸架性能試驗(yàn)、蛇形試驗(yàn)及平順性試驗(yàn),分析并驗(yàn)證車身穩(wěn)定系統(tǒng)的有效性。

    2.2車身穩(wěn)定系統(tǒng)參數(shù)匹配設(shè)計(jì)

    車身穩(wěn)定系統(tǒng)提供的側(cè)傾剛度與蓄能器預(yù)充壓力、預(yù)充體積以及工作壓力、液壓缸大小共同決定。

    針對(duì)車輛性能相關(guān)參數(shù)的設(shè)計(jì)過(guò)程,定義設(shè)計(jì)變量

    。 (1)

    式中,P0,V0分別為蓄能器預(yù)充氮?dú)鈮毫腕w積;P1為蓄能器工作壓力;D,d分別為液壓缸內(nèi)徑和活塞桿直徑。

    根據(jù)上述分析可知,參數(shù)匹配的評(píng)估目標(biāo)應(yīng)該為設(shè)計(jì)過(guò)程中獲得的整車性能指標(biāo)的集合,如懸架系統(tǒng)的側(cè)傾剛度和側(cè)傾固有頻率,整車操穩(wěn)性中車身側(cè)傾角,整車平順性試驗(yàn)中車身質(zhì)心處的加速度響應(yīng)。定義參數(shù)匹配目標(biāo)

    。 (2)

    。 (3)

    。 (4)

    。 (5)

    式中,K_φ,K ?_φ,ω_nroll,ω ?_nroll,φ,φ ?,a_ω,a ?_ω分別為安裝車身穩(wěn)定系統(tǒng)車輛和原車的側(cè)傾角剛度,側(cè)傾固有頻率,車身側(cè)傾角和車身垂向加速度響應(yīng);ei為整車性能各目標(biāo)的范圍。

    車身穩(wěn)定系統(tǒng)的設(shè)計(jì)必須處理整車性能的折衷,使得Ji(i=1, 2, 3, 4)均達(dá)到理想中最小值,從而得到實(shí)現(xiàn)整車性能最優(yōu)的最佳車身穩(wěn)定系統(tǒng)參數(shù)組合。

    定義車身穩(wěn)定系統(tǒng)設(shè)計(jì)變量x的約束范圍

    。 (6)

    。 (7)

    式中,C為系數(shù)矩陣。

    2.2.1車身穩(wěn)定系統(tǒng)設(shè)計(jì)參數(shù)初始值

    車身穩(wěn)定系統(tǒng)由安裝于車輪(或車軸)與車身之間的單向或者雙向作用液壓缸、連接各液壓缸之間的油管以及安裝于油路中間的蓄能器組成,油管的連接方式取決于車輛需要解決的動(dòng)力學(xué)問(wèn)題[13]。

    綜合考慮車輛的操縱穩(wěn)定性以及行駛平順性,結(jié)合車輛系統(tǒng)動(dòng)力學(xué)分析,車身穩(wěn)定系統(tǒng)作用原理圖如圖3所示。

    圖3 車身穩(wěn)定系統(tǒng)原理圖

    車身穩(wěn)定系統(tǒng)中蓄能器需要保證補(bǔ)充和儲(chǔ)存系統(tǒng)中的油液,蓄能器體積V0的約束范圍如式(8)所示,懸架最大行程為Smax=152 mm,

    。 (8)

    式中, A_T=πD^2/4,A_B=π〖(D〗^2-d^2)/4; AT,AB分別為液壓缸無(wú)桿腔和有桿腔截面積。

    蓄能器預(yù)充氮?dú)鈮毫0的約束范圍如式(9)所示,蓄能器壓縮比為1:8,

    P_0 〖≤P〗_1≤8P_0 。 (9)

    液壓缸密封件往復(fù)運(yùn)動(dòng)最高壓力P_max=20 MPa。系統(tǒng)工作壓力越高對(duì)液壓系統(tǒng)中密封件的要求也越高。

    液壓缸無(wú)桿腔和有桿腔截面積經(jīng)驗(yàn)比A_T?(A_B=0.8)。根據(jù)原越野車前懸架的設(shè)計(jì)幾何關(guān)系,液壓缸外徑D_max≤0.065m。

    根據(jù)車身穩(wěn)定系統(tǒng)在穩(wěn)態(tài)條件下能夠達(dá)到的整車性能目標(biāo)初步確定車身穩(wěn)定系統(tǒng)變量的初始值x0=[x10, x20, x30, x40, x50]T=[AT0, AB0, P00, V00, P10]T 。車身穩(wěn)定系統(tǒng)設(shè)計(jì)參數(shù)的初始值和約束范圍見(jiàn)表2。

    表2車身穩(wěn)定系統(tǒng)設(shè)計(jì)參數(shù)變量

    初始值

    (x0) 下限

    (xmin) 上限

    (xmax)

    液壓缸內(nèi)徑 D/m 0.05 0.04 0.05

    活塞桿直徑 d/m 0.014 0.014 0.02

    蓄能器體積 V0/L 0.75 0.403 0.8

    蓄能器預(yù)充壓力 P0/MPa 1 0.5 2.5

    系統(tǒng)工作壓力 P1/MPa 2.5 1 4

    車身穩(wěn)定系統(tǒng)穩(wěn)定狀態(tài)時(shí),根據(jù)理想氣體多變狀態(tài)的規(guī)律,可得

    P_0 V_0=P_1 V_1 。 (10)

    考慮當(dāng)車輛左、右輪同時(shí)反向跳動(dòng)時(shí)所激起的車輛側(cè)傾模態(tài)振動(dòng),引起車身的側(cè)傾角為φ。此時(shí)A回路蓄能器中的油液將流進(jìn)左側(cè)液壓缸上腔和右側(cè)液壓缸下腔,B回路的右側(cè)液壓缸上腔和左側(cè)液壓缸下腔的油液流進(jìn)蓄能器中,引起液壓回路A和B的壓力變化。根據(jù)理想氣體多變狀態(tài)的規(guī)律,可以得到

    P_A=P_1 〖V_1〗^γ/〖V_A〗^γ 。 (11)

    P_B=P_1 〖V_1〗^γ/〖V_B〗^γ 。 (12)

    式中,γ為氣體多變指數(shù),取值為1.4;PA、PB分別為回路A、B的液體壓力;VA、VB分別為回路A、B的蓄能器氣體體積;V1為蓄能器的初始工作體積。

    根據(jù)車身姿態(tài)變化可以得出回路A、B中油液的體積變化?V_A、?V_B為

    ?V_A=-[(A_T1+A_B1)l_f+(A_T2+A_B2)l_r]sinφ。 (13)

    ?V_B=[(A_T1+A_B1)l_f+(A_T2+A_B2)l_r]sinφ。 (14)

    式中,D,d分別為液壓缸內(nèi)徑和活塞桿直徑。

    從而得到回路A、B中的蓄能器氣體體積VA和VB為

    V_A=V_1-?V_A 。 (15)

    V_B=V_1-?V_B 。 (16)

    則車身穩(wěn)定系統(tǒng)提供的抗側(cè)傾力矩M_φ為

    M_φ=(P_B-P_A )[(A_T1+A_B1 ) l_f+(A_T2+A_B2 ) l_r] 。(17)

    定義,λ= A_T2?A_T1 , λ_1=A_B1?A_T1 , λ_2=A_B2?A_T2 , α_1=P_1 A_T1 [(1+λ_1 ) l_f+(1+λ_2)λl_r],

    K_φ=(dM_φ)?dφ,可以得到提供的側(cè)傾剛度K_φ與側(cè)傾角φ之間的關(guān)系如下

    K_φ=γP_0^γ V_0^γ α_1^2 cosφ[1?〖(P_0 V_0+α_1 sinφ)〗^(γ+1) +( 1)?〖(P_0 V_0-α_1 sinφ)〗^(γ+1) ]。 (18)

    當(dāng)車身發(fā)生垂向位移Z_V時(shí),根據(jù)車身姿態(tài)變化可以得出回路A、B中油液的體積變化?V_A、?V_B為

    ?V_A=?V_B=(-A_T1+A_B1-A_T2+A_B2)Z_V 。 (19)

    則車身穩(wěn)定系統(tǒng)提供的抗垂向振動(dòng)的力F為

    F=(P_A+P_B)(A_B1-A_T1+A_B2-A_T2) 。 (20)

    定義 ,K_V=dF?(dZ_V ),α_2=P_1 A_T1 [(1-λ_1)+λ(1-λ_2)] , 可以得到垂向剛度K_V和垂向位移Z_V的關(guān)系如下

    K_V=(2γP_0^γ V_0^γ α_2^2)?〖(P_0 V_0-α_2 Z_V)〗^(γ+1) 。 (21)

    取車身穩(wěn)定系統(tǒng)設(shè)計(jì)參數(shù)初始值,系統(tǒng)工作壓力分別為1.5 MPa、2 MPa、2.5 MPa時(shí),仿真得到的車身穩(wěn)定系統(tǒng)所提供的側(cè)傾剛度增加的百分比隨車輛側(cè)傾角的變化如圖4所示,車身穩(wěn)定系統(tǒng)提供垂向剛度增加的百分比隨車身垂向位移的變化曲線如圖5所示。

    圖4車身穩(wěn)定系統(tǒng)提供側(cè)傾剛度增加的百分比

    圖5車身穩(wěn)定系統(tǒng)提供垂向剛度增加的百分比

    安裝車身穩(wěn)定系統(tǒng)之后,拆掉原車橫向穩(wěn)定桿,車輛側(cè)傾剛度比原車的側(cè)傾剛度增加最少為40%,即?K_φ≥40%。車輛平順性相比原車不能夠變差,即?K_V=0。

    取所增加的側(cè)傾剛度誤差為e(?K_φ)=0.02,取所增加的垂向剛度誤差為e(?K_V)=0.02。

    得到可用的參數(shù)組合見(jiàn)表3。

    表3車身穩(wěn)定系統(tǒng)參數(shù)組合

    序號(hào) 液壓缸內(nèi)徑 活塞桿直徑 蓄能器體積 蓄能器預(yù)充壓力 系統(tǒng)工作壓力

    D/m d/m V0/L P0/MPa P1/MPa

    1 0.042 0.014 0.4 1.3 2.1

    2 0.044 0.014 0.5 1.1 2

    3 0.046 0.014 0.7 0.9 2

    4 0.048 0.014 0.6 1.1 1.8

    5 0.048 0.016 0.7 1.1 2

    6 0.050 0.016 0.8 1.3 2.1

    2.2.2車身動(dòng)態(tài)穩(wěn)定系統(tǒng)參數(shù)匹配評(píng)估

    根據(jù)牛頓第二定律得到車輛機(jī)械系統(tǒng)動(dòng)力學(xué)方程

    mz ?_M (t)+cz ?_M (t)+kz_M (t)=F_H+F_M。 (22)

    式中,m,c,k分別為車輛的質(zhì)量、阻尼和剛度系數(shù);z_M (t),z ?_M (t),z ?_M (t)分別為位移、速度和加速度;F_H,F(xiàn)_M分別為液壓系統(tǒng)和外界的作用力。

    液壓缸為機(jī)械系統(tǒng)與液壓系統(tǒng)的邊界條件,根據(jù)流體基本公式可以得到無(wú)桿腔和有桿腔的流體壓縮量為

    Q_(c(T))=Z ?_(s-u) (t) A_T-Q_T (t)=(V_T P ?_T (t))?β 。 (23)

    Q_(c(B))=Q_B (t)-Z ?_(s-u) (t) A_B=(V_B P ?_B (t))?β。 (24)

    式中,Q_T (t),Q_B (t)分別為無(wú)桿腔和有桿腔的流量變化;VT,VB分別為液壓缸無(wú)桿腔和有桿腔體積變化;Z ?_(s-u)為懸架位移的變化量;P ?_T (t),P ?_B (t)分別為液壓缸無(wú)桿腔和有桿腔的壓力變化量;β為液壓介質(zhì)的體積彈性模量。

    液壓子系統(tǒng)的動(dòng)力學(xué)方程建立過(guò)程中,蓄能器的壓力變化量與流量和壓力的關(guān)系為

    P ?_a=(γQ_a P_a 〖(P_a?P_1 )〗^(1?γ))?V_1 。 (25)

    式中,P_a分別為蓄能器的瞬態(tài)工作壓力;Qa為蓄能器中的流量。

    將液壓管路分成多段進(jìn)行分析,每一段的動(dòng)力學(xué)方程為

    Q ?_i=(A_i (P_i1-P_i2))?(ρl_i ) 。 (26)

    式中,Ai管路的截面面積;Pi1,Pi2分別為其中某一段管路的兩端壓力;ρ為液壓介質(zhì)密度;li為液壓管路的長(zhǎng)度。

    引入狀態(tài)變量

    z=[z_s θ φ z_u1 z_u2 z_u3 z_u4

    z ?_s θ ? φ ? z ?_u1 z ?_u2 z ?_u3 z ?_u4

    P_T^C1 P_B^C1 P_T^C2 P_B^C2 P_T^C3 P_B^C3 P_T^C4 P_B^C4 P_a^1 P_a^2

    Q_T^C1 Q_B^C1 Q_T^C2 Q_B^C2 Q_T^C3 Q_B^C3 Q_T^C4 Q_B^C4 ]^T

    式中,zs,θ,φ,zu1,zu2,zu3,zu4分別為車身垂向位移,俯仰角,側(cè)傾角,4個(gè)車輪的垂向位移;P_T^Ci,P_B^Ci (i=1, 2, 3, 4)分別為液壓缸無(wú)桿腔和有桿腔的壓力;P_a^i (i=1, 2)為蓄能器的壓力;Q_T^Ci,Q_B^Ci(i=1,2,3,4)分別為液壓缸無(wú)桿腔和有桿腔的流量。從而得到狀態(tài)方程

    z ?=Az+B 。 (27)

    式中,

    A=[■(■( 0_(7×7) I_(7×7) 0_(7×8) 0_(7×2) 0_(7×8) @ -M^(-1) K -M^(-1) C 〖A1〗_(7×8) 0_(7×2) 0_(7×8) )@■( 0_(8×7) 〖A2〗_(8×7) 0_(8×8) 0_(8×2) 〖A3〗_(8×8) @ 0_(2×7) 0_(2×7) 0_(2×8) 0_(2×2) 〖A4〗_(2×8) )@ 0_(8×7) 0_(8×7) 〖A5〗_(8×8) 〖A6〗_(8×2) 0_(8×8) )]

    式中,A1為機(jī)械系統(tǒng)加速度與液壓系統(tǒng)的液壓缸進(jìn)出口壓力之間的關(guān)系矩陣;A2為液壓系統(tǒng)液壓缸進(jìn)出口壓力的變化量與機(jī)械系統(tǒng)的速度之間的關(guān)系矩陣;A3為液壓系統(tǒng)液壓缸進(jìn)出口壓力的變化量與液壓系統(tǒng)液壓缸進(jìn)出口流量之間的關(guān)系矩陣;A4液壓系統(tǒng)蓄能器壓力的變化量與液壓缸進(jìn)出口流量之間的關(guān)系矩陣;A5液壓系統(tǒng)液壓缸進(jìn)出口流量的變化量與液壓系統(tǒng)液壓缸進(jìn)出口壓力之間的關(guān)系矩陣;A6液壓系統(tǒng)液壓缸進(jìn)出口流量的變化量與蓄能器的壓力之間的關(guān)系矩陣。狀態(tài)矩陣A由液壓缸大小、蓄能器體積、蓄能器預(yù)充壓力以及系統(tǒng)工作壓力等參數(shù)決定。

    B=[0_(1×7) (4F_0)?m_S (2F_0 (b-a))?I_yy (m_s a_y (t))?I_xx

    ((k_t1 z_g1+F_0 ))?m_u1 ((k_t2 z_g2+F_0 ))?m_u2

    ((k_t3 z_g3+F_0))?m_u3 ((k_t4 z_g4+F_0))?m_u4 0_(1×18) ]^T。(28)

    式中,F(xiàn)0車身穩(wěn)定系統(tǒng)液壓缸對(duì)車身的初始作用力;ay(t)為車身側(cè)向加速度;kti(i=1, 2, 3, 4)為輪胎剛度;zgi(i=1, 2, 3, 4)為路面輸入激勵(lì),mui(i=1, 2, 3, 4)為車輛簧下質(zhì)量。

    顯然,由式(25)確定的機(jī)械液壓的動(dòng)力學(xué)方程具有時(shí)域特性,通過(guò)給定的車身側(cè)傾角和垂向加速度目標(biāo),對(duì)車身穩(wěn)定系統(tǒng)設(shè)計(jì)參數(shù)進(jìn)行優(yōu)化評(píng)估。

    2.3車身動(dòng)態(tài)穩(wěn)定系統(tǒng)參數(shù)匹配結(jié)果

    整車操縱穩(wěn)定性和平順性仿真測(cè)試中,車輛操縱穩(wěn)定性仿真試驗(yàn)輸入激勵(lì)為原越野車蛇行試驗(yàn)時(shí)所采集的加速度如圖6所示,仿真中車身穩(wěn)定系統(tǒng)的設(shè)計(jì)參數(shù)設(shè)置見(jiàn)表3,得到不同組合的車輛側(cè)傾角的時(shí)間歷程曲線如圖8所示。其中安裝車身穩(wěn)定系統(tǒng)之后車身側(cè)傾角雖然有輕微波動(dòng)但是獲得的車身側(cè)傾角皆小于裝有橫向穩(wěn)定桿車輛。

    車輛平順性仿真試驗(yàn)輸入激勵(lì)為C級(jí)路面隨機(jī)輸入激勵(lì)如圖7所示,仿真中車身穩(wěn)定系統(tǒng)的設(shè)置參數(shù)見(jiàn)表3,得到不同參數(shù)組合的車輛垂向加速度時(shí)間歷程曲線如圖9所示。

    圖6車輛操穩(wěn)性仿真測(cè)試側(cè)向加速度輸入

    圖7車輛平順性仿真測(cè)試路面隨機(jī)輸入

    圖8整車側(cè)傾角曲線

    圖9在路面隨機(jī)輸入下質(zhì)心處加速度響應(yīng)

    根據(jù)車輛懸架性能試驗(yàn)和整車性能試驗(yàn)評(píng)價(jià)標(biāo)準(zhǔn)和試驗(yàn)方法進(jìn)行仿真測(cè)試,以此來(lái)驗(yàn)證車身穩(wěn)定系統(tǒng)參數(shù)匹配的有效性和可行性。綜合考慮懸架設(shè)計(jì)中操縱穩(wěn)定性和平順性的性能指標(biāo),以及液壓系統(tǒng)設(shè)計(jì)標(biāo)準(zhǔn)和經(jīng)驗(yàn),最終得到仿真的整車操縱穩(wěn)定性和平順性的結(jié)果如圖10和圖11所示。

    圖10車速50 km/h整車側(cè)傾角曲線

    從汽車操縱穩(wěn)定性試驗(yàn)仿真結(jié)果圖10可以看出,比較車輛的側(cè)傾角大小,裝有車身穩(wěn)定系統(tǒng)車輛的側(cè)傾角要比裝有橫向穩(wěn)定桿車輛的要小,車身穩(wěn)定系統(tǒng)匹配的設(shè)計(jì)參數(shù)使車輛的側(cè)傾角大小在整車性能的目標(biāo)范圍內(nèi)。

    圖11在路面隨機(jī)輸入下質(zhì)心處加速度響應(yīng)

    從汽車平順性試驗(yàn)仿真結(jié)果圖11可以看出,裝有車身穩(wěn)定系統(tǒng)車輛的加速度響應(yīng)相比于原車的差值非常小,車身穩(wěn)定系統(tǒng)匹配的設(shè)計(jì)參數(shù)使車輛的質(zhì)心垂向加速度響應(yīng)大小在整車性能的目標(biāo)范圍內(nèi)。

    根據(jù)液壓缸和蓄能器的設(shè)計(jì)標(biāo)準(zhǔn),結(jié)合仿真試驗(yàn)分析,最終確定匹配的車身穩(wěn)定系統(tǒng)參數(shù)見(jiàn)表4。

    表4車身穩(wěn)定系統(tǒng)部分參數(shù)

    參數(shù)

    液壓缸內(nèi)徑 D/m 0.04

    活塞桿直徑 d/m 0.018

    蓄能器體積 V0/L 0.5

    蓄能器預(yù)充壓力 P0/MPa 1

    系統(tǒng)工作壓力/MPa P1 2

    2.4車身穩(wěn)定系統(tǒng)與整車匹配結(jié)構(gòu)的工程實(shí)現(xiàn)

    車輛在粗糙不平的壞路上行駛時(shí),由于不平路面的激勵(lì)造成車輪上下劇烈跳動(dòng),嚴(yán)重時(shí)會(huì)出現(xiàn)車軸撞擊緩沖塊的現(xiàn)象。如果車身穩(wěn)定系統(tǒng)的液壓缸行程設(shè)計(jì)不合理就會(huì)導(dǎo)致在遇到惡劣路面時(shí),液壓缸活塞撞擊液壓缸頂部端蓋現(xiàn)象,對(duì)于懸架構(gòu)件和車架的強(qiáng)度、疲勞壽命等造成影響[1]。

    運(yùn)用CATIA分析越野車后懸鋼板彈簧變形運(yùn)動(dòng),以此進(jìn)行車身穩(wěn)定系統(tǒng)與整車懸架系統(tǒng)的匹配設(shè)計(jì)。使得車身穩(wěn)定系統(tǒng)液壓缸行程在原車懸架系統(tǒng)的運(yùn)動(dòng)行程范圍之內(nèi)。分析得到板簧變形運(yùn)動(dòng)的結(jié)果如圖12所示,從圖中可以得到后懸架板簧動(dòng)載軸管中心到車架下表面最小距離為149.3 mm,板簧動(dòng)載軸管中心到車架下表面最大距離為329.4 mm。

    圖12 板簧變形運(yùn)動(dòng)分析

    支撐構(gòu)件設(shè)計(jì)需滿足強(qiáng)度要求,圖13(a)為所設(shè)計(jì)的液壓缸安裝支架。

    在強(qiáng)度分析的框架中,需要確定負(fù)荷作用,材料參數(shù)以及邊界條件,運(yùn)用有限元方法(Finite Element Method,F(xiàn)EM)將拉、壓等載荷作用下的局部、整體范圍內(nèi)的最大值與材料參數(shù)進(jìn)行比較,得出失效可能性以及材料是否充分應(yīng)用和優(yōu)化潛力的結(jié)論[1]。通過(guò)UG三維建模,二維CAD加工圖紙?jiān)O(shè)計(jì)以及有限元分析,圖13(b)所示為液壓缸安裝支架應(yīng)力分布云圖。

    (a) 安裝支架三維建模圖

    (b)安裝支架應(yīng)力分布云圖

    圖13液壓缸安裝支架

    有限元仿真分析拉、壓載荷情況下,安裝支架FEM仿真結(jié)果如圖13(b)所示,安裝支架節(jié)點(diǎn)位移和單元節(jié)點(diǎn)的最大值V_max,最小值V_min結(jié)果匯總見(jiàn)表5和表6。

    表5安裝支架節(jié)點(diǎn)位移仿真結(jié)果

    位移-節(jié)點(diǎn)VEL /mm

    X Y Z 幅值

    V_max 3.35e-003 1.76e-003 9.93e-004 1.38e-002

    V_min -4.21e003 -2.40e003 -1.30e-002 0

    表6安裝支架單元節(jié)點(diǎn)仿真結(jié)果

    應(yīng)力-單元節(jié)點(diǎn)σ/kPa

    Von

    Mises 最小

    主應(yīng)力 最大

    主應(yīng)力 最大剪切

    V_max 2.54e+004 2.35e+003 2.30e+004 1.45e+004

    V_min 2.47e-007 -2.57e+004 -3.8e+003 1.41e-007

    仿真結(jié)果表明液壓缸安裝支架在所選用材料范圍內(nèi)滿足材料的性能要求。

    根據(jù)工程設(shè)計(jì)方法與經(jīng)驗(yàn),車身穩(wěn)定系統(tǒng)與車輛的匹配設(shè)計(jì)結(jié)果如圖14所示。橡膠金屬連接件對(duì)于固體聲的傳播非常重要,橡膠金屬對(duì)于在聲源與聲的接受者(人體)之間的具有隔斷作用[1]。越野車車身穩(wěn)定系統(tǒng)工程實(shí)際安裝中,在液壓缸安裝支架與車身之間增加橡膠金屬連接件如圖14(b)所示。

    (a) 蓄能器安裝圖 (b) 液壓缸安裝圖

    圖14 車身穩(wěn)定系統(tǒng)與整車匹配完成局部圖

    3試驗(yàn)結(jié)果及分析

    為驗(yàn)證車身穩(wěn)定系統(tǒng)的有效性和可行性,根據(jù)車輛懸架性能試驗(yàn)和整車性能試驗(yàn)評(píng)價(jià)標(biāo)準(zhǔn)和試驗(yàn)方法進(jìn)行測(cè)試。依據(jù)車身穩(wěn)定系統(tǒng)與整車的匹配性設(shè)計(jì)與研究結(jié)果,在越野車上安裝車身穩(wěn)定系統(tǒng)(HIS),取代原車橫向穩(wěn)定桿(ARB),得到改裝后的試驗(yàn)樣車并進(jìn)行實(shí)車性能測(cè)試試驗(yàn),驗(yàn)證車身穩(wěn)定系統(tǒng)在改善車輛操縱穩(wěn)定性上的有效性。試驗(yàn)數(shù)據(jù)分析結(jié)果見(jiàn)表7~9,試驗(yàn)曲線見(jiàn)圖15~17。

    表7 懸架性能實(shí)驗(yàn)結(jié)果與分析

    HIS ARB E

    ω_nroll /Hz 1.657 1.373 20.68%

    Κ_roll /(kN﹒m/rad) 193.3 132.7 45.65%

    圖15 車身側(cè)傾角頻率響應(yīng)

    從懸架性能試驗(yàn)結(jié)果表7和圖15可以看出,安裝車身穩(wěn)定系統(tǒng)后,車輛的側(cè)傾固有頻率比原車的側(cè)傾固有頻率增加的百分比E為20.68%,車輛側(cè)傾剛度比原車整車側(cè)傾剛度增加45.65%。說(shuō)明車身穩(wěn)定系統(tǒng)能夠有效提高車輛的側(cè)傾穩(wěn)定性。

    表8蛇行試驗(yàn)結(jié)果與分析

    HIS ARB E

    φ_max /( ) 2.984 3.510 14.99%

    /( ) 1.797 4 2.834 0 36.58%

    圖16 蛇行實(shí)驗(yàn)—車速50km/h整車側(cè)傾角曲線

    從蛇行試驗(yàn)結(jié)果表8和圖16可以看出,比較裝有車身穩(wěn)定系統(tǒng)和橫向穩(wěn)定桿車輛的側(cè)傾角大小,可以看出安裝車身穩(wěn)定系統(tǒng)后,相同的試驗(yàn)工況下,車速為50 km/h時(shí),最大側(cè)傾角減小了14.99%,平均側(cè)傾角減小了36.58%。說(shuō)明車身穩(wěn)定系統(tǒng)能夠有效降低車身側(cè)傾角,提高車輛的側(cè)傾穩(wěn)定性,改善乘員的舒適性和安全性。

    表9 平順性隨機(jī)輸入行駛試驗(yàn)結(jié)果分析

    試驗(yàn)

    測(cè)量參數(shù) 結(jié)果分析:各測(cè)點(diǎn)的加速度均方根值

    差值分析E

    速度v/(km·h-1) 30 40 50

    駕駛員座椅處 0.002 3 -0.006 4 -0.027 4

    質(zhì)心 0.000 9 0.013 2 -0.006 5

    后右座椅處 0.004 5 0.026 8 0.005 7

    后左座椅處 0.056 5 0.095 9 0.042 2

    圖17在路面隨機(jī)輸入下質(zhì)心處垂向加速度響應(yīng)

    從汽車平順性試驗(yàn)結(jié)果表9和圖17可以看出,裝有車身穩(wěn)定系統(tǒng)車輛的加速度均方根值相比于原車的差值非常小,不超過(guò)0.1。說(shuō)明車身穩(wěn)定系統(tǒng)對(duì)平順性的影響非常小。

    4結(jié)論

    本文對(duì)某越野車所設(shè)計(jì)的車身穩(wěn)定系統(tǒng),考慮了該車的底盤結(jié)構(gòu),進(jìn)行了強(qiáng)度分析。同時(shí)對(duì)改裝后的越野車進(jìn)行了懸架及整車性能試驗(yàn),從試驗(yàn)結(jié)果中可以得到如下結(jié)論:

    車身穩(wěn)定系統(tǒng)關(guān)鍵物理參數(shù)對(duì)車輛性能的影響進(jìn)行優(yōu)化分析,得到了與整車匹配的設(shè)計(jì)參數(shù)。匹配設(shè)計(jì)兼顧了車輛懸架設(shè)計(jì)中平順性和操縱穩(wěn)定性,取得了較好的實(shí)際應(yīng)用效果。

    安裝車身穩(wěn)定系統(tǒng)的該越野車的懸架性能得到優(yōu)化,車輛側(cè)傾剛度增加。

    安裝車身穩(wěn)定系統(tǒng)的車輛操縱穩(wěn)定性得到提高,車身側(cè)傾角減小,乘員的乘坐舒適性和安全性得到改善,性能優(yōu)于通常所采用的橫向穩(wěn)定桿。

    車身穩(wěn)定系統(tǒng)在兼顧車輛操縱穩(wěn)定性和行駛平順性方面的有效性和可行性得到試驗(yàn)驗(yàn)證。

    參考文獻(xiàn)

    [1] 日本自動(dòng)車技術(shù)會(huì)(中國(guó)汽車工程學(xué)會(huì)組譯). 汽車工程手冊(cè). 5底盤設(shè)計(jì)篇 [M]. 北京:北京理工大學(xué)出版社,2010: 31-33.

    Society of Automotive Engineering of Japan, Society of Automotive Engineering of China. Automotive Technology Handbook 5.Chassis Design [M]. Beijing: Beijing Institute of Technology Press, 2010: 31-33. (in Chinese)

    [2] 劉偉,史文庫(kù),方德廣,等. 針對(duì)汽車操縱穩(wěn)定性的懸架系統(tǒng)控制的研究 [J]. 汽車工程, 2011, 33(6): 502-511.

    Liu Wei, Shi Wenku, Fang Deguang, et al. A Research on Suspension System Control for Vehicle Handling and Stability [J]. Automotive Engineering, 2011, 33(6): 502-511. (in Chinese)

    [3] 劉偉,史文庫(kù),方德廣,等. 汽車動(dòng)力學(xué)分析及懸架子系統(tǒng)優(yōu)化設(shè)計(jì) [J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào), 2012, 44(3): 96-101.

    Liu Wei, Shi Wenku, Fang Deguang, et al. Vehicle Dynamics Analysis and Optimization of Suspension Design [J]. Journal of Harbin Institute of Technology, 2012, 44(3): 96-101. (in Chinese)

    [4] 劉偉,史文庫(kù),桂龍明,等. 基于平順性與操縱穩(wěn)定性的懸架系統(tǒng)多目標(biāo)優(yōu)化 [J]. 吉林大學(xué)學(xué)報(bào), 2011, 41(5): 1199-1204 .

    Liu Wei, Shi Wenku, Gui Longming, et al. Multi-objective Optimization of Suspension System Based on Vehicle Ride Comfort and Handling Stability [J]. Journal of Jilin University, 2011, 41(5): 1199-1204. (in Chinese)

    [5] 黃東明. 某越野車動(dòng)力學(xué)仿真研究 [D]. 南京: 南京理工大學(xué), 2011.

    Huang Dongming. Off-road Vehicle Dynamic Simulation Research [D]. Nanjing: Nanjing University, 2011. (in Chinese)

    [6] 蘇小平. 依維柯汽車多體動(dòng)力學(xué)仿真分析、優(yōu)化研究及工程實(shí)現(xiàn) [D]. 南京: 南京理工大學(xué), 2004.

    Su Xiaoping. The Muti-body Dynamic Simulation Analysis, Optimization Design and Engineering Application of IVECO Vehicle [D]. Nanjing: Nanjing University, 2004. (in Chinese)

    [7] ZHANG N, SMITH W, JEYAKUMARAN J. Hydraulically Interconnected Vehicle Suspension: Background and Modeling [J],.Vehicle System Dynamics, 2010, 48(1): 17-40.

    [8] 丁飛,張農(nóng),韓旭. 安裝液壓互聯(lián)懸架貨車的機(jī)械液壓多體系統(tǒng)建模及模態(tài)分析 [J]. 機(jī)械工程學(xué)報(bào), 2012, 48(6): 116-123.

    Ding Fei, Zhang Nong, Han Xu. Modeling and Modal Analysis of Multi-body Truck System Fitted with Hydraulically Interconnected Suspension [J]. Journal of Mechanical Engineering, 2012, 48(6): 116-123. (in Chinese)

    [9] SMITH W, ZHANG N, JEYAKUMARAN J. Hydraulically Interconnected Vehicle Suspension: Theoretical and Experimental Ride Analysis [J]. Vehicle System Dynamics, 2010, 48(1): 41-64.

    [10] SMITH W, ZHANG N, HU W. Hydraulically Interconnected Vehicle Suspension: Handling Performance [J]. Vehicle System Dynamics, 2010,48(1):1-20.

    [11] ZHANG N, DONG G, DU H P. Investigation into Unstripped Rollover of Light Vehicles in the Modified Fishhook and the Sine Maneuvers Part I: Vehicle Modelling Roll and Yaw Instability [J]. Vehicle System Dynamics, 2008, 46(4): 271-293.

    [12] DONG G, ZHANG N, DU Hai-Ping. Investigation into Unstripped Rollover of Light Vehicles in the Modified Fishhook and the Sine Maneuvers PartⅡ: Vehicle Modelling, Roll and Yaw Instability [J]. Vehicle System Dynamics, 2011, 49(6): 949-968.

    [13]WANG L, ZHANG N, DU H. Modeling, Parameter Estimation and Testing of a Vehicle with Anti-roll Systems [C]//. 6th Australasian Congress on Applied Mechanics, Perth, Australia,2010.

    作者簡(jiǎn)介:

    責(zé)任作者:方德廣(1965—),男,江蘇射陽(yáng)人。博士,高級(jí)工程師,主要從事整車總布置及匹配的研究。

    電話:13505177248

    郵箱:fangdeguang2013@126.com

    通訊作者:張雨(1990—),女,湖北隨州人。碩士研究生,主要研究方向?yàn)閰?shù)識(shí)別及液壓互聯(lián)懸架。

    電話:15211096747

    郵箱:zhangyu_jane@126.com

    猜你喜歡
    試驗(yàn)設(shè)計(jì)
    何為設(shè)計(jì)的守護(hù)之道?
    《豐收的喜悅展示設(shè)計(jì)》
    流行色(2020年1期)2020-04-28 11:16:38
    瞞天過(guò)?!律O(shè)計(jì)萌到家
    CS95
    世界汽車(2017年8期)2017-08-12 04:39:15
    510
    世界汽車(2017年8期)2017-08-12 04:26:42
    馭勝S330
    世界汽車(2017年8期)2017-08-12 04:17:18
    設(shè)計(jì)秀
    海峽姐妹(2017年7期)2017-07-31 19:08:17
    有種設(shè)計(jì)叫而專
    Coco薇(2017年5期)2017-06-05 08:53:16
    C-NCAP 2016年第八號(hào)試驗(yàn)發(fā)布
    汽車與安全(2016年5期)2016-12-01 05:22:16
    試驗(yàn)
    太空探索(2016年12期)2016-07-18 11:13:43
    色综合色国产| 国产一区二区三区av在线 | 国产v大片淫在线免费观看| 一进一出好大好爽视频| 天天躁夜夜躁狠狠久久av| av天堂中文字幕网| 成人高潮视频无遮挡免费网站| 97超级碰碰碰精品色视频在线观看| 久久午夜福利片| 男女下面进入的视频免费午夜| 国产精品亚洲一级av第二区| 人妻制服诱惑在线中文字幕| 亚洲,欧美,日韩| 尾随美女入室| 伦精品一区二区三区| 精品少妇黑人巨大在线播放 | 色视频www国产| 亚洲一区高清亚洲精品| 最近的中文字幕免费完整| 国产精品美女特级片免费视频播放器| 在线免费观看不下载黄p国产| 麻豆精品久久久久久蜜桃| 啦啦啦啦在线视频资源| 在线免费十八禁| 99九九线精品视频在线观看视频| 噜噜噜噜噜久久久久久91| 亚洲自偷自拍三级| 国产精品日韩av在线免费观看| 亚洲av美国av| 日韩中字成人| 日韩av不卡免费在线播放| 我的老师免费观看完整版| 在线观看66精品国产| 少妇熟女aⅴ在线视频| 丝袜喷水一区| 久久精品国产亚洲av天美| av女优亚洲男人天堂| 简卡轻食公司| 国产精品久久久久久久久免| 亚洲图色成人| 中文字幕人妻熟人妻熟丝袜美| 99热全是精品| 99热全是精品| 久久精品国产亚洲av天美| 桃色一区二区三区在线观看| 一区二区三区免费毛片| 精品午夜福利在线看| 亚洲国产欧美人成| 久久久精品94久久精品| av中文乱码字幕在线| 大型黄色视频在线免费观看| av在线观看视频网站免费| 久久精品人妻少妇| 一个人看的www免费观看视频| 人妻丰满熟妇av一区二区三区| 欧美另类亚洲清纯唯美| 免费在线观看成人毛片| 18禁裸乳无遮挡免费网站照片| 日本免费a在线| 国产一区亚洲一区在线观看| 色哟哟·www| 欧美日韩乱码在线| 色播亚洲综合网| 欧美性感艳星| 欧美精品国产亚洲| 国产私拍福利视频在线观看| 日韩精品有码人妻一区| av天堂在线播放| 又粗又爽又猛毛片免费看| 国产免费一级a男人的天堂| av.在线天堂| 日韩国内少妇激情av| 成年版毛片免费区| 日本成人三级电影网站| 成人高潮视频无遮挡免费网站| 在线国产一区二区在线| 小蜜桃在线观看免费完整版高清| 国产成人一区二区在线| 久久久久国产精品人妻aⅴ院| 少妇的逼水好多| 日本一本二区三区精品| 免费黄网站久久成人精品| 99热6这里只有精品| 国产片特级美女逼逼视频| 欧美日韩乱码在线| 美女内射精品一级片tv| 欧美成人a在线观看| 最近最新中文字幕大全电影3| 欧美一区二区精品小视频在线| 亚洲人成网站在线播放欧美日韩| 精品一区二区三区视频在线观看免费| 免费观看在线日韩| 国产午夜福利久久久久久| 日本黄大片高清| 亚洲第一区二区三区不卡| 在线天堂最新版资源| 精品不卡国产一区二区三区| 国产精品亚洲一级av第二区| 一级a爱片免费观看的视频| АⅤ资源中文在线天堂| 欧美在线一区亚洲| 欧美xxxx性猛交bbbb| 99热这里只有是精品50| 国产探花在线观看一区二区| 啦啦啦韩国在线观看视频| 国产久久久一区二区三区| 日本三级黄在线观看| 日本三级黄在线观看| 国产91av在线免费观看| 欧美成人精品欧美一级黄| 精品人妻一区二区三区麻豆 | 亚洲美女搞黄在线观看 | 18禁黄网站禁片免费观看直播| 久久精品国产自在天天线| 丰满的人妻完整版| 久久精品影院6| 国产精品女同一区二区软件| 黄色视频,在线免费观看| 超碰av人人做人人爽久久| 亚洲av熟女| 女人十人毛片免费观看3o分钟| 欧美性感艳星| 久久久久国内视频| 女人十人毛片免费观看3o分钟| 一a级毛片在线观看| av天堂在线播放| 成人毛片a级毛片在线播放| 国产精品99久久久久久久久| 国产乱人偷精品视频| 干丝袜人妻中文字幕| 亚洲av电影不卡..在线观看| 直男gayav资源| 啦啦啦观看免费观看视频高清| 日韩欧美三级三区| 久久亚洲精品不卡| 97超视频在线观看视频| 日本a在线网址| 免费搜索国产男女视频| 久久天躁狠狠躁夜夜2o2o| 国产成人影院久久av| 国产精品一二三区在线看| 亚洲国产精品久久男人天堂| 中文资源天堂在线| av在线蜜桃| 精品人妻视频免费看| 日本 av在线| 久久久精品94久久精品| 狠狠狠狠99中文字幕| 久久国内精品自在自线图片| 国产人妻一区二区三区在| aaaaa片日本免费| 国产乱人偷精品视频| 日本在线视频免费播放| 午夜激情欧美在线| 国产精品永久免费网站| 国产成人a区在线观看| 青春草视频在线免费观看| 97超碰精品成人国产| 精品福利观看| 老熟妇仑乱视频hdxx| 深夜a级毛片| 久久婷婷人人爽人人干人人爱| 精品午夜福利在线看| 天美传媒精品一区二区| 国产精华一区二区三区| 欧美绝顶高潮抽搐喷水| 午夜日韩欧美国产| 尾随美女入室| 不卡视频在线观看欧美| av在线蜜桃| 国内精品美女久久久久久| 91久久精品电影网| 国产一区二区在线观看日韩| 成人国产麻豆网| 亚洲18禁久久av| 校园人妻丝袜中文字幕| 高清毛片免费观看视频网站| 国产成人aa在线观看| 久久人人爽人人片av| 老女人水多毛片| 亚洲人成网站在线观看播放| 日韩精品中文字幕看吧| 日韩欧美 国产精品| 美女免费视频网站| 国产不卡一卡二| 久久精品综合一区二区三区| 美女被艹到高潮喷水动态| 国产精品99久久久久久久久| 免费在线观看影片大全网站| www.色视频.com| 精品乱码久久久久久99久播| 我要搜黄色片| 一本一本综合久久| 老司机福利观看| 免费人成在线观看视频色| 国产午夜精品论理片| 国产视频一区二区在线看| 亚洲欧美日韩无卡精品| 十八禁网站免费在线| 男插女下体视频免费在线播放| 精品人妻一区二区三区麻豆 | 国产一区二区在线观看日韩| 在线观看午夜福利视频| 日本一本二区三区精品| 国产爱豆传媒在线观看| 最近视频中文字幕2019在线8| 午夜精品国产一区二区电影 | 亚洲精品456在线播放app| 色5月婷婷丁香| 无遮挡黄片免费观看| 亚洲av第一区精品v没综合| 精品久久久久久成人av| 婷婷亚洲欧美| 久久精品国产亚洲av香蕉五月| 午夜福利在线观看吧| 特大巨黑吊av在线直播| 黄色欧美视频在线观看| 亚洲人成网站在线观看播放| 国产亚洲精品久久久久久毛片| 国产三级在线视频| av免费在线看不卡| 欧美日韩在线观看h| 亚洲国产欧洲综合997久久,| 久久久久久久久大av| 欧美日韩乱码在线| 久久精品综合一区二区三区| 好男人在线观看高清免费视频| 一区二区三区免费毛片| 一级av片app| 免费在线观看成人毛片| 精品久久久噜噜| 真实男女啪啪啪动态图| 成人特级黄色片久久久久久久| 婷婷精品国产亚洲av| 五月伊人婷婷丁香| 五月伊人婷婷丁香| 岛国在线免费视频观看| 日本精品一区二区三区蜜桃| 女人十人毛片免费观看3o分钟| 色播亚洲综合网| 熟女人妻精品中文字幕| 国产中年淑女户外野战色| 久久草成人影院| 国产精品伦人一区二区| 免费人成在线观看视频色| 大又大粗又爽又黄少妇毛片口| 亚洲经典国产精华液单| 联通29元200g的流量卡| 精品一区二区三区视频在线观看免费| 你懂的网址亚洲精品在线观看 | 性插视频无遮挡在线免费观看| 狂野欧美白嫩少妇大欣赏| 亚洲一区二区三区色噜噜| 国产男靠女视频免费网站| 日日啪夜夜撸| 女生性感内裤真人,穿戴方法视频| 深夜精品福利| 国产亚洲精品久久久com| 国产成人影院久久av| 又黄又爽又刺激的免费视频.| 久久人妻av系列| 黄色日韩在线| eeuss影院久久| 99在线人妻在线中文字幕| 免费高清视频大片| 久久久精品94久久精品| 国产亚洲精品久久久com| 一边摸一边抽搐一进一小说| 亚洲av不卡在线观看| 亚洲精品乱码久久久v下载方式| 国产成人一区二区在线| 中文资源天堂在线| 午夜福利在线在线| 色在线成人网| 天堂√8在线中文| 女生性感内裤真人,穿戴方法视频| 亚洲四区av| 熟女人妻精品中文字幕| 国内精品久久久久精免费| 日韩大尺度精品在线看网址| 精品不卡国产一区二区三区| 三级国产精品欧美在线观看| 久久99热这里只有精品18| 国产精品国产三级国产av玫瑰| videossex国产| 男人和女人高潮做爰伦理| 亚洲av成人精品一区久久| 看片在线看免费视频| 中文字幕av成人在线电影| 国产私拍福利视频在线观看| 18禁在线无遮挡免费观看视频 | 日日摸夜夜添夜夜添av毛片| 亚洲欧美日韩东京热| 一进一出抽搐gif免费好疼| 精品人妻偷拍中文字幕| 国产男人的电影天堂91| 女生性感内裤真人,穿戴方法视频| 插阴视频在线观看视频| 国产精品美女特级片免费视频播放器| 一本精品99久久精品77| 天堂网av新在线| 亚洲第一区二区三区不卡| 97碰自拍视频| 天堂√8在线中文| 欧美一区二区亚洲| 久久精品夜色国产| 内地一区二区视频在线| 97超碰精品成人国产| 国产精品1区2区在线观看.| 日韩一本色道免费dvd| 97人妻精品一区二区三区麻豆| 亚洲精品色激情综合| 亚洲性夜色夜夜综合| 国产欧美日韩一区二区精品| 性欧美人与动物交配| 成人欧美大片| 大香蕉久久网| 91午夜精品亚洲一区二区三区| 色噜噜av男人的天堂激情| 国产国拍精品亚洲av在线观看| 国产成人a∨麻豆精品| 国产淫片久久久久久久久| 亚洲第一电影网av| a级一级毛片免费在线观看| 国产精品1区2区在线观看.| 99久久久亚洲精品蜜臀av| 日本 av在线| 国产精品三级大全| 国产亚洲av嫩草精品影院| 97超级碰碰碰精品色视频在线观看| 久久久久久久久久久丰满| 欧美潮喷喷水| 婷婷精品国产亚洲av| 中文字幕av成人在线电影| 久久精品久久久久久噜噜老黄 | 99久国产av精品| 别揉我奶头 嗯啊视频| 精品无人区乱码1区二区| 国产在视频线在精品| 99热这里只有是精品在线观看| 国产久久久一区二区三区| 99热6这里只有精品| 韩国av在线不卡| 最近在线观看免费完整版| 午夜福利视频1000在线观看| av视频在线观看入口| 欧美一区二区亚洲| 少妇的逼好多水| 国产高清视频在线观看网站| 看非洲黑人一级黄片| 午夜福利18| 国产女主播在线喷水免费视频网站 | 一级毛片久久久久久久久女| 亚洲五月天丁香| 俄罗斯特黄特色一大片| 国产亚洲欧美98| 真实男女啪啪啪动态图| 免费无遮挡裸体视频| 亚洲精品亚洲一区二区| 欧美一区二区精品小视频在线| 亚洲内射少妇av| 久99久视频精品免费| 午夜视频国产福利| 国产av在哪里看| 噜噜噜噜噜久久久久久91| 秋霞在线观看毛片| 亚洲无线观看免费| 深夜a级毛片| 99国产精品一区二区蜜桃av| 中出人妻视频一区二区| 搡老妇女老女人老熟妇| 亚洲国产欧美人成| 精品午夜福利视频在线观看一区| 两个人的视频大全免费| 日韩强制内射视频| 内地一区二区视频在线| 国产极品精品免费视频能看的| 婷婷色综合大香蕉| 美女免费视频网站| 69av精品久久久久久| 国产 一区 欧美 日韩| 97碰自拍视频| 六月丁香七月| 黄色欧美视频在线观看| 欧美日韩在线观看h| 在线观看午夜福利视频| 淫妇啪啪啪对白视频| ponron亚洲| 99九九线精品视频在线观看视频| 亚洲av不卡在线观看| 成人亚洲欧美一区二区av| 精品午夜福利在线看| 最近2019中文字幕mv第一页| 免费av观看视频| 亚洲美女黄片视频| 我要搜黄色片| 别揉我奶头 嗯啊视频| 日本-黄色视频高清免费观看| 露出奶头的视频| 麻豆乱淫一区二区| 日本熟妇午夜| 久久婷婷人人爽人人干人人爱| 免费观看的影片在线观看| 国产在视频线在精品| 91久久精品电影网| 91在线精品国自产拍蜜月| 亚洲丝袜综合中文字幕| 成人国产麻豆网| 狂野欧美白嫩少妇大欣赏| 插逼视频在线观看| 我要搜黄色片| 老师上课跳d突然被开到最大视频| 此物有八面人人有两片| 夜夜爽天天搞| 亚洲丝袜综合中文字幕| eeuss影院久久| 午夜老司机福利剧场| 久久亚洲国产成人精品v| 97超视频在线观看视频| 国产69精品久久久久777片| 亚洲电影在线观看av| 两个人视频免费观看高清| 国产精品野战在线观看| 日韩制服骚丝袜av| 午夜激情欧美在线| 亚洲电影在线观看av| 两个人视频免费观看高清| 亚洲人成网站在线观看播放| 五月伊人婷婷丁香| 最后的刺客免费高清国语| 日韩在线高清观看一区二区三区| 黄片wwwwww| 99在线视频只有这里精品首页| 国产免费男女视频| 高清日韩中文字幕在线| 天天躁日日操中文字幕| 亚洲人成网站高清观看| 国产精华一区二区三区| 国产精品女同一区二区软件| 欧美3d第一页| 国产老妇女一区| 偷拍熟女少妇极品色| 亚洲天堂国产精品一区在线| 国产精品乱码一区二三区的特点| 精品少妇黑人巨大在线播放 | 久久久色成人| 欧美+日韩+精品| 色综合亚洲欧美另类图片| 啦啦啦啦在线视频资源| 黄色日韩在线| 亚洲无线在线观看| 国内久久婷婷六月综合欲色啪| 日韩av不卡免费在线播放| 嫩草影院精品99| 久久久久精品国产欧美久久久| 亚洲激情五月婷婷啪啪| 欧美日韩乱码在线| 乱码一卡2卡4卡精品| 天天躁夜夜躁狠狠久久av| 欧美日韩乱码在线| av黄色大香蕉| 欧美+亚洲+日韩+国产| 日日摸夜夜添夜夜爱| 日本一本二区三区精品| 国产av在哪里看| 中国美女看黄片| 精品久久久久久久人妻蜜臀av| av国产免费在线观看| 久久久色成人| 亚洲专区国产一区二区| 偷拍熟女少妇极品色| 99热6这里只有精品| 国产熟女欧美一区二区| 乱人视频在线观看| 最近最新中文字幕大全电影3| av中文乱码字幕在线| 18禁裸乳无遮挡免费网站照片| 日韩在线高清观看一区二区三区| 搡女人真爽免费视频火全软件 | 亚洲av不卡在线观看| 日韩欧美 国产精品| 夜夜爽天天搞| 久久亚洲国产成人精品v| 国产精品,欧美在线| 欧美中文日本在线观看视频| 欧美一区二区精品小视频在线| 深爱激情五月婷婷| 免费黄网站久久成人精品| 日本与韩国留学比较| 免费搜索国产男女视频| 成人美女网站在线观看视频| 亚洲成人精品中文字幕电影| 91久久精品国产一区二区成人| 日韩大尺度精品在线看网址| 最好的美女福利视频网| 两个人视频免费观看高清| 在线国产一区二区在线| 插逼视频在线观看| 亚洲熟妇中文字幕五十中出| 亚洲色图av天堂| 国产毛片a区久久久久| 日韩三级伦理在线观看| 毛片女人毛片| 欧美日本视频| 成年女人永久免费观看视频| 成人二区视频| 亚洲中文字幕日韩| 最近视频中文字幕2019在线8| 欧美另类亚洲清纯唯美| 在线观看av片永久免费下载| 亚洲av熟女| 成人美女网站在线观看视频| 热99在线观看视频| 精品久久久久久久末码| 久久久精品大字幕| 久久精品国产清高在天天线| 精品人妻熟女av久视频| 欧美绝顶高潮抽搐喷水| 国产一区二区三区在线臀色熟女| 成熟少妇高潮喷水视频| 欧美三级亚洲精品| 舔av片在线| 久久久久精品国产欧美久久久| 亚洲成人精品中文字幕电影| 国产视频内射| 亚洲中文字幕一区二区三区有码在线看| 免费搜索国产男女视频| 99热这里只有精品一区| 中文字幕av在线有码专区| 亚洲熟妇中文字幕五十中出| 国产单亲对白刺激| 精品熟女少妇av免费看| 久久久久九九精品影院| 日韩欧美在线乱码| 国产成人91sexporn| 97超级碰碰碰精品色视频在线观看| 三级男女做爰猛烈吃奶摸视频| a级毛片a级免费在线| 精品久久久久久久久久免费视频| 免费看a级黄色片| 在线免费十八禁| 一个人看的www免费观看视频| 97在线视频观看| 可以在线观看毛片的网站| 日本 av在线| 麻豆精品久久久久久蜜桃| 99国产极品粉嫩在线观看| 亚洲欧美清纯卡通| 天天一区二区日本电影三级| 欧美日韩乱码在线| 国产精品电影一区二区三区| 一a级毛片在线观看| 成人亚洲精品av一区二区| 嫩草影视91久久| 99riav亚洲国产免费| 波野结衣二区三区在线| 在线观看美女被高潮喷水网站| 啦啦啦韩国在线观看视频| 国产精品,欧美在线| 欧美另类亚洲清纯唯美| 波多野结衣高清无吗| 国产欧美日韩一区二区精品| 九九在线视频观看精品| 国产高清有码在线观看视频| 少妇的逼水好多| 国产精品三级大全| 国产高清三级在线| 热99在线观看视频| av黄色大香蕉| 看非洲黑人一级黄片| 久久99热6这里只有精品| 亚洲精品在线观看二区| 精品一区二区三区视频在线| 女同久久另类99精品国产91| 国产成人a∨麻豆精品| 免费看光身美女| 乱人视频在线观看| 99久久精品一区二区三区| 精品久久国产蜜桃| 国产高清三级在线| 精品一区二区免费观看| 久久人妻av系列| 国产伦精品一区二区三区视频9| 成年女人毛片免费观看观看9| 色在线成人网| 成人鲁丝片一二三区免费| 中文字幕av成人在线电影| 1024手机看黄色片| 悠悠久久av| 亚洲av免费高清在线观看| 在线观看免费视频日本深夜| 美女cb高潮喷水在线观看| .国产精品久久| 国产精品久久久久久av不卡| 成年版毛片免费区| 亚洲国产精品成人综合色| 两性午夜刺激爽爽歪歪视频在线观看| 激情 狠狠 欧美| 亚洲性久久影院| 国产一区二区亚洲精品在线观看| 国产探花在线观看一区二区| 精品久久国产蜜桃| 国产大屁股一区二区在线视频| 波多野结衣高清无吗| 18禁在线无遮挡免费观看视频 | 婷婷精品国产亚洲av| or卡值多少钱| 国内久久婷婷六月综合欲色啪| 国产高清视频在线播放一区| 69人妻影院| 中文字幕av在线有码专区| 亚洲欧美成人综合另类久久久 | 别揉我奶头~嗯~啊~动态视频| 黄色日韩在线| 大香蕉久久网|