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

    突加集中荷載作用下飽和土-深埋圓形隧道襯砌系統(tǒng)的非軸對稱動力響應

    2016-12-21 03:31:50陸劍駿楊驍
    上海大學學報(自然科學版) 2016年5期
    關鍵詞:圓形徑向孔隙

    陸劍駿,楊驍

    (上海大學土木工程系,上?!?00072)

    突加集中荷載作用下飽和土-深埋圓形隧道襯砌系統(tǒng)的非軸對稱動力響應

    陸劍駿,楊驍

    (上海大學土木工程系,上海200072)

    考慮飽和土與深埋圓形襯砌的相互作用,研究了突加集中荷載作用下飽和土-襯砌系統(tǒng)的非軸對稱動力響應.基于Biot理論和彈性理論,采用Laplace變換和Fourier級數(shù),考慮襯砌邊界條件以及襯砌與飽和土交界面處的連續(xù)性條件,在Laplace變換域內求得突加集中荷載作用下飽和土-彈性襯砌耦合系統(tǒng)的位移、應力和孔隙水壓力等的解析表達式.利用Laplace逆變換Crump數(shù)值反演方法得到飽和土-襯砌系統(tǒng)動力響應的數(shù)值解,并分析了土體和襯砌系統(tǒng)的力學、幾何等參數(shù)對系統(tǒng)動力響應的影響.結果表明:5倍隧道襯砌半徑以外處土體的動力響應遠小于隧道附近土體的動力響應;襯砌剛度和厚度對土體位移和應力影響顯著,但對孔隙水壓力影響較??;孔隙水的可壓縮性對土體位移的幅值影響不大,但對應力幅值的影響較為顯著.

    飽和土;彈性襯砌;集中荷載;動力響應;Laplace變換

    隨著城市建設的發(fā)展,土地資源日益緊張,地下空間的開發(fā)利用顯得尤為重要.地下隧道是地下工程的常見結構形式之一,常用于交通、人防和管線.隧道在全壽命運營周期內可能受到動荷載的作用,如機車突加荷載和振動、碰撞以及爆炸等.瞬間的動力作用會對隧道結構和周圍土體產(chǎn)生較大影響,因此研究地下隧道在突加荷載作用下的瞬態(tài)響應不僅具有重要的理論意義,也能夠為隧道安全性的評估提供指導.

    近年來,基于飽和土的Biot理論[1-2]和多孔介質理論[3-4],飽和土中隧道的動力特性和動力響應已有較豐富的研究成果.除深埋隧道動力特性的研究成果[5-10]外,Senjuntichai等[11]研究了不同荷載下飽和彈性介質中圓形隧道的動力響應,楊峻等[12]利用積分變換法給出了飽和土體中圓柱形孔洞動力響應解.將襯砌等效為彈性Flügge殼體,Zakout等[13]研究了彈性土體中圓形隧道動力響應;丁伯陽等[14]基于Green函數(shù)給出了飽和土中矩形隧道承受集中荷載的振動位移;劉干斌等[15]研究了飽和黏彈性土體中隧道開挖的土體應力和位移變化.考慮隧道襯砌的影響,Lu等[16]將襯砌視為片狀彈性體,采用微分求積法數(shù)值研究了飽和彈性土中隧道的動力響應;Gao等[17]解析研究了不同荷載作用下彈性襯砌深埋圓形隧道飽和彈性土的動力響應.將襯砌等效為Flügge殼體,蔡袁強等[18]研究了簡化爆炸荷載作用下圓形隧道飽和土的動力響應.這些工作揭示了飽和土和襯砌系統(tǒng)的動力學行為.然而,上述研究成果主要針對深埋圓形隧道飽和土-襯砌系統(tǒng)的軸對稱動力響應,而關于飽和土-深埋隧道襯砌系統(tǒng)的非軸對稱動力響應研究成果則較少[19],其中非軸對稱變形的動力特性,如振動沿環(huán)向的傳播規(guī)律等有待深入研究.

    基于飽和土的Biot模型,本工作研究了飽和彈性土-深埋圓形隧道襯砌耦合系統(tǒng)在突加集中荷載作用下的非軸對稱動力響應.首先,將問題視為平面應變問題,基于飽和土的Biot理論,通過引入勢函數(shù),利用分離變量法和Laplace變換,得到飽和土在Laplace變換域中的動力響應通解;其次,將隧道襯砌變形分別視為平面應變問題和Flügge殼體,得到襯砌的動力響應通解.在此基礎上,利用隧道襯砌和飽和土接觸面處的連續(xù)性條件及邊界條件,得到Laplace變換域中飽和土-襯砌耦合系統(tǒng)非軸對稱的位移、應力和孔隙水壓力解析解,并利用Laplace逆變換的Crump方法,得到時間域中飽和土-襯砌系統(tǒng)的瞬態(tài)位移、應力和孔隙水壓力響應.在驗證本工作理論和算法的基礎上,數(shù)值研究了飽和彈性土-深埋圓形隧道襯砌耦合系統(tǒng)在突加集中荷載作用下的動力響應,分析了不同隧道模型和材料參數(shù)等對飽和土動力響應的影響.

    1 問題描述及飽和土控制方程求解

    圖1為飽和土中深埋圓形隧道的示意圖.設無限飽和彈性土體中有襯砌厚度為h的無限長圓形隧道承受徑向集中壓力f(θ,t)=p(t)δ(θ)的作用,其中δ(θ)為delta函數(shù).襯砌的內外半徑分別為R1和R2,平均半徑R=(R1+R2)/2,h=R2?R1.假設飽和土-襯砌系統(tǒng)的動力響應為小變形,且土體與彈性襯砌完全接觸,二者間無相對滑移,接觸面處不滲透.根據(jù)圓形隧道幾何和荷載分布特征,可以將飽和土-襯砌系統(tǒng)的動力響應視為非軸對稱平面應變問題.

    圖1 飽和土中的深埋圓形隧道Fig.1 Circular tunnel deeply embedded in saturated soil

    記飽和土平面應變的徑向和環(huán)向位移分別為uSr和uSθ,且uS=(uSr,uSθ,0),而孔隙水相對于土體骨架的徑向和環(huán)向位移分別為wr和wθ,且w=(wr,wθ,0),則根據(jù)飽和土的Biot理論[1-2],飽和彈性土的動力方程為

    式中,λS和μS為飽和土體的表觀Lame常數(shù);α和M分別為反映土骨架和孔隙流體壓縮性的常數(shù);ρf為孔隙水的真實密度,飽和土體的表觀密度為ρ=(1?n)ρS+nρf,其中n為孔隙率, ρS為土骨架的真實密度;m=ρf/n;b=η0/kd為流固相互作用系數(shù),其中η0和kd分別為流體的黏滯系數(shù)和動力滲透系數(shù).

    應力-應變本構關系可以表示為

    式中,σSij(i,j=r,θ)為土體的總應力;δij為Kronecker符號;e和ξ分別為土體和孔隙流體的體積應變,e=divuS,ξ=?divw;εSij(i,j=r,θ)為土骨架的應變;pf為孔隙水壓力.

    并引入如下無量綱參數(shù):

    可得極坐標中Laplace變換域內飽和彈性土的無量綱控制方程為

    引入如下位移勢函數(shù):

    則方程(3)可化為

    方程(5)可進一步解耦為

    由于隧道承受的徑向集中壓力f(θ,t)在極軸上,因此可設位移勢函數(shù)為

    式中,Kn(·)為n階第二類變形Bessel函數(shù);Bn,Dn,En為待定常數(shù);

    利用方程(2)和(4)可求得Laplace變換域中飽和土的位移和應力以及孔隙水的相對位移和壓力:

    2 襯砌的控制方程及求解

    2.1平面應變模型

    將彈性襯砌視為平面應變問題,設襯砌的徑向和環(huán)向位移分別為uLr和uLθ,且uL= (uLr,uLθ,0),則襯砌的動力方程為

    式中,λL和μL為襯砌的Lame常數(shù),ρL為襯砌密度.

    襯砌的應力-應變本構關系為

    式中,σLij(i,j=r,θ),εLij(i,j=r,θ)分別為襯砌的應力和應變.

    考慮初始條件

    利用Laplace變換,并引入如下無量綱參數(shù):

    可得襯砌在Laplace變換域中的無量綱控制方程為

    類似地,引入襯砌位移勢函數(shù)

    由方程(13)可得

    假設彈性襯砌的位移勢函數(shù)為

    式中,In(·)為n階第一類變形Bessel函數(shù);Fn,Gn,Mn,Nn為待定常數(shù).

    因此,由襯砌的應力-應變關系以及位移勢函數(shù)可求得彈性襯砌在Laplace變換域內的位移和應力分別為

    2.2Flügge殼體模型

    若將襯砌等效為薄壁圓柱殼體,在無扭矩條件下,根據(jù)Flügge殼體理論[13],襯砌的運動方程可表示為

    并對方程(19)進行Laplace變換,可得

    3 邊界條件

    隧道內部承受的徑向突加集中常荷載p(t)=p0H(t)關于θ的Fourier展式及其Laplace變換分別為

    3.1飽和土-平面應變模型襯砌系統(tǒng)

    設飽和土與彈性襯砌完全接觸,則襯砌與土體交界面處的位移和應力均連續(xù).由于地下隧道通常需保持干燥無水的運營環(huán)境,因此不考慮襯砌的透水性,即襯砌與土體接觸面處不滲透,由此可得Laplace變換域中土體和襯砌接觸面(η=η2)處的連續(xù)性條件為

    由式(7)~(11),(15)~(18)和(21)~(23)可得,待定常數(shù)Bn,Dn,En,Fn,Gn,Mn,Nn(n=0, 1,…)滿足線性方程組

    3.2飽和土-Flügge殼體模型襯砌系統(tǒng)

    若將襯砌結構視為圓柱殼體,當襯砌厚度h遠小于隧道中曲面半徑R時,可以將襯砌的中曲面(r=R)視為襯砌和土體的接觸面,即忽略h/2[13].在飽和土與彈性襯砌完全接觸且接觸面處不滲透的條件下,Laplace變換域中襯砌與土體接觸面(η=ηm=(1+η2)/2)上的連續(xù)性條件和應力協(xié)調條件為

    3.3無襯砌隧道

    若不考慮隧道襯砌的影響,即隧道無襯砌,有R1=R2,則Laplace變換域中飽和土動力響應的邊界條件為

    類似地,可得到確定待定常數(shù)Bn,Dn和En的線性代數(shù)方程組,從而確定待定常數(shù)Bn,Dn和En.

    為求得飽和土-襯砌系統(tǒng)在時間域中的動力響應,需對Laplace變換域中的解進行Laplace逆變換,然而通常情況下很難求得相應解Laplace逆變換的解析表達式.因此,本工作采用Laplace逆變換的Crump數(shù)值算法[20],得到飽和土-襯砌耦合系統(tǒng)在時間域中的動力響應結果.

    4 數(shù)值結果和分析

    4.1對比驗證

    為了驗證本工作中模型和Crump數(shù)值逆變換計算結果的正確性和可靠性,對均布突加荷載作用下無襯砌隧道的軸對稱動力響應[11]進行對比分析.令n=0,h?=0,ρL?=1, λS?=λL?和μS?=μL?,并按文獻[11]取

    圖2 本工作數(shù)值結果與文獻[11]結果的比較Fig.2 Result comparisons between numerical results presented in this paper and those in Ref.[11]

    圖2給出了無襯砌隧道邊界(η=1)處土體無量綱徑向位移USr、無量綱環(huán)向應力ΣSθ和無量綱孔隙水壓力P在均布突加荷載作用下隨無量綱時間τ的響應.由圖可見,當τ<5時,采用本工作中的方法得到的數(shù)值結果與文獻[11]中的結果吻合良好;當τ>5時,二者略微有所區(qū)別,原因是本工作采用的Laplace數(shù)值逆變換方法與文獻[11]中的方法有所不同,且Crump反演方法的結果具有更明顯的波動性質,計算結果更加準確[21].

    4.2參數(shù)分析

    在通常情況下,可采用3種模型研究隧道附近飽和土的力學性能,即無襯砌飽和土隧道模型、將襯砌視為二維平面問題的飽和土-襯砌耦合系統(tǒng)模型和將襯砌視作無扭矩薄壁圓柱殼體的飽和土-殼體襯砌耦合系統(tǒng)模型.

    采用表1給出的參數(shù),分析飽和土-隧道襯砌系統(tǒng)在突加集中荷載作用下的動力響應.由于采用Fourier級數(shù)進行分析,因此所得數(shù)值結果依賴于Fourier級數(shù)的項數(shù)選取.通過數(shù)值計算發(fā)現(xiàn),當選取Fourier級數(shù)項數(shù)n≥15時,響應增加的幅值小于10-3,故本工作取n=15進行相關的數(shù)值計算.

    表1 計算參數(shù)Table 1 Parameters for calculation

    圖3給出了3種隧道模型中θ=0°處USr,P以及θ=90°處ΣSθ隨τ的動力響應.可見,飽和土-襯砌耦合系統(tǒng)模型與無襯砌飽和土隧道模型的徑向位移、環(huán)向應力和孔隙水壓力隨時間的響應性態(tài)基本一致,而飽和土-殼體襯砌耦合系統(tǒng)模型的響應小于飽和土-平面應變襯砌耦合系統(tǒng)模型的響應,并且徑向位移響應差值較大,襯砌可明顯降低土體的徑向位移和環(huán)向應力,但對孔隙水影響很小.另外,當τ≈3時,徑向位移和環(huán)向應力達到峰值,孔隙水壓力最小;而在突加荷載作用的瞬間,孔隙水壓力最大.而后,隨著時間的延長,由于孔隙水和土體顆粒相互作用的阻尼性質,土體位移和應力逐漸趨為常數(shù),孔隙水壓力逐漸消散,土骨架承擔荷載.

    圖3 集中突加常荷載作用下飽和土的動力響應Fig.3 Dynamic responses of saturated soil under the constant step load

    圖4給出了飽和土-平面應變襯砌耦合系統(tǒng)在τ值不同時土體的USr,ΣSθ和P沿η= r/R1的分布.可見,土體響應具有振動性質.當τ較小時,較遠位置η處的響應為0,表明突加集中荷載引起的響應尚未傳播至該處,并且離隧道較遠(η≥5)處的土體動力響應遠小于隧道附近土體的動力響應.

    圖4 突加集中常荷載作用下飽和土的動力響應沿徑向的分布Fig.4 Dynamic responses of saturated soil under the constant step load along the radial position

    下面以飽和土-平面應變襯砌耦合系統(tǒng)為例進行相關的參數(shù)分析.圖5給出了當I= μL/μS取不同值(其他參數(shù)按表1選取)時,θ=0°處USr,P以及θ=90°處ΣSθ隨τ的動力響應.可見,彈性襯砌剛度對土體的位移和應力響應有顯著的影響,但對孔隙水壓力影響較小.隨著參數(shù)I的增大,土體位移、應力以及孔隙水壓力幅值減小,振動減弱,即適當增加襯砌剛度可以有效減小土體的動力響應.當參數(shù)I較大時,飽和土的徑向位移、環(huán)向應力和孔隙水壓力幾乎為0,此時外荷載基本由襯砌承擔.

    圖6給出了h?取不同值(其他參數(shù)按表1選取)時,θ=0°處USr,P以及θ=90°處ΣSθ隨τ的動力響應.可見,類似于襯砌剛度對土體的影響,h?對土體位移和應力響應影響顯著,但孔隙水壓力變化較小.隨著h?的增大,土體的位移和應力幅值減小,襯砌承擔的外荷載越來越大,土體變形逐漸減小.

    圖7給出了M?取不同值(其他參數(shù)仍按表1選取)時,θ=0°處USr,P以及θ=90°處ΣSθ隨τ的動力響應.可見,隨著M?的增大,土體位移、應力和孔隙水壓力振動幅度和周期增加,但M?對土體位移幅值的影響相對較小.

    圖5 I取不同值時飽和土-襯砌系統(tǒng)的動力響應Fig.5 Dynamic responses of saturated soil-lining system with different I

    圖6 襯砌厚度h?不同時飽和土-襯砌系統(tǒng)的動力響應Fig.6 Dynamic responses of saturated soil-lining system with different thickness of the lining h?

    圖7 M?不同時飽和土-襯砌系統(tǒng)的動力響應Fig.7 Dynamic responses of saturated soil-lining system with different M?

    5 結論

    本工作基于飽和土的Biot理論和彈性理論,采用Laplace變換及其數(shù)值逆變換研究了飽和土-隧道襯砌耦合系統(tǒng)在突加集中荷載作用下的動力響應,得到了Laplace變換域中的解析解,重點分析了土體和襯砌力學和幾何參數(shù)對飽和土動力響應的影響.研究結果表明:①襯砌可明顯降低土體的徑向位移和環(huán)向應力,而且飽和土-殼體襯砌耦合系統(tǒng)模型的響應小于飽和土-平面應變襯砌耦合系統(tǒng)模型的響應;②突加集中荷載作用下土體的動力響應具有明顯的波傳播特征,且離隧道較遠處的土體動力響應遠小于隧道附近土體的動力響應;③襯砌剛度和厚度對土體位移和應力響應影響顯著,但對孔隙水壓力影響較小,適當增加襯砌剛度和厚度可以有效減小土體的動力響應;④孔隙水可壓縮性對土體位移的幅值影響不大,但對應力幅值的影響較為顯著,且隨著孔隙水壓溶性的減小,土體位移、應力和孔隙水壓力振動幅度和周期增加.

    [1]BIOT M A.Theory of propagation of elastic waves in a fluid-saturated porous solid.Ⅰ.Lowfrequency range[J].Journal of Acoustical Society of American,1956,28(2):168-178.

    [2]BIOT M A.Mechanics of deformation and acoustic propagation in porous media[J].Journal of Applied Physics,1962,33(4):1482-1498.

    [3]BOwEN R M.Compressible porous media models by use of the theory of mixtures[J].International Journal of Engineering Science,1982,20(6):697-735.

    [4]DEBOER R.Theory of porous media:highlights in historical development and current state[M]. Berlin:Springer-Verlag,2000.

    [5]XIE K H,LIU G B,SHI Z Y.Dynamic response of partially sealed circular tunnel in viscoelastic saturated soil[J].Soil Dynamic and Earthquake Engineering,2004,24(12):1003-1011.

    [6]劉干斌,謝康和,施祖元.黏彈性飽和多孔介質中圓柱孔洞的頻域響應[J].力學學報,2004,36(5):557-563.

    [7]LIU G B,XIE K H,LIU X H.Dynamic response of a partially sealed tunnel in porous rock under inner water pressure[J].Tunnelling and Underground Space Technology,2010,25(4):407-414.

    [8]楊驍,聞敏杰.飽和分數(shù)導數(shù)型粘彈性土-深埋圓形隧洞襯砌系統(tǒng)的動力特性[J].工程力學,2012, 29(12):248-255.

    [9]楊驍,聞敏杰.深埋圓形隧洞飽和土-襯砌簡諧振動的解析解[J].上海大學學報(自然科學版),2012, 18(5):525-530.

    [10]高華喜,聞敏杰.內水壓作用下粘彈性飽和土-隧洞襯砌相互作用[J].工程力學,2013,30(3):289-296.

    [11]SENjUNTIcHAI T,RAjApAKSE R K N D.Transient response of a circular cavity in a poroelastic medium[J].International Journal for Numerical and Analytical Method in Geomechanics,1993, 17(6):357-383.

    [12]楊峻,宮全美,吳世明,等.飽和土體中圓柱形孔洞的動力分析[J].上海力學,1996,17(1):37-45.

    [13]ZAKOUT U,AKKAS N.Transient response of a cylindrical cavity with and without a bonded shell in an infinite elastic medium[J].International Journal of Engineering Science,1997,35(12):1203-1220.

    [14]丁伯陽,宋新初,袁金華.飽和土隧道內集中荷載作用下振動位移反應的Green函數(shù)解答[J].工程力學,2009,26(6):153-157.

    [15]劉干斌,謝康和,施祖元.粘彈性土體中深埋圓形隧道的應力和位移分析[J].工程力學,2004,21(5):132-138.

    [16]LU J F,JENG D S,LEE T L.Dynamic response of a piecewise circular tunnel embedded in a poroelastic medium[J].Soil Dynamics and Earthquake Engineering,2007,27(9):875-891.

    [17]GAO M,WANG Y,GAO G Y,et al.An analytical solution for the transient response of a cylindrical lined cavity in a poroelastic medium[J].Soil Dynamics and Earthquake Engineering, 2013,46(3):30-40.

    [18]蔡袁強,陳成振,孫宏磊.爆炸荷載作用下飽和土中隧道的瞬態(tài)動力響應[J].巖土工程學報,2011, 33(3):361-367.

    [19]劉干斌,謝康和,施祖元.黏彈性飽和土體中圓形隧洞動力相互作用[J].浙江大學學報(工學版), 2005,39(10):1576-1581.

    [20]CRUMp K S.Numerical inversion of Laplace transforms using a Fourier series approximation[J]. Journal of the ACM,1976,23(1):89-96.

    [21]HONIG G,HIRDES U.A method for the numerical inversion of Laplace transforms[J].Journal of Computational and Applied Mathematics,1984,10(1):113-132.

    本文彩色版可登陸本刊網(wǎng)站查詢:http://www.journal.shu.edu.cn

    Non-axisymmetrical dynamic response of saturated soil and lining system with deeply embedded circular tunnel under step concentrated load

    LU Jianjun,YANG Xiao
    (Department of Civil Engineering,Shanghai University,Shanghai 200072,China)

    Considering the interaction between soil and lining in a deeply embedded circular tunnel,the non-axisymmetrical dynamic response of a saturated soil-lining coupled system under step concentrated load was investigated.Based on the Biot’s theory and elasticity,using the Laplace transform and Fourier series,analytical expressions of displacements,stresses and pore water pressure of the saturated soil-lining system subject to step concentrated load were obtained in the Laplace transform domain under the boundary conditions of the lining and the continuity conditions on the interface between the saturated soil and lining.Numerical solutions of the dynamic responses of the saturated soil and lining system were obtained with a Crump method of the inverse Laplace transform.The influences of mechanical and geometric parameters of the soil and lining on the dynamic response of the system were analyzed.It was shown that dynamic response of the soil in adistance from the tunnel center more than 5 times of the radius was much less than that of the soil in the vicinity of tunnel.Influences of stiffness and thickness on displacements and stresses of the soil were significant,while the influences on pore water pressure were trivial. Furthermore,compressibility of pore water had stronger influence on the amplitude of the stress than that of the displacement of the soil.

    saturated soil;elastic lining;concentrated load;dynamic response;Laplace transform

    O 327;TU 435

    A

    1007-2861(2016)05-0665-15

    10.3969/j.issn.1007-2861.2015.01.003

    2015-02-11

    國家自然科學基金資助項目(10872124)

    楊驍(1965—),男,教授,博士生導師,博士,研究方向為非線性固體力學、多孔介質理論等. E-mail:xyang@shu.edu.cn

    猜你喜歡
    圓形徑向孔隙
    淺探徑向連接體的圓周運動
    RN上一類Kirchhoff型方程徑向對稱正解的存在性
    基于PID+前饋的3MN徑向鍛造機控制系統(tǒng)的研究
    重型機械(2020年3期)2020-08-24 08:31:40
    二氧化碳在高嶺石孔隙中吸附的分子模擬
    一類無窮下級整函數(shù)的Julia集的徑向分布
    Preparation of bimodal grain size 7075 aviation aluminum alloys and the ir corrosion properties
    為什么窨井蓋大多都是圓形的
    肥皂泡為什么是圓形?
    圓形題
    圓形變身喵星人
    国产亚洲精品久久久久5区| 香蕉国产在线看| 国产成人欧美在线观看| 成人一区二区视频在线观看| 欧美精品亚洲一区二区| 国产精品香港三级国产av潘金莲| 亚洲aⅴ乱码一区二区在线播放 | 婷婷精品国产亚洲av在线| 国产一区在线观看成人免费| 女生性感内裤真人,穿戴方法视频| www.www免费av| tocl精华| 亚洲av电影不卡..在线观看| 久久久久国产精品人妻aⅴ院| 精品电影一区二区在线| 亚洲熟妇熟女久久| 亚洲欧美日韩高清专用| 免费观看精品视频网站| 国产高清视频在线播放一区| 亚洲精品色激情综合| 制服丝袜大香蕉在线| 国产成人av教育| 啦啦啦观看免费观看视频高清| 日本成人三级电影网站| 男女做爰动态图高潮gif福利片| 精品免费久久久久久久清纯| 国产精品久久电影中文字幕| 黄色成人免费大全| 美女免费视频网站| 亚洲成人国产一区在线观看| 一级毛片精品| 亚洲人成电影免费在线| 女警被强在线播放| 岛国视频午夜一区免费看| 波多野结衣高清作品| 亚洲熟妇熟女久久| 亚洲av电影在线进入| 国产免费av片在线观看野外av| 夜夜看夜夜爽夜夜摸| 久久精品国产综合久久久| svipshipincom国产片| 亚洲精品国产一区二区精华液| 怎么达到女性高潮| 真人做人爱边吃奶动态| 好看av亚洲va欧美ⅴa在| 在线观看舔阴道视频| 一级毛片女人18水好多| 国产午夜精品久久久久久| 国产精品免费一区二区三区在线| 欧美黑人精品巨大| avwww免费| 91麻豆精品激情在线观看国产| 久久香蕉国产精品| 欧美日韩乱码在线| 男插女下体视频免费在线播放| 两性夫妻黄色片| 亚洲人成伊人成综合网2020| 成人亚洲精品av一区二区| 精品国内亚洲2022精品成人| 久久久久免费精品人妻一区二区| 一二三四在线观看免费中文在| 男人舔女人的私密视频| 久久亚洲精品不卡| 一个人免费在线观看的高清视频| 在线观看午夜福利视频| 欧美日韩一级在线毛片| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩精品网址| 久久天堂一区二区三区四区| 91国产中文字幕| 日本一本二区三区精品| 国产一区二区三区视频了| 免费av毛片视频| 久久久久久久精品吃奶| 亚洲国产精品久久男人天堂| 熟女电影av网| 中文字幕人妻丝袜一区二区| 麻豆成人午夜福利视频| 一区二区三区国产精品乱码| 午夜福利欧美成人| 一本一本综合久久| 看片在线看免费视频| 欧美+亚洲+日韩+国产| 亚洲无线在线观看| 久久久久国产一级毛片高清牌| 日本熟妇午夜| 国产av一区二区精品久久| 夜夜爽天天搞| 亚洲av成人精品一区久久| 国产99白浆流出| 色噜噜av男人的天堂激情| 51午夜福利影视在线观看| 亚洲成a人片在线一区二区| 欧美高清成人免费视频www| www.自偷自拍.com| 欧美一区二区国产精品久久精品 | 国内精品久久久久精免费| 久久久久久免费高清国产稀缺| 久久中文看片网| 欧美色视频一区免费| 91老司机精品| 99国产极品粉嫩在线观看| 亚洲性夜色夜夜综合| 首页视频小说图片口味搜索| 少妇熟女aⅴ在线视频| 精品久久久久久,| 亚洲av熟女| 日韩欧美在线二视频| 国产精品亚洲一级av第二区| 国产精品九九99| 一进一出抽搐动态| www日本在线高清视频| 一a级毛片在线观看| 亚洲五月天丁香| 天堂影院成人在线观看| 大型黄色视频在线免费观看| 国产一区二区三区在线臀色熟女| 精品欧美一区二区三区在线| 在线观看www视频免费| 国产av不卡久久| 无人区码免费观看不卡| 婷婷亚洲欧美| 免费电影在线观看免费观看| 非洲黑人性xxxx精品又粗又长| 欧美绝顶高潮抽搐喷水| 国产精品免费一区二区三区在线| 别揉我奶头~嗯~啊~动态视频| 99国产精品99久久久久| 色噜噜av男人的天堂激情| 久久久久国内视频| 男男h啪啪无遮挡| 久久人人精品亚洲av| 两个人视频免费观看高清| 欧洲精品卡2卡3卡4卡5卡区| √禁漫天堂资源中文www| 在线十欧美十亚洲十日本专区| 久久亚洲精品不卡| 三级毛片av免费| 不卡av一区二区三区| 亚洲片人在线观看| 亚洲人成电影免费在线| 国语自产精品视频在线第100页| 在线永久观看黄色视频| 一级作爱视频免费观看| 亚洲最大成人中文| 国产av不卡久久| 成年免费大片在线观看| 中文亚洲av片在线观看爽| 久久久水蜜桃国产精品网| 国内精品久久久久久久电影| 小说图片视频综合网站| 久久久久久大精品| svipshipincom国产片| 12—13女人毛片做爰片一| 免费看美女性在线毛片视频| 正在播放国产对白刺激| 午夜免费成人在线视频| 母亲3免费完整高清在线观看| 精品久久久久久久久久免费视频| 村上凉子中文字幕在线| 午夜福利视频1000在线观看| 好看av亚洲va欧美ⅴa在| 国产精品久久久人人做人人爽| 亚洲色图 男人天堂 中文字幕| 亚洲欧美日韩高清在线视频| 日韩国内少妇激情av| 成人av一区二区三区在线看| 欧美又色又爽又黄视频| 国产一区二区在线观看日韩 | 美女扒开内裤让男人捅视频| 特级一级黄色大片| 久久久久久人人人人人| 国产成年人精品一区二区| 亚洲五月天丁香| www.精华液| 91九色精品人成在线观看| 97碰自拍视频| 亚洲精品国产精品久久久不卡| 国产成人精品久久二区二区91| 亚洲免费av在线视频| 啦啦啦韩国在线观看视频| 18禁观看日本| 又大又爽又粗| 久久精品国产综合久久久| 国产蜜桃级精品一区二区三区| 日韩有码中文字幕| 成人高潮视频无遮挡免费网站| 日韩欧美三级三区| 亚洲精品国产一区二区精华液| 人妻夜夜爽99麻豆av| 欧美日韩乱码在线| 亚洲av中文字字幕乱码综合| 国产真人三级小视频在线观看| 国产亚洲av嫩草精品影院| 国产亚洲精品一区二区www| 婷婷六月久久综合丁香| 亚洲人成电影免费在线| av免费在线观看网站| 老汉色av国产亚洲站长工具| www.精华液| 日本免费一区二区三区高清不卡| а√天堂www在线а√下载| 久久热在线av| 亚洲国产看品久久| 三级毛片av免费| 国产黄片美女视频| 亚洲性夜色夜夜综合| 精品久久久久久久久久久久久| 午夜视频精品福利| 欧美成人一区二区免费高清观看 | 午夜福利在线在线| 成人手机av| 女人爽到高潮嗷嗷叫在线视频| 久久婷婷人人爽人人干人人爱| 免费在线观看黄色视频的| 白带黄色成豆腐渣| 国产精品精品国产色婷婷| 伦理电影免费视频| 97碰自拍视频| 国产熟女午夜一区二区三区| av免费在线观看网站| 露出奶头的视频| 国产精品永久免费网站| av天堂在线播放| 中国美女看黄片| 色综合婷婷激情| 国产免费av片在线观看野外av| 免费在线观看成人毛片| 亚洲一区二区三区色噜噜| 一个人免费在线观看电影 | 国产99白浆流出| 亚洲精品中文字幕一二三四区| 夜夜躁狠狠躁天天躁| 97人妻精品一区二区三区麻豆| 两性夫妻黄色片| 天天一区二区日本电影三级| 欧美性猛交╳xxx乱大交人| svipshipincom国产片| 国产av在哪里看| 国产精华一区二区三区| 欧美性长视频在线观看| av在线天堂中文字幕| 观看免费一级毛片| 国产av一区二区精品久久| e午夜精品久久久久久久| 全区人妻精品视频| 亚洲黑人精品在线| 亚洲av片天天在线观看| 亚洲男人天堂网一区| 国产单亲对白刺激| 午夜福利高清视频| 青草久久国产| 久热爱精品视频在线9| 亚洲狠狠婷婷综合久久图片| 国产99久久九九免费精品| 久久婷婷人人爽人人干人人爱| 两人在一起打扑克的视频| 999精品在线视频| 啪啪无遮挡十八禁网站| 又黄又粗又硬又大视频| 听说在线观看完整版免费高清| xxx96com| 国产免费av片在线观看野外av| 久久久久久亚洲精品国产蜜桃av| 久久久久亚洲av毛片大全| 嫁个100分男人电影在线观看| 一本大道久久a久久精品| 天堂√8在线中文| 在线播放国产精品三级| 亚洲片人在线观看| 欧美丝袜亚洲另类 | 日韩三级视频一区二区三区| 亚洲在线自拍视频| 亚洲,欧美精品.| 国产v大片淫在线免费观看| 久久久久久人人人人人| 麻豆国产av国片精品| 国产精品99久久99久久久不卡| 中文字幕熟女人妻在线| 国产av麻豆久久久久久久| 国产精品,欧美在线| 老司机午夜福利在线观看视频| 国产三级中文精品| 亚洲人与动物交配视频| 日韩欧美 国产精品| 欧美黑人巨大hd| 又黄又粗又硬又大视频| 亚洲九九香蕉| 人人妻人人看人人澡| 五月玫瑰六月丁香| 91老司机精品| 深夜精品福利| 成人18禁高潮啪啪吃奶动态图| 亚洲成人中文字幕在线播放| 免费看日本二区| 国产精品影院久久| 国内精品久久久久久久电影| 久久欧美精品欧美久久欧美| 在线观看午夜福利视频| 国产熟女xx| 日韩三级视频一区二区三区| 欧美一级a爱片免费观看看 | 色综合婷婷激情| 天堂√8在线中文| 观看免费一级毛片| 久久精品国产清高在天天线| 国产精品久久电影中文字幕| 中文字幕人成人乱码亚洲影| 精品欧美一区二区三区在线| 日韩成人在线观看一区二区三区| 欧美最黄视频在线播放免费| 一进一出好大好爽视频| svipshipincom国产片| 啦啦啦免费观看视频1| 桃色一区二区三区在线观看| 2021天堂中文幕一二区在线观| 狂野欧美白嫩少妇大欣赏| 国产亚洲精品久久久久久毛片| 国语自产精品视频在线第100页| 天天躁狠狠躁夜夜躁狠狠躁| 搡老妇女老女人老熟妇| 欧美一级毛片孕妇| 亚洲一区中文字幕在线| av国产免费在线观看| 亚洲国产精品成人综合色| 精品第一国产精品| 国产野战对白在线观看| 中国美女看黄片| 成人精品一区二区免费| av国产免费在线观看| a在线观看视频网站| 亚洲av中文字字幕乱码综合| 日韩精品免费视频一区二区三区| 成人手机av| 一本综合久久免费| 国产一区二区在线观看日韩 | 国产伦在线观看视频一区| 18禁国产床啪视频网站| 可以在线观看的亚洲视频| 成年版毛片免费区| 亚洲av成人av| 夜夜躁狠狠躁天天躁| 男女床上黄色一级片免费看| 国产aⅴ精品一区二区三区波| 日本 av在线| 欧美丝袜亚洲另类 | 老汉色av国产亚洲站长工具| 亚洲欧美日韩高清在线视频| 午夜免费观看网址| 亚洲人成网站高清观看| 亚洲精品色激情综合| 成人18禁在线播放| 宅男免费午夜| 精品第一国产精品| 久久精品影院6| 黄色a级毛片大全视频| 一进一出好大好爽视频| 黑人欧美特级aaaaaa片| 亚洲精品美女久久av网站| 国产在线精品亚洲第一网站| 亚洲熟妇熟女久久| 中文字幕av在线有码专区| 久久香蕉国产精品| 亚洲午夜理论影院| 久久久精品欧美日韩精品| 波多野结衣巨乳人妻| 黄色视频不卡| 亚洲午夜理论影院| 90打野战视频偷拍视频| 久久精品国产综合久久久| 黄色视频,在线免费观看| 中文资源天堂在线| 免费观看人在逋| 老司机深夜福利视频在线观看| 亚洲狠狠婷婷综合久久图片| 两性午夜刺激爽爽歪歪视频在线观看 | 国产高清有码在线观看视频 | 91成年电影在线观看| 脱女人内裤的视频| 久久精品综合一区二区三区| 亚洲 欧美一区二区三区| 国产精品九九99| 成人精品一区二区免费| 丝袜美腿诱惑在线| 色老头精品视频在线观看| 国产精品久久久久久人妻精品电影| 老司机在亚洲福利影院| 午夜免费激情av| 亚洲人成网站在线播放欧美日韩| 国产精华一区二区三区| 久久久久久九九精品二区国产 | 在线观看一区二区三区| 日韩欧美国产在线观看| 天堂√8在线中文| 免费一级毛片在线播放高清视频| 黄色视频不卡| 好看av亚洲va欧美ⅴa在| 午夜亚洲福利在线播放| 超碰成人久久| 99久久无色码亚洲精品果冻| 精品少妇一区二区三区视频日本电影| tocl精华| 久久久久亚洲av毛片大全| 久久这里只有精品19| 亚洲 国产 在线| 日本撒尿小便嘘嘘汇集6| 亚洲欧美日韩无卡精品| 熟女少妇亚洲综合色aaa.| 日本 av在线| 亚洲天堂国产精品一区在线| 国产黄a三级三级三级人| 男男h啪啪无遮挡| 丝袜美腿诱惑在线| 亚洲国产欧美网| 一个人免费在线观看电影 | 两个人免费观看高清视频| 中文字幕最新亚洲高清| 床上黄色一级片| 看片在线看免费视频| 国产精品精品国产色婷婷| 不卡av一区二区三区| 母亲3免费完整高清在线观看| 一边摸一边做爽爽视频免费| 在线观看免费午夜福利视频| 无遮挡黄片免费观看| 老司机午夜福利在线观看视频| 亚洲精品久久成人aⅴ小说| 中文字幕最新亚洲高清| 黄色视频不卡| 在线观看舔阴道视频| 日日干狠狠操夜夜爽| 成年版毛片免费区| 午夜精品在线福利| 国产三级中文精品| 成人三级做爰电影| 丰满人妻一区二区三区视频av | av福利片在线| 天天躁狠狠躁夜夜躁狠狠躁| 在线播放国产精品三级| 亚洲熟妇熟女久久| 欧美乱色亚洲激情| 老鸭窝网址在线观看| 在线a可以看的网站| 好男人在线观看高清免费视频| 男人舔奶头视频| 亚洲第一欧美日韩一区二区三区| 久久精品国产清高在天天线| 国产激情偷乱视频一区二区| 性色av乱码一区二区三区2| 激情在线观看视频在线高清| 成年人黄色毛片网站| 嫩草影视91久久| 首页视频小说图片口味搜索| 日韩欧美一区二区三区在线观看| 亚洲av成人精品一区久久| 欧美一区二区国产精品久久精品 | videosex国产| 99国产精品一区二区三区| 国产av麻豆久久久久久久| 在线观看一区二区三区| 久久久精品大字幕| 一个人免费在线观看的高清视频| 欧美日韩国产亚洲二区| 老鸭窝网址在线观看| 欧美日韩黄片免| 美女 人体艺术 gogo| 男人舔奶头视频| 中出人妻视频一区二区| 欧美日本亚洲视频在线播放| 免费在线观看成人毛片| 亚洲午夜理论影院| 午夜两性在线视频| 一级片免费观看大全| 国内精品一区二区在线观看| 88av欧美| 亚洲avbb在线观看| 国产野战对白在线观看| 欧美日韩福利视频一区二区| 日韩欧美 国产精品| 欧美人与性动交α欧美精品济南到| 久久香蕉精品热| 搡老妇女老女人老熟妇| 中文资源天堂在线| av在线天堂中文字幕| 亚洲专区字幕在线| 国内毛片毛片毛片毛片毛片| 精品国内亚洲2022精品成人| 久久天堂一区二区三区四区| 国产片内射在线| 免费av毛片视频| 最新美女视频免费是黄的| 中文字幕人成人乱码亚洲影| 人妻久久中文字幕网| 日韩欧美精品v在线| 欧美黄色淫秽网站| 国产亚洲精品第一综合不卡| 欧美乱妇无乱码| 中文亚洲av片在线观看爽| 两个人视频免费观看高清| 特级一级黄色大片| 一夜夜www| 亚洲精品av麻豆狂野| АⅤ资源中文在线天堂| 观看免费一级毛片| 久久久久久久久久黄片| 日日夜夜操网爽| 亚洲精品美女久久av网站| 国产av麻豆久久久久久久| 一个人免费在线观看的高清视频| 麻豆一二三区av精品| 九色成人免费人妻av| 亚洲欧美日韩高清在线视频| 黄片小视频在线播放| √禁漫天堂资源中文www| 1024香蕉在线观看| 久久久久久久久免费视频了| 美女黄网站色视频| 午夜久久久久精精品| 日韩成人在线观看一区二区三区| 午夜激情福利司机影院| 午夜影院日韩av| 这个男人来自地球电影免费观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲一卡2卡3卡4卡5卡精品中文| 嫩草影视91久久| 一进一出抽搐动态| 波多野结衣高清无吗| 欧美一级毛片孕妇| 国产欧美日韩一区二区三| 中文字幕最新亚洲高清| 国产黄a三级三级三级人| 搡老熟女国产l中国老女人| 成人午夜高清在线视频| 国产精品1区2区在线观看.| 久久天躁狠狠躁夜夜2o2o| 在线观看免费视频日本深夜| 色精品久久人妻99蜜桃| 人妻夜夜爽99麻豆av| 三级毛片av免费| 午夜影院日韩av| √禁漫天堂资源中文www| 88av欧美| 成人精品一区二区免费| 国产黄色小视频在线观看| 在线a可以看的网站| 国产真实乱freesex| 国产三级中文精品| 精品午夜福利视频在线观看一区| 色综合亚洲欧美另类图片| 成人一区二区视频在线观看| 亚洲av电影在线进入| 亚洲自偷自拍图片 自拍| 中亚洲国语对白在线视频| 日本黄大片高清| 国产成人一区二区三区免费视频网站| 一级a爱片免费观看的视频| 国产人伦9x9x在线观看| 亚洲熟妇中文字幕五十中出| 韩国av一区二区三区四区| 99国产精品一区二区蜜桃av| 国产精品免费视频内射| 制服诱惑二区| 伊人久久大香线蕉亚洲五| 欧美3d第一页| 亚洲国产精品999在线| 婷婷丁香在线五月| 国产精品av视频在线免费观看| 亚洲国产精品久久男人天堂| 草草在线视频免费看| 成人18禁在线播放| 国产精品,欧美在线| 久久久精品大字幕| 99热只有精品国产| 国产人伦9x9x在线观看| 在线免费观看的www视频| 免费一级毛片在线播放高清视频| 精品国产乱码久久久久久男人| 久热爱精品视频在线9| 高清毛片免费观看视频网站| 一区二区三区高清视频在线| 亚洲免费av在线视频| 欧美黑人精品巨大| 脱女人内裤的视频| 久久人妻福利社区极品人妻图片| 一二三四在线观看免费中文在| 午夜免费激情av| 我的老师免费观看完整版| 欧美日韩一级在线毛片| 色哟哟哟哟哟哟| 久久久久九九精品影院| 久久久久久人人人人人| 成人国语在线视频| 亚洲,欧美精品.| 欧美高清成人免费视频www| 99国产精品一区二区蜜桃av| 麻豆成人av在线观看| 91国产中文字幕| 亚洲专区中文字幕在线| 亚洲一区高清亚洲精品| 99久久精品热视频| 中文字幕av在线有码专区| 欧美黑人巨大hd| 18禁观看日本| 成熟少妇高潮喷水视频| 久久天堂一区二区三区四区| 国产精品电影一区二区三区| 日本a在线网址| 婷婷亚洲欧美| 精品欧美国产一区二区三| 成人国产一区最新在线观看| 精品久久久久久成人av| 久久久久久久久久黄片| 久久久久免费精品人妻一区二区| 1024视频免费在线观看| 国产午夜福利久久久久久|