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

    基于模糊廣義預(yù)測控制算法的ATO系統(tǒng)仿真研究

    2014-04-07 08:01:33顧桂梅蒲松芝
    城市軌道交通研究 2014年1期
    關(guān)鍵詞:坡道列車運(yùn)行模糊控制

    顧桂梅蒲松芝

    (1.蘭州交通大學(xué)自動化與電氣工程學(xué)院,730070,蘭州;2.卡斯柯信號有限公司,200070,上?!蔚谝蛔髡?,副教授)

    基于模糊廣義預(yù)測控制算法的ATO系統(tǒng)仿真研究

    顧桂梅1蒲松芝2

    (1.蘭州交通大學(xué)自動化與電氣工程學(xué)院,730070,蘭州;2.卡斯柯信號有限公司,200070,上?!蔚谝蛔髡?,副教授)

    針對ATO(列車自動運(yùn)行)系統(tǒng)的特點(diǎn),設(shè)計(jì)了基于模糊廣義預(yù)測控制算法的速度控制器,并結(jié)合實(shí)際線路建立了線路模型,生成了優(yōu)化的目標(biāo)曲線。利用Matlab分別對基于模糊廣義預(yù)測控制算法和基于隱式廣義預(yù)測控制算法的速度控制器進(jìn)行了仿真。仿真結(jié)果表明,基于模糊廣義預(yù)測控制算法的速度控制器不僅能使ATO系統(tǒng)的性能得到優(yōu)化,而且在一定程度上削弱了坡道對列車運(yùn)行速度造成的影響,具備良好的可移植性。

    列車自動運(yùn)行;模糊廣義預(yù)測控制;仿真

    First-author's address School of Automation and Electrical Engineering,Lanzhou Jiaotong University,730070,Lanzhou,China

    列車自動運(yùn)行(ATO)系統(tǒng)不僅能自動控制列車運(yùn)行,減輕司機(jī)負(fù)擔(dān),而且能在保證列車安全運(yùn)行的前提下,提高列車運(yùn)行的各項(xiàng)性能,因此得到廣泛關(guān)注,針對ATO系統(tǒng)的研究也取得了不少成果。然而,由于列車運(yùn)行過程十分復(fù)雜且極易受環(huán)境因素影響,因此很難找到精確的數(shù)學(xué)模型。模糊廣義預(yù)測控制(FGPC)算法是一種將模糊邏輯和廣義預(yù)測控制(GPC)算法相集成的控制策略。該算法對被控過程模型要求不高,且能進(jìn)一步提高系統(tǒng)的穩(wěn)態(tài)精度,減小動態(tài)誤差[1]。為此,本文將FGPC算法應(yīng)用到ATO系統(tǒng)中。Matlab仿真結(jié)果表明,該算法能減弱或消除坡道對列車運(yùn)行速度造成的影響。

    1 ATO系統(tǒng)的設(shè)計(jì)

    1.1 速度控制器的結(jié)構(gòu)

    ATO系統(tǒng)的性能指標(biāo)包括運(yùn)行的安全性、準(zhǔn)時(shí)性,以及停車精準(zhǔn)性、節(jié)能性和乘客舒適性。指標(biāo)的好壞與系統(tǒng)對速度的控制密不可分,因此,對ATO仿真系統(tǒng)的設(shè)計(jì),核心是對速度控制器的設(shè)計(jì)[2]。

    根據(jù)系統(tǒng)功能、FGPC算法理論及機(jī)車實(shí)際特點(diǎn),設(shè)計(jì)列車速度控制器(如圖1所示)。首先,速度控制器對目標(biāo)曲線進(jìn)行柔化處理,得到參考軌跡并將其作為控制器的輸入;其次,將廣義預(yù)測控制器的輸出作為列車運(yùn)動模型的輸入,并控制列車運(yùn)行速度;其三,將當(dāng)前時(shí)刻列車實(shí)際速度與目標(biāo)速度進(jìn)行比較,構(gòu)成偏差向量;最后,將偏差向量及偏差變化量進(jìn)行模糊推理并對預(yù)測誤差加以補(bǔ)償,即作用在列車運(yùn)動模型上的控制信號u(t)由廣義預(yù)測控制量u1(t)和模糊控制量u2(t)兩部分組成,使列車實(shí)現(xiàn)無誤差運(yùn)行。這種將模糊推理和GPC算法組合的方法不僅提高了系統(tǒng)的魯棒性,而且降低了對預(yù)測模型精確度的要求,對列車運(yùn)行這樣的復(fù)雜系統(tǒng)具有實(shí)際價(jià)值。

    為提高ATO系統(tǒng)的適應(yīng)性,在圖1中,將線路坡道對列車運(yùn)行速度的影響作為擾動量加入控制器中。因此,反饋的實(shí)際運(yùn)行速度包括列車在GPC算法控制下對參考軌跡的追蹤速度和環(huán)境因素對列車造成的干擾速度。模糊推理不僅對GPC誤差進(jìn)行補(bǔ)償,還對環(huán)境引起的干擾速度進(jìn)行補(bǔ)償,從而降低了環(huán)境對列車運(yùn)行速度的影響。

    圖1 速度控制器的設(shè)計(jì)

    1.2 列車運(yùn)動模型

    廣義預(yù)測控制是以受控自回歸積分滑動平均模型(CARIMA)為基礎(chǔ)的一種控制算法。設(shè)被控對象的數(shù)學(xué)模型采用如下的CARIMA方程描述:

    A(z-1)y(k)=B(z-1)Δu(k-1)+ξ(k)/δ (1)其中

    式中:

    u(k),y(k)---分別為被控對象的輸入和輸出;

    δ---差分算子,δ=1-z-1;

    ξ(k)---均值為零、方差為σ2的白噪聲。

    為了使控制效果更加良好,在實(shí)際確定目標(biāo)函數(shù)時(shí),需考慮當(dāng)前時(shí)刻的控制對系統(tǒng)未來時(shí)刻的干擾。目標(biāo)函數(shù)為

    式中:

    n---最大預(yù)測時(shí)域;

    m---控制時(shí)域;

    λ---控制加權(quán)系數(shù)。

    Δu(k+j-1)=0表示在m步后控制量不再變化。

    GPC算法的控制目標(biāo)是使被控對象輸出值y(k)盡量靠近預(yù)定值yr(k)。在實(shí)際中,控制目標(biāo)不是使y(k)直接跟隨yr(k),而是追蹤參考軌跡w(k+j)。參考軌跡為:

    式中:

    α-柔化系數(shù),0<α<1。

    由式(2)和式(3)可知,GPC算法的最終目的是通過Δu(k),Δu(k+1),…,Δu(k+m-1)的求解使得目標(biāo)函數(shù)達(dá)到最小值。

    要完成對ATO系統(tǒng)的仿真,首先要有列車運(yùn)動模型作為控制對象。在此,對文獻(xiàn)[3]中列車運(yùn)動模型進(jìn)行離散變換,并考慮坡道對列車運(yùn)行的影響,確定速度控制器的輸出差分方程為:

    式(4)中,ξ(k-1)表示k-1時(shí)刻坡道對列車速度的干擾量。將列車運(yùn)動模型的差分方程轉(zhuǎn)換為廣義預(yù)測CARIMA模型為:

    其中:

    1.3 線路模型

    在圖1中,將坡道對列車的影響作為干擾速度疊加在控制器上。因此,設(shè)計(jì)線路模型的最終目的是將坡道對列車產(chǎn)生的附加阻力轉(zhuǎn)換為速度形式。

    1.3.1 線路參數(shù)的化簡

    城市軌道交通線路由各種不同的線路條件組合而成,為減少工作量,在實(shí)際計(jì)算中需對坡道進(jìn)行化簡?;喥碌揽捎墒剑?)求出[4]。

    式中:

    i1,i2,…,in---各實(shí)際坡道坡度千分?jǐn)?shù);

    l1,l2,…,ln---各實(shí)際坡道長度;

    H1,H2---化簡后坡道的始點(diǎn)和終點(diǎn)標(biāo)高;

    ih---化簡后坡道坡度千分?jǐn)?shù);

    lh---化簡后坡道長度。

    本文以上海軌道交通16號線惠南站到惠南東站上行為例,結(jié)合式(5)對坡道進(jìn)行化簡。

    1.3.2 坡道附加阻力及其轉(zhuǎn)換

    坡道附加阻力是列車在坡道上運(yùn)行時(shí),其重力在沿下坡道方向的分力引起的。根據(jù)文獻(xiàn)[4],可得列車的坡道附加阻力為:

    式中:

    Wi---坡道附加阻力;

    i---坡道坡度千分?jǐn)?shù);

    g---重力加速度;

    M---列車質(zhì)量;

    G1---列車所受重力;

    θ---坡道夾角。

    列車運(yùn)行時(shí)的動能主要包括列車線性運(yùn)動的動能和列車旋轉(zhuǎn)部分的轉(zhuǎn)動動能。

    式中:

    Ek---列車動能;

    I---列車旋轉(zhuǎn)部分轉(zhuǎn)動慣量;

    Rh---旋轉(zhuǎn)部分換算半徑;

    w---列車旋轉(zhuǎn)部分角速度,w=v/Rh;

    γ---回轉(zhuǎn)質(zhì)量系數(shù),是列車旋轉(zhuǎn)部分的動能折算質(zhì)量與列車全部質(zhì)量的比值。

    列車在坡道上運(yùn)行,可將坡道產(chǎn)生的附加阻力視為獨(dú)自作用在勻速運(yùn)行列車上的合力。根據(jù)牛頓第二定律,結(jié)合式(6)和式(7),可得坡道附加阻力對列車運(yùn)行產(chǎn)生的干擾速度為:

    式中:

    si---當(dāng)前列車所在里程;

    si0---列車所處坡道開始里程;

    a---坡道附加阻力產(chǎn)生的加速度。

    根據(jù)式(8),可得化簡后坡道產(chǎn)生的速度干擾曲線,如圖2所示。

    圖2 坡道產(chǎn)生的干擾速度曲線

    分析圖2可知,由于單質(zhì)點(diǎn)列車模型忽略了列車的長度和車輛間的相互作用力,因此當(dāng)列車跨越換坡點(diǎn)時(shí),列車受力是瞬時(shí)變化的,從而導(dǎo)致了干擾速度的瞬時(shí)變化。

    1.4 目標(biāo)曲線的生成

    速度控制器的功能是確保列車能很好地跟隨目標(biāo)曲線運(yùn)行。當(dāng)提供給速度控制器的目標(biāo)曲線較差時(shí),無論設(shè)計(jì)的速度控制器性能多好,也不能使ATO系統(tǒng)獲得良好性能。這就要求供給速度控制器的目標(biāo)曲線必須是一條最優(yōu)的曲線。為提高系統(tǒng)節(jié)能性,本文在生成目標(biāo)曲線時(shí)采用了節(jié)能策略。即:加速階段以最大牽引力加速,中間階段盡可能地延長列車勻速運(yùn)行或惰行時(shí)間;進(jìn)站階段以最大制動力制動[5]。

    在線路模型中,由于條件限制,區(qū)間被分為120 km/h和100 km/h兩個(gè)限速不同的區(qū)段,進(jìn)、出站限速為80 km/h,線路總里程為5 899 m。根據(jù)列車運(yùn)動模型、線路模型及操縱策略得到列車運(yùn)行過程理想的v-s目標(biāo)曲線如圖3所示。

    圖3 v-s目標(biāo)曲線圖

    列車跟隨目標(biāo)曲線可能存在超調(diào)量,在設(shè)計(jì)目標(biāo)曲線時(shí),取區(qū)間最大仿真運(yùn)行速度分別為32 m/s(區(qū)間最高限速)和26.5 m/s(曲線區(qū)段限速),進(jìn)、出站最大運(yùn)行速度為20.8 m/s,稍低于ATP(列車自動保護(hù))系統(tǒng)給出的線路限速33.33 m/s、27.78 m/s和22.2 m/s。

    1.5 參數(shù)的確定

    GPC算法是基于參數(shù)模型的控制算法,模型參數(shù)的選取對速度控制器的控制性能會產(chǎn)生極大影響。對于模糊控制器來說,輸入、輸出變量的模糊論域及量化因子、比例因子等參數(shù)對系統(tǒng)的控制性能有著十分重要的影響,因此需要對這些參數(shù)進(jìn)行說明。

    1.5.1 廣義預(yù)測控制參數(shù)的選取

    根據(jù)文獻(xiàn)[6],在選擇預(yù)測時(shí)域長度n時(shí),必須使其超過受控對象的時(shí)滯和慣性部分。n的取值越大,系統(tǒng)的抗干擾和抗模型失配能力越好,但n取值太大會增加計(jì)算量,故取n=9??刂茣r(shí)域長度m取值越大則跟蹤性能越好,但取值太大會導(dǎo)致系統(tǒng)不穩(wěn)定,此處取m=2??刂萍訖?quán)系數(shù)λ用來限制控制增量Δu的劇烈變化,取λ=2。柔化系數(shù)α取值越大,系統(tǒng)的魯棒性越好,但取值太大會導(dǎo)致系統(tǒng)的快速性降低,故取α=0.92。采樣周期T過大會使受控模型不準(zhǔn),控制效果下降,T太小則會造成控制器延時(shí),影響系統(tǒng)的穩(wěn)定性,本文結(jié)合ATO系統(tǒng)的特點(diǎn),取T=5。需特別說明的是,在選擇采樣周期時(shí),為改善啟動階段的估計(jì)和控制,可首先測出對象的單位階躍響應(yīng)系數(shù),以此作為初始參數(shù)。

    1.5.2 模糊控制參數(shù)的選取

    1.5.2.1 模糊語言變量的語言值分檔

    本文采用二維模糊控制結(jié)構(gòu),選取預(yù)測輸出與參考軌跡的偏差e、偏差變化量ec作為模糊輸入量,模糊控制量u2作為輸出量。

    模糊輸入量的物理論域上下限分別為預(yù)測輸出與參考軌跡的誤差、誤差變化量中的最大值和最小值;模糊控制量的物理論域取為u2=[-10,10]。兩個(gè)模糊輸入變量、輸出控制量的模糊論域分別為[-6,6]和[-9.333,9.333]。模糊控制器的輸入變量e、ec及輸出變量u2的模糊子集均為{NB,NM,NS,ZO,PS,PM,PB},含義依次為負(fù)大、負(fù)中、負(fù)小、零、正小、正中、正大。

    1.5.2.2 模糊子集隸屬函數(shù)的確定

    根據(jù)系統(tǒng)特性,選用均勻分布的三角形函數(shù)作為模糊輸入隸屬函數(shù),非均勻分布的三角形函數(shù)作為模糊輸出隸屬函數(shù),且隸屬度函數(shù)的中心值分別為[-6,-4,-2,0,2,4,6]和[-6,-2.667,-0.667,0.667,2.667,6]。圖4、圖5分別為輸入變量模糊子集隸屬函數(shù)和輸出變量模糊子集隸屬函數(shù)。

    圖4 輸入變量e、ec模糊子集隸屬函數(shù)

    圖5 輸出變量u2模糊子集隸屬函數(shù)

    1.5.2.3 模糊控制規(guī)則的確定

    為了使系統(tǒng)更好地利用模糊控制器的反饋補(bǔ)償作用,必須認(rèn)真分析模糊控制器的輸入變量,并結(jié)合列車操作者的工作經(jīng)驗(yàn)建立模糊規(guī)則庫,本文以表1為模糊規(guī)則。

    表1 模糊控制規(guī)則表

    2 仿真結(jié)果

    2.1 未考慮坡道時(shí)的仿真結(jié)果

    圖6為不考慮坡道的情況,且列車分別在GPC隱式算法和FGPC算法的控制下跟蹤v-t目標(biāo)曲線運(yùn)行時(shí)的仿真結(jié)果。GPC算法參數(shù)已在前文做介紹。模糊控制的量化因子為ke=0.714,kec= 5.12,比例因子為ku=2.14。

    圖6 追蹤v-t目標(biāo)曲線仿真圖

    分析圖6可知,前50 s,列車在GPC隱式算法的控制下運(yùn)行時(shí),列車速度變化緩慢,不能快速地跟隨目標(biāo)曲線,跟蹤誤差較大;經(jīng)過模糊控制量的補(bǔ)償作用后,速度控制器不僅能快速跟蹤v-t目標(biāo)曲線,減小跟蹤誤差,還能確保列車準(zhǔn)時(shí)到站;同時(shí),當(dāng)列車在區(qū)間低限速區(qū)運(yùn)行時(shí),基于FGPC算法的速度控制器能更好地保證列車運(yùn)行的安全性。由此可知,基于FGPC算法的速度控制器能很好地優(yōu)化ATO系統(tǒng)的性能。

    圖7為不考慮坡道的情況,且列車分別在GPC隱式算法和FGPC算法的控制下跟蹤a-t目標(biāo)曲線運(yùn)行時(shí)的仿真結(jié)果。其中,GPC算法參數(shù)與跟蹤vt目標(biāo)曲線時(shí)相同,模糊控制的量化因子為ke= 7.13,kec=6.67,比例因子為ku=2.14。

    圖7 追蹤a-t目標(biāo)曲線仿真圖

    由圖7可知,在50~100 s、180~200 s的時(shí)間段內(nèi),列車在基于GPC隱式算法的速度控制器下運(yùn)行時(shí)跟蹤誤差基本為0;但在第25 s和220 s時(shí),列車經(jīng)過模糊控制量的補(bǔ)償后能快速跟蹤目標(biāo)曲線,減小了跟蹤誤差,增強(qiáng)了系統(tǒng)的動態(tài)性和實(shí)時(shí)性。因此,基于FGPC算法的ATO系統(tǒng)具有更好的性能。

    2.2 考慮坡道的仿真結(jié)果

    以圖2中化簡后的第二個(gè)坡道為例進(jìn)行仿真,驗(yàn)證本文設(shè)計(jì)的速度控制器在有坡道的線路上也能滿足ATO系統(tǒng)的各項(xiàng)性能指標(biāo)。化簡后的第二個(gè)坡道為下坡道,坡長為1 350 m,坡度為4‰,坡道的起始位移為240.3 m。圖8為列車分別在有坡道與無坡道情況下運(yùn)行時(shí)的仿真曲線圖。由圖8可知,當(dāng)列車運(yùn)行至240.3 m處時(shí),由于下坡道的存在,列車運(yùn)行速度迅速增加,不僅超出了目標(biāo)曲線,而且超越了限速曲線。因此,需對坡道產(chǎn)生的速度影響進(jìn)行補(bǔ)償控制。

    圖8 存在坡道干擾時(shí)v-s曲線圖

    圖9為經(jīng)過模糊補(bǔ)償之后的仿真曲線圖。在模糊控制量的補(bǔ)償作用下,仿真曲線滿足了低于限速曲線的要求,保證了列車的安全運(yùn)行,使系統(tǒng)的性能得到了很大提高。由圖9可以看出,應(yīng)用FGPC算法設(shè)計(jì)的速度控制器能適應(yīng)任何線路情況,能保證列車的安全運(yùn)行且具有可靠的跟隨性。

    圖9 有坡道時(shí)FGPC算法仿真曲線圖

    圖10為速度控制器在FGPC算法控制下,列車分別運(yùn)行在無坡道和有坡道線路上的仿真曲線圖。當(dāng)列車在有坡道的線路上運(yùn)行時(shí),列車運(yùn)行速度變化較慢,但整體的跟隨效果與在無坡道的線路上運(yùn)行時(shí)基本相同。雖然其犧牲了ATO系統(tǒng)部分的快速性,降低了系統(tǒng)的響應(yīng)速度,但確保了列車運(yùn)行的安全性、平穩(wěn)性和舒適性,體現(xiàn)了基于FGPC算法的速度控制器的良好性能。

    圖10 有坡道、無坡道時(shí)速度跟隨曲線

    3 結(jié)語

    本文根據(jù)ATO系統(tǒng)的性能設(shè)計(jì)了基于FGPC算法的ATO系統(tǒng)速度控制器,并根據(jù)列車運(yùn)動模型及線路模型生成優(yōu)化的目標(biāo)曲線。通過系統(tǒng)仿真,并與GPC隱式算法的控制效果進(jìn)行比較,證明了將FGPC算法應(yīng)用于ATO系統(tǒng)是可行的,且該算法不僅能改善ATO系統(tǒng)的各項(xiàng)性能指標(biāo),還能改善坡道對列車運(yùn)行速度的影響,具有很好的移植性和適應(yīng)性。

    [1] 宮赤坤,華澤釗.基于模糊補(bǔ)償?shù)膹V義預(yù)測控制[J].吉林大學(xué)學(xué)報(bào),2002,32(2):78.

    [2] 黃秀玲,王長林.廣義預(yù)測控制算法在ATO速度控制器中的研究[J].軌道交通縱橫,2009,1(5):13.

    [3] 劉智勇.高速磁浮列車安全速度曲線算法及緊急制動控制研究[D].杭州:浙江大學(xué),2008.

    [4] 石紅國.列車運(yùn)行過程仿真及優(yōu)化研究[D].成都:西南交通大學(xué),2006.

    [5] 崔世文,馮曉云.列車優(yōu)化操縱與自動駕駛模式的研究與仿真[J].鐵道機(jī)車車輛,2005,1(10):1.

    [6] 李國勇.智能控制及其MATLAB實(shí)現(xiàn)[M].北京:電子工業(yè)出版社,2005:289.

    [7] 宗明,張國侯.軌道交通列車自動運(yùn)行仿真平臺的開發(fā)研究[J].城市軌道交通研究,2012(6):52.

    Simulation of ATO System Based on Fuzzy Generalized Predictive Control Algorithm

    Gu Guimei,Pu Songzhi

    According to the characteristics of automatic train operation(ATO)system,a speed controller of ATO system based on fuzzy generalized predictive control algorithms(FGPCA)is designed,a line model combined with the actual lines is established and an optimized target curveis generated.Then the speed controller based on FGPCA and implicit generalized predictive control algorithms with MATLAB issimulated.The result proves that the speed controller based on FGPCA could not only optimize the properties of ATO system,but could also weaken the impact of ramps on the train speed to some extent.Therefore,this newly designed speed controller has a good transportability.

    automatic train operation(ATO);fuzzy generalized predictive control;simulation

    U 231.6

    2012-06-18)

    猜你喜歡
    坡道列車運(yùn)行模糊控制
    AMT坡道起步輔助控制策略
    平朔東露天礦工作幫移動坡道優(yōu)化研究
    改善地鐵列車運(yùn)行舒適度方案探討
    T-S模糊控制綜述與展望
    列車運(yùn)行控制系統(tǒng)技術(shù)發(fā)展趨勢分析
    基于模糊控制的PLC在溫度控制中的應(yīng)用
    電子制作(2016年21期)2016-05-17 03:53:01
    基于模糊控制的恒壓供水系統(tǒng)的研究
    用于光伏MPPT中的模糊控制占空比擾動法
    相同徑路的高速列車運(yùn)行圖編制方法
    基于傾角傳感器的坡道角度識別研究
    50天的宝宝边吃奶边哭怎么回事| 男女那种视频在线观看| 免费看日本二区| tocl精华| 又大又爽又粗| 久久天躁狠狠躁夜夜2o2o| 成在线人永久免费视频| 黄片大片在线免费观看| 精品欧美一区二区三区在线| 中文字幕高清在线视频| 亚洲va日本ⅴa欧美va伊人久久| 天天一区二区日本电影三级| 999久久久国产精品视频| 国产三级在线视频| АⅤ资源中文在线天堂| 又大又爽又粗| 日本在线视频免费播放| 国产精品乱码一区二三区的特点| 嫩草影视91久久| 国产精品av视频在线免费观看| 亚洲无线在线观看| 久久久久九九精品影院| www.精华液| 精品久久久久久久人妻蜜臀av| 老司机福利观看| 一级毛片精品| 精品一区二区三区四区五区乱码| 久久精品国产综合久久久| 国产91精品成人一区二区三区| 国产精品电影一区二区三区| 国产午夜精品论理片| 亚洲激情在线av| 久久精品影院6| 国产成人啪精品午夜网站| 18禁黄网站禁片免费观看直播| 美女 人体艺术 gogo| 黄色毛片三级朝国网站| 熟女电影av网| 国模一区二区三区四区视频 | 我的老师免费观看完整版| 香蕉国产在线看| 99精品久久久久人妻精品| АⅤ资源中文在线天堂| 欧美一级毛片孕妇| 男人舔女人的私密视频| 免费看美女性在线毛片视频| 麻豆成人av在线观看| 亚洲欧美日韩高清专用| 午夜福利18| 夜夜爽天天搞| 熟女电影av网| 亚洲免费av在线视频| 日日夜夜操网爽| 精品不卡国产一区二区三区| 国产蜜桃级精品一区二区三区| 在线观看美女被高潮喷水网站 | 亚洲自偷自拍图片 自拍| 欧美+亚洲+日韩+国产| 亚洲精品美女久久av网站| 成人手机av| 观看免费一级毛片| 老汉色∧v一级毛片| 99久久精品热视频| 麻豆av在线久日| 无人区码免费观看不卡| av超薄肉色丝袜交足视频| 中文字幕精品亚洲无线码一区| 一二三四社区在线视频社区8| www.自偷自拍.com| 婷婷精品国产亚洲av| 国产精品一区二区精品视频观看| 50天的宝宝边吃奶边哭怎么回事| 脱女人内裤的视频| 人妻久久中文字幕网| 夜夜爽天天搞| 丝袜美腿诱惑在线| 欧美成人免费av一区二区三区| 免费在线观看视频国产中文字幕亚洲| 国产三级在线视频| 久久亚洲真实| 亚洲精品国产一区二区精华液| 在线十欧美十亚洲十日本专区| 亚洲国产精品sss在线观看| 俄罗斯特黄特色一大片| 国产精品亚洲av一区麻豆| 91在线观看av| 97人妻精品一区二区三区麻豆| 久久中文字幕一级| 欧美极品一区二区三区四区| 亚洲人成网站在线播放欧美日韩| 五月玫瑰六月丁香| 久久久久性生活片| 亚洲男人的天堂狠狠| 成人国语在线视频| www国产在线视频色| 欧美成狂野欧美在线观看| 久久中文字幕人妻熟女| 国产91精品成人一区二区三区| 成人三级黄色视频| 老熟妇仑乱视频hdxx| 精品久久久久久成人av| 香蕉av资源在线| 精品久久久久久,| 少妇粗大呻吟视频| 嫩草影视91久久| 亚洲真实伦在线观看| 人人妻,人人澡人人爽秒播| 国产精品亚洲美女久久久| 怎么达到女性高潮| 成熟少妇高潮喷水视频| 麻豆av在线久日| 亚洲18禁久久av| 亚洲国产精品999在线| 亚洲 国产 在线| 一进一出抽搐gif免费好疼| 一夜夜www| 国产三级在线视频| 精品一区二区三区av网在线观看| 免费搜索国产男女视频| 日韩三级视频一区二区三区| 精品电影一区二区在线| 欧美日韩瑟瑟在线播放| 看片在线看免费视频| bbb黄色大片| 日本a在线网址| 日韩欧美一区二区三区在线观看| 国产精品香港三级国产av潘金莲| 国内少妇人妻偷人精品xxx网站 | 亚洲精品在线观看二区| 99精品在免费线老司机午夜| 精品日产1卡2卡| 麻豆成人午夜福利视频| 久久人妻福利社区极品人妻图片| 变态另类丝袜制服| 波多野结衣高清无吗| 亚洲欧美日韩东京热| 在线播放国产精品三级| 午夜精品久久久久久毛片777| 淫秽高清视频在线观看| 国产精品久久久av美女十八| 欧美精品亚洲一区二区| 制服人妻中文乱码| 国产视频内射| 性色av乱码一区二区三区2| 美女大奶头视频| 一区二区三区激情视频| 黑人操中国人逼视频| 成人亚洲精品av一区二区| 久久久久国内视频| av有码第一页| 美女午夜性视频免费| 老汉色∧v一级毛片| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产精品久久男人天堂| 一级作爱视频免费观看| 91av网站免费观看| 18禁国产床啪视频网站| 亚洲av成人不卡在线观看播放网| 欧美激情久久久久久爽电影| 国产激情欧美一区二区| 欧美另类亚洲清纯唯美| 免费在线观看黄色视频的| 亚洲一码二码三码区别大吗| 国产av不卡久久| 在线观看一区二区三区| 一区二区三区高清视频在线| 99国产精品一区二区三区| 久久久久精品国产欧美久久久| 亚洲色图 男人天堂 中文字幕| 香蕉丝袜av| 日本精品一区二区三区蜜桃| 亚洲成人精品中文字幕电影| 999精品在线视频| 白带黄色成豆腐渣| 国产成人精品久久二区二区免费| 正在播放国产对白刺激| 国产视频内射| 99在线视频只有这里精品首页| 搡老妇女老女人老熟妇| 国产一区二区在线av高清观看| 国产一区在线观看成人免费| 日韩欧美国产一区二区入口| 日韩欧美国产一区二区入口| 日本免费一区二区三区高清不卡| 老熟妇乱子伦视频在线观看| 成人18禁在线播放| 亚洲第一欧美日韩一区二区三区| 美女 人体艺术 gogo| 精品一区二区三区四区五区乱码| 亚洲av电影不卡..在线观看| 老熟妇仑乱视频hdxx| 国产亚洲精品久久久久5区| 亚洲第一电影网av| 久久欧美精品欧美久久欧美| 在线观看午夜福利视频| 国产精品亚洲美女久久久| 在线永久观看黄色视频| 国内精品久久久久久久电影| 亚洲人成网站在线播放欧美日韩| 少妇裸体淫交视频免费看高清 | 夜夜看夜夜爽夜夜摸| videosex国产| 亚洲天堂国产精品一区在线| 少妇人妻一区二区三区视频| 欧美人与性动交α欧美精品济南到| 国产伦一二天堂av在线观看| 欧美黑人欧美精品刺激| 婷婷丁香在线五月| 18禁观看日本| 丝袜人妻中文字幕| 国模一区二区三区四区视频 | 999久久久国产精品视频| 91av网站免费观看| 日韩大尺度精品在线看网址| 日韩 欧美 亚洲 中文字幕| 国产av又大| 超碰成人久久| 91成年电影在线观看| 久久伊人香网站| 999久久久国产精品视频| 亚洲黑人精品在线| 天堂av国产一区二区熟女人妻 | 1024手机看黄色片| 亚洲欧美日韩高清在线视频| 精品午夜福利视频在线观看一区| 级片在线观看| 麻豆成人av在线观看| 99热只有精品国产| 日本 欧美在线| www.精华液| 国产乱人伦免费视频| 91麻豆精品激情在线观看国产| 又大又爽又粗| 中文在线观看免费www的网站 | 久久中文字幕一级| 久久精品国产99精品国产亚洲性色| 欧美 亚洲 国产 日韩一| 又粗又爽又猛毛片免费看| 国产激情欧美一区二区| 日本免费一区二区三区高清不卡| 久久久水蜜桃国产精品网| 黄色视频,在线免费观看| 亚洲自拍偷在线| 亚洲狠狠婷婷综合久久图片| 黄频高清免费视频| 性欧美人与动物交配| 欧美性长视频在线观看| 男人舔奶头视频| 在线观看美女被高潮喷水网站 | 亚洲精华国产精华精| 亚洲片人在线观看| 露出奶头的视频| 日韩有码中文字幕| 欧美不卡视频在线免费观看 | 50天的宝宝边吃奶边哭怎么回事| 国产欧美日韩一区二区三| 免费在线观看成人毛片| 麻豆久久精品国产亚洲av| 亚洲人成77777在线视频| 国产成人一区二区三区免费视频网站| 三级男女做爰猛烈吃奶摸视频| 19禁男女啪啪无遮挡网站| 国产成人啪精品午夜网站| 12—13女人毛片做爰片一| 久久这里只有精品中国| 黄色视频不卡| 国产精品免费视频内射| 白带黄色成豆腐渣| 亚洲自拍偷在线| 国产伦人伦偷精品视频| 国产成人欧美在线观看| 午夜精品一区二区三区免费看| 亚洲国产精品合色在线| 18禁黄网站禁片免费观看直播| 99国产综合亚洲精品| 中文字幕精品亚洲无线码一区| 人妻丰满熟妇av一区二区三区| 成人高潮视频无遮挡免费网站| 婷婷精品国产亚洲av在线| 99精品欧美一区二区三区四区| 亚洲人成77777在线视频| 精品久久久久久久毛片微露脸| 亚洲五月婷婷丁香| 午夜久久久久精精品| 老司机深夜福利视频在线观看| 久久久久久九九精品二区国产 | 久久性视频一级片| 国产成人aa在线观看| 日日摸夜夜添夜夜添小说| 老司机靠b影院| 亚洲精品色激情综合| 亚洲国产欧美一区二区综合| 亚洲一卡2卡3卡4卡5卡精品中文| 国产av一区在线观看免费| 99国产精品一区二区三区| 日韩欧美 国产精品| 长腿黑丝高跟| 亚洲欧美日韩高清专用| 日日夜夜操网爽| 人人妻人人澡欧美一区二区| 免费观看人在逋| 日本三级黄在线观看| 亚洲中文字幕一区二区三区有码在线看 | 精品久久久久久久久久久久久| 久久精品91蜜桃| 国产午夜精品久久久久久| 成年免费大片在线观看| 国内精品一区二区在线观看| 久久亚洲精品不卡| 黄片大片在线免费观看| 天堂av国产一区二区熟女人妻 | 岛国在线免费视频观看| 国产又黄又爽又无遮挡在线| 18禁观看日本| 久久人妻福利社区极品人妻图片| www.精华液| 悠悠久久av| 在线观看美女被高潮喷水网站 | 国产久久久一区二区三区| 在线观看www视频免费| 亚洲人成77777在线视频| 熟女少妇亚洲综合色aaa.| 欧美三级亚洲精品| 天天添夜夜摸| www.精华液| 老司机靠b影院| 黄频高清免费视频| 制服人妻中文乱码| 国产精品久久久久久亚洲av鲁大| 男女之事视频高清在线观看| 天堂影院成人在线观看| 黄色 视频免费看| 夜夜爽天天搞| 国产aⅴ精品一区二区三区波| 国产一级毛片七仙女欲春2| or卡值多少钱| 欧洲精品卡2卡3卡4卡5卡区| 老司机靠b影院| 少妇裸体淫交视频免费看高清 | 国产精品亚洲一级av第二区| 视频区欧美日本亚洲| 91麻豆av在线| 日韩中文字幕欧美一区二区| 精品久久久久久久毛片微露脸| 精品久久久久久成人av| 两个人的视频大全免费| 中文字幕熟女人妻在线| 欧美日韩亚洲国产一区二区在线观看| 国产一级毛片七仙女欲春2| 久久人妻福利社区极品人妻图片| 亚洲欧洲精品一区二区精品久久久| 国产精品综合久久久久久久免费| 日日干狠狠操夜夜爽| 精品第一国产精品| 2021天堂中文幕一二区在线观| 亚洲国产欧美一区二区综合| 国产1区2区3区精品| 久久亚洲精品不卡| 国产成人aa在线观看| 黄色成人免费大全| 亚洲熟妇熟女久久| 久久精品国产99精品国产亚洲性色| 搞女人的毛片| 国产av一区二区精品久久| 国产成人aa在线观看| 白带黄色成豆腐渣| 国产亚洲欧美在线一区二区| 老司机福利观看| 老熟妇仑乱视频hdxx| 亚洲国产欧美一区二区综合| 久久热在线av| 伊人久久大香线蕉亚洲五| 久久天躁狠狠躁夜夜2o2o| 久久国产精品人妻蜜桃| 亚洲色图av天堂| 久久这里只有精品中国| 在线观看日韩欧美| 欧美绝顶高潮抽搐喷水| 一二三四社区在线视频社区8| 免费在线观看成人毛片| 又紧又爽又黄一区二区| 午夜福利免费观看在线| 亚洲国产精品久久男人天堂| 亚洲av成人不卡在线观看播放网| 久久精品成人免费网站| netflix在线观看网站| 久久中文看片网| 女人被狂操c到高潮| 悠悠久久av| 国产欧美日韩一区二区精品| 一进一出抽搐gif免费好疼| 操出白浆在线播放| 女人高潮潮喷娇喘18禁视频| 欧美一级毛片孕妇| 国产精品久久久人人做人人爽| 午夜成年电影在线免费观看| 国产精品,欧美在线| 最近最新中文字幕大全电影3| 国产91精品成人一区二区三区| 精品第一国产精品| 国产精品九九99| 精品高清国产在线一区| 午夜亚洲福利在线播放| 国产精品99久久99久久久不卡| 1024视频免费在线观看| 亚洲中文字幕一区二区三区有码在线看 | 久久久久国内视频| 日本黄色视频三级网站网址| 欧美日韩福利视频一区二区| 女人高潮潮喷娇喘18禁视频| 午夜福利成人在线免费观看| 少妇被粗大的猛进出69影院| 黄色a级毛片大全视频| 十八禁网站免费在线| 免费一级毛片在线播放高清视频| a级毛片a级免费在线| 亚洲国产欧美网| 日韩欧美免费精品| 国产精品国产高清国产av| av超薄肉色丝袜交足视频| 一边摸一边做爽爽视频免费| 别揉我奶头~嗯~啊~动态视频| 嫩草影视91久久| 国产亚洲欧美98| 欧美在线黄色| 成人18禁高潮啪啪吃奶动态图| 88av欧美| 欧美高清成人免费视频www| 少妇裸体淫交视频免费看高清 | 亚洲av电影在线进入| 中文资源天堂在线| 最近最新中文字幕大全免费视频| 国产视频一区二区在线看| 精品欧美国产一区二区三| 大型黄色视频在线免费观看| 狂野欧美激情性xxxx| 国产不卡一卡二| 欧美国产日韩亚洲一区| 男女床上黄色一级片免费看| 欧美高清成人免费视频www| 黄频高清免费视频| 小说图片视频综合网站| av福利片在线| 亚洲欧美日韩无卡精品| 波多野结衣巨乳人妻| 女生性感内裤真人,穿戴方法视频| 人人妻人人看人人澡| 国产黄色小视频在线观看| 国产乱人伦免费视频| 禁无遮挡网站| 91字幕亚洲| 夜夜夜夜夜久久久久| avwww免费| 黑人欧美特级aaaaaa片| 久久久久久大精品| 亚洲欧美日韩东京热| 亚洲精品久久国产高清桃花| 免费在线观看影片大全网站| 国产区一区二久久| 国产免费av片在线观看野外av| 国产精品一区二区免费欧美| 午夜福利免费观看在线| 国产精品一区二区三区四区免费观看 | 久久这里只有精品中国| 人妻丰满熟妇av一区二区三区| or卡值多少钱| 国产乱人伦免费视频| 两性夫妻黄色片| 国产高清视频在线播放一区| 欧美3d第一页| 1024香蕉在线观看| 色综合婷婷激情| 国产成人精品久久二区二区91| 全区人妻精品视频| bbb黄色大片| 精华霜和精华液先用哪个| 午夜久久久久精精品| 亚洲精品国产精品久久久不卡| 91在线观看av| 久久精品91无色码中文字幕| 久久久久久九九精品二区国产 | 国产精品,欧美在线| 国产激情偷乱视频一区二区| 搡老妇女老女人老熟妇| 国产视频一区二区在线看| 亚洲成人久久爱视频| 97人妻精品一区二区三区麻豆| 级片在线观看| 露出奶头的视频| 亚洲午夜精品一区,二区,三区| 又大又爽又粗| 亚洲 国产 在线| 成人三级做爰电影| 天堂动漫精品| 亚洲一区二区三区不卡视频| 国产精品香港三级国产av潘金莲| 一级毛片女人18水好多| 精品国产乱码久久久久久男人| 黄色成人免费大全| 国产av又大| 又粗又爽又猛毛片免费看| 亚洲欧美激情综合另类| av中文乱码字幕在线| 国产成人av教育| 97碰自拍视频| 亚洲人成伊人成综合网2020| 久久久久久人人人人人| 欧美日韩一级在线毛片| 神马国产精品三级电影在线观看 | 岛国视频午夜一区免费看| 中文字幕高清在线视频| 一级片免费观看大全| 在线观看www视频免费| 成年免费大片在线观看| 亚洲,欧美精品.| 91九色精品人成在线观看| 欧美黑人巨大hd| 12—13女人毛片做爰片一| 91麻豆av在线| 亚洲精品美女久久av网站| 亚洲欧美日韩东京热| 国产av又大| 国产在线观看jvid| 欧美久久黑人一区二区| 极品教师在线免费播放| 日韩中文字幕欧美一区二区| 国产69精品久久久久777片 | 亚洲欧美日韩高清专用| 五月伊人婷婷丁香| 国产精品电影一区二区三区| 免费人成视频x8x8入口观看| 91字幕亚洲| svipshipincom国产片| 女生性感内裤真人,穿戴方法视频| 欧美激情久久久久久爽电影| 久99久视频精品免费| av欧美777| 窝窝影院91人妻| 色播亚洲综合网| 老熟妇乱子伦视频在线观看| 日韩av在线大香蕉| 国产精品免费视频内射| 一边摸一边抽搐一进一小说| 麻豆成人av在线观看| 麻豆av在线久日| 久热爱精品视频在线9| 国模一区二区三区四区视频 | 国产三级中文精品| 欧美3d第一页| 久久精品aⅴ一区二区三区四区| 天堂影院成人在线观看| 国产精品一区二区三区四区久久| 午夜精品一区二区三区免费看| 床上黄色一级片| 亚洲一区二区三区不卡视频| 成年人黄色毛片网站| 少妇人妻一区二区三区视频| 中文字幕久久专区| 免费搜索国产男女视频| 久久久久亚洲av毛片大全| 观看免费一级毛片| 亚洲av第一区精品v没综合| 日本黄大片高清| 校园春色视频在线观看| 亚洲中文av在线| 一级黄色大片毛片| 久久久久久亚洲精品国产蜜桃av| 亚洲18禁久久av| 岛国在线观看网站| tocl精华| 久久精品夜夜夜夜夜久久蜜豆 | 免费在线观看完整版高清| 香蕉久久夜色| 久久久国产成人免费| 日本 av在线| 男男h啪啪无遮挡| 俺也久久电影网| 精华霜和精华液先用哪个| 黄色毛片三级朝国网站| 欧美成人一区二区免费高清观看 | 一区二区三区国产精品乱码| 亚洲美女黄片视频| 一卡2卡三卡四卡精品乱码亚洲| 精华霜和精华液先用哪个| 色av中文字幕| 日韩欧美一区二区三区在线观看| 99热这里只有是精品50| 淫秽高清视频在线观看| 日韩国内少妇激情av| 欧美成人一区二区免费高清观看 | 亚洲五月婷婷丁香| 777久久人妻少妇嫩草av网站| 久久中文字幕一级| 亚洲av中文字字幕乱码综合| 亚洲国产欧洲综合997久久,| 午夜成年电影在线免费观看| 国产一区二区在线观看日韩 | 国产高清videossex| 久久精品aⅴ一区二区三区四区| 99热这里只有是精品50| 精品久久久久久久毛片微露脸| 亚洲午夜精品一区,二区,三区| 日本一二三区视频观看| 欧美又色又爽又黄视频| 精品午夜福利视频在线观看一区| 在线免费观看的www视频| 色综合婷婷激情| 99久久无色码亚洲精品果冻| 国产99白浆流出| 亚洲欧美激情综合另类| 国产主播在线观看一区二区| 看免费av毛片| 日本免费a在线| 精品不卡国产一区二区三区| 亚洲人与动物交配视频|