• <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一区在线观看免费| 国产97色在线日韩免费| 99在线视频只有这里精品首页| 国产av又大| 久久九九热精品免费| 1024手机看黄色片| 操出白浆在线播放| 亚洲成国产人片在线观看| 国产精品一区二区免费欧美| 亚洲第一欧美日韩一区二区三区| 99在线视频只有这里精品首页| 国产精品久久久久久人妻精品电影| 亚洲国产精品成人综合色| 欧美激情高清一区二区三区| 国产精品免费一区二区三区在线| 久久天堂一区二区三区四区| 哪里可以看免费的av片| 国产亚洲欧美在线一区二区| 搡老妇女老女人老熟妇| 国产伦人伦偷精品视频| 亚洲国产精品成人综合色| 亚洲av熟女| 亚洲色图 男人天堂 中文字幕| 嫁个100分男人电影在线观看| 一级作爱视频免费观看| 99riav亚洲国产免费| 国产成人av教育| 国产精品 欧美亚洲| 99国产精品99久久久久| 欧美日本视频| 国产伦在线观看视频一区| 色av中文字幕| 国产av一区在线观看免费| 久久人妻av系列| 成人亚洲精品av一区二区| 欧美一级毛片孕妇| 久久久久国产一级毛片高清牌| 亚洲国产精品久久男人天堂| 特大巨黑吊av在线直播 | 又黄又爽又免费观看的视频| 真人做人爱边吃奶动态| 一级毛片女人18水好多| 最近最新中文字幕大全电影3 | 国产精品国产高清国产av| 超碰成人久久| 看片在线看免费视频| 精品国产国语对白av| 一本综合久久免费| 亚洲国产日韩欧美精品在线观看 | 亚洲一区二区三区色噜噜| 国产精品综合久久久久久久免费| 午夜福利在线在线| 午夜福利一区二区在线看| 国产免费av片在线观看野外av| 久久久久久大精品| 禁无遮挡网站| 日本在线视频免费播放| www日本在线高清视频| 精品人妻1区二区| 国产三级黄色录像| 真人一进一出gif抽搐免费| 级片在线观看| 免费无遮挡裸体视频| 亚洲国产精品sss在线观看| 99国产综合亚洲精品| 国产麻豆成人av免费视频| 高潮久久久久久久久久久不卡| 亚洲中文字幕日韩| 国产免费av片在线观看野外av| 国产成人精品无人区| 成人国语在线视频| 亚洲欧美精品综合久久99| 正在播放国产对白刺激| 久久久久免费精品人妻一区二区 | 黄色女人牲交| 免费av毛片视频| 国内毛片毛片毛片毛片毛片| 亚洲欧美一区二区三区黑人| 俺也久久电影网| 亚洲成a人片在线一区二区| 欧美大码av| 精品高清国产在线一区| 日韩中文字幕欧美一区二区| 侵犯人妻中文字幕一二三四区| 又紧又爽又黄一区二区| 久久婷婷成人综合色麻豆| 国产午夜福利久久久久久| 悠悠久久av| 国产免费av片在线观看野外av| 久久久久久大精品| 午夜亚洲福利在线播放| 18禁美女被吸乳视频| 一夜夜www| 在线观看免费日韩欧美大片| 日韩一卡2卡3卡4卡2021年| 人人妻人人看人人澡| 欧美日本亚洲视频在线播放| 91av网站免费观看| 欧美日韩亚洲国产一区二区在线观看| 一级a爱视频在线免费观看| 男人舔奶头视频| 中文字幕另类日韩欧美亚洲嫩草| 在线观看舔阴道视频| 成人一区二区视频在线观看| 亚洲第一电影网av| 日日爽夜夜爽网站| 白带黄色成豆腐渣| av中文乱码字幕在线| 成在线人永久免费视频| 国产成人欧美在线观看| 国产精品98久久久久久宅男小说| 久久精品国产清高在天天线| 别揉我奶头~嗯~啊~动态视频| 成人精品一区二区免费| 亚洲国产毛片av蜜桃av| 18美女黄网站色大片免费观看| 精品第一国产精品| 亚洲全国av大片| 久久久国产精品麻豆| e午夜精品久久久久久久| 日韩国内少妇激情av| 久久精品国产亚洲av香蕉五月| 久久亚洲真实| 中文字幕最新亚洲高清| 99久久久亚洲精品蜜臀av| 视频区欧美日本亚洲| 性欧美人与动物交配| 天天一区二区日本电影三级| 久久中文看片网| 中文字幕人妻熟女乱码| 天天一区二区日本电影三级| 免费在线观看完整版高清| 亚洲人成电影免费在线| 超碰成人久久| 岛国视频午夜一区免费看| 在线观看日韩欧美| 欧美精品啪啪一区二区三区| www日本黄色视频网| 欧美成人一区二区免费高清观看 | 可以在线观看的亚洲视频| 黑人操中国人逼视频| 18禁黄网站禁片免费观看直播| 色精品久久人妻99蜜桃| 日日夜夜操网爽| 免费看a级黄色片| 亚洲七黄色美女视频| 国产高清videossex| 身体一侧抽搐| 国产精品av久久久久免费| 成人亚洲精品av一区二区| 久久精品91蜜桃| 国产一区二区三区在线臀色熟女| 精品午夜福利视频在线观看一区| 黄网站色视频无遮挡免费观看| 美女高潮到喷水免费观看| 少妇的丰满在线观看| 在线观看免费视频日本深夜| 色综合站精品国产| 老司机福利观看| 麻豆成人av在线观看| 99re在线观看精品视频| 精品免费久久久久久久清纯| 日韩欧美在线二视频| 国产一级毛片七仙女欲春2 | 美国免费a级毛片| 久久精品夜夜夜夜夜久久蜜豆 | 桃色一区二区三区在线观看| 久久香蕉精品热| 天堂动漫精品| 国产单亲对白刺激| 妹子高潮喷水视频| 久久中文字幕一级| 国产aⅴ精品一区二区三区波| 欧美一区二区精品小视频在线| 男女之事视频高清在线观看| 两个人看的免费小视频| 波多野结衣av一区二区av| 国产精品久久久久久人妻精品电影| 黄色a级毛片大全视频| av天堂在线播放| 久久天躁狠狠躁夜夜2o2o| 窝窝影院91人妻| 免费看a级黄色片| 中国美女看黄片| 男男h啪啪无遮挡| 国产又爽黄色视频| 午夜精品久久久久久毛片777| 亚洲中文日韩欧美视频| 久热这里只有精品99| 最近在线观看免费完整版| 欧美 亚洲 国产 日韩一| 午夜免费观看网址| 欧美黄色淫秽网站| 精品一区二区三区视频在线观看免费| 精品免费久久久久久久清纯| 高潮久久久久久久久久久不卡| 波多野结衣av一区二区av| 禁无遮挡网站| 天天一区二区日本电影三级| 国产精品一区二区三区四区久久 | 国内少妇人妻偷人精品xxx网站 | 51午夜福利影视在线观看| svipshipincom国产片| 久久精品人妻少妇| 久久精品91无色码中文字幕| 搡老熟女国产l中国老女人| 香蕉av资源在线| 两人在一起打扑克的视频| 在线观看免费视频日本深夜| 日本免费a在线| 欧美久久黑人一区二区| 女人高潮潮喷娇喘18禁视频| 国产av又大| 女性被躁到高潮视频| www.999成人在线观看| 成人国语在线视频| ponron亚洲| 亚洲aⅴ乱码一区二区在线播放 | 国产av在哪里看| 午夜免费激情av| 韩国av一区二区三区四区| 久久精品国产清高在天天线| 久久草成人影院| 妹子高潮喷水视频| 色播亚洲综合网| 中文字幕av电影在线播放| av在线播放免费不卡| 操出白浆在线播放| 黄网站色视频无遮挡免费观看| 制服丝袜大香蕉在线| 50天的宝宝边吃奶边哭怎么回事| 男女下面进入的视频免费午夜 | 国产亚洲av高清不卡| 天堂√8在线中文| 欧美人与性动交α欧美精品济南到| 男女之事视频高清在线观看| 成年免费大片在线观看| 国产伦在线观看视频一区| 国产又爽黄色视频| 午夜福利视频1000在线观看| 日本在线视频免费播放| 欧美日韩精品网址| 99精品欧美一区二区三区四区| 99在线视频只有这里精品首页| 美国免费a级毛片| 他把我摸到了高潮在线观看| 欧美乱妇无乱码| 黄色毛片三级朝国网站| 女人高潮潮喷娇喘18禁视频| 亚洲av熟女| 怎么达到女性高潮| 成人三级黄色视频| 大香蕉久久成人网| 色播在线永久视频| 久久人人精品亚洲av| 人成视频在线观看免费观看| 久久伊人香网站| 国内少妇人妻偷人精品xxx网站 | 午夜日韩欧美国产| 免费看a级黄色片| 亚洲一码二码三码区别大吗| 日本在线视频免费播放| 熟女电影av网| 亚洲熟妇中文字幕五十中出| 久久久国产成人免费| 一级毛片高清免费大全| 日本免费一区二区三区高清不卡| 中国美女看黄片| 欧美 亚洲 国产 日韩一| 亚洲精品美女久久av网站| 亚洲欧美激情综合另类| 国产伦在线观看视频一区| 啦啦啦免费观看视频1| 欧美在线一区亚洲| 国产片内射在线| 动漫黄色视频在线观看| 18禁国产床啪视频网站| 久久性视频一级片| 国产激情欧美一区二区| 成人国产一区最新在线观看| 听说在线观看完整版免费高清| 精品久久蜜臀av无| 亚洲五月色婷婷综合| 午夜福利免费观看在线| 一级毛片女人18水好多| 88av欧美| 国产又黄又爽又无遮挡在线| 在线永久观看黄色视频| 日韩欧美在线二视频| 色在线成人网| 桃色一区二区三区在线观看| 久久欧美精品欧美久久欧美| 久9热在线精品视频| 麻豆成人av在线观看| av视频在线观看入口| 亚洲午夜理论影院| 女性生殖器流出的白浆| 黄色a级毛片大全视频| 欧美午夜高清在线| 欧美国产日韩亚洲一区| www日本黄色视频网| xxxwww97欧美| 此物有八面人人有两片| bbb黄色大片| 妹子高潮喷水视频| 给我免费播放毛片高清在线观看| 黑人巨大精品欧美一区二区mp4| av视频在线观看入口| 91九色精品人成在线观看| 哪里可以看免费的av片| 精品国内亚洲2022精品成人| 一二三四社区在线视频社区8| 中文在线观看免费www的网站 | 免费女性裸体啪啪无遮挡网站| 亚洲av成人一区二区三| 精品久久久久久久毛片微露脸| 夜夜躁狠狠躁天天躁| 欧美午夜高清在线| 日韩欧美在线二视频| 国产真人三级小视频在线观看| av片东京热男人的天堂| 国产亚洲精品第一综合不卡| 久久久久久亚洲精品国产蜜桃av| 伊人久久大香线蕉亚洲五| 久久久国产成人免费| 久久久久久免费高清国产稀缺| 久久精品成人免费网站| 亚洲人成网站在线播放欧美日韩| 中亚洲国语对白在线视频| 又黄又爽又免费观看的视频| 日韩av在线大香蕉| 成人国产综合亚洲| 69av精品久久久久久| 无限看片的www在线观看| 黄色视频,在线免费观看| 国产成人精品久久二区二区免费| 国产一区二区三区在线臀色熟女| 在线观看舔阴道视频| 亚洲av电影在线进入| 国产精品久久电影中文字幕| 在线观看www视频免费| 国产精品自产拍在线观看55亚洲| 一级毛片高清免费大全| 国产一区二区三区在线臀色熟女| 亚洲国产日韩欧美精品在线观看 | 午夜激情av网站| 欧美zozozo另类| 757午夜福利合集在线观看| 男女之事视频高清在线观看| 欧美日韩瑟瑟在线播放| 99国产精品一区二区蜜桃av| 欧美不卡视频在线免费观看 | 午夜成年电影在线免费观看| 精品国产美女av久久久久小说| 国产激情久久老熟女| 久久国产精品男人的天堂亚洲| 少妇的丰满在线观看| 无遮挡黄片免费观看| 欧美黄色淫秽网站| 国产真人三级小视频在线观看| 91麻豆av在线| 日韩欧美一区视频在线观看| 日本在线视频免费播放| 亚洲第一青青草原| 午夜久久久久精精品| 色综合亚洲欧美另类图片| 国产精品 欧美亚洲| 国产av一区二区精品久久| 丰满人妻熟妇乱又伦精品不卡| 午夜福利一区二区在线看| 国产av一区在线观看免费| 日韩欧美一区二区三区在线观看| 九色国产91popny在线| 国产91精品成人一区二区三区| 午夜久久久久精精品| 国产精品,欧美在线| 最新美女视频免费是黄的| 欧美黄色片欧美黄色片| 可以在线观看毛片的网站| 久久天堂一区二区三区四区| 国产97色在线日韩免费| 国产成人欧美在线观看| 无人区码免费观看不卡| 欧美激情久久久久久爽电影| АⅤ资源中文在线天堂| 观看免费一级毛片| 搡老岳熟女国产| 亚洲欧美激情综合另类| 国产亚洲av高清不卡| 听说在线观看完整版免费高清| 99精品在免费线老司机午夜| 久久精品国产清高在天天线| 日韩欧美在线二视频| 美女大奶头视频| 国产日本99.免费观看| 999久久久国产精品视频| 欧美色视频一区免费| 久久精品91蜜桃| 色精品久久人妻99蜜桃| 老熟妇乱子伦视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 怎么达到女性高潮| 国产高清激情床上av| 亚洲成国产人片在线观看| 韩国精品一区二区三区| 男人的好看免费观看在线视频 | 久久中文字幕一级| 观看免费一级毛片| 国产精品香港三级国产av潘金莲| 亚洲欧美激情综合另类| 妹子高潮喷水视频| 欧美成人性av电影在线观看| 免费在线观看日本一区| www日本黄色视频网| 精品熟女少妇八av免费久了| 曰老女人黄片| 天天添夜夜摸| 高潮久久久久久久久久久不卡| 午夜福利一区二区在线看| www国产在线视频色| 神马国产精品三级电影在线观看 | 久久久久免费精品人妻一区二区 | 在线观看舔阴道视频| 性色av乱码一区二区三区2| 国产精品国产高清国产av| 亚洲av熟女| 欧美日韩精品网址| 给我免费播放毛片高清在线观看| www.自偷自拍.com| 婷婷丁香在线五月| 久久天堂一区二区三区四区| 日日夜夜操网爽| 黄色视频,在线免费观看| 夜夜夜夜夜久久久久| 久久久国产欧美日韩av| 丰满人妻熟妇乱又伦精品不卡| 99精品欧美一区二区三区四区| 免费一级毛片在线播放高清视频| 日本免费a在线| 午夜激情av网站| 成人手机av| 久久久久久九九精品二区国产 | 久久久久免费精品人妻一区二区 | 国产伦在线观看视频一区| 69av精品久久久久久| 国产亚洲精品久久久久久毛片| 国产熟女午夜一区二区三区| 国产亚洲精品第一综合不卡| 日本在线视频免费播放| 黄色视频不卡| 成人永久免费在线观看视频| 欧美黑人巨大hd| 亚洲精品国产区一区二| 桃色一区二区三区在线观看| 男女做爰动态图高潮gif福利片| 淫妇啪啪啪对白视频| 色综合站精品国产| 欧美激情久久久久久爽电影| 给我免费播放毛片高清在线观看| 午夜免费成人在线视频| 亚洲国产精品久久男人天堂| 啦啦啦 在线观看视频| av电影中文网址| 日本 欧美在线| 妹子高潮喷水视频| 成年免费大片在线观看| 国产精品久久久人人做人人爽| 久久青草综合色| 欧美性猛交╳xxx乱大交人| 免费无遮挡裸体视频| 午夜精品在线福利| 色在线成人网| 久久午夜综合久久蜜桃| 色综合站精品国产| 欧美不卡视频在线免费观看 | 啦啦啦韩国在线观看视频| av在线天堂中文字幕| 国产精品99久久99久久久不卡| 精品欧美国产一区二区三| 久久久久亚洲av毛片大全| 99久久综合精品五月天人人| 男女床上黄色一级片免费看| 91九色精品人成在线观看| 色精品久久人妻99蜜桃| 久久久久久人人人人人| 狂野欧美激情性xxxx| 国产精品免费一区二区三区在线| 精品国产国语对白av| 欧美黑人巨大hd| 亚洲成人久久性| 最新在线观看一区二区三区| 国产91精品成人一区二区三区| 成人精品一区二区免费| 男女视频在线观看网站免费 | 日韩高清综合在线| 69av精品久久久久久| 又黄又爽又免费观看的视频| 成人午夜高清在线视频 | 久久久久久久午夜电影| 亚洲午夜理论影院| 亚洲av日韩精品久久久久久密| 日日夜夜操网爽| 国产精品日韩av在线免费观看| 中文字幕高清在线视频| 91成年电影在线观看| 国产亚洲精品综合一区在线观看 | 日韩av在线大香蕉| 久热爱精品视频在线9| 亚洲国产欧美一区二区综合| 午夜两性在线视频| 国产av又大| 韩国精品一区二区三区| 欧美色欧美亚洲另类二区| 精品人妻1区二区| 日本撒尿小便嘘嘘汇集6| 自线自在国产av| 午夜精品久久久久久毛片777| 动漫黄色视频在线观看| 在线观看午夜福利视频| 亚洲成人免费电影在线观看| www.www免费av| 亚洲成人免费电影在线观看| 午夜福利视频1000在线观看| 亚洲中文字幕一区二区三区有码在线看 | e午夜精品久久久久久久| 欧美日本亚洲视频在线播放| 99久久精品国产亚洲精品| 亚洲美女黄片视频| 精品午夜福利视频在线观看一区| 中文字幕人成人乱码亚洲影| 色av中文字幕| 日本 av在线| 精品久久久久久久久久久久久 | 久久午夜亚洲精品久久| 不卡一级毛片| 黄色片一级片一级黄色片| 日韩 欧美 亚洲 中文字幕| 国产精品电影一区二区三区| 免费搜索国产男女视频| 韩国精品一区二区三区| 精品久久久久久久久久久久久 | 可以在线观看毛片的网站| 色播亚洲综合网| www.自偷自拍.com| 亚洲天堂国产精品一区在线| 亚洲性夜色夜夜综合| 一区二区三区激情视频| 国产单亲对白刺激| 哪里可以看免费的av片| 国产乱人伦免费视频| 人人妻人人澡人人看| 一区二区日韩欧美中文字幕| 亚洲午夜理论影院| 午夜福利视频1000在线观看| 国产精品久久久av美女十八| 制服人妻中文乱码| 日韩欧美免费精品| 亚洲中文字幕一区二区三区有码在线看 | 丝袜美腿诱惑在线| 动漫黄色视频在线观看| 老司机福利观看| 久久久久久久午夜电影| 侵犯人妻中文字幕一二三四区| x7x7x7水蜜桃| 久久精品国产清高在天天线| 亚洲自偷自拍图片 自拍| 久久香蕉激情| 久久精品aⅴ一区二区三区四区| 精品卡一卡二卡四卡免费| 高清毛片免费观看视频网站| 丁香欧美五月| 可以在线观看毛片的网站| 国产精品久久视频播放| 村上凉子中文字幕在线| 国产99久久九九免费精品| 亚洲五月天丁香| 人人妻,人人澡人人爽秒播| 最近最新中文字幕大全电影3 | 亚洲av成人一区二区三| 亚洲五月婷婷丁香| 午夜福利欧美成人| 成人精品一区二区免费| 男人的好看免费观看在线视频 | 夜夜爽天天搞| 亚洲真实伦在线观看| 午夜免费鲁丝| 亚洲精华国产精华精| 午夜福利18| 成人欧美大片| 18禁裸乳无遮挡免费网站照片 | 18禁黄网站禁片午夜丰满| 黄色女人牲交| 搞女人的毛片| 精品熟女少妇八av免费久了| 久久精品国产清高在天天线| 亚洲va日本ⅴa欧美va伊人久久| 99精品久久久久人妻精品| 免费av毛片视频| 亚洲va日本ⅴa欧美va伊人久久| bbb黄色大片| 欧美黄色片欧美黄色片| 曰老女人黄片| 国产97色在线日韩免费| 黄色视频不卡| 午夜两性在线视频| 亚洲av电影不卡..在线观看| 成在线人永久免费视频| 12—13女人毛片做爰片一| 天堂影院成人在线观看| 性色av乱码一区二区三区2| 自线自在国产av|