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

    基坑開(kāi)挖引起下臥地鐵隧道變形的簡(jiǎn)化計(jì)算方法

    2021-04-02 05:59:32邱居濤歐孝奪龍團(tuán)元侯凱文賴(lài)增任
    鐵道學(xué)報(bào) 2021年10期
    關(guān)鍵詞:有限元變形模型

    江 杰,邱居濤,歐孝奪,龍團(tuán)元,侯凱文,賴(lài)增任

    (1.廣西大學(xué) 土木建筑工程學(xué)院, 廣西 南寧 530004;2.廣西大學(xué) 工程防災(zāi)與結(jié)構(gòu)安全教育部重點(diǎn)實(shí)驗(yàn)室,廣西 南寧 530004;3.廣西大學(xué) 廣西防災(zāi)減災(zāi)與工程安全重點(diǎn)實(shí)驗(yàn)室,廣西 南寧 530004)

    近年來(lái),隨著城市地鐵建設(shè)的大規(guī)模發(fā)展,在既有地鐵線(xiàn)路上方建造的商業(yè)建筑也隨之增多?;邮┕ねǔO扔谥黧w結(jié)構(gòu)施工,基坑開(kāi)挖難免會(huì)對(duì)周?chē)馏w產(chǎn)生擾動(dòng),導(dǎo)致下臥隧道縱向不均勻變形、地下水滲漏、軌道脫離等。因此,確保地鐵隧道結(jié)構(gòu)穩(wěn)定和運(yùn)營(yíng)安全尤為重要。

    為了研究下臥隧道在基坑開(kāi)挖卸荷下的變形響應(yīng),在試驗(yàn)研究方面,張玉偉等[1]、姚愛(ài)軍等[2]利用模型試驗(yàn)研究了上方基坑開(kāi)挖卸荷-加載作用下地鐵盾構(gòu)隧道的變形特征及圍土壓力分布規(guī)律。對(duì)于計(jì)算分析方法而言通常有兩種,第一種是數(shù)值模擬,Dolezalova[3]、Xiao等[4]、Huang等[5]和鄭剛等[6]使用FEM模型分析基坑開(kāi)挖引起的隧道變形規(guī)律,并針對(duì)不同的影響因素進(jìn)行了模擬。但該方法建模復(fù)雜且計(jì)算耗時(shí)。第二種為兩階段分析法,第一階段根據(jù)解析解得到基坑開(kāi)挖卸荷引起的土體附加應(yīng)力,第二階段假設(shè)基于彈性均質(zhì)地基模型計(jì)算由于附加應(yīng)力導(dǎo)致的隧道變形,相關(guān)學(xué)者對(duì)此已經(jīng)進(jìn)行了大量的研究[7-12]。但既有研究均假定施工場(chǎng)地為均質(zhì)土體,而在實(shí)際施工過(guò)程中,不同土體之間差異性較大,土體的成層性不可忽略。此外,土體并非由一系列獨(dú)立彈簧所組成,在外荷載作用下,土層相互之間存在剪切變形,以往研究將隧道擱置于Winkler地基模型中,僅考慮地基基床系數(shù)對(duì)隧道變形的影響是存在缺陷的。因此,提出一種不但能體現(xiàn)土體分層特性,而且能反映土體單元之間剪切作用的簡(jiǎn)化計(jì)算方法更具現(xiàn)實(shí)意義。

    本文基于彈性層狀體系基本解,提出了考慮分層效應(yīng)的基坑開(kāi)挖擾動(dòng)地層附加應(yīng)力的計(jì)算方法,改變了過(guò)去研究該問(wèn)題時(shí)考慮均質(zhì)地基的單一現(xiàn)狀。首先以彈性層狀體系應(yīng)力解取代傳統(tǒng)的Mindlin應(yīng)力解[13]計(jì)算開(kāi)挖卸荷土體擾動(dòng)作用于隧道的附加應(yīng)力,通過(guò)引入剪切參數(shù)求解出Winkler-Pasternak地基模型中下臥隧道的縱向變形響應(yīng)。最后將有限元數(shù)值模擬與工程監(jiān)測(cè)數(shù)據(jù)進(jìn)行對(duì)比,證明本文所提出簡(jiǎn)化方法的正確性。

    1 兩階段分析法

    1.1 彈性層狀體系應(yīng)力與位移基本解

    n層土體構(gòu)成的彈性地基模型受軸對(duì)稱(chēng)垂直荷載的作用簡(jiǎn)圖見(jiàn)圖1。圖1中,p(r)為軸對(duì)稱(chēng)荷載;δ為荷載半徑;Δhi(i=1,2,…,n)為第i層體系厚度;Ei、μi分別為第i層體系的彈性模量和泊松比。

    圖1 軸對(duì)稱(chēng)荷載施加于多層地基表面示意圖

    參照Burmister彈性層狀理論解[14-16],推導(dǎo)多層彈性體系中任意深度z處的應(yīng)力和位移解表達(dá)式

    (1)

    (2)

    (3)

    (4)

    (5)

    矩陣L和K定義為

    (6)

    式中:z為計(jì)算點(diǎn)深度;hi為第i層土體底部到地表面的距離,Φ(ξ,z)稱(chēng)為傳遞矩陣,定義為

    (7)

    傳遞矩陣Φ各元素表達(dá)式為

    (8)

    式中:E和μ為對(duì)應(yīng)于層狀體系中每一層土的彈性模量和泊松比。

    1.1.1 表面受集中荷載作用的解析解

    當(dāng)多層地基表面作用一集中荷載p時(shí),p(r)的表達(dá)式為

    (9)

    式中:δ(r)為Dirac函數(shù),定義為

    (10)

    結(jié)合式(5)、式(9)和式(10),可以得到p(r)經(jīng)過(guò)Hankel積分變換后的表達(dá)式為

    (11)

    將式(11)代入式(2)和式(4)中,可得集中荷載作用下,與集中荷載水平距離為r處的地基表面的位移和地基內(nèi)任意深度z的豎向附加應(yīng)力,表達(dá)式為

    (12)

    (13)

    1.1.2 表面受矩形均布荷載作用的解析解

    設(shè)基坑的長(zhǎng)度和寬度分別為a和b,作用在坑底的垂直均布荷載為p0,見(jiàn)圖2。

    圖2 層狀地基矩形荷載分布

    (14)

    p0可按下式計(jì)算

    (15)

    式中:γi、hi對(duì)應(yīng)于第i層土底部以上土層的單位重度和厚度。

    利用式(12)、式(13)對(duì)矩形面積積分,可得到矩形均布垂直荷載作用下,荷載作用面任意點(diǎn)B(x0,y0,0)處的位移和荷載作用面以下任意點(diǎn)B'(x0,y0,z)的豎向附加應(yīng)力表達(dá)式

    (16)

    (17)

    式中:Ω為矩形均布垂直荷載的面積。

    1.2 隧道在開(kāi)挖卸荷作用下的變形分析

    基坑與盾構(gòu)隧道相互作用的簡(jiǎn)化模型見(jiàn)圖3,將盾構(gòu)隧道簡(jiǎn)化成Euler-Bernoulli長(zhǎng)梁,假定隧道與周?chē)馏w完全接觸并且土層為各向同性的彈性體,地基土體的非線(xiàn)性變化不納入計(jì)算考慮范圍。Winkler和Winkler-Pasternak地基模型的隧道變形的平衡微分方程表達(dá)式為

    圖3 盾構(gòu)隧道受基坑開(kāi)挖影響示意圖

    (18)

    式中:EI為盾構(gòu)隧道的抗彎剛度;k為地基基床系數(shù);Gp為土體的剪切模量;δ(x)為隧道的撓度;D為盾構(gòu)隧道的外徑。

    因微分方程中存在四階導(dǎo)數(shù),根據(jù)有限差分定理Δmfk=hmf(m)(ξ)可對(duì)導(dǎo)數(shù)進(jìn)行降階,表達(dá)式為

    (19)

    將隧道劃分為n+5個(gè)節(jié)點(diǎn),每個(gè)節(jié)點(diǎn)長(zhǎng)為l,所受附加應(yīng)力為σz(i),豎向位移為δi,將式(19)代入式(18)并整理成矩陣形式為

    (20)

    式中:Q為作用于隧道的附加應(yīng)力并可由式(17)獲得;δ為隧道的縱向位移;ψ1、ψ2、ψ3為系數(shù)矩陣。

    假設(shè)隧道足夠長(zhǎng)且10倍基坑開(kāi)挖尺寸之外隧道兩端剪力和彎矩為零[18],由此可知邊界條件表達(dá)式為

    (21)

    根據(jù)上述邊界條件,進(jìn)一步推導(dǎo)出系數(shù)矩陣ψ1、ψ2、ψ3的表達(dá)式為

    (22)

    (23)

    ψ2為(n+1)×(n+1)階單位矩陣。

    將ψ1、ψ2、ψ3的表達(dá)式代入式(20)并用Matlab進(jìn)行編程求解,便可得到Winkler和Winkler-Pasternak地基模型下隧道的隆起值δ,分別對(duì)δ進(jìn)行二階、三階求導(dǎo)便能求出隧道的彎矩M與剪力Fs,公式為

    (24)

    選取合適的基床系數(shù)k與剪切模量Gp是進(jìn)行合理計(jì)算的前提。根據(jù)Klar等[19]的建議,本文計(jì)算采用2倍Vesic基床系數(shù)值

    (25)

    式中:Es為隧道所處土層彈性模量。

    對(duì)于剪切模量Gp的取值,通常采用簡(jiǎn)化彈性空間法[20],即

    (26)

    式中:Gs為土層的剪切模量;Hp為剪切層的厚度并取[21]Hp=2.5D。

    2 算例驗(yàn)證

    2.1 與均質(zhì)地基中附加應(yīng)力的對(duì)比

    某基坑開(kāi)挖平面尺寸為32 m×16 m,開(kāi)挖深度為8 m,場(chǎng)地共有四層土,具體參數(shù)見(jiàn)表1。

    表1 施工場(chǎng)地土體參數(shù)

    為分析地基土體的非均質(zhì)性對(duì)基坑開(kāi)挖擾動(dòng)地層附加應(yīng)力的影響,選取基坑底面中心以下15 m處平行于基坑短邊的縱向不同節(jié)點(diǎn),將均質(zhì)地基與分層地基兩種情況土體豎向附加應(yīng)力進(jìn)行了對(duì)比,其中,均質(zhì)地基的土體模量按Poulos等[22]提出的加權(quán)平均值計(jì)算,經(jīng)折算后均質(zhì)土體彈性模量取101 MPa,泊松比取0.26。

    兩種不同地基下的附加應(yīng)力分布見(jiàn)圖4,規(guī)律大體相似。兩者的附加應(yīng)力均在x=0處達(dá)到最大,在距基坑中心12 m外均質(zhì)地基中附加應(yīng)力值要大于分層地基的附加應(yīng)力值,而0~12 m內(nèi)分層地基附加應(yīng)力值大于均質(zhì)地基附加應(yīng)力值。由于土體彈性模量自上而下依次增大,土質(zhì)由軟變硬將導(dǎo)致地基產(chǎn)生應(yīng)力集中現(xiàn)象,相反,土質(zhì)由硬變軟將產(chǎn)生應(yīng)力擴(kuò)散現(xiàn)象。本節(jié)土體附加應(yīng)力曲線(xiàn)分布規(guī)律體現(xiàn)了這一土力學(xué)結(jié)論。由于開(kāi)挖卸荷具有對(duì)稱(chēng)性,以基坑中心為劃分點(diǎn),選取一側(cè)的土體附加應(yīng)力值進(jìn)行比較分析,見(jiàn)表2。通過(guò)比較可知,兩者之間最小應(yīng)力偏差為2.4%,最大應(yīng)力偏差可達(dá)67.2%,其中x=0處應(yīng)力偏差為19.8%。目前計(jì)算外荷載作用下土體附加應(yīng)力值主要采用基于均質(zhì)地基的Boussinesq解和Mindlin解,當(dāng)土層參數(shù)相差不大時(shí),采用Boussinesq解和Mindlin解計(jì)算的土體附加應(yīng)力可近似等同于彈性層狀體系解;當(dāng)土層參數(shù)相差較大時(shí),Boussinesq解和Mindlin解的計(jì)算結(jié)果則與彈性層狀體系解明顯不同。因此在理論分析中,考慮地基土體彈性模量和泊松比的變化對(duì)于更準(zhǔn)確計(jì)算土體附加應(yīng)力是具有實(shí)際意義的。

    圖4 均質(zhì)地基與分層地基附加應(yīng)力對(duì)比

    表2 不同地基模型附加應(yīng)力比較

    2.2 與有限元模型的對(duì)比

    某基坑開(kāi)挖施工恰好位于既有盾構(gòu)隧道上方,基坑開(kāi)挖尺寸為8 m×8 m×6 m(長(zhǎng)×寬×深)。施工場(chǎng)地共有四層土,各層土體的具體參數(shù)見(jiàn)表3。盾構(gòu)隧道位于地表以下20 m且隧道縱軸平行于與基坑長(zhǎng)邊,外徑為6.2 m。隧道等效抗彎剛度(EI)eq=7.8×107kN·m2,基床系數(shù)k=4 260 kN/m3,剪切模量Gp=4.257 3×104kN/m。

    表3 模型計(jì)算參數(shù)

    算例采用有限差分軟件FLAC3D進(jìn)行驗(yàn)證??紤]有限元計(jì)算中邊界效應(yīng)的影響,將模型尺寸定義為150 m×150 m×60 m(長(zhǎng)×寬×高),沿著側(cè)壁設(shè)置四道地連墻用來(lái)抵抗基坑側(cè)壁土體的水平應(yīng)力,每道地連墻厚1 m,埋深10 m。為了與理論計(jì)算中涉及的參數(shù)條件保持一致,在有限元模型中隧道采用梁?jiǎn)卧M(jìn)行模擬,周?chē)馏w采用線(xiàn)彈性模型進(jìn)行模擬,并假設(shè)隧道與土層滿(mǎn)足變形協(xié)調(diào)條件。有限元網(wǎng)格劃分見(jiàn)圖5,在整個(gè)模型區(qū)域中,網(wǎng)格個(gè)數(shù)共計(jì)135 218個(gè),節(jié)點(diǎn)共有148 800個(gè),以基坑中心為原點(diǎn),截取位于區(qū)間x=[-75 m,75 m]的隧道進(jìn)行研究。

    圖5 有限元網(wǎng)格劃分示意圖

    采用Winkler模型、Winkler-Pasternak模型與有限元數(shù)值模擬得到的隧道隆起量對(duì)比見(jiàn)圖6。由圖6可見(jiàn),三種方法中采用Winkler方法得到的數(shù)值解遠(yuǎn)大于其余兩者,原因在于Winkler模型忽略了地基彈簧之間的相互作用。此外,在距離基坑中心約50 m處,隧道開(kāi)始產(chǎn)生縱向位移。Winkler模型、解析理論得到的隧道隆起值與有限元解的對(duì)比見(jiàn)表4。采用解析理論計(jì)算得到的隧道縱向最大位移為2.02 mm,數(shù)值模擬為2.25 mm。在開(kāi)挖影響范圍內(nèi),解析理論與有限元結(jié)果的位移偏差都在40%以?xún)?nèi),遠(yuǎn)低于Winkler模型。兩種方法的計(jì)算結(jié)果分布規(guī)律基本一致,數(shù)值上存在微小差別可能原因在于基坑開(kāi)挖土體擾動(dòng)導(dǎo)致既有隧道相對(duì)剛度的減小,本文理論算法中并沒(méi)有考慮這方面的影響。

    圖6 隧道隆起值對(duì)比曲線(xiàn)

    表4 Winkler模型、解析理論與有限元結(jié)果比較(隧道隆起)

    隧道在附加應(yīng)力作用下采用三種計(jì)算方法得到的彎矩分布圖見(jiàn)圖7。由圖7可知,正彎矩最大值位于基坑中心正下方處,此時(shí)襯砌管片處于最不利受彎拉狀態(tài),最大負(fù)彎矩產(chǎn)生于距基坑中心兩側(cè)約25 m的位置。盾構(gòu)隧道是由一系列襯砌管片拼裝而成的結(jié)構(gòu)體,當(dāng)作用于其上的彎矩超出隧道極限抗彎能力時(shí),襯砌管片環(huán)縫會(huì)大幅度張開(kāi)從而影響隧道結(jié)構(gòu)的穩(wěn)定。Winkler模型、解析理論得到的隧道彎矩值與有限元解的對(duì)比見(jiàn)表5。采用解析理論得到的彎矩偏差均值為57.4%,小于Winkler模型得到的彎矩偏差均值148.8%,與有限元計(jì)算結(jié)果相比,Winkler模型計(jì)算得到的最大正、負(fù)彎矩值高出79.7%、105.3%,總體來(lái)說(shuō),采用本文方法與有限元算法得到的彎矩值吻合較好,曲線(xiàn)分布規(guī)律也基本一致。

    圖7 隧道彎矩對(duì)比曲線(xiàn)

    表5 Winkler模型、解析理論與有限元結(jié)果比較(隧道彎矩)

    采用本文理論解、Winkler模型解以及有限元模擬得到的隧道縱向剪力對(duì)比曲線(xiàn)見(jiàn)圖8。其中本文計(jì)算結(jié)果與有限元計(jì)算結(jié)果數(shù)值吻合較好,剪力最大(小)值皆發(fā)生于x=±9 m處。為了維護(hù)隧道的結(jié)構(gòu)穩(wěn)定,可在剪力最大處采取加固接頭螺栓的辦法。Winkler模型、解析理論得到的隧道剪力值與有限元解的對(duì)比見(jiàn)表6。在x=±9 m處Winkler模型計(jì)算出的最大(小)剪力值要超出數(shù)值模擬的82.7%。與有限元結(jié)果相比,采用解析理論得到的剪力偏差均值為62.3%,小于Winkler模型得到的剪力偏差均值108.1%。

    表6 Winkler模型、解析理論與有限元結(jié)果比較(隧道剪力)

    圖8 隧道縱向剪力對(duì)比曲線(xiàn)

    根據(jù)上述分析可知,在預(yù)測(cè)由于開(kāi)挖卸荷引起下臥隧道變形中,本文提供了一種快速而簡(jiǎn)便的方法。

    2.3 與監(jiān)測(cè)數(shù)據(jù)對(duì)比

    上海東方路地下過(guò)街通道采用明挖法施工,上海軌道交通2號(hào)線(xiàn)的一段隧道正好位于基坑正下方并與基坑成45°夾角,見(jiàn)圖9。為了簡(jiǎn)化分析過(guò)程,將基坑視為26 m×18 m(長(zhǎng)×寬)的矩形,各層土體參數(shù)見(jiàn)表7。基坑開(kāi)挖深度為6.5 m且隧道頂部距離坑底僅為2.76 m,地鐵隧道使用盾構(gòu)法施工,采用鋼螺栓連接的預(yù)制節(jié)段環(huán)作為永久襯砌,厚度為0.35 m。隧道軸線(xiàn)位于地表以下12.36 m,地鐵2號(hào)線(xiàn)隧道外徑為6.2 m。隧道等效抗彎剛度取1.087×108kN·m2?;拥撞繚仓?.88 m厚的混凝土底板。考慮地層改良后地基剛度的增加,理論計(jì)算中隧道所處土層彈性模量取16.01 MPa。

    圖9 基坑與隧道位置關(guān)系示意(單位:m)

    表7 場(chǎng)地土層參數(shù)

    為了驗(yàn)證本文兩階段方法的適用性,將理論計(jì)算結(jié)果與現(xiàn)場(chǎng)監(jiān)測(cè)結(jié)果進(jìn)行了對(duì)比分析,本文理論解計(jì)算得到的隧道隆起值見(jiàn)圖10,并與文獻(xiàn)[23]中現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)以及Winkler模型解進(jìn)行了對(duì)比分析。

    從圖10可知,采用本文方法計(jì)算的隧道隆起值與實(shí)測(cè)數(shù)據(jù)分布規(guī)律基本一致,并且比大多數(shù)實(shí)測(cè)隆起值要大一些,原因在于實(shí)際施工過(guò)程會(huì)對(duì)基坑采取各種加固措施如分層分區(qū)開(kāi)挖、底板澆筑、圍護(hù)樁的打設(shè)等,然而本文理論計(jì)算沒(méi)有考慮這些因素。隧道彎矩對(duì)比曲線(xiàn)見(jiàn)圖11,由圖10、圖11可知,采用本文解析理論得到的隧道隆起和彎矩曲線(xiàn)低于Winkler理論模型得到的曲線(xiàn),意味著如用Winkler模型進(jìn)行隧道縱向變形預(yù)測(cè),需采取過(guò)多的保護(hù)措施來(lái)減小隧道隆起量。此外,本文理論解與隧道實(shí)測(cè)隆起量都滿(mǎn)足地鐵隧道的規(guī)范標(biāo)準(zhǔn),即隧道最大隆起值不超過(guò)20 mm。

    圖10 隧道隆起值對(duì)比曲線(xiàn)

    圖11 隧道彎矩對(duì)比曲線(xiàn)

    3 結(jié)論

    (1)與以往基于均質(zhì)地基的兩階段法不同,在計(jì)算基坑開(kāi)挖卸荷作用于下臥隧道的附加荷載時(shí),本文建立的集中荷載與矩形荷載下彈性層狀體系應(yīng)力解不僅考慮了彈性模量與泊松比隨土體分層變化而不同的情況,還避免了均質(zhì)地基中應(yīng)力擴(kuò)散能力過(guò)大的弊端,與工程實(shí)際情況符合,具有更好的計(jì)算精度。

    (2)將地鐵隧道簡(jiǎn)化為Euler-Bernoulli長(zhǎng)梁架臥于Winkler-Pasternak地基模型中,從而建立梁的撓度平衡微分方程,通過(guò)有限差分法對(duì)方程進(jìn)行降階處理,最終得到隧道在附加應(yīng)力作用下的縱向變形解答。

    (3)在預(yù)測(cè)隧道縱向變形響應(yīng)的方法中,采用Winkler地基模型通常導(dǎo)致計(jì)算結(jié)果偏大,Winkler-Pasternak模型充分考慮了土體彈簧之間的相互作用,與工程實(shí)際更貼近,計(jì)算結(jié)果更具參考性與精確性。此外,基坑中心下方區(qū)域?qū)е滤淼喇a(chǎn)生最大縱向位移且使襯砌管片處于最不利受彎拉狀態(tài),施工中應(yīng)加強(qiáng)這部分區(qū)域的變形監(jiān)測(cè)與保護(hù)力度。

    (4)本文未考慮基坑側(cè)壁土體卸荷作用以及隧道與周?chē)馏w之間的非線(xiàn)性作用,在將來(lái)研究中可進(jìn)一步完善。

    猜你喜歡
    有限元變形模型
    一半模型
    重要模型『一線(xiàn)三等角』
    談詩(shī)的變形
    重尾非線(xiàn)性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    “我”的變形計(jì)
    例談拼圖與整式變形
    會(huì)變形的餅
    3D打印中的模型分割與打包
    磨削淬硬殘余應(yīng)力的有限元分析
    基于SolidWorks的吸嘴支撐臂有限元分析
    国产高清视频在线播放一区| av视频免费观看在线观看| 亚洲av第一区精品v没综合| 大型黄色视频在线免费观看| 国产日韩一区二区三区精品不卡| 十分钟在线观看高清视频www| 日韩 欧美 亚洲 中文字幕| 又黄又粗又硬又大视频| 国产精品免费大片| 欧美日韩一级在线毛片| 精品国产乱码久久久久久男人| 欧美黑人精品巨大| 另类精品久久| 在线观看舔阴道视频| 久久久久久久精品吃奶| 国产视频一区二区在线看| 久久久久久久久免费视频了| 国产黄色免费在线视频| 国产精品二区激情视频| 日韩欧美一区二区三区在线观看 | 18禁国产床啪视频网站| 天堂动漫精品| 国产主播在线观看一区二区| 亚洲第一av免费看| 亚洲欧美色中文字幕在线| 亚洲精品久久成人aⅴ小说| 97在线人人人人妻| 热re99久久精品国产66热6| 黑人巨大精品欧美一区二区mp4| 午夜福利影视在线免费观看| 久久久久久久久免费视频了| 国产单亲对白刺激| 69精品国产乱码久久久| 超色免费av| 亚洲人成电影观看| 这个男人来自地球电影免费观看| 精品国产一区二区三区久久久樱花| 亚洲成av片中文字幕在线观看| 国产精品免费一区二区三区在线 | 少妇 在线观看| 老司机亚洲免费影院| 久久人人爽av亚洲精品天堂| 国产单亲对白刺激| 狠狠狠狠99中文字幕| 国产伦理片在线播放av一区| 黄色 视频免费看| 怎么达到女性高潮| 无人区码免费观看不卡 | 欧美大码av| av福利片在线| 香蕉丝袜av| 久久毛片免费看一区二区三区| 9热在线视频观看99| 精品第一国产精品| 欧美成狂野欧美在线观看| 亚洲人成77777在线视频| 男女无遮挡免费网站观看| 国产视频一区二区在线看| 久久影院123| 精品久久蜜臀av无| 母亲3免费完整高清在线观看| 一进一出好大好爽视频| 免费黄频网站在线观看国产| 欧美 日韩 精品 国产| 日韩免费高清中文字幕av| 久久ye,这里只有精品| 一本色道久久久久久精品综合| 日韩人妻精品一区2区三区| 日韩熟女老妇一区二区性免费视频| 久久中文看片网| 精品亚洲成国产av| 亚洲精品在线美女| 老司机亚洲免费影院| 丰满少妇做爰视频| 精品少妇内射三级| 国产精品一区二区免费欧美| 免费观看a级毛片全部| 最近最新中文字幕大全免费视频| 成人免费观看视频高清| 国产不卡一卡二| 自线自在国产av| 在线观看免费视频日本深夜| 精品午夜福利视频在线观看一区 | 亚洲九九香蕉| 人人妻人人澡人人爽人人夜夜| 大陆偷拍与自拍| 一夜夜www| 国产成人精品无人区| 在线看a的网站| 久久精品成人免费网站| 十八禁人妻一区二区| 满18在线观看网站| 最近最新中文字幕大全免费视频| 露出奶头的视频| 久久国产精品大桥未久av| 黄色视频,在线免费观看| 久久亚洲精品不卡| h视频一区二区三区| h视频一区二区三区| 日韩 欧美 亚洲 中文字幕| 亚洲五月色婷婷综合| 国产成人精品久久二区二区免费| 99国产精品99久久久久| 亚洲午夜理论影院| 99国产精品一区二区三区| 岛国在线观看网站| 久久久久精品国产欧美久久久| tocl精华| 色婷婷av一区二区三区视频| 亚洲国产av影院在线观看| 国产成人av教育| aaaaa片日本免费| 桃花免费在线播放| 欧美成人免费av一区二区三区 | 日本黄色日本黄色录像| 亚洲精品国产区一区二| 欧美日韩精品网址| 蜜桃国产av成人99| 欧美精品啪啪一区二区三区| 亚洲专区字幕在线| 午夜91福利影院| 国产精品 欧美亚洲| www.精华液| 大片电影免费在线观看免费| 国产欧美日韩综合在线一区二区| 日本黄色日本黄色录像| 免费少妇av软件| 热99久久久久精品小说推荐| 欧美激情 高清一区二区三区| 18禁国产床啪视频网站| 国产在线视频一区二区| 国产免费现黄频在线看| 免费观看av网站的网址| 久久久久国内视频| 日韩一卡2卡3卡4卡2021年| av片东京热男人的天堂| 99精品久久久久人妻精品| 嫁个100分男人电影在线观看| 欧美亚洲 丝袜 人妻 在线| 制服诱惑二区| 精品欧美一区二区三区在线| 美女扒开内裤让男人捅视频| 青青草视频在线视频观看| 亚洲国产欧美在线一区| 午夜久久久在线观看| 久久人人爽av亚洲精品天堂| 淫妇啪啪啪对白视频| 少妇的丰满在线观看| 亚洲欧洲日产国产| 久久热在线av| 精品国内亚洲2022精品成人 | 午夜福利,免费看| 91精品三级在线观看| 亚洲av日韩精品久久久久久密| 狠狠婷婷综合久久久久久88av| 亚洲中文av在线| 99国产精品一区二区蜜桃av | 国产一区二区激情短视频| 欧美+亚洲+日韩+国产| 欧美日韩精品网址| 一边摸一边抽搐一进一出视频| 欧美 亚洲 国产 日韩一| 999久久久精品免费观看国产| 乱人伦中国视频| 丝瓜视频免费看黄片| tube8黄色片| av欧美777| 啦啦啦免费观看视频1| 女人被躁到高潮嗷嗷叫费观| 国产精品免费视频内射| 建设人人有责人人尽责人人享有的| 亚洲国产欧美日韩在线播放| 国产精品香港三级国产av潘金莲| 精品一品国产午夜福利视频| 99国产精品免费福利视频| 一二三四社区在线视频社区8| 久久精品人人爽人人爽视色| 亚洲精品美女久久久久99蜜臀| 成年动漫av网址| 亚洲专区中文字幕在线| 欧美 日韩 精品 国产| 精品亚洲成a人片在线观看| 久久精品熟女亚洲av麻豆精品| 欧美老熟妇乱子伦牲交| 性高湖久久久久久久久免费观看| 国产伦理片在线播放av一区| 国产一卡二卡三卡精品| 真人做人爱边吃奶动态| 一本大道久久a久久精品| 久久国产精品人妻蜜桃| 啦啦啦视频在线资源免费观看| 亚洲av成人一区二区三| 欧美日韩国产mv在线观看视频| 午夜免费鲁丝| 精品国产亚洲在线| 国产精品久久电影中文字幕 | 热99久久久久精品小说推荐| 亚洲av美国av| 久久精品人人爽人人爽视色| 黄频高清免费视频| 久久人妻熟女aⅴ| 少妇裸体淫交视频免费看高清 | 十分钟在线观看高清视频www| 亚洲av欧美aⅴ国产| 国产精品九九99| 久久av网站| 国产伦理片在线播放av一区| 欧美在线黄色| 成人国产一区最新在线观看| 亚洲五月色婷婷综合| 激情在线观看视频在线高清 | 色婷婷av一区二区三区视频| 免费av中文字幕在线| 亚洲av国产av综合av卡| 国产一区二区三区综合在线观看| 亚洲中文av在线| av有码第一页| 男女免费视频国产| 精品免费久久久久久久清纯 | 一边摸一边做爽爽视频免费| 亚洲专区国产一区二区| 午夜福利影视在线免费观看| 夜夜骑夜夜射夜夜干| 国产在视频线精品| 色播在线永久视频| 国产精品一区二区精品视频观看| 少妇猛男粗大的猛烈进出视频| 久久午夜综合久久蜜桃| 美女扒开内裤让男人捅视频| 黑人欧美特级aaaaaa片| 一本大道久久a久久精品| 高清欧美精品videossex| 欧美精品啪啪一区二区三区| 一级a爱视频在线免费观看| 每晚都被弄得嗷嗷叫到高潮| videos熟女内射| 男女无遮挡免费网站观看| 18禁国产床啪视频网站| 欧美激情 高清一区二区三区| 国产精品久久久人人做人人爽| 久久久国产欧美日韩av| 韩国精品一区二区三区| 亚洲国产欧美日韩在线播放| 另类精品久久| 男人舔女人的私密视频| 欧美精品高潮呻吟av久久| av在线播放免费不卡| 国产日韩欧美在线精品| 国产又爽黄色视频| 亚洲精品一二三| 十八禁网站网址无遮挡| 纵有疾风起免费观看全集完整版| 精品国产乱码久久久久久男人| 一夜夜www| 午夜精品久久久久久毛片777| 欧美日韩亚洲国产一区二区在线观看 | 男人舔女人的私密视频| 亚洲av成人不卡在线观看播放网| 国产在线观看jvid| 亚洲精品一二三| 欧美日韩中文字幕国产精品一区二区三区 | 嫁个100分男人电影在线观看| 一级片免费观看大全| 久久人妻福利社区极品人妻图片| 久久精品亚洲av国产电影网| 亚洲天堂av无毛| 大码成人一级视频| 成人18禁高潮啪啪吃奶动态图| 香蕉久久夜色| 久久人人97超碰香蕉20202| 久久av网站| 母亲3免费完整高清在线观看| 99热国产这里只有精品6| 婷婷丁香在线五月| 中文字幕高清在线视频| 精品少妇一区二区三区视频日本电影| 精品国产一区二区三区四区第35| 精品久久久久久久毛片微露脸| 又大又爽又粗| 桃花免费在线播放| 男女免费视频国产| tocl精华| 一级片免费观看大全| 国产av又大| 久久午夜综合久久蜜桃| 曰老女人黄片| 欧美在线黄色| 精品国产乱码久久久久久男人| 国产精品偷伦视频观看了| 黑人猛操日本美女一级片| 中文字幕av电影在线播放| a级毛片黄视频| 久久久久久人人人人人| 午夜福利一区二区在线看| 中亚洲国语对白在线视频| 亚洲国产中文字幕在线视频| 我要看黄色一级片免费的| 美女视频免费永久观看网站| 日韩欧美一区二区三区在线观看 | 亚洲自偷自拍图片 自拍| 老熟女久久久| 中文字幕人妻丝袜一区二区| 最近最新免费中文字幕在线| 成人三级做爰电影| av电影中文网址| 99久久人妻综合| 性高湖久久久久久久久免费观看| 国产精品美女特级片免费视频播放器 | xxxhd国产人妻xxx| 激情在线观看视频在线高清 | 精品少妇一区二区三区视频日本电影| 欧美在线黄色| 国内毛片毛片毛片毛片毛片| bbb黄色大片| 最近最新中文字幕大全电影3 | 每晚都被弄得嗷嗷叫到高潮| 老汉色∧v一级毛片| 精品一区二区三卡| 天天添夜夜摸| 欧美精品一区二区大全| av片东京热男人的天堂| 国产精品 欧美亚洲| 丁香六月欧美| 麻豆乱淫一区二区| 久久av网站| 久久人妻熟女aⅴ| 欧美日韩一级在线毛片| 午夜成年电影在线免费观看| 精品卡一卡二卡四卡免费| 国产日韩欧美视频二区| 久久毛片免费看一区二区三区| 夜夜爽天天搞| 最新在线观看一区二区三区| 美女福利国产在线| 每晚都被弄得嗷嗷叫到高潮| 亚洲视频免费观看视频| 天天操日日干夜夜撸| 一级毛片精品| 久久毛片免费看一区二区三区| 男人操女人黄网站| 美女福利国产在线| 亚洲伊人色综图| 精品国产一区二区三区久久久樱花| 免费在线观看黄色视频的| 免费一级毛片在线播放高清视频 | 亚洲一区中文字幕在线| 午夜福利欧美成人| 国产亚洲精品久久久久5区| 最黄视频免费看| xxxhd国产人妻xxx| 超碰97精品在线观看| a级片在线免费高清观看视频| 桃红色精品国产亚洲av| 久久午夜综合久久蜜桃| 中文字幕最新亚洲高清| 久久影院123| 无遮挡黄片免费观看| 老司机影院毛片| 免费在线观看黄色视频的| 色视频在线一区二区三区| 国产视频一区二区在线看| 在线观看免费视频网站a站| 性色av乱码一区二区三区2| 老司机午夜十八禁免费视频| 在线 av 中文字幕| 久久久欧美国产精品| 一区二区三区精品91| 免费不卡黄色视频| 黄色毛片三级朝国网站| 精品国产超薄肉色丝袜足j| 亚洲视频免费观看视频| 久久亚洲真实| 欧美日韩亚洲国产一区二区在线观看 | 成年版毛片免费区| 国产又爽黄色视频| 日本av免费视频播放| 亚洲va日本ⅴa欧美va伊人久久| 精品人妻在线不人妻| 99精品久久久久人妻精品| 高清欧美精品videossex| 精品人妻熟女毛片av久久网站| 国产区一区二久久| 水蜜桃什么品种好| 国产一区二区三区综合在线观看| 在线av久久热| 日韩一区二区三区影片| 在线永久观看黄色视频| 日本av免费视频播放| 欧美日韩国产mv在线观看视频| 免费在线观看黄色视频的| 国产免费现黄频在线看| 精品一区二区三区四区五区乱码| 久久久久久久久免费视频了| 9热在线视频观看99| 日韩欧美一区二区三区在线观看 | 激情在线观看视频在线高清 | 1024香蕉在线观看| 18禁黄网站禁片午夜丰满| 欧美激情久久久久久爽电影 | 国产精品影院久久| 狠狠婷婷综合久久久久久88av| 一级片'在线观看视频| 一二三四在线观看免费中文在| 高清视频免费观看一区二区| 亚洲人成伊人成综合网2020| 国产成+人综合+亚洲专区| 99re在线观看精品视频| 嫁个100分男人电影在线观看| 久久 成人 亚洲| 国产精品久久久av美女十八| 精品福利永久在线观看| 啦啦啦视频在线资源免费观看| 丝袜喷水一区| 久久久久精品人妻al黑| 69精品国产乱码久久久| 亚洲国产欧美在线一区| 女同久久另类99精品国产91| 十八禁网站免费在线| 国产精品亚洲av一区麻豆| 在线观看免费视频日本深夜| 久久精品人人爽人人爽视色| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久免费高清国产稀缺| 国产精品.久久久| 桃花免费在线播放| 高清欧美精品videossex| 18禁国产床啪视频网站| 精品一区二区三区四区五区乱码| 欧美激情久久久久久爽电影 | cao死你这个sao货| 50天的宝宝边吃奶边哭怎么回事| 不卡一级毛片| 国产一区二区三区在线臀色熟女 | 在线亚洲精品国产二区图片欧美| 极品教师在线免费播放| 色尼玛亚洲综合影院| 美女视频免费永久观看网站| 美女高潮到喷水免费观看| 十八禁人妻一区二区| 久久婷婷人人爽人人干人人爱| 亚洲五月婷婷丁香| 欧美日本视频| 香蕉av资源在线| 欧美色欧美亚洲另类二区| 亚洲欧美日韩东京热| av天堂在线播放| 日本a在线网址| 欧美日韩亚洲国产一区二区在线观看| 在线观看66精品国产| 国产高清视频在线观看网站| 三级毛片av免费| 两性午夜刺激爽爽歪歪视频在线观看| 久久亚洲精品不卡| 午夜精品在线福利| 欧美高清成人免费视频www| 一级a爱片免费观看的视频| 欧美xxxx黑人xx丫x性爽| 国内精品久久久久精免费| av天堂在线播放| 亚洲国产欧美一区二区综合| 亚洲在线观看片| 婷婷精品国产亚洲av在线| 免费人成视频x8x8入口观看| 观看免费一级毛片| 国产爱豆传媒在线观看| 51午夜福利影视在线观看| 天堂影院成人在线观看| 悠悠久久av| 成人午夜高清在线视频| 免费观看精品视频网站| 国产又色又爽无遮挡免费看| 久久性视频一级片| 特级一级黄色大片| 麻豆国产97在线/欧美| 国产成人一区二区三区免费视频网站| 欧美日韩精品网址| 成年女人永久免费观看视频| 国产成人av教育| 成人高潮视频无遮挡免费网站| 99精品欧美一区二区三区四区| 很黄的视频免费| 日本熟妇午夜| 九色成人免费人妻av| 亚洲avbb在线观看| 亚洲av成人不卡在线观看播放网| 中国美女看黄片| 真人做人爱边吃奶动态| 久久精品aⅴ一区二区三区四区| 操出白浆在线播放| 国产亚洲精品久久久久久毛片| 搡老妇女老女人老熟妇| 亚洲欧美日韩高清在线视频| 色噜噜av男人的天堂激情| 女警被强在线播放| 小说图片视频综合网站| 亚洲av电影在线进入| 国产久久久一区二区三区| 日本免费一区二区三区高清不卡| 亚洲av熟女| 成人一区二区视频在线观看| 国产高清激情床上av| 欧美3d第一页| 精品电影一区二区在线| 一本久久中文字幕| 免费看日本二区| 每晚都被弄得嗷嗷叫到高潮| 黄片大片在线免费观看| 又黄又爽又免费观看的视频| 少妇裸体淫交视频免费看高清| 国产综合懂色| 亚洲真实伦在线观看| 久久天堂一区二区三区四区| 精品久久蜜臀av无| 久久精品91蜜桃| 国产亚洲精品久久久com| 99久久久亚洲精品蜜臀av| АⅤ资源中文在线天堂| 最新在线观看一区二区三区| 一区二区三区高清视频在线| 成在线人永久免费视频| 国产久久久一区二区三区| 黄色 视频免费看| 亚洲aⅴ乱码一区二区在线播放| 中文字幕av在线有码专区| 精华霜和精华液先用哪个| av女优亚洲男人天堂 | 九九热线精品视视频播放| 看免费av毛片| 久久精品综合一区二区三区| 熟女人妻精品中文字幕| 村上凉子中文字幕在线| 九色国产91popny在线| 极品教师在线免费播放| 国产精品久久久久久久电影 | 少妇的逼水好多| 夜夜爽天天搞| 亚洲 国产 在线| 欧美中文综合在线视频| 19禁男女啪啪无遮挡网站| 国产精品国产高清国产av| 一本一本综合久久| 精华霜和精华液先用哪个| 精品国产乱子伦一区二区三区| 一夜夜www| 国内精品美女久久久久久| 精品久久蜜臀av无| 国产av不卡久久| 又黄又爽又免费观看的视频| 成人av在线播放网站| 少妇丰满av| 美女被艹到高潮喷水动态| 国产精品久久久久久人妻精品电影| 三级国产精品欧美在线观看 | 1024香蕉在线观看| 国产亚洲精品一区二区www| 午夜免费激情av| 国产男靠女视频免费网站| 国产 一区 欧美 日韩| 香蕉丝袜av| 国产一区二区三区视频了| 99久国产av精品| 99在线人妻在线中文字幕| 婷婷亚洲欧美| 一进一出抽搐动态| 久久这里只有精品中国| 可以在线观看毛片的网站| 夜夜爽天天搞| 亚洲精品在线观看二区| 国产av一区在线观看免费| 黄色视频,在线免费观看| 全区人妻精品视频| 少妇的丰满在线观看| 亚洲18禁久久av| 岛国在线观看网站| 黄色片一级片一级黄色片| 他把我摸到了高潮在线观看| 亚洲成av人片在线播放无| 男人的好看免费观看在线视频| 欧美日本亚洲视频在线播放| 久久久精品大字幕| 亚洲七黄色美女视频| 欧美一区二区精品小视频在线| 国产久久久一区二区三区| 亚洲狠狠婷婷综合久久图片| 日本黄色视频三级网站网址| 99久久综合精品五月天人人| 亚洲中文av在线| 精品不卡国产一区二区三区| 欧美黄色淫秽网站| 波多野结衣高清作品| 婷婷亚洲欧美| 国内精品久久久久久久电影| 久久精品国产清高在天天线| 免费在线观看视频国产中文字幕亚洲| 午夜精品久久久久久毛片777| 久久这里只有精品中国| 国产熟女xx| 麻豆国产av国片精品| 久久这里只有精品中国| 国产熟女xx| 久9热在线精品视频| 国内揄拍国产精品人妻在线| 在线观看舔阴道视频| 青草久久国产| 日韩 欧美 亚洲 中文字幕| 免费观看精品视频网站| 色精品久久人妻99蜜桃| 亚洲美女黄片视频| 欧美最黄视频在线播放免费| 色精品久久人妻99蜜桃| 无遮挡黄片免费观看| 一级毛片精品| 俺也久久电影网| 18禁黄网站禁片免费观看直播| 一区二区三区国产精品乱码|