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

    基于橫向加速度監(jiān)測的高速鐵路橋梁動力性能異常預(yù)警方法

    2016-10-18 11:55:52耿方方丁幼亮趙瀚瑋吳來義岳青
    鐵道建筑 2016年9期
    關(guān)鍵詞:方根值主梁大橋

    耿方方,丁幼亮,趙瀚瑋,吳來義,岳青

    (1.東南大學(xué)混凝土及預(yù)應(yīng)力混凝土結(jié)構(gòu)教育部重點(diǎn)實(shí)驗(yàn)室,江蘇南京210096;2.南京工程學(xué)院,江蘇南京211167;3.河北省大型結(jié)構(gòu)健康診斷與控制實(shí)驗(yàn)室,河北石家莊050043;4.中鐵大橋(南京)橋隧診治有限公司,江蘇南京210032)

    基于橫向加速度監(jiān)測的高速鐵路橋梁動力性能異常預(yù)警方法

    耿方方1,2,丁幼亮1,3,趙瀚瑋1,吳來義4,岳青4

    (1.東南大學(xué)混凝土及預(yù)應(yīng)力混凝土結(jié)構(gòu)教育部重點(diǎn)實(shí)驗(yàn)室,江蘇南京210096;2.南京工程學(xué)院,江蘇南京211167;3.河北省大型結(jié)構(gòu)健康診斷與控制實(shí)驗(yàn)室,河北石家莊050043;4.中鐵大橋(南京)橋隧診治有限公司,江蘇南京210032)

    基于大勝關(guān)大橋的主梁橫向加速度監(jiān)測數(shù)據(jù),采用主成分分析方法建立了高速鐵路橋梁動力性能異常預(yù)警方法。采用主梁橫向加速度的均方根值作為反映橋梁動力性能的監(jiān)測參數(shù),以其互相關(guān)性模型作為橋梁動力性能評價模型,并采用主成分分析方法有效地分離了環(huán)境因素變化對橋梁動力性能評價模型的影響,建立了橋梁動力性能異常預(yù)警指標(biāo)。分析結(jié)果表明:大勝關(guān)大橋主梁不同測點(diǎn)之間均方根值的互相關(guān)性可以采用二次多項(xiàng)式擬合;所建立的預(yù)警方法可以敏感地發(fā)現(xiàn)主梁橫向加速度的異常變化,從而對橋梁動力性能及時做出預(yù)警。

    高速鐵路橋梁;振動監(jiān)測;主成分分析;環(huán)境因素;實(shí)時預(yù)警

    京滬高鐵南京大勝關(guān)長江大橋?yàn)榫咚勹F路和滬漢蓉鐵路共用的跨江通道,是世界首座6線鐵路大橋。大勝關(guān)大橋由于跨度大、荷載大、行車速度高,列車在橋上運(yùn)行時可能產(chǎn)生比較大的振動,對高速列車運(yùn)行安全的影響必須予以重視。特別是大勝關(guān)大橋在長期服役過程中,由于受環(huán)境侵蝕、材料老化、荷載的長期效應(yīng)與疲勞效應(yīng)等因素影響,導(dǎo)致了橋梁服役性能特別是動力性能發(fā)生退化,需要采取必要的手段和措施保障大勝關(guān)大橋的列車走行安全性。

    大勝關(guān)大橋結(jié)構(gòu)健康監(jiān)測系統(tǒng)為該橋動力性能監(jiān)測與異常預(yù)警提供了新的契機(jī)。該橋結(jié)構(gòu)健康監(jiān)測系統(tǒng)主要監(jiān)測結(jié)構(gòu)振動、應(yīng)力、變形響應(yīng)以及與橋址環(huán)境相關(guān)的風(fēng)速、溫度等參數(shù)[1]。結(jié)構(gòu)健康監(jiān)測系統(tǒng)通過分析實(shí)測數(shù)據(jù)建立大橋的連續(xù)評價體系,并預(yù)測大橋的未來變化趨勢,實(shí)現(xiàn)以健康監(jiān)測系統(tǒng)為平臺的大橋養(yǎng)護(hù)管理評價系統(tǒng)[2]。本文依據(jù)京滬高鐵南京大勝關(guān)大橋2014年的主梁實(shí)測橫向加速度數(shù)據(jù),研究基于主成分分析(Principal Component Analysis,PCA)的動力性能監(jiān)測與異常預(yù)警方法,重點(diǎn)討論了橋梁動力性能監(jiān)測參數(shù)的選擇、動力性能評價模型的建立、動力性能異常預(yù)警指標(biāo)的構(gòu)建等,通過長期監(jiān)測數(shù)據(jù)分析探討了該方法的可行性與有效性。

    1 基于橫向加速度監(jiān)測的橋梁動力性能異常預(yù)警原理

    一般采用脫軌系數(shù)和輪重減載率來判斷列車走行安全性。文獻(xiàn)[3]指出現(xiàn)行脫軌系數(shù)和輪重減載率的限值不能防止脫軌發(fā)生,列車脫軌的根本原因是列車橋梁(軌道)時變系統(tǒng)的橫向振動喪失穩(wěn)定。本文將文獻(xiàn)[3]建立的列車脫軌能量隨機(jī)分析理論和鐵路橋梁結(jié)構(gòu)健康監(jiān)測系統(tǒng)相結(jié)合,建立基于橫向加速度監(jiān)測的橋梁動力性能異常預(yù)警方法。該方法的基本原理是將列車和橋梁視為一個整體系統(tǒng),對高速列車通過時主梁的橫向振動響應(yīng)進(jìn)行分析,及時發(fā)現(xiàn)橋梁因結(jié)構(gòu)動力性能退化造成的主梁橫向振動響應(yīng)異常變化,并作出預(yù)警。

    基于橫向加速度監(jiān)測的橋梁動力性能異常預(yù)警方法包括3項(xiàng)研究內(nèi)容:

    1)提取大橋動力性能監(jiān)測參數(shù)。本文將考察橫向加速度峰值和均方根值這2類監(jiān)測參數(shù)的適用性。

    2)構(gòu)建大橋動力性能評價模型。由于高速列車通過時列車-軌道-橋梁系統(tǒng)的橫向振動響應(yīng)具有顯著的隨機(jī)性,故需要提取動力性能監(jiān)測參數(shù),構(gòu)建動力性能評價模型。

    3)建立大橋動力性能異常預(yù)警方法。采用統(tǒng)計(jì)模式識別的方法對大橋動力性能評價模型進(jìn)行長期監(jiān)測并建立異常預(yù)警指標(biāo),從而較好地識別出大橋動力性能是否發(fā)生異常變化。

    2 大勝關(guān)大橋動力性能監(jiān)測與異常預(yù)警方法

    2.1主梁橫向振動監(jiān)測系統(tǒng)

    大勝關(guān)大橋主橋?yàn)橐宦?lián)(108+192+2×336+192 +108)m鋼桁拱橋,主橋立面見圖1(a)。為了監(jiān)測主梁2主跨在高速列車作用下的橫向加速度響應(yīng),分別在2主跨跨中斷面處各安裝1個橫向加速度傳感器,分別為傳感器JSD-11-04和JSD-15-06,傳感器布置見圖1(b)。本文依據(jù)2014年中310 d的橫向加速度數(shù)據(jù)研究大勝關(guān)大橋動力性能監(jiān)測與異常預(yù)警方法。

    圖1 主橋立面及橫向加速度傳感器布置

    2.2主梁橫向加速度監(jiān)測結(jié)果與分析

    高速列車通過時的橫向加速度時程曲線見圖2。從圖2可以看出,高速列車通過時2主跨跨中的橫向加速度具有顯著的非平穩(wěn)振動特性,并且主梁第1跨跨中的橫向加速度響應(yīng)明顯小于第2跨跨中。說明雖然大勝關(guān)大橋2主跨的結(jié)構(gòu)布置相同,但是由于軌道橫向不平順等原因,2主跨在列車-軌道-橋梁系統(tǒng)的橫向振動特性上存在較大的差異,有必要針對大勝關(guān)大橋主梁橫向加速度響應(yīng)進(jìn)行長期監(jiān)測,以監(jiān)控高速列車通過時橋梁的橫向振動狀態(tài)。

    本文采用加速度幅值作為大勝關(guān)大橋主梁動力性能的監(jiān)測參數(shù)。加速度幅值一般指加速度響應(yīng)的某種最大值或某種意義下的有效值。因此,考察加速度峰值和均方根值2種參數(shù)在主梁橫向振動監(jiān)測中的適用性。2014年2月1日傳感器JSD-11-04的橫向加速度幅值與列車速度的相關(guān)性見圖3。可知,隨著列車車速的增加,加速度峰值和均方根值2種參數(shù)均有增大的趨勢,但二者的變化幅度較大,與列車速度的相關(guān)性較差。這主要是由每次行車的輪軌接觸狀態(tài)不同引起的。上述分析結(jié)果表明,加速度幅值與列車速度之間難以建立有效的相關(guān)性模型用以表征高速列車通過時的主梁橫向振動狀態(tài)。此外,高速列車通過時列車-軌道-橋梁系統(tǒng)橫向振動具有顯著的隨機(jī)性,主梁橫向加速度幅值的波動性較大,難以根據(jù)某一測點(diǎn)的橫向加速度響應(yīng)表征橋梁的橫向動力性能。

    圖2 高速列車通過時橫向加速度時程曲線

    圖3 橫向加速度幅值與列車速度的相關(guān)性

    2014年2月1日不同測點(diǎn)之間橫向加速度幅值的互相關(guān)性見圖4??芍?,采用均方根值表征的加速度幅值的互相關(guān)性要明顯優(yōu)于加速度峰值。根據(jù)圖4(b)建立的二次多項(xiàng)式擬合公式分析第2跨跨中加速度均方根的實(shí)測值與擬合值的相關(guān)性,見圖5(a)。其中,第2跨跨中加速度均方根的擬合值由第1跨跨中加速度均方根值代入二次多項(xiàng)式擬合公式計(jì)算得到,實(shí)測值與擬合值的相關(guān)系數(shù)R為0.939 1。同理根據(jù)2014年310 d每天橫向振動數(shù)據(jù)計(jì)算得到的相關(guān)系數(shù),見圖5(b),均在0.90以上。上述分析結(jié)果表明,高速列車通過時主梁2主跨跨中的橫向加速度響應(yīng)的均方根值具有良好的互相關(guān)性。這是因?yàn)橹髁簷M向加速度的均方根值表征了列車-軌道-橋梁系統(tǒng)橫向振動的能量,不同測點(diǎn)之間均方根值的互相關(guān)性可以定量表征橫向振動能量的空間分布特性。

    圖4 不同測點(diǎn)之間橫向加速度幅值的互相關(guān)性

    圖5 加速度均方根值互相關(guān)性的二次多項(xiàng)式擬合效果

    2.3基于橫向加速度監(jiān)測的大橋動力性能異常預(yù)警方法

    根據(jù)2.2節(jié)的分析結(jié)果,橫向加速度的均方根值可以作為反映橋梁動力性能的監(jiān)測參數(shù),并且主梁不同測點(diǎn)之間均方根值的互相關(guān)性模型可以作為橋梁動力性能評價模型。在此基礎(chǔ)上討論大勝關(guān)大橋動力性能異常預(yù)警方法。首先分析均方根值互相關(guān)性模型中二次多項(xiàng)式擬合系數(shù)的長期監(jiān)測結(jié)果。主梁2主跨跨中加速度均方根值互相關(guān)性的二次多項(xiàng)式擬合系數(shù)見圖6。從圖6可以看出,雖然每天的橫向加速度監(jiān)測數(shù)據(jù)可以擬合相關(guān)性良好的二次多項(xiàng)式模型,但是每天的二次多項(xiàng)式擬合系數(shù)仍存在較大的波動。這是因?yàn)?,在溫度等環(huán)境條件的變化影響下列車-軌道-橋梁系統(tǒng)的橫向振動性能會在較寬的范圍內(nèi)波動,直接根據(jù)二次多項(xiàng)式模型的3種擬合系數(shù)難以反映橫向動力性能的異常變化。由于環(huán)境因素的多樣性以及對橋梁振動影響的復(fù)雜性,目前難以全面準(zhǔn)確地測量各種環(huán)境因素,無法有效地分離出各因素的影響量[4-7]。因此,本文針對二次多項(xiàng)式模型的3種擬合系數(shù)研究了基于主成分分析的“環(huán)境條件歸一化”方法。主成分分析方法理論上不需要監(jiān)測環(huán)境因素變量,具有較好的工程應(yīng)用前景[6-7]。

    圖6 加速度均方根值互相關(guān)性的二次多項(xiàng)式擬合系數(shù)

    主成分分析方法是一種將多個相關(guān)變量簡化為幾個不相關(guān)變量線性組合的多元統(tǒng)計(jì)分析方法[8]。它可以在保證數(shù)據(jù)信息丟失最少的情況下,實(shí)現(xiàn)對高維變量空間的降維處理,從而揭示數(shù)據(jù)中所含的關(guān)鍵信息。該方法特別適合于在多變量因子分析中探索系統(tǒng)中某些不可觀測的潛在特征。PCA方法的基本原理是對一個n維向量u=[u1u2…un]T進(jìn)行降維變換,變換矩陣為A,并得到原模式向量v=[v1v2…vm]T(m<n),即

    式中,m個分量v1,v2,…,vm稱為主成分。

    所謂主成分分析,就是找到可以解釋總體方差80%~90%以上的m個分量,組成新的主成分向量v。v中包含有原模式向量中大部分信息,實(shí)現(xiàn)了原數(shù)據(jù)的信息凝聚。對于本文而言,對二次多項(xiàng)式模型的3種擬合系數(shù)組成的3維向量進(jìn)行降維變換得到1維向量。計(jì)算結(jié)果表明,當(dāng)樣本個數(shù)足夠時第1個主成分分量的累積方差可以達(dá)到90%以上。

    根據(jù)上述PCA方法的基本原理,大勝關(guān)大橋動力性能異常預(yù)警方法的主要步驟如下:

    1)正常列車走行狀態(tài)下計(jì)算每天橫向加速度均方根值的二次多項(xiàng)式擬合系數(shù),將二次多項(xiàng)式模型的3種擬合系數(shù)組成原模式向量u。

    2)將k天采集的原模式向量u采用主成分分析進(jìn)行降維變換,變換矩陣為A,得到主成分向量v。變換矩陣A可以較好地反映實(shí)際環(huán)境影響特征,從而有效地消除環(huán)境因素的影響,變換矩陣A的計(jì)算方法見文獻(xiàn)[8]。

    3)根據(jù)主成分向量v可以得到橫向加速度均方根值的重構(gòu)值,定義動力性能異常預(yù)警指標(biāo)D為

    式中:vm(k)為第k天橫向加速度均方根值的實(shí)測日平均值;ve(k)為主成分重構(gòu)后的第k天橫向加速度均方根值的日平均值。當(dāng)監(jiān)測天數(shù)k足夠時,異常預(yù)警指標(biāo)D趨近于0。

    4)將未知列車走行狀態(tài)下的橫向加速度均方根值采用上述正常狀態(tài)下的變換矩陣A進(jìn)行主成分重構(gòu),并計(jì)算異常預(yù)警指標(biāo)D。比較正常列車走行狀態(tài)和未知列車走行狀態(tài)下異常預(yù)警指標(biāo)D是否存在顯著性差異,即可判斷主梁橫向振動是否發(fā)生異常而及時作出預(yù)警。

    以2014年310 d中前260 d采集數(shù)據(jù)作為訓(xùn)練樣本,以后50 d采集數(shù)據(jù)作為檢驗(yàn)樣本。基于PCA方法的異常預(yù)警指標(biāo)訓(xùn)練效果見圖7??芍?dāng)訓(xùn)練樣本天數(shù)較少時,D在-0.08~-0.04波動;隨著監(jiān)測天數(shù)的增多,D逐漸趨近于0;當(dāng)監(jiān)測天數(shù)達(dá)到160 d后D保持穩(wěn)定。在此基礎(chǔ)上將2014年310 d中后50 d的監(jiān)測數(shù)據(jù)進(jìn)行主成分重構(gòu),需要指出的是PCA重構(gòu)時采用的變換矩陣A是根據(jù)前260 d采集數(shù)據(jù)的PCA分析得到的。160~310 d的異常預(yù)警指標(biāo)曲線見圖8。可以看出,無論是訓(xùn)練樣本或是檢驗(yàn)樣本,D均在-0.002~0.002波動,說明PCA方法可以較好地考慮環(huán)境因素的影響,具有較強(qiáng)的重構(gòu)和預(yù)測能力。

    在此基礎(chǔ)上研究本文方法應(yīng)用于大橋動力性能異常預(yù)警的效果,對2014后50 d的主梁第2跨跨中橫向加速度均方根值放大1.05倍,模擬結(jié)構(gòu)動力性能退化對橫向加速度均方根值的影響。主梁第2跨跨中橫向加速度均方根值日平均值的實(shí)測曲線和模擬退化曲線見圖9??芍?,直接由橫向加速度均方根值難以反映出橫向加速度的異常變化,采用異常預(yù)警指標(biāo)可以敏感地發(fā)現(xiàn)橫向加速度的異常變化,從而對大橋動力性能及時作出異常預(yù)警。

    圖7 基于PCA方法的異常預(yù)警指標(biāo)訓(xùn)練效果

    圖8 基于PCA方法的異常預(yù)警指標(biāo)曲線

    圖9 橫向加速度均方根值的日平均值實(shí)測曲線和模擬退化曲線

    3 結(jié)論

    1)高速列車通過時大勝關(guān)大橋主梁的橫向振動具有顯著的隨機(jī)性,根據(jù)主梁某一測點(diǎn)的橫向加速度響應(yīng)及其與列車速度之間的相關(guān)性均難以表征橋梁橫向動力性能。

    2)高速列車通過時主梁2主跨跨中橫向加速度響應(yīng)的均方根值具有良好的互相關(guān)性,通過主梁橫向振動能量的空間分布特性可以表征高速列車通過時橋梁的橫向動力性能。

    3)采用PCA方法可以有效地消除環(huán)境因素對實(shí)測橫向加速度均方根值互相關(guān)性的影響,在此基礎(chǔ)上建立的異常預(yù)警指標(biāo)能夠識別出橫向加速度均方根值的異常變化,可以應(yīng)用于大勝關(guān)大橋動力性能的實(shí)時監(jiān)測與預(yù)警。

    [1]ZHANG N,XIA H,GUO W W.Vehicle-Bridge Interaction Analysis Under High-Speed Trains[J].Journal of Sound and Vibration,2008,309(3/5):407-425.

    [2]岳青,吳來義,朱利明,等.高速鐵路大跨度橋梁運(yùn)營監(jiān)測系統(tǒng)研究[J].鐵道建筑,2015(12):1-6.

    [3]曾慶元,向俊,婁平,等.列車脫軌的力學(xué)機(jī)理與防止脫軌理論[J].鐵道科學(xué)與工程學(xué)報,2004,1(1):19-31.

    [4]HUA X G,NI Y Q,KO J M,et al.Modeling of Temperature-Frequency Correlation Using Combined Principal Component Analysis and Support Vector Regression Technique[J].Journal of Computing in Civil Engineering,2007,21(2):122-135.

    [5]ZHOU H F,Ni Y Q,KO J M.Modeling of Wind and Temperature Effects on Modal Frequencies and Analysis of Relative Strength of Effect[J].Wind and Structures,2008,11(1):35-50.

    [6]DERAEMAEKER A,REYNDERS E,ROECK G D,et al.Vibration-Based StructuralHealthMonitoringUsingOutput-only Measurements Under Changing Environment[J].Mechanical Systems and Signal Processing,2008,22(1):34-56.

    [7]丁幼亮,李愛群,耿方方.考慮環(huán)境因素影響的懸索橋整體狀態(tài)預(yù)警方法研究[J].東南大學(xué)學(xué)報(自然科學(xué)版),2010,40(5):1052-1056.

    [8]梅長林,周家良.實(shí)用統(tǒng)計(jì)方法[M].北京:科學(xué)出版社,2002.

    Early Warning Method of Abnormal Dynamic Performance of High Speed Railway Bridge Based on Transverse Acceleration Monitoring

    GENG Fangfang1,2,DING Youliang1,3,ZHAO Hanwei1,WU Laiyi4,YUE Qing4
    (1.Key Laboratory of Concrete and Prestressed Concrete Structures of Ministry of Education,Southeast University,Nanjing Jiangsu 210096,China;2.Nanjing Institute of Technology,Nanjing Jiangsu 211167,China;3.Key Laboratory for Health Monitoring and Control of Large Structures,Shijiazhuang Hebei 050043,China;4.China Railway Major Bridge(Nanjing)Bridge and Tunnel Inspect&Retrofit Co.,Ltd.,Nanjing Jiangsu 210032,China)

    M aking use of transverse acceleration monitoring data of the main girder of Dashengguan bridge,the monitoring and warning method of abnormal dynamic performance was formulated based on principal component analysis.T he root mean square value of transverse acceleration of the girder was used as the monitoring parameter reflecting the dynamic performance.T he correlation model of root mean square values was established as the assessment model of the dynamic performance.T he influences of environmental conditions on the assessment model were eliminated by using the principal component analysis,the corresponding warning index was proposed.Analysis results show that the correlation of root mean square values measured at different points can be expressed as quadratic polynomials.And the proposed waning method provides a good capability for detecting the abnormal changes of transverse accelerations of the maingirder.It is suitable for real-time warning of the dynamic performance.

    High speedrailwaybridge;Vibrationmonitoring;Principal componentanalysis;Environmental condition;Real-time warning

    U448.13

    ADOI:10.3969/j.issn.1003-1995.2016.09.01

    1003-1995(2016)09-0001-05

    (責(zé)任審編鄭冰)

    2016-03-31;

    2016-08-01

    國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)(2015CB060000);國家自然科學(xué)基金(51438002,51578138);河北省大型結(jié)構(gòu)健康診斷與控制重點(diǎn)實(shí)驗(yàn)室開放課題(201502)

    耿方方(1983—),女,講師,博士。

    丁幼亮(1979—),男,研究員,博士。

    猜你喜歡
    方根值主梁大橋
    磁流變彈性體減振單元動力學(xué)分析
    龍門起重機(jī)主梁腹板拓?fù)鋬?yōu)化
    一圖看懂藏木大橋
    客車平順性仿真及優(yōu)化
    數(shù)鴨子
    上天入地——港珠澳大橋
    超寬組合梁斜拉橋主梁剪力滯效應(yīng)研究
    上海公路(2017年4期)2017-03-01 07:04:16
    斜拉橋主梁施工牽鎖掛籃結(jié)構(gòu)設(shè)計(jì)
    1000噸龍門吊主梁跨塢安裝
    Studyof velocityfluctuations in the plenum of a 3/4open jet automotive wind tunnel
    精品午夜福利在线看| 亚洲av国产av综合av卡| 欧美黑人精品巨大| www.自偷自拍.com| a 毛片基地| 在现免费观看毛片| 国产免费一区二区三区四区乱码| 国产男人的电影天堂91| 又大又爽又粗| 亚洲国产精品一区二区三区在线| 欧美人与善性xxx| 久久久欧美国产精品| 国产亚洲av片在线观看秒播厂| 久久久久国产精品人妻一区二区| 久久韩国三级中文字幕| 国产片内射在线| 国产精品三级大全| 久久久久视频综合| 国产探花极品一区二区| 久久人人97超碰香蕉20202| 少妇 在线观看| 美女国产高潮福利片在线看| 大片免费播放器 马上看| 黄片播放在线免费| 久久久久久久精品精品| 欧美国产精品一级二级三级| 美女扒开内裤让男人捅视频| 亚洲国产欧美在线一区| 街头女战士在线观看网站| 最近中文字幕2019免费版| 午夜福利免费观看在线| 波野结衣二区三区在线| 欧美激情 高清一区二区三区| 久热这里只有精品99| 欧美 日韩 精品 国产| 国产亚洲一区二区精品| 亚洲精品第二区| 人妻人人澡人人爽人人| 人人澡人人妻人| 国产av国产精品国产| 啦啦啦在线免费观看视频4| 亚洲av电影在线观看一区二区三区| 免费观看a级毛片全部| 国产高清不卡午夜福利| 男女午夜视频在线观看| 亚洲国产成人一精品久久久| 欧美成人精品欧美一级黄| 晚上一个人看的免费电影| 国产精品久久久久久精品古装| 天堂俺去俺来也www色官网| 久久久久久免费高清国产稀缺| 国产亚洲一区二区精品| 国产精品久久久久久精品古装| 亚洲精品国产av蜜桃| 免费看av在线观看网站| 青春草国产在线视频| 亚洲三区欧美一区| 一边摸一边做爽爽视频免费| 久久久久国产精品人妻一区二区| 亚洲五月色婷婷综合| 国产亚洲av高清不卡| 麻豆精品久久久久久蜜桃| 妹子高潮喷水视频| 九色亚洲精品在线播放| 亚洲欧美中文字幕日韩二区| 一级a爱视频在线免费观看| 久热爱精品视频在线9| 国产精品国产三级国产专区5o| 美女视频免费永久观看网站| 三上悠亚av全集在线观看| 亚洲精品中文字幕在线视频| 伊人久久国产一区二区| 精品国产一区二区三区四区第35| 欧美黑人精品巨大| 精品人妻熟女毛片av久久网站| 国产毛片在线视频| 成人免费观看视频高清| 满18在线观看网站| 久久99精品国语久久久| 亚洲av欧美aⅴ国产| 亚洲精品,欧美精品| 亚洲伊人久久精品综合| 国产成人欧美在线观看 | 午夜福利视频精品| 99久久99久久久精品蜜桃| 亚洲四区av| 亚洲欧美激情在线| 日韩中文字幕欧美一区二区 | 国产精品国产三级国产专区5o| 热re99久久精品国产66热6| 999久久久国产精品视频| 涩涩av久久男人的天堂| 欧美97在线视频| 1024香蕉在线观看| 日韩,欧美,国产一区二区三区| 国产一区二区三区综合在线观看| 老汉色av国产亚洲站长工具| 天天添夜夜摸| 波野结衣二区三区在线| 国产欧美亚洲国产| 亚洲精品美女久久av网站| 亚洲成色77777| 免费女性裸体啪啪无遮挡网站| 国产免费福利视频在线观看| 久久久久久久久久久久大奶| 亚洲av欧美aⅴ国产| av网站免费在线观看视频| 婷婷色综合大香蕉| 性少妇av在线| 亚洲av男天堂| 亚洲一级一片aⅴ在线观看| 熟女少妇亚洲综合色aaa.| 搡老乐熟女国产| 国产1区2区3区精品| 一本大道久久a久久精品| 美女国产高潮福利片在线看| 精品少妇一区二区三区视频日本电影 | 国产一区二区三区综合在线观看| 免费高清在线观看日韩| 精品国产露脸久久av麻豆| 国产熟女欧美一区二区| 满18在线观看网站| 免费人妻精品一区二区三区视频| 只有这里有精品99| bbb黄色大片| 中文精品一卡2卡3卡4更新| 熟女少妇亚洲综合色aaa.| 欧美精品av麻豆av| 久久韩国三级中文字幕| 亚洲欧美中文字幕日韩二区| 精品国产超薄肉色丝袜足j| 晚上一个人看的免费电影| 国产精品久久久久久精品古装| 嫩草影视91久久| 天天影视国产精品| 欧美变态另类bdsm刘玥| 久久久久精品人妻al黑| 久久精品人人爽人人爽视色| 精品免费久久久久久久清纯 | 久久久久久久久免费视频了| 日日爽夜夜爽网站| 午夜福利影视在线免费观看| 亚洲精品国产av成人精品| 亚洲精品日本国产第一区| 亚洲精品成人av观看孕妇| 精品人妻一区二区三区麻豆| 一个人免费看片子| 精品人妻熟女毛片av久久网站| 精品国产一区二区三区久久久樱花| 国产成人精品久久二区二区91 | 日韩一卡2卡3卡4卡2021年| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av在线观看美女高潮| 国产爽快片一区二区三区| 亚洲一区中文字幕在线| 天堂8中文在线网| 久久青草综合色| 国产成人a∨麻豆精品| 一边摸一边做爽爽视频免费| 热re99久久精品国产66热6| 天堂俺去俺来也www色官网| 婷婷色av中文字幕| 亚洲第一av免费看| 国产淫语在线视频| av视频免费观看在线观看| 国产成人午夜福利电影在线观看| 国产一卡二卡三卡精品 | 精品国产露脸久久av麻豆| 亚洲av电影在线进入| 亚洲,欧美,日韩| 国产不卡av网站在线观看| 18禁国产床啪视频网站| 欧美人与善性xxx| 亚洲成色77777| 精品国产乱码久久久久久男人| 1024视频免费在线观看| 国产成人91sexporn| 国产成人91sexporn| 免费av中文字幕在线| √禁漫天堂资源中文www| svipshipincom国产片| 久久热在线av| 日本av手机在线免费观看| 国语对白做爰xxxⅹ性视频网站| 韩国av在线不卡| 男女午夜视频在线观看| 久久国产精品男人的天堂亚洲| 男女无遮挡免费网站观看| 精品一品国产午夜福利视频| a级片在线免费高清观看视频| 中文字幕人妻丝袜一区二区 | 国产亚洲最大av| 美女扒开内裤让男人捅视频| 成人影院久久| 亚洲精品久久午夜乱码| 亚洲欧美中文字幕日韩二区| 韩国av在线不卡| 国产精品蜜桃在线观看| 欧美成人精品欧美一级黄| 国产亚洲av片在线观看秒播厂| 亚洲精品久久成人aⅴ小说| 亚洲国产精品999| 欧美日韩亚洲综合一区二区三区_| 亚洲精品一二三| 高清视频免费观看一区二区| 美女高潮到喷水免费观看| 国产女主播在线喷水免费视频网站| 在线 av 中文字幕| 亚洲综合精品二区| 久久久精品区二区三区| 欧美97在线视频| 深夜精品福利| 性高湖久久久久久久久免费观看| 制服丝袜香蕉在线| 亚洲精品,欧美精品| 国产欧美日韩综合在线一区二区| 国产黄频视频在线观看| 一区二区三区精品91| 亚洲国产av新网站| 中文精品一卡2卡3卡4更新| 久久热在线av| 啦啦啦视频在线资源免费观看| av国产精品久久久久影院| 午夜福利视频精品| 少妇人妻精品综合一区二区| 免费黄频网站在线观看国产| 蜜桃在线观看..| 精品国产乱码久久久久久男人| 欧美精品人与动牲交sv欧美| 母亲3免费完整高清在线观看| 欧美日韩视频精品一区| 亚洲av国产av综合av卡| 狂野欧美激情性xxxx| 亚洲精品久久午夜乱码| 欧美国产精品va在线观看不卡| 最近2019中文字幕mv第一页| 亚洲欧美一区二区三区国产| 热re99久久国产66热| 又粗又硬又长又爽又黄的视频| 操美女的视频在线观看| 欧美国产精品一级二级三级| 久久久久网色| 嫩草影视91久久| 精品国产一区二区三区四区第35| 国产成人精品在线电影| 满18在线观看网站| 午夜av观看不卡| 熟女少妇亚洲综合色aaa.| 丰满迷人的少妇在线观看| 老司机影院成人| 男女无遮挡免费网站观看| 十八禁网站网址无遮挡| 99精国产麻豆久久婷婷| av网站免费在线观看视频| 日韩制服骚丝袜av| 欧美另类一区| av国产久精品久网站免费入址| 成年av动漫网址| 另类亚洲欧美激情| videos熟女内射| 青春草视频在线免费观看| 1024香蕉在线观看| 男的添女的下面高潮视频| 哪个播放器可以免费观看大片| 欧美日韩亚洲国产一区二区在线观看 | 一级a爱视频在线免费观看| 十八禁高潮呻吟视频| 少妇被粗大的猛进出69影院| 精品国产一区二区三区久久久樱花| 国产精品偷伦视频观看了| av一本久久久久| 无限看片的www在线观看| 久久精品久久久久久久性| 久久久久久人妻| 亚洲精品美女久久av网站| 男女免费视频国产| 久久久久精品久久久久真实原创| 啦啦啦在线观看免费高清www| 高清欧美精品videossex| 天天躁夜夜躁狠狠躁躁| 欧美黄色片欧美黄色片| 最新的欧美精品一区二区| 可以免费在线观看a视频的电影网站 | 性少妇av在线| 欧美日韩视频高清一区二区三区二| 一本久久精品| 欧美av亚洲av综合av国产av | 国产精品久久久久久人妻精品电影 | 一二三四在线观看免费中文在| 久久天堂一区二区三区四区| 视频在线观看一区二区三区| 日韩中文字幕视频在线看片| 最近2019中文字幕mv第一页| 国产乱人偷精品视频| 精品视频人人做人人爽| 国产av码专区亚洲av| 国产成人精品久久久久久| 两性夫妻黄色片| 色综合欧美亚洲国产小说| 国精品久久久久久国模美| 丝袜脚勾引网站| 国产成人系列免费观看| 性高湖久久久久久久久免费观看| 国产成人精品无人区| 色婷婷av一区二区三区视频| 国产人伦9x9x在线观看| 欧美乱码精品一区二区三区| 麻豆精品久久久久久蜜桃| 欧美变态另类bdsm刘玥| 久久99精品国语久久久| 成人国产av品久久久| 精品亚洲成a人片在线观看| 一本一本久久a久久精品综合妖精| 一级毛片黄色毛片免费观看视频| 性少妇av在线| 丰满迷人的少妇在线观看| 欧美精品亚洲一区二区| 亚洲图色成人| 黄色一级大片看看| 午夜日本视频在线| 亚洲成国产人片在线观看| 91精品国产国语对白视频| 国产欧美日韩一区二区三区在线| 久久天堂一区二区三区四区| 国产福利在线免费观看视频| 操美女的视频在线观看| 日本欧美国产在线视频| 久久久久久久大尺度免费视频| 精品一品国产午夜福利视频| 午夜福利免费观看在线| 七月丁香在线播放| 91精品国产国语对白视频| 另类精品久久| 在线天堂中文资源库| 男的添女的下面高潮视频| 日韩不卡一区二区三区视频在线| 亚洲国产精品一区三区| 亚洲国产精品成人久久小说| 日本91视频免费播放| 久久国产精品大桥未久av| 人妻人人澡人人爽人人| 自线自在国产av| 成年人免费黄色播放视频| 99国产综合亚洲精品| 精品福利永久在线观看| 七月丁香在线播放| 欧美国产精品va在线观看不卡| 精品人妻一区二区三区麻豆| videosex国产| 久久久久国产一级毛片高清牌| 99九九在线精品视频| 狂野欧美激情性bbbbbb| 欧美日本中文国产一区发布| 日韩大码丰满熟妇| 久久热在线av| 一个人免费看片子| 久久久久精品人妻al黑| 波多野结衣一区麻豆| 我要看黄色一级片免费的| 巨乳人妻的诱惑在线观看| 丝袜美腿诱惑在线| 欧美少妇被猛烈插入视频| 超碰97精品在线观看| 欧美变态另类bdsm刘玥| 免费观看a级毛片全部| 久久久久久久精品精品| 国产精品蜜桃在线观看| 精品国产乱码久久久久久男人| av片东京热男人的天堂| 一本久久精品| 在线观看一区二区三区激情| 最近中文字幕2019免费版| 一本—道久久a久久精品蜜桃钙片| 看免费av毛片| 丝瓜视频免费看黄片| 日本黄色日本黄色录像| 久久久精品国产亚洲av高清涩受| 久久久久久久久免费视频了| 精品少妇一区二区三区视频日本电影 | 欧美国产精品va在线观看不卡| 18禁观看日本| 国产精品.久久久| 女人精品久久久久毛片| 啦啦啦在线免费观看视频4| 国产精品一二三区在线看| 亚洲第一青青草原| 操出白浆在线播放| 亚洲精品av麻豆狂野| 国产精品久久久久久久久免| 久久精品国产亚洲av涩爱| 亚洲国产av新网站| 亚洲成人手机| 国产乱人偷精品视频| 1024视频免费在线观看| 伦理电影免费视频| 侵犯人妻中文字幕一二三四区| 精品一区二区三卡| 国产精品偷伦视频观看了| 两个人看的免费小视频| 9热在线视频观看99| 另类精品久久| 亚洲国产精品国产精品| 国产成人欧美| 999久久久国产精品视频| 超碰97精品在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 日本91视频免费播放| 免费黄网站久久成人精品| 中文乱码字字幕精品一区二区三区| 人人妻人人添人人爽欧美一区卜| 亚洲欧美激情在线| 色94色欧美一区二区| 七月丁香在线播放| 国产老妇伦熟女老妇高清| 精品人妻在线不人妻| 人人妻人人澡人人看| 99久国产av精品国产电影| 人妻 亚洲 视频| 中文字幕人妻熟女乱码| av卡一久久| 欧美在线黄色| 亚洲七黄色美女视频| 午夜激情久久久久久久| 女性被躁到高潮视频| 精品国产一区二区三区四区第35| 成人黄色视频免费在线看| 热re99久久国产66热| 美女午夜性视频免费| 亚洲人成电影观看| av在线播放精品| 女人高潮潮喷娇喘18禁视频| 婷婷色麻豆天堂久久| av在线播放精品| 精品人妻一区二区三区麻豆| 1024视频免费在线观看| 看非洲黑人一级黄片| 精品久久久精品久久久| 免费黄网站久久成人精品| 精品免费久久久久久久清纯 | 国产精品熟女久久久久浪| 午夜精品国产一区二区电影| av.在线天堂| www.精华液| 三上悠亚av全集在线观看| 91成人精品电影| 婷婷色麻豆天堂久久| 午夜免费鲁丝| 精品人妻一区二区三区麻豆| 久久精品亚洲av国产电影网| 中文字幕高清在线视频| 久久久久久人人人人人| 亚洲精品视频女| 看免费av毛片| 极品少妇高潮喷水抽搐| 成人手机av| 亚洲 欧美一区二区三区| 精品亚洲乱码少妇综合久久| 男女之事视频高清在线观看 | 欧美成人精品欧美一级黄| 一级片免费观看大全| 国产成人av激情在线播放| 伊人久久国产一区二区| 久久性视频一级片| 久久97久久精品| avwww免费| 七月丁香在线播放| 亚洲欧美日韩另类电影网站| 看免费成人av毛片| 多毛熟女@视频| 丝袜在线中文字幕| 看十八女毛片水多多多| 高清av免费在线| 一边摸一边做爽爽视频免费| 欧美激情极品国产一区二区三区| 日韩一卡2卡3卡4卡2021年| 亚洲成人一二三区av| 看免费成人av毛片| 最近最新中文字幕免费大全7| 亚洲成人一二三区av| 欧美日韩亚洲综合一区二区三区_| 777米奇影视久久| 一级黄片播放器| 久久久久久久久久久久大奶| 亚洲精品国产色婷婷电影| 无限看片的www在线观看| 亚洲国产日韩一区二区| 99热网站在线观看| 老司机亚洲免费影院| 国产色婷婷99| 免费在线观看黄色视频的| 人妻人人澡人人爽人人| 在线观看免费视频网站a站| 亚洲精品,欧美精品| 一边摸一边抽搐一进一出视频| 尾随美女入室| 中国三级夫妇交换| 久久久久久久国产电影| 制服丝袜香蕉在线| av国产精品久久久久影院| 国产伦理片在线播放av一区| 国产麻豆69| 日韩一区二区三区影片| 电影成人av| 日日爽夜夜爽网站| 国产免费又黄又爽又色| 色吧在线观看| 亚洲精品成人av观看孕妇| 色94色欧美一区二区| 校园人妻丝袜中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 国产极品天堂在线| 国产精品久久久久久精品电影小说| 久久鲁丝午夜福利片| 日韩伦理黄色片| 激情视频va一区二区三区| 久久国产精品大桥未久av| 天堂8中文在线网| 巨乳人妻的诱惑在线观看| 亚洲国产欧美一区二区综合| 久久国产亚洲av麻豆专区| 久久久久精品人妻al黑| 欧美日韩亚洲国产一区二区在线观看 | 十分钟在线观看高清视频www| 九色亚洲精品在线播放| 国产欧美亚洲国产| 久久这里只有精品19| 免费看av在线观看网站| 秋霞伦理黄片| 国产欧美日韩综合在线一区二区| 亚洲精品在线美女| 叶爱在线成人免费视频播放| 最近中文字幕高清免费大全6| 国产伦理片在线播放av一区| 色吧在线观看| 欧美亚洲日本最大视频资源| 视频在线观看一区二区三区| 国产在视频线精品| 亚洲国产日韩一区二区| 看免费av毛片| 两个人免费观看高清视频| av在线观看视频网站免费| 97精品久久久久久久久久精品| videosex国产| www.精华液| 亚洲第一av免费看| 91国产中文字幕| 女人爽到高潮嗷嗷叫在线视频| 亚洲人成77777在线视频| 国产精品女同一区二区软件| 高清在线视频一区二区三区| 日韩精品有码人妻一区| 99久国产av精品国产电影| videosex国产| 国产1区2区3区精品| 人妻 亚洲 视频| 免费不卡黄色视频| 亚洲国产精品国产精品| 91aial.com中文字幕在线观看| 天天躁日日躁夜夜躁夜夜| 国产精品av久久久久免费| 我的亚洲天堂| 国产日韩欧美在线精品| 免费观看a级毛片全部| 97人妻天天添夜夜摸| 美女大奶头黄色视频| 啦啦啦视频在线资源免费观看| 人妻 亚洲 视频| 伦理电影免费视频| 亚洲成人免费av在线播放| 99精品久久久久人妻精品| 性色av一级| 另类亚洲欧美激情| 国产一区亚洲一区在线观看| 男人舔女人的私密视频| 人人妻人人澡人人看| 亚洲三区欧美一区| 日韩一本色道免费dvd| 亚洲欧美日韩另类电影网站| 中文字幕人妻丝袜制服| 久久 成人 亚洲| 看十八女毛片水多多多| 青春草亚洲视频在线观看| 考比视频在线观看| 制服人妻中文乱码| 美女午夜性视频免费| 一区二区av电影网| 亚洲精品在线美女| 久久免费观看电影| 99热全是精品| 免费女性裸体啪啪无遮挡网站| 亚洲,一卡二卡三卡| 欧美人与善性xxx| 美女国产高潮福利片在线看| 波野结衣二区三区在线| av电影中文网址| 亚洲国产日韩一区二区| 久久午夜综合久久蜜桃| 1024香蕉在线观看| 亚洲精品国产一区二区精华液| 亚洲人成网站在线观看播放| 大码成人一级视频| tube8黄色片| 熟女少妇亚洲综合色aaa.| 欧美在线黄色| 色婷婷av一区二区三区视频| 国产精品久久久久久久久免| 只有这里有精品99| 欧美 日韩 精品 国产| 丝袜美腿诱惑在线| 免费久久久久久久精品成人欧美视频| 十八禁高潮呻吟视频| 美女主播在线视频| 婷婷色av中文字幕| 久久精品久久久久久噜噜老黄| 99热全是精品|