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

超大跨度斜拉橋合理抗震結構體系研究

2017-04-21 06:19:35張新軍許江江浙江工業大學建筑工程學院浙江杭州310014
浙江工業大學學報 2017年2期
關鍵詞:結構

張新軍,許江江(浙江工業大學 建筑工程學院,浙江 杭州 310014)

超大跨度斜拉橋合理抗震結構體系研究

張新軍,許江江
(浙江工業大學 建筑工程學院,浙江 杭州 310014)

為研究超大跨度斜拉橋的抗震性能,以主跨1 400 m的超大跨度斜拉橋設計方案為工程背景,采用多振型地震反應譜方法進行水平和豎向地震作用的結構反應分析,揭示超大跨度地震反應的特點,同時分析了主梁的高度和寬度、索塔結構型式、索塔高跨比、邊中跨比、邊跨輔助墩以及斜拉索錨固體系等設計參數對超大跨度斜拉橋地震反應的影響,探討了其合理的抗震結構體系.結果表明:水平地震作用下橋塔和主梁的地震反應顯著,橋塔的塔底截面和主梁在塔梁交接處截面的地震內力非常大,應特別重視這些截面的抗震設計;超大跨度斜拉橋當采用A型橋塔并降低其高度、短邊跨并設置輔助墩以及斜拉索采用部分地錨方式時,結構地震反應小,抗震性能良好,是其合理的抗震結構體系.

超大跨度斜拉橋;抗震性能;多振型地震反應譜分析;結構設計參數

本世紀初,蘇通長江大橋(主跨1 088 m)和香港昂船洲大橋(主跨1 018 m)的相繼建成使得斜拉橋跨入了千米級橋梁的行列,更大跨徑的俄羅斯Russky Island大橋(主跨1 104 m)也于2012年建成.為避免跨海聯島工程建設中超深水基礎施工困難并滿足日益增長的通航要求,斜拉橋的跨徑將會進一步增大,如規劃的韓國東南部連接馬山市和Geoje島的1 200 m主跨的斜拉橋方案以及日本本四聯絡線的1 400 m主跨的斜拉橋方案等.較之于懸索橋,斜拉橋結構剛度更大,抗風安全性更好,拉索容易更換,施工簡便以及不需要大體積錨碇,因而在近年來的國際跨海大橋設計方案競賽中屢次被提出和采用,斜拉橋已經成為現代大跨度橋梁的主流橋型[1].已有關于斜拉橋的極限跨徑研究也表明:主跨1 200 m以下是斜拉橋比較經濟合理的主跨,主跨在1 200~1 500 m范圍內的斜拉橋仍具有競爭力[1-2].

隨著斜拉橋跨度的進一步增大,結構剛度將進一步降低,結構的抗震安全性已成為其設計考慮的重要問題.迄今為止,國內學者對千米級主跨斜拉橋的抗震性能及減震方法開展了分析研究.武芳文等基于隨機振動理論,對蘇通長江公路大橋在隨機地震荷載作用下的動力響應展開研究,分析了結構彈性模量、樁-土-結構動力相互作用、地震動模型、斜拉索振動效應以及模態組合數量對結構內力和位移的影響[3].王蒂等通過結構非線性地震分析研究了蘇通長江公路大橋塔梁連接處分別設置彈性連接裝置和阻尼器的減震效果[4].葉愛君等以蘇通長江公路大橋為研究對象,采用非線性時程分析方法分析了3種邊墩和梁橫向約束體系對超大跨度斜拉橋地震反應的影響以及阻尼器的合理設置方式及其設計參數[5].任亮以主跨1 090 m的鋼索鋼主梁斜拉橋設計方案為原型擬定了一座相同跨度的CFRP索RPC梁斜拉橋設計方案,分別采用反應譜法和時程分析法對比分析了原方案和新方案的彈性地震響應,并對所提方案進行了抗震性能評價[6].韓振峰等考慮幾何非線性對兩座千米級斜拉橋進行成橋狀態的地震動力時程分析,分析幾何非線性對千米級斜拉橋的地震反應特性的影響[7].韓振峰等針對漂浮型的千米級斜拉橋,分析研究在塔梁之間分別采用彈性連接裝置、流體黏滯阻尼器以及彈性連接裝置與流體黏滯阻尼器組合的3種連接方式的縱向減震效果及其參數設置[8].以上研究主要針對主跨略超千米斜拉橋的抗震性能尤其是減震措施,對于更大跨度斜拉橋的抗震性能及其合理抗震結構體系方面的研究則沒有涉及.為此,筆者以主跨1 400 m的斜拉橋設計方案為工程背景,采用多振型地震反應譜分析方法進行水平和豎向地震作用下的結構反應分析,闡述超大跨度斜拉橋地震反應的特點.在此基礎上,分析主梁高度和寬度、索塔結構型式、索塔高跨比、邊中跨比、輔助墩設置及斜拉索錨固體系等結構設計參數對超大跨度斜拉橋地震反應的影響,并探討其合理的抗震結構體系.

1 橋梁簡介

圖1為一座主跨1 400 m的超大跨度斜拉橋設計方案[9],縱橋向結構對稱布置,兩側邊跨長度均為680 m,為提高結構整體剛度在兩側邊跨靠近錨固墩位置設置間距為100 m的3個輔助墩.主梁采用流線型鋼箱梁,梁寬35 m,梁高3.5 m,并在橋塔兩側各80 m(Xu)范圍內對主梁截面進行加強.橋塔采用全鋼結構,橫橋向A型結構型式,索塔總高度約327 m,其中橋面以上高度約287 m,索塔高跨比為1/5.斜拉索采用全自錨體系,在主梁上的錨固間距為20 m,在索塔頂部的錨固間距為4 m,全橋共設有4×34對斜拉索.主梁、索塔和斜拉索的截面特性見表1.

圖1 主跨1 400 m斜拉橋設計方案立面布置圖(單位:m)Fig.1 Vertical layout of a cable-stayed bridge design scheme with 1 400 m main span(unit: m)

構件A/m2Ix/m4Iy/m4Iz/m4W/(kN·m-1)主梁1.761(2.046)8.330(9.739)3.939(4.434)193.200(261.1)258.400(280.4)索塔1.76039.27030.67040.320189.600斜拉索0.038~0.0092.926~0.693

注:1)A為截面積;Ix為抗扭慣性矩;Iy為豎向抗彎慣性矩;Iz為側向抗彎慣性矩;W為構件單位長度重量;括號內數值為加強截面的特性值.

2 結構三維有限元分析模型

地震反應分析時,基于MIDAS/CIVIL有限元分析軟件,將該橋離散為三維桿系結構有限元模型,如圖2所示.主梁采用魚骨梁計算模型,主梁、索塔及其橫梁等均采用空間梁單元模擬,并考慮其P—⊿效應;斜拉索采用空間桿單元模擬,并考慮其初始軸向力對單元剛度的影響;斜拉索與主梁之間連接關系采用剛臂單元模擬.主梁縱橋向可自由位移,與錨固墩及輔助墩之間保持沿橫橋向和豎向位移以及繞縱軸轉動三個自由度的主從關系,其余自由度均可自由位移;在塔梁交接處,主梁豎向無支撐,但其橫橋向位移則受橋塔約束.

圖2 結構三維有限元分析模型Fig.2 Structural 3D finite element model

3 超大跨度斜拉橋地震反應特性分析

3.1 地震動參數

依據《公路橋梁抗震設計細則(JTG/T B02-01—2008)》[10]并結合橋梁類型和橋址處場地土特性,該橋為A類橋梁,按地震烈度7度設防;水平向設計基本地震動加速度峰值為0.15 g;場地類型為Ⅱ類,場地土特征周期為0.40 s;結構阻尼比為3%.圖3為E1地震作用下的水平設計加速度反應譜,其豎向設計加速度反應譜取為水平設計加速度反應譜的65%.

圖3 E1地震作用下水平設計加速度反應譜Fig.3 The horizontal seismic design acceleration response spectrum under E1 seismic action

3.2 地震反應分析

基于MIDAS/CIVIL有限元分析軟件,采用多振型地震反應譜分析方法,對該方案橋進行縱向、橫向和豎向地震作用的結構反應分析,各振型的地震反應采用CQC方法組合.限于篇幅,僅給出結構地震反應的峰值,如表2所示.

表2 E1地震動作用下結構地震反應1)Table 2 Structural seismic response under E1 seismic action

注:1) 縱向和豎向地震作用下,索塔彎矩和剪力均為縱向,主梁彎矩和剪力均為豎向;橫向地震作用下,索塔和主梁的彎矩和剪力均為橫向.

結合結構的地震反應分布圖和表1的結構地震反應峰值進行綜合分析,超大跨度斜拉橋的地震反應特點總結如下:

1) 縱向地震作用下,橋塔縱向振動,主梁表現為縱向和豎向耦合振動.橋塔塔頂處的縱向位移最大,并在塔底截面產生最大縱橋向彎矩、剪力和軸力;主梁沿橋軸向的縱向位移基本一致,并在中跨兩四分點處和邊跨靠近橋塔的1/3左右處產生最大豎向位移.主梁的最大彎矩和剪力均出現在邊跨靠近橋塔的第一個輔助墩處.相比較而言,橋塔的地震反應遠大于主梁,縱向地震作用對橋塔受力不利,并應特別重視塔底截面的抗震設計.

2) 橫向地震作用下,橋塔和主梁同時橫橋向振動.橋塔橫橋向彎曲,塔頂處發生最大橫向位移,并在塔底產生最大橫向彎矩、剪力和軸向力;主梁的最大橫向位移出現在跨中,橫向彎矩和剪力的最大值則出現在塔梁交接處截面.相比較而言,主梁的橫向地震反應更為顯著,橫向地震作用對主梁受力不利,同時應特別重視塔梁交接處主梁截面的抗震設計.

3) 豎向地震作用下,橋塔和主梁的振動形式與縱向地震作用相同.豎向地震作用下的主梁的慣性荷載通過斜拉索傳遞到橋塔中去,橋塔所受的軸向力因而顯著增大,伴隨著橋塔的縱向彎曲,在塔底截面產生了較大的縱向彎矩和剪力;主梁在跨中處的豎向位移及彎矩均達到最大值.

4) 縱向和橫向地震作用下結構的地震反應均顯著大于豎向地震作用效應,因此應重視水平地震作用下超大跨度斜拉橋的抗震性能.同時,在縱向、橫向和豎向地震綜合作用下,橋塔的塔底截面和塔梁交接處主梁截面的地震反應均非常大,應特別重視這些截面的抗震設計.

4 超大跨度斜拉橋地震反應的設計參數影響分析

為全面地把握超大跨度斜拉橋的抗震性能,對斜拉橋主要設計參數如主梁的高度和寬度、邊中跨比、索塔高跨比、橋塔結構形式、邊跨輔助墩設置以及斜拉索的錨固體系等對超大跨度斜拉橋地震反應的影響進行了分析,并探討其合理的抗震結構體系.限于篇幅,以下的各項內力和位移值按縱向、橫向和豎向地震單獨作用產生的峰值采用平方和開根的計算方法得到.

4.1 主梁高度

主梁高度是影響斜拉橋主梁豎向彎曲剛度的重要參數,增加主梁高度可以增大結構的豎向剛度,減小結構變形改善其靜力性能[11-12].基于設計方案橋,分別設計了主梁高度為3,4 m的兩座對比方案橋并進行地震反應分析,主梁高度對索塔和主梁地震反應的影響分別如表3,4所示.

可以看出:主梁高度對斜拉橋索塔和主梁的地震反應影響比較小,總體而言隨著主梁高度的增加,結構的地震位移略有減小,但地震內力則有所增大.因此,主梁高度是影響斜拉橋地震反應的非敏感參數,可以根據結構靜力性能需求確定.

4.2 主梁寬度

基于設計方案橋,將主梁寬度分別調整為28,32 m設計了相應的兩座對比方案橋,并進行結構地震反應分析,主梁寬度對索塔和主梁地震反應的影響分別如表5,6所示.

表3 主梁高度對索塔地震反應的影響Table 3 Effect of the girder depth on the seismic response of tower

表4 主梁高度對主梁地震反應的影響Table 4 Effect of the girder depth on the seismic response of girder

表5 主梁寬度對索塔地震反應的影響Table 5 Effect of the girder width on the seismic response of tower

表6 主梁寬度對主梁地震反應的影響Table 6 Effect of the girder width on the seismic response of girder

同主梁高度參數一樣,主梁寬度的變化對斜拉橋索塔和主梁的地震反應影響也很小.總體上看,減小主梁寬度有利于減小結構地震位移和內力,改善其抗震性能.因此,在滿足設計交通量需求的前提下,采用較小的主梁寬度有利于改善斜拉橋結構的抗震性能.

4.3 索塔結構型式

此處的索塔結構型式主要指其橫橋向結構布置型式,它對斜拉橋結構的橫向和扭轉剛度影響比較顯著[11-12].為揭示橫橋向索塔結構型式對超大跨度斜拉橋抗震性能的影響,基于設計方案橋,將橋塔斜拉索錨固區段的兩根塔柱并攏設計了橋塔橫橋向倒Y型的方案橋,并進行地震反應分析,索塔結構型式對索塔和主梁地震反應的影響分別如表7,8所示.

表7 索塔結構型式對索塔地震反應的影響Table 7 Effect of the tower structural style on the seismic response of tower

表8 索塔結構型式對主梁地震反應的影響Table 8 Effect of the tower structural style on the seismic response of girder

可以看出:與A型索塔結構型式相比,斜拉橋采用倒Y型橋塔后,索塔的主梁的地震內力雖略有減小,但橋塔的縱向位移、主梁的橫向和豎向位移都顯著增大,說明結構的整體剛度有所降低.因此,從抗震性能角度而言,超大跨度斜拉橋適宜采用A型塔.

4.4 索塔高跨比

索塔高跨比是影響斜拉橋的整體剛度和經濟性的重要設計參數,比值一般居于1/4與1/7之間,而以1/5附近取值較常見[11-12].為揭示索塔高跨比對超大跨度斜拉橋抗震性能的影響,通過調整橋塔高度建立索塔高跨比分別為1/6和1/4的兩座對比方案橋,并進行地震反應分析,索塔高跨比對索塔和主梁地震反應的影響分別如表9,10所示.

表9 索塔高跨比對索塔地震反應的影響Table 9 Effect of the tower height-to-span ratio on the seismic response of tower

表10 索塔高跨比對主梁地震反應的影響Table 10 Effect of the tower height-to-span ratio on the seismic response of girder

隨著索塔高跨比的減小,斜拉橋索塔和主梁的各向位移都顯著減小,其所受的地震內力隨之降低,結構的抗震性能因而得到明顯改善.此外,降低索塔高度也有利于提高經濟性.因此,在滿足靜力性能要求的前提下,盡可能采用較小的索塔高跨比.

4.5 邊中跨比

邊中跨比對斜拉橋的力學和經濟性能有重要影響,取值一般在0.25到0.5之間[11-12].為揭示邊中跨比對超大跨度斜拉橋抗震性能的影響,保持中跨長度不變,改變邊跨斜拉索的錨固間距分別擬定了邊跨為408,544 m(邊中跨比分別為0.29,0.39)兩座對比方案橋,并進行地震反應分析,邊中跨比對索塔和主梁地震反應的影響分別如表11,12所示.

表11 邊中跨比對索塔地震反應的影響Table 11 Effect of the side-to-main span ratio on the seismic response of tower

表12 邊中跨比對主梁地震反應的影響Table 12 Effect of the side-to-main span ratio on the seismic response of girder

可以看出:邊跨長度對橫橋向地震內力和位移峰值影響很小,但對縱向和豎向地震反應影響比較顯著.縮短邊跨長度后,橋塔和主梁的縱向和豎向地震反應隨之明顯減小,結構的抗震性能因而提高.因此,再次從抗震性能角度說明了斜拉橋適宜采用短邊跨的橋跨布置方式.

4.6 輔助墩設置

為增強結構的豎向剛度并提高施工安全性,大跨度斜拉橋設計時通常會在邊跨設置若干數量的輔助墩[11-12].為揭示邊跨輔助墩設置對超大跨度斜拉橋抗震性能的影響,在兩側邊跨分別設置1個和2個輔助墩擬定兩座對比方案橋,并進行地震反應分析,邊跨輔助墩設置數量對索塔和主梁地震反應的影響分別如表13,14所示.

表13 輔助墩數量對索塔地震反應的影響Table 13 Effect of the auxiliary piers on the seismic response of tower

表14 輔助墩數量對主梁地震反應的影響Table 14 Effect of the auxiliary piers on the seismic response of girder

可以看出:隨著邊跨輔助墩設置數量的增加,斜拉橋索塔和主梁的地震反應都隨之單調遞減,說明邊跨輔助墩的設置有利于提高超大跨度斜拉橋的抗震性能.應當指出的是,邊跨輔助墩的最優數量需要結合經濟性、結構剛度以及施工安全性綜合確定.

4.7 斜拉索錨固體系

目前,大跨度斜拉橋基本都屬于自錨式斜拉橋,即斜拉索全部在橋面主梁上錨固.為揭示斜拉索錨固體系對超大跨度斜拉橋抗震性能的影響,將設計方案橋邊跨兩側索面的最外側5根拉索由自錨方式改為地錨方式設計了部分地錨式斜拉橋對比方案[11-12],并進行地震反應分析,斜拉索錨固體系對索塔和主梁地震反應的影響分別如表15,16所示.

表15 斜拉索錨固體系對索塔地震反應的影響Table 15 Effect of the anchorage system of stay cables on the seismic response of tower

表16 斜拉索錨固體系對主梁地震反應的影響Table 16 Effect of the anchorage system of stay cables on the seismic response of girder

可以看出:與斜拉索全自錨方式相比,斜拉橋邊跨部分斜拉索改為地錨后,橋塔的地震內力雖有所增大,但其地震位移則明顯減小,說明斜拉索對橋塔的約束剛度得到增強;主梁的橫向地震內力略有增加,豎向地震內力則有所減小,主梁的各向位移都顯著減小,說明斜拉索對主梁的支承剛度得到增強.因此,相比于全自錨式斜拉橋,部分地錨式斜拉橋的抗震性能更加優越.

5 結 論

基于MIDAS/CIVIL有限元分析軟件,采用多振型地震反應譜方法,對1 400 m主跨超大跨度斜拉橋進行E1地震作用下的結構地震反應特性和主要結構設計參數的影響分析,并得到了兩點結論:1) 超大跨度斜拉橋在縱向和橫向地震作用下的地震反應都顯著大于豎向地震作用,因此應重視水平地震作用下結構的抗震性能.同時,在縱向、橫向和豎向地震綜合作用下,橋塔的塔底和塔梁結合處主梁截面的地震內力都非常大,應特別重視這些截面的抗震設計.2) 超大跨度斜拉橋當采用A型橋塔并降低其高度、短邊跨布置并設置若干輔助墩以及邊跨部分拉索采用地錨方式時,結構地震反應小,抗震性能良好,是其合理的抗震結構體系.

[1] 項海帆.世界大橋的未來趨勢——2011年倫敦國際橋協會議的啟示[J].橋梁,2012(3):12-16.

[2] GIMSING N J, GEORGAKIS C T. Cable-supported bridges: concept and design[M]. Chichester: Wiley,2012.

[3] 武芳文,徐超,趙雷.超大跨度斜拉橋隨機地震響應參數敏感性分析[J].鐵道學報,2014,36(6):107-113.

[4] 王蒂,黃平明.超大跨度斜拉橋縱向減震耗能塔、梁連接裝置研究[J].鄭州大學學報(工學版),2008,29(4):112-115.

[5] 葉愛君,范立礎.超大跨度斜拉橋的橫向約束體系[J].中國公路學報,2007,20(2):63-67.

[6] 任亮.基于高性能材料的千米級跨徑混凝土斜拉橋力學性能研究[D].長沙:湖南大學,2013.

[7] 韓振峰,葉愛君,范立礎.千米級斜拉橋的動力幾何非線性分析[J].土木工程學報,2010,43(6):67-73.

[8] 韓振峰,葉愛君.千米級斜拉橋的縱向減震體系研究[J].地震工程與工程振動,2015,35(6):64-70.

[9] NAGAI M, FUJINO Y, YAMAGUCHI H, IWASAKI E. Feasibility of a 1400 m span steel cable-stayed bridge[J]. Journal of bridge engineering, ASCE,2004,9(5):444-452.

[10] 中華人民共和國行業標準.公路橋梁抗震設計細則:JTG/T B02-01—2008[S].北京:人民交通出版社,2008.

[11] 張新軍,孫海凌.超大跨度斜拉橋空氣靜力和動力穩定性研究[J].中國工程科學,2014,16(3):50-58.

[12] 張新軍,虞周均,孫海凌.超大跨度斜拉橋空氣靜力穩定性研究[J].浙江工業大學學報,2014,42(2):182-189.

(責任編輯:陳石平)

Study on the favorable earthquake-resistant structural system of super long-span cable-stayed bridges

ZHANG Xinjun, XU Jiangjiang
(College of Civil Engineering and Architecture, Zhejiang University of Technology, Hangzhou 310014, China)

By taking the design scheme of a super long-span cable-stayed bridge with a main span of 1400 m as an example, the seismic response of the bridge under horizontal and vertical seismic excitations is investigated through multimode response spectrum analysis and the seismic performance of the bridge is revealed. The influence of design parameters, including the depth and width of the girder, the structural style and height-to-span ratio of the tower, the side-to-main span ratio, the auxiliary piers in the side span and the anchorage system of stay cables, on the seismic response of the bridge is analyzed numerically and its favorable earthquake-resistant structural system is proposed. The results show that the horizontal seismic excitation produces significant seismic response of the girder and tower, there exists great seismic response for the sections at the tower bottom and the girder near the tower, and therefore more attentions should be paid to the seismic design of these sections. For A-shaped towers with smaller tower height, shorter side span, several auxiliary piers installed in the side span and several stay cables partially earth-anchored, a smaller seismic response and a better seismic performance are achieved for super long-span cable-stayed bridges.

super long-span cable-stayed bridge; seismic performance; multimode response spectrum analysis; design parameter

2016-06-15

浙江省自然科學基金資助項目(LY15E080020)

張新軍(1971—),男,浙江東陽人,教授,主要從事橋梁抗風抗震研究,E-mail:xjzhang@zjut.edu.cn.

U448.27

A

1006-4303(2017)02-0230-07

猜你喜歡
結構
DNA結構的發現
《形而上學》△卷的結構和位置
哲學評論(2021年2期)2021-08-22 01:53:34
論結構
中華詩詞(2019年7期)2019-11-25 01:43:04
新型平衡塊結構的應用
模具制造(2019年3期)2019-06-06 02:10:54
循環結構謹防“死循環”
論《日出》的結構
縱向結構
縱向結構
我國社會結構的重建
人間(2015年21期)2015-03-11 15:23:21
創新治理結構促進中小企業持續成長
現代企業(2015年9期)2015-02-28 18:56:50
主站蜘蛛池模板: 人妻丰满熟妇AV无码区| 日韩第一页在线| 亚洲无码视频一区二区三区| 无码AV高清毛片中国一级毛片| 在线亚洲精品自拍| 亚洲国产精品日韩av专区| 天天综合色天天综合网| 免费Aⅴ片在线观看蜜芽Tⅴ| 99精品在线看| 国产微拍精品| 97久久人人超碰国产精品| 久久久久亚洲av成人网人人软件 | 国产精品妖精视频| www亚洲天堂| 97在线碰| 国产精品林美惠子在线观看| 日日噜噜夜夜狠狠视频| 欧美在线网| 日本伊人色综合网| 伊人久综合| 欧美爱爱网| 91无码人妻精品一区二区蜜桃| 特级欧美视频aaaaaa| 国产香蕉一区二区在线网站| 久久黄色小视频| 一本无码在线观看| 久久一本日韩精品中文字幕屁孩| 亚洲国产亚综合在线区| 99精品一区二区免费视频| 激情网址在线观看| 亚洲色欲色欲www网| 久久久久88色偷偷| 找国产毛片看| 亚洲欧洲日韩国产综合在线二区| 99久久国产综合精品2023| 国产小视频免费观看| 国产SUV精品一区二区| a级高清毛片| 一级爆乳无码av| 欧美啪啪一区| 极品国产一区二区三区| 任我操在线视频| 六月婷婷激情综合| 天天躁日日躁狠狠躁中文字幕| 久久伊伊香蕉综合精品| 久久国产免费观看| 2021最新国产精品网站| 亚洲成人精品久久| 欧美亚洲香蕉| 狂欢视频在线观看不卡| 一本大道东京热无码av| 黄色污网站在线观看| 久久国产精品国产自线拍| 中文字幕佐山爱一区二区免费| 91久久精品日日躁夜夜躁欧美| 久久五月视频| 国产精彩视频在线观看| 99爱视频精品免视看| 99热这里只有精品免费| 999精品在线视频| 国产欧美成人不卡视频| 色婷婷狠狠干| 草草影院国产第一页| 2021亚洲精品不卡a| 国产精品美女自慰喷水| 国产高清不卡视频| 亚洲精品天堂自在久久77| 日韩精品无码免费一区二区三区 | 国产h视频在线观看视频| 久久综合丝袜日本网| 欧美亚洲中文精品三区| 啊嗯不日本网站| 亚洲国产第一区二区香蕉| 亚洲日韩精品伊甸| 一级毛片高清| 97av视频在线观看| 直接黄91麻豆网站| 国产乱子伦一区二区=| 国产91无码福利在线| 日韩成人在线视频| 亚洲av日韩综合一区尤物| 在线观看av永久|