999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

近斷層脈沖效應對大跨度結合梁斜拉橋地震響應的影響

2019-11-13 11:58:52侯烈張龍奇師新虎
鐵道科學與工程學報 2019年10期
關鍵詞:橋梁

侯烈,張龍奇,師新虎

近斷層脈沖效應對大跨度結合梁斜拉橋地震響應的影響

侯烈1,張龍奇1,師新虎2

(1. 四川交通職業技術學院,四川 成都 611130;2. 西南交通大學 土木工程學院,四川 成都 610031)

為了研究近斷層脈沖效應對大跨度結合梁斜拉橋地震響應的影響,明確近斷層脈沖地震動長周期脈沖和高頻分量對大跨度結合梁斜拉橋地震響應的貢獻。首先采用基于小波分析的近斷層脈沖地震動識別方法從集集地震中選取相同場地下的脈沖/非脈沖2類近斷層地震記錄,并對脈沖型地震動的長周期速度脈沖進行提取,研究長周期、高幅值速度脈沖對近斷層脈沖地震動反應譜的放大效應。最后以某一大跨度結合梁斜拉橋為例,研究近斷層速度脈沖效應對其地震響應的影響。研究結果表明:由于大跨度結合梁斜拉橋結構體系較柔,自振周期較長,近斷層脈沖效應對其地震響應的影響顯著,與非脈沖地震作用相比,近斷層脈沖效應使得主梁和橋塔的地震響應都顯著增大,最大改變率可達177.9%;主梁和橋塔的地震位移響應主要由長周期脈沖所決定,高頻成分的影響相對較小,但是對橋塔內力響應而言,其受近斷層地震動高頻成分的影響較大,在基于內力控制的結構抗震設計中應該充分考慮長周期脈沖和高頻成分2部分的影響。

近斷層脈沖;結合梁斜拉橋;地震響應;小波分析

大跨度結合梁斜拉橋由于其造型美觀,受力合理,能夠充分發揮鋼材和混凝土2種材料的性能優勢,因此在許多大型的跨江、跨海工程中都得到了廣泛應用[1]。上世紀建成的南浦大橋、楊浦大橋、徐浦大橋以及后來的武漢二七長江大橋等一大批橋梁都采用了這種方案。對于這類大跨度結合梁斜拉橋而言,其跨度較大,結構體系較柔,自振周期較長,對地震動的長周期成分較為敏感[2]。前場方向性效應(Forward Directivity effect)和滑沖效應(Fling Step effect)引起的2類近斷層脈沖地震動是最為典型的長周期地震動,其中除了含有豐富的高頻成分外,還包含有顯著的長周期脈沖[3]。近30年來發生的幾次大地震,美國Northridge地震、日本Kobe地震、臺灣集集地震以及我國2008年的汶川地震等多次地震都集中在淺源近場區[4],這類近斷層長周期地震動對橋梁結構的不利影響已得到充分體現,國內外很多專家學者都對其進行了深入的研究。陳令坤等[5]研究了近斷層脈沖效應對多跨簡支梁橋非線性地震響應的影響,研究結果表明:與遠場地震作用相比,近斷層方向脈沖效應會顯著增大橋梁結構的變形;夏修身等[6]以斷層距作為近斷層地震動的識別參數,研究了近斷層地震動對高墩鐵路橋梁高墩搖擺反應的影響,研究結果表明近斷層地震動會顯著增大高墩墩頂的搖擺位移。孫穎等[7]以某五跨一聯的隔震連續梁橋為例,研究具有明顯脈沖的近、遠場長周期地震動作用下鉛芯橡膠支座(LRB)和摩擦擺(FPB)的減震效果,研究結果表明長周期地震作用下隔震支座的隔震效果會嚴重下降。目前關于近斷層長周期地震動對結構抗震性能影響的研究主要是以斷層距作為選擇近/遠場地震記錄的依據,進而通過結構的地震響應分析研究脈沖地震動對結構抗震性能的影響。近斷層脈沖地震動的長周期和高頻分量對橋梁地震響應的影響并不明確。此外,目前的研究對象主要集中在簡支梁和連續梁等常規橋梁結構,通常這類結構的剛度較大,自振周期較小,脈沖周期T與結構自振周期的比值T/較大,地震作用下橋梁結構的地震反應主要以一階模態的貢獻為主。對于大跨度結合梁斜拉橋而言,其結構較為輕柔,地震作用下,高階振型的參與較大,自振周期與脈沖周期T更為接近,結構發生共振的可能性更大,近斷層脈沖地震動對結構的不利影響將會更加嚴重。鑒于此,本文以在建的某一大跨度結合梁斜拉橋為例,基于ANSYS建立其三維空間有限元模型,對其進行地震響應分析,研究近斷層脈沖效應對大跨度結合梁斜拉橋地震響應的影響,從而為近斷層地區大跨度結合梁斜拉橋的抗震設計提供參考。

1 近斷層脈沖地震動的選取及頻譜特性分析

1.1 近斷層脈沖地震動的選取

目前,關于近斷層地震動的定義還存在較大爭議,不少學者都采用斷層距作為近斷層地震動的劃分標準[8?9],根據CALTRANS的采用的定義方式,斷層距小于10英里(約15 km)的地震動都為近斷層地震動[8],而Stewart等[9]考慮到近場效應受到震級和場地條件等各種因素的影響,將斷層距20~60 km的區域作為界定近場和遠場的依據。Iervolino等[10]的研究表明并不是所有的近斷層地震記錄都為脈沖地震記錄,很大一部分近斷層地震記錄中并不包含脈沖成分,其特性和遠場地震動并無差別。Baker[11]通過連續小波變換對近斷層脈沖地震動進行識別,從而實現了對近斷層脈沖地震動和非脈沖地震動的有效區分,并得到了普遍的認同,因此,本文采用該方法選擇近斷層脈沖地震記錄和非脈沖地震記錄,從而為近斷層作用下大跨度結合梁斜拉橋的地震響應提供依據。

除了地震動的長周期脈沖外,場地效應也是影響橋梁結構地震響應的一個重要因素。為了對場地效應的影響進行考慮,本文采用30 m土層平均剪切波速30作為場地條件的表征,分別從太平洋地震中心PEER數據庫選擇了相同場地下的脈沖地震記錄TCU036-EW (簡稱TCU036)和非脈沖地震記錄TCU070-EW (簡稱TCU070)作為分析依據。

表1 地震動詳細信息

圖1 TCU036-EW地震記錄

表1分別為TCU036和TCU070 2條地震記錄的詳細信息,圖1和圖2為2條地震動的速度時程和按照連續小波變換提取的脈沖時程。其中表1中的脈沖指標P即為按照Baker[11]所提出的脈沖地震動識別方法計算出的脈沖性指標,當P大于0.85,該地震動即被判定為脈沖型地震動,當P小于0.15即被判定為非脈沖地震記錄。從表1和圖1~2中可以看出,雖然TCU036和TCU070 2條地震動均來自于同一次地震,并且斷層距和場地條件都相同,但是由于地震動測站與斷層相對位置的隨機性,2條地震記錄顯示出完全不同的脈沖性,TCU036速度時程中包含有一個顯著的長周期高幅值脈沖,屬于典型的脈沖型地震動,而根據Baker[11]所提出的脈沖地震動識別方法,其脈沖指標P為0.999 9,按照脈沖地震動分類標準,其也屬于脈沖型地震動,由此也證明了Baker所提出的脈沖地震動識別方法的有效性。

圖2 TCU070-EW地震記錄

1.2 近斷層地震動的反應譜分析

鑒于TCU036和TCU070 2條地震脈沖性的顯著差異,有必要對其進行頻譜特性的分析。圖3(a)和圖3(b)分別為阻尼比為5%時TCU036和TCU070 2條地震記錄的擬加速度反應譜和擬速度反應譜(PGA統一調整到0.3,調幅系數見表1)。

從圖3中可以看出:對于線彈性的單自由度體系而言,當周期小于1 s時非脈沖地震動記錄TCU070的反應譜明顯要高于脈沖型地震動TCU036的譜值,但是當周期大于1 s后,脈沖地震記錄TCU036的反應譜明顯要高于非脈沖地震動TCU070的反應譜,而且在周期為4 s位置處附近,二者的差距達到最大。由此說明了近斷層脈沖效應對地震動的反應譜存在顯著的放大效應,當結構的自振周期恰好落在這一附近時,其地震需求將顯著的增大。

(a) 加速度反應譜;(b) 擬速度反應譜

2 近斷層速度脈沖對單自由度系統(SDOF)位移響應的影響

由于在基于位移性能的橋梁抗震設計中,一般以最大位移作為工程需求參數,因此這里研究近斷層速度脈沖效應對單自由度系統(SDOF)位移響應的影響。該單自由度系統滿足:1) 自振周期位于0.02~10 s,大部分工程結構的自振周期都位于該周期范圍內;2) 延性系數為1,即該系統為線彈性系統;3) 阻尼比=5%。

圖5分別為TCU036,TCU070以及按照連續小波變換[11]從TCU036速度時程中提取的脈沖作用下該單自由度系統的位移響應。從圖5中可以看出,脈沖地震記錄TCU036作用下,在整個周期內單自由度系統的位移響應都隨著結構自振周期的增大而增大。而非脈沖地震記錄TCU070作用下,雖然在周期5 s位置處出現了一個顯著的峰值,但是整體而言,隨著自振周期的增大,單自由度系統的響應也呈現出逐漸增大的趨勢。在TCU036的長周期脈沖時程作用下,隨著結構周期的增大,單自由度系統的響應呈現出現增大后減小的趨勢,并且當結構周期在8 s位置處達到了最大值。這主要是由于該長周期脈沖中高頻成分都已經被過濾,僅含有長周期成分,其頻率成分相對較為單一所導致的。

圖4 單自由度體系計算模型

圖5 單自由度體系位移響應

此外,從圖5中還可以看出,對于長周期結構而言,TCU036地震記錄以及其對應的長周期脈沖作用下,單自由度系統的位移響應十分接近,由此說明了對于長周期結構,其地震響應主要是由速度時程中含有的長周期高幅值脈沖所決定的,高頻成分對其貢獻很小,因此在其抗震性能分析時應該考慮近斷層脈沖效應的影響,否則將嚴重低估結構的地震響應。

3 研究對象與橋梁有限元模型

3.1 工程背景及橋梁結構建模

本文以某大跨度結合梁斜拉橋為例,研究近斷層脈沖效應對大跨度結合梁斜拉橋地震響應的影響。該橋主跨720 m,全長1 355 m,跨徑布置為70 m+93 m+142 m+720 m+240 m+90 m。主梁采用鋼?混結合梁,梁高3.8 m,橋面寬度36.5 m。橋塔采用H型橋塔,塔高238 m,上塔柱高69.33 m材料為C55預應力混凝土,斜拉索采用抗拉強度1 670 MPa高強平行鋼絲索,標準間距12 m,全橋共設116對(232根),索面呈扇形布置,支撐體系為半漂浮體系,整個橋梁結構布置如圖6所示。

單位:m

圖7 結合梁斜拉橋動力計算模型

基于ANSYS有限平臺建立橋梁結構動力計算三維有限元模型,整個橋梁結構包括主梁、拉索、橋塔、邊墩以及輔助墩5個部分。橋塔、邊墩及輔助墩均采用BEAM44單元進行模擬,拉索采用LINK8單元進行模擬,主梁的鋼箱部分采用BEAM44單元模擬,混凝土橋面板采用SHELL181單元進行模擬。主梁與橋塔、邊墩和輔助墩均采用節點耦合的方式設置主從約束,墩底和塔底均采用固結,即不考慮土?結構相互作用。整個橋梁結構空間有限元模型如圖7所示。

3.2 橋梁結構動力特性

基于上述有限元動力計算模型,采用Lanczos方法對該算例橋梁進行動力特性分析,限于文章篇幅,表2僅列出了前5階振型的頻率和對應的振型特征。從表2中可以看出,該橋梁1階模態為全橋縱向漂移,自振周期為9.07 s,符合該類橋型的特點,屬于典型的長周期結構。由于主梁采用的是鋼?混結合梁斜拉橋,截面類似與π型截面,橫向剛度較小,因此第2階振型表現為對稱側彎,且早于了主梁對稱豎彎出現;橋塔的側彎振型從第4階開始出現,并且最早出現的是右側橋塔的側彎,這主要是由于右側只有一個輔助墩,而左側有2個輔助墩,因此左側橋墩產生的約束更強,橋塔的振動頻率也就更高。

表2 橋梁動力特性

4 大跨度結合梁斜拉橋地震響應分析

為了研究近斷層脈沖效應對大跨度結合梁斜拉橋地震響應的影響,將TCU036和TCU070 2條地震記錄以及按照連續小波變化提取的TCU036長周期脈沖分別乘以調幅系數,使其峰值加速度達到0.3,然后將其按照縱向輸入計算模型,通過時程分析獲得橋梁結構的地震響應。

4.1 橋塔地震響應

圖8~9分別為TCU036和TCU070以及TCU036長周期脈沖作用下該橋主塔縱向位移和剪力的峰值響應包絡曲線。從圖8和圖9中可以看出:

圖8 橋塔縱向位移響應峰值

圖9 橋塔縱向剪力響應峰值

1) 隨著橋塔高墩的增加,橋塔位移呈逐漸增大的趨勢,在塔頂位置處,位移達到最大,而主塔內力總體呈逐漸下降的趨勢,在下塔柱位置處,主塔的內力增長最為明顯,并且在塔底位置處,橋塔內力達到了最大值,這符合斜拉橋的受力特點;

2) 近斷層脈沖效應對大跨度結合梁斜拉橋地震響應有顯著的影響。以塔頂位移為例,TCU036作用下,塔頂位移為188.65 cm,而TCU070作用下,塔頂位移為105.72 cm。與非脈沖地震作用相比,近斷層脈沖效應使得塔頂位移增長了78.44%。同樣,對于塔底剪力而言,與非脈沖地震記錄TCU070作用相比,TCU036作用下,其增長了48.01%。由此說明了脈沖型地震對大跨度結合梁斜拉橋的地震響應有著不利的影響,在近斷層地區的橋梁抗震設計時,應該充分考慮近斷層脈沖效應的影響;

3) TCU036原始地震記錄以及其對應的長周期脈沖地震作用下,橋塔的位移基本完全重合,而上塔柱的內力也基本完全相等,但是下塔柱的內力存在較大差異。由此說明了對于橋塔位移而言,其地震響應主要是長周期脈沖的貢獻,但是對于橋塔內力而言,高頻成分的影響也非常大,在以內力控制為主的結構抗震設計時,應該充分考慮近斷層的長周期脈沖和高頻成分2部分對結構抗震性能的 影響。

4.2 主梁地震響應

圖10為TCU036和TCU070以及TCU036長周期脈沖作用下該算例橋梁主梁的豎向位移響應峰值包絡曲線。從圖10中可以看出:主梁的豎向位移整體呈兩邊小、中間大的趨勢,在跨中偏右的地方達到了最大值,并且左右兩側的變形顯示出明顯的不對稱。這主要是由于跨中跨度最大,剛度較小,地震作用下主梁的豎向位移最大,但是由于左右兩側結構不對稱,因此主梁豎向位移最大值出現在了跨中偏右的位置。此外,從圖10中還可以看出,與主塔地震響應類似,近斷層脈沖效應對主梁的地震響應產生了顯著的影響,以中跨豎向位移最大值為例,非脈沖地震TCU070作用下,主梁豎向位移最大值僅為68.79 cm,而脈沖地震動作用下TCU036作用下,主梁的豎向位移最大值達到了191.17 cm,近斷層脈沖效應使得主梁的豎向位移增大了177.9%。

對比TCU036及其對應長周期脈沖作用下主梁的豎向位移可以發現,TCU036原始地震記錄下主梁的豎向位移與其脈沖時程作用下主梁的豎向位移非常接近,特別是左右兩側的邊跨,二者基本完全重合。由此說明了主梁的豎向位移主要是有近斷層脈沖地震動的長周期脈沖成分控制的,高頻成分對其影響較小。

圖10 主梁豎向位移響應峰值

綜上所述可以得出,近斷層脈沖效應對大跨度結合梁斜拉橋的地震響應有著顯著的影響,近斷層脈沖效應使得大跨度斜拉橋的橋塔和主梁的變形都顯著增大,在近斷層地區的橋梁抗震設計時應該充分考慮近斷層脈沖效應的不利影響。

5 結論

1) 由于地震動測站和斷層空間相對位置的隨機性,近斷層地震記錄可能顯示出不同的脈沖特性,部分近斷層地震記錄中可能并不含有長周期、高幅值脈沖。

2) 與非脈沖地震作用相比,近斷層脈沖地震作用下,大跨度結合梁斜拉橋的地震響應顯著增大,主梁和橋塔的地震響應改變率最大分別可達到177.9%和78.44%,近斷層和跨斷層的大跨度結合梁斜拉橋抗震設計時應當充分考慮近斷層脈沖效應的不利影響。

3) 近斷層脈沖作用下大跨度結合梁斜拉橋的主塔位移響應主要是長周期成分的貢獻,高頻成分的影響基本可以忽略不計,但是對其內力響應而言,高頻成分將激起橋塔的高階振型,從而對其內力產生顯著的影響,橋塔的抗震設計應當同時考慮長周期脈沖和高頻成分2部分的影響。

[1] 劉克蛟. 三跨結合梁斜拉橋的靜力性能研究[D]. 蘭州:蘭州交通大學, 2015. LIU Kejiao. Research on static performance of the three-span composite beam cable-stayed bridge[D]. Lanzhou: Lanzhou Jiaotong University, 2015.

[2] LI S, ZHANG F, WANG J, et al. Effects of near-fault motions and artificial pulse-type ground motions on super-span cable-stayed bridge systems[J]. Journal of Bridge Engineering, 2016, 22(3): 04016128.

[3] ZHANG S, WANG G. Effects of near-fault and far-fault ground motions on nonlinear dynamic response and seismic damage of concrete gravity dams[J]. Soil Dynamics and Earthquake Engineering, 2013, 53: 217? 229.

[4] 江輝, 楚芹, 崔禹婷. 地震動脈沖參數對近斷層區深水橋墩動力響應的影響分析[J]. 鐵道學報, 2015, 37(8): 80?90. JIANG Hui, CHU Qin, CUI Yuting. Influence of pulse parameters on dynamic response of bridge pier in deep water excited by near-fault earthquake[J]. Journal of the China Railway Society, 2015, 37(8): 80?90.

[5] 陳令坤, 張楠, 胡超, 等. 近斷層地震方向脈沖效應對高速鐵路橋梁彈塑性反應的影響[J]. 振動與沖擊, 2013, 32(15): 149?155. CHEN Lingkun, ZHANG Nan, HU Chao, et al. Effects of near-fault directivity pulse-like ground motion on elastic-elastic seismic response of high-speed railway bridge[J]. Journal of Vibration and Shock, 2013, 32(15): 149?155.

[6] 夏修身, 李建中. 近場地震動對樁基礎高墩搖擺反應的影響[J]. 哈爾濱工業大學學報, 2014, 46(4): 82?86. XIA Xiushen, LI Jiangzhong. Effect of near-field ground motion on the rocking response of tall pier with pile foundations[J]. Journal of Harbin Institute of Technology, 2014, 46(4): 82?86.

[7] 孫穎, 陳天海, 卓衛東, 等. 長周期地震動作用下隔震連續梁橋地震反應特性研究[J]. 工程力學, 2016, 33(增1): 244?250. SUN Ying, CHEN Tianhai, ZHUO Weidong, et al. Study on seismic response characteristics of isolated continuous bridges under long-period ground motion[J]. Engineering Mechanics, 2016, 33(Suppl 1): 244?250.

[8] Billah A H M M, Alam M S, Bhuiyan M A R. Fragility analysis of retrofitted multicolumn bridge bent subjected to near-fault and far-field ground motion[J]. Journal of Bridge Engineering, 2012, 18(10): 992?1004.

[9] Stewart J P, Chiou S J, Bray J D, et al. Ground motion evaluation procedures for performance-based design[J]. Soil Dynamics and Earthquake Engineering, 2002, 22(9/12): 765?772.

[10] Iervolino I, Cornell C A. Probability of occurrence of velocity pulses in near-source ground motions[J]. Bulletin of the Seismological Society of America, 2008, 98(5): 2262?2277.

[11] Baker J W. Quantitative classification of near-fault ground motions using wavelet analysis[J]. Bulletin of the Seismological Society of America, 2007, 97(97): 1486? 1501.

Influence of near-fault pulse effect on seismic response of long-span composite girder cable-stayed bridge

HOU Lie1, ZHANG Longqi1, SHI Xinhu2

(1. Sichuan Vocational and Technical College of Communications, Chengdu 611130, China; 2. School of Civil Engineering, Southwest Jiaotong University, Chengdu 610031, China)

In order to study the effect of near fault pulse effect on seismic response of long-span combined beam cable-stayed bridge cable-stayed bridge, determine the contribution of long-period impulse and high-frequency component of near-fault impulse ground motion to the seismic response of long-span composite beam cable-stayed bridge. Firstly, two types of near-fault seismic records, namely, pulse/non-pulse, were selected from the Chi-Chi earthquake by using the wavelet analysis-based method, and the long period velocity pulse extracted. The amplification effect of long-period and high-amplitude velocity pulses on seismic response spectrum was studied. Then the influence of the near-fault pulse effect on the seismic response of a cable-stayed bridge with a large span combined beam was compared. The results from this research show that: Because of the flexible structure system of the long-span combined beam cable-stayed bridge, the self-vibration period is longer, and the near fault pulse effect has a significant influence on its seismic response. Compared with the non-pulse earthquake,thenear-fault impulse effect makes the seismic response of the main beam and the bridge tower increase significantly, and the maximum change rate is up to 177.9%; The seismic displacement response of girder and tower is mainly determined by the long period pulse, the influence of high frequency components is relatively small. As for the internal force response of the bridge tower, it is greatly influenced by the high frequency component of the near fault ground motion. In the seismic design based on internal force control, the effect of two parts of the long period pulse and the high frequency component should be fully considered.

near-fault impulsive; composite cable-stayed bridges; seismic response; wavelet analysis

U442.5

A

1672 ? 7029(2019)10? 2514 ? 07

10.19713/j.cnki.43?1423/u.2019.10.018

2019?01?02

國家自然科學基金資助項目(U1434205);教育部春暉計劃資助項目(Z2014032);中國鐵路總公司科技開發項目(2015G002-A)

侯烈(1978?),男,四川達州人,副教授,從事橋梁設計與研究工作;E?mail:10619975@qq.com

(編輯 陽麗霞)

猜你喜歡
橋梁
一種橋梁伸縮縫防滲水裝置
工程與建設(2019年4期)2019-10-10 01:45:56
手拉手 共搭愛的橋梁
句子也需要橋梁
加固技術創新,為橋梁健康保駕護航
中國公路(2017年11期)2017-07-31 17:56:30
無人機在橋梁檢測中的應用
中國公路(2017年10期)2017-07-21 14:02:37
高性能砼在橋梁中的應用
現代鋼橋制造對橋梁鋼的更高要求
焊接(2016年8期)2016-02-27 13:05:15
城鄉建設一體化要注重橋梁的建筑設計
南昌54座橋梁進行兩個月的夏季體檢
橋梁伸縮縫損壞因素與加固
主站蜘蛛池模板: 亚洲国产成人自拍| 国产乱人视频免费观看| 精品国产一二三区| 最新精品久久精品| 亚洲欧美在线看片AI| 中文字幕在线欧美| 99re在线视频观看| jizz在线观看| 国产香蕉在线视频| 国产精品视频第一专区| 欧美在线观看不卡| 广东一级毛片| 亚洲九九视频| 日本亚洲欧美在线| 精品国产Av电影无码久久久| 丁香综合在线| 久操线在视频在线观看| 久久国产精品夜色| 国产女人在线视频| 91久久性奴调教国产免费| 亚洲va在线观看| 一级一毛片a级毛片| 亚洲综合片| 成人精品午夜福利在线播放| 久久精品国产精品青草app| 免费精品一区二区h| 亚洲欧美成人在线视频| 丁香五月婷婷激情基地| 中文字幕 日韩 欧美| 啪啪啪亚洲无码| 日韩无码视频网站| 一本色道久久88亚洲综合| 国产jizzjizz视频| 在线a视频免费观看| AV老司机AV天堂| 亚洲欧洲日产国产无码AV| 国产女同自拍视频| 久久香蕉国产线看精品| 77777亚洲午夜久久多人| 午夜激情福利视频| 无码视频国产精品一区二区| 国产人成乱码视频免费观看| 伊人成人在线| 狼友视频一区二区三区| a天堂视频在线| 日本高清在线看免费观看| 无码乱人伦一区二区亚洲一| 国内嫩模私拍精品视频| 国产亚洲精品自在久久不卡 | 国产男人的天堂| 秋霞国产在线| 青青草国产免费国产| 91黄视频在线观看| 伊人AV天堂| 国产白浆在线观看| 亚洲Av综合日韩精品久久久| 国产毛片不卡| 亚洲精品少妇熟女| 久久永久视频| 精品少妇三级亚洲| 欧美成人精品一级在线观看| 91人妻日韩人妻无码专区精品| 国产精选自拍| 青青青草国产| 波多野结衣中文字幕一区二区| 四虎国产在线观看| 青青草原国产| 99这里只有精品免费视频| www亚洲天堂| 丁香婷婷久久| 丝袜无码一区二区三区| 精品国产免费观看一区| 国产永久免费视频m3u8| 永久免费无码日韩视频| 亚洲欧美精品在线| 99精品伊人久久久大香线蕉| 中文字幕啪啪| 秋霞午夜国产精品成人片| 国产成人精品18| 日韩乱码免费一区二区三区| 久久精品最新免费国产成人| 国产成人精品男人的天堂下载 |