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

    基于WIM的大跨鋼箱梁懸索橋疲勞荷載譜研究

    2022-12-03 10:18:14賀偉誠(chéng)
    關(guān)鍵詞:車(chē)重軸重概率密度

    賀偉誠(chéng),王 超

    (湖北工業(yè)大學(xué)土木建筑與環(huán)境學(xué)院,湖北 武漢 430068)

    大跨度橋梁主梁廣泛采用正交異性鋼橋面板鋼箱梁結(jié)構(gòu),這種結(jié)構(gòu)具有自重輕、抗扭及橫向抗彎剛度大、整體性強(qiáng)、經(jīng)濟(jì)性好等優(yōu)點(diǎn)[1]。然而,由于車(chē)輪荷載的反復(fù)作用,加之焊接殘余應(yīng)力及初始缺陷的影響,橋梁極易產(chǎn)生疲勞損傷累積,嚴(yán)重時(shí)甚至造成結(jié)構(gòu)破壞。針對(duì)鋼箱梁疲勞損傷問(wèn)題,研究建立疲勞車(chē)輛模型并分析各車(chē)型對(duì)橋梁疲勞損傷貢獻(xiàn)率是重要的方向之一,國(guó)內(nèi)外學(xué)者在該領(lǐng)域開(kāi)展了相關(guān)的研究工作。宗周紅等[2]研究了京滬高速公路的重車(chē)疲勞荷載譜,發(fā)現(xiàn)重車(chē)占比較高,且六軸二類(lèi)車(chē)型為等效標(biāo)準(zhǔn)疲勞車(chē)。胡森林等[3]提出了針對(duì)廣西地區(qū)疲勞設(shè)計(jì)的疲勞損傷修正系數(shù)建議值。王元清等[4]針對(duì)建立的車(chē)輛荷載模型進(jìn)行分析,得出42 d樣本量可以用來(lái)代表整體分布特征的結(jié)論。鄧揚(yáng)等[5]、Zhang等[6]研究了分時(shí)段的疲勞荷載譜,前者發(fā)現(xiàn)晚上超車(chē)道的車(chē)重均值比白天高,后者發(fā)現(xiàn)汀九橋的疲勞壽命在120 a的設(shè)計(jì)限值內(nèi)。Yan等[7]研究了中、小跨徑橋梁的疲勞應(yīng)力譜,發(fā)現(xiàn)超載控制對(duì)于疲勞應(yīng)力譜中高幅應(yīng)力的概率密度有很大的影響。翟慕賽等[8]建立了不同等級(jí)公路的疲勞荷載譜,發(fā)現(xiàn)不同等級(jí)公路其疲勞荷載模型也不同。孫秀雅等[9]基于疲勞荷載譜針對(duì)壩陵河大橋疲勞易損部位進(jìn)行了疲勞壽命分析,得出橋上疲勞細(xì)節(jié)的疲勞壽命高于設(shè)計(jì)值的結(jié)論。陳斌等[10]研究了目標(biāo)基準(zhǔn)期內(nèi)車(chē)輛荷載的極值分布,發(fā)現(xiàn)車(chē)輛荷載存在離散性大的特點(diǎn)。羅媛等[11]提出了橋梁上基于隨機(jī)重載車(chē)流的疲勞可靠度分析方法,發(fā)現(xiàn)車(chē)輛軸數(shù)多少和橋梁等效疲勞應(yīng)力幅的高低存在負(fù)相關(guān)性。劉揚(yáng)等[12]研究了關(guān)于隨機(jī)車(chē)流作用下簡(jiǎn)支梁橋設(shè)計(jì)基準(zhǔn)期內(nèi)荷載效應(yīng)極值概率分布的問(wèn)題。Deng等[13]建立了疲勞荷載譜并單獨(dú)討論了六軸車(chē)的車(chē)輛模型,提出了一種針對(duì)疲勞車(chē)輛模型研究的新方法。MulyadiEndang等[14]進(jìn)行了對(duì)比S-N曲線(xiàn)與標(biāo)稱(chēng)疲勞抗力的研究,得出S-N曲線(xiàn)上50 a壽命依然低于疲勞極限的結(jié)論。Guo等[15]利用車(chē)輛疲勞荷載譜對(duì)橋梁的可靠度進(jìn)行了分析,發(fā)現(xiàn)隨著時(shí)間的推移和車(chē)流量增加,可靠度會(huì)加速降低。

    我國(guó)幅員遼闊,各地環(huán)境不同、橋型各異、交通流量狀況復(fù)雜多變,對(duì)于不同地區(qū)橋梁特別是特大跨度鋼箱梁橋的疲勞問(wèn)題,有必要針對(duì)其實(shí)際運(yùn)營(yíng)環(huán)境進(jìn)行專(zhuān)門(mén)的研究分析[16]。另外,橋梁的實(shí)際交通狀況會(huì)受工作日、節(jié)假日、季節(jié)變化影響發(fā)生改變,為了確保統(tǒng)計(jì)得到的疲勞荷載譜更加準(zhǔn)確和具有代表性,有必要采用較長(zhǎng)時(shí)間段的數(shù)據(jù)進(jìn)行分析。因此,本文將針對(duì)某大跨度正交異性橋面板鋼箱梁懸索橋進(jìn)行疲勞分析,利用動(dòng)態(tài)稱(chēng)重系統(tǒng)(WIM)長(zhǎng)期采集的實(shí)橋車(chē)輛荷載數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析,引入確定模型最佳分量數(shù)的兩種準(zhǔn)則,采用混合高斯分布擬合車(chē)重?cái)?shù)據(jù)的分布規(guī)律,得到車(chē)重概率密度曲線(xiàn)并建立疲勞荷載譜,進(jìn)而計(jì)算出各車(chē)型對(duì)該橋的疲勞損傷貢獻(xiàn)率。

    1 車(chē)輛荷載數(shù)據(jù)統(tǒng)計(jì)

    1.1 工程概況

    某大跨度鋼箱梁懸索橋是繞城公路的重要組成部分,設(shè)計(jì)速度為120 km/h,鋼箱梁主梁全寬38.5 m,設(shè)雙向六車(chē)道,主梁橫斷面布置圖如圖1所示。為了保障大橋的安全運(yùn)營(yíng),大橋建設(shè)了橋梁健康監(jiān)測(cè)系統(tǒng),實(shí)時(shí)采集包括車(chē)輛信息在內(nèi)的各類(lèi)橋梁運(yùn)營(yíng)期的關(guān)鍵參數(shù)。由于此橋地處城市外環(huán),是連接市內(nèi)外的樞紐,交通量大,車(chē)型多樣,交通狀況復(fù)雜,不同車(chē)型會(huì)對(duì)橋梁造成不同程度的疲勞損傷,因此,這里將對(duì)該橋WIM系統(tǒng)采集的長(zhǎng)期車(chē)輛荷載數(shù)據(jù)進(jìn)行分析研究。

    圖1 某鋼箱梁懸索橋橫斷面圖 mm

    1.2 車(chē)型分類(lèi)

    由于實(shí)際車(chē)輛類(lèi)型龐雜,沒(méi)有必要對(duì)每一種車(chē)型進(jìn)行統(tǒng)計(jì)分析,為了分析簡(jiǎn)易性和有效性,需要對(duì)其進(jìn)行分類(lèi)。由于結(jié)構(gòu)疲勞損傷主要由應(yīng)力幅和循環(huán)次數(shù)控制,而正交異性橋面板鋼箱梁應(yīng)力幅循環(huán)次數(shù)主要受軸數(shù)影響,這里依據(jù)軸數(shù)將車(chē)輛進(jìn)行分類(lèi),特別是實(shí)際交通流中大量存在的各種兩軸車(chē)都?xì)w為一類(lèi)。另外,車(chē)輛荷載與橋梁疲勞壽命成負(fù)相關(guān)性,過(guò)小的車(chē)輛荷載對(duì)疲勞壽命的影響可忽略不計(jì),為了簡(jiǎn)化分析,可以先對(duì)數(shù)據(jù)進(jìn)行預(yù)處理,剔除車(chē)重過(guò)小(小于3 t)的車(chē)輛[17],再進(jìn)行統(tǒng)計(jì)分析,預(yù)處理后得到的各典型車(chē)型交通量如表1所示。

    表1 各典型車(chē)型交通量

    從統(tǒng)計(jì)結(jié)果可以看出該橋車(chē)流量以V2型車(chē)即兩軸車(chē)為主,占到了44.5%,六軸大貨車(chē)僅次于兩軸車(chē),相對(duì)占比也比較高,達(dá)到33.14%,其他四種車(chē)型占比均低于10%。考慮到相對(duì)于其他大貨車(chē)來(lái)說(shuō),兩軸車(chē)載重和軸數(shù)相對(duì)較小,多軸大貨車(chē)的疲勞影響也不容小視,下文將對(duì)此進(jìn)行詳細(xì)分析。

    1.3 交通量統(tǒng)計(jì)

    基于表1的典型車(chē)輛分類(lèi),現(xiàn)統(tǒng)計(jì)出每周各車(chē)型的車(chē)流量情況,由于少部分?jǐn)?shù)據(jù)缺失,故共分析了24周的車(chē)輛數(shù)據(jù),統(tǒng)計(jì)的各車(chē)型車(chē)流量變化如圖2所示。從圖2a可知各周的車(chē)流量變化都比較明顯,其中第22周的變化尤為顯著,推測(cè)可能是橋面檢測(cè)或維修等原因?qū)е拢粡膱D2b可以直觀看出各周車(chē)流量的構(gòu)成也存在不同程度的變化,兩軸車(chē)和六軸車(chē)占比變化更為突出,并且兩類(lèi)車(chē)占比趨勢(shì)基本相反。

    圖2 交通量變化

    由于不同車(chē)型及車(chē)重的車(chē)流量每周變化明顯,需要采用較長(zhǎng)時(shí)間的監(jiān)測(cè)數(shù)據(jù)統(tǒng)計(jì)分析車(chē)流荷載分布特征,才能減少因偶然性而導(dǎo)致數(shù)據(jù)出現(xiàn)過(guò)大或過(guò)小的情況發(fā)生,從而確保建立的疲勞荷載譜更加準(zhǔn)確并且具有代表性。

    2 模型車(chē)輛疲勞荷載譜

    2.1 車(chē)重統(tǒng)計(jì)

    基于表1的交通量統(tǒng)計(jì),再依據(jù)各類(lèi)型車(chē)車(chē)重?cái)?shù)據(jù)即可統(tǒng)計(jì)得到各車(chē)型的車(chē)重分布特征。以V2車(chē)型為例,根據(jù)實(shí)測(cè)數(shù)據(jù)統(tǒng)計(jì)其車(chē)重分布如圖3所示。這里將根據(jù)車(chē)重分布圖擬合其概率密度曲線(xiàn),由圖3可以看出,車(chē)重分布圖呈現(xiàn)出多峰的形式,普通正態(tài)分布(圖中虛線(xiàn))很難較好地?cái)M合其分布特征,擬采用混合高斯分布擬合該分布規(guī)律,其概率密度函數(shù)如式(1):

    圖3 二軸車(chē)的車(chē)重-概率密度及擬合曲線(xiàn)

    (1)

    式中:M為分量數(shù),描述分布峰值個(gè)數(shù);Pi為各分布型對(duì)應(yīng)的系數(shù);fi(x)為各高斯分布的概率密度函數(shù)。

    采用混合高斯分布擬合需要確定擬合曲線(xiàn)的高斯分量數(shù),當(dāng)分量數(shù)過(guò)小時(shí),其擬合效果較差,分量數(shù)過(guò)多又會(huì)使問(wèn)題變復(fù)雜,對(duì)于解決實(shí)際工程問(wèn)題收效不大。為了更加準(zhǔn)確描述概率分布情況,這里引入赤池信息準(zhǔn)則(AIC)和貝葉斯信息準(zhǔn)則(BIC)來(lái)優(yōu)化確定最佳分量數(shù),AIC和BIC分別定義為式(2)和(3):

    AIC=-2lnL+2M

    (2)

    BIC=-2lnL+MlnN

    (3)

    (4)

    式中:M為分量數(shù);lnL為模型的對(duì)數(shù)極大似然函數(shù)值;N為模型中數(shù)據(jù)的個(gè)數(shù);SSE為殘差平方和,具體見(jiàn)式(5)。

    模型最佳分量數(shù)的選取應(yīng)優(yōu)先考慮AIC值,取AIC最小時(shí)對(duì)應(yīng)的M值;而當(dāng)AIC值相差不大時(shí),則引入BIC,最優(yōu)分量數(shù)依舊取其值最小時(shí)對(duì)應(yīng)的M值。

    對(duì)于概率密度函數(shù)擬合的效果,采用擬合函數(shù)的SSE(和方差)以及R2(確定系數(shù))進(jìn)行評(píng)估,其定義為:

    (5)

    (6)

    擬合曲線(xiàn)的和方差值越接近0,確定系數(shù)越接近1,則擬合效果越好。

    針對(duì)圖3所示的V2型車(chē),計(jì)算出采用不同分量數(shù)的混合高斯分布擬合車(chē)重概率分布時(shí)對(duì)應(yīng)的AIC和BIC值,如圖4所示。

    圖4 V2型車(chē)混合分布模型信息準(zhǔn)則值變化曲線(xiàn)

    由圖可見(jiàn),當(dāng)分量數(shù)為3時(shí),AIC和BIC的值均達(dá)到最小,故V2型車(chē)多峰分布模型的分量數(shù)M取3。相應(yīng)的混合高斯分布擬合曲線(xiàn)如圖3中紅線(xiàn)所示。

    為了檢驗(yàn)實(shí)測(cè)值的統(tǒng)計(jì)規(guī)律和擬合曲線(xiàn)之間的契合度,可求解出擬合函數(shù)的SSE(和方差)以及R2(確定系數(shù)),其擬合優(yōu)度如表2所示。

    從統(tǒng)計(jì)結(jié)果可以看出,V2型車(chē)數(shù)量最多,車(chē)型復(fù)雜,但整體而言車(chē)重量都不大,基本在6~13 t之間。采用同樣的方法,對(duì)其他四種車(chē)型的車(chē)重進(jìn)行統(tǒng)計(jì),具體統(tǒng)計(jì)結(jié)果如圖5所示,可以發(fā)現(xiàn)三軸車(chē)和四軸車(chē)的車(chē)重分別為10.64 t和14.55 t時(shí)概率密度均達(dá)到最大峰值,且三軸車(chē)和四軸車(chē)分布規(guī)律類(lèi)似,均出現(xiàn)峰與峰之間距離過(guò)大的情況,最右的一個(gè)峰其峰值較低,此峰值表明存在部分車(chē)輛滿(mǎn)載或超載;五軸車(chē)和六軸車(chē)均呈現(xiàn)為雙峰,兩類(lèi)車(chē)型的車(chē)重分別為16.39 t和17.98 t時(shí)概率密度達(dá)到最大峰值。而用于擬合各車(chē)型車(chē)重概率分布的分量數(shù)、擬合曲線(xiàn)的參數(shù)以及擬合優(yōu)度等參數(shù)可見(jiàn)表2。

    圖5 各車(chē)型車(chē)重概率密度模型

    表2 各類(lèi)車(chē)型車(chē)重概率密度曲線(xiàn)擬合參數(shù)及擬合優(yōu)度

    由表2各車(chē)型概率密度擬合曲線(xiàn)的和方差以及確定系數(shù)可知,五種車(chē)型混合高斯分布的和方差均接近于0,確定系數(shù)均大于0.96,由此可以確定混合高斯分布能較好擬合各車(chē)型的概率密度曲線(xiàn)。

    2.2 建立疲勞荷載譜

    依據(jù)Miner線(xiàn)性積傷率和等效疲勞損傷原理,當(dāng)累計(jì)疲勞損傷相同時(shí),可以將i個(gè)不同車(chē)重作用1次視為1個(gè)等效車(chē)重作用i次,等效車(chē)重計(jì)算式如下:

    (7)

    式中:W(eq)k為第k種車(chē)型的等效車(chē)重;pi為第k種車(chē)型中i車(chē)出現(xiàn)的頻率;Wi為該車(chē)車(chē)重。

    將原始數(shù)據(jù)中每種車(chē)型的各軸重分別統(tǒng)計(jì),計(jì)算各軸重在車(chē)重的占比,即可得到各車(chē)型軸重分配系數(shù)γk,計(jì)算方法如式(8)所示,表3為各車(chē)型軸重分配系數(shù)。

    表3 各車(chē)型軸重分配系數(shù)

    (8)

    式中:γk為第k種車(chē)型的軸重分配系數(shù);ω(eq)k,j為第k種車(chē)型第j軸的等效軸重;W(eq)k為第k種車(chē)型的等效車(chē)重。

    運(yùn)用式(7)計(jì)算出等效車(chē)重后,再結(jié)合軸重分配系數(shù)即可得到各車(chē)型的等效軸重?cái)?shù)據(jù),然后建立車(chē)輛的疲勞荷載譜,見(jiàn)表4。

    表4 疲勞荷載譜

    3 疲勞損傷貢獻(xiàn)率

    橋梁結(jié)構(gòu)的疲勞損傷性能可以由不同車(chē)輛對(duì)橋梁產(chǎn)生的疲勞損傷貢獻(xiàn)率來(lái)反映,各車(chē)型對(duì)結(jié)構(gòu)產(chǎn)生的損傷值和疲勞損傷貢獻(xiàn)率根據(jù)上文統(tǒng)計(jì)的車(chē)輛疲勞荷載譜來(lái)計(jì)算。對(duì)于正交異性橋面板鋼箱梁,橋面板鋼結(jié)構(gòu)在車(chē)輛局部荷載作用下應(yīng)力影響線(xiàn)很短[18],局部效應(yīng)比較明顯,所以結(jié)構(gòu)疲勞損傷與軸重直接相關(guān),可以通過(guò)計(jì)算等效軸重的疲勞貢獻(xiàn)率再結(jié)合軸數(shù)來(lái)得到等效車(chē)重的疲勞貢獻(xiàn)率。具體計(jì)算方法如下:

    第k種車(chē)型第j軸等效軸重造成的橋梁損傷值D(eq)k,j為:

    (9)

    式中,c為結(jié)構(gòu)或材料的性能參數(shù);β為軸重引起的應(yīng)力幅值系數(shù);m為S-N曲線(xiàn)的斜率負(fù)倒數(shù),取m=3;nk,j為第k種車(chē)型第j軸的總軸數(shù);ω(eq)k,j為第k種車(chē)型第j軸的等效軸重。

    各軸產(chǎn)生的總疲勞損傷值D為:

    (10)

    式中,kmax為車(chē)型總數(shù);jmax為第k種車(chē)型的最大軸數(shù)。

    第k種車(chē)型第j軸等效軸重造成的疲勞損傷貢獻(xiàn)率為:

    (11)

    第k種車(chē)型的疲勞損傷貢獻(xiàn)率λk定義為:

    (12)

    由此可以計(jì)算得到各軸及各車(chē)型對(duì)橋梁造成的疲勞損傷貢獻(xiàn)率,如表5所示。

    表5 各軸及各車(chē)型的疲勞貢獻(xiàn)率

    從表5數(shù)據(jù)可以看出:1)兩軸車(chē)的頻率在各類(lèi)型車(chē)中最大,但六軸車(chē)等效車(chē)重和各軸等效軸重產(chǎn)生的疲勞貢獻(xiàn)率最高,原因是在疲勞損傷貢獻(xiàn)率的計(jì)算中等效軸重的權(quán)重高于車(chē)輛數(shù),并且六軸車(chē)的軸數(shù)最多、該車(chē)型各軸等效軸重相比于其它車(chē)型也較高;2)兩軸車(chē)頻率遠(yuǎn)遠(yuǎn)高于三軸車(chē),但等效車(chē)重產(chǎn)生的疲勞損傷貢獻(xiàn)率卻略低于三軸車(chē),原因是雖然這兩類(lèi)車(chē)型的軸1和軸2等效軸重相差不大,但三軸車(chē)的第3軸的等效軸重遠(yuǎn)高于兩軸車(chē)的所有軸重;3)四軸車(chē)和五軸車(chē)的頻率本身較低,并且各軸等效軸重都不高,最終產(chǎn)生的疲勞損傷貢獻(xiàn)率都較低。

    4 結(jié)論

    通過(guò)對(duì)某正交異性橋面板鋼箱梁橋動(dòng)態(tài)稱(chēng)重系統(tǒng)長(zhǎng)期監(jiān)測(cè)的車(chē)輛荷載數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析,采用混合高斯模型擬合車(chē)重荷載分布概率密度,建立了該橋的疲勞車(chē)輛荷載譜,最后分析研究了各車(chē)型對(duì)該橋的疲勞損傷貢獻(xiàn)率,分析結(jié)果表明:

    1)基于長(zhǎng)時(shí)間的數(shù)據(jù)采集和分析,發(fā)現(xiàn)各周的車(chē)流量變化顯著,每周車(chē)流量的車(chē)型構(gòu)成也有很大變化,其中兩軸車(chē)和六軸車(chē)占比變化最明顯,且這兩類(lèi)型車(chē)的占比趨勢(shì)基本相反。長(zhǎng)期的統(tǒng)計(jì)數(shù)據(jù)可以減小偶然因素帶來(lái)的影響,從而使研究更具有代表性和普遍性,所以有必要對(duì)長(zhǎng)期車(chē)流情況進(jìn)行研究。

    2)各車(chē)型的車(chē)重概率密度分布均服從混合高斯分布,且由赤池信息準(zhǔn)則和貝葉斯信息準(zhǔn)則確定高斯分量數(shù)后的擬合曲線(xiàn)其和方差均接近于0,確定系數(shù)均大于0.96,可知擬合效果較好。

    3)通過(guò)建立各車(chē)型疲勞荷載譜可知,V3型車(chē)第3軸的等效軸重與各類(lèi)車(chē)型的等效軸重相比都最大,除此之外的各車(chē)型等效軸重相差都不明顯,但是隨著軸數(shù)的增加,等效車(chē)重也隨之增大,所以V2型車(chē)等效車(chē)重最小,V6型車(chē)最大。

    4)六軸車(chē)各等效軸重產(chǎn)生的總疲勞損傷貢獻(xiàn)率最高,三軸車(chē)其次,兩軸車(chē)僅次于三軸車(chē),四軸車(chē)和五軸車(chē)最低。原因是六軸車(chē)頻率較高,各軸的等效軸重也較大,且軸數(shù)多,再加上等效軸重相比于車(chē)流量更加影響疲勞損傷貢獻(xiàn)率的計(jì)算,導(dǎo)致疲勞損傷貢獻(xiàn)率最高和最不利車(chē)型都是六軸車(chē)。

    猜你喜歡
    車(chē)重軸重概率密度
    基于三軸加速度計(jì)最小二乘法質(zhì)量和坡度估計(jì)
    基于發(fā)動(dòng)機(jī)輔助制動(dòng)功能的智能車(chē)重預(yù)估方法
    連續(xù)型隨機(jī)變量函數(shù)的概率密度公式
    考慮橫豎向車(chē)橋耦合曲線(xiàn)橋的動(dòng)力響應(yīng)影響
    20t軸重米軌轉(zhuǎn)向架裝用搖枕的研制
    基于車(chē)橋耦合振動(dòng)的鋼管混凝土系桿拱橋動(dòng)力響應(yīng)分析
    Hunt過(guò)程在Girsanov變換下的轉(zhuǎn)移概率密度的表示公式
    32.5t軸重貨車(chē)作用下重載鐵路軌道的合理剛度
    隨機(jī)變量線(xiàn)性組合的分布的一個(gè)算法
    隨機(jī)結(jié)構(gòu)-TMD優(yōu)化設(shè)計(jì)與概率密度演化研究
    热99re8久久精品国产| 久久久精品免费免费高清| 成人永久免费在线观看视频 | 50天的宝宝边吃奶边哭怎么回事| 国产熟女午夜一区二区三区| 亚洲第一av免费看| 日本黄色视频三级网站网址 | 午夜激情久久久久久久| 一区在线观看完整版| 国产精品一区二区在线观看99| 久久久精品免费免费高清| 三上悠亚av全集在线观看| 精品卡一卡二卡四卡免费| 丰满迷人的少妇在线观看| e午夜精品久久久久久久| 国产成人影院久久av| videosex国产| 亚洲久久久国产精品| 亚洲国产毛片av蜜桃av| 悠悠久久av| 久久久精品区二区三区| 老司机深夜福利视频在线观看| 亚洲欧美日韩高清在线视频 | 12—13女人毛片做爰片一| 岛国在线观看网站| 99久久99久久久精品蜜桃| 久久久久久久久免费视频了| av网站在线播放免费| e午夜精品久久久久久久| 免费观看a级毛片全部| 久久久久久久久久久久大奶| 国产av精品麻豆| 成人亚洲精品一区在线观看| 国产在线视频一区二区| 日韩精品免费视频一区二区三区| 久久人人97超碰香蕉20202| 另类精品久久| 亚洲精品一卡2卡三卡4卡5卡| 欧美人与性动交α欧美软件| 看免费av毛片| av国产精品久久久久影院| 91大片在线观看| 亚洲国产欧美在线一区| 极品人妻少妇av视频| 日本av免费视频播放| 电影成人av| 精品免费久久久久久久清纯 | 国产精品二区激情视频| avwww免费| 新久久久久国产一级毛片| 男人舔女人的私密视频| 天堂俺去俺来也www色官网| 亚洲一区中文字幕在线| 久久精品亚洲精品国产色婷小说| 啦啦啦在线免费观看视频4| 两个人看的免费小视频| 国产精品久久电影中文字幕 | 19禁男女啪啪无遮挡网站| 午夜福利欧美成人| 亚洲 国产 在线| 免费在线观看影片大全网站| 亚洲成人免费电影在线观看| 欧美国产精品va在线观看不卡| 中文字幕高清在线视频| 一区二区三区国产精品乱码| 91老司机精品| 久久国产精品人妻蜜桃| 成人三级做爰电影| 桃花免费在线播放| 精品国产一区二区久久| 精品福利永久在线观看| 亚洲精品中文字幕一二三四区 | 久久国产亚洲av麻豆专区| 精品少妇久久久久久888优播| 久久午夜综合久久蜜桃| 欧美国产精品一级二级三级| 法律面前人人平等表现在哪些方面| e午夜精品久久久久久久| 最新在线观看一区二区三区| 人人妻人人澡人人看| 欧美亚洲 丝袜 人妻 在线| 亚洲av欧美aⅴ国产| 精品福利永久在线观看| 国产成人精品无人区| 久久中文字幕一级| 国产精品一区二区免费欧美| 久久精品91无色码中文字幕| a级毛片在线看网站| 午夜福利在线观看吧| 欧美精品人与动牲交sv欧美| 免费观看av网站的网址| 精品久久久久久久毛片微露脸| 成人免费观看视频高清| 国产精品自产拍在线观看55亚洲 | 午夜福利免费观看在线| 十八禁网站免费在线| 亚洲成av片中文字幕在线观看| 日韩欧美一区二区三区在线观看 | 成人永久免费在线观看视频 | 中文字幕人妻丝袜一区二区| 欧美激情 高清一区二区三区| 蜜桃在线观看..| 亚洲男人天堂网一区| 国产一区二区三区综合在线观看| 女人久久www免费人成看片| 欧美激情极品国产一区二区三区| 男女床上黄色一级片免费看| 欧美午夜高清在线| 18禁国产床啪视频网站| 国产午夜精品久久久久久| 欧美激情极品国产一区二区三区| 日韩三级视频一区二区三区| 美女高潮到喷水免费观看| 俄罗斯特黄特色一大片| 肉色欧美久久久久久久蜜桃| 中文字幕另类日韩欧美亚洲嫩草| 久久中文字幕一级| 久久久久久久精品吃奶| 最近最新免费中文字幕在线| 麻豆av在线久日| 日韩大码丰满熟妇| 国产av又大| 欧美日韩一级在线毛片| 日本wwww免费看| 国产成人免费观看mmmm| 国产成人一区二区三区免费视频网站| 亚洲人成电影免费在线| 18禁黄网站禁片午夜丰满| 18禁黄网站禁片午夜丰满| 亚洲成a人片在线一区二区| 国产成人av教育| 最新在线观看一区二区三区| 久久久久久亚洲精品国产蜜桃av| 老司机午夜十八禁免费视频| 免费不卡黄色视频| 一本—道久久a久久精品蜜桃钙片| 久久久国产一区二区| av一本久久久久| 久久热在线av| 窝窝影院91人妻| 91老司机精品| 久久精品国产综合久久久| 两性午夜刺激爽爽歪歪视频在线观看 | 电影成人av| 亚洲av欧美aⅴ国产| 精品福利永久在线观看| 天天操日日干夜夜撸| 视频在线观看一区二区三区| 如日韩欧美国产精品一区二区三区| 国产国语露脸激情在线看| 精品国产一区二区久久| 中文字幕色久视频| 国产高清videossex| 一级毛片精品| 一本综合久久免费| 麻豆乱淫一区二区| 亚洲熟女毛片儿| 精品国产超薄肉色丝袜足j| 欧美 日韩 精品 国产| 亚洲欧美色中文字幕在线| 精品一品国产午夜福利视频| 天堂动漫精品| 免费观看人在逋| 电影成人av| 最黄视频免费看| 国产一区二区在线观看av| 法律面前人人平等表现在哪些方面| 最近最新中文字幕大全电影3 | 999久久久国产精品视频| 91成人精品电影| 人妻 亚洲 视频| 美女高潮喷水抽搐中文字幕| 女性生殖器流出的白浆| 最近最新中文字幕大全免费视频| 国产av又大| av又黄又爽大尺度在线免费看| 国产成人啪精品午夜网站| 啦啦啦在线免费观看视频4| 午夜福利视频精品| 国产免费视频播放在线视频| 无限看片的www在线观看| 麻豆成人av在线观看| 亚洲专区中文字幕在线| 天天躁狠狠躁夜夜躁狠狠躁| 日韩欧美三级三区| 大香蕉久久网| 香蕉丝袜av| 男女边摸边吃奶| tocl精华| 老司机午夜十八禁免费视频| 亚洲国产中文字幕在线视频| 久久热在线av| 久久精品亚洲熟妇少妇任你| 天堂俺去俺来也www色官网| 99国产极品粉嫩在线观看| 国产淫语在线视频| 欧美日韩中文字幕国产精品一区二区三区 | 成人三级做爰电影| 黄色怎么调成土黄色| www日本在线高清视频| 日韩一区二区三区影片| 女人高潮潮喷娇喘18禁视频| 99国产精品免费福利视频| 肉色欧美久久久久久久蜜桃| 麻豆国产av国片精品| 熟女少妇亚洲综合色aaa.| 精品一区二区三区av网在线观看 | 精品人妻1区二区| 午夜福利乱码中文字幕| 五月开心婷婷网| 啪啪无遮挡十八禁网站| 国产精品美女特级片免费视频播放器 | 免费在线观看日本一区| 香蕉久久夜色| 亚洲成人免费电影在线观看| 十八禁高潮呻吟视频| 亚洲成av片中文字幕在线观看| 这个男人来自地球电影免费观看| 性少妇av在线| 精品久久久精品久久久| 国产无遮挡羞羞视频在线观看| 国产一区有黄有色的免费视频| 国产欧美亚洲国产| 国产成人精品久久二区二区91| 欧美精品一区二区免费开放| 久久久久久亚洲精品国产蜜桃av| 成人18禁在线播放| 久久99热这里只频精品6学生| netflix在线观看网站| 亚洲色图av天堂| 美女视频免费永久观看网站| 韩国精品一区二区三区| 搡老岳熟女国产| 欧美黑人欧美精品刺激| 国产亚洲精品久久久久5区| 日韩人妻精品一区2区三区| av福利片在线| 亚洲国产中文字幕在线视频| 99国产精品99久久久久| 一级黄色大片毛片| 国产精品影院久久| 午夜福利一区二区在线看| 亚洲精品粉嫩美女一区| cao死你这个sao货| 久久午夜亚洲精品久久| 国产单亲对白刺激| 久久国产精品人妻蜜桃| 色播在线永久视频| 动漫黄色视频在线观看| 久久久国产欧美日韩av| 99re6热这里在线精品视频| 在线看a的网站| 国产av又大| 可以免费在线观看a视频的电影网站| 黄色a级毛片大全视频| 精品福利观看| 欧美 日韩 精品 国产| 考比视频在线观看| 亚洲av日韩精品久久久久久密| 久久久久国内视频| 亚洲国产av新网站| 最近最新中文字幕大全电影3 | 色尼玛亚洲综合影院| 大片免费播放器 马上看| 麻豆乱淫一区二区| 亚洲伊人色综图| 国产高清激情床上av| 国产一区二区三区在线臀色熟女 | 看免费av毛片| 亚洲色图av天堂| 日韩 欧美 亚洲 中文字幕| 国产一区二区三区综合在线观看| 欧美在线一区亚洲| 精品国产一区二区久久| 精品少妇久久久久久888优播| av免费在线观看网站| 无遮挡黄片免费观看| 国产欧美日韩精品亚洲av| 丝袜人妻中文字幕| 波多野结衣一区麻豆| 在线av久久热| 精品少妇内射三级| 精品国产超薄肉色丝袜足j| 美国免费a级毛片| 欧美精品亚洲一区二区| 国产一区有黄有色的免费视频| 国产精品.久久久| 午夜老司机福利片| 脱女人内裤的视频| 亚洲熟女毛片儿| 久久久精品免费免费高清| 久久久久久久精品吃奶| 欧美一级毛片孕妇| 成人永久免费在线观看视频 | 99国产精品一区二区蜜桃av | 女警被强在线播放| av国产精品久久久久影院| 欧美亚洲日本最大视频资源| 十八禁网站免费在线| 日韩大码丰满熟妇| 国产片内射在线| 女同久久另类99精品国产91| 18禁国产床啪视频网站| 电影成人av| 色综合欧美亚洲国产小说| 午夜91福利影院| 精品亚洲乱码少妇综合久久| 老司机影院毛片| 免费观看a级毛片全部| 天堂8中文在线网| 动漫黄色视频在线观看| 免费在线观看完整版高清| 99国产精品一区二区蜜桃av | 婷婷丁香在线五月| 欧美国产精品va在线观看不卡| 日本欧美视频一区| 亚洲性夜色夜夜综合| 黄色片一级片一级黄色片| av国产精品久久久久影院| 每晚都被弄得嗷嗷叫到高潮| 国产一区二区三区综合在线观看| 色精品久久人妻99蜜桃| 高清在线国产一区| 黑人巨大精品欧美一区二区mp4| 一区二区日韩欧美中文字幕| 男女免费视频国产| 久久久久国产一级毛片高清牌| 黑人操中国人逼视频| 国产男靠女视频免费网站| 肉色欧美久久久久久久蜜桃| 极品人妻少妇av视频| 色播在线永久视频| 精品人妻1区二区| 香蕉丝袜av| 熟女少妇亚洲综合色aaa.| 777久久人妻少妇嫩草av网站| 午夜福利欧美成人| 国产成人精品久久二区二区91| 欧美成人免费av一区二区三区 | 国产成人一区二区三区免费视频网站| 亚洲av成人一区二区三| 男女午夜视频在线观看| 亚洲性夜色夜夜综合| av片东京热男人的天堂| 成人永久免费在线观看视频 | 日本一区二区免费在线视频| 啦啦啦在线免费观看视频4| 天天躁夜夜躁狠狠躁躁| 免费高清在线观看日韩| 国产97色在线日韩免费| 日韩免费高清中文字幕av| 久久国产精品男人的天堂亚洲| 国产成人精品无人区| 一级毛片电影观看| 亚洲精品在线美女| 亚洲精品美女久久av网站| bbb黄色大片| 搡老熟女国产l中国老女人| 亚洲成人免费av在线播放| cao死你这个sao货| 久久人妻av系列| 丝袜在线中文字幕| 国产精品av久久久久免费| 欧美日韩亚洲国产一区二区在线观看 | 天天躁狠狠躁夜夜躁狠狠躁| 国产免费福利视频在线观看| www.精华液| 亚洲国产欧美日韩在线播放| 午夜福利在线免费观看网站| 可以免费在线观看a视频的电影网站| 中文字幕高清在线视频| 精品亚洲乱码少妇综合久久| 黄频高清免费视频| 亚洲国产中文字幕在线视频| 美女视频免费永久观看网站| 国产av精品麻豆| 国产成人啪精品午夜网站| 建设人人有责人人尽责人人享有的| 亚洲熟妇熟女久久| 操美女的视频在线观看| 亚洲自偷自拍图片 自拍| 国产又色又爽无遮挡免费看| 视频在线观看一区二区三区| 日本av手机在线免费观看| 老汉色∧v一级毛片| 成人免费观看视频高清| 国产区一区二久久| 成年人午夜在线观看视频| a级片在线免费高清观看视频| 18在线观看网站| 青青草视频在线视频观看| 国产又色又爽无遮挡免费看| 一区在线观看完整版| 国产视频一区二区在线看| 人成视频在线观看免费观看| 免费一级毛片在线播放高清视频 | 日韩欧美一区视频在线观看| 如日韩欧美国产精品一区二区三区| 色播在线永久视频| 久久亚洲精品不卡| 日本av免费视频播放| 在线观看舔阴道视频| 国产主播在线观看一区二区| 十分钟在线观看高清视频www| 日韩免费高清中文字幕av| 日韩欧美一区视频在线观看| 欧美乱妇无乱码| 精品久久久精品久久久| 久久精品aⅴ一区二区三区四区| 亚洲熟女精品中文字幕| 一本综合久久免费| 精品乱码久久久久久99久播| 日本wwww免费看| 黄色视频在线播放观看不卡| 一二三四社区在线视频社区8| 日韩一卡2卡3卡4卡2021年| av免费在线观看网站| videosex国产| 国产精品久久电影中文字幕 | 久久人妻av系列| 亚洲av成人不卡在线观看播放网| 欧美在线一区亚洲| 国产av国产精品国产| 精品视频人人做人人爽| 午夜精品久久久久久毛片777| 色婷婷久久久亚洲欧美| 高清视频免费观看一区二区| 亚洲精品美女久久av网站| 纵有疾风起免费观看全集完整版| 两个人看的免费小视频| 精品人妻1区二区| 国产亚洲精品久久久久5区| 91大片在线观看| 乱人伦中国视频| 亚洲国产欧美一区二区综合| 国产成人免费无遮挡视频| 成年动漫av网址| 操出白浆在线播放| 黄色视频在线播放观看不卡| 成年人免费黄色播放视频| 91麻豆av在线| 三上悠亚av全集在线观看| 老鸭窝网址在线观看| 精品卡一卡二卡四卡免费| 久久天堂一区二区三区四区| 日本vs欧美在线观看视频| 另类精品久久| 一夜夜www| 人人澡人人妻人| 丝袜喷水一区| 一区二区三区精品91| 精品高清国产在线一区| 欧美变态另类bdsm刘玥| 亚洲人成伊人成综合网2020| 九色亚洲精品在线播放| 欧美日韩中文字幕国产精品一区二区三区 | 日韩一区二区三区影片| 黄色视频不卡| 亚洲国产中文字幕在线视频| 一边摸一边抽搐一进一出视频| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久国内视频| 一级毛片电影观看| 欧美成人午夜精品| 成在线人永久免费视频| 久久ye,这里只有精品| 日本五十路高清| 超色免费av| 中文字幕精品免费在线观看视频| 一二三四社区在线视频社区8| 亚洲国产中文字幕在线视频| 黑人巨大精品欧美一区二区蜜桃| 久久精品成人免费网站| 久久精品亚洲精品国产色婷小说| 亚洲成av片中文字幕在线观看| 可以免费在线观看a视频的电影网站| 国产一区二区三区视频了| 亚洲七黄色美女视频| 丁香六月欧美| 国产精品国产av在线观看| 国产97色在线日韩免费| 国产成人欧美| 一级,二级,三级黄色视频| 欧美日韩国产mv在线观看视频| 久久天堂一区二区三区四区| 精品久久久久久电影网| 日本黄色视频三级网站网址 | 香蕉丝袜av| 在线观看免费高清a一片| 中文字幕人妻丝袜一区二区| 美女扒开内裤让男人捅视频| 乱人伦中国视频| 国产av又大| 日本黄色视频三级网站网址 | 国产成人欧美在线观看 | 久久毛片免费看一区二区三区| 亚洲国产欧美在线一区| 亚洲美女黄片视频| 国产成人免费观看mmmm| 精品亚洲成国产av| 久久精品国产99精品国产亚洲性色 | 真人做人爱边吃奶动态| 视频区欧美日本亚洲| 咕卡用的链子| 91av网站免费观看| 国产精品久久久久久精品电影小说| 久久亚洲真实| 女同久久另类99精品国产91| 十八禁人妻一区二区| 精品一区二区三区视频在线观看免费 | 大香蕉久久成人网| 丝袜在线中文字幕| 国产国语露脸激情在线看| 久久精品国产a三级三级三级| 无人区码免费观看不卡 | 最新在线观看一区二区三区| 免费看a级黄色片| 日本撒尿小便嘘嘘汇集6| 午夜福利视频在线观看免费| 免费不卡黄色视频| 天天影视国产精品| 99热网站在线观看| 一个人免费看片子| 亚洲综合色网址| 两个人免费观看高清视频| 国产日韩一区二区三区精品不卡| 国产精品国产av在线观看| √禁漫天堂资源中文www| 妹子高潮喷水视频| 交换朋友夫妻互换小说| 9热在线视频观看99| av视频免费观看在线观看| 熟女少妇亚洲综合色aaa.| 在线观看一区二区三区激情| 国产成人一区二区三区免费视频网站| 妹子高潮喷水视频| 电影成人av| 欧美精品啪啪一区二区三区| 国产有黄有色有爽视频| 黄片播放在线免费| 国产xxxxx性猛交| 操美女的视频在线观看| 老司机靠b影院| 亚洲av欧美aⅴ国产| 老熟妇乱子伦视频在线观看| 超色免费av| 欧美 日韩 精品 国产| 精品国产一区二区久久| 亚洲国产毛片av蜜桃av| www.精华液| 亚洲男人天堂网一区| 一区二区三区国产精品乱码| 国产亚洲精品久久久久5区| 国产精品免费视频内射| 亚洲人成电影观看| 每晚都被弄得嗷嗷叫到高潮| 精品久久久久久电影网| 搡老岳熟女国产| 91大片在线观看| 国产精品免费一区二区三区在线 | 99riav亚洲国产免费| 97人妻天天添夜夜摸| 国产福利在线免费观看视频| 不卡av一区二区三区| 精品福利永久在线观看| 国产视频一区二区在线看| 一区在线观看完整版| 国产男靠女视频免费网站| 美女扒开内裤让男人捅视频| 91字幕亚洲| 桃花免费在线播放| 人妻 亚洲 视频| 日本av手机在线免费观看| 国产一区二区三区在线臀色熟女 | 男女免费视频国产| 黄色视频在线播放观看不卡| 男女边摸边吃奶| 热re99久久精品国产66热6| 女警被强在线播放| 91成人精品电影| 免费久久久久久久精品成人欧美视频| 1024香蕉在线观看| 国精品久久久久久国模美| 国产精品av久久久久免费| 国产精品久久久久久精品电影小说| 欧美日韩成人在线一区二区| 纯流量卡能插随身wifi吗| 99热网站在线观看| 成年动漫av网址| 人人澡人人妻人| 亚洲av电影在线进入| 欧美激情久久久久久爽电影 | 制服人妻中文乱码| 欧美激情极品国产一区二区三区| av电影中文网址| 久久亚洲真实| 亚洲国产欧美日韩在线播放| 99久久精品国产亚洲精品| 国产主播在线观看一区二区| 青青草视频在线视频观看| 国产男女超爽视频在线观看| 建设人人有责人人尽责人人享有的| 婷婷丁香在线五月| 91成年电影在线观看| 久久久久久免费高清国产稀缺| 熟女少妇亚洲综合色aaa.| 黄色视频,在线免费观看| 老司机靠b影院| 丰满人妻熟妇乱又伦精品不卡| 丝袜美腿诱惑在线| 悠悠久久av| 中文字幕最新亚洲高清| 亚洲欧美一区二区三区黑人| 国产深夜福利视频在线观看| 黄色a级毛片大全视频|