• <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ì)與概率密度演化研究
    多毛熟女@视频| 每晚都被弄得嗷嗷叫到高潮| 欧美精品av麻豆av| 女人久久www免费人成看片| 久久久久久免费高清国产稀缺| 99九九在线精品视频| 十八禁人妻一区二区| 纯流量卡能插随身wifi吗| 久久热在线av| 高清在线国产一区| 国产成人一区二区三区免费视频网站| 亚洲人成电影观看| 这个男人来自地球电影免费观看| 久久99一区二区三区| 大陆偷拍与自拍| 国产欧美日韩一区二区精品| 久热这里只有精品99| 男人添女人高潮全过程视频| 热re99久久精品国产66热6| 久久久久国产精品人妻一区二区| 亚洲精品国产av蜜桃| 久久精品国产a三级三级三级| 在线永久观看黄色视频| 十八禁网站免费在线| 亚洲精品国产精品久久久不卡| 亚洲精品美女久久av网站| 黑人巨大精品欧美一区二区蜜桃| 久久久久久久国产电影| 亚洲av成人不卡在线观看播放网 | 日本黄色日本黄色录像| 国产精品自产拍在线观看55亚洲 | 欧美大码av| 亚洲 欧美一区二区三区| xxxhd国产人妻xxx| 91老司机精品| 啦啦啦在线免费观看视频4| 精品少妇久久久久久888优播| 天天躁日日躁夜夜躁夜夜| 搡老熟女国产l中国老女人| 一本大道久久a久久精品| 不卡av一区二区三区| 久久国产精品影院| 欧美国产精品va在线观看不卡| 美女福利国产在线| 一本综合久久免费| 一本久久精品| 手机成人av网站| 久久人人爽人人片av| 777久久人妻少妇嫩草av网站| 欧美 亚洲 国产 日韩一| 亚洲第一欧美日韩一区二区三区 | 精品卡一卡二卡四卡免费| 首页视频小说图片口味搜索| 这个男人来自地球电影免费观看| 黑丝袜美女国产一区| 他把我摸到了高潮在线观看 | 免费少妇av软件| av在线app专区| 成年av动漫网址| 天堂俺去俺来也www色官网| 亚洲国产欧美日韩在线播放| 一二三四在线观看免费中文在| 久久久久久免费高清国产稀缺| 黑人巨大精品欧美一区二区mp4| 国产1区2区3区精品| 亚洲欧洲精品一区二区精品久久久| 首页视频小说图片口味搜索| 成人三级做爰电影| 国产激情久久老熟女| 国产人伦9x9x在线观看| 亚洲av电影在线观看一区二区三区| 在线观看舔阴道视频| 免费一级毛片在线播放高清视频 | 在线观看免费视频网站a站| 午夜福利一区二区在线看| 80岁老熟妇乱子伦牲交| 欧美午夜高清在线| 国产精品国产av在线观看| 国产精品国产av在线观看| 日韩 亚洲 欧美在线| av在线播放精品| 亚洲欧洲精品一区二区精品久久久| 汤姆久久久久久久影院中文字幕| 成人国语在线视频| 大型av网站在线播放| 日韩有码中文字幕| 亚洲中文av在线| 日韩人妻精品一区2区三区| 80岁老熟妇乱子伦牲交| 国产欧美日韩精品亚洲av| 亚洲 国产 在线| 欧美日韩av久久| 国产精品久久久久久精品电影小说| 午夜久久久在线观看| 久久影院123| 两个人免费观看高清视频| 国产一区二区激情短视频 | 亚洲欧美成人综合另类久久久| 亚洲人成电影观看| 交换朋友夫妻互换小说| 免费黄频网站在线观看国产| 国产精品影院久久| 一区二区三区激情视频| 免费在线观看日本一区| 各种免费的搞黄视频| tocl精华| 成人亚洲精品一区在线观看| 80岁老熟妇乱子伦牲交| 啦啦啦免费观看视频1| 欧美精品亚洲一区二区| 女人爽到高潮嗷嗷叫在线视频| 精品国内亚洲2022精品成人 | 男女国产视频网站| 日本av免费视频播放| 成人黄色视频免费在线看| 国产亚洲欧美在线一区二区| 欧美97在线视频| 日韩精品免费视频一区二区三区| 欧美成人午夜精品| 夜夜骑夜夜射夜夜干| 中文字幕av电影在线播放| 操美女的视频在线观看| 91老司机精品| 99热网站在线观看| 久久99热这里只频精品6学生| 日本撒尿小便嘘嘘汇集6| 免费人妻精品一区二区三区视频| 亚洲精品久久久久久婷婷小说| 十八禁网站网址无遮挡| 欧美xxⅹ黑人| 亚洲男人天堂网一区| 精品人妻1区二区| av在线app专区| 亚洲精品一卡2卡三卡4卡5卡 | 美女国产高潮福利片在线看| 一级毛片女人18水好多| 一本久久精品| 久久午夜综合久久蜜桃| 亚洲 欧美一区二区三区| 精品久久蜜臀av无| 久久久精品国产亚洲av高清涩受| 日本撒尿小便嘘嘘汇集6| 国产精品 国内视频| 一本—道久久a久久精品蜜桃钙片| 亚洲av欧美aⅴ国产| 亚洲第一欧美日韩一区二区三区 | 丝袜喷水一区| 精品一品国产午夜福利视频| 捣出白浆h1v1| 成人黄色视频免费在线看| 伦理电影免费视频| a 毛片基地| 老熟女久久久| 蜜桃在线观看..| 国产精品影院久久| 精品亚洲乱码少妇综合久久| 丰满饥渴人妻一区二区三| 美国免费a级毛片| 久热爱精品视频在线9| 啦啦啦视频在线资源免费观看| 女人精品久久久久毛片| 久久久国产一区二区| av福利片在线| 亚洲国产精品一区三区| 十分钟在线观看高清视频www| 大片电影免费在线观看免费| 麻豆国产av国片精品| 亚洲美女黄色视频免费看| 在线观看免费日韩欧美大片| 亚洲国产欧美一区二区综合| 亚洲精品日韩在线中文字幕| 色94色欧美一区二区| 久久99一区二区三区| 日韩三级视频一区二区三区| 久久天堂一区二区三区四区| 日韩熟女老妇一区二区性免费视频| 考比视频在线观看| 91av网站免费观看| 91精品三级在线观看| 男女床上黄色一级片免费看| 日韩欧美一区二区三区在线观看 | 日韩熟女老妇一区二区性免费视频| 成人亚洲精品一区在线观看| 亚洲视频免费观看视频| 捣出白浆h1v1| 久久精品亚洲av国产电影网| 视频区欧美日本亚洲| 人妻一区二区av| 国产一区二区三区在线臀色熟女 | 国产在线视频一区二区| 男女之事视频高清在线观看| 免费高清在线观看视频在线观看| 两人在一起打扑克的视频| 多毛熟女@视频| 咕卡用的链子| 少妇粗大呻吟视频| 国产精品久久久av美女十八| www.av在线官网国产| 国产精品 国内视频| 欧美中文综合在线视频| 免费日韩欧美在线观看| 99久久人妻综合| 视频区欧美日本亚洲| 国产极品粉嫩免费观看在线| 精品欧美一区二区三区在线| 久久国产亚洲av麻豆专区| 国产无遮挡羞羞视频在线观看| 国产主播在线观看一区二区| 男女下面插进去视频免费观看| 午夜激情久久久久久久| 日韩免费高清中文字幕av| 亚洲国产成人一精品久久久| 少妇被粗大的猛进出69影院| av超薄肉色丝袜交足视频| 欧美日韩亚洲国产一区二区在线观看 | 天堂8中文在线网| 王馨瑶露胸无遮挡在线观看| 色婷婷久久久亚洲欧美| 嫩草影视91久久| 久久天堂一区二区三区四区| 欧美精品人与动牲交sv欧美| 欧美另类亚洲清纯唯美| 精品卡一卡二卡四卡免费| 女人高潮潮喷娇喘18禁视频| 在线精品无人区一区二区三| 日韩人妻精品一区2区三区| 人人妻人人澡人人看| 亚洲avbb在线观看| 久久热在线av| 9191精品国产免费久久| a 毛片基地| 悠悠久久av| 在线观看一区二区三区激情| 黄色视频不卡| av网站免费在线观看视频| 日韩一卡2卡3卡4卡2021年| 国产亚洲午夜精品一区二区久久| 欧美性长视频在线观看| 亚洲成人免费av在线播放| 国产一区二区三区在线臀色熟女 | 首页视频小说图片口味搜索| 久久精品国产综合久久久| 成人av一区二区三区在线看 | 搡老乐熟女国产| 亚洲中文字幕日韩| 黄片小视频在线播放| 麻豆av在线久日| 少妇的丰满在线观看| 日日夜夜操网爽| 欧美亚洲日本最大视频资源| 久久天堂一区二区三区四区| 少妇裸体淫交视频免费看高清 | 久久久久久久久免费视频了| 亚洲欧美日韩高清在线视频 | 免费日韩欧美在线观看| 亚洲成人免费av在线播放| 俄罗斯特黄特色一大片| 最新在线观看一区二区三区| 精品国产乱子伦一区二区三区 | 成人免费观看视频高清| 香蕉国产在线看| 亚洲第一欧美日韩一区二区三区 | 欧美亚洲 丝袜 人妻 在线| 一级片'在线观看视频| 精品少妇久久久久久888优播| 国产精品一区二区精品视频观看| 国产在线观看jvid| 狠狠狠狠99中文字幕| 中文字幕色久视频| 亚洲va日本ⅴa欧美va伊人久久 | 久久人妻福利社区极品人妻图片| 国产成人免费观看mmmm| 欧美日韩国产mv在线观看视频| 日韩欧美国产一区二区入口| 一本一本久久a久久精品综合妖精| 99国产精品一区二区蜜桃av | 国产一级毛片在线| 亚洲成av片中文字幕在线观看| 亚洲精品自拍成人| 法律面前人人平等表现在哪些方面 | 国产精品久久久久成人av| 伊人久久大香线蕉亚洲五| 在线天堂中文资源库| 男人操女人黄网站| 91老司机精品| 天堂俺去俺来也www色官网| 18禁国产床啪视频网站| 国产色视频综合| 国产精品自产拍在线观看55亚洲 | 欧美国产精品va在线观看不卡| 男女无遮挡免费网站观看| 久久久欧美国产精品| 男人操女人黄网站| 一本色道久久久久久精品综合| 久久久久久久久免费视频了| 国产日韩欧美亚洲二区| 国产真人三级小视频在线观看| 交换朋友夫妻互换小说| 国产国语露脸激情在线看| 成年女人毛片免费观看观看9 | 黄片小视频在线播放| 国产亚洲欧美在线一区二区| 欧美人与性动交α欧美精品济南到| 一级毛片女人18水好多| 咕卡用的链子| 日韩欧美一区二区三区在线观看 | 桃红色精品国产亚洲av| 久久人人爽av亚洲精品天堂| 午夜日韩欧美国产| 好男人电影高清在线观看| 亚洲第一av免费看| 美女脱内裤让男人舔精品视频| 欧美亚洲 丝袜 人妻 在线| 精品一区二区三区四区五区乱码| 国产亚洲午夜精品一区二区久久| 国产成人免费无遮挡视频| 女性被躁到高潮视频| 欧美日韩亚洲高清精品| 国产xxxxx性猛交| 亚洲性夜色夜夜综合| 麻豆av在线久日| 国产男人的电影天堂91| 天堂8中文在线网| 曰老女人黄片| 在线观看人妻少妇| 成人国产一区最新在线观看| videos熟女内射| av欧美777| 亚洲中文av在线| 国产一卡二卡三卡精品| 大陆偷拍与自拍| 亚洲久久久国产精品| 69av精品久久久久久 | 久久免费观看电影| 精品久久久久久电影网| 中文字幕最新亚洲高清| 精品久久久精品久久久| 男人添女人高潮全过程视频| 在线精品无人区一区二区三| 操美女的视频在线观看| 中文字幕最新亚洲高清| 99热全是精品| 美女中出高潮动态图| 高清视频免费观看一区二区| 国产精品自产拍在线观看55亚洲 | 国产精品熟女久久久久浪| 欧美午夜高清在线| 美女脱内裤让男人舔精品视频| 各种免费的搞黄视频| 久久久久久久国产电影| 国产精品久久久人人做人人爽| 中文字幕人妻丝袜制服| 99香蕉大伊视频| 夫妻午夜视频| 欧美精品啪啪一区二区三区 | 精品福利永久在线观看| 黄色片一级片一级黄色片| 亚洲色图综合在线观看| 一级黄色大片毛片| 免费黄频网站在线观看国产| 久久精品人人爽人人爽视色| 亚洲国产欧美一区二区综合| 亚洲av片天天在线观看| 韩国高清视频一区二区三区| 精品少妇黑人巨大在线播放| 最近最新中文字幕大全免费视频| 夜夜夜夜夜久久久久| 免费在线观看日本一区| av免费在线观看网站| 天天影视国产精品| 欧美日韩中文字幕国产精品一区二区三区 | 国产福利在线免费观看视频| 亚洲国产av新网站| 人人妻人人添人人爽欧美一区卜| 久热爱精品视频在线9| 一个人免费看片子| 日韩欧美免费精品| 免费在线观看黄色视频的| 美国免费a级毛片| 午夜视频精品福利| av又黄又爽大尺度在线免费看| 国产亚洲av高清不卡| 黄色视频,在线免费观看| 美女扒开内裤让男人捅视频| 午夜精品国产一区二区电影| 51午夜福利影视在线观看| 久久久久久久久免费视频了| tube8黄色片| 国产极品粉嫩免费观看在线| 亚洲中文字幕日韩| 亚洲av电影在线进入| 大码成人一级视频| 99久久99久久久精品蜜桃| 黄色 视频免费看| 国产一区二区激情短视频 | 激情视频va一区二区三区| 久久精品熟女亚洲av麻豆精品| 少妇裸体淫交视频免费看高清 | 一本大道久久a久久精品| 男人舔女人的私密视频| 久久精品亚洲av国产电影网| 中文字幕最新亚洲高清| 老熟妇乱子伦视频在线观看 | 成在线人永久免费视频| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品成人av观看孕妇| 国产免费av片在线观看野外av| 69av精品久久久久久 | 亚洲第一欧美日韩一区二区三区 | 日韩制服骚丝袜av| 久久国产亚洲av麻豆专区| 亚洲专区中文字幕在线| 久久久水蜜桃国产精品网| 欧美中文综合在线视频| 色老头精品视频在线观看| 国产在线观看jvid| 日本一区二区免费在线视频| 一级毛片女人18水好多| 日本五十路高清| 满18在线观看网站| 精品久久蜜臀av无| 久久久久久久国产电影| 最近中文字幕2019免费版| 久久精品久久久久久噜噜老黄| 青青草视频在线视频观看| 日韩 欧美 亚洲 中文字幕| 少妇被粗大的猛进出69影院| 久久狼人影院| 国产成人免费观看mmmm| 亚洲全国av大片| 成人手机av| 侵犯人妻中文字幕一二三四区| 免费女性裸体啪啪无遮挡网站| 中文字幕av电影在线播放| 人人妻人人爽人人添夜夜欢视频| 亚洲欧洲日产国产| 纯流量卡能插随身wifi吗| 亚洲国产欧美一区二区综合| 欧美精品啪啪一区二区三区 | 日韩制服骚丝袜av| 久久久久久免费高清国产稀缺| 亚洲欧美一区二区三区黑人| 99久久99久久久精品蜜桃| 亚洲男人天堂网一区| 欧美国产精品一级二级三级| 国产真人三级小视频在线观看| 性少妇av在线| 国产亚洲一区二区精品| 国产精品久久久久久人妻精品电影 | 精品久久蜜臀av无| 丝瓜视频免费看黄片| 大香蕉久久成人网| 亚洲va日本ⅴa欧美va伊人久久 | 免费观看a级毛片全部| 一级a爱视频在线免费观看| 美女脱内裤让男人舔精品视频| 国产av国产精品国产| 99国产精品免费福利视频| 亚洲国产av影院在线观看| 国产视频一区二区在线看| 国产免费福利视频在线观看| 午夜福利影视在线免费观看| 亚洲精品粉嫩美女一区| 午夜免费鲁丝| 国产91精品成人一区二区三区 | 中文字幕最新亚洲高清| 国产日韩欧美在线精品| 精品熟女少妇八av免费久了| www日本在线高清视频| 男女之事视频高清在线观看| 亚洲av成人一区二区三| 19禁男女啪啪无遮挡网站| av天堂久久9| 美女中出高潮动态图| 成人国语在线视频| 国产深夜福利视频在线观看| 免费人妻精品一区二区三区视频| 国产精品99久久99久久久不卡| 中文字幕人妻丝袜制服| 中文字幕av电影在线播放| 日韩欧美一区视频在线观看| 亚洲国产av影院在线观看| 美女视频免费永久观看网站| 91国产中文字幕| 人妻人人澡人人爽人人| 亚洲精品一二三| 国产男人的电影天堂91| 婷婷成人精品国产| 大型av网站在线播放| 日日爽夜夜爽网站| 午夜免费观看性视频| 女警被强在线播放| 亚洲精品成人av观看孕妇| 黑丝袜美女国产一区| 欧美久久黑人一区二区| tube8黄色片| 18禁黄网站禁片午夜丰满| 在线观看www视频免费| a级毛片在线看网站| 欧美人与性动交α欧美软件| 黄色视频不卡| 国产日韩一区二区三区精品不卡| 午夜激情av网站| 99精国产麻豆久久婷婷| 精品国内亚洲2022精品成人 | 久久精品亚洲av国产电影网| av免费在线观看网站| 男女免费视频国产| 精品国产国语对白av| 十八禁网站免费在线| 国产精品免费大片| 国产精品久久久av美女十八| 少妇人妻久久综合中文| 欧美 日韩 精品 国产| 中国美女看黄片| 亚洲国产精品成人久久小说| 亚洲性夜色夜夜综合| 看免费av毛片| h视频一区二区三区| 桃花免费在线播放| 99久久人妻综合| 精品国产乱码久久久久久小说| 999精品在线视频| 麻豆av在线久日| 色老头精品视频在线观看| 亚洲一码二码三码区别大吗| 又黄又粗又硬又大视频| 大陆偷拍与自拍| 91精品国产国语对白视频| 亚洲成人国产一区在线观看| 多毛熟女@视频| 久久精品熟女亚洲av麻豆精品| 欧美黄色淫秽网站| 精品亚洲乱码少妇综合久久| 一进一出抽搐动态| 女人被躁到高潮嗷嗷叫费观| 日本av免费视频播放| 国产精品1区2区在线观看. | 久久久精品国产亚洲av高清涩受| 少妇粗大呻吟视频| 日本av免费视频播放| 99热全是精品| 18在线观看网站| 午夜激情av网站| 成年美女黄网站色视频大全免费| 精品国产国语对白av| 国产一区二区激情短视频 | 99精品欧美一区二区三区四区| 欧美少妇被猛烈插入视频| 精品一区二区三卡| 亚洲综合色网址| 日韩欧美一区视频在线观看| 国产高清国产精品国产三级| 久久毛片免费看一区二区三区| 搡老乐熟女国产| 另类亚洲欧美激情| 久久久精品94久久精品| 美女午夜性视频免费| 桃花免费在线播放| 涩涩av久久男人的天堂| 国产成人一区二区三区免费视频网站| 黑人巨大精品欧美一区二区蜜桃| 亚洲欧美成人综合另类久久久| 99久久人妻综合| 性少妇av在线| 国产日韩欧美视频二区| 中文字幕另类日韩欧美亚洲嫩草| 美国免费a级毛片| 精品少妇黑人巨大在线播放| 久久久国产一区二区| 男人添女人高潮全过程视频| 大码成人一级视频| 欧美在线黄色| 美女扒开内裤让男人捅视频| 日韩视频在线欧美| 大片电影免费在线观看免费| 精品一区二区三卡| 91麻豆av在线| 在线观看免费视频网站a站| 欧美另类一区| 亚洲视频免费观看视频| 两个人看的免费小视频| www.999成人在线观看| 久久99热这里只频精品6学生| 十八禁网站网址无遮挡| 黑人欧美特级aaaaaa片| 一个人免费看片子| 高清av免费在线| 亚洲精品久久成人aⅴ小说| 欧美少妇被猛烈插入视频| 午夜福利乱码中文字幕| 亚洲欧洲日产国产| 欧美变态另类bdsm刘玥| 精品一区在线观看国产| 国产欧美日韩一区二区三 | 久久久久网色| 永久免费av网站大全| 水蜜桃什么品种好| 18在线观看网站| 天堂俺去俺来也www色官网| 亚洲熟女精品中文字幕| 飞空精品影院首页| 午夜精品久久久久久毛片777| 免费在线观看影片大全网站| 18禁黄网站禁片午夜丰满| 人人妻,人人澡人人爽秒播| 女人精品久久久久毛片| 国产日韩欧美在线精品| 欧美中文综合在线视频| 老汉色av国产亚洲站长工具| 一级毛片女人18水好多| 国产一区二区三区综合在线观看| 色视频在线一区二区三区| 最近中文字幕2019免费版| 久久精品国产综合久久久|