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

新意法在野豬山隧道施工中的應(yīng)用研究

2017-06-28 16:22:44王克忠朱海斌孔令民梁其東張海軍
關(guān)鍵詞:圍巖變形施工

王克忠,朱海斌,孔令民,梁其東,張海軍

(1.浙江工業(yè)大學(xué) 建筑工程學(xué)院,浙江 杭州 310014;2.日照市供水管理處,山東 日照 276800)

新意法在野豬山隧道施工中的應(yīng)用研究

王克忠1,朱海斌1,孔令民2,梁其東2,張海軍1

(1.浙江工業(yè)大學(xué) 建筑工程學(xué)院,浙江 杭州 310014;2.日照市供水管理處,山東 日照 276800)

為探討新意法在復(fù)雜地質(zhì)環(huán)境下隧道施工中的可行性,依托具體工程,對山嶺隧洞出口段進行了數(shù)值仿真分析,研究了新意法和臺階法開挖時掌子面的縱向內(nèi)空變形和應(yīng)力變化規(guī)律,預(yù)測了沿隧道軸線方向拱頂沉降的時空變化規(guī)律,對比分析了隧道不同開挖、臨時加固方式對掌子面大變形、失穩(wěn)的影響,得出了掌子面中部易失穩(wěn),新意法施工可有效抑制掌子面變形的結(jié)論,對同類隧道工程研究具有一定的理論意義,并對工程實踐具有重要參考價值.

新意法;臺階法;掌子面;數(shù)值模擬

山嶺隧道出口段圍巖破碎,在隧道開挖后,掌子面應(yīng)力解除,圍巖應(yīng)力重分布[1],隧道出口段坍塌、冒頂和掌子面失穩(wěn)等現(xiàn)象時常發(fā)生.多年來,國內(nèi)外一些專家、學(xué)者采用解析法、模型試驗法和經(jīng)驗公式法等[2]對隧道出口段開展了大量研究.周藝等[3]以現(xiàn)場試驗結(jié)果為依據(jù),結(jié)合數(shù)值仿真模擬,對比分析了隧道在不同開挖方法下的圍巖變形特征、初期支護和二襯的內(nèi)力變化規(guī)律.趙錄學(xué)[4]依據(jù)新意法原理,對特殊工程環(huán)境下隧道選線、支護參數(shù)的優(yōu)化和施工方案進行了探討,為新意法在我國的推廣應(yīng)用提供了借鑒.國內(nèi)外參考文獻表明,國內(nèi)學(xué)者利用新意法的理論,對超前預(yù)加固在軟巖淺埋隧道施工的作用進行了研究[5-7],但未能全面闡釋施工中圍巖力學(xué)與變形規(guī)律.筆者依托野豬山隧道工程,采用有限元軟件Midas/GTS對新意法和臺階法的施工過程進行了數(shù)值仿真計算,并結(jié)合監(jiān)測數(shù)據(jù)對比分析了兩種施工方法在淺埋軟弱地層中的適用性,為今后同類工程的設(shè)計和施工提供了借鑒.

1 掌子面失穩(wěn)判據(jù)

巖石強度理論是與隧道掌子面的失穩(wěn)機理緊密相關(guān)的,反映了巖石破壞時的應(yīng)力狀態(tài)或變形狀態(tài).目前廣泛應(yīng)用的強度理論,多數(shù)是從應(yīng)力觀點來考察材料破壞的,其中庫倫-莫爾強度理論就是應(yīng)力觀點強度理論的繼續(xù)和發(fā)展.

庫倫在1773年提出了,當(dāng)材料某截面上的剪應(yīng)力達到黏聚力與內(nèi)摩擦力之和時,材料將發(fā)生剪切破壞.土體剪切破壞準(zhǔn)則的庫倫準(zhǔn)則表達式[8]為

τn=c-σntanφ

(1)

式中:τn為極限抗剪強度;σn為剪切面上的法向應(yīng)力;c,φ分別為土的黏聚力和內(nèi)摩擦角.

摩爾(1882)將庫倫破壞條件表示為與最大主應(yīng)力圓相切的一條直線,當(dāng)σ1≥σ2≥σ3時,由圖1可推導(dǎo)各主應(yīng)力關(guān)系式為

(2)

屈服面方程為

(3)

圖1 摩爾-庫侖破壞條件Fig.1 Mohr-Coulomb failure conditions

如何判別掌子面是否處于整體失穩(wěn)狀態(tài)是隧道有限元穩(wěn)定分析的一個關(guān)鍵問題.按照巖土力學(xué)中的摩爾-庫倫屈服準(zhǔn)則,隧道掌子面失穩(wěn)時破壞面上各點應(yīng)力都達到極限平衡狀態(tài),極限應(yīng)力圓與屈服面相切.因此,有限元分析時以掌子面在開挖過程中應(yīng)力狀態(tài)與屈服面的關(guān)系作為失穩(wěn)判據(jù)是比較合理的.

2 數(shù)值仿真模擬

2.1 工程簡介

野豬山隧道為分離式隧道,兩隧洞相距約35m.其出口段分布殘坡積、洪積含碎石粉質(zhì)黏土和含黏性土碎石,褐黃色,碎石較多為強風(fēng)化巖塊,厚度較大,埋深為10~20m.兩隧洞出口段分別采用新意法和臺階法兩種施工方法,兩種工法初期支護一致,采用25cm厚C25噴射混凝土、錨桿、鋼筋網(wǎng)與鋼架聯(lián)合支護,鋼拱架為18號工字鋼,縱向間距為0.5m.

2.2 模型及邊界條件

采用Midas/GTS建立三維有限元模型,由于野豬山兩隧洞相距較遠,相互影響可忽略不計,故對單洞施工過程進行模擬.選取野豬山隧道出口段50m進行數(shù)值仿真分析.以模型正面隧道底部中心為原點,取隧道軸線在水平面上的投影線為y軸,其正方向由洞口指向隧洞內(nèi)部,x軸垂直隧洞軸線,z軸的正方向為豎直向上.模型范圍:-40 m≤x≤40 m,0≤y≤50 m,-35 m≤z≤地面.該模型上邊界設(shè)為自由邊界.以隧道開挖跨度的3~5倍[9]選取模型左右邊界和底部邊界,其中模型底部邊界設(shè)置豎向約束,其余各邊界分別設(shè)置該方向上的水平約束,計算模型如圖2所示.

圖2 計算模型Fig.2 Calculation model

2.3 模型參數(shù)選定

計算本構(gòu)模型采用Mohr-coulomb屈服準(zhǔn)則,選取彈塑性模型進行計算.采用板單元模擬噴射混凝土,用植入桁架單元模擬錨桿,鋼拱架采用等效計算的方法[10]并參考型鋼型號將彈性模量折算并附加到混凝土噴層中,其中折算噴混層的縱向長度為0.5 m,計算式為

E=E0+Sg·Eg/Sc

(4)

式中:E為折算后等效初期支護的彈性模量;E0為原混凝土材料的彈性模量;Sg為鋼拱架的截面積;Sc為混凝土截面積;Eg為鋼材的彈性模量.

參照野豬山隧道工程地質(zhì)勘察報告中的地層物理參數(shù),土體劃分為3層.具體物理力學(xué)參數(shù)如表1,2所示.

表1 巖體計算參數(shù)

表2 支護計算參數(shù)

2.4 施工模擬

為比較臺階法和新意法施工的優(yōu)缺點及其對掌子面穩(wěn)定性的影響,對兩種工法下隧道的開挖及支護過程進行了對比分析研究.臺階法施工:采用三臺階法開挖,先后開挖第一步和第二步臺階,各開挖8 m和4 m并支護,再以1 m為一個循環(huán)進尺,每一個循環(huán)完成后及時施做初期支護,施工工序如圖3所示.

新意法施工:采用全斷面開挖,以1 m為一個循環(huán)進尺,開挖前掌子面植入玻璃纖維錨桿并注漿進行超前預(yù)加固,保證錨桿嵌入掌子面前方至少12 m,搭接6 m,掌子面錨桿布置間距為120 mm,如圖4所示.

圖3 臺階法開挖工序Fig.3 Process of bench cut method

圖4 新意法掌子面錨桿布置Fig.4 Arrangement of anchors on excavated face of ADECO-RS approach

3 計算結(jié)果及分析

3.1 掌子面擠出位移

巖土控制變形分析法依據(jù)超前核心圍巖體系的應(yīng)力—應(yīng)變狀態(tài),把掌子面穩(wěn)定形態(tài)分為3類[11-12]:穩(wěn)定狀態(tài)、短期穩(wěn)定狀態(tài)和不穩(wěn)定狀態(tài).而掌子面超前核心巖土體穩(wěn)定性將直接影響到掌子面的擠出位移.掌子面擠出位移超過某個界限就表明掌子面穩(wěn)定狀態(tài)轉(zhuǎn)變?yōu)槎唐诜€(wěn)定或不穩(wěn)定狀態(tài).

三臺階法開挖時,上、中、下三個臺階的掌子面沿隧道縱向產(chǎn)生不同程度的擠出位移,土體出現(xiàn)鼓出現(xiàn)象,其中最大位移出現(xiàn)在下臺階,達到15 mm左右,如圖5所示(單位為m).

新意法施工時,掌子面同樣出現(xiàn)擠出位移,最大位移出現(xiàn)在掌子面中心附近,達到8 mm左右,相比于臺階法開挖縮小了近一半,如圖6所示(單位為m).新意法施工改變了開挖面處超前核心土的強度和剛度,一定程度上抑制了掌子面過大變形.

圖5 臺階法掌子面縱向位移圖Fig.5 Longitudinal displacement on excavated face of bench cut method

圖6 新意法掌子面縱向位移圖Fig.6 Longitudinal displacement on excavated face of ADECO-RS approach

3.2 掌子面應(yīng)力狀態(tài)

根據(jù)掌子面的擠出位移可知:在施工過程中掌子面穩(wěn)定性最差的部位在掌子面中心附近.下面選取臺階法上、中、下三個臺階各開挖到30 m處時,與三個臺階相對應(yīng)的掌子面中心附近的節(jié)點作為研究對象,同時新意法三個測點的位置與臺階法的三個節(jié)點相對應(yīng),然后通過偏應(yīng)力和球應(yīng)力的關(guān)系曲線,分析以上節(jié)點在開挖過程中的屈服軌跡和應(yīng)力變化趨勢.其中,圍巖采用摩爾-庫倫準(zhǔn)則[13],以其Mohr形式作為屈服函數(shù),即

(5)

式中:σ1,σ3分別為巖土體剪切面上的第一、第三主應(yīng)力;c,φ分別為巖土體的黏聚力和內(nèi)摩擦角.

兩種施工方法在開挖過程中,上部掌子面中心點的應(yīng)力變化情況如圖7所示.兩種方法在隨著開挖面接近上部中心監(jiān)控位置,球應(yīng)力都逐漸增大,同時偏應(yīng)力也逐漸增大.其中,臺階法的應(yīng)力點向屈服面接近,但并沒有達到屈服線,掌子面處于穩(wěn)定狀態(tài)和短期穩(wěn)定狀態(tài).而新意法由于施工時對掌子面進行植入錨桿并注漿處理,改善了開挖面前方的巖土體物理力學(xué)性質(zhì),應(yīng)力點則超過原有的屈服面,但應(yīng)力變化并處于彈塑性階段,即掌子面處于穩(wěn)定狀態(tài)和短期穩(wěn)定狀態(tài).

圖7 上部掌子面應(yīng)力變化Fig.7 Face stress of upper face

中部掌子面隨著開挖,其中心區(qū)域應(yīng)力變化情況如圖8所示.中部掌子面的監(jiān)控點應(yīng)力隨著開挖面的推進,兩種開挖方法的球應(yīng)力都逐漸減小,而剪應(yīng)力變化出現(xiàn)差異.臺階法在開挖面距監(jiān)控點4m時,剪應(yīng)力幾乎不隨球應(yīng)力的減小發(fā)生變化,當(dāng)開挖面距監(jiān)控點2m時,圍巖應(yīng)力達到M-C屈服線而進入塑性狀態(tài),即掌子面失穩(wěn).而新意法開挖過程中,其應(yīng)力點超過了M-C屈服線,但圍巖應(yīng)力并沒有發(fā)生明顯的屈服現(xiàn)象.計算表明,由于隧道所處地層軟弱,開挖面前方2m范圍內(nèi)的土體難以維持平衡狀態(tài),而通過對開挖面的超前核心土進行加固,可以增強開挖面前方2m內(nèi)的巖土體的穩(wěn)定性.

圖8 中部掌子面應(yīng)力變化Fig.8 Face stress of middle face

下部掌子面隨著開挖,其中心區(qū)域應(yīng)力變化情況如圖9所示.臺階法的監(jiān)控點應(yīng)力變化情況和中部掌子面類似,在開挖面距監(jiān)控點4m時,剪應(yīng)力幾乎不隨球應(yīng)力的減小發(fā)生變化,當(dāng)開挖面距監(jiān)控點2m時,圍巖應(yīng)力達到M-C屈服線而進入塑性狀態(tài),即掌子面失穩(wěn).而新意法開挖過程中,球應(yīng)力和剪應(yīng)力都在逐漸增大,其應(yīng)力點超過了M-C屈服線,但圍巖應(yīng)力點并沒有完全進入塑性階段.由圖9可知,新意法施工時掌子面更穩(wěn)定.

3.3 拱頂沉降分析

沿隧道軸線方向每2m選取一個截面,各截面拱頂?shù)淖罱K豎向位移如圖10所示,同時選取沿隧道軸線方向上的拱頂特征點的沉降值進行對比分析如表3所示.

圖9 下部掌子面應(yīng)力變化Fig.9 Face stress of lower face

圖10 拱頂沉降曲線Fig.10 Result of vault settlement

施工工法拱頂特征點位置/mmY=0mY=6mY=12mY=18mY=26mY=28mY=30m臺階法-19.3-17.38-15.65-14.45-12.53-12.40-10.88新意法-13.10-11.88-10.66-10.31-8.81-8.36-3.62減小比率/%32.132.131.928.729.732.666.7

由圖10和表3可知:兩種工法開挖隧道完工后,其拱頂豎向位移量均沿著隧道軸線方向逐漸減小.但新意法施工各斷面的拱頂沉降量均比臺階法開挖時小,且變化趨勢更平緩.此外,新意法開挖時拱頂?shù)淖畲筘Q向位移比臺階法的最大豎向位移小了32.1%,且新意法掌子面拱頂處的彈性收斂較臺階法減小了66.7%.由此可見,新意法采用全斷面開挖,減少了施工工序,使隧道圍巖更快地形成閉合圈,同時結(jié)合其超前加固工藝技術(shù),有效地抑制了隧道拱頂?shù)某两?

3.4 掌子面加固區(qū)長度的影響

以新意法掌子面中部擠出位移作為指標(biāo)對加固區(qū)效應(yīng)進行分析,掌子面中部的擠出位移在不同加固區(qū)長度下的變化情況如圖11所示.

由圖11可知:1) 隨著加固區(qū)長度的增加,掌子面的縱向位移逐漸減小,且其變化幅度也在逐漸減小,其中加固區(qū)長度在8 m以內(nèi)時,縱向位移變化幅度較大,當(dāng)超過12 m時,縱向位移曲線變化趨于平緩.對于本隧道出口段,加固區(qū)比較適宜的長度是1倍開挖跨度(12 m),當(dāng)超過這個距離時,加固區(qū)的擴大對于掌子面穩(wěn)定性影響相對較弱;2) 野豬山隧道出口段處于較為軟弱的地層,開挖過程中掌子面易發(fā)生失穩(wěn)破壞,而通過植入全粘結(jié)型玻璃纖維錨桿并加壓注漿的施工工藝,改變了開挖面處超前核心土的物理力學(xué)性質(zhì),在一定程度上增強了掌子面的穩(wěn)定性.

圖11 加固區(qū)長度對縱向位移影響Fig.11 The effect of reinforcement length on the face stability

4 實測數(shù)據(jù)與模擬結(jié)果對比分析

新意法施工現(xiàn)場YK60+280斷面(Y=5 m)和YK60+285斷面(Y=0 m)實測的拱頂沉降隨開挖時間變化曲線如圖12所示.由圖12可知:由于280斷面初始監(jiān)測晚,使得拱頂沉降實測值受到較大的影響,280斷面拱頂沉降的實測值比285斷面的實測值小了近一半.由此可知,初始監(jiān)測時刻對監(jiān)測數(shù)據(jù)量值有較大的影響,且285斷面拱頂沉降的實測值更接近實際收斂變形.其中,實測拱頂沉降值在開挖后5~10 d內(nèi)出現(xiàn)回縮現(xiàn)象是因為掌子面停滯施工時降雨導(dǎo)致淺埋的土體軟化、掌子面剝落造成部分圍巖應(yīng)力釋放,而YK60+280斷面在開挖25 d左右時拱頂沉降值出現(xiàn)回縮現(xiàn)象,是由于施工現(xiàn)場監(jiān)測的后視棱鏡松動導(dǎo)致的.

圖12 拱頂沉降監(jiān)測結(jié)果Fig.12 Monitoring result of vault settlement

隧道圍巖的豎向變形包括了預(yù)收斂和收斂變形,該監(jiān)測項目只能對圍巖的收斂變形數(shù)據(jù)進行采集.由285斷面(Y=0 m)的監(jiān)測結(jié)果可知,隧道最終拱頂沉降約為11.88 mm.由數(shù)值仿真模擬的計算結(jié)果可知,掌子面植入錨桿并注漿加固后,掌子面處拱頂沉降為3.62 mm,285斷面(Y=0 m)的模擬沉降為13.10 mm,其中圍巖預(yù)收斂變形為3.62 mm,收斂變形為9.48 mm左右,模擬值與實測值相差較小,模擬值比實測值小20.2%.由于現(xiàn)場施工過程影響因素較多,模擬值與實測值雖存在一定偏差,但基本吻合.現(xiàn)場監(jiān)測數(shù)據(jù)較好地驗證了該數(shù)值仿真模型的可行性.

5 結(jié) 論

通過野豬山隧道出口段數(shù)值模擬,對比分析新意法和臺階法下隧道掌子面應(yīng)力變化和拱頂?shù)某两担Y(jié)果表明:野豬山隧道臺階法施工時,盡管開挖循環(huán)進尺較小,但掌子面擠出位移較大,其中部極易發(fā)生塑性變形,進行圍巖周邊加固的單一技術(shù)措施并不能有效控制掌子面的變形;臺階法開挖時,掌子面前方2 m范圍內(nèi)的土體穩(wěn)定性問題突出,而新意法的超前加固措施使掌子面前方土體強度發(fā)生改變,有效抑制了掌子面失穩(wěn);新意法通過對掌子面植入錨桿并注漿加固,改變了超前核心土的強度和剛度,其中加固區(qū)比較適宜的長度是1倍開挖跨度;通過對比掌子面穩(wěn)定性、擠出位移和拱頂沉降量,表明在淺埋軟弱地層等特殊工程環(huán)境下新意法施工相比于臺階法更為合理.

[1] 姚軍,王國才,陳祥林.地應(yīng)力釋放對隧道圍巖穩(wěn)定性影響的研究[J].浙江工業(yè)大學(xué)學(xué)報,2010,38(6):629-632.

[2] 王克忠,程青云,王玉培,等.粉質(zhì)砂性土地基中盾構(gòu)區(qū)間隧道開挖過程數(shù)值計算研究[J].浙江工業(yè)大學(xué)學(xué)報,2012,40(5):587-590.

[3] 周藝,何川,鄒育麟,等.破碎千枚巖隧道施工工法比選試驗研究[J].巖石力學(xué)與工程學(xué)報,2013(3):537-548.

[4] 趙錄學(xué).關(guān)于新意法隧道設(shè)計的幾點建議[J].現(xiàn)代隧道技術(shù),2012(1):50-52.

[5] 喻軍,劉松玉.隧道洞口邊坡變形控制與數(shù)值分析[J].浙江工業(yè)大學(xué)學(xué)報,2012,40(1):101-105.

[6] TONON F. ADECO full-face tunnel excavation of two 260 m2tubes in clays with sub-horizontal jet-grouting under minimal urban cover[J]. Tunnelling and underground space technology, 2011, 26(2): 253-266.

[7] TONON F. Sequential excavation, NATM and ADECO: what they have in common and how they differ[J]. Tunnelling and underground space technology, 2010, 25(3): 245-265.

[8] 潘昌實.隧道力學(xué)數(shù)值方法[M].北京:中國鐵道出版社,1995.

[9] 羅祿森,王明年,郭軍.淺埋黃土隧道破壞模式的探討[J].現(xiàn)代隧道技術(shù),2008(4):28-31.

[10] 扈世民,張頂立,李鵬飛,等.大斷面黃土隧道初期支護適應(yīng)性研究[J].巖土力學(xué),2011(S2):660-665.

[11] 周藝.隧道掌子面穩(wěn)定性分析及其控制技術(shù)研究[D].成都:西南交通大學(xué),2010.

[12] YOO C, SHIN H K. Deformation behaviour of tunnel face reinforced with longitudinal pipes-laboratory and numerical investigation[J]. Tunnelling and underground space technology, 2003, 18(4): 303-319.

[13] 陳明祥.彈塑性力學(xué)[M].北京:科學(xué)出版社,2007.

(責(zé)任編輯:劉 巖)

An application study on the ADECO-RS approach in the construction of Yezhu mountain tunnel

WANG Kezhong1, ZHU Haibin1, KONG Lingmin2, LIANG Qidong2, ZHANG Haijun1

(1.College of Civil Engineering and Architecture, Zhejiang University of Technology, Hangzhou 310014, China; 2.Watersupply Management Office of Rizhao, Rizhao 276800, China)

In order to study the applicability of the ADECO-RS approach in tunnel construction under complicated geological conditions, numerical simulations are used to study the variations of the longitudinal deformation and the stress of faces. The spatial-temporal variations of the vault subsidence along the axial direction are predicted with the ADECO-RS approach and the bench method for the outlet section in the mountain tunnel. The effects of the excessive displacement and instability on the faces between different methods of excavation and temporary reinforcement are analyzed. This indicates that the middle part of the face is liable to lose stability and the ADECO-RS approach can control the deformations effectively. This study is of theoretical and practical significance for the construction of similar tunnels.

ADECO-RS approach; bench method; heading face; numerical modeling

2016-09-27

國家自然科學(xué)基金面上資助項目(51679215);清華大學(xué)水沙科學(xué)與水利水電工程國家重點實驗室開放基金(sklhse-2014-C-02)

王克忠(1965—),男,山東冠縣人,教授,研究方向為隧道工程,E-mail:wkz@zjut.edu.cn.

U459.2

A

1006-4303(2017)03-0237-06

猜你喜歡
圍巖變形施工
土木工程施工技術(shù)創(chuàng)新探討
談詩的變形
中華詩詞(2020年1期)2020-09-21 09:24:52
隧道開挖圍巖穩(wěn)定性分析
“我”的變形計
軟弱破碎圍巖隧道初期支護大變形治理技術(shù)
江西建材(2018年4期)2018-04-10 12:37:22
例談拼圖與整式變形
會變形的餅
土木工程施工實習(xí)的探討與實踐
扶貧村里施工忙
河南電力(2016年5期)2016-02-06 02:11:34
采空側(cè)巷道圍巖加固與巷道底臌的防治
主站蜘蛛池模板: 久久精品国产91久久综合麻豆自制| 日韩欧美在线观看| 波多野结衣无码视频在线观看| 亚洲精品动漫| 色亚洲成人| 婷婷99视频精品全部在线观看| 老司机午夜精品视频你懂的| 日本尹人综合香蕉在线观看| 本亚洲精品网站| 特黄日韩免费一区二区三区| 欧美另类精品一区二区三区| 欧美a在线看| 日本高清免费不卡视频| 亚洲另类国产欧美一区二区| 久久精品无码一区二区日韩免费| 一本大道视频精品人妻| 色婷婷综合激情视频免费看 | 日韩精品成人在线| 99免费在线观看视频| 人妻丰满熟妇av五码区| 在线看片中文字幕| 一级成人a毛片免费播放| 五月天久久婷婷| 国产精品一区在线观看你懂的| 国产色爱av资源综合区| 国产人免费人成免费视频| 亚洲人网站| 国产精品亚欧美一区二区| 国产一区二区三区精品欧美日韩| a级毛片免费在线观看| 国产成人av一区二区三区| www.91中文字幕| 国产精品免费p区| 伊人久久大香线蕉影院| 国产高清在线精品一区二区三区| 国内精品视频| 扒开粉嫩的小缝隙喷白浆视频| 2020久久国产综合精品swag| 国产麻豆另类AV| 日本三级欧美三级| 久草中文网| 欧美日韩国产系列在线观看| 亚洲青涩在线| 国产免费久久精品44| 试看120秒男女啪啪免费| 国产91精品久久| 九九香蕉视频| 国产一区二区精品福利| 免费人成视频在线观看网站| 动漫精品啪啪一区二区三区| 午夜无码一区二区三区在线app| 又爽又大又光又色的午夜视频| 在线精品亚洲国产| 久久99精品久久久久纯品| 福利姬国产精品一区在线| yy6080理论大片一级久久| 天堂在线视频精品| 国产精品第三页在线看| 国产亚洲精品自在线| 久久女人网| 国产成人91精品免费网址在线| 国产精品3p视频| 国产 日韩 欧美 第二页| 亚洲精品免费网站| 亚洲国产日韩欧美在线| 福利在线一区| 国产成人91精品| 欧美精品啪啪一区二区三区| 久久特级毛片| 这里只有精品在线| 国产午夜一级毛片| 日韩福利视频导航| 久久久亚洲色| 国内精品视频| 国产精品香蕉在线观看不卡| 亚洲国产亚洲综合在线尤物| 久久无码高潮喷水| 国产区在线观看视频| 五月婷婷导航| 99视频精品全国免费品| 久久午夜夜伦鲁鲁片不卡| 免费三A级毛片视频|