• <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)行圖編制方法
    基于傾角傳感器的坡道角度識別研究
    两人在一起打扑克的视频| www.自偷自拍.com| 色噜噜av男人的天堂激情| 亚洲av美国av| 国产精品久久久久久人妻精品电影| 亚洲精品久久国产高清桃花| 久久久成人免费电影| 久久精品国产99精品国产亚洲性色| av视频在线观看入口| 亚洲熟妇中文字幕五十中出| 美女黄网站色视频| 久久久久久久久久黄片| 国产一区在线观看成人免费| 在线观看一区二区三区| 在线视频色国产色| 午夜福利欧美成人| 亚洲精品在线观看二区| 男人舔奶头视频| 国内精品美女久久久久久| 成人av一区二区三区在线看| 国产三级中文精品| 国产91精品成人一区二区三区| 一个人免费在线观看电影 | 淫妇啪啪啪对白视频| 欧美激情久久久久久爽电影| 久久欧美精品欧美久久欧美| av在线天堂中文字幕| 真人一进一出gif抽搐免费| 法律面前人人平等表现在哪些方面| 亚洲av片天天在线观看| av中文乱码字幕在线| 国产一区二区在线av高清观看| 男人的好看免费观看在线视频| 欧洲精品卡2卡3卡4卡5卡区| 三级男女做爰猛烈吃奶摸视频| 巨乳人妻的诱惑在线观看| 国产精品精品国产色婷婷| 精品一区二区三区视频在线 | 国产精品电影一区二区三区| 日本黄色视频三级网站网址| 男人舔女人下体高潮全视频| 亚洲18禁久久av| 亚洲精华国产精华精| 欧美黄色淫秽网站| 国产伦在线观看视频一区| 日韩欧美国产一区二区入口| 免费看美女性在线毛片视频| 欧美一级毛片孕妇| 网址你懂的国产日韩在线| 天天躁狠狠躁夜夜躁狠狠躁| 嫩草影院精品99| 国产美女午夜福利| 在线国产一区二区在线| 老熟妇仑乱视频hdxx| 深夜精品福利| 亚洲成人久久爱视频| 天堂av国产一区二区熟女人妻| 日本黄大片高清| 在线免费观看的www视频| 麻豆av在线久日| 国产精品 欧美亚洲| 岛国在线观看网站| 免费av毛片视频| 最近最新免费中文字幕在线| 国产成人啪精品午夜网站| 757午夜福利合集在线观看| 五月伊人婷婷丁香| 亚洲国产欧美一区二区综合| xxx96com| 12—13女人毛片做爰片一| 国产精品一区二区三区四区久久| 国产真实乱freesex| 精品国产乱子伦一区二区三区| 999精品在线视频| 天堂动漫精品| 久久这里只有精品19| 麻豆一二三区av精品| 国产精品综合久久久久久久免费| 国产精品日韩av在线免费观看| 精品熟女少妇八av免费久了| 好看av亚洲va欧美ⅴa在| av天堂在线播放| 国产精品日韩av在线免费观看| 亚洲人与动物交配视频| 国产伦一二天堂av在线观看| 国产真人三级小视频在线观看| cao死你这个sao货| 国产又黄又爽又无遮挡在线| 日韩人妻高清精品专区| 偷拍熟女少妇极品色| 黄色片一级片一级黄色片| 亚洲av中文字字幕乱码综合| 在线视频色国产色| 桃红色精品国产亚洲av| 亚洲 欧美一区二区三区| 免费在线观看影片大全网站| 国产黄a三级三级三级人| 国产真实乱freesex| 精品久久久久久久久久免费视频| 国产高清视频在线播放一区| a级毛片a级免费在线| 一进一出抽搐动态| 欧美日韩瑟瑟在线播放| 午夜免费成人在线视频| 国产精品自产拍在线观看55亚洲| 精品久久久久久成人av| 国产精品1区2区在线观看.| 国产免费av片在线观看野外av| 母亲3免费完整高清在线观看| 国产成人欧美在线观看| 免费观看精品视频网站| 欧美日韩国产亚洲二区| 久久精品综合一区二区三区| 国内毛片毛片毛片毛片毛片| 久久亚洲真实| 老鸭窝网址在线观看| 国产精品久久视频播放| 国产欧美日韩一区二区精品| 国产单亲对白刺激| 脱女人内裤的视频| 天堂av国产一区二区熟女人妻| 特级一级黄色大片| 中文亚洲av片在线观看爽| 波多野结衣高清无吗| 久久亚洲精品不卡| 亚洲av成人av| 两个人视频免费观看高清| 麻豆国产av国片精品| 精品99又大又爽又粗少妇毛片 | 亚洲黑人精品在线| av女优亚洲男人天堂 | 国产乱人视频| 国产精品免费一区二区三区在线| 欧美性猛交黑人性爽| 亚洲精品久久国产高清桃花| 亚洲国产精品久久男人天堂| 精品国产乱子伦一区二区三区| 亚洲国产高清在线一区二区三| 免费av毛片视频| 岛国视频午夜一区免费看| 夜夜爽天天搞| 亚洲国产精品成人综合色| 欧美日韩一级在线毛片| av黄色大香蕉| 国产精品1区2区在线观看.| 久久精品夜夜夜夜夜久久蜜豆| 少妇裸体淫交视频免费看高清| 成人性生交大片免费视频hd| 国产精品久久电影中文字幕| 国产 一区 欧美 日韩| 搡老熟女国产l中国老女人| 亚洲熟女毛片儿| 黄色 视频免费看| 亚洲欧美日韩无卡精品| 五月玫瑰六月丁香| 观看美女的网站| 99久久综合精品五月天人人| 一本综合久久免费| 两个人的视频大全免费| 免费电影在线观看免费观看| 国内毛片毛片毛片毛片毛片| 国产淫片久久久久久久久 | 精品不卡国产一区二区三区| 97碰自拍视频| 亚洲国产欧洲综合997久久,| 18禁美女被吸乳视频| 国产精华一区二区三区| 老司机福利观看| 啦啦啦观看免费观看视频高清| 熟妇人妻久久中文字幕3abv| 久久久久久九九精品二区国产| 一进一出抽搐gif免费好疼| 美女黄网站色视频| 日本免费a在线| 丝袜人妻中文字幕| 亚洲人成电影免费在线| 国产综合懂色| 日韩免费av在线播放| 日本a在线网址| 亚洲欧洲精品一区二区精品久久久| 午夜福利在线观看免费完整高清在 | 亚洲欧美日韩高清专用| 在线十欧美十亚洲十日本专区| 草草在线视频免费看| 真人一进一出gif抽搐免费| 国产亚洲精品一区二区www| 一区福利在线观看| 给我免费播放毛片高清在线观看| 又黄又粗又硬又大视频| 十八禁网站免费在线| 国内精品久久久久久久电影| 精品福利观看| 90打野战视频偷拍视频| 国产v大片淫在线免费观看| 曰老女人黄片| 国产成年人精品一区二区| 美女 人体艺术 gogo| 特级一级黄色大片| 成年女人永久免费观看视频| 久久伊人香网站| 久久久久久人人人人人| 亚洲精华国产精华精| 国产av麻豆久久久久久久| 日韩 欧美 亚洲 中文字幕| 国产亚洲av嫩草精品影院| 亚洲国产日韩欧美精品在线观看 | 日日干狠狠操夜夜爽| 亚洲在线自拍视频| 中文字幕人成人乱码亚洲影| 国产日本99.免费观看| 在线观看免费视频日本深夜| 国产黄色小视频在线观看| 制服人妻中文乱码| 成人三级黄色视频| 午夜激情福利司机影院| 又黄又爽又免费观看的视频| 久久久久久国产a免费观看| 成年免费大片在线观看| 三级国产精品欧美在线观看 | 日本成人三级电影网站| 国产精品免费一区二区三区在线| 男人和女人高潮做爰伦理| 亚洲熟妇中文字幕五十中出| 一区二区三区激情视频| 久久久国产成人精品二区| 91老司机精品| 日本成人三级电影网站| 久久精品国产亚洲av香蕉五月| 日本a在线网址| 精品久久蜜臀av无| 好看av亚洲va欧美ⅴa在| 国产精品久久久久久亚洲av鲁大| 亚洲精品美女久久久久99蜜臀| 这个男人来自地球电影免费观看| 久久午夜综合久久蜜桃| 神马国产精品三级电影在线观看| 午夜福利视频1000在线观看| 成人国产综合亚洲| 免费观看人在逋| 一级毛片女人18水好多| 九九久久精品国产亚洲av麻豆 | 成人一区二区视频在线观看| 日韩欧美一区二区三区在线观看| 国产成人欧美在线观看| 久久久久久人人人人人| 给我免费播放毛片高清在线观看| 亚洲专区国产一区二区| 极品教师在线免费播放| 啦啦啦韩国在线观看视频| 少妇丰满av| 熟女人妻精品中文字幕| 熟女电影av网| 免费人成视频x8x8入口观看| 欧美乱色亚洲激情| 久久久久久国产a免费观看| 久久久久精品国产欧美久久久| 欧美又色又爽又黄视频| 亚洲五月婷婷丁香| 免费av不卡在线播放| 午夜福利免费观看在线| 91老司机精品| 色吧在线观看| 日韩欧美在线二视频| 国产男靠女视频免费网站| 真人一进一出gif抽搐免费| 亚洲色图av天堂| 一个人免费在线观看的高清视频| 精品一区二区三区av网在线观看| 午夜免费观看网址| 午夜福利成人在线免费观看| 亚洲自拍偷在线| 亚洲成av人片免费观看| 三级男女做爰猛烈吃奶摸视频| 在线观看美女被高潮喷水网站 | 精品一区二区三区视频在线 | 在线永久观看黄色视频| 婷婷丁香在线五月| 午夜福利18| 99在线视频只有这里精品首页| 国产精品1区2区在线观看.| 精品国产乱码久久久久久男人| 淫妇啪啪啪对白视频| 日本 av在线| 免费在线观看亚洲国产| 亚洲 国产 在线| 免费观看的影片在线观看| 在线十欧美十亚洲十日本专区| 国产亚洲欧美在线一区二区| 99精品在免费线老司机午夜| 国产精品免费一区二区三区在线| 国产av一区在线观看免费| 日本与韩国留学比较| 久99久视频精品免费| 亚洲精品一卡2卡三卡4卡5卡| 极品教师在线免费播放| 少妇的逼水好多| 久久中文看片网| 国产v大片淫在线免费观看| 97碰自拍视频| 午夜激情欧美在线| 欧美三级亚洲精品| 国产精华一区二区三区| 久久精品91蜜桃| cao死你这个sao货| 国产毛片a区久久久久| 亚洲性夜色夜夜综合| 国产欧美日韩精品一区二区| 波多野结衣高清作品| 亚洲欧洲精品一区二区精品久久久| 精华霜和精华液先用哪个| 久久久久性生活片| 久久精品人妻少妇| 国产美女午夜福利| 久久精品人妻少妇| 国产成人av教育| 999精品在线视频| 又大又爽又粗| 一级黄色大片毛片| 久久久久国内视频| 在线a可以看的网站| 69av精品久久久久久| 欧美日本视频| 精品久久久久久久毛片微露脸| 成人18禁在线播放| 欧美日韩黄片免| 国产精品99久久久久久久久| 黑人操中国人逼视频| 一a级毛片在线观看| 国产亚洲av嫩草精品影院| 精品99又大又爽又粗少妇毛片 | 成人特级av手机在线观看| 禁无遮挡网站| netflix在线观看网站| 国产成年人精品一区二区| 欧美3d第一页| 欧美一区二区精品小视频在线| 最近最新中文字幕大全免费视频| 怎么达到女性高潮| 婷婷精品国产亚洲av在线| 午夜亚洲福利在线播放| 亚洲国产精品久久男人天堂| ponron亚洲| 免费电影在线观看免费观看| 真人一进一出gif抽搐免费| 亚洲无线观看免费| 亚洲国产精品sss在线观看| 亚洲五月婷婷丁香| 久久久久免费精品人妻一区二区| 99精品欧美一区二区三区四区| 狂野欧美激情性xxxx| 小说图片视频综合网站| 久久久精品欧美日韩精品| 少妇的丰满在线观看| 国产午夜福利久久久久久| 18禁美女被吸乳视频| 窝窝影院91人妻| 国产精品av视频在线免费观看| 成年人黄色毛片网站| 中文资源天堂在线| 日本精品一区二区三区蜜桃| av片东京热男人的天堂| 成人国产综合亚洲| 久久久久久大精品| 最近最新中文字幕大全电影3| 91字幕亚洲| 亚洲国产欧美一区二区综合| 在线国产一区二区在线| 久久精品国产亚洲av香蕉五月| 好男人在线观看高清免费视频| 无限看片的www在线观看| 国产精品久久久人人做人人爽| 国产精品美女特级片免费视频播放器 | 丁香六月欧美| 此物有八面人人有两片| 色播亚洲综合网| 人人妻人人看人人澡| 99国产精品一区二区蜜桃av| 亚洲一区二区三区不卡视频| 日韩欧美免费精品| 757午夜福利合集在线观看| а√天堂www在线а√下载| 超碰成人久久| 两个人视频免费观看高清| 日韩精品中文字幕看吧| 午夜福利在线观看吧| 白带黄色成豆腐渣| 精品久久久久久久末码| 99视频精品全部免费 在线 | 日本黄色视频三级网站网址| 国模一区二区三区四区视频 | 99在线视频只有这里精品首页| 999精品在线视频| 午夜精品久久久久久毛片777| 亚洲电影在线观看av| 欧美成狂野欧美在线观看| 老汉色∧v一级毛片| 午夜精品久久久久久毛片777| 老熟妇仑乱视频hdxx| 国产成人aa在线观看| 午夜福利高清视频| 国产成人一区二区三区免费视频网站| 午夜福利18| 亚洲第一欧美日韩一区二区三区| 久久精品aⅴ一区二区三区四区| 色尼玛亚洲综合影院| 国产精品久久久人人做人人爽| 天堂影院成人在线观看| 一进一出好大好爽视频| 亚洲欧美日韩高清在线视频| 亚洲无线在线观看| 亚洲一区二区三区不卡视频| 亚洲欧美日韩卡通动漫| 国产精品98久久久久久宅男小说| 久久久久国内视频| 12—13女人毛片做爰片一| 女警被强在线播放| 999久久久国产精品视频| 香蕉国产在线看| 三级国产精品欧美在线观看 | 99re在线观看精品视频| 日本三级黄在线观看| 国产私拍福利视频在线观看| 亚洲在线自拍视频| 久久久久久大精品| 日本熟妇午夜| 国产精品 国内视频| 亚洲中文字幕日韩| 男人的好看免费观看在线视频| 三级国产精品欧美在线观看 | 亚洲性夜色夜夜综合| 日韩欧美在线乱码| 日日干狠狠操夜夜爽| 国产精品98久久久久久宅男小说| 国产亚洲av嫩草精品影院| 久久久久久人人人人人| 美女高潮的动态| 国产不卡一卡二| av黄色大香蕉| 国产精品九九99| 国产欧美日韩精品一区二区| 久久中文字幕人妻熟女| 在线十欧美十亚洲十日本专区| 母亲3免费完整高清在线观看| 中文亚洲av片在线观看爽| 亚洲欧美日韩高清专用| 18禁观看日本| 老司机午夜福利在线观看视频| 国产精品久久久av美女十八| 法律面前人人平等表现在哪些方面| 国产91精品成人一区二区三区| 中文字幕久久专区| 国产精品女同一区二区软件 | 欧美zozozo另类| 小说图片视频综合网站| 成人性生交大片免费视频hd| 亚洲成人精品中文字幕电影| 国产亚洲欧美98| 少妇丰满av| 久久人妻av系列| 久久精品影院6| 黑人操中国人逼视频| 亚洲av熟女| 在线观看66精品国产| 欧美成狂野欧美在线观看| www.自偷自拍.com| 午夜福利成人在线免费观看| 美女高潮喷水抽搐中文字幕| 亚洲午夜理论影院| 母亲3免费完整高清在线观看| 亚洲精品在线美女| 午夜激情福利司机影院| 国产精品香港三级国产av潘金莲| 久久精品91蜜桃| 国产伦人伦偷精品视频| 国产视频一区二区在线看| 亚洲国产色片| 美女大奶头视频| 黄色丝袜av网址大全| 日韩三级视频一区二区三区| 又黄又粗又硬又大视频| a在线观看视频网站| 日韩人妻高清精品专区| 成人三级做爰电影| 88av欧美| 亚洲一区二区三区不卡视频| 免费看十八禁软件| www.自偷自拍.com| 蜜桃久久精品国产亚洲av| 国产精品一区二区三区四区久久| 好看av亚洲va欧美ⅴa在| 国产成人精品久久二区二区91| 全区人妻精品视频| aaaaa片日本免费| 99re在线观看精品视频| 亚洲av成人精品一区久久| 舔av片在线| 国产精品99久久99久久久不卡| 最近在线观看免费完整版| 亚洲精品在线观看二区| 丰满人妻熟妇乱又伦精品不卡| 九色成人免费人妻av| netflix在线观看网站| 老鸭窝网址在线观看| 成人av一区二区三区在线看| 99久久综合精品五月天人人| 日本精品一区二区三区蜜桃| 欧美日韩国产亚洲二区| 午夜免费激情av| 国产乱人视频| 狂野欧美激情性xxxx| 国产精品,欧美在线| 亚洲av成人一区二区三| a级毛片在线看网站| 亚洲国产高清在线一区二区三| 国产日本99.免费观看| 99久久精品热视频| 美女被艹到高潮喷水动态| 亚洲精品中文字幕一二三四区| 久久精品91蜜桃| 国产亚洲欧美98| 国产91精品成人一区二区三区| 色综合婷婷激情| a级毛片在线看网站| 999久久久国产精品视频| 欧美日韩亚洲国产一区二区在线观看| 国产成年人精品一区二区| 中文资源天堂在线| 十八禁人妻一区二区| 午夜免费成人在线视频| 日韩欧美国产在线观看| 成在线人永久免费视频| 中文亚洲av片在线观看爽| 在线看三级毛片| 精品国产美女av久久久久小说| 久久久久久久久中文| 最近最新免费中文字幕在线| 国产三级中文精品| 伦理电影免费视频| 日韩三级视频一区二区三区| 亚洲午夜精品一区,二区,三区| 国产真实乱freesex| 国产一区二区在线观看日韩 | 一级毛片高清免费大全| 日本 欧美在线| 国产三级中文精品| 国产高清videossex| 嫩草影院精品99| 亚洲欧美精品综合一区二区三区| 真人一进一出gif抽搐免费| 亚洲成人精品中文字幕电影| 淫妇啪啪啪对白视频| 性色av乱码一区二区三区2| 日本黄大片高清| 亚洲欧美精品综合久久99| 在线免费观看的www视频| 国产成+人综合+亚洲专区| 18禁裸乳无遮挡免费网站照片| 亚洲真实伦在线观看| 国产 一区 欧美 日韩| 日韩欧美免费精品| 亚洲中文字幕一区二区三区有码在线看 | 亚洲精品一卡2卡三卡4卡5卡| 午夜激情欧美在线| 又紧又爽又黄一区二区| 又大又爽又粗| 国产午夜福利久久久久久| 久久久久久大精品| 国产人伦9x9x在线观看| 国产精品久久久人人做人人爽| 黄色丝袜av网址大全| 亚洲,欧美精品.| 精品久久久久久久久久免费视频| 黄色女人牲交| 国产综合懂色| 亚洲国产欧洲综合997久久,| 亚洲av片天天在线观看| 国产伦人伦偷精品视频| 日韩欧美国产在线观看| 国产三级在线视频| 美女高潮的动态| 2021天堂中文幕一二区在线观| 桃红色精品国产亚洲av| 亚洲av成人不卡在线观看播放网| 亚洲精品456在线播放app | 色噜噜av男人的天堂激情| 国内久久婷婷六月综合欲色啪| 香蕉av资源在线| 国产探花在线观看一区二区| 舔av片在线| 亚洲自偷自拍图片 自拍| 国产成人影院久久av| 麻豆久久精品国产亚洲av| 成人av在线播放网站| 国产野战对白在线观看| 99久久成人亚洲精品观看| 亚洲,欧美精品.| 亚洲在线自拍视频| or卡值多少钱| 婷婷亚洲欧美| 久久久色成人| 久久欧美精品欧美久久欧美| 亚洲欧美激情综合另类| 久久人人精品亚洲av| 丰满人妻熟妇乱又伦精品不卡| 一二三四社区在线视频社区8| 99国产精品99久久久久| 亚洲中文av在线| 亚洲av成人一区二区三| av视频在线观看入口| 亚洲熟女毛片儿| 精品一区二区三区视频在线 | 精品人妻1区二区| 一卡2卡三卡四卡精品乱码亚洲| 97超视频在线观看视频| 真人做人爱边吃奶动态|