亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        波浪循環(huán)荷載作用下濱海軟土水平向動(dòng)力特性試驗(yàn)研究

        2016-04-06 07:06:38郝斌趙玉成劉珍巖李占嶺盧二巖
        土木建筑與環(huán)境工程 2016年1期

        郝斌 趙玉成 劉珍巖 李占嶺 盧二巖

        摘要:為研究海上風(fēng)電樁基在波浪荷載作用下,產(chǎn)生水平向循環(huán)荷載對樁基周圍土體動(dòng)力特性的影響,以唐山地區(qū)濱海軟土為研究對象,通過室內(nèi)動(dòng)三軸試驗(yàn),研究不同圍壓、動(dòng)應(yīng)力幅值和振動(dòng)次數(shù)條件下對軟土水平向動(dòng)力特性的影響。結(jié)果表明:軟土水平向動(dòng)強(qiáng)度隨圍壓增加而增加,隨振動(dòng)次數(shù)增加而減??;動(dòng)應(yīng)力幅值增大,破壞振次減小;水平向動(dòng)應(yīng)變?chǔ)興隨振動(dòng)次數(shù)增加變大,且動(dòng)應(yīng)力幅值越大,增速越明顯,變化規(guī)律遵循Monismith模型;動(dòng)應(yīng)力幅值改變時(shí),軟土水平向動(dòng)模量變化明顯,當(dāng)圍壓減小,動(dòng)彈性模量減?。徊苠檐浲了较蜷g具有明顯的結(jié)構(gòu)性,不同圍壓條件下,隨動(dòng)應(yīng)力幅值增加動(dòng)阻尼比均表現(xiàn)增大趨勢。

        關(guān)鍵詞:濱海軟土;循環(huán)荷載;水平向動(dòng)力特性;動(dòng)三軸試驗(yàn)

        中圖分類號(hào):TU411.8

        文獻(xiàn)標(biāo)志碼:A文章編號(hào):16744764(2016)01011606

        Abstract:

        A series of dynamic triaxle test on soft soil in Tangshan Binhai area was carried out to understand the horizontal dynamic characteristics of the soil around the offshore wind power pile foundation under wave load. The influence of confining pressure, dynamic stress amplitude and vibration frequency on horizontal dynamic characteristics of soft soil were studied. The results show that the horizontal dynamic strength of soft soil increases with the increase of confining pressure increases,and decreases with the increase of vibration. When the dynamic stress amplitude increases, the vibration times decreases. The horizontal dynamic strain εd increases with the increase of vibration frequency, and the larger the dynamic stress amplitude is, the more significant the growth is. The change agrees well with that of the Monismith model. The change of the horizontal dynamic modulus of soft soil is significantly influenced by the dynamic stress amplitude. When confining pressure decreases, dynamic modulus decreases. The Caofeidian soft soil has distinct structural horizontal dynamic characteristics. Under different confining pressures, the damping ratio shows an increasing trend when dynamic stress amplitude increases.

        Keywords:coastal soft soil; cyclic loading; horizontal dynamic characteristics; dynamic triaxial test

        風(fēng)電產(chǎn)業(yè)迅速發(fā)展,受陸地條件制約,地廣人稀的海岸灘涂和海上風(fēng)電產(chǎn)業(yè)得到重視,海上風(fēng)電場建設(shè)步入發(fā)展的快車道,但沿海地區(qū)廣泛分布著工程性質(zhì)較差的軟土。風(fēng)電設(shè)施建成后樁基周圍土體受到常見的動(dòng)力循環(huán)荷載——波浪荷載[1]長期作用,在波浪荷載作用下土體產(chǎn)生的累積變形、強(qiáng)度及穩(wěn)定性問題成為眾多學(xué)者關(guān)注的重點(diǎn)[2]。波浪荷載作為一種作用時(shí)間長、振動(dòng)較為規(guī)則的循環(huán)荷載,需考慮循環(huán)效應(yīng)對樁基周圍土體動(dòng)力特性的影響。應(yīng)用GDS室內(nèi)動(dòng)三軸試驗(yàn)系統(tǒng)研究波浪荷載作用下軟土水平向的動(dòng)力特性。

        對于海上風(fēng)電建設(shè)的研究,最為核心的是樁基及周圍土體在波浪長期循環(huán)荷載作用下分析土體動(dòng)力特性變化規(guī)律和運(yùn)用PY曲線[3]研究樁基側(cè)向位移。近年來,關(guān)于循環(huán)動(dòng)荷載樁的研究受到關(guān)注,Basack[4]通過建立波浪產(chǎn)生的水平向循環(huán)荷載對于土體強(qiáng)度和剛度的衰減函數(shù)模型,應(yīng)用于樁基承載力分析,提出基于該地區(qū)波浪循環(huán)荷載作用下軟土動(dòng)力特性的樁基設(shè)計(jì)參考; Basak[5]在軟土地基上進(jìn)行了一系列室內(nèi)試驗(yàn),研究海洋中波浪荷載作用下軟土地基中樁基周圍水平向循環(huán)荷載對樁基承載力影響,以樁的變形為研究重點(diǎn),缺少波浪產(chǎn)生的水平向循環(huán)荷載次數(shù)、動(dòng)應(yīng)力幅值及頻率變化對土體動(dòng)力特性深入的分析;Banerjee等[6]通過室內(nèi)動(dòng)三軸試驗(yàn)獲得波浪水平向循環(huán)荷載作用下軟土剛度衰減模型,應(yīng)用ABAQUS有效模擬復(fù)雜樁土系統(tǒng)進(jìn)一步了解樁基周圍土體動(dòng)力特性對樁的影響。中國眾多學(xué)者對飽和軟土在交通循環(huán)荷載及地震循環(huán)荷載作用下產(chǎn)生的動(dòng)力特性進(jìn)行了大量試驗(yàn)與學(xué)術(shù)研究,形成了相對完備的理論[710],而對于土體在波浪循環(huán)荷載長期作用下水平向動(dòng)力特性的研究較少涉及。曾向軍等[7]以洞庭湖區(qū)沉積軟土為研究對象,應(yīng)用動(dòng)三軸試驗(yàn),研究不同動(dòng)應(yīng)力幅值、振動(dòng)頻率和圍壓下湖湘軟土在豎向荷載作用時(shí)動(dòng)力特性,表明軸向累積應(yīng)變受動(dòng)應(yīng)力幅值、圍壓及振動(dòng)頻率影響呈穩(wěn)定型增長,動(dòng)孔隙水壓受動(dòng)應(yīng)力幅值影響;張向東等[8]以營口軟土為研究對象,通過動(dòng)三軸試驗(yàn)研究豎向荷載作用下動(dòng)力特性變化及由此導(dǎo)致的地基承載力不足等問題,試驗(yàn)結(jié)果表明營口地區(qū)軟土具有明顯的結(jié)構(gòu)性,其動(dòng)應(yīng)變振次曲線和動(dòng)骨干曲線存在轉(zhuǎn)折點(diǎn);曹勇等[9]對海積結(jié)構(gòu)軟土動(dòng)力特性與變形特征進(jìn)行分析,在3種波形循環(huán)荷載作用下考慮振動(dòng)次數(shù)、幅值及頻率的影響,并類比不同循環(huán)荷載波形下軟土剛度軟化過程,得出軟化指數(shù)隨振動(dòng)次數(shù)變化規(guī)律;曹洋等[10]以杭州某工地基坑結(jié)構(gòu)性軟粘土為研究對象,應(yīng)用空心圓柱扭剪儀,進(jìn)行不同循環(huán)應(yīng)力比和不同頻率下的主應(yīng)力軸旋轉(zhuǎn)循環(huán)剪切試驗(yàn),模擬波浪荷載作用下原狀軟粘土孔隙水壓累積、模量衰減等動(dòng)力特性,確立土體應(yīng)變破壞標(biāo)準(zhǔn)及土體強(qiáng)度變化規(guī)律;郭玉樹等[11]應(yīng)用循環(huán)三軸試驗(yàn)分析海上風(fēng)力發(fā)電機(jī)單樁基礎(chǔ)側(cè)向位移,針對砂性土壤進(jìn)行室內(nèi)循環(huán)三軸試驗(yàn),得到塑性應(yīng)變增量推求土壤割線衰減剛度,導(dǎo)入三維有限元數(shù)值模型計(jì)算可得樁身在循環(huán)側(cè)向力下的位移;蔡袁強(qiáng)等[12]以蕭山正常固結(jié)飽和軟粘土為研究對象,進(jìn)行豎向應(yīng)力控制的循環(huán)三軸試驗(yàn),表明初始偏應(yīng)力對動(dòng)彈模量及阻尼比有較大影響,在試驗(yàn)數(shù)據(jù)的基礎(chǔ)上建立了蕭山地區(qū)飽和軟粘土剛度軟化規(guī)律經(jīng)驗(yàn)擬合公式。

        筆者針對波浪循環(huán)荷載作用下樁基周圍土體,以唐山地區(qū)曹妃甸周邊軟粘土為研究對象進(jìn)行水平向動(dòng)三軸試驗(yàn),研究不同圍壓、密度和動(dòng)應(yīng)力幅值條件下軟土的水平向動(dòng)強(qiáng)度和動(dòng)模量變化,分析圍壓、動(dòng)應(yīng)力幅值和振動(dòng)次數(shù)對軟土動(dòng)強(qiáng)度、動(dòng)變形、動(dòng)模量和動(dòng)阻尼比特性的影響,將試驗(yàn)結(jié)果與現(xiàn)有的研究成果進(jìn)行對比分析,總結(jié)唐山地區(qū)濱海軟土在波浪荷載作用下引起的軟土水平向累積應(yīng)變規(guī)律和動(dòng)力特性,獲得唐山地區(qū)濱海軟土工程性質(zhì)。

        1唐山濱海地區(qū)軟土水平向動(dòng)三軸

        試驗(yàn)

        1.1動(dòng)三軸試驗(yàn)系統(tǒng)與土樣制備

        試驗(yàn)應(yīng)用GDSDYNTTS電機(jī)控制的動(dòng)三軸試驗(yàn)系統(tǒng),可對直徑39.1 mm,高度80 mm的三軸試樣進(jìn)行試驗(yàn)。

        試樣取自唐山地區(qū)曹妃甸周邊的軟粘土,對原狀土樣進(jìn)行常規(guī)室內(nèi)土工試驗(yàn),得到軟土的基本物理力學(xué)指標(biāo),見表1。

        1.2動(dòng)三軸試驗(yàn)方案

        對原狀土樣水平向取樣,制備12個(gè)原狀試驗(yàn)土樣,采用GDS動(dòng)三軸系統(tǒng),進(jìn)行軟土水平向動(dòng)力特性研究。試驗(yàn)方案見表2,不同圍壓反映不同埋深。試驗(yàn)前充分飽和試樣,應(yīng)用GDS動(dòng)三軸系統(tǒng)中BCheck模塊對土樣進(jìn)行檢測,飽和度均達(dá)到97%以上。試驗(yàn)固結(jié)比取Kc=1.0,振動(dòng)頻率為f=1 Hz(每組循環(huán)采集10個(gè)點(diǎn)),采用全幅應(yīng)變=5%作為原狀土樣破壞標(biāo)準(zhǔn)。

        2試驗(yàn)結(jié)果分析與數(shù)據(jù)處理

        2.1軟土水平向動(dòng)強(qiáng)度特性分析

        根據(jù)εd達(dá)到5%時(shí)作為土樣破壞標(biāo)準(zhǔn),繪制出軟土水平向抗剪強(qiáng)度τd和破壞振次Nf的關(guān)系曲線,即τdlgN曲線,研究長期循環(huán)荷載作用下軟土水平向抗剪強(qiáng)度變化規(guī)律。

        圖1所示為相同固結(jié)比、振動(dòng)頻率時(shí)不同圍壓和動(dòng)應(yīng)力幅值條件下唐山濱海地區(qū)軟土水平向的動(dòng)剪強(qiáng)度曲線。表明曹妃甸軟土水平向動(dòng)剪應(yīng)力隨振動(dòng)次數(shù)增加而減小,對試驗(yàn)數(shù)據(jù)進(jìn)行擬合分析,二者關(guān)系可用負(fù)冪次函數(shù)表示

        τd=aN-b

        (1)

        即曹妃甸軟土水平向動(dòng)強(qiáng)度特性與張向東等[8]得到的營口軟土動(dòng)強(qiáng)度特性中動(dòng)剪應(yīng)力振動(dòng)次數(shù)函數(shù)關(guān)系一致。式(1)中τd為水平向動(dòng)剪應(yīng)力;N為振動(dòng)次數(shù);a和b為試驗(yàn)條件下相關(guān)參數(shù)。

        由圖1可知,軟土水平向動(dòng)抗剪強(qiáng)度隨圍壓增加而增大,原因在于軟土顆粒之間孔隙隨圍壓增大而被擠密。此外,圍壓大小影響土顆粒之間相互作用,隨圍壓減小,土顆粒間相互作用減弱,動(dòng)荷載對土體破壞增強(qiáng),土體抗剪強(qiáng)度降低。

        進(jìn)行擬合分析得表3,擬合度均達(dá)0.99以上。隨a值降低,動(dòng)剪應(yīng)力降低;隨b值增大,動(dòng)剪應(yīng)力減小。為合理評價(jià)循環(huán)荷載對軟土水平向動(dòng)強(qiáng)度影響,需進(jìn)一步考慮不同固結(jié)比、振動(dòng)頻率及初始偏應(yīng)力條件下對水平向動(dòng)剪應(yīng)力的影響。

        2.2軟土水平向累積動(dòng)應(yīng)變特性分析

        圖2為不同圍壓條件下濱海地區(qū)的軟土水平向累積動(dòng)應(yīng)變與振動(dòng)次數(shù)關(guān)系曲線。

        當(dāng)動(dòng)應(yīng)力較小時(shí),對于應(yīng)變模型的研究,許多學(xué)者結(jié)合試驗(yàn)數(shù)據(jù)擬合得到的關(guān)系曲線,提出了相關(guān)經(jīng)驗(yàn)公式,如Lentz[13]的對數(shù)關(guān)系式、Monismith等[14]的指數(shù)關(guān)系式等。模型中應(yīng)用最為廣泛的為Monismith指數(shù)模型,即

        圖3表明曹妃甸軟土水平向動(dòng)彈性模量受動(dòng)應(yīng)變及圍壓影響顯著。動(dòng)應(yīng)變增加,水平向動(dòng)彈性模量減??;圍壓減小,水平向動(dòng)彈性模量隨之減小。不同圍壓條件下,隨動(dòng)應(yīng)力幅值增加,動(dòng)應(yīng)變增大,動(dòng)彈性模量整體趨勢減小,在于動(dòng)應(yīng)力幅值的增加對土體結(jié)構(gòu)破壞加強(qiáng),使土體水平向動(dòng)應(yīng)變增大。

        2.4軟土水平向動(dòng)阻尼比特性分析

        圖4為不同試驗(yàn)條件下曹妃甸軟土動(dòng)阻尼比試驗(yàn)曲線。動(dòng)三軸試驗(yàn)測定的動(dòng)阻尼比λd表征為每振動(dòng)一周中土體能量的耗散,定義為實(shí)際阻尼系數(shù)C與臨界阻尼系數(shù)Ccr之比

        圖4(a)、(b)、(c)表明當(dāng)圍壓取值范圍在50~150 kPa時(shí),曹妃甸軟土水平向動(dòng)阻尼比介于0.2~06之間,隨動(dòng)應(yīng)力幅值增大而減小。由圖4(d)可知,圍壓450 kPa時(shí),阻尼比介于0.3~0.8之間,可知隨圍壓大幅度增加阻尼比數(shù)值增大,動(dòng)阻尼比增大即荷載作用過程中軟土能量的損耗增加。圍壓增大時(shí),軟土顆粒之間的聯(lián)結(jié)作用加強(qiáng),對波浪循環(huán)荷載振動(dòng)形式形成更大干擾,其振幅隨作用時(shí)間的增長而衰減,且衰減越快,阻尼越大。不同圍壓條件下,動(dòng)應(yīng)力幅值的增加,動(dòng)阻尼比均表現(xiàn)為增大趨勢。

        2.5波浪荷載大小的變化分析

        綜合研究表征波浪荷載的動(dòng)荷因素:動(dòng)荷波形、動(dòng)荷頻率、動(dòng)荷幅值,可獲得動(dòng)彈性模量隨波浪循環(huán)荷載變化規(guī)律。如圖5所示。

        圖5表明隨波浪荷載大小的改變對不同圍壓、密度條件下軟土的動(dòng)彈性模量影響顯著。分析由表5得到的最大彈性模量Edmax和最大剪切彈性模量Gdmax表明,隨波浪荷載改變,不同圍壓條件下,隨波浪荷載振幅增大,Edmax呈增大趨勢。

        3結(jié)論

        1)曹妃甸軟土水平向動(dòng)強(qiáng)度隨圍壓的增加而增大,隨振動(dòng)次數(shù)增加而減小。當(dāng)動(dòng)應(yīng)力幅值減小,破壞振次增加。波浪循環(huán)動(dòng)荷載長期作用下,動(dòng)荷載頻率成為影響軟土動(dòng)強(qiáng)度的主要因素,為更合理評價(jià)軟土水平動(dòng)力特性,需進(jìn)一步考慮不同固結(jié)比、振動(dòng)頻率等因素的影響。

        2)曹妃甸軟土的水平向動(dòng)應(yīng)變隨振動(dòng)次數(shù)增加而增大。初始階段,動(dòng)應(yīng)變增速較小隨后逐漸增大,直至土體發(fā)生破壞;受動(dòng)應(yīng)力幅值影響,幅值越大增速越明顯。對試驗(yàn)數(shù)據(jù)進(jìn)行回歸擬合分析,該地區(qū)軟土水平向累積變形規(guī)律適用于Monismith指數(shù)模型。

        3)曹妃甸軟土水平向動(dòng)彈模量受動(dòng)應(yīng)力幅值的影響。隨動(dòng)應(yīng)力幅值增大,水平向動(dòng)應(yīng)變增大,動(dòng)彈模量減?。浑S圍壓減小,動(dòng)彈性模量隨之減小。對試驗(yàn)數(shù)據(jù)進(jìn)行回歸擬合分析,水平向動(dòng)彈模量變化可近似表示為雙曲線模型。

        4)阻尼比λd與動(dòng)應(yīng)變關(guān)系曲線表明曹妃甸軟土水平向具有顯著的結(jié)構(gòu)性,圍壓較大時(shí)(450 kPa)隨土顆粒擠密孔隙減小,土顆粒間聯(lián)結(jié)作用增強(qiáng),能量消耗增加,阻尼比較大。隨圍壓減小,軟土水平向間作用力減小,阻尼比降低。不同圍壓條件下,隨動(dòng)應(yīng)力幅值的增加,阻尼比均表現(xiàn)增大的趨勢。

        5)綜合分析不同波浪荷載形式作用下,土體Edmax和Gdmax受波浪荷載的變化影響顯著,隨波浪荷載增大,Edmax和Gdmax呈增大趨勢。

        參考文獻(xiàn):

        [1] 葛川, 何炎平, 葉宇,等. 海上風(fēng)電場的發(fā)展、構(gòu)成和基礎(chǔ)形式[J]. 中國海洋平臺(tái), 2008, 23(6): 3135.

        Ge C, He Y P, Ye Y, et al. Development composing and foundation of offshore windfarm [J]. China Offshore Platform, 2008, 23(6): 3135. (in Chinese)

        [2] 畢雪梅, 趙瑞斌, 曲國勝. 循環(huán)荷載下天津吹填土動(dòng)應(yīng)力應(yīng)變關(guān)系研究[J]. 工程勘察, 2014(1): 15.

        Bi X M, Zhao R B, Qu G S. Experimental study on dynamic stressstrain relation of reclaimed soil in Tianjin subjected to cyclic loads [J].Geotechnical Investigation & Surveying, 2014 (1): 15. (in Chinese)

        [3] 胡安峰, 張光建, 賈玉帥, 等. 剛度衰減模型在大直徑樁累積側(cè)向位移分析中的應(yīng)用[J]. 浙江大學(xué)學(xué)報(bào):工學(xué)版, 2014, 48(4): 721726.

        Hu A F, Zhang G J, Jia Y S, et al. Application of degradation stiffness model in analysis of cumulative lateral displacement of monopile foundation [J]. Journal of Zhejiang University: Engineering Science, 2014, 48(4): 721726. (in Chinese)

        [4] Basack S. Design recommendations for pile subjected to cyclic load [J]. Marine Georesources & Geotechnology, 2013, 33(4): 356360.

        [5] Basack S. Response of vertical pile group subjected to horizontal cyclic load in soft clay [J]. Latin American Journal of Solids and Structures, 2010(7): 91103.

        [6] Banerjee S, Shirole O N. Numerical analysis of piles under cyclic lateral load [J]. Indian Geotechnical Journal, 2014, 44(4): 436448.

        [7] 曾向軍, 鄧宗偉, 高乾豐, 等. 循環(huán)荷載下湖相軟土動(dòng)力特性研究[J]. 工業(yè)建筑, 2015, 45(3): 109114.

        Zeng X J, Deng Z W, Gao Q F, et al. Study of dynamic characteristics of lacustrine soft soil under cyclic loading [J]. Industrial Construction, 2015, 45(3): 109114. (in Chinese)

        [8] 張向東, 劉家順, 張虎偉. 循環(huán)荷載作用下軟土動(dòng)力特性試驗(yàn)研究[J]. 公路交通科技, 2014, 31(5): 17.

        Zhang X D, Liu J S, Zhang H W. Experimental study on dynamic characteristic of soft soil under cyclic loading [J]. Journal of Highway and Transportation Research and Development, 2014, 31(5): 17. (in Chinese)

        [9] 曹勇, 孔令偉, 楊愛武. 海積結(jié)構(gòu)性軟土動(dòng)力性狀的循環(huán)荷載波形效應(yīng)與剛度軟化特征[J]. 巖土工程學(xué)報(bào), 2013, 35(3): 583589.

        Cao Y, Kong L W, Yang A W. Waveform effect of cyclic loading of dynamic character and stiffness degradation characteristics of marine deposited natural soft clay [J]. Chinese Journal of Geotechnical Engineering, 2013, 35(3): 583589. (in Chinese)

        [10] 曹洋, 周建, 嚴(yán)佳佳. 原狀軟黏土主應(yīng)力軸動(dòng)態(tài)旋轉(zhuǎn)試驗(yàn)研究[J]. 中南大學(xué)學(xué)報(bào):自然科學(xué)版, 2014, 45(3): 902909.

        Cao Y, Zhou J, Yan J J. Dynamic principal stress rotation test on intact soft clay [J]. Journal of Central South University: Science and Technology, 2014, 45(3): 902909. (in Chinese)

        [11] 郭玉樹, 亞克慕斯·馬丁, 阿布達(dá)雷赫曼·哈里. 用循環(huán)三軸試驗(yàn)分析海上風(fēng)力發(fā)電機(jī)單樁基礎(chǔ)側(cè)向位移[J]. 巖土工程學(xué)報(bào), 2009, 31(11): 17291734.

        Kuo Y S, Achmus M, AbdelRahman K. Estimation of lateral deformation of monopile foundations by use of cyclic triaxial tests [J]. Chinese Journal of Geotechnical Engineering, 2009, 31(11): 17291734. (in Chinese)

        [12] 蔡袁強(qiáng), 王軍, 徐長節(jié). 初始偏應(yīng)力作用對蕭山軟粘土動(dòng)彈模量與阻尼影響試驗(yàn)研究[J]. 巖土力學(xué), 2007, 28(11): 22912302.

        Cai Y Q, Wang J, Xu C J. Experimental study on dynamic elastic and damping ratio of Xiaoshan saturated soft clay considering initial deviator stress [J]. Rock and Soil Mechanics, 2007, 28(11): 22912302. (in Chinese)

        [13] Lentz R W. Permanent deformation of a cohesionless subgrade material under cyclic loading [D]. East Lansing:Michigan State University, 1979.

        [14] Monismith C L, Ogawa N, Freeme C R. Permanent deformation characteristics of subsoil due to repeated loading [J]. Transportation Research Record, 1975, (537): 117.

        [15] 魏星, 黃茂松. 交通荷載作用下公路軟土地基長期沉降的計(jì)算[J]. 巖土力學(xué), 2009, 30(11): 33423346.

        Wei X, Huang M S. A simple method to predict trafficloadinduced permanent settlement of road on soft subsoil [J]. Rock and Soil Mechanics, 2009, 30(11): 33423346. (in Chinese)

        (編輯胡英奎)

        领导边摸边吃奶边做爽在线观看 | 丝袜美腿久久亚洲一区| 最新在线观看免费的a站国产| 在线观看的网站| 高中生粉嫩无套第一次| 亚洲欧美日韩高清一区二区三区| 亚洲永久精品日韩成人av| 激情综合婷婷色五月蜜桃| 国产一区二区精品久久| 久久久久国产一级毛片高清版A| av最新版天堂在资源在线| 美丽的小蜜桃在线观看| 亚洲热线99精品视频| 国产内射在线激情一区| 久久精品成人91一区二区| 加勒比一区二区三区av| 欧美老妇牲交videos| 国自产精品手机在线观看视频| 岛国AV一区二区三区在线观看| 亚洲AV无码中文AV日韩A| 81久久免费精品国产色夜| 无码人妻久久一区二区三区蜜桃| 精品亚洲欧美无人区乱码| 人妻少妇看A偷人无码电影| 亚洲永久精品日韩成人av| 亚洲av无码一区二区三区网址| 亚洲人成绝费网站色www| 亚洲都市校园激情另类| 国产精品麻豆一区二区三区| 亚洲理论电影在线观看| 精品熟女少妇av免费观看| 亚洲AV手机专区久久精品| 一区二区三区国产内射 | 99久久精品人妻一区| 久久精品夜色噜噜亚洲a∨| 国产熟人av一二三区| 一本一本久久久久a久久综合激情| 少妇人妻出水中文字幕乱码| 国产综合精品久久99之一| 青楼妓女禁脔道具调教sm| 欧美日韩亚洲一区二区精品|