999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

深基坑雙排樁結(jié)構(gòu)受力變形機(jī)理及其計(jì)算模型研究

2016-08-01 00:49:30薛光橋
關(guān)鍵詞:有限元法

薛光橋

(1.中鐵第四勘察設(shè)計(jì)院集團(tuán)有限公司,武漢 430063;2.水下隧道技術(shù)湖北省工程實(shí)驗(yàn)室,武漢 430063)

?

深基坑雙排樁結(jié)構(gòu)受力變形機(jī)理及其計(jì)算模型研究

薛光橋1,2

(1.中鐵第四勘察設(shè)計(jì)院集團(tuán)有限公司,武漢430063;2.水下隧道技術(shù)湖北省工程實(shí)驗(yàn)室,武漢430063)

摘要:近年來,為了適應(yīng)國內(nèi)城市地下工程建設(shè)蓬勃發(fā)展帶來的地下工程密集化、復(fù)雜化等方面的需要,深基坑開挖時雙排樁圍護(hù)結(jié)構(gòu)因具有側(cè)向剛度大、控制變形能力強(qiáng)、施工快速等特點(diǎn)在各種大型深基坑工程中得到了較為廣泛的推廣和應(yīng)用。但實(shí)際設(shè)計(jì)過程中,由于該結(jié)構(gòu)的變形機(jī)理與計(jì)算模型尚未有統(tǒng)一的理論,其設(shè)計(jì)方法比較雜亂,存在一定程度的保守設(shè)計(jì)。依托南京梅子洲青奧城地下空間工程及其現(xiàn)場監(jiān)測數(shù)據(jù),采用有限元計(jì)算模型對近年來在地下工程中逐漸得到推廣使用的雙排樁圍護(hù)結(jié)構(gòu)的變形機(jī)理、土壓力分布規(guī)律等進(jìn)行探討和研究;在此基礎(chǔ)上,結(jié)合與目前廣為應(yīng)用的理正雙排樁模型和新規(guī)范雙排樁模型的對比分析,研究并提出更吻合實(shí)際的、精確反映結(jié)構(gòu)變形和力學(xué)特性的雙排樁圍護(hù)結(jié)構(gòu)的修正設(shè)計(jì)計(jì)算模型,該模型架構(gòu)清晰,準(zhǔn)確性高,貼近實(shí)際工程設(shè)計(jì),可以大力推廣應(yīng)用于雙排樁結(jié)構(gòu)的計(jì)算分析。

關(guān)鍵詞:深基坑,雙排樁,變形機(jī)理,計(jì)算模型,有限元法

1概述

近年來,為了適應(yīng)國內(nèi)城市地下工程建設(shè)蓬勃發(fā)展帶來的地下工程密集化、復(fù)雜化等方面的需要,深基坑開挖時雙排樁圍護(hù)結(jié)構(gòu)因具有側(cè)向剛度大、控制變形能力強(qiáng)、施工快速等特點(diǎn)在各種大型深基坑工程中得到了較為廣泛的推廣和應(yīng)用[1-2]。

國內(nèi)從20世紀(jì)90年代初即開始對雙排樁圍護(hù)結(jié)構(gòu)進(jìn)行了初步的探索,隨后大量的研究人員結(jié)合模型試驗(yàn)、工程實(shí)測理論推導(dǎo)等對雙排樁結(jié)構(gòu)的設(shè)計(jì)理論、土壓力分配方法、計(jì)算模型、工程應(yīng)用等進(jìn)行了研究,例如,張弘、林棟[3-4]等人考慮連梁與前、后排樁的鉸接情況,近似地將樁間土體視為受側(cè)向約束的無限長彈性土體,且橫向應(yīng)變?yōu)榱悖瑧?yīng)用經(jīng)驗(yàn)系數(shù)調(diào)整前、后排樁所受的側(cè)向土壓力;何頤華[5]根據(jù)雙排樁前、后排樁之間的滑動土體占樁后滑動土體總量的體積比例來確定前、后排樁所受側(cè)向壓力;熊巨華[6]首次在文獻(xiàn)中提出對于前、后排距小于4倍樁徑的雙排樁支護(hù)結(jié)構(gòu),給出等效抗彎剛度的計(jì)算方法;鄭剛[7]等人提出將雙排樁之間的土視為薄壓縮層,并以水平向彈簧模擬,可以考慮兩排樁之間的土層分布變化、壓縮性、樁間土加固等對雙排樁相互作用的影響,避免對前、后排樁土壓力做出人為的分配。此外,還有彈性地基梁模型[8]、破裂面假設(shè)模型[9-10]和土拱理論模型[11-12]等。

上述研究成果對雙排樁圍護(hù)結(jié)構(gòu)的理論完善和工程應(yīng)用起到了較大的促進(jìn)作用,但是這些現(xiàn)有的不同計(jì)算模型和方法即使在相同計(jì)算參數(shù)的情形下,對同一工程項(xiàng)目得到的計(jì)算結(jié)果仍有較大的差異,不同的計(jì)算模型對雙排樁結(jié)構(gòu)的變形和受力機(jī)理也存在諸多的解釋,給工程設(shè)計(jì)和工程建設(shè)帶來了諸多困惑與困難。

以南京梅子洲基坑工程為背景依托,結(jié)合監(jiān)測數(shù)據(jù)和有限元模型探討了雙排樁結(jié)構(gòu)的變形和受力機(jī)理,在較廣泛應(yīng)用的理正雙排樁模型和新規(guī)范雙排樁模型的基礎(chǔ)上,研究并提出了更吻合實(shí)際的、適用于雙排樁圍護(hù)結(jié)構(gòu)的修正設(shè)計(jì)計(jì)算模型,以更好地服務(wù)于工程實(shí)際。

圖1 雙排樁縱向剖面(單位:mm)

2工程概況

南京梅子洲基坑工程周圍建筑物眾多,為嚴(yán)格控制變形,部分區(qū)域采用雙排樁圍護(hù)結(jié)構(gòu)。如圖1所示,前樁長17.5 m,后樁長13 m,前、后排間距13 m。鋼筋混凝土冠梁尺寸為1.4 m×1.0 m,鋼筋混凝土連梁尺寸為0.8 m×1.0 m,前排樁后采用三軸攪拌樁止水帷幕。基坑懸臂開挖7.9 m,開挖前基坑降水至開挖面以下1 m處。

樁體結(jié)構(gòu)主要涉及4層土,分別為雜填土①、淤泥質(zhì)粉質(zhì)黏土②、粉砂③以及粉細(xì)砂④。場地內(nèi)地表以下10~13 m均為淤泥質(zhì)粉質(zhì)黏土,局部厚度達(dá)23 m。該地層以下以粉細(xì)砂為主,深度達(dá)50 m,之后進(jìn)入強(qiáng)、中風(fēng)化的泥巖地層。

3雙排樁變形機(jī)理和結(jié)構(gòu)受力研究

3.1有限元模型

基于連續(xù)介質(zhì)理論的有限元法是目前解決工程領(lǐng)域數(shù)值計(jì)算最流行的計(jì)算方法[13-14],在合理選擇計(jì)算參數(shù)和本構(gòu)模型的前提下,有限元法得到的計(jì)算結(jié)果具有很好的研究和參考價值。本文采用PLAXIS有限元軟件,建立二維有限元模型,其中前樁和后樁按照抗彎剛度等效為地連墻,連梁等效為連板。有限元模型尺寸為100 m×50 m,左右邊界水平位移約束,底邊水平、豎直雙向位移約束。

土體本構(gòu)模型選用Harding-Soil(HS)模型,相較于Mohr-Coulomb模型,HS模型采用雙曲線應(yīng)力-應(yīng)變關(guān)系,同時剛度模量不斷變化的特性適用于具有典型卸載特色的基坑開挖過程。

結(jié)合該計(jì)算斷面的現(xiàn)場監(jiān)測位移數(shù)據(jù),通過反演分析獲得了表1所示的較合理的有限元地層計(jì)算參數(shù)。利用表1參數(shù)得到的有限元計(jì)算結(jié)果可以對雙排樁結(jié)構(gòu)的計(jì)算參數(shù)作進(jìn)一步分析研究。

表1 土體參數(shù)計(jì)算值

3.2變形機(jī)理

為分析不同間距下雙排樁的受力變形機(jī)理,分別計(jì)算前、后排樁間距為2、4、6 m和8 m的有限元模型,計(jì)算得出的前樁變形如表2所示。可知,雙排樁整體側(cè)向位移較小,控制變形能力強(qiáng)。樁間距為2~6 m時,樁頂位移即為樁身最大位移;樁間距為8 m時,最大位移位于樁頂以下,連梁的錨拉作用較為明顯。綜上,雙排樁間距小時,整體性明顯;樁間距大時,錨拉性明顯。

表2 前樁位移值

雙排樁結(jié)構(gòu)主要受到樁身側(cè)向土壓力作用,通過有限元模型中的接觸單元可計(jì)算出樁身側(cè)向土壓力值,為雙排樁荷載結(jié)構(gòu)模型的土壓力荷載取值提供參考依據(jù)。

由于開挖面以上土壓力對結(jié)構(gòu)變形影響大,本文土壓力研究主要針對開挖面以上樁體部分。前排樁有限元土壓力計(jì)算值如表3所示。后排樁側(cè)向土壓力計(jì)算值如表4所示。

表3 前排樁外側(cè)側(cè)向土壓力 kPa

表4 后排樁側(cè)向土壓力 kPa

由表3可以看出,不同樁間距下前排樁側(cè)向土壓力值相差較小,說明前、后樁間距對前樁土壓力影響很小。樁身上部位移大,土壓力偏向主動土壓力;樁身中下部位移小,土壓力偏向靜止土壓力。側(cè)向土壓力總體處于主動和靜止土壓力之間。

由表4可以看出,樁身上部位移大,土壓力偏于主動,樁身中下部位移小,土壓力偏于靜止。內(nèi)側(cè)與外側(cè)側(cè)向土壓力值相差不大,總體內(nèi)側(cè)值略大。表明后排樁通過樁間土體將側(cè)向土壓力傳遞給前排樁。

4與常用雙排樁計(jì)算模型的對比分析研究

4.1理正模型

該模型根據(jù)土體積比例系數(shù)法[5,15]理論進(jìn)行分析,即前排樁土壓力按開挖面確定,后排樁土壓力按等效開挖面確定。該方法在理正軟件中得到了實(shí)現(xiàn)和推廣應(yīng)用,如圖2所示,此處暫將其稱為理正模型。

圖2 理正模型

理正模型表明前、后排樁共同承擔(dān)土壓力,并通過土壓力分配系數(shù)進(jìn)行分配。開挖面和等效開挖面以下采用接地彈簧模擬土體抗力,其中:

前樁土壓力

(1)

后樁土壓力

(2)

(3)

(4)

式中,ea為作用于樁身土壓力,kPa。

水平地基基床系數(shù)

(5)

式中m——水平基床系數(shù)比例系數(shù),kN/ m4;

z——計(jì)算點(diǎn)距開挖面距離,m。

4.2新規(guī)范模型

圖3 新規(guī)范模型

《建筑基坑技術(shù)規(guī)程》(JGJ120—2012)中對雙排樁結(jié)構(gòu)建議采用如圖3所示的平面剛架模型計(jì)算[1]。

后排樁外側(cè)作用主動土壓力;

前排樁嵌固段土反力

(6)

前、后排樁的樁間土體對樁側(cè)壓力

(7)

(8)

式中ps0——初始土反力強(qiáng)度,kPa;

Δv——前、后樁相對位移,減小為正值,增大為零,mm;

Sy——雙排樁排距,m;

d——樁徑,m。

4.3計(jì)算結(jié)果對比分析

本次研究分別令前、后排樁間距為2、4、6 m和8 m,理正模型、新規(guī)范模型和有限元模型的計(jì)算結(jié)果見圖4~圖6。

圖4 前排樁樁頂水平位移

圖5 前排樁坑底以上最大彎矩

圖6 前排樁坑底以下最大彎矩

圖4表明,理正模型與有限元模型的位移計(jì)算結(jié)果變化規(guī)律一致,但理正結(jié)果整體偏小。新規(guī)范模型中位移隨樁間距增大而增大,與實(shí)際不太符合。

圖5、圖6表明理正模型計(jì)算的結(jié)構(gòu)內(nèi)力與有限元模型計(jì)算結(jié)果變化趨勢一致,但樁間距為2~4 m時偏差較大。新規(guī)范模型內(nèi)力計(jì)算結(jié)果與有限元相差較大。

綜上,從變形機(jī)理和結(jié)構(gòu)受力的角度而言,理正模型與基于監(jiān)測數(shù)據(jù)的有限元模型計(jì)算結(jié)果比較接近,可針對理正模型進(jìn)行修正。

5修正的雙排樁設(shè)計(jì)計(jì)算模型

5.1土壓力修正

通過有限元模型計(jì)算可知(表3、表4),前樁承擔(dān)所有主動土壓力,后樁內(nèi)側(cè)土壓力較外側(cè)大。雙排樁荷載結(jié)構(gòu)模型中可將土壓力全部作用于前樁,后樁不作用土壓力荷載,本次研究據(jù)此結(jié)論對土壓力進(jìn)行修正。

5.2嵌固段彈簧剛度修正

基坑圍護(hù)結(jié)構(gòu)嵌固段水平地基基床系數(shù)K常采用“m”法計(jì)算。其中最為常見的兩種剛度取值是三角形分布與梯形分布。

三角形

(9)

梯形

(10)

沿海地區(qū)土質(zhì)較軟,根據(jù)多個工程實(shí)際經(jīng)驗(yàn),梯形分布臨界深度約5 m左右。

以樁間距4 m的雙排樁為例,計(jì)算上述兩種剛度取值模式下的結(jié)構(gòu)變形和內(nèi)力,如表5所示。

表5 樁間距4 m時計(jì)算值

由表5可知,從變形和內(nèi)力的角度而言,梯形彈簧剛度取值法與有限元模型計(jì)算結(jié)果更加吻合;三角形取值法會使得嵌固段土體剛度較大,結(jié)構(gòu)變形較小。

5.3后樁內(nèi)側(cè)彈簧剛度折減修正

理正模型采用的接地彈簧邊界與新規(guī)范模型采用的樁間連接彈簧是兩種模型間最本質(zhì)的區(qū)別。前面的計(jì)算結(jié)果表明采用接地彈簧更符合結(jié)構(gòu)變形和受力特性,但樁間土體本質(zhì)上是有限空間土體,采用接地彈簧邊界模擬對后排樁的限制作用會造成土抗力過大,與實(shí)際不符。可針對不同前、后排樁間距雙排樁模型,計(jì)算出合理的后排樁內(nèi)側(cè)彈簧剛度折減系數(shù)ω。

經(jīng)過多次計(jì)算,合理的后排樁彈簧剛度折減系數(shù)ω如表6所示。

表6 不同樁間距的ω

5.4雙排樁修正模型

結(jié)合前述研究,基于理正計(jì)算模型,本文對土壓力、嵌固段彈簧剛度取值以及后樁內(nèi)側(cè)彈簧剛度折減等進(jìn)行了修正,提出了相應(yīng)的更符合實(shí)際的“雙排樁修正模型”。

雙排樁修正模型、理正模型、有限元模型的計(jì)算值對比如圖7~圖9所示。結(jié)果表明,雙排樁修正模型的計(jì)算結(jié)果從變形規(guī)律和數(shù)值大小方面相比于理正模型,更加吻合基于監(jiān)測值的有限元模型,適用于普通間距的雙排樁圍護(hù)結(jié)構(gòu)。

圖7 前排樁樁頂水平位移

圖8 前排樁坑底以上最大彎矩

圖9 后排樁樁身最大彎矩

6結(jié)論

本文首先結(jié)合監(jiān)測數(shù)據(jù),通過有限元模型計(jì)算對雙排樁結(jié)構(gòu)受力和變形機(jī)理進(jìn)行研究;在此基礎(chǔ)上,結(jié)合理正模型、新規(guī)范模型和有限元模型的對比分析研究,提出了新的雙排樁修正計(jì)算模型。主要得出以下結(jié)論。

(1)有限元模型計(jì)算結(jié)果表明:雙排樁控制變形能力強(qiáng)。樁間距較小時,體現(xiàn)整體性;樁間距較大時,體現(xiàn)錨拉性。前樁外側(cè)土壓力整體介于主動土壓力和靜止土壓力之間;后樁內(nèi)側(cè)土壓力較外側(cè)稍大。

(2)理正模型、新規(guī)范模型和有限元模型計(jì)算結(jié)果表明:理正模型的變形和內(nèi)力變化趨勢和有限元較為接近,新規(guī)范模型部分規(guī)律不符合有限元模型,可針對理正模型進(jìn)行修正。

(3)將理正模型中土壓力進(jìn)行修正,即前樁開挖面以上外側(cè)承擔(dān)全部主動土壓力,開挖面以下為矩形土壓力荷載;后樁外側(cè)不承擔(dān)土壓力。通過計(jì)算表明,修正土壓力符合雙排樁受力變形特性。

(4)嵌固段彈簧剛度采用“m”法計(jì)算時,梯形分布比三角形分布的計(jì)算結(jié)果更符合雙排樁結(jié)構(gòu)特性。

(5)后樁內(nèi)側(cè)彈簧采取剛度折減時,可使計(jì)算結(jié)果更為準(zhǔn)確。一般前、后樁間距為2~8 m時,折減系數(shù)為0.2~0.5。

(6)雙排樁修正模型能夠較好地反映結(jié)構(gòu)變形和力學(xué)特性,模型簡單,易于實(shí)際工程設(shè)計(jì)應(yīng)用。

參考文獻(xiàn):

[1]中華人民共和國住房和城鄉(xiāng)建設(shè)部.JGJ120—2012建筑基坑支護(hù)技術(shù)規(guī)程[S].北京:中國建筑工業(yè)出版社,2012.

[2]張勉,薛光橋,蔡永昌.拉錨式雙排樁荷載結(jié)構(gòu)模型計(jì)算研究[C]∥中國建筑業(yè)協(xié)會深基礎(chǔ)施工分會,中國土木工程學(xué)會土力學(xué)及巖土工程分會樁基礎(chǔ)學(xué)術(shù)委員會,中國工程機(jī)械學(xué)會樁工機(jī)械分會.第五屆全國深基礎(chǔ)工程發(fā)展論壇會議論文集,2015:41-46.

[3]張弘.深基坑開挖中雙排支護(hù)結(jié)構(gòu)的應(yīng)用與探討[J].地基處理,1993,4(3):42-47.

[4]林棟.用雙排灌注樁作深基坑圍護(hù)結(jié)構(gòu)的嘗試[J].建筑施工,1994(4):18-19.

[5]何頤華,楊斌,金寶森,等.雙排護(hù)坡樁試驗(yàn)與計(jì)算的研究[J].建筑結(jié)構(gòu)學(xué)報,1996(2):58-66.

[6]熊巨華.一類雙排樁支護(hù)結(jié)構(gòu)的簡化計(jì)算方法[J].勘察科學(xué)技術(shù),1999(2):32-34.

[7]鄭剛,李欣,劉暢,等.考慮樁土相互作用的雙排樁分析[J].建筑結(jié)構(gòu)學(xué)報,2004(1):99-106.

[8]劉釗.雙排支護(hù)樁結(jié)構(gòu)的分析及試驗(yàn)研究[J].巖土工程學(xué)報,1992(5):76-80.

[9]程知言,裘慰倫,張可能,等.雙排樁支護(hù)結(jié)構(gòu)設(shè)計(jì)計(jì)算方法探討[J].地質(zhì)與勘探,2001(2):88-90,93.

[10]戴智敏,陽凱凱.深基坑雙排樁支護(hù)結(jié)構(gòu)體系受力分析與計(jì)算[J].信陽師范學(xué)院學(xué)報:自然科學(xué)版,2002(3):348-352.

[11]蔡海兵,黃以春,李陽.基于隨機(jī)介質(zhì)理論的土體融沉預(yù)測及其參數(shù)敏感性分析[J].鐵道標(biāo)準(zhǔn)設(shè)計(jì),2015(8):107-111.

[12]侯新宇,劉松玉,童立元.主、被動土壓區(qū)加固對地鐵深基坑變形影響的數(shù)值分析[J].鐵道標(biāo)準(zhǔn)設(shè)計(jì),2012(4):94-97.

[13]焦蒼,范鵬,徐成家.深基坑支護(hù)結(jié)構(gòu)變形預(yù)測研究與應(yīng)用[J].鐵道標(biāo)準(zhǔn)設(shè)計(jì),2005(1):81-84.

[14]孫毅,戴李春,齊俊,等.深淺式異形基坑開挖工法比選與地表沉降特點(diǎn)研究[J].鐵道標(biāo)準(zhǔn)設(shè)計(jì),2015(9):125-129.

[15]曠慶華.同深基坑開挖引起緊鄰地鐵車站變形特性研究[J].鐵道標(biāo)準(zhǔn)設(shè)計(jì),2015(5):130-132.

收稿日期:2015-11-11; 修回日期:2016-01-04

作者簡介:薛光橋(1979—),男,高級工程師,2005年畢業(yè)于武漢理工大學(xué)結(jié)構(gòu)工程系,工學(xué)碩士,E-mail:21987532@qq.com。

文章編號:1004-2954(2016)07-0116-05

中圖分類號:U452

文獻(xiàn)標(biāo)識碼:A

DOI:10.13238/j.issn.1004-2954.2016.07.027

Research on Mechanism and Computation Model of Deep Foundation Pit Double-row Pile Structure

XUE Guang-qiao1,2

(1.China Railway Siyuan Survey and Design Group Co., Ltd., Wuhan 430000, China;2.Hubei Provincial Engineering Laboratory for Underwater Tunnel, Wuhan 430000, China)

Abstract:In recent years, the double-row piles structure is widely promoted and applied in deep foundation pit to meet the needs of intensive and complicated underground projects due to its high lateral stiffness, less deformation and quick construction and the vigorous development of urban underground engineering construction. However, in practical design, there is no unified design method and there exists a certain degree of conservative due to the lack of standard theory to match the structure deformation mechanism with calculation models. With Nanjing Meizizhou Qingao city underground space engineering program and field monitoring data, the deformation mechanism and soil pressure distribution law of double-row pile retaining structure are studied, which is gradually promoted and used in underground engineering field in recent years with finite element computation model. On the basis of this, a practical and reasonable corrected computation model of double-row pile retaining structure is recommended after the comparison and analysis of the widely used Lizheng and Xinguifang double-row pile models. This model structure is featured by clearness, high accuracy, and practical application in engineering design. The model can be easily applied in computation and analysis of double-row pile structure.

Key words:Deep foundation pit; Double-row pile; Deformation mechanism; Computation model; Finite element method

猜你喜歡
有限元法
正交各向異性材料裂紋疲勞擴(kuò)展的擴(kuò)展有限元法研究
基于有限元法的高頻變壓器繞組損耗研究
基于有限元法副發(fā)動機(jī)托架輕量化設(shè)計(jì)
專用汽車(2016年8期)2016-03-01 04:16:43
傳遞矩陣法與有限元法計(jì)算電機(jī)轉(zhuǎn)子臨界轉(zhuǎn)速的對比分析
Sine-Gordon方程H1-Galerkin非協(xié)調(diào)混合有限元法的誤差分析
三維有限元法在口腔正畸生物力學(xué)研究中發(fā)揮的作用
RKDG有限元法求解一維拉格朗日形式的Euler方程
集成對稱模糊數(shù)及有限元法的切削力預(yù)測
有限元法在機(jī)械設(shè)計(jì)方向中的教學(xué)實(shí)踐
基于HCSR和CSR-OT的油船疲勞有限元法對比分析
船海工程(2013年6期)2013-03-11 18:57:25
主站蜘蛛池模板: 午夜福利在线观看成人| 久久青青草原亚洲av无码| 91蜜芽尤物福利在线观看| 欧美亚洲第一页| 五月婷婷综合网| 人妻精品久久久无码区色视| 国产精品久久久久久久久久98 | 亚洲成网站| 干中文字幕| 午夜a视频| 99精品欧美一区| 国产在线91在线电影| a天堂视频| 538国产在线| 午夜人性色福利无码视频在线观看 | 国产91导航| 色亚洲成人| 黄色在线网| 自拍中文字幕| 国产成人麻豆精品| 最新国产高清在线| 亚洲天堂精品视频| 高清无码不卡视频| 韩日午夜在线资源一区二区| 免费看a毛片| 国产丝袜啪啪| 国产v精品成人免费视频71pao| 91久久国产综合精品女同我| 99久久人妻精品免费二区| 国模视频一区二区| 免费在线成人网| 国产精品一区在线观看你懂的| 亚洲成在线观看 | 亚洲欧美在线综合一区二区三区| 日韩第九页| 久久久久九九精品影院| 亚洲女同欧美在线| 久久国产精品影院| 国内精品九九久久久精品| 午夜日韩久久影院| 亚洲啪啪网| 丰满人妻久久中文字幕| 免费精品一区二区h| 精品国产香蕉在线播出| www欧美在线观看| 欧美成人精品在线| 在线色综合| 伊人久久综在合线亚洲2019| 尤物午夜福利视频| 亚洲精品欧美重口| 91精品国产一区自在线拍| 亚洲系列无码专区偷窥无码| 青草精品视频| 免费人成又黄又爽的视频网站| 伊人成人在线| 免费看一级毛片波多结衣| 亚洲国产综合第一精品小说| 91久久国产综合精品女同我| 婷婷久久综合九色综合88| 亚洲狼网站狼狼鲁亚洲下载| 久久亚洲精少妇毛片午夜无码| 国产高清不卡| 亚洲香蕉伊综合在人在线| 无码高清专区| 91啪在线| 亚洲一区免费看| 在线看免费无码av天堂的| 日本久久久久久免费网络| 欧美无遮挡国产欧美另类| www成人国产在线观看网站| 免费看黄片一区二区三区| 91免费国产高清观看| 国产婬乱a一级毛片多女| 丁香六月激情婷婷| 啦啦啦网站在线观看a毛片 | 亚洲无线观看| 91热爆在线| 亚洲欧美日韩动漫| 成人在线观看一区| 久久久久九九精品影院| 成人国产精品视频频| 99999久久久久久亚洲|