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

    雙向八車道連拱隧道圍巖力學參數(shù)反演分析

    2011-11-06 13:13:42崔志盛
    鉆探工程 2011年5期
    關鍵詞:導坑力學反演

    崔志盛,金 磊,趙 凱

    (1.山西平陽高速公路建設管理處,山西陽泉 045100;2.同濟大學地下建筑與工程系,上海 200092)

    雙向八車道連拱隧道圍巖力學參數(shù)反演分析

    崔志盛1,金 磊2,趙 凱1

    (1.山西平陽高速公路建設管理處,山西陽泉 045100;2.同濟大學地下建筑與工程系,上海 200092)

    在采用雙側(cè)壁導坑法施工的雙向八車道特大斷面連拱隧道中,施工步繁多,臨時支護設置的時間長,隧道全斷面的變形量測只能在臨時支護拆除后進行。由于全斷面變形數(shù)據(jù)獲取得較晚,故較難將其用于圍巖力學參數(shù)的反演。將有限元計算和BP神經(jīng)網(wǎng)絡技術相結(jié)合,并在有限元計算過程中考慮實際的施工步,建立起所有臨時支護拆除之前這一施工狀態(tài)下導坑的變形量與圍巖力學參數(shù)之間的非線性映射關系,并通過對應狀態(tài)下實測的導坑變形值反演了圍巖的力學參數(shù)。將反演的結(jié)果用于正分析驗算,驗證了該方法是可行的。

    隧道工程;雙連拱隧道;位移反分析;有限元;BP神經(jīng)網(wǎng)絡

    0 引言

    由于巖土介質(zhì)的復雜性和不可預見性,巖土材料的物理力學參數(shù)的確定一直是該領域的難題。數(shù)值計算方法的發(fā)展大大提高了各種巖土工程應力、位移和破壞分析的效率,但是在此之前除了要確定合理的計算模型和本構(gòu)關系外,還必須確定合理的輸入?yún)?shù),否則數(shù)值計算的結(jié)果會與實際有較大偏差甚至是產(chǎn)生錯誤的結(jié)果。

    通過現(xiàn)場實測資料反分析巖土材料參數(shù)往往比室內(nèi)的試驗結(jié)果更為可靠,因為室內(nèi)試驗加卸載的過程不可能與巖土體的形成和開挖過程相同,所得的參數(shù)在用于穩(wěn)定性評價的時候存在應力歷史和途徑的差別[1]。結(jié)合監(jiān)控量測,通過反分析技術獲取巖土體的物理力學參數(shù),進而將其應用于正分析,評價圍巖穩(wěn)定性和確定合理的開挖支護方案,并對工程結(jié)構(gòu)進行更準確的安全性評價和分析,就顯得十分重要。位移反分析是目前巖土工程反分析問題中應用最多的一種方法。在隧道工程中,它是利用在隧道開挖和支護過程中量測到的圍巖位移反演圍巖物理力學參數(shù)。由于在位移反分析中輸入的圍巖位移值是圍巖支護體系在開挖影響的空間域和時間域中力學變化最直接、最明顯、綜合性的宏觀體現(xiàn)[2],因此位移反分析所得的參數(shù)從整體上反應了巖體復雜性、施工方法以及隧道工程結(jié)構(gòu)的影響。

    位移反分析的求解方法主要可分為 3類:

    (1)位移逆反分析法,即直接建立基本控制方程的逆方程,通過最小二乘解由量測位移反求待定參數(shù)。Kavanagh[3]在有限元的基礎上,給出了根據(jù)實測的位移或應變逆反分析巖土材料力學參數(shù)的方法; S.Sakurai和 K.Takeuchi[4]以有限元計算為基礎,將巖土介質(zhì)視為線彈性和各向同性材料,給出了位移逆反分析初始地應力和楊氏模量的計算方法。但是由于量測數(shù)據(jù)往往具有離散型,計算可能無法收斂。

    (2)位移正反分析法,即在可能的上下界范圍內(nèi),首先給定一個參數(shù)的試探值,通過迭代運算和誤差函數(shù)的優(yōu)化,求得反演參數(shù)的最佳值。賈超等[5]應用位移正反分析了某隧道圍巖力學參數(shù);朱合華等[6]應用位移正反分析法,用上一個施工步監(jiān)測的位移反演所需計算參數(shù),并以此作為輸入?yún)?shù)用于預測下一個施工步產(chǎn)生的圍巖位移。但是對于大型工程而言,計算范圍大,節(jié)點多,本構(gòu)關系復雜,計算量巨大,正反分析的迭代過程往往較為耗時[6]。

    (3)基于統(tǒng)計理論考慮先驗誤差信息的“貝葉斯方法”和“卡曼濾波法”。蔣樹屏等采用卡爾曼濾波有限元法[7],根據(jù)量測的位移等數(shù)據(jù)反演得到了圍巖內(nèi)部剪應變和塑性區(qū)的歷時變化估計。

    將數(shù)值計算與BP神經(jīng)網(wǎng)絡技術耦合的位移反分析法,實質(zhì)上是采用了位移的正反分析的思路。但是,有別于一般的位移正反分析,該方法是利用BP神經(jīng)網(wǎng)絡強大的學習、存儲和計算能力以及較強的容錯特性[8],通過較少的數(shù)值計算的結(jié)果,建立起位移和待反演參數(shù)之間的高度非線性的映射關系,去代替剩余的數(shù)值計算,以減少計算量,節(jié)約時間。文獻[2、9~11]利用該方法反演了圍巖的力學參數(shù),但是在數(shù)值計算的過程中均未進行詳細的施工步模擬;而被用于反演計算的隧道實測變形是與具體的施工狀態(tài)相對應的。另外,目前鮮有關于采用雙側(cè)壁導坑法施工的雙向八車道特大斷面連拱隧道的圍巖力學參數(shù)的反演的報導。由于施工步繁多,臨時支護設置的時間長,隧道全斷面的變形量測只能在臨時支護拆除后進行,變形數(shù)據(jù)獲取得較晚,故較難將其用于圍巖力學參數(shù)的反演。因此,利用臨時支護拆除前導坑的最終變形量測數(shù)據(jù)進行反演,并在數(shù)值計算過程中進行詳細的施工步模擬,是一項有意義的探索工作。

    1 BP神經(jīng)網(wǎng)絡的結(jié)構(gòu)和工作原理

    由 Rumelhart和 McCelland[12]首先提出的 BP (Back Propagation)神經(jīng)網(wǎng)絡,是一種按誤差逆?zhèn)鞑ニ惴ㄓ柧毜亩鄬忧梆伨W(wǎng)絡,是目前應用最廣泛的神經(jīng)網(wǎng)絡模型之一。BP神經(jīng)網(wǎng)絡模型拓撲結(jié)構(gòu)包括輸入層、隱含層和輸出層。

    Robert Hecht2nielson[13]證明了對任何在閉區(qū)間內(nèi)的連續(xù)函數(shù),都可以用含一個隱含層的 BP網(wǎng)絡來逼近,因而一個三層的BP網(wǎng)絡可以完成任意的 n維到m維的映照。圖 1是三層BP神經(jīng)網(wǎng)絡的結(jié)構(gòu)示意圖。其中的輸入、輸出變換層是為了進行歸一化處理而添加的兩層。

    圖 1 BP神經(jīng)網(wǎng)絡結(jié)構(gòu)示意圖

    BP網(wǎng)絡模型處理信息的基本原理是:輸入信號Xi通過隱含層節(jié)點作用于輸出層節(jié)點,經(jīng)過非線性變換,產(chǎn)生輸出信號 Yk。其過程可表示為:

    式中:Xi——輸入層輸入量,i=1,2,…,n;Hj——隱含層輸出量,j=1,2,…,r,r的確定原則見式 (5); Yk——輸出層輸出量,k=1,2,…,m;Wij——輸入層節(jié)點與隱含層節(jié)點的聯(lián)接權(quán)值;Tjk——隱含層節(jié)點與輸出層節(jié)點的聯(lián)接權(quán)值;θi、θj——神經(jīng)單元閾; f——傳遞函數(shù)。常用的傳遞函數(shù)類型有:Log-sig2 moid型函數(shù)的輸入值可取任意值,輸出值在 0和 1之間;tan-sigmod型傳遞函數(shù)的輸入值可取任意值,輸出值在 -1到 +1之間;線性傳遞函數(shù) purelin的輸入輸出值可取任意值。

    通過調(diào)整Wij、Tjk和θi、θj的值,使誤差沿梯度方向下降,經(jīng)過反復學習訓練,確定與最小誤差相對應的網(wǎng)絡參數(shù)(權(quán)值和閾值),訓練即告停止。此時經(jīng)過訓練的BP神經(jīng)網(wǎng)絡就能對類似樣本的輸入信息自行處理,并輸出使誤差最小的輸出信息。

    2 圍巖力學參數(shù)反演

    2.1 工程概況

    某隧道全長 248 m。隧道設計為雙向八車道特大斷面連拱隧道,單洞寬度為 18.2 m。進口洞底標高 17.311 m,出口洞底標高 18.854 m。單洞標準斷面內(nèi)輪廓面積 (路面以上)左洞為 143.16 m2,右洞為 140.65 m2,含仰拱面積為 171.06 m2。隧道中墻采用 3層復合式曲中墻,整體開挖跨度 39.62 m,高度 11.84 m。

    隧道位于低山丘陵區(qū),地表覆蓋薄層殘坡積土,進出口基巖風化層較厚,場區(qū)植被發(fā)育。地形稍起伏,洞身最高點海拔 60.0 m。隧道進出口自然斜坡穩(wěn)定,天然坡度為 10°~30°。隧道圍巖為花崗巖、花崗斑巖,為較硬巖~堅硬巖,局部為閃長玢巖巖脈,隧道區(qū)巖石風化較強烈,且受巖脈侵入的影響,巖體較破碎。隧道進出口以Ⅴ級圍巖為主,洞身以Ⅲ~Ⅳ級為主。

    2.2 模型建立

    以該隧道位于Ⅴ級圍巖中的 Z24+660斷面為例進行圍巖力學參數(shù)的反演。依據(jù)隧道開挖的影響范圍,模型尺寸為:上邊界至地表自由面;下邊界至洞底距離不小于 3倍洞高;左右邊界距洞口距離不小于 6倍隧道跨度。因此取模型長度為 300 m,模型下邊界距隧道底部 40 m,上邊界為地表面。

    由于隧道圍巖應力場復雜,其在施工過程中的破壞形式也具有多樣性,如受拉破裂、脆性剪切破壞、塑性破壞等,這些隧道圍巖破壞形式一般認為是塑性破壞。因此,在有限元計算中選用理想彈塑性本構(gòu)模型、德魯克 -普拉格 (Drucker-Prager)屈服準則。隧道施工過程通過軟件提供的單元“開挖”、“回填”、“復制線性材料”等命令來實現(xiàn)。模型單元采用四節(jié)點四邊形單元,共劃分 6032個單元。

    邊界約束條件:采用二維平面應變模型,模型的左、右邊界均受到水平方向的位移約束,下部邊界受豎向位移約束,地表則為自由邊界。模型網(wǎng)格劃分情況如圖 2、圖 3所示。

    圖2 Z24+660斷面有限元計算模型

    圖3 Z24+660斷面有限元計算模型局部詳圖

    此外,在建模過程中,圍巖、中墻采用實體單元模擬;初期支護采用梁單元模擬。

    由于只進行圍巖力學參數(shù)(位移模量 E、泊松比μ、粘聚力 c、內(nèi)摩擦角φ)的反演,因此圍巖的物理參數(shù)(即重度γ)按照《羅漢山隧道工程地質(zhì)勘察報告》選取,Ⅴ級圍巖重度γ=22 kN/m3。支護結(jié)構(gòu)計算參數(shù)見表1。

    表1 Z24+660斷面支護結(jié)構(gòu)參數(shù)表

    2.3 施工過程模擬方案

    Ⅴ級圍巖采用雙側(cè)壁導坑法施工。由于在導坑的位移量測過程中不澆注二襯與仰拱,因此模擬開挖與支護的過程分為如下 11步 (見圖 4):(1)開挖中導坑,中導坑初期支護、澆筑中墻;(2)開挖左(右)洞外側(cè)導坑上臺階,施作初期支護、臨時支護; (3)開挖左(右)洞外側(cè)導坑下臺階,施作初期支護、臨時支護;(4)開挖左洞內(nèi)側(cè)導坑上臺階,施作初期支護、臨時支護;(5)開挖左洞內(nèi)側(cè)導坑下臺階,施作初期支護、臨時支護;(6)開挖右洞內(nèi)側(cè)導坑上臺階,施作初期支護、臨時支護;(7)開挖右洞內(nèi)側(cè)導坑下臺階,施作初期支護、臨時支護;(8)開挖左洞中部導坑上臺階,施作初期支護、臨時支護;(9)開挖左洞中部導坑下臺階,施作初期支護;(10)開挖右洞中部導坑上臺階,施作初期支護、臨時支護; (11)開挖右洞中部導坑下臺階,施作初期支護。圖4中加粗的部分表示初期支護或臨時支護。

    另外,通過正分析試算,確定出應力釋放系數(shù)為:開挖后至初期支護、臨時支護施作前為 0.3,初期支護、臨時支護后為 0.7。

    圖4 雙側(cè)壁導坑法施工步

    限于篇幅,僅給出第 4步和第 11步完成后的有限元模型局部詳圖,見圖 5。

    圖5 有限元模型局部詳圖(施工步4、11)

    2.4 導坑測點與測線選取

    圖 6、圖 7是雙洞八車道特大斷面連拱隧道位移量測的測點與測線的布置情況。考慮到左右洞均為外側(cè)導坑上臺階施工最早,故選擇以這兩個導坑的拱頂沉降和水平收斂,即 G1、G3垂直沉降和B1C1、B3C3水平收斂的實測數(shù)據(jù)進行反演。

    圖6 拱頂沉降測點布置圖

    圖7 周邊收斂測線布置圖

    2.5 BP神經(jīng)網(wǎng)絡訓練樣本確定

    2.5.1 圍巖力學參數(shù)取值水平

    根據(jù)《公路隧道設計規(guī)范》(JTG D70-2004)與《錨桿噴射混凝土支護技術規(guī)范》(GB 50086-2001)推薦的各級圍巖力學參數(shù)的取值范圍,確定出圍巖力學參數(shù)取值水平如表 2。

    表2 Ⅴ級圍巖力學參數(shù)取值水平

    2.5.2 BP神經(jīng)網(wǎng)絡訓練樣本

    在已知圍巖的力學參數(shù)取值水平的情況下,結(jié)合四因素四水平正交表 L16(44)確定出不同圍巖力學參數(shù)的組合方式,將其輸入建好的模型,可得到不同圍巖力學參數(shù)組合情況下,第 11個施工步完成后相應的導坑位移值,見表 3。從而各組圍巖力學參數(shù)和相應的導坑位移即組成了 BP神經(jīng)網(wǎng)絡的訓練樣本,其中以導坑位移值作為輸入樣本,圍巖力學參數(shù)作為輸出樣本,即可進行神經(jīng)網(wǎng)絡訓練。

    表3 BP神經(jīng)網(wǎng)絡訓練樣本

    2.6 BP神經(jīng)網(wǎng)絡反演圍巖力學參數(shù)

    神經(jīng)網(wǎng)絡的結(jié)構(gòu)計算參數(shù)對最終的訓練結(jié)果有很大的影響。以表 3數(shù)值計算的導坑位移值為實際輸入,以與之相應的圍巖力學參數(shù)為實際輸出,分別對其進行歸一化處理并將處理結(jié)果轉(zhuǎn)置后分別作為輸入層與輸出層,這是因為考慮到 BP網(wǎng)絡輸入層節(jié)點數(shù)不至于過多以及輸出矩陣和輸入矩陣應具有相同的列數(shù)。其中歸一化過程按照下式計算:

    式中:u——各圍巖力學參數(shù)組合情況下導坑位移值;umin、umax——分別為各圍巖力學參數(shù)組合情況下導坑位移最小值、最大值。

    經(jīng)過這樣的處理后神經(jīng)網(wǎng)絡輸入層、輸出層節(jié)點單元數(shù)均為 4,每個節(jié)點 (神經(jīng)元)為 1×16的矩陣。采用 newff函數(shù)進行訓練,輸入層與隱含層、隱含層與輸出層之間的傳遞用 tansig函數(shù)與 logsig函數(shù)。初步擬定訓練的誤差精度 E=1×10-15。隱含層節(jié)點數(shù)按照下式確定:

    式中:r——隱含層節(jié)點數(shù),r≥log2n;n、m——分別為輸入層與輸出層節(jié)點數(shù);C——1~10之間的常數(shù)。

    經(jīng)過反復調(diào)試,當取隱含層節(jié)點數(shù) n1=13、學習速率η=0.4、訓練步數(shù) t=4000時,訓練的效果較好。訓練結(jié)束的時候,共經(jīng)歷了 4000步運算,此時訓練誤差為 5.29×10-15,符合要求。訓練過程的誤差曲線見圖8。整個過程借助于MATLAB7.0完成。網(wǎng)絡訓練好以后,以表 4中實測的位移值(考慮了 30%的測前位移損失[3])為實際輸入,經(jīng)過歸一化后即可反演Ⅴ級圍巖的力學參數(shù),反演結(jié)果見表 5。再將反演的結(jié)果輸入建立的模型進行正分析計算,則可得導坑位移值,繼而得到其與實測值的相對誤差,見表 4。

    圖8 BP神經(jīng)網(wǎng)絡訓練誤差曲線

    表4 導坑位移實測值與模擬值比較表

    表5 Ⅴ級圍巖力學參數(shù)反演結(jié)果

    由表 4可知,將反演的結(jié)果用于正分析時收到了較好的效果,模擬位移值與實測位移值的平均相對誤差為 8.46%。因此反演的結(jié)果可以用于數(shù)值計算。

    3 結(jié)語

    利用有限元數(shù)值計算和BP神經(jīng)網(wǎng)絡相結(jié)合的位移反分析法,在有限元計算過程中進行了詳細的施工步模擬,根據(jù)雙側(cè)壁導坑法施工過程中臨時支護拆除前這一施工狀態(tài)下導坑的變形量測數(shù)據(jù),反演了雙向八車道特大斷面連拱隧道圍巖的力學參數(shù),并將反演的結(jié)果用于數(shù)值正分析驗算,驗證了反演結(jié)果符合實際情況,其精度能夠滿足工程要求。因此,對于施工步復雜、施工時間跨度長的特大斷面連拱隧道,根據(jù)導坑的變形量測數(shù)據(jù)反演圍巖的力學參數(shù)是可行的,它克服了由于全斷面變形量測數(shù)據(jù)獲取得較晚而難以被用于反演的缺點,能及時反饋信息,指導下一步設計和施工。

    [1] 張勇,張子新,華安增.高速公路隧道的圍巖位移特性[J].中南公路工程,2001,26(3):12-14.

    [2] 李曉紅,靳曉光,亢會明,等.隧道位移智能化反分析及其應用[J].地下空間,2001,21(4):299-304.

    [3] Kavanagh K.Experiment versus analysis:Computational tech2 niques for the description of static material response[J].Interna2 tional Journal forNumericalMethods in Engineering,1973,5:503 -515.

    [4] S.Sakurai,K.Takeuchi.Back analysis of measured displacement of tunnels[J].RockMechanics and Rock Engineering,1983,16: 173-180.

    [5] 賈超,劉寧,肖樹芳.洞室?guī)r體參數(shù)的位移正演反分析[J].巖土力學,2003,24(3):450-454.

    [6] 朱合華,張晨明,王建秀,等.龍山雙連拱隧道動態(tài)位移反分析與預測[J].巖石力學與工程學報,2006,25(1):67-73.

    [7] 蔣樹屏,趙陽.復雜地質(zhì)條件下公路隧道圍巖監(jiān)控量測與非確定性反分析研究 [J].巖石力學與工程學報,2004,23(20): 3460-3464.

    [8] 樊琨,劉宇敏,張艷華.基于人工神經(jīng)網(wǎng)絡的巖土工程力學參數(shù)反分析[J].河海大學學報 (自然科學版),1998,26(4):98 -102.

    [9] 郝哲,萬明富,劉斌,等.韓家?guī)X隧道圍巖物理力學參數(shù)反分析[J].東北大學學報(自然科學版),2005,26(3):300-303.

    [10] 易小明,陳衛(wèi)忠,李術才,等.BP神經(jīng)網(wǎng)絡在分岔隧道位移反分析中的應用[J].巖石力學與工程學報,2006,25(S2):3927 -3932.

    [11] 王向剛,李術才,王剛,等.公路隧道位移智能反演與穩(wěn)定性研究[J].巖土力學,2006,27(S1):322-326.

    [12] RumelhartD.E.,McClelland J.L.Parallel distributed process2 ing[M].M IT press,Cambridge,1986.

    [13] Robert Hecht2nielsen.Theory of the back propagation neural net2 work[A].Proceedings of the International Joint Conference on NeuralNetworks[C].IEEE press,1989.121-125.

    Back Analysis on M echan ical Parameters of Surround ing Rock in a Biarch Tunnel with Eight Lanes

    CU I Zhi-sheng1,JIN Lei2,ZHAO Kai1(1.Office of Pingyang Highway Construction and Management,Yangquan Shanxi 045100, China;2.Department of Geotechnical Engineering,TongjiUniversity,Shanghai 200092,China)

    In a biarch tunnelwith eight lanes using double side driftmethod,it is the deformation of the drifts rather than the whole tunnel section that ismeasured due to the restriction from the temporary linings.Only by the time when all the temporary linings are removed can the back analysis be implemented if we choose to use the whole section’s deformation data to get the mechanical parameters of surrounding rock,which is time consuming.This paper discusses a method which involves both FEM and BP neural network to establish a nonlinear mapping relationship between surrounding rock’s mechanical parameters and drifts’final deformation rather than the whole section’s defor mation with construction sequences taking into consideration,followed by a back analysis on the mechanical parameters using the drifts’in2situ monitoring de for mation data on the basis of the above relationship.There’s a good agreement be tween the in2situ defor mation data and the numerical defor mation data achieved from a normal analysis using the back analyzed parameters.

    tunneling engineering;biarch tunnel;displacement back analysis;FE M;BP neural network

    TU443

    A

    1672-7428(2011)05-0065-05

    2011-01-24

    國家自然科學基金項目(40972178)

    崔志盛(1962-),男(漢族),山西盂縣人,山西平陽高速公路建設管理處總工程師、高級工程師,土木工程專業(yè),從事高速公路建設的工程技術和管理工作,山西省陽泉市,fzhzk@163.com。

    猜你喜歡
    導坑力學反演
    力學
    反演對稱變換在解決平面幾何問題中的應用
    高速公路隧道洞口淺埋段施工中的雙側(cè)壁導坑法
    雙側(cè)壁導坑法隧道下穿既有車站施工方法研究
    弟子規(guī)·余力學文(十)
    快樂語文(2021年11期)2021-07-20 07:41:32
    弟子規(guī)·余力學文(四)
    快樂語文(2020年30期)2021-01-14 01:05:28
    臺階法和中隔壁導坑法開挖在公路隧道施工中的運用
    江西建材(2018年1期)2018-04-04 05:26:24
    基于低頻軟約束的疊前AVA稀疏層反演
    基于自適應遺傳算法的CSAMT一維反演
    力學 等
    av免费观看日本| 国产高清三级在线| 九色成人免费人妻av| 婷婷色综合大香蕉| 3wmmmm亚洲av在线观看| 亚洲欧美成人综合另类久久久| 国产极品天堂在线| 性高湖久久久久久久久免费观看| 国产精品久久久久久久电影| 国内少妇人妻偷人精品xxx网站| 国产免费一级a男人的天堂| 少妇高潮的动态图| 成人毛片60女人毛片免费| 久久精品人妻少妇| 久久精品久久久久久噜噜老黄| 99久久精品热视频| 日本vs欧美在线观看视频 | 亚洲av免费高清在线观看| 狂野欧美白嫩少妇大欣赏| 超碰97精品在线观看| 干丝袜人妻中文字幕| 精品亚洲成a人片在线观看 | 毛片一级片免费看久久久久| 99热这里只有精品一区| 亚洲av二区三区四区| av免费观看日本| 99九九线精品视频在线观看视频| 成人二区视频| 精品少妇黑人巨大在线播放| 少妇高潮的动态图| 边亲边吃奶的免费视频| 成人午夜精彩视频在线观看| 黄色视频在线播放观看不卡| 日本黄大片高清| 成人亚洲精品一区在线观看 | 亚洲精品一二三| 国产精品99久久久久久久久| 国产成人一区二区在线| 精品酒店卫生间| 女性生殖器流出的白浆| 欧美日韩一区二区视频在线观看视频在线| 午夜福利在线在线| 男人舔奶头视频| 性色av一级| 七月丁香在线播放| 亚洲人成网站在线观看播放| 一级毛片我不卡| 久久久久久久久大av| 男女免费视频国产| 人妻一区二区av| 精品视频人人做人人爽| 嘟嘟电影网在线观看| 久久国产亚洲av麻豆专区| 国产精品成人在线| 成人国产麻豆网| 亚洲精品色激情综合| 国产精品精品国产色婷婷| 最后的刺客免费高清国语| 久久精品国产亚洲网站| 国产在视频线精品| 久久国产亚洲av麻豆专区| 久久久久国产精品人妻一区二区| 国产男女超爽视频在线观看| 在线观看av片永久免费下载| 免费高清在线观看视频在线观看| 搡女人真爽免费视频火全软件| 久久久久久久久大av| 国产精品熟女久久久久浪| 人妻一区二区av| av线在线观看网站| 亚洲国产av新网站| av天堂中文字幕网| 99精国产麻豆久久婷婷| 久久99蜜桃精品久久| 美女国产视频在线观看| av一本久久久久| 美女福利国产在线 | 妹子高潮喷水视频| 亚洲第一av免费看| 成人高潮视频无遮挡免费网站| 午夜福利在线观看免费完整高清在| 国产v大片淫在线免费观看| 亚洲精品456在线播放app| 国产一区有黄有色的免费视频| 亚洲精品一区蜜桃| 80岁老熟妇乱子伦牲交| 国产色爽女视频免费观看| 美女中出高潮动态图| 日韩欧美一区视频在线观看 | 少妇的逼水好多| 91精品一卡2卡3卡4卡| 人体艺术视频欧美日本| 亚洲精品久久久久久婷婷小说| 中国三级夫妇交换| 夫妻午夜视频| 99久久中文字幕三级久久日本| 国产国拍精品亚洲av在线观看| 18禁在线播放成人免费| 国产欧美日韩精品一区二区| 精品一区二区三卡| av专区在线播放| 欧美+日韩+精品| 我的老师免费观看完整版| 国产av码专区亚洲av| 26uuu在线亚洲综合色| 日本与韩国留学比较| 精品国产三级普通话版| 中国美白少妇内射xxxbb| 久久亚洲国产成人精品v| 欧美日韩综合久久久久久| 一本久久精品| 九色成人免费人妻av| 在线精品无人区一区二区三 | 国产精品久久久久久久久免| 日韩一区二区三区影片| 亚洲av不卡在线观看| 精品亚洲成国产av| 九草在线视频观看| 波野结衣二区三区在线| www.色视频.com| 人人妻人人添人人爽欧美一区卜 | 国产又色又爽无遮挡免| 成人毛片60女人毛片免费| 国产精品.久久久| 高清黄色对白视频在线免费看 | 亚洲精品自拍成人| 国产欧美另类精品又又久久亚洲欧美| 久久 成人 亚洲| 国产亚洲精品久久久com| 91精品国产九色| 少妇的逼水好多| av福利片在线观看| 午夜日本视频在线| 美女国产视频在线观看| 日日摸夜夜添夜夜爱| 色视频www国产| 欧美变态另类bdsm刘玥| 国产免费福利视频在线观看| 亚洲精品,欧美精品| 亚洲精品色激情综合| 各种免费的搞黄视频| 欧美日韩视频高清一区二区三区二| 嫩草影院入口| 亚洲精品视频女| 欧美日韩精品成人综合77777| 黑人高潮一二区| 超碰av人人做人人爽久久| 春色校园在线视频观看| 99久久综合免费| 免费在线观看成人毛片| 日韩亚洲欧美综合| 亚洲成人手机| 亚洲欧美日韩另类电影网站 | 我要看日韩黄色一级片| 26uuu在线亚洲综合色| 日韩av免费高清视频| 亚洲欧美清纯卡通| 天堂俺去俺来也www色官网| 蜜桃久久精品国产亚洲av| 国产 一区 欧美 日韩| 久久女婷五月综合色啪小说| 亚洲精品,欧美精品| 国产成人免费观看mmmm| 成人漫画全彩无遮挡| 久久6这里有精品| 亚洲精品日本国产第一区| 国产在线一区二区三区精| 亚洲精华国产精华液的使用体验| 国产黄片美女视频| 精品酒店卫生间| 久久人妻熟女aⅴ| 在现免费观看毛片| av国产久精品久网站免费入址| 极品教师在线视频| 亚洲国产日韩一区二区| 交换朋友夫妻互换小说| 亚洲av福利一区| 99久久综合免费| 啦啦啦中文免费视频观看日本| 18禁裸乳无遮挡免费网站照片| 免费在线观看成人毛片| 丝袜喷水一区| 大香蕉久久网| 国产乱来视频区| 老师上课跳d突然被开到最大视频| 日本黄色日本黄色录像| 男的添女的下面高潮视频| 中文欧美无线码| 国产精品人妻久久久影院| 欧美成人a在线观看| 麻豆乱淫一区二区| 国产 一区 欧美 日韩| 岛国毛片在线播放| 中文字幕av成人在线电影| 中文精品一卡2卡3卡4更新| 18禁在线无遮挡免费观看视频| 国产v大片淫在线免费观看| 亚洲久久久国产精品| 亚洲精品乱码久久久久久按摩| 午夜激情久久久久久久| 美女cb高潮喷水在线观看| av在线app专区| .国产精品久久| 成人国产麻豆网| 国产成人免费无遮挡视频| 亚洲av综合色区一区| 亚洲无线观看免费| 国产精品麻豆人妻色哟哟久久| 免费av不卡在线播放| 99热国产这里只有精品6| 日本av手机在线免费观看| 国产视频内射| 日韩电影二区| 人妻制服诱惑在线中文字幕| 国产精品无大码| 有码 亚洲区| 男人和女人高潮做爰伦理| av黄色大香蕉| 国产成人freesex在线| 日本av免费视频播放| 青春草亚洲视频在线观看| 国产美女午夜福利| 亚洲av成人精品一区久久| 免费观看在线日韩| 成人国产av品久久久| 久久久久视频综合| 99久久综合免费| 老师上课跳d突然被开到最大视频| 精品酒店卫生间| 国产爽快片一区二区三区| 丝袜喷水一区| 亚洲中文av在线| 久久久久精品久久久久真实原创| 久久精品久久精品一区二区三区| 国产精品久久久久久久久免| 亚洲中文av在线| 草草在线视频免费看| 男男h啪啪无遮挡| 成年av动漫网址| 免费大片黄手机在线观看| 观看av在线不卡| 91久久精品国产一区二区成人| 国产91av在线免费观看| 一个人看的www免费观看视频| 看非洲黑人一级黄片| 久久国产亚洲av麻豆专区| 精品人妻一区二区三区麻豆| 国产av码专区亚洲av| 国产成人精品婷婷| 51国产日韩欧美| 国产精品不卡视频一区二区| 日本午夜av视频| 亚洲综合精品二区| 免费看不卡的av| 国内揄拍国产精品人妻在线| 久久影院123| 在线亚洲精品国产二区图片欧美 | 人妻少妇偷人精品九色| h视频一区二区三区| 久久99精品国语久久久| 中国国产av一级| 日本欧美国产在线视频| 国产综合精华液| 免费看日本二区| 久久久久国产网址| 自拍欧美九色日韩亚洲蝌蚪91 | 18禁动态无遮挡网站| 久久久久久久久久人人人人人人| 一级片'在线观看视频| 妹子高潮喷水视频| 全区人妻精品视频| 18禁在线播放成人免费| 亚洲精品乱码久久久v下载方式| 国产成人精品久久久久久| 97在线视频观看| 热re99久久精品国产66热6| 人体艺术视频欧美日本| 国产精品免费大片| 99热这里只有精品一区| 天堂中文最新版在线下载| 成人特级av手机在线观看| 欧美3d第一页| 高清午夜精品一区二区三区| 亚洲精品国产色婷婷电影| 99久久中文字幕三级久久日本| 简卡轻食公司| 亚洲精华国产精华液的使用体验| 天堂中文最新版在线下载| 夫妻午夜视频| 久久久久久久国产电影| 你懂的网址亚洲精品在线观看| 国产永久视频网站| 国国产精品蜜臀av免费| 在现免费观看毛片| 97在线人人人人妻| 婷婷色麻豆天堂久久| 久久 成人 亚洲| 国产成人a区在线观看| 在线观看av片永久免费下载| 午夜免费鲁丝| 乱码一卡2卡4卡精品| 免费高清在线观看视频在线观看| 嫩草影院入口| 小蜜桃在线观看免费完整版高清| 久久99精品国语久久久| 欧美另类一区| 婷婷色av中文字幕| 三级经典国产精品| 777米奇影视久久| 能在线免费看毛片的网站| 伊人久久国产一区二区| 黄色一级大片看看| 少妇丰满av| 99热6这里只有精品| 黄色日韩在线| 2018国产大陆天天弄谢| 啦啦啦啦在线视频资源| 人体艺术视频欧美日本| 最近中文字幕2019免费版| 精品久久久久久久久亚洲| 久久国产亚洲av麻豆专区| 亚洲自偷自拍三级| 麻豆成人午夜福利视频| 91狼人影院| 久久99热这里只有精品18| 日韩亚洲欧美综合| 人人妻人人爽人人添夜夜欢视频 | 久久久久久久久久久免费av| 五月玫瑰六月丁香| 国产视频内射| 大陆偷拍与自拍| 日韩av免费高清视频| 国产精品久久久久久久久免| 国产有黄有色有爽视频| 舔av片在线| 极品少妇高潮喷水抽搐| 中文字幕av成人在线电影| 卡戴珊不雅视频在线播放| 青春草视频在线免费观看| 日本与韩国留学比较| 美女cb高潮喷水在线观看| 欧美激情国产日韩精品一区| 黑人高潮一二区| 91精品国产九色| 国产免费福利视频在线观看| 黄色怎么调成土黄色| 亚洲国产毛片av蜜桃av| 2022亚洲国产成人精品| 亚洲欧美精品自产自拍| 有码 亚洲区| 亚洲精品一区蜜桃| 日韩人妻高清精品专区| 麻豆国产97在线/欧美| 成年av动漫网址| 日产精品乱码卡一卡2卡三| 女性被躁到高潮视频| 大陆偷拍与自拍| 国产高清有码在线观看视频| 久久毛片免费看一区二区三区| 国产男人的电影天堂91| 久久人人爽人人片av| 色综合色国产| 国产综合精华液| 国产伦精品一区二区三区四那| 欧美日韩精品成人综合77777| 国产精品熟女久久久久浪| 精品少妇黑人巨大在线播放| 街头女战士在线观看网站| 中文天堂在线官网| 国产黄片视频在线免费观看| 精品亚洲成国产av| 汤姆久久久久久久影院中文字幕| 久久午夜福利片| 精品一区二区三卡| 免费播放大片免费观看视频在线观看| 婷婷色综合www| 亚洲国产日韩一区二区| 美女xxoo啪啪120秒动态图| 日韩电影二区| 亚洲综合色惰| a级毛色黄片| 一级毛片黄色毛片免费观看视频| 精品一区在线观看国产| 国产深夜福利视频在线观看| 交换朋友夫妻互换小说| 天天躁日日操中文字幕| 成人午夜精彩视频在线观看| 十分钟在线观看高清视频www | 97在线视频观看| 又爽又黄a免费视频| 亚洲欧美日韩无卡精品| 亚洲精品日本国产第一区| 网址你懂的国产日韩在线| 国产精品精品国产色婷婷| 青春草视频在线免费观看| 久久久久精品久久久久真实原创| 99热6这里只有精品| 夜夜看夜夜爽夜夜摸| 深爱激情五月婷婷| 成人综合一区亚洲| 日韩制服骚丝袜av| 国产91av在线免费观看| 91aial.com中文字幕在线观看| 51国产日韩欧美| 肉色欧美久久久久久久蜜桃| 国产免费一区二区三区四区乱码| 人妻制服诱惑在线中文字幕| 亚洲av电影在线观看一区二区三区| 亚洲av成人精品一区久久| 亚洲美女视频黄频| 久久久久性生活片| 久久女婷五月综合色啪小说| 久久综合国产亚洲精品| 男人添女人高潮全过程视频| 视频区图区小说| 国产视频内射| 中文字幕av成人在线电影| 性色avwww在线观看| 成人亚洲精品一区在线观看 | 亚洲精品国产色婷婷电影| 免费黄频网站在线观看国产| 日韩中字成人| 麻豆乱淫一区二区| 午夜福利网站1000一区二区三区| 最后的刺客免费高清国语| 亚洲最大成人中文| 国产精品爽爽va在线观看网站| 人妻 亚洲 视频| 丰满人妻一区二区三区视频av| 亚洲精品国产av蜜桃| 欧美精品一区二区免费开放| 久久人人爽人人片av| 国内少妇人妻偷人精品xxx网站| av不卡在线播放| 女性生殖器流出的白浆| 在线观看美女被高潮喷水网站| 久久婷婷青草| 老师上课跳d突然被开到最大视频| 日韩精品有码人妻一区| 欧美日韩亚洲高清精品| 日韩中字成人| 少妇人妻精品综合一区二区| 看十八女毛片水多多多| 亚洲av福利一区| 天堂俺去俺来也www色官网| 网址你懂的国产日韩在线| 精品午夜福利在线看| 亚洲不卡免费看| 亚洲怡红院男人天堂| 18禁动态无遮挡网站| av免费观看日本| 嘟嘟电影网在线观看| 黄片无遮挡物在线观看| 男人爽女人下面视频在线观看| 国产高清不卡午夜福利| 一区二区三区乱码不卡18| 免费大片黄手机在线观看| 亚洲内射少妇av| 熟妇人妻不卡中文字幕| 亚洲精品aⅴ在线观看| 免费观看无遮挡的男女| 成人一区二区视频在线观看| 日本午夜av视频| 中文乱码字字幕精品一区二区三区| 99re6热这里在线精品视频| 搡女人真爽免费视频火全软件| 成年人午夜在线观看视频| 一个人看视频在线观看www免费| 精品久久久久久电影网| 国产v大片淫在线免费观看| 精品午夜福利在线看| 亚洲成人手机| 国产91av在线免费观看| 少妇丰满av| 国产欧美日韩一区二区三区在线 | 成人二区视频| 久久热精品热| 高清不卡的av网站| 中文字幕久久专区| 婷婷色综合大香蕉| 少妇熟女欧美另类| 色网站视频免费| 丰满迷人的少妇在线观看| videos熟女内射| 亚洲精品自拍成人| 久久午夜福利片| 国产精品不卡视频一区二区| 国产成人免费观看mmmm| 国产精品一二三区在线看| 国产精品一区二区在线观看99| 国产精品免费大片| 成年女人在线观看亚洲视频| 亚州av有码| 久久精品国产a三级三级三级| 男女啪啪激烈高潮av片| 国产免费一区二区三区四区乱码| 亚洲人成网站在线观看播放| 草草在线视频免费看| 免费在线观看成人毛片| 全区人妻精品视频| 少妇裸体淫交视频免费看高清| 在线观看一区二区三区激情| 亚洲熟女精品中文字幕| 久久久久网色| 亚洲成色77777| 一级二级三级毛片免费看| 亚洲av成人精品一区久久| 日韩中字成人| 午夜激情福利司机影院| 成人亚洲精品一区在线观看 | 日韩亚洲欧美综合| 成人毛片60女人毛片免费| 大香蕉97超碰在线| 精品国产一区二区三区久久久樱花 | 国产爱豆传媒在线观看| 久久人人爽人人片av| 亚洲欧美日韩无卡精品| 亚洲国产成人一精品久久久| 亚洲美女视频黄频| 久久久a久久爽久久v久久| 黑人高潮一二区| 黄色欧美视频在线观看| 丰满乱子伦码专区| 国产av一区二区精品久久 | 亚洲欧美成人综合另类久久久| 美女cb高潮喷水在线观看| 日本黄大片高清| 午夜激情福利司机影院| 国产精品偷伦视频观看了| 成年免费大片在线观看| 亚洲成人中文字幕在线播放| 久久久久久人妻| 99久久中文字幕三级久久日本| 国产精品麻豆人妻色哟哟久久| xxx大片免费视频| 久久毛片免费看一区二区三区| 观看美女的网站| 老司机影院毛片| 亚洲国产色片| 中文欧美无线码| 老熟女久久久| 国产亚洲91精品色在线| 国产伦理片在线播放av一区| 视频中文字幕在线观看| 亚洲无线观看免费| 欧美zozozo另类| 免费av不卡在线播放| 五月天丁香电影| 交换朋友夫妻互换小说| 久久精品久久久久久久性| 欧美亚洲 丝袜 人妻 在线| 国产69精品久久久久777片| 国产精品偷伦视频观看了| 韩国高清视频一区二区三区| 国产精品福利在线免费观看| 国产精品一二三区在线看| av国产免费在线观看| 纯流量卡能插随身wifi吗| 最新中文字幕久久久久| 欧美成人午夜免费资源| 一本—道久久a久久精品蜜桃钙片| av专区在线播放| 777米奇影视久久| 亚洲三级黄色毛片| 丰满人妻一区二区三区视频av| 日日啪夜夜爽| 亚洲成人中文字幕在线播放| 色视频www国产| 日本与韩国留学比较| 99热全是精品| 午夜福利在线在线| 最近2019中文字幕mv第一页| 我的老师免费观看完整版| 亚洲四区av| 精品熟女少妇av免费看| 亚洲精品国产av成人精品| 老司机影院毛片| av国产久精品久网站免费入址| 五月开心婷婷网| 亚洲欧美日韩卡通动漫| 精品熟女少妇av免费看| 久久久久性生活片| 亚洲av成人精品一二三区| 成人美女网站在线观看视频| 国产高清国产精品国产三级 | 国产免费福利视频在线观看| 亚洲成人手机| 人妻制服诱惑在线中文字幕| 97热精品久久久久久| 午夜福利高清视频| 国产中年淑女户外野战色| 亚洲天堂av无毛| 一级爰片在线观看| 欧美一区二区亚洲| 亚洲成人av在线免费| 我要看黄色一级片免费的| 91午夜精品亚洲一区二区三区| 欧美成人精品欧美一级黄| 一级片'在线观看视频| 国产欧美亚洲国产| 97热精品久久久久久| 日韩一本色道免费dvd| 精品99又大又爽又粗少妇毛片| 日韩av免费高清视频| 精品酒店卫生间| av女优亚洲男人天堂| 少妇丰满av| 日韩一区二区视频免费看| 91久久精品国产一区二区成人| 三级国产精品欧美在线观看| 免费黄色在线免费观看| 日韩精品有码人妻一区| 亚洲第一av免费看| 日韩成人av中文字幕在线观看| 欧美区成人在线视频|