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

    局部動(dòng)荷載作用下土拱效應(yīng)的離散元分析

    2022-03-18 02:27:52李格燁沈盼盼張興亞羅敏敏
    長江科學(xué)院院報(bào) 2022年3期
    關(guān)鍵詞:模型試驗(yàn)填料幅值

    李格燁,徐 超,沈盼盼,張興亞,羅敏敏

    (1.同濟(jì)大學(xué) 土木工程學(xué)院,上海 200092; 2.中國長江三峽集團(tuán)有限公司 上??睖y(cè)設(shè)計(jì)研究院有限公司,上海 200434; 3.新城控股集團(tuán)股份有限公司,上海 200092; 4.浙江大學(xué)建筑設(shè)計(jì)研究院有限公司,杭州 310027)

    1 研究背景

    土拱效應(yīng)是土工結(jié)構(gòu)中常見的現(xiàn)象。Terzaghi[1]通過Trapdoor試驗(yàn)證實(shí)了土拱效應(yīng)的存在,并提出了“土拱效應(yīng)”的概念。為了評(píng)價(jià)土拱效應(yīng)的強(qiáng)弱,McNulty[2]定義了土拱率ρ,其表達(dá)式為

    (1)

    式中:σv為活動(dòng)土體底部作用的平均豎向應(yīng)力;γH為活動(dòng)土體的自重應(yīng)力,γ為活動(dòng)土體的重度,H為活動(dòng)土體的高度;q為土體表面作用的附加荷載。

    自重荷載和表面靜荷載作用下土拱效應(yīng)的研究已有較多的成果,涉及土拱的變形特性、應(yīng)力重分布規(guī)律及理論計(jì)算等。然而,在實(shí)際工程中,土拱不可避免地受到動(dòng)荷載的影響。越來越多的學(xué)者對(duì)循環(huán)動(dòng)荷載作用下的土拱效應(yīng)開展研究。在模型試驗(yàn)方面,韓高孝等[3]以某樁-網(wǎng)結(jié)構(gòu)斷面為原型,開展了三維模型試驗(yàn)研究,結(jié)果表明在循環(huán)正弦波荷載作用下,先前由樁所承擔(dān)的部分動(dòng)應(yīng)力會(huì)轉(zhuǎn)移至樁間土,這意味著動(dòng)荷載作用下土拱效應(yīng)會(huì)發(fā)生退化。Xu等[4-5]利用自制的平面應(yīng)變Trapdoor模型試驗(yàn)箱進(jìn)行了一系列試驗(yàn),結(jié)果表明,局部靜動(dòng)荷載均會(huì)弱化土拱效應(yīng),動(dòng)荷載對(duì)土拱效應(yīng)的退化程度大于靜荷載。Bi等[6]采用砂土開展了Trapdoor試驗(yàn),研究了均布循環(huán)荷載下土拱的結(jié)構(gòu)形態(tài)、位移場(chǎng)、豎向應(yīng)力的演變規(guī)律,再現(xiàn)了土拱從形成到完全塌陷的過程。Aqoub等[7]研究了循環(huán)荷載下低填方樁承式路堤的性能,表明土拱的破壞在荷載的前幾個(gè)周期內(nèi)發(fā)生,但是,路堤材料的振動(dòng)密實(shí)以及路基軟土的變形導(dǎo)致土體強(qiáng)度恢復(fù),進(jìn)而致使土拱效應(yīng)部分或全部恢復(fù)。在數(shù)值計(jì)算方面,賴漢江等[8]通過離散元數(shù)值模擬,分析了動(dòng)荷載作用下樁承式路堤中的土拱效應(yīng),結(jié)果表明,循環(huán)荷載作用下,低填方路堤中土拱效應(yīng)先逐漸減弱,最終趨于穩(wěn)定。梁自立等[9]采用 ABAQUS 有限元分析軟件,分析了列車動(dòng)荷載(循環(huán)正弦波荷載)對(duì)樁網(wǎng)結(jié)構(gòu)路基土拱效應(yīng)的影響,結(jié)果表明動(dòng)荷載對(duì)土拱效應(yīng)的影響與路基高度、樁間土性質(zhì)密切相關(guān);當(dāng)路基內(nèi)部形成完整土拱時(shí),動(dòng)荷載對(duì)土拱效應(yīng)基本無影響,但會(huì)引起低矮路基中的非完整土拱發(fā)生失穩(wěn)、退化。Pham等[10]利用有限元方法建立了樁承式加筋路堤三維模型,結(jié)果表明,循環(huán)荷載次數(shù)不同,土拱效應(yīng)退化程度不同;交通荷載的速度越大,即荷載頻率越大,土拱效應(yīng)退化更明顯。

    上述研究成果集中于動(dòng)荷載作用下土拱效應(yīng)的退化甚至是破壞現(xiàn)象的認(rèn)識(shí),在作用機(jī)理、定量評(píng)價(jià)等方面有待進(jìn)一步深入。本文基于Xu等(2019)[4-5]開展的室內(nèi)Trapdoor模型試驗(yàn),采用PFC2D建立了相應(yīng)的離散元模型,研究了內(nèi)摩擦角、荷載幅值和荷載頻率對(duì)土拱效應(yīng)的影響。

    2 離散元模型

    2.1 Trapdoor模型試驗(yàn)

    Xu等(2019)[4-5]采用自制的Trapdoor模型試驗(yàn)箱開展了試驗(yàn),其中填料采用橫斷面為圓形的鋁棒相似土。鋁棒相似土的單位重度為22.4 kN/m3,孔隙率約為17%?;顒?dòng)門的寬度B=128 mm。試驗(yàn)過程中先下移活動(dòng)門至總位移達(dá)到8 mm,形成初始土拱,接著在填料表面施加動(dòng)荷載。其他信息參見Xu等(2019)[4-5]。

    2.2 離散元模型建立

    2.2.1 試樣制備方法

    試樣制備方法及過程決定了初始狀態(tài)下試樣中顆粒的孔隙及受力分布是否合理,因而對(duì)試驗(yàn)的結(jié)果有重要的影響。常用的試樣制備方法有:半徑擴(kuò)大和重力沉積法(Bhandari等[11]、Han等[12])、改進(jìn)的多層壓縮法(Improved Multi-layer Compaction Method,IMCM)(Lai等[13])等。半徑擴(kuò)大和重力沉積法生成的試樣內(nèi)顆粒分布很不均勻,從而導(dǎo)致整個(gè)試樣內(nèi)都會(huì)出現(xiàn)較高的初始地應(yīng)力。更關(guān)鍵的是,試樣中的應(yīng)力并不隨高度呈現(xiàn)梯度變化,這種現(xiàn)象并不符合實(shí)際情況(Rui等[14])。IMCM試樣制備的過程與模型試驗(yàn)過程相似,因此初始應(yīng)力分布與實(shí)際大致相同。本文數(shù)值模擬基于PFC2D5.0軟件,借鑒了IMCM的基本原理,提出了計(jì)算效率更高的多層收斂法,如圖1所示。其主要過程如下:

    圖1 試樣制備過程Fig.1 Flow chart for generating sample

    (1)在指定區(qū)域內(nèi)按照填料的孔隙率和質(zhì)量比生成Ball顆粒層,如圖1(a)所示,其顆粒接觸如圖1(b)所示,并對(duì)其進(jìn)行求解直至收斂,結(jié)果如圖1(c)所示。需要指出的是,在試樣達(dá)到指定高度前,顆粒重力加速度g和表面摩擦系數(shù)μ均設(shè)置為0。

    (2)在已經(jīng)形成的顆粒層上方創(chuàng)建新的墻體形成新的空間,按照步驟(1)生成新的Ball顆粒層,如圖1(d)所示,這時(shí)顆粒接觸如圖1(e)所示。刪除2層顆粒間的墻體并再次計(jì)算至收斂,此時(shí)顆粒間接觸消失,如圖1(f)所示。

    (3)重復(fù)第(2)步,直至達(dá)到試樣的目標(biāo)高度。

    (4)對(duì)所有顆粒賦予重力加速度g及合理的表面摩擦系數(shù)μ,并計(jì)算至收斂。

    2.2.2 Trapdoor 離散元模型

    基于Trapdoor模型試驗(yàn)[4],建立了如圖2所示的Trapdoor離散元模型,其中,模型箱和活動(dòng)門(B1—B7)均采用Wall單元模擬,填料利用2.2.1節(jié)中多層收斂法生成的Ball顆粒模擬。每次生成厚度為20 mm的填料(與模型試驗(yàn)相同),直至試樣高度達(dá)到2B。在活動(dòng)門矩形區(qū)域上方有均勻分布的Ball顆粒來模擬物理模型中的墊層,墊層厚度h=22.4 mm。采用Clump單元模擬加載板,在Clump上施加動(dòng)荷載模擬表面荷載。

    圖2 Trapdoor離散元模型Fig.2 Trapdoor discrete model

    2.3 細(xì)觀參數(shù)標(biāo)定

    填料的細(xì)觀特性通過數(shù)值雙軸試驗(yàn)標(biāo)定,其細(xì)觀力學(xué)參數(shù)列于表1。與物理試驗(yàn)相同,數(shù)值雙軸試驗(yàn)的試樣是將3、4、5 mm三種直徑的顆粒以質(zhì)量比1∶1∶1混合而成,采用50、100、150 kPa三級(jí)圍壓,試驗(yàn)所得偏應(yīng)力及體積應(yīng)變隨軸向應(yīng)變的變化關(guān)系如圖3所示。數(shù)值雙軸試驗(yàn)結(jié)果驗(yàn)證了文中顆粒填料細(xì)觀參數(shù)取值的合理性,求解得出顆粒填料的宏觀內(nèi)摩擦角為24.8°,與實(shí)際雙軸試驗(yàn)得出的結(jié)果25°十分接近。

    表1 填料的細(xì)觀力學(xué)參數(shù)Table 1 Micromechanics properties of fillings

    圖3 數(shù)值及物理雙軸試驗(yàn)中偏應(yīng)力和體積應(yīng)變隨軸向應(yīng)變的曲線Fig.3 Curves of deviatoric stress and volumetric strain versus axial strain obtained from numerical simulation and laboratory biaxial tests

    3 離散元模擬結(jié)果與分析

    選取Xu等(2019)[4]中的正弦波CT4開展數(shù)值模擬,即填料高度H= 2B,加載板的長度L= 2B。數(shù)值試驗(yàn)過程與模型試驗(yàn)過程相同。

    3.1 自重荷載下Trapdoor試驗(yàn)

    在數(shù)值試驗(yàn)中,通過Fish語言編程記錄活動(dòng)門B4每下移0.1 mm時(shí)作用在活動(dòng)門上的應(yīng)力,繪制土拱率ρ隨Trapdoor歸一化位移(Trapdoor位移與Trapdoor寬度的比值)的變化關(guān)系曲線,即地基反應(yīng)曲線(Ground Reaction Curve,GRC)。圖4給出了模型試驗(yàn)和數(shù)值模擬對(duì)應(yīng)的GRC。整體而言,模型試驗(yàn)與數(shù)值模擬中GRC的重合度較高,均分為初始拱、最大拱、應(yīng)力恢復(fù)及極限狀態(tài)4個(gè)階段。最大拱對(duì)應(yīng)歸一化位移為1.25%,模型試驗(yàn)中的最大拱對(duì)應(yīng)的土拱率,即最小土拱率是0.31,數(shù)值模擬中該值為0.32;極限狀態(tài)開始階段對(duì)應(yīng)的歸一化位移均>10%。因此,自重條件下數(shù)值模型較好地還原了物理模型。

    圖4 模型試驗(yàn)和數(shù)值模擬的GRCFig.4 GRCs obtained from model tests and numerical simulation

    3.2 局部動(dòng)荷載作用下Trapdoor試驗(yàn)

    以幅值為8 kPa、頻率為0.5 Hz,加載時(shí)間為200 s的正弦波(CT4)作為動(dòng)荷載,施加在數(shù)值模型的填料表面。表2中,8 mm位移是指活動(dòng)門B4下降8 mm,形成初始土拱,并未施加表面動(dòng)荷載的階段。8 kPa峰值是指在表面荷載8 kPa作用下,最后一個(gè)循環(huán)周期峰值所對(duì)應(yīng)的土拱。8 kPa谷值是指最后一個(gè)循環(huán)周期谷值所對(duì)應(yīng)的土拱。從表2可以看出,在局部動(dòng)荷載加載階段(荷載峰值)作用下,土拱效應(yīng)有一定程度的退化。在卸載階段(荷載谷值),土拱效應(yīng)有一定程度的恢復(fù)。對(duì)比模型試驗(yàn)和數(shù)值模擬的結(jié)果,可以發(fā)現(xiàn)各階段土拱率均很接近。從應(yīng)力角度來看,數(shù)值模型對(duì)物理模型的還原度較好。

    表2 模型試驗(yàn)及數(shù)值模擬對(duì)應(yīng)的土拱率Table 2 Soil arching ratios obtained from numerical simulation and model tests

    圖5(a)顯示的是8 mm階段所對(duì)應(yīng)的模型試驗(yàn)和數(shù)值模擬的填料豎向位移云圖。其填料的變形模式與塔形變形模式(Rui等[14])相同。圖5(b)和圖5(c)分別是8 kPa峰值和8 kPa谷值所對(duì)應(yīng)的填料位移云圖。如圖5所示,模型試驗(yàn)和數(shù)值模擬的變形形態(tài)及變形大小二者均非常接近。從變形角度來看,數(shù)值模型很好地還原了物理模型。

    圖5 模型試驗(yàn)和數(shù)值模擬填料的豎向位移云圖Fig.5 Vertical displacement contours of the backfill of model tests and numerical simulation

    為了定量衡量動(dòng)荷載作用下土拱效應(yīng)的退化程度,定義了土拱退化參數(shù)α,如式(2)所示。

    (2)

    式中:ρload代表動(dòng)荷載作用結(jié)束后的土拱率;ρinitial代表動(dòng)荷載施加前的土拱率。當(dāng)α=0時(shí)表明土拱效應(yīng)在動(dòng)荷載作用下沒有退化;當(dāng)α=100% 時(shí)表明土拱效應(yīng)消失,即ρload=1。在8 kPa動(dòng)荷載作用下,土拱退化參數(shù)α=64%,表明土拱效應(yīng)退化了64%,即動(dòng)荷載作用下土拱的發(fā)揮程度是初始土拱的36%。

    填料顆粒間接觸力是土拱效應(yīng)的細(xì)觀表現(xiàn)。圖6展示的是加載前后接觸力的變化規(guī)律,接觸力線條顏色越深代表接觸力越大。圖6(a)是8 mm位移階段接觸力的分布情況,可以看出,穩(wěn)定木塊上方(B3和B5)的接觸力明顯大于活動(dòng)門(B4),且力鏈沿高度發(fā)生偏轉(zhuǎn),形成明顯的拱形結(jié)構(gòu)。圖6(b)展示的是8 kPa谷值所對(duì)應(yīng)的接觸力分布,穩(wěn)定木塊上方(B3和B5)的接觸力與活動(dòng)門(B4)差別不明顯,拱形結(jié)構(gòu)基本消失,表明土拱效應(yīng)退化明顯。

    圖6 接觸力分布Fig.6 Distribution of contact forces

    4 參數(shù)分析

    4.1 內(nèi)摩擦角

    通過改變表1中顆粒表面摩擦研究試樣內(nèi)摩擦角對(duì)土拱發(fā)揮程度的影響。顆粒表面摩擦系數(shù)μ=0.5、0.8、1.2對(duì)應(yīng)填料的內(nèi)摩擦角分別為24.8°、29.7°、32.6°。數(shù)值模型中填料高度H= 2B,荷載幅值為8 kPa,荷載頻率為0.5 Hz,加載范圍L= 2B。

    圖7展示了無荷載、動(dòng)荷載峰值對(duì)應(yīng)的土拱率與顆粒表面摩擦系數(shù)間的關(guān)系。填料的土拱率均隨試樣內(nèi)摩擦角的增大而減小,表明土拱效應(yīng)增強(qiáng)。但是,動(dòng)荷載作用下土拱率的值均大于無荷載情況,說明在動(dòng)荷載作用下土拱有一定程度的退化。土拱退化參數(shù)α與顆粒表面摩擦系數(shù)的關(guān)系如圖7所示,隨著顆粒表面摩擦系數(shù)的增大,土拱退化參數(shù)α逐漸減小,說明土拱效應(yīng)的退化程度降低,抵御動(dòng)荷載的能力增強(qiáng)。

    圖7 土拱率、土拱退化參數(shù)與顆粒表面摩擦系數(shù)間的關(guān)系Fig.7 Variations of soil arching ratio and degradation parameter of soil arching against surface frictional coefficient of particle

    4.2 荷載幅值

    數(shù)值模型中填料的高度H=2B和H=3B,動(dòng)荷載加載頻率為0.5 Hz,加載范圍L=2B,顆粒表面摩擦系數(shù)μ=0.8。動(dòng)荷載作用下H=2B和H=3B的填料土拱率隨荷載幅值的變化(見圖8)。荷載幅值為0對(duì)應(yīng)活動(dòng)門B4下降8 mm,形成初始土拱,并未施加表面動(dòng)荷載的階段。從圖8可以看出,動(dòng)荷載作用下填料的土拱率均隨荷載幅值的增加而增大。在相同大小的動(dòng)荷載作用下,高度較高(即H=3B)填料的土拱率明顯小于高度較低(即H=2B)的土拱率,說明填料高度越高,土拱效應(yīng)更明顯。如圖8所示,隨著動(dòng)荷載幅值的增大,土拱效應(yīng)的退化越來越明顯。相同荷載幅值作用下,高度較高(H=3B)的土拱退化參數(shù)α均小于高度較低(H=2B)的情況,如荷載幅值為12 kPa時(shí),填料(H=2B)的土拱退化參數(shù)α為77%,H=3B的土拱退化參數(shù)α約為47%,前者是后者的1.6倍,說明增加填料高度可以明顯減少動(dòng)荷載作用導(dǎo)致的土拱效應(yīng)退化。

    圖8 土拱率、土拱退化參數(shù)隨動(dòng)荷載幅值的變化Fig.8 Variations of soil arching ratio and degradation parameter of soil arching with cyclic loading amplitudes

    4.3 荷載頻率

    選用4種頻率的動(dòng)荷載,分別為0.1、0.5、2.0、10.0 Hz,并與靜荷載(頻率為0 Hz)的結(jié)果進(jìn)行了對(duì)比。數(shù)值模型填料的高度H=2B,顆粒的表面摩擦系數(shù)μ=0.8,荷載幅值為8 kPa,加載范圍L=2B。

    如圖9所示,填料的土拱率隨動(dòng)荷載頻率的增大而增大,表明局部荷載作用下土拱效應(yīng)明顯退化,并且動(dòng)荷載作用下土拱效應(yīng)的退化程度大于靜荷載;當(dāng)荷載頻率為0時(shí),土拱退化參數(shù)α值最小也說明了這一點(diǎn)。此外,圖9表明隨著荷載頻率的增大,土拱退化程度越來越明顯,在荷載頻率為10 Hz時(shí),土拱效退化參數(shù)α為72%,約是靜荷載的2倍。

    圖9 土拱率、土拱退化參數(shù)隨動(dòng)荷載頻率的變化Fig.9 Variations of soil arching ratio and degradationparameter of soil arching with cyclic loading frequency

    5 結(jié) 論

    (1)土拱形成后,在局部動(dòng)荷載作用下土拱率增大,表明土拱效應(yīng)有一定的退化;在細(xì)觀上表現(xiàn)為接觸力的拱形結(jié)構(gòu)消失。

    (2)定義了土拱退化參數(shù)α以定量衡量土拱效應(yīng)的退化程度。隨著內(nèi)摩擦角的增大,土拱退化參數(shù)α的值逐漸減小,說明土拱效應(yīng)的退化程度降低,土拱抵御動(dòng)荷載的能力增強(qiáng)。

    (3)土拱退化參數(shù)α均隨荷載幅值、荷載頻率的增加而增大,表明土拱效應(yīng)的退化程度越來越嚴(yán)重。相同荷載幅值和頻率作用下,填料高度越高(即H=3B),土拱抵抗動(dòng)荷載導(dǎo)致退化的能力越強(qiáng)。

    猜你喜歡
    模型試驗(yàn)填料幅值
    反推力裝置模型試驗(yàn)臺(tái)的研制及驗(yàn)證
    基于S變換的交流電網(wǎng)幅值檢測(cè)系統(tǒng)計(jì)算機(jī)仿真研究
    電子制作(2017年7期)2017-06-05 09:36:13
    聚合物型色譜填料的研究進(jìn)展
    臺(tái)階式短加筋土擋墻行為特征的離心模型試驗(yàn)
    正序電壓幅值檢測(cè)及諧波抑制的改進(jìn)
    填料填裝工藝技術(shù)創(chuàng)新與應(yīng)用
    巨厚堅(jiān)硬巖漿巖不同配比的模型試驗(yàn)研究
    錯(cuò)流旋轉(zhuǎn)填料床脫除細(xì)顆粒物研究
    低壓電力線信道脈沖噪聲的幅值與寬度特征
    基于零序電壓幅值增量的消弧線圈調(diào)諧新方法
    小蜜桃在线观看免费完整版高清| 国内精品一区二区在线观看| 亚洲欧美精品自产自拍| 国产 一区精品| 国产黄a三级三级三级人| 六月丁香七月| 国产精品久久电影中文字幕| 白带黄色成豆腐渣| 少妇人妻一区二区三区视频| 国产精品女同一区二区软件| 免费在线观看成人毛片| 国产精品人妻久久久久久| 免费av不卡在线播放| 亚洲精品亚洲一区二区| 免费电影在线观看免费观看| av视频在线观看入口| 国产中年淑女户外野战色| 欧美国产日韩亚洲一区| 国产黄片美女视频| 久久这里只有精品中国| 欧美成人a在线观看| 啦啦啦啦在线视频资源| 国产亚洲精品综合一区在线观看| 日韩人妻高清精品专区| 亚州av有码| 午夜亚洲福利在线播放| 欧美3d第一页| 欧美激情国产日韩精品一区| 日韩精品青青久久久久久| 天堂动漫精品| 在线免费十八禁| 一本精品99久久精品77| 久久精品91蜜桃| 日韩制服骚丝袜av| 国产精品久久久久久av不卡| 亚洲美女搞黄在线观看 | 美女xxoo啪啪120秒动态图| 精品一区二区三区视频在线| 人妻丰满熟妇av一区二区三区| av天堂中文字幕网| 国产精品一及| 亚洲高清免费不卡视频| 精品人妻一区二区三区麻豆 | 国产av不卡久久| 亚洲七黄色美女视频| 国产视频内射| 国产精品99久久久久久久久| 日产精品乱码卡一卡2卡三| 国产av不卡久久| 美女cb高潮喷水在线观看| 99久久成人亚洲精品观看| 大香蕉久久网| 九九久久精品国产亚洲av麻豆| 男女之事视频高清在线观看| 精品久久久久久成人av| 国产精品一区二区三区四区免费观看 | 亚洲不卡免费看| 九九热线精品视视频播放| 亚洲欧美日韩无卡精品| 免费看a级黄色片| 人人妻,人人澡人人爽秒播| 2021天堂中文幕一二区在线观| 在线观看av片永久免费下载| av中文乱码字幕在线| 日韩一本色道免费dvd| 男女啪啪激烈高潮av片| 亚洲美女黄片视频| eeuss影院久久| 成人国产麻豆网| 免费电影在线观看免费观看| 国产成人福利小说| 简卡轻食公司| 亚洲性久久影院| 嫩草影视91久久| 99热只有精品国产| 偷拍熟女少妇极品色| 亚洲经典国产精华液单| 亚洲最大成人中文| 国产成人精品久久久久久| 久久精品91蜜桃| 亚洲av中文av极速乱| 精品乱码久久久久久99久播| 亚洲18禁久久av| 午夜福利高清视频| 欧美三级亚洲精品| 国产午夜福利久久久久久| 99视频精品全部免费 在线| 白带黄色成豆腐渣| 国产亚洲精品综合一区在线观看| 真人做人爱边吃奶动态| 国产色婷婷99| 日日摸夜夜添夜夜添av毛片| 国产精品亚洲一级av第二区| 国内少妇人妻偷人精品xxx网站| 国产蜜桃级精品一区二区三区| 亚洲精品在线观看二区| 日韩欧美三级三区| 精品乱码久久久久久99久播| 免费人成视频x8x8入口观看| 日韩 亚洲 欧美在线| 99久久成人亚洲精品观看| 国产亚洲精品av在线| 欧美人与善性xxx| av天堂中文字幕网| av在线蜜桃| 亚洲成人中文字幕在线播放| 国产 一区精品| 国产精品一及| 久久久久久国产a免费观看| 亚洲av熟女| 人妻少妇偷人精品九色| 国产aⅴ精品一区二区三区波| 成人欧美大片| 成人高潮视频无遮挡免费网站| 国产成人freesex在线 | 搡女人真爽免费视频火全软件 | 精品免费久久久久久久清纯| 99国产精品一区二区蜜桃av| 免费搜索国产男女视频| 日韩欧美精品免费久久| 最近在线观看免费完整版| 韩国av在线不卡| 亚洲天堂国产精品一区在线| 三级经典国产精品| 麻豆国产97在线/欧美| 熟女人妻精品中文字幕| 欧美日本视频| 精品一区二区免费观看| 亚洲精品国产av成人精品 | 狂野欧美激情性xxxx在线观看| 非洲黑人性xxxx精品又粗又长| 午夜免费男女啪啪视频观看 | 日本三级黄在线观看| 国产老妇女一区| av免费在线看不卡| 秋霞在线观看毛片| 国产精品一区二区三区四区免费观看 | 99在线视频只有这里精品首页| 国产探花在线观看一区二区| 美女大奶头视频| 夜夜爽天天搞| 最近最新中文字幕大全电影3| 欧美三级亚洲精品| 美女xxoo啪啪120秒动态图| 精品日产1卡2卡| 黄色日韩在线| 村上凉子中文字幕在线| 97超视频在线观看视频| 国产精品女同一区二区软件| 国产真实乱freesex| 麻豆精品久久久久久蜜桃| 国产精品永久免费网站| 波多野结衣巨乳人妻| 一个人看视频在线观看www免费| 夜夜看夜夜爽夜夜摸| 一区二区三区高清视频在线| 少妇丰满av| 国产成人福利小说| videossex国产| 韩国av在线不卡| 国产高清三级在线| 内射极品少妇av片p| 国产成人精品久久久久久| 久久婷婷人人爽人人干人人爱| 露出奶头的视频| 乱系列少妇在线播放| av在线亚洲专区| 亚洲va在线va天堂va国产| 女的被弄到高潮叫床怎么办| 亚洲成人av在线免费| 国产精品人妻久久久影院| 人妻少妇偷人精品九色| 精品福利观看| 久久人人爽人人爽人人片va| 极品教师在线视频| 悠悠久久av| 美女xxoo啪啪120秒动态图| 免费在线观看成人毛片| 看片在线看免费视频| 少妇高潮的动态图| 国产三级在线视频| 亚洲国产欧洲综合997久久,| av视频在线观看入口| 国产成人aa在线观看| 女的被弄到高潮叫床怎么办| 波多野结衣巨乳人妻| 成年女人看的毛片在线观看| 国产精品福利在线免费观看| АⅤ资源中文在线天堂| 91久久精品国产一区二区成人| 最近在线观看免费完整版| 亚洲人成网站在线播| av专区在线播放| 少妇高潮的动态图| 国产精品无大码| 亚洲无线在线观看| 麻豆乱淫一区二区| 一级毛片aaaaaa免费看小| 一区二区三区高清视频在线| 国产一区二区亚洲精品在线观看| 午夜福利视频1000在线观看| 午夜久久久久精精品| 午夜视频国产福利| 欧美区成人在线视频| 村上凉子中文字幕在线| 亚洲av成人av| 亚洲精品影视一区二区三区av| 日本黄色片子视频| 中文字幕av在线有码专区| 国产亚洲精品久久久久久毛片| 亚洲成人精品中文字幕电影| 亚洲熟妇熟女久久| 国产成人影院久久av| 中出人妻视频一区二区| 我的老师免费观看完整版| 久久久精品欧美日韩精品| 亚洲美女搞黄在线观看 | 韩国av在线不卡| 三级经典国产精品| 色哟哟·www| 99久国产av精品| 国产一级毛片七仙女欲春2| 欧美性感艳星| 美女被艹到高潮喷水动态| 久久精品国产自在天天线| 久久久色成人| 变态另类成人亚洲欧美熟女| 日韩欧美免费精品| a级一级毛片免费在线观看| 在线观看一区二区三区| 男插女下体视频免费在线播放| 国产在视频线在精品| 国产精品亚洲美女久久久| 精品一区二区三区视频在线观看免费| 3wmmmm亚洲av在线观看| 色综合站精品国产| 春色校园在线视频观看| 高清毛片免费看| 网址你懂的国产日韩在线| 日韩精品有码人妻一区| 久久天躁狠狠躁夜夜2o2o| 欧美中文日本在线观看视频| 午夜免费男女啪啪视频观看 | 一级毛片我不卡| 色av中文字幕| 欧美日本视频| 日韩人妻高清精品专区| 99热只有精品国产| 99久久无色码亚洲精品果冻| 麻豆av噜噜一区二区三区| 一级毛片aaaaaa免费看小| 国产午夜福利久久久久久| 免费av毛片视频| 亚洲欧美日韩卡通动漫| 亚洲成a人片在线一区二区| 精品久久久久久久人妻蜜臀av| 久久久久久九九精品二区国产| 日韩欧美免费精品| 亚洲人成网站在线观看播放| 亚洲综合色惰| 成年av动漫网址| 日韩精品青青久久久久久| 亚洲va在线va天堂va国产| 欧美日韩国产亚洲二区| 欧美性猛交黑人性爽| 人妻久久中文字幕网| 欧美最新免费一区二区三区| 3wmmmm亚洲av在线观看| 日本免费a在线| 欧美成人免费av一区二区三区| 国产精品人妻久久久久久| 日韩欧美 国产精品| 亚州av有码| 日本欧美国产在线视频| 国产精品一区二区免费欧美| 黄色配什么色好看| 午夜影院日韩av| 久久久精品94久久精品| 狂野欧美白嫩少妇大欣赏| 成人欧美大片| 一进一出抽搐gif免费好疼| 99热只有精品国产| 久久久色成人| 免费观看精品视频网站| 国产精品久久久久久亚洲av鲁大| 深爱激情五月婷婷| 精品午夜福利视频在线观看一区| 日韩精品中文字幕看吧| 国产精品国产三级国产av玫瑰| av女优亚洲男人天堂| 一本精品99久久精品77| 51国产日韩欧美| 国产精品电影一区二区三区| 午夜福利18| 99riav亚洲国产免费| 日韩一区二区视频免费看| 午夜影院日韩av| 插阴视频在线观看视频| 啦啦啦韩国在线观看视频| 国产三级在线视频| 国产精品一区二区性色av| 亚洲无线在线观看| 一级av片app| 欧美日韩综合久久久久久| 久久精品影院6| 久久精品国产99精品国产亚洲性色| 永久网站在线| 亚洲熟妇熟女久久| 亚洲av电影不卡..在线观看| av在线蜜桃| 最近中文字幕高清免费大全6| 日韩欧美国产在线观看| 久久热精品热| 久久久久久伊人网av| 美女cb高潮喷水在线观看| 成人一区二区视频在线观看| 人妻夜夜爽99麻豆av| 亚洲一级一片aⅴ在线观看| 日韩高清综合在线| 国产成人a区在线观看| 国产高潮美女av| 舔av片在线| aaaaa片日本免费| 免费看光身美女| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲第一电影网av| 寂寞人妻少妇视频99o| 欧美性猛交黑人性爽| 久久久国产成人精品二区| 久久鲁丝午夜福利片| 久久精品国产清高在天天线| 高清毛片免费看| 日本色播在线视频| 亚洲精品456在线播放app| 日产精品乱码卡一卡2卡三| 乱系列少妇在线播放| 国产亚洲精品久久久久久毛片| 一个人观看的视频www高清免费观看| 成人美女网站在线观看视频| 搡老岳熟女国产| 热99re8久久精品国产| 两个人的视频大全免费| 麻豆成人午夜福利视频| 国产精品综合久久久久久久免费| 精品一区二区三区视频在线| 国产精品综合久久久久久久免费| 可以在线观看的亚洲视频| 精品一区二区三区人妻视频| 国产淫片久久久久久久久| 国产成人精品久久久久久| 一级毛片我不卡| 波野结衣二区三区在线| 美女内射精品一级片tv| 精品不卡国产一区二区三区| 欧美xxxx黑人xx丫x性爽| 国产黄片美女视频| 国模一区二区三区四区视频| 在线观看免费视频日本深夜| 午夜免费激情av| 国产老妇女一区| 国产精品久久视频播放| 插逼视频在线观看| 99久久精品国产国产毛片| 国产亚洲精品久久久久久毛片| 欧美xxxx黑人xx丫x性爽| 婷婷精品国产亚洲av| 亚洲aⅴ乱码一区二区在线播放| 国产精品乱码一区二三区的特点| 一夜夜www| 国产免费男女视频| 俄罗斯特黄特色一大片| av在线亚洲专区| 成人漫画全彩无遮挡| 亚洲成av人片在线播放无| 97在线视频观看| 亚洲av中文字字幕乱码综合| 久久久精品欧美日韩精品| 97超碰精品成人国产| 国国产精品蜜臀av免费| 日韩欧美在线乱码| 观看美女的网站| 91狼人影院| 综合色丁香网| 少妇熟女aⅴ在线视频| а√天堂www在线а√下载| 欧美绝顶高潮抽搐喷水| 校园人妻丝袜中文字幕| 亚洲av五月六月丁香网| 日韩欧美精品免费久久| 不卡一级毛片| 最近2019中文字幕mv第一页| 91午夜精品亚洲一区二区三区| 三级毛片av免费| 波多野结衣高清无吗| 在线看三级毛片| 欧美日本亚洲视频在线播放| 中出人妻视频一区二区| 男女做爰动态图高潮gif福利片| 亚洲不卡免费看| 99久久精品热视频| 又黄又爽又刺激的免费视频.| 亚洲自偷自拍三级| 真人做人爱边吃奶动态| 偷拍熟女少妇极品色| 搞女人的毛片| 一个人观看的视频www高清免费观看| h日本视频在线播放| 中国美女看黄片| 九九在线视频观看精品| 麻豆成人午夜福利视频| 婷婷精品国产亚洲av在线| 97超级碰碰碰精品色视频在线观看| 午夜影院日韩av| 老司机午夜福利在线观看视频| 黄片wwwwww| 精品久久久久久成人av| 久久中文看片网| 黄色欧美视频在线观看| 欧美一区二区国产精品久久精品| 亚洲性久久影院| 日本a在线网址| 伦理电影大哥的女人| 欧美xxxx黑人xx丫x性爽| av.在线天堂| 日本欧美国产在线视频| 欧美+亚洲+日韩+国产| 国产视频内射| 最新中文字幕久久久久| 久久久a久久爽久久v久久| 91麻豆精品激情在线观看国产| 老司机午夜福利在线观看视频| 黄色一级大片看看| 小说图片视频综合网站| 淫秽高清视频在线观看| 欧美国产日韩亚洲一区| 2021天堂中文幕一二区在线观| 免费搜索国产男女视频| 搡老熟女国产l中国老女人| 亚州av有码| 又爽又黄无遮挡网站| 日本-黄色视频高清免费观看| 免费看美女性在线毛片视频| 一级毛片aaaaaa免费看小| 91麻豆精品激情在线观看国产| 51国产日韩欧美| 亚洲人成网站在线播放欧美日韩| 丰满人妻一区二区三区视频av| 成人三级黄色视频| 男人和女人高潮做爰伦理| 老熟妇乱子伦视频在线观看| 免费av毛片视频| 亚洲成a人片在线一区二区| .国产精品久久| 观看美女的网站| 午夜精品一区二区三区免费看| 久久99热这里只有精品18| 久久人人爽人人片av| 网址你懂的国产日韩在线| 国产精品精品国产色婷婷| 蜜桃久久精品国产亚洲av| 最好的美女福利视频网| 日本欧美国产在线视频| 免费无遮挡裸体视频| 欧美在线一区亚洲| 午夜日韩欧美国产| 国产亚洲91精品色在线| 99热这里只有是精品在线观看| 1000部很黄的大片| 97热精品久久久久久| 亚洲最大成人中文| 激情 狠狠 欧美| 亚洲人成网站高清观看| 国产伦精品一区二区三区视频9| 少妇熟女欧美另类| 国产女主播在线喷水免费视频网站 | 精品一区二区三区视频在线| 全区人妻精品视频| 99在线视频只有这里精品首页| 亚洲精品在线观看二区| 欧洲精品卡2卡3卡4卡5卡区| 久久精品影院6| 中文资源天堂在线| 国产精品久久视频播放| 国产精品综合久久久久久久免费| 欧美3d第一页| 日日撸夜夜添| 99久久精品一区二区三区| 天美传媒精品一区二区| 全区人妻精品视频| 久久久久性生活片| 久久婷婷人人爽人人干人人爱| 中文字幕久久专区| 日韩成人av中文字幕在线观看 | 国产真实乱freesex| 午夜福利18| 亚洲第一电影网av| 久久精品综合一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 日韩高清综合在线| 久久久午夜欧美精品| 亚洲av电影不卡..在线观看| 亚洲熟妇中文字幕五十中出| 中文资源天堂在线| 在线观看午夜福利视频| 精品人妻偷拍中文字幕| 校园春色视频在线观看| 哪里可以看免费的av片| 国产免费男女视频| 久久国内精品自在自线图片| 午夜福利在线在线| 国产av一区在线观看免费| 精品久久国产蜜桃| 黄色一级大片看看| 99在线视频只有这里精品首页| 午夜福利18| 老师上课跳d突然被开到最大视频| 97在线视频观看| 国产精品一区二区三区四区免费观看 | 欧美+日韩+精品| 嫩草影院入口| 亚洲激情五月婷婷啪啪| 午夜福利成人在线免费观看| 哪里可以看免费的av片| 看免费成人av毛片| 久久久国产成人免费| 亚洲美女黄片视频| 丝袜喷水一区| 超碰av人人做人人爽久久| 成人亚洲欧美一区二区av| 午夜爱爱视频在线播放| avwww免费| 成人一区二区视频在线观看| av视频在线观看入口| 国产精品一区二区性色av| 免费看日本二区| 97超视频在线观看视频| 欧美人与善性xxx| 99久久精品国产国产毛片| 亚洲乱码一区二区免费版| 久久久久久久久久久丰满| 九九热线精品视视频播放| 国语自产精品视频在线第100页| 乱系列少妇在线播放| 国内揄拍国产精品人妻在线| 欧美最黄视频在线播放免费| 欧美zozozo另类| 国产女主播在线喷水免费视频网站 | 成年女人看的毛片在线观看| 干丝袜人妻中文字幕| 免费不卡的大黄色大毛片视频在线观看 | 日韩欧美一区二区三区在线观看| 中文字幕av在线有码专区| 国产精品久久视频播放| 国产av在哪里看| 精品人妻一区二区三区麻豆 | 国产色爽女视频免费观看| 国产精品女同一区二区软件| 日本-黄色视频高清免费观看| 成年免费大片在线观看| 国产一区亚洲一区在线观看| 亚洲性久久影院| 99久国产av精品国产电影| 97人妻精品一区二区三区麻豆| 成人综合一区亚洲| 丝袜美腿在线中文| 嫩草影视91久久| 级片在线观看| 中国国产av一级| 欧美日韩一区二区视频在线观看视频在线 | 亚洲国产精品成人综合色| 色哟哟哟哟哟哟| 欧美极品一区二区三区四区| 国产精品日韩av在线免费观看| 99久久九九国产精品国产免费| 亚洲成人久久爱视频| 九九久久精品国产亚洲av麻豆| 欧美高清成人免费视频www| 18禁在线无遮挡免费观看视频 | 高清毛片免费看| 日韩强制内射视频| av.在线天堂| 欧美成人精品欧美一级黄| 国产黄a三级三级三级人| 亚洲经典国产精华液单| 麻豆精品久久久久久蜜桃| 亚洲av免费在线观看| 亚洲国产精品合色在线| 51国产日韩欧美| 国产私拍福利视频在线观看| 一个人看视频在线观看www免费| 国产淫片久久久久久久久| 国产精品av视频在线免费观看| 国产在线男女| 高清毛片免费观看视频网站| 国产一区二区三区av在线 | 不卡一级毛片| 精品99又大又爽又粗少妇毛片| 综合色av麻豆| 在线天堂最新版资源| 成年女人毛片免费观看观看9| 有码 亚洲区| 在线观看免费视频日本深夜| 成年女人毛片免费观看观看9| 国产一区二区在线av高清观看| 免费无遮挡裸体视频| 国产视频内射| 久久精品影院6| 欧美最新免费一区二区三区| 欧美3d第一页| 国产一区二区在线av高清观看| 精品不卡国产一区二区三区| 老司机午夜福利在线观看视频| 男女边吃奶边做爰视频| 熟女电影av网| 亚洲av电影不卡..在线观看|