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

    四級(jí)公路貨車路側(cè)事故概率預(yù)測(cè)與多發(fā)點(diǎn)判別方法研究

    2024-01-03 04:43:14張麗莉靳鋼睿程國(guó)柱徐亮
    關(guān)鍵詞:縱坡坡度貨車

    張麗莉,靳鋼睿,程國(guó)柱,徐亮

    (1.東北林業(yè)大學(xué) 交通學(xué)院,黑龍江 哈爾濱 150040;2.北京汽車研究總院有限公司,北京 101300; 3.長(zhǎng)春工程學(xué)院 土木學(xué)院,吉林 長(zhǎng)春 130012)

    隨著我國(guó)公路貨運(yùn)的快速發(fā)展,涉及貨車的路側(cè)事故占比也逐年增大。美國(guó)的路側(cè)安全研究項(xiàng)目數(shù)據(jù)顯示,由交通事故造成的死亡人數(shù)里一半以上都是由路側(cè)事故造成的[1]。中國(guó)《道路交通事故統(tǒng)計(jì)年報(bào)》顯示:2018年中國(guó)貨運(yùn)車輛發(fā)生了約4萬起的道路交通事故,占全年總事故數(shù)的21%,這些事故造成了約2萬的死亡人數(shù),占全年事故死亡人數(shù)的22%,導(dǎo)致約5億元的直接經(jīng)濟(jì)損失[2]。截止到2021年末,四級(jí)公路的公路里程占全國(guó)公路總里程的76.4%。由于四級(jí)公路的路面狀況和圓曲線半徑、縱坡等線形條件較差,以及貨車車身長(zhǎng)、駕駛性能差、重量大等相關(guān)因素,貨車參與的事故中重大、特大事故較多,而且路側(cè)事故在發(fā)生的事故中所占比例較大、造成的死亡率較高。

    路側(cè)事故致因的相關(guān)研究表明,Lord等[3]認(rèn)為平曲線路段側(cè)翻事故發(fā)生率隨路肩寬度的增大而減小。Rusli等[4]認(rèn)為雨天和縱坡坡度較大容易導(dǎo)致路側(cè)事故的發(fā)生,設(shè)置公路輪廓線和寬路肩能明顯降低路側(cè)事故概率。程國(guó)柱等[5]認(rèn)為車速對(duì)平曲線路段發(fā)生路側(cè)事故的影響很大。Cheng等[6]通過分析路側(cè)事故發(fā)生頻率與護(hù)欄、路側(cè)固定物間距、邊坡坡度等之間的關(guān)系來研究風(fēng)險(xiǎn)因素。Butsick等[7]用負(fù)二項(xiàng)式回歸建立了兩種低等級(jí)公路的安全性能函數(shù),使用經(jīng)驗(yàn)貝葉斯調(diào)整后提出了一種適用于特定風(fēng)險(xiǎn)的交通事故多發(fā)點(diǎn)鑒別方法。Briz-Redon等[8]利用不同碰撞類型的事故風(fēng)險(xiǎn)和風(fēng)險(xiǎn)概率,提出一種適用于不同碰撞類型的事故高風(fēng)險(xiǎn)區(qū)域識(shí)別方法。袁泉等[9]構(gòu)建了鑒別事故風(fēng)險(xiǎn)點(diǎn)的模型,對(duì)道路交通事故風(fēng)險(xiǎn)點(diǎn)進(jìn)行多層面鑒別分析。孟祥海等[10]利用統(tǒng)計(jì)和假設(shè)檢驗(yàn)的方法識(shí)別事故顯著性風(fēng)險(xiǎn)因素,鑒別分析了高速公路路段上的事故多發(fā)點(diǎn)。

    早期的研究者采用對(duì)數(shù)正態(tài)回歸模型對(duì)路側(cè)事故頻率進(jìn)行研究。然而,在事故預(yù)測(cè)方面,這種傳統(tǒng)線性回歸模型的預(yù)測(cè)精度較低。隨著后繼交通安全科研人員的不斷研究,事故預(yù)測(cè)中存在大量零值的問題可由泊松和負(fù)二項(xiàng)回歸模型解決[11]。馬壯林等[12]認(rèn)為用零堆積負(fù)二項(xiàng)回歸模型可以更好地預(yù)測(cè)路側(cè)事故數(shù)。盡管有多種統(tǒng)計(jì)回歸方法可以對(duì)路側(cè)事故數(shù)量或頻率進(jìn)行預(yù)測(cè),但對(duì)預(yù)測(cè)事故概率卻并不適用。考慮到與事故頻率和事故數(shù)相比,事故概率更能體現(xiàn)路段事故的高發(fā)水平,因而對(duì)事故概率預(yù)測(cè)進(jìn)行研究是有必要的。綜上所述,本文使用PC-Crash軟件進(jìn)行仿真試驗(yàn)來獲取研究數(shù)據(jù),提出了四級(jí)公路貨車路側(cè)事故概率預(yù)測(cè)模型和事故多發(fā)點(diǎn)判別方法。研究成果可為交通管理部門鑒別路側(cè)事故風(fēng)險(xiǎn)要素提供依據(jù),有助于減少公路貨車路側(cè)事故。

    1 交通事故模擬技術(shù)

    大部分學(xué)者對(duì)事故進(jìn)行統(tǒng)計(jì)分析和研究時(shí)常依賴歷史事故數(shù)據(jù)。但是隨著時(shí)間的變化,事故現(xiàn)場(chǎng)的道路交通環(huán)境會(huì)發(fā)生改變,再加上事故資料的殘缺,研究成果的精準(zhǔn)性和有效性會(huì)受這些情況的影響。因此,近年來交通事故模擬技術(shù)被廣泛地應(yīng)用于交通事故數(shù)據(jù)獲取。

    目前常用的事故模擬技術(shù)有駕駛模擬器、MADYMO、LS—DYNA、Carsim、PC-Crash等。本文采用PC-Crash軟件來收集路側(cè)交通事故數(shù)據(jù)。相關(guān)文獻(xiàn)中的國(guó)內(nèi)外學(xué)者使用PC-Crash軟件進(jìn)行路側(cè)事故致因分析和事故嚴(yán)重度研究,并進(jìn)行了案例驗(yàn)證,證明了該軟件在模擬單車側(cè)翻事故和與固定物體的碰撞事故方面具有高度的可靠性和獨(dú)特的優(yōu)越性[5-6,13]。此外,與收集的事故數(shù)據(jù)相比,利用PC-Crash軟件獲取的數(shù)據(jù)具有以下優(yōu)點(diǎn):

    (1)事故現(xiàn)場(chǎng)的信息更加完善。

    (2)事故數(shù)據(jù)不會(huì)因?yàn)闀r(shí)間和交通等條件變化而受到影響。

    (3)區(qū)域特性不會(huì)對(duì)事故數(shù)據(jù)造成影響,獲得的研究成果具有普遍適用性。

    (4)研究成本低。

    PC-Crash仿真軟件進(jìn)行數(shù)字化重構(gòu)事故過程時(shí)使用的是優(yōu)化計(jì)算的方式。使用者首先調(diào)入所需的車型和道路模型,隨后進(jìn)行模型的參數(shù)設(shè)置,包括道路模型的幾何條件、車輛模型的尺寸和位置等參數(shù)估計(jì)值。將優(yōu)化變量設(shè)置為所需車型的終止位置或制動(dòng)軌跡,將車輛行駛的中間位置作為約束條件,經(jīng)不斷迭代計(jì)算,反復(fù)調(diào)整事故車輛運(yùn)動(dòng)參數(shù),如碰撞角度、加速度、減速度、初始速度等初始值,以及碰撞點(diǎn)位置、道路附著系數(shù)等參數(shù),一直到優(yōu)化所得結(jié)果與實(shí)際事故狀態(tài)基本一致,由此可得到最趨近于實(shí)際事故狀況的車輛運(yùn)動(dòng)參數(shù)。

    PC-Crash軟件通過調(diào)查、收集和分析研究事故信息、資料和數(shù)據(jù)來實(shí)現(xiàn)交通事故再現(xiàn)與分析,將事故車輛從發(fā)生事故后的終點(diǎn)位置反向推演到事故發(fā)生的整個(gè)過程,再將其反向推演到事故發(fā)生之前的車輛行駛狀況,從而對(duì)事故的成因進(jìn)行分析。目前,我國(guó)在交通事故分析再現(xiàn)和事故責(zé)任認(rèn)定方面主要依賴于人工的估算和分析判斷,但是人的主觀因素會(huì)對(duì)其造成很大的影響,通過事故仿真與重現(xiàn)技術(shù)對(duì)道路交通事故的再現(xiàn)模擬,既可為事故責(zé)任認(rèn)定提供依據(jù),又可輔助研究事故的一般成因,具有很高的研究?jī)r(jià)值和現(xiàn)實(shí)意義。

    2 數(shù)據(jù)收集

    根據(jù)路側(cè)事故調(diào)查經(jīng)驗(yàn),并結(jié)合相關(guān)研究成果,本文將車輛駛離行車道和未駛離行車道作為因變量,將圓曲線半徑、車速、縱坡坡度、超高橫坡度、路面附著系數(shù)、車型6個(gè)風(fēng)險(xiǎn)要素作為自變量,其中車輛駛離行車道代表出現(xiàn)路側(cè)事故,未駛離行車道代表未出現(xiàn)路側(cè)事故。

    國(guó)家統(tǒng)計(jì)局的統(tǒng)計(jì)結(jié)果顯示,2020年輕型民用載貨汽車和重型民用載貨汽車的擁有量占民用載貨汽車擁有量的96%,其中輕型民用載貨汽車載重量大于1 800 kg同時(shí)小于等于6 000 kg,重型民用載貨汽車載重量大于14 000 kg。因此本文車型選取重型貨車(重型民用載貨汽車)和輕型貨車(輕型民用載貨汽車)。

    本文使用PC-Crash仿真軟件構(gòu)建了3種公路模型。第一種是30 km/h的設(shè)計(jì)速度、3.25 m的行車道,0.5 m的土路肩寬度、4 m的路基高度,邊坡坡度為1∶1.5的雙向兩車道公路模型;第二種是20 km/h的設(shè)計(jì)速度、3 m的行車道,0.25 m的土路肩寬度、4 m的路基高度,邊坡坡度為1∶1.5的雙向兩車道公路模型;第三種是20 km/h的設(shè)計(jì)速度、3.5 m的行車道,0.5 m的土路肩寬度、4 m的路基高度,邊坡坡度為1∶1.5的單向單車道公路模型。從軟件中調(diào)取MAN-19.362 作為輕型貨車的代表車型,調(diào)取ASCHERSLEBEN KAROSS作為重型貨車的代表車型,見圖1。

    圖1 試驗(yàn)場(chǎng)景

    按照我國(guó)《公路路線設(shè)計(jì)規(guī)范》(JTG D20-2017)[14]有關(guān)規(guī)定,第一種公路模型中每個(gè)風(fēng)險(xiǎn)要素分別取為:2%、4%、6%的縱坡坡度;4%、6%、8%的超高橫坡度;30、35、40、65、100、200 m的圓曲線半徑;0.2、0.4、0.6、0.8的路面附著系數(shù);20、30、40、60、80 km/h的車速。第二種和第三種公路模型中每個(gè)風(fēng)險(xiǎn)要素分別取為:2%、4%、6%的縱坡坡度;4%、6%、8%的超高橫坡度;15、20、30、65、100、200 m的圓曲線半徑;0.2、0.4、0.6、0.8的路面附著系數(shù);20、30、40、60、80 km/h的車速。在仿真試驗(yàn)中需要對(duì)圓曲線半徑進(jìn)行加寬,根據(jù)規(guī)范[14],貨車的加寬值見表1,其中單車道公路路面加寬值為雙車道公路路面加寬值的一半。由于駕駛?cè)嗽谟龅綇澋缆范螘r(shí)會(huì)進(jìn)行轉(zhuǎn)向操作,所以在每次開展試驗(yàn)前,根據(jù)不同的圓曲線半徑,提前設(shè)置車輛的轉(zhuǎn)彎半徑,完成車輛轉(zhuǎn)向輪的偏轉(zhuǎn)。

    表1 雙車道路面加寬值 m

    根據(jù)每種風(fēng)險(xiǎn)要素的取值,在3種公路模型中分別建立了180種道路情況,對(duì)車輛設(shè)定5種初始速度,并對(duì)圓曲線進(jìn)行加寬,對(duì)重型貨車和輕型貨車分別進(jìn)行仿真試驗(yàn),分別得到1 800組仿真數(shù)據(jù)。

    第一種公路模型中未出現(xiàn)路側(cè)事故1 491組,占82.8%,出現(xiàn)路側(cè)事故309組,占17.2%。第二種公路模型中未出現(xiàn)路側(cè)事故1 473組,占81.8%,出現(xiàn)路側(cè)事故327組,占18.2%。第三種公路模型中未出現(xiàn)路側(cè)事故1 490組,占82.8%,出現(xiàn)路側(cè)事故310組,占17.2%。

    3 四級(jí)公路貨車路側(cè)事故概率預(yù)測(cè)

    3.1 Logistic回歸分析

    Logistic回歸屬于廣義線性回歸模型,廣泛用于收集分析數(shù)據(jù),對(duì)給定事件進(jìn)行預(yù)估等范疇。由于未出現(xiàn)路側(cè)事故和出現(xiàn)路側(cè)事故這兩個(gè)因變量是二分類,所以對(duì)路側(cè)事故概率進(jìn)行預(yù)測(cè)時(shí)選取二元Logistic回歸模型。假定因變量與n個(gè)自變量相關(guān),記X=(x1,x2,...,xn),則二元Logistic回歸模型[15]為:

    (1)

    式中:a1為偏回歸系數(shù);a0為常數(shù)項(xiàng);P為事件發(fā)生概率。

    因模型中各自變量的單位不同,偏回歸系數(shù)不能直接體現(xiàn)其對(duì)因變量的重要程度大小,故采用標(biāo)準(zhǔn)回歸系數(shù)來檢驗(yàn)自變量對(duì)因變量的重要程度大小。標(biāo)準(zhǔn)回歸系數(shù)αι′的計(jì)算公式[16]為:

    αι′=αi(Si/S)

    (2)

    將研究數(shù)據(jù)調(diào)入二元Logistic回歸模型中,分析結(jié)果見表2。可以看出,每個(gè)變量的p值都小于0.05,表明四級(jí)公路貨車路側(cè)事故是否發(fā)生與這些變量顯著相關(guān)。

    將表2中的偏回歸系數(shù)和標(biāo)準(zhǔn)差代入式(2), 可求得各個(gè)風(fēng)險(xiǎn)要素的標(biāo)準(zhǔn)回歸系數(shù)。由表2可以看出顯著性風(fēng)險(xiǎn)要素圓曲線半徑對(duì)雙車道四級(jí)公路貨車路側(cè)事故的重要程度最大,其次為車速、路面附著系數(shù)、車型、縱坡坡度,超高橫坡度的重要程度最小。由表2可以看出顯著性風(fēng)險(xiǎn)要素圓曲線半徑對(duì)單車道四級(jí)公路貨車路側(cè)事故的重要程度最大,其次為車速、路面附著系數(shù)、車型、超高橫坡度,縱坡坡度的重要程度最小。

    表2 三種公路模型分析結(jié)果

    表2(續(xù))

    根據(jù)表2中各個(gè)變量的偏回歸系數(shù),可構(gòu)建第一種公路模型的貨車路側(cè)事故概率預(yù)測(cè)模型為:

    (3)

    式中:v為車速,km/h;R為圓曲線半徑,m;μ為路面附著系數(shù);i1為縱坡坡度,%(下坡為正);i2為超高橫坡度,%(當(dāng)試驗(yàn)段橫斷面的外部高度大于內(nèi)部高度時(shí),其值為正);η為車型,輕型貨車取0,重型貨車取1。

    根據(jù)表2中各個(gè)變量的偏回歸系數(shù),可構(gòu)建第二種公路模型的貨車路側(cè)事故概率預(yù)測(cè)模型為:

    (4)

    根據(jù)表2中各個(gè)變量的偏回歸系數(shù),可構(gòu)建第三種公路模型的貨車路側(cè)事故概率預(yù)測(cè)模型為:

    (5)

    3.2 模型檢驗(yàn)

    本文對(duì)模型進(jìn)行檢驗(yàn)時(shí)使用的是ROC曲線,此方法用來檢驗(yàn)?zāi)P蛿M合效果時(shí)通常依據(jù)的是曲線下方面積。模型的擬合效果在置信度為95%時(shí),隨AUC值的增大而增強(qiáng)。AUC的值通常在0.5~1.0范圍內(nèi)。在預(yù)測(cè)模型結(jié)果的基礎(chǔ)上,使用SPSS軟件生成了3種模型的ROC曲線,見圖2。

    (a) 第一種

    (b) 第二種

    (c) 第三種圖2 3種模型的ROC曲線

    圖2(a)中的曲線極其接近左上方,AUC的值為0.992(大于0.9),因此第一種公路模型的路側(cè)事故概率預(yù)測(cè)模型在置信度為95%時(shí)比較精準(zhǔn)。

    圖2(b)和圖2(c)中的曲線極其接近左上方,AUC的值為0.994(大于0.9),因此第二種和第三種公路模型的路側(cè)事故概率預(yù)測(cè)模型在置信度為95%時(shí)比較精準(zhǔn)。

    根據(jù)上文得到的概率預(yù)測(cè)模型,綜合考慮路面附著系數(shù)、縱坡坡度、車型等風(fēng)險(xiǎn)要素,從而獲取各個(gè)路段單元可能發(fā)生的貨車路側(cè)事故概率。其中四級(jí)公路各個(gè)路段單元的運(yùn)行速度可由JTG B05—2015[17]中公路運(yùn)行速度計(jì)算模型測(cè)算得到。

    4 四級(jí)公路貨車路側(cè)事故多發(fā)點(diǎn)判別

    根據(jù)上文3種模型的預(yù)測(cè)結(jié)果,可以看出:第一種模型中1 500組概率值小于50%,占總數(shù)的83.3%,300組概率值大于50%,占總數(shù)的16.7%。第二種模型中1 479組概率值小于50%,占總數(shù)的82.2%,321組概率值大于50%,占總數(shù)的17.8%。第三種模型中1497組概率值小于50%,占總數(shù)的83.2%,303組概率值大于50%,占總數(shù)的16.8%。因此,在概率預(yù)測(cè)總結(jié)果中,路側(cè)事故低概率事件占比較大,路側(cè)事故高概率事件占比較小。

    根據(jù)上述分析,本文使用了累計(jì)頻率法去識(shí)別四級(jí)公路貨車路側(cè)事故多發(fā)路段。此方法認(rèn)為路側(cè)事故低概率事件的路段占大部分,路側(cè)事故高概率事件的路段占少部分,并且路側(cè)事故發(fā)生概率越高的路段所占比例越小。

    對(duì)3種模型中的概率值進(jìn)行由小到大的排序,計(jì)算得到每個(gè)概率值出現(xiàn)的累計(jì)頻率,由此可分離出累計(jì)頻率很小,但事故風(fēng)險(xiǎn)很高的概率閾值。因此將縱坐標(biāo)設(shè)置為對(duì)應(yīng)路側(cè)事故概率的累計(jì)頻率,將橫坐標(biāo)設(shè)置為路側(cè)事故概率。繪制了累計(jì)頻率關(guān)于概率值的變化曲線,見圖3。

    (a) 第一種

    (b) 第二種

    (c) 第三種圖3 3種模型累計(jì)頻率曲線

    圖3(a)中累計(jì)頻率隨著概率值的增大而增大的趨勢(shì)逐漸平緩。通過SPSS軟件,再結(jié)合曲線的變化趨勢(shì),擬合的曲線方程為:

    y=-1.145x2+1.948x+0.131

    (6)

    式中:x為事故概率;y為累計(jì)頻率。

    曲線方程的擬合度R2=0.99。依據(jù)“累計(jì)頻率曲線法”鑒別規(guī)則[18],經(jīng)過曲線方程求導(dǎo)得到曲線拐點(diǎn)處的概率值為85%。即當(dāng)設(shè)計(jì)速度為30 km/h的雙車道四級(jí)公路貨車路側(cè)事故概率∈[85%, 100%]時(shí),判定路段為路側(cè)事故潛在多發(fā)點(diǎn)。

    圖3(b)中累計(jì)頻率隨著概率值的增大而增大的趨勢(shì)逐漸平緩。通過SPSS軟件,再結(jié)合曲線的變化趨勢(shì),擬合的曲線方程為:

    y=-1.065x2+1.893x+0.115

    (7)

    曲線方程的擬合度R2=0.995。曲線拐點(diǎn)處的概率值為89%。即當(dāng)設(shè)計(jì)速度為20 km/h的雙車道四級(jí)公路貨車路側(cè)事故概率∈[89%, 100%]時(shí),判定路段為路側(cè)事故潛在多發(fā)點(diǎn)。

    圖3(c)中累計(jì)頻率隨著概率值的增大而增大的趨勢(shì)逐漸平緩。通過SPSS軟件,再結(jié)合曲線的變化趨勢(shì),擬合的曲線方程為:

    y=-1.121x2+1.975x+0.105

    (8)

    曲線方程的擬合度R2=0.996。曲線拐點(diǎn)處的概率值為88%。即當(dāng)設(shè)計(jì)速度為20 km/h的單車道四級(jí)公路貨車路側(cè)事故概率∈[88%, 100%]時(shí),判定路段為路側(cè)事故潛在多發(fā)點(diǎn)。

    5 結(jié)論

    本文得到的主要結(jié)論如下:

    (1)基于二元 Logistic 回歸方法,分析識(shí)別了四級(jí)公路貨車路側(cè)事故的顯著性風(fēng)險(xiǎn)要素,探討了各個(gè)顯著性風(fēng)險(xiǎn)要素對(duì)路側(cè)事故的重要程度,并建立了四級(jí)公路貨車路側(cè)事故概率預(yù)測(cè)模型,該模型可以計(jì)算各個(gè)顯著性風(fēng)險(xiǎn)要素條件下的路側(cè)事故概率,對(duì)路側(cè)事故概率實(shí)現(xiàn)了定量分析。

    (2)對(duì)雙車道四級(jí)公路貨車路側(cè)事故重要程度影響最大的是圓曲線半徑,其次為車速、路面附著系數(shù)、車型、縱坡坡度,超高橫坡度的重要程度最小。對(duì)單車道公路貨車路側(cè)事故重要程度影響最大的也是圓曲線半徑,其次為車速、路面附著系數(shù)、車型、超高橫坡度,縱坡坡度的重要程度最小。

    (3)基于四級(jí)公路貨車路側(cè)事故概率預(yù)測(cè)模型,得到3種四級(jí)公路的貨車路側(cè)事故潛在多發(fā)點(diǎn)概率閾值,分別為85%、89%和88%,提出了四級(jí)公路貨車路側(cè)事故多發(fā)點(diǎn)判別方法。

    本文對(duì)路側(cè)事故致因研究分析后發(fā)現(xiàn),除文中得到的顯著性風(fēng)險(xiǎn)因素外,還存在其他對(duì)事故有顯著影響的風(fēng)險(xiǎn)因素,如駕駛?cè)似诔潭?、光線情況、天氣情況等,后續(xù)研究可將多源數(shù)據(jù)融合,提出范圍更廣、適用性更強(qiáng)的事故概率預(yù)測(cè)模型,并完成實(shí)例驗(yàn)證。

    猜你喜歡
    縱坡坡度貨車
    貨車
    幼兒畫刊(2023年12期)2024-01-15 07:06:14
    基于交通事故多發(fā)位置的區(qū)間平均縱坡控制指標(biāo)研究
    半柔性復(fù)合路面在重交通長(zhǎng)大縱坡路段的應(yīng)用
    石油瀝青(2020年1期)2020-05-25 06:54:08
    道路最小設(shè)計(jì)縱坡探討
    關(guān)于公路超高漸變段合成坡度解析與應(yīng)用
    揭示車路協(xié)同矛盾 破解長(zhǎng)大縱坡難題
    貨車也便捷之ETC新時(shí)代!——看高速公路貨車ETC如何實(shí)現(xiàn)
    學(xué)與玩(2017年6期)2017-02-16 07:07:24
    基于圖像處理的定位器坡度計(jì)算
    電氣化鐵道(2016年4期)2016-04-16 05:59:46
    治超新規(guī)實(shí)施在即 深究貨車非法改裝亂象
    專用汽車(2016年9期)2016-03-01 04:16:52
    一本—道久久a久久精品蜜桃钙片| 老女人水多毛片| 日韩成人伦理影院| 日韩一本色道免费dvd| av免费在线看不卡| 亚洲va在线va天堂va国产| 国产69精品久久久久777片| 久久久精品免费免费高清| 欧美日韩一区二区视频在线观看视频在线| 精华霜和精华液先用哪个| 噜噜噜噜噜久久久久久91| 最近手机中文字幕大全| 亚洲av在线观看美女高潮| 久久鲁丝午夜福利片| 伊人亚洲综合成人网| 久久久久久伊人网av| 啦啦啦中文免费视频观看日本| 国产日韩欧美视频二区| 亚洲国产精品999| 亚洲激情五月婷婷啪啪| 欧美另类一区| 国产视频内射| 人人妻人人爽人人添夜夜欢视频 | 日韩三级伦理在线观看| 99久国产av精品国产电影| 精品人妻一区二区三区麻豆| 女性生殖器流出的白浆| 亚洲怡红院男人天堂| av国产精品久久久久影院| 欧美3d第一页| 九草在线视频观看| 国产一区二区在线观看av| 免费观看性生交大片5| 韩国高清视频一区二区三区| 国产伦在线观看视频一区| 汤姆久久久久久久影院中文字幕| 成人漫画全彩无遮挡| 在线免费观看不下载黄p国产| 美女中出高潮动态图| 寂寞人妻少妇视频99o| 97超视频在线观看视频| 麻豆成人午夜福利视频| 18+在线观看网站| 91午夜精品亚洲一区二区三区| 一级片'在线观看视频| av不卡在线播放| 观看免费一级毛片| 国产精品一区二区性色av| 黑人高潮一二区| 我要看黄色一级片免费的| 亚洲国产精品国产精品| 日本-黄色视频高清免费观看| 亚洲av.av天堂| 夜夜骑夜夜射夜夜干| 亚洲国产精品成人久久小说| kizo精华| 国产熟女午夜一区二区三区 | 另类精品久久| 国产精品.久久久| 26uuu在线亚洲综合色| 一区二区三区四区激情视频| 又爽又黄a免费视频| av天堂中文字幕网| 午夜激情久久久久久久| 五月天丁香电影| 成人综合一区亚洲| 免费高清在线观看视频在线观看| 免费观看性生交大片5| 欧美日韩亚洲高清精品| 成年人免费黄色播放视频 | 久久久久国产精品人妻一区二区| 蜜桃在线观看..| 一区二区三区精品91| 成人18禁高潮啪啪吃奶动态图 | 中文字幕精品免费在线观看视频 | 亚洲国产成人一精品久久久| 日日啪夜夜爽| 中文字幕久久专区| 国产av一区二区精品久久| 国产探花极品一区二区| 黄色配什么色好看| 国产免费福利视频在线观看| 黄片无遮挡物在线观看| 亚洲天堂av无毛| 最近最新中文字幕免费大全7| 亚洲精品国产色婷婷电影| 国产精品久久久久久av不卡| 亚洲人与动物交配视频| 亚洲久久久国产精品| av福利片在线| 2021少妇久久久久久久久久久| 丝袜喷水一区| 2022亚洲国产成人精品| 卡戴珊不雅视频在线播放| 99热这里只有精品一区| 国产成人精品久久久久久| 免费久久久久久久精品成人欧美视频 | 午夜免费男女啪啪视频观看| 五月伊人婷婷丁香| 高清不卡的av网站| 涩涩av久久男人的天堂| 成人毛片60女人毛片免费| 国产亚洲午夜精品一区二区久久| 亚洲精品456在线播放app| 国产成人精品福利久久| 爱豆传媒免费全集在线观看| 狂野欧美白嫩少妇大欣赏| 日韩一区二区视频免费看| 亚洲一区二区三区欧美精品| 日日摸夜夜添夜夜爱| 久久精品国产鲁丝片午夜精品| 美女脱内裤让男人舔精品视频| 97超视频在线观看视频| 欧美日韩av久久| 亚洲精华国产精华液的使用体验| 婷婷色综合大香蕉| 国产免费一区二区三区四区乱码| 精品亚洲成国产av| 久久热精品热| 亚洲欧美日韩卡通动漫| 在线观看一区二区三区激情| 国产精品偷伦视频观看了| 激情五月婷婷亚洲| 午夜免费鲁丝| 美女脱内裤让男人舔精品视频| 免费观看av网站的网址| 最新中文字幕久久久久| 午夜激情久久久久久久| 男男h啪啪无遮挡| 在线观看www视频免费| 日韩一区二区视频免费看| 国产免费一级a男人的天堂| 99久久中文字幕三级久久日本| 亚洲,一卡二卡三卡| 欧美精品亚洲一区二区| 观看免费一级毛片| 看十八女毛片水多多多| 69精品国产乱码久久久| 人人澡人人妻人| 国产成人aa在线观看| 亚洲国产精品国产精品| 另类亚洲欧美激情| 大话2 男鬼变身卡| 菩萨蛮人人尽说江南好唐韦庄| 国产日韩一区二区三区精品不卡 | 欧美三级亚洲精品| 国产精品无大码| 一边亲一边摸免费视频| 下体分泌物呈黄色| 国产日韩欧美亚洲二区| 少妇高潮的动态图| 欧美日韩视频精品一区| 一级a做视频免费观看| 99久久综合免费| 青春草国产在线视频| 一个人免费看片子| 三级经典国产精品| 婷婷色综合大香蕉| 成人18禁高潮啪啪吃奶动态图 | 国产成人91sexporn| 亚洲av电影在线观看一区二区三区| 2021少妇久久久久久久久久久| 中文字幕制服av| 欧美高清成人免费视频www| 国产av精品麻豆| 国产 一区精品| 成人影院久久| 五月天丁香电影| 亚洲三级黄色毛片| 日本黄大片高清| 色吧在线观看| 在线观看免费高清a一片| 亚洲欧美日韩另类电影网站| 成人美女网站在线观看视频| 国产精品99久久久久久久久| 天天操日日干夜夜撸| 最近最新中文字幕免费大全7| √禁漫天堂资源中文www| 好男人视频免费观看在线| 久久久久久久国产电影| 两个人免费观看高清视频 | 男女边吃奶边做爰视频| 中文字幕精品免费在线观看视频 | 极品人妻少妇av视频| 亚洲国产精品一区三区| 丰满饥渴人妻一区二区三| 久久久久久久久久成人| 亚洲精品自拍成人| 精品一区二区免费观看| 中国国产av一级| a 毛片基地| 亚洲欧美精品自产自拍| 91午夜精品亚洲一区二区三区| 777米奇影视久久| 最后的刺客免费高清国语| 丝瓜视频免费看黄片| 精品国产露脸久久av麻豆| 精品久久久久久电影网| 人妻人人澡人人爽人人| 国产乱人偷精品视频| 中文资源天堂在线| 高清午夜精品一区二区三区| 少妇人妻一区二区三区视频| 在线观看三级黄色| 少妇裸体淫交视频免费看高清| 久久国内精品自在自线图片| 国产一级毛片在线| 在线观看人妻少妇| 亚洲国产精品999| 中文字幕人妻丝袜制服| av播播在线观看一区| 精品一区二区免费观看| 国产精品女同一区二区软件| 亚洲欧美清纯卡通| 国产精品久久久久久精品古装| 青春草视频在线免费观看| 久久鲁丝午夜福利片| 免费观看在线日韩| 26uuu在线亚洲综合色| 亚洲一级一片aⅴ在线观看| 校园人妻丝袜中文字幕| 欧美日韩av久久| 国产高清三级在线| 国产伦在线观看视频一区| 国产精品一区二区三区四区免费观看| 99视频精品全部免费 在线| 蜜臀久久99精品久久宅男| 一本色道久久久久久精品综合| av国产精品久久久久影院| 麻豆成人午夜福利视频| 久久99热6这里只有精品| 新久久久久国产一级毛片| 欧美高清成人免费视频www| 性高湖久久久久久久久免费观看| 日韩av不卡免费在线播放| 国产成人91sexporn| 视频中文字幕在线观看| 麻豆成人午夜福利视频| 久久久国产欧美日韩av| 国产精品99久久久久久久久| 91aial.com中文字幕在线观看| 国产成人精品婷婷| 黄色视频在线播放观看不卡| 中文资源天堂在线| 日韩欧美 国产精品| 日本wwww免费看| 91久久精品国产一区二区成人| 国产视频内射| 国产又色又爽无遮挡免| 交换朋友夫妻互换小说| 啦啦啦在线观看免费高清www| 91久久精品国产一区二区成人| 韩国av在线不卡| 在线观看免费视频网站a站| 肉色欧美久久久久久久蜜桃| 久久久久久久久大av| 国模一区二区三区四区视频| 狠狠精品人妻久久久久久综合| 久久av网站| 久久久久久久久久久久大奶| 亚洲国产欧美在线一区| 婷婷色综合www| 丁香六月天网| 欧美区成人在线视频| 日韩不卡一区二区三区视频在线| 久久ye,这里只有精品| 国产片特级美女逼逼视频| av在线app专区| 国产男女超爽视频在线观看| 精品久久久久久电影网| 国产精品久久久久久av不卡| 美女视频免费永久观看网站| 欧美变态另类bdsm刘玥| 91久久精品电影网| 久久久久视频综合| 2021少妇久久久久久久久久久| 一级爰片在线观看| 欧美三级亚洲精品| 久久国产乱子免费精品| 一级毛片久久久久久久久女| 99热这里只有是精品在线观看| 熟女av电影| 亚洲精品久久久久久婷婷小说| 男女边吃奶边做爰视频| 国产成人aa在线观看| 91成人精品电影| 精品酒店卫生间| 免费播放大片免费观看视频在线观看| 久久97久久精品| 久久久久久人妻| 亚洲成色77777| 精品一区二区免费观看| 性色avwww在线观看| 精品人妻一区二区三区麻豆| 内射极品少妇av片p| 国产精品福利在线免费观看| 免费观看性生交大片5| 亚洲成人一二三区av| 人人妻人人爽人人添夜夜欢视频 | 欧美性感艳星| 97在线人人人人妻| 我的老师免费观看完整版| 午夜免费鲁丝| 久久久a久久爽久久v久久| 国产黄色免费在线视频| 久久99热这里只频精品6学生| 国产乱人偷精品视频| 交换朋友夫妻互换小说| 亚洲国产精品一区三区| 国产乱人偷精品视频| 日本-黄色视频高清免费观看| 亚洲av电影在线观看一区二区三区| 老熟女久久久| 国产熟女欧美一区二区| 亚洲av.av天堂| 一级黄片播放器| 精品酒店卫生间| 一级片'在线观看视频| 桃花免费在线播放| 97精品久久久久久久久久精品| 国产在线免费精品| 日韩大片免费观看网站| 亚洲国产最新在线播放| 日韩视频在线欧美| 国产成人精品无人区| 精品一区在线观看国产| 欧美3d第一页| av有码第一页| av又黄又爽大尺度在线免费看| 在线免费观看不下载黄p国产| 看免费成人av毛片| 丰满人妻一区二区三区视频av| 97精品久久久久久久久久精品| 岛国毛片在线播放| 欧美精品国产亚洲| 男人狂女人下面高潮的视频| 夜夜看夜夜爽夜夜摸| 人人澡人人妻人| 国产亚洲av片在线观看秒播厂| 久久热精品热| 高清不卡的av网站| 中文字幕久久专区| 日韩中字成人| a 毛片基地| 国产精品一区二区性色av| 国内精品宾馆在线| 欧美日韩一区二区视频在线观看视频在线| 国模一区二区三区四区视频| 国产精品久久久久久久电影| 免费观看的影片在线观看| 在线天堂最新版资源| 精品久久久精品久久久| .国产精品久久| 卡戴珊不雅视频在线播放| 欧美激情极品国产一区二区三区 | 一边亲一边摸免费视频| 少妇被粗大的猛进出69影院 | 亚洲av国产av综合av卡| 建设人人有责人人尽责人人享有的| 男人舔奶头视频| 少妇人妻精品综合一区二区| av在线观看视频网站免费| 一级黄片播放器| freevideosex欧美| 国产无遮挡羞羞视频在线观看| 欧美亚洲 丝袜 人妻 在线| 一级黄片播放器| 汤姆久久久久久久影院中文字幕| 免费黄色在线免费观看| 最新的欧美精品一区二区| 韩国av在线不卡| 国产亚洲5aaaaa淫片| 建设人人有责人人尽责人人享有的| 国产极品天堂在线| 国产无遮挡羞羞视频在线观看| 亚洲美女黄色视频免费看| 日本黄色片子视频| 久久精品国产自在天天线| 亚洲国产成人一精品久久久| 99久久精品热视频| 国产一区有黄有色的免费视频| 国产欧美另类精品又又久久亚洲欧美| 日韩精品免费视频一区二区三区 | 亚洲精品视频女| 97精品久久久久久久久久精品| 男女边摸边吃奶| 51国产日韩欧美| 亚洲一区二区三区欧美精品| 超碰97精品在线观看| 哪个播放器可以免费观看大片| 亚洲国产成人一精品久久久| 99热网站在线观看| 免费少妇av软件| 高清av免费在线| 成年女人在线观看亚洲视频| 伦理电影免费视频| 高清视频免费观看一区二区| 日日摸夜夜添夜夜爱| 亚洲成色77777| 精品久久久精品久久久| 两个人的视频大全免费| 99热6这里只有精品| 国产黄频视频在线观看| 欧美成人精品欧美一级黄| 精品国产乱码久久久久久小说| 高清毛片免费看| 一级毛片 在线播放| 国产日韩欧美在线精品| 激情五月婷婷亚洲| 99re6热这里在线精品视频| av国产精品久久久久影院| 国产精品.久久久| 成人无遮挡网站| 中文字幕人妻丝袜制服| 日韩 亚洲 欧美在线| 国产精品国产三级专区第一集| 韩国高清视频一区二区三区| 久久久久精品久久久久真实原创| 精品亚洲乱码少妇综合久久| 中文精品一卡2卡3卡4更新| 午夜免费观看性视频| 久久影院123| 亚洲,欧美,日韩| 日日爽夜夜爽网站| 色网站视频免费| 在现免费观看毛片| 亚洲精品亚洲一区二区| 欧美 日韩 精品 国产| 99久久精品热视频| 大又大粗又爽又黄少妇毛片口| 免费av中文字幕在线| 国产精品不卡视频一区二区| 一二三四中文在线观看免费高清| 插逼视频在线观看| 水蜜桃什么品种好| 亚洲天堂av无毛| 高清毛片免费看| 久久久久视频综合| 97超碰精品成人国产| 最近中文字幕高清免费大全6| 亚洲经典国产精华液单| 国产片特级美女逼逼视频| 婷婷色麻豆天堂久久| 国产精品久久久久久久久免| 我的女老师完整版在线观看| 日韩一区二区视频免费看| 国产日韩一区二区三区精品不卡 | 男女免费视频国产| 麻豆精品久久久久久蜜桃| 日韩中字成人| 国产一区有黄有色的免费视频| 日本欧美国产在线视频| 成年av动漫网址| 日本午夜av视频| 啦啦啦啦在线视频资源| 亚洲精品亚洲一区二区| 亚洲精品成人av观看孕妇| 99视频精品全部免费 在线| 卡戴珊不雅视频在线播放| 女性生殖器流出的白浆| 亚洲色图综合在线观看| 黄色配什么色好看| 国精品久久久久久国模美| 欧美日本中文国产一区发布| 观看美女的网站| 97超碰精品成人国产| 国产高清国产精品国产三级| 男人添女人高潮全过程视频| 国产精品国产三级国产av玫瑰| 观看免费一级毛片| 日韩三级伦理在线观看| 国产高清三级在线| 亚洲自偷自拍三级| 精品99又大又爽又粗少妇毛片| 又大又黄又爽视频免费| 啦啦啦啦在线视频资源| 只有这里有精品99| 人妻少妇偷人精品九色| 国产成人精品久久久久久| 夜夜骑夜夜射夜夜干| 久久久久久久久久人人人人人人| 久久97久久精品| 97超碰精品成人国产| 中国国产av一级| 欧美激情国产日韩精品一区| 免费大片18禁| 婷婷色综合www| 看十八女毛片水多多多| 免费观看性生交大片5| 国产一区二区三区综合在线观看 | 在线天堂最新版资源| 国产亚洲精品久久久com| 久久人妻熟女aⅴ| 国产精品久久久久久av不卡| 久久久国产精品麻豆| 国产精品久久久久久av不卡| 女人久久www免费人成看片| 久久精品久久精品一区二区三区| 国产淫语在线视频| 男女国产视频网站| 91精品一卡2卡3卡4卡| 777米奇影视久久| 在线观看一区二区三区激情| 我的女老师完整版在线观看| 久久久国产欧美日韩av| 久久久久久久久大av| 精品亚洲成a人片在线观看| 十分钟在线观看高清视频www | 免费观看在线日韩| 日韩亚洲欧美综合| 九九在线视频观看精品| 六月丁香七月| 久久女婷五月综合色啪小说| 好男人视频免费观看在线| 国产高清不卡午夜福利| 99久久综合免费| 国产在线视频一区二区| 国产精品女同一区二区软件| 人妻一区二区av| 嫩草影院新地址| 成年人免费黄色播放视频 | 老熟女久久久| 一级毛片 在线播放| 日本91视频免费播放| 丰满人妻一区二区三区视频av| 老熟女久久久| 少妇人妻久久综合中文| 日日啪夜夜爽| 亚洲精品色激情综合| 少妇精品久久久久久久| 日韩制服骚丝袜av| 性色av一级| 日本av免费视频播放| 亚洲av综合色区一区| 精品久久久久久久久亚洲| 熟妇人妻不卡中文字幕| 最新中文字幕久久久久| 三上悠亚av全集在线观看 | 亚洲成人一二三区av| 只有这里有精品99| 国产视频内射| 精品人妻偷拍中文字幕| 乱系列少妇在线播放| av在线观看视频网站免费| 久久99热6这里只有精品| 一区二区三区精品91| 日韩 亚洲 欧美在线| 国产免费一级a男人的天堂| 免费观看的影片在线观看| 看十八女毛片水多多多| 天美传媒精品一区二区| 亚洲,欧美,日韩| 日日摸夜夜添夜夜爱| 国产男女超爽视频在线观看| 亚洲婷婷狠狠爱综合网| 成人免费观看视频高清| 91在线精品国自产拍蜜月| 欧美日本中文国产一区发布| 久久99热这里只频精品6学生| 国产伦理片在线播放av一区| 日韩强制内射视频| 97超碰精品成人国产| 黑人猛操日本美女一级片| 在线观看免费视频网站a站| 久久人人爽人人爽人人片va| 免费人成在线观看视频色| 免费久久久久久久精品成人欧美视频 | 人妻一区二区av| 高清不卡的av网站| av福利片在线观看| 国产老妇伦熟女老妇高清| 久久免费观看电影| 另类精品久久| 26uuu在线亚洲综合色| 新久久久久国产一级毛片| 高清av免费在线| 少妇被粗大的猛进出69影院 | 香蕉精品网在线| 久久精品国产自在天天线| av在线app专区| 午夜激情久久久久久久| 少妇的逼水好多| 久久久久久久精品精品| 亚洲经典国产精华液单| 婷婷色综合www| 在线观看av片永久免费下载| 99热全是精品| 久久精品国产自在天天线| 国产色爽女视频免费观看| 夫妻午夜视频| 亚洲精品一二三| 国产精品久久久久久av不卡| 久久99一区二区三区| 熟女人妻精品中文字幕| 两个人的视频大全免费| 亚洲精品日本国产第一区| 国产永久视频网站| 国产成人aa在线观看| 天天躁夜夜躁狠狠久久av| 三级国产精品片| 成人黄色视频免费在线看| 免费黄网站久久成人精品| 丰满人妻一区二区三区视频av| 日韩精品免费视频一区二区三区 | 香蕉精品网在线| 九九久久精品国产亚洲av麻豆| 又爽又黄a免费视频| 日韩中文字幕视频在线看片| 黄色欧美视频在线观看| 少妇的逼好多水| 赤兔流量卡办理| 9色porny在线观看| 欧美激情极品国产一区二区三区 | 欧美日韩综合久久久久久| 日韩亚洲欧美综合|