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

    賽車懸架轉(zhuǎn)向桿系的優(yōu)化設(shè)計(jì)

    2016-09-27 02:35:31楊永旺張代勝
    關(guān)鍵詞:約束方程桿系橫臂

    楊永旺, 張代勝

    (合肥工業(yè)大學(xué) 機(jī)械與汽車工程學(xué)院,安徽 合肥 230009)

    ?

    賽車懸架轉(zhuǎn)向桿系的優(yōu)化設(shè)計(jì)

    楊永旺,張代勝

    (合肥工業(yè)大學(xué) 機(jī)械與汽車工程學(xué)院,安徽 合肥230009)

    文章運(yùn)用自然坐標(biāo)法建立了賽車懸架轉(zhuǎn)向桿系的運(yùn)動(dòng)學(xué)模型,推導(dǎo)了約束方程,給出了賽車轉(zhuǎn)向時(shí)車輪偏轉(zhuǎn)角度和輪胎上下跳動(dòng)時(shí)前束值的計(jì)算方法;采用加權(quán)合并法,建立了綜合考慮賽車輪胎磨損、運(yùn)動(dòng)協(xié)調(diào)性和機(jī)動(dòng)性的非線性優(yōu)化目標(biāo)函數(shù),并提出了一種GA-PSO(Genetic Algorithm-Particle Swarm Optimization)混合算法進(jìn)行求解。針對(duì)特定賽車進(jìn)行了優(yōu)化設(shè)計(jì),優(yōu)化結(jié)果表明了該方法是有效的。

    自然坐標(biāo)法;GA-PSO混合算法;懸架轉(zhuǎn)向;優(yōu)化設(shè)計(jì);賽車

    0 引  言

    賽車懸架轉(zhuǎn)向桿系設(shè)計(jì)合理與否直接關(guān)系到賽車的操縱性能和輪胎的使用壽命。目前多數(shù)賽車采用雙橫臂獨(dú)立懸架和齒輪齒條機(jī)械式轉(zhuǎn)向器,雙橫臂獨(dú)立懸架與轉(zhuǎn)向桿系構(gòu)成復(fù)雜的空間機(jī)構(gòu),因此對(duì)賽車的運(yùn)動(dòng)分析變得較為困難。

    針對(duì)汽車懸架轉(zhuǎn)向系統(tǒng)的優(yōu)化設(shè)計(jì)問(wèn)題,國(guó)內(nèi)外學(xué)者已經(jīng)進(jìn)行了許多研究。文獻(xiàn)[1]對(duì)某整體式轉(zhuǎn)向機(jī)構(gòu)進(jìn)行了穩(wěn)健設(shè)計(jì),減少了不確定性運(yùn)動(dòng)誤差對(duì)汽車運(yùn)動(dòng)軌跡的影響;文獻(xiàn)[2-3]借助ADAMS軟件對(duì)懸架桿系與轉(zhuǎn)向桿系的運(yùn)行協(xié)調(diào)性進(jìn)行了優(yōu)化分析;文獻(xiàn)[4]運(yùn)用空間RSSR四桿機(jī)構(gòu)的旋轉(zhuǎn)矢量法建立了轉(zhuǎn)向桿系空間結(jié)構(gòu)模型并對(duì)車輪轉(zhuǎn)角進(jìn)行了分析;文獻(xiàn)[5-6]應(yīng)用空間機(jī)構(gòu)運(yùn)動(dòng)學(xué)原理對(duì)雙橫臂獨(dú)立懸架轉(zhuǎn)向梯形機(jī)構(gòu)進(jìn)行了優(yōu)化設(shè)計(jì);文獻(xiàn)[7]運(yùn)用多體動(dòng)力學(xué)中坐標(biāo)轉(zhuǎn)換原理實(shí)現(xiàn)了雙橫臂獨(dú)立懸架機(jī)構(gòu)的運(yùn)動(dòng)特性分析;文獻(xiàn)[8-9]分別采用分組實(shí)數(shù)遺傳算法和改進(jìn)粒子群優(yōu)化算法對(duì)整體式轉(zhuǎn)向梯形機(jī)構(gòu)的優(yōu)化模型進(jìn)行了求解;文獻(xiàn)[10]應(yīng)用目標(biāo)級(jí)聯(lián)方法對(duì)獨(dú)立懸架轉(zhuǎn)向機(jī)構(gòu)進(jìn)行了運(yùn)動(dòng)學(xué)優(yōu)化。

    本文對(duì)賽車懸架轉(zhuǎn)向桿系進(jìn)行優(yōu)化設(shè)計(jì),首先作出如下假設(shè):① 所有桿件均為剛體;② 忽略鉸接處的間隙。

    1 賽車懸架轉(zhuǎn)向桿系運(yùn)動(dòng)學(xué)分析

    1.1自然坐標(biāo)

    某賽車單側(cè)懸架轉(zhuǎn)向桿系運(yùn)動(dòng)學(xué)模型如圖1所示。其中,慣性基原點(diǎn)為O,X軸位于賽車縱向?qū)ΨQ面內(nèi)且指向車尾,Z軸垂直地面向上,X、Y、Z軸滿足右手法則;A點(diǎn)和B點(diǎn)為上橫臂與車架鉸接點(diǎn),C點(diǎn)和D點(diǎn)為下橫臂與車架鉸接點(diǎn),上、下橫臂分別繞AB軸和CD軸轉(zhuǎn)動(dòng);E、F分別為轉(zhuǎn)向節(jié)上、下球銷中心;MN為轉(zhuǎn)向橫拉桿,M點(diǎn)為轉(zhuǎn)向節(jié)臂連接點(diǎn),N點(diǎn)為轉(zhuǎn)向梯形斷開點(diǎn),且M、N兩處均為球鉸副;G點(diǎn)為車輪中心點(diǎn),L點(diǎn)為輪軸與主銷軸線交點(diǎn)。

    圖1 懸架轉(zhuǎn)向桿系運(yùn)動(dòng)學(xué)模型

    1.2約束方程

    上橫臂的自然坐標(biāo)為B和E2點(diǎn)的坐標(biāo),首先根據(jù)剛性約束給出3個(gè)約束方程為:

    (1)

    (2)

    (3)

    由于上橫臂只有1個(gè)自由度,因此根據(jù)旋轉(zhuǎn)副補(bǔ)充2個(gè)約束方程[11],即

    (4)

    其中,B0為B點(diǎn)固結(jié)在慣性基上的點(diǎn),其坐標(biāo)值為常量。兩向量叉乘得到3個(gè)代數(shù)方程,但其中只有2個(gè)是獨(dú)立的,因此(4)式實(shí)際包含2個(gè)約束方程。

    同理,下橫臂共有5個(gè)約束方程,即

    (5)

    (6)

    (7)

    (8)

    其中,D0為D點(diǎn)固結(jié)在慣性基上的點(diǎn),其坐標(biāo)值為常量。

    轉(zhuǎn)向橫拉桿有1個(gè)剛性約束方程,即

    (9)

    轉(zhuǎn)向節(jié)包含E、F、G、L、M5個(gè)基本點(diǎn),需要9個(gè)剛性約束方程,即

    (10)

    (11)

    (12)

    (13)

    (14)

    (15)

    (16)

    (17)

    (18)

    因此,該運(yùn)動(dòng)學(xué)模型共有20個(gè)約束方程,剩余4個(gè)自由度,這4個(gè)自由度分別為斷開點(diǎn)N的三維坐標(biāo)和輪心G點(diǎn)的Z軸坐標(biāo)。

    1.3運(yùn)動(dòng)學(xué)參數(shù)

    (19)

    其中,L0、G0分別為L(zhǎng)、G的初始坐標(biāo)。

    由于賽車轉(zhuǎn)向器與車架剛性連接,在轉(zhuǎn)向過(guò)程中點(diǎn)N只會(huì)沿Y軸移動(dòng),即只有Y方向的坐標(biāo)值發(fā)生變化,因此易得到方向盤轉(zhuǎn)角θSW與N點(diǎn)坐標(biāo)的函數(shù)關(guān)系,記為:

    (20)

    聯(lián)立(19)式和(20)式即可得到方向盤轉(zhuǎn)角與車輪轉(zhuǎn)角的關(guān)系為:

    (21)

    (2) 當(dāng)N點(diǎn)固定不動(dòng),車輪發(fā)生跳動(dòng)即GZ發(fā)生變化時(shí),同樣可求得1組基本點(diǎn)的對(duì)應(yīng)坐標(biāo)值。此時(shí),對(duì)應(yīng)的車輪前束值為:

    (22)

    其中,D為車輪直徑。

    2 賽車懸架轉(zhuǎn)向桿系非線性優(yōu)化模型

    2.1目標(biāo)函數(shù)

    為了綜合考慮賽車輪胎磨損、懸架轉(zhuǎn)向運(yùn)動(dòng)的協(xié)調(diào)性以及機(jī)動(dòng)性,目標(biāo)函數(shù)涉及以下3個(gè)方面。

    (1) 賽車在轉(zhuǎn)向過(guò)程中,為減少輪胎磨損,應(yīng)保證正確的內(nèi)外輪轉(zhuǎn)角關(guān)系。但由于輪胎側(cè)偏角的存在,特別是賽車高速行駛時(shí),100%的Ackerman不一定滿足實(shí)際要求。為修正輪胎側(cè)偏角對(duì)轉(zhuǎn)角誤差的影響,平均百分比Ackerman校正率η被應(yīng)用于優(yōu)化目標(biāo)函數(shù)中,即

    (23)

    其中,θi和θo分別為內(nèi)側(cè)車輪轉(zhuǎn)角和外側(cè)車輪轉(zhuǎn)角;θo′為根據(jù)Ackerman原理計(jì)算得到的外側(cè)車輪轉(zhuǎn)角理論值[5];w0(θo)為權(quán)函數(shù),表征不同轉(zhuǎn)角的使用情況,按照文獻(xiàn)[5]選取。

    (2) 賽車懸架轉(zhuǎn)向桿系的匹配設(shè)計(jì)應(yīng)處理好懸架桿系與轉(zhuǎn)向桿系的干涉問(wèn)題,即應(yīng)實(shí)現(xiàn)車輪跳動(dòng)過(guò)程中前束變化量最小,因此目標(biāo)函數(shù)如下:

    (24)

    其中,Δ0為初始時(shí)車輪的前束值。

    (3) 受比賽場(chǎng)地的限制,賽車需滿足最小轉(zhuǎn)彎半徑的要求,即應(yīng)滿足方向盤最大轉(zhuǎn)角一定時(shí)外側(cè)車輪最大轉(zhuǎn)角盡可能接近目標(biāo)值θC,因此目標(biāo)函數(shù)如下:

    (25)

    綜合以上3個(gè)目標(biāo)函數(shù),采用加權(quán)合并法建立的總目標(biāo)函數(shù)如下:

    F(X)=min[w1F1(X)+

    (26)

    其中,w1、w2和w3均為加權(quán)系數(shù),為了不失一般性,需滿足w1+w2+w3=1。

    2.2設(shè)計(jì)變量

    考慮到賽車懸架轉(zhuǎn)向桿系關(guān)于XOZ平面對(duì)稱,選取一側(cè)轉(zhuǎn)向橫拉桿兩端點(diǎn)坐標(biāo)作為優(yōu)化變量,即

    (27)

    2.3約束條件

    根據(jù)轉(zhuǎn)向器的安裝位置以及避免與懸架、車架、輪輞發(fā)生干涉等幾個(gè)方面來(lái)確定設(shè)計(jì)變量的取值范圍,記為:

    (28)

    3 求解非線性優(yōu)化模型

    3.1GA-PSO混合算法

    GA-PSO(Genetic Algorithm-Particle Swarm Optimization)混合算法是一種融合了遺傳算法全局性與粒子群算法記憶性的智能優(yōu)化算法。

    傳統(tǒng)遺傳算法產(chǎn)生初始種群后,首先計(jì)算適應(yīng)度,再進(jìn)行選擇、交叉、變異一系列操作,形成新種群[12]。該過(guò)程中父代的優(yōu)秀個(gè)體沒有得到保護(hù),即父代中的優(yōu)秀個(gè)體可能被交叉產(chǎn)生的劣質(zhì)后代所取代。交叉產(chǎn)生的優(yōu)良個(gè)體也沒有得到保留,可能在變異環(huán)節(jié)被破壞,因此不容易快速收斂。粒子群算法能夠?qū)崟r(shí)記憶個(gè)體和全局最優(yōu)信息,通過(guò)追隨當(dāng)前的最優(yōu)粒子在解空間進(jìn)行搜索,具有加快收斂的能力。

    針對(duì)遺傳算法和粒子群算法各自的優(yōu)缺點(diǎn),本文提出一種GA-PSO混合算法,其流程如圖2所示。

    首先對(duì)第n代群體P進(jìn)行適應(yīng)度計(jì)算,按優(yōu)劣排序得到群體P′,對(duì)P′中的前w4%較優(yōu)個(gè)體進(jìn)行保留,并按照基本粒子群算法的速度和位置更新公式[13]進(jìn)行加速,加速后的粒子群Q1進(jìn)入第n+1代群體Q;其次,按照聯(lián)賽選擇法,從P′中選出w5%的個(gè)體進(jìn)行實(shí)數(shù)中間重組得到Q2,以確保子代種群的多樣性,避免陷于局部最優(yōu);最后對(duì)P′中排序靠后的w6%個(gè)體進(jìn)行變異操作得到Q3,以保證局部隨機(jī)搜索的能力。需要注意的是,為了保證較好的全局尋優(yōu)能力,w5取值通常比w4和w6大。

    圖2 一種GA-PSO混合算法流程

    3.2優(yōu)化步驟

    針對(duì)本文建立的非線性優(yōu)化模型,采用GA-PSO混合算法對(duì)其進(jìn)行求解,具體步驟如下:

    (1) 采用實(shí)數(shù)編碼,按(28)式所給的變量區(qū)間,通過(guò)隨機(jī)方法產(chǎn)生初始種群P,設(shè)置種群數(shù)為50,最大迭代次數(shù)T=50。

    (2) 由于目標(biāo)函數(shù)(26)式始終為正值,且為最小化問(wèn)題,所以確定適應(yīng)度函數(shù)為1/F(X)。計(jì)算種群適應(yīng)度,按適應(yīng)度值大小排序得到種群P′。

    (3) 按PSO算法[13]更新P′中前20%的個(gè)體,得到種群Q1。其中,取慣性因子為0.8,取加速因子c1=1.494 45,c2=1.515 54。

    (4) 采用聯(lián)賽選擇法從P′中選出70%的個(gè)體,以100%的交叉概率,根據(jù)如下方程進(jìn)行交叉得到Q2,即

    (29)

    其中,λ=rand(0,1)。

    (5) 對(duì)P′中剩余10%的最差個(gè)體全部執(zhí)行變異得到Q3。

    (6) 至此得到完整的子代群體Q,檢查Q中是否有個(gè)體已跳出解空間,即是否仍然滿足約束條件,若滿足則接受,否則用變量對(duì)應(yīng)元素邊界值替換該個(gè)體的變量值。

    (7) 檢驗(yàn)是否達(dá)到迭代次數(shù),滿足則停止迭代,否則轉(zhuǎn)到步驟(2)。

    4 實(shí)例分析

    某賽車懸架轉(zhuǎn)向桿系各關(guān)鍵點(diǎn)初始坐標(biāo)值見表1所列。

    設(shè)計(jì)變量的上、下限為:

    Xmin=(-30,220,-80,-78,550,-70),

    Xmax=(-10,250,-68,-50,575,0)。

    賽車方向盤最大轉(zhuǎn)角θSWmax為90°,車輪上、下跳動(dòng)區(qū)間為(-25,25),取平均百分比Ackerman校正率η=80%,加權(quán)系數(shù)w1=w3=0.4,w2=0.2。

    采用GA-PSO混合算法,最終得到優(yōu)化結(jié)果如下:

    X=(-10.00,222.88,-68.00,-52.75,

    575.00,-47.75)。

    適應(yīng)度隨迭代次數(shù)的變化曲線如圖3所示。由圖3可知,采用的GA-PSO混合算法能夠快速有效收斂。

    表1 關(guān)鍵點(diǎn)初始坐標(biāo)值 mm

    圖3 迭代適應(yīng)度曲線

    優(yōu)化前、后賽車各目標(biāo)函數(shù)值變化如圖4所示。

    由圖4a可看出,優(yōu)化后賽車左、右車輪實(shí)際轉(zhuǎn)角更接近目標(biāo)值曲線,能夠有效減少輪胎磨損;優(yōu)化前的轉(zhuǎn)角曲線反映出賽車具有無(wú)法實(shí)現(xiàn)外側(cè)車輪較大轉(zhuǎn)角的問(wèn)題;優(yōu)化后的曲線表明,在同樣最大方向盤轉(zhuǎn)角的情況下,優(yōu)化后的賽車可以獲得理想的外側(cè)車輪轉(zhuǎn)角,增強(qiáng)了賽車的過(guò)彎能力,提高了賽車機(jī)動(dòng)性。由圖4b可以看出,賽車經(jīng)過(guò)優(yōu)化后,由車輪上下跳動(dòng)所引起的前束變化值大幅度減小,懸架與轉(zhuǎn)向桿系間的運(yùn)動(dòng)協(xié)調(diào)性明顯得到提高。

    圖4 優(yōu)化前、后賽車各目標(biāo)函數(shù)值變化曲線

    5 結(jié)束語(yǔ)

    綜合考慮賽車輪胎磨損、懸架轉(zhuǎn)向運(yùn)動(dòng)協(xié)調(diào)性以及機(jī)動(dòng)性能夠更全面地反映賽車的設(shè)計(jì)要求。優(yōu)化結(jié)果表明本文所提出的目標(biāo)函數(shù)具有一定的工程實(shí)用價(jià)值。

    自然坐標(biāo)法可以實(shí)現(xiàn)快速建模,由于其約束方程均為線性或二次方程,因此求解方便,易于編程,本文的實(shí)例分析說(shuō)明了該方法用于復(fù)雜空間機(jī)構(gòu)問(wèn)題求解的有效性和方便性,可應(yīng)用于麥弗遜懸架、雙前橋牽引車轉(zhuǎn)向系等系統(tǒng)的建模。

    針對(duì)非線性多目標(biāo)優(yōu)化問(wèn)題,本文提出一種GA-PSO混合算法。通過(guò)采取保護(hù)父代中和交叉產(chǎn)生的優(yōu)良個(gè)體,并對(duì)優(yōu)良個(gè)體進(jìn)行PSO加速等措施,使進(jìn)化具有較好的方向性且收斂速度加快。同時(shí)采取100%的交叉概率有效保證了子代種群的多樣性。優(yōu)化結(jié)果表明,GA-PSO混合算法能夠較好地求解非線性優(yōu)化問(wèn)題,具有一定的通用性。

    [1]張蕾,董恩國(guó).基于多目標(biāo)優(yōu)化的車輛轉(zhuǎn)向機(jī)構(gòu)穩(wěn)健設(shè)計(jì)[J].機(jī)械設(shè)計(jì),2015,32(3):51-55.

    [2]PALIWAL G,SUKUMAR N,GUPTA U,et al.Investigation and optimization of front suspension and steering geometrical compatibility:SAE Technical Paper 2015-01-0492[R].SAE,2015.DOI:10.4271/2015-01-0492.

    [3]向鐵明,沈理真.基于ADAMS的某賽車前懸架桿系優(yōu)化[J].汽車工程,2014,36(4):486-490.

    [4]古玉鋒,呂彭民,單增海,等.轉(zhuǎn)向桿系空間結(jié)構(gòu)非線性建模與分析[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2014,45(10):7-14.

    [5]楊永旺.基于空間機(jī)構(gòu)學(xué)的FSAE賽車轉(zhuǎn)向梯形機(jī)構(gòu)優(yōu)化設(shè)計(jì)[J].農(nóng)業(yè)裝備與車輛工程,2014,52(8):27-30.

    [6]方錫邦,王其東,關(guān)柯.雙橫臂獨(dú)立懸架轉(zhuǎn)向傳動(dòng)機(jī)構(gòu)的優(yōu)化設(shè)計(jì)[J].合肥工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2001,24(5):950-954.

    [7]王其東,趙韓,李巖,等.汽車雙橫臂式獨(dú)立懸架機(jī)構(gòu)運(yùn)動(dòng)特性分析[J].合肥工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2001,24(6):1066-1071.

    [8]劉嶺,閆光榮,雷毅,等.基于改進(jìn)粒子群算法的車輛轉(zhuǎn)向梯形機(jī)構(gòu)優(yōu)化[J].農(nóng)業(yè)工程學(xué)報(bào),2013,29(10):76-82.

    [9]張勇,王福林,玄登影,等.分組實(shí)數(shù)遺傳算法優(yōu)化拖拉機(jī)轉(zhuǎn)向梯形問(wèn)題研究[J].農(nóng)機(jī)化研究,2014(12):92-95.

    [10]BLANCO J,MUNOZ L.Target cascading optimization of the kinematics of a steering system coupled to a double wishbone suspension system of a hybrid off-road vehicle:SAE Technical Paper 2013-01-0967[R].SAE,2013.DOI:10.4271/2013-01-0967.

    [11]DE JAL N J G,BAYO E.Kinematic and dynamic simulation of multibodt stytems[M].New York:Springer-Verlag,1994:47-56.

    [12]王小平,曹立明.遺傳算法:理論、應(yīng)用及軟件實(shí)現(xiàn)[M].西安:西安交通大學(xué)出版社,2002:28-38.

    [13]劉波.粒子群優(yōu)化算法及其工程應(yīng)用[M].北京:機(jī)械工業(yè)出版社,2010:34-35.

    (責(zé)任編輯胡亞敏)

    Optimization design of suspension and steering linkages in racing car

    YANG Yongwang, ZHANG Daisheng

    (School of Machinery and Automobile Engineering, Hefei University of Technology, Hefei 230009, China)

    By using the natural coordinates method, a kinematical model of suspension and steering linkages in racing car was established. The constraint equations of this system were deduced. The calculation method of both the wheel steering angle during cornering and the toe caused by wheel bouncing was given. By using the weighted combination method, a nonlinear optimization object function was established, which took the tire wear, motion compatibility and mobility into account. Then a hybrid algorithm of genetic algorithm-particle swarm optimization (GA-PSO) was put forward to solve the problem. Finally, the effectiveness of the proposed method was verified by an example of a racing car.

    natural coordinates method; hybrid algorithm of genetic algorithm-particle swarm optimization (GA-PSO); suspension and steering; optimization design; racing car

    2015-03-20;

    2015-05-06

    合肥工業(yè)大學(xué)產(chǎn)學(xué)研校企合作資助項(xiàng)目(W2014JSKF0458)

    楊永旺(1989-),男,山西大同人,合肥工業(yè)大學(xué)碩士生;

    張代勝(1962-),男,安徽青陽(yáng)人,合肥工業(yè)大學(xué)教授,碩士生導(dǎo)師.

    10.3969/j.issn.1003-5060.2016.07.005

    U463.45

    A

    1003-5060(2016)07-0887-05

    猜你喜歡
    約束方程桿系橫臂
    移動(dòng)機(jī)器人動(dòng)力學(xué)方程的約束違約穩(wěn)定方法
    雙橫臂懸架導(dǎo)向機(jī)構(gòu)運(yùn)動(dòng)學(xué)分析及優(yōu)化
    用轉(zhuǎn)置矩陣求解靜不定桁架的內(nèi)力
    含剛性斜桿的平面有側(cè)移剛架內(nèi)力計(jì)算1)
    基于GPU的桿系離散元并行算法在大型工程結(jié)構(gòu)中的應(yīng)用
    礦井巷道三維建模方法探討
    GEVO16柴油機(jī)氣門鎖夾逃逸原因分析
    基于ADAMS的車輛轉(zhuǎn)向桿系動(dòng)態(tài)仿真分析
    多體系統(tǒng)指標(biāo)2運(yùn)動(dòng)方程HHT方法違約校正1)
    空間桿系結(jié)構(gòu)動(dòng)力失穩(wěn)區(qū)域研究
    99久国产av精品国产电影| 久久亚洲国产成人精品v| 韩国av在线不卡| 欧美日韩av久久| 精品一区二区三区视频在线| 丝袜美足系列| 欧美另类一区| 99久久中文字幕三级久久日本| 最新中文字幕久久久久| 免费人妻精品一区二区三区视频| 大话2 男鬼变身卡| 亚洲精品aⅴ在线观看| 90打野战视频偷拍视频| 乱人伦中国视频| 国产亚洲精品第一综合不卡 | 黄片无遮挡物在线观看| 亚洲av中文av极速乱| 欧美bdsm另类| 一区二区av电影网| 国产精品国产三级国产av玫瑰| 极品人妻少妇av视频| 少妇人妻 视频| 亚洲国产av新网站| 国产一区二区在线观看av| 亚洲国产av影院在线观看| 亚洲成国产人片在线观看| 2018国产大陆天天弄谢| 午夜久久久在线观看| 啦啦啦在线观看免费高清www| 国产一区二区在线观看日韩| 亚洲欧美成人综合另类久久久| 久久久久精品久久久久真实原创| 久久精品国产自在天天线| 99热国产这里只有精品6| 国产国语露脸激情在线看| 亚洲成国产人片在线观看| 黑人高潮一二区| 各种免费的搞黄视频| av片东京热男人的天堂| 精品国产一区二区久久| 国产欧美日韩综合在线一区二区| 国产欧美另类精品又又久久亚洲欧美| 欧美少妇被猛烈插入视频| 视频中文字幕在线观看| 国产av国产精品国产| 国语对白做爰xxxⅹ性视频网站| 亚洲av免费高清在线观看| 亚洲欧美日韩另类电影网站| 在线观看www视频免费| 久久精品国产亚洲av天美| 在线观看一区二区三区激情| 99re6热这里在线精品视频| 欧美日韩综合久久久久久| 中文天堂在线官网| 久久免费观看电影| 亚洲四区av| 18禁观看日本| 99热6这里只有精品| 一级毛片我不卡| 午夜91福利影院| 亚洲 欧美一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 国产极品天堂在线| 日韩一本色道免费dvd| 精品人妻偷拍中文字幕| 日本黄大片高清| 色网站视频免费| 国产一区二区激情短视频 | 国产av码专区亚洲av| 99久久精品国产国产毛片| 激情五月婷婷亚洲| 国产麻豆69| 日韩av免费高清视频| 国产成人精品婷婷| 大片电影免费在线观看免费| 亚洲综合色网址| 国产男人的电影天堂91| 久久热在线av| 丝袜美足系列| 人妻少妇偷人精品九色| 天堂俺去俺来也www色官网| 国产爽快片一区二区三区| 日韩视频在线欧美| 欧美+日韩+精品| 亚洲成色77777| 欧美激情极品国产一区二区三区 | 国产精品三级大全| 黑人高潮一二区| 看免费av毛片| 最近手机中文字幕大全| 老司机亚洲免费影院| 亚洲精品国产av成人精品| 老司机亚洲免费影院| 韩国高清视频一区二区三区| 亚洲国产精品国产精品| 80岁老熟妇乱子伦牲交| 久久av网站| 国产高清不卡午夜福利| 午夜福利网站1000一区二区三区| av播播在线观看一区| 两性夫妻黄色片 | 激情五月婷婷亚洲| 91在线精品国自产拍蜜月| 99久久精品国产国产毛片| 国产精品99久久99久久久不卡 | 亚洲中文av在线| 美女脱内裤让男人舔精品视频| 免费观看在线日韩| 永久网站在线| 永久免费av网站大全| 免费人妻精品一区二区三区视频| 免费观看a级毛片全部| 久久精品国产自在天天线| 国产欧美另类精品又又久久亚洲欧美| √禁漫天堂资源中文www| 免费大片18禁| 毛片一级片免费看久久久久| 成人无遮挡网站| av不卡在线播放| 女性被躁到高潮视频| 少妇人妻精品综合一区二区| 女人被躁到高潮嗷嗷叫费观| av黄色大香蕉| 久久久欧美国产精品| 波野结衣二区三区在线| 日韩av不卡免费在线播放| 看非洲黑人一级黄片| 中文字幕人妻熟女乱码| 日韩伦理黄色片| 少妇人妻久久综合中文| 一本—道久久a久久精品蜜桃钙片| 日韩中文字幕视频在线看片| 十分钟在线观看高清视频www| 一区在线观看完整版| 久久这里只有精品19| 欧美精品国产亚洲| 亚洲欧洲精品一区二区精品久久久 | 香蕉国产在线看| 国产精品一区二区在线观看99| 久久精品熟女亚洲av麻豆精品| 狠狠精品人妻久久久久久综合| 亚洲欧美日韩另类电影网站| 一本久久精品| 美女福利国产在线| 国产成人免费无遮挡视频| 91在线精品国自产拍蜜月| 老司机亚洲免费影院| 亚洲美女搞黄在线观看| 久久婷婷青草| 天堂中文最新版在线下载| 精品国产乱码久久久久久小说| 国产高清国产精品国产三级| 婷婷色综合www| 18禁观看日本| 欧美人与善性xxx| 街头女战士在线观看网站| 亚洲av福利一区| 日日摸夜夜添夜夜爱| 狂野欧美激情性bbbbbb| 久久精品国产亚洲av天美| 国产成人一区二区在线| 一级片免费观看大全| 26uuu在线亚洲综合色| 久久综合国产亚洲精品| 热99国产精品久久久久久7| 日本91视频免费播放| 精品酒店卫生间| 欧美日韩精品成人综合77777| 国产精品免费大片| 久久 成人 亚洲| 欧美国产精品va在线观看不卡| 国产精品蜜桃在线观看| 午夜免费男女啪啪视频观看| 中文字幕另类日韩欧美亚洲嫩草| 99香蕉大伊视频| 制服人妻中文乱码| 亚洲综合精品二区| 欧美成人午夜免费资源| 蜜桃在线观看..| 寂寞人妻少妇视频99o| 十分钟在线观看高清视频www| av在线播放精品| 亚洲精品第二区| 午夜福利视频在线观看免费| 亚洲精品美女久久av网站| 久久免费观看电影| 国产欧美另类精品又又久久亚洲欧美| 9191精品国产免费久久| 久热久热在线精品观看| 在线观看免费视频网站a站| 全区人妻精品视频| 欧美人与性动交α欧美精品济南到 | 飞空精品影院首页| 亚洲欧美色中文字幕在线| 全区人妻精品视频| 国产欧美日韩综合在线一区二区| 亚洲欧洲日产国产| 春色校园在线视频观看| 欧美日韩亚洲高清精品| 两个人免费观看高清视频| 久久久久国产网址| 国产成人精品无人区| 亚洲欧洲国产日韩| 丰满乱子伦码专区| 一区二区三区乱码不卡18| 乱人伦中国视频| 午夜福利,免费看| 日本黄大片高清| 久久精品久久久久久噜噜老黄| 丰满乱子伦码专区| 国产精品嫩草影院av在线观看| 久久人人97超碰香蕉20202| 人人妻人人爽人人添夜夜欢视频| 国产精品 国内视频| 自拍欧美九色日韩亚洲蝌蚪91| 高清不卡的av网站| 一边摸一边做爽爽视频免费| 男女无遮挡免费网站观看| 天天操日日干夜夜撸| 黑人猛操日本美女一级片| 夜夜骑夜夜射夜夜干| av黄色大香蕉| 精品人妻熟女毛片av久久网站| 国产精品 国内视频| 精品久久蜜臀av无| 美女内射精品一级片tv| 国产精品99久久99久久久不卡 | xxx大片免费视频| 一二三四中文在线观看免费高清| 亚洲精品中文字幕在线视频| 国产毛片在线视频| 一级,二级,三级黄色视频| 极品人妻少妇av视频| 丰满乱子伦码专区| 日韩 亚洲 欧美在线| 秋霞伦理黄片| 日本-黄色视频高清免费观看| 大香蕉久久成人网| 国产国拍精品亚洲av在线观看| 美女福利国产在线| 亚洲成人手机| 一边摸一边做爽爽视频免费| 亚洲精品久久成人aⅴ小说| 人成视频在线观看免费观看| 欧美精品国产亚洲| 五月天丁香电影| 少妇的逼水好多| 制服丝袜香蕉在线| 大码成人一级视频| 精品久久国产蜜桃| 黄色怎么调成土黄色| 久久精品久久精品一区二区三区| 午夜激情av网站| 久久精品熟女亚洲av麻豆精品| 亚洲欧美精品自产自拍| 国产一区二区在线观看av| 精品卡一卡二卡四卡免费| 精品99又大又爽又粗少妇毛片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 免费看av在线观看网站| 美女xxoo啪啪120秒动态图| 国产精品久久久久久精品电影小说| 国产淫语在线视频| 91午夜精品亚洲一区二区三区| 国产亚洲午夜精品一区二区久久| 大香蕉97超碰在线| 免费少妇av软件| 精品福利永久在线观看| 国国产精品蜜臀av免费| 国产日韩欧美视频二区| 国产欧美亚洲国产| www.色视频.com| 人妻人人澡人人爽人人| 亚洲美女搞黄在线观看| 欧美激情国产日韩精品一区| 国产极品天堂在线| 色视频在线一区二区三区| 999精品在线视频| 亚洲精品美女久久av网站| 亚洲激情五月婷婷啪啪| xxx大片免费视频| 国产精品久久久久久精品电影小说| 久久精品国产亚洲av天美| 免费av不卡在线播放| 五月伊人婷婷丁香| 国产亚洲欧美精品永久| 亚洲四区av| 日本欧美视频一区| 最近最新中文字幕大全免费视频 | 欧美少妇被猛烈插入视频| av电影中文网址| 免费黄频网站在线观看国产| 亚洲欧美一区二区三区黑人 | 国产男女内射视频| 最近2019中文字幕mv第一页| 亚洲精华国产精华液的使用体验| 婷婷色麻豆天堂久久| 天堂8中文在线网| 亚洲国产欧美日韩在线播放| 九草在线视频观看| 国产一区二区在线观看av| 大片免费播放器 马上看| 国产深夜福利视频在线观看| 日本-黄色视频高清免费观看| 色哟哟·www| 国产黄色视频一区二区在线观看| 巨乳人妻的诱惑在线观看| 永久网站在线| 日韩三级伦理在线观看| 中文字幕精品免费在线观看视频 | 日本vs欧美在线观看视频| 咕卡用的链子| 午夜影院在线不卡| 2022亚洲国产成人精品| 久久久久久久精品精品| 亚洲精品一区蜜桃| 亚洲精品自拍成人| 青青草视频在线视频观看| 人妻系列 视频| 久久鲁丝午夜福利片| 各种免费的搞黄视频| 亚洲五月色婷婷综合| 欧美老熟妇乱子伦牲交| 亚洲精品乱码久久久久久按摩| 啦啦啦在线观看免费高清www| 中文乱码字字幕精品一区二区三区| 赤兔流量卡办理| 国产精品一国产av| 久久精品久久精品一区二区三区| 欧美成人午夜免费资源| 成年美女黄网站色视频大全免费| 午夜激情久久久久久久| a 毛片基地| 亚洲精品久久成人aⅴ小说| 欧美日韩视频精品一区| 少妇猛男粗大的猛烈进出视频| 涩涩av久久男人的天堂| 寂寞人妻少妇视频99o| 综合色丁香网| 欧美成人午夜免费资源| 高清黄色对白视频在线免费看| 91久久精品国产一区二区三区| 亚洲精品中文字幕在线视频| 精品人妻熟女毛片av久久网站| 一级毛片 在线播放| 国产精品一区www在线观看| 99国产精品免费福利视频| 日韩精品有码人妻一区| 欧美 亚洲 国产 日韩一| 亚洲精品久久成人aⅴ小说| 国产精品无大码| 亚洲情色 制服丝袜| 飞空精品影院首页| 中文字幕另类日韩欧美亚洲嫩草| 自拍欧美九色日韩亚洲蝌蚪91| 青春草国产在线视频| av有码第一页| 少妇的丰满在线观看| 一二三四在线观看免费中文在 | 汤姆久久久久久久影院中文字幕| 天堂8中文在线网| 久久精品国产a三级三级三级| 男女下面插进去视频免费观看 | 高清视频免费观看一区二区| 少妇人妻 视频| 美女国产高潮福利片在线看| 国产成人午夜福利电影在线观看| 99热国产这里只有精品6| 春色校园在线视频观看| 国产在线一区二区三区精| 亚洲av欧美aⅴ国产| 综合色丁香网| 亚洲伊人色综图| 久久久久久久久久久久大奶| 亚洲色图综合在线观看| 久久精品国产鲁丝片午夜精品| 成年美女黄网站色视频大全免费| 国产精品三级大全| 在线亚洲精品国产二区图片欧美| 精品一品国产午夜福利视频| 成人午夜精彩视频在线观看| 国产xxxxx性猛交| 亚洲av综合色区一区| 中文字幕精品免费在线观看视频 | 女性被躁到高潮视频| 校园人妻丝袜中文字幕| 一二三四中文在线观看免费高清| 国产成人aa在线观看| 亚洲一区二区三区欧美精品| 日日爽夜夜爽网站| 人妻 亚洲 视频| 一区在线观看完整版| 婷婷色综合www| 久久久久久久久久成人| 亚洲欧美成人综合另类久久久| 最近的中文字幕免费完整| 看十八女毛片水多多多| 中国三级夫妇交换| 亚洲国产精品成人久久小说| 久久久久精品人妻al黑| 亚洲天堂av无毛| 丝袜美足系列| 亚洲丝袜综合中文字幕| 免费观看a级毛片全部| 春色校园在线视频观看| 丁香六月天网| 亚洲精品,欧美精品| 亚洲成人手机| 久久这里有精品视频免费| 久久精品久久久久久噜噜老黄| 2021少妇久久久久久久久久久| 免费在线观看完整版高清| 七月丁香在线播放| 久久人人爽人人片av| 性高湖久久久久久久久免费观看| 一本—道久久a久久精品蜜桃钙片| 久久人人爽人人片av| 夫妻性生交免费视频一级片| 国产精品欧美亚洲77777| 免费黄网站久久成人精品| av网站免费在线观看视频| 国产成人av激情在线播放| 国产xxxxx性猛交| 日本色播在线视频| 精品少妇内射三级| 久久这里有精品视频免费| 亚洲欧洲国产日韩| 少妇人妻精品综合一区二区| 99视频精品全部免费 在线| 成人18禁高潮啪啪吃奶动态图| 精品亚洲成国产av| 街头女战士在线观看网站| 欧美日韩国产mv在线观看视频| 热99久久久久精品小说推荐| 精品一区二区免费观看| 99热6这里只有精品| 久久精品人人爽人人爽视色| av一本久久久久| 国产xxxxx性猛交| 色婷婷av一区二区三区视频| 啦啦啦啦在线视频资源| 国产女主播在线喷水免费视频网站| 亚洲色图综合在线观看| 久久精品aⅴ一区二区三区四区 | 人妻一区二区av| 国产精品人妻久久久影院| 亚洲精品中文字幕在线视频| 国产永久视频网站| 欧美性感艳星| 免费少妇av软件| 久久精品国产综合久久久 | 国产免费视频播放在线视频| av在线app专区| 亚洲精品aⅴ在线观看| 美女内射精品一级片tv| 亚洲精品美女久久久久99蜜臀 | 国产成人a∨麻豆精品| 日韩在线高清观看一区二区三区| 日本vs欧美在线观看视频| 亚洲av综合色区一区| 久久免费观看电影| 春色校园在线视频观看| 大香蕉久久网| 狂野欧美激情性xxxx在线观看| 久久久a久久爽久久v久久| 免费观看a级毛片全部| a级毛色黄片| 天天操日日干夜夜撸| videossex国产| 亚洲国产日韩一区二区| 精品久久国产蜜桃| 亚洲伊人色综图| 看非洲黑人一级黄片| 成人二区视频| 亚洲精品自拍成人| 五月天丁香电影| 国产熟女欧美一区二区| 国产精品久久久久久av不卡| 宅男免费午夜| 色哟哟·www| 男女无遮挡免费网站观看| 色吧在线观看| 天天躁夜夜躁狠狠躁躁| 91aial.com中文字幕在线观看| 51国产日韩欧美| 欧美日韩一区二区视频在线观看视频在线| 亚洲色图综合在线观看| a级片在线免费高清观看视频| 欧美少妇被猛烈插入视频| 日本欧美视频一区| 街头女战士在线观看网站| 免费观看无遮挡的男女| 色5月婷婷丁香| 伦理电影免费视频| 午夜福利乱码中文字幕| xxx大片免费视频| 国产亚洲av片在线观看秒播厂| 免费在线观看黄色视频的| 黄网站色视频无遮挡免费观看| 满18在线观看网站| 欧美 亚洲 国产 日韩一| 免费av不卡在线播放| 男男h啪啪无遮挡| 捣出白浆h1v1| 午夜免费鲁丝| 国产精品不卡视频一区二区| 夫妻午夜视频| 国产xxxxx性猛交| 18禁裸乳无遮挡动漫免费视频| 亚洲精品一区蜜桃| 精品国产一区二区久久| 亚洲精品456在线播放app| 日本av免费视频播放| 韩国av在线不卡| 交换朋友夫妻互换小说| 久久国内精品自在自线图片| 亚洲在久久综合| 人人妻人人添人人爽欧美一区卜| 国产又色又爽无遮挡免| 午夜激情久久久久久久| 国产精品麻豆人妻色哟哟久久| 国产成人免费观看mmmm| av在线观看视频网站免费| 亚洲五月色婷婷综合| 国产精品国产三级专区第一集| 麻豆乱淫一区二区| 18在线观看网站| 国产精品国产三级国产专区5o| 国产成人精品福利久久| 久久久久久久精品精品| 精品少妇黑人巨大在线播放| 丰满迷人的少妇在线观看| 国产精品不卡视频一区二区| 成人漫画全彩无遮挡| 午夜福利,免费看| 欧美精品av麻豆av| 在线观看免费视频网站a站| 精品午夜福利在线看| 乱码一卡2卡4卡精品| 日韩不卡一区二区三区视频在线| 自线自在国产av| 日本与韩国留学比较| a级毛色黄片| 亚洲精品,欧美精品| 在线免费观看不下载黄p国产| 交换朋友夫妻互换小说| 午夜激情av网站| 成人毛片a级毛片在线播放| 日本欧美国产在线视频| 永久免费av网站大全| a 毛片基地| 精品人妻偷拍中文字幕| 久久久久国产精品人妻一区二区| 国产成人精品久久久久久| 在现免费观看毛片| 欧美精品人与动牲交sv欧美| 亚洲国产精品成人久久小说| 日本免费在线观看一区| 女性被躁到高潮视频| 纵有疾风起免费观看全集完整版| 少妇人妻精品综合一区二区| 免费人妻精品一区二区三区视频| 国产黄色视频一区二区在线观看| 天天躁夜夜躁狠狠久久av| 亚洲欧美一区二区三区国产| 日韩一本色道免费dvd| 一边摸一边做爽爽视频免费| 91aial.com中文字幕在线观看| 欧美成人精品欧美一级黄| 日韩伦理黄色片| 精品人妻一区二区三区麻豆| 女人精品久久久久毛片| 亚洲国产毛片av蜜桃av| 国产男人的电影天堂91| 男女无遮挡免费网站观看| 午夜久久久在线观看| 亚洲人成网站在线观看播放| 成人漫画全彩无遮挡| 在线精品无人区一区二区三| 国产av国产精品国产| 免费大片黄手机在线观看| 亚洲,一卡二卡三卡| 女性生殖器流出的白浆| 91久久精品国产一区二区三区| 久久亚洲国产成人精品v| 成人18禁高潮啪啪吃奶动态图| 亚洲国产精品999| 十八禁高潮呻吟视频| 七月丁香在线播放| 免费少妇av软件| 日本猛色少妇xxxxx猛交久久| 国产精品久久久av美女十八| 国产精品人妻久久久影院| 免费人成在线观看视频色| 男女国产视频网站| kizo精华| 91在线精品国自产拍蜜月| 免费高清在线观看日韩| 99久国产av精品国产电影| 亚洲激情五月婷婷啪啪| 建设人人有责人人尽责人人享有的| 超色免费av| 18禁观看日本| 免费不卡的大黄色大毛片视频在线观看| 巨乳人妻的诱惑在线观看| 精品一区在线观看国产| 免费看光身美女| 黄色毛片三级朝国网站| 中文字幕另类日韩欧美亚洲嫩草| 亚洲综合精品二区| 母亲3免费完整高清在线观看 | 亚洲av福利一区| 亚洲国产看品久久| 久久午夜综合久久蜜桃| 九草在线视频观看|