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

    多目標(biāo)粒子群算法在地鐵速度曲線優(yōu)化中的應(yīng)用

    2022-09-06 08:55:10黃江平鐘曉靜
    關(guān)鍵詞:質(zhì)點(diǎn)阻力能耗

    黃江平,鐘曉靜

    (華東交通大學(xué)電氣與自動(dòng)化工程學(xué)院,南昌 330013)

    隨著城軌列車運(yùn)營里程的增加,城市軌道交通系統(tǒng)面臨的壓力也越來越大,要求最大程度地減少能源消耗,降低其運(yùn)營成本和對環(huán)境影響[1]。截止2019年底,全國城市軌道交通運(yùn)營里程從2064 km增長到6333.3 km。據(jù)統(tǒng)計(jì),我國城市軌道交通每年用電量約15.26億kW·h[2],同時(shí),根據(jù)現(xiàn)有城市軌道交通能耗數(shù)據(jù)顯示,列車牽引是列車耗能的主要因素,約為總能耗的50%[3]。因此,研究降低列車牽引能耗方法,對城市軌道交通建設(shè)具有重要的價(jià)值和意義。

    HOWLETT等[4-5]用最大值原理證明了在直線道路上最佳駕駛策略是4種模式的組合,包括牽引、巡航、惰行和制動(dòng);MIYATAKE等[6]討論了再生制動(dòng)的列車節(jié)能策略,比較了梯度法、二次規(guī)劃法的優(yōu)劣。智能算法的迅速發(fā)展也為列車運(yùn)行優(yōu)化提供了更好的選擇方案,段玉瓊等[7]提出基于人工蜂群算法與操縱工況序列相結(jié)合的方法來優(yōu)化計(jì)算,確定高速列車操縱工況關(guān)鍵轉(zhuǎn)換點(diǎn)最優(yōu)位置和速度。與此同時(shí),李旭陽等[8]應(yīng)用“四階段”最優(yōu)運(yùn)行策略,以牽引能耗值和運(yùn)行時(shí)間誤差的權(quán)重作為目標(biāo)函數(shù),采用遺傳算法進(jìn)行優(yōu)化;SHANG GUAN等[9]將時(shí)間和能耗作為2個(gè)優(yōu)化目標(biāo),混合差分和模擬退火2種進(jìn)化算法求解最優(yōu)運(yùn)行曲線。國內(nèi)外學(xué)者雖然都對列車駕駛策略和運(yùn)行曲線優(yōu)化進(jìn)行了深入研究,但在其選擇上還有一些不足。比如,駕駛策略通常采用長時(shí)間巡航運(yùn)行,以至在復(fù)雜路況下為保持勻速,列車需在“牽引-制動(dòng)”下反復(fù)切換,導(dǎo)致能耗過大;而常規(guī)處理多目標(biāo)問題的算法,通常是對多個(gè)目標(biāo)進(jìn)行加權(quán)求和,將其轉(zhuǎn)換為單目標(biāo)問題進(jìn)行處理,這可能導(dǎo)致最終結(jié)果波動(dòng)較大;同時(shí),為便于分析計(jì)算往往采用單質(zhì)點(diǎn)模型,不符合實(shí)際線路阻力變化情況。

    綜上所述,針對單質(zhì)點(diǎn)和多質(zhì)點(diǎn)建模的區(qū)別,同時(shí)考慮惰行工況對列車牽引能耗的影響,設(shè)置不同的控制策略加以比較,以列車運(yùn)行時(shí)間和能耗作為優(yōu)化目標(biāo),采用基于自適應(yīng)網(wǎng)格的多目標(biāo)粒子群算法,并對慣性權(quán)重采用線性微分遞減法處理,獲得Pareto[10]最優(yōu)解集后,再結(jié)合模糊隸屬度函數(shù)篩選最優(yōu)方案,并與常規(guī)粒子群算法進(jìn)行對比。最后,采用北京某地鐵線路數(shù)據(jù)對算法進(jìn)行驗(yàn)證。

    1 列車模型構(gòu)建

    1.1 模型分析

    為簡化分析,通常將列車看成有質(zhì)量無長度的單質(zhì)點(diǎn)[11],而單質(zhì)點(diǎn)模型在經(jīng)過變坡點(diǎn)或變曲率點(diǎn)時(shí)列車所受附加阻力會(huì)發(fā)生突變,這與實(shí)際的附加阻力變化情況不符。因此,在單質(zhì)點(diǎn)模型基礎(chǔ)上進(jìn)行改進(jìn),將列車視為多個(gè)質(zhì)點(diǎn)構(gòu)成的質(zhì)點(diǎn)鏈,將每節(jié)車廂視為一個(gè)質(zhì)點(diǎn),同時(shí)每節(jié)車廂之間為剛性連接[12],每節(jié)車廂的運(yùn)行狀態(tài)與整個(gè)列車的運(yùn)行狀態(tài)相同。

    由于多質(zhì)點(diǎn)模型未忽略列車長度,故在實(shí)際線路中列車阻力變化情況與單質(zhì)點(diǎn)模型不同。列車所受附加阻力中,隧道附加阻力只與隧道長度有關(guān),模型影響不大,因此,主要分析不同模型下坡道及曲線附加阻力的變化情況。圖1分別給出列車在不同模型下附加阻力的變化,仿真時(shí)以列車前端作為參照點(diǎn)。

    圖1 附加阻力變化

    通過仿真結(jié)果可知,列車在經(jīng)過變坡點(diǎn)及變曲率點(diǎn)的過程,坡道和曲線附加阻力是一個(gè)逐漸變化的過程,與單質(zhì)點(diǎn)模型的突變相比多質(zhì)點(diǎn)模型更貼合實(shí)際,減少誤差。

    1.2 列車運(yùn)動(dòng)方程

    對列車的運(yùn)動(dòng)進(jìn)行受力分析,通過牛頓第二定律推導(dǎo),得出列車運(yùn)動(dòng)方程。

    (1)

    (2)

    W(v)=(w1+w2+w3+w0)×G×M/1 000

    (3)

    w0=a+bv+cv2

    (4)

    式中,x為列車的位置;F牽為牽引力;B0為制動(dòng)力;W(v)為列車運(yùn)行阻力,w0為單位基本阻力,其中,a、b和c是與列車外觀和機(jī)械結(jié)構(gòu)相關(guān)的常數(shù),w1為單位坡道附加阻力,w2為單位曲線附加阻力,w3為單位隧道附加阻力;ηq為牽引力使用系數(shù);ηz為制動(dòng)力使用系數(shù);M為列車質(zhì)量;γ為回轉(zhuǎn)質(zhì)量系數(shù);G為重力加速度,通常取9.8 m/s2。

    2 優(yōu)化模型構(gòu)建

    列車優(yōu)化問題是研究列車在運(yùn)行過程中滿足運(yùn)行位移、限制速度、工況轉(zhuǎn)換等約束條件下,搜尋最佳工況切換點(diǎn)。通過優(yōu)化算法對建立的優(yōu)化模型進(jìn)行處理,獲得最優(yōu)的速度-距離曲線。

    2.1 目標(biāo)方程

    列車運(yùn)行狀態(tài)通常有4種:牽引,巡航,惰行,制動(dòng)。在此基礎(chǔ)上對牽引總能耗E和運(yùn)行時(shí)間T兩個(gè)優(yōu)化目標(biāo)建立關(guān)系式。

    (5)

    (6)

    式中,s為兩站間的區(qū)間長度;x為列車運(yùn)行時(shí)的位置;δ為再生制動(dòng)反饋能量的效率。

    2.2 評價(jià)指標(biāo)

    在保證列車準(zhǔn)時(shí)性的同時(shí),力求降低列車運(yùn)行過程牽引能耗,將fE和fT作為評價(jià)指標(biāo)。

    fE=E

    (7)

    fT=T

    (8)

    式中,E為列車在整個(gè)區(qū)間的牽引能耗;T為列車實(shí)際運(yùn)行時(shí)間。

    2.3 目標(biāo)及約束條件

    根據(jù)上述目標(biāo)方程和評價(jià)指標(biāo)公式,建立多目標(biāo)優(yōu)化數(shù)學(xué)模型如下。

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

    F=min(fE,fT)

    (9)

    約束條件

    (10)

    式中,V0和Vd分別為列車起始速度及終止速度;Vlim為區(qū)間線路的限制速度;X為列車最終運(yùn)行位移,x為線路實(shí)際長度;T0為列車計(jì)劃運(yùn)行時(shí)間。

    2.4 駕駛策略選擇

    傳統(tǒng)的ATO駕駛策略在牽引過程后,列車通常將進(jìn)入巡航工況,而在復(fù)雜的短距離路段中,長時(shí)間保持勻速會(huì)使得列車在牽引和制動(dòng)之間反復(fù)切換,損耗很大[13]。因此,用牽引和惰行工況來代替巡航工況可有效降低列車牽引能耗。如圖2所示,將駕駛策略分為未加入惰行、加入一次惰行及加入多次惰行3種控制策略。

    圖2 列車控制策略

    3 速度曲線優(yōu)化算法

    采用一種基于自適應(yīng)網(wǎng)格的多目標(biāo)粒子群算法(MOPSO based on adaptive grid algorithm,簡稱AGA-MOPSO)[14],將Pareto解集的概念引入到算法中,對每次迭代產(chǎn)生的Pareto最優(yōu)粒子進(jìn)行存檔[15]。與將多目標(biāo)加權(quán)求和轉(zhuǎn)化為單目標(biāo)的傳統(tǒng)粒子群算法相比,該算法無論是在可行解的多樣性及穩(wěn)定性上,還是最后的優(yōu)化效果上都更加出色。

    3.1 基于自適應(yīng)網(wǎng)格的MOPSO算法

    算法求解的具體步驟如下。

    (1)設(shè)置算法參數(shù)。種群數(shù)量M,外部存儲(chǔ)數(shù)量N,加速度常量C1和C2,慣性權(quán)重Wmax,Wmin及粒子的維度(即工況轉(zhuǎn)換點(diǎn)的個(gè)數(shù))。

    (2)初始化粒子種群,包括粒子群的位置pop[i](列車運(yùn)行可能的工況轉(zhuǎn)換點(diǎn)),速度V[i](工況轉(zhuǎn)換點(diǎn)位置改變的量)。

    (3)計(jì)算粒子目標(biāo)函數(shù)值,確定粒子個(gè)體最優(yōu)位置pbest[i]。構(gòu)建自適應(yīng)網(wǎng)格,通過Archive集中粒子的密度信息來選擇出全局最優(yōu)粒子gbest[i]。

    (4)更新粒子群中粒子位置和速度公式。

    V[i]=W×V[i]+rand( )×C1×(pbest[i]-

    pop[i])+rand( )×C2×(gbest[i]-pop[i])

    (11)

    pop[i]=pop[i]+V[i]

    (12)

    式中,rand( )為范圍[0,1]內(nèi)的隨機(jī)數(shù)。

    (5)計(jì)算各個(gè)粒子的評價(jià)指標(biāo),利用支配關(guān)系更新個(gè)體最優(yōu)粒子pbest[i]。

    (6)對Archive集進(jìn)行更新,同時(shí)對Archive集進(jìn)行截?cái)?,保留N個(gè)最優(yōu)解。

    (7)判斷是否滿足設(shè)定的迭代次數(shù),若否,則跳出程序,否則跳至第三步。若是,則結(jié)束迭代,輸出Archive解集。

    3.1.1 自適應(yīng)網(wǎng)格的構(gòu)建

    把目標(biāo)空間E-T(能耗-時(shí)間)以M×M的方式劃分為M2個(gè)網(wǎng)格,每個(gè)網(wǎng)格的長寬計(jì)算如下

    (13)

    (14)

    式中,M為Archive解集在當(dāng)前迭代中解的個(gè)數(shù);maxEt,minEt分別為迭代t次時(shí)Archive解集中列車運(yùn)行能耗的最大和最小值;maxTt和minTt分別為迭代t次時(shí)列車運(yùn)行時(shí)間的最大和最小值。

    統(tǒng)計(jì)每個(gè)網(wǎng)格中粒子數(shù)量,將一個(gè)網(wǎng)格內(nèi)解的個(gè)數(shù)作為粒子分布密度。圖3給出了網(wǎng)格構(gòu)建示意。

    圖3 自適應(yīng)網(wǎng)格構(gòu)建

    3.1.2 Archive集的截?cái)?/p>

    為維持Archive集中解的數(shù)量,保證解分布的多樣性及算法的全局搜索能力,要對外部存檔進(jìn)行刪除。通過式(15)來計(jì)算粒子數(shù)量大于1的網(wǎng)格中要隨機(jī)刪除的粒子個(gè)數(shù)P。

    (15)

    式中,Int( )為取整函數(shù);A為Archive集中的粒子數(shù);N為要保留最優(yōu)解的個(gè)數(shù);Di為網(wǎng)格i的粒子個(gè)數(shù)。

    3.2 慣性權(quán)重W

    慣性權(quán)重W的變化對算法搜索能力有很大影響。與傳統(tǒng)的線性遞減法相比,采用線性微分遞減法可以快速地收斂逼近于全局最優(yōu)值,提高運(yùn)算效率。W計(jì)算公式如下

    (16)

    式中,Wmax為慣性權(quán)重最大值;Wmin為慣性權(quán)重最小值;ta為總迭代次數(shù);t為當(dāng)前迭代次數(shù)。

    3.3 Pareto解集中最優(yōu)方案選擇

    采用模糊隸屬度函數(shù)[16]對得到的Pareto解集中的每個(gè)解進(jìn)行評價(jià),以此來選出所有解中的最優(yōu)方案。模糊隸屬度函數(shù)如下

    (17)

    (18)

    式中,UE,UT分別為能耗和準(zhǔn)點(diǎn)誤差的隸屬度值;Emax和Emin分別為Pareto解集中能耗的最大、最小值;Ei為第i個(gè)解的能耗;T0為計(jì)劃運(yùn)行時(shí)間;Ti為第i個(gè)解的運(yùn)行時(shí)間;Te為距離計(jì)劃時(shí)間的誤差。最優(yōu)駕駛方案為UE與UT的和為最小值所對應(yīng)的Pareto解。

    4 仿真與分析

    4.1 仿真參數(shù)設(shè)置

    為驗(yàn)證算法的有效性,以北京地鐵常用的B型列車為研究對象[17],選取北京某段地鐵線路進(jìn)行仿真研究。算法參數(shù)和線路基本參數(shù)分別如表1、表2所示,區(qū)間長2 810 m,最大坡度7‰,最大限速80 km/h,列車計(jì)劃運(yùn)行時(shí)間為190 s,列車實(shí)際運(yùn)行時(shí)間為180~200 s,誤差為±10 s。

    表1 算法參數(shù)

    4.2 兩種算法分析對比

    將傳統(tǒng)PSO算法和AGA-MOPSO算法在相同駕駛策略下進(jìn)行仿真實(shí)驗(yàn)。兩種算法獲得的解如圖4所示。其中,PSO算法通過加權(quán)求和,將其轉(zhuǎn)換為單目標(biāo)優(yōu)化問題,獲得唯一最優(yōu)解;而AGA-MOPSO引入Pareto解原理,獲得一組由時(shí)間-能耗平衡的Pareto解集。

    表2 線路參數(shù)

    通過模糊隸屬度函數(shù)對Pareto解集中的所有解進(jìn)行評價(jià),選出最優(yōu)方案,如表3所示。與PSO相比,AGA-MOPSO獲得的最優(yōu)解在更加準(zhǔn)時(shí)的同時(shí)可節(jié)能達(dá)7.4%。

    圖4 兩種算法的解

    表3 最優(yōu)方案對比

    分別將兩種算法獨(dú)立運(yùn)行20次,得到了每次優(yōu)化獲得的最優(yōu)解的時(shí)間-能耗,如圖5所示。

    為比較兩組數(shù)據(jù)的穩(wěn)定程度,采用平均值、標(biāo)準(zhǔn)差及方差作為評價(jià)指標(biāo),分別對兩種算法下的時(shí)間誤差-能耗兩個(gè)變量進(jìn)行比較。

    如表4所示,兩種算法時(shí)間誤差波動(dòng)范圍為2.12 s和2.76 s,相對于計(jì)劃運(yùn)行時(shí)間190 s來說影響不大,因此,時(shí)間誤差波動(dòng)可忽略,主要對能耗變量的波動(dòng)進(jìn)行描述。其中,PSO算法得到的能耗波動(dòng)范圍為2.02 kW·h,平均值為13.89 kW·h,方差為0.497;與之相比AGA-MOPSO算法的能耗波動(dòng)范圍更小為1.40 kW·h,平均能耗為12.48 kW·h,方差為0.137。因此,AGA-MOPSO算法獲得的最終解波動(dòng)范圍更小,且更加穩(wěn)定。

    圖5 算法穩(wěn)定性比較

    表4 兩個(gè)變量評價(jià)指標(biāo)描述統(tǒng)計(jì)

    綜上所述,對比兩種算法,AGA-MOPSO算法優(yōu)化得到最優(yōu)解的效果比PSO算法更好,更加準(zhǔn)時(shí)且能耗更小。同時(shí),無論是解的多樣性還是穩(wěn)定性等方面都優(yōu)于PSO算法。

    4.3 不同控制策略對比

    針對未加入惰行、加入一次惰行以及加入多次惰行3種不同控制策略,在利用AGA-MOPSO算法獲得各自的Pareto解集后,采用上述模糊隸屬度函數(shù)篩選出最優(yōu)方案。表5為各種駕駛策略對應(yīng)的最優(yōu)駕駛方案。

    表5 不同控制策略下對應(yīng)的最優(yōu)方案

    3種控制策略的最優(yōu)方案所對應(yīng)的速度-距離曲線如圖6所示。

    圖6 不同策略下優(yōu)化曲線

    在滿足列車準(zhǔn)點(diǎn)要求的前提下,多次加入惰行的控制策略運(yùn)行時(shí)間為190.3 s,能耗12.47 kW·h為3種駕駛策略中耗能最低,與策略未加入惰行相比節(jié)能18.3%,與策略加入一次惰行相比節(jié)能13.2%。綜合以上仿真結(jié)果可以得出,策略加入多次惰行的節(jié)能準(zhǔn)點(diǎn)的優(yōu)化效果最好。

    5 結(jié)論

    針對城市軌道交通中列車速度曲線優(yōu)化問題,建立了多質(zhì)點(diǎn)列車模型,設(shè)計(jì)以降低牽引能耗和提高列車準(zhǔn)點(diǎn)的多目標(biāo)優(yōu)化模型,采用AGA-MOPSO多目標(biāo)優(yōu)化算法對不同駕駛策略進(jìn)行優(yōu)化,生成一組能耗-時(shí)間Pareto解,且通過模糊隸屬度函數(shù)進(jìn)行篩選,獲得最優(yōu)運(yùn)行方案,得出以下結(jié)論。

    (1)多質(zhì)點(diǎn)模型下附加阻力的變化更加符合實(shí)際情況,與單質(zhì)點(diǎn)模型相比減小了列車在變坡點(diǎn)和變曲率點(diǎn)的附加阻力誤差。

    (2)對比傳統(tǒng)PSO算法處理多目標(biāo)問題,AGA-MOPSO多目標(biāo)優(yōu)化算法在解方案的多樣性上優(yōu)于PSO算法,優(yōu)化得到的最優(yōu)方案效果更好,在準(zhǔn)點(diǎn)誤差0.3s的同時(shí),節(jié)能可達(dá)7.4%。同時(shí),在最優(yōu)解穩(wěn)定性上也優(yōu)于PSO算法。

    (3)對于3種不同的駕駛策略,采用模糊隸屬函數(shù)從Pareto解集中選出3種駕駛策略所對應(yīng)的最優(yōu)方案。對比可知,在滿足列車準(zhǔn)點(diǎn)的前提下,加入多次惰行控制的駕駛策略優(yōu)化效果更好。

    本次優(yōu)化結(jié)果為地鐵列車運(yùn)行方案提供了一定思路,但只考慮了列車準(zhǔn)點(diǎn)以及能耗兩個(gè)目標(biāo)變量,下一步研究可考慮加入舒適度、電分相等其他因素。

    猜你喜歡
    質(zhì)點(diǎn)阻力能耗
    120t轉(zhuǎn)爐降低工序能耗生產(chǎn)實(shí)踐
    昆鋼科技(2022年2期)2022-07-08 06:36:14
    能耗雙控下,漲價(jià)潮再度來襲!
    鼻阻力測定在兒童OSA診療中的臨床作用
    巧用“搬運(yùn)法”解決連續(xù)質(zhì)點(diǎn)模型的做功問題
    探討如何設(shè)計(jì)零能耗住宅
    零阻力
    英語文摘(2020年10期)2020-11-26 08:12:12
    日本先進(jìn)的“零能耗住宅”
    別讓摩擦成為學(xué)習(xí)的阻力
    質(zhì)點(diǎn)的直線運(yùn)動(dòng)
    質(zhì)點(diǎn)的直線運(yùn)動(dòng)
    午夜福利乱码中文字幕| 一区二区日韩欧美中文字幕| av欧美777| 飞空精品影院首页| 欧美在线黄色| 国产三级黄色录像| 久久九九热精品免费| 久久久精品国产亚洲av高清涩受| 精品国产乱子伦一区二区三区 | 精品少妇黑人巨大在线播放| 国产成人精品无人区| 人人妻人人澡人人看| 制服人妻中文乱码| 亚洲欧美精品综合一区二区三区| av在线app专区| 亚洲久久久国产精品| 亚洲七黄色美女视频| 大香蕉久久成人网| 欧美精品av麻豆av| 久久精品国产综合久久久| 老司机在亚洲福利影院| 人妻 亚洲 视频| 丝袜美腿诱惑在线| 1024视频免费在线观看| 国产亚洲午夜精品一区二区久久| 久久 成人 亚洲| 国产精品一区二区精品视频观看| 天天影视国产精品| 岛国毛片在线播放| av福利片在线| 精品国内亚洲2022精品成人 | 搡老乐熟女国产| 黑人操中国人逼视频| 亚洲熟女毛片儿| 丁香六月天网| 亚洲国产精品一区二区三区在线| 欧美日韩一级在线毛片| 一区二区三区四区激情视频| 在线亚洲精品国产二区图片欧美| 亚洲成人国产一区在线观看| 亚洲五月婷婷丁香| 国产精品一区二区在线不卡| 亚洲中文字幕日韩| 久久天堂一区二区三区四区| 电影成人av| 亚洲成国产人片在线观看| 色播在线永久视频| 亚洲少妇的诱惑av| 国产免费现黄频在线看| 国产国语露脸激情在线看| 欧美大码av| 亚洲一区中文字幕在线| 亚洲国产看品久久| 一级,二级,三级黄色视频| 老司机午夜福利在线观看视频 | 在线 av 中文字幕| 国产日韩欧美亚洲二区| 黄色视频,在线免费观看| 国产欧美日韩综合在线一区二区| 日日摸夜夜添夜夜添小说| 亚洲精品成人av观看孕妇| 啦啦啦免费观看视频1| 叶爱在线成人免费视频播放| 极品人妻少妇av视频| 人人妻人人添人人爽欧美一区卜| 国产无遮挡羞羞视频在线观看| 国产精品香港三级国产av潘金莲| 亚洲一码二码三码区别大吗| 精品久久蜜臀av无| 国产精品一区二区精品视频观看| 久久久欧美国产精品| 国产一区二区三区在线臀色熟女 | 久久久国产欧美日韩av| 99精品欧美一区二区三区四区| 久久天躁狠狠躁夜夜2o2o| 日韩大片免费观看网站| 欧美日韩国产mv在线观看视频| 国产精品国产三级国产专区5o| 久久亚洲精品不卡| 国产精品 国内视频| 国产激情久久老熟女| 欧美性长视频在线观看| h视频一区二区三区| 亚洲九九香蕉| 久久久久视频综合| 在线 av 中文字幕| 考比视频在线观看| 美女大奶头黄色视频| 成人国语在线视频| 91精品伊人久久大香线蕉| 久热爱精品视频在线9| 成年动漫av网址| 大型av网站在线播放| 亚洲综合色网址| 在线十欧美十亚洲十日本专区| 丁香六月天网| 在线看a的网站| 亚洲伊人久久精品综合| 欧美日韩国产mv在线观看视频| 女人被躁到高潮嗷嗷叫费观| 国产亚洲av片在线观看秒播厂| 亚洲人成77777在线视频| 十八禁人妻一区二区| 亚洲欧洲精品一区二区精品久久久| 亚洲精品一卡2卡三卡4卡5卡 | 久久久久久亚洲精品国产蜜桃av| 色播在线永久视频| 男人添女人高潮全过程视频| 日韩大码丰满熟妇| 黄色a级毛片大全视频| 啦啦啦啦在线视频资源| 精品人妻熟女毛片av久久网站| 中亚洲国语对白在线视频| 国产精品欧美亚洲77777| 9热在线视频观看99| 久久精品国产亚洲av香蕉五月 | 久久女婷五月综合色啪小说| 午夜老司机福利片| 亚洲国产精品成人久久小说| 午夜日韩欧美国产| 韩国精品一区二区三区| 久久中文字幕一级| 首页视频小说图片口味搜索| 欧美日韩精品网址| 18禁国产床啪视频网站| 久久久国产成人免费| 国产成人免费无遮挡视频| 自线自在国产av| 成人亚洲精品一区在线观看| 午夜激情久久久久久久| 国产精品一二三区在线看| 天堂俺去俺来也www色官网| 99国产极品粉嫩在线观看| 亚洲国产欧美日韩在线播放| av在线老鸭窝| 男女床上黄色一级片免费看| 国产色视频综合| 人人妻人人添人人爽欧美一区卜| 日韩制服丝袜自拍偷拍| 叶爱在线成人免费视频播放| 午夜久久久在线观看| 亚洲欧美日韩另类电影网站| 亚洲精品国产精品久久久不卡| 女警被强在线播放| 亚洲va日本ⅴa欧美va伊人久久 | 十分钟在线观看高清视频www| 丰满少妇做爰视频| 国产日韩欧美在线精品| 麻豆av在线久日| 亚洲国产精品一区三区| 国产成人欧美| 99精国产麻豆久久婷婷| 亚洲精品国产一区二区精华液| 午夜激情av网站| 欧美午夜高清在线| 亚洲国产av新网站| 一边摸一边做爽爽视频免费| 男女床上黄色一级片免费看| 少妇人妻久久综合中文| 在线十欧美十亚洲十日本专区| 精品人妻一区二区三区麻豆| 99国产精品一区二区蜜桃av | 国产成人免费观看mmmm| 9热在线视频观看99| 他把我摸到了高潮在线观看 | 久久久精品94久久精品| 777久久人妻少妇嫩草av网站| 一区二区日韩欧美中文字幕| 波多野结衣一区麻豆| 丝袜美足系列| 亚洲国产看品久久| 麻豆国产av国片精品| 中文字幕高清在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 麻豆av在线久日| 国产精品一区二区免费欧美 | 国产精品亚洲av一区麻豆| 一区二区三区乱码不卡18| 精品亚洲成国产av| 一级,二级,三级黄色视频| 91老司机精品| 丁香六月天网| 国产精品偷伦视频观看了| 久久久久久久精品精品| 日本欧美视频一区| 免费高清在线观看视频在线观看| 少妇猛男粗大的猛烈进出视频| 国产野战对白在线观看| 国产成人精品无人区| 人人澡人人妻人| 国产亚洲av片在线观看秒播厂| 欧美精品一区二区大全| 亚洲色图综合在线观看| 乱人伦中国视频| 51午夜福利影视在线观看| 免费久久久久久久精品成人欧美视频| 又紧又爽又黄一区二区| 色94色欧美一区二区| 婷婷丁香在线五月| 水蜜桃什么品种好| 老司机影院毛片| 久久久久国产一级毛片高清牌| www日本在线高清视频| 在线观看一区二区三区激情| 1024视频免费在线观看| 脱女人内裤的视频| 国产精品久久久av美女十八| 波多野结衣一区麻豆| netflix在线观看网站| 99国产极品粉嫩在线观看| 最近最新中文字幕大全免费视频| 国产片内射在线| 欧美精品av麻豆av| 久久久精品94久久精品| 一进一出抽搐动态| 国产免费av片在线观看野外av| 中文字幕人妻丝袜一区二区| 视频区欧美日本亚洲| 天天操日日干夜夜撸| 久久精品久久久久久噜噜老黄| 99精品久久久久人妻精品| 桃红色精品国产亚洲av| 少妇人妻久久综合中文| 亚洲美女黄色视频免费看| 久久午夜综合久久蜜桃| av超薄肉色丝袜交足视频| 久久人人爽人人片av| 午夜激情av网站| 午夜福利视频精品| 午夜福利视频精品| 99久久综合免费| 中文字幕高清在线视频| 国产精品一区二区免费欧美 | 午夜视频精品福利| 日日夜夜操网爽| 国产成人精品久久二区二区91| 黄色 视频免费看| 丰满人妻熟妇乱又伦精品不卡| 国产色视频综合| 90打野战视频偷拍视频| 国产又色又爽无遮挡免| 午夜福利一区二区在线看| 国产欧美日韩一区二区三区在线| 又紧又爽又黄一区二区| 亚洲精品中文字幕在线视频| 自线自在国产av| avwww免费| 久久av网站| 亚洲欧美精品自产自拍| 亚洲精品日韩在线中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 成人18禁高潮啪啪吃奶动态图| 亚洲国产精品一区三区| 欧美人与性动交α欧美软件| 精品卡一卡二卡四卡免费| 日本五十路高清| videosex国产| 国产成人欧美在线观看 | 91老司机精品| 搡老岳熟女国产| 亚洲激情五月婷婷啪啪| 久久青草综合色| 丁香六月天网| 两性夫妻黄色片| 9191精品国产免费久久| 久久久久精品人妻al黑| 男女高潮啪啪啪动态图| 国产日韩欧美在线精品| 久久这里只有精品19| 亚洲五月色婷婷综合| 精品国产一区二区三区久久久樱花| 丰满人妻熟妇乱又伦精品不卡| 男女午夜视频在线观看| 欧美日韩福利视频一区二区| 亚洲中文字幕日韩| 伊人亚洲综合成人网| 亚洲中文日韩欧美视频| 欧美黑人欧美精品刺激| 男女高潮啪啪啪动态图| 日韩精品免费视频一区二区三区| 女人精品久久久久毛片| 18禁观看日本| 制服诱惑二区| 秋霞在线观看毛片| 国产成人精品无人区| 免费一级毛片在线播放高清视频 | 最近最新免费中文字幕在线| tocl精华| 欧美xxⅹ黑人| 丝袜脚勾引网站| 久久久久久久久免费视频了| av片东京热男人的天堂| 久久精品久久久久久噜噜老黄| 国产成人免费观看mmmm| 婷婷丁香在线五月| 人人妻人人澡人人爽人人夜夜| 欧美激情极品国产一区二区三区| 一个人免费看片子| 久久久久久免费高清国产稀缺| 久久人人爽人人片av| 国产免费av片在线观看野外av| 国产国语露脸激情在线看| 欧美乱码精品一区二区三区| 老汉色av国产亚洲站长工具| 中国国产av一级| 国产又爽黄色视频| 久久久久久亚洲精品国产蜜桃av| 欧美+亚洲+日韩+国产| 欧美日韩成人在线一区二区| 亚洲国产精品成人久久小说| 男女之事视频高清在线观看| 高清欧美精品videossex| 成年人免费黄色播放视频| 成年女人毛片免费观看观看9 | 首页视频小说图片口味搜索| 久久精品国产a三级三级三级| 亚洲五月婷婷丁香| 亚洲精品成人av观看孕妇| 中文字幕另类日韩欧美亚洲嫩草| 色播在线永久视频| 国产亚洲欧美在线一区二区| 精品国产乱子伦一区二区三区 | 啦啦啦中文免费视频观看日本| 考比视频在线观看| 午夜激情久久久久久久| 国产精品久久久av美女十八| 国产精品亚洲av一区麻豆| 亚洲国产成人一精品久久久| 制服诱惑二区| 丝袜在线中文字幕| 欧美日韩视频精品一区| 精品国产超薄肉色丝袜足j| 日韩人妻精品一区2区三区| 欧美人与性动交α欧美软件| 性色av乱码一区二区三区2| 丰满少妇做爰视频| 五月开心婷婷网| 淫妇啪啪啪对白视频 | 狂野欧美激情性xxxx| 好男人电影高清在线观看| 欧美成狂野欧美在线观看| 一级a爱视频在线免费观看| 亚洲欧美一区二区三区黑人| 亚洲专区字幕在线| 欧美xxⅹ黑人| 女警被强在线播放| 久久精品成人免费网站| 亚洲精品国产av蜜桃| 亚洲伊人色综图| 少妇精品久久久久久久| 丝袜喷水一区| 国产精品久久久av美女十八| 最近中文字幕2019免费版| 国产极品粉嫩免费观看在线| 欧美精品av麻豆av| 黄色视频在线播放观看不卡| 男女无遮挡免费网站观看| 99国产综合亚洲精品| 成人国语在线视频| 亚洲va日本ⅴa欧美va伊人久久 | 视频区欧美日本亚洲| 麻豆国产av国片精品| 久久久久精品国产欧美久久久 | 俄罗斯特黄特色一大片| 99九九在线精品视频| 久久久水蜜桃国产精品网| 丰满迷人的少妇在线观看| 国产精品香港三级国产av潘金莲| av网站在线播放免费| 国产又爽黄色视频| 国产成人啪精品午夜网站| 美女午夜性视频免费| 日韩中文字幕视频在线看片| 国产成人系列免费观看| 91麻豆av在线| 精品视频人人做人人爽| 欧美中文综合在线视频| 人妻 亚洲 视频| 日本猛色少妇xxxxx猛交久久| 一区福利在线观看| 国产精品 欧美亚洲| 日韩欧美一区视频在线观看| 国产精品欧美亚洲77777| 男女无遮挡免费网站观看| 熟女少妇亚洲综合色aaa.| 巨乳人妻的诱惑在线观看| 亚洲精品乱久久久久久| 啦啦啦免费观看视频1| 各种免费的搞黄视频| 丁香六月天网| 久久久国产一区二区| 啦啦啦免费观看视频1| 成人三级做爰电影| 欧美黑人精品巨大| 亚洲精品国产色婷婷电影| 国产免费av片在线观看野外av| av一本久久久久| 日韩欧美一区二区三区在线观看 | 亚洲精品日韩在线中文字幕| 啦啦啦啦在线视频资源| 免费观看a级毛片全部| 亚洲国产成人一精品久久久| 少妇人妻久久综合中文| 国产又爽黄色视频| 国产片内射在线| 欧美午夜高清在线| 9色porny在线观看| av有码第一页| 无限看片的www在线观看| 51午夜福利影视在线观看| 免费在线观看视频国产中文字幕亚洲 | 最新的欧美精品一区二区| 97精品久久久久久久久久精品| 中国美女看黄片| 秋霞在线观看毛片| 久久久国产成人免费| 欧美日韩亚洲高清精品| 热99re8久久精品国产| 日韩欧美国产一区二区入口| 日韩中文字幕视频在线看片| av视频免费观看在线观看| 黑人猛操日本美女一级片| 考比视频在线观看| 亚洲精品久久午夜乱码| 国产亚洲欧美在线一区二区| 夜夜骑夜夜射夜夜干| 性色av乱码一区二区三区2| 一边摸一边抽搐一进一出视频| 欧美精品一区二区大全| av免费在线观看网站| 亚洲精品自拍成人| 啦啦啦啦在线视频资源| 国产精品秋霞免费鲁丝片| 亚洲少妇的诱惑av| 91麻豆av在线| 欧美变态另类bdsm刘玥| 久久久久久久久免费视频了| 久久精品人人爽人人爽视色| 国产成+人综合+亚洲专区| 69av精品久久久久久 | 国产主播在线观看一区二区| 搡老乐熟女国产| 国产精品二区激情视频| 在线观看免费高清a一片| 国产男女超爽视频在线观看| 午夜福利影视在线免费观看| 亚洲欧美精品自产自拍| 国产精品一区二区精品视频观看| 黄色 视频免费看| 18禁观看日本| 操出白浆在线播放| 亚洲熟女毛片儿| 交换朋友夫妻互换小说| 啦啦啦啦在线视频资源| 麻豆国产av国片精品| 国产日韩欧美视频二区| 麻豆av在线久日| 动漫黄色视频在线观看| 桃红色精品国产亚洲av| 自线自在国产av| 成年女人毛片免费观看观看9 | 亚洲精品成人av观看孕妇| 麻豆国产av国片精品| 精品人妻1区二区| 2018国产大陆天天弄谢| 国产不卡av网站在线观看| 午夜久久久在线观看| 精品国产一区二区三区四区第35| 国产欧美日韩综合在线一区二区| 精品福利观看| 成年美女黄网站色视频大全免费| 人妻人人澡人人爽人人| 国产精品香港三级国产av潘金莲| 国产精品成人在线| 精品卡一卡二卡四卡免费| 欧美午夜高清在线| 免费久久久久久久精品成人欧美视频| 丝袜在线中文字幕| 黄片大片在线免费观看| 久久 成人 亚洲| 国产淫语在线视频| 香蕉国产在线看| 大片免费播放器 马上看| 男女无遮挡免费网站观看| 成人手机av| 欧美日韩国产mv在线观看视频| h视频一区二区三区| 国产极品粉嫩免费观看在线| a 毛片基地| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲午夜精品一区二区久久| 一个人免费看片子| 亚洲成国产人片在线观看| 女人高潮潮喷娇喘18禁视频| av福利片在线| 18禁观看日本| 一区二区av电影网| 成年av动漫网址| 人人妻人人澡人人看| 18禁观看日本| 热99re8久久精品国产| 国产精品影院久久| 久久久精品国产亚洲av高清涩受| 日韩欧美国产一区二区入口| 菩萨蛮人人尽说江南好唐韦庄| 极品少妇高潮喷水抽搐| 久久九九热精品免费| av电影中文网址| 美女脱内裤让男人舔精品视频| 正在播放国产对白刺激| 女人高潮潮喷娇喘18禁视频| 亚洲五月婷婷丁香| 五月天丁香电影| 中文字幕另类日韩欧美亚洲嫩草| netflix在线观看网站| 老汉色av国产亚洲站长工具| 欧美国产精品一级二级三级| 日韩电影二区| av又黄又爽大尺度在线免费看| 国产精品麻豆人妻色哟哟久久| 亚洲成国产人片在线观看| 欧美精品高潮呻吟av久久| 最新在线观看一区二区三区| 丝袜喷水一区| 欧美一级毛片孕妇| 国产精品久久久av美女十八| 大片电影免费在线观看免费| 精品视频人人做人人爽| 制服人妻中文乱码| 久久天堂一区二区三区四区| 成人18禁高潮啪啪吃奶动态图| 深夜精品福利| 午夜免费鲁丝| 亚洲精品美女久久久久99蜜臀| 91老司机精品| 国产欧美亚洲国产| netflix在线观看网站| 香蕉丝袜av| 亚洲一区中文字幕在线| 人成视频在线观看免费观看| 亚洲国产精品成人久久小说| 一级,二级,三级黄色视频| 日韩一区二区三区影片| 两性夫妻黄色片| 亚洲少妇的诱惑av| 伦理电影免费视频| 一个人免费在线观看的高清视频 | 香蕉丝袜av| 超碰97精品在线观看| 一本—道久久a久久精品蜜桃钙片| 国产一卡二卡三卡精品| 亚洲专区字幕在线| 一区在线观看完整版| 搡老岳熟女国产| 啦啦啦中文免费视频观看日本| 人成视频在线观看免费观看| 伊人久久大香线蕉亚洲五| 国产激情久久老熟女| 久久久国产精品麻豆| 日韩 亚洲 欧美在线| 少妇人妻久久综合中文| 人妻一区二区av| 中文字幕精品免费在线观看视频| 亚洲熟女精品中文字幕| 在线观看免费视频网站a站| 日日爽夜夜爽网站| 亚洲精品自拍成人| 亚洲成人免费av在线播放| 国产片内射在线| 香蕉丝袜av| 亚洲国产欧美日韩在线播放| 一级片'在线观看视频| 日本黄色日本黄色录像| 五月天丁香电影| 一进一出抽搐动态| 午夜福利一区二区在线看| 大片免费播放器 马上看| 狠狠狠狠99中文字幕| 一边摸一边做爽爽视频免费| 日韩制服丝袜自拍偷拍| 又大又爽又粗| 男人爽女人下面视频在线观看| 中文字幕人妻丝袜制服| 欧美国产精品va在线观看不卡| 免费女性裸体啪啪无遮挡网站| 国产精品一区二区精品视频观看| 国产成人影院久久av| 99热全是精品| 国产色视频综合| 自线自在国产av| 久久久久久亚洲精品国产蜜桃av| 青春草视频在线免费观看| av天堂久久9| 日韩三级视频一区二区三区| 成人亚洲精品一区在线观看| 黑人猛操日本美女一级片| 他把我摸到了高潮在线观看 | tocl精华| 久久久久久久国产电影| 欧美黄色淫秽网站| 在线观看免费视频网站a站| 国产精品一区二区免费欧美 | 国产成人精品久久二区二区免费| 亚洲欧美一区二区三区黑人| 国产亚洲一区二区精品| 国产一区二区三区av在线| 黑人巨大精品欧美一区二区蜜桃| av片东京热男人的天堂| 日韩欧美一区视频在线观看| av天堂久久9| 捣出白浆h1v1| 夜夜骑夜夜射夜夜干| 电影成人av| 婷婷成人精品国产| 精品一区二区三区四区五区乱码|